GapFilling with Cobrapy.

[Bello Suleiman]

Hello all,

Im new to the group. I recently started using cobrapy for my genome scale metabolic models analysis. 

I have a few models that have quiet a number of blocked metabolites. 

The Cobra toolbox for Matlab has functions that can detect theses metabolites in a metabolic model network. 

Functions such as:

1. findRootNPmets (Finds both no consumption and production metabolites)
2. gapFind (Finds any blocked metabolites downstream of a gap)

My question is are there any functions similar to the ones mentioned above that are available in cobrapy. I read through the cobrapy documentations and packages.

The only form of GapFilling functions i came across was "GrowMatch" and "SMILEY" which try to to make the minimal number of changes to a model and allow it to simulate growth. 

I do not want to rely on Flux Balance Analysis as a means to determine what reactions may be suitable candidates to fill in gaps. 

Thank you in advance.