PubMed Health⌕ Search

Biomedical subjects

Hamid Bolouri

Publications and source records attributed to Hamid Bolouri.

At least 19 recordsLinked to original sources

Dual feedback loops in the GAL regulon suppress cellular heterogeneity in yeast.

Transcriptional noise is known to be an important cause of cellular heterogeneity and phenotypic variation. The extent to which molecular interaction networks may have evolved to either filter or exploit transcriptional noise is a much debated question. The yeast genetic network regulating galactose metabolism involves two proteins, Gal3p and Gal80p, that feed back positively and negatively, respectively, on GAL gene expression. Using kinetic modeling and experimental validation, we demonstrate that these feedback interactions together are important for (i) controlling the cell-to-cell variability of GAL gene expression and (ii) ensuring that cells rapidly switch to an induced state for galactose uptake.

Computer Simulation↗

Systems biology approaches identify ATF3 as a negative regulator of Toll-like receptor 4.

The innate immune system is absolutely required for host defence, but, uncontrolled, it leads to inflammatory disease. This control is mediated, in part, by cytokines that are secreted by macrophages. Immune regulation is extraordinarily complex, and can be best investigated with systems approaches (that is, using computational tools to predict regulatory networks arising from global, high-throughput data sets). Here we use cluster analysis of a comprehensive set of transcriptomic data derived from Toll-like receptor (TLR)-activated macrophages to identify a prominent group of genes that appear to be regulated by activating transcription factor 3 (ATF3), a member of the CREB/ATF family of transcription factors. Network analysis predicted that ATF3 is part of a transcriptional complex that also contains members of the nuclear factor (NF)-kappaB family of transcription factors. Promoter analysis of the putative ATF3-regulated gene cluster demonstrated an over-representation of closely apposed ATF3 and NF-kappaB binding sites, which was verified by chromatin immunoprecipitation and hybridization to a DNA microarray. This cluster included important cytokines such as interleukin (IL)-6 and IL-12b. ATF3 and Rel (a component of NF-kappaB) were shown to bind to the regulatory regions of these genes upon macrophage activation. A kinetic model of Il6 and Il12b messenger RNA expression as a function of ATF3 and NF-kappaB promoter binding predicted that ATF3 is a negative regulator of Il6 and Il12b transcription, and this hypothesis was validated using Atf3-null mice. ATF3 seems to inhibit Il6 and Il12b transcription by altering chromatin structure, thereby restricting access to transcription factors. Because ATF3 is itself induced by lipopolysaccharide, it seems to regulate TLR-stimulated inflammatory responses as part of a negative-feedback loop.

Activating Transcription Factor 3↗

Expression and functional profiling reveal distinct gene classes involved in fatty acid metabolism.

Cells respond to fatty acid exposure by metabolic reorganization and proliferation of peroxisomes. Described here is the development and application of a genome-wide screen to identify nonessential yeast genes necessary for efficient metabolism of myristic and oleic acids. Comparison of the resultant fitness data set with an integrated data set of genes transcriptionally responsive to fatty acids revealed very little overlap between the data sets. Furthermore, the fitness data set enriched for genes involved in peroxisome biogenesis and other processes related to cell morphology, whereas the expression data set enriched for genes related to metabolism. These data suggest that in response to fatty acid exposure, transcriptional control is biased towards metabolic reorganization, and structural changes tend to be controlled post-transcriptionally. They also suggest that fatty acid responsive metabolic networks are more robust than those related to cell structure. Statistical analyses of these and other global data sets suggest that the utilization of distinct control mechanisms for the execution of morphological versus metabolic responses is widespread.

Fatty Acids↗

A data integration methodology for systems biology.

Different experimental technologies measure different aspects of a system and to differing depth and breadth. High-throughput assays have inherently high false-positive and false-negative rates. Moreover, each technology includes systematic biases of a different nature. These differences make network reconstruction from multiple data sets difficult and error-prone. Additionally, because of the rapid rate of progress in biotechnology, there is usually no curated exemplar data set from which one might estimate data integration parameters. To address these concerns, we have developed data integration methods that can handle multiple data sets differing in statistical power, type, size, and network coverage without requiring a curated training data set. Our methodology is general in purpose and may be applied to integrate data from any existing and future technologies. Here we outline our methods and then demonstrate their performance by applying them to simulated data sets. The results show that these methods select true-positive data elements much more accurately than classical approaches. In an accompanying companion paper, we demonstrate the applicability of our approach to biological data. We have integrated our methodology into a free open source software package named POINTILLIST.

Informatics↗

A data integration methodology for systems biology: experimental verification.

The integration of data from multiple global assays is essential to understanding dynamic spatiotemporal interactions within cells. In a companion paper, we reported a data integration methodology, designated Pointillist, that can handle multiple data types from technologies with different noise characteristics. Here we demonstrate its application to the integration of 18 data sets relating to galactose utilization in yeast. These data include global changes in mRNA and protein abundance, genome-wide protein-DNA interaction data, database information, and computational predictions of protein-DNA and protein-protein interactions. We divided the integration task to determine three network components: key system elements (genes and proteins), protein-protein interactions, and protein-DNA interactions. Results indicate that the reconstructed network efficiently focuses on and recapitulates the known biology of galactose utilization. It also provided new insights, some of which were verified experimentally. The methodology described here, addresses a critical need across all domains of molecular and cell biology, to effectively integrate large and disparate data sets.

Chromatin Immunoprecipitation↗

Computational representation of developmental genetic regulatory networks.

Developmental genetic regulatory networks (GRNs) have unique architectural characteristics. They are typically large-scale, multi-layered, and organized in a nested, modular hierarchy of regulatory network kernels, function-specific building blocks, and structural gene batteries. They are also inherently multicellular and involve changing topological relationships among a growing number of cells. Reconstruction of developmental GRNs requires unique computational tools that support the above representational requirements. In addition, we argue that DNA-centered network modeling, separate descriptions of network organization and network behavior, and support for network documentation and annotation are essential requirements for computational modeling of developmental GRNs. Based on these observations, we have developed a freely available, platform-independent, open source software package (BioTapestry) which supports both the process of model construction and also model visualization, analysis, documentation, and dissemination. We provide an overview of the main features of BioTapestry. The BioTapestry software and additional documents are available from http://www.biotapestry.org. We recommend BioTapestry as the substrate for further co-development for and by the developmental biology community.

Computer Simulation↗

Is the regulation of galactose 1-phosphate tuned against gene expression noise?

The average number of mRNA molecules per active gene in yeast can be remarkably low. Consequently, the relative number of copies of each transcript per cell can vary greatly from moment to moment. When these transcripts are encoding metabolic enzymes, how do the resulting variations in enzyme concentrations affect the regulation of metabolic intermediates? Using a kinetic model of galactose utilization in yeast, we analysed the transmission of noise from transcription and translation on metabolic intermediate regulation. In particular, the effect of the kinetic properties of the galactose-1-phosphate uridylyltransferase reaction on the transmission of noise was analysed.

Galactose↗

Dissecting innate immune responses with the tools of systems biology.

Systems biology strives to derive accurate predictive descriptions of complex systems such as innate immunity. The innate immune system is essential for host defense, yet the resulting inflammatory response must be tightly regulated. Current understanding indicates that this system is controlled by complex regulatory networks, which maintain homoeostasis while accurately distinguishing pathogenic infections from harmless exposures. Recent studies have used high throughput technologies and computational techniques that presage predictive models and will be the foundation of a systems level understanding of innate immunity.

Animals↗

Dizzy: stochastic simulation of large-scale genetic regulatory networks.

We describe Dizzy, a software tool for stochastically and deterministically modeling the spatially homogeneous kinetics of integrated large-scale genetic, metabolic, and signaling networks. Notable features include a modular simulation framework, reusable modeling elements, complex kinetic rate laws, multi-step reaction processes, steady-state noise estimation, and spatial compartmentalization.

Algorithms↗

Control of internal and external noise in genetic regulatory networks.

Positive and negative feedback loops, for example, where a protein regulates its own transcription, play an important role in many genetic regulatory networks. Such systems will be subject to internal noise, which occurs due to the small number of molecules taking part in some reactions. This paper examines the effect of feedback loops on noise levels. Error growth techniques from nonlinear dynamics are used to estimate the variance of a system around a steady-state attractor. It is shown that variablity due to intrinsic stochasticity is directly linked to the stability of the steady state, and therefore to the system's resistance to external perturbations. The methods are demonstrated for a number of simple systems, including a genetic switch with homo-dimerizing regulatory protein, and an oscillator.

Animals↗

A method for estimating stochastic noise in large genetic regulatory networks.

MOTIVATION: Genetic regulatory networks are often affected by stochastic noise, due to the low number of molecules taking part in certain reactions. The networks can be simulated using stochastic techniques that model each reaction as a stochastic event. As models become increasingly large and sophisticated, however, the solution time can become excessive; particularly if one wishes to determine the effect on noise of changes to a series of parameters, or the model structure. Methods are therefore required to rapidly estimate stochastic noise. RESULTS: This paper presents an algorithm, based on error growth techniques from non-linear dynamics, to rapidly estimate the noise characteristics of genetic networks of arbitrary size. The method can also be used to determine analytical solutions for simple sub-systems. It is demonstrated on a number of cases, including a prototype model of the galactose regulatory pathway in yeast. AVAILABILITY: A software tool which incorporates the algorithm is available for use as part of the stochastic simulation package Dizzy. It is available for download at http://labs.systemsbiology.net/bolouri/software/Dizzy/ CONTACT: dorrell@systemsbiology.org SUPPLEMENTARY INFORMATION: A conceptual model of the regulatory part of the galactose utilization pathway in yeast, used as an example in the paper, is available at http://labs.systemsbiology.net/bolouri/models/galconcept.dizzy

Algorithms↗

Transcriptional regulatory cascades in development: initial rates, not steady state, determine network kinetics.

A model was built to examine the kinetics of regulatory cascades such as occur in developmental gene networks. The model relates occupancy of cis-regulatory target sites to transcriptional initiation rate, and thence to RNA and protein output. The model was used to simulate regulatory cascades in which genes encoding transcription factors are successively activated. Using realistic parameter ranges based on extensive earlier measurements in sea urchin embryos, we find that transitions of regulatory states occur sharply in these simulations, with respect to time or changing transcription factor concentrations. As is often observed in developing systems, the simulated regulatory cascades display a succession of gene activations separated by delays of some hours. The most important causes of this behavior are cooperativity in the assembly of cis-regulatory complexes and the high specificity of transcription factors for their target sites. Successive transitions in state occur long in advance of the approach to steady-state levels of the molecules that drive the process. The kinetics of such developmental systems thus depend mainly on the initial output rates of genes activated in response to the advent of new transcription factors.

Animals↗

Next generation simulation tools: the Systems Biology Workbench and BioSPICE integration.

Researchers in quantitative systems biology make use of a large number of different software packages for modelling, analysis, visualization, and general data manipulation. In this paper, we describe the Systems Biology Workbench (SBW), a software framework that allows heterogeneous application components--written in diverse programming languages and running on different platforms--to communicate and use each others' capabilities via a fast binary encoded-message system. Our goal was to create a simple, high performance, opensource software infrastructure which is easy to implement and understand. SBW enables applications (potentially running on separate, distributed computers) to communicate via a simple network protocol. The interfaces to the system are encapsulated in client-side libraries that we provide for different programming languages. We describe in this paper the SBW architecture, a selection of current modules, including Jarnac, JDesigner, and SBWMeta-tool, and the close integration of SBW into BioSPICE, which enables both frameworks to share tools and compliment and strengthen each others capabilities.

Biochemical Phenomena↗

Modeling DNA sequence-based cis-regulatory gene networks.

Gene network analysis requires computationally based models which represent the functional architecture of regulatory interactions, and which provide directly testable predictions. The type of model that is useful is constrained by the particular features of developmentally active cis-regulatory systems. These systems function by processing diverse regulatory inputs, generating novel regulatory outputs. A computational model which explicitly accommodates this basic concept was developed earlier for the cis-regulatory system of the endo16 gene of the sea urchin. This model represents the genetically mandated logic functions that the system executes, but also shows how time-varying kinetic inputs are processed in different circumstances into particular kinetic outputs. The same basic design features can be utilized to construct models that connect the large number of cis-regulatory elements constituting developmental gene networks. The ultimate aim of the network models discussed here is to represent the regulatory relationships among the genomic control systems of the genes in the network, and to state their functional meaning. The target site sequences of the cis-regulatory elements of these genes constitute the physical basis of the network architecture. Useful models for developmental regulatory networks must represent the genetic logic by which the system operates, but must also be capable of explaining the real time dynamics of cis-regulatory response as kinetic input and output data become available. Most importantly, however, such models must display in a direct and transparent manner fundamental network design features such as intra- and intercellular feedback circuitry; the sources of parallel inputs into each cis-regulatory element; gene battery organization; and use of repressive spatial inputs in specification and boundary formation. Successful network models lead to direct tests of key architectural features by targeted cis-regulatory analysis.

DNA↗

New computational approaches for analysis of cis-regulatory networks.

The investigation and modeling of gene regulatory networks requires computational tools specific to the task. We present several locally developed software tools that have been used in support of our ongoing research into the embryogenesis of the sea urchin. These tools are especially well suited to iterative refinement of models through experimental and computational investigation. They include: BioArray, a macroarray spot processing program; SUGAR, a system to display and correlate large-BAC sequence analyses; SeqComp and FamilyRelations, programs for comparative sequence analysis; and NetBuilder, an environment for creating and analyzing models of gene networks. We also present an overview of the process used to build our model of the Strongylocentrotus purpuratus endomesoderm gene network. Several of the tools discussed in this paper are still in active development and some are available as open source.

Chromosomes, Artificial, Bacterial↗

A provisional regulatory gene network for specification of endomesoderm in the sea urchin embryo.

We present the current form of a provisional DNA sequence-based regulatory gene network that explains in outline how endomesodermal specification in the sea urchin embryo is controlled. The model of the network is in a continuous process of revision and growth as new genes are added and new experimental results become available; see http://www.its.caltech.edu/~mirsky/endomeso.htm (End-mes Gene Network Update) for the latest version. The network contains over 40 genes at present, many newly uncovered in the course of this work, and most encoding DNA-binding transcriptional regulatory factors. The architecture of the network was approached initially by construction of a logic model that integrated the extensive experimental evidence now available on endomesoderm specification. The internal linkages between genes in the network have been determined functionally, by measurement of the effects of regulatory perturbations on the expression of all relevant genes in the network. Five kinds of perturbation have been applied: (1) use of morpholino antisense oligonucleotides targeted to many of the key regulatory genes in the network; (2) transformation of other regulatory factors into dominant repressors by construction of Engrailed repressor domain fusions; (3) ectopic expression of given regulatory factors, from genetic expression constructs and from injected mRNAs; (4) blockade of the beta-catenin/Tcf pathway by introduction of mRNA encoding the intracellular domain of cadherin; and (5) blockade of the Notch signaling pathway by introduction of mRNA encoding the extracellular domain of the Notch receptor. The network model predicts the cis-regulatory inputs that link each gene into the network. Therefore, its architecture is testable by cis-regulatory analysis. Strongylocentrotus purpuratus and Lytechinus variegatus genomic BAC recombinants that include a large number of the genes in the network have been sequenced and annotated. Tests of the cis-regulatory predictions of the model are greatly facilitated by interspecific computational sequence comparison, which affords a rapid identification of likely cis-regulatory elements in advance of experimental analysis. The network specifies genomically encoded regulatory processes between early cleavage and gastrula stages. These control the specification of the micromere lineage and of the initial veg(2) endomesodermal domain; the blastula-stage separation of the central veg(2) mesodermal domain (i.e., the secondary mesenchyme progenitor field) from the peripheral veg(2) endodermal domain; the stabilization of specification state within these domains; and activation of some downstream differentiation genes. Each of the temporal-spatial phases of specification is represented in a subelement of the network model, that treats regulatory events within the relevant embryonic nuclei at particular stages.

Animals↗

Robustness as a measure of plausibility in models of biochemical networks.

Theory, experiment, and observation suggest that biochemical networks which are conserved across species are robust to variations in concentrations and kinetic parameters. Here, we exploit this expectation to propose an approach to model building and selection. We represent a model as a mapping from parameter space to behavior space, and utilize bifurcation analysis to study the robustness of each region of steady-state behavior to parameter variations. The hypothesis that potential errors in models will result in parameter sensitivities is tested by analysis of two models of the biochemical oscillator underlying the Xenopus cell cycle. Our analysis successfully identifies known weaknesses in the older model and suggests areas for further investigation in the more recent, more plausible model. It also correctly highlights why the more recent model is more plausible.

Animals↗