Rational design of hetero-ring-expanded guanine analogs with enhanced properties for modified DNA building blocks

The properties and modes of recognition of physiological DNAs associated with the four natural nucleobases might be extended, in principle, by the design of non-natural nucleobase derivatives. The goal is an expansion of the genetic alphabet, with the possible outcome of producing new DNAs with improved physical or biological properties. In this work, a new series of hetero-ring-expanded guanine analogs are proposed, and their relevant structural characteristics and electronic properties are determined by density functional theory. The stabilities of the decamer DNA duplexes (dn.dC)10 (where n represents the corresponding expanded guanine analog designed here) are also examined, using molecular dynamics. The simulations show that the designed motifs can form stable DNA-like structures. We determined the pairing energies for the Watson-Crick (WC) hydrogen-bonded dimers between the expanded G-analogs and the natural C, and found that the pairing energies are close to those of the natural GC pair. The calculated adiabatic ionization potentials (IPs) of the size-expanded guanine analogs and their base pairs, and the corresponding vertical ionization potentials, show that some are distinctly smaller than the corresponding natural versions. The HOMO-LUMO energy gaps for most of the size-expanded guanine analogs and their WC base pairs are considerably lower than those of the corresponding natural base and base pairs. Thus, the expanded G bases may be considered as DNA genetic motifs, and they may serve as building blocks for potential biological applications and the development of molecular electronic devices.     

How do size-expanded DNA nucleobases enhance duplex stability? Computational analysis of the hydrogen-bonding and stacking ability of xDNA bases

Computational chemistry (B3LYP, MP2) is used to study the properties of size-expanded DNA nucleobases generated by inserting a benzene spacer into the natural nucleobases. Although the addition of the spacer does not significantly affect the hydrogen-bonding properties of natural nucleobases, the orientation of the base about the glycosidic bond necessary for Watson-Crick binding is destabilized, which could have implications for the selectivity of expanded bases, as well as the stability of expanded duplexes. Consideration of the (stacked) binding energies in the preferred relative orientation of natural and expanded nucleobases aligned according to their centers of mass reveals that the stacking within natural dimers can be increased by up to 50% upon expansion of one nucleobase and up to 90% upon expansion of two nucleobases. The implications of these findings to the stability of expanded duplexes were revealed by considering simplified models of natural and mixed duplexes composed of four nucleobases. Although intra- and interstrand interactions within double helices are typically less than those predicted when nucleobases are stacked according to their centers of mass, some nucleobases utilize their full stacking potential within double helices, where both intra- and interstrand interactions can be significant. Most importantly, increasing the size of nucleobases within the duplex significantly increases both intra- and interstrand stacking interactions. Specifically, some interactions are double the magnitude of the corresponding intrastrand interactions in natural helices, and even greater increases in interstrand interactions are sometimes found. Thus, our work suggests that mixed duplexes composed of natural bases hydrogen bound to expanded bases may exploit the increase in the inherent stacking ability of the expanded bases in more than one way and thereby afford duplexes with greater stability than natural DNA.     

Novel benzopyrimidines as widened analogues of DNA bases

We report on the synthesis, stacking, and pairing properties of a new structural class of size-expanded pyrimidine nucleosides, abbreviated dyT and dyC. Their bases are benzo-homologated variants of thymine and cytosine and have a design that is distinct from a previously described class of size-expanded (xDNA) pyrimidines, with a different vector of expansion relative to the sugar. We term this new base geometry "yDNA" (a mnemonic for "wide DNA"). Both C-glycosides were prepared using Pd-mediated coupling of iodinated base derivatives with a deoxyribose precursor. As free deoxynucleosides, both dyT and dyC displayed robust fluorescence, with emission maxima at 375 and 390 nm, respectively. Both widened pyrimidines could be incorporated readily as protected phosphoramidite derivatives into synthetic oligonucleotides. Experiments in "dangling end" DNA contexts revealed that both yT and yC stack more favorably than their natural counterparts. When opposite natural bases in the context of Watson-Crick DNA were paired, the yT nucleotide formed a pair with A that was equally stable as a T-A pair, despite the mismatch in size with the neighboring natural pairs. The yC nucleotide (paired opposite G) was destabilizing by a small amount in the same context. Despite the large size of the pairs, both yT and yC were selective for their Watson-Crick complementary partners A and G, respectively. The pairing properties and fluorescence of yDNA nucleotides may lead to useful applications in the study of steric effects in DNA-protein interactions. In addition, the compounds may serve as building blocks for a large-sized artificial genetic system.     

Exploring the limits of DNA size: naphtho-homologated DNA bases and pairs

A new design for DNA bases and base pairs is described in which the pyrimidine bases are widened by naphtho-homologation. Two naphtho-homologated deoxyribosides, dyyT (1) and dyyC (2), were synthesized and could be incorporated into oligonucleotides as suitably protected phosphoramidite derivatives. The deoxyribosides were found to be fluorescent, with emission maxima at 446 and 433 nm, respectively. Studies with single substitutions of 1 and 2 in the natural DNA context revealed exceptionally strong base stacking propensity for both. Sequences containing multiple substitutions of 1 and 2 paired opposite adenine and guanine were subsequently mixed and studied by several analytical methods. Data from UV mixing experiments, FRET measurements, fluorescence quenching experiments, and hybridizations on beads suggest that complementary "doublewide DNA" (yyDNA) strands may self-assemble into helical complexes with 1:1 stoichiometry. Data from thermal denaturation plots and CD spectra were less conclusive. Control experiments in one sequence context gave evidence that yyDNA helices, if formed, are preferentially antiparallel and are sequence selective. Hypothesized base pairing schemes are analogous to Watson-Crick pairing, but with glycosidic C1'-C1' distances widened by over 45%, to ca. 15.2 A. The possible self-assembly of the double-wide DNA helix establishes a new limit for the size of information-encoding, DNA-like molecules, and the fluorescence of yyDNA bases suggests uses as reporters in monomeric and oligomeric forms.     

yDNA versus xDNA pyrimidine nucleobases: computational evidence for dependence of duplex stability on spacer location

The structural and binding properties of the natural and x- and y-pyrimidines were compared using computational methods. Our calculations show that although the x-pyrimidines favor different orientations about the glycosidic bond compared to the natural pyrimidines, which could have implications for the formation and resulting stability of xDNA duplexes and jeopardize the selectivity of expanded nucleobases, y-pyrimidines have rotational profiles more similar to the natural bases. Increasing the pyrimidine size using a benzene spacer leads to relatively minor changes in the hydrogen-bond strength of isolated Watson-Crick base pairs. However, differences in the anomeric carbon distances in pairs composed of x- or y-pyrimidines suggest yDNA may yield a more optimal expanded structure. By stacking two monomers via their centers of mass, we find that the expanded nucleobases stack much stronger than the natural bases. Additionally, although replacing xT by yT changes the stacking energy by less than 5 kJ mol (-1), replacing xC by yC significantly strengthens complexes with the natural nucleobases (by up to 30%). Calculations on larger duplex models composed of four nucleobases reveal that x- and y-pyrimidines can increase duplex stability of natural helices by strengthening both the intra and interstrand stacking interactions. Furthermore, when the total stability (sum of all hydrogen-bonding and (intrastrand and interstrand) stacking interactions) of the larger models is considered, y-pyrimidines lead to more stable complexes than x-pyrimidines for all but three duplex sequences. Thus, through analysis of a variety of properties, our calculations suggest that the location of the benzene spacer affects the properties of expanded nucleobases and the stability of expanded duplexes, and therefore should be carefully considered when designing future expanded analogues.     

Structural, dynamical, and electronic transport properties of modified DNA duplexes containing size-expanded nucleobases

Among the distinct strategies proposed to expand the genetic alphabet, size-expanded nucleobases are promising for the development of modified DNA duplexes with improved biotechnological properties. In particular, duplexes built up by replacing canonical bases with the corresponding benzo-fused counterparts could be valuable as molecular nanowires. In this context, this study reports the results of classical molecular dynamics simulations carried out to examine the structural and dynamical features of size-expanded DNAs, including both hybrid duplexes containing mixed pairs of natural and benzo-fused bases (xDNA) and pure size-expanded (xxDNA) duplexes. Furthermore, the electronic structure of both natural and size-expanded duplexes is examined by means of density functional computations. The results confirm that the structural and flexibility properties of the canonical DNA are globally little affected by the presence of benzo-fused bases. The most relevant differences are found in the enhanced size of the grooves, and the reduction in the twist. However, the analysis also reveals subtle structural effects related to the nature and sequence of benzo-fused bases in the duplex. On the other hand, electronic structure calculations performed for xxDNAs confirm the reduction in the HOMO-LUMO gap predicted from the analysis of the natural bases and their size-expanded counterparts, which suggests that pure size-expanded DNAs can be good conductors. A more complex situation is found for xDNAs, where fluctuations in the electrostatic interaction between base pairs exerts a decisive influence on the modulation of the energy gap.     

Cis Watson-Crick/Watson-Crick Base Pairs in DNA: A Density Functional Theory(DFT) Study

The four nucleic acid DNA bases(adenine, thymine, guanine, cytosine) and ten cis Watson-Crick/Watson-Crick(cis WC/WC) DNA base pairs were investigated by density functional theory(DFT) quantum chemical calculations. Geometry optimizations were carried out on the four bases and ten base pairs at the B3LYP level with 6-31G~(**) basis set. All the optimizations were performed within Cs symmetry. The optimum structures for the four bases and seven cis WC/WC base pairs were obtained, and Natural Bond Orbital analysis(NBO) was based on these structures. The possibilities of matches between any two of the four bases through their Watson-Crick(WC) edges were discussed. The structures of seven cis WC/WC base pairs change to a certain extent relative to these of the four bases due to the formation of hydrogen bonds. These base pairs existing in DNA have an important influence on the structural stability of the double helix. The analysis of the electronic structures and molecular orbitals for seven cis WC/WC base pairs can provide significant information about the relationship between charge transfer along the hydrogen bond and the Frontier orbitals of these base pairs.     

A computational proposal for the experimentally observed discriminatory behavior of hypoxanthine, a weak universal nucleobase

A computational model composed of six nucleobases was used to investigate why hypoxanthine does not yield duplexes of equal stability when paired opposite each of the natural DNA nucleobases. The magnitudes of all nearest-neighbor interactions in a DNA helix were calculated, including hydrogen-bonding, intra- and interstrand stacking interactions, as well as 1-3 intrastrand stacking interactions. Although the stacking interactions in DNA relevant arrangements are significant and account for at least one third of the total stabilization energy in our nucleobase complexes, the trends in the magnitude of the stacking interactions cannot explain the relative experimental melting temperatures previously reported in the literature. Furthermore, although the total hydrogen-bonding interactions explain why hypoxanthine preferentially pairs with cytosine, the experimental trend for the remaining nucleobases (A, T, G) is not explained. In fact, the calculated pairing preference of hypoxanthine matches that determined experimentally only when the sum of all types of nearest-neighbor interactions is considered. This finding highlights a strong correlation between the relative magnitude of the total nucleobase-nucleobase interactions and measured melting temperatures for DNA strands containing hypoxanthine despite the potential role of other factors (including hydration, temperature, sugar-phosphate backbone). By considering a large range of sequence combinations, we reveal that the binding preference of hypoxanthine is strongly dependent on the nucleobase sequence, which may explain the varied ability of hypoxanthine to universally bind to the natural bases. As a result, we propose that future work should closely examine the interplay between the dominant nucleobase-nucleobase interactions and the overall strand stability to fully understand how sequence context affects the universal binding properties of modified bases and to aid the design of new molecules with ambiguous pairing properties.     

Ab initio study of naphtho-homologated DNA bases

Naphtho-homologated DNA bases have been recently used to build a new type of size-expanded DNA known as yyDNA. We have used theoretical techniques to investigate the structure, tautomeric preferences, base-pairing ability, stacking interactions, and HOMO-LUMO gaps of the naphtho-bases. The structure of these bases is found to be similar to that of the benzo-fused predecessors (y-bases) with respect to the planarity of the aromatic rings and amino groups. Tautomeric studies reveal that the canonical-like forms of naphtho-thymine (yyT) and naphtho-adenine (yyA) are the most stable tautomers, leading to hydrogen-bonded dimers with the corresponding natural nucleobases that mimic the Watson-Crick pairing. However, the canonical-like species of naphtho-guanine (yyG) and naphtho-cytosine (yyC) are not the most stable tautomers, and the most favorable hydrogen-bonded dimers involve wobble-like pairings. The expanded size of the naphtho-bases leads to stacking interactions notably larger than those found for the natural bases, and they should presumably play a dominant contribution in modulating the structure of yyDNA duplexes. Finally, the HOMO-LUMO gap of the naphtho-bases is smaller than that of their benzo-base counterparts, indicating that size-expansion of DNA bases is an efficient way of reducing their HOMO-LUMO gap. These results are examined in light of the available experimental evidence reported for yyT and yyC.     

A new four-base genetic helix, yDNA, composed of widened benzopyrimidine-purine pairs

We describe the properties of stable DNA-like self-assembled helices composed entirely of base pairs involving two new size-expanded pyrimidines. We term this new helix geometry "yDNA" (an abbreviation of "wide DNA"). The new pyrimidine analogues, yT and yC, are increased in size by benzo-homologation and have a geometry that is distinct from previous size-expanded pyrimidines. The yT and yC deoxyribosides were incorporated into oligodeoxynucleotides designed to form four pairs: yT-A, A-yT, yC-G, and G-yC. Helices were characterized by thermal denaturation, mixing data, and circular dichroism spectra. Results showed that highly stable double-stranded helices were formed in several sequence contexts. The data further showed that yT and yC could be segregated onto one strand and used to bind to natural strands of DNA with high sequence selectivity. The combination of yC, yT, G, and A make up a new selective, self-assembling four-base genetic pairing system that functions in many respects like natural DNA, but which is structurally distinct. The results establish that multiple variants of size-expanded DNA-like helices are feasible and suggest the possibility of a future eight-base genetic system based on the yDNA geometry. Finally, the high binding selectivity, affinity, and fluorescence of yDNA strands may yield useful applications in detection of nucleic acid sequences.     

Fluorous base-pairing effects in a DNA polymerase active site

We describe selective "fluorous" effects in the active site of a DNA polymerase, by using nucleotide analogues whose pairing edges are perfluorinated. The 5'-triphosphate deoxynucleotide derivatives of DNA base analogues 2,3,4,5-tetrafluorobenzene ((F)B) and 4,5,6,7-tetrafluoroindole ((F)I), as well as hydrocarbon controls benzene (B) and indole (I), were synthesized and studied as substrates for the DNA Polymerase I Klenow fragment (KF exo-). Modified nucleotides were present in the DNA template or were supplied as nucleoside triphosphates in studies of the steady-state kinetics of single nucleotide insertion. When supplied opposite the non-natural bases in the template strand, the hydrophobic nucleoside triphosphates were incorporated by up to two orders of magnitude more efficiently than the natural deoxynucleoside triphosphates. The purine-like fluorinated indole nucleotide ((F)I) was the most efficiently inserted of the four hydrophobic analogues, with the most effective incorporation occurring opposite the pyrimidine-like tetrafluorobenzene ((F)B). In all cases, the polyfluorinated base pairs were more efficiently processed than the analogous hydrocarbon pairs. A preliminary test of polymerase extension beyond these pairs showed that only the (F)B base is appreciably extended; the inefficient extension is consistent with recently published data regarding other nonpolar base pairs. These results suggest the importance of hydrophobicity, stacking, and steric interactions in the polymerase-mediated replication of DNA base pairs that lack hydrogen bonds. These findings further suggest that the enhanced hydrophobicity of polyfluoroaromatic bases could be employed in the design of new, selective base pairs that are orthogonal to the natural Watson-Crick pairs used in replication.     

Toward a designed, functioning genetic system with expanded-size base pairs: solution structure of the eight-base xDNA double helix

We describe the NMR-derived solution structure of the double-helical form of a designed eight-base genetic pairing system, termed xDNA. The benzo-homologous xDNA design contains base pairs that are wider than natural DNA pairs by ca. 2.4 A (the width of a benzene ring). The eight component bases of this xDNA helix are A, C, G, T, xA, xT, xC, and xG. The structure was solved in aqueous buffer using 1D and 2D NMR methods combined with restrained molecular dynamics. The data show that the decamer duplex is right-handed and antiparallel, and hydrogen-bonded in a way analogous to that of Watson-Crick DNA. The sugar-phosphate backbone adopts a regular conformation similar to that of B-form DNA, with small dihedral adjustments due to the larger circumference of the helix. The grooves are much wider and more shallow than those of B-form DNA, and the helix turn is slower, with ca. 12 base pairs per 360 degrees turn. There is an extensive intra- and interstrand base stacking surface area, providing an explanation for the greater stability of xDNA relative to natural DNA. There is also evidence for greater motion in this structure compared to a previous two-base-expanded helix; possible chemical and structural reasons for this are discussed. The results confirm paired self-assembly of the designed xDNA system. This suggests the possibility that other genetic system structures besides the natural one might be functional in encoding information and transferring it to new complementary strands.     

Size-expanded analogues of dG and dC: synthesis and pairing properties in DNA

We describe the completion of the set of four benzo-fused expanded DNA (xDNA) nucleoside analogues. We previously reported the development of benzo-fused analogues of dA and dT and their inclusion in an exceptionally stable new four-base genetic system, termed xDNA, in which the base pairs were expanded in size. Here we describe the preparation and properties of the second half of this nucleotide set: namely, the previously unknown dxC and dxG nucleosides. The dxC analogue was prepared from the previously reported dxT nucleoside in three steps and 57% yield. The large-sized deoxyguanosine analogue was prepared from an intermediate in the synthesis of dxA, yielding dxG in 14 steps overall (2.4%). Suitably protected versions of the deoxynucleosides were prepared for oligonucleotide synthesis following standard procedures, and they were readily incorporated into DNA by automated synthesizer. "Dangling-end" measurements revealed that the benzo-fused homologues stack considerably more strongly on neighboring DNA sequences than do their natural counterparts. Base pairing experiments with xC or xG bases showed that they pair selectively with their Watson-Crick partners, but with mild destabilization, due apparently to their larger size. Overall, the data suggest that the fluorescent xG and xC bases may be useful probes of steric effects in the study of biological nucleotide recognition mechanisms. In addition, the completion of the xDNA nucleoside set makes it possible in the future to construct full four-base xDNA strands that can target any sequence of natural DNA and RNA.     

On the aromatic character of the heterocyclic bases of DNA and RNA

Studies based on ab initio optimized geometries (at B3LYP/6-311+G** and MP2/6-311+G** levels) and on experimental structures retrieved from the Cambridge Structural Database (CSD) reveal that the nucleobases constituting DNA and RNA differ significantly in their aromatic character, as shown by the geometry-based index of aromaticity HOMA that ranges from 0.466 for thymine to 0.917 for adenine, based on B3LYP/6-311+G** calculations, and 0.495-0.926, respectively, if based on the MP2/6-311+G** level. Aromaticity of the bases decreases markedly with an increase of the number of double-bond C=X (X = N, O) substituents at the rings. H-bonds involving C=O groups in Watson-Crick pairs cause an increase of the aromatic character of the rings.     

Photophysical characters of rationally designed hetero-ring-expanded guanine analogues and effect of cytosine pairing

We present the results of the CIS and TDB3LYP calculations of the optical absorption and emission spectra of some newly designed guanine (G) analogues and their Watson-Crick base pairs. Compared with natural G, the onset absorption peaks of these newly designed analogues are red-shifted, while the fluorescence peaks are blue-shifted. In general, the first excited singlet states (pipi*) of these analogues are nonplanar for all bases considered here. But, the Stokes shifts for the designed G-analogues are much smaller than that of natural G, suggesting that they have stronger molecular rigidity and higher fluorescence quantum yields than those of natural G. The first excited states of these Watson-Crick base pairs essentially originate from those of their isolated purine moieties, as demonstrated from the S1 geometries of their Watson-Crick base pairs. For G and its analogues, A1 and A2 (they are ring-expanded with one-bond intercalation at the C5 site), the pairing with cytosine reduces the oscillator strengths of both the first absorption peak (by 27%-60%) and the fluorescent emission (by 19%-23%), while for the analogues A3, A4, and xG in which G is ring-expanded with a two-bond intercalation at the C5 site, the pairing, in contrast, increases the oscillator strengths of both the first absorption peak (by 11%-15%) and the fluorescent emission (by 3%-20%). These observations indicate that the pairing with cytosine can quench the fluorescence for G, A1, and A2 but enhance the fluorescence quantum yields for A3, A4, and xG. The significant shifts induced by ring-expansion in the ring-expanded G with a two-bond intercalation at the C5 site reveal a possibility for their fluorescent detections.     

A series of nonpolar thymidine analogues of increasing size: DNA base pairing and stacking properties

[reaction: see text] We describe the properties in DNA of a set of five nonpolar nucleoside mimics in which shape is similar but size is increased gradually. The compounds vary in the size of their exocyclic substituents, which range from hydrogen to iodine, and are designed to test the steric effects of nucleosides, nucleotides, and DNA in biological systems in a systematic way. We describe the conversion of toluene, 2,4-difluorotoluene, 2,4-dichlorotoluene, 2,4-dibromotoluene, and 2,4-diiodotoluene deoxyribosides into suitably protected phosphoramidite derivatives and their incorporation into synthetic DNAs. Studies of their behavior in the context of hexamer and dodecamer duplexes were carried out, with comparison to natural thymine. Thermal melting data with compounds in 5' dangling positions showed that all five compounds stack more strongly than thymine, and all the dihalo-substituted cases stack more strongly than the unsubstituted toluene case. Stacking correlated with surface area and hydrophobicity, both of which increase across the series. In base-pairing studies, all five compounds showed destabilized pairing opposite natural bases (relative to thymine-adenine pairing), as expected. Notably, pairing among the nonpolar base analogues was considerably more stable, and some of the pairs involving the largest analogues showed stability equal to that of a natural thymine-adenine pair. The results establish the base pairing properties of a potentially useful new series of biochemical probes for DNA-protein interactions and also identify a set of new, stable hydrophobic base pairs for designed genetic pairing systems.     

Quadracyclic Adenine: A Non‐Perturbing Fluorescent Adenine Analogue

Fluorescent‐base analogues (FBAs) comprise a group of increasingly important molecules for the investigation of nucleic acid structure and dynamics as well as of interactions between nucleic acids and other molecules. Here, we report on the synthesis, detailed spectroscopic characterisation and base‐pairing properties of a new environment‐sensitive fluorescent adenine analogue, quadracyclic adenine (qA). After developing an efficient route of synthesis for the phosphoramidite of qA it was incorporated into DNA in high yield by using standard solid‐phase synthesis procedures. In DNA qA serves as an adenine analogue that preserves the B‐form and, in contrast to most currently available FBAs, maintains or even increases the stability of the duplex. We demonstrate that, unlike fluorescent adenine analogues, such as the most commonly used one, 2‐aminopurine, and the recently developed triazole adenine, qA shows highly specific base‐pairing with thymine. Moreover, qA has an absorption band outside the absorption of the natural nucleobases (>300 nm) and can thus be selectively excited. Upon excitation the qA monomer displays a fluorescence quantum yield of 6.8 % with an emission maximum at 456 nm. More importantly, upon incorporation into DNA the fluorescence of qA is significantly less quenched than most FBAs. This results in quantum yields that in some sequences reach values that are up to fourfold higher than maximum values reported for 2‐aminopurine. To facilitate future utilisation of qA in biochemical and biophysical studies we investigated its fluorescence properties in greater detail and resolved its absorption band outside the DNA absorption region into distinct transition dipole moments. In conclusion, the unique combination of properties of qA make it a promising alternative to current fluorescent adenine analogues for future detailed studies of nucleic acid‐containing systems.     

Autonomously pairing cysteinyl-linked nucleotide analogues with a unique architecture

We report the efficient pairing in water of the first representative of oligonucleotide analogues in which the backbone is replaced by linking elements between the nucleobases. The architecture of the new analogue demonstrates that the structural differentiation of oligonucleotides into a contiguous backbone and nucleobases, as embodied by the natural nucleic acids and all nucleotide analogues analyzed to date, is not a prerequisite for pairing. UV and circular dichroism analyses of self-complementary and non-self-complementary octanucleotide analogues strongly suggest the fully reversible, sequence-specific association of our new analogues to form a left-handed double helix with an antiparallel strand orientation that is characterized by melting temperatures and free enthalpies higher than those of natural RNA and DNA of the same sequence. The linking element incorporates an L-cysteine moiety that allows a short and efficient synthesis of the monomeric building blocks and, through the choice of either L- or D-cysteine, gives access to either one of the enantiomeric oligomers and thus to left- or right-handed helices.