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.
Diffusive dynamics on multidimensional rough free energy surfaces
Rent this article for


Image of FIG. 1.
FIG. 1.

(a) Coordinate system for CO rebinding in Mb. The center of mass of the CO molecule is at position with a distance to the center of mass of the four pyrrole nitrogen atoms which is at (0,0,0). The projection of defines the angle . The vector , which is used in the definition of , is given by the direction of the projection of and the length of . Note that . (b) is the free energy surface from the present simulations (see also Fig. 2) with the position of the Fe atom, the state, and the Xe4 pocket labeled. (c) shows the projection of the position of the CO center of mass onto the heme plane drawn together with a wire frame of the heme unit from previous MD simulations of unbound CO (Ref. 20). The projection in (c) is rotated to illustrate that representations (b) and (c) cover a very similar phase space. The earlier simulations were too short to sample the second state Xe4(2) (see Fig. 2).

Image of FIG. 2.
FIG. 2.

(Color) (a) Contour map of FES for CO rebinding to native myoglobin. The state is the bound state, whereas DP (distal pocket, state), Xe4(1), and Xe4(2) correspond to unbound CO. The ruggedness of the FES can be clearly seen. Contour lines are given in kcal/mol. (b) Minimum energy path along the yellow arrow in panel (a). (c) Probability distribution of CO for starting from (centered at locations “X”). At , the bound state is already populated, whereas much of the population still resides in the Xe4 pocket. Red contours are low probability, blue contours are high probability, and green contour is intermediate. (d) Active site region of MbCO with CO in the state (1), at the transition between DP and Xe4 (2), in Xe4(1) (3), and in Xe4(2) (4). Residues that form the top and bottom of the Xe4 cavity are colored blue and gold, respectively.

Image of FIG. 3.
FIG. 3.

Free energy surface as a function of (number of native contacts) and (the radius of gyration in Å). (a) original FES; (b) FES with a modified barrier between the unfolded and the folded state. Contours are separated by .

Image of FIG. 4.
FIG. 4.

(Color online) Folding time (ps) for protein G depending on the transition barrier and the diffusion coefficient : the green curve (solid black line) corresponds to , the blue line (dot dashed black) corresponds to the -dependent diffusion coefficient, the red line (solid grey) corresponds to , and the black line to . The circles are calculated folding times and the continuous lines are the fit. Dashed lines mark the experimentally observed folding times: exponential phase with a time scale of and the slower, denaturant-dependent time constant (Ref. 32).


Generic image for table
Table I.

The positions and energies of the minima in the unbound state on the FES for , the corresponding center-of-mass distance (Å), and the values of the forward and backward free energy barriers (in kcal/mol) between the states indicated.

Generic image for table
Table II.

The rebinding times (ns) for CO depending on the initial distribution or for different values of (kcal/mol).


Article metrics loading...


Full text loading...

This is a required field
Please enter a valid email address
752b84549af89a08dbdd7fdb8b9568b5 journal.articlezxybnytfddd
Scitation: Diffusive dynamics on multidimensional rough free energy surfaces