A statistical mechanical model for hydrogen exchange in globular proteins

Proteins' Evolution upon Point Mutations

As the reader must be already aware, state-of-the-art protein folding prediction methods have reached a smashing success in their goal of accurately determining the three-dimensional structures of proteins. Yet, a solution to simple problems such as the effects of protein point mutations on their (i) native conformation; (ii) marginal stability; (iii) ensemble of high-energy nativelike conformations; and (iv) metamorphism propensity and, hence, their evolvability, remains as an unsolved problem. As a plausible solution to the latter, some properties of the amide hydrogen-deuterium exchange, a highly sensitive probe of the structure, stability, and folding of proteins, are assessed from a new perspective. The preliminary results indicate that the protein marginal stability change upon point mutations provides the necessary and sufficient information to estimate, through a Boltzmann factor, the evolution of the amide hydrogen exchange protection factors and, consequently, that of the ensemble of folded conformations coexisting with the native state. This work contributes to our general understanding of the effects of point mutations on proteins and may spur significant progress in our efforts to develop methods to determine the appearance of new folds and functions accurately.

Protecting thermodynamic stability of protein: The basic paradigm against stress and unfolded protein response by osmolytes

Organic osmolytes are known to play important role in stress protection by stabilizing macromolecules and suppressing harmful effects on functional activity. There is existence of several reports in the literature regarding their effects on structural, functional and thermodynamic aspects of many enzymes and the interaction parameters with proteins have been explored. Osmolytes are compatible with enzyme function and therefore, can be accumulated up to several millimolar concentrations. From the thermodynamic point of view, osmolyte raises mid-point of thermal denaturation (T_m) of proteins while having no significant effect on ΔG_D° (free energy change at physiological condition). Unfavorable interaction with the peptide backbone due to preferential hydration is the major driving force for folding of unfolded polypeptide in presence of osmolyte. However, the thermodynamic basis of stress protection and origin of compatibility paradigm has been a debatable issue. In the present manuscript, we attempt to elaborate the origin of stress protection and compatibility paradigm of osmolytes based on the effect on thermodynamic stability of proteins. We also infer that protective effects of osmolytes on ΔG_D° (of proteins) could also indicate its potential involvement in unfolded protein response and overall stress biology on macromolecular level.

Interpreting Hydrogen-Deuterium Exchange Events in Proteins Using Atomistic Simulations: Case Studies on Regulators of G-Protein Signaling Proteins

Hydrogen-deuterium exchange (HDX) experiments are widely used in studies of protein dynamics. To predict the propensity of amide hydrogens for exchange with deuterium, several models have been reported in which computations of amide-hydrogen protection factors are carried out using molecular dynamics (MD) simulations. Given significant variation in the criteria used in different models, the robustness and broader applicability of these models to other proteins, especially homologous proteins showing distinct amide-exchange patterns, remains unknown. The sensitivity of the predictions when MD simulations are conducted with different force-fields is yet to tested and quantified. Using MD simulations and experimental HDX data on three homologous signaling proteins, we report detailed studies quantifying the performance of seven previously reported models (M1-M7) of two general types: empirical and fractional-population models. We find that the empirical models show inconsistent predictions but predictions of the fractional population models are robust. Contrary to previously reported work, we find that the solvent-accessible surface area of amide hydrogens is a useful metric when combined with a new metric defining the distances of amide hydrogens from the first polar atoms in proteins. On the basis of this, we report two new models, one empirical (M8) and one population-based (M9). We find strong protection of amide hydrogens from solvent exchange both within the stable helical motifs and also in the interhelical loops. We further observe that the exchange-competent states of amide hydrogens occur on the sub 100 ps time-scale via localized fluctuations, and such states among amides of a given protein do not appear to show any cooperativity or allosteric coupling.

Structural Dynamics of the GW182 Silencing Domain Including its RNA Recognition motif (RRM) Revealed by Hydrogen-Deuterium Exchange Mass Spectrometry

The human GW182 protein plays an essential role in micro(mi)RNA-dependent gene silencing. miRNA silencing is mediated, in part, by a GW182 C-terminal region called the silencing domain, which interacts with the poly(A) binding protein and the CCR4-NOT deadenylase complex to repress protein synthesis. Structural studies of this GW182 fragment are challenging due to its predicted intrinsically disordered character, except for its RRM domain. However, detailed insights into the properties of proteins containing disordered regions can be provided by hydrogen-deuterium exchange mass spectrometry (HDX/MS). In this work, we applied HDX/MS to define the structural state of the GW182 silencing domain. HDX/MS analysis revealed that this domain is clearly divided into a natively unstructured part, including the CCR4-NOT interacting motif 1, and a distinct RRM domain. The GW182 RRM has a very dynamic structure, since water molecules can penetrate the whole domain in 2 h. The finding of this high structural dynamics sheds new light on the RRM structure. Though this domain is one of the most frequently occurring canonical protein domains in eukaryotes, these results are - to our knowledge - the first HDX/MS characteristics of an RRM. The HDX/MS studies show also that the α2 helix of the RRM can display EX1 behavior after a freezing-thawing cycle. This means that the RRM structure is sensitive to environmental conditions and can change its conformation, which suggests that the state of the RRM containing proteins should be checked by HDX/MS in regard of the conformational uniformity. Graphical Abstract.

Foldamer hypothesis for the growth and sequence differentiation of prebiotic polymers

It is not known how life originated. It is thought that prebiotic processes were able to synthesize short random polymers. However, then, how do short-chain molecules spontaneously grow longer? Also, how would random chains grow more informational and become autocatalytic (i.e., increasing their own concentrations)? We study the folding and binding of random sequences of hydrophobic ([Formula: see text]) and polar ([Formula: see text]) monomers in a computational model. We find that even short hydrophobic polar (HP) chains can collapse into relatively compact structures, exposing hydrophobic surfaces. In this way, they act as primitive versions of today's protein catalysts, elongating other such HP polymers as ribosomes would now do. Such foldamer catalysts are shown to form an autocatalytic set, through which short chains grow into longer chains that have particular sequences. An attractive feature of this model is that it does not overconverge to a single solution; it gives ensembles that could further evolve under selection. This mechanism describes how specific sequences and conformations could contribute to the chemistry-to-biology (CTB) transition.

Non-equilibrium hydrogen exchange for determination of H-bond strength and water accessibility in solid proteins

We demonstrate measurement of non-equilibrium backbone amide hydrogen-deuterium exchange rates (HDX) for solid proteins. The target of this study are the slowly exchanging residues in solid samples, which are associated with stable secondary-structural elements of proteins. These hydrogen exchange processes escape methods measuring equilibrium exchange rates of faster processes. The method was applied to a micro-crystalline preparation of the SH3 domain of chicken α-spectrin. Therefore, from a 100% back-exchanged micro-crystalline protein preparation, the supernatant buffer was exchanged by a partially deuterated buffer to reach a final protonation level of approximately 20% before packing the sample in a 1.3 mm rotor. Tracking of the HN peak intensities for 2 weeks reports on site-specific hydrogen bond strength and also likely reflects water accessibility in a qualitative manner. H/D exchange can be directly determined for hydrogen-bonded amides using ¹H detection under fast magic angle spinning. This approach complements existing methods and provides the means to elucidate interesting site-specific characteristics for protein functionality in the solid state.

Estimation of Hydrogen-Exchange Protection Factors from MD Simulation Based on Amide Hydrogen Bonding Analysis

Hydrogen exchange (HX) studies have provided critical insight into our understanding of protein folding, structure, and dynamics. More recently, hydrogen exchange mass spectrometry (HX-MS) has become a widely applicable tool for HX studies. The interpretation of the wealth of data generated by HX-MS experiments as well as other HX methods would greatly benefit from the availability of exchange predictions derived from structures or models for comparison with experiment. Most reported computational HX modeling studies have employed solvent-accessible-surface-area based metrics in attempts to interpret HX data on the basis of structures or models. In this study, a computational HX-MS prediction method based on classification of the amide hydrogen bonding modes mimicking the local unfolding model is demonstrated. Analysis of the NH bonding configurations from molecular dynamics (MD) simulation snapshots is used to determine partitioning over bonded and nonbonded NH states and is directly mapped into a protection factor (PF) using a logistics growth function. Predicted PFs are then used for calculating deuteration values of peptides and compared with experimental data. Hydrogen exchange MS data for fatty acid synthase thioesterase (FAS-TE) collected for a range of pHs and temperatures was used for detailed evaluation of the approach. High correlation between prediction and experiment for observable fragment peptides is observed in the FAS-TE and additional benchmarking systems that included various apo/holo proteins for which literature data were available. In addition, it is shown that HX modeling can improve experimental resolution through decomposition of in-exchange curves into rate classes, which correlate with prediction from MD. Successful rate class decompositions provide further evidence that the presented approach captures the underlying physical processes correctly at the single residue level. This assessment is further strengthened in a comparison of residue resolved protection factor predictions for staphylococcal nuclease with NMR data, which was also used to compare prediction performance with other algorithms described in the literature. The demonstrated transferable and scalable MD based HX prediction approach adds significantly to the available tools for HX-MS data interpretation based on available structures and models.

Prediction of native-state hydrogen exchange from perfectly funneled energy landscapes

Simulations based on perfectly funneled energy landscapes often capture many of the kinetic features of protein folding. We examined whether simulations based on funneled energy functions can also describe fluctuations in native-state protein ensembles. We quantitatively compared the site-specific local stability determined from structure-based folding simulations, with hydrogen exchange protection factors measured experimentally for ubiquitin, chymotrypsin inhibitor 2, and staphylococcal nuclease. Different structural definitions for the open and closed states based on the number of native contacts for each residue, as well as the hydrogen-bonding state, or a combination of both criteria were evaluated. The predicted exchange patterns agree with the experiments under native conditions, indicating that protein topology indeed has a dominant effect on the exchange kinetics. Insights into the simplest mechanistic interpretation of the amide exchange process were thus obtained.

Strategies for the thermodynamic characterization of linked binding/local folding reactions within the native state application to the LID domain of adenylate kinase from Escherichia coli

Conformational fluctuations in proteins have emerged as an important aspect of biological function, having been linked to processes ranging from molecular recognition and catalysis to allostery and signal transduction. In spite of the realization of their importance, however, the connections between fluctuations and function have largely been empirical, even when they have been quantitative. Part of the problem in understanding the role of fluctuations in function is the fact that the mere existence of fluctuations complicates the interpretation of classic mutagenesis approaches. Namely, mutagenesis, which is typically targeted to an internal position (to elicit an effect), will change the fluctuations as well as the structure of the native state. Decoupling these effects is essential to an unambiguous understanding of the role of fluctuations in function. Here, we use a mutation strategy that targets surface-exposed sites in flexible parts of the molecule for mutation to glycine. Such mutations leave the ground-state structure unaffected. As a result, we can assess the nature of the fluctuations, develop a quantitative model relating fluctuations to function (in this case, molecular recognition), and unambiguously resolve the probabilities of the fluctuating states. We show that when this approach is applied to Escherichia coli adenylate kinase (AK), unique thermodynamic and structural insights are obtained, even when classic mutagenesis approaches targeted to the same region yield ambiguous results.

Neutralizing positive charges at the surface of a protein lowers its rate of amide hydrogen exchange without altering its structure or increasing its thermostability

This paper combines two techniques--mass spectrometry and protein charge ladders--to examine the relationship between the surface charge and hydrophobicity of a representative globular protein (bovine carbonic anhydrase II; BCA II) and its rate of amide hydrogen-deuterium (H/D) exchange. Mass spectrometric analysis indicated that the sequential acetylation of surface lysine-ε-NH3(+) groups--a type of modification that increases the net negative charge and hydrophobicity of the surface of BCA II without affecting its secondary or tertiary structure--resulted in a linear decrease in the aggregate rate of amide H/D exchange at pD 7.4, 15 °C. According to analysis with MS, the acetylation of each additional lysine generated between 1.4 and 0.9 additional hydrogens that are protected from H/D exchange during the 2 h exchange experiment at 15 °C, pD 7.4. NMR spectroscopy demonstrated that none of the hydrogen atoms which became protected upon acetylation were located on the side chain of the acetylated lysine residues (i.e., lys-ε-NHCOCH3) but were instead located on amide NHCO moieties in the backbone. The decrease in rate of exchange associated with acetylation paralleled a decrease in thermostability: the most slowly exchanging rungs of the charge ladder were the least thermostable (as measured by differential scanning calorimetry). This observation--that faster rates of exchange are associated with slower rates of denaturation--is contrary to the usual assumptions in protein chemistry. The fact that the rates of H/D exchange were similar for perbutyrated BCA II (e.g., [lys-ε-NHCO(CH2)2CH3]18) and peracetylated BCA II (e.g., [lys-ε-NHCOCH3]18) suggests that the electrostatic charge is more important than the hydrophobicity of surface groups in determining the rate of H/D exchange. These electrostatic effects on the kinetics of H/D exchange could complicate (or aid) the interpretation of experiments in which H/D exchange methods are used to probe the structural effects of non-isoelectric perturbations to proteins (i.e., phosphorylation, acetylation, or the binding of the protein to an oligonucleotide or to another charged ligand or protein).

Analysis of protein conformation and dynamics by hydrogen/deuterium exchange MS

Understanding as much as possible about proteins in the shortest amount of time has long been a goal of hydrogen exchange (HX) MS. Recent technological advances have led to improvements in the technique, but has this goal yet been achieved? (To listen to a podcast about this Feature, please go to the Analytical Chemistry Web site at pubs.acs.org/journal/ancham.).

NMR analysis of native-state protein conformational flexibility by hydrogen exchange

The rate of hydrogen exchange for the most protected amides of a protein is widely used to provide an estimate of global conformational stability by analyzing the exchange kinetics in the unfolded state in terms of model peptide exchange rates. The exchange behavior of the other amides of the protein which do not exchange via a global unfolding mechanism can provide insight into the smaller-scale conformational transitions that facilitate access to solvent as required for the exchange reaction. However, since the residual tertiary structure in the exchange-competent conformation can modulate the chemistry of the exchange reaction, equilibrium values estimated from normalization with model peptide rates are open to question. To overcome this limitation, the most robust approaches utilize differential analyses as a function of experimental variables such as denaturant concentration, temperature, pH, and mutational variation. Practical aspects of these various differential analysis techniques are considered with illustrations drawn from the literature.

The extremely slow-exchanging core and acid-denatured state of green fluorescent protein

Green fluorescent protein (GFP) is a large protein with a complex eleven-stranded beta-barrel structure. Previous studies have shown that it has a complex energy landscape for folding on which there are several intermediate states and a denatured state with significant residual structure. Here, we use two different types of HD exchange measurement and nuclear magnetic resonance (NMR) techniques to probe the energy landscape for folding of GFP in further detail. HD exchange experiments were performed over a wide range of conditions including different concentrations of denaturant. Results show that the penetration model dominates the exchange mechanism, consistent with the known stability and slow unfolding kinetics of GFP. HD exchange experiments at high pH establish that there is an extremely slow-exchanging superstable core of amide protons in GFP that are clustered and located in beta-strands 1, 2, 4, 5, and 6. These residues form part of a mini-beta-sheet which we propose constitutes a folding nucleus. Using a pulsed-labeling strategy, the acid-denatured state has been investigated and the residual structure observed in earlier studies shown to locate to beta-strands 1 and 3. There is some evidence that this residual structure is stabilized by a localized hydrophobic collapse of the polypeptide chain.

Denaturation of jack-bean urease by sodium n-dodecyl sulphate: A kinetic study below the critical micelle concentration

Kinetics of urease denaturation by anionic surfactant (sodium n-dodecyl sulphate, SDS) at concentrations below the critical micelle concentration (CMC) is investigated spectrophotometrically at neutral pH and the corresponding two-phase kinetic parameters of the process are estimated from a three-state reversible process using a binomial exponential relation based on the relaxation time method as: Using a prepared computer program, the experimental data are properly fitted into a binomial exponential relation, considering a two-phase denaturation pathway including a kinetically stable folded intermediate formed at SDS concentration of 1.1 mM. Forward and backward rate constants are estimated as: k(1)=0.2141+/-4.5 x 10(-3), k(2)=5.173 x 10(-3)+/-8.3 x 10(-5), k(-1)=0.09432+/-3.6 x 10(-4) and k(-2)=2.079 x 10(-3)+/-5.6 x 10(-5)s(-1) for the proposed mechanism. The rate-limiting step as well as the reaction coordinates in the denaturation mechanism are established. The mechanism involves formation of a kinetically stable folded native like intermediate through the electrostatic interactions. The intermediate was found to be more stable even than the native form (by about 9 kJmol(-1)) and still hexamer, because no loss of amplitude was observed. Electrophoresis experiments on the native and surfactant/urease complexes indicated a higher mobility for the kinetically folded native like intermediate.

Effects of pH, temperature, and sucrose on benzyl alcohol-induced aggregation of recombinant human granulocyte colony stimulating factor

Antimicrobial preservatives (e.g., benzyl alcohol), which are required in multidose formulations, can induce protein aggregation. In this study, the mechanism of benzyl alcohol-induced aggregation of recombinant human granulocyte colony-stimulating factor (rhGCSF) was investigated by determining the effects of temperature, pH, and sucrose on this process. rhGCSF was incubated at 25 and 37 degrees C and at pH 7.0 (phosphate-buffered saline, PBS) and pH 3.5 (HCl). Benzyl alcohol (0.9% w/v) accelerated aggregation of rhGCSF at pH 7.0, an effect that was much greater at 37 degrees C than at 25 degrees C and partially counteracted by 1.0 M sucrose. At pH 3.5, benzyl alcohol did not induce aggregation of rhGCSF. Spectroscopic studies showed that 0.9% benzyl alcohol altered the tertiary structure of rhGCSF at both pH, without detectably altering secondary structure. Structural perturbation was greater at 37 degrees C than at 25 degrees C. At both pH 7.0 and 3.5, the hydrogen-deuterium (H-D) exchange rate for rhGCSF was increased by 0.9% benzyl alcohol. Sucrose (1.0 M) partially counteracted the benzyl alcohol-induced perturbation of tertiary structure and the increase in H-D exchange rate. Thus, benzyl alcohol accelerates aggregation of rhGCSF at pH 7.0, because it favors partially unfolded aggregation-prone conformations of the protein. Sucrose partially counteracts benzyl alcohol-induced rhGCSF aggregation by shifting the molecular population away from these species and towards more compact conformations. We postulate that the absence of aggregation at pH 3.5, even with benzyl alcohol-induced structural perturbation, is due to the unfavorable energetics of intermolecular interactions (i.e., colloidal stability) between rhGCSF molecules at this pH.

Cooperative α-helix unfolding in a protein-DNA complex from hydrogen-deuterium exchange

We present experimental evidence for a cooperative unfolding transition of an alpha-helix in the lac repressor headpiece bound to a symmetric variant of the lac operator, as inferred from hydrogen-deuterium (H-D) exchange experiments monitored by NMR spectroscopy. In the EX1 limit, observed exchange rates become pH-independent and exclusively sensitive to local structure fluctuations that expose the amide proton HN to exchange. Close to this regime, we measured decay rates of individual backbone HN signals in D2O, and of their mutual HN-HN NOE by time-resolved two-dimensional (2D) NMR experiments. The data revealed correlated exchange at the center of the lac headpiece recognition helix, Val20-Val23, and suggested that the correlation breaks down at Val24, at the C terminus of the helix. A lower degree of correlation was observed for the exchange of Val9 and Ala10 at the center of helix 1, while no correlation was observed for Val38 and Glu39 at the center of helix 3. We conclude that HN exchange in the recognition helix and, to some extent, in helix 1 is a cooperative event involving the unfolding of these helices, whereas the HN exchange in helix 3 is dominated by random local structure fluctuations.

Structural interpretation of hydrogen exchange protection factors in proteins: characterization of the native state fluctuations of CI2

Protection factors obtained from equilibrium hydrogen exchange experiments are an important source of structural information on both native and nonnative states of proteins. We present a method for determining ensembles of protein structures by using hydrogen exchange data as restraints in molecular dynamics simulations in conjunction with an empirical force-field. The method is applied to determine the ensemble of structures representing the native state of chymotrypsin inhibitor 2 (CI2), including the rare, large fluctuations responsible for hydrogen exchange.

Hydrogen exchange mass spectrometry for the analysis of protein dynamics

Hydrogen exchange coupled to mass spectrometry (MS) has become a valuable analytical tool for the study of protein dynamics. By combining information about protein dynamics with more classical functional data, a more thorough understanding of protein function can be obtained. In many cases, protein dynamics are directly related to specific protein functions such as conformational changes during enzyme activation or protein movements during binding. The method is made possible because labile backbone hydrogens in a protein will exchange with deuterium atoms when the protein is placed in a D2O solution. The subsequent increase in protein mass over time is measured with high-resolution MS. The location of the deuterium incorporation is determined by monitoring deuterium incorporation in peptic fragments that are produced after the labeling reaction. In this review, we will summarize the general principles of the method, discuss the latest variations on the experimental protocol that probe different types of protein movements, and review other recent work and improvements in the field.

Hydrogen-deuterium exchange in bovine serum albumin protein monitored by Fourier transform infrared spectroscopy, part II: kinetic studies

The set of infrared spectra recorded at different levels of hydrogen-deuterium (H/D) exchange as a function of time were processed using spectral decomposition. The most precise information about H/D exchange of the NH groups of the protein backbone was retrieved by observing the intensity change of the nearest C=O stretching vibration. The H/D exchange at the protein backbone begins with NH groups bonded to C=O with a characteristic frequency of 1683 cm(-1). These amide groups were initially free and are the first to accept H-bonds from water during the hydration process. The NH groups, which are connected to C=O groups with a characteristic band at 1657 cm(-1), exchange at a slower rate. For both populations the pairs of comparable exchange rates were calculated with rate constants of 10(-3) min(-1), 0.014 min(-1), 0.0046 min(-1), and 0.09 min(-1). The appearance of two different exchange rates for each population is a consequence of the distinct exposure of particular molecular groups to the solvent. Two additional bands sensitive to exchange are attributed to NH bending modes in the side chains and are located at 1610 cm(-1) and 1585 cm(-1). These NH groups undergo H/D exchange at the beginning of the exchange with exchange rates of 0.019 min(-1) and 0.17 min(-1), respectively. The calculated exchange rate for hydrating water molecules is 0.037 min(-1). The distribution of water's exchange rate is extremely broad and covers almost the entire interval of the time-dependent experiment. From the efficiency of the exchange (96%) it is evident that some parts of the protein are completely forbidden to water molecules.

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.

Explicit-chain model of native-state hydrogen exchange: Implications for event ordering and cooperativity in protein folding

Native-state hydrogen exchange experiments on several proteins have revealed partially unfolded conformations with diverse stabilities. These equilibrium observations have been used to support kinetic arguments that folding proceeds via a sequential "pathway." This interpretative logic is evaluated here by analyzing the relationship between thermodynamic behavior and folding kinetics in a class of simplified lattice protein models. The chain models studied have varying degrees of cooperative interplay (coupling) between local helical conformational preference and favorable nonlocal interactions. When model cooperativity is high, as native conditions are weakened, "isotherms" of free energy of exchange for residues belonging to the same helix merge together before global unfolding. The point of merger depends on the model energetic favorability of the helix. This trend is similar to the corresponding experimental observations. Kinetically, we find that the ordering of helix formation in the very last stage of native core assembly tends to follow the stabilities of their converged isotherms. In a majority (but not all) of folding trajectories, the final assembly of helices that are thermodynamically more stable against exchange precedes that of helices that are less stable against exchange. These model features are in partial agreement with common experimental interpretations. However, the model results also underscore the ensemble nature of the folding process: the kinetics of helix formation is not a discrete, strictly "all-or-none" process as that envisioned by certain non-explicit-chain models. Helices generally undergo many cycles of partial formation and dissolution before their conformations are fixed in the final assembly stage of folding, a kinetic stage that takes up only approximately 2% of the average folding time in the present model; and the ordering of the helices' final assembly in some trajectories can be different from the dominant ordering stipulated by the exchange isotherms.

Structural effects of drying and rehydration for enzymes in soils: a kinetics-FTIR analysis of α-chymotrypsin adsorbed on montmorillonite

The effects of desiccation and rehydration cycles encountered by extracellular enzymes in soils are studied on -chymotrypsin adsorbed on montmorillonite. The controlled hygrometric FTIR cell used in this study enables to monitor drying and rehydration steps undergone by the -chymotrypsin-montmorillonite suspension or by the enzyme alone. Relative humidity (RH) determines the amount of deuterated water in the FTIR cell atmosphere. The molar water/protein ratio (W/P) as well as the conformational and solvation states of the enzyme have been determined using H/D exchange monitored by FTIR-transmission spectroscopy. When the W/P ratio decreases from 3500 to approximately 400, unfolding of beta-secondary structure in three different domains involves about 8% of the polypeptide backbone with respect to the most solvated states. Desiccation induces beta-unfolding, which opens channels allowing free vapor water molecules to diffuse into the enzyme at 15% RH. On drying to 0% RH, displacements of internal water (H2O) in the enzyme are demonstrated by reverse peptide isotopic exchanges (COND ==> CONH). Specific beta-structures, only formed in highly solvated states, sequester around 20 internal H2O molecules. Indeed, most of the unfolded secondary structures during the drying step are refolded at W/P approximately 1000 during rehydration. However, self-association hinders the recovery of the initial closed tertiary structure. The pD-dependent structural changes controlling inward and outward water diffusion are suppressed, whether the protein is initially in an adsorbed state or in solution. Changes in secondary structures encountered during desiccation/rehydration cycle are similar for the protein either free or in the adsorbed state. Thus domains that are unfolded by adsorption are not concerned by the desiccation/rehydration cycle.

Rare Fluctuations of Native Proteins Sampled by Equilibrium Hydrogen Exchange

We present a method for determining the ensembles of native protein structures that result from the large fluctuations of low probability revealed by hydrogen-exchange experiments. The measured protection factors are used to bias Monte Carlo simulations to sample the structures of the exchange competent species. The approach is illustrated by its application to the case of alpha-lactalbumin.

Effects of sucrose on conformational equilibria and fluctuations within the native-state ensemble of proteins

Osmolytes increase the thermodynamic conformational stability of proteins, shifting the equilibrium between native and denatured states to favor the native state. However, their effects on conformational equilibria within native-state ensembles of proteins remain controversial. We investigated the effects of sucrose, a model osmolyte, on conformational equilibria and fluctuations within the native-state ensembles of bovine pancreatic ribonuclease A and S and horse heart cytochrome c. In the presence of sucrose, the far- and near-UV circular dichroism spectra of all three native proteins were slightly altered and indicated that the sugar shifted the native-state ensemble toward species with more ordered, compact conformations, without detectable changes in secondary structural contents. Thermodynamic stability of the proteins, as measured by guanidine HCl-induced unfolding, increased in proportion to sucrose concentration. Native-state hydrogen exchange (HX) studies monitored by infrared spectroscopy showed that addition of 1 M sucrose reduced average HX rate constants at all degrees of exchange of the proteins, for which comparison could be made in the presence and absence of sucrose. Sucrose also increased the exchange-resistant core regions of the proteins. A coupling factor analysis relating the free energy of HX to the free energy of unfolding showed that sucrose had greater effects on large-scale than on small-scale fluctuations. These results indicate that the presence of sucrose shifts the conformational equilibria toward the most compact protein species within native-state ensembles, which can be explained by preferential exclusion of sucrose from the protein surface.

Protein-folding kinetics and mechanisms studied by pulse-labeling and mass spectrometry

The "protein-folding problem" refers to the question of how and why a denatured polypeptide chain can spontaneously fold into a compact and highly ordered conformation. The classical description of this process in terms of reaction pathways has been complemented by models that describe folding as a biased conformational diffusion on a multidimensional energy landscape. The identification and characterization of short-lived intermediates provide important insights into the mechanism of folding. Pulsed hydrogen/deuterium exchange (HDX) methods are among the most powerful tools for studying the properties of kinetic intermediates. Analysis of pulse-labeled proteins by mass spectrometry (MS) provides information that is complementary to that obtained in nuclear magnetic resonance (NMR) studies; NMR data represent an average of entire protein ensembles, whereas MS can detect co-existing protein species. MS-based pulse-labeling experiments can distinguish between folding scenarios that involve parallel pathways, and those where folding is channeled through obligatory intermediates. The proteolytic digestion/MS technique provides spatially resolved information on the HDX pattern of folding intermediates. This method is especially important for proteins that are too large to be studied by NMR. Although traditional pulsed HDX protocols are based on quench-flow techniques, it is also possible to use electrospray (ESI) MS to analyze the reaction mixture on-line and "quasi-instantaneously" after labeling. This approach allows short-lived protein conformations to be studied by their HDX level, their ESI charge-state distribution, and their ligand-binding state. Covalent labeling of free cysteinyl residues provides an alternative approach to pulsed HDX experiments. Another promising development is the use of synchrotron X-rays to induce oxidation at specific sites within a protein for studying their solvent accessibility during folding.

Protein hydrogen exchange mechanism: local fluctuations

Experiments were done to study the dynamic structural motions that determine protein hydrogen exchange (HX) behavior. The replacement of a solvent-exposed lysine residue with glycine (Lys8Gly) in a helix of recombinant cytochrome c does not perturb the native structure, but it entropically potentiates main-chain flexibility and thus can promote local distortional motions and large-scale unfolding. The mutation accelerates amide hydrogen exchange of the mutated residue by about 50-fold, neighboring residues in the same helix by less, and residues elsewhere in the protein not at all, except for Leu98, which registers the change in global stability. The pattern of HX changes shows that the coupled structural distortions that dominate exchange can be several residues in extent, but they expose to exchange only one amide NH at a time. This "local fluctuation" mode of hydrogen exchange may be generally recognized by disparate near-neighbor rates and a low dependence on destabilants (denaturant, temperature, pressure). In contrast, concerted unfolding reactions expose multiple neighboring amide NHs with very similar computed protection factors, and they show marked destabilant sensitivity. In both modes, ionic hydrogen exchange catalysts attack from the bulk solvent without diffusing through the protein matrix.

Thermodynamic and kinetic exploration of the energy landscape of Borrelia burgdorferi OspA by native-state hydrogen exchange

We report a native-state hydrogen-exchange (HX) method to simultaneously obtain both thermodynamic and kinetic information on the formation of multiple excited states in a folding energy landscape. Our method exploits the inherent dispersion and pH dependence of the intrinsic HX rates to cover both the EX2 (thermodynamic) and EX1 (kinetic) regimes. At each concentration of denaturant, HX measurements are performed over a range of pH values. Using this strategy, we dissected Borrelia burgdorferi OspA, a predominantly beta-sheet protein containing a unique single-layer beta-sheet, into five cooperative units and postulated excited states predominantly responsible for HX. More importantly, we determined the interconversion rates between these excited states and the native state. The use of both thermodynamic and kinetic information from native-state HX enabled us to construct a folding landscape of this 28kDa protein, including local minima and maxima, and to discriminate on-pathway and off-pathway intermediates. This method, which we term EX2/EX1 HX, should be a powerful tool for characterizing the complex folding mechanisms exhibited by the majority of proteins.

Fast-folding protein kinetics, hidden intermediates, and the sequential stabilization model

Do two-state proteins fold by pathways or funnels? Native-state hydrogen exchange experiments show discrete nonnative structures in equilibrium with the native state. These could be called hidden intermediates (HI) because their populations are small at equilibrium, and they are not detected in kinetic experiments. HIs have been invoked as disproof of funnel models, because funnel pictures appear to indicate (1) no specific sequences of events in folding; (2) a continuum, rather than a discrete ladder, of structures; and (3) smooth landscapes. In the present study, we solve the exact dynamics of a simple model. We find, instead, that the present microscopic model is indeed consistent with HIs and transition states, but such states occur in parallel, rather than along the single pathway predicted by the sequential stabilization model. At the microscopic level, we observe a huge multiplicity of trajectories. But at the macroscopic level, we observe two pathways of specific sequences of events that are relatively traditional except that they are in parallel, so there is not a single reaction coordinate. Using singular value decomposition, we show an accurate representation of the shapes of the model energy landscapes. They are highly complex funnels.

Propagation of a single destabilizing mutation throughout the Escherichia coli ribonuclease HI native state

A point mutation (I53A) in the core of Escherichia coli RNase H* is known to destabilize both the native conformation (DeltaG(UN)) and the kinetic intermediate (DeltaG(UI)) by 2 kcal/mole. Here, we have used native-state hydrogen deuterium exchange to ask how this destabilization is propagated throughout the molecule. Stability parameters were obtained for individual residues in I53A and compared with those from the wild-type protein. A destabilization of 2 kcal/mole was observed in residues in the core but was not detected in the periphery of the molecule. These results are consistent with the localized destabilization of the core observed in the early intermediate of the kinetic folding pathway, supporting the resemblance of this kinetic intermediate to the partially unfolded form detected in the native state at equilibrium. A thermodynamic cycle also shows no interaction between Ile 53 and a residue in the periphery. There is, however, an increase in the number of denaturant-independent exchange events in the periphery of I53A, showing that effects of the point mutation are communicated to regions outside the core, although perhaps not through changes in stability. In sum, this work shows that localized regions within a protein can be destabilized independently. Furthermore, it implies a correspondence between the kinetic intermediate and the equilibrium PUF, as the magnitude and localization of the destabilization are the same in both.

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.

Detecting structural changes in viral capsids by hydrogen exchange and mass spectrometry

Amide hydrogen exchange and mass spectrometry have been used to study the pH-induced structural changes in the capsid of brome mosaic virus (BMV). Capsid protein was labeled in a structurally sensitive way by incubating intact viral particles in D(2)O at pH 5.4 and 7.3. Deuterium levels in the intact coat protein and its proteolytic fragments were determined by mass spectrometry. The largest deuterium increases induced by structural alteration occurred in the regions around the quasi-threefold axes, which are located at the center of the asymmetric unit. The increased levels of deuterium indicate loosening of structure in these regions. This observation confirms the previously proposed swelling model for BMV and cowpea chlorotic mottle virus (CCMV) and is consistent with the structure of swollen CCMV recently determined by cryo-electron microscopy and image reconstruction. Structural changes in the extended N- and C-terminal arms were also detected and compared with the results obtained with other swollen plant viruses. This study demonstrates that protein fragmentation/amide hydrogen exchange is a useful tool for probing structural changes in viral capsids.

Changes in protein conformational mobility upon activation of extracellular regulated protein kinase-2 as detected by hydrogen exchange

Changes in protein mobility accompany changes in conformation during the trans-activation of enzymes; however, few studies exist that validate or characterize this behavior. In this study, amide hydrogen/deuterium exchange/mass spectrometry was used to probe the conformational flexibility of extracellular signal-regulated protein kinase-2 before and after activation by phosphorylation. The exchange data indicated that extracellular regulated protein kinase-2 activation caused altered backbone flexibility in addition to the conformational changes previously established by x-ray crystallography. The changes in flexibility occurred in regions involved in substrate binding and turnover, suggesting their importance in enzyme regulation.

Modeling protein density of states: additive hydrophobic effects are insufficient for calorimetric two-state cooperativity

A well-established experimental criterion for two-state thermodynamic cooperativity in protein folding is that the van't Hoff enthalpy DeltaH(vH) around the transition midpoint is equal, or very nearly so, to the calorimetric enthalpy DeltaH(cal) of the entire transition. This condition is satisfied by many small proteins. We use simple lattice models to provide a statistical mechanical framework to elucidate how this calorimetric two-state picture may be reconciled with the hierarchical multistate scenario emerging from recent hydrogen exchange experiments. We investigate the feasibility of using inverse Laplace transforms to recover the underlying density of states (i.e., enthalpy distribution) from calorimetric data. We find that the constraint imposed by DeltaH(vH)/DeltaH(cal) approximately 1 on densities of states of proteins is often more stringent than other "two-state" criteria proposed in recent theoretical studies. In conjunction with reasonable assumptions, the calorimetric two-state condition implies a narrow distribution of denatured-state enthalpies relative to the overall enthalpy difference between the native and the denatured conformations. This requirement does not always correlate with simple definitions of "sharpness" of a transition and has important ramifications for theoretical modeling. We find that protein models that assume capillarity cooperativity can exhibit overall calorimetric two-state-like behaviors. However, common heteropolymer models based on additive hydrophobic-like interactions, including highly specific two-dimensional Gō models, fail to produce proteinlike DeltaH(vH)/DeltaH(cal) approximately 1. A simple model is constructed to illustrate a proposed scenario in which physically plausible local and nonlocal cooperative terms, which mimic helical cooperativity and environment-dependent hydrogen bonding strength, can lead to thermodynamic behaviors closer to experiment. Our results suggest that proteinlike thermodynamic cooperativity may require a cooperative interplay between local and nonlocal interactions. The prospect of using calorimetric data to constrain Z-scores of knowledge-based potentials is discussed.

Ensemble modulation as an origin of denaturant-independent hydrogen exchange in proteins

Native state hydrogen exchange (HX) has become a powerful tool for the analysis of conformational states that exist under native conditions. However, the interpretation of HX data in terms of conformational fluctuations is still controversial. In particular, it has been shown that many residues display exchange behavior that is independent of denaturant concentration. It has been postulated that this lack of denaturant dependence results from local fluctuations that do not expose appreciable amounts of buried surface area. Here, we use a general thermodynamic description of HX to explore the different possibilities for this behavior. We find that the denaturant dependence seen in HX experiments under native conditions is not a de facto indication of the amount of surface area exposure required for exchange. Instead, this behavior results from the relatively homogenous character of the conformational ensemble that exists under native conditions and the non-specific nature of denaturant effects. Furthermore, a comparison of the HX behavior from a stabilized mutant of Staphylococcal nuclease (SNase) with that predicted for the wild-type SNase from the COREX algorithm suggests that denaturant-independent exchange of many residues is consistent with significant (approximately 10 %) surface area exposure for this protein.

Rate and equilibrium constants for protein unfolding and refolding determined by hydrogen exchange-mass spectrometry

Studies of protein unfolding and refolding may help us understand the more general problem of protein folding. Recent studies from this laboratory demonstrated that the unfolding and refolding of a large protein, rabbit muscle aldolase (M(r) 157 kDa), can be studied by combining amide hydrogen exchange and mass spectrometry. Results of these studies indicated that aldolase has three unfolding domains which likely unfold sequentially. Urea was used to increase the populations of partially unfolded states which were labeled with deuterium following a brief exposure to D(2)O. Electrospray ionization mass spectra of both the intact protein and its peptic fragments had multiple envelopes of isotope peaks from which the populations of unfolded forms were determined. The present study extends the previous investigations to include different urea concentrations and kinetic modeling of data taken as the system approaches equilibrium. Analysis of these results gives rate and equilibrium constants describing the unfolding and refolding processes characteristic of aldolase destabilized in urea. The change in solvent-accessible surface, which has been used as a reaction coordinate for protein folding, is estimated from the dependence of the equilibrium and rate constants on the concentration of urea.

Hydrogen exchange demonstrates three domains in aldolase unfold sequentially

Rabbit muscle aldolase is a homotetramer in which the subunits have a classical alpha/beta-barrel structure and Mr 39,212 Da. We have previously reported that aldolase incubated in 3 M urea has three unfolding domains distinguished by their different unfolding rates. The unfolding rates of these domains were determined from isotope patterns in the mass spectra of peptic fragments derived from aldolase incubated in 3 M urea and pulse labeled in (2)H2O. The present study extends this investigation to more thoroughly characterize the structures of these unfolding intermediates. Mass spectra of intact monomers had four envelopes of isotope peaks corresponding to four structural forms of aldolase. Analysis of the present results suggests that these structural forms consist of native aldolase and three forms that have one to three domains unfolded. The molecular masses of these four structural forms indicate that there are 107 residues in each of the three unfolding domains of aldolase. Present results also show that aldolase remains a tetramer in 4 M urea, even though hydrogen exchange and circular dichroism indicate that it has lost most of its secondary and tertiary structure. The abundances of unfolded domains, which were determined from mass spectra of either intact aldolase or its peptic fragments, were used to determine the abundances of specific, partially unfolded forms of aldolase. Kinetic modeling of the abundances of these structures suggests that these structures are formed sequentially as aldolase unfolds in urea.

The hydrogen exchange core and protein folding

A database of hydrogen-deuterium exchange results has been compiled for proteins for which there are published rates of out-exchange in the native state, protection against exchange during folding, and out-exchange in partially folded forms. The question of whether the slow exchange core is the folding core (Woodward C, 1993, Trends Biochem Sci 18:359-360) is reexamined in a detailed comparison of the specific amide protons (NHs) and the elements of secondary structure on which they are located. For each pulsed exchange or competition experiment, probe NHs are shown explicitly; the large number and broad distribution of probe NHs support the validity of comparing out-exchange with pulsed-exchange/competition experiments. There is a strong tendency for the same elements of secondary structure to carry NHs most protected in the native state, NHs first protected during folding, and NHs most protected in partially folded species. There is not a one-to-one correspondence of individual NHs. Proteins for which there are published data for native state out-exchange and theta values are also reviewed. The elements of secondary structure containing the slowest exchanging NHs in native proteins tend to contain side chains with high theta values or be connected to a turn/loop with high theta values. A definition for a protein core is proposed, and the implications for protein folding are discussed. Apparently, during folding and in the native state, nonlocal interactions between core sequences are favored more than other possible nonlocal interactions. Other studies of partially folded bovine pancreatic trypsin inhibitor (Barbar E, Barany G, Woodward C, 1995, Biochemistry 34:11423-11434; Barber E, Hare M, Daragan V, Barany G, Woodward C, 1998, Biochemistry 37:7822-7833), suggest that developing cores have site-specific energy barriers between microstates, one disordered, and the other(s) more ordered.

Comparison of continuous and pulsed labeling amide hydrogen exchange/mass spectrometry for studies of protein dynamics

In contrast to the rigid structures portrayed by X-ray diffraction, proteins in solution display constant motion which leads to populations that are momentarily unfolded. To begin to understand protein dynamics, we must have experimental methods for determining rates of folding and unfolding, as well as for identifying structures of folding and unfolding intermediates. Amide hydrogen exchange has become an important tool for such measurements. When urea is used to stabilize unfolded forms of proteins, the refolding rates may become slower than the rates of isotope exchange. In such cases, the intermolecular distribution of deuterium among the entire population of molecules may become bimodal, giving rise to a bimodal distribution of isotope peaks in mass spectra of the protein or its peptic fragments. When the protein is exposed continuously to D2O, the relative intensities of the two envelopes of isotope peaks give an integrated account of populations participating in the folding/unfolding process. However, when the protein is exposed only briefly to D2O, the relative intensities of the two envelopes of isotope peaks give an instantaneous measure of the folded/unfolded populations. Application of these two labeling methods to a large protein, aldolase, is described along with a discussion of specific parameters required to optimize these experiments.

Modeling deuterium exchange behavior of ERK2 using pepsin mapping to probe secondary structure

Recently, mass spectrometry has been applied to studies of hydrogen exchange of backbone amides, allowing analysis of large proteins at physiological concentrations. Low resolution spatial information is obtained by digesting proteins after exchange into D2O, using electrospray ionization liquid chromatography/mass spectrometry (ESI-LC/MS) to measure deuteration by mass increases of resulting peptides. This study develops modeling paradigms to increase resolution, using the signal transduction kinase ERK2 as a prototype for larger, less stable proteins. In-exchange data for peptides were analyzed by nonlinear least squares and a maximum entropy method, distinguishing amides into fast, intermediate, slow, and nonexchanging classes. Analysis of completely nonexchanging or in-exchanging peptides and peptides with sequence overlaps showed that nonexchanging amides were generally hydrogen bonded and sterically constrained or buried > or = 2.2 A from the protein surface, while fast exchanging hydrogens were generally exposed at the protein surface. In order to more fully understand the intermediate and slow exchanging classes, an empirical model was developed by analyzing published exchange rates in cytochrome c. The model correlated protection factors with a combined dependency on surface accessibility, hydrogen bond length, and position of residues from alpha helix ends. Together with analysis of partial proteolytic products, the derived rules for exchange allowed modeling of exchange behavior of peptides. Substantial deviation from the predicted rates in some cases suggested a role for conformational freedom in regulating fast and intermediate exchanging amides.

Experimental study of the protein folding landscape: unfolding reactions in cytochrome c

Hydrogen exchange results for cytochrome c have been interpreted in terms of transient hydrogen bond-breaking reactions that include large unfolding reactions and small fluctuational distortions. The differential sensitivity of these opening reactions to denaturant, temperature, and protein stability makes it possible to distinguish the different opening reactions and to characterize their structural and thermodynamic parameters. The partially unfolded forms (PUFs) observed are few and discrete, evidently because they are produced by the reversible unfolding of the protein's several intrinsically cooperative secondary structural elements. The PUFs are robust, evidently because the structural elements do not change over a wide range of conditions. The discrete nature of the PUFs and their small number is as expected for classical folding intermediates but not for theoretically derived folding models apparently because the simplified non-protein models usually analyzed in theoretical studies encompass only a single cooperative unit rather than multiple separable units.

Polymer principles and protein folding

This paper surveys the emerging role of statistical mechanics and polymer theory in protein folding. In the polymer perspective, the folding code is more a solvation code than a code of local phipsi propensities. The polymer perspective resolves two classic puzzles: (1) the Blind Watchmaker's Paradox that biological proteins could not have originated from random sequences, and (2) Levinthal's Paradox that the folded state of a protein cannot be found by random search. Both paradoxes are traditionally framed in terms of random unguided searches through vast spaces, and vastness is equated with impossibility. But both processes are partly guided. The searches are more akin to balls rolling down funnels than balls rolling aimlessly on flat surfaces. In both cases, the vastness of the search is largely irrelevant to the search time and success. These ideas are captured by energy and fitness landscapes. Energy landscapes give a language for bridging between microscopics and macroscopics, for relating folding kinetics to equilibrium fluctuations, and for developing new and faster computational search strategies.

Comparison of SH3 and SH2 domain dynamics when expressed alone or in an SH(3+2) construct: the role of protein dynamics in functional regulation

Protein dynamics play an important role in protein function and regulation of enzymatic activity. To determine how additional interactions with surrounding structure affects local protein dynamics, we have used hydrogen exchange and mass spectrometry to investigate the SH2 and SH3 domains of the protein tyrosine kinase Hck. Exchange rates of isolated Hck SH3 and SH2 domains were compared with rates for the same domains when part of a larger SH(3+2) construct. Increased deuterium incorporation was observed for the SH3 domain in the joint construct, particularly near the SH2 interface and the short sequence that connects SH3 to SH2, implying greater flexibility of SH3 when it is part of SH(3+2). Slow cooperative unfolding of the SH3 domain occurred at the same rate in isolated SH3 as in the SH(3+2) construct, suggesting a functional significance for this unfolding. The SH2 domain displayed relatively smaller changes in flexibility when part of the SH(3+2) construct. These results suggest that the domains influence each other. Further, our results imply a link between functional regulation and structural dynamics of SH3 and SH2 domains.

Native protein fluctuations: the conformational-motion temperature and the inverse correlation of protein flexibility with protein stability

We study the fluctuations of native proteins by exact enumeration using the HP lattice model. The model fluctuations increase with temperature. We observe a low-temperature point, below which large fluctuations are frozen out. This prediction is consistent with the observation by Tilton et al. [R. F. Tilton, Jr., J. C. Dewan, and G. A. Petsko, Biochemistry 31, 2469 (1992)], that the thermal motions of ribonuclease A increase sharply above about 200 K. We also explore protein "flexibility" as defined by Debye-Waller-like factors and solvent accessibilities of core residues to hydrogen exchange. We find that proteins having greater stability tend to have fewer large fluctuations, and hence lower flexibilities. If flexibility is necessary for enzyme catalysis, this could explain why proteins from thermophilic organisms, which are exceptionally stable, may be catalytically inactive at normal temperatures.

Structural dynamics of the Streptomyces lividans K+ channel (SKC1): secondary structure characterization from FTIR spectroscopy

Fourier transform infrared (FTIR) spectroscopy was used to probe the secondary structure, orientation, and the kinetics of amide hydrogen-deuterium exchange (HX) of the small K+ channel from Streptomyces lividans. Frequency component analysis of the amide I band showed that SKC1 is composed of 44-46% alpha-helix, 21-24% beta-sheet, 10-12% turns and 18-20% unordered structures. The order parameter S of the helical component of SKC1 was between 0.60 and 0.69. Close to 80% of SKC1 amide protons exchange within approximately 3 h of D2O exposure, suggesting that the channel is largely accessible to solvent exchange. These results are consistent with a model of SKC1 in which helices slightly tilted from the membrane normal line the water-filled vestibules that flank the K+ selectivity filter.