Determination of 13Cα relaxation times in uniformly 13C/15N-enriched proteins

Summary Relaxation times of ¹³C carbons of uniformly ¹³C/¹⁵N-enriched probes have been investigated. The relaxation behaviour was analyzed in terms of a multispin system. Pulse sequences for the determination of T₁, T₂ and the heteronuclear NOE of ¹³C in uniformly ¹³C/¹⁵N-enriched ribonuclease T1 are presented. The experiments performed in order to obtain T₁ and the heteronuclear NOE were similar to those of the corresponding ¹⁵N experiments published previously. The determination of T₂ for the C-carbon in a completely labeled protein is more complicated, since the magnetization transfer during the T₂ evolution period owing to the scalar coupling of C–C must be suppressed. Various different pulse sequences for the T₂ evolution period were simulated in order to optimize the bandwidth for which reliable T₂ relaxation times can be obtained. A proof for the quality of these pulse sequences is given by fitting the intensity decay of individual ¹H–¹³C cross peaks, in a series of (¹H, ¹³C)-ct-HSQC spectra with a modified CPMG sequence as well as a T_1p sequence for the transverse relaxation time, to a single exponential using a simplex algorithm.     

Automated prediction of 15N, 13Cα, 13Cβ and 13C′ chemical shifts in proteins using a density functional database

A database of peptide chemical shifts, computed at the density functional level, has been used to develop an algorithm for prediction of ¹⁵N and ¹³C shifts in proteins from their structure; the method is incorporated into a program called SHIFTS (version 4.0). The database was built from the calculated chemical shift patterns of 1335 peptides whose backbone torsion angles are limited to areas of the Ramachandran map around helical and sheet configurations. For each tripeptide in these regions of regular secondary structure (which constitute about 40% of residues in globular proteins) SHIFTS also consults the database for information about sidechain torsion angle effects for the residue of interest and for the preceding residue, and estimates hydrogen bonding effects through an empirical formula that is also based on density functional calculations on peptides. The program optionally searches for alternate side-chain torsion angles that could significantly improve agreement between calculated and observed shifts. The application of the program on 20 proteins shows good consistency with experimental data, with correlation coefficients of 0.92, 0.98, 0.99 and 0.90 and r.m.s. deviations of 1.94, 0.97, 1.05, and 1.08 ppm for ¹⁵N, ¹³C, ¹³C and ¹³C, respectively. Reference shifts fit to protein data are in good agreement with `random-coil' values derived from experimental measurements on peptides. This prediction algorithm should be helpful in NMR assignment, crystal and solution structure comparison, and structure refinement.     

MQ-HCN-based pulse sequences for the measurement of 13C1′-1H1′, 13C1′-15N, 1H1′-15N, 13C1′-13C2′, 1H1′-13C2′, 13C6/8-1H6/8, 13C6/8-15N, 1H6/8-15N, 13C6-13C5, 1H6-13C5 dipolar couplings in 13C, 15N-labeled DNA (and RNA)

A suite of multiple quantum (MQ) HCN-based pulse sequences has been developed for the purpose of collecting dipolar coupling data in labeled nucleic acids. All the pulse sequences are based on the robust MQ-HCN experiment which has been utilized for assignment purposes in labeled nucleic acids for a number of years and provides much-needed resolution for the dipolar coupling measurements. We have attempted to collect multiple couplings centered on the 13C1' and 13C6/8 positions. Six pulse sequences are described, one each for measurement of one-bond 13C1'-1H1' and 13C6/8-1H6/8 couplings, one for measurement of one-bond 13C1'-15N and two-bond 1H1'-15N couplings, one for measurement of one-bond 13C6/8-15N and two-bond 1H6/8-15N couplings, one for measurement of one-bond 13C1'- 13C2' and two-bond 1H1'-13C2' couplings, and one for measurement of one-bond 13C6-13C5 and two-bond 1H6-13C5 couplings in the bases of C and T. These sequences are demonstrated for a labeled 18 bp DNA duplex in a 47 kDa ternary complex of DNA, CBFbeta, and the CBFalpha Runt domain, thus clearly demonstrating the robustness of the pulse sequences even for a very large complex.     

Identification of Zinc-ligated Cysteine Residues Based on 13Cα and 13Cβ Chemical Shift Data

Although a significant number of proteins include bound metals as part of their structure, the identification of amino acid residues coordinated to non-paramagnetic metals by NMR remains a challenge. Metal ligands can stabilize the native structure and/or play critical catalytic roles in the underlying biochemistry. An atom’s chemical shift is exquisitely sensitive to its electronic environment. Chemical shift data can provide valuable insights into structural features, including metal ligation. In this study, we demonstrate that overlapped ¹³Cβ chemical shift distributions of Zn-ligated and non-metal-ligated cysteine residues are largely resolved by the inclusion of the corresponding ¹³Cα chemical shift information, together with secondary structural information. We demonstrate this with a bivariate distribution plot, and statistically with a multivariate analysis of variance (MANOVA) and hierarchical logistic regression analysis. Using 287 ¹³Cα/¹³Cβ shift pairs from 79 proteins with known three-dimensional structures, including 86 ¹³Cα and¹³Cβ shifts for 43 Zn-ligated cysteine residues, along with corresponding oxidation state and secondary structure information, we have built a logistic regression model that distinguishes between oxidized cystines, reduced (non-metal ligated) cysteines, and Zn-ligated cysteines. Classifying cysteines/cystines with a statisical model incorporating all three phenomena resulted in a predictor of Zn ligation with a recall, precision and F-measure of 83.7%, and an accuracy of 95.1%. This model was applied in the analysis of Bacillus subtilis IscU, a protein involved in iron–sulfur cluster assembly. The model predicts that all three cysteines of IscU are metal ligands. We confirmed these results by (i) examining the effect of metal chelation on the NMR spectrum of IscU, and (ii) inductively coupled plasma mass spectrometry analysis. To gain further insight into the frequency of occurrence of non-cysteine Zn ligands, we analyzed the Protein Data Bank and found that 78% of the Zn ligands are histidine and cysteine (with nearly identical frequencies), and 18% are acidic residues aspartate and glutamate. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

13C–13C NOESY spectra of a 480 kDa protein: solution NMR of ferritin

Molecular size has limited solution NMR analyses of proteins. We report ¹³C–¹³C NOESY experiments on a 480 kDa protein, the multi-subunit ferritin nanocage with gated pores. By exploiting ¹³C-resonance-specific chemical shifts and spin diffusion effects, we identified 75% of the amino acids, with intraresidue C–C connectivities between nuclei separated by 1–4 bonds. These results show the potential of ¹³C–¹³C NOESY for solution studies of molecular assemblies >100 kDa.     

Measurement of 13Cα-13CO cross-relaxation rates in 15N-/13C-labelled proteins

Summary Internal motions play an important role in the biological function of proteins and NMR relaxation studies may characterize them over a wide range of frequencies. An experimental pulse scheme is proposed for the measurement of the ¹³CO-¹³C cross-relaxation rate. For sensitivity reasons, this measurement is performed in an indirect manner through several coherence transfer steps, which should thus be calibrated independently. Contributions of other relaxation pathways can be eliminated by the determination of the initial slope of the buildup curve. The cross-relaxation rates have been determined on a ¹⁵N-/¹³C-labelled 116-residue protein and the significant variations along the sequence have been interpreted as evidence of an increased amount of fast local motion.     

13Cα and 13Cβ chemical shifts as a tool to delineate β-hairpin structures in peptides

Unravelling the factors that contribute to the formation and the stability of -sheet structure in peptides is a subject of great current interest. A -hairpin, the smallest -sheet motif, consists of two antiparallel hydrogen-bonded -strands linked by a loop region. We have performed a statistical analysis on protein -hairpins showing that the most abundant types of -hairpins, 2:2, 3:5 and 4:4, have characteristic patterns of ¹³C and ¹³C conformational shifts, as expected on the basis of their and angles. This fact strongly supports the potential value of ¹³C and ¹³C conformational shifts as a means to identify -hairpin motifs in peptides. Their usefulness was confirmed by analysing the patterns of ¹³C and ¹³C conformational shifts in 13 short peptides, 10–15 residues long, that adopt -hairpin structures in aqueous solution. Furthermore, we have investigated their potential as a method to quantify -hairpin populations in peptides.     

13C代谢流量分析发展30年

~(13)C代谢流量分析(~(13)C metabolic flux analysis,~(13)C-MFA),是通过标记实验分析蛋白氨基酸或胞内代谢物同位素标记异构体的分布情况,从而准确定量胞内反应速率。该技术在系统理解细胞代谢特性、指导代谢工程改造和揭示病理生理学等方面起着重要作用,引起研究者的广泛重视。文中重点综述了代谢流分析30年的发展历程,尤其在工业生物技术和生物医药领域的应用,并对未来的发展方向进行展望。     

Out-and-back 13C–13C scalar transfers in protein resonance assignment by proton-detected solid-state NMR under ultra-fast MAS

We present here ¹H-detected triple-resonance H/N/C experiments that incorporate CO–CA and CA–CB out-and-back scalar-transfer blocks optimized for robust resonance assignment in biosolids under ultra-fast magic-angle spinning (MAS). The first experiment, (H)(CO)CA(CO)NH, yields ¹H-detected inter-residue correlations, in which we record the chemical shifts of the CA spins in the first indirect dimension while during the scalar-transfer delays the coherences are present only on the longer-lived CO spins. The second experiment, (H)(CA)CB(CA)NH, correlates the side-chain CB chemical shifts with the NH of the same residue. These high sensitivity experiments are demonstrated on both fully-protonated and 100 %-H^N back-protonated perdeuterated microcrystalline samples of Acinetobacter phage 205 (AP205) capsids at 60 kHz MAS.     

TROSY-based HNCO pulse sequences for the measurement of 1HN-15N, 15N-13CO, 1HN-13CO, 13CO-13Cα and 1HN-13Cα dipolar couplings in 15N, 13C, 2H-labeled proteins

HNCO-based 3D pulse schemes are presented for measuring ¹HN-¹⁵N,¹⁵N-¹³CO, ¹HN-¹³CO,¹³CO-¹³C and ¹HN-¹³C dipolar couplings in ¹⁵N,¹³C,²-labeled proteins. The experiments are based on recently developed TROSY methodology for improving spectral resolution and sensitivity. Data sets recorded on a complex of Val, Leu, Ile (1 only) methyl protonated ¹⁵N,¹³C,²H-labeled maltose binding protein and -cyclodextrin as well as ¹⁵N,¹³C,²H-labeled human carbonic anhydrase II demonstrate that precise dipolar couplings can be obtained on proteins in the 30–40 kDa molecular weight range. These couplings will serve as powerful restraints for obtaining global folds of highly deuterated proteins.     

Quantitative J correlation methods for the accurate measurement of 13C′-13Cαdipolar couplings in proteins

Methods are described for the precise and accurate measurement of one-bond dipolar (13)C'-(13)C(alpha) couplings in weakly aligned proteins. The experiments are based on the principle of quantitative J correlation, where (1)J(C'C(alpha)) (or (1)J(C'C(alpha)) + 1D(C'C(alpha)) is measured from the relative intensity of two interleaved 3D TROSY-HN(CO)CA or 3DTROSY-HNCO spectra recorded with dephasing intervals of zero (reference spectrum) and approximately 3/(2(1)J(C'C(alpha)) (attenuated spectrum). In analogy to other quantitative J correlation techniques, the random error in the measured (1)J(C'C(alpha)) value is inversely proportional to the signal-to-noise ratio in the reference spectrum. It is shown that for weakly aligned proteins, with the magnitude of the alignment tensor of D(a)(NH) < or = 10-15 Hz, the systematic errors are typically negligible. The methods are demonstrated for the third IgG-binding domain of protein G (GB3) and a-synuclein in complex with a detergent micelle, where errors in (1)D(C'C(alpha)) of less than 0.1 Hz and ca. 0.2 Hz,respectively, are estimated. Remarkably, the dipolar couplings determined for GB3 are in even better agreement with the recently refined 1.1-angstroms X-ray structure than the input (13)C'-(13)C(alpha) couplings used for the refinement.     

A 13C NMR spectrometric method for the determination of intramolecular δ13C values in fructose from plant sucrose samples

Recent developments in (13) C NMR spectrometry have allowed the determination of intramolecular (13) C/(12) C ratios with high precision. However, the analysis of carbohydrates requires their derivatization to constrain the anomeric carbon. Fructose has proved to be particularly problematic because of a byproduct occurring during derivatization and the complexity of the NMR spectrum of the derivative. Here, we describe a method to determine the intramolecular (13) C/(12) C ratios in fructose by (13) C NMR analysis of the acetyl-isopropylidene derivative. We have applied this method to measure the intramolecular (13) C/(12) C distribution in the fructosyl moiety of sucrose and have compared this with that in the glucosyl moiety. Three prominent features stand out. First, in sucrose from both C(3) and C(4) plants, the C-1 and C-2 positions of the glucosyl and fructosyl moieties are markedly different. Second, these positions in C(3) and C(4) plants show a similar profile. Third, the glucosyl and fructosyl moieties of sucrose from Crassulacean acid metabolism (CAM) metabolism have a different profile. These contrasting values can be interpreted as a result of the isotopic selectivity of enzymes that break or make covalent bonds in glucose metabolism, whereas the distinctive (13) C pattern in CAM sucrose probably indicates a substantial contribution of gluconeogenesis to glucose synthesis.     

Selective diagonal-free 13C,13C-edited aliphatic–aromatic NOESY experiment with non-uniform sampling

A band-selective aromatic–aliphatic C,C-edited four-dimensional NOESY experiment is proposed here. Its key advantage is the absence of auto-correlation signals which makes it very attractive for joint use with non-uniform sampling. It is demonstrated here that the sensitivity of the experiment is not significantly affected by utilization of selective pulses (for either aromatic-13C or aliphatic-13C spins). The method was applied to the sample of E32Q mutant of human S100A1 protein, a homodimer of total molecular mass ~20 kDa. High-resolution 4D spectra were obtained from ~1.5 % of sampling points required conventionally. It is shown that superior resolution facilitates unambiguous assignment of observed aliphatic–aromatic cross-peaks. Additionally, the addition of aliphatic-13C dimension enables to resolve peaks with degenerated aliphatic ¹H chemical shifts. All observed cross-peaks were validated against previously determined 3D structure of E32Q mutant of S100A1 protein (PDB 2LHL). The increased reliability of structural constraints obtained from the proposed high-resolution 4D 13C(ali),13C(aro)-edited NOESY can be exploited in the automated protocols of structure determination of proteins.     

Synthesis and SAR of C12–C13-oxazoline derivatives of epothilone A

The SAR of a series of new epothilone A derivatives with a 2-substituted-1,3-oxazoline moiety trans-fused to the C12–C13 bond of the deoxy macrocycle have been investigated with regard to tubulin polymerization induction and cancer cell growth inhibition. Significant differences in antiproliferative activity were observed between different analogs, depending on the nature of the substituent at the 2-position of the oxazoline ring. The most potent compounds showed comparable activity with the natural product epothilone A. Modeling studies provide a preliminary rationale for the observed SAR.     

Study of protein dynamics in solution by measurement of 13Cα-13CO NOE and 13CO longitudinal relaxation

Summary ¹³C-¹³CO homonuclear NOE and ¹³CO T₁ relaxation were measured for a 20 kDa protein using tripleresonance pulse sequences. The experiments were sufficiently sensitive to obtain statistically significant differences in relaxation parameters over the molecule. The ¹³C-¹³CO cross-relaxation rate, obtained from these data, is directly proportional to an order parameter describing local motion and it is largely independent of the local correlation time. It is therefore a relatively straightforward observable for the identification of local dynamics.     

δ13C values of some succulent plants from Madagascar

Summary ¹³C values were determined in 20 succulents from Madagascar. The values were indicative of Crassulacean Acid Metabolism in 10 species of the Didiereaceae, 4 species of the Euphorbiaceae, 2 species of the Crassulaceae and 1 species of the Cucurbitaceae. The Didiereaceae and Euphorbiaceae studied are major components of a high biomass xerophytic flora in the semi-arid southwest and south of Madagascar. Three species of the Euphorbiaceae with succulent stems and non-succulent leaves, which were cultivated outdoors in the Tananarive Botanic Garden, showed C₃ like ¹³C values for both leaves and stems. ¹³C values of leaf and stem material from a similar species, collected in the south of Madagascar, indicated Crassulacean Acid Metabolism.Abbreviations CAM Crassulacean Acid Metabolism     

Does the 13C of foliage‐respired CO2 and biochemical pools reflect the 13C of recently assimilated carbon?

Isotopic labelling experiments were conducted to assess relationships among ¹³C of recently assimilated carbon ( δC _A), foliage respiration ( δC _F), soluble carbohydrate ( δC _SC), leaf waxes ( δC _LW) and bulk organic matter ( δC _OM). Slash pine, sweetgum and maize were grown under ¹³C depleted CO₂ to label biomass and then placed under ambient conditions to monitor the loss of label. In pine and sweetgum, δC _F of labelled plants (∼−44 and −35‰, respectively) rapidly approached control values but remained depleted by ∼4–6‰ after 3–4 months. For these tree species, no or minimal label was lost from δC _SC, δC _LW and δC _OM during the observation periods. δC _F and δC _SC of labelled maize plants rapidly changed and were indistinguishable from controls after 1 month, while δC _LW and δC _OM more slowly approached control values and remained depleted by 2–6‰. Changes in δC _F in slash pine and sweetgum fit a two-pool exponential model, with the fast turnover metabolic pool (∼3–4 d half-life) constituting only 1–2% of the total. In maize, change in δC _F fits a single pool model with a half-life of 6.4 d. The ¹³C of foliage respiration and biochemical pools reflect temporally integrated values of δC _A, with change in isotopic composition dampened by the size of metabolic carbon reserves and turnover rates.     

Cα and Cβ Carbon-13 Chemical Shifts in Proteins From an Empirical Database

We have constructed an extensive database of 13C Cα and Cβ chemical shifts in proteins of solution, for proteins of which a high-resolution crystal structure exists, and for which the crystal structure has been shown to be essentially identical to the solution structure. There is no systematic effect of temperature, reference compound, or pH on reported shifts, but there appear to be differences in reported shifts arising from referencing differences of up to 4.2 ppm. The major factor affecting chemical shifts is the backbone geometry, which causes differences of ca. 4 ppm between typical α- helix and β-sheet geometries for Cα, and of ca. 2 ppm for Cβ. The side-chain dihedral angle χ1 has an effect of up to 0.5 ppm on the Cα shift, particularly for amino acids with branched side-chains at Cβ. Hydrogen bonding to main-chain atoms has an effect of up to 0.9 ppm, which depends on the main- chain conformation. The sequence of the protein and ring-current shifts from aromatic rings have an insignificant effect (except for residues following proline). There are significant differences between different amino acid types in the backbone geometry dependence; the amino acids can be grouped together into five different groups with different φ,ψ shielding surfaces. The overall fit of individual residues to a single non-residue-specific surface, incorporating the effects of hydrogen bonding and χ1 angle, is 0.96 ppm for both Cα and Cβ. The results from this study are broadly similar to those from ab initio studies, but there are some differences which could merit further attention.     

13C NMR在马利筋族植物C21甾类成分结构研究中的应用

本文介绍了＾１３ＣＮＭＲ在以列筋族植物Ｃ２１甾类成分甙元类型工基的确定、糖的种类、连接顺序的判断等方面的应用。