Dual conformational recognition by Z-DNA binding protein is important for the B-Z transition process.

Left-handed Z-DNA is radically different from the most common right-handed B-DNA and can be stabilized by interactions with the Zα domain, which is found in a group of proteins, such as human ADAR1 and viral E3L proteins. It is well-known that most Zα domains bind to Z-DNA in a conformation-specific manner and induce rapid B-Z transition in physiological conditions. Although many structural and biochemical studies have identified the detailed interactions between the Zα domain and Z-DNA, little is known about the molecular basis of the B-Z transition process. In this study, we successfully converted the B-Z transition-defective Zα domain, vvZαE3L, into a B-Z converter by improving B-DNA binding ability, suggesting that B-DNA binding is involved in the B-Z transition. In addition, we engineered the canonical B-DNA binding protein GH5 into a Zα-like protein having both Z-DNA binding and B-Z transition activities by introducing Z-DNA interacting residues. Crystal structures of these mutants of vvZαE3L and GH5 complexed with Z-DNA confirmed the significance of conserved Z-DNA binding interactions. Altogether, our results provide molecular insight into how Zα domains obtain unusual conformational specificity and induce the B-Z transition.

Environmental Effects on Guanine-Thymine Mispair Tautomerization Explored with Quantum Mechanical/Molecular Mechanical Free Energy Simulations.

DNA bases can adopt energetically unfavorable tautomeric forms that enable the formation of Watson-Crick-like (WC-like) mispairs, which have been proposed to give rise to spontaneous mutations in DNA and misincorporation errors in DNA replication and translation. Previous NMR and computational studies have indicated that the population of WC-like guanine-thymine (G-T) mispairs depends on the environment, such as the local nucleic acid sequence and solvation. To investigate these environmental effects, herein G-T mispair tautomerization processes are studied computationally in aqueous solution, in A-form and B-form DNA duplexes, and within the active site of a DNA polymerase λ variant. The wobble G-T (wG-T), WC-like G-T*, and WC-like G*-T forms are considered, where * indicates the enol tautomer of the base. The minimum free energy paths for the tautomerization from the wG-T to the WC-like G-T* and from the WC-like G-T* to the WC-like G*-T are computed with mixed quantum mechanical/molecular mechanical (QM/MM) free energy simulations. The reaction free energies and free energy barriers are found to be significantly influenced by the environment. The wG-T→G-T* tautomerization is predicted to be endoergic in aqueous solution and the DNA duplexes but slightly exoergic in the polymerase, with Arg517 and Asn513 providing electrostatic stabilization of G-T*. The G-T*→G*-T tautomerization is also predicted to be slightly more thermodynamically favorable in the polymerase relative to these DNA duplexes. These simulations are consistent with an experimentally driven kinetic misincorporation model suggesting that G-T mispair tautomerization occurs in the ajar polymerase conformation or concertedly with the transition from the ajar to the closed polymerase conformation. Furthermore, the order of the associated two proton transfer reactions is predicted to be different in the polymerase than in aqueous solution and the DNA duplexes. These studies highlight the impact of the environment on the thermodynamics, kinetics, and fundamental mechanisms of G-T mispair tautomerization, which plays a role in a wide range of biochemically important processes.

The many lives of type IA topoisomerases.

The double-helical structure of genomic DNA is both elegant and functional in that it serves both to protect vulnerable DNA bases and to facilitate DNA replication and compaction. However, these design advantages come at the cost of having to evolve and maintain a cellular machinery that can manipulate a long polymeric molecule that readily becomes topologically entangled whenever it has to be opened for translation, replication, or repair. If such a machinery fails to eliminate detrimental topological entanglements, utilization of the information stored in the DNA double helix is compromised. As a consequence, the use of B-form DNA as the carrier of genetic information must have co-evolved with a means to manipulate its complex topology. This duty is performed by DNA topoisomerases, which therefore are, unsurprisingly, ubiquitous in all kingdoms of life. In this review, we focus on how DNA topoisomerases catalyze their impressive range of DNA-conjuring tricks, with a particular emphasis on DNA topoisomerase III (TOP3). Once thought to be the most unremarkable of topoisomerases, the many lives of these type IA topoisomerases are now being progressively revealed. This research interest is driven by a realization that their substrate versatility and their ability to engage in intimate collaborations with translocases and other DNA-processing enzymes are far more extensive and impressive than was thought hitherto. This, coupled with the recent associations of TOP3s with developmental and neurological pathologies in humans, is clearly making us reconsider their undeserved reputation as being unexceptional enzymes.

The complex formed between a synthetic RNA aptamer and the transcription repressor TetR is a structural and functional twin of the operator DNA-TetR regulator complex.

RNAs play major roles in the regulation of gene expression. Hence, designer RNA molecules are increasingly explored as regulatory switches in synthetic biology. Among these, the TetR-binding RNA aptamer was selected by its ability to compete with operator DNA for binding to the bacterial repressor TetR. A fortuitous finding was that induction of TetR by tetracycline abolishes both RNA aptamer and operator DNA binding in TetR. This enabled numerous applications exploiting both the specificity of the RNA aptamer and the efficient gene repressor properties of TetR. Here, we present the crystal structure of the TetR-RNA aptamer complex at 2.7 Å resolution together with a comprehensive characterization of the TetR-RNA aptamer versus TetR-operator DNA interaction using site-directed mutagenesis, size exclusion chromatography, electrophoretic mobility shift assays and isothermal titration calorimetry. The fold of the RNA aptamer bears no resemblance to regular B-DNA, and neither does the thermodynamic characterization of the complex formation reaction. Nevertheless, the functional aptamer-binding epitope of TetR is fully contained within its DNA-binding epitope. In the RNA aptamer complex, TetR adopts the well-characterized DNA-binding-competent conformation of TetR, thus revealing how the synthetic TetR-binding aptamer strikes the chords of the bimodal allosteric behaviour of TetR to function as a synthetic regulator.

Enrichment technique to allow early detection and monitor emergence of KRAS mutation in response to treatment.

Sensitivity of cell-free circulating tumour DNA (ctDNA) assays is often hampered by the limited quantity of intact mutant nucleotide fragments. To overcome the issue of substrate limitation in clinical applications, we developed an enrichment method utilizing pyrrole-imidazole (PI) polyamides and their ability to bind the minor groove of B-DNA. We present here a proof-of-concept experiment to enrich specific mutant KRAS alleles with biotinylated PI polyamides. We investigated the clinical feasibility of incorporating PI polyamides to detect KRAS mutations in ctDNA from 40 colorectal cancer (CRC) patients, of whom 17 carried mutations in KRAS. After enriching ctDNA with those polyamides, we used digital PCR to detect several common KRAS codon 12 mutations. Enrichment by biotinylated PI polyamides improved the sensitivity of ctDNA analysis (88.9% vs. 11.1%, P < 0.01) in 9 non-metastatic mutation-positive patients. We observed no differences in performance for the 8 metastatic subjects (100% vs. 75%, P = 0.47). In the remaining 23/40 patients with wild type KRAS codon 12, no mutant alleles were detected with or without polyamide-facilitated enrichment. Enriching B-form of ctDNA with PI polyamides significantly improved the assay sensitivity in detecting KRAS mutations in non-metastatic CRC patient samples.

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 MgCl2 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.

Structural Investigation into the Threading Intercalation of a Chiral Dinuclear Ruthenium(II) Polypyridyl Complex through a B-DNA Oligonucleotide.

Herein we report the separation of the three stereoisomers of the DNA light-switch compound [{Ru(bpy)2}2(tpphz)]4+ (tpphz = tetrapyrido[3,2-a:2',3'-c:3″,2″-h:2‴,3‴-j]phenazine) by column chromatography and the characterization of each stereoisomer by X-ray crystallography. The interaction of these compounds with a DNA octanucleotide d(GCATATCG).d(CGATATGC) has been studied using NMR techniques. Selective deuteration of the bipyridyl rings was needed to provide sufficient spectral resolution to characterize structures. NMR-derived structures for these complexes show a threading intercalation binding mode with slow and chirality-dependent rates. This represents the first solution structure of an intercalated bis-ruthenium ligand. Intriguingly, we find that the binding site selectivity is dependent on the nature of the stereoisomer employed, with Λ RuII centers showing a better intercalation fit.

DNA Base Pair Mismatches Induce Structural Changes and Alter the Free-Energy Landscape of Base Flip.

Double-stranded DNA may contain mismatched base pairs beyond the Watson-Crick pairs, guanine-cytosine and adenine-thymine. Such mismatches bear adverse consequences for human health. We utilize molecular dynamics and metadynamics computer simulations to study the structure and dynamics of both matched and mismatched base pairs. We discover significant differences between the matched and mismatched pairs in structure and base flip work profiles. Mismatched pairs shift more in the plane normal to the DNA strand and exhibit more noncanonical structures, including the e-motif. We discuss the potential implications on the mismatch repair enzymes' detection of DNA mismatches.

Sequence preference and structural heterogeneity of BZ junctions.

BZ junctions, which connect B-DNA to Z-DNA, are necessary for local transformation of B-DNA to Z-DNA in the genome. However, the limited information on the junction-forming sequences and junction structures has led to a lack of understanding of the structural diversity and sequence preferences of BZ junctions. We determined three crystal structures of BZ junctions with diverse sequences followed by spectroscopic validation of DNA conformation. The structural features of the BZ junctions were well conserved regardless of sequences via the continuous base stacking through B-to-Z DNA with A-T base extrusion. However, the sequence-dependent structural heterogeneity of the junctions was also observed in base step parameters that are correlated with steric constraints imposed during Z-DNA formation. Further, circular dichroism and fluorescence-based analysis of BZ junctions revealed that a base extrusion was only found at the A-T base pair present next to a stable dinucleotide Z-DNA unit. Our findings suggest that Z-DNA formation in the genome is influenced by the sequence preference for BZ junctions.

Toward a Universal Structural and Energetic Model for Prokaryotic Promoters.

With almost no consensus promoter sequence in prokaryotes, recruitment of RNA polymerase (RNAP) to precise transcriptional start sites (TSSs) has remained an unsolved puzzle. Uncovering the underlying mechanism is critical for understanding the principle of gene regulation. We attempted to search the hidden code in ∼16,500 promoters of 12 prokaryotes representing two kingdoms in their structure and energetics. Twenty-eight fundamental parameters of DNA structure including backbone angles, basepair axis, and interbasepair and intrabasepair parameters were used, and information was extracted from x-ray crystallography data. Three parameters (solvation energy, hydrogen-bond energy, and stacking energy) were selected for creating energetics profiles using in-house programs. DNA of promoter regions was found to be inherently designed to undergo a change in every parameter undertaken for the study, in all prokaryotes. The change starts from some distance upstream of TSSs and continues past some distance from TSS, hence giving a signature state to promoter regions. These signature states might be the universal hidden codes recognized by RNAP. This observation was reiterated when randomly selected promoter sequences (with little sequence conservation) were subjected to structure generation; all developed into very similar three-dimensional structures quite distinct from those of conventional B-DNA and coding sequences. Fine structural details at important motifs (viz. -11, -35, and -75 positions relative to TSS) of promoters reveal novel to our knowledge and pointed insights for RNAP interaction at these locations; it could be correlated with how some particular structural changes at the -11 region may allow insertion of RNAP amino acids in interbasepair space as well as facilitate the flipping out of bases from the DNA duplex.

Noncanonical secondary structures arising from non-B DNA motifs are determinants of mutagenesis.

Somatic mutations show variation in density across cancer genomes. Previous studies have shown that chromatin organization and replication time domains are correlated with, and thus predictive of, this variation. Here, we analyze 1809 whole-genome sequences from 10 cancer types to show that a subset of repetitive DNA sequences, called non-B motifs that predict noncanonical secondary structure formation can independently account for variation in mutation density. Combined with epigenetic factors and replication timing, the variance explained can be improved to 43%-76%. Approximately twofold mutation enrichment is observed directly within non-B motifs, is focused on exposed structural components, and is dependent on physical properties that are optimal for secondary structure formation. Therefore, there is mounting evidence that secondary structures arising from non-B motifs are not simply associated with increased mutation density-they are possibly causally implicated. Our results suggest that they are determinants of mutagenesis and increase the likelihood of recurrent mutations in the genome. This analysis calls for caution in the interpretation of recurrent mutations and highlights the importance of taking non-B motifs that can simply be inferred from the reference sequence into consideration in background models of mutability henceforth.

6S RNA Mimics B-Form DNA to Regulate Escherichia coli RNA Polymerase.

Noncoding RNAs (ncRNAs) regulate gene expression in all organisms. Bacterial 6S RNAs globally regulate transcription by binding RNA polymerase (RNAP) holoenzyme and competing with promoter DNA. Escherichia coli (Eco) 6S RNA interacts specifically with the housekeeping σ70-holoenzyme (Eσ70) and plays a key role in the transcriptional reprogramming upon shifts between exponential and stationary phase. Inhibition is relieved upon 6S RNA-templated RNA synthesis. We report here the 3.8 Å resolution structure of a complex between 6S RNA and Eσ70 determined by single-particle cryo-electron microscopy and validation of the structure using footprinting and crosslinking approaches. Duplex RNA segments have A-form C3' endo sugar puckers but widened major groove widths, giving the RNA an overall architecture that mimics B-form promoter DNA. Our results help explain the specificity of Eco 6S RNA for Eσ70 and show how an ncRNA can mimic B-form DNA to directly regulate transcription by the DNA-dependent RNAP.

A B-Z junction induced by an A A mismatch in GAC repeats in the gene for cartilage oligomeric matrix protein promotes binding with the hZαADAR1 protein.

GAC repeat expansion from five to seven in the exonic region of the gene for cartilage oligomeric matrix protein (COMP) leads to pseudoachondroplasia, a skeletal abnormality. However, the molecular mechanism by which GAC expansions in the COMP gene lead to skeletal dysplasias is poorly understood. Here we used molecular dynamics simulations, which indicate that an A … A mismatch in a d(GAC)6·d(GAC)6 duplex induces negative supercoiling, leading to a local B-to-Z DNA transition. This transition facilitates the binding of d(GAC)7·d(GAC)7 with the Zα-binding domain of human adenosine deaminase acting on RNA 1 (ADAR1, hZαADAR1), as confirmed by CD, NMR, and microscale thermophoresis studies. The CD results indicated that hZαADAR1 recognizes the zigzag backbone of d(GAC)7·d(GAC)7 at the B-Z junction and subsequently converts it into Z-DNA via the so-called passive mechanism. Molecular dynamics simulations carried out for the modeled hZαADAR1-d(GAC)6d(GAC)6 complex confirmed the retention of previously reported important interactions between the two molecules. These findings suggest that hZαADAR1 binding with the GAC hairpin stem in COMP can lead to a non-genetic, RNA editing-mediated substitution in COMP that may then play a crucial role in the development of pseudoachondroplasia.

Effect of gold nanoparticle on stability of the DNA molecule: A study of molecular dynamics simulation.

An understanding of the mechanism of DNA interactions with gold nanoparticles is useful in today medicine applications. We have performed a molecular dynamics simulation on a B-DNA duplex (CCTCAGGCCTCC) in the vicinity of a gold nanoparticle with a truncated octahedron structure composed of 201 gold atoms (diameter ∼1.8 nm) to investigate gold nanoparticle (GNP) effects on the stability of DNA. During simulation, the nanoparticle is closed to DNA and phosphate groups direct the particles into the major grooves of the DNA molecule. Because of peeling and untwisting states that are occur at end of DNA, the nucleotide base lies flat on the surface of GNP. The configuration entropy is estimated using the covariance matrix of atom-positional fluctuations for different bases. The results show that when a gold nanoparticle has interaction with DNA, entropy increases. The results of conformational energy and the hydrogen bond numbers for DNA indicated that DNA becomes unstable in the vicinity of a gold nanoparticle. The radial distribution function was calculated for water hydrogen-phosphate oxygen pairs. Almost for all nucleotide, the presence of a nanoparticle around DNA caused water molecules to be released from the DNA duplex and cations were close to the DNA.

Surface-assisted DNA self-assembly: An enzyme-free strategy towards formation of branched DNA lattice.

DNA based self-assembled nanostructures and DNA origami has proven useful for organizing nanomaterials with firm precision. However, for advanced applications like nanoelectronics and photonics, large-scale organization of self-assembled branched DNA (bDNA) into periodic lattices is desired. In this communication for the first time we report a facile method of self-assembly of Y-shaped bDNA nanostructures on the cationic surface of Aluminum (Al) foil to prepare periodic two dimensional (2D) bDNA lattice. Particularly those Y-shaped bDNA structures having smaller overhangs and unable to self-assemble in solution, they are easily assembled on the surface of Al foil in the absence of ligase. Field emission scanning electron microscopy (FESEM) analysis shows homogenous distribution of two-dimensional bDNA lattices across the Al foil. When the assembled bDNA structures were recovered from the Al foil and electrophoresed in nPAGE only higher order polymeric bDNA structures were observed without a trace of monomeric structures which confirms the stability and high yield of the bDNA lattices. Therefore, this enzyme-free economic and efficient strategy for developing bDNA lattices can be utilized in assembling various nanomaterials for functional molecular components towards development of DNA based self-assembled nanodevices.

Influence of BII Backbone Substates on DNA Twist: A Unified View and Comparison of Simulation and Experiment for All 136 Distinct Tetranucleotide Sequences.

Reliable representation of the B-DNA base-pair step twist is one of the crucial requirements for theoretical modeling of DNA supercoiling and other biologically relevant phenomena in B-DNA. It has long been suspected that the twist is inaccurately described by current empirical force fields. Unfortunately, comparison of simulation results with experiments is not straightforward because of the presence of BII backbone substates, whose populations may differ in experimental and simulation ensembles. In this work, we provide a comprehensive view of the effect of BII substates on the overall B-DNA helix twist and show how to reliably compare twist values from experiment and simulation in two scenarios. First, for longer DNA segments freely moving in solution, we show that sequence-averaged twists of different BI/BII ensembles can be compared directly because of approximate cancellation of the opposing BII effects. Second, for sequence-specific data, such as a particular base-pair step or tetranucleotide twist, can be compared only for a clearly defined BI/BII backbone conformation. For the purpose of force field testing, we designed a compact set of fourteen 22-base-pair B-DNA duplexes (Set 14) containing all 136 distinct tetranucleotide sequences and carried out a total of 84 µs of molecular dynamics simulations, primarily with the OL15 force field. Our results show that the ff99bsc0εζOL1χOL4, parmbsc1, and OL15 force fields model the B-DNA helical twist in good agreement with X-ray and minicircle ligation experiments. The comprehensive understanding obtained regarding the effect of BII substates on the base-pair step geometry should aid meaningful comparisons of various conformational ensembles in future research.

2'-Alkynylnucleotides: A Sequence- and Spin Label-Flexible Strategy for EPR Spectroscopy in DNA.

Electron paramagnetic resonance (EPR) spectroscopy is a powerful method to elucidate molecular structure through the measurement of distances between conformationally well-defined spin labels. Here we report a sequence-flexible approach to the synthesis of double spin-labeled DNA duplexes, where 2'-alkynylnucleosides are incorporated at terminal and internal positions on complementary strands. Post-DNA synthesis copper-catalyzed azide-alkyne cycloaddition (CuAAC) reactions with a variety of spin labels enable the use of double electron-electron resonance experiments to measure a number of distances on the duplex, affording a high level of detailed structural information.

Synthesis, base pairing and structure studies of geranylated RNA.

Natural RNAs utilize extensive chemical modifications to diversify their structures and functions. 2-Thiouridine geranylation is a special hydrophobic tRNA modification that has been discovered very recently in several bacteria, such as Escherichia coli, Enterobacter aerogenes, Pseudomonas aeruginosa and Salmonella Typhimurium The geranylated residues are located in the first anticodon position of tRNAs specific for lysine, glutamine and glutamic acid. This big hydrophobic terpene functional group affects the codon recognition patterns and reduces frameshifting errors during translation. We aimed to systematically study the structure, function and biosynthesis mechanism of this geranylation pathway, as well as answer the question of why nature uses such a hydrophobic modification in hydrophilic RNA systems. Recently, we have synthesized the deoxy-analog of S-geranyluridine and showed the geranylated T-G pair is much stronger than the geranylated T-A pair and other mismatched pairs in the B-form DNA duplex context, which is consistent with the observation that the geranylated tRNA(Glu) UUC recognizes GAG more efficiently than GAA. In this manuscript we report the synthesis and base pairing specificity studies of geranylated RNA oligos. We also report extensive molecular simulation studies to explore the structural features of the geranyl group in the context of A-form RNA and its effect on codon-anticodon interaction during ribosome binding.

Translocation and deletion breakpoints in cancer genomes are associated with potential non-B DNA-forming sequences.

Gross chromosomal rearrangements (including translocations, deletions, insertions and duplications) are a hallmark of cancer genomes and often create oncogenic fusion genes. An obligate step in the generation of such gross rearrangements is the formation of DNA double-strand breaks (DSBs). Since the genomic distribution of rearrangement breakpoints is non-random, intrinsic cellular factors may predispose certain genomic regions to breakage. Notably, certain DNA sequences with the potential to fold into secondary structures [potential non-B DNA structures (PONDS); e.g. triplexes, quadruplexes, hairpin/cruciforms, Z-DNA and single-stranded looped-out structures with implications in DNA replication and transcription] can stimulate the formation of DNA DSBs. Here, we tested the postulate that these DNA sequences might be found at, or in close proximity to, rearrangement breakpoints. By analyzing the distribution of PONDS-forming sequences within ±500 bases of 19 947 translocation and 46 365 sequence-characterized deletion breakpoints in cancer genomes, we find significant association between PONDS-forming repeats and cancer breakpoints. Specifically, (AT)n, (GAA)n and (GAAA)n constitute the most frequent repeats at translocation breakpoints, whereas A-tracts occur preferentially at deletion breakpoints. Translocation breakpoints near PONDS-forming repeats also recur in different individuals and patient tumor samples. Hence, PONDS-forming sequences represent an intrinsic risk factor for genomic rearrangements in cancer genomes.

Premeltons in DNA.

Premeltons are examples of emergent-structures (i.e., structural-solitons) that arise spontaneously in DNA due to the presence of nonlinear-excitations in its structure. They are of two kinds: B-B (or A-A) premeltons form at specific DNA-regions to nucleate site-specific DNA melting. These are stationary and, being globally-nontopological, undergo breather-motions that allow drugs and dyes to intercalate into DNA. B-A (or A-B) premeltons, on the other hand, are mobile, and being globally-topological, act as phase-boundaries transforming B- into A-DNA during the structural phase-transition. They are not expected to undergo breather motions. A key feature of both types of premeltons is the presence of an intermediate structural-form in their central regions (proposed as being a transition-state intermediate in DNA-melting and in the B- to A-transition), which differs from either A- or B-DNA. Called beta-DNA, this is both metastable and hyperflexible--and contains an alternating sugar-puckering pattern along the polymer backbone combined with the partial unstacking (in its lower energy-forms) of every-other base-pair. Beta-DNA is connected to either B- or to A-DNA on either side by boundaries possessing a gradation of nonlinear structural-change, these being called the kink and the antikink regions. The presence of premeltons in DNA leads to a unifying theory to understand much of DNA physical chemistry and molecular biology. In particular, premeltons are predicted to define the 5' and 3' ends of genes in naked-DNA and DNA in active-chromatin, this having important implications for understanding physical aspects of the initiation, elongation and termination of RNA-synthesis during transcription. For these and other reasons, the model will be of broader interest to the general-audience working in these areas. The model explains a wide variety of data, and carries with it a number of experimental predictions--all readily testable--as will be described in this review.