Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2010 Jun 28;132(24):244101.
doi: 10.1063/1.3435332.

Enhanced sampling and applications in protein folding in explicit solvent

Cheng Zhang  1 , Jianpeng Ma
Affiliations

Enhanced sampling and applications in protein folding in explicit solvent

Cheng Zhang et al. J Chem Phys. .

Abstract

We report a single-copy tempering method for simulating large complex systems. In a generalized ensemble, the method uses runtime estimate of the thermal average energy computed from a novel integral identity to guide a continuous temperature-space random walk. We first validated the method in a two-dimensional Ising model and a Lennard-Jones liquid system. It was then applied to folding of three small proteins, trpzip2, trp-cage, and villin headpiece in explicit solvent. Within 0.5-1 microsecond, all three systems were reversibly folded into atomic accuracy: the alpha carbon root mean square deviations of the best folded conformations from the native states were 0.2, 0.4, and 0.4 A, for trpzip2, trp-cage, and villin headpiece, respectively.

PubMed Disclaimer

Figures

Figure 1
Figure 1
(a) Schematic illustration of auxiliary functions ϕ(β), ϕs(β) and ϕt(β), which are used in the integral identity for estimating the thermal average energy E˜(β). The estimate computed in this way uses statistics from a large temperature window (β−,β+) instead of a single bin (βi,βi+1) and also avoids systematic bias. ϕ(β) is a combination of a smooth function ϕs(β) and a function ϕt(β) localized at (βi,βi+1). ϕs(β) is controlled by two parameters a+ and a− that satisfy a++a−=1. (b) Schematic illustration of ϕs′(β) and −ϕt′(β). ϕs′(β) (shaded) spans over the whole temperature window (β−,β+) while −ϕt′(β) is localized in (βi,βi+1).
Figure 2
Figure 2
Thermodynamic quantities as functions of β for the 32×32 Ising model. Results from the method introduced in Sec. 2 are labeled as method 1 (thin solid lines, cross for errors), and those from the method introduced in Appendix B are labeled as method 2 (dashed lines, circles for errors); (a) the partition function, (b) the heat capacity, (c) the average energy, and (d) E˜(β).
Figure 3
Figure 3
Reconstructed energy distribution at a few temperatures using Eq. 20 with a window size Δβ=0.02 for the 32×32 Ising model. For comparison, we also show the resulting distributions from averaging energy distributions of two adjacent temperature bins, with the bin size δβ=0.0002. The energy distributions constructed from Eq. 20 with a larger window are more precise than those from simple averages of two adjacent bins.
Figure 4
Figure 4
(a) E˜(β) of an 864-particle Lennard-Jones system. As a comparison, the values of the average energy from several constant temperature simulations are shown as dots. (b) The reconstructed radial distribution functions of an 864-particle Lennard-Jones system at two selected temperatures T=1.0 and 2.0 are shown as points. As a comparison, the corresponding radial distribution functions from independent constant temperature simulations are shown as lines.
Figure 5
Figure 5
Trpzip2: (a) the initial fully extended conformation, (b) a typical folded structure from trajectory 2. The Cα-RMSD and heavy atom RMSD are 0.25 and 1.08 Å, respectively. Gray: reference structure (PDB ID: 1LE1).
Figure 6
Figure 6
Trpzip2: quantities along two independent simulation trajectories. Left: the first 200 ns of a fast-folding trajectory (trajectory 1). Right: 2 μs of trajectory 3. Panels from top to bottom: (a) Cα-RMSD from native structure, (b) Cα radius of gyration, (c) temperature, and (d) potential energy.
Figure 7
Figure 7
Trpzip2: (a) distribution along the RMSD from the native structure at three temperatures 300, 400, and 500 K, calculated from trajectory 3 and (b) fraction P of the folded state vs the temperature. Inset: linear fitting of log(P0∕P−1) vs β according to the two-state model.
Figure 8
Figure 8
Trpzip2: the heat capacity, computed from four independent trajectories. Bold solid line: empirical formula CV≈4.2×103∕T.
Figure 9
Figure 9
Trpzip2: joint distributions of (a) the radius of gyration vs Cα-RMSD and (b) the potential energy vs temperature. A brighter color represents a higher population density.
Figure 10
Figure 10
Trp-cage: (a) the initial fully extended structure, (b) a typical folded structure. The Cα-RMSD and the all heavy atom RMSD are 0.44 and 1.54 Å, respectively. Gray: reference structure (PDB ID: 1L2Y).
Figure 11
Figure 11
Trp-cage: three independent trajectories (a) Cα-RMSD from native structure, (b) Cα radius of gyration, (c) temperature, and (d) potential energy.
Figure 12
Figure 12
Villin headpiece: (a) the initial fully extended structure, (b) a typical folded structure compared with an NMR reference structure (gray, PDB ID: 1VII, Cα-RMSD: 1.15 Å), (c) a typical folded structure compared with an x-ray structure (gray, PDB ID: 1YRF, Cα-RMSD: 0.47 Å), and the N-terminal is not shown due to the sequence difference.
Figure 13
Figure 13
Villin headpiece: joint distributions of (a) Cα-RMSD and Cα-RMSDcore (from residues 9∼32) using the NMR structure as the reference and (b) Cα-RMSD and Cα-RMSDcore, using the x-ray structure as the reference. Statistics from the two trajectories were combined.
See this image and copyright information in PMC

References

    1. Lyubartsev A. P., Martsinovski A. A., Shevkunov S. V., and Vorontsovvelyaminov P. N., J. Chem. Phys. JCPSA6 96, 1776 (1992)10.1063/1.462133 - DOI
    2. Marinari E. and Parisi G., Europhys. Lett. EULEEJ 19, 451 (1992)10.1209/0295-5075/19/6/002 - DOI
    3. Zhang C. and Ma J., Phys. Rev. E PLEEE8 76, 036708 (2007)10.1103/PhysRevE.76.036708 - DOI - PMC - PubMed
    4. Li H., Fajer M., and Yang W., J. Chem. Phys. JCPSA6 126, 024106 (2007)10.1063/1.2424700 - DOI - PubMed
    5. Gao Y. Q., J. Chem. Phys. JCPSA6 128, 064105 (2008);10.1063/1.2825614 - DOI - PubMed
    6. Kim J., Straub J. E., and Keyes T., Phys. Rev. Lett. PRLTAO 97, 050601 (2006).10.1003/PhysRevLett.97.050601 - DOI - PubMed
    1. Zhang C. and Ma J., J. Chem. Phys. JCPSA6 130, 194112 (2009).10.1063/1.3139192 - DOI - PMC - PubMed
    1. Geyer C. J., Proceedings of the 23rd Symposium on the Interface (American Statistical Association, New York, 1991)
    2. Hansmann U. H. E., Chem. Phys. Lett. CHPLBC 281, 140 (1997)10.1016/S0009-2614(97)01198-6 - DOI
    3. Hukushima K. and Nemoto K., J. Phys. Soc. Jpn. JUPSAU 65, 1604 (1996)10.1143/JPSJ.65.1604 - DOI
    4. Swendsen R. H. and Wang J. S., Phys. Rev. Lett. PRLTAO 57, 2607 (1986);10.1103/PhysRevLett.57.2607 - DOI - PubMed
    5. Kim J., Straub J. E., and Keyes T., J. Chem. Phys. JCPSA6 132, 224107 (2010).10.1063/1.3432176 - DOI - PMC - PubMed
    1. Darve E. and Pohorille A., J. Chem. Phys. JCPSA6 115, 9169 (2001)10.1063/1.1410978 - DOI
    2. Fasnacht M., Swendsen R. H., and Rosenberg J. M., Phys. Rev. E PLEEE8 69, 056704 (2004).10.1103/PhysRevE.69.056704 - DOI - PubMed
    1. Adib A. B. and Jarzynski C., J. Chem. Phys. JCPSA6 122, 014114 (2005).10.1063/1.1829631 - DOI - PubMed

Publication types