Hydrogen exchange in native and alcohol forms of ubiquitin.

Ubiquitin adopts a non-native folded structure in 60% methanol solution at low pH. Two-dimensional nuclear magnetic resonance (2D NMR) was used to measure the hydrogen-exchange rates of backbone amide protons of ubiquitin in both native and methanol forms, and to characterize the structure of ubiquitin in the methanol state. Protection factors (the ratios of experimentally determined exchange rates to the rates calculated for an unfolded polypeptide) for protons in the native form of ubiquitin range from less than 10 to greater than 10(5). Most of the protons that are protected from exchange are located in regions of hydrogen-bonded secondary structure. The most strongly protected backbone amide protons are those of residues comprising the hydrophobic core. Hydrogen exchange from ubiquitin in methanol solution was too rapid to measure directly by 2D NMR, so a labeling scheme was employed, in which exchange with solvent occurred while the protein was in methanol solution. Exchange was quenched by dilution with aqueous buffer after the desired labeling time, and proton occupancies were measured by 1H NMR of the native form of the protein. Protection factors for protons in the methanol form of ubiquitin range from 2.6 to 42, with all protected protons located in hydrogen-bonded structure in the native form. Again, the most strongly protected protons are those of residues in the hydrophobic core. Comparison of the patterns of the hydrogen-exchange rates in the native and methanol forms indicates that almost all of the native secondary structure persists in the methanol form, but that it is almost uniformly destabilized by 4-6 kcal/mol.(ABSTRACT TRUNCATED AT 250 WORDS)     

Three-dimensional solution structure and stability of thioredoxin m from spinach.

Proton NMR spectral resonances of thioredoxin m from spinach have been assigned, and its solution structure has been determined on the basis of 1156 nuclear Overhauser effect- (NOE-) derived distance constraints by using restrained molecular dynamics calculations. The average pairwise root-mean-square deviation (RMSD) for the 25 best NMR structures for the backbone was 1.0 +/- 0.1, when the structurally well-defined residues were considered. The N- and C-terminal segments (1-13 and 118-119) and residues 41-49, comprising the active site, are highly disordered. At the time of concluding this work, a crystal structure of this protein was reported, in which thioredoxin m was found to crystallize as noncovalent dimers. Although the solution and crystal structures are very similar, no evidence was found about the existence of dimers in solution, thus confirming that dimerization is not needed for the regulatory activity of thioredoxin m. The spinach thioredoxin m does not unfold by heat in the range 25-85 degrees C, as revealed by thermal circular dichroic (CD) measurements. However, its unfolding free energy (9.1 +/- 0.8 kcal mol(-1), at pH 5.3 and 25 degrees C) could be determined by extrapolating the free energy values obtained at different concentrations of guanidinium chloride (GdmCl). The folding-unfolding process is two-state as indicated by the coincidence of the CD denaturation curves obtained at far and near UV. The H/D exchange behavior of backbone amide protons was analyzed. The slowest-exchanging protons, requiring a global-unfolding mechanism in order to exchange, are those from beta2, beta3, and beta4, the central strands of the beta-sheet, which constitute the main element of the core of the protein. The free energies obtained from exchange measurements of protons belonging to the alpha-helices are lower than those derived from GdmCl denaturation studies, indicating that those protons exchange by local-unfolding mechanisms.     

Apparent local stability of the secondary structure of Azotobacter vinelandii holoflavodoxin II as probed by hydrogen exchange: implications for redox potential regulation and flavodoxin folding.

As a first step to determine the folding pathway of a protein with an alpha/beta doubly wound topology, the 1H, 13C, and 15N backbone chemical shifts of Azotobacter vinelandii holoflavodoxin II (179 residues) have been determined using multidimensional NMR spectroscopy. Its secondary structure is shown to contain a five-stranded parallel beta-sheet (beta2-beta1-beta3-beta4-beta5) and five alpha-helices. Exchange rates for the individual amide protons of holoflavodoxin were determined using the hydrogen exchange method. The amide protons of 65 residues distributed throughout the structure of holoflavodoxin exchange slowly at pH* 6.2 [kex < 10(-5) s(-1)] and can be used as probes in future folding studies. Measured exchange rates relate to apparent local free energies for transient opening. We propose that the amide protons in the core of holoflavodoxin only exchange by global unfolding of the apo state of the protein. The results obtained are discussed with respect to their implications for flavodoxin folding and for modulation of the flavin redox potential by the apoprotein. We do not find any evidence that A. vinelandii holoflavodoxin II is divided into two subdomains based on its amide proton exchange rates, as opposed to what is found for the structurally but not sequentially homologous alpha/beta doubly wound protein Che Y.     

Dihydrofolate reductase: sequential resonance assignments using 2D and 3D NMR and secondary structure determination in solution

Three-dimensional (3D) heteronuclear NMR techniques have been used to make sequential 1H and 15N resonance assignments for most of the residues of Lactobacillus casei dihydrofolate reductase (DHFR), a monomeric protein of molecular mass 18,300 Da. A uniformly 15N-labeled sample of the protein was prepared and its complex with methotrexate (MTX) studied by 3D 15N/1H nuclear Overhauser-heteronuclear multiple quantum coherence (NOESY-HMQC), Hartmann-Hahn-heteronuclear multiple quantum coherence (HOHAHA-HMQC), and HMQC-NOESY-HMQC experiments. These experiments overcame most of the spectral overlap problems caused by chemical shift degeneracies in 2D spectra and allowed the 1H-1H through-space and through-bond connectivities to be identified unambiguously, leading to the resonance assignments. The novel HMQC-NOESY-HMQC experiment allows NOE cross peaks to be detected between NH protons even when their 1H chemical shifts are degenerate as long as the amide 15N chemical shifts are nondegenerate. The 3D experiments, in combination with conventional 2D NOESY, COSY, and HOHAHA experiments on unlabelled and selectively deuterated DHFR, provide backbone assignments for 146 of the 162 residues and side-chain assignments for 104 residues of the protein. Data from the NOE-based experiments and identification of the slowly exchanging amide protons provide detailed information about the secondary structure of the binary complex of the protein with methotrexate. Sequential NHi-NHi+1 NOEs define four regions with helical structure. Two of these regions, residues 44-49 and 79-89, correspond to within one amino acid to helices C and E in the crystal structure of the DHFR.methotrexate.NADPH complex [Bolin et al. (1982) J. Biol. Chem. 257, 13650-13662], while the NMR-determined helix formed by residues 26-35 is about one turn shorter at the N-terminus than helix B in the crystal structure, which spans residues 23-34. Similarly, the NMR-determined helical region comprising residues 102-110 is somewhat offset from the crystal structure's helix F, which encompasses residues 97-107. Regions of beta-sheet structure were characterized in the binary complex by strong alpha CHi-NHi+1 NOEs and by slowly exchanging amide protons. In addition, several long-range NOEs were identified linking together these stretches to form a beta-sheet. These elements align perfectly with corresponding elements in the crystal structure of the DHFR.methotrexate.NADPH complex, which contains an eight-stranded beta-sheet, indicating that the main body of the beta-sheet is preserved in the binary complex in solution.     

Study of amide proton exchange in 15N-enriched cryptogein using pH-dependent off-resonance ROESY-HSQC experiments

A new pH-dependent off-resonance ROESY-HSQC experiment has been used to characterize the degree of protection of the amide protons of cryptogein, a protein of the elicitin family, against solvent exchange. The study of the pH dependence of solvent-shielded amide protons in this protein reveals that the helices have different levels of stability. Two of the five helices exhibit strong protection of amide hydrogens against exchange with the solvent. By contrast, greater flexibility is observed in the other three helices, particularly in the C-terminal helix. These results provide information on the dynamic features of the protein and are consistent with the RMSD for the backbone atoms of residues involved in helical structures. In addition, the question of the flexibility in a hydrophobic cavity made of conserved residues, which represent a plausible binding site, is addressed by this method. Received: 30 July 1997 / Accepted: 11 November 1997     

Insight into the factors influencing the backbone dynamics of three homologous proteins,dendrotoxins I and K,and BPTI: FTIR and time-resolved fluorescence investigations

Attenuated total reflection Fourier transform infrared (ATR-FTIR) spectroscopy, combined with hydrogen/deuterium exchange technique and time-resolved fluorescence spectroscopy, has been used to investigate the changes in structure and dynamics that underlie the thermodynamic stability differences observed for three closely homologous proteins: dendrotoxins I and K, and bovine pancreatic trypsin inhibitor (BPTI). The experiments were performed on proteins under their native state and a modified form, obtained by selective reduction of a disulfide bond at the surface of the molecule, increasing slightly the backbone flexibility without changing the average structure. The data confirmed the high local as well as global rigidity of BPTI. In protein K, the exchange process was slow during the first 2 h of exchange, presumably reflecting a compact three-dimensional conformation, and then increased rapidly, the internal amide protons of the beta-strands exchanging 10-fold faster than in BPTI or protein I. The most probable destabilizing element was identified as Pro32, in the core of the beta-sheet. Protein I was found to present a 10% more expanded volume than protein K or BPTI, and there is a possible correlation between the resulting increased flexibility of the molecule and the lower thermodynamic stability observed for this protein. Interestingly, the interior amide protons of the beta-sheet structure were found to be as protected against exchange in protein I as in BPTI, suggesting that, although globally more flexible than that of Toxin K or BPTI, the structure of Toxin I could be locally quite rigid. The structural factors suspected to be responsible for the differences in internal flexibility of the two toxins could play a significant role in determining their functional properties.     

Nuclear magnetic resonance studies of the internal dynamics in Apo, (Cd2+)1 and (Ca2+)2 calbindin D9k. The rates of amide proton exchange with solvent.

The backbone dynamics of the EF-hand Ca(2+)-binding protein, calbindin D9k, has been investigated in the apo, (Cd2+)1 and (Ca2+)2 states by measuring the rate constants for amide proton exchange with solvent. 15N-1H correlation spectroscopy was utilized to follow direct 1H-->2H exchange of the slowly exchanging amide protons and to follow indirect proton exchange via saturation transfer from water to the rapidly exchanging amide protons. Plots of experimental rate constants versus intrinsic rate constants have been analyzed to give qualitative insight into the opening modes of the protein that lead to exchange. These results have been interpreted within the context of a progressive unfolding model, wherein hydrophobic interactions and metal chelation serve to anchor portions of the protein, thereby damping fluctuations and retarding amide proton exchange. The addition of Ca2+ or Cd2+ was found to retard the exchange of many amide protons observed to be in hydrogen-bonding environments in the crystal structure of the (Ca2+)2 state, but not of those amide protons that were not involved in hydrogen bonds. The largest changes in rate constant occur for residues in the ion-binding loops, with substantial effects also found for the adjacent residues in helices I, II and III, but not helix IV. The results are consistent with a reorganization of the hydrogen-bonding networks in the metal ion-binding loops, accompanied by a change in the conformation of helix IV, as metal ions are chelated. Further analysis of the results obtained for the three states of metal occupancy provides insight into the nature of the changes in conformational fluctuations induced by ion binding.     

Determination of the secondary structure of interleukin-8 by nuclear magnetic resonance spectroscopy

The solution conformation of interleukin-8 (IL-8), a small protein of 72 residues with a wide range of proinflammatory activities, has been investigated by two-dimensional NMR spectroscopy. The 1H-NMR spectrum of IL-8 is assigned in a sequential manner and regular elements of secondary structure are identified on the basis of a qualitative interpretation of the nuclear Overhauser, coupling constant and amide exchange data. The IL-8 monomer contains a triple stranded anti-parallel beta-sheet arranged in a Greek key and a long C-terminal helix (residues 57-72). It is shown that IL-8 is a dimer in solution in which the interface is principally formed by six backbone hydrogen bonds between residues 25, 27, and 29 of one monomer and residues 29, 27, and 25, respectively, of the other. As a result, the two units of the dimer form a contiguous six-stranded anti-parallel beta-sheet. The secondary structure of IL-8 is similar to that found in the crystal structure of the sequence related protein platelet factor 4.     

The three-dimensional structure of guanine-specific ribonuclease F1 in solution determined by NMR spectroscopy and distance geometry.

Two-dimensional 1H-NMR studies have been performed on ribonuclease F1 (RNase F1), which contains 106 amino acid residues. Sequence-specific resonance assignments were accomplished for the backbone protons of 99 amino acid residues and for most of their side-chain protons. The three-dimensional structures were constructed on the basis of 820 interproton-distance restraints derived from NOE, 64 distance restraints for 32 hydrogen bonds and 33 phi torsion-angle restraints. A total of 40 structures were obtained by distance geometry and simulated-annealing calculations. The average root-mean-square deviation (residues 1-106) between the 40 converged structures and the mean structure obtained by averaging their coordinates was 0.116 +/- 0.018 nm for the backbone atoms and 0.182 +/- 0.015 nm for all atoms including the hydrogen atoms. RNase F1 was determined to be an alpha/beta-type protein. A well-defined structure constitutes the core region, which consists of a small N-terminal beta-sheet (beta 1, beta 2) and a central five-stranded beta-sheet (beta 3-beta 7) packed on a long helix. The structure of RNase F1 has been compared with that of RNase T1, which was determined by X-ray crystallography. Both belong to the same family of microbial ribonucleases. The polypeptide backbone fold of RNase F1 is basically identical to that of RNase T1. The conformation-dependent chemical shifts of the C alpha protons are well conserved between RNase F1 and RNase T1. The residues implicated in catalysis are all located on the central beta-sheet in a geometry similar to that of RNase T1.     

A stabilizing alpha/beta-hydrophobic core greatly contributes to hyperthermostability of archaeal [P62A]Ssh10b

The hyperthermophilic Ssh10b from Sulfolobus shibatae is a member of the Sac10b family, which has been postulated to play a role in chromosomal organization in Archaea. Ssh10b is capable of significantly constraining negative DNA supercoils at elevated temperatures. In this study, the solution structure of the dimeric P62A mutant Ssh10b ([P62A]Ssh10b) was determined by multidimensional NMR spectroscopy. The backbone 15N dynamics, H/D exchange with and without the denaturant GdmSCN, and chemical and thermal denaturation experiments were performed to investigate the molecular basis of high thermostability of [P62A]Ssh10b. Data analysis has revealed an alpha/beta-hydrophobic core consisting of two alpha-helices and one beta-sheet which are stabilized by cooperative hydrophobic and hydrogen-bonding interactions. This stabilizing alpha/beta-hydrophobic core of [P62A]Ssh10b exhibiting highly restricted internal motions is composed of residues having highly protected amide protons which exchange with solvent mostly by means of a global unfolding process. The K40N mutation greatly destabilizes the mutant [P62A]Ssh10b because this mutation disturbs the packing of alpha-helix against the beta-sheet reducing the stability of the alpha/beta-hydrophobic core in the mutant protein. In comparison with homologous mesophilic and thermophilic proteins, it can be presumed that the stabilizing alpha/beta-hydrophobic core in the [P62A]Ssh10b structure greatly contributes to the high thermostability of the protein.     

Investigation of the structural stability of cardiotoxin analogue III from the Taiwan cobra by hydrogen-deuterium exchange kinetics.

The conformational stability of a small ( approximately 7 kDa), all beta-sheet protein, cardiotoxin analogue III (CTX III), from the venom of the Taiwan cobra has been investigated by hydrogen-deuterium (H/D) exchange using two-dimensional NMR spectroscopy. The H/D exchange kinetics of backbone amide protons in CTX III has been monitored at pD 3.6 and 6.6 (at 25 degrees C), for over 5000 h. Examination of H/D exchange kinetics in the protein showed that a number of slowly exchanging residues are in the hydrophobic core of the protein. The average protection factor of the amide protons of residues belonging to the triple-stranded beta-sheet domain is about 20 times greater than that of those in the double-stranded beta-sheet segment. The residues in the C-terminal tail of the molecule, though structureless, have been found to exhibit significant protection against H/D exchange. Comparison of the quenched-flow H/D exchange data on CTX III with those obtained in the present study reveals that the most slowly exchanging portion constitutes the folding core of the protein.     

1H NMR resonance assignments, secondary structure, and global fold of Apo bovine calbindin D9k

The solution structure and dynamics of apo bovine calbindin D9k have been studied by a wide range of two-dimensional 1H nuclear magnetic resonance experiments. Due to the presence of conformational heterogeneity in the wild-type protein, the sequential resonance assignment was carried out on a Pro43----Gly mutant. By use of a combination of scalar correlation experiments acquired from H2O solution, 61 of the 76 1H spin systems could be assigned to particular amino acid types. The remaining resonances were assigned by a parallel series of experiments acquired from 2H2O solution. These spin system assignments provided a basis for complete sequential resonance assignments from interresidue backbone nuclear Overhauser effects (NOEs). Elements of secondary structure were identified from sequential and medium-range NOEs, backbone spin-spin coupling constants, and slowly exchanging amide protons. Four sections of helix are delineated, together with a short antiparallel beta-sheet interaction between the peptide loops involved in Ca2+ binding. The global fold is provided by combining these elements of secondary structure with a subset of the long-range, interhelix NOEs. Comparison with similar studies on the Ca2(+)-saturated protein indicates that at this crude level the structures are very similar. However, removal of the Ca2+ does dramatically affect the dynamics of the protein, as judged by amide proton exchange rates and aromatic ring rotation. This is particularly evident in the increased flexibility of the residues in the hydrophobic core.     

Investigation of the four cooperative unfolding units existing in the MICAL-1 CH domain

The solution structure of human MICAL-1 calpolnin homology (CH) domain is composed of six alpha helices and one 3(10) helix. To study the unfolding of this domain, we carry out native-state hydrogen exchange, intrinsic fluorescence and far-UV circular dichroism experiments. The free energy of unfolding, DeltaG(H2O), is calculated to be 7.11+/-0.58 kcal mol(-1) from GuHCl denaturation at pH 6.5. Four cooperative unfolding units are found using native-state hydrogen exchange experiment. Forty-seven slow-exchange residues can be studied by native-state hydrogen exchange experiments. From the concentration dependence of exchange rates, free energy of amide hydrogen with solvent, DeltaG(HX) and m-value (sensitivity of exposure to denaturant) are obtained, which reveal four cooperative unfolding units. The slowest exchanging protons are distributed throughout the whole hydrophobic core of the protein, which might be the folding core. These results will help us understand the structure of MICAL-1 CH domain more deeply.     

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

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

Exchange behavior of the H-bonded amide protons in the 3 to 13 helix of ribonuclease S

The preceding article shows that there are eight highly protected amide protons in the S-peptide moiety of RNAase S at pH 5, 0 degrees C. The residues with protected NH protons are 7 to 13, whose amide protons are H-bonded in the 3 to 13 alpha-helix, and Asp 14, whose NH proton is H-bonded to the CO group of Val47. We describe here the exchange behavior of these eight protected protons as a function of pH. Exchange rates of the individual NH protons are measured by 1H nuclear magnetic resonance in D2O. A procedure is used for specifically labeling with 1H only these eight NH protons. The resonance assignments of the eight protons are made chiefly by partial exchange, through correlating the resonance intensities in spectra taken when the peptide is bound and when it is dissociated from S-protein in 3.5 M-urea-d4, in D2O, pH 2.3, -4 degrees C. The two remaining assignments are made and some other assignments are checked by measurements of the nuclear Overhauser effect between adjacent NH protons of the alpha-helix. There is a transition in exchange behavior between pH 3, where the helix is weakly protected against exchange, and pH 5 where the helix is much more stable. At pH 3.1, 20 degrees C, exchange rates are uniform within the helix within a factor of two, after correction for different intrinsic exchange rates. The degree of protection within the helix is only 10 to 20-fold at this pH. At pH 5.1, 20 degrees C, the helix is more stable by two orders of magnitude and exchange occurs preferentially from the N-terminal end. At both pH values the NH proton of Asp 14, which is just outside the helix, is less protected by an order of magnitude than the adjacent NH protons inside the helix. Opening of the helix can be observed below pH 3.7 by changes in chemical shifts of the NH protons in the helix. At pH 2.4 the changes are 25% of those expected for complete opening. Helix opening is a fast reaction on the n.m.r. time scale (tau much less than 1 ms) unlike the generalized unfolding of RNAase S which is a slow reaction. Dissociation of S-peptide from S-protein in native RNAase S at pH 3.0 also is a slow reaction. Opening of the helix below pH 3.7 is a two-state reaction, as judged by comparing chemical shifts with exchange rates. The exchange rates at pH 3.1 are predicted correctly from the changes in chemical shift by assuming that helix opening is a two-state reaction. At pH values above 3.7, the nature of the helix opening reaction changes. These results indicate that at least one partially unfolded state of RNAase S is populated in the low pH unfolding transition.     

Microscopic stability of cold shock protein A examined by NMR native state hydrogen exchange as a function of urea and trimethylamine N-oxide

Native state hydrogen exchange of cold shock protein A (CspA) has been characterized as a function of the denaturant urea and of the stabilizing agent trimethylamine N-oxide (TMAO). The structure of CspA has five strands of beta-sheet. Strands beta1-beta4 have strongly protected amide protons that, based on experiments as a function of urea, exchange through a simple all-or-none global unfolding mechanism. By contrast, the protection of amide protons from strand beta5 is too weak to measure in water. Strand beta5 is hydrogen bonded to strands beta3 and beta4, both of which afford strong protection from solvent exchange. Gaussian network model (GNM) simulations, which assume that the degree of protection depends on tertiary contact density in the native structure, accurately predict the strong protection observed in strands beta1-beta4 but fail to account for the weak protection in strand beta5. The most conspicuous feature of strand beta5 is its low sequence hydrophobicity. In the presence of TMAO, there is an increase in the protection of strands beta1-beta4, and protection extends to amide protons in more hydrophilic segments of the protein, including strand beta5 and the loops connecting the beta-strands. TMAO stabilizes proteins by raising the free energy of the denatured state, due to highly unfavorable interactions between TMAO and the exposed peptide backbone. As such, the stabilizing effects of TMAO are expected to be relatively independent of sequence hydrophobicity. The present results suggest that the magnitude of solvent exchange protection depends more on solvent accessibility in the ensemble of exchange susceptible conformations than on the strength of hydrogen-bonding interactions in the native structure.     

Influence of pH on NMR structure and stability of the human prion protein globular domain

The NMR structure of the globular domain of the human prion protein (hPrP) with residues 121-230 at pH 7.0 shows the same global fold as the previously published structure determined at pH 4.5. It contains three alpha-helices, comprising residues 144-156, 174-194, and 200-228, and a short anti-parallel beta-sheet, comprising residues 128-131 and 161-164. There are slight, strictly localized, conformational changes at neutral pH when compared with acidic solution conditions: helix alpha1 is elongated at the C-terminal end with residues 153-156 forming a 310-helix, and the population of helical structure in the C-terminal two turns of helix alpha 2 is increased. The protonation of His155 and His187 presumably contributes to these structural changes. Thermal unfolding monitored by far UV CD indicates that hPrP-(121-230) is significantly more stable at neutral pH. Measurements of amide proton protection factors map local differences in protein stability within residues 154-157 at the C-terminal end of helix alpha 1 and residues 161-164 of beta-strand 2. These two segments appear to form a separate domain that at acidic pH has a larger tendency to unfold than the overall protein structure. This domain could provide a "starting point" for pH-induced unfolding and thus may be implicated in endosomic PrPC to PrPSc conformational transition resulting in transmissible spongiform encephalopathies.