PubMed Health⌕ Search
Home/Archive/PubMed

PubMed · 7710478

How does a protein fold?

A Sali·E Shakhnovich·M Karplus
View original source ↗Export metadata ↓

Abstract

The number of all possible conformations of a polypeptide chain is too large to be sampled exhaustively. Nevertheless, protein sequences do fold into unique native states in seconds (the Levinthal paradox). To determine how the Levinthal paradox is resolved, we use a lattice Monte Carlo model in which the global minimum (native state) is known. The necessary and sufficient condition for folding in this model is that the native state be a pronounced global minimum on the potential surface. This guarantees thermodynamic stability of the native state at a temperature where the chain does not get trapped in local minima. Folding starts by a rapid collapse from a random-coil state to a random semi-compact globule. It then proceeds by a slow, rate-determining search through the semi-compact states to find a transition state from which the chain folds rapidly to the native state. The elements of the folding mechanism that lead to the resolution of the Levinthal paradox are the reduced number of conformations that need to be searched in the semi-compact globule (approximately 10(10) versus approximately 10(16) for the random coil) and the existence of many (approximately 10(3)) transition states. The results have evolutionary implications and suggest principles for the folding of real proteins.

Explore related subjects

Models, ChemicalMonte Carlo MethodProtein FoldingThermodynamicsJournal ArticleResearch Support, Non-U.S. Gov'tResearch Support, U.S. Gov't, Non-P.H.S.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

A Sali, E Shakhnovich, M Karplus. 1994-05-19. How does a protein fold?. https://doi.org/10.1038/369248a0

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

PubMed1999

Folding pathway of a lattice model for proteins.

The folding of a protein-like heteropolymer is studied by using direct simulation of a lattice model that folds rapidly to a well-defined "native" structure. The details of each molecular folding event depend on the random initial conformation as well as the random thermal fluctuations of the polymer. By analyzing the statistical properties of hundreds of folding events, a classical folding "pathway" for such a polymer is found that includes partially folded, on-pathway intermediates that are shown to be metastable equilibrium states of the polymer. These results are discussed in the context of the "classical" and "new" views of folding.

Models, Chemical↗
PubMed1998

Polymorphism and disorder of poly(dA).poly(dT) in fibers.

The absence of proportionality between fiber length and interbase axial distance reveals a high degree of disorder in fibers of poly(dA).poly(dT) when the relative humidity (r.h.) varies. According to X-ray patterns, the increase of fiber length does not correspond to any conformational change. The B' structure is present alone and in a better-organized form at higher r.h. We noted a dependence of length values on the direction of variation of the r.h. An unusual minimum of fiber length is observed around 65% r.h. and explained as been due to an amorphous part of the fiber, essentially composed of single strands of poly(dA) and poly(dT). At a temperature higher than 30 degrees C, the transition B'(alpha)-B* takes place but is associated with an important decrease of fiber length in spite of constancy of the helical rise per base pair. It is shown that this behavior is also due to the amorphous part of the fiber and not to a conformational transition.

Models, Chemical↗
PubMed1998

Conformational interconversions in peptide beta-turns: analysis of turns in proteins and computational estimates of barriers.

The two most important beta-turn features in peptides and proteins are the type I and type II turns, which differ mainly in the orientation of the central peptide unit. Facile conformational interconversion is possible, in principle, by a flip of the central peptide unit. Homologous crystal structures afford an opportunity to structurally characterize both possible conformational states, thus allowing identification of sites that are potentially stereochemically mobile. A representative data set of 250 high-resolution (</=2.0 A), non-homologous protein crystal structures and corresponding variant and homologous entries, obtained from the Brookhaven Protein Data Bank, was examined to identify turns that are assigned different conformational types (type I/type II) in related structures. A total of 55 examples of beta-turns were identified as possible candidates for a stereochemically mobile site. Of the 55 examples, 45 could be classified as a potential site for interconversion between type I and type II beta-turns, while ten correspond to flips from type I' to type II' structures. As a further check, the temperature factors of the central peptide unit carbonyl oxygen atom of the 55 examples were examined. The analysis reveals that the turn assignments are indeed reliable. Examination of the secondary structures at the flanking positions of the flippable beta-turns reveals that seven examples occur in the loop region of beta-hairpins, indicating that the formation of ordered secondary structures on either side of the beta-turn does not preclude local conformational variations. In these beta-turns, Pro (11 examples), Lys (nine examples) and Ser (seven examples) were most often found at the i+1 position. Glycine was found to occur overwhelmingly at position i+2 (28 examples), while Ser (seven examples) and Asn (six examples) were amongst the most frequent residues. Activation energy barriers for the interconversion between type I and type II beta-turns were computed using the peptide models Ac-Pro-Aib-NHMe and Ac-Pro-Gly-NHMe within the framework of the AM1 semi-empirical molecular orbital procedure. In order to have a uniform basis for comparison and to eliminate the distracting influence of the deviation of backbone dihedral angles from that expected for ideal beta-turns, the dihedral angles phii+1 and psii+2 were fixed at the ideal values (phii+1=-60 degrees and psii+2=0 degrees). The other two angles (psii+1 and phii+2) were varied systematically to go from type II to type I beta-turn structures. The computational results suggest that there exists one stereospecific, concerted flip of the central peptide unit involving correlated single bond rotation that can occur with an activation barrier of the order of 3 kcal/mol. The results presented here suggest that conformational variations in beta-turns are observed in protein crystal structures and such changes may be an important dynamic feature in solution.

Models, Chemical↗