Measurement of nucleobase pKa values in model mononucleotides shows RNA-RNA duplexes to be more stable than DNA-DNA duplexes

Chemical Synthesis of Modified RNA

Ribonucleic acids (RNAs) play a vital role in living organisms. Many of their cellular functions depend critically on chemical modification. Methods to modify RNA in a controlled manner-both in vitro and in vivo-are thus essential to evaluate and understand RNA biology at the molecular and mechanistic levels. The diversity of modifications, combined with the size and uniformity of RNA (made up of only 4 nucleotides) makes its site-specific modification a challenging task that needs to be addressed by complementary approaches. One such approach is solid-phase RNA synthesis. We discuss recent developments in this field, starting with new protection concepts in the ongoing effort to overcome current size limitations. We continue with selected modifications that have posed significant challenges for their incorporation into RNA. These include deazapurine bases required for atomic mutagenesis to elucidate mechanistic aspects of catalytic RNAs, and RNA containing xanthosine, N4-acetylcytidine, 5-hydroxymethylcytidine, 3-methylcytidine, 2'-OCF3, and 2'-N3 ribose modifications. We also discuss the all-chemical synthesis of 5'-capped mRNAs and the enzymatic ligation of chemically synthesized oligoribonucleotides to obtain long RNA with multiple distinct modifications, such as those needed for single-molecule FRET studies. Finally, we highlight promising developments in RNA-catalyzed RNA modification using cofactors that transfer bioorthogonal functionalities.

NMR measurements of transient low-populated tautomeric and anionic Watson-Crick-like G·T/U in RNA:DNA hybrids: implications for the fidelity of transcription and CRISPR/Cas9 gene editing

Many biochemical processes use the Watson-Crick geometry to distinguish correct from incorrect base pairing. However, on rare occasions, mismatches such as G·T/U can transiently adopt Watson-Crick-like conformations through tautomerization or ionization of the bases, giving rise to replicative and translational errors. The propensities to form Watson-Crick-like mismatches in RNA:DNA hybrids remain unknown, making it unclear whether they can also contribute to errors during processes such as transcription and CRISPR/Cas editing. Here, using NMR R1ρ experiments, we show that dG·rU and dT·rG mismatches in two RNA:DNA hybrids transiently form tautomeric (Genol·T/U $ \mathbin{\lower.3ex\hbox{$\buildrel\textstyle\rightarrow\over {\smash{\leftarrow}\vphantom{_{\vbox to.5ex{\vss}}}}$}}$ G·Tenol/Uenol) and anionic (G·T-/U-) Watson-Crick-like conformations. The tautomerization dynamics were like those measured in A-RNA and B-DNA duplexes. However, anionic dG·rU- formed with a ten-fold higher propensity relative to dT-·rG and dG·dT- and this could be attributed to the lower pKa (ΔpKa ∼0.4-0.9) of U versus T. Our findings suggest plausible roles for Watson-Crick-like G·T/U mismatches in transcriptional errors and CRISPR/Cas9 off-target gene editing, uncover a crucial difference between the chemical dynamics of G·U versus G·T, and indicate that anionic Watson-Crick-like G·U- could play a significant role evading Watson-Crick fidelity checkpoints in RNA:DNA hybrids and RNA duplexes.

NMR measurements of transient low-populated tautomeric and anionic Watson-Crick-like G·T/U in RNA:DNA hybrids: Implications for the fidelity of transcription and CRISPR/Cas9 gene editing

Many biochemical processes use the Watson-Crick geometry to distinguish correct from incorrect base pairing. However, on rare occasions, mismatches such as G•T/U can transiently adopt Watson-Crick-like conformations through tautomerization or ionization of the bases, giving rise to replicative and translational errors. The propensities to form Watson-Crick-like mismatches in RNA:DNA hybrids remain unknown, making it unclear whether they can also contribute to errors during processes such as transcription and CRISPR/Cas editing. Here, using NMR R 1ρ experiments, we show that dG•rU and dT•rG mismatches in two RNA:DNA hybrids transiently form tautomeric (G enol •T/U ⇄G•T enol /U enol ) and anionic (G•T - /U - ) Watson-Crick-like conformations. The tautomerization dynamics were like those measured in A-RNA and B-DNA duplexes. However, anionic dG•rU - formed with a ten-fold higher propensity relative to dT - •rG and dG•dT - and this could be attributed to the lower pK a (Δ pK a ∼0.4-0.9) of U versus T. Our findings suggest plausible roles for Watson-Crick-like G•T/U mismatches in transcriptional errors and CRISPR/Cas9 off-target gene editing, uncover a crucial difference between the chemical dynamics of G•U versus G•T, and indicate that anionic Watson-Crick-like G•U - could play a significant role evading Watson-Crick fidelity checkpoints in RNA:DNA hybrids and RNA duplexes.

Calculated pKa Values for a Series of Aza- and Deaza-Modified Nucleobases

A variety of synthetic modified nucleobases have been used to investigate the structure and function of RNA and DNA or act as enzyme inhibitors. A set of these modifications involves the addition or removal of a nitrogen atom in the ring. These aza and deaza modifications have garnered interest as useful biochemical tools, but information on some of their physical characteristics is lacking. In this study, the B3LYP density functional with the 6-31+G(d,p) basis set and an implicit-explicit solvent model was used to perform ab initio quantum mechanical studies to estimate pK_a values of aza- and deaza-modified nucleobases. A comparison between theoretical and known experimental pK_a values was carried out, and adjustment factors were applied to 57 pK_a values in the purine and pyrimidine data sets.

1-Deazaguanosine-Modified RNA: The Missing Piece for Functional RNA Atomic Mutagenesis

Atomic mutagenesis is the key to advance our understanding of RNA recognition and RNA catalysis. To this end, deazanucleosides are utilized to evaluate the participation of specific atoms in these processes. One of the remaining challenges is access to RNA-containing 1-deazaguanosine (c¹G). Here, we present the synthesis of this nucleoside and its phosphoramidite, allowing first time access to c¹G-modified RNA. Thermodynamic analyses revealed the base pairing parameters for c¹G-modified RNA. Furthermore, by NMR spectroscopy, a c¹G-triggered switch of Watson-Crick into Hoogsteen pairing in HIV-2 TAR RNA was identified. Additionally, using X-ray structure analysis, a guanine–phosphate backbone interaction affecting RNA fold stability was characterized, and finally, the critical impact of an active-site guanine in twister ribozyme on the phosphodiester cleavage was revealed. Taken together, our study lays the synthetic basis for c¹G-modified RNA and demonstrates the power of the completed deazanucleoside toolbox for RNA atomic mutagenesis needed to achieve in-depth understanding of RNA recognition and catalysis.

Calculations of pKa Values for a Series of Naturally Occurring Modified Nucleobases

Modified nucleobases are found in functionally important regions of RNA and are often responsible for essential structural roles. Many of these nucleobase modifications are dynamically regulated in nature, with each modification having a different biological role in RNA. Despite the high abundance of modifications, many of their characteristics are still poorly understood. One important property of a nucleobase is its pK_a value, which has been widely studied for unmodified nucleobases, but not for the modified versions. In this study, the pK_a values of modified nucleobases were determined by performing ab initio quantum mechanical calculations using a B3LYP density functional with the 6-31+G(d,p) basis set and a combination of implicit-explicit solvation systems. This method, which was previously employed to determine the pK_a values of unmodified nucleobases, is applicable to a variety of modified nucleobases. Comparisons of the pK_a values of modified nucleobases give insight into their structural and energetic impacts within nucleic acids.

Acid-base and lipophilic properties of peptide nucleic acid derivatives

The first combined experimental and theoretical study on the ionization and lipophilic properties of peptide nucleic acid (PNA) derivatives, including eleven PNA monomers and two PNA decamers, is described. The acidity constants (pK_a) of individual acidic and basic centers of PNA monomers were measured by automated potentiometric pH titrations in water/methanol solution, and these values were found to be in agreement with those obtained by MoKa software. These results indicate that single nucleobases do not change their pK_a values when included in PNA monomers and oligomers. In addition, immobilized artificial membrane chromatography was employed to evaluate the lipophilic properties of PNA monomers and oligomers, which showed the PNA derivatives had poor affinity towards membrane phospholipids, and confirmed their scarce cell penetrating ability. Overall, our study not only is of potential relevance to evaluate the pharmacokinetic properties of PNA, but also constitutes a reliable basis to properly modify PNA to obtain mimics with enhanced cell penetration properties.

•

The first study on acid-base and lipophilic properties of peptide nucleic acids (PNA).

•

pK_a of acid-base centers of PNA evaluated by potentiometric method and MoKa prediction.

•

NMR experiments provide additional information on the protonation of PNA monomers.

•

Lipophilicity of PNA monomers and oligomers is investigated by IAM chromatography.

•

This study can lay the basis of evaluating the pharmacokinetic properties of PNA.

Impact of 3-deazapurine nucleobases on RNA properties

Deazapurine nucleosides such as 3-deazaadenosine (c3A) are crucial for atomic mutagenesis studies of functional RNAs. They were the key for our current mechanistic understanding of ribosomal peptide bond formation and of phosphodiester cleavage in recently discovered small ribozymes, such as twister and pistol RNAs. Here, we present a comprehensive study on the impact of c3A and the thus far underinvestigated 3-deazaguanosine (c3G) on RNA properties. We found that these nucleosides can decrease thermodynamic stability of base pairing to a significant extent. The effects are much more pronounced for 3-deazapurine nucleosides compared to their constitutional isomers of 7-deazapurine nucleosides (c7G, c7A). We furthermore investigated base pair opening dynamics by solution NMR spectroscopy and revealed significantly enhanced imino proton exchange rates. Additionally, we solved the X-ray structure of a c3A-modified RNA and visualized the hydration pattern of the minor groove. Importantly, the characteristic water molecule that is hydrogen-bonded to the purine N3 atom and always observed in a natural double helix is lacking in the 3-deazapurine-modified counterpart. Both, the findings by NMR and X-ray crystallographic methods hence provide a rationale for the reduced pairing strength. Taken together, our comparative study is a first major step towards a comprehensive understanding of this important class of nucleoside modifications.

The effect of adenine protonation on RNA phosphodiester backbone bond cleavage elucidated by deaza-nucleobase modifications and mass spectrometry

The catalytic strategies of small self-cleaving ribozymes often involve interactions between nucleobases and the ribonucleic acid (RNA) backbone. Here we show that multiply protonated, gaseous RNA has an intrinsic preference for the formation of ionic hydrogen bonds between adenine protonated at N3 and the phosphodiester backbone moiety on its 5'-side that facilitates preferential phosphodiester backbone bond cleavage upon vibrational excitation by low-energy collisionally activated dissociation. Removal of the basic N3 site by deaza-modification of adenine was found to abrogate preferential phosphodiester backbone bond cleavage. No such effects were observed for N1 or N7 of adenine. Importantly, we found that the pH of the solution used for generation of the multiply protonated, gaseous RNA ions by electrospray ionization affects phosphodiester backbone bond cleavage next to adenine, which implies that the protonation patterns in solution are at least in part preserved during and after transfer into the gas phase. Our study suggests that interactions between protonated adenine and phosphodiester moieties of RNA may play a more important mechanistic role in biological processes than considered until now.

General Recognition of U-G, U-A, and C-G Pairs by Double-Stranded RNA-Binding PNAs Incorporated with an Artificial Nucleobase

Chemically modified peptide nucleic acids (PNAs) show great promise in the recognition of RNA duplexes by major-groove PNA·RNA-RNA triplex formation. Triplex formation is favored for RNA duplexes with a purine tract within one of the RNA duplex strands, and is severely destabilized if the purine tract is interrupted by pyrimidine residues. Here, we report the synthesis of a PNA monomer incorporated with an artificial nucleobase S, followed by the binding studies of a series of S-modified PNAs. Our data suggest that an S residue incorporated into short 8-mer dsRNA-binding PNAs (dbPNAs) can recognize internal Watson-Crick C-G and U-A, and wobble U-G base pairs (but not G-C, A-U, and G-U pairs) in RNA duplexes. The short S-modified PNAs show no appreciable binding to DNA duplexes or single-stranded RNAs. Interestingly, replacement of the C residue in an S·C-G triple with a 5-methyl C results in the disruption of the triplex, probably due to a steric clash between S and 5-methyl C. Previously reported PNA E base shows recognition of U-A and A-U pairs, but not a U-G pair. Thus, S-modified dbPNAs may be uniquely useful for the general recognition of RNA U-G, U-A, and C-G pairs. Shortening the succinyl linker of our PNA S monomer by one carbon atom to have a malonyl linker causes a severe destabilization of triplex formation. Our experimental and modeling data indicate that part of the succinyl moiety in a PNA S monomer may serve to expand the S base forming stacking interactions with adjacent PNA bases.

The influence of substituents and the environment on the NMR shielding constants of supramolecular complexes based on A-T and A-U base pairs

In the present study, we have theoretically analyzed supramolecular complexes based on the Watson-Crick A-T and A-U base pairs using dispersion-corrected density functional theory (DFT). Hydrogen atoms H8 and/or H6 in the natural adenine and thymine/uracil bases were replaced, respectively, by substituents X8, Y6 = NH^-, NH₂, NH₃⁺ (N series), O^-, OH, OH₂⁺ (O series), F, Cl or Br (halogen series). We examined the effect of the substituents on the hydrogen-bond lengths, strength and bonding mechanism, and the NMR shielding constants of the C2-adenine and C2-thymine/uracil atoms in the base pairs. The general belief in the literature that there is a direct connection between changes in the hydrogen-bond strength and the C2-adenine shielding constant is conclusively rejected by our computations.

How RNase HI (Escherichia coli) promoted site-selective hydrolysis works on RNA in duplex with carba-LNA and LNA substituted antisense strands in an antisense strategy context?

A detailed kinetic study of 36 single modified AON-RNA heteroduplexes shows that substitution of a single native nucleotide in the antisense strand (AON) by locked nucleic acid (LNA) or by diastereomerically pure carba-LNA results in site-dependent modulation of RNase H promoted cleavage of complementary mRNA strands by 2 to 5 fold at 5'-GpN-3' cleavage sites, giving up to 70% of the RNA cleavage products. The experiments have been performed using RNase HI of Escherichia coli. The 2nd best cleavage site, being the 5'-ApN-3' sites, cleaves up to 23%, depending upon the substitution site in 36 isosequential complementary AONs. A comparison of the modified AON promoted RNA cleavage rates with that of the native AON shows that sequence-specificity is considerably enhanced as a result of modification. Clearly, relatively weaker 5'-purine (Pu)-pyrimidine (Py)-3' stacking in the complementary RNA strand is preferred (giving ∼90% of total cleavage products), which plays an important role in RNase H promoted RNA cleavage. A plausible mechanism of RNase H mediated cleavage of the RNA has been proposed to be two-fold, dictated by the balancing effect of the aromatic character of the purine aglycone: first, the locally formed 9-guanylate ion (pK_a 9.3, ∼18-20% N1 ionized at pH 8) alters the adjoining sugar-phosphate backbone around the scissile phosphate, transforming its sugar N/S conformational equilibrium, to preferential S-type, causing preferential cleavage at 5'-GpN-3' sites around the center of 20 mer complementary mRNA. Second, the weaker nearest-neighbor strength of 5'-Pu-p-Py-3' stacking promotes preferential 5'-GpN-3' and 5'-ApN-3' cleavage, providing ∼90% of the total products, compared to ∼50% in that of the native one, because of the cLNA/LNA substituent effect on the neighboring 5'-Pu-p-Py-3' sites, providing both local steric flexibility and additional hydration. This facilitates both the water and water/Mg²⁺ ion availability at the cleavage site causing sequence-specific hydrolysis of the phosphodiester bond of scissile phosphate. The enhancement of the total rate of cleavage of the complementary mRNA strand by up to 25%, presented in this work, provides opportunities to engineer a single modification site in appropriately substituted AONs to design an effective antisense strategy based on the nucleolytic stability of the AON strand versus RNase H capability to cleave the complementary RNA strand.

On the mechanism of RNA phosphodiester backbone cleavage in the absence of solvent

Ribonucleic acid (RNA) modifications play an important role in the regulation of gene expression and the development of RNA-based therapeutics, but their identification, localization and relative quantitation by conventional biochemical methods can be quite challenging. As a promising alternative, mass spectrometry (MS) based approaches that involve RNA dissociation in ‘top-down’ strategies are currently being developed. For this purpose, it is essential to understand the dissociation mechanisms of unmodified and posttranscriptionally or synthetically modified RNA. Here, we have studied the effect of select nucleobase, ribose and backbone modifications on phosphodiester bond cleavage in collisionally activated dissociation (CAD) of positively and negatively charged RNA. We found that CAD of RNA is a stepwise reaction that is facilitated by, but does not require, the presence of positive charge. Preferred backbone cleavage next to adenosine and guanosine in CAD of (M+nH)ⁿ⁺ and (M−nH)ⁿ⁻ ions, respectively, is based on hydrogen bonding between nucleobase and phosphodiester moieties. Moreover, CAD of RNA involves an intermediate that is sufficiently stable to survive extension of the RNA structure and intramolecular proton redistribution according to simple Coulombic repulsion prior to backbone cleavage into c and y ions from phosphodiester bond cleavage.

CH-π hydrogen bonds in biological macromolecules

This is a sequel to the previous Perspective "The CH-π hydrogen bond in chemistry. Conformation, supramolecules, optical resolution and interactions involving carbohydrates", which featured in a PCCP themed issue on "Weak Hydrogen Bonds - Strong Effects?": Phys. Chem. Chem. Phys., 2011, 13, 13873-13900. Evidence that weak hydrogen bonds play an enormously important role in chemistry and biochemistry has now accumulated to an extent that the rigid classical concept of hydrogen bonds formulated by Pauling needs to be seriously revised and extended. The concept of a more generalized hydrogen bond definition is indispensable for understanding the folding mechanisms of proteins. The CH-π hydrogen bond, a weak molecular force occurring between a soft acid CH and a soft base π-electron system, among all is one of the most important and plays a functional role in defining the conformation and stability of 3D structures as well as in many molecular recognition events. This concept is also valuable in structure-based drug design efforts. Despite their frequent occurrence in organic molecules and bio-molecules, the importance of CH-π hydrogen bonds is still largely unknown to many chemists and biochemists. Here we present a review that deals with the evidence, nature, characteristics and consequences of the CH-π hydrogen bond in biological macromolecules (proteins, nucleic acids, lipids and polysaccharides). It is hoped that the present Perspective will show the importance of CH-π hydrogen bonds and stimulate interest in the interactions of biological macromolecules, one of the most fascinating fields in bioorganic chemistry. Implication of this concept is enormous and valuable in the scientific community.

Calculations of pKa's and redox potentials of nucleobases with explicit waters and polarizable continuum solvation

The SMD implicit solvation model augmented with one and four explicit water molecules was used to calculate pKa's and redox potentials of N-methyl-substituted nucleic acid bases guanine, adenine, cytosine, thymine, and uracil. Calculations were carried out with the B3LYP/6-31+G(d,p) level of theory. The same numbers of water molecules were hydrogen bonded to the neutral, protonated, and deprotonated nucleobases in their unoxidized and oxidized forms. The improvement in pKa1 involving neutrals and cations was modest. By contrast, the improvement in pKa2 involving neutrals and anions was quite significant, reducing the mean absolute error from 4.6 pKa units with no waters, to 2.6 with one water and 1.7 with four waters. For the oxidation of nucleobases, adding explicit waters did little to improve E(X(•),H(+)/XH), possibly because both species in the redox couple are neutral molecules at pH 7.

Incorporation of thio-pseudoisocytosine into triplex-forming peptide nucleic acids for enhanced recognition of RNA duplexes

Peptide nucleic acids (PNAs) have been developed for applications in biotechnology and therapeutics. There is great potential in the development of chemically modified PNAs or other triplex-forming ligands that selectively bind to RNA duplexes, but not single-stranded regions, at near-physiological conditions. Here, we report on a convenient synthesis route to a modified PNA monomer, thio-pseudoisocytosine (L), and binding studies of PNAs incorporating the monomer L. Thermal melting and gel electrophoresis studies reveal that L-incorporated 8-mer PNAs have superior affinity and specificity in recognizing the duplex region of a model RNA hairpin to form a pyrimidine motif major-groove RNA₂–PNA triplex, without appreciable binding to single-stranded regions to form an RNA–PNA duplex or, via strand invasion, forming an RNA–PNA₂ triplex at near-physiological buffer condition. In addition, an L-incorporated 8-mer PNA shows essentially no binding to single-stranded or double-stranded DNA. Furthermore, an L-modified 6-mer PNA, but not pseudoisocytosine (J) modified or unmodified PNA, binds to the HIV-1 programmed −1 ribosomal frameshift stimulatory RNA hairpin at near-physiological buffer conditions. The stabilization of an RNA₂–PNA triplex by L modification is facilitated by enhanced van der Waals contacts, base stacking, hydrogen bonding and reduced dehydration energy. The destabilization of RNA–PNA and DNA–PNA duplexes by L modification is due to the steric clash and loss of two hydrogen bonds in a Watson–Crick-like G–L pair. An RNA₂–PNA triplex is significantly more stable than a DNA₂–PNA triplex, probably because the RNA duplex major groove provides geometry compatibility and favorable backbone–backbone interactions with PNA. Thus, L-modified triplex-forming PNAs may be utilized for sequence-specifically targeting duplex regions in RNAs for biological and therapeutic applications.

Ultra-sensitive pH control of supramolecular polymers and hydrogels: pKa matching of biomimetic monomers

pK_a-matched monomers self-assemble into supramolecular polymers at pH 7, producing hydrogels with high pH sensitivity.

Factors influencing the electrokinetic injection of oligonucleotides in capillary gel electrophoresis when using laser-induced fluorescence detection

Capillary gel electrophoresis (CGE) is a powerful tool for the analysis of oligonucleotides owing to its extraordinary resolving power. However, the only feasible injection mode for CGE, electrokinetic injection, can cause bias of the injected amount and thus reproducibility issues for CGE methods. Although the source of the bias in electrokinetic injection for analysis of small molecules by capillary zone electrophoresis has long been identified, there are very few studies on electrokinetic injection issues for biological molecules analyzed by CGE. In this study, we report three issues related to electrokinetic injection for oligonucleotides. First, the relationship between the injection amount and the sample solution resistance is not always linear for oligonucleotides, as has been observed for small molecules. Second, the injecting water prior to an oligonucleotide sample dramatically improves the reproducibility of both the injected amount and resolution through a 'stacking-like' mechanism. Third, optimizing the gel concentration dramatically increases the amount of oligonucleotide that is injected into the column.

Intrinsic acid-base properties of a hexa-2'-deoxynucleoside pentaphosphate, d(ApGpGpCpCpT): neighboring effects and isomeric equilibria

The intrinsic acid-base properties of the hexa-2'-deoxynucleoside pentaphosphate, d(ApGpGpCpCpT) [=(A1∙G2∙G3∙C4∙C5∙T6)=(HNPP)⁵⁻] have been determined by ¹H NMR shift experiments. The pKa values of the individual sites of the adenosine (A), guanosine (G), cytidine (C), and thymidine (T) residues were measured in water under single-strand conditions (i.e., 10% D₂O, 47 °C, I=0.1 M, NaClO₄). These results quantify the release of H⁺ from the two (N7)H⁺ (G∙G), the two (N3)H⁺ (C∙C), and the (N1)H⁺ (A) units, as well as from the two (N1)H (G∙G) and the (N3)H (T) sites. Based on measurements with 2'-deoxynucleosides at 25 °C and 47 °C, they were transferred to pKa values valid in water at 25 °C and I=0.1 M. Intramolecular stacks between the nucleobases A1 and G2 as well as most likely also between G2 and G3 are formed. For HNPP three pKa clusters occur, that is those encompassing the pKa values of 2.44, 2.97, and 3.71 of G2(N7)H⁺, G3(N7)H⁺, and A1(N1)H⁺, respectively, with overlapping buffer regions. The tautomer populations were estimated, giving for the release of a single proton from five-fold protonated H₅(HNPP)(±) , the tautomers (G2)N7, (G3)N7, and (A1)N1 with formation degrees of about 74, 22, and 4%, respectively. Tautomer distributions reveal pathways for proton-donating as well as for proton-accepting reactions both being expected to be fast and to occur practically at no "cost". The eight pKa values for H₅(HNPP)(±) are compared with data for nucleosides and nucleotides, revealing that the nucleoside residues are in part affected very differently by their neighbors. In addition, the intrinsic acidity constants for the RNA derivative r(A1∙G2∙G3∙C4∙C5∙U6), where U=uridine, were calculated. Finally, the effect of metal ions on the pKa values of nucleobase sites is briefly discussed because in this way deprotonation reactions can easily be shifted to the physiological pH range.

Quantum-mechanical analysis of the energetic contributions to π stacking in nucleic acids versus rise, twist, and slide

Symmetry-adapted perturbation theory (SAPT) is applied to pairs of hydrogen-bonded nucleobases to obtain the energetic components of base stacking (electrostatic, exchange-repulsion, induction/polarization, and London dispersion interactions) and how they vary as a function of the helical parameters Rise, Twist, and Slide. Computed average values of Rise and Twist agree well with experimental data for B-form DNA from the Nucleic Acids Database, even though the model computations omitted the backbone atoms (suggesting that the backbone in B-form DNA is compatible with having the bases adopt their ideal stacking geometries). London dispersion forces are the most important attractive component in base stacking, followed by electrostatic interactions. At values of Rise typical of those in DNA (3.36 Å), the electrostatic contribution is nearly always attractive, providing further evidence for the importance of charge-penetration effects in π-π interactions (a term neglected in classical force fields). Comparison of the computed stacking energies with those from model complexes made of the "parent" nucleobases purine and 2-pyrimidone indicates that chemical substituents in DNA and RNA account for 20-40% of the base-stacking energy. A lack of correspondence between the SAPT results and experiment for Slide in RNA base-pair steps suggests that the backbone plays a larger role in determining stacking geometries in RNA than in B-form DNA. In comparisons of base-pair steps with thymine versus uracil, the thymine methyl group tends to enhance the strength of the stacking interaction through a combination of dispersion and electrosatic interactions.

Base-pairing preferences, physicochemical properties and mutational behaviour of the DNA lesion 8-nitroguanine

8-Nitro-2'-deoxyguanosine (8-nitrodG) is a relatively unstable, mutagenic lesion of DNA that is increasingly believed to be associated with tissue inflammation. Due to the lability of the glycosidic bond, 8-nitrodG cannot be incorporated into oligodeoxynucleotides (ODNs) by chemical DNA synthesis and thus very little is known about its physicochemical properties and base-pairing preferences. Here we describe the synthesis of 8-nitro-2'-O-methylguanosine, a ribonucleoside analogue of this lesion, which is sufficiently stable to be incorporated into ODNs. Physicochemical studies demonstrated that 8-nitro-2'-O-methylguanosine adopts a syn conformation about the glycosidic bond; thermal melting studies and molecular modelling suggest a relatively stable syn-8-nitroG·anti-G base pair. Interestingly, when this lesion analogue was placed in a primer-template system, extension of the primer by either avian myeloblastosis virus reverse transcriptase (AMV-RT) or human DNA polymerase β (pol β), was significantly impaired, but where incorporation opposite 8-nitroguanine did occur, pol β showed a 2:1 preference to insert dA over dC, while AMV-RT incorporated predominantly dC. The fact that no 8-nitroG·G base pairing is seen in the primer extension products suggests that the polymerases may discriminate against this pairing system on the basis of its poor geometric match to a Watson-Crick pair.

Acid dissociation constants of uridine-5'-diphosphate compounds determined by 31phosphorus nuclear magnetic resonance spectroscopy and internal pH referencing

The acid dissociation constant (pK(a)) of small, biological molecules is an important physical property used for investigating enzyme mechanisms and inhibitor design. For phosphorus-containing molecules, the (31)P nuclear magnetic resonance (NMR) chemical shift is sensitive to the local chemical environment, particularly to changes in the electronic state of the molecule. Taking advantage of this property, we present a (31)P NMR approach that uses inorganic phosphate buffer as an internal pH reference to determine the pK(a) values of the imide and second diphosphate of uridine-5'-diphosphate compounds, including the first reported values for UDP-GlcNAc and UDP-S-GlcNAc. New methods for using inorganic phosphate buffer as an internal pH reference, involving mathematical correction factors and careful control of the chemical shift reference sample, are illustrated. A comparison of the newly determined imide and diphosphate pK(a) values of UDP, UDP-GlcNAc, and UDP-S-GlcNAc with other nucleotide phosphate and thio-analogs reveals the significance of the monosaccharide and sulfur position on the pK(a) values.

Recovery Based Nanowire Field-Effect Transistor Detection of Pathogenic Avian Influenza DNA

Fast and accurate diagnosis is critical in infectious disease surveillance and management. We proposed a DNA recovery system that can easily be adapted to DNA chip or DNA biosensor for fast identification and confirmation of target DNA. This method was based on the re-hybridization of DNA target with a recovery DNA to free the DNA probe. Functionalized silicon nanowire field-effect transistor (SiNW FET) was demonstrated to monitor such specific DNA–DNA interaction using high pathogenic strain virus hemagglutinin 1 (H1) DNA of avian influenza (AI) as target. Specific electric changes were observed in real-time for AI virus DNA sensing and device recovery when nanowire surface of SiNW FET was modified with complementary captured DNA probe. The recovery based SiNW FET biosensor can be further developed for fast identification and further confirmation of a variety of influenza virus strains and other infectious diseases.

Nucleobase recognition at alkaline pH and apparent pKa of single DNA bases immobilised within a biological nanopore

The four DNA bases are recognized in immobilized DNA strands at high alkaline pH by nanopore current recordings. Ionic currents through the biological nanopores are also employed to measure the apparent pK(a) values of single nucleobases within the immobilised DNA strands.

Phase behavior of crowded like-charged mixed polyelectrolytes in a cell-sized sphere

We studied the phase behavior of a mixture of two semiflexible negatively charged polyelectrolytes, giant DNA and alginate, under crowded condition in a cell-sized sphere (5-40 μm in diameter), where the persistence length of giant DNA is 50 nm and that of alginate is 5 nm. Through microscopic observation, we found that the polymer mixture exhibits a unique phase behavior, which depends on the size of the sphere, whereas the mixture remains homogeneous and isotropic in bulk solution. When the sphere is small, DNA is completely depleted on the surface. When the sphere is medium sized, a portion of the DNA is depleted on the surface, and the remainder stays within the sphere. By introducing a curvature-dependent term for the interaction between DNA and the surface into the Flory-Huggins model, we interpret the observed characteristics of the phase behavior in terms of the relative importance of the surface-to-volume effect.

The solvation, partitioning, hydrogen bonding, and dimerization of nucleotide bases: a multifaceted challenge for quantum chemistry

We present M06-2X density functional calculations of the chloroform/water partition coefficients of cytosine, thymine, uracil, adenine, and guanine and calculations of the free energies of association of selected unsubstituted and alkylated nucleotide base pairs in chloroform and water. Both hydrogen bonding and π-π stacking interactions are considered. Solvation effects are treated using the continuum solvent models SM8, SM8AD, and SMD, including geometry optimization in solution. Comparison of theoretical results with available experimental data indicates that all three of these solvation models predict the chloroform-water partition coefficients for the studied nucleobases qualitatively well, with mean unsigned errors in the range of 0.4-1.3 log units. All three models correctly predict the preference for hydrogen bonding over stacking for nucleobase pairs solvated in chloroform, and SM8, SM8AD, and SMD show similar accuracy in predicting the corresponding free energies of association. The agreement between theory and experiment for the association free energies of the dimers in water is more difficult to assess, as the relevant experimental data are indirect. Theory predicts that the stacking interaction of nucleobases in water is more favorable than hydrogen bonding for only two out of three tested hetero-dimers.

Carba-LNA-5MeC/A/G/T modified oligos show nucleobase-specific modulation of 3'-exonuclease activity, thermodynamic stability, RNA selectivity, and RNase H elicitation: synthesis and biochemistry

Using the intramolecular 5-exo-5-hexenyl radical as a key cyclization step, we previously reported an unambiguous synthesis of carba-LNA thymine (cLNA-T), which we subsequently incorporated in antisense oligonucleotides (AON) and investigated their biochemical properties [J. Am. Chem. Soc.2007, 129 (26), 8362-8379]. These cLNA-T incorporated oligos showed specific RNA affinity of +3.5-5 °C/modification for AON:RNA heteroduplexes, which is comparable to what is found for those of LNAs (Locked Nucleic Acids). These modified oligos however showed significantly enhanced nuclease stability (ca. 100 times more) in the blood serum compared to those of the LNA modified counterparts without compromising any RNase H recruitment capability. We herein report the synthesis of 5-methylcytosine-1-yl ((Me)C), 9-adeninyl (A), and 9-guaninyl (G) derivatives of cLNA and their oligonucleotides and report their biochemical properties as potential RNA-directed inhibitors. In a series of isosequential carba-LNA modified AONs, we herein show that all the cLNA modified AONs are found to be RNA-selective, but the magnitude of RNA-selectivity of 7'-R-Me-cLNA-G (cLNA-G) (ΔT(m) = 2.9 °C/modification) and intractable isomeric mixtures of 7'-(S/R)-Me-cLNA-T (cLNA-T, ΔT(m) = 2.2 °C/modification) was found to be better than diastereomeric mixtures of 7'-(S/R)-Me-cLNA-(Me)C with trace of cENA-(Me)C (cLNA-(Me)C, ΔT(m) = 1.8 °C/modification) and 7'-R-Me-cLNA-A (cLNA-A, ΔT(m) = 0.9 °C/modification). cLNA-(Me)C modified AONs however exhibited the best nuclease stability, which is 4-, 7-, and 20-fold better, respectively, than cLNA-T, cLNA-A, and cLNA-G modified counterparts, which in turn was more than 100 times stable than that of the native. When the modification sites are appropriately chosen in the AONs, the cLNA-A, -G, and -(Me)C modified sites in the AON:RNA hybrids can be easily recognized by RNase H, and the RNA strand of the hybrid is degraded in a specific manner, which is important for the design of oligos for therapeutic purposes. The cLNA-(Me)C modified AON/RNA, however, has been found to be degraded 4 times faster than cLNA-A and G modified counterparts. By appropriately choosing the carba-LNA modification sites in AON strands, the digestion of AON:RNA can be either totally repressed or be limited to cleavage at specific sites or at a single site only (similar to that of catalytic RNAzyme or DNAzyme). Considering all physico- and biochemical aspects of cLNA modified oligos, the work suggests that the cLNA modified antisense oligos have the potential of being a promising therapeutic candidate due to their (i) higher nucleobase-specific RNA affinity and RNA selectivity, (ii) greatly improved nuclease stability, and (iii) efficient RNase H recruitment capability, which can induce target RNA cleavage in a very specific manner at multiple or at a single site, in a designed manner.

A convenient protection for 4-oxopyrimidine moieties in nucleosides by the pivaloyl group

Application of the pivaloyl group as a protection for the N3 position of thymidine and uridine was investigated. Pivaloylation of thymidine is a very rapid reaction proceeding under mild conditions with excellent regioselectivity for sugar or thymine moiety, depending on the amines used. Several pivaloylated thymidine derivatives were obtained by treatment of unprotected thymidine with pivaloyl chloride under various experimental conditions. Stability of the N3-pivaloyl protecting group under basic and acidic conditions was evaluated and the conditions for its selective removal were found.

DNA strands from denatured duplexes are translocated through engineered protein nanopores at alkaline pH

Nanopores are under development for the detection of a variety of analytes and the investigation of chemical reactions at the single molecule level. In particular, the analysis of nucleic acid molecules is under intense investigation, including the development of systems for rapid, low-cost DNA sequencing. Here, we show that DNA can be translocated through an engineered alphaHL protein pore at pH 11.7, a value at which dsDNA is denatured. Therefore, the alphaHL pore is sufficiently stable to entertain the possibility of direct nanopore sequencing of genomic dsDNA samples, which are more readily obtained and handled than ssDNA.

An analysis of the different behavior of DNA and RNA through the study of the mutual relationship between stacking and hydrogen bonding

The mutual relationship between stacking and hydrogen bonding and the possible influence of stacking in the different behavior of DNA and RNA base pairs have been studied through complete DFT optimization of different structures of A-U and A-T dimers (i.e., A-U/U-A and A-T/T-A), using some functionals developed by the group of Truhlar. The results obtained in this work clearly show that stacking and hydrogen bonding are deeply connected. The different behavior of DNA and RNA when replacing uracil by thymine can be interpreted through the formation of a stabilizing CH/pi interaction between the methyl group of thymine and the five-member ring of adenine.

DNA with branched internal side chains: synthesis of 5-tripropargylamine-dU and conjugation by an azide-alkyne double click reaction

5-Tripropargylamine-2'-deoxyuridine (1 a) containing two terminal triple bonds was synthesized by a Pd-assisted Sonogashira cross-coupling reaction and was subsequently converted into the corresponding phosphoramidite building block (9) and employed in solid-phase oligonucleotide synthesis. T(m) experiments demonstrate that the presence of covalently attached branched tripropargylamine residues has a positive effect on the base pair stability. The two terminal C[triple chemical bond]C bonds of modified DNA were functionalized by means of Cu(I)-mediated 1,3-dipolar cycloaddition reactions (click chemistry) with azides such as 3-azido-7-hydroxycoumarin or 3'-azido-3'-deoxythymidine (AZT) both in solution and on solid support. In particular, with the nonfluorescent 3-azido-7-hydroxycoumarin a strongly fluorescent oligonucleotide bis-dye conjugate was generated. For comparison, the N(3)-propargylated 2'-deoxyuridine 2 was prepared from 2'-deoxyuridine and propargyl bromide and incorporated into DNA. The two terminal triple bonds of 1 a allow the simultaneous post-modification of DNA by two reporter molecules and can be applied to almost any azido derivatives (oligonucleotides, proteins, polysaccharides etc.) including those forming dendrimeric side chains.

Comparison of the acid-base properties of ribose and 2'-deoxyribose nucleotides

The extent to which the replacement of a ribose unit by a 2'-deoxyribose unit influences the acid-base properties of nucleotides has not hitherto been determined in detail. In this study, by potentiometric pH titrations in aqueous solution, we have measured the acidity constants of the 5'-di- and 5'-triphosphates of 2'-deoxyguanosine [i.e., of H(2)(dGDP)(-) and H(2)(dGTP)(2-)] as well as of the 5'-mono-, 5'-di-, and 5'-triphosphates of 2'-deoxyadenosine [i.e., of H(2)(dAMP)(+/-), H(2)(dADP)(-), and H(2)(dATP)(2-)]. These 12 acidity constants (of the 56 that are listed) are compared with those of the corresponding ribose derivatives (published data) measured under the same experimental conditions. The results show that all protonation sites in the 2'-deoxynucleotides are more basic than those in their ribose counterparts. The influence of the 2'-OH group is dependent on the number of 5'-phosphate groups as well as on the nature of the purine nucleobase. The basicity of N7 in guanine nucleotides is most significantly enhanced (by about 0.2 pK units), while the effect on the phosphate groups and the N1H or N1H(+) sites is less pronounced but clearly present. In addition, (1)H NMR chemical shift change studies in dependence on pD in D(2)O have been carried out for the dAMP, dADP, and dATP systems, which confirmed the results from the potentiometric pH titrations and showed the nucleotides to be in their anti conformations. Overall, our results are not only of relevance for metal ion binding to nucleotides or nucleic acids, but also constitute an exact basis for the calculation, determination, and understanding of perturbed pK(a) values in DNAzymes and ribozymes, as needed for the delineation of acid-base mechanisms in catalysis.

Impact of chemical and structural anisotropy on the electrophoretic mobility of spherical soft multilayer particles: the case of bacteriophage MS2

We report a theoretical investigation of the electrohydrodynamic properties of spherical soft particles composed of permeable concentric layers that differ in thickness, soft material density, chemical composition, and flow penetration degree. Starting from a recent numerical scheme developed for the computation of the direct-current electrophoretic mobility (mu) of diffuse soft bioparticles, the dependence of mu on the electrolyte concentration and solution pH is evaluated taking the known three-layered structure of bacteriophage MS2 as a supporting model system (bulk RNA, RNA-protein bound layer, and coat protein). The electrokinetic results are discussed for various layer thicknesses, hydrodynamic flow penetration degrees, and chemical compositions, and are discussed on the basis of the equilibrium electrostatic potential and hydrodynamic flow field profiles that develop within and around the structured particle. This study allows for identifying the cases where the electrophoretic mobility is a function of the inner structural and chemical specificity of the particle and not only of its outer surface properties. Along these lines, we demonstrate the general inapplicability of the notions of zeta potential (zeta) and surface charge for quantitatively interpreting electrokinetic data collected for such systems. We further shed some light on the physical meaning of the isoelectric point. In particular, numerical and analytical simulations performed on structured soft layers in indifferent electrolytic solution demonstrate that the isoelectric point is a complex ionic strength-dependent signature of the flow permeation properties and of the chemical and structural details of the particle. Finally, the electrophoretic mobilities of the MS2 virus measured at various ionic strength levels and pH values are interpreted on the basis of the theoretical formalism aforementioned. It is shown that the electrokinetic features of MS2 are to a large extent determined not only by the external proteic capsid but also by the chemical composition and hydrodynamic flow permeation of/within the inner RNA-protein bound layer and bulk RNA part of the bacteriophage. The impact of virus aggregation, as revealed by decreasing diffusion coefficients for decreasing pH values, is also discussed.

Self-assembly of nucleoamphiphiles: investigating nucleosides effect and the mechanism of micrometric helix formation

A new family of self-assembling systems based on nucleoamphiphiles is described. Nano to micrometric left-handed helix formation in aqueous solution was induced simply by complexing a GMP or an AMP with a nonchiral monocationic amphiphile. The assembling behavior such as micellar formation, monolayer at air-water interface, as well as the aggregates in solution of these nucleoamphiphiles are strongly influenced by the presence of nucleosides in solution. The observed effects depend on the properties of complexed nucleotides and nucleosides with a complex mixture of pi stacking, hydrophobicity of the bases, and hydrogen bonding.

Modification of DNA with octadiynyl side chains: synthesis, base pairing, and formation of fluorescent coumarin dye conjugates of four nucleobases by the alkyne--azide "click" reaction

Oligonucleotides incorporating 5-(octa-1,7-diynyl)-2'-deoxycytidine 1a, 5-(octa-1,7-diynyl)-2'-deoxyuridine 2a and 7-deaza-7-(octa-1,7-diynyl)-2'-deoxyguanosine 3a, 7-deaza-7-(octa-1,7-diynyl)-2'-deoxyadenosine 4a were prepared. For this, the phosphoramidites 7, 10, and 13 were synthesized and employed in solid-phase oligonucleotide synthesis. The octa-1,7-diynyl nucleosides 1a- 4a were obtained from their corresponding iodo derivatives using the palladium-assisted Sonogashira cross-coupling reaction. The Tm values demonstrated that DNA duplexes containing octa-1,7-diynyl nucleosides show a positive influence on the DNA duplex stability when they are introduced at the 5-position of pyrimidines or at the 7-position of 7-deazapurines. The terminal alkyne residue of oligonucleotides were selectively conjugated to the azide residue of the nonfluorescent 3-azido-7-hydroxycoumarin ( 38) using the protocol of copper(I)-catalyzed [3 + 2] Huisgen--Sharpless--Meldal cycloaddition "click chemistry" resulting in the formation of strongly fluorescent 1,2,3-triazole conjugates. The fluorescence properties of oligonucleotides with covalently linked coumarin--nucleobase assemblies were investigated. Among the four modified bases, the 7-deazapurines show stronger fluorescence quenching than that of pyrimidines.

Acid-base and metal-ion binding properties of the RNA dinucleotide uridylyl-(5'-->3')-[5']uridylate (pUpU3-)

It is well known that Mg2+ and other divalent metal ions bind to the phosphate groups of nucleic acids. Subtle differences in the coordination properties of these metal ions to RNA, especially to ribozymes, determine whether they either promote or inhibit catalytic activity. The ability of metal ions to coordinate simultaneously with two neighboring phosphate groups is important for ribozyme structure and activity. However, such an interaction has not yet been quantified. Here, we have performed potentiometric pH titrations to determine the acidity constants of the protonated dinucleotide H2(pUpU)-, as well as the binding properties of pUpU3- towards Mg2+, Mn2+, Cd2+, Zn2+, and Pb2+. Whereas Mg2+, Mn2+, and Cd2+ only bind to the more basic 5'-terminal phosphate group, Pb2+, and to a certain extent also Zn2+, show a remarkably enhanced stability of the [M(pUpU)]- complex. This can be attributed to the formation of a macrochelate by bridging the two phosphate groups within this dinucleotide by these metal ions. Such a macrochelate is also possible in an oligonucleotide, because the basic structural units are the same, despite the difference in charge. The formation degrees of the macrochelated species of [Zn(pUpU)]- and [Pb(pUpU)]- amount to around 25 and 90 %, respectively. These findings are important in the context of ribozyme and DNAzyme catalysis, and explain, for example, why the leadzyme could be selected in the first place, and why this artificial ribozyme is inhibited by other divalent metal ions, such as Mg2+.

Influence of dT20 and [d(AT)10]2 on cisplatin hydrolysis studied by two-dimensional [1H,15N] HMQC NMR spectroscopy

The influence of the presence of DNA on the kinetics of cisplatin (cis-[PtCl2(NH3)2]) aquation (replacement of Cl- by H2O) and anation (replacement of H2O by Cl-) involved in the hydrolysis of cisplatin have been determined by two-dimensional [1H,15N] HMQC NMR spectroscopy. Single-stranded dT20 and double-stranded [d(AT)10]2 oligonucleotides were used as DNA models, avoiding guanines which are known to react rapidly with aquated cisplatin forms. Reactions starting from cis-[PtCl2(15NH3)2], or from a stoichiometric mixture of cis-[Pt(15NH3)2(H2O)2]2+ and Cl- (all 0.5 mM Pt(II); in ionic strength, adjusted to 0.095 M or 0.011 M with NaClO4, pH between 3.0 and 4.0) were followed in an NMR tube in both the absence and presence of 0.7 mM dT20 or [d(AT)10]2. In the presence of dT20, we observed a slight and ionic-strength-independent decrease (15-20 %) of the first aquation rate constant, and a more significant decrease of the second anation rate constant. The latter was more important at low ionic strength, and can be explained by efficient condensation of cis-[Pt(15NH3)2(H2O)2]2+ on the surface of single-stranded DNA, in a region depleted of chloride anions. At low ionic strength, we observed an additional set of [1H,15N] HMQC spectral signals indicative of an asymmetric species of PtN2O2 coordination, and we assigned them to phosphate-bound monoadducts of cis-[Pt(15NH3)2(H2O)2]2+. Double-stranded [d(AT)10]2 slowed down the first aquation step also by approximately 15 %; however, we could not determine the influence on the second hydrolysis step because of a significant background reaction with cis-[Pt(NH3)2(H2O)2]2+.

Computational and empirical trans-hydrogen bond deuterium isotope shifts suggest that N1-N3 A:U hydrogen bonds of RNA are shorter than those of A:T hydrogen bonds of DNA

Density functional theory calculations of isolated Watson-Crick A:U and A:T base pairs predict that adenine 13C2 trans-hydrogen bond deuterium isotope shifts due to isotopic substitution at the pyrimidine H3, (2h)Delta13C2, are sensitive to the hydrogen-bond distance between the N1 of adenine and the N3 of uracil or thymine, which supports the notion that (2h)Delta13C2 is sensitive to hydrogen-bond strength. Calculated (2h)Delta13C2 values at a given N1-N3 distance are the same for isolated A:U and A:T base pairs. Replacing uridine residues in RNA with 5-methyl uridine and substituting deoxythymidines in DNA with deoxyuridines do not statistically shift empirical (2h)Delta13C2 values. Thus, we show experimentally and computationally that the C7 methyl group of thymine has no measurable affect on (2h)Delta13C2 values. Furthermore, (2h)Delta13C2 values of modified and unmodified RNA are more negative than those of modified and unmodified DNA, which supports our hypothesis that RNA hydrogen bonds are stronger than those of DNA. It is also shown here that (2h)Delta13C2 is context dependent and that this dependence is similar for RNA and DNA.

The chemical nature of the 2'-substituent in the pentose-sugar dictates the pseudoaromatic character of the nucleobase (pKa) in DNA/RNA

We here show that the pKa (error limit: 0.01 to 0.03 pKa unit) of a nucleobase in a nucleotide can be modulated by the chemical nature of the 2'-substituent at the sugar moiety. This has been evidenced by the measurement of nucleobase pKa in 47 different model nucleoside 3',5'-bis- and 3'-mono-ethylphosphates. The fact that the electronic character of each of the 2'-substituents (Fig. 1) alters the chemical shift of the H2' sugar proton, and also alters the pKa of the nucleobase in the nucleotides has been evidenced by a correlation plot of pKa of N3 of pyrimidine (T/C/U) or pKa of N7 of 9-guaninyl with the corresponding deltaH2' chemical shifts at the neutral pH, which shows linear correlation with high Pearson's correlation coefficients (R = 0.85-0.97). That this modulation of the pKa of the nucleobase by a 2'-substituent is a through-bond as well as through-space effect has been proven by ab initio determined pKa estimation. Interestingly, experimental pKas of nucleobases from NMR titration and the calculated pKas (by ab initio calculations utilizing closed shell HF 6-31G** basis set) are linearly correlated with R = 0.98. It has also been observed that the difference of ground and protonated/de-protonated HOMO orbital energies (DeltaHOMO, a.u.) for the nucleobases (A/G/C/T/U) are well correlated with their pK(a)s in different 2'-substituted 3',5'-bis-ethylphosphate analogs suggesting that only the orbital energy of HOMO can be successfully used to predict the modulation of the chemical reactivity of the nucleobase by the 2'-substituent. It has also been demonstrated that pKa values of nucleobases in 3',5'-bis-ethylphosphates (Table 1) are well correlated with the change in dipole moment for the respective nucleobases after protonation or de-protonation. This work thus unambiguously shows that alteration of the thermodynamic stability (Tm) of the donor-acceptor complexes [ref. 20], as found with various 2'-modified duplexes in the antisense, siRNA or in triplexes by many workers in the field, is a result of alteration of the pseudoaromatic character of the nucleobases engineered by alteration of the chemical nature of the 2'-substitution.