Electrostatic effects in macromolecules: Fundamental concepts and practical modeling

Current Opinion in Structural Biology

Volumn 8, Issue 2, 1998, Pages 211-217

  Warshel, Arieh ^a   Papazyan, Arno ^a  

^a UNIVERSITY OF SOUTHERN CALIFORNIA   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

DIELECTRIC CONSTANT; ELECTRIC FIELD; ELECTRICITY; MACROMOLECULE; MODEL; OXIDATION REDUCTION POTENTIAL; PKA; PRIORITY JOURNAL; REVIEW;

EID: 0032054517     PISSN: 0959440X     EISSN: None     Source Type: Journal    
DOI: 10.1016/S0959-440X(98)80041-9     Document Type: Article

References (85)

1. Warshel A, Russell ST. Calculations of electrostatic interactions in biological systems and in solutions. Q Rev Biophys. 17:1984;283-422.
2. Nakamura H. Roles of electrostatic interaction in proteins. Q Rev Biophys. 29:1996;1-90.
3. Warshel A. Electrostatic basis of structure - function correlation in proteins. Accounts Chem Res. 14:1981;284-290.
4. Sharp KA, Honig B. Electrostatic interactions in macromolecules: theory and applications. Annu Rev Biophys Biophys Chem. 19:1990;301-332.
5. Warshel A, Åqvist J. Electrostatic energy and macromolecular function. Annu Rev Biophys Biophys Chem. 20:1991;267-298.
6. Gilson MK. Theory of electrostatic interactions in macromolecules. Curr Opin Struct Biol. 5:1995;216-223.
7. Schaefer M, Sommer M, Karplus M. pH-dependence of protein stability: absolute electrostatic free energy differences between conformations. J Phys Chem B. 101:1997;1663-1683.
8. Dufner H, Kast SM, Brickmann J, Schlenkrich M. Ewald summation versus direct summation of shifted-force potentials for the calculation of electrostatic interactions in solids - a quantitative study. J Comput Chem. 18:1997;660-676.
9. Figueirido F, Delbuono GS, Levy RM. On finite-size corrections to the free-energy of ionic hydration. J Phys Chem B. 101:1997;5622-5623.
10. Hummer G, Pratt LR, Garcia AE. Free energy of ionic hydration. J Phys Chem. 100:1996;1206-1215.
11. Lee FS, Warshel A. A local reaction field method for fast evaluation of long-range electrostatic interactions in molecular simulations. J Chem Phys. 97:1992;3100-3107.
12. Aqvist J. Calculation of absolute binding free-energies for charged ligands and effects of long-range electrostatic interactions. of special interest J Comput Chem. 17:1996;1587-1597 The absolute binding free energies of charged ligands are evaluated using an innovative variant of the LRA treatment of [15] and the LRF treatment [11]. This LRA variant neglects the average over the configuration generated with the uncharged solute and estimates the hydrophobic contribution by scaling the average van der Waals contribution (this contribution is evaluated in [15] by the PDLD surface area estimate). This paper pays significant attention to the importance of the proper treatment of long-range effects, clearly illustrating the superiority of the LRF treatment over cut-off methods. The problems associated with FEP approaches are also illustrated clearly.
13. Warshel A, Sussman F, King G. Free energy of charges in solvated proteins: microscopic calculations using a reversible charging process. Biochemistry. 25:1986;8368-8372.
14. Kollman P. Free-energy calculations - applications to chemical and biochemical phenomena. Chem Rev. 93:1993;2395-2417.
15. Lee FS, Chu ZT, Bolger MB, Warshel A. Calculations of antibody antigen interactions - microscopic and semimicroscopic evaluation of the free-energies of binding of phosphorylcholine analogs to McPC603. Protein Eng. 5:1992;215-228.
16. as of ionizable residues in proteins - an explicit solvent calculation for lysozyme. Proteins. 20:1994;85-97.
17. Aqvist J, Medina C, Samuelsson J-E. New method for predicting binding affinity in computer-aided drug design. Protein Eng. 7:1994;385-391.
18. Liu Y, Ichiye T. The static dielectric-constant of the soft sticky dipole model of liquid water - Monte-Carlo simulation. Chem Phys Lett. 256:1996;334-340.
19. Coalson RD, Duncan A. Statistical mechanics of a multipolar gas: a lattice field theory approach. J Phys Chem. 100:1996;2612-2620.
20. Papazyan A, Warshel A. Continuum and dipole-lattice models of solvation. of special interest J Phys Chem B. 101:1997;11254-11264 The connection between discrete and continuum representations of solvation is made in a concrete manner. New insights into the meaning and domain of applicability of the Clausius - Mossotti equation are revealed. The insufficiency of capturing the macroscopic dielectric constant for solvation behavior is demonstrated, with implications for cut-off schemes in molecular simulations. Further validation of the Langevin dipoles model is provided by its essentially exact agreement with the more 'realistic' Brownian dipole model at a given dielectric constant. A novel way of obtaining the 'microscopic' dielectric constants is devised, which involves a pointwise comparison of microscopic electric polarization and potential with the numerical solution of the Poisson equation. The fact that the continuum dielectric is the infinite dipole density limit of a more general dipole-lattice description is shown by recovering the continuum results with dipole lattices of high number density. The linearity of the continuum models is explained as a direct consequence of being the infinite density limit of a dipolar model.
21. eff for protein interiors.
22. Papazyan A, Warshel A. A stringent test of the cavity concept in continuum dielectrics. J Chem Phys. 107:1997;7975-7978.
23. Warshel A, Naray-Szabo G, Sussman F, Hwang J-K. How do serine proteases work? Biochemistry. 28:1989;3629-3638.
24. Warshel A, Russell ST, Churg AK. Macroscopic models for studies of electrostatic interactions in proteins: limitations and applicability. Proc Natl Acad Sci USA. 81:1984;4785-4789.
25. You TJ, Bashford D. Conformation and hydrogen-ion titration of proteins - a continuum electrostatic model with conformational flexibility. Biophys J. 69:1995;1721-1733.
26. eff. At any rate, the null model is not really a predictive model in the sense that it only 'predicts' the trivial cases in which the protein does very little; as argued in the text and [24], these are not the interesting cases.
27. St WC II, Tempozyk A, Hawley RC, Hendrickson T. Semianalytical treatment of solvation for molecular mechanics and dynamics. J Am Chem Soc. 112:1990;6127-6129.
28. Edinger SR, Cortis C, Shenkin PS, Friesner RA. Solvation free energies of peptides: comparison of approximate continuum solvation models with accurate solution of the Poisson - Boltzmann equation. of special interest J Phys Chem B. 101:1997;1190-1197 An attempt at a systematic comparison of several simplified models of solvation. Readers must be warned, however, about a couple of shortcomings: it singles out the GB model as an approximation for the PB equation or a solution of the multicavity continuum model, implying that the other models are not. In fact, none of the compared models are rigorously derived from the PB equation, and all of them can be regarded as an 'approximation' to the PB equation. However, a good solvation model does not necessarily aim to approximate PB equation which is itself an approximation to the real solvation behavior. GB is the mathematical statement of the empirically observed behavior of a pair of charges in homogeneous solutions. As such, it is a good solvation model within its domain of applicability, namely uniform dielectrics. Despite the implication in this paper, GB is no more valid than a distance-dependent dielectric (DDD) model, when the latter is properly parametrized based on empirical data (just as GB is parametrized). GB and DDD models have been shown to be equivalent [1,41,85].
29. Warshel A. Calculations of chemical processes in solutions. J Phys Chem. 83:1979;1640-16542.
30. Figueirido F, Delbuono GS, Levy RM. On finite-size effects in computer simulations using the Ewald potential. J Chem Phys. 103:1995;6133-6142.
31. Åqvist J, Hansson T. On the validity of electrostatic linear-response in polar solvents. J Phys Chem. 100:1996;9512-9521.
32. Orozco M, Luque FJ. Generalized linear-response approximation in discrete methods. Chem Phys Lett. 265:1997;473-480.
33. Florian J, Warshel A. Langevin dipoles model for ab initio calculations of chemical processes in solution: parametrization and application to hydration free energies of neutral and ionic solutes and conformational analysis in aqueous solution. J Phys Chem B. 101:1997;5583-5595.
34. Malcolm NOJ, McDouall JJW. Assessment of the Langevin dipoles solvation model for Hartree - Fock wave-functions. J Mol Struct. 366:1996;1-9.
35. Marten B, Kim K, Cortis C, Friesner RA, Murphy RB, Ringnalda MN, Sitkoff D, Honig B. New model for calculation of solvation free energies: correction of self-consistent reaction field continuum dielectric theory for short-range hydrogen-bonding effects. J Phys Chem. 100:1996;11775-11788.
36. Tomasi J, Mennucci B, Cammi R, Cossi M. Quantum mechanical models for reactions in solution. Náray-Szabó G, Warshel A. Computational Approaches to Biochemical Reactivity. 1997;1-102 Kluwer Academic Publishers, Dordrecht, The Netherlands.
37. Hawkins GD, Cramer CJ, Truhlar DG. Parametrized models of aqueous free-energies of solvation based on pairwise descreening of solute atomic charges from a dielectric medium. J Phys Chem. 100:1996;19824-19839.
38. Fothergill M, Goodman MF, Petruska J, Warshel A. Structure-energy analysis of the role of metal ion in phosphodiester bond hydrolysis by DNA polymerase I. J Am Chem Soc. 117:1995;11619-11627.
39. Åqvist J, Fothergill M. Computer simulation of the trisephosphate isomerase catalyzed reaction. J Biol Chem. 271:1996;10010-10016.
40. Saito M. Molecular dynamics free-energy study of a protein in solution with all degrees of freedom and long-range Coulomb interactions. J Phys Chem. 99:1995;17043-17048.
41. Mehler EL. Self-consistent, free energy based approximation to calculate pH dependent electrostatic effects in proteins. of special interest J Phys Chem. 100:1996;16006-16018 This paper outlines the relationship between the Lorentz - Debye - Sack theory of solvation and a DDD model. It shows that the latter is equivalent to the GB method.
42. as, but this is due to the fact that the method results in a large effective dielectric constant for charge - charge interactions. In this sense, the use of a large dielectric constant is not much different than the approach used in [24]; however, here and in many related studies, the crucial self-energy (i.e. solvation) contribution is overlooked and charge - charge interactions are mistakenly assumed to be the main electrostatic effects. This along with the assumption that insensitivity to ionic strength implies a diminished importance of electrostatics are two important misconceptions that still surface in current literature.
43. Lamm G, Pack GR. Calculation of dielectric constants near polyelectrolytes in solution. J Phys Chem B. 101:1997;959-965.
44. Svensson B, Jonsson B. An efficient simulation technique for electrostatic free-energies with applications to azurin. J Comput Chem. 16:1995;370-377.
45. a values of lysine residues in calbindin d-9k. J Mol Biol. 259:1996;828-839.
46. a calculations.
47. Zhou H-X. Control of reduction potential by protein matrix: lesson from a spherical protein model. J Biol Inorg Chem. 2:1997;109-113.
48. Bertini I, Gori-Savellini G, Luchinat C. Are unit charges always negligible? J Biol Inorg Chem. 2:1997;114-118.
49. Mauk AG, Moore GR. Control of metalloprotein redox potentials: what does site-directed mutagenesis of hemoproteins tell us? J Biol Inorg Chem. 2:1997;119-125.
50. Gunner MR, Alexov E, Torres E, Lipovaca S. The importance of the protein in controlling the electrochemistry of heme metalloproteins: methods of calculation and analysis. J Biol Inorg Chem. 2:1997;126-134.
51. Náray-Szabó G. Electrostatic modulation of electron transfer in the active site of heme peroxidases. J Biol Inorg Chem. 2:1997;135-138.
52. Armstrong FA. Evaluations of reduction potential data in relation to coupling, kinetics and function. J Biol Inorg Chem. 2:1997;139-142.
53. Warshel A, Papazyan A, Muegge I. Microscopic and semimacroscopic redox calculations: what can and cannot be learned from continuum models. J Biol Inorg Chem. 2:1997;143-152.
54. Muegge I, Qi PX, Wand AJ, Chu ZT, Warshel A. The reorganization energy of cytochrome-c revisited. of special interest J Phys Chem B. 101:1997;825-836 The experimentally observed structural changes upon oxidation of cytochrome c are coverted to reorganization energy using different microscopic and macroscopic models. The results provide a structure based estimate of the activation barrier for electron transfer between two cytochromes, and establish that different dielectric constants are associated with different models.
55. Kassner RJ. Effects of nonpolar environments on the redox potentials of heme proteins. Proc Natl Acad Sci USA. 69:1972;2263-2267.
56. Churg AK, Warshel A. Control of redox potential of cytochrome c and microscopic dielectric effects in proteins. Biochemistry. 25:1986;1675-1681.
57. Stephens PJ, Jollie DR, Warshel A. Protein control of redox potentials of iron-sulfur proteins. Chem Rev. 96:1996;2491-2513.
58. Swartz PD, Beck BW, Ichiye T. Structural origins of redox potentials in Fe-S proteins - electrostatic potentials of crystal structures. Biophys J. 71:1996;2958-2969.
59. Alden RG, Parson WW, Chu ZT, Warshel A. Calculations of electrostatic energies in photosynthetic reaction centers. J Am Chem Soc. 117:1995;12284-12298.
60. Yelle RB, Ichiye T. Solvation free-energy reaction curves for electron-transfer in aqueous solution - theory and simulation. J Phys Chem B. 101:1997;4127-4135.
61. Merz KM, Banci L. Binding of azide to human carbonic anydrase II. The role electrostatic complementarity plays in selecting the preferred resonance structure of azide. J Phys Chem. 100:1996;17414-17420.
62. Madura JD, Nakajima Y, Hamilton RM, Wierzbicki A, Warshel A. Calculations of the electrostatic free-energy contributions to the binding free-energy of sulfonamides to carbonic-anhydrase. Struct Chem. 7:1996;131-138.
63. Muegge I, Schweins T, Langer R, Warshel A. Electrostatic control of GTP and GDP binding in the oncoprotein P21(ras). Structure. 4:1996;475-489.
64. Muegge I, Tao H, Warshel A. A fast estimate of electrostatic group contributions to the free energy of protein-inhibitor binding. Protein Eng. 10:1998;1363-1372.
65. Warshel A, Tao H, Fothergill M, Chu ZT. Effective methods for estimation of binding energies in computer-aided drug design. Isr J Chem. 34:1994;253-256.
66. Froloff N, Windemuth A, Honig B. On the calculation of binding free-energies using continuum methods - application to MHC class-1 protein - peptide interactions. of special interest Protein Sci. 6:1997;1293-1301 Binding free energies of peptides to proteins are evaluated by the Poisson - Boltzmann method, using the charging cycle introduced in earlier PDLD/S studies [15]; however, current results are still not satisfactory both in terms of the observed trends and absolute values.
67. Novotny J, Bruccoleri RE, Davis M, Sharp KA. Empirical free energy calculations - a blind test and further improvements to the method. J Mol Biol. 268:1997;401-411.
68. Muegge I, Schweins T, Warshel A. Electrostatic contributions to protein - protein binding affinities: application to Rap/Raf interaction. Proteins. 30:1998;407-423.
69. Mccoy AJ, Epa VC, Colman PM. Electrostatic complementarity at protein/protein interfaces. J Mol Biol. 268:1997;570-584.
70. Náray-Szabó G, Fuxreiter M, Warshel A. Electrostatic basis of enzyme catalysis. Náray-Szabó G, Warshel A. Computational Approaches to Biochemical Reactivity. 1997;237-294 Kluwer Academic Publishers, Dordrecht, The Netherlands.
71. Van Beek J, Callender R, Gunner MR. The contribution of electrostatic and van-der-Waals interactions to the stereospecificity of the reaction catalyzed by lactate-dehydrogenase. Biophys J. 72:1997;619-626.
72. Fisher CL, Cabelli DE, Hallewell RA, Beroza P, Lo TP, Getzoff ED, Tainer JA. Computational, pulse-radiolytic, and structural investigations of lysine-136 and its role in the electrostatic triad of human Cu, Zn superoxide dismutase. Proteins. 29:1997;103-112.
73. Antosiewicz J, McCammon JA, Wlodek ST, Gilson MK. Simulation of charge-mutant acetylcholinesterases. Biochemistry. 34:1995;4211-4219.
74. Fuxreiter M, Warshel A. On the origin of the catalytic power of acetylcholinesterase: computer simulation studies. J Am Chem Soc. 120:1997;183-194.
75. + through the gramicidin channel. This is done by comparing the calculated and observed conductance. It is found that the barriers calculated by the approach of Roux and Karplus [76] drastically overestimate the actual barrier whereas the profile obtained by Åqvist and Warshel [77] is more realistic. Although the authors are hesitant to criticize [76] and point out that a scaling factor of ~0.3 gives 'excellent' agreement with experiments, the fact remains that seemingly rigorous potential-of-mean-force force treatments, when used with an improper treatment of long-range interactions, produces incorrect penetration barriers (which are the most important factors in controlling ion conductance).
76. Roux B, Karplus M. Ion transport in the gramicidin channel: free energy of the solvated right-handed dimer in a model membrane. J Am Chem Soc. 115:1993;3250-3262.
77. + by gramicidin A. Biophys J. 56:1989;171-182.
78. Dorman V, Partenskii MB, Jordan PC. A semimicroscopic Monte Carlo study of permeation energetics in a gramicidin-like channel - the origin of selectivity. Biophys J. 70:1996;121-134.
79. Ben Tal N, Honig B, Peitzsch RM, Denisov G. Binding of small basic peptides to membranes containing acidic lipids - theoretical models and experimental results. Biophys J. 71:1996;561-575.
80. Sheng Q, Schulten K, Pidgeon C. Molecular dynamics simulation of immobilized artificial membranes. J Phys Chem. 99:1995;11018-11027.
81. King G, Lee FS, Warshel A. Microscopic simulations of macroscopic dielectric constants of solvated proteins. J Chem Phys. 95:1991;4366-4377.
82. Simonson T, Brooks CL. Charge screening and the dielectric-constant of proteins - insights from molecular dynamics. J Am Chem Soc. 118:1996;8452-8458.
83. Sham YY, Muegge I, Warshel A. The effect of protein relaxation on charge - charge interactions and dielectric constants in proteins. Biophys J. 1998;. in press.
84. Alexov EG, Gunner MR. Incorporating protein conformational flexibility into the calculation of pH-dependent protein properties. Biophys J. 72:1997;2075-2093.
85. Luzhkov V, Warshel A. Microscopic models for quantum mechanical calculations of chemical processes in solution. J Comput Chem. 13:1992;199-213.