Monte Carlo vs molecular dynamics for conformational sampling

被引:125
作者
Jorgensen, WL
TiradoRives, J
机构
[1] Department of Chemistry, Yale University, New Haven
关键词
D O I
10.1021/jp960880x
中图分类号
O64 [物理化学(理论化学)、化学物理学];
学科分类号
070304 ; 081704 ;
摘要
A comparison study has been carried out to test the relative efficiency of Metropolis Monte Carlo and molecular dynamics simulations for conformational sampling. The test case that has been examined is liquid hexane. OPLS potential functions have been used with sampling of external motions and internal angle bending and torsions. The BOSS 3.6 and AMBER 4.0 programs were used to perform the internal-coordinate Monte Carlo and Cartesian molecular dynamics simulations, respectively, for samples of 267 molecules in the NPT ensemble at 25 degrees C and 1 atm with periodic boundary conditions. The initial configurations of the system were equilibrated with the-hexane molecules in the all trans configuration. The major findings are (1) the Monte Carlo and molecular dynamics results for thermodynamic properties and conformer populations are in accord, (2) the conformer populations reach equilibrium in ca. 10(7) Monte Carlo configurations or 100 ps of molecular dynamics, and (3) the molecular dynamics calculations require 1.6-3.8 times more computer time to achieve the same level of convergence as the Monte Carlo simulations. The viability of internal-coordinate Monte Carlo methodology contrasts with earlier notions on the inefficiency of Cartesian Monte Carlo methods. Monte Carlo sampling in related contexts such as for side-chain torsional motion in proteins appears promising.
引用
收藏
页码:14508 / 14513
页数:6
相关论文
共 19 条
[1]  
Allen M. P., 1987, J COMPUTER SIMULATIO, DOI DOI 10.2307/2938686
[2]   MOLECULAR-DYNAMICS WITH COUPLING TO AN EXTERNAL BATH [J].
BERENDSEN, HJC ;
POSTMA, JPM ;
VANGUNSTEREN, WF ;
DINOLA, A ;
HAAK, JR .
JOURNAL OF CHEMICAL PHYSICS, 1984, 81 (08) :3684-3690
[3]  
BROOKS CL, 1988, ADV CHEM PHYS, V71, P1
[4]  
Deem MW, 1996, MOL PHYS, V87, P1245, DOI 10.1080/00268979650026785
[5]   A CONCERTED ROTATION ALGORITHM FOR ATOMISTIC MONTE-CARLO SIMULATION OF POLYMER MELTS AND GLASSES [J].
DODD, LR ;
BOONE, TD ;
THEODOROU, DN .
MOLECULAR PHYSICS, 1993, 78 (04) :961-996
[6]   OPTIMIZED INTERMOLECULAR POTENTIAL FUNCTIONS FOR LIQUID HYDROCARBONS [J].
JORGENSEN, WL ;
MADURA, JD ;
SWENSON, CJ .
JOURNAL OF THE AMERICAN CHEMICAL SOCIETY, 1984, 106 (22) :6638-6646
[7]   THEORETICAL-STUDIES OF MEDIUM EFFECTS ON CONFORMATIONAL EQUILIBRIA [J].
JORGENSEN, WL .
JOURNAL OF PHYSICAL CHEMISTRY, 1983, 87 (26) :5304-5314
[8]  
JORGENSEN WL, 1995, BOSS VERSION 3 6
[9]   FREE-ENERGY CALCULATIONS - APPLICATIONS TO CHEMICAL AND BIOCHEMICAL PHENOMENA [J].
KOLLMAN, P .
CHEMICAL REVIEWS, 1993, 93 (07) :2395-2417
[10]   A COMPREHENSIVE STUDY OF THE ROTATIONAL ENERGY PROFILES OF ORGANIC-SYSTEMS BY AB-INITIO MO THEORY, FORMING A BASIS FOR PEPTIDE TORSIONAL PARAMETERS [J].
MAXWELL, DS ;
TIRADORIVES, J ;
JORGENSEN, WL .
JOURNAL OF COMPUTATIONAL CHEMISTRY, 1995, 16 (08) :984-1010