蛋白质折叠机理的理论研究

简要综述了近年来蛋白质折叠机理的理论研究。首先回顾了蛋白质折叠理论的发展历程,然后对折叠中间体的研究现状作了较详细的介绍。同时,对折叠机理理论研究中的几种理论模型和模拟算法作了细致评述,分析了其现状和存在的问题。最后,总结和讨论了折叠机理理论研究的现存问题及研究热点,并展望了该领域研究的发展趋势。     

Gatekeepers in the ribosomal protein s6: thermodynamics, kinetics, and folding pathways revealed by a minimalist protein model

We investigate the effect of structural gatekeepers on the folding of the ribosomal protein S6. Folding thermodynamics and early refolding kinetics are studied for this system utilizing computer simulations of a minimalist protein model. When gatekeepers are eliminated, the thermodynamic signature of a folding intermediate emerges, and a marked decrease in folding efficiency is observed. We explain the prerequisites that determine the "strength" of a given gatekeeper. The investigated gatekeepers are found to have distinct functions, and to guide the folding and time-dependent packing of non-overlapping secondary structure elements in the protein. Gatekeepers avoid kinetic traps during folding by favoring the formation of "productive topologies" on the way to the native state. The trends in folding rates in the presence/absence of gatekeepers observed for our minimalist model of S6 are in very good agreement with experimental data on this protein.     

Hierarchical protein folding pathways: a computational study of protein fragments

We have previously presented a building block folding model. The model postulates that protein folding is a hierarchical top-down process. The basic unit from which a fold is constructed, referred to as a hydrophobic folding unit, is the outcome of combinatorial assembly of a set of "building blocks." Results obtained by the computational cutting procedure yield fragments that are in agreement with those obtained experimentally by limited proteolysis. Here we show that as expected, proteins from the same family give very similar building blocks. However, different proteins can also give building blocks that are similar in structure. In such cases the building blocks differ in sequence, stability, contacts with other building blocks, and in their 3D locations in the protein structure. This result, which we have repeatedly observed in many cases, leads us to conclude that while a building block is influenced by its environment, nevertheless, it can be viewed as a stand-alone unit. For small-sized building blocks existing in multiple conformations, interactions with sister building blocks in the protein will increase the population time of the native conformer. With this conclusion in hand, it is possible to develop an algorithm that predicts the building block assignment of a protein sequence whose structure is unknown. Toward this goal, we have created sequentially nonredundant databases of building block sequences. A protein sequence can be aligned against these, in order to be matched to a set of potential building blocks.     

Rapid folding of the prion protein captured by pressure-jump

The conversion of the cellular form of the prion protein (PrP^C) to an altered disease state, generally denoted as scrapie isoform (PrP^Sc), appears to be a crucial molecular event in prion diseases. The details of this conformational transition are not fully understood, but it is perceived that they are associated with misfolding of PrP or its incapacity to maintain the native fold during its cell cycle. Here we present a tryptophan mutant of PrP (F198W), which has enhanced fluorescence sensitivity to unfolding/refolding transitions. Equilibrium folding was studied by circular dichroism and fluorescence. Pressure-jump experiments were successfully applied to reveal rapid submillisecond folding events of PrP at temperatures not accessed before. D. C. Jenkins and D. S. Pearson contributed equally.     

Statistical analyses of protein folding rates from the view of quantum transition

Understanding protein folding rate is the primary key to unlock the fundamental physics underlying protein structure and its folding mechanism.Especially,the temperature dependence of the folding rate remains unsolved in the literature.Starting from the assumption that protein folding is an event of quantum transition between molecular conformations,we calculated the folding rate for all two-state proteins in a database and studied their temperature dependencies.The non-Arrhenius temperature relation for 16 proteins,whose experimental data had previously been available,was successfully interpreted by comparing the Arrhenius plot with the first-principle calculation.A statistical formula for the prediction of two-state protein folding rate was proposed based on quantum folding theory.The statistical comparisons of the folding rates for 65 two-state proteins were carried out,and the theoretical vs.experimental correlation coefficient was 0.73.Moreover,the maximum and the minimum folding rates given by the theory were consistent with the experimental results.     

Entropic criteria for protein folding derived from recurrences: six residues patch as the basic protein word

Some research has suggested that patches of six constitute an important amino acid window length in proteins for conveying information. We present database evidence that supports this conjecture, as well as additional recurrence-based data that characterization and quantification of these words affect the folding/aggregation features of proteins. Other indirect evidence is presented and discussed.     

Mechanisms of protein folding

Understanding the mechanism by which a polypeptide chain folds into its native structure is a central problem of modern biophysics. The collaborative efforts of experimental and theoretical studies recently raised the tantalizing possibility to define a unifying mechanism for protein folding. In this review we summarize some of these intriguing advances and analyze them together with a discussion on the new findings concerning the so-called downhill folding.     

Membrane protein folding: beyond the two stage model

The folding of alpha-helical membrane proteins has previously been described using the two stage model, in which the membrane insertion of independently stable alpha-helices is followed by their mutual interactions within the membrane to give higher order folding and oligomerization. Given recent advances in our understanding of membrane protein structure it has become apparent that in some cases the model may not fully represent the folding process. Here we present a three stage model which gives considerations to ligand binding, folding of extramembranous loops, insertion of peripheral domains and the formation of quaternary structure.     

Structural features of protein folding nuclei

A crucial event of protein folding is the formation of a folding nucleus. We demonstrate the presence of a considerable coincidence between the location of folding nuclei and the location of so-called "root structural motifs", which have unique overall folds and handedness. In the case of proteins with a single root structural motif, the involvement in the formation of a folding nucleus is in average significantly higher for amino acids residues that are in root structural motifs, compared to residues in other parts of the protein. The tests carried out revealed that the observed difference is statistically reliable. Thus, a structural feature that corresponds to the protein folding nucleus is now found.     

Probing the energy landscape of protein folding/unfolding transition states

Previous molecular dynamics (MD) simulations of the thermal denaturation of chymotrypsin inhibitor 2 (CI2) have provided atomic-resolution models of the transition state ensemble that is well supported by experimental studies. Here, we use simulations to further investigate the energy landscape around the transition state region. Nine structures within approximately 35 ps and 3 A C(alpha) RMSD of the transition state ensemble identified in a previous 498 K thermal denaturation simulation were quenched under the quasi-native conditions of 335 K and neutral pH. All of the structures underwent hydrophobically driven collapse in response to the drop in temperature. Structures less denatured than the transition state became structurally more native-like, while structures that were more denatured than the transition state tended to show additional loss of native structure. The structures in the immediate region of the transition state fluctuated between becoming more and less native-like. All of the starting structures had the same native-like topology and were quite similar (within 3.5 A C(alpha) RMSD). That the structures all shared native-like topology, yet diverged into either more or less native-like structures depending on which side of the transition state they occupied on the unfolding trajectory, indicates that topology alone does not dictate protein folding. Instead, our results suggest that a detailed interplay of packing interactions and interactions with water determine whether a partially denatured protein will become more native-like under refolding conditions.     

Evolutionary computer programming of protein folding and structure predictions

In order to understand the mechanism of protein folding and to assist the rational de-novo design of fast-folding, non-aggregating and stable artificial enzymes it is very helpful to be able to simulate protein folding reactions and to predict the structures of proteins and other biomacromolecules. Here, we use a method of computer programming called "evolutionary computer programming" in which a program evolves depending on the evolutionary pressure exerted on the program. In the case of the presented application of this method on a computer program for folding simulations, the evolutionary pressure exerted was towards faster finding deep minima in the energy landscape of protein folding. Already after 20 evolution steps, the evolved program was able to find deep minima in the energy landscape more than 10 times faster than the original program prior to the evolution process.     

Balancing oxidative protein folding: The influences of reducing pathways on disulfide bond formation

Oxidative protein folding is confined to few compartments, including the endoplasmic reticulum, the mitochondrial intermembrane space and the bacterial periplasm. Conversely, in compartments in which proteins are translated such as the cytosol, the mitochondrial matrix and the chloroplast stroma proteins are kept reduced by the thioredoxin and glutaredoxin systems that functionally overlap. The highly reducing NADPH pool thereby serves as electron donor that enables glutathione reductase and thioredoxin reductase to keep glutathione pools and thioredoxins in their reduced redox state, respectively. Notably, also compartments containing oxidizing machineries are linked to these reducing pathways. Reducing pathways aid in proofreading of disulfide bond formation by isomerization or they provide reducing equivalents for the reduction of disulfides prior to degradation. In addition, they contribute to the thiol-dependent regulation of protein activities, and they help to counteract oxidative stress. The existence of oxidizing and reducing pathways in the same compartment poses a potential problem as the cell has to avoid futile cycles of oxidation and subsequent reduction reactions. Thus, compartments that contain oxidizing machineries have developed sophisticated ways to spatiotemporally balance and regulate oxidation and reduction. In this review, we discuss oxidizing and reducing pathways in the endoplasmic reticulum, the periplasm and the mitochondrial intermembrane space and highlight the role of glutathione especially in the endoplasmic reticulum and the intermembrane space. This article is part of a Special Issue entitled: Thiol-Based Redox Processes.     

Folding barrier in horse cytochrome c: support for a classical folding pathway

Native-state structures and conformations of ferrocytochrome c, nitrosylcytochrome c, and carbonmonoxycytochrome c are very similar. They are, however, immensely different from each other in terms of thermodynamic stability. The dramatic destabilization of ferrocytochrome c to the extent of 12 kcal mol(-1) produces no effect on the folding rate, and this is so in spite of the fact that all three test-tube variants fold in an apparent two-state manner. For all three proteins the folding barrier is early in time, sizable in energy, and is of the same magnitude (approximately 6.5 kcal mol(-1)). These results raise some challenges to the "new view" of protein folding. An early transition state, the search for which consumes most of the observed folding time, is suggested.     

Thermal folding and mechanical unfolding pathways of protein secondary structures

Mechanical stretching of secondary structures is studied through molecular dynamics simulations of a Go-like model. Force versus displacement curves are studied as a function of the stiffness and velocity of the pulling device. The succession of stretching events, as measured by the order in which contacts are ruptured, is compared to the sequencing of events during thermal folding and unfolding. Opposite cross-correlations are found for an alpha-helix and a beta-hairpin structure. In a tandem of two alpha-helices, the two constituent helices unravel nearly simultaneously. A simple condition for simultaneous versus sequential unraveling of repeat units is presented.     

The folding mechanism of a two-domain protein: folding kinetics and domain docking of the gene-3 protein of phage fd

The gene-3 protein (G3P) of filamentous phages is essential for the infection of Escherichia coli. The carboxy-terminal domain anchors this protein in the phage coat, whereas the two amino-terminal domains N1 and N2 protrude from the phage surface. We analyzed the folding mechanism of the two-domain fragment N1-N2 of G3P (G3P(*)) and the interplay between folding and domain assembly. For this analysis, a variant of G3P(*) was used that contained four stabilizing mutations (IIHY-G3P(*)). The observed refolding kinetics extend from 10 ms to several hours. Domain N1 refolds very rapidly (with a time constant of 9.4 ms at 0.5 M guanidinium chloride, 25 degrees C) both as a part of IIHY-G3P(*) and as an isolated protein fragment. The refolding of domain N2 is slower and involves two reactions with time constants of seven seconds and 42 seconds. These folding reactions of the individual domains are followed by a very slow, spectroscopically silent docking process, which shows a time constant of 6200 seconds. This reaction was detected by a kinetic unfolding assay for native molecules. Before docking, N1 and N2 unfold fast and independently, after docking they unfold slowly in a correlated fashion. A high energy barrier is thus created by domain docking, which protects G3P kinetically against unfolding. The slow domain docking is possibly important for the infection of E.coli by the phage. Upon binding to the F pilus, the N2 domain separates from N1 and the binding site for TolA on domain N1 is exposed. Since domain reassembly is so slow, this binding site remains accessible until pilus retraction has brought N1 close to TolA on the bacterial surface.     

Alpha helical crossovers favor right‐handed supersecondary structures by kinetic trapping: The phone cord effect in protein folding

The remarkable predominance of right‐handedness in beta‐alpha‐beta helical crossovers has been previously explained in terms of thermodynamic stability and kinetic accessibility, but a different kinetic trapping mechanism may also play a role. If the beta‐sheet contacts are made before the crossover helix is fully formed, and if the backbone angles of the folding helix follows the energetic pathway of least resistance, then the helix would impart a torque on the ends of the two strands. Such a torque would tear apart a left‐handed conformation but hold together a right‐handed one. Right‐handed helical crossovers predominate even in all‐alpha proteins, where previous explanations based on the preferred twist of the beta sheet do not apply. Using simple molecular simulations, we can reproduce the right‐handed preference in beta‐alpha‐beta units, without imposing specific beta‐strand geometry. The new kinetic trapping mechanism is dubbed the “phone cord effect” because it is reminiscent of the way a helical phone cord forms superhelices to relieve torsional stress. Kinetic trapping explains the presence of a right‐handed superhelical preference in alpha helical crossovers and provides a possible folding mechanism for knotted proteins.     

The effects of nonnative interactions on protein folding rates: theory and simulation

Proteins are minimally frustrated polymers. However, for realistic protein models, nonnative interactions must be taken into account. In this paper, we analyze the effect of nonnative interactions on the folding rate and on the folding free energy barrier. We present an analytic theory to account for the modification on the free energy landscape upon introduction of nonnative contacts, added as a perturbation to the strong native interactions driving folding. Our theory predicts a rate-enhancement regime at fixed temperature, under the introduction of weak, nonnative interactions. We have thoroughly tested this theoretical prediction with simulations of a coarse-grained protein model, by using an off-lattice C(alpha)model of the src-SH3 domain. The strong agreement between results from simulations and theory confirm the nontrivial result that a relatively small amount of nonnative interaction energy can actually assist the folding to the native structure.     

Multiple folding mechanisms of protein ubiquitin

Based on the C(alpha) Go-type model, the folding kinetics and mechanisms of protein ubiquitin with mixed alpha/beta topology are studied by molecular dynamics simulations. The relaxation kinetics shows that there are three phases, namely the major phase, the intermediate phase and the slowest minor phase. The existence of these three phases are relevant to the phenomenon found in experiments. According to our simulations, the folding at high temperatures around the folding transition temperature T(f) is of a two-state process, and the folding nucleus is consisted of contacts between the front end of alpha-helix and the turn(4). The folding at low temperature (approximately T = 0.8) is also studied, where an A-state like structure is found lying on the major folding pathway. The appearance of this structure is related to the stability of the first part (residue 1-51) of protein ubiquitin. As the temperature decreases, the formation of secondary structures, tertiary structures and collapse of the protein are found to be decoupled gradually and the folding mechanism changes from the nucleation-condensation to the diffusion-collision. This feature indicates a unifying common folding mechanism for proteins. The intermediate phase is also studied and is found to represent a folding process via a long-lived intermediate state which is stabilized by strong interactions between the beta(1) and the beta(5) strand. These strong interactions are important for the function of protein ubiquitin as a molecular chaperone. Thus the intermediate phase is assumed as a byproduct of the requirement of protein function. In addition, the validity of the current Go-model is also investigated, and a lower limited temperature for protein ubiquitin T(limit) = 0.8 is proposed. At temperatures higher than this value, the kinetic traps due to glass dynamics cannot be significantly populated and the intermediate states can be reliably identified although there is slight chevron rollover in the folding rates. At temperature lower than T(limit), however, the traps due to glass dynamics become dominant and may be mistaken for real intermediate states. This limitation of valid temperature range prevents us to reveal the burst phase intermediate in the major folding phase since it might only be stabilized at temperatures lower than T(limit), according to experiments. Our works show that caution must be taken when studying low-temperature intermediate states by using the C(alpha) Go-models.