Automatic Design of Boolean Networks for Modelling Differentiation Trees

Real living cell is a complex system governed by many process which are not yet fully understood: the process of cell differentiation is one of these. In this thesis work we make use of a cell differentiation model to develop gene regulatory networks (Boolean networks) with desired differentiation dynamics. To accomplish this task we have introduced techniques of automatic design and we have performed experiments using various differentiation trees. The results obtained have shown that the developed algorithms, except the Random algorithm, are able to generate Boolean networks with interesting differentiation dynamics. Moreover, we have presented some possible future applications and developments of the cell differentiation model in robotics and in medical research. Understanding the mechanisms involved in biological cells can gives us the possibility to explain some not yet understood dangerous disease, i.e the cancer. Le cellula è un sistema complesso governato da molti processi ancora non pienamente compresi: il differenziamento cellulare è uno di questi. In questa tesi utilizziamo un modello di differenziamento cellulare per sviluppare reti di regolazione genica (reti Booleane) con dinamiche di differenziamento desiderate. Per svolgere questo compito abbiamo introdotto tecniche di progettazione automatica e abbiamo eseguito esperimenti utilizzando vari alberi di differenziamento. I risultati ottenuti hanno mostrato che gli algoritmi sviluppati, eccetto l'algoritmo Random, sono in grado di poter generare reti Booleane con dinamiche di differenziamento interessanti. Inoltre, abbiamo presentato alcune possibili applicazioni e sviluppi futuri del modello di differenziamento in robotica e nella ricerca medica. Capire i meccanismi alla base del funzionamento cellulare può fornirci la possibilità di spiegare patologie ancora oggi non comprese, come il cancro.

Extended kalman filter for estimation of parameters in nonlinear state-space models of biochemical networks.

It is system dynamics that determines the function of cells, tissues and organisms. To develop mathematical models and estimate their parameters are an essential issue for studying dynamic behaviors of biological systems which include metabolic networks, genetic regulatory networks and signal transduction pathways, under perturbation of external stimuli. In general, biological dynamic systems are partially observed. Therefore, a natural way to model dynamic biological systems is to employ nonlinear state-space equations. Although statistical methods for parameter estimation of linear models in biological dynamic systems have been developed intensively in the recent years, the estimation of both states and parameters of nonlinear dynamic systems remains a challenging task. In this report, we apply extended Kalman Filter (EKF) to the estimation of both states and parameters of nonlinear state-space models. To evaluate the performance of the EKF for parameter estimation, we apply the EKF to a simulation dataset and two real datasets: JAK-STAT signal transduction pathway and Ras/Raf/MEK/ERK signaling transduction pathways datasets. The preliminary results show that EKF can accurately estimate the parameters and predict states in nonlinear state-space equations for modeling dynamic biochemical networks.

Isolation, characterization and transcriptional regulation of {\it Xenopus myod\/}

I have cloned cDNAs corresponding to two distinct genes, Xlmf1 and Xlmf25, which encode skeletal muscle-specific, transcriptional regulatory proteins. These proteins are members of the helix-loop-helix family of DNA binding factors, and are most homologous to MyoD1. These two genes have disparate temporal expression patterns during early embryogenesis; although, both transcripts are present exclusively in skeletal muscle of the adult. Xlmf1 is first detected 7 hours after fertilization, shortly after the midblastula transition. Xlmf25 is detected in maternal stores of mRNA, during early cleavage stages of the embryo and throughout later development. Both Xlmf1 and Xlmf25 transcripts are detected prior to the expression of other, previously characterized, muscle-specific genes. The ability of Xlmf1 and Xlmf25 to convert mouse 10T1/2 fibroblasts to a myogenic phenotype demonstrates their activity as myogenic regulatory factors. Additionally, Xlmf1 and Xlmf25 can directly transactivate a reporter gene linked to the muscle-specific, muscle creatine kinase (MCK) enhancer. The functional properties of Xlmf1 and Xlmf25 proteins were further explored by investigating their interactions with the binding site in the MCK enhancer. Analysis of dissociation rates revealed that Xlmf25-E12 dimers had a two-fold lower avidity for this site than did Xlmf1-E12 dimers. Clones containing genomic sequence of Xlmf1 and Xlmf25 have been isolated. Reporter gene constructs containing a lac-z gene driven by Xlmf1 regulatory sequences were analyzed by embryo injections and transfections into cultured muscle cells. Elements within $-$200 bp of the transcription start site can promote high levels of muscle specific expression. Embryo injections show that 3500 bp of upstream sequence is sufficient to drive somite specific expression. EMSAs and DNAse I footprint analysis has shown the discrete interaction of factors with several cis-elements within 200 bp of the transcription start site. Mutation of several of these elements shows a positive requirement for two CCAAT boxes and two E boxes. It is evident from the work performed with this promoter that Xlmf1 is tightly regulated during muscle cell differentiation. This is not surprising given the fact that its gene product is crucial to the determination of cell fate choices. ^

The ribosome as novel target for small regulatory RNAs

Small non-protein-coding RNA (ncRNA) molecules have been recognized recently as major contributors to regulatory networks in controlling gene expression in a highly efficient manner. While the list of validated ncRNAs that regulate crucial cellular processes grows steadily, not a single ncRNA has been identified that directly interacts and regulates the ribosome during protein biosynthesis (with the notable exceptions of 7SL RNA and tmRNA). All of the recently discovered regulatory ncRNAs that act on translation (e.g. microRNAs, siRNAs or antisense RNAs) target the mRNA rather than the ribosome. This is unexpected, given the central position the ribosome plays during gene expression. Furthermore it is strongly assumed that the primordial translation system in the ‘RNA world’ most likely received direct regulatory input from ncRNA-like cofactors. The fundamental question that we would like to ask is: Does the ‘RNA world still communicate’ with the ribosome? To address this question, we have analyzed the small ncRNA interactomes of ribosomes of prokaryotic (H. volcanii, S. aureus) and unicellular eukaryotic model organisms. Deep-sequencing and subsequent bioinformatic analyses revealed thousands of putative ribosome-associated ncRNAs. For a subset of these ncRNA candidates we have gathered experimental evidence that they are expressed in a stress-dependent manner and indeed directly target the ribosome. In the archaeon H. volcanii a tRNA-derived fragment was identified to target the small ribosomal subunit upon alkaline stress in vitro and in vivo. As a consequence of ribosome binding, this tRNA-fragment reduces protein synthesis by interfering with the peptidyl transferase activity. Our data reveal the ribosome as a novel target for small regulatory ncRNAs in all domains of life. Ribosome-bound ncRNAs are capable of fine tuning translation and might represent a so far largely unexplored class of regulatory sRNAs.

The ribosome as novel target for small regulatory RNAs

Small non-protein-coding RNA (ncRNA) molecules have been recognized recently as major contributors to regulatory networks in controlling gene expression in a highly efficient manner. While the list of validated ncRNAs that regulate crucial cellular processes grows steadily, not a single ncRNA has been identified that directly interacts and regulates the ribosome during protein biosynthesis (with the notable exceptions of 7SL RNA and tmRNA). All of the recently discovered regulatory ncRNAs that act on translation (e.g. microRNAs, siRNAs or antisense RNAs) target the mRNA rather than the ribosome. This is unexpected, given the central position the ribosome plays during gene expression. Furthermore it is strongly assumed that the primordial translation system in the ‘RNA world’ most likely received direct regulatory input from ncRNA-like cofactors. The fundamental question that we would like to ask is: Does the ‘RNA world still communicate’ with the ribosome? To address this question, we have analyzed the small ncRNA interactomes of ribosomes of organisms from all three domains of life. Deep-sequencing and subsequent bioinformatic analyses revealed thousands of putative ribosome-associated ncRNAs.1,2 For a subset of these ncRNA candidates we have gathered experimental evidence that they are expressed in a stress-dependent manner and indeed directly target the ribosome. We show that some of these ribosome-bound small ncRNAs are capable of fine tuning protein synthesis in vitro and in vivo. Our data therefore reveal the ribosome as a novel target for small regulatory ncRNAs in all domains of life and suggest the existence of a so far largely unexplored mechanism of translation regulation.

Ribosome-associated ncRNAs (rancRNAs) An emerging class of translation regulators

Small non-protein-coding RNA (ncRNA) molecules represent major contributors to regulatory networks in controlling gene expression in a highly efficient manner. Most of the recently discovered regulatory ncRNAs acting on translation target the mRNA rather than the ribosome (e.g.: miRNAs, siRNAs, antisense RNAs). To address the question, whether ncRNA regulators exist that are capable of modulating the rate of protein production by directly interacting with the ribosome, we have analyzed the small ncRNA interactomes of ribosomes. Deep-sequencing analyses revealed thousands of putative rancRNAs in various model organisms (1,2). For a subset of these ncRNA candidates we have gathered experimental evidence that they associate with ribosomes in a stress-dependent manner and fine-tune the rate of protein biosynthesis (3,4). Many of the investigated rancRNAs appear to be processing products of larger functional RNAs, such as tRNAs (2,3), mRNAs (3), or snoRNAs (2). Post-transcriptional cleavage of RNA to generate smaller fragments is a widespread mechanism that enlarges the structural and functional complexity of cellular RNomes. Our data disclose the ribosome as target for small regulatory RNAs. rancRNAs are found in all domains of life and represent a prevalent but so far largely unexplored class of regulatory molecules (5). Ongoing work in our lab revealed first insight into rancRNA processing and mechanism of this emerging class of translation regulators.

Biology of tRNA halves in Trypanosoma brucei

Although T. brucei has to challenge tremendous environment changes, e.g. switch from the bloodstream form in mammalian hosts to the mid gut form present in tsetse flies, there is no evidence for differential regulation of RNA Pol II transcription. Instead, constitutive transcription appears to occur. This observation indicates that protein levels have to be regulated by post-transcriptional mechanisms. It has been shown that non-protein coding RNAs (ncRNAs) are crucial in regulatory networks (e.g. chromosome remodelling; RNA polymerase activity; mRNA turnover; etc.), but all of the recently discovered ncRNAs involved in translation regulation target the mRNA rather than the ribosome. This is unexpected, since the ribosome has a central role during gene expression and due to the assumption that the primordial translation system most likely received direct regulatory input from small molecules including ncRNA cofactors. In our lab, it has been discovered that ncRNAs are able to directly bind to the ribosome, therefore influencing the translation rate in Haloferax volcanii and Saccharomyces cerevisiae. In order to extend this idea of ribosome-binding ncRNAs in mammalian parasites, we want to investigate this mechanism in T. brucei. Accordingly, we performed a genomic screen for small ribosome-associated RNAs followed by functional analyses of possible candidates. With the help of this genomic screen, we found tRNAs that are alternated and tRNA halves that are differentially expressed upon nutritional stress.

Biology of tRNA halves in Trypanosoma brucei

Although T. brucei has to challenge tremendous environment changes, e.g. switch from the bloodstream form in mammalian hosts to the mid gut form present in tsetse flies, there is no evidence for differential regulation of RNA Pol II transcription. Instead, constitutive transcription appears to occur. This observation indicates that protein levels have to be regulated by post-transcriptional mechanisms. It has been shown that non-protein coding RNAs (ncRNAs) are crucial in regulatory networks (e.g. chromosome remodelling; RNA polymerase activity; mRNA turnover; etc.), but all of the recently discovered ncRNAs involved in translation regulation target the mRNA rather than the ribosome. This is unexpected, since the ribosome has a central role during gene expression and due to the assumption that the primordial translation system most likely received direct regulatory input from small molecules including ncRNA cofactors. In our lab, it has been discovered that ncRNAs are able to directly bind to the ribosome, therefore influencing the translation rate in Haloferax volcanii and Saccharomyces cerevisiae. In order to extend this idea of ribosome-binding ncRNAs in mammalian parasites, we want to investigate this mechanism in T. brucei. Accordingly, we performed a genomic screen for small ribosome-associated RNAs followed by functional analyses of possible candidates. With the help of this genomic screen, we found tRNAs that are alternated and tRNA halves that are differentially expressed upon nutritional stress.

Systems Biology Approaches to Probe Gene Regulation in Bacteria

Mechanisms that allow pathogens to colonize the host are not the product of isolated genes, but instead emerge from the concerted operation of regulatory networks. Therefore, identifying components and the systemic behavior of networks is necessary to a better understanding of gene regulation and pathogenesis. To this end, I have developed systems biology approaches to study transcriptional and post-transcriptional gene regulation in bacteria, with an emphasis in the human pathogen Mycobacterium tuberculosis (Mtb). First, I developed a network response method to identify parts of the Mtb global transcriptional regulatory network utilized by the pathogen to counteract phagosomal stresses and survive within resting macrophages. As a result, the method unveiled transcriptional regulators and associated regulons utilized by Mtb to establish a successful infection of macrophages throughout the first 14 days of infection. Additionally, this network-based analysis identified the production of Fe-S proteins coupled to lipid metabolism through the alkane hydroxylase complex as a possible strategy employed by Mtb to survive in the host. Second, I developed a network inference method to infer the small non-coding RNA (sRNA) regulatory network in Mtb. The method identifies sRNA-mRNA interactions by integrating a priori knowledge of possible binding sites with structure-driven identification of binding sites. The reconstructed network was useful to predict functional roles for the multitude of sRNAs recently discovered in the pathogen, being that several sRNAs were postulated to be involved in virulence-related processes. Finally, I applied a combined experimental and computational approach to study post-transcriptional repression mediated by small non-coding RNAs in bacteria. Specifically, a probabilistic ranking methodology termed rank-conciliation was developed to infer sRNA-mRNA interactions based on multiple types of data. The method was shown to improve target prediction in Escherichia coli, and therefore is useful to prioritize candidate targets for experimental validation.

REGULATION OF ALTERNATIVE CARBON METABOLISM IN CANDIDA ALBICANS

Candida albicans is the most important fungal pathogen of humans. Transcript profiling studies show that upon phagocytosis by macrophages, C. albicans undergoes a massive metabolic reorganization activating genes involved in alternative carbon metabolism, including the glyoxylate cycle, β-oxidation and gluconeogenesis. Mutations in key enzymes such as ICL1 (glyoxylate cycle) and FOX2 (fatty acid β-oxidation) revealed that alternative carbon metabolic pathways are required for full virulence in C. albicans. These studies indicate C. albicans uses non-preferred carbon sources allowing its adaptation to microenvironments were nutrients are scarce. It has become apparent that the regulatory networks required for regulation of alternative carbon metabolism in C. albicans are considerably different from the Saccharomyces cerevisiae paradigm and appear more analogous to the Aspergillus nidulans systems. Well-characterized transcription factors in S. cerevisiae have no apparent phenotype or are missing in C. albicans. CTF1 was found to be a single functional homolog of the A. nidulans FarA/FarB proteins, which are transcription factors required for fatty acid utilization. Both FOX2 and ICL1 were found to be part of a large CTF1 regulon. To increase our understanding of how CTF1 regulates its target genes, including whether regulation is direct or indirect, the FOX2 and ICL1 promoter regions were analyzed using a combination of bioinformatics and promoter deletion analysis. To begin characterizing the FOX2 and ICL1 promoters, 5’ rapid amplification of cDNA ends (5’RACE) was used to identify two transcriptional initiation sites in FOX2 and one in ICL1. GFP reporter assays show FOX2 and ICL1 are rapidly expressed in the presence of alternative carbon sources. Both FOX2 and ICL1 harbor the CCTCGG sequence known to be bound by the Far proteins, hence rendering the motif as a putative CTF1 DNA binding element. In this study, the CCTCGG sequence was found to be essential for FOX2 regulation. However, this motif does not appear to be equally important for the regulation of ICL1. This study supports the notion that although C. albicans has diverged from the paradigms of model fungi, C. albicans has made specific adaptations to its transcription-based regulatory network that may contribute to its metabolic flexibility.

Caenorhabditis elegans Mediator complexes are required for developmental-specific transcriptional activation

Mediator proteins are required for transcriptional regulation of most genes in yeast. Mammalian Mediator homologs also function as transcriptional coactivators in vitro; however, their physiological role in gene-specific transcription is not yet known. To determine the role of Mediator proteins in the development of complex organisms, we purified putative Mediator complexes from Caenorhabditis elegans and analyzed their phenotypes in vivo. C. elegans Mediator homologs were assembled into two multiprotein complexes. RNA interference assays showed that the CeMed6, CeMed7, and CeMed10/CeNut2 gene products are required for the expression of developmentally regulated genes, but are dispensable for expression of the ubiquitously expressed genes tested in this study. Therefore, the gene-specific function of Mediator as an integrator of transcriptional regulatory signals is evolutionarily conserved and is essential for C. elegans development.

Rational design of a bacterial transcriptional cascade for amplifying gene expression capacity

Cascade regulatory circuits have been described that control numerous cell processes, and may provide models for the design of artificial circuits with novel properties. Here we describe the design of a transcriptional regulatory cascade to amplify the cell response to a given signal. We used the salicylate-responsive activators of Pseudomonas putida NahR of the naphthalene degradation plasmid NAH7 and XylS2, a mutant regulator of the TOL plasmid for catabolism of m-xylene and their respective cognate promoters Psal and Pm. Control of the expression of xylS2 with the nahR/Psal system permitted either their selective activation with specific effectors for each protein or the simultaneous activation of both of them with salicylate. When cells face the common effector of the two regulators, both the increase in XylS2 concentration and the stimulation of its activity act synergistically on the Pm promoter, amplifying the gene expression capacity by at least one order of magnitude with respect to the individual systems. By changing the hierarchy of regulators, we showed that the specific features of the downstream regulator were crucial for the amplification effect. Directed changes in the effector profile of the regulators allowed the extension of the amplifying system to other molecular signals.

Id1 regulation of cellular senescence through transcriptional repression of p16/Ink4a

The Id family of helix–loop–helix (HLH) transcriptional regulatory proteins does not possess a basic DNA-binding domain and functions as a negative regulator of basic HLH transcription factors. Id proteins coordinate cell growth and differentiation pathways within mammalian cells and have been shown to regulate G1-S cell-cycle transitions. Although much recent data has implicated Id1 in playing a critical role in modulating cellular senescence, no direct genetic evidence has been reported to substantiate such work. Here we show that Id1-null primary mouse embryo fibroblasts undergo premature senescence despite normal growth profiles at early passage. These cells possess increased expression of the tumor-suppressor protein p16/Ink4a but not p19/ARF, and have decreased cyclin-dependent kinase (cdk) 2 and cdk4 kinase activity. We also show that Id1 is able to directly inhibit p16/Ink4a but not p19/ARF promoter activity via its HLH domain, and that Id1inhibits transcriptional activation at E-boxes within the p16/Ink4a promoter. Our data provide, to our knowledge, the first genetic evidence for a role for Id1 as an inhibitor of cellular senescence and suggest that Id1 functions to delay cellular senescence through repression of p16/Ink4a. Because epigenetic and genetic abrogation of p16/Ink4a function has been implicated in the evolution of several human malignancies, we propose that transcriptional regulation of p16/Ink4a may also provide a mechanism for the dysregulation of normal cellular growth controls during the evolution of human malignancies.

Transient Activation of Meox1 Is an Early Component of the Gene Regulatory Network Downstream of Hoxa2

Hox genes encode transcription factors that regulate morphogenesis in all animals with bilateral symmetry. Although Hox genes have been extensively studied, their molecular function is not clear in vertebrates, and only a limited number of genes regulated by Hox transcription factors have been identified. Hoxa2 is required for correct development of the second branchial arch, its major domain of expression. We now show that Meox1 is genetically downstream from Hoxa2 and is a direct target. Meox1 expression is downregulated in the second arch of Hoxa2 mouse mutant embryos. In chromatin immunoprecipitation (ChIP), Hoxa2 binds to the Meox1 proximal promoter. Two highly conserved binding sites contained in this sequence are required for Hoxa2-dependent activation of the Meox1 promoter. Remarkably, in the absence of Meox1 and its close homolog Meox2, the second branchial arch develops abnormally and two of the three skeletal elements patterned by Hoxa2 are malformed. Finally, we show that Meox1 can specifically bind the DNA sequences recognized by Hoxa2 on its functional target genes. These results provide new insight into the Hoxa2 regulatory network that controls branchial arch identity.

Bistability and switching in the lysis/lysogeny genetic regulatory network of bacteriophage lambda

Bistability and switching are two important aspects of the genetic regulatory network of phage. Positive and negative feedbacks are key regulatory mechanisms in this network. By the introduction of threshold values, the developmental pathway of A phage is divided into different stages. If the protein level reaches a threshold value, positive or negative feedback will be effective and regulate the process of development. Using this regulatory mechanism, we present a quantitative model to realize bistability and switching of phage based on experimental data. This model gives descriptions of decisive mechanisms for different pathways in induction. A stochastic model is also introduced for describing statistical properties of switching in induction. A stochastic degradation rate is used to represent intrinsic noise in induction for switching the system from the lysogenic pathway to the lysis pathway. The approach in this paper represents an attempt to describe the regulatory mechanism in genetic regulatory network under the influence of intrinsic noise in the framework of continuous models. (C) 2003 Elsevier Ltd. All rights reserved.