Systems biology relies heavily on building mathematical models to help understand and make predictions of biological processes. Specialized software to assist in building models has been developed since the arrival of the first digital computers. [1] [2] [3] [4] The following list gives the currently supported software applications available to researchers.
The vast majority of modern systems biology modeling software support SBML, which is the de facto standard for exchanging models of biological cellular processes. Some tools also support CellML, a standard used for representing physiological processes. The advantage of using standard formats is that even though a particular software application may eventually become unsupported and even unusable, the models developed by that application can be easily transferred to more modern equivalents. This allows scientific research to be reproducible long after the original publication of the work.
To obtain more information about a particular tool, click on the name of the tool. This will direct you either to a peer-reviewed publication or, in some rare cases, to a dedicated Wikipedia page.
When an entry in the SBML column states "Yes, but only for reactions.", it means that the tool only supports the reaction component of SBML. For example, rules, events, etc. are not supported.
Name | Description/Notability | OS | License | Site | SBML Support |
---|---|---|---|---|---|
iBioSim | iBioSim [5] [6] is a computer-aided design (CAD) tool for the modeling, analysis, and design of genetic circuits. | multiplatform (Java/C++) | Apache License | Yes | |
CompuCell3D | GUI/Scripting tool [7] for building and simulating multicellular models. | multiplatform (C++/Python) | MIT | Yes, but only for reactions. | |
COPASI | GUI tool [8] [9] for analyzing and simulating SBML models. | multiplatform (C++) | Artistic License | Yes | |
Cytosim | Spatial simulator for flexible cytoskeletal filaments and motor proteins [10] | Mac, Linux, Cygwin (C++) | GPL3 | Not applicable | |
libroadrunner | High-performance software library for simulation and analysis of SBML models [11] [12] | multiplatform (C/C++) | Apache License | Yes | |
massPy | Simulation tool [13] [14] that can work with COBRApy [15] | multiplatform (Python) | MIT | Yes | |
MCell | GUI tool for particle-based spatial stochastic simulation with individual molecules [16] [17] [18] | multiplatform | MIT and GPLv2 | Not applicable | |
OpenCOR | A cross-platform modelling environment, which is aimed at organizing, editing, simulating, and analysing CellML files on Windows, Linux and macOS. | multiplatform (C++/Python) | GPLv3 | Uses CellML | |
PhysiBoSS | A specialized form of the PhysiCell agent-based modeling platform that directly integrates Boolean signaling networks into cell Agents [19] | multiplatform (C++) | BSD-3 | Yes, but only for reactions | |
PhysiCell | A agent-based [20] modeling framework for multicellular systems biology. | multiplatform (C++) | BSD-3 | Yes, but only for reactions | |
PySCeS | Python tool for modeling and analyzing SBML models [21] [22] [23] | multiplatform (Python) | BSD-3 | Yes | |
pySB | Python-based [24] platform with specialization in rule-based models. | multiplatform (Python) | BSD-3 | Partial | |
ReaDDy | Particle-based spatial simulator with intermolecular potentials [25] | Linux and Mac | Custom | Not applicable | |
SBSCL | Java library [26] [27] with efficient and exhaustive support for SBML | multiplatform (Java) | LGPL | Yes | |
SBW (software) | A distributed workbench [28] [29] that includes many modeling tools | multiplatform (C/C++) | BSD-3 | Yes | |
Smoldyn | Particle-based simulator for spatial stochastic simulations with individual molecules [30] [31] [32] [33] | multiplatform (C/C++/Python) | LGPL | Not applicable | |
Spatiocyte | Spatial modeling software that uses a fine lattice with up to one molecule per site [34] [35] | multiplatform | Unknown | Not applicable | |
SpringSaLaD | Particle-based spatial simulator in which molecules are spheres that are linked by springs [36] | multiplatform | Unknown | Not applicable | |
STEPS | Stochastic reaction-diffusion and membrane potential solver on distributed meshes [37] [38] [39] [40] | multiplatform (C++/Python) | GPLv2 | Partial | |
Tellurium (software) | Simulation environment, [41] [42] that packages multiple libraries into one platform. | multiplatform (Python) | Apache License | Yes | |
URDME | Stochastic reaction-diffusion simulation on unstructured meshes [43] | MatLab on Mac, Linux | GPL3 | Not applicable | |
VCell | Comprehensive modeling platform [44] [45] for non-spatial, spatial, deterministic and stochastic simulations, including both reaction networks and reaction rules. | multiplatform (Java) | MIT | Yes | |
The following table lists specialist tools that cannot be grouped with the modeling tools.
Name | Description/Notability | OS | License | Site |
---|---|---|---|---|
PySCeSToolbox | PySCeSToolbox [46] is a set of metabolic model analysis tools. Among other features, it can be used to generate the control analysis equations that relate the elasticities to the control coefficients. The package is cross-platform and requires PySCeS and Maxima to operate. | multiplatform (C++/Python) | BSD-3 | |
Name | ODE | Constraint based | Stochastic | Logical | Agent based | Spatial (particle) | Spatial (continuous) |
---|---|---|---|---|---|---|---|
iBioSim | Yes | No | Yes | No | Limited | No | No |
CompuCell3D | Yes | No | No | No | Yes | No | Yes |
COPASI | Yes | No | Yes | No | No | No | No |
Cytosim | No | No | Yes | No | ? | Yes | ? |
libroadrunner | Yes | No | Yes | No | No | No | No |
massPy | Uses libroadrunner | Uses COBRApy | No | No | No | No | |
MCell | No | No | Yes | No | No | Yes | No |
OpenCOR | Yes | No | No | No | No | No | No |
PhysiBoSS | |||||||
PhysiCell | Uses libroadrunner | No | No | No | Yes | ? | Yes |
PySCeS | Yes | No | ? | No | No | No | No |
pySB | Yes | No | No | No | No | No | No |
ReaDDy | |||||||
SBSCL | Yes | ? | ? | No | No | No | No |
SBW | Yes | No | Yes | No | No | No | No |
Smoldyn | No | No | Yes | No | No | Yes | No |
Spatiocyte | |||||||
SpringSaLaD | |||||||
STEPS | |||||||
Tellurium (software) | Uses libroadrunner | ||||||
URDME | |||||||
VCell | Yes | No | ? | No | No | No | Single Cell |
Name | Non-stiff solver | Stiff solver | Steady-state solver | Steady-state sensitivities | Time-dependent sensitivities | Bifurcation Analysis |
---|---|---|---|---|---|---|
iBioSim | Yes | Yes | No | No | ? | No |
CompuCell3D | Uses libroadrunner | NA | ||||
COPASI | Yes | Yes | Yes | Yes | ? | Limited |
libroadrunner | Yes | Yes | Yes | Yes | Yes | via AUTO2000 plugin |
masspy | Uses libroadrunner | |||||
OpenCOR | Yes | Yes | ? | ? | ? | No |
PhysiBoSS | ||||||
PhysiCell | Uses libroadrunner | |||||
PySCeS | Yes | Yes | Yes | Yes | ? | Limited+ |
pySB | Yes | No | No | No | No | No |
SBSCL | ||||||
SBW | Uses C# edition of roadrunner | Yes | ||||
Tellurium (software) | Uses libroadrunner | |||||
VCell | Yes | Yes | No | No | No | No |
Name | Import | Export | Primary Interface | Network visualization (editing) |
---|---|---|---|---|
iBioSim | SBML | SBML | GUI | Yes (Yes) |
CompuCell3D | Native XML specification format and SBML | Native XML | GUI/Python scripting | No |
COPASI | Native XML specification format and SBML | Native XML and SBML | GUI | Yes (No) |
libroadrunner | SBML | SBML | Python scripting | No |
masspy | SBML | SBML | Python scripting | No |
Name | Stoichiometry matrix | Reduced stoich matrix | Conserved moiety analysis | Jacobian | MCA |
---|---|---|---|---|---|
COPASI | Yes | Yes | Yes | Yes | Yes |
libroadrunner | Yes | Yes | Yes | Yes | Yes |
masspy | via libroadrunner | ||||
PySCeS | Yes | Yes | Yes | Yes | Yes |
VCell | ? | ? | ? | ? | Limited |
Name | Parameter Estimation | DAE support | Units support |
---|---|---|---|
iBioSim | No | ? | ? |
ComputeCell3D | NA | NA | ? |
COPASI | Yes | Limited | Yes |
libroadrunner | via Python packages | Limited | Yes |
masspy | via Python packages | Limited | Yes |
Particle based simulators treat each molecule of interest as an individual particle in continuous space, simulating molecular diffusion, molecule-membrane interactions and chemical reactions. [47]
The following list compares the features for several particle-based simulators. This table is edited from a version that was originally published in the Encyclopedia of Computational Neuroscience. [48] System boundaries codes: R = reflecting, A = absorbing, T = transmitting, P = periodic, and I = interacting. * Algorithm is exact but software produced incorrect results at the time of original table compilation. † These benchmark run times are not comparable with others due to differing levels of detail.
Feature | MCell | Smoldyn | eGFRD | SpringSaLaD | ReaDDy |
---|---|---|---|---|---|
Time steps | ~1 us | ns to ms | event-based | ~10 ns | ~0.1 ns to us |
Molecules | points | points, spheres | spheres | multi-spheres | multi-spheres |
Dimensions | 2,3 | 1,2,3 | 3 | 3 | 3 |
System boundaries | R,A,P,T | R,A,P,T | P | R | P,I |
Surfaces | triangle mesh | many primitives | - | 1 flat surface | plane, sphere |
Surface molecules | 1/tile, 2 states | unlimited, 4 states | - | unlimited, 3 states | - |
Excluded volume | - | excellent | exact | good | excellent |
Multimers | states only | rule-based model | - | explicit | explicit |
Allostery | - | yes | - | yes | - |
Reaction accuracy | very good | excellent | exact* | excellent | excellent |
Dissociation products | stochastic | fixed separation | adjacent | adjacent | adjacent |
Molecule-surface interactions | good | excellent | - | to sites only | potentials |
Long-range interactions | - | yes | - | - | yes |
Benchmark run time | 67 s | 22 s | 13 days† | 9.1 months† | 13 minutes |
Distribution | executable | executable | self-compile | Java file | self-compile |
User interface | GUI, text | text | text | GUI | Python script |
Graphical output | excellent | good | partial support | partial support | good |
Library interface | Python | C/C++, Python | - | - | Python |
References | [52] [53] | [54] [55] [56] | [57] | [58] |
Model calibration is a key activity when developing systems biology models. This table highlights some of the current model calibration tools available to systems biology modelers. The first table list tools that are SBML compatible.
Tool | PEtab Compatible | P1 | P2 |
---|---|---|---|
pyPESTO [59] | Yes | NA | NA |
COPASI | Yes | NA | NA |
PEtab [60] is a community standard for specifying model calibration runs.
The following list some very early software for modeling biochemical systems that were developed pre-1980s There are listed for historical interest.
Name | Description/Notability | Language | Terminus ante quem [61] |
---|---|---|---|
BIOSIM [62] | The first ever recorded digital simulator of biochemical networks (by David Garfinkel) | FORTRAN IV | 1968 |
KDF 9 [63] | First simulator to support MCA. Developed by the late Jim Burns in Edinburgh | Early form of FORTRAN | 1968 |
METASIM [64] | Early simulator by Park and Wright | PL/1 | 1973 |
The following list shows some of the software modeling applications that were developed in the 1980s and 1990s. There are listed for historical interest.
Name | Description/Notability | Language | SBML Support | Terminus ante quem [65] |
---|---|---|---|---|
COR [66] | First public CellML-based environment. | Object Pascal | Uses CellML | 2010 |
DBsolve [67] | Early GUI based simulation platform. | C/C++ | No | 1999 |
E-Cell [68] | One of the earliest attempts at a whole-cell modeling platform. | C/C++ | No | 1999 |
Gepasi [69] | First GUI application that supported metabolic control analysis and parameter estimation. | C/C++ | Yes | 1993 |
Jarnac [70] | First GUI based application to support scripting in systems biology modeling. | Object Pascal | Yes | 2000 |
JSim [71] | First Java-based systems biology modeling platform | Java | Yes | 2003 |
MetaMod [72] | One of the first PC-based systems biology simulators | BBC Micro | No | 1986 |
MetaModel [73] | Early PC-based systems biology simulator | Turbo Pascal 5.0 | No | 1991 |
MIST [74] | GUI based simulator | Borland Pascal 7.0 | No | 1995 |
SCAMP [75] | First application to support metabolic control analysis and simulation on a PC | Pascal, later in C | No | 1985 (Thesis) |
Bioinformatics is an interdisciplinary field of science that develops methods and software tools for understanding biological data, especially when the data sets are large and complex. Bioinformatics uses biology, chemistry, physics, computer science, computer programming, information engineering, mathematics and statistics to analyze and interpret biological data. The subsequent process of analyzing and interpreting data is referred to as computational biology.
Computational biology refers to the use of data analysis, mathematical modeling and computational simulations to understand biological systems and relationships. An intersection of computer science, biology, and big data, the field also has foundations in applied mathematics, chemistry, and genetics. It differs from biological computing, a subfield of computer science and engineering which uses bioengineering to build computers.
A generegulatory network (GRN) is a collection of molecular regulators that interact with each other and with other substances in the cell to govern the gene expression levels of mRNA and proteins which, in turn, determine the function of the cell. GRN also play a central role in morphogenesis, the creation of body structures, which in turn is central to evolutionary developmental biology (evo-devo).
Systems biology is the computational and mathematical analysis and modeling of complex biological systems. It is a biology-based interdisciplinary field of study that focuses on complex interactions within biological systems, using a holistic approach to biological research.
In academia, computational immunology is a field of science that encompasses high-throughput genomic and bioinformatics approaches to immunology. The field's main aim is to convert immunological data into computational problems, solve these problems using mathematical and computational approaches and then convert these results into immunologically meaningful interpretations.
Metabolic network modelling, also known as metabolic network reconstruction or metabolic pathway analysis, allows for an in-depth insight into the molecular mechanisms of a particular organism. In particular, these models correlate the genome with molecular physiology. A reconstruction breaks down metabolic pathways into their respective reactions and enzymes, and analyzes them within the perspective of the entire network. In simplified terms, a reconstruction collects all of the relevant metabolic information of an organism and compiles it in a mathematical model. Validation and analysis of reconstructions can allow identification of key features of metabolism such as growth yield, resource distribution, network robustness, and gene essentiality. This knowledge can then be applied to create novel biotechnology.
Systems immunology is a research field under systems biology that uses mathematical approaches and computational methods to examine the interactions within cellular and molecular networks of the immune system. The immune system has been thoroughly analyzed as regards to its components and function by using a "reductionist" approach, but its overall function can't be easily predicted by studying the characteristics of its isolated components because they strongly rely on the interactions among these numerous constituents. It focuses on in silico experiments rather than in vivo.
Fluxomics describes the various approaches that seek to determine the rates of metabolic reactions within a biological entity. While metabolomics can provide instantaneous information on the metabolites in a biological sample, metabolism is a dynamic process. The significance of fluxomics is that metabolic fluxes determine the cellular phenotype. It has the added advantage of being based on the metabolome which has fewer components than the genome or proteome.
In enzymology, a guanylate kinase is an enzyme that catalyzes the chemical reaction
Anders Krogh is a bioinformatician at the University of Copenhagen, where he leads the university's bioinformatics center. He is known for his pioneering work on the use of hidden Markov models in bioinformatics, and is co-author of a widely used textbook in bioinformatics. In addition, he also co-authored one of the early textbooks on neural networks. His current research interests include promoter analysis, non-coding RNA, gene prediction and protein structure prediction.
The ISCB Overton Prize is a computational biology prize awarded annually for outstanding accomplishment by a scientist in the early to mid stage of his or her career. Laureates have made significant contribution to the field of computational biology either through research, education, service, or a combination of the three.
The ISCB Accomplishment by a Senior Scientist Award is an annual prize awarded by the International Society for Computational Biology for contributions to the field of computational biology.
Trey Ideker is a professor of medicine and bioengineering at UC San Diego. He is the Director of the National Resource for Network Biology, the San Diego Center for Systems Biology, and the Cancer Cell Map Initiative. He uses genome-scale measurements to construct network models of cellular processes and disease.
Multi-state modeling of biomolecules refers to a series of techniques used to represent and compute the behaviour of biological molecules or complexes that can adopt a large number of possible functional states.
Curtis Huttenhower is a Professor of Computational Biology and Bioinformatics in the Department of Biostatistics, School of Public Health, Harvard University.
Nicolas Le Novère is a British and French biologist. His research focuses on modeling signaling pathways and developing tools to share mathematical models.
Cell-based models are mathematical models that represent biological cells as discrete entities. Within the field of computational biology they are often simply called agent-based models of which they are a specific application and they are used for simulating the biomechanics of multicellular structures such as tissues. to study the influence of these behaviors on how tissues are organised in time and space. Their main advantage is the easy integration of cell level processes such as cell division, intracellular processes and single-cell variability within a cell population.
C-ImmSim started, in 1995, as the C-language "version" of IMMSIM, the IMMune system SIMulator, a program written back in 1991 in APL-2 by the astrophysicist Phil E. Seiden together with the immunologist Franco Celada to implement the Celada-Seiden model. The porting was mainly conducted and further developed by Filippo Castiglione with the help of few other people.
libRoadRunner is a C/C++ software library that supports simulation of SBML based models.. It uses LLVM to generate extremely high-performance code and is the fastest SBML-based simulator currently available. Its main purpose is for use as a reusable library that can be hosted by other applications, particularly on large compute clusters for doing parameter optimization where performance is critical. It also has a set of Python bindings that allow it to be easily used from Python.
Smoldyn is an open-source software application for cell-scale biochemical simulations. It uses particle-based simulation, meaning that it simulates each molecule of interest individually, in order to capture natural stochasticity and yield nanometer-scale spatial resolution. Simulated molecules diffuse, react, are confined by surfaces, and bind to membranes in similar manners as in real biochemical systems.
{{cite journal}}
: Cite journal requires |journal=
(help)