Viscosity Prediction of High-Concentration Antibody Solutions with Atomistic Simulations.

The computational prediction of the viscosity of dense protein solutions is highly desirable, for example, in the early development phase of high-concentration biopharmaceutical formulations where the material needed for experimental determination is typically limited. Here, we use large-scale atomistic molecular dynamics (MD) simulations with explicit solvation to de novo predict the dynamic viscosities of solutions of a monoclonal IgG1 antibody (mAb) from the pressure fluctuations using a Green-Kubo approach. The viscosities at simulated mAb concentrations of 200 and 250 mg/mL are compared to the experimental values, which we measured with rotational rheometry. The computational viscosity of 24 mPa·s at the mAb concentration of 250 mg/mL matches the experimental value of 23 mPa·s obtained at a concentration of 213 mg/mL, indicating slightly different effective concentrations (or activities) in the MD simulations and in the experiments. This difference is assigned to a slight underestimation of the effective mAb-mAb interactions in the simulations, leading to a too loose dynamic mAb network that governs the viscosity. Taken together, this study demonstrates the feasibility of all-atom MD simulations for predicting the properties of dense mAb solutions and provides detailed microscopic insights into the underlying molecular interactions. At the same time, it also shows that there is room for further improvements and highlights challenges, such as the massive sampling required for computing collective properties of dense biomolecular solutions in the high-viscosity regime with reasonable statistical precision.

Exploring molecular determinants and pharmacokinetic properties of IgG1-scFv bispecific antibodies.

Bispecific antibodies (BsAbs) capable of recognizing two distinct epitopes or antigens offer promising therapeutic options for various diseases by targeting multiple pathways. The favorable pharmacokinetic (PK) properties of monoclonal antibodies (mAbs) are crucial, as they directly influence patient safety and therapeutic efficacy. For numerous mAb therapeutics, optimization of neonatal Fc receptor (FcRn) interactions and elimination of unfavorable molecular properties have led to improved PK properties. However, many BsAbs exhibit unfavorable PK, which has precluded their development as drugs. In this report, we present studies on the molecular determinants underlying the distinct PK profiles of three IgG1-scFv BsAbs. Our study indicated that high levels of nonspecific interactions, elevated isoelectric point (pI), and increased number of positively charged patches contributed to the fast clearance of IgG1-scFv. FcRn chromatography results revealed specific scFv-FcRn interactions that are unique to the IgG1-scFv, which was further supported by molecular dynamics (MD) simulation. These interactions likely stabilize the BsAb FcRn interaction at physiological pH, which in turn could disrupt FcRn-mediated BsAb recycling. In addition to the empirical observations, we also evaluated the impact of in silico properties, including pI differential between the Fab and scFv and the ratio of dipole moment to hydrophobic moment (RM) and their correlation with the observed clearance. These findings highlight that the PK properties of BsAbs may be governed by novel determinants, owing to their increased structural complexity compared to immunoglobulin G (IgG) 1 antibodies.

Entropy Tug-of-War Determines Solvent Effects in the Liquid-Liquid Phase Separation of a Globular Protein.

Liquid-liquid phase separation (LLPS) plays a key role in the compartmentalization of cells via the formation of biomolecular condensates. Here, we combined atomistic molecular dynamics (MD) simulations and terahertz (THz) spectroscopy to determine the solvent entropy contribution to the formation of condensates of the human eye lens protein γD-Crystallin. The MD simulations reveal an entropy tug-of-war between water molecules that are released from the protein droplets and those that are retained within the condensates, two categories of water molecules that were also assigned spectroscopically. A recently developed THz-calorimetry method enables quantitative comparison of the experimental and computational entropy changes of the released water molecules. The strong correlation mutually validates the two approaches and opens the way to a detailed atomic-level understanding of the different driving forces underlying the LLPS.