Добавил:
kiopkiopkiop18@yandex.ru t.me/Prokururor I Вовсе не секретарь, но почту проверяю Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Ординатура / Хирургия / Библиотека им академика М.И. Перельмана / Книга_5395_Библиотеки_им_академика_М_И_Перельмана

.pdf
Скачиваний:
0
Добавлен:
15.09.2026
Размер:
15 Мб
Скачать
☆
Overview of Scheduling Methods for Pharmaceutical Production 367
https://t.me/medicina_free
time might constrain the production, and incorporating these constraints is essential to generate an implementable schedule. In some chemical or biopharmaceutical production, products are generated in batches, and different batches of the same products/intermediates cannot be mixed, which demands multiple storage vessels for storing these batches. Furthermore, if batching decisions are considered, the size of the storage vessels will influence batch sizes. If enough vessels of sufficient sizes are available, then we can assume unlimited storage policy; else, the storage policy is limited, and the multistage production schedule is influenced by both the number and size of the available storage vessels. Another important factor is storage time, which is often important for pharmaceutical processes due to the low shelf life of the intermediate products. Accounting for the waiting time in the processing units is also important and might be used instead of an actual storage vessel in case of limited availability. For further information on modelling of the storage constraints in a multistage facility, the interested reader can refer to Sundaramoorthy et al. [19]
5 Illustrative Example
In this section, we demonstrate the applicability of the abovementioned discrete­time scheduling models through some case studies. The instances are designed to demonstrate the modelling of changeovers and utility (resource) constraints and the efficacy of discrete-time models. We have solved all the instances using the FICO Xpress Optimization Suite and the “mmxprs” module version 2.8.1 on an Intel Core i7 (2.3 GHz) unit with 32GB RAM [43], and the model characteristics for each of these instances are provided in Table 1.
For the first two instances, we consider a production facility involving two production stages and four units: two units per stage. Ten batches need to be processed in this facility, and unlimited storage policy is assumed. Data for this example is taken from Sundaramoorthy et al. [19] and modified accordingly (Figs.
2 and 3).
Two different instances ofthis example are considered. In thefirst instance, wedo not consider changeover times, while in the second, we find the optimum schedule incorporating changeover times. The objective is to minimize the makespan. The resultant optimal schedules for these two instances are shown in Fig. 2a and b.
Table 1 Model statistics for the instances
Objective Makespan minimization Makespan minimization Cost minimization Constraints 604 10,659 51,831 Continuous variables 1 1 99 Binary variables 890 890 80,661 Computational time (s) 0.1 3 2309
Instance 1 Instance 2 Instance 3
368 S. Misra and C. T. Maravelias
https://t.me/medicina_free
Fig. 2 Gantt chart of optimal solutions for Instances 1 and 2
Fig. 3 Gantt chart and utility profiles of optimal solutions for Instance 3. (a) Optimum schedule;
(b) consumption and availability of utility (EL) required for batch execution; (c) consumption of utility (CW) required for changeover
Overview of Scheduling Methods for Pharmaceutical Production 369
https://t.me/medicina_free
We next present a large size instance, in which 20 batches are to be scheduled. Each batch is to be processed in two stages, and each stage has two dedicated units in which the batches can be processed. We also consider utility requirements and illustrate the effect of limited utility on the optimum schedule. We assume that the processing of each batch requires electricity (EL) in both stages. Furthermore, it is assumed that each changeover among batches incurs cost and requires the consumption of cleaning water (CW), which is assumed to be a nonrenewable utility with initial inventory of 130 tons. The cost of CW is assumed to be included in the changeover cost. The objective is to minimize the total (processing +changeover + utility) cost.
Figure 3 shows the Gantt chart and utility profiles of the optimum solution for the third instance. Figure 3b depicts the engagements of EL utility that coincide with the processing of batches, whereas Fig. 3c illustrates the consumption profile of CW utility, which is triggered by the changeovers. This example demonstrates how the intricacies related to changeover time cost and utility requirements can be modeled. Nevertheless, as the computational requirements increase with problem size, advanced solution methods may be necessary to solve large industrial-scale instances. A set of such solution methods, including tightening methodologies and reformulation techniques applicable to multistage, multiproduct production environment, are described in Merchan et al. [20] and Lee and Maravelias [21, 22].
6 Conclusion
This chapter provides an overview of discrete-time batch scheduling models that can be used to schedule chemical and biopharmaceutical manufacturing environments. Concepts and techniques regarding the formulation of such models are discussed, starting from the single-unit problem and building up to problems in multiproduct multistage facilities. The structure of the discrete-time models remains almost the same from a single-unit to multiunit multistage problems, and complex decisions related to batching, changeover, and shared resources can be easily modeled. Furthermore, as the time grid is explicit, time-sensitive nuances such as time varying resource availability, intermediate orders can also be easily incorporated into the model. We showed, through a medium-size case study of a multiproduct multistage facility, how the aforementioned features can be incorporated.
References
1. Marques, C. M., Moniz, S., de Sousa, J. P., Barbosa-Povoa, A. P. & Reklaitis, G. Decision­support challenges in the chemical-pharmaceutical industry: Findings and future research directions. Comput. Chem. Eng.134, 106672 (2020).
370 S. Misra and C. T. Maravelias
https://t.me/medicina_free
2. Lee, S. L. et al. Modernizing Pharmaceutical Manufacturing: from Batch to Continuous Production. J. Pharm. Innov.10, 191–199 (2015).
3. Georgiadis, G. P., Elekidis, A. P. & Georgiadis, M. C. Optimization-Based Scheduling for the Process Industries: From Theory to Real-Life Industrial Applications. Processes vol. 7 (2019).
4. Reklaitis, G. V. Overview of Scheduling and Planning of Batch Process Operations BT - Batch Processing Systems Engineering. in (eds. Reklaitis, G. V, Sunol, A. K., Rippin, D. W. T. & Hortaçsu, Ö.) 660–705 (Springer Berlin Heidelberg, 1996).
5. Kallrath, J. Planning and scheduling in the process industry. OR Spectr.24, 219–250 (2002).
6. Floudas, C. a. & Lin, X. Continuous-time versus discrete-time approaches for scheduling of chemical processes: A review. Comput. Chem. Eng.28, 2109–2129 (2004).
7. Méndez, C. a., Cerdá, J., Grossmann, I. E., Harjunkoski, I. & Fahl, M. State-of-the-art review of optimization methods for short-term scheduling of batch processes. Comput. Chem. Eng.30, 913–946 (2006).
8. Maravelias, C. T. General framework and modeling approach classification for chemical production scheduling. AIChE J.58, 1812–1828 (2012).
9. Harjunkoski, I. et al. Scope for industrial applications of production scheduling models and solution methods. Comput. Chem. Eng.62, 161–193 (2014).
10. Méndez, C. A., Henning, G. P. & Cerdá, J. Optimal scheduling of batch plants satisfying multiple product orders with different due-dates. Comput. Chem. Eng.24, 2223–2245 (2000).
11. Méndez, C. A., Henning, G. P. & Cerdá, J. An MILP continuous-time approach to short­term scheduling of resource-constrained multistage flowshop batch facilities. Comput. Chem. Eng.25, 701–711 (2001).
12. Gupta, S. & Karimi, I. A. An Improved MILP Formulation for Scheduling Multiproduct, Multistage Batch Plants. Ind. Eng. Chem. Res.42, 2365–2380 (2003).
13. Sundaramoorthy, A. & Maravelias, C. T. Modeling of Storage in Batching and Scheduling of Multistage Processes. Ind. Eng. Chem. Res.47, 6648–6660 (2008).
14. Sundaramoorthy, A. & Maravelias, C. T. Simultaneous Batching and Scheduling in Multistage Multiproduct Processes. Ind. Eng. Chem. Res.47, 1546–1555 (2008).
15. Kopanos, G. M., Laínez, J. M. & Puigjaner, L. An Efficient Mixed-Integer Linear Programming Scheduling Framework for Addressing Sequence-Dependent Setup Issues in Batch Plants. Ind. Eng. Chem. Res.48, 6346–6357 (2009).
16. Cerdá, J., Henning, G. P. & Grossmann, I. E. A Mixed-Integer Linear Programming Model for Short-Term Scheduling of Single-Stage Multiproduct Batch Plants with Parallel Lines. Ind. Eng. Chem. Res.36, 1695–1707 (1997).
17. Kondili, E., Pantelides, C. C. & Sargent, R. W. H. A general algorithm for short-term scheduling of batch operations—I. MILP formulation. Comput. Chem. Eng.17, 211–227 (1993).
18. Shah, N., Pantelides, C. C. & Sargent, R. W. H. A general algorithm for short-term scheduling of batch operations—II. Computational issues. Comput. Chem. Eng.17, 229–244 (1993).
19. Sundaramoorthy, A., Maravelias, C. T. & Prasad, P. Scheduling of Multistage Batch Processes under Utility Constraints. Ind. Eng. Chem. Res.48, 6050–6058 (2009).
20. Merchan, A. F., Lee, H. & Maravelias, C. T. Discrete-time mixed-integer programming models and solution methods for production scheduling in multistagefacilities. Comput. Chem. Eng.94, 387–410 (2016).
21. Lee, H. & Maravelias, C. T. Discrete-time mixed-integer programming models for short-term scheduling in multipurpose environments. Comput. Chem. Eng.107, 171–183 (2017).
22. Lee, H. & Maravelias, C. T. Mixed-integer programming models for simultaneous batching and scheduling in multipurpose batch plants. Comput. Chem. Eng.106, 621–644 (2017).
23. Schilling, G. & Pantelides, C. C. A simple continuous-time process scheduling formulation and a novel solution algorithm. Comput. Chem. Eng.20, S1221–S1226 (1996).
24. Castro, P. M., Barbosa-Póvoa, A. P., Matos, H. A. & Novais, A. Q. Simple Continuous-Time Formulation for Short-Term Scheduling of Batch and Continuous Processes. Ind. Eng. Chem. Res.43, 105–118 (2004).
Overview of Scheduling Methods for Pharmaceutical Production 371
https://t.me/medicina_free
25. Maravelias, C. T. & Grossmann, I. E. New General Continuous-Time State−Task Network Formulation for Short-Term Schedulingof Multipurpose Batch Plants. Ind. Eng. Chem. Res.42, 3056–3074 (2003).
26. Janak, S. L., Floudas, C. A., Kallrath, J. & Vormbrock, N. Production scheduling of a large­scale industrial batch plant. I. Short-term and medium-term scheduling. Ind. Eng. Chem. Res.45, 8234–8252 (2006).
27. Shaik, M. A. & Floudas, C. A. Unit-specific event-based continuous-time approach for short­term scheduling of batch plants using RTN framework. Comput. Chem. Eng.32, 260–274 (2008).
28. Mostafaei, H. & Harjunkoski, I. Continuous-time scheduling formulation for multipurpose batch plants. AIChE J.66, e16804 (2020).
29. Ierapetritou, M. G. & Floudas, C. A. Effective Continuous-Time Formulation for Short-Term Scheduling. 1. Multipurpose Batch Processes. Ind. Eng. Chem. Res.37, 4341–4359 (1998).
30. Maravelias, C. T. Mixed-Time Representation for State-TaskNetwork Models. Ind. Eng. Chem. Res.44, 9129–9145 (2005).
31. Westerlund, J., Hästbacka, M., Forssell, S. & Westerlund, T. Mixed-Time Mixed-Integer Linear Programming Scheduling Model. Ind. Eng. Chem. Res.46, 2781–2796 (2007).
32. Lee, H. & Maravelias, C. T. Combining the advantages of discrete- and continuous-time scheduling models: Part 1. Framework and mathematical formulations. Comput. Chem. Eng.116, 176–190 (2018).
33. Castro, P. M., Harjunkoski, I. & Grossmann, I. E. Optimal Short-Term Scheduling of Large­Scale Multistage Batch Plants. Ind. Eng. Chem. Res.48, 11002–11016 (2009).
34. Kopanos, G. M., Méndez, C. A. & Puigjaner, L. MIP-based decomposition strategies for large­scale scheduling problems in multiproduct multistage batch plants: A benchmark scheduling problem of the pharmaceutical industry. Eur. J. Oper. Res.207, 644–655 (2010).
35. Stefansson, H., Sigmarsdottir, S., Jensson, P. & Shah, N. Discrete and continuous time representations and mathematical models for large production scheduling problems: A case study from the pharmaceutical industry. Eur. J. Oper. Res.215, 383–392 (2011).
36. Lakhdar, K., Zhou, Y., Savery, J., Titchener-Hooker, N. J. & Papageorgiou, L.G. Medium Term Planning of Biopharmaceutical Manufacture using Mathematical Programming. Biotechnol. Prog.21, 1478–1489 (2005).
37. Kabra, S., Shaik, M. A. & Rathore, A. S. Multi-period scheduling of a multistage multiproduct bio-pharmaceutical process. Comput. Chem. Eng.57, 95–103 (2013).
38. Liu, S., Yahia, A. & Papageorgiou, L. G. Optimal Production and Maintenance Planning of Biopharmaceutical Manufacturing under Performance Decay. Ind. Eng. Chem. Res.53, 17075– 17091 (2014).
39. Moniz, S., Barbosa-Póvoa, A. P., de Sousa, J. P. & Duarte, P. Solution Methodology for Scheduling Problems in Batch Plants. Ind. Eng. Chem. Res.53, 19265–19281 (2014).
40. Eberle, L. et al. Rigorous approach to scheduling of sterile drug product manufacturing. Comput. Chem. Eng.94, 221–234 (2016).
41. Zyngier, D. & Kelly, J. D. Multi-Product Inventory Logistics Modeling in the Process Indus­tries BT - Optimization and Logistics Challenges in the Enterprise. in (eds. Chaovalitwongse, W., Furman, K. C. & Pardalos, P. M.) 61–95 (Springer US, 2009). doi: https://doi.org/10.1007/
978-0-387-88617-6_2.
42. Velez, S. & Maravelias, C. T. Advances in Mixed-Integer Programming Methods for Chemical Production Scheduling. Annu. Rev. Chem. Biomol. Eng.5, 97–121 (2014).
43. Guéret, C., Prins, C., Sevaux, M. & Heipcke, S. (Dash Optimization Limited, 2002).
Applications of Optimization with Xpress-MP.
Model-Based Risk Assessment of mAb
https://t.me/medicina_free
Developability
M. Karlberg, A. Kizhedath, and J. Glassey
Monoclonal antibodies were already one of the fastest growing sectors of bio­pharmaceutical industry [1]. Recent research on the significant benefits of various antibodies in reducing the risks of fatality or reducing the symptoms of COVID­19, e.g., tocilizumab and sarilumab (Cortegiani et al. 2021), inevitably increases the importance of the rapid discovery of mAbs and the development of efficient manufacturing processes. Research reports frequently concentrated on the rapid discovery of new mAbs, but the developability and the manufacturability of mAbs were less explored until recently [1–4]. Thischapter adopts the framework of quality by design (QbD), concentrating particularly on the model-based risk assessment of mAb developability. A case study highlights specific areas where advanced modelling approaches can contribute to speeding up the manufacturability and developability of mAbs and demonstrates the benefits and the challenges of this approach.
1 Quality by Design
The QbD paradigm was introduced in 2004 and is a systematic approach that aligns with the Process Analytical Technology (PAT) principles and aims to build quality into the product through product and process understanding (US FDA 2004). The framework is especially useful for the process development of mAbs, which consists of many different steps (unit operations). A typical mAb process can be divided into two parts: the upstream (USP) or the cell culture processing where the mAbs are expressed and the downstream (DSP) or purification where the mAbs are isolated
M. Karlberg · A. Kizhedath · J. Glassey () School of Engineering, Newcastle University, Newcastle upon Tyne, UK e-mail: jarka.glassey@ncl.ac.uk
© The Author(s), under exclusive license to Springer Nature Switzerland AG 2022 A. Fytopoulos et al. (eds.), Optimization of Pharmaceutical Processes, Springer Optimization and Its Applications 189, https://doi.org/10.1007/978-3-030-90924-6_14
373
374 M. Karlberg et al.
https://t.me/medicina_free
Fig. 1 General outline of the QbD methodology. The process design space is shown as the dashed box where the effects of process parameters and raw material input (blue box) on the product quality are characterised.Steps highlighted in red indicate risk assessment of either product quality attributes, process parameters or raw materials. The green box indicates the availability of clinical data which can be used to better define the QTPP
and contaminants removed. Typically, a mAb process will consist of between 15 and 20 different unit operations which must be individually characterised, e.g. finding the optimal process parameters that deliver consistent drug quality and safety of use [5]. The major steps of a QbD implementation to develop and characterise a product/process are illustrated in Fig. 1. It is important to note that these are a subset of steps involved in the whole QbD process. For more details on QbD, including important aspects such as Critical Material Attributes (CMAs), which will not be discussed in this contribution in detail, please refer to Rathore [6] and Lawrence et al. [7].
1.1 Quality Target Product Profile
The implementation of QbD starts by defining the Quality Target Product Profile (QTPP) which forms the basis of the design for the development and contains information about the drug quality criteria such as delivery mechanisms, intended use and route of administration for the intended product to ensure clinical safety and efficacy. The QTPP is generated from knowledge that is based on literature research, clinical trials and existing experience from industry or academia [8, 9]. For mAbs, the QTPP relates to the product’s intended use and properties that can affect patients and needs to be clearly stated to avoid adverse effects in patients.
Model-Based Risk Assessment of mAb Developability 375
https://t.me/medicina_free
Issues with the lack of efficacy, adverse effects and/or high manufacturing costs have led to the withdrawal or discontinuation of several mAb products (Kizhedath et al. 2017). It is thus important that QTPPs include antigen binding, pharmacokinetics, effector function, stability and half-life of the mAb [6, 10]. However, much of this information does not become available until later when clinical data has been obtained. Thus, instead, many aspects of the QTPP are based on prior knowledge in early process development of a mAb. Recently, the computational prediction and simulation of the mAb structure have been shown to be a valuable tool for mAb design due to their ability to provide estimates of behaviour and protein stability which can aid in more accurate QTPP specification [11, 12].
1.2 Critical Quality Attributes and Risk Assessment
QTPP provides the basis for the identification of Critical Quality Attributes (CQAs) from a list of Quality Attributes (QAs) using risk-based analysis. This is in accordance with the ICH Q9 guideline to investigate properties that might affect product quality. The CQAs are physical, chemical or biological properties of the drug product that need to be within appropriate ranges to ensure the desired product quality. These ranges, similar to the generation of QTPP, are obtained through literature research, clinical data and previous experience, but they are also updated during process development as new information from characterisation studies becomes available. The most frequently used method for risk assessment in industries is Failure Mode and Effect Analysis (FMEA) where the impact of different unit operations in the process on the QAs is listed. Each effect is ranked according to a Severity rating (S), an Occurrence rating (O) and a Detectability rating (D). A final Risk Priority Number (RPN) is calculated by multiplying the ratings which are then ranked to identify the effects that potentially affect the product quality and efficacy [13, 14]. Tailored risk assessment methods have also been proposed for biopharmaceuticals by Zalai et al. (2013). In their work, the authors argued that traditional methods do not consider the “complexity” of how a process might affect the product or the “uncertainty” which includes the quality of the input material as a possible source of risk and which need to be added as additional factors to the risk assessment.
Once the CQAs have been selected, a process design space is defined by screening process parameters (PPs) for each of the unit operations that have a significant effect on the CQAs. PPs that have a significant impact on the CQAs are called Critical Process Parameters (CPPs) and are identified and controlled using the following steps:
1. Similar to identification of CQAs, risk analysis methods, such as FMEA, are used
to reduce the large number of PPs to those that may affect CQAs.
376 M. Karlberg et al.
https://t.me/medicina_free
2. Systematic experimental studies using Design of Experiments (DoE) over a
range of PP settings are carried out in small scale to obtain experimental data
for process characterisation to identify CPPs and their optimal ranges. This is
referred to as the control space.
3. Multivariate data analysis (MVDA) is used for the implementation of appropriate
real-time monitoring and control strategies needed for the defined CQAs and
CPPs to ensure product quality. Movement outside of the defined control space
would cause the product quality to drop below that of the desired quality stated
in the QTPP.
The use of statistical DoE is preferred over univariate analysis in process development of pharmaceuticals as it can generate qualitative and quantitative information about important process parameters and their impact on the product quality [15]. Response Surface Modelling (RSM) and leverage plots are often used to analyse the DoE data to investigate the significance of PPs on the explored CQAs as well as to define allowed ranges for the identified CPPs [16]. Various experimental designs exist, and selecting an appropriate design is critical in order to maximise the information gained from the experiments. Kumar et al. [17] compared different experimental designs for the DoE of downstream unit operations to demonstrate how these affect the response surface of each unit operation [17]. Tai et a l. [ 18] showed that a well-chosen experimental design can lead to diverse and informative data about the system and when combined with high-throughput experimentation techniques, it can be a powerful tool in defining the process design space.
The fundamental principle of the QbD framework is to increase process under­standing in terms of the effect that PPs have on product quality. Zurdo [19] suggested that the QbD framework needed to be extended to incorporate product understanding in terms of the developability of the pharmaceutical which would include manufacturability, safety, pharmacology and biological activity. The author argued that by using in silico risk assessment tools based on structural features of the mAbs and historic development data, predictions concerning manufacturability of an mAb could be made. In a later publication, two case studies were presented where structural properties of mAbs were successfully linked to CQAs related to aggregation and half-life [4]. Such tools can add great value to early process development of mAbs when implementing the QbD framework where little is known about both the process and product. Thus, a more in-depth investigation of Quantitative Structure-Activity Relationship (QSAR) framework and its potential benefits for QbD integration is explored here.
2 Advanced Modelling Approaches in QbD
The QSAR framework relates structural properties and features (also known as descriptors) of a compound to biological or physicochemical activity [20, 21].
Model-Based Risk Assessment of mAb Developability 377
https://t.me/medicina_free
This methodology was first introduced by Hammet in the 1930s and was later refined by Hansch and Fujita and has become a standard tool for small-molecule drug discovery [22]. A method derived from QSAR, referred to as Quantitative Sequence-Activity Modelling (QSAM), has been introduced in recent years and focuses on relating structural descriptors of proteins, peptides and nucleic acids to activity [23]. The only difference between QSAR and QSAM is the development of descriptors, whereas the guidelines for the predictive model development remain the same. Given the complex protein structure of mAbs, the QSAM methodology for descriptor generation will be of more relevance, and the workflow described below will therefore focus more on protein-based rather than small molecule-based QSAR.
2.1 Descriptor Generation
One of the most important steps in QSAR is the numerical representation of structures of the pharmaceuticals so that they can be used in correlation studies with prediction outputs of interest. For proteins, such as mAbs, two approaches to generate descriptors are discussed here: (1) descriptors generated from the amino acid primary sequence and (2) descriptors generated from three-dimensional models of the mAbs. It has been shown that a combination of both physicochemical and 3D structure descriptors works best and also ensures that the model is not overly reliant on a single type of a descriptor [24].
2.2 Amino Acid Composition-Based Descriptor Generation
Extensive research has been carried out to develop new informative descriptors for peptides and proteins generated from their primary sequence [25]. This was first introduced by Sneath [26] who derived amino acid descriptors for the 20 naturally occurring amino acids from qualitative data. Later on, Kidera et al. [27] used 188 properties of the 20 naturally occurring amino acids, which were converted into ten orthogonal new descriptors to describe the amino acids. Later, the Z-scale, which consists of three new amino acid descriptors derived by applying PCA to 29 physiochemical properties [28, 29], was introduced. Other amino acid scales, which were also derived through PCA, include the extended Z-scale and T-scale [30, 31]. Other descriptors include the so-called isotropic surface area (ISA) and the electronic charge index (ECI), which are derived from the 3D structures of the amino acids [32]. All these descriptors were tested and performed well in respective studies on small peptides. In a two-part review by van Westen et al. [33, 34], many of the existing amino acid scales were benchmarked and compared. The authors demonstrated that the different scales described different physiochemical and topological properties which is useful when deciding on which scales to use [33,
34]. Doytchinova et al. [35] applied the Z-scale’s descriptors to successfully predict