Question About MMT Code for Building modelStorage

[Nicholas Y]

Hi,

I have a question regarding the `buildModelStorage` function in the Microbiome Modeling Toolbox code for `mgPipe`.

My question concerns the part in the file:

cobratoolbox/src/analysis/multiSpecies/microbiomeModelingToolbox/mgPipe
/buildModelStorage.m

in lines 47-94, when getting the list of active exchange metabolites.

I noticed that in line 48, the COBRA model is loaded, and at line 80, FVA is run to determine which exchange reactions can carry flux. However, when I look at the code, I do not see any diet constraints being applied to the model before FVA is run.

Therefore, I was wondering, is there a specific diet that is used to run this FVA to determine the active exchange metabolites? Or does the code intentionally use whatever exchange bounds the model happens to have when it is loaded?

I noticed that some of the AGORA2 models have exchange bounds of (-1000, 1000) by default. However, there is nothing that guarantees the input model would have all its exchange bounds set to (-1000, 1000). For example, if I were to use a custom model where the exchange bounds were all set to (-0.1, 1000) in the .mat file, would it run FVA with the bounds set to (-0.1, 1000)?

If that is the case, is that the intended behavior or is this a potential bug? If I am using custom models that were saved with custom diet constraints applied to them, should I ensure all exchange bounds are set to (-1000, 1000) in all models in the model directory before running mgPipe?

Thank you for your time and insights!

Best regards,

Nicholas Yousefi