Helix stability and hydrophobicity in the folding mechanism of the bacterial immunity protein Im9

Recent models suggest that the mechanism of protein folding is determined by the balance between the stability of secondary structural elements and the hydrophobicity of the sequence. Here we determine the role of these factors in the folding kinetics of Im9* by altering the secondary structure propensity or hydrophobicity of helices I, II or IV by the substitution of residues at solvent exposed sites. The folding kinetics of each variant were measured at pH 7.0 and 10 degrees C, under which conditions wild-type Im9* folds with two-state kinetics. We show that increasing the helicity of these sequences in regions known to be structured in the folding intermediate of Im7*, switches the folding of Im9* from a two- to three-state mechanism. By contrast, increasing the hydrophobicity of helices I or IV has no effect on the kinetic folding mechanism. Interestingly, however, increasing the hydrophobicity of solvent-exposed residues in helix II stabilizes the folding intermediate and the rate-limiting transition state, consistent with the view that this helix makes significant non-native interactions during folding. The results highlight the generic importance of intermediates in folding and show that such species can be populated by increasing helical propensity or by stabilizing inter-helix contacts through non-native interactions.     

Characterisation of transition state structures for protein folding using 'high', 'medium' and 'low' {Phi}-values

It has been suggested that Phi-values, which allow structural information about transition states (TSs) for protein folding to be obtained, are most reliably interpreted when divided into three classes (high, medium and low). High Phi-values indicate almost completely folded regions in the TS, intermediate Phi-values regions with a detectable amount of structure and low Phi-values indicate mostly unstructured regions. To explore the extent to which this classification can be used to characterise in detail the structure of TSs for protein folding, we used Phi-values divided into these classes as restraints in molecular dynamics simulations. This type of procedure is related to that used in NMR spectroscopy to define the structure of native proteins from the measurement of inter-proton distances derived from nuclear Overhauser effects. We illustrate this approach by determining the TS ensembles of five proteins and by showing that the results are similar to those obtained by using as restraints the actual numerical Phi-values measured experimentally. Our results indicate that the simultaneous consideration of a set of low-resolution Phi-values can provide sufficient information for characterising the architecture of a TS for folding of a protein.     

A simple approach for protein structure discrimination based on the network pattern of conserved hydrophobic residues

Evolutionarily conserved hydrophobic residues at the core of protein structures are generally assumed to play a structural role in protein folding and stability. Recent studies have implicated that their importance to protein structures is uneven, with a few of them being crucial and the rest of them being secondary. In this work, we explored the possibility of employing this feature of native structures for discriminating non-native structures from native ones. First, we developed a network tool to quantitatively measure the structural contributions of individual amino acid residues. We systematically applied this method to diverse fold-type sets of native proteins. It was confirmed that this method could grasp the essential structural features of native proteins. Next, we applied it to a number of decoy sets of proteins. The results indicate that such an approach indeed identified non-native structures in most test cases. This finding should be of help for the investigation of the fundamental problem of protein structure prediction.     

Antagonism, non-native interactions and non-two-state folding in S6 revealed by double-mutant cycle analysis

When folding to the native state N in the presence of salt, the apparent two-state folder S6 transiently forms a transient off-pathway state C with substantial secondary and tertiary structure. Fifteen double mutant cycles were analysed to compare side-chain interaction energies DeltaDeltaG(int) in C, N and TS (the transition state between N and the denatured state). The kinetic signatures of these destabilizing mutants suggest folding scenarios involving unfolding intermediates and even alternative unfolding pathways. However, restricting the kinetic data to linear parts of the chevron plot allows reliable extrapolation to zero molar denaturant of rate constants of folding, unfolding and misfolding. Side-chain interactions appear to contribute to the stability of C, but in a substantially non-native environment, as shown by changes in the sign of DeltaDeltaG(int) between C and N. Remarkably, there appear to be significant (0.7-2 kcal/mol) antagonistic interactions between the two residues Leu30 and Leu75 in N and TS, which may be linked to subtle structural changes seen in the crystal structures of the mutants. A small number of overlapping residues are involved in these kinds of antagonistic interactions in N, TS and C, suggesting that repulsive interactions are coded into the protein topology whether the protein folds or misfolds. Destabilizing double mutants indicate that apparent two-state folders can be induced to behave in more complex ways provided that the native state is suitably destabilized.     

Fluorescence resonance energy transfer analysis of the folding pathway of Engrailed Homeodomain

The Engrailed Homeodomain folds on the microsecond time scale via an intermediate that is experimentally well characterised using structural Engrailed-Homeodomain mimics. Here, we analysed directly the changes in distance between key residues during the kinetics of unfolding and at equilibrium using fluorescence resonance energy transfer (FRET). Trp was the donor and 5-(((acetylamino)ethyl)amino) naphthalene-1-sulphate, the acceptor, substituted in positions that caused little change in stability. Distances calculated for the native state were in good agreement with those derived from the NMR structure. The distances between the N- and C-termini of Helix I and of Helix III increased, then decreased and finally increased again with increasing GdmCl concentration on equilibrium denaturation. This behaviour implied that there was a folding intermediate on the folding pathway and that this intermediate was populated at low concentrations of GdmCl concentration ( approximately 1 M). We analysed the changes in distance during temperature-jump relaxation kinetics, using a qualitative and very conservative procedure that drew conclusions only when changes in fluorescence of mutants containing either the donor or the acceptor alone would not obscure the change in the FRET signal when both donor and acceptor were present. The distance changes obtained under equilibrium and kinetic measurements were self-consistent and also consistent with the known high-resolution structures of the mimics of the folding intermediates. We showed that for analysing distances in disordered ensembles, it is important to use FRET probes with a critical distance close to the average separation in the ensemble. Otherwise, average distances could be over or underestimated.     

Yeast hexokinase isoenzyme ScHxk2: stability of a two-domain protein with discontinuous domains

The hexokinase isoenzyme 2 of Saccharomyces cerevisiae (ScHxk2) represents an archetype of a two-domain protein with the active site located in a cleft between the two domains. Binding of the substrate glucose results in a rigid body movement of the two domains leading to a cleft closure of the active site. Both domains of this enzyme are composed of discontinuous peptide sequences. This structural feature is reflected in the stability and folding of the ScHxk2 protein. Structural transitions induced by urea treatment resulted in the population of a thermodynamically stable folding intermediate, which, however, does not correspond to a molecule with one domain folded and the other unfolded. As demonstrated by different spectroscopic techniques, both domains are structurally affected by the partial denaturation. The intermediate possesses only 40% of the native secondary structural content and a substantial increase in the Stokes radius as judged by circular dichroism and dynamic light scattering analyses. One-dimensional 1H NMR data prove that all tryptophan residues are in a non-native environment in the intermediate, indicating substantial changes in the tertiary structure. Still, the intermediate possesses quite a high stability for a transition intermediate of about ΔG = -22 kJ mol?1.     

Protein Folding and Misfolding on Surfaces

Protein folding, misfolding and aggregation, as well as the way misfolded and aggregated proteins affects cell viability are emerging as key themes in molecular and structural biology and in molecular medicine. Recent advances in the knowledge of the biophysical basis of protein folding have led to propose the energy landscape theory which provides a consistent framework to better understand how a protein folds rapidly and efficiently to the compact, biologically active structure. The increased knowledge on protein folding has highlighted its strict relation to protein misfolding and aggregation, either process being in close competition with the other, both relying on the same physicochemical basis. The theory has also provided information to better understand the structural and environmental factors affecting protein folding resulting in protein misfolding and aggregation into ordered or disordered polymeric assemblies. Among these, particular importance is given to the effects of surfaces. The latter, in some cases make possible rapid and efficient protein folding but most often recruit proteins/peptides increasing their local concentration thus favouring misfolding and accelerating the rate of nucleation. It is also emerging that surfaces can modify the path of protein misfolding and aggregation generating oligomers and polymers structurally different from those arising in the bulk solution and endowed with different physical properties and cytotoxicities.     

RNA folding pathways and the self-assembly of ribosomes

Many RNAs do not directly code proteins but are nonetheless indispensable to cellular function. These strands fold into intricate three-dimensional shapes that are essential structures in protein synthesis, splicing, and many other processes of gene regulation and expression. A variety of biophysical and biochemical methods are now showing, in real time, how ribosomal subunits and other ribonucleoprotein complexes assemble from their molecular components. Footprinting methods are particularly useful for studying the folding of long RNAs: they provide quantitative information about the conformational state of each residue and require little material. Data from footprinting complement the global information available from small-angle X-ray scattering or cryo-electron microscopy, as well as the dynamic information derived from single-molecule F?rster resonance energy transfer (FRET) and NMR methods. In this Account, I discuss how we have used hydroxyl radical footprinting and other experimental methods to study pathways of RNA folding and 30S ribosome assembly. Hydroxyl radical footprinting probes the solvent accessibility of the RNA backbone at each residue in as little as 10 ms, providing detailed views of RNA folding pathways in real time. In conjunction with other methods such as solution scattering and single-molecule FRET, time-resolved footprinting of ribozymes showed that stable domains of RNA tertiary structure fold in less than 1 s. However, the free energy landscapes for RNA folding are rugged, and individual molecules kinetically partition into folding pathways that lead through metastable intermediates, stalling the folding or assembly process. Time-resolved footprinting was used to follow the formation of tertiary structure and protein interactions in the 16S ribosomal RNA (rRNA) during the assembly of 30S ribosomes. As previously observed in much simpler ribozymes, assembly occurs in stages, with individual molecules taking different routes to the final complex. Interactions occur concurrently in all domains of the 16S rRNA, and multistage protection of binding sites of individual proteins suggests that initial encounter complexes between the rRNA and ribosomal proteins are remodeled during assembly. Equilibrium footprinting experiments showed that one primary binding protein was sufficient to stabilize the tertiary structure of the entire 16S 5'-domain. The rich detail available from the footprinting data showed that the secondary assembly protein S16 suppresses non-native structures in the 16S 5'-domain. In doing so, S16 enables a conformational switch distant from its own binding site, which may play a role in establishing interactions with other domains of the 30S subunit. Together, the footprinting results show how protein-induced changes in RNA structure are communicated over long distances, ensuring cooperative assembly of even very large RNA-protein complexes such as the ribosome.     

Contrasting disease and nondisease protein aggregation by molecular simulation

[Figurre: see text]. Protein aggregation can be defined as the sacrifice of stabilizing intrachain contacts of the functional state that are replaced with interchain contacts to form non-functional states. The resulting aggregate morphologies range from amorphous structures without long-range order typical of nondisease proteins involved in inclusion bodies to highly structured fibril assemblies typical of amyloid disease proteins. In this Account, we describe the development and application of computational models for the investigation of nondisease and disease protein aggregation as illustrated for the proteins L and G and the Alzheimer's Abeta systems. In each case, we validate the models against relevant experimental observables and then expand on the experimental window to better elucidate the link between molecular properties and aggregation outcomes. Our studies show that each class of protein exhibits distinct aggregation mechanisms that are dependent on protein sequence, protein concentration, and solution conditions. Nondisease proteins can have native structural elements in the denatured state ensemble or rapidly form early folding intermediates, which offers avenues of protection against aggregation even at relatively high concentrations. The possibility that early folding intermediates may be evolutionarily selected for their protective role against unwanted aggregation could be a useful strategy for reengineering sequences to slow aggregation and increase folding yield in industrial protein production. The observed oligomeric aggregates that we see for nondisease proteins L and G may represent the nuclei for larger aggregates, not just for large amorphous inclusion bodies, but potentially as the seeds of ordered fibrillar aggregates, since most nondisease proteins can form amyloid fibrils under conditions that destabilize the native state. By contrast, amyloidogenic protein sequences such as Abeta 1-40,42 and the familial Alzheimer's disease (FAD) mutants favor aggregation into ordered fibrils once the free-energy barrier for forming a critical nucleus is crossed. However, the structural characteristics and oligomer size of the soluble nucleation species have yet to be determined experimentally for any disease peptide sequence, and the molecular mechanism of polymerization that eventually delineates a mature fibril is unknown. This is in part due to the limited experimental access to very low peptide concentrations that are required to characterize these early aggregation events, providing an opportunity for theoretical studies to bridge the gap between the monomer and fibril end points and to develop testable hypotheses. Our model shows that Abeta 1-40 requires as few as 6-10 monomer chains (depending on sequence) to begin manifesting the cross-beta order that is a signature of formation of amyloid filaments or fibrils assessed in dye-binding kinetic assays. The richness of the oligomeric structures and viable filament and fibril polymorphs that we observe may offer structural clues to disease virulence variations that are seen for the WT and hereditary mutants.     

Sequence and structural analysis of two designed proteins with 88% identity adopting different folds

Protein folding is a natural phenomenon by which a sequence of amino acids folds into a unique functional three-dimensional structure. Although the sequence code that governs folding remains a mystery, one can identify key inter-residue contacts responsible for a given topology. In nature, there are many pairs of proteins of a given length that share little or no sequence identity. Similarly, there are many proteins that share a common topology but lack significant evidence of homology. In order to tackle this problem, protein engineering studies have been used to determine the minimal number of amino acid residues that codes for a particular fold. In recent years, the coupling of theoretical models and experiments in the study of protein folding has resulted in providing some fruitful clues. He et al. have designed two proteins with 88% sequence identity, which adopt different folds and functions. In this work, we have systematically analysed these two proteins by performing pentapeptide search, secondary structure predictions, variation in inter-residue interactions and residue-residue pair preferences, surrounding hydrophobicity computations, conformational switching and energy computations. We conclude that the local secondary structural preference of the two designed proteins at the Nand C-terminal ends to adopt either coil or strand conformation may be a crucial factor in adopting the different folds. Early on during the process of folding, both proteins may choose different energetically favourable pathways to attain the different folds.     

Probing invisible, low-populated States of protein molecules by relaxation dispersion NMR spectroscopy: an application to protein folding

Biological function depends on molecular dynamics that lead to excursions from highly populated ground states to much less populated excited states. The low populations and the transient formation of such excited states render them invisible to the conventional methods of structural biology. Thus, while detailed pictures of ground-state structures of biomolecules have emerged over the years, largely through X-ray diffraction and solution nuclear magnetic resonance (NMR) spectroscopy studies, much less structural data has been accumulated on the conformational properties of the invisible excited states that are necessary to fully explain function. NMR spectroscopy is a powerful tool for studying conformational dynamics because it is sensitive to dynamics over a wide range of time scales, extending from picoseconds to seconds and because information is, in principle, available at nearly every position in the molecule. Here an NMR method for quantifying millisecond time scale dynamics that involve transitions between different molecular conformations is described. The basic experimental approach, termed relaxation dispersion NMR spectroscopy, is outlined to provide the reader with an intuitive feel for the technology. A variety of different experiments that probe conformational exchange at different sites in proteins are described, including a brief summary of data-fitting procedures to extract both the kinetic and thermodynamic properties of the exchange process and the structural features of the invisible excited states along the exchange pathway. It is shown that the methodology facilitates detection of intermediates and other excited states that are populated at low levels, 0.5% or higher, that cannot be observed directly in spectra, so long as they exchange with the observable ground state of the protein on the millisecond time scale. The power of the methodology is illustrated by a detailed application to the study of protein folding of the small modular SH3 domain. The kinetics and thermodynamics that describe the folding of this domain have been characterized through the effects of temperature, pressure, side-chain deuteration, and mutation, and the structural features of a low-populated folding intermediate have been assessed. Despite the fact that many previous studies have shown that SH3 domains fold via a two-state mechanism, the NMR methods presented unequivocally establish the presence of an on-pathway folding intermediate. The unique capabilities of NMR relaxation dispersion follow from the fact that large numbers of residues can be probed individually in a single experiment. By contrast, many other forms of spectroscopy monitor properties that are averaged over all residues in the molecule or that make use of only one or two reporters. The NMR methodology is not limited to protein folding, and applications to enzymatic catalysis, binding, and molecular recognition are beginning to emerge.     

Close-range electrostatic interactions in proteins

Two types of noncovalent bonding interactions are present in protein structures, specific and nonspecific. Nonspecific interactions are mostly hydrophobic and van der Waals. Specific interactions are largely electrostatic. While the hydrophobic effect is the major driving force in protein folding, electrostatic interactions are important in protein folding, stability, flexibility, and function. Here we review the role of close-range electrostatic interactions (salt bridges) and their networks in proteins. Salt bridges are formed by spatially proximal pairs of oppositely charged residues in native protein structures. Often salt-bridging residues are also close in the protein sequence and fall in the same secondary structural element, building block, autonomous folding unit, domain, or subunit, consistent with the hierarchical model for protein folding. Recent evidence also suggests that charged and polar residues in largely hydrophobic interfaces may act as hot spots for binding. Salt bridges are rarely found across protein parts which are joined by flexible hinges, a fact suggesting that salt bridges constrain flexibility and motion. While conventional chemical intuition expects that salt bridges contribute favorably to protein stability, recent computational and experimental evidence shows that salt bridges can be stabilizing or destabilizing. Due to systemic protein flexibility, reflected in small-scale side-chain and backbone atom motions, salt bridges and their stabilities fluctuate in proteins. At the same time, genome-wide, amino acid sequence composition, structural, and thermodynamic comparisons of thermophilic and mesophilic proteins indicate that specific interactions, such as salt bridges, may contribute significantly towards the thermophilic-mesophilic protein stability differential.     

Kinetic folding mechanism of PDZ2 from PTP-BL

PDZ domains represent a large family of protein-interaction modules associated with a variety of unrelated proteins with different functions. We report a complete characterization of the kinetic folding mechanism of a fluorescent variant of PDZ2 from PTP-BL, investigated under a variety of different experimental conditions. For this purpose, we engineered a fluorescent variant of this protein Y43W (called pseudo-wild-type, pWT43). The results suggest the presence of a high-energy intermediate in the folding of PDZ2, as revealed by a pronounced non-linear dependence of the unfolding rate constant on denaturant concentration. Such an intermediate may or may not be detectable depending on the experimental conditions, giving rise to apparent two-state folding under stabilizing conditions (e.g. in the presence of sodium sulfate). Interestingly, even under these conditions, three-state folding can be restored by selectively destabilizing the native-like rate-limiting barrier by one specific mutation (V44A). Finally, we show that data taken on pWT43 under different experimental conditions (e.g. different pH values from 2.1 to 8.0 or in the presence of a stabilizing salt) and also data on a site-directed conservative mutant can be rationalized in terms of a simple reaction scheme involving a single set of intermediates and transition states.     

The Role of Disordered Ribosomal Protein Extensions in the Early Steps of Eubacterial 50 S Ribosomal Subunit Assembly

Although during the past decade research has shown the functional importance of disorder in proteins, many of the structural and dynamics properties of intrinsically unstructured proteins (IUPs) remain to be elucidated. This review is focused on the role of the extensions of the ribosomal proteins in the early steps of the assembly of the eubacterial 50 S subunit. The recent crystallographic structures of the ribosomal particles have revealed the picture of a complex assembly pathway that condenses the rRNA and the ribosomal proteins into active ribosomes. However, little is know about the molecular mechanisms of this process. It is thought that the long basic r-protein extensions that penetrate deeply into the subunit cores play a key role through disorder-order transitions and/or co-folding mechanisms. A current view is that such structural transitions may facilitate the proper rRNA folding. In this paper, the structures of the proteins L3, L4, L13, L20, L22 and L24 that have been experimentally found to be essential for the first steps of ribosome assembly have been compared. On the basis of their structural and dynamics properties, three categories of extensions have been identified. Each of them seems to play a distinct function. Among them, only the coil-helix transition that occurs in a phylogenetically conserved cluster of basic residues of the L20 extension appears to be strictly required for the large subunit assembly in eubacteria. The role of α helix-coil transitions in 23 S RNA folding is discussed in the light of the calcium binding protein calmodulin that shares many structural and dynamics properties with L20.     

Coarse-Grained and Atomistic MD Simulations of RNA and DNA Folding

Although the main features of the protein folding problem are coming into clearer focus, the microscopic viewpoint of nucleic acid folding mechanisms is only just beginning to be addressed. Experiments, theory, and simulations are pointing to complex thermodynamic and kinetic mechanisms. As is the case for proteins, molecular dynamics (MD) simulations continue to be indispensable tools for providing a molecular basis for nucleic acid folding mechanisms. In this review, we provide an overview of biomolecular folding mechanisms focusing on nucleic acids. We outline the important interactions that are likely to be the main determinants of nucleic acid folding energy landscapes. We discuss recent MD simulation studies of empirical force field and Go-type MD simulations of RNA and DNA folding mechanisms to outline recent successes and the theoretical and computational challenges that lie ahead.     

Intermediates in the protein folding process: a computational model

The paper presents a model for simulating the protein folding process in silico. The two-step model (which consists of the early stage-ES and the late stage-LS) is verified using two proteins, one of which is treated (according to experimental observations) as the early stage and the second as an example of the LS step. The early stage is based solely on backbone structural preferences, while the LS model takes into account the water environment, treated as an external hydrophobic force field and represented by a 3D Gauss function. The characteristics of 1ZTR (the ES intermediate, as compared with 1ENH, which is the LS intermediate) confirm the link between the gradual disappearance of ES characteristics in LS structural forms and the simultaneous emergence of LS properties in the 1ENH protein. Positive verification of ES and LS characteristics in these two proteins (1ZTR and 1ENH respectively) suggest potential applicability of the presented model to in silico protein folding simulations.     

Is Protein Folding a Thermodynamically Unfavorable,Active, Energy-Dependent Process?

The prevailing current view of protein folding is the thermodynamic hypothesis, under which the native folded conformation of a protein corresponds to the global minimum of Gibbs free energy G. We question this concept and show that the empirical evidence behind the thermodynamic hypothesis of folding is far from strong. Furthermore, physical theory-based approaches to the prediction of protein folds and their folding pathways so far have invariably failed except for some very small proteins, despite decades of intensive theory development and the enormous increase of computer power. The recent spectacular successes in protein structure prediction owe to evolutionary modeling of amino acid sequence substitutions enhanced by deep learning methods, but even these breakthroughs provide no information on the protein folding mechanisms and pathways. We discuss an alternative view of protein folding, under which the native state of most proteins does not occupy the global free energy minimum, but rather, a local minimum on a fluctuating free energy landscape. We further argue that ΔG of folding is likely to be positive for the majority of proteins, which therefore fold into their native conformations only through interactions with the energy-dependent molecular machinery of living cells, in particular, the translation system and chaperones. Accordingly, protein folding should be modeled as it occurs in vivo, that is, as a non-equilibrium, active, energy-dependent process.