Tautomeric Equilibria of Nucleobases in the Hachimoji Expanded Genetic Alphabet

Electronic and Nuclear Quantum Effects on Proton Transfer Reactions of Guanine-Thymine (G-T) Mispairs Using Combined Quantum Mechanical/Molecular Mechanical and Machine Learning Potentials

Rare tautomeric forms of nucleobases can lead to Watson-Crick-like (WC-like) mispairs in DNA, but the process of proton transfer is fast and difficult to detect experimentally. NMR studies show evidence for the existence of short-time WC-like guanine-thymine (G-T) mispairs; however, the mechanism of proton transfer and the degree to which nuclear quantum effects play a role are unclear. We use a B-DNA helix exhibiting a wGT mispair as a model system to study tautomerization reactions. We perform ab initio (PBE0/6-31G*) quantum mechanical/molecular mechanical (QM/MM) simulations to examine the free energy surface for tautomerization. We demonstrate that while the ab initio QM/MM simulations are accurate, considerable sampling is required to achieve high precision in the free energy barriers. To address this problem, we develop a QM/MM machine learning potential correction (QM/MM-ΔMLP) that is able to improve the computational efficiency, greatly extend the accessible time scales of the simulations, and enable practical application of path integral molecular dynamics to examine nuclear quantum effects. We find that the inclusion of nuclear quantum effects has only a modest effect on the mechanistic pathway but leads to a considerable lowering of the free energy barrier for the GT*⇌G*T equilibrium. Our results enable a rationalization of observed experimental data and the prediction of populations of rare tautomeric forms of nucleobases and rates of their interconversion in B-DNA.

Investigation of the stability of D5SIC-DNAM-incorporated DNA duplex in Taq polymerase binary system: a systematic classical MD approach

DNA polymerases are fundamental enzymes that play a crucial role in processing DNA with high fidelity and accuracy ensuring the faithful transmission of genetic information. The recognition of unnatural base pairs (UBPs) by polymerases, enabling their replication, represents a significant and groundbreaking discovery with profound implications for genetic expansion. Romesberg et al. examined the impact of DNA containing 2,6-dimethyl-2H-isoquiniline-1-thione: D5SIC (DS) and 2-methoxy-3-methylnaphthalene: DNAM (DN) UBPs bound to T. aquaticus DNA polymerase (Taq) through crystal structure analysis. Here, we have used polarizable and nonpolarizable classical molecular dynamics (MD) simulations to investigate the structural aspects and stability of Taq in complex with a DNA duplex including a DS-DN pair in the terminal 3' and 5' positions. Our results suggest that the flexibility of UBP-incorporated DNA in the terminal position is arrested by the polymerase, thus preventing fraying and mispairing. Our investigation also reveals that the UBP remains in an intercalated conformation inside the active site, exhibiting two distinct orientations in agreement with experimental findings. Our analysis pinpoints particular residues responsible for favorable interactions with the UBP, with some relying on van der Waals interactions while other on Coulombic forces.

Investigation of dynamical flexibility of D5SIC-DNAM inside DNA duplex in aqueous solution: a systematic classical MD approach

Incorporation of artificial 3rd base pairs (unnatural base pairs, UBPs) has emerged as a fundamental technique in pursuit of expanding the genetic alphabet. 2,6-Dimethyl-2H-isoquiniline-1-thione: D5SIC (DS) and 2-methoxy-3-methylnaphthalene: DNAM (DN), a potential unnatural base pair (UBP) developed by Romesberg and colleagues, has been shown to have remarkable capability for replication within DNA. Crystal structures of a Taq polymerase/double-stranded DNA (ds-DNA) complex containing a DS-DN pair in the 3' terminus showed a parallelly stacked geometry for the pre-insertion, and an intercalated geometry for the post-insertion structure. Unconventional orientations of DS-DN inside a DNA duplex have inspired scientists to investigate the conformational orientations and structural properties of UBP-incorporated DNA. In recent years, computational simulations have been used to investigate the geometry of DS-DN within the DNA duplex; nevertheless, unresolved questions persist owing to inconclusive findings. In this work, we investigate the structural and dynamical properties of DS and DN inside a ds-DNA strand in aqueous solution considering both short and long DNA templates using polarizable, and non-polarizable classical MD simulations. Flexible conformational change of UBP with major populations of Watson-Crick-Franklin (WCF) and three distinct non-Watson-Crick-Franklin (nWCFP1, nWCFP2, nWCFO) conformations through intra and inter-strand flipping have been observed. Our results suggest that a dynamical conformational change leads to the production of diffierent conformational distribution for the systems. Simulations with a short ds-DNA duplex suggest nWCF (P1 and O) as the predominant structures, whereas long ds-DNA duplex simulations indicate almost equal populations of WCF, nWCFP1, nWCFO. DS-DN in the terminal position is found to be more flexible with occasional mispairing and fraying. Overall, these results suggest flexibility and dynamical conformational change of the UBP as well as indicate varied conformational distribution irrespective of starting orientation of the UBP and length og DNA strand.

Recent developments and applications of reference interaction site model self-consistent field with constrained spatial electron density (RISM-SCF-cSED): A hybrid model of quantum chemistry and integral equation theory of molecular liquids

The significance of solvent effects in electronic structure calculations has long been noted, and various methods have been developed to consider this effect. The reference interaction site model self-consistent field with constrained spatial electron density (RISM-SCF-cSED) is a hybrid model that combines the integral equation theory of molecular liquids with quantum chemistry. This method can consider the statistically convergent solvent distribution at a significantly lower cost than molecular dynamics simulations. Because the RISM theory explicitly considers the solvent structure, it performs well for systems where hydrogen bonds are formed between the solute and solvent molecules, which is a challenge for continuum solvent models. Taking advantage of being founded on the variational principle, theoretical developments have been made in calculating various properties and incorporating electron correlation effects. In this review, we organize the theoretical aspects of RISM-SCF-cSED and its distinctions from other hybrid methods involving integral equation theories. Furthermore, we carefully present its progress in terms of theoretical developments and recent applications.

A unified Watson-Crick geometry drives transcription of six-letter expanded DNA alphabets by E. coli RNA polymerase

Artificially Expanded Genetic Information Systems (AEGIS) add independently replicable unnatural nucleotide pairs to the natural G:C and A:T/U pairs found in native DNA, joining the unnatural pairs through alternative modes of hydrogen bonding. Whether and how AEGIS pairs are recognized and processed by multi-subunit cellular RNA polymerases (RNAPs) remains unknown. Here, we show that E. coli RNAP selectively recognizes unnatural nucleobases in a six-letter expanded genetic system. High-resolution cryo-EM structures of three RNAP elongation complexes containing template-substrate UBPs reveal the shared principles behind the recognition of AEGIS and natural base pairs. In these structures, RNAPs are captured in an active state, poised to perform the chemistry step. At this point, the unnatural base pair adopts a Watson-Crick geometry, and the trigger loop is folded into an active conformation, indicating that the mechanistic principles underlying recognition and incorporation of natural base pairs also apply to AEGIS unnatural base pairs. These data validate the design philosophy of AEGIS unnatural basepairs. Further, we provide structural evidence supporting a long-standing hypothesis that pair mismatch during transcription occurs via tautomerization. Together, our work highlights the importance of Watson-Crick complementarity underlying the design principles of AEGIS base pair recognition.

Tautomeric equilibrium and spectroscopic properties of 8-azaguanine revealed by quantum chemistry methods

8-azaguanine is a triazolopyrimidine nucleobase analog possessing potent antibacterial and antitumor activities, and it has been implicated as a lead molecule in cancer and malaria therapy. Its intrinsic fluorescence properties can be utilized for monitoring its interactions with biological polymers like proteins or nucleic acids. In order to better understand these interactions, it is important to know the tautomeric equilibrium of this compound. In this work, the tautomeric equilibrium of all natural neutral and anionic compound forms (except highly improbable imino-enol tautomers) as well as their methyl derivatives and ribosides was revealed by quantum chemistry methods. It was shown that, as expected, tautomers protonated at positions 1 and 9 dominate neutral forms both in gas phase and in aqueous solution. 8-azaguanines methylated at any position of the triazole ring are protonated at position 1. The computed vertical absorption and emission energies are in very good agreement with the experimental data. They confirm the validity of the assumption that replacing the proton with the methyl group does not significantly change the positions of absorption and fluorescence peaks.

How proton transfer impacts hachimoji DNA

Hachimoji DNA is a synthetic nucleic acid extension of DNA, formed by an additional four bases, Z, P, S, and B, that can encode information and sustain Darwinian evolution. In this paper, we aim to look into the properties of hachimoji DNA and investigate the probability of proton transfer between the bases, resulting in base mismatch under replication. First, we present a proton transfer mechanism for hachimoji DNA, analogous to the one presented by Löwdin years prior. Then, we use density functional theory to calculate proton transfer rates, tunnelling factors and the kinetic isotope effect in hachimoji DNA. We determined that the reaction barriers are sufficiently low that proton transfer is likely to occur even at biological temperatures. Furthermore, the rates of proton transfer of hachimoji DNA are much faster than in Watson-Crick DNA due to the barrier for Z-P and S-B being 30% lower than in G-C and A-T. Suggesting that proton transfer occurs more frequently in hachimoji DNA than canonical DNA, potentially leading to a higher mutation rate.

QDπ: A Quantum Deep Potential Interaction Model for Drug Discovery

We report QDπ-v1.0 for modeling the internal energy of drug molecules containing H, C, N, and O atoms. The QDπ model is in the form of a quantum mechanical/machine learning potential correction (QM/Δ-MLP) that uses a fast third-order self-consistent density-functional tight-binding (DFTB3/3OB) model that is corrected to a quantitatively high-level of accuracy through a deep-learning potential (DeepPot-SE). The model has the advantage that it is able to properly treat electrostatic interactions and handle changes in charge/protonation states. The model is trained against reference data computed at the ωB97X/6-31G* level (as in the ANI-1x data set) and compared to several other approximate semiempirical and machine learning potentials (ANI-1x, ANI-2x, DFTB3, MNDO/d, AM1, PM6, GFN1-xTB, and GFN2-xTB). The QDπ model is demonstrated to be accurate for a wide range of intra- and intermolecular interactions (despite its intended use as an internal energy model) and has shown to perform exceptionally well for relative protonation/deprotonation energies and tautomers. An example application to model reactions involved in RNA strand cleavage catalyzed by protein and nucleic acid enzymes illustrates QDπ has average errors less than 0.5 kcal/mol, whereas the other models compared have errors over an order of magnitude greater. Taken together, this makes QDπ highly attractive as a potential force field model for drug discovery.

Rethinking nucleic acids from their origins to their applications

Reviewed are three decades of synthetic biology research in our laboratory that has generated alternatives to standard DNA and RNA as possible informational systems to support Darwinian evolution, and therefore life, and to understand their natural history, on Earth and throughout the cosmos. From this, we have learned that: • the core structure of nucleic acids appears to be a natural outcome of non-biological chemical processes probably in constrained, intermittently irrigated, sub-aerial aquifers on the surfaces of rocky planets like Earth and/or Mars approximately 4.36 ± 0.05 billion years ago; • however, this core is not unique. Synthetic biology has generated many different molecular systems able to support the evolution of molecular information; • these alternatives to standard DNA and RNA support biotechnology, including DNA synthesis, human diagnostics, biomedical research and medicine; • in particular, they support laboratory in vitro evolution (LIVE) with performance to generate catalysts at least 104-105 fold better than standard DNA libraries, enhancing access to receptors and catalysts on demand. Coupling nanostructures to the products of LIVE with expanded DNA offers new approaches for disease therapy; and • nevertheless, a polyelectrolyte structure and size regular building blocks are required for any informational polymer to support Darwinian evolution. These features serve as universal and agnostic biosignatures, useful for seeking life throughout the Solar System. This article is part of the theme issue 'Reactivity and mechanism in chemical and synthetic biology'.

The Base Pairs of Isoguanine and 8-Aza-7-deazaisoguanine with 5-Methylisocytosine as Targets for DNA Functionalization

The isoguanine-isocytosine base pair (isoG-isoC) represents an important expansion of the DNA coding system. The base pair is more stable than the canonical adenine-thymine or guanine-cytosine pairs. However, nothing is known on the functionalization of the noncanonical isoG-isoC pair at the isoguanine site. In this work, functionalization of the isoG-isoC and the isosteric base pair that contains 8-aza-7-deazaisoguanine in place of isoguanine is studied. Short ethynyl, more space demanding octadiynyl, and dendritic tripropargylamine residues attached to the isoG-isoC base pairs were introduced to oligonucleotides. 12-mer duplexes were formed by hybridization with single base pair modification. The use of the two modified nucleobases gave us the freedom to shift nucleobase substituents within the major groove of double helical DNA. Clickable side chains at position-7 stabilize the base pair, whereas 8-substituents reduce its stability strongly. The weak isoguanine-thymine or 8-aza-7-deazaisoguanine-thymine base pairs show a similar sensitivity to the position of nucleobase functionalization as base pair matches formed with 5-methylisocytosine. CD spectra of all modified duplexes display the typical shape of a B-DNA with only marginal changes. Fluorescent pyrene labeled DNA with long, short, and branched linkers was generated using click chemistry. Pyrene click adducts with long linkers are essential to maintain or to increase base pair stability. Labeled duplexes are more fluorescent than corresponding single strands. For the dendritic linker excimer emission was observed for single strands but only monomer emission in duplexes.

Photophysical Characterization of Isoguanine in a Prebiotic-Like Environment

It is intriguing how a mixture of organic molecules survived the prebiotic UV fluxes and evolved into the actual genetic building blocks. Scientists are trying to shed light on this issue by synthesizing nucleic acid monomers and their analogues under prebiotic Era-like conditions and by exploring their excited state dynamics. To further add to this important body of knowledge, this study discloses new insights into the photophysical properties of protonated isoguanine, an isomorph of guanine, using steady-state and femtosecond broadband transient absorption spectroscopies, and quantum mechanical calculations. Protonated isoguanine decays in ultrafast time scales following 292 nm excitation, consistently with the barrierless paths connecting the bright S₁ (ππ*) state with different internal conversion funnels. Complementary calculations for neutral isoguanine predict similar photophysical properties. These results demonstrate that protonated isoguanine can be considered photostable in contrast to protonated guanine, which exhibits 40-fold longer excited state lifetimes.

Resolving Ultrafast Photoinitiated Dynamics of the Hachimoji 5-Aza-7-Deazaguanine Nucleobase: Impact of Synthetically Expanding the Genetic Alphabet

The guanine derivative, 5-aza-7-deazaguanine (^5N7C G) has recently been proposed as one of four unnatural bases, termed Hachimoji (8-letter) to expand the genetic code. We apply steady-state and time-resolved spectroscopy to investigate its electronic relaxation mechanism and probe the effect of atom substitution on the relaxation mechanism in polar protic and polar aprotic solvents. Mapping of the excited state potential energy surfaces is performed, from which the critical points are optimized by using the state-of-art Extended Multi-State Complete Active Space Second-Order Perturbation Theory. It is demonstrated that excitation to the lowest energy ¹ ππ* state of ^5N7C G results in complex dynamics leading to ca. 10 to 30-fold slower relaxation (depending on solvent) compared to guanine. A significant conformational change occurs at the S₁ minimum, resulting in a 10-fold greater fluorescence quantum yield compared to guanine. The fluorescence quantum yield and S₁ decay lifetime increase going from water to acetonitrile to propanol. The solvent-dependent results are supported by the quantum chemical calculations showing an increase in the energy barrier between the S₁ minimum and the S₁ /S₀ conical intersection going from water to propanol. The longer lifetimes might make ^5N7C G more photochemical active to adjacent nucleobases than guanine or other nucleobases within DNA.

Investigations Into Chemically Stabilized Four-Letter DNA for DNA-Encoded Chemistry

DNA-encoded libraries are a prime technology for target-based small molecule screening. Native DNA used as genetic compound barcode is chemically vulnerable under many reaction conditions. DNA barcodes that are composed of pyrimidine nucleobases, 7-deazaadenine, and 7-deaza-8-azaguanine have been investigated for their suitability for encoded chemistry both experimentally and computationally. These four-letter barcodes were readily ligated by T4 ligation, amplifiable by Taq polymerase, and the resultant amplicons were correctly sequenced. Chemical stability profiling showed a superior chemical stability compared to native DNA, though higher susceptibility to depurination than a three-letter code based on pyrimidine DNA and 7-deazaadenine.

Oxidation Enhances Binding of Extrahelical 5-Methyl-Cytosines by Thymine DNA Glycosylase

The DNA repair protein thymine DNA glycosylase (TDG) removes mispaired or damaged bases, such as oxidized methyl-cytosine, from DNA by cleavage of the glycosidic bond between the sugar and the target base flipped into the enzyme's active site. The enzyme is active against formyl-cytosine and carboxyl-cytosine, whereas the lower oxidized hydroxymethyl-cytosine and methyl-cytosine itself are not processed by the enzyme. Molecular dynamics simulations with thermodynamic integration of TDG complexed to DNA carrying one of four different (oxidized) methyl-cytosine bases in extrahelcial conformation, methyl-cytosine (mC), hydroxymethyl-cytosine (hmC), formyl-cytosine (fC), or carboxyl-cytosine (caC), show a more favorable binding affinity of the higher oxidized forms, fC and caC, than the nonsubstrate bases hmC and mC. Despite rather comparable, reaction-competent conformations of the flipped bases in the active site of the enzyme, more and stronger interactions with active site residues account for the preferred binding of the higher oxidized bases. Binding of the negatively charged caC and the neutral fC are strengthened by interactions with positively charged His151. Our calculated proton affinities find this protonation state of His151 the preferred one in the presence of caC and conceivable in the presence of fC as well as increasing the binding affinity toward the two bases. Discrimination of the substrate bases is further achieved by the backbone of Tyr152 that forms a strong hydrogen bond to the carboxyl and formyl oxygen atoms of caC and fC, respectively, a contact that is completely lacking in mC and much weaker in hmC. Overall, our computational results indicate that the enzyme discriminates the different oxidation forms of methyl-cytosine already at the formation of the extrahelical complexes.

Rare Tautomers of Artificially Expanded Genetic Letters and their Effects on the Base pair Stabilities

To expand the existing genetic letters, it is necessary to design robust nucleotides that can function naturally in living cells. Therefore, it is desirable to examine the roles of recently proposed second-generation artificially expanded genetic letters in producing stable duplex DNA. Here, a reliable dispersion-corrected density functional theory method is used to understand the electronic structures and properties of different rare tautomers of proposed expanded genetic letters and their effects on the base pair stabilities in the duplex DNA. It is found that the rare tautomers are not only stable in the aqueous medium but can also base pair with natural bases to produce stable mispairs. Except for J and V, all the artificial genetic letters are found to produce mispairs that are about 1-7 kcal/mol more stable than their complementary counterparts. They are also appreciably more stable than the naturally occurring G:C, A:T, and G:T pairs. The higher base pair stabilities are found to be mainly because of the polarity of monomers and attractive electrostatic interactions.

New chemistries and enzymes for synthetic genetics

Beyond the natural nucleic acids DNA and RNA, nucleic acid chemistry has unlocked a whole universe of modifications to their canonical chemical structure, which can in various ways modify and enhance nucleic acid function and utility for applications in biotechnology and medicine. Unlike the natural modifications of tRNA and rRNA or the epigenetic modifications in mRNA and genomic DNA, these altered chemistries are not found in nature and therefore these molecules are referred to as xeno-nucleic acids (XNAs). In this review we aim to focus specifically on recent progress in a subsection of this vast field-synthetic genetics-concerned with encoded synthesis, reverse transcription, and evolution of XNAs.

Interaction of Thymine DNA Glycosylase with Oxidised 5-Methyl-cytosines in Their Amino- and Imino-Forms

Thymine DNA Glycosylase (TDG) is an enzyme of the base excision repair mechanism and removes damaged or mispaired bases from DNA via hydrolysis of the glycosidic bond. Specificity is of high importance for such a glycosylase, so as to avoid the damage of intact DNA. Among the substrates reported for TDG are mispaired uracil and thymine but also formyl-cytosine and carboxyl-cytosine. Methyl-cytosine and hydroxylmethyl-cytosine are, in contrast, not processed by the TDG enzyme. We have in this work employed molecular dynamics simulations to explore the conformational dynamics of DNA carrying a formyl-cytosine or carboxyl-cytosine and compared those to DNA with the non-cognate bases methyl-cytosine and hydroxylmethyl-cytosine, as amino and imino tautomers. Whereas for the mispairs a wobble conformation is likely decisive for recognition, all amino tautomers of formyl-cytosine and carboxyl-cytosine exhibit the same Watson–Crick conformation, but all imino tautomers indeed form wobble pairs. The conformational dynamics of the amino tautomers in free DNA do not exhibit differences that could be exploited for recognition, and also complexation to the TDG enzyme does not induce any alteration that would indicate preferable binding to one or the other oxidised methyl-cytosine. The imino tautomers, in contrast, undergo a shift in the equilibrium between a closed and a more open, partially flipped state, towards the more open form upon complexation to the TDG enzyme. This stabilisation of the more open conformation is most pronounced for the non-cognate bases methyl-cytosine and hydroxyl-cytosine and is thus not a likely mode for recognition. Moreover, calculated binding affinities for the different forms indicate the imino forms to be less likely in the complexed DNA. These findings, together with the low probability of imino tautomers in free DNA and the indifference of the complexed amino tautomers, suggest that discrimination of the oxidised methyl-cytosines does not take place in the initial complex formation.

Chemically Stabilized DNA Barcodes for DNA-Encoded Chemistry

DNA‐encoded compound libraries are a widely used small molecule screening technology. One important aim in library design is the coverage of chemical space through structurally diverse molecules. Yet, the chemical reactivity of native DNA barcodes limits the toolbox of reactions for library design. Substituting the chemically vulnerable purines by 7‐deazaadenine, which exhibits tautomerization stability similar to natural adenine with respect to the formation of stable Watson–Crick pairs, yielded ligation‐competent, amplifiable, and readable DNA barcodes for encoded chemistry with enhanced stability against protic acid‐ and metal ion‐promoted depurination. The barcode stability allowed for straightforward translation of 16 exemplary reactions that included isocyanide multicomponent reactions, acid‐promoted Pictet–Spengler and Biginelli reactions, and metal‐promoted pyrazole syntheses on controlled pore glass‐coupled barcodes for diverse DEL design. The Boc protective group of reaction products offered a convenient handle for encoded compound purification.

Quantum-mechanical property prediction of solvated drug molecules: what have we learned from a decade of SAMPL blind prediction challenges?

Joint academic-industrial projects supporting drug discovery are frequently pursued to deploy and benchmark cutting-edge methodical developments from academia in a real-world industrial environment at different scales. The dimensionality of tasks ranges from small molecule physicochemical property assessment over protein-ligand interaction up to statistical analyses of biological data. This way, method development and usability both benefit from insights gained at both ends, when predictiveness and readiness of novel approaches are confirmed, but the pharmaceutical drug makers get early access to novel tools for the quality of drug products and benefit of patients. Quantum-mechanical and simulation methods particularly fall into this group of methods, as they require skills and expense in their development but also significant resources in their application, thus are comparatively slowly dripping into the realm of industrial use. Nevertheless, these physics-based methods are becoming more and more useful. Starting with a general overview of these and in particular quantum-mechanical methods for drug discovery we review a decade-long and ongoing collaboration between Sanofi and the Kast group focused on the application of the embedded cluster reference interaction site model (EC-RISM), a solvation model for quantum chemistry, to study small molecule chemistry in the context of joint participation in several SAMPL (Statistical Assessment of Modeling of Proteins and Ligands) blind prediction challenges. Starting with early application to tautomer equilibria in water (SAMPL2) the methodology was further developed to allow for challenge contributions related to predictions of distribution coefficients (SAMPL5) and acidity constants (SAMPL6) over the years. Particular emphasis is put on a frequently overlooked aspect of measuring the quality of models, namely the retrospective analysis of earlier datasets and predictions in light of more recent and advanced developments. We therefore demonstrate the performance of the current methodical state of the art as developed and optimized for the SAMPL6 pKa and octanol-water log P challenges when re-applied to the earlier SAMPL5 cyclohexane-water log D and SAMPL2 tautomer equilibria datasets. Systematic improvement is not consistently found throughout despite the similarity of the problem class, i.e. protonation reactions and phase distribution. Hence, it is possible to learn about hidden bias in model assessment, as results derived from more elaborate methods do not necessarily improve quantitative agreement. This indicates the role of chance or coincidence for model development on the one hand which allows for the identification of systematic error and opportunities toward improvement and reveals possible sources of experimental uncertainty on the other. These insights are particularly useful for further academia-industry collaborations, as both partners are then enabled to optimize both the computational and experimental settings for data generation.

Unnatural bases for recognition of noncoding nucleic acid interfaces

The notion of using synthetic heterocycles instead of the native bases to interface with DNA and RNA has been explored for nearly 60 years. Unnatural bases compatible with the DNA/RNA coding interface have the potential to expand the genetic code and co-opt the machinery of biology to access new macromolecular function; accordingly, this body of research is core to synthetic biology. While much of the literature on artificial bases focuses on code expansion, there is a significant and growing effort on docking synthetic heterocycles to noncoding nucleic acid interfaces; this approach seeks to illuminate major processes of nucleic acids, including regulation of transcription, translation, transport, and transcript lifetimes. These major avenues of research at the coding and noncoding interfaces have in common fundamental principles in molecular recognition. Herein, we provide an overview of foundational literature in biophysics of base recognition and unnatural bases in coding to provide context for the developing area of targeting noncoding nucleic acid interfaces with synthetic bases, with a focus on systems developed through iterative design and biophysical study.

Building better polymerases: Engineering the replication of expanded genetic alphabets

DNA polymerases are today used throughout scientific research, biotechnology, and medicine, in part for their ability to interact with unnatural forms of DNA created by synthetic biologists. Here especially, natural DNA polymerases often do not have the "performance specifications" needed for transformative technologies. This creates a need for science-guided rational (or semi-rational) engineering to identify variants that replicate unnatural base pairs (UBPs), unnatural backbones, tags, or other evolutionarily novel features of unnatural DNA. In this review, we provide a brief overview of the chemistry and properties of replicative DNA polymerases and their evolved variants, focusing on the Klenow fragment of Taq DNA polymerase (Klentaq). We describe comparative structural, enzymatic, and molecular dynamics studies of WT and Klentaq variants, complexed with natural or noncanonical substrates. Combining these methods provides insight into how specific amino acid substitutions distant from the active site in a Klentaq DNA polymerase variant (ZP Klentaq) contribute to its ability to replicate UBPs with improved efficiency compared with Klentaq. This approach can therefore serve to guide any future rational engineering of replicative DNA polymerases.