Sequence Specific Interaction between Promoter DNA and Escherichia coli RNA Polymerase: Comparative Thermodynamic Analysis with One Immobilized Partner

Sequence specific interaction between DNA and protein molecules has been a subject of active investigation for decades now. Here, we have chosen single promoter containing bacteriophage Delta D-III T7 DNA and Escherichia coli RNA polymerase and followed their recognition at the air-water interface by using the surface plasmon resonance (SPR) technique, where the movement of one of the reacting species is restricted by way of arraying them on an immobilized support. For the Langmuir monolayer studies, we used a RNA polymerase with a histidine tag attached to one of its subunits, thus making it an xcellent substrate for Ni(II) ions, while the SPR Studies were done using biotin-labeled DNA immobilized on a streptavidin-coated chip. Detailed analysis of the thermodynamic parameters as a function of concentration and temperature revealed that the interaction of RNA polymerase with T7 DNA is largely entropy driven (83 (+/- 12) kcal mol(-1)) with a positive enthalpy of 13.6 (+/- 3.6) kcal mol(-1), The free energy of reaction determined by SPR and Langmuir-Blodgett technique was -11 (+/- 2) and -15.6 kcal mol(-1), respectively. The ability of these methods to retain the specificity of the recognition process was also established.

Compound heterozygous mutations in RIPPLY2 associated with vertebral segmentation defects

Segmentation defects of the vertebrae (SDV) are caused by aberrant somite formation during embryogenesis and result in irregular formation of the vertebrae and ribs. The Notch signal transduction pathway plays a critical role in somite formation and patterning in model vertebrates. In humans, mutations in several genes involved in the Notch pathway are associated with SDV, with both autosomal recessive (MESP2, DLL3, LFNG, HES7) and autosomal dominant (TBX6) inheritance. However, many individuals with SDV do not carry mutations in these genes. Using whole-exome capture and massive parallel sequencing, we identified compound heterozygous mutations in RIPPLY2 in two brothers with multiple regional SDV, with appropriate familial segregation. One novel mutation (c.A238T:p.Arg80*) introduces a premature stop codon. In transiently transfected C2C12 mouse myoblasts, the RIPPLY2 mutant protein demonstrated impaired transcriptional repression activity compared with wild-type RIPPLY2 despite similar levels of expression. The other mutation (c.240-4T>G), with minor allele frequency <0.002, lies in the highly conserved splice site consensus sequence 5' to the terminal exon. Ripply2 has a well-established role in somitogenesis and vertebral column formation, interacting at both gene and protein levels with SDV-associated Mesp2 and Tbx6. We conclude that compound heterozygous mutations in RIPPLY2 are associated with SDV, a new gene for this condition. © The Author 2014.

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

Next generation sequencing identifies novel CACNA1A gene mutations in Episodic Ataxia type 2

Episodic Ataxia type 2 (EA2) is a rare autosomal dominantly inherited neurological disorder characterized by recurrent disabling imbalance, vertigo and episodes of ataxia lasting minutes to hours. EA2 is caused most often by loss of function mutations of the calcium channel gene CACNA1A. In addition to EA2, mutations in CACNA1A are responsible for two other allelic disorders: familial hemiplegic migraine type1 (FHM1) and spinocerebellar ataxia type 6 (SCA6). Herein, we have utilised Next Generation Sequencing (NGS) to screen the coding sequence, exon-intron boundaries and UTRs of five genes where mutation is known to produce symptoms related to EA2, including CACNA1A. We performed this screening in a group of 31 unrelated patients with EA2 symptoms. Both novel and known mutations were detected through NGS technology, and confirmed through Sanger sequencing. Genetic testing showed in total 15 mutation bearing patients (48%), of which 9 were novel mutations (6 missense and 3 small frameshift deletion mutations) and six known mutations (4 missense and 2 nonsense).These results demonstrate the efficiency of our NGS-panel for detecting known and novel mutations for EA2 in the CACNA1A gene, also identifying a novel missense mutation in ATP1A2 which is not a normal target for EA2 screening.

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.