Lattice simulations of aggregation funnels for protein folding.

A computer model of protein aggregation competing with productive folding is proposed. Our model adapts techniques from lattice Monte Carlo studies of protein folding to the problem of aggregation. However, rather than starting with a single string of residues, we allow independently folding strings to undergo collisions and consider their interactions in different orientations. We first present some background into the nature and significance of protein aggregation and the use of lattice Monte Carlo simulations in understanding other aspects of protein folding. The results of a series of simulation experiments involving simple versions of the model illustrate the importance of considering aggregation in simulations of protein folding and provide some preliminary understanding of the characteristics of the model. Finally, we discuss the value of the model in general and of our particular design decisions and experiments. We conclude that computer simulation techniques developed to study protein folding can provide insights into protein aggregation, and that a better understanding of aggregation may in turn provide new insights into and constraints on the more general protein folding problem.     

Native disulfide bond formation in proteins

Native disulfide bond formation is critical for the proper folding of many proteins. Recent studies using newly identified protein oxidants, folding catalysts, and mutant cells provide insight into the mechanism of oxidative protein folding in vivo. This insight promises new strategies for more efficient protein production.     

Combining experiment and simulation in protein folding: closing the gap for small model systems

All-atom molecular dynamics (MD) simulations on increasingly powerful computers have been combined with experiments to characterize protein folding in detail over wider time ranges. The folding of small ultrafast folding proteins is being simulated on micros timescales, leading to improved structural predictions and folding rates. To what extent is 'closing the gap' between simulation and experiment for such systems providing insights into general mechanisms of protein folding?     

Distinguishing between sequential and nonsequentially folded proteins: implications for folding and misfolding.

We describe here an algorithm for distinguishing sequential from nonsequentially folding proteins. Several experiments have recently suggested that most of the proteins that are synthesized in the eukaryotic cell may fold sequentially. This proposed folding mechanism in vivo is particularly advantageous to the organism. In the absence of chaperones, the probability that a sequentially folding protein will misfold is reduced significantly. The problem we address here is devising a procedure that would differentiate between the two types of folding patterns. Footprints of sequential folding may be found in structures where consecutive fragments of the chain interact with each other. In such cases, the folding complexity may be viewed as being lower. On the other hand, higher folding complexity suggests that at least a portion of the polypeptide backbone folds back upon itself to form three-dimensional (3D) interactions with noncontiguous portion(s) of the chain. Hence, we look at the mechanism of folding of the molecule via analysis of its complexity, that is, through the 3D interactions formed by contiguous segments on the polypeptide chain. To computationally splice the structure into consecutively interacting fragments, we either cut it into compact hydrophobic folding units or into a set of hypothetical, transient, highly populated, contiguous fragments ("building blocks" of the structure). In sequential folding, successive building blocks interact with each other from the amino to the carboxy terminus of the polypeptide chain. Consequently, the results of the parsing differentiate between sequentially vs. nonsequentially folded chains. The automated assessment of the folding complexity provides insight into both the likelihood of misfolding and the kinetic folding rate of the given protein. In terms of the funnel free energy landscape theory, a protein that truly follows the mechanism of sequential folding, in principle, encounters smoother free energy barriers. A simple sequentially folded protein should, therefore, be less error prone and fold faster than a protein with a complex folding pattern.     

Rapid protein-folding assay using green fluorescent protein.

Formation of the chromophore of green fluorescent protein (GFP) depends on the correct folding of the protein. We constructed a "folding reporter" vector, in which a test protein is expressed as an N-terminal fusion with GFP. Using a test panel of 20 proteins, we demonstrated that the fluorescence of Escherichia coli cells expressing such GFP fusions is related to the productive folding of the upstream protein domains expressed alone. We used this fluorescent indicator of protein folding to evolve proteins that are normally prone to aggregation during expression in E. coli into closely related proteins that fold robustly and are fully soluble and functional. This approach to improving protein folding does not require functional assays for the protein of interest and provides a simple route to improving protein folding and expression by directed evolution.     

Cyclophilin 20 is involved in mitochondrial protein folding in cooperation with molecular chaperones Hsp70 and Hsp60.

We studied the role of mitochondrial cyclophilin 20 (CyP20), a peptidyl-prolyl cis-trans isomerase, in preprotein translocation across the mitochondrial membranes and protein folding inside the organelle. The inhibitory drug cyclosporin A did not impair membrane translocation of preproteins, but it delayed the folding of an imported protein in wild-type mitochondria. Similarly, Neurospora crassa mitochondria lacking CyP20 efficiently imported preproteins into the matrix, but folding of an imported protein was significantly delayed, indicating that CyP20 is involved in protein folding in the matrix. The slow folding in the mutant mitochondria was not inhibited by cyclosporin A. Folding intermediates of precursor molecules reversibly accumulated at the molecular chaperones Hsp70 and Hsp60 in the matrix. We conclude that CyP20 is a component of the mitochondrial protein folding machinery and that it cooperates with Hsp70 and Hsp60. It is speculated that peptidyl-prolyl cis-trans isomerases in other cellular compartments may similarly promote protein folding in cooperation with chaperone proteins.     

The Dominant Folding Route Minimizes Backbone Distortion in SH3

Energetic frustration in protein folding is minimized by evolution to create a smooth and robust energy landscape. As a result the geometry of the native structure provides key constraints that shape protein folding mechanisms. Chain connectivity in particular has been identified as an essential component for realistic behavior of protein folding models. We study the quantitative balance of energetic and geometrical influences on the folding of SH3 in a structure-based model with minimal energetic frustration. A decomposition of the two-dimensional free energy landscape for the folding reaction into relevant energy and entropy contributions reveals that the entropy of the chain is not responsible for the folding mechanism. Instead the preferred folding route through the transition state arises from a cooperative energetic effect. Off-pathway structures are penalized by excess distortion in local backbone configurations and contact pair distances. This energy cost is a new ingredient in the malleable balance of interactions that controls the choice of routes during protein folding.     

Förster resonance energy transfer as a probe of membrane protein folding

The folding reaction of a β-barrel membrane protein, outer membrane protein A (OmpA), is probed with F?rster resonance energy transfer (FRET) experiments. Four mutants of OmpA were generated in which the donor fluorophore, tryptophan, and acceptor molecule, a naphthalene derivative, are placed in various locations on the protein to report the evolution of distances across the bilayer and across the protein pore during a folding event. Analysis of the FRET efficiencies reveals three timescales for tertiary structure changes associated with insertion and folding into a synthetic bilayer. A narrow pore forms during the initial stage of insertion, followed by bilayer traversal. Finally, a long-time component is attributed to equilibration and relaxation, and may involve global changes such as pore expansion and strand extension. These results augment the existing models that describe concerted insertion and folding events, and highlight the ability of FRET to provide insight into the complex mechanisms of membrane protein folding. This article is part of a Special Issue entitled: Membrane protein structure and function.     

Regulation of mRNA translation by protein folding in the endoplasmic reticulum

The process of protein secretion is intimately linked to the rate and potential of proper folding and assembly of secretory proteins. The efficiency of protein folding is communicated to the cytoplasm via several signal transduction pathways. This regulates the rate of polypeptide chain synthesis and induction of genes encoding functions that reduce protein-folding load on the endoplasmic reticulum (ER). This review summarizes recent insights into the mechanisms that couple protein translation with protein folding in the ER.     

Conversion of two-state to multi-state folding kinetics on fusion of two protein foldons

Chymotrypsin inhibitor 2 (CI2) is the archetypal single-foldon protein that folds in simple two-state kinetics without the accumulation of a folding intermediate. To model the effects of fusion of single foldons to give a multi-foldon protein, we engineered a "double-CI2" protein, in which another CI2 polypeptide was inserted into the loop region of the parent CI2. CD and HSQC spectra demonstrated that while the double-CI2 protein adopted two kinds of native conformations, CI2-like structure was almost preserved in both the domains of double-CI2. In the folding kinetic studies, double-CI2 exhibited a remarkable rollover of the observed folding rates at low denaturant concentrations, indicating that double-CI2 accumulated a kinetic folding intermediate. The different folding mechanisms between WT-CI2 and double-CI2 support the present view that protein size or number of domains is an important determinant for formation of folding intermediates.     

Energy landscape remodeling mechanism of Hsp70-chaperone-accelerated protein folding

Hsp70 chaperone is one of the key protein machines responsible for the quality control of protein production in cells. Facilitating in vivo protein folding by counteracting misfolding and aggregation is the essence of its biological function. Although the allosteric cycle during its functional actions has been well characterized both experimentally and computationally, the mechanism by which Hsp70 assists protein folding is still not fully understood. In this work, we studied the Hsp70-mediated folding of model proteins with rugged energy landscape by using molecular simulations. Different from the canonical scenario of Hsp70 functioning, which assumes that folding of substrate proteins occurs spontaneously after releasing from chaperones, our results showed that the substrate protein remains in contacts with the chaperone during its folding process. The direct chaperone-substrate interactions in the open conformation of Hsp70 tend to shield the substrate sites prone to form non-native contacts, which therefore avoids the frustrated folding pathway, leading to a higher folding rate and less probability of misfolding. Our results suggest that in addition to the unfoldase and holdase functions widely addressed in previous studies, Hsp70 can facilitate the folding of its substrate proteins by remodeling the folding energy landscape and directing the folding processes, demonstrating the foldase scenario. These findings add new, to our knowledge, insights into the general molecular mechanisms of chaperone-mediated protein folding.     

Predicting protein folding pathways

A structured folding pathway, which is a time ordered sequence of folding events, plays an important role in the protein folding process and hence, in the conformational search. Pathway prediction, thus gives more insight into the folding process and is a valuable guiding tool to search the conformation space. In this paper, we propose a novel 'unfolding' approach to predict the folding pathway. We apply graph-based methods on a weighted secondary structure graph of a protein to predict the sequence of unfolding events. When viewed in reverse this yields the folding pathway. We demonstrate the success of our approach on several proteins whose pathway is partially known.     

Malleability of the folding mechanism of the outer membrane protein PagP: parallel pathways and the effect of membrane elasticity

Understanding the interactions between membrane proteins and the lipid bilayer is key to increasing our ability to predict and tailor the folding mechanism, structure and stability of membrane proteins. Here, we have investigated the effects of changing the membrane composition and the relative concentrations of protein and lipid on the folding mechanism of the bacterial outer membrane protein PagP. The folding pathway, monitored by tryptophan fluorescence, was found to be characterized by a burst phase, representing PagP adsorption to the liposome surface, followed by a time course that reflects the folding and insertion of the protein into the membrane. In 1,2-dilauroyl-sn-glycero-3-phosphocholine (diC(12:0)PC) liposomes, the post-adsorption time course fits well to a single exponential at high lipid-to-protein ratios (LPRs), but at low LPRs, a second exponential phase with a slower folding rate constant is observed. Interrupted refolding assays demonstrated that the two exponential phases reflect the presence of parallel folding pathways. Partitioning between these pathways was found to be modulated by the elastic properties of the membrane. Folding into mixed 1,2-dilauroyl-sn-glycero-3-phosphoethanolamine:diC(12:0)PC liposomes resulted in a decrease in PagP adsorption to the liposomes and a switch to the slower folding pathway. By contrast, inclusion of 1,2-dilauroyl-sn-glycero-3-phosphoserine into diC(12:0)PC liposomes resulted in a decrease in the folding rate of the fast pathway. The results highlight the effect of lipid composition in tailoring the folding mechanism of a membrane protein, revealing that membrane proteins have access to multiple, competing folding routes to a unique native structure.     

Identification and characterization of folding inhibitors of hen egg lysozyme: an example of a new paradigm of drug design

Studies of protein folding indicate the presence of native contacts in the denatured state, giving rise to folding elements which contribute to the accomplishment of the native state. The possibility of finding molecules which can interact with specific folding elements of a target protein preventing it from reaching its native state, and hence from becoming biologically active, is particularly attractive. The notion that folding elements not only provide molecular recognition directing the folding process, but also have conserved sequence, implies that targeting such elements will make protein folding inhibitors less susceptible to mutations which, in many cases, abrogate drug effects. The folding-inhibition strategy can lead to a truly novel and rational approach to drug design, aside from providing new insight into folding. This is illustrated in the case of hen egg lysozyme.     

An adaptable standard for protein export from the endoplasmic reticulum

To provide an integrated view of endoplasmic reticulum (ER) function in protein export, we have described the interdependence of protein folding energetics and the adaptable biology of cellular protein folding and transport through the exocytic pathway. A simplified treatment of the protein homeostasis network and a formalism for how this network of competing pathways interprets protein folding kinetics and thermodynamics provides a framework for understanding cellular protein trafficking. We illustrate how folding and misfolding energetics, in concert with the adjustable biological capacities of the folding, degradation, and export pathways, collectively dictate an adaptable standard for protein export from the ER. A model of folding for export (FoldEx) establishes that no single feature dictates folding and transport efficiency. Instead, a network view provides insight into the basis for cellular diversity, disease origins, and protein homeostasis, and predicts strategies for restoring protein homeostasis in protein-misfolding diseases.     

What determines protein folding type? An investigation of intrinsic structural properties and its implications for understanding folding mechanisms

Protein folding experiments demonstrate that the folding behaviors of many proteins can be roughly classified into two types: two-state kinetics and multi-state kinetics. Although the two types of protein folding kinetics have been observed for a long time, what determines the folding type of a protein is still largely unclear. The present work performed a comparative study based on a dataset of 43 two-state and 42 multi-state folders at different levels of proteins' intrinsic properties from the simplest sequence length to native structure topology. The results show that protein's amino acids composition and the long-range interaction-based topological complexity rather than secondary structure contents are the major determinants of protein folding type. Furthermore, a sequence-based folding type prediction achieved an accuracy of more than 80%. These findings implicate that there is no clear boundary between secondary and tertiary structure formation during the protein folding process and support the existence of a continuum of folding mechanism between the two ends of hierarchic and nucleation folding scenarios.     

Matching simulation and experiment: a new simplified model for simulating protein folding.

Simulations of simplified protein folding models have provided much insight into solving the protein folding problem. We propose here a new off-lattice bead model, capable of simulating several different fold classes of small proteins. We present the sequence for an alpha/beta protein resembling the IgG-binding proteins L and G. The thermodynamics of the folding process for this model are characterized using the multiple multihistogram method combined with constant-temperature Langevin simulations. The folding is shown to be highly cooperative, with chain collapse nearly accompanying folding. Two parallel folding pathways are shown to exist on the folding free energy landscape. One pathway contains an intermediate--similar to experiments on protein G, and one pathway contains no intermediates-similar to experiments on protein L. The folding kinetics are characterized by tabulating mean-first passage times, and we show that the onset of glasslike kinetics occurs at much lower temperatures than the folding temperature. This model is expected to be useful in many future contexts: investigating questions of the role of local versus nonlocal interactions in various fold classes, addressing the effect of sequence mutations affecting secondary structure propensities, and providing a computationally feasible model for studying the role of solvation forces in protein folding.     

Critical role of beta-hairpin formation in protein G folding

Comparison of the folding mechanisms of proteins with similar structures but very different sequences can provide fundamental insights into the determinants of protein folding mechanisms. Despite very little sequence similarity, the approximately 60 residue IgG binding domains of protein G and protein L both consist of a single helix packed against a four-stranded sheet formed by two symmetrically disposed beta-hairpins. We demonstrate that, as in the case of protein L, one of the two beta-turns of protein G is formed and the other disrupted in the folding transition state. Unlike protein L, however, in protein G it is the second beta-turn that is formed in the folding transition state ensemble. Substitution of an Asp residue by Ala in protein G that eliminates an i,i+2 side chain-main chain hydrogen bond in the second beta-turn slows the folding rate approximately 20-fold but has virtually no effect on the unfolding rate. Taken together with previous results, these findings suggest that the presence of an intact beta-turn in the folding transition state is a consequence of the overall topology of protein L and protein G, but the particular hairpin that is formed is determined by the detailed interatomic interactions that determine the free energies of formation of the isolated beta-hairpins.     

The dynamical contact order: Protein folding rate parameters based on quantum conformational transitions

Protein folding is regarded as a quantum transition between the torsion states of a polypeptide chain. According to the quantum theory of conformational dynamics, we propose the dynamical contact order (DCO) defined as a characteristic of the contact described by the moment of inertia and the torsion potential energy of the polypeptide chain between contact residues. Consequently, the protein folding rate can be quantitatively studied from the point of view of dynamics. By comparing theoretical calculations and experimental data on the folding rate of 80 proteins, we successfully validate the view that protein folding is a quantum conformational transition. We conclude that (i) a correlation between the protein folding rate and the contact inertial moment exists; (ii) multi-state protein folding can be regarded as a quantum conformational transition similar to that of two-state proteins but with an intermediate delay. We have estimated the order of magnitude of the time delay; (iii) folding can be classified into two types, exergonic and endergonic. Most of the two-state proteins with higher folding rate are exergonic and most of the multi-state proteins with low folding rate are endergonic. The folding speed limit is determined by exergonic folding.     

Direct correlation between proteins' folding rates and their amino acid compositions: an ab initio folding rate prediction

Discovering the mechanism of protein folding, in molecular biology, is a great challenge. A key step to this end is to find factors that correlate with protein folding rates. Over the past few years, many empirical parameters, such as contact order, long-range order, total contact distance, secondary structure contents, have been developed to reflect the correlation between folding rates and protein tertiary or secondary structures. However, the correlation between proteins' folding rates and their amino acid compositions has not been explored. In the present work, we examined systematically the correlation between proteins' folding rates and their amino acid compositions for two-state and multistate folders and found that different amino acids contributed differently to the folding progress. The relation between the amino acids' molecular weight and degeneracy and the folding rates was examined, and the role of hydrophobicity in the protein folding process was also inspected. As a consequence, a new indicator called composition index was derived, which takes no structure factors into account and is merely determined by the amino acid composition of a protein. Such an indicator is found to be highly correlated with the protein's folding rate (r > 0.7). From the results of this work, three points of concluding remarks are evident. (1) Two-state folders and multistate folders have different rate-determining amino acids. (2) The main determining information of a protein's folding rate is largely reflected in its amino acid composition. (3) Composition index may be the best predictor for an ab initio protein folding rate prediction directly from protein sequence from the standpoint of practical application.