Combined Yamamoto approach for simultaneous estimation of adsorption isotherm and kinetic parameters in ion-exchange chromatography

Realization of digital twin for dynamic control toward sample variation of ion exchange chromatography in antibody separation

Digital twin (DT) is a virtual and digital representation of physical objects or processes. In this paper, this concept is applied to dynamic control of the collection window in the ion exchange chromatography (IEC) toward sample variations. A possible structure of a feedforward model-based control DT system was proposed. Initially, a precise IEC mechanistic model was established through experiments, model fitting, and validation. The average root mean square error (RMSE) of fitting and validation was 8.1% and 7.4%, respectively. Then a model-based gradient optimization was performed, resulting in a 70.0% yield with a remarkable 11.2% increase. Subsequently, the DT was established by systematically integrating the model, chromatography system, online high-performance liquid chromatography, and a server computer. The DT was validated under varying load conditions. The results demonstrated that the DT could offer an accurate control with acidic variants proportion and yield difference of less than 2% compared to the offline analysis. The embedding mechanistic model also showed a positive predictive performance with an average RMSE of 11.7% during the DT test under >10% sample variation. Practical scenario tests indicated that tightening the control target could further enhance the DT robustness, achieving over 98% success rate with an average yield of 72.7%. The results demonstrated that the constructed DT could accurately mimic real-world situations and perform an automated and flexible pooling in IEC. Additionally, a detailed methodology for applying DT was summarized.

Mechanistic model-based characterization of size-exclusion-mixed-mode resins for removal of monoclonal antibody fragments

Although antibody fragments are a critical impurity to remove from process streams, few platformable purification techniques have been developed to this end. In this work, a novel size-exclusion-mixed-mode (SEMM) resin was characterized with respect to its efficacy in mAb fragment removal. Inverse size-exclusion chromatography showed that the silica-based resin had a narrow pore size distribution and a median pore radius of roughly 6.2 nm. Model-based characterization was carried out with Chromatography Analysis and Design Toolkit (CADET), using the general rate model and the multicomponent Langmuir isotherm. Model parameters were obtained from fitting breakthrough curves, performed at multiple residence times, for a mixture of mAb, aggregates, and an array of fragments (varying in size). Accurate fits were obtained to the frontal chromatographic data across a range of residence times. Model validation was then performed with a scaled-up column, altering residence time and feed composition from the calibration run. Accurate predictions were obtained, thereby illustrating the model's interpolative and extrapolative capabilities. Additionally, the SEMM resin achieved 90% mAb yield, 37% aggregate removal, 29% [Formula: see text] removal, 54% Fab/Fc removal, 100% Fc fragments removal, and a productivity of 72.3 g mAbL×h. Model predictions for these statistics were all within 5%. Simulated batch uptake experiments showed that resin penetration depth was directly related to protein size, with the exception of the aggregate species, and that separation was governed by differential pore diffusion rates. Additional simulations were performed to characterize the dependence of fragment removal on column dimension, load density, and feed composition. Fragment removal was found to be highly dependent on column load density, where optimal purification was achieved below 100 mg protein/mL column. Furthermore, fragment removal was dependent on column volume (constant load mass), but agnostic to whether column length or diameter was changed. Lastly, the dependence on feed composition was shown to be complex. While fragment removal was inversely related to fragment mass fraction in the feed, the extent depended on fragment size. Overall, the results from this study illustrated the efficacy of the SEMM resin in fragment and aggregate removal and elucidated relationships with key operational parameters through model-based characterization.

Application of mechanistic modelling in membrane and fiber chromatography for purification of biotherapeutics - A review

Mechanistic modelling is a simulation tool which has been effectively applied in downstream bioprocessing to model resin chromatography. Membrane and fiber chromatography are newer approaches that offer higher rates of mass transfer and consequently higher flow rates and reduced processing times. This review describes the key considerations in the development of mechanistic models for these unit operations. Mass transfer is less complex than in resin columns, but internal housing volumes can make modelling difficult, particularly for laboratory-scale devices. Flow paths are often non-linear and the dead volume is often a larger fraction of the overall volume, which may require more complex hydrodynamic models to capture residence time distributions accurately. In this respect, the combination of computational fluid dynamics with appropriate protein binding models is emerging as an ideal approach.

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 cation exchange chromatography scale-up considering packing inhomogeneities

In process development and characterization, the scale-up of chromatographic steps is a crucial part and brings a number of challenges. Usually, scale-down models are used to represent the process step, and constant column properties are assumed. The scaling is then typically based on the concept of linear scale-up. In this work, a mechanistic model describing an anti-Langmuirian to Langmuirian elution behavior of a polypeptide, calibrated with a pre-packed 1 ml column, is applied to demonstrate the scalability to larger column volumes up to 28.2 ml. Using individual column parameters for each column size, scaling to similar eluting salt concentrations, peak heights, and shapes is experimentally demonstrated by considering the model's relationship between the normalized gradient slope and the eluting salt concentration. Further scale-up simulations show improved model predictions when radial inhomogeneities in packing quality are considered.

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.

Separation of closely related monoclonal antibody charge variant impurities using poly(ethylenimine)-grafted cation-exchange chromatography resin

The removal of protein charge variants due to complex chemical and enzymatic modifications like glycosylation, fragmentation and deamidation presents a significant challenge in the purification of monoclonal antibodies (mAb) and complicates downstream processing. These protein modifications occur either in vivo or during fermentation and downstream processing. The presence of charge variants can lead to diminished biological activity, differences in pharmacokinetics, pharmacodynamics, stability and efficacy. Therefore, these different product variants should be appropriately controlled for the consistency of product quality and to ensure patient safety. This investigation focuses on the development of a chromatography step for the removal of the charge variants from a recombinant single-chain variable antibody fragment (scFv-Fc-Ab). Poly(ethyleneimine)-grafted cation-exchange resins (Poly CSX and Poly ABX) were evaluated and compared to traditional macroporous cation-exchange and tentacle cation-exchange resins. Linear salt gradient experiments were conducted to study the separation efficiency of scFv-Fc-Ab variants using different resins. A classical thermodynamic model was used to develop a mechanistic understanding of the differences in charge variant retention behaviour of different resins. High selectivity in separation of scFv-Fc-Ab charge variants is obtained in the Poly CSX resin.

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.

Robust mechanistic modeling of protein ion-exchange chromatography

Mechanistic models for ion-exchange chromatography of proteins are well-established and a broad consensus exists on most aspects of the detailed mathematical and physical description. A variety of specializations of these models can typically capture the general locations of elution peaks, but discrepancies are often observed in peak position and shape, especially if the column load level is in the non-linear range. These discrepancies may prevent the use of models for high-fidelity predictive applications such as process characterization and development of high-purity and -productivity process steps. Our objective is to develop a sufficiently robust mechanistic framework to make both conventional and anomalous phenomena more readily predictable using model parameters that can be evaluated based on independent measurements or well-accepted correlations. This work demonstrates the implementation of this approach for industry-relevant case studies using both a model protein, lysozyme, and biopharmaceutical product monoclonal antibodies, using cation-exchange resins with a variety of architectures (SP Sepharose FF, Fractogel EMD SO₃^-, Capto S and Toyopearl SP650M). The modeling employs the general rate model with the extension of the surface diffusivity to be variable, as a function of ionic strength or binding affinity. A colloidal isotherm that accounts for protein-surface and protein-protein interactions independently was used, with each characterized by a parameter determined as a function of ionic strength and pH. Both of these isotherm parameters, along with the variable surface diffusivity, were successfully estimated using breakthrough data at different ionic strengths and pH. The model developed was used to predict overloads and elution curves with high accuracy for a wide variety of gradients and different flow rates and protein loads. The in-silico methodology used in this work for parameter estimation, along with a minimal amount of experimental data, can help the industry adopt model-based optimization and control of preparative ion-exchange chromatography with high accuracy.

Advanced score system and automated search strategies for parameter estimation in mechanistic chromatography modeling

Least squares estimation of unknown parameters from measurement data is a well-established standard method in chromatography modeling but can suffer from critical disadvantages. The description of real-world systems is generally prone to unaccounted mechanisms, such as dispersion in external holdup volumes, and systematic measurement errors, such as caused by pump delays. In this scenario, matching the shape between simulated and measured chromatograms has been found to be more important than the exact peak positions. We have therefore developed a new score system that separately accounts for the shape, position and height of individual peaks. A genetic algorithm is used for optimizing these multiple objectives. Even for non-conflicting objectives, this approach shows superior convergence in comparison to single-objective gradient search, while conflicting objectives indicate incomplete models or inconsistent data. In the latter case, Pareto optima provide important information for understanding the system and improving experiments. The proposed method is demonstrated with synthetic and experimental case studies of increasing complexity. All software is freely available as open source code (https://github.com/modsim/CADET-Match).

Characterization of a macroporous epoxy-polymer based resin for the ion-exchange chromatography of therapeutic proteins

This study investigated the adsorption capacity and mass transfer properties of a novel macroporous epoxy-polymer-based anion-exchanger, MPR Q, for the efficient separation of therapeutic proteins. MPR Q resin was prepared by phase separation based on spinodal decomposition followed by dextran grafting and ligand conjugation. Under static conditions, MPR Q exhibited a binding capacity of 49.8 mg-IgG/cm³-resin at pH 10, whereas the fastest adsorption was observed among the anion-exchanger resins tested. Inverse size-exclusion chromatography (iSEC) experiments revealed that the apparent pore diameter of MPR Q was approximately 90 nm, which was sufficiently large for the penetration of human IgG and bovine IgM. Moreover, the reduced height equivalent to a theoretical plate, h, of human IgG, determined using the linear gradient elution method was 65.8 and was not significantly changed in the range of linear velocities from 20.37 to 50.93 cm/min. The dynamic binding capacity at 10% breakthrough of MPR Q, determined by frontal analysis, exhibited a capacity of 43.8 mg/cm³ at 5.09 cm/min and 58% of DBC_10% was maintained even though the linear velocity was increased to 50.93 cm/min. Furthermore, a resolution for separation of IgG and BSA by MPR Q was 1.06 at 5.09 cm/min, while it was higher than that for the conventional resin at all linear velocities from 5.09 cm/min to 50.93 cm/min. Thus, it was suggested that the MPR Q developed in this study is a promising resin that can efficiently separate large biomacromolecules such as human IgG at higher velocities.

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.

Correlation between protein desorption behavior and its adsorption enthalpy change in polymer grafted anion exchange chromatography

Thermodynamic studies on protein adsorption onto chromatographic surfaces mainly focus on the molecular level interaction between proteins and ligands. Yet, not much attention is given to the study of polymer grafted ligand architecture effect on thermodynamic parameters, nor to the relation between chromatographic parameters and the directly obtained thermodynamic parameters. These relations are needed in order to confer meaning and to ease future data interpretation of thermodynamic studies of protein adsorption. In this study, the adsorption of bovine serum albumin monomer (BSA_m) onto chromatographic surfaces with grafted ligands was studied from a thermodynamic point of view together with chromatographic data. Isothermal titration calorimetry (ITC) results showed that BSA_m adsorption is exothermic (ΔH¯^ads < 0) when adsorbs onto Toyopearl GigaCapQ 650 M, Toyopearl Q600AR, and Q Sepharose XL, but endothermic (ΔH¯^ads > 0) when adsorbs onto Toyopearl SuperQ and a conventional resin (Q Sepharose Fast Flow), showing clear differences in the driving forces of adsorption caused by different ligand architectures. In addition, we found a new relation between the salt required for protein elution and the change in adsorption enthalpy (ΔH¯^ads) directly measured with ITC, intrinsically connecting both adsorption and desorption mechanisms.

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.

Mechanistic Modeling of Preparative Column Chromatography for Biotherapeutics

Chromatography has long been, and remains, the workhorse of downstream processing in the production of biopharmaceuticals. As bioprocessing has matured, there has been a growing trend toward seeking a detailed fundamental understanding of the relevant unit operations, which for some operations include the use of mechanistic modeling in a way similar to its use in the conventional chemical process industries. Mechanistic models of chromatography have been developed for almost a century, but although the essential features are generally understood, the specialization of such models to biopharmaceutical processing includes several areas that require further elucidation. This review outlines the overall approaches used in such modeling and emphasizes current needs, specifically in the context of typical uses of such models; these include selection and improvement of isotherm models and methods to estimate isotherm and transport parameters independently. Further insights are likely to be aided by molecular-level modeling, as well as by the copious amounts of empirical data available for existing processes.