Why does substitution of thymine by 6-ethynylpyridone increase the thermostability of DNA double helices?

Caging of a Strongly Pairing Fluorescent Thymidine Analog with Soft Nucleophiles

Controlling the pairing strength of nucleobases in DNA through reactions with compounds found inside the cell is a formidable challenge. Here we report how a thiazolyl substituent turns a strongly pairing ethynylpyridone C-nucleoside into a reactive residue in oligonucleotides. The thiazolyl-bearing pyridone reacts with soft nucleophiles, such as glutathione, but not with hard nucleophiles like hydroxide or carbonate. The addition products pair much more weakly with adenine in a complementary strand than the starting material, and also change their fluorescence. This makes oligonucleotides containing the new deoxynucleoside interesting for controlled release. Due to its reactivity toward N, P, S, and Se-nucleophiles, and the visual signal accompanying chemical conversion, the fluorescent nucleotide reported here may also have applications in chemical biology, sensing and diagnostics.

High Fidelity Enzyme-Free Primer Extension with an Ethynylpyridone Thymidine Analog

High fidelity base pairing is important for the transmission of genetic information. Weak base pairs can lower fidelity, complicating sequencing, amplification and replication of DNA. Thymidine 5'-monophosphate (TMP) is the most weakly pairing nucleotide among the canonical deoxynucleotides, causing high errors rates in enzyme-free primer extension. Here we report the synthesis of an ethynylpyridone C-nucleoside analog of 3'-amino-2',3'-dideoxythymidine monophosphate and its incorporation in a growing strand by enzyme-free primer extension. The ethynylpyridone C-nucleotide accelerates extension more than five-fold, reduces misincorporation and readily displaces TMP in competition experiments. The results bode well for the use of the C-nucleoside as replacements for thymidine in practical applications.

Replacing thymine with a strongly pairing fifth Base: A combined quantum mechanics and molecular dynamics study

The non-natural ethynylmethylpyridone C-nucleoside (W), a thymidine (T) analogue that can be incorporated in oligonucleotides by automated synthesis, has recently been reported to form a high fidelity base pair with adenosine (A) and to be well accommodated in B-DNA duplexes. The enhanced binding affinity for A of W, as compared to T, makes it an ideal modification for biotechnological applications, such as efficient probe hybridization for the parallel detection of multiple DNA strands. In order to complement the experimental study and rationalize the impact of the non-natural W nucleoside on the structure, stability and dynamics of DNA structures, we performed quantum mechanics (QM) calculations along with molecular dynamics (MD) simulations. Consistently with the experimental study, our QM calculations show that the A:W base pair has an increased stability as compared to the natural A:T pair, due to an additional CH-π interaction. Furthermore, we show that mispairing between W and guanine (G) causes a distortion in the planarity of the base pair, thus explaining the destabilization of DNA duplexes featuring a G:W pair. MD simulations show that incorporation of single or multiple consecutive A:W pairs in DNA duplexes causes minor changes to the intra- and inter-base geometrical parameters, while a moderate widening/shrinking of the major/minor groove of the duplexes is observed. QM calculations applied to selected stacks from the MD simulations also show an increased stacking energy for W, over T, with the neighboring bases.

An AZT Analog with Strongly Pairing Ethynylpyridone Nucleobase and Its Antiviral Activity against HSV1

Challenges resulting from novel viruses or new strains of known viruses call for new antiviral agents. Nucleoside analogs that act as inhibitors of viral polymerases are an attractive class of antivirals. For nucleosides containing thymine, base pairing is weak, making it desirable to identify nucleobase analogs that pair more strongly with adenine, in order to compete successfully with the natural substrate. We have recently described a new class of strongly binding thymidine analogs that contain an ethynylmethylpyridone as base and a C-nucleosidic linkage to the deoxyribose. Here we report the synthesis of the 3'-azido-2',3'-deoxyribose derivative of this compound, dubbed AZW, both as free nucleoside and as ProTide phosphoramidate. As a proof of principle, we studied the activity against Herpes simplex virus type 1 (HSV1). Whereas the ProTide phosphoramidate suffered from low solubility, the free nucleoside showed a stronger inhibitory effect than that of AZT in a plaque reduction assay. This suggests that strongly pairing C-nucleoside analogs of pyrimidines have the potential to become active pharmaceutical ingredients with antiviral activity.

Estimating Strengths of Individual Hydrogen Bonds in RNA Base Pairs: Toward a Consensus between Different Computational Approaches

Noncoding RNA molecules are composed of a large variety of noncanonical base pairs that shape up their functionally competent folded structures. Each base pair is composed of at least two interbase hydrogen bonds (H-bonds). It is expected that the characteristic geometry and stability of different noncanonical base pairs are determined collectively by the properties of these interbase H-bonds. We have studied the ground-state electronic properties [using density functional theory (DFT) and DFT-D3-based methods] of all the 118 normal base pairs and 36 modified base pairs, belonging to 12 different geometric families (cis and trans of WW, WH, HH, WS, HS, and SS) that occur in a nonredundant set of high-resolution RNA crystal structures. Having addressed some of the limitations of the earlier approaches, we provide here a comprehensive compilation of the average energies of different types of interbase H-bonds (E _HB). We have also characterized each interbase H-bond using 13 different parameters that describe its geometry, charge distribution at its bond critical point (BCP), and n → σ*-type charge transfer from filled π orbitals of the H-bond acceptor to the empty antibonding orbital of the H-bond donor. On the basis of the extent of their linear correlation with the H-bonding energy, we have shortlisted five parameters to model linear equations for predicting E _HB values. They are (i) electron density at the BCP: ρ, (ii) its Laplacian: ∇²ρ, (iii) stabilization energy due to n → σ*-type charge transfer: E(2), (iv) donor-hydrogen distance, and (v) hydrogen-acceptor distance. We have performed single variable and multivariable linear regression analysis over the normal base pairs and have modeled sets of linear relationships between these five parameters and E _HB. Performance testing of our model over the set of modified base pairs shows promising results, at least for the moderately strong H-bonds.

Theoretical characterization of the conformational features of unnatural oligonucleotides containing a six nucleotide genetic alphabet

The addition of the unnatural P:Z base pair to the four naturally occurring DNA bases expands the genetic alphabet and yields an artificially expanded genetic information system (AEGIS). Herein, the structural feature of oligonucleotides containing a novel unnatural P:Z base pair is characterized using both molecular dynamics and quantum chemistry. The results show that the incorporation of the novel artificial base pair (P:Z) preserves the global conformational feature of duplex DNA except for some local structures. The Z-nitro group imparts new properties to the groove width, which widens the major groove. The unnatural oligonucleotides containing mismatched base pairs exhibit low stability. This ensures efficient and high-fidelity replication. In general, the incorporation of the P:Z pair strengthens the stability of the corresponding DNA duplex. The calculated results also show that the thermostability originates from both hydrogen interaction and stacking interaction. The Z-nitro group plays an important role in enhancing the stability of the H-bonds and stacking strength of the P:Z pair. Overall, the present results provide theoretical insights in the exploration of artificially expanded genetic information systems.

Consequences of Mg2+ binding on the geometry and stability of RNA base pairs

Quantum chemical calculations reveal the role of magnesium in stabilizing the geometries of intrinsically unstable RNA base pairs.

How Does Mg2+ Modulate the RNA Folding Mechanism: A Case Study of the G:C W:W Trans Basepair

Reverse Watson-Crick G:C basepairs (G:C W:W Trans) occur frequently in different functional RNAs. This is one of the few basepairs whose gas-phase-optimized isolated geometry is inconsistent with the corresponding experimental geometry. Several earlier studies indicate that through post-transcriptional modification, direct protonation, or coordination with Mg²⁺, accumulation of positive charge near N7 of guanine can stabilize the experimental geometry. Interestingly, recent studies reveal significant variation in the position of putatively bound Mg²⁺. This, in conjunction with recently raised doubts regarding some of the Mg²⁺ assignments near the imino nitrogen of guanine, is suggestive of the existence of multiple Mg²⁺ binding modes for this basepair. Our detailed investigation of Mg²⁺-bound G:C W:W Trans pairs occurring in high-resolution RNA crystal structures shows that they are found in 14 different contexts, eight of which display Mg²⁺ binding at the Hoogsteen edge of guanine. Further examination of occurrences in these eight contexts led to the characterization of three different Mg²⁺ binding modes: 1) direct binding via N7 coordination, 2) direct binding via O6 coordination, and 3) binding via hydrogen-bonding interaction with the first-shell water molecules. In the crystal structures, the latter two modes are associated with a buckled and propeller-twisted geometry of the basepair. Interestingly, respective optimized geometries of these different Mg²⁺ binding modes (optimized using six different DFT functionals) are consistent with their corresponding experimental geometries. Subsequent interaction energy calculations at the MP2 level, and decomposition of its components, suggest that for G:C W:W Trans , Mg²⁺ binding can fine tune the basepair geometries without compromising with their stability. Our results, therefore, underline the importance of the mode of binding of Mg²⁺ ions in shaping RNA structure, folding and function.

Computational studies on the interactions of glycine amino acid with graphene, h-BN and h-SiC monolayers

In the present work, the adsorption of glycine amino acid and its zwitterionic form onto three different hexagonal sheets, namely graphene, boron-nitride (h-BN) and silicon carbide (h-SiC), has been investigated within the framework of density functional theory (DFT) calculations.

Feasibility of occurrence of different types of protonated base pairs in RNA: a quantum chemical study

Protonated nucleobases have significant roles in facilitating catalytic functions of RNA, and in stabilizing different structural motifs. Reported pKa values of nucleobase protonation suggest that the population of neutral nucleobases is 10(3)-10(4) times higher than that of protonated nucleobases under physiological conditions (pH ∼ 7.4). Therefore, a molecular level understanding of various putative roles of protonated nucleobases cannot be achieved without addressing the question of how their occurrence propensities and stabilities are related to the free energy costs associated with the process of protonation under physiological conditions. With water as the proton donor, we use advanced QM methods to evaluate the site specific protonation propensities of nucleobases in terms of their associated free energy changes (ΔGprot). Quantitative follow up on the energetics of base pair formation and database search for evaluating their occurrence frequencies, reveal a lack of correlation between base pair stability and occurrence propensities on the one hand, and ease of protonation on the other. For example, although N7 protonated adenine (ΔGprot = 40.0 kcal mol(-1)) is found to participate in stable base pairing, base pairs involving N7 protonated guanine (ΔGprot = 36.8 kcal mol(-1)), on geometry optimization, converge to a minima where guanine transfers its extra proton to its partner base. Such observations, along with examples of weak base pairs involving N3 protonation of cytosine (ΔGprot = 37.0 kcal mol(-1)) are rationalized by analysing the protonation induced charge redistributions which are found to significantly influence, both positively and negatively, the hydrogen bonding potentials of different functional sites of individual nucleobases. Protonation induced charge redistribution is also found to strongly influence (i) the aromatic character of the rings of the participating bases and (ii) hydrogen bonding potential of the free edges of the protonated base pair. Comprehensive analysis of a non-redundant RNA crystal structure dataset further reveals that, while availability of stabilization possibilities determine the feasibility of occurrence of protonated bases, their occurrence context and specific functional roles are important factors determining their occurrence propensities.

The role of N7 protonation of guanine in determining the structure, stability and function of RNA base pairs

Ab initio computations and bioinformatics studies reveal that stabilization of some important RNA structural motifs might involve N7 protonation of guanine.