Optimization of unnatural base pair packing for polymerase recognition

As part of an effort to expand the genetic alphabet, we have been examining the ability of predominately hydrophobic nucleobase analogues to pair in duplex DNA and during polymerase-mediated replication. We previously reported the synthesis and thermal stability of unnatural base pairs formed between nucleotides bearing simple methyl-substituted phenyl ring nucleobase analogues. Several of these pairs are virtually as stable and selective as natural base pairs in the same sequence context. Here, we report the characterization of polymerase-mediated replication of the same unnatural base pairs. We find that every facet of replication, including correct and incorrect base pair synthesis, as well as continued primer extension beyond the unnatural base pair, is sensitive to the specific methyl substitution pattern of the nucleobase analogue. The results demonstrate that neither hydrogen bonding nor large aromatic surface area is required for polymerase recognition, and that interstrand interactions between small aromatic rings may be optimized for replication. Combined with our previous results, these studies suggest that appropriately derivatized phenyl nucleobase analogues represent a promising approach toward developing a third base pair and expanding the genetic alphabet.     

Major groove substituents and polymerase recognition of a class of predominantly hydrophobic unnatural base pairs

Expansion of the genetic alphabet with an unnatural base pair is a long‐standing goal of synthetic biology. We have developed a class of unnatural base pairs, formed between d 5SICS and analogues of d MMO2 that are efficiently and selectively replicated by the Klenow fragment (Kf) DNA polymerase. In an effort to further characterize and optimize replication, we report the synthesis of five new d MMO2 analogues bearing different substituents designed to be oriented into the developing major groove and an analysis of their insertion opposite d 5SICS by Kf and Thermus aquaticus DNA polymerase I (Taq). We also expand the analysis of the previously optimized pair, d NaM –d 5SICS , to include replication by Taq. Finally, the efficiency and fidelity of PCR amplification of the base pairs by Taq or Deep Vent polymerases was examined. The resulting structure–activity relationship data suggest that the major determinants of efficient replication are the minimization of desolvation effects and the introduction of favorable hydrophobic packing, and that Taq is more sensitive than Kf to structural changes. In addition, we identify an analogue (d NMO1 ) that is a better partner for d 5SICS than any of the previously identified d MMO2 analogues with the exception of d NaM . We also found that d NaM –d 5SICS is replicated by both Kf and Taq with rates approaching those of a natural base pair.     

An efficiently extended class of unnatural base pairs

A third DNA base pair, which is synthesized efficiently and selectively, would have wide ranging applications from synthetic organisms to nucleic acids biotechnology. Hydrophobic unnatural nucleobases offer a promising route to such a pair, but are often limited by inefficient extension, defined as synthesis immediately following the unnatural pair. Here, we describe a simple screen which enables the characterization of large numbers of previously uncharacterized hetero base pairs. From this screen, we identified a class of unnatural base pairs which are extended more efficiently than any unnatural base pair reported to date. Screening, when complemented by further kinetic analysis, can improve the understanding of the determinants of efficient extension as well as identify viable hetero base pairs.     

Minor groove hydrogen bonds and the replication of unnatural base pairs

As part of an effort to expand the genetic alphabet, we examined the synthesis of DNA with six different unnatural nucleotides bearing methoxy-derivatized nucleobase analogues. Different nucleobase substitution patterns were used to systematically alter the nucleobase electronics, sterics, and hydrogen-bonding potential. We determined the ability of the Klenow fragment of E. coli DNA polymerase I to synthesize and extend the different unnatural base pairs and mispairs under steady-state conditions. Unlike other hydrogen-bond acceptors examined in the past, the methoxy groups do not facilitate mispairing, implying that they are not recognized by any of the hydrogen-bond donors of the natural nucleobases; however, they do facilitate replication. The more efficient replication results largely from an increase in the rate of extension of primers terminating at the unnatural base pair and, interestingly, requires that the methoxy group be at the ortho position where it is positioned in the developing minor groove and can form a functionally important hydrogen bond with the polymerase. Thus, ortho methoxy groups should be generally useful for the effort to expand the genetic alphabet.     

Configurational specificity of stacking interactions in DNA base pairs: A computational analysis

A theoretical study of stacking patterns of various hydrogen-bonded base pair complexes has been undertaken. Modified Rayleigh-Schrodinger perturbation theory for intermediate range interactions, has been employed to evaluate the stacking interactions using multicentered-multipole expansion method. Net atomic charge and corresponding dipole components located at each of the atomic centers have been computed by CNDO/2 method. An analysis of the intermolecular forces involved in the stable formation of the various base pair complexes, has been presented and the results have been discussed in the light of experimental as well as other theoretical observations. The possibility of relative preference of the left-handed configuration for alternating sequences has been quantitatively explored.     

Preopening of the DNA base pairs

The lower‐energy stable structures of the A?T base pair are revealed under a search of its potential energy surface in the vicinity of its Watson–Crick configuration performed at the PM3 computational level. Their properties and the mutual position of the nucleic acid bases A and T in these structures allow to partition them into three classes: partially preopened, stretched, and fully preopened. The preferable monohydration sites of the preopened, stretched, and fully preopened pairs are also determined. It is demonstrated, first, that the monohydration of the A?T pair at particular sites favors a base pair preopeness and, second, that a binding of the water molecule to the preopened A?T base pair on the major groove side enhances its stabilization. It is also shown that water molecule placing in the vicinity of the central H bond of the A?T pair significantly facilitates its preopening. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 82: 193–204, 2001     

Fluorescent hydrophobic zippers inside duplex DNA: interstrand stacking of perylene-3,4:9,10-tetracarboxylic acid bisimides as artificial DNA base dyes

Perylene-3,4:9,10-tetracarboxylic acid bisimides (PBs) were incorporated synthetically into oligonucleotides by using automated DNA building-block chemistry. The 2'-deoxyribofuranoside of the natural nucleosides was replaced by (S)-aminopropan-2,3-diol as an acyclic linker between the phosphodiester bridges that is tethered to one of the imide nitrogen atoms of the PB dye. The S configuration of this linker was chosen to mimic the stereochemical situation at the 3'-position of the natural 2'-deoxyribofuranosides. By using this strategy, up to six PB dyes were incorporated in the middle of 18-mer DNA duplexes by using interstrand alternating sequences of PBs with thymines or an abasic site analogue. Both PB dimers and PB hexamers as artificial base substitutions inside the duplexes yield characteristic excimer-type fluorescence. The stacking properties of the PB chromophores are modulated by the presence or absence of thymines opposite the PB modification site in the counterstrand. The interstrand PB dimers can be regarded as hydrophobically interacting base pairs, which display a characteristic fluorescence readout signal. Hence, for the PB hexamers, we proposed a zipperlike recognition motif that is formed inside duplex DNA. The PB zipper shows characteristic excimer-type emission as a fluorescence readout signal for the pairing interaction.     

Halogen bonding in DNA base pairs

Halogen bonding (R-X···Y) is a qualitative analogue of hydrogen bonding that may prove useful in the rational design of artificial proteins and nucleotides. We explore halogen-bonded DNA base pairs containing modified guanine, cytosine, adenine and thymine nucleosides. The structures and stabilities of the halogenated systems are compared to the normal hydrogen bonded base pairs. In most cases, energetically stable, coplanar structures are identified. In the most favorable cases, halogenated base pair stabilities are within 2 kcal mol(-1) of the hydrogen bonded analogues. Among the halogens X = Cl, Br, and I, bromine is best suited for inclusion in these biological systems because it possesses the best combination of polarizability and steric suitability. We find that the most stable structures result from a single substitution of a hydrogen bond for a halogen bond in dA:dT and dG:dC base pairs, which allows 1 or 2 hydrogen bonds, respectively, to complement the halogen bond.     

Electronic properties of metal-modified DNA base pairs

The electronic properties of several metal-modified Watson-Crick guanine-cytosine base pairs are investigated by means of first-principle density functional theory calculations. Focus is placed on a new structure recently proposed as a plausible model for building an antiparallel duplex with Zn-guanine-cytosine pairs, but we also inspect several other conformations and the incorporation of Ag and Cu ions. We analyze the effects induced by the incorporation of one metal cation per base pair by comparing the structures and the electronic properties of the metalated pairs to those of the natural guanine-cytosine pair, particularly for what concerns the modifications of energy levels and charge density distributions of the frontier orbitals. Our results reveal the establishment of covalent bonding between the metal cation and the nucleobases, identified in the presence of hybrid metal-guanine and metal-cytosine orbitals. Attachment of the cation can occur either at the N1 or the N7 site of guanine and is compatible with altering or not altering the H-bond pattern of the natural pair. Cu(II) strongly contributes to the hybridization of the orbitals around the band gap, whereas Ag(I) and Zn(II) give hybrid states farther from the band gap. Most metalated pairs have smaller band gaps than the natural guanine-cytosine pair. The band gap shrinking along with the metal-base coupling suggests interesting consequences for electron transfer through DNA double helices.     

Symmetry selection in artificial DNA base pairs

We report the results of density functional theory (DFT) studies on the formation of the complex H1--Cu2+-H1- consisting of two deprotonated hydroxypyridone ligands (H1-) and a Cu2+ ion. We compare the total energies of three possible structures with different symmetries and show that the structure with plane reflection symmetry has the lowest energy. The electronic structure of the periodic extended DNA-like double helix consisting of stacked H1--Cu2+-H1- units is then calculated within the density functional method, and the double helix is found to be an insulating ferromagnet.     

Chimeric GNA/DNA metal-mediated base pairs

DNA double helices comprising chimeric GNA/DNA metal-mediated base pairs have been synthesized and characterized (GNA = glycol nucleic acid). The possibility to combine different nucleic acid backbones within one metal-mediated base pair expands the applicability of metal-functionalized nucleic acids.     

Efforts toward expansion of the genetic alphabet: structure and replication of unnatural base pairs

Expansion of the genetic alphabet has been a long-time goal of chemical biology. A third DNA base pair that is stable and replicable would have a great number of practical applications and would also lay the foundation for a semisynthetic organism. We have reported that DNA base pairs formed between deoxyribonucleotides with large aromatic, predominantly hydrophobic nucleobase analogues, such as propynylisocarbostyril (dPICS), are stable and efficiently synthesized by DNA polymerases. However, once incorporated into the primer, these analogues inhibit continued primer elongation. More recently, we have found that DNA base pairs formed between nucleobase analogues that have minimal aromatic surface area in addition to little or no hydrogen-bonding potential, such as 3-fluorobenzene (d3FB), are synthesized and extended by DNA polymerases with greatly increased efficiency. Here we show that the rate of synthesis and extension of the self-pair formed between two d3FB analogues is sufficient for in vitro DNA replication. To better understand the origins of efficient replication, we examined the structure of DNA duplexes containing either the d3FB or dPICS self-pairs. We find that the large aromatic rings of dPICS pair in an intercalative manner within duplex DNA, while the d3FB nucleobases interact in an edge-on manner, much closer in structure to natural base pairs. We also synthesized duplexes containing the 5-methyl-substituted derivatives of d3FB (d5Me3FB) paired opposite d3FB or the unsubstituted analogue (dBEN). In all, the data suggest that the structure, electrostatics, and dynamics can all contribute to the extension of unnatural primer termini. The results also help explain the replication properties of many previously examined unnatural base pairs and should help design unnatural base pairs that are better replicated.     

Dynamic behavior of DNA base pairs containing 8-oxoguanine

The process by which DNA repair enzymes recognize and selectively excise damaged bases in duplex DNA is fundamental to our mechanistic understanding of these critical biological reactions. 8-Oxoguanine (8-oxoG) is the most common form of oxidative DNA damage; unrepaired, this lesion generates a G:C-->T:A mutation. Central to the recognition and repair of DNA damage is base extrusion, a process in which the damaged base lesion or, in some cases, its partner disengages from the helix and is bound to the enzyme's active site where base excision takes place. The conformation adopted by 8-oxoG in duplex DNA is affected by the base positioned opposite this lesion; conformational changes may also take place when the damaged base binds to its cognate repair enzyme. We performed unrestrained molecular dynamics simulations for several 13-mer DNA duplexes. Oligomers containing G:C and 8oxoG:C pairs adopted Watson-Crick geometries in stable B-form duplexes; 8oxoG showed increased local and global flexibility and a reduced barrier to base extrusion. Duplexes containing the G:A mismatch showed much larger structural fluctuations and failed to adopt a well-defined structure. For the 8oxoG:A mismatch that is recognized by the DNA glycosylase MutY, the damaged nucleoside underwent spontaneous and reproducible anti-->syn transitions. The syn conformation is thermodynamically preferred. Steric hindrance and unfavorable electrostatics associated with the 8oxoG O8 atom in the anti conformation were the major driving forces for this transition. Transition events follow two qualitatively different pathways. The overall anti-->syn transition rate and relative probability of the two transition paths were dependent on local sequence context. These simulations indicate that both the dynamic and equilibrium behavior of the duplex change as a result of oxidation; these differences may provide valuable new insight into the selective action of enzymes on damaged DNA.