Problems and Algorithms for Sequence Segmentations

The analysis of sequential data is required in many diverse areas such as telecommunications, stock market analysis, and bioinformatics. A basic problem related to the analysis of sequential data is the sequence segmentation problem. A sequence segmentation is a partition of the sequence into a number of non-overlapping segments that cover all data points, such that each segment is as homogeneous as possible. This problem can be solved optimally using a standard dynamic programming algorithm. In the first part of the thesis, we present a new approximation algorithm for the sequence segmentation problem. This algorithm has smaller running time than the optimal dynamic programming algorithm, while it has bounded approximation ratio. The basic idea is to divide the input sequence into subsequences, solve the problem optimally in each subsequence, and then appropriately combine the solutions to the subproblems into one final solution. In the second part of the thesis, we study alternative segmentation models that are devised to better fit the data. More specifically, we focus on clustered segmentations and segmentations with rearrangements. While in the standard segmentation of a multidimensional sequence all dimensions share the same segment boundaries, in a clustered segmentation the multidimensional sequence is segmented in such a way that dimensions are allowed to form clusters. Each cluster of dimensions is then segmented separately. We formally define the problem of clustered segmentations and we experimentally show that segmenting sequences using this segmentation model, leads to solutions with smaller error for the same model cost. Segmentation with rearrangements is a novel variation to the segmentation problem: in addition to partitioning the sequence we also seek to apply a limited amount of reordering, so that the overall representation error is minimized. We formulate the problem of segmentation with rearrangements and we show that it is an NP-hard problem to solve or even to approximate. We devise effective algorithms for the proposed problem, combining ideas from dynamic programming and outlier detection algorithms in sequences. In the final part of the thesis, we discuss the problem of aggregating results of segmentation algorithms on the same set of data points. In this case, we are interested in producing a partitioning of the data that agrees as much as possible with the input partitions. We show that this problem can be solved optimally in polynomial time using dynamic programming. Furthermore, we show that not all data points are candidates for segment boundaries in the optimal solution.

Evidence of an intermediate phase in ternary Ge7Se93-xSbx glasses

The compositional dependence of thermal properties, such as glass transition temperature (T-g), non-reversing enthalpy change (Delta H-NR) and the specific heat capacity change (Delta C-p) of melt quenched Ge7Se93-xSbx (21 a parts per thousand currency sign x a parts per thousand currency sign 31) glasses, has been studied using alternating differential scanning calorimetry (ADSC) which is analogous to modulated differential scanning calorimetry (MDSC). The glass transition temperature, T-g, which is a measure of global connectivity of the glass, has been found to increase with the addition of Sb. In addition, a change in slope has been observed in the composition dependence of T-g at an average coordination aOE (c) r > = 2.40. The experimentally observed compositional variation of glass transition temperature, has been compared with the theoretical predictions from the stochastic agglomeration theory (SAT) and has been found to be consistent. Further, a narrow thermally reversing window is seen in the compositional variation of the relaxation enthalpy (Delta H-NR), which is centered around aOE (c) r > = 2.40. The change in specific heat capacity (Delta C-p) at T-g is also found to exhibit a distinct minima at aOE (c) r > = 2.40, suggesting that the structural rearrangements for the liquid in the glass transition region are minimized around aOE (c) r > = 2.4.

Characterization of genomic diversity in extraintestinal pathogenic Escherichia coli (ExPEC) and development of a diagnostic DNA microarray for the differentiation of ExPEC isolates causing urinary tract infections

Extraintestinal pathogenic Escherichia coli (ExPEC) represent a diverse group of strains of E. coli, which infect extraintestinal sites, such as the urinary tract, the bloodstream, the meninges, the peritoneal cavity, and the lungs. Urinary tract infections (UTIs) caused by uropathogenic E. coli (UPEC), the major subgroup of ExPEC, are among the most prevalent microbial diseases world wide and a substantial burden for public health care systems. UTIs are responsible for serious morbidity and mortality in the elderly, in young children, and in immune-compromised and hospitalized patients. ExPEC strains are different, both from genetic and clinical perspectives, from commensal E. coli strains belonging to the normal intestinal flora and from intestinal pathogenic E. coli strains causing diarrhea. ExPEC strains are characterized by a broad range of alternate virulence factors, such as adhesins, toxins, and iron accumulation systems. Unlike diarrheagenic E. coli, whose distinctive virulence determinants evoke characteristic diarrheagenic symptoms and signs, ExPEC strains are exceedingly heterogeneous and are known to possess no specific virulence factors or a set of factors, which are obligatory for the infection of a certain extraintestinal site (e. g. the urinary tract). The ExPEC genomes are highly diverse mosaic structures in permanent flux. These strains have obtained a significant amount of DNA (predictably up to 25% of the genomes) through acquisition of foreign DNA from diverse related or non-related donor species by lateral transfer of mobile genetic elements, including pathogenicity islands (PAIs), plasmids, phages, transposons, and insertion elements. The ability of ExPEC strains to cause disease is mainly derived from this horizontally acquired gene pool; the extragenous DNA facilitates rapid adaptation of the pathogen to changing conditions and hence the extent of the spectrum of sites that can be infected. However, neither the amount of unique DNA in different ExPEC strains (or UPEC strains) nor the mechanisms lying behind the observed genomic mobility are known. Due to this extreme heterogeneity of the UPEC and ExPEC populations in general, the routine surveillance of ExPEC is exceedingly difficult. In this project, we presented a novel virulence gene algorithm (VGA) for the estimation of the extraintestinal virulence potential (VP, pathogenicity risk) of clinically relevant ExPECs and fecal E. coli isolates. The VGA was based on a DNA microarray specific for the ExPEC phenotype (ExPEC pathoarray). This array contained 77 DNA probes homologous with known (e.g. adhesion factors, iron accumulation systems, and toxins) and putative (e.g. genes predictably involved in adhesion, iron uptake, or in metabolic functions) ExPEC virulence determinants. In total, 25 of DNA probes homologous with known virulence factors and 36 of DNA probes representing putative extraintestinal virulence determinants were found at significantly higher frequency in virulent ExPEC isolates than in commensal E. coli strains. We showed that the ExPEC pathoarray and the VGA could be readily used for the differentiation of highly virulent ExPECs both from less virulent ExPEC clones and from commensal E. coli strains as well. Implementing the VGA in a group of unknown ExPECs (n=53) and fecal E. coli isolates (n=37), 83% of strains were correctly identified as extraintestinal virulent or commensal E. coli. Conversely, 15% of clinical ExPECs and 19% of fecal E. coli strains failed to raster into their respective pathogenic and non-pathogenic groups. Clinical data and virulence gene profiles of these strains warranted the estimated VPs; UPEC strains with atypically low risk-ratios were largely isolated from patients with certain medical history, including diabetes mellitus or catheterization, or from elderly patients. In addition, fecal E. coli strains with VPs characteristic for ExPEC were shown to represent the diagnostically important fraction of resident strains of the gut flora with a high potential of causing extraintestinal infections. Interestingly, a large fraction of DNA probes associated with the ExPEC phenotype corresponded to novel DNA sequences without any known function in UTIs and thus represented new genetic markers for the extraintestinal virulence. These DNA probes included unknown DNA sequences originating from the genomic subtractions of four clinical ExPEC isolates as well as from five novel cosmid sequences identified in the UPEC strains HE300 and JS299. The characterized cosmid sequences (pJS332, pJS448, pJS666, pJS700, and pJS706) revealed complex modular DNA structures with known and unknown DNA fragments arranged in a puzzle-like manner and integrated into the common E. coli genomic backbone. Furthermore, cosmid pJS332 of the UPEC strain HE300, which carried a chromosomal virulence gene cluster (iroBCDEN) encoding the salmochelin siderophore system, was shown to be part of a transmissible plasmid of Salmonella enterica. Taken together, the results of this project pointed towards the assumptions that first, (i) homologous recombination, even within coding genes, contributes to the observed mosaicism of ExPEC genomes and secondly, (ii) besides en block transfer of large DNA regions (e.g. chromosomal PAIs) also rearrangements of small DNA modules provide a means of genomic plasticity. The data presented in this project supplemented previous whole genome sequencing projects of E. coli and indicated that each E. coli genome displays a unique assemblage of individual mosaic structures, which enable these strains to successfully colonize and infect different anatomical sites.

Plasmids and aromatic degradation in Sphingomonas for bioremediation : Aromatic ring cleavage genes in soil and rhizosphere

Microbial degradation pathways play a key role in the detoxification and the mineralization of polyaromatic hydrocarbons (PAHs), which are widespread pollutants in soil and constituents of petroleum hydrocarbons. In microbiology the aromatic degradation pathways are traditionally studied from single bacterial strains with capacity to degrade certain pollutant. In soil the degradation of aromatics is performed by a diverse community of micro-organisms. The aim of this thesis was to study biodegradation on different levels starting from a versatile aromatic degrader Sphingobium sp. HV3 and its megaplasmid, extending to revelation of diversity of key catabolic enzymes in the environment and finally studying birch rhizoremediation in PAH-polluted soil. To understand biodegradation of aromatics on bacterial species level, the aromatic degradation capacity of Sphingobium sp. HV3 and the role of the plasmid pSKY4, was studied. Toluene, m-xylene, biphenyl, fluorene, phenanthrene were detected as carbon and energy sources of the HV3 strain. Tn5 transposon mutagenesis linked the degradation capacity of toluene, m-xylene, biphenyl and naphthalene to the pSKY4 plasmid and qPCR expression analysis showed that plasmid extradiol dioxygenases genes (bphC and xylE) are inducted by phenanthrene, m-xylene and biphenyl whereas the 2,4-dichlorophenoxyacetic acid herbicide induced the chlorocatechol 1,2-dioxygenase gene (tfdC) from the ortho-pathway. A method to study upper meta-pathway extradiol dioxygenase gene diversity in soil was developed. The extradiol dioxygenases catalyse cleavage of the aromatic ring between a hydroxylated carbon and an adjacent non-hydroxylated carbon (meta-cleavage). A high diversity of extradiol dioxygenases were detected from polluted soils. The detected extradiol dioxygenases showed sequence similarity to known catabolic genes of Alpha-, Beta-, and Gammaproteobacteria. Five groups of extradiol dioxygenases contained sequences with no close homologues in the database, representing novel genes. In rhizoremediation experiment with birch (Betula pendula) treatment specific changes of extradiol dioxygenase communities were shown. PAH pollution changed the bulk soil extradiol dioxygenase community structure and birch rhizosphere contained a more diverse extradiol dioxygenase community than the bulk soil showing a rhizosphere effect. The degradation of pyrene in soil was enhanced with birch seedlings compared to soil without birch. The complete 280,923 kb nucleotide sequence of pSKY4 plasmid was determined. The open reading frames of pSKY4 were divided into putative conjugative transfer, aromatic degradation, replication/maintaining and transposition/integration function-encoding proteins. Aromatic degradation orfs shared high similarity to corresponding genes in pNL1, a plasmid from the deep subsurface strain Novosphingobium aromaticivorans F199. The plasmid backbones were considerably more divergent with lower similarity, which suggests that the aromatic pathway has functioned as a plasmid independent mobile genetic element. The functional diversity of microbial communities in soil is still largely unknown. Several novel clusters of extradiol dioxygenases representing catabolic bacteria, whose function, biodegradation pathways and phylogenetic position is not known were amplified with single primer pair from polluted soils. These extradiol dioxygenase communities were shown to change upon PAH pollution, which indicates that their hosts function in PAH biodegradation in soil. Although the degradation pathways of specific bacterial species are substantially better depicted than pathways in situ, the evolution of degradation pathways for the xenobiotic compounds is largely unknown. The pSKY4 plasmid contains aromatic degradation genes in putative mobile genetic element causing flexibility/instability to the pathway. The localisation of the aromatic biodegradation pathway in mobile genetic elements suggests that gene transfer and rearrangements are a competetive advantage for Sphingomonas bacteria in the environment.

Studies in the Sandalwood Oil Series. III. Chain Effect on Terpene Transformations

The carboxyl chain of some molecules has been found to be responsible for causing rearrangements and controlling their course. This chain effect, which operates during reactions involving carbonium ions, is illustrated with examples from Sandalwood oil chemistry.

Cytosines, but Not Purines, Determine Recombination Activating Gene (RAG)-induced Breaks on Heteroduplex DNA Structures

The sequence specificity of the recombination activating gene (RAG) complex during V(D)J recombination has been well studied. RAGs can also act as structure-specific nuclease; however, little is known about the mechanism of its action. Here, we show that in addition to DNA structure, sequence dictates the pattern and efficiency of RAG cleavage on altered DNA structures. Cytosine nucleotides are preferentially nicked by RAGs when present at single-stranded regions of heteroduplex DNA. Although unpaired thymine nucleotides are also nicked, the efficiency is many fold weaker. Induction of single- or double-strand breaks by RAGs depends on the position of cytosines and whether it is present on one or both of the strands. Interestingly, RAGs are unable to induce breaks when adenine or guanine nucleotides are present at single-strand regions. The nucleotide present immediately next to the bubble sequence could also affect RAG cleavage. Hence, we propose “C(d)C(S)C(S)” (d, double-stranded; s, single-stranded) as a consensus sequence for RAG-induced breaks at single-/double-strand DNA transitions. Such a consensus sequence motif is useful for explaining RAG cleavage on other types of DNA structures described in the literature. Therefore, the mechanism of RAG cleavage described here could explain facets of chromosomal rearrangements specific to lymphoid tissues leading to genomic instability.

Anti-apoptotic Protein BCL2 Down-regulates DNA End Joining in Cancer Cells

Cancer cells are often associated with secondary chromosomal rearrangements, such as deletions, inversions, and translocations, which could be the consequence of unrepaired/misrepaired DNA double strand breaks (DSBs). Nonhomologous DNA end joining is one of the most common pathways to repair DSBs in higher eukaryotes. By using oligomeric DNA substrates mimicking various endogenous DSBs in a cell-free system, we studied end joining (EJ) in different cancer cell lines. We found that the efficiency of EJ varies among cancer cells; however, there was no remarkable difference in the mechanism and expression of EJ proteins. Interestingly, cancer cells with lower levels of EJ possessed elevated expression of BCL2 and vice versa. Removal of BCL2 by immunoprecipitation or protein fractionation led to elevated EJ. More importantly, we show that overexpression of BCL2 or the addition of purified BCL2 led to the down-regulation of EJ. Further, we found that BCL2 interacts with KU proteins both in vitro and in vivo. Hence, our results suggest that EJ in cancer cells could be negatively regulated by the anti-apoptotic protein, BCL2, and this may contribute toward increased chromosomal abnormalities in cancer.

Parallel Tempering Simulations of Liquid-Phase Adsorption of n-Alkane Mixtures in Zeolite LTA-5A

Adsorption of n-alkane mixtures in the zeolite LTA-5A under liquid-phase conditions has been studied using grand canonical Monte Carlo (GCMC) simulations combined with parallel tempering. Normal GCMC techniques fail for some of these systems due to the preference of linear molecules to coil within a single cage in the zeolite. The narrow zeolite windows severerly restrict interactions of the molecules, making it difficult to simulate cooperative rearrangements necessary to explore configuration space. Because of these reasons, normal GCMC simulations results show poor reproducibility in some cases. These problems were overcome with parallel tempering techniques. Even with parallel tempering, these are very challenging systems for molecular simulation. Similar problems may arise for other zeolites such as CHA, AFX, ERI, KFI, and RHO having cages connected by narrow windows. The simulations capture the complex selectivity behavior observed in experiments such as selectivity inversion and azeotrope formation.

Effect of SDS on human hair: Study on the molecular structure and morphology

This paper presents a model study to understand the effect of surfactants on the physicochemical properties of human hair. FT-IR ATR spectroscopy has been employed to understand the chemical changes induced by sodium dodecyl sulfate (SDS) on human scalp hair. In particular, the SDS induced changes in the secondary structure of protein present in the outer protective layer of hair, i.e. cuticle, have been investigated. Conformational changes in the secondary structure of protein were studied by curve fitting of the amide I band after every phase of SDS treatment. It has been found that SDS brings rearrangements in the protein backbone conformations by transforming beta-sheet structure to random coil and beta-turn. Additionally, AFM and SEM studies were carried out to understand the morphological changes induced on the hair surface. SEM and AFM images demonstrated the rupture and partial erosion of cuticle sublayers.

Studies on adsorption of guar gum onto biotite mica

The interaction of guar gum with biotite mica has been investigated through adsorption, flotation and electrokinetic measurements. The adsorption densities of guar gum increase with increase of pH and the isotherms exhibit Langmuirian behaviour. Pretreatment of mica with a complexing agent such as EDTA results in a decrease in the adsorption density, highlighting the contribution of metal ions to the adsorption process. An increase in the surface face-to-edge ratio lends to an increase in the adsorption density. The flotation recoveries decrease as a function of pH, complementing the adsorption results. However, polymer depressant ability is reduced in the case of EDTA treated mica, consequent to reduction of metallic sites. Electrokinetic measurements portray conformational rearrangements of macromolecules with the loading, resulting in the shift of the shear plane, further away from the interface. Dissolution experiments indicate release of metal ions from mica, while co-precipitation tests confirm polymer-metal ion interaction in the bulk solution. The adsorption process is governed by hydrogen bonding as well as chemical interaction between guar gum and the surface metal hydroxide groups of mica. (C) 1997 Published by Elsevier Science Ltd.

Probing the Allosteric Mechanism in Pyrrolysyl-tRNA Synthetase Using Energy-Weighted Network Formalism

Pyrrolysyl-tRNA synthetase (PyIRS) is an atypical enzyme responsible for charging tRNA(Pyl) with pyrrolysine, despite lacking precise tRNA anticodon recognition. This dimeric protein exhibits allosteric regulation of function, like any other tRNA synthetases. In this study we examine the paths of allosteric communication at the atomic level, through energy-weighted networks of Desulfitobacterium hafniense PyIRS (DhPyIRS) and its complexes with tRNA(Pyl) and activated pyrrolysine. We performed molecular dynamics simulations of the structures of these complexes to obtain an ensemble conformation-population perspective. Weighted graph parameters relevant to identifying key players and ties in the context of social networks such as edge/node betweenness, closeness index, and the concept of funneling are explored in identifying key residues and interactions leading to shortest paths of communication in the structure networks of DhPylRS. Further, the changes in the status of important residues and connections and the costs of communication due to ligand induced perturbations are evaluated. The optimal, suboptimal, and preexisting paths are also investigated. Many of these parameters have exhibited an enhanced asymmetry between the two subunits of the dimeric protein, especially in the pretransfer complex, leading us to conclude that encoding of function goes beyond the sequence/structure of proteins. The local and global perturbations mediated by appropriate ligands and their influence on the equilibrium ensemble of conformations also have a significant role to play in the functioning of proteins. Taking a comprehensive view of these observations, we propose that the origin of many functional aspects (allostery rand half-sites reactivity in the case of DhPyIRS) lies in subtle rearrangements of interactions and dynamics at a global level.

Effect of aging on swelling and swell-shrink behavior of a compacted expansive soil

Effect of aging on swelling and swell-shrink behavior of a compacted expansive soil is investigated in this paper. An expansive soil having a liquid limit of 100% is used for this purpose. Compacted specimens were prepared and aged for a predetermined number of days (7, 15, 30, and 90 days) to study their swelling and swell-shrink behavior. It has been shown that aging improves the resistance to compression of compacted specimens. The swelling potentials of specimens also decreased with aging. The dominant factors that influence the aging effects are the water content and degree of saturation at the beginning of the aging process. The changed behavior of aged specimens is attributed to particle rearrangements and formation of bonds, which affect the surface area absorbing water during swelling. The cyclic swell-shrink tests on aged specimens indicated that the differences in vertical displacement during the first swelling were eliminated in the subsequent cycles when specimens were shrunk more, but the aging effect was found to persist with cycles for specimens subjected to lower shrinkage magnitudes.

How does DNA break during chromosomal translocations?

Chromosomal translocations are one of the most common types of genetic rearrangements and are molecular signatures for many types of cancers. They are considered as primary causes for cancers, especially lymphoma and leukemia. Although many translocations have been reported in the last four decades, the mechanism by which chromosomes break during a translocation remains largely unknown. In this review, we summarize recent advances made in understanding the molecular mechanism of chromosomal translocations.

Mechanism of Fragility at BCL2 Gene Minor Breakpoint Cluster Region during t(14;18) Chromosomal Translocation

The t(14;18) translocation in follicular lymphoma is one of the most common chromosomal translocations. Breaks in chromosome 18 are localized at the 3'-UTR of BCL2 gene or downstream and are mainly clustered in either the major breakpoint region or the minor breakpoint cluster region (mcr). The recombination activating gene (RAG) complex induces breaks at IgH locus of chromosome 14, whereas the mechanism of fragility at BCL2 mcr remains unclear. Here, for the first time, we show that RAGs can nick mcr; however, the mechanism is unique. Three independent nicks of equal efficiency are generated, when both Mg2+ and Mn2+ are present, unlike a single nick during V(D)J recombination. Further, we demonstrate that RAG binding and nicking at the mcr are independent of nonamer, whereas a CCACCTCT motif plays a critical role in its fragility, as shown by sequential mutagenesis. More importantly, we recapitulate the BCL2 mcr translocation and find that mcr can undergo synapsis with a standard recombination signal sequence within the cells, in a RAG-dependent manner. Further, mutation to the CCACCTCT motif abolishes recombination within the cells, indicating its vital role. Hence, our data suggest a novel, physiologically relevant, nonamer-independent mechanism of RAG nicking at mcr, which may be important for generation of chromosomal translocations in humans.

Potential G-quadruplex formation at breakpoint regions of chromosomal translocations in cancer may explain their fragility

Genetic alterations like point mutations, insertions, deletions, inversions and translocations are frequently found in cancers. Chromosomal translocations are one of the most common genomic aberrations associated with nearly all types of cancers especially leukemia and lymphoma. Recent studies have shown the role of non-B DNA structures in generation of translocations. In the present study, using various bioinformatic tools, we show the propensity of formation of different types of altered DNA structures near translocation breakpoint regions. In particular, we find close association between occurrence of G-quadruplex forming motifs and fragile regions in almost 70% of genes involved in rearrangements in lymphoid cancers. However, such an analysis did not provide any evidence for the occurrence of G-quadruplexes at the close vicinity of translocation breakpoint regions in nonlymphoid cancers. Overall, this study will help in the identification of novel non-B DNA targets that may be responsible for generation of chromosomal translocations in cancer. (C) 2012 Elsevier Inc. All rights reserved.