Wolfram Liebermeister

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.

A comprehensive library of fluorescent transcriptional reporters for Escherichia coli

E. coli is widely used for systems biology research; there exists a need, however, for tools that can be used to accurately and comprehensively measure expression dynamics in individual living cells. To address this we present a library of transcriptional fusions of gfp to each of about 2,000 different promoters in E. coli K12, covering the great majority of the promoters in the organism. Each promoter fusion is expressed from a low-copy plasmid. We demonstrate that this library can be used to obtain highly accurate dynamic measurements of promoter activity on a genomic scale, in a glucose-lactose diauxic shift experiment. The library allowed detection of about 80 previously uncharacterized transcription units in E. coli, including putative internal promoters within previously known operons, such as the lac operon. This library can serve as a tool for accurate, high-resolution analysis of transcription networks in living E. coli cells.

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.

Glycolytic strategy as a tradeoff between energy yield and protein cost

Contrary to the textbook portrayal of glycolysis as a single pathway conserved across all domains of life, not all sugar-consuming organisms use the canonical Embden–Meyerhoff–Parnass (EMP) glycolytic pathway. Prokaryotic glucose metabolism is particularly diverse, including several alternative glycolytic pathways, the most common of which is the Entner–Doudoroff (ED) pathway. The prevalence of the ED pathway is puzzling as it produces only one ATP per glucose—half as much as the EMP pathway. We argue that the diversity of prokaryotic glucose metabolism may reflect a tradeoff between a pathway’s energy (ATP) yield and the amount of enzymatic protein required to catalyze pathway flux. We introduce methods for analyzing pathways in terms of thermodynamics and kinetics and show that the ED pathway is expected to require several-fold less enzymatic protein to achieve the same glucose conversion rate as the EMP pathway. Through genomic analysis, we further show that prokaryotes use different glycolytic pathways depending on their energy supply. Specifically, energy-deprived anaerobes overwhelmingly rely upon the higher ATP yield of the EMP pathway, whereas the ED pathway is common among facultative anaerobes and even more common among aerobes. In addition to demonstrating how protein costs can explain the use of alternative metabolic strategies, this study illustrates a direct connection between an organism’s environment and the thermodynamic and biochemical properties of the metabolic pathways it employs.

Integrating quantitative proteomics and metabolomics with a genome-scale metabolic network model

Motivation: The availability of modern sequencing techniques has led to a rapid increase in the amount of reconstructed metabolic networks. Using these models as a platform for the analysis of high throughput transcriptomic, proteomic and metabolomic data can provide valuable insight into conditional changes in the metabolic activity of an organism. While transcriptomics and proteomics provide important insights into the hierarchical regulation of metabolic flux, metabolomics shed light on the actual enzyme activity through metabolic regulation and mass action effects. Here we introduce a new method, termed integrative omics-metabolic analysis (IOMA) that quantitatively integrates proteomic and metabolomic data with genome-scale metabolic models, to more accurately predict metabolic flux distributions. The method is formulated as a quadratic programming (QP) problem that seeks a steady-state flux distribution in which flux through reactions with measured proteomic and metabolomic data, is as consistent as possible with kinetically derived flux estimations. Results: IOMA is shown to successfully predict the metabolic state of human erythrocytes (compared to kinetic model simulations), showing a significant advantage over the commonly used methods flux balance analysis and minimization of metabolic adjustment. Thereafter, IOMA is shown to correctly predict metabolic fluxes in Escherichia coli under different gene knockouts for which both metabolomic and proteomic data is available, achieving higher prediction accuracy over the extant methods. Considering the lack of high-throughput flux measurements, while high-throughput metabolomic and proteomic data are becoming readily available, we expect IOMA to significantly contribute to future research of cellular metabolism. Contacts: kerenyiz@post.tau.ac.il; tomersh@cs.technion.ac.il

Spanning high-dimensional expression space using ribosome-binding site combinatorics

Protein levels are a dominant factor shaping natural and synthetic biological systems. Although proper functioning of metabolic pathways relies on precise control of enzyme levels, the experimental ability to balance the levels of many genes in parallel is a major outstanding challenge. Here, we introduce a rapid and modular method to span the expression space of several proteins in parallel. By combinatorially pairing genes with a compact set of ribosome-binding sites, we modulate protein abundance by several orders of magnitude. We demonstrate our strategy by using a synthetic operon containing fluorescent proteins to span a 3D color space. Using the same approach, we modulate a recombinant carotenoid biosynthesis pathway in Escherichia coli to reveal a diversity of phenotypes, each characterized by a distinct carotenoid accumulation profile. In a single combinatorial assembly, we achieve a yield of the industrially valuable compound astaxanthin 4-fold higher than previously reported. The methodology presented here provides an efficient tool for exploring a high-dimensional expression space to locate desirable phenotypes.

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.

Systems biology standards—the community speaks

VOLUME 25 NUMBER 4 APRIL 2007 NATURE BIOTECHNOLOGY available filters and dramatically decreases clogging and breakage. This allows undesired components to be easily eliminated during blood diagnosis. Third, microscale electrophoresis using bacterial cellulose in a low-viscosity polymer matrix results in a dramatic improvement in separation and detection of analytes compared with current techniques. And finally, the unique light amplification properties of Nata de Coco6,7 provide a much greater level of detection sensitivity compared with other light amplification systems. All these properties allow optimized detection of genetic variation by virtue of better separation resolution and higher detection intensity of biomolecules7. Our in vitro results using the methylthiazolydiphenyl-tetrazolium bromide (MTT) assay in the human leukemia cell line THP-1 also indicate that bacterial cellulose does not upregulate CD86 or CD54 expression, suggesting it is unlikely to be allergenic (Fig. 1e). Indeed, bacterial cellulose is already under clinical development as a scaffold for tissueengineered products9. Nata de Coco is just one of numerous nanomaterials found in nature that could be applied in innovative nanotechnology products and devices. Cellulose is the most abundant material in nature and, apart from bacteria, it can be derived from sources as diverse as wood, cotton, animals and plant cell walls. Other abundant natural materials that may mined for nanomaterials include silk, collagen, chitin, chitosan and polylactic acid. Recently, much of the emphasis in nanotechnology has been on synthetic materials, such as carbon nanotubes, nanoparticles (colloidal gold, quantum dots, latex and so on) and inorganic nanomaterials (ZnO, TiO2, silica). For these synthetic nanomaterials to be applied as drug carriers, medical treatments, implants, tissue engineering constructs and cosmetics, extensive investigation will be necessary to establish their biocompatibility, safety, immunogenicity and allergenicity. Although advances in these areas are exciting and should continue to be funded and supported, we should not forget that many materials are available directly from the natural world with interesting properties on the nanoscale. These materials have a proven safety record in humans, biodegradable properties that make them environmentally friendly and in certain cases economic advantages over more sophisticated and expensive products and technologies. Bacterial cellulose is just one attractive source of nanomaterials for use in medical and food applications. But the potential for naturally occurring nanomaterials is boundless. As we begin to use naturally occurring nanomaterials like Nata de Coco in applications beyond medical diagnostics, we will see that these natural nanomaterials are not just an alternative; rather, they may well be the preferred alternative.

Modular rate laws for enzymatic reactions

MOTIVATION Standard rate laws are a key requisite for systematically turning metabolic networks into kinetic models. They should provide simple, general and biochemically plausible formulae for reaction velocities and reaction elasticities. At the same time, they need to respect thermodynamic relations between the kinetic constants and the metabolic fluxes and concentrations. RESULTS We present a family of reversible rate laws for reactions with arbitrary stoichiometries and various types of regulation, including mass-action, Michaelis-Menten and uni-uni reversible Hill kinetics as special cases. With a thermodynamically safe parameterization of these rate laws, parameter sets obtained by model fitting, sampling or optimization are guaranteed to lead to consistent chemical equilibrium states. A reformulation using saturation values yields simple formulae for rates and elasticities, which can be easily adjusted to the given stationary flux distributions. Furthermore, this formulation highlights the role of chemical potential differences as thermodynamic driving forces. We compare the modular rate laws to the thermodynamic-kinetic modelling formalism and discuss a simplified rate law in which the reaction rate directly depends on the reaction affinity. For automatic handling of modular rate laws, we propose a standard syntax and semantic annotations for the Systems Biology Markup Language. AVAILABILITY An online tool for inserting the rate laws into SBML models is freely available at www.semanticsbml.org. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

Annotation and merging of SBML models with semanticSBML

SUMMARY Systems Biology Markup Language (SBML) is the leading exchange format for mathematical models in Systems Biology. Semantic annotations link model elements with external knowledge via unique database identifiers and ontology terms, enabling software to check and process models by their biochemical meaning. Such information is essential for model merging, one of the key steps towards the construction of large kinetic models. SemanticSBML is a tool that helps users to check and edit MIRIAM annotations and SBO terms in SBML models. Using a large collection of biochemical names and database identifiers, it supports modellers in finding the right annotations and in merging existing models. Initially, an element matching is derived from the MIRIAM annotations and conflicting element attributes are categorized and highlighted. Conflicts can then be resolved automatically or manually, allowing the user to control the merging process in detail. AVAILABILITY SemanticSBML comes as a free software written in Python and released under the GPL 3. A Debian package, a source package for other Linux distributions, a Windows installer and an online version of semanticSBML with limited functionality are available at http://www.semanticsbml.org. A preinstalled version can be found on the Linux live DVD SB.OS, available at http://www.sbos.eu. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.