Cross-modulation of physicochemical character of aglycones in dinucleoside (3'-->5') monophosphates by the nearest neighbor interaction in the stacked state

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.

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.

Determination of pKa values for deprotonable nucleobases in short model oligonucleotides

The deprotonation of ionizable nucleobases centrally placed in short model oligonucleotides was examined under different physical conditions, using UV absorption spectroscopy. The oligonucleotide sequences were designed so that only the central base would be ionized over the pH range examined. pKa values of 9.90±0.01 and 9.34±0.04 were determined for the guanine group in the oligomer d-ACAGCAC and 2'-deoxyguanosine, respectively, both at 25°C and 0.1M NaCl. Lengthening the oligonucleotide up to the tridecamer stage further increases the pKa of the central guanine moiety. Electrolyte concentration, temperature, and mixed water-ethanol solvents affect the acidity of the central base. Changes in the sequence surrounding the central guanine can also have a significant effect, especially in the case of strongly stacking sequences. The pKa values were also determined for the hepta(2'-O-methyl)ribonucleotide and the heptamer PNA of identical sequence, as well as for oligodeoxyribonucleotides with different deprotonable bases, viz. thymine, uracil, or hypoxanthine, in the central position. The results are interpreted in terms of the electric-field effect exerted on the departing proton by the negative electric charges located on the internucleotide phosphate groups, and calculations show this effect to approximately explain the magnitude of the pKa difference observed between the deoxyriboheptanucleotide and its electroneutral PNA analogue.

Enzymatic polymerization of phosphonate nucleosides

5'-O-phosphonomethyl-2'-deoxyadenosine (PMdA) proved to be a good substrate of the Therminator polymerase. In this article, we investigated whether the A, C, T and U analogues of this phosphonate nucleoside (PMdN) series can function as substrates of natural DNA polymerases. PMdT and PMdU could only be polymerized enzymatically to a limited extent. Nevertheless, PMdA and PMdC could be incorporated into a DNA duplex with complete chain elongation by all the DNA polymerases tested. A mixed sequence of four nucleotides containing modified C, T and A residues could be obtained with the Vent(exo(-)) and Therminator polymerases. The kinetic values for the incorporation of PMdA by Vent(exo(-)) polymerase were determined; a reduced K(M) value was found for the incorporation of PMdA compared to the natural substrate. Future polymerase directed evolution studies will allow us to select an enzyme with a heightened capacity to process these modified DNA building blocks into modified strands.

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.

Non-identical electronic characters of the internucleotidic phosphates in RNA modulate the chemical reactivity of the phosphodiester bonds

We here show that the electronic properties and the chemical reactivities of the internucleotidic phosphates in the heptameric ssRNAs are dissimilar in a sequence-specific manner because of their non-identical microenvironments, in contrast with the corresponding isosequential ssDNAs. This has been evidenced by monitoring the delta H8(G) shifts upon pH-dependent ionization (pK(a1)) of the central 9-guaninyl (G) to the 9-guanylate ion (G-), and its electrostatic effect on each of the internucleotidic phosphate anions, as measured from the resultant delta 31P shifts (pKa2) in the isosequential heptameric ssRNAs vis-à-vis ssDNAs: [d/r(5'-Cp1Ap2Q1p3Gp4Q2p5Ap6C-3'): Q1 = Q2 = A (5a/5b) or C (8a/8b), Q(1) = A, Q(2) = C (6a/6b), Q1 = C, Q2 = A (7a/7b)]. These oligos with single ionizable G in the centre are chosen because of the fact that the pseudoaromatic character of G can be easily modulated in a pH-dependent manner by its transformation to G- (the 2'-OH to 2-O- ionization effect is not detectable below pH 11.6 as evident from the N(1-Me)-G analog), thereby modulating/titrating the nature of the electrostatic interactions of G to G- with the phosphates, which therefore constitute simple models to interrogate how the variable pseudoaromatic characters of nucleobases under different sequence context (J. Am. Chem. Soc., 2004, 126, 8674-8681) can actually influence the reactivity of the internucleotide phosphates as a result of modulation of sequence context-specific electrostatic interactions. In order to better understand the impact of the electrostatic effect of the G to G- on the tunability of the electronic character of internucleotidic phosphates in the heptameric ssRNAs 5b, 6b, 7b and 8b, we have also performed their alkaline hydrolysis at pH 12.5 at 20 degrees C, and have identified the preferences of the cleavage sites at various phosphates, which are p2, p3 and p4 (Fig.3). The results of these alkaline hydrolysis studies have been compared with the hydrolysis of analogous N(1-Me)-G heptameric ssRNA sequences 5c, 7c and 8c under identical conditions in order to establish the role of the electrostatic effect of the 9-guanylate ion (and the 2'-OH to 2-O- ionization) on the internucleotidic phosphate. It turned out that the relative alkaline hydrolysis rate at those particular phosphates (p2, p3 and p4) in the N(1-Me)-G heptamers was reduced from 16-78% compared to those in the native counterparts [Fig. 4, and ESI 2 (Fig. S11)]. Thus, these physico-chemical studies have shown that those p2, p3 and p4 phosphates in the native heptameric RNAs, which show pKa2 as well as more deshielding (owing to weaker 31P screening) in the alkaline pH compared to those at the neutral pH, are more prone to the alkaline hydrolysis because of their relatively enhanced electrophilic character resulting from weaker 31P screening. This screening effect originates as a result of the systematic charge repulsion effect between the electron cloud in the outermost orbitals of phosphorus and the central guanylate ion, leading to delocalization of the phosphorus p(pi) charge into its dpi orbitals. It is thus likely that, just as in the non-enzymatic hydrolysis, the enzymatic hydrolysis of a specific phosphate in RNA by general base-catalysis in RNA-cleaving proteins (RNase A, RNA phosphodiesterase or nuclease) can potentially be electrostatically influenced by tuning the transient charge on the nucleobase in the steric proximity or as a result of specific sequence context owing to nearest-neighbor interactions.

Single-stranded adenine-rich DNA and RNA retain structural characteristics of their respective double-stranded conformations and show directional differences in stacking pattern

The structural preorganization of isosequential ssDNA and ssRNA hexamers d/r(GAAAAC)(1) [J. Am. Chem. Soc. 2003, 125, 9948] have been investigated by NMR and molecular dynamics simulations. Analysis of the nuclear Overhauser effect spectrometry (NOESY) footprints in the aqueous solution has shown that there is a substantial population of ordered right-handed helical structure in both hexameric single-stranded DNA and RNA, which are reminiscent of their respective right-handed helical duplex form, despite the fact these single-stranded molecules are devoid of any intermolecular hydrogen bonds. The NMR-constrained molecular dynamics (1.5 ns) derived geometries of the adenine-adenine overlaps at each dinucleotide step of the hexameric ssDNA (1a) and ssRNA (1b) show that the relatively electron-rich imidazole stacks above the electron-deficient pyrimidine in 5' to 3' direction in ssDNA (1a) while, in contradistinction, the pyrimidine stacks above the imidazole in the 5' to 3' direction in ssRNA (1b). This also means that the pi-frame of the 5'-pyrimidine can interact with the relatively positively charged imino and amino protons in the 3' direction in ssRNA and in the 5' direction in ssDNA, thereby stabilizing the twist and slide observed in the stacked oligonucleotides. The differently preferred stacking geometries in ssDNA and ssRNA have direct physicochemical implications for self-assembly and pK(a) modulation by the nearest-neighbor interactions, as well as for the dangling-end stabilization effects and imino-proton reactivity.

Significant pKa perturbation of nucleobases is an intrinsic property of the sequence context in DNA and RNA

The pH titration and NMR studies (pH 6.6-12.5) in the heptameric isosequential ssDNA and ssRNA molecules, [d/r(5'-CAQ1GQ2AC-3', with variable Q1/Q2)], show that the pKa of the central G residue within the heptameric ssDNAs (DeltapKa = 0.67 +/- 0.03) and ssRNAs (DeltapKa = 0.49 +/- 0.02) is sequence-dependent. This variable pKa of the G clearly shows that its pseudoaromatic character, hence, its chemical reactivity, is strongly modulated and tuned by its sequence context. In contradistinction to the ssDNAs, the electrostatic transmission of the pKa of the G moiety to the neighboring A or C residues in the heptameric ssRNAs (as observed by the response of the aromatic marker protons of As or Cs) is found to be uniquely dependent upon the sequence composition. This demonstrates that the neighboring As or Cs in ssRNAs have variable electrostatic efficiency to interact with the central G/G-, which is owing to the variable pseudoaromatic characters (giving variable chemical reactivities) of the flanking As or Cs compared to those of the isosequential ssDNAs. The sequence-dependent variation of pKa of the central G and the modulation of its pKa transmission through the nearest-neighbors by variable electrostatic interaction is owing to the electronically coupled nature of the constituent nucleobases across the single strand, which demonstrates the unique chemical basis of the sequence context specificity of DNA or RNA in dictating the biological interaction, recognition, and function with any specific ligand.

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

To understand why the RNA-RNA duplexes in general has a higher thermodynamic stability over the corresponding DNA-DNA duplexes, we have measured the pK(a) values of both nucleoside 3',5'-bis-ethyl phosphates [Etp(d/rN)pEt] and nucleoside 3'-ethyl phosphates [(d/rN)pEt] (N = A, G, C, or T/U), modeling as donors and acceptors of base pairs in duplexes. While the 3',5'-bis-phosphates, Etp(d/rN)pEt, mimic the internucleotidic monomeric units of DNA and RNA, in which the stacking contribution is completely absent, the 3'-ethyl phosphates, (d/rN)pEt, mimic the nucleotide at the 5'-end. The pK(a) values of the nucleobase in each of these model nucleoside phosphates have been determined with low pK(a) error (sigma = +/-0.01 to 0.02) by (1)H NMR (at 500 MHz) with 20-33 different pH measurements for each compound. This study has led us to show the following: (1) All monomeric DNA nucleobases are more basic than the corresponding RNA nucleobases in their respective Etp(d/rN)pEt and (d/rN)pEt. (2) The pK(a) values of the monomeric nucleotide blocks as well as Delta pK(a) values between the donor and acceptor can be used to understand the relative base-pairing strength in the oligomeric duplexes in the RNA and DNA series. (3) The Delta G*(pKa) of the donor and acceptor of the base pair in duplexes enables a qualitative dissection of the relative strength of the base-pairing and stacking in the RNA-RNA over the DNA-DNA duplexes. (4) It is also found that the relative contribution of base-pairing strength and nucleobase stacking in RNA-RNA over DNA-DNA is mutually compensating as the % A-T/U content increases or decreases. This interdependency of stacking and hydrogen bonding can be potentially important in the molecular design of the base-pair mimics to expand the alphabet of the genetic code.

Cross-modulation of the pKa of nucleobases in a single-stranded hexameric-RNA due to tandem electrostatic nearest-neighbor interactions

The pH titration studies (pH 6.7-12.1) in a series of dimeric, trimeric, tetrameric, pentameric, and hexameric oligo-RNA molecules [GpA (2a), GpC (3a), GpApC (5), GpA(1)pA(2)pC (6), GpA(1)pA(2)pA(3)pC (7), GpA(1)pA(2)pA(3)pA(4)pC (8)] have shown that the pK(a) of N(1)-H of 9-guaninyl could be measured not only from its own deltaH8G, but also from the aromatic marker protons of other constituent nucleobases. The relative chemical shift differences [Deltadelta((N)(-)(D))] between the protons in various nucleotide residues in the oligo-RNAs at the neutral (N) and deprotonated (D) states of the guanine moiety show that the generation of the 5'-(9-guanylate ion) in oligo-RNAs 2-8 reduces the stability of the stacked helical RNA conformation owing to the destabilizing anion(G(-))-pi/dipole(Im(delta)(-)) interaction. This destabilizing effect in the deprotonated RNA is, however, opposed by the electrostatically attractive atom-pisigma (major) as well as the anion(G(-))-pi/dipole(Py(delta)(+)) (minor) interactions. Our studies have demonstrated that the electrostatically repulsive anion(G(-))-pi/dipole(Im(delta)(-)) interaction propagates from the first to the third nucleobase quite strongly in the oligo-RNAs 6-8, causing destacking of the helix, and then its effect is gradually reduced, although it is clearly NMR detectable along the RNA chain. Thus, such specific generation of a charge at a single nucleobase moiety allows us to explore the relative strength of stacking within a single-stranded helix. The pK(a) of 5'-Gp residue from its own deltaH8G in the hexameric RNA 8 is found to be 9.76 +/- 0.01; it, however, varies from 9.65 +/- 0.01 to 10.5 +/- 0.07 along the RNA chain as measured from the other marker protons (H2, H8, H5, and H6) of 9-adeninyl and 1-cytosinyl residues. This nucleobase-dependent modulation of pK(a)s (DeltapK(a) +/- 0.9) of 9-guaninyl obtained from other nucleobases in the hexameric RNA 8 represents a difference of ca. 5.1 kJ mol(-)(1), which has been attributed to the variable strength of electrostatic interactions between the electron densities of the involved atoms in the offset stacked nucleobases as well as with that of the phosphates. The chemical implication of this variable pK(a) for guanin-9-yl deprotonation as obtained from all other marker protons of each nucleotide residue within a ssRNA molecule is that it enables us to experimentally understand the variation of the electronic microenvironment around each constituent nucleobase along the RNA chain in a stepwise manner with very high accuracy without having to make any assumption. This means that the pseudoaromaticity of neighboring 9-adeninyl and next-neighbor nucleobases within a polyanionic sugar-phosphate backbone of a ssRNA can vary from one case to another due to cross-modulation of an electronically coupled pi system by a neighboring nucleobase. This modulation may depend on the sequence context, spatial proximity of the negatively charged phosphates, as well as whether the offset stacking is ON or OFF. The net outcome of this electrostatic interaction between the neighbors is creation of new sequence-dependent hybrid nucleobases in an oligo- or polynucleotide whose properties are unlike the monomeric counterpart, which may have considerable biological implications.

Tandem electrostatic effect from the first to the third aglycon in the trimeric RNA owing to the nearest-neighbor interaction

We here show an electrostatic polar-pi interaction from the first to the third aglycon, via the second aglycon, in the ground state in two single stranded trimeric RNAs, 5'-GpA(1)pA(2)-3' (3) and 5'-GpApC-3' (4), as a result of intramolecular nearest neighbor offset-stacking. The experimental evidence in support of this conclusion has been obtained by comparing the pK(a)s of each aglycone in the two trimers with those of guanosine 3'-ethyl phosphate, GpEt (1) and 5'-GpA-3' (2): Thus, the pK(a) of N(1)-H of guanin-9-yl of 5'-GpA(1)pA(2)-3' (3) could be measured by pH titration (pH 7.3-11.6) of its own deltaH8G (pK(a) 9.75 +/- 0.02) as well as from deltaH8A(1) (pK(a) 9.72 +/- 0.02) and deltaH2A(1) (pK(a) 9.83 +/- 0.04) of the neighboring pA(1)p moiety and the deltaH8A(2) (pK(a) 9.83 +/- 0.02) of the terminal pA(2) moiety. Similarly, the pH titration of GpApC (4) shows the pK(a) of N(1)-H of guanin-9-yl from its own deltaH8G (pK(a) 9.88 +/- 0.03) as well as from deltaH8A (pK(a) 9.87 +/- 0.01) of the neighboring pAp moiety, and deltaH5/H6C (pK(a) 9.88 +/- 0.01 and 9.90 +/- 0.01 respectively) of the 3'-terminal cytosin-1-yl. This intramolecular nearest neighbor electrostatic interaction in the single-stranded RNA modulates the pseudoaromaticity of the nearest neighbors by almost total transmission of because they constitute an extended array of offset-stacked coupled aromatic heterocycles within a polyanionic sugar-phosphate backbone at the ground state. The enhanced basicity of Gp residue by ca. 0.6 pK(a) unit in the trimers compared to that of the dimer is a result of the change in the electrostatic microenvironment owing to the nearest neighbors in the former (the nucleobases as well as the phosphates). Thus, the from the 5'-guanylate ion to the 3'-end aglycon via the central adenin-9-yl is 55 to 56 kJ mol(-)(1) in each step through a distance spanning approximately 6.8 A in an unfolded state. As a result, the pK(a) of guanin-9-yl moiety has become 9.25 +/- 0.02 in GpEt (1), 9.17 +/- 0.02 in GpA (2), 9.75 +/- 0.02 in GpApA (3), and 9.88 +/- 0.03 in GpApC (4). This means that guanin-9-yl moiety of trimers 3 and 4 is more basic than in the monomer or the dimer. The net outcome of this electrostatic cross-talk between the two neighboring heterocycles is creation of new hybrid aglycones in an oligo or polynucleotide, whose physicochemical property and the pseudoaromatic character are completely dependent both upon the nearest neighbors, and whether they are stacked or unstacked. Thus, this tunable physicochemical property of an aglycon (an array of the extended genetic code) may have considerable implication in our understanding of the specific ligand binding ability of an aptamer, the pK(a) and the hydrogen bonding ability of nucleic acids in a microenvironment, or in the triplet usage by the anticodon-codon interaction in the protein biosynthesis in the ribosome.