CONTRABASS is a tool for the detection of vulnerabilities in metabolic models. The main purpose of the tool is to compute chokepoint and essential reactions by taking into account both the topology and the dynamic information of the model. In addition to the detection of vulnerabilities, CONTRABASS can compute essential genes, compute and remove dead-end metabolites, compute different sets of growth-dependent reactions, and update the flux bounds of the reactions according to the results of Flux Variability Analysis.
CONTRABASS takes as input the SBML file of a metabolic model, and provides as output a spreadsheet file and an html file reporting the obtained results. CONTRABASS accounts for the following sets of reactions and metabolites:
- Chokepoint reactions: A reaction is a chokepoint if it is the unique consumer or the only producer of a given metabolite.
- Essential Reactions: A reaction is essential if its deletion, or equivalently, restricting its flux to zero, causes a significant decrease in the objective function (e.g. cellular growth).
- Dead-End Metabolites (DEM): A metabolite is a dead-end metabolite if it is not produced or not consumed by any reaction.
- Essential reactions for optimal growth: A reaction is essential for optimal growth if its deletion, or equivalently, restricting its flux to zero, causes a decrease in the objective function.
- Dead reactions: A reaction is dead is its upper flux bound and its lower flux bound are equal to zero.
- Blocked reactions: A reaction is blocked if its flux is necessarily zero at any possible steady state of the model.
- Reversible reactions: A reaction is reversible if its upper flux bound is strictly positive and its lower flux bound is strictly negative.
- Non-reversible reactions: A reaction is non-reversible if it is not dead and not reversible.
- Essential genes: A gene is essential if the objetive function (e.g. cellular growth) is zero when it is knocked down.
Figure: In the network, reactions Reac_2, Reac_5, Reac_6, Reac_7 and Reac_8 are chokepoint reactions.
Metabolites Lysine and Glutamate are dead-end metabolites.
The computation of vulnerabilities can also be exploited programmatically via the Low Level API which is based on COBRApy. You can also try the Low Level API with Binder.
Table of Contents
Citation
[Paper]
@article{10.1093/bioinformatics/btad053,
author = {Oarga, Alexandru and Bannerman, Bridget P and Júlvez, Jorge},
title = "{CONTRABASS: exploiting flux constraints in genome-scale models for the detection of vulnerabilities}",
journal = {Bioinformatics},
volume = {39},
number = {2},
pages = {btad053},
year = {2023},
month = {01},
issn = {1367-4811},
doi = {10.1093/bioinformatics/btad053},
url = {https://doi.org/10.1093/bioinformatics/btad053},
eprint = {https://academic.oup.com/bioinformatics/article-pdf/39/2/btad053/49123769/btad053.pdf},
}
License
CONTRABASS is released under GPLv3 license.
Online Web
CONTRABASS can also be executed through the CONTRABASS Online Web App (http://contrabass.unizar.es).
Install
CONTRABASS can be installed via pip package manager:
$ pip install contrabass
Quickstart
Generate report on vulnerabilities on input model MODEL.xml
$ contrabass report critical-reactions MODEL.xml
Generate report on growth-dependent reactions on input model MODEL.xml
$ contrabass report growth-dependent-reactions MODEL.xml
Documentation
Documentation is available at readthedocs and can also be downloaded.
The previous links include examples and descriptions of the operations that can be performed with the tool.
The next flowchart provides a graphical description of the available operations that can be performed with CONTRABASS and their respective commands:
More information about the parameters of the tool can be obtained by executing contrabass -h.
For a detailed description of the operations see the documentation.
$ contrabass
Usage: contrabass [OPTIONS] COMMAND [ARGS]...
Compute vulnerabilities on constraint-based models
Options:
-h, --help Show this message and exit.
-V, --version Show the version and exit.
Commands:
new-model Export refined constraint-based model.
report Compute vulnerabilities on constraint-based models.
Low Level API
The computation of vulnerabilities can also be exploited via the COBRApy based low level API. For further information see the Low Level API documentation.
You can also try it with Binder.
Example of network refinement and chokepoint computation:
from contrabass.core import CobraMetabolicModel
model = CobraMetabolicModel("aureus.xml")
# update flux bounds with FVA
model.fva(update_flux=True)
# compute chokepoints
model.find_chokepoints()
# get chokepoints
model.chokepoints()
Maintainers
@alexOarga
Contributing
Feel free to dive in! Open an issue or submit PRs.
Standard Readme follows the Contributor Covenant Code of Conduct.
评论