Specific nucleus as the transition state for protein folding: evidence from the lattice model

Small Proteins Fold Through Transition States With Native-like Topologies

The folding pathway of common-type acyl phosphatase (ctAcP) is characterized using psi-analysis, which identifies specific chain-chain contacts using bi-histidine (biHis) metal-ion binding sites. In the transition state ensemble (TSE), the majority of the protein is structured with a near-native topology, only lacking one beta-strand and an alpha-helix. psi-Values are zero or unity for all sites except one at the amino terminus of helix H2. This fractional psi-value remains unchanged when three metal ions of differing coordination geometries are used, indicating this end of the helix experiences microscopic heterogeneity through fraying in the TSE. Ubiquitin, the other globular protein characterized using psi-analysis, also exhibits a single consensus TSE structure. Hence, the TSE of both proteins have converged to a single configuration, albeit one that contains some fraying at the periphery. Models of the TSE of both proteins are created using all-atom Langevin dynamics simulations using distance constraints derived from the experimental psi-values. For both proteins, the relative contact order of the TS models is approximately 80% of the native value. This shared value viewed in the context of the known correlation between contact order and folding rates, suggests that other proteins will have a similarly high fraction of the native contact order. This constraint greatly limits the range of possible configurations at the rate-limiting step.

Testing simplified proteins models of the hPin1 WW domain

The WW domain of the human Pin1 protein for its simple topology and large amount of experimental data is an ideal candidate to assess theoretical approaches to protein folding. The purpose of this work is to compare the reliability of the chemically based Sorenson/Head-Gordon (SHG) model and a standard native centric model in reproducing, through molecular dynamics simulations, some of the well known features of the folding transition of this small domain. Our results show that the Gō model correctly reproduces the cooperative, two-state, folding mechanism of the WW-domain, while the SHG model predicts a transition occurring in two stages: a collapse, followed by a structural rearrangement. The lack of a cooperative folding in the SHG simulations appears to be related to the nonfunnel shape of the energy landscape featuring a partitioning of the native valley in subbasins corresponding to different chain chiralities. However, the SHG approach remains more reliable in estimating the phi-values with respect to Gō-like description. This may suggest that the WW-domain folding process is stirred by energetic and topological factors as well, and it highlights the better suitability of chemically based models in simulating mutations.

Hydrophobic collapse in (in silico) protein folding

A model of hydrophobic collapse, which is treated as the driving force for protein folding, is presented. This model is the superposition of three models commonly used in protein structure prediction: (1) 'oil-drop' model introduced by Kauzmann, (2) a lattice model introduced to decrease the number of degrees of freedom for structural changes and (3) a model of the formation of hydrophobic core as a key feature in driving the folding of proteins. These three models together helped to develop the idea of a fuzzy-oil-drop as a model for an external force field of hydrophobic character mimicking the hydrophobicity-differentiated environment for hydrophobic collapse. All amino acids in the polypeptide interact pair-wise during the folding process (energy minimization procedure) and interact with the external hydrophobic force field defined by a three-dimensional Gaussian function. The value of the Gaussian function usually interpreted as a probability distribution is treated as a normalized hydrophobicity distribution, with its maximum in the center of the ellipsoid and decreasing proportionally with the distance versus the center. The fuzzy-oil-drop is elastic and changes its shape and size during the simulated folding procedure.

Fuzzy-oil-drop hydrophobic force field--a model to represent late-stage folding (in silico) of lysozyme

A model of hydrophobic collapse (in silico), which is generally considered to be the driving force for protein folding, is presented in this work. The model introduces the external field in the form of a fuzzy-oil-drop assumed to represent the environment. The drop is expressed in the form of a three-dimensional Gauss function. The usual probability value is assumed to represent the hydrophobicity distribution in the three-dimensional space of the virtual environment. The differences between this idealized hydrophobicity distribution and the one represented by the folded polypeptide chain is the parameter to be minimized in the structure optimization procedure. The size of fuzzy-oil-drop is critical for the folding process. A strong correlation between protein length and the dimension of the native and early-stage molecular form was found on the basis of single-domain proteins analysis. A previously presented early-stage folding (in silico) model was used to create the starting structure for the procedure of late-stage folding of lysozyme. The results of simulation were found to be promising, although additional improvements for the formation of beta-structure and disulfide bonds as well as the participation of natural ligand in folding process seem to be necessary.

Hydrophobic collapse in late-stage folding (in silico) of bovine pancreatic trypsin inhibitor

Hydrophobic collapse is commonly considered as a process of significance for protein folding. Many models have been applied for description of this event. A model introducing an external force field mimicking the hydrophobic environment and simultaneously the driving force for the folding process is presented in this paper. Bovine pancreatic trypsin inhibitor (BPTI) was taken as a test protein. An early-stage folding (in silico) model presented elsewhere was used to create the starting structure for hydrophobic collapse. The resulting structure was energy-refined using molecular dynamics simulation in an explicit solvent. The similarity versus the crystal structure of BPTI is estimated using visual analysis, residue-residue contacts, phi, psi angle distributions, RMSD, accessible solvent area, radii of gyration and hydrodynamic radii. A program allowing creation of early-stage folding structural forms to be created for any protein is available from http://bioinformatics.cm-uj.krakow.pl/earlystage. The program for late-stage folding simulation is available on request.

Identification of the minimal protein-folding nucleus through loop-entropy perturbations

To explore the plasticity and structural constraints of the protein-folding nucleus we have constructed through circular permutation four topological variants of the ribosomal protein S6. In effect, these topological variants represent entropy mutants with maintained spatial contacts. The proteins were characterized at two complementary levels of detail: by phi-value analysis estimating the extent of contact formation in the transition-state ensemble and by Hammond analysis measuring the site-specific growth of the folding nucleus. The results show that, although the loop-entropy alterations markedly influence the appearance and structural location of the folding nucleus, it retains a common motif of one helix docking against two strands. This nucleation motif is built around a shared subset of side chains in the center of the hydrophobic core but extends in different directions of the S6 structure following the permutant-specific differences in local loop entropies. The adjustment of the critical folding nucleus to alterations in loop entropies is reflected by a direct correlation between the phi-value change and the accompanying change in local sequence separation.

Effective potentials for folding proteins

A coarse-grained off-lattice model that is not biased in any way to the native state is proposed to fold proteins. To predict the native structure in a reasonable time, the model has included the essential effects of water in an effective potential. Two new ingredients, the dipole-dipole interaction and the local hydrophobic interaction, are introduced and are shown to be as crucial as the hydrogen bonding. The model allows successful folding of the wild-type sequence of protein G and may have provided important hints to the study of protein folding.

Physical origins of protein superfamilies

In this work, we discovered a fundamental connection between selection for protein stability and emergence of preferred structures of proteins. Using a standard exact three-dimensional lattice model we evolve sequences starting from random ones and determine the exact native structure after each mutation. Acceptance of mutations is biased to select for stable proteins. We found that certain structures, "wonderfolds", are independently discovered numerous times as native states of stable proteins in many unrelated runs of selection. The strong dependence of lattice fold usage on the structural determinant of designability quantitatively reproduces uneven fold usage in natural proteins. Diversity of sequences that fold into wonderfold structures gives rise to superfamilies, i.e. sets of dissimilar sequences that fold into the same or very similar structures. The present work establishes a model of pre-biotic structure selection, which identifies dominant structural patterns emerging upon optimization of proteins for survival in a hot environment. Convergently discovered pre-biotic initial superfamilies with wonderfold structures could have served as a seed for subsequent biological evolution involving gene duplications and divergence.

Entropy capacity determines protein folding

Search and study of the general principles that govern kinetics and thermodynamics of protein folding generate a new insight into the factors controlling this process. Here, based on the known experimental data and using theoretical modeling of protein folding, we demonstrate that there exists an optimal relationship between the average conformational entropy and the average energy of contacts per residue-that is, an entropy capacity-for fast protein folding. Statistical analysis of conformational entropy and number of contacts per residue for 5829 protein structures from four general structural classes (all-alpha, all-beta, alpha/beta, alpha+beta) demonstrates that each class of proteins has its own class-specific average number of contacts (class alpha/beta has the largest number of contacts) and average conformational entropy per residue (class all-alpha has the largest number of rotatable angles phi, psi, and chi per residue). These class-specific features determine the folding rates: alpha proteins are the fastest folding proteins, then follow beta and alpha+beta proteins, and finally alpha/beta proteins are the slowest ones. Our result is in agreement with the experimental folding rates for 60 proteins. This suggests that structural and sequence properties are important determinants of protein folding rates.

Design of HIV-1-PR inhibitors that do not create resistance: blocking the folding of single monomers

The main problems found in designing drugs are those of optimizing the drug-target interaction and of avoiding the insurgence of resistance. We suggest a scheme for the design of inhibitors that can be used as leads for the development of a drug and that do not face either of these problems, and then apply it to the case of HIV-1-PR. It is based on the knowledge that the folding of single-domain proteins, such as each of the monomers forming the HIV-1-PR homodimer, is controlled by local elementary structures (LES), stabilized by local contacts among hydrophobic, strongly interacting, and highly conserved amino acids that play a central role in the folding process. Because LES have evolved over many generations to recognize and strongly interact with each other so as to make the protein fold fast and avoid aggregation with other proteins, highly specific (and thus little toxic) as well as effective folding-inhibitor molecules suggest themselves: short peptides (or eventually their mimetic molecules) displaying the same amino acid sequence of that of LES (p-LES). Aside from being specific and efficient, these inhibitors are expected not to induce resistance; in fact, mutations in HIV-1-PR that successfully avoid the action of p-LES imply the destabilization of one or more LES and thus should lead to protein denaturation. Making use of Monte Carlo simulations, we first identify the LES of the HIV-1-PR and then show that the corresponding p-LES peptides act as effective inhibitors of the folding of the protease.

A folding inhibitor of the HIV-1 protease

Because the human immunodeficiency virus type 1 protease (HIV-1-PR) is an essential enzyme in the viral life cycle, its inhibition can control AIDS. The folding of single-domain proteins, like each of the monomers forming the HIV-1-PR homodimer, is controlled by local elementary structures (LES, folding units stabilized by strongly interacting, highly conserved, as a rule hydrophobic, amino acids). These LES have evolved over myriad generations to recognize and strongly attract each other, so as to make the protein fold fast and be stable in its native conformation. Consequently, peptides displaying a sequence identical to those segments of the monomers associated with LES are expected to act as competitive inhibitors and thus destabilize the native structure of the enzyme. These inhibitors are unlikely to lead to escape mutants as they bind to the protease monomers through highly conserved amino acids, which play an essential role in the folding process. The properties of one of the most promising inhibitors of the folding of the HIV-1-PR monomers found among these peptides are demonstrated with the help of spectrophotometric assays and circular dichroism spectroscopy.

Effects of disulfide bonds on folding behavior and mechanism of the beta-sheet protein tendamistat

Tendamistat, a small disulfide-bonded beta-sheet protein, and its three single/double-disulfide mutants are investigated by using a modified Gō-like model, aiming to understand the folding mechanism of disulfide-bonded protein as well as the effects of removal of disulfide bond on the folding process. Our simulations show that tendamistat and its two single-disulfide mutants are all two-state folders, consistent with the experimental observations. It is found that the disulfide bonds as well as three hydrogen bonds between the N-terminal loop-0 and strand-6 are of significant importance for the folding of tendamistat. Without these interactions, their two-state behaviors become unstable and the predictions of the model are inconsistent with experiments. In addition, the effect of disulfide bonds on the folding process are studied by comparing the wild-type tendamistat and its two mutants; it is found that the removal of either of the C11-C27 or C45-C73 disulfide bond leads to a large decrease in the thermodynamical stability and loss of structure in the unfolded state, and the effect of the former is stronger than that of the later. These simulation results are in good agreement with experiments and, thus, validate our model. Based on the same model, the detailed folding pathways of the wild-type tendamistat and two mutants are studied, and the effect of disulfide bonds on the folding kinetics are discussed. The obtained results provide a detailed folding picture of these proteins and complement experimental findings. Finally, the folding nuclei predicted to be existent in this protein tendamistat as well as its mutants are firstly identified in this work. The positions of the nucleus are consistent with those argued in experimental studies. Therefore, a nucleation/growth folding mechanism that can explain the two-state folding manner is clearly characterized. Moreover, the effect by the removal of each disulfide bond on the folding thermodynamics and dynamics can also be well interpreted from their influence on the folding nucleus. The implementation of this work indicates that the modified Gō-like model really describes the folding behavior of protein tendamistat and could be used to study the folding of other disulfide-bonded proteins.

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.

A critical assessment of the topomer search model of protein folding using a continuum explicit-chain model with extensive conformational sampling

Recently, a series of closely related theoretical constructs termed the "topomer search model" (TSM) has been proposed for the folding mechanism of small, single-domain proteins. A basic assumption of the proposed scenarios is that the rate-limiting step in folding is an essentially unbiased, diffusive search for a conformational state called the native topomer defined by an overall native-like topological pattern. Successes in correlating TSM-predicted folding rates with that of real proteins have been interpreted as experimental support for the model. To better delineate the physics entailed, key TSM concepts are examined here using extensive Langevin dynamics simulations of continuum C(alpha) chain models. The theoretical native topomers of four experimentally well-studied two-state proteins are characterized. Consistent with the TSM perspective, we found that the sizes of the native topomers increase with experimental folding rate. However, a careful determination of the corresponding probabilities that the native topomers are populated during a random search fails to reproduce the previously predicted folding rates. Instead, our results indicate that an unbiased TSM search for the native topomer amounts to a Levinthal-like process that would take an impossibly long average time to complete. Furthermore, intraprotein contacts in all four native topomers considered exhibit no apparent correlation with the experimental phi-values determined from the folding kinetics of these proteins. Thus, the present findings suggest that certain basic, generic yet essential energetic features in protein folding are not accounted for by TSM scenarios to date.

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.

Unexpectedly strong energy stabilization inside the hydrophobic core of small protein rubredoxin mediated by aromatic residues: correlated ab initio quantum chemical calculations

The formation of a hydrophobic core of globular proteins is believed to be the consequence of exterior hydrophobic forces of entropic nature. This, together with the low occurrence of hydrogen bonds in the protein core, leads to the opinion that the energy contribution of core formation to protein folding and stability is negligible. We show that stabilization inside the hydrophobic core of a small protein, rubredoxin, determined by means of high-level correlated ab initio calculations (complete basis set limit of MP2 stabilization energy + CCSD(T) correction term), amounted to approximately 50 kcal/mol. These results clearly demonstrate strong attraction inside a hydrophobic core. This finding may lead to substantial changes in the current view of protein folding. We also point out the inability of the DFT/B3LYP method to describe a strong attraction between studied amino acids.

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.

Extending the folding nucleus of ubiquitin with an independently folding beta-hairpin finger: hurdles to rapid folding arising from the stabilisation of local interactions

The N-terminal beta-hairpin sequence of ubiquitin has been implicated as a folding nucleation site. To extend and stabilise the ubiquitin folding nucleus, we have inserted an autonomously folding 14-residue peptide sequence beta4 which in isolation forms a highly populated beta-hairpin (>70%) stabilised by local interactions. NMR structural analysis of the ubiquitin mutant (Ubeta4) shows that the hairpin finger is fully structured and stabilises ubiquitin by approximately 8kJmol(-1). Protein engineering and kinetic (phi(F)-value) analysis of a series of Ubeta4 mutants shows that the hairpin extension of Ubeta4 is also significantly populated in the transition state (phi(F)-values >0.7) and has the effect of templating the formation of native contacts in the folding nucleus of ubiquitin. However, at low denaturant concentrations the chevron plot of Ubeta4 shows a small deviation from linearity (roll-over effect), indicative of the population of a compact collapsed state, which appears to arise from over-stabilisation of local interactions. Destabilising mutations within the native hairpin sequence and within the engineered hairpin extension, but not elsewhere, eliminate this non-linearity and restore apparent two-state behaviour. The pitfall to stabilising local interactions is to present hurdles to the rapid and efficient folding of small proteins down a smooth folding funnel by trapping partially folded or misfolded states that must unfold or rearrange before refolding.

Multiple folding pathways of the SH3 domain

Experimental observations suggest that proteins follow different folding pathways under different environmental conditions. We perform molecular dynamics simulations of a model of the c-Crk SH3 domain over a broad range of temperatures, and identify distinct pathways in the folding transition. We determine the kinetic partition temperature-the temperature for which the c-Crk SH3 domain undergoes a rapid folding transition with minimal kinetic barriers-and observe that below this temperature the model protein may undergo a folding transition by multiple folding pathways via only one or two intermediates. Our findings suggest the hypothesis that the SH3 domain, a protein fold for which only two-state folding kinetics was observed in previous experiments, may exhibit intermediate states under conditions that strongly stabilize the native state.

Folding Trp-cage to NMR resolution native structure using a coarse-grained protein model

We develop a coarse-grained protein model with a simplified amino acid interaction potential. Using this model, we perform discrete molecular dynamics folding simulations of a small 20-residue protein--Trp-cage--from a fully extended conformation. We demonstrate the ability of the Trp-cage model to consistently reach conformations within 2-angstroms backbone root-mean-square distance from the corresponding NMR structures. The minimum root-mean-square distance of Trp-cage conformations in simulations can be <1 angstroms. Our findings suggest that, at least in the case of Trp-cage, a detailed all-atom protein model with a molecular mechanics force field is not necessary to reach the native state of a protein. Our results also suggest that the success of folding Trp-cage in our simulations and in the reported all-atom molecular mechanics simulation studies may be mainly due to the special stabilizing features specific to this miniprotein.

Identification of residual structure within denatured antichymotrypsin: implications for serpin folding and misfolding

The native serpin fold is metastable and possesses the inherent ability to convert into more stable, but inactive, conformations. In order to understand why serpins attain the native fold instead of other more thermodynamically favourable folds we have investigated the presence of residual structure within denatured antichymotrypsin (ACT). Through mutagenesis we created a single tryptophan variant of ACT in which a Trp residue (276) is situated on the H-helix, located within a region known as the B/C barrel. The presence of residual structure around Trp 276 in 5 M guanidine hydrochloride (GdnHCl) was shown by fluorescence and circular dichroism spectroscopy and fluorescence lifetime experiments. The residual structure was disrupted in the presence of 5 M guanidine thiocyanate (GdnSCN). Protein refolding studies showed that significant refolding could be achieved from the GdnHCl denatured state but not the GdnSCN denatured form. The implications of these data on the folding and misfolding of the serpin superfamily are discussed.

Protein folding rates correlate with heterogeneity of folding mechanism

By observing trends in the folding kinetics of experimental 2-state proteins at their transition midpoints, and by observing trends in the barrier heights of numerous simulations of coarse-grained, C(alpha) model Go proteins, we show that folding rates correlate with the degree of heterogeneity in the formation of native contacts. Statistically significant correlations are observed between folding rates and measures of heterogeneity inherent in the native topology, as well as between rates and the variance in the distribution of either experimentally measured or simulated phi values.

Intermediates and the folding of proteins L and G

We use a minimalist protein model, in combination with a sequence design strategy, to determine differences in primary structure for proteins L and G, which are responsible for the two proteins folding through distinctly different folding mechanisms. We find that the folding of proteins L and G are consistent with a nucleation-condensation mechanism, each of which is described as helix-assisted beta-1 and beta-2 hairpin formation, respectively. We determine that the model for protein G exhibits an early intermediate that precedes the rate-limiting barrier of folding, and which draws together misaligned secondary structure elements that are stabilized by hydrophobic core contacts involving the third beta-strand, and presages the later transition state in which the correct strand alignment of these same secondary structure elements is restored. Finally, the validity of the targeted intermediate ensemble for protein G was analyzed by fitting the kinetic data to a two-step first-order reversible reaction, proving that protein G folding involves an on-pathway early intermediate, and should be populated and therefore observable by experiment.

Kinetics of protein-DNA interaction: facilitated target location in sequence-dependent potential

Recognition and binding of specific sites on DNA by proteins is central for many cellular functions such as transcription, replication, and recombination. In the process of recognition, a protein rapidly searches for its specific site on a long DNA molecule and then strongly binds this site. Here we aim to find a mechanism that can provide both a fast search (1-10 s) and high stability of the specific protein-DNA complex (Kd=10(-15)-10(-8) M). Earlier studies have suggested that rapid search involves sliding of the protein along the DNA. Here we consider sliding as a one-dimensional diffusion in a sequence-dependent rough energy landscape. We demonstrate that, despite the landscape's roughness, rapid search can be achieved if one-dimensional sliding is accompanied by three-dimensional diffusion. We estimate the range of the specific and nonspecific DNA-binding energy required for rapid search and suggest experiments that can test our mechanism. We show that optimal search requires a protein to spend half of its time sliding along the DNA and the other half diffusing in three dimensions. We also establish that, paradoxically, realistic energy functions cannot provide both rapid search and strong binding of a rigid protein. To reconcile these two fundamental requirements we propose a search-and-fold mechanism that involves the coupling of protein binding and partial protein folding. The proposed mechanism has several important biological implications for search in the presence of other proteins and nucleosomes, simultaneous search by several proteins, etc. The proposed mechanism also provides a new framework for interpretation of experimental and structural data on protein-DNA interactions.

What can one learn from experiments about the elusive transition state?

We present the results of an exact analysis of a model energy landscape of a protein to clarify the idea of the transition state and the physical meaning of the phi values determined in protein engineering experiments. We benchmark our findings to various theoretical approaches proposed in the literature for the identification and characterization of the transition state.

Understanding the determinants of stability and folding of small globular proteins from their energetics

The results of minimal model calculations indicate that the stability and the kinetic accessibility of the native state of small globular proteins are controlled by few "hot" sites. By means of molecular dynamics simulations around the native conformation, which describe the protein and the surrounding solvent at the all-atom level, an accurate and compact energetic map of the native state of the protein is generated. This map is further simplified by means of an eigenvalue decomposition. The components of the eigenvector associated with the lowest eigenvalue indicate which hot sites are likely to be responsible for the stability and for the rapid folding of the protein. The comparison of the results of the model with the findings of mutagenesis experiments performed for four small proteins show that the eigenvalue decomposition method is able to identify between 60% and 80% of these (hot) sites.

Importance of hydrophobic cluster formation through long-range contacts in the folding transition state of two-state proteins

Understanding the folding pathways of proteins is a challenging task. The Phi value approach provides a detailed understanding of transition-state structures of folded proteins. In this work, we have computed the hydrophobicity associated with each residue in the folded state of 16 two-state proteins and compared the Phi values of each mutant residue. We found that most of the residues with high Phi value coincide with local maximum in surrounding hydrophobicity, or have nearby residues that show such maximum in hydrophobicity, indicating the importance of hydrophobic interactions in the transition state. We have tested our approach to different structural classes of proteins, such as alpha-helical, SH3 domains of all-beta proteins, beta-sandwich, and alpha/beta proteins, and we observed a good agreement with experimental results. Further, we have proposed a hydrophobic contact network pattern to relate the Phi values with long-range contacts, which will be helpful to understand the transition-state structures of folded proteins. The present approach could be used to identify potential hydrophobic clusters that may form through long-range contacts during the transition state.

Scattered Hammond plots reveal second level of site-specific information in protein folding: phi' (beta++)

Site-specific information about structural heterogeneities of the protein-folding transition-state ensemble is commonly derived from the scatter of the Brønsted plot through the individual values of = Delta logk(f)/Delta logK(D-N). Here, we provide a second level of site-specific detail in the transition-state analysis by demonstrating that the scatter of the Hammond plot is related to heterogeneities in the -value growth. That is, the extent of transition-state movement (Delta beta(++)) is proportional to the free-energy gradient of the mutational perturbation across the top of the activation barrier, '(beta(++)) proportional, variant Delta logK(D-N). The analysis is applied to the two-state protein L23 where the site-specific free-energy gradients are used to identify the interactions that show the highest degree of consolidation after crossing the barrier top. These interactions are distributed as a shell around the high- initiation point and denote the side-chain contacts that add criticality to the folding nucleus.

Simulation, experiment, and evolution: understanding nucleation in protein S6 folding

In this study, we explore nucleation and the transition state ensemble of the ribosomal protein S6 using a Monte Carlo (MC) Go model in conjunction with restraints from experiment. The results are analyzed in the context of extensive experimental and evolutionary data. The roles of individual residues in the folding nucleus are identified, and the order of events in the S6 folding mechanism is explored in detail. Interpretation of our results agrees with, and extends the utility of, experiments that shift phi-values by modulating denaturant concentration and presents strong evidence for the realism of the mechanistic details in our MC Go model and the structural interpretation of experimental phi-values. We also observe plasticity in the contacts of the hydrophobic core that support the specific nucleus. For S6, which binds to RNA and protein after folding, this plasticity may result from the conformational flexibility required to achieve biological function. These results present a theoretical and conceptual picture that is relevant in understanding the mechanism of nucleation in protein folding.

Folding probabilities: A novel approach to folding transitions and the two-dimensional Ising-model

The theoretical concept of folding probability, p(fold), has proven to be a useful means to characterize the kinetics of protein folding. Here, we illustrate the practical importance of p(fold) and demonstrate how it can be determined theoretically. We derive a general analytical expression for p(fold) and show how it can be estimated from simulations for systems where the transition rates between the relevant microstates are not known. By analyzing the Ising model we are able to determine the scaling behavior of the numerical error in the p(fold) estimate as function of the number of analyzed Monte Carlo runs. We apply our method to a simple, newly developed protein folding model for the formation of alpha helices. It is demonstrated that our technique highly parallelizes the calculation of p(fold) and that it is orders of magnitude more efficient than conventional approaches.

Outlining folding nuclei in globular proteins

Our theoretical approach for prediction of folding/unfolding nuclei in three-dimensional protein structures is based on a search for free energy saddle points on networks of protein unfolding pathways. Under some approximations, this search is performed rapidly by dynamic programming and results in prediction of Phi values, which can be compared with those found experimentally. In this study, we optimize some details of the model (specifically, hydrogen atoms are taken into account in addition to heavy atoms), and compare the theoretically obtained and experimental Phi values (which characterize involvement of residues in folding nuclei) for all 17 proteins, where Phi values are now known for many residues. We show that the model provides good Phi value predictions for proteins whose structures have been determined by X-ray analysis (the average correlation coefficient is 0.65), with a more limited success for proteins whose structures have been determined by NMR techniques only (the average correlation coefficient is 0.34), and that the transition state free energies computed from the same model are in a good anticorrelation with logarithms of experimentally measured folding rates at mid-transition (the correlation coefficient is -0.73).

Imprint of evolution on protein structures

We attempt to understand the evolutionary origin of protein folds by simulating their divergent evolution with a three-dimensional lattice model. Starting from an initial seed lattice structure, evolution of model proteins progresses by sequence duplication and subsequent point mutations. A new gene's ability to fold into a stable and unique structure is tested each time through direct kinetic folding simulations. Where possible, the algorithm accepts the new sequence and structure and thus a "new protein structure" is born. During the course of each run, this model evolutionary algorithm provides several thousand new proteins with diverse structures. Analysis of evolved structures shows that later evolved structures are more designable than seed structures as judged by recently developed structural determinant of protein designability, as well as direct estimate of designability for selected structures by thermodynamic sampling of their sequence space. We test the significance of this trend predicted on lattice models on real proteins and show that protein domains that are found in eukaryotic organisms only feature statistically significant higher designability than their prokaryotic counterparts. These results present a fundamental view on protein evolution highlighting the relative roles of structural selection and evolutionary dynamics on genesis of modern proteins.

Discerning the Structure and Energy of Multiple Transition States in Protein Folding using ψ-Analysis

We quantify the degree to which folding occurs along a complex landscape with structurally distinct pathways using psi-analysis in combination with a protein engineering method that identifies native, non-covalent polypeptide interactions and their relative populations at the rate-limiting step. By probing the proximity of two specific partners, this method is extremely well-suited for comparison to theoretical simulations. Using ubiquitin as a model system, we detect individual pathways with site-resolved resolution, demonstrating that the protein folds through a native-like transition state ensemble with a common nucleus that contains heterogeneous features on its periphery. The consensus transition state topology has part of the major helix docked against four properly aligned beta-strands. However, structural heterogeneity exists in the transition state ensemble, wherein peripheral regions are differentially populated according to their relative stability. Pathway diversity reflects the variable order of formation of these peripheral elements, which radiate outward from the common nucleus. These results, which show only moderate agreement with traditional mutational phi-analysis, provide an extraordinarily detailed and quantitative description of protein folding.

The importance of explicit chain representation in protein folding models: an examination of Ising-like models

A class of models that represents a protein chain as a sequence of "folded" and "unfolded" residues has recently been used to correlate rates and mechanisms of protein folding with the protein native structure. In order to better understand the conditions under which these "Ising-like" models apply, we compare results from this model to those obtained from an off-lattice model which uses the same potential function. We find that Ising-like models by construction impose folding via a highly sequential nucleation-condensation mechanism, which in turn leads to more rugged energy landscapes, fewer "pathways" to the native state, and in the specific case examined here, the cold shock protein A from Escherichia coli, a qualitative difference in the most likely order of events in folding.

Uncovering Network Systems Within Protein Structures

Traditionally, proteins have been viewed as a construct based on elements of secondary structure and their arrangement in three-dimensional space. In a departure from this perspective we show that protein structures can be modelled as network systems that exhibit small-world, single-scale, and to some degree, scale-free properties. The phenomenological network concept of degrees of separation is applied to three-dimensional protein structure networks and reveals how amino acid residues can be connected to each other within six degrees of separation. This work also illuminates the unique features of protein networks in comparison to other networks currently studied. Recognising that proteins are networks provides a means of rationalising the robustness in the overall three-dimensional fold of a protein against random mutations and suggests an alternative avenue to investigate the determinants of protein structure, function and folding.

Folding of Cu, Zn Superoxide Dismutase and Familial Amyotrophic Lateral Sclerosis

Cu, Zn superoxide dismutase (SOD1) has been implicated in the familial form of the neurodegenerative disease amyotrophic lateral sclerosis (ALS). It has been suggested that mutant mediated SOD1 misfolding/aggregation is an integral part of the pathology of ALS. We study the folding thermodynamics and kinetics of SOD1 using a hybrid molecular dynamics approach. We reproduce the experimentally observed SOD1 folding thermodynamics and find that the residues which contribute the most to SOD1 thermal stability are also crucial for apparent two-state folding kinetics. Surprisingly, we find that these residues are located on the surface of the protein and not in the hydrophobic core. Mutations in some of the identified residues are found in patients with the disease. We argue that the identified residues may play an important role in aggregation. To further characterize the folding of SOD1, we study the role of cysteine residues in folding and find that non-native disulfide bond formation may significantly alter SOD1 folding dynamics and aggregation propensity.

Structural determinants of the rate of protein folding

To understand the mechanism of protein folding and to assist rational design of fast-folding, non-aggregating and stable artificial enzymes, it is essential to determine the structural parameters which govern the rate constants of folding, kf. It has been found that -logkf is a linear function of the so-called chain topology parameter (CTP) within the range of 10(-1)s(-1)< or = kf < or =10(8)s(-1). The correlation between -logkf and CTP is much improved than using previously published contact order (CO) method. It has been further suggested that short sequence separations may be preferred for the establishment of stable interactions for the design of novel artificial enzymes and the modification of slow-folding proteins with aggregating intermediates.

Side-chain dynamics and protein folding

The processes by which protein side chains reach equilibrium during a folding reaction are investigated using both lattice and all-atom simulations. We find that rates of side-chain relaxation exhibit a distribution over the protein structure, with the fastest relaxing side chains located in positions kinetically important for folding. Traversal of the major folding transition state corresponds to the freezing of a small number of side chains, belonging to the folding nucleus, whereas the rest of the protein proceeds toward equilibrium via backbone fluctuations around the native fold. The postnucleation processes by which side chains relax are characterized by very slow dynamics and many barrier crossings, and thus resemble the behavior of a glass.

Master equation approach to finding the rate-limiting steps in biopolymer folding

A master equation approach is developed to find the rate-limiting steps in biopolymer folding, where the folding kinetics is described as a linear combination of basic kinetic modes determined from the eigenvalues and eigenvectors of the rate matrix. Because the passage of a rate-limiting step is intrinsically related to the folding speed, it is possible to probe and to identify the rate-limiting steps through the folding from different unfolded initial conformations. In a master equation approach, slow and fast folding speeds are directly correlated to the large and small contributions of the (rate-limiting) slow kinetic modes. Because the contributions from the slow modes can be computed from the corresponding eigenvectors, the rate-limiting steps can be identified from the eigenvectors of the slow modes. Our rate-limiting searching method has been tested for a simplified hairpin folding kinetics model, and it may provide a general transition state searching method for biopolymer folding.