GC skew in protein-coding genes between the leading and lagging strands in bacterial genomes: new substitution models incorporating strand bias

The DNA strands in most prokaryotic genomes experience strand-biased spontaneous mutation, especially C→T mutations produced by deamination that occur preferentially in the leading strand. This has often been invoked to account for the asymmetry in nucleotide composition, typically measured by GC skew, between the leading and the lagging strand. Casting such strand asymmetry in the framework of a nucleotide substitution model is important for understanding genomic evolution and phylogenetic reconstruction. We present a substitution model showing that the increased C→T mutation will lead to positive GC skew in one strand but negative GC skew in the other, with greater C→T mutation pressure associated with greater differences in GC skew between the leading and the lagging strand. However, the model based on mutation bias alone does not predict any positive correlation in GC skew between the leading and lagging strands. We computed GC skew for coding sequences collinear with the leading and lagging strands across 339 prokaryotic genomes and found a strong and positive correlation in GC skew between the two strands. We show that the observed positive correlation can be satisfactorily explained by an improved substitution model with one additional parameter incorporating a general trend of C avoidance.     

Modelling the probability distribution of the number of DNA double-strand breaks due to sporadic alkylation of nucleotide bases

Metabolites and certain chemical agents (for example methyl methanesulfonate) can induce nucleotide bases on chromosomal strands to become alkylated. These alkylated sites have the potential to become single-strand chromosomal breaks, a form of DNA damage, if they are exposed to a sufficient temperature in vitro. It has been proposed that a single-strand break (SSB) sufficiently close to another SSB on the opposite chromosomal strand will form a double-strand break (DSB). DNA repair mechanisms are less able to repair DSBs compared to SSBs. Because of the complex three-dimensional structure of DNA, some chromosomal regions are more susceptible to alkylation than others. A question of interest is therefore whether these alkylated bases are randomly distributed or tend to be clustered. Pulsed-field gel electrophoresis allows the number of DNA fragments (and hence the number of DSBs) to be observed directly. The randomness of alkylation events can therefore be tested using the standard statistical hypothesis-testing framework. Under the null hypothesis, that the SSBs are randomly distributed on each of the strands, we can calculate the probability of observing a number of DSBs at least as large as that observed and hence the associated p-value. Previously, the probability distribution of the number of DSBs has been determined by Monte Carlo simulations; when considering the whole genome this can be very time consuming. In this paper, we theoretically derive an approximation to the distribution enabling appropriate probabilities to be calculated quickly. Based on previous findings we assume that the number of breaks on each strand is small compared to the number of nucleotide bases. We show that our method can give the correct probability distribution when alkylation events are relatively rare, discuss how rare these events have to be and suggest potential extensions to the model when a greater proportion of bases are alkylated.     

A model for cooperative ligand binding at complementary sites of DNA

The process of cooperative binding of ligands to DNA has been classified into different modes. An additional mode of cooperative interaction amongst ligands binding at sites on complementary strands has been emphasised. A statistical mechanical method has been applied to obtain an analytical expression for the fraction of nucleotide sites bound. Theoretical Scatchard plots have been drawn and analysed.     

Coding capacity of complementary DNA strands.

A Fortran computer algorithm has been used to analyze the nucleotide sequence of several structural genes. The analysis performed on both coding and complementary DNA strands shows that whereas open reading frames shorter than 100 codons are randomly distributed on both DNA strands, open reading frames longer than 100 codons ("virtual genes") are significantly more frequent on the complementary DNA strand than on the coding one. These "virtual genes" were further investigated by looking at intron sequences, splicing points, signal sequences and by analyzing gene mutations. On the basis of this analysis coding and complementary DNA strands of several eukaryotic structural genes cannot be distinguished. In particular we suggest that the complementary DNA strand of the human epsilon-globin gene might indeed code for a protein.     

Elucidating a key intermediate in homologous DNA strand exchange: structural characterization of the RecA-triple-stranded DNA complex using fluorescence resonance energy transfer

The RecA protein of Escherichia coli plays essential roles in homologous recombination and restarting stalled DNA replication forks. In vitro, the protein mediates DNA strand exchange between single-stranded (ssDNA) and homologous double-stranded DNA (dsDNA) molecules that serves as a model system for the in vivo processes. To date, no high-resolution structure of the key intermediate, comprised of three DNA strands simultaneously bound to a RecA filament (RecA-tsDNA complex), has been reported. We present a systematic characterization of the helical geometries of the three DNA strands of the RecA-tsDNA complex using fluorescence resonance energy transfer (FRET) under physiologically relevant solution conditions. FRET donor and acceptor dyes were used to label different DNA strands, and the interfluorophore distances were inferred from energy transfer efficiencies measured as a function of the base-pair separation between the two dyes. The energy transfer efficiencies were first measured on a control RecA-dsDNA complex, and the calculated helical parameters (h approximately 5 A, Omega(h) approximately 20 degrees ) were consistent with structural conclusions derived from electron microscopy (EM) and other classic biochemical methods. Measurements of the helical parameters for the RecA-tsDNA complex revealed that all three DNA strands adopt extended and unwound conformations similar to those of RecA-bound dsDNA. The structural data are consistent with the hypothesis that this complex is a late, post-strand-exchange intermediate with the outgoing strand shifted by about three base-pairs with respect to its registry with the incoming and complementary strands. Furthermore, the bases of the incoming and complementary strands are displaced away from the helix axis toward the minor groove of the heteroduplex, and the bases of the outgoing strand lie in the major groove of the heteroduplex. We present a model for the strand exchange intermediate in which homologous contacts preceding strand exchange arise in the minor groove of the substrate dsDNA.     

Rifampicin-resistant initiation of DNA synthesis on the isolated strands of ColE plasmid DNA.

The opposite strands of the ColE1 and ColE3 plasmids were isolated as circular single-stranded DNA molecules. These molecules were compared with M13 and phi X174 viral DNA with respect to their capacity to function as templates for in vitro DNA synthesis by a replication enzyme fraction from Escherichia coli. It was found for both ColE plasmids that the conversion of H as well as L strands to duplex DNA molecules closely resembles phi X174 complementary strand synthesis and occurs by a rifampicin-resistant priming mechanism involving the dnaB, dnaC, and dnaG gene products. Restriction analysis of partially double-stranded intermediates indicates that preferred start sites for DNA synthesis are present on both strands of the ColE1 HaeII-C fragment. Inspection of the nucleotide sequence of this region reveals structural similarities with the origin of phi X174 complementary strand synthesis. We propose that the rifampicin-resistant initiation site (rri) in the ColE1 L strand is required for the priming of discontinuous lagging strand synthesis during vegetative replication and that the rri site in the H strand is involved in the initiation of L strand synthesis during conjugative transfer.     

The crystal structure of the RNA/DNA hybrid r(GAAGAGAAGC). d(GCTTCTCTTC) shows significant differences to that found in solution.

The crystal structure of the RNA/DNA hybrid r(GAAGAGAAGC). d(GCTTCTCTTC) has been solved and refined at 2.5 A resolution. The refinement procedure converged at R = 0.181 for all reflections in the range 20.0-2.5 A. In the crystal, the RNA/DNA hybrid duplex has an A' conformation with all but one of the nucleotide sugar moieties adopting a C3'- endo (N) conformation. Both strands in the double helix adopt a global conformation close to the A-form and the width of the minor groove is typical of that found in the crystal structures of other A-form duplexes. However, differences are observed between the RNA and DNA strands that make up the hybrid at the local level. In the central portion of the duplex, the RNA strand has backbone alpha, beta and gamma torsion angles that alternate between the normal gauche -/ trans / gauche + conformation and an unusual trans / trans / trans conformation. Coupled with this so-called 'alpha/gamma flipping' of the backbone torsion angles, the distance between adjacent phosphorous atoms on the RNA strand systematically varies. Neither of these phenomena are observed on the DNA strand. The structure of the RNA/DNA hybrid presented here differs significantly from that found in solution for this and other sequences. Possible reasons for these differences and their implications for the current model of RNase H activity are discussed.     

Mutation detection using nucleotide analogs that alter electrophoretic mobility.

A simple primer extension assay has been developed to distinguish homologous DNA segments differing by as little as a single nucleotide. DNA strands are synthesized with one of the four natural nucleotides replaced with an analog that affects electrophoretic mobility. DNAs that are the same length but differ in the number of analog molecules per strand exhibit different mobilities on a sequencing gel. In combination with the polymerase chain reaction (PCR; 1, 2), this method has been used to distinguish mutant and normal alleles of the human insulin receptor gene that differ by a single-base substitution. The method appears to be generally applicable to the detection of any nucleotide polymorphism in any segment of DNA.     

Substrate preferences of human placental DNA methyltransferase investigated with synthetic polydeoxynucleotides

A DNA methyltransferase partly purified from human placenta has been tested on a variety of synthetic polydeoxynucleotides. The results showed that: the enzyme is most active as a 'maintenance' or 'hemi-' methylase but also has some de novo methylating activity; the presence or absence of A or T in the substrate strand has little influence on maintenance or de novo activity, while polymers containing C but not G in the same strand are poor de novo substrates and bind poorly to the enzyme; single-stranded polymers are about as good substrates as double-stranded ones, and the effects of nucleotide composition (particularly G and mC content) on enzyme activity with single strands are similar to those with double-stranded polymers; strands in which all the cytosines are methylated bind the enzyme well. A mechanism is suggested involving two different sites on the enzyme that recognize CG and mCG, and which rationalizes the activity of eukaryotic DNA methyltransferases towards single-stranded DNA.     

Fidelity of replication of the leading and the lagging DNA strands opposite N-methyl-N-nitrosourea-induced DNA damage in human cells.

Semi-conservative replication of double-stranded DNA in eukaryotic cells is an asymmetric process involving leading and lagging strand synthesis and different DNA polymerases. We report a study to analyze the effect of these asymmetries when the replication machinery encounters alkylation-induced DNA adducts. The model system is an EBV-derived shuttle vector which replicates in synchrony with the host human cells and carries as marker gene the bacterial gpt gene. A preferential distribution of N-methyl-N-nitrosourea (MNU)-induced mutations in the non transcribed DNA strand of the shuttle vector pF1-EBV was previously reported. The hypermutated strand was the leading strand. To test whether the different fidelity of DNA polymerases synthesizing the leading and the lagging strands might contribute to MNU-induced mutation distribution the mutagenesis study was repeated on the shuttle vector pTF-EBV which contains the gpt gene in the inverted orientation. We show that the base substitution error rates on an alkylated substrate are similar for the replication of the leading and lagging strands. Moreover, we present evidence that the fidelity of replication opposite O6-methylguanine adducts of both the leading and lagging strands is not affected by the 3' flanking base. The preferential targeting of mutations after replication of alkylated DNA is mainly driven by the base at the 5' side of the G residues.     

Deoxyribonuclease I produces staggered cuts in the DNA of chromatin

The relationship of cuts made by deoxyribonuclease I (DNase I, EC. 3.1.4.5) on the two strands of DNA of chromatin has been investigated. DNA was extracted from a DNase I digest of rat liver nuclei and incubated with the large fragment of DNA polyrnerase I. Analysis of the products of this incubation indicates the cuts made by DNase I on opposite strands are staggered with respect to one another. A cut on one strand is about two bases in the 3′ direction or eight bases in the 5′ direction from the position on its own strand which is directly across from the cut on the other strand. A different result is obtained when a DNase I digest of native DNA is analyzed. Current models for the organization of DNA in the nucleosome are discussed with respect to these results.     

The stretched DNA geometry of recombination and repair nucleoprotein filaments

The RecA protein of Escherichia coli plays essential roles in homologous recombination and restarting stalled DNA replication forks. In vitro, the protein mediates DNA strand exchange between single-stranded (ssDNA) and homologous double-stranded DNA (dsDNA) molecules that serves as a model system for the in vivo processes. To date, no high-resolution structure of the key intermediate, comprised of three DNA strands simultaneously bound to a RecA filament (RecA x tsDNA complex), has been elucidated by classical methods. Here we review the systematic characterization of the helical geometries of the three DNA strands of the RecA x tsDNA complex using fluorescence resonance energy transfer (FRET) under physiologically relevant solution conditions. Measurements of the helical parameters for the RecA x tsDNA complex are consistent with the hypothesis that this complex is a late, poststrand-exchange intermediate with the outgoing strand shifted by about three base pairs with respect to its registry with the incoming and complementary strands. All three strands in the RecA x tsDNA complex adopt extended and unwound conformations similar to those of RecA-bound ssDNA and dsDNA.     

Analysis of a circular code model

A circular code has been identified in the protein (coding) genes of both eukaryotes and prokaryotes by using a statistical method called trinucleotide frequency (TF) method [Arquès & Michel (1996). J. theor. Biol. 182, 45-58]. Recently, a probabilistic model based on the nucleotide frequencies with a hypothesis of absence of correlation between successive bases on a DNA strand, has been proposed by Koch & Lehmann [(1997). J. theor. Biol. 189, 171-174] for constructing some particular circular codes. Their interesting method which we call here nucleotide frequency (NF) method, reveals several limits for constructing the circular code observed with protein genes.     

A 189-bp fragment of Crithidia fasciculata maxicircle DNA confers autonomous replication in Saccharomyces cerevisiae

A 189-bp fragment capable of promoting high-frequency transformation in Saccharomyces cerevisiae has been isolated from the maxicircle of the insect trypanosomatid Crithidia fasciculata. Chimeric plasmids containing this autonomously replicating sequence (ars) are maintained as extrachromosomal elements in S. cerevisiae. The nucleotide sequence of the maxicircle fragment, termed ars189, has been determined, and its position has been mapped in the maxicircle. The ars189 fragment has an A + T content of 79.4% and shows a large asymmetry in the distribution of adenine and thymine residues between the two strands. In one strand (the T strand) thymine accounts for 118 out of 189 nucleotides while adenine accounts for only 32 nucleotides. The ars189 DNA does not hybridize with minicircles, and its sequence appears to be unique in the C. fasciculata maxicircle genome. This sequence also shows extensive homology to a sequence within a 2.6-kb ars fragment of the Leishmania tarentolae maxicircle. In addition, ars189 contains two copies of a yeast consensus ars sequence (A/T)TTTATPuTTT(T/A).     

DNA Bearing Bulky Fluorescent and Photoreactive Damage in Both Strands as Substrates of the Nucleotide Excision Repair System

Model DNA molecules that contain bulky lesions in both strands have been created, and their properties as substrates of the nucleotide excision repair (NER) system have been analyzed. The modified nucleoside, 5-[3-(4-azido-2,3,5,6-tetrafluorobenzamido)-1-propoxypropyl]-2′-deoxycytidine (dC^FAB), or the nonnucleoside fragment, N-[6-(9-anthracenylcarbamoyl)hexanoyl]-3-amino-1,2-propanediol (nAnt), have been inserted as damage in certain positions of the first DNA strand (“0”). The position of N-[6-5(6)- fluoresceinylcarbamoyl]hexanoyl]-3-amino-1,2-propanediol (nFlu) has been varied within the second DNA strand. This residue has been located opposite the removable damaging fragment of the first strand at positions–20,–10,–4, 0, +3, and +8 relative to the first lesion). It has been demonstrated that the presence of nFlu at the–4, 0, or +3 position of the second strand significantly reduces the thermostability of DNA duplexes, especially in the case of nAnt-DNA and completely excludes the possibility of NER-catalyzed excision from dC^FAB- and nAnt-containing 137-meric DNA with the second lesion at these positions. The introduction of nFlu at positions–20,–10, or +8 differently affects the excision efficiency of dC^FAB- and nAntcontaining fragments from the first strand. The excision efficiency of dC^FAB-containing fragments from extended double-damaged DNA is as high as from DNA that contains a single dC^FAB damage, while the excision of nAnt-containing fragments occurs with 80–90% lower efficiency from double-damaged DNA occurs from DNA that contains the single nAnt insert. The nFlu insert differently affects the interaction of the sensory XPC-HR23B dimer with dC^FAB- and nAnt-containing DNAs, although in all cases, this interaction occurs with increased efficiency compared to that with single-damaged DNAs. No direct correlation between the thermostability of the DNA duplex and XPC-DNA affinity on the one hand, and the excision efficiency of lesions on the other hand has been shown. The absence of the correlation may be caused by both functional features of variable multiprotein complexes involved in the recognition and verification of damage during NER and the sensitivity of the complexes to the structure of the damage and damage-surrounding DNA. The results are important for understanding the NER mechanism of elimination of bulky damage located in both DNA strands.     

In vitro synthesis of uniform poly(dG)-poly(dC) by Klenow exo- fragment of polymerase I

In this paper, we describe a production procedure of the one-to-one double helical complex of poly(dG)–poly(dC), characterized by a well-defined length (up to 10 kb) and narrow size distribution of molecules. Direct evidence of strands slippage during poly(dG)–poly(dC) synthesis by Klenow exo⁻ fragment of polymerase I is obtained by fluorescence resonance energy transfer (FRET). We show that the polymer extension results in an increase in the separation distance between fluorescent dyes attached to 5′ ends of the strands in time and, as a result, losing communication between the dyes via FRET. Analysis of the products of the early steps of the synthesis by high-performance liquid chromatography and mass spectroscopy suggest that only one nucleotide is added to each of the strand composing poly(dG)–poly(dC) in the elementary step of the polymer extension. We show that proper pairing of a base at the 3′ end of the primer strand with a base in sequence of the template strand is required for initiation of the synthesis. If the 3′ end nucleotide in either poly(dG) or poly(dC) strand is substituted for A, the polymer does not grow. Introduction of the T-nucleotide into the complementary strand to permit pairing with A-nucleotide results in the restoration of the synthesis. The data reported here correspond with a slippage model of replication, which includes the formation of loops on the 3′ ends of both strands composing poly(dG)–poly(dC) and their migration over long-molecular distances (μm) to 5′ ends of the strands.