Higher order structural effects stabilizing the reverse Watson-Crick Guanine-Cytosine base pair in functional RNAs

Unveiling structural and energetic characterization of the emissive RNA alphabet anchored in the methylthieno[3,4-d]pyrimidine heterocycle core

This study presents a comprehensive theoretical exploration of the fluorescent non-natural emissive nucleobases- mthA, mthG, mthC, and mthU derived from the methylthieno[3,4-d]pyrimidine heterocycle. Our calculations, aligning with experimental findings, reveal that these non-natural bases exert minimal influence on the geometry of classical Watson-Crick base pairs within an RNA duplex, maintaining H-bonding akin to natural bases. In terms of energy, the impact of the modified bases, but for mthG, is also found to be little significant. We delved into an in-depth analysis of the photophysical properties of these non-natural bases. This investigation unveiled a correlation between their absorption/emission peaks and the substantial impact of the modification on the energy levels of the highest unoccupied molecular orbitals (HOMO) and the lowest unoccupied molecular orbital (LUMO). Notably, this alteration in energy levels resulted in a significant reduction of the HOMO-LUMO gap, from approximately 5.4-5.5 eV in the natural bases, to roughly 3.9-4.7 eV in the modified bases. This shift led to a consequential change in absorption and emission spectra towards longer wavelengths, elucidating their bathochromic shift.

An Atlas of the base inter-RNA stacks involved in bacterial translation

Nucleobase-specific noncovalent interactions play a crucial role in translation. Herein, we provide a comprehensive analysis of the stacks between different RNA components in the crystal structures of the bacterial ribosome caught at different translation stages. Analysis of tRNA||rRNA stacks reveals distinct behaviour; both the A-and E-site tRNAs exhibit unique stacking patterns with 23S rRNA bases, while P-site tRNAs stack with 16S rRNA bases. Furthermore, E-site stacks exhibit diverse face orientations and ring topologies-rare for inter-chain RNA interactions-with higher average interaction energies than A or P-site stacks. This suggests that stacking may be essential for stabilizing tRNA progression through the E-site. Additionally, mRNA||rRNA stacks reveal other geometries, which depend on the tRNA binding site, whereas 16S rRNA||23S rRNA stacks highlight the importance of specific bases in maintaining the integrity of the translational complex by linking the two rRNAs. Furthermore, tRNA||mRNA stacks exhibit distinct geometries and energetics at the E-site, indicating their significance during tRNA translocation and elimination. Overall, both A and E-sites display a more diverse distribution of inter-RNA stacks compared to the P-site. Stacking interactions in the active ribosome are not simply accidental byproducts of biochemistry but are likely invoked to compensate and support the integrity and dynamics of translation.

Designing a multi-epitope candidate vaccine by employing immunoinformatics approaches to control African swine fever spread

The African swine fever virus has been circulating for decades and is highly infectious, often fatal to farmed and wild pigs. There is currently no approved vaccine or treatment for the disease, making prevention even more difficult. Therefore, vaccine development is necessary and urgent to limit the consequences of ASF and ensure the food chain and sustainability of the swine industry. This research study was conducted to design a multi-epitope vaccine for controlling veterinary diseases caused by the African swine fever virus. We employed the immunoinformatics approaches to reveal 37 epitopes from different viral proteins of ASFV. These epitopes were linked to adjuvants and linkers to form a full-fledged immunogenic vaccine construct. The tertiary structure of the final vaccine was predicted using a deep-learning approach. The molecular docking and molecular dynamics predicted stable interactions between the vaccine and immune receptor TLR5 of Sus scrofa (Pig). The MD simulation studies reflect that the calculated parameters like RMSD, RMSF, number of hydrogen bonds, and finally, the buried interface surface area for the complex remained stable throughout the simulation time. This analysis suggests the stability of interface interactions between the TLR5 and the multi-epitope vaccine construct. Further, the physiochemical analysis demonstrated that our designed vaccine construct was expected to have high stability and prolonged half-life time in mammalian cells. Traditional vaccine design experiments require significant time and financial input from the development stage to the final product. Studies like this can assist in accelerating vaccine development while minimizing the cost.Communicated by Ramaswamy H. Sarma.

Dynamical characterization and multiple unbinding paths of two PreQ1 ligands in one pocket

Riboswitches naturally regulate gene expression in bacteria by binding to specific small molecules. Class 1 preQ1 riboswitch aptamer is an important model not only for RNA folding but also as a target for designing small molecule antibiotics due to its well-known minimal aptamer domain. Here, we ran a total of 62.4 μs conventional and enhanced-sampling molecular dynamics (MD) simulations to characterize the determinants underlying the binding of the preQ1-II riboswitch aptamer to two preQ1 ligands in one binding pocket. Decomposition of binding free energy suggested that preQ1 ligands at α and β sites interact with four nucleotides (G5, C17, C18, and A30) and two nucleotides (A12 and C31), respectively. Mg2+ ions play a crucial role in both stabilizing the binding pocket and facilitating ligand binding. The flexible preQ1 ligand at the β site leads to the top of the binding pocket loosening and thus pre-organizes the riboswitch for ligand entry. Enhanced sampling simulations further revealed that the preQ1 ligand at the α site unbinds through two orthogonal pathways, which are dependent on whether or not a β site preQ1 ligand is present. One of the two preQ1 ligands has been identified in the binding pocket, which will aid to identify the second preQ1 Ligand. Our work provides new information for designing robust ligands.

Structural and Energetic Features of Base-Base Stacking Contacts in RNA

Nucleobase π-π stacking is one of the crucial organizing interactions within three-dimensional (3D) RNA architectures. Characterizing the structural variability of these contacts in RNA crystal structures will help delineate their subtleties and their role in determining function. This analysis of different stacking geometries found in RNA X-ray crystal structures is the largest such survey to date; coupled with quantum-mechanical calculations on typical representatives of each possible stacking arrangement, we determined the distribution of stacking interaction energies. A total of 1,735,481 stacking contacts, spanning 359 of the 384 theoretically possible distinct stacking geometries, were identified. Our analysis reveals preferential occurrences of specific consecutive stacking arrangements in certain regions of RNA architectures. Quantum chemical calculations suggest that 88 of the 359 contacts possess intrinsically stable stacking geometries, whereas the remaining stacks require the RNA backbone or surrounding macromolecular environment to force their formation and maintain their stability. Our systematic analysis of π-π stacks in RNA highlights trends in the occurrence and localization of these noncovalent interactions and may help better understand the structural intricacies of functional RNA-based molecular architectures.

Occurrence and stability of anion-π interactions between phosphate and nucleobases in functional RNA molecules

We present a systematic structural and energetic characterization of phosphate(OP)-nucleobase anion…π stacking interactions in RNAs. We observed OP-nucleobase stacking contacts in a variety of structural motifs other than regular helices and spanning broadly diverse sequence distances. Apart from the stacking between a phosphate and a guanine or a uracil two-residue upstream in specific U-turns, such interactions in RNA have been scarcely characterized to date. Our QM calculations showed an energy minimum at a distance between the OP atom and the nucleobase plane centroid slightly below 3 Å for all the nucleobases. By sliding the OP atom over the nucleobase plane we localized the optimal mutual positioning of the stacked moieties, corresponding to an energy minimum below -6 kcal•mol-1, for all the nucleobases, consistently with the projections of the OP atoms over the different π-rings we observed in experimental occurrences. We also found that the strength of the interaction clearly correlates with its electrostatic component, pointing to it as the most relevant contribution. Finally, as OP-uracil and OP-guanine interactions represent together 86% of the instances we detected, we also proved their stability under dynamic conditions in model systems simulated by state-of-the art DFT-MD calculations.

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

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

Halogen-Bonded Guanine Base Pairs, Quartets and Ribbons

Halogen bonding is studied in different structures consisting of halogenated guanine DNA bases, including the Hoogsteen guanine–guanine base pair, two different types of guanine ribbons (R-I and R-II) consisting of two or three monomers, and guanine quartets. In the halogenated base pairs (except the Cl-base pair, which has a very non-planar structure with no halogen bonds) and R-I ribbons (except the At trimer), the potential N-X•••O interaction is sacrificed to optimise the N-X•••N halogen bond. In the At trimer, the astatines originally bonded to N1 in the halogen bond donating guanines have moved to the adjacent O6 atom, enabling O-At•••N, N-At•••O, and N-At•••At halogen bonds. The brominated and chlorinated R-II trimers contain two N-X•••N and two N-X•••O halogen bonds, whereas in the iodinated and astatinated trimers, one of the N-X•••N halogen bonds is lost. The corresponding R-II dimers keep the same halogen bond patterns. The G-quartets display a rich diversity of symmetries and halogen bond patterns, including N-X•••N, N-X•••O, N-X•••X, O-X•••X, and O-X•••O halogen bonds (the latter two facilitated by the transfer of halogens from N1 to O6). In general, halogenation decreases the stability of the structures. However, the stability increases with the increasing atomic number of the halogen, and the At-doped R-I trimer and the three most stable At-doped quartets are more stable than their hydrogenated counterparts. Significant deviations from linearity are found for some of the halogen bonds (with halogen bond angles around 150°).

Occurrence and stability of lone pair-π and OH-π interactions between water and nucleobases in functional RNAs

We identified over 1000 instances of water-nucleobase stacking contacts in a variety of RNA molecules from a non-redundant set of crystal structures with resolution ≤3.0 Å. Such contacts may be of either the lone pair-π (lp-π) or the OH-π type, in nature. The distribution of the distances of the water oxygen from the nucleobase plane peaks at 3.5 Å for A, G and C, and approximately at 3.1-3.2 Å for U. Quantum mechanics (QM) calculations confirm, as expected, that the optimal energy is reached at a shorter distance for the lp-π interaction as compared to the OH-π one (3.0 versus 3.5 Å). The preference of each nucleobase for either type of interaction closely correlates with its electrostatic potential map. Furthermore, QM calculations show that for all the nucleobases a favorable interaction, of either the lp-π or the OH-π type, can be established at virtually any position of the water molecule above the nucleobase skeleton, which is consistent with the uniform projection of the OW atoms over the nucleobases ring we observed in the experimental occurrences. Finally, molecular dynamics simulations of a model system for the characterization of water-nucleobase stacking contacts confirm the stability of these interactions also under dynamic conditions.

Pairing interactions between nucleobases and ligands in aptamer:ligand complexes of riboswitches: crystal structure analysis, classification, optimal structures, and accurate interaction energies

In the present work, 67 crystal structures of the aptamer domains of RNA riboswitches are chosen for analysis of the structure and strength of hydrogen bonding (pairing) interactions between nucleobases constituting the aptamer binding pockets and the bound ligands. A total of 80 unique base:ligand hydrogen-bonded pairs containing at least two hydrogen bonds were identified through visual inspection. Classification of these contacts in terms of the interacting edge of the aptamer nucleobase revealed that interactions involving the Watson-Crick edge are the most common, followed by the sugar edge of purines and the Hoogsteen edge of uracil. Alternatively, classification in terms of the chemical constitution of the ligand yields five unique classes of base:ligand pairs: base:base, base:amino acid, base:sugar, base:phosphate, and base:other. Further, quantum mechanical (QM) geometry optimizations revealed that 67 out of 80 pairs exhibit stable geometries and optimal deviations from their macromolecular crystal occurrences. This indicates that these contacts are well-defined RNA aptamer:ligand interaction motifs. QM calculated interaction energies of base:ligand pairs reveal a rich hydrogen bonding landscape, ranging from weak interactions (base:other, -3 kcal/mol) to strong (base:phosphate, -48 kcal/mol) contacts. The analysis was further extended to study the biological importance of base:ligand interactions in the binding pocket of the tetrahydrofolate riboswitch and thiamine pyrophosphate riboswitch. Overall, our study helps in understanding the structural and energetic features of base:ligand pairs in riboswitches, which could aid in developing meaningful hypotheses in the context of RNA:ligand recognition. This can, in turn, contribute toward current efforts to develop antimicrobials that target RNAs.

Structural and Energetic Impact of Non-natural 7-Deaza-8-azaguanine, 7-Deaza-8-azaisoguanine, and Their 7-Substituted Derivatives on Hydrogen-Bond Pairing with Cytosine and Isocytosine

The impact of 7-deaza-8-azaguanine (DAG) and 7-deaza-8-azaisoguanine (DAiG) modifications on the geometry and stability of the G:C Watson-Crick (cWW) base pair and the G:iC and iG:C reverse Watson-Crick (tWW) base pairs has been characterized theoretically. In addition, the effect on the same base pairs of seven C7-substituted DAG and DAiG derivatives, some of which have been previously experimentally characterized, has been investigated. Calculations indicate that all of these modifications have a negligible impact on the geometry of the above base pairs, and that modification of the heterocycle skeleton has a small impact on the base-pair interaction energies. Instead, base-pair interaction energies are dependent on the nature of the C7 substituent. For the 7-substituted DAG-C cWW systems, a linear correlation between the base-pair interaction energy and the Hammett constant of the 7-substituent is found, with higher interaction energies corresponding to more electron-withdrawing substituents. Therefore, the explored modifications are expected to be accommodated in both parallel and antiparallel nucleic acid duplexes without perturbing their geometry, while the strength of a base pair (and duplex) featuring a DAG modification can, in principle, be tuned by incorporating different substituents at the C7 position.

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.

Type-II tRNAs and Evolution of Translation Systems and the Genetic Code

Because tRNA is the core biological intellectual property that was necessary to evolve translation systems, tRNAomes, ribosomes, aminoacyl-tRNA synthetases, and the genetic code, the evolution of tRNA is the core story in evolution of life on earth. We have previously described the evolution of type-I tRNAs. Here, we use the same model to describe the evolution of type-II tRNAs, with expanded V loops. The models are strongly supported by inspection of typical tRNA diagrams, measuring lengths of V loop expansions, and analyzing the homology of V loop sequences to tRNA acceptor stems. Models for tRNA evolution provide a pathway for the inanimate-to-animate transition and for the evolution of translation systems, the genetic code, and cellular life.

An intricate balance of hydrogen bonding, ion atmosphere and dynamics facilitates a seamless uracil to cytosine substitution in the U-turn of the neomycin-sensing riboswitch

The neomycin sensing riboswitch is the smallest biologically functional RNA riboswitch, forming a hairpin capped with a U-turn loop-a well-known RNA motif containing a conserved uracil. It was shown previously that a U→C substitution of the eponymous conserved uracil does not alter the riboswitch structure due to C protonation at N3. Furthermore, cytosine is evolutionary permitted to replace uracil in other U-turns. Here, we use molecular dynamics simulations to study the molecular basis of this substitution in the neomycin sensing riboswitch and show that a structure-stabilizing monovalent cation-binding site in the wild-type RNA is the main reason for its negligible structural effect. We then use NMR spectroscopy to confirm the existence of this cation-binding site and to demonstrate its effects on RNA stability. Lastly, using quantum chemical calculations, we show that the cation-binding site is altering the electronic environment of the wild-type U-turn so that it is more similar to the cytosine mutant. The study reveals an amazingly complex and delicate interplay between various energy contributions shaping up the 3D structure and evolution of nucleic acids.

Theoretical characterization of sulfur-to-selenium substitution in an emissive RNA alphabet: impact on H-bonding potential and photophysical properties

We employ density functional theory (DFT) and time-dependent DFT (TDDFT) calculations to investigate the structural, energetic and optical properties of a new computationally designed RNA alphabet, where the nucleobases, ^tsA, ^tsG, ^tsC, and ^tsU (ts-bases), have been derived by replacing sulfur with selenium in the previously reported tz-bases, based on the isothiazolo[4,3-d]pyrimidine heterocycle core. We find out that the modeled non-natural bases have minimal impact on the geometry and energetics of the classical Watson-Crick base pairs, thus potentially mimicking the natural bases in a RNA duplex in terms of H-bonding. In contrast, our calculations indicate that H-bonded base pairs involving the Hoogsteen edge of purines are destabilized as compared to their natural counterparts. We also focus on the photophysical properties of the non-natural bases and correlate their absorption/emission peaks to the strong impact of the modification on the energy of the lowest unoccupied molecular orbital. It is indeed stabilized by roughly 1.1-1.6 eV as compared to the natural analogues, resulting in a reduction of the gap between the highest occupied and the lowest unoccupied molecular orbital from 5.3-5.5 eV in the natural bases to 3.9-4.2 eV in the modified ones, with a consequent bathochromic shift in the absorption and emission spectra. Overall, our analysis clearly indicates that the newly modelled ts-bases are expected to exhibit better fluorescent properties as compared to the previously reported tz-bases, while retaining similar H-bonding properties. In addition, we show that a new RNA alphabet based on size-extended benzo-homologated ts-bases can also form stable Watson-Crick base pairs with the natural complementary nucleobases.

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.

Energetics and dynamics of the non-natural fluorescent 4AP:DAP base pair

Quantum mechanics and molecular dynamics methods are used to compare the non-natural 4AP–DAP base pair to natural base pairs.

Occurrence and stability of lone pair-π stacking interactions between ribose and nucleobases in functional RNAs

The specific folding pattern and function of RNA molecules lies in various weak interactions, in addition to the strong base-base pairing and stacking. One of these relatively weak interactions, characterized by the stacking of the O4' atom of a ribose on top of the heterocycle ring of a nucleobase, has been known to occur but has largely been ignored in the description of RNA structures. We identified 2015 ribose-base stacking interactions in a high-resolution set of non-redundant RNA crystal structures. They are widespread in structured RNA molecules and are located in structural motifs other than regular stems. Over 50% of them involve an adenine, as we found ribose-adenine contacts to be recurring elements in A-minor motifs. Fewer than 50% of the interactions involve a ribose and a base of neighboring residues, while approximately 30% of them involve a ribose and a nucleobase at least four residues apart. Some of them establish inter-domain or inter-molecular contacts and often implicate functionally relevant nucleotides. In vacuo ribose-nucleobase stacking interaction energies were calculated by quantum mechanics methods. Finally, we found that lone pair-π stacking interactions also occur between ribose and aromatic amino acids in RNA-protein complexes.

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.

Higher order structures involving post transcriptionally modified nucleobases in RNA

Quantum chemical studies are carried out to understand the structures and stabilities of higher order structures involving post-transcriptionally modified nucleobases in RNA.

How Mg2+ ion and water network affect the stability and structure of non-Watson-Crick base pairs in E. coli loop E of 5S rRNA: a molecular dynamics and reference interaction site model (RISM) study

The non-Watson-Crick (non-WC) base pairs of Escherichia coli loop E of 5S rRNA are stabilized by Mg²⁺ ions through water-mediated interaction. It is important to know the synergic role of Mg²⁺ and the water network surrounding Mg²⁺ in stabilizing the non-WC base pairs of RNA. For this purpose, free energy change of the system is calculated using molecular dynamics (MD) simulation as Mg²⁺ is pulled from RNA, which causes disturbance of the water network. It was found that Mg²⁺ remains hexahydrated unless it is close to or far from RNA. In the pentahydrated form, Mg²⁺ interacts directly with RNA. Water network has been identified by two complimentary methods; MD followed by a density-based clustering algorithm and three-dimensional-reference interaction site model. These two methods gave similar results. Identification of water network around Mg²⁺ and non-WC base pairs gives a clue to the strong effect of water network on the stability of this RNA. Based on sequence analysis of all Eubacteria 5s rRNA, we propose that hexahydrated Mg²⁺ is an integral part of this RNA and geometry of base pairs surrounding it adjust to accommodate the [Formula: see text]. Overall the findings from this work can help in understanding the basis of the complex structure and stability of RNA with non-WC base pairs.

Structural and energetic characterization of the emissive RNA alphabet based on the isothiazolo[4,3-d]pyrimidine heterocycle core

We present theoretical characterization of fluorescent non-natural nucleobases, (tz)A, (tz)G, (tz)C, and (tz)U, derived from the isothiazolo[4,3-d]pyrimidine heterocycle. Consistent with the experimental evidence, our calculations show that the non-natural bases have minimal impact on the geometry and stability of the classical Watson-Crick base pairs, allowing them to accurately mimic natural bases in a RNA duplex, in terms of H-bonding. In contrast, our calculations indicate that H-bonded base pairs involving the Hoogsteen edge are destabilized relative to their natural counterparts. Analysis of the photophysical properties of the non-natural bases allowed us to correlate their absorption/emission peaks to the strong impact of the modification on the energy of the lowest unoccupied molecular orbital, LUMO, which is stabilized by roughly 1.0-1.2 eV relative to the natural analogues, while the highest occupied molecular orbital, HOMO, is not substantially affected. As a result, the HOMO-LUMO gap is reduced from 5.3-5.5 eV in the natural bases to 4.0-4.4 eV in the modified ones, with a consequent bathochromic shift in the absorption and emission spectra.

Structural and Energetic Impact of Non-Natural 7-Deaza-8-Azaadenine and Its 7-Substituted Derivatives on H-Bonding Potential with Uracil in RNA Molecules

Non-natural (synthetic) nucleobases, including 7-ethynyl- and 7-triazolyl-8-aza-7-deazaadenine, have been introduced in RNA molecules for targeted applications, and have been characterized experimentally. However, no theoretical characterization of the impact of these modifications on the structure and energetics of the corresponding H-bonded base pair is available. To fill this gap, we performed quantum mechanics calculations, starting with the analysis of the impact of the 8-aza-7-deaza modification of the adenine skeleton, and we moved then to analyze the impact of the specific substituents on the modified 8-aza-7-deazaadenine. Our analysis indicates that, despite of these severe structural modifications, the H-bonding properties of the modified base pair gratifyingly replicate those of the unmodified base pair. Similar behavior is predicted when the same skeleton modifications are applied to guanine when paired to cytosine. To stress further the H-bonding pairing in the modified adenine-uracil base pair, we explored the impact of strong electron donor and electron withdrawing substituents on the C7 position. Also in this case we found minimal impact on the base pair geometry and energy, confirming the validity of this modification strategy to functionalize RNAs without perturbing its stability and biological functionality.

Quantum chemical benchmark study on 46 RNA backbone families using a dinucleotide unit

We have created a benchmark set of quantum chemical structure-energy data denoted as UpU46, which consists of 46 uracil dinucleotides (UpU), representing all known 46 RNA backbone conformational families. Penalty-function-based restrained optimizations with COSMO TPSS-D3/def2-TZVP ensure a balance between keeping the target conformation and geometry relaxation. The backbone geometries are close to the clustering-means of their respective RNA bioinformatics family classification. High-level wave function methods (DLPNO-CCSD(T) as reference) and a wide-range of dispersion-corrected or inclusive DFT methods (DFT-D3, VV10, LC-BOP-LRD, M06-2X, M11, and more) are used to evaluate the conformational energies. The results are compared to the Amber RNA bsc0χOL3 force field. Most dispersion-corrected DFT methods surpass the Amber force field significantly in accuracy and yield mean absolute deviations (MADs) for relative conformational energies of ∼0.4-0.6 kcal/mol. Double-hybrid density functionals represent the most accurate class of density functionals. Low-cost quantum chemical methods such as PM6-D3H+, HF-3c, DFTB3-D3, as well as small basis set calculations corrected for basis set superposition errors (BSSEs) by the gCP procedure are also tested. Unfortunately, the presently available low-cost methods are struggling to describe the UpU conformational energies with satisfactory accuracy. The UpU46 benchmark is an ideal test for benchmarking and development of fast methods to describe nucleic acids, including force fields.

An atlas of RNA base pairs involving modified nucleobases with optimal geometries and accurate energies

Posttranscriptional modifications greatly enhance the chemical information of RNA molecules, contributing to explain the diversity of their structures and functions. A significant fraction of RNA experimental structures available to date present modified nucleobases, with half of them being involved in H-bonding interactions with other bases, i.e. ‘modified base pairs’. Herein we present a systematic investigation of modified base pairs, in the context of experimental RNA structures. To this end, we first compiled an atlas of experimentally observed modified base pairs, for which we recorded occurrences and structural context. Then, for each base pair, we selected a representative for subsequent quantum mechanics calculations, to find out its optimal geometry and interaction energy. Our structural analyses show that most of the modified base pairs are non Watson–Crick like and are involved in RNA tertiary structure motifs. In addition, quantum mechanics calculations quantify and provide a rationale for the impact of the different modifications on the geometry and stability of the base pairs they participate in.

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

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

Structural stability, acidity, and halide selectivity of the fluoride riboswitch recognition site

Using static and dynamics DFT methods we show that the Mg(2+)/F(-)/phosphate/water cluster at the center of the fluoride riboswitch is stable by its own and, once assembled, does not rely on any additional factor from the overall RNA fold. Further, we predict that the pKa of the water molecule bridging two Mg cations is around 8.4. We also demonstrate that the halide selectivity of the fluoride riboswitch is determined by the stronger Mg-F bond, which is capable of keeping together the cluster. Replacing F(-) with Cl(-) results in a cluster that is unstable under dynamic conditions. Similar conclusions on the structure and energetics of the cluster in the binding pocket of fluoride-inhibited pyrophosphatase suggest that the peculiarity of fluoride is in its ability to establish much stronger metal-halide bonds.

Towards biochemically relevant QM computations on nucleic acids: controlled electronic structure geometry optimization of nucleic acid structural motifs using penalty restraint functions

Recent developments in dispersion-corrected density functional theory methods allow for the first time the description of large fragments of nucleic acids (hundreds of atoms) with an accuracy clearly surpassing the accuracy of common biomolecular force fields. Such calculations can significantly improve the description of the potential energy surface of nucleic acid molecules, which may be useful for studies of molecular interactions and conformational preferences of nucleic acids, as well as verification and parameterization of other methods. The first of such studies, however, demonstrated that successful applications of accurate QM calculations to larger nucleic acid building blocks are hampered by difficulties in obtaining geometries that are biochemically relevant and are not biased by non-native structural features. We present an approach that can greatly facilitate large-scale QM studies on nucleic acids, namely electronic structure geometry optimization of nucleic acid fragments utilizing a penalty function to restrain key internal coordinates with a specific focus on the torsional backbone angles. This work explores the viability of these restraint optimizations for DFT-D3, PM6-D3H and HF-3c optimizations on a set of examples (a UpA dinucleotide, a DNA G-quadruplex and a B-DNA fragment). Evaluation of different penalty function strengths reveals only a minor system-dependency and reasonable restraint values range from 0.01 to 0.05 Eh rad(-2) for the backbone torsions. Restraints are crucial to perform the QM calculations on biochemically relevant conformations in implicit solvation and gas phase geometry optimizations. The reasons for using restrained instead of constrained or unconstrained optimizations are explained and an open-source external optimizer is provided.

Dynamics of riboswitches: Molecular simulations

Riboswitch RNAs play key roles in bacterial metabolism and represent a promising new class of antibiotic targets for treatment of infectious disease. While many studies of riboswitches have been performed, the exact mechanism of riboswitch operation is still not fully understood at the atomistic level of detail. Molecular dynamics simulations are useful for interpreting existing experimental data and producing predictions for new experiments. Here, a wide range of computational studies on riboswitches is reviewed. By elucidating the key principles of riboswitch operation, computation may aid in the effort to design more specific antibiotics with affinities greater than those of the native ligand. Such a detailed understanding may be required to improve efficacy and reduce side effects. These studies are laying the groundwork for understanding the action mechanism of new compounds that inhibit riboswitch activity. Future directions such as magnesium effects, large-scale conformational changes, expression platforms and co-transcriptional folding are also discussed. This article is part of a Special Issue entitled: Riboswitches.