A Criterion That Determines Fast Folding of Proteins: A Model Study

We consider the statistical mechanics of a full set of two-dimensional protein-like heteropolymers, whose thermodynamics is characterized by the coil-to-globular ($T_\theta$) and the folding ($T_f$) transition temperatures. For our model, the typical time scale for reaching the unique native conformation is shown to scale as $\tau_f\sim F(M)\exp(\sigma/\sigma_0)$, where $\sigma=1-T_f/T_\theta$, $M$ is the number of residues, and $F(M)$ scales algebraically with $M$. We argue that $T_f$ scales linearly with the inverse of entropy of low energy non-native states, whereas $T_\theta$ is almost independent of it. As $\sigma\rightarrow 0$, non-productive intermediates decrease, and the initial rapid collapse of the protein leads to structures resembling the native state. Based solely on {\it accessible} information, $\sigma$ can be used to predict sequences that fold rapidly.Comment: 10 pages, latex, figures upon reques

Reply to Comment on "Criterion that Determines the Foldability of Proteins"

We point out that the correlation between folding times and $\sigma = (T_{\theta } - T_{f})/T_{\theta }$ in protein-like heteropolymer models where $T_{\theta }$ and $T_{f}$ are the collapse and folding transition temperatures was already established in 1993 before the other presumed equivalent criterion (folding times correlating with $T_{f}$ alone) was suggested. We argue that the folding times for these models show no useful correlation with the energy gap even if restricted to the ensemble of compact structures as suggested by Karplus and Shakhnovich (cond-mat/9606037).Comment: 6 pages, Latex, 2 Postscript figures. Plots explicitly showing the lack of correlation between folding time and energy gap are adde

Twist and writhe dynamics of stiff filaments

This letter considers the dynamics of a stiff filament, in particular the coupling of twist and bend via writhe. The time dependence of the writhe of a filament is $W_r^2\sim L t^{1/4}$ for a linear filament and $W_r^2\sim t^{1/2} / L$ for a curved filament. Simulations are used to study the relative importance of crankshaft motion and tube like motion in twist dynamics. Fuller's theorem, and its relation with the Berry phase, is reconsidered for open filamentsComment: 7 Pages with 2 figure

The effect of local thermal fluctuations on the folding kinetics: a study from the perspective of the nonextensive statistical mechanics

Protein folding is a universal process, very fast and accurate, which works consistently (as it should be) in a wide range of physiological conditions. The present work is based on three premises, namely: ($i$) folding reaction is a process with two consecutive and independent stages, namely the search mechanism and the overall productive stabilization; ($ii$) the folding kinetics results from a mechanism as fast as can be; and ($iii$) at nanoscale dimensions, local thermal fluctuations may have important role on the folding kinetics. Here the first stage of folding process (search mechanism) is focused exclusively. The effects and consequences of local thermal fluctuations on the configurational kinetics, treated here in the context of non extensive statistical mechanics, is analyzed in detail through the dependence of the characteristic time of folding ($\tau$) on the temperature $T$ and on the nonextensive parameter $q$.The model used consists of effective residues forming a chain of 27 beads, which occupy different sites of a $3-$D infinite lattice, representing a single protein chain in solution. The configurational evolution, treated by Monte Carlo simulation, is driven mainly by the change in free energy of transfer between consecutive configurations. ...Comment: 19 pages, 3 figures, 1 tabl

Sequence Dependence of Self-Interacting Random Chains

We study the thermodynamic behavior of the random chain model proposed by Iori, Marinari and Parisi, and how this depends on the actual sequence of interactions along the chain. The properties of randomly chosen sequences are compared to those of designed ones, obtained through a simulated annealing procedure in sequence space. We show that the transition to the folded phase takes place at a smaller strength of the quenched disorder for designed sequences. As a result, folding can be relatively fast for these sequences.Comment: 14 pages, uuencoded compressed postscript fil

Soliton concepts and the protein structure

Structural classification shows that the number of different protein folds is surprisingly small. It also appears that proteins are built in a modular fashion, from a relatively small number of components. Here we propose to identify the modular building blocks of proteins with the dark soliton solution of a generalized discrete nonlinear Schrodinger equation. For this we show that practically all protein loops can be obtained simply by scaling the size and by joining together a number of copies of the soliton, one after another. The soliton has only two loop specific parameters and we identify their possible values in Protein Data Bank. We show that with a collection of 200 sets of parameters, each determining a soliton profile that describes a different short loop, we cover over 90% of all proteins with experimental accuracy. We also present two examples that describe how the loop library can be employed both to model and to analyze the structure of folded proteins.Comment: 7 pages 6 fig

Entropic Barriers, Frustration and Order: Basic Ingredients in Protein Folding

We solve a model that takes into account entropic barriers, frustration, and the organization of a protein-like molecule. For a chain of size $M$, there is an effective folding transition to an ordered structure. Without frustration, this state is reached in a time that scales as $M^{\lambda}$, with $\lambda\simeq 3$. This scaling is limited by the amount of frustration which leads to the dynamical selectivity of proteins: foldable proteins are limited to $\sim 300$ monomers; and they are stable in {\it one} range of temperatures, independent of size and structure. These predictions explain generic properties of {\it in vivo} proteins.Comment: 4 pages, 4 Figures appended as postscript fil

Exploring the Levinthal limit in protein folding

According to the thermodynamic hypothesis, the native state of proteins is uniquely defined by their amino acid sequence. On the other hand, according to Levinthal, the native state is just a local minimum of the free energy and a given amino acid sequence, in the same thermodynamic conditions, can assume many, very different structures that are as thermodynamically stable as the native state. This is the Levinthal limit explored in this work. Using computer simulations, we compare the interactions that stabilize the native state of four different proteins with those that stabilize three non-native states of each protein and find that the nature of the interactions is very similar for all such 16 conformers. Furthermore, an enhancement of the degree of fluctuation of the non-native conformers can be explained by an insufficient relaxation to their local free energy minimum. These results favor Levinthal's hypothesis that protein folding is a kinetic non-equilibrium process.FCT - Foundation for Science and Technology, Portugal [UID/Multi/04326/2013]; Fundacao de Amparo a Pesquisa do Estado de Sao Paulo (FAPESP); Conselho Nacional de Desenvolvimento Cientia co e Tecnologico (CNPq

Role of Secondary Motifs in Fast Folding Polymers: A Dynamical Variational Principle

A fascinating and open question challenging biochemistry, physics and even geometry is the presence of highly regular motifs such as alpha-helices in the folded state of biopolymers and proteins. Stimulating explanations ranging from chemical propensity to simple geometrical reasoning have been invoked to rationalize the existence of such secondary structures. We formulate a dynamical variational principle for selection in conformation space based on the requirement that the backbone of the native state of biologically viable polymers be rapidly accessible from the denatured state. The variational principle is shown to result in the emergence of helical order in compact structures.Comment: 4 pages, RevTex, 4 eps figure