This package provides importers to load reaction networks into
Catalyst.jl
ReactionSystems
from several file formats. Currently it supports loading networks in the
following formats:
- A subset of the BioNetGen .net file format.
- Networks represented by dense or sparse substrate and product stoichiometric matrices.
- Networks represented by dense or sparse complex stoichiometric and incidence matrices.
SBMLToolkit.jl provides an alternative for loading SBML files into Catalyst models, offering a much broader set of supported features. It allows the import of models that include features such as constant species, boundary condition species, events, constraint equations and more. SBML files can be generated from many standard modeling tools, including BioNetGen, COPASI, and Virtual Cell.
For information on using the package, see the stable documentation. Use the in-development documentation for the version of the documentation which contains the unreleased features.
A simple network from the builtin BioNetGen bngl examples is the
repressilator. The generate_network
command in the bngl file outputs a reduced network description, i.e. a
.net file, which can be loaded into a
Catalyst ReactionSystem as:
using ReactionNetworkImporters
fname = "PATH/TO/Repressilator.net"
prnbng = loadrxnetwork(BNGNetwork(), fname)Here BNGNetwork is a type specifying the file format that is being loaded.
prnbng is a ParsedReactionNetwork structure with the following fields:
rn, a CatalystReactionSystemu0, aDictmapping initial condition symbolic variables to numeric values and/or symbolic expressions.p, aDictmapping parameter symbolic variables to numeric values and/or symbolic expressions.varstonames, aDictmapping the internal symbolic variable of a species used in the generatedReactionSystemto aStringgenerated from the name in the .net file. This is necessary as BioNetGen can generate exceptionally long species names, involving characters that lead to malformed species names when used withCatalyst.groupstosyms, aDictmapping theStrings representing names for any groups defined in the BioNetGen file to the corresponding symbolic variable representing theModelingToolkitsymbolic observable associated with the group.
Given prnbng, we can construct and solve the corresponding ODE model for the
reaction system by
using OrdinaryDiffEq, Catalyst
rn = prnbng.rn
tf = 100000.0
oprob = ODEProblem(rn, Float64[], (0.0, tf), Float64[])
sol = solve(oprob, Tsit5(), saveat = tf / 1000.0)Note that we specify empty parameter and initial condition vectors as these are
already stored in the generated ReactionSystem, rn. A Dict mapping each
symbolic species and parameter to its initial value or symbolic expression can
be obtained using ModelingToolkit.defaults(rn).
See the Catalyst documentation for how to generate ODE, SDE, jump and other types of models.
Catalyst ReactionSystems can also be constructed from
- substrate and product stoichiometric matrices.
- complex stoichiometric and incidence matrices.
For example, here we both directly build a Catalyst
network using the @reaction_network macro, and then show how to build the same
network from these matrices using ReactionNetworkImporters:
# Catalyst network from the macro:
rs = @reaction_network testnetwork begin
k1, 2A --> B
k2, B --> 2A
k3, A + B --> C
k4, C --> A + B
k5, 3C --> 3A
end
# network from basic stoichiometry using ReactionNetworkImporters
@parameters k1 k2 k3 k4 k5
@variables t
@species A(t) B(t) C(t)
species = [A, B, C]
pars = [k1, k2, k3, k4, k5]
substoich = [2 0 1 0 0;
0 1 1 0 0;
0 0 0 1 3]
prodstoich = [0 2 0 1 3;
1 0 0 1 0;
0 0 1 0 0]
mn = MatrixNetwork(pars, substoich, prodstoich; species = species,
params = pars) # a matrix network
prn = loadrxnetwork(mn; name = :testnetwork) # dense version
# test the two networks are the same
@assert rs == prn.rn
# network from reaction complex stoichiometry
stoichmat = [2 0 1 0 0 3;
0 1 1 0 0 0;
0 0 0 1 3 0]
incidencemat = [-1 1 0 0 0;
1 -1 0 0 0;
0 0 -1 1 0;
0 0 1 -1 0;
0 0 0 0 -1;
0 0 0 0 1]
cmn = ComplexMatrixNetwork(pars, stoichmat, incidencemat; species = species,
params = pars) # a complex matrix network
prn = loadrxnetwork(cmn)
# test the two networks are the same
@assert rs == prn.rnThe basic usages are
mn = MatrixNetwork(rateexprs, substoich, prodstoich; species = Any[],
params = Any[], t = nothing)
prn = loadrxnetwork(mn::MatrixNetwork)
cmn = ComplexMatrixNetwork(rateexprs, stoichmat, incidencemat; species = Any[],
params = Any[], t = nothing)
prn = loadrxnetwork(cmn::ComplexMatrixNetwork)Here MatrixNetwork and ComplexMatrixNetwork are the types, which select that
we are constructing a substrate/product stoichiometric matrix-based or a
reaction complex matrix-based stoichiometric representation as input. See the
Catalyst.jl API for more
discussion on these matrix representations, and how Catalyst handles symbolic
reaction rate expressions. These two types have the following fields:
-
rateexprs, any valid Symbolics.jl expression for the rates, or any basic number type. This can be a hardcoded rate constant like1.0, a parameter likek1above, or an general Symbolics expression involving parameters and species likek*A. -
matrix inputs
-
For
MatrixNetworksubstoich, a number of species by number of reactions matrix with entry(i,j)giving the stoichiometric coefficient of speciesias a substrate in reactionj.prodstoich, a number of species by number of reactions matrix with entry(i,j)giving the stoichiometric coefficient of speciesias a product in reactionj.
-
For
ComplexMatrixNetworkstoichmat, the complex stoichiometry matrix defined here.incidencemat, the complex incidence matrix defined here.
-
-
species, an optional vector of symbolic variables representing each species in the network. Can be constructed using the Catalyst.jl@speciesmacro. Each species should be dependent on the same time variable (tin the example above). -
parameters, a vector of symbolic variables representing each parameter in the network. Can be constructed with the ModelingToolkit.jl@parametersmacro. If no parameters are used it is an optional keyword. -
t, an optional Symbolics.jl variable representing time as the independent variable of the reaction network. If not providedCatalyst.DEFAULT_IVis used to determine the default time variable.
For both input types, loadrxnetwork returns a ParsedReactionNetwork, prn,
with only the field, prn.rn, filled in. prn.rn corresponds to the generated
Catalyst.jl
ReactionSystem
that represents the network.
Dispatches are added if substoich and prodstoich both have the type
SparseMatrixCSCin case of MatrixNetwork (or stoichmat and incidencemat
both have the type SparseMatrixCSC in case of ComplexMatrixNetwork), in
which case they are efficiently iterated through using the SparseArrays
interface.
If the keyword argument species is not set, the resulting reaction network
will simply name the species S1, S2,..., SN for a system with N total
species. params defaults to an empty vector, so that it does not need to be
set for systems with no parameters.
![dependabot[bot] avatar](https://avatars.githubusercontent.com/in/29110?v=4 "dependabot[bot]")
![github-actions[bot] avatar](https://avatars.githubusercontent.com/in/15368?v=4 "github-actions[bot]")
React
Typescript
Django
Laravel
Facebook
Microsoft
Google
Alibaba
D3
Tencent