Are isolated nucleic acid bases really planar? A Car-Parrinello molecular dynamics study

Car-Parrinello molecular dynamics simulations of the flexibility of isolated DNA bases have been carried out. The comparison of lowest ring out-of-plane vibrations calculated by using MP2/cc-pvdz and BLYP/PW methods reveals that the DFT method with the plane wave basis set reasonably reproduces out-of-plane deformability of the pyrimidine ring in nucleic acid bases and could be used for reliable modeling of conformational flexibility of nucleobases. The conformational phase space of pyrimidine rings in thymine, cytosine, guanine, and adenine has been investigated by using the ab initio Car-Parrinello molecular dynamics method. It is demonstrated that all nucleic acid bases are highly flexible molecules and possess a nonplanar effective conformation of the pyrimidine ring despite the fact that the planar geometry corresponds to a minimum on the potential energy surface. The population of the planar geometry of the pyrimidine ring does not exceed 30%. Among the nonplanar conformations of the pyrimidine rings, the boat-like and half-chair conformations are the most populated.     

How flexible are fleximer nucleobases? A computational study

Density functional theory was used to study the potential energy surface for rotation about the carbon-carbon bonds in a variety of guanosine, adenosine, and inosine fleximers, which are modified purines with the imidazole and pyrimidine rings separated by a single carbon-carbon bond. Various connectivities between C4 or C5 of the imidazole ring and C5' or C6' of the pyrimidine ring were considered. Calculations on fleximer nucleobases in the absence of the ribose moiety suggest that a planar relative arrangement of the imidazole and pyrimidine rings is favored, and that all fleximers are indeed very flexible with regards to rotation about the carbon-carbon bond, where calculated barriers are generally less than 40 kJ mol(-1). Furthermore, calculated binding energies of fleximer-pyrimidine pairs indicate that the hydrogen-bonding properties of these modified nucleobases mimic those of the corresponding natural purine. Inclusion of the sugar moiety often leads to a favored nonplanar orientation of the two rings, and either a reduction in the rotational barrier height or small changes in the rotational surface depending on the connectivity and nucleobase considered. It is concluded that several connectivities may have favorable properties for biochemical applications where flexible nucleobases would be beneficial.     

Intrinsic lifetimes of the excited state of DNA and RNA bases

The lifetimes of the excited state of free nucleobases were measured in the gas phase for the first time. They are, respectively, 1.0 and 0.8 ps for the purine bases adenine (shown above) and guanine and 3.2, 2.4, and 6.4 ps for the pyrimidine bases cytosine, uracil, and thymine at 267 nm. The longer lifetimes of the pyrimidine bases may be associated with their higher propensity toward photodegradation, especially in the case of thymine. The ultrashort lifetime of nucleobases conventionally known in solution was found to be an intrinsic molecular property due to extremely facile internal conversion, and therefore the lifetime should be largely independent of the medium at this energy, that is, whether in vacuo, in solution, or in vivo. The evolutionary selection of nucleobases as the durable carriers of genetic information is suggested to be due to their inherent immunity from photochemical reactions.     

Effects of nucleobase metalation on frontier molecular orbitals: potential implications for pi-stacking interactions with tryptophan

Biochemical recognition processes mediated through pi-stacking interactions are a potential target for rational drug synthesis. A combination of electrostatic, hydrophobic, solvation, charge-transfer, induction, and dispersion interactions has been used to account for the three-dimensional arrangements observed in such motifs. A principal example involves the interaction of purine and pyrimidine rings of nucleic acids with aromatic amino-acid residues such as tryptophan, phenylalanine, and tyrosine. Protonation, alkylation, or coordination of a metal ion such as Pd(II) or Pt(II) to a nucleobase strengthens this interaction by lowering the energy of the lowest unoccupied molecular orbital (LUMO) of the modified nucleobase and improving overlap with the highest occupied molecular orbital (HOMO) in N-acetyl tryptophan. The relative energy difference between the frontier orbitals of isolated molecules, obtained using Density Functional Theory (DFT), is explored as a predictive tool for the strength of the pi-stacking interaction of the nucleobase/tryptophan pair. From the optimized structures of these species, evaluation of the donor-acceptor HOMO-LUMO gap (Deltaepsilon d-->a) suggests that this parameter is a promising predictor of pi-stacking strength for the donor-acceptor pairs presented in this study. The analysis correlates well with experimental association constants, measured by fluorescence spectroscopy, of metallated and alkylated nucleobases with tryptophan in comparison to free nucleobases.     

On the photoproduction of DNA/RNA cyclobutane pyrimidine dimers

The UV photoreactivity of different pyrimidine DNA/RNA nucleobases along the singlet manifold leading to the formation of cyclobutane pyrimidine dimers has been studied by using the CASPT2 level of theory. The initially irradiated singlet state promotes the formation of excimers between pairs of properly oriented nucleobases through the overlap between the ?? structures of two stacked nucleobases. The system evolves then to the formation of cyclobutane pyrimidine dimers via a shearing-type conical intersection activating a [2?+?2] photocycloaddition mechanism. The relative location of stable excimer conformations or alternative decay channels with respect to the reactive degeneracy region explains the differences in the photoproduction efficiency observed in the experiments for different nucleobases sequences. A comparative analysis of the main structural parameters and energetic profiles in the singlet manifold is carried out for thymine, uracil, cytosine, and 5-methylcytosine homodimers. Thymine and uracil dimers display the most favorable paths, in contrast to cytosine. Methylation of the nucleobases seems to increase the probability for dimerization.     

Double Functionalization of Carbon Nanotubes with Purine and Pyrimidine Derivatives

Herein, we have developed a synthetic strategy for the covalent double functionalization of single‐walled carbon nanotubes (SWCNTs) with a combination of purine–pyrimidine and purine–purine nucleobase systems. The nucleobases were introduced on the sidewall of oxidized SWCNTs through 1,3‐dipolar cycloaddition and by amidation of the carboxylic acids located at the tips and defect sites of the nanotubes. The new nanohybrids were characterized by transmission electron microscopy, thermogravimetric analysis, FTIR and Raman spectroscopy, magic‐angle spinning NMR spectroscopy, and Kaiser test. The nucleobase/SWCNT conjugates can be envisaged for the modulation of the interactions with nucleic acids by means of base pairing, thereby opening new possibilities in the development of DNA/CNT nanobioconjugates.     

Nucleotide recognition by protonated aminoglycosides

Interactions of protonated forms of kanamycin A with nucleotides and several simple phosphate anions have been studied by potentiometric and NMR titrations. The affinity of kanamycin A to anions is comparable to that observed with other aliphatic polyammonium receptors of similar charge, but it discriminates triphosphate nucleotides with different nucleobases with binding constants following the order GTP?CTP ≈ ATP. Kanamycin A also binds the respective uncharged nucleosides with the same selectivity. Binding of ATP is exothermic with a negative entropic contribution in contrast to what is expected for simple ion pairing. Other tested aminoglycosides, amikacin and streptomycin, bind ATP less efficiently than kanamycin A. Models of structures of kanamycin A complexes with ATP and GTP obtained by molecular mechanics (OPLS-2005) calculations based on ¹H and ³¹P NMR data confirm the possibility of nucleotide discrimination by simultaneous ion pairing of terminal nucleotide phosphate groups with ammonium sites of rings B and C and hydrogen bonding of the nucleobase at the ring A of the aminoglycoside.     

alpha-L-LNA (alpha-L-ribo configured locked nucleic acid) recognition of DNA: an NMR spectroscopic study

We have used NMR and CD spectroscopy to study and characterise two alpha-L-LNA:DNA duplexes, a nonamer that incorporates three alpha-L-LNA nucleotides and a decamer that incorporates four alpha-L-LNA nucleotides, in which alpha-L-LNA is alpha-L-ribo-configured locked nucleic acid. Both duplexes adopt right-handed helical conformations and form normal Watson-Crick base pairing with all nucleobases in the anti conformation. Deoxyribose conformations were determined from measurements of scalar coupling constants in the sugar rings, and for the decamer duplex, distance information was derived from 1H-1H NOE measurements. In general, the deoxyriboses in both of the alpha-L-LNA:DNA duplexes adopt S-type (B-type structure) sugar puckers, that is the inclusion of the modified alpha-L-LNA nucleotides does not perturb the local native B-like double-stranded DNA (dsDNA) structure. The CD spectra of the duplexes confirm these findings, as these display B-type characteristic features that allow us to characterise the overall duplex type as B-like. The 1H-1H NOE distances which were determined for the decamer duplex were employed in a simulated annealing protocol to generate a model structure for this duplex, thus allowing a more detailed inspection of the impact of the alpha-L-ribo-configured nucleotides. In this structure, it is evident that the malleable DNA backbone rearranges in the vicinity of the modified nucleotides in order to accommodate them and present their nucleobases in a geometry suitable for Watson-Crick base pairing.     

X-ray absorption spectroscopy of the nucleotide bases at the carbon, nitrogen, and oxygen K-edges

Near-edge X-ray absorption fine structure spectra of three pyrimidine (viz., cytosine, uracil, and thymine) and two purine (viz., adenine and guanine) nucleobases, which are the key constituents of DNA and RNA, were measured at the C, N, and O K-edges using the self-absorption-free partial electron yield mode. The nucleobase samples were prepared as highly pure native polycrystalline powder films. The spectra are analyzed in terms of the electronic structure of the nucleobases. Subtle chemical effects related to the molecular structures of these heterocyclic compounds with extended pi-electron systems are considered and discussed.     

Three-dimensional organization of helices: design principles for nucleobase-functionalized beta-peptides

The construction and molecular recognition of various three-dimensional biomimetic structures is based on the predictable de novo design of artificial molecules. In this regard beta-peptides are especially interesting, since stable secondary structures are obtained already with short sequences; one of them is the 14-helix in which every third residue has the same orientation. The covalent functionalization of every third 14-helix side chain with nucleobases was used for a reversible organization of two helices based on nucleobase pairing. A series of beta-peptides with various nucleobase sequences was synthesized and the stability of double strand formation was investigated. As few as four nucleobases are sufficient for considerable duplex stability. The stability of base pairing was examined by temperature-dependent UV spectroscopy and the formation of the 14-helix was confirmed by circular dichroism (CD) spectroscopy. The preferred strand orientation of complementary-nucleobase-modified beta-peptide helices was investigated as well as the influence of helix content on the duplex stability. The preorganization of a 14-helix in regard to double-strand recognition was tuned by the sequential order of polar beta-amino acids or by the amount of 2-aminocyclohexanecarboxylic acid units incorporated, which are known to facilitate 14-helix formation, respectively.     

Radiolysis and radioracemization of RNA ribonucleosides: implications for the origins of life

The RNA nucleosides, namely adenosine, cytidine, guanosine and uridine were γ-radioyzed in solid state and in vacuo at room temperature to a total dose of 3.2 MGy. Through electronic absorption spectroscopy, differential scanning calorimetry and through polarimetry and optical rotatory dispersion spectroscopy, it was found that the purine-based nucleosides (adenosine and guanosine) show a much higher radiolysis resistance than the pyrimidine-based nucleosides (cytidine and uridine). In an astrochemical/astrobiological context, these results may explain why purine nucleobases are found in practically all carbonaceous chondrite meteorites while the pyrimidine nucleobases are absent or below the detection limits of the current analytical techniques. In the hypothesis that both purines and pyrimidines nucleobases were present in certain bodies at the beginning of the solar system 4.6?×?10⁹ years ago, the radiolysis due to radionuclides decay has destroyed more easily and completely the pyrimidine bases due to their much lower radiolysis resistance than the purine bases.     

Cyclohexyl "base pairs" stabilize duplexes and intensify pyrene fluorescence by shielding it from natural base pairs

In this study, we investigated the stability and structure of artificial base pairs that contain cyclohexyl rings. The introduction of a single pair of isopropylcyclohexanes into the middle of DNA slightly destabilized the duplex. Interestingly, as the number of the "base pairs" increased, the duplex was remarkably stabilized. A duplex with six base pairs was even more stable than one containing six A-T pairs. Thermodynamic analysis revealed that changes in entropy and not enthalpy contributed to duplex stability, demonstrating that hydrophobic interactions between isopropyl groups facilitated the base pairing, and thus stabilized the duplex. NOESY of a duplex containing an isopropylcyclohexane-methylcyclohexane pair unambiguously demonstrated its "pairing" in the duplex because distinct NOEs between the protons of cyclohexyl moieties and imino protons of both of the neighboring natural base pairs were observed. CD spectra of duplexes tethering cyclohexyl moieties also showed a positive-negative couplet that is characteristic of the B-form DNA duplex. Taken together, these results showed that cyclohexyl moieties formed base pairs in the DNA duplex without severely disturbing the helical structure of natural DNA. Next, we introduced cyclohexyl base pairs between pyrene and nucleobases as an "insulator" that suppresses electron transfer between them. We found a massive increase in the quantum yield of pyrene due to the efficient shielding of pyrene from nucleobases. The cyclohexyl base pairs reported here have the potential to prepare highly fluorescent labeling agents by multiplying fluorophores and insulators alternately into DNA duplexes.     

Specific Purine‐Purine Base Pairing in Linear Alanyl‐Peptide Nucleic Acids

Purine‐purine base pairing with guanine, isoguanine, 2,6‐diaminopurine, and xanthine is investigated within the topology of alanyl‐PNA. Alanyl‐PNA is based on a regular peptide backbone with alternating configuration of the amino acids. The nucleobases are covalently linked as side chains. Their distance in peptides with β‐sheet conformation is similar to the favored base‐pair stacking distance. Therefore, alanyl‐PNA provides self‐pairing linear double‐strands. The linear double‐strand topology does not restrict base‐pair size and geometry. The favored base pairs are formed mostly dependent on recognition by H‐bonding. The synthesis of the nucleo‐amino acids with unnatural nucleobases and their oligomerization is described. Hexamers and a tetramer based on 2,6‐diaminopurine‐xanthine and guanine‐isoguanine base pairs were observed with very high stabilities. For xanthine‐xanthine self‐pairing, an unusual tridentate reverse Watson‐Crick pairing mode is indicated, that is only possible with xanthines pairing in different tautomeric forms. To investigate the nature of xanthine‐xanthine base pairs in more detail, quantum‐chemical calculations were performed. They establish the easier tautomerization of xanthine compared to uracil and indicate that, in the AO basis‐set limit, the tridentate pairing mode with mixed tautomers is favored.     

Theoretical molecular double-core-hole spectroscopy of nucleobases

Double-core-hole (DCH) spectra have been investigated for pyrimidine, purine, the RNA/DNA nucleobases, and formamide, using the density functional theory (DFT) method. DCH spectra of formamide were also examined by the complete-active-space self-consistent-field (CASSCF) method. All possible single- and two-site DCH (ssDCH and tsDCH) states of the nucleobases were calculated. The generalized relaxation energy and interatomic generalized relaxation energy were evaluated from the energy differences between ssDCH and single-core-hole (SCH) states and between tsDCH and SCH states, respectively. The generalized relaxation energy is correlated to natural bond orbital charge, whereas the interatomic generalized relaxation energy is correlated to the interatomic distance between the core holes at two sites. The present analysis using DCH spectroscopy demonstrates that the method is useful for the chemical analysis of large molecular systems.     

Direct Chemoselective Synthesis of N-3-Substituted Pyrimidinones in a Microwave-Assisted Method

Synthesis of selectively N-3-substituted pyrimidine nucleobases or pyrimidinones has always been a challenge because of poor regioselectivity and chemoselectivity. In this article we demonstrate a single-step, de novo synthesis of selectively N-3-substituted modified pyrimidinones. We have developed a microwave-assisted methodology for direct, chemoselective alkylation, benzylation, and arylation of C-5 and C-6 substituted pyrimidine nucleobases selectively at the N-3 position. The reactions were found to proceed, with high efficiency, without the requirement of solvent and were complete within 10–15 min of irradiation. The efficiency of the method was further improved by addition of a Lewis acid, which not only increases the yield significantly but also accelerates the reaction rate.     

Solvent- and environment-dependent fluorescence of modified nucleobases

Fluorescent nucleosides with modified nucleobases are useful tools for detecting nucleic acids and probing their structures and functions. Nucleobases are suitable for modification because 1) intrinsically light-absorbing nucleobases can be converted to fluorescent chromophores by simple chemical modification, 2) attaching substituents to nucleobases at appropriately selected positions does not inhibit base pairing or duplex formation, and 3) duplex formation and protein interactions affect the environment of nucleobases, causing changes in their fluorescence intensities and/or wavelengths. This review summarizes recent fluorescent nucleosides and their photophysical properties, such as absorption wavelength, emission wavelength, and fluorescence quantum yield together with their solvent dependency.     

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.     

On heteroaromaticity of nucleobases. Bond lengths as multidimensional phenomena

Three hundred and nine carbon-carbon, carbon-nitrogen, and carbon-oxygen pi-bond lengths in high precision crystal structures of 31 purine and pyrimidine nucleobases were related to the Pauling pi-bond order, its analogues corrected to crystal packing effects, the numbers of non-hydrogen atoms around the bond, and the sum of atomic numbers of the bond atoms. Principal Component Analysis (PCA) and Hierarchical Cluster Analysis (HCA) demonstrated that the bond lengths in the nucleobases are three-dimensional phenomenon, characterized by nine distinct classes of bonds. Bond lengths predicted by Linear Regression models, Pauling Harmonic Potential Curves, Multiple Linear Regression, Principal Component, and Partial Least Squares Regression were compared to those calculated by molecular mechanics, semiempirical, and ab initio methods using PCA-HCA procedure on the calculated bond lengths, statistical parameters, and structural aromaticity indices. Incorporation of crystal packing effects into bond orders makes multivariate models to be competitive to semiempirical results, while further improvement of quantum chemical calculations can be achieved by geometry optimization of molecular clusters.     

Mapping the sites for selective oxidation of guanines in DNA

The effective energy of a positive charge when it is localized at a specific guanine nucleobase in DNA was calculated using density functional theory. The results demonstrate that the efficiency of a guanine to act as a hole-trap in DNA strongly depends on the nature of the flanking nucleobases. The presence of a pyrimidine base at the 3' position adjacent to a guanine significantly increases the localization energy of the positive charge. The calculated distributions of a positive charge in sequences of two or three adjacent guanines, flanked by other nucleobases, provide an explanation for experimental literature data on the site-selective oxidation of DNA.     

A highly efficient protocol for the synthesis of N-aryl nucleobases using zinc oxide in ionic liquids

A highly efficient and simple method for the synthesis of N-aryl derivatives of pyrimidine and purine nucleobases via Narylation of nucleobases using zinc oxide in 1-butyl-3-methylimidazolium bromide ([bmim]Br) under microwave as well as thermal conditions is described. In both conditions, the title compounds were produced in high to excellent yields and in short reaction times.