michaelwitting/wormLipidPredictR
Lipid reaction prediction based on WormJam

Package index
Search the michaelwitting/wormLipidPredictR package
Vignettes
Functions
34
Source code
16
Man pages
13
Home
/
GitHub
/
michaelwitting/wormLipidPredictR
/

In michaelwitting/wormLipidPredictR: Lipid reaction prediction based on WormJam 

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>",
  warning = FALSE
)

This vignette describes the complete workflow for the prediction of Caenorhabditis elegans (C. elegans) lipids performed by Witting et al. (in preparation). All template reactions at the current stage have been curated from the WormJam model.

# load required library
library(LipidNetworkPredictR)
library(dplyr)
library(rgoslin)

In order to store all the results a data frame for the reactions is created. Lipids abbreviations are just appended.

## create data frame to store reactions
lipid_reactions <- data.frame()

## create vector to store shorthand notations of generated lipids
##lipid_abbrev <- vector()
lipid_df <- tibble()

The following abbreviations are used (in alphabetical order):

Fatty acids used for biosynthesis

First input for the prediction are fatty acids that were shown to be present in C. elegans by different groups. A consensus list was generated based on fatty acids from WormJam and detected in C. elegans. A data frame of these fatty acids is supplied with this package. Fatty acids used in glycero- and glycerophospholipids are different from the ones used in sphingolipid biosynthesis. The pools are only partially overlapping between the two. Therefore, two different files for the two different pools are supplied. The following lines read the fatty acids into a data frame and isolates the shorthand notation of the fatty acids, which are further used for downstream processing.

## read data frame with fatty acids detected in C. elegans
f <- system.file("extdata", "celegans_fatty_acids_GL_GP.txt",
    package = "LipidNetworkPredictR")

fatty_acid_gl_gp_df <- read.table(f, sep = "\t", header = TRUE, 
        comment.char = "", stringsAsFactors = FALSE) |> 
    data.frame()

## get short hand notations for fatty acids
FA_GPGL <- fatty_acid_gl_gp_df$shorthand

## take only a subset of those to speed things up
FA_GPGL <- FA_GPGL[FA_GPGL %in%
    c("FA(12:0)", "FA(14:0)", "FA(15:0)", "FA(16:0)", "FA(17:0)", "FA(18:0)", "FA(19:0)", "FA(20:0)")]
lipid_info <- rgoslin::parseLipidNames(FA_GPGL)

## add lipids to data frame
lipid_df <- dplyr::bind_rows(lipid_df, lipid_info)

The lines below perform the same action for fatty acids used in sphingolipid biosynthesis.

## read data frame with fatty acids detected in C. elegans
f <- system.file("extdata", "celegans_fatty_acids_SPH.txt",
    package = "LipidNetworkPredictR")

fatty_acid_sph_df <- read.table(f, sep = "\t", header = TRUE,
    comment.char = "", stringsAsFactors = FALSE) |>
    data.frame()

## get short hand notations for fatty acids
FA_SPH <- fatty_acid_sph_df$shorthand

## take only a subset of those to speed things up
FA_SPH <- FA_SPH[FA_SPH %in%
    c("FA(12:0)", "FA(14:0)", "FA(15:0)", "FA(16:0)", "FA(17:0)", "FA(18:0)", "FA(19:0)", "FA(20:0)" )]
lipid_info <- rgoslin::parseLipidNames(FA_SPH)

## add lipids to data frame
lipid_df <- dplyr::bind_rows(lipid_df, lipid_info)

GPGL pool

Block 0: Preparation of Acyl-CoA

Biosynthesis of CoA from FA (RHEA:15421)

Acyl groups are transfered in their activate form of acyl CoA (CoA). Their synthesis requires ATP. C. elegans harbors different Acyl-CoA synthetases and their substrate specificity is not established so far. In order to be able to use all defined acyls they are converted to acyl CoA. The reaction converts fatty acids (FA) to acyl-CoA.

Block 1: Biosynthesis and Degradation of Phosphaditic acid and Diacylglycerol

Phosphaditic acid (PA) and Diacylgylcerols (DG) are an important intermediate for the biosynthesis of different lipids. The represent a branching point between triacylglycerols and glycerophospholipids. They are synthesized from lysosphosphaditic and phosphaditic acid by successive acylation of glycerol 3-phosphate. The first block of reactions handles the generation of PA and DG as well as their degradation via monoacylglycerols to fatty acids.

Biosynthesis of LPA (RHEA:15325)

The first step to obtain DG is the biosynthesis of lysophophaditic acid from acyl CoA_GPGL and glycerol 3-phosphate. This reaction takes place in the ER and the mitochondria. In the ER the reaction is catalyzed by acl-5 (WBGene00011106) or acl-6 (WBGene00017261), while in the mitochondria acl-4 (WBGene00018657) is used. No subcellular specificity is currently used, therefore the assignments are merged together to yield a single reaction.

Biosynthesis of PA (RHEA:19709)

In the next step a second acyl group is transferred to the LPA to produce PA. This reaction is catalyzed by acl-9 (WBGene00022646) or acl-8 (WBGene00020264) or mboa-3 (WBGene00016934) or acl-1 (WBGene00010339) or acl-2 (WBGene00011543) in the ER. For the mitochondrial counterpart no genes are known so far. acl-3 (WBGene00006491) might be a good candidate since it is predicted to locate in the mitochondria.

Biosynthesis of DG (RHEA:27429)

Diacylglycerols (DG) are produced from PA by dephosphorylation. DG are important intermediates for biosynthesis of glycerophospholipids or triacylglycerols (DG). A prominent enzyme catalyzing this reactions is LPIN-1. Plasma membrane or ER WBGene00008749 or WBGene00018756 or WBGene00020895 or WBGene00010425 nucleus WBGene00010425

DG Kinase (RHEA:10272)

DG can be phosphorylated to create PA, which is the reversal of the reaction above. Under normal conditions this reactions is inactive, but can be activated yielding PA. Normally, DG is used for synthesis of TGs and phospho- lipids, but phosphorylation yields a switch in DG to PA, which both have also different signaling functions. In this case a lipid, which has been generated previously by a different reactions is produced.

Degradation of DG (RHEA:44712, RHEA:33275, sn2mg_to_sn1mg)

Production of 1-MG (RHEA:44712)

DG are degraded by the activity of different lipases yielding mono acylglycerols. In the case of this reaction the fatty acid at sn1 position is removed yieling an sn2-MG. sn2-MG can be converted to sn1-MG, which in turn can be re-acylated at a later stage.

Production of 2-MG (RHEA:33275)

DG are degraded by the activity of different lipases yielding mono acylglycerols. In the case of this reaction the fatty acid at sn1 position is removed yieling an sn2-MG. sn2-MG can be converted to sn1-MG, which in turn can be re-acylated at a later stage.

Acyl chain migration (sn2mg_to_sn1mg)

Acyl chain migration can happen to produced 1-Acyl-MGs from 2-Acyl-MGs. Both species are in equilibrium to each other. However, 1-Acyl-MGs can be reacylated to yields DGs.

Degradation of MG (RHEA:34019, RHEA:32871)

Both, sn1- and sn2-MG, can degraded to free fatty acids and glycerol, which can be used for energy production upon conversion to Acyl-CoAs and import into mitochondria for beta-oxidation. Other possibilities include the re-use in lipid biosynthesis.

Phosphorylation and reacylation of MG (RHEA:38463, RHEA:32947, RHEA:33747)

Similar to the degradation, all MG can be reacylated to yields DG again.

Additionally, 1-Acyl-MG can be also phosphorylated to yield LPA, which is then fed into lipid biosynthesis.

Block 2: Biosynthesis and degradation of TG

Triacylglycerols (TG) are storage lipids in C. elegans and other organisms. They are stored in lipid droplets. In C. elegans they are also found in yolk for the production of eggs.

Biosynthesis of TG (RHEA:10868)

TG are synthesized by acylation of DG in the ER and sequestered into lipid droplets. Since TGs harbor three acyl chains their combinatorial spaces is large.

Degradation of TG (RHEA:33271)

TG are "symmetric" molecules, therefore from each TG two different DG can produced.

Block 3: Biosynthesis of Glycerophospholipids (PC, PE, PS)

The glycerophospholipids phosphatidylcholine, -ethanolamine and serine are synthesized from DG by addition of headgroups. Additionally they can be interconverted

Biosynthesis of PC (RHEA:32939)

Biosynthesis of PE (RHEA:32943)

Biosynthesis of PS (RHEA:45088, RHEA:27606)

PS can be formed from PC (RHEA:45088) and from PE (RHEA:27606).

Decarboxylation of PS to PE (RHEA:20828)

PE can be formed from PS species by decarboxylation. The reaction takes place only in mitochondria.

Block 4: Remodeling and Degradation of Glycerophospholipids

Phospholipase action on PC (RHEA:18689, RHEA:15801, RHEA:14445)

PLA1 cleaves of the acyl at the sn1 position producing an sn2LPC (RHEA:18689_.

PLA2 on PC (RHEA:15801).

PLD on PC (RHEA:14445).

LPC to FA_GPGL (RHEA:15177, RHEA:44696)

sn1LPC to FA_GPGL.

sn2LPC to FA_GPGL.

Reacylation of LPC to PC (Remodeling, RHEA:12937)

Only sn1LPC get re-acylated --> constraining? = remodeling!

Phospholipase action on PE (RHEA:44408)

PLA1 on PE (RHEA:44408).

PLA2 on PE (RHEA:44604).

PLD on PE (pe_to_pa).

LPE to FA_GPGL (RHEA:32967)

sn1LPE to FA_GPGL (RHEA:32967).

sn2LPE to FA_GPGL (sn2lpe_to_fa)

Reacylation of LPE to PE (Remodeling, RHEA:32995)

Only sn1LPE get re-acylated --> constraining? = remodeling!

Block 5: Biosynthesis of Glycerophospholipids (PG, PI, CL)

PG, PI and CL are synthesized based on the common precursor CDP-DG

Biosynthesis of CDP-DG (RHEA:16229)

Biosynthesis of PI (RHEA:11580)

Biosynthesis of PGP and PG (RHEA:12593, RHEA:33751, RHEA:32931)

This reaction is taking place only in the mitochondria.

This reaction has no gene assignment at the moment (RHEA:33751).

The reaction from PG to CL is only taking place in mitochondria (RHEA:32931).

Block 6: Biosynthesis of Etherlipids

Generation of fatty alcohols (RHEA:52716)

Some etherlipids... mostly PE, small amounts in PC. This reaction is constrained to 16:0, 18:0, 20:0!

Generation of Acyl-DHAP (RHEA:17657)

Generation of Alkyl-DHAP (RHEA:36171)

Reduction of Alkyl-DHAP to LPA-O (RHEA:36175)

Enzyme not known so far in C. elegans!!! EC 1.1.1.101

Biosynthesis of PA-O (RHEA:36235)

Same genes as for LPA to PA --> substrate specificity not known yet.

Biosynthesis of DG-O (RHEA:36239)

Biosnythesis of PE-O (RHEA:36187)

Headgroup transfer via cept-1 and cept-2.

Biosynthesis of PE-P (RHEA:22956)

reductase is unknown so far.

Biosynthesis of PC-O (RHEA:36179)

same as for PC!

Block 7: Remodeling of Etherlipids

PLA2 on PE-O (peo_to_lpeo)

Reacylation of PE-O (lpeo_to_peo)

PLA2 on PE-P (RHEA:36195)

Reacylation of PE-P (RHEA:16245)

PLA2 on PC-O (RHEA:36231)

Remodeling of PC-O???

Degradation of LPEP

Block Y: Biosynthesis and Degradation of N-Acyl-Ethanolamines

This biosynthetic pathway is putative. Only a few enzymes have been found so far.

Biosynthesis of N-Acyl-phosphatidylethanolamine (NAPE) from PE and PC (pe_to_nape_sn1, peo_to_napeo_sn1, pep_to_napep_sn1, pe_to_nape_sn2, peo_to_napeo_sn2, pep_to_napep_sn2)

N-acyl-phosphatidylethanolamines are produced by the transfer of an acyl group from PC to the nitrogen of PE. In C. elegans this reaction is so far not known. It is unknown if the sn1 or sn2 position is transferred. Since C. elegans contains about equal amounts of PE, PE-O and PE-P it might be possible that the acyl group is transferred to all of them producing NAPE, NAPE-O and NAPE-P.

transfer of Sn1 fatty acyl to PE (pe_to_nape_sn1)

transfer of Sn1 fatty acyl to PE-O (peo_to_napeo_sn1)

transfer of Sn1 fatty acyl to PE-P (pep_to_napep_sn1)

transfer of Sn2 fatty acyl to PE (pe_to_nape_sn2)

transfer of Sn2 fatty acyl to PE-O (peo_to_napeo_sn2)

transfer of Sn2 fatty acyl to PE-P (pep_to_napep_sn2)

Biosynthesis of N-Acyl-ethanolamine (nape_to_nae, napeo_to_nae)

From NAPE (nape_to_nae)

From NAPEO (napeo_to_nae)

Biosynthesis of from NAPE (nape_to_pnae)

From NAPE.

Generation of N-Acyl-Ethanolamines from N-Acyl-Phosphoethanolamines

Generation of LysoNAPE from NAPE (nape_to_lnape)

Generation of GPNAE from LysoNAPE (RHEA:45420)

Degradation of N-Acyl-Ethanolamines (nae_to_fa)

Run the function

reactions <- rbind(
    ## Block 0
    c(1, "RHEA:15421", "M_ATP + M_CoA + M_FA = M_PPi + M_AMP + M_AcylCoA"), ## FA_GPGL + CoA = PPi + AMP + CoA_GPGL
    ## Block 1
    c(2, "RHEA:15325", "M_Glycerol-3-P + M_AcylCoA = M_CoA + M_LPA"), ## G3P + CoA_GPGL = CoA + LPA
    c(3, "RHEA:19709", "M_LPA + M_AcylCoA = M_CoA + M_PA"), ## LPA + CoA_GPGL = CoA + PA
    c(4, "RHEA:27429", "M_H2O + M_PA = M_Pi + M_1,2-DG"), ## H2O + PA = Pi + DG
    c(5, "RHEA:10272", "M_ATP + M_1,2-DG = M_H+ + M_ADP + M_PA"), ## ATP + DG = H+ + ADP + PA
    c(6, "RHEA:44712", "M_H2O + M_1,2-DG = M_H+ + M_1-MG + M_FA"), ## H2O + DG = H+ + sn1MG + FA
    c(7, "RHEA:33275", "M_H2O + M_1,2-DG = M_H+ + M_2-MG + M_FA"), ## H2O + DG = H+ + sn2MG + FA
    c(8, "sn2mg_to_sn1mg", "M_2-MG <=> M_1-MG"), ## sn1MG <=> sn2MG
    c(9, "RHEA:34019", "M_H2O + M_1-MG = M_Glycerol + M_H+ + M_FA"), ## H2O + sn1MG = Glycerol + H+ + FA
    c(10, "RHEA:32871", "M_H2O + M_2-MG = M_Glycerol + M_H+ + M_FA"), ## H2O + sn2MG = Glycerol + H+ + FA
    c(11, "RHEA:38463", "M_1-MG + M_AcylCoA = M_CoA + M_1,2-DG"), ## sn1MG + CoA_GPGL = CoA + DG
    c(12, "RHEA:32947", "M_AcylCoA + M_1-MG = M_CoA + M_1,2-DG"), #################### CoA_GPGL + sn2MG = CoA + DG #################### should be 2-MG???
    c(13, "RHEA:33747", "M_ATP + M_1-MG = M_H+ + M_ADP + M_LPA"), ## ATP + 1-MG = H+ + ADP + LPA
    ## Block 2
    c(14, "RHEA:10868", "M_AcylCoA + M_1,2-DG = M_CoA + M_TG"), ## CoA_GPGL + DG = CoA + TG
    c(15, "RHEA:33271", "M_H2O + M_TG = M_H+ + M_FA + M_1,2-DG"), ## H2O + TG = H+ + FA + DG
    ## Block 3
    c(16, "RHEA:32939", "M_CDP-Choline + M_1,2-DG = M_H+ + M_CMP + M_PC"), ## CDP-Choline + DG = H+ + CMP + PC
    c(17, "RHEA:32943", "M_CDP-Ethanolamine + M_1,2-DG = M_H+ + M_CMP + M_PE"), ## CDP-Ethanolamine + DG = H+ + CMP + PE
    c(18, "RHEA:45088", "M_L-Serine + M_PC = M_Choline + M_PS"), ## L-Serine + PC = Choline + PS
    c(19, "RHEA:27606", "M_L-Serine + M_PE = M_Ethanolamine + M_PS"), ## L-Serine + PE = Ethanolamine + PS
    c(20, "RHEA:20828", "M_H+ + M_PS = M_CO2 + M_PE"), ## H+ + PS = CO2 + PE
    ## Block 4
    c(21, "RHEA:18689", "M_H2O_c + M_PC_c = M_2-LPC_c + M_FA_c"), ## H2O + PC = sn2LPC + FA_GPGL
    c(22, "RHEA:15801", "M_H2O + M_PC = M_1-LPC + M_FA"), ## H2O + PC = sn1LPC + FA_GPGL
    c(23, "RHEA:14445", "M_H2O + M_PC = M_Choline + M_PA"), ## H2O + PC = Choline + PA
    c(24, "RHEA:15177", "M_H2O + M_1-LPC = M_Glycerophosphocholine + M_H+ + M_FA"), ## H2O + sn1LPC = Glycerophosphocholine + H+ + FA_GPGL 
    c(25, "RHEA:44696", "M_H2O + M_2-LPC = M_Glycerophosphocholine + M_H+ + M_FA"), ## H2O + sn2LPC = Glycerophosphocholine + H+ + FA_GPGL
    c(26, "RHEA:12937", "M_1-LPC_c + M_AcylCoA_c = M_PC_c + M_CoA_c"), ## sn1LPC + CoA_GPGL = PC + CoA_GPGL
    c(27, "RHEA:44408", "M_H2O + M_PE = M_H+ + M_FA + M_2-LPE"), ## H2O + PE = H+ + FA_GPGL + sn2LPE
    c(28, "RHEA:44604", "M_H2O + M_PE = M_H+ + M_FA + M_1-LPE"), ## H2O + PE = H+ + FA_GPGL + sn1LPE
    c(29, "pe_to_pa", "M_H2O + M_PE <=> M_Ethanolamine + M_H+ + M_PA"), ## H2O + PE <=> Ethanolamine + H+ + PA
    c(30, "RHEA:32967", "M_H2O + M_1-LPE = M_H+ + M_FA + M_Glycerophosphoethanolamine"), ## H2O + sn1LPE = H+ + FA_GPGL + Glycerophosphoethanolamine
    c(31, "sn2lpe_to_fa", "M_H2O + M_2-LPE = M_H+ + M_FA + M_Glycerophosphoethanolamine"), ## H2O + sn2LPE = H+ + FA_GPGL + Glycerophosphoethanolamine
    c(32, "RHEA:32995", "M_AcylCoA + M_1-LPE = M_CoA + M_PE"), ## CoA_GPGL + sn1LPE = CoA + PE
    ## Block 5
    c(33, "RHEA:16229", "M_CTP + M_PA = M_PPi + M_CDP-DG"), ## CTP + PA = PPi + CDP-DG
    c(34, "RHEA:11580", "M_myo-Inositol + M_CDP-DG = M_H+ + M_CMP + M_PI"), ## myo-Inositol + CDP-DG = H+ + CMP + PI
    c(35, "RHEA:12593", "M_Glycerol-3-P + M_CDP-DG = M_H+ + M_CMP + M_PGP"), ## Glycerol-3-P + CDP-DG = H+ + CMP + PGP
    c(36, "RHEA:33751", "M_H2O + M_PGP = M_Pi + M_PG"), ## H2O + PGP = Pi + PG
    c(37, "RHEA:32931", "M_CDP-DG + M_PG = M_H+ + M_CMP + M_CL"), ## CDP-DG + PG = H+ + CMP + CL
    ## Block 6
    c(38, "RHEA:52716", "M_AcylCoA + 2 M_NADPH + 2 M_H+ = M_FAOH + 2 M_NADP + M_CoA"), ## AcylCoA_GPGL + 2 NADPH + 2 H+ = FAO + 2 NADP + CoA
    c(39, "RHEA:17657", "M_Dihydroxyacetone-P + M_AcylCoA = M_CoA + M_AcylDHAP"), ## Dihydroxyacetone-P + CoA_GPGL = CoA + AcylDHAP 
    c(40, "RHEA:36171", "M_AcylDHAP + M_FAO = M_H+ + M_FA + M_AlkylDHAP"), ## AcylDHAP + FAO = H+ + FA_GPGL + M_AlkylDHAP
    c(41, "RHEA:36175", "M_H+ + M_NADPH + M_AlkylDHAP = M_LPA-O + M_NADP"), ## H+ + NADPH + AlkylDHAP = LPAO + NADP
    c(42, "RHEA:36235", "M_LPA-O + M_AcylCoA = M_PA-O + M_CoA"), ## LPAO + CoA_GPGL = PAO + CoA
    c(43, "RHEA:36239", "M_H2O + M_PA-O = M_Pi + M_DG-O"), ## H2O + PAO = Pi + DGO
    c(44, "RHEA:36187", "M_CDP-Ethanolamine + M_DG-O = M_H+ + M_CMP + M_PE-O"), ## CDP-Ethanolamine + DGO = H+ + CMP + PEO
    c(45, "RHEA:22956", "M_PE-O + M_Fe2+cytochrome_b5 + 2 M_H+ + M_O2 = M_PE-P + M_Fe3+cytochrome_b5 + 2 M_H2O"), ## PEO + Fe2+cytochrome_b5 + 2 H+ + O2 = PEP + Fe3+cytochrome_b5 + 2 H2O
    c(46, "RHEA:36179", "M_CDP-Choline + M_DG-O = M_H+ + M_CMP + M_PC-O"), ## CDP-Choline + DGO = H+ + CMP + PCO
    ## Block 7
    c(47, "peo_to_lpeo", "M_H2O + M_PE-O = M_H+ + M_LPE-O+ M_FA"), ## H2O + PEO = H+ + LPEO + FA_GPGL
    c(48, "lpeo_to_peo", "M_AcylCoA + M_LPE-O = M_CoA + M_PE-O"), ## CoA_GPGL + LPEO = CoA + M_PEO
    c(49, "RHEA:36195", "M_H2O + M_PE-P = M_H+ + M_LPE-P + M_FA"), ## H2O + PEP = H+ + LPEP + FA
    c(50, "RHEA:16245", "M_AcylCoA + M_LPE-P = M_CoA + M_PE-P"), ## CoA_GPGL + LPEP = CoA + PEP
    c(51, "RHEA:36231", "M_H2O + M_PC-O = M_H+ + M_LPC-O + M_FA"), ## H2O + PCO = H+ + LPCO + FA_GPGL
    ## Block Y
    c(52, "pe_to_nape_sn1", "M_PE + M_PC = M_NAPE + M_2-LPC"), ## PE + PC = NAPE + sn2LPC"
    c(53, "peo_to_napeo_sn1", "M_PE-O + M_PC = M_NAPEO + M_2-LPC"), ## PEO + PC = NAPEO + sn2LPC
    c(54, "pep_to_napep_sn1", "M_PE-P + M_PC = M_NAPEP + M_2-LPC"), ## PEP + PC = NAPEP + sn2LPC
    c(55, "pe_to_nape_sn2", "M_PE + M_PC = M_NAPE + M_1-LPC"), ## PE + PC = NAPE + sn1LPC
    c(56, "peo_to_napeo_sn2", "M_PE-O + M_PC = M_NAPEO + M_1-LPC"), ## PEO + PC = NAPEO + sn1LPC
    c(57, "pep_to_napep_sn2", "M_PE-P + M_PC = M_NAPEP + M_1-LPC"), ## PEP + PC = NAPEP + sn1LPC
    c(58, "nape_to_nae", "M_NAPE + M_H2O <=> M_NAE + M_PA"), ## NAPE + H2O <=> NAE + PA
    c(59, "napeo_to_nae", "M_NAPEO + M_H2O <=> M_NAE + M_PA-O"), ## akac2nape + H2O <=> NAE + PAO
    c(60, "nape_to_pnae", "M_LNAPE + M_H2O <=> M_pnae + M_1,2-DG"), ## LNAPE + H2O <=> PNAE + DG
    c(61, "nape_to_lnape", "M_LNAPE + M_H2O <=> M_LNAPE + M_FA"), ## LNAPE + H2O <=> LNAPE + FA_GPGL
    c(62, "RHEA:45420", "M_LNAPE + M_H2O = M_GPNAE + M_FA"), ## LNAPE + H2O = GPNAE + FA_GPGL
    c(63, "nae_to_fa", "M_H2O + M_NAE = M_Ethanolamine + M_H+ + M_FA") ## H2O + nae = Ethanolamine + H+ + FA_GPGL
)

reactions <- data.frame(order = as.numeric(reactions[, 1]), RHEA = reactions[, 2],
    reactions = reactions[, 3])

## run the function
reactions_GPGL <- create_reactions(substrates = list(FA = FA_GPGL), reactions = reactions)

## show the first entry of reactions_GPGL
reactions_GPGL[[1]]

The object reactions_GPGL is a list that contains 63 entries, each consisting of another two list entries. The first entry will contain the products of each reaction and the second entry will contain a data.frame with the following information on each reaction:

The function create_reaction_adjacency_matrix will create an adjacency matrix from the reactions_GPGL. This adjacency matrix may be used in subsequent analysis, e.g. using the igraph package. 

adj_GPGL <- create_reaction_adjacency_matrix(reactions_GPGL)

SPH pool

Block 0: Preparation of Acyl-CoA

Biosynthesis from FA to CoA (RHEA:15421)

Likewise, the fatty acyls for the synthesis of sphingolipids are required as acyl-CoA. The main differences are that fatty acyls in sphingolipids are generally saturated and longer compared to glycero- and glycerophospholipids and often contain a hydroxygroup at the second position of the acyl.

Block X: Biosynthesis and degradation of sphingolipids

Biosynthesis of Dihydroceramides (RHEA:53424)

Biosynthesis of Dihydrosphingomyelins (RHEA:44620)

Degradation of Dihydrosphingomyelins (DhSM, RHEA:45300)

Biosynthesis of Ceramides (Cer, dhcer_to_cer)

Biosynthesis of Sphingomyelins (SM, RHEA:18765)

Degradation of Sphingomyelins (SM) to Ceramides (Cer) (sm_to_cer)

Biosynthesis of Ceramide-1-Phosphates (CerP) from Ceramides (Cer) (RHEA:17929)

Degradation of Ceramide-1-Phosphates (CerP) to Ceramides (Cer) (CerP_to_cer)

Gene is unknown.

Biosynthesis of Glucosylceramides (RHEA:12088)

Run the function

reactions <- rbind(
    c(1, "RHEA:15421", "M_ATP + M_CoA + M_FA = M_PPi + M_AMP + M_AcylCoA"), ## ATP + CoA + FA_SPH = PPi + AMP + CoA_SPH
    ## Block X
    c(2, "RHEA:53424", "M_AcylCoA + M_Sphinganine = M_CoA + M_DhCer"), ## CoA_SPH + c17isosphgn <=> CoA + DHCer
    c(3, "RHEA:44620", "M_Phosphocholine + M_DhCer = M_1,2-DG + M_DhSM"), ## Phosphocholine + DhCer = DG + DhSM
    c(4, "RHEA:45300", "M_H2O + M_DhSM = M_Phosphocholine + M_H+ + M_DhCer"), ## H2O + DhSM = Phosphocholine + H+ + DhCer
    c(5, "dhcer_to_cer", "M_H+ + M_NADH + M_O2 + M_DhCer = 2 M_H2O + M_NAD + M_Cer"), ## H+ + NADH + O2 + DhCer = 2 H2O + NAD + Cer
    c(6, "RHEA:18765", "M_Phosphocholine + M_Cer = M_1,2-DG + M_SM"), ## Phosphocholine + Cer = 1,2-DG + SM
    c(7, "sm_to_cer", "M_H2O + M_SM = M_Phosphocholine + M_H+ + M_Cer"), ## H2O + SM = Phosphocholine + H+ + Cer
    c(8, "RHEA:17929", "M_ATP + M_Cer = M_H+ + M_ADP + M_CerP"), ## ATP + Cer = H+ + ADP + CerP
    c(9, "cerp_to_cer", "M_H2O + M_CerP = M_Pi + M_Cer"), ## H2O + CerP = Pi + Cer
    c(10, "RHEA:12088", "M_UDP-Glucose + M_Cer = M_H+ + M_UDP + GlcCer") ## UDP-Glucose + Cer = H+ + UDP + GlcCer
)

reactions <- data.frame(order = as.numeric(reactions[, 1]), RHEA = reactions[, 2],
    reactions = reactions[, 3])

## run the function
reactions_SPH <- create_reactions(substrates = list(FA = FA_SPH), reactions = reactions)

## show the first entry of reactions_SPH
reactions_SPH[[1]]

The object reactions_SPH is a list that contains 10 entries, each consisting of another two list entries. The first entry will contain the products of each reaction and the second entry will contain a data.frame with the following information on each reaction:

The function create_reaction_adjacency_matrix will create an adjacency matrix from the reactions_GPGL. This adjacency matrix may be used in subsequent analysis, e.g. using the igraph package. 

adj_GPGL <- create_reaction_adjacency_matrix(reactions_SPH)

List of implemented reactions

Please find here a list of the implemented reactions. The list is ordered according to their RHEA ids.

cerp_to_cer: - "cerp_to_cer": "M_H2O + M_CerP <=> M_Pi + M_Cer"

dhcer_to_cer - "dhcer_to_cer": "M_H+ + M_NADH + M_O2 + M_DhCer <=> 2 M_H2O + M_NAD + M_Cer"

sn2lpe_to_fa: - "sn2lpe_to_fa": "M_H2O + M_2-LPE <=> M_H+ + M_FA + M_Glycerophosphoethanolamine"

lpeo_to_peo: - "lpeo_to_peo": "M_AcylCoA + M_ak2lgpe <=> M_CoA + M_PE-O"

lpep_to_pep: - "RHEA:16245": "M_AcylCoA + M_LPE-P = M_CoA + M_PE-P" - "RHEA:16246": "M_AcylCoA + M_LPE-P => M_CoA + M_PE-P" - "RHEA:16247": "M_AcylCoA + M_LPE-P <= M_CoA + M_PE-P" - "RHEA:16248": "M_AcylCoA + M_LPE-P <=> M_CoA + M_PE-P"

sn2mg_to_sn1mg: - "sn2mg_to_sn1mg": "M_2-MG <=> M_1-MG"

nae_to_fa: - "nae_to_fa": "M_H2O + M_NAE <=> M_Ethanolamine + M_H+ + M_FA"

nape_to_lnape
- "nape_to_lnape": "M_NAPE + M_H2O <=> M_LNAPE + M_FA"

nape_to_nae
- "nape_to_nae": "M_NAPE + M_H2O <=> M_NAE + M_PA"

nape_to_pnae
- "nape_to_pnae": "M_NAPE + M_H2O <=> M_PNAE + M_1,2-DG"

napeo_to_nae
- "napeo_to_nae": "M_NAPEO + M_H2O <=> M_NAE + M_PA-O"

pco_to_lpco:
- "RHEA:36231": "M_H2O + M_PC-O <=> M_H+ + M_LPC-O + M_FA" - "RHEA:36232": "M_H2O + M_PC-O <=> M_H+ + M_LPC-O + M_FA" - "RHEA:36233": "M_H2O + M_PC-O <=> M_H+ + M_LPC-O + M_FA" - "RHEA:36234": "M_H2O + M_PC-O <=> M_H+ + M_LPC-O + M_FA"

pe_to_dg:
- "pe_to_dg": "M_H2O + M_PE <=> M_P-Ethanolamine + M_1,2-DG"

pe_to_nape_sn1:
- "pe_to_nape_sn1": "M_PE + M_PC <=> M_NAPE + M_2-LPC"

pe_to_nape_sn2:
- "pe_to_nape_sn2": "M_PE + M_PC <=> M_NAPE + M_1-LPC"

pe_to_pa - "pe_to_pa": "M_H2O + M_PE <=> M_Ethanolamine + M_H+ + M_PA"

peo_to_lpeo: - "peo_to_lpeo": "M_H2O + M_PE-O <=> M_H+ + M_LPE-O + M_FA"

peo_to_napeo_sn1
- "peo_to_napeo_sn1": "M_PE-O + M_PC <=> M_NAPEO + M_2-LPC"

peo_to_napeo_sn2
- "peo_to_napeo_sn2": "M_PE-O + M_PC <=> M_NAPEO + M_1-LPC"

pep_to_lpep:
- "RHEA:36195": "M_H2O + M_PE-P = M_H+ + M_LPE-P + M_FA" - "RHEA:36196": "M_H2O + M_PE-P => M_H+ + M_LPE-P + M_FA" - "RHEA:36197": "M_H2O + M_PE-P <= M_H+ + M_LPE-P + M_FA" - "RHEA:36198": "M_H2O + M_PE-P <=> M_H+ + M_LPE-P + M_FA"

pep_to_napep_sn1:
- "pep_to_napep_sn1": "M_PE-P + M_PC <=> M_NAPEP + M_2-LPC"

pep_to_napep_sn2:
- "pep_to_napep_sn2": "M_PE-P + M_PC <=> M_NAPEP + M_1-LPC"

dg_to_pa: - "RHEA:10272": "M_ATP + M_1,2-DG = M_H+ + M_ADP + M_PA" - "RHEA:10273": "M_ATP + M_1,2-DG => M_H+ + M_ADP + M_PA" - "RHEA:10274": "M_ATP + M_1,2-DG <= M_H+ + M_ADP + M_PA" - "RHEA:10275": "M_ATP + M_1,2-DG <=> M_H+ + M_ADP + M_PA"

pc_to_dg - "RHEA:10604": "M_H2O + M_PC = M_Phosphocholine + M_1,2-DG" - "RHEA:10605": "M_H2O + M_PC => M_Phosphocholine + M_1,2-DG" - "RHEA:10606": "M_H2O + M_PC <= M_Phosphocholine + M_1,2-DG" - "RHEA:10607": "M_H2O + M_PC <=> M_Phosphocholine + M_1,2-DG"

dg_to_tg: - "RHEA:10868": "M_AcylCoA + M_1,2-DG = M_CoA + M_TG" - "RHEA:10869": "M_AcylCoA + M_1,2-DG => M_CoA + M_TG" - "RHEA:10870": "M_AcylCoA + M_1,2-DG <= M_CoA + M_TG" - "RHEA:10871": "M_AcylCoA + M_1,2-DG <=> M_CoA + M_TG"

cdpdg_to_pi: - "RHEA:11580": "M_myo-Inositol + M_CDP-DG = M_H+ + M_CMP + M_PI" - "RHEA:11581": "M_myo-Inositol + M_CDP-DG => M_H+ + M_CMP + M_PI" - "RHEA:11582": "M_myo-Inositol + M_CDP-DG <= M_H+ + M_CMP + M_PI" - "RHEA:11583": "M_myo-Inositol + M_CDP-DG <=> M_H+ + M_CMP + M_PI"

cer_to_glccer: - "RHEA:12088": "M_UDP-Glucose + M_Cer <=> M_H+ + M_UDP + M_GlcCer"

cdpdg_to_pgp: - "RHEA:12593": "M_Glycerol-3-P + M_CDP-DG = M_H+ + M_CMP + M_pgp" - "RHEA:12594": "M_Glycerol-3-P + M_CDP-DG => M_H+ + M_CMP + M_PGP" - "RHEA:12595": "M_Glycerol-3-P + M_CDP-DG <= M_H+ + M_CMP + M_PGP" - "RHEA:12596": "M_Glycerol-3-P + M_CDP-DG <=> M_H+ + M_CMP + M_PGP"

sn1lpc_to_pc: - "RHEA:12937": "M_1-LPC + M_AcylCoA = M_PC + M_CoA" - "RHEA:12938": "M_1-LPC + M_AcylCoA => M_PC + M_CoA" - "RHEA:12939": "M_1-LPC + M_AcylCoA <= M_PC + M_CoA" - "RHEA:12940": "M_1-LPC + M_AcylCoA <=> M_PC + M_CoA"

pc_to_pa - "RHEA:14445": "M_H2O + M_PC = M_Choline + M_PA" - "RHEA:14446": "M_H2O + M_PC => M_Choline + M_PA" - "RHEA:14447": "M_H2O + M_PC <= M_Choline + M_PA" - "RHEA:14448": "M_H2O + M_PC <=> M_Choline + M_PA"

sn1lpc_to_fa: - "RHEA:15177": "M_H2O + M_1-LPC = M_Glycerophosphocholine + M_H+ + M_FA" - "RHEA:15178": "M_H2O + M_1-LPC => M_Glycerophosphocholine + M_H+ + M_FA" - "RHEA:15179": "M_H2O + M_1-LPC <= M_Glycerophosphocholine + M_H+ + M_FA" - "RHEA:15180": "M_H2O + M_1-LPC <=> M_Glycerophosphocholine + M_H+ + M_FA"

coa_to_lpa: - "RHEA:15325": "M_Glycerol-3-P + M_AcylCoA = M_CoA + M_LPA" - "RHEA:15326": "M_Glycerol-3-P + M_AcylCoA => M_CoA + M_LPA" - "RHEA:15327": "M_Glycerol-3-P + M_AcylCoA <= M_CoA + M_LPA" - "RHEA:15328": "M_Glycerol-3-P + M_AcylCoA <=> M_CoA + M_LPA"

fa_to_coa - "RHEA:15421"/"RHEA:38883": "M_ATP + M_CoA + M_FA = M_PPi + M_AMP + M_AcylCoA" - "RHEA:15422"/"RHEA:38884": "M_ATP + M_CoA + M_FA => M_PPi + M_AMP + M_AcylCoA" - "RHEA:15423"/"RHEA:38885": "M_ATP + M_CoA + M_FA <= M_PPi + M_AMP + M_AcylCoA" -"RHEA:15424"/"RHEA:38886": "M_ATP + M_CoA + M_FA <=> M_PPi + M_AMP + M_AcylCoA"

pc_to_sn1lpc - "RHEA:15801": "M_H2O + M_PC = M_1-LPC + M_FA" - "RHEA:15802": "M_H2O + M_PC => M_1-LPC + M_FA" - "RHEA:15803": "M_H2O + M_PC <= M_1-LPC + M_FA" - "RHEA:15804": "M_H2O + M_PC <=> M_1-LPC + M_FA"

pa_to_cdpdg: - "RHEA:16229": "M_CTP + M_PA = M_PPi + M_CDP-DG" - "RHEA:16230": "M_CTP + M_PA => M_PPi + M_CDP-DG" - "RHEA:16231": "M_CTP + M_PA <= M_PPi + M_CDP-DG" - "RHEA:16232": "M_CTP + M_PA <=> M_PPi + M_CDP-DG"

sn2mg_to_dg: - "RHEA:16741"/"RHEA:32947": "M_AcylCoA + M_1-MG = M_CoA + M_1,2-DG" - "RHEA:16742"/"RHEA:32948": "M_AcylCoA + M_1-MG => M_CoA + M_1,2-DG" - "RHEA:16743"/"RHEA:32949": "M_AcylCoA + M_1-MG <= M_CoA + M_1,2-DG" - "RHEA:16744"/"RHEA:32950": "M_AcylCoA + M_1-MG <=> M_CoA + M_1,2-DG"

lpep_to_fal: - RHEA:16905: "H2O + 1-LPE-P = FAL + Glycerophosphoethanolamine" - RHEA:16906: "H2O + 1-LPE-P => FAL + Glycerophosphoethanolamine" - RHEA:16907: "H2O + 1-LPE-P <= FAL + Glycerophosphoethanolamine" - RHEA:16908: "H2O + 1-LPE-P <=> FAL + Glycerophosphoethanolamine"

coa_to_acdhap: - "RHEA:17657": "M_Dihydroxyacetone-P + M_AcylCoA = M_CoA + M_AcylDHAP" - "RHEA:17658": "M_Dihydroxyacetone-P + M_AcylCoA => M_CoA + M_AcylDHAP" - "RHEA:17659": "M_Dihydroxyacetone-P + M_AcylCoA <= M_CoA + M_AcylDHAP" - "RHEA:17660": "M_Dihydroxyacetone-P + M_AcylCoA <=> M_CoA + M_AcylDHAP"

cer_to_cerp: - "RHEA:17929": "M_ATP + M_Cer <=> M_H+ + M_ADP + M_CerP"

pi_to_sn1lpi: - "RHEA:18001": "M_H2O + M_PI = M_1-LPI + M_H+ + M_FA" - "RHEA:18002": "M_H2O + M_PI => M_1-LPI + M_H+ + M_FA" - "RHEA:18003": "M_H2O + M_PI <= M_1-LPI + M_H+ + M_FA" - "RHEA:18004": "M_H2O + M_PI <=> M_1-LPI + M_H+ + M_FA"

pc_to_sn2lpc - "RHEA:18689": "M_H2O + M_PC = M_2-LPC + M_FA" - "RHEA:18690": "M_H2O + M_PC => M_2-LPC + M_FA" - "RHEA:18691": "M_H2O + M_PC <= M_2-LPC + M_FA" - "RHEA:18692": "M_H2O + M_PC <=> M_2-LPC + M_FA"

cer_to_sm: - "RHEA:18765": "M_PC + M_Cer <=> M_1,2-DG + M_SM"

lpa_to_pa: - "RHEA:19709": "M_LPA + M_AcylCoA = M_CoA + M_PA" - "RHEA:19710": "M_LPA + M_AcylCoA => M_CoA + M_PA" - "RHEA:19711": "M_LPA + M_AcylCoA <= M_CoA + M_PA" - "RHEA:19712": "M_LPA + M_AcylCoA <=> M_CoA + M_PA"

ps_to_pe: - "RHEA:20828": "M_H+ + M_PS = M_CO2 + M_PE" - "RHEA:20829": "M_H+ + M_PS => M_CO2 + M_PE" - "RHEA:20830": "M_H+ + M_PS <= M_CO2 + M_PE" - "RHEA:20831": "M_H+ + M_PS <=> M_CO2 + M_PE"

peo_to_pep: - "RHEA:22956": "M_PE-O + M_Fe2+-cytochrome_b5 + 2 M_H+ + M_O2 = M_PE-P + M_Fe3+-cytochrome_b5 + 2 M_H2O" - "RHEA:22957": "M_PE-O + M_Fe2+-cytochrome_b5 + 2 M_H+ + M_O2 => M_PE-P + M_Fe3+-cytochrome_b5 + 2 M_H2O" - "RHEA:22958": "M_PE-O + M_Fe2+-cytochrome_b5 + 2 M_H+ + M_O2 <= M_PE-P + M_Fe3+-cytochrome_b5 + 2 M_H2O" - "RHEA:22959": "M_PE-O + M_Fe2+-cytochrome_b5 + 2 M_H+ + M_O2 <=> M_PE-P + M_Fe3+-cytochrome_b5 + 2 M_H2O"

lpco_to_pco: - "RHEA:23992": "1-LPC-O + AcylCoA = PC-O + CoA" - "RHEA:23993": "1-LPC-O + AcylCoA => PC-O + CoA" - "RHEA:23994": "1-LPC-O + AcylCoA <= PC-O + CoA" - "RHEA:23995": "1-LPC-O + AcylCoA <=> PC-O + CoA"

pa_to_dg: - "RHEA:27429": "M_H2O + M_PA = M_Pi + M_1,2-DG" - "RHEA:27430": "M_H2O + M_PA => M_Pi + M_1,2-DG" - "RHEA:27431": "M_H2O + M_PA <= M_Pi + M_1,2-DG" - "RHEA:27432": "M_H2O + M_PA <=> M_Pi + M_1,2-DG"

pe_to_ps: - "RHEA:27606": "M_L-Serine + M_PE = M_Ethanolamine + M_PS" - "RHEA:27607": "M_L-Serine + M_PE => M_Ethanolamine + M_PS" - "RHEA:27608": "M_L-Serine + M_PE <= M_Ethanolamine + M_PS" - "RHEA:27609": "M_L-Serine + M_PE <=> M_Ethanolamine + M_PS"

sn2mg_to_fa: - "RHEA:32871": "M_H2O + M_2-MG = M_Glycerol + M_H+ + M_FA" - "RHEA:32872": "M_H2O + M_2-MG => M_Glycerol + M_H+ + M_FA" - "RHEA:32873": "M_H2O + M_2-MG <= M_Glycerol + M_H+ + M_FA" - "RHEA:32874": "M_H2O + M_2-MG <=> M_Glycerol + M_H+ + M_FA"

pg_to_cl: - "RHEA:32931": "M_CDP-DG + M_PG = M_H+ + M_CMP + M_CL" - "RHEA:32932": "M_CDP-DG + M_PG => M_H+ + M_CMP + M_CL" - "RHEA:32933": "M_CDP-DG + M_PG <= M_H+ + M_CMP + M_CL" - "RHEA:32934": "M_CDP-DG + M_PG <=> M_H+ + M_CMP + M_CL"

cl_to_lcl: - "RHEA:32935": "M_H2O + M_CL = M_1,2,4-LCL + M_H+ + M_FA" - "RHEA:32936": "M_H2O + M_CL => M_1,2,4-LCL + M_H+ + M_FA" - "RHEA:32937": "M_H2O + M_CL <= M_1,2,4-LCL + M_H+ + M_FA" - "RHEA:32938": "M_H2O + M_CL <=> M_1,2,4-LCL + M_H+ + M_FA"

dg_to_pc - "RHEA:32939": "M_CDP-Choline + M_1,2-DG = M_H+ + M_CMP + M_PC" - "RHEA:32940": "M_CDP-Choline + M_1,2-DG => M_H+ + M_CMP + M_PC" - "RHEA:32941": "M_CDP-Choline + M_1,2-DG <= M_H+ + M_CMP + M_PC" - "RHEA:32942": "M_CDP-Choline + M_1,2-DG <=> M_H+ + M_CMP + M_PC"

dg_to_pe: - "RHEA:32943": "M_CDP-Ethanolamine + M_1,2-DG = M_H+ + M_CMP + M_PE" - "RHEA:32944": "M_CDP-Ethanolamine + M_1,2-DG => M_H+ + M_CMP + M_PE" - "RHEA:32945": "M_CDP-Ethanolamine + M_1,2-DG <= M_H+ + M_CMP + M_PE" - "RHEA:32946": "M_CDP-Ethanolamine + M_1,2-DG <=> M_H+ + M_CMP + M_PE"

sn1lpe_to_fa: - "RHEA:32967": "M_H2O + M_1-LPE = M_H+ + M_FA + M_Glycerophosphoethanolamine" - "RHEA:32968": "M_H2O + M_1-LPE => M_H+ + M_FA + M_Glycerophosphoethanolamine" - "RHEA:32969": "M_H2O + M_1-LPE <= M_H+ + M_FA + M_Glycerophosphoethanolamine" - "RHEA:32970": "M_H2O + M_1-LPE <=> M_H+ + M_FA + M_Glycerophosphoethanolamine"

sn1lpe_to_pe: - "RHEA:32995": "M_AcylCoA + M_1-LPE = M_CoA + M_PE" - "RHEA:32996": "M_AcylCoA + M_1-LPE => M_CoA + M_PE" - "RHEA:32997": "M_AcylCoA + M_1-LPE <= M_CoA + M_PE" - "RHEA:32998": "M_AcylCoA + M_1-LPE <=> M_CoA + M_PE"

sn1lpi_to_pi: - "RHEA:33195": "M_1-LPI + M_AcylCoA = M_PI + M_CoA" - "RHEA:33196": "M_1-LPI + M_AcylCoA => M_PI + M_CoA" - "RHEA:33197": "M_1-LPI + M_AcylCoA <= M_PI + M_CoA" - "RHEA:33198": "M_1-LPI + M_AcylCoA <=> M_PI + M_CoA"

tg_to_dg: - "RHEA:33271"/"RHEA:44864": "M_H2O + M_TG = M_H+ + M_FA + M_1,2-DG" - "RHEA:33272"/"RHEA:44865": "M_H2O + M_TG => M_H+ + M_FA + M_1,2-DG" - "RHEA:33273"/"RHEA:44866": "M_H2O + M_TG <= M_H+ + M_FA + M_1,2-DG" - "RHEA:33274"/"RHEA:44867": "M_H2O + M_TG <=> M_H+ + M_FA + M_1,2-DG"

dg_to_sn2mg - "RHEA:33275": "M_H2O + M_1,2-DG = M_H+ + M_2-MG + M_FA" - "RHEA:33276": "M_H2O + M_1,2-DG => M_H+ + M_2-MG + M_FA" - "RHEA:33277": "M_H2O + M_1,2-DG <= M_H+ + M_2-MG + M_FA" - "RHEA:33278": "M_H2O + M_1,2-DG <=> M_H+ + M_2-MG + M_FA"

sn1mg_to_lpa - "RHEA:33747": "M_ATP + M_1-MG = M_H+ + M_ADP + M_LPA" - "RHEA:33748": "M_ATP + M_1-MG => M_H+ + M_ADP + M_LPA" - "RHEA:33749": "M_ATP + M_1-MG <= M_H+ + M_ADP + M_LPA" - "RHEA:33750": "M_ATP + M_1-MG <=> M_H+ + M_ADP + M_LPA"

pgp_to_pg: - "RHEA:33751": "M_H2O + M_pgp = M_Pi + M_PG" - "RHEA:33752": "M_H2O + M_pgp => M_Pi + M_PG" - "RHEA:33753": "M_H2O + M_pgp <= M_Pi + M_PG" - "RHEA:33754": "M_H2O + M_pgp <=> M_Pi + M_PG"

sn1mg_to_fa: - "RHEA:34019": "M_H2O + M_1-MG = M_Glycerol + M_H+ + M_FA" - "RHEA:34020": "M_H2O + M_1-MG => M_Glycerol + M_H+ + M_FA" - "RHEA:34021": "M_H2O + M_1-MG <= M_Glycerol + M_H+ + M_FA" - "RHEA:34022": "M_H2O + M_1-MG <=> M_Glycerol + M_H+ + M_FA"

dg_to_sn1mg: - "RHEA:35663"/"RHEA:44712": "M_H2O + M_1,2-DG = M_H+ + M_1-MG + M_FA" - "RHEA:35664"/"RHEA:44713": "M_H2O + M_1,2-DG => M_H+ + M_1-MG + M_FA" - "RHEA:35665"/"RHEA:44714": "M_H2O + M_1,2-DG <= M_H+ + M_1-MG + M_FA" - "RHEA:35666"/"RHEA:44715": "M_H2O + M_1,2-DG <=> M_H+ + M_1-MG + M_FA"

lcl_to_cl: - "RHEA:35839": "M_1,2,4-LCL + M_AcylCoA = M_CL + M_CoA" - "RHEA:35840": "M_1,2,4-LCL + M_AcylCoA => M_CL + M_CoA" - "RHEA:35841": "M_1,2,4-LCL + M_AcylCoA <= M_CL + M_CoA" - "RHEA:35842": "M_1,2,4-LCL + M_AcylCoA <=> M_CL + M_CoA"

lpco_to_mgo: - "RHEA:36083": "H2O + 1-LPC-O = H+ + Phosphocholine + 1-MG-O" - "RHEA:36084": "H2O + 1-LPC-O => H+ + Phosphocholine + 1-MG-O" - "RHEA:36085": "H2O + 1-LPC-O <= H+ + Phosphocholine + 1-MG-O" - "RHEA:36086": "H2O + 1-LPC-O <=> H+ + Phosphocholine + 1-MG-O"

acyldhap_to_alkyldhap: - "RHEA:36171": "M_AcylDHAP + M_FAO = M_H+ + M_FA + M_AlkylDHAP" - "RHEA:36172": "M_AcylDHAP + M_FAO => M_H+ + M_FA + M_AlkylDHAP" - "RHEA:36173": "M_AcylDHAP + M_FAO <= M_H+ + M_FA + M_AlkylDHAP" - "RHEA:36174": "M_AcylDHAP + M_FAO <=> M_H+ + M_FA + M_AlkylDHAP"

alkyldhap_to_lpao: - "RHEA:36175": "M_H+ + M_NADPH + M_AlkylDHAP = M_LPA-O + M_NADP" - "RHEA:36176": "M_H+ + M_NADPH + M_AlkylDHAP => M_LPA-O + M_NADP" - "RHEA:36177": "M_H+ + M_NADPH + M_AlkylDHAP <= M_LPA-O + M_NADP" - "RHEA:36178": "M_H+ + M_NADPH + M_AlkylDHAP <=> M_LPA-O + M_NADP"

dgo_to_pco - "RHEA:36179": "M_CDP-Choline + M_DG-O <=> M_H+ + M_CMP + M_PC-O" - "RHEA:36180": "M_CDP-Choline + M_DG-O <=> M_H+ + M_CMP + M_PC-O" - "RHEA:36181": "M_CDP-Choline + M_DG-O <=> M_H+ + M_CMP + M_PC-O" - "RHEA:36182": "M_CDP-Choline + M_DG-O <=> M_H+ + M_CMP + M_PC-O"

dgo_to_peo: - "RHEA:36187": "M_CDP-Ethanolamine + M_DG-O = M_H+ + M_CMP + M_PE-O" - "RHEA:36188": "M_CDP-Ethanolamine + M_DG-O => M_H+ + M_CMP + M_PE-O" - "RHEA:36189": "M_CDP-Ethanolamine + M_DG-O <= M_H+ + M_CMP + M_PE-O" - "RHEA:36190": "M_CDP-Ethanolamine + M_DG-O <=> M_H+ + M_CMP + M_PE-O"

lpep_to_lpap: - "RHEA:36203": "H2O + 1-LPE-P = LPA-P + H+ + Ethanolamine" - "RHEA:36204": "H2O + 1-LPE-P => LPA-P + H+ + Ethanolamine" - "RHEA:36205": "H2O + 1-LPE-P <= LPA-P + H+ + Ethanolamine" - "RHEA:36206": "H2O + 1-LPE-P <=> LPA-P + H+ + Ethanolamine"

lpao_to_pao: - "RHEA:36235": "M_LPA-O + M_AcylCoA = M_PA-O + M_CoA" - "RHEA:36236": "M_LPA-O + M_AcylCoA => M_PA-O + M_CoA" - "RHEA:36237": "M_LPA-O + M_AcylCoA <= M_PA-O + M_CoA" - "RHEA:36238": "M_LPA-O + M_AcylCoA <=> M_PA-O + M_CoA"

lpep_to_mgp: - "RHEA:36199": "H2O + 1-LPE-P = Phosphoethanolamine + H+ + 1-MG-P" - "RHEA:36200": "H2O + 1-LPE-P => Phosphoethanolamine + H+ + 1-MG-P" - "RHEA:36201": "H2O + 1-LPE-P <= Phosphoethanolamine + H+ + 1-MG-P" - "RHEA:36202": "H2O + 1-LPE-P <=> Phosphoethanolamine + H+ + 1-MG-P"

pao_to_dgo - "RHEA:36239": "M_H2O + M_PA-O = M_Pi + M_DG-O" - "RHEA:36240": "M_H2O + M_PA-O => M_Pi + M_DG-O" - "RHEA:36241": "M_H2O + M_PA-O <= M_Pi + M_DG-O" - "RHEA:36242": "M_H2O + M_PA-O <=> M_Pi + M_DG-O"

sn1mg_to_dg: - "RHEA:38463"/"RHEA:39943": "M_1-MG + M_AcylCoA = M_CoA + M_1,2-DG" - "RHEA:38464"/"RHEA:39944": "M_1-MG + M_AcylCoA => M_CoA + M_1,2-DG" - "RHEA:38465"/"RHEA:39945": "M_1-MG + M_AcylCoA <= M_CoA + M_1,2-DG" - "RHEA:38466"/"RHEA:39946": "M_1-MG + M_AcylCoA <=> M_CoA + M_1,2-DG"

lpco_to_lpao: - "RHEA:39927": "M_H2O + M_1-LPC-O = M_1-LPA-O + M_H+ + M_Choline" - "RHEA:39928": "M_H2O + M_1-LPC-O => M_1-LPA-O + M_H+ + M_Choline" - "RHEA:39929": "M_H2O + M_1-LPC-O <= M_1-LPA-O + M_H+ + M_Choline" - "RHEA:39930": "M_H2O + M_1-LPC-O <=> M_1-LPA-O + M_H+ + M_Choline"

pi_to_dg: - "RHEA:43484": "M_H2O + M_PI = M_myo-Inositol-1-P + M_H+ + M_1,2-DG" - "RHEA:43485": "M_H2O + M_PI => M_myo-Inositol-1-P + M_H+ + M_1,2-DG" - "RHEA:43486": "M_H2O + M_PI <= M_myo-Inositol-1-P + M_H+ + M_1,2-DG" - "RHEA:43487": "M_H2O + M_PI <=> M_myo-Inositol-1-P + M_H+ + M_1,2-DG"

pe_to_sn2lpe: - "RHEA:44408": "M_H2O + M_PE = M_H+ + M_FA + M_2-LPE" - "RHEA:44409": "M_H2O + M_PE => M_H+ + M_FA + M_2-LPE" - "RHEA:44410": "M_H2O + M_PE <= M_H+ + M_FA + M_2-LPE" - "RHEA:44411": "M_H2O + M_PE <=> M_H+ + M_FA + M_2-LPE"

pe_to_sn1lpe: - "RHEA:44604": "M_H2O + M_PE = M_H+ + M_FA + M_1-LPE" - "RHEA:44605": "M_H2O + M_PE => M_H+ + M_FA + M_1-LPE" - "RHEA:44606": "M_H2O + M_PE <= M_H+ + M_FA + M_1-LPE" - "RHEA:44607": "M_H2O + M_PE <=> M_H+ + M_FA + M_1-LPE"

dhcer_to_dhsm - "RHEA:44620": "M_PC + M_DhCer <=> M_1,2-DG + M_DhSM"

sn2lpc_to_fa: - "RHEA:44696": "M_H2O + M_2-LPC = M_Glycerophosphocholine + M_H+ + M_FA" - "RHEA:44697": "M_H2O + M_2-LPC => M_Glycerophosphocholine + M_H+ + M_FA" - "RHEA:44698": "M_H2O + M_2-LPC <= M_Glycerophosphocholine + M_H+ + M_FA" - "RHEA:44699": "M_H2O + M_2-LPC <=> M_Glycerophosphocholine + M_H+ + M_FA"

pc_to_ps - "RHEA:45088": "M_L-Serine + M_PC = M_Choline + M_PS" - "RHEA:45089": "M_L-Serine + M_PC => M_Choline + M_PS" - "RHEA:45090": "M_L-Serine + M_PC <= M_Choline + M_PS" - "RHEA:45091": "M_L-Serine + M_PC <=> M_Choline + M_PS"

dhsm_to_dhcer - "RHEA:45300": "M_H2O + M_DhSM = M_Phosphocholine + M_H+ + M_DhCer" - "RHEA:45301": "M_H2O + M_DhSM => M_Phosphocholine + M_H+ + M_DhCer" - "RHEA:45302": "M_H2O + M_DhSM <= M_Phosphocholine + M_H+ + M_DhCer" - "RHEA:45303": "M_H2O + M_DhSM <=> M_Phosphocholine + M_H+ + M_DhCer"

lnape_to_gpnae: - "RHEA:45420": "M_LNAPE + M_H2O <=> M_gpnae + M_FA"

coa_to_FAO: - "RHEA:52716": "M_AcylCoA + 2 M_NADPH + 2 M_H+ = M_FAO + 2 M_NADP + M_CoA" - "RHEA:52717": "M_AcylCoA + 2 M_NADPH + 2 M_H+ => M_FAO + 2 M_NADP + M_CoA" - "RHEA:52718": "M_AcylCoA + 2 M_NADPH + 2 M_H+ <= M_FAO + 2 M_NADP + M_CoA" - "RHEA:52719": "M_AcylCoA + 2 M_NADPH + 2 M_H+ <=> M_FAO + 2 M_NADP + M_CoA"

sm_to_cer:
- "sm_to_cer": "M_H2O + M_SM <=> M_Phosphocholine + M_H+ + M_Cer"

sphinga_to_dhcer:
- "RHEA:53424": "M_AcylCoA + M_Sphinganine = M_CoA + M_DhCer" - "RHEA:53425": "M_AcylCoA + M_Sphinganine => M_CoA + M_DhCer" - "RHEA:53426": "M_AcylCoA + M_Sphinganine <= M_CoA + M_DhCer" - "RHEA:53427": "M_AcylCoA + M_Sphinganine <=> M_CoA + M_DhCer"



michaelwitting/wormLipidPredictR documentation built on July 1, 2023, 9:20 p.m.

R Package Documentation

rdrr.io home R language documentation Run R code online

Browse R Packages

CRAN packages Bioconductor packages R-Forge packages GitHub packages

We want your feedback!

Note that we can't provide technical support on individual packages. You should contact the package authors for that.
GitHub issue tracker
ian@mutexlabs.com
Personal blog

 
Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.

Close