Crystal Structures of Non-Natural Nucleobase Pairs in A- and B-DNA\xe2\x80\xa0
The extent to which synthetic biology can be used to expand genetic information systems compatible with natural enzymes and cells will depend on the extent to which multiple and contiguous non-natural nucleobase pairs fit within the standard double helical conformations of DNA. Toward this goal, two non-standard nucleobases (Z, 6-amino-5-nitro-2(1H)-pyridone and P, 2-amino-imidazo[1,2-a]-1,3,5-triazin-4(8H)one) were designed to form a Z:P pair with a standard \xe2\x80\x9cedge on\xe2\x80\x9d Watson-Crick geometry, but with rearranged hydrogen bond donor and acceptor groups. Here, we present the crystal structures of two self-complementary 16-mer oligonucleotides containing Z:P pairs. The first contained two consecutive Z:P nucleobase pairs and was found to crystallize within a host-guest complex in B-form. The second contained six consecutive Z:P pairs; it was found to crystallize as an A-form DNA duplex, although it can adopt B-form in solution as inferred from circular dichroism spectra. Although Z:P pairs have some structural properties that are similar to those of G:C pairs, unique features include stacking of the nitro group on Z with the adjacent heterocyclic nucleobase ring in A-DNA. In both B-and A-DNA, major groove widths associated with the Z:P pairs are approximately 1 \xc3\x85 wider than those of comparable G:C pairs potentially due to the presence of the nitro group in Z. Thus, our structural studies suggest that multiple and consecutive Z:P pairs are readily accommodated in DNA duplex structures recognized by natural polymerases, and therefore the GACTZP synthetic genetic system has the requisite properties to expand sequence space.
INTRODUCTION<!>Synthesis and purification of Z-P containing oligonucleotides<!>Crystallization of ZP containing and other oligonucleotides<!>Data collection, structure determination, and crystallographic refinement<!>Circular dichroism experiments<!>Crystallization and structure determinations of ZP-containing oligonucleotides<!>Characterization of the helical properties of ZP-containing oligonucleotides<!>Structural properties of ZP nucleobase pairs in A- and B-DNA
<p>Originally defined as simply the field seeking to create artificial life,1 "synthetic biology" became more narrowly construed in the 1970s to mean the use of emerging recombinant DNA technologies to rearrange and modify natural genetic pieces.2 Today, synthetic biologists are no longer constrained to exploiting or manipulating the natural nucleotides and amino acids found in modern terran life, which are merely products of prebiotic chemistry, adaptation, accident, and evolutionary "drift". Instead, several kinds of nucleotide analogues have been reported that form non-standard nucleobase pairs "orthogonal" to the standard T:A and C:G pairs when they are incorporated into DNA.3-14 Further, in many cases, these extra deoxyribonucleotides can direct the synthesis of non-standard RNA containing extra ribonucleotides thereby increasing the number of codons that, in turn, can encode proteins containing additional amino acids.15,16 These developments have led to a new view of synthetic biology, adumbrated by Eric Kool,9 which seeks to create the properties that we value in life (including reproduction, adaptation, and evolution) using molecular platforms that combine natural and non-natural synthetic components.</p><p>In some cases, however, molecular behaviors peculiar to the added synthetic genetic components can limit that combination. For example, Romesberg has shown that two nonstandard nucleobases, d5SCIS and dNaM,8,17,18 designed to pair by hydrophobic and geometric complementarity, can be maintained in duplex DNA both in vitro by PCR polymerases and in living Escherichia coli cells for as long as nine hours.19 However, a potentially undesirable property of d5SCIS:dNaM within duplex DNA is that it intrinsically does not pair in the designed "edge-on" geometry, but rather by inter-strand "stacking" of the two hydrophobic species, which introduces significant distortion within the duplex DNA backbone.20 Only when bound within the active site of an engineered DNA polymerase is the d5SCIS:dNaM pair forced to resemble a the geometrical properties of a Watson-Crick pair.21 One might therefore predict that inclusion of d5SCIS:dNaM pairs at multiple or consecutive sites would have negatively impact the overall DNA conformation, especially given the known toxicity of intercalative DNA binding agents, which distort the DNA structure and prevent protein binding with deleterious consequences for gene regulation and subsequent transcriptional events.</p><p>A class of artificially expanded genetic information systems (AEGIS) does not differ so greatly from those found in natural biology.13 Rather than pairing via hydrophobic interactions, these retain the hydrogen bonding that joins the standard A:T and G:C pairs. The pairing is "orthogonal" to the pairing of standard bases because the hydrogen bonding units are shuffled. For example, the small and 6-amino-5-nitro-2(1H)-pyridone heterocycle (trivially named Z) presents a donor-donor-acceptor pattern of hydrogen bonding units to a complementary large 2-amino-imidazo[1,2-a]-1,3,5-triazin-4(8H)one heterocycle (trivially named P), which presents an acceptor-acceptor-donor pattern of hydrogen bonding units (Figure 1).22 Further, Z brings to nucleic acids a functional moiety (the nitro group) that is not found in any natural encoded biopolymer, a functional group that may contribute to the ability of GACTZP DNA to display especially effective binding power following in vitro selection.23 Superficially, Z:P resembles the C-G nucleobase pair,24 but nevertheless it has been shown that DNA containing up to four consecutive Z:P units can be amplified successfully with high fidelity.24,25</p><p>Although X-ray crystallography has in the past been used to examine DNA containing single non-standard pairs,18,20,21,26 no data inform us about how the double-helical structure of DNA might be perturbed by the presence multiple non-natural nucleobase pairs, particularly DNA containing contiguous non-natural pairs. We therefore sought to obtain a series of X-ray structures for 16-bp oligonucleotide duplexes containing contiguous Z:P paired nucleotides. Our initial strategy was to employ a "host-guest" system in which the Z:P-containing oligonucleotides were co-crystallized with the N-terminal fragment of the Moloney murine leukemia virus reverse transcriptase (MMLV-RT).27-29 This system has most recently been used to obtain structures of double-helical DNA containing a thymine dimer30 and DNA bound to a Co(III)-substituted form of bleomycin.31</p><p>Here, we present the 1.8 Å crystal structure of a 16-bp, self-complementary oligonucleotide containing only two consecutive Z:P pairs in B-form DNA. For comparative analysis, we also determined the crystal structures at 1.78 Å and 1.68 Å of host-guest complexes including oligonucleotides with consecutive G:C and A:T pairs in place of the Z:P pairs, respectively. A 16-bp, self-complementary oligonucleotide containing six consecutive Z:P nucleobase pairs was found to crystallize in the absence of the host; its 1.98 Å crystal structure is reported here as determined by Br SAD phasing. These structural studies provide a basis for analyzing the structural properties of Z:P nucleobase pairs in the context of both A- and B-DNA and their resultant helical properties.</p><!><p>Standard phosphoramidites (Bz-dA, Ac-dC, dmf-dG, dT, and 5-Br-dU-CE) and controlled pore glass (CPG) having standard nucleosides were purchased from Glen Research (Sterling, VA). AEGIS phosphoramidites (dZ and dP) were obtained from Firebird Biomolecular Sciences LLC (Alachua, FL). All oligonucleotides containing Z and P (2P, 3/6ZP, 3/6ZP Br1 and 3/6ZP Br2, see Figure 1 and Table 1) were synthesized on an ABI 394 DNA Synthesizer following standard phosphoramidite chemistry, as previously reported.32 The CPGs carrying the synthetic oligonucleotides were treated with 1 M DBU in anhydrous acetonitrile (2.0 mL) at room temperature for 24 hours to remove the NPE group from the Z nucleobase. Then, the CPGs were filtered and dried. The CPGs carrying 2P and 3/6ZP (without 5-Br-dU) were treated with concentrated ammonium hydroxide at 55 °C for 16 hours, while the CPGs carrying 3/6ZP Br1 and 3/6ZP Br2 (with 5-Br-dU) were treated with concentrated ammonium hydroxide at room temperature for 24 hours. After removal of ammonium hydroxide, the oligonucleotides containing Z and P were purified on ion-exchange HPLC, and then desalted using Sep-Pac® Plus C18 cartridges (Waters). Fully standard oligonucleotides were purchased from Midland Certified Reagent Co. (Midland, Texas) in desalted form and used without further purification.</p><!><p>Sixteen base pair self-complementary oligonucleotides containing two or six consecutive Z:P nucleobase pairs were designed to be compatible with our host-guest crystallization. The host in this system is the N-terminal fragment of Moloney murine leukemia virus reverse transcriptase, which was purified as previously described,33 and the guests are the various oligonucleotides. Oligonucleotide sequences screened for crystallization are compiled in Table 1. The oligonucleotides were resuspended in buffer containing 10 mM HEPES pH 7.0 and 10 mM MgCl2 to give a final concentration of 2.5 mM duplex DNA and then annealed by heating to 70 °C followed by slow cooling to room temperature prior to crystallization. A 2.9 mM stock solution of the protein was diluted to 1.4 mM using 50 mM MES pH 6.0 and 0.3 M NaCl. This 1.4 mM sub-stock was then further diluted to 0.65 mM in 100 mM HEPES pH 7.5 and 0.3 M NaCl.</p><p>Each oligonucleotide was pre-complexed with the N-terminal fragment of Moloney murine leukemia virus reverse transcriptase with final concentrations of 0.71 and 0.46 mM, respectively, as previously described, and then subjected to self-nucleation or microseeding experiments. For microseeding, vapor diffusion hanging drops included 1 μL each of protein-DNA microseeds diluted in reservoir solution containing 7 % PEG 4000, 5 mM magnesium acetate and 50 mM N-(2-acetamido)iminodiacetic acid (ADA) at pH 6.5 and 1 μL of the protein-DNA complex solution, suspended over the reservoir solution. Using a similar strategy to that described above, self-complementary 16-mer duplex DNA oligonucleotides including either G:C or A:T pairs replacing the Z:P pairs found in the 3/6ZP oligonucleotide (Table 1) were crystallized as host-guest complexes.</p><p>The oligonucleotide including 6 consecutive Z:P pairs (3/6 Z:P) did not crystallize under our normal host-guest complex conditions and was then subjected to further screening of crystallization conditions, specifically the Natrix screening kit (Hampton Research, Inc), by using the ARI Gryphon crystallization robot. For this experiment, the protein and DNA concentrations were reduced to 0.23 mM and 0.35 mM, respectively. Promising crystals were obtained for this oligonucleotide; however, we noted a large amount of precipitant present in conditions that included 10 mM magnesium acetate, 50 mM MES pH 5.6, and 2.5 M ammonium sulfate. Subsequent crystallizations including just the oligonucleotide, diluted to 0.35 mM in the same buffer that was used to dilute the protein described above, produced crystals under similar conditions. Two different oligonucleotides including 5-bromouracil replacing Ts in the sequence were synthesized (Table 1) for phasing purposes. Both of these oligonucleotides crystallized under similar conditions that included 10 mM magnesium acetate, 50 mM MES pH 5.6, and 1.7-2.0 M ammonium sulfate. The crystal used for the structure determination was that of the Br1 sequence, which crystallized in 10 mM magnesium acetate, 50 mM MES pH 5.6, and 1.7 M ammonium sulfate and was cryo-cooled in a solution including the reservoir with 20% glycerol added.</p><!><p>A Br single wavelength anomalous dispersion experiment was performed to obtain phases for the 3/6 ZP-Br structure. Approximately 10-fold redundant data were collected to 1.98 Å at the APS SBC 19-ID beamline (Table 2) for a crystal of the Br1 sequence including two Br atoms from 5-bromouracil replacing thymines in the sequence as shown in Table 1. Initial Patterson searching and phasing calculations were done by using HKL3000.34 Specifically, a single Br site was identified in SHELXD35 and refined in SHELXE36. Initial phases were then calculated to 2.45 Å using MLPHARE, figure of merit 0.56, and improved and extended to 2.25 Å by solvent flattening using DM.37-39</p><p>This initial experimentally phased electron density map was of high quality and allowed automated building of the DNA structure by NAUTILUS,40 albeit with A:T pairs modeled in place of Z:P pairs. Following initial refinement of the structure generated by NAUTILUS in REFMAC5,41 positive peaks were observed in the Fo-Fc electron density maps consistent with missing chemical moieties present within the Z and P nucleotides. While maintaining the coordinate system, additional functional groups were added to the core NAUTILUS generated structure using the Molefacture Plug-in for Visual Molecular Dynamics (VMD) Graphics Viewer.42 Upon generation of the complete coordinates, each Z or P nucleobase was optimized for ten steps using GAMESS ab initio molecular quantum chemistry,43 with a 6-31G44 basis set. Each idealized nucleobase analog was then placed back into the model to form the completed DNA structure, including the novel functionality of the Z and P bases. Parameter files and linking statements were created for refinement in PHENIX.45 Addition of solvent molecules and manual model adjustment was done iteratively in COOT46 followed by refinement in PHENIX with maximum likelihood targets and isotropic B-factors. Addition of solvent molecules and manual model adjustment was done iteratively in COOT followed by refinement in PHENIX with maximum likelihood targets and isotropic B-factors.</p><p>Data were collected for the 2P host-guest complex crystal at the APS SBC 19-BM beamline to Bragg spacings of 1.8 Å and for the GC and AT host-guest complex crystals at the APS GM/CA 23-ID-D beamline to Bragg spacings of 1.78 Å and 1.7 Å, respectively (Table 2). The host-guest crystal structures were determined by molecular replacement as implemented in PHASER47 within the CCP4 suite of programs48 using the model of the N-terminal fragment of MMLV RT as the search model. This approach provides unbiased electron density for the DNA complexed to the protein. The protein model and associated water molecules were first adjusted in COOT and then refined initially in REFMAC and later in PHENIX to improve the electron density, and the DNA was subsequently modeled into Fo-Fc electron density maps. In building the DNA model, the first three nucleobase pairs were modeled and refined, then the next two pairs, and finally the remaining three pairs.</p><p>For the 2P structure, the Z:P nucleobase pairs were initially modeled as G:C pairs and then subsequently replaced by superimposing the base coordinates for either Z or P generated for the 3/6 ZP structure on the common atoms manually in COOT. The parameter and linking files were created for refinement in PHENIX. As the asymmetric unit includes one protein molecule and half of the DNA molecule, the DNA can be modeled as 8 pairs of duplex DNA or a single 16-mer strand of DNA. To ensure that the phosphodiester bond between bases 8 and 9 is appropriately connected, the DNA was modeled as a single 16-mer strand in the final round of the refinement, and the 8-mer duplex was then regenerated by symmetry.</p><p>Coordinates have been deposited for the 2P, GC, AT, and 3/6 ZP structures (Table 1) with PDB identifiers, 4XO0, 4XPE, 4XPC, and 4XNO, respectively.</p><!><p>The DNA sequences including 3/6 ZP (5'- CTTATPPPZZZATAAG -3'), AT sequence (5' CTTATAAATTTATAAG 3') and GC sequence (5' CTTATGGGCCCATAAG 3') were analyzed by circular dichroism to determine the helical form of the oligonucleotides buffered at neutral pH in low salt. For CD analysis, the 2.5 mM stocks of these DNA sequences were diluted to 5 μM with buffer containing 10 mM HEPES pH 7.0 and 10 mM MgCl2. The CD spectra for DNA sequences were collected on a Jasco J-810 CD instrument at 25 °C, at a rate of 50 nm/min and a wavelength increment of 0.1 nm. Ellipticity, Ø (mdegrees) was recorded for the DNA sequences from a wavelength of 320 to 220 nm. The spectra were the average of five scans corrected for ellipticity readings obtained for buffer (10 mM HEPES pH 7.0, 10 mM MgCl2) by itself.</p><!><p>Expansion of genetic alphabets that combine standard and non-natural nucleobase pairs in duplex DNA will be an essential element in the evolution of efforts in synthetic biology. The need to maintain (or at least prevent radical divergence from) existing canonical structures, however, represents an important constraint upon the incorporation of synthetic nucleotides into DNA if existing nucleic acid binding proteins or enzymes are to act on them. Of course, even canonical structures have considerable diversity, and the special properties of stacked A:T versus G:C nucleobase pairs in duplex DNA are well documented.49 Thus, characterizing the structural characteristics of DNA containing multiple and/or contiguous non-natural nucleobase pairs as compared to natural nucleobase pairs has considerable importance.</p><p>The focus of this study is the structural properties of stacked Z:P nucleobase pairs (Figure 1), which provide novel hydrogen-bonding arrangements between novel purine and pyrimidine-like bases. Accordingly, 16-mer self-complementary oligonucleotides containing either 4 or 6 Z:P pairs in total, with 2 or 6 consecutive Z:P pairs, respectively, were subjected to crystallization trials (Table 1). Initial crystallization screening involved the use of the host-guest system that the Georgiadis laboratory developed for the crystallization and analysis of novel DNA structures and DNA-ligand complexes. In this system, the N-terminal fragment of Moloney murine leukemia virus reverse transcriptase (MMLV-RT) serves as the host and a self-complementary 16-mer duplex DNA oligonucleotide as the guest. The host interacts with the three terminal nucleobase pairs on either end through minor groove hydrogen-bonding while the central 10 nucleobase pairs are free of interactions with either protein or other DNA molecules. The crystal lattice results from packing of complexes that include two protein molecules and one 16-mer oligonucleotide, while the asymmetric unit, the unique repeating unit in the crystal lattice, includes only one protein molecule and 8 nucleobase pairs of DNA.</p><p>The oligonucleotide including two consecutive Z:P nucleobase pairs (2P) but not the oligonucleotide including six contiguous Z:P pairs (3/6 ZP) crystallized in the host-guest system. As compared to other oligonucleotides crystallized as host-guest complexes, the 2P host-guest complex produced much smaller crystals that took longer to grow. The 2P host-guest crystal was phased by molecular replacement using the protein structure as the search model and refined to 1.8 Å resolution. The 3/6 ZP oligonucleotide was subjected to a high-throughput crystallization screen as a complex with the N-terminal fragment of MMLV-RT. Crystals obtained for the oligonucleotide were found to grow without the host protein under conditions including relatively high salt concentrations (~ 2 M ammonium sulfate). As no structural model was available for molecular replacement phasing of the 3/6 ZP structure, oligonucleotides including 5-bromouracil in place of thymine were synthesized and crystallized (Table 1) for experimental phasing purposes. The crystals obtained for the brominated oligonucleotides were actually larger and diffracted to higher resolution than the non-brominated crystals. A Br SAD phasing experiment was performed for the 3/6 ZP Br1 crystals producing a 2.25 Å experimental electron density map of excellent quality as shown in Figure 2, and the structure was refined to 1.98 Å.</p><p>The host-guest system readily accommodates B-form duplex DNA oligonucleotides that are 16 bps in length. It was therefore of interest to determine whether oligonucleotides analogous to 3/6 ZP containing either AT or GC pairs in place of the ZP pairs would crystallize as host-guest complexes (Table 1). Both the AT and GC oligonucleotides crystallized as host-guest complexes; the GC crystals were much smaller and took longer to grow than the AT oligonucleotides, which grew rapidly and to large size (0.2-0.3 mm in each dimension). Both the AT and GC structures were determined to 1.68 Å and 1.78 Å, respectively, by molecular replacement. Final 2Fo-Fc electron density maps are shown in Figure 3 for the DNA in each structure.</p><!><p>The 3/6 ZP oligonucleotide crystallized as A-form DNA in a lattice involving molecular packing interactions typically observed in other A-DNA structures.50 In this case, one end of the duplex forms hydrogen bonds in the minor groove with Z:P pairs located in the middle of another DNA duplex. Specifically, O2 of C1 hydrogen bonds to N2 of P7 and N2 of G16 hydrogen bonds to O2 of Z10 (see Figure 2 for numbering scheme). As analyzed by the software package 3DNA,51 the helical form of this structure is classified as A-form excluding the first three dinucleotide steps and the 13th dinucleotide step, which includes the 5-bromouracil. All six contiguous Z:P-containing dinucleotide steps are therefore classified as A-DNA.</p><p>Other oligonucleotides known to crystallize under high salt conditions in A-form are GC-rich and involve similar hydrogen-bonding interactions between the ends and G:C pairs in the middle of the oligonucleotides. The fact that 3/6 ZP contains 8 A:T pairs and yet still crystallized as A-form suggests that Z:P pairs may have an even higher propensity to adopt A-form DNA than GC-rich sequences. In a survey of A-DNA structures, the oligonucleotide sequences were G:C-rich with most containing fewer than 2 A:T nucleobase pairs. In fact, the structure most similar to 3/6 ZP in terms of length is that of a GC-rich 14 base pair oligonucleotides (PDB ID: 4OKL), which includes only 2 central A:T nucleobase pairs.</p><p>Given that 3/6 ZP crystallized as A-DNA in a high salt condition, it was of interest to determine its helical form in solution under low salt conditions. The CD spectra of all three sequences (AT, GC and 3/6 ZP) exhibit right-handed B-like spectra as shown in Figure 4, with subtle differences in the peak positions and heights pertaining to differences in the primary sequence of DNA. The spectra for the A:T-rich sequence exhibits a negative peak at 248 nm and a positive long wavelength peak at about 279 nm, typical of right handed B-DNA. The GC-rich sequence on the other hand has a broad negative peak centered at approximately 245 nm. The positive peak for the GC-rich sequence shifts to 270 nm instead of 279 nm observed for the AT rich sequence. The 3/6 ZP sequence behaves more similarly to the G:C-rich sequence with its broad negative peak at 241 nm and a positive peak at 273 nm. However, the spectral features of 3/6 ZP differ slightly from the control GC sequence. The positive peak for the 3/6 ZP sequence never reaches zero ellipticity, and at wavelengths greater than 290 nm, it attains another shorter positive peak extending from 290-320 nm. Overall, 3/6 ZP does appear to exist as B-DNA in solution, with CD spectral features most similar to GC.</p><p>In contrast to the 3/6 ZP, 2P, AT, and GC oligonucleotides (Table 1) crystallized in the host-guest system and exhibit B-form throughout, excluding the terminal dinucleotide steps, as analyzed by 3DNA.51 The AT and GC oligonucleotides are analogous to 3/6 ZP in that the 5'-P3Z3 sequence is replaced with 5'-G3C3 or 5'-A3T3, maintaining the positioning of purine and pyrimidine-like bases within the oligonucleotides. In theory, given that the 3/6 ZP oligonucleotide like the AT and GC oligonucleotides appears to adopt B-form under low salt conditions in solution as assessed by CD, it should have been possible to obtain host-guest crystals. The fact that it does not crystallize as a host-guest complex may be consistent with altered dynamic or other properties associate with the Z:P pairs.</p><!><p>Of particular interest are the structural properties associated with the unique NO2 group in Z in the context of A- or B-DNA. This nitro group was introduced initially to manage the chemical properties of the system presenting a hydrogen bond donor-donor-acceptor pattern (from the major to the minor groove).22 Subsequently, the nitro group was found to confer binding potential on GACTZP libraries that appears to be absent in standard GACT libraries.23 As anticipated, N4, N3, and O2 of Z hydrogen bond to O6, N1, and N2 of P, respectively, with distances of 2.7-3.0 Å in 2P and 3/6 ZP structures.</p><p>In the A-DNA structure, the Z nucleobase, including its NO2 group, is planar; this appears to facilitate stacking interactions between the NO2 of one Z nucleobase and the pyrimidine (or purine) ring of the adjacent nucleobase (Figure 5A). One oxygen atom of the nitro group in these stacked planar Z nucleobases is within hydrogen-bonding distance (~2.7 Å) of N4 of the pyrimidine-like ring. In contrast, in B-DNA, the NO2 group does not stack with adjacent bases but is nearly planar in both Z nucleobases with similar hydrogen-bonding between a nitro oxygen and N4 to that observed in 3/6ZP (Figure 5B). Accordingly, the electron density for the Z-NO2 groups in B-DNA is not as well ordered as in A-DNA (Figure 2), consistent with the more constrained conformation observed in the stacking interactions in A-form DNA. Potentially, the favorable stacking interactions of the Z-NO2 groups contribute to the ability of the 16-mer oligonucleotide including only six Z:P pairs to crystallize as A-form DNA. The fact that Z is a C-glycoside with a carbon-carbon linkage to its deoxyribose sugar rather than nitrogen-carbon linkage found in natural bases may also be a contributing factor.</p><p>To characterize the DNA parameters of Z:P pairs in A-DNA, we analyzed the local base pair step parameters, local base pair helical parameters, and the groove widths for the central Z:P pairs in our structure and compared them to those of the GC-rich 4OKL structure (5'-CCCCGGTACCGGGG) using the program 3DNA. Average values were calculated for five Z:P or P:Z dinucleotide steps for our structure and 10 GC or CG dinucleotide steps from 4OKL. As previously noted, A:T pairs are found infrequently in A-DNA structures and thus are not available for comparison. As shown in Table 3, the average rise, roll, and twist values for Z:P dinucleotide steps are similar to the values observed for GC dinucleotide steps and the average parameters for A-DNA reported by Lu et al.51 However, the average slide value for Z:P pairs of −2.08 Å in A-DNA is slightly larger than observed in G:C pairs in 4OKL (−1.95 Å) and on average in A-DNA structures −1.53 +/− 0.34 Å. This finding is consistent with the preferential stacking of the NO2 group with the adjacent pyrimidine or purine ring.</p><p>Otherwise, local base pair helical parameters including H-rise, inclination, H-twist, and X-displacement are all within the average range observed in A-DNA structures and are similar to those observed for G:C pairs in 4OKL. Within the central Z:P pairs of our structure, the average major groove width of 18.92 Å is approximately 1 Å wider the groove width associated with the G:C pairs (18.0 Å) in the 4OKL structure. Widening of the major groove may result from the presence of the Z-NO2 in the major groove. The minor groove widths for both Z:P and G:C pairs in A-DNA are more similar than the major groove widths with average values of 16.52 and 16.92 Å, respectively.</p><p>Similarly, the properties of the Z:P nucleobase pairs in B-DNA were assessed and compared to A:T and G:C pairs at the same positions in the DNA sequence in our host-guest structures. There are four Z:P pairs within six nucleotides (5'-PPTAZZ) in the 16-mer sequence crystallized; however, only half of this sequence and its complementary sequence are unique as the other half is related by crystallographic symmetry. These sequences are referred to as B-ZP, B-GC, and B-AT in contrast to those in the A-DNA structures, which are A-ZP and A-GC. In this case, the three dinucleotide steps containing ZP pairs and the analogous dinucleotide steps from the AT and GC oligonucleotide structures were analyzed. The average local base pair step parameters and local base pair helical parameters for B-ZP all fall within the range observed for B-DNA structures. Average values for each of these parameters for B-ZP are more similar to those calculated for B-GC than for B-AT. As was true for A-ZP, the major groove width of B-ZP of 18.7 Å is significantly (~ 0.7 Å) wider than that of B-GC (17.97 Å). However, it is slightly narrower than B-AT, which has an average major groove width of 19.07 Å. The minor groove width for B-ZP of 12.67 Å is similar to that of B-GC (12.43 Å) and significantly wider than that of B-AT (9.67 Å). AT-rich sequences are known to exhibit deep, narrow minor grooves, which bind to a number of small molecules.</p><p>Local base pair parameters including shear, stretch, stagger, buckle, propeller, and opening were also analyzed using 3DNA for A-ZP pairs and found to be fairly uniform in value. Propeller values of −10.76 ° and −13.68 ° are the largest observed values and occur at the junctions between the PZ or ZP and natural nucleobase pairs. For B-ZP, the pair Z6P11 exhibits an unusually large buckle angle of −14.25 ° (Figure 5C), which is much larger than buckle angles of −1.4 and 5.38 ° exhibited for equivalent A:T or G:C base pairs, respectively. This large buckle angle for Z6P11 may contribute to the offset hydrogen-bonding interactions observed for this base pair. In the major groove, the Z:P pair provides novel major groove hydrogen bonding opportunities with four electronegative atoms, two provided by the NO2 group, as compared to three present in either A:T or G:C pairs (Figure 6A-C). In the minor groove, Z:P pairs present the same hydrogen bonding pattern of electronegative atoms as the GC base pair (Figure 6D-F), with all three nucleobase pairs including O2 and N3 atoms.</p><p>Thus, our structural analysis indicates that Z:P nucleobase pairs are accommodated in both A- and B-DNA and overall exhibit DNA parameters that are quite similar to those observed for G:C pairs in the same helical form. The most significant impact of Z:P pairs on the structure of the DNA duplexes is widening of the major groove in both A and B-DNA as compared to values calculated for similar GC regions. This finding suggests that the NO2 group which faces the major groove has a widening effect. We speculate that stacked Z:P nucleobase pairs have a higher propensity to form A-DNA than equivalent sequences with G:C pairs based on the inability to crystallize the 3/6 ZP sequence in our host-guest system, which has a strong selection for B-form oligonucleotides, and the fact that the AT-rich 3/6 ZP sequence still crystallized in A-form while other examples of A-DNA are more G:C-rich. However, the CD spectrum of the 3/6 ZP oligonucleotide in a low salt buffering solution has features consistent with B-DNA and is similar overall to that obtained for the equivalent GC sequence. In conclusion, our analysis suggests that Z and P nucleotides will be useful in expanding the genetic alphabet, as contiguous Z:P pairs are well tolerated in standard helical forms of DNA.</p>