Calculation of collective variable-based PMF by combining WHAM with umbrella sampling

Wei Xin Xu, Yang Li, John Zhang

Research output: Contribution to journalArticle

Abstract

Potential of mean force (PMF) with respect to localized reaction coordinates (RCs) such as distance is often applied to evaluate the free energy profile along the reaction pathway for complex molecular systems. However, calculation of PMF as a function of global RCs is still a challenging and important problem in computational biology. We examine the combined use of the weighted histogram analysis method and the umbrella sampling method for the calculation of PMF as a function of a global RC from the coarse-grained Langevin dynamics simulations for a model protein. The method yields the folding free energy profile projected onto a global RC, which is in accord with benchmark results. With this method rare global events would be sufficiently sampled because the biased potential can be used for restricting the global conformation to specific regions during free energy calculations. The strategy presented can also be utilized in calculating the global intra- and intermolecular PMF at more detailed levels.

Original languageEnglish (US)
Article number068702
JournalChinese Physics Letters
Volume29
Issue number6
DOIs
StatePublished - Jun 2012

Fingerprint

sampling
free energy
profiles
biology
histograms
folding
proteins
simulation

ASJC Scopus subject areas

  • Physics and Astronomy(all)

Cite this

Calculation of collective variable-based PMF by combining WHAM with umbrella sampling. / Xu, Wei Xin; Li, Yang; Zhang, John.

In: Chinese Physics Letters, Vol. 29, No. 6, 068702, 06.2012.

Research output: Contribution to journalArticle

@article{ba0a01f76acf498bb63667aff4bb6b4e,
title = "Calculation of collective variable-based PMF by combining WHAM with umbrella sampling",
abstract = "Potential of mean force (PMF) with respect to localized reaction coordinates (RCs) such as distance is often applied to evaluate the free energy profile along the reaction pathway for complex molecular systems. However, calculation of PMF as a function of global RCs is still a challenging and important problem in computational biology. We examine the combined use of the weighted histogram analysis method and the umbrella sampling method for the calculation of PMF as a function of a global RC from the coarse-grained Langevin dynamics simulations for a model protein. The method yields the folding free energy profile projected onto a global RC, which is in accord with benchmark results. With this method rare global events would be sufficiently sampled because the biased potential can be used for restricting the global conformation to specific regions during free energy calculations. The strategy presented can also be utilized in calculating the global intra- and intermolecular PMF at more detailed levels.",
author = "Xu, {Wei Xin} and Yang Li and John Zhang",
year = "2012",
month = "6",
doi = "10.1088/0256-307X/29/6/068702",
language = "English (US)",
volume = "29",
journal = "Chinese Physics Letters",
issn = "0256-307X",
publisher = "IOP Publishing Ltd.",
number = "6",

}

TY - JOUR

T1 - Calculation of collective variable-based PMF by combining WHAM with umbrella sampling

AU - Xu, Wei Xin

AU - Li, Yang

AU - Zhang, John

PY - 2012/6

Y1 - 2012/6

N2 - Potential of mean force (PMF) with respect to localized reaction coordinates (RCs) such as distance is often applied to evaluate the free energy profile along the reaction pathway for complex molecular systems. However, calculation of PMF as a function of global RCs is still a challenging and important problem in computational biology. We examine the combined use of the weighted histogram analysis method and the umbrella sampling method for the calculation of PMF as a function of a global RC from the coarse-grained Langevin dynamics simulations for a model protein. The method yields the folding free energy profile projected onto a global RC, which is in accord with benchmark results. With this method rare global events would be sufficiently sampled because the biased potential can be used for restricting the global conformation to specific regions during free energy calculations. The strategy presented can also be utilized in calculating the global intra- and intermolecular PMF at more detailed levels.

AB - Potential of mean force (PMF) with respect to localized reaction coordinates (RCs) such as distance is often applied to evaluate the free energy profile along the reaction pathway for complex molecular systems. However, calculation of PMF as a function of global RCs is still a challenging and important problem in computational biology. We examine the combined use of the weighted histogram analysis method and the umbrella sampling method for the calculation of PMF as a function of a global RC from the coarse-grained Langevin dynamics simulations for a model protein. The method yields the folding free energy profile projected onto a global RC, which is in accord with benchmark results. With this method rare global events would be sufficiently sampled because the biased potential can be used for restricting the global conformation to specific regions during free energy calculations. The strategy presented can also be utilized in calculating the global intra- and intermolecular PMF at more detailed levels.

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

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

U2 - 10.1088/0256-307X/29/6/068702

DO - 10.1088/0256-307X/29/6/068702

M3 - Article

AN - SCOPUS:84862689513

VL - 29

JO - Chinese Physics Letters

JF - Chinese Physics Letters

SN - 0256-307X

IS - 6

M1 - 068702

ER -