Добавил:
Sekretar
kiopkiopkiop18@yandex.ru
t.me/Prokururor I Вовсе не секретарь, но почту проверяю
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:Ординатура / Хирургия / Библиотека им академика М.И. Перельмана / Книга_5395_Библиотеки_им_академика_М_И_Перельмана
.pdf
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 discretetime 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. Decisionsupport 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 shortterm 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 largescale 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 shortterm 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 LargeScale 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 largescale 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 Industries 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 biopharmaceutical industry [1]. Recent research on the significant benefits of various
antibodies in reducing the risks of fatality or reducing the symptoms of COVID19, 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 understanding 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
Соседние файлы в папке Библиотека им академика М.И. Перельмана
