Given the following input data:
EC numbers for the enzymes involved + compound ids of reactants involved in each reaction (we actually have to do conversion steps from gene names and arbitrary compound ids in a SBML file before in order to generate these, but I am positive that this is the easiest step).
Does anybody know or has anybody already implemented a programmatic way of determining whether or not a reaction in KEGG or maybe ecocyc metacyc is reversible or not, and would like to share the principle or code, in whatever language?
The background: We are trying to evaluate and quality check a genome scale metabolic model (in SBML) and we suspect that it erroneously has be added a single (hopefully not more) reversible reaction where the reaction in reality is irreversible. So it would be possible but tedious process to check manually with the kegg maps if the reversibility in our model conforms with the reaction in KEGG.
I hope this has become clear, if not please comment for more information.