Toward an Expanded Genome: Structural and Computational Characterization of an Artificially Expanded Genetic Information System

Solid-state synthesis of polyfunctionalized 2-pyridones and conjugated dienes

Functionalized 2-pyridones are important biologically active compounds, DNA base analogues and synthetic intermediates. Herein, we report a simple, green, solid-state synthesis of differently substituted 2-pyridones. It starts from commercially available amines and activated alkynes, uses silica gel (15%Cs2CO3/SiO2) as the solid phase and a reaction vial as the only equipment. If necessary, heating is performed in a laboratory oven. Since most reactions are completed within a few hours, no additional energy consumption is required. The syntheses do not require solvents and other reagents and are easily monitored by standard analytical techniques. The atom economy is high, since all atoms of reactants are present in the products and EtOH is the only by-product. The syntheses produce polyfunctionalized conjugated dienes as the only intermediates, which are also important building blocks.

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.

Theoretical investigation of tautomerism of 2- and 4-pyridones: origin, substituent and solvent effects

Computational investigation at the BHandHLYP/6-311+G(d,p) level of theory of the gas-phase tautomerism of 2- and 4-pyridones confirmed the slight prevalence of lactim in the case of the former, but its dominance in the case of the latter, as shown previously. Examination of aromaticity by using HOMA, EDDB, NBOdel, NICS and AICD led to the conclusion that tautomerization of 4-pyridone results in greater aromaticity gain. It is also driven by the Pauli repulsion relief, which was revealed by the tautomerization energy decomposition analysis. By contrast, in the case of 2-pyridone, lactim is favoured by orbital and electrostatic interactions and disfavoured by the Pauli repulsion. Aromaticity gain in this case is smaller. The position of the tautomeric equilibrium can be modulated by substituent inductive effects (Cl and F), inductive and resonance effects (NH2 and NO2), hydrogen bonding (NO2), and medium polarity, the increase of which increases lactam population.

Visible Helicity Induction and Memory in Polyallene toward Circularly Polarized Luminescence, Helicity Discrimination, and Enantiomer Separation

Inspired by biological helices (e.g., DNA), artificial helical polymers have attracted intense attention. However, precise synthesis of one-handed helices from achiral materials remains a formidable challenge. Herein, a series of achiral poly(biphenyl allene)s with controlled molar mass and low dispersity were prepared and induced into one-handed helices using chiral amines and alcohols. The induced one-handed helix was simultaneously memorized, even after the chiral inducer was removed. The switchable induction processes were visible to naked eye; the achiral polymers exhibited blue emission (irradiated at 365 nm), whereas the induced one-handed helices exhibited cyan emission with clear circularly polarized luminescence. The induced helices formed stable gels in various solvents with helicity discrimination ability: the same-handed helix gels were self-healing, whereas the gels of opposite-handed helicity were self-sorted. Moreover, the induced helices could separate enantiomers via enantioselective crystallization with high efficiency and switchable enantioselectivity.

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.

Crystal structures of 'ALternative Isoinformational ENgineered' DNA in B-form

The first structural model of duplex DNA reported in 1953 by Watson & Crick presented the double helix in B-form, the form that genomic DNA exists in much of the time. Thus, artificial DNA seeking to mimic the properties of natural DNA should also be able to adopt B-form. Using a host-guest system in which Moloney murine leukemia virus reverse transcriptase serves as the host and DNA as the guests, we determined high-resolution crystal structures of three complexes including 5'-CTTBPPBBSSZZSAAG, 5'-CTTSSPBZPSZBBAAG and 5'-CTTZZPBSBSZPPAAG with 10 consecutive unnatural nucleobase pairs in B-form within self-complementary 16 bp duplex oligonucleotides. We refer to this ALternative Isoinformational ENgineered (ALIEN) genetic system containing two nucleobase pairs (P:Z, pairing 2-amino-imidazo-[1,2-a]-1,3,5-triazin-(8H)-4-one with 6-amino-5-nitro-(1H)-pyridin-2-one, and B:S, 6-amino-4-hydroxy-5-(1H)-purin-2-one with 3-methyl-6-amino-pyrimidin-2-one) as ALIEN DNA. We characterized both position- and sequence-specific helical, nucleobase pair and dinucleotide step parameters of P:Z and B:S pairs in the context of B-form DNA. We conclude that ALIEN DNA exhibits structural features that vary with sequence. Further, Z can participate in alternative stacking modes within a similar sequence context as captured in two different structures. This finding suggests that ALIEN DNA may have a larger repertoire of B-form structures than natural DNA. This article is part of the theme issue 'Reactivity and mechanism in chemical and synthetic biology'.

Macromolecular Information Transfer

Macromolecular information transfer can be defined as the process by which a coded monomer sequence is communicated from one macromolecule to another. In such a transfer process, the information sequence can be kept identical, transformed into a complementary sequence or even translated into a different molecular language. Such mechanisms are crucial in biology and take place in DNA→DNA replication, DNA→RNA transcription and RNA→protein translation. In fact, there would be no life on Earth without macromolecular information transfer. Mimicking such processes with synthetic macromolecules would also be of major scientific relevance because it would open up new avenues for technological applications (e.g. data storage and processing) but also for the creation of artificial life. In this important context, this minireview summarizes recent research about information transfer in synthetic oligomers and polymers. Medium- and long-term perspectives are also discussed.

Visualizing "Alternative Isoinformational Engineered" DNA in A- and B-Forms at High Resolution

A fundamental property of DNA built from four informational nucleotide units (GCAT) is its ability to adopt different helical forms within the context of the Watson-Crick pair. Well-characterized examples include A-, B-, and Z-DNA. For this study, we created an isoinformational biomimetic polymer, built (like standard DNA) from four informational "letters", but with the building blocks being artificial. This ALternative Isoinformational ENgineered (ALIEN) DNA was hypothesized to support two nucleobase pairs, the P:Z pair matching 2-amino-imidazo-[1,2a]-1,3,5-triazin-[8H]-4-one with 6-amino-3-5-nitro-1H-pyridin-2-one and the B:S pair matching 6-amino-4-hydroxy-5-1H-purin-2-one with 3-methyl-6-amino-pyrimidin-2-one. We report two structures of ALIEN DNA duplexes at 1.2 Å resolution and a third at 1.65 Å. All of these are built from a single self-complementary sequence (5'-CTSZZPBSBSZPPBAG) that includes 12 consecutive ALIEN nucleotides. We characterized the helical, nucleobase pair, and dinucleotide step parameters of ALIEN DNA in these structures. In addition to showing that ALIEN pairs retain basic Watson-Crick pairing geometry, two of the ALIEN DNA structures are characterized as A-form DNA and one as B-form DNA. We identified parameters that map differences effecting the transition between the two helical forms; these same parameters distinguish helical forms of isoinformational natural DNA. Collectively, our analyses suggest that ALIEN DNA retains essential structural features of natural DNA, not only its information density and Watson-Crick pairing but also its ability to adopt two canonical forms.

The role of nitro group on the excited-state relaxation mechanism of P-Z base pair

DNAs' photostability is significant to the normal function of organisms. P-Z is a hydrogen bonded artificial DNA base pair, where P and Z represent 2-amino-imidazo[1,2-a]-1,3,5-triazin-4(8H)one and 6-amino-5nitro-2(1H)-pyridone, respectively. The excited-state relaxation mechanism of P-Z pair is investigated using static TDDFT calculations combined with the non-adiabatic dynamic simulations at TDDFT level. The roles of nitro rotation, nitro out-of-plane deformation, and single proton transfer (SPT) along hydrogen bond are revealed. The results of potential energy profile calculations demonstrate that the SPT processes along the hydrogen bonds are unfavorable to occur statically, which is in great contrast to the natural base pair. The non-adiabatic dynamic simulations show that the excited-state nitro rotation and nitro out-of-plane deformation are the two important relaxation channels which lead to the fast internal conversion to S₀ state. The SPT from Z to P is also observed, followed by distortion on P, inducing the fast internal conversion to S₀ state. However, this channel (decay via SPT process) plays minor roles on the excited-state relaxation mechanism statistically. This work shows the great differences of the excited-state relaxation mechanism between the natural base pairs and artificial base pair, also sheds new light into the role of hydrogen bond and nitro group in P-Z base pair.

Theoretical investigation of conformational deviation of the human parallel telomeric G-quadruplex DNA in the presence of different salt concentrations and temperatures under confinement

Various experimental reports address the stability of G-quadruplex DNA inside a close confinement such as α-hemolysin, nanocavity water pool and different metal-organic-frameworks (MOFs). To understand the conformational change of G-quadruplex DNA at the atomistic level, we have carried out a total of 40 μs simulation run under both non-polar and polar confinement conditions. To investigate the dynamics, we have considered two different KCl salt concentrations, i.e., 0.47 M (minimal salt concentration) and higher than 2 M (higher salt concentration), at two distinct temperatures, 300 K and 350 K. Here, we have observed that the human telomeric G-quadruplex DNA deviates more from its crystal structure at minimal salt concentration under both non-polar and polar confinement conditions. Besides, the loop regions deviate and fluctuate more compared to the other regions, i.e., sugar-phosphate backbone and tetrad regions. The presence of K+ ions is found to be primarily responsible for this phenomenon. From the spatial density function (SDF) plots, a higher density of K+ ions is observed in the backbone region. Furthermore, from the residue-wise first solvation shell estimation, we have noticed that the K+ ions mainly accumulate in the tetrad region under both non-polar and polar confinement conditions due to which the tetrad regions are more rigid than the loop regions. Higher salt concentration results in increased rigidity of the G-quadruplex DNA. Our study provides valuable insight into the conformational deviation of the G-quadruplex DNA under nanoconfinement conditions.

An in silico investigation of the binding modes and pathway of APTO-253 on c-KIT G-quadruplex DNA

The stability of c-KIT G-quadruplex DNA via ligands has been a significant concern in the growing field of cancer therapy. Thus, it is very important to understand the mechanism behind the high binding affinity of the small drug molecules on the c-KIT G-quadruplex DNA. In this study, we have investigated the binding mode and pathway of the APTO-253 ligand on the c-KIT G-quadruplex DNA employing a total of 10 μs all atom molecular dynamics simulations and further 8.82 μs simulations via the umbrella sampling method using both OL15 and BSC1 latest force fields for DNA structures. From the cluster structure analysis, mainly three binding pathways i.e., top, bottom and side loop stacking modes are identified. Moreover, RMSD, RMSF and 2D-RMSD values indicate that the c-KIT G-quadruplex DNA and APTO-253 molecules are stable throughout the simulation run. Furthermore, the number of hydrogen bonds in each tetrad and the distance between the two central K⁺ cations confirm that the c-KIT G-quadruplex DNA maintains its conformation in the process of complex formation with the APTO-253 ligand. The binding free energies and the minimum values in the potential of mean forces suggest that the binding processes are energetically favorable. Furthermore, we have found that the bottom stacking mode is the most favorable binding mode among all the three modes for the OL15 force field. However, for the BSC1 force field, both the top and bottom binding modes of the APTO-253 ligand in c-KIT G-quadruplex DNA are comparable to each other. To investigate the driving force for the complex formation, we have noticed that the van der Waals (vdW) and π-π stacking interactions are mainly responsible. Our detailed studies provide useful information for the discovery of novel drugs in the field of stabilization of G-quadruplex DNAs.

Understanding The Role of Reline, a Natural DES, on Temperature-Induced Conformational Changes of C-Kit G-Quadruplex DNA: A Molecular Dynamics Study

The noncanonical guanine-rich DNAs have drawn particular attention to the scientific world due to their controllable diverse and polymorphic structures. Apart from biological and medical significance, G-quadruplex DNAs are widely used in various fields such as nanotechnology, nanomachine, biosensors, and biocatalyst. So far, the applications of the G-quadruplex DNA are mainly limited in the water medium. Recently, a new generation of solvent named deep eutectic solvent (DES) has become very popular and has been widely used as a reaction medium of biocatalytic reactions and long-term storage medium for nucleic acids, even at high temperature. Hence, it is essential to understand the role of DES on temperature-induced conformational changes of a G-quadruplex DNA. In this research work, we have explored the temperature-mediated conformational dynamics of c-kit oncogene promoter G-quadruplex DNA in reline medium in the temperature range of 300-500 K, using a total of 10 μs unbiased all-atom molecular dynamics simulation. Here, from RMSD, RMSF, R_g and principal component analyses, we notice that the c-kit G-quadruplex DNA is stable up to 450 K in reline medium. However, it unfolds in water medium at 450 K. It is found that the hydrogen bonding interactions between c-kit G-quadruplex DNA and reline play a key role in the stabilization of the G-quadruplex DNA even at high temperature. Furthermore, in this work we have observed a very interesting and distinctive phenomenon of the central cation of the G-quadruplex DNA. Its position was seen to fluctuate between the two tetrad cores, that is, the region between tetrad-1 and tetrad-2 and that between tetrad-2 and tetrad-3 and vice versa at 450 and 500 K in reline medium which is absent in water medium at 450 K. Moreover, the rate of its oscillation is increased when temperature is increased.

Building better enzymes: Molecular basis of improved non-natural nucleobase incorporation by an evolved DNA polymerase

Obtaining semisynthetic microorganisms that exploit the information density of "hachimoji" DNA requires access to engineered DNA polymerases. A KlenTaq variant has been reported that incorporates the "hachimoji" P:Z nucleobase pair with a similar efficiency to that seen for Watson-Crick nucleobase incorporation by the wild type (WT) KlenTaq DNA polymerase. The variant polymerase differs from WT KlenTaq by only four amino acid substitutions, none of which are located within the active site. We now report molecular dynamics (MD) simulations on a series of binary complexes aimed at elucidating the contributions of the four amino acid substitutions to altered catalytic activity. These simulations suggest that WT KlenTaq is insufficiently flexible to be able to bind AEGIS DNA correctly, leading to the loss of key protein/DNA interactions needed to position the binary complex for efficient incorporation of the "hachimoji" Z nucleobase. In addition, we test literature hypotheses about the functional roles of each amino acid substitution and provide a molecular description of how individual residue changes contribute to the improved activity of the KlenTaq variant. We demonstrate that MD simulations have a clear role to play in systematically screening DNA polymerase variants capable of incorporating different types of nonnatural nucleobases thereby limiting the number that need to be characterized by experiment. It is now possible to build DNA molecules containing nonnatural nucleobase pairs in addition to A:T and G:C. Exploiting this development in synthetic biology requires engineered DNA polymerases that can replicate nonnatural nucleobase pairs. Computational studies on a DNA polymerase variant reveal how amino acid substitutions outside of the active site yield an enzyme that replicates nonnatural nucleobase pairs with high efficiency. This work will facilitate efforts to obtain bacteria possessing an expanded genetic alphabet.

Snapshots of an evolved DNA polymerase pre- and post-incorporation of an unnatural nucleotide

The next challenge in synthetic biology is to be able to replicate synthetic nucleic acid sequences efficiently. The synthetic pair, 2-amino-8-(1-beta-d-2′- deoxyribofuranosyl) imidazo [1,2-a]-1,3,5-triazin-[8H]-4-one (trivially designated P) with 6-amino-3-(2′-deoxyribofuranosyl)-5-nitro-1H-pyridin-2-one (trivially designated Z), is replicated by certain Family A polymerases, albeit with lower efficiency. Through directed evolution, we identified a variant KlenTaq polymerase (M444V, P527A, D551E, E832V) that incorporates dZTP opposite P more efficiently than the wild-type enzyme. Here, we report two crystal structures of this variant KlenTaq, a post-incorporation complex that includes a template-primer with P:Z trapped in the active site (binary complex) and a pre-incorporation complex with dZTP paired to template P in the active site (ternary complex). In forming the ternary complex, the fingers domain exhibits a larger closure angle than in natural complexes but engages the template-primer and incoming dNTP through similar interactions. In the binary complex, although many of the interactions found in the natural complexes are retained, there is increased relative motion of the thumb domain. Collectively, our analyses suggest that it is the post-incorporation complex for unnatural substrates that presents a challenge to the natural enzyme and that more efficient replication of P:Z pairs requires a more flexible polymerase.

Aptamer chemistry

Aptamers are single-stranded DNA or RNA molecules capable of tightly binding to specific targets. These functional nucleic acids are obtained by an in vitro Darwinian evolution method coined SELEX (Systematic Evolution of Ligands by EXponential enrichment). Compared to their proteinaceous counterparts, aptamers offer a number of advantages including a low immunogenicity, a relative ease of large-scale synthesis at affordable costs with little or no batch-to-batch variation, physical stability, and facile chemical modification. These alluring properties have propelled aptamers into the forefront of numerous practical applications such as the development of therapeutic and diagnostic agents as well as the construction of biosensing platforms. However, commercial success of aptamers still proceeds at a weak pace. The main factors responsible for this delay are the susceptibility of aptamers to degradation by nucleases, their rapid renal filtration, suboptimal thermal stability, and the lack of functional group diversity. Here, we describe the different chemical methods available to mitigate these shortcomings. Particularly, we describe the chemical post-SELEX processing of aptamers to include functional groups as well as the inclusion of modified nucleoside triphosphates into the SELEX protocol. These methods will be illustrated with successful examples of chemically modified aptamers used as drug delivery systems, in therapeutic applications, and as biosensing devices.