Evidence for a new Z-type left-handed DNA helix: properties of Z(WC)-DNA

Extensive Bioinformatics Analyses Reveal a Phylogenetically Conserved Winged Helix (WH) Domain (Zτ) of Topoisomerase IIα, Elucidating Its Very High Affinity for Left-Handed Z-DNA and Suggesting Novel Putative Functions

The dynamic processes operating on genomic DNA, such as gene expression and cellular division, lead inexorably to topological challenges in the form of entanglements, catenanes, knots, "bubbles", R-loops, and other outcomes of supercoiling and helical disruption. The resolution of toxic topological stress is the function attributed to DNA topoisomerases. A prominent example is the negative supercoiling (nsc) trailing processive enzymes such as DNA and RNA polymerases. The multiple equilibrium states that nscDNA can adopt by redistribution of helical twist and writhe include the left-handed double-helical conformation known as Z-DNA. Thirty years ago, one of our labs isolated a protein from Drosophila cells and embryos with a 100-fold greater affinity for Z-DNA than for B-DNA, and identified it as topoisomerase II (gene Top2, orthologous to the human UniProt proteins TOP2A and TOP2B). GTP increased the affinity and selectivity for Z-DNA even further and also led to inhibition of the isomerase enzymatic activity. An allosteric mechanism was proposed, in which topoII acts as a Z-DNA-binding protein (ZBP) to stabilize given states of topological (sub)domains and associated multiprotein complexes. We have now explored this possibility by comprehensive bioinformatic analyses of the available protein sequences of topoII representing organisms covering the whole tree of life. Multiple alignment of these sequences revealed an extremely high level of evolutionary conservation, including a winged-helix protein segment, here denoted as Zτ, constituting the putative structural homolog of Zα, the canonical Z-DNA/Z-RNA binding domain previously identified in the interferon-inducible RNA Adenosine-to-Inosine-editing deaminase, ADAR1p150. In contrast to Zα, which is separate from the protein segment responsible for catalysis, Zτ encompasses the active site tyrosine of topoII; a GTP-binding site and a GxxG sequence motif are in close proximity. Quantitative Zτ-Zα similarity comparisons and molecular docking with interaction scoring further supported the "B-Z-topoII hypothesis" and has led to an expanded mechanism for topoII function incorporating the recognition of Z-DNA segments ("Z-flipons") as an inherent and essential element. We further propose that the two Zτ domains of the topoII homodimer exhibit a single-turnover "conformase" activity on given G(ate) B-DNA segments ("Z-flipins"), inducing their transition to the left-handed Z-conformation. Inasmuch as the topoII-Z-DNA complexes are isomerase inactive, we infer that they fulfill important structural roles in key processes such as mitosis. Topoisomerases are preeminent targets of anti-cancer drug discovery, and we anticipate that detailed elucidation of their structural-functional interactions with Z-DNA and GTP will facilitate the design of novel, more potent and selective anti-cancer chemotherapeutic agents.

The Origin of Left-Handed Poly[d(G-C)]

The discovery of a reversible transition in the helical sense of a double-helical DNA was initiated by the first synthesis in 1967 of the alternating sequence poly[d(G-C)]. In 1968, exposure to high salt concentration led to a cooperative isomerization of the double helix manifested by an inversion in the CD spectrum in the 240-310 nm range and in an altered absorption spectrum. The tentative interpretation, reported in 1970 and then in detailed form in a 1972 publication by Pohl and Jovin, was that the conventional right-handed B-DNA structure (R) of poly[d(G-C)] transforms at high salt concentration into a novel, alternative left-handed (L) conformation. The historical course of this development and its aftermath, culminating in the first crystal structure of left-handed Z-DNA in 1979, is described in detail. The research conducted by Pohl and Jovin after 1979 is summarized, ending with an assessment of "unfinished business": condensed Z*-DNA; topoisomerase IIα (TOP2A) as an allosteric ZBP (Z-DNA-binding protein); B-Z transitions of phosphorothioate-modified DNAs; and parallel-stranded poly[d(G-A)], a double helix with high stability under physiological conditions and potentially also left-handed.

Ion distributions around left- and right-handed DNA and RNA duplexes: a comparative study

The ion atmosphere around nucleic acids is an integral part of their solvated structure. However, detailed aspects of the ionic distribution are difficult to probe experimentally, and comparative studies for different structures of the same sequence are almost non-existent. Here, we have used large-scale molecular dynamics simulations to perform a comparative study of the ion distribution around (5'-CGCGCGCGCGCG-3')2 dodecamers in solution in B-DNA, A-RNA, Z-DNA and Z-RNA forms. The CG sequence is very sensitive to ionic strength and it allows the comparison with the rare but important left-handed forms. The ions investigated include Na(+), K(+) and Mg(2 +), with various concentrations of their chloride salts. Our results quantitatively describe the characteristics of the ionic distributions for different structures at varying ionic strengths, tracing these differences to nucleic acid structure and ion type. Several binding pockets with rather long ion residence times are described, both for the monovalent ions and for the hexahydrated Mg[(H2O)6](2+) ion. The conformations of these binding pockets include direct binding through desolvated ion bridges in the GpC steps in B-DNA and A-RNA; direct binding to backbone oxygens; binding of Mg[(H2O)6](2+) to distant phosphates, resulting in acute bending of A-RNA; tight 'ion traps' in Z-RNA between C-O2 and the C-O2' atoms in GpC steps; and others.

Reaction path ensemble of the B-Z-DNA transition: a comprehensive atomistic study

Since its discovery in 1979, left-handed Z-DNA has evolved from an in vitro curiosity to a challenging DNA structure with crucial roles in gene expression, regulation and recombination. A fundamental question that has puzzled researchers for decades is how the transition from B-DNA, the prevalent right-handed form of DNA, to Z-DNA is accomplished. Due to the complexity of the B–Z-DNA transition, experimental and computational studies have resulted in several different, apparently contradictory models. Here, we use molecular dynamics simulations coupled with state-of-the-art enhanced sampling techniques operating through non-conventional reaction coordinates, to investigate the B–Z-DNA transition at the atomic level. Our results show a complex free energy landscape, where several phenomena such as over-stretching, unpeeling, base pair extrusion and base pair flipping are observed resulting in interconversions between different DNA conformations such as B-DNA, Z-DNA and S-DNA. In particular, different minimum free energy paths allow for the coexistence of different mechanisms (such as zipper and stretch–collapse mechanisms) that previously had been proposed as independent, disconnected models. We find that the B–Z-DNA transition—in absence of other molecular partners—can encompass more than one mechanism of comparable free energy, and is therefore better described in terms of a reaction path ensemble.

AN OVERVIEW OF THE STRETCHED INTERMEDIATE MODEL OF B–Z DNA TRANSITION

Recently, the stretched intermediate model was proposed for the B–Z deoxyribonucleic acid (DNA) transition based on simulation results carried out using the Stochastic Difference Equation (SDE) that showed unwinding and elongation of the oligomer during the transition. This model has proven to be successful in resolving the steric dilemma in short oligomers. However, extending the simulation method to larger DNA strands may prove to be computationally challenging. Such difficulty has led us to adopt a mean field approach using phenomenological interaction potentials to simulate the transition. Like the atomistic approach, the SDE simulations based on the mean field approach, also suggest the presence of a stretched intermediate during the transition.

Z-DNA crystallography

We review the effect of sequence on the structure of left-handed Z-DNA in single crystals. The various substituent groups that define a nucleotide base as guanine, cytosine,thymine, or adenine affect both the DNA conformation and the organization of solvent around the duplex. These are discussed in terms of their effect on the ability of sequences to adopt the unusual Z-DNA structure. In addition, the experimental and theoretical methods used to treat DNA hydration are discussed as they relate to the stability of Z-DNA . Finally, we argue that Z-DNA , as defined by the crystal conformation, is sufficient in itself to account for the physical properties of left-handed conformations observed in polymers and in genomic sequences

Solitary excitations in B-Z DNA transition: a theoretical and numerical study

The molecular mechanism of B-Z DNA transition remains elusive since the elucidation of the left-handed Z-DNA structure using atomic resolution crystallographic study. Numerous proposals for the molecular mechanism have been advanced, but none has provided a satisfactory explanation for the process. A nonlinear DNA model is proposed which enables one to derive various hypothesized molecular mechanisms, namely the Harvey model, Zang and Olson model, and the stretched intermediate model, by imposing certain constraints and conditions on the model. These constraints raise the need to reevaluate experimental investigations on B-Z DNA transition.

Intrinsic fluorescence of B and Z forms of poly d(G-m5C)·poly d(G-m5C), a synthetic double-stranded DNA: spectra and lifetimes by the maximum entropy method

A study has been made of the fluorescence of poly d(G-m5C).poly d(G-m5C), a synthetic double-stranded DNA, in buffered neutral aqueous solution at room temperature, excited by synchrotron radiation at 280 nm and 250 nm and by a frequency-doubled pulse dye laser at 290 nm. Exciting at 280 nm, the B form shows a uni-modal UV spectrum with lambdaf(max) approximately 340 nm. The Z form has in addition a visible emission lambdaf(max) at 450 nm. The spectral positions remain unchanged on exciting at 250 nm but the relative intensities change considerably. Decay profiles have been obtained at 360 nm and 450 nm for both the B and Z forms and have been analyzed by fitting to a pseudo-continuous distribution of 100 (and occasionally 200) exponentials, ranging from 10 ps to 20 ns, by optimizing the 'entropy' of the signal (the method of maximum entropy). We find the mean lifetimes for both wavelengths of emission and for both structural forms fall into three well-separated regions in the ranges indicated tau1 approximately 0.04-0.21 ns, tau2 approximately 0.9-1.26 ns, and tau3 approximately 5.1-6.5 ns. The UV emission, from its spectral position and half-width, correlates with monomeric emission from m5C (and from C for poly d(G-C)). However the lifetime tau1 is approximately 2 orders of magnitude longer than the monomers and points to an involvement of protonated guanosine (GH+, tauf approximately 200 ps) in the overall absorption/emission sequence. In the UV the tau3 emission is predominant, with fractional time-integrated emission approximately 86% for B DNA and approximately 64% for Z. We suggest it results from exciton (stacked) absorption followed by dissociative emission. For Z DNA the visible (450 nm) emission is dominated by a tau3 species (approximately 91%) with a lifetime of 6.5 ns and we suggest it represents a hetero-excimer emission consequent upon absorption by the strongly overlapped base-stacking, which differs from that in B DNA. The weak emission corresponding to tau2 is made more apparent by scanned gated detection of the emission from laser excitation (290 nm) of single-crystal d(m5C-G)3. A central role is attributed to the tight stacking of the bases in the Z form which correlates with enhanced hypochromism at 250 nm vs. 280 nm and with the reversal of the fluorescence intensity ratios UV-visible between these wavelengths.

Applying the stochastic difference equation to DNA conformational transitions: a study of B-Z and B-A DNA transitions

Despite the existence of numerous models to account for the B-Z DNA transition, experimenters have not yet arrived at a conclusive answer to the structural and dynamical features of the B-Z transition. By applying the stochastic difference equation to simulate the B-Z DNA transition, we have shown that the stretched intermediate model of the B-Z transition is more probable than other B-Z transition models such as the Harvey model. This is accomplished by comparing potential energy profiles of various B-Z DNA transition models and calculating relative probabilities based on the stochastic difference equation with respect to length (SDEL) formalism. The results garnered in this article allow for new approaches in determining the structural transition of B-DNA to Z-DNA experimentally. We have also simulated the B-A DNA transition using the stochastic difference equation. Unlike the B-Z DNA transition, the mechanism for the B-A DNA transition is well established. The variation in the pseudorotation angle during the transition is in good agreement with experimental results. Qualitative features of the simulated B-A transition also agree well with experimental data. The SDEL approach is thus a suitable numerical technique to compute long-time molecular dynamics trajectory for DNA molecules.

The stretched intermediate model of B-Z DNA transition

There have been numerous attempts to describe the mechanism of B-Z transition. Our simulations based on the stochastic difference equation with length algorithm show that a short DNA oligomer will tend to unwind and overstretch during the transition. The overstretching of DNA is then understood from the Zhou, Zhang, and Ou-Yang model. Unlike the Harvey model, the stretched intermediate model does not pose any steric dilemma; we are able to show that the chain sense reversal progresses spontaneously using the stretched intermediate model. A nonlinear DNA model is used to describe the origins and mechanism of base rotation in the stretched intermediate state of DNA. We also propose an experiment that can verify the existence of a stretched intermediate state during B-Z transition, thus opening up fresh grounds for experimentation in this field.

The halting arrival of left-handed Z-DNA

Forty-nine years ago Watson and Crick proposed a double-stranded (ds-) model for DNA. This double helix has become an icon of molecular biology. Twenty-six years later, Rich accidently discovered Z-DNA, an exotic left-handed nucleic acid. For many years thereafter, this left-handed DNA was thought to be an artifact. DNA is no longer looked upon as a static molecule but rather an extremely dynamic structure in which different conformations are in equilibrium with each other. Many researchers have spent the last two decades characterizing this novel left-handed DNA structure. Now many investigators are beginning to accept the possibility that this novel ds-DNA conformation may play a significant in vivo role within eukaryotic and prokaryotic cells. However, more research needs to be performed before it is absolutely accepted by all in the scientific community.

Structure and conformation of helical nucleic acids: analysis program (SCHNAaP)

We present a new versatile program, SCHNAaP, for the analysis of double-helical nucleic acid structures. The program uses mathematically rigorous and fully reversible procedures for calculating the structural parameters: the Cambridge University Engineering Department Helix computation Scheme (CEHS) is used to determine the local helical parameters and an analogous procedure is used to determine the global helical parameters. These parameters form a complete set that conforms to the "Cambridge Accord" on definitions and nomenclature of nucleic acid structure parameters. In addition to the two standard Watson-Crick base-pairs, the program handles mismatched base-pairs and chemically modified bases. An analysis of the sugar-phosphate backbone conformation is included. Standardized base-stacking diagrams of each dinucleotide step with reference to the mid-step triad are generated. Structures are classified as one of the four polymorphic families, A/B, Z, W or R, although W- and R-DNA (two types of hypothetical structure) have yet to be observed experimentally.

Structure and conformation of helical nucleic acids: rebuilding program (SCHNArP)

We present a program, SCHNArP, for rebuilding double-helical nucleic acid structures from a set of helical parameters. The parameter sets are based on mathematically reversible schemes that allow direct comparison of data from experimental X-ray crystal structures analyzed using the analysis program, SCHNAaP (see accompanying paper), and structures built using the rebuilding program, SCHNArP. The program uses either local CEHS helical parameters or global helical parameters. A number of standard parameter sets from the literature are included that allow comparison of oligomer and polymer structures generated using different models for sequence-dependent DNA bending. Exact atomic models are provided for the bases. Schematic models that trace the path of the backbone and use rectangular blocks for the bases can be generated.

Refinement of the solution structure of a branched DNA three-way junction

We have refined the structure of the DNA Three-Way Junction complex, TWJ-TC, described in the companion paper by quantitative analysis of two 2D NOESY spectra (mixing times 60 and 200 ms) obtained in D2O solution. NOESY crosspeak intensities extracted from these spectra were used in two kinds of refinement procedure: 1) distance-restrained energy minimization (EM) and molecular dynamics (MD) and 2) full relaxation matrix back calculation refinement. The global geometry of the refined model is very similar to that of a published, preliminary model (Leontis, 1993). Two of the helical arms of the junction are stacked. These are Helix 1, defined by basepairs S1-G1/S3-C12 through S1-C5/S3-G8 and Helix 2, which comprises basepairs S1-C6/S2-G5 through S1-G10/S2-G1. The third helical arm (Helix 3), comprised of basepairs S2-C6/S3-G5 through S2-C10/S3-G1 extends almost perpendicularly from the axis defined by Helices 1 and 2. The bases S1-C5 and S1-C6 of Strand 1 are continuously stacked across the junction region. The conformation of this strand is close to that of B-form DNA along its entire length, including the S1-C5 to S1-C6 dinucleotide step at the junction. The two unpaired bases S3-T6 and S3-C7 lie outside of the junction along the minor groove of Helix 1 and largely exposed to solvent. Analysis of the refined structure reveals that the glycosidic bond of S3-T6 exists in the syn conformation, allowing the methyl group of this residue to contact the hydrophobic surface of the minor groove of Helix 1, at S3-G11. The helical parameters of the three helical arms of the structure exhibit only weak deviations from typical values for right-handed B-form DNA. Unusual dihedral angles are only observed for the sugarphosphate backbone joining the "hinge" residues, S2-G5 and S2-C6, and S3-G5 through S3-G8. The glycosidic bond of S3-G8also lies within the syn range, allowing favorable Watson-Crick base-pairing interactions with Si -C5. The stability of this structure was checked in 39 ps molecular dynamic simulation at 330 K in water. The structure of TWJ-TC retained the geometrical features mentioned above at the end of the simulation period. The final R(1/6)-factor of the refined structure is 5%.