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

All-atom modeling of methacrylate-based multi-modal chromatography resins for Langmuir constant prediction of peptides

In downstream processing, the intricate nature of the interactions between biomolecules and adsorbent materials presents a significant challenge in the prediction of their binding and elution behaviors. This complexity is further heightened in multi-modal chromatography (MMC), which employs two distinct binding mechanisms. To gain a deeper understanding of the involved interactions, simulating the adsorption of biomolecules on resin surfaces is a focal point of ongoing research. However, previous studies often simplified the adsorbent surface, modeling it as a flat or slightly curved plane without including a realistic backbone structure. Here, we introduce and validate two novel workflows aimed at predicting peptide binding behaviors in MMC, specifically targeting methacrylate-based resins. Our first achievement was the development of an all-atom model of a commercial MMC resin surface, incorporating its polymethacrylic backbone. Furthermore, we established and tested a workflow for rapid calculations of binding free energies (ΔG) with 10 linear peptides as target molecules. These ΔG calculations were effectively used to predict Langmuir constants, achieving a high coefficient of determination (R²) of 0.96. In subsequent benchmarking tests, our model outperformed established, simpler resin surface models in terms of predictive capabilities.

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.