Hydrogen bonding of nucleotide base pairs: Application of the PM3 method

The ability of the pm3 semiempirical quantum mechanical method to reproduce hydrogen bonding in nucleotide base pairs was assessed. Results of pm3 calculations on the nucleotides 2′-deoxyadenosine 5′-monophosphate (pdA), 2′-deoxyguanosine 5′-monophosphate (pdG), 2′-deoxycytidine 5′-monophosphate (pdC), and 2′-deoxythymidine 5′-monophosphate (pdT) and the base pairs pdA–pdT, pdG–pdC, and pdG(syn)–pdC are presented and discussed. The pm3 method is the first of the parameterized nddo quantum mechanical models with any ability to reproduce hydrogen bonding between nucleotide base pairs. Intermolecular hydrogen bond lengths between nucleotides displaying Watson–Crick base pairing are 0.1–0.2 Å less than experimental results. Nucleotide bond distances, bond angles, and torsion angles about the glycosyl bond (χ), the C4′C5′ bond (γ), and the C5′O5′ bond (β) agree with experimental results. There are many possible conformations of nucleotides. pm3 calculations reveal that many of the most stable conformations are stabilized by intramolecular CHO hydrogen bonds. These interactions disrupt the usual sugar puckering. The stacking interactions of a dT–pdA duplex are examined at different levels of gradient optimization. The intramolecular hydrogen bonds found in the nucleotide base pairs disappear in the duplex, as a result of the additional constraints on the phosphate group when part of a DNA backbone. Sugar puckering is reproduced by the pm3 method for the four bases in the dT–pdA duplex. pm3 underestimates the attractive stacking interactions of base pairs in a B-DNA helical conformation. The performance of the pm3 method implemented in SPARTAN is contrasted with that implemented in MOPAC. At present, accurate ab initio calculations are too timeconsuming to be of practical use, and molecular mechanics methods cannot be used to determine quantum mechanical properties such as reaction-path calculations, transition-state structures, and activation energies. The pm3 method should be used with extreme caution for examination of small DNA systems. Future parameterizations of semiempirical methods should incorporate base stacking interactions into the parameterization data set to enhance the ability of these methods.

Caractérisation structurale et thermodynamique de la reconnaissance du substrat par le ribozyme VS de Neurospora

Les interactions ARN/ARN de type kissing-loop sont des éléments de structure tertiaire qui jouent souvent des rôles clés chez les ARN, tant au niveau fonctionnel que structural. En effet, ce type d’interaction est crucial pour plusieurs processus dépendant des ARN, notamment pour l’initiation de la traduction, la reconnaissance des ARN antisens et la dimérisation de génome rétroviral. Les interactions kissing-loop sont également importantes pour le repliement des ARN, puisqu’elles permettent d’établir des contacts à longue distance entre différents ARN ou encore entre les domaines éloignés d’un même ARN. Ce type d’interaction stabilise aussi les structures complexes des ARN fonctionnels tels que les ARNt, les riborégulateurs et les ribozymes. Comme d’autres ARN fonctionnels, le ribozyme VS de Neurospora contient une interaction kissing-loop importante. Celle-ci est impliquée dans la reconnaissance du substrat et se forme entre la tige-boucle I (stem-loop I, SLI) du substrat et la tige-boucle V (stem-loop V, SLV) du domaine catalytique. Des études biochimiques ont démontré que l’interaction kissing-loop I/V, dépendante du magnésium, implique trois paires de bases Watson-Crick (W-C). De plus, cette interaction est associée à un réarrangement de la structure du substrat, le faisant passer d’une conformation inactive dite unshifted à une conformation active dite shifted. Les travaux présentés dans cette thèse consistent en une caractérisation structurale et thermodynamique de l’interaction kissing-loop I/V du ribozyme VS, laquelle est formée de fragments d’ARN représentant les tige-boucles I et V dérivées du ribozyme VS (SLI et SLV). Cette caractérisation a été réalisée principalement par spectroscopie de résonance magnétique nucléaire (RMN) et par titrage calorimétrique isotherme (isothermal titration calorimetry, ITC) en utilisant différents complexes SLI/SLV dans lesquels l’ARN SLV est commun à tous les complexes, alors que différentes variations de l’ARN SLI ont été utilisées, soit en conformation shiftable ou preshifted. Les données d’ITC ont permis de démontrer qu’en présence d’une concentration saturante de magnésium, l’affinité d’un substrat SLI preshifted pour SLV est extrêmement élevée, rendant cette interaction plus stable que ce qui est prédit pour un duplexe d’ARN équivalent. De plus, l’étude effectuée par ITC montre que des ARN SLI preshifted présentent une meilleure affinité pour SLV que des ARN SLI shiftable, ce qui a permis de calculer le coût énergétique associé au réarrangement de structure du substrat. En plus de confirmer la formation des trois paires de bases W-C prédites à la jonction I/V, les études de RMN ont permis d’obtenir une preuve structurale directe du réarrangement structural des substrats SLI shiftable en présence de magnésium et de l’ARN SLV. La structure RMN d’un complexe SLI/SLV de grande affinité démontre que les boucles terminales de SLI et SLV forment chacune un motif U-turn, ce qui facilite l’appariement W-C intermoléculaire. Plusieurs autres interactions ont été définies à l’interface I/V, notamment des triplets de bases, ainsi que des empilements de bases. Ces interactions contribuent d’ailleurs à la création d’une structure présentant un empilement continu, c’est-à-dire qui se propage du centre de l’interaction jusqu’aux bouts des tiges de SLI et SLV. Ces études de RMN permettent donc de mieux comprendre la stabilité exceptionnelle de l’interaction kissing-loop I/V au niveau structural et mènent à l’élaboration d’un modèle cinétique de l’activation du substrat par le ribozyme VS. En considérant l’ensemble des données d’ITC et de RMN, l’étonnante stabilité de l’interaction I/V s’explique probablement par une combinaison de facteurs, dont les motifs U-turn, la présence d’un nucléotide exclu de la boucle de SLV (U700), la liaison de cations magnésium et l’empilement de bases continu à la jonction I/V.

Long-lived excited-state dynamics of i-motif structures probed by time-resolved infrared spectroscopy

UV-generated excited states of cytosine (C) nucleobases are precursors to mutagenic photoproduct formation. The i-motif formed from C-rich sequences is known to exhibit high yields of long-lived excited states following UV absorption. Here the excited states of several i-motif structures have been characterized following 267 nm laser excitation using time-resolved infrared spectroscopy (TRIR). All structures possess a long-lived excited state of ~300 ps and notably in some cases decays greater than 1 ns are observed. These unusually long-lived lifetimes are attributed to the interdigitated DNA structure which prevents direct base stacking overlap.

From Rational Design of Drug Crystals to Understanding of Nucleic Acid Structures: Lamivudine Duplex

A DNA-like duplex of nucleosides is probable to exist even without the 5`-phosphate groups needed to assemble the chain backbone. However, double-stranded helical structures of nucleosides are unknown. Here, we report a duplex of nucleoside analogs that is spontaneously assembled due to stacking of the neutral and protonated molecules of lamivudine, a nucleoside reverse transcriptase inhibitor (NTRI) widely used in anti-HIV drug combinatory medication. The left-handed lamivudine duplex has features similar to those of i-motif DNA, as the face-to-face base stacking and the helix rise per base pair. Furthermore, the protonation pattern on alternate bases expected for it DNA-like duplex stabilized by pairing of neutral and protonated cytosine fragments was observed for the first time in the lamivudine double-stranded helix. This structure demonstrates that hydrogen bonds can substitute for covalent phosphodiester linkage in the stabilization of the duplex backbone. This interesting example of spontaneous molecular self-organization indicates that the 5`-phosphate group could not be a requirement for duplex assembly.

Further Quantum Mechanical Evidence that Difluorotoluene does not Hydrogen Bond

Calculations were run on the methylated DNA base pairs adenine:thymine and adenine:difluorotoluene to further investigate the hydrogen-bonding properties of difluorotoluene (F). Geometries were optimized using hybrid density functional theory. Single-point calculations at the MP2(full) level were performed to obtain more rigorous energies. The functional counterpoise method was used to correct for the basis set superposition error (BSSE), and the interaction energies were also corrected for fragment relaxation. These corrections brought the B3LYP and MP2 interaction energies into excellent agreement. In the gas phase, the Gibbs free energies calculated at the B3LYP and MP2 levels of theory predict that A and T will spontaneously form an A:T pair while A:F spontaneously dissociates into A and F. Solvation effects on the pairing of the bases were explored using implicit solvent models for water and chloroform. In aqueous solution, both A:T and A:F are predicted to dissociate into their component monomers. Semiempirical calculations were performed on small sections of B-form DNA containing the two pairs, and the results provide support for the concept that base stacking is more important than hydrogen bonding for the stability of the A:F pair within a DNA helix.

Molecular Dynamics Simulation of a PNA·DNA·PNA Triple Helix in Aqueous Solution

Molecular dynamics simulations have been used to explore the conformational flexibility of a PNA·DNA·PNA triple helix in aqueous solution. Three 1.05 ns trajectories starting from different but reasonable conformations have been generated and analyzed in detail. All three trajectories converge within about 300 ps to produce stable and very similar conformational ensembles, which resemble the crystal structure conformation in many details. However, in contrast to the crystal structure, there is a tendency for the direct hydrogen-bonds observed between the amide hydrogens of the Hoogsteen-binding PNA strand and the phosphate oxygens of the DNA strand to be replaced by water-mediated hydrogen bonds, which also involve pyrimidine O2 atoms. This structural transition does not appear to weaken the triplex structure but alters groove widths and so may relate to the potential for recognition of such structures by other ligands (small molecules or proteins). Energetic analysis leads us to conclude that the reason that the hybrid PNA/DNA triplex has quite different helical characteristics from the all-DNA triplex is not because the additional flexibility imparted by the replacement of sugar−phosphate by PNA backbones allows motions to improve base-stacking but rather that base-stacking interactions are very similar in both types of triplex and the driving force comes from weak but definate conformational preferences of the PNA strands.

Warum Pentose- und nicht Hexose-Nucleinsäuren??. Teil V. (Purin-Purin)-Basenpaarung in der homo-DNS-Reihe: Guanin, Isoguanin, 2,6-Diaminopurin und Xanthin

Why Pentose- and Not Hexose-Nucleic Acids? Purine-Purine Pairing in homo-DNA: Guanine,Isoguanine, 2,6-Diaminopurine, and Xanthine This paper concludes the series of reports in this journal [1–4] on the chemistry of homo-DNA, the constitutionally simplifie dmodel system of hexopyranosyl-(6′ → 4′)-oligonucleotide systems stidued in our laboratory as potentially natural-nucleic-acid alternatives in the context of a chemical aetiology of nucleic-acid structure. The report describes the synthesis and pairing properties of homo-DNA oligonucleotides which contain as nucleobases exclusively purines, and gives, together with part III of the series [3], a survey of what we know today about purine-purine pairingin homo-DNA. In addition, the paper discusses those aspects of the chemistry of homo-DNA which, we think, influence the way how some of the structural features of DNA (and RNA) are to be interpreted on a qualitative level. Purine-purine pairing occurs in the homo-DNA domain in great variety. Most prominent is a novel tridentate Watson-Crick pair between guanine and isoguanine, as well as one between 2,6-diaminopurine and xanthinone, both giving rise to very stable duplexes containing the all-purine strands in antiparallel orientation. For the guanine-isoguanine pair, constitutional assignment is based on temperature-dependent UV and CD spectroscopy of various guanine- and isoguanine-containg duplexes in comparison with duplexes known to be paired in the reverse guanine is replaced by 7-carbauguanine. Isoguanine and 2,6-diaminopurine also have the capability of self-pariring in the reverse-Hoogsteen mode, as previously observed for adenine and guanine [3]. In this type of pairing, the interchangeably. Fig. 36 provides an overall survey of the relative strength of pairing in all possible purine-purine combinations. Watson-Crick pairing of isoguanine with guanine demands the former to participate in its 3H-tautomeric form; hitherto this specific tautomer had not been considered in the pairing chemistry of isoguanine. Whereas (cumulative) purine-purine pairing in DNA (reverse-Hoogsten or Hoogsteen) seems to occur in triplexes and tetrapalexes only, its occurrence in duplexes in a characteristic feature of homo-DNA chemistry. The occurrence of purine-purine Watson-Crick base pairs is probably a consequence of homo-DNA's quasi-linear ladder structure [1][4]. In a double helix, the distance between the two sugar C-atoms, on which a base pair is anchored, is expected to be constrained by the dimensions of the helix; in a linear duplex, however, there would be no restrictions with regard to base-pair length. Homo-DNA's ladder-like model also allows one to recognize one of the reasons why nucleic-acid duplexes prefer to pair in antiparallel, rather than parallel strand orientation: in homo-DNA duplexes, (averaged) backbone and base pair axes are strongly inclined toward one another [4]; the stronger this inclination, the higher the preference for antiparallel strand orientation is expected to be (Fig. 16). In retrospect, homo-DNA turns out to be one of the first artificial oligonucleotide systems (cf. Footnote 65) to demonstrate in a comprehensive way that informational base pairing involving purines and pyrimidines is not a capability unique to ribofuranosyl systems. Stability and helical shape of pairing complexes are not necessary conditions of one another; it is the potential for extensive conformational cooperativity of hte backbone structure with respect to the constellational demands of base pairing and base stacking that determines whether or nor a given type of base-carrying backbone structure is an informational pairing system. From the viewpoint of the chemical aetiology of nucleic-acid structure, which inspired our investigations on hexopyranosyl-(6′ → 4′)-oligonucleotide systems in the first place, the work on homo-DNA is only an extensive model study, because homo-DNA is not to be considered a potential natural-nucleic-acid altenratie. In retrospect, it seems fortunate that the model study was carried out, because without it we could hardly have comprehended the pairing behavior of the proper nucleic-acid alternatives which we have studied later and which will be discussed in Part VI of this series. The English footnotes to Fig. 1–49 provide an extension of this summary.

Structure of the histone mRNA hairpin required for cell cycle regulation of histone gene expression.

Expression of replication-dependent histone genes requires a conserved hairpin RNA element in the 3' untranslated regions of poly(A)-less histone mRNAs. The 3' hairpin element is recognized by the hairpin-binding protein or stem-loop-binding protein (HBP/SLBP). This protein-RNA interaction is important for the endonucleolytic cleavage generating the mature mRNA 3' end. The 3' hairpin and presumably HBP/SLBP are also required for nucleocytoplasmic transport, translation, and stability of histone mRNAs. RNA 3' processing and mRNA stability are both regulated during the cell cycle. Here, we have determined the three-dimensional structure of a 24-mer RNA comprising a mammalian histone RNA hairpin using heteronuclear multidimensional NMR spectroscopy. The hairpin adopts a novel UUUC tetraloop conformation that is stabilized by base stacking involving the first and third loop uridines and a closing U-A base pair, and by hydrogen bonding between the first and third uridines in the tetraloop. The HBP interaction of hairpin RNA variants was analyzed in band shift experiments. Particularly important interactions for HBP recognition are mediated by the closing U-A base pair and the first and third loop uridines, whose Watson-Crick functional groups are exposed towards the major groove of the RNA hairpin. The results obtained provide novel structural insight into the interaction of the histone 3' hairpin with HBP, and thus the regulation of histone mRNA metabolism.

Induction of topoisomerase I cleavage complexes by 1-β-d-arabinofuranosylcytosine (ara-C) in vitro and in ara-C-treated cells

1-β-d-Arabinofuranosylcytosine (Ara-C) is a nucleoside analog commonly used in the treatment of leukemias. Ara-C inhibits DNA polymerases and can be incorporated into DNA. Its mechanism of cytotoxicity is not fully understood. Using oligonucleotides and purified human topoisomerase I (top1), we found a 4- to 6-fold enhancement of top1 cleavage complexes when ara-C was incorporated at the +1 position (immediately 3′) relative to a unique top1 cleavage site. This enhancement was primarily due to a reversible inhibition of top1-mediated DNA religation. Because ara-C incorporation is known to alter base stacking and sugar puckering at the misincorporation site and at the neighboring base pairs, the observed inhibition of religation at the ara-C site suggests the importance of the alignment of the 5′-hydroxyl end for religation with the phosphate group of the top1 phosphotyrosine bond. This study also demonstrates that ara-C treatment and DNA incorporation trap top1 cleavage complexes in human leukemia cells. Finally, we report that camptothecin-resistant mouse P388/CPT45 cells with no detectable top1 are crossresistant to ara-C, which suggests that top1 poisoning is a potential mechanism for ara-C cytotoxicity.

A nicked duplex decamer DNA with a PEG6 tether

A dumbbell double-stranded DNA decamer tethered with a hexaethylene glycol linker moiety (DDSDPEG), with a nick in the centre of one strand, has been synthesised. The standard NMR methods, E.COSY, TOCSY, NOESY and HMQC, were used to measure 1H, 31P and T1 spectral parameters. Molecular modelling using rMD-simulated annealing was used to compute the structure. Scalar couplings and dipolar contacts show that the molecule adopts a right-handed B-DNA helix in 38 mM phosphate buffer at pH 7. Its high melting temperature confirms the good base stacking and stability of the duplex. This is partly attributed to the presence of the PEG6 linker at both ends of the duplex that restricts the dynamics of the stem pentamers and thus stabilises the oligonucleotide. The inspection of the global parameters shows that the linker does not distort the B-DNA geometry. The computed structure suggests that the presence of the nick is not disturbing the overall tertiary structure, base pair geometry or duplex base pairing to a substantial extent. The nick has, however, a noticeable impact on the local geometry at the nick site, indicated clearly by NMR analysis and reflected in the conformational parameters of the computed structure. The 1H spectra also show much sharper resonances in the presence of K+ indicating that conformational heterogeneity of DDSDPEG is reduced in the presence of potassium as compared to sodium or caesium ions. At the same time the 1H resonances have longer T1 times. This parameter is suggested as a sensitive gauge of stabilisation.

Inter-strand C-H...O hydrogen bonds stabilizing four-stranded intercalated molecules: stereoelectronic effects of O4' in cytosine-rich DNA.

DNA fragments with stretches of cytosine residues can fold into four-stranded structures in which two parallel duplexes, held together by hemiprotonated cytosine.cytosine+ (C.C+) base pairs, intercalate into each other with opposite polarity. The structural details of this intercalated DNA quadruplex have been assessed by solution NMR and single crystal x-ray diffraction studies of cytosine-rich sequences, including those present in metazoan telomeres. A conserved feature of these structures is the absence of stabilizing stacking interactions between the aromatic ring systems of adjacent C.C+ base pairs from intercalated duplexes. Effective stacking involves only the exocyclic keto groups and amino groups of the cytidine bases. The apparent absence of stability provided by stacking interactions between the bases in this intercalated DNA has prompted us to examine the available structures in detail, in particular with regard to unusual features that could compensate for the lack of base stacking. In addition to base-on-deoxyribose stacking and intra-cytidine C-H...O hydrogen bonds, this analysis reveals the presence of a hitherto unobserved, systematic intermolecular C-H...O hydrogen bonding network between the deoxyribose sugar moieties of antiparallel backbones in the four-stranded molecule.

A molecular anchor for stabilizing triple-helical DNA.

Molecular modeling has been used to predict that 2,6-disubstituted amidoanthraquinones, and not the 1,4 series, should preferentially interact with and stabilize triple-stranded DNA structures over duplex DNA. This is due to marked differences in the nature of chromophore-base stacking and groove accessibility for the two series. A DNA foot-printing method that monitors the extent of protection from DNase I cleavage on triplex formation has been used to examine the effects of a number of synthetic isomer compounds in the 1,4 and 2,6 series. The experimental results are in accord with the predicted behavior and confirm that the 1,4 series bind preferentially to double- rather than triple-stranded DNA, whereas the isomeric 2,6 derivatives markedly favor binding to triplex DNA.

DNA with Hydrophobic Base Substitutes: A Stable, Zipperlike Recognition Motif Based On Interstrand-Stacking Interactions

No Abstract

Synthesis, characterization, and structure of a new cobalt(II) Schiff-base complex with quinoxaline- 2-carboxalidine-2-amino-5-methylphenol

The mononuclear cobalt(II) complex [CoL2] H2O (where HL is quinoxaline-2-carboxalidine- 2-amino-5-methylphenol) has been prepared and characterized by elemental analysis, conductivity measurement, IR, UV-Vis spectroscopy, TG-DTA, and X-ray structure determination. The crystallographic study shows that cobalt(II) is distorted octahedral with each tridentate NNO Schiff base in a cis arrangement. The crystal exhibits a 2-D polymeric structure parallel to [010] plane, formed by O-H...N and O-H... O intermolecular hydrogen bonds and pye stacking interactions, as a racemic mixture of optical enantiomers. The ligand is a Schiff base derived from quinoxaline-2-carboxaldehyde