1887
banner image
No data available.
Please log in to see this content.
You have no subscription access to this content.
No metrics data to plot.
The attempt to load metrics for this article has failed.
The attempt to plot a graph for these metrics has failed.
An experimentally guided umbrella sampling protocol for biomolecules
Rent:
Rent this article for
USD
10.1063/1.2976440
/content/aip/journal/jcp/129/11/10.1063/1.2976440
http://aip.metastore.ingenta.com/content/aip/journal/jcp/129/11/10.1063/1.2976440
View: Figures

Figures

Image of FIG. 1.
FIG. 1.

Schematic examples of restraining potentials for the pedagogical case of a double well potential, showing the effect of subtracting the PMF from the restraints. (a) The original double well potential. (b) A set of uneducated-guess (poor) starting potentials (dashed lines) with the negative of the PMF (solid line) and (c) the results of subtracting the PMF from these harmonic potentials (dashed lines) overlaid with the original PMF. (d) A nearly optimal set of restraining potentials and (e) the results of subtracting the PMF from these optimal harmonic potentials, overlaid with the original PMF. The potentials in (e) resulting from the subtraction (dotted lines) produce biasing restraints similar to the type we suggest to use in our guided umbrella sampling protocol.

Image of FIG. 2.
FIG. 2.

(a) “Experimental” potential of mean force for the end-to-end distance of pentapeptide with fitted biasing potentials. (b) Second derivative of the PMF calculated using a finite-difference scheme and by eliminating, with a running average filter of width , changes in the curvature occurring on length scales smaller than the thermally accessible range for the harmonic potentials. The plot has five flat regions; correspondingly, five guiding umbrella potentials are chosen with the locations of their minima within the five regions placed according to a root-mean-square best fit to the experimental PMF.

Image of FIG. 3.
FIG. 3.

Biased histograms for the pentapeptide in the case of (a) the unguided 17 window umbrella sampling, (b) 5 window unguided umbrella sampling, and (c) 5 window guided umbrella sampling after per window for each. Note the poor overlap for the histograms in panel (b) relative to (a), and the recovery of the overlap in (c).

Image of FIG. 4.
FIG. 4.

K-S test of similarity for the pentapeptide model. The cumulative probability distribution functions for the distance between the sulfur of cysteine and the center of mass of the aromatic ring of tryptophan of the pentapeptide for the three umbrella sampling methods at different simulation times were compared to the exhaustive simulations at . The K-S values for the 17 window umbrella sampling, 5 window umbrella sampling, and the 5 window guided umbrella sampling are shown. The K-S statistics for the guided umbrella sampling at various times as compared to the 17 window umbrella sampling at is also shown. The lines are meant as guides to the eye. A lower K-S value indicates higher likelihood that the distributions represent the same data.

Image of FIG. 5.
FIG. 5.

(a) Potentials of mean force for the end-to-end distance of the pentapeptide for the guiding extensive sampling at and for the umbrella sampling methods at the time they most resemble the extensive run: i.e., at for the 17 window unguided run and at for the guided run. The 5 window unguided umbrella sampling never fully resembles the extensive simulation, not even after . (b) PMFs for the 17 window unguided umbrella sampling after the full and the guided umbrella sampling after the full compared to of extensive sampling.

Image of FIG. 6.
FIG. 6.

Comparison of probability curves of pentapeptide end-to-end distance at various simulation times. (a) Exhaustive sampling at , guided umbrella sampling at ( per window), and the 17 window unguided umbrella sampling at ( per window). (b) Exhaustive sampling at , guided umbrella sampling at ( per window), and 17 window unguided umbrella sampling at ( per window). It can be seen that the umbrella sampling probabilities are beginning to shift to the right. (c) Probability distribution for the unguided 17 window umbrella sampling is shown at ( per window) with the guided umbrella sampling at 50 and . There is very little difference in the guided umbrella sampling between 50 and .

Image of FIG. 7.
FIG. 7.

Comparison of the exhaustive and guided umbrella sampling PMFs for the pentapeptide end-to-end distance as a function of simulation time. (a) Plot of exhaustive sampling times vs the guided umbrella sampling times with a similar PMF; exponential fit yields . (b) Guided umbrella sampling PMF as a function of simulation time. The time evolution of the potential of mean force of the reaction coordinate shows that the free energy well begins to narrow after about and does not converge until about . (c) Overlay of the PMFs from the exhaustive sampling (purple) and the guided umbrella sampling (red). The time of the guided umbrella sampling has been transformed using the exponential from panel (a) in order to make it line up with the exhaustive run results. The two simulations follow the same relative time path with respect to the reaction coordinate.

Image of FIG. 8.
FIG. 8.

Free energy profiles for “perpendicular” degrees of freedon, i.e., the pentapeptide dihedral angles. Exhaustive sampling distribution at is shown in green, 17 window umbrella sampling in blue, and 5 window guided umbrella sampling in red.

Image of FIG. 9.
FIG. 9.

(A) Experimental PMF and a selection of guiding umbrella sampling potentials used in the experimentally guided umbrella sampling simulations of unfolding of titin I27. As the simulations progess toward higher energy states, the force constant increases. (B) PMF for guided and unguided umbrella sampling simulations. The two methods produce substantially different curves; in contrast to the unguided runs, the guided sampling converges significantly faster, and to a profile consistent with experimental data (see text).

Image of FIG. 10.
FIG. 10.

(A) Equilibrated structure of the titin I27 domain. (B) Guided umbrella sampling titin structure for . (C) Unguided umbrella sampling structure for . The dotted lines mark the distance between carbons for residues 3 to 26 and 5 to 24, which are hydrogen bonded in the equilibrium structure. The dissociation of the A strand (in red), residues 3 to 7, is the first step in the forced unfolding pathway, which is correctly reproduced by the guided sampling, but not by the unguided one, even for twice the simulation time (see text).

Loading

Article metrics loading...

/content/aip/journal/jcp/129/11/10.1063/1.2976440
2008-09-15
2014-04-25
Loading

Full text loading...

This is a required field
Please enter a valid email address
752b84549af89a08dbdd7fdb8b9568b5 journal.articlezxybnytfddd
Scitation: An experimentally guided umbrella sampling protocol for biomolecules
http://aip.metastore.ingenta.com/content/aip/journal/jcp/129/11/10.1063/1.2976440
10.1063/1.2976440
SEARCH_EXPAND_ITEM