Genomic Instability of G-Quadruplex Sequences in Escherichia coli: Roles of DinG, RecG, and RecQ Helicases.

Guanine-rich DNA can fold into highly stable four-stranded DNA structures called G-quadruplexes (G4). Originally identified in sequences from telomeres and oncogene promoters, they can alter DNA metabolism. Indeed, G4-forming sequences represent obstacles for the DNA polymerase, with important consequences for cell life as they may lead to genomic instability. To understand their role in bacterial genomic instability, different G-quadruplex-forming repeats were cloned into an Escherichia coli genetic system that reports frameshifts and complete or partial deletions of the repeat when the G-tract comprises either the leading or lagging template strand during replication. These repeats formed stable G-quadruplexes in single-stranded DNA but not naturally supercoiled double-stranded DNA. Nevertheless, transcription promoted G-quadruplex formation in the resulting R-loop for (G3T)4 and (G3T)8 repeats. Depending on genetic background and sequence propensity for structure formation, mutation rates varied by five orders of magnitude. Furthermore, while in vitro approaches have shown that bacterial helicases can resolve G4, it is still unclear whether G4 unwinding is important in vivo. Here, we show that a mutation in recG decreased mutation rates, while deficiencies in the structure-specific helicases DinG and RecQ increased mutation rates. These results suggest that G-quadruplex formation promotes genetic instability in bacteria and that helicases play an important role in controlling this process in vivo.

NACDDB: Nucleic Acid Circular Dichroism Database.

The Nucleic Acid Circular Dichroism Database (NACDDB) is a public repository that archives and freely distributes circular dichroism (CD) and synchrotron radiation CD (SRCD) spectral data about nucleic acids, and the associated experimental metadata, structural models, and links to literature. NACDDB covers CD data for various nucleic acid molecules, including DNA, RNA, DNA/RNA hybrids, and various nucleic acid derivatives. The entries are linked to primary sequence and experimental structural data, as well as to the literature. Additionally, for all entries, 3D structure models are provided. All entries undergo expert validation and curation procedures to ensure completeness, consistency, and quality of the data included. The NACDDB is open for submission of the CD data for nucleic acids. NACDDB is available at: https://genesilico.pl/nacddb/.

Crucial Role of the C-Terminal Domain of Hfq Protein in Genomic Instability.

G-rich DNA repeats that can form G-quadruplex structures are prevalent in bacterial genomes and are frequently associated with regulatory regions of genes involved in virulence, antigenic variation, and antibiotic resistance. These sequences are also inherently mutagenic and can lead to changes affecting cell survival and adaptation. Transcription of the G-quadruplex-forming repeat (G3T)n in E. coli, when mRNA comprised the G-rich strand, promotes G-quadruplex formation in DNA and increases rates of deletion of G-quadruplex-forming sequences. The genomic instability of G-quadruplex repeats may be a source of genetic variability that can influence alterations and evolution of bacteria. The DNA chaperone Hfq is involved in the genetic instability of these G-quadruplex sequences. Inactivation of the hfq gene decreases the genetic instability of G-quadruplex, demonstrating that the genomic instability of this regulatory element can be influenced by the E. coli highly pleiotropic Hfq protein, which is involved in small noncoding RNA regulation pathways, and DNA organization and packaging. We have shown previously that the protein binds to and stabilizes these sequences, increasing rates of their genomic instability. Here, we extend this analysis to characterize the role of the C-terminal domain of Hfq protein in interaction with G-quadruplex structures. This allows to better understand the function of this specific region of the Hfq protein in genomic instability.

Role of Hfq in Genome Evolution: Instability of G-Quadruplex Sequences in E. coli.

Certain G-rich DNA repeats can form quadruplex in bacterial chromatin that can present blocks to DNA replication and, if not properly resolved, may lead to mutations. To understand the participation of quadruplex DNA in genomic instability in Escherichia coli (E. coli), mutation rates were measured for quadruplex-forming DNA repeats, including (G3T)4, (G3T)8, and a RET oncogene sequence, cloned as the template or nontemplate strand. We evidence that these alternative structures strongly influence mutagenesis rates. Precisely, our results suggest that G-quadruplexes form in E. coli cells, especially during transcription when the G-rich strand can be displaced by R-loop formation. Structure formation may then facilitate replication misalignment, presumably associated with replication fork blockage, promoting genomic instability. Furthermore, our results also evidence that the nucleoid-associated protein Hfq is involved in the genetic instability associated with these sequences. Hfq binds and stabilizes G-quadruplex structure in vitro and likely in cells. Collectively, our results thus implicate quadruplexes structures and Hfq nucleoid protein in the potential for genetic change that may drive evolution or alterations of bacterial gene expression.

Novel computational study on π-stacking to understand mechanistic interactions of Tryptanthrin analogues with DNA.

Based on recently published initial experimental results on the intercalation of a class of broad spectrum antiparasitic compounds, we present a purely theoretical approach for determining if these compounds may preferentially intercalate with guanosine/cytosine (GC)-rich or adenosine/thymidine (TA)-rich regions of DNA. The predictive model presented herein is based upon utilization of density functional theory (DFT) to determine a priori how the best intercalator may energetically and sterically interact with each of the nucleoside base pairs. A potential new method using electrostatic potential maps (EPMs) to visually select the best poses is introduced and compared to the existing brute-force center of mass (COM) approach. The EPM and COM predictions are in agreement with each other, but the EPM method is potentially much more efficient. We report that 4-azatryptantrin, the best intercalator, is predicted to favor π-stacking with GC over that of TA by approximately 2-4 kcal/mol. This represents a significant difference if one takes into account the Boltzmann distribution at physiological temperature. This theoretical method will be utilized to guide future experimental studies on the elucidation of possible mechanism(s) for the action of these antiparasitic compounds at the molecular level.

Replication fork stalling and checkpoint activation by a PKD1 locus mirror repeat polypurine-polypyrimidine (Pu-Py) tract.

DNA sequences prone to forming noncanonical structures (hairpins, triplexes, G-quadruplexes) cause DNA replication fork stalling, activate DNA damage responses, and represent hotspots of genomic instability associated with human disease. The 88-bp asymmetric polypurine-polypyrimidine (Pu-Py) mirror repeat tract from the human polycystic kidney disease (PKD1) intron 21 forms non-B DNA secondary structures in vitro. We show that the PKD1 mirror repeat also causes orientation-dependent fork stalling during replication in vitro and in vivo. When integrated alongside the c-myc replicator at an ectopic chromosomal site in the HeLa genome, the Pu-Py mirror repeat tract elicits a polar replication fork barrier. Increased replication protein A (RPA), Rad9, and ataxia telangiectasia- and Rad3-related (ATR) checkpoint protein binding near the mirror repeat sequence suggests that the DNA damage response is activated upon replication fork stalling. Moreover, the proximal c-myc origin of replication was not required to cause orientation-dependent checkpoint activation. Cells expressing the replication fork barrier display constitutive Chk1 phosphorylation and continued growth, i.e. checkpoint adaptation. Excision of the Pu-Py mirror repeat tract abrogates the DNA damage response. Adaptation to Chk1 phosphorylation in cells expressing the replication fork barrier may allow the accumulation of mutations that would otherwise be remediated by the DNA damage response.

On an unbiased and consistent estimator for mutation rates.

Spontaneous mutations are stochastic events. The mutation rate, defined as mutations per genome per replication, is generally very low, and it is widely accepted that spontaneous mutations occur at defined, but different, rates in bacteriophage and in bacterial, insect, and mammalian cells. The calculation of mutation rates has proved to be a significant problem. Mutation rates can be calculated by following mutant accumulation during growth or from the distribution of mutants obtained in parallel cultures. As Luria and Delbrück described in 1943, the number of mutants in parallel populations of bacterial cells varies widely depending on when a spontaneous mutation occurs during growth of the culture. Since 1943, many mathematical refinements to estimating rates, called estimators, have been described to facilitate determination of the mutation rate from the distribution or frequency of mutants detected following growth of parallel cultures. We present a rigorous mathematical solution to the mutation rate problem using an unbiased and consistent estimator. Using this estimator we demonstrate experimentally that mutation rates can be easily calculated by determining mutant accumulation, that is, from the number of mutants measured in two successive generations. Moreover, to verify the consistency of our estimator we conduct a series of simulation trials that show a surprisingly rapid convergence to the targeted mutation rate (reached between 25th and 30th generations).

Replication-dependent instability at (CTG) x (CAG) repeat hairpins in human cells.

Instability of (CTG) x (CAG) microsatellite trinucleotide repeat (TNR) sequences is responsible for more than a dozen neurological or neuromuscular diseases. TNR instability during DNA synthesis is thought to involve slipped-strand or hairpin structures in template or nascent DNA strands, although direct evidence for hairpin formation in human cells is lacking. We have used targeted recombination to create a series of isogenic HeLa cell lines in which (CTG) x (CAG) repeats are replicated from an ectopic copy of the Myc (also known as c-myc) replication origin. In this system, the tendency of chromosomal (CTG) x (CAG) tracts to expand or contract was affected by origin location and the leading or lagging strand replication orientation of the repeats, and instability was enhanced by prolonged cell culture, increased TNR length and replication inhibition. Hairpin cleavage by synthetic zinc finger nucleases in these cells has provided the first direct evidence for the formation of hairpin structures during replication in vivo.

Antimicrobial activity of tryptanthrins in Escherichia coli.

Tryptanthrins have potential therapeutic activity against a wide variety of pathogenic organisms, although little is known about their mechanism. Activity against Escherichia coli, however, has not been examined. The effects of tryptanthrin (indolo[2,1-b]quinazolin-6,12-dione) and nine derivatives on growth, survival, and mutagenesis in E. coli were examined. Analogues with a nitrogen atom at the 4-position of tryptanthrin stopped log phase growth of E. coli cultures at concentrations as low as 5 microM. Tryptanthrins decreased viability during incubation with cells in buffer by factors of 10(-2) to <10(-6) at 0.2-40 microM. Derivatives with an oxime group at the 6-position exhibited the greatest bactericidal activity. Most tryptanthrins were not mutagenic in several independent assays, although the 4-aza and 4 aza-8-fluoro derivatives increased frameshift mutations about 22- and 4-fold, respectively. Given the structure of trypanthrins, binding to DNA may occur by intercalation. From analysis using a sensitive linking number assay, several tryptanthrins, especially the 4-aza and 6-oximo derivatives, intercalate into DNA.

Genetic instabilities of (CCTG).(CAGG) and (ATTCT).(AGAAT) disease-associated repeats reveal multiple pathways for repeat deletion.

The DNA repeats (CTG).(CAG), (CGG).(CCG), (GAA).(TTC), (ATTCT).(AGAAT), and (CCTG).(CAGG), undergo expansion in humans leading to neurodegenerative disease. A genetic assay for repeat instability has revealed that the activities of RecA and RecB during replication restart are involved in a high rate of deletion of (CTG).(CAG) repeats in E. coli. This assay has been applied to (CCTG).(CAGG) repeats associated with myotonic dystrophy type 2 (DM2) that expand to 11 000 copies and to spinocerebellar ataxia type 10 (SCA10) (ATTCT).(AGAAT) repeats that expand to 4500 copies in affected individuals. DM2 (CCTG).(CAGG) repeats show a moderate rate of instability, less than that observed for the myotonic dystrophy type 1 (CTG).(CAG) repeats, while the SCA10 (ATTCT).(AGAAT) repeats were remarkably stable in E. coli. In contrast to (CTG).(CAG) repeats, deletions of the DM2 and SCA10 repeats were not dependent on RecA and RecB, suggesting that replication restart may not be a predominant mechanism by which these repeats undergo deletion. These results suggest that different molecular mechanisms, or pathways, are responsible for the instability of different disease-associated DNA repeats in E. coli. These pathways involve simple replication slippage and various sister strand exchange events leading to deletions or expansions, often associated with plasmid dimerization. The differences in the mechanisms of repeat deletion may result from the differential propensity of these repeats to form various DNA secondary structures and their differential proclivity for primer-template misalignment during replication.

A Z-DNA sequence reduces slipped-strand structure formation in the myotonic dystrophy type 2 (CCTG) x (CAGG) repeat.

All DNA repeats known to undergo expansion leading to human neurodegenerative disease can form one, or several, alternative conformations, including hairpin, slipped strand, triplex, quadruplex, or unwound DNA structures. These alternative structures may interfere with the normal cellular processes of transcription, DNA repair, replication initiation, or polymerase elongation and thereby contribute to the genetic instability of these repeat tracts. We show that (CCTG) x (CAGG) repeats, in the first intron of the ZNF9 gene associated with myotonic dystrophy type 2, form slipped-strand DNA structures in a length-dependent fashion upon reduplexing. The threshold for structure formation on reduplexing is between 36 and 42 repeats in length. Alternative DNA structures also form in (CCTG)(58) x (CAGG)(58) and larger repeat tracts in plasmids at physiological superhelical densities. This represents an example of a sequence that forms slipped-strand DNA from the energy of DNA supercoiling. Moreover, Z-DNA forms in a (TG) x (CA) tract within the complex repeat sequence 5' of the (CCTG)(n) x (CAGG)(n) repeat in the ZNF9 gene. Upon reduplexing, the presence of the flanking sequence containing the Z-DNA-forming tract reduced the extent of slipped-strand DNA formation by 62% for (CCTG)(57) x (CAGG)(57) compared with 58 pure repeats without the flanking sequence. This finding suggests that the Z-DNA-forming sequence in the DM2 gene locus may have a protective effect of reducing the potential for slipped-strand DNA formation in (CCTG)(n) x (CAGG)(n) repeats.

Unstable spinocerebellar ataxia type 10 (ATTCT*(AGAAT) repeats are associated with aberrant replication at the ATX10 locus and replication origin-dependent expansion at an ectopic site in human cells.

Spinocerebellar ataxia type 10 (SCA10) is associated with expansion of (ATTCT)n repeats (where n is the number of repeats) within the ataxin 10 (ATX10/E46L) gene. The demonstration that (ATTCT)n tracts can act as DNA unwinding elements (DUEs) in vitro has suggested that aberrant replication origin activity occurs at expanded (ATTCT)n tracts and may lead to their instability. Here, we confirm these predictions. The wild-type ATX10 locus displays inefficient origin activity, but origin activity is elevated at the expanded ATX10 loci in patient-derived cells. To test whether (ATTCT)n tracts can potentiate origin activity, cell lines were constructed that contain ectopic copies of the c-myc replicator in which the essential DUE was replaced by ATX10 DUEs with (ATTCT)n. ATX10 DUEs containing (ATTCT)27 or (ATTCT)48, but not (ATTCT)8 or (ATTCT)13, could substitute functionally for the c-myc DUE, but (ATTCT)48 could not act as an autonomous replicator. Significantly, chimeric c-myc replicators containing ATX10 DUEs displayed length-dependent (ATTCT)n instability. By 250 population doublings, dramatic two- and fourfold length expansions were observed for (ATTCT)27 and (ATTCT)48 but not for (ATTCT)8 or (ATTCT)13. These results implicate replication origin activity as one molecular mechanism associated with the instability of (ATTCT)n tracts that are longer than normal length.

Slipped strand DNA structures.

Slipped strand DNA structures are formed when complementary strands comprising direct repeats pair in a misaligned, or slipped, fashion along the DNA helix axis. Although slipped strand DNA may form in almost any direct repeat, to date, these structures have only been detected in short DNA repeats, termed unstable DNA repeats, in which expansion is associated with many neurodegenerative diseases. This alternative DNA structure, or a similar slipped intermediate DNA that may form during DNA replication or repair, may be a causative factor in the instability of the DNA sequences that can form these structures.

Target DNA structure plays a critical role in RAG transposition.

Antigen receptor gene rearrangements are initiated by the RAG1/2 protein complex, which recognizes specific DNA sequences termed RSS (recombination signal sequences). The RAG recombinase can also catalyze transposition: integration of a DNA segment bounded by RSS into an unrelated DNA target. For reasons that remain poorly understood, such events occur readily in vitro, but are rarely detected in vivo. Previous work showed that non-B DNA structures, particularly hairpins, stimulate transposition. Here we show that the sequence of the four nucleotides at a hairpin tip modulates transposition efficiency over a surprisingly wide (>100-fold) range. Some hairpin targets stimulate extraordinarily efficient transposition (up to 15%); one serves as a potent and specific transposition inhibitor, blocking capture of targets and destabilizing preformed target capture complexes. These findings suggest novel regulatory possibilities and may provide insight into the activities of other transposases.

(CAG)*(CTG) repeats associated with neurodegenerative diseases are stable in the Escherichia coli chromosome.

(CAG)(n)*(CTG)(n) expansion is associated with many neurodegenerative diseases. Repeat instability has been extensively studied in bacterial plasmids, where repeats undergo deletion at high rates. We report an assay for (CAG)(n)*(CTG)(n) deletion from the chloramphenicol acetyltransferase gene integrated into the Escherichia coli chromosome. In strain AB1157, deletion rates for 25-60 (CAG) x (CTG) repeats integrated in the chromosome ranged from 6.88 x 10(-9) to 1.33 x 10(-10), or approximately 6,300 to 660,000-fold lower than in plasmid pBR325. In contrast to the situation in plasmids, deletions occur at a higher rate when (CTG)(43), rather than (CAG)(43), comprised the leading template strand, and complete rather than partial deletions were the predominant mutation observed. Repeats were also stable on long term growth following multiple passages through exponential and stationary phase. Mutations in priA and recG increased or decreased deletion rates, but repeats were still greatly stabilized in the chromosome. The remarkable stability of (CAG)(n) x (CTG)(n) repeats in the E. coli chromosome may result from the differences in the mechanisms for replication or the probability for recombination afforded by a high plasmid copy number. The integration of (CAG)(n) x (CTG)(n) repeats into the chromosome provides a model system in which the inherent stability of these repeats reflects that in the human genome more closely.

Replication restart: a pathway for (CTG).(CAG) repeat deletion in Escherichia coli.

(CTG)n.(CAG)n repeats undergo deletion at a high rate in plasmids in Escherichia coli in a process that involves RecA and RecB. In addition, DNA replication fork progression can be blocked during synthesis of (CTG)n.(CAG)n repeats. Replication forks stalled at (CTG)n.(CAG)n repeats may be rescued by replication restart that involves recombination as well as enzymes involved in replication and DNA repair, and this process may be responsible for the high rate of repeat deletion in E. coli. To test this hypothesis (CAG)n.(CTG)n deletion rates were measured in several E. coli strains carrying mutations involved in replication restart. (CAG)n.(CTG)n deletion rates were decreased, relative to the rates in wild type cells, in strains containing mutations in priA, recG, ruvAB, and recO. Mutations in priB and priC resulted in small reductions in deletion rates. In a recF strain, rates were decreased when (CAG)n comprised the leading template strand, but rates were increased when (CTG)n comprised the leading template. Deletion rates were increased slightly in a recJ strain. The mutational spectra for most mutant strains were altered relative to those in parental strains. In addition, purified PriA and RecG proteins showed unexpected binding to single-stranded, duplex, and forked DNAs containing (CAG)n and/or (CTG)n loop-outs in various positions. The results presented are consistent with an interpretation that the high rates of trinucleotide repeat instability observed in E. coli result from the attempted restart of replication forks stalled at (CAG)n.(CTG)n repeats.

DNA strand arrangement within the SfiI-DNA complex: atomic force microscopy analysis.

The SfiI restriction enzyme binds to DNA as a tetramer holding two usually distant DNA recognition sites together before cleavage of the four DNA strands. To elucidate structural properties of the SfiI-DNA complex, atomic force microscopy (AFM) imaging of the complexes under noncleaving conditions (Ca2+ instead of Mg2+ in the reaction buffer) was performed. Intramolecular complexes formed by protein interaction between two binding sites in one DNA molecule (cis interaction) as well as complexes formed by the interaction of two sites in different molecules (trans interaction) were analyzed. Complexes were identified unambiguously by the presence of a tall spherical blob at the DNA intersections. To characterize the path of DNA within the complex, the angles between the DNA helices in the proximity of the complex were systematically analyzed. All the data show clear-cut bimodal distributions centered around peak values corresponding to 60 degrees and 120 degrees. To unambiguously distinguish between the crossed and bent models for the DNA orientation within the complex, DNA molecules with different arm lengths flanking the SfiI binding site were designed. The analysis of the AFM images for complexes of this type led to the conclusion that the DNA recognition sites within the complex are crossed. The angles of 60 degrees or 120 degrees between the DNA helices correspond to a complex in which one of the helices is flipped with respect to the orientation of the other. Complexes formed by five different recognition sequences (5'-GGCCNNNNNGGCC-3'), with different central base pairs, were also analyzed. Our results showed that complexes containing the two possible orientations of the helices were formed almost equally. This suggests no preferential orientation of the DNA cognate site within the complex, suggesting that the central part of the DNA binding site does not form strong sequence specific contacts with the protein.

Duplications between direct repeats stabilized by DNA secondary structure occur preferentially in the leading strand during DNA replication.

To ascertain a leading or lagging strand preference for duplication mutations, several short DNA sequences, i.e. mutation inserts, were designed that should demonstrate an asymmetric propensity for duplication mutations in the two complementary DNA strands during replication. The design of the mutation insert involved a 7-bp quasi inverted repeat that forms a remarkably stable hairpin in one DNA strand, but not the other. The inverted repeat is asymmetrically placed between flanking direct repeats. This sequence was cloned into a modified chloramphenicol acetyltransferase (CAT) gene containing a -1 frameshift mutation. Duplication of the mutation insert restores the reading frame of the CAT gene resulting in a chloramphenicol resistant phenotype. The mutation insert showed greater than a 200-fold preference for duplication mutations during leading strand, compared with lagging strand, replication. This result suggests that misalignment stabilized by DNA secondary structure, leading to duplication between direct repeats, occurred preferentially during leading strand synthesis.