Characteristics of available models for <i>Phaeodactylum tricornutum</i>.

<p>Characteristics of available models for <i>Phaeodactylum tricornutum</i>.</p

Photosynthetic electron flow constraints as determined by Bailleul <i>et al</i>. [83].

<p>Photosynthetic electron flow constraints as determined by Bailleul <i>et al</i>. [<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0155038#pone.0155038.ref083" target="_blank">83</a>].</p

Chloroplastic ornithine cycle as revealed by the model.

<p>Metabolic network usage of a chloroplastic ornithine cycle is shown under a saturating photon constraint of 16 mM allowing maximum carbon uptake. Minor reactants and products are omitted for visual clarity (i.e., water, protons and phosphate). Metabolite and reaction abbreviation suffixes indicate cellular compartment; c, cytosol; h, chloroplast; m, mitochondria. Reversible reactions are indicated by arrowheads at both ends. The filled arrowhead indicates the direction in which the reaction is running, i.e. from substrate (open arrowhead) to product (filled arrowhead). Abbreviations used: ACOAT, acetylornithine transaminase; AGK, acetylglutamate kinase; AGPR, N-acetyl-δ-glutamyl-phosphate reductase; GACT, glutamate N-acetyltransferase; GLNA, glutamine synthase; GLTS, glutamate synthase (ferredoxin dependent); GLUDH2, glutamine dehydrogenase (NAD dependent); GLUSA, glutamate semialdehyde degradation (spontaneous); OAT, ornithine aminotransferase; P5CDH, 1-pyrroline-5-carboxylate dehydrogenase; acorn, N-acetylornithine; acglu, N-acetyl-L-glutamate; acg5p, N-acetyl-L-glutamate 5-phosphate; acg5sa, N-Acetyl-L-glutamate 5-semialdehyde; adp, ADP; akg, α-ketoglutarate; atp, ATP; fdxox, ferredoxin (oxidized); fdxrd, ferredoxin (reduced); gln__L, L-glutamine; glu__L, L-glutamate; glu5sa, L-glutamate 5-semialdehyde; nad, NAD⁺; nadh, NADH; nadp, NADP⁺; nadph, NADPH; nh4, ammonium ion; orn, ornithine; 1pyr5c, (S)-1-Pyrroline-5-carboxylate.</p

Genes in reconstruction over predicted genes in genome against genome size for selected eukaryotic metabolic reconstructions.

<p>The three reconstructions with the highest ratio of genes in reconstruction per genes in genome are highlighted. bna572+ has a comparable ratio as <i>i</i>LB1025 and <i>i</i>LB1027_lipid, <i>i</i>TO977 has a higher ratio. Compared to <i>i</i>TO977 and bna572+, <i>i</i>LB1025 and <i>i</i>LB1027_lipid contain more reactions. The number of reactions in the respective reconstructions is used to scale the circle diameters. Note the discontinuous x-axis. Abbreviations: AraGEM: <i>Arabidopsis thaliana</i> [<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0155038#pone.0155038.ref073" target="_blank">73</a>]; bna572+: <i>Brassica napus</i> [<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0155038#pone.0155038.ref074" target="_blank">74</a>]; AlgaGEM: <i>Chlamydomonas reinhardtii</i> [<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0155038#pone.0155038.ref075" target="_blank">75</a>]; <i>i</i>RC1080: <i>Chlamydomonas reinhardtii</i> [<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0155038#pone.0155038.ref039" target="_blank">39</a>]; <i>i</i>RS1563: <i>Zea mays</i> [<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0155038#pone.0155038.ref076" target="_blank">76</a>]; <i>i</i>LB1025 and <i>i</i>LB1027_lipid: <i>Phaeodactylum tricornutum</i>, this study; <i>i</i>TO977: <i>Saccharomyces cerevisiae</i> Sc288 [<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0155038#pone.0155038.ref077" target="_blank">77</a>]; Recon2: <i>Homo sapiens</i> [<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0155038#pone.0155038.ref078" target="_blank">78</a>]; <i>i</i>MM1415: <i>Mus musculus</i> [<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0155038#pone.0155038.ref079" target="_blank">79</a>].</p

Light-dependent carbon partitioning.

<p>(A) Simulations indicated as photon uptake exceeds carbon uptake, excess redox potential is stored in triacylglycerol. The saturation of carbon uptake is shown in black. (B) Percent of carbon fixed in TAG against percent of metabolite flow through NADHOR (vNADHOR; EC 1.6.5.3,1.6.99.3) over metabolite flow through PSI (vPSI; EC 1.97.1.2) at a super-saturating photon uptake of 22 mM. According to our simulations TAG accumulation is inversely proportional to energetic coupling. TAG accumulation is prohibited when at least 35% of photosynthetically fixed electrons are redirected to the mitochondria.</p

Metabolic network reconstruction workflow.

<p>In step one we obtained a draft reconstruction based on <i>P</i>. <i>tricornutum</i>’s genome annotation and reference reconstructions. This draft reconstruction was manually curated using several resources such as an improved genome annotation, subcellular localization predictions and external databases. All reactions were elementally and charge balanced, QC/QA was performed and a biomass objective function was defined before transforming the reconstruction into a computational model. In an iterative process, the <i>in silico</i> predictions are compared with experimental observations to validate and improve the metabolic model.</p

Reconstruction characteristics <i>i</i>LB1027_lipid.

<p>(A) Reactions per subsystem. Most reactions are involved in lipid metabolism. Our FTIR measurements underline the fact that the lipids make up the highest fraction of biomass. Due to the presence of multiple compartments and the fact that many pathways are split among compartments, many reactions are attributed to intracellular transport. The modeling subsystem contains ATP maintenance, biomass, demand, sink, and exchange reactions. (B) Percent reactions and metabolites per compartment. Most reactions and metabolites are present in the cytosol, followed by chloroplast and mitochondria in the case of reactions and mitochondria and chloroplast for metabolites. Peroxisome, extracellular space, and thylakoid contain less than 5% and 8% of all reactions and metabolites in the reconstruction, respectively.</p

Subcellular localization prediction pipeline.

<p>Schematic representation of the implemented subcellular localization prediction pipeline for <i>Phaeodactylum tricornutum</i> adapted from previous work [<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0155038#pone.0155038.ref057" target="_blank">57</a>]. Subcellular compartments are given in ellipses and bioinformatics programs are displayed in rectangles. Our added steps are highlighted in gray. The ER retention signal is (K/D)-(D/E)-E-L in the protein C-terminal region. A protein is categorized as peroxisomal if the signal (S/A/C)-(K/R/H)-(L/M) or S-S-L is found in the C-terminal region.</p