Non-natural nucleotides as probes for the mechanism and fidelity of DNA polymerases

DNA synthesis technologies to close the gene writing gap

Synthetic DNA is of increasing demand across many sectors of research and commercial activities. Engineering biology, therapy, data storage and nanotechnology are set for rapid developments if DNA can be provided at scale and low cost. Stimulated by successes in next generation sequencing and gene editing technologies, DNA synthesis is already a burgeoning industry. However, the synthesis of >200 bp sequences remains unaffordable. To overcome these limitations and start writing DNA as effectively as it is read, alternative technologies have been developed including molecular assembly and cloning methods, template-independent enzymatic synthesis, microarray and rolling circle amplification techniques. Here, we review the progress in developing and commercializing these technologies, which are exemplified by innovations from leading companies. We discuss pros and cons of each technology, the need for oversight and regulatory policies for DNA synthesis as a whole and give an overview of DNA synthesis business models.

Pd-Catalyzed Ring-Opening/Arylation/Cyclization of 2-Aminothiazole Derivatives Provides Modular Access to Isocytosine Analogues

We report a Pd-catalyzed ring-opening/arylation/cyclization reaction sequence between 2-aminothiazoles and aryl (pseudo)halides that provides modular access to isocytosine analogues. The scope of the reaction is broad with respect to both coupling partners and a robustness test demonstrated the functional group tolerance of the methodology. Visual kinetic analysis revealed that the product may inhibit catalyst turnover for some substrates.

Completing the β,γ-CXY-dNTP Stereochemical Probe Toolkit: Synthetic Access to the dCTP Diastereomers and 31P and 19F NMR Correlations with Absolute Configurations

Nucleoside 5'-triphosphate (dNTP) analogues in which the β,γ-oxygen is mimicked by a CXY group (β,γ-CXY-dNTPs) have provided information about DNA polymerase catalysis and fidelity. Definition of CXY stereochemistry is important to elucidate precise binding modes. We previously reported the (R)- and (S)-β,γ-CHX-dGTP diastereomers (X = F, Cl), prepared via P,C-dimorpholinamide CHCl (6a, 6b) and CHF (7a, 7b) bisphosphonates (BPs) equipped with an (R)-mandelic acid as a chiral auxiliary, with final deprotection using H₂/Pd. This method also affords the β,γ-CHCl-dTTP (11a, 11b), β,γ-CHF (12a, 12b), and β,γ-CHCl (13a, 13b) dATP diastereomers as documented here, but the reductive deprotection step is not compatible with dCTP or the bromo substituent in β,γ-CHBr-dNTP analogues. To complete assembly of the toolkit, we describe an alternative synthetic strategy featuring ethylbenzylamine or phenylglycine-derived chiral BP synthons incorporating a photolabile protecting group. After acid-catalyzed removal of the (R)-(+)-α-ethylbenzylamine auxiliary, coupling with activated dCMP and photochemical deprotection, the individual diastereomers of β,γ-CHBr- (33a, 33b), β,γ-CHCl- (34a, 34b), β,γ-CHF-dCTP (35a, 35b) were obtained. The β,γ-CH(CH₃)-dATPs (44a, 44b) were obtained using a methyl (R)-(-)-phenylglycinate auxiliary. ³¹P and ¹⁹F NMR Δδ values are correlated with CXY stereochemistry and pK_a2-4 values for 13 CXY-bisphosphonic acids and imidodiphosphonic acid are tabulated.

Chemical Modification of Aptamers for Increased Binding Affinity in Diagnostic Applications: Current Status and Future Prospects

Aptamers are short single stranded DNA or RNA oligonucleotides that can recognize analytes with extraordinary target selectivity and affinity. Despite their promising properties and diagnostic potential, the number of commercial applications remains scarce. In order to endow them with novel recognition motifs and enhanced properties, chemical modification of aptamers has been pursued. This review focuses on chemical modifications, aimed at increasing the binding affinity for the aptamer's target either in a non-covalent or covalent fashion, hereby improving their application potential in a diagnostic context. An overview of current methodologies will be given, thereby distinguishing between pre- and post-SELEX (Systematic Evolution of Ligands by Exponential Enrichment) modifications.

Base modifications affecting RNA polymerase and reverse transcriptase fidelity

Ribonucleic acid (RNA) is capable of hosting a variety of chemically diverse modifications, in both naturally-occurring post-transcriptional modifications and artificial chemical modifications used to expand the functionality of RNA. However, few studies have addressed how base modifications affect RNA polymerase and reverse transcriptase activity and fidelity. Here, we describe the fidelity of RNA synthesis and reverse transcription of modified ribonucleotides using an assay based on Pacific Biosciences Single Molecule Real-Time sequencing. Several modified bases, including methylated (m⁶A, m⁵C and m⁵U), hydroxymethylated (hm⁵U) and isomeric bases (pseudouridine), were examined. By comparing each modified base to the equivalent unmodified RNA base, we can determine how the modification affected cumulative RNA polymerase and reverse transcriptase fidelity. 5-hydroxymethyluridine and N⁶-methyladenosine both increased the combined error rate of T7 RNA polymerase and reverse transcriptases, while pseudouridine specifically increased the error rate of RNA synthesis by T7 RNA polymerase. In addition, we examined the frequency, mutational spectrum and sequence context of reverse transcription errors on DNA templates from an analysis of second strand DNA synthesis.

Structural and Energetic Impact of Non-natural 7-Deaza-8-azaguanine, 7-Deaza-8-azaisoguanine, and Their 7-Substituted Derivatives on Hydrogen-Bond Pairing with Cytosine and Isocytosine

The impact of 7-deaza-8-azaguanine (DAG) and 7-deaza-8-azaisoguanine (DAiG) modifications on the geometry and stability of the G:C Watson-Crick (cWW) base pair and the G:iC and iG:C reverse Watson-Crick (tWW) base pairs has been characterized theoretically. In addition, the effect on the same base pairs of seven C7-substituted DAG and DAiG derivatives, some of which have been previously experimentally characterized, has been investigated. Calculations indicate that all of these modifications have a negligible impact on the geometry of the above base pairs, and that modification of the heterocycle skeleton has a small impact on the base-pair interaction energies. Instead, base-pair interaction energies are dependent on the nature of the C7 substituent. For the 7-substituted DAG-C cWW systems, a linear correlation between the base-pair interaction energy and the Hammett constant of the 7-substituent is found, with higher interaction energies corresponding to more electron-withdrawing substituents. Therefore, the explored modifications are expected to be accommodated in both parallel and antiparallel nucleic acid duplexes without perturbing their geometry, while the strength of a base pair (and duplex) featuring a DAG modification can, in principle, be tuned by incorporating different substituents at the C7 position.

Rolling Circle Transcription-Amplified Hierarchically Structured Organic-Inorganic Hybrid RNA Flowers for Enzyme Immobilization

Programmable nucleic acids have emerged as powerful building blocks for the bottom-up fabrication of two- or three-dimensional nano- and microsized constructs. Here we describe the construction of organic-inorganic hybrid RNA flowers (hRNFs) via rolling circle transcription (RCT), an enzyme-catalyzed nucleic acid amplification reaction. These hRNFs are highly adaptive structures with controlled sizes, specific nucleic acid sequences, and a highly porous nature. We demonstrated that hRNFs are applicable as potential biological platforms, where the hRNF scaffold can be engineered for versatile surface functionalization and the inorganic component (magnesium ions) can serve as an enzyme cofactor. For surface functionalization, we proposed robust and straightforward approaches including in situ synthesis of functional hRNFs and postfunctionalization of hRNFs that enable facile conjugation with various biomolecules and nanomaterials (i.e., proteins, enzymes, organic dyes, inorganic nanoparticles) using selective chemistries (i.e., avidin-biotin interaction, copper-free click reaction). In particular, we showed that hRNFs can serve as soft scaffolds for β-galactosidase immobilization and greatly enhance enzymatic activity and stability. Therefore, the proposed concepts and methodologies are not only fundamentally interesting when designing RNA scaffolds or RNA bionanomaterials assembled with enzymes but also have significant implications on their future utilization in biomedical applications ranging from enzyme cascades to biosensing and drug delivery.

2'-Deoxyribonucleoside 5'-triphosphates bearing 4-phenyl and 4-pyrimidinyl imidazoles as DNA polymerase substrates

We developed a versatile access to a series of 4-substituted imidazole 2'-deoxynucleoside triphosphate bearing functionalized phenyl or pyrimidinyl rings. 4-Iodo-1H-imidazole was enzymatically converted into the corresponding 2'-deoxynucleoside, which was then chemically derived into its 5'-triphosphate, followed by 4-arylation via Suzuki-Miyaura coupling using (hetero)arylboronic acids. Both KF (exo-) and Deep Vent (exo-) DNA polymerases incorporated these modified nucleotides in primer-extension assays, adenine being the preferred pairing partner in the template. The 4-(3-aminophenyl)imidazole derivative (3APh) was the most efficiently inserted opposite A by KF (exo-) with only a 37-fold lower efficiency (Vmax/KM) than that of the correct dTTP. No further extension occurred after the incorporation of a single aryl-imidazole nucleotide. Interestingly, the aryl-imidazole dNTPs were found to undergo successive incorporation by calf thymus terminal deoxynucleotidyl transferase with different tailing efficiencies among this series and with a marked preference for 2APyr polymerization.

Mass-spectrometry analysis of modifications at DNA termini induced by DNA polymerases

Non-natural nucleotide substrates are widely used in the enzymatic synthesis of modified DNA. The terminal activity of polymerases in the presence of modified nucleotides is an important, but poorly characterized, aspect of enzymatic DNA synthesis. Here, we studied different types of polymerase activity at sequence ends using extendable and non-extendable synthetic models in the presence of the Cy5-dUTP analog Y. In primer extension reactions with selected exonuclease-deficient polymerases, nucleotide Y appeared to be a preferential substrate for non-templated 3'-tailing, as determined by MALDI mass-spectrometry and gel-electrophoresis. This result was further confirmed by the 3'-tailing of a non-extendable hairpin oligonucleotide model. Additionally, DNA polymerases induce an exchange of the 3' terminal thymidine for a non-natural nucleotide via pyrophosphorolysis in the presence of inorganic pyrophosphate. In primer extension reactions, the proofreading polymerases Vent, Pfu, and Phusion did not support the synthesis of Y-modified primer strand. Nevertheless, Pfu and Phusion polymerases were shown to initiate terminal nucleotide exchange at the template. Unlike non-proofreading polymerases, these two enzymes recruit 3'-5' exonuclease functions to cleave the 3' terminal thymidine in the absence of pyrophosphate.

The Toolbox for Modified Aptamers

Aptamers are nucleic acid-based scaffolds that can bind with high affinity to a variety of biological targets. Aptamers are identified from large DNA or RNA libraries through a process of directed molecular evolution (SELEX). Chemical modification of nucleic acids considerably increases the functional and structural diversity of aptamer libraries and substantially increases the affinity of the aptamers. Additionally, modified aptamers exhibit much greater resistance to biodegradation. The evolutionary selection of modified aptamers is conditioned by the possibility of the enzymatic synthesis and replication of non-natural nucleic acids. Wild-type or mutant polymerases and their non-natural nucleotide substrates that can support SELEX are highlighted in the present review. A focus is made on the efforts to find the most suitable type of nucleotide modifications and the engineering of new polymerases. Post-SELEX modification as a complementary method will be briefly considered as well.

C(α) torsion angles as a flexible criterion to extract secrets from a molecular dynamics simulation

Given the increasing complexity of simulated molecular systems, and the fact that simulation times have now reached milliseconds to seconds, immense amounts of data (in the gigabyte to terabyte range) are produced in current molecular dynamics simulations. Manual analysis of these data is a very time-consuming task, and important events that lead from one intermediate structure to another can become occluded in the noise resulting from random thermal fluctuations. To overcome these problems and facilitate a semi-automated data analysis, we introduce in this work a measure based on C(α) torsion angles: torsion angles formed by four consecutive C(α) atoms. This measure describes changes in the backbones of large systems on a residual length scale (i.e., a small number of residues at a time). Cluster analysis of individual C(α) torsion angles and its fuzzification led to continuous time patches representing (meta)stable conformations and to the identification of events acting as transitions between these conformations. The importance of a change in torsion angle to structural integrity is assessed by comparing this change to the average fluctuations in the same torsion angle over the complete simulation. Using this novel measure in combination with other measures such as the root mean square deviation (RMSD) and time series of distance measures, we performed an in-depth analysis of a simulation of the open form of DNA polymerase I. The times at which major conformational changes occur and the most important parts of the molecule and their interrelations were pinpointed in this analysis. The simultaneous determination of the time points and localizations of major events is a significant advantage of the new bottom-up approach presented here, as compared to many other (top-down) approaches in which only the similarity of the complete structure is analyzed.

Conformational changes of the phenyl and naphthyl isocyanate-DNA adducts during DNA replication and by minor groove binding molecules

DNA lesions produced by aromatic isocyanates have an extra bulky group on the nucleotide bases, with the capability of forming stacking interaction within a DNA helix. In this work, we investigated the conformation of the 2′-deoxyadenosine and 2′-deoxycytidine derivatives tethering a phenyl or naphthyl group, introduced in a DNA duplex. The chemical modification experiments using KMnO₄ and 1-cyclohexyl-3 -(2-morpholinoethyl) carbodiimide metho-p-toluenesulfonate have shown that the 2′-deoxycytidine lesions form the base pair with guanine while the 2′-deoxyadenosine lesions have less ability of forming the base pair with thymine in solution. Nevertheless, the kinetic analysis shows that these DNA lesions are compatible with DNA ligase and DNA polymerase reactions, as much as natural DNA bases. We suggest that the adduct lesions have a capability of adopting dual conformations, depending on the difference in their interaction energies between stacking of the attached aromatic group and base pairing through hydrogen bonds. It is also presented that the attached aromatic groups change their orientation by interacting with the minor groove binding netropsin, distamycin and synthetic polyamide. The nucleotide derivatives would be useful for enhancing the phenotypic diversity of DNA molecules and for exploring new non-natural nucleotides.

Tolerance of base pair size and shape in postlesion DNA synthesis

The influence of base pair size and shape on the fidelity of DNA polymerase-mediated extension past lesion-containing mispairs was examined. Primer extension analysis was performed with synthetic nucleosides paired opposite the pro-mutagenic DNA lesion O(6)-benzylguanine (O(6)-BnG). These data indicate that the error-prone DNA polymerase IV (Dpo4) inefficiently extended past the larger Peri:O(6)-BnG base pair, and in contrast, error-free extension was observed for the smaller BIM:O(6)-BnG base pair. Steady-state kinetic analysis revealed that Dpo4 catalytic efficiency was strongly influenced by the primer:template base pair. Compared to the C:G pair, a 1.9- and 79,000-fold reduction in Dpo4 efficiency was observed for terminal C:O(6)-BnG and BIM:G base pairs respectively. These results demonstrate the impact of geometrical size and shape on polymerase-mediated mispair extension.

New pyrosequencing method to analyze the function of the Klenow fragment (EXO-) for unnatural nucleic acids: pyrophosphorolysis and incorporation efficiency

The incorporation of deoxynucleoside triphosphates (dNTPs) catalyzed by polymerases is conventionally examined using gel electrophoresis autoradiography. Here, we studied an alternative method, pyrosequencing, to verify the incorporation of dNTPs containing unnatural nucleotides by polymerases. We found that the pyrosequencing method more rapidly and easily confirmed the incorporation of dNTPs than the conventional method, especially in the presence of low-efficiency dNTP polymerases. Furthermore, the method can detect the pyrophosphorolysis reaction just before the position of the unnatural nucleic acid, and the efficiency of incorporation just after it.

The steric hypothesis for DNA replication and fluorine hydrogen bonding revisited in light of structural data

In DNA, bases pair in a molecular interaction that is both highly predictable and exquisitely specific. Therefore researchers have generally believed that the insertion of the matching nucleotide opposite a template base by DNA polymerases (pols) required Watson-Crick (W-C) hydrogen bond formation. However pioneering work by Kool and co-workers using hydrophobic base analogs such as the thymine (T) isostere 2,4-difluorotoluene (F) showed that shape rather than H-bonding served as the primary source of specificity in DNA replication by certain pols. This steric hypothesis for DNA replication has gained popularity, perhaps discouraging further experimental studies to address potential limitations of this new idea. The idea that shape trumps H-bonding in terms of pol selectivity largely hinges on the belief that fluorine is a poor H-bond acceptor. However, the shape complementarity model was embraced in the absence of any detailed structural data for match (F:A) and mismatch pairs (F:G, F:C, F:T) in DNA duplexes or at active sites of pols. Although the F and T nucleosides are roughly isosteric, it is unclear whether F:A and T:A pairs exhibit similar geometries. If the F:A pair is devoid of H-bonding, it will be notably wider than a T:A pair. Because shape and size and H-bonding are intimately related, it may not be possible to separate these two properties. Thus the geometries of an isolated F:A pair in water may differ considerably from an F:A pair embedded in a stretch of duplex DNA, at the tight active site of an A-family replicative pol, or within the spacious active site of a Y-family translesion pol. The shape complementarity model may have more significance for pol accuracy than efficiency: this model appears to be most relevant for replicative pols that use specific residues to probe the identity of the nascent base pair from the minor groove side. However, researchers have not fully considered the importance of such interactions that include H-bonds compared with W-C H-bonds in terms of pol fidelity and the shape complementarity model. This Account revisits the steric hypothesis for DNA replication in light of recent structural data and discusses the role of fluorine as an H-bond acceptor. Over the last 5 years, crystal structures have emerged for nucleic acid duplexes with F paired opposite to natural bases or located at the active sites of DNA pols. These data permit a more nuanced understanding of the role of shape in DNA replication and the capacity of fluorine to form H-bonds. These studies and additional research involving RNA or other fluorine-containing nucleoside analogs within duplexes indicate that fluorine engages in H-bonding in many cases. Although T and F are isosteric at the nucleoside level, replacement of a natural base by F in pairs often changes their shapes and sizes, and dF in DNA behaves differently from rF in RNA. Similarly, the pairing geometries observed for F and T opposite dATP, dGTP, dTTP, or dCTP and their H-bonding patterns at the active site of a replicative pol differ considerably.

Synthetic nucleotides as probes of DNA polymerase specificity

The genetic code is continuously expanding with new nucleobases designed to suit specific research needs. These synthetic nucleotides are used to study DNA polymerase dynamics and specificity and may even inhibit DNA polymerase activity. The availability of an increasing chemical diversity of nucleotides allows questions of utilization by different DNA polymerases to be addressed. Much of the work in this area deals with the A family DNA polymerases, for example, Escherichia coli DNA polymerase I, which are DNA polymerases involved in replication and whose fidelity is relatively high, but more recent work includes other families of polymerases, including the Y family, whose members are known to be error prone. This paper focuses on the ability of DNA polymerases to utilize nonnatural nucleotides in DNA templates or as the incoming nucleoside triphosphates. Beyond the utility of nonnatural nucleotides as probes of DNA polymerase specificity, such entities can also provide insight into the functions of DNA polymerases when encountering DNA that is damaged by natural agents. Thus, synthetic nucleotides provide insight into how polymerases deal with nonnatural nucleotides as well as into the mutagenic potential of nonnatural nucleotides.

A computational proposal for the experimentally observed discriminatory behavior of hypoxanthine, a weak universal nucleobase

A computational model composed of six nucleobases was used to investigate why hypoxanthine does not yield duplexes of equal stability when paired opposite each of the natural DNA nucleobases. The magnitudes of all nearest-neighbor interactions in a DNA helix were calculated, including hydrogen-bonding, intra- and interstrand stacking interactions, as well as 1-3 intrastrand stacking interactions. Although the stacking interactions in DNA relevant arrangements are significant and account for at least one third of the total stabilization energy in our nucleobase complexes, the trends in the magnitude of the stacking interactions cannot explain the relative experimental melting temperatures previously reported in the literature. Furthermore, although the total hydrogen-bonding interactions explain why hypoxanthine preferentially pairs with cytosine, the experimental trend for the remaining nucleobases (A, T, G) is not explained. In fact, the calculated pairing preference of hypoxanthine matches that determined experimentally only when the sum of all types of nearest-neighbor interactions is considered. This finding highlights a strong correlation between the relative magnitude of the total nucleobase-nucleobase interactions and measured melting temperatures for DNA strands containing hypoxanthine despite the potential role of other factors (including hydration, temperature, sugar-phosphate backbone). By considering a large range of sequence combinations, we reveal that the binding preference of hypoxanthine is strongly dependent on the nucleobase sequence, which may explain the varied ability of hypoxanthine to universally bind to the natural bases. As a result, we propose that future work should closely examine the interplay between the dominant nucleobase-nucleobase interactions and the overall strand stability to fully understand how sequence context affects the universal binding properties of modified bases and to aid the design of new molecules with ambiguous pairing properties.

Nucleotide Analogues as Probes for DNA and RNA Polymerases

Nucleotide analogues represent a major class of anti-cancer and anti-viral drugs, and provide an extremely powerful tool for dissecting the mechanisms of DNA and RNA polymerases. While the basic assays themselves are relatively straight-forward, a key issue is to appropriately design the studies to answer the mechanistic question of interest. This article addresses the major issues involved in designing these studies, and some of the potential difficulties that arise in interpreting the data. Examples are given both of the type of analogues typically used, the experimental approaches with different polymerases, and issues with data interpretation.