Quantum mechanical calculations of NMR chemical shifts in nucleic acids

Structure and Dynamics of ATP and the ATP-Zn2+ Complex in Solution

Despite the crucial role of ATP in life and artificial life-like applications, fundamental aspects relevant to its function, such as its conformational properties and its interaction with water and ions, remain unclear. Here, by integrating linear and two-dimensional infrared spectroscopy with ab initio molecular dynamics, we provide a detailed characterization of the vibrational spectra of the phosphate groups in ATP and in its complex with Zn2+ in water. Our study highlights the role of conformational disorder and solvation dynamics, beyond the harmonic normal-mode analysis, and reveals a complex scenario in which electronic and environmental effects tune the coupling between phosphate vibrations. We identify βγ-bidentate and αβγ-tridentate modes as the preferential coordination modes of Zn2+, as was proposed in the literature for Mg2+, although this conclusion is reached by a different spectral interpretation.

NMR chemical shift assignments of RNA oligonucleotides to expand the RNA chemical shift database

RNAs play myriad functional and regulatory roles in the cell. Despite their significance, three-dimensional structure elucidation of RNA molecules lags significantly behind that of proteins. NMR-based studies are often rate-limited by the assignment of chemical shifts. Automation of the chemical shift assignment process can greatly facilitate structural studies, however, accurate chemical shift predictions rely on a robust and complete chemical shift database for training. We searched the Biological Magnetic Resonance Data Bank (BMRB) to identify sequences that had no (or limited) chemical shift information. Here, we report the chemical shift assignments for 12 RNA hairpins designed specifically to help populate the BMRB.

Classification of RNA backbone conformations into rotamers using 13C' chemical shifts: exploring how far we can go

The conformational space of the ribose-phosphate backbone is very complex as it is defined in terms of six torsional angles. To help delimit the RNA backbone conformational preferences, 46 rotamers have been defined in terms of these torsional angles. In the present work, we use the ribose experimental and theoretical ¹³C′ chemical shifts data and machine learning methods to classify RNA backbone conformations into rotamers and families of rotamers. We show to what extent the experimental ¹³C′ chemical shifts can be used to identify rotamers and discuss some problem with the theoretical computations of ¹³C′ chemical shifts.

A simple and fast approach for predicting (1)H and (13)C chemical shifts: toward chemical shift-guided simulations of RNA

We introduce a simple and fast approach for predicting RNA chemical shifts from interatomic distances that performs with an accuracy similar to existing predictors and enables the first chemical shift-restrained simulations of RNA to be carried out. Our analysis demonstrates that the applied restraints can effectively guide conformational sampling toward regions of space that are more consistent with chemical shifts than the initial coordinates used for the simulations. As such, our approach should be widely applicable in mapping the conformational landscape of RNAs via chemical shift-guided molecular dynamics simulations. The simplicity and demonstrated sensitivity to three-dimensional structure should also allow our method to be used in chemical shift-based RNA structure prediction, validation, and refinement.

Spectroscopic study of proflavine adsorption on the carbon nanotube surface

Despite the fact that non-covalent interactions between various aromatic compounds and carbon nanotubes are being extensively investigated now, there is still a lack of understanding about the nature of such interactions. The present paper sheds light on one of the possible mechanisms of interaction between the typical aromatic dye proflavine and the carbon nanotube surface, namely, π-stacking between aromatic rings of these compounds. To investigate such a complexation, a qualitative analysis was performed by means of ultraviolet visible, infrared, and nuclear magnetic resonance spectroscopy. The data obtained suggest that π-stacking brings the major contribution to the stabilization of the complex between proflavine and the carbon nanotube.

Prediction of RNA 1H and 13C chemical shifts: a structure based approach

The use of NMR-derived chemical shifts in protein structure determination and prediction has received much attention, and, as such, many methods have been developed to predict protein chemical shifts from three-dimensional (3D) coordinates. In contrast, little attention has been paid to predicting chemical shifts from RNA coordinates. Using the random forest machine learning approach, we developed RAMSEY, which is capable of predicting both (1)H and protonated (13)C chemical shifts from RNA coordinates. In this report, we introduce RAMSEY, assess its accuracy, and demonstrate the sensitivity of RAMSEY-predicted chemical shifts to RNA 3D structure.

Nucleic acid helix structure determination from NMR proton chemical shifts

We present a method for de novo derivation of the three-dimensional helix structure of nucleic acids using non-exchangeable proton chemical shifts as sole source of experimental restraints. The method is called chemical shift de novo structure derivation protocol employing singular value decomposition (CHEOPS) and uses iterative singular value decomposition to optimize the structure in helix parameter space. The correct performance of CHEOPS and its range of application are established via an extensive set of structure derivations using either simulated or experimental chemical shifts as input. The simulated input data are used to assess in a defined manner the effect of errors or limitations in the input data on the derived structures. We find that the RNA helix parameters can be determined with high accuracy. We finally demonstrate via three deposited RNA structures that experimental proton chemical shifts suffice to derive RNA helix structures with high precision and accuracy. CHEOPS provides, subject to further development, new directions for high-resolution NMR structure determination of nucleic acids.

Chemical shifts in nucleic acids studied by density functional theory calculations and comparison with experiment

NMR chemical shifts are highly sensitive probes of local molecular conformation and environment and form an important source of structural information. In this study, the relationship between the NMR chemical shifts of nucleic acids and the glycosidic torsion angle, χ, has been investigated for the two commonly occurring sugar conformations. We have calculated by means of DFT the chemical shifts of all atoms in the eight DNA and RNA mono-nucleosides as a function of these two variables. From the DFT calculations, structures and potential energy surfaces were determined by using constrained geometry optimizations at the BP86/TZ2P level of theory. The NMR parameters were subsequently calculated by single-point calculations at the SAOP/TZ2P level of theory. Comparison of the (1)H and (13)C NMR shifts calculated for the mono-nucleosides with the shifts determined by NMR spectroscopy for nucleic acids demonstrates that the theoretical shifts are valuable for the characterization of nucleic acid conformation. For example, a clear distinction can be made between χ angles in the anti and syn domains. Furthermore, a quantitative determination of the χ angle in the syn domain is possible, in particular when (13)C and (1)H chemical shift data are combined. The approximate linear dependence of the C1' shift on the χ angle in the anti domain provides a good estimate of the angle in this region. It is also possible to derive the sugar conformation from the chemical shift information. The DFT calculations reported herein were performed on mono-nucleosides, but examples are also provided to estimate intramolecularly induced shifts as a result of hydrogen bonding, polarization effects, or ring-current effects.

Structural analysis of the binding of the diquaternary pyridophenazine derivative dqdppn to B-DNA oligonucleotides

The interaction of the ethylene-bipyridyldiylium-naphthaphenazine dication, dqdppn, with several hexa- and octanucleotide duplexes has been studied using CD and NMR. Taken together, these studies reveal that with the hexanucleotide, dqdppn intercalates into the terminal base pair, and causes a large twisting of the terminal base pair. In contrast, with all three octanucleotides, dqdppn intercalates more centrally within the sequence. The NMR-derived structures of two of the binding complexes demonstrate that dqdppn intercalates from the major groove in an unusual 'side-on' geometry, rather than threading through the helix. An analysis of these results indicates that the preferred binding site is not sequence-specific, but primarily at the most conformationally flexible DNA step.

Structure-specific binding of [Co(phen)(2)(HPIP)](3+) to a DNA duplex containing sheared G:A mismatch base pairs

The binding of a Co(III) complex to the decanucleotide d(CCGAATGAGG)(2) containing two pairs of G:A mismatches was studied by 2D-NMR, UV absorption, and molecular modeling. NMR investigations indicate that racemic [Co(phen)(2)(HPIP)]Cl(3) [HPIP=2-(2-hydroxyphenyl) imidazo [4,5-f][1,10] phenanthroline] binds the decanucleotide by intercalation: the HPIP ligand selectively inserts between the stacked bases from the minor groove at the terminal regions and from the major groove at the sheared region. Further, molecular modeling revealed that the recognition shows strong enantioselectivity: the Lambda-isomer preferentially intercalates into the T(6)G(7):A(5)A(4) region from the DNA major groove, while Delta-isomer favors the terminal C(1)C(2):G(10)G(9) region and intercalates from the minor groove. Detailed energy analysis suggests that the steric interaction, especially the electrostatic effect, is the primary determinants of the recognition event. Melting experiments indicate that binding stabilizes the DNA duplex and increases the melting temperature by 9.5 degrees C. The intrinsic binding constant of the complex to the mismatched duplex was determined to be 3.5x105M(-1).

Recent advances in RNA structure determination by NMR

Despite recent advances in the solution of NMR structures of RNA and RNA-ligand complexes, the rate limiting step remains the gathering of a large number of NOE and torsion restraints. Additional sources of information for structure determination of larger RNA molecules have recently become available, and it is possible to supplement NOE and J-coupling data with the measurement of dipolar couplings and cross-correlated relaxation rates in high-resolution NMR spectroscopy.

Effect of methylation on the pyrimidine-pyrimidine stacking interaction studied by (1)H NMR chemical shift

Mutually induced proton chemical shift changes were measured for the mixed solutions of pyrimidine and its methylated forms in deuterium oxide at 35 degrees C. The chemical shift vs. concentration profiles were analyzed using a three-state decomposition model based on competitive self- and hetero-association dimer equilibria. The equilibrium constants show an increasing association tendency within the series pyrimidine-5-methyl-pyrimidine (0.23 +/- 0.02 M(-1)) < pyrimidine-4,6-dimethyl-pyrimidine (0.32 +/- 0.04 M(-1)) < 5-methyl-pyrimidine-4,6-dimethyl-pyrimidine (0.51 +/- 0.04 M(-1)). The upfield dimer shifts suggest an offset stacked geometry for the structure of associations between the parent molecule of the pyrimidine nucleobases and its methylated derivatives in aqueous solution.

DSHIFT: a web server for predicting DNA chemical shifts

DSHIFT is a web server for predicting chemical shifts of DNA sequences in random coil form or double helical B-form. The prediction methods are based on sets of published reference chemical shift values and correction factors which account for shielding or deshielding effects from neighboring nucleotides. Proton, carbon and phosphorus chemical shift predictions are available for random coil DNAs. For double helical B-DNA, only proton chemical shift prediction is available. Results from these predictions will be useful for facilitating NMR resonance assignments and investigating structural features of solution DNA molecules. The URL of this server is: http://www.chem.cuhk.edu.hk/DSHIFT.

Identification of the alpha and beta anomers of 1-(2-deoxy-D-erythro-pentofuranosyl)-oxaluric acid at the site of riboflavin-mediated photooxidation of guanine in 2'-deoxyguanosine and thymidylyl-(3'-5')-2'-deoxyguanosine

Products of riboflavin-mediated photosensitization of 2'-deoxyguanosine (dG) and thymidylyl-(3'-5')-2'-deoxyguanosine (TpdG) by 350-nm light in oxygen-saturated aqueous solution have been isolated and identified as 1-(2-deoxy-beta-D-erythro-pentofuranosyl) oxaluric acid (beta-dOx) and thymidylyl-(3'-5')-1-(2-deoxy-beta-D-erythro-pentofuranosyl) oxaluric acid (Tpbeta-dOx), respectively. In aqueous solution the modified beta-deoxyribonucleoside is slowly converted to the alpha-anomer, generating alpha-dOx and Tpalpha-dOx. These modified nucleosides and dinucleoside monophosphates have been isolated by HPLC and characterized by proton and carbon NMR spectroscopy, fast atom bombardment mass spectrometry, and enzymatic analyses. Both alpha-dOx and Tpalpha-dOx slowly convert back into the modified beta-deoxyribonucleoside, indicating that the furanosidic anomers are in dynamic equilibrium. Relative to TpdG, the rate of hydrolysis of Tpbeta-dOx and Tpalpha-dOx by spleen phosphodiesterase is greatly reduced. Hot piperidine (1.0 M, 90 degrees C, 30 min) destroys Tpbeta-dOx and Tpalpha-dOx. Riboflavin-mediated photosensitization of TpdG in D2O instead of H2O has no detectable effect on the yield of Tpbeta-dOx, suggesting that oxaluric acid is generated through a Type-I reaction mechanism, likely through the intermediary on initially generated 8-oxo-7,8-dihydro-2'-deoxyguanosine.

Hetero-association of anticancer antibiotics in aqueous solution: NMR and molecular mechanics analysis

In order to investigate the effect on combinations of aromatic antibiotics used in chemotherapy, the hetero-association of the antitumour antibiotics actinomycin D (AMD) with daunomycin (DAU) or novatrone (NOV) has been studied by the methods of 1D- and 2D 500 MHz 1H-NMR spectroscopy and molecular mechanics calculations. The experimental concentration and temperature dependences of the proton chemical shifts of mixtures of the aromatic drugs have been analyzed in terms of a modified statistical-thermodynamical model of hetero-association to give the equilibrium reaction constants, the thermodynamical parameters (deltaH, deltaS) of hetero-association of AMD with DAU or NOV and the limiting values of proton chemical shifts of the molecules in the hetero-complexes. The most favorable averaged structures of the 1:1 DAU-AMD and NOV-AMD hetero-association complexes have been determined using both the limiting values of proton chemical shifts of the molecules and molecular mechanics methods (X-PLOR software). The results show that intermolecular complexes between DAU-AMD and NOV-AMD are mainly stabilized by stacking interactions of the aromatic chromophores, although the DAU-AMD hetero-complex has additional stabilization, which may be explained by an intermolecular hydrogen bond between a carbonyl group of ring C of DAU and the NH group of D-Val of the pentapeptide side chain ring of AMD. The relative content of each type of molecular complex in the mixed solution has been calculated at different values of the ratio (r) of the initial concentrations of DAU and AMD. It is found that the contributions of hetero-complexes to the general equilibrium in solution are predominant at quite different values of r, viz. at r>12 for AMD with NOV and at r>2 for AMD with DAU, compared to r>0.3 for the DAU-NOV system observed previously. It is concluded that anticancer drugs have quite different affinities for formation of hetero-complexes with other aromatic antibiotics in aqueous solution, which may need to be taken into consideration for their use in combination chemotherapy.

Random coil phosphorus chemical shift of deoxyribonucleic acids

Random coil phosphorus chemical shift has been studied using 16 17-nucleotide DNA sequences. Due to the presence of residual base stacking in these sequences, the temperature and sequence effects were investigated at 50 and 55 degrees C. The phosphorus chemical shifts of random coil DNA sequences have been found to be independent of temperature. Sequence effect analysis shows that the phosphorus chemical shift of a nucleotide in a random coil DNA sequence depends on both its 5'- and 3'-nearest neighbors. A trimer model has been used to establish the random coil 31P chemical shift prediction protocol which shows an accuracy of 0.02 ppm.

Calculation of 13C chemical shifts in rna nucleosides: structure-13C chemical shift relationships

Isotropic 13C chemical shifts of the ribose sugar in model RNA nucleosides are calculated using SCF and DFT-GIAO ab initio methods for different combinations of ribose sugar pucker, exocyclic torsion angle, and glycosidic torsion angle. Idealized conformations were obtained using structures that were fully optimized by ab initio DFT methods starting with averaged parameters from a collection of crystallographic data. Solid-state coordinates of accurate crystal or neutron diffraction structures were also examined directly without optimization. The resulting 13C chemical shifts for the two sets of calculations are then compared. The GIAO-DFT method overestimates the shifts by an average of 5 ppm while the GIAO-SCF underestimates the shifts by the same amount. However, in the majority of cases the errors appear to be systematic, as the slope of a plot of calculated vs experimental shifts is very close to unity, with minimal scatter. The values of the 13C NMR shifts of the ribose sugar are therefore sufficiently precise to allow for statistical separation of sugar puckering modes and exocyclic torsion angle conformers, based on the canonical equation model formulated in a previous paper.

2-Deoxy-beta-D-erythro-pentofuranose: hydroxymethyl group conformation and substituent effects on molecular structure, ring geometry, and NMR spin-spin coupling constants from quantum chemical calculations

The effect of hydroxymethyl conformation (gg, gt, and tg rotamers about the C4-C5 bond) on the conformational energies and structural parameters (bond lengths, bond angles, bond torsions) of the 10 envelope forms of the biologically relevant aldopentofuranose, 2-deoxy-beta-D-erythro-pentofuranose (2-deoxy-D-ribofuranose) 2, has been investigated by ab initio molecular orbital calculations at the HF/6-31G level of theory. C4-C5 bond rotation induces significant changes in the conformational energy profile of 2 (2gt and 2tg exhibit one global energy minimum, whereas 2gg exhibits two nearly equivalent energy minima), and structural changes, especially those in bond lengths, are consistent with predictions based on previously reported vicinal, 1,3- and 1,4-oxygen lone pair effects. HF/6-31G-optimized envelope geometries of 2gg were re-optimized using density functional theory (DFT, B3LYP/6-31G), and the resulting structures were used in DFT calculations of NMR spin-spin coupling constants involving 13C (i.e., J(CH) and J(CC) over one, two, and three bonds) in 2gg according to methods described previously. The computed J-couplings were compared to those reported previously in 2gt to assess the effect of C4-C5 bond rotation on scalar couplings within the furanose ring and hydroxymethyl side chain. The results confirm prior predictions of correlations between 2J(CH), 3J(CH), 2J(CC) and 3J(CC), and ring conformation, and verify the usefulness of a concerted application of these couplings (both their magnitudes and signs) in assigning preferred ring and C4-C5 bond conformations in aldopentofuranosyl rings. The new calculated J-couplings in 2gg have particular relevance to related J-couplings in DNA (and RNA indirectly), where the gg rotamer, rather than the gt rotamer, is observed in most native structures. The effects of two additional structural perturbations on 2 were also studied, namely, deoxygenation at C5 (yielding 2,5-dideoxy-beta-D-erythro-pentofuranose 4) and methyl glycosidation at O1 (yielding methyl 2-deoxy-beta-D-erythro-pentofuranoside 5) at the HF/6-31G level. The conformational energy profile of 4 resembles that found for 2gt, not 2gg, indicating that 4 is an inappropriate structural mimic of the furanose ring in DNA. Glycosidation failed to induce differential stabilization of ring conformations containing an axial C1-O1 bond (anomeric effect), contrary to experimental data. The latter discrepancy indicates that either the magnitude of this differential stabilization depends on ring configuration or that solvent effects, which are neglected in these calculations, play a role in promoting this stabilization.

Dependence of 13C NMR chemical shifts on conformations of rna nucleosides and nucleotides

Cross-polarization magic-angle spinning solid-state NMR spectroscopy has been used to investigate the dependence of (13)C sugar chemical shifts on specific conformational parameters of a variety of ribonucleotides and ribonucleosides. Solid-state NMR is a valuable tool for nucleoside and nucleotide structural studies since it provides the means to acquire spectra that correspond to single conformations, as opposed to (13)C solution NMR methods. The distinct effects of sugar puckering on the C1', C4', and C5' resonances of C2' endo (S type) and C3' endo (N type) furanoid conformations allow us to separate them into two groups. Further analysis of each group reveals an additional dependence of the C1' and C5' resonances on the glycosidic and C4'-C5' exocyclic torsion angles, respectively. However, it is found that the glycosidic conformation cannot independently be determined from sugar (13)C chemical shift data. The statistical methods of exploratory data analysis and discriminant analysis are used to construct two canonical coordinates-linear combinations of chemical shifts which give the statistically optimal determination of the conformation from the NMR data.

A-like guanine-guanine stacking in the aqueous DNA duplex of d(GGGGCCCC)

We have used CD spectroscopy, NMR spectroscopy and unrestrained molecular dynamics to study conformational properties of a DNA duplex formed by the self-complementary octamer d(GGGGCCCC). Its unusual CD spectrum contains features indicating A-like stacking of half of the bases, whereas the other half stack in a B-like fashion. Unrestrained molecular dynamics simulations converged to a stable B-like double-helix of d(GGGGCCCC). However, the double-helix contained a central hole whose size was half of that occurring in structure A. In the canonical structure B, the hole does not exist at all because the base-pairs cross the double-helix centre. The cytosine bases were stacked in the duplex of d(GGGGCCCC) as in structure B, while stacking of the guanine bases displayed features characteristic for structure A. NMR spectroscopy revealed that the A-like guanine-guanine stacking was accompanied by an increased tendency of the deoxyribose rings attached to the guanine bases to be puckered in an A-like fashion. Otherwise, the duplex of d(GGGGCCCC) showed no clash, no bend and no other significant deviation from structure B. The present analysis demonstrates a remarkable propensity of the guanine runs to stack in an A-like fashion even within the B-DNA framework. This property explains why the oligo(dG). oligo(dC) tracts switch into structure A so easily. Secondly, this property may influence replication, because structure A is replicated more faithfully than structure B. Thirdly, the oligo(dG) runs might have played an important role in early evolution, when DNA took on functions that originally evolved on RNA. Fourthly, the present study extends the vocabulary of DNA secondary structures by the heteronomous duplex of d(GGGGCCCC) in which the B-like strand of oligo(dC) is bound to the A-like strand of oligo(dG).

The C-A mismatch base pair and the single-strand terminus in the E. coli initiator tRNA(fMet) acceptor stem adopt unusual conformations

Acceptor stem variants of tRNA(fMet) (Escherichia coli) have been characterized by nuclear magnetic resonance. The wild type contains a C1-A72 mismatch pair which is crucial for its biological function. For comparison, the mismatch was replaced by regular pairs U1-A72 and C1-G72. Further variants contain an altered discriminator base, G73, or a G1-C72/U73 combination. The stems of variants U1-A72/A73 and C1-G72/A73 have A-RNA geometry, which extends essentially to the single-strand terminus. C1-A72/G73 variant and wild type are structurally almost identical. C1 and A72 adopt peculiar conformations with C1 being largely destacked with respect to G2, while A73 stacks upon C1. The unique arrangement of the mismatch causes a distinctly different orientation of the single-strand terminus compared to variants with regular 1-72 base pairs, and to formyltransferase-complexed tRNA(fMet).

An A-type double helix of DNA having B-type puckering of the deoxyribose rings

DNA usually adopts structure B in aqueous solution, while structure A is preferred in mixtures of trifluoroethanol (TFE) with water. However, the octamer d(CCCCGGGG) and other d(C(n)G(n)) fragments of DNA provide CD spectra that suggest that the base-pairs are stacked in an A-like fashion even in aqueous solution. Yet, d(CCCCGGGG) undergoes a cooperative TFE-induced transition into structure A, indicating that an important part of the aqueous duplex retains structure B. NMR spectroscopy shows that puckering of the deoxyribose rings is of the B-type. Hence, combination of the information provided by CD spectroscopy and NMR spectroscopy suggests an unprecedented double helix of DNA in which A-like base stacking is combined with B-type puckering of the deoxyribose rings. In order to determine whether this combination is possible, we used molecular dynamics to simulate the duplex of d(CCCCGGGG). Remarkably, the simulations, completely unrestrained by the experimental data, provided a very stable double helix of DNA, exhibiting just the intermediate B/A features described above. The double helix contained well-stacked guanine bases but almost unstacked cytosine bases. This generated a hole in the double helix center, which is a property characteristic for A-DNA, but absent from B-DNA. The minor groove was narrow at the double helix ends but wide at the central CG step where the Watson-Crick base-pairs were buckled in opposite directions. The base-pairs stacked tightly at the ends but stacking was loose in the duplex center. The present double helix, in which A-like base stacking is combined with B-type sugar puckering, is relevant to replication and transcription because both of these phenomena involve a local B-to-A transition.

NMR investigation of the complexation of daunomycin with deoxytetranucleotides of different base sequence in aqueous solution

500 MHz NMR spectroscopy has been used to investigate the complexation of the anthracycline antibiotic daunomycin (DAU) with self-complementary deoxytetranucleotides, 5'-d(CGCG), 5'-d(GCGC), 5'-d(TGCA), 5'-d(ACGT) and 5'-d(AGCT), of different base sequence in aqueous salt solution. 2D homonuclear 1H NMR spectroscopy (TOCSY and NOESY) and heteronuclear 1H - 31P NMR spectroscopy (HMBC) have been used for complete assignment of the non-exchangeable protons and the phosphorus resonance signals, respectively, and for a qualitative determination of the preferred binding sites of the drug. Analysis shows that DAU intercalates preferentially into the terminal sites of each of the tetranucleotides and that the aminosugar of the antibiotic is situated in the minor groove of the tetramer duplex, partly eclipsing the third base pair. A quantitative determination of the complexation of DAU with the deoxytetranucleotides has been made using the experimental concentration and temperature dependences of the drug proton chemical shifts; these have been analysed in terms of the equilibrium reaction constants, limiting proton chemical shifts and thermodynamical parameters (enthalpies deltaH, entropies deltaS) of different drug-DNA complexes (1:1, 1:2, 2:1, 2:2) in aqueous solution. It is found that DAU interacts with sites containing three adjacent base pairs but does not show any significant sequence specificity of binding with either single or double-stranded tetranucleotides, in contrast with other intercalating drugs such as proflavine, ethidium bromide and actinomycin D. The most favourable structures of the 1:2 complexes have been derived from the induced limiting proton chemical shifts of the drug in the intercalated complexes with the tetranucleotide duplex, in conjunction with 2D NOE data. It has been found that the conformational parameters of the double helix and the orientation of the DAU chromophore in the intercalated complexes depend on base sequence at the binding site of the tetramer duplexes in aqueous solution.

Detection of N-H...N hydrogen bonding in RNA via scalar couplings in the absence of observable imino proton resonances

Hydrogen bond networks stabilize RNA secondary and tertiary structure and are thus essentially important for protein recognition. During structure refinements using either NMR or X-ray techniques, hydrogen bonds were usually inferred indirectly from the proximity of donor and acceptor functional groups. Recently, quantitative heteronuclear J(N,N)-HNN COSY NMR experiments were introduced that allowed the direct identification of donor and acceptor nitrogen atoms involved in hydrogen bonds. However, protons involved in base pairing interactions in nucleic acids are often not observable due to exchange processes. The application of a modified quantitative J(N,N)-HNN COSY pulse scheme permits observation of(2h)J(N,N) couplings via non-exchangeable protons. This approach allowed the unambiguous identification of the A27.U23 reverse Hoogsteen base pair involved in a U-A.U base triple in the HIV-2 transactivation response element-argininamide complex. Despite a wealth of NOE information, direct evidence for this interaction was lacking due to the rapid exchange of the U23 imino proton. The ability to directly observe hydrogen bonds, even in D(2)O and in the presence of rapid exchange, should facilitate structural studies of RNA.