Comparison of the pH- and thermally-induced fluctuations of a therapeutic antibody Fab fragment by molecular dynamics simulation

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

Computational methods including machine learning and molecular dynamics simulations have strong potential to characterize, understand, and ultimately predict the properties of proteins relevant to their stability and function as therapeutics. Such methods would streamline the development pathway by minimizing the current experimental testing required for many protein variants and formulations. The molecular understanding of thermostability and aggregation propensity has advanced significantly along with predictive algorithms based on the sequence-level or structural-level information on a protein. However, these approaches focus largely on a comparison of protein sequence variations to correlate the properties of proteins to their stability, solubility, and aggregation propensity. For therapeutic protein development, it is of equal importance to take into account the impact of the formulation conditions to elucidate and predict the stability of the antibody drugs. At the macroscopic level, changing temperature, pH, ionic strength, and the addition of excipients can significantly alter the kinetics of protein aggregation. The mechanisms controlling aggregation kinetics have been traced back to a combination of molecular features, including conformational stability, partial unfolding to aggregation-prone states, and the colloidal stability governed by surface charges and hydrophobicity. However, very little has been done to evaluate these features in the context of protein dynamics in different formulations. In this work, we have combined a range of molecular features calculated from the Fab A33 protein sequence and molecular dynamics simulations. Using the power of advanced, yet interpretable, statistical tools, it has been possible to uncover greater insights into the mechanisms behind protein stability, validating previous findings, and also develop models that can predict the aggregation kinetics within a range of 49 different solution conditions.

Stabilization of adalimumab Fab through the introduction of disulfide bonds between the variable and constant domains

Fab is a promising format for antibody drug. Therefore, efforts have been made to improve its thermal stability for therapeutic and commercial use. So far, we have attempted to introduce a disulfide bond into the Fab fragment to improve its thermal stability and demonstrated that it is possible to do this without sacrificing its biochemical function. In this study, to develop a novel stabilization strategy for Fab, we attempted to introduce a disulfide bond between the variable and constant domains and prepared three variants of Fab; H:G10C + H:P210C, L:P40C + L:E165C, and H:G10C + H:P210C + L:P40C + L:E165C. Differential scanning calorimetry measurements showed that each of these variants had improved thermal stability. In addition, the variants with two disulfide bonds demonstrated a 6.5 °C increase in their denaturation temperatures compared to wild-type Fab. The introduction of disulfide bonds was confirmed by X-ray crystallography, and the variants retained their antigen-binding activity. The variants were also found to be less aggregative than the wild type. Our results demonstrate that the introduction of a disulfide bond between the variable and constant domains significantly improves the thermal stability of Fab.

Analysis of thermostability for seven Phe to Ala and six Pro to Gly mutants in the Fab constant region of adalimumab

To identify amino acids that play important roles in the structural stability of Fab, seven phenylalanine residues in the Fab constant region of the therapeutic antibody adalimumab were subjected to alanine mutagenesis. Six Fab mutants, H:F130A, H:F154A, H:F174A, L:F118A, L:F139A and L:F209A, showed decreased thermostability compared with wild-type Fab. In contrast, the Tm for the L:F116A mutant was 1.7°C higher than that of wild-type Fab, indicating that the F116 residue was unfavorable for Fab thermostability. Six proline mutants, H:P131G, H:P155G, H:P175G, L:P119G, L:P120G and L:P141G, were also prepared to investigate the effect of proline residues adjacent to mutated phenylalanine residues. The thermostability of the H:P155G and L:P141G mutants in particular was significantly reduced, with decreases in Tm of 5.0 and 3.0°C, respectively, compared with wild-type Fab. The H:P155 and L:P141 residues have a cis conformation, whereas the other mutated proline residues have a trans conformation. H:P155 and L:P141 had stacking interactions with the H:F154 and L:Y140, respectively, at the interface between the variable and constant regions. It is suggested that the interactions of the aromatic ring with a cis-form proline at the interface between the variable and constant regions is important for stability of Fab.

Crystal structures and molecular dynamics simulations of a humanised antibody fragment at acidic to basic pH

Antibody-fragment (Fab) therapy development has the potential to be accelerated by computational modelling and simulations that predict their target binding, stability, formulation, manufacturability, and the impact of further protein engineering. Such approaches are currently predicated on starting with good crystal structures that closely represent those found under the solution conditions to be simulated. A33 Fab, is an undeveloped immunotherapeutic antibody candidate that was targeted to the human A33 antigen homogeneously expressed in 95% cases of primary and metastatic colorectal cancers. It is now used as a very well characterised testing ground for developing analytics, formulation and protein engineering strategies, and to gain a deeper understanding of mechanisms of destabilisation, representative of the wider therapeutic Fab platform. In this article, we report the structure of A33 Fab in two different crystal forms obtained at acidic and basic pH. The structures overlapped with RMSD of 1.33 Å overall, yet only 0.5 Å and 0.76 Å for the variable- and constant regions alone. While most of the differences were within experimental error, the switch linker between the variable and the constant regions showed some small differences between the two pHs. The two structures then enabled a direct evaluation of the impact of initial crystal structure selection on the outcomes of molecular dynamics simulations under different conditions, and their subsequent use for determining best fit solution structures using previously obtained small-angle x-ray scattering (SAXS) data. The differences in the two structures did not have a major impact on MD simulations regardless of the pH, other than a slight persistence of structure affecting the solvent accessibility of one of the predicted APR regions of A33 Fab. Interestingly, despite being obtained at pH 4 and pH 9, the two crystal structures were more similar to the SAXS solution structures obtained at pH 7, than to those at pH 4 or pH 9. Furthermore, the P65 crystal structure from pH 4 was also a better representation of the solution structures at any other pH, than was the P1 structure obtained at pH 9. Thus, while obtained at different pH, the two crystal structures may represent highly (P65) and lesser (P1) populated species that both exist at pH 7 in solution. These results now lay the foundation for confident MD simulations of A33 Fab that rationalise or predict behaviours in a range of conditions.

Getting Smaller by Denaturation: Acid-Induced Compaction of Antibodies

Protein denaturation is a ubiquitous process that occurs both in vitro and in vivo. While our molecular understanding of the denatured structures of proteins is limited, it is commonly accepted that the loss of unique intramolecular contacts makes proteins larger. Herein, we report compaction of the immunoglobulin G1 (IgG1) protein upon acid denaturation. Small-angle X-ray scattering coupled with size exclusion chromatography revealed that IgG1 radii of gyration at pH 2 were ∼75% of those at a neutral pH. Scattering profiles showed a compact globular shape, supported by analytical ultracentrifugation. The acid denaturation of proteins with a decrease in size is energetically costly, and acid-induced compaction requires an attractive force for domain reorientation. Such intramolecular aggregation may be widespread in immunoglobulin proteins as noncanonical structures. Herein, we discuss the potential biological significance of these noncanonical structures of antibodies.

Protein Engineering and HDX Identify Structural Regions of G-CSF Critical to Its Stability and Aggregation

The protein engineering and formulation of therapeutic proteins for prolonged shelf-life remain a major challenge in the biopharmaceutical industry. Understanding the influence of mutations and formulations on the protein structure and dynamics could lead to more predictive approaches to their improvement. Previous intrinsic fluorescence analysis of the chemically denatured granulocyte colony-stimulating factor (G-CSF) suggested that loop AB could subtly reorganize to form an aggregation-prone intermediate state. Hydrogen deuterium exchange mass spectrometry (HDX-MS) has also revealed that excipient binding increased the thermal unfolding transition midpoint (T_m) by stabilizing loop AB. Here, we have combined protein engineering with biophysical analyses and HDX-MS to reveal that increased exchange in a core region of the G-CSF comprising loop AB (ABI, a small helix, ABII) and loop CD packed onto helix B and the beginning of loop BC leads to a decrease in T_m and higher aggregation rates. Furthermore, some mutations can increase the population of the aggregation-prone conformation within the native ensemble, as measured by the greater local exchange within this core region.