The effect of noncollinearity of 15N-1H dipolar and 15N CSA tensors and rotational anisotropy on 15N relaxation, CSA/dipolar cross correlation, and TROSY

Current approaches to 15N relaxation in proteins assume that the 15N-1H dipolar and 15N CSA tensors are collinear. We show theoretically that, when there is significant anisotropy of molecular rotation, different orientations of the two tensors, experimentally observed in proteins, nucleic acids, and small peptides, will result in differences in site- specific correlation functions and spectral densities. The standard treatments of the rates of longitudinal and transverse relaxation of amide 15N nuclei, of the 15N CSA/15N-1H dipolar cross correlation, and of the TROSY experiment are extended to account for the effect of noncollinearity of the 15N-1H dipolar and 15N CSA (chemical shift anisotropy) tensors. This effect, proportional to the degree of anisotropy of the overall motion, (D/D–1), is sensitive to the relative orientation of the two tensors and to the orientation of the peptide plane with respect to the diffusion coordinate frame. The effect is negligible at small degrees of anisotropy, but is predicted to become significant for D/D1.5, and at high magnetic fields. The effect of noncollinearity of 15N CSA and 15N-1H dipolar interaction is sensitive to both gross (hydrodynamic) properties and atomic-level details of protein structure. Incorporation of this effect into relaxation data analysis is likely to improve both precision and accuracy of the derived characteristics of protein dynamics, especially at high magnetic fields and for molecules with a high degree of anisotropy of the overall motion. The effect will also make TROSY efficiency dependent on local orientation in moderately anisotropic systems.     

Cross-correlation between a carbonyl C' chemical shift anisotropy and a long-range dipolar C'HA coupling in proteins using symmetrical reconversion

A new sequence is described to measure the cross-correlation rates between the chemical shift anisotropy of the carbonyl carbon-13 nucleus and the dipole-dipole interaction between this carbonyl and the alpha-proton in proteins. The sequence is based on the symmetrical reconversion principle and is insensitive to experimental errors and to violations of the secular approximation. The cross-correlation rate depends on the backbone angle . The advantages and limitations of the sequence are discussed.     

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

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

Rapid and accurate calculation of protein 1H, 13C and 15N chemical shifts

A computer program (SHIFTX) is described which rapidly and accurately calculates the diamagnetic ¹H, ¹³C and ¹⁵N chemical shifts of both backbone and sidechain atoms in proteins. The program uses a hybrid predictive approach that employs pre-calculated, empirically derived chemical shift hypersurfaces in combination with classical or semi-classical equations (for ring current, electric field, hydrogen bond and solvent effects) to calculate ¹H, ¹³C and ¹⁵N chemical shifts from atomic coordinates. The chemical shift hypersurfaces capture dihedral angle, sidechain orientation, secondary structure and nearest neighbor effects that cannot easily be translated to analytical formulae or predicted via classical means. The chemical shift hypersurfaces were generated using a database of IUPAC-referenced protein chemical shifts – RefDB (Zhang et al., 2003), and a corresponding set of high resolution (<2.1 Å) X-ray structures. Data mining techniques were used to extract the largest pairwise contributors (from a list of 20 derived geometric, sequential and structural parameters) to generate the necessary hypersurfaces. SHIFTX is rapid (< 1 CPU second for a complete shift calculation of 100 residues) and accurate. Overall, the program was able to attain a correlation coefficient (r) between observed and calculated shifts of 0.911 (¹H), 0.980 (¹³C), 0.996 (¹³C), 0.863 (¹³CO), 0.909 (¹⁵N), 0.741 (¹HN), and 0.907 (sidechain ¹H) with RMS errors of 0.23, 0.98, 1.10, 1.16, 2.43, 0.49, and 0.30 ppm, respectively on test data sets. We further show that the agreement between observed and SHIFTX calculated chemical shifts can be an extremely sensitive measure of the quality of protein structures. Our results suggest that if NMR-derived structures could be refined using heteronuclear chemical shifts calculated by SHIFTX, their precision could approach that of the highest resolution X-ray structures. SHIFTX is freely available as a web server at http://redpoll.pharmacy.ualberta.ca.     

Overall structure and sugar dynamics of a DNA dodecamer from homo- and heteronuclear dipolar couplings and 31P chemical shift anisotropy

The solution structure of d(CGCGAATTCGCG)₂ has been determined on the basis of an exceptionally large set of residual dipolar couplings. In addition to the heteronuclear ¹³C-¹H and ¹⁵N-¹H and qualitative homonuclear ¹H-¹H dipolar couplings, previously measured in bicelle medium, more than 300 quantitative ¹H-¹H and 22 ³¹P-¹H dipolar restraints were obtained in liquid crystalline Pf1 medium, and 22 ³¹P chemical shift anisotropy restraints. High quality DNA structures can be obtained solely on the basis of these new restraints, and these structures are in close agreement with those calculated previously on the basis of ¹³C-¹H and ¹⁵N-¹H dipolar couplings. In the newly calculated structures, ³¹P-¹H dipolar and ³Jsub H3 P sub couplings and ³¹P CSA data restrain the phosphodiester backbone torsion angles. The final structure represents a quite regular B-form helix with a modest bending of 10°, which is essentially independent of whether or not electrostatic terms are used in the calculation. Combined, the number of homo- and heteronuclear dipolar couplings significantly exceeds the number of degrees of freedom in the system. Results indicate that the dipolar coupling data cannot be fit by a single structure, but are compatible with the presence of rapid equilibria between C2-endo and C3-endo deoxyribose puckers (sugar switching). The C2-H2/H2 dipolar couplings in B-form DNA are particularly sensitive to sugar pucker and yield the largest discrepancies when fit to a single structure. To resolve these discrepancies, we suggest a simplified dipolar coupling analysis that yields N/S equilibria for the ribose sugar puckers, which are in good agreement with previous analyses of NMR J_HH couplings, with a population of the minor C3-endo form higher for pyrimidines than for purines.     

1H, 13C and 15N chemical shift referencing in biomolecular NMR

Summary A considerable degree of variability exists in the way that ¹H, ¹³C and ¹⁵N chemical shifts are reported and referenced for biomolecules. In this article we explore some of the reasons for this situation and propose guidelines for future chemical shift referencing and for conversion from many common ¹H, ¹³C and ¹⁵N chemical shift standards, now used in biomolecular NMR, to those proposed here.Abbreviations TMS tetramethylsilane - TSP 3-(trimethylsilyl)-propionate, sodium salt - DSS 2,2-dimethyl-2-silapentane-5-sulfonate, sodium salt - TFE 2,2,2-trifluoroethanol - DMSO dimethyl sulfoxide     

13C solid-state NMR chemical shift anisotropy analysis of the anomeric carbon in carbohydrates

(13)C NMR solid-state structural analysis of the anomeric center in carbohydrates was performed on six monosaccharides: glucose (Glc), mannose (Man), galactose (Gal), galactosamine hydrochloride (GalN), glucosamine hydrochloride (GlcN), and N-acetyl-glucosamine (GlcNAc). In the 1D (13)C cross-polarization/magic-angle spinning (CP/MAS) spectrum, the anomeric center C-1 of these carbohydrates revealed two well resolved resonances shifted by 3-5ppm, which were readily assigned to the anomeric alpha and beta forms. From this experiment, we also extracted the (13)C chemical shift anisotropy (CSA) tensor elements of the two forms from their spinning sideband intensities, respectively. It was found out that the chemical shift tensor for the alpha anomer was more axially symmetrical than that of the beta form. A strong linear correlation was obtained when the ratio of the axial asymmetry of the (13)C chemical shift tensors of the two anomeric forms was plotted in a semilogarithmic plot against the relative population of the two anomers. Finally, we applied REDOR spectroscopy to discern whether or not there were any differences in the sugar ring conformation between the anomers. Identical two-bond distances of 2.57A (2.48A) were deduced for both the alpha and beta forms in GlcNAc (GlcN), suggesting that the two anomers have essentially identical sugar ring scaffolds in these sugars. In light of these REDOR distance measurements and the strong correlation observed between the ratio of the axial asymmetry parameters of the (13)C chemical shift tensors and the relative population between the two anomeric forms, we concluded that the anomeric effect arises principally from interaction of the electron charge clouds between the C-1-O-5 and the C-1-O-1 bonds in these monosaccharides.     

Direct structure refinement of high molecular weight proteins against residual dipolar couplings and carbonyl chemical shift changes upon alignment: an application to maltose binding protein

The global fold of maltose binding protein in complex with -cyclodextrin has been determined using a CNS-based torsion angle molecular dynamics protocol involving direct refinement against dipolar couplings and carbonyl chemical shift changes that occur upon alignment. The shift changes have been included as structural restraints using a new module, CANI, that has been incorporated into CNS. Force constants and timesteps have been determined that are particularly effective in structure refinement applications involving high molecular weight proteins with small to moderate numbers of NOE restraints. Solution structures of the N- and C-domains of MBP calculated with this new protocol are within 2 Å of the X-ray conformation.     

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     

蛋白质/核酸相互作用研究方法进展

蛋白质和核酸是构成生命体最为重要的两类生物大分子,蛋白质与核酸的相互作用是分子生物学研究的中心问题之一,它是许多生命活动的重要组成部分。研究蛋白质/核酸相互作用近期采用的新技术有:核酸适体技术、生物信息学方法、蛋白质芯片技术以及纳米技术等。本文就这些新的研究方法进行综述。     

Measurement of methyl 13C-1H cross-correlation in uniformly 13C-, 15N-, labeled proteins

An understanding of side chain motions in protein is of great interest since side chains often play an important role in protein folding and intermolecular interactions. A novel method for measuring the dynamics of methyl groups in uniformly ¹³C-, ¹⁵N-labeled proteins has been developed by our group. The method relies on the difference in peak intensities of ¹³C quartet components of methyl groups, in a spectrum recording the free evolution of ¹³C under proton coupling in a constant-time period. Cross-correlated relaxation rates between ¹³C-¹H dipoles can be easily measured from the intensities of the multiplet components. The degree of the methyl restrictions (S ²) can be estimated from the cross-correlated relaxation rate. The method is demonstrated on a sample of human fatty acid binding protein in the absence of fatty acid. We obtained relaxation data for 33 out of 46 residues having methyl groups in apo-IFABP. It has been found that the magnitude of the CSA tensor of spin ¹³C in a methyl group could be estimated from the intensities of the ¹³C multiplet components.     

Empirical correlation between protein backbone <Superscript>15</Superscript>N and <Superscript>13</Superscript>C secondary chemical shifts and its application to nitrogen chemical shift re-referencing

The linear analysis of chemical shifts (LACS) has provided a robust method for identifying and correcting ¹³C chemical shift referencing problems in data from protein NMR spectroscopy. Unlike other approaches, LACS does not require prior knowledge of the three-dimensional structure or inference of the secondary structure of the protein. It also does not require extensive assignment of the NMR data. We report here a way of extending the LACS approach to ¹⁵N NMR data from proteins, so as to enable the detection and correction of inconsistencies in chemical shift referencing for this nucleus. The approach is based on our finding that the secondary ¹⁵N chemical shift of the backbone nitrogen atom of residue i is strongly correlated with the secondary chemical shift difference (experimental minus random coil) between the alpha and beta carbons of residue i − 1. Thus once alpha and beta ¹³C chemical shifts are available (their difference is referencing error-free), the ¹⁵N referencing can be validated, and an appropriate offset correction can be derived. This approach can be implemented prior to a structure determination and can be used to analyze potential referencing problems in database data not associated with three-dimensional structure. Application of the LACS algorithm to the current BMRB protein chemical shift database, revealed that nearly 35% of the BMRB entries have δ ¹⁵N values mis-referenced by over 0.7 ppm and over 25% of them have δ ¹H^N values mis-referenced by over 0.12 ppm. One implication of the findings reported here is that a backbone ¹⁵N chemical shift provides a better indicator of the conformation of the preceding residue than of the residue itself. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Differentiation between unlabeled and very-low-level fully 15N,13C-labeled nucleotides for resonance assignments in nucleic acids

Based on different characteristics between unlabeled and fully ¹⁵N,¹³C-labeled nucleotides, we develop a method for unambiguous resonance assignments in nucleic acids following site-specific fully ¹⁵N,¹³C isotope incorporation at very low levels¹. The J-couplings between heteronuclei provide for distinction between the NMR signals of the fully labeled nucleotides and those of the natural abundance nucleotides. The method is demonstrated for DNA oligonucleotides², in the dimeric G-quadruplex [d(GGGTTCAGG)]₂and in the 22-nucleotide human telomeric fragment d[AG₃(TTAG₃)₃]. We expect this approach to be useful for selective monitoring of important functional domains and of their interactions in large nucleic acids.     

1H, 13C and 15N assignments and chemical shift-derived secondary structure of intestinal fatty acid-binding protein

Summary Sequence-specific ¹H, ¹³C and ¹⁵N resonance assignments have been established for rat intestinal fatty acid-binding protein complexed with palmitate (15.4 kDa) at pH 7.2 and 37°C. The resonance assignment strategy involved the concerted use of seven 3D triple-resonance expriments (CC-TOCSY, HCCH-TOCSY, HNCO, HNCA, ¹⁵N-TOCSY-HMQC, HCACO and HCA(CO)N). A central feature of this strategy was the concurrent assignment of both backbone and side-chain aliphatic atoms, which was critical for overcoming ambiguities in the assignment process. The CC-TOCSY experiment provided the unambiguous links between the side-chain spin systems observed in HCCH-TOCSY and the backbone correlations observed in the other experiments. Assignments were established for 124 of the 131 residues, although 6 of the 124 had missing amide ¹H resonances, presumably due to rapid exchange with solvent under these experimental conditions. The assignment database was used to determine the solution secondary structure of the complex, based on chemical shift indices for the ¹H, ¹³C, ¹³C and ¹³CO atoms. Overall, the secondary structure agreed well with that determined by X-ray crystallography [Sacchettini et al. (1989) J. Mol. Biol., 208, 327–339], although minor differences were observed at the edges of secondary structure elements.