Electrostatics calculations: latest methodological advances

Current Opinion in Structural Biology

Volumn 16, Issue 2, 2006, Pages 142-151

  Koehl, Patrice ^a  

^a Davis   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

SOLVENT;

CALCULATION; DIELECTRIC CONSTANT; ELECTRICITY; MATHEMATICAL COMPUTING; METHODOLOGY; MOLECULAR DYNAMICS; PRIORITY JOURNAL; REVIEW; THEORETICAL MODEL;

ALGORITHMS; COMPUTATIONAL BIOLOGY; COMPUTER SIMULATION; ELECTROSTATICS; MOLECULAR CONFORMATION; PROTEINS; SOLVENTS; THERMODYNAMICS;

EID: 33645981305     PISSN: 0959440X     EISSN: None     Source Type: Journal    
DOI: 10.1016/j.sbi.2006.03.001     Document Type: Review

References (87)

1. Kauzmann W. Some factors in the interpretation of protein denaturation. Adv Protein Chem 14 (1959) 1-63
2. Perutz M.F. Electrostatic effects in proteins. Science 201 (1978) 1187-1191
3. Neves-Petersen M.T., and Petersen S.B. Protein electrostatics: a review of the equations and methods used to model electrostatic equations in biomolecules-applications in biotechnology. Biotechnol Annu Rev 9 (2003) 315-395
4. Tsui V., and Case D.A. Theory and applications of the generalized Born solvation model in macromolecular simulations. Biopolymers 56 (2000) 275-291
5. Bashford D., and Case D.A. Generalized born models of macromolecular solvation effects. Annu Rev Phys Chem 51 (2000) 129-152
6. Simonson T. Macromolecular electrostatics: continuum models and their growing pains. Curr Opin Struct Biol 11 (2001) 243-252
7. Feig M., and Brooks III C.L. Recent advances in the development and application of implicit solvent models in biomolecule simulations. Curr Opin Struct Biol 14 (2004) 217-224
8. Baker N.A. Poisson-Boltzmann methods for biomolecular electrostatics. Methods Enzymol 383 (2004) 94-118
9. Baker N.A. Improving implicit solvent simulations: a Poisson-centric view. Curr Opin Struct Biol 15 (2005) 137-143. This review of implicit solvent models surveys recent methodological advances in solving the PB equation.
10. Sagui C., and Darden T.A. Molecular dynamics simulations of biomolecules: long-range electrostatic effects. Annu Rev Biophys Biomol Struct 28 (1999) 155-179
11. Ponder J.W., and Case D.A. Force fields for protein simulations. Adv Protein Chem 66 (2003) 27-85
12. Mackerell Jr. A.D. Empirical force fields for biological macromolecules: overview and issues. J Comput Chem 25 (2004) 1584-1604
13. Tobias D.J. Electrostatics calculations: recent methodological advances and applications to membranes. Curr Opin Struct Biol 11 (2001) 253-261
14. Brannigan G., Lin L.C., and Brown F.L. Implicit solvent simulation models for biomembranes. Eur Biophys J 35 (2005) 104-124. This review of implicit solvent models focuses on calculations relevant to membrane systems. It surveys coarse-grained representations of membranes and describes simulation techniques adapted to these approximate representations, in which the solvent is considered implicitly.
15. Gibas C.J., Jambeck P., and Subramaniam S. Continuum electrostatic methods applied to pH-dependent properties of antibody-antigen association. Methods 20 (2000) 292-309
16. Sheinerman F.B., Norel R., and Honig B. Electrostatic aspects of protein-protein interactions. Curr Opin Struct Biol 10 (2000) 153-159
17. Wang W., Donini O., Reyes C.M., and Kollman P.A. Biomolecular simulations: recent developments in force fields, simulations of enzyme catalysis, protein-ligand, protein-protein, and protein-nucleic acid noncovalent interactions. Annu Rev Biophys Biomol Struct 30 (2001) 211-243
18. Sinha N., and Smith-Gill S.J. Electrostatics in protein binding and function. Curr Protein Pept Sci 3 (2002) 601-614
19. Vizcarra C.L., and Mayo S.L. Electrostatics in computational protein design. Curr Opin Chem Biol 9 (2005) 622-626. The complexity and computational cost of in silico protein design have led to the development of many approximate methods for computing electrostatics, focusing on pairwise decomposable models. These methods are reviewed in this paper, with emphasis on implicit solvent models.
20. Kumar S., and Nussinov R. Close-range electrostatic interactions in proteins. ChemBioChem 3 (2002) 604-617
21. Kumar S., and Nussinov R. Different roles of electrostatics in heat and in cold: adaptation by citrate synthase. ChemBioChem 5 (2004) 280-290
22. Donchev A.G., Ozrin V.D., Subbotin M.V., Tarasov O.V., and Tarasov V.I. A quantum mechanical polarizable force field for biomolecular interactions. Proc Natl Acad Sci USA 102 (2005) 7829-7834. This paper introduces a quantum mechanical polarizable force-field that is optimized for drug design. It is shown to agree well with experimental data in simulations of methane and methanol in gas phases. This force-field was developed with the simulation of large biomolecular systems in mind.
23. Friesner R.A. Ab initio quantum chemistry: methodology and applications. Proc Natl Acad Sci USA 102 (2005) 6648-6653. This perspective article provides an overview of the state of the art in ab initio quantum chemical methodology and applications.
24. Van Der Spoel D., Lindahl E., Hess B., Groenhof G., Mark A.E., and Berendsen H.J. GROMACS: fast, flexible, and free. J Comput Chem 26 (2005) 1701-1718. This paper presents an overview of the GROMACS package for molecular simulations. GROMACS includes many tricks to improve speed. Its main strength, however, lies in its flexibility, and in the fact that it is in the public domain and distributed under the GNU General Public License.
25. Shen M.Y., and Freed K.F. A simple method for faster nonbonded force evaluations. J Comput Chem 26 (2005) 691-698. This paper introduces approximations for the evaluation of inverse distances, such that square root or division operations are not required.
26. Shaw D.E. A fast, scalable method for the parallel evaluation of distance-limited pairwise particle interactions. J Comput Chem 26 (2005) 1318-1328. This paper introduces a new method for the parallel evaluation of pairwise particle interactions. Its originality lies in the reduction of the amount of data transferred between processors, taking advantage of the fact that particles only interact with neighboring particles within a given cutoff.
27. Loncharich R.J., and Brooks B.R. The effects of truncating long-range forces on protein dynamics. Proteins 6 (1989) 32-45
28. Schreiber H., and Steinhauser O. Cutoff size does strongly influence molecular dynamics results on solvated polypeptides. Biochemistry 31 (1992) 5856-5860
29. Schreiber H., and Steinhauser O. Taming cut-off induced artifacts in molecular dynamics studies of solvated polypeptides. The reaction field method. J Mol Biol 228 (1992) 909-923
30. Schreiber H., and Steinhauser O. Molecular dynamics studies of solvated polypeptides: why cutoff scheme does not work. Chem Phys 168 (1992) 75-89
31. Beck D.A., Armen R.S., and Daggett V. Cutoff size need not strongly influence molecular dynamics results for solvated polypeptides. Biochemistry 44 (2005) 609-616. The authors revisit the problem of artifacts introduced by cutoff schemes in MD simulation. They show that, by using a proper force-shifted spherical cutoff, it is possible to correctly model peptide dynamics in solution.
32. Norberg J., and Nilsson L. On the truncation of long-range electrostatic interactions in DNA. Biophys J 79 (2000) 1537-1553
33. Garemyr R., and Elofsson A. Study of the electrostatics treatment in molecular dynamics simulations. Proteins 37 (1999) 417-428
34. Varekova R.D., Koca J., and Zhan C.G. Complexity and convergence of electrostatic and van der Waals energies within PME and cutoff methods. Int J Mol Sci 5 (2004) 154-173
35. Allen M.P., and Tildesley D.J. Computer Simulation of Liquids (1989), Oxford University Press, New York
36. Frenkel D., and Smit B. Understanding Molecular Simulation: From Algorithms To Applications (1996), Academic Press, San Diego
37. Darden T., York D., and Pedersen L. Particle mesh Ewald: an n log n method for Ewald sums in large systems. J Chem Phys 99 (1993) 10089-10092
38. Sagui C., and Darden T. Multigrid methods for classical molecular dynamics simulations of biomolecules. J Chem Phys 114 (2001) 6578-6591
39. Skeel R.D., Tezcan I., and Hardy D.J. Multiple grid methods for classical molecular dynamics. J Comput Chem 23 (2002) 673-684
40. Shan Y., Klepeis J.L., Eastwood M.P., Dror R.O., and Shaw D.E. Gaussian split Ewald: a fast Ewald mesh method for molecular simulation. J Chem Phys 122 (2005) 54101. The authors propose a new Ewald mesh method for computing the electrostatic interactions between particles in the case of explicit solvent with PBCs. This method involves two stages of charge spreading on the mesh; it can be used with both a real-space and a Fourier-space solver of the Poisson equation. Noteworthy in this paper is the introduction of a uniform framework that describes the differences between existing Ewald mesh methods.
41. Essmann U., Perera L., Berkowitz M.L., Darden T., Lee H., and Pedersen L.G. A smooth particle mesh Ewald method. J Chem Phys 103 (1995) 8577-8593
42. Crocker M.S., Hampton S.S., Matthey T., and Izaguirre J.A. MDSIMAID: automatic parameter optimization in fast electrostatic algorithms. J Comput Chem 26 (2005) 1021-1031. MDSimAID is a recommender system that optimizes parallel fast electrostatics solvers and presents recommended parameters to the user. It is part of a general effort (see [43]) to facilitate the use and development of new methods for computing electrostatic interactions in explicit solvent models.
43. Matthey T., Cickovski T., Hampton S., Ko A., Ma Q., Nyerges M., Raeder T., Slabach T., and Izaguirre J.A. ProtoMol, an object-oriented framework for prototyping novel algorithms for molecular dynamics. ACM Trans Math Softw 20 (2004) 237-265
44. Tironi I.G., Sperb R., Smith P.E., and van Gunsteren W.F. A generalized reaction field method for molecular dynamics simulations. J Chem Phys 102 (1995) 5451-5459
45. Heinz T.N., and Hunenberger P.H. Combining the lattice-sum and reaction-field approaches for evaluating long-range electrostatic interactions in molecular simulations. J Chem Phys 123 (2005) 34107. The authors describe a new scheme for computing electrostatic interactions in explicit solvent with PBCs that combines a mesh or lattice-sum approach and a reaction-field approach. They show that there is continuity between the lattice-sum and reaction-field schemes, and unify them into a common framework.
46. Wu X., and Brooks B.R. Isotropic periodic sum: a method for the calculation of long-range interactions. J Chem Phys 122 (2005) 44107. This paper describes a new approach to the calculation of long-range interactions in explicit solvent models. It defines a local region around each particle and describes the remaining region as images of the local region, distributed statistically in an isotropic and periodic way.
47. Petrella R.J., and Karplus M. Electrostatic energies and forces computed without explicit interparticle interactions: a linear time complexity formulation. J Comput Chem 26 (2005) 755-787. The authors describe a new, rapid method for the calculation of electrostatic interactions between particles. It does not require a cutoff distance and the computing time scales linearly with the number of particles or charges. The main idea presented is to replace the double sum in the Coulomb summation with the product of two simple sums, using tricks introduced to compute Euler sums. It does require the inverse of the distance to be approximated with a polynomial.
48. Bailey D.H., Borwein J.M., and Girgensohn R. Experimental evaluation of Euler sums. Exp Math 3 (1994) 17-28
49. Hunenberger P.H., and McCammon J.A. Effect of artificial periodicity in simulations of biomolecules under Ewald boundary conditions: a continuum electrostatics study. Biophys Chem 78 (1999) 69-88
50. Hunenberger P.H., and McCammon J.A. Ewald artifacts in computer simulations of ionic solvation and ion-ion interaction: a continuum electrostatics study. J Chem Phys 110 (1999) 1856-1872
51. Weber W., Hunenberger P.H., and McCammon J.A. Molecular dynamics simulations of a polyalanine octapeptide under Ewald boundary conditions: influence of artificial periodicity on peptide conformation. J Phys Chem B104 (2000) 3668-3675
52. Villarreal M.A., and Montich G.G. On the Ewald artifacts in computer simulations. The test-case of the octaalanine peptide with charged termini. J Biomol Struct Dyn 23 (2005) 135-142. The authors re-examine the problem of artifacts detected in molecular simulations based on Ewald summation. Using simulations of octa-alanine peptides, they show that the artifacts are more a consequence of poor sampling than of the Ewald summation technique.
53. Monticelli L., and Colombo G. The influence of simulation conditions in molecular dynamics investigations of model beta-sheet peptides. Theoretical Chemical Accounts 3 (2004) 145-157
54. Halgren T.A., and Damm W. Polarizable force fields. Curr Opin Struct Biol 11 (2001) 236-242
55. Jorgensen W.L., and Tirado-Rives J. Potential energy functions for atomic-level simulations of water and organic and biomolecular systems. Proc Natl Acad Sci USA 102 (2005) 6665-6670. An excellent review of potential energy functions for molecular simulations.
56. Caldwell J.W., and Kollman P.A. Structure and properties of neat liquids using non additive molecular mechanics: water, methanol, and N-methylacetamide. J Phys Chem 99 (1995) 6208-6219
57. Cubero E., Luque F.J., and Orozco M. Is polarization important in cation-pi interactions?. Proc Natl Acad Sci USA 95 (1998) 5976-5980
58. Rick S.W., and Stuart S.J. Potentials and algorithms for incorporating polarizability in computer simulations. Rev Comp Chem 18 (2002) 89-146
59. Kaminski G.A., Stern H.A., Berne B.J., Friesner R.A., Cao Y.X., Murphy R.B., Zhou R., and Halgren T.A. Development of a polarizable force field for proteins via ab initio quantum chemistry: first generation model and gas phase tests. J Comput Chem 23 (2002) 1515-1531
60. Patel S., Mackerell Jr. A.D., and Brooks III C.L. CHARMM fluctuating charge force field for proteins: II protein/solvent properties from molecular dynamics simulations using a nonadditive electrostatic model. J Comput Chem 25 (2004) 1504-1514
61. Baucom J., Transue T., Fuentes-Cabrera M., Krahn J.M., Darden T.A., and Sagui C. Molecular dynamics simulations of the d(CCAACGTTGG)2 decamer in crystal environment: comparison of atomic point-charge, extra-point, and polarizable force fields. J Chem Phys 121 (2004) 6998-7008
62. Anisimov V.M., Lamoureux G., Vorobyov I.G., Huang N., Roux B., and Mackerell Jr. A.D. Determination of electrostatic parameters for a polarizable force field based on the classical Drude oscillator. J Chem Theory Comput 1 (2005) 153-168. The authors have developed a polarizable empirical force-field based on the classical Drude oscillator model. They incorporated it into CHARMM and have tested it in a 1000 ps simulation of a DNA octamer in a box of water with sodium counterions, again proving the feasibility of using polarizable force-fields (see also [59-61]).
63. Wisz M.S., and Hellinga H.W. An empirical model for the electrostatic interactions in proteins incorporating multiple geometry-dependent dielectric constants. Proteins 51 (2003) 360-377
64. Schwarzl S.M., Huang D., Smith J.C., and Fisher S. Nonuniform charge scaling (NUCS): a practical approximation of solvent electrostatic screening in proteins. J Comput Chem 26 (2005) 1359-1371. The authors have developed a pragmatic approach to implicit solvation that approximates solvent screening by individually scaling the partial charges of the atoms, so as to reproduce electrostatic interactions obtained from a PB analysis. The scaled charges are consequently maintained as fixed during molecular simulations.
65. Feig M., Onufriev A., Lee M.S., Im W., Case D.A., and Brooks III C.L. Performance comparison of generalized Born and Poisson methods in the calculation of electrostatic solvation energies for protein structures. J Comput Chem 25 (2004) 265-284
66. David L., Luo R., and Gilson M. Comparison of generalized Born and Poisson models: energetics and dynamics of HIV protease. J Comput Chem 21 (2000) 295-309
67. Lu B.Z., Chen W.Z., Wang C.X., and Xu X.J. Protein molecular dynamics with electrostatic force entirely determined by a single Poisson-Boltzmann calculation. Proteins 48 (2002) 497-504
68. Luo R., David L., and Gilson M.K. Accelerated Poisson-Boltzmann calculations for static and dynamic systems. J Comput Chem 23 (2002) 1244-1253
69. Fogolari F., Brigo A., and Molinari H. Protocol for MM/PBSA molecular dynamics simulations of proteins. Biophys J 85 (2003) 159-166
70. Prabhu N.V., Zhu P., and Sharp K.A. Implementation and testing of stable, fast implicit solvation in molecular dynamics using the smooth-permittivity finite difference Poisson-Boltzmann method. J Comput Chem 25 (2004) 2049-2064
71. Marshall S.A., Vizcarra C.L., and Mayo S.L. One and two-body decomposable Poisson-Boltzmann methods for protein design calculations. Protein Sci 14 (2005) 1293-1304. One-body and two-body decomposable PB methods are proposed for computing electrostatic interactions. These approximate methods were specifically developed for protein design.
72. Lu B., Zhang D., and McCammon J.A. Computation of electrostatic forces between solvated molecules determined by the Poisson-Boltzmann equation using a boundary element method. J Chem Phys 122 (2005) 214102. The authors propose a new approach for calculating electrostatics forces among solvated molecules when the electric field is obtained by solving the linearized PB equation. They use the variational principle in the context of BEM. This approach does not require the calculation of the stress tensor on the surface of the molecule and therefore avoids the problem of hypersingularity.
73. ••] are still at the developmental stage.
74. Hofinger S. Solving the Poisson-Boltzmann equation with the specialized computer chip MD-GRAPE-2. J Comput Chem 26 (2005) 1148-1154. A new way of solving the PB equation on a special-purpose hardware chip is presented. The author reports acceleration factors of up to 40-fold with no loss of accuracy compared with conventional approaches.
75. Still W.C., Tempczyk A., Hawley R.C., and Hendrickson T. Semi-analytical treatment of solvation for molecular mechanics and dynamics. J Am Chem Soc 112 (1990) 6127-6129
76. Zhu J., Alexov E., and Honig B. Comparative study of generalized Born models: Born radii and peptide folding. J Phys Chem B 109 (2005) 3008-3022. Four GB models are compared in terms of their ability to reproduce the results of PB calculations and their performance in the ab initio folding of two peptides. Surprisingly, the method that is most successful in reproducing PB results is not the most effective in the folding experiments.
77. ••] are compared in terms of their performance in retaining native properties of protein structures during MD simulations. The highlight of this paper is the development of a statistical framework to verify that the results of MD simulations are meaningful.
78. as and related characteristics of biomolecules. This server is designed to facilitate access to electrostatics analysis by a large community, including non-experts.
79. •] are two independent servers providing similar services.
80. Bashford D. Macroscopic electrostatic models for protonation states in proteins. Front Biosci 9 (2004) 1082-1099
81. Dolinsky T.J., Nielsen J.E., McCammon J.A., and Baker N.A. PDB2PQR: an automated pipeline for the setup of Poisson-Boltzmann electrostatics calculations. Nucleic Acids Res 32 (2004) W665-W667
82. Gilson M.K., Sharp K., and Honig B. Calculating the electrostatic potential of molecules in solution: method and error assessment. Proteins 4 (1988) 7-18
83. Baker N.A., Sept D., Joseph S., Holst M.J., and McCammon J.A. Electrostatics of nanosystems: application to microtubules and the ribosome. Proc Natl Acad Sci USA 98 (2001) 10037-10041
84. Grant J.A., Pickup B.T., and Nicholls A. A smooth permittivity function for Poisson Boltzmann solvation methods. J Comput Chem 22 (2001) 608-640
85. Sayyed-Ahmad A., Tuncay K., and Ortoleva P. Efficient solution technique for solving the Poisson-Boltzmann equations. J Comput Chem 25 (2004) 1068-1074
86. Nicholls A., Sharp K., and Honig B. Protein folding and association: insights from the interfacial and thermodynamic properties of hydrocarbons. Proteins 11 (1991) 281-296
87. Petrey D., and Honig B. GRASP2: visualization, surface properties, and electrostatics of macromolecular structures and sequences. Methods Enzymol 374 (2003) 492-509