CPMG relaxation dispersion NMR experiments measuring glycine 1Hα and 13Cα chemical shifts in the ‘invisible’ excited states of proteins

Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion NMR experiments are extremely powerful for characterizing millisecond time-scale conformational exchange processes in biomolecules. A large number of such CPMG experiments have now emerged for measuring protein backbone chemical shifts of sparsely populated (>0.5%), excited state conformers that cannot be directly detected in NMR spectra and that are invisible to most other biophysical methods as well. A notable deficiency is, however, the absence of CPMG experiments for measurement of ¹H^α and ¹³C^α chemical shifts of glycine residues in the excited state that reflects the fact that in this case the ¹H^α, ¹³C^α spins form a three-spin system that is more complex than the AX ¹H^α–¹³C^α spin systems in the other amino acids. Here pulse sequences for recording ¹H^α and ¹³C^α CPMG relaxation dispersion profiles derived from glycine residues are presented that provide information from which ¹H^α, ¹³C^α chemical shifts can be obtained. The utility of these experiments is demonstrated by an application to a mutant of T4 lysozyme that undergoes a millisecond time-scale exchange process facilitating the binding of hydrophobic ligands to an internal cavity in the protein.     

Predicting <Superscript>13</Superscript>C<Superscript>α</Superscript> chemical shifts for validation of protein structures

The ¹³C^α chemical shifts for 16,299 residues from 213 conformations of four proteins (experimentally determined by X-ray crystallography and Nuclear Magnetic Resonance methods) were computed by using a combination of approaches that includes, but is not limited to, the use of density functional theory. Initially, a validation test of this methodology was carried out by a detailed examination of the correlation between computed and observed ¹³C^α chemical shifts of 10,564 (of the 16,299) residues from 139 conformations of the human protein ubiquitin. The results of this validation test on ubiquitin show agreement with conclusions derived from computation of the chemical shifts at the ab initio Hartree–Fock level. Further, application of this methodology to 5,735 residues from 74 conformations of the three remaining proteins that differ in their number of amino acid residues, sequence and three-dimensional structure, together with a new scoring function, namely the conformationally averaged root-mean-square-deviation, enables us to: (a) offer a criterion for an accurate assessment of the quality of NMR-derived protein conformations; (b) examine whether X-ray or NMR-solved structures are better representations of the observed ¹³C^α chemical shifts in solution; (c) provide evidence indicating that the proposed methodology is more accurate than automated predictors for validation of protein structures; (d) shed light as to whether the agreement between computed and observed ¹³C^α chemical shifts is influenced by the identity of an amino acid residue or its location in the sequence; and (e) provide evidence confirming the presence of dynamics for proteins in solution, and hence showing that an ensemble of conformations is a better representation of the structure in solution than any single conformation. Electronic Supplementary Material The online version of this article (doi: ) contains supplementary material, which is available to authorized users.     

Sequential nearest-neighbor effects on computed <Superscript>13</Superscript>C<Superscript>α</Superscript> chemical shifts

To evaluate sequential nearest-neighbor effects on quantum-chemical calculations of ¹³C^α chemical shifts, we selected the structure of the nucleic acid binding (NAB) protein from the SARS coronavirus determined by NMR in solution (PDB id 2K87). NAB is a 116-residue α/β protein, which contains 9 prolines and has 50% of its residues located in loops and turns. Overall, the results presented here show that sizeable nearest-neighbor effects are seen only for residues preceding proline, where Pro introduces an overestimation, on average, of 1.73 ppm in the computed ¹³C^α chemical shifts. A new ensemble of 20 conformers representing the NMR structure of the NAB, which was calculated with an input containing backbone torsion angle constraints derived from the theoretical ¹³C^α chemical shifts as supplementary data to the NOE distance constraints, exhibits very similar topology and comparable agreement with the NOE constraints as the published NMR structure. However, the two structures differ in the patterns of differences between observed and computed ¹³C^α chemical shifts, Δ _ca,i , for the individual residues along the sequence. This indicates that the Δ _ca,i -values for the NAB protein are primarily a consequence of the limited sampling by the bundles of 20 conformers used, as in common practice, to represent the two NMR structures, rather than of local flaws in the structures.     

Improved accuracy in measuring one-bond and two-bond <Superscript>15</Superscript>N,<Superscript>13</Superscript>C<Superscript>α</Superscript> coupling constants in proteins by double-inphase/antiphase (DIPAP) spectroscopy

An extension to HN(CO-α/β-N,C^α-J)-TROSY (Permi and Annila in J Biomol NMR 16:221–227, 2000) is proposed that permits the simultaneous determination of the four coupling constants ¹ J _{N′(i)Cα(i)}, ² J _HN(i)Cα(i), ² J _{Cα(i−1)N′(i)}, and ³ J _{Cα(i−1)HN(i)} in ¹⁵N,¹³C-labeled proteins. Contrasting the original scheme, in which two separate subspectra exhibit the ² J _CαN′ coupling as inphase and antiphase splitting (IPAP), we here record four subspectra that exhibit all combinations of inphase and antiphase splittings possible with respect to both ² J _CαN′ and ¹ J _N′Cα (DIPAP). Complementary sign patterns in the different spectrum constituents overdetermine the coupling constants which can thus be extracted at higher accuracy than is possible with the original experiment. Fully exploiting data redundance, simultaneous 2D lineshape fitting of the E.COSY multiplet tilts in all four subspectra provides all coupling constants at ultimate precision. Cross-correlation and differential-relaxation effects were taken into account in the evaluation procedure. By applying a four-point Fourier transform, the set of spectra is reversibly interconverted between DIPAP and spin-state representations. Methods are exemplified using proteins of various size.     

Measuring the signs of the methyl <Superscript>1</Superscript>H chemical shift differences between major and ‘invisible’ minor protein conformational states using methyl <Superscript>1</Superscript>H multi-quantum spectroscopy

Carr–Purcell–Meiboom–Gill (CPMG) type relaxation dispersion experiments are now routinely used to characterise protein conformational dynamics that occurs on the μs to millisecond (ms) timescale between a visible major state and ‘invisible’ minor states. The exchange rate(s) (\( k_{{{\text{ex}}}} \)), population(s) of the minor state(s) and the absolute value of the chemical shift difference \(|{\Delta \varpi }|\) (ppm) between different exchanging states can be extracted from the CPMG data. However the sign of \({\Delta \varpi }\) that is required to reconstruct the spectrum of the ‘invisible’ minor state(s) cannot be obtained from CPMG data alone. Building upon the recently developed triple quantum (TQ) methyl \( ^{1} {\text{H}} \) CPMG experiment (Yuwen in Angew Chem 55:11490–11494, 2016) we have developed pulse sequences that use carbon detection to generate and evolve single quantum (SQ), double quantum (DQ) and TQ coherences from methyl protons in the indirect dimension to measure the chemical exchange-induced shifts of the SQ, DQ and TQ coherences from which the sign of \({\Delta \varpi }\) is readily obtained for two state exchange. Further a combined analysis of the CPMG data and the difference in exchange induced shifts between the SQ and DQ resonances and between the SQ and TQ resonances improves the estimates of exchange parameters like the population of the minor state. We demonstrate the use of these experiments on two proteins undergoing exchange: (1) the ~ 18 kDa cavity mutant of T4 Lysozyme (\( k_{{{\text{ex}}}} \sim\,3500{\text{ s}}^{{ - 1}} \)) and (2) the \(\sim\,4.7\) kDa Peripheral Sub-unit Binding Domain (PSBD) from the acetyl transferase of Bacillus stearothermophilus (\(k_{ex} \sim\,13,000\hbox { s}^{-1}\)).     

Identification of C-terminal neighbours of amino acid residues without an aliphatic <Superscript>13</Superscript>C<Superscript>γ</Superscript> as an aid to NMR assignments in proteins

We propose a methodology that uses GFT (3,2)D CB(CACO)NNH experiment to rapidly collect the data and readily identify six amino acid residue types (Ala, Asn, Asp, Cys, Gly and Ser) in any given protein. Further, the experiment can distinguish the redox state of Cys residues. The proposed experiment in its two forms will have wide range of applications in resonance assignment strategies and structure determination of proteins.     

A <Superscript>13</Superscript>C-detected <Superscript>15</Superscript>N double-quantum NMR experiment to probe arginine side-chain guanidinium <Superscript>15</Superscript>N<Superscript>η</Superscript> chemical shifts

Arginine side-chains are often key for enzyme catalysis, protein–ligand and protein–protein interactions. The importance of arginine stems from the ability of the terminal guanidinium group to form many key interactions, such as hydrogen bonds and salt bridges, as well as its perpetual positive charge. We present here an arginine ¹³C^ζ-detected NMR experiment in which a double-quantum coherence involving the two ¹⁵N^η nuclei is evolved during the indirect chemical shift evolution period. As the precession frequency of the double-quantum coherence is insensitive to exchange of the two ¹⁵N^η; this new approach is shown to eliminate the previously deleterious line broadenings of ¹⁵N^η resonances caused by the partially restricted rotation about the C^ζ–N^ε bond. Consequently, sharp and well-resolved ¹⁵N^η resonances can be observed. The utility of the presented method is demonstrated on the L99A mutant of the 19 kDa protein T4 lysozyme, where the measurement of small chemical shift perturbations, such as one-bond deuterium isotope shifts, of the arginine amine ¹⁵N^η nuclei becomes possible using the double-quantum experiment.     

Rapid measurement of long-range distances in proteins by multidimensional <Superscript>13</Superscript>C–<Superscript>19</Superscript>F REDOR NMR under fast magic-angle spinning

The ability to simultaneously measure many long-range distances is critical to efficient and accurate determination of protein structures by solid-state NMR (SSNMR). So far, the most common distance constraints for proteins are ¹³C–¹⁵N distances, which are usually measured using the rotational-echo double-resonance (REDOR) technique. However, these measurements are restricted to distances of up to ~?5 Å due to the low gyromagnetic ratios of ¹⁵N and ¹³C. Here we present a robust 2D ¹³C–¹⁹F REDOR experiment to measure multiple distances to ~?10 Å. The technique targets proteins that contain a small number of recombinantly or synthetically incorporated fluorines. The ¹³C–¹⁹F REDOR sequence is combined with 2D ¹³C–¹³C correlation to resolve multiple distances in highly ¹³C-labeled proteins. We show that, at the high magnetic fields which are important for obtaining well resolved ¹³C spectra, the deleterious effect of the large ¹⁹F chemical shift anisotropy for REDOR is ameliorated by fast magic-angle spinning and is further taken into account in numerical simulations. We demonstrate this 2D ¹³C–¹³C resolved ¹³C–¹⁹F REDOR technique on ¹³C, ¹⁵N-labeled GB1. A 5-¹⁹F-Trp tagged GB1 sample shows the extraction of distances to a single fluorine atom, while a 3-¹⁹F-Tyr labeled GB1 sample allows us to evaluate the effects of multi-spin coupling and statistical ¹⁹F labeling on distance measurement. Finally, we apply this 2D REDOR experiment to membrane-bound influenza B M2 transmembrane peptide, and show that the distance between the proton-selective histidine residue and the gating tryptophan residue differs from the distances in the solution NMR structure of detergent-bound BM2. This 2D ¹³C–¹⁹F REDOR technique should facilitate SSNMR-based protein structure determination by increasing the measurable distances to the ~?10 Å range.     

<Superscript>1</Superscript>H–<Superscript>15</Superscript>N correlation spectroscopy of nanocrystalline proteins

The limits of resolution that can be obtained in ¹H–¹⁵N 2D NMR spectroscopy of isotopically enriched nanocrystalline proteins are explored. Combinations of frequency switched Lee–Goldburg (FSLG) decoupling, fast magic angle sample spinning (MAS), and isotopic dilution via deuteration are investigated as methods for narrowing the amide ¹H resonances. Heteronuclear decoupling of ¹⁵N from the ¹H resonances is also studied. Using human ubiquitin as a model system, the best resolution is most easily obtained with uniformly ²H and ¹⁵N enriched protein where the amides have been exchanged in normal water, MAS at 20 kHz, and WALTZ-16 decoupling of the ¹⁵N nuclei. The combination of these techniques results in average ¹H lines of only 0.26 ppm full width at half maximum. Techniques for optimizing instrument stability and ¹⁵N decoupling are described for achieving the best possible performance in these experiments.     

Assessment of chemical exchange in tryptophan–albumin solution through <Superscript>19</Superscript>F multicomponent transverse relaxation dispersion analysis

A number of NMR methods possess the capability of probing chemical exchange dynamics in solution. However, certain drawbacks limit the applications of these NMR approaches, particularly, to a complex system. Here, we propose a procedure that integrates the regularized nonnegative least squares (NNLS) analysis of multiexponential T₂ relaxation into Carr–Purcell–Meiboom–Gill (CPMG) relaxation dispersion experiments to probe chemical exchange in a multicompartmental system. The proposed procedure was validated through analysis of ¹⁹F T₂ relaxation data of 6-fluoro-DL-tryptophan in a two-compartment solution with and without bovine serum albumin. Given the regularized NNLS analysis of a T₂ relaxation curve acquired, for example, at the CPMG frequency υ _CPMG  = 125, the nature of two distinct peaks in the associated T₂ distribution spectrum indicated 6-fluoro-DL-tryptophan either retaining the free state, with geometric mean */multiplicative standard deviation (MSD) = 1851.2 ms */1.51, or undergoing free/albumin-bound interconversion, with geometric mean */MSD = 236.8 ms */1.54, in the two-compartment system. Quantities of the individual tryptophan species were accurately reflected by the associated T₂ peak areas, with an interconversion state-to-free state ratio of 0.45 ± 0.11. Furthermore, the CPMG relaxation dispersion analysis estimated the exchange rate between the free and albumin-bound states in this fluorinated tryptophan analog and the corresponding dissociation constant of the fluorinated tryptophan–albumin complex in the chemical-exchanging, two-compartment system.     

Effects of side-chain orientation on the <Superscript>13</Superscript>C chemical shifts of antiparallel β-sheet model peptides

The dependence of the ¹³C chemical shift on side-chain orientation was investigated at the density functional level for a two-strand antiparallel β-sheet model peptide represented by the amino acid sequence Ac-(Ala)₃-X-(Ala)₁₂-NH₂ where X represents any of the 17 naturally occurring amino acids, i.e., not including alanine, glycine and proline. The dihedral angles adopted for the backbone were taken from, and fixed at, observed experimental values of an antiparallel β-sheet. We carried out a cluster analysis of the ensembles of conformations generated by considering the side-chain dihedral angles for each residue X as variables, and use them to compute the ¹³C chemical shifts at the density functional theory level. It is shown that the adoption of the locally-dense basis set approach for the quantum chemical calculations enabled us to reduce the length of the chemical-shift calculations while maintaining good accuracy of the results. For the 17 naturally occurring amino acids in an antiparallel β-sheet, there is (i) good agreement between computed and observed ¹³C^α and ¹³C^β chemical shifts, with correlation coefficients of 0.95 and 0.99, respectively; (ii) significant variability of the computed ¹³C^α and ¹³C^β chemical shifts as a function of χ¹ for all amino acid residues except Ser; and (iii) a smaller, although significant, dependence of the computed ¹³C^α chemical shifts on χ^ξ (with ξ ≥ 2) compared to χ¹ for eleven out of seventeen residues. Our results suggest that predicted ¹³C^α and ¹³C^β chemical shifts, based only on backbone (φ,ψ) dihedral angles from high-resolution X-ray structure data or from NMR-derived models, may differ significantly from those observed in solution if the dihedral-angle preferences for the side chains are not taken into account. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Major groove width variations in RNA structures determined by NMR and impact of <Superscript>13</Superscript>C residual chemical shift anisotropy and <Superscript>1</Superscript>H–<Superscript>13</Superscript>C residual dipolar coupling on refinement

Ribonucleic acid structure determination by NMR spectroscopy relies primarily on local structural restraints provided by ¹H– ¹H NOEs and J-couplings. When employed loosely, these restraints are broadly compatible with A- and B-like helical geometries and give rise to calculated structures that are highly sensitive to the force fields employed during refinement. A survey of recently reported NMR structures reveals significant variations in helical parameters, particularly the major groove width. Although helical parameters observed in high-resolution X-ray crystal structures of isolated A-form RNA helices are sensitive to crystal packing effects, variations among the published X-ray structures are significantly smaller than those observed in NMR structures. Here we show that restraints derived from aromatic ¹H– ¹³C residual dipolar couplings (RDCs) and residual chemical shift anisotropies (RCSAs) can overcome NMR restraint and force field deficiencies and afford structures with helical properties similar to those observed in high-resolution X-ray structures.     

<Superscript>2</Superscript>H–<Superscript>13</Superscript>C correlation solid-state NMR for investigating dynamics and water accessibilities of proteins and carbohydrates

Site-specific determination of molecular motion and water accessibility by indirect detection of ²H NMR spectra has advantages over dipolar-coupling based techniques due to the large quadrupolar couplings and the ensuing high angular resolution. Recently, a Rotor Echo Short Pulse IRrAdiaTION mediated cross polarization (^RESPIRATIONCP) technique was developed, which allowed efficient transfer of ²H magnetization to ¹³C at moderate ²H radiofrequency field strengths available on most commercial MAS probes. In this work, we investigate the ²H–¹³C magnetization transfer characteristics of one-bond perdeuterated CD _n spin systems and two-bond H/D exchanged C–(O)–D and C–(N)–D spin systems in carbohydrates and proteins. Our results show that multi-bond, broadband ²H–¹³C polarization transfer can be achieved using ²H radiofrequency fields of ~50 kHz, relatively short contact times of 1.3–1.7 ms, and with sufficiently high sensitivity to enable 2D ²H–¹³C correlation experiments with undistorted ²H spectra in the indirect dimension. To demonstrate the utility of this ²H–¹³C technique for studying molecular motion, we show ²H–¹³C correlation spectra of perdeuterated bacterial cellulose, whose surface glucan chains exhibit a motionally averaged C6 ²H quadrupolar coupling that indicates fast trans-gauche isomerization about the C5–C6 bond. In comparison, the interior chains in the microfibril core are fully immobilized. Application of the ²H–¹³C correlation experiment to H/D exchanged Arabidopsis primary cell walls show that the O–D quadrupolar spectra of the highest polysaccharide peaks can be fit to a two-component model, in which 74% of the spectral intensity, assigned to cellulose, has a near-rigid-limit coupling, while 26% of the intensity, assigned to matrix polysaccharides, has a weakened coupling of 50 kHz. The latter O–D quadrupolar order parameter of 0.22 is significantly smaller than previously reported C–D dipolar order parameters of 0.46–0.55 for pectins, suggesting that additional motions exist at the C–O bonds in the wall polysaccharides. ²H–¹³C polarization transfer profiles are also compared between statistically deuterated and H/D exchanged GB1.     

Comprehensive determination of <Superscript>3</Superscript>J<Subscript>HNHα</Subscript> for unfolded proteins using <Superscript>13</Superscript>C′-resolved spin-echo difference spectroscopy

An experiment is presented to determine ³J_HNHα coupling constants, with significant advantages for applications to unfolded proteins. The determination of coupling constants for the peptide chain using 1D ¹H, or 2D and 3D ¹H-¹⁵N correlation spectroscopy is often hampered by extensive resonance overlap when dealing with flexible, disordered proteins. In the experiment detailed here, the overlap problem is largely circumvented by recording ¹H-¹³C′ correlation spectra, which demonstrate superior resolution for unfolded proteins. J-coupling constants are extracted from the peak intensities in a pair of 2D spin-echo difference experiments, affording rapid acquisition of the coupling data. In an application to the cytoplasmic domain of human neuroligin-3 (hNlg3cyt) data were obtained for 78 residues, compared to 54 coupling constants obtained from a 3D HNHA experiment. The coupling constants suggest that hNlg3cyt is intrinsically disordered, with little propensity for structure.     

Fractional <Superscript>13</Superscript>C enrichment of isolated carbons using [1-<Superscript>13</Superscript>C]- or [2-<Superscript>13</Superscript>C]-glucose facilitates the accurate measurement of dynamics at backbone C<Superscript>α</Superscript> and side-chain methyl positions in proteins

A simple labeling approach is presented based on protein expression in [1-¹³C]- or [2-¹³C]-glucose containing media that produces molecules enriched at methyl carbon positions or backbone C^α sites, respectively. All of the methyl groups, with the exception of Thr and Ile(δ1) are produced with isolated ¹³C spins (i.e., no ¹³C–¹³C one bond couplings), facilitating studies of dynamics through the use of spin-spin relaxation experiments without artifacts introduced by evolution due to large homonuclear scalar couplings. Carbon-α sites are labeled without concomitant labeling at C^β positions for 17 of the common 20 amino acids and there are no cases for which ¹³C^α−¹³CO spin pairs are observed. A large number of probes are thus available for the study of protein dynamics with the results obtained complimenting those from more traditional backbone ¹⁵N studies. The utility of the labeling is established by recording ¹³C R _1ρ and CPMG-based experiments on a number of different protein systems.     

Impact of two-bond <Superscript>15</Superscript>N–<Superscript>15</Superscript>N scalar couplings on <Superscript>15</Superscript>N transverse relaxation measurements for arginine side chains of proteins

NMR relaxation of arginine (Arg) ¹⁵N_ε nuclei is useful for studying side-chain dynamics of proteins. In this work, we studied the impact of two geminal ¹⁵N–¹⁵N scalar couplings on measurements of transverse relaxation rates (R ₂ ) for Arg side-chain ¹⁵N_ε nuclei. For 12 Arg side chains of the DNA-binding domain of the Antp protein, we measured the geminal ¹⁵N–¹⁵N couplings ( ² J _NN ) of the ¹⁵N_ε nuclei and found that the magnitudes of the ² J _NN coupling constants were virtually uniform with an average of 1.2 Hz. Our simulations, assuming ideal 180° rotations for all ¹⁵N nuclei, suggested that the two ² J _NN couplings of this magnitude could in principle cause significant modulation in signal intensities during the Carr–Purcell-Meiboom–Gill (CPMG) scheme for Arg ¹⁵N_ε R ₂ measurements. However, our experimental data show that the expected modulation via two ² J _NN couplings vanishes during the ¹⁵N CPMG scheme. This quenching of J modulation can be explained by the mechanism described in Dittmer and Bodenhausen (Chemphyschem 7:831–836, 2006). This effect allows for accurate measurements of R ₂ relaxation rates for Arg side-chain ¹⁵N_ε nuclei despite the presence of two ² J _NN couplings. Although the so-called recoupling conditions may cause overestimate of R ₂ rates for very mobile Arg side chains, such conditions can readily be avoided through appropriate experimental settings.     

H<Superscript>N</Superscript>, N,C<Superscript>α</Superscript> and C<Superscript>β</Superscript> assignments of the two periplasmic domains of <Emphasis Type="Italic">Neisseria meningitidis</Emphasis> DsbD

DsbD is a disulfide bond reductase present in the inner membrane of many Gamma-Proteobacteria. In the human pathogen Neisseria meningitidis, DsbD is required for viability and represents a potential target for the development of antibiotics. Here we report the chemical shift assignments (H^N, N, C^α and C^β) for the reduced and oxidized forms of the two periplasmic domains of N. meningitidis DsbD, n-NmDsbD and c-NmDsbD. The backbone amide resonances in all four forms were completely assigned, and the secondary structures for the core regions of the proteins were calculated using ¹³C^αβ shifts. The reduced and oxidized forms of each domain have similar secondary shifts suggesting they retain the same fold. We anticipate that these data will provide an important basis for studying the interaction between n-NmDsbD and c-NmDsbD, which is required for electron transfer across the bacterial cytoplasmic membrane.     

Facile measurement of <Superscript>1</Superscript>H–<Superscript>15</Superscript>N residual dipolar couplings in larger perdeuterated proteins

We present a simple method, ARTSY, for extracting ¹J_NH couplings and ¹H–¹⁵N RDCs from an interleaved set of two-dimensional ¹H–¹⁵N TROSY-HSQC spectra, based on the principle of quantitative J correlation. The primary advantage of the ARTSY method over other methods is the ability to measure couplings without scaling peak positions or altering the narrow line widths characteristic of TROSY spectra. Accuracy of the method is demonstrated for the model system GB3. Application to the catalytic core domain of HIV integrase, a 36 kDa homodimer with unfavorable spectral characteristics, demonstrates its practical utility. Precision of the RDC measurement is limited by the signal-to-noise ratio, S/N, achievable in the 2D TROSY-HSQC spectrum, and is approximately given by 30/(S/N) Hz.