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.     

CACA-TOCSY with alternate <Superscript>13</Superscript>C–<Superscript>12</Superscript>C labeling: a <Superscript>13</Superscript>C<Superscript>α</Superscript> direct detection experiment for mainchain resonance assignment,dihedral angle information,and amino acid type identification

We present a ¹³C direct detection CACA-TOCSY experiment for samples with alternate ¹³C–¹²C labeling. It provides inter-residue correlations between ¹³C^α resonances of residue i and adjacent C^αs at positions i − 1 and i + 1. Furthermore, longer mixing times yield correlations to C^α nuclei separated by more than one residue. The experiment also provides C^α-to-sidechain correlations, some amino acid type identifications and estimates for ψ dihedral angles. The power of the experiment derives from the alternate ¹³C–¹²C labeling with [1,3-¹³C] glycerol or [2-¹³C] glycerol, which allows utilizing the small scalar ³J_CC couplings that are masked by strong ¹J_CC couplings in uniformly ¹³C labeled samples.     

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.     

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.     

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 <Superscript>13</Superscript>C<Superscript>β</Superscript> chemical shifts of invisible excited states in proteins by relaxation dispersion NMR spectroscopy

A labeling scheme is introduced that facilitates the measurement of accurate ¹³C^β chemical shifts of invisible, excited states of proteins by relaxation dispersion NMR spectroscopy. The approach makes use of protein over-expression in a strain of E. coli in which the TCA cycle enzyme succinate dehydrogenase is knocked out, leading to the production of samples with high levels of ¹³C enrichment (30–40%) at C^β side-chain carbon positions for 15 of the amino acids with little ¹³C label at positions one bond removed (≈5%). A pair of samples are produced using [1-¹³C]-glucose/NaH¹²CO₃ or [2-¹³C]-glucose as carbon sources with isolated and enriched (>30%) ¹³C^β positions for 11 and 4 residues, respectively. The efficacy of the labeling procedure is established by NMR spectroscopy. The utility of such samples for measurement of ¹³C^β chemical shifts of invisible, excited states in exchange with visible, ground conformations is confirmed by relaxation dispersion studies of a protein–ligand binding exchange reaction in which the extracted chemical shift differences from dispersion profiles compare favorably with those obtained directly from measurements on ligand free and fully bound protein samples.     

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.     

Estimation of denitrification potential in a karst aquifer using the <Superscript>15</Superscript>N and <Superscript>18</Superscript>O isotopes of NO<Stack><Subscript>3</Subscript><Superscript>−</Superscript></Stack>

A confined aquifer in the Malm Karst of the Franconian Alb, South Germany was investigated in order to understand the role of the vadose zone in denitrifiaction processes. The concentrations of chemical tracers Sr²⁺ and Cl^– and concentrations of stable isotope ¹⁸O were measured in spring water and precipitation during storm events. Based on these measurements a conceptual model for runoff was constructed. The results indicate that pre-event water, already stored in the system at the beginning of the event, flows downslope on vertical and lateral preferential flow paths. Chemical tracers used in a mixing model for hydrograph separation have shown that the pre-event water contribution is up to 30%. Applying this information to a conceptual runoff generation model, the values of ¹⁵N and ¹⁸O in nitrate could be calculated. Field observations showed the occurence of significant microbial denitrification processes above the soil/bedrock interface before nitrate percolates through to the deeper horizon of the vadose zone. The source of nitrate could be determined and denitrification processes were calculated. Assuming that the nitrate reduction follows a Rayleigh process one could approximate a nitrate input concentration of about 170 mg/l and a residual nitrate concentration of only about 15%. The results of the chemical and isotopic tracers postulate fertilizers as nitrate source with some influence of atmospheric nitrate. The combined application of hydrograph separation and determination of isotope values in ¹⁵N and ¹⁸O of nitrate lead to an improved understanding of microbial processes (nitrification, denitrification) in dynamic systems.     

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.     

Carbonyl carbon label selective (CCLS) <Superscript>1</Superscript>H–<Superscript>15</Superscript>N HSQC experiment for improved detection of backbone <Superscript>13</Superscript>C–<Superscript>15</Superscript>N cross peaks in larger proteins

We present a highly sensitive pulse sequence, carbonyl carbon label selective ¹H–¹⁵N HSQC (CCLS-HSQC) for the detection of signals from ¹H–¹⁵N units involved in ¹³C′–¹⁵N linkages. The CCLS-HSQC pulse sequence utilizes a modified ¹⁵N CT evolution period equal to 1/( ) (∼33 ms) to select for ¹³C′–¹⁵N pairs. By collecting CCLS-HSQC and HNCO data for two proteins (8 kDa ubiquitin and 20 kDa HscB) at various temperatures (5–40°C) in order to vary correlation times, we demonstrate the superiority of the CCLS-HSQC pulse sequence for proteins with long correlation times (i.e. higher molecular weight). We then show that the CCLS-HSQC experiment yields assignments in the case of a 41 kDa protein incorporating pairs of ¹⁵N- and ¹³C′-labeled amino acids, where a TROSY 2D-HN(CO) had failed. Although the approach requires that the ¹H–¹⁵N HSQC cross peaks be observable, it does not require deuteration of the protein. The method is suitable for larger proteins and is less affected by conformational exchange than HNCO experiments, which require a longer period of transverse ¹⁵N magnetization. The method also is tolerant to the partial loss of signal from isotopic dilution (scrambling). This approach will be applicable to families of proteins that have been resistant to NMR structural and dynamic analysis, such as large enzymes, and partially folded or unfolded proteins.     

Simultaneous α/β spin-state selection for <Superscript>13</Superscript>C and <Superscript>15</Superscript>N from a time-shared HSQC-IPAP experiment

Two novel HSQC-IPAP approaches are proposed to achieve α/β spin-state editing simultaneously for ¹³C and ¹⁵N in a single NMR experiment. The pulse schemes are based on a time-shared (TS) 2D ¹H,¹³C/¹H,¹⁵N-HSQC correlation experiment that combines concatenated echo elements for simultaneous J(CH) and J(NH) coupling constants evolution, TS evolution of ¹³C and ¹⁵N chemical shifts in the indirect dimension and heteronuclear α/β-spin-state selection by means of the IPAP principle. Heteronuclear α/β-editing for all CH _n (n = 1–3) and NH _n (1–2) multiplicities can be achieved in the detected F2 dimension of a single TS-HSQC-F2-IPAP experiment. On the other hand, an alternative TS-HSQC-F1-IPAP experiment is also proposed to achieve α/β-editing in the indirect F1 dimension. Experimental and simulated data is provided to evaluate these principles in terms of sensitivity and performance simultaneously on backbone and side-chain CH, CH₂, CH₃, NH, and NH₂ spin systems in uniformly ¹³C/¹⁵N-labeled proteins and in small natural-abundance peptides.     

<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.     

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.     

<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.     

Variability and directionality of temporal changes in δ<Superscript>13</Superscript>C and δ<Superscript>15</Superscript>N of aquatic invertebrate primary consumers

Seasonal oscillations in the carbon (δ¹³C) and nitrogen (δ¹⁵N) isotope signatures of aquatic algae can cause seasonal enrichment–depletion cycles in the isotopic composition of planktonic invertebrates (e.g., copepods). Yet, there is growing evidence that seasonal enrichment–depletion cycles also occur in the isotope signatures of larger invertebrate consumers, taxa used to define reference points in isotope-based trophic models (e.g., trophic baselines). To evaluate the general assumption of temporal stability in non-zooplankton aquatic invertebrates, δ¹³C and δ¹⁵N time series data from the literature were analyzed for seasonality and the influence of biotic (feeding group) and abiotic (trophic state, climate regime) factors on isotope temporal patterns. The amplitude of δ¹³C and δ¹⁵N enrichment–depletion cycles was negatively related to body size, although all size-classes of invertebrates displayed a winter-to-summer enrichment in δ¹³C and depletion in δ¹⁵N. Among feeding groups, periphytic grazers were more variable and displayed larger temporal changes in δ¹³C than detritivores. For nitrogen, temporal variability and magnitude of directional change of δ¹⁵N was most strongly related to ecosystem trophic state (eutrophic > mesotrophic, oligotrophic). This study provides evidence of seasonality in the isotopic composition of aquatic invertebrates across very broad geographical and ecological gradients as well as identifying factors that are likely to modulate the strength and variability of seasonality. These results emphasize the need for researchers to recognize the likelihood of temporal changes in non-zooplankton aquatic invertebrate consumers at time scales relevant to seasonal studies and, if present, to account for temporal dynamics in isotope trophic models.     

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.     

Sequence specific <Superscript>1</Superscript>H, <Superscript>13</Superscript>C and <Superscript>15</Superscript>N resonance assignments of a cataract-related variant G57W of human γS-crystallin

γS-crystallin is a major structural component of the human eye lens, which maintains its stability over the lifetime of an organism with negligible turnover. The G57W mutant of human γS-crystallin (abbreviated hereafter as γS-G57W) is associated with dominant congenital cataracts. In order to provide a structural basis for the ability of γS-G57W causing cataract, we have cloned, overexpressed, isolated and purified the protein. The 2D [¹⁵N–¹H]-HSQC spectrum recorded with uniformly ¹³C/¹⁵N-labelled γS-G57W was highly dispersed indicating the protein to adopt an ordered conformation. In this paper, we report almost complete sequence-specific ¹H, ¹³C and ¹⁵N resonance assignments of γS-G57W using a suite of heteronuclear 3D NMR experiments.     

Spatial variation in flux, δ<Superscript>13</Superscript>C and δ<Superscript>2</Superscript>H of methane in a small Arctic lake with fringing wetland in western Greenland

Small lakes in northern latitudes represent a significant source of CH₄ to the atmosphere that is predicted to increase with warming in the Arctic. Yet, whole-lake CH₄ budgets are lacking as are measurements of δ¹³C-CH₄ and δ²H-CH₄. In this study, we quantify spatial variability of diffusive and ebullitive fluxes of CH₄ and corresponding δ¹³C-CH₄ and δ²H-CH₄ in a small, Arctic lake system with fringing wetland in southwestern Greenland during summer. Net CH₄ flux was highly variable, ranging from an average flux of 7 mg CH₄ m^?2 d^?1 in the deep-water zone to 154 mg CH₄ m^?2 d^?1 along the lake margin. Diffusive flux accounted for ~8.5 % of mean net CH₄ flux, with plant-mediated and ebullitive flux accounting for the balance of the total net flux. Methane content of emitted ebullition was low (mean ± SD 10 ± 17 %) compared to previous studies from boreal lakes and wetlands. Isotopic composition of net CH₄ emissions varied widely throughout the system, with δ¹³C-CH₄ ranging from ?66.2 to ?55.5 ‰, and δ²H-CH₄ ranging from ?345 to ?258 ‰. Carbon isotope composition of CH₄ in ebullitive flux showed wider variation compared to net flux, ranging from ?69.2 to ?49.2 ‰. Dissolved CH₄ concentrations were highest in the sediment and decreased up the water column. Higher concentrations of CH₄ in the hypoxic deep water coincided with decreasing dissolved O₂ concentrations, while methanotrophic oxidation dominated in the epilimnion based upon decreasing concentrations and increasing values of δ¹³C-CH₄ and δ²H-CH₄. The most depleted ¹³C- and ²H-isotopic values were observed in profundal bottom waters and in subsurface profundal sediments. Based upon paired δ¹³C and δ²H observations of CH₄, acetate fermentation was likely the dominant production pathway throughout the system. However, isotopic ratios of CH₄ in deeper sediments were consistent with mixing/transition between CH₄ production pathways, indicating a higher contribution of the CO₂ reduction pathway. The large spatial variability in fluxes of CH₄ and in isotopic composition of CH₄ throughout a single lake system indicates that the underlying mechanisms controlling CH₄ cycling (production, consumption and transport) are spatially heterogeneous. Net flux along the lake margin dominated whole-lake flux, suggesting the nearshore littoral area dominates CH₄ emissions in these systems. Future studies of whole-lake CH₄ budgets should consider this significant spatial heterogeneity.