I got a list of proteins (UniProt ID) from a proteomics experiment and i am trying to build a model containing all the routes where my proteins are implicated. To achieve this, my first idea is to extract the Reactome route identifier from all the Uniprot IDs and then obtain the SBML for each one of those routes.
Joining the models (SBMLs) is a problem for future me, but i actually dont know how to obtain the SBML from Reactome with a list of identifiers like this (R-HSA-112409). I obtained the list of identifiers with a small python script.
I'm trying to build the model, not to get the routes more represented within my list of proteins. Reactome provides this type of analysis, but the SBML resulting from the top routes does not contain the same information that those accesed from the UniProt identifiers.
This is my first contact with SBMLs, so i would like to ask for advice for constructing the model (other databases providing SBMLs given a list of proteins for example), and if it were possible, for a way to obtain the SBMLs from Reactome using the route identifiers!
Thank you for your answers!