TY - JOUR
T1 - Enhanced sampling by multiple molecular dynamics trajectories
T2 - Carbonmonoxy myoglobin 10 μs A0 → A1-3 transition from ten 400 picosecond simulations
AU - Loccisano, Anne E.
AU - Acevedo, Orlando
AU - DeChancie, Jason
AU - Schulze, Brita G.
AU - Evanseck, Jeffrey D.
PY - 2004/5/1
Y1 - 2004/5/1
N2 - The utility of multiple trajectories to extend the time scale of molecular dynamics simulations is reported for the spectroscopic A-states of carbonmonoxy myoglobin (MbCO). Experimentally, the A0→A1-3 transition has been observed to be 10μs at 300K, which is beyond the time scale of standard molecular dynamics simulations. To simulate this transition, 10 short (400ps) and two longer time (1.2ns) molecular dynamics trajectories, starting from five different crystallographic and solution phase structures with random initial velocities centered in a 37Å radius sphere of water, have been used to sample the native-fold of MbCO. Analysis of the ensemble of structures gathered over the cumulative 5.6ns reveals two biomolecular motions involving the side chains of His64 and Arg45 to explain the spectroscopic states of MbCO. The 10μs A0→A1-3 transition involves the motion of His64, where distance between His64 and CO is found to vary up to 8.8±1.0Å during the transition of His64 from the ligand (A 1-3) to bulk solvent (A0). The His64 motion occurs within a single trajectory only once, however the multiple trajectories populate the spectroscopic A-states fully. Consequently, multiple independent molecular dynamics simulations have been found to extend biomolecular motion from 5ns of total simulation to experimental phenomena on the microsecond time scale.
AB - The utility of multiple trajectories to extend the time scale of molecular dynamics simulations is reported for the spectroscopic A-states of carbonmonoxy myoglobin (MbCO). Experimentally, the A0→A1-3 transition has been observed to be 10μs at 300K, which is beyond the time scale of standard molecular dynamics simulations. To simulate this transition, 10 short (400ps) and two longer time (1.2ns) molecular dynamics trajectories, starting from five different crystallographic and solution phase structures with random initial velocities centered in a 37Å radius sphere of water, have been used to sample the native-fold of MbCO. Analysis of the ensemble of structures gathered over the cumulative 5.6ns reveals two biomolecular motions involving the side chains of His64 and Arg45 to explain the spectroscopic states of MbCO. The 10μs A0→A1-3 transition involves the motion of His64, where distance between His64 and CO is found to vary up to 8.8±1.0Å during the transition of His64 from the ligand (A 1-3) to bulk solvent (A0). The His64 motion occurs within a single trajectory only once, however the multiple trajectories populate the spectroscopic A-states fully. Consequently, multiple independent molecular dynamics simulations have been found to extend biomolecular motion from 5ns of total simulation to experimental phenomena on the microsecond time scale.
UR - http://www.scopus.com/inward/record.url?scp=1942487728&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=1942487728&partnerID=8YFLogxK
U2 - 10.1016/j.jmgm.2003.12.004
DO - 10.1016/j.jmgm.2003.12.004
M3 - Article
C2 - 15099833
AN - SCOPUS:1942487728
VL - 22
SP - 369
EP - 376
JO - Journal of Molecular Graphics and Modelling
JF - Journal of Molecular Graphics and Modelling
SN - 1093-3263
IS - 5
ER -