Electrostatics in the Self-assembly of DNA and RNA
Turner and his coworkers were the first to realize the role of the nearest neighbor interactions in the self-organization and stabilization of the DNA and RNa duplexes. They also elucidated the free-energy contributions in the homo and hetero duplex stability owing to the the nearest neighbor interactions. Readers are cited to the Turner Group's home page for their pioneering work.(A) Electrostatic cross-modulation of the pseudoaromatic character in single-stranded RNA by nearest-neighbor interactions
We have however recently given concrete experimental evidences through pKa measurements (PDF 338, PDF 339, PDF 341, PDF 343, PDF 355, PDF 361) that the electrostatic interactions through the nearesst-neighbors play an important role in steering both intra and intermolecular forces culminating into the self-assembly of both single-stranded and homo and hetero duplexes of DNA and RNA. We indeed could show that because of cross-modulation of the aromatic character (evidenced by a change of pKa ) amongst the nearest-neighbors as a result of electrostatic cross-talk, a clear pKa modulation take place. The generality of this pKa modulation by electrostatic cross-talk was also clearly shown in a much simpler coupled system consisting of pyridine-benzene, covalently linked by linker, which was transformed to pyridinium-benzene by simple change of the pH. The change of the aromatic character from pyridine to pyridinium was electrostatically registered by the phenyl group in the steric proximity by showing the pKa modulation [PDF 343]
Basis of the above pKa modulation: Such electrostatic interactions is also expected to steer recognition and interactions of the ligand by the nucleic acid by the complementary electric field. The overall stability of nucleic acids, as well as their complexes with a ligand, is considered a sum of contributions of intermolecular electrostatic interactions, hydrogen bonding, van der Waals and hydrophobic interactions. ( PDF 338, PDF 339, PDF 341, PDF 343, PDF 355, PDF 361) Measurable entity of the different electrostatic/stereoelectronic effects in nucleic acids and protein folding is pKa perturbation. Determination and interpretation of pKa of an ionizable group (a basic center for protonation or a dissociable group) in biomolecules highlights the molecular mechanism of its biological function. The shift in pKa values is an important source of information about neighboring charges, electrostatics, structural perturbation as well as partial charge distribution over the whole molecule, and differential hydration of the microenvironment. Structural studies have revealed that in a large number of RNAs (and in certain cases in DNAs), the pKavalues of nucleobases (particularly adenine and cytosine) are significantly perturbed relative to that of mononucleotides. Thus, the most effective acid-base catalysis can be performed with ionizable groups having pKa near physiological pH of 7, thereby accounting for widespread use of histidines in protein enzyme active sites.
For reviews see Chapter 6 of Dr. Parag Acharya's Ph.D. Thesis. and Chapter 6 of Dr. Sandipta Acharya's Ph.D. thesis
The generation of a single anionic or cationic center at an alkaline or acidic pH in a given molecule presents a unique opportunity to examine the electrostatic make-up of these molecules both at the neutral or ionic state. The generation of a single cationic center in the phenyl-nicotinamide system provided new straightforward evidence showing that the charge density of the electron-deficient pyridinium was actually enhanced by the donation of the charge from the electron-rich phenyl group (i.e., the pyridinyl became more basic by ca. 0.5 pKa unit compared to an analogous system where phenyl was absent) owing to the electrostatic interactions between these two moieties. On the other hand, the generation of the 5'-guanylate ion in, for example, the hexameric single-strand (ss) RNA [5'-GAAAAC-3'], in comparison with the constituent trimeric [PDF 341], tetrameric, pentameric, and hexameric [PDF 344, 355]-ssRNAs, has unequivocally shown how far the electrostatic cross-talk (as an interplay of Coulombic attractive or repulsive forces) in this electronically coupled system propagates through the intervening pAp nucleotide steps until the terminal pC-3' residue in comparison with the neutral counterpart. The footprint of the propagation of this electrostatic cross-talk among the neighboring nucleobases is evident by measurement of pKas from the marker protons of ionization point (i.e., of G) as well as from the neighboring marker protons (i.e., of A or C) in the vicinity, as well as from the change of the chemical environment (i.e., chemical shifts) around their aromatic marker protons owing to a change of the stacking-destacking equilibrium as a function of pH.
Details are provided in References
as well as in the Ph.D. Thesises of Dr.
Parag Acharya and Dr.
B) Propagation of electrostatic interplay across the single strand (Papers 344 and 347)
Perturbed pKa is not repulsive destacking. What I meant in my papers is that the influence of guanylate ion is most felt on the stacked nucleobases by electrostatics because they are stacked.and in proximity..... with what is written here..... "It has been observed that the propagation of the interplay of various electrostatic forces across the ssDNA chain is considerably weak, since the perturbed pKa is only NMR detectable effectively up to the fourth nucleobase moiety in the hexameric ssDNA (15c). In comparison, this perturbed pKa is observed, by NMR titration experiment, for the corresponding isosequential ssRNA up to the sixth nucleobase residue, which measures upto ~21 Å in the unfolded state. This suggests that the intramoleclar attractive stacking is weaker in the isosequential ssDNA than those of the ssRNA. Here, we are focussing our attention to the relative penaly of stacking in ssDNA versus ssRNA. That is to compare the N-type sugar/ with pseudoaxial nucleobase in ssRNA versus S-type sugar with pseudoequatorial nucleobsae in ssDNA, in one hand, visavis more hydrophobic environment in ssDNA versus hydrophilic 2'-OH group in ssRNA compared to ssDNA. In this competition, if we assume that generation of a gradient of negative charge in guanylate ion will electrostatically influence by an influx of negative charges to those nuclebases which are most stacked and meaning they are in the orbital-overlapping proximity. Then, the stronger stacking in ssRNA, as evident by its most perturbed pKas of the constituent nucleobases, will clearly mean that the N-type sugar with pseudoaxial nucleobase in ssRNA is more stabilizing than the constituent hydrophilic 2'-OH group in each pentofuranose sugar [PDF 355]
Figure 1. Electrostatic Interplay across GpA1pA2pA3pA4pC
Figure 2. pKa perturbations
from nearest-neighbor interactions
(C) Base pairing versus stacking contribution in overall RNA-RNA and DNA-DNA duplex stabilization - a qualitative understanding from pKa measurement of model monomeric nucleotides (Paper 347 )
Figure 3. Watson-Crick Base Pairing
Because we can successfully measure the pKa of guanine-9-yl from either
of the aglycones in the RNA timers 3 and 4,it shows that the aglycones
in the trimeric RNAs constitute a coupled heterocyclic system right across
the pH range, 6.7 to 11.5 owing to both 3'->5' and 5'->3' two-way cross-modulation
by electrostatic interaction. This may be the reason a trimeric RNA sequence
constitutes a single codon signal in recognition and function in the protein
synthesis machinery. The magnitude of the chemical shift change in any
of the aromatic protons in either of the two coupled aglycones differs
in a variable manner depending upon the geometry of stacking, partial charge
of the heteroatom as well as the sequence
(compare GpApA and GpApC), which is evident from relative chemical shift change.
(1) Increasing content of A-T base pairs weakens the stability of DNA-DNA
duplex over the corresponding RNA-RNA duplexes
(2) Strength of stacking of A-T rich DNA-DNA sequence [due to CH3 (1-Thyminyl) -p (of 9-adeninyl) interaction] increases compared to A-U rich sequence in RNA-RNA duplexes
Conclusions (Paper 347 ):
(1) The nucleobases of the monomeric DNA are uniformly more basic than the corresponding RNA counterparts.
(2) The strength of the base-pairing based on the pKa difference (DpKa) of the monomeric donor and acceptor can be used to understand the relative base-pairing strength of larger oligomeric DNA-DNA and RNA-RNA duplexes.
(3) The use of DpKa among
the monomer blocks, modeling the A-T/U and G-C base pairs, allows us to
understand the base-pairing contributions in the free energy of DNA-DNA
and RNA-RNA duplex stability, which is evident from a high correlation
coefficient found between
DG°37 of the helix stability and the sum of the DpKa values of donor/acceptor in the base-pair formation.
(4) The high correlation of DG°37
of the helix stability and the sum of DpKa
in a duplex showed that a simple subtraction of the base-pairing contribution
of RNA-RNA over DNA-DNA ([DG°bp
]RR-DD) from the free energy of the total helix stability
gives us the relative stacking contribution ([DG°stacking
]RR-DD) in a qualitative manner.
Figure 4. Free energy stabilization of the base-pairing ([DG°bp ]RR-DD) [red triangles, R = 0.84] as well as that for stacking ([[DG°bp ]RR-DD) [blue squares, R = 0.85] in RNA-RNA (RR) over the DNA-DNA (DD) duplex as a function of %A-T/U bp content in the duplexes. (R - correlation coefficient).
(5) A comparison of these two linear plots ([DG°bp
]RR-DD) versus ([DG°stacking
]RR-DD) as a function of % A-T/U bp content, Figure
with opposite slope shows that with the increasing content of A-T base pairs the stability of DNA-DNA duplex weakens over the corresponding RNA-RNA duplexes ([DG°bp ]RR-DD) while the strength of stacking ([DG°stacking ]RR-DD) of A-T rich DNA-DNA sequence increases in comparison with A-U rich sequence in RNA-RNA duplexes. This increased stacking contribution from T compared to U, in DNA-DNA over RNA-RNA duplex, comes from favorable electrostatic CH/s interaction between the 5-methyl group of T with the nearest-neighbor A in the AT rich sequence.
See References [PDF
electronic characters of the internucleotidic phosphates in RNA modulate
the chemical reactivity of the phosphodiester bonds (Paper
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. [ PDF 361]
(E) Extended Genetic Code (Paper 341 )
The net result of this electrostatic cross-talk between two neighboring
aglycones as a result of base-base stacking is creation of a unique set
of aglycones in an oligo or polynucleotide, whose physico-chemical properties
are completely dependent upon the nearest neighbor electrostatic interactions.
This has considerable implication in the specific ligand binding ability,
aptamer recognition, RNA catalysis, and most probably in codon-anticodon
interaction. The poorer pseudoaromatic cross-modulation and stronger stacking
in ssDNA compared to ssRNA (as manifested in their respective pKas)
may have considerable implications in serving the purpose of DNA as a carrier
of the genetic code (its almost error-free replication and transcription).
On the other hand, the conformational flexibility of ssRNA allows it to
create different scaffolds and nascent folded states with the characteristic
tunable dielectrics, giving variable microenvironments (resulting in to
larger pKa variation), which is manifested in its its dual
biological role in general, as we witness in the translation machninery and catalysis.
Thus in a RNA sequence, P1 Q1
N Q2 P2 , the actual physico-chemical
integrity of N is dictated by the pseudoaromatic
character of both neighboring Q1 and
Q2 , whose properties are further tuned by the electronic
nature of P1 and P2.
Hence, the relative stacking-destacking in any two adjacent nucleotides
will actually set the ON and OFF
switch for the tunability of the pseudoaromatic character of a particular
nucleobase, N. Thus, the pseudoaromatic
character of N can have at least
of variations, depending upon the chemical nature of the neighboring
and Q2, which therefore implies that a given nucleobase
sequence in a polynucleotide chain constitutes an unique extended genetic
code, which can be turned
ON and OFF
depending upon the intrinsic dynamics of folding and unfolding within the
molecule owing to the sequence context or interaction
with an external ligand.
cross-modulation of the pseudoaromatic character in single-stranded RNA
by nearest-neighbor interactions.
Pure and Applied Chemistry, 77 (1), 291-311 (2005).
349. Acharya, S.; Barman, J.; Cheruku, P.; Chatterjee, S.; Acharya, P.; Isaksson, J. and Chattopadhyaya, J. Significant pKa perturbation of Nucleobases Is an Intrinsic Property of the Sequence Context in DNA and RNA. J. Am. Chem. Soc. 2004, 126, 8674-8681.
347. Measurement of Nucleobase p Ka Values in Model Mononucleotides Shows RNA-RNA Duplexes To Be More Stable than DNA-DNA Duplexes. Acharya, P.; Cheruku, P.; Chatterjee, S.; Acharya, S. and Chattopadhyaya, J. J. Am. Chem. Soc; 126, 2862-2869 (2004).
344. Acharya, P.; Acharya, S.; Amirkhanov, N. V.; Cheruku, P.; Földesi, A. and Chattopadhyaya, J. Cross-modulation of the pKa of Nucleobases in a Single Stranded Hexameric-RNA Due to Tandem Electrostatic Nearest-neighbor Interactions. J. Am. Chem. Soc. 2003, 125, 9948-9961.
341. Acharya, P.; Acharya, S.; Földesi, A. and Chattopadhyaya, J. Tandem Electrostatic Effect From the First to the Third Aglycon in the Trimeric RNA Owing to the Nearest-neighbor Stacking. J. Am. Chem. Soc. 2003, 125, 2094-2100.
343. Acharya, P.; Plashkevych, O.; Morita, C.; Yamada, S. and Chattopadhyaya, J. A Repertoire of Pyridinium-Phenyl-Methyl Cross-Talk through a Cascade of Intramolecular Electrostatic Interactions. J. Org. Chem. 2003, 68, 1529-1538.
339. Acharya, S.; Acharya, P.; Földesi, A. and Chattopadhyaya, J. Cross-Modulation of Physicochemical Character of Aglycones in Dinucleoside (3'->5') monophosphates by the Nearest Neighbor Interaction in the Stacked State. J. Am. Chem. Soc. 2002, 124, 13722-13730.
331. Acharya, P. and Chattopadhyaya, J. The Hydrogen Bonding and Hydration of 2'-OH in Adenosine and Adenosine 3'-ethylphosphate. J. Org. Chem. 2002, 67, 1852-1865.
301. Acharya, P.; Trifonova, A.; Thibaudeau, C.; Földesi, A. and Chattopadhyaya, J. The Transmission of the Electronic Character of Guanin-9-yl Drives the Sugar-phosphate Backbone Torsions in Guanosine 3',5'-bisphosphate. Angew Chem. Int. Ed. 1999, 38, 3645-3650.
Chemical Characters of the Internucleotidic Phosphates can Modulate the
Non-enzymatic Reactivity of the Phosphodiester Bonds in RNA. J. Barman,
S. Acharya, S. Chatterjee, Å. Engström and J. Chattopadhyaya,
Org. Biomol. Chem, 4, 928-941 (2006).