Predicting pKa for proteins using COSMO-RS

We have used the COSMO-RS implicit solvation method to calculate the equilibrium constants, pKa, for deprotonation of the acidic residues of the ovomucoid inhibitor protein, OMTKY3. The root mean square error for comparison with experimental data is only 0.5 pH units and the maximum error 0.8 pH units. The results show that the accuracy of pKa prediction using COSMO-RS is as good for large biomolecules as it is for smaller inorganic and organic acids and that the method compares very well to previous pKa predictions of the OMTKY3 protein using Quantum Mechanics/Molecular Mechanics. Our approach works well for systems of about 1000 atoms or less, which makes it useful for small proteins as well as for investigating portions of larger proteins such as active sites in enzymes.

[1]  Wei Wang,et al.  External Factors Affecting Protein Aggregation , 2010 .

[2]  Sarah L. Williams,et al.  Progress in the prediction of pKa values in proteins , 2011, Proteins.

[3]  M. Olsson,et al.  Improving the desolvation penalty in empirical protein pKa modeling , 2012, Journal of Molecular Modeling.

[4]  Hans W. Horn,et al.  ELECTRONIC STRUCTURE CALCULATIONS ON WORKSTATION COMPUTERS: THE PROGRAM SYSTEM TURBOMOLE , 1989 .

[5]  Eamonn F. Healy,et al.  Development and use of quantum mechanical molecular models. 76. AM1: a new general purpose quantum mechanical molecular model , 1985 .

[6]  A. Bonincontro,et al.  Lysozyme binds onto functionalized carbon nanotubes. , 2013, Colloids and surfaces. B, Biointerfaces.

[7]  A. Klamt,et al.  Fast solvent screening via quantum chemistry: COSMO‐RS approach , 2002 .

[8]  P. Harbury,et al.  Tanford-Kirkwood electrostatics for protein modeling. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[9]  F Guarnieri,et al.  A self-consistent, microenvironment modulated screened coulomb potential approximation to calculate pH-dependent electrostatic effects in proteins. , 1999, Biophysical journal.

[10]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[11]  A. D. Robertson,et al.  pH, ionic strength, and temperature dependences of ionization equilibria for the carboxyl groups in turkey ovomucoid third domain. , 1995, Biochemistry.

[12]  G. Vriend,et al.  Optimizing the hydrogen‐bond network in Poisson–Boltzmann equation‐based pKa calculations , 2001, Proteins.

[13]  G Vriend,et al.  WHAT IF: a molecular modeling and drug design program. , 1990, Journal of molecular graphics.

[14]  A. Warshel,et al.  Electrostatic contributions to protein–protein binding affinities: Application to Rap/Raf interaction , 1998, Proteins.

[15]  J. Perdew,et al.  Density-functional approximation for the correlation energy of the inhomogeneous electron gas. , 1986, Physical review. B, Condensed matter.

[16]  F. Weigend,et al.  Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. , 2005, Physical chemistry chemical physics : PCCP.

[17]  Martin Korth,et al.  Third-Generation Hydrogen-Bonding Corrections for Semiempirical QM Methods and Force Fields , 2010 .

[18]  M K Gilson,et al.  Theoretical and experimental analysis of ionization equilibria in ovomucoid third domain. , 1998, Biochemistry.

[19]  Andreas Klamt,et al.  Accurate prediction of basicity in aqueous solution with COSMO‐RS , 2006, J. Comput. Chem..

[20]  Martin Korth,et al.  Comparison of molecular mechanics, semi-empirical quantum mechanical, and density functional theory methods for scoring protein-ligand interactions. , 2013, The journal of physical chemistry. B.

[21]  Jan H. Jensen,et al.  The determinants of carboxyl pKa values in turkey ovomucoid third domain , 2004, Proteins.

[22]  R. Norel,et al.  Electrostatic aspects of protein-protein interactions. , 2000, Current opinion in structural biology.

[23]  S. Hovmöller,et al.  The average atomic volume and density of proteins , 1998 .

[24]  Nathan A. Baker,et al.  PROGRESS IN THE PREDICTION OF pK a VALUES IN PROTEINS , 2011 .

[25]  R. Huber,et al.  X‐ray crystal structure of the complex of human leukocyte elastase (PMN elastase) and the third domain of the turkey ovomucoid inhibitor. , 1986, The EMBO journal.

[26]  Andreas Klamt,et al.  First Principles Calculations of Aqueous pKa Values for Organic and Inorganic Acids Using COSMO-RS Reveal an Inconsistency in the Slope of the pKa Scale. , 2003, The journal of physical chemistry. A.

[27]  Hans W. Horn,et al.  Fully optimized contracted Gaussian basis sets for atoms Li to Kr , 1992 .