Scripts provided for the Paper "The intertwined metabolism during symbiotic nitrogen fixation elucidated by metabolic modelling"
By Thomas Pfau, Nils Christian, Shyam K Masakapalli, Lee J. Sweetlove, Mark G. Poolman and Oliver Ebenhöh
Cobra Toolbox (available here) and IBM CPLEX (Academic licenses are available free of charge via the IBM academic initiative). The solver used in the paper was CPLEX v12.6.1.
All Matlab Scripts are contained in the Matlab folder.
The relevant functions to plot the different Figures are:
importMedicago.m - Imports the Medicago truncatula network from the xml file into a MATLAB COBRA structure and adjusts some encodings
BuildTissueModel.m - Builds the Two tissue model. You can supply a previously loaded Model, otherwise it is read from the xml file.
BuildSymbioticModel.m - Build the Symbiotic model and the Non Symbiotic Model wth adjusted Asparagine Biomass fraction. You can provide the TissueModel generated from BuildTissueModel to this function. Otherwise it will be generated from the xml.
AmmoniumScan.m - Generates the Ammonium scan plot (underlying Figure 2), if the Tissue Model is not supplied it will be loaded from the xml.
GrowthRateComparison.m - Generates the comparison plot between Growth Rates of the symbiotic and non symbiotic system (Figure 4). You can supply the outputs of the BuildSymbioticModel function, otherwise the respective models will be built from the xml file.
OxWeightScan.m - Generates the plot for the Oxygen scan underlying Figure 5. You can supply the SymbioticModel generated by BuildSymbioticModel to this function. Otherwise it will be constructed from the xml.
It is recommended to use importMedicago, BuildTissueModel and BuildSymbioticModel once, since otherwise extensive computations will have to be run every time the plots are generated.
The ScrumPy (available here) to Cplex interface is provided with the File ScrumPyQP.py in the ScrumPy directory.
Pfau, T., Christian, N., Masakapalli, S. K., Sweetlove, L. J., Poolman, M. G., & Ebenhöh, O. (2018). The intertwined metabolism during symbiotic nitrogen fixation elucidated by metabolic modelling. Scientific reports, 8(1), 12504.