Straightforward method for calibration of mechanistic cation exchange chromatography models for industrial applications

An examination of process models and model risk frameworks for pharmaceutical manufacturing

Process models are a growing tool for pharmaceutical manufacturing process design and control. The Industry 4.0 paradigm promises to increase the amount of data available to understand manufacturing processes. Tools such as Artificial Intelligence (AI) might accelerate process development and allow better predictions of process trajectories. Several examples of process improvements realized through the application of process models have been shown in lyophilization, chromatography, fluid bed drying, bioreactor control, continuous direct compression, and wet granulation. An important consideration of implementing a process model is determining the impact of the model on the quality of the product and the risks associated with model maintenance over the product lifecycle. Several regulatory documents address risk-based considerations for process models. This work discusses existing risk-based frameworks for model validation and lifecycle maintenance that could aid the adoption of process models in pharmaceutical manufacturing. Hypothetical case studies illustrate the implications of applying a model risk framework to facilitate model validation and lifecycle maintenance in the manufacture of pharmaceuticals and biological products.

From protein structure to an optimized chromatographic capture step using multiscale modeling

Optimizing a biopharmaceutical chromatographic purification process is currently the greatest challenge during process development. A lack of process understanding calls for extensive experimental efforts in pursuit of an optimal process. In silico techniques, such as mechanistic or data driven modeling, enhance the understanding, allowing more cost-effective and time efficient process optimization. This work presents a modeling strategy integrating quantitative structure property relationship (QSPR) models and chromatographic mechanistic models (MM) to optimize a cation exchange (CEX) capture step, limiting experiments. In QSPR, structural characteristics obtained from the protein structure are used to describe physicochemical behavior. This QSPR information can be applied in MM to predict the chromatogram and optimize the entire process. To validate this approach, retention profiles of six proteins were determined experimentally from mixtures, at different pH (3.5, 4.3, 5.0, and 7.0). Four proteins at different pH's were used to train QSPR models predicting the retention volumes and characteristic charge, subsequently the equilibrium constant was determined. For an unseen protein knowing only the protein structure, the retention peak difference between the modeled and experimental peaks was 0.2% relative to the gradient length (60 column volume). Next, the CEX capture step was optimized, demonstrating a consistent result in both the experimental and QSPR-based methods. The impact of model parameter confidence on the final optimization revealed two viable process conditions, one of which is similar to the optimization achieved using experimentally obtained parameters. The multiscale modeling approach reduces the required experimental effort by identification of initial process conditions, which can be optimized.

Real-time monitoring of gradient chromatography using dual Kalman-filters

Real-time state estimation in chromatography is a useful tool to improve monitoring of biopharmaceutical downstream processes, combining mechanistic model predictions with real-time data acquisition to obtain an estimation that surpasses that of either approach individually. One common technique for real-time state estimation is Kalman filtering. However, non-linear adsorption isotherms pose a significant challenge to Kalman filters, which are dependent on fast algorithm execution to function. In this work, we apply Kalman filtering of non-constant elution conditions using a non-linear adsorption isotherm using a novel approach where dual Kalman filters are used to estimate the states of the adsorption modifier, salt, and the components to be separated. We performed offline tuning of the Kalman filters on real chromatogram data from a linear gradient, ion-exchange separation of two proteins. The tuning was then validated by running the Kalman filters in parallel with a chromatographic separation in real time. The resulting, tuned, dual Kalman filters improved the L2 norm by 53 % over the open-loop model prediction, when compared to the true elution profiles. The Kalman filters were also applicable in real-time with a signal sampling frequency of 5 s, enabling accurate and robust estimation and paving the way for future applications beyond monitoring, such as real-time optimal pooling control.

Modeling multi-component separation in hydrophobic interaction chromatography with improved parameter-by-parameter estimation method

Mechanistic models are powerful tools for chromatographic process development and optimization. However, hydrophobic interaction chromatography (HIC) mechanistic models lack an effective and logical parameter estimation method, especially for multi-component system. In this study, a parameter-by-parameter method for multi-component system (called as mPbP-HIC) was derived based on the retention mechanism to estimate the six parameters of the Mollerup isotherm for HIC. The linear parameters (ks,i and keq,i) and nonlinear parameters (ni and qmax,i) of the isotherm can be estimated by the linear regression (LR) and the linear approximation (LA) steps, respectively. The remaining two parameters (kp,i and kkin,i) are obtained by the inverse method (IM). The proposed method was verified with a two-component model system. The results showed that the model could accurately predict the protein elution at a loading of 10 g/L. However, the elution curve fitting was unsatisfactory for high loadings (12 g/L and 14 g/L), which is mainly attributed to the demanding experimental conditions of the LA step and the potential large estimation error of the parameter qmax. Therefore, the inverse method was introduced to further calibrate the parameter qmax, thereby reducing the estimation error and improving the curve fitting. Moreover, the simplified linear approximation (SLA) was proposed by reasonable assumption, which provides the initial guess of qmax without solving any complex matrix and avoids the problem of matrix unsolvable. In the improved mPbP-HIC method, qmax would be initialized by the SLA and finally determined by the inverse method, and this strategy was named as SLA+IM. The experimental validation showed that the improved mPbP-HIC method has a better curve fitting, and the use of SLA+IM reduces the error accumulation effect. In process optimization, the parameters estimated by the improved mPbP-HIC method provided the model with excellent predictive ability and reasonable extrapolation. In conclusion, the SLA+IM strategy makes the improved mPbP-HIC method more rational and can be easily applied to the practical separation of protein mixture, which would accelerate the process development for HIC in downstream of biopharmaceuticals.

Methodology for fast development of digital solutions in integrated continuous downstream processing

The methodology for production of biologics is going through a paradigm shift from batch-wise operation to continuous production. Lot of efforts are focused on integration, intensification, and continuous operation for decreased foot-print, material, equipment, and increased productivity and product quality. These integrated continuous processes with on-line analytics become complex processes, which requires automation, monitoring, and control of the operation, even unmanned or remote, which means bioprocesses with high level of automation or even autonomous capabilities. The development of these digital solutions becomes an important part of the process development and needs to be assessed early in the development chain. This work discusses a platform that allows fast development, advanced studies, and validation of digital solutions for integrated continuous downstream processes. It uses an open, flexible, and extendable real-time supervisory controller, called Orbit, developed in Python. Orbit makes it possible to communicate with a set of different physical setups and on the same time perform real-time execution. Integrated continuous processing often implies parallel operation of several setups and network of Orbit controllers makes it possible to synchronize complex process system. Data handling, storage, and analysis are important properties for handling heterogeneous and asynchronous data generated in complex downstream systems. Digital twin applications, such as advanced model-based and plant-wide monitoring and control, are exemplified using computational extensions in Orbit, exploiting data and models. Examples of novel digital solutions in integrated downstream processes are automatic operation parameter optimization, Kalman filter monitoring, and model-based batch-to-batch control.

Simultaneous prediction of 16 quality attributes during protein A chromatography using machine learning based Raman spectroscopy models

Several key technologies for advancing biopharmaceutical manufacturing depend on the successful implementation of process analytical technologies that can monitor multiple product quality attributes in a continuous in-line setting. Raman spectroscopy is an emerging technology in the biopharma industry that promises to fit this strategic need, yet its application is not widespread due to limited success for predicting a meaningful number of quality attributes. In this study, we addressed this very problem by demonstrating new capabilities for preprocessing Raman spectra using a series of Butterworth filters. The resulting increase in the number of spectral features is paired with a machine learning algorithm and laboratory automation hardware to drive the automated collection and training of a calibration model that allows for the prediction of 16 different product quality attributes in an in-line mode. The demonstrated ability to generate these Raman-based models for in-process product quality monitoring is the breakthrough to increase process understanding by delivering product quality data in a continuous manner. The implementation of this multiattribute in-line technology will create new workflows within process development, characterization, validation, and control.

Automatic procedure for modelling, calibration, and optimization of a three-component chromatographic separation

The current landscape of biopharmaceutical production necessitates an ever-growing set of tools to meet the demands for shorter development times and lower production costs. One path towards meeting these demands is the implementation of digital tools in the development stages. Mathematical modelling of process chromatography, one of the key unit operations in the biopharmaceutical downstream process, is one such tool. However, obtaining parameter values for such models is a time-consuming task that grows in complexity with the number of compounds in the mixture being purified. In this study, we tackle this issue by developing an automated model calibration procedure for purification of a multi-component mixture by linear gradient ion exchange chromatography. The procedure was implemented using the Orbit software (Lund University, Department of Chemical Engineering), which both generates a mathematical model structure and performs the experiments necessary to obtain data for model calibration. The procedure was extended to suggest operating points for the purification of one of the components in the mixture by means of multi-objective optimization using three different objectives. The procedure was tested on a three-component protein mixture and was able to generate a calibrated model capable of reproducing the experimental chromatograms to a satisfactory degree, using a total of six assays. An additional seventh experiment was performed to validate the model response under one of the suggested optimum conditions, respecting a 95 % purity requirement. All of the above was automated and set in motion by the push of a button. With these results, we have taken a step towards fully automating model calibration and thus accelerating digitalization in the development stages of new biopharmaceuticals.

Application of Raman spectroscopy during pharmaceutical process development for determination of critical quality attributes in Protein A chromatography

Raman spectroscopy is considered a Process Analytical Technology (PAT) tool in biopharmaceutical downstream processes. In the past decade, researchers have shown Raman spectroscopy's feasibility in determining Critical Quality Attributes (CQAs) in bioprocessing. This study verifies the feasibility of implementing a Raman-based PAT tool in Protein A chromatography as a CQA monitoring technique, for the purpose of accelerating process development and achieving real-time release in manufacturing. A system connecting Raman to a Tecan liquid handling station enables high-throughput model calibration. One calibration experiment collects Raman spectra of 183 samples with 8 CQAs within 25 h. After applying Butterworth high-pass filters and k-nearest neighbor (KNN) regression for model training, the model showed high predictive accuracy for fragments (Q2 = 0.965) and strong predictability for target protein concentration, aggregates, as well as charge variants (Q2≥ 0.922). The model's robustness was confirmed by varying the elution pH, load density, and residence time using 19 external validation preparative Protein A chromatography runs. The model can deliver elution profiles of multiple CQAs within a set point ± 0.3 pH range. The CQA readouts were presented as continuous chromatograms with a resolution of every 28 s for enhanced process understanding. In external validation datasets, the model maintained strong predictability especially for target protein concentration (Q2 = 0.956) and basic charge variants (Q2 = 0.943), except for overpredicted HCP (Q2 = 0.539). This study demonstrates a rapid, effective method for implementing Raman spectroscopy for in-line CQA monitoring in process development and biomanufacturing, eliminating the need for labor-intensive sample pooling and handling.

Parameter-by-parameter estimation method for adsorption isotherm in hydrophobic interaction chromatography

Hydrophobic interaction chromatography (HIC) is used as a critical polishing step in the downstream processing of biopharmaceuticals. Normally the process development of HIC is a cumbersome and time-consuming task, and the mechanical models can provide a powerful tool to characterize the process, assist process design and accelerate process development. However, the current estimation of model parameters relies on the inverse method, which lacks an efficient and logical parameter estimation strategy. In this study, a parameter-by-parameter (PbP) method based on the theoretical derivation and simplifying assumptions was proposed to estimate the Mollerup isotherm parameters for HIC. The method involves three key steps: (1) linear regression (LR) to estimate the salt-protein interaction parameter and the equilibrium constant; (2) linear approximation (LA) to estimate the stoichiometric parameter and the maximum binding capacity; and (3) inverse method to estimate the protein-protein interaction parameter and the kinetic coefficient. The results indicated that the LR step should be used for dilution condition (loading factor below 5%), while the LA step should be conducted when the isotherm is in the transition or nonlinear regions. Six numerical experiments were conducted to implement the PbP method. The results demonstrated that the PbP method developed allows for the systematic estimation of HIC parameters one-by-one, effectively reducing the number of parameters required for inverse method estimation from six to two. This helps prevent non-identifiability of structural parameters. The feasibility of the PbP-HIC method was further validated by real-world experiments. Moreover, the PbP method enhances the mechanistic understanding of adsorption behavior of HIC and shows a promising application to other stoichiometric displacement model-derived isotherms.

Mechanistic modeling of empty-full separation in recombinant adeno-associated virus production using anion-exchange membrane chromatography

Recombinant adeno-associated viral vectors (rAAVs) have become an industry-standard technology in the field of gene therapy, but there are still challenges to be addressed in their biomanufacturing. One of the biggest challenges is the removal of capsid species other than that which contains the gene of interest. In this work, we develop a mechanistic model for the removal of empty capsids-those that contain no genetic material-and enrichment of full rAAV using anion-exchange membrane chromatography. The mechanistic model was calibrated using linear gradient experiments, resulting in good agreement with the experimental data. The model was then applied to optimize the purification process through maximization of yield studying the impact of mobile phase salt concentration and pH, isocratic wash and elution length, flow rate, percent full (purity) requirement, loading density (challenge), and the use of single-step or two-step elution modes. A solution from the optimization with purity of 90% and recovery yield of 84% was selected and successfully validated, as the model could predict the recovery yield with remarkable fidelity and was able to find process conditions that led to significant enrichment. This is, to the best of our knowledge, the first case study of the application of de novo mechanistic modeling for the enrichment of full capsids in rAAV manufacturing, and it serves as demonstration of the potential of mechanistic modeling in rAAV process development.

Model based process optimization of an industrial chromatographic process for separation of lactoferrin from bovine milk

Model based process development using predictive mechanistic models is a powerful tool for in-silico downstream process development. It allows to obtain a thorough understanding of the process reducing experimental effort. While in pharma industry, mechanistic modeling becomes more common in the last years, it is rarely applied in food industry. This case study investigates risk ranking and possible optimization of the industrial process of purifying lactoferrin from bovine milk using SP Sepharose Big Beads with a resin particle diameter of 200 µm, based on a minimal number of lab-scale experiments combining traditional scale-down experiments with mechanistic modeling. Depending on the location and season, process water pH and the composition of raw milk can vary, posing a challenge for highly efficient process development. A predictive model based on the general rate model with steric mass action binding, extended for pH dependence, was calibrated to describe the elution behavior of lactoferrin and main impurities. The gained model was evaluated against changes in flow rate, step elution conditions, and higher loading and showed excellent agreement with the observed experimental data. The model was then used to investigate the critical process parameters, such as water pH, conductivity of elution steps, and flow rate, on process performance and purity. It was found that the elution behavior of lactoferrin is relatively consistent over the pH range of 5.5 to 7.6, while the elution behavior of the main impurities varies greatly with elution pH. As a result, a significant loss in lactoferrin is unavoidable to achieve desired purities at pH levels below pH 6.0. Optimal process parameters were identified to reduce water and salt consumption and increase purity, depending on water pH and raw milk composition. The optimal conductivity for impurity removal in a low conductivity elution step was found to be 43 mS/cm, while a conductivity of 95 mS/cm leads to the lowest overall salt usage during lactoferrin elution. Further increasing the conductivity during lactoferrin elution can only slightly lower the elution volume thus can also lead to higher total salt usage. Low flow rates during elution of 0.2 column volume per minute are beneficial compared to higher flow rates of 1 column volume per minute. The, on lab-scale, calibrated model allows predicting elution volume and impurity removal for large-scale experiments in a commercial plant processing over 106 liters of milk per day. The successful model extrapolation was possible without recalibration or detailed knowledge of the manufacturing plant. This study therefore provides a possible pathway for rapid process development of chromatographic purification in the food industries combining traditional scale-down experiments with mechanistic modeling.

Understanding adsorption behavior of antiviral labyrinthopeptin peptides in anion exchange chromatography

Lantipeptides from bacterial sources are increasingly important as biopharmaceuticals because of their broad range of applications. However, the availability of most lantipeptides is low, and systematic approaches for downstream processing of this group of peptides is still lacking. Model-based development for chromatographic separations has proven to be a useful tool for developing reliable purification processes. One important compound of such a model is the adsorption behavior of the components of interest. In ion-exchange chromatography, the adsorption equilibrium between salt and proteins can be described using the steric mass action (SMA) formalism. Beyond, the model parameters may be related to the lanthipeptides physico-chemical properties. In this study, the antiviral lantipeptides labyrinthopeptin A1 and A2, purified from Actinomadura namibiensis culture broth, were characterized for their adsorption behavior in anion-exchange chromatography in the range from pH 5.0-7.4. The experiments necessary to determine the three SMA parameters were chosen in a way to limit the amount of peptides needed. Linear gradient elution was applied successfully to separate A1 and A2 and to determine the characteristic charge ν_i and the equilibrium constant [Formula: see text] . Batch adsorption experiments using a robotic workstation for high throughput and accuracy provided non-linear adsorption isotherms and the steric factor σ_i. Labyrinthopeptin A1 and A2 show a very different adsorption behavior even though the fundamental structure of the two peptides is similar. k_eq of A1 ranging from 0.18 to 0.88 are approximately one order of magnitude smaller than that of A2 ranging from 3.44 to 9.73 indicating the higher affinity of A2 to the stationary phase. At pH 7.0 σ was 1.12 and 0.60 for A1 and A2, respectively which was expected based on the molecular weight of the peptides. The characteristic charge for both peptides was also theoretically estimated from the amino acids involved in electrostatic interactions which was in good agreement with experimental data. Thereby, this work provides an useful approach to estimate SMA parameters based on simple structural information that can be applied early in chromatographic ion-exchange process development for peptides and may help adapting the processes for future designed lanthipeptides.

Standardized method for mechanistic modeling of multimodal anion exchange chromatography in flow through operation

Multimodal chromatography offers an increased selectivity compared to unimodal chromatographic methods and is often employed for challenging separation tasks in industrial downstream processing (DSP). Unfortunately, the implementation of multimodal polishing into a generic downstream platform can be hampered by non-robust platform conditions leading to a time and cost intensive process development. Mechanistic modeling can assist experimental process development but readily applicable and easy to calibrate multimodal chromatography models are lacking. In this work, we present a mechanistic modeling aided approach that paves the way for an accelerated development of anionic mixed-mode chromatography (MMC) for biopharmaceutical purification. A modified multimodal isotherm model was calibrated using only three chromatographic experiments and was employed in the retention prediction of four antibody formats including a Fab, a bispecific, as well as an IgG1 and IgG4 antibody subtype at pH 5.0 and 6.0. The chromatographic experiments were conducted using the anionic mixed-mode resin Capto adhere at industrial relevant process conditions to enable flow through purification. An existing multimodal isotherm model was reduced to hydrophobic interactions in the linear range of the adsorption isotherm and successfully employed in the simulation of six chromatographic experiments per molecule in concert with the transport dispersive model (TDM). The model reduction to only three parameters did prevent structural parameter non-identifiability and enabled an analytical isotherm parameter determination that was further refined by incorporation of size exclusion effects of the selected multimodal resin. During the model calibration, three linear salt gradient elution experiments were performed for each molecule followed by an isotherm parameter uncertainty assessment. Lastly, each model was validated with a set of step and isocratic elution experiments. This standardized modeling approach facilitates the implementation of multimodal chromatography as a key unit operation for the biopharmaceutical downstream platform, while increasing the mechanistic insight to the multimodal adsorption behavior of complex biologics.

Parameter-by-parameter method for steric mass action model of ion exchange chromatography: Simplified estimation for steric shielding factor

Mechanistic models play a crucial role in the process development and optimization of ion-exchange chromatography (IEC). Recent researches in steric mass action (SMA) model have heightened the need for better estimation of nonlinear parameter, steric shielding factor σ. In this work, a straightforward approach combination of simplified linear approximation (SLA) and inverse method (IM) was proposed to initialize and further determine σ, respectively. An existed, unique, and positive σ can be derived from SLA. Compared with linear approximation (LA) developed in our previous study, σ of the multi-component system can be calculated easily without solving the complex system of linear equations, leading to a time complexity reduction from O(n³) to O(n). The proposed method was verified first in numerical experiments about the separation of three charge variants. The calculated σ was more reasonable than that of LA, and the error of elution profiles with the parameters estimated by SLA+IM was only one-sixth of that by LA in numerical experiments. Moreover, the error accumulation effect could also be reduced. The proposed method was further confirmed in real-world experiments about the separation of monomer-dimer mixtures of monoclonal antibody. The results gave a lower error and better physical understanding compared to LA. In conclusion, SLA+IM developed in the present work provides a novel and straightforward way to determine σ. This simplification would help to save the effort of calibration experiments and accelerate the process development for the multi-component IEC separation.

Mechanistic modeling, simulation, and optimization of mixed-mode chromatography for an antibody polishing step

Mixed-mode chromatography combines features of ion-exchange chromatography and hydrophobic interaction chromatography and is increasingly used in antibody purification. As a replacement for flow-through operations on traditional unmixed resins or as a pH-controlled bind-and-elute step, the use of both interaction modes promises a better removal of product-specific impurities. However, the combination of the functionalities makes industrial process development significantly more complex, in particular the identification of the often small elution window that delivers the desired selectivity. Mechanistic modeling has proven that even difficult separation problems can be solved in a computer-optimized manner once the process dynamics have been modeled. The adsorption models described in the literature are also very complex, which makes model calibration difficult. In this work, we approach this problem with a newly constructed model that describes the adsorber saturation with the help of the surface coverage function of the colloidal particle adsorption model for ion-exchange chromatography. In a case study, a model for a pH-controlled antibody polishing step was created from six experiments. The behavior of fragments, aggregates, and host cell proteins was described with the help of offline analysis. After in silico optimization, a validation experiment confirmed an improved process performance in comparison to the historical process set point. In addition to these good results, the work also shows that the high dynamics of mixed-mode chromatography can produce unexpected results if process parameters deviate too far from tried and tested conditions.

The use of predictive models to develop chromatography-based purification processes

Chromatography is the workhorse of biopharmaceutical downstream processing because it can selectively enrich a target product while removing impurities from complex feed streams. This is achieved by exploiting differences in molecular properties, such as size, charge and hydrophobicity (alone or in different combinations). Accordingly, many parameters must be tested during process development in order to maximize product purity and recovery, including resin and ligand types, conductivity, pH, gradient profiles, and the sequence of separation operations. The number of possible experimental conditions quickly becomes unmanageable. Although the range of suitable conditions can be narrowed based on experience, the time and cost of the work remain high even when using high-throughput laboratory automation. In contrast, chromatography modeling using inexpensive, parallelized computer hardware can provide expert knowledge, predicting conditions that achieve high purity and efficient recovery. The prediction of suitable conditions in silico reduces the number of empirical tests required and provides in-depth process understanding, which is recommended by regulatory authorities. In this article, we discuss the benefits and specific challenges of chromatography modeling. We describe the experimental characterization of chromatography devices and settings prior to modeling, such as the determination of column porosity. We also consider the challenges that must be overcome when models are set up and calibrated, including the cross-validation and verification of data-driven and hybrid (combined data-driven and mechanistic) models. This review will therefore support researchers intending to establish a chromatography modeling workflow in their laboratory.

An accelerated approach for mechanistic model based prediction of linear gradient elution ion-exchange chromatography of proteins

With growing demands for therapeutic monoclonal antibodies, in silico downstream process development based on mechanistic modeling of chromatography separation process is being increasingly used for process optimization and process characterization. Application of mechanistic modeling in biopharmaceutical industry has been sparse due to the significant investment of time and resources that are required for performing model calibration. Mechanistic modeling of the chromatography process involves a large number of mass transport and binding parameters and their initial input values are required for simulations. These input values of column parameters can be easily obtained either from experiments or from empirical correlations available in literature. On the other hand, obtaining the model input valves for binding kinetic parameters is usually a cumbersome process as it involves performing batch experiments which are not only tedious but also require significant quantities of purely isolated main product and its related impurities, which is challenging as the product related impurities are typically present in smaller quantities and hence are difficult to obtain as pure species. In the present work, a mechanistic model that is based on the general rate model coupled with extended Langmuir binding model has been used for prediction of linear gradient elution peaks of monoclonal antibody on cation exchanger chromatography. The present work describes an accelerated approach for obtaining the input values for binding kinetic parameters in the extended Langmuir binding model from the two Yamamoto coefficient A and B values obtained by Yamamoto method directly from the model calibration linear gradient elution runs of different gradient slopes and at low to moderate protein loadings. The equations that can relate the two coefficients to the extended Langmuir model equation binding kinetic parameters were derived. Therefore, once A and B are determined, the binding kinetic parameter values were determined straightforward, thereby avoiding the problem of multiple solutions for the model parameters. The estimated binding parameters were successfully validated from isocratic elution experiments performed at low loading. What we demonstrate is that the proposed approach allows us to estimate binding kinetic parameters in a significantly more efficient and accelerated manner than presently used approaches, thereby accelerating development and implementation of mechanistic modeling for process chromatography.

Parameter-by-parameter method for steric mass action model of ion exchange chromatography: Theoretical considerations and experimental verification

Ion exchange chromatography (IEC) is one of the most widely-used techniques for protein separation and has been characterized by mechanistic models. However, the time-consuming and cumbersome model calibration hinders the application of mechanistic models for process development. A new methodology called "parameter-by-parameter method (PbP)" was proposed with mechanistic derivations of the steric mass action (SMA) model of IEC. The protocol includes four steps: (1) first linear regression (LR1) for characteristic charge; (2) second linear regression (LR2) for equilibrium coefficient; (3) linear approximation (LA) for shielding factor; (4) inverse method (IM) for kinetic coefficient. Four SMA parameters could be one-by-one determined in sequence, reducing the number of unknown parameters per species from four to one, and predicting almost consistent retention. Numerical single-component experiments were investigated firstly, and the PbP method showed excellent agreement between experiments and simulations. The effects of loadings on the PbP and Yamamoto methods were compared. It was found that the PbP method had higher accuracy and robustness than the Yamamoto method. Moreover, a five-experiment strategy was suggested to implement the PbP method, which is straightforward to reduce the cost of calibration experiments. Finally, a real-world multi-component separation was challenged and further confirmed the feasibility of the PbP method. In general, the proposed method can not only reliably estimate the SMA parameters with comprehensive physical understanding but also accurately predict retention over a wide range of loading conditions.

Automation of Modeling and Calibration of Integrated Preparative Protein Chromatography Systems

With the increasing global demand for precise and efficient pharmaceuticals and the biopharma industry moving towards Industry 4.0, the need for advanced process integration, automation, and modeling has increased as well. In this work, a method for automatic modeling and calibration of an integrated preparative chromatographic system for pharmaceutical development and production is presented. Based on a user-defined system description, a system model was automatically generated and then calibrated using a sequence of experiments. The system description and model was implemented in the Python-based preparative chromatography control software Orbit.

Holistic Process Models: A Bayesian Predictive Ensemble Method for Single and Coupled Unit Operation Models

The coupling of individual models in terms of end-to-end calculations for unit operations in manufacturing processes is a challenging task. We present a probability distribution-based approach for the combined outcomes of parametric and non-parametric models. With this so-called Bayesian predictive ensemble, the statistical moments such as mean value and standard deviation can be accurately computed without any further approximation. It is shown that the ensemble of different model predictions leads to an uninformed prior distribution, which can be transformed into a predictive posterior distribution using Bayesian inference and numerical Markov Chain Monte Carlo calculations. We demonstrate the advantages of our method using several numerical examples. Our approach is not restricted to certain unit operations, and can also be used for the more robust interpretation and assessment of model predictions in general.

Mechanistic modeling and simulation of a complex low and high loading elution behavior of a polypeptide in cation exchange chromatography

The mechanistic modeling of preparative liquid chromatography is still a challenging task. Non-ideal thermodynamic conditions may require activity coefficients for the mechanistic description of preparative chromatography. In this work, a chromatographic cation exchange step with a polypeptide having a complex elution behavior in low and high loading situations is modeled. Model calibration in the linear range of the isotherm is done by applying counterion-induced linear gradient elution experiments between pH 3.3 and pH 4.3. Inverse fitting with column loads up to 25 mg/mL_CV is performed for parameter estimation in the non-linear range. The polypeptide elution peak shows an anti-Langmuirian behavior with fronting under low loading conditions and a switch to a Langmuirian behavior with increasing load. This unusual elution behavior could be described with an extended version of the sigmoidal Self-Association isotherm, including two activity coefficients for the polypeptide and counterion in solution. The activity coefficient of the solute polypeptide shows a strong influence on the model parameters and is crucial in the linear and non-linear range of the isotherm. The modeling procedure results in a unique and robust model parameter set that is sufficient to describe the complex elution behavior and allows modeling over the full isotherm range. This article is protected by copyright. All rights reserved.

Behavior of weakly adsorbing protein impurities in flow-through ion-exchange chromatography

Flow-through ion-exchange chromatography is frequently used in polishing biotherapeutics, but the factors that contribute to impurity persistence are incompletely understood. A large number of dilute impurities may be encountered that exhibit physicochemical diversity, making the flow-through separation performance highly sensitive to process conditions. The analysis presented in this work develops two novel correlations that offer transferable insights into the chromatographic behavior of weakly adsorbing impurities. The first, based on column simulations and validated experimentally, delineates the relative contributions of thermodynamic, transport, and geometric properties in dictating the initial breakthrough volumes of dilute species. The Graetz number for mass transfer was found to generalize the transport contributions, enabling estimation of a threshold in the equilibrium constant below which impurity persistence is expected. Impurity adsorption equilibria are needed to use this correlation, but such data are not typically available. The second relationship presented in this work may be used to reduce the experimental burden of estimating adsorption equilibria as a function of ionic strength. A correlation between stoichiometric displacement model parameters was found by consolidating isocratic retention data for over 200 protein-pH-resin combinations from the extant literature. Coupled with Yamamoto's analysis of linear gradient elution data, this correlation may be used to estimate retentivity approximately from a single experimental measurement, which could prove useful in predicting host-cell protein chromatographic behavior.

Stress-dependent flexibility of a full-length human monoclonal antibody: Insights from molecular dynamics to support biopharmaceutical development

After several decades of advancements in drug discovery, product development of biopharmaceuticals remains a time- and resource-consuming endeavor. One of the main reasons is associated to the lack of fundamental understanding of conformational dynamics of such biologic entities, and how they respond to various stresses encountered during manufacturing. In this work, we have studied the conformational dynamics of human IgG₁κ b12 monoclonal antibody (mAb) using molecular dynamics simulations. The hundreds of nanoseconds long trajectories reveal that b12 mAb is highly flexible. Its variable domains show greater conformational fluctuations than the constant domains. Additionally, it collapses towards a more globular shape in response to thermal stress, leading to decrease in the total solvent exposed surface area and radius of gyration. This behavior is more pronounced for the deglycosylated b12 mAb, and it appears to correlate with increase in inter-domain contacts between specific regions of the antibody. Conformational fluctuations also cause temporary formation and disruption of hydrophobic and charged patches on the antibody surface, which is particularly important for the prediction of CMC properties during development phases of antibody-based biotherapeutics. The insights gained through these simulations may help the development of biologic drugs, especially with regards to manufacturing processes where antibodies may undergo significant thermal stress.

Generic and specific recurrent neural network models: Applications for large and small scale biopharmaceutical upstream processes

The calculation of temporally varying upstream process outcomes is a challenging task. Over the last years, several parametric, semi-parametric as well as non-parametric approaches were developed to provide reliable estimates for key process parameters. We present generic and product-specific recurrent neural network (RNN) models for the computation and study of growth and metabolite-related upstream process parameters as well as their temporal evolution. Our approach can be used for the control and study of single product-specific large-scale manufacturing runs as well as generic small-scale evaluations for combined processes and products at development stage. The computational results for the product titer as well as various major upstream outcomes in addition to relevant process parameters show a high degree of accuracy when compared to experimental data and, accordingly, a reasonable predictive capability of the RNN models. The calculated values for the root-mean squared errors of prediction are significantly smaller than the experimental standard deviation for the considered process run ensembles, which highlights the broad applicability of our approach. As a specific benefit for platform processes, the generic RNN model is also used to simulate process outcomes for different temperatures in good agreement with experimental results. The high level of accuracy and the straightforward usage of the approach without sophisticated parameterization and recalibration procedures highlight the benefits of the RNN models, which can be regarded as promising alternatives to existing parametric and semi-parametric methods.

In silico process characterization for biopharmaceutical development following the quality by design concept

With the quality by design (QbD) initiative, regulatory authorities demand a consistent drug quality originating from a well-understood manufacturing process. This study demonstrates the application of a previously published mechanistic chromatography model to the in silico process characterization (PCS) of a monoclonal antibody polishing step. The proposed modeling workflow covered the main tasks of traditional PCS studies following the QbD principles, including criticality assessment of 11 process parameters and establishment of their proven acceptable ranges of operation. Analyzing effects of multi-variate sampling of process parameters on the purification outcome allowed identification of the edge-of-failure. Experimental validation of in silico results demanded approximately 75% less experiments compared to a purely wet-lab based PCS study. Stochastic simulation, considering the measured variances of process parameters and loading material composition, was used to estimate the capability of the process to meet the acceptance criteria for critical quality attributes and key performance indicators. The proposed workflow enables the implementation of digital process twins as QbD tool for improved development of biopharmaceutical manufacturing processes.

Analysis of complex protein elution behavior in preparative ion exchange processes using a colloidal particle adsorption model

A fundamental understanding of the protein retention mechanism in preparative ion exchange (IEX) chromatography columns is essential for a model-based process development approach. For the past three decades, the mechanistic description of protein retention has been based predominantly on the steric mass action (SMA) model. In recent years, however, retention profiles of proteins have been reported more frequently for preparative processes that are not consistent with the mechanistic understanding relying on the SMA model. In this work, complex elution behavior of proteins in preparative IEX processes is analyzed using a colloidal particle adsorption (CPA) model. The CPA model is found to be capable of reproducing elution profiles that cannot be described by the traditional SMA model. According to the CPA model, the reported complex behavior can be ascribed to a strong compression and concentration of the elution front in the lower unsaturated part of the chromatography column. As the unsaturated part of the column decreases with increasing protein load density, exceeding a critical load density can lead to the formation of a shoulder in the peak front. The general applicability of the model in describing preparative IEX processes is demonstrated using several industrial case studies including multiple monoclonal antibodies on different IEX adsorber systems. In this context, the work covers both salt controlled and pH-controlled protein elution.

Protein adsorption on ion exchange adsorbers: A comparison of a stoichiometric and non-stoichiometric modeling approach

For mechanistic modeling of ion exchange (IEX) processes, a profound understanding of the adsorption mechanism is important. While the description of protein adsorption in IEX processes has been dominated by stoichiometric models like the steric mass action (SMA) model, discrepancies between experimental data and model results suggest that the conceptually simple stoichiometric description of protein adsorption provides not always an accurate representation of nonlinear adsorption behavior. In this work an alternative colloidal particle adsorption (CPA) model is introduced. Based on the colloidal nature of proteins, the CPA model provides a non-stoichiometric description of electrostatic interactions within IEX columns. Steric hindrance at the adsorber surface is considered by hard-body interactions between proteins using the scaled-particle theory. The model's capability of describing nonlinear protein adsorption is demonstrated by simulating adsorption isotherms of a monoclonal antibody (mAb) over a wide range of ionic strength and pH. A comparison of the CPA model with the SMA model shows comparable model results in the linear adsorption range, but significant differences in the nonlinear adsorption range due to the different mechanistic interpretation of steric hindrance in both models. The results suggest that nonlinear adsorption effects can be overestimated by the stoichiometric formalism of the SMA model and are generally better reproduced by the CPA model.

About Model Validation in Bioprocessing

In bioprocess engineering the Qualtiy by Design (QbD) initiative encourages the use of models to define design spaces. However, clear guidelines on how models for QbD are validated are still missing. In this review we provide a comprehensive overview of the validation methods, mathematical approaches, and metrics currently applied in bioprocess modeling. The methods cover analytics for data used for modeling, model training and selection, measures for predictiveness, and model uncertainties. We point out the general issues in model validation and calibration for different types of models and put this into the context of existing health authority recommendations. This review provides a starting point for developing a guide for model validation approaches. There is no one-fits-all approach, but this review should help to identify the best fitting validation method, or combination of methods, for the specific task and the type of bioprocess model that is being developed.

Validation of Novel Lattice Boltzmann Large Eddy Simulations (LB LES) for Equipment Characterization in Biopharma

Detailed process and equipment knowledge is crucial for the successful production of biopharmaceuticals. An essential part is the characterization of equipment for which Computational Fluid Dynamics (CFD) is an important tool. While the steady, Reynolds-averaged Navier–Stokes (RANS) k − ε approach has been extensively reviewed in the literature and may be used for fast equipment characterization in terms of power number determination, transient schemes have to be further investigated and validated to gain more detailed insights into flow patterns because they are the method of choice for mixing time simulations. Due to the availability of commercial solvers, such as M-Star CFD, Lattice Boltzmann simulations have recently become popular in the industry, as they are easy to set up and require relatively low computing power. However, extensive validation studies for transient Lattice Boltzmann Large Eddy Simulations (LB LES) are still missing. In this study, transient LB LES were applied to simulate a 3 L bioreactor system. The results were compared to novel 4D particle tracking (4D PTV) experiments, which resolve the motion of thousands of passive tracer particles on their journey through the bioreactor. Steady simulations for the determination of the power number followed a structured workflow, including grid studies and rotating reference frame volume studies, resulting in high prediction accuracy with less than 11% deviation, compared to experimental data. Likewise, deviations for the transient simulations were less than 10% after computational demand was reduced as a result of prior grid studies. The time averaged flow fields from LB LES were in good accordance with the novel 4D PTV data. Moreover, 4D PTV data enabled the validation of transient flow structures by analyzing Lagrangian particle trajectories. This enables a more detailed determination of mixing times and mass transfer as well as local exposure times of local velocity and shear stress peaks. For the purpose of standardization of common industry CFD models, steady RANS simulations for the 3 L vessel were included in this study as well.

Modeling the impact of amino acid substitution in a monoclonal antibody on cation exchange chromatography

A vital part of biopharmaceutical research is decision making around which lead candidate should be progressed in early-phase development. When multiple antibody candidates show similar biological activity, developability aspects are taken into account to ease the challenges of manufacturing the potential drug candidate. While current strategies for developability assessment mainly focus on drug product stability, only limited information is available on how antibody candidates with minimal differences in their primary structure behave during downstream processing. With increasing time-to-market pressure and an abundance of monoclonal antibodies (mAbs) in development pipelines, developability assessments should also consider the ability of mAbs to integrate into the downstream platform. This study investigates the influence of amino acid substitutions in the complementarity-determining region (CDR) of a full-length IgG1 mAb on the elution behavior in preparative cation exchange chromatography. Single amino acid substitutions within the investigated mAb resulted in an additional positive charge in the light chain (L) and heavy chain (H) CDR, respectively. The mAb variants showed an increased retention volume in linear gradient elution compared with the wild-type antibody. Furthermore, the substitution of tryptophan with lysine in the H-CDR3 increased charge heterogeneity of the product. A multiscale in silico analysis, consisting of homology modeling, protein surface analysis, and mechanistic chromatography modeling increased understanding of the adsorption mechanism. The results reveal the potential effects of lead optimization during antibody drug discovery on downstream processing.

Cross-scale quality assessment of a mechanistic cation exchange chromatography model

Cation exchange chromatography (CEX) is an essential part of most monoclonal antibody (mAb) purification platforms. Process characterization and root cause investigation of chromatographic unit operations are performed using scale down models (SDM). SDM chromatography columns typically have the identical bed height as the respective manufacturing-scale, but a significantly reduced inner diameter. While SDMs enable process development demanding less material and time, their comparability to manufacturing-scale can be affected by variability in feed composition, mobile phase and resin properties, or dispersion effects depending on the chromatography system at hand. Mechanistic models can help to close gaps between scales and reduce experimental efforts compared to experimental SDM applications. In this study, a multicomponent steric mass-action (SMA) adsorption model was applied to the scale-up of a CEX polishing step. Based on chromatograms and elution pool data ranging from laboratory- to manufacturing-scale, the proposed modeling workflow enabled early identification of differences between scales, for example, system dispersion effects or ionic capacity variability. A multistage model qualification approach was introduced to measure the model quality and to understand the model's limitations across scales. The experimental SDM and the in silico model were qualified against large-scale data using the identical state of the art equivalence testing procedure. The mechanistic chromatography model avoided limitations of the SDM by capturing effects of bed height, loading density, feed composition, and mobile phase properties. The results demonstrate the applicability of mechanistic chromatography models as a possible alternative to conventional SDM approaches.