Herbert M. Sauro

The systems biology markup language (SBML) : a medium for representation and exchange of biochemical network models

MOTIVATION Molecular biotechnology now makes it possible to build elaborate systems models, but the systems biology community needs information standards if models are to be shared, evaluated and developed cooperatively. RESULTS We summarize the Systems Biology Markup Language (SBML) Level 1, a free, open, XML-based format for representing biochemical reaction networks. SBML is a software-independent language for describing models common to research in many areas of computational biology, including cell signaling pathways, metabolic pathways, gene regulation, and others. AVAILABILITY The specification of SBML Level 1 is freely available from http://www.sbml.org/

BioModels Database: a free, centralized database of curated, published, quantitative kinetic models of biochemical and cellular systems

BioModels Database (), part of the international initiative BioModels.net, provides access to published, peer-reviewed, quantitative models of biochemical and cellular systems. Each model is carefully curated to verify that it corresponds to the reference publication and gives the proper numerical results. Curators also annotate the components of the models with terms from controlled vocabularies and links to other relevant data resources. This allows the users to search accurately for the models they need. The models can currently be retrieved in the SBML format, and import/export facilities are being developed to extend the spectrum of formats supported by the resource.

The Systems Biology Graphical Notation

Circuit diagrams and Unified Modeling Language diagrams are just two examples of standard visual languages that help accelerate work by promoting regularity, removing ambiguity and enabling software tool support for communication of complex information. Ironically, despite having one of the highest ratios of graphical to textual information, biology still lacks standard graphical notations. The recent deluge of biological knowledge makes addressing this deficit a pressing concern. Toward this goal, we present the Systems Biology Graphical Notation (SBGN), a visual language developed by a community of biochemists, modelers and computer scientists. SBGN consists of three complementary languages: process diagram, entity relationship diagram and activity flow diagram. Together they enable scientists to represent networks of biochemical interactions in a standard, unambiguous way. We believe that SBGN will foster efficient and accurate representation, visualization, storage, exchange and reuse of information on all kinds of biological knowledge, from gene regulation, to metabolism, to cellular signaling.

Transcriptional dynamics of the embryonic stem cell switch

Recent ChIP experiments of human and mouse embryonic stem cells have elucidated the architecture of the transcriptional regulatory circuitry responsible for cell determination, which involves the transcription factors OCT4, SOX2, and NANOG. In addition to regulating each other through feedback loops, these genes also regulate downstream target genes involved in the maintenance and differentiation of embryonic stem cells. A search for the OCT4–SOX2–NANOG network motif in other species reveals that it is unique to mammals. With a kinetic modeling approach, we ascribe function to the observed OCT4–SOX2–NANOG network by making plausible assumptions about the interactions between the transcription factors at the gene promoter binding sites and RNA polymerase (RNAP), at each of the three genes as well as at the target genes. We identify a bistable switch in the network, which arises due to several positive feedback loops, and is switched on/off by input environmental signals. The switch stabilizes the expression levels of the three genes, and through their regulatory roles on the downstream target genes, leads to a binary decision: when OCT4, SOX2, and NANOG are expressed and the switch is on, the self-renewal genes are on and the differentiation genes are off. The opposite holds when the switch is off. The model is extremely robust to parameter changes. In addition to providing a self-consistent picture of the transcriptional circuit, the model generates several predictions. Increasing the binding strength of NANOG to OCT4 and SOX2, or increasing its basal transcriptional rate, leads to an irreversible bistable switch: the switch remains on even when the activating signal is removed. Hence, the stem cell can be manipulated to be self-renewing without the requirement of input signals. We also suggest tests that could discriminate between a variety of feedforward regulation architectures of the target genes by OCT4, SOX2, and NANOG.

Sensitivity analysis of stoichiometric networks: an extension of metabolic control analysis to non-steady state trajectories

A sensitivity analysis of general stoichiometric networks is considered. The results are presented as a generalization of Metabolic Control Analysis, which has been concerned primarily with system sensitivities at steady state. An expression for time-varying sensitivity coefficients is given and the Summation and Connectivity Theorems are generalized. The results are compared to previous treatments. The analysis is accompanied by a discussion of the computation of the sensitivity coefficients and an application to a model of phototransduction.

The Synthetic Biology Open Language (SBOL) provides a community standard for communicating designs in synthetic biology

The re-use of previously validated designs is critical to the evolution of synthetic biology from a research discipline to an engineering practice. Here we describe the Synthetic Biology Open Language (SBOL), a proposed data standard for exchanging designs within the synthetic biology community. SBOL represents synthetic biology designs in a community-driven, formalized format for exchange between software tools, research groups and commercial service providers. The SBOL Developers Group has implemented SBOL as an XML/RDF serialization and provides software libraries and specification documentation to help developers implement SBOL in their own software. We describe early successes, including a demonstration of the utility of SBOL for information exchange between several different software tools and repositories from both academic and industrial partners. As a community-driven standard, SBOL will be updated as synthetic biology evolves to provide specific capabilities for different aspects of the synthetic biology workflow.

SCAMP : a general - purpose simulator and metabolic control analysis program

SCAMP is a general-purpose simulator of metabolic and chemical networks. The program is written in C and is portable to all computer systems that support an ANSI C compiler. SCAMP accepts metabolic models described in a biochemical language, and this enables novice as well as experienced users rapidly to build and simulate metabolic systems. The language is sufficiently flexible to enable other types of model to be built, e.g. chemostat or ecological models. The language offers many facilities, including: the ability to describe metabolic pathways of any structure and possessing any kinetics using normal chemical notation; optionally build models directly from the differential equations; differing compartment volumes; access to flux, concentration and rate of change information; detection of conserved cycles; access to all coefficients and elasticities of metabolic control analysis; user-defined forcing functions at the model boundaries; user-defined monitoring functions; user-configurable output of any quantity. From the model description SCAMP can either generate C code for later compilation to produce fast executable stand-alone models or run-time code for input to a run-time interpreter for immediate execution. The simulator also incorporates an inbuilt symbolic differentiator for evaluating the Jacobian and elasticity matrices.

In-Fusion BioBrick assembly and re-engineering

Genetic circuits can be assembled from standardized biological parts called BioBricks. Examples of BioBricks include promoters, ribosome-binding sites, coding sequences and transcriptional terminators. Standard BioBrick assembly normally involves restriction enzyme digestion and ligation of two BioBricks at a time. The method described here is an alternative assembly strategy that allows for two or more PCR-amplified BioBricks to be quickly assembled and re-engineered using the Clontech In-Fusion PCR Cloning Kit. This method allows for a large number of parallel assemblies to be performed and is a flexible way to mix and match BioBricks. In-Fusion assembly can be semi-standardized by the use of simple primer design rules that minimize the time involved in planning assembly reactions. We describe the success rate and mutation rate of In-Fusion assembled genetic circuits using various homology and primer lengths. We also demonstrate the success and flexibility of this method with six specific examples of BioBrick assembly and re-engineering. These examples include assembly of two basic parts, part swapping, a deletion, an insertion, and three-way In-Fusion assemblies.

Designing and engineering evolutionary robust genetic circuits

BackgroundOne problem with engineered genetic circuits in synthetic microbes is their stability over evolutionary time in the absence of selective pressure. Since design of a selective environment for maintaining function of a circuit will be unique to every circuit, general design principles are needed for engineering evolutionary robust circuits that permit the long-term study or applied use of synthetic circuits.ResultsWe first measured the stability of two BioBrick-assembled genetic circuits propagated in Escherichia coli over multiple generations and the mutations that caused their loss-of-function. The first circuit, T9002, loses function in less than 20 generations and the mutation that repeatedly causes its loss-of-function is a deletion between two homologous transcriptional terminators. To measure the effect between transcriptional terminator homology levels and evolutionary stability, we re-engineered six versions of T9002 with a different transcriptional terminator at the end of the circuit. When there is no homology between terminators, the evolutionary half-life of this circuit is significantly improved over 2-fold and is independent of the expression level. Removing homology between terminators and decreasing expression level 4-fold increases the evolutionary half-life over 17-fold. The second circuit, I7101, loses function in less than 50 generations due to a deletion between repeated operator sequences in the promoter. This circuit was re-engineered with different promoters from a promoter library and using a kanamycin resistance gene (kanR) within the circuit to put a selective pressure on the promoter. The evolutionary stability dynamics and loss-of-function mutations in all these circuits are described. We also found that on average, evolutionary half-life exponentially decreases with increasing expression levels.ConclusionsA wide variety of loss-of-function mutations are observed in BioBrick-assembled genetic circuits including point mutations, small insertions and deletions, large deletions, and insertion sequence (IS) element insertions that often occur in the scar sequence between parts. Promoter mutations are selected for more than any other biological part. Genetic circuits can be re-engineered to be more evolutionary robust with a few simple design principles: high expression of genetic circuits comes with the cost of low evolutionary stability, avoid repeated sequences, and the use of inducible promoters increases stability. Inclusion of an antibiotic resistance gene within the circuit does not ensure evolutionary stability.

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.