We report the use of molecular modeling to predict the oxidation propensity of methionine residues in proteins. Oxidation of methionine to the sulfoxide form is one of the major degradation pathways for therapeutic proteins. Oxidation can occur during production, formulation, or storage of pharmaceuticals and it often reduces or eliminates biological activity. We use a molecular model based on atomistic simulations called 2-shell water coordination number to predict the oxidation rates for several model proteins and therapeutic candidates. In addition, we implement models that are based on static and simulation average of the solvent-accessible area (SAA) for either the side chain or the sulfur atom in the methionine residue. We then compare the results from the different models against the experimentally measured relative rates of methionine oxidation. We find that both the 2-shell model and the simulation-averaged SAA models are accurate in predicting the oxidation propensity of methionine residues for the proteins tested. We also find the appropriate parameter ranges where the models are most accurate. These models have significant predictive power and can be used to enable further protein engineering or to guide formulation approaches in stabilizing the unstable methionine residues.
Keywords: biophysical models; computer aided drug design; in silico modeling; molecular dynamics; molecular modeling; oxidation; physicochemical properties; protein formulation; protein structure; stabilization.
© 2015 Wiley Periodicals, Inc. and the American Pharmacists Association.