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.

How hydrophobicity and the glycosylation site of glycans affect protein folding and stability: a molecular dynamics simulation.

Design of Glyco-Linkers at Multiple Structural Levels to Modulate Protein Stability.

Sialylation-induced stabilization of dynamic glycoprotein conformations unveiled by time-aligned parallel unfolding and glycan release mass spectrometry

Solvation and Desolvation Effects in Protein Folding: Native Flexibility, Kinetic Cooperativity and Enthalpic Barriers under Isostability Conditions

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

Conformational Changes of the Protein Domains Upon Binding with Carbon Nanotubes Studied by Molecular Dynamics Simulations

Glycosylation increases active site rigidity leading to improved enzyme stability and turnover

Chevron Behavior and Isostable Enthalpic Barriers in Protein Folding: Successes and Limitations of Simple Gō-like Modeling.

Interplay Between Drying and Stability of a TIM Barrel Protein: a Combined Simulation-Experimental Study.

Molecular Dynamics Simulation of Protein Conformation Transition in Solution

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

Structural Transitions of Confined Model Proteins: Molecular Dynamics Simulation and Experimental Validation

Configurational Entropy Modulates The Mechanical Stability Of Protein Gb1

Effects of Turn Stability and Side-Chain Hydrophobicity on the Folding of Β-Structures

Interactions of sucrose and trehalose with lysozyme in different media: a perspective from atomistic and coarse-grained molecular dynamics simulations

N-Glycosylation as a Modulator of Protein Conformation and Assembly in Disease

Dissecting the Conformational Stability of a Glycan Hairpin

Computational impact of hydrophobicity in protein stability

Entropy-Enthalpy Compensations Fold Proteins in Precise Ways

Oscillatory molecular driving force for protein folding at high concentration: a molecular simulation.

Interactions between amino acid side chains in cylindrical hydrophobic nanopores with applications to peptide stability