INDRA (Integrated Network and Dynamical Reasoning Assembler) is an automated model assembly system, originally developed for molecular systems biology and then generalized to other domains (see INDRA World). INDRA draws on natural language processing systems and structured databases to collect mechanistic and causal assertions, represents them in a standardized form (INDRA Statements), and assembles them into various modeling formalisms including causal graphs and dynamical models.
At the core of INDRA are its knowledge-level assembly procedures, allowing sources to be assembled into coherent models, a process that involves correcting systematic input errors, finding and resolving redundancies, inferring missing information, filtering to a relevant scope and assessing the reliability of causal information.
The detailed INDRA documentation is available at http://indra.readthedocs.io.
INDRA is currently integrated with the following natural language processing
systems and structured databases. These input modules (available in
indra.sources) all produce INDRA Statements.
Biological pathway databases:
|Database / Exchange format||Module||Reference|
|PathwayCommons / BioPax||http://pathwaycommons.org/ |
|Large Corpus / BEL||https://github.com/pybel/pybel |
|Target Affinity Spectrum||https://doi.org/10.1101/358978|
Custom knowledge bases:
|Database / Exchange format||Module||Reference|
|NDEx / CX||http://ndexbio.org|
|INDRA DB / INDRA Statements||https://github.com/indralab/indra_db|
INDRA also provides several model output assemblers that take INDRA Statements as input. The most sophisticated model assembler is the PySB Assembler, which implements a policy-guided automated assembly procedure of a rule-based executable model (that can then be further compiled into other formats such as SBML, Kappa, BNGL and SBGN to connect to a vast ecosystem of downstream tools). Several other model assembly modules target various network formats for visualization, and graph/structural analysis (PyBEL, CyJS, Graphviz, SBGN, CX, SIF, etc.) and curation (HTML, TSV, IndexCards). Finally, the English Assembler produces English language descriptions of a set of INDRA Statements.
INDRA also supports extension by outside model assembly tools which take INDRA Statements as input and produce models. One such example is Delphi (https://github.com/ml4ai/delphi), which is a Dynamic Bayesian Network model assembler. Similarly, outside tools that support INDRA Statements can implement custom visualization methods, such as CauseMos, developed by Uncharted Software (https://uncharted.software/).
Assemblers aimed at model-driven discovery and analysis:
|Modeling formalism / Exchange format||Purpose||Module||Reference|
|PySB (-> SBML, SBGN, BNGL, Kappa, etc.)||Detailed, mechanistic modeling, simulation, analysis||http://pysb.org|
|PyBEL||Causal analysis, visualization||https://github.com/pybel/pybel |
|IndraNet||Causal analysis, signed and unsigned|
|SIF||Network analysis, logic modeling, visualization||SIF format|
|Figaro||Bayesian network inference||https://github.com/p2t2/figaro/|
|KAMI||Knowledge aggregation of protein sites/states and Kappa modeling||https://github.com/Kappa-Dev/KAMI|
Assemblers primarily aimed at visualization:
|Network / Exchange format||Purpose||Module||Reference|
|Causal Analysis Graph||General causal graph visualization|
|CX||Network browsing, versioning on NDEx||http://ndexbio.org|
|Cytoscape JS||Interactive Cytoscape JS network to embed in websites||http://js.cytoscape.org/|
|Graphviz||Static PDF/PNG visualization with powerful automated layout using Graphviz||https://www.graphviz.org/|
|SBGN||Visualization with Systems Biology Graphical Notation||http://sbgn.org|
Assemblers primarily aimed at expert curation and browsing:
|English language||Human-readable descriptions, reports, dialogue|
|HTML||Web-based browsing, linking out to provenance, curation||Curation tutorial|
|TSV (Tab/Comma Separated Values)||Spreadsheet-based browsing and curation|
|Index Cards||Custom JSON format for curating biological mechanisms|
A key feature of INDRA is providing internal knowledge-assembly modules that operate on INDRA Statements and perform the following tasks:
The internal assembly steps of INDRA including the ones listed above, and also a large collection of filters (filter by source, belief, the presence of grounding information, semantic filters by entity role, etc.) are exposed in the indra.tools.assemble_corpus submodule. This submodule contains functions that take Statements as input and produce processed Statements as output. They can be composed to form an assembly pipeline connecting knowledge collected from sources with an output model.
This diagram illustrates the assembly pipeline process.
The choice of assembly functions can vary depending on the domain (i.e, biology or world modeling), the modeling goal (i.e., the type of model that will be assembled and how that model will be used), desired features, and confidence (e.g., filter to human genes only or apply a belief cutoff), and any other user preferences.
An example of a typical assembly pipeline for biology statements is as follows. Some of the below steps can be removed, rearranged, and other steps added to change the assembly pipeline.
from indra.tools import assemble_corpus as ac stmts = <the collection of all raw statements to use> stmts = ac.filter_no_hypothesis(stmts) # Filter out hypothetical statements stmts = ac.map_grounding(stmts) # Map grounding stmts = ac.filter_grounded_only(stmts) # Filter out ungrounded agents stmts = ac.filter_human_only(stmts) # Filter out non-human genes stmts = ac.map_sequence(stmts) # Map sequence stmts = ac.run_preassembly(stmts, # Run preassembly return_toplevel=False) stmts = ac.filter_belief(stmts, 0.8) # Apply belief cutoff of 0.8
An example of an assembly pipeline for statements in the world modeling domain
is as follows (note how biology-specific functions are not used, and a custom
belief_scorer and ontology is passed to
run_preassembly here, while the
biology pipeline used default values). Note that this example requires
indra_world package to be installed.
from indra.tools import assemble_corpus as ac from indra_world.belief.wm_scorer import get_eidos_scorer from indra_world.ontology.world import world_ontology stmts = <the collection of all raw statements to use> stmts = ac.filter_grounded_only(stmts) # Filter out ungrounded agents belief_scorer = get_eidos_scorer() stmts = ac.run_preassembly(stmts, # Run preassembly return_toplevel=False, belief_scorer=belief_scorer, ontology=world_ontology, normalize_equivalences=True, # Optional: rewrite equivalent groundings to one standard normalize_opposites=True, # Optional: rewrite opposite groundings to one standard normalize_ns='WM') # Use 'WM' namespace to normalize equivalences and opposites stmts = ac.filter_belief(stmts, 0.8) # Apply belief cutoff of e.g., 0.8
Assembled statements returned after running the assembly pipeline can be passed into any of the output model assemblers.
INDRA also contains modules to access literature content (e.g., PubMed, Elsevier), available in
access ontological information and convert between identifiers (e.g., UniProt,
HGNC), available in
A full list of further INDRA modules is available in the
Gyori B.M., Bachman J.A., Subramanian K., Muhlich J.L., Galescu L., Sorger P.K. From word models to executable models of signaling networks using automated assembly (2017), Molecular Systems Biology, 13, 954.
For detailed installation instructions, see the documentation.
INDRA currently supports Python 3.6+. The last release of INDRA compatible with Python 2.7 is 1.10, and the last release fully compatible with Python 3.5 is 1.17.
The preferred way to install INDRA is by pointing pip to the source repository as
$ pip install git+https://github.com/sorgerlab/indra.git
Releases of INDRA are also available on PyPI, you can install the latest release as
$ pip install indra
However, releases will usually be behind the latest code available in this repository.
INDRA depends on a few standard Python packages. These packages are installed by pip during setup. For certain modules and use cases, other "extra" dependencies may be needed, which are described in detail in the documentation.
A REST API for INDRA is available at http://api.indra.bio:8000. Note that the REST API is ideal for prototyping and for building light-weight web apps, but should not be used for large reading and assembly workflows.
INDRA is available as a Docker image on Dockerhub and can be pulled as
docker pull labsyspharm/indra
You can run the INDRA REST API using the container as
docker run -id -p 8080:8080 --entrypoint python labsyspharm/indra /sw/indra/rest_api/api.py
The Dockerfile to build the image locally is available in a separate repository at https://github.com/indralab/indra_docker.
In this example INDRA assembles a PySB model from the natural language description of a mechanism via the TRIPS reading web service.
from indra.sources import trips from indra.assemblers.pysb import PysbAssembler pa = PysbAssembler() # Process a natural language description of a mechanism trips_processor = trips.process_text('MEK2 phosphorylates ERK1 at Thr-202 and Tyr-204') # Collect extracted mechanisms in PysbAssembler pa.add_statements(trips_processor.statements) # Assemble the model model = pa.make_model(policies='two_step')
INDRA also provides an interface for the
REACH natural language
processor. In this example, a full paper from PubMed
Central is processed. The paper's PMC ID is
PMC3717945. The example
assumest that a REACH server is running locally (see documentation at
Note that REACH takes about 8 minutes to read this full-text paper.
from indra.sources import reach reach_processor = reach.process_pmc('PMC3717945', url=reach.local_nxml_url)
At this point,
reach_processor.statements contains a list of INDRA statements
extracted from the PMC paper.
Next we look at an example of reading the 10 most recent PubMed abstracts on BRAF and collecting the results in INDRA statements.
from indra.sources import reach from indra.literature import pubmed_client # Search for 10 most recent abstracts in PubMed on 'BRAF' pmids = pubmed_client.get_ids('BRAF', retmax=10) all_statements =  for pmid in pmids: abs = pubmed_client.get_abstract(pmid) if abs is not None: reach_processor = reach.process_text(abs, url=reach.local_text_url) if reach_processor is not None: all_statements += reach_processor.statements
At this point, the
all_statements list contains all the statements
extracted from the 10 abstracts.
The next example shows querying the BEL large corpus network for a neighborhood of a given list of proteins using their HGNC gene names. This example performs the query via PyBEL.
from indra.sources import bel # Process the neighborhood of BRAF and MAP2K1 bel_processor = bel.process_pybel_neighborhood(['BRAF', 'MAP2K1'])
At this point,
bel_processor.statements contains a list of INDRA statements
extracted from the neighborhood query.
Next, we look at an example of querying the Pathway Commons database for paths between two lists of proteins.
from indra.sources import biopax # Process the neighborhood of BRAF and MAP2K1 biopax_processor = biopax.process_pc_pathsfromto(['BRAF', 'RAF1'], ['MAP2K1', 'MAP2K2'])
At this point,
biopax_processor.statements contains a list of INDRA
Statements extracted from the paths-from-to query.
The development of INDRA has been funded from the following sources:
|DARPA Big Mechanism||W911NF-14-1-0397|
|DARPA World Modelers||W911NF-18-1-0014|
|DARPA Communicating with Computers||W911NF-15-1-0544|
|DARPA Automated Scientific Discovery Framework||W911NF-18-1-0124|
|DARPA Automating Scientific Knowledge Extraction||HR00111990009|
|DARPA Young Faculty Award||W911NF-20-1-0255|