Prebiotic N-(2-Aminoethyl)-Glycine (AEG)-Assisted Synthesis of Proto-RNA?

Quantum mechanical calculations are used to explore the thermodynamics of possible prebiotic synthesis of the building blocks of nucleic acids. Different combinations of D-ribofuranose (Ribf) and N-(2-aminoethyl)-glycine (AEG) (trifunctional connectors (TCs)); the nature of the Ribf, its anomeric form, and its ring puckering (conformation); and the nature of the nucleobases (recognition units (RUs)) are considered. The combinatorial explosion of possible nucleosides has been drastically reduced on physicochemical grounds followed by a detailed thermodynamic evaluation of alternative synthetic pathways. The synthesis of nucleosides containing N-(2-aminoethyl)-glycine (AEG) is predicted to be thermodynamically favored suggesting a possible role of AEG as a component of an ancestral proto-RNA that may have preceded today's nucleic acids. A new pathway for the building of free nucleotides (exemplified by 5'-uridine monophosphate (UMP)) and of AEG dipeptides is proposed. This new pathway leads to a spontaneous formation of free UMP assisted by an AEG nucleoside in an aqueous environment. This appears to be a workaround to the "water problem" that prohibits the synthesis of nucleotides in water.

Accelerated Screening of Alternative DNA Base-Organic Molecule-Base Architectures via Integrated Theory and Experiment

Organic molecule-mediated noncanonical DNA self-assembly expands the standard DNA base-pairing alphabets. However, only a very limited number of small molecules have been recognized as mediators because of the tedious and complicated experiments like crystallization and microscopy imaging. Here we present an integrative screening protocol incorporating molecular dynamics (MD) for fast theoretical simulation and native polyacrylamide gel electrophoresis for convenient experimental validation. Melamine, the molecule that was confirmed mediating noncanonical DNA base-pairing, and 38 other candidate molecules were applied to demonstrate the feasibility of this protocol. We successfully identified seven stable noncanonical DNA duplex structures, and another eight novel structures with sub-stability. In addition, we discovered that hairpins at both ends can significantly stabilize the noncanonical DNA structures, providing a guideline to design small organic molecule-incorporated DNA structures. Such an efficient screening protocol will accelerate the design of alternative DNA-molecule architectures beyond Watson-Crick pairs. Considering the wide range of potential mediators, it will also facilitate applications such as noncovalent, highly dense loading of drug molecules in DNA-based delivery system and probe design for sensitive detection of certain molecules.

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.

Guanidinium–amino acid hydrogen-bonding interactions in protein crystal structures: implications for guanidinium-induced protein denaturation

Structural analysis of guanidinium–amino acid interaction pairs in protein crystal structures is coupled with an effective scheme for classifying the optimized pairs, to gain understanding of the guanidinium:protein hydrogen bonding modes.

On the prebiotic selection of nucleotide anomers: A computational study

Present-day known predominance of the β- over the α-anomers in nucleosides and nucleotides emerges from a thermodynamic analysis of their assembly from their components, i.e. bases, sugars, and a phosphate group. Furthermore, the incorporation of uracil into RNA and thymine into DNA rather than the other way around is also predicted from the calculations. An interplay of kinetics and thermodynamics must have driven evolutionary selection of life's building blocks. In this work, based on quantum chemical calculations, we focus on the latter control as a tool for “natural selection”.

Formations of bimolecular barbituric acid complexes through hydrogen bonding interactions: DFT analyses of structural and electronic features

Formations of bimolecular barbituric acid (BA) complexes through hydrogen-bonding (HB) interactions were investigated in this work. BA has been known as a starting compound of pharmaceutical compounds developments, in which the molecular and atomic features of parent BA in homo-paring with another BA molecule were investigated here. The models were optimized to reach the stabilized structures and their properties were evaluated at the molecular and atomic scales. Density functional theory (DFT) calculations were performed to provide required information for achieving the goal of this work. Six dimer models were obtained finally according to examining all possible starting dimers configurations for involving in optimization calculations. N-H . . . O and C-H . . . O interactions were also involved in dimers formations besides participation of the X-center of parent BA in interaction. Molecular and atomic scales features were evaluated for characterizing the dimers formations. As a consequence, several configurations of BA dimers were obtained showing the importance of performing such structural analyses for developing further compounds from BA.

Appraising the potency of small molecule inhibitors and their graphene surface-mediated organizational attributes on uric acid-melamine clusters

Uric acid (UA) and melamine (MM) crystallization in humans is associated with adverse medical conditions, including the germination of kidney stones, because of their low solubility. The growth of kidney stones, usually formed on renal papillary facades, is accomplished on the matrix-coated surface by the aggregation of preformed crystals or secondary crystal nucleation. Therefore, the effects of inhibitors such as theobromine (TB) and allopurinol (AP) on MM-UA aggregation are investigated by employing classical molecular dynamics simulations on a graphene surface. This impersonates the exact essence of the precipitation of kidney stones. The interaction between MM-UA is very intense and, thus, large clusters are formed on the surface. The presence of TB and AP will, however, substantially inhibit their aggregation. TB and AP significantly impede UA aggregation in particular. Therefore, lower order UA clusters are formed. These smaller UA clusters then pull a lower number of MM towards themselves, resulting in a smaller order UA-MM cluster. MM and UA aggregation on a 2D graphene surface is found to be spontaneous. There is no difference in these molecules' adsorption with a change in the force field parameters (i.e., GAFF and OPLS-AA) for graphene. Moreover, the greater the surface area of graphene, the more molecules are absorbed. The solute-surface van der Waals interaction energy plays a driving force in the adsorption of solute molecules on the surface. In addition, interactions like hydrogen bonding and π-stacking over the graphene surface involve binding all like molecules. These aggregated solute molecules strongly attract more like molecules until all solute molecules are adsorbed on the graphene surface, as estimated by enhanced sampling. The molecular origin of graphene exfoliation by MM is also described here. The present work helps to design novel kidney stone inhibitors.

Investigating Functionalization Impacts on Structural Features of Barbituric Acid Derivatives: DFT Approach

Density functional theory (DFT) calculations were performed to investigate electronic and structural properties of barbituric acid (BA) and sixtheen of its derivatives to show impacts of structural functionalization on the features of parent BA. The models were optimized and the minimum energy structures were confirmed by frequency calculations. Molecular and atomic descriptors were evaluated for the optimized models, in which the results of formation binding strength and molecular orbital features indicated significance of such functionalization processes on the observed properties. The highest occupied and the lowest unoccupied molecular orbitals (HOMO and LUMO) and their related parameters all indicated magnitudes of changes from one molecule to another one. Furthermore, atomic scale quadrupole coupling constants (Cq) were evaluated for the nitrogen and oxygen atoms of BA compounds showing significance of structural functionalization impacts on the atomic properties of parent BA. As a consequence, such structural analyses of BA compounds could show their characteristic features for further developments especially for their efficient pharmaceutical applications.

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.

Influence of the Number, Nature and Position of Methyl Posttranscriptional Modifications on Nucleobase Stacking in RNA

DFT calculations are employed to quantify the influence of the presence, number, nature, and position of posttranscriptional methylation on stacking strength of RNA bases. We carry out detailed potential energy scans of the variation in stacking energies with characteristic geometrical parameters in three categories of forty stacked dimers - canonical base homodimers (N||N), methylated base homodimers (mN||mN) and heterodimers of canonical bases and methylated counterparts (N||mN). Our analysis reveals that neutral methylation invariably enhances the stacking of bases. Further, N||mN stacking is stronger than mN||mN stacking and charged N||mN exhibit strongest stacking among all dimers. This indicates that methylations greatly enhance stacking when dispersed in RNA sequences containing identical bases. Comparison of stacks involving singly- and doubly-methylated purines reveal that incremental methylation enhances the stacking in neutral dimers. Although methylation at the carbon position of neutral pyrimidine dimers greatly enhances the stacking, methylation on the 5-membered ring imparts better stacking compared to methylation on the 6-membered ring in adenine dimers. However, methylation at the ring nitrogen (N¹ ) provides better stacking than the amino group (N² ) in guanine dimers. Our results thus highlight subtle structural effects of methylation on RNA base stacking and will enhance our understanding of the physicochemical principles of RNA structure and dynamics.

Investigation of corona discharge ionization of barbituric acid using ion mobility spectrometry along with quantum chemical calculations

This study aimed at examining atmospheric-pressure chemical ionization of barbituric acid through the corona discharge ion mobility spectrometry (CD-IMS) and the quantum chemical calculations. The results indicated two product ion peaks in the IMS spectrum of barbituric acid. The thermal decomposition of the barbituric acid sample was investigated by scanning the temperature of the injection port and analyzing the temporal evolution of the IMS peaks over elapsed time. It was found that the barbituric acid sample was not thermally decomposed in the injection port of the instrument. Experimental evidences were collected by changing the reactant ions, concentration of barbituric acid sample, and IMS cell temperature. The two observed peaks were then assigned to cationic form and oxygen protonated isomers of barbituric acid. The positions of the product ion peaks were explicated considering the dipole moments of the product ions.

Can modified DNA base pairs with chalcogen bonding expand the genetic alphabet? A combined quantum chemical and molecular dynamics simulation study

A comprehensive (DFT and MD) computational study is presented with the goal to design and analyze model chalcogen-bonded modified nucleobase pairs that replace one (i.e., A_XY:T, G:C_XY, G_XY:C) or two (G_XY:C_X'Y', X/X' = S, Se and Y/Y' = F, Cl, Br) Watson-Crick (WC) hydrogen bonds of the canonical A:T or G:C pair with chalcogen bond(s). DFT calculations on 18 base pair combinations that replace one WC hydrogen bond with a chalcogen bond reveal that the bases favorably interact in the gas phase (binding strengths up to -140 kJ mol^-1) and water (up to -85 kJ mol^-1). Although the remaining hydrogen bond(s) exhibits similar characteristics to those in the canonical base pairs, the structural features of the (Y-XO) chalcogen bond(s) change significantly with the identity of X and Y. The 36 doubly-substituted (G_XY:C_X'Y') base pairs have structural deviations from canonical G:C similar to those of the singly-substituted modifications (G:C_XY or G_XY:C). Furthermore, despite the replacement of two strong hydrogen bonds with chalcogen bonds, some G_XY:C_X'Y' pairs possess comparable binding energies (up to -132 kJ mol^-1 in the gas phase and up to -92 kJ mol^-1 in water) to the most stable G:C_XY or G_XY:C pairs, as well as canonical G:C. More importantly, G:C-modified pairs containing X = Se (high polarizability) and Y = F (high electronegativity) are the most stable, with comparable or slightly larger (by up to 13 kJ mol^-1) binding energies than G:C. Further characterization of the chalcogen bonding in all modified base pairs (AIM, NBO and NCI analyses) reveals that the differences in the binding energies of modified base pairs are mainly dictated by the differences in the strengths of their chalcogen bonds. Finally, MD simulations on DNA oligonucleotides containing the most stable chalcogen-bonded base pair from each of the four classifications (A_XY:T, G:C_XY, G_XY:C and G_XY:C_X'Y') reveal that the singly-modified G:C pairs best retain the local helical structure and pairing stability to a greater extent than the modified A:T pair. Overall, our study identifies two (G:C_SeF and G_SeF:C) promising pairs that retain chalcogen bonding in DNA and should be synthesized and further explored in terms of their potential to expand the genetic alphabet.

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.

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.

The Power Without the Glory: Multiple Roles of Hydrogen Peroxide in Mediating the Origin of Life

The hydrogen peroxide (HP) crucible hypothesis proposed here holds that life began in a localized environment on Earth that was perfused with a flow of hydrogen peroxide from a sustained external source, which powered and mediated molecular evolution and the protocellular RNA world. In this article, we consolidate and review recent evidence, both circumstantial and tested in simulation in our work and in the laboratory in others' work, for its multiple roles in the evolution of the first living systems: (1) it provides a periodic power source as the thiosulfate-hydrogen peroxide (THP) redox oscillator, (2) it may act as an agent of molecular change and evolution and mediator of homochirality, and (3) the THP oscillator, subject to Brownian input perturbations, produces a weighted distribution of output thermal fluctuations that favor polymerization and chemical diversification over chemical degradation and simplification. The hypothesis can help to clarify the hero and villain roles of hydrogen peroxide in cell function, and on the singularity of life: of necessity, life evolved early an armory of catalases, the continuing, and all-pervasive presence of which prevents hydrogen peroxide from accumulating anywhere in sufficient quantities to host a second origin. The HP crucible hypothesis is radical, but based on well-known chemistry and physics, it is eminently testable in the laboratory, and many of our simulations provide recipes for such experiments.

The A·T(rWC)/A·T(H)/A·T(rH) ↔ A·T*(rwWC)/A·T*(wH)/A·T*(rwH) mutagenic tautomerization via sequential proton transfer: a QM/QTAIM study

In this study for the first time we have revealed by QM and QTAIM calculations at the MP2/aug-cc-pVDZ//B3LYP/6-311++G(d,p) level of QM theory the novel routes of the mutagenic tautomerization of three biologically important A·T DNA base pairs – reverse Watson–Crick A·T(rWC), Hoogsteen A·T(H) and reverse Hoogsteen A·T(rH) – followed by their rebuilding into the wobble (w) A·T*(rw_WC), A·T*(w_H) and A·T*(rw_H) base mispairs by the participation of the mutagenic tautomers of the DNA bases (denoted by asterisk) and vice versa, thus complementing the physico-chemical property of the canonical A·T(WC) Watson–Crick DNA base pair reported earlier (Brovarets' et al., RSC Adv., 2015, 5, 99594–99605). These non-dissociative tautomeric transformations in the classical A·T(rWC), A·T(H) and A·T(rH) DNA base pairs proceed similarly to the canonical A·T(WC) DNA base pair via the intrapair sequential proton transfer with shifting towards major or minor grooves of DNA followed by further double proton transfer along the intermolecular H-bonds and are controlled by the plane symmetric and highly stable transition states – tight ion pairs formed by the A⁺ nucleobase, protonated by the N1/N7 nitrogen atoms, and T⁻ nucleobase, deprotonated by the N3H imino group. Comparison of the estimated populations of the tautomerised states (10⁻²¹ to 10⁻¹⁴) with similar characteristics for the canonical A·T(WC) DNA base pair (10⁻⁸ to 10⁻⁷) leads authors to the conclusion, that only a base pair with WC architecture can be a building block of the DNA macromolecule as a genetic material, which is able for the evolutionary self-development. Among all four classical DNA base pairs, only A·T(WC) DNA base pair can ensure the proper rate of the spontaneous point errors of replication in DNA.

We discovered tautomeric wobbling of the classical A·T DNA base pairs. This data evidence, that only a base pair with Watson–Crick architecture can be a building block of the DNA macromolecule as a genetic material, which is able for the evolutionary self-development.