High resolution <Superscript>13</Superscript>C-detected solid-state NMR spectroscopy of a deuterated protein

High resolution ¹³C-detected solid-state NMR spectra of the deuterated beta-1 immunoglobulin binding domain of the protein G (GB1) have been collected to show that all ¹⁵N, ¹³C′, ¹³Cα and ¹³Cβ sites are resolved in ¹³C–¹³C and ¹⁵N–¹³C spectra, with significant improvement in T ₂ relaxation times and resolution at high magnetic field (750 MHz). The comparison of echo T ₂ values between deuterated and protonated GB1 at various spinning rates and under different decoupling schemes indicates that ¹³Cα T ₂′ times increase by almost a factor of two upon deuteration at all spinning rates and under moderate decoupling strength, and thus the deuteration enables application of scalar-based correlation experiments that are challenging from the standpoint of transverse relaxation, with moderate proton decoupling. Additionally, deuteration in large proteins is a useful strategy to selectively detect polar residues that are often important for protein function and protein–protein interactions.     

Temperature dependence of fast carbonyl backbone dynamics in chicken villin headpiece subdomain

Temperature-dependence of protein dynamics can provide information on details of the free energy landscape by probing the characteristics of the potential responsible for the fluctuations. We have investigated the temperature-dependence of picosecond to nanosecond backbone dynamics at carbonyl carbon sites in chicken villin headpiece subdomain protein using a combination of three NMR relaxation rates: ¹³C′ longitudinal rate, and two cross-correlated rates involving dipolar and chemical shift anisotropy (CSA) relaxation mechanisms, ¹³C′/¹³C′-¹³C^α CSA/dipolar and ¹³C′/¹³C′–¹⁵N CSA/dipolar. Order parameters have been extracted using the Lipari-Szabo model-free approach assuming a separation of the time scales of internal and molecular motions in the 2–16°C temperature range. There is a gradual deviation from this assumption from lower to higher temperatures, such that above 16°C the separation of the time scales is inconsistent with the experimental data and, thus, the Lipari-Szabo formalism can not be applied. While there are variations among the residues, on the average the order parameters indicate a markedly steeper temperature dependence at backbone carbonyl carbons compared to that probed at amide nitrogens in an earlier study. This strongly advocates for probing sites other than amide nitrogen for accurate characterization of the potential and other thermodynamics characteristics of protein backbone.     

Backbone dynamics of bacteriorhodopsin as studied by <Superscript>13</Superscript>C solid-state NMR spectroscopy

The surface dynamics of bacteriorhodopsin was examined by measurements of site-specific ¹³C–¹H dipolar couplings in [3-¹³C]Ala-labeled bacteriorhodopsin. Motions of slow or intermediate frequency (correlation time <50 µs) scale down ¹³C–¹H dipolar couplings according to the motional amplitude. The two-dimensional dipolar and chemical shift (DIPSHIFT) correlation technique was utilized to obtain the dipolar coupling strength for each resolved peak in the ¹³C MAS solid-state NMR spectrum, providing the molecular order parameter of the respective site. In addition to the rotation of the Ala methyl group, which scales the dipolar coupling to 1/3 of the rigid limit value, fluctuations of the C–C vector result in additional motional averaging. Typical order parameters measured for mobile sites in bacteriorhodopsin are between 0.25 and 0.29. These can be assigned to Ala103 of the C–D loop and Ala235 at the C-terminal -helix protruded from the membrane surface, and Ala196 of the F–G loop, as well as to Ala228 and Ala233 of the C-terminal -helix and Ala51 from the transmembrane -helix. Such order parameters departing significantly from the value of 0.33 for rotating methyl groups are obviously direct evidence for the presence of fluctuation motions of the Ala C–C vectors of intact preparations of fully hydrated, wild-type bacteriorhodopsin at ambient temperature. The order parameter for Ala160 from the expectantly more flexible E–F loop, however, is unavailable under highest-field NMR conditions, probably because increased chemical shift anisotropy together with intrinsic fluctuation motions result in an unresolved ¹³C NMR signal.     

Cross-correlated spin relaxation effects in methyl <Superscript>1</Superscript>H CPMG-based relaxation dispersion experiments: Complications and a simple solution

Artifacts associated with the measurement of methyl ¹H single quantum CPMG-based relaxation dispersion profiles are described. These artifacts arise due to the combination of cross-correlated spin relaxation effects involving intra-methyl ¹H–¹H dipolar interactions and imperfections in ¹H refocusing pulses that are applied during CPMG intervals that quantify the effects of chemical exchange on measured transverse relaxation rates. As a result substantial errors in extracted exchange parameters can be obtained. A simple work-around is presented where the ¹H chemical shift difference between the exchanging states is extracted from a combination of ¹³C single quantum and ¹³C–¹H multiple quantum dispersion profiles. The approach is demonstrated with an application to a folding/unfolding reaction involving a G48M mutant Fyn SH3 domain.     

A suite of Mathematica notebooks for the analysis of protein main chain 15N NMR relaxation data

A suite of Mathematica notebooks has been designed to ease the analysis of protein main chain ¹⁵N NMR relaxation data collected at a single magnetic field strength. Individual notebooks were developed to perform the following tasks: nonlinear fitting of ¹⁵N-T ₁ and -T ₂ relaxation decays to a two parameter exponential decay, calculation of the principal components of the inertia tensor from protein structural coordinates, nonlinear optimization of the principal components and orientation of the axially symmetric rotational diffusion tensor, model-free analysis of ¹⁵N-T ₁, -T ₂, and {¹H}–¹⁵N NOE data, and reduced spectral density analysis of the relaxation data. The principle features of the notebooks include use of a minimal number of input files, integrated notebook data management, ease of use, cross-platform compatibility, automatic visualization of results and generation of high-quality graphics, and output of analyses in text format.L. Spyracopoulos is an AHFMR Medical Research Senior Scholar     

Residue Specific Ribose and Nucleobase Dynamics of the cUUCGg RNA Tetraloop Motif by MNMR 13C Relaxation

The dynamics of the nucleobase and the ribose moieties in a 14-nt RNA cUUCGg hairpin-loop uniformly labeled with ¹³C and ¹⁵N were studied by ¹³C spin relaxation experiments. R₁, R_1ρ and the ¹³C-{¹H} steady-state NOE of C₆ and C_1′ in pyrimidine and C₈ and C_1′ in purine residues were obtained at 298 K. The relaxation data were analyzed by the model-free formalism to yield dynamic information on timescales of pico-, nano- and milli-seconds. An axially symmetric diffusion tensor with an overall rotational correlation time τ_c of 2.31±0.13 ns and an axial ratio of 1.35±0.02 were determined. Both findings are in agreement with hydrodynamic calculations. For the nucleobase carbons, the validity of different reported ¹³C chemical shift anisotropy values (Stueber, D. and Grant, D. M., 2002 J. Am. Chem. Soc. 124, 10539–10551; Fiala et al., 2000 J. Biomol. NMR 16, 291–302; Sitkoff, D. and Case, D. A., 1998 Prog. NMR Spectroscopy 32, 165–190) is discussed. The resulting dynamics are in agreement with the structural features of the cUUCGg motif in that all residues are mostly rigid (0.82 < S² < 0.96) in both the nucleobase and the ribose moiety except for the nucleobase of U7, which is protruding into solution (S² = 0.76). In general, ribose mobility follows nucleobase dynamics, but is less pronounced. Nucleobase dynamics resulting from the analysis of ¹³C relaxation rates were found to be in agreement with ¹⁵N relaxation data derived dynamic information (Akke et al., 1997 RNA 3, 702–709). Electronic supplementary material Electronic supplementary material is available for this article at and accessible for authorised users.     

Comparison of fast backbone dynamics at amide nitrogen and carbonyl sites in dematin headpiece C-terminal domain and its S74E mutant

We perform a detailed comparison of fast backbone dynamics probed at amide nitrogen versus carbonyl carbon sites for dematin headpiece C-terminal domain (DHP) and its S74E mutant (DHPS74E). Carbonyl dynamics is probed via auto-correlated longitudinal rates and transverse C′/C′-C^α CSA/dipolar and C′/C′–N CSA/dipolar cross-correlated rates, while ¹⁵N data are taken from a previous study. Resulting values of effective order parameters and internal correlation times support the conclusion that C′ relaxation reports on a different subset of fast motions compared to those probed at N–H bond vectors in the same peptide planes. ¹³C′ order parameters are on the average 0.08 lower than ¹⁵N order parameters with the exception of the flexible loop region in DHP. The reduction of mobility in the loop region upon the S74E mutation can be seen from the ¹⁵N order parameters but not from the ¹³C order parameters. Internal correlation times at ¹³C′ sites are on the average an order of magnitude longer than those at ¹⁵N sites for the well-structured C-terminal subdomains, while the more flexible N-terminal subdomains have more comparable average internal correlation times.     

Structure determination of proteins in <Superscript>2</Superscript>H<Subscript>2</Subscript>O solution aided by a deuterium-decoupled 3D HCA(N)CO experiment

We developed an NMR pulse sequence, 3D HCA(N)CO, to correlate the chemical shifts of protein backbone ¹Hα and ¹³Cα to those of ¹³C′ in the preceding residue. By applying ²H decoupling, the experiment was accomplished with high sensitivity comparable to that of HCA(CO)N. When combined with HCACO, HCAN and HCA(CO)N, the HCA(N)CO sequence allows the sequential assignment using backbone ¹³C′ and amide ¹⁵N chemical shifts without resort to backbone amide protons. This assignment strategy was demonstrated for ¹³C/¹⁵N-labeled GB1 dissolved in ²H₂O. The quality of the GB1 structure determined in ²H₂O was similar to that determined in H₂O in spite of significantly smaller number of NOE correlations. Thus this strategy enables the determination of protein structures in ²H₂O or H₂O at high pH values.     

Broadband <Superscript>15</Superscript>N–<Superscript>13</Superscript>C dipolar recoupling via symmetry-based RF pulse schemes at high MAS frequencies

An approach for generating efficient RN_n^{n_S, n_k} {\rm{RN}}_{n}^{\nu_{\rm{S}}, {\nu_{\rm{k}}}} symmetry-based dual channel RF pulse schemes for γ-encoded broadband ¹⁵N–¹³C dipolar recoupling at high magic angle spinning frequencies is presented. The method involves the numerical optimisation of the RF phase-modulation profile of the basic “R” element so as to obtain heteronuclear double quantum dipolar recoupling sequences with satisfactory magnetisation transfer characteristics. The basic “R” element was implemented as a sandwich of a small number of short pulses of equal duration with each pulse characterised by a RF phase and amplitude values. The performance characteristics of the sequences were evaluated via numerical simulations and ¹⁵N–¹³C chemical shift correlation experiments. Employing such ¹³C–¹⁵N double-quantum recoupling sequences and the multiple receiver capabilities available in the current generation of NMR spectrometers, the possibility to simultaneously acquire 3D NCC and CNH chemical shift correlation spectra is also demonstrated.     

Strong coupling effects during <Emphasis Type="Italic">X</Emphasis>-pulse CPMG experiments recorded on heteronuclear <Emphasis Type="Italic">ABX</Emphasis> spin systems: artifacts and a simple solution

Simulation and experiment have been used to establish that significant artifacts can be generated in X-pulse CPMG relaxation dispersion experiments recorded on heteronuclear ABX spin-systems, such as ¹³C _i –¹³C _j –¹H, where ¹³C _i and ¹³C _j are strongly coupled. A qualitative explanation of the origin of these artifacts is presented along with a simple method to significantly reduce them. An application to the measurement of ¹H CPMG relaxation dispersion profiles in an HIV-2 TAR RNA molecule where all ribose sugars are protonated at the 2′ position, deuterated at all other sugar positions and ¹³C labeled at all sugar carbons is presented to illustrate the problems that strong ¹³C–¹³C coupling introduces and a simple solution is proposed.     

HNCO-based measurement of one-bond amide <Superscript>15</Superscript>N-<Superscript>1</Superscript>H couplings with optimized precision

A pair of 3D HNCO-based experiments have been developed with the aim of optimizing the precision of measurement of ¹J_NH couplings. Both pulse sequences record ¹J_NH coupling evolution during the entire constant time interval that ¹⁵N magnetization is dephasing or rephasing with respect to the directly bonded ¹³C′ nucleus, with ¹⁵N¹³C′ multiple quantum coherence maintained during the ¹³C′ evolution period. The first experiment, designed for smaller proteins, produces an apparent doubling of the ¹J_NH coupling without any accompanying increases in line width. The second experiment is a J-scaled TROSY-HNCO experiment in which the ¹J_NH coupling is measured by frequency difference between resonances offset symmetrically about the position of the downfield component of the ¹⁵N doublet (i.e. the TROSY resonance). This experiment delivers significant gains in precision of ¹J_NH coupling measurement compared to existing J-scaled TROSY-HNCO experiments. With the proper choice of acquisition parameters and sufficient sensitivity to acquire a 3D TROSY-HNCO experiment, it is shown that ¹J_NH couplings can be measured with a precision which approaches or exceeds the precision of measurement with which the frequency of the TROSY resonance itself can be determined.     

Addressing the overlap problem in the quantitative analysis of two dimensional NMR spectra: Application to <Superscript>15</Superscript>N relaxation measurements

A quantitative analysis of 2D ¹H-¹⁵N spectra is often complicated by resonance overlap. Here a simple method is presented for resolving overlapped correlations by recording 2D projection planes from HNCO data sets. Applications are presented involving the measurement of ¹⁵N T₁ relaxation rates in a high molecular weight protein, malate synthase G, and in a system that exchanges between folded and unfolded states, the drkN SH3 domain. By supplementing relaxation data recorded in the conventional way as a series of 2D ¹H-¹⁵N data sets with a series of a pair of projection planes the number of dynamics probes is increased significantly for both systems studied.     

Determination of <Superscript>13</Superscript>C CSA Tensors: Extension of the Model-independent Approach to an RNA Kissing Complex Undergoing Anisotropic Rotational Diffusion in Solution

Chemical shift anisotropy (CSA) tensor parameters have been determined for the protonated carbons of the purine bases in an RNA kissing complex in solution by extending the model-independent approach [Fushman, D., Cowburn, D. (1998) J. Am. Chem. Soc. 120, 7109–7110]. A strategy for determining CSA tensor parameters of heteronuclei in isolated X–H two-spin systems (X = ¹³C or ¹⁵N) in molecules undergoing anisotropic rotational diffusion is presented. The original method relies on the fact that the ratio κ₂=R₂^auto/R₂^cross of the transverse auto- and cross-correlated relaxation rates involving the X CSA and the X–H dipolar interaction is independent of parameters related to molecular motion, provided rotational diffusion is isotropic. However, if the overall motion is anisotropic κ₂ depends on the anisotropy D_||/D_⊥ of rotational diffusion. In this paper, the field dependence of both κ₂ and its longitudinal counterpart κ₁=R₁^auto/R₁^cross are determined. For anisotropic rotational diffusion, our calculations show that the average κ_av = 1/2 (κ₁+κ₂), of the ratios is largely independent of the anisotropy parameter D_||/D_⊥. The field dependence of the average ratio κ_av may thus be utilized to determine CSA tensor parameters by a generalized model-independent approach in the case of molecules with an overall motion described by an axially symmetric rotational diffusion tensor.     

Assignment of paramagnetic <Superscript>15</Superscript>N-HSQC spectra by heteronuclear exchange spectroscopy

Paramagnetic metal ions in proteins provide a rich source of structural information, but the resonance assignments required to extract the information can be challenging. Here we demonstrate that paramagnetically shifted ¹⁵N-HSQC cross-peaks can be assigned using N_Z-exchange spectroscopy under conditions in which the paramagnetic form of the protein is in dynamic equilibrium with its diamagnetic form. Even slow exchange of specifically bound metal ions may be detected within the long lifetime of ¹⁵N longitudinal magnetization of large proteins at high magnetic fields. Alternatively, the exchange can be accelerated using an excess of metal ions. In the resulting exchange spectra, paramagnetic ¹⁵N resonances become visible for residues that are not directly observed in a conventional ¹⁵N-HSQC spectrum due to paramagnetic ¹H^N broadening. The experiments are illustrated by the 30 kDa lanthanide-binding ɛ186/θ complex of DNA polymerase III in the presence of sub-stoichiometric amounts of Dy³⁺ or a mixture of Dy³⁺ and La³⁺.     

<Superscript>13</Superscript>C-<Superscript>13</Superscript>C NOESY: A constructive use of <Superscript>13</Superscript>C-<Superscript>13</Superscript>C spin-diffusion

¹³C-¹³C NOESY experiments were performed under long mixing time conditions on reduced human superoxide dismutase (32 kDa, ¹⁵N, ^13C and 70% ²H labeled). ¹³C-¹³C couplings were successfully eliminated through post-processing of in-phase-anti-phase (IPAP) data. It appears that at mixing time _m of 3.0 s the spin diffusion mechanism allows the detection of 96% of the two-bond correlations involving C and C. The interpretation was confirmed by simulations. This approach broadens the range of applicability of ¹³C-¹³C NOESY spectroscopy.     

An investigation of the dynamics of spermine bound to duplex and quadruplex DNA by <Superscript>13</Superscript>C NMR spectroscopy

A detailed analysis of the ¹³C relaxation of ¹³C-labelled spermine bound to duplex and quadruplex DNA is presented. T₁, T₂ and heteronuclear NOE data were collected at four ¹³C frequencies (75.4, 125.7, 150.9 and 201.2 MHz). The data were analyzed in terms of a frequency-dependent order parameter, S ²(ω), to estimate the generalized order parameter and the contributions to the relaxation from different motional frequencies in the picosecond–nanosecond timescale and from any exchange processes that may be occurring on the microsecond–millisecond timescale. The relaxation data was surprisingly similar for spermine bound to two different duplexes and a linear parallel quadruplex. Analysis of the relaxation data from these complexes confirmed the conclusions of previous studies that the dominant motion of spermine is independent of the macroscopic tumbling of the DNA and has an effective correlation time of ∼50 ps. In contrast, spermine bound to a folded antiparallel quadruplex had faster relaxation rates, especially R ₂. As with the other complexes, a fast internal motion of the order of 50 ps makes a substantial contribution to the relaxation. The generalized order parameter for spermine bound to duplex DNA and the linear quadruplex is small but is larger for spermine bound to the folded quadruplex. In the latter case, there is evidence for exchange between at least two populations of spermine occurring on the microsecond–millisecond timescale. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Selective backbone labeling of proteins using {1,2-<Superscript>13</Superscript>C<Subscript>2</Subscript>}-pyruvate as carbon source

A simple isotope labeling approach for selective ¹³C/¹⁵N backbone labeling of proteins is described. Using {1,2-¹³C₂}-pyruvate as the sole carbon source in bacterial growth media, selective incorporation of ¹³C^α-¹³CO spin-pairs into the backbones of protein molecules with medium-to-high levels of ¹³C-enrichment is possible for a subset of 12 amino acids. The isotope labeling scheme has been tested on a pair of proteins—a 7-kDa immunoglobulin binding domain B1 of streptococcal protein G and an 82-kDa enzyme malate synthase G. A number of protein NMR applications are expected to benefit from the {1,2-¹³C₂}-pyruvate based protein production.     

Pressure-dependent <Superscript>13</Superscript>C chemical shifts in proteins: origins and applications

Pressure-dependent ¹³C chemical shifts have been measured for aliphatic carbons in barnase and Protein G. Up to 200 MPa (2 kbar), most shift changes are linear, demonstrating pressure-independent compressibilities. CH₃, CH₂ and CH carbon shifts change on average by +0.23, −0.09 and −0.18 ppm, respectively, due to a combination of bond shortening and changes in bond angles, the latter matching one explanation for the γ-gauche effect. In addition, there is a residue-specific component, arising from both local compression and conformational change. To assess the relative magnitudes of these effects, residue-specific shift changes for protein G were converted into structural restraints and used to calculate the change in structure with pressure, using a genetic algorithm to convert shift changes into dihedral angle restraints. The results demonstrate that residual ¹³Cα shifts are dominated by dihedral angle changes and can be used to calculate structural change, whereas ¹³Cβ shifts retain significant dependence on local compression, making them less useful as structural restraints.     

Quantifying Lipari–Szabo modelfree parameters from 13CO NMR relaxation experiments

It is proposed to obtain effective Lipari–Szabo order parameters and local correlation times for relaxation vectors of protein ¹³CO nuclei by carrying out a ¹³CO-R₁ auto relaxation experiment, a transverse CSA/dipolar cross correlation and a transverse ¹³CO CSA/¹³CO–¹⁵N CSA/dipolar cross correlation experiment. Given the global rotational correlation time from ¹⁵N relaxation experiments, a new program COMFORD (CO-Modelfree Fitting Of Relaxation Data) is presented to fit the ¹³CO data to an effective order parameter , an effective local correlation time and the orientation of the CSA tensor with respect to the molecular frame. It is shown that the effective is least sensitive to rotational fluctuations about an imaginary axis and most sensitive to rotational fluctuations about an imaginary axis parallel to the NH bond direction. As such, the information is fully complementary to the ¹⁵N relaxation order parameter, which is least sensitive to fluctuations about the NH axis and most sensitive to fluctuations about the axis. The new paradigm is applied on data of Ca²⁺ saturated Calmodulin, and on available literature data for Ubiquitin. Our data indicate that the order parameters rapport on slower, and sometimes different, motions than the ¹⁵N relaxation order parameters. The CO local correlation times correlate well with the calmodulin’s secondary structure. Electronic Supplementary Material Supplementary material is available to authorized users in the online version of this article at .     

An intraresidual i(HCA)CO(CA)NH experiment for the assignment of main-chain resonances in <Superscript>15</Superscript>N, <Superscript>13</Superscript>C labeled proteins

An improved pulse sequence, intraresidual i(HCA)CO(CA)NH, is described for establishing solely ¹³C′(i), ¹⁵N(i), ¹H^N(i) connectivities in uniformly ¹⁵N/¹³C-labeled proteins. In comparison to the “out-and-back” style intra-HN(CA)CO experiment, the new pulse sequence offers at least two-fold higher experimental resolution in the ¹³C′ dimension and on average 1.6 times higher sensitivity especially for residues in α-helices. Performance of the new experiment was tested on a small globular protein ubiquitin and an intrinsically unfolded 110-residue cancer/testis antigen CT16/PAGE5. Use of intraresidual i(HCA)CO(CA)NH experiment in combination with the established HNCO experiment was crucial for the assignment of highly disordered CT16.