initial state
Short overview of the workflow and its principles:
- You have a model/set of reactions
- You need SMILES notation of all metabolites present in the reactions
- You calculate INCHi-keys for all metabolites
- You create SMILES strings for all reactions, each in iits own directory, including the model-IDs and InCHi-keys for the used metabolites
- You apply RDT (could be replaced by another automatic atom mapping tool) on the reaction
- You get a .rxn file as result
- you parse the individual molecules out that file
- you calculate the InCHi-key for all molecules
- you identify the model ID for the molecule sequence of the .rxn file
- you extract the mapping information
- you calculate InCHi for all individual molecules of the .rxn file, with auxillary information (containing a mapping of which InCHi-ordered atom is equivalent to which input atom
- You adapt the mapping information, so now all mappings identify atoms with their InCHi-order in the molecule
- combine the InCHi-ordered mappings from all reactions