Edda Klipp

Condition-Dependent Transcriptome Reveals High-Level Regulatory Architecture in Bacillus subtilis

Outside In Acquisition and analysis of large data sets promises to move us toward a greater understanding of the mechanisms by which biological systems are dynamically regulated to respond to external cues. Now, two papers explore the responses of a bacterium to changing nutritional conditions (see the Perspective by Chalancon et al.). Nicolas et al. (p. 1103) measured transcriptional regulation for more than 100 different conditions. Greater amounts of antisense RNA were generated than expected and appeared to be produced by alternative RNA polymerase targeting subunits called sigma factors. One transition, from malate to glucose as the primary nutrient, was studied in more detail by Buescher et al. (p. 1099) who monitored RNA abundance, promoter activity in live cells, protein abundance, and absolute concentrations of intracellular and extracellular metabolites. In this case, the bacteria responded rapidly and largely without transcriptional changes to life on malate, but only slowly adapted to use glucose, a shift that required changes in nearly half the transcription network. These data offer an initial understanding of why certain regulatory strategies may be favored during evolution of dynamic control systems. A horizontal analysis reveals the breadth of genes turned on and off as nutrients change. Bacteria adapt to environmental stimuli by adjusting their transcriptomes in a complex manner, the full potential of which has yet to be established for any individual bacterial species. Here, we report the transcriptomes of Bacillus subtilis exposed to a wide range of environmental and nutritional conditions that the organism might encounter in nature. We comprehensively mapped transcription units (TUs) and grouped 2935 promoters into regulons controlled by various RNA polymerase sigma factors, accounting for ~66% of the observed variance in transcriptional activity. This global classification of promoters and detailed description of TUs revealed that a large proportion of the detected antisense RNAs arose from potentially spurious transcription initiation by alternative sigma factors and from imperfect control of transcription termination.

Integrative model of the response of yeast to osmotic shock

Integration of experimental studies with mathematical modeling allows insight into systems properties, prediction of perturbation effects and generation of hypotheses for further research. We present a comprehensive mathematical description of the cellular response of yeast to hyperosmotic shock. The model integrates a biochemical reaction network comprising receptor stimulation, mitogen-activated protein kinase cascade dynamics, activation of gene expression and adaptation of cellular metabolism with a thermodynamic description of volume regulation and osmotic pressure. Simulations agree well with experimental results obtained under different stress conditions or with specific mutants. The model is predictive since it suggests previously unrecognized features of the system with respect to osmolyte accumulation and feedback control, as confirmed with experiments. The mathematical description presented is a valuable tool for future studies on osmoregulation in yeast and—with appropriate modifications—other organisms. It also serves as a starting point for a comprehensive description of cellular signaling.

Dynamic rerouting of the carbohydrate flux is key to counteracting oxidative stress

Background Eukaryotic cells have evolved various response mechanisms to counteract the deleterious consequences of oxidative stress. Among these processes, metabolic alterations seem to play an important role. Results We recently discovered that yeast cells with reduced activity of the key glycolytic enzyme triosephosphate isomerase exhibit an increased resistance to the thiol-oxidizing reagent diamide. Here we show that this phenotype is conserved in Caenorhabditis elegans and that the underlying mechanism is based on a redirection of the metabolic flux from glycolysis to the pentose phosphate pathway, altering the redox equilibrium of the cytoplasmic NADP(H) pool. Remarkably, another key glycolytic enzyme, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), is known to be inactivated in response to various oxidant treatments, and we show that this provokes a similar redirection of the metabolic flux. Conclusion The naturally occurring inactivation of GAPDH functions as a metabolic switch for rerouting the carbohydrate flux to counteract oxidative stress. As a consequence, altering the homoeostasis of cytoplasmic metabolites is a fundamental mechanism for balancing the redox state of eukaryotic cells under stress conditions.

Global Network Reorganization During Dynamic Adaptations of Bacillus subtilis Metabolism

Outside In Acquisition and analysis of large data sets promises to move us toward a greater understanding of the mechanisms by which biological systems are dynamically regulated to respond to external cues. Now, two papers explore the responses of a bacterium to changing nutritional conditions (see the Perspective by Chalancon et al.). Nicolas et al. (p. 1103) measured transcriptional regulation for more than 100 different conditions. Greater amounts of antisense RNA were generated than expected and appeared to be produced by alternative RNA polymerase targeting subunits called sigma factors. One transition, from malate to glucose as the primary nutrient, was studied in more detail by Buescher et al. (p. 1099) who monitored RNA abundance, promoter activity in live cells, protein abundance, and absolute concentrations of intracellular and extracellular metabolites. In this case, the bacteria responded rapidly and largely without transcriptional changes to life on malate, but only slowly adapted to use glucose, a shift that required changes in nearly half the transcription network. These data offer an initial understanding of why certain regulatory strategies may be favored during evolution of dynamic control systems. A vertical analysis reveals that a simple switch of one food for another evokes changes at many levels. Adaptation of cells to environmental changes requires dynamic interactions between metabolic and regulatory networks, but studies typically address only one or a few layers of regulation. For nutritional shifts between two preferred carbon sources of Bacillus subtilis, we combined statistical and model-based data analyses of dynamic transcript, protein, and metabolite abundances and promoter activities. Adaptation to malate was rapid and primarily controlled posttranscriptionally compared with the slow, mainly transcriptionally controlled adaptation to glucose that entailed nearly half of the known transcription regulation network. Interactions across multiple levels of regulation were involved in adaptive changes that could also be achieved by controlling single genes. Our analysis suggests that global trade-offs and evolutionary constraints provide incentives to favor complex control programs.

Triple-negative breast cancer: Present challenges and new perspectives

Triple‐negative breast cancers (TNBC), characterized by absence of estrogen receptor (ER), progesterone receptor (PR) and lack of overexpression of human epidermal growth factor receptor 2 (HER2), are typically associated with poor prognosis, due to aggressive tumor phenotype(s), only partial response to chemotherapy and present lack of clinically established targeted therapies. Advances in the design of individualized strategies for treatment of TNBC patients require further elucidation, by combined ‘omics’ approaches, of the molecular mechanisms underlying TNBC phenotypic heterogeneity, and the still poorly understood association of TNBC with BRCA1 mutations. An overview is here presented on TNBC profiling in terms of expression signatures, within the functional genomic breast tumor classification, and ongoing efforts toward identification of new therapy targets and bioimaging markers. Due to the complexity of aberrant molecular patterns involved in expression, pathological progression and biological/clinical heterogeneity, the search for novel TNBC biomarkers and therapy targets requires collection of multi‐dimensional data sets, use of robust multivariate data analysis techniques and development of innovative systems biology approaches.

Selective benefits of damage partitioning in unicellular systems and its effects on aging

Cytokinesis in unicellular organisms sometimes entails a division of labor between cells leading to lineage-specific aging. To investigate the potential benefits of asymmetrical cytokinesis, we created a mathematical model to simulate the robustness and fitness of dividing systems displaying different degrees of damage segregation and size asymmetries. The model suggests that systems dividing asymmetrically (size-wise) or displaying damage segregation can withstand higher degrees of damage before entering clonal senescence. When considering population fitness, a system producing different-sized progeny like budding yeast is predicted to benefit from damage retention only at high damage propagation rates. In contrast, the fitness of a system of equal-sized progeny is enhanced by damage segregation regardless of damage propagation rates, suggesting that damage partitioning may also provide an evolutionary advantage in systems dividing by binary fission. Indeed, by using Schizosaccharomyces pombe as a model, we experimentally demonstrate that damaged proteins are unevenly partitioned during cytokinesis and the damage-enriched sibling suffers from a prolonged generation time and accelerated aging. This damage retention in S. pombe is, like in Saccharomyces cerevisiae, Sir2p- and cytoskeleton-dependent, suggesting this to be an evolutionarily conserved mechanism. We suggest that sibling-specific aging may be a result of the strong selective advantage of damage segregation, which may be more common in nature than previously anticipated.

Modelling the dynamics of the yeast pheromone pathway

We present a mathematical model of the dynamics of the pheromone pathways in haploid yeast cells of mating type MATa after stimulation with pheromone α‐factor. The model consists of a set of differential equations and describes the dynamics of signal transduction from the receptor via several steps, including a G protein and a scaffold MAP kinase cascade, up to changes in the gene expression after pheromone stimulation in terms of biochemical changes (complex formations, phosphorylations, etc.). The parameters entering the models have been taken from the literature or adapted to observed time courses or behaviour. Using this model we can follow the time course of the various complex formation processes and of the phosphorylation states of the proteins involved. Furthermore, we can explain the phenotype of more than a dozen well‐characterized mutants and also the graded response of yeast cells to varying concentrations of the stimulating pheromone. Copyright

Cell Size at S Phase Initiation: An Emergent Property of the G1/S Network

The eukaryotic cell cycle is the repeated sequence of events that enable the division of a cell into two daughter cells. It is divided into four phases: G1, S, G2, and M. Passage through the cell cycle is strictly regulated by a molecular interaction network, which involves the periodic synthesis and destruction of cyclins that bind and activate cyclin-dependent kinases that are present in nonlimiting amounts. Cyclin-dependent kinase inhibitors contribute to cell cycle control. Budding yeast is an established model organism for cell cycle studies, and several mathematical models have been proposed for its cell cycle. An area of major relevance in cell cycle control is the G1 to S transition. In any given growth condition, it is characterized by the requirement of a specific, critical cell size, PS, to enter S phase. The molecular basis of this control is still under discussion. The authors report a mathematical model of the G1 to S network that newly takes into account nucleo/cytoplasmic localization, the role of the cyclin-dependent kinase Sic1 in facilitating nuclear import of its cognate Cdk1-Clb5, Whi5 control, and carbon source regulation of Sic1 and Sic1-containing complexes. The model was implemented by a set of ordinary differential equations that describe the temporal change of the concentration of the involved proteins and protein complexes. The model was tested by simulation in several genetic and nutritional setups and was found to be neatly consistent with experimental data. To estimate PS, the authors developed a hybrid model including a probabilistic component for firing of DNA replication origins. Sensitivity analysis of PS provides a novel relevant conclusion: PS is an emergent property of the G1 to S network that strongly depends on growth rate.

Minimum Information About a Simulation Experiment (MIASE).

Reproducibility of experiments is a basic requirement for science. Minimum Information (MI) guidelines have proved a helpful means of enabling reuse of existing work in modern biology. The Minimum Information Required in the Annotation of Models (MIRIAM) guidelines promote the exchange and reuse of biochemical computational models. However, information about a model alone is not sufficient to enable its efficient reuse in a computational setting. Advanced numerical algorithms and complex modeling workflows used in modern computational biology make reproduction of simulations difficult. It is therefore essential to define the core information necessary to perform simulations of those models. The Minimum Information About a Simulation Experiment (MIASE, Glossary in Box 1) describes the minimal set of information that must be provided to make the description of a simulation experiment available to others. It includes the list of models to use and their modifications, all the simulation procedures to apply and in which order, the processing of the raw numerical results, and the description of the final output. MIASE allows for the reproduction of any simulation experiment. The provision of this information, along with a set of required models, guarantees that the simulation experiment represents the intention of the original authors. Following MIASE guidelines will thus improve the quality of scientific reporting, and will also allow collaborative, more distributed efforts in computational modeling and simulation of biological processes.

Mathematical modeling of intracellular signaling pathways

Dynamic modeling and simulation of signal transduction pathways is an important topic in systems biology and is obtaining growing attention from researchers with experimental or theoretical background. Here we review attempts to analyze and model specific signaling systems. We review the structure of recurrent building blocks of signaling pathways and their integration into more comprehensive models, which enables the understanding of complex cellular processes. The variety of mechanisms found and modeling techniques used are illustrated with models of different signaling pathways. Focusing on the close interplay between experimental investigation of pathways and the mathematical representations of cellular dynamics, we discuss challenges and perspectives that emerge in studies of signaling systems.