Enhancement of intrinsic guanine fluorescence by protonation in DNA of various structures

Understanding the diversity of DNA structure and functions in biology requires tools to study this biomolecule selectively and thoroughly. Fluorescence methods are powerful technique for non-invasive research. Due to the low quantum yield, the intrinsic fluorescence of nucleotides has not been considered for use in the detection and differentiation of nucleic acid bases. Here, we have studied the influence of protonation of nucleotides on their fluorescence properties. We show that protonation of ATP and GTP leads to enhanced intrinsic fluorescence. Fluorescence enhancement at acidic pH has been observed for double-stranded DNA and single-stranded oligonucleotides. The formation of G4 secondary structures apparently protected certain nucleotides from protonation, resulting in less pronounced fluorescence enhancement. Furthermore, acid-induced depurination under protonation was less noticeable in G4 structures than in double-stranded and single-stranded DNA. We show that changes in the intrinsic fluorescence of guanine can be used as a sensitive sensor for changes in the structure of the DNA and for the protonation of specific nucleotides.

Scalable multimodal weak anion exchange chromatographic purification for stable mRNA drug substance

Messenger RNA (mRNA) has emerged as a modality with immense therapeutic potential. Recent innovations in production process of mRNA call for procedures to isolate pure mRNA drug substance (DS) with high yield, high capacity, scalability, and compatibility with GMP production systems. Novel RNA modalities, such as circular RNA (circRNA), have further driven the need for non-affinity capture possibilities which are already widely used in the biopharmaceutical industry, for example, in monoclonal antibody processing. The principle that multimodal ion exchange/hydrogen bonding chromatography can be used to separate mRNA from in vitro transcription components has recently been demonstrated. Here, we apply and refine this approach to be suitable for scalable purification of multiple mRNA constructs with sufficient yields, purity, and stability, for use in mRNA production process. Binding capacity of the PrimaS-modified monolithic chromatographic column for mRNA enabled up to 7 mg/mL product isolation in a single chromatographic run, with 98% recovery and room temperature stability of the eGFP mRNA demonstrated for up to 28 days. This approach is independent of construct size or the presence of polyadenylic acid tail and is applicable for capture of a wide variety of RNAs, including mRNA, self-amplifying RNA, circRNA, and with optimization also smaller RNAs such as transfer RNA and others.

Tracking Topological and Electronic Effects on the Folding and Stability of Guanine-Deficient RNA G-Quadruplexes, Engineered with a New Computational Tool for De Novo Quadruplex Folding

The initial aim of this work was to elucidate the mutual influence of different single-stranded segments (loops and caps) on the thermodynamic stability of RNA G-quadruplexes. To this end, we used a new NAB-GQ-builder software program, to construct dozens of two-tetrad G-quadruplex topologies, based on a designed library of sequences. Then, to probe the sequence–morphology–stability relationships of the designed topologies, we performed molecular dynamics simulations. Their results provide guidance for the design of G-quadruplexes with balanced structures, and in turn programmable physicochemical properties for applications as biomaterials. Moreover, by comparative examinations of the single-stranded segments of three oncogene promoter G-quadruplexes, we assess their druggability potential for future therapeutic strategies. Finally, on the basis of a thorough analysis at the quantum mechanical level of theory on a series of guanine assemblies, we demonstrate how a valence tautomerism, triggered by a coordination of cations, initiates the process of G-quadruplex folding, and we propose a sequential folding mechanism, otherwise dictated by the cancellation of the dipole moments on guanines.

Multidimensional insights into the repeated electromagnetic field stimulation and biosystems interaction in aging and age-related diseases

We provide a multidimensional sequence of events that describe the electromagnetic field (EMF) stimulation and biological system interaction. We describe this process from the quantum to the molecular, cellular, and organismal levels. We hypothesized that the sequence of events of these interactions starts with the oscillatory effect of the repeated electromagnetic stimulation (REMFS). These oscillations affect the interfacial water of an RNA causing changes at the quantum and molecular levels that release protons by quantum tunneling. Then protonation of RNA produces conformational changes that allow it to bind and activate Heat Shock Transcription Factor 1 (HSF1). Activated HSF1 binds to the DNA expressing chaperones that help regulate autophagy and degradation of abnormal proteins. This action helps to prevent and treat diseases such as Alzheimer’s and Parkinson’s disease (PD) by increasing clearance of pathologic proteins. This framework is based on multiple mathematical models, computer simulations, biophysical experiments, and cellular and animal studies. Results of the literature review and our research point towards the capacity of REMFS to manipulate various networks altered in aging (Reale et al. PloS one 9, e104973, 2014), including delay of cellular senescence (Perez et al. 2008, Exp Gerontol 43, 307-316) and reduction in levels of amyloid-β peptides (Aβ) (Perez et al. 2021, Sci Rep 11, 621). Results of these experiments using REMFS at low frequencies can be applied to the treatment of patients with age-related diseases. The use of EMF as a non-invasive therapeutic modality for Alzheimer’s disease, specifically, holds promise. It is also necessary to consider the complicated and interconnected genetic and epigenetic effects of the REMFS-biological system’s interaction while avoiding any possible adverse effects.

Understanding the role of non-Watson-Crick base pairs in DNA-protein recognition: Structural and energetic aspects using crystallographic database analysis and quantum chemical calculation

Specific recognition of DNA base sequences by proteins is vital for life-cycles of all organisms. In a large number of crystal structures of protein-DNA complexes, DNA conformation significantly deviates from the canonical B-DNA structure. A key question is whether such alternate conformations exist prior to protein binding and one is selected for complexation or the structure observed is induced by protein binding. Non-canonical base pairs, such as Hoogsteen base pairs, are often observed in crystal structures of protein-DNA complexes. We decided to explore whether the occurrence of such non-canonical base pairs in protein-DNA complexes is induced by the protein or is selected from pre-existing conformations. Detailed quantum chemical calculations with dispersion-corrected density functional theory (DFT-D) indicated that most of the non-canonical base pairs with DNA bases are stable even in the absence of the interacting amino acids. However, the G:G Hoogsteen base pair, which also appears in the telomere structure, appears to be unstable in the absence of other stabilizing agents, such as positively charged amino acids. Thus, the stability of many of the non-canonical base pair containing duplexes may be close to the canonical B-DNA structure and hence energetically accessible in the ground state; suggesting that the selection from pre-existing conformations may be an important mechanism for observed non-canonical base pairs in protein-DNA complexes.

Collision-induced dissociation of homodimeric and heterodimeric radical cations of 9-methylguanine and 9-methyl-8-oxoguanine: correlation between intra-base pair proton transfer originating from the N1-H at a Watson-Crick edge and non-statistical dissociation

It has been shown previously in protonated, deprotonated and ionized guanine-cytosine base pairs that intra-base pair proton transfer from the N1-H at the Watson-Crick edge of guanine to the complementary nucleobase prompts non-statistical dissociation of the base-pair system, and the dissociation of a proton-transferred base-pair structure is kinetically more favored than that of the starting, conventional base-pair structure. However, the fundamental chemistry underlying this anomalous and intriguing kinetics has not been completely revealed, which warrants the examination of more base-pair systems in different structural contexts in order to derive a generalized base-pair structure-kinetics correlation. The purpose of the present work is to expand the investigation to the non-canonical homodimeric and heterodimeric radical cations of 9-methylguanine (9MG) and 9-methyl-8-oxoguanine (9MOG), i.e., [9MG·9MG]˙⁺, [9MOG·9MG]˙⁺ and [9MOG·9MOG]˙⁺. Experimentally, collision-induced dissociation tandem mass spectrometry coupled with an electrospray ionization (ESI) source was used for the formation of base-pair radical cations, followed by detection of dissociation product ions and cross sections in the collisions with Xe gas under single ion-molecule collision conditions and as a function of the center-of-mass collision energy. Computationally, density functional theory and coupled cluster theory were used to calculate and identify probable base-pair structures and intra-base pair proton transfer and hydrogen transfer reactions, followed by kinetics modeling to explore the properties of dissociation transition states and kinetic factors. The significance of this work is twofold: it provides insight into base-pair opening kinetics in three biologically-important, non-canonical systems upon oxidative and ionization damage; and it links non-statistical dissociation to intra-base pair proton-transfer originating from the N1-H at the Watson-Crick edge of 8-oxoguanine, enhancing understanding towards the base-pair fragmentation assisted by proton transfer.

The Effects of Potential and pH on the Adsorption of Guanine on the Au(111) Electrode

The adsorption of nucleobase at a gold electrode has been a model system to study the interaction between biomolecule and metal, which is relevant to the development of sensors and molecular electronics. The current study has employed in situ scanning tunneling microscopy (STM) and voltammetry to investigate the adsorption configuration and spatial structure of guanine (G) on a well-defined Au(111) electrode in perchloric acid (HClO₄) and neutral phosphate buffer solution (PBS) containing 50 μM G. Potential control had a profound effect on the adsorption of the G molecule on the Au(111) electrode. No adsorption of G was observed at a potential more negative than 0 V in HClO₄ and -0.2 V (versus Ag/AgCl) in PBS; shifting potential positively triggered a rapid adsorption of G to yield a well-ordered G array. Different spatial structures of G admolecules were imaged with STM in HClO₄ and PBS, suggesting that ions in the electrolyte were important in this adsorption event. Shifting potential positively caused a more compact G adlayer with molecules adopting a tilted orientation. Meanwhile, G molecules continued to deposit on the Au(111) electrode leading to a multilayer G film. These processes were reversible to the potential modulation. G admolecules on the Au(111) electrode could be irreversibly oxidized in 0.1 M PBS, which resulted in a prominent peak at 0.74 V in the voltammogram. This oxidation process could be used to analyze the G molecule in a sample.

Systematic Optimization of Ambient Ionization Ion Mobility Mass Spectrometry for Rapid Separation of Isomers

Current methods typically used for metabolite screening and disease diagnosis often require extensive sample preparation, which increases analysis time and associated costs. While ambient ionization techniques enable the analysis of various samples in complex matrices with little or no sample preparation in a short time (typically within a minute), their reduced selectivity, even when coupled with high-resolution mass spectrometers, limits their application in certain fields. In this study, we have optimized the coupling of paper spray (PS) and leaf spray (LS) ambient ionization techniques with a commercially available ion mobility mass spectrometer (IM-MS) and demonstrated the separation of geometric and constitutional isomers. Ambient ionization techniques allow simultaneous introduction and ionization of samples, while background noise and matrix interference from paper and leaf substrates are filtered out by IM separation, resulting in high sensitivity and selectivity of the PS-IM-MS and LS-IM-MS workflows. In addition, we introduced a novel approach to perform single-field collision cross section (CCS) measurements, which resulted in CCS values that differ by 0.15% and 0.25% from traditional stepped-field and single-field methods, respectively. In addition, we used advanced computational tools to confidently identify analyte structures by comparing CCS values from experimental IM measurements and theoretical calculations. These results suggest that the coupling of ambient ionization methods with ion mobility techniques enables rapid, sensitive, and highly selective analysis that can be used in different fields, such as agrochemical screening and disease diagnostics.

Tackling a Curious Case: Generation of Charge-Tagged Guanosine Radicals by Gas-Phase Electron Transfer and Their Characterization by UV-vis Photodissociation Action Spectroscopy and Theory

We report the generation of gas-phase riboguanosine radicals that were tagged at ribose with a fixed-charge 6-(trimethylammonium)hexane-1-aminocarbonyl group. The radical generation relied on electron transfer from fluoranthene anion to noncovalent dibenzocrown-ether dication complexes which formed nucleoside cation radicals upon one-electron reduction and crown-ether ligand loss. The cation radicals were characterized by collision-induced dissociation (CID), photodissociation (UVPD), and UV-vis action spectroscopy. Identification of charge-tagged guanosine radicals was challenging because of spontaneous dissociations by loss of a hydrogen atom and guanine that occurred upon storing the ions in the ion trap without further excitation. The loss of H proceeded from an exchangeable position on N-7 in guanine that was established by deuterium labeling and was the lowest energy dissociation of the guanosine radicals according to transition-state energy calculations. Rate constant measurements revealed an inverse isotope effect on the loss of either hydrogen or deuterium with rate constants k_H = 0.25-0.26 s^-1 and k_D = 0.39-0.54 s^-1. We used time-dependent density functional theory calculations, including thermal vibronic effects, to predict the absorption spectra of several protomeric radical isomers. The calculated spectra of low-energy N-7-H guanine-radical tautomers closely matched the action spectra. Transition-state-theory calculations of the rate constants for the loss of H-7 and guanine agreed with the experimental rate constants for a narrow range of ion effective temperatures. Our calculations suggest that the observed inverse isotope effect does not arise from the isotope-dependent differences in the transition-state energies. Instead, it may be caused by the dynamics of post-transition-state complexes preceding the product separation.

Frequency and hydrogen bonding of nucleobase homopairs in small molecule crystals

We used the high resolution and accuracy of the Cambridge Structural Database (CSD) to provide detailed information regarding base pairing interactions of selected nucleobases. We searched for base pairs in which nucleobases interact with each other through two or more hydrogen bonds and form more or less planar structures. The investigated compounds were either free forms or derivatives of adenine, guanine, hypoxanthine, thymine, uracil and cytosine. We divided our findings into categories including types of pairs, protonation patterns and whether they are formed by free bases or substituted ones. We found base pair types that are exclusive to small molecule crystal structures, some that can be found only in RNA containing crystal structures and many that are native to both environments. With a few exceptions, nucleobase protonation generally followed a standard pattern governed by pKa values. The lengths of hydrogen bonds did not depend on whether the nucleobases forming a base pair were charged or not. The reasons why particular nucleobases formed base pairs in a certain way varied significantly.

An IRMPD spectroscopic and computational study of protonated guanine-containing mismatched base pairs in the gas phase

Infrared multiple photon dissociation spectroscopy has been used to probe the structures of the three protonated base-pair mismatches containing 9-ethylguanine (9eG) in the gas phase. Some of these protonated base-pairs have been identified in RNA.

Guanine–adenine interactions in DNA tetranucleotide cation radicals revealed by UV/vis photodissociation action spectroscopy and theory

Hydrogen-rich cation radicals (GATT + 2H)⁺˙ and (AGTT + 2H)⁺˙ represent oligonucleotide models of charged hydrogen atom adducts to DNA.

Single-stranded DNA structural diversity: TAGGGT from monomers to dimers to tetramer formation

RATIONALE

DNA quadruplex structures have emerged as novel drug targets due to their role in preventing abnormal gene transcription and maintaining telomere stability. Trapped Ion Mobility Spectrometry-Mass Spectrometry (TIMS-MS), combined with theoretical modeling, is a powerful tool for studying the kinetic intermediates of DNA complexes formed in solution and interrogated in the gas phase after desolvation.

METHODS

A TAGGGT ssDNA sequence was purchased and studied in 10 mM ammonium acetate using nanospray electrospray ionization (nESI)-TIMS-MS in positive and negative ion mode. Collisional cross section (CCS) profiles were measured using internal calibration (Tune Mix). Theoretical structures were proposed based on molecular dynamics, charge location and geometry optimization for the most intense IMS bands based on the number of TAGGGT units, adduct form and charge states.

RESULTS

A distribution of monomeric, dimeric and tetrameric TAGGGT structures were formed in solution and separated in the gas phase based on their mobility and m/z value (e.g., [M + 2H]⁺² , [2M + 3H]⁺³ , [M - 2H]^-2 , [2M - 3H]^-3 , [4M + 4H]⁺⁴ , [4M + 3H + NH₄ ]⁺⁴ , [4M + 2H + 2NH₄ ]⁺⁴ and [4M + H + 3NH₄ ]⁺⁴ ). The high mobility resolution of the TIMS-MS analyzer permitted the observation of multiple CCS bands per molecular ion form. Comparison with theoretical candidate structures suggests that monomeric TAGGGT species are stabilized by A-T and G⁺ -G interactions, with the size of the conformer influenced by the proton location. In the case of the TAGGGT quadruplex, the protonated species displayed a broad CCS distribution, while six discrete conformers were stabilized by the presence of ammonium ions (n = 1-3).

CONCLUSIONS

This is the first observation of multiple conformations of TAGGGT complexes (n = 1, 2 and 4) in 10 mM ammonium acetate. Candidate structures with intramolecular interactions of the form of G⁺ -G and traditional A-T base pairing agreed with the experimental trends. Our results demonstrate the structural diversity of TAGGGT monomers, dimers and tetramers in the gas phase beyond the previously reported solution structure, using 10 mM ammonium acetate to replicate biological conditions.

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

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

Evidence for Hidden Involvement of N3-Protonated Guanine in RNA Structure and Function

Charged nucleobases have been found to occur in several known RNA molecules and are considered essential for their structure and function. The mechanism of their involvement is however not yet fully understood. Revelation of the role of N7-protonated guanine, in modulating the geometry and stability of noncanonical base pairs formed through its unprotonated edges [Watson-Crick (WC) and sugar], has triggered the need to evaluate the feasibility of similar roles of other protonated nucleobases [Halder et al., Phys Chem Chem Phys, 2015, 17, 26249]. In this context, N3 protonation of guanine makes an interesting case as its influence on the charge distribution of the WC edge is similar to that of N7 protonation, though its thermodynamic cost of protonation is significantly higher. In this work, we have carried out structural bioinformatics analyses and quantum mechanics-based calculations to show that N3 protonation of guanine may take place in a cellular environment, at least in the G:C W:W Trans and G:G W:H Cis base pairs. Our results provide a reasonable starting point for future investigations in order to address the larger mechanistic question.

Hydrogen-Rich Cation Radicals of DNA Dinucleotides: Generation and Structure Elucidation by UV-Vis Action Spectroscopy

Hydrogen-rich DNA dinucleotide cation radicals (dGG + 2H)^+•, (dCG + 2H)^+•, and (dGC + 2H)^+• represent transient species comprising protonated and hydrogen atom adducted nucleobase rings that serve as models for proton and radical migrations in ionized DNA. These DNA cation radicals were generated in the gas phase by electron-transfer dissociation of dinucleotide dication-crown-ether complexes and characterized by UV-vis photodissociation action spectra, ab initio calculations of structures and relative energies, and time-dependent density functional theory calculations of UV-vis absorption spectra. Theoretical calculations indicate that (dGG + 2H)^+• cation radicals formed by electron transfer underwent an exothermic conformational collapse that was accompanied by guanine ring stacking and facile internucleobase hydrogen atom transfer, forming 3'-guanine C-8-H radicals. In contrast, exothermic hydrogen transfer from the 5'-cytosine radical onto the guanine ring in (dCG + 2H)^+• was kinetically hampered, resulting in the formation of a mixture of 5'-cytosine and 3'-guanine radicals. Conformational folding and nucleobase stacking were energetically unfavorable in (dGC + 2H)^+• that retained its structure of a 3'-cytosine radical, as formed by one-electron reduction of the dication. Hydrogen-rich guanine (G + H)^• and cytosine (C + H)^• radicals were calculated to have vastly different basicities in water, as illustrated by the respective p K_a values of 20.0 and 4.6, which is pertinent to their different abilities to undergo proton-transfer reactions in solution.

Protonated nucleobases are not fully ionized in their chloride salt crystals and form metastable base pairs further stabilized by the surrounding anions

This paper presents experimental charge-density studies of cytosinium chloride, adeninium chloride hemihydrate and guaninium dichloride crystals based on ultra-high-resolution X-ray diffraction data and extensive theoretical calculations. The results confirm that the cohesive energies of the studied systems are dominated by contributions from intermolecular electrostatic interactions, as expected for ionic crystals. Electrostatic interaction energies (Ees) usually constitute 95% of the total interaction energy. The Ees energies in this study were several times larger in absolute value when compared, for example, with dimers of neutral nucleobases. However, they were not as large as some theoretical calculations have predicted. This was because the molecules appeared not to be fully ionized in the studied crystals. Apart from charge transfer from chlorine to the protonated nucleobases, small but visible charge redistribution within the nucleobase cations was observed. Some dimers of singly protonated bases in the studied crystals, namely a cytosinium-cytosinium trans sugar/sugar edge pair and an adeninium-adeninium trans Hoogsteen/Hoogsteen edge pair, exhibited attractive interactions (negative values of Ees) or unusually low repulsion despite identical molecular charges. The pairs are metastable as a result of strong hydrogen bonding between bases which overcompensates the overall cation-cation repulsion, the latter being weakened due to charge transfer and molecular charge-density polarization.

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

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

Structural landscape of base pairs containing post-transcriptional modifications in RNA

Base pairs involving post-transcriptionally modified nucleobases are believed to play important roles in a wide variety of functional RNAs. Here we present our attempts toward understanding the structural and functional role of naturally occurring modified base pairs using a combination of X-ray crystal structure database analysis, sequence analysis, and advanced quantum chemical methods. Our bioinformatics analysis reveals that despite their presence in all major secondary structural elements, modified base pairs are most prevalent in tRNA crystal structures and most commonly involve guanine or uridine modifications. Further, analysis of tRNA sequences reveals additional examples of modified base pairs at structurally conserved tRNA regions and highlights the conservation patterns of these base pairs in three domains of life. Comparison of structures and binding energies of modified base pairs with their unmodified counterparts, using quantum chemical methods, allowed us to classify the base modifications in terms of the nature of their electronic structure effects on base-pairing. Analysis of specific structural contexts of modified base pairs in RNA crystal structures revealed several interesting scenarios, including those at the tRNA:rRNA interface, antibiotic-binding sites on the ribosome, and the three-way junctions within tRNA. These scenarios, when analyzed in the context of available experimental data, allowed us to correlate the occurrence and strength of modified base pairs with their specific functional roles. Overall, our study highlights the structural importance of modified base pairs in RNA and points toward the need for greater appreciation of the role of modified bases and their interactions, in the context of many biological processes involving RNA.

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

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