- Research article
- Open Access
Reconstruction of a genome-scale metabolic model for Actinobacillus succinogenes 130Z
BMC Systems Biology volume 12, Article number: 61 (2018)
Actinobacillus succinogenes is a promising bacterial catalyst for the bioproduction of succinic acid from low-cost raw materials. In this work, a genome-scale metabolic model was reconstructed and used to assess the metabolic capabilities of this microorganism under producing conditions.
The model, iBP722, was reconstructed based on the functional reannotation of the complete genome sequence of A. succinogenes 130Z and manual inspection of metabolic pathways, covering 1072 enzymatic reactions associated with 722 metabolic genes that involve 713 metabolites. The highly curated model was effective in capturing the growth of A. succinogenes on various carbon sources, as well as the SA production under various growth conditions with fair agreement between experimental and predicted data. Calculated flux distributions under different conditions show that a number of metabolic pathways are affected by the activity of some metabolic enzymes at key nodes in metabolism, including the transport mechanism of carbon sources and the ability to fix carbon dioxide.
The established genome-scale metabolic model can be used for model-driven strain design and medium alteration to improve succinic acid yields.
Actinobacillus succinogenes is a gram-negative facultative anaerobic bacterium and is one of the major natural producers of succinic acid (SA). It can grow on a broad range of substrates, including arabinose, cellobiose, fructose, galactose, glucose, lactose, maltose, mannitol, mannose, sucrose and xylose, producing a mixture of by-products (e.g., SA, formic acid (FA), acetic acid (AA), and ethanol (EtOH)) as the main by-products [1, 2]. Its tolerance to high sugar concentrations (up to 160 g.L− 1 of glucose ) and high levels of organic acids , as well as its capnophilic nature , make this microorganism potentially interesting for the production of SA at the industrial scale. High-titer succinate production using low-cost feedstocks like cane molasses or corn straw has been obtained [2, 6,7,8,9,10,11,12,13,14,15]; however, significant amounts of other organic acids were also produced, increasing the downstream processing costs, which makes this bioprocess less competitive. Currently, SA is mainly produced from petrochemical feedstocks through the hydrogenation of maleic acid or maleic anhydride . However, the bio-based production using low pH yeast fermentation [17,18,19] or anaerobic fermentation using bacteria [20,21,22,23] has been successfully implemented by companies like Myriant , BASF  or BioAmber , offering economically and ecologically attractive alternatives to the conventional petro-based SA production [27,28,29]. Some examples of SA producing systems are given in Table 1, including naturally-producing bacteria like Basfia succiniciproducens and Mannheimia succiniproducens and genetically engineered organisms such as E. coli or S. cerevisiae. So far, natural producers appear to outperform most engineered strains, but developments in strain design and fermentative processes are expected to promote the production of bio-based SA by metabolically engineered microorganisms. For instance, M. succiniproducens has been metabolically engineered by removing competing pathways, resulting in an increase in the SA yield from 0.45 to 0.76 g of SA per g of glucose .
Under optimized conditions, the A. succinogenes wild-type strain is able to produce up to 98 g.L− 1 of SA with an approximate yield of 90% (w/w) on glucose . The optimization of bioprocesses has proven to further increase SA production by using high concentrations of carbon dioxide (CO2) and/or hydrogen (H2) [32, 33]. Other studies have shown that the redox state of the fermentation broth affects SA production, which can be improved by manipulating the supplementation of oxidant and reducing agents [34, 35]. Yet, due to the accumulation of other fermentative by-products, SA yields are still far below the maximum theoretical yield of 1.12 g.g− 1 of glucose consumed (YSA/Glc) . A comprehensive understanding of the metabolism and the phenotypic responses to environmental perturbations is a major step for developing efficient bioprocesses for SA production.
Genome-scale metabolic models (GSMMs) have proven to be powerful tools for understanding and re-designing the metabolism of microbial strains. For instance, the optimization of SA production in E. coli or S. cerevisiae has been achieved by applying metabolic engineering strategies supported by in silico modelling of metabolic networks [37,38,39,40]. There are three main pathways for SA biosynthesis, including the tricarboxylic acid (TCA) cycle in the oxidative direction, the glyoxylate shunt and the reductive TCA pathway . Typically, under aerobic conditions, either the oxidative TCA cycle or the glyoxylate shunt can be used for SA production and several studies have used both S. cerevisiae and E. coli to exploit these metabolic pathways [38, 40]. The redirection of the carbon flux through the glyoxylate shunt provides some advantages over the oxidative TCA cycle, mainly because the decarboxylation of isocitrate to succinyl-CoA leads to carbon loss  and, in the case of S. cerevisiae, the SA channelling from the mitochondria to cytosol is avoided as TCA cycle enzymes are located in mitochondria. However, if a reductive TCA pathway is used, a 2-fold maximum theoretical yield (2 mol.mol− 1 of glucose) can be achieved compared to the oxidative route (1 mol.mol− 1 of glucose) . Many organisms, such as E. coli and S. cerevisiae have been tested for SA production under anaerobic conditions using the reductive branch of the TCA cycle [38, 42]. However, reducing power limitations (i.e. NADH) or redox balance issues have shown to have an impact on the final SA yields . Metabolic engineering strategies driven by in silico modelling may allow to overcome these disadvantages, both under aerobic and anaerobic conditions. Strategies for increasing energy and/or cofactor pools [43,44,45], to overcome enzyme limitations  or to decrease by-product generation  are just a few examples that have been exploited .
A. succinogenes produces SA anaerobically through the reductive branch of the TCA cycle (i.e., C4 pathway) using fumarate as the final electron acceptor, which makes this metabolic branch highly dependent on the redox state of cultures. Phosphoenolpyruvate (PEP) node controls the amount of flux that is directed towards the C4 and C3 pathways, adjusting the level of fermentative products generated by each pathway. Other metabolic nodes like oxaloacetate (OAA) and malate (MAL) have shown to link the C4 and C3 pathways via decarboxylating enzymes, i.e. oxaloacetate decarboxylase and the NADPH-dependent malic enzyme, respectively. The split of carbon flux at these branching points is largely influenced by various factors, such as the availability of CO2, pH or carbon sources. It has been previously shown that an increase in the concentration of dissolved CO2 in the fermentation broth through the supplementation of magnesium carbonate  or sodium bicarbonate  can promote an increase in the carbon flow toward the C4 pathway, thus increasing the production of SA. The presence of carbonic anhydrases (putative coding gene, Asuc_1199), which interconverts CO2 and bicarbonate (HCO3−), may also contribute to increase CO2 fixation [12, 49].
Here, a GSMM (named iBP722) representing a wide range of metabolic capabilities of A. succinogenes 130Z is proposed. The model allows predicting and analysing the impact of stoichiometric and physiological constraints known to apply at steady-state conditions. Although the central carbon metabolism of A. succinogenes 130Z has been comprehensively described [32, 33, 50, 51], the overall representation of metabolic pathways associated with various catabolic and anabolic capabilities of the organism are now made available. The biosynthetic pathways for vitamins, cofactors and other biomass building blocks are described, as well as respiratory and energy consuming assimilatory pathway. This model provides a detailed insight on the metabolism of A. succinogenes that can be systematically explored to improve the bioproduction of SA.
Metabolic functional annotation
The complete genome sequence of A. succinogenes 130Z (GenBank accession number NC_009655)  was used for the functional annotation of genes based on homology searching methods. The annotated genes with potential metabolic roles were manually inspected and associated with the corresponding coding enzyme(s) and biochemical reaction(s). An internally developed platform was used to compute, assign and curate gene metabolic functions. This platform couples automated annotation tools with manual curation procedures that allows the assignment of metabolic functions to coding sequences (CDSs) of a particular genome. The pipeline consists in five main steps (see Fig. 1):
The application of homology search tools like BLAST  and HMMER  against sequence databases, such as UniProt  to find the best alignment between sequences. A ranked list of hits with the most significant matches to each query is obtained with the respective information, including scores and protein features that may contain Enzyme Commission (EC) numbers.
The computation of “functional scores” based on BLAST and HMMER scores. As the range of enzymatic functions attributed to each CDS can vary from tool to tool, a second score was computed, i.e. so-called functional score, in order to propose the best candidate metabolic functions. Besides BLAST and HMMER scores, a functional score was computed similarly to the method used by Merlin , which is based on the frequency of the EC number in the homology hits and on the taxonomy distance between the target strain and other strains within hit results. This weighed functional scores range between 0 and 1 (1 corresponding to a high confidence score).
Assignment of putative metabolic function(s). EC numbers associated with the highest functional score are automatically attributed to each CDS, as well as the corresponding metabolic reaction(s) from an internal reactions database.
The assignment of metabolic functions to each CDS is manually revised, as automatic assignments may fail when more than one high scoring EC number is found and/or EC number(s) are incorrectly associated in databases. Therefore the user is allowed to inspect all putative assignments and select the most appropriate or modify the functional assignment, which can be based on previous knowledge or other assumptions defined by the user. Furthermore, more than one CDS can be associated to the same EC number, which in many cases consists of subunits of the same multimeric enzyme. The classification of multimeric or monomeric subunits is also defined at this stage using an internal database.
After validation of functional assignments and enzyme subunits, the association of metabolic reactions is carried out based on the associated EC number(s) and/or previous knowledge using as a reference an internal reactions database. This step is perhaps the most critical during the reconstruction of the metabolic network, as it will define the set of stoichiometric reactions that characterize a specific organism.
To note that this annotation pipeline is flexible enough, such that unassigned CDSs in step (3) can be later reviewed, particularly during the gap filling process.
Construction of the metabolic network
The GSMM was initially constructed by compiling the annotated metabolic genes and their corresponding coding enzyme(s) and biochemical reaction(s). Additionally, spontaneous and transport reactions from databases like KEGG  and Transporter Classification Database (TCDB)  were added. Some metabolic reactions, although not associated with genes, were also included due to evidences found in literature. The preliminary draft model was then processed to:
Identify metabolic gaps (or missing reactions) that either consume or produce isolated (or dead-end) metabolites within the metabolic network. Then, each dead-end metabolite was inspected to search for metabolic reactions that consume, produce or transport this metabolite. Typically, MetaCyc or KEGG databases were used for gap-filling, i.e. to identify the sets of biochemical reactions that link each dead-end to a metabolite in the network. When several alternatives are found, a manual inspection is required and sequence-based homology searches using one or more amino acid sequences collected from potential candidates are used to find the most likely reactions set in the metabolic network of A. succinogenes 130Z.
Infer and correct the mass and charge balance of biochemical reactions. Stoichiometric coefficients of compounds in reactions are corrected, such that the reaction is balanced for mass and charge, usually by adding missing protons or water molecules.
Identify and correct the reversibility of reactions based on their thermodynamic properties and information found in literature. Databases like MetaCyc  and tools like eQuilibrator  were used.
Include a biomass reaction representing the basic macromolecular composition of A. succinogenes in terms of proteins, DNA, RNA, lipopolysaccharide (LPS), phospholipids, peptidoglycan, glycogen and cofactors and vitamins (CAV). The synthesis of each macromolecule was also represented by individual reactions considering the building blocks molar composition. For instance, the synthesis of one gram of protein was calculated based on the average amino acids composition using the set of encoded proteins of the A. succinogenes 130Z genome, according to the methodology proposed in . The synthesis of other cellular components, like CAV, was calculated assuming that each small molecule is equally present in one gram of CAV. Detailed information on the biomass composition can be found in Additional file 1.
Constraints-based flux analysis
Basic stoichiometric modelling methods, such as parsimonious flux balance analysis (pFBA) and flux variability analysis (FVA), were used to interrogate the metabolic properties and capabilities of the reconstructed GSMM for A. succinogenes 130Z under varying environmental conditions. Phenotype simulations were performed by maximizing the biomass reaction assuming growth under defined conditions, i.e. defined minimum media containing basic components required for biomass synthesis, such as vitamins, minerals and trace metals and explicit carbon, nitrogen and sulphur sources. Maximum theoretical product yields were calculated by maximizing the target product instead, ignoring the formation of biomass and ATP maintenance requirements, such that the costs of product biosynthesis in terms of carbon, energy and reducing equivalents were properly evaluated.
A. succinogenes genome-scale metabolic network
The construction of the GSMM for A. succinogenes 130Z was carried out in three different phases: (1) first, metabolic functions were assigned to genes; (2) then, biochemical reactions and enzymatic complexes (assigning proper gene-reaction relationships) were compiled to build a draft metabolic model, (3) which was thereafter completed and corrected by defining a biomass reaction, identifying network gaps and correcting inconsistencies when comparing with reported information.
The final iBP722 model consists of 722 unique genes (open reading frame (ORF) coverage − 35%), 1072 reactions and 713 unique metabolites. The model is available as a Systems Biology Markup Language (SBML) file at http://darwin.di.uminho.pt/models and BioModels database  assigned with the identifier MODEL1804130001 and detailed information on the curated metabolic network can be found in Additional file 2.
Predicted growth on different carbon sources
The iBP722 model was inspected for the ability to simulate the A. succinogenes growth on different conditions. Model simulations were performed using the OptFlux software  applying FBA-based methods that maximize the biomass reaction under defined conditions, i.e. exchange fluxes were constrained to specific values (usually experimentally measured fluxes) that allow testing growth under defined environmental conditions, such as sole carbon sources. Some exchange fluxes like those associated with the exchange of CO2, NH4+, Pi, H+, vitamins and trace elements were kept unconstrained (i.e. unlimited uptake rates) to provide unlimited basic nutrients for biomass synthesis.
A. succinogenes is auxotrophic for three amino acids: L-glutamate, L-cysteine and L-methionine. Glutamate auxotrophy is due to the organism’s inability to synthesize α-ketoglutarate, since the genes encoding for isocitrate dehydrogenase and α-ketoglutarate oxidoreductase enzymes are absent in the genome. A. succinogenes possesses most of the genes encoding enzymes associated with the cysteine biosynthetic pathway, but the absence of an adenylsulfate kinase to assimilate sulfate prevents the synthesis of hydrogen sulfide required for the synthesis of L-cysteine. This organism also lacks several genes for the biosynthesis of methionine as identified during the functional annotation process. Given that, exchange fluxes associated with amino acid auxotrophies in A. succinogenes, such as L-cysteine and L-methionine, were maintained unconstrained, except for L-glutamate that was limited to a minimum flux value to support growth requirements and avoid glutamate consumption as an additional carbon source. Further flux constraints were introduced when inspecting carbon flux distributions in the central carbon metabolism, and are further detailed in Table SI 11 (Additional file 2). For instance, succinyl-CoA ligase reaction was limited to a zero flux to avoid the formation of succinic acid from succinyl-CoA, in order to be consistent with in vivo observations . Assuming these flux constraints, growth predictions were computed and compared with in vivo observations [4, 27].
The model predicted accurately growth on more than 90% of the carbon sources, with only two not supporting in silico growth (FN = 2). In silico growth on beta-gentiobiose and D-arabitol was not predicted, as catabolic and transport reactions were not identified in A. succinogenes. Further information on transport activities included in the model is given in Additional file 4.
Predicted yields for fermentation products
The iBP722 model was further validated by predicting anaerobic production yields and comparing with experimental data from batch and/or chemostat cultures of A. succinogenes 130Z growing on glucose or xylose at different initial concentrations. [32, 51, 63] (Fig. 2). For each condition, carbon uptake rates were defined based on experimental values, except for condition C that was set to 8 mmol.gDCW− 1.h− 1, and then predicted production rates were used to calculate minimum and maximum FVA yields for biomass, SA, AA, FA and EtOH (Yx/S, YSA/S, YAA/S, YFA/S, YEtOH/S, respectively). FVA spans are given by the difference between the maximum and minimum predicted yields while maintaining 95% of the maximum biomass formation.
FVA predictions indicate higher yields for C4 by-products (i.e. SA) compared to C3 by-products (i.e. AA, EtOH and FA), which is in good agreement with experimental data; however predicted mass ratios between fermentative by-products, specifically SA/AA and SA/FA, are higher compared with in vivo observations, especially when considering maximum FVA ratios. Yet, FVA spans indicate a significant flexibility for these ratios, which might explain variations in the in vivo observations.
Improving model predictions
An FVA analysis was performed to elucidate these discrepancies and the metabolic flexibility associated with C3 and C4 metabolic pathways in A. succinogenes was examined. Given a set of flux constraints (e.g. substrate uptake rates, q S ), the predicted flux spans for the main reactions in the central carbon metabolism were calculated (Fig. 3a). Flux spans are given by the difference between the maximum and minimum predicted flux values of each metabolic reaction while maintaining 95% of the maximum biomass formation. As shown in Fig. 3b, the largest flux spans were associated with reactions around the PEP node (e.g. PEP carboxykinase (PPCK) and pyruvate kinase (PYK)), while reactions associated to the Embden-Meyerhoff-Parnas pathway (glucose-6-phosphate isomerase (PGI), phosphofructokinase (PFK) and enolase (ENO)) presented the lowest flux spans. Metabolic flux data estimated from 13C-labeling experiments  was further used to assess the accuracy of predicted flux spans. Most in vivo measurements were between estimated flux ranges, except for the CO2 exchange and formate dehydrogenase (FDHmq) reactions, indicating that the NADH-producing FDHmq reaction should be active and CO2 exchange flux should be lower. Changes in flux constraints associated with CO2 uptake (from unlimited to a maximum of 4 mmol.gDCW− 1.h− 1) altered predicted flux spans (Fig. 3c), especially for FDHmq, pyruvate-formate lyase (PFL) and the FA exchange reaction (EX_FA), indicating a higher flexibility in metabolic activities linked to FA accumulation.
The production of SA in A. succinogenes is influenced by several factors, namely the utilized carbon sources  or the availability of CO2 [5, 33, 48]. The iBP722 model was investigated for predicting the metabolic flexibility associated with the production of SA when changing CO2 availability or carbon sources, as well as the reversibility of metabolic reactions like malic enzyme (ME2) (Fig. 4). FVA yields for maximum growth simulations show that the production of reduced by-products, particularly EtOH, changes with the carbon source. As presented in Fig. 4a, the minimum and maximum FVA yields for ALCD2x under D-sorbitol growth conditions (0.62 and 0.69 mol.mol− 1) were higher compared to glucose (0.14 and 0.22 mol.mol− 1, respectively), with the consequent accumulation of higher amounts of EtOH. Under glucose conditions most of the carbon flux through the C3 branch is redirected toward the production of AA instead, via acetate kinase (ACKr) with the production of ATP (minimum and maximum FVA yields of 0.47 and 0.55 mol.mol− 1, respectively). Interestingly, however, is that predicted FVA yields for SA production hardly change between glucose and D-sorbitol growth conditions when considering the same transport mechanism, i.e. PEP:sugar phosphotransferase system (PEP:PTS). Yet, assuming that both transport mechanisms could be active (i.e., symport and PEP:PTS) for glucose uptake would increase SA minimum and maximum FVA yields from 0.63 and 0.67 to 0.82 and 0.96 mol.mol− 1, respectively.
The flux exchange between the C3 and C4 branches has been also investigated as a major factor affecting the metabolic flexibility of sugars fermentation in A. succinogenes . The decarboxylation of L-malate to pyruvate (PYR) reducing NADP+ to NADPH by ME2, may have a major role in the fermentative metabolism. Although the thermodynamics of this reaction is not conclusive regarding its reversibility, the metabolic network was tested using both the forward and the reverse directions. FVA results (Fig. 4b) show that SA production should increase when ME2 occurs in the reverse direction, especially for D-sorbitol conditions, redirecting most of the PYR pool toward the C4 branch. Consequently, a greater metabolic flexibility in SA production, especially under glucose growth conditions, was predicted.
The availability of CO2 was also shown to influence the production of SA (Fig. 4c). Model predictions indicated that maximum SA yields on glucose can decrease nearly 20% when decreasing the maximum CO2 uptake rates by 50%. A shift in carbon flux distributions is observed when decreasing CO2 uptake levels, redirecting most of the carbon flux toward the production of C3 by-products like FA and EtOH.
GSMMs are powerful tools to explore the metabolism of biological systems. In this work, the iBP722 model of A. succinogenes 130Z was reconstructed and used as a platform for the in silico analysis of the metabolic behaviour of this organism during anaerobic growth. The major end-product is SA, but significant amounts of other by-products such as EtOH, FA and AA are also accumulated. The possibility to predict and adjust the fermentative metabolism of A. succinogenes 130Z under different conditions brings new opportunities to exploit this host as a platform for the industrial production of SA and other reduced by-products.
In silico simulations were carried out using pFBA and FVA methods predicting growth behaviour under chemically defined medium. Model predictions were validated using reported physiological data and flux distributions from 13C experiments found in literature [32, 33, 51]. The iBP722 model supports growth predictions on 19 carbon sources under anaerobic conditions (Additional file 3), including C6 and C5 sugars that were found to enter the cell mainly through active transport systems (Additional file 4). The only two carbon sources with incorrect growth predictions were β-gentiobiose and D-arabitol, due to the absence of transport and catabolic pathways in the model, since no coding proteins were found. On the other hand, no-growth predictions for glycerol were correctly identified, though transport and catabolic reactions for glycerol consumption are present in the model. According to in vivo growth experiments , no cellular growth on glycerol as a sole carbon source is observed under anaerobic conditions, possibly caused by redox imbalance under anaerobiosis. However, the addition of external electron acceptors like dimethylsulfoxide (DMSO) has shown to recover cellular growth , which was also confirmed by in silico analysis (Additional file 5).
Further validations included the comparison of in vivo measurements and in silico FVA predictions for minimum and maximum yields under various conditions (Fig. 2). Despite being relatively variable (from 0.20 to 0.23 g.g− 1 under glucose conditions), experimental values for biomass yields were used to validate model predictions. Variations in experimental conditions, particularly associated with culture media that is often supplemented with yeast extract, interfering with carbon yields, or the initial sugar concentration that may affect bacterial growth due to substrate inhibition , may explain these differences. Nevertheless, biomass composition represented in the model is based on the work of McKinlay and co-workers , which is expected to provide accurate in silico predictions, especially regarding metabolic requirements to generate biomass contents per unit of substrate (i.e. biomass yields, Yx/s). Additional reaction constraints were included in the model to improve model predictions. For instance, flux constraints of L-glutamate uptake and succinyl-CoA synthase were changed to improve carbon-to-nitrogen ratios according to experimental measurements . Moreover, the energetic requirement for non-growth associated maintenance, i.e. the amount of ATP spent for cellular maintenance without growth, was adjusted to improve predicted biomass yields.
Fermentative products ratios were also compared, showing some consistency between experimental measurements and FVA predictions. Predicted mass ratios for SA/AA and SA/FA are higher than for SA/FA, which is consistent with in vivo data; but FVA spans for SA/AA and SA/FA suggest a huge flexibility in these ratios. This suggests that carbon flux distributions in the C3 branch are rather challenging to predict. In fact, FVA flux spans of fermentative pathways (Fig. 3b) indicate a high metabolic flexibility in the accumulation of by-products which, in some cases are inconsistent with experimental measurements. However, when using data from 13C-labeling experiments  to constrain fluxes, in particular the CO2 uptake rate (Fig. 3c), model predictions improved, especially for reactions associated with FA accumulation (PFL, FDHmq and EX_FA).
Overall, the iBP722 model allow us to evaluate the production of SA, showing that maximum theoretical yield for SA (1.1 g per g of glucose, assuming a symport system and no ATP requirements for maintenance) is comparable to those predicted using E. coli or S. cerevisiae models (1.1 and 0.8 g per g of glucose for anaerobic conditions using iJO1366  and iMM904 , respectively, under the same previous assumptions). Thus, metabolic capabilities seem equivalent to other organisms being exploited for SA production. It also allows describing the impact of growth conditions on the production of C3 and C4 fermentative by-products. The carbon split between C4 and C3 pathways has been investigated and showed to be influenced by various factors like the available reducing power (i.e. NADH/NAD+ ratio) or CO2 availability, therefore affecting the SA production in A. succinogenes growing cultures [64, 69]. In silico predicted SA yields are higher with more reduced carbon sources (e.g. minimum and maximum FVA yields of 0.44 and 0.46 g.g− 1 on sorbitol compared to 0.41 and 0.44 g.g− 1 on glucose, correspondingly, assuming the same sugar transport mechanism) and higher CO2 availability, favouring carbon flux through the C4 branch, as a consequence of higher reducing power and higher PPCK carboxylation activity. Moreover, flexibility in sugar transport mechanisms or enzymes reversibility (e.g. ME2) may lead to increased levels of SA, as carbon flux partitioning between C3 and C4 pathways was shown to be largely affected (Fig. 4). The reverse activity of ME2 has shown to increase SA yields, redirecting part of the carbon flux from the C3 toward the C4 branch through the carboxylation of pyruvate to L-malate, with the simultaneous production of reducing power (NADPH). On the other hand, transport activities limited to PEP:PTS-based systems decrease SA yields, since PEP is used as the energy source for sugar uptake generating pyruvate, which is necessarily consumed via the C3 branch.
In this work, the GSMM of A. succinogenes 130Z (iBP722) was reconstructed and validated using different sets of experimental data from literature. The reconstruction of the model included the compilation of functionally annotated metabolic genes and the corresponding coding proteins, as well as associated biochemical reactions. The model was complemented with a biomass equation and spontaneous and transport reactions. It was further amended after a gap filling process, including the identification of genetic evidences based on homology searches. Model accuracy to predict growth phenotypes and the production of fermentative by-products was evaluated using FVA and pFBA simulation methods. The ability to predict changes in carbon flux distributions due to environmental perturbations like CO2 limitations or alterations in the redox state was also tested. Predicted SA yields were in good agreement with experimental data, suggesting that the model is able to characterize the fermentative metabolism under various conditions. The increase in CO2 availability showed to have a positive impact in SA yields, which is consistent with reported data [48, 49]. As such, optimal conditions for increased SA yields may include increased CO2 availability, the use of more reduced carbon sources like sorbitol or the use of external energy source like H2 . Besides improving process conditions for the production of SA, the design of microbial strains by metabolic engineering to increase the flux through the C4 branch has been attempted, albeit with limited success . Model predictions show that changes in sugar transport systems, CO2 fixation activities or the reversibility of the malic enzyme may have an impact in the SA yield. Therefore, it is expected that modifications in the metabolism that would include these activities would further improve the SA yield.
Overall, the iBP722 model enables a better understanding of the metabolic behaviour and capabilities of this organism, which can be explored to further improve SA productivity. The capacity to consume a wide range of C5 and C6 sugars, as well as other low-cost carbon sources (e.g. glycerol or lactose) and its metabolic flexibility may provide some advantages over other SA-producing strains, like recombinant Escherichia coli or Mannheimia succiniciproducens .
- AAex :
Extracellular acetic acid
- EtOHex :
- FAex :
Extracellular formic acid
- FUMex :
- Glcex :
- SAex :
Extracellular succinic acid
- EX_CO2 :
Carbon dioxide exchange reaction
Acetic acid exchange reaction
Ethanol exchange reaction
Formic acid exchange reaction
Fumarate exchange reaction
Glucose exchange reaction
Succinic acid exchange reaction
Glucose 6-phosphate dehydrogenase of the oxidative pentose phosphate pathway
Van der Werf MJ, Guettler MV, Jain MK, Zeikus JG. Environmental and physiological factors affecting the succinate product ratio during carbohydrate fermentation by Actinobacillus sp. 130Z. Arch Microbiol. 1997;167:332–42.
Liu YP, Zheng P, Sun ZH, Ni Y, Dong JJ, Zhu LL. Economical succinic acid production from cane molasses by Actinobacillus succinogenes. Bioresour Technol. 2008;99:1736–42.
Urbance SE, Pometto AL, DiSpirito A a., Denli Y: Evaluation of succinic acid continuous and repeat-batch biofilm fermentation by Actinobacillus succinogenes using plastic composite support bioreactors. Appl Microbiol Biotechnol 2004, 65:664–670.
Guettler MV, Rumler D, Jainf MK. Actinobacillus succinogenes sp. nov., a novel succinic-acid-producing strain from the bovine rumen. Int J Syst Bacteriol. 1999;49(1 999):207–16.
Gunnarsson IB, Alvarado-morales M, Angelidaki I. Utilization of CO2 fixating bacterium Actinobacillus succinogenes 130Z for simultaneous biogas upgrading and biosuccinic acid production. Environ Sci Technol. 2014;48:12464–8.
Shen N, Qin Y, Wang Q, Liao S, Zhu J, Zhu Q, Mi H, Adhikari B, Wei Y, Huang R. Production of succinic acid from sugarcane molasses supplemented with a mixture of corn steep liquor powder and peanut meal as nitrogen sources by Actinobacillus succinogenes. Lett Appl Microbiol. 2015;60:544–51.
Zheng P, Dong JJ, Sun ZH, Ni Y, Fang L. Fermentative production of succinic acid from straw hydrolysate by Actinobacillus succinogenes. Bioresour Technol. 2009;100:2425–9.
Wan C, Li Y, Shahbazi A, Xiu S. Succinic acid production from cheese whey using Actinobacillus succinogenes 130 Z. Appl Biochem Biotechnol. 2008;145:111–9.
Chen K, Jiang M, Wei P, Yao J, Wu H. Succinic acid production from acid hydrolysate of corn fiber by Actinobacillus succinogenes. Appl Biochem Biotechnol. 2010;160:477–85.
Jiang M, Chen K, Liu Z, Wei P, Ying H, Chang H. Succinic acid production by Actinobacillus succinogenes using spent brewer’s yeast hydrolysate as a nitrogen source. Appl Biochem Biotechnol. 2010;160:244–54.
Xi YL, Chen KQ, Dai WY, Ma JF, Zhang M, Jiang M, Wei P, Ouyang PK. Succinic acid production by Actinobacillus succinogenes NJ113 using corn steep liquor powder as nitrogen source. Bioresour Technol. 2013;136:775–9.
Li Q, Yang M, Wang D, Li W, Wu Y, Zhang Y, Xing J, Su Z. Efficient conversion of crop stalk wastes into succinic acid production by Actinobacillus succinogenes. Bioresour Technol. 2010;101:3292–4.
Chen K, Zhang H, Miao Y, Jiang M, Chen J. Succinic acid production from enzymatic hydrolysate of sake lees using Actinobacillus succinogenes 130Z. Enzym Microb Technol. 2010;47:236–40.
Chen KQ, Li J, Ma JF, Jiang M, Wei P, Liu ZM, Ying HJ. Succinic acid production by Actinobacillus succinogenes using hydrolysates of spent yeast cells and corn fiber. Bioresour Technol. 2011;102:1704–8.
Yu J, Li Z, Ye Q, Yang Y, Chen S. Development of succinic acid production from corncob hydrolysate by Actinobacillus succinogenes. J Ind Microbiol Biotechnol. 2010;37:1033–40.
Zeikus JG, Jain MK, Elankovan P. Biotechnology of succinic acid production and markets for derived industrial products. Appl Microbiol Biotechnol. 1999;51:545–52.
Kamzolova SV, Vinokurova NG, Shemshura ON, Bekmakhanova NE, Lunina JN, Samoilenko VA, Morgunov IG. The production of succinic acid by yeast Yarrowia lipolytica through a two-step process. Appl Microbiol Biotechnol. 2014;98:7959–69.
Raab AM, Gebhardt G, Bolotina N, Weuster-Botz D, Lang C. Metabolic engineering of Saccharomyces cerevisiae for the biotechnological production of succinic acid. Metab Eng. 2010;12:518–25.
Yan D, Wang C, Zhou J, Liu Y, Yang M, Xing J. Construction of reductive pathway in Saccharomyces cerevisiae for effective succinic acid fermentation at low pH value. Bioresour Technol. 2014;156:232–9.
Liu Y, Wu H, Li Q, Tang X, Li Z, Ye Q. Process development of succinic acid production by Escherichia coli NZN111 using acetate as an aerobic carbon source. Enzym Microb Technol. 2011;49:459–64.
Scholten E, Renz T, Thomas J. Continuous cultivation approach for fermentative succinic acid production from crude glycerol by Basfia succiniciproducens DD1. Biotechnol Lett. 2009;31:1947–51.
Kim DY, Yim SC, Lee PC, Lee WG, Lee SY, Chang HN. Batch and continuous fermentation of succinic acid from wood hydrolysate by Mannheimia succiniciproducens MBEL55E. In Enzyme Microb Technol Volume. 2004;35:648–53.
Lee SY, Kim JM, Song H, Lee JW, Kim TY, Jang Y-S. From genome sequence to integrated bioprocess for succinic acid production by Mannheimia succiniciproducens. Appl Microbiol Biotechnol. 2008;79:11–22.
Cok B, Tsiropoulos I, Roes AL, Patel MK. Succinic acid production derived from carbohydrates: an energy and greenhouse gas assessment of a platform chemical toward a bio-based economy. Biofuels Bioprod Biorefin. 2014;8:16–29.
Choi S, Song CW, Shin JH, Lee SY. Biorefineries for the production of top building block chemicals and their derivatives. Metab Eng. 2015:223–39.
Jansen MLA, van Gulik WM. Towards large scale fermentative production of succinic acid. Curr Opin Biotechnol. 2014;30:190–7.
Lee SJ, Song H, Lee SY. Genome-based metabolic engineering of Mannheimia succiniciproducens for succinic acid production. Appl Environ Microbiol. 2006;72:1939–48.
Yan Q, Zheng P, Dong JJ, Sun ZH. A fibrous bed bioreactor to improve the productivity of succinic acid by Actinobacillus succinogenes. J Chem Technol Biotechnol. 2014;89:1760–6.
McKinlay JB, Zeikus JG, Vieille C. Insights into Actinobacillus succinogenes fermentative metabolism in a chemically defined growth medium. Appl Environ Microbiol. 2005;71:6651–6.
McKinlay JB, Vieille C. 13C-metabolic flux analysis of Actinobacillus succinogenes fermentative metabolism at different NaHCO3 and H2 concentrations. Metab Eng. 2008;10:55–68.
Li J, Jiang M, Chen K-QQ, Ye Q, Shang L-AA, Wei P, Ying H-JJ, Chang H-NN. Effect of redox potential regulation on succinic acid production by Actinobacillus succinogenes. Bioprocess Biosyst Eng. 2010;33:911–20.
Carvalho M, Matos M, Roca C, Reis M a M. Succinic acid production from glycerol by Actinobacillus succinogenes using dimethylsulfoxide as electron acceptor. New Biotechnol. 2014;31:133–9.
van Heerden CD, Nicol W. Continuous succinic acid fermentation by Actinobacillus succinogenes. Biochem Eng J. 2013;73:5–11.
Otero JM, Cimini D, Patil KR, Poulsen SG, Olsson L, Nielsen J. Industrial systems biology of Saccharomyces cerevisiae enables novel succinic acid cell factory. PLoS One. 2013;8
Sajo Mienda B, Shahir Shamsir M. Model-driven in silico glpC gene knockout predicts increased succinate production from glycerol in Escherichia coli. AIMS Bioeng. 2015;2:40–8.
Xu G, Zou W, Chen X, Xu N, Liu L, Chen J. Fumaric acid production in Saccharomyces cerevisiae by in silico aided metabolic engineering. PLoS One. 2012;7
Raab AM, Lang C. Oxidative versus reductive succinic acid production in the yeast Saccharomyces cerevisiae. Bioeng Bugs. 2011;2
Cheng KK, Wang GY, Zeng J, Zhang JA. Improved succinate production by metabolic engineering. Biomed Res Int. 2013;2013
Lin H, Bennett GN, San KY. Fed-batch culture of a metabolically engineered Escherichia coli strain designed for high-level succinate production and yield under aerobic conditions. Biotechnol Bioeng. 2005;90:775–9.
Liu R, Liang L, Chen K, Ma J, Jiang M, Wei P, Ouyang P. Fermentation of xylose to succinate by enhancement of ATP supply in metabolically engineered Escherichia coli. Appl Microbiol Biotechnol. 2012;94:959–68.
Zhang X, Jantama K, Moore JC, Jarboe LR, Shanmugam KT, Ingram LO. Metabolic evolution of energy-conserving pathways for succinate production in Escherichia coli. Proc Natl Acad Sci U S A. 2009;106:20180–5.
ya LL, ming LR, feng MJ, quan CK, Jiang M, Wei P. Increased production of succinic acid in Escherichia coli by overexpression of malate dehydrogenase. Biotechnol Lett. 2011;33:2439–44.
Yang J, Wang Z, Zhu N, Wang B, Chen T, Zhao X. Metabolic engineering of Escherichia coli and in silico comparing of carboxylation pathways for high succinate productivity under aerobic conditions. Microbiol Res. 2014;169:432–40.
Jiang M, Ma J, Wu M, Liu R, Liang L, Xin F, Zhang W, Jia H, Dong W. Progress of succinic acid production from renewable resources: metabolic and fermentative strategies. Bioresour Technol. 2017;
Zou W, Zhu L-W, Li H-M, Tang Y-J. Significance of CO2 donor on the production of succinic acid by Actinobacillus succinogenes ATCC 55618. Microb Cell Factories. 2011;10:87.
Xi Y, Chen K, Li J, Fang X, Zheng X, Sui S, Jiang M, Wei P. Optimization of culture conditions in CO2 fixation for succinic acid production using Actinobacillus succinogenes. J Ind Microbiol Biotechnol. 2011;38:1605–12.
McKinlay JB, Laivenieks M, Schindler BD, McKinlay A a, Siddaramappa S, Challacombe JF, Lowry SR, Clum A, Lapidus AL, Burkhart KB, Harkins V, Vieille C. A genomic perspective on the potential of Actinobacillus succinogenes for industrial succinate production. BMC Genomics 2010, 11:680.
McKinlay JB, Shachar-Hill Y, Zeikus JG, Vieille C. Determining Actinobacillus succinogenes metabolic pathways and fluxes by NMR and GC-MS analyses of 13C-labeled metabolic product isotopomers. Metab Eng. 2007;9:177–92.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.
Eddy SR. Profile hidden Markov models. Bioinformatics. 1998;14(9):755–63.
The UniProt Consortium. UniProt: the universal protein knowledgebase. Nucleic Acids Res. 2018;46(5):2699.
Dias O, Rocha M, Ferreira EC, Rocha I. Reconstructing genome-scale metabolic models with merlin. Nucleic Acids Research. 2015;43(8):3899–910.
Kanehisa M, Goto S. KEGG: Kyoto Encyclopaedia of genes and genomes. Nucl Acids Res. 2000;28:27–30.
Saier MH, Reddy VS, Tamang DG, Västermark Å. The transporter classification database. Nucleic Acids Res. 2014;42
Caspi R, Billington R, Ferrer L, Foerster H, Fulcher CA, Keseler IM, Kothari A, Krummenacker M, Latendresse M, Mueller LA, Ong Q, Paley S, Subhraveti P, Weaver DS, Karp PD. The MetaCyc database of metabolic pathways and enzymes and the BioCyc collection of pathway/genome databases. Nucleic Acids Res. 2016;44:D471–80.
Flamholz A, Noor E, Bar-even A, Milo R. EQuilibrator - the biochemical thermodynamics calculator. Nucleic Acids Res. 2012;40
Santos S, Rocha I. A computation tool for the estimation of biomass composition from genomic and transcriptomic information. Adv. Intell. Syst. Comput. 2016:161–9.
Juty N, Ali R, Glont M, Keating S, Rodriguez N, Swat MJ, Wimalaratne SM, Hermjakob H, Le Novère N, Laibe C, Chelliah V. BioModels: content, features, functionality, and use. CPT Pharmacometrics Syst Pharmacol. 2015:55–68.
Rocha I, Maia P, Evangelista P, Vilaça P, Soares S, Pinto JP, Nielsen J, Patil KR, Ferreira EC, Rocha M. OptFlux: an open-source software platform for in silico metabolic engineering. BMC Syst Biol. 2010;4:45.
Bradfield MFA, Nicol W. Continuous succinic acid production from xylose by Actinobacillus succinogenes. Bioprocess Biosyst Eng. 2016;39:233–44.
Li J, Jiang M, Chen K, Shang L, Wei P, Ying H, Ye Q, Ouyang P, Chang H. Enhanced production of succinic acid by Actinobacillus succinogenes with reductive carbon source. Process Biochem. 2010;45:980–5.
Schindler BD, Joshi RV, Vieille C. Respiratory glycerol metabolism of Actinobacillus succinogenes 130Z for succinate production. J Ind Microbiol Biotechnol. 2014;41:1339–52.
Lin SKC, Du C, Koutinas A, Wang R, Webb C. Substrate and product inhibition kinetics in succinic acid production by Actinobacillus succinogenes. Biochem Eng J. 2008;41:128–35.
Orth JD, Conrad TM, Na J, Lerman JA, Nam H, Feist AM, Palsson BØ. A comprehensive genome-scale reconstruction of Escherichia coli metabolism—2011. Mol Syst Biol. 2011;
Mo ML, Palsson BO, Herrgård MJ. Connecting extracellular metabolomic measurements to intracellular flux states in yeast. BMC Syst Biol. 2009;3:37.
Pateraki C, Patsalou M, Vlysidis A, Kopsahelis N, Webb C, Koutinas AA, Koutinas M. Actinobacillus succinogenes: advances on succinic acid production and prospects for development of integrated biorefineries. Biochem Eng J. 2016;112:285–303.
Guarnieri MT, Chou YC, Salvachúa D, Mohagheghi A, St. John PC, Peterson DJ, Bomble YJ, Beckham GT. Metabolic engineering of Actinobacillus succinogenes provides insights into succinic acid biosynthesis. Appl Environ Microbiol. 2017;83
Song H, Lee JW, Choi S, You JK, Hong WH, Lee SY. Effects of dissolved CO2 levels on the growth of Mannheimia succiniciproducens and succinic acid production. Biotechnol Bioeng. 2007;98:1296–304.
Brink HG, Nicol W. Succinic acid production with Actinobacillus succinogenes: rate and yield analysis of chemostat and biofilm cultures. Microb Cell Factories. 2014;13:111.
Lin H, Bennett GN, San KY. Metabolic engineering of aerobic succinate production systems in Escherichia coli to improve process productivity and achieve the maximum theoretical succinate yield. Metab Eng. 2005;7:116–27.
The authors want to thank Apostolis Koutinas and co-workers from Agricultural University of Athens for the helpful discussions.
Financially supported by BRIGIT (KBBE-2012-6-311935, FP7 project Contract nr 311935) and by the Portuguese Foundation for Science and Technology (FCT) under the scope of the strategic funding of UID/BIO/04469/2013 unit and COMPETE 2020 (POCI-01-0145-FEDER-006684), in addition to the BioTecNorte operation (NORTE-01-0145-FEDER-000004) funded by European Regional Development Fund under the scope of Norte2020 - Programa Operacional Regional do Norte.
Availability of data and materials
The GSMM generated during the current study is available at http://darwin.di.uminho.pt/models and at BioModels with the identifier MODEL1804130001. A. succinogenes 130Z complete genome - GenBank accession number NC_009655.
Ethics approval and consent to participate
The authors declare that they have no competing interests
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Biomass composition of A. succinogenes 130Z. (DOCX 43 kb)
Details on the iBP722 model. (DOCX 24 kb)
Testing A. succinogenes growth on different carbon sources. (DOCX 22 kb)
Details on transport mechanisms included in the iBP722 model. (DOCX 37 kb)
Growth predictions on glycerol with and without DMSO. (DOCX 21 kb)
About this article
Cite this article
Pereira, B., Miguel, J., Vilaça, P. et al. Reconstruction of a genome-scale metabolic model for Actinobacillus succinogenes 130Z. BMC Syst Biol 12, 61 (2018). https://0-doi-org.brum.beds.ac.uk/10.1186/s12918-018-0585-7
- Genome-scale metabolic reconstruction
- Constraints-based flux analysis
- Succinic acid fermentation