Microscopic theory of critical folding nuclei and reconfiguration activation barriers in folding proteins

Comparing a simple theoretical model for protein folding with all-atom molecular dynamics simulations

Advances in computing have enabled microsecond all-atom molecular dynamics trajectories of protein folding that can be used to compare with and test critical assumptions of theoretical models. We show that recent simulations by the Shaw group (10, 11, 14, 15) are consistent with a key assumption of an Ising-like theoretical model that native structure grows in only a few regions of the amino acid sequence as folding progresses. The distribution of mechanisms predicted by simulating the master equation of this native-centric model for the benchmark villin subdomain, with only two adjustable thermodynamic parameters and one temperature-dependent kinetic parameter, is remarkably similar to the distribution in the molecular dynamics trajectories.

Variationally determined free energy profiles for structural models of proteins: characteristic temperatures for folding and trapping

Characterizing the phase diagram for proteins is important both for laboratory studies and for the development of structure prediction algorithms. Using a variational scheme, we calculated the generic features of the protein thermostability over a large range of temperatures for a set of more than 50 different proteins using a model based on native structure alone. Focusing on a specific system, protein G, we further examined, using a more realistic model that includes the nonnative interaction, the thermostability of both the native state and a collection of trap structures. By surveying the native structures for many proteins and by paying closer attention to the various trap structures of protein G, we obtained an overall understanding of the folding dynamics far from the conditions usually focused on; namely, those near the folding temperature alone. Two characteristic temperatures (shown to scale with folding temperature in general) signal drastic changes in the folding mechanism. The variational calculations suggest that most proteins would, indeed, fold in a barrierless manner below a critical temperature analogous to a spinodal in crystallization. For fixed interaction strengths, this temperature, however, seems to be generally very low, approximately 50% of the equilibrium folding temperature. Likewise, native proteins, in general, would unfold in a completely barrierless way at a temperature 25% above folding temperature according to these variational calculations. We also studied the distribution of free energy profiles for escape from a set of trap structures generated by simulations.

Origins of barriers and barrierless folding in BBL

In contrast to classical chemical phenomenology, theory suggests that proteins may undergo downhill folding without an activation barrier under certain thermodynamic conditions. Recently, the BBL protein was proposed to fold by such a downhill scenario, but discrepancies between experimental results found in different groups argue against this. After briefly reviewing the major experimental studies of the BBL folding mechanism, we show that simulations of both coarse-grained and atomistic models can reconcile the seemingly conflicting observations.

Thermodynamics of Gō-type models for protein folding

Go-type potentials, based on the inter-residue contacts present in the native structure of a protein, are frequently used to predict dynamic and structural features of the folding pathways through computer simulations. However, the mathematical form used to define the model interactions includes several arbitrary choices, whose consequences are not usually analyzed. In this work, we use a simple off-lattice protein model and a parallel tempering Monte Carlo simulation technique to carry out such analysis, centered in the thermodynamic characteristics of the folding transition. We show how the definition of a native contact has a deep impact on the presence of simple or complex transitions, with or without thermodynamic intermediates. In addition, we have checked that the width of the attractive wells has a profound effect on the free-energy barrier between the folded and unfolded states, mainly through its influence on the entropy of the denatured state.

Localization and freezing of a Gaussian chain in a quenched random potential

The Gaussian chain in a quenched random potential (which is characterized by the disorder strength Delta) is investigated in the d-dimensional space by the replicated variational method. The general expression for the free energy within so-called one-step-replica symmetry breaking (1-RSB) scenario has been systematically derived. We have shown that the replica symmetrical (RS) limit of this expression can describe the chain center-of-mass localization and collapse. The critical disorder when the chain becomes localized scales as Delta(c) approximately b(d)N(-2+d/2) (where b is the length of the Kuhn segment length and N is the chain length) whereas the chain gyration radius R(g) approximately b(b(d)/Delta)(1/(4-d)). The freezing of the internal degrees of freedom follows to the 1-RSB-scenario and is characterized by the beads localization length D(2). It was demonstrated that the solution for D(2) appears as a metastable state at Delta=Delta(A) and behaves similarly to the corresponding frozen states in heteropolymers or in p-spin random spherical model.

Embedding a native state into a random heteropolymer model: the dynamic approach

We study a random heteropolymer model with Langevin dynamics, in the supersymmetric formulation. Employing a procedure similar to one that has been used in static calculations, we construct an ensemble in which the affinity of the system for a native state is controlled by a "selection temperature" T0. In the limit of high T0, the model reduces to a random heteropolymer, while for T0-->0 the system is forced into the native state. Within the Gaussian variational approach that we employed previously for the random heteropolymer, we explore the phases of the system for high and low T0. For high T0, the system exhibits a (dynamical) spin-glass phase, like that found for the random heteropolymer, below a temperature T(g). For low T0, we find an ordered phase, characterized by a nonzero overlap with the native state, below a temperature T(n) proportional to 1/T(0)>T(g). However, the random-globule phase remains locally stable below T(n), down to the dynamical glass transition at T(g). Thus, in this model, folding is rapid for temperatures between T(g) and T(n), but below T(g) the system can get trapped in conformations uncorrelated with the native state. At a lower temperature, the ordered phase can also undergo a dynamical glass transition, splitting into substates separated by large barriers.

Universality classes in folding times of proteins

Molecular dynamics simulations in simplified models allow one to study the scaling properties of folding times for many proteins together under a controlled setting. We consider three variants of the Go models with different contact potentials and demonstrate scaling described by power laws and no correlation with the relative contact order parameter. We demonstrate existence of at least three kinetic universality classes that are correlated with the types of structure: the alpha-, alpha-beta-, and beta- proteins have the scaling exponents of approximately 1.7, 2.5, and 3.2, respectively. The three classes merge into one when the contact range is truncated at a reasonable value. We elucidate the role of the potential associated with the chirality of a protein.

Self-generated disorder and structural glass formation in homopolymer globules

We have investigated the interrelation between spin glasses and structural glasses. Spin glasses in this case are p-spin interaction spin glasses (at p>2) or Potts glasses that contain quenched disorder, whereas the structural glasses are here exemplified by a homopolymeric globule, which can be viewed as a liquid of connected units on a nanoscale. It is argued that the homopolymeric globule problem can be mapped onto a disorder field theoretical model whose effective Hamiltonian resembles the corresponding one for the spin glass model. In this sense the disorder in the globule is self-generated (in contrast to spin glasses) and can be related to competing interactions (virial coefficients of different signs) and the chain connectivity. The work is aimed at giving a quantitative description of this analogy. We have investigated in the mean-field approximation the phase diagram of the homopolymeric globule where the transition line from the liquid to glassy globule is treated in terms of the replica symmetry breaking paradigm. The configurational entropy temperature dependence is also discussed.

Kinetic nonoptimality and vibrational stability of proteins

Scaling of folding times in Go models of proteins and of decoy structures with the Lennard-Jones potentials in the native contacts reveal power law trends when studied under optimal folding conditions. The power law exponent depends on the type of native geometry. Its value indicates lack of kinetic optimality in the model proteins. In proteins, mechanical and thermodynamic stabilities are correlated.

How does a beta -hairpin fold/unfold? competition between topology and heterogeneity in a solvable model

We study the competition between topological effects and sequence inhomogeneities in determining the thermodynamics and the un/folding kinetics of a beta-hairpin. Our work utilizes a new exactly solvable model that allows for arbitrary configurations of native contacts. In general, the competition between heterogeneity and topology results in a crossover of the dominant transition state. Interestingly, near this crossover, the single reaction coordinate picture can be seriously misleading. Our results also suggest that inferring the folding pathway from unfolding simulations is not always justified.

Two state behavior in a solvable model of beta-hairpin folding

Understanding the mechanism of protein secondary structure formation is an essential part of the protein-folding puzzle. Here we describe a simple model for the formation of a beta hairpin, motivated by the fact that folding of a beta hairpin captures much of the basic physics of protein folding. The modeled hairpin is composed of two interacting Gaussian chains with one pairwise (two-body) and two many-body interactions. We show that these many-body interactions, arising from side chain packing effects, are responsible for producing an "all-or-none" folding transition. We also estimate the (single exponential) folding/unfolding rate via calculating the thermodynamic weight of the "critical" droplet/bubble.

Limited internal friction in the rate-limiting step of a two-state protein folding reaction

Small, single-domain proteins typically fold via a compact transition-state ensemble in a process well fitted by a simple, two-state model. To characterize the rate-limiting conformational changes that underlie two-state folding, we have investigated experimentally the effects of changing solvent viscosity on the refolding of the IgG binding domain of protein L. In conjunction with numerical simulations, our results indicate that the rate-limiting conformational changes of the folding of this domain are strongly coupled to solvent viscosity and lack any significant "internal friction" arising from intrachain collisions. When compared with the previously determined solvent viscosity dependencies of other, more restricted conformational changes, our results suggest that the rate-limiting folding transition involves conformational fluctuations that displace considerable amounts of solvent. Reconciling evidence that the folding transition state ensemble is comprised of highly collapsed species with these and similar, previously reported results should provide a significant constraint for theoretical models of the folding process.