This project contains scripts allowing to add a "virus biomass function" to a cell metabolic model (a SBML model), and then to perform an analysis of this "Host-Virus Model" (HVM).

The scripts allowing to create the HVM and to analyse it are stored in the modelAnalysis module. The analyse_model.py script in the root directory of this project consists in a command-line interface to this module. An example of invokation (using the Recon2.2 model) would be

python analyse_model.py SARS_CoV_2.txt recon2.xml biomass_reaction

where "SARS_CoV_2.txt" is the genbank file for SARS_CoV_2, "recon2.xml" is the SBML file for the Recon2.2 model and "biomass_reaction" is the id of the biomass reaction in the Recon2.2 model. For informations about the other options of the command-line interface of analyse_model.py, refer to its documentation.

The minimalMediaCreation directory contains a script to determine the composition of the minimal medium allowing growth for a given model.

The scripts create_reaction.py and get_tissue-specific_stoichiometry.sh in the root directory are used to produce modified versions of the Recon2.2 model where the biomass_protein reaction is modified such that its stoichiometry reflects the RNA expression data from a specific tissue, from the human protein atlas project (https://www.proteinatlas.org/). In order to use it, use the get_tissue-specific_stoichiometry.sh script.

The scripts have been used in the following article: hadrien delattre, Kalesh Sasidharan, orkun soyer. Inhibiting the reproduction of COVID-19-causing SARS-CoV-2 through perturbations in human cell metabolic network. Authorea. May 12, 2020. Preprint available here.

The analysis codes provided here are “AS IS" and covered under a BSD2 Licence. They are intended solely for non-commercial, academic use in the hope that they will be useful for analysing other host-virus pairs in the context of infection or viral ecology. If you would like to explore application of the provided code and analysis approach in a commercial setting, please get in touch with the project developers or University Warwick tech transfer office at “[email protected]” quoting reference ‘host-virus metabolic modelling