Showing papers in "Molecular Systems Biology in 2010"

A side effect resource to capture phenotypic effects of drugs

Abstract: The molecular understanding of phenotypes caused by drugs in humans is essential for elucidating mechanisms of action and for developing personalized medicines. Side effects of drugs (also known as adverse drug reactions) are an important source of human phenotypic information, but so far research on this topic has been hampered by insufficient accessibility of data. Consequently, we have developed a public, computer-readable side effect resource (SIDER) that connects 888 drugs to 1450 side effect terms. It contains information on frequency in patients for one-third of the drug–side effect pairs. For 199 drugs, the side effect frequency of placebo administration could also be extracted. We illustrate the potential of SIDER with a number of analyses. The resource is freely available for academic research at http://sideeffects.embl.de.

Feedback between p21 and reactive oxygen production is necessary for cell senescence

Abstract: Cellular senescence—the permanent arrest of cycling in normally proliferating cells such as fibroblasts—contributes both to age-related loss of mammalian tissue homeostasis and acts as a tumour suppressor mechanism. The pathways leading to establishment of senescence are proving to be more complex than was previously envisaged. Combining in-silico interactome analysis and functional target gene inhibition, stochastic modelling and live cell microscopy, we show here that there exists a dynamic feedback loop that is triggered by a DNA damage response (DDR) and, which after a delay of several days, locks the cell into an actively maintained state of ‘deep' cellular senescence. The essential feature of the loop is that long-term activation of the checkpoint gene CDKN1A (p21) induces mitochondrial dysfunction and production of reactive oxygen species (ROS) through serial signalling through GADD45-MAPK14(p38MAPK)-GRB2-TGFBR2-TGFβ. These ROS in turn replenish short-lived DNA damage foci and maintain an ongoing DDR. We show that this loop is both necessary and sufficient for the stability of growth arrest during the establishment of the senescent phenotype.

Omic data from evolved E. coli are consistent with computed optimal growth from genome‐scale models

Abstract: After hundreds of generations of adaptive evolution at exponential growth, Escherichia coli grows as predicted using flux balance analysis (FBA) on genome-scale metabolic models (GEMs). However, it is not known whether the predicted pathway usage in FBA solutions is consistent with gene and protein expression in the wild-type and evolved strains. Here, we report that >98% of active reactions from FBA optimal growth solutions are supported by transcriptomic and proteomic data. Moreover, when E. coli adapts to growth rate selective pressure, the evolved strains upregulate genes within the optimal growth predictions, and downregulate genes outside of the optimal growth solutions. In addition, bottlenecks from dosage limitations of computationally predicted essential genes are overcome in the evolved strains. We also identify regulatory processes that may contribute to the development of the optimal growth phenotype in the evolved strains, such as the downregulation of known regulons and stringent response suppression. Thus, differential gene and protein expression from wild-type and adaptively evolved strains supports observed growth phenotype changes, and is consistent with GEM-computed optimal growth states.

Sequence signatures and mRNA concentration can explain two-thirds of protein abundance variation in a human cell line

Abstract: Transcription, mRNA decay, translation and protein degradation are essential processes during eukaryotic gene expression, but their relative global contributions to steady-state protein concentrations in multi-cellular eukaryotes are largely unknown. Using measurements of absolute protein and mRNA abundances in cellular lysate from the human Daoy medulloblastoma cell line, we quantitatively evaluate the impact of mRNA concentration and sequence features implicated in translation and protein degradation on protein expression. Sequence features related to translation and protein degradation have an impact similar to that of mRNA abundance, and their combined contribution explains two-thirds of protein abundance variation. mRNA sequence lengths, amino-acid properties, upstream open reading frames and secondary structures in the 5′ untranslated region (UTR) were the strongest individual correlates of protein concentrations. In a combined model, characteristics of the coding region and the 3′UTR explained a larger proportion of protein abundance variation than characteristics of the 5′UTR. The absolute protein and mRNA concentration measurements for >1000 human genes described here represent one of the largest datasets currently available, and reveal both general trends and specific examples of post-transcriptional regulation.

The endocannabinoid system links gut microbiota to adipogenesis

Abstract: Obesity is characterised by altered gut microbiota, low-grade inflammation and increased endocannabinoid (eCB) system tone; however, a clear connection between gut microbiota and eCB signalling has yet to be confirmed. Here, we report that gut microbiota modulate the intestinal eCB system tone, which in turn regulates gut permeability and plasma lipopolysaccharide (LPS) levels. The impact of the increased plasma LPS levels and eCB system tone found in obesity on adipose tissue metabolism (e.g. differentiation and lipogenesis) remains unknown. By interfering with the eCB system using CB1 agonist and antagonist in lean and obese mouse models, we found that the eCB system controls gut permeability and adipogenesis. We also show that LPS acts as a master switch to control adipose tissue metabolism both in vivo and ex vivo by blocking cannabinoid-driven adipogenesis. These data indicate that gut microbiota determine adipose tissue physiology through LPS-eCB system regulatory loops and may have critical functions in adipose tissue plasticity during obesity.

Metabolomic and transcriptomic stress response of Escherichia coli

Abstract: Environmental fluctuations lead to a rapid adjustment of the physiology of Escherichia coli, necessitating changes on every level of the underlying cellular and molecular network. Thus far, the majority of global analyses of E. coli stress responses have been limited to just one level, gene expression. Here, we incorporate the metabolite composition together with gene expression data to provide a more comprehensive insight on system level stress adjustments by describing detailed time-resolved E. coli response to five different perturbations (cold, heat, oxidative stress, lactose diauxie, and stationary phase). The metabolite response is more specific as compared with the general response observed on the transcript level and is reflected by much higher specificity during the early stress adaptation phase and when comparing the stationary phase response to other perturbations. Despite these differences, the response on both levels still follows the same dynamics and general strategy of energy conservation as reflected by rapid decrease of central carbon metabolism intermediates coinciding with downregulation of genes related to cell growth. Application of co-clustering and canonical correlation analysis on combined metabolite and transcript data identified a number of significant condition-dependent associations between metabolites and transcripts. The results confirm and extend existing models about co-regulation between gene expression and metabolites demonstrating the power of integrated systems oriented analysis.

Computational reconstruction of tissue-specific metabolic models: application to human liver metabolism

Abstract: The computational study of human metabolism has been advanced with the advent of the first generic (non-tissue specific) stoichiometric model of human metabolism. In this study, we present a new algorithm for rapid reconstruction of tissue-specific genome-scale models of human metabolism. The algorithm generates a tissue-specific model from the generic human model by integrating a variety of tissue-specific molecular data sources, including literature-based knowledge, transcriptomic, proteomic, metabolomic and phenotypic data. Applying the algorithm, we constructed the first genome-scale stoichiometric model of hepatic metabolism. The model is verified using standard cross-validation procedures, and through its ability to carry out hepatic metabolic functions. The model's flux predictions correlate with flux measurements across a variety of hormonal and dietary conditions, and improve upon the predictive performance obtained using the original, generic human model (prediction accuracy of 0.67 versus 0.46). Finally, the model better predicts biomarker changes in genetic metabolic disorders than the generic human model (accuracy of 0.67 versus 0.59). The approach presented can be used to construct other human tissue-specific models, and be applied to other organisms.

Defining the transcriptome and proteome in three functionally different human cell lines

Abstract: An essential question in human biology is how cells and tissues differ in gene and protein expression and how these differences delineate specific biological function. Here, we have performed a global analysis of both mRNA and protein levels based on sequence-based transcriptome analysis (RNA-seq), SILAC-based mass spectrometry analysis and antibody-based confocal microscopy. The study was performed in three functionally different human cell lines and based on the global analysis, we estimated the fractions of mRNA and protein that are cell specific or expressed at similar/different levels in the cell lines. A highly ubiquitous RNA expression was found with >60% of the gene products detected in all cells. The changes of mRNA and protein levels in the cell lines using SILAC and RNA ratios show high correlations, even though the genome-wide dynamic range is substantially higher for the proteins as compared with the transcripts. Large general differences in abundance for proteins from various functional classes are observed and, in general, the cell-type specific proteins are low abundant and highly enriched for cell-surface proteins. Thus, this study shows a path to characterize the transcriptome and proteome in human cells from different origins.

A human B‐cell interactome identifies MYB and FOXM1 as master regulators of proliferation in germinal centers

Abstract: Assembly of a transcriptional and post-translational molecular interaction network in B cells, the human B-cell interactome (HBCI), reveals a hierarchical, transcriptional control module, where MYB and FOXM1 act as synergistic master regulators of proliferation in the germinal center (GC). Eighty percent of genes jointly regulated by these transcription factors are activated in the GC, including those encoding proteins in a complex regulating DNA pre-replication, replication, and mitosis. These results indicate that the HBCI analysis can be used for the identification of determinants of major human cell phenotypes and provides a paradigm of general applicability to normal and pathologic tissues.

Target mRNA abundance dilutes microRNA and siRNA activity

Abstract: Post-transcriptional regulation by microRNAs and siRNAs depends not only on characteristics of individual binding sites in target mRNA molecules, but also on system-level properties such as overall molecular concentrations. We hypothesize that an intracellular pool of microRNAs/siRNAs faced with a larger number of available predicted target transcripts will downregulate each individual target gene to a lesser extent. To test this hypothesis, we analyzed mRNA expression change from 178 microRNA and siRNA transfection experiments in two cell lines. Wefind that downregulation of particular genes mediated by microRNAs and siRNAs indeed varies with the total concentration of available target transcripts. We conclude that to interpret and design experiments involving gene regulation bysmall RNAs,global properties,such as target mRNA abundance,need to be considered in addition to local determinants. We propose that analysis of microRNA/siRNA targeting would benefit from a more quantitative definition, rather than simple categorization of genes as ‘target’ or ‘not a target.’ Our results are important for understanding microRNA regulation and may also have implications for siRNA design and small RNA therapeutics. Molecular Systems Biology 6: 363; published online 20 April 2010; doi:10.1038/msb.2010.24 Subject Categories: functional genomics; RNA

Emergent cooperation in microbial metabolism.

Abstract: Mixed microbial communities exhibit emergent biochemical properties not found in clonal monocultures. We report a new type of synthetic genetic interaction, synthetic mutualism in trans (SMIT), in which certain pairs of auxotrophic Escherichia coli mutants complement one another's growth by cross-feeding essential metabolites. We find significant metabolic synergy in 17% of 1035 such pairs tested, with SMIT partners identified throughout the metabolic network. Cooperative phenotypes show more growth on average by aiding the proliferation of their conjugate partner, thereby expanding the source of their own essential metabolites. We construct a quantitative, predictive, framework for describing SMIT interactions as governed by stoichiometric models of the metabolic networks of the interacting strains.

Targeted interactomics reveals a complex core cell cycle machinery in Arabidopsis thaliana

Abstract: Cell proliferation is the main driving force for plant growth. Although genome sequence analysis revealed a high number of cell cycle genes in plants, little is known about the molecular complexes steering cell division. In a targeted proteomics approach, we mapped the core complex machinery at the heart of the Arabidopsis thaliana cell cycle control. Besides a central regulatory network of core complexes, we distinguished a peripheral network that links the core machinery to up- and downstream pathways. Over 100 new candidate cell cycle proteins were predicted and an in-depth biological interpretation demonstrated the hypothesis-generating power of the interaction data. The data set provided a comprehensive view on heterodimeric cyclin-dependent kinase (CDK)–cyclin complexes in plants. For the first time, inhibitory proteins of plant-specific B-type CDKs were discovered and the anaphase-promoting complex was characterized and extended. Important conclusions were that mitotic A- and B-type cyclins form complexes with the plant-specific B-type CDKs and not with CDKA;1, and that D-type cyclins and S-phase-specific A-type cyclins seem to be associated exclusively with CDKA;1. Furthermore, we could show that plants have evolved a combinatorial toolkit consisting of at least 92 different CDK–cyclin complex variants, which strongly underscores the functional diversification among the large family of cyclins and reflects the pivotal role of cell cycle regulation in the developmental plasticity of plants.

Coupling governs entrainment range of circadian clocks.

Abstract: Circadian clocks are endogenous oscillators driving daily rhythms in physiology and behavior. Synchronization of these timers to environmental light–dark cycles (‘entrainment’) is crucial for an organism’s fitness. Little is known about which oscillator qualities determine entrainment, i.e., entrainment range, phase and amplitude. In a systematic theoretical and experimental study, we uncovered these qualities for circadian oscillators in the suprachiasmatic nucleus (SCN—the master clock in mammals) and the lung (a peripheral clock): (i) the ratio between stimulus (zeitgeber) strength and oscillator amplitude and (ii) the rigidity of the oscillatory system (relaxation rate upon perturbation) determine entrainment properties. Coupling among oscillators affects both qualities resulting in increased amplitude and rigidity. These principles explain our experimental findings that lung clocks entrain to extreme zeitgeber cycles, whereas SCN clocks do not. We confirmed our theoretical predictions by showing that pharmacological inhibition of coupling in the SCN leads to larger ranges of entrainment. These differences between master and the peripheral clocks suggest that coupling-induced rigidity in the SCN filters environmental noise to create a robust circadian system.

Evolution, genomic analysis, and reconstruction of isobutanol tolerance in Escherichia coli

Abstract: Escherichia coli has been engineered to produce isobutanol, with titers reaching greater than the toxicity level. However, the specific effects of isobutanol on the cell have never been fully understood. Here, we aim to identify genotype–phenotype relationships in isobutanol response. An isobutanol-tolerant mutant was isolated with serial transfers. Using whole-genome sequencing followed by gene repair and knockout, we identified five mutations (acrA, gatY, tnaA, yhbJ, and marCRAB) that were primarily responsible for the increased isobutanol tolerance. We successfully reconstructed the tolerance phenotype by combining deletions of these five loci, and identified glucosamine-6-phosphate as an important metabolite for isobutanol tolerance, which presumably enhanced membrane synthesis. The isobutanol-tolerant mutants also show increased tolerance to n-butanol and 2-methyl-1-butanol, but showed no improvement in ethanol tolerance and higher sensitivity to hexane and chloramphenicol than the parental strain. These results suggest that C4, C5 alcohol stress impacts the cell differently compared with the general solvent or antibiotic stresses. Interestingly, improved isobutanol tolerance did not increase the final titer of isobutanol production.

HepatoNet1: a comprehensive metabolic reconstruction of the human hepatocyte for the analysis of liver physiology

Abstract: We present HepatoNet1, the first reconstruction of a comprehensive metabolic network of the human hepatocyte that is shown to accomplish a large canon of known metabolic liver functions. The network comprises 777 metabolites in six intracellular and two extracellular compartments and 2539 reactions, including 1466 transport reactions. It is based on the manual evaluation of >1500 original scientific research publications to warrant a high-quality evidence-based model. The final network is the result of an iterative process of data compilation and rigorous computational testing of network functionality by means of constraint-based modeling techniques. Taking the hepatic detoxification of ammonia as an example, we show how the availability of nutrients and oxygen may modulate the interplay of various metabolic pathways to allow an efficient response of the liver to perturbations of the homeostasis of blood compounds.

Metabonomic, transcriptomic, and genomic variation of a population cohort

Abstract: Comprehensive characterization of human tissues promises novel insights into the biological architecture of human diseases and traits. We assessed metabonomic, transcriptomic, and genomic variation for a large population-based cohort from the capital region of Finland. Network analyses identified a set of highly correlated genes, the lipid‐leukocyte (LL) module, as having a prominent role in over 80 serum metabolites (of 134 measures quantified), including lipoprotein subclasses, lipids, and amino acids. Concurrent association with immune response markers suggested the LL module as a possible link between inflammation, metabolism, and adiposity. Further, genomic variation was used to generate a directed network and infer LL module’s largely reactive nature to metabolites. Finally, gene co-expression in circulating leukocytes was shown to be dependent on serum metabolite concentrations, providing evidence for the hypothesis that the coherence of molecular networks themselves is conditional on environmental factors. These findings show the importance and opportunityof systematic molecular investigation of human population samples.To facilitate and encourage this investigation, the metabonomic, transcriptomic, and genomic data used in this study have been made available as a resource for the research community. Molecular Systems Biology 6: 441; published online 21 December 2010; doi:10.1038/msb.2010.93 Subject Categories: functional genomics; metabolic and regulatory networks

Insight into human alveolar macrophage and M. tuberculosis interactions via metabolic reconstructions

Abstract: Metabolic coupling of Mycobacterium tuberculosis to its host is foundational to its pathogenesis. Computational genome-scale metabolic models have shown utility in integrating -omic as well as physiologic data for systemic, mechanistic analysis of metabolism. To date, integrative analysis of host-pathogen interactions using in silico mass-balanced, genome-scale models has not been performed. We, therefore, constructed a cell-specific alveolar macrophage model, iAB-AMO-1410, from the global human metabolic reconstruction, Recon 1. The model successfully predicted experimentally verified ATP and nitric oxide production rates in macrophages. This model was then integrated with an M. tuberculosis H37Rv model, iNJ661, to build an integrated host-pathogen genome-scale reconstruction, iAB-AMO-1410-Mt-661. The integrated host-pathogen network enables simulation of the metabolic changes during infection. The resulting reaction activity and gene essentiality targets of the integrated model represent an altered infectious state. High-throughput data from infected macrophages were mapped onto the host-pathogen network and were able to describe three distinct pathological states. Integrated host-pathogen reconstructions thus form a foundation upon which understanding the biology and pathophysiology of infections can be developed.

Bacterial adaptation through distributed sensing of metabolic fluxes

Abstract: The recognition of carbon sources and the regulatory adjustments to recognized changes are of particular importance for bacterial survival in fluctuating environments. Despite a thorough knowledge base of Escherichia coli's central metabolism and its regulation, fundamental aspects of the employed sensing and regulatory adjustment mechanisms remain unclear. In this paper, using a differential equation model that couples enzymatic and transcriptional regulation of E. coli's central metabolism, we show that the interplay of known interactions explains in molecular-level detail the system-wide adjustments of metabolic operation between glycolytic and gluconeogenic carbon sources. We show that these adaptations are enabled by an indirect recognition of carbon sources through a mechanism we termed distributed sensing of intracellular metabolic fluxes. This mechanism uses two general motifs to establish flux-signaling metabolites, whose bindings to transcription factors form flux sensors. As these sensors are embedded in global feedback loop architectures, closed-loop self-regulation can emerge within metabolism itself and therefore, metabolic operation may adapt itself autonomously (not requiring upstream sensing and signaling) to fluctuating carbon sources.

A modular gradient-sensing network for chemotaxis in Escherichia coli revealed by responses to time-varying stimuli.

Abstract: The Escherichia coli chemotaxis-signaling pathway computes time derivatives of chemoeffector concentrations. This network features modules for signal reception/amplification and robust adaptation, with sensing of chemoeffector gradients determined by the way in which these modules are coupled in vivo. We characterized these modules and their coupling by using fluorescence resonance energy transfer to measure intracellular responses to time-varying stimuli. Receptor sensitivity was characterized by step stimuli, the gradient sensitivity by exponential ramp stimuli, and the frequency response by exponential sine-wave stimuli. Analysis of these data revealed the structure of the feedback transfer function linking the amplification and adaptation modules. Feedback near steady state was found to be weak, consistent with strong fluctuations and slow recovery from small perturbations. Gradient sensitivity and frequency response both depended strongly on temperature. We found that time derivatives can be computed by the chemotaxis system for input frequencies below 0.006Hz at 221C and below 0.018Hz at 321C. Our results show how dynamic input‐output measurements, time honored in physiology, can serve as powerful tools in deciphering cell-signaling mechanisms. Molecular Systems Biology 6: 382; published online 22 June 2010; doi:10.1038/msb.2010.37 Subject Categories: signal transduction

A comprehensive map of the mTOR signaling network.

Abstract: The mammalian target of rapamycin (mTOR) is a central regulator of cell growth and proliferation. mTOR signaling is frequently dysregulated in oncogenic cells, and thus an attractive target for anticancer therapy. Using CellDesigner, a modeling support software for graphical notation, we present herein a comprehensive map of the mTOR signaling network, which includes 964 species connected by 777 reactions. The map complies with both the systems biology markup language (SBML) and graphical notation (SBGN) for computational analysis and graphical representation, respectively. As captured in the mTOR map, we review and discuss our current understanding of the mTOR signaling network and highlight the impact of mTOR feedback and crosstalk regulations on drug-based cancer therapy. This map is available on the Payao platform, a Web 2.0 based community-wide interactive process for creating more accurate and information-rich databases. Thus, this comprehensive map of the mTOR network will serve as a tool to facilitate systems-level study of up-to-date mTOR network components and signaling events toward the discovery of novel regulatory processes and therapeutic strategies for cancer.

Single-cell quantification of IL-2 response by effector and regulatory T cells reveals critical plasticity in immune response.

Abstract: Understanding how the immune system decides between tolerance and activation by antigens requires addressing cytokine regulation as a highly dynamic process. We quantified the dynamics of interleukin-2 (IL-2) signaling in a population of T cells during an immune response by combining in silico modeling and single-cell measurements in vitro. We demonstrate that IL-2 receptor expression levels vary widely among T cells creating a large variability in the ability of the individual cells to consume, produce and participate in IL-2 signaling within the population. Our model reveals that at the population level, these heterogeneous cells are engaged in a tug-of-war for IL-2 between regulatory (Treg) and effector (Teff) T cells, whereby access to IL-2 can either increase the survival of Teff cells or the suppressive capacity of Treg cells. This tug-of-war is the mechanism enforcing, at the systems level, a core function of Treg cells, namely the specific suppression of survival signals for weakly activated Teff cells but not for strongly activated cells. Our integrated model yields quantitative, experimentally validated predictions for the manipulation of Treg suppression.

An atlas of gene regulatory networks reveals multiple three-gene mechanisms for interpreting morphogen gradients

Abstract: The interpretation of morphogen gradients is a pivotal concept in developmental biology, and several mechanisms have been proposed to explain how gene regulatory networks (GRNs) achieve concentration-dependent responses. However, the number of different mechanisms that may exist for cells to interpret morphogens, and the importance of design features such as feedback or local cell–cell communication, is unclear. A complete understanding of such systems will require going beyond a case-by-case analysis of real morphogen interpretation mechanisms and mapping out a complete GRN ‘design space.’ Here, we generate a first atlas of design space for GRNs capable of patterning a homogeneous field of cells into discrete gene expression domains by interpreting a fixed morphogen gradient. We uncover multiple very distinct mechanisms distributed discretely across the atlas, thereby expanding the repertoire of morphogen interpretation network motifs. Analyzing this diverse collection of mechanisms also allows us to predict that local cell–cell communication will rarely be responsible for the basic dose-dependent response of morphogen

microRNA‐122 as a regulator of mitochondrial metabolic gene network in hepatocellular carcinoma

Abstract: Tumorigenesis involves multistep genetic alterations. To elucidate the microRNA (miRNA)–gene interaction network in carcinogenesis, we examined their genome-wide expression profiles in 96 pairs of tumor/non-tumor tissues from hepatocellular carcinoma (HCC). Comprehensive analysis of the coordinate expression of miRNAs and mRNAs reveals that miR-122 is under-expressed in HCC and that increased expression of miR-122 seed-matched genes leads to a loss of mitochondrial metabolic function. Furthermore, the miR-122 secondary targets, which decrease in expression, are good prognostic markers for HCC. Transcriptome profiling data from additional 180 HCC and 40 liver cirrhotic patients in the same cohort were used to confirm the anti-correlation of miR-122 primary and secondary target gene sets. The HCC findings can be recapitulated in mouse liver by silencing miR-122 with antagomir treatment followed by gene-expression microarray analysis. In vitro miR-122 data further provided a direct link between induction of miR-122-controlled genes and impairment of mitochondrial metabolism. In conclusion, miR-122 regulates mitochondrial metabolism and its loss may be detrimental to sustaining critical liver function and contribute to morbidity and mortality of liver cancer patients.

Integrative model of genomic factors for determining binding site selection by estrogen receptor-α.

Abstract: A major question in transcription factor (TF) biology is why a TF binds to only a small fraction of motif eligible binding sites in the genome. Using the estrogen receptor-α as a model system, we sought to explicitly define parameters that determine TF-binding site selection. By examining 12 genetic and epigenetic parameters, we find that an energetically favorable estrogen response element (ERE) motif sequence, co-occupancy by the TF FOXA1, the presence of the H3K4me1 mark and an open chromatin configuration in the pre-ligand state provide specificity for ER binding. These factors can model estrogen-induced ER binding with high accuracy (ROC-AUC=0.95 and 0.88 using different genomic backgrounds). Moreover, when assessed in another estrogen-responsive cell line, this model was highly predictive for ERα binding (ROC-AUC=0.86). Variance in binding site selection between MCF-7 and T47D resides in sites with suboptimal ERE motifs, but modulated by the chromatin configuration. These results suggest a definable interplay between sequence motifs and local chromatin in selecting TF binding.

Data assimilation constrains new connections and components in a complex, eukaryotic circadian clock model

Abstract: Circadian clocks generate 24-h rhythms that are entrained by the day/night cycle. Clock circuits include several light inputs and interlocked feedback loops, with complex dynamics. Multiple biological components can contribute to each part of the circuit in higher organisms. Mechanistic models with morning, evening and central feedback loops have provided a heuristic framework for the clock in plants, but were based on transcriptional control. Here, we model observed, posttranscriptional and post-translational regulation and constrain many parameter values based on experimental data. The model’s feedback circuit is revised and now includes PSEUDO-RESPONSE REGULATOR 7 (PRR7) andZEITLUPE. The revised model matches data in varying environments and mutants, and gains robustness to parameter variation. Our results suggest that the activation of important morning-expressed genes follows their release from a night inhibitor (NI). Experiments inspired by the new model support the predicted NI function and show that the PRR5 gene contributes to the NI. The multiple PRR genes of Arabidopsis uncouple events in the late night from light-driven responses in the day, increasing the flexibility of rhythmic regulation. Molecular Systems Biology 6: 416; published online 21 September 2010; doi:10.1038/msb.2010.69 Subject Categories: metabolic and regulatory networks; simulation and data analysis

Clustering phenotype populations by genome-wide RNAi and multiparametric imaging.

Abstract: Genetic screens for phenotypic similarity have made key contributions to associating genes with biological processes. With RNA interference (RNAi), highly parallel phenotyping of loss-of-function effects in cells has become feasible. One of the current challenges however is the computational categorization of visual phenotypes and the prediction of biological function and processes. In this study, we describe a combined computational and experimental approach to discover novel gene functions and explore functional relationships. We performed a genome-wide RNAi screen in human cells and used quantitative descriptors derived from high-throughput imaging to generate multiparametric phenotypic profiles. We show that profiles predicted functions of genes by phenotypic similarity. Specifically, we examined several candidates including the largely uncharacterized gene DONSON, which shared phenotype similarity with known factors of DNA damage response (DDR) and genomic integrity. Experimental evidence supports that DONSON is a novel centrosomal protein required for DDR signalling and genomic integrity. Multiparametric phenotyping by automated imaging and computational annotation is a powerful method for functional discovery and mapping the landscape of phenotypic responses to cellular perturbations.

Cell-type specific analysis of translating RNAs in developing flowers reveals new levels of control

Abstract: Determining both the expression levels of mRNA and the regulation of its translation is important in understanding specialized cell functions. In this study, we describe both the expression profiles of cells within spatiotemporal domains of the Arabidopsis thaliana flower and the post-transcriptional regulation of these mRNAs, at nucleotide resolution. We express a tagged ribosomal protein under the promoters of three master regulators of flower development. By precipitating tagged polysomes, we isolated cell type-specific mRNAs that are probably translating, and quantified those mRNAs through deep sequencing. Cell type comparisons identified known cell-specific transcripts and uncovered many new ones, from which we inferred cell type-specific hormone responses, promoter motifs and coexpressed cognate binding factor candidates, and splicing isoforms. By comparing translating mRNAs with steady-state overall transcripts, we found evidence for widespread post-transcriptional regulation at both the intron splicing and translational stages. Sequence analyses identified structural features associated with each step. Finally, we identified a new class of noncoding RNAs associated with polysomes. Findings from our profiling lead to new hypotheses in the understanding of flower development.

Cross-species chemogenomic profiling reveals evolutionarily conserved drug mode of action

Abstract: We present a cross-species chemogenomic screening platform using libraries of haploid deletion mutants from two yeast species, Saccharomyces cerevisiae and Schizosaccharomyces pombe .W e screened a set of compounds of known and unknown mode of action (MoA) and derived quantitative drug scores (or D-scores), identifying mutants that are either sensitive or resistant to particular compounds. We found that compound–functional module relationships are more conserved than individual compound–gene interactions between these two species. Furthermore, we observed that combining data from both species allows for more accurate prediction of MoA. Finally, using this platform, we identified a novel small molecule that acts as a DNA damaging agent and demonstrate

Tradeoff between enzyme and metabolite efficiency maintains metabolic homeostasis upon perturbations in enzyme capacity

Abstract: What is the relationship between enzymes and metabolites, the two major constituents of metabolic networks? We propose three alternative relationships between enzyme capacity and metabolite concentration alterations based on a Michaelis–Menten kinetic; that is enzyme capacities, metabolite concentrations, or both could limit the metabolic reaction rates. These relationships imply different correlations between changes in enzyme capacity and metabolite concentration, which we tested by quantifying metabolite, transcript, and enzyme abundances upon local (single-enzyme modulation) and global (GCR2 transcription factor mutant) perturbations in Saccharomyces cerevisiae. Our results reveal an inverse relationship between fold-changes in substrate metabolites and their catalyzing enzymes. These data provide evidence for the hypothesis that reaction rates are jointly limited by enzyme capacity and metabolite concentration. Hence, alteration in one network constituent can be efficiently buffered by converse alterations in the other constituent, implying a passive mechanism to maintain metabolic homeostasis upon perturbations in enzyme capacity.

Regulatory and metabolic rewiring during laboratory evolution of ethanol tolerance in E. coli

Abstract: Understanding the genetic basis of adaptation is a central problem in biology. However, revealing the underlying molecular mechanisms has been challenging as changes in fitness may result from perturbations to many pathways, any of which may contribute relatively little. We have developed a combined experimental/computational framework to address this problem and used it to understand the genetic basis of ethanol tolerance in Escherichia coli. We used fitness profiling to measure the consequences of single-locus perturbations in the context of ethanol exposure. A module-level computational analysis was then used to reveal the organization of the contributing loci into cellular processes and regulatory pathways (e.g. osmoregulation and cell-wall biogenesis) whose modifications significantly affect ethanol tolerance. Strikingly, we discovered that a dominant component of adaptation involves metabolic rewiring that boosts intracellular ethanol degradation and assimilation. Through phenotypic and metabolomic analysis of laboratory-evolved ethanol-tolerant strains, we investigated naturally accessible pathways of ethanol tolerance. Remarkably, these laboratory-evolved strains, by and large, follow the same adaptive paths as inferred from our coarse-grained search of the fitness landscape.