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

High-throughput in silico workflow for optimization and characterization of multimodal chromatographic processes

While high-throughput (HT) experimentation and mechanistic modeling have long been employed in chromatographic process development, it remains unclear how these techniques should be used in concert within development workflows. In this work, a process development workflow based on HT experiments and mechanistic modeling was constructed. The integration of HT and modeling approaches offers improved workflow efficiency and speed. This high-throughput in silico (HT-IS) workflow was employed to develop a Capto MMC polishing step for mAb aggregate removal. High-throughput batch isotherm data was first generated over a range of mobile phase conditions and a suite of analytics were employed. Parameters for the extended steric mass action (SMA) isotherm were regressed for the multicomponent system. Model validation was performed using the extended SMA isotherm in concert with the general rate model of chromatography using the CADET modeling software. Here, step elution profiles were predicted for eight RoboColumn runs across a range of ionic strength, pH, and load density. Optimized processes were generated through minimization of a complex objective function based on key process metrics. Processes were evaluated at lab-scale using two feedstocks, differing in composition. The results confirmed that both processes obtained high monomer yield (>85%) and removed∼ 50 % $$ \sim 50\% $$ of aggregate species. Column simulations were then carried out to determine sensitivity to a wide range of process inputs. Elution buffer pH was found to be the most critical process parameter, followed by resin ionic capacity. Overall, this study demonstrated the utility of the HT-IS workflow for rapid process development and characterization.

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.

Automated calibration and in-line measurement of product quality during therapeutic monoclonal antibody purification using Raman spectroscopy

Current manufacturing and development processes for therapeutic monoclonal antibodies demand increasing volumes of analytical testing for both real-time process controls and high-throughput process development. The feasibility of using Raman spectroscopy as an in-line product quality measuring tool has been recently demonstrated and promises to relieve this analytical bottleneck. Here, we resolve time-consuming calibration process that requires fractionation and preparative experiments covering variations of product quality attributes (PQAs) by engineering an automation system capable of collecting Raman spectra on the order of hundreds of calibration points from two to three stock seed solutions differing in protein concentration and aggregate level using controlled mixing. We used this automated system to calibrate multi-PQA models that accurately measured product concentration and aggregation every 9.3 s using an in-line flow-cell. We demonstrate the application of a nonlinear calibration model for monitoring product quality in real-time during a biopharmaceutical purification process intended for clinical and commercial manufacturing. These results demonstrate potential feasibility to implement quality monitoring during GGMP manufacturing as well as to increase chemistry, manufacturing, and controls understanding during process development, ultimately leading to more robust and controlled manufacturing processes.

Model-assisted process development, characterization and design of continuous chromatography for antibody separation

Continuous manufacturing in monoclonal antibody production has generated increased interest due to its consistent quality, high productivity, high equipment utilization, and low cost. One of the major challenges in realizing continuous biological manufacturing lies in implementing continuous chromatography. Given the complex operation mode and various operation parameters, it is challenging to develop a continuous process. Due to the process parameters being mainly determined by the breakthrough curves and elution behaviors, chromatographic modeling has gradually been used to assist in process development and characterization. Model-assisted approaches could realize multi-parameter interaction investigation and multi-objective optimization by integrating continuous process models. These approaches could reduce time and resource consumption while achieving a comprehensive and systematic understanding of the process. This paper reviews the application of modeling tools in continuous chromatography process development, characterization and design. Model-assisted process development approaches for continuous capture and polishing steps are introduced and summarized. The challenges and potential of model-assisted process characterization are discussed, emphasizing the need for further research on the design space determination strategy and parameter robustness analysis method. Additionally, some model applications for process design were highlighted to promote the establishment of the process optimization and process simulation platform.

Isotherm model discrimination for multimodal chromatography using mechanistic models derived from high-throughput batch isotherm data

In this work, we have examined an array of isotherm formalisms and characterized them based on their relative complexities and predictive abilities with multimodal chromatography. The set of isotherm models studied were all based on the stoichiometric displacement framework, with considerations for electrostatic interactions, hydrophobic interactions, and thermodynamic activities. Isotherm parameters for each model were first determined through twenty repeated fits to a set of mAb - Capto MMC batch isotherm data spanning a range of loading, ionic strength, and pH as well as a set of mAb - Capto Adhere batch data at constant pH. The batch isotherm data were used in two ways-spanning the full range of loading or consisting of only the high concentration data points. Predictive ability was defined through the model's capacity to capture prominent changes in salt gradient elution behavior with respect to pH for Capto MMC or unique elution patterns and yield losses with respect to gradient slope for Capto Adhere. In both cases, model performance was quantified using a scoring metric based on agreement in peak characteristics for column predictions and accuracy of fit for the batch data. These scores were evaluated for all twenty isotherm fits and their corresponding column predictions, thereby producing a statistical distribution of model performances. Model complexity (number of isotherm parameters) was then considered through use of the Akaike information criterion (AIC) calculated from the score distributions. While model performance for Capto MMC benefitted substantially from removal of low protein concentration data, this was not the case for Capto Adhere; this difference was likely due to the qualitatively different shapes of the isotherms between the two resins. Surprisingly, the top-performing (high accuracy with minimal number of parameters) isotherm model was the same for both resins. The extended steric mass action (SMA) isotherm (containing both protein-salt and protein-protein activity terms) accurately captured both the pH-dependent elution behavior for Capto MMC as well as loss in protein recovery with increasing gradient slope for Capto Adhere. In addition, this isotherm model achieved the highest median score in both resin systems, despite it lacking any explicit hydrophobic stoichiometric terms. The more complex isotherm models, which explicitly accounted for both electrostatic and hydrophobic interaction stoichiometries, were ill-suited for Capto MMC and had lower AIC model likelihoods for Capto Adhere due to their increased complexity. Interestingly, the ability of the extended SMA isotherm to predict the Capto Adhere results was largely due to the protein-salt activity coefficient, as determined via isotherm parameter sensitivity analyses. Further, parametric studies on this parameter demonstrated that it had a major impact on both binding affinity and elution behavior, therein fully capturing the impact of hydrophobic interactions. In summary, we were able to determine the isotherm formalisms most capable of consistently predicting a wide range of column behavior for both a multimodal cation-exchange and multimodal anion-exchange resin with high accuracy, while containing a minimized set of model parameters.

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.

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.