Comparison of the transition state ensembles for folding of Im7 and Im9 determined using all-atom molecular dynamics simulations with ϕ value restraints

From covalent transition states in chemistry to noncovalent in biology: from β- to Φ-value analysis of protein folding

Solving the mechanism of a chemical reaction requires determining the structures of all the ground states on the pathway and the elusive transition states linking them. 2024 is the centenary of Brønsted's landmark paper that introduced the β-value and structure-activity studies as the only experimental means to infer the structures of transition states. It involves making systematic small changes in the covalent structure of the reactants and analysing changes in activation and equilibrium-free energies. Protein engineering was introduced for an analogous procedure, Φ-value analysis, to analyse the noncovalent interactions in proteins central to biological chemistry. The methodology was developed first by analysing noncovalent interactions in transition states in enzyme catalysis. The mature procedure was then applied to study transition states in the pathway of protein folding - 'part (b) of the protein folding problem'. This review describes the development of Φ-value analysis of transition states and compares and contrasts the interpretation of β- and Φ-values and their limitations. Φ-analysis afforded the first description of transition states in protein folding at the level of individual residues. It revealed the nucleation-condensation folding mechanism of protein domains with the transition state as an expanded, distorted native structure, containing little fully formed secondary structure but many weak tertiary interactions. A spectrum of transition states with various degrees of structural polarisation was then uncovered that spanned from nucleation-condensation to the framework mechanism of fully formed secondary structure. Φ-analysis revealed how movement of the expanded transition state on an energy landscape accommodates the transition from framework to nucleation-condensation mechanisms with a malleability of structure as a unifying feature of folding mechanisms. Such movement follows the rubric of analysis of classical covalent chemical mechanisms that began with Brønsted. Φ-values are used to benchmark computer simulation, and Φ and simulation combine to describe folding pathways at atomic resolution.

Some useful ideas for multistate protein design: Effect of amino acid substitutions on the multistate proteins stability and the rate of protein structure formation

The design of new protein variants is usually confined to slightly “fixing” an already existing protein, adapting it to certain conditions or to a new substrate. This is relatively easy to do if the fragment of the protein to be affected, such as the active site of the protein, is known. But what if you need to “fix” the stability of a protein or the rate of its native or intermediate state formation? Having studied a large number of protein mutant forms, we have established the effect of various amino acid substitutions on the energy landscape of the protein. As a result, we have revealed a number of patterns to help researchers identify amino acid residues that determine the folding rate and the stability of globular proteins states and design a mutant form of a protein with desired properties.

Nonnative Energetic Frustrations in Protein Folding at Residual Level: A Simulation Study of Homologous Immunoglobulin-like β-Sandwich Proteins

Nonnative interactions cause energetic frustrations in protein folding and were found to dominate key events in folding intermediates. However, systematically characterizing energetic frustrations that are caused by nonnative intra-residue interactions at residual resolution is still lacking. Recently, we studied the folding of a set of homologous all-α proteins and found that nonnative-contact-based energetic frustrations are highly correlated to topology of the protein native-contact network. Here, we studied the folding of nine homologous immunoglobulin-like (Ig-like) β-sandwich proteins, and examined nonnative-contact-based energetic frustrations Gō-like model. Our calculations showed that nonnative-interaction-based energetic frustrations in β-sandwich proteins are much more complicated than those in all-α proteins, and they exhibit highly heterogeneous effects on the folding of secondary structures. Further, the nonnative interactions introduced distinct correlations in the folding of different folding-patches of β-sandwich proteins. Taken together, a strong interplay might exist between nonnative-interaction energetic frustrations and the protein native-contact networks, which ensures that β-sandwich domains adopt a common folding mechanism.

Molecular recognition and packing frustration in a helical protein

Biomolecular recognition entails attractive forces for the functional native states and discrimination against potential nonnative interactions that favor alternate stable configurations. The challenge posed by the competition of nonnative stabilization against native-centric forces is conceptualized as frustration. Experiment indicates that frustration is often minimal in evolved biological systems although nonnative possibilities are intuitively abundant. Much of the physical basis of minimal frustration in protein folding thus remains to be elucidated. Here we make progress by studying the colicin immunity protein Im9. To assess the energetic favorability of nonnative versus native interactions, we compute free energies of association of various combinations of the four helices in Im9 (referred to as H1, H2, H3, and H4) by extensive explicit-water molecular dynamics simulations (total simulated time > 300 μs), focusing primarily on the pairs with the largest native contact surfaces, H1-H2 and H1-H4. Frustration is detected in H1-H2 packing in that a nonnative packing orientation is significantly stabilized relative to native, whereas such a prominent nonnative effect is not observed for H1-H4 packing. However, in contrast to the favored nonnative H1-H2 packing in isolation, the native H1-H2 packing orientation is stabilized by H3 and loop residues surrounding H4. Taken together, these results showcase the contextual nature of molecular recognition, and suggest further that nonnative effects in H1-H2 packing may be largely avoided by the experimentally inferred Im9 folding transition state with native packing most developed at the H1-H4 rather than the H1-H2 interface.

Biomolecules need to recognize one another with high specificity: promoting “native” functional intermolecular binding events while avoiding detrimental “nonnative” bound configurations; i.e., “frustration”—the tendency for nonnative interactions—has to be minimized. Folding of globular proteins entails a similar discrimination. To gain physical insight, we computed the binding affinities of helical structures of the protein Im9 in various native or nonnative configurations by atomic simulations, discovering that partial packing of the Im9 core is frustrated. This frustration is overcome when the entire core of the protein is assembled, consistent with experiment indicating no significant kinetic trapping in Im9 folding. Our systematic analysis thus reveals a subtle, contextual aspect of biomolecular recognition and provides a general approach to characterize folding frustration.

Folding Mechanism of Proteins Im7 and Im9: Insight from All-Atom Simulations in Implicit and Explicit Solvent

Im7 and Im9 are evolutionary related proteins with almost identical native structures. In spite of their structural similarity, experiments show that Im7 folds through a long-lived on-pathway intermediate, while Im9 folds according to two-state kinetics. In this work, we use a recently developed enhanced path sampling method to generate many folding trajectories for these proteins, using realistic atomistic force fields, in both implicit and explicit solvent. Overall, our results are in good agreement with the experimental ϕ values and with the result of ϕ-value-restrained molecular dynamics (MD) simulations. However, our implicit solvent simulations fail to predict a qualitative difference in the folding pathways of Im7 and Im9. In contrast, our simulations in explicit solvent correctly reproduce the fact that only protein Im7 folds through a on-pathway intermediate. By analyzing our atomistic trajectories, we provide a physical picture which explains the observed difference in the folding kinetics of these chains.

Native contact density and nonnative hydrophobic effects in the folding of bacterial immunity proteins

The bacterial colicin-immunity proteins Im7 and Im9 fold by different mechanisms. Experimentally, at pH 7.0 and 10°C, Im7 folds in a three-state manner via an intermediate but Im9 folding is two-state-like. Accordingly, Im7 exhibits a chevron rollover, whereas the chevron arm for Im9 folding is linear. Here we address the biophysical basis of their different behaviors by using native-centric models with and without additional transferrable, sequence-dependent energies. The Im7 chevron rollover is not captured by either a pure native-centric model or a model augmented by nonnative hydrophobic interactions with a uniform strength irrespective of residue type. By contrast, a more realistic nonnative interaction scheme that accounts for the difference in hydrophobicity among residues leads simultaneously to a chevron rollover for Im7 and an essentially linear folding chevron arm for Im9. Hydrophobic residues identified by published experiments to be involved in nonnative interactions during Im7 folding are found to participate in the strongest nonnative contacts in this model. Thus our observations support the experimental perspective that the Im7 folding intermediate is largely underpinned by nonnative interactions involving large hydrophobics. Our simulation suggests further that nonnative effects in Im7 are facilitated by a lower local native contact density relative to that of Im9. In a one-dimensional diffusion picture of Im7 folding with a coordinate- and stability-dependent diffusion coefficient, a significant chevron rollover is consistent with a diffusion coefficient that depends strongly on native stability at the conformational position of the folding intermediate.

In order to fold correctly, a globular protein must avoid being trapped in wrong, i.e., nonnative conformations. Thus a biophysical account of how attractive nonnative interactions are bypassed by some amino acid sequences but not others is key to deciphering protein structure and function. We examine two closely related bacterial immunity proteins, Im7 and Im9, that are experimentally known to fold very differently: Whereas Im9 folds directly, Im7 folds through a mispacked conformational intermediate. A simple model we developed accounts for their intriguingly different folding kinetics in terms of a balance between the density of native-promoting contacts and the hydrophobicity of local amino acid sequences. This emergent principle is extensible to other biomolecular recognition processes.

Energetic frustrations in protein folding at residue resolution: a homologous simulation study of Im9 proteins

Energetic frustration is becoming an important topic for understanding the mechanisms of protein folding, which is a long-standing big biological problem usually investigated by the free energy landscape theory. Despite the significant advances in probing the effects of folding frustrations on the overall features of protein folding pathways and folding intermediates, detailed characterizations of folding frustrations at an atomic or residue level are still lacking. In addition, how and to what extent folding frustrations interact with protein topology in determining folding mechanisms remains unclear. In this paper, we tried to understand energetic frustrations in the context of protein topology structures or native-contact networks by comparing the energetic frustrations of five homologous Im9 alpha-helix proteins that share very similar topology structures but have a single hydrophilic-to-hydrophobic mutual mutation. The folding simulations were performed using a coarse-grained Gō-like model, while non-native hydrophobic interactions were introduced as energetic frustrations using a Lennard-Jones potential function. Energetic frustrations were then examined at residue level based on φ-value analyses of the transition state ensemble structures and mapped back to native-contact networks. Our calculations show that energetic frustrations have highly heterogeneous influences on the folding of the four helices of the examined structures depending on the local environment of the frustration centers. Also, the closer the introduced frustration is to the center of the native-contact network, the larger the changes in the protein folding. Our findings add a new dimension to the understanding of protein folding the topology determination in that energetic frustrations works closely with native-contact networks to affect the protein folding.

Assessing and refining molecular dynamics simulations of proteins with nuclear magnetic resonance data

The sophistication of the force fields, algorithms and hardware used for molecular dynamics (MD) simulations of proteins is continuously increasing. No matter how advanced the methodology, however, it is essential to evaluate the appropriateness of the structures sampled in a simulation by comparison with quantitative experimental data. Solution nuclear magnetic resonance (NMR) data are particularly useful for checking the quality of protein simulations, as they provide both structural and dynamic information on a variety of temporal and spatial scales. Here, various features and implications of using NMR data to validate and bias MD simulations are outlined, including an overview of the different types of NMR data that report directly on structural properties and of relevant simulation techniques. The focus throughout is on how to properly account for conformational averaging, particularly within the context of the assumptions inherent in the relationships that link NMR data to structural properties.

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.

Nuclease colicins and their immunity proteins

It is more than 80 years since Gratia first described 'a remarkable antagonism between two strains of Escherichia coli'. Shown subsequently to be due to the action of proteins (or peptides) produced by one bacterium to kill closely related species with which it might be cohabiting, such bacteriocins have since been shown to be commonplace in the internecine warfare between bacteria. Bacteriocins have been studied primarily from the twin perspectives of how they shape microbial communities and how they penetrate bacteria to kill them. Here, we review the modes of action of a family of bacteriocins that cleave nucleic acid substrates in E. coli, known collectively as nuclease colicins, and the specific immunity (inhibitor) proteins that colicin-producing organisms make in order to avoid committing suicide. In a process akin to targeting in mitochondria, nuclease colicins engage in a variety of cellular associations in order to translocate their cytotoxic domains through the cell envelope to the cytoplasm. As well as informing on the process itself, the study of nuclease colicin import has also illuminated functional aspects of the host proteins they parasitize. We also review recent studies where nuclease colicins and their immunity proteins have been used as model systems for addressing fundamental problems in protein folding and protein-protein interactions, areas of biophysics that are intimately linked to the role of colicins in bacterial competition and to the import process itself.

NMR characterisation of the relationship between frustration and the excited state of Im7

Previous work shows that Im9 folds in a two-state transition while its homologue Im7 folds in a three-state transition via an on-pathway kinetic intermediate state (KIS), with this difference being related to frustration in the structure of Im7. We have used NMR spectroscopy to study conformational dynamics connected to the frustration. A combination of equilibrium peptide N(1)H/N(2)H exchange, model-free analyses of backbone NH relaxation data and relaxation dispersion (RD)-NMR shows that the native state of Im7 is in equilibrium with an intermediate state that is lowly populated [equilibrium intermediate state (EIS)]. Comparison of kinetic and thermodynamic parameters describing the EIS native-state equilibrium obtained by RD-NMR with previously reported parameters describing the KIS native-state equilibrium obtained from stopped-flow fluorescence studies of refolding His-tagged Im7 shows that the KIS and the EIS are the same species. (15)N chemical shifts of the EIS obtained from the RD-NMR analysis show that residues forming helix III in the native state are unstructured in the EIS while other residues experiencing frustration in the native state are in structured regions of the EIS. We show that binding of Im7 and its L53A/I54A variant (which resembles the EIS as shown in previous work) to the cognate partner for Im7, the DNase domain of colicin E7, causes the dynamic processes associated with the frustration to be dampened.

Nonequilibrium NMR Methods for Monitoring Protein and RNA Folding

The review introduces to time-resolved NMR spectroscopic investigations of the kinetics of protein and RNA folding. The description of the experimental investigations is discussed in the context of possible kinetic folding pathways showing the extent of information that can be gained from the various kinetic experiments. The review introduces to four different methods to initiate folding reactions in connection with time-resolved NMR experiments and discusses examples of refolding of the model proteinα-lactalbumin and of bistable RNAs.

Simulation of fluorescence resonance energy transfer experiments: effect of the dyes on protein folding

Fluorescence resonance energy transfer is a powerful technique which is often used to probe the properties of proteins and complex macromolecules. The technique relies on relatively large fluorescent dyes which are engineered into the molecule of interest. In the case of small proteins, these dyes may affect the stability of the protein, and modify the folding kinetics and the folding mechanisms which are being probed. Here we use atomistic simulation to investigate the effect that commonly used fluorescent dyes have on the folding of a four-helix bundle protein. We show that, depending on where the dyes are attached, their effect on the kinetic and thermodynamic properties of the protein may be significant. We find that, while the overall folding mechanism is not affected by the dyes, they can destabilize, or even stabilize, intermediate states.

The nucleation mechanism of protein folding: a survey of computer simulation studies

The nucleation mechanism of protein folding, originally proposed by Baldwin in the early 1970s, was firstly observed by Shakhnovich and co-workers two decades later in the context of Monte Carlo simulations of a simple lattice model. At about the same time the extensive use of φ-value analysis provided the first experimental evidence that the folding of Chymotrypsin-inhibitor 2, a small single-domain protein, which folds with two-state kinetics, is also driven by a nucleation mechanism. Since then, the nucleation mechanism is generally considered the most common form of folding mechanism amongst two-state proteins. However, recent experimental data has put forward the idea that this may not necessarily be so, since the accuracy of the experimentally determined φ values, which are used to identify the critical (i.e. nucleating) residues, is typically poor. Here, we provide a survey of in silico results on the nucleation mechanism, ranging from simple lattice Monte Carlo to more sophisticated off-lattice molecular dynamics simulations, and discuss them in light of experimental data.

Probing the free energy landscape of the FBP28WW domain using multiple techniques

The free-energy landscape of a small protein, the FBP 28 WW domain, has been explored using molecular dynamics (MD) simulations with alternative descriptions of the molecule. The molecular models used range from coarse-grained to all-atom with either an implicit or explicit treatment of the solvent. Sampling of conformation space was performed using both conventional and temperature-replica exchange MD simulations. Experimental chemical shifts and NOEs were used to validate the simulations, and experimental phi values both for validation and as restraints. This combination of different approaches has provided insight into the free energy landscape and barriers encountered by the protein during folding and enabled the characterization of native, denatured and transition states which are compatible with the available experimental data. All the molecular models used stabilize well defined native and denatured basins; however, the degree of agreement with the available experimental data varies. While the most detailed, explicit solvent model predicts the data reasonably accurately, it does not fold despite a simulation time 10 times that of the experimental folding time. The less detailed models performed poorly relative to the explicit solvent model: an implicit solvent model stabilizes a ground state which differs from the experimental native state, and a structure-based model underestimates the size of the barrier between the two states. The use of experimental phi values both as restraints, and to extract structures from unfolding simulations, result in conformations which, although not necessarily true transition states, appear to share the geometrical characteristics of transition state structures. In addition to characterizing the native, transition and denatured states of this particular system in this work, the advantages and limitations of using varying levels of representation are discussed.

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.

The Effect of Increasing the Stability of Non-native Interactions on the Folding Landscape of the Bacterial Immunity Protein Im9

How stabilising non-native interactions influence protein folding energy landscapes is currently not well understood: such interactions could speed folding by reducing the conformational search to the native state, or could slow folding by increasing ruggedness. Here, we examine the influence of non-native interactions in the folding process of the bacterial immunity protein Im9, by exploiting our ability to manipulate the stability of the intermediate and rate-limiting transition state (TS) in the folding of this protein by minor alteration of its sequence or changes in solvent conditions. By analysing the properties of these species using Phi-value analysis, and exploration of the structural properties of the TS ensemble using molecular dynamics simulations, we demonstrate the importance of non-native interactions in immunity protein folding and demonstrate that the rate-limiting step involves partial reorganisation of these interactions as the TS ensemble is traversed. Moreover, we show that increasing the contribution to stability made by non-native interactions results in an increase in Phi-values of the TS ensemble without altering its structural properties or solvent-accessible surface area. The data suggest that the immunity proteins fold on multiple, but closely related, micropathways, resulting in a heterogeneous TS ensemble that responds subtly to mutation or changes in the solvent conditions. Thus, altering the relative strength of native and non-native interactions influences the search to the native state by restricting the pathways through the folding energy landscape.

Folding of tandem-linked domains

One of the factors, which influences protein folding in vivo, is a linkage of protein domains into multidomain tandems. However, relatively little is known about the impact of domain connectivity on protein folding mechanisms. In this article, we use coarse grained models of proteins to explore folding of tandem-linked domains (TLD). We found TLD folding to follow two scenarios. In the first, the tandem connectivity produces relatively minor impact on folding and the mechanisms of folding of tandem-linked and single domains remain similar. The second scenario involves qualitative changes in folding mechanism because of tandem linkage. As a result, protein domains, which fold via two-state mechanism as single isolated domains, may form new stable intermediates when inserted into tandems. The new intermediates are created by topological constraints imposed by the linkers between domains. In both cases tandem linkage slows down folding. We propose that the impact of tandem connectivity can be minimized, if the terminal secondary structure elements (SSEs) are flexible. In particular, two factors appear to facilitate TLD folding: (1) the interactions between terminal SSE are poorly ordered in the folding transition state, whereas nonterminal SSE are better structured, (2) the interactions between terminal SSE are weak in the native state. We apply these findings to wild-type proteins by examining experimental phi-value data and by performing all-atom molecular dynamics simulations. We show that immunoglobulin-like domains appear to utilize the factors, which minimize the impact of tandem connectivity on their folding. Several single domain proteins, which are likely to misfold in tandems, are also identified.

Temperature-dependent folding pathways of Pin1 WW domain: an all-atom molecular dynamics simulation of a Gō model

We study the folding thermodynamics and kinetics of the Pin1 WW domain, a three-stranded beta-sheet protein, by using all-atom (except nonpolar hydrogens) discontinuous molecular dynamics simulations at various temperatures with a Gō model. The protein exhibits a two-state folding kinetics near the folding transition temperature. A good agreement between our simulations and the experimental measurements by the Gruebele group has been found, and the simulation sheds new insights into the structure of transition state, which is hard to be straightforwardly captured in experiments. The simulation also reveals that the folding pathways at approximately the transition temperature and at low temperatures are much different, and an intermediate state at a low temperature is predicted. The transition state of this small beta-protein at its folding transition temperature has a well-established hairpin 1 made of beta1 and beta2 strands while its low-temperature kinetic intermediate has a formed hairpin 2 composed of beta2 and beta3 strands. Theoretical results are compared with other simulation results as well as available experimental data. This study confirms that specific side-chain packing in an all-atom Gō model can yield a reasonable prediction of specific folding kinetics for a given protein. Different folding behaviors at different temperatures are interpreted in terms of the interplay of entropy and enthalpy in folding process.

Site-specific collapse dynamics guide the formation of the cytochrome c' four-helix bundle

The evolution of tryptophan-to-heme (W/heme) distance distributions extracted from analysis of fluorescence energy transfer kinetics during the refolding of Rhodopseudomonas palustris cytochrome c' reveals dramatic differences between two variants [W32 (Q1A/F32W/W72F) and W72 (Q1A)]. Both W32/heme and W72/heme distance distributions measured at the earliest time point attainable with a continuous-flow mixer (150 mus) confirm that the polypeptide ensemble is not uniformly collapsed and that native structure is not formed. Time-resolved fluorescence spectra indicate that W32 is sequestered from the aqueous solution during the first 700 mus of folding, whereas W72 remains exposed to solvent. The first moment of the W32/heme distance distribution evolves to its native value faster than that of W72, suggesting that the approach of W32 to the heme precedes that of W72.

Characterisation of the Conformational Properties of Urea-unfolded Im7: Implications for the Early Stages of Protein Folding

The colicin immunity protein Im7 folds from its unfolded state in 6 M urea to its native four-helix structure through an on-pathway intermediate that lacks one of the helices of the native structure (helix III). In order to further characterize the folding mechanism of Im7, we have studied the conformational properties of the protein unfolded in 6 M urea in detail using heteronuclear NMR. Triple-resonance experiments with 13C/15N-labelled Im7 in 6 M urea provided almost complete resonance assignments for the backbone nuclei, and measurement of backbone 15N relaxation parameters allowed dynamic ordering of the unfolded polypeptide chain to be investigated. Reduced spectral density mapping and fitting backbone R2 relaxation rates to a polymer dynamics model identified four clusters of interacting residues, each predicted by the average area buried upon folding for each residue. Chemical shift analyses and measurement of NOEs detected with a long mixing-time 1H-1H-15N NOESY-HSQC spectrum confirmed the formation of four clusters. Each cluster of interacting side-chains in urea-unfolded Im7 occurs in a region of the protein that forms a helix in the protein, with the largest clusters being associated with the three long helices that are formed in the on-pathway folding intermediate, whilst the smallest cluster forms a helix only in the native state. NMR studies of a Phe15Ala Im7 variant and a protein in which residues 51-56 are replaced by three glycine residues (H3G3 Im7*), indicated that the clusters do not interact with each other, possibly because they are solvated by urea, as indicated by analysis of NOEs between the protein and the solvent. Based on these data, we suggest that dilution of the chaotrope to initiate refolding will result in collapse of the clusters, leading to the formation of persistent helical structure and the generation of the three-helix folding intermediate.

Structural interpretation of hydrogen exchange protection factors in proteins: characterization of the native state fluctuations of CI2

Protection factors obtained from equilibrium hydrogen exchange experiments are an important source of structural information on both native and nonnative states of proteins. We present a method for determining ensembles of protein structures by using hydrogen exchange data as restraints in molecular dynamics simulations in conjunction with an empirical force-field. The method is applied to determine the ensemble of structures representing the native state of chymotrypsin inhibitor 2 (CI2), including the rare, large fluctuations responsible for hydrogen exchange.

Structural comparison of the two alternative transition states for folding of TI I27

TI I27, a beta-sandwich domain from the human muscle protein titin, has been shown to fold via two alternative pathways, which correspond to a change in the folding mechanism. Under physiological conditions, TI I27 folds by a classical nucleation-condensation mechanism (diffuse transition state), whereas at extreme conditions of temperature and denaturant it switches to having a polarized transition state. We have used experimental Phi-values as restraints in ensemble-averaged molecular dynamics simulations to determine the ensembles of structures representing the two transition states. The comparison of these ensembles indicates that when native interactions are substantially weakened, a protein may still be able to fold if it can access an alternative transition state characterized by a much larger entropic contribution. Analysis of the probability distribution of Phi-values derived from ensemble averaged simulations, enables us to identify residues that form contacts in some members of the ensemble but not in others illustrating that many interactions present in transition states are not strictly required for the successful completion of the folding process.

Engineering Stabilising β-Sheet Interactions into a Conformationally Flexible Region of the Folding Transition State of Ubiquitin

Protein engineering studies suggest that the transition state for the folding of ubiquitin is highly polarised towards the N-terminal part of the sequence and involves a nucleus of residues within the beta-hairpin (residues 1-17) and main alpha-helix (residues 23-34). In contrast, the observation of small phi-values for residues in the C-terminal portion of the sequence (residues 35-76), coupled with a folding topology that results in a much higher contact order, suggests that fast folding of ubiquitin is dependent upon configurational flexibility in the C-terminal part of the polypeptide chain to ensure passage down a relatively smooth folding funnel to the native state. We show that the introduction of a small mini-hairpin motif as an extension of the native 43-50 hairpin stabilises local interactions in the C-terminal part of the sequence, resulting largely in a deceleration of the unfolding kinetics without perturbing the apparent two-state folding mechanism. However, a single-point Leu-->Phe substitution within the engineered hairpin sequence leads to the premature collapse of the denatured ensemble through the stabilisation of non-native interactions and the population of a compact intermediate. Non-linear effects in the kinetic data at low concentrations of denaturant suggest that the collapsed state, which is further stabilised in the presence of cosmotropic salts, may subsequently fold directly to the native state through a "triangular" reaction scheme involving internal rearrangement rather than unfolding and refolding.

Transition State Contact Orders Correlate with Protein Folding Rates

We have used molecular dynamics simulations restrained by experimental phi values derived from protein engineering experiments to determine the structures of the transition state ensembles of ten proteins that fold with two-state kinetics. For each of these proteins we then calculated the average contact order in the transition state ensemble and compared it with the corresponding experimental folding rate. The resulting correlation coefficient is similar to that computed for the contact orders of the native structures, supporting the use of native state contact orders for predicting folding rates. The native contacts in the transition state also correlate with those of the native state but are found to be about 30% lower. These results show that, despite the high levels of heterogeneity in the transition state ensemble, the large majority of contributing structures have native-like topologies and that the native state contact order captures this phenomenon.

Determination of the folding transition states of barnase by using PhiI-value-restrained simulations validated by double mutant PhiIJ-values

The protein barnase folds from the denatured state into its native conformation via a high-energy intermediate. Using PhiI-values determined experimentally from single-point mutations as restraints in all-atom molecular dynamics simulations, we have determined ensembles of structures corresponding to the transition states for the formation of the folding intermediate and its conversion into the native state. We have also introduced a stringent validation of the approach used to calculate such structures by predicting interaction PhiIJ-values determined experimentally from a series of double-mutant cycles. The ensembles that we have obtained illustrate the heterogeneity in the nucleation-condensation process by which barnase folds. Obligatory steps of this process include the sequential formation of two folding nuclei, which correspond to the two main hydrophobic cores of the protein. Nonobligatory steps include the elongation of the strand beta1 and the formation of the helix alpha2. The results confirm that the use of experimental observables such as PhiI-values as restraints in molecular dynamics simulations is a powerful general strategy to characterize the relatively heterogeneous structural ensembles that populate nonnative regions of the energy landscapes of proteins.

Balancing energy and entropy: a minimalist model for the characterization of protein folding landscapes

Coarse-grained models have been extremely valuable in promoting our understanding of protein folding. However, the quantitative accuracy of existing simplified models is strongly hindered either from the complete removal of frustration (as in the widely used Gō-like models) or from the compromise with the minimal frustration principle and/or realistic protein geometry (as in the simple on-lattice models). We present a coarse-grained model that "naturally" incorporates sequence details and energetic frustration into an overall minimally frustrated folding landscape. The model is coupled with an optimization procedure to design the parameters of the protein Hamiltonian to fold into a desired native structure. The application to the study of src-Src homology 3 domain shows that this coarse-grained model contains the main physical-chemical ingredients that are responsible for shaping the folding landscape of this protein. The results illustrate the importance of nonnative interactions and energetic heterogeneity for a quantitative characterization of folding mechanisms.

Towards complete descriptions of the free-energy landscapes of proteins

In recent years increasingly detailed information about the structures and dynamics of protein molecules has been obtained by innovative applications of experimental techniques, in particular nuclear magnetic resonance spectroscopy and protein engineering, and theoretical methods, notably molecular dynamics simulations. In this article we discuss how such approaches can be combined by incorporating a wide range of different types of experimental data as restraints in computer simulations to provide unprecedented detail about the ensembles of structures that describe proteins in a wide variety of states from the native structure to highly unfolded species. Knowledge of these ensembles is beginning to enable the complete free-energy landscapes of individual proteins to be defined at atomic resolution. This strategy has provided new insights into the mechanism by which proteins are able to fold into their native states, or by which they fail to do so and give rise to harmful aggregates that are associated with a wide range of debilitating human diseases.

Different secondary structure elements as scaffolds for protein folding transition states of two homologous four-helix bundles

Comparison of the folding processes for homologue proteins can provide valuable information about details in the interactions leading to the formation of the folding transition state. Here the folding kinetics of 18 variants of yACBP and 3 variants of bACBP have been studied by Phi-value analysis. In combination with Phi-values from previous work, detailed insight into the transition states for folding of both yACBP and bACBP has been obtained. Of the 16 sequence positions that have been studied in both yACBP and bACBP, 5 (V12, I/L27, Y73, V77, and L80) have high Phi-values and appear to be important for the transition state formation in both homologues. Y31, A34, and A69 have high Phi-values only in yACBP, while F5, A9, and I74 have high Phi-values only in bACBP. Thus, additional interactions between helices A2 and A4 appear to be important for the transition state of yACBP, whereas additional interactions between helices A1 and A4 appear to be important for the transition state of bACBP. To examine whether these differences could be assigned to different packing of the residues in the native state, a solution structure of yACBP was determined by NMR. Small changes in the packing of the hydrophobic side-chains, which strengthen the interactions between helices A2 and A4, are observed in yACBP relative to bACBP. It is suggested that different structure elements serve as scaffolds for the folding of the 2 ACBP homologues.

Protein folding and the organization of the protein topology universe

The mechanism by which proteins fold to their native states has been the focus of intense research in recent years. The rate-limiting event in the folding reaction is the formation of a conformation in a set known as the transition-state ensemble. The structural features present within such ensembles have now been analysed for a series of proteins using data from a combination of biochemical and biophysical experiments together with computer-simulation methods. These studies show that the topology of the transition state is determined by a set of interactions involving a small number of key residues and, in addition, that the topology of the transition state is closer to that of the native state than to that of any other fold in the protein universe. Here, we review the evidence for these conclusions and suggest a molecular mechanism that rationalizes these findings by presenting a view of protein folds that is based on the topological features of the polypeptide backbone, rather than the conventional view that depends on the arrangement of different types of secondary-structure elements. By linking the folding process to the organization of the protein structure universe, we propose an explanation for the overwhelming importance of topology in the transition states for protein folding.

Comparison of the transition states for folding of two Ig-like proteins from different superfamilies

In the "fold approach" proteins with a similar fold but different sequences are compared in order to investigate the relationship between native state structure and folding behaviour. Here we compare the properties of the transition states for folding of TI I27, the 27th immunoglobulin domain from human cardiac titin, and that of TNfn3, the third fibronectin type III domain from human tenascin. Experimental phi-values were used as restraints in molecular dynamics simulations to determine the structures that make up the transition state ensembles (TSEs) for folding of the two proteins. The restrained simulations that we present allow a detailed structural comparison of the two TSEs to be made. Further calculations show explicitly that for both proteins the formation of the interactions involving the residues in the folding nucleus is sufficient for the establishment of the topology of the Ig-like fold. We found that, although the folding nuclei of the two proteins are similar, the packing of the folding nucleus of TI I27 is much tighter than that of TNfn3, reflecting the higher experimental phi-values and beta(T) (Tanford Beta) of TI I27. These results suggest that the folding nucleus can be significantly deformed to accommodate extensive sequence variation while conserving the same folding mechanism.

The folding of spectrin domains II: phi-value analysis of R16

Studies on the folding of helical proteins have shown a wide range of different mechanisms and highlighted the importance of helical propensity as a factor in determining folding mechanism. Here, we contribute to this interesting field with the protein engineering phi-value analysis of the 16th domain of chicken brain alpha-spectrin, R16. The fortuitous curvature seen in the unfolding arm of the chevron plot allows us to investigate both early and late events in folding. R16 is the first two-state helical protein for which this has been possible.

Optimal combination of theory and experiment for the characterization of the protein folding landscape of S6: how far can a minimalist model go?

The detailed characterization of the overall free energy landscape associated with the folding process of a protein is the ultimate goal in protein folding studies. Modern experimental techniques provide accurate thermodynamic and kinetic measurements on restricted regions of a protein landscape. Although simplified protein models can access larger regions of the landscape, they are oftentimes built on assumptions and approximations that affect the accuracy of the results. We present a new methodology that allows to combine the complementary strengths of theory and experiment for a more complete characterization of a protein folding landscape. We prove that this new procedure allows a simplified protein model to reproduce remarkably well (correlation coefficient > 0.9) all experimental data available on free energies differences upon single mutations for S6 ribosomal protein and two circular permutants. Our results confirm and quantify the hypothesis, recently formulated on the basis of experimental data, that the folding landscape of protein S6 is strongly affected by an atypical distribution of contact energies.

Trapping the on-pathway folding intermediate of Im7 at equilibrium

The four-helical protein Im7 folds via a rapidly formed on-pathway intermediate (k(UI)=3000 s(-1) at pH 7.0, 10 degrees C) that contains three (helices I, II and IV) of the four native alpha-helices. The relatively slow (k(IN)=300 s(-1)) conversion of this intermediate into the native structure is driven by the folding and docking of the six residue helix III onto the developing hydrophobic core. Here, we describe the structural properties of four Im7* variants designed to trap the protein in the intermediate state by disrupting the stabilising interactions formed between helix III and the rest of the protein structure. In two of these variants (I54A and L53AI54A), hydrophobic residues within helix III have been mutated to alanine, whilst in the other two mutants the sequence encompassing the native helix III was replaced by a glycine linker, three (H3G3) or six (H3G6) residues in length. All four variants were shown to be monomeric, as judged by analytical ultracentrifugation, and highly helical as measured by far-UV CD. In addition, all the variants denature co-operatively and have a stability (DeltaG(UF)) and buried hydrophobic surface area (M(UF)) similar to those of the on-pathway kinetic intermediate. Structural characterisation of these variants using 1-anilino-8-napthalene sulphonic acid (ANS) binding, near-UV CD and 1D (1)H NMR demonstrate further that the trapped intermediate ensemble is highly structured with little exposed hydrophobic surface area. Interestingly, however, the structural properties of the variants I54A and L53AI54A differ in detail from those of H3G3 and H3G6. In particular, the single tryptophan residue, located near the end of helix IV, and distant from helix III, is in a distinct environment in the two sets of mutants as judged by fluorescence, near-UV CD and the sensitivity of tryptophan fluorescence to iodide quenching. Overall, the results confirm previous kinetic analysis that demonstrated the hierarchical folding of Im7 via an on-pathway intermediate, and show that this species is a highly helical ensemble with a well-formed hydrophobic core. By contrast with the native state, however, the intermediate ensemble is flexible enough to change in response to mutation, its structural properties being tailored by residues in the sequence encompassing the native helix III.

Switching Two-state to Three-state Kinetics in the Helical Protein Im9 via the Optimisation of Stabilising Non-native Interactions by Design

The four-helix protein Im7 folds through an on-pathway intermediate at pH 7.0 and 10 degrees C. By contrast, under these conditions there is no evidence for a populated intermediate in the folding of its more stable homologue, Im9, even in the presence of 0.4 M sodium sulphate. Previous studies using phi-value analysis have shown that the Im7 intermediate is misfolded, in that three of its four native helices are formed, but are docked in a non-native manner. Using knowledge of the structure of the intermediate of Im7, we have used rational design to stabilise an intermediate formed during the folding of Im9 by the introduction of specific stabilising interactions at positions known to stabilise the Im7 folding intermediate through non-native interactions. We show that the redesigned Im9 sequence folds with three-state kinetics at pH 7.0 and have used phi-value analysis to demonstrate that this species resembles the misfolded intermediate populated during Im7 folding. The redesigned Im9 sequence folds 20-fold faster than the wild-type protein under conditions in which folding is two-state. The data show that intermediate formation is an important feature of folding, even for small proteins such as Im9 for which these partially folded states do not become significantly populated. In addition, they show that the introduction of stabilising interactions can lead to rapid refolding, even when the contacts introduced are non-native.