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.

Solvent effects on the NMR shieldings of stacked DNA base pairs

Stacking effects are among the most important effects in DNA. We have recently studied their influence in fragments of DNA through the analysis of NMR magnetic shieldings, firstly in vacuo. As a continuation of this line of research we show here the influence of solvent effects on the shieldings through the application of both explicit and implicit models. We found that the explicit solvent model is more appropriate for consideration due to the results matching better in general with experiments, as well as providing clear knowledge of the electronic origin of the value of the shieldings. Our study is grounded on a recently developed theoretical model of our own, by which we are able to learn about the magnetic effects of given fragments of DNA molecules on selected base pairs. We use the shieldings of the atoms of a central base pair (guanine-cytosine) of a selected fragment of DNA molecules as descriptors of physical effects, like π-stacking and solvent effects. They can be taken separately and altogether. The effect of π-stacking is introduced through the addition of some pairs above and below of the central base pair, and now, the solvent effect is considered including a network of water molecules that consist of two solvation layers, which were fixed in the calculations performed in all fragments. We show that the solvent effects enhance the stacking effects on the magnetic shieldings of atoms that belong to the external N-H bonds. The net effect is of deshielding on both atoms. There is also a deshielding effect on the carbon atoms that belong to CO bonds, for which the oxygen atom has an explicit hydrogen bond (HB) with a solvent water molecule. Solvent effects are found to be no higher than a few percent of the total value of the shieldings (between 1% and 5%) for most atoms, although there are few for which such an effect can be higher. There is one nitrogen atom, the acceptor of the HB between guanine and cytosine, that is more highly shielded (around 15 ppm or 10%) when the explicit solvent is considered. In a similar manner, the most external nitrogen atom of cytosine and the hydrogen atom that is bonded to it are highly deshielded (around 10 ppm for nitrogen and around 3 ppm for hydrogen).

Role of Stacking Interactions in the Stability of Primitive Genetics: A Quantum Chemical View

The origin of genetic material on earth is an age-old, entangled mystery that lacks a unanimous explanation. Recent studies have suggested that noncanonical bases such as barbituric acid (BA), melamine (MM), cyanuric acid (CA), and 2,4,6-triaminopyrimidine (TAP) may have undergone molecular selection within the "prebiotic soup" to spontaneously form supramolecular assemblies, which then covalently assembled into an RNA-like polymer (preRNA). However, information on the role of intrinsic interactions of these candidate heterocycles in their molecular selection as the components of preRNA, and the subsequent transition from preRNA to RNA, is currently missing in the literature. To fill this gap in our knowledge on the origin and evolution of primitive genetics, the present work employs density functional theory (B3LYP-D3) to evaluate and compare the stacking propensities of dimers containing prebiotic noncanonical (BA, MM, CA, and TAP) and/or canonical RNA bases (A, C, G, and U). Our detailed analysis of the variation in stacking strength with respect to four characteristic geometrical parameters between the monomers [i.e., the vertical distance, the angle of rotation, and (two) displacements in the x and y directions] reveals that stacking between nonidentical bases is preferred over identical bases for both prebiotic-prebiotic and canonical-canonical dimers. This not only underscores the similarity between the fundamental chemical properties of preRNA and RNA constituents but also supports the likelihood of the evolution of modern (RNA) genetics from primitive (preRNA) genetics. Furthermore, greater average stacking stabilization of canonical dimers than that of dimers containing one canonical and one preRNA nucleobase (by ∼5 kJ mol^-1) or dimers solely containing preRNA nucleobases (by ∼12 kJ mol^-1) indicates that enhanced stacking is an important factor that may have spurred the evolution of preRNA to an intermediate informational polymer to RNA. More importantly, our study identifies the central roles of CA, BA, and TAP in stacking stabilization within the preRNA and of BA in stacking interactions within the intermediate polymers and suggests that these heterocycles may have played distinct roles in various stages during the evolution from preRNA to RNA. Overall, our results highlight the significance of stacking interactions in the selection of nucleobase components of preRNA.

Fundamental photophysics of isomorphic and expanded fluorescent nucleoside analogues

Fluorescent nucleoside analogues (FNAs) are structurally diverse mimics of the natural essentially non-fluorescent nucleosides which have found numerous applications in probing the structure and dynamics of nucleic acids as well as their interactions with various biomolecules. In order to minimize disturbance in the labelled nucleic acid sequences, the FNA chromophoric groups should resemble the natural nucleobases in size and hydrogen-bonding patterns. Isomorphic and expanded FNAs are the two groups that best meet the criteria of non-perturbing fluorescent labels for DNA and RNA. Significant progress has been made over the past decades in understanding the fundamental photophysics that governs the spectroscopic and environmentally sensitive properties of these FNAs. Herein, we review recent advances in the spectroscopic and computational studies of selected isomorphic and expanded FNAs. We also show how this information can be used as a rational basis to design new FNAs, select appropriate sequences for optimal spectroscopic response and interpret fluorescence data in FNA applications.

Computational insight into the halogen bonded self-assembly of hexa-coordinated metalloporphyrins

We demonstrate herein a computational study probing the influence of metalloporphyrins on intermolecular halogen bonding (XB) during supramolecular self-assembly. The results demonstrate that porphyrin aromatic rings can activate or deactivate halogen bonding interactions, especially those on axial ligands, and further influence the preference type of halogen···halogen bonding during the supramolecular self-assembly. Calculations show that the halogen atom present at the equatorial position has a higher sigma hole potential (V_S,max) than that at the axial position. The computational analysis and our observations from the X-ray structure analysis are in good agreement. From structural analysis it is clear that equatorial halogen atoms prefer to participate in Type-II XB interactions whereas the axial halogen atoms either participate in Type-I XB interaction or reluctant to participate in XB interactions due to the decrease of their sigma hole potential. Thus, we demonstrate, herein, for the first time a computational study probing the direct influence of the porphyrin's ring current on the sigma hole potential (V_S,max) of the halogen atoms and subsequently the effects of the supramolecular self-assembly.

Structural Patterns and Stabilities of Hydrogen-Bonded Pairs Involving Ribonucleotide Bases and Arginine, Glutamic Acid, or Glutamine Residues of Proteins from Quantum Mechanical Calculations

Ribonucleotide:protein interactions play crucial roles in a number of biological processes. Unlike the RNA:protein interface where van der Waals contacts are prevalent, the recognition of a single ribonucleotide such as ATP by a protein occurs predominantly through hydrogen-bonding interactions. As a first step toward understanding the role of hydrogen bonding in ribonucleotide:protein recognition, the present work employs density functional theory to provide a detailed quantum-mechanical analysis of the structural and energetic characteristics of 18 unique hydrogen-bonded pairs involving the nucleobase/nucleoside moiety of four canonical ribonucleotides and the side chains of three polar amino-acid residues (arginine, glutamine, and glutamic acid) of proteins. In addition, we model five new pairs that are till now not observed in crystallographically identified ribonucleotide:protein complexes but may be identified in complexes crystallized in the future. We critically examine the characteristics of each pair in its ribonucleotide:protein crystal structure occurrence and (gas phase and water phase) optimized intrinsic structure. We further evaluated the interaction energy of each pair and characterized the associated hydrogen bonds using a number of quantum mechanics-based relationships including natural bond orbital analysis, quantum theory atoms in molecules analysis, Iogansen relationships, Nikolaienko-Bulavin-Hovorun relationships, and noncovalent interaction-reduced density gradient analysis. Our analyses reveal rich variability in hydrogen bonds in the crystallographic as well as intrinsic structure of each pair, which includes conventional O/N-H···N/O and C-H···O hydrogen bonds as well as donor/acceptor-bifurcated hydrogen bonds. Further, we identify five combinations of nucleobase and amino acid moieties; each of which exhibits at least two alternate (i.e., multimodal) structures that interact through the same nucleobase edge. In fact, one such pair exhibits four multimodal structures; one of which possesses unconventional "amino-acceptor" hydrogen bonding with comparable (-9.4 kcal mol^-1) strength to the corresponding conventional (i.e., amino:donor) structure (-9.2 kcal mol^-1). This points to the importance of amino-acceptor hydrogen bonds in RNA:protein interactions and suggests that such interactions must be considered in the future while studying the dynamics in the context of molecular recognition. Overall, our study provides preliminary insights into the intrinsic features of ribonucleotide:amino acid interactions, which may help frame a clearer picture of the molecular basis of RNA:protein recognition and further appreciate the role of such contacts in biology.

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.

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.

How do hydrophobic nucleobases differ from natural DNA nucleobases? Comparison of structural features and duplex properties from QM calculations and MD simulations

Computational (DFT and MD simulation) methods are employed to systematically characterize the structural and energetic properties of five hydrophobic nucleobases (FEMO, MMO2, NaM, 5SICS and TPT3) that constitute four unnatural base pairs (FEMO:5SICS, MMO2:5SICS, NaM:5SICS and TPT3:NaM). These hydrophobic bases have been recently shown to be replicated when present between natural bases in DNA duplexes, with the highest replication fidelity and efficiency occuring for the TPT3:NaM pair. Our QM calculations suggest that the preferred (anti) glycosidic orientations of nucleosides containing hydrophobic bases are similar to the natural DNA nucleosides despite differences in their chemical structures. However, due to the inability to form interbase hydrogen bonds, hydrophobic base pairs intrinsically prefer nonplanar, distorted geometries, many of which are stabilized through π-π stacking interactions. Furthermore, the intrinsic stacking potential between a hydrophobic and a natural base is similar to that between two natural bases, indicating that the strength of stacking interactions in DNA duplexes containing hydrophobic bases is likely comparable to natural DNA. However, in contrast to the isolated base-pair geometries, our MD simulations suggest that the hydrophobic base pairs adopt variable geometries within DNA, which range from stacked (5SICS:FEMO) to nearly planar (5SICS:NaM and SICS:MMO2) to planar (TPT3:NaM). As a result, the duplex structural features at the site of modification depend on the identity of the hydrophobic base pair, where the TPT3:NaM pair causes the least structural changes compared to natural DNA. Overall, the structural insight obtained from our calculations on DNA containing hydrophobic base pairs explains the experimentally-observed higher fidelity and efficiency during replication of TPT3:NaM compared to other hydrophobic nucleobase pairs. By providing valuable structural information that explains the intrinsic and duplex properties of this class of unnatural nucleobases, the present work may aid the future design of improved hydrophobic analogues.

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.

Structural and electronic properties of barbituric acid and melamine-containing ribonucleosides as plausible components of prebiotic RNA: implications for prebiotic self-assembly

DFT calculations are employed to assess the fitness of barbituric acid and melamine to act as components of prebiotic RNA.

Nucleotides containing variously modified sugars: energetics, structure, and mechanical properties

The influence of various sugar residue modifications on intrinsic energetic, conformational, and mechanical properties of 2'-deoxyribonucleotide-5'-monophosphates (dNs) was comprehensively investigated using modern quantum chemical approaches. In total, fourteen sugar modifications, including double bonds and heteroatoms (S and N) inside the sugar ring, as well as fluorination in various positions, were analyzed. Among hundreds of possible conformational states of dNs, only two - AI and BI, corresponding to the most biologically significant forms of a double-helical DNA, were considered for each dN. It was established that the most of the studied modifications tend to strongly stabilize either AI or BI conformation of dNs both in the gas phase and in aqueous solution (modelled by implicit solvent models). Therefore, some of these modifications can be used as a tool for reducing structural polymorphism of nucleic acids in solution as well as for designing oligonucleotides with specific structural features. The evaluation of relaxed force constants (RFC) for glycosidic bonds suggests that the majority of the studied modifications of the sugar residue yield increased strengths of glycosidic bonds in dNs, and can therefore be used for designing modified nucleic acids with an increased resistance to abasic lesions. The most significant reinforcement of the glycosidic bond occurs in dNs containing the CF2 group instead of the O4' oxygen and the fluorine atom at the 2'-α-position. The calculation of the RFC and vibrational root-mean-square (VRMS) deviations for conformational degrees of freedom revealed a strong dependence between mechanical properties of dNs and their energetic characteristics. In particular, electronic energies of AI and BI conformers of dNs calculated in vacuo are closely connected with the values of relaxed force constants (RFC) for the δ angle: the higher RFC(δ) values correspond to more energetically favorable conformers.

New size-expanded RNA nucleobase analogs: a detailed theoretical study

Fluorescent nucleobase analogs have attracted much attention in recent years due to their potential applications in nucleic acids research. In this work, four new size-expanded RNA base analogs were computationally designed and their structural, electronic, and optical properties are investigated by means of DFT calculations. The results indicate that these analogs can form stable Watson-Crick base pairs with natural counterparts and they have smaller ionization potentials and HOMO-LUMO gaps than natural ones. Particularly, the electronic absorption spectra and fluorescent emission spectra are calculated. The calculated excitation maxima are greatly red-shifted compared with their parental and natural bases, allowing them to be selectively excited. In gas phase, the fluorescence from them would be expected to occur around 526, 489, 510, and 462 nm, respectively. The influences of water solution and base pairing on the relevant absorption spectra of these base analogs are also examined.