Biomolecular dynamics at long timesteps: Bridging the timescale gap between simulation and experimentation

Tamar Schlick, Eric Barth, Margaret Mandziuk

Research output: Contribution to journalArticle

Abstract

Innovative algorithms have been developed during the past decade for simulating Newtonian physics for macromolecules. A major goal is alleviation of the severe requirement that the integration timestep be small enough to resolve the fastest components of the motion and thus guarantee numerical stability. This timestep problem is challenging if strictly faster methods with the same all-atom resolution at small timesteps are sought. Mathematical techniques that have worked well in other multiple-timescale contexts-where the fast motions are rapidly decaying or largely decoupled from others-have not been as successful for biomolecules, where vibrational coupling is strong. This review examines general issues that limit the timestep and describes available methods (constrained, reduced-variable, implicit, symplectic, multiple-timestep, and normal-mode-based schemes). A section compares results of selected integrators for a model dipeptide, assessing physical and numerical performance. Included is our dual timestep method LN, which relies on an approximate linearization of the equations of motion every Δt interval (5 fs or less), the solution of which is obtained by explicit integration at the inner timestep Δτ (e.g., 0.5 fs). LN is computationally competitive, providing 4-5 speedup factors, and results are in good agreement, in comparison to 0.5 fs trajectories. These collective algorithmic efforts help fill the gap between the time range that can be simulated and the timespans of major biological interest (milliseconds and longer). Still, only a hierarchy of models and methods, along with experimentational improvements, will ultimately give theoretical modeling the status of partner with experiment.

Original languageEnglish (US)
Pages (from-to)181-222
Number of pages42
JournalAnnual Review of Biophysics and Biomolecular Structure
Volume26
DOIs
StatePublished - 1997

Fingerprint

Dipeptides
Convergence of numerical methods
Biomolecules
Macromolecules
Linearization
Equations of motion
Physics
Trajectories
Atoms
Experiments

Keywords

  • Linearized models
  • Molecular and Langevin dynamics
  • Multiple timesteps
  • Normal modes
  • Numerical integration

ASJC Scopus subject areas

  • Biophysics
  • Structural Biology

Cite this

Biomolecular dynamics at long timesteps : Bridging the timescale gap between simulation and experimentation. / Schlick, Tamar; Barth, Eric; Mandziuk, Margaret.

In: Annual Review of Biophysics and Biomolecular Structure, Vol. 26, 1997, p. 181-222.

Research output: Contribution to journalArticle

@article{d466490df6804ed8bb9ea0b7ac5a2f1c,
title = "Biomolecular dynamics at long timesteps: Bridging the timescale gap between simulation and experimentation",
abstract = "Innovative algorithms have been developed during the past decade for simulating Newtonian physics for macromolecules. A major goal is alleviation of the severe requirement that the integration timestep be small enough to resolve the fastest components of the motion and thus guarantee numerical stability. This timestep problem is challenging if strictly faster methods with the same all-atom resolution at small timesteps are sought. Mathematical techniques that have worked well in other multiple-timescale contexts-where the fast motions are rapidly decaying or largely decoupled from others-have not been as successful for biomolecules, where vibrational coupling is strong. This review examines general issues that limit the timestep and describes available methods (constrained, reduced-variable, implicit, symplectic, multiple-timestep, and normal-mode-based schemes). A section compares results of selected integrators for a model dipeptide, assessing physical and numerical performance. Included is our dual timestep method LN, which relies on an approximate linearization of the equations of motion every Δt interval (5 fs or less), the solution of which is obtained by explicit integration at the inner timestep Δτ (e.g., 0.5 fs). LN is computationally competitive, providing 4-5 speedup factors, and results are in good agreement, in comparison to 0.5 fs trajectories. These collective algorithmic efforts help fill the gap between the time range that can be simulated and the timespans of major biological interest (milliseconds and longer). Still, only a hierarchy of models and methods, along with experimentational improvements, will ultimately give theoretical modeling the status of partner with experiment.",
keywords = "Linearized models, Molecular and Langevin dynamics, Multiple timesteps, Normal modes, Numerical integration",
author = "Tamar Schlick and Eric Barth and Margaret Mandziuk",
year = "1997",
doi = "10.1146/annurev.biophys.26.1.181",
language = "English (US)",
volume = "26",
pages = "181--222",
journal = "Annual Review of Biophysics",
issn = "1936-122X",
publisher = "Annual Reviews Inc.",

}

TY - JOUR

T1 - Biomolecular dynamics at long timesteps

T2 - Bridging the timescale gap between simulation and experimentation

AU - Schlick, Tamar

AU - Barth, Eric

AU - Mandziuk, Margaret

PY - 1997

Y1 - 1997

N2 - Innovative algorithms have been developed during the past decade for simulating Newtonian physics for macromolecules. A major goal is alleviation of the severe requirement that the integration timestep be small enough to resolve the fastest components of the motion and thus guarantee numerical stability. This timestep problem is challenging if strictly faster methods with the same all-atom resolution at small timesteps are sought. Mathematical techniques that have worked well in other multiple-timescale contexts-where the fast motions are rapidly decaying or largely decoupled from others-have not been as successful for biomolecules, where vibrational coupling is strong. This review examines general issues that limit the timestep and describes available methods (constrained, reduced-variable, implicit, symplectic, multiple-timestep, and normal-mode-based schemes). A section compares results of selected integrators for a model dipeptide, assessing physical and numerical performance. Included is our dual timestep method LN, which relies on an approximate linearization of the equations of motion every Δt interval (5 fs or less), the solution of which is obtained by explicit integration at the inner timestep Δτ (e.g., 0.5 fs). LN is computationally competitive, providing 4-5 speedup factors, and results are in good agreement, in comparison to 0.5 fs trajectories. These collective algorithmic efforts help fill the gap between the time range that can be simulated and the timespans of major biological interest (milliseconds and longer). Still, only a hierarchy of models and methods, along with experimentational improvements, will ultimately give theoretical modeling the status of partner with experiment.

AB - Innovative algorithms have been developed during the past decade for simulating Newtonian physics for macromolecules. A major goal is alleviation of the severe requirement that the integration timestep be small enough to resolve the fastest components of the motion and thus guarantee numerical stability. This timestep problem is challenging if strictly faster methods with the same all-atom resolution at small timesteps are sought. Mathematical techniques that have worked well in other multiple-timescale contexts-where the fast motions are rapidly decaying or largely decoupled from others-have not been as successful for biomolecules, where vibrational coupling is strong. This review examines general issues that limit the timestep and describes available methods (constrained, reduced-variable, implicit, symplectic, multiple-timestep, and normal-mode-based schemes). A section compares results of selected integrators for a model dipeptide, assessing physical and numerical performance. Included is our dual timestep method LN, which relies on an approximate linearization of the equations of motion every Δt interval (5 fs or less), the solution of which is obtained by explicit integration at the inner timestep Δτ (e.g., 0.5 fs). LN is computationally competitive, providing 4-5 speedup factors, and results are in good agreement, in comparison to 0.5 fs trajectories. These collective algorithmic efforts help fill the gap between the time range that can be simulated and the timespans of major biological interest (milliseconds and longer). Still, only a hierarchy of models and methods, along with experimentational improvements, will ultimately give theoretical modeling the status of partner with experiment.

KW - Linearized models

KW - Molecular and Langevin dynamics

KW - Multiple timesteps

KW - Normal modes

KW - Numerical integration

UR - http://www.scopus.com/inward/record.url?scp=0030918284&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0030918284&partnerID=8YFLogxK

U2 - 10.1146/annurev.biophys.26.1.181

DO - 10.1146/annurev.biophys.26.1.181

M3 - Article

VL - 26

SP - 181

EP - 222

JO - Annual Review of Biophysics

JF - Annual Review of Biophysics

SN - 1936-122X

ER -