Novel methods of sampling phase space in the simulation of biological systems

Current Opinion in Structural Biology

Volumn 7, Issue 2, 1997, Pages 181-189

  Berne, Bruce J ^a   Straub, John E ^b  

^a Och Spine at New York Presbyterian Hospitals   (United States)

^b Boston University   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

ALGORITHM; CHEMICAL STRUCTURE; COMPARATIVE STUDY; CONFORMATION; ELECTRICITY; MOLECULAR BIOLOGY; MOLECULAR DYNAMICS; PRIORITY JOURNAL; REVIEW; SAMPLING; SIMULATION; THERMODYNAMICS;

EID: 0030904524     PISSN: 0959440X     EISSN: None     Source Type: Journal    
DOI: 10.1016/S0959-440X(97)80023-1     Document Type: Article

References (76)

1. Thirumalai D, Mountain RD, Kirkpatrick TR: Ergodic behavior in supercooled liquids and glasses. Phys Rev A 1989, 39:3563-3574.
2. Straub JE, Rashkin A, Thirumalai D: Dynamics in rugged energy landscapes with applications to the S-peptide and ribonuclease A. J Am Chem Soc 1994, 116:2049-2063.
3. Hodel A, Simonsen T, Fox RO, Brünger AT: Conformational substates and uncertainty in macromolecular free energy calculations. J Phys Chem 1993, 97:3409-3417.
4. Allen MP, Tildesley, DJ: Computer Simulation of Liquids. Oxford: Oxford University Press; 1991.
5. Duane S, Kennedy AD, Pendleton BJ, Roweth D: Hybrid Monte Carlo. Phys Lett B 1987, 195:216-221.
6. Bennett CH: Mass tensor molecular dynamics. J Comp Phys 1975, 19:267-279.
7. Berne BJ: Molecular dynamics and Monte Carlo simulations of rare events. In Multiple Time Scales. Edited by Brackbill JU, Cohen BI. New York: Academic Press; 1985:419-436.
8. Tuckerman ME, Martyna GJ, Berne BJ: Reversible multiple time scale molecular dynamics. J Chem Phys 1992, 97:1990-2001.
9. Stuart SJ, Zhou R, Berne BJ: Molecular dynamics with multiple timescales: The selection of efficient reference system propagators. J Chem Phys 1996, 105:1426-1436. Various strategies for factorizing the propagator in r-RESPA are discussed. This paper shows which strategies to avoid and also how to subdivide the Fourier space contribution to the force in Ewald simulations into fast and slow components and how to incorporate this subdivision into an efficient implementation of r-RESPA.
10. Humphreys D, Friesner RA, Berne BJ: A multiple time scale molecular dynamics algorithm for macromolecules. J Phys Chem 1994, 98:6885-6892.
11. Humphreys D, Friesner RA, Berne BJ: Simulated annealing of a protein in a continuum solvent by multiple-time-step molecular dynamics. J Phys Chem 1995, 99:10674-10685. Building on work in [10], r-RESPA is applied to proteins using the generalized Born solvation model.
12. Zhou R, Berne BJ: A new molecular dynamics method combining the reference system propagator algorithm with a fast multipole method for simulating proteins and other complex systems. J Chem Phys 1995, 103:9444-9458. This paper combines a new top-down FMM with r-RESPA for simulating proteins. The combination of these two powerful methodologies leads to better than an order of magnitude speed-up for large proteins such as the photosynthetic reaction center.
13. Watanabe M, Karplus M: Dynamics of molecules with internal degrees of freedom by multiple time-step methods. J Chem Phys 1993, 99:8063-8074.
14. Procacci P, Marchi M: Taming the Ewald sum in molecular dynamics simulations of solvated proteins via a multiple time step algorithm. J Chem Phys 1996, 104:3003-3012. In this paper, standard Ewald is combined with r-RESPA. The whole of the Fourier space contribution to the force is included in an outer loop of r-RESPA. Because this sum contains a rapidly varying short time contribution, this is not the best strategy for making the r-RESPA break-up. Despite this, this paper reports a significant speed-up over non-r-RESPA implementations.
15. Greengard L: The Rapid Evaluation of Potential Fields in Particle Systems. Cambridge, MA: MIT Press; 1988.
16. Greengard L, Rokhlin V: On the evaluation of electrostatic interactions in molecular modeling. Phys Scr 1989, 29:139-144.
17. Board JA, Causey JW, Leathrum JF, Windemuth A, Schulten K: Accelerated molecular dynamics simulation with the parallel fast multipole algorithm. Chem Phys Lett 1992, 198:89-94.
18. Shimada J, Kaneko H, Takada T: Performance of Fast Multipole Methods for calculating electrostatic interactions in biomacromolecular simulations. J Comp Chem 1994, 15:28-43.
19. White CA, Head-Gordon M: Derivation and efficient implementation of the fast multipole method. J Chem Phys 1994, 101:6593-6605.
20. Ding H-Q, Karasawa N, Goddard WA III: Atomic level simulations on a million particles: the cell multipole method for Coulomb and London nonbonded interactions. J Chem Phys 1992, 97:4309-4315.
21. Lee FS, Warshel A: A local reaction field method for fast evaluation of long-range electrostatic interactions in molecular simulations. J Chem Phys 1992, 97:3100-3107.
22. Stote RH, States DJ, Karplus M: On the treatment of electrostatic interactions in biomolecular simulation. J Chem Phys 1991, 88:2419-2433.
23. Saito M: Molecular dynamics simulations of proteins in water without the truncation of long-range Coulomb interactions. Mol Simulat 1992, 8:321-333.
24. Shimada J, Kaneko H, Takada T: Efficient calculations of Coulombic interactions in biomolecular simulations with periodic boundary conditions. J Comp Chem 1993, 14:867-878.
25. Mathiowetz AM, Jain A, Karasawa N, Goddard WA III: Protein simulations using techniques suitable for very large systems: the Cell Multipole Method for nonbonded interactions and the Newton-Euler Inverse Mass Operator Method for internal coordinate dynamics. Proteins 1994, 20:227-247.
26. Lee MA, Schmidt KE: Implementing the Fast Multipole Method in three dimensions. J Stat Phys 1991, 63:1223-1235.
27. Darden TA, York DM, Pedersen LG: Particle Mesh Ewald: an N log(N) method for Ewald sums in large systems. J Chem Phys 1993, 98:10089-10092.
28. Petersen HG: Accuracy and efficiency of the Particle Mesh Ewald method. J Chem Phys 1995, 103:3668-3679.
29. Essman U, Perera L, Berkowitz ML, Darden T, Lee H, Pedersen LG: A smooth Particle Mesh Ewald method. J Chem Phys 1995, 103:8577-8593.
30. Bernardo DN, Ding Y, Krogh-Jespersen K, Levy RM: An anisotropic polarizable water model: incorporation of all-atom polarizabilities into molecular mechanics force fields. J Phys Chem 1994, 98:4180-4187.
31. Stuart SJ, Berne BJ: Effects of polarizability on the hydration of the chloride ion. J Phys Chem 1996, 100:11934-11943. This paper shows how to incorporate ionic polarizability into a simulation of ionic solubility in polarizable water using the fluctuating charge model and a Drude oscillator polarizability model for the ion. r-RESPA is used to deal with the fast vibrations of the Drude oscillator.
32. Cao J, Berne BJ: Theory and simulation of polar and non-polar polarizable fluids. J Chem Phys 1993, 99:6998-7011.
33. Rick SW, Stuart S, Berne BJ: Dynamical fluctuating charge force fields: application to liquid water. J Chem Phys 1994, 101:6141-6161.
34. Rick SW, Berne BJ: The aqueous solvation of water: a comparison of molecular and continuum methods. J Am Chem Soc 1994, 116:3949-3954.
35. Rick S, Stuart S, Bader JS, Berne BJ: Fluctuating charge force fields for aqueous solutions. J Mol Liq 1995, 65/66:31-40.
36. Rick S, Berne BJ: Dynamical fluctuating charge force fields: the aqueous solvation of amides. J Am Chem Soc 1996, 118:672-679.
37. Northup SH, McCammon JA: Simulation methods for protein structure fluctuations. Biopolymers 1980, 19:1001-1016.
38. Hannsmann UHE, Okamoto Y, Eisenmenger F: Molecular dynamics, Langevin and hybrid Monte Carlo simulations in a mutticanonical ensemble. Chem Phys Lett 1996, 259:321-330. The paper presents a number of simulation methods based on sampling the multicanonical probability distribution using a variety of algorithms. There is no critical discussion of the relative efficiency of the various algorithms.
39. Gromov DG, De Pablo JJ: Structure of binary polymer blends: multiple time step hybrid Monte Carlo simulations and self-consistent integral equation theory. J Chem Phys 1995, 103:8247-8256. HMC was first invented in [5]. The need for such a method has an interesting history. Parrinello and Rahman first used MD to simulate quantum degrees of freedom. Calloway and Rahman then showed how MD can be used for lattice-gauge simulations. Hall and Berne pointed out that MD could not often be used by itself to simulate path integrals because the dynamics would often be in the KAM regime where they would not be ergodic. Berne and Hall suggested that one should then punctuate the MD evolution by sampling the momenta from the Maxwell distribution. They then showed that the MD simulation would agree with MC simulations. This gives rise to BGK dynamics, a type of stochastic dynamics used by Montgomery, Chandler and Berne to simulate activated barrier crossing and later by Andersen in his canonical MD method. In BGK dynamics, exact Hamiltonian flow between collisions gives the canonical distribution, as first shown by Lebowitz and Bergmann. Unfortunately, when the equations of motion are integrated numerically, one cannot prove that the states sample the canonical distribution. Duane et al. [5] then introduced the Metropolis criterion for accepting or rejecting the momentum refreshed MD step. Their HMC does rigorously give the correct canonical distribution and is thus the method of choice.
40. Forrest BM, Suter UW: Hybrid Monte Carlo simulations of dense polymer system. J Chem Phys 1994, 101:2616-2029.
41. Irbäck A: Hybrid Monte Carlo simulation of polymer chains. J Chem Phys 1994, 101:1661-1667.
42. Neal RM: An improved acceptance procedure for the Hybrid Monte Carlo algorithm. J Comp Phys 1994, 111:194-203.
43. Gupta S, Irbäck A, Karsch F, Petersson B: The acceptance probability in the Hybrid Monte Carlo method. Phys Lett B 1990, 242:437-443.
44. Mackenzie PB: An improved Hybrid Monte Carlo method. Phys Lett B 1989, 226:369-371.
45. Nosé S: A unified formulation of the constant temperature molecular dynamic methods. J Chem Phys 1984, 81:511-519.
46. Tuckerman M, Berne BJ, Martyna G, Klein M: Efficient molecular dynamics and hybrid Monte Carlo algorithm for path integrals. J Chem Phys 1993, 99:2796-2808.
47. Frantz DD, Freeman DL, Doll JD: Reducing quasi-ergodic behavior in Monte Carlo simulation by J-walking: applications to atomic clusters. J Chem Phys 1990, 93:2769-2784.
48. Frantz DD, Freeman DL, Doll JD: Extending J-walking to quantum systems: applications to atomic clusters. J Chem Phys 1992, 97:5713-5731.
49. Matro A, Freeman DL, Topper RQ: Computational study of the structures and thermodynamic properties of ammonium chloride clusters using a parallel jump-walking approach. J Chem Phys 1996, 104:8690-8702.
50. Freeman DL, Doll JD: Computational studies of clusters. Methods and results. Annu Rev Phys Chem 1996, 47:43-80.
51. Strozak MA, Lopez GE, Freeman DL: Gibbs free-energy changes for the growth of argon clusters adsorbed on graphite. J Chem Phys 1992, 97:4445-4452.
52. Lopez GE, Freeman DL: A study of low temperature heat capacity anomalies in bimetallic alloy clusters using J-walking Monte Carlo methods. J Chem Phys 1993, 98:1428-1435.
53. Liu Z, Berne BJ: Methods for accelerating chain folding and mixing. J Chem Phys 1993, 99:6071-6077.
54. Berg BA, Neuhaus T: Multicanonical algorithms for first order phase transitions. Phys Lett B 1991, B267:249-253.
55. Hansmann UHE, Okamoto Y: Comparative study of multicanonical and simulated annealing algorithms in the protein folding problem. Physica A 1994, 212:415-437.
56. Okamoto Y, Hansmann UHE: Thermodynamics of helix-coil transitions studied by multicanonical algorithms. J Phys Chem 1995, 99:11276-11287 An application of multicanonical sampling to examine helix-coil transitions in isolated polypeptides ranging from 10-20 amino acids.
57. Hao M-H, Scheraga HA: Monte Carlo simulation of a first-order transition for protein folding. J Phys Chem 1994, 98:4940-4948.
58. Tsallis C: Possible generalization of Boltzmann-Gibbs statistics. J Stat Phys 1988, 52:479-487.
59. Andricioaei I, Straub JE: Generalized simulated annealing algorithms using Tsallis statistics: application to conformational optimization of a tetrapeptide. Phys Rev E 1996, 53:R3055-R3058. A generalized MC algorithm that samples equilibrium Tsallis statistical distributions is presented. The relative efficiency of a simulated annealing protocol employing the generalized Monte Carlo algorithm is shown to be significantly more effective than MD or HMC based annealing algorithms in isolating the global energy minimum of a tetrapeptide.
60. Guarnieri F, Still WC: A rapidly convergent simulation method: Mixed Monte Carlo/stochastic dynamics. J Comp Chem 1994, 11:1302-1310.
61. Senderowitz H, Guarnieri F, Still, WC: A smart Monte Carlo technique for free energy simulations of multiconformational molecules. Direct calculations of the conformational populations of organic molecules. J Am Chem Soc 1995, 117:8211-8219.
62. Stillinger FH, Stillinger DK: Cluster optimization simplified by interaction modification. J Chem Phys 1990, 93:6106-6107.
63. Head-Gordon T, Stillinger FH: Predicting polypeptide and protein structures from amino acid sequence: Antlion method applied to melittin. Biopolymers 1993, 33:293-303.
64. Piela L, Kostrowicki J, Scheraga HA: The multiple-minima problem in the conformational analysis of molecules. Deformation of the potential energy hypersurface by the Diffusion Equation Method. J Phys Chem 1989, 93:3339-3346.
65. Kostrowicki J, Scheraga HA: Application of the Diffusion Equation Method for global optimization to oligopeptides. J Phys Chem 1992, 96:7442-7449.
66. Kostrowicki J, Scheraga HA: Some approaches to the multiple-minima problem in protein folding. In Global Minimization of Nonconvex Energy Functions: Molecular Conformation and Protein Folding. Edited by Pardalos PM, Shalloway D, Xue G. Providence: American Mathematical Society; 1996:123-132. A discussion of methods for extending the DEM to internal coordinates (dihedral-angle space) is presented.
67. Shalloway D: Recent advances in global optimization. Edited by Floudas CA, Pardalos PM. Princeton: Princeton University Press; 1992:433-477.
68. Church BW, Orešič M, Shalloway D: Tracking metastable states to free-energy global minima. In Global Minimization of Nonconvex Energy Functions: Molecular Conformation and Protein Folding. Edited by Pardalos PM, Shalloway D, Xue G. American Mathematical Society; Providence; 1996:41-64. This paper provides and excellent introduction to the PA method. The PA method is related to a reduced Smoluchowski diffusion equation making an intriguing connection to diffusive kinetics between conformational substates. Application of the PA method to Met-enkephalin is described as a means of gaining insight into the pathways available to peptides and proteins during annealing or folding.
69. Amara P, Hsu D, Straub JE: Global energy minimum searches using an approximate solution of the imaginary time Schrödinger equation. J Phys Chem 1993, 97:6715-6721.
70. Ma J, Hsu D, Straub JE: Approximate solution of the classical Liouville equation using Gaussian phase packet dynamics: application to enhanced equilibrium averaging and global optimization. J Chem Phys 1993, 49:4024-4035.
71. Straub JE: Optimization techniques with applications to proteins. In Recent Developments in Theoretical Studies of Proteins. Edited by Elber R. Singapore: World Scientific; 1996:137-196. This review presents a complete overview of potential-smoothing algorithms applied to global energy minimization of biomolecular systems.
72. Amara P, Straub JE: Folding model proteins using kinetic and thermodynamic annealing of the classical density distribution. J Phys Chem 1995, 99:14840-14852. This paper presents a study of global energy minimization for a series of model proteins using simulated annealing and potential-smoothing algorithms. A critical evaluation of Gaussian potential-smoothing methods and a comparison with MD simulated annealing is presented.
73. Roitberg A, Elber R: Modeling side chains in peptides and proteins: application of the locally enhanced sampling and the simulated annealing methods to find minimum energy conformations. J Chem Phys 1991, 95:9277-9287.
74. Verkhivker G, Elber R, Nowak W: Locally enhanced sampling in free energy calculations: application of mean field approximation to accurate calculation of free energy differences. J Chem Phys 1992, 97:7838-7841.
75. Finnila AB, Gomez MA, Sebenik C, Stenson C, Doll JD: Quantum annealing: a new method for minimizing multidimensional functions. Chem Phys Lett 1994, 219:343-348.
76. Andricioaei I, Straub JE: Finding the needle in the haystack: algorithms for conformational optimization. Comput Phys 1996, 10:449-454. This tutorial introduction to molecular conformational optimization presents an overview of simulated annealing, Gaussian potential smoothing, quantum mechanical annealing and the Method of Bad Derivatives.