Abstract:High drug dosages required for biotherapeutics, such as monoclonal antibodies (mAbs), and the small volumes that can be administered to patients via subcutaneous injections pose challenges due to high-concentration formulations. The addition of excipients such as arginine to high-concentration protein formulations can increase solubility and reduce the tendency of protein particle formation. Studying highconcentration mAb systems with molecular dynamics (MD) simulations can provide microscopic insights into the mode of action of excipients but requires large system sizes and long time scales that are currently out of reach at the fully atomistic level. Computationally efficient coarse grained models such as the Martini 3 force field can tackle this challenge but require careful parametrization, testing, and validation. This study extends the popular Martini 3 force field towards realistic protein-excipient interactions of arginine and glutamate excipients, using the Fab domains of the therapeutic mAbs trastuzumab and omalizumab as model systems. A novel Martini 3 mapping of the amino acid excipients is introduced, which explicitly captures the zwitterionic character of the backbone. The Fab-excipient interactions of arginine and glutamate are characterized at the single-residue level concerning molecular contacts. The Martini 3 simulations are compared with results from all-atom simulations as a reference. Our findings reveal an overestimation of Fab-excipient contacts in the original Martini 3 force field, suggesting a too strong attraction between protein surface residues and excipients. Therefore, we reparametrized the protein-excipient interaction parameters in Martini 3 against all-atom simulations. The excipient interactions obtained with the new Martini 3 mapping and Lennard-Jones (LJ) interaction parameters, coined Martini 3-exc, agree closely with the all-atom reference data. This work presents an improved parameter set for mAb-arginine and mAb-glutamate interactions in the Martini 3 coarse grained force field, enabling large-scale simulations of high-concentration mAb formulations and the stabilizing effect of the excipients.

Computational Procedure for Predicting Excipient Effects on Protein-Protein Affinities

Modeling the protein partitioning in aqueous polymer two-phase systems: Influence of polymer concentration and molecular weight

A Discontinuous Potential Model for Protein-Protein Interactions.

Computational prediction of protein–protein binding affinities

Residue-resolved insights into the stabilization of therapeutic proteins by excipients: A case study of two monoclonal antibodies with arginine and glutamate

Optimized Protein-Excipient Interactions in the Martini 3 Force Field

Complementary evaluation of computational methods for predicting single residue effects on peptide binding specificities

Impact of Co-Excipient Selection on Hydrophobic Polymer Folding: Insights for Optimal Formulation Design

Predicting Colloidal Stability of High-Concentration Monoclonal Antibody Formulations in Common Pharmaceutical Buffers Using Improved Polyethylene Glycol Induced Protein Precipitation Assay.

Application of a thermodynamic model to the prediction of phase separations in freeze-concentrated formulations for protein lyophilization

Towards safer and efficient formulations: Machine learning approaches to predict drug-excipient compatibility

Residue interactions guide translational diffusion of proteins

Computational Prediction of Amino Acid Preferences of Potentially Multispecific Peptide-Binding Domains Involved in Protein-Protein Interactions

Prediction of protein–ligand binding affinity by free energy simulations: assumptions, pitfalls and expectations

Predictive Model Building for Aggregation Kinetics Based on Molecular Dynamics Simulations of an Antibody Fragment

Learning effective amino acid interactions through iterative stochastic techniques

Computational Structure Prediction for Antibody-Antigen Complexes From Hydrogen-Deuterium Exchange Mass Spectrometry: Challenges and Outlook

Explainable Machine Learning Models to Predict Gibbs–Donnan Effect During Ultrafiltration and Diafiltration of High‐Concentration Monoclonal Antibody Formulations

Molecular Dynamics Study of Protein Aggregation at Moving Interfaces

On achieving high accuracy and reliability in the calculation of relative protein–ligand binding affinities

Massively parallel free energy calculations for affinity maturation of designed miniproteins