Abstract:The design of stable formulations remains a major challenge for protein therapeutics, particularly the need to minimize aggregation. Experimental formulation screens are typically based on thermal transition midpoints (Tm), and forced degradation studies at elevated temperatures. Both approaches give limited predictions of long-term storage stability, particularly at low temperatures. Better understanding of the mechanisms of action for formulation of excipients and buffers could lead to improved strategies for formulation design. Here, we identified a complex impact of glycine concentration on the experimentally determined stability of an antibody Fab fragment and then used molecular dynamics simulations to reveal mechanisms that underpin these complex behaviors. Tm values increased monotonically with glycine concentration, but associated ΔSvh measurements revealed more complex changes in the native ensemble dynamics, which reached a maximum at 30 mg/mL. The aggregation kinetics at 65 °C were similar at 0 and 20 mg/mL glycine, but then significantly slower at 50 mg/mL. These complex behaviors indicated changes in the dominant stabilizing mechanisms as the glycine concentration was increased. MD revealed a complex balance of glycine self-interaction, and differentially preferred interactions of glycine with the Fab as it displaced hydration-shell water, and surface-bound water and citrate buffer molecules. As a result, glycine binding to the Fab surface had different effects at different concentrations, and led from preferential interactions at low concentrations to preferential exclusion at higher concentrations. During preferential interaction, glycine displaced water from the Fab hydration shell, and a small number of water and citrate molecules from the Fab surface, which reduced the protein dynamics as measured by root-mean-square fluctuation (RMSF) on the short time scales of MD. By contrast, the native ensemble dynamics increased according to ΔSvh, suggesting increased conformational changes on longer time scales. The aggregation kinetics did not change at low glycine concentrations, and so the opposing dynamics effects either canceled out or were not directly relevant to aggregation. During preferential exclusion at higher glycine concentrations, glycine could only bind to the Fab surface through the displacement of citrate buffer molecules already favorably bound on the Fab surface. Displacement of citrate increased the flexibility (RMSF) of the Fab, as glycine formed fewer bridging hydrogen bonds to the Fab surface. Overall, the slowing of aggregation kinetics coincided with reduced flexibility in the Fab ensemble at the very highest glycine concentrations, as determined by both RMSF and ΔSvh, and occurred at a point where glycine binding displaced neither water nor citrate. These final interactions with the Fab surface were driven by mass action and were the least favorable, leading to a macromolecular crowding effect under the regime of preferential exclusion that stabilized the dynamics of Fab.

Molecular Dynamics Simulation of an Iron(III) Binding Site on the Fc Domain of IgG1 Relevant for Visible Light-Induced Protein Fragmentation.

Molecular Recognition of Fc-Specific Ligands Binding Onto the Consensus Binding Site of Igg: Insights from Molecular Simulation

Molecular Simulation of the Interactions between 4-Mercaptoethyl-pyridine Ligand and IgG

Molecular Insight into Protein Binding Orientations and Interaction Modes on Hydrophobic Charge-Induction Resin.

Visible Light Induces Site-Specific Oxidative Heavy Chain Fragmentation of a Monoclonal Antibody (IgG1) Mediated by an Iron(III)-Containing Histidine Buffer

Fragment Density Functional Theory Calculation of NMR Chemical Shifts for Proteins with Implicit Solvation.

Photoinduced Tyrosine Side Chain Fragmentation in IgG4-Fc: Mechanisms and Solvent Isotope Effects

Effects of Ligand Binding on the Mechanical Stability of Protein GB1 Studied by Steered Molecular Dynamics Simulation

The Solvent at Antigen-Binding Site Regulated C3d-CR2 Interactions Through the C-terminal Tail of C3d at Different Ion Strengths: Insights from Molecular Dynamics Simulation.

Comparison of the Ph- and Thermally-Induced Fluctuations of a Therapeutic Antibody Fab Fragment by Molecular Dynamics Simulation

Computational Design to Reduce Conformational Flexibility and Aggregation Rates of an Antibody Fab Fragment.

Insight into the photodynamic mechanism and protein binding of a nitrosyl iron-sulfur [Fe2S2(NO)4]2- cluster

Stabilization of the Ferryl=oxoheme Form of Staphylococcus aureus IsdG by Electron Transfer from a Second-Sphere Tryptophan

Electrostatic Affinities and Binding Kinetics among $α_2\textrm{I}$ Integrin Domains, Divalent Cations and 21-mer Collagen Fragment

Crystal Structure, Theoretical Analysis, and Protein/DNA Binding Activity of Iron(III) Complex Containing Differently Protonated Pyridoxal- S-Methyl-Isothiosemicarbazone Ligands

Molecular Dynamics Simulation and Kinetic Study of Fluoride Binding to V21C/V66C Myoglobin with a Cytoglobin-like Disulfide Bond

Iron insertion into coproporphyrin III‐ferrochelatase complex: evidence for an intermediate distorted catalytic species

Molecular Dynamics Simulations Reveal How Competing Protein-Surface Interactions for Glycine, Citrate, and Water Modulate Stability in Antibody Fragment Formulations

Ligand Binding Rate Constants in Heme Proteins Using Markov State Models and Molecular Dynamics Simulations

Sequential Proton Loss Electron Transfer in Deactivation of Iron(IV) Binding Protein by Tyrosine Based Food Components

Structure and dynamics of a dizinc metalloprotein: Effect of charge transfer and polarization