Stabilization of Z DNA in vivo by localized supercoiling

Contrasting effects of mismatch locations on Z-DNA formation under bending force

Z-DNA is a non-canonical, left-handed helical structure that plays crucial roles in various cellular processes. DNA mismatches, which involve the incorporation of incorrect Watson-Crick base pairs, are present in all living organisms and contribute to the mechanism of Z-DNA formation. However, the impact of mismatches on Z-DNA formation remains poorly understood. Moreover, the combined effect of DNA mismatches and bending, a common biological phenomenon observed in vivo, has not yet been explored due to technological limitations. Here, using single-molecule FRET, we show that a mismatch inside the Z-DNA region, i.e., the CG repeat region, hinders Z-DNA formation. In stark contrast, however, a mismatch in the B-Z junction facilitates Z-DNA formation. When the bending force is applied on double stranded DNA, a mismatch in the B-Z junction releases the bending stress more effectively than one in the CG repeat region. These findings provide mechanical insights into the role of DNA mismatches and bending forces in regulating Z-DNA formation, whether promoting or inhibiting it in biological environments.

Global identification and functional characterization of Z-DNA in rice

Z-DNA is a left-handed double helix form of DNA that is believed to be involved in various DNA transactions. However, comprehensive investigations aimed at global profiling of Z-DNA landscapes are still missing in both humans and plants. We here report the development of two techniques: anti-Z-DNA antibody-based immunoprecipitation followed by sequencing (ZIP-seq), and cleavage under targets and tagmentation (CUT&TAG) for characterizing Z-DNA in nipponbare rice (Oryza sativa L., Japonica). We found that Z-DNA-IP+ (Z-DNA recognized by the antibody) exhibits distinct genomic features as compared to Z-DNA-IP- (Z-DNA not recognized by the antibody). The concomitant presence of G-quadruplexes (G4s) and i-motifs (iMs) may promote Z-DNA formation. DNA modifications such as DNA-6mA/-4acC generally disfavours Z-DNA formation, while modifications like DNA-5mC (CHH) and 8-oxodG promote it, highlighting the distinct roles of DNA base modifications in modulating Z-DNA formation. Importantly, Z-DNA located at transcription start sites (TSSs) enhances gene expression, whereas Z-DNA in genic regions represses it, underscoring its dual roles in regulating the expression of genes involved in fundamental biological functions and responses to salt stress. Furthermore, Z-DNA may play a role in transcriptional initiation and termination rather than in transcriptional elongation. Finally, the presence of Z-DNA in promoters is correlated with the coevolution of overlapping genes, thereby regulating gene domestication. Consequently, our study represents as a pivotal point and a solid foundation for reliably launching genome-wide investigations of Z-DNA, thereby advancing the understanding of Z-DNA biology in both plants and non-plant systems.

Z-nucleic acid sensor ZBP1 in sterile inflammation

Z-DNA binding protein 1 (ZBP1), a cytosolic nucleic acid sensor for Z-form nucleic acids (Z-NA), can detect both exogenous and endogenous nucleic acids. Upon sensing of self Z-NA or exposure to diverse noxious stimuli, ZBP1 regulates inflammation by activating nuclear factor kappa B and interferon regulating factor 3 signaling pathways. In addition, ZBP1 promotes the assembly of ZBP1 PANoptosome, which initiates caspase 3-mediated apoptosis, mixed lineage kinase domain like pseudokinase (MLKL)-mediated necroptosis, and gasdermin D (GSDMD)-mediated pyroptosis (PANoptosis), leading to the release of various damage-associated molecular patterns. Thereby, ZBP1 is implicated in the development and progression of diverse sterile inflammatory diseases. This review outlines the expression, structure, and function of ZBP1, along with its dual roles in controlling inflammation and cell death to orchestrate innate immunity in sterile inflammation, especially autoimmune diseases, and cancers. ZBP1 has emerged as an attractive therapeutic target for various sterile inflammatory diseases.

Extracellular G-quadruplexes and Z-DNA protect biofilms from DNase I, and G-quadruplexes form a DNAzyme with peroxidase activity

Many bacteria form biofilms to protect themselves from predators or stressful environmental conditions. In the biofilm, bacteria are embedded in a protective extracellular matrix composed of polysaccharides, proteins and extracellular DNA (eDNA). eDNA most often is released from lysed bacteria or host mammalian cells, and it is the only matrix component most biofilms appear to have in common. However, little is known about the form DNA takes in the extracellular space, and how different non-canonical DNA structures such as Z-DNA or G-quadruplexes might contribute to its function in the biofilm. The aim of this study was to determine if non-canonical DNA structures form in eDNA-rich staphylococcal biofilms, and if these structures protect the biofilm from degradation by nucleases. We grew Staphylococcus epidermidis biofilms in laboratory media supplemented with hemin and NaCl to stabilize secondary DNA structures and visualized their location by immunolabelling and fluorescence microscopy. We furthermore visualized the macroscopic biofilm structure by optical coherence tomography. We developed assays to quantify degradation of Z-DNA and G-quadruplex DNA oligos by different nucleases, and subsequently investigated how these enzymes affected eDNA in the biofilms. Z-DNA and G-quadruplex DNA were abundant in the biofilm matrix, and were often present in a web-like structures. In vitro, the structures did not form in the absence of NaCl or mechanical shaking during biofilm growth, or in bacterial strains deficient in eDNA or exopolysaccharide production. We thus infer that eDNA and polysaccharides interact, leading to non-canonical DNA structures under mechanical stress when stabilized by salt. We also confirmed that G-quadruplex DNA and Z-DNA was present in biofilms from infected implants in a murine implant-associated osteomyelitis model. Mammalian DNase I lacked activity against Z-DNA and G-quadruplex DNA, while Micrococcal nuclease could degrade G-quadruplex DNA and S1 Aspergillus nuclease could degrade Z-DNA. Micrococcal nuclease, which originates from Staphylococcus aureus, may thus be key for dispersal of biofilm in staphylococci. In addition to its structural role, we show for the first time that the eDNA in biofilms forms a DNAzyme with peroxidase-like activity in the presence of hemin. While peroxidases are part of host defenses against pathogens, we now show that biofilms can possess intrinsic peroxidase activity in the extracellular matrix.

Dynamic alternative DNA structures in biology and disease

Repetitive elements in the human genome, once considered 'junk DNA', are now known to adopt more than a dozen alternative (that is, non-B) DNA structures, such as self-annealed hairpins, left-handed Z-DNA, three-stranded triplexes (H-DNA) or four-stranded guanine quadruplex structures (G4 DNA). These dynamic conformations can act as functional genomic elements involved in DNA replication and transcription, chromatin organization and genome stability. In addition, recent studies have revealed a role for these alternative structures in triggering error-generating DNA repair processes, thereby actively enabling genome plasticity. As a driving force for genetic variation, non-B DNA structures thus contribute to both disease aetiology and evolution.

Machine learning shows torsion angle preferences in left-handed and right-handed quadruplex DNAs

Left-handed G quadruplexes (LHG4) have been recently discovered as a new class of G quadruplexes. The biological functions of LHG4s are still unknown, but they share a striking resemblance to Z-DNA in their helicity and jagged phosphate backbone. To further understand structural features of the LHG4s that define their left handedness, we have employed human-interpretable machine-learning methods to classify right- and left-handed G4s purely based on torsional angle analysis. Our results reveal the importance of the α, β, δ, and χ angles in left-handed structuring across both Z-DNAs and LHG4s. Our analysis may serve as the first step to understanding the conditions of formation for LHG4s and their potential biological relevance.

Z-nucleic acids: Uncovering the functions from past to present

Since Z-nucleic acid was identified in the 1970s, much is still unknown about its biological functions and nature in vivo. Recent studies on adenosine deaminase acting on RNA 1 (ADAR1) and Z-DNA-binding protein 1 (ZBP1) have highlighted its function in immune responses. Specifically, Z-RNAs, either endogenous or induced by viral infection, are sensed by ZBP1 and activate necroptosis. Z-RNAs act as the stimuli that induce innate immune responses through various pathways, including melanoma differentiation-associated protein 5 (MAD5)-mitochondrial antiviral-signaling protein (MAVS)-mediated type I IFN activation and proteinase kinase R (PKR)-dependent integrated stress response, and their immunostimulatory potential is curtailed by RNA editing conducted by ADAR1. Aberrant immune responses induced by Z-RNAs are associated with human diseases. They also induce pathogenesis in mice. Unlike Z-RNAs, the biological functions of Z-DNAs were barely studied, especially in mammals. Moreover, the origin or sequence preference of Z-nucleic acids requires further investigation. Such knowledge will expand our understanding of Z-nucleic acids, including from which genomic loci and under which circumstances they form, and the mechanisms by which they participate in the physiological activities. In this review, we provide insights in Z-nucleic acid research and highlight the unsolved puzzles.

DNA Bending Force Facilitates Z-DNA Formation under Physiological Salt Conditions

Z-DNA, a noncanonical helical structure of double-stranded DNA (dsDNA), plays pivotal roles in various biological processes, including transcription regulation. Mechanical stresses on dsDNA, such as twisting and stretching, help to form Z-DNA. However, the effect of DNA bending, one of the most common dsDNA deformations, on Z-DNA formation is utterly unknown. Here, we show that DNA bending induces the formation of Z-DNA, that is, more Z-DNA is formed as the bending force becomes stronger. We regulated the bending force on dsDNA by using D-shaped DNA nanostructures. The B-Z transition was observed by single-molecule fluorescence resonance energy transfer. We found that as the bending force became stronger, Z-DNA was formed at lower Mg²⁺ concentrations. When dsDNA contained cytosine methylations, the B-Z transition occurred at 78 mM Mg²⁺ (midpoint) in the absence of the bending force. However, the B-Z transition occurred at a 28-fold lower Mg²⁺ concentration (2.8 mM) in the presence of the bending force. Monte Carlo simulation suggested that the B-Z transition stabilizes the bent form via the formation of the B-Z junction with base extrusion, which effectively releases the bending stress on DNA. Our results clearly show that the bending force facilitates the B-Z transition under physiological salt conditions.

Nonalternating purine pyrimidine sequences can form stable left-handed DNA duplex by strong topological constraint

In vivo, left-handed DNA duplex (usually refers to Z-DNA) is mainly formed in the region of DNA with alternating purine pyrimidine (APP) sequence and plays significant biological roles. It is well known that d(CG)_n sequence can form Z-DNA most easily under negative supercoil conditions, but its essence has not been well clarified. The study on sequence dependence of Z-DNA stability is very difficult without modification or inducers. Here, by the strong topological constraint caused by hybridization of two complementary short circular ssDNAs, left-handed duplex part was generated for various sequences, and their characteristics were investigated by using gel-shift after binding to specific proteins, CD and T_m analysis, and restriction enzyme cleavage. Under the strong topological constraint, non-APP sequences can also form left-handed DNA duplex as stable as that of APP sequences. As compared with non-APP sequences, the thermal stability difference for APP sequences between Z-form and B-form is smaller, which may be the reason that Z-DNA forms preferentially for APP ones. This result can help us to understand why nature selected APP sequences to regulate gene expression by transient Z-DNA formation, as well as why polymer with chirality can usually form both duplexes with left- or right-handed helix.

Interaction of the Escherichia coli HU Protein with Various Topological Forms of DNA

E. coli histone-like protein HU has been shown to interact with different topological forms of DNA. Using radiolabeled HU, we examine the effects of DNA supercoiling on HU-DNA interactions. We show that HU binds preferentially to negatively supercoiled DNA and that the affinity of HU for DNA increases with increases in the negative superhelical density of DNA. Binding of HU to DNA is most sensitively influenced by DNA supercoiling within a narrow but physiologically relevant range of superhelicity (σ = -0.06-0). Under stoichiometric binding conditions, the affinity of HU for negatively supercoiled DNA (σ = -0.06) is more than 10 times higher than that for relaxed DNA at physiologically relevant HU/DNA mass ratios (e.g., 1:10). This binding preference, however, becomes negligible at HU/DNA mass ratios higher than 1:2. At saturation, HU binds both negatively supercoiled and relaxed DNA with similar stoichiometries, i.e., 5-6 base pairs per HU dimer. In our chemical crosslinking studies, we demonstrate that HU molecules bound to negatively supercoiled DNA are more readily crosslinked than those bound to linear DNA. At in vivo HU/DNA ratios, HU appears to exist predominantly in a tetrameric form on negatively supercoiled DNA and in a dimeric form on linear DNA. Using a DNA ligase-mediated nick closure assay, we show that approximately 20 HU dimers are required to constrain one negative supercoil on relaxed DNA. Although fewer HU dimers may be needed to constrain one negative supercoil on negatively supercoiled DNA, our results and estimates of the cellular level of HU argue against a major role for HU in constraining supercoils in vivo. We discuss our data within the context of the dynamic distribution of the HU protein in cells, where temporal and local changes of DNA supercoiling are known to take place.

Decoding the Structural Dynamics and Conformational Alternations of DNA Secondary Structures by Single-Molecule FRET Microspectroscopy

In addition to the canonical double helix form, DNA is known to be extrapolated into several other secondary structural patterns involving themselves in inter- and intramolecular type hydrogen bonding. The secondary structures of nucleic acids go through several stages of multiple, complex, and interconvertible heterogeneous conformations. The journey of DNA through these conformers has significant importance and has been monitored thoroughly to establish qualitative and quantitative information about the transition between the unfolded, folded, misfolded, and partially folded states. During this structural interconversion, there always exist specific populations of intermediates, which are short-lived or sometimes even do not accumulate within a heterogeneous population and are challenging to characterize using conventional ensemble techniques. The single-molecule FRET(sm-FRET) microspectroscopic method has the advantages to overcome these limitations and monitors biological phenomena transpiring at a measurable high rate and balanced stochastically over time. Thus, tracing the time trajectory of a particular molecule enables direct measurement of the rate constant of each transition step, including the intermediates that are hidden in the ensemble level due to their low concentrations. This review is focused on the advantages of the employment of single-molecule Forster's resonance energy transfer (sm-FRET), which is worthwhile to access the dynamic architecture and structural transition of various secondary structures that DNA adopts, without letting the donor of one molecule to cross-talk with the acceptor of any other. We have emphasized the studies performed to explore the states of folding and unfolding of several nucleic acid secondary structures, for example, the DNA hairpin, Holliday junction, G-quadruplex, and i-motif.

Dynamic regulation of Z-DNA in the mouse prefrontal cortex by the RNA-editing enzyme Adar1 is required for fear extinction

DNA forms conformational states beyond the right-handed double-helix; however, the functional relevance of these non-canonical structures in the brain remains unknown. We show that, in the prefrontal cortex of mice, the formation of one such structure, Z-DNA, is involved in the regulation of extinction memory. Z-DNA is formed during fear learning, and reduced during extinction learning, which is mediated, in part, by a direct interaction between Z-DNA and the RNA editing enzyme Adar1. Adar1 binds to Z-DNA during fear extinction learning which leads to a reduction in Z-DNA at sites where Adar1 is recruited. Knockdown of Adar1 leads to an inability to modify a previously acquired fear memory and blocks activity-dependent changes in DNA structure and RNA state; effects that are fully rescued by the introduction of full-length Adar1. These findings suggest a novel mechanism of learning-induced gene regulation dependent on both proteins which recognize DNA structure, and the state.

Chromosomal Organization and Regulation of Genetic Function in Escherichia coli Integrates the DNA Analog and Digital Information

In this article, we summarize our current understanding of the bacterial genetic regulation brought about by decades of studies using the Escherichia coli model. It became increasingly evident that the cellular genetic regulation system is organizationally closed, and a major challenge is to describe its circular operation in quantitative terms. We argue that integration of the DNA analog information (i.e., the probability distribution of the thermodynamic stability of base steps) and digital information (i.e., the probability distribution of unique triplets) in the genome provides a key to understanding the organizational logic of genetic control. During bacterial growth and adaptation, this integration is mediated by changes of DNA supercoiling contingent on environmentally induced shifts in intracellular ionic strength and energy charge. More specifically, coupling of dynamic alterations of the local intrinsic helical repeat in the structurally heterogeneous DNA polymer with structural-compositional changes of RNA polymerase holoenzyme emerges as a fundamental organizational principle of the genetic regulation system. We present a model of genetic regulation integrating the genomic pattern of DNA thermodynamic stability with the gene order and function along the chromosomal OriC-Ter axis, which acts as a principal coordinate system organizing the regulatory interactions in the genome.

Unveiling the pathway to Z-DNA in the protein-induced B-Z transition

Left-handed Z-DNA is an extraordinary conformation of DNA, which can form by special sequences under specific biological, chemical or physical conditions. Human ADAR1, prototypic Z-DNA binding protein (ZBP), binds to Z-DNA with high affinity. Utilizing single-molecule FRET assays for Z-DNA forming sequences embedded in a long inactive DNA, we measure thermodynamic populations of ADAR1-bound DNA conformations in both GC and TG repeat sequences. Based on a statistical physics model, we determined quantitatively the affinities of ADAR1 to both Z-form and B-form of these sequences. We also reported what pathways it takes to induce the B–Z transition in those sequences. Due to the high junction energy, an intermediate B* state has to accumulate prior to the B–Z transition. Our study showing the stable B* state supports the active picture for the protein-induced B–Z transition that occurs under a physiological setting.

Computational Approaches to Predict the Non-canonical DNAs

Background:

Although most nucleotides in the genome form canonical double-stranded B-DNA, many repeated sequences transiently present as non-canonical conformations (non-B DNA) such as triplexes, quadruplexes, Z-DNA, cruciforms, and slipped/hairpins. Those noncanonical DNAs (ncDNAs) are not only associated with many genetic events such as replication, transcription, and recombination, but are also related to the genetic instability that results in the predisposition to disease. Due to the crucial roles of ncDNAs in cellular and genetic functions, various computational methods have been implemented to predict sequence motifs that generate ncDNA.

Objective:

Here, we review strategies for the identification of ncDNA motifs across the whole genome, which is necessary for further understanding and investigation of the structure and function of ncDNAs.

Conclusion:

There is a great demand for computational prediction of non-canonical DNAs that play key functional roles in gene expression and genome biology. In this study, we review the currently available computational methods for predicting the non-canonical DNAs in the genome. Current studies not only provide an insight into the computational methods for predicting the secondary structures of DNA but also increase our understanding of the roles of non-canonical DNA in the genome.

Topologically Constrained Formation of Stable Z-DNA from Normal Sequence under Physiological Conditions

Z-DNA, a left-handed duplex, has been shown to form in vivo and regulate expression of the corresponding gene. However, its biological roles have not been satisfactorily understood, mainly because Z-DNA is easily converted to the thermodynamically favorable B-DNA. Here we present a new idea to form stable Z-DNA under normal physiological conditions and achieve detailed analysis on its fundamental features. Simply by mixing two complementary minicircles of single-stranded DNA with no chemical modification, the hybridization spontaneously induces topological constraint which twines one-half of the double-stranded DNA into stable Z-DNA. The formation of Z-conformation with high stability has been proved by using circular dichroism spectroscopy, Z-DNA-specific antibody binding assay, nuclease digestion, etc. Even at a concentration of MgCl₂ as low as 0.5 mM, Z-DNA was successfully obtained, avoiding the use of high salt conditions, limited sequences, ancillary additives, or chemical modifications, criteria which have hampered Z-DNA research. The resultant Z-DNA has the potential to be used as a canonical standard sample in Z-DNA research. By using this approach, further developments of Z-DNA science and its applications become highly promising.

Free Energy Landscape and Dynamics of Supercoiled DNA by High-Speed Atomic Force Microscopy

DNA supercoiling fundamentally constrains and regulates the storage and use of genetic information. While the equilibrium properties of supercoiled DNA are relatively well understood, the dynamics of supercoils are much harder to probe. Here we use atomic force microscopy (AFM) imaging to demonstrate that positively supercoiled DNA plasmids, in contrast to their negatively supercoiled counterparts, preserve their plectonemic geometry upon adsorption under conditions that allow for dynamics and equilibration on the surface. Our results are in quantitative agreement with a physical polymer model for supercoiled plasmids that takes into account the known mechanical properties and torque-induced melting of DNA. We directly probe supercoil dynamics using high-speed AFM imaging with subsecond time and ∼nanometer spatial resolution. From our recordings we quantify self-diffusion, branch point flexibility, and slithering dynamics and demonstrate that reconfiguration of molecular extensions is predominantly governed by the bending flexibility of plectoneme arms. We expect that our methodology can be an asset to probe protein-DNA interactions and topochemical reactions on physiological relevant DNA length and supercoiling scales by high-resolution AFM imaging.

DNA Conformation Regulates Gene Expression: The MYC Promoter and Beyond

Emerging evidence suggests that DNA topology plays an instructive role in cell fate control through regulation of gene expression. Transcription produces torsional stress, and the resultant supercoiling of the DNA molecule generates an array of secondary structures. In turn, local DNA architecture is harnessed by the cell, acting within sensory feedback mechanisms to mediate transcriptional output. MYC is a potent oncogene, which is upregulated in the majority of cancers; thus numerous studies have focused on detailed understanding of its regulation. Dissection of regulatory regions within the MYC promoter provided the first hint that intimate feedback between DNA topology and associated DNA remodeling proteins is critical for moderating transcription. As evidence of such regulation is also found in the context of many other genes, here we expand on the prototypical example of the MYC promoter, and also explore DNA architecture in a genome-wide context as a global mechanism of transcriptional control.

Osmium tetroxide as a probe of RNA structure

Structured RNAs have a central role in cellular function. The capability of structured RNAs to adopt fixed architectural structures or undergo dynamic conformational changes contributes to their diverse role in the regulation of gene expression. Although numerous biophysical and biochemical tools have been developed to study structured RNAs, there is a continuing need for the development of new methods for the investigation of RNA structures, especially methods that allow RNA structure to be studied in solution close to its native cellular conditions. Here we use osmium tetroxide (OsO₄) as a chemical probe of RNA structure. In this method, we have used fluorescence-based sequencing technologies to detect OsO₄ modified RNA. We characterized the requirements for OsO₄ modification of RNA by investigating three known structured RNAs: the M-box, glycine riboswitch RNAs, and tRNA^asp Our results show that OsO₄ predominantly modifies RNA at uracils that are conformationally exposed on the surface of the RNA. We also show that changes in OsO₄ reactivity at flexible positions in the RNA correlate with ligand-driven conformational changes in the RNA structure. Osmium tetroxide modification of RNA will provide insights into the structural features of RNAs that are relevant to their underlying biological functions.

Opening the Strands of Replication Origins-Still an Open Question

The local separation of duplex DNA strands (strand opening) is necessary for initiating basic transactions on DNA such as transcription, replication, and homologous recombination. Strand opening is commonly a stage at which these processes are regulated. Many different mechanisms are used to open the DNA duplex, the details of which are of great current interest. In this review, we focus on a few well-studied cases of DNA replication origin opening in bacteria. In particular, we discuss the opening of origins that support the theta (θ) mode of replication, which is used by all chromosomal origins and many extra-chromosomal elements such as plasmids and phages. Although the details of opening can vary among different origins, a common theme is binding of the initiator to multiple sites at the origin, causing stress that opens an adjacent and intrinsically unstable A+T rich region. The initiator stabilizes the opening by capturing one of the open strands. How the initiator binding energy is harnessed for strand opening remains to be understood.

Biosynthesis and Regulation of the Branched-Chain Amino Acids†

This review focuses on more recent studies concerning the systems biology of branched-chain amino acid biosynthesis, that is, the pathway-specific and global metabolic and genetic regulatory networks that enable the cell to adjust branched-chain amino acid synthesis rates to changing nutritional and environmental conditions. It begins with an overview of the enzymatic steps and metabolic regulatory mechanisms of the pathways and descriptions of the genetic regulatory mechanisms of the individual operons of the isoleucine-leucine-valine (ilv) regulon. This is followed by more-detailed discussions of recent evidence that global control mechanisms that coordinate the expression of the operons of this regulon with one another and the growth conditions of the cell are mediated by changes in DNA supercoiling that occur in response to changes in cellular energy charge levels that, in turn, are modulated by nutrient and environmental signals. Since the parallel pathways for isoleucine and valine biosynthesis are catalyzed by a single set of enzymes, and because the AHAS-catalyzed reaction is the first step specific for valine biosynthesis but the second step of isoleucine biosynthesis, valine inhibition of a single enzyme for this enzymatic step might compromise the cell for isoleucine or result in the accumulation of toxic intermediates. The operon-specific regulatory mechanisms of the operons of the ilv regulon are discussed in the review followed by a consideration and brief review of global regulatory proteins such as integration host factor (IHF), Lrp, and CAP (CRP) that affect the expression of these operons.

A SERS/fluorescence dual-mode nanosensor based on the human telomeric G-quadruplex DNA: Application to mercury (II) detection

DNA-metal nanoparticle conjugates have been increasingly exploited for sensing purposes over the past decades. However, most of the existing strategies are operated with canonical DNA structures, such as single-stranded forms, stem-loop structures, and double helix structures. There is intense interest in the development of nano-system based on high order DNA secondary structures. Herein, we propose a SERS/fluorescence dual-mode nanosensor, where the signal transduction mechanism is based on the conformational switching of the human telomeric G-quadruplex DNA. The nanosensor exhibits excellent SERS/fluorescence responses to the complementary strands of G-quadruplexes. Based on T-Hg(2+)-T coordination chemistry, this sensor is effectively applied to determination of Hg(2+) in buffer solution and real samples. It achieves a limit of detection (LOD) as low as 1ppt, which is ~100 times more sensitive than conventional optical sensors. We anticipate that the proposed G-quadruplex-based nanosensor could be applied to the analysis of other metal ions and small molecules in environmental samples and biological systems.

Conformation-dependent DNA attraction

Understanding how DNA molecules interact with other biomolecules is related to how they utilize their functions and is therefore critical for understanding their structure-function relationships. For a long time, the existence of Z-form DNA (a left-handed double helical version of DNA, instead of the common right-handed B-form) has puzzled the scientists, and the definitive biological significance of Z-DNA has not yet been clarified. In this study, the effects of DNA conformation in DNA-DNA interactions are explored by molecular dynamics simulations. Using umbrella sampling, we find that for both B- and Z-form DNA, surrounding Mg(2+) ions always exert themselves to screen the Coulomb repulsion between DNA phosphates, resulting in very weak attractive force. On the contrary, a tight and stable bound state is discovered for Z-DNA in the presence of Mg(2+) or Na(+), benefiting from their hydrophobic nature. Based on the contact surface and a dewetting process analysis, a two-stage binding process of Z-DNA is outlined: two Z-DNA first attract each other through charge screening and Mg(2+) bridges to phosphate groups in the same way as that of B-DNA, after which hydrophobic contacts of the deoxyribose groups are formed via a dewetting effect, resulting in stable attraction between two Z-DNA molecules. The highlighted hydrophobic nature of Z-DNA interaction from the current study may help to understand the biological functions of Z-DNA in gene transcription.

Distinct Z-DNA binding mode of a PKR-like protein kinase containing a Z-DNA binding domain (PKZ)

Double-stranded ribonucleic acid-activated protein kinase (PKR) downregulates translation as a defense mechanism against viral infection. In fish species, PKZ, a PKR-like protein kinase containing left-handed deoxyribonucleic acid (Z-DNA) binding domains, performs a similar role in the antiviral response. To understand the role of PKZ in Z-DNA recognition and innate immune response, we performed structural and functional studies of the Z-DNA binding domain (Zα) of PKZ from Carassius auratus (caZα_PKZ). The 1.7-Å resolution crystal structure of caZα_PKZ:Z-DNA revealed that caZα_PKZ shares the overall fold with other Zα, but has discrete structural features that differentiate its DNA binding mode from others. Functional analyses of caZα_PKZ and its mutants revealed that caZα_PKZ mediates the fastest B-to-Z transition of DNA among Zα, and the minimal interaction for Z-DNA recognition is mediated by three backbone phosphates and six residues of caZα_PKZ. Structure-based mutagenesis and B-to-Z transition assays confirmed that Lys56 located in the β-wing contributes to its fast B-to-Z transition kinetics. Investigation of the DNA binding kinetics of caZα_PKZ further revealed that the B-to-Z transition rate is positively correlated with the association rate constant. Taking these results together, we conclude that the positive charge in the β-wing largely affects fast B-to-Z transition activity by enhancing the DNA binding rate.

Energetics of Z-DNA binding protein-mediated helicity reversals in DNA, RNA, and DNA-RNA duplexes

Z-DNA binding proteins (ZBPs) specifically recognize and stabilize left-handed double helices, including Z-DNA and Z-RNA. However, the energetics of Z-form stabilization by ZBPs have never been characterized due to the technical limitations of bulk studies, resulting in an unclear understanding of the ZBP operational mechanism at the molecular level. Here, we use single-molecule fluorescence resonance energy transfer (FRET) to determine the energetics of Z-form stabilization by ZBP for DNA, RNA, and DNA-RNA duplexes, revealing that the formation of B-Z or A-Z junctions dominates the thermodynamics and kinetics of Z-form stabilization. Furthermore, in contrast to general assumptions, the Z-form is most efficiently and most rapidly formed in the DNA-RNA hybrid duplex due to the greatly reduced junction energy in the DNA-RNA hybrid.

Competitive superhelical transitions involving cruciform extrusion

A DNA molecule under negative superhelical stress becomes susceptible to transitions to alternate structures. The accessible alternate conformations depend on base sequence and compete for occupancy. We have developed a method to calculate equilibrium distributions among the states available to such systems, as well as their average thermodynamic properties. Here we extend this approach to include superhelical cruciform extrusion at both perfect and imperfect inverted repeat (IR) sequences. We find that short IRs do not extrude cruciforms, even in the absence of competition. But as the length of an IR increases, its extrusion can come to dominate both strand separation and B-Z transitions. Although many IRs are present in human genomic DNA, we find that extrusion-susceptible ones occur infrequently. Moreover, their avoidance of transcription start sites in eukaryotes suggests that cruciform formation is rarely involved in mechanisms of gene regulation. We examine a set of clinically important chromosomal translocation breakpoints that occur at long IRs, whose rearrangement has been proposed to be driven by cruciform extrusion. Our results show that the susceptibilities of these IRs to cruciform formation correspond closely with their observed translocation frequencies.

DNA secondary structures: stability and function of G-quadruplex structures

In addition to the canonical double helix, DNA can fold into various other inter- and intramolecular secondary structures. Although many such structures were long thought to be in vitro artefacts, bioinformatics demonstrates that DNA sequences capable of forming these structures are conserved throughout evolution, suggesting the existence of non-B-form DNA in vivo. In addition, genes whose products promote formation or resolution of these structures are found in diverse organisms, and a growing body of work suggests that the resolution of DNA secondary structures is critical for genome integrity. This Review focuses on emerging evidence relating to the characteristics of G-quadruplex structures and the possible influence of such structures on genomic stability and cellular processes, such as transcription.

Theoretical analysis of competing conformational transitions in superhelical DNA

We develop a statistical mechanical model to analyze the competitive behavior of transitions to multiple alternate conformations in a negatively supercoiled DNA molecule of kilobase length and specified base sequence. Since DNA superhelicity topologically couples together the transition behaviors of all base pairs, a unified model is required to analyze all the transitions to which the DNA sequence is susceptible. Here we present a first model of this type. Our numerical approach generalizes the strategy of previously developed algorithms, which studied superhelical transitions to a single alternate conformation. We apply our multi-state model to study the competition between strand separation and B-Z transitions in superhelical DNA. We show this competition to be highly sensitive to temperature and to the imposed level of supercoiling. Comparison of our results with experimental data shows that, when the energetics appropriate to the experimental conditions are used, the competition between these two transitions is accurately captured by our algorithm. We analyze the superhelical competition between B-Z transitions and denaturation around the c-myc oncogene, where both transitions are known to occur when this gene is transcribing. We apply our model to explore the correlation between stress-induced transitions and transcriptional activity in various organisms. In higher eukaryotes we find a strong enhancement of Z-forming regions immediately 5′ to their transcription start sites (TSS), and a depletion of strand separating sites in a broad region around the TSS. The opposite patterns occur around transcript end locations. We also show that susceptibility to each type of transition is different in eukaryotes and prokaryotes. By analyzing a set of untranscribed pseudogenes we show that the Z-susceptibility just downstream of the TSS is not preserved, suggesting it may be under selection pressure.

The stresses imposed on DNA within organisms can drive the molecule from its standard B-form double-helical structure into other conformations at susceptible sites within the sequence. We present a theoretical method to calculate this transition behavior due to stresses induced by supercoiling. We also develop a numerical algorithm that calculates the transformation probability of each base pair in a user-specified DNA sequence under stress. We apply this method to analyze the competition between transitions to strand separated and left-handed Z-form structures. We find that these two conformations are both competitive under physiological environmental conditions, and that this competition is especially sensitive to temperature. By comparing its results to experimental data we also show that the algorithm properly describes the competition between melting and Z-DNA formation. Analysis of large gene sets from various organisms shows a correlation between sites of stress-induced transitions and locations that are involved in regulating gene expression.

Cooperative kinking at distant sites in mechanically stressed DNA

In cells, DNA is routinely subjected to significant levels of bending and twisting. In some cases, such as under physiological levels of supercoiling, DNA can be so highly strained, that it transitions into non-canonical structural conformations that are capable of relieving mechanical stress within the template. DNA minicircles offer a robust model system to study stress-induced DNA structures. Using DNA minicircles on the order of 100 bp in size, we have been able to control the bending and torsional stresses within a looped DNA construct. Through a combination of cryo-EM image reconstructions, Bal31 sensitivity assays and Brownian dynamics simulations, we have been able to analyze the effects of biologically relevant underwinding-induced kinks in DNA on the overall shape of DNA minicircles. Our results indicate that strongly underwound DNA minicircles, which mimic the physical behavior of small regulatory DNA loops, minimize their free energy by undergoing sequential, cooperative kinking at two sites that are located about 180° apart along the periphery of the minicircle. This novel form of structural cooperativity in DNA demonstrates that bending strain can localize hyperflexible kinks within the DNA template, which in turn reduces the energetic cost to tightly loop DNA.

The torsional state of DNA within the chromosome

Virtually all processes of the genome biology affect or are affected by the torsional state of DNA. Torsional energy associated with an altered twist facilitates or hinders the melting of the double helix, its molecular interactions, and its spatial folding in the form of supercoils. Yet, understanding how the torsional state of DNA is modulated remains a challenging task due to the multiplicity of cellular factors involved in the generation, transmission, and dissipation of DNA twisting forces. Here, an overview of the implication of DNA topoisomerases, DNA revolving motors, and other DNA interactions that determine local levels of torsional stress in bacterial and eukaryotic chromosomes is provided. Particular emphasis is made on the experimental approaches being developed to assess the torsional state of intracellular DNA and its organization into topological domains.

Theoretical analysis of the stress induced B-Z transition in superhelical DNA

We present a method to calculate the propensities of regions within a DNA molecule to transition from B-form to Z-form under negative superhelical stresses. We use statistical mechanics to analyze the competition that occurs among all susceptible Z-forming regions at thermodynamic equilibrium in a superhelically stressed DNA of specified sequence. This method, which we call SIBZ, is similar to the SIDD algorithm that was previously developed to analyze superhelical duplex destabilization. A state of the system is determined by assigning to each base pair either the B- or the Z-conformation, accounting for the dinucleotide repeat unit of Z-DNA. The free energy of a state is comprised of the nucleation energy, the sequence-dependent B-Z transition energy, and the energy associated with the residual superhelicity remaining after the change of twist due to transition. Using this information, SIBZ calculates the equilibrium B-Z transition probability of each base pair in the sequence. This can be done at any physiologically reasonable level of negative superhelicity. We use SIBZ to analyze a variety of representative genomic DNA sequences. We show that the dominant Z-DNA forming regions in a sequence can compete in highly complex ways as the superhelicity level changes. Despite having no tunable parameters, the predictions of SIBZ agree precisely with experimental results, both for the onset of transition in plasmids containing introduced Z-forming sequences and for the locations of Z-forming regions in genomic sequences. We calculate the transition profiles of 5 kb regions taken from each of 12,841 mouse genes and centered on the transcription start site (TSS). We find a substantial increase in the frequency of Z-forming regions immediately upstream from the TSS. The approach developed here has the potential to illuminate the occurrence of Z-form regions in vivo, and the possible roles this transition may play in biological processes.

Adenosine-to-inosine RNA editing

Ribonucleic acid (RNA) editing is a mechanism that generates RNA and protein diversity, which is not directly encoded in the genome. The most common type of RNA editing in vertebrates is the conversion of adenosine to inosine in double-stranded RNA which occurs in the higher eukaryotes. This editing is carried out by the family of adenosine deaminase acting on RNA (ADAR) proteins. The most-studied substrates of ADAR proteins undergo editing which is very consistent, highly conserved, and functionally important. However, editing causes changes in protein-coding regions only at a small proportion of all editing sites. The vast majority of editing sites are in noncoding sequences. This includes microRNAs, as well as the introns and 3' untranslated regions of messenger RNAs, which play important roles in the RNA-mediated regulation of gene expression.

Minute negative superhelicity is sufficient to induce the B-Z transition in the presence of low tension

Left-handed Z-DNA has fascinated biological scientists for decades by its extraordinary structure and potential involvement in biological phenomena. Despite its instability relative to B-DNA, Z-DNA is stabilized in vivo by negative supercoiling. A detailed understanding of Z-DNA formation is, however, still lacking. In this study, we have examined the B-Z transition in a short guanine/cytosine (GC) repeat in the presence of controlled tension and superhelicity via a hybrid technique of single-molecule FRET and magnetic tweezers. The hybrid scheme enabled us to identify the states of the specific GC region under mechanical control and trace conformational changes synchronously at local and global scales. Intriguingly, minute negative superhelicity can facilitate the B-Z transition at low tension, indicating that tension, as well as torsion, plays a pivotal role in the transition. Dynamic interconversions between the states at elevated temperatures yielded thermodynamic and kinetic constants of the transition. Our single-molecule studies shed light on the understanding of Z-DNA formation by highlighting the highly cooperative and dynamic nature of the B-Z transition.

Methods to determine DNA structural alterations and genetic instability

Chromosomal DNA is a dynamic structure that can adopt a variety of non-canonical (i.e., non-B) conformations. In this regard, at least 10 different forms of non-B DNA conformations have been identified; many of them have been found to be mutagenic, and associated with human disease development. Despite the importance of non-B DNA structures in genetic instability and DNA metabolic processes, mechanisms by which instability occurs remain largely undefined. The purpose of this review is to summarize current methodologies that are used to address questions in the field of non-B DNA structure-induced genetic instability. Advantages and disadvantages of each method will be discussed. A focused effort to further elucidate the mechanisms of non-B DNA-induced genetic instability will lead to a better understanding of how these structure-forming sequences contribute to the development of human disease.

Comparative analyses of distributions and functions of Z-DNA in Arabidopsis and rice

Left-handed Z-DNA is an energetically unfavorable DNA structure that could form mostly under certain physiological conditions and was known to be involved in a number of cellular activities such as transcription regulation. We have compared the distributions and functions of Z-DNA in the genomes of Arabidopsis and rice, and observed that Z-DNA occurs in rice at least 9 times more often than in Arabidopsis; similar observations hold for other monocots and dicots. In addition, Z-DNA is significantly enriched in the coding regions of Arabidopsis, and in the high-GC-content regions of rice. Based on our analyses, we speculate that Z-DNA may play a role in regulating the expression of transcription factors, inhibitors, translation repressors, succinate dehydrogenases and glutathione-disulfide reductases in Arabidopsis, and it may affect the expression of vesicle and nucleosome genes and genes involved in alcohol transporter activity, stem cell maintenance, meristem development and reproductive structure development in rice.

Effect of sequence-directed nucleosome disruption on cell-type-specific repression by alpha2/Mcm1 in the yeast genome

In Saccharomyces cerevisiae, a-cell-specific genes are repressed in MATalpha cells by alpha2/Mcm1, acting in concert with the Ssn6-Tup1 corepressors and the Isw2 chromatin remodeling complex, and nucleosome positioning has been proposed as one mechanism of repression. However, prior studies showed that nucleosome positioning is not essential for repression by alpha2/Mcm1 in artificial reporter plasmids, and the importance of the nucleosome positioning remains questionable. We have tested the function of positioned nucleosomes through alteration of genomic chromatin at the a-cell-specific gene BAR1. We report here that a positioned nucleosome in the BAR1 promoter is disrupted in cis by the insertion of diverse DNA sequences such as poly(dA) . poly(dT) and poly(dC-dG) . poly(dC-dG), leading to inappropriate partial derepression of BAR1. Also, we show that isw2 mutation causes loss of nucleosome positioning in BAR1 in MATalpha cells as well as partial disruption of repression. Thus, nucleosome positioning is required for full repression, but loss of nucleosome positioning is not sufficient to relieve repression completely. Even though disruption of nucleosome positioning by the cis- and trans-acting modulators of chromatin has a modest effect on the level of transcription, it causes significant degradation of the alpha-mating pheromone in MATalpha cells, thereby affecting its cell type identity. Our results illustrate a useful paradigm for analysis of chromatin structural effects at genomic loci.

Non-B DNA conformations formed by long repeating tracts of myotonic dystrophy type 1, myotonic dystrophy type 2, and Friedreich's ataxia genes, not the sequences per se, promote mutagenesis in flanking regions

The expansions of long repeating tracts of CTG.CAG, CCTG.CAGG, and GAA.TTC are integral to the etiology of myotonic dystrophy type 1 (DM1), myotonic dystrophy type 2 (DM2), and Friedreich's ataxia (FRDA). Essentially all studies on the molecular mechanisms of this expansion process invoke an important role for non-B DNA conformations which may be adopted by these repeat sequences. We have directly evaluated the role(s) of the repeating sequences per se, or of the non-B DNA conformations formed by these sequences, in the mutagenic process. Studies in Escherichia coli and three types of mammalian (COS-7, CV-1, and HEK-293) fibroblast-like cells revealed that conditions which promoted the formation of the non-B DNA structures enhanced the genetic instabilities, both within the repeat sequences and in the flanking sequences of up to approximately 4 kbp. The three strategies utilized included: the in vivo modulation of global negative supercoil density using topA and gyrB mutant E. coli strains; the in vivo cleavage of hairpin loops, which are an obligate consequence of slipped-strand structures, cruciforms, and intramolecular triplexes, by inactivation of the SbcC protein; and by genetic instability studies with plasmids containing long repeating sequence inserts that do, and do not, adopt non-B DNA structures in vitro. Hence, non-B DNA conformations are critical for these mutagenesis mechanisms.

Sticky DNA: in vivo formation in E. coli and in vitro association of long GAA*TTC tracts to generate two independent supercoiled domains

The expanded GAA*TTC repeat sequence associated with Friedreich's ataxia (FRDA) adopts non-B DNA structures, (triplexes and sticky DNA). Sticky DNA is formed in plasmids by the association of two long GAA*TTC tracts at lengths that are found in the sequence of the frataxin gene in patients. Most FRDA patients have expanded GAA*TTC repeats (up to 1700 triplets), which inhibit the transcription of the gene, thus diminishing the synthesis of frataxin, a mitochondrial protein involved in iron-sulfur cluster biogenesis. Negative supercoiling and MgCl(2) (or MnCl(2)) are required to stabilize sticky DNA (a dumbbell-shaped structure) in plasmids with a pair of repeat tracts where n> or =60 in the direct repeat orientation in vitro. Since the triplet repeat sequences (TRS) were symmetrically positioned in the plasmids and because a number of unique restriction sites were present in the vector, studies were conducted to evaluate the influence of selectively linearizing one or the other supercoiled domains created by the DNA*DNA associated region, i.e. the stable complex at the pair of TRS's. The two domains behave independently, thus confirming the association of the two tracts and the dumbbell-shaped plasmid in our model for sticky DNA. Linking number investigations were performed on a family of plasmids harboring different lengths (30, 60, or 176 repeats), orientations and number of tracts (one or two) of a GAA*TTC repeat in Escherichia coli to evaluate the in vivo role, if any, of sticky DNA. Unexpectedly, this non-B DNA conformation elicited the formation of a TRS-length dependent change in the global topology of the plasmids, indicative of an apparent compression of the primary helices. Thus, linking number determinations confirm that sticky DNA has an important consequence in vivo.

Analysis of Non‐B DNA Structure at Chromosomal Sites in the Mammalian Genome

Changes at sites of genetic instability ultimately involve DNA repair pathways. Some sites of genetic instability in the mammalian genome appear to be unstable because they adopt a non-B DNA conformation. We describe two structural approaches for determination of whether a genomic region is configured in a non-B DNA conformation. Our studies indicate that at least some chromosomal fragile sites can be explained by such altered DNA conformations. One of the methods that we describe is called the bisulfite modification assay. This is a powerful assay because it provides information on individual DNA molecules. The second approach uses preexisting DNA structural reagents, but describes our specific application of them to analysis of DNA in vivo.

Thermodynamic basis for antibody binding to Z-DNA: comparison of a monoclonal antibody and its recombinant derivatives

Antibody engineering represents a promising area in biotechnology. Recombinant antibodies can be easily manipulated generating new ligand and effector activities that can be used as prototype magic bullets. On the other hand, an extensive knowledge of recombinant antibody binding and stability features are essential for an efficient substitution. In this study, we compared the stability and protein binding properties of two recombinant antibody fragments with their parental monoclonal antibody. The recombinant fragments were a monomeric scFv and a dimeric one, harboring human IgG1 CH2-CH3 domains. We have used fluorescence titration quenching to determine the thermodynamics of the interaction between an anti-Z-DNA monoclonal antibody and its recombinant antibody fragments with Z-DNA. All the antibody fragments seemed to bind DNA similarly, in peculiar two-affinity states. Enthalpy-entropy compensation was observed for both affinity states, but a marked entropy difference was observed for the monomeric scFv antibody fragment, mainly for the high affinity binding. In addition, we compared the stability of the dimeric antibody fragment and found differences favoring the monoclonal antibody. These differences seem to derive from the heterologous expression system used.

Increased negative superhelical density in vivo enhances the genetic instability of triplet repeat sequences

The influence of negative superhelical density on the genetic instabilities of long GAA.TTC, CGG.CCG, and CTG.CAG repeat sequences was studied in vivo in topologically constrained plasmids in Escherichia coli. These repeat tracts are involved in the etiologies of Friedreich ataxia, fragile X syndrome, and myotonic dystrophy type 1, respectively. The capacity of these DNA tracts to undergo deletions-expansions was explored with three genetic-biochemical approaches including first, the utilization of topoisomerase I and/or DNA gyrase mutants, second, the specific inhibition of DNA gyrase by novobiocin, and third, the genetic removal of the HU protein, thus lowering the negative supercoil density (-sigma). All three strategies revealed that higher -sigma in vivo enhanced the formation of deleted repeat sequences. The effects were most pronounced for the Friedreich ataxia and the fragile X triplet repeat sequences. Higher levels of -sigma stabilize non-B DNA conformations (i.e. triplexes, sticky DNA, flexible and writhed DNA, slipped structures) at appropriate repeat tracts; also, numerous prior genetic instability investigations invoke a role for these structures in promoting the slippage of the DNA complementary strands. Thus, we propose that the in vivo modulation of the DNA structure, localized to the repeat tracts, is responsible for these behaviors. Presuming that these interrelationships are also found in humans, dynamic alterations in the chromosomal nuclear matrix may modulate the -sigma of certain DNA regions and, thus, stabilize/destabilize certain non-B conformations which regulate the genetic expansions-deletions responsible for the diseases.

Advances in mechanisms of genetic instability related to hereditary neurological diseases

Substantial progress has been realized in the past several years in our understanding of the molecular mechanisms responsible for the expansions and deletions (genetic instabilities) of repeating tri-, tetra- and pentanucleotide repeating sequences associated with a number of hereditary neurological diseases. These instabilities occur by replication, recombination and repair processes, probably acting in concert, due to slippage of the DNA complementary strands relative to each other. The biophysical properties of the folded-back repeating sequence strands play a critical role in these instabilities. Non-B DNA structural elements (hairpins and slipped structures, DNA unwinding elements, tetraplexes, triplexes and sticky DNA) are described. The replication mechanisms are influenced by pausing of the replication fork, orientation of the repeat strands, location of the repeat sequences relative to replication origins and the flap endonuclease. Methyl-directed mismatch repair, nucleotide excision repair, and repair of damage caused by mutagens are discussed. Genetic recombination and double-strand break repair advances in Escherichia coli, yeast and mammalian models are reviewed. Furthermore, the newly discovered capacities of certain triplet repeat sequences to cause gross chromosomal rearrangements are discussed.

Topological insulators inhibit diffusion of transcription-induced positive supercoils in the chromosome of Escherichia coli

The double helical nature of DNA implies that progression of transcription machinery that cannot rotate easily around the DNA axis creates waves of positive supercoils ahead of it and negative supercoils behind it. Using topological reporters that detect local variations in DNA supercoiling, we have characterized the diffusion of transcription-induced (TI) positive supercoils in plasmids or in the chromosome of wild type Escherichia coli cells. Transcription-induced positive supercoils were able to diffuse and affect local supercoiling several kilobases away from the site of origin. By testing the effect of various DNA sequences, these reporters enabled us to identify elements that impede supercoil diffusion, i.e. behave as topological insulators. All the elements tested correspond to DNA gyrase catalytic targets. These results correlate the ability of a DNA sequence to be cleaved by DNA gyrase with topological insulator activity. Implications of the asymmetry in supercoil diffusion for the control of DNA topology are discussed.

The analysis of stress-induced duplex destabilization in long genomic DNA sequences

We present a method for calculating predicted locations and extents of stress-induced DNA duplex destabilization (SIDD) as functions of base sequence and stress level in long DNA molecules. The base pair denaturation energies are assigned individually, so the influences of near neighbors, methylated bases, adducts, or lesions can be included. Sample calculations indicate that copolymeric energetics give results that are close to those derived when full near-neighbor energetics are used; small but potentially informative differences occur only in the calculated SIDD properties of moderately destabilized regions. The method presented here for analyzing long sequences calculates the destabilization properties within windows of fixed length N, with successive windows displaced by an offset distance d(o). The final values of the relevant destabilization parameters for each base pair are calculated as weighted averages of the values computed for each window in which that base pair appears. This approach implicitly assumes that the strength of the direct coupling between remote base pairs that is induced by the imposed stress attenuates with their separation distance. This strategy enables calculations of the destabilization properties of DNA sequences of any length, up to and including complete chromosomes. We illustrate its utility by calculating the destabilization properties of the entire E. coli genomic DNA sequence. A preliminary analysis of the results shows that promoters are associated with SIDD regions in a highly statistically significant manner, suggesting that SIDD attributes may prove useful in the computational prediction of promoter locations in prokaryotes.

Distributions of Z-DNA and nuclear factor I in human chromosome 22: a model for coupled transcriptional regulation

An analysis of the human chromosome 22 genomic sequence shows that both Z-DNA forming regions (ZDRs) and promoter sites for nuclear factor-I (NFI) are correlated with the locations of known and predicted genes across the chromosome and accumulate around the transcriptional start sites of the known genes. Thus, the occurrence of Z-DNA across human genomic sequences mirrors that of a known eukaryotic transcription factor. In addition, 43 of the 383 fully annotated chromosomal genes have ZDRs within 2 nucleosomes upstream of strong NFIs. This suggests a distinct class of human genes that may potentially be transcriptionally regulated by a mechanism that couples Z-DNA with NFI activation, similar to the mechanism previously elucidated for the human colony stimulation factor-I promoter [Liu et al. (2001) Cell, 106, 309-318]. The results from this study will facilitate the design of experimental studies to test the generality of this mechanism for other genes in the cell.

Hairpin Structure-forming Propensity of the (CCTG·CAGG) Tetranucleotide Repeats Contributes to the Genetic Instability Associated with Myotonic Dystrophy Type 2

The genetic instabilities of (CCTG.CAGG)(n) tetranucleotide repeats were investigated to evaluate the molecular mechanisms responsible for the massive expansions found in myotonic dystrophy type 2 (DM2) patients. DM2 is caused by an expansion of the repeat from the normal allele of 26 to as many as 11,000 repeats. Genetic expansions and deletions were monitored in an African green monkey kidney cell culture system (COS-7 cells) as a function of the length (30, 114, or 200 repeats), orientation, or proximity of the repeat tracts to the origin (SV40) of replication. As found for CTG.CAG repeats related to DM1, the instabilities were greater for the longer tetranucleotide repeat tracts. Also, the expansions and deletions predominated when cloned in orientation II (CAGG on the leading strand template) rather than I and when cloned proximal rather than distal to the replication origin. Biochemical studies on synthetic d(CAGG)(26) and d(CCTG)(26) as models of unpaired regions of the replication fork revealed that d(CAGG)(26) has a marked propensity to adopt a defined base paired hairpin structure, whereas the complementary d(CCTG)(26) lacks this capacity. The effect of orientation described above differs from all previous results with three triplet repeat sequences (including CTG.CAG), which are also involved in the etiologies of other hereditary neurological diseases. However, similar to the triplet repeat sequences, the ability of one of the two strands to form a more stable folded structure, in our case the CAGG strand, explains this unorthodox "reversed" behavior.

Interaction of the Zalpha domain of human ADAR1 with a negatively supercoiled plasmid visualized by atomic force microscopy

Interest to the left-handed DNA conformation has been recently boosted by the findings that a number of proteins contain the Zalpha domain, which has been shown to specifically recognize Z-DNA. The biological function of Zalpha is presently unknown, but it has been suggested that it may specifically direct protein regions of Z-DNA induced by negative supercoiling in actively transcribing genes. Many studies, including a crystal structure in complex with Z-DNA, have focused on the human ADAR1 Zalpha domain in isolation. We have hypothesized that the recognition of a Z-DNA sequence by the Zalpha(ADAR1) domain is context specific, occurring under energetic conditions, which favor Z-DNA formation. To test this hypothesis, we have applied atomic force microscopy to image Zalpha(ADAR1) complexed with supercoiled plasmid DNAs. We have demonstrated that the Zalpha(ADAR1) binds specifically to Z-DNA and preferentially to d(CG)(n) inserts, which require less energy for Z-DNA induction compared to other sequences. A notable finding is that site-specific Zalpha binding to d(GC)(13) or d(GC)(2)C(GC)(10) inserts is observed when DNA supercoiling is insufficient to induce Z-DNA formation. These results indicate that Zalpha(ADAR1) binding facilities the B-to-Z transition and provides additional support to the model that Z-DNA binding proteins may regulate biological processes through structure-specific recognition.