Using DelPhi capabilities to mimic protein's conformational reorganization with amino acid specific dielectric constants.

Many molecular events are associated with small or large conformational changes occurring in the corresponding proteins. Modeling such changes is a challenge and requires significant amount of computing time. From point of view of electrostatics, these changes can be viewed as a reorganization of local charges and dipoles in response to the changes of the electrostatic field, if the cause is insertion or deletion of a charged amino acid. Here we report a large scale investigation of modeling the changes of the folding energy due to single mutations involving charged group. This allows the changes of the folding energy to be considered mostly electrostatics in origin and to be calculated with DelPhi assigning residue-specific value of the internal dielectric constant of protein. The predicted energy changes are benchmarked against experimentally measured changes of the folding energy on a set of 257 single mutations. The best fit between experimental values and predicted changes is used to find out the effective value of the internal dielectric constant for each type of amino acid. The predicted folding free energy changes with the optimal, amino acid specific, dielectric constants are within RMSD=0.86 kcal/mol from experimentally measured changes.

[1]  C. Kieslich,et al.  An evaluation of Poisson-Boltzmann electrostatic free energy calculations through comparison with experimental mutagenesis data. , 2011, Biopolymers.

[2]  L. R. Scott,et al.  Electrostatics and diffusion of molecules in solution: simulations with the University of Houston Brownian dynamics program , 1995 .

[3]  Xueyu Song,et al.  An Inhomogeneous Model of Protein Dielectric Properties: Intrinsic Polarizabilities of Amino Acids , 2002 .

[4]  Zhe Zhang,et al.  In Silico and In Vitro Investigations of the Mutability of Disease-Causing Missense Mutation Sites in Spermine Synthase , 2011, PloS one.

[5]  Emil Alexov,et al.  Poisson-Boltzmann calculations of nonspecific salt effects on protein-protein binding free energies. , 2007, Biophysical journal.

[6]  K A Dill,et al.  Explicit-water molecular dynamics study of a short-chain 3,3 ionene in solutions with sodium halides. , 2009, Journal of Chemical Physics.

[7]  S J Wodak,et al.  Calculations of electrostatic properties in proteins. Analysis of contributions from induced protein dipoles. , 1987, Journal of molecular biology.

[8]  An-Suei Yang,et al.  Electrostatic contributions to the binding free energy of the lambdacI repressor to DNA. , 1998, Biophysical journal.

[9]  Wilfred F van Gunsteren,et al.  Explicit-solvent molecular dynamics simulations of the polysaccharide schizophyllan in water. , 2007, Biophysical journal.

[10]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[11]  Emil Alexov,et al.  Rapid grid‐based construction of the molecular surface and the use of induced surface charge to calculate reaction field energies: Applications to the molecular systems and geometric objects , 2002, J. Comput. Chem..

[12]  Akinori Sarai,et al.  ProTherm, version 4.0: thermodynamic database for proteins and mutants , 2004, Nucleic Acids Res..

[13]  Arieh Warshel,et al.  A surface constrained all‐atom solvent model for effective simulations of polar solutions , 1989 .

[14]  E. Padan,et al.  Multiconformation continuum electrostatics analysis of the NhaA Na+/H+ antiporter of Escherichia coli with functional implications. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[15]  W. C. Still,et al.  Semianalytical treatment of solvation for molecular mechanics and dynamics , 1990 .

[16]  B Honig,et al.  On the pH dependence of protein stability. , 1993, Journal of molecular biology.

[17]  M K Gilson,et al.  The dielectric constant of a folded protein , 1986, Biopolymers.

[18]  M. Sanner,et al.  Reduced surface: an efficient way to compute molecular surfaces. , 1996, Biopolymers.

[19]  Barry Honig,et al.  Extending the Applicability of the Nonlinear Poisson−Boltzmann Equation: Multiple Dielectric Constants and Multivalent Ions† , 2001 .

[20]  Junmei Wang,et al.  Junmei Wang, Romain M. Wolf, James W. Caldwell, Peter A. Kollman, and David A. Case, "Development and testing of a general amber force field"Journal of Computational Chemistry(2004) 25(9) 1157–1174 , 2005, J. Comput. Chem..

[21]  F M Richards,et al.  Areas, volumes, packing and protein structure. , 1977, Annual review of biophysics and bioengineering.

[22]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

[23]  Emil Alexov,et al.  Role of the protein side‐chain fluctuations on the strength of pair‐wise electrostatic interactions: Comparing experimental with computed pKas , 2002, Proteins.

[24]  A. Karshikoff,et al.  A model of a local dielectric constant in proteins , 1998 .

[25]  J. Schlessman,et al.  Electrostatic effects in a network of polar and ionizable groups in staphylococcal nuclease. , 2008, Journal of molecular biology.

[26]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[27]  A. Wada,et al.  A theoretical study of the dielectric constant of protein. , 1988, Protein engineering.

[28]  R Nussinov,et al.  Explicit and implicit water simulations of a β‐hairpin peptide , 1999, Proteins.

[29]  Justin R. Spaeth,et al.  A comparison of implicit- and explicit-solvent simulations of self-assembly in block copolymer and solute systems. , 2011, The Journal of chemical physics.

[30]  Huan-Xiang Zhou,et al.  Comparison of calculation and experiment implicates significant electrostatic contributions to the binding stability of barnase and barstar. , 2003, Biophysical journal.

[31]  M. Gilson,et al.  Prediction of pH-dependent properties of proteins. , 1994, Journal of molecular biology.

[32]  Huan-Xiang Zhou,et al.  Electrostatic contributions to T4 lysozyme stability: solvent-exposed charges versus semi-buried salt bridges. , 2002, Biophysical journal.

[33]  J. B. Matthew.  Electrostatic effects in proteins. , 1985, Annual review of biophysics and biophysical chemistry.

[34]  M. Perutz.  Electrostatic effects in proteins. , 1978, Science.

[35]  Kelly K. Lee,et al.  Distance dependence and salt sensitivity of pairwise, coulombic interactions in a protein , 2002, Protein science : a publication of the Protein Society.

[36]  Martin Karplus,et al.  pH-Dependence of Protein Stability: Absolute Electrostatic Free Energy Differences between Conformations† , 1997 .

[37]  Mark A Olson,et al.  Modeling loop reorganization free energies of acetylcholinesterase: A comparison of explicit and implicit solvent models , 2004, Proteins.

[38]  A. Warshel,et al.  Calculations of electrostatic interactions in biological systems and in solutions , 1984, Quarterly Reviews of Biophysics.

[39]  B. Honig,et al.  Classical electrostatics in biology and chemistry. , 1995, Science.

[40]  Electrostatic contributions to the binding free energy of the lambdacI repressor to DNA. , 1998, Biophysical journal.

[41]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[42]  D. Porschke.  Electrostatics and electrodynamics of bacteriorhodopsin. , 1996, Biophysical journal.

[43]  Emil Alexov,et al.  On the electrostatic component of protein-protein binding free energy , 2008, PMC biophysics.

[44]  Charles R. Cantor,et al.  Annual Review of Biophysics and Bioengineering , 1972 .

[45]  Emil Alexov,et al.  Numerical calculations of the pH of maximal protein stability. The effect of the sequence composition and three-dimensional structure. , 2003, European journal of biochemistry.

[46]  D. D. Yue,et al.  Theory of Electric Polarization , 1974 .

[47]  E. Alexov,et al.  Calculated protein and proton motions coupled to electron transfer: electron transfer from QA- to QB in bacterial photosynthetic reaction centers. , 1999, Biochemistry.

[48]  K. Sharp,et al.  Electrostatic interactions in macromolecules: theory and applications. , 1990, Annual review of biophysics and biophysical chemistry.

[49]  S. Hassan,et al.  A critical analysis of continuum electrostatics: The screened Coulomb potential–implicit solvent model and the study of the alanine dipeptide and discrimination of misfolded structures of proteins , 2002, Proteins.

[50]  V. Barone,et al.  A hybrid explicit/implicit solvation method for first-principle molecular dynamics simulations. , 2008, The Journal of chemical physics.

[51]  J. Trylska,et al.  Continuum molecular electrostatics, salt effects, and counterion binding—A review of the Poisson–Boltzmann theory and its modifications , 2008, Biopolymers.

[52]  Huan-Xiang Zhou,et al.  Do electrostatic interactions destabilize protein-nucleic acid binding? , 2007, Biopolymers.

[53]  Emil Alexov,et al.  A missense mutation in CLIC2 associated with intellectual disability is predicted by in silico modeling to affect protein stability and dynamics , 2011, Proteins.

[54]  Mark A Olson,et al.  An efficient hybrid explicit/implicit solvent method for biomolecular simulations , 2004, J. Comput. Chem..

[55]  Alan E. Mark,et al.  Dielectric properties of trypsin inhibitor and lysozyme calculated from molecular dynamics simulations , 1993 .

[56]  D. Bashford,et al.  Electrostatic coupling to pH‐titrating sites as a source of cooperativity in protein‐ligand binding , 1998, Protein science : a publication of the Protein Society.

[57]  Thomas Simonson,et al.  Dielectric relaxation in proteins: Microscopic and macroscopic models , 1999 .

[58]  Zhe Zhang,et al.  Computational analysis of missense mutations causing Snyder‐Robinson syndrome , 2010, Human mutation.

[59]  Arieh Warshel,et al.  Microscopic simulations of macroscopic dielectric constants of solvated proteins , 1991 .