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.

Computational Mutagenesis of Antibody Fragments: Disentangling Side Chains from ΔΔ G Predictions

Locuaz: an in-silico platform for antibody fragments optimization

Improving the Stability of an Antibody Variable Fragment by a Combination of Knowledge-based Approaches: Validation and Mechanisms

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

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

Identification of Strategic Residues at the Interface of Antigen–Antibody Interactions by in Silico Mutagenesis

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

Antibodies: Computer-Aided Prediction of Structure and Design of Function

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

Affinity maturation of antibody fragments: A review encompassing the development from random approaches to computational rational optimization

Binding symmetry and surface flexibility mediate antibody self-association

Assessing and Engineering Antibody Stability Using Experimental and Computational Methods.

Predicting deamidation and isomerization sites in therapeutic antibodies using structure-based in silico approaches

Experimental Structures of Antibody/MHC-I Complexes Reveal Details of Epitopes Overlooked by Computational Prediction

A Comparison of Antibody-Antigen Complex Sequence-to-Structure Prediction Methods and their Systematic Biases

Antibody complementarity-determining region design using AlphaFold2 and DDG predictor

In silico design of stable single-domain antibodies with high affinity

Facile Affinity Maturation of Antibody Variable Domains Using Natural Diversity Mutagenesis

In Silico Methods In Antibody Design

Computational design of developable therapeutic antibodies: efficient traversal of binder landscapes and rescue of escape mutations

Antigen–Antibody Interactions and Structural Flexibility of a Femtomolar-Affinity Antibody