A multiscale modeling method for therapeutic antibodies in ion exchange chromatography

High throughput screening for rapid and reliable prediction of monovalent antibody binding behavior in flowthrough mode

Flowthrough (FT) anion exchange (AEX) chromatography is a widely used polishing step for the purification of monoclonal antibody (mAb) formats. To accelerate downstream process development, high throughput screening (HTS) tools have proven useful. In this study, the binding behavior of six monovalent mAbs (mvAbs) was investigated by HTS in batch binding mode on different AEX and mixed-mode resins at process-relevant pH and NaCl concentrations. The HTS entailed the evaluation of mvAb partition coefficients (Kp) and visualization of results in surface-response models. Interestingly, the HTS data grouped the mvAbs into either a strong-binding group or a weak-binding/FT group independent of theoretical Isoelectric point. Mapping the charged and hydrophobic patches by in silico protein surface property analyses revealed that the distribution of patches play a major role in predicting FT behavior. Importantly, the conditions identified by HTS were successfully verified by 1 mL on-column experiments. Finally, employing the optimal FT conditions (7-9 mS/cm and pH 7.0) at a mini-pilot scale (CV = 259 mL) resulted in 99% yield and a 21-23-fold reduction of host cell protein to <100 ppm, depending on the varying host cell protein (HCP) levels in the load. This work opens the possibility of using HTS in FT mode to accelerate downstream process development for mvAb candidates in early research.

Current state of implementation of in silico tools in the biopharmaceutical industry-Proceedings of the 5th modeling workshop

The fifth modeling workshop (5MW) was held in June 2023 at Favrholm, Denmark and sponsored by Recovery of Biological Products Conference Series. The goal of the workshop was to assemble modeling practitioners to review and discuss the current state, progress since the last fourth mini modeling workshop (4MMW), gaps and opportunities for development, deployment and maintenance of models in bioprocess applications. Areas of focus were four categories: biophysics and molecular modeling, mechanistic modeling, computational fluid dynamics (CFD) and plant modeling. Highlights of the workshop included significant advancements in biophysical/molecular modeling to novel protein constructs, mechanistic models for filtration and initial forays into modeling of multiphase systems using CFD for a bioreactor and mapped strategically to cell line selection/facility fit. A significant impediment to more fully quantitative and calibrated models for biophysics is the lack of large, anonymized datasets. A potential solution would be the use of specific descriptors in a database that would allow for detailed analyzes without sharing proprietary information. Another gap identified was the lack of a consistent framework for use of models that are included or support a regulatory filing beyond the high-level guidance in ICH Q8-Q11. One perspective is that modeling can be viewed as a component or precursor of machine learning (ML) and artificial intelligence (AI). Another outcome was alignment on a key definition for "mechanistic modeling." Feedback from participants was that there was progression in all of the fields of modeling within scope of the conference. Some areas (e.g., biophysics and molecular modeling) have opportunities for significant research investment to realize full impact. However, the need for ongoing research and development for all model types does not preclude the application to support process development, manufacturing and use in regulatory filings. Analogous to ML and AI, given the current state of the four modeling types, a prospective investment in educating inter-disciplinary subject matter experts (e.g., data science, chromatography) is essential to advancing the modeling community.

Predicting multimodal chromatography of therapeutic antibodies using multiscale modeling

Multimodal chromatography has emerged as a powerful method for the purification of therapeutic antibodies. However, process development of this separation technique remains challenging because of an intricate and molecule-specific interaction towards multimodal ligands, leading to time-consuming and costly experimental optimization. This study presents a multiscale modeling approach to predict the multimodal chromatographic behavior of therapeutic antibodies based on their sequence information. Linear gradient elution (LGE) experiments were performed on an anionic multimodal resin for 59 full-length antibodies, including five different antibody formats at pH 5.0, 6.0, and 7.0 that were used for parameter determination of a linear adsorption model at low loading density conditions. Quantitative structure-property relationship (QSPR) modeling was utilized to correlate the adsorption parameters with up to 1374 global and local physicochemical descriptors calculated from antibody homology models. The final QSPR models employed less than eight descriptors per model and demonstrated high training accuracy (R² > 0.93) and reasonable test set prediction accuracy (Q² > 0.83) for the adsorption parameters. Model evaluation revealed the significance of electrostatic interaction and hydrophobicity in determining the chromatographic behavior of antibodies, as well as the importance of the HFR3 region in antibody binding to the multimodal resin. Chromatographic simulations using the predicted adsorption parameters showed good agreement with the experimental data for the vast majority of antibodies not employed during the model training. The results of this study demonstrate the potential of sequence-based prediction for determining chromatographic behavior in therapeutic antibody purification. This approach leads to more efficient and cost-effective process development, providing a valuable tool for the biopharmaceutical industry.

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.

Antibody sequence-based prediction of pH gradient elution in multimodal chromatography

Multimodal chromatography has emerged as a promising technique for antibody purification, owing to its capacity to selectively capture and separate target molecules. However, the optimization of chromatography parameters remains a challenge due to the intricate nature of protein-ligand interactions. To tackle this issue, efficient predictive tools are essential for the development and optimization of multimodal chromatography processes. In this study, we introduce a methodology that predicts the elution behavior of antibodies in multimodal chromatography based on their amino acid sequences. We analyzed a total of 64 full-length antibodies, including IgG1, IgG4, and IgG-like multispecific formats, which were eluted using linear pH gradients from pH 9.0 to 4.0 on the anionic mixed-mode resin Capto adhere. Homology models were constructed, and 1312 antibody-specific physicochemical descriptors were calculated for each molecule. Our analysis identified six key structural features of the multimodal antibody interaction, which were correlated with the elution behavior, emphasizing the antibody variable region. The results show that our methodology can predict pH gradient elution for a diverse range of antibodies and antibody formats, with a test set R² of 0.898. The developed model can inform process development by predicting initial conditions for multimodal elution, thereby reducing trial and error during process optimization. Furthermore, the model holds the potential to enable an in silico manufacturability assessment by screening target antibodies that adhere to standardized purification conditions. In conclusion, this study highlights the feasibility of using structure-based prediction to enhance antibody purification in the biopharmaceutical industry. This approach can lead to more efficient and cost-effective process development while increasing process understanding.