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.

Thermal Stability Estimation of Single Domain Antibodies Using Molecular Dynamics Simulations

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

AbMelt: Learning antibody thermostability from molecular dynamics

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

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

Molecular Origin of the Stability Difference in Four Shark IgNAR Constant Domains.

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

Structural and mutational analysis of a monomeric and dimeric form of a single domain antibody with implications for protein misfolding

Thermodynamic versus kinetic basis for the high conformational stability of nanobodies for therapeutic applications

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

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

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

Evaluation and applications of antibody variable stability

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

Probing Thermal Stability of Proteins with Temperature Scanning Viscometer.

TEMPRO: nanobody melting temperature estimation model using protein embeddings

Assessing and Engineering Antibody Stability Using Experimental and Computational Methods.

High‐throughput system for the thermostability analysis of proteins

Evaluating molecular modeling tools for thermal stability using an independently generated dataset

Molecular Dynamics Perspective on the Protein Thermal Stability: A Case Study Using SAICAR Synthetase

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