Abstract
Over the past decade, therapeutic monoclonal antibodies (mAbs) have established their role as valuable agents in the treatment of various diseases ranging from cancers to infectious, cardiovascular and autoimmune diseases. Reactive groups of the amino acids within these proteins make them susceptible to many kinds of chemical modifications during manufacturing, storage and in vivo circulation. Among these reactions, the oxidation of methionine residues to their sulfoxide form is a commonly observed chemical modification in mAbs. When the oxidized methionine is in the complementarity-determining region (CDR), this modification can affect antigen binding and thus abrogate biological activity. For these reasons, it is essential to identify oxidation liabilities during the antibody discovery and development phases. Here, we present an in silico method, based on protein modeling and molecular dynamics simulations, to predict the oxidation-liable residues in the variable region of therapeutic antibodies. Previous studies have used the 2-shell water coordination number descriptor (WCN) to identify methionine residues susceptible to oxidation. Although the WCN descriptor successfully predicted oxidation liabilities when the residue was solvent exposed, the method was much less accurate for partially buried methionine residues. Consequently, we introduce a new descriptor, WCN-OH, that improves the accuracy of prediction of methionine oxidation susceptibility by extending the theoretical framework of the water coordination number to incorporate the effects of polar amino acids side chains in close proximity to the methionine of interest.