Abstract:Computationally guided semirational design has significant potential for improving the aggregation kinetics of protein biopharmaceuticals. While improvement in the global conformational stability can stabilize proteins to aggregation under some conditions, previous studies suggest that such an approach is limited, because thermal transition temperatures ( Tm) and the fraction of protein unfolded ( fT) tend to only correlate with aggregation kinetics where the protein is incubated at temperatures approaching the Tm. This is because under these conditions, aggregation from globally unfolded protein becomes dominant. However, under native conditions, the aggregation kinetics are presumed to be dependent on local structural fluctuations or partial unfolding of the native state, which reveal regions of high propensity to form protein-protein interactions that lead to aggregation. In this work, we have targeted the design of stabilizing mutations to regions of the A33 Fab surface structure, which were predicted to be more flexible. This Fab already has high global stability, and global unfolding is not the main cause of aggregation under most conditions. Therefore, the aim was to reduce the conformational flexibility and entropy of the native protein at various locations and thus identify which of those regions has the greatest influence on the aggregation kinetics. Highly dynamic regions of structure were identified through both molecular dynamics simulation and B-factor analysis of related X-ray crystal structures. The most flexible residues were mutated into more stable variants, as predicted by Rosetta, which evaluates the ΔΔ GND for each potential point mutation. Additional destabilizing variants were prepared as controls to evaluate the prediction accuracy and also to assess the general influence of conformational stability on aggregation kinetics. The thermal conformational stability, and aggregation rates of 18 variants at 65 °C, were each examined at pH 4, 200 mM ionic strength, under which conditions the initial wild-type protein was <5% unfolded. Variants with decreased Tm values led to more rapid aggregation due to an increase in the fraction of protein unfolded under the conditions studied. As expected, no significant improvements were observed in the global conformational stability as measured by Tm. However, 6 of the 12 stable variants led to an increase in the cooperativity of unfolding, consistent with lower conformational flexibility and entropy in the native ensemble. Three of these had 5-11% lower aggregation rates, and their structural clustering indicated that the local dynamics of the C-terminus of the heavy chain had a role in influencing the aggregation rate.

AbMelt: Learning antibody thermostability from molecular dynamics

Thermal stability of single-domain antibodies estimated by molecular dynamics simulations

Can template-based protein models guide the design of sequence fitness for enhanced thermal stability of single domain antibodies?

Thermal Stability Estimation of Single Domain Antibodies Using Molecular Dynamics Simulations

TEMPRO: nanobody melting temperature estimation model using protein embeddings

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

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

Formulation development of therapeutic monoclonal antibodies using high-throughput fluorescence and static light scattering techniques: role of conformational and colloidal stability.

DSC Derived (Ea & ΔG) Energetics and Aggregation Predictions for mAbs

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

Prediction of protein biophysical traits from limited data: a case study on nanobody thermostability through NanoMelt

ThermoPCD: a database of molecular dynamics trajectories of antibody–antigen complexes at physiologic and fever-range temperatures

Designing Robust Monoclonal Antibody Drug Products: Pitfalls of Simplistic Approaches for Stability Prediction

Delicate balance among thermal stability, binding affinity, and conformational space explored by single-domain VHH antibodies

Thermal Denaturing of Mutant Lysozyme with Both the OPLSAA and the CHARMM Force Fields

Toward enhancement of antibody thermostability and affinity by computational design in the absence of antigen

Structural mechanism of Fab domain dissociation as a measure of interface stability

Increasing the melting temperature of VHH with the in silico free energy score

Learning the heterogeneous hypermutation landscape of immunoglobulins from high-throughput repertoire data

Investigating Substitutions in Antibody–Antigen Complexes Using Molecular Dynamics: A Case Study with Broad-spectrum, Influenza A Antibodies

An Efficient Method to Predict Protein Thermostability in Alanine Mutation