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.     

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.     

Kinetic consequences of native state optimization of surface-exposed electrostatic interactions in the Fyn SH3 domain

Optimization of surface exposed charge-charge interactions in the native state has emerged as an effective means to enhance protein stability; but the effect of electrostatic interactions on the kinetics of protein folding is not well understood. To investigate the kinetic consequences of surface charge optimization, we characterized the folding kinetics of a Fyn SH3 domain variant containing five amino acid substitutions that was computationally designed to optimize surface charge-charge interactions. Our results demonstrate that this optimized Fyn SH3 domain is stabilized primarily through an eight-fold acceleration in the folding rate. Analyses of the constituent single amino acid substitutions indicate that the effects of optimization of charge-charge interactions on folding rate are additive. This is in contrast to the trend seen in folded state stability, and suggests that electrostatic interactions are less specific in the transition state compared to the folded state. Simulations of the transition state using a coarse-grained chain model show that native electrostatic contacts are weakly formed, thereby making the transition state conducive to nonspecific, or even nonnative, electrostatic interactions. Because folding from the unfolded state to the folding transition state for small proteins is accompanied by an increase in charge density, nonspecific electrostatic interactions, that is, generic charge density effects can have a significant contribution to the kinetics of protein folding. Thus, the interpretation of the effects of amino acid substitutions at surface charged positions may be complicated and consideration of only native-state interactions may fail to provide an adequate picture.     

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.     

Molecular dynamics simulations from putative transition states of alpha-spectrin SH3 domain

A series of molecular dynamics simulations in explicit solvent were started from nine structural models of the transition state of the SH3 domain of alpha-spectrin, which were generated by Lindorff-Larsen et al. (Nat Struct Mol Biol 2004;11:443-449) using molecular dynamics simulations in which experimental Phi - values were incorporated as restraints. Two of the nine models were simulated 10 times for 200 ns and the remaining models simulated two times for 200 ns. Complete folding was observed in one case, while in the other simulations partial folding or unfolding events were observed, which were characterized by a regularization of elements of secondary structure. These results are consistent with recent experimental evidence that the folding of SH3 domains involves low populated intermediate states.     

Protein folding kinetics beyond the phi value: using multiple amino acid substitutions to investigate the structure of the SH3 domain folding transition state

The SH3 domain folding transition state structure contains two well-ordered turn regions, known as the diverging turn and the distal loop. In the Src SH3 domain transition state, these regions are stabilized by a hydrogen bond between Glu30 in the diverging turn and Ser47 in the distal loop. We have examined the effects on folding kinetics of amino acid substitutions at the homologous positions (Glu24 and Ser41) in the Fyn SH3 domain. In contrast to most other folding kinetics studies which have focused primarily on non-disruptive substitutions with Ala or Gly, here we have examined the effects of substitutions with diverse amino acid residues. Using this approach, we demonstrate that the transition state structure is generally tolerant to amino acid substitutions. We also uncover a unique role for Ser at position 41 in facilitating folding of the distal loop, which can only be replicated by Asp at the same position. Both these residues appear to accelerate folding through the formation of short-range side-chain to backbone hydrogen bonds. The folding of the diverging turn region is shown to be driven primarily by local interactions. The diverging turn and distal loop regions are found to interact in the transition state structure, but only in the context of particular mutant backgrounds. This work demonstrates that studying the effects of a variety of amino acid substitutions on protein folding kinetics can provide unique insights into folding mechanisms which cannot be obtained by standard Phi value analysis.     

Probing the energy landscape of protein folding/unfolding transition states

Previous molecular dynamics (MD) simulations of the thermal denaturation of chymotrypsin inhibitor 2 (CI2) have provided atomic-resolution models of the transition state ensemble that is well supported by experimental studies. Here, we use simulations to further investigate the energy landscape around the transition state region. Nine structures within approximately 35 ps and 3 A C(alpha) RMSD of the transition state ensemble identified in a previous 498 K thermal denaturation simulation were quenched under the quasi-native conditions of 335 K and neutral pH. All of the structures underwent hydrophobically driven collapse in response to the drop in temperature. Structures less denatured than the transition state became structurally more native-like, while structures that were more denatured than the transition state tended to show additional loss of native structure. The structures in the immediate region of the transition state fluctuated between becoming more and less native-like. All of the starting structures had the same native-like topology and were quite similar (within 3.5 A C(alpha) RMSD). That the structures all shared native-like topology, yet diverged into either more or less native-like structures depending on which side of the transition state they occupied on the unfolding trajectory, indicates that topology alone does not dictate protein folding. Instead, our results suggest that a detailed interplay of packing interactions and interactions with water determine whether a partially denatured protein will become more native-like under refolding conditions.     

A thermodynamic and kinetic analysis of the folding pathway of an SH3 domain entropically stabilised by a redesigned hydrophobic core

The folding thermodynamics and kinetics of the alpha-spectrin SH3 domain with a redesigned hydrophobic core have been studied. The introduction of five replacements, A11V, V23L, M25V, V44I and V58L, resulted in an increase of 16% in the overall volume of the side-chains forming the hydrophobic core but caused no remarkable changes to the positions of the backbone atoms. Judging by the scanning calorimetry data, the increased stability of the folded structure of the new SH3-variant is caused by entropic factors, since the changes in heat capacity and enthalpy upon the unfolding of the wild-type and mutant proteins were identical at 298 K. It appears that the design process resulted in an increase in burying both the hydrophobic and hydrophilic surfaces, which resulted in a compensatory effect upon the changes in heat capacity and enthalpy. Kinetic analysis shows that both the folding and unfolding rate constants are higher for the new variant, suggesting that its transition state becomes more stable compared to the folded and unfolded states. The phi(double dagger-U) values found for a number of side-chains are slightly lower than those of the wild-type protein, indicating that although the transition state ensemble (TSE) did not change overall, it has moved towards a more denatured conformation, in accordance with Hammond's postulate. Thus, the acceleration of the folding-unfolding reactions is caused mainly by an improvement in the specific and/or non-specific hydrophobic interactions within the TSE rather than by changes in the contact order. Experimental evidence showing that the TSE changes globally according to its hydrophobic content suggests that hydrophobicity may modulate the kinetic behaviour and also the folding pathway of a protein.     

Cooperative interactions and a non-native buried Trp in the unfolded state of an SH3 domain

The presence of residual structure in the unfolded state of the N-terminal SH3 domain of Drosophila drk (drkN SH3 domain) has been investigated using far- and near-UV circular dichroism (CD), fluorescence, and NMR spectroscopy. The unfolded (U(exch)) state of the drkN SH3 domain is significantly populated and exists in equilibrium with the folded (F(exch)) state under non-denaturing conditions near physiological pH. Denaturation experiments have been performed on the drkN SH3 domain in order to monitor the change in ellipticity, fluorescence intensity, and chemical shift between the U(exch) state and chemically or thermally denatured states. Differences between the unfolded and chemically or thermally denatured states highlight specific areas of residual structure in the unfolded state that are cooperatively disrupted upon denaturation. Results provide evidence for cooperative interactions in the unfolded state involving residues of the central beta-sheet, particularly the beta4 strand. Denaturation as well as hydrogen-exchange experiments demonstrate a non-native burial of the Trp ring within this "cooperative" core of the unfolded state. These findings support the presence of non-native hydrophobic clusters, organised by Trp rings, within disordered states.     

The osmolyte trimethylamine‐N‐oxide stabilizes the Fyn SH3 domain without altering the structure of its folding transition state

Trimethylamine‐N‐oxide (TMAO) is a naturally occurring osmolyte that stabilizes proteins against denaturation. Although the impact of TMAO on the folding thermodynamics of many proteins has been well characterized, far fewer studies have investigated its effects on protein folding kinetics. In particular, no previous studies have used Φ‐value analysis to determine whether TMAO may alter the structure of the folding transition state. Here we have measured the effects on folding kinetics of 16 different amino acid substitutions distributed across the structure of the Fyn SH3 domain both in the presence and absence of TMAO. The folding and unfolding rates in TMAO, on average, improved to equivalent degrees, with a twofold increase in the protein folding rate accompanied by a twofold decrease in the unfolding rate. Importantly, TMAO caused little alteration to the Φ‐values of the mutants tested, implying that this compound minimally perturbs the folding transition state structure. Furthermore, the solvent accessibility of the transition state was not altered as reflected in an absence of a TMAO‐induced change in the denaturant β factors. Through TMAO‐induced folding studies, a β factor of 0.5 was calculated for this compound, suggesting that the protein backbone, which is the target of action of TMAO, is 50% exposed in the transition state as compared to the native state. This finding is consistent with the equivalent effects of TMAO on the folding and unfolding rates. Through thermodynamic analysis of mutants, we also discovered that the stabilizing effect of TMAO is lessened with increasing temperature.     

The identification of conserved interactions within the SH3 domain by alignment of sequences and structures

The SH3 domain, comprised of approximately 60 residues, is found within a wide variety of proteins, and is a mediator of protein-protein interactions. Due to the large number of SH3 domain sequences and structures in the databases, this domain provides one of the best available systems for the examination of sequence and structural conservation within a protein family. In this study, a large and diverse alignment of SH3 domain sequences was constructed, and the pattern of conservation within this alignment was compared to conserved structural features, as deduced from analysis of eighteen different SH3 domain structures. Seventeen SH3 domain structures solved in the presence of bound peptide were also examined to identify positions that are consistently most important in mediating the peptide-binding function of this domain. Although residues at the two most conserved positions in the alignment are directly involved in peptide binding, residues at most other conserved positions play structural roles, such as stabilizing turns or comprising the hydrophobic core. Surprisingly, several highly conserved side-chain to main-chain hydrogen bonds were observed in the functionally crucial RT-Src loop between residues with little direct involvement in peptide binding. These hydrogen bonds may be important for maintaining this region in the precise conformation necessary for specific peptide recognition. In addition, a previously unrecognized yet highly conserved beta-bulge was identified in the second beta-strand of the domain, which appears to provide a necessary kink in this strand, allowing it to hydrogen bond to both sheets comprising the fold.     

Native state energetics of the Src SH2 domain: evidence for a partially structured state in the denatured ensemble

We have defined the free-energy profile of the Src SH2 domain using a variety of biophysical techniques. Equilibrium and kinetic experiments monitored by tryptophan fluorescence show that Src SH2 is quite stable and folds rapidly by a two-state mechanism, without populating any intermediates. Native state hydrogen-deuterium exchange confirms this two-state behavior; we detect no cooperative partially unfolded forms in equilibrium with the native conformation under any conditions. Interestingly, the apparent stability of the protein from hydrogen exchange is 2 kcal/mol greater than the stability determined by both equilibrium and kinetic studies followed by fluorescence. Native-state proteolysis demonstrates that this "super protection" does not result from a deviation from the linear extrapolation model used to fit the fluorescence data. Instead, it likely arises from a notable compaction in the unfolded state under native conditions, resulting in an ensemble of conformations with substantial solvent exposure of side chains and flexible regions sensitive to proteolysis, but backbone amides that exchange with solvent approximately 30-fold slower than would be expected for a random coil. The apparently simple behavior of Src SH2 in traditional unfolding studies masks the significant complexity present in the denatured-state ensemble.     

The design of a hyperstable mutant of the Abp1p SH3 domain by sequence alignment analysis

We have characterized the thermodynamic stability of the SH3 domain from the Saccharomyces cerevisiae Abp1p protein and found it to be relatively low compared to most other SH3 domains, with a Tm of 60 degrees C and a deltaGu of 3.08 kcal/mol. Analysis of a large alignment of SH3 domains led to the identification of atypical residues at eight positions in the wild-type Abp1p SH3 domain sequence that were subsequently replaced by the residue seen most frequently at that position in the alignment. Three of the eight mutants constructed in this way displayed increases in Tm ranging from 8 to 15 degrees C with concomitant increases in deltaGu of up to 1.4 kcal/mol. The effects of these substitutions on folding thermodynamics and kinetics were entirely additive, and a mutant containing all three was dramatically stabilized with a Tm greater than 90 degrees C and a deltaGu more than double that of the wild-type domain. The folding rate of this hyperstable mutant was 10-fold faster than wild-type, while its unfolding rate was fivefold slower. All of the stabilized mutants were still able to bind a target peptide with wild-type affinity. We have analyzed the stabilizing amino acid substitutions isolated in this study and several other similar sequence alignment based studies. In approximately 25% of cases, increased stability can be explained by enhanced propensity of the substituted residue for the local backbone conformation at the mutagenized site.     

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.     

Design and structural thermodynamic studies of the chimeric protein derived from spectrin SH3 domain

A number of the chimeric constructs with spectrin SH3 domain were designed for structural and thermodynamic studies of protein self-assembly and protein-ligand interactions. SH3 domains, components of many regulatory proteins, operate through weak interactions with proline-rich regions of polypeptide chains. The recombinant construct (WT-CIIA) studied in this work was constructed by linking the peptide ligand PPPVPPYSAG to the domain C-terminus via a long 12-residue linker to increase the affinity of this ligand for the spectrin domain, thereby ensuring a stable positioning of the polyproline helix to the conserved ligand-binding site in orientation II, which is regarded as untypical of the interaction between this domain and oligopeptides. A comparison of fluorescence spectra of the initial domain and the recombinant protein WT-CIIA suggests that the ligand sticks to the conservative binding site. However, analysis of the equilibrium urea-induced unfolding has demonstrated that this is an unstable contact, which leads to a two-stage unfolding of the chimeric protein. The protein WT-CIIA was crystallized; a set of X-ray diffraction data with a resolution of 1.75 Å was recorded from individual crystals. A preliminary analysis of these diffraction data has demonstrated that the crystals belong to space group P32 with the following unit cell parameters: a = b = 36.39, c = 112.17 Å, a = β = 90.0, and γ = 120.0.     

Thermodynamics and kinetics of non-native interactions in protein folding: a single point mutant significantly stabilizes the N-terminal domain of L9 by modulating non-native interactions in the denatured state

Comparatively little is known about the role of non-native interactions in protein folding and their role in both folding and stability is controversial. We demonstrate that non-native electrostatic interactions involving specific residues in the denatured state can have a significant effect upon protein stability and can persist in the transition state for folding. Mutation of a single surface exposed residue, Lys12 to Met, in the N-terminal domain of the ribosomal protein L9 (NTL9), significantly increased the stability of the protein and led to faster folding. Structural and energetic studies of the wild-type and K12M mutant show that the 1.9 kcal mol(-1) increase in stability is not due to native state effects, but rather is caused by modulation of specific non-native electrostatic interactions in the denatured state. pH dependent stability measurements confirm that the increased stability of the K12M is due to the elimination of favorable non-native interactions in the denatured state. Kinetic studies show that the non-native electrostatic interactions involving K12 persist in the transition state. The analysis demonstrates that canonical Phi-values can arise from the disruption of non-native interactions as well as from the development of native interactions.     

Artificial proteins as allosteric modulators of PDZ3 and SH3 in two‐domain constructs: A computational characterization of novel chimeric proteins

Artificial multidomain proteins with enhanced structural and functional properties can be utilized in a broad spectrum of applications. The design of chimeric fusion proteins utilizing protein domains or one‐domain miniproteins as building blocks is an important advancement for the creation of new biomolecules for biotechnology and medical applications. However, computational studies to describe in detail the dynamics and geometry properties of two‐domain constructs made from structurally and functionally different proteins are lacking. Here, we tested an in silico design strategy using all‐atom explicit solvent molecular dynamics simulations. The well‐characterized PDZ3 and SH3 domains of human zonula occludens (ZO‐1) (3TSZ), along with 5 artificial domains and 2 types of molecular linkers, were selected to construct chimeric two‐domain molecules. The influence of the artificial domains on the structure and dynamics of the PDZ3 and SH3 domains was determined using a range of analyses. We conclude that the artificial domains can function as allosteric modulators of the PDZ3 and SH3 domains. Proteins 2016; 84:1358–1374. © 2016 Wiley Periodicals, Inc.     

Hydrogen exchange and ligand binding: ligand-dependent and ligand-independent protection in the Src SH3 domain

Amide hydrogen-deuterium exchange has proven to be a powerful tool for detecting and characterizing high-energy conformations in protein ensembles. Since interactions with ligands can modulate these high-energy conformations, hydrogen exchange appears to be an ideal experimental probe of the physical mechanisms underlying processes like allosteric regulation. The chemical mechanism of hydrogen exchange, however, can complicate such studies. Here, we examine hydrogen exchange rates in a simple model system, the c-Src SH3 domain interacting with a short peptide ligand. Addition of ligand slows the rates of hydrogen exchange at nearly every amide for which we can obtain data. Careful analysis, however, reveals that this slowing is due primarily to a reduction in the population of free protein in the system, and not to any specific property of the complex. We present a method to separate the contributions of free and bound protein to the exchange kinetics that has allowed us to identify the subset of amides where exchange arises directly from the complex. These results demonstrate that the slowing of hydrogen exchange induced by ligand interactions should be interpreted with caution, and more extensive experiments are required to correlate changes in hydrogen exchange with changes in structure or internal dynamics.     

Structural and kinetic characterization of the simplified SH3 domain FP1

The simplified SH3 domain sequence, FP1, obtained in phage display selection experiments has an amino acid composition that is 95% Ile, Lys, Glu, Ala, Gly. Here we use NMR to investigate the tertiary structure of FP1. We find that the overall topology of FP1 resembles that of the src SH3 domain, the hydrogen-deuterium exchange and chemical shift perturbation profiles are similar to those of naturally occurring SH3 domains, and the (15)N relaxation rates are in the range of naturally occurring small proteins. Guided by the structure, we further simplify the FP1 sequence and compare the effects on folding kinetics of point mutations in FP1 and the wild-type src SH3 domain. The results suggest that the folding transition state of FP1 is similar to but somewhat less polarized than that of the wild-type src SH3 domain.