Native state EX2 and EX1 hydrogen exchange of Escherichia coli CspA, a small beta-sheet protein

Escherichia coli CspA is a small all-beta-sheet protein that folds fast (tau = 4 ms) via an apparent two-state mechanism. Our previous studies have shown that a large aromatic cluster on the surface of the protein participates in the rate-limiting step of folding and thus may be part of the folding nucleus of this protein. To obtain a more detailed picture of molecular events at the peptide backbone during unfolding and folding of CspA, we used native state hydrogen exchange and nuclear magnetic resonance spectroscopy (NMR). The experiments with native CspA were performed over a range of pH values from low pH, where exchange is governed by a rapid equilibrium before chemical exchange (EX2 exchange), to high pH, where exchange is dictated by the rate of unfolding (EX1 exchange). Rates of folding and unfolding were determined for 11 residues. The distribution of rates of folding within the structure of CspA suggests that hairpin turns, including one near the aromatic cluster, may nucleate the folding of CspA.     

Equilibrium amide hydrogen exchange and protein folding kinetics

The classical Linderstrøm-Lang hydrogen exchange (HX) model is extended to describe the relationship between the HX behaviors (EX1 and EX2) and protein folding kinetics for the amide protons that can only exchange by global unfolding in a three-state system including native (N), intermediate (I), and unfolded (U) states. For these slowly exchanging amide protons, it is shown that the existence of an intermediate (I) has no effect on the HX behavior in an off-pathway three-state system (IUN). On the other hand, in an on-pathway three-state system (UIN), the existence of a stable folding intermediate has profound effect on the HX behavior. It is shown that fast refolding from the unfolded state to the stable intermediate state alone does not guarantee EX2 behavior. The rate of refolding from the intermediate state to the native state also plays a crucial role in determining whether EX1 or EX2 behavior should occur. This is mainly due to the fact that only amide protons in the native state are observed in the hydrogen exchange experiment. These new concepts suggest that caution needs to be taken if one tries to derive the kinetic events of protein folding from equilibrium hydrogen exchange experiments.     

Kinetics of unfolding and folding from amide hydrogen exchange in native ubiquitin

Amide hydrogen (NH) exchange is one of the few experimental techniques with the potential for determining the thermodynamics and kinetics of conformational motions at nearly every residue in native proteins. Quantitative interpretation of NH exchange in terms of molecular motions relies on a simple two-state kinetic model: at any given slowly exchanging NH, a closed or exchange-incompetent conformation is in equilibrium with an open or exchange-competent conformation. Previous studies have demonstrated the accuracy of this model in measuring conformational equilibria by comparing exchange data with the thermodynamics of protein unfolding. We report here a test of the accuracy of the model in determining the kinetics of conformational changes in native proteins. The kinetics of folding and unfolding for ubiquitin have been measured by conventional methods and compared with those derived from a comprehensive analysis of the pH dependence of exchange in native ubiquitin. Rate constants for folding and unfolding from these two very different types of experiments show good agreement. The simple model for NH exchange thus appears to be a robust framework for obtaining quantitative information about molecular motions in native proteins.     

Partially folded states of staphylococcal nuclease highlight the conserved structural hierarchy of OB-fold proteins

We have been interested in whether three proteins that share a five-stranded beta-barrel "OB-fold" structural motif but no detectable sequence homology fold by similar mechanisms. Here we describe native-state hydrogen exchange experiments as a function of urea for SN (staphylococcal nuclease), a protein with an OB-fold motif and additional nonconserved elements of structure. The regions of structure with the largest stability and unfolding cooperativity are contained within the conserved OB-fold portion of SN, consistent with previous results for CspA (cold shock protein A) and LysN (anticodon binding domain of lysyl tRNA synthetase). The OB-fold also has the subset of residues with the slowest unfolding rates in the three proteins, as determined by hydrogen exchange experiments in the EX1 limit. Although the protein folding hierarchy is maintained at the level of supersecondary structure, it is not evident for individual residues as might be expected if folding depended on obligatory nucleation sites. Rather, the site-specific stability profiles appear to be linked to sequence hydrophobicity and to the density of long-range contacts at each site in the three-dimensional structures of the proteins. We discuss the implications of the correlation between stability to unfolding and conservation of structure for mechanisms of protein structure evolution.     

Beyond the EX1 limit: probing the structure of high-energy states in protein unfolding

Hydrogen exchange kinetics in native solvent conditions have been used to explore the conformational fluctuations of an immunoglobulin domain (CD2.domain1). The global folding/unfolding kinetics of the protein are unaltered between pH 4.5 and pH 9.5, allowing us to use the pH-dependence of amide hydrogen/deuterium exchange to characterise conformational states with energies up to 7.2kcal/mol higher than the folded ground state. The study was intended to search for discreet unfolding intermediates in this region of the energy spectrum, their presence being revealed by the concerted exchange behaviour of subsets of amide groups that become accessible at a given free energy, i.e. the spectrum would contain discreet groupings. Protection factors for 58 amide groups were measured across the pH range and the hydrogen-exchange energy profile is described.More interestingly, exchange behaviour could be grouped into three categories; the first two unremarkable, the third unexpected. (1) In 33 cases, amide exchange was dominated by rapid fluctuation, i.e. the free energy difference between the ground state and the rapidly accessed open state is sufficiently low that the contribution from crossing the unfolding barrier is negligible. (2) In 18 cases exchange is dominated by the global folding transition barrier across the whole pH range measured. The relationship between hydroxyl ion concentration and observed exchange rate is hyperbolic, with the limiting rate being that for global unfolding; the so-called EX1 limit. For these, the free energy difference between the folded ground state and any rapidly-accessed open state is too great for the proton to be exchanged through such fluctuations, even at the highest pH employed in this study. (3) For the third group, comprising five cases, we observe a behaviour that has not been described. In this group, as in category 2, the rate of exchange reaches a plateau; the EX1 limit. However, as the intrinsic exchange rate (k(int)) is increased, this limit is breached and the rate begins to rise again. This unintuitive behaviour does not result from pH instability, rather it is a consequence of amide groups experiencing two processes; rapid fluctuation of structure and crossing the global barrier for unfolding. The boundary at which the EX1 limit is overcome is determined by the equilibrium distribution of the fluctuating open and closed states (K(O/C)) and the rate constant for unfolding (k(u)). This critical boundary is reached when k(int)K(O/C)=k(u). Given that, in a simple transition state formalism: k(u)=K(#)k' (where K(#) describes the equilibrium distribution between the transition and ground state and k' describes the rate of a barrierless rearrangement), it follows that if the pH is raised to a level where k(int)=k', then the entire free energy spectrum from ground state to transition state could be sampled.     

Mapping protein energy landscapes with amide hydrogen exchange and mass spectrometry: I. A generalized model for a two-state protein and comparison with experiment

Protein amide hydrogen exchange (HDX) is a convoluted process, whose kinetics is determined by both dynamics of the protein and the intrinsic exchange rate of labile hydrogen atoms fully exposed to solvent. Both processes are influenced by a variety of intrinsic and extrinsic factors. A mathematical formalism initially developed to rationalize exchange kinetics of individual amide hydrogen atoms is now often used to interpret global exchange kinetics (e.g., as measured in HDX MS experiments). One particularly important advantage of HDX MS is direct visualization of various protein states by observing distinct protein ion populations with different levels of isotope labeling under conditions favoring correlated exchange (the so-called EX1 exchange mechanism). However, mildly denaturing conditions often lead to a situation where the overall HDX kinetics cannot be clearly classified as either EX1 or EX2. The goal of this work is to develop a framework for a generalized exchange model that takes into account multiple processes leading to amide hydrogen exchange, and does not require that the exchange proceed strictly via EX1 or EX2 kinetics. To achieve this goal, we use a probabilistic approach that assigns a transition probability and a residual protection to each equilibrium state of the protein. When applied to a small protein chymotrypsin inhibitor 2, the algorithm allows complex HDX patterns observed experimentally to be modeled with remarkably good fidelity. On the basis of the model we are now in a position to begin to extract quantitative dynamic information from convoluted exchange kinetics.     

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.     

Equilibrium hydrogen exchange reveals extensive hydrogen bonded secondary structure in the on-pathway intermediate of Im7

The four-helical immunity protein Im7 folds through an on-pathway intermediate that has a specific, but partially misfolded, hydrophobic core. In order to gain further insight into the structure of this species, we have identified the backbone hydrogen bonds formed in the ensemble by measuring the amide exchange rates (under EX2 conditions) of the wild-type protein and a variant, I72V. In this mutant the intermediate is significantly destabilised relative to the unfolded state (deltadeltaG(ui) = 4.4 kJ/mol) but the native state is only slightly destabilised (deltadeltaG(nu) = 1.8 kJ/mol) at 10 degrees C in 2H2O, pH* 7.0 containing 0.4 M Na2SO4, consistent with the view that this residue forms significant non-native stabilising interactions in the intermediate state. Comparison of the hydrogen exchange rates of the two proteins, therefore, enables the state from which hydrogen exchange occurs to be identified. The data show that amides in helices I, II and IV in both proteins exchange slowly with a free energy similar to that associated with global unfolding, suggesting that these helices form highly protected hydrogen-bonded helical structure in the intermediate. By contrast, amides in helix III exchange rapidly in both proteins. Importantly, the rate of exchange of amides in helix III are slowed substantially in the Im7* variant, I72V, compared with the wild-type protein, whilst other amides exchange more rapidly in the mutant protein, in accord with the kinetics of folding/unfolding measured using chevron analysis. These data demonstrate, therefore, that local fluctuations do not dominate the exchange mechanism and confirm that helix III does not form stable secondary structure in the intermediate. By combining these results with previously obtained Phi-values, we show that the on-pathway folding intermediate of Im7 contains extensive, stable hydrogen-bonded structure in helices I, II and IV, and that this structure is stabilised by both native and non-native interactions involving amino acid side-chains in these helices.     

Energetics of the native energy landscape of a two-domain calcium sensor protein: distinct folding features of the two domains

The protein folding energy landscape allows a thorough understanding of the protein folding problem which in turn helps in understanding various aspects of biological functions. Characterizing the cooperative unfolding units and the intermediates along the folding funnel of a protein is a challenging task. In this paper, we investigated the native energy landscape of EhCaBP, a calcium sensor, belonging to the same EF-hand superfamily as calmodulin. EhCaBP is a two-domain EF-hand protein consisting of two EF-hands in each domain and binding to four Ca2+ cations. Native-state hydrogen exchange (HX) was used to assess the folding features of the landscape and also to throw light on the structure-folding function paradigm of calcium sensor proteins. HX measurements under the EX2 regime provided the thermodynamic information about the protein folding events under native conditions. HX studies revealed that the unfolding of EhCaBP is not a two-state process. Instead, it proceeds through cooperative units. The C-terminal domain exhibits less denaturant dependence than the N-terminal domain, suggesting that the former is dominated by local fluctuations. It is interesting to note that the N- and C-terminal domains of EhCaBP have distinct folding features. In fact, these observed differences can regulate the domain-dependent target recognition of two-domain Ca2+ sensor proteins.     

Rapid cooperative two-state folding of a miniature alpha-beta protein and design of a thermostable variant

There is a great deal of interest in developing small stably folded miniature proteins. A limited number of these molecules have been described, however they typically have not been characterized in depth. In particular, almost no detailed studies of the thermodynamics and folding kinetics of these proteins have been reported. Here we describe detailed studies of the thermodynamics and kinetics of folding of a 39 residue mixed alpha-beta protein (NTL9(1-39)) derived from the N-terminal domain of the ribosomal protein L9. The protein folds cooperatively and rapidly in a two-state fashion to a native state typical of those found for normal globular proteins. At pH 5.4 in 20mM sodium acetate, 100mM NaCl the temperature of maximum stability is 6 degrees C, the t(m) is 65.3 degrees C, deltaH degrees (t(m)) is between 24.6 kcalmol(-1) and 26.3 kcalmol(-1), and deltaC(p) degrees is 0.38 kcalmol(-1)deg(-1). The thermodynamic parameters are in the range expected on the basis of per residue values determined from databases of globular proteins. H/2H exchange measurements reveal a set of amides that exchange via global unfolding, exactly as expected for a normal cooperatively folded globular protein. Kinetic measurements show that folding is two-state folding. The folding rate is 640 s(-1) and the value of deltaG degrees calculated from the folding and unfolding rates is in excellent agreement with the equilibrium value. A designed thermostable variant, generated by mutating K12 to M, was characterized and found to have a t(m) of 82 degrees C. Equilibrium and kinetic measurements demonstrate that its folding is cooperative and two-state.     

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.     

Stabilization of native protein fold by intein-mediated covalent cyclization

A mutant version of the N-terminal domain of Escherichia coli DnaB helicase was used as a model system to assess the stabilization against unfolding gained by covalent cyclization. Cyclization was achieved in vivo by formation of an amide bond between the N and C termini with the help of a split mini-intein. Linear and circular proteins were constructed to be identical in amino acid sequence. Mutagenesis of Phe102 to Glu rendered the protein monomeric even at high concentration. A difference in free energy of unfolding, DeltaDeltaG, between circular and linear protein of 2.3(+/-0.5) kcal mol(-1) was measured at 10 degrees C by circular dichroism. A theoretical estimate of the difference in conformational entropy of linear and circular random chains in a three-dimensional cubic lattice model predicted DeltaDeltaG=2.3 kcal mol(-1), suggesting that stabilization by protein cyclization is driven by the reduced conformational entropy of the unfolded state. Amide-proton exchange rates measured by NMR spectroscopy and mass spectrometry showed a uniform, approximately tenfold decrease of the exchange rates of the most slowly exchanging amide protons, demonstrating that cyclization globally decreases the unfolding rate of the protein. The amide proton exchange was found to follow EX1 kinetics at near-neutral pH, in agreement with an unusually slow refolding rate of less than 4 min(-1) measured by stopped-flow circular dichroism. The linear and circular proteins differed more in their unfolding than in their folding rates. Global unfolding of the N-terminal domain of E.coli DnaB is thus promoted strongly by spatial separation of the N and C termini, whereas their proximity is much less important for folding.     

Absence of kinetic thermal stabilization in a hyperthermophile rubredoxin indicated by 40 microsecond folding in the presence of irreversible denaturation

The striking kinetic stability of many proteins derived from hyperthermophilic organisms has led to the proposal that such stability may result from a heightened activation barrier for unfolding independent of a corresponding increase in the thermodynamic stability. This in turn implies a corresponding retardation of the folding reaction. A commonly cited model for kinetic thermal stabilization is the rubredoxin from Pyrococcus furiosus (Pf), which exhibits an irreversible denaturation lifetime at 100 degrees C of nearly a week. Utilizing protein resonances shifted well outside of the random coil chemical shift envelope, nuclear magnetic resonance (NMR) chemical exchange measurements on Pf rubredoxin as well as on the mesophile Clostridium pasteurianum (Cp) rubredoxin demonstrate reversible thermal transition temperatures of 144 degrees C (137 degrees C for the N-terminal modified A2K variant) and 104 degrees C, respectively, with similar (un)folding rates of approximately 25,000 s(-1), only modestly slower than the diffusion controlled rate. The absence of a substantial activation barrier to rubredoxin folding as well as the similar folding kinetics of the mesophile protein indicate that kinetic stabilization has not been utilized by the hyperthermophile rubredoxin in achieving its extreme thermal stability. The two-state folding kinetics observed for Pf rubredoxin contradict a previous assertion of multiphasic folding based on hydrogen exchange data extrapolated to an estimated midpoint of transition temperature (T(m)) of nearly 200 degrees C. This discrepancy is resolved by the observation that the base-catalyzed hydrogen exchange of the model dipeptide (N-acetyl-L-cysteine-N-methylamide)4-Cd2+ is 23-fold slower than that of the free cysteine model dipeptide used to normalize the Pf rubredoxin hydrogen exchange data.     

On the nature of conformational openings: native and unfolded-state hydrogen and thiol-disulfide exchange studies of ferric aquomyoglobin.

Native-state amide hydrogen exchange (HX) of proteins in the presence of denaturant has provided valuable details on the structures of equilibrium folding intermediates. Here, we extend HX theory to model thiol group exchange (SX) in single cysteine-containing variants of sperm whale ferric aquomyoglobin. SX is complementary to HX in that it monitors conformational opening events that expose side-chains, rather than the main chain, to solvent. A simple two-process model, consisting of EX2-limited local structural fluctuations and EX1-limited global unfolding, adequately accounts for all HX data. SX is described by the same model except at very low denaturant concentrations and when the bulky labeling reagent 5,5'-dithiobis (2-nitrobenzoic acid) (DTNB) is used. Under these conditions SX can occur by a novel denaturant-dependent process. This anomalous behavior is not observed when the smaller labeling reagent methyl methanethiosulfonate is employed, suggesting that it reflects a denaturant-induced increase in the amplitudes of local structural fluctuations. It also is not seen in heme-free apomyoglobin, which may indicate that local openings are sufficiently large in the absence of denaturant to allow DTNB unhindered access. Differences in SX kinetics obtained using the two labeling reagents provide estimates of the sizes of local opening reactions at different sites in the protein. At all sequence positions examined except for position 73, the same opening event appears to facilitate exchange of both backbone amide and side-chain thiol groups. The C73 thiol group is exposed by a low-energy fluctuation that does not expose its amide group to exchange.     

Conformational transitions in the membrane scaffold protein of phospholipid bilayer nanodiscs

Phospholipid bilayer nanodiscs are model membrane systems that provide an environment where membrane proteins are highly stable and monodisperse without the use of detergents or liposomes. Nanodiscs consist of a discoidal phospholipid bilayer encircled by two copies of an amphipathic alpha helical membrane scaffold protein, which is modeled from apolipoprotein A-1. Hydrogen exchange mass spectrometry was used to probe the structure and dynamics of the scaffold protein in the presence and absence of lipid. On nanodisc self-assembly, the entire scaffold protein gained significant protection from exchange, consistent with a large, protein-wide, structural rearrangement. This protection was short-lived and the scaffold protein was highly deuterated within 2 h. Several regions of the scaffold protein, in both the lipid-free and lipid-associated states, displayed EX1 unfolding kinetics. The rapid deuteration of the scaffold protein and the presence of correlated unfolding events both indicate that nanodiscs are dynamic rather than rigid bodies in solution. This work provides a catalog of the expected scaffold protein peptic peptides in a nanodisc-hydrogen exchange mass spectrometry experiment and their deuterium uptake signatures, data that can be used as a benchmark to verify correct assembly and nanodisc structure. Such reference data will be useful control data for all hydrogen exchange mass spectrometry experiments involving nanodiscs in which transmembrane or lipid-associated proteins are the primary molecule(s) of interest.     

Conformational dynamics of partially denatured myoglobin studied by time-resolved electrospray mass spectrometry with online hydrogen-deuterium exchange

This study demonstrates the use of electrospray mass spectrometry in conjunction with rapid online mixing ("time-resolved" ESI-MS) for monitoring protein conformational dynamics under equilibrium conditions. The hydrogen/deuterium exchange (HDX) kinetics of mildly denatured myoglobin (Mb) at pD 9.3, in the presence of 27% acetonitrile, were studied with millisecond time resolution. Analytical ultracentrifugation indicates that the average protein compactness under these solvent conditions is similar to that of native holomyoglobin (hMb). The mass spectrum shows protein ions in a wide array of charge and heme binding states, indicating the presence of multiple coexisting conformations. The experimental approach used allows the HDX kinetics of all of these species to be monitored separately. A combination of EX1 and EX2 behavior was observed for hMb ions in charge states 7+ to 9+, which predominantly represent nativelike hMb in solution. The EX1 kinetics are biphasic, indicating the presence of two protein populations that undergo conformational opening events with different rate constants. The EX2 kinetics observed for nativelike hMb are biphasic as well. All other charge and heme binding states represent non-native protein conformations that are involved in rapid interconversion processes, thus leading to monoexponential EX2 kinetics with a common rate constant. Burst phase labeling for these non-native proteins occurs at 125 sites. In contrast, the nativelike protein conformation shows burst phase labeling only for 88 sites. A kinetic model is developed which is based on the assumption of three distinct (un)folding units in Mb. The model implies that the free energy landscape of the protein exhibits a major barrier. The crossing of this barrier is most likely associated with slow, cooperative opening/closing events of the heme binding pocket. Rapid conformational fluctuations on either side of the barrier give rise to the observed EX2 kinetics. Simulated HDX kinetics based on this model are in excellent agreement with the experimental data.     

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.     

Dynamics of the Tec‐family tyrosine kinase SH3 domains

The Src Homology 3 (SH3) domain is an important regulatory domain found in many signaling proteins. X‐ray crystallography and NMR structures of SH3 domains are generally conserved but other studies indicate that protein flexibility and dynamics are not. We previously reported that based on hydrogen exchange mass spectrometry (HX MS) studies, there is variable flexibility and dynamics among the SH3 domains of the Src‐family tyrosine kinases and related proteins. Here we have extended our studies to the SH3 domains of the Tec family tyrosine kinases (Itk, Btk, Tec, Txk, Bmx). The SH3 domains of members of this family augment the variety in dynamics observed in previous SH3 domains. Txk and Bmx SH3 were found to be highly dynamic in solution by HX MS and Bmx was unstructured by NMR. Itk and Btk SH3 underwent a clear EX1 cooperative unfolding event, which was localized using pepsin digestion and mass spectrometry after hydrogen exchange labeling. The unfolding was localized to peptide regions that had been previously identified in the Src‐family and related protein SH3 domains, yet the kinetics of unfolding were not. Sequence alignment does not provide an easy explanation for the observed dynamics behavior, yet the similarity of location of EX1 unfolding suggests that higher‐order structural properties may play a role. While the exact reason for such dynamics is not clear, such motions can be exploited in intra‐ and intermolecular binding assays of proteins containing the domains.     

Conformational changes in chemically modified Escherichia coli thioredoxin monitored by H/D exchange and electrospray ionization mass spectrometry

Hydrogen/deuterium (H/D) exchange in combination with electrospray ionization mass spectrometry and near-ultraviolet (UV) circular dichroism (CD) was used to study the conformational properties and thermal unfolding of Escherichia coli thioredoxin and its Cys32-alkylated derivatives in 1% acetic acid (pH 2.7). Thermal unfolding of oxidized (Oxi) and reduced (Red) -thioredoxin (TRX) and Cys-32-ethylglutathionyl (GS-ethyl-TRX) and Cys-32-ethylcysteinyl (Cys-ethyl-TRX), which are derivatives of Red-TRX, follow apparent EX1 kinetics as charge-state envelopes, H/D mass spectral exchange profiles, and near-UV CD appear to support a two-state folding/unfolding model. Minor mass peaks in the H/D exchange profiles and nonsuperimposable MS- and CD-derived melting curves, however, suggest the participation of unfolding intermediates leading to the conclusion that the two-state model is an oversimplification of the process. The relative stabilities as measured by melting temperatures by both CD and mass spectral charge states are, Oxi-TRX, GS-ethyl-TRX, Cys-ethyl-TRX, and Red-TRX. The introduction of the Cys-32-ethylglutathionyl group provides extra stabilization that results from additional hydrogen bonding interactions between the ethylglutathionyl group and the protein. Near-UV CD data show that the local environment near the active site is perturbed to almost an identical degree regardless of whether alkylation at Cys-32 is by the ethylglutathionyl group, or the smaller, nonhydrogen-bonding ethylcysteinyl group. Mass spectral data, however, indicate a tighter structure for GS-ethyl-TRX.     

Towards a consistent modeling of protein thermodynamic and kinetic cooperativity: how applicable is the transition state picture to folding and unfolding?

To what extent do general features of folding/unfolding kinetics of small globular proteins follow from their thermodynamic properties? To address this question, we investigate a new simplified protein chain model that embodies a cooperative interplay between local conformational preferences and hydrophobic burial. The present four-helix-bundle 55mer model exhibits protein-like calorimetric two-state cooperativity. It rationalizes native-state hydrogen exchange observations. Our analysis indicates that a coherent, self-consistent physical account of both the thermodynamic and kinetic properties of the model leads naturally to the concept of a native state ensemble that encompasses considerable conformational fluctuations. Such a multiple-conformation native state is seen to involve conformational states similar to those revealed by native-state hydrogen exchange. Many of these conformational states are predicted to lie below native baselines commonly used in interpreting calorimetric data. Folding and unfolding kinetics are studied under a range of intrachain interaction strengths as in experimental chevron plots. Kinetically determined transition midpoints match well with their thermodynamic counterparts. Kinetic relaxations are found to be essentially single-exponential over an extended range of model interaction strengths. This includes the entire unfolding regime and a significant part of a folding regime with a chevron rollover, as has been observed for real proteins that fold with non-two-state kinetics. The transition state picture of protein folding and unfolding is evaluated by comparing thermodynamic free energy profiles with actual kinetic rates. These analyses suggest that some chevron rollovers may arise from an internal frictional effect that increasingly impedes chain motions with more native conditions, rather than being caused by discrete deadtime folding intermediates or shifts of the transition state peak as previously posited.