Determination of the glycosidic bond angle chi in RNA from cross-correlated relaxation of CH dipolar coupling and N chemical shift anisotropy

A new heteronuclear NMR pulse sequence, the quantitative Gamma(HCN) experiment, for the determination of the glycosidic torsion angle chi in (13)C,(15)N-labeled oligonucleotides is described. The Gamma(HCN) experiment allows measurement of CH dipole-dipole, N chemical shift anisotropy cross-correlated relaxation rates (Gamma(C1'H1',N1)(DD,CSA) and Gamma(C2'H2',N9)(DD,CSA) for pyrimidines Gamma(C1'H1'N9)(DD,CSA) and Gamma(C2'H2',N9)(DD,CSA) for purines). A nucleotide-specific parametrization for the dependence of these Gamma-rates on chi based on (15)N chemical shift tensors determined by solid-state NMR experiments on mononucleosides (Stueber, D.; Grant, D. M. J. Am. Chem. Soc. 2002, 124, 10539-10551) is presented. For a 14-mer and a 30-mer RNA of known structures, it is found that the Gamma(HCN) experiment offers a very sensitive parameter for changes in the angle chi and allows restraining of chi with an accuracy of around 10 degrees for residues which do not undergo conformational averaging. Therefore, the Gamma(HCN) experiment can be used for the determination of chi in addition to data derived from (3)J(C,H)-coupling constants. As shown for the 30-mer RNA, the derived torsion angle information can be incorporated as additional restraint, improving RNA structure calculations.     

Chemical shift anisotropy tensors of carbonyl, nitrogen, and amide proton nuclei in proteins through cross-correlated relaxation in NMR spectroscopy

The principal components and orientations of the chemical shift anisotropy (CSA) tensors of the carbonyl (C'), nitrogen (N), and amide proton (H(N)) nuclei of 64 distinct amide bonds in human ubiquitin have been determined in isotropic solution by a set of 14 complementary auto- and cross-correlated relaxation rates involving the CSA interactions of the nuclei of interest and several dipole-dipole (DD) interactions. The CSA parameters thus obtained depend to some degree on the models used for local motions. Three cases have been considered: restricted isotropic diffusion, three-dimensional Gaussian axial fluctuations (3D-GAF), and independent out-of-plane motions of the NH(N) vectors with respect to the peptide planes.     

31P chemical shift tensors for canonical and non-canonical conformations of nucleic acids: a DFT study and NMR implications

31P chemical shift anisotropy (CSA) tensors have been calculated for a set of selected DNA and RNA backbone conformations using density functional theory. The set includes canonical A-RNA, A-DNA, BI-DNA, BII-DNA, ZI-DNA, and ZII-DNA as well as four A-RNA-type, seven non-A-RNA-type, and three non-canonical DNA conformations. Hexahydrated dimethyl phosphate has been employed as a model. The 31P chemical shift tensors obtained are discussed in terms of similarities in the behavior observed for gauche-gauche (gg) and gauche-trans (gt) conformations around the P-O bonds. We show that torsion angles alpha and zeta are major determinants of the isotropic chemical shift deltaiso and of the deltaCSA11 component of the traceless chemical shift tensor, which is revealed in separate ranges of both deltaiso and deltaCSA11 for gg- and gt-conformers, respectively. A clear distinction between the two conformation types has not been found for the deltaCSA22 and deltaCSA33 components, which is attributed to their different directional properties. The 31P CSA tensors exhibit considerable variations resulting in large spans of approximately 16 ppm for deltaCSA11 and approximately 22 ppm for deltaCSA22 and deltaCSA33. We examine the consequences of the CSA variations for predicting the chemical shift changes upon partial alignment deltacsa and for the values of CSA order parameters extracted from the analysis of 31P NMR relaxation data. The theoretical 31P CSA tensors as well as the experimental 31P CSA tensor of barium diethyl phosphate (BDEP) are used to calculate deltacsa for two eclipsed orientations of the CSA and molecular alignment tensors. Percentage differences between the CSA order parameters obtained using the theoretical 31P CSA tensors and the experimental 31P CSA tensor of BDEP, respectively, are also determined.     

Eta(z)/kappa: a transverse relaxation optimized spectroscopy NMR experiment measuring longitudinal relaxation interference

NMR spin relaxation experiments provide a powerful tool for the measurement of global and local biomolecular rotational dynamics at subnanosecond time scales. Technical limitations restrict most spin relaxation studies to biomolecules weighing less than 10 kDa, considerably smaller than the average protein molecular weight of 30 kDa. In particular, experiments measuring eta(z), the longitudinal (1)H(N)-(15)N dipole-dipole (DD)/(15)N chemical shift anisotropy (CSA) cross-correlated relaxation rate, are among those least suitable for use with larger biosystems. This is unfortunate because these experiments yield valuable insight into the variability of the (15)N CSA tensor over the polypeptide backbone, and this knowledge is critical to the correct interpretation of most (15)N-NMR backbone relaxation experiments, including R(2) and R(1). In order to remedy this situation, we present a new (1)H(N)-(15)N transverse relaxation optimized spectroscopy experiment measuring eta(z) suitable for applications with larger proteins (up to at least 30 kDa). The presented experiment also yields kappa, the site-specific rate of longitudinal (1)H(N)-(1)H(') DD cross relaxation. We describe the eta(z)/kappa experiment's performance in protonated human ubiquitin at 30.0 degrees C and in protonated calcium-saturated calmodulin/peptide complex at 20.0 degrees C, and demonstrate preliminary experimental results for a deuterated E. coli DnaK ATPase domain construct at 34 degrees C.     

Carbon-13 chemical shift anisotropy in DNA bases from field dependence of solution NMR relaxation rates

Knowledge of (13)C chemical shift anisotropy (CSA) in nucleotide bases is important for the interpretation of solution-state NMR relaxation data in terms of local dynamic properties of DNA and RNA. Accurate knowledge of the CSA becomes particularly important at high magnetic fields, prerequisite for adequate spectral resolution in larger oligonucleotides. Measurement of (13)C relaxation rates of protonated carbons in the bases of the so-called Dickerson dodecamer, d(CGCGAATTCGCG)(2), at 500 and 800 MHz (1)H frequency, together with the previously characterized structure and diffusion tensor yields CSA values for C5 in C, C6 in C and T, C8 in A and G, and C2 in A that are closest to values previously reported on the basis of solid-state FIREMAT NMR measurements, and mostly larger than values obtained by in vacuo DFT calculations. Owing to the noncollinearity of dipolar and CSA interactions, interpretation of the NMR relaxation rates is particularly sensitive to anisotropy of rotational diffusion, and use of isotropic diffusion models can result in considerable errors.     

Chemical shift tensors of protonated base carbons in helical RNA and DNA from NMR relaxation and liquid crystal measurements

Knowledge of (13)C chemical shift anisotropy (CSA) tensors in nucleotide bases is important for interpretation of NMR relaxation data in terms of local dynamic properties of nucleic acids and for analysis of residual chemical shift anisotropy (RCSA) resulting from weak alignment. CSA tensors for protonated nucleic acid base carbons have been derived from measurements on a uniformly (13)C-enriched helical A-form RNA segment and a helical B-form DNA dodecamer at natural (13)C abundance. The magnitudes of the derived CSA principal values are tightly restricted by the magnetic field dependencies of the (13)C transverse relaxation rates, whereas the tensor orientation and asymmetry follow from quantitative measurements of interference between (13)C-{(1)H} dipolar and (13)C CSA relaxation mechanisms. Changes in the chemical shift between the isotropic and aligned states, Deltadelta, complement these measurements and permit cross-validation. The CSA tensors are determined from the experimental Deltadelta values and relaxation rates, under the assumption that the CSA tensor of any specific carbon in a given type of base is independent of the base position in either the RNA or DNA helix. However, the experimental data indicate that for pyrimidine C(6) carbons in A-form RNA the CSA magnitude is considerably larger than in B-form DNA. This result is supported by quantum chemical calculations and is attributed in part to the close proximity between intranucleotide C(6)H and O(5)' atoms in RNA. The magnitudes of the measured CSA tensors, on average, agree better with previous solid-state NMR results obtained on powdered nucleosides than with prior results from quantum chemical calculations on isolated bases, which depend rather strongly on the level of theory at which the calculations are carried out. In contrast, previously computed orientations of the chemical shift tensors agree well with the present experimental results and exhibit less dependence on the level of theory at which the computations are performed.     

Site-specific variations of carbonyl chemical shift anisotropies in proteins

By combining molecular dynamics (MD) simulation with a novel extended density functional theory method, we calculate site-specific carbonyl chemical shift tensors in the SMN Tudor domain. We formulate a simple model for the C' chemical shift anisotropy (CSA) based solely on the isotropic chemical shift. Using this simple chemical shift tensor model and the MD simulation, an accurate prediction of transverse C'/N-H cross-correlated relaxation rates can be obtained.     

New probes of ligand flexibility in drug design: transferred (13)C CSA-dipolar cross-correlated relaxation at natural abundance

Understanding the impact of molecular flexibility remains an important outstanding problem in rational drug design. Toward this end, we present new NMR relaxation methods that describe ligand flexibility at the atomic level. Specifically, we measure natural abundance (13)C cross-correlated relaxation parameters for ligands in rapid exchange between the free and receptor-bound states. The rapid exchange transfers the bound state relaxation parameters to the free state, such that a comparison of relaxation rates in the absence and presence of protein receptor yields site-specific information concerning the bound ligand flexibility. We perform these measurements for aromatic carbons, which are highly prevalent in drug-like molecules and demonstrate significant cross-correlated relaxation between the (13)C-(1)H dipole-dipole (DD) and (13)C chemical shift anisotropy (CSA) relaxation mechanisms. Our use of natural abundance measurements addresses the practical difficulties of obtaining isotope-labeled ligands in pharmaceutical research settings. We demonstrate our methods on a small ligand of the 42 kDa kinase domain of the p38 MAP kinase. We show that exchange-transferred cross-correlated relaxation measurements are not only sensitive probes of bound ligand flexibility but also offer complementary advantages over standard R(1) = 1/T(1) and R(2) = 1/T(2) measurements. The ligand flexibility profiles obtained from the relaxation data can help assess the influence of dynamics on ligand potency or pharmacokinetic properties or both, and thereby include inherent molecular flexibility in drug design.     

Near neighbor approximation in nuclear magnetic resonance

A new way to deal with the excitation by multiple effective RF fields with interference is presented using the coherent averaging theory. It significantly simplifies the calculation of the effect of RF interference that occurs in the excitations by periodic pulses and phase-incremented pulses (PIPs). This approach shows that each neighboring RF field contributes to an excitation profile an offset shift, which is termed the Bloch-Siegert offset shift (BSOS). The BSOS depends not only on the strengths of both RF fields that interfere with each other but also on their relative phase between the two RF fields. Consequently, it can be positive, negative, and zero. In addition, the BSOS is also inversely proportional to the frequency separation of the two RF fields. Therefore, only a few near neighbors need to be taken into account in most cases, resulting in a near neighbor approximation (NNA). The BSOS for two multiband excitation profiles, one by a periodic pulse and the other by a PIP, are calculated using the NNA. The results are in good agreement with the computer simulated ones.     

Residue-specific 13C' CSA tensor principal components for ubiquitin: correlation between tensor components and hydrogen bonding

The residue-specific 13C' CSA tensor principal components, sigma(11), sigma(22), sigma(33), and the tensor orientation defined by the rotation angles beta and gamma have been determined by solution NMR for uniformly labeled ubiquitin partially aligned in four different media. Spurious chemical shift deviations due to solvent effects were corrected with an offset calculated by linear regression of the residual dipolar couplings and chemical shifts at increasing alignment strengths. Analysis of this effect revealed no obvious correlation to solvent exposure. Data obtained in solution from a protein offer a better sampling of 13C' CSA for different amino acid types in a complex heterogeneous environment, thereby allowing for the evaluation of structural variables that would be challenging to achieve by other methods. The 13C' CSA principal components cluster about the average values previously determined, and experimental correlations observed between sigma(11), sigma(22) tensorial components and C'O...H(N) hydrogen bonding are discussed. The inverse association of sigma(11) and sigma(22) exemplify the calculated and solid-state NMR observed effect on the tensor components by hydrogen bonding. We also show that 13C' CSA tensors are sensitive to hydrogen-bond length but not hydrogen-bond angle. This differentiation was previously unavailable. Similarly, hydrogen bonding to the conjugated NH of the same peptide plane has no detectable effect. Importantly, the observed weak correlations signify the presence of confounding influences such as nearest-neighbor effects, side-chain conformation, electrostatics, and other long-range factors to the 13C' CSA tensor. These analyses hold future potential for exploration provided that more accurate data from a larger number of proteins and alignments become available.     

Ultrafast relaxation of zinc protoporphyrin encapsulated within apomyoglobin in buffer solutions

The relaxation dynamics of a zinc protoporphyrin (ZnPP) in THF, KPi buffer, and encapsulated within apomyoglobin (apoMb) was investigated in its excited state using femtosecond fluorescence up-conversion spectroscopy with S2 excitation (lambda(ex) = 430 nm). The S2 --> S1 internal conversion of ZnPP is ultrafast (tau < 100 fs), and the hot S1 ZnPP species are produced promptly after excitation. The relaxation dynamics of ZnPP in THF solution showed a dominant offset component (tau = 2.0 ns), but it disappeared completely when ZnPP formed aggregates in KPi buffer solution. When ZnPP was reconstituted into the heme pocket of apoMb to form a complex in KPi buffer solution, the fluorescence transients exhibited a biphasic decay feature with the signal approaching an asymptotic offset: at lambda(em) = 600 nm, the rapid component decayed in 710 fs and the slow one in 27 ps; at lambda(em) = 680 nm, the two time constants were 950 fs and 40 ps. We conclude that (1) the fast-decay component pertains to an efficient transfer of energy from the hot S1 ZnPP species to apoMb through a dative bond between zinc and proximal histidine of the protein; (2) the slow-decay component arises from the water-induced vibrational relaxation of the hot S1 ZnPP species; and (3) the offset component is due to S1 --> T1 intersystem crossing of the surviving cold S1 ZnPP species. The transfer of energy through bonds might lead the dative bond to break, which explains our observation of the degradation of ZnPP-Mb samples in UV-vis and CD spectra upon protracted excitation.     

Relaxations Studies of 2IZSZ Spin Order with Signal Enhancement by Coherence Transfer

The pulse sequence for generating coherence transfer (or polarization transfer) is invoked to enhance the signal of heteronuclear two spin order (e.g. 2I_ZS_Z) in a spin system with CH moiety. This allows the observation of selective conversion of Zeeman order, in a sample with natural abundant ¹³C nuclei, into two spin order for measuring cross-correlation of chemical shift anisotropy and dipole-dipole interactions. The molecular reorientational correlation time and the orientation of the C–H axis with respect to the principal axes of carboxyl CSA tensor may be determined simultaneously in the relaxation profile of two spin order.     

A study of homonuclear dipolar recoupling pulse sequences in solid-state nuclear magnetic resonance

Dipolar recoupling pulse sequences are of great importance in magic angle spinning solid-state NMR. Recoupling sequences are used for excitation of double-quantum coherence, which, in turn, is employed in experiments to estimate internuclear distances and molecular torsion angles. Much effort is spent on the design of recoupling sequences that are able to produce double-quantum coherence with high efficiency in demanding spin systems, i.e., spin systems with small dipole-dipole couplings and large chemical-shift anisotropies (CSAs). The sequence should perform robustly under a variety of experimental conditions. This paper presents experiments and computer calculations that extend the theory of double-quantum coherence preparation from the strong coupling/small CSA limit to the weak coupling limit. The performance of several popular dipole-dipole recoupling sequences-DRAWS, POST-C7, SPC-5, R1, and R2-are compared. It is found that the optimum performance for several of these sequences, in the weak coupling/large CSA limit, varies dramatically, with respect to the sample spinning speed, the magnitude and orientation of the CSAs, and the magnitude of dipole-dipole couplings. It is found that the efficiency of double-quantum coherence preparation by gamma-encoded sequences departs from the predictions of first-order theory. The discussion is supported by density-matrix calculations.     

The anharmonic vibrational potential and relaxation pathways of the amide I and II modes of N-methylacetamide

We investigate the influence of isotopic substitution and solvation of N-methylacetamide (NMA) on anharmonic vibrational coupling and vibrational relaxation of the amide I and amide II modes. Differences in the anharmonic potential of isotopic derivatives of NMA in D2O and DMSO-d6 are quantified by extraction of the anharmonic parameters and the transition dipole moment angles from cross-peaks in the two-dimensional infrared (2D-IR) spectra. To interpret the effects of isotopic substitution and solvent interaction on the anharmonic potential, density functional theory and potential energy distribution calculations are performed. It is shown that the origin of anharmonic variation arises from differing local mode contributions to the normal modes of the NMA isotopologues, particularly in amide II. The time domain manifestation of the coupling is the coherent exchange of excitation between amide modes seen as the quantum beats in femtosecond pump-probes. The biphasic behavior of population relaxation of the pump-probe and 2D-IR experiments can be understood by the rapid exchange of strongly coupled modes within the peptide backbone, followed by picosecond dissipation into weakly coupled modes of the bath.     

High field 207Pb spin-lattice relaxation in solid lead nitrate and lead molybdate

Spin-lattice relaxation rates of lead have been measured at 17.6 T (156.9 MHz) as a function of temperature in polycrystalline lead nitrate and lead molybdate. Comparing the results with relaxation rates measured at lower fields, it is found that at high fields and low temperature, chemical shift anisotropy (CSA) makes small but observable contributions to lead relaxation in both materials. At 17.6 T and 200 K, CSA accounts for about 15% of the observed relaxation rate. Above 300 K, the dominant relaxation mechanism even at 17.6 T is an indirect Raman process involving modulation of the (207)Pb spin-rotation tensor, as first proposed by Grutzner et al. [J. Am. Chem. Soc. 123, 7094 (2001)] and later treated theoretically in more detail by Vega et al. [Phys. Rev. B 74, 214420 (2006)]. The improved signal to noise ratio at high fields makes it possible to quantify relaxation time anisotropy by analyzing saturation-recovery functions for individual frequencies on the powder pattern line shape. No orientation dependence is found for the spin-lattice relaxation rate of either material. It is argued from examination of the appropriate theoretical expressions, derived here for the first time, that the lack of observable relaxation time anisotropy is probably a general feature of this indirect Raman mechanism.     

Nuclear spin relaxation due to chemical shift anisotropy of gas-phase 129Xe

Nuclear spin relaxation provides detailed dynamical information on molecular systems and materials. Here, first-principles modeling of the chemical shift anisotropy (CSA) relaxation time for the prototypic monoatomic (129)Xe gas is carried out, both complementing and predicting the results of NMR measurements. Our approach is based on molecular dynamics simulations combined with pre-parametrized ab initio binary nuclear shielding tensors, an "NMR force field". By using the Redfield relaxation formalism, the simulated CSA time correlation functions lead to spectral density functions that, for the first time, quantitatively determine the experimental spin-lattice relaxation times T(1). The quality requirements on both the Xe-Xe interaction potential and binary shielding tensor are investigated in the context of CSA T(1). Persistent dimers Xe(2) are found to be responsible for the CSA relaxation mechanism in the low-density limit of the gas, completely in line with the earlier experimental findings.     

Site-specific 13C chemical shift anisotropy measurements in a uniformly 15N,13C-labeled microcrystalline protein by 3D magic-angle spinning NMR spectroscopy

In this Communication, we introduce a 3D magic-angle spinning recoupling experiment that correlates chemical shift anisotropy (CSA) powder line shapes with two dimensions of site-resolved isotropic chemical shifts. The principal tensor elements from 127 ROCSA line shapes are reported, constraining 102 unique backbone and side-chain 13C sites in a microcrystalline protein (the 56 residue beta1 immunoglobulin binding domain of protein G). The tensor elements, determined by fitting to numerical simulations, agree well with quantum chemical predictions. The experiments, therefore, validate calculations of CSAs in a protein of known structure. The data will be useful for the development of side-chain CSA quantum calculations and will aid in the design and interpretation of solution NMR experiments that utilize CSA-dipole cross-correlation to constrain torsion angles or to enhance resolution and sensitivity (such as in TROSY). Furthermore, the methodology described here will enable databases of CSA data to be generated with higher efficiency, for purposes of direct protein structure refinement.     

Cross-correlated relaxation between H1' chemical shift anisotropy and H1'-H2' dipolar relaxation mechanisms in ribonucleosides: application to the characterization of their anomeric configuration

Cross-correlated nuclear spin relaxation between 1H chemical shift anisotropy (CSA) and 1H-1H dipolar relaxation mechanisms in ribonucleosides in solution phase are observed and used to identify their anomeric configuration. Only alpha-ribonucleosides showed the presence of cross-correlated spin relaxation through differential spin-lattice relaxation (T1) of the H1' doublet. Dependence of the magnitude and the orientation of the H1' CSA tensor values on the glycosidic torsion angle and the fast time-scale internal motions present in the ribose moiety play a significant role in the characterization of the anomeric configuration of the nucleosides via cross-correlated relaxation.     

Determination of calpha chemical shift tensor orientation in peptides by dipolar-modulated chemical shift recoupling NMR spectroscopy

We present a new method for determining the orientation of chemical shift tensors in polycrystalline solids with site resolution and demonstrate its application to the determination of the Calpha chemical shift tensor orientation in a model peptide with beta-sheet torsion angles. The tensor orientation is obtained under magic angle spinning by modulating a recoupled chemical shift anisotropy (CSA) pattern with various dipolar couplings. These dipolar-modulated chemical shift patterns constitute the indirect dimension of a 2D spectrum and are resolved according to the isotropic chemical shifts of different sites in the direct dimension. These dipolar-modulated CSA spectra are equivalent to the projection of a 2D static separated-local-field spectrum onto its chemical shift dimension, except that its dipolar dimension is multiplied with a modulation function. Both (13)C-(1)H and (13)C-(15)N dipolar couplings can modulate the CSA spectra of the Calpha site in an amino acid and yield the relative orientations of the chemical shift principal axes to the C-H and C-N bonds. We demonstrate the C-H and C-N modulated CSA experiments on methylmalonic acid and N-tBoc-glycine, respectively. The MAS results agree well with the results of the 2D separated-local-field spectra, thus confirming the validity of this MAS dipolar-modulation approach. Using this technique, we measured the Val Calpha tensor orientation in N-acetylvaline, which has beta-sheet torsion angles. The sigma(11) axis is oriented at 158 degrees (or 22 degrees) from the C-H bond, while the sigma(22) axis is tilted by 144 degrees (or 36 degrees) from the C-N bond. Both the orientations and the magnitude of this chemical shift tensor are in excellent agreement with quantum chemical calculations.     

Femtosecond evolution of the pyrrole molecule excited in the near part of its UV spectrum

The evolution of the isolated pyrrole molecule has been followed after excitation in the 265-217 nm range by using femtosecond time delayed ionization. The transients collected in the whole excitation range show the vanishing of the ionization signal in the femtosecond time scale, caused by the relaxation along a πσ(?) type state (3s a(1)←π 1a(2)), which is the lowest excited electronic state of the molecule. This surface is dissociative along the NH bond, yielding a 15 ± 3 fs lifetime that reflects the loss of the ionization cross-section induced by the ultrafast wavepacket motion. Although a weak πσ(?) absorption is detected, the state is mainly reached through internal conversion of the higher bright ππ(?) transitions, which occurs with a 19 ± 3 fs lifetime. In addition to its resonant excitation, the intense ππ(?) absorption extending in the 220-190 nm interval is also out-of-resonance populated at energies far to the red from its absorption onset. This coherent adiabatic excitation of the ππ(?) transition should follow the excitation pulse (coherent population return effect), but instead the system relaxes toward the lower πσ(?) surface through a conical intersection during the interaction time, leading to the population of πσ(?) state at wavelengths as long as 265 nm. According to the observed behavior, the time evolution of the system in the full excitation range studied is modeled by a coherent treatment that provides key insights on the photophysical properties of the molecule.