SBML Level 3: an extensible format for the exchange and reuse of biological models.

Systems biology has experienced dramatic growth in the number, size, and complexity of computational models. To reproduce simulation results and reuse models, researchers must exchange unambiguous model descriptions. We review the latest edition of the Systems Biology Markup Language (SBML), a format designed for this purpose. A community of modelers and software authors developed SBML Level 3 over the past decade. Its modular form consists of a core suited to representing reaction-based models and packages that extend the core with features suited to other model types including constraint-based models, reaction-diffusion models, logical network models, and rule-based models. The format leverages two decades of SBML and a rich software ecosystem that transformed how systems biologists build and interact with models. More recently, the rise of multiscale models of whole cells and organs, and new data sources such as single-cell measurements and live imaging, has precipitated new ways of integrating data with models. We provide our perspectives on the challenges presented by these developments and how SBML Level 3 provides the foundation needed to support this evolution.

BioModels: expanding horizons to include more modelling approaches and formats.

BioModels serves as a central repository of mathematical models representing biological processes. It offers a platform to make mathematical models easily shareable across the systems modelling community, thereby supporting model reuse. To facilitate hosting a broader range of model formats derived from diverse modelling approaches and tools, a new infrastructure for BioModels has been developed that is available at http://www.ebi.ac.uk/biomodels. This new system allows submitting and sharing of a wide range of models with improved support for formats other than SBML. It also offers a version-control backed environment in which authors and curators can work collaboratively to curate models. This article summarises the features available in the current system and discusses the potential benefit they offer to the users over the previous system. In summary, the new portal broadens the scope of models accepted in BioModels and supports collaborative model curation which is crucial for model reproducibility and sharing.

ProbOnto: ontology and knowledge base of probability distributions.

MOTIVATION: Probability distributions play a central role in mathematical and statistical modelling. The encoding, annotation and exchange of such models could be greatly simplified by a resource providing a common reference for the definition of probability distributions. Although some resources exist, no suitably detailed and complex ontology exists nor any database allowing programmatic access. RESULTS: ProbOnto, is an ontology-based knowledge base of probability distributions, featuring more than 80 uni- and multivariate distributions with their defining functions, characteristics, relationships and re-parameterization formulas. It can be used for model annotation and facilitates the encoding of distribution-based models, related functions and quantities. AVAILABILITY AND IMPLEMENTATION: http://probonto.org CONTACT: mjswat@ebi.ac.uk SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Filopodial-Tension Model of Convergent-Extension of Tissues.

In convergent-extension (CE), a planar-polarized epithelial tissue elongates (extends) in-plane in one direction while shortening (converging) in the perpendicular in-plane direction, with the cells both elongating and intercalating along the converging axis. CE occurs during the development of most multicellular organisms. Current CE models assume cell or tissue asymmetry, but neglect the preferential filopodial activity along the convergent axis observed in many tissues. We propose a cell-based CE model based on asymmetric filopodial tension forces between cells and investigate how cell-level filopodial interactions drive tissue-level CE. The final tissue geometry depends on the balance between external rounding forces and cell-intercalation traction. Filopodial-tension CE is robust to relatively high levels of planar cell polarity misalignment and to the presence of non-active cells. Addition of a simple mechanical feedback between cells fully rescues and even improves CE of tissues with high levels of polarity misalignments. Our model extends easily to three dimensions, with either one converging and two extending axes, or two converging and one extending axes, producing distinct tissue morphologies, as observed in vivo.

Systems oncology: towards patient-specific treatment regimes informed by multiscale mathematical modelling.

The multiscale complexity of cancer as a disease necessitates a corresponding multiscale modelling approach to produce truly predictive mathematical models capable of improving existing treatment protocols. To capture all the dynamics of solid tumour growth and its progression, mathematical modellers need to couple biological processes occurring at various spatial and temporal scales (from genes to tissues). Because effectiveness of cancer therapy is considerably affected by intracellular and extracellular heterogeneities as well as by the dynamical changes in the tissue microenvironment, any model attempt to optimise existing protocols must consider these factors ultimately leading to improved multimodal treatment regimes. By improving existing and building new mathematical models of cancer, modellers can play important role in preventing the use of potentially sub-optimal treatment combinations. In this paper, we analyse a multiscale computational mathematical model for cancer growth and spread, incorporating the multiple effects of radiation therapy and chemotherapy in the patient survival probability and implement the model using two different cell based modelling techniques. We show that the insights provided by such multiscale modelling approaches can ultimately help in designing optimal patient-specific multi-modality treatment protocols that may increase patients quality of life.

libRoadRunner: a high performance SBML simulation and analysis library.

MOTIVATION: This article presents libRoadRunner, an extensible, high-performance, cross-platform, open-source software library for the simulation and analysis of models expressed using Systems Biology Markup Language (SBML). SBML is the most widely used standard for representing dynamic networks, especially biochemical networks. libRoadRunner is fast enough to support large-scale problems such as tissue models, studies that require large numbers of repeated runs and interactive simulations. RESULTS: libRoadRunner is a self-contained library, able to run both as a component inside other tools via its C++ and C bindings, and interactively through its Python interface. Its Python Application Programming Interface (API) is similar to the APIs of MATLAB ( WWWMATHWORKSCOM: ) and SciPy ( HTTP//WWWSCIPYORG/: ), making it fast and easy to learn. libRoadRunner uses a custom Just-In-Time (JIT) compiler built on the widely used LLVM JIT compiler framework. It compiles SBML-specified models directly into native machine code for a variety of processors, making it appropriate for solving extremely large models or repeated runs. libRoadRunner is flexible, supporting the bulk of the SBML specification (except for delay and non-linear algebraic equations) including several SBML extensions (composition and distributions). It offers multiple deterministic and stochastic integrators, as well as tools for steady-state analysis, stability analysis and structural analysis of the stoichiometric matrix. AVAILABILITY AND IMPLEMENTATION: libRoadRunner binary distributions are available for Mac OS X, Linux and Windows. The library is licensed under Apache License Version 2.0. libRoadRunner is also available for ARM-based computers such as the Raspberry Pi. http://www.libroadrunner.org provides online documentation, full build instructions, binaries and a git source repository. CONTACTS: hsauro@u.washington.edu or somogyie@indiana.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Bystander effects and their implications for clinical radiation therapy: Insights from multiscale in silico experiments.

Radiotherapy is a commonly used treatment for cancer and is usually given in varying doses. At low radiation doses relatively few cells die as a direct response to radiation but secondary radiation effects, such as DNA mutation or bystander phenomena, may affect many cells. Consequently it is at low radiation levels where an understanding of bystander effects is essential in designing novel therapies with superior clinical outcomes. In this paper, we use a hybrid multiscale mathematical model to study the direct effects of radiation as well as radiation-induced bystander effects on both tumour cells and normal cells. We show that bystander responses play a major role in mediating radiation damage to cells at low-doses of radiotherapy, doing more damage than that due to direct radiation. The survival curves derived from our computational simulations showed an area of hyper-radiosensitivity at low-doses that are not obtained using a traditional radiobiological model.

The cell behavior ontology: describing the intrinsic biological behaviors of real and model cells seen as active agents.

MOTIVATION: Currently, there are no ontologies capable of describing both the spatial organization of groups of cells and the behaviors of those cells. The lack of a formalized method for describing the spatiality and intrinsic biological behaviors of cells makes it difficult to adequately describe cells, tissues and organs as spatial objects in living tissues, in vitro assays and in computational models of tissues. RESULTS: We have developed an OWL-2 ontology to describe the intrinsic physical and biological characteristics of cells and tissues. The Cell Behavior Ontology (CBO) provides a basis for describing the spatial and observable behaviors of cells and extracellular components suitable for describing in vivo, in vitro and in silico multicell systems. Using the CBO, a modeler can create a meta-model of a simulation of a biological model and link that meta-model to experiment or simulation results. Annotation of a multicell model and its computational representation, using the CBO, makes the statement of the underlying biology explicit. The formal representation of such biological abstraction facilitates the validation, falsification, discovery, sharing and reuse of both models and experimental data. AVAILABILITY AND IMPLEMENTATION: The CBO, developed using Protégé 4, is available at http://cbo.biocomplexity.indiana.edu/cbo/ and at BioPortal (http://bioportal.bioontology.org/ontologies/CBO).

Effects of a hypercaloric diet on ß-cell responsivity in lean healthy men.

OBJECTIVE: Insulin resistance and hyperinsulinaemia precede the onset of obesity-induced DM2. The early adaptation of the ß-cell during the initial phase of overfeeding and weight gain has only been partly elucidated. We studied the early changes in insulin clearance and ß-cell responsivity during a positive and negative energy balance in lean healthy men. DESIGN: We studied in nine healthy lean men [age, 37 (27-43) years; BMI, 23·6 (20·6-25·6) kg/m(2) ] insulin sensitivity, insulin clearance, insulin secretion and static and dynamic ß-cell responsivity at baseline and after the hypercaloric and subsequent hypocaloric diet. RESULTS: Participants gained 7 [5·1-7·6]% of their initial body weight on the hypercaloric diet. Compared to baseline, insulin sensitivity and insulin clearance decreased, while glucose-stimulated insulin secretion was higher. The GLP-1 response to oral glucose did not change. The dynamic ß-cell responsivity index increased but the basal and static responsivity indexes did not change. Total and static disposition indexes (DIs) in the hypercaloric state showed a trend towards a decrease. During the hypocaloric diet, insulin sensitivity, glucose-stimulated insulin secretion and insulin clearance returned to baseline. The responsivity and the DIs were not different in the hypocaloric phase compared to baseline. CONCLUSION: A positive energy balance resulting in weight gain in lean men induces hyperinsulinaemia, which is explained by a combined effect on insulin clearance and insulin secretion. Increased insulin secretion was related to insulin resistance-induced higher glucose concentrations but also to increased dynamic ß-cell responsivity. Glucose sensitivity of the ß-cell did not change. These early adaptations are completely reversible during a negative energy balance after loss of the gained weight.

Adhesion failures determine the pattern of choroidal neovascularization in the eye: a computer simulation study.

Choroidal neovascularization (CNV) of the macular area of the retina is the major cause of severe vision loss in adults. In CNV, after choriocapillaries initially penetrate Bruch's membrane (BrM), invading vessels may regress or expand (CNV initiation). Next, during Early and Late CNV, the expanding vasculature usually spreads in one of three distinct patterns: in a layer between BrM and the retinal pigment epithelium (sub-RPE or Type 1 CNV), in a layer between the RPE and the photoreceptors (sub-retinal or Type 2 CNV) or in both loci simultaneously (combined pattern or Type 3 CNV). While most studies hypothesize that CNV primarily results from growth-factor effects or holes in BrM, our three-dimensional simulations of multi-cell model of the normal and pathological maculae recapitulate the three growth patterns, under the hypothesis that CNV results from combinations of impairment of: 1) RPE-RPE epithelial junctional adhesion, 2) Adhesion of the RPE basement membrane complex to BrM (RPE-BrM adhesion), and 3) Adhesion of the RPE to the photoreceptor outer segments (RPE-POS adhesion). Our key findings are that when an endothelial tip cell penetrates BrM: 1) RPE with normal epithelial junctions, basal attachment to BrM and apical attachment to POS resists CNV. 2) Small holes in BrM do not, by themselves, initiate CNV. 3) RPE with normal epithelial junctions and normal apical RPE-POS adhesion, but weak adhesion to BrM (e.g. due to lipid accumulation in BrM) results in Early sub-RPE CNV. 4) Normal adhesion of RBaM to BrM, but reduced apical RPE-POS or epithelial RPE-RPE adhesion (e.g. due to inflammation) results in Early sub-retinal CNV. 5) Simultaneous reduction in RPE-RPE epithelial binding and RPE-BrM adhesion results in either sub-RPE or sub-retinal CNV which often progresses to combined pattern CNV. These findings suggest that defects in adhesion dominate CNV initiation and progression.

Learn-confirm in model-informed drug development: Assessing an immunogenicity quantitative systems pharmacology platform.

Immunogenicity against therapeutic proteins frequently causes attrition owing to its potential impact on pharmacokinetics, pharmacodynamics, efficacy, and safety. Predicting immunogenicity is complex because of its multifactorial drivers, including compound properties, subject characteristics, and treatment parameters. To integrate these, the Immunogenicity Simulator was developed using published, predominantly late-stage trial data from 15 therapeutic proteins. This single-blinded evaluation with subject-level data from 10 further monoclonals assesses the Immunogenicity Simulator's credibility for application during the drug development process.

Computer simulation of cellular patterning within the Drosophila pupal eye.

We present a computer simulation and associated experimental validation of assembly of glial-like support cells into the interweaving hexagonal lattice that spans the Drosophila pupal eye. This process of cell movements organizes the ommatidial array into a functional pattern. Unlike earlier simulations that focused on the arrangements of cells within individual ommatidia, here we examine the local movements that lead to large-scale organization of the emerging eye field. Simulations based on our experimental observations of cell adhesion, cell death, and cell movement successfully patterned a tracing of an emerging wild-type pupal eye. Surprisingly, altering cell adhesion had only a mild effect on patterning, contradicting our previous hypothesis that the patterning was primarily the result of preferential adhesion between IRM-class surface proteins. Instead, our simulations highlighted the importance of programmed cell death (PCD) as well as a previously unappreciated variable: the expansion of cells' apical surface areas, which promoted rearrangement of neighboring cells. We tested this prediction experimentally by preventing expansion in the apical area of individual cells: patterning was disrupted in a manner predicted by our simulations. Our work demonstrates the value of combining computer simulation with in vivo experiments to uncover novel mechanisms that are perpetuated throughout the eye field. It also demonstrates the utility of the Glazier-Graner-Hogeweg model (GGH) for modeling the links between local cellular interactions and emergent properties of developing epithelia as well as predicting unanticipated results in vivo.

Development of a coupled simulation toolkit for computational radiation biology based on Geant4 and CompuCell3D.

Understanding and designing clinical radiation therapy is one of the most important areas of state-of-the-art oncological treatment regimens. Decades of research have gone into developing sophisticated treatment devices and optimization protocols for schedules and dosages. In this paper, we presented a comprehensive computational platform that facilitates building of the sophisticated multi-cell-based model of how radiation affects the biology of living tissue. We designed and implemented a coupled simulation method, including a radiation transport model, and a cell biology model, to simulate the tumor response after irradiation. The radiation transport simulation was implemented through Geant4 which is an open-source Monte Carlo simulation platform that provides many flexibilities for users, as well as low energy DNA damage simulation physics, Geant4-DNA. The cell biology simulation was implemented using CompuCell3D (CC3D) which is a cell biology simulation platform. In order to couple Geant4 solver with CC3D, we developed a 'bridging' module, RADCELL, that extracts tumor cellular geometry of the CC3D simulation (including specification of the individual cells) and ported it to the Geant4 for radiation transport simulation. The cell dose and cell DNA damage distribution in multicellular system were obtained using Geant4. The tumor response was simulated using cell-based tissue models based on CC3D, and the cell dose and cell DNA damage information were fed back through RADCELL to CC3D for updating the cell properties. By merging two powerful and widely used modeling platforms, CC3D and Geant4, we delivered a novel tool that can give us the ability to simulate the dynamics of biological tissue in the presence of ionizing radiation, which provides a framework for quantifying the biological consequences of radiation therapy. In this introductory methods paper, we described our modeling platform in detail and showed how it can be applied to study the application of radiotherapy to a vascularized tumor.

Cycle Network Model of Prostaglandin H Synthase-1.

The kinetic model of Prostaglandin H Synthase-1 (PGHS-1) was developed to investigate its complex network kinetics and non-steroidal anti-inflammatory drugs (NSAIDs) efficacy in different in vitro and in vivo conditions. To correctly describe the complex mechanism of PGHS-1 catalysis, we developed a microscopic approach to modelling of intricate network dynamics of 35 intraenzyme reactions among 24 intermediate states of the enzyme. The developed model quantitatively describes interconnection between cyclooxygenase and peroxidase enzyme activities; substrate (arachidonic acid, AA) and reducing cosubstrate competitive consumption; enzyme self-inactivation; autocatalytic role of AA; enzyme activation threshold; and synthesis of intermediate prostaglandin G2 (PGG2) and final prostaglandin H2 (PGH2) products under wide experimental conditions. In the paper, we provide a detailed description of the enzyme catalytic cycle, model calibration based on a series of in vitro kinetic data, and model validation using experimental data on the regulatory properties of PGHS-1. The validated model of PGHS-1 with a unified set of kinetic parameters is applicable for in silico screening and prediction of the inhibition effects of NSAIDs and their combination on the balance of pro-thrombotic (thromboxane) and anti-thrombotic (prostacyclin) prostaglandin biosynthesis in platelets and endothelial cells expressing PGHS-1.

Multicell simulations of development and disease using the CompuCell3D simulation environment.

Mathematical modeling and computer simulation have become crucial to biological fields from genomics to ecology. However, multicell, tissue-level simulations of development and disease have lagged behind other areas because they are mathematically more complex and lack easy-to-use software tools that allow building and running in silico experiments without requiring in-depth knowledge of programming. This tutorial introduces Glazier-Graner-Hogeweg (GGH) multicell simulations and CompuCell3D, a simulation framework that allows users to build, test, and run GGH simulations.

Front instabilities and invasiveness of simulated avascular tumors.

We study the interface morphology of a 2D simulation of an avascular tumor composed of identical cells growing in an homogeneous healthy tissue matrix (TM), in order to understand the origin of the morphological changes often observed during real tumor growth. We use the Glazier-Graner-Hogeweg model, which treats tumor cells as extended, deformable objects, to study the effects of two parameters: a dimensionless diffusion-limitation parameter defined as the ratio of the tumor consumption rate to the substrate transport rate, and the tumor-TM surface tension. We model TM as a nondiffusing field, neglecting the TM pressure and haptotactic repulsion acting on a real growing tumor; thus, our model is appropriate for studying tumors with highly motile cells, e.g., gliomas. We show that the diffusion-limitation parameter determines whether the growing tumor develops a smooth (noninvasive) or fingered (invasive) interface, and that the sensitivity of tumor morphology to tumor-TM surface tension increases with the size of the dimensionless diffusion-limitation parameter. For large diffusion-limitation parameters, we find a transition (missed in previous work) between dendritic structures, produced when tumor-TM surface tension is high, and seaweed-like structures, produced when tumor-TM surface tension is low. This observation leads to a direct analogy between the mathematics and dynamics of tumors and those observed in nonbiological directional solidification. Our results are also consistent with the biological observation that hypoxia promotes invasive growth of tumor cells by inducing higher levels of receptors for scatter factors that weaken cell-cell adhesion and increase cell motility. These findings suggest that tumor morphology may have value in predicting the efficiency of antiangiogenic therapy in individual patients.

Development of computational model for cell dose and DNA damage quantification of multicellular system.

Purpose: The aim of this study is to build a computational model to investigate the cell dose and cell DNA damage distribution of a multicellular tissue system under the irradiation.Materials and methods: In this work, we developed a computational model for quantifying cell dose and double strand break (DSB) number in a multicellular system by simulating the radiation transport in 2D and 3D cell culture. The model was based on an open-source radiation transport package, Geant4 with Geant4-DNA physics. First, the computational multicellular system was created using a developed program, CelllMaker. Second, the radiation transport simulation for cells was conducted using Geant4 package with the Geant4-DNA physics to obtain the cellular dose and cellular DSB yield.Results: Using the method described in this work, it is possible to obtain the cellular dose and DNA damage simultaneously. The developed model provides a solution for quantifying the cellular dose and cellular DNA damage which are not easily determined in a radiobiological experiment.Conclusions: With limited validation data for the model, this preliminary study provides a roadmap for building a comprehensive toolkit for simulating cellular dose and DNA damage of multicellular tissue systems.

Computational model of wound healing: EGF secreted by fibroblasts promotes delayed re-epithelialization of epithelial keratinocytes.

It is widely agreed that keratinocyte migration plays a crucial role in wound re-epithelialization. Defects in this function contribute to wound reoccurrence causing significant clinical problems. Several in vitro studies have shown that the speed of migrating keratinocytes can be regulated by epidermal growth factor (EGF) which affects keratinocyte's integrin expression. The relationship between integrin expression (through cell-matrix adhesion) stimulated by EGF and keratinocyte migration speed is not linear since increased adhesion, due to increased integrin expression, has been experimentally shown to slow down cell migration due to the biphasic dependence of cell speed on adhesion. In our previous work we showed that keratinocytes that were co-cultured with EGF-enhanced fibroblasts formed an asymmetric migration pattern, where, the cumulative distances of keratinocytes migrating toward fibroblasts were smaller than those migrating away from fibroblasts. This asymmetric pattern is thought to be provoked by high EGF concentration secreted by fibroblasts. The EGF stimulates the expression of integrin receptors on the surface of keratinocytes migrating toward fibroblasts via paracrine signaling. In this paper, we present a computational model of keratinocyte migration that is controlled by EGF secreted by fibroblasts using the Cellular Potts Model (CPM). Our computational simulation results confirm the asymmetric pattern observed in experiments. These results provide a deeper insight into our understanding of the complexity of keratinocyte migration in the presence of growth factor gradients and may explain re-epithelialization failure in impaired wound healing.

Cytoskeletal remodeling in lymphocyte activation.

The formerly distinct fields of lymphocyte signal transduction and cytoskeletal remodeling have recently become linked, as proteins involved in transducing signals downstream of lymphocyte antigen receptors have also been implicated in actin cytoskeleton remodeling, microtubule dynamics and regulation of cell polarity. These discoveries have fuelled interest in understanding both the role of the actin cytoskeleton as an integral component of lymphocyte activation and the interplay between lymphoid cell-cell contact sites (immunological synapse), retractile pole structures (uropod, distal pole complex), and Rho-family GTPases (Rac, Rho, Cdc42), their upstream activators (Dbl-family guanine nucleotide exchange factors) and their downstream effectors (WASp, Arp2/3, ADAP). To understand how these complex regulatory networks are wired, a new breed of computational biologists uses mathematical language to reproduce and simulate signaling circuits 'in silico'.