Coordinate-dependent diffusion in protein folding

The origin of nonmonotonic complex behavior and the effects of nonnative interactions on the diffusive properties of protein folding

We present a method for calculating the configurational-dependent diffusion coefficient of a globular protein as a function of the global folding process. Using a coarse-grained structure-based model, we determined the diffusion coefficient, in reaction coordinate space, as a function of the fraction of native contacts formed Q for the cold shock protein (TmCSP). We find nonmonotonic behavior for the diffusion coefficient, with high values for the folded and unfolded ensembles and a lower range of values in the transition state ensemble. We also characterized the folding landscape associated with an energetically frustrated variant of the model. We find that a low-level of frustration can actually stabilize the native ensemble and increase the associated diffusion coefficient. These findings can be understood from a mechanistic standpoint, in that the transition state ensemble has a more homogeneous structural content when frustration is present. Additionally, these findings are consistent with earlier calculations based on lattice models of protein folding and more recent single-molecule fluorescence measurements.

Protein folding pathways and state transitions described by classical equations of motion of an elastic network model

Protein topology defined by the matrix of residue contacts has proved to be a fruitful basis for the study of protein dynamics. The widely implemented coarse-grained elastic network model of backbone fluctuations has been used to describe crystallographic temperature factors, allosteric couplings, and some aspects of the folding pathway. In the present study, we develop a model of protein dynamics based on the classical equations of motion of a damped network model (DNM) that describes the folding path from a completely unfolded state to the native conformation through a single-well potential derived purely from the native conformation. The kinetic energy gained through the collapse of the protein chain is dissipated through a friction term in the equations of motion that models the water bath. This approach is completely general and sufficiently fast that it can be applied to large proteins. Folding pathways for various proteins of different classes are described and shown to correlate with experimental observations and molecular dynamics and Monte Carlo simulations. Allosteric transitions between alternative protein structures are also modeled within the DNM through an asymmetric double-well potential.

Dynamics of protein folding: probing the kinetic network of folding-unfolding transitions with experiment and theory

The problem of spontaneous folding of amino acid chains into highly organized, biologically functional three-dimensional protein structures continues to challenge the modern science. Understanding how proteins fold requires characterization of the underlying energy landscapes as well as the dynamics of the polypeptide chains in all stages of the folding process. In recent years, important advances toward these goals have been achieved owing to the rapidly growing interdisciplinary interest and significant progress in both experimental techniques and theoretical methods. Improvements in the experimental time resolution led to determination of the timescales of the important elementary events in folding, such as formation of secondary structure and tertiary contacts. Sensitive single molecule methods made possible probing the distributions of the unfolded and folded states and following the folding reaction of individual protein molecules. Discovery of proteins that fold in microseconds opened the possibility of atomic-level theoretical simulations of folding and their direct comparisons with experimental data, as well as of direct experimental observation of the barrier-less folding transition. The ultra-fast folding also brought new questions, concerning the intrinsic limits of the folding rates and experimental signatures of barrier-less "downhill" folding. These problems will require novel approaches for even more detailed experimental investigations of the folding dynamics as well as for the analysis of the folding kinetic data. For theoretical simulations of folding, a main challenge is how to extract the relevant information from overwhelmingly detailed atomistic trajectories. New theoretical methods have been devised to allow a systematic approach towards a quantitative analysis of the kinetic network of folding-unfolding transitions between various configuration states of a protein, revealing the transition states and the associated folding pathways at multiple levels, from atomistic to coarse-grained representations. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Balance between alpha and beta structures in ab initio protein folding

Despite initial successes in folding of proteins by molecular simulation, it is becoming increasingly evident that current energy functions (force fields) tend to favor either alpha or beta secondary structure, such that the choice of force field is governed by the structural class of the protein. Here, we study the folding of peptides with either predominantly alpha (Trp cage) or beta (GB1 hairpin) structure with a modified version of the Amber ff03 force field, optimized to reproduce structural propensity in a helix-forming peptide. Using extensive replica exchange molecular dynamics simulations starting from completely unfolded configurations, we obtain the correct folded structure for each peptide, in close agreement with the experimental native structure (<1.5 A all-atom root-mean-square deviation). We obtain converged equilibrium distributions, with folded populations at standard conditions (approximately 300 K), in remarkable accord with experiment. Further comparison to experimental data from NMR spectroscopy and FRET suggests that although the folded structures are accurately reproduced, the unfolded state remains too structured and compact. Our results suggest that the backbone correction results in a force field that is transferable to the folding of proteins from different structural classes.

How the diffusivity profile reduces the arbitrariness of protein folding free energies

The concept of a protein diffusing in its free-energy folding landscape has been fruitful for both theory and experiment. Yet the choice of the reaction coordinate (RC) introduces an undesirable degree of arbitrariness into the problem. We analyze extensive simulation data of an alpha-helix in explicit water solvent as it stochastically folds and unfolds. The free-energy profiles for different RCs exhibit significant variations, some having an activation barrier, while others not. We show that this variation has little effect on the predicted folding kinetics if the diffusivity profiles are properly taken into account. This kinetic quasi-universality is rationalized by an RC rescaling, which, due to the reparameterization invariance of the Fokker-Planck equation, allows the combination of free-energy and diffusivity effects into a single function, the rescaled free-energy profile. This rescaled free energy indeed shows less variation among different RCs than the bare free energy and diffusivity profiles separately do, if we properly distinguish between RCs that contain knowledge of the native state and those that are purely geometric in nature. Our method for extracting diffusivity profiles is easily applied to experimental single molecule time series data and might help to reconcile conflicts that arise when comparing results from different experimental probes for the same protein.

Coordinate and time-dependent diffusion dynamics in protein folding

We developed both analytical and simulation methods to explore the diffusion dynamics in protein folding. We found the diffusion as a quantitative measure of escape from local traps along the protein folding funnel with chosen reaction coordinates has two remarkable effects on kinetics. At a fixed coordinate, local escape time depends on the distribution of barriers around it, therefore the diffusion is often time distributed. On the other hand, the environments (local escape barriers) change along the coordinates, therefore diffusion is coordinate dependent. The effects of time-dependent diffusion on folding can lead to non-exponential kinetics and non-Poisson statistics of folding time distribution. The effects of coordinate dependent diffusion on folding can lead to the change of the kinetic barrier height as well as the position of the corresponding transition state and therefore modify the folding kinetic rates as well as the kinetic routes. Our analytical models for folding are based on a generalized Fokker-Planck diffusion equation with diffusion coefficient both dependent on coordinate and time. Our simulation for folding are based on structure-based folding models with a specific fast folding protein CspTm studied experimentally on diffusion and folding with single molecules. The coordinate and time-dependent diffusion are especially important to be considered in fast folding and single molecule studies, when there is a small or no free energy barrier and kinetics is controlled by diffusion while underlying statistics of kinetics become important. Including the coordinate dependence of diffusion will challenge the transition state theory of protein folding. The classical transition state theory will have to be modified to be consistent. The more detailed folding mechanistic studies involving phi value analysis based on the classical transition state theory will also have to be quantitatively modified. Complex kinetics with multiple time scales may allow us not only to explore the folding kinetics but also probe the local landscape and barrier height distribution with single-molecule experiments.

A natural missing link between activated and downhill protein folding scenarios

We propose protein PTB1 : 4W as a good candidate for engineering into a downhill folder. PTB1 : 4W has a probe-dependent thermal unfolding curve and sub-millisecond T-jump relaxation kinetics on more than one time scale. Its refolding rate in denaturant is a non-linear function of denaturant concentration (curved chevron plot). Yet at high denaturant concentration its unfolding is probe-independent, and the folding kinetics can be fitted to a single exponential decay. The domain appears to fold via a mechanism between downhill folding and activated folding over several small barriers, and when denaturant is added, one of these barriers greatly increases and simplifies the observed folding to apparent two-state kinetics. We predict the simplest free energy function consistent with the thermal denaturation and kinetics experiments by using the singular value Smoluchowski dynamics (SVSD) model. PTB1 : 4W is a natural 'missing link' between downhill and activated folding. We suggest mutations that could move the protein into the downhill folding limit.