Protein folding in high-dimensional spaces:hypergutters and the role of non-native interactions

The prediction of the three-dimensional native structure of proteins from the knowledge of their amino acid sequence, known as the protein folding problem, is one of the most important yet unsolved issues of modern science. Since the conformational behaviour of flexible molecules is nothing more than a complex physical problem, increasingly more physicists are moving into the study of protein systems, bringing with them powerful mathematical and computational tools, as well a...

The energy landscapes of proteins have evolved to be different from most random heteropolymers. Many studies have concluded that evolutionary selection for rapid and reliable folding to a given structure that is stable at biological temperatures leads to energy landscapes having a single dominant basin and an overall funnel topography. We show here that, while such a landscape topography is indeed a sufficient condition for folding, another possibility also exists, giving a n...

Perturbing a Go model towards a realistic protein Hamiltonian by adding non-native interactions, we find that the folding rate is in general enhanced as ruggedness is initially increased, as long as the protein is sufficiently large and flexible. Eventually the rate drops rapidly towards zero when ruggedness significantly slows conformational transitions. Energy landscape arguments for thermodynamics and kinetics are coupled with a treatment of non-native collapse to elucidat...

A central goal of protein-folding theory is to predict the stochastic dynamics of transition paths --- the rare trajectories that transit between the folded and unfolded ensembles --- using only thermodynamic information, such as a low-dimensional equilibrium free-energy landscape. However, commonly used one-dimensional landscapes typically fall short of this aim, because an empirical coordinate-dependent diffusion coefficient has to be fit to transition-path trajectory data ...

We use a three dimensional cubic lattice model of proteins to study their properties that determine folding to the native state. The protein chain is modeled as a sequence of $N$ beads. The interactions between beads are taken from a Gaussian distribution of energies. We studied 56 sequences with unique ground states for $N = 15$ and $27$. Thermodynamic and kinetic properties were determined using Monte Carlo simulations and exhaustive enumeration. For all sequences we find c...

We present a simple model of protein folding dynamics that captures key qualitative elements recently seen in all-atom simulations. The goals of this theory are to serve as a simple formalism for gaining deeper insight into the physical properties seen in detailed simulations as well as to serve as a model to easily compare why these simulations suggest a different kinetic mechanism than previous simple models. Specifically, we find that non-native contacts play a key role in...

To what extent do general features of folding/unfolding kinetics of small globular proteins follow from their thermodynamic properties? To address this question, we investigate a new simplifed protein chain model that embodies a cooperative interplay between local conformational preferences and hydrophobic burial. The present four-helix-bundle 55mer model exhibits proteinlike calorimetric two-state cooperativity. It rationalizes native-state hydrogen exchange observations. Ou...

By observing trends in the folding kinetics of experimental 2-state proteins at their transition midpoints, and by observing trends in the barrier heights of numerous simulations of coarse grained, C-alpha model, Go proteins, we show that folding rates correlate with the degree of heterogeneity in the formation of native contacts. Statistically significant correlations are observed between folding rates and measures of heterogeneity inherent in the native topology, as well as...

With the help of a simple 20 letters, lattice model of heteropolymers, we investigate the energy landscape in the space of designed good-folder sequences. Low-energy sequences form clusters, interconnected via neutral networks, in the space of sequences. Residues which play a key role in the foldability of the chain and in the stability of the native state are highly conserved, even among the chains belonging to different clusters. If, according to the interaction matrix, som...

The number of protein structures is far less than the number of sequences. By imposing simple generic features of proteins (low energy and compaction) on all possible sequences we show that the structure space is sparse compared to the sequence space. Even though the sequence space grows exponentially with N (the number of amino acids) we conjecture that the number of low energy compact structures only scales as ln N. This implies that many sequences must map onto countable n...