DFT analysis of NMR scalar interactions across the glycosidic bond in DNA

The relationship between the glycosidic torsion angle chi, the three-bond couplings (3)J(C2/4-H1') and (3)J(C6/8-H1'), and the one-bond coupling (1)J(C1'-H1') in deoxyribonucleosides and a number of uracil cyclo-nucleosides has been analyzed using density functional theory. The influence of the sugar pucker and the hydroxymethyl conformation has also been considered. The parameters of the Karplus relationships between the three-bond couplings and chi depend strongly on the aromatic base. (3)J(C2/4-H1') reveals different behavior for deoxyadenosine, deoxyguanosine, and deoxycytidine as compared to deoxythymidine and deoxyuridine. In the case of (3)J(C6/8-H1'), an opposite trans to cis ratio of couplings is obtained for pyrimidine nucleosides in contrast to purine nucleosides. The extremes of the Karplus curves are shifted by ca. 10 degrees with respect to syn and anti-periplanar orientations of the coupled nuclei. The change in the sugar pucker from S to N decreases (3)J(C2/4-H1') and (3)J(C6/8-H1'), while increasing (1)J(C1'-H1') for the syn rotamers, whereas all of the trends are reversed for the anti rotamers. The influence of the sugar pucker on (1)J(C1'-H1') is interpreted in terms of interactions between the n(O4'), sigma*(C1'-H1') orbitals. The (1)J(C1'-H1') are related to chi through a generalized Karplus relationship, which combines cos(chi) and cos(2)(chi) functions with mutually different phase shifts that implicitly accounts for a significant portion of the related sugar pucker effects. Most of theoretical (3)J(C2/4-H1') and (3)J(C6/8-H1') for uracil cyclo-nucleosides compare well with available experimental data. (3)J(C6/8-H1') couplings for all C2-bridged nucleosides are up to 3 Hz smaller than in the genuine nucleosides with the corresponding chi, revealing a nonlocal aspect of the spin-spin interactions across the glycosidic bond. Theoretical (1)J(C1'-H1') are underestimated with respect to the experiment by ca. 10% but reproduce the trends in (1)J(C1'-H1') vs chi.     

Determination of DNA structure in solution: enzymatic deuteration of the ribose 2' carbon

An enzymatic solution to the problem of obtaining 13C/15N-labeled nucleotides that are deuterated uniquely at the H2' ' position within the ribose ring is presented. Selective deuteration occurs with an overall yield of >80%. The deuteron at the H2' ' position allows measurement of the scalar and residual dipolar couplings for the bond vectors attached to the C2' carbon of each ribose sugar. These data allow the accurate determination of sugar conformation. Interesting DNA double helices of 2-3 turns are now within the reach of solution NMR spectroscopy. As an example, these labeled nucleotides are incorporated uniquely at positions 6-14 in a 20-bp DNA sequence containing the adenovirus major late promoter.     

Measurement of ribose carbon chemical shift tensors for A-form RNA by liquid crystal NMR spectroscopy

Incomplete motional averaging of chemical shift anisotropy upon weak alignment of nucleic acids and proteins in a magnetic field results in small changes in chemical shift. Knowledge of nucleus-specific chemical shift (CS) tensor magnitudes and orientations is necessary to take full advantage of these measurements in biomolecular structure determination. We report the determination by liquid crystal NMR of the CS tensors for all ribose carbons in A-form helical RNA, using a series of novel 3D NMR pulse sequences for accurate and resolved measurement of the ribose (13)C chemical shifts. The orientation of the riboses relative to the rhombic alignment tensor of the molecule studied, a stem-loop sequence corresponding to helix-35 of 23S rRNA, is known from an extensive set of residual dipolar couplings (RDC), previously used to refine its structure. Singular-value-decomposition fits of the chemical shift changes to this structure, or alternatively to a database of helical RNA X-ray structures, provide the CS tensor for each type of carbon. Quantum chemical calculations complement the experimental results and confirm that the most shielded tensor component lies approximately along the local carbon-oxygen bond axis in all cases and that shielding anisotropy for C3' and C4' is much larger than for C1' and C2', with C5' being intermediate.     

Simultaneous NMR study of protein structure and dynamics using conservative mutagenesis

A novel iterative procedure is described that allows both the orientation and dynamics of internuclear bond vectors to be determined from direct interpretation of NMR dipolar couplings, measured under at least three orthogonal alignment conditions. If five orthogonal alignments are available, the approach also yields information on the degree of motional anisotropy and the direction in which the largest amplitude internal motion of each bond vector takes place. The method is demonstrated for the backbone (15)N-(1)H, (13)C(alpha)-(1)H(alpha), and (13)C(alpha)-13C' interactions in the previously well-studied protein domain GB3, dissolved in a liquid crystalline suspension of filamentous phage Pf1. Alignment variation is achieved by using conservative mutations of charged surface residues. Results indicate remarkably uniform backbone dynamics, with amplitudes that agree well with those of previous (15)N relaxation studies for most residues involved in elements of secondary structure, but larger amplitude dynamics than those found by (15)N relaxation for residues in loop and turn regions. In agreement with a previous analysis of dipolar couplings, the N-H bonds in the second beta-strand, which is involved in antibody recognition, show elevated dynamics with largest amplitudes orthogonal to the chain direction.     

1H-1H dipolar couplings provide a unique probe of RNA backbone structure

A NMR method is described that permits simultaneous measurement of the geminal 2JH1H2 + 2DH1H2 splitting and the sum of the 1JCH1 + 1DCH1 + 1JCH2 + 1DCH2 couplings for methylene groups, where 2DH1H2 and 1DCH are residual dipolar couplings, occurring when molecules are weakly oriented relative to the magnetic field. By suppressing either the upfield or downfield half of the 1H-1H geminal doublet, the experiment yields improved resolution relative to regular two-dimensional 1H-13C correlation spectra, making it applicable to systems of considerable complexity. The method is demonstrated for measurement of all 2DH5'H5' couplings in a 24-nucleotide 13C-enriched RNA stem loop structure, weakly aligned in liquid crystalline Pf1. The method is equally applicable to methylene groups in 13C-labeled proteins and to natural abundance samples of smaller molecules.     

Transition-state analysis of S. pneumoniae 5'-methylthioadenosine nucleosidase

Kinetic isotope effects (KIEs) and computer modeling are used to approximate the transition state of S. pneumoniae 5'-methylthioadenosine/S-adenosylhomocysteine nucleosidase (MTAN). Experimental KIEs were measured and corrected for a small forward commitment factor. Intrinsic KIEs were obtained for [1'-3H], [1'-14C], [2'-3H], [4'-3H], [5'-3H(2)], [9-15N] and [Me-3H(3)] MTAs. The intrinsic KIEs suggest an SN1 transition state with no covalent participation of the adenine or the water nucleophile. The transition state was modeled as a stable ribooxacarbenium ion intermediate and was constrained to fit the intrinsic KIEs. The isotope effects predicted a 3-endo conformation for the ribosyl oxacarbenium-ion corresponding to H1'-C1'-C2'-H2' dihedral angle of 70 degrees. Ab initio Hartree-Fock and DFT calculations were performed to study the effect of polarization of ribosyl hydroxyls, torsional angles, and the effect of base orientation on isotope effects. Calculations suggest that the 4'-3H KIE arises from hyperconjugation between the lonepair (n(p)) of O4' and the sigma* (C4'-H4') antibonding orbital owing to polarization of the 3'-hydroxyl by Glu174. A [methyl-3H(3)] KIE is due to hyperconjugation between np of sulfur and sigma* of methyl C-H bonds. The van der Waal contacts increase the 1'-3H KIE because of induced dipole-dipole interactions. The 1'-3H KIE is also influenced by the torsion angles of adjacent atoms and by polarization of the 2'-hydroxyl. Changing the virtual solvent (dielectric constant) does not influence the isotope effects. Unlike most N-ribosyltransferases, N7 of the leaving group adenine is not protonated at the transition state of S. pneumoniae MTAN. This feature differentiates the S. pneumoniae and E. coli transition states and explains the 10(3)-fold decrease in the catalytic efficiency of S. pneumoniae MTAN relative to that from E. coli.     

Sugar pucker modulates the cross-correlated relaxation rates across the glycosidic bond in DNA

The dependence of N1/9 and C1' chemical shielding (CS) tensors on the glycosidic bond orientation (chi) and sugar pucker (P) in the DNA nucleosides 2'-deoxyadenosine, 2'-deoxyguanosine, 2'-deoxycytidine, and 2'-deoxythymidine was studied using the calculation methods of quantum chemistry. The results indicate that these CS-tensors exhibit a significant degree of conformational dependence on chi and P structural parameters. The presented data test underlying assumptions of currently established methods for interpretation of cross-correlated relaxation rates between the N1/9 chemical shielding tensor and C1'-H1' dipole-dipole (Ravindranathan et al. J. Biomol. NMR 2003, 27, 365-75. Duchardt et al. J. Am. Chem. Soc. 2004, 126, 1962-70) and highlight possible limitations of these methods when applied to DNA.     

A ribose sugar conformational switch in the LTR-retrotransposon Ty3 polypurine tract-containing RNA/DNA hybrid

To probe structural features of a polypurine tract (PPT) that mediate its specific recognition and processing, a model 20 bp RNA/DNA hybrid duplex, which includes the full PPT sequence of the Saccharomyces cerevisiae LTR-retrotransposon Ty3, has been investigated using solution NMR spectroscopy. While homonuclear NOESY and DQF-COSY analyses indicate that this PPT-containing RNA/DNA hybrid adopts an overall A-form-like helical geometry, an unexpected sugar pucker switch has been detected for the ribose at position +1, relative to the cleavage site, on the RNA strand. A model of the conformational changes induced by the A- to B-type sugar pucker switch shows a significant change in the backbone trajectory of the RNA strand, which alters the presentation of backbone phosphate and 2' hydroxyl groups 3' of this residue. This observation implies that part of the mechanism governing RNase H fidelity may be through distortion of the RNA/DNA helix one base ahead of the scissile bond.     

Determination of natural abundance 15N-1H and 13C-1H dipolar couplings of molecules in a strongly orienting media using two-dimensional inverse experiments

NMR spectra of molecules oriented in liquid crystals provide homo- and heteronuclear dipolar couplings and thereby the geometry of the molecules. Several inequivalent dilute spins such as 13C and 15N coupled to protons form different coupled spin systems in their natural abundance and appear as satellites in the proton spectra. Identification of transitions belonging to each spin system is essential to determine heteronuclear dipolar couplings, which is a formidable task. In the present study, using 15N-1H and 13C-1H HSQC, and HMQC experiments we have selectively detected spectra of each rare spin coupled to protons. The 15N-1H and 13C-1H dipolar couplings have been determined in the natural abundance of 13C and 15N for the molecules pyrazine, pyrimidine and pyridazine oriented in a thermotropic liquid crystal.     

Theoretical study of the scalar coupling constants across the noncovalent contacts in RNA base pairs: the cis- and trans-watson-crick/sugar edge base pair family

The structure and function of RNA molecules are substantially affected by non-Watson-Crick base pairs actively utilizing the 2'-hydroxyl group of ribose. Here we correlate scalar coupling constants across the noncovalent contacts calculated for the cis- and trans-WC/SE (Watson-Crick/sugar edge) RNA base pairs with the geometry of base to base and sugar to base hydrogen bond(s). 23 RNA base pairs from the 32 investigated were found in RNA crystal structures, and the calculated scalar couplings are therefore experimentally relevant with regard to the binding patterns occurring in this class of RNA base pairs. The intermolecular scalar couplings 1hJ(N,H), 2hJ(N,N), 2hJ(C,H), and 3hJ(C,N) were calculated for the N-H...N and N-H...O=C base to base contacts and various noncovalent links between the sugar hydroxyl and RNA base. Also, the intramolecular 1J(N,H) and 2J(C,H) couplings were calculated for the amino or imino group of RNA base and the ribose 2'-hydroxyl group involved in the noncovalent interactions. The calculated scalar couplings have implications for validation of local geometry, show specificity for the amino and imino groups of RNA base involved in the linkage, and can be used for discrimination between the cis- and trans-WC/SE base pairs. The RNA base pairs within an isosteric subclass of the WC/SE binding patterns can be further sorted according to the scalar couplings calculated across different local noncovalent contacts. The effect of explicit water inserted in the RNA base pairs on the magnitude of the scalar couplings was calculated, and the data for discrimination between the water-inserted and direct RNA base pairs are presented. The calculated NMR data are significant for structural interpretation of the scalar couplings in the noncanonical RNA base pairs.     

NMR Structure Determination of Saccharose and Raffinose by Means of Homo‐ and Heteronuclear Dipolar Couplings

Residual dipolar couplings have dramatically improved the accuracy and precision of high‐resolution NMR structures during the last years. This was first demonstrated for proteins. In this article, we describe, with raffinose and saccharose as examples, that dipolar couplings improve the precision of structures of carbohydrates for which usually very few structural parameters are available. The relative orientation as well as the dynamics of the monosaccharide moieties with respect to each other can be determined with the help of ¹³C,¹H and ¹H,¹H dipolar couplings, which can easily be measured. Significant differences between the solution and the X‐ray crystal structure exist. These results indicate that residual dipolar‐coupling data may provide a more complete and dynamic model of carbohydrates in particular, and small molecules in general.     

Insights into the mobility of methyl-bearing side chains in proteins from (3)J(CC) and (3)J(CN) couplings

Side-chain dynamics in proteins can be characterized by the NMR measurement of (13)C and (2)H relaxation rates. Evaluation of the corresponding spectral densities limits the slowest motions that can be studied quantitatively to the time scale on which the overall molecular tumbling takes place. A different measure for the degree of side-chain order about the C(alpha)-C(beta) bond (chi(1) angle) can be derived from (3)J(C)(')(-)(C)(gamma) and (3)J(N)(-)(C)(gamma) couplings. These couplings can be measured at high accuracy, in particular for Thr, Ile, and Val residues. In conjunction with the known backbone structures of ubiquitin and the third IgG-binding domain of protein G, and an extensive set of (13)C-(1)H side-chain dipolar coupling measurements in oriented media, these (3)J couplings were used to parametrize empirical Karplus relationships for (3)J(C)(')(-)(C)(gamma) and (3)J(N)(-)(C)(gamma). These Karplus curves agree well with results from DFT calculations, including an unusual phase shift, which causes the maximum (3)J(CC) and (3)J(CN) couplings to occur for dihedral angles slightly smaller than 180 degrees, particularly noticeable in Thr residues. The new Karplus curves permit determination of rotamer populations for the chi(1) torsion angles. Similar rotamer populations can be derived from side-chain dipolar couplings. Conversion of these rotamer populations into generalized order parameters, S(J)(2) and S(D)(2), provides a view of side-chain dynamics that is complementary to that obtained from (13)C and (2)H relaxation. On average, results agree well with literature values for (2)H-relaxation-derived S(rel)(2) values in ubiquitin and HIV protease, but also identify a fraction of residues for which S(J,D)(2) < S(rel)(2). This indicates that some of the rotameric averaging occurs on a time scale too slow to be observable in traditional relaxation measurements.     

A new approach in 1D and 2D 13C high-resolution solid-state NMR spectroscopy of paramagnetic organometallic complexes by very fast magic-angle spinning

Novel 1D and multidimensional solid-state NMR (SSNMR) methods using very fast magic-angle spinning (VFMAS) (spinning speed > 20 kHz) for performing 13C high-resolution SSNMR of paramagnetic organometallic complexes are discussed. VFMAS removes a majority of 13C-1H and 1H-1H dipolar couplings, which are often difficult to remove by RF pulse techniques in paramagnetic complexes because of large paramagnetic shifts. In the first systematic approach using the unique feature of VFMAS for paramagnetic complexes, we demonstrate a means of obtaining well-resolved 1D and multidimensional 13C SSNMR spectra, sensitivity enhancements via cross polarization, and signal assignments, and applications of dipolar recoupling methods for nonlabeled paramagnetic organometallic complexes of moderate paramagnetic shifts ( approximately 800 ppm). Experimental results for powder samples of small nonlabeled coordination complexes at 1H frequencies of 400.2-400.3 MHz show that highly resolved 13C SSNMR spectra can be obtained under VFMAS, without requirements of 1H decoupling. Sensitivity enhancement in 13C SSNMR via cross polarization from 1H spins was demonstrated with an amplitude-sweep high-power CP sequence using strong RF fields ( approximately 100 kHz) available in the VFMAS probe. 13C CPMAS spectra of nonlabeled Cu(II)(dl-alanine)2.(H2O) and V(III)(acetylacetonate)3 (V(acac)3) show that it is possible to obtain high-resolution spectra for a small quantity ( approximately 15 mg) of nonlabeled paramagnetic organometal complexes within a few minutes under VFMAS. Experiments on Cu(II)(dl-alanine)2.(H2O) demonstrated that 1H-13C dipolar recoupling for paramagnetic organometal complexes can be performed under VFMAS by application of rotor-synchronous pi-pulses to 1H and 13C spins. The results also showed that signal assignments for 13CH, 13CH3, and 13CO groups in paramagnetic complexes are possible on the basis of the amount of 13C-1H dipolar dephasing induced by dipolar recoupling. Furthermore, the experimental 2D 13C/1H chemical-shift correlation NMR spectrum obtained for nonlabeled V(acac)3 exhibits well-resolved lines, which overlap in 1D 13C and 1H spectra. Signals for different chemical groups in the 2D spectrum are distinguished by the 13C-1H dipolar dephasing method combined with the 2D 13C/1H correlation NMR. The assignments offer information on the existence of nonequivalent ligands in the coordination complex in solids, without requiring a single-crystal sample.     

An alternative method for pucker determination in carbohydrates from residual dipolar couplings: a solution NMR study of the fructofuranosyl ring of sucrose

A simple solution NMR method is presented for pucker determination of five-membered rings using only residual dipolar couplings obtained in a single liquid crystalline medium, DMPC/DHPC bicelles (DMPC = dimyristoylphosphatidylcholine; DHPC = dihexanoylphosphatidylcholine). The method was applied to determine the pucker of the fructofuranosyl ring of sucrose. The results indicate a fructofuranosyl pucker phase in the 20 degrees - 70 degrees range. The pucker phases are in agreement with those from previous NMR and optical spectroscopic studies and, importantly, do not rely on empirically parametrized Karplus curves. Furthermore, the analysis implies more than one stable pucker phase and rapid ring interconversion in this range. The present results suggest that using residual dipolar couplings alone can reveal multiple conformations present in solution. Hence, when a sufficient number of residual dipolar coupling constants is measured, the outcome is a robust, reliable, and independent route for determining carbohydrate and nucleic acid structure by NMR spectroscopy.     

Selective measurement of heteronuclear 1H-13C dipolar couplings in motionally heterogeneous semicrystalline polymer systems

A pulse sequence for the selective recoupling of heteronuclear dipolar interactions in mobile amorphous phase of powdered semicrystalline polymers is described. 1H-13C dipolar interactions are selectively measured by PISEMA-type sequence. Selection of 13C magnetization originating from amorphous phase is achieved by a train of saturation pulses followed by a short delay and a direct excitation pulse on 13C spins. The development of undesired net 13C magnetization during the recoupling sequence is prevented by the efficient "reverse" 13C --> 1H cross-polarization. The efficacy of the 2D method to measure 1H-13C dipolar couplings selectively for mobile components is demonstrated on powdered crystalline L-alanine, semicrystalline polyethylene, and nanocomposite polyamide-6/montmorillonite.     

Comparison of methyl rotation axis order parameters derived from model-free analyses of (2)H and (13)C longitudinal and transverse relaxation rates measured in the same protein sample.

Recombinant HIV-1 protease was obtained from bacteria grown on a 98% D(2)O medium containing 3-(13)C pyruvic acid as the sole source of (13)C and (1)H. The purified protein is highly deuterated at non-methyl carbons, but contains significant populations of (13)CHD(2) and (13)CH(2)D methyl isotopomers. This pattern of isotope labeling permitted measurements of (1)H and (13)C relaxation rates of (13)CHD(2) isotopomers and (2)H (D) relaxation rates of (13)CH(2)D isotopomers using a single sample. The order parameters S(axis)(2), which characterize the motions of the methyl rotation axes, were derived from model-free analyses of R(1) and R(2) data sets measured for (13)C and (2)H spins. Our primary goal was to compare the S(axis)(2) values derived from the two independent types of data sets to test our understanding of the relaxation mechanisms involved. However, S(axis)(2) values derived from the analyses depend strongly on the geometry of the methyl group, the sizes of the quadrupolar and dipolar couplings, and the effects of bond vibrations and librations on these couplings. Therefore uncertainties in these basic physical parameters complicate comparison of the order parameters. This problem was circumvented by using an experimental relationship, between the methyl quadrupolar, (13)C-(13)C and (13)C-(1)H dipolar couplings, derived from independent measurements of residual static couplings of weakly aligned proteins by Ottiger and Bax (J. Am. Chem. Soc. 1999, 121, 4690-4695) and Mittermaier and Kay (J. Am. Chem. Soc. 1999, 121, 10608-10613). This approach placed a tight experimental restraint on the values of the (2)H quadrupolar and (13)C-(1)H dipolar interactions and greatly facilitated the accurate comparison of the relative values of the order parameters. When applied to our data this approach yielded satisfactory agreement between the S(axis)(2) values derived from the (13)C and (2)H data sets.     

Measurements of side-chain 13C-13C residual dipolar couplings in uniformly deuterated proteins

13C-only spectroscopy was used to measure multiple residual (13)C-(13)C dipolar couplings (RDCs) in uniformly deuterated and (13)C-labeled proteins. We demonstrate that (13)C-start and (13)C-observe spectra can be routinely used to measure an extensive set of the side-chain residual (13)C-(13)C dipolar couplings upon partial alignment of human ubiquitin in the presence of bacteriophages Pf1. We establish that, among different broadband polarization transfer schemes, the FLOPSY family can be used to exchange magnetization between a J coupled network of spins while largely decoupling dipolar interactions between these spins. An excellent correlation between measured RDCs and the 3D structure of the protein was observed, indicating a potential use of the (13)C-(13)C RDCs in the structure determination of perdeuterated proteins.     

Measurement of 13C-1H dipolar couplings in solids by using ultra-fast magic-angle spinning NMR spectroscopy with symmetry-based sequences

We show that (13)C-(1)H dipolar couplings in fully protonated organic solids can be measured by applying a Symmetry-based Resonance-Echo DOuble-Resonance (S-REDOR) experiment at ultra-fast Magic-Angle Spinning (MAS). The (13)C-(1)H dipolar couplings are recovered by using the R12 recoupling scheme, while the interference of (1)H-(1)H dipolar couplings are suppressed by the symmetry properties of this sequence and the use of high MAS frequency (65 kHz). The R12 method is especially advantageous for large (13)C-(1)H dipolar interactions, since the dipolar recoupling time can be incremented by steps as short as one rotor period. This allows a fine sampling for the rising part of the dipolar dephasing curve. We demonstrate experimentally that one-bond (13)C-(1)H dipolar coupling in the order of 22 kHz can be accurately determined. Furthermore, the proposed method allows a rapid evaluation of the dipolar coupling by fitting the S-REDOR dipolar dephasing curve with an analytical expression.