Local and non-local native topologies reveal the underlying folding landscape of proteins

Due to Plaxco, Simons, Baker and others, it is now well known that the two-state single domain protein folding rate is fairly well predicted from knowledge of the topology of the native structure. Plaxco et al found that the folding rates of two-state proteins correlate with the average degree to which native contacts are 'local' within the chain sequence: fast-folders usually have mostly local structures. Here, we dissected the native topology further by focusing on non-local and local contacts using lower and upper bounds of allowable sequence separation in computing the average contact order. We analyzed non-local and local contacts of 82 two-state proteins whose experimental folding rates span over six orders of magnitude. We observed that both the number of non-local contacts and the average sequence separation of non-local contacts (non-local CO) are both negatively correlated with the folding rate, showing that the non-local contacts dominate the barrier-crossing process. Surprisingly, the local contact orders of the proteins also correlate with the folding rates. However, this correlation shows a strong positive trend indicating the role of a diffusive search in the denatured basin.     

Analysis of the folding pathway of chymotrypsin inhibitor by correlation of phi-values with inter-residue contacts

The transition state for folding of chymotrypsin inhibitor 2 (CI2) is investigated by correlating Phi-values with inter-residue contacts. In agreement with former work, the strongest consolidation of secondary structure is found in the alpha-helix. There are correlations for tertiary structure interactions between the residues Leu49, Ile57 and the helix which have been suggested to represent the main components of the nucleation site of CI2 folding. However, correlations for tertiary structure interactions of comparable magnitude are also found in the helix-strand2-strand1-motif and between strand3and strand4. Copyright 1999 Academic Press.     

Recent successes of the energy landscape theory of protein folding and function

Protein folding and binding can be understood using energy landscape theory. When seeming deviations from the predictions of the funnel hypothesis are found, landscape theory helps us locate the cause. Sometimes the deviation reflects symmetry effects, allowing extra degeneracies to occur. Such effects seem to explain some kinetic anomalies in helical bundles. When binding processes were found to use apparently non-funneled landscapes this was traced to an inadequate understanding of biomolecular forces. The discrepancy allowed the discovery of new water-mediated forces - some of which act between hydrophilic residues. Introducing such forces into the algorithms greatly improves the quality of structure predictions.     

Navigating the RNA folding landscape

Analysis of individual RNA folding reactions reveals that, as in proteins, cooperative interactions selectively drive RNA toward its biologically active, native conformation. This new work establishes a platform for future investigations of the physical principles underlying the assembly of large RNA enzymes.     

Energetics of the native energy landscape of a two-domain calcium sensor protein: distinct folding features of the two domains

The protein folding energy landscape allows a thorough understanding of the protein folding problem which in turn helps in understanding various aspects of biological functions. Characterizing the cooperative unfolding units and the intermediates along the folding funnel of a protein is a challenging task. In this paper, we investigated the native energy landscape of EhCaBP, a calcium sensor, belonging to the same EF-hand superfamily as calmodulin. EhCaBP is a two-domain EF-hand protein consisting of two EF-hands in each domain and binding to four Ca2+ cations. Native-state hydrogen exchange (HX) was used to assess the folding features of the landscape and also to throw light on the structure-folding function paradigm of calcium sensor proteins. HX measurements under the EX2 regime provided the thermodynamic information about the protein folding events under native conditions. HX studies revealed that the unfolding of EhCaBP is not a two-state process. Instead, it proceeds through cooperative units. The C-terminal domain exhibits less denaturant dependence than the N-terminal domain, suggesting that the former is dominated by local fluctuations. It is interesting to note that the N- and C-terminal domains of EhCaBP have distinct folding features. In fact, these observed differences can regulate the domain-dependent target recognition of two-domain Ca2+ sensor proteins.     

Probing the cytochrome c' folding landscape

The folding kinetics of R. palustris cytochrome c' (cyt c') have been monitored by heme absorption and native Trp72 fluorescence at pH 5. The Trp72 fluorescence burst signal suggests early compaction of the polypeptide ensemble. Analysis of heme transient absorption spectra reveals deviations from two-state behavior, including a prominent slow phase that is accelerated by the prolyl isomerase cyclophilin. A nonnative proline configuration (Pro21) likely interferes with the formation of the helical bundle surrounding the heme.     

Morphogenesis of a protein: folding pathways and the energy landscape

Current knowledge on the reaction whereby a protein acquires its native three-dimensional structure was obtained by and large through characterization of the folding mechanism of simple systems. Given the multiplicity of amino acid sequences and unique folds, it is not so easy, however, to draw general rules by comparing folding pathways of different proteins. In fact, quantitative comparison may be jeopardized not only because of the vast repertoire of sequences but also in view of a multiplicity of structures of the native and denatured states. We have tackled the problem of the relationships between the sequence information and the folding pathway of a protein, using a combination of kinetics, protein engineering and computational methods, applied to relatively simple systems. Our strategy has been to investigate the folding mechanism determinants using two complementary approaches, i.e. (i) the study of members of the same family characterized by a common fold, but substantial differences in amino acid sequence, or (ii) heteromorphic pairs characterized by largely identical sequences but with different folds. We discuss some recent data on protein-folding mechanisms by presenting experiments on different members of the PDZ domain family and their circularly permuted variants. Characterization of the energetics and structures of intermediates and TSs (transition states), obtained by Φ-value analysis and restrained MD (molecular dynamics) simulations, provides a glimpse of the malleability of the dynamic states and of the role of the topology of the native states and of the denatured states in dictating folding and misfolding pathways.     

Pathway of oxidative folding of bovine alpha-interferon: predominance of native disulfide-bonded folding intermediates

Bovine alpha-interferon (BoINF-alpha) is a single polypeptide protein containing 166 amino acids, two disulfide bonds (Cys1-Cys99 and Cys29-Cys138), and five stretches of alpha-helical structure. The pathway of oxidative folding of BoINF-alpha has been investigated here. Of the eight possible one- and two-disulfide isomers, only two nativelike one-disulfide isomers, BoINF-alpha (Cys1-Cys99) and BoINF-alpha (Cys29-Cys138), predominate as intermediates along the folding pathway. More strikingly, alpha-helical structures formed almost quantitatively before any detectable formation of a disulfide bond. This is demonstrated by the observation that fully reduced BoINF-alpha (starting material of oxidative folding) and reduced carboxymethylated BoINF-alpha both exhibit alpha-helical structure content indistinguishable form that of native BoINF-alpha. The folding mechanism of BoINF-alpha appears to be compatible with the framework model, in which secondary structures fold first, followed by docking (compaction) of preformed secondary structural elements yielding the native structure.     

Evolution of the folding landscape of effector caspases

Caspases are a family of cysteinyl proteases that control programmed cell death and maintain homeostasis in multicellular organisms. The caspase family is an excellent model to study protein evolution because all caspases are produced as zymogens (procaspases [PCPs]) that must be activated to gain full activity; the protein structures are conserved through hundreds of millions of years of evolution; and some allosteric features arose with the early ancestor, whereas others are more recent evolutionary events. The apoptotic caspases evolved from a common ancestor (CA) into two distinct subfamilies: monomers (initiator caspases) or dimers (effector caspases). Differences in activation mechanisms of the two subfamilies, and their oligomeric forms, play a central role in the regulation of apoptosis. Here, we examine changes in the folding landscape by characterizing human effector caspases and their CA. The results show that the effector caspases unfold by a minimum three-state equilibrium model at pH 7.5, where the native dimer is in equilibrium with a partially folded monomeric (PCP-7, CA) or dimeric (PCP-6) intermediate. In comparison, the unfolding pathway of PCP-3 contains both oligomeric forms of the intermediate. Overall, the data show that the folding landscape was first established with the CA and was retained for >650 million years. Partially folded monomeric or dimeric intermediates in the ancestral ensemble provide mechanisms for evolutionary changes that affect stability of extant caspases. The conserved folding landscape allows for the fine-tuning of enzyme stability in a species-dependent manner while retaining the overall caspase–hemoglobinase fold.     

Interdomain contacts control folding of transcription factor RfaH

Escherichia coli RfaH activates gene expression by tethering the elongating RNA polymerase to the ribosome. This bridging action requires a complete refolding of the RfaH C-terminal domain (CTD) from an α-helical hairpin, which binds to the N-terminal domain (NTD) in the free protein, to a β-barrel, which interacts with the ribosomal protein S10 following RfaH recruitment to its target operons. The CTD forms a β-barrel when expressed alone or proteolytically separated from the NTD, indicating that the α-helical state is trapped by the NTD, perhaps co-translationally. Alternatively, the interdomain contacts may be sufficient to drive the formation of the α-helical form. Here, we use functional and NMR analyses to show that the denatured RfaH refolds into the native state and that RfaH in which the order of the domains is reversed is fully functional in vitro and in vivo. Our results indicate that all information necessary to determine its fold is encoded within RfaH itself, whereas accessory factors or sequential folding of NTD and CTD during translation are dispensable. These findings suggest that universally conserved RfaH homologs may change folds to accommodate diverse interaction partners and that context-dependent protein refolding may be widespread in nature.     

Protein folding in the landscape perspective: Chevron plots and non-arrhenius kinetics

We use two simple models and the energy landscape perspective to study protein folding kinetics. A major challenge has been to use the landscape perspective to interpret experimental data, which requires ensemble averaging over the microscopic trajectories usually observed in such models. Here, because of the simplicity of the model, this can be achieved. The kinetics of protein folding falls into two classes: multiple-exponential and two-state (single-exponential) kinetics. Experiments show that two-state relaxation times have “chevron plot” dependences on denaturant and non-Arrhenius dependences on temperature. We find that HP and HP+ models can account for these behaviors. The HP model often gives bumpy landscapes with many kinetic traps and multiple-exponental behavior, whereas the HP+ model gives more smooth funnels and two-state behavior. Multiple-exponential kinetics often involves fast collapse into kinetic traps and slower barrier climbing out of the traps. Two-state kinetics often involves entropic barriers where conformational searching limits the folding speed. Transition states and activation barriers need not define a single conformation; they can involve a broad ensemble of the conformations searched on the way to the native state. We find that unfolding is not always a direct reversal of the folding process. Proteins 30:2–33, 1998. © 1998 Wiley-Liss, Inc.     

Nucleation‐based prediction of the protein folding rate and its correlation with the folding nucleus size

The problem of protein self‐organization is in the focus of current molecular biology studies. Although the general principles are understood, many details remain unclear. Specifically, protein folding rates are of interest because they dictate the rate of protein aggregation which underlies many human diseases. Here we offer predictions of protein folding rates and their correlation with folding nucleus sizes. We calculated free energies of the transition state and sizes of folding nuclei for 84 proteins and peptides whose other parameters were measured at the point of thermodynamic equilibrium between their unfolded and native states. We used the dynamic programming method where each residue was considered to be either as folded as in its native state or completely disordered. The calculated and measured folding rates showed a good correlation at the temperature mid‐transition point (the correlation coefficient was 0.75). Also, we pioneered in demonstrating a moderate (‐0.57) correlation coefficient between the calculated sizes of folding nuclei and the folding rates. Predictions made by different methods were compared. The established good correlation between the estimated free energy barrier and the experimentally found folding rate of each studied protein/peptide indicates that our model gives reliable results for the considered data set. Proteins 2012; © 2012 Wiley Periodicals, Inc.     

The folding funnel landscape for the peptide Met-enkephalin

We study the free energy landscape of the small peptide Met-enkephalin. Our data were obtained from a generalized-ensemble Monte Carlo simulation taking the interactions among all atoms into account. We show that the free energy landscape resembles that of a funnel, indicating that this peptide is a good folder. Our work demonstrates that the energy landscape picture and folding concept, developed in the context of simplified protein models, can also be used to describe the folding in more realistic models.     

Finding the path in an RNA folding landscape

In this issue of Structure, Reymond et?al. (2010) combine molecular and computational biology approaches to provide structural details for intermediates in the folding pathway of the hepatitis delta virus ribozyme.     

Probing the self-association, intermolecular contacts, and folding propensity of amelogenin

Amelogenins are an intrinsically disordered protein family that plays a major role in the development of tooth enamel, one of the most highly mineralized materials in nature. Monomeric porcine amelogenin possesses random coil and residual secondary structures, but it is not known which sequence regions would be conformationally attractive to potential enamel matrix targets such as other amelogenins (self-assembly), other matrix proteins, cell surfaces, or biominerals. To address this further, we investigated recombinant porcine amelogenin (rP172) using "solvent engineering" techniques to simultaneously promote native-like structure and induce amelogenin oligomerization in a manner that allows identification of intermolecular contacts between amelogenin molecules. We discovered that in the presence of 2,2,2-trifluoroethanol (TFE) significant folding transitions and stabilization occurred primarily within the N- and C-termini, while the polyproline Type II central domain was largely resistant to conformational transitions. Seven Pro residues (P2, P127, P130, P139, P154, P157, P162) exhibited conformational response to TFE, and this indicates these Pro residues act as folding enhancers in rP172. The remaining Pro residues resisted TFE perturbations and thus act as conformational stabilizers. We also noted that TFE induced rP172 self-association via the formation of intermolecular contacts involving P4-H6, V19-P33, and E40-T58 regions of the N-terminus. Collectively, these results confirm that the N- and C-termini of amelogenin are conformationally responsive and represent potential interactive sites for amelogenin-target interactions during enamel matrix mineralization. Conversely, the Pro, Gln central domain is resistant to folding and this may have important functional significance for amelogenin.     

A minimal physically realistic protein-like lattice model: designing an energy landscape that ensures all-or-none folding to a unique native state

A simple protein model restricted to the face-centered cubic lattice has been studied. The model interaction scheme includes attractive interactions between hydrophobic (H) residues, repulsive interactions between hydrophobic and polar (P) residues, and orientation-dependent P-P interactions. Additionally, there is a potential that favors extended beta-type conformations. A sequence has been designed that adopts a native structure, consisting of an antiparallel, six-member Greek-key beta-barrel with protein-like structural degeneracy. It has been shown that the proposed model is a minimal one, i.e., all the above listed types of interactions are necessary for cooperative (all-or-none) type folding to the native state. Simulations were performed via the Replica Exchange Monte Carlo method and the numerical data analyzed via a multihistogram method.     

Energy landscape of a native protein: Jumping-among-minima model

We have investigated energy landscape of human lysozyme in its native state by using principal component analysis and a model, jumping-among-minima (JAM) model. These analyses are applied to 1 nsec molecular dynamics trajectory of the protein in water. An assumption embodied in the JAM model allows us to divide protein motions into intra-substate and inter-substate motions. By examining intra-substate motions, it is shown that energy surfaces of individual conformational substates are nearly harmonic and mutually similar. As a result of principal component analysis and JAM model analysis, protein motions are shown to consist of three types of collective modes, multiply hierarchical modes, singly hierarchical modes, and harmonic modes. Multiply hierarchical modes, the number of which accounts only for 0.5% of all modes, dominate contributions to total mean-square atomic fluctuation. Inter-substate motions are observed only in a small-dimensional subspace spanned by the axes of multiplyhierarchical and singly hierarchical modes. Inter-substate motions have two notable time components: faster component seen within 200 psec and slower component. The former involves transitions among the conformational substates of the low-level hierarchy, whereas the latter involves transitions of the higher level substates observed along the first four multiply hierarchical modes. We also discuss dependence of the subspace, which contains conformational substates, on time duration of simulation. Proteins 33:496–517, 1998. © 1998 Wiley-Liss, Inc.     

Measuring the folding landscape of a harmonically constrained biopolymer

Pioneering studies have shown that the probability distribution of opening length for a DNA hairpin, recorded under constant force using an optical trap, can be used to reconstruct the energy landscape of the transition. However, measurements made under constant force are subject to some limitations. Under constant force a system with a sufficiently high energy barrier spends most of its time in the closed or open conformation, with relatively few statistics collected in the transition state region. We describe a measurement scheme in which the system is driven progressively through the transition by an optical trap and an algorithm is used to extract the energy landscape of the transition from the fluctuations recorded during this process. We illustrate this technique in simulations and demonstrate its effectiveness in experiments on a DNA hairpin. We find that the combination of this technique with the use of short DNA handles facilitates a high-resolution measurement of the hairpin's folding landscape with a very short measurement time.