NMR investigation of the dynamics of tryptophan side-chains in hemoglobins

NMR relaxation measurements of 15N spin-lattice relaxation rate (R(1)), spin-spin relaxation rate (R(2)), and heteronuclear nuclear Overhauser effect (NOE) have been carried out at 11.7T and 14.1T as a function of temperature for the side-chains of the tryptophan residues of 15N-labeled and/or (2H,15N)-labeled recombinant human normal adult hemoglobin (Hb A) and three recombinant mutant hemoglobins, rHb Kempsey (betaD99N), rHb (alphaY42D/betaD99N), and rHb (alphaV96W), in the carbonmonoxy and the deoxy forms as well as in the presence and in the absence of an allosteric effector, inositol hexaphosphate (IHP). There are three Trp residues (alpha14, beta15, and beta37) in Hb A for each alphabeta dimer. These Trp residues are located in important regions of the Hb molecule, i.e. alpha14Trp and beta15Trp are located in the alpha(1)beta(1) subunit interface and beta37Trp is located in the alpha(1)beta(2) subunit interface. The relaxation experiments show that amino acid substitutions in the alpha(1)beta(2) subunit interface can alter the dynamics of beta37Trp. The transverse relaxation rate (R(2)) for beta37Trp can serve as a marker for the dynamics of the alpha(1)beta(2) subunit interface. The relaxation parameters of deoxy-rHb Kemspey (betaD99N), which is a naturally occurring abnormal human hemoglobin with high oxygen affinity and very low cooperativity, are quite different from those of deoxy-Hb A, even in the presence of IHP. The relaxation parameters for rHb (alphaY42D/betaD99N), which is a compensatory mutant of rHb Kempsey, are more similar to those of Hb A. In addition, TROSY-CPMG experiments have been used to investigate conformational exchange in the Trp residues of Hb A and the three mutant rHbs. Experimental results indicate that the side-chain of beta37Trp is involved in a relatively slow conformational exchange on the micro- to millisecond time-scale under certain experimental conditions. The present results provide new dynamic insights into the structure-function relationship in hemoglobin.     

Backbone dynamics and energetics of a calmodulin domain mutant exchanging between closed and open conformations.

Previous studies have suggested that the Ca2+-saturated E140Q mutant of the C-terminal domain of calmodulin exhibits equilibrium exchange between "open" and "closed" conformations similar to those of the Ca2+-free and Ca2+-saturated states of wild-type calmodulin. The backbone dynamics of this mutant were studied using15N spin relaxation experiments at three different temperatures. Measurements at each temperature of the15N rate constants for longitudinal and transverse auto-relaxation, longitudinal and transverse cross-correlation relaxation, and the1H-15N cross-relaxation afforded unequivocal identification of conformational exchange processes on microsecond to millisecond time-scales, and characterization of fast fluctuations on picosecond to nanosecond time-scales using model-free approaches. The results show that essentially all residues of the protein are involved in conformational exchange. Generalized order parameters of the fast internal motions indicate that the conformational substates are well folded, and exclude the possibility that the exchange involves a significant population of unfolded or disordered species. The temperature dependence of the order parameters offers qualitative estimates of the contribution to the heat capacity from fast fluctuations of the protein backbone, revealing significant variation between the well-ordered secondary structure elements and the more flexible regions. The temperature dependence of the conformational exchange contributions to the transverse auto-relaxation rate constants directly demonstrates that the microscopic exchange rate constants are greater than 2.7x10(3)s-1at 291 K. The conformational exchange contributions correlate with the chemical shift differences between the Ca2+-free and Ca2+-saturated states of the wild-type protein, thereby substantiating that the conformational substates are similar to the open and closed states of wild-type calmodulin. Taking the wild-type chemical shifts to represent the conformational substates of the mutant and populations estimated previously, the microscopic exchange rate constants could be estimated as 2x10(4)to 3x10(4)s-1at 291 K for a subset of residues. The temperature depen dence of the exchange allows the characterization of apparent energy barriers of the conformational transition, with results suggesting a complex process that does not correspond to a single global transition between substates.     

Measurement of 15N relaxation in deuterated amide groups in proteins using direct nitrogen detection

15N chemical shielding tensors contain useful structural information, and their knowledge is essential for accurate analysis of protein backbone dynamics. The anisotropic component (CSA) of 15N chemical shielding can be obtained from 15N relaxation measurements in solution. However, the predominant contribution to nitrogen relaxation from 15N-(1)H dipolar coupling in amide groups limits the sensitivity of these measurements to the actual CSA values. Here we present nitrogen-detected NMR experiments for measuring 15N relaxation in deuterated amide groups in proteins, where the dipolar contribution to 15N relaxation is significantly reduced by the deuteration. Under these conditions nitrogen spin relaxation becomes a sensitive probe for variations in 15N chemical shielding tensors. Using the nitrogen direct-detection experiments we measured the rates of longitudinal and transverse 15N relaxation for backbone amides in protein G in D(2)O at 11.7 T. The measured relaxation rates are validated by comparing the overall rotational diffusion tensor obtained from these data with that from the conventional 15N relaxation measurements in H(2)O. This analysis revealed a 17-24 degree angle between the NH-bond and the unique axis of the 15N chemical shielding tensor.     

Improved sensitivity in indirect monitoring of chemical shifts of proton-heteronuclear spin pairs (1H-13C and 1H-15N) in 3D and 4D NMR spectroscopy

In three-dimensional and four-dimensional experiments on doubly labelled proteins not only heteronuclear (¹³C or ¹⁵N) but also proton (¹H) frequencies are often indirectly monitored, rather than being directly observed. In this communication we show how in these experiments by overlaying ¹H and heteronuclear evolutions one can obtain decreased apparent relaxation rates of ¹H signals, yielding improved sensitivity. The new method applies to spin pairs like ¹H-¹⁵N, as in amide groups, or ¹H-¹³C, as in methine groups of alpha or aromatic systems.     

Local mobility of 15N labeled biomolecules characterized through cross-correlation rates: Applications to paramagnetic proteins*

The mobility of 15N labeled proteins can be characterized by measuring the cross-correlation rates N,NI that govern the conversion of Zeeman order Nz of an amide 15N nucleus into longitudinal two-spin order 2NzIz involving the amide 15N and 1H nuclei. This represents an alternative to the measurement of 15N self-relaxation rates 1/T1 and 1/T2 or 1/T1. The rate of interconversion between Nz and 2NzIz is due to cross-correlation between fluctuations of different interactions and is not affected by a variety of relaxation mechanisms that contribute to the self-relaxation rates 1/T1, 1/T2 and 1/T1. Spin diffusion among protons, which affects the measurements, can be quenched by various means that are evaluated by experiments and simulations. By applying an off-resonance radio-frequency (RF) field in the vicinity of the nitrogen resonance, the spectral density function J() can be determined at the frequency origin and at the nitrogen Larmor frequency. The methods are applied to the paramagnetic High-Potential Iron-Sulfur Protein iso I (HiPIP I) from E. halophila in its reduced state.     

Protein global fold determination using site-directed spin and isotope labeling

We describe a simple experimental approach for the rapid determination of protein global folds. This strategy utilizes site-directed spin labeling (SDSL) in combination with isotope enrichment to determine long-range distance restraints between amide protons and the unpaired electron of a nitroxide spin label using the paramagnetic effect on relaxation rates. The precision and accuracy of calculating a protein global fold from only paramagnetic effects have been demonstrated on barnase, a well-characterized protein. Two monocysteine derivatives of barnase, (H102C) and (H102A/Q15C), were 15N enriched, and the paramagnetic nitroxide spin label, MTSSL, attached to the single Cys residue of each. Measurement of amide 1H longitudinal relaxation times, in both the oxidized and reduced states, allowed the determination of the paramagnetic contribution to the relaxation processes. Correlation times were obtained from the frequency dependence of these relaxation processes at 800, 600, and 500 MHz. Distances in the range of 8 to 35 A were calculated from the magnitude of the paramagnetic contribution to the relaxation processes and individual amide 1H correlation times. Distance restraints from the nitroxide spin to amide protons were used as restraints in structure calculations. Using nitroxide to amide 1H distances as long-range restraints and known secondary structure restraints, barnase global folds were calculated having backbone RMSDs <3 A from the crystal structure. This approach makes it possible to rapidly obtain the overall topology of a protein using a limited number of paramagnetic distance restraints.     

An evaluation of detergents for NMR structural studies of membrane proteins

Structural information on membrane proteins lags far behind that on soluble proteins, in large part due to difficulties producing homogeneous, stable, structurally relevant samples in a membrane-like environment. In this study 25 membrane mimetics were screened using 2D (1)H-(15)N heteronuclear single quantum correlation NMR experiments to establish sample homogeneity and predict fitness for structure determination. A single detergent, 1-palmitoyl-2-hydroxy-sn-glycero-3-[phospho-RAC-(1-glycerol)] (LPPG), yielded high quality NMR spectra with sample lifetimes greater than one month for the five proteins tested - R. sphaeroides LH1 alpha and beta subunits, E. coli and B. pseudofirmus OF4 ATP synthase c subunits, and S. aureus small multidrug resistance transporter - with 1, 2, or 4 membrane spanning alpha-helices, respectively. Site-specific spin labeling established interhelical distances in the drug transporter and genetically fused dimers of c subunits in LPPG consistent with in vivo distances. Optical spectroscopy showed that LH1 beta subunits form native-like complexes with bacteriochlorophyll a in LPPG. All the protein/micelle complexes were estimated to exceed 100 kDaltons by translational diffusion measurements. However, analysis of (15)N transverse, longitudinal and (15)N[(1)H] nuclear Overhauser effect relaxation measurements yielded overall rotational correlation times of 8 to 12 nsec, similar to a 15-20 kDalton protein tumbling isotropically in solution, and consistent with the high quality NMR data observed.     

Signal assignments and chemical-shift structural analysis of uniformly 13C, 15N-labeled peptide, mastoparan-X, by multidimensional solid-state NMR under magic-angle spinning

Carbon-13 and nitrogen-15 signals of fully isotope-labeled 15-residue peptide, glycinated mastoparan-X, in a solid state were assigned by two- and three-dimensional NMR experiments under magic-angle spinning conditions. Intra-residue spin connectivities were obtained with multidimensional correlation experiments for C'-C(alpha)-C(beta) and N-C(alpha)-C(beta). Sequence specific assignments were performed with inter-residue C(alpha)-C(alpha) and N-C(alpha)C(beta) correlation experiments. Pulse sequences for these experiments have mixing periods under recoupled zero- and double-quantum (13)C-(13)C and (15)N-(13)C dipolar interactions. These correlation spectra allowed the complete assignments of (13)C and (15)N backbone and (13)C(beta) signals. Chemical shift analysis of the (13)C and (15)N signals based on empirical and quantum chemical databases for proteins indicated that the backbone between residues 3 and 14 forms alpha-helix and residue 2 has extended conformation in the solid state. This structure was compared with the G-protein- and membrane-bound structures of mastoparan-X.     

Probing Slow Backbone Dynamics in Proteins Using TROSY-Based Experiments to Detect Cross-Correlated Time-Modulation of Isotropic Chemical Shifts

The difference in the relaxation rates of zero-quantum (ZQ) and double-quantum (DQ) coherences is the result of three principal mechanisms. These include the cross-correlation between the chemical shift anisotropies of the two participating nuclei, dipolar interactions with remote protons as well as interference effects due to the time-modulation of their isotropic chemical shifts as a consequence of slow micros-ms dynamics. The last effect when present, dominates the others resulting in large differences between the relaxation rates of ZQ and DQ coherences. We present here four sets of TROSY-based (Salzmann et al., 1998) experiments that measure this effect for several pairs of backbone nuclei including (15)N, (13)C(alpha) and (13)C'. These experiments allow the detection of the presence of slow dynamic processes in the protein backbone including correlated motion over two and three bonds. Further, we define a new parameter chi which represents the extent of correlated motion on the slow (micros-ms) timescale. This methodology has been applied to (15)N,(13)C,REDPRO-(2)H-labeled (Shekhtman et al., 2002) human ubiquitin. The ubiquitin backbone is seen to exhibit extensive dynamics on the slow timescale. This is most pronounced in several residues in the N-terminal region of the alpha-helix and in the loop connecting the strands beta(4) and beta(5). These residues which include Glu24, Asn25, Glu51 and Asp 52 form a continuous surface. As an additional benefit, the measured rates confirm the dependence of the (13)C(alpha) chemical shift tensor on local secondary structure of the protein backbone.     

Enhancement of the steady-state magnetization in TROSY experiments

Under the condition that the longitudinal relaxation time of spin I is shorter than the longitudinal relaxation time of spin S the steady-state magnetization in [S,I]-TROSY-type experiments can be enhanced by intermediate storage of a part of the steady-state magnetization of spin I on spin S with a pulse sequence element during the relaxation delay. It is demonstrated with samples ranging in size from the 1 kDa cyclosporin to the 110 kDa ¹⁵N,²H-labeled dihydroneopterin Aldolase that intermediate storage of steady-state magnetization in a [¹⁵N,¹H]-TROSY experiment yields a signal gain of 10–25%. The method proposed here for intermediate storage of steady-state magnetization can be implemented in any [¹⁵N,¹H]-TROSY-type experiments.     

Direct measurement of the 15N CSA/dipolar relaxation interference from coupled HSQC spectra

Here we propose a method for the measurement of the ¹⁵N CSA/dipolar relaxation interference based on direct comparison of the ¹⁵N doublet components observed in a ¹H-coupled ¹H-¹⁵N HSQC-type spectrum. This allows the determination of the cross-correlation rates with no need for correction factors associated with other methods. The signal overlap problem of coupled HSQC spectra is addressed here by using the IPAP scheme (Ottiger et al., 1998). The approach is applied to the B3 domain of protein G to show that the method provides accurate measurements of the ¹⁵N CSA/dipolar cross-correlation rates.     

Mapping of the spectral densities of N-H bond motions in eglin c using heteronuclear relaxation experiments.

A new strategy is used for studying the internal motions of proteins based on measurements of NMR relaxation parameters. The strategy yields values of the so-called spectral density functions J(omega) for N-H bond vectors. The spectral density functions are related to the distribution of frequencies contained in the rotational (overall and internal) motions of these NH bond vectors. No a priori model assumptions about the dynamics are required in this approach. The method involves measurements of six relaxation parameters consisting of 15N longitudinal relaxation rates, transverse relaxation rates of in-phase and antiphase coherence, the relaxation rates of heteronuclear 1H-15N two-spin order, the heteronuclear 1H-15N nuclear Overhauser effects, and longitudinal relaxation rates of the amide protons. The values of the spectral density functions at the five frequencies 0, omega N, omega H + omega N, omega H, and omega H - omega N are determined from the relaxation parameters using analytical relations derived previously [Peng & Wagner (1992) J. Magn. Reson. 98, 308-332]. Here, the method is applied to characterize the backbone dynamics of the 15N-enriched proteinase inhibitor eglin c, a protein of 70 residues. The values for J(0) and J(omega N = 50 MHz) vary significantly with the amino acid sequence, whereas the spectral densities at higher frequencies, J(450 MHz), J(500 MHz), and J(550 MHz), are typically much smaller and show no significant variation with the sequence. The collective behavior of the J(omega) values indicate greater internal motion for the proteinase binding loop residues and the first eight N-terminal residues. The additional internal motion in these regions is in the rate range below 450 MHz. The values of J(omega) are also compared with root mean square deviations (rmsds) of backbone atoms as obtained in NMR structure determinations. Low values of J(0) and J(omega N) are correlated with high rmsds. Spectral densities at higher frequencies, J(450 MHz), J(500 MHz), and J(550 MHz), are small and show no correlation with rmsds. A comparison with the spectral density functions obtained by fitting the experimental data to the functional dependence of the Lipari and Szabo formalism [Lipari & Szabo (1982a) J. Am. Chem. Soc. 104, 4546-4559] is made.     

An NMR and circular dichroism study of the interaction of thiocyanate with human and cross-linked hemoglobin: identification of Lys-alpha-99 as a possible dissociation linked binding site

The interaction of thiocyanate with human native and cross-linked oxyhemoglobin (oxyHb), and methemoglobin (metHb) has been investigated by optical spectroscopy, circular dichroism (CD) and nuclear spin lattice relaxation rate measurements. The interaction of thiocyanate anion with human hemoglobin has been investigated by NMR measurements of the nuclear spin lattice relaxation rate of N(15) labeled thiocyanate in the presence of cyanomethemoglobin and cross-linked cyanomethemoglobin. Results show that thiocyanate is located approximately 8.9 and 6.2 A away from the heme group in cyanomethemoglobin and cross-linked cyanomethemoblobin, respectively. These results are consistent with the binding of SCN(-) at the lys-alpha-99 in the unmodified hemoglobin. Since this site is blocked in the cross-linked hemoglobin, the binding site is different. Results show that one mole of SCN(-) is binding to one mole of oxyhemoglobin suggesting that binding at the lys-alpha-99 is linked to dissociation of the hemoglobin tetramer into dimers due to its location at the alpha(1)beta(2) interface. Circular dichroism studies show that the interaction of thiocyanate with oxyHb decreases the optical rotation at 240 nm indicating a conformational change of the protein, which influences the electronic transitions of a number of peptide bonds or (and) a few aromatic side chains.     

Electron-electron spin-spin interaction in spin-labeled low-spin methemoglobin.

Nitroxyl free radical electron spin relaxation times for spin-labeled low-spin methemoglobins were measured between 6 and 120 K by two-pulse electron spin echo spectroscopy and by saturation recovery electron paramagnetic resonance (EPR). Spin-lattice relaxation times for cyano-methemoglobin and imidazole-methemoglobin were measured between 8 and 25 K by saturation recovery and between 4.2 and 20 K by electron spin echo. At low temperature the iron electron spin relaxation rates are slow relative to the iron-nitroxyl electron-electron spin-spin splitting. As temperature is increased, the relaxation rates for the Fe(III) become comparable to and then greater than the spin-spin splitting, which collapses the splitting in the continuous wave EPR spectra and causes an increase and then a decrease in the nitroxyl electron spin echo decay rate. Throughout the temperature range examined, interaction with the Fe(III) increases the spin lattice relaxation rate (1/T1) for the nitroxyl. The measured relaxation times for the Fe(III) were used to analyze the temperature-dependent changes in the spin echo decays and in the saturation recovery (T1) data for the interacting nitroxyl and to determine the interspin distance, r. The values of r for three spin-labeled methemoglobins were between 15 and 15.5 A, with good agreement between values obtained by electron spin echo and saturation recovery. Analysis of the nitroxyl spin echo and saturation recovery data also provides values of the iron relaxation rates at temperatures where the iron relaxation rates are too fast to measure directly by saturation recovery or electron spin echo spectroscopy. These results demonstrate the power of using time-domain EPR measurements to probe the distance between a slowly relaxing spin and a relatively rapidly relaxing metal in a protein.     

Cross correlations between 13C-1H dipolar interactions and 15N chemical shift anisotropy in nucleic acids

Two sets of cross-correlated relaxation rates involving chemical shift anisotropy and dipolar interactions have been measured in an RNA kissing complex. In one case, both the CSA and dipolar interaction tensors are located on the same nucleotide base and are rigidly fixed with respect to each other. In the other case, the CSA tensor is located on the nucleotide base whereas the dipolar interaction is located on the adjoining ribose unit. Analysis of the measured rates in terms of isotropic or anisotropic rotational diffusion has been carried out for both cases. A marked difference between the two models is observed for the cross-correlation rates involving rigidly fixed spin interactions. The influence of internal motions about the glycosidic linkage between the nucleotide base and the ribose unit on cross-correlated relaxation rates has been estimated by applying a model of restricted rotational diffusion. Local motions seem to have a more pronounced effect on cross-correlated relaxation rates when the two spin interactions are not rigidly fixed with respect to each other.     

Quantitative estimation of magnitude and orientation of the CSA tensor from field dependence of longitudinal NMR relaxation rates

A method is presented that makes it possible to estimate both the orientation and the magnitude of the chemical shift anisotropy (CSA) tensor in molecules with a pair of spin 1/2 nuclei, typically ¹³C-¹H or ¹⁵ N-¹H. The method relies on the fact that the longitudinal cross-correlation rate as well as a linear combination of the autorelaxation rates of longitudinal heterospin magnetization, longitudinal two-spin order and longitudinal proton magnetization are proportional to the spectral density at the Larmor frequency of the heterospin. Therefore the ratio between the cross-correlation rate and the above linear combination is independent of the dynamics. From the field dependence of the ratio both the magnitude and the orientation of the CSA tensor can be estimated. The method is applicable to molecules in all motional regimes and is not limited to molecules in extreme narrowing or slow tumbling, nor is it sensitive to chemical exchange broadening. It is tested on the 22 amino acid residue peptide motilin, selectively ¹³ C labeled in the ortho positions in the ring of the single tyrosine residue. In the approximation of an axially symmetric ¹³C CSA tensor, the symmetry axis of the CSA tensor makes an angle of 23° ± 1° to the ¹³ C-¹H bond vector, and has a magnitude of 156 ± 5 ppm. This is in close agreement with solid-state NMR data on tyrosine powder [Frydman et al. (1992) Isr. J. Chem., 32, 161–164].     

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.