BioMed Research International
Volume 2016 (2016), Article ID 7863706, 9 pages
http://dx.doi.org/10.1155/2016/7863706

Research Article

Reconstruction of the Fatty Acid Biosynthetic Pathway of Exiguobacterium antarcticum B7 Based on Genomic and Bibliomic Data

Regiane Kawasaki,¹ Rafael A. Baraúna,¹ Artur Silva,¹ Marta S. P. Carepo,² Rui Oliveira,² Rodolfo Marques,² Rommel T. J. Ramos,¹ and Maria P. C. Schneider¹

¹Genomics and Systems Biology Center, Institute of Biological Sciences, Federal University of Pará, 66075-110 Belém, PA, Brazil
²Rede de Química e Tecnologia/Centro de Química Fina e Biológica, Chemistry Department, Universidade Nova de Lisboa, 2829-516 Costa da Caparica, Portugal

Received 10 November 2015; Accepted 16 June 2016

Academic Editor: Yongsheng Bai

Copyright © 2016 Regiane Kawasaki et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Exiguobacterium antarcticum B7 is extremophile Gram-positive bacteria able to survive in cold environments. A key factor to understanding cold adaptation processes is related to the modification of fatty acids composing the cell membranes of psychrotrophic bacteria. In our study we show the in silico reconstruction of the fatty acid biosynthesis pathway of E. antarcticum B7. To build the stoichiometric model, a semiautomatic procedure was applied, which integrates genome information using KEGG and RAST/SEED. Constraint-based methods, namely, Flux Balance Analysis (FBA) and elementary modes (EM), were applied. FBA was implemented in the sense of hexadecenoic acid production maximization. To evaluate the influence of the gene expression in the fluxome analysis, FBA was also calculated using the values obtained in the transcriptome analysis at 0°C and 37°C. The fatty acid biosynthesis pathway showed a total of 13 elementary flux modes, four of which showed routes for the production of hexadecenoic acid. The reconstructed pathway demonstrated the capacity of E. antarcticum B7 to de novo produce fatty acid molecules. Under the influence of the transcriptome, the fluxome was altered, promoting the production of short-chain fatty acids. The calculated models contribute to better understanding of the bacterial adaptation at cold environments.

1. Introduction

Bacteria are increasingly used in industrial processes to produce chemicals, foods, and drugs, among other products [1]. The main biochemical pathways of bacteria may be manipulated and optimized to more efficiently produce compounds of industrial interest in various areas; for example, metabolic pathways of Corynebacterium glutamicum are rationally engineered to produce L-amino acids on an industrial scale [2]. To accomplish this task, specific tools are used such as FMM (from metabolite to metabolite) [3] Cytoscape [4], CellDesigner [5], SBW (Systems Biology Workbench) [6], COPASI (COmplex PAthway SImulator) [7], and COBRA (COnstraints Based Reconstruction and Analysis) toolbox [8].

The genomes of several bacterial strains have been sequenced and annotated and have been used in combination with biochemical and physiological data to reconstruct metabolic networks at a genome-scale [9]. Recently, genomic models were reconstructed for some bacterial species aiming to increase the amount and quality of data that has been annotated in either the literature or databases [10–12]. The draft network generated from the annotated genome still requires significant manual curation for a comprehensive and accurate metabolic representation of the organism [13].

The need to develop automatic or at least semiautomatic methods to reconstruct metabolic networks from genome annotation is increasing because the number of complete genome sequences available is growing fast. Recent studies [13, 14] have highlighted the problems associated with genome annotations and databases, which perform automatic reconstructions and, thus, require manual assessment. The currently available 96-step protocol, proposed by Thiele and Palsson [15], is a well-established process for the assembly, curation, and validation of metabolic reconstruction. This protocol is combined with computational tools, including the visualization and numerical calculation software package, MATLAB® (MathWorks, USA).

Constraint-based modeling is frequently used as final validation step of the reconstructed network. This step is extremely useful for simulating the phenotypic behavior under different physiological environments [16–18] thereby allowing assessing if the reconstructed network represents well the in vivo cellular system. Microbial adaptation to cold environments is one of the applications of these methods.

Psychrotrophic microorganisms have an optimal growth temperature higher than 15°C but are also able to grow and adapt to extremely cold environments, with temperatures of approximately 0°C [19]. Thus, the unique physicochemical characteristics of their habitat and the biological apparatus developed by these microorganisms to survive under these conditions render these organisms valuable sources of biotechnological processes. The cellular response to cold by psychrotrophic bacteria may be studied from a general standpoint with the advent of omics methods. Recently, the B7 strain of Exiguobacterium antarcticum was isolated from Ginger Lake sediments located in the Antarctic Peninsula Region (69°30′S, 65°W). This lake was formed due to the warming in the region, which led to partial melting of ice caps [20]. The genome of this strain was sequenced [21], and its response to cold was evaluated through differential expression of its genome at 37°C and 0°C using omics methods [22]. One of the mechanisms of cold adaptation of all psychrophilic or psychrotrophic organisms is the change in the chemical structures of the membrane phospholipids. The fatty acid chains become shorter and unsaturated at low temperatures. Accordingly, the fluidity of the membrane is kept intact [19, 23].

Bacterial de novo synthesis of fatty acids is regulated by the protein FapR [24], which is responsible for activating/disabling a regulon consisting of four operons in E. antarcticum B7. In cold, two of these operons are repressed (fabH1-fabF and fabI), and the expression levels of the other two remain unaltered [22]. The chemical components, which are included in this regulon, must be reconstructed and evaluated and then associated with their respective genetic elements to further understand this metabolic pathway and to reach more complete conclusions about its importance for adaptation to cold [23, 24]. Bioinformatics methods are used for in silico reconstruction of metabolic pathways [25].

In this work we present the in silico reconstruction of the fatty acid biosynthesis pathway of the Exiguobacterium antarcticum B7, based on linear programming (FBA) and convex cone method (elementary modes). The influence of transcriptome in FBA calculation was also evaluated.

2. Materials and Methods

2.1. Data Collection

The genomic data of E. antarcticum B7 in the formats  .gbk and  .fasta were collected from NCBI under accession number NC_018665. The metabolic pathway of the fatty acid biosynthesis was initially evaluated in the KEGG (Kyoto Encyclopedia of Genes and Genomes) database [26]. When necessary, the visualization of the genome of the bacterium was performed using Artemis software [27].

2.2. Preliminary Reconstruction

This step was performed following two methods: one semiautomatic and the other automatic. The tools within the KEGG databases were used essentially in the semiautomatic method. The.fasta file with the E. antarcticum B7 genome was submitted to the online tool KAAS (KEGG Automatic Annotation Server) [28], available at http://www.genome.jp/kaas-bin/kaas_main. The parameters chosen to run this software were as follows: (a) bidirectional Best Hit (BBH) Method, recommended for complete genomes, performs the search for orthologous genes between a specific group of organisms, and (b) prokaryote, the set of genes chosen, should be representative of the target organism, in this case, the bacterium E. antarcticum B7. Following the processing, a text file was generated (query.ko). Each line of this file is formed by two parameters: the first consists of the sequence identified (gene), and the second, when present, consists of the KO assignment, termed number. This value indicates orthologous groups encoding the same enzymatic activity. Afterwards, the file generated is passed through a filter, an auxiliary computer software program (script) Python developed for the present study, which only selects numbers and individually and increasingly commands per line into a new file (new_query.ko). This file was used as entry in the option User Data Mapping of the Pathway Mapping tool of the KEGG.

The automatic method essentially consisted of submitting the  .fasta file of the E. antarcticum B7 genome to the online tool RAST (Rapid Annotation using Subsystems Technology) [29], available at http://rast.nmpdr.org/, to generate the drafts of the metabolic network and of the fatty acid biosynthesis pathway of the target microorganism. The final draft of this step was generated from the combination of the resulting pathways of the semiautomatic and automatic models. The common pathways were maintained, while surplus compounds, enzymes, and reactions, that is, present in some, but absent in others, were not directly excluded but were instead reserved for the curated step.

2.3. Manual Curation

The following steps were completed in this manual curation stage, following the protocol explained above. (i) Draft refinement: this phase began with the analysis of enzymes and reactions, components of the fatty acid biosynthesis pathway, by reading books and articles specifically on the subject. The objective was to diagnose the absence or presence of more than one element of the study pathway. The online databases KEGG, ENZYME [30], and SEED [31] were consulted to ratify the enzymes and the structures of the reactions.

(ii) Assessment of the stoichiometry and reversibility of the reactions: in this step, all model reactions were assessed and stoichiometrically corrected, if necessary. The biochemical data on the organism are very important to determine the reversibility of the reaction. For this purpose, the databases (KEGG, SEED, and ENZYME) and the tool eQuilibrator [32] were used to analyze the reversibility of the reactions. Thus, the thermodynamic constraints were respected.

(iii) Addition of gene data and reaction location: the Artemis tool was used to identify the genes of the reactions (enzymes) from their locations in the genome assessed using the draft generated.

(iv) Assessment of Gene-Protein-Reaction (GPR) associations: in this step, the function of each gene is indicated. GPR associations were identified using databases of the organism and specific literature.

(v) Definition of constraints: some constraints were defined in the model in this manual curation step, including stoichiometric and thermodynamic constraints (through the reversibility and irreversibility of fluxes).

2.4. Metabolic Model Design

The metabolic model designed and refined following the manual curation step was converted into a mathematical representation, termed a stoichiometric matrix, which encouraged the development of a wide variety of computational tools to analyze network properties.

The constraints of capacity, which are the upper and lower limits defining the maximum and minimum fluxes allowed for the reactions, were added in this step. The inputs of the stoichiometric matrix are the coefficients of the metabolites in the reactions with negative values for consumed metabolites (substrates) and positive values when the metabolites are produced or secreted (products) (Additional File 1 see Supplementary Material available online at http://dx.doi.org/10.1155/2016/7863706).

2.5. Metabolic Pathway Validation

The computational model sought to examine the metabolic capabilities and to evaluate the system properties they may perform under the constraints imposed on the cell. Thus, the final step in the reconstruction process consisted of assessing, evaluating, and validating the fatty acid biosynthesis pathway of E. antarcticum B7. The validation of that metabolic model was performed using simulation and flux analysis. The fatty acid biosynthesis pathway is well described in the literature because it is a highly conserved process among organisms, which facilitated its complete definition. Thus, most gaps had already been filled during the manual curation process.

3. Results and Discussion

The expectation to understand the relationship between the genome and the physiology of a particular organism was a key incentive for reconstructing metabolic networks. Protocol adaptations using semiautomatic and automatic methods are necessary to reconstruct the metabolic networks of organisms with few reported data on their metabolic capabilities, including E. antarcticum B7.

3.1. Pathway Reconstruction Using the Semiautomatic Method

The draft of the metabolic network of E. antarcticum B7 was retrieved from the KEGG database [26]. The KEGG Metabolic Pathway tool was used to highlight the fatty acid biosynthesis pathway from the resulting draft of the metabolic network. The genes annotated and identified using KEGG and their respective enzymes are shown in Figure 1(a) in green, and the others are listed in white boxes. Table 1 shows the genes, locus tags, and enzymes identified using the KEGG Metabolic Pathway tool. A total of 11 locus tags associated with their respective genes and Enzyme Commission (EC) numbers were identified; only the locus tag Eab7_2235 has no added gene associated with it.

Table 1: Genes, locus tags, and EC numbers identified in the draft of the fatty acid biosynthesis pathway of E. antarcticum B7. The features displayed were generated using the methods: semiautomatic and automatic.

Semiautomatic method
Automatic method

Gene
Locus tag
EC number
Gene
Locus tag
EC number

accA
Eab7_2059
6.4.1.2
accA
Eab7_2059
6.4.1.2

accB
Eab7_0870
6.4.1.2
accB
Eab7_0870
6.4.1.14

accC
Eab7_0871
6.4.1.2
accC
Eab7_0871
6.4.1.2

accD
Eab7_2060
6.4.1.2
accD
Eab7_2060
6.4.1.14

fabD
Eab7_1760
2.3.1.39
fabD
Eab7_1760
2.3.1.39

fabH1
Eab7_1911
2.3.1.180
fabH1
Eab7_1911
2.3.1.180

fabF
Eab7_1910
2.3.1.179
fabF
Eab7_1910
2.3.1.179

fabG
Eab7_1795
1.1.1.100
fabG
Eab7_1795
1.1.1.100

fabZ
Eab7_2463
4.2.1.59
fabI
Eab7_1885
1.3.1.10

fabI
Eab7_1885
1.3.1.10

fabK
Eab7_0377
1.3.1.9

—
Eab7_2235
1.14.19.2

(a)

(b)

Figure 1: Drafts of the fatty acid biosynthesis pathway of E. antarcticum B7 bacteria. The drafts were designed using the following methods: (a) semiautomatic method, generated by the KEGG database, and (b) automatic method, generated by online tool RAST. Colored boxes indicate the possibility of the presence of enzymes in the pathway.

3.2. Pathway Reconstruction Using the Automatic Method

The RAST/SEED tool does not provide graphic display of the metabolic map draft as KEGG; for this purpose, it uses a standard table to list the 247 metabolic pathways that compose the network, regardless of whether they were identified in the genome of the microorganism. RAST identified 25 reactions, 40 compounds, and 20 EC numbers in the fatty acid biosynthesis pathway of E. antarcticum B7. Table 1 outlines the genes, locus tags, and enzymes identified at this step. Annotated genes identified by RAST and their respective enzymes are colored in purple, and the others are listed in white boxes (Figure 1(b)).

The analysis of both fatty acid biosynthetic pathway drafts shows that the draft generated using KEGG apparently has the most complete flux, except for enzyme 6.3.4.14, which is exclusively present in the draft resulting from the RAST tool. The draft generated using RAST has a gap in which the enzymes FabA and FabB are not included in the pathway elongation process. The flux for the production of hexadecenoic acid is also absent from the pathway generated using RAST.

Artemis software was used to confirm the presence of all genes selected through the automatic and semiautomatic methods. The genes accABCD, fabD, fabH1, fabF, fabG, fabZ, and fabI and the locus tag Eab7_2235 were described in the genome of E. antarcticum B7, except fabK gene, which was detected only by the automatic method.

The KGML file produced by KEGG was submitted to the software KEGGtranslator [33] to be converted into a SBML (System Biology Markup Language) file [34]. This file was converted into an Excel spreadsheet using MATLAB functions. The files in SBML format and the Excel spreadsheets are the most used formats in metabolic reconstructions. The reactions and metabolites of the preliminary network generated using KEGG could be visualized in the spreadsheet.

The data generated using the RAST/SEED tool were analyzed and added to the first step of the process, supplementing the data collected using KEGG. The files generated with both platforms were used to manage the manual curation data.

The larger number of genes identified using KEGG (12) compared to those found using RAST/SEED (9) may be explained because the former uses orthology (KEGG Orthology (KO)) through protein homology to identify the so-called metabolite genes [35] in a genome, which facilitates finding gene-protein-reaction (GPR) associations.

3.3. Manual Curation of the Metabolic Pathway of De Novo Fatty Acid Synthesis

The reactions of the fatty acid synthesis pathway were annotated and refined. The metabolites were organized into two compartments (cytoplasm and extracellular compartment) based on the location of the enzymes associated with each pathway. The cofactors and the reversibility of the reactions were compiled from the data published in the literature and online tools (ENZYME and BRENDA). The EC number was noted, and the genes were identified. A summary of those results is shown in Table 2. Thermodynamic analysis of the reactions revealed that malonyl-CoA synthesis from acetyl-CoA (AcCoA) is an irreversible process; similar to the process regarding the fabF gene, the Eab7_2235 locus tag, and the extracellular metabolites the other processes are reversible.

Table 2: Relationships between components of the E. antarcticum B7 fatty acid biosynthesis pathway following curation. The signals and indicate irreversible and reversible reactions, respectively.

Gene
Locus tag
EC number
Enzyme
Reaction
Fold change ()

accA
Eab7_2059
6.4.1.2
Acetyl-CoA carboxylase carboxyl transferase alpha subunit
ATP + acetyl-CoA + ADP + orthophosphate + malonyl-CoA
0.4562

accB
Eab7_0870
6.4.1.2
Acetyl-CoA carboxylase biotin-carboxyl carrier protein
ATP + acetyl-CoA + ADP + orthophosphate + malonyl-CoA
−0.05773

accC
Eab7_0871
6.4.1.2
6.4.1.14
Acetyl-CoA carboxylase, biotin carboxylase subunit
ATP + acetyl-CoA + ADP + orthophosphate + malonyl-CoA
−0.5623

accD
Eab7_2060
6.4.1.2
Acetyl-CoA carboxylase carboxyl transferase beta subunit
ATP + acetyl-CoA + ADP + orthophosphate + malonyl-CoA
−0.2811

fabD
Eab7_1760
2.3.1.39
ACP S-malonyl transferase
Malonyl-CoA + ACP CoA + malonyl-(acp)
0.9448

fabH1
Eab7_1911
2.3.1.180
3-Oxoacyl-ACP synthase III
Acetyl-CoA + malonyl-(acp) acetoacetyl-(acp) + CoA + CO₂
0.8512

fabF
Eab7_1910
2.3.1.179
3-Oxoacyl-ACP synthase II
Acetyl-(acp) + malonyl-(acp) acetoacetyl-(acp) + CO₂ + ACP
0.6942

fabG
Eab7_1895
1.1.1.100
3-Oxoacyl-ACP reductase
Acetoacetyl-(acp) + NADPH + H⁺   (R)-3-hydroxybutanoyl-(acp) + NADP⁺
0.8523

fabZ
Eab7_2463
4.2.1.59
3-Hydroxyacyl-ACP dehydratase
(R)-3-hydroxybutanoyl-(acp) but-2-enoyl-(acp) + H₂O
0.12902

fabI
Eab7_1885
1.3.1.9
1.3.1.10
Enoyl-ACP reductase I
But-2-enoyl-(acp) + NADH + H⁺   butyryl-(acp) + NAD⁺
0.2969

—
Eab7_2235
1.14.19.2
Acyl-ACP desaturase
Hexadecanoyl-(acp) + acceptor_reduced + O₂   hexadecenoyl-(acp) + acceptor + 2H₂O
1.3768

It is very important to assess the quality of the annotated genome submitted to the online tools during curation. The literature categorically states that the quality of the reconstructed network directly depends on the annotated genome of the organism. The rule is to use the latest updated version of the annotated genome [36–38].

3.4. Curated Metabolic Network

The previous step added constraints regarding the stoichiometry (chemical balance), thermodynamics (reversibility of reactions), and physiology (cofactors used) of the study model. The result was a system of equations that describes the cell metabolism according to the metabolites of interest. The mathematical representation of this model essentially consisted of describing the performance of the fatty acid biosynthesis pathway using a stoichiometric matrix. This data structure consists of 54 metabolites and 59 reactions, resulting in a 54 × 59 stoichiometric matrix (Additional File 1). Additional File 2 shows the list of biochemical reactions identified in the curated model.

The reconstructed pathway model was converted into SBML in the used MATLAB toolbox. The SBML file was validated using the tool SBML validator and was then submitted to the tool Cytoscape, which generated the network of Figure 2. The gene-protein-reaction (GPR) representation therein describes the degree of connectivity of each enzyme in the pathway. Vertices with few connections are in green, the vertices with regular numbers of connections are in yellow, and the vertices with large numbers of connections are in red. The network connectivity obeys a scale-free model [39].

Figure 2: Layout of the fatty acid biosynthesis pathway generated using Cytoscape. Green vertices: fewer connections. Yellow vertices: regular number of connections. Red vertices: large number of connections. This network is a free model scale.

3.5. Flux Balance Analysis (FBA)

The FBA was coded in MATLAB implementing a constrained linear program using the GLPK (GNU Linear Programming Kit) linear optimization library [8]. All fluxes were calculated in percentage of the input flux of AcCoA (reaction 39), which was fixed to 100. Hexadecenoic acid is the key metabolic product; thus the respective flux (reaction 41) was set as the maximization target for FBA. To improve convergence, upper and lower bounds were for irreversible reactions and for reversible reaction. The final optimized fluxes are shown in Additional File 3. The target maximum, reaction R41 in Additional File 3, was 7.69, which may be read as 7.69 moles of hexadecenoic acid produced for every 100 moles of AcCoA consumed per unit time per unit cell mass.

Figure 3(a) shows the generated flux plot, which shows the variation occurring between the response fluxes, with the majority, approximately 37, showing positive values smaller than 20, while 15 are above that range.

(a)

(b)

Figure 3: Flux graph generated using the software MATLAB, which uses methods based on FBA constraints. The vertical axis represents the fluxes calculated from the input of 100 moles of acetyl-CoA (AcCoA). The horizontal axis represents the reactions participating in the fatty acid biosynthesis pathway. The blue bars determine the output percentage for each pathway reaction. (a) FBA for maximizing the production of hexadecenoic acid. (b) FBA calculated with the values obtained from the transcriptome.

3.6. Influence of Transcriptome in FBA Calculation

Log base 2-fold change values () obtained in vitro by Dall’Agnol and colleagues [22] were used to evaluate the influence of differential expression in the FBA calculation. These values were obtained by comparison of RPKM (reads per kilobase per million reads sequenced) generated in the transcriptome of the bacterium at 0°C and 37°C. The of genes that composes the fatty acid biosynthesis pathway is shown in Table 2.

As presented by Dall’Agnol and colleagues [22], the aerobic energetic metabolism of E. antarcticum B7 at 0°C is repressed, and a fraction of the acetyl-CoA is probably used as a substrate to synthesize short-chain fatty acids in cold. The synthesis begins with the conversion of acetyl-CoA to malonyl-CoA catalyzed by the multimeric enzyme encoded by the genes accABCD. In FBA analysis, only half of the input of acetyl-CoA (100 mM) is converted to malonyl-CoA, which binds to the acyl carrier protein (ACP) at 0°C (Figure 3(b)). As stated earlier, the other half of acetyl-CoA is probably used for energy generation.

The remaining route is cyclical, being the reactions catalyzed by enzymes encoded by fabF, fabG, fabZ, and fabI genes. At each round, two carbons are added to the growing chain of fatty acid. In these reactions, the flow of metabolites remains unchanged until the fatty acid molecule reaches a size of eight carbons (octadecanoyl-ACP in reaction 38) where the percentage of the flow amount decreases (Figure 3(b)). These results are consistent with the previously published data which affirm that bacteria decrease their fatty acid chains to survive in cold. These short fatty acid molecules will be converted into membrane phospholipids in order to maintain the fluidity of this biological barrier.

3.7. Elementary Flux Modes

The calculation of elementary modes was performed in MATLAB using the Metatool toolbox [40] (modo_elementar.m code in the Supplemental Information). A total of 13 elementary flux modes were found for the fatty acid biosynthesis pathway of E. antarcticum B7 (Additional File 4). Of these, only 4 elementary modes (2, 5, 8, and 11) have a positive nonzero coefficient for reaction 41, which indicates that the target product hexadecenoic acid may only be generated by one of these four elementary modes. The routes identified in EM2 and EM5 begin at the second reaction (R2) which is catalyzed by the enzyme FabD. The value of R2 for both elementary modes indicates a considerable production of malonyl-CoA. The remaining reactions of EM2 and EM5 are presented in Additional File 4. In EM8 and EM11, the routes begin from R1. The value 1 for this reaction indicates a lower activity, reflecting in a lower production of hexadecenoic acid. Regarding reaction 41, the values of elementary modes 2 and 5, in this case 1, are higher than the values of 8 and 11 (0.142857 each), indicating that both elementary modes 2 and 5 produce 1 mol hexadecenoic acid when they are active, while elementary modes 8 and 11 produce 0.142857 moles.

4. Conclusions

The first metabolic pathway of E. antarcticum B7, reconstructed following the steps defined in this work, suggests that the protocol used is a suitable tool for further metabolic reconstruction studies. Almost all the first steps of the process were automated; however, manual curation was, as usual, laborious because it required an intensive search for available data.

The metabolic pathway of fatty acid biosynthesis was representative and consistent under the limits and boundary conditions set. The FBA and elementary mode methods were used to examine the hexadecenoic acid production potential of the reconstructed pathway. The application of constraint-based modeling revealed being very useful to assess network operation plasticity, even if the intracellular kinetics are unknown. The in silico analysis performed using FBA enabled a quantitative assessment of hexadecenoic acid production potential.

Finally, a key issue involves deciding when to stop the process and to consider the reconstruction finalized, at least temporarily. This decision is usually based on the reconstruction purpose. The most complete metabolic model currently available is the E. coli model, which has been researched and refined for over 10 years [41–45]. Other studies constantly updating their models are Homo sapiens, with three reconstructions [46–48], and S. cerevisiae, with more than a dozen reconstructions, including two in 2013 [49, 50]. The protocol reported in the present study may be used to compile several data pieces available in the literature aimed at proposing possible metabolic pathways, thereby enabling deeper research of the metabolism under study.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This study was supported by the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES).

References

J. W. Lee, D. Na, J. M. Park, J. Lee, S. Choi, and S. Y. Lee, “Systems metabolic engineering of microorganisms for natural and non-natural chemicals,” Nature Chemical Biology, vol. 8, no. 6, pp. 536–546, 2012.

J. Becker and C. Wittmann, “Systems and synthetic metabolic engineering for amino acid production—the heartbeat of industrial strain development,” Current Opinion in Biotechnology, vol. 23, no. 5, pp. 718–726, 2012.

C.-H. Chou, W.-C. Chang, C.-M. Chiu, C.-C. Huang, and H.-D. Huang, “FMM: a web server for metabolic pathway reconstruction and comparative analysis,” Nucleic Acids Research, vol. 37, no. 2, pp. W129–W134, 2009.

P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software Environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, pp. 2498–2504, 2003.

A. Funahashi, Y. Matsuoka, A. Jouraku, M. Morohashi, N. Kikuchi, and H. Kitano, “CellDesigner 3.5: a versatile modeling tool for biochemical networks,” Proceedings of the IEEE, vol. 96, no. 8, pp. 1254–1265, 2008.

F. T. Bergmann and H. M. Sauro, “SBW—a modular framework for systems biology,” in Proceedings of the 37th Winter Simulation Conference (WSC '06), pp. 1637–1645, December 2006.

S. Hoops, S. Sahle, R. Gauges et al., “COPASI—a COmplex PAthway SImulator,” Bioinformatics, vol. 22, no. 24, pp. 3067–3074, 2006.

J. Schellenberger, R. Que, R. M. T. Fleming et al., “Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0,” Nature Protocols, vol. 6, no. 9, pp. 1290–1307, 2011.

M. Terzer, N. D. Maynard, M. W. Covert, and J. Stelling, “Genome-scale metabolic networks,” Wiley Interdisciplinary Reviews: Systems Biology and Medicine, vol. 1, no. 3, pp. 285–297, 2009.

D. McCloskey, B. Ø. Palsson, and A. M. Feist, “Basic and applied uses of genome-scale metabolic network reconstructions of Escherichia coli,” Molecular Systems Biology, vol. 9, article 661, 2013.

M. Fondi, I. Maida, E. Perrin et al., “Genome-scale metabolic reconstruction and constraint-based modelling of the Antarctic bacterium Pseudoalteromonas haloplanktis TAC125,” Environmental Microbiology, vol. 17, no. 3, pp. 751–766, 2015.

D. McCloskey, B. Ø. Palsson, and A. M. Feist, “Basic and applied uses of genome‐scale metabolic network reconstructions of Escherichia coli,” Molecular Systems Biology, vol. 9, article 661, 2013.

A. M. Feist, M. J. Herrgård, I. Thiele, J. L. Reed, and B. Ø. Palsson, “Reconstruction of biochemical networks in microorganisms,” Nature Reviews Microbiology, vol. 7, no. 2, pp. 129–143, 2009.

J. L. Reed, I. Famili, I. Thiele, and B. Ø. Palsson, “Towards multidimensional genome annotation,” Nature Reviews Genetics, vol. 7, no. 2, pp. 130–141, 2006.

I. Thiele and B. Ø. Palsson, “A protocol for generating a high-quality genome-scale metabolic reconstruction,” Nature Protocols, vol. 5, no. 1, pp. 93–121, 2010.

J. S. Edwards, R. U. Ibarra, and B. O. Palsson, “In silico predictions of Escherichia coli metabolic capabilities are consistent with experimental data,” Nature Biotechnology, vol. 19, no. 2, pp. 125–130, 2001.

R. Mahadevan, B. Ø. Palsson, and D. R. Lovley, “In situ to in silico and back: elucidating the physiology and ecology of Geobacter spp. using genome-scale modelling,” Nature Reviews Microbiology, vol. 9, no. 1, pp. 39–50, 2011.

J. D. Trawick and C. H. Schilling, “Use of constraint-based modeling for the prediction and validation of antimicrobial targets,” Biochemical Pharmacology, vol. 71, no. 7, pp. 1026–1035, 2006.

R. Y. Morita, “Psychrophilic bacteria,” Bacteriological Reviews, vol. 39, no. 2, pp. 144–167, 1975.

R. Mulvaney, N. J. Abram, R. C. A. Hindmarsh et al., “Recent Antarctic Peninsula warming relative to Holocene climate and ice-shelf history,” Nature, vol. 489, no. 7414, pp. 141–144, 2012.

A. R. Carneiro, R. T. J. Ramos, H. Dall'Agnol et al., “Genome sequence of Exiguobacterium antarcticum B7, isolated from a biofilm in Ginger Lake, King George Island, Antarctica,” Journal of Bacteriology, vol. 194, no. 23, pp. 6689–6690, 2012.

H. Dall'Agnol, R. A. Baraúna, P. de Sá et al., “Omics profiles used to evaluate the gene expression of Exiguobacterium antarcticum B7 during cold adaptation,” BMC Genomics, vol. 15, no. 1, article 986, 2014.

A. Casanueva, M. Tuffin, C. Cary, and D. A. Cowan, “Molecular adaptations to psychrophily: the impact of ‘omic’ technologies,” Trends in Microbiology, vol. 18, no. 8, pp. 374–381, 2010.

G. E. Schujman, L. Paoletti, A. D. Grossman, and D. de Mendoza, “FapR, a bacterial transcription factor involved in global regulation of membrane lipid biosynthesis,” Developmental Cell, vol. 4, no. 5, pp. 663–672, 2003.

E. Pitkänen, J. Rousu, and E. Ukkonen, “Computational methods for metabolic reconstruction,” Current Opinion in Biotechnology, vol. 21, no. 1, pp. 70–77, 2010.

H. Ogata, S. Goto, K. Sato, W. Fujibuchi, H. Bono, and M. Kanehisa, “KEGG: kyoto encyclopedia of genes and genomes,” Nucleic Acids Research, vol. 27, no. 1, pp. 29–34, 1999.

T. Carver, S. R. Harris, M. Berriman, J. Parkhill, and J. A. McQuillan, “Artemis: an integrated platform for visualization and analysis of high-throughput sequence-based experimental data,” Bioinformatics, vol. 28, no. 4, pp. 464–469, 2012.

Y. Moriya, M. Itoh, S. Okuda, A. C. Yoshizawa, and M. Kanehisa, “KAAS: an automatic genome annotation and pathway reconstruction server,” Nucleic Acids Research, vol. 35, no. 2, pp. W182–W185, 2007.

R. K. Aziz, D. Bartels, A. Best et al., “The RAST Server: rapid annotations using subsystems technology,” BMC Genomics, vol. 9, article 75, 2008.

A. Bairoch, “The ENZYME database in 2000,” Nucleic Acids Research, vol. 28, no. 1, pp. 304–305, 2000.

R. Overbeek, R. Olson, G. D. Pusch et al., “The SEED and the rapid annotation of microbial genomes using Subsystems Technology (RAST),” Nucleic Acids Research, vol. 42, no. 1, pp. D206–D214, 2014.

A. Flamholz, E. Noor, A. Bar-Even, and R. Milo, “EQuilibrator—the biochemical thermodynamics calculator,” Nucleic Acids Research, vol. 40, no. 1, pp. D770–D775, 2012.

C. Wrzodek, A. Dräger, and A. Zell, “KEGGtranslator: visualizing and converting the KEGG PATHWAY database to various formats,” Bioinformatics, vol. 27, no. 16, pp. 2314–2315, 2011.

M. Hucka, A. Finney, H. M. Sauro et al., “The systems biology markup language (SBML): a medium for representation and exchange of biochemical network models,” Bioinformatics, vol. 19, no. 4, pp. 524–531, 2003.

M. Ashburner, C. A. Ball, J. A. Blake et al., “Gene ontology: tool for the unification of biology,” Nature Genetics, vol. 25, no. 1, pp. 25–29, 2000.

J. J. Hamilton and J. L. Reed, “Software platforms to facilitate reconstructing genome-scale metabolic networks,” Environmental Microbiology, vol. 16, no. 1, pp. 49–59, 2014.

G. J. E. Baart and D. E. Martens, “Genome-scale metabolic models: reconstruction and analysis,” Methods in Molecular Biology, vol. 799, pp. 107–126, 2012.

V. Lacroix, L. Cottret, P. Thébault, and M.-F. Sagot, “An introduction to metabolic networks and their structural analysis,” IEEE/ACM Transactions on Computational Biology and Bioinformatics, vol. 5, no. 4, pp. 594–617, 2008.

A.-L. Barabási and Z. N. Öltvai, “Network biology: understanding the cell's functional organization,” Nature Reviews Genetics, vol. 5, no. 2, pp. 101–113, 2004.

T. Pfeiffer, I. Sánchez-Valdenebro, J. C. Nuño, F. Montero, and S. Schuster, “METATOOL: for studying metabolic networks,” Bioinformatics, vol. 15, no. 3, pp. 251–257, 1999.

E. Almaas, B. Kovács, T. Vicsek, Z. N. Oltvai, and A.-L. Barabási, “Global organization of metabolic fluxes in the bacterium Escherichia coli,” Nature, vol. 427, no. 6977, pp. 839–843, 2004.

J. S. Edwards and B. Ø. Palsson, “The Escherichia coli MG1655 in silico metabolic genotype: its definition, characteristics, and capabilities,” Proceedings of the National Academy of Sciences of the United States of America, vol. 97, no. 10, pp. 5528–5533, 2000.

I. M. Keseler, J. Collado-Vides, S. Gama-Castro et al., “EcoCyc: a comprehensive database resource for Escherichia coli,” Nucleic Acids Research, vol. 33, pp. D334–D337, 2005.

J. D. Orth, T. M. Conrad, J. Na et al., “A comprehensive genome-scale reconstruction of Escherichia coli metabolism-2011,” Molecular Systems Biology, vol. 7, article 535, 2011.

A. M. Feist and B. Ø. Palsson, “The growing scope of applications of genome-scale metabolic reconstructions using Escherichia coli,” Nature Biotechnology, vol. 26, no. 6, pp. 659–667, 2008.

N. C. Duarte, S. A. Becker, N. Jamshidi et al., “Global reconstruction of the human metabolic network based on genomic and bibliomic data,” Proceedings of the National Academy of Sciences of the United States of America, vol. 104, no. 6, pp. 1777–1782, 2007.

T. Y. Kim, S. B. Sohn, Y. B. Kim, W. J. Kim, and S. Y. Lee, “Recent advances in reconstruction and applications of genome-scale metabolic models,” Current Opinion in Biotechnology, vol. 23, no. 4, pp. 617–623, 2012.

I. Thiele, N. Swainston, R. M. Fleming et al., “A community-driven global reconstruction of human metabolism,” Nature Biotechnology, vol. 31, no. 5, pp. 419–425, 2013.

B. D. Heavner, K. Smallbone, N. D. Price, and L. P. Walker, “Version 6 of the consensus yeast metabolic network refines biochemical coverage and improves model performance,” Database, vol. 2013, Article ID bat059, 2013.

T. Österlund, I. Nookaew, S. Bordel, and J. Nielsen, “Mapping condition-dependent regulation of metabolism in yeast through genome-scale modeling,” BMC Systems Biology, vol. 7, article 36, 2013.