Browse
Submit Data
Databases
API
Help

Dataset Information

31 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

How the folding rate constant of simple, single-domain proteins depends on the number of native contacts.

ABSTRACT: Experiments have shown that the folding rate constants of two dozen structurally unrelated, small, single-domain proteins can be expressed in terms of one quantity (the contact order) that depends exclusively on the topology of the folded state. Such dependence is unique in chemical kinetics. Here we investigate its physical origin and derive the approximate formula ln(k) = ln(N) + a + bN, were N is the number of contacts in the folded state, and a and b are constants whose physical meaning is understood. This formula fits well the experimentally determined folding rate constants of the 24 proteins, with single values for a and b.

SUBMITTER: Makarov DE

PROVIDER: S-EPMC122558 | biostudies-literature | 2002 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

How the folding rate constant of simple, single-domain proteins depends on the number of native contacts.

Makarov Dmitrii E DE Keller Craig A CA Plaxco Kevin W KW Metiu Horia H

Proceedings of the National Academy of Sciences of the United States of America 20020301 6

Experiments have shown that the folding rate constants of two dozen structurally unrelated, small, single-domain proteins can be expressed in terms of one quantity (the contact order) that depends exclusively on the topology of the folded state. Such dependence is unique in chemical kinetics. Here we investigate its physical origin and derive the approximate formula ln(k) = ln(N) + a + bN, were N is the number of contacts in the folded state, and a and b are constants whose physical meaning is u ...[more]

PMID: 11904417

Similar Datasets

Native contacts determine protein folding mechanisms in atomistic simulations.

Project description:The recent availability of long equilibrium simulations of protein folding in atomistic detail for more than 10 proteins allows us to identify the key interactions driving folding. We find that the collective fraction of native amino acid contacts, Q, captures remarkably well the transition states for all the proteins with a folding free energy barrier. Going beyond this global picture, we devise two different measures to quantify the importance of individual interresidue contacts in the folding mechanism: (i) the log-ratio of lifetimes of contacts during folding transition paths and in the unfolded state and (ii) a Bayesian measure of how predictive the formation of each contact is for being on a transition path. Both of these measures indicate that native, or near-native, contacts are important for determining mechanism, as might be expected. More remarkably, however, we found that for almost all the proteins, with the designed protein α3D being a notable exception, nonnative contacts play no significant part in determining folding mechanisms.

| S-EPMC3816414 | biostudies-literature

Rate constant and reaction coordinate of Trp-cage folding in explicit water.

Project description:We report rate constant calculations and a reaction coordinate analysis of the rate-limiting folding and unfolding process of the Trp-cage mini-protein in explicit solvent using transition interface sampling. Previous transition path sampling simulations revealed that in this (un)folding process the protein maintains its compact configuration, while a (de)increase of secondary structure is observed. The calculated folding rate agrees reasonably with experiment, while the unfolding rate is 10 times higher. We discuss possible origins for this mismatch. We recomputed the rates with the forward flux sampling method, and found a discrepancy of four orders of magnitude, probably caused by the method's higher sensitivity to the choice of order parameter with respect to transition interface sampling. Finally, we used the previously computed transition path-sampling ensemble to screen combinations of many order parameters for the best model of the reaction coordinate by employing likelihood maximization. We found that a combination of the root mean-square deviation of the helix and of the entire protein was, of the set of tried order parameters, the one that best describes the reaction coordination.

| S-EPMC2567936 | biostudies-literature

Competition between native topology and nonnative interactions in simple and complex folding kinetics of natural and designed proteins.

Project description:We compared folding properties of designed protein Top7 and natural protein S6 by using coarse-grained chain models with a mainly native-centric construct that accounted also for nonnative hydrophobic interactions and desolvation barriers. Top7 and S6 have similar secondary structure elements and are approximately equal in length and hydrophobic composition. Yet their experimental folding kinetics were drastically different. Consistent with experiment, our simulated folding chevron arm for Top7 exhibited a severe rollover, whereas that for S6 was essentially linear, and Top7 model kinetic relaxation was multiphasic under strongly folding conditions. The peculiar behavior of Top7 was associated with several classes of kinetic traps in our model. Significantly, the amino acid residues participating in nonnative interactions in trapped conformations in our Top7 model overlapped with those deduced experimentally. These affirmations suggest that the simple ingredients of native topology plus sequence-dependent nonnative interactions are sufficient to account for some key features of protein folding kinetics. Notably, when nonnative interactions were absent in the model, Top7 chevron rollover was not correctly predicted. In contrast, nonnative interactions had little effect on the quasi linearity of the model folding chevron arm for S6. This intriguing distinction indicates that folding cooperativity is governed by a subtle interplay between the sequence-dependent driving forces for native topology and the locations of favorable nonnative interactions entailed by the same sequence. Constructed with a capability to mimic this interplay, our simple modeling approach should be useful in general for assessing a designed sequence's potential to fold cooperatively.

| S-EPMC2840274 | biostudies-literature

Coevolution of function and the folding landscape: correlation with density of native contacts.

Project description:The relationship between the folding landscape and function of evolved proteins is explored by comparison of the folding mechanisms for members of the flavodoxin fold. CheY, Spo0F, and NtrC have unrelated functions and low sequence homology but share an identical topology. Recent coarse-grained simulations show that their folding landscapes are uniquely tuned to properly suit their respective biological functions. Enhanced packing in Spo0F and its limited conformational dynamics compared to CheY or NtrC lead to frustration in its folding landscape. Simulation as well as experimental results correlate with the local density of native contacts for these and a sample of other proteins. In particular, protein regions of low contact density are observed to become structured late in folding; concomitantly, these dynamic regions are often involved in binding or conformational rearrangements of functional importance. These observations help to explain the widespread success of Gō-like coarse-grained models in reproducing protein dynamics.

| S-EPMC2567933 | biostudies-literature

Folding and Stabilization of Native-Sequence-Reversed Proteins.

Project description:Though the problem of sequence-reversed protein folding is largely unexplored, one might speculate that reversed native protein sequences should be significantly more foldable than purely random heteropolymer sequences. In this article, we investigate how the reverse-sequences of native proteins might fold by examining a series of small proteins of increasing structural complexity (?-helix, ?-hairpin, ?-helix bundle, and ?/?-protein). Employing a tandem protein structure prediction algorithmic and molecular dynamics simulation approach, we find that the ability of reverse sequences to adopt native-like folds is strongly influenced by protein size and the flexibility of the native hydrophobic core. For ?-hairpins with reverse-sequences that fail to fold, we employ a simple mutational strategy for guiding stable hairpin formation that involves the insertion of amino acids into the ?-turn region. This systematic look at reverse sequence duality sheds new light on the problem of protein sequence-structure mapping and may serve to inspire new protein design and protein structure prediction protocols.

| S-EPMC4844985 | biostudies-literature

A protein folding potential that places the native states of a large number of proteins near a local minimum.

Project description:BackgroundWe present a simple method to train a potential function for the protein folding problem which, even though trained using a small number of proteins, is able to place a significantly large number of native conformations near a local minimum. The training relies on generating decoys by energy minimization of the native conformations using the current potential and using a physically meaningful objective function (derivative of energy with respect to torsion angles at the native conformation) during the quadratic programming to place the native conformation near a local minimum.ResultsWe also compare the performance of three different types of energy functions and find that while the pairwise energy function is trainable, a solvation energy function by itself is untrainable if decoys are generated by minimizing the current potential starting at the native conformation. The best results are obtained when a pairwise interaction energy function is used with solvation energy function.ConclusionsWe are able to train a potential function using six proteins which places a total of 42 native conformations within approximately 4 A rmsd and 71 native conformations within approximately 6 A rmsd of a local minimum out of a total of 91 proteins. Furthermore, the threading test using the same 91 proteins ranks 89 native conformations to be first and the other two as second.

| S-EPMC126205 | biostudies-literature

Native states of fast-folding proteins are kinetic traps.

Project description:It has been suggested that the native state of a protein acts as a kinetic hub that can facilitate transitions between nonnative states. Using recently developed tools to quantify mediation probabilities ("hub scores"), we quantify hub-like behavior in atomic resolution trajectories for the first time. We use a data set of trajectory ensembles for 12 fast-folding proteins previously published by D. E. Shaw Research (Lindorff-Larsen, K.; et al. How Fast-Folding Proteins Fold. Science2011, 334, 517) with an aggregate simulation time of over 8.2 ms. We visualize the free-energy landscape of each molecule using configuration space networks, and show that dynamic quantities can be qualitatively understood from visual inspection of the networks. Modularity optimization is used to provide a parameter-free means of tessellating the network into a group of communities. Using hub scores, we find that the percentage of trajectories that are mediated by the native state is 31% when averaged over all molecules, and reaches a maximum of 52% for the homeodomain and chignolin. Furthermore, for these mediated transitions, we use Markov models to determine whether the native state acts as a facilitator for the transition, or as a trap (i.e., an off-pathway detour). Although instances of facilitation are found in 4 of the 12 molecules, we conclude that the native state acts primarily as a trap, which is consistent with the idea of a funnel-like landscape.

| S-EPMC3619186 | biostudies-literature

Dynamics of equilibrium folding and unfolding transitions of titin immunoglobulin domain under constant forces.

Project description:The mechanical stability of force-bearing proteins is crucial for their functions. However, slow transition rates of complex protein domains have made it challenging to investigate their equilibrium force-dependent structural transitions. Using ultra stable magnetic tweezers, we report the first equilibrium single-molecule force manipulation study of the classic titin I27 immunoglobulin domain. We found that individual I27 in a tandem repeat unfold/fold independently. We obtained the force-dependent free energy difference between unfolded and folded I27 and determined the critical force (?5.4 pN) at which unfolding and folding have equal probability. We also determined the force-dependent free energy landscape of unfolding/folding transitions based on measurement of the free energy cost of unfolding. In addition to providing insights into the force-dependent structural transitions of titin I27, our results suggest that the conformations of titin immunoglobulin domains can be significantly altered during low force, long duration muscle stretching.

| S-EPMC4446711 | biostudies-literature

Folding pathway of a multidomain protein depends on its topology of domain connectivity.

Project description:How do the folding mechanisms of multidomain proteins depend on protein topology? We addressed this question by developing an Ising-like structure-based model and applying it for the analysis of free-energy landscapes and folding kinetics of an example protein, Escherichia coli dihydrofolate reductase (DHFR). DHFR has two domains, one comprising discontinuous N- and C-terminal parts and the other comprising a continuous middle part of the chain. The simulated folding pathway of DHFR is a sequential process during which the continuous domain folds first, followed by the discontinuous domain, thereby avoiding the rapid decrease in conformation entropy caused by the association of the N- and C-terminal parts during the early phase of folding. Our simulated results consistently explain the observed experimental data on folding kinetics and predict an off-pathway structural fluctuation at equilibrium. For a circular permutant for which the topological complexity of wild-type DHFR is resolved, the balance between energy and entropy is modulated, resulting in the coexistence of the two folding pathways. This coexistence of pathways should account for the experimentally observed complex folding behavior of the circular permutant.

| S-EPMC4234588 | biostudies-literature

Simple MD-based model for oxidative folding of peptides and proteins.

Project description:Significant strides have been recently made to fold peptides and small proteins in silico using MD simulations. However, facilities are currently lacking to include disulfide bonding in the MD models of protein folding. To address this problem, we have developed a simple empirical protocol to model formation of disulfides, which is perturbation-free, retains the same speed as conventional MD simulations and allows one to control the reaction rate. The new protocol has been tested on 15-aminoacid peptide guanylin containing four cysteine residues; the net simulation time using Amber ff14SB force field was 61??s. The resulting isomer distribution is in qualitative agreement with experiment, suggesting that oxidative folding of guanylin in vitro occurs under kinetic control. The highly stable conformation of the so-called isomer 2(B) has been obtained for full-length guanylin, which is significantly different from the poorly ordered structure of the truncated peptide PDB ID 1GNB. In addition, we have simulated oxidative folding of guanylin within the 94-aminoacid prohormone proguanylin. The obtained structure is in good agreement with the NMR coordinates 1O8R. The proposed modeling strategy can help to explore certain fundamental aspects of protein folding and is potentially relevant for manufacturing of synthetic peptides and recombinant proteins.

| S-EPMC5570944 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data