In silico selection of therapeutic antibodies for development: Viscosity, clearance, and chemical stability

Significance mAbs are increasingly being used for treatment of chronic diseases wherein the subcutaneous delivery route is preferred to enable self-administration and at-home use. To deliver high doses (several hundred milligrams) through a small volume (∼1 mL) into the subcutaneous space, mAb solutions need to have low viscosity. Concomitantly, acceptable chemical stability is required for adequate shelf life, and normal in vivo clearance is needed for less frequent dosing. We propose in silico tools that provide rapid assessment of atypical behavior of mAbs (high viscosity, chemical degradation, and fast plasma clearance), which are simply predicted from sequence and/or structure-derived parameters. Such analysis will greatly improve the probability of success to move mAb-based therapeutics efficiently into clinical development and ultimately benefit patients. For mAbs to be viable therapeutics, they must be formulated to have low viscosity, be chemically stable, and have normal in vivo clearance rates. We explored these properties by observing correlations of up to 60 different antibodies of the IgG1 isotype. Unexpectedly, we observe significant correlations with simple physical properties obtainable from antibody sequences and by molecular dynamics simulations of individual antibody molecules. mAbs viscosities increase strongly with hydrophobicity and charge dipole distribution and decrease with net charge. Fast clearance correlates with high hydrophobicities of certain complementarity determining regions and with high positive or high negative net charge. Chemical degradation from tryptophan oxidation correlates with the average solvent exposure time of tryptophan residues. Aspartic acid isomerization rates can be predicted from solvent exposure and flexibility as determined by molecular dynamics simulations. These studies should aid in more rapid screening and selection of mAb candidates during early discovery.

[1]  Yaning Wang,et al.  Accurate determination of succinimide degradation products using high fidelity trypsin digestion peptide map analysis. , 2011, Analytical chemistry.

[2]  A. Klibanov,et al.  Hydrophobic salts markedly diminish viscosity of concentrated protein solutions , 2011, Biotechnology and bioengineering.

[3]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[4]  Steven J Shire,et al.  Challenges in the development of high protein concentration formulations. , 2004, Journal of pharmaceutical sciences.

[5]  M. Eisenstein Something new under the skin , 2011, Nature Biotechnology.

[6]  Vincent A. Voelz,et al.  Probing antibody internal dynamics with fluorescence anisotropy and molecular dynamics simulations , 2013, mAbs.

[7]  T. Blundell,et al.  Comparative protein modelling by satisfaction of spatial restraints. , 1993, Journal of molecular biology.

[8]  K N Houk,et al.  Neighboring side chain effects on asparaginyl and aspartyl degradation: an ab initio study of the relationship between peptide conformation and backbone NH acidity. , 2001, Journal of the American Chemical Society.

[9]  W Wang,et al.  Monoclonal Antibody Pharmacokinetics and Pharmacodynamics , 2008, Clinical pharmacology and therapeutics.

[10]  Y John Wang,et al.  Methionine, tryptophan, and histidine oxidation in a model protein, PTH: mechanisms and stabilization. , 2009, Journal of pharmaceutical sciences.

[11]  Collaborative Computational,et al.  The CCP4 suite: programs for protein crystallography. , 1994, Acta crystallographica. Section D, Biological crystallography.

[12]  T. Igawa,et al.  Reduced elimination of IgG antibodies by engineering the variable region. , 2010, Protein engineering, design & selection : PEDS.

[13]  Apollon Papadimitriou,et al.  Structure-Based Prediction of Asparagine and Aspartate Degradation Sites in Antibody Variable Regions , 2014, PloS one.

[14]  Chris Sander,et al.  The double cubic lattice method: Efficient approaches to numerical integration of surface area and volume and to dot surface contouring of molecular assemblies , 1995, J. Comput. Chem..

[15]  Li Li,et al.  Concentration Dependent Viscosity of Monoclonal Antibody Solutions: Explaining Experimental Behavior in Terms of Molecular Properties , 2014, Pharmaceutical Research.

[16]  Sandeep Yadav,et al.  Specific interactions in high concentration antibody solutions resulting in high viscosity. , 2010, Journal of pharmaceutical sciences.

[17]  Steven J Shire,et al.  Reversible self-association increases the viscosity of a concentrated monoclonal antibody in aqueous solution. , 2005, Journal of pharmaceutical sciences.

[18]  Steven J Shire,et al.  Reversible self-association of a concentrated monoclonal antibody solution mediated by Fab-Fab interaction that impacts solution viscosity. , 2008, Journal of pharmaceutical sciences.

[19]  Steven J. Shire,et al.  Establishing a Link Between Amino Acid Sequences and Self-Associating and Viscoelastic Behavior of Two Closely Related Monoclonal Antibodies , 2011, Pharmaceutical Research.

[20]  Oliver F. Lange,et al.  Full correlation analysis of conformational protein dynamics , 2007, Proteins.

[21]  Bernhardt L Trout,et al.  Computational methods to predict therapeutic protein aggregation. , 2012, Methods in molecular biology.

[22]  T. Singer,et al.  Minipig as a potential translatable model for monoclonal antibody pharmacokinetics after intravenous and subcutaneous administration , 2012, mAbs.

[23]  Sandeep Yadav,et al.  Weak interactions govern the viscosity of concentrated antibody solutions: high-throughput analysis using the diffusion interaction parameter. , 2012, Biophysical journal.

[24]  Nia Beckley,et al.  Isomerization of Asp-Asp motif in model peptides and a monoclonal antibody Fab fragment. , 2013, Journal of pharmaceutical sciences.

[25]  Randal R. Ketchem,et al.  Entanglement model of antibody viscosity. , 2014, The journal of physical chemistry. B.

[26]  Ronald T Borchardt,et al.  Aspartate isomerization in the complementarity-determining regions of two closely related monoclonal antibodies. , 2007, Biochemistry.

[27]  Herren Wu,et al.  Development of motavizumab, an ultra-potent antibody for the prevention of respiratory syncytial virus infection in the upper and lower respiratory tract. , 2007, Journal of molecular biology.

[28]  D. Eisenberg,et al.  Analysis of membrane and surface protein sequences with the hydrophobic moment plot. , 1984, Journal of molecular biology.

[29]  Lisa J. Bernstein,et al.  A strategy for risk mitigation of antibodies with fast clearance , 2012, mAbs.

[30]  Bernhardt L. Trout,et al.  Design of therapeutic proteins with enhanced stability , 2009, Proceedings of the National Academy of Sciences.