Premium
Multidimensional adaptive umbrella sampling: Applications to main chain and side chain peptide conformations
Author(s) -
Bartels Christian,
Karplus Martin
Publication year - 1997
Publication title -
journal of computational chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.907
H-Index - 188
eISSN - 1096-987X
pISSN - 0192-8651
DOI - 10.1002/(sici)1096-987x(199709)18:12<1450::aid-jcc3>3.0.co;2-i
Subject(s) - umbrella sampling , extrapolation , molecular dynamics , sampling (signal processing) , dihedral angle , degrees of freedom (physics and chemistry) , dipeptide , potential of mean force , statistical physics , chemistry , algorithm , computational chemistry , computer science , mathematics , physics , thermodynamics , mathematical analysis , molecule , organic chemistry , hydrogen bond , biochemistry , filter (signal processing) , amino acid , computer vision
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