A thorough literature search was performed to obtain an in-depth knowledge of all the major metabolic pathways known to occur in an epithelial cell of the small intestine. We then retrieved the corresponding reactions and genes from the human metabolic reconstruction (13 (
link)), which is accessible through the BiGG database (70 (
link)), to compile a draft reconstruction. Missing transport and metabolic reactions for peptides and for dietary fibers were added to the initial draft reconstruction upon detailed manual gap analysis and further review of the corresponding literature. Genome annotations from the EntrezGene database (71 (
link)) as well as protein information from the UniProt (72 (
link)) and BRENDA database (73 (
link)) were used in addition to the information retrieved from the scientific literature to assign GPR associations to the reactions not present in Recon 1. For the reactions that were extracted from the Recon 1, GPR associations were kept as reported in Recon 1, since no comprehensive transcriptomic data are available for sIECs. The sIEC metabolic reconstruction was assembled and converted to a mathematical model using rBioNet as a reconstruction environment (74 (
link)) and an established protocol (75 ).
We used the global human metabolic network, Recon 1 (13 (
link)), as reaction database, but adjusted sub-cellular and extracellular location, reaction stoichiometry and directionality according to literature evidence. Only those reactions and pathways with literature evidence for their occurrence in human small intestinal enterocytes were incorporated into hs_sIEC611 from the global human metabolic reconstruction Recon 1, which captures metabolic capabilities known to occur in any human cell. Moreover, we added 262 transport and 50 metabolic reactions, which were not present in Recon 1, but for which supporting information for their presence in sIECs could be found in the scientific literature (Fig.
1E,
Supplementary Material, Table S2). These reactions included many transport systems specific for enterocytes and metabolic pathways for sulfo-cysteine metabolism, dietary fiber metabolism, di- and tri-peptide degradation and cholesterol-ester synthesis (Fig.
1E,
Supplementary Material, Table S2). In addition to these reactions, 73 reactions were added from our recently published acylcarnitine/fatty acid oxidation module for the human metabolic reconstruction (20 (
link)). We added further 95 reactions, which were present in Recon 1, but for which the compartment was adjusted by placing them into the lumen compartment. The stoichiometry of the reactions catalyzed by the glucose 6-phosphate dehydrogenase (E.C. 1.1.1.49), the 6-phosphogluconolactonase (E.C. 3.1.1.31) and the phosphogluconate dehydrogenase (E.C. 1.1.1.44) was changed to three, since they were required to generate three molecules each of 6-phospho-
d-glucono-1,5-lactone, 6-Phospho-
d-gluconate and ribulose-5-phosphate simultaneously (76 ). The directionality of ATP requiring fatty acid activation reactions catalyzed by the fatty acyl-CoA ligase (E.C. 6.2.1.3) was changed in agreement with a recent report (76 ). Also, the cofactor requirement and sub-cellular localization of reactions included in the cholesterol synthesis pathway, which are catalyzed by the desmosterol reductase (E.C. 1.3.1.72) and HMG-CoA reductase (E.C. 2.3.3.10) reactions, were updated in accordance with the current literature evidence (76 ).
Sahoo S, & Thiele I. (2013). Predicting the impact of diet and enzymopathies on human small intestinal epithelial cells. Human Molecular Genetics, 22(13), 2705-2722.