Relative stability of DNA as a generic criterion for promoter prediction: whole genome annotation of microbial genomes with varying nucleotide base composition

The rapid increase in genome sequence information has necessitated the annotation of their functional elements, particularly those occurring in the non-coding regions, in the genomic context. Promoter region is the key regulatory region, which enables the gene to be transcribed or repressed, but it is difficult to determine experimentally. Hence an in silico identification of promoters is crucial in order to guide experimental work and to pin point the key region that controls the transcription initiation of a gene. In this analysis, we demonstrate that while the promoter regions are in general less stable than the flanking regions, their average free energy varies depending on the GC composition of the flanking genomic sequence. We have therefore obtained a set of free energy threshold values, for genomic DNA with varying GC content and used them as generic criteria for predicting promoter regions in several microbial genomes, using an in-house developed tool `PromPredict'. On applying it to predict promoter regions corresponding to the 1144 and 612 experimentally validated TSSs in E. coli (50.8% GC) and B. subtilis (43.5% GC) sensitivity of 99% and 95% and precision values of 58% and 60%, respectively, were achieved. For the limited data set of 81 TSSs available for M. tuberculosis (65.6% GC) a sensitivity of 100% and precision of 49% was obtained.

DNA- and chromatin-condensing properties of rat testes H1a and H1t compared to those of rat liver H1bdec: H1t is a poor condenser of chromatin

Histones H1a and H1t are two major linker histone variants present at the pachytene interval of mammalian spermatogenesis. The DNA- and chromatin-condensing properties of these two variants isolated from rat testes were studied and compared with those from rat liver. For this purpose, the histone H1 subtypes were purified from the respective tissues using bath acid and salt extraction procedures, Circular dichroism studies revealed that acid exposure during isolation affects the alpha-helical structure of both the globular domain (in the presence of 1 M NaCl) and the C-terminal lambda-tail (in the presence of 60% trifluoroethanol). The condensation of rat oligonucleosomal DNA, as measured by circular dichroism spectroscopy, by the salt-extracted histone H1 was at least 10 times more efficient than condensation by the acid-extracted histone H1. A site size of 16-20 base pairs was calculated for the salt-extracted histone H1. Among the different histone H1 subtypes, somatic histone H1bdec had the highest DNA-condensing property, followed by histone H1a and histone H1t. All the salt-extracted histones condensed rat oligonucleosomal DNA more efficiently than linear pBR-322 DNA, Histones H1bdec and H1a condensed histone H1-depleted chromatin, prepared from rat liver nuclei, with relatively equal efficiency. On the other hand, there was no condensation of histone H1-depleted chromatin with the testes specific histone H1t. A comparison of the amino acid sequences of histone H1d (rat) and histone H1t (rat) revealed several interesting differences in the occurrence of DNA-binding motifs at the C-terminus. A striking observation is the presence of a direct repeat of an octapeptide motif K(A)T(S)PKKA(S)K(T)K(A) in histone H1d that is absent in histone H1t.

DNA-condensing and chromatin-condensing properties of rat testes hla and hit compared to those of rat-liver hlbdec - hlt is a poor condenser of chromatin

Histones H1a and H1t are two major linker histone variants present at the pachytene interval of mammalian spermatogenesis. The DNA- and chromatin-condensing properties of these two variants isolated from rat testes were studied and compared with those from rat liver. For this purpose, the histone H1 subtypes were purified from the respective tissues using bath acid and salt extraction procedures, Circular dichroism studies revealed that acid exposure during isolation affects the alpha-helical structure of both the globular domain (in the presence of 1 M NaCl) and the C-terminal lambda-tail (in the presence of 60% trifluoroethanol). The condensation of rat oligonucleosomal DNA, as measured by circular dichroism spectroscopy, by the salt-extracted histone H1 was at least 10 times more efficient than condensation by the acid-extracted histone H1. A site size of 16-20 base pairs was calculated for the salt-extracted histone H1. Among the different histone H1 subtypes, somatic histone H1bdec had the highest DNA-condensing property, followed by histone H1a and histone H1t. All the salt-extracted histones condensed rat oligonucleosomal DNA more efficiently than linear pBR-322 DNA, Histones H1bdec and H1a condensed histone H1-depleted chromatin, prepared from rat liver nuclei, with relatively equal efficiency. On the other hand, there was no condensation of histone H1-depleted chromatin with the testes specific histone H1t. A comparison of the amino acid sequences of histone H1d (rat) and histone H1t (rat) revealed several interesting differences in the occurrence of DNA-binding motifs at the C-terminus. A striking observation is the presence of a direct repeat of an octapeptide motif K(A)T(S)PKKA(S)K(T)K(A) in histone H1d that is absent in histone H1t.

Mutations as molecular tools : The metabolic-rate dependent molecular clock and DNA barcoding of allied species

Mutation and recombination are the fundamental processes leading to genetic variation in natural populations. This variation forms the raw material for evolution through natural selection and drift. Therefore, studying mutation rates may reveal information about evolutionary histories as well as phylogenetic interrelationships of organisms. In this thesis two molecular tools, DNA barcoding and the molecular clock were examined. In the first part, the efficiency of mutations to delineate closely related species was tested and the implications for conservation practices were assessed. The second part investigated the proposition that a constant mutation rate exists within invertebrates, in form of a metabolic-rate dependent molecular clock, which can be applied to accurately date speciation events. DNA barcoding aspires to be an efficient technique to not only distinguish between species but also reveal population-level variation solely relying on mutations found on a short stretch of a single gene. In this thesis barcoding was applied to discriminate between Hylochares populations from Russian Karelia and new Hylochares findings from the greater Helsinki region in Finland. Although barcoding failed to delineate the two reproductively isolated groups, their distinct morphological features and differing life-history traits led to their classification as two closely related, although separate species. The lack of genetic differentiation appears to be due to a recent divergence event not yet reflected in the beetles molecular make-up. Thus, the Russian Hylochares was described as a new species. The Finnish species, previously considered as locally extinct, was recognized as endangered. Even if, due to their identical genetic make-up, the populations had been regarded as conspecific, conservation strategies based on prior knowledge from Russia would not have guaranteed the survival of the Finnish beetle. Therefore, new conservation actions based on detailed studies of the biology and life-history of the Finnish Hylochares were conducted to protect this endemic rarity in Finland. The idea behind the strict molecular clock is that mutation rates are constant over evolutionary time and may thus be used to infer species divergence dates. However, one of the most recent theories argues that a strict clock does not tick per unit of time but that it has a constant substitution rate per unit of mass-specific metabolic energy. Therefore, according to this hypothesis, molecular clocks have to be recalibrated taking body size and temperature into account. This thesis tested the temperature effect on mutation rates in equally sized invertebrates. For the first dataset (family Eucnemidae, Coleoptera) the phylogenetic interrelationships and evolutionary history of the genus Arrhipis had to be inferred before the influence of temperature on substitution rates could be studied. Further, a second, larger invertebrate dataset (family Syrphidae, Diptera) was employed. Several methodological approaches, a number of genes and multiple molecular clock models revealed that there was no consistent relationship between temperature and mutation rate for the taxa under study. Thus, the body size effect, observed in vertebrates but controversial for invertebrates, rather than temperature may be the underlying driving force behind the metabolic-rate dependent molecular clock. Therefore, the metabolic-rate dependent molecular clock does not hold for the here studied invertebrate groups. This thesis emphasizes that molecular techniques relying on mutation rates have to be applied with caution. Whereas they may work satisfactorily under certain conditions for specific taxa, they may fail for others. The molecular clock as well as DNA barcoding should incorporate all the information and data available to obtain comprehensive estimations of the existing biodiversity and its evolutionary history.

Transcriptional activation of the Escherichia coli bgl operon: negative regulation by DNA structural elements near the promoter

The bgl operon of Escherichia coil is transcriptionally inactive in wild-type cells. DNA insertion sequences (IS) constitute a major class of spontaneous mutations that activate the cryptic bgl promoter. In an attempt to study the molecular mechanism of activation mediated by insertion sequences, transcription of the bgl promoter was carried out in vitro. Stimulation of transcription is observed when a plasmid containing an insertionally activated bgl promoter is used as a template in the absence of proteins other than RNA polymerase. Deletions that remove sequences upstream of the bgl promoter, and insertion of a 1.2 kb DNA fragment encoding resistance to kanamycin, activate the promoter. Point mutations within a region of dyad symmetry upstream of the promoter, which has the potential to extrude into a cruciform structure under torsional stress, also lead to activation, Introduction of a sequence with dyad symmetry, upstream of an activated bgl promoter carrying a deletion of upstream sequences, results in a fourfold reduction in transcription, These results suggest that the cryptic nature of the bgl promoter is because of the presence of DNA structural elements near the promoter that negatively affect transcription.

Nonspecifically bound proteins spin while diffusing along DNA

It is known that DNA-binding proteins can slide along the DNA helix while searching for specific binding sites, but their path of motion remains obscure. Do these proteins undergo simple one-dimensional (1D) translational diffusion, or do they rotate to maintain a specific orientation with respect to the DNA helix? We measured 1D diffusion constants as a function of protein size while maintaining the DNA-protein interface. Using bootstrap analysis of single-molecule diffusion data, we compared the results to theoretical predictions for pure translational motion and rotation-coupled sliding along the DNA. The data indicate that DNA-binding proteins undergo rotation-coupled sliding along the DNA helix and can be described by a model of diffusion along the DNA helix on a rugged free-energy landscape. A similar analysis including the 1D diffusion constants of eight proteins of varying size shows that rotation-coupled sliding is a general phenomenon. The average free-energy barrier for sliding along the DNA was 1.1 +/- 0.2 k(B)T. Such small barriers facilitate rapid search for binding sites.

The Mycobacterial MsDps2 Protein Is a Nucleoid-Forming DNA Binding Protein Regulated by Sigma Factors sigma(A) and sigma(B)

The Dps (DNA-binding protein from starved cells) proteins from Mycobacterium smegmatis MsDps1 and MsDps2 are both DNA-binding proteins with some differences. While MsDps1 has two oligomeric states, with one of them responsible for DNA binding, MsDps2 has only one DNA-binding oligomeric state. Both the proteins however, show iron-binding activity. The MsDps1 protein has been shown previously to be induced under conditions of starvation and osmotic stress and is regulated by the extra cellular sigma factors sigma(H) and sigma(F). We show here, that the second Dps homologue in M. smegmatis, namely MsDps2, is purified in a DNA-bound form and exhibits nucleoid-like structures under the atomic force microscope. It appears that the N-terminal sequence of Dps2 plays a role in nucleoid formation. MsDps2, unlike MsDps1, does not show elevated expression in nutritionally starved or stationary phase conditions; rather its promoter is recognized by RNA polymerase containing sigma(A) or sigma(B), under in vitro conditions. We propose that due to the nucleoid-condensing ability, the expression of MsDps2 is tightly regulated inside the cells.

O6-Methylguanine Repair by O6-Alkylguanine-DNA Alkyltransferase

O6-Alkylguanine-DNA alkyltransferase (AGT) repairs O6-methylguanine (O6mG) in DNA that is known to cause Mutation and cancer. On the basis of Calculations performed using density functional theory involving the active site of AGT, a mechanism for catalytic demethylation of O6mG to guanine has been proposed. In this mechanism, roles of six amino acids, i.e., Cys145, His 146, Glu172, Tyr114, Lys165, and Ser159 in catalytic demethylation of O6mG are involved. This mechanism has three steps as follows. At the first step, Cys145 in the Cys145-water-His146-Glu172 tetrad is converted to cysteine thiolate anion while at the second step, abstraction of the Tyr114 proton by the N3 site of O6mG occurs in a barrierless manner. In the third step, abstraction of Lys165 proton by deprotonated Tyr114 and transfer of the methyl group of O6mG to the thiolate group of Cys145 anion Occur simultaneously. As AGT is a major target in cancer therapy, identification of the roles of the different amino acids in demethylation of O6mG is expected to be useful in designing efficient AGT inhibitors.

Mutations of mitochondrial DNA polymerase gamma: an important cause of neurological disorders

Thesis focuses on mutations of POLG1 gene encoding catalytic subunit polγ-α of mitochondrial DNA polymerase gamma holoenzyme (polG) and the association of mutations with different clinical phenotypes. In addition, particular defective mutant variants of the protein were characterized biochemically in vitro. PolG-holoenzyme is the sole DNA polymerase found in mitochondria. It is involved in replication and repair of the mitochondrial genome, mtDNA. Holoenzyme also includes the accessory subunit polγ-β, which is required for the enhanced processivity of polγ-α. Defective polγ-α causes accumulation of secondary mutations on mtDNA, which leads to a defective oxidative phosphorylation system. The clinical consequences of such mutations are variable, affecting nervous system, skeletal muscles, liver and other post-mitotic tissues. The aims of the studies included: 1) Determination of the role of POLG1 mutations in neurological syndromes with features of mitochondrial dysfunction and an unknown molecular cause. 2) Development and set up of diagnostic tests for routine clinical purposes. 3) Biochemical characterization of the functional consequences of the identified polγ-α variants. Studies describe new neurological phenotypes in addition to PEO caused by POLG1 mutations, including parkinsonism, premature amenorrhea, ataxia and Parkinson s disease (PD). POLG1 mutations and polymorphisms are both common and/or potential genetic risk factors at least among the Finnish population. The major findings and applications reported here are: 1) POLG1 mutations cause parkinsonism and premature menopause in PEO families in either a recessive or a dominant manner. 2) A common recessive POLG1 mutations (A467T and W748S) in the homozygous state causes severe adult or juvenile-onset ataxia without muscular symptoms or histological or mtDNA abnormalities in muscles. 3) A common recessive pathogenic change A467T can also cause a mild dominant disease in heterozygote carriers. 4) The A467T variant shows reduced polymerase activity due to defective template binding. 5) Rare polyglutamine tract length variants of POLG1 are significantly enriched in Finnish idiopathic Parkinson s disease patients. 6) Dominant mutations are clearly restricted to the highly conserved polymerase domain motifs, whereas recessive ones are more evenly distributed along the protein. The present results highlight and confirm the new role of mitochondria in parkinsonism/Parkinson s disease and describe a new mitochondrial ataxia. Based on these results, a POLG1 diagnostic routine has been set up in Helsinki University Central Hospital (HUSLAB).

A novel DNA intercalator, butylamino-pyrimido[4',5':4,5]selenolo(2,3-b)quinoline, induces cell cycle arrest and apoptosis in leukemic cells

DNA intercalators are one of the most commonly used chemotherapeutic agents. Novel intercalating compounds of pyrimido[4',5':4,5]selenolo(2,3-b)quinoline series having a butylamino or piperazino group at fourth position (BPSQ and PPSQ, respectively) are studied. Our results showed that BPSQ induced cytotoxicity whereas PPSQ was cytostatic. The cytotoxicity induced by BPSQ was concentration- and time-dependent. Cell cycle analysis and tritiated thymidine assay revealed that BPSQ affects the cell cycle progression by arresting at S phase. The absence of p-histone H3 and reduction in the levels of PCNA in the cells treated with BPSQ further confirmed the cell cycle arrest. Further, annexin V staining, DNA fragmentation, nuclear condensation and changes in the expression levels of BCL2/BAD confirmed the activation of apoptosis. Activation of caspase 8 and lack of cleavage of caspase 9, caspase 3 and PARP suggest the possibility of BPSQ triggering extrinsic pathway for induction of apoptosis, which is discussed. Hence, we have identified a novel compound which would have clinical relevance in cancer chemotherapeutics.

Mutational analysis of active-site residues in the Mycobacterium leprae RecA intein, a LAGLIDADG homing endonuclease: Asp(122) and Asp(193) are crucial to the double-stranded DNA cleavage activity whereas Asp(218) is not

Mycobacterium leprae recA harbors an in-frame insertion sequence that encodes an intein homing endonuclease (PI-MleI). Most inteins (intein endonucleases) possess two conserved LAGLIDADG (DOD) motifs at their ctive center. A common feature of LAGLIDADG-type homing endonucleases is that they recognize and cleave the same or very similar DNA sequences. However, PI-MleI is distinctive from other members of the family of LAGLIDADG-type HEases for its modular structure with functionally separable domains for DNA-binding and cleavage, each with distinct sequence preferences. Sequence alignment analyses of PI-MleI revealed three putative LAGLIDADG motifs; however, there is conflicting bioinformatics data in regard to their identity and specific location within the intein polypeptide. To resolve this conflict and to determine the active-site residues essential for DNA target site recognition and double-stranded DNA cleavage, we performed site-directed mutagenesis of presumptive catalytic residues in the LAGLIDADG motifs. Analysis of target DNA recognition and kinetic parameters of the wild-type PI-MleI and its variants disclosed that the two amino acid residues, Asp(122) (in Block C) and Asp(193) (in functional Block E), are crucial to the double-stranded DNA endonuclease activity, whereas Asp(218) (in pseudo-Block E) is not. However, despite the reduced catalytic activity, the PI-MleI variants, like the wild-type PI-MleI, generated a footprint of the same length around the insertion site. The D122T variant showed significantly reduced catalytic activity, and D122A and D193A mutations although failed to affect their DNA-binding affinities, but abolished the double-stranded DNA cleavage activity. On the other hand, D122C variant showed approximately twofold higher double-stranded DNA cleavage activity, compared with the wild-type PI-MleI. These results provide compelling evidence that Asp(122) and Asp(193) in DOD motif I and II, respectively, are bona fide active-site residues essential for DNA cleavage activity. The implications of these results are discussed in this report.

Anomalous reaction-diffusion as a model of nonexponential DNA escape kinetics

We show that data from recent experiments carried out on the kinetics of DNA escape from alpha-hemolysin nanopores [M. Wiggin, C. Tropini, C. T. Cossa, N. N. Jetha, and A. Marziali, Biophys. J. 95, 5317 (2008)] may be rationalized by a model of chain dynamics based on the anomalous diffusion of a particle moving in a harmonic well in the presence of a delta function sink. The experiments of Wiggin found, among other things, that the occasional occurrence of unusually long escape times in the distribution of chain trapping events led to nonexponential decays in the survival probability, S(t), of the DNA molecules within the nanopore. Wiggin ascribed this nonexponentiality to the existence of a distribution of trapping potentials, which they suggested was theresult of stochastic interactions between the bases of the DNA and the amino acids located on the surface of the nanopore. Based on this idea, they showed that the experimentally determined S(t) could be well fit in both the short and long time regimes by a function of the form (1+t/tau)(-alpha) (the so called Becquerel function). In our model, S(t) is found to be given by a Mittag-Leffler function at short times and by a generalized Mittag-Leffler function at long times. By suitable choice of certain parameter values, these functions are found to fit the experimental S(t) even better than the Becquerel function. Anomalous diffusion of DNA within the trap prior to escape over a barrier of fixed height may therefore provide a second, plausible explanation of the data, and may offer fresh perspectives on similar trapping and escape problems.

Photocytotoxic Oxovanadium(IV) Complexes Showing Light-Induced DNA and Protein Cleavage Activity

Oxovanadium(IV) complexes [VO(L)(B)]Cl-2 (1-3), where L is bis(2-benzimidazolylmethyl)amine and B is 1,10-phenanthroline(phen),dipyrido[3,2-d:2',3'-f]quinoxaline(dpq) or dipyrido[3,2-a:2',3'-c]phenazine (dppz), have been prepared, characterized, and their photo-induced DNA and protein cleavage activity studied. The photocytotoxicity of complex 3 has been studied using adenocarcinoma A549 cells, The phen complex 1, structurally characterized by single-crystal X-ray crystallography, shows the presence of a vanadyl group in six-coordinate VON5 coordination geometry. The ligands L and phen display tridentate and bidentate N-donor chelating binding modes, respectively. The complexes exhibit a d-d band near 740 nm in 15% DMF-Tris-HCl buffer (pH 7.2). The phen and dpq complexes display an irreversible cathodic cyclic voltammetric response near -0.8 V in 20% DMF-Tris-HCl buffer having 0.1 M KCl as supporting electrolyte. The dppz complex 3 exhibits a quasi-reversible voltammogram near -0.6 V (vs SCE) that is assignable to the V(IV)-V(III)couple. The complexes bind to calf thymus DNA giving binding constant values in the range of 6.6 x 10(4)-2.9 x 10(5) M-1. The binding site size, thermal melting and viscosity binding data suggest DNA surface and/or groove binding nature of the complexes. The complexes show poor ``chemical nuclease'' activity in dark in the presence of 3-mercaptopropionic acid or hydrogen peroxide. The dpq and dppz complexes are efficient photocleavers of plasmid DNA in UV-A light of 365 nm via a mechanistic pathway that involves formation of both singlet oxygen and hydroxyl radicals. The complexes show significant photocleavage of DNA in near-IR light (>750 nm) via hydroxyl radical pathway. Among the three complexes, the dppz complex 3 shows significant BSA and lysozyme protein cleavage activity in UV-A light of 365 nm via hydroxyl radical pathway. The dppz complex 3 also exhibits photocytotoxicity in non-small cell lung carcinoma/human lung adenocarcinoma A549 cells giving IC50 value of 17 mu M in visible light(IC50 = 175 mu M in dark).

Characterization of DNA Strand Transfer Promoted by Mycobacterium smegmatis RecA Reveals Functional Diversity with Mycobacterium tuberculosis RecA†

The RecA-like proteins constitute a group of DNA strand transfer proteins ubiquitous in eubacteria, eukarya, and archaea. However, the functional relationship among RecA proteins is poorly understood. For instance, Mycobacterium tuberculosis RecA is synthesized as a large precursor, which undergoes an unusual protein-splicing reaction to generate an active form. Whereas the precursor was inactive, the active form promoted DNA strand transfer less efficiently compared to EcRecA. Furthermore, gene disruption studies have indicated that the frequencies of allele exchange are relatively lower in Mycobacterium tuberculosis compared to Mycobacterium smegmatis. The mechanistic basis and the factors that contribute to differences in allele exchange remain to be understood. Here, we show that the extent of DNA strand transfer promoted by the M. smegmatis RecA in vitro differs significantly from that of M. tuberculosis RecA. Importantly, M. smegmatis RecA by itself was unable to promote strand transfer, but cognate or noncognate SSBs rendered it efficient even when added prior to RecA. In the presence of SSB, MsRecA or MtRecA catalyzed strand transfer between ssDNA and varying lengths of linear duplex DNA with distinctly different pH profiles. The factors that were able to suppress the formation of DNA networks greatly stimulated strand transfer reactions promoted by MsRecA or MtRecA. Although the rate and pH profiles of dATP hydrolysis catalyzed by MtRecA and MsRecA were similar, only MsRecA was able to couple dATP hydrolysis to DNA strand transfer. Together, these results provide insights into the functional diversity in DNA strand transfer promoted by RecA proteins of pathogenic and nonpathogenic species of mycobacteria.

The RecA intein of Mycobacterium tuberculosis promotes cleavage of ectopic DNA sites. Implications for the dispersal of inteins in natural populations

The RecA intein of Mycobacterium tuberculosis, a novel double-stranded DNA endonuclease, requires both Mn(2+) and ATP for efficient cleavage of the inteinless recA allele. In this study, we show that Mg(2+) alone was sufficient to stimulate PI-MtuI to cleave double-stranded DNA at ectopic sites. In the absence of Mg(2+), PI-MtuI formed complexes with topologically different forms of DNA containing ectopic recognition sequences with equal affinity but failed to cleave DNA. We observed that PI-MtuI was able to inflict double-strand breaks robustly within the ectopic recognition sequence to generate either a blunt end or 1-2-nucleotide 3'-hydroxyl overhangs. Mutational analyses of the presumptive metal ion-binding ligands (Asp(122), Asp(222), and Glu(220)) together with immunoprecipitation assays provided compelling evidence to link both the Mg(2+)- and Mn(2+) and ATP-dependent endonuclease activities to PI-MtuI. The kinetic mechanism of PI-MtuI promoted cleavage of ectopic DNA sites proceeded through a sequential mechanism with transient accumulation of nicked circular duplex DNA as an intermediate. Together, these data suggest that PI-MtuI, like group II introns, might mediate ectopic DNA transposition and hence its lateral transfer in natural populations.