Markovian milestoning with Voronoi tessellations

Eric Vanden Eijnden, Maddalena Venturoli

Research output: Contribution to journalArticle

Abstract

A new milestoning procedure using Voronoi tessellations is proposed. In the new procedure, the edges of Voronoi cells are used as milestones, and the necessary kinetic information about the transitions between the milestones is calculated by running molecular dynamics (MD) simulations restricted to these cells. Like the traditional milestoning technique, the new procedure offers a reduced description of the original dynamics and permits to efficiently compute the various quantities necessary in this description. However, unlike traditional milestoning, the new procedure does not require to reinitialize trajectories from the milestones, and thereby it avoids the approximation made in traditional milestoning that the distribution for reinitialization is the equilibrium one. In this paper we concentrate on Markovian milestoning, which we show to be valid under suitable assumptions, and we explain how to estimate the rate matrix of transitions between the milestones from data collected from the MD trajectories in the Voronoi cells. The rate matrix can then be used to compute mean first passage times between milestones and reaction rates. The procedure is first illustrated on test-case examples in two dimensions and then applied to study the kinetics of protein insertion into a lipid bilayer by means of a coarse-grained model.

Original languageEnglish (US)
Article number194101
JournalJournal of Chemical Physics
Volume130
Issue number19
DOIs
StatePublished - 2009

Fingerprint

Molecular dynamics
Trajectories
Lipid bilayers
Kinetics
cells
trajectories
molecular dynamics
Reaction rates
Molecular Dynamics Simulation
kinetics
lipids
insertion
Computer simulation
reaction kinetics
proteins
Lipid Bilayers
Proteins
estimates
matrices
approximation

ASJC Scopus subject areas

  • Physics and Astronomy(all)
  • Physical and Theoretical Chemistry
  • Medicine(all)

Cite this

Markovian milestoning with Voronoi tessellations. / Vanden Eijnden, Eric; Venturoli, Maddalena.

In: Journal of Chemical Physics, Vol. 130, No. 19, 194101, 2009.

Research output: Contribution to journalArticle

@article{45d6c34c624e452a90fdcd83f04d5fc5,
title = "Markovian milestoning with Voronoi tessellations",
abstract = "A new milestoning procedure using Voronoi tessellations is proposed. In the new procedure, the edges of Voronoi cells are used as milestones, and the necessary kinetic information about the transitions between the milestones is calculated by running molecular dynamics (MD) simulations restricted to these cells. Like the traditional milestoning technique, the new procedure offers a reduced description of the original dynamics and permits to efficiently compute the various quantities necessary in this description. However, unlike traditional milestoning, the new procedure does not require to reinitialize trajectories from the milestones, and thereby it avoids the approximation made in traditional milestoning that the distribution for reinitialization is the equilibrium one. In this paper we concentrate on Markovian milestoning, which we show to be valid under suitable assumptions, and we explain how to estimate the rate matrix of transitions between the milestones from data collected from the MD trajectories in the Voronoi cells. The rate matrix can then be used to compute mean first passage times between milestones and reaction rates. The procedure is first illustrated on test-case examples in two dimensions and then applied to study the kinetics of protein insertion into a lipid bilayer by means of a coarse-grained model.",
author = "{Vanden Eijnden}, Eric and Maddalena Venturoli",
year = "2009",
doi = "10.1063/1.3129843",
language = "English (US)",
volume = "130",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "19",

}

TY - JOUR

T1 - Markovian milestoning with Voronoi tessellations

AU - Vanden Eijnden, Eric

AU - Venturoli, Maddalena

PY - 2009

Y1 - 2009

N2 - A new milestoning procedure using Voronoi tessellations is proposed. In the new procedure, the edges of Voronoi cells are used as milestones, and the necessary kinetic information about the transitions between the milestones is calculated by running molecular dynamics (MD) simulations restricted to these cells. Like the traditional milestoning technique, the new procedure offers a reduced description of the original dynamics and permits to efficiently compute the various quantities necessary in this description. However, unlike traditional milestoning, the new procedure does not require to reinitialize trajectories from the milestones, and thereby it avoids the approximation made in traditional milestoning that the distribution for reinitialization is the equilibrium one. In this paper we concentrate on Markovian milestoning, which we show to be valid under suitable assumptions, and we explain how to estimate the rate matrix of transitions between the milestones from data collected from the MD trajectories in the Voronoi cells. The rate matrix can then be used to compute mean first passage times between milestones and reaction rates. The procedure is first illustrated on test-case examples in two dimensions and then applied to study the kinetics of protein insertion into a lipid bilayer by means of a coarse-grained model.

AB - A new milestoning procedure using Voronoi tessellations is proposed. In the new procedure, the edges of Voronoi cells are used as milestones, and the necessary kinetic information about the transitions between the milestones is calculated by running molecular dynamics (MD) simulations restricted to these cells. Like the traditional milestoning technique, the new procedure offers a reduced description of the original dynamics and permits to efficiently compute the various quantities necessary in this description. However, unlike traditional milestoning, the new procedure does not require to reinitialize trajectories from the milestones, and thereby it avoids the approximation made in traditional milestoning that the distribution for reinitialization is the equilibrium one. In this paper we concentrate on Markovian milestoning, which we show to be valid under suitable assumptions, and we explain how to estimate the rate matrix of transitions between the milestones from data collected from the MD trajectories in the Voronoi cells. The rate matrix can then be used to compute mean first passage times between milestones and reaction rates. The procedure is first illustrated on test-case examples in two dimensions and then applied to study the kinetics of protein insertion into a lipid bilayer by means of a coarse-grained model.

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

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

U2 - 10.1063/1.3129843

DO - 10.1063/1.3129843

M3 - Article

C2 - 19466815

AN - SCOPUS:66049102550

VL - 130

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 19

M1 - 194101

ER -