Posttransition state desolvation of the hydrophobic core of the src-SH3 protein domain

Conserved patterns and interactions in the unfolding transition state across SH3 domain structural homologues

Proteins with similar structures are generally assumed to arise from similar sequences. However, there are more cases than not where this is not true. The dogma is that sequence determines structure; how, then, can very different sequences fold to the same structure? Here, we employ high temperature unfolding simulations to probe the pathways and specific interactions that direct the folding and unfolding of the SH3 domain. The SH3 metafold in the Dynameomics Database consists of 753 proteins with the same structure, but varied sequences and functions. To investigate the relationship between sequence and structure, we selected 17 targets from the SH3 metafold with high sequence variability. Six unfolding simulations were performed for each target, transition states were identified, revealing two general folding/unfolding pathways at the transition state. Transition states were also expressed as mathematical graphs of connected chemical nodes, and it was found that three positions within the structure, independent of sequence, were consistently more connected within the graph than any other nearby positions in the sequence. These positions represent a hub connecting different portions of the structure. Multiple sequence alignment and covariation analyses also revealed certain positions that were more conserved due to packing constraints and stabilizing long-range contacts. This study demonstrates that members of the SH3 domain with different sequences can unfold through two main pathways, but certain characteristics are conserved regardless of the sequence or unfolding pathway. While sequence determines structure, we show that disparate sequences can provide similar interactions that influence folding and lead to similar structures.

Thermodynamics of Protein Folding Studied by Umbrella Sampling along a Reaction Coordinate of Native Contacts

Spontaneous transitions between the native and non-native protein conformations are normally rare events that hardly take place in typical unbiased molecular dynamics simulations. It was recently demonstrated that such transitions can be well described by a reaction coordinate, Q, that represents the collective fraction of the native contacts between the protein atoms. Here we attempt to use this reaction coordinate to enhance the conformational sampling. We perform umbrella sampling simulations with biasing potentials on Q for two model proteins, Trp-Cage and BBA, using the CHARMM force field. Hamiltonian replica exchange is implemented in these simulations to further facilitate the sampling. The simulations appear to have reached satisfactory convergence, resulting in unbiased free energies as a function of Q. In addition to the native structure, multiple folded conformations are identified in the reconstructed equilibrium ensemble. Some conformations without any native contacts nonetheless have rather compact geometries and are stabilized by hydrogen bonds not present in the native structure. Whereas the enhanced sampling along Q reasonably reproduces the equilibrium conformational space, we also find that the folding of an α-helix in Trp-Cage is a slow degree of freedom orthogonal to Q and therefore cannot be accelerated by biasing the reaction coordinate. Overall, we conclude that whereas Q is an excellent parameter to analyze the simulations, it is not necessarily a perfect reaction coordinate for enhanced sampling, and better incorporation of other slow degrees of freedom may further improve this reaction coordinate.

Crystallographic studies on protein misfolding: Domain swapping and amyloid formation in the SH3 domain

Oligomerization by 3D domain swapping is found in a variety of proteins of diverse size, fold and function. In the early 1960s this phenomenon was postulated for the oligomers of ribonuclease A, but it was not until the 1990s that X-ray diffraction provided the first experimental evidence of this special manner of oligomerization. Nowadays, structural information has allowed the identification of these swapped oligomers in over one hundred proteins. Although the functional relevance of this phenomenon is not clear, this alternative folding of protomers into intertwined oligomers has been related to amyloid formation. Studies on proteins that develop 3D domain swapping might provide some clues on the early stages of amyloid formation. The SH3 domain is a small modular domain that has been used as a model to study the basis of protein folding. Among SH3 domains, the c-Src-SH3 domain emerges as a helpful model to study 3D domain swapping and amyloid formation.

Electrostatic effects in the folding of the SH3 domain of the c-Src tyrosine kinase: pH-dependence in 3D-domain swapping and amyloid formation

The SH3 domain of the c-Src tyrosine kinase (c-Src-SH3) aggregates to form intertwined dimers and amyloid fibrils at mild acid pHs. In this work, we show that a single mutation of residue Gln128 of this SH3 domain has a significant effect on: (i) its thermal stability; and (ii) its propensity to form amyloid fibrils. The Gln128Glu mutant forms amyloid fibrils at neutral pH but not at mild acid pH, while Gln128Lys and Gln128Arg mutants do not form these aggregates under any of the conditions assayed. We have also solved the crystallographic structures of the wild-type (WT) and Gln128Glu, Gln128Lys and Gln128Arg mutants from crystals obtained at different pHs. At pH 5.0, crystals belong to the hexagonal space group P6₅22 and the asymmetric unit is formed by one chain of the protomer of the c-Src-SH3 domain in an open conformation. At pH 7.0, crystals belong to the orthorhombic space group P2₁2₁2₁, with two molecules at the asymmetric unit showing the characteristic fold of the SH3 domain. Analysis of these crystallographic structures shows that the residue at position 128 is connected to Glu106 at the diverging β-turn through a cluster of water molecules. Changes in this hydrogen-bond network lead to the displacement of the c-Src-SH3 distal loop, resulting also in conformational changes of Leu100 that might be related to the binding of proline rich motifs. Our findings show that electrostatic interactions and solvation of residues close to the folding nucleation site of the c-Src-SH3 domain might play an important role during the folding reaction and the amyloid fibril formation.

Mapping the Native Conformational Ensemble of Proteins from a Combination of Simulations and Experiments: New Insight into the src-SH3 Domain

The biological function of a protein is strongly tied to the ensemble of three-dimensional conformations populated at physiological temperature, and dynamically transforming into each other. Experimental techniques such as nuclear magnetic resonance spectroscopy (NMR) provide a wealth of structural and dynamical information, which, in combination with an accurate atomic-level computational modeling, can disclose the details of protein behavior. We here propose a fast and efficient protocol employing molecular dynamics (MD) simulations and NMR chemical shifts, which allows one to reconstruct the detailed conformational ensemble of small globular proteins. In the case of the well-studied src-SH3 domain, we are able to obtain new important insight including the existence of a helical state in the RT loop and a pathway for single-file water diffusion in and out of the core.

Hydrated and dehydrated tertiary interactions--opening and closing--of a four-helix bundle peptide

The structural heterogeneity and thermal denaturation of a dansyl-labeled four-helix bundle homodimeric peptide was studied with steady-state and time-resolved fluorescence spectroscopy and with circular dichroism (CD). At room temperature the fluorescence decay of the polarity-sensitive dansyl, located in the hydrophobic core region, can be described by a broad distribution of fluorescence lifetimes, reflecting the heterogeneous microenvironment. However, the lifetime distribution is nearly bimodal, which we ascribe to the presence of two major conformational subgroups. Since the fluorescence lifetime reflects the water content of the four-helix bundle conformations, we can use the lifetime analysis to monitor the change in hydration state of the hydrophobic core of the four-helix bundle. Increasing the temperature from 9 degrees C to 23 degrees C leads to an increased population of molten-globule-like conformations with a less ordered helical backbone structure. The fluorescence emission maximum remains constant in this temperature interval, and the hydrophobic core is not strongly affected. Above 30 degrees C the structural dynamics involve transient openings of the four-helix bundle structure, as evidenced by the emergence of a water-quenched component and less negative CD. Above 60 degrees C the homodimer starts to dissociate, as shown by the increasing loss of CD and narrow, short-lived fluorescence lifetime distributions.

Probing the role of hydration in the unfolding transitions of carbonmonoxy myoglobin and apomyoglobin

We show that the equilibrium unfolding transition of horse carbonmonoxy myoglobin monitored by the stretching vibration of the CO ligand, a local environmental probe, is very sharp and, thus, quite different from those measured by global conformational reporters. In addition, the denatured protein exhibits an A(0)-like CO band. We hypothesize that this sharp transition reports penetration of water into the heme pocket of the protein. Parallel experiments on horse apomyoglobin, wherein an environment-sensitive fluorescent probe, nile red, was used, also reveals a similar putative hydration event. Given the importance of dehydration in protein folding and also the recent debate over the interpretation of probe-dependent unfolding transitions, these results have strong implications on the mechanism of protein folding.

Relaxation dispersion NMR spectroscopy as a tool for detailed studies of protein folding

Characterization of the mechanisms by which proteins fold into their native conformations is important not only for protein structure prediction and design but also because protein misfolding intermediates may play critical roles in fibril formation that are commonplace in neurodegenerative disorders. In practice, the study of folding pathways is complicated by the fact that for the most part intermediates are low-populated and short-lived so that biophysical studies are difficult. Due to recent methodological advances, relaxation dispersion NMR spectroscopy has emerged as a particularly powerful tool to obtain high-resolution structural information about protein folding events on the millisecond timescale. Applications of the methodology to study the folding of SH3 domains have shown that folding proceeds via previously undetected on-pathway intermediates, sometimes stabilized by nonnative long-range interactions. The relaxation dispersion approach provides a detailed kinetic and thermodynamic description of the folding process as well as the promise of obtaining an atomic level structural description of intermediate states. We review the concerted application of a variety of recently developed NMR relaxation dispersion experiments to obtain a "high-resolution" picture of the folding pathway of the A39V/N53P/V55L Fyn SH3 domain.

Role of water in mediating the assembly of Alzheimer amyloid-beta Abeta16-22 protofilaments

The role of water in promoting the formation of protofilaments (the basic building blocks of amyloid fibrils) is investigated using fully atomic molecular dynamics simulations. Our model protofilament consists of two parallel beta-sheets of Alzheimer Amyloid-beta 16-22 peptides (Ac-K(16)-L(17)-V(18)-F(19)-F(20)-A(21)-E(22)-NH2). Each sheet presents a distinct hydrophobic and hydrophilic face and together self-assemble to a stable protofilament with a core consisting of purely hydrophobic residues (L(17), F(19), A(21)), with the two charged residues (K(16), E(22)) pointing to the solvent. Our simulations reveal a subtle interplay between a water mediated assembly and one driven by favorable energetic interactions between specific residues forming the interior of the protofilament. A dewetting transition, in which water expulsion precedes hydrophobic collapse, is observed for some, but not all molecular dynamics trajectories. In the trajectories in which no dewetting is observed, water expulsion and hydrophobic collapse occur simultaneously, with protofilament assembly driven by direct interactions between the hydrophobic side chains of the peptides (particularly between F-F residues). For those same trajectories, a small increase in the temperature of the simulation (on the order of 20 K) or a modest reduction in the peptide-water van der Waals attraction (on the order of 10%) is sufficient to induce a dewetting transition, suggesting that the existence of a dewetting transition in simulation might be sensitive to the details of the force field parametrization.

Hydrophobicity of protein surfaces: Separating geometry from chemistry

To better understand the role of surface chemical heterogeneity in natural nanoscale hydration, we study via molecular dynamics simulation the structure and thermodynamics of water confined between two protein-like surfaces. Each surface is constructed to have interactions with water corresponding to those of the putative hydrophobic surface of a melittin dimer, but is flattened rather than having its native "cupped" configuration. Furthermore, peripheral charged groups are removed. Thus, the role of a rough surface topography is removed, and results can be productively compared with those previously observed for idealized, atomically smooth hydrophilic and hydrophobic flat surfaces. The results indicate that the protein surface is less hydrophobic than the idealized counterpart. The density and compressibility of water adjacent to a melittin dimer is intermediate between that observed adjacent to idealized hydrophobic or hydrophilic surfaces. We find that solvent evacuation of the hydrophobic gap (cavitation) between dimers is observed when the gap has closed to sterically permit a single water layer. This cavitation occurs at smaller pressures and separations than in the case of idealized hydrophobic flat surfaces. The vapor phase between the melittin dimers occupies a much smaller lateral region than in the case of the idealized surfaces; cavitation is localized in a narrow central region between the dimers, where an apolar amino acid is located. When that amino acid is replaced by a polar residue, cavitation is no longer observed.

Phi-value analysis of a three-state protein folding pathway by NMR relaxation dispersion spectroscopy

Experimental studies of protein folding frequently are consistent with two-state folding kinetics. However, recent NMR relaxation dispersion studies of several fast-folding mutants of the Fyn Src homology 3 (SH3) domain have established that folding proceeds through a low-populated on-pathway intermediate, which could not be detected with stopped-flow experiments. The dispersion experiments provide precise kinetic and thermodynamic parameters that describe the folding pathway, along with a detailed site-specific structural characterization of both the intermediate and unfolded states from the NMR chemical shifts that are extracted. Here we describe NMR relaxation dispersion Phi-value analysis of the A39V/N53P/V55L Fyn SH3 domain, where the effects of suitable point mutations on the energy landscape are quantified, providing additional insight into the structure of the folding intermediate along with per-residue structural information of both rate-limiting transition states that was not available from previous studies. In addition to the advantage of delineating the full three-state folding pathway, the use of NMR relaxation dispersion as opposed to stopped-flow kinetics to quantify Phi values facilitates their interpretation because the obtained chemical shifts monitor any potential structural changes along the folding pathway that might be introduced by mutation, a significant concern in their analysis. Phi-Value analysis of several point mutations of A39V/N53P/V55L Fyn SH3 establishes that the beta(3)-beta(4)-hairpin already is formed in the first transition state, whereas strand beta(1), which forms nonnative interactions in the intermediate, does not fully adopt its native conformation until after the final transition state. The results further support the notion that on-pathway intermediates can be stabilized by nonnative contacts.

Hydration of a hydrophobic cavity and its functional role: a simulation study of human interleukin-1beta

Molecular dynamics simulations reveal that the hydrophobic cavity in human cytokine Interleukin-1beta is hydrated and can dynamically accommodate between one and four water molecules. These waters have residence times >> 500 ps and can give rise to detectable NOEs, in agreement with NMR observations of Ernst et al. (Science 1995; 267:1813-1817). The waters also display high positional disorder within the cavity, which explains why they have not been resolved crystallographically. The average distribution of water molecules over time within the cavity matches well the low resolution electron density extracted by Yu et al. (Proc Natl Acad Sci 1999; 96:103-108). The water molecules hydrate the hydrophobic cavity preferentially as complex clusters. These clusters result from a combination of hydrogen bonds between the waters and stabilizing interactions between the waters and aromatic rings forming the cavity. Free energy estimates suggest that it takes 4-waters to hydrate the cavity in a thermodynamically stable manner leading to a gain in free energy of transfer from bulk of approximately approximately 3.6 kcal/mol. This arises from the existence of the water clusters in multiple hydrogen bonded states. In addition, the waters are found to migrate either individually or as clusters out of the cavity through several pathways. The upper limit for one-dimensional diffusion of the waters within the protein matrix is 4 A/ps (relative to 6 A/ps for bulk). Simulations reveal pathways in addition to those identified crystallographically, with motions controlled by the rotations of sidechains. We find that only when the hydrophobic cavity is hydrated, do correlated motions couple distant sites with the sites that make contact with the receptor and this data partly offers an explanation of experimental mutagenesis data. Simulations, together with recent observations based on mutagenesis by Heidary et al. (J Mol Biol 2005; 353:1187-1198) that hydrogen bond networks couple motions across long distances in interleukin-1beta, lead us to hypothesize that the hydration of the cavity (conserved across mammals) can thermodynamically enhance hydrogen bond networks to enable coupling across long distances by acting as a plug and this in turn enables a kinetic control of the rate of transmission of signals.

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.

Identification of a Collapsed Intermediate with Non-native Long-range Interactions on the Folding Pathway of a Pair of Fyn SH3 Domain Mutants by NMR Relaxation Dispersion Spectroscopy

Recent 15N and 13C spin-relaxation dispersion studies of fast-folding mutants of the Fyn SH3 domain have established that folding proceeds through a low-populated on-pathway intermediate (I) where the central beta-sheet is at least partially formed, but without interactions between the NH2- and COOH-terminal beta-strands that exist in the folded state (F). Initial studies focused on mutants where Gly48 is replaced; in an effort to establish whether this intermediate is a general feature of Fyn SH3 folding a series of 15N relaxation experiments monitoring the folding of Fyn SH3 mutants N53P/V55L and A39V/N53P/V55L are reported here. For these mutants as well, folding proceeds through an on-pathway intermediate with similar features to those observed for G48M and G48V Fyn SH3 domains. However, the 15N chemical shifts extracted for the intermediate indicate pronounced non-native contacts between the NH2 and COOH-terminal regions not observed previously. The kinetic parameters extracted for the folding of A39V/N53P/V55L Fyn SH3 from the three-state folding model F<-->I<-->U are in good agreement with folding and unfolding rates extrapolated to zero denaturant obtained from stopped-flow experiments analyzed in terms of a simplified two-state folding reaction. The folding of the triple mutant was studied over a wide range of temperatures, establishing that there is no difference in heat capacities between F and I states. This confirms a compact folding intermediate structure, which is supported by the 15N chemical shifts of the I state extracted from the dispersion data. The temperature-dependent relaxation data simplifies data analysis because at low temperatures (< 25 degrees C) the unfolded state (U) is negligibly populated relative to I and F. A comparison between parameters extracted at low temperatures where the F<-->I exchange model is appropriate with those from the more complex, three-state model at higher temperatures has been used to validate the protocol for analysis of three-site exchange relaxation data.

Transition-states in protein folding kinetics: the structural interpretation of Phi values

Phi values are experimental measures of the effects of mutations on the folding kinetics of a protein. A central question is what structural information Phi values give about the transition-state of folding. Traditionally, a Phi value is interpreted as representing the "nativeness" of a mutated residue in the transition-state. However, this interpretation is often problematic. We present here a better structural interpretation of Phi values for mutations within a given helix. Our interpretation is based on a simple physical model that distinguishes between secondary and tertiary free energy contributions of helical residues. From a linear fit of the model to experimental data, we obtain two structural parameters: the extent of helix formation in the transition-state, and the nativeness of tertiary interactions in the transition-state. We apply the model to all proteins with well-characterized helices for which more than 10 Phi values are available: protein A, CI2, and protein L. The model is simple to apply to experimental data, captures nonclassical Phi values <0 or >1 in these helices, and explains how different mutations at a given site can lead to different Phi values.

Sampling the multiple folding mechanisms of Trp-cage in explicit solvent

We investigate the kinetic pathways of folding and unfolding of the designed miniprotein Trp- cage in explicit solvent. Straightforward molecular dynamics and replica exchange methods both have severe convergence problems, whereas transition path sampling allows us to sample unbiased dynamical pathways between folded and unfolded states and leads to deeper understanding of the mechanisms of (un)folding. In contrast to previous predictions employing an implicit solvent, we find that Trp-cage folds primarily (80% of the paths) via a pathway forming the tertiary contacts and the salt bridge, before helix formation. The remaining 20% of the paths occur in the opposite order, by first forming the helix. The transition states of the rate-limiting steps are solvated native-like structures. Water expulsion is found to be the last step upon folding for each route. Committor analysis suggests that the dynamics of the solvent is not part of the reaction coordinate. Nevertheless, during the transition, specific water molecules are strongly bound and can play a structural role in the folding.

GroEL-mediated protein folding: making the impossible, possible

Protein folding is a spontaneous process that is essential for life, yet the concentrated and complex interior of a cell is an inherently hostile environment for the efficient folding of many proteins. Some proteins-constrained by sequence, topology, size, and function-simply cannot fold by themselves and are instead prone to misfolding and aggregation. This problem is so deeply entrenched that a specialized family of proteins, known as molecular chaperones, evolved to assist in protein folding. Here we examine one essential class of molecular chaperones, the large, oligomeric, and energy utilizing chaperonins or Hsp60s. The bacterial chaperonin GroEL, along with its co-chaperonin GroES, is probably the best-studied example of this family of protein-folding machine. In this review, we examine some of the general properties of proteins that do not fold well in the absence of GroEL and then consider how folding of these proteins is enhanced by GroEL and GroES. Recent experimental and theoretical studies suggest that chaperonins like GroEL and GroES employ a combination of protein isolation, unfolding, and conformational restriction to drive protein folding under conditions where it is otherwise not possible.

WATER MEDIATION IN PROTEIN FOLDING AND MOLECULAR RECOGNITION

Water is essential for life in many ways, and without it biomolecules might no longer truly be biomolecules. In particular, water is important to the structure, stability, dynamics, and function of biological macromolecules. In protein folding, water mediates the collapse of the chain and the search for the native topology through a funneled energy landscape. Water actively participates in molecular recognition by mediating the interactions between binding partners and contributes to either enthalpic or entropic stabilization. Accordingly, water must be included in recognition and structure prediction codes to capture specificity. Thus water should not be treated as an inert environment, but rather as an integral and active component of biomolecular systems, where it has both dynamic and structural roles. Focusing on water sheds light on the physics and function of biological machinery and self-assembly and may advance our understanding of the natural design of proteins and nucleic acids.

Hydration of the folding transition state ensemble of a protein

A complete description of the mechanisms of protein folding requires knowledge of the structural and physical character of the folding transition state ensembles (TSEs). A key question concerning the role of hydration of the hydrophobic core in determining folding mechanisms remains. To address this, we probed the state of hydration of the TSE of staphylococcal nuclease (SNase) by examining the fluorescence-detected pressure-jump relaxation behavior of six SNase variants in which a residue in the hydrophobic core, Val-66, was replaced with polar or ionizable residues (Lys, Arg, His, Asp, Glu, and Asn). Because of a large positive activation volume for folding, the major effect of pressure on the wild-type protein is to decrease the folding rate. By the time wild-type SNase reaches the folding transition state, most water has already been expelled from its hydrophobic core. In contrast, the major effect of pressure on the variant proteins is an increase in the unfolding rate due to a large negative activation volume for unfolding. This results from a significant increase in the level of hydration of the TSE when an internal ionizable group is present. These data confirm that the role of water in the folding reaction can differ from protein to protein and that even a single substitution in a critical position can modulate significantly the properties of the TSE.

Solvation of the folding-transition state in Pseudomonas aeruginosa azurin is modulated by metal: Solvation of azurin's folding nucleus

The role of water in protein folding, specifically its presence or not in the transition-state structure, is an unsolved question. There are two common classes of folding-transition states: diffuse transition states, in which almost all side chains have similar, rather low phi (phi) values, and polarized transition states, which instead display distinct substructures with very high phi-values. Apo-and zinc-forms of Pseudomonas aeruginosa azurin both fold in two-state equilibrium and kinetic reactions; while the apo-form exhibits a polarized transition state, the zinc form entails a diffuse, moving transition state. To examine the presence of water in these two types of folding-transition states, we probed the equilibrium and kinetic consequences of replacing core valines with isosteric threonines at six positions in azurin. In contrast to regular hydrophobic-to-alanine phi-value analysis, valine-to-threonine mutations do not disrupt the core packing but stabilize the unfolded state and can be used to assess the degree of solvation in the folding-transition state upon combination with regular phi-values. We find that the transition state for folding of apo-azurin appears completely dry, while that for zinc-azurin involves partially formed interactions that engage water molecules. This distinct difference between the apo-and holo-folding nuclei can be rationalized in terms of the shape of the free-energy barrier.

The Physics of the Interactions Governing Folding and Association of Proteins

The review discusses the molecular origins of the forces and free energies that determine several things about proteins, and how experiment and theory reveal this information. The first subject is the stability of the folded, native structures. The second is the range of molecular mechanisms by which proteins find their way to those folded structures in laboratory environments. The third is the much more complex problem of how folding occurs in the cellular environment. This topic includes a discussion of crowding and of the roles of chaperone molecules. The review concludes with a discussion of protein aggregation and fibril formation and of misfolding and therapies associated with it.

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.

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.

Formation of the folding nucleus of an SH3 domain investigated by loosely coupled molecular dynamics simulations

The experimentally well-established folding mechanism of the src-SH3 domain, and in particular the phi-value analysis of its transition state, represents a sort of testing table for computational investigations of protein folding. Here, parallel molecular dynamics simulations of the src-SH3 domain have been performed starting from denatured conformations. By rescuing and restarting only trajectories approaching the folding transition state, an ensemble of conformations was obtained with a completely structured central beta-sheet and a native-like packing of residues Ile-110, Ala-121, and Ile-132. An analysis of the trajectories shows that there are several pathways leading to the formation of the central beta-sheet whereas its two hairpins form in a different but consistent way.

Sparsely populated folding intermediates of the Fyn SH3 domain: matching native-centric essential dynamics and experiment

A complete description of how a protein folds requires the characterization of intermediate conformations traversed during the folding transition. We have calculated dynamics trajectories of a simplified model of the Fyn SH3 domain with a native-centric potential energy function. Analysis of the resulting site-resolved energy trajectory identifies an ensemble of intermediate conformations for folding and another for unfolding. The model's folding intermediate is structured in the three beta-strands that make up the protein's core and is strikingly similar to intermediates detected in a recent NMR study of Fyn SH3 folding and to folding transition states elucidated in mutagenesis studies of SH3 domains. The unfolding intermediate is formed by dissociation of the folded protein's two terminal beta-strands from its core. The presence of such an intermediate is consistent with the results of a protein-engineering study on the src SH3 domain showing that these strands separate before the rate-limiting step of unfolding. Despite the presence of these conformations intermediate between the native and fully unfolded states, the computed heat capacity vs. temperature profile of the model protein indicates that its thermodynamics satisfies the usual calorimetric criterion for two-state folding. This observation highlights the fact that, if not properly interpreted, methods such as calorimetry that do not probe multiple sites in a molecule can lead to an oversimplified view of folding. The close agreement between results from this simplified model and experimental work underscores the important contributions that computational methods can make in providing insights into protein folding.

Temperature dependence of the free energy landscape of the src-SH3 protein domain

The temperature dependence of the free energy landscape of the src-SH3 protein domain is investigated through fully atomic simulations in explicit solvent. Simulations are performed above and below the folding transition temperature, enabling an analysis of both protein folding and unfolding. The transition state for folding and unfolding, identified from the free energy surfaces, is found to be very similar, with structure in the central hydrophobic sheet and little structure throughout the rest of the protein. This is a result of a polarized folding (unfolding) mechanism involving early formation (late loss) of the central hydrophobic sheet at the transition state. Unfolding simulations map qualitatively well onto low-temperature free energy surfaces but appear, however, to miss important features observed in folding simulations. In particular, details of the folding mechanism involving the opening and closing of the hydrophobic core are not captured by unfolding simulations performed under strongly denaturing conditions. In addition, free energy surfaces at high temperatures do not display a desolvation barrier found at lower temperatures, involving the expulsion of water molecules from the hydrophobic core.

Transition states for protein folding have native topologies despite high structural variability

We present a structural analysis of the folding transition states of three SH3 domains. Our results reveal that the secondary structure is not yet fully formed at this stage of folding and that the solvent is only partially excluded from the interior of the protein. Comparison of the members of the transition state ensemble with a database of native folds shows that, despite substantial local variability, the transition state structures can all be classified as having the topology characteristic of an SH3 domain. Our results suggest a mechanism for folding in which the formation of a network of interactions among a subset of hydrophobic residues ensures that the native topology is generated. Such a mechanism enables high fidelity in folding while minimizing the need to establish a large number of specific interactions in the conformational search.