| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, People's Republic of China; and The Department of Chemistry and Department of Physics, State University of New York, Stony Brook, New York
Correspondence: Address reprint requests to Jin Wang, E-mail: jinwang{at}sprynet.com.
| ABSTRACT |
|---|
|
|
|---|
6)), the kinetics is an exponential Poissonian process again. | INTRODUCTION |
|---|
|
|
|---|
Both theoretical and experimental investigations on folding and reaction kinetics show complex kinetics in different ranges of temperature (Bryngelson and Wolynes, 1989
; Bryngelson et al., 1995
; Chan and Dill, 1994
; Abkevich et al., 1994
; Wang et al., 1996; Gutin et al., 1996
; Cieplak et al., 1999
; Seno et al., 1998
; Klimov and Thirumalai, 1998
; Kaya and Chan, 2000
, 2002
; Itzhaki et al., 1995
; Schuler et al., 2002
; Lipman et al., 2003
; Sabelko et al., 1999
; Nguyen et al., 2003
; Frauenfelder et al., 1988
, 1991
; Yang and Xie, 2002a
,b
). By varying the temperature, the underlying energy landscape structures can be probed in different levels, from the global to the local detail perspectives (Frauenfelder et al., 1988
, 1991
). The relationship between dynamics and functions of the biomolecules can be revealed. Although different theoretical approaches explain kinetic behavior within specific temperature ranges, the unified picture of kinetics of the whole temperature range seems still lacking. This is the purpose of the current study. Lately, the diffusive dynamics of folding is shown to have complex kinetics (Nguyen et al., 2003
; Lee et al., 2003
; Zhou et al., 2003
; V. B. P. Leite, J. N. Onuchic, G. Stell, and J. Wang, unpublished results).
In this article, we study the kinetics in the whole temperature range. We show that the Poisson(exponential)-non-Poisson(non-exponential)-Poisson (exponential) kinetics emerge from high to low temperatures. This phenomena seems to be universal not only for protein folding, but might also exist in other biomolecular folding, biomolecular binding and reaction systems, electron transfer, viscous liquid, and glassy materials.
The current results of the study are also relevant to the single molecule studies where the mean of the observables is often unreliable due to the large statistical fluctuations (which are not smeared out by the number of molecules as in the bulk case), and cannot be used to accurately characterize the system. In general, fluctuations of the observables intrinsic for characterizing the system are obtained from the information on moments and distributions (Wang and Wolynes, 1995
, 1999
; Onuchic et al., 1999
; Wang, 2003
). The connections of the theory and simulations with the single molecule kinetic experiments can be made through the analysis of the long-time dynamic trajectories or multiple short-time runs where information on the mean, high-order moments and distributions, or histograms of the important observables can be extracted (Lu et al., 1998
; Schenter et al., 1999
; Moerner, 1996
; Zhuang et al., 2000
, 2002
; Jia et al., 1999
).
| MATERIALS AND METHODS |
|---|
|
|
|---|
as an order parameter to represent the folding progress. The system is assumed to be in quasiequilibrium with respect to
, and the states are kinetically locally connected. This is a good approximation when folding has a definite free energy barrier between non-native and native states. It might not be as good for the downhill case, where there is no barrier. The results we obtain in this article, however, seem not to be influenced qualitatively by the kinetic connectivity assumption. It has been shown that in the global kinetic connectivity case, similar kinetic behavior occurs (Saven et al., 1994
In this model there are N residues in a polypeptide chain. For each residue there are
+ 1 available conformational states, one being the native state. A simplified version of the polypeptide chain energy is expressed as
![]() | (1) |
i is the state of ith residue. The three terms represent the one-body potential, two-body interactions for nearest-neighbor residues in sequence, and interactions for residues close in space but not in sequence, respectively. Due to the sequence heterogeneity, the energies and interactions can be approximated by random variables of Gaussian distributions (Derrida, 1981
![]() | (2) |
) is the average free energy for the polypeptide chain. T is a scaled temperature,
+ 1 is the number of conformational states of each residue, and 
and
L are energy differences between the native and average non-native states for one- and two-body interactions (energy gap biased to native state), respectively. 
and
L are energy spreads of one- and two-body non-native interactions (roughness of the landscape). Note that the two-body energies
L and
L include contributions from the second and third terms in Eq. 1. The last term in F is the configurational entropy contribution.
The kinetic process along the above free-energy landscape is approximated via the use of Metropolis rate dynamics. Using continuous-time random walks, the generalized Fokker-Planck diffusion equation in the Laplace-transformed space can be obtained (Bryngelson and Wolynes, 1989
; Lee et al., 2003
),
![]() | (3) |
, s)
F(
)/T + log [D(
, s)/D(
, 0)]. In Eq. 3, s is the Laplace transform variable over time
.
is the Laplace transform of G(
,
), the probability density function. G(
,
)d
gives the probability for a polypeptide chain to stay between
and
+ d
at time
. The value ni(
) is the initial condition for G(
,
). D(
, s) is the frequency-dependent diffusion coefficient (Bryngelson and Wolynes, 1989
![]() | (4) |
(
)
1/
+ (1 1/
)
. The average
...
R is taken over P(R,
), the probability distribution function of transition rate R from one state with order parameter
to its neighboring states, which may have order parameters equal to
,
, or
The explicit expression of P(R,
) can be found in Bryngelson and Wolynes (1989)
=
i and an absorbing one at
=
f. The choice of an absorbing boundary condition at
=
f facilitates our calculation for the first passage time and its distribution.
One can rewrite Eq. 3 in its integral-equation representation by integrating it twice over
:
![]() | (5) |
In this work we mainly study the behavior of the first passage time (FPT) for the order parameter to reach
f. This FPT characterizes the folding time. By taking the derivatives with respect to s in Eq. 5 and taking the limit of s = 0, we can iteratively obtain the moments of the first passage time.
![]() |
) is the distribution of the first passage time. When n = 1, the mean first passage time is given as
![]() |
, and by the observation that the distribution of the first passage time
, where
and
are Laplace transforms of PFPT(
) and
(
) (
), respectively, we can obtain the information of PFPT(
) by studying the behavior of
Due to the fact that Eq. 5 is linear in G(
, s), we can solve it with the numerical matrix-inversion technique. | RESULTS |
|---|
|
|
|---|
= 10 to match the physical scales. For simplicity we assume 
=
L and 
=
L. The ratio of the energy gap between the native state and the average of non-native states over the spread of non-native states, 
/
, representing the energy bias or slope toward the native state relative to the spread (variance) or the fluctuations in energy of the polypeptide chain on the folding landscape, can be shown to be the controlling parameter in this problem. The large values of 
/
imply a steep folding funnel and small values of 
/
imply a rough folding landscape. We set the initial distribution to be ni(
) =
(
i), where
i is set to be 0.05. In our calculations we set
f = 0.95.
The mean first passage time (MFPT) 

for the folding process versus a scaled inverse temperature, T0/T (T0 is defined as the temperature of the minimum or optimal (fastest) mean first passage time), is plotted in Fig. 1 for several settings of the parameters 
/
, ranging from the more funneled energy landscape to the more rough energy landscape. Notice that the energy gap 
and roughness of the landscape 
are dependent on both internal sequence compositions of the protein and external environments such as solvents, denaturants, pressure, etc. We have a U-shape curve for each fixed 
/
, and the MFPT reaches its minimum at temperature T0. At high temperatures, the MFPT is large although the diffusion process itself is fast (i.e., D(
, s) is large). This long-time folding behavior is due to the instability of the native state. The MFPT is also large at low temperature, which indicates that the polypeptide chain is trapped in low-energy non-native states. The kinetic diffusion is very slow. This is in agreement with simulation studies and the experiments (this chevron rollover phenomena was first investigated and explained by Miller et al. (1992)
, Socci and Onuchic (1994
, 1995
), Socci et al. (1996)
, Gutin et al. (1996)
, Itzhaki et al. (1995)
, Cieplak et al. (1999)
, Seno et al. (1998)
, Klimov and Thirumalai (1998)
, Kaya and Chan (2000
, 2002
, 2003
), Chan et al. (2004)
, Plotkin and Wolynes (1998)
, Plotkin and Onuchic (2002a
,b
), Zhou et al. (2003)
, and Nguyen et al. (2003)
.
|
|
|

2
/

2 versus inverse kinetic minimum (fastest folding) scaled temperature. We can see that at high temperatures (T0/T < 1 in this case), the ratio is a constant close to 2 (for example, the ratio is equal to 2.7 for 
/
= 2.8). Notice that for the Poissonian process, 
n
= n!

n. This implies an approximately underlying Poissonian statistics. The kinetic process can thus be shown as exponential (Wang, 2003
6, the range of values is due to different landscape parameters: gap/roughness ratios), the second-order moment ratios decrease rapidly to values close to 2. This implies an approximate Poissonian process and exponential kinetics again.
|

2
/

2 versus inverse folding scaled temperature Tf/T. We can see above the folding transition temperature, the second moment ratio is close to 2. This implies Poissonian exponential kinetics. Below Tf and above Tk, the second moment ratio is still close to 2. It implies again Poissonian exponential kinetics. The non-exponential kinetics only emerges when the temperature is below Tk. As the temperature drops even lower (Tf/T > 4.5
7.5), the second-order moment ratios drop rapidly to values close to 2, implying Poissonian exponential kinetics again.
|

/
). We found that in all the cases, T0 < Tk < Tf. Since Tk is the marker for complex kinetics, this leads to exponential kinetics at temperatures T > Tk. The non-exponential kinetics quickly emerges at temperatures T0 < T < Tk (and continue to temperatures at (1/6
1/4) T0 < T < T0). As the gap/roughness ratio increases (
/
), the ratio of folding temperature to kinetic transition temperature Tk (where kinetics is switched from exponential to non-exponential) increases. This indicates, as the underlying folding energy landscape is more funneled toward native state relative to the local traps, that there is a wider temperature window at Tk < T < Tf for the exponential kinetics. In other words, it is relatively easier for folding (increasing Tf) and relatively harder for kinetic transition (decreasing Tk) to occur. Notice that Tk/T0 is almost a constant. T0 represents the onset of the rollover behavior in kinetics and Tk represents the onset of the non-exponential kineticsthe transition from a "smoother" appearing folding landscape to a "rougher" or "bumpier" landscape, so that basins of attraction of part of the landscape start to be partially frozen or become traps. Notice that the ratio of Tk and T0 is almost a constant. The kinetic behavior around (Tk) is analogous to the glassy material at TA above the glass transition temperature where partial freezing occurs (Kirkpatrick and Wolynes, 1987
|

/
= 4.0). The fact the curve looks like a straight line implies that the distribution of first passage time in Laplace space is approximately a stretched exponential (
) which is the form of Laplace transform of the Lévy distribution in the time space. So we have
![]() |
for large
, approaching power law distribution at long times. The power law exponent is linear with the slope of the above figure. It can be shown that power law exponent is a monotonically increasing function of the temperature below T0, implying that the tail of the FPT distribution becomes fattier (u is smaller and the power law decay is slower) as the temperature drops lower.
|
| DISCUSSIONS |
|---|
|
|
|---|
6 in this case), more and more traps become important. One can expand the Gaussian-like density of states near a frozen one, resulting in the linearization in energy of the exponential. The density of states in this temperature range thus becomes exponentially distributed. Since the first passage time is exponentially dependent on the barrier (the energy), then the distribution of first passage time follows a power law (
) for certain low temperature regimes at long times. The tail of the distribution of the first passage time becomes fattier as the temperature decreases. The rare kinetic events can play an important role. Specific kinetic path might be dominating. At very low temperature (lower then T0/T > 4
6 in this case), only very few states are kinetically accessible. The kinetics therefore is dominated with the energy barrier in the deepest valley the system is trapped into, resulting to exponential process and Poissonian statistics again.
Four temperature scales appear in this study, all representing different scales or levels of the landscape. Tf is the folding transition temperature. The average kinetics characterized by the mean first passage time has a U-shape dependence on temperature, with the fastest time at temperature T0. The fluctuations of the kinetics measured by the ratios of the second order moment to the square of the first-order moment has a bell-shape dependence on temperature, with the turning point at high temperature side at Tk. At very low temperatures, T0/T > 4
6, the system is frozen (to single traps). It is important to notice that the point of the chevron rollover (the fast folding time), T0 is related but different in value from the onset of the complex kinetics (exponential-non-exponential transition), Tk. In other words, although the turning point of the average kinetics is at T0, the turning point of the fluctuations of the kinetics or the transition from exponential to non-exponential kinetics is at Tk (on the high temperature side). Since Tk > T0 (in fact Tf > Tk > T0 > T0/(4
6)), the fluctuations in kinetics first sense the traps or bumps of the landscape. In other words, the fluctuations are more sensitive toward the shape of the landscape. They can be used to probe the underlying structure.
It is worth mentioning that although we focus on the study of the protein folding problem in this article, the approach we use here is very general for treating problems with barrier crossings on a complex energy landscape. In fact, several experiments on folding (Sabelko et al., 1999
; Nguyen et al., 2003
), binding (Frauenfelder et al., 1988
, 1991
), and reaction-conformational dynamics (Yang and Xie, 2002a
,b
) already show the existence of complex kinetics in different temperature ranges. With the rapid advances in the experiments as well as the computational power, the dynamic trajectories (long-time or multiple short-time) can be obtained and analyzed relatively easily compared with those done before (V. B. P. Leite, J. N. Onuchic, G. Stell, and J. Wang, unpublished results). It will be interesting to see the comparisons with the analytical results obtained here in the wide temperature ranges to reveal the intrinsic features and topology of the underlying energy landscape.
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Submitted on March 12, 2004; accepted for publication June 28, 2004.
| REFERENCES |
|---|
|
|
|---|
Bryngelson, J. D., and P. G. Wolynes. 1989. Intermediates and barrier crossing in a random energy model (with applications to protein folding). J. Phys. Chem. 93:69026915.[CrossRef]
Bryngelson, J. D., J. O. Onuchic, N. D. Socci, and P. G. Wolynes. 1995. Funnels, pathways, and the energy landscape of protein-foldinga synthesis. Proteins Struct. Funct. Genet. 21:167195.[CrossRef][Medline]
Chan, H. S., and K. A. Dill. 1994. Transition states and folding dynamics of proteins and heteropolymers. J. Chem. Phys. 100:92389257.[CrossRef]
Chan, H. S., S. Shimizu, and H. Kaya. 2004. Cooperativity principles in protein folding. Meth. Enzymol. 380:350379.[Medline]
Cieplak, M., T. X. Hoang, and M. S. Li. 1999. Scaling of folding properties in simple models of proteins. Phys. Rev. Lett. 83:16841687.[CrossRef]
Derrida, B. 1981. Random-energy model: an exactly solvable model of disordered systems. Phys. Rev. B. 24:26132626.[CrossRef]
Frauenfelder, H., F. Parak, and R. D. Young. 1988. Conformational substrates in proteins. Annu. Rev. Biophys. Biophys. Chem. 17:451479.[CrossRef][Medline]
Frauenfelder, H., S. G. Sligar, and P. G. Wolynes. 1991. The energy landscapes and motions of proteins. Science. 254:15981603.
Gutin, A. M., V. I. Abkevich, and E. I. Shakhnovich. 1996. Chain-length scaling of protein folding time. Phys. Rev. Lett. 77:54335436.[CrossRef][Medline]
Itzhaki, L. S., D. E. Otzen, and A. R. Fersht. 1995. The structure of the transition state for folding of chymotrypsin inhibitor-2 analyzed by protein engineering methodsevidence for a nucleation-condensation mechanism for protein-folding. J. Mol. Biol. 254:260288.[CrossRef][Medline]
Jia, Y. W., D. S. Talaga, W. L. Lau, H. S. M. Lu, W. F. DeGrado, and R. M. Hochstrasser. 1999. Folding dynamics of single GCN-4 peptides by fluorescence resonant energy transfer confocal microscopy. Chem. Phys. 247:6983.[CrossRef]
Kaya, H., and H. S. Chan. 2000. Energetic components of cooperative protein folding. Phys. Rev. Lett. 85:48234826.[CrossRef][Medline]
Kaya, H., and H. S. Chan. 2002. Towards a consistent modeling of protein thermodynamic and kinetic cooperativity: how applicable is the transition state picture to folding and unfolding? J. Mol. Biol. 315:899909.[CrossRef][Medline]
Kaya, H., and H. S. Chan. 2003. Contact order dependent protein folding rates: kinetic consequences of a cooperative interplay between favorable nonlocal interactions and local conformational preferences. Proteins. 52:524533.[CrossRef][Medline]
Kirkpatrick, T. R., and P. G. Wolynes. 1987. Stable and metastable states in mean-field Potts and structural glasses. Phys. Rev. B. 36:85528564.[CrossRef]
Kirkpatrick, T. R., D. Thirumalai, and P. G. Wolynes. 1989. Scaling concepts for the dynamics of viscous liquids near an ideal glassy state. Phys. Rev. A. 40:10451054.[CrossRef][Medline]
Klimov, D. K., and D. Thirumalai. 1998. Linking rates of folding in lattice models of proteins with underlying thermodynamic characteristics. J. Chem. Phys. 109:41194125.[CrossRef]
Lee, C. L., C. T. Lin, G. Stell, and J. Wang. 2003. Diffusion dynamics, moments, and distribution of first-passage time on the protein-folding energy landscape, with applications to single molecules. Phys. Rev. E. 67:041905.[CrossRef]
Lee, C. L., G. Stell, and J. Wang. 2003. First-passage time distribution and non-Markovian diffusion dynamics of protein folding. J. Chem. Phys. 118:959968.[CrossRef]
Levinthal, C. 1969. Proceedings in Mössbauer Spectroscopy in Biological Systems. P. DeBrunner, J. Tsibris, and E. Munck, editors. University of Illinois Press, Urbana, IL. p22.
Lipman, E. A., B. Schuler, O. Bakajin, and W. A. Eaton. 2003. Single-molecule measurement of protein folding kinetics. Science. 301:12331235.
Lu, H. P., L. Y. Xun, and X. S. Xie. 1998. Single-molecule enzymatic dynamics. Science. 282:18771882.
Moerner, W. E. 1996. High-resolution optical spectroscopy of single molecules in solids. Acc. Chem. Res. 29:563571.[CrossRef]
Miller, R., C. A. Danko, M. J. Fasolka, A. C. Balazs, H. S. Chan, and K. A. Dill. 1992. Folding kinetics of proteins and copolymers. J. Chem. Phys. 96:768780.[CrossRef]
Nguyen, H., M. Jäger, A. Moretto, M. Gruebele, and J. W. Kelly. 2003. Tuning the free-energy landscape of a WW domain by temperature, mutation, and truncation. Proc. Natl. Acad. Sci. USA. 100:39483953.
Onuchic, J. N., J. Wang, and P. G. Wolynes. 1999. Analyzing single molecule trajectories on complex energy landscapes using replica correlation functions. Chem. Phys. 247:175184.[CrossRef]
Plotkin, S. S., and P. G. Wolynes. 1998. Non-Markovian configurational diffusion and reaction coordinates for protein folding. Phys. Rev. Lett. 80:50155018.[CrossRef]
Plotkin, S. S., and J. N. Onuchic. 2002a. Understanding protein folding with energy landscape theory. I. Basic concepts. Quar. Rev. Biophys. 35:111167.[CrossRef][Medline]
Plotkin, S. S., and J. N. Onuchic. 2002b. Understanding protein folding with energy landscape theory. II. Quantitative aspects. Quar. Rev. Biophys. 35:205286.[CrossRef][Medline]
Sabelko, J., J. Ervin, and M. Gruebele. 1999. Stretching lattice models of protein folding. Proc. Natl. Acad. Sci. USA. 96:60316036.
Saven, J. G., J. Wang, and P. G. Wolynes. 1994. Kinetics of protein folding: the dynamics of globally connected rough energy landscapes with biases. J. Chem. Phys. 101:1103711043.[CrossRef]
Seno, F., C. Micheletti, A. Maritan, and J. R. Banavar. 1998. Variational approach to protein design and extraction of interaction potentials. Phys. Rev. Lett. 81:21722175.[CrossRef]
Schenter, G. K., H. P. Lu, and X. S. Xie. 1999. Statistical analyses and theoretical models of single-molecule enzymatic dynamics. J. Phys. Chem. A. 103:1047710488.[CrossRef]
Schuler, B., E. A. Lipman, and W. A. Eaton. 2002. Probing the free-energy surface for protein folding with single-molecule fluorescence spectroscopy. Nature. 419:743747.[CrossRef][Medline]
Socci, N. D., and J. N. Onuchic. 1994. Folding kinetics of proteinlike heteropolymers. J. Chem. Phys. 101:15191528.[CrossRef]
Socci, N. D., and J. N. Onuchic. 1995. Kinetic and thermodynamic analysis of proteinlike heteropolymers: Monte Carlo histogram technique. J. Chem. Phys. 103:47324744.[CrossRef]
Socci, N. D., J. N. Onuchic, and P. G. Wolynes. 1996. Diffusive dynamics of the reaction coordinate for protein folding funnels. J. Chem. Phys. 104:58605868.[CrossRef]
Wang, J., and P. G. Wolynes. 1995. Intermittency of single molecule reaction dynamics in fluctuating environments. Phys. Rev. Lett. 74:43174320.[CrossRef][Medline]
Wang, J., J. Onuchic, and P. Wolynes. 1996a. Statistics of kinetic pathways on biased rough energy landscapes with applications to protein folding. Phys. Rev. Lett. 76:48614864.[CrossRef][Medline]
Wang, J., J. G. Saven, and P. G. Wolynes. 1996b. Kinetics in a globally connected, correlated random energy model. J. Chem. Phys. 105:1127611284.[CrossRef]
Wang, J., and P. G. Wolynes. 1999. Intermittency of activated events in single molecules: the reaction diffusion description. J. Chem. Phys. 110:48124819.[CrossRef]
Wang, J. 2003. Statistics, pathways and dynamics of single molecule protein folding. J. Chem. Phys. 118:952958.[CrossRef]
Yang, H., and X. S. Xie. 2002a. Statistical approaches for probing single-molecule dynamics photon by photon. Chem. Phys. 284:423437.[CrossRef]
Yang, H., and X. S. Xie. 2002b. Probing single-molecule dynamics photon-by-photon. J. Chem. Phys. 117:1096510979.[CrossRef]
Zhou, Y., C. Zhang, G. Stell, and J. Wang. 2003. Temperature dependence of the distribution of the first passage time: results from discontinuous molecular dynamics simulations of an all-atom model of the second-hairpin fragment of protein G. J. Am. Chem. Soc. 125:63006305.[CrossRef][Medline]
Zhuang, X. W., L. E. Bartley, H. P. Babcock, R. Russell, T. Ha, D. Herschlag, and S. Chu. 2000. A single-molecule study of RNA catalysis and folding. Science. 288:20482051.
Zhuang, X. W., H. Kim, M. J. B. Pereira, H. P. Babcock, N. G. Walter, and S. Chu. 2002. Correlating structural dynamics and function in single ribozyme molecules. Science. 296:14731476.
This article has been cited by other articles:
![]() |
G. Baldini, F. Cannone, G. Chirico, M. Collini, B. Campanini, S. Bettati, and A. Mozzarelli Evidence of Discrete Substates and Unfolding Pathways in Green Fluorescent Protein Biophys. J., March 1, 2007; 92(5): 1724 - 1731. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Kaya, Z. Liu, and H. S. Chan Chevron Behavior and Isostable Enthalpic Barriers in Protein Folding: Successes and Limitations of Simple Go-like Modeling Biophys. J., July 1, 2005; 89(1): 520 - 535. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. B. P. Leite, J. N. Onuchic, G. Stell, and J. Wang Probing the Kinetics of Single Molecule Protein Folding Biophys. J., December 1, 2004; 87(6): 3633 - 3641. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |