Abstract
Asparagine (Asn) deamidation and aspartic acid (Asp) isomerization are common degradation pathways that affect the stability of therapeutic antibodies. These modifications can pose a significant challenge in the development of biopharmaceuticals. As such, the early engineering and selection of chemically stable monoclonal antibodies (mAbs) can substantially mitigate the risk of subsequent failure. In this study, we introduce a novel in silico approach for predicting deamidation and isomerization sites in therapeutic antibodies by analyzing the structural environment surrounding asparagine and aspartate residues. The resulting quantitative structure-activity relationship (QSAR) model was trained using previously published forced degradation data from 57 clinical-stage mAbs. The predictive accuracy of the model was evaluated for four different states of the protein structure: (1) static homology models, (2) enhancing low-frequency vibrational modes during short molecular dynamics (MD) runs, (3) a combination of (2) with a protonation state reassignment, and (4) conventional full-atomistic MD simulations. The most effective QSAR model considered the accessible surface area (ASA) of the residue, the pKa value of the backbone amide, and the root mean square deviations of both the alpha carbon and the side chain. The accuracy was further enhanced by incorporating the QSAR model into a decision tree, which also includes empirical information about the sequential successor and the position in the protein. The resulting model has been implemented as a plugin named "Forecasting Reactivity of Isomerization and Deamidation in Antibodies" in MOE software, completed with a user-friendly graphical interface to facilitate its use.