On the role of Hoogsteen:Hoogsteen interactions in RNA: ab initio investigations of structures and energies

Decoding the enigma of RNA-protein recognition: quantum chemical insights into arginine fork motifs

Arginine (Arg) forks are noncovalent recognition motifs wherein an Arg interacts with the phosphates and guanine nucleobases of RNA, providing extraordinary specific RNA:protein recognition. In this work, we carried out an in-depth DFT based quantum mechanical investigation on all known classes of Arg forks to estimate their intrinsic structural stabilities and interaction energies. The optimized structures closely mimic the structural characteristics of Arg forks and this close match between experimental and optimized geometries suggests that Arg forks are intrinsically stable and do not require additional support from other RNA or protein components. Both hydrogen-bonding and cation-π interactions are important for the intrinsic stability of Arg forks, providing an average interaction energy of -36.7 kcal mol-1. Furthermore, we found a direct correlation between Arg forks' interaction energies and the number of phosphates involved, which is more delicately modulated by other factors, like the types of hydrogen bonds and cation-π interactions that constitute the Arg fork. Additionally, we observed a positive correlation between the average interaction energies of Arg forks and the frequency of their occurrence in available crystal structures. At the broader level, this work establishes the groundwork for more precise modeling and understanding of RNA-protein interfaces, which could have potential implications in advancing the knowledge of biomolecular recognition patterns.

The assessment of molecular dynamics results of three-dimensional RNA aptamer structure prediction

Aptamers are single-stranded DNA or RNA that bind to specific targets such as proteins, thus having similar characteristics to antibodies. It can be synthesized at a lower cost, with no batch-to-batch variations, and is easier to modify chemically than antibodies, thus potentially being used as therapeutic and biosensing agents. The current method for RNA aptamer identification in vitro uses the SELEX method, which is considered inefficient due to its complex process. Computational models of aptamers have been used to predict and study the molecular interaction of modified aptamers to improve affinity. In this study, we generated three-dimensional models of five RNA aptamers from their sequence using mFold, RNAComposer web server, and molecular dynamics simulation. The model structures were then evaluated and compared with the experimentally determined structures. This study showed that the combination of mFold, RNAComposer, and molecular dynamics simulation could generate 14-16, 28, or 29 nucleotides length of 3D RNA aptamer with similar geometry and topology to the experimentally determined structures. The non-canonical basepair structure of the aptamer loop was formed through the MD simulation, which also improved the three-dimensional RNA aptamers model. Clustering analysis was recommended to choose the more representative model.

Structural Properties of Hachimoji Nucleic Acids and Their Building Blocks: Comparison of Genetic Systems with Four, Six and Eight Alphabets

Expansion of the genetic alphabet is an ambitious goal. A recent breakthrough has led to the eight-base (hachimoji) genetics having canonical and unnatural bases. However, very little is known on the molecular-level features that facilitate the candidature of unnatural bases as genetic alphabets. Here we amalgamated DFT calculations and MD simulations to analyse the properties of the constituents of hachimoji DNA and RNA. DFT reveals the dominant syn conformation for isolated unnatural deoxyribonucleosides and at the 5'-end of oligonucleotides, although an anti/syn mixture is predicted at the nonterminal and 3'-terminal positions. However, isolated ribonucleotides prefer an anti/syn mixture, but mostly prefer anti conformation at the nonterminal positions. Further, the canonical base pairing combinations reveals significant strength, which may facilitate replication of hachimoji DNA. We also identify noncanonical base pairs that can better tolerate the substitution of unnatural pairs in RNA. Stacking strengths of 51 dimers reveals higher average stacking stabilization of dimers of hachimoji bases than canonical bases, which provides clues for choosing energetically stable sequences. A total of 14.4 μs MD simulations reveal the influence of solvent on the properties of hachimoji oligonucleotides and point to the likely fidelity of replication of hachimoji DNA. Our results pinpoint the features that explain the experimentally observed stability of hachimoji DNA.

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.

RNA folding using quantum computers

The 3-dimensional fold of an RNA molecule is largely determined by patterns of intramolecular hydrogen bonds between bases. Predicting the base pairing network from the sequence, also referred to as RNA secondary structure prediction or RNA folding, is a nondeterministic polynomial-time (NP)-complete computational problem. The structure of the molecule is strongly predictive of its functions and biochemical properties, and therefore the ability to accurately predict the structure is a crucial tool for biochemists. Many methods have been proposed to efficiently sample possible secondary structure patterns. Classic approaches employ dynamic programming, and recent studies have explored approaches inspired by evolutionary and machine learning algorithms. This work demonstrates leveraging quantum computing hardware to predict the secondary structure of RNA. A Hamiltonian written in the form of a Binary Quadratic Model (BQM) is derived to drive the system toward maximizing the number of consecutive base pairs while jointly maximizing the average length of the stems. A Quantum Annealer (QA) is compared to a Replica Exchange Monte Carlo (REMC) algorithm programmed with the same objective function, with the QA being shown to be highly competitive at rapidly identifying low energy solutions. The method proposed in this study was compared to three algorithms from literature and, despite its simplicity, was found to be competitive on a test set containing known structures with pseudoknots.

The recent FDA approval of mRNA-based vaccines has increased public interest in synthetically designed RNA molecules. RNA molecules fold into complex secondary structures which determine their molecular properties and in part their efficacy. Determining the folded structure of an RNA molecule is a computationally challenging task with exponential scaling that is intractable to solve exactly, and therefore approximate methods are used. Quantum computing technology offers a new approach to finding approximate solutions to problems with exponential scaling. We formulate a simplistic, yet effective, model of RNA folding that can easily be mapped to quantum computers and we show that currently available quantum computing hardware is competitive with classical methods.

PolyPurine Reverse Hoogsteen Hairpins Work as RNA Species for Gene Silencing

PolyPurine Reverse Hoogsteen Hairpins (PPRHs) are gene-silencing DNA-oligonucleotides developed in our laboratory that are formed by two antiparallel polypurine mirror repeat domains bound intramolecularly by Hoogsteen bonds. The aim of this work was to explore the feasibility of using viral vectors to deliver PPRHs as a gene therapy tool. After treatment with synthetic RNA, plasmid transfection, or viral infection targeting the survivin gene, viability was determined by the MTT assay, mRNA was determined by RT-qPCR, and protein levels were determined by Western blot. We showed that the RNA-PPRH induced a decrease in cell viability in a dose-dependent manner and an increase in apoptosis in PC-3 and HeLa cells. Both synthetic RNA-PPRH and RNA-PPRH intracellularly generated upon the transfection of a plasmid vector were able to reduce survivin mRNA and protein levels in PC-3 cells. An adenovirus type-5 vector encoding the PPRH against survivin was also able to decrease survivin mRNA and protein levels, leading to a reduction in HeLa cell viability. In this work, we demonstrated that PPRHs can also work as RNA species, either chemically synthesized, transcribed from a plasmid construct, or transcribed from viral vectors. Therefore, all these results are the proof of principle that viral vectors could be considered as a delivery system for PPRHs.

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.

Molecular dynamics study of ways of RNA base-pair formation

Base pairing is a fundamental phenomenon in RNA structure and function. However, although there have been considerable recent advances, some important aspects of base-pair formation are still unknown, including the ways of base-pair formation and detailed roles of metal ions. Here we show that base pairs can form through four different ways: stabilizing, bridging, rotating, and shifting. Among them the stabilizing and bridging ways involve direct binding of metal ions while the rotating and shifting ways do not in most cases. Furthermore, we find that the formations of base pairs in different positions of the hairpin stem may adopt different ways.

A quantum chemical view of the interaction of RNA nucleobases and base pairs with the side chains of polar amino acids

Hydrogen bonding between amino acids and nucleobases is important for RNA-protein recognition. As a first step toward understanding the physicochemical features of these contacts, the present work employs density functional theory calculations to critically analyze the intrinsic structures and strength of all theoretically possible model hydrogen-bonded complexes involving RNA nucleobase edges and polar amino acid side chains. Our geometry optimizations uncover a number of unique complexes that involve variable hydrogen-bonding characteristics, including conventional donor-acceptor interactions, bifurcated interactions and single hydrogen-bonded contacts. Further, significant strength of these complexes in the gas phase (-27 kJ mol^-1 to -226 kJ mol^-1) and solvent phase (-19 kJ mol^-1 to -78 kJ mol^-1) points toward the ability of associated contacts to provide stability to RNA-protein complexes. More importantly, for the first time, our study uncovers the features of complexes involving protonated nucleobases, as well as those involving the weakly polar cysteine side chain, and thereby highlights their potential importance in biological processes that involve RNA-protein interactions. Additional analysis on select base pair-amino acid complexes uncovers the ability of amino acid side chain to simultaneously interact with both nucleobases of the base pair, and highlights the greater strength of such interactions compared to base-amino acid interactions. Overall, our analysis provides a basic physicochemical framework for understanding the molecular basis of nucleic acid-protein interactions. Further, our quantum chemical data can be used to design better algorithms for automated search of these contacts at the RNA-protein interface.Communicated by Ramaswamy H. Sarma.

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.

Can Cyanuric Acid and 2,4,6-Triaminopyrimidine Containing Ribonucleosides be Components of Prebiotic RNA? Insights from QM Calculations and MD Simulations

As a step toward assessing their fitness as pre-RNA nucleobases, we employ DFT and MD simulations to analyze the noncovalent interactions of cyanuric acid (CA) and 2,4,6-triaminopyrimidine (TAP), and the structural properties of the associated ribonucleosides (rNs) and oligonucleotides. Our calculations reveal that the TAP : CA pair has a comparable hydrogen-bond strength to the canonical A : U pair. This strengthens the candidature of CA and TAP as prebiotic nucleobases. Further, the stacking between two canonical nucleobases is stronger than those between TAP or CA and a canonical base, as well as those between two TAP and/or CA, which indicates that enhanced stacking may have served as a driving force for the evolution from prebiotic to canonical nucleobases. Similarities in the DFT-derived anti/syn rotational barriers and MD-derived (anti) glycosidic conformation of the CA and TAP rNs and canonical rNs further substantiate their candidature as pre-RNA components. Greater deglycosylation barriers (as obtained by DFT calculations) for TAP rNs compared to canonical rNs suggest TAP rNs indicate higher resistance to environmental factors, while lower barriers indicate that CA rNs were likely more suitable for less-challenging locations. Finally, the tight packing in narrow CA:TAP-containing helices suggests that the prebiotic polymers were shielded from water, which would aid their evolution into self-replicating systems. Our calculations thus support proposals that CA and TAP can act as nucleobases of pre-RNA.

Going beyond base-pairs: topology-based characterization of base-multiplets in RNA

Identification and characterization of base-multiplets, which are essentially mediated by base-pairing interactions, can provide insights into the diversity in the structure and dynamics of complex functional RNAs, and thus facilitate hypothesis driven biological research. The necessary nomenclature scheme, an extension of the geometric classification scheme for base-pairs by Leontis and Westhof, is however available only for base-triplets. In the absence of information on topology, this scheme is not applicable to quartets and higher order multiplets. Here we propose a topology-based classification scheme which, in conjunction with a graph-based algorithm, can be used for the automated identification and characterization of higher order base-multiplets in RNA structures. Here, the RNA structure is represented as a graph, where nodes represent nucleotides and edges represent base-pairing connectivity. Sets of connected components (of n nodes) within these graphs constitute subgraphs representing multiplets of "n" nucleotides. The different topological variants of the RNA multiplets thus correspond to different nonisomorphic forms of these subgraphs. To annotate RNA base-multiplets unambiguously, we propose a set of topology-based nomenclature rules for quartets, which are extendable to higher multiplets. We also demonstrate the utility of our approach toward the identification and annotation of higher order RNA multiplets, by investigating the occurrence contexts of selected examples in order to gain insights regarding their probable functional roles.

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.

The H/ACA complex disrupts triplex in hTR precursor to permit processing by RRP6 and PARN

Human telomerase RNA (hTR) is transcribed as a precursor that is then posttranscriptionally modified and processed. A fraction of the transcripts is oligoadenylated by TRAMP and either processed into the mature hTR or degraded by the exosome. Here, we characterize the processing of 3′ extended forms of varying length by PARN and RRP6. We show that tertiary RNA interactions unique to the longer transcripts favor RNA degradation, whereas H/ACA RNP assembly stimulates productive processing. Interestingly, the H/ACA complex actively promotes processing in addition to protecting the mature 3′ end. Processing occurs in two steps with longer forms first being trimmed by RRP6 and shorter forms then being processed by PARN. These results reveal how RNA structure and RNP assembly affect the kinetics of processing and degradation and ultimately determine the amount of functional telomerase produced in cells.

Telomerase RNA (hTR) is transcribed as a 3′-extended precursor. Here the authors examine the processing of hTR precursors of various lengths and show that processing occurs in distinct steps involving different nucleases PARN and RRP6.

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.

B-DNA model systems in non-terran bio-solvents: implications for structure, stability and replication

We have computationally analyzed a comprehensive series of Watson–Crick and mismatched B-DNA base pairs, in the gas phase and in several solvents, including toluene, chloroform, ammonia, methanol and water, using dispersion-corrected density functional theory and implicit solvation.

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.

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.

RNABP COGEST: a resource for investigating functional RNAs

Structural bioinformatics of RNA has evolved mainly in response to the rapidly accumulating evidence that non-(protein)-coding RNAs (ncRNAs) play critical roles in gene regulation and development. The structures and functions of most ncRNAs are however still unknown. Most of the available RNA structural databases rely heavily on known 3D structures, and contextually correlate base pairing geometry with actual 3D RNA structures. None of the databases provide any direct information about stabilization energies. However, the intrinsic interaction energies of constituent base pairs can provide significant insights into their roles in the overall dynamics of RNA motifs and structures. Quantum mechanical (QM) computations provide the only approach toward their accurate quantification and characterization. ‘RNA Base Pair Count, Geometry and Stability’ (http://bioinf.iiit.ac.in/RNABPCOGEST) brings together information, extracted from literature data, regarding occurrence frequency, experimental and quantum chemically optimized geometries, and computed interaction energies, for non-canonical base pairs observed in a non-redundant dataset of functional RNA structures. The database is designed to enable the QM community, on the one hand, to identify appropriate biologically relevant model systems and also enable the biology community to easily sift through diverse computational results to gain theoretical insights which could promote hypothesis driven biological research.

Database URL:http://bioinf.iiit.ac.in/RNABPCOGEST

Why does substitution of thymine by 6-ethynylpyridone increase the thermostability of DNA double helices?

Efficiency of 6-ethynylpyridone (E), a potential thymine (T) analogue, which forms high-fidelity base pairs with adenine (A) and gives rise to stabler DNA duplexes, with stability comparable to those containing canonical cytosine(C):guanine(G) base pairs, has been reported recently. Estimates of the interaction energies, involving geometry optimization at the DFT level (including middle range dispersion interactions) followed by single point energy calculation at MP2 level, in excellent correlation with the experimentally observed trends, show that E binds more strongly and more discriminately with A than T does. Detailed analysis reveals that the increase in base-base interaction arises out of conjugation of acetylenic π electrons with the ring π system of E, which results in not only an extra stabilizing C-H···π interaction in the EA pair, but also a strengthening of the conventional hydrogen bonds. However, the computed base-base interaction energy for the EA pair was found to be much less than that of the canonical CG pair, implying that the difference in the TA versus EA base pairing interaction alone cannot explain the large experimentally observed increase in the thermostability of DNA duplexes, where a TA pair is replaced with an EA pair. Our computations show that the conjugation of acetylenic π electrons with the ring π system also possibly plays a role in increasing the stacking potential of the EA pair, which in turn can explain its marked influence in the enhancement of duplex stability.

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

The G:C reverse Watson-Crick (W:W trans) base pair, also known as Levitt base pair in the context of tRNAs, is a structurally and functionally important base pair that contributes to tertiary interactions joining distant domains in functional RNA molecules and also participates in metabolite binding in riboswitches. We previously indicated that the isolated G:C W:W trans base pair is a rather unstable geometry, and that dicationic metal binding to the Guanine base or posttranscriptional modification of the Guanine can increase its stability. Herein, we extend our survey and report on other H-bonding interactions that can increase the stability of this base pair. To this aim, we performed a bioinformatics search of the PDB to locate all the occurencies of G:C trans base pairs. Interestingly, 66% of the G:C trans base pairs in the PDB are engaged in additional H-bonding interactions with other bases, the RNA backbone or structured water molecules. High level quantum mechanical calculations on a data set of representative crystal structures were performed to shed light on the structural stability and energetics of the various crystallographic motifs. This analysis was extended to the binding of the preQ1 metabolite to a preQ1-II riboswitch.

Configuration of wobble base pairs having pyrimidines as anticodon wobble bases: significance for codon degeneracy

Degeneracy of the genetic code was attributed by Crick to imprecise hydrogen-bonded base-pairing at the wobble position during codon-anticodon pairing. The Crick wobble rules define but do not explain the RNA base pair combinations allowed at this position. We select six pyrimidine bases functioning as anticodon wobble bases (AWBs) to study their H-bonded pairing properties with the four major RNA bases using density functional theory at the B3LYP/6-31G(d,p) level. This is done to assess the extent to which the configuration of a solitary RNA wobble base pair may in itself determine specificity and degeneracy of the genetic code by allowing or disallowing the given base pair during codon-anticodon pairing. Calculated values of select configuration markers for the base pairs screen well between allowed and disallowed base pairs for most cases examined here, where the base pair width emerges as an important factor. A few allowed wobble pairs invoke the involvement of RNA nucleoside conformation, as well as involvement of the exocyclic substituent in H-bonding. This study, however, cannot explain the disallowed status of the Ura⋯Gua wobble pair on the basis of configuration alone. Explanation of the allowed status of the V⋯Ura pair requires further study on the mediatory role of water molecules. Apart from these two cases, these computational results are sufficient, on the basis of base pair configuration alone, to account for the specificity and degeneracy of the genetic code for all known cases of codon-anticodon pairing which involve the pyrimidine AWBs studied here.

Role of wobble base pair geometry for codon degeneracy: purine-type bases at the anticodon wobble position

Codon degeneracy is a key feature of the genetic code, explained by Crick (J Mol Biol 19:548-555, 1966) in terms of imprecision of base pairing at the codon third position (the wobble position) of the codon-anticodon duplex. The Crick wobble rules define, but do not explain, which base pairs are allowed/disallowed at the wobble position of this duplex. This work examines whether the H-bonded configurations of solitary RNA base pairs can in themselves help decide which base pairs are allowed at the wobble position during codon-anticodon pairing. Taking the purine-type bases guanine, hypoxanthine, queuine and adenine as anticodon wobble bases, H-bonded pairing energies and optimized configurations of numerous RNA base pairs are calculated in gas and modeled aqueous phase at the B3LYP/6-31 G(d,p) level. Calculated descriptors of alignment of these solitary base pairs are able to screen between allowed and disallowed base pairs for all cases studied here, except two cases which invoke base-sugar interactions in the codon wobble nucleoside. The exclusion of adenine from the anticodon wobble position cannot be explained on the basis of pairing facility or base pair geometry. These DFT results thus account for the specificity and degeneracy of the genetic code for all cases involving guanine, hypoxanthine and queuine as anticodon wobble bases.

A-minor tertiary interactions in RNA kink-turns. Molecular dynamics and quantum chemical analysis

The RNA kink-turn is an important recurrent RNA motif, an internal loop with characteristic consensus sequence forming highly conserved three-dimensional structure. Functional arrangement of RNA kink-turns shows a sharp bend in the phosphodiester backbone. Among other signature interactions, kink-turns form A-minor interaction between their two stems. Most kink-turns possess extended A-minor I (A-I) interaction where adenine of the second A•G base pair of the NC-stem interacts with the first canonical pair of the C-stem (i.e., the receptor pair) via trans-sugar-edge/sugar-edge (tSS) and cis-sugar-edge/sugar-edge (cSS) interactions. The remaining kink-turns have less compact A-minor 0 (A-0) interaction with just one tSS contact. We show that kink-turns with A-I in ribosomal X-ray structures keep G═C receptor base pair during evolution while the inverted pair (C═G) is not realized. In contrast, kink-turns with A-0 in the observed structures alternate G═C and C═G base pairs in sequences. We carried out an extended set (~5 μs) of explicit-solvent molecular dynamics simulations of kink-turns to rationalize this structural/evolutionary pattern. The simulations were done using a net-neutral Na(+) cation atmosphere (with ~0.25 M cation concentration) supplemented by simulations with either excess salt KCl atmosphere or inclusion of Mg(2+). The results do not seem to depend on the treatment of ions. The simulations started with X-ray structures of several kink-turns while we tested the response of the simulated system to base substitutions, modest structural perturbations and constraints. The trends seen in the simulations reveal that the A-I/G═C arrangement is preferred over all three other structures. The A-I/C═G triple appears structurally entirely unstable, consistent with the covariation patterns seen during the evolution. The A-0 arrangements tend to shift toward the A-I pattern in simulations, which suggests that formation of the A-0 interaction is likely supported by the surrounding protein and RNA molecules. A-0 may also be stabilized by additional kink-turn nucleotides not belonging to the kink-turn consensus, as shown for the kink-turn from ribosomal Helix 15. Quantum-chemical calculations on all four A-minor triples suggest that there is a different balance of electrostatic and dispersion stabilization in the A-I/G═C and A-I/C═G triples, which may explain different behavior of these otherwise isosteric triples in the context of kink-turns.

Noncanonical hydrogen bonding in nucleic acids. Benchmark evaluation of key base-phosphate interactions in folded RNA molecules using quantum-chemical calculations and molecular dynamics simulations

RNA molecules are stabilized by a wide range of noncanonical interactions that are not present in DNA. Among them, the recently classified base-phosphate (BPh) interactions belong to the most important ones. Twelve percent of nucleotides in the ribosomal crystal structures are involved in BPh interactions. BPh interactions are highly conserved and provide major constraints on RNA sequence evolution. Here we provide assessment of the energetics of BPh interactions using MP2 computations extrapolated to the complete basis set of atomic orbitals and corrected for higher-order electron correlation effects. The reference computations are compared with DFT-D and DFT-D3 approaches, the SAPT method, and the molecular mechanics force field. The computations, besides providing the basic benchmark for the BPh interactions, allow some refinements of the original classification, including identification of some potential doubly bonded BPh patterns. The reference computations are followed by analysis of some larger RNA fragments that consider the context of the BPh interactions. The computations demonstrate the complexity of interaction patterns utilizing the BPh interactions in real RNA structures. The BPh interactions are often involved in intricate interaction networks. We studied BPh interactions of protonated adenine that can contribute to catalysis of hairpin ribozyme, the key BPh interaction in the S-turn motif of the sarcin-ricin loop, which may predetermine the S-turn topology and complex BPh patterns from the glmS riboswitch. Finally, the structural stability of BPh interactions in explicit solvent molecular dynamics simulations is assessed. The simulations well preserve key BPh interactions and allow dissection of structurally/functionally important water-meditated BPh bridges, which could not be considered in earlier bioinformatics classification of BPh interactions.

Protonation of base pairs in RNA: context analysis and quantum chemical investigations of their geometries and stabilities

Base pairs involving protonated nucleobases play important roles in mediating global macromolecular conformational changes and in facilitation of catalysis in a variety of functional RNA molecules. Here we present our attempts at understanding the role of such base pairs by detecting possible protonated base pairs in the available RNA crystal structures using BPFind software, in their specific structural contexts, and by the characterization of their geometries, interaction energies, and stabilities using advanced quantum chemical computations. We report occurrences of 18 distinct protonated base pair combinations from a representative data set of RNA crystal structures and propose a theoretical model for one putative base pair combination. Optimization of base pair geometries was carried out at the B3LYP/cc-pVTZ level, and the BSSE corrected interaction energies were calculated at the MP2/aug-cc-pVDZ level of theory. The geometries for each of the base pairs were characterized in terms of H-bonding patterns observed, rmsd values observed on optimization, and base pair geometrical parameters. In addition, the intermolecular interaction in these complexes was also analyzed using Morokuma energy decomposition. The gas phase interaction energies of the base pairs range from -24 to -49 kcal/mol and reveal the dominance of Hartree-Fock component of interaction energy constituting 73% to 98% of the total interaction energy values. On the basis of our combined bioinformatics and quantum chemical analysis of different protonated base pairs, we suggest resolution of structural ambiguities and correlate their geometric and energetic features with their structural and functional roles. In addition, we also examine the suitability of specific base pairs as key elements in molecular switches and as nucleators for higher order structures such as base triplets and quartets.

Quantum chemical studies of nucleic acids: can we construct a bridge to the RNA structural biology and bioinformatics communities?

In this feature article, we provide a side-by-side introduction for two research fields: quantum chemical calculations of molecular interaction in nucleic acids and RNA structural bioinformatics. Our main aim is to demonstrate that these research areas, while largely separated in contemporary literature, have substantial potential to complement each other that could significantly contribute to our understanding of the exciting world of nucleic acids. We identify research questions amenable to the combined application of modern ab initio methods and bioinformatics analysis of experimental structures while also assessing the limitations of these approaches. The ultimate aim is to attain valuable physicochemical insights regarding the nature of the fundamental molecular interactions and how they shape RNA structures, dynamics, function, and evolution.