A simple method to adjust inconsistently referenced <Superscript>13</Superscript>C and <Superscript>15</Superscript>N chemical shift assignments of proteins

Inconsistent ¹³C and ¹⁵N chemical shift referencing is a continuing problem associated with protein chemical shift assignments deposited in BioMagResBank (BMRB). Here we describe a simple and robust approach that can quantitatively determine the ¹³C and ¹⁵N referencing offsets solely from chemical shift assignment data and independently of 3D coordinate data. This novel structure-independent approach permitted the assessment and determination of ¹³C and ¹⁵N reference offsets for all protein entries deposited in the BMRB. Tests on 452 proteins with known 3D structures show that this structure-independent approach yields ¹³C and ¹⁵N referencing offsets that exhibit excellent agreement with those calculated on the basis of 3D structures. Furthermore, this protocol appears to improve the accuracy of chemical shift-derived secondary structural identification, and has been formally incorporated into a computer program called PSSI (http//www.pronmr.com).Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1007/s10858-004-7441-3     

Conformational signatures of <Superscript>13</Superscript>C chemical shifts in RNA ribose

The conformational dependence of ¹³C chemical shift values of RNA riboses determined by liquid-state NMR spectroscopy was evaluated using data deposited for RNA structures in the RCSD and BMRB data bases. Results derived support the applicability of the canonical coordinates approach of Rossi and Harbison (J Magn Reson 151:1–8, 2001) in liquid-state NMR to assess the sugar pucker of ribose units in RNA. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Deuterium isotope effects on <Superscript>15</Superscript>N backbone chemical shifts in proteins

Quantum mechanical calculations are presented that predict that one-bond deuterium isotope effects on the ¹⁵N chemical shift of backbone amides of proteins, ¹Δ¹⁵N(D), are sensitive to backbone conformation and hydrogen bonding. A quantitative empirical model for ¹Δ¹⁵N(D) including the backbone dihedral angles, Φ and Ψ, and the hydrogen bonding geometry is presented for glycine and amino acid residues with aliphatic side chains. The effect of hydrogen bonding is rationalized in part as an electric-field effect on the first derivative of the nuclear shielding with respect to N–H bond length. Another contributing factor is the effect of increased anharmonicity of the N–H stretching vibrational state upon hydrogen bonding, which results in an altered N–H/N–D equilibrium bond length ratio. The N–H stretching anharmonicity contribution falls off with the cosine of the N–H···O bond angle. For residues with uncharged side chains a very good prediction of isotope effects can be made. Thus, for proteins with known secondary structures, ¹Δ¹⁵N(D) can provide insights into hydrogen bonding geometries. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

RefDB: a database of uniformly referenced protein chemical shifts

RefDB is a secondary database of reference-corrected protein chemical shifts derived from the BioMagResBank (BMRB). The database was assembled by using a recently developed program (SHIFTX) to predict protein ¹H, ¹³C and ¹⁵N chemical shifts from X-ray or NMR coordinate data of previously assigned proteins. The predicted shifts were then compared with the corresponding observed shifts and a variety of statistical evaluations performed. In this way, potential mis-assignments, typographical errors and chemical referencing errors could be identified and, in many cases, corrected. This approach allows for an unbiased, instrument-independent solution to the problem of retrospectively re-referencing published protein chemical shifts. Results from this study indicate that nearly 25% of BMRB entries with ¹³C protein assignments and 27% of BMRB entries with ¹⁵N protein assignments required significant chemical shift reference readjustments. Additionally, nearly 40% of protein entries deposited in the BioMagResBank appear to have at least one assignment error. From this study it evident that protein NMR spectroscopists are increasingly adhering to recommended IUPAC ¹³C and ¹⁵N chemical shift referencing conventions, however, approximately 20% of newly deposited protein entries in the BMRB are still being incorrectly referenced. This is cause for some concern. However, the utilization of RefDB and its companion programs may help mitigate this ongoing problem. RefDB is updated weekly and the database, along with its associated software, is freely available at http://redpoll.pharmacy.ualberta.ca and the BMRB website.     

Linear analysis of carbon-13 chemical shift differences and its application to the detection and correction of errors in referencing and spin system identifications

Statistical analysis reveals that the set of differences between the secondary shifts of the α- and β-carbons for residues i of a protein (Δδ¹³C^α_i- Δδ¹³C^β_i) provides the means to detect and correct referencing errors for ¹H and ¹³C nuclei within a given dataset. In a correctly referenced protein dataset, linear regression plots of Δδ¹³C^α_i,Δδ¹³C^β_i, or Δδ¹H^α_i vs. (Δδ¹³C^α_i- Δδ¹³C^β_i) pass through the origin from two directions, the helix-to-coil and strand-to-coil directions. Thus, linear analysis of chemical shifts (LACS) can be used to detect referencing errors and to recalibrate the ¹H and ¹³C chemical shift scales if needed. The analysis requires only that the signals be identified with distinct residue types (intra-residue spin systems). LACS allows errors in calibration to be detected and corrected in advance of sequence-specific assignments and secondary structure determinations. Signals that do not fit the linear model (outliers) deserve scrutiny since they could represent errors in identifying signals with a particular residue, or interesting features such as a cis-peptide bond. LACS provides the basis for the automated detection of such features and for testing reassignment hypotheses. Early detection and correction of errors in referencing and spin system identifications can improve the speed and accuracy of chemical shift assignments and secondary structure determinations. We have used LACS to create a database of offset-corrected chemical shifts corresponding to nearly 1800 BMRB entries: 300 with and 1500 without corresponding three-dimensional (3D) structures. This database can serve as a resource for future analysis of the effects of amino acid sequence and protein secondary and tertiary structure on NMR chemical shifts.Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1007/s10858-005-1717-0     

TALOS+: a hybrid method for predicting protein backbone torsion angles from NMR chemical shifts

NMR chemical shifts in proteins depend strongly on local structure. The program TALOS establishes an empirical relation between ¹³C, ¹⁵N and ¹H chemical shifts and backbone torsion angles ϕ and ψ (Cornilescu et al. J Biomol NMR 13 289–302, 1999). Extension of the original 20-protein database to 200 proteins increased the fraction of residues for which backbone angles could be predicted from 65 to 74%, while reducing the error rate from 3 to 2.5%. Addition of a two-layer neural network filter to the database fragment selection process forms the basis for a new program, TALOS+, which further enhances the prediction rate to 88.5%, without increasing the error rate. Excluding the 2.5% of residues for which TALOS+ makes predictions that strongly differ from those observed in the crystalline state, the accuracy of predicted ϕ and ψ angles, equals ±13°. Large discrepancies between predictions and crystal structures are primarily limited to loop regions, and for the few cases where multiple X-ray structures are available such residues are often found in different states in the different structures. The TALOS+ output includes predictions for individual residues with missing chemical shifts, and the neural network component of the program also predicts secondary structure with good accuracy. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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.     

Automatic <Superscript>13</Superscript>C chemical shift reference correction for unassigned protein NMR spectra

Poor chemical shift referencing, especially for ¹³C in protein Nuclear Magnetic Resonance (NMR) experiments, fundamentally limits and even prevents effective study of biomacromolecules via NMR, including protein structure determination and analysis of protein dynamics. To solve this problem, we constructed a Bayesian probabilistic framework that circumvents the limitations of previous reference correction methods that required protein resonance assignment and/or three-dimensional protein structure. Our algorithm named Bayesian Model Optimized Reference Correction (BaMORC) can detect and correct ¹³C chemical shift referencing errors before the protein resonance assignment step of analysis and without three-dimensional structure. By combining the BaMORC methodology with a new intra-peaklist grouping algorithm, we created a combined method called Unassigned BaMORC that utilizes only unassigned experimental peak lists and the amino acid sequence. Unassigned BaMORC kept all experimental three-dimensional HN(CO)CACB-type peak lists tested within ±?0.4 ppm of the correct ¹³C reference value. On a much larger unassigned chemical shift test set, the base method kept ¹³C chemical shift referencing errors to within ±?0.45 ppm at a 90% confidence interval. With chemical shift assignments, Assigned BaMORC can detect and correct ¹³C chemical shift referencing errors to within ±?0.22 at a 90% confidence interval. Therefore, Unassigned BaMORC can correct ¹³C chemical shift referencing errors when it will have the most impact, right before protein resonance assignment and other downstream analyses are started. After assignment, chemical shift reference correction can be further refined with Assigned BaMORC. These new methods will allow non-NMR experts to detect and correct ¹³C referencing error at critical early data analysis steps, lowering the bar of NMR expertise required for effective protein NMR analysis.     

Automated amino acid side-chain NMR assignment of proteins using <Superscript>13</Superscript>C- and <Superscript>15</Superscript>N-resolved 3D [<Superscript>1</Superscript>H,<Superscript>1</Superscript>H]-NOESY

ASCAN is a new algorithm for automatic sequence-specific NMR assignment of amino acid side-chains in proteins, which uses as input the primary structure of the protein, chemical shift lists of (1)H(N), (15)N, (13)C(alpha), (13)C(beta) and possibly (1)H(alpha) from the previous polypeptide backbone assignment, and one or several 3D (13)C- or (15)N-resolved [(1)H,(1)H]-NOESY spectra. ASCAN has also been laid out for the use of TOCSY-type data sets as supplementary input. The program assigns new resonances based on comparison of the NMR signals expected from the chemical structure with the experimentally observed NOESY peak patterns. The core parts of the algorithm are a procedure for generating expected peak positions, which is based on variable combinations of assigned and unassigned resonances that arise for the different amino acid types during the assignment procedure, and a corresponding set of acceptance criteria for assignments based on the NMR experiments used. Expected patterns of NOESY cross peaks involving unassigned resonances are generated using the list of previously assigned resonances, and tentative chemical shift values for the unassigned signals taken from the BMRB statistics for globular proteins. Use of this approach with the 101-amino acid residue protein FimD(25-125) resulted in 84% of the hydrogen atoms and their covalently bound heavy atoms being assigned with a correctness rate of 90%. Use of these side-chain assignments as input for automated NOE assignment and structure calculation with the ATNOS/CANDID/DYANA program suite yielded structure bundles of comparable quality, in terms of precision and accuracy of the atomic coordinates, as those of a reference structure determined with interactive assignment procedures. A rationale for the high quality of the ASCAN-based structure determination results from an analysis of the distribution of the assigned side chains, which revealed near-complete assignments in the core of the protein, with most of the incompletely assigned residues located at or near the protein surface.     

Applications of Variable-angle Sample Spinning Experiments to the Measurement of Scaled Residual Dipolar Couplings and <Superscript>15</Superscript>N CSA in Soluble Proteins

NMR spectra of ubiquitin in the presence of bicelles at a concentration of 32% w/v have been recorded at 700 MHz under sample spinning conditions at the magic angle (54.7°) and at an angle of 45.5°. At the magic angle, the ¹H–¹⁵N HSQC spectrum of ubiquitin in bicelles is virtually indistinguishable from the one recorded on the protein in solution. Spinning the sample at the magic angle creates an isotropic environment with no preferred bicelle orientations, thus allowing the determination of scalar coupling constants. For an angle of rotation of 45.5°, the bicelles orient with their normal perpendicular to the spinning axis leading to the observation of strong residual dipolar couplings and chemical shift variations of the ¹⁵N resonances. Electronic supplementary material Electronic supplementary material is available for this article at and accessible for authorised users.     

Effects of growth and tissue type on the kinetics of <Superscript>13</Superscript>C and <Superscript>15</Superscript>N incorporation in a rapidly growing ectotherm

The use of stable isotopes to investigate animal diets, habitat use, and trophic level requires understanding the rate at which animals incorporate the ¹³C and ¹⁵N from their diets and the factors that determine the magnitude of the difference in isotopic composition between the animal’s diet and that of its tissues. We determined the contribution of growth and catabolic turnover to the rate of ¹³C and ¹⁵N incorporation into several tissues that can be sampled non-invasively (skin, scute, whole blood, red blood cells, and plasma solutes) in two age classes of a rapidly growing ectotherm (loggerhead turtles, Caretta caretta). We found significant differences in C and N incorporation rates and isotopic discrimination factors (Δ¹³C = δ¹³C_tissues − δ¹³C_diet and Δ¹⁵N = δ¹⁵N_tissues − δ¹⁵N_diet) among tissues and between age classes. Growth explained from 26 to 100% of the total rate of incorporation in hatchling turtles and from 15 to 52% of the total rate of incorporation in juvenile turtles. Because growth contributed significantly to the rate of isotopic incorporation, variation in rates among tissues was lower than reported in previous studies. The contribution of growth can homogenize the rate of isotopic incorporation and limit the application of stable isotopes to identify dietary changes at contrasting time scales and to determine the timing of diet shifts. The isotopic discrimination factor of nitrogen ranged from −0.64 to 1.77‰ in the turtles’ tissues. These values are lower than the commonly assumed average 3.4‰ discrimination factors reported for whole body and muscle isotopic analyses. The increasing reliance on non-invasive and non-destructive sampling in animal isotopic ecology requires that we recognize and understand why different tissues differ in isotopic discrimination factors.     

Backbone <Superscript>1</Superscript>H, <Superscript>13</Superscript>C and <Superscript>15</Superscript>N chemical shift assignment of full-length human uracil DNA glycosylase UNG2

Human uracil N-glycosylase isoform 2—UNG2 consists of an N-terminal intrinsically disordered regulatory domain (UNG2 residues 1–92, 9.3 kDa) and a C-terminal structured catalytic domain (UNG2 residues 93–313, 25.1 kDa). Here, we report the backbone ¹H, ¹³C, and ¹⁵N chemical shift assignment as well as secondary structure analysis of the N-and C-terminal domains of UNG2 representing the full-length UNG2 protein.     

Fish scale δ<Superscript>15</Superscript>N and δ<Superscript>13</Superscript>C values provide biogeochemical tags of fish comparable in performance to element concentrations in scales and otoliths

Natural variability in stable isotope ratios and element concentrations in calcified structures of fish (e.g. scales and otoliths) has provided biogeochemical ‘tags’ for studying origins and movements of marine species, but has been little used in freshwater studies. We examine whether variability in scale δ¹⁵N and δ¹³C values of Salmo trutta L., could provide a tag of fish over small spatial scales in a small river catchment (River Dee, U.K.) and compared their performance as tags with that of scale/otolith element concentrations. Whole scale δ¹⁵N and δ¹³C values differed among six study sites and fish could be classified to their site of origin with a high degree of accuracy. Classifying fish to their site of capture was marginally superior using scale δ¹⁵N and δ¹³C values compared to that achieved using Sr, Mn, Ba and Mg in scale hydroxyapatite or otolith aragonite. Scale δ¹⁵N and δ¹³C values could therefore provide non-lethally collectable biogeochemical tags superior in performance to element concentrations in otoliths and scales. A comprehensive study of δ¹⁵N and δ¹³C values within freshwater systems would develop our understanding of factors influencing geographical variability in baseline δ¹⁵N and δ¹³C signatures.     

<Superscript>1</Superscript>H, <Superscript>13</Superscript>C and <Superscript>15</Superscript>N resonance assignment of human guanylate kinase

Human guanylate kinase (hGMPK) is a critical enzyme that, in addition to phosphorylating its physiological substrate (d)GMP, catalyzes the second phosphorylation step in the conversion of anti-viral and anti-cancer nucleoside analogs to their corresponding active nucleoside analog triphosphates. Until now, a high-resolution structure of hGMPK is unavailable and thus, we studied free hGMPK by NMR and assigned the chemical shift resonances of backbone and side chain ¹H, ¹³C, and ¹⁵N nuclei as a first step towards the enzyme’s structural and mechanistic analysis with atomic resolution.     

<Superscript>1</Superscript>H, <Superscript>13</Superscript>C and <Superscript>15</Superscript>N backbone and side-chain resonance assignments of the ZnF-UBP domain of USP20/VDU2

Deubiquitinase USP20/VDU2 has been identified as a regulator of multiple proteins including hypoxia-inducible factor (HIF)-1α, β2-adrenergic receptor, and tumor necrosis factor receptor associated factor 6 etc. It contains four structural domains, including an N-terminal zinc-finger ubiquitin binding domain (ZnF-UBP) that potentially helps USP20 to recruit its ubiquitin substrates. Here we report the ¹H, ¹³C and ¹⁵N backbone and side-chain resonance assignments of the ZnF-UBP domain of USP20/VDU2. The BMRB accession number is 26901. The secondary structural elements predicted from the NMR data reveal a global fold consisting of three α-helices and four β-strands. The complete assignments can be used to explore the protein dynamics of the USP20 ZnF-UBP and its interactions with monoubiquitin and ubiquitin chains.     

Foliar and fungal <Superscript>15</Superscript> N:<Superscript>14</Superscript> N ratios reflect development of mycorrhizae and nitrogen supply during primary succession: testing analytical models

Nitrogen isotopes (¹⁵N/¹⁴N ratios, expressed as δ¹⁵N values) are useful markers of the mycorrhizal role in plant nitrogen supply because discrimination against ¹⁵N during creation of transfer compounds within mycorrhizal fungi decreases the ¹⁵N/¹⁴N in plants (low δ¹⁵N) and increases the ¹⁵N/¹⁴N of the fungi (high δ¹⁵N). Analytical models of ¹⁵N distribution would be helpful in interpreting δ¹⁵N patterns in fungi and plants. To compare different analytical models, we measured nitrogen isotope patterns in soils, saprotrophic fungi, ectomycorrhizal fungi, and plants with different mycorrhizal habits on a glacier foreland exposed during the last 100 years of glacial retreat and on adjacent non-glaciated terrain. Since plants during early primary succession may have only limited access to propagules of mycorrhizal fungi, we hypothesized that mycorrhizal plants would initially be similar to nonmycorrhizal plants in δ¹⁵N and then decrease, if mycorrhizal colonization were an important factor influencing plant δ¹⁵N. As hypothesized, plants with different mycorrhizal habits initially showed similar δ¹⁵N values (−4 to −6‰ relative to the standard of atmospheric N₂ at 0‰), corresponding to low mycorrhizal colonization in all plant species and an absence of ectomycorrhizal sporocarps. In later successional stages where ectomycorrhizal sporocarps were present, most ectomycorrhizal and ericoid mycorrhizal plants declined by 5–6‰ in δ¹⁵N, suggesting transfer of ¹⁵N-depleted N from fungi to plants. The values recorded (−8 to −11‰) are among the lowest yet observed in vascular plants. In contrast, the δ¹⁵N of nonmycorrhizal plants and arbuscular mycorrhizal plants declined only slightly or not at all. On the forefront, most ectomycorrhizal and saprotrophic fungi were similar in δ¹⁵N (−1 to −3‰), but the host-specific ectomycorrhizal fungus Cortinarius tenebricus had values of up to 7‰. Plants, fungi and soil were at least 4‰ higher in δ¹⁵N from the mature site than in recently exposed sites. On both the forefront and the mature site, host-specific ectomycorrhizal fungi had higher δ¹⁵N values than ectomycorrhizal fungi with a broad host range. From these isotopic patterns, we conclude:(1) large enrichments in ¹⁵N of many ectomycorrhizal fungi relative to co-occurring ectomycorrhizal plants are best explained by treating the plant-fungal-soil system as a closed system with a discrimination against ¹⁵N of 8–10‰ during transfer from fungi to plants, (2) based on models of ¹⁵N mass balance, ericoid and ectomycorrhizal fungi retain up to two-thirds of the N in the plant-mycorrhizal system under the N-limited conditions at forefront sites, (3) sporocarps are probably enriched in ¹⁵N by an additional 3‰ relative to available nitrogen, and (4) host-specific ectomycorrhizal fungi may transfer more N to plant hosts than non-host-specific ectomycorrhizal fungi. Our study confirms that nitrogen isotopes are a powerful tool for probing nitrogen dynamics between mycorrhizal fungi and associated plants.     

<Superscript>1</Superscript>H, <Superscript>15</Superscript>N, <Superscript>13</Superscript>C backbone resonance assignment of the C-terminal domain of enzyme I from <Emphasis Type="Italic">Thermoanaerobacter tengcongensis</Emphasis>

Phosphoenolpyruvate binding to the C-terminal domain (EIC) of enzyme I of the bacterial phosphotransferase system (PTS) initiates a phosphorylation cascade that results in sugar translocation across the cell membrane and controls a large number of essential pathways in bacterial metabolism. EIC undergoes an expanded to compact conformational equilibrium that is regulated by ligand binding and determines the phosphorylation state of the overall PTS. Here, we report the backbone ¹H, ¹⁵N and ¹³C chemical shift assignments of the 70 kDa EIC dimer from the thermophilic bacterium Thermoanaerobacter tengcongensis. Assignments were obtained at 70 °C by heteronuclear multidimensional NMR spectroscopy. In total, 90% of all backbone resonances were assigned, with 264 out of a possible 299 residues assigned in the ¹H–¹⁵N TROSY spectrum. The secondary structure predicted from the assigned backbone resonance using the program TALOS+ is in good agreement with the X-ray crystal structure of T. tengcongensis EIC. The reported assignments will allow detailed structural and thermodynamic investigations on the coupling between ligand binding and conformational dynamics in EIC.     

Foliar <Superscript>15</Superscript>N natural abundance indicates phosphorus limitation of bog species

Foliar δ¹⁵N, %N and %P in the dominant woody and herbaceous species across nutrient gradients in New Zealand restiad (family Restionaceae) raised bogs revealed marked differences in plant δ¹⁵N correlations with P. The two heath shrubs, Leptospermum scoparium (Myrtaceae) and Dracophyllum scoparium (Epacridaceae), showed considerable isotopic variation (−2.03 to −15.55‰, and −0.39 to −12.06‰, respectively) across the bogs, with foliar δ¹⁵N strongly and positively correlated with P concentrations in foliage and peat, and negatively correlated with foliar N:P ratios. For L. scoparium, the isotopic gradient was not linked to ectomycorrhizal (ECM) fractionation as ECMs occurred only on higher nutrient marginal peats where ¹⁵N depletion was least. In strong contrast, restiad species (Empodisma minus Sporadanthus ferrugineus, S. traversii) showed little isotopic variation across the same nutrient gradients. Empodisma minus and S. traversii had δ¹⁵N levels consistently around 0‰ (means of −0.12‰ and +0.15‰ respectively), and S. ferrugineus, which co-habited with E. minus, was more depleted (mean −4.97‰). The isotopic differences between heath shrubs and restiads were similar in floristically dissimilar bogs and may be linked to contrasting nutrient demands, acquisition mechanisms, and root morphology. Leptospermum scoparium shrubs on low nutrient peats were stunted, with low tissue P concentrations, and high N:P ratios, suggesting they were P-limited, which was probably exacerbated by markedly reduced mycorrhizal colonisations. The coupling of δ¹⁵N depletion and %P in heath shrubs suggests that N fractionation is promoted by P limitation. In contrast, the constancy in δ¹⁵N of the restiad species through the N and P gradients suggests that these are not suffering from P limitation.     

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.