pH dependence of the hydrogen exchange in the SH3 domain of alpha-spectrin

Using nuclear magnetic resonance we have measured the hydrogen exchange (HX) in the Src homology region 3 (SH3) domain of alpha-spectrin as a function of pH*. At very acidic pH* values the exchange of most residues appears to occur via global unfolding, although several residues show abnormally large Gibbs energies of exchange, suggesting the presence of some residual structure in the unfolded state. At higher pH* HX occurs mainly via local or partial unfoldings. We have been able to characterize the coupling between the electrostatic interactions in this domain and the conformational fluctuations occurring under native conditions by analyzing the dependence upon pH* of the Gibbs energy of exchange. The SH3 domain seems to be composed of a central core, which requires large structural disruptions to become exposed to the solvent, surrounded by smaller subdomains, which fluctuate independently.     

Identification of rare partially unfolded states in equilibrium with the native conformation in an all beta-barrel protein

Human acidic fibroblast growth factor 1 (hFGF-1) is an all beta-barrel protein, and the secondary structural elements in the protein include 12 antiparallel beta-strands arranged into a beta-trefoil fold. In the present study, we investigate the stability of hFGF-1 by hydrogen-deuterium exchange as a function of urea concentration. Urea-induced equilibrium unfolding of hFGF-1 monitored by fluorescence and CD spectroscopy suggests that the protein unfolds by a two-state (native to denatured) mechanism. Hydrogen exchange in hFGF-1, under the experimental conditions used, occurs by the EX2 mechanism. In contrast to the equilibrium unfolding events monitored by optical probes, native state hydrogen exchange data show that the beta-trefoil architecture of hFGF-1 does not behave as a single cooperative unit. There are at least two structurally independent units with differing stabilities in hFGF-1. Beta-strands I, II, III, VI, VII, X, XI, and XII fit into the global unfolding isotherm. By contrast, residues in beta-strands IV, V, VIII, and IX exchange by the subfolding isotherm and could be responsible for the occurrence of high-energy partially unfolded state(s) in hFGF-1. There appears to be a broad continuum of stabilities among the four beta-strands (beta-strands IV, V, VIII, and IX) constituting the subglobal folding unit. The slow exchanging residues in hFGF-1 do not represent the folding nucleus of the protein.     

Predicting the equilibrium protein folding pathway: Structure-based analysis of staphylococcal nuclease

The equilibrium folding pathway of staphylococcal nucleas (SNase) has been approximated using a statistical thermodynamic formalism that utilizes the high-resolution structure of the native state as a template to generate a large ensemble of partially folded states. Close to 400,000 different states ranging from the native to the completely unfolded states were included in the analysis. The probability of each state was estimated using an empirical structural parametrization of the folding energetics. It is shown that this formalism predicts accurately the stability of the protein, the cooperativity of the folding/unfolding transition observed by differential scanning calorimetry (DSC) or urea denaturation and the thermodynamic parameters for unfolding. More importantly, this formalism provides a quantitative account of the experimental hydrogen exchange protection factors measured under native conditions for SNase. These results suggest that the computer-generated distribution of states approximates well the ensemble of conformations existing in solution. Furthermore, this formalism represents the first model capable of quantitatively predicting within a unified framework the probability distribution of states seen under native conditions and its change upon unfolding. © 1997 Wiley-Liss, Inc.     

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.     

Detection of a hidden folding intermediate of the third domain of PDZ

The folding pathway of the third domain of PDZ from the synaptic protein PSD-95 was characterized using kinetic and equilibrium methods by monitoring the fluorescence signal from a Trp residue that is incorporated at a near-surface position. Kinetic folding of this domain showed multiple exponential phases, whereas unfolding showed a single exponential phase. The slow kinetic phases were attributed to isomerization of proline residues, since there are five proline residues in this domain. We found that the logarithms of the rate constants for the fast phase of folding and unfolding are linearly dependent on the concentrations of denaturant. The unfolding free energy derived from these rate constants at zero denaturant was close to the value measured using the equilibrium method, suggesting the absence of detectable sub-millisecond folding intermediates. However, native-state hydrogen exchange experiments detected a partially unfolded intermediate under native conditions. It was further confirmed by a protein engineering study. These data suggest that a hidden intermediate exists after the rate-limiting step in the folding of the third domain of PDZ.     

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

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

Hierarchy of structure loss in MD simulations of src SH3 domain unfolding.

To complement experimental studies of the src SH3 domain folding, we studied 30 independent, high-temperature, molecular dynamics simulations of src SH3 domain unfolding. These trajectories were observed to differ widely from each other. Thus, rather than analyzing individual trajectories, we sought to identify the recurrent features of the high-temperature unfolding process. The conformations from all simulations were combined and then divided into groups based on the number of native contacts. Average occupancies of each side-chain hydrophobic contact and hydrogen bond in the protein were then determined. In the symmetric funnel limit, the occupancies of all contacts should decrease in concert with the loss in total number of native contacts. If there is a lack of symmetry or hierarchy to the unfolding process, the occupancies of some contacts should decrease more slowly, and others more rapidly. Despite the heterogeneity of the individual trajectories, the ensemble averaging revealed an order to the unfolding process: contacts between the N and C-terminal strands are the first to disappear, whereas contacts within the distal beta-hairpin and a hydrogen-bonding network involving the distal loop beta-turn and the diverging turn persist well after the majority of the native contacts are lost. This hierarchy of events resembles but is somewhat less pronounced than that observed in our experimental studies of the folding of src SH3 domain.     

Hydrogen exchange of the tetramerization domain of the human tumour suppressor p53 probed by denaturants and temperature.

We have analysed the hydrogen/deuterium exchange of the tetramerization domain of human tumour suppressor p53 under mild chemical denaturation conditions, and at different temperatures. Exchange behaviour has been measured for 16 amide protons in the chemical-denaturation studies and for seven protons in the temperature-denaturation studies. The exchange rates are in the range observed for other proteins with similar elements of secondary structure. The slowest-exchange core includes contributions from residues in the alpha helix and the beta sheet. However, only some of the slowest-exchanging protons correspond to residues involved in native interactions in the transient intermediate detected during the folding of this domain. The guanidinium-chloride denaturation curves of all residues seem to merge together, although they are well below the main isotherm of global unfolding. Thus, there is no evidence for several subglobal unfolding units. The activation parameters obtained from the temperature-denaturation experiments are similar to those obtained for monomeric proteins, and well below the global unfolding enthalpy obtained by circular dichroism measurements. Thus, the exchange studies at different denaturant concentrations and temperatures indicate that no particular folding intermediate is populated under those conditions.     

Stable intermediate states and high energy barriers in the unfolding of GFP

We present a study of the denaturation of a truncated, cycle3 variant of green fluorescent protein (GFP). Chemical denaturation is used to unfold the protein, with changes in structure being monitored by the green fluorescence, tyrosine fluorescence and far-UV circular dichroism. The results show that the denaturation behaviour of GFP is complex compared to many small proteins: equilibrium is established only very slowly, over the time course of weeks, suggesting that there are high folding/unfolding energy barriers. Unfolding kinetics confirm that the rates of unfolding at low concentrations of denaturant are very low, consistent with the slow establishment of the equilibrium. In addition, we find that GFP significantly populates an intermediate state under equilibrium conditions, which is compact and stable with respect to the unfolded state (m(IU)=4.6 kcal mol(-1) M(-1) and Delta G(IU)=12.5 kcal mol(-1)). The global and local stability of GFP was probed further by measuring the hydrogen/deuterium (H/D) NMR exchange rates of more than 157 assigned amide protons. Analysis at two different values of pH showed that amide protons within the beta-barrel structure exchange at the EX2 limit, consequently, free energies of exchange could be calculated and compared to those obtained from the denaturation-curve studies providing further support for the three-state model and the existence of a stable intermediate state. Analysis reveals that amide protons in beta-strands 7, 8, 9 and 10 have, on average, higher exchange rates than others in the beta-barrel, suggesting that there is greater flexibility in this region of the protein. Forty or so amide protons were found which do not undergo significant exchange even after several months and these are clustered into a core region encompassing most of the beta-strands, at least at one end of the barrel structure. It is likely that these residues play an important role in stabilizing the structure of the intermediate state. The intermediate state observed in the chemical denaturation studies described here, is similar to that observed at pH 4 in other studies.     

Mapping the interactions present in the transition state for unfolding/folding of FKBP12.

The structure of the transition state for folding/unfolding of the immunophilin FKBP12 has been characterised using a combination of protein engineering techniques, unfolding kinetics, and molecular dynamics simulations. A total of 34 mutations were made at sites throughout the protein to probe the extent of secondary and tertiary structure in the transition state. The transition state for folding is compact compared with the unfolded state, with an approximately 30 % increase in the native solvent-accessible surface area. All of the interactions are substantially weaker in the transition state, as probed by both experiment and molecular dynamics simulations. In contrast to some other proteins of this size, no element of structure is fully formed in the transition state; instead, the transition state is similar to that found for smaller, single-domain proteins, such as chymotrypsin inhibitor 2 and the SH3 domain from alpha-spectrin. For FKBP12, the central three strands of the beta-sheet, beta-strand 2, beta-strand 4 and beta-strand 5, comprise the most structured region of the transition state. In particular Val101, which is one of the most highly buried residues and located in the middle of the central beta-strand, makes approximately 60 % of its native interactions. The outer beta-strands and the ends of the central beta-strands are formed to a lesser degree. The short alpha-helix is largely unstructured in the transition state, as are the loops. The data are consistent with a nucleation-condensation model of folding, the nucleus of which is formed by side-chains within beta-strands 2, 4 and 5, and the C terminus of the alpha-helix. The precise residues involved in the nucleus differ in the two simulated transition state ensembles, but the interacting regions of the protein are conserved. These residues are distant in the primary sequence, demonstrating the importance of tertiary interactions in the transition state. The two independently derived transition state ensembles are structurally similar, which is consistent with a Bronsted analysis confirming that the transition state is an ensemble of states close in structure.     

Cooperative propagation of local stability changes from low-stability and high-stability regions in a SH3 domain

Site-directed mutagenesis has been used to produce local stability changes at two regions of the binding site surface of the alpha-spectrin SH3 domain (Spc-SH3) differing in their intrinsic stability. Mutations were made at residue 56, located at the solvent-exposed side of the short 3(10) helix, and at residue 21 in the tip of the flexible RT-loop. NMR chemical-shift analysis and X-ray crystallography indicated negligible changes produced by the mutations in the native structure limited to subtle rearrangements near the mutated residue and at flexible loops. Additionally, mutations do not alter importantly the SH3 binding site structure, although produce significant changes in its affinity for a proline-rich decapeptide. The changes in global stability measured by differential scanning calorimetry are consistent the local energy changes predicted by theoretical models, with the most significant effects observed for the Ala-Gly mutations. Propagation of the local stability changes throughout the domain structure has been studied at a per-residue level of resolution by NMR-detected amide hydrogen-deuterium exchange (HX). Stability propagation is remarkably efficient in this small domain, apparently due to its intrinsically low stability. Nevertheless, the HX-core of the domain is not fully cooperative, indicating the existence of co-operative subunits within the core, which is markedly polarized. An equilibrium phi-analysis of the changes in the apparent Gibbs energies of HX per residue produced by the mutations has allowed us to characterize structurally the conformational states leading to HX. Some of these states resemble notably the folding transition state of the Spc-SH3 domain, suggesting a great potential of this approach to explore the folding energy landscape of proteins. An energy perturbation propagates more effectively from a flexible region to the core than in the opposite direction, because the former affects a broader region of the energy landscape than the latter. This might be of importance in understanding the special thermodynamic signature of the SH3-peptide interaction and the relevance of the dual character of SH3 binding sites.     

Altered dynamics in Lck SH3 upon binding to the LBD1 domain of Herpesvirus saimiri Tip

The Tip protein from Herpesvirus saimiri interacts with the SH3 domain from the Src-family kinase Lck via a proline-containing sequence termed LBD1. Src-family kinase SH3 domains related to Lck have been shown to be dynamic in solution and partially unfold under physiological conditions. The rate of such partial unfolding is reduced by viral protein binding. To determine if the Lck SH3 domain displayed similar behavior, the domain was investigated with hydrogen exchange and mass spectrometry. Lck SH3 was found to be highly dynamic in solution. While other SH3 domains require as much as 10,000 sec to become totally deuterated, Lck SH3 became almost completely labeled within 200 sec. A partial unfolding event involving 8-10 residues was observed with a half-life of approximately 10 sec. Tip LBD1 binding did not cause gross structural changes in Lck SH3 but globally stabilized the domain and reduced the rate of partial unfolding by a factor of five. The region of partial unfolding in Lck SH3 was found to be similar to that identified for other SH3 domains that partially unfold. Although the sequence conservation between Lck SH3 and other closely related SH3 domains is high, the dynamics do not appear to be conserved.     

Secondary structures and structural fluctuation in a dimeric protein, Streptomyces subtilisin inhibitor.

Based on the nuclear magnetic resonance assignments of a dimeric protein, Streptomyces subtilisin inhibitor (SSI), microscopic details of secondary structures in solution have been elucidated. The chemical shift index of C(alpha) signals, together with information on the hydrogen exchange rates of the backbone amide protons, were used to identify secondary structures. The locations of these secondary structures were found to be different in some critical points from those determined earlier by X-ray crystallography of the crystal. Notably, the beta3 strand is completely missing and the alpha2 helix is extended toward the C-terminus. Furthermore, hydrogen exchange experiments of individual peptide NH protons under strongly folding conditions revealed mechanisms of global and local structural fluctuation within the dimeric structure. It has been suggested that the global fluctuation of the monomeric unit occurs without affecting the accompanying monomer, in contrast to the equilibrium thermal unfolding, which is cooperative. Higher protection against hydrogen exchange for residues in part of the beta4 strand implies that this region might serve as a folding core.     

Partially unfolded forms and non-two-state folding of a beta-sandwich: FHA domain from Arabidopsis receptor kinase-associated protein phosphatase

FHA domains adopt a beta-sandwich fold with 11 strands. The first evidence of partially unfolded forms of a beta-sandwich is derived from native-state hydrogen exchange (NHX) of the forkhead-associated (FHA) domain from kinase-associated protein phosphatase from Arabidopsis. The folding kinetics of this FHA domain indicate that EX2 behavior prevails at pH 6.3. In the chevron plot, rollover in the folding arm and bends in the unfolding arm suggest folding intermediates. NHX of this FHA domain suggests a core of six most stable beta-strands and two loops, characterized by rare global unfolding events. Flanking this stable core are beta-strands and recognition loops with less stability, termed subglobal motifs. These suggest partially unfolded forms (near-native intermediates) with two levels of stability. The spatial separation of the subglobal motifs on the flanks suggests possible parallelism in their folding as additional beta-strands align with the stable core of six strands. Intermediates may contribute to differences in stabilities and m-values suggested by NHX or kinetics relative to chemical denaturation. Residual structure in the unfolded regime is suggested by superprotection of beta-strand 6 and by GdmCl-dependence of adjustments in amide NMR spectra and residual optical signal. The global folding stability depends strongly on pH, with at least 3 kcal/mol more stability at pH 7.3 than at pH 6.3. This FHA domain is hypothesized to fold progressively with initial hydrophobic collapse of its stable six-stranded core followed by addition of less stable flanking beta-strands and ordering of recognition loops.     

A thermodynamic analysis of a family of small globular proteins: SH3 domains.

The stability and folding thermodynamics of two SH3-domains, belonging to Fyn and Abl proteins, have been studied by scanning calorimetry and urea-induced unfolding. They undergo an essentially two-state unfolding with parameters similar to those of the previously studied alpha-spectrin SH3 domain. The correlations between the thermodynamic parameters (heat capacity increment, delta Cp,U, the proportionality factor, m, and the Gibbs energy, delta Gw298) of unfolding and some integral structural parameters, such as polar and non-polar areas exposed upon domain denaturation, have been analyzed. The experimental data on delta Cp,U and the m-factor of the linear extrapolation model (LEM) obey the simple empirical correlations deduced elsewhere. The Gibbs energies calculated from the DSC data were compared with those found by fitting urea-unfolding curves to the LEM and the denaturant-binding model (DBM). The delta Gw298 values found with DBM correlate better with the DSC data, while those obtained with LEM are systematically smaller. The systematic difference between the parameters calculated with LEM and DBM are explained by an inherent imperfection of the LEM.     

A single mutation induces amyloid aggregation in the alpha-spectrin SH3 domain: analysis of the early stages of fibril formation

The Src-homology region 3 domain of chicken alpha-spectrin (Spc-SH3) is a small two-state folding protein, which has never been described to form amyloid fibrils under any condition investigated so far. We show here that the mutation of asparagine 47 to alanine at the distal loop, which destabilises similarly the native and folding transition states of the domain, induces the formation of amyloid fibrils under mild acid conditions. Amyloid aggregation of the mutant is enhanced by the increase in temperature, protein concentration and NaCl concentration. The early stages of amyloid formation have been monitored as a function of time and temperature using a variety of biophysical methods. Differential scanning calorimetry experiments under conditions of amyloid formation have allowed the identification of different thermal transitions corresponding to conformational and aggregation processes as well as to the high-temperature disaggregation and unfolding of the amyloid fibrils. Aggregation is preceded by a rapid conformational change in the monomeric domain involving about 40% of the global unfolding enthalpy, considerable change in secondary structure, large loss of tertiary structure and exposure of hydrophobic patches to the solvent. The conformational change is followed by formation of a majority of oligomeric species with apparent hydrodynamic radius between 2.5 nm and 10nm, depending on temperature, together with the appearance and progressive growth of protofibrillar aggregates. After these early aggregation stages, long and curved fibrils of up to several micrometers start to develop by elongation of the protofibrils. The calorimetric data indicate that the specific enthalpy of fibril disaggregation and unfolding is relatively low, suggesting a low density of interactions within the fibril structure as compared to the native protein and a main entropy contribution to the stability of the amyloid fibrils.     

Order of steps in the cytochrome C folding pathway: evidence for a sequential stabilization mechanism

Previous work used hydrogen exchange (HX) experiments in kinetic and equilibrium modes to study the reversible unfolding and refolding of cytochrome c (Cyt c) under native conditions. Accumulated results now show that Cyt c is composed of five individually cooperative folding units, called foldons, which unfold and refold as concerted units in a stepwise pathway sequence. The first three steps of the folding pathway are linear and sequential. The ordering of the last two steps has been unclear because the fast HX of the amino acid residues in these foldons has made measurement difficult. New HX experiments done under slower exchange conditions show that the final two foldons do not unfold and refold in an obligatory sequence. They unfold separately and neither unfolding obligately contains the other, as indicated by their similar unfolding surface exposure and the specific effects of destabilizing and stabilizing mutations, pH change, and oxidation state. These results taken together support a sequential stabilization mechanism in which folding occurs in the native context with prior native-like structure serving to template the stepwise formation of subsequent native-like foldon units. Where the native structure of Cyt c requires sequential folding, in the first three steps, this is found. Where structural determination is ambiguous, in the final two steps, alternative parallel folding is found.     

Rapid folding and unfolding of Apaf-1 CARD

Caspase recruitment domains (CARDs) are members of the death domain superfamily and contain six antiparallel helices in an alpha-helical Greek key topology. We have examined the equilibrium and kinetic folding of the CARD of Apaf-1 (apoptotic protease activating factor 1), which consists of 97 amino acid residues, at pH 6 and pH 8. The results showed that an apparent two state equilibrium mechanism is not adequate to describe the folding of Apaf-1 CARD at either pH, suggesting the presence of intermediates in equilibrium unfolding. Interestingly, the results showed that the secondary structure is less stable than the tertiary structure, based on the transition mid-points for unfolding. Single mixing and sequential mixing stopped-flow studies showed that Apaf-1 CARD folds and unfolds rapidly and suggest a folding mechanism that contains parallel channels with two unfolded conformations folding to the native conformation. Kinetic simulations show that a slow folding phase is described by a third conformation in the unfolded ensemble that interconverts with one or both unfolded species. Overall, the native ensemble is formed rapidly upon refolding. This is in contrast to other CARDs in which folding appears to be dominated by formation of kinetic traps.     

Dramatic stabilization of an SH3 domain by a single substitution: roles of the folded and unfolded states

The N-terminal SH3 domain of the Drosophila drk protein (drkN SH3) exists in equilibrium between folded and unfolded states under non-denaturing buffer conditions. In order to examine the origins of this instability, we have made mutations in the domain and characterized the thermodynamics and kinetics of folding. Results of substitutions of negatively charged residues to neutral amino acid residues suggest that the large electrostatic potential of the domain does not play a dominant role in the instability of the domain. Sequence alignment of a large number of SH3 domains reveals that the drkN SH3 domain has a threonine (T22) at a position corresponding to an otherwise highly conserved glycine residue in the diverging beta-turn connecting the beta3 and beta4 strands. Mutation of T22 to glycine results in significant stabilization of the drkN SH3 domain by 2.5 kcal/mole. To further characterize the basis for the stabilization of the T22 mutant relative to wild-type, we made additional mutant proteins with substitutions of residue T22. A strong correlation is seen between protein stability or folding rate and propensity for native beta-turn structure at this position. Correlation of folding rates with AGADIR predictions of non-native helical structure in the diverging turn region, along with our previous NMR evidence for non-native structure in this region of the unfolded state of the drkN SH3 domain, suggests that the free energy of the unfolded state also plays a role in stability. This result highlights the importance of both folded and unfolded states for understanding protein stability.     

Exploring subdomain cooperativity in T4 lysozyme II: uncovering the C-terminal subdomain as a hidden intermediate in the kinetic folding pathway

Intermediates along a protein's folding pathway can play an important role in its biology. Previous kinetics studies have revealed an early folding intermediate for T4 lysozyme, a small, well-characterized protein composed of an N-terminal and a C-terminal subdomain. Pulse-labeling hydrogen exchange studies suggest that residues from both subdomains contribute to the structure of this intermediate. On the other hand, equilibrium native state hydrogen experiments have revealed a high-energy, partially unfolded form of the protein that has an unstructured N-terminal subdomain and a structured C-terminal subdomain. To resolve this discrepancy between kinetics and equilibrium data, we performed detailed kinetics analyses of the folding and unfolding pathways of T4 lysozyme, as well as several point mutants and large-scale variants. The data support the argument for the presence of two distinct intermediates, one present on each side of the rate-limiting transition state barrier. The effects of circular permutation and site-specific mutations in the wild-type and circular permutant background, as well as a fragment containing just the C-terminal subdomain, support a model for the unfolding intermediate with an unfolded N-terminal and a folded C-terminal subdomain. Our results suggest that the partially unfolded form identified by native state hydrogen exchange resides on the folded side of the rate-limiting transition state and is, therefore, under most conditions, a "hidden" intermediate.