Structural characterization of protein folding intermediates by proton magnetic resonance and hydrogen exchange

High Pressure CPMG and CEST Reveal That Cavity Position Dictates Distinct Dynamic Disorder in the PP32 Repeat Protein

Given the central role of conformational dynamics in protein function, it is essential to characterize the time scales and structures associated with these transitions. High pressure (HP) perturbation favors transitions to excited states because they typically occupy a smaller molar volume, thus facilitating characterization of conformational dynamics. Repeat proteins, with their straightforward architecture, provide good models for probing the sequence dependence of protein conformational dynamics. Investigations of chemical exchange by 15N CPMG relaxation dispersion analysis revealed that introduction of a cavity via substitution of isoleucine 7 by alanine in the N-terminal capping motif of the pp32 leucine-rich repeat protein leads to pressure-dependent conformational exchange detected on the 500 μs-2 ms CPMG time scale. Exchange amplitude decreased from the N- to C-terminus, revealing a gradient of conformational exchange across the protein. In contrast, introduction of a cavity in the central core of pp32 via the L60A mutation led to pressure-induced exchange on a slower (>2 ms) time scale detected by 15N-CEST analysis. Excited state 15N chemical shifts indicated that in the excited state detected by HP CEST, the N-terminal region is mostly unfolded, while the core retains native-like structure. These HP chemical exchange measurements reveal that cavity position dictates exchange on distinct time scales, highlighting the subtle, yet central role of sequence in determining protein conformational dynamics.

Complex Folding Landscape of Apomyoglobin at Acidic pH Revealed by Ultrafast Kinetic Analysis of Core Mutants

Under mildly acidic conditions (pH 4-4.5) apomyoglobin (apoMb) adopts a partially structured equilibrium state ( M-state) that structurally resembles a kinetic intermediate encountered at a late stage of folding to the native structure at neutral pH. We have previously reported that the M-state is formed rapidly (<1 ms) via a multistate process and thus offers a unique opportunity for exploring early stages of folding by both experimental and computational techniques. In order to gain structural insight into intermediates and barriers at the residue level, we studied the folding/unfolding kinetics of 12 apoMb mutants at pH 4.2 using fluorescence-detected ultrafast mixing techniques. Global analysis of the submillisecond folding/unfolding kinetics vs urea concentration for each variant, based on a sequential four-state mechanism ( U ⇔ I ⇔ L ⇔ M), allowed us to determine elementary rate constants and their dependence on urea concentration for most transitions. Comparison of the free energy diagrams constructed from the kinetic data of the mutants with that of wild-type apoMb yielded quantitative information on the effects of mutations on the free energy (ΔΔ G) of both intermediates and the first two kinetic barriers encountered during folding. Truncation of conserved aliphatic side chains on helices A, G, and H gives rise to a stepwise increase in ΔΔ G as the protein advances from U toward M, consistent with progressive stabilization of native-like contacts within the primary core of apoMb. Helix-helix contacts in the primary core contribute little to the first folding barrier ( U ⇔ I) and thus are not required for folding initiation but are critical for the stability of the late intermediate, L, and the M-state. Alanine substitution of hydrophobic residues at more peripheral helix-helix contact sites of the native structure, which are still absent or unstable in the M-state, shows both positive (destabilizing) and negative (stabilizing) ΔΔ G, indicating that non-native contacts are formed initially and weakened or lost as a result of subsequent structural rearrangement steps.

Interrupted Pressure-Jump NMR Experiments Reveal Resonances of On-Pathway Protein Folding Intermediate

Previous pressure-jump NMR experiments on a pressure-sensitized double mutant of ubiquitin showed evidence that its folding occurs via two parallel, comparably efficient pathways: a single barrier and a two-barrier pathway. An interrupted folding NMR experiment is introduced, where for a brief period the pressure is dropped to atmospheric conditions (1 bar), followed by a jump back to high pressure for signal detection. Conventional, forward sampling of the indirect dimension during the low-pressure period correlates the ¹⁵N or ¹³C' chemical shifts of the unfolded protein at 1 bar to the ¹H frequencies of both the unfolded and folded proteins at high pressure. Remarkably, sampling the data of the same experiment in the reverse direction yields the frequencies of proteins present at the end of the low-pressure interval, which include unfolded, intermediate, and folded species. Although the folding intermediate ¹⁵N shifts differ strongly from natively folded protein, its ¹³C' chemical shifts, which are more sensitive probes for secondary structure, closely match those of the folded protein and indicate that the folding intermediate must have a structure that is quite similar to the native state.

High-Resolution Mapping of a Repeat Protein Folding Free Energy Landscape

A complete description of the pathways and mechanisms of protein folding requires a detailed structural and energetic characterization of the conformational ensemble along the entire folding reaction coordinate. Simulations can provide this level of insight for small proteins. In contrast, with the exception of hydrogen exchange, which does not monitor folding directly, experimental studies of protein folding have not yielded such structural and energetic detail. NMR can provide residue specific atomic level structural information, but its implementation in protein folding studies using chemical or temperature perturbation is problematic. Here we present a highly detailed structural and energetic map of the entire folding landscape of the leucine-rich repeat protein, pp32 (Anp32), obtained by combining pressure-dependent site-specific ¹H-¹⁵N HSQC data with coarse-grained molecular dynamics simulations. The results obtained using this equilibrium approach demonstrate that the main barrier to folding of pp32 is quite broad and lies near the unfolded state, with structure apparent only in the C-terminal region. Significant deviation from two-state unfolding under pressure reveals an intermediate on the folded side of the main barrier in which the N-terminal region is disordered. A nonlinear temperature dependence of the population of this intermediate suggests a large heat capacity change associated with its formation. The combination of pressure, which favors the population of folding intermediates relative to chemical denaturants; NMR, which allows their observation; and constrained structure-based simulations yield unparalleled insight into protein folding mechanisms.

Conformational dynamics of insulin

We have exploited a prandial insulin analog to elucidate the underlying structure and dynamics of insulin as a monomer in solution. A model was provided by insulin lispro (the active component of Humalog(®); Eli Lilly and Co.). Whereas NMR-based modeling recapitulated structural relationships of insulin crystals (T-state protomers), dynamic anomalies were revealed by amide-proton exchange kinetics in D(2)O. Surprisingly, the majority of hydrogen bonds observed in crystal structures are only transiently maintained in solution, including key T-state-specific inter-chain contacts. Long-lived hydrogen bonds (as defined by global exchange kinetics) exist only at a subset of four α-helical sites (two per chain) flanking an internal disulfide bridge (cystine A20-B19); these sites map within the proposed folding nucleus of proinsulin. The anomalous flexibility of insulin otherwise spans its active surface and may facilitate receptor binding. Because conformational fluctuations promote the degradation of pharmaceutical formulations, we envisage that "dynamic re-engineering" of insulin may enable design of ultra-stable formulations for humanitarian use in the developing world.

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.

Complexity of aromatic ring-flip motions in proteins: Y97 ring dynamics in cytochrome c observed by cross-relaxation suppressed exchange NMR spectroscopy

Dynamics of large-amplitude conformational motions in proteins are complex and less understood, although these processes are intimately associated with structure, folding, stability, and function of proteins. Here, we use a large set of spectra obtained by cross-relaxation suppressed exchange NMR spectroscopy (EXSY) to study the 180 degrees flipping motion of the Y97 ring of horse ferricytochrome c as a function of near-physiological temperature in the 288-308 K range. With rising temperature, the ring-flip rate constant makes a continuous transition from Arrhenius to anti-Arrhenius behavior through a narrow Arrhenius-like zone. This behavior is seen not only for the native state of the protein, but also for native-like states generated by adding subdenaturing amounts of guanidine deuterochloride (GdnDCl). Moderately destabilizing concentrations of the denaturant (1.5 M GdnDCl) completely removes the Arrhenius-like feature from the temperature window employed. The Arrhenius to anti-Arrhenius transition can be explained by the heat capacity model where temperature strengthens ground state interactions, perhaps hydrophobic in nature. The effect of the denaturant may appear to arise from direct protein-denaturant interactions that are structure-stabilizing under subdenaturing conditions. The temperature distribution of rate constants under different stability conditions also suggests that the prefactor in Arrhenius-like relations is temperature dependent. Although the use of the transition state theory (TST) offers several challenges associated with data interpretation, the present results and a consideration of others published earlier provide evidence for complexity of ring-flip dynamics in proteins.

Destabilization of transthyretin by pathogenic mutations in the DE loop

Transthyretin single-amino-acid variants are responsible for familial amyloidotic polyneuropathy, in which transthyretin variants accumulate extracellularly in the form of fibrillar aggregates. We studied the structural stabilities of four transthyretin variants (L58H, L58R, T59K, and E61K), in which a positively charged amino acid is introduced in a loop region between the D- and E-strands. In addition to being located in the DE-loop, L58 and T59 are involved in the core of the transthyretin monomer. The L58H, L58R, and T59K substitutions destabilized transthyretin more than the E61K mutation did, indicating that transthyretin is substantially destabilized by the substitution of residues located in both the DE-loop and the monomer core. By utilizing hydrogen-deuterium exchange and nuclear magnetic resonance, we demonstrated that residues in the G-strand and the loop between the A- and B-strands were destabilized by these pathogenic mutations in the DE loop. At the quaternary structural level, the DE-loop mutations destabilized the dimer-dimer contact area, which may lead to transient dissociation into a dimer. Our results suggest that the destabilization of the dimer-dimer interface and the monomer core is important for the amyloidogenesis of transthyretin.

GuHCl and NaCl-dependent hydrogen exchange in MerP reveals a well-defined core with an unusual exchange pattern

We have analysed hydrogen exchange at amide groups to characterise the energy landscape of the 72 amino acid residue protein MerP. From the guanidine hydrochloride (GuHCl) dependence of exchange in the pre-transitional region we have determined free energy values of exchange (DeltaG(HX)) and corresponding m-values for individual amide protons. Detailed analysis of the exchange patterns indicates that for one set of amide protons there is a weak dependence on denaturant, indicating that the exchange is dominated by local fluctuations. For another set of amide protons a linear, but much stronger, denaturant dependence is observed. Notably, the plots of free energy of exchange versus [GuHCl] for 16 amide protons show pronounced upward curvature, and a close inspection of the structure shows that these residues form a well-defined core in the protein. The hydrogen exchange that was measured at various concentrations of NaCl shows an apparent selective stabilisation of this core. Detailed analysis of this exchange pattern indicates that it may originate from selective destabilisation of the unfolded state by guanidinium ions and/or selective stabilisation of the core in the native state by chloride ions.

Effects of co-operative ligand binding on protein amide NH hydrogen exchange

Amide protection factors have been determined from NMR measurements of hydrogen/deuterium amide NH exchange rates measured on assigned signals from Lactobacillus casei apo-DHFR and its binary and ternary complexes with trimethoprim (TMP), folinic acid and coenzymes (NADPH/NADP(+)). The substantial sizes of the residue-specific DeltaH and TDeltaS values for the opening/closing events in NH exchange for most of the measurable residues in apo-DHFR indicate that sub-global or global rather than local exchange mechanisms are usually involved. The amide groups of residues in helices and sheets are those most protected in apo-DHFR and its complexes, and the protection factors are generally related to the tightness of ligand binding. The effects of ligand binding that lead to changes in amide protection are not localised to specific binding sites but are spread throughout the structure via a network of intramolecular interactions. Although the increase in protein stability in the DHFR.TMP.NADPH complex involves increased ordering in the protein structure (requiring TDeltaS energy) this is recovered, to a large extent, by the stronger binding (enthalpic DeltaH) interactions made possible by the reduced motion in the protein. The ligand-induced protection effects in the ternary complexes DHFR.TMP.NADPH (large positive binding co-operativity) and DHFR.folinic acid.NADPH (large negative binding co-operativity) mirror the co-operative effects seen in the ligand binding. For the DHFR.TMP.NADPH complex, the ligand-induced protection factors result in DeltaDeltaG(o) values for many residues being larger than the DeltaDeltaG(o) values in the corresponding binary complexes. In contrast, for DHFR.folinic acid.NADPH, the DeltaDeltaG(o) values are generally smaller than many of those in the corresponding binary complexes. The results indicate that changes in protein conformational flexibility on formation of the ligand complex play an important role in determining the co-operativity in the ligand binding.

Protein conformations, interactions, and H/D exchange

Modern mass spectrometry (MS) is well known for its exquisite sensitivity in probing the covalent structure of macromolecules, and for that reason, it has become the major tool used to identify individual proteins in proteomics studies. This use of MS is now widespread and routine. In addition to this application of MS, a handful of laboratories are developing and using a methodology by which MS can be used to probe protein conformation and dynamics. This application involves using MS to analyze amide hydrogen/deuterium (H/D) content from exchange experiments. Introduced by Linderstøm-Lang in the 1950s, H/D exchange involves using (2)H labeling to probe the rate at which protein backbone amide protons undergo chemical exchange with the protons of water. With the advent of highly sensitive electrospray ionization (ESI)-MS, a powerful new technique for measuring H/D exchange in proteins at unprecedented sensitivity levels also became available. Although it is still not routine, over the past decade the methodology has been developed and successfully applied to study various proteins and it has contributed to an understanding of the functional dynamics of those proteins.

Refolding processes of cytochrome P450cam from ferric and ferrous acid forms to the native conformation. Formations of folding intermediates with non-native heme coordination state

Changes in heme coordination state and protein conformation of cytochrome P450(cam) (P450(cam)), a b-type heme protein, were investigated by employing pH jump experiments coupled with time-resolved optical absorption, fluorescence, circular dichroism, and resonance Raman techniques. We found a partially unfolded form (acid form) of ferric P450(cam) at pH 2.5, in which a Cys(-)-heme coordination bond in the native conformation was ruptured. When the pH was raised to pH 7.5, the acid form refolded to the native conformation through a distinctive intermediate. Formations of similar acid and intermediate forms were also observed for ferrous P450(cam). Both the ferric and ferrous forms of the intermediate were found to have an unidentified axial ligand of the heme at the 6th coordination sphere, which is vacant in the high spin ferric and ferrous forms at the native conformation. For the ferrous form, it was also indicated that the 5th axial ligand is different from the native cysteinate. The folding intermediates identified in this study demonstrate occurrences of non-native coordination state of heme during the refolding processes of the large b-type heme protein, being akin to the well known folding intermediates of cytochromes c, in which c-type heme is covalently attached to a smaller protein.

NMR-detected hydrogen exchange and molecular dynamics simulations provide structural insight into fibril formation of prion protein fragment 106-126

PrP106-126, a peptide corresponding to residues 107-127 of the human prion protein, induces neuronal cell death by apoptosis and causes proliferation and hypertrophy of glia, reproducing the main neuropathological features of prion-related transmissible spongiform encephalopathies, such as bovine spongiform encephalopathy and Creutzfeldt-Jakob disease. Although PrP106-126 has been shown to form amyloid-like fibrils in vitro, their structural properties have not been elucidated. Here, we investigate the conformational characteristics of a fibril-forming fragment of the mouse prion protein, MoPrP106-126, by using electron microscopy, CD spectroscopy, NMR-detected hydrogen-deuterium exchange measurements, and molecular dynamics simulations. The fibrils contain approximately 50% beta-sheet structure, and strong amide exchange protection is limited to the central portion of the peptide spanning the palindromic sequence VAGAAAAGAV. Molecular dynamics simulations indicate that MoPrP106-126 in water assumes a stable structure consisting of two four-stranded parallel beta-sheets that are tightly packed against each other by methyl-methyl interactions. Fibril formation involving polyalanine stacking is consistent with the experimental observations.

Sequential vs. parallel protein-folding mechanisms: experimental tests for complex folding reactions

The recent emphasis on rough energy landscapes for protein folding reactions by theoreticians, and the many observations of complex folding kinetics by experimentalists provide a rationale for a brief literature survey of various empirical approaches for validating the underlying mechanisms. The determination of the folding mechanism is a key step in defining the energy surface on which the folding reactions occurs and in interpreting the effects of amino acid replacements on this reaction. Case studies that illustrate methods for differentiating between sequential and parallel channel folding mechanisms are presented. The ultimate goal of such efforts is to understand how the one-dimensional information contained in the amino acid sequence is rapidly and efficiently translated into three-dimensional structure.

Partly folded states of members of the lysozyme/lactalbumin superfamily: a comparative study by circular dichroism spectroscopy and limited proteolysis

The partly folded states of protein members of the lysozyme (LYS)/alpha-lactalbumin (LA) superfamily have been analyzed by circular dichroism (CD) measurements and limited proteolysis experiments. Hen, horse, dog, and pigeon LYSs and bovine LA were used in the present study. These are related proteins of 123- to 129-amino-acid residues with similar three-dimensional structures but low similarity in amino acid sequences. Moreover, notable differences among them reside in their calcium-binding properties and capability to adopt partly folded states or molten globules in acid solution (A-state) or on depletion of calcium at neutral pH (apo-state). Far- and near-UV CD measurements revealed that although the structures of hen and dog LYS are rather stable in acid at pH 2.0 or at neutral pH in the absence of calcium, conformational transitions to various extents occur with all other LYS/LA proteins herewith investigated. The most significant perturbation of tertiary structure in acid was observed with bovine LA and LYS from horse milk and pigeon egg-white. Pepsin and proteinase K were used as proteolytic probes, because these proteases show broad substrate specificity, and therefore, their sites of proteolysis are dictated not by the specific amino acid sequence of the protein substrate but by its overall structure and dynamics. Although hen LYS at pH 2.0 was fully resistant to proteolysis by pepsin, the other members of the LYS/LA superfamily were cleaved at different rates at few sites of the polypeptide chain and thus producing rather large protein fragments. The apo-form of bovine LA, horse LYS, and pigeon LYS were attacked by proteinase K at pH 8.3, whereas dog and hen LYSs were resistant to proteolysis when reacted under identical experimental conditions. Briefly, it has been found that the proteolysis data correlate well with the extent of conformational transitions inferred from CD spectra and with existing structural informations regarding the proteins herewith investigated, mainly derived from NMR and hydrogen exchange measurements. The sites of initial proteolytic cleavages in the LYS variants occur at the level of the beta-subdomain (approximately chain region 34-57), in analogy to those observed with bovine LA. Proteolysis data are in agreement with the current view that the molten globule of the LYS/LA proteins is characterized by a structured alpha-domain and a largely disrupted beta-subdomain. Our results underscore the utility of the limited proteolysis approach for analyzing structure and dynamics of proteins, even if adopting an ensemble of dynamic states as in the molten globule.

High-symmetry coordination cages via self-assembly

We provide a summary of our results in three-dimensional, coordination-driven self-assembly based on the directional-bonding methodology, in which the stoichiometric mixing of complementary building blocks, with appropriate, predefined geometries, leads to targeted, nanoscopic cages. Using this motif, we have synthesized high-symmetry ensembles resembling the Platonic solids, such as dodecahedra, and the Archimedean solids, such as truncated tetrahedra and cuboctahedra, as well as other cages, like trigonal bipyramids, adamantanoids, and trigonal prisms. The synthesis and characterization of these compounds is discussed, as is some host-guest chemistry.

Studies of biomolecular conformations and conformational dynamics by mass spectrometry

In the post-genomic era, a wealth of structural information has been amassed for proteins from NMR and crystallography. However, static protein structures alone are not a sufficient description: knowledge of the dynamic nature of proteins is essential to understand their wide range of functions and behavior during the life cycle from synthesis to degradation. Furthermore, few proteins have the ability to act alone in the crowded cellular environment. Assemblies of multiple proteins governed by complex signaling pathways are often required for the tasks of target recognition, binding, transport, and function. Mass spectrometry has emerged over the past several years as a powerful tool to address many of these questions. Recent improvements in "soft" ionization techniques have enabled researchers to study proteins and biomolecular complexes, both directly and indirectly. Likewise, continuous improvements in instrumental design in recent years have resulted in a dramatic expansion of the m/z range and resolution, enabling observation of large multi-protein assemblies whose structures are retained in the gas phase. In this article, we discuss some of the mass spectrometric techniques applied to investigate the nature of the conformations and dynamical properties that govern protein function.

A kinetic folding intermediate probed by native state hydrogen exchange

Stopped-flow fluorescence studies on the N-terminal domain of rat CD2 (CD2.d1) have demonstrated that folding from the fully denatured state (U) proceeds via the transient accumulation of an apparent intermediate (I) in a so-called burst phase that precedes the rate-limiting transition leading to the native state (N). A previous pH-dependent equilibrium hydrogen exchange (HX) study identified a subset of amides in CD2.d1 which, under EX2 conditions, exchange from N with free energies greater than or equal to the free energy difference between the N and I states calculated from the stopped-flow data. Under EX1 conditions the rates of HX for these amides tend towards an asymptote that matches the global unfolding rate calculated from the stopped-flow data, suggesting that exchange for these amides requires traversing the N-to-I transition state barrier. Exchange for these amides presumably occurs from exchange-competent forms comprising the kinetic burst phase therefore. To explore this idea further, native state HX (NHX) data have been collected for CD2.d1 under EX2 conditions using denaturant concentrations which span either side of the denaturant concentration where, according to the stopped-flow data, the apparent U and I states are iso-energetic. The data fit to a two-component, sub-global (sg)/global (g) NHX mechanism, yielding Delta G and m value parameters (where the m value is a measure of hydrocarbon solvation). Regression analysis demonstrates that the (m(sg), Delta G(sg)) and (m(g), Delta G(g)) values calculated for this subset of amides correspond with those describing the kinetic burst phase transition. This result confirms the ability of the NHX technique to explore the structural and energetic properties of kinetic folding intermediates.

An amino acid code for protein folding

Direct structural information obtained for many proteins supports the following conclusions. The amino acid sequences of proteins can stabilize not only the final native state but also a small set of discrete partially folded native-like intermediates. Intermediates are formed in steps that use as units the cooperative secondary structural elements of the native protein. Earlier intermediates guide the addition of subsequent units in a process of sequential stabilization mediated by native-like tertiary interactions. The resulting stepwise self-assembly process automatically constructs a folding pathway, whether linear or branched. These conclusions are drawn mainly from hydrogen exchange-based methods, which can depict the structure of infinitesimally populated folding intermediates at equilibrium and kinetic intermediates with subsecond lifetimes. Other kinetic studies show that the polypeptide chain enters the folding pathway after an initial free-energy-uphill conformational search. The search culminates by finding a native-like topology that can support forward (native-like) folding in a free-energy-downhill manner. This condition automatically defines an initial transition state, the search for which sets the maximum possible (two-state) folding rate. It also extends the sequential stabilization strategy, which depends on a native-like context, to the first step in the folding process. Thus the native structure naturally generates its own folding pathway. The same amino acid code that translates into the final equilibrium native structure-by virtue of propensities, patterning, secondary structural cueing, and tertiary context-also produces its kinetic accessibility.

Structural equilibrium fluctuations in mesophilic and thermophilic alpha-amylase

By comparing a mesophilic alpha-amylase with its thermophilic homolog, we investigated the relationship between thermal stability and internal equilibrium fluctuations. Fourier transform infrared spectroscopy monitoring hydrogen/deuterium (H/D) exchange kinetics and incoherent neutron scattering measuring picosecond dynamics were used to study dynamic features of the folded state at room temperature. Fairly similar rates of slowly exchanging amide protons indicate about the same free energy of stabilization DeltaG(stab) for both enzymes at room temperature. With respect to motions on shorter time scales, the thermophilic enzyme is characterized by an unexpected higher structural flexibility as compared to the mesophilic counterpart. In particular, the picosecond dynamics revealed a higher degree of conformational freedom for the thermophilic alpha-amylase. The mechanism proposed for increasing thermal stability in the present case is characterized by entropic stabilization and by flattening of the curvature of DeltaG(stab) as a function of temperature.

Enhanced prediction accuracy of protein secondary structure using hydrogen exchange Fourier transform infrared spectroscopy

A novel equilibrium hydrogen exchange Fourier transform IR (HX-FTIR) spectroscopy method for predicting secondary structure content was employed using spectra obtained for a training set of 23 globular proteins. The IR bandshape and frequency changes resulting from controlled levels of H-D exchange were observed to be protein-dependent. Their analysis revealed these variations to be partly correlated to secondary structure. For each protein, a set of 6 spectra was measured with a systematic variation of the solvent H-D ratio and was subjected to factor analysis. The most significant component spectra for each protein, representing independent aspects of the spectral response to deuteration, were each subjected to a second factor analysis over the entire training set. Restricted multiple regression (RMR) analysis using the loadings of the principal components from 19 of these H-D analyses revealed an improvement in prediction accuracy compared with conventional bandshape-based analyses of FTIR data. Nearly a factor of 2 reduction in error for prediction of helix fractions was found using s1, the average spectral response for the H-D set. In some cases, significant error reduction for prediction of minor components was found using higher factors. Using the same analytical methods, prediction errors with this new deuteration-response-FTIR method were shown to be even better than those obtained by use of electronic circular dichroism (ECD) data for helix predictions and to be significantly lower for ECD-based sheet prediction, making these the best secondary structure predictions obtained with the RMR method. Tests of a limited variable selection scheme showed further improvements, consistent with previous results of this approach using ECD data.

Microsecond to minute dynamics revealed by EX1-type hydrogen exchange at nearly every backbone hydrogen bond in a native protein

A previous comprehensive analysis of the pH dependence of native-state amide hydrogen (NH) exchange in turkey ovomucoid third domain (OMTKY3) yielded apparent opening and closing rate constants (k(op) and k(cl)) at 14 NH groups involved in global conformational changes. This analysis has been extended to 18 additional slowly exchanging NH groups. Quench-flow experiments were performed to monitor NH exchange in native OMTKY3 from neutral to very alkaline pH ( approximately 12) conditions. Above pH 10 the mechanism of exchange switched from one governed by a rapid equilibrium preceding the chemistry of exchange (i.e. EX2 exchange), to one where exchange was limited by the rate of opening (i.e. EX1 exchange). Kinetics of solvent exposure are now known for nearly all backbone NH groups in native OMTKY3, yielding rate constants that span five orders of magnitude, 0.004 to 200 s(-1).

Hydrogen exchange in ribonuclease A and ribonuclease S: evidence for residual structure in the unfolded state under native conditions

Two-dimensional NMR spectroscopy has been used to monitor the exchange of backbone amide protons in ribonuclease A (RNase A) and its subtilisin-cleaved form, ribonuclease S (RNase S). Exchange measurements at two different pH values (5.4 and 6.0) show that the exchange process occurs according to the conditions of the EX2 limit. Differential scanning calorimetry measurements have been carried out in 2H2O under conditions analogous to those used in the NMR experiments in order to determine the values of DeltaCp, DeltaHu and Tm, corresponding to the thermal denaturation of both proteins. For the amide protons of a large number of residues in RNase A, the free energies at 25 degreesC for exchange competent unfolding processes are much lower than the calorimetric denaturation free energies, thus showing that exchange occurs through local fluctuations in the native state. For 20 other protons, the cleavage reaction had approximately the same effect on the exchange rate constants than on the equilibrium constant for unfolding, indicating that those protons exchange by global unfolding. There is a good agreement between the residues to which these protons belong and those involved in the putative folding nucleation site identified by quench-flow NMR studies. The unfolding free energies of the slowest exchanging protons, DeltaGex, as evaluated from exchange data, are much larger than the calorimetric free energies of unfolding, DeltaGu. Given the agreement between DeltaDeltaGex(A-S), the difference in free energy from exchange for a given proton of the two proteins, and DeltaDeltaGu(A-S), the difference in the calorimetric free energy of the two proteins, the discrepancy indicates that the intrinsic exchange rates in the unfolded state of those protons cannot be approximated by those measured in short unstructured peptides and, consequently, exchange for those protons in RNase A and S must occur through a rather structured denatured state.

Quench-flow experiments combined with mass spectrometry show apomyoglobin folds through and obligatory intermediate

Folding of apomyoglobin is characterized by formation of a compact intermediate that contains substantial helicity. To determine whether this intermediate is obligatory or whether the protein can fold directly into the native state via an alternate parallel pathway, we have combined quench-flow hydrogen-exchange pulse labeling techniques with electrospray ionization mass spectrometry. The mass spectra of apomyoglobin obtained at various refolding times suggest that apomyoglobin indeed folds through a single pathway containing an obligatory intermediate with a significant hydrogen-bonded secondary structure content.

The structural distribution of cooperative interactions in proteins: analysis of the native state ensemble

Cooperative interactions link the behavior of different amino acid residues within a protein molecule. As a result, the effects of chemical or physical perturbations to any given residue are propagated to other residues by an intricate network of interactions. Very often, amino acids "sense" the effects of perturbations occurring at very distant locations in the protein molecule. In these studies, we have investigated by computer simulation the structural distribution of those interactions. We show here that cooperative interactions are not intrinsically bi-directional and that different residues play different roles within the intricate network of interactions existing in a protein. The effect of a perturbation to residue j on residue k is not necessarily equal to the effect of the same perturbation to residue k on residue j. In this paper, we introduce a computer algorithm aimed at mapping the network of cooperative interactions within a protein. This algorithm exhaustively performs single site thermodynamic mutations to each residue in the protein and examines the effects of those mutations on the distribution of conformational states. The algorithm has been applied to three different proteins (lambda repressor fragment 6-85, chymotrypsin inhibitor 2, and barnase). This algorithm accounts well for the observed behavior of these proteins.

Evidence for barrier-limited protein folding kinetics on the microsecond time scale

Although important structural events in protein folding are known to occur on the submillisecond time scale, the limited time resolution of conventional kinetic methods has precluded direct observation of the initial collapse of the polypeptide chain. A continuous-flow capillary mixing method recently developed by us made it possible to account for the entire fluorescence change associated with refolding of cytochrome c from approximately 5-10(-5)-10(2) s, including the previously unresolved quenching of Trp 59 fluorescence (burst phase) indicative of the formation of compact states. The kinetics of folding exhibits a major exponential process with a time constant of approximately 50 micros, independent of initial conditions and heme ligation state, indicating that a common free energy barrier is encountered during the initial collapse of the polypeptide chain. The resulting loosely packed intermediate accumulates prior to the rate-limiting formation of specific tertiary interactions, confirming previous indications that folding involves at least two distinct stages.

A continuous-flow capillary mixing method to monitor reactions on the microsecond time scale

A continuous-flow capillary mixing apparatus, based on the original design of Regenfuss et al. (Regenfuss, P., R. M. Clegg, M. J. Fulwyler, F. J. Barrantes, and T. M. Jovin. 1985. Rev. Sci. Instrum. 56:283-290), has been developed with significant advances in mixer design, detection method and data analysis. To overcome the problems associated with the free-flowing jet used for observation in the original design (instability, optical artifacts due to scattering, poor definition of the geometry), the solution emerging from the capillary is injected directly into a flow-cell joined to the tip of the outer capillary via a ground-glass joint. The reaction kinetics are followed by measuring fluorescence versus distance downstream from the mixer, using an Hg(Xe) arc lamp for excitation and a digital camera with a UV-sensitized CCD detector for detection. Test reactions involving fluorescent dyes indicate that mixing is completed within 15 micros of its initiation and that the dead time of the measurement is 45 +/- 5 micros, which represents a >30-fold improvement in time resolution over conventional stopped-flow instruments. The high sensitivity and linearity of the CCD camera have been instrumental in obtaining artifact-free kinetic data over the time window from approximately 45 micros to a few milliseconds with signal-to-noise levels comparable to those of conventional methods. The scope of the method is discussed and illustrated with an example of a protein folding reaction.

Characterisation of urea-denatured states of an immunoglobulin superfamily domain by heteronuclear NMR

The structural and dynamic properties of an immunoglobulin superfamily domain (IgSF), Ig 18', have been characterised by NMR at 285 K, in the presence of 4.2 M and 6.0 M urea, respectively. Analysis of chemical shift deviations, 3JHNHalpha coupling constants, sequential NOE pattern, and 15N relaxation data reveals that although the two urea-denatured states are highly disordered, some local turn-like residual structures do exist. Moreover, some distinct differences between the properties of the two denatured states are observed. In 4.2 M urea-denatured Ig 18', regions 80-83 and 86-92 adopt turn-like conformations, furthermore, region 84-93 is involved in slow exchange processes that occur on a micro- to millisecond time-scale. In the 6.0 M urea-denatured state, these turn-like conformations are less occupied, and chemical exchange processes in region 84-93 are largely reduced. In contrast, region 32-36 has persistent turn-like structures in both urea-denatured states. Some correlation between the spectral density function at 0 frequency, Jeff(0), for the urea-denatured states and the secondary structure elements of the folded state have been observed. Except for the terminal regions, residues corresponding to beta-strands have higher Jeff(0) values compared to residues corresponding to loops. The characterisation and comparison of the two urea-denatured states highlight residues that possess properties that may be crucial for the initiation of folding of this domain.

Definition of amide protection factors for early kinetic intermediates in protein folding

Hydrogen-deuterium exchange experiments have been used previously to investigate the structures of well defined states of a given protein. These include the native state, the unfolded state, and any intermediates that can be stably populated at equilibrium. More recently, the hydrogen-deuterium exchange technique has been applied in kinetic labeling experiments to probe the structures of transiently formed intermediates on the kinetic folding pathway of a given protein. From these equilibrium and nonequilibrium studies, protection factors are usually obtained. These protection factors are defined as the ratio of the rate of exchange of a given backbone amide when it is in a fully solvent-exposed state (usually obtained from model peptides) to the rate of exchange of that amide in some state of the protein or in some intermediate on the folding pathway of the protein. This definition is straightforward for the case of equilibrium studies; however, it is less clear-cut for the case of transient kinetic intermediates. To clarify the concept for the case of burst-phase intermediates, we have introduced and mathematically defined two different types of protection factors: one is P struc, which is more related to the structure of the intermediate, and the other is P app, which is more related to the stability of the intermediate. Kinetic hydrogen-deuterium exchange data from disulfide-intact ribonuclease A and from cytochrome c are discussed to explain the use and implications of these two definitions.

Folding kinetics of the SH3 domain of PI3 kinase by real-time NMR combined with optical spectroscopy

The refolding kinetics of the chemically denatured SH3 domain of phosphatidylinositol 3'-kinase (PI3-SH3) have been monitored by real-time one-dimensional 1H NMR coupled with a variety of other biophysical techniques. These experiments indicate that the refolding kinetics of PI3-SH3 are biphasic. The slow phase (27 (+/- 8)% amplitude) is due to a population of substantially unfolded molecules with an incorrectly configured cis proline residue. The fast phase (73 (+/- 8)% amplitude) corresponds to the folding of protein molecules with proline residues in a trans configuration in the unfolded state. NMR experiments indicate that the first species populated after the initiation of folding exhibit poor chemical shift dispersion and have spectra very similar to that of the denatured protein in 8 M guanidine hydrochloride. Linear combinations of the first spectrum and of the spectrum of the native protein accurately reconstruct all of the spectra acquired during refolding. Consistent with this, native side-chain and backbone H alpha atom packing (NMR), secondary structure (far-UV circular dichroism), burial of aromatic residues (near-UV circular dichroism), intrinsic fluorescence and peptide binding activity are all recovered with effectively identical kinetics. Equilibrium unfolding and folding/unfolding kinetics yield, within experimental error, identical values for the free energy of unfolding (delta Gu-H2O = 3.38 kcal mol-1) and for the slope of the free energy of unfolding versus denaturant concentration (meq = 2.33 kcal mol-1 M-1). Together, these data provide strong evidence that PI3-SH3 folds without significant population of kinetic well-structured intermediates. That PI3-SH3 folds slowly (time constant 2.8 seconds in H2O at 20 degrees C) indicates that slow refolding is not always a consequence of kinetic traps but may be observed even when a protein appears to fold via a simple, two-state mechanism.

Hydrogen exchange and protein folding

Amide hydrogen-deuterium exchange is a sensitive probe of the structure, stability and dynamics of proteins. The significant increase in the number of small, model proteins that have been studied has allowed a better understanding of the structural fluctuations that lead to hydrogen exchange. Recent technical advances enable the methodology to be applied to the study of protein-protein interactions in much larger, more complex systems.

A biophysical study of integral membrane protein folding

In order to characterize the thermodynamic constraints on the process of integral membrane protein folding and assembly, we have conducted a biophysical dissection of the structure of bacteriorhodopsin (BR), a prototypical alpha-helical integral membrane protein. Seven polypeptides were synthesized, corresponding to each of the seven transmembrane alpha-helices in BR, and the structure of each individual polypeptide was characterized in reconstituted phospholipid vesicles. Five of the seven polypeptides form stable transmembrane alpha-helices in isolation from the remainder of the tertiary structure of BR. However, using our reconstitution protocols, the polypeptide corresponding to the F helix in BR does not form any stable secondary structure in reconstituted vesicles, and the polypeptide corresponding to the G helix forms a hyperstable beta-sheet structure with its strands oriented perpendicular to the plane of the membrane. [The polypeptide corresponding to the C helix spontaneously equilibrates in a pH-dependent manner between a transmembrane alpha-helical conformation, a peripherally bound nonhelical conformation, and a fully water soluble conformation; the conformational properties of this polypeptide are the subject of the accompanying paper: Hunt et al. (1997) Biochemistry 36, 15177-15192.] Our observations suggest that the folding of alpha-helical integral membrane proteins may proceed spontaneously. However, the preference for a non-native conformation exhibited by two of the polypeptides suggests that the formation of some transmembrane substructures could require external constraints such as the links between the helices, interactions with the rest of the protein, or the involvement of cellular chaperones or translocases. Our results also suggest a strategy for improving the thermodynamic stability of alpha-helical integral membrane proteins, a goal that could facilitate attempts to overexpress and/or refold them.

Conformational properties of the A-state of cytochrome c studied by hydrogen/deuterium exchange and electrospray mass spectrometry

Hydrogen/deuterium (H/D) exchange studies that were monitored by liquid chromatography-electrospray ionization mass spectrometry (LC-ESIMS) were used to obtain a structural description of the compact acid-denatured state of ferricytochrome c (A-state). Due to the very different solvent conditions necessary to generate the nonnative states, it was essential that after deuterium labeling the nonnative states were refolded to the native state to insure high reproducibility during sample preparation and LC-ESIMS analysis. Approximately 30% lower deuterium was found incorporated in the A-state compared to the acid-denatured (UA) state. The analysis of the width of the mass peak suggests that the distribution of conformers sampled in the A-state was relatively narrow and that the compactness of the A-state was much closer to that of the native state than to the acid-denatured state. The LC-ESIMS study of partially deuterium-labeled peptic fragments derived from the A-state conformer generated under H/D quenching conditions were interpreted in terms of a significant loss of structural integrity within amino acid region 22-46.

Protein folding and intermediates

With the exception of the discovery of the rate of formation of the earliest intermediates, there have been no major conceptual leaps in our understanding of protein folding reactions over the past two years. Rather, this period has seen an extension of two established techniques: first, mutational analysis combined with a kinetic definition of the energy landscape of the reaction; and second, the use of hydrogen/deuterium exchange of backbone amide groups combined with NMR. Owing to the application of these methods to a wider range of proteins, it is now possible to draw some general conclusions about the physical processes that direct a protein to its native fold.

Absence of a stable intermediate on the folding pathway of protein A

The B-domain of protein A has one of the simplest protein topologies, a three-helix bundle. Its folding has been studied as a model for elementary steps in the folding of larger proteins. Earlier studies suggested that folding might occur by way of a helical hairpin intermediate. Equilibrium hydrogen exchange measurements indicate that the C-terminal helical hairpin could be a potential folding intermediate. Kinetic refolding experiments were performed using stopped-flow circular dichroism and NMR hydrogen-deuterium exchange pulse labeling. Folding of the entire molecule is essentially complete within the 6 ms dead time of the quench-flow apparatus, indicating that the intermediate, if formed, progresses rapidly to the final folded state. Site-directed mutagenesis of the isoleucine residue at position 16 was used to generate a variant protein containing tryptophan (the 116 W mutant). The formation of the putative folding intermediate was expected to be favored in this mutant at the expense of the native folded form, due to predicted unfavorable steric interactions of the bulky tryptophan side chain in the folded state. The 116 W mutant refolds completely within the dead time of a stopped-flow fluorescence experiment. No partly folded intermediate could be detected by either kinetic or equilibrium measurements. Studies of peptide fragments suggest that the protein A sequence has an intrinsic propensity to form a helix II/helix III hairpin. However, its stability appears to be marginal (of the order of 1/2 kT) and it could not be an obligatory intermediate on a defined folding pathway. These results explicitly demonstrate that the protein A B domain folds extremely rapidly by an apparent two-state mechanism without formation of stable partly folded intermediates. Similar mechanisms may also be involved in the rapid folding of subdomains of larger proteins to form the compact molten globule intermediates that often accumulate during the folding process.

Kinetic role of early intermediates in protein folding

The traditional view that partly folded intermediates are important for directing a protein toward the native state has been challenged by the notion that proteins can intrinsically fold rapidly in a single step if kinetic complications due to slow conformational events are avoided. Intermediates that accumulate within the first few milliseconds of folding are, however, a common observation even for small single-domain proteins. Recent spectroscopic studies, coupled with quantitative kinetic analysis, suggest that folding is facilitated by the rapid formation of compact intermediates with some native-like structural features.

Insights into protein folding from NMR

NMR has emerged as an important tool for studies of protein folding because of the unique structural insights it can provide into many aspects of the folding process. Applications include measurements of kinetic folding events and structural characterization of folding intermediates, partly folded states, and unfolded states. Kinetic information on a time scale of milliseconds or longer can be obtained by real-time NMR experiments and by quench-flow hydrogen-exchange pulse labeling. Although NMR cannot provide direct information on the very rapid processes occurring during the earliest stages of protein folding, studies of isolated peptide fragments provide insights into likely protein folding initiation events. Multidimensional NMR techniques are providing new information on the structure and dynamics of protein folding intermediates and both partly folded and unfolded states.

Influence of N-cap mutations on the structure and stability of Escherichia coli HPr

This paper describes the effect of N-capping substitutions on the structure and stability of histidine-containing protein (HPr). We have used NMR spectroscopy and conformational stability studies to quantify changes in local and global free energy due to mutagenesis at Ser46, the N-cap for helix B in HPr. Previous NMR studies suggested that helix B of Escherichia coli HPr is dynamic as judged by the rate of exchange of amide protons with solvent. Ser46 was chosen because it is the site of regulatory phosphorylation in HPrs from Gram-positive bacteria, and mutation of this residue to an aspartic acid (S46D) in E. coli HPr (Gram-negative) also makes it a poor substrate in the bacterial phosphoenolpyruvate: sugar phosphotransferase system. Therefore, to understand the mechanism of inactivation of E. coli S46D HPr, as well as the effect of mutagenesis on protein stability, we have characterized three mutants of E. coli HPr: Ser46 has been mutated to an Asp, Asn, and Ala in S46D, S46N, and S46A HPrs, respectively. The results indicate that these N-cap replacements have a marked influence on helix B stability. The effect of mutagenesis on local stability is correlated to global unfolding of HPr. The ability of amino acids to stabilize helix B is Asp > Asn > Ser > Ala. In addition, since there are neither large-scale conformational changes nor detectable changes in the active site of S46D HPr, it is proposed that the loss of phosphotransfer activity of S46D HPr is due to unfavorable steric and/or electrostatic interactions of the Asp with enzyme I of the PTS.

The carboxy-terminal region of human interleukin-5 is essential for maintenance of tertiary structure but not for dimerization

The C-terminal region of interleukin-5 has previously been suggested to be important for biological activity [Mackenzie et al., (1991), Mol. Immunol. 28, 155-158; Kodama et al. (1991), Biochem. Biophys. Res. Commun. 178, 514-519]. We have investigated this region by making a series of truncation mutants. The proteins were expressed in Escherichia coli, purified from inclusion bodies, and were able to refold with the disulfide homodimeric topology typical of interleukin-5. Analysis of the truncated carboxy-terminal proteins in an interleukin-5-dependent proliferation assay on TF-1 cells showed a rapid loss of activity as the C-terminal was shortened by more than two amino acids. This loss of biological activity correlated with a drop in binding affinity to both the alpha chain of the receptor and the high-affinity complex consisting of the alpha and beta subunits. Analysis of the proteins by 1H-NMR showed that the truncated mutants have higher exchange rates with solvent, indicating a less rigid structure. The carboxy-terminal region is therefore necessary to maintain the stability of the four-helix bundle and to orient correctly the important residues of the fourth helix. Inspection of the structure determined by X-ray crystallography shows that Trp-110 acts as the major residue in anchoring the fourth helix.

Refolding and unfolding kinetics of the equilibrium folding intermediate of apomyoglobin

Little is known about the kinetic process in which stable intermediates in protein folding are formed: whether their folding is highly cooperative (two-state) or weakly cooperative is controversial. We report here that the folding and unfolding kinetics of the pH 4-stable intermediate (I1) of apomyoglobin are measurable, in the millisecond time range, when monitored by stopped-flow measurements of tryptophan fluorescence. The kinetics confirm that folding of I1 is strongly cooperative, but there is a burst phase (missing amplitude) in unfolding. If the faster steps in unfolding of I1 can be measured directly by suitable fast-reaction methods, they will give information about the nature of the folding transition.

Hydrogen exchange in the carbon monoxide complex of soybean leghemoglobin

Hydrogen/deuterium exchange rates for individual amide protons have been measured for the carbon monoxide complex of soybean leghemoglobin. Fast two-dimensional NOESY experiments were performed, with 5.2-min data-collection time for each spectrum, which made possible the measurement of NOE cross-peaks of relatively rapidly exchanging amide protons at early time points. Exchange rates were measured for 61 backbone amides, the protection factors were calculated to provide information on the packing and local stability of the protein. The data are consistent with the presence of transient cooperative local unfolding of helical segments. The B-, E-, G- and H-helices have extensive regions of slow-, medium- and fast-exchanging amide protons. For each of these helices, there is a progressive decrease in protection on moving from the helix center to the termini. This is consistent with a stable helix center, with dynamic fraying at the ends. Amide exchange from the A-helix and C-helix is rapid except in small local regions. The F-helix, which is located on the proximal side of the heme pocket and is well formed in solution as demonstrated by characteristic medium range NOE connectivities [Morikis, D. Lepre, C.A. & Wright, P.E. (1994) Eur. J. Biochem. 219, 611-626], exhibits fast exchange for all amide protons. The implied flexibility and low stability of the F-helix may be functionally important in facilitating movement of the helix upon ligand binding. Fast exchange has also been observed for all amide protons in the CE-loop and in turns, as expected for flexible or solvent exposed regions. A strong tertiary contact has been established between the A-, G- and H-helices by the presence of a slowly exchanging indole N epsilon H of Trp129.

Watching protein folding unfold

Direct NMR observation of a transient folding intermediate provides new evidence for the importance of molten globules as general intermediates in protein folding.

Stopped-flow NMR spectroscopy: real-time unfolding studies of 6-19F-tryptophan-labeled Escherichia coli dihydrofolate reductase

Escherichia coli dihydrofolate reductase (DHFR; EC 1.5.1.3) contains five tryptophan residues that have been replaced with 6-19F-tryptophan. The 19F NMR assignments are known in the native, unliganded form and the unfolded form. We have used these assignments with stopped-flow 19F NMR spectroscopy to investigate the behavior of specific regions of the protein in real time during urea-induced unfolding. The NMR data show that within 1.5 sec most of the intensities of the native 19F resonances of the protein are lost but only a fraction (approximately 20%) of the intensities of the unfolded resonances appears. We postulate that the early disappearance of the native resonances indicates that most of the protein rapidly forms an intermediate in which the side chains have considerable mobility. Stopped-flow far-UV circular dichroism measurements indicate that this intermediate retains native-like secondary structure. Eighty percent of the intensities of the NMR resonances assigned to the individual tryptophans in the unfolded state appear with similar rate constants (k approximately 0.14 sec-1), consistent with the major phase of unfolding observed by stopped-flow circular dichroism (representing 80% of total amplitude). These data imply that after formation of the intermediate, which appears to represent an expanded structural form, all regions of the protein unfold at the same rate. Stopped-flow measurements of the fluorescence and circular dichroism changes associated with the urea-induced unfolding show a fast phase (half-time of about 1 sec) representing 20% of the total amplitude in addition to the slow phase mentioned above. The NMR data show that approximately 20% of the total intensity for each of the unfolded tryptophan resonances is present at 1.5 sec, indicating that these two phases may represent the complete unfolding of the two different populations of the native protein.