Chapter 3. Problems associated with measurement of coupling constants in proteins

Traditional one-dimensional 1H NMR spectroscopy is sufficient to provide information on chemical shifts and coupling constants in small organic molecules. Owing to the long transverse relaxation times, proton resonances have very narrow line widths in comparison with the inhomogeneity of the magnetic field. Hence, relatively small coupling constants can often be measured with high accuracy from a well-dispersed spectrum.

In the case of biological macromolecules, the situation is very different for two reasons: increased spectral overlap and faster transverse relaxation of protons. The number of proton resonances in the same rather limited chemical shift range increases rapidly as the protein size increases. For example, a 10-kDa protein has approximately 100 amide proton resonances in the chemical shift range of 3 to 4 ppm. In the case of a 30-kDa protein, the number of resonances is on the order of 300. In addition, as the isotropic molecular rotational correlation time, τc, increases with increasing molecular weight, the line width is dictated by the shorter transverse relaxation time of the corresponding spin instead of magnetic field inhomogeneity. Ultimately, the line width is larger than the coupling constant of interest. These two dilemmas make it impossible to obtain structurally important vicinal proton-proton coupling constants simply from a 1D proton spectrum.

An obvious solution to the overlap problem is to increase the dimensionality of the spectrum. For example, a two-dimensional phase-sensitive COSY experiment (Marion et al. 1983) has been successfully used with small proteins. The dispersion of the proton signals in the COSY spectrum is to the power of two better than in the 1D proton spectrum. However, applicability of the phase-sensitive COSY is limited to small proteins, because in the case of determination of 3JHNHα, the cross-peak intensity between 1HN and 1Hα is low due to self-cancellation of antiphase splitting in the presence of broad lines. Furthermore, the magnitude of the coupling constant is not accurate because the antiphase splitting greatly overestimates the true coupling constant (Neuhaus et al. 1985). For this reason, several line-fitting techniques have been developed for the extraction of the true coupling constant in the phase-sensitive COSY experiment (Neuhaus et al. 1985).

One major limitation of homonuclear experiments stems from the rapid transverse relaxation of protons. This reduces their usability for protein structure determination due to coherence transfer inefficiency in multi-dimensional experiments. Proton relaxation is dominated by the large dipole-dipole (DD) interaction between protons in the protein main-chain and side-chains. Consequently, the proton line width increases rapidly with increasing rotational correlation time. As can be seen from Figure 1a, the proton line width is close to 10 Hz when the correlation time is of the order 10 ns. In contrast, the line width of the backbone amide nitrogen is less than 5 Hz (Figure 1b). It is then obvious that the development of heteronuclear correlation experiments, which allow a greater dispersion of NMR signals with higher coherence transfer efficiency, is indispensable for larger proteins.

Figure 1. Protein resonance line widths as a function of isotropic molecular rotational correlation time τc. (a) Plots are shown for (--) 1H spins, (…) 1H spins directly bonded to 13C, and (- - -) 1H spins directly bonded to 15N spin. (b) Line widths for 13C and 15N spins. Plots are shown for (---) 13C, (- - -) 13C1Hz, (- . - . - . -) 15N, and (…) 15N1Hz coherences (Cavanagh et al. 1994).

Introduction of isotopic labeling in protein NMR spectroscopy in the late 1980’s improved the sensitivity and resolution of protein NMR spectra. Thus, magnetically inactive 12C, and spin-1 14N nuclei are replaced by active spin- nuclei 13C and 15N. In addition, the isotope enrichment facilitated the measurement of several, also heteronuclear, coupling constants in the protein main-chain. More importantly, new techniques were devised to determine coupling constants in general. It is now possible to measure finite couplings, without the need to resolve them directly, by utilizing for instance, E.COSY (Griesinger et al. 1986), multiple-quantum coherence principles (Braunschweiler et al. 1983) and quantitative J-correlation (Blake et al. 1992). The 15N, 13C double labeling of protein samples has increased the size limit of NMR-feasible proteins from 10 kDa to the 20-30 kDa limit. Table 1 represents typical transverse relaxation times found for 1H, 15N, and 13C nuclei in a protein backbone at a 30-kDa size regime.

Table 1. Transverse relaxation times for different spins in the 29 kDa, uniformly 15N/13C and 15N/13C/2H labeled, HCA II protein at 600 1H frequency. (Farmer & Venters 1999).

NucleusaT2,exp (ms)bT2,calc (ms)
13Ca124.018.0
13C’47.045.0
15N52.049.0
15N(Hz)43.0-
1HN (α-helix)24.021.0
1HN (β-sheet)29.017.0
aMeasured T2 values for perdeuterated HCA II, bCalculated T2 values for protonated HCA II

To increase the protein size further, rapid transverse relaxation of heteronuclei must be taken into account as well. Transverse relaxation times for aliphatic carbons with proteins at the 20-30 kDa regime have been shown to be as short as 10-20 ms (Gardner & Kay 1998) (Table 1). It is inevitable that this seriously limits sensitivity and resolution in the experiments that record chemical shift or relay magnetization via aliphatic carbons. A very efficient solution to this problem is perdeuteration, that is, replacement of non-exchangeable protons with deuterium. This decreases the carbon relaxation rate over 10-fold, since the dipole-dipole relaxation between 13C and its directly bound 2H is much smaller than that between 13C and 1H, owing to the 6.5 times smaller magnetogyric ratio of 2H. Perdeuteration is beneficial for 15N, 1H correlation experiments as well, since the relaxation rate of the amide proton decreases approximately twofold, due to the disappearance of homonuclear dipolar relaxation caused by aliphatic protons. It has been shown that 40% of the transverse relaxation rate of the amide proton arises from the dipolar contribution with the aliphatic protons (Markus et al. 1994). Amide nitrogens and carbonyl carbons are, however, much less affected by perdeuteration (Table 1).

Unfortunately, the improved sensitivity and resolution does not come without a cost. As mentioned earlier, NOE information available from non-exchangeable protons is reduced or completely lost (in the case of fully perdeuterated proteins), which results in low-resolution protein structures (Venters et al. 1995; Gardner & Kay 1998). Problems associated with a phenomenon referred to as spin diffusion are, however, less serious because of the decreased proton density (Gardner & Kay 1998; Farmer & Venters 1999). Thus, a separation between protons i and j can be determined more accurately thanks to the absence of magnetization relay via mutual spin k. Furthermore, it permits the use of longer NOE mixing times, allowing the measurement of larger distances than would be possible in protonated systems (Venters et al. 1995).

Chemical shift anisotropy (CSA) is a significant source of T2 relaxation for 15N and 13C’ spins in proteins. As the CSA interaction increases with increasing magnetic field, it reduces the gain in sensitivity and resolution obtained by the use of the highest magnetic fields. The second important source of relaxation is dipole-dipole (DD) interaction between spins. It has been commonly known for several years that the interference between chemical shift anisotropy (CSA) and dipole-dipole (DD) relaxation mechanisms results in differently relaxing 15N-1H multiplet components (Goldman 1984). Pervushin and co-workers recently showed that at proton frequencies near 1.1 GHz, almost complete cancellation of transverse relaxation effects within a 15N-1H moiety is expected for one of the four multiplet components in fully 15N-1H coupled 2D experiment. By using a certain spin-state edited experiment for the selection of the most slowly relaxing 15N-1H multiplet component (TROSY, transverse relaxation optimized spectroscopy), significant improvement in sensitivity and resolution in 15N-1H detected experiments, especially with larger perdeuterated protein samples, can be achieved (Pervushin et al. 1997). Interestingly, in addition to 1HN spin relaxation, use of perdeuterated samples with TROSY selection also dramatically decreases the relaxation rate of 15N spin. Remote protons have been shown to correspond with 75% of the residual T2 of 15N relaxation (Pervushin et al. 1997). Therefore, TROSY selection concomitantly with perdeuteration is most effective for large β-sheet proteins, where there are only a few routes to amide nitrogen relaxation. For α-helical proteins, the TROSY effect is somewhat less dramatic because of the proximity of 1HN(i-1) and 1HN(i+1) spins to the 15N(i) spin. In highly perdeuterated β-sheet proteins, it may also be advantageous to prevent the formation of NxCz antiphase coherence by coherent carbon decoupling, since longitudinal 13Cz spin relaxation forms a significant source of residual relaxation for 15N T2. Table 2 shows the average T2 relaxation times of the TROSY peak for the backbone 15N spin in the uniformly 15N/13C-labeled protein EIN, and the corresponding T2 (15N) for the 15N/13C/2H labeled sample (Kontaxis et al. 2000).

Table 2. Transverse relaxation times for the TROSY and anti-TROSY components of backbone 15N in 8.6 kDa 15N/13C and 15N/13C/2H labeled ubiquitin, and in 30 kDa 15N/13C and 15N/13C/2H labeled EIN at 600 and 800 MHz 1H frequency (Kontaxis et al. 2000).

ProteinaT2 (ms)bT2 (ms)cT2 (ms)

dT2 (ms)

EIN-15N/13C/2H (800 MHz)131.027.0

51.0

52.0
EIN-15N/13C/2H (600 MHz)118.033.0

60.0

57.0
EIN-15N (800 MHz)79.022.0

52.0

35.0
EIN-15N (600 MHz)72.028.0

57.0

39.0
UBI-15N/13C/2H (800 MHz)200.045.0

74.0

74.0
UBI-15N/13C/2H (600 MHz)185.055.0

93.0

93.0
UBI-15N/13C (800 MHz)111.040.0

73.0

59.0
UBI-15N/13C (600 MHz)104.046.0

87.0

67.0
adownfield TROSY component, bupfield anti-TROSY component, caverage time for α and β components in the refocused INEPT experiment, daverage time for α and β components in the usual HSQC experiment

According to Table 2, it can be seen that the transverse relaxation time of the 15N spin shortens with respect to increasing magnetic field strength due to increasing cross-correlation between CSA and DD in the usual 15N-HSQC experiment, where slow and fast relaxing spin-states are interchanged by a 180°(1H) pulse applied during t1. Consequently, the 15N line width increases with increasing magnetic field strength. On the other hand, line width of the TROSY component becomes narrower as a function of polarizing magnetic field strength. Interestingly, the transverse relaxation of 15N is largely affected by remote 1H spins in the regular HSQC experiment, resulting in a broader line for the 15N spin than for that in the corresponding refocused INEPT experiment. In contrast, the transverse relaxation of the 15N spin in perdeuterated proteins does not receive any major contribution from the 1H spin flips, yielding comparable line widths between the regular HSQC and the refocused INEPT experiments. It can also be seen that the 15N transverse relaxation time in perdeuterated proteins can be almost twice as long as that in protonated samples due to cancellation of remote, aliphatic proton interaction with the 15N nucleus. The so-called TROSY implementation (Pervushin et al. 1997; Andersson et al. 1998; Meissner et al. 1998; Salzmann et al. 1998; Yang & Kay 1999; Rance et al. 1999), together with deuterium labeling and utilization of directional information available through dipolar couplings, will increase the size limit of proteins applicable for structure determination much further than was predicted earlier.