Multidimensional adaptive umbrella sampling: Applications to main chain and side chain peptide conformations

A new adaptive umbrella sampling technique for molecular dynamics simulations is described. The high efficiency of the technique renders multidimensional adaptive umbrella sampling possible and thereby enables uniform sampling of the conformational space spanned by several degrees of freedom. The efficiency is achieved by using the weighted histogram analysis method to combine the results from different simulations, by a suitable extrapolation scheme to define the umbrella potential for regions that have not been sampled, and by a criterion to identify simulations during which the system was not in equilibrium. The technique is applied to two test systems, the alanine dipeptide and the threonine dipeptide, to sample the configurational space spanned by one or two dihedral angles. The umbrella potentials applied at the end of each adaptive umbrella sampling run are equal to the negative of the corresponding potentials of mean force. The trajectories obtained in the simulations can be used to calculate dynamical variables that are of interest. An example is the distribution of the distance between the HN and the Hβ proton that can be important for the interpretation of NMR experiments. Factors influencing the accuracy of the calculated quantities are discussed. © 1997 John Wiley & Sons, Inc. J Comput Chem 18: 1450–1462, 1997

[1]  M. Mezei.  Adaptive umbrella sampling: Self-consistent determination of the non-Boltzmann bias , 1987 .

[2]  Yuko Okamoto,et al.  Comparative Study of Multicanonical and Simulated Annealing Algorithms in the Protein Folding Problem , 1994 .

[3]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[4]  C. Brooks,et al.  Constant-temperature free energy surfaces for physical and chemical processes , 1993 .

[5]  Shankar Kumar,et al.  Multidimensional free‐energy calculations using the weighted histogram analysis method , 1995, J. Comput. Chem..

[6]  J. Mccammon,et al.  Dynamics of Proteins and Nucleic Acids , 2018 .

[7]  Charles L. Brooks,et al.  Reaction paths and free energy profiles for conformational transitions: An internal coordinate approach , 1991 .

[8]  J Hermans,et al.  Microfolding: Conformational probability map for the alanine dipeptide in water from molecular dynamics simulations , 1988, Proteins.

[9]  M Karplus,et al.  Dynamical theory of activated processes in globular proteins. , 1982, Proceedings of the National Academy of Sciences of the United States of America.

[10]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

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

[12]  H. Scheraga,et al.  Prediction of the native conformation of a polypeptide by a statistical‐mechanical procedure. I. Backbone structure of enkephalin , 1985, Biopolymers.

[13]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[14]  Charles L. Brooks,et al.  Conformational equilibrium in the alanine dipeptide in the gas phase and aqueous solution : a comparison of theoretical results , 1992 .

[15]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[16]  W. F. Gunsteren,et al.  Motion and Conformation of Side Chains in Peptides. A Comparison of 2D Umbrella-Sampling Molecular Dynamics and NMR Results , 1996 .

[17]  Yuko Okamoto,et al.  Prediction of peptide conformation by multicanonical algorithm: New approach to the multiple‐minima problem , 1993, J. Comput. Chem..

[18]  Umbrella sampling along linear combinations of generalized coordinates. Theory and application to a glycine dipeptide , 1995 .

[19]  A. Lyubartsev,et al.  New approach to Monte Carlo calculation of the free energy: Method of expanded ensembles , 1992 .

[20]  Berg,et al.  Multicanonical ensemble: A new approach to simulate first-order phase transitions. , 1992, Physical review letters.

[21]  Lee,et al.  New Monte Carlo algorithm: Entropic sampling. , 1993, Physical review letters.

[22]  Roland L. Dunbrack,et al.  Backbone-dependent rotamer library for proteins. Application to side-chain prediction. , 1993, Journal of molecular biology.

[23]  R. Hooft,et al.  An adaptive umbrella sampling procedure in conformational analysis using molecular dynamics and its application to glycol , 1992 .

[24]  A. Kidera,et al.  Enhanced conformational sampling in Monte Carlo simulations of proteins: application to a constrained peptide. , 1995, Proceedings of the National Academy of Sciences of the United States of America.

[25]  I N Bronstein,et al.  Taschenbuch der Mathematik , 1966 .

[26]  Harold A. Scheraga,et al.  MONTE CARLO SIMULATION OF A FIRST-ORDER TRANSITION FOR PROTEIN FOLDING , 1994 .

[27]  U. Krengel.  Einführung in die Wahrscheinlichkeitstheorie und Statistik , 1988 .

[28]  Shankar Kumar,et al.  Method for free‐energy calculations using iterative techniques , 1996 .

[29]  M. Karplus,et al.  Proteins: A Theoretical Perspective of Dynamics, Structure, and Thermodynamics , 1988 .