In silico protein fragmentation reveals the importance of critical nuclei on domain reassembly

Protein complementation assays (PCAs) based on split protein fragments have become powerful tools that facilitate the study and engineering of intracellular protein-protein interactions. These assays are based on the observation that a given protein can be split into two inactive fragments and these fragments can reassemble into the original properly folded and functional structure. However, one experimentally observed limitation of PCA systems is that the folding of a protein from its fragments is dramatically slower relative to that of the unsplit parent protein. This is due in part to a poor understanding of how PCA design parameters such as split site position in the primary sequence and size of the resulting fragments contribute to the efficiency of protein reassembly. We used a minimalist on-lattice model to analyze how the dynamics of the reassembly process for two model proteins was affected by the location of the split site. Our results demonstrate that the balanced distribution of the “folding nucleus,” a subset of residues that are critical to the formation of the transition state leading to productive folding, between protein fragments is key to their reassembly.     

Thermodynamics and folding kinetics analysis of the SH3 domain form discrete molecular dynamics

We perform a detailed analysis of the thermodynamics and folding kinetics of the SH3 domain fold with discrete molecular dynamic simulations. We propose a protein model that reproduces some of the experimentally observed thermodynamic and folding kinetic properties of proteins. Specifically, we use our model to study the transition state ensemble of the SH3 fold family of proteins, a set of unstable conformations that fold to the protein native state with probability 1/2. We analyze the participation of each secondary structure element formed at the transition state ensemble. We also identify the folding nucleus of the SH3 fold and test extensively its importance for folding kinetics. We predict that a set of amino acid contacts between the RT-loop and the distal hairpin are the critical folding nucleus of the SH3 fold and propose a hypothesis that explains this result.     

Influence of denatured and intermediate states of folding on protein aggregation

We simulate the aggregation thermodynamics and kinetics of proteins L and G, each of which self-assembles to the same alpha/beta [corrected] topology through distinct folding mechanisms. We find that the aggregation kinetics of both proteins at an experimentally relevant concentration exhibit both fast and slow aggregation pathways, although a greater proportion of protein G aggregation events are slow relative to those of found for protein L. These kinetic differences are correlated with the amount and distribution of intrachain contacts formed in the denatured state ensemble (DSE), or an intermediate state ensemble (ISE) if it exists, as well as the folding timescales of the two proteins. Protein G aggregates more slowly than protein L due to its rapidly formed folding intermediate, which exhibits native intrachain contacts spread across the protein, suggesting that certain early folding intermediates may be selected for by evolution due to their protective role against unwanted aggregation. Protein L shows only localized native structure in the DSE with timescales of folding that are commensurate with the aggregation timescale, leaving it vulnerable to domain swapping or nonnative interactions with other chains that increase the aggregation rate. Folding experiments that characterize the structural signatures of the DSE, ISE, or the transition state ensemble (TSE) under nonaggregating conditions should be able to predict regions where interchain contacts will be made in the aggregate, and to predict slower aggregation rates for proteins with contacts that are dispersed across the fold. Since proteins L and G can both form amyloid fibrils, this work also provides mechanistic and structural insight into the formation of prefibrillar species.     

Stability and folding kinetics of a ubiquitin mutant with a strong propensity for nonnative beta-hairpin conformation in the unfolded state

A F45W mutant of yeast ubiquitin has been used as a model system to examine the effects of nonnative local interactions on protein folding and stability. Mutating the native TLTGK G-bulged type I turn in the N-terminal beta-hairpin to NPDG stabilizes a nonnative beta-strand alignment in the isolated peptide fragment. However, NMR structural analysis of the native and mutant proteins shows that the NPDG mutant is forced to adopt the native beta-strand alignment and an unfavorable type I NPDG turn. The mutant is significantly less stable (approximately 9 kJ mol(-1)) and folds 30 times slower than the native sequence, demonstrating that local interactions can modulate protein stability and that attainment of a nativelike beta-hairpin conformation in the transition state ensemble is frustrated by the turn mutations. Surprising, alcoholic cosolvents [5-10% (v/v) TFE] are shown to accelerate the folding rate of the NPDG mutant. We conclude, backed-up by NMR data on the peptide fragments, that even though nonnative states in the denatured ensemble are highly populated and their stability further enhanced in the presence of cosolvents, the simultaneous increase in the proportion of nativelike secondary structure (hairpin or helix), in rapid equilibrium with nonnative states, is sufficient to accelerate the folding process. It is evident that modulating local interactions and increasing nonnative secondary structure propensities can change protein stability and folding kinetics. However, nonlocal contacts formed in the global cooperative folding event appear to determine structural specificity.     

Reconstruction of the src-SH3 protein domain transition state ensemble using multiscale molecular dynamics simulations

We use an integrated computational approach to reconstruct accurately the transition state ensemble (TSE) for folding of the src-SH3 protein domain. We first identify putative TSE conformations from free energy surfaces generated by importance sampling molecular dynamics for a fully atomic, solvated model of the src-SH3 protein domain. These putative TSE conformations are then subjected to a folding analysis using a coarse-grained representation of the protein and rapid discrete molecular dynamics simulations. Those conformations that fold to the native conformation with a probability (P(fold)) of approximately 0.5, constitute the true transition state. Approximately 20% of the putative TSE structures were found to have a P(fold) near 0.5, indicating that, although correct TSE conformations are populated at the free energy barrier, there is a critical need to refine this ensemble. Our simulations indicate that the true TSE conformations are compact, with a well-defined central beta sheet, in good agreement with previous experimental and theoretical studies. A structured central beta sheet was found to be present in a number of pre-TSE conformations, however, indicating that this element, although required in the transition state, does not define it uniquely. An additional tight cluster of contacts between highly conserved residues belonging to the diverging turn and second beta-sheet of the protein emerged as being critical elements of the folding nucleus. A number of commonly used order parameters to identify the transition state for folding were investigated, with the number of native Cbeta contacts displaying the most satisfactory correlation with P(fold) values.     

Identifying folding nucleus based on residue contact networks of proteins

In the native structure of a protein, all the residues are tightly parked together in a specific order following its folding and every residue contacts with some spatially neighbor residues. A residue contact network can be constructed by defining the residues as nodes and the native contacts as edges. During the folding of small single-domain proteins, there is a set of contacts (or bonds), defined as the folding nucleus (FN), which is formed around the transition state, i.e., a rate-limiting barrier located at about the middle between the unfolded states and the native state on the free energy landscape. Such a FN plays an essential role in the folding dynamics and the residues, which form the related contacts called as folding nucleus residues (FNRs). In this work, the FNRs in proteins are identified by using quantities which characterize the topology of residue contact networks of proteins. By comparing the specificities of residues with the network quantities K(R), L(R), and D(R), up to 90% FNRs of six typical proteins found experimentally are identified. It is found that the FNRs behave the full-closeness centrals rather than degree or closeness centers in the residue contact network, implying that they are important to the folding cooperativity of proteins. Our study shows that the FNRs can be identified solely from the native structures of proteins based on the analysis of residue contact network without any knowledge of the transition state ensemble.     

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.     

Commitment and nucleation in the protein G transition state

An accurate characterization of the transition state ensemble (TSE) is central to furthering our understanding of the protein folding reaction. We have extensively tested a recently reported method for studying a protein's TSE, utilizing phi-value data from protein engineering experiments and computational studies as restraints in all-atom Monte Carlo (MC) simulations. The validity of interpreting experimental phi-values as the fraction of native contacts made by a residue in the TSE was explored, revealing that this definition is unable to uniquely specify a TSE. The identification of protein G's second hairpin, in both pre and post-transition conformations demonstrates that high experimental phi-values do not guarantee a residue's importance in the TSE. An analysis of simulations based on structures restrained by experimental phi-values is necessary to yield this result, which is not obvious from a simplistic interpretation of individual phi-values. The TSE that we obtain corresponds to a single, specific nucleation event, characterized by six residues common to all three observed, convergent folding pathways. The same specific nucleus was independently identified from computational and experimental data, and "Conservation of Conservation" analysis in the protein G fold. When associated strictly with complete nucleus formation and concomitant chain collapse, folding is a well-defined two state event. Once the nucleus has formed, the folding reaction enters a slow relaxation process associated with side-chain packing and small, local backbone rearrangements. A detailed analysis of phi-values and their relationship to the transition state ensemble allows us to construct a unified theoretical model of protein G folding.     

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.     

The Origin of Nonmonotonic Complex Behavior and the Effects of Nonnative Interactions on the Diffusive Properties of Protein Folding

We present a method for calculating the configurational-dependent diffusion coefficient of a globular protein as a function of the global folding process. Using a coarse-grained structure-based model, we determined the diffusion coefficient, in reaction coordinate space, as a function of the fraction of native contacts formed Q for the cold shock protein (TmCSP). We find nonmonotonic behavior for the diffusion coefficient, with high values for the folded and unfolded ensembles and a lower range of values in the transition state ensemble. We also characterized the folding landscape associated with an energetically frustrated variant of the model. We find that a low-level of frustration can actually stabilize the native ensemble and increase the associated diffusion coefficient. These findings can be understood from a mechanistic standpoint, in that the transition state ensemble has a more homogeneous structural content when frustration is present. Additionally, these findings are consistent with earlier calculations based on lattice models of protein folding and more recent single-molecule fluorescence measurements.     

Folding pathway dependence on energetic frustration and interaction heterogeneity for a three-dimensional hydrophobic protein model

Monte Carlo simulations of a hydrophobic protein model of 40 monomers in the cubic lattice are used to explore the effect of energetic frustration and interaction heterogeneity on its folding pathway. The folding pathway is described by the dependence of relevant conformational averages on an appropriate reaction coordinate, pfold, defined as the probability for a given conformation to reach the native structure before unfolding. We compare the energetically frustrated and heterogeneous hydrophobic potential, according to which individual monomers have a higher or lower tendency to form contacts unspecifically depending on their hydrophobicities, to an unfrustrated homogeneous Go-type potential with uniformly attractive native interactions and neutral non-native interactions (called Go1 in this study), and to an unfrustrated heterogeneous potential with neutral non-native interactions and native interactions having the same energy as the hydrophobic potential (called Go2 in this study). Folding kinetics are slowed down dramatically when energetic frustration increases, as expected and previously observed in a two-dimensional model. Contrary to our previous results in two dimensions, however, it appears that the folding pathway and transition state ensemble can be significantly dependent on the energy function used to stabilize the native structure. The sequence of events along the reaction coordinate, or the order along this coordinate in which different regions of the native conformation become structured, turns out to be similar for the hydrophobic and Go2 potentials, but with analogous events tending to occur at lower pfold values in the first case. In particular, the transition state obtained from the ensemble around pfold = 0.5 is more structured for the hydrophobic potential. For Go1, not only the transition state ensemble but the order of events itself is modified, suggesting that interaction heterogeneity, in addition to energetic frustration, can have significant effects on the folding mechanism, most likely by modifying the probability of different contacts in the unfolded state, the starting point for the folding reaction. Although based on a simple model, these results provide interesting insight into how sequence-dependent switching between folding pathways might occur in real proteins.     

Thermodynamics and Kinetics of Single-Chain Monellin Folding with Structural Insights into Specific Collapse in the Denatured State Ensemble

Proteins, which behave as random coils in high denaturant concentrations undergo collapse transition similar to polymers on denaturant dilution. We study collapse in the denatured ensemble of single-chain monellin (MNEI) using a coarse-grained protein model and molecular dynamics simulations. The model is validated by quantitatively comparing the computed guanidinium chloride and pH-dependent thermodynamic properties of MNEI folding with the experiments. The computed properties such as the fraction of the protein in the folded state and radius of gyration (R_g) as function of [GuHCl] are in good agreement with the experiments. The folded state of MNEI is destabilized with an increase in pH due to the deprotonation of the residues Glu24 and Cys42. On decreasing [GuHCl], the protein in the unfolded ensemble showed specific compaction. The R_g of the protein decreased steadily with [GuHCl] dilution due to increase in the number of native contacts in all the secondary structural elements present in the protein. MNEI folding kinetics is complex with multiple folding pathways and transiently stable intermediates are populated in these pathways. In strong stabilizing conditions, the protein in the unfolded ensemble showed transition to a more compact unfolded state where R_g decreased by ≈ 17% due to the formation of specific native contacts in the protein. The intermediate populated in the dominant MNEI folding pathway satisfies the structural features of the dry molten globule inferred from experiments.     

Prediction of folding transition-state position (betaT) of small, two-state proteins from local secondary structure content

Folding kinetics of proteins is governed by the free energy and position of transition states. But attempts to predict the position of folding transition state on reaction pathway from protein structure have been met with only limited success, unlike the folding-rate prediction. Here, we find that the folding transition-state position is related to the secondary structure content of native two-state proteins. We present a simple method for predicting the transition-state position from their alpha-helix, turn and polyproline secondary structures. The method achieves 81% correlation with experiment over 24 small, two-state proteins, suggesting that the local secondary structure content, especially for content of alpha-helix, is a determinant of the solvent accessibility of the transition state ensemble and size of folding nucleus.     

Non-native interactions,effective contact order,and protein folding: a mutational investigation with the energetically frustrated hydrophobic model

By Monte Carlo simulations, we explored the effect of single mutations on the thermodynamics and kinetics of the folding of a two-dimensional, energetically frustrated, hydrophobic protein model. Phi-Value analysis, corroborated by simulations beginning from given sets of judiciously chosen initial contacts, suggests that the transition state of the model consists of a limited region of the native structure, that is, a folding nucleus. It seems that the most important contacts in the transition state (large and positive Phi) are not the ones with the highest contact order, because in this case the entropic cost of their formation would be too high, but exactly the ones that decrease the entropic cost of difficult contacts, reducing their effective contact order. Mutations of internal monomers involved in high-order contacts were actually the ones resulting in the fastest kinetics (and Phi < 0), indicating they tend to make low order, non-native contacts of low entropic cost that stabilize the unfolded state with respect to the transition state. Folding acceleration by other non-native interactions was also observed and a simple general mechanism is proposed according to which non-native contacts can act indirectly over the folding nucleus, "chelating" out potentially harmful contacts. The polymer graph of our model, which facilitates the visualization of effective contact orders, successfully suggests the relative kinetic importance of different contacts and is reasonably consistent with analogous graphs for the well characterized family of SH3 domains.     

Nucleation and the transition state of the SH3 domain

We present a verified computational model of the SH3 domain transition state (TS) ensemble. This model was built for three separate SH3 domains using experimental phi-values as structural constraints in all-atom protein folding simulations. While averaging over all conformations incorrectly considers non-TS conformations as transition states, quantifying structures as pre-TS, TS, and post-TS by measurement of their transmission coefficient ("probability to fold", or p(fold)) allows for rigorous conclusions regarding the structure of the folding nucleus and a full mechanistic analysis of the folding process. Through analysis of the TS, we observe a highly polarized nucleus in which many residues are solvent-exposed. Mechanistic analysis suggests the hydrophobic core forms largely after an early nucleation step. SH3 presents an ideal system for studying the nucleation-condensation mechanism and highlights the synergistic relationship between experiment and simulation in the study of protein folding.     

Experiment and theory highlight role of native state topology in SH3 folding.

We use a combination of experiments, computer simulations and simple model calculations to characterize, first, the folding transition state ensemble of the src SH3 domain, and second, the features of the protein that determine its folding mechanism. Kinetic analysis of mutations at 52 of the 57 residues in the src SH3 domain revealed that the transition state ensemble is even more polarized than suspected earlier: no single alanine substitution in the N-terminal 15 residues or the C-terminal 9 residues has more than a two-fold effect on the folding rate, while such substitutions at 15 sites in the central three-stranded beta-sheet cause significant decreases in the folding rate. Molecular dynamics (MD) unfolding simulations and ab initio folding simulations on the src SH3 domain exhibit a hierarchy of folding similar to that observed in the experiments. The similarity in folding mechanism of different SH3 domains and the similar hierarchy of structure formation observed in the experiments and the simulations can be largely accounted for by a simple native state topology-based model of protein folding energy landscapes.     

Amino Acid Insertion Reveals a Necessary Three-Helical Intermediate in the Folding Pathway of the Colicin E7 Immunity Protein Im7

The small (87-residue) α-helical protein Im7 (an inhibitor protein for colicin E7 that provides immunity to cells producing colicin E7) folds via a three-state mechanism involving an on-pathway intermediate. This kinetic intermediate contains three of four native helices that are oriented in a non-native manner so as to minimise exposed hydrophobic surface area at this point in folding. The short (6-residue) helix III has been shown to be unstructured in the intermediate ensemble and does not dock onto the developing hydrophobic core until after the rate-limiting transition state has been traversed. After helix III has docked, it adopts an α-helical secondary structure, and the side chains of residues within this region provide contacts that are crucial to native-state stability. In order to probe further the role of helix III in the folding mechanism of Im7, we created a variant that contains an eight-amino-acid polyalanine-like helix stabilised by a Glu-Arg salt bridge and an Asn-Pro-Gly capping motif, juxtaposed C-terminal to the natural 6-residue helix III. The effect of this insertion on the structure of the native protein and its folding mechanism were studied using NMR and ?-value analysis, respectively. The results reveal a robust native structure that is not perturbed by the presence of the extended helix III. Mutational analysis performed to probe the folding mechanism of the redesigned protein revealed a conserved mechanism involving the canonical three-helical intermediate. The results suggest that folding via a three-helical species stabilised by both native and non-native interactions is an essential feature of Im7 folding, independent of the helical propensity of helix III.     

Stiffness of the distal loop restricts the structural heterogeneity of the transition state ensemble in SH3 domains

Protein engineering experiments and Phi(F)-value analysis of SH3 domains reveal that their transition state ensemble (TSE) is conformationally restricted, i.e. the fluctuations in the transition state (TS) structures are small. In the TS of src SH3 and alpha-spectrin SH3 the distal loop and the associated hairpin are fully structured, while the rest of the protein is relatively disordered. If native structure predominantly determines the folding mechanism, the findings for SH3 folds raise the question: What are the features of the native topology that determine the nature of the TSE? We propose that the presence of stiff loops in the native state that connect local structural elements (such as the distal hairpin in SH3 domains) conformationally restricts TSE. We validate this hypothesis using the simulations of a "control" system (16 residue beta-hairpin forming C-terminal fragment of the GBl protein) and its variants. In these fragments the role of bending rigidity in determining the nature of the TSE can be directly examined without complications arising from interactions with the rest of the protein. The TSE structures in the beta-hairpins are determined computationally using cluster analysis and limited Phi(F)-value analysis. Both techniques prove that the conformational heterogeneity decreases as the bending rigidity of the loop increases. To extend this finding to SH3 domains a measure of bending rigidity based on loop curvature, which utilizes native structures in the Protein Data Bank (PDB), is introduced. Using this measure we show that, with few exceptions, the ordering of stiffness of the distal, n-src, and RT loops in the 29 PDB structures of SH3 domains is conserved. Combining the simulation results for beta-hairpins and the analysis of PDB structures for SH3 domains, we propose that the stiff distal loop restricts the conformational fluctuations in the TSE. We also predict that constraining the distal loop to be preformed in the denatured ensemble should not alter the nature of TSE. On the other hand, if the amino and carboxy terminals are cross-linked to form a circular polypeptide chain, the pathways and TSs are altered. These contrasting scenarios are illustrated using simulations of cross-linked WT beta-hairpin fragments. Computations of bending rigidities for immunoglobulin-like domain proteins reveal no clear separation in the stiffness of their loops. In the beta-sandwich proteins, which have large fractions of non-local native contacts, the nature of the TSE cannot be apparently determined using purely local structural characteristics. Nevertheless, the measure of loop stiffness still provides qualitative predictions of the ordered regions in the TSE of Ig27 and TenFn3.     

Pathway diversity and concertedness in protein folding: an ab-initio approach

Making use of an ab-initio folding simulator, we generate in vitro pathways leading to the native fold in moderate size single- domain proteins. The assessment of pathway diversity is not biased by any a priori information on the native fold. We focus on two study cases, hyperthermophile variant of protein G domain (1gb4) and ubiquitin (1ubi), with the same topology but different context dependence in their native folds. We demonstrate that a quenching of structural fluctuations is achieved once the proteins find a stationary plateau maximizing the number of highly protected hydrogen bonds. This enables us to identify the folding nucleus and show that folding does not become expeditious until a concerted event takes place generating a topology able to prevent water attack on a maximal number of hydrogen bonds. This result is consistent with the standard nucleation mechanism postulated for two-state folders. Pathway diversity is correlated with the extent of conflict between local structural propensity and large-scale context, rather than with contact order: In highly context-dependent proteins, the success of folding cannot rely on a single fortuitous event in which local propensity is overruled by large-scale effects. We predict mutational Pi values on individual pathways, compute ensemble averages and predict extent of surface burial and percentage of hydrogen bonding on each component of the transition state ensemble, thus deconvoluting individual folding-route contributions to the averaged two-state kinetic picture. Our predicted kinetic isotopic effects find experimental support and lead to further probes. Finally, the molecular redesign potentiality of the method, aimed at increasing folding expediency, is explored.