Identical RNA-Protein Interactions in Vivo and in Vitro and a Scheme of Folding the Newly Synthesized Proteins by Ribosomes

A distinct three-dimensional shape of rRNA inside the ribosome is required for the peptidyl transfer activity of its peptidyltransferase center (PTC). In contrast, even the in vitro transcribed PTC RNA interacts with unfolded protein(s) at about five sites to let them attain their native states. We found that the same set of conserved nucleotides in the PTC interact identically with nascent and chemically unfolded proteins in vivo and in vitro, respectively. The time course of this interaction, difficult to follow in vivo, was observed in vitro. It suggested nucleation of folding of cytosolic globular proteins vectorially from hydrophilic N to hydrophobic C termini, consistent with our discovery of a regular arrangement of cumulative hydrophobic indices of the peptide segments of cytosolic proteins from N to C termini. Based on this observation, we propose a model here for the nucleation of folding of the nascent protein chain by the PTC.     

Folding of β-barrel membrane proteins in lipid bilayers — Unassisted and assisted folding and insertion

In cells, β-barrel membrane proteins are transported in unfolded form to an outer membrane into which they fold and insert. Model systems have been established to investigate the mechanisms of insertion and folding of these versatile proteins into detergent micelles, lipid bilayers and even synthetic amphipathic polymers. In these experiments, insertion into lipid membranes is initiated from unfolded forms that do not display residual β-sheet secondary structure. These studies therefore have allowed the investigation of membrane protein folding and insertion in great detail. Folding of β-barrel membrane proteins into lipid bilayers has been monitored from unfolded forms by dilution of chaotropic denaturants that keep the protein unfolded as well as from unfolded forms present in complexes with molecular chaperones from cells. This review is aimed to provide an overview of the principles and mechanisms observed for the folding of β-barrel transmembrane proteins into lipid bilayers, the importance of lipid–protein interactions and the function of molecular chaperones and folding assistants. This article is part of a Special Issue entitled: Lipid–protein interactions.     

Direct Observation of Parallel Folding Pathways Revealed Using a Symmetric Repeat Protein System

Although progress has been made to determine the native fold of a polypeptide from its primary structure, the diversity of pathways that connect the unfolded and folded states has not been adequately explored. Theoretical and computational studies predict that proteins fold through parallel pathways on funneled energy landscapes, although experimental detection of pathway diversity has been challenging. Here, we exploit the high translational symmetry and the direct length variation afforded by linear repeat proteins to directly detect folding through parallel pathways. By comparing folding rates of consensus ankyrin repeat proteins (CARPs), we find a clear increase in folding rates with increasing size and repeat number, although the size of the transition states (estimated from denaturant sensitivity) remains unchanged. The increase in folding rate with chain length, as opposed to a decrease expected from typical models for globular proteins, is a clear demonstration of parallel pathways. This conclusion is not dependent on extensive curve-fitting or structural perturbation of protein structure. By globally fitting a simple parallel-Ising pathway model, we have directly measured nucleation and propagation rates in protein folding, and have quantified the fluxes along each path, providing a detailed energy landscape for folding. This finding of parallel pathways differs from results from kinetic studies of repeat-proteins composed of sequence-variable repeats, where modest repeat-to-repeat energy variation coalesces folding into a single, dominant channel. Thus, for globular proteins, which have much higher variation in local structure and topology, parallel pathways are expected to be the exception rather than the rule.     

Structural Characteristic of the Initial Unfolded State on Refolding Determines Catalytic Efficiency of the Folded Protein in Presence of Osmolytes

Osmolytes are low molecular weight organic molecules accumulated by organisms to assist proper protein folding, and to provide protection to the structural integrity of proteins under denaturing stress conditions. It is known that osmolyte-induced protein folding is brought by unfavorable interaction of osmolytes with the denatured/unfolded states. The interaction of osmolyte with the native state does not significantly contribute to the osmolyte-induced protein folding. We have therefore investigated if different denatured states of a protein (generated by different denaturing agents) interact differently with the osmolytes to induce protein folding. We observed that osmolyte-assisted refolding of protein obtained from heat-induced denatured state produces native molecules with higher enzyme activity than those initiated from GdmCl- or urea-induced denatured state indicating that the structural property of the initial denatured state during refolding by osmolytes determines the catalytic efficiency of the folded protein molecule. These conclusions have been reached from the systematic measurements of enzymatic kinetic parameters (K _m and k _cat), thermodynamic stability (T _m and ΔH _m) and secondary and tertiary structures of the folded native proteins obtained from refolding of various denatured states (due to heat-, urea- and GdmCl-induced denaturation) of RNase-A in the presence of various osmolytes.     

Role of Protein Stabilizers on the Conformation of the Unfolded State of Cytochrome c and Its Early Folding Kinetics: INVESTIGATION AT SINGLE MOLECULAR RESOLUTION*

An insight into the conformation and dynamics of unfolded and early intermediate states of a protein is essential to understand the mechanism of its aggregation and to design potent inhibitor molecules. Fluorescence correlation spectroscopy has been used to study the effects of several model protein stabilizers on the conformation of the unfolded state and early folding dynamics of tetramethyl rhodamine-labeled cytochrome c from Saccharomyces cerevisiae at single molecular resolution. Special attention has been given to arginine, which is a widely used stabilizer for improving refolding yield of different proteins. The value of the hydrodynamic radius (r_H) obtained by analyzing the intensity fluctuations of the diffusing molecules has been found to increase in a two-state manner as the protein is unfolded by urea. The results further show that the presence of arginine and other protein stabilizers favors a relatively structured conformation of the unfolded states (r_H of 29 Å) over an extended one (r_H of 40 Å), which forms in their absence. Also, the time constant of a kinetic component (τ_R) of about 30 μs has been observed by analyzing the correlation functions, which represents formation of a collapsed state. This time constant varies with urea concentration representing an inverted Chevron plot that shows a roll-over and behavior in the absence of arginine. To the best of our knowledge, this is one of the first applications of fluorescence correlation spectroscopy to study direct folding kinetics of a protein.     

Sequential folding of transfer RNA: A nuclear magnetic resonance study of successively longer tRNA fragments with a common 5′ end

Most folding studies on proteins and nucleic acids have been addressed to the transition between the folded and unfolded states of an intact molecule, where an entire residue sequence is present during the folding event. However, since these polymers are synthesized sequentially from one terminus to the other in vivo, their folding pathways may be influenced greatly by the sequential appearance of the residues as a function of time.The three-dimensional structure of yeast tRNA^Phe in the crystalline state is correlated with 360 MHz proton nuclear magnetic resonances from three fragments plus an intact molecule of the tRNA that share a common 5′ end and are in a solution condition similar to that of the crystal structure. This has allowed identification of folded structures present in the fragments and presumably present in the growing tRNA molecule as it is being synthesized from the 5′ end. The experiments show that only the correct stems are formed in the fragments; no additional or competing helical region is produced. This suggests that in the biosynthesis of this tRNA, correct folding of helical stems occurs before the entire molecule is formed. Further, some of the tertiary interactions (hydrogen bonds) found in the crystal structure are also probably present before the synthesis is completed. These findings are generalized to consider the precursor of the tRNA as well as other tRNAs.     

Single-Molecule Studies of the Im7 Folding Landscape

Under appropriate conditions, the four-helical Im7 (immunity protein 7) folds from an ensemble of unfolded conformers to a highly compact native state via an on-pathway intermediate. Here, we investigate the unfolded, intermediate, and native states populated during folding using diffusion single-pair fluorescence resonance energy transfer by measuring the efficiency of energy transfer (or proximity or P ratio) between pairs of fluorophores introduced into the side chains of cysteine residues placed in the center of helices 1 and 4, 1 and 3, or 2 and 4. We show that while the native states of each variant give rise to a single narrow distribution with high P values, the distributions of the intermediates trapped at equilibrium (denoted I^eqm) are fitted by two Gaussian distributions. Modulation of the folding conditions from those that stabilize the intermediate to those that destabilize the intermediate enabled the distribution of lower P value to be assigned to the population of the unfolded ensemble in equilibrium with the intermediate state. The reduced stability of the I^eqm variants allowed analysis of the effect of denaturant concentration on the compaction and breadth of the unfolded state ensemble to be quantified from 0 to 6 M urea. Significant compaction is observed as the concentration of urea is decreased in both the presence and absence of sodium sulfate, as previously reported for a variety of proteins. In the presence of Na₂SO₄ in 0 M urea, the P value of the unfolded state ensemble approaches that of the native state. Concurrent with compaction, the ensemble displays increased peak width of P values, possibly reflecting a reduction in the rate of conformational exchange among iso-energetic unfolded, but compact conformations. The results provide new insights into the initial stages of folding of Im7 and suggest that the unfolded state is highly conformationally constrained at the outset of folding.     

Molecular basis of cooperativity in protein folding IV. Core: A general cooperative folding model

The cooperative nature of the protein folding process is independent of the characteristic fold and the specific secondary structure attributes of a globular protein. A general folding/unfolding model should, therefore, be based upon structural features that transcend the peculiarities of α-helices, β-sheets, and other structural motifs found in proteins. The studies presented in this paper suggest that a single structural characteristic common to all globular proteins is essential for cooperative folding. The formation of a partly folded state from the native state results in the exposure to solvent of two distinct regions: (1) the portions of the protein that are unfolded; and (2) the “complementary surfaces,” located in the regions of the protein that remain folded. The cooperative character of the folding/unfolding transition is determined largely by the energetics of exposing complementary surface regions to the solvent. By definition, complementary regions are present only in partly folded states; they are absent from the native and unfolded states. An unfavorable free energy lowers the probability of partly folded states and increases the cooperativity of the transition. In this paper we present a mathematical formulation of this behavior and develop a general cooperative folding/unfolding model, termed the “complementary region” (CORE) model. This model successfully reproduces the main properties of folding/unfolding transitions without limiting the number of partly folded states accessible to the protein, thereby permitting a systematic examination of the structural and solvent conditions under which intermediates become populated. It is shown that the CORE model predicts two-state folding/unfolding behavior, even though the two-state character is not assumed in the model. © 1993 Wiley-Liss, Inc.     

Understanding how small helical proteins fold: conformational dynamics of Im proteins relevant to their folding landscapes

Understanding the mechanism of folding of small proteins requires characterization of their starting unfolded states and any partially unfolded states populated during folding. Here, we review what is known from NMR about these states of Im7, a 4-helix bundle protein that folds via an on-pathway intermediate, and show that there is an alignment of non-native structure in urea-unfolded Im7 with the helices of native Im7 that is a consequence of hydrophobic helix-promoting residues also promoting cluster-formation in the unfolded protein. We suggest that this kind of alignment is present in other proteins and is relevant to how native state topology determines folding rates.     

A novel approach for large-scale polypeptide folding based on elastic networks using continuous optimization

We present a new computationally efficient method for large-scale polypeptide folding using coarse-grained elastic networks and gradient-based continuous optimization techniques. The folding is governed by minimization of energy based on Miyazawa-Jernigan contact potentials. Using this method we are able to substantially reduce the computation time on ordinary desktop computers for simulation of polypeptide folding starting from a fully unfolded state. We compare our results with available native state structures from Protein Data Bank (PDB) for a few de-novo proteins and two natural proteins, Ubiquitin and Lysozyme. Based on our simulations we are able to draw the energy landscape for a small de-novo protein, Chignolin. We also use two well known protein structure prediction software, MODELLER and GROMACS to compare our results. In the end, we show how a modification of normal elastic network model can lead to higher accuracy and lower time required for simulation.     

Folding network of villin headpiece subdomain

Protein folding is a complex multidimensional process that is difficult to illustrate by the traditional analyses based on one- or two-dimensional profiles. Analyses based on transition networks have become an alternative approach that has the potential to reveal detailed features of protein folding dynamics. However, due to the lack of successful reversible folding of proteins from conventional molecular-dynamics simulations, this approach has rarely been utilized. Here, we analyzed the folding network from several 10 μs conventional molecular-dynamics reversible folding trajectories of villin headpiece subdomain (HP35). The folding network revealed more complexity than the traditional two-dimensional map and demonstrated a variety of conformations in the unfolded state, intermediate states, and the native state. Of note, deep enthalpic traps at the unfolded state were observed on the folding landscape. Furthermore, in contrast to the clear separation of the native state and the primary intermediate state shown on the two-dimensional map, the two states were mingled on the folding network, and prevalent interstate transitions were observed between these two states. A more complete picture of the folding mechanism of HP35 emerged when the traditional and network analyses were considered together.     

Non-native hydrophobic interactions detected in unfolded apoflavodoxin by paramagnetic relaxation enhancement

Transient structures in unfolded proteins are important in elucidating the molecular details of initiation of protein folding. Recently, native and non-native secondary structure have been discovered in unfolded A. vinelandii flavodoxin. These structured elements transiently interact and subsequently form the ordered core of an off-pathway folding intermediate, which is extensively formed during folding of this α–β parallel protein. Here, site-directed spin-labelling and paramagnetic relaxation enhancement are used to investigate long-range interactions in unfolded apoflavodoxin. For this purpose, glutamine-48, which resides in a non-native α-helix of unfolded apoflavodoxin, is replaced by cysteine. This replacement enables covalent attachment of nitroxide spin-labels MTSL and CMTSL. Substitution of Gln-48 by Cys-48 destabilises native apoflavodoxin and reduces flexibility of the ordered regions in unfolded apoflavodoxin in 3.4 M GuHCl, because of increased hydrophobic interactions in the unfolded protein. Here, we report that in the study of the conformational and dynamic properties of unfolded proteins interpretation of spin-label data can be complicated. The covalently attached spin-label to Cys-48 (or Cys-69 of wild-type apoflavodoxin) perturbs the unfolded protein, because hydrophobic interactions occur between the label and hydrophobic patches of unfolded apoflavodoxin. Concomitant hydrophobic free energy changes of the unfolded protein (and possibly of the off-pathway intermediate) reduce the stability of native spin-labelled protein against unfolding. In addition, attachment of MTSL or CMTSL to Cys-48 induces the presence of distinct states in unfolded apoflavodoxin. Despite these difficulties, the spin-label data obtained here show that non-native contacts exist between transiently ordered structured elements in unfolded apoflavodoxin.     

Structural Interpretation of Paramagnetic Relaxation Enhancement-Derived Distances for Disordered Protein States

Paramagnetic relaxation enhancement (PRE) is a powerful technique for studying transient tertiary organizations of unfolded and partially folded proteins. The heterogeneous and dynamic nature of disordered protein states, together with the r⁻⁶ dependence of PRE, presents significant challenges for reliable structural interpretation of PRE-derived distances. Without additional knowledge of accessible conformational substates, ensemble-simulation-based protocols have been used to calculate structure ensembles that appear to be consistent with the PRE distance restraints imposed on the ensemble level with the proper r⁻⁶ weighting. However, rigorous assessment of the reliability of such protocols has been difficult without intimate knowledge of the true nature of disordered protein states. Here we utilize sets of theoretical PRE distances derived from simulated structure ensembles that represent the folded, partially folded and unfolded states of a small protein to investigate the efficacy of ensemble-simulation-based structural interpretation of PRE distances. The results confirm a critical limitation that, due to r⁻⁶ weighting, only one or a few members need to satisfy the distance restraints and the rest of the ensemble are essentially unrestrained. Consequently, calculated structure ensembles will appear artificially heterogeneous no matter whether the PRE distances are derived from the folded, partially unfolded or unfolded state. Furthermore, the nature of the heterogeneous ensembles is largely determined by the protein model employed in structure calculation and reflects little on the true nature of the underlying disordered state. These findings suggest that PRE measurements on disordered protein states alone generally do not contain enough information for a reliable structural interpretation and that the latter will require additional knowledge of accessible conformational substates. Interestingly, when a very large number of PRE measurements is available, faithful structural interpretation might be possible with intermediate ensemble sizes under ideal conditions.     

Structural analysis of kinetic folding intermediates for a TIM barrel protein, indole-3-glycerol phosphate synthase, by hydrogen exchange mass spectrometry and Gō model simulation

The structures of partially folded states appearing during the folding of a (βα)₈ TIM barrel protein, the indole-3-glycerol phosphate synthase from Sulfolobus solfataricus (sIGPS), was assessed by hydrogen exchange mass spectrometry (HX-MS) and Gō model simulations. HX-MS analysis of the peptic peptides derived from the pulse-labeled product of the sub-millisecond folding reaction from the urea-denatured state revealed strong protection in the (βα)₄ region, modest protection in the neighboring (βα)_1-3 and (βα)₅β₆ segments and no significant protection in the remaining N and C-terminal segments. These results demonstrate that this species is not a collapsed form of the unfolded state under native-favoring conditions nor is it the native state formed via fast-track folding. However, the striking contrast of these results with the strong protection observed in the (βα)_2-5β₆ region after 5 s of folding demonstrates that these species represent kinetically distinct folding intermediates that are not identical as previously thought. A re-examination of the kinetic folding mechanism by chevron analysis of fluorescence data confirmed distinct roles for these two species: the burst-phase intermediate is predicted to be a misfolded, off-pathway intermediate, while the subsequent 5 s intermediate corresponds to an on-pathway equilibrium intermediate. Comparison with the predictions using a C^α Gō model simulation of the kinetic folding reaction for sIGPS shows good agreement with the core of the structure offering protection against exchange in the on-pathway intermediate(s). Because the native-centric Gō model simulations do not explicitly include sequence-specific information, the simulation results support the hypothesis that the topology of TIM barrel proteins is a primary determinant of the folding free energy surface for the productive folding reaction. The early misfolding reaction must involve aspects of non-native structure not detected by the Gō model simulation.     

Tryptophan-tryptophan energy migration as a tool to follow apoflavodoxin folding

Submolecular details of Azotobacter vinelandii apoflavodoxin (apoFD) (un)folding are revealed by time-resolved fluorescence anisotropy using wild-type protein and variants lacking one or two of apoFD's three tryptophans. ApoFD equilibrium (un)folding by guanidine hydrochloride follows a three-state model: native ↔ unfolded ↔ intermediate. In native protein, W128 is a sink for Förster resonance energy transfer (FRET). Consequently, unidirectional FRET with a 50-ps transfer correlation time occurs from W167 to W128. FRET from W74 to W167 is much slower (6.9 ns). In the intermediate, W128 and W167 have native-like geometry because the 50-ps transfer time is observed. However, non-native structure exists between W74 and W167 because instead of 6.9 ns the transfer correlation time is 2.0 ns. In unfolded apoFD this 2.0-ns transfer correlation time is also detected. This decrease in transfer correlation time is a result of W74 and W167 becoming solvent accessible and randomly oriented toward one another. Apparently W74 and W167 are near-natively separated in the folding intermediate and in unfolded apoFD. Both tryptophans may actually be slightly closer in space than in the native state, even though apoFD's radius increases substantially upon unfolding. In unfolded apoFD the 50-ps transfer time observed for native and intermediate folding states becomes 200 ps as W128 and W167 are marginally further separated than in the native state. Apparently, apoFD's unfolded state is not a featureless statistical coil but contains well-defined substructures. The approach presented is a powerful tool to study protein folding.     

Folding Factors and Partners for the Intrinsically Disordered Protein Micro-Exon Gene 14 (MEG-14)

The micro-exon genes (MEG) of Schistosoma mansoni, a parasite responsible for the second most widely spread tropical disease, code for small secreted proteins with sequences unique to the Schistosoma genera. Bioinformatics analyses suggest the soluble domain of the MEG-14 protein will be largely disordered, and using synchrotron radiation circular dichroism spectroscopy, its secondary structure was shown to be essentially completely unfolded in aqueous solution. It does, however, show a strong propensity to fold into more ordered structures under a wide range of conditions. Partial folding was produced by increasing temperature (in a reversible process), contrary to the behavior of most soluble proteins. Furthermore, significant folding was observed in the presence of negatively charged lipids and detergents, but not in zwitterionic or neutral lipids or detergents. Absorption onto a surface followed by dehydration stimulated it to fold into a helical structure, as it did when the aqueous solution was replaced by nonaqueous solvents. Hydration of the dehydrated folded protein was accompanied by complete unfolding. These results support the identification of MEG-14 as a classic intrinsically disordered protein, and open the possibility of its interaction/folding with different partners and factors being related to multifunctional roles and states within the host.     

Characterization of the unfolded state of bovine alpha-lactalbumin and comparison with unfolded states of homologous proteins

The unfolded states of three homologous proteins with a very similar fold have been investigated by heteronuclear NMR spectroscopy. Secondary structure propensities as derived from interpretation of chemical shifts and motional restrictions as evidenced by heteronuclear (15)N relaxation rates have been analyzed in the reduced unfolded states of hen lysozyme and the calcium-binding proteins bovine alpha-lactalbumin and human alpha-lactalbumin. For all three proteins, significant deviations from random-coil predictions can be identified; in addition, the unfolded states also differ from each other, despite the fact that they possess very similar structures in their native states. Deviations from random-coil motional properties are observed in the alpha- and the beta-domain in bovine alpha-lactalbumin and lysozyme, while only regions within the alpha-domain deviate in human alpha-lactalbumin. The motional restrictions and residual secondary structure are determined both by the amino acid sequence of the protein and by residual long-range interactions. Even a conservative single point mutation from I to L in a highly conserved region between the two alpha-lactalbumins results in considerable differences in the motional properties. Given the differences in oxidative folding between hen lysozyme and alpha-lactalbumin, the results obtained on the unfolded states suggest that residual long-range interactions, i.e., those between the alpha- and the beta-domain of lysozyme, may act as nucleation sites for protein folding, while this property of residual structure is replaced by the calcium-binding site between the domains in alpha-lactalbumin.     

On the Prediction of Folding Nuclei in Globular Proteins

An approach to predicting folding nuclei in globular proteins with known three-dimensional structures is proposed. This approach is based on the pinpointing of the lowest saddle points on the barrier between the unfolded state and native structure on the free-energy landscape of a protein chain; the proposed technique uses the dynamic programming method. A comparison of calculation results with experimental data on the folding nuclei of 21 proteins shows that the model provides good Φ value predictions for protein structures determined by X-ray analysis and, less successfully, in structures determined by nuclear magnetic resonance. Consideration of the whole ensemble of transition states provides a better prediction of folding nuclei than consideration of only transition states with lowest free energies. In addition, we predict the location of folding nuclei in three-dimensional structures of some proteins whose folding kinetics is being studied, but there is no experimental evidence concerning their folding nuclei.     

The Folding Trajectory of RNase H Is Dominated by Its Topology and Not Local Stability: A Protein Engineering Study of Variants that Fold via Two-State and Three-State Mechanisms

Proteins can sample a variety of partially folded conformations during the transition between the unfolded and native states. Some proteins never significantly populate these high-energy states and fold by an apparently two-state process. However, many proteins populate detectable, partially folded forms during the folding process. The role of such intermediates is a matter of considerable debate. A single amino acid change can convert Escherichia coli ribonuclease H from a three-state folder that populates a kinetic intermediate to one that folds in an apparent two-state fashion. We have compared the folding trajectories of the three-state RNase H and the two-state RNase H, proteins with the same native-state topology but altered regional stability, using a protein engineering approach. Our data suggest that both versions of RNase H fold through a similar trajectory with similar high-energy conformations. Mutations in the core and the periphery of the protein affect similar aspects of folding for both variants, suggesting a common trajectory with folding of the core region followed by the folding of the periphery. Our results suggest that formation of specific partially folded conformations may be a general feature of protein folding that can promote, rather than hinder, efficient folding.