Skip to main content

Metabolic modelling of polyhydroxyalkanoate copolymers production by mixed microbial cultures



This paper presents a metabolic model describing the production of polyhydroxyalkanoate (PHA) copolymers in mixed microbial cultures, using mixtures of acetic and propionic acid as carbon source material. Material and energetic balances were established on the basis of previously elucidated metabolic pathways. Equations were derived for the theoretical yields for cell growth and PHA production on mixtures of acetic and propionic acid as functions of the oxidative phosphorylation efficiency, P/O ratio. The oxidative phosphorylation efficiency was estimated from rate measurements, which in turn allowed the estimation of the theoretical yield coefficients.


The model was validated with experimental data collected in a sequencing batch reactor (SBR) operated under varying feeding conditions: feeding of acetic and propionic acid separately (control experiments), and the feeding of acetic and propionic acid simultaneously. Two different feast and famine culture enrichment strategies were studied: (i) either with acetate or (ii) with propionate as carbon source material. Metabolic flux analysis (MFA) was performed for the different feeding conditions and culture enrichment strategies. Flux balance analysis (FBA) was used to calculate optimal feeding scenarios for high quality PHA polymers production, where it was found that a suitable polymer would be obtained when acetate is fed in excess and the feeding rate of propionate is limited to ~0.17 C-mol/(C-mol.h). The results were compared with published pure culture metabolic studies.


Acetate was more conducive toward the enrichment of a microbial culture with higher PHA storage fluxes and yields as compared to propionate. The P/O ratio was not only influenced by the selected microbial culture, but also by the carbon substrate fed to each culture, where higher P/O ratio values were consistently observed for acetate than propionate. MFA studies suggest that when mixtures of acetate and propionate are fed to the cultures, the catabolic activity is primarily guaranteed through acetate uptake, and the characteristic P/O ratio of acetate prevails over that of propionate. This study suggests that the PHA production process by mixed microbial cultures has the potential to be comparable or even more favourable than pure cultures.


Polyhydroxyalkanoates are biopolymers with physico-chemical properties similar to polypropylene but with the advantage of being biodegradable and biocompatible. Industrial PHA production technology is currently based on bacteria cultivation using pure cultures grown in well-defined nutrient deficient synthetic media with single substrates [1]. The relatively high production costs constitute presently the most important barrier to PHA becoming a commodity polymer in direct competition with oil-based polymers. The cost of the substrate and of the equipment required for aseptic operation is responsible for about 40% of the total PHA production cost [1].

Potential strategies for reducing both the operational and capital expenses of PHA production are the use of open mixed cultures of microorganisms instead of pure cultures, and the use of waste materials as the carbon substrate. Activated sludge is known to accumulate up to 65% PHA per cell dry weight [24]. There are currently many possibilities of industrial and agricultural wastes that could be directed to the production of PHA on the basis of open mixed microbial systems [5, 6].

Mixed microbial cultures cannot directly metabolize sugars into PHA. A pre-fermentation step is required in order to convert sugars into volatile fatty acids (VFA) such as acetic, propionic, valeric and butyric acid. This pre-fermentation step can, however, be accomplished with high yield [5, 7]. The final PHA molecular structure is highly dependent on the VFA composition used as carbon source. The most well-studied case is the production of polyhydroxybutyrate (PHB), a homopolymer, obtained when acetic acid is used as the sole carbon source. The metabolism for this transformation (included in Figure 1) can be defined by 6 basic reactions [8, 9]: acetate uptake (R1), growth on acetyl-CoA (R4), catabolism (R6), oxidative phosphorylation (R8), maintenance (R7) and PHB production (RPHB). A metabolic model based on these reactions was comprehensively validated with experimental data [10]. This study demonstrated an excellent agreement between the model predictions and experimental data.

Figure 1
figure 1

PHA metabolic network. Figure shows schematic representation of the metabolic reactions, substrates, products (both marked by full lines) and metabolites (marked by dashed lines) involved in PHA production by mixed microbial cultures.

The properties and quality of PHA are highly dependent on the monomeric composition, which, in turn, is dependent on the VFA mixture adopted as carbon source. PHB is a relatively low quality polymer since it is too brittle for many applications [1, 11]. In terms of PHA processing, the production of copolymers appears to be of higher commercial interest [1214]. The production of PHB homopolymer and copolymers of 3-hydroxybutyrate (3HB), 3-hydroxyvalerate (3HV) and 3-hydroxy-2-methylvalerate (3H2MV) from mixtures of acetic and propionic acid has been studied by several researchers [11, 1517]. There are, however, no metabolic models currently available describing the production of such copolymers. Metabolic models can be used to facilitate process optimization, and also serve as a reference basis in the interpretation of data arising from the study of biological processes. The comparison of experimentally-determined stoichiometry with the theoretical model predictions allows a better understanding of the processes under study. The main objective of this work was to extend existing PHB models [8, 10] to the production of PHA copolymers from mixtures of VFA. This study will be focused on the PHA production phase, and is confined to mixtures of acetic and propionic acid, which are seen as the most promising for the synthesis of high quality PHAs [11, 18].

Metabolic model

A schematic of the metabolic model developed in this study is shown in Figure 1. Acetate and propionate are taken up inside the cell, and used for PHA production and biomass growth. Energy is generated within the cells through oxidative phosphorylation, while a portion of this energy is also necessary for cell maintenance. A detailed description of the biochemical reactions, material and energy balances, kinetic reactions and metabolic fluxes is provided below.

Basic reactions

The detailed biochemical reactions involved in the PHA production process using acetate and propionate carbon sources under aerobic conditions are shown in Table 1. All reactions are expressed on a carbon-mole basis. In the feast phase, acetate and propionate are taken up inside the cells by means of active transport, requiring one mole of ATP per mole of carbon source. They are then converted into acetyl-CoA and propionyl-CoA, respectively, at the cost of another mole of ATP per mole of acetate or propionate. The net reactions are shown in R1 and R2 [19].

Table 1 Metabolic model of PHA production from mixtures of acetate and propionate by mixed microbial cultures

While acetate is converted to acetyl-CoA and propionyl-CoA is produced from propionate, it has also been observed that a portion of the propionyl-CoA produced from propionate uptake is converted into acetyl-CoA [11]. This transformation may potentially proceed through 5 different biochemical pathways [11], however, the net material and reducing power transformations are identical through each pathway. The difference between the pathways lies only in the amount of energy generated. For the purposes of this model, it was assumed that propionyl-CoA was first converted to succinyl-CoA via the methylmalonyl-CoA pathway, then converted to oxaloacetate, pyruvate and acetyl-CoA. The net reaction is described in R3.

In PHA production, acetyl-CoA and propionyl-CoA are reduced and condensed to produce a polymer consisting of various 3-hydroxyalkanoates (3HAs) monomers. 3HB is formed from 2 acetyl-CoA monomers, one acetyl-CoA and one propionyl-CoA combine to form either 3HV or 3-hydroxy-2-methylbutyrate (3H2MB) (which are isomers of each other, and henceforth expressed only as 3HV, indicating the sum of these 2 compounds), while 2 propionyl-CoA monomers form 3H2MV. The fraction of each PHA monomer produced by the cells depends on the fluxes of acetyl-CoA and propionyl-CoA available for polymer production. It also depends on whether the cells tend to reduce and condense acetyl-CoA and propionyl-CoA in a random fashion, or whether preferential binding occurs between e.g. acetyl-CoA and propionyl-CoA. For the purposes of metabolic model development, however, it is necessary to express each reaction in relation to the overall mass, energy and redox transformations. Thus, it is more convenient to express the various potential PHA fractions produced by the cells in terms of the reduced and condensed monomers produced from acetyl-CoA or propionyl-CoA. As shown in R9 and R10, these monomers are represented as acetyl-CoA* and propionyl-CoA*, respectively, in a similar fashion as expressed by a previous study [20]. The PHA formation process from acetyl-CoA* and propionyl-CoA* has no further demand for energy or reducing power.

Biomass growth occurs from acetyl-CoA and propionyl-CoA, respectively. The energy required for the production of 1 C-mol of biomass from acetyl-CoA is represented as K1, which has been estimated as 1.7 mol ATP [21]. The reaction of biomass growth from acetyl-CoA was determined through stoichiometric and redox balancing, and is shown in R4. The amount of energy required for biomass synthesis from propionyl-CoA is represented as K2, which has been estimated to be 1.38 mol ATP per C-mol of biomass formed [22]. Similarly, it was assumed that propionyl-CoA is converted to succinyl-CoA for biomass synthesis [22]. The net reaction is shown in R5. The biomass formula was assumed to be CH1.4O0.4N0.2 [23].

Catabolism of acetyl-CoA and propionyl-CoA also occurs in the microbial cells. Acetyl-CoA is converted to CO2 through the tricarboxylic acid cycle (TCA) as described in R6 [19]. The catabolism of propionyl-CoA is assumed to occur via acetyl-CoA, where R3 is followed by R6.

Energy in the form of ATP is produced from NADH2 through oxidative phosphorylation. The amount of ATP generated per mole of NADH2 oxidized is expressed by the P/O ratio, δ, which represents the efficiency of oxidative phosphorylation [24]. This reaction is expressed in R8.

In addition to the energetic transformations described above, energy is also required for cell maintenance. The rate of ATP consumption for maintenance purposes is described in R7 as mATP.

Pseudo steady-state material balancing of intracellular intermediates

The metabolic network of Figure 1 has q = 13 metabolic reactions, m = 6 intracellular metabolites (marked by dashed lines), 4 input substrates (acetic and propionic acid, ammonia and oxygen) and 5 end-products [biomass, PHB, poly(3-hydroxyvalerate) (PHV), poly(2-methyl-3-hydroxyvalerate) (PH2MV) and CO2]. The steady state material balances to the m intracellular intermediates (acetyl-CoA, acetyl-CoA*, propionyl-CoA, propionyl-CoA*, NADH2 and ATP) are:

Ac-CoA: R1 + R3 - 1.27·R4 - R6 - R9= 0

Ac-CoA*: R9 - R PHB - 0.4·R PHV = 0

Prop-CoA: R2 - 1.5·R3 - 1.06·R5 - R10 = 0

Prop-CoA*: R10 - 0.6·R PHV - RPH 2MV= 0

NADH 2 : 1.5·R3 + 0.53·R4 + 0.47·R5 + 2·R6 - R8 - 0.25·R9 - 0.17·R10 = 0

ATP: -R1 - 0.67·R2 - K1·R4 - K2·R5 + 0.5·R6 - R7 + δ·R8 = 0

Constraints to the metabolic network

From acetyl-CoA* and propionyl-CoA*, the resulting PHA polymer can be comprised of either PHB (2 acetyl-CoA* molecules), PHV (one acetyl-CoA* and one propionyl-CoA*) and PH2MV (2 propionyl-CoA*). The PHA composition depends on whether selective or random condensation of acetyl-CoA* and propionyl-CoA* takes place. In the propionate enriched culture selective condensation was observed, where acetyl-CoA* preferentially bounds with propionyl-CoA*, forming PHV [11]. Since acetyl-CoA* was generated in higher abundance than propionyl-CoA*, the remaining acetyl-CoA* condensed to form PHB, while PH2MV was not produced by the sludge. This selective condensation is expressed by the following three additional equations:

R P H B = R 9 2 3 R 10 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOuai1aaSbaaSqaaiabdcfaqjabdIeaijabdkeacbqabaGccqGH9aqpcqWGsbGudaWgaaWcbaGaeGyoaKdabeaakiabgkHiTKqbaoaalaaabaGaeGOmaidabaGaeG4mamdaaOGaeyyXICTaemOuai1aaSbaaSqaaiabigdaXiabicdaWaqabaaaaa@3D3A@
R P H V = 5 3 R 10 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOuai1aaSbaaSqaaiabdcfaqjabdIeaijabdAfawbqabaGccqGH9aqpjuaGdaWcaaqaaiabiwda1aqaaiabiodaZaaakiabgwSixlabdkfasnaaBaaaleaacqaIXaqmcqaIWaamaeqaaaaa@3A18@

RPH 2MV= 0 (2c)

since RPH 2MV= 0 (Eq. 2c), Eqs. (2a–b) are linearly dependent to Eqs. (1a–f), only a single additional constraint is introduced with RPH2MV = 0.

In the acetate enriched culture, the PHA composition consisted of a lower PHV fraction and a higher PH2MV and PHB fraction, as is expected by microorganisms performing random condensation of acetyl-CoA* and propionyl-CoA*. In this latter case the following three equations are applied [20]:

R P H B = 2 R 9 2 2 R 9 + 3 R 10 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOuai1aaSbaaSqaaiabdcfaqjabdIeaijabdkeacbqabaGccqGH9aqpjuaGdaWcaaqaaiabikdaYiabgwSixlabdkfasnaaDaaabaGaeGyoaKdabaGaeGOmaidaaaqaaiabikdaYiabgwSixlabdkfasnaaBaaabaGaeGyoaKdabeaacqGHRaWkcqaIZaWmcqGHflY1cqWGsbGudaWgaaqaaiabigdaXiabicdaWaqabaaaaaaa@45CC@
R P H V = 5 R 9 R 10 2 R 9 + 3 R 10 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOuai1aaSbaaSqaaiabdcfaqjabdIeaijabdAfawbqabaGccqGH9aqpjuaGdaWcaaqaaiabiwda1iabgwSixlabdkfasnaaBaaabaGaeGyoaKdabeaacqGHflY1cqWGsbGudaWgaaqaaiabigdaXiabicdaWaqabaaabaGaeGOmaiJaeyyXICTaemOuai1aaSbaaeaacqaI5aqoaeqaaiabgUcaRiabiodaZiabgwSixlabdkfasnaaBaaabaGaeGymaeJaeGimaadabeaaaaaaaa@4A7D@
R P H 2 M V = 3 R 10 2 2 R 9 + 3 R 10 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOuai1aaSbaaSqaaiabdcfaqjabdIeaijabikdaYiabd2eanjabdAfawbqabaGccqGH9aqpjuaGdaWcaaqaaiabiodaZiabgwSixlabdkfasnaaDaaabaGaeGymaeJaeGimaadabaGaeGOmaidaaaqaaiabikdaYiabgwSixlabdkfasnaaBaaabaGaeGyoaKdabeaacqGHRaWkcqaIZaWmcqGHflY1cqWGsbGudaWgaaqaaiabigdaXiabicdaWaqabaaaaaaa@48E9@

The R9 and R10 fluxes are sufficient to calculate the fluxes of R PHB , R PHV and RPH 2MV.

On the other hand, constraints (3a–c) automatically obey the material balances (1a–f) (see [20] for details), thus also a single additional constraint can also be added to the material balances here (1a–f).

Dynamic material balancing of substrates and end-products

The transient material balances in a batch system of input substrates (acetate, propionate, active biomass, ammonia) and intracellular contents of PHB, PHV and PH2MV are the following:

d A c d t = R 1 X MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqWGKbazcqWGbbqqcqWGJbWyaeaacqWGKbazcqWG0baDaaGccqGH9aqpcqGHsislcqWGsbGudaWgaaWcbaGaeGymaedabeaakiabgwSixlabdIfaybaa@3B01@
d Pr o p d t = R 2 X MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqWGKbazcyGGqbaucqGGYbGCcqWGVbWBcqWGWbaCaeaacqWGKbazcqWG0baDaaGccqGH9aqpcqGHsislcqWGsbGudaWgaaWcbaGaeGOmaidabeaakiabgwSixlabdIfaybaa@3E0F@
d N d t = 0.2 ( R 4 + R 5 ) X MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqWGKbazcqWGobGtaeaacqWGKbazcqWG0baDaaGccqGH9aqpcqGHsislcqaIWaamcqGGUaGlcqaIYaGmcqGHflY1cqGGOaakcqWGsbGudaWgaaWcbaGaeGinaqdabeaakiabgUcaRiabdkfasnaaBaaaleaacqaI1aqnaeqaaOGaeiykaKIaeyyXICTaemiwaGfaaa@43CF@
d X d t = ( R 4 + R 5 ) X MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqWGKbazcqWGybawaeaacqWGKbazcqWG0baDaaGccqGH9aqpcqGGOaakcqWGsbGudaWgaaWcbaGaeGinaqdabeaakiabgUcaRiabdkfasnaaBaaaleaacqaI1aqnaeqaaOGaeiykaKIaeyyXICTaemiwaGfaaa@3DE8@
d f P H B d t = R P H B ( R 4 + R 5 ) f P H B MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqWGKbazcqWGMbGzdaWgaaqaaiabdcfaqjabdIeaijabdkeacbqabaaabaGaemizaqMaemiDaqhaaOGaeyypa0JaemOuai1aaSbaaSqaaiabdcfaqjabdIeaijabdkeacbqabaGccqGHsislcqGGOaakcqWGsbGudaWgaaWcbaGaeGinaqdabeaakiabgUcaRiabdkfasnaaBaaaleaacqaI1aqnaeqaaOGaeiykaKIaeyyXICTaemOzay2aaSbaaSqaaiabdcfaqjabdIeaijabdkeacbqabaaaaa@4AAA@
d f P H V d t = R P H V ( R 4 + R 5 ) f P H V MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqWGKbazcqWGMbGzdaWgaaqaaiabdcfaqjabdIeaijabdAfawbqabaaabaGaemizaqMaemiDaqhaaOGaeyypa0JaemOuai1aaSbaaSqaaiabdcfaqjabdIeaijabdAfawbqabaGccqGHsislcqGGOaakcqWGsbGudaWgaaWcbaGaeGinaqdabeaakiabgUcaRiabdkfasnaaBaaaleaacqaI1aqnaeqaaOGaeiykaKIaeyyXICTaemOzay2aaSbaaSqaaiabdcfaqjabdIeaijabdAfawbqabaaaaa@4B22@
d f P H 2 M V d t = R P H 2 M V ( R 4 + R 5 ) f P H 2 M V MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqWGKbazcqWGMbGzdaWgaaqaaiabdcfaqjabdIeaijabikdaYiabd2eanjabdAfawbqabaaabaGaemizaqMaemiDaqhaaOGaeyypa0JaemOuai1aaSbaaSqaaiabdcfaqjabdIeaijabikdaYiabd2eanjabdAfawbqabaGccqGHsislcqGGOaakcqWGsbGudaWgaaWcbaGaeGinaqdabeaakiabgUcaRiabdkfasnaaBaaaleaacqaI1aqnaeqaaOGaeiykaKIaeyyXICTaemOzay2aaSbaaSqaaiabdcfaqjabdIeaijabikdaYiabd2eanjabdAfawbqabaaaaa@5161@

Theoretical yields

The theoretical yields and VFA and oxygen maintenance coefficients were derived analytically from the material balance equations (1a–f). The expressions are compiled in Table 2. The theoretical yield coefficients are functions of δ, K1, K2 and y, the ratio of propionate flux to total carbon flux (see Eq. 5):

Table 2 Theoretical yields and maintenance coefficients
y = R 2 R 1 + R 2 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyEaKNaeyypa0tcfa4aaSaaaeaacqWGsbGudaWgaaqaaiabikdaYaqabaaabaGaemOuai1aaSbaaeaacqaIXaqmaeqaaiabgUcaRiabdkfasnaaBaaabaGaeGOmaidabeaaaaaaaa@36E2@

Kinetic model

As stated above, the metabolic network has q = 13 fluxes and a total number of 7 constraints. Therefore, at least 6 fluxes must be defined a priori in order to calculate the network fluxes. Table 3 compiles the kinetic equations used in this work to model the missing fluxes. This model is an extension of a previous acetate to PHB metabolic model [10]. Briefly, the model defines the following kinetic rates:

Table 3 Kinetic model
  • Growth on acetyl-CoA (R4) is limited by the concentrations of acetate and ammonia.

  • Growth on propionyl-CoA (R5) is limited by the concentrations of propionate and ammonia.

  • Catabolism (R6) is limited by total VFA concentration. The kinetic equation for catabolism is only used for mixtures of acetate and propionate to fulfil the 6 fluxes required to perform MFA.

  • Maintenance (R7) is limited by total VFA concentration.

  • Acetyl-CoA* synthesis (R9) is limited by total VFA concentration and inhibited by the intracellular PHA content. Intracellular PHA inhibition has been validated experimentally by several authors [810, 2527].

  • Propionyl-CoA* synthesis (R10) is limited by propionate concentration and inhibited by the intracellular PHA content.


Reactors operation

In this study, two SBR for the production of PHA, were operated over four years. The sludge in each reactor was adapted to either acetate or propionate as the sole carbon source. The reactors working volume was 1 litre and the total SBR cycle duration was 12 h, consisting of 10.5 h of aerobiosis, 1 h of settling (agitation and air sparging switched off) and 0.5 h to withdraw half of the volume, which was replaced by the same volume of fresh medium during the first 15 min at the beginning of the next cycle. The hydraulic retention time (HRT) was therefore 1 day. At the end of each cycle, before settling, a defined volume of biomass was removed to keep the sludge retention time (SRT) at 10 days. Oxygen was supplied by an air compressor through a ceramic membrane disperser introduced inside the reactor at an airflow rate of 1.0 vvm (volume air/(volume reactor.min)), allowing for the dissolved oxygen (DO) concentration to be around 80% of the saturation value. The reactors were operated without pH control but its value was monitored on-line and ranged between 8.0 and 9.2; the temperature was controlled at 22°C and the stirring rate at 250 rpm.

Culture medium

The standard medium used in the SBRs was composed of (per litre of distilled water): 1.269 g CH3CH2COOH or 4.0796 g of CH3COONa.3H2O, 600 mg MgSO4.7H2O, 160 mg NH4Cl, 100 mg EDTA, 92 mg K2HPO4, 45 mg KH2PO4, 70 mg CaCl2.2H2O and 2 ml of trace elements solution. The trace solution consisted of (per litre of distilled water): 1500 mg FeCl3.6H2O, 150 mg H3BO3, 150 mg CoCl2.6H2O, 120 mg MnCl2.4H2O, 120 mg ZnSO4.7H2O, 60 mg Na2MoO4.2H2O, 30 mg CuSO4.5H2O and 30 mg of KI. Thiourea (10 mg/l) was added to inhibit nitrification. The pH of the salt solution was adjusted to 7.2 and then sterilized, where the phosphorus components of the solution was sterilized separately. After sterilization, the two solutions were allowed to cool, and were then mixed together.

In the batch experiments, different concentrations of acetate, propionate and ammonia were tested for the two systems and different feeding regimens were used. For the acetate reactor: the acetate concentrations tested with 1.4 N-mmol/l of ammonia were 15 C-mmol/l, 30 C-mmol/l and 60 C-mmol/l. The ammonia concentrations with 30 Cmmol/l of acetate: 0.7 N-mmol/l, 1.4 N-mmol/l and 2.8 N-mmol/l. The feeding regimen was tested by supplying multiple pulses of 60 C-mmol/l of acetate (3 tests with 3 pulses and 1 test with 4 pulses), where 0.7 N-mmol/l of ammonia was also added in the first pulse for these tests. Two more assays were performed by supplying a 30 C-mmol/l of propionate and a mixture of acetate and propionate with 15 C-mmol/l each. The batch tests performed in the propionate system analyzed the effect of different concentrations of propionate (30 C-mmol/l, 60 C-mmol/l, 90 C-mmol/l and 120 C-mmol/l) with 1.4 N-mmol/l of ammonia and the use of acetate (in one pulse of 30 C-mmol/l and mixed with propionate, 15 C-mmol/l each). In total, 18 batch tests were performed with the two systems.

Analytical techniques

Cell dry weight was determined as volatile suspended solids (VSS), according to Standard Methods [28]. Acetate and propionate were analyzed by HPLC using a BioRad Aminex HPX-87H column, with 0.01 N sulphuric acid as eluent, an elution rate of 0.6 ml/min and an operating temperature of 50°C. A UV detector (Merck) set at 210 nm was used. Prior to injection, samples were filtered using a 0.2 μm membrane. PHA were determined by GC after acidic estherification (see [4, 11] for details). Ammonia was analyzed using an ammonia gas sensing combination electrode (ThermoOrion 9512).

Parameter estimation

The kinetic parameters were estimated by non-linear weighted least-squares (MATLAB's lsqnonlin function) using the Levenberg-Marquardt algorithm. The program minimized the root mean squared error:

r m s e = e T e ( n p ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOCaiNaemyBa0Maem4CamNaemyzauMaeyypa0ZaaOaaaKqbagaadaWcaaqaaiabhwgaLnaaCaaabeqaaiabdsfaubaacqGHflY1cqWHLbqzaeaadaqadaqaaiabd6gaUjabgkHiTiabdchaWbGaayjkaiaawMcaaaaaaSqabaaaaa@3F03@

with e the vector of residuals scaled by their maximum values, n the number of measurements and p the number of parameters to estimate.

The confidence bounds of the parameters were estimated by the approximation of the Hessian matrix, H, by the Jacobian matrix, J, at the minimum root mean squared error, rmse.

H = JT·J·rmse (7)

Finally, confidence bounds, CBp, were obtained by the estimate of standard deviations for a level of confidence of 95%:

CB = diag(H)·t(1 - 0.95, n - p) (8)

with t the t-student distribution.

For the calculation of residuals, Eqs. (4a–g) were integrated using a 4th/5th order Runge-Kutta solver (MATLAB's ode45 function), with the obtained concentrations of acetate, propionate, PHB, PHV, PH2MV and oxygen subtracted from the respective measured values to give the residuals.

Metabolic Flux Analysis

Metabolic Flux Analysis (MFA) is a methodology that allows the determination of a set of unknown metabolic fluxes, v n , from a set of known fluxes, v b [29] by steady-state material balancing (here using Eqs. (1a–f)). In the present case, v b is composed of the modelled fluxes:

v b = [R4, R5, R6, R7, R9, R10]T

The remaining fluxes were calculated by solving the central MFA equation [29]:

v n = -A n -1·A b ·v b

with v n the vector of unknown fluxes,

v n = [R1, R2, R3, R8, R PHB , R PHV , RPH 2MV]T

A n and A b are the corresponding stoichiometric matrices of the unknown and known fluxes respectively. Note that in this case A n is a square (7 × 7) matrix.

For single substrate feeding studies, either R1 or R2 equals zero, thus only 5 modelled

fluxes had to be included in v b , being R6 moved to v n .

Flux Balance Analysis

Flux Balance Analysis (FBA) is a method to optimize a given metabolic objective function under the constraints of steady-state material balancing of intracellular intermediates and other known biological constraints [30]. FBA was applied in this study to maximize the flux of PHA synthesis with a desired 24% (C-mol/C-mol) propionyl-CoA* composition. Mathematically, the problem can be stated in the following way.

max R 1 , R 2 ( R P H B + R P H V + R P H 2 M V ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaWaaCbeaeaacyGGTbqBcqGGHbqycqGG4baEaSqaaiabdkfasnaaBaaameaacqaIXaqmaeqaaSGaeiilaWIaemOuai1aaSbaaWqaaiabikdaYaqabaaaleqaaOWaaeWaaeaacqWGsbGudaWgaaWcbaGaemiuaaLaemisaGKaemOqaieabeaakiabgUcaRiabdkfasnaaBaaaleaacqWGqbaucqWGibascqWGwbGvaeqaaOGaey4kaSIaemOuai1aaSbaaSqaaiabdcfaqjabdIeaijabikdaYiabd2eanjabdAfawbqabaaakiaawIcacaGLPaaaaaa@49E0@

under the following constraints:

0 = A n ·v n + A b ·v b (12b)

0 = R PHB +R PHV + RPH 2MV- 0.76·R9 - 0.24·R10 (12c)

R4 = R5 = 0 (12d)

R7 = 0.02 (12e)

R 3 R 2 = 0.43 R 2 R 1 + R 2 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqWGsbGudaWgaaqaaiabiodaZaqabaaabaGaemOuai1aaSbaaeaacqaIYaGmaeqaaaaakiabg2da9iabicdaWiabc6caUiabisda0iabiodaZKqbaoaalaaabaGaemOuai1aaSbaaeaacqaIYaGmaeqaaaqaaiabdkfasnaaBaaabaGaeGymaedabeaacqGHRaWkcqWGsbGudaWgaaqaaiabikdaYaqabaaaaaaa@3E4D@

Constraint (12c) defines the desired monomeric composition of the final polymer. Constraint (12d) states that the cells are not growing during the PHA synthesis phase. The maintenance flux (Eq. 12e) was set to a known flux (see Table 4). Constraint (12f) was established on the basis of experimental observations and discussed later in the results section. FBA problems are normally solved through linear programming using the simplex algorithm. In our case, constraint (12f) is a nonlinear equation, thus a nonlinear function solver was adopted. The MATLAB function 'fmincon' based on the quasi-Newton optimization algorithm was adopted to solve the system (12a–f).

Table 4 Parameters estimation results of acetate and propionate enriched cultures for the different feeding conditions

Results and discussion

Estimation of the kinetic parameters

The data of 18 batch experiments were used for model calibration and validation. In mixed microbial cultures there is a diversity of organisms with different phenotypes competing for the nutrients in the medium. The observed metabolic activity is thus an "average" metabolic activity over all organisms present in the culture. This metabolic activity is obviously highly dependent on the culture enrichment regimen and also on the VFA feeding strategy. As such, the parameter estimation was done separately for the following 6 distinct experimental conditions:

  1. 1.

    Acetate enriched culture fed with acetate.

  2. 2.

    Acetate enriched culture fed with propionate.

  3. 3.

    Acetate enriched culture fed with mixtures of acetate and propionate.

  4. 4.

    Propionate enriched culture fed with acetate.

  5. 5.

    Propionate enriched culture fed with propionate.

  6. 6.

    Propionate enriched culture fed with mixtures of acetate and propionate.

The estimated parameter values and respective 95% confidence limits are compiled in Table 4. Note that the estimated confidence bounds are generally quite low, denoting the high sensitivity of residuals to parameters and the high statistical confidence of the estimated parameter values.

Figures 2, 3, 4, 5, 6, 7 compare model predictions and respective measurements for the data of the 18 experiments used in this study. The modelling results for acetate enriched cultures fed with acetate (7 calibration experiments and 3 validation experiments) and propionate enriched culture fed with propionate (3 calibration experiments and 1 validation experiment) are shown in Figures 2 and 3, respectively. Full symbols represent calibration points while open symbols denote validation points. These results show that the model was able to consistently describe the experimental data for a wide range of acetate and propionate feeding conditions, with regression coefficients, r2, always above 0.94. Figures 4, 5, 6, 7 show the modelling results obtained for the remaining experiments, where both cultures were fed with substrates other than those used in the enrichment phase. In these Figures, the experimental data of the acetate, propionate, PHB, PHV, PH2MV and oxygen is represented by full symbols and the model prediction by full lines. A good agreement between experimental data and model predictions was obtained in all cases.

Figure 2
figure 2

Modelling results. Figure shows the results obtained for acetate enriched cultures under acetate feeding conditions. Full symbols represent parameter estimation results and open symbols the model validation results. Dashed lines are the 95% confidence limits for the parameter estimation and model validation results.

Figure 3
figure 3

Modelling results. Figure shows the results obtained for propionate enriched cultures under propionate feeding conditions. Full symbols represent parameter estimation results and open symbols the model validation results. Dashed lines are the 95% confidence limits for the parameter estimation and model validation results.

Figure 4
figure 4

Modelling results. Figure shows the results obtained for acetate enriched cultures under propionate feeding conditions. Full symbols represent the experimental data and full lines the modelling results.

Figure 5
figure 5

Modelling results. Figure shows the results obtained for acetate enriched cultures under acetate and propionate feeding conditions. Full symbols represent the experimental data and full lines the modelling results.

Figure 6
figure 6

Modelling results. Figure shows the results obtained for propionate enriched cultures under acetate feeding conditions. Full symbols represent the experimental data and full lines the modelling results.

Figure 7
figure 7

Modelling results. Figure shows the results obtained for propionate enriched cultures under acetate and propionate feeding conditions. Full symbols represent the experimental data and full lines the modelling results.

P/O ratio

The P/O ratio is a measure of the efficiency of ATP synthesis coupled to cell respiration, indicative of the efficiency of catabolism. In acetate enriched cultures, the P/O ratio was significantly higher than in propionate enriched cultures (Table 4). Moreover, the acetate feast and famine enrichment strategy yielded an extremely energetically effective population with a P/O ratio close to the theoretical maximum, which is between 2–3 mol-ATP/mol-NADH2 in bacteria [9, 10, 3133]. This result may suggest that the feast and famine strategy may induce not only the accumulation of intracellular reserves but also the optimization of the global energetic efficiency of the final selected culture.

The P/O ratio results in Table 4 also show that after a short term swap in substrate, the P/O ratio is highly affected. This behaviour is common to both cultures used in the present study. When propionate is the main source of energy, the P/O ratio is consistently much lower (~1 in both acetate enriched and propionate enriched cultures) than when acetate is the main source of energy (2.9 and 1.8 in acetate and propionate enriched cultures, respectively). This result is also consistent with the generally lower energetic efficiency associated with the propionate enriched culture. Interestingly, when mixtures of acetate and propionate are fed to the culture, the characteristic P/O ratio of acetate prevails over that of propionate, which suggests that catabolic activity and respiration is preferentially executed through acetate metabolism (see discussion below for more detail).

A detailed study about the P/O ratio for microbial cells using different substrates was performed by Stouthamer [34]. In this study, a range of P/O ratios between 2.25 and 3 were found for different substrates, whereas 2.25 was the value calculated for acetate. On the other hand, Gottschalk [19] stated that P/O ratio may also vary with bacterial species and more specifically with the number of phospholylation sites they contain. In mixed cultures containing enrichments of polyphosphate and glycogen-accumulating organisms (PAO and GAO), different P/O ratios have been reported depending on the carbon source used for culture enrichment, namely acetate and propionate. The P/O ratio reported for acetate enriched cultures was 1.85 for PAO [35] and 1.73 for GAO [22], while in propionate enriched cultures the P/O ratio was reported as 1.37 for PAO and 1.29 for GAO [36]. Furthermore, it was observed that the mixed microbial culture fed with acetate exhibited a significantly lower P/O ratio (1.38) at the beginning of the reactor operational period, before the acetate culture had adapted to a feast and famine regimen (data not shown). A very high P/O ratio (2.88) was also observed for a pure culture of Cupriavidus necator (formerly Ralstonia eutropha) when imposed to the feast and famine regimen for PHA production [33]. These statements support the hypothesis that the carbon source affects the P/O ratio and the feast and famine regimen may enable the selection of bacteria with a higher P/O ratio, mainly in acetate enriched cultures.

Higher P/O ratios are normally associated with lower maintenance coefficients and higher biomass/substrate and product/substrate yields. Table 5 compiles the theoretical yields for the different culture enrichment scenarios. As expected, for the same feeding conditions a higher P/O ratio led to higher yields and lower maintenance coefficients. Within this frame, it is clear that the acetate feast and famine culture enrichment strategy produces the highest yields of cell growth and PHA formation.

Table 5 Yields and maintenance coefficients for acetate and propionate enriched cultures for the different feeding conditions

Metabolic flux distribution

The calculated metabolic flux distribution (MFD) was obtained in all cases in conditions of excess of carbon source material and negligible cell growth conditions (low ammonia concentration in the medium). The results are shown in Figures 8a–d. When acetate is fed as the sole carbon source to acetate or propionate enriched cultures, the overall metabolic activity is significantly higher in the former case than in the latter case (Figure 8a). All fluxes are at least two-threefold higher in the acetate enriched culture as compared to the propionate enriched culture.

Figure 8
figure 8

MFD results. Figure shows MFD results for acetate (a), propionate (b) and acetate + propionate (c) feeding using acetate (■) and propionate (□) enriched cultures. Figure 8d presents MFD results of the main substrates uptake, PHA and oxidative phosphorylation fluxes for the two enriched cultures.

On the other hand, if propionate is fed to both cultures, all fluxes, without exception, are higher in the propionate enriched culture when compared to the acetate culture (Figure 8b). These results confirm that the substrate used for culture enrichment is more effectively metabolized by the selected culture. However, the differences between the fluxes of both cultures in Figure (8b) are markedly lower than those of Figure (8a). In particular, the total carbon uptake and corresponding PHA storage are much higher in the acetate enriched culture fed with acetate. It seems that acetate was a much more effective substrate than propionate for selective enrichment of PHA producing cultures, leading to higher quantities of PHA stored.

The simultaneous feeding of acetate and propionate to the enriched cultures is shown in Figure 8c. Acetate and propionate are taken up at similar fluxes in each culture: 0.33 C-mol/(C-mol.h) of acetate and 0.34 C-mol/(C-mol.h) of propionate was taken up in acetate enriched cultures and 0.22 C-mol/(C-mol.h) of acetate and 0.23 C-mol/(C-mol.h) of propionate was taken up in propionate enriched cultures. The total carbon flux is 33% higher in acetate enriched cultures when compared to the propionate enriched cultures. The acetate enriched culture produced a total PHA flux of 0.52 C-mol/(C-mol.h), whereas the PHA flux in the propionate enriched culture is 0.31 C-mol/(C-mol.h) (40% lower). However, the catabolism flux, R6, is about the same in both cultures, suggesting that the percentage of carbon source spent for maintenance is higher for propionate enriched cultures due to the much lower P/O ratio in propionate enriched cultures. This is further confirmed by the higher oxidative phosphorylation flux, R8, for the propionate enriched culture.

It is also clear from Figures 8b–c that PH2MV was produced in the acetate culture and not in the propionate culture, due to the fact that acetyl-CoA* and propionyl-CoA* tended to condense randomly in the acetate culture and selectively in the propionate culture. A possible explanation for this result is that selective condensation of propionyl-CoA* with acetyl-CoA* may be a property of cells acclimatized to a propionate substrate, and that the acetate culture, which was never previously exposed to propionate, may have lacked the necessary enzymes or microbial population that is responsible for selective condensation.

Figure 8d shows a comparison of the fluxes for some of the main reactions in the two enriched cultures during the different feeding conditions. In acetate enriched cultures, the total VFA uptake capacity obtained with acetate feeding (0.73 C-mol/(C-mol.h)) was approximately equal to the sum of acetate and propionate uptake, when the substrates were fed simultaneously. Propionate enriched cultures showed a completely different behaviour. The total carbon uptake observed when both substrates were fed to this culture was approximately double that of the propionate or acetate uptake when the substrates were fed individually. A similar trend can be observed for the oxygen consumption flux and total PHA production flux, where the total flux in the propionate culture is approximately equal to the individual acetate and propionate fluxes combined, while the mixed substrate feed exhibited similar fluxes to the case of acetate feeding in the acetate culture. The mixed substrate feed had a synergistic effect on the individual metabolism of both substrates: both studied cultures were more able to take up propionate and convert it into PHA with a combined acetate/propionate feeding. This increase in VFA uptake was also observed in a culture adapted to a mixture of acetate, propionate and lactic acid [16]. A possible explanation for this result is that the conversion of propionyl-CoA to acetyl-CoA is the rate limiting step for PHV synthesis using only propionate. The effect of this limitation on VFA uptake rate is attenuated when acetate is also fed, because the requirement of propionyl-CoA driven to acetyl-CoA is lower and more propionyl-CoA can be driven for propionyl-CoA* synthesis without losing carbon in reaction, R3.

Since the fraction of propionate per total VFA in the feed had a substantial influence on the amount of propionate uptake driven to acetyl-CoA, a relationship was established based on these two fractions, for each mixed culture under either a propionate or a mixed propionate-acetate feed (see Eq. 12f). The linear regression coefficient, r2, for this equation was found to be 0.98, based on the data presented in Table 4. It is clear that an increase in the propionate fraction fed to the mixed cultures leads to an increase in the amount of propionate converted through propionyl-CoA to acetyl-CoA. Thus, there was a higher requirement for acetyl-CoA production when the relative acetate uptake rate is lower. This supports the hypothesis discussed above, where it was proposed that energy generation through the catabolic activity of the cells are preferentially executed through acetyl-CoA.

Overall, the results from the tests with the simultaneous feeding of acetate and propionate again corroborate the hypothesis that acetate was a much more effective substrate than propionate for the selection of an optimal PHA producing culture through enrichment via the feast and famine regimen.

Flux balance analysis

From the results presented above, it is clear that acetate is a superior substrate to propionate for the enrichment of a microbial culture performing PHA production under the feast and famine regimen, due to the higher oxidative phosphorylation efficiency with this carbon source, resulting in higher PHA productivity. It is also clear that the feeding of both acetate and propionate carbon sources simultaneously resulted in either a similar or higher total PHA productivity, with a higher diversity of 3HA monomers as compared to the feeding of a single substrate. Thus, in PHA producing systems where the microbial culture enrichment and PHA production phases are separated, it is desirable to enrich the microbial culture using acetate, while feeding a combination of acetate and propionate in the PHA production phase.

During the PHA production phase, it is possible to control the final PHA polymer composition through manipulating the acetate and propionate feeding fractions, using the metabolic model presented in this study. It was stated that a PHA copolymer with 15–20 mol % of propionyl-CoA* (18–24 C-mol %) is desired for improved toughness of the final polymer [18]. Assuming this propionyl-CoA* content as a constraint, the optimal feeding strategy required to obtain this copolymer was defined by performing FBA on the acetate enriched culture, where Eq. 12f was also used as a constraint. The FBA results are presented in Figure 9. These results show that the maximum PHA production rate with the desired propionyl-CoA* content is obtained when the propionate uptake rate is below the maximum calculated from parameter estimation (0.34 C-mol/(C-mol.h)). The limitation of the propionate uptake rate can be achieved by ensuring that the feeding rate of propionate is equal to the desired propionate uptake rate during reactor operation. When the process was operated at the maximum propionate uptake rate, the obtained propionyl-CoA* content in the PHA polymer was about 42%.

Figure 9
figure 9

FBA results. Figure shows FBA results obtained by maximizing the PHA flux (RPHA*) with a propionyl-CoA* content of 24% (C-mol basis).

Metabolic comparison of PHA production by pure and mixed cultures

The process efficiency for PHA production by mixed microbial cultures was evaluated through a comparison of the metabolic fluxes for the PHA production process by pure cultures of microorganisms using the same carbon substrates (i.e. acetate and propionate). In the present study using mixed microbial cultures, the biomass growth was assumed to be negligible based on the very low feed concentration of ammonia. In the pure culture studies [18, 37], the maximum carbon fraction used for biomass growth was less than 15%, providing a good basis for comparison with the present study. MFA was performed for a pure culture of Cupriavidus necator using acetate as the sole carbon source [37]. As shown in Table 6, the acetate uptake rate is higher for mixed microbial cultures enriched on acetate (0.73 C-mol/(C-mol.h)) than for pure cultures (0.44 C-mol/(C-mol.h)). However, the acetate uptake rate in the cultures enriched on propionate is half of the rate observed in pure cultures (0.21 C-mol/(C-mol.h)). The fluxes for CO2 formation and O2 consumption are slightly higher for pure cultures. These differences could be explained by the additional energetic consumption for cell growth in pure cultures that is negligible in both mixed cultures. The fraction of carbon driven for PHB is significantly higher for both mixed cultures. This is likely due to the very high P/O ratio observed, especially in the mixed culture fed with acetate and, hence, a more efficient utilization of the carbon supply is expected.

Table 6 Fluxes distribution for acetate and acetate + propionate feeding between mixed microbial cultures and Cupriavidus necator

For the case of a Cupriavidus necator culture fed with a mixture of acetate and propionate [18], the VFA uptake kinetics were observed to be different from the mixed microbial cultures enriched with either acetate or propionate, despite the fact that the initial fraction of acetate and propionate were the same in all cases (50:50 on a C-mol basis). The VFA uptake rate observed in the acetate enriched culture is much higher than in pure cultures (about twofold), whereas, the propionate enriched culture has only a slightly higher VFA uptake rate as compared to the pure cultures. The PHA production fluxes were observed to be much higher for the cases of the mixed cultures (78.0 and 68.6%), as compared to the pure culture (41.2%), during the feeding of a mixture of acetate and propionate. Additionally, despite the fact that the relative propionate uptake rate is higher in pure cultures (63%) than in mixed cultures (51%), the fraction of propionyl-CoA driven to propionyl-CoA* is significantly lower (86% lower). This is due to the fact that a large portion of the propionate fed to Cupriavidus necator was first converted to acetyl-CoA prior to PHA production, unlike the mixed culture cases. This metabolic pathway causes a large loss of CO2 by the pure culture, thus lowering the PHA production efficiency. Furthermore, a large fraction of this propionyl-CoA converted to acetyl-CoA was then used for cell catabolism, contributing to the much higher carbon flux observed in the case of the pure culture as compared to the mixed cultures.

It should be noted that the carbon source used in pure culture studies [18, 37] was neither acetate nor propionate, but mixtures of yeast extract, meat extract and peptone. This could be one explanation for the inferior VFA uptake and PHA production results observed in these studies as compared to the present study, since the culture may require an adaptation period to synthesize the necessary enzymes for the metabolism of acetate and propionate carbon sources. Alternatively, it may be that the microorganisms present in the mixed cultures are more efficient PHA producers when compared to the pure cultures of Cupriavidus necator, due to the feast and famine regimen imposed. An increased VFA uptake rate was observed for a pure culture of Amaricoccus kaplicensis when the length of the famine period was increased [38].


• This paper presents a metabolic model of PHA copolymers production in mixed microbial cultures. The model was applied to two different cultures, obtained through two distinct enrichment protocols (selected with either acetate or propionate), under different feeding conditions (fed with either a single substrate or with mixtures of the two substrates). With this model, intracellular flux distributions for the different cases were calculated. Also, feeding scenarios were optimized by FBA targeting maximal productivity of a copolymer with a desired monomeric composition. These results were benchmarked with published results of pure cultures of Cupriavidus necator. From these studies the following main conclusions may be highlighted:

• The substrate used for culture enrichment by feast and famine feeding has a high selective pressure on the organisms that compose the final culture in the sense that the organisms selected are those that most effectively metabolize the adopted substrate.

• Sludge enrichment by acetate feeding through feast and famine regimen is much more selective towards a culture with high PHA storage fluxes and yields than when propionate is used.

• The P/O ratio is highly dependent on the substrate and the microbial culture selected. Acetate metabolism has a consistently higher P/O ratio (close to the theoretical maximum of 3 mol-ATP/mol-NADH2) than propionate metabolism.

• MFA studies suggest that when mixtures of acetate and propionate are fed to the cultures, the catabolic activity is primarily guaranteed through acetate uptake. Coherently, the acetate P/O ratio tends to prevail over that of propionate when both substrates are fed simultaneously.

• The application of FBA targeting the optimization of the PHA synthesis flux with a 24% (C-mol/C-mol) propionyl-CoA* composition revealed that acetate should be fed in excess whereas the feeding of propionate should be limited to ~0.17 C-mol/(C-mol.h). Under excess of both substrates the final propionyl-CoA* content is around 42% (C-mol/C-mol).

• By comparing metabolic fluxes between mixed and pure cultures of Cupriavidus necator it can be concluded that the PHA production process by mixed microbial cultures has the potential to be comparable or even more favourable than what is achieved by pure cultures.

• Although the metabolic characterization of mixed cultures was shown in this study to be highly favourable and promising, it should be noted that these results were obtained under ammonia limitation, thus, with conditions where cells growth is negligible. The compatibility of high cell growth rates with high PHA synthesis fluxes in mixed cultures is still not undoubtedly demonstrated in the literature.


Acronyms :

ASM3: activated sludge model No. 3

DO :

dissolved oxygen


flux balance analysis


glycogen-accumulating organisms

3HA :


3HB :


3H2MB :


3H2MV :


3HV :



hydraulic retention time


metabolic flux analysis


metabolic flux distribution


polyphosphate-accumulating organisms









P/O ratio :

ATP synthesis/oxygen consumption ratio

rmse :

root mean squared error


sequencing batch reactor


sludge retention time


tricarboxylic acid cycle


volatile fatty acids


volatile suspended solids

Symbols: Ac:

acetate/acetate concentration (C-mmol/l)

f i :

intracellular PHA, PHB, PHV and PH2MV contents (C-mol/C-mol)

K1, K2:

energy requirements for biomass synthesis (mol-ATP/C-mol)

K S :

acetate and propionate half-saturation constants (C-mmol/l)

K N :

ammonia half-saturation constant (N-mmol/l)

m j :

maintenance on component j [C-mol/(C-mol.h)]

m ATP :

maintenance coefficient on ATP [mol-ATP/(C-mol.h)]

N :

ammonia/ammonia concentration (N-mmol/l)

O 2 :

oxygen/oxygen concentration (mmol/l)

Prop :

propionate/propionate concentration (C-mmol/l)

R i :

specific rate of reaction on compound i [C-mol/(C-mol.h)]

S :

total VFA/total VFA concentration (C-mmol/l)

t :

culture runtime (h)

X :

active biomass/active biomass concentration (C-mmol/l)

y :

fraction of propionate uptake rate in the total VFA uptake rate

Y i/j :

yield of component i on component j (C-mol/C-mol)

α :

PHB production saturation order constant (dimensionless)

δ :

efficiency of oxidative phosphorylation (mol-ATP/mol-NADH2).


  1. Dias JML, Lemos PC, Serafim LS, Oliveira C, Eiroa M, Albuquerque MGE, Ramos AM, Oliveira R, Reis MAM: Recent advances in polyhydroxyalkanoate production by mixed aerobic cultures: From the substrate to the final product. Macromolecular Bioscience. 2006, 6 (11): 885-906. 10.1002/mabi.200600112

    Article  CAS  PubMed  Google Scholar 

  2. Beccari M, Majone M, Massanisso P, Ramadori R: A bulking sludge with high storage response selected under intermittent feeding. Water Research. 1998, 32 (11): 3403-3413. 10.1016/S0043-1354(98)00100-6.

    Article  CAS  Google Scholar 

  3. Satoh H, Iwamoto Y, Mino T, Matsuo T: Activated sludge as a possible source of biodegradable plastic. Water Science and Technology. 1998, 38 (2): 103-109. 10.1016/S0273-1223(98)00435-1.

    Article  CAS  Google Scholar 

  4. Serafim LS, Lemos PC, Oliveira R, Reis MAM: Optimization of polyhydroxybutyrate production by mixed cultures submitted to aerobic dynamic feeding conditions. Biotechnol Bioeng. 2004, 87 (2): 145-160. 10.1002/bit.20085

    Article  CAS  PubMed  Google Scholar 

  5. Albuquerque MGE, Eiroa M, Torres C, Nunes BR, Reis MAM: Strategies for the development of a side stream process for polyhydroxyalkanoate (PHA) production from sugar cane molasses. Journal of Biotechnology. 2007, 130 (4): 411-421. 10.1016/j.jbiotec.2007.05.011

    Article  CAS  PubMed  Google Scholar 

  6. Bengtsson S, Werker A, Christensson M, Welander T: Production of polyhydroxyalkanoates by activated sludge treating a paper mill wastewater. Bioresour Technol. 2008, 99 (3): 509-516. 10.1016/j.biortech.2007.01.020

    Article  CAS  PubMed  Google Scholar 

  7. Temudo MF, Kleerebezem R, van Loosdrecht M: Influence of the pH on (open) mixed culture fermentation of glucose: A chemostat study. Biotechnol Bioeng. 2007, 98 (1): 69-79. 10.1002/bit.21412

    Article  CAS  PubMed  Google Scholar 

  8. van Aalst-van Leeuwen MA, Pot MA, van Loosdrecht MCM, Heijnen JJ: Kinetic modeling of poly(beta-hydroxybutyrate) production and consumption by Paracoccus pantotrophus under dynamic substrate supply. Biotechnol Bioeng. 1997, 55 (5): 773-782. 10.1002/(SICI)1097-0290(19970905)55:5<773::AID-BIT7>3.0.CO;2-8.

    Article  CAS  PubMed  Google Scholar 

  9. Third KA, Newland M, Cord-Ruwisch R: The effect of dissolved oxygen on PHB accumulation in activated sludge cultures. Biotechnol Bioeng. 2003, 82 (2): 238-250. 10.1002/bit.10564

    Article  CAS  PubMed  Google Scholar 

  10. Dias JML, Serafim LS, Lemos PC, Reis MAM, Oliveira R: Mathematical modelling of a mixed culture cultivation process for the production of polyhydroxybutyrate. Biotechnol Bioeng. 2005, 92 (2): 209-222. 10.1002/bit.20598

    Article  CAS  PubMed  Google Scholar 

  11. Lemos PC, Serafim LS, Reis MAM: Synthesis of polyhydroxyalkanoates from different short-chain fatty acids by mixed cultures submitted to aerobic dynamic feeding. Journal of Biotechnology. 2006, 122 (2): 226-238. 10.1016/j.jbiotec.2005.09.006

    Article  CAS  PubMed  Google Scholar 

  12. Lafferty M, Korsatko B, Korsatko W: Microbial Production of Poly-β-hydroxybutyric acid. 1988, 6b: Weinheim: VCH Publishers

    Google Scholar 

  13. Anderson AJ, Dawes EA: Occurrence, Metabolism, Metabolic Role, and Industrial Uses of Bacterial Polyhydroxyalkanoates. Microbiological Reviews. 1990, 54 (4): 450-472.

    PubMed Central  CAS  PubMed  Google Scholar 

  14. Lee SY: Plastic bacteria? Progress and prospects for polyhydroxyalkanoate production in bacteria. Trends in Biotechnology. 1996, 14 (11): 431-438. 10.1016/0167-7799(96)10061-5.

    Article  CAS  Google Scholar 

  15. Rhu DH, Lee WH, Kim JY, Choi E: Polyhydroxyalkanoate (PHA) production from waste. Water Sci Technol. 2003, 48 (8): 221-228.

    CAS  PubMed  Google Scholar 

  16. Dionisi D, Majone M, Papa V, Beccari M: Biodegradable polymers from organic acids by using activated sludge enriched by aerobic periodic feeding. Biotechnol Bioeng. 2004, 85 (6): 569-579. 10.1002/bit.10910

    Article  CAS  PubMed  Google Scholar 

  17. Dionisi D, Carucci G, Papini MP, Riccardi C, Majone M, Carrasco F: Olive oil mill effluents as a feedstock for production of biodegradable polymers. Water Research. 2005, 39 (10): 2076-2084. 10.1016/j.watres.2005.03.011

    Article  CAS  PubMed  Google Scholar 

  18. Yu J, Si YT: Metabolic carbon fluxes and biosynthesis of polyhydroxyalkanoates in Ralstonia eutropha on short chain fatty acids. Biotechnol Prog. 2004, 20 (4): 1015-1024. 10.1021/bp034380e

    Article  CAS  PubMed  Google Scholar 

  19. Gottschalk : Bacterial metabolism. 1986, New York: Springer-Verlag, Second

    Chapter  Google Scholar 

  20. Filipe CDM, Daigger GT, Grady CPL: A metabolic model for acetate uptake under anaerobic conditions by glycogen accumulating organisms: Stoichiometry, kinetics, and the effect of pH. Biotechnol Bioeng. 2001, 76 (1): 17-31. 10.1002/bit.1022

    Article  CAS  PubMed  Google Scholar 

  21. Gommers PJF, Vanschie BJ, Vandijken JP, Kuenen JG: Biochemical Limits to Microbial-Growth Yields – an Analysis of Mixed Substrate Utilization. Biotechnol Bioeng. 1988, 32 (1): 86-94. 10.1002/bit.260320112

    Article  CAS  PubMed  Google Scholar 

  22. Zeng RJ, van Loosdrecht MCM, Yuan ZG, Keller J: Metabolic model for glycogen-accumulating organisms in anaerobic/aerobic activated sludge systems. Biotechnol Bioeng. 2003, 81 (1): 92-105. 10.1002/bit.10455

    Article  CAS  PubMed  Google Scholar 

  23. Henze M, Harremoes P, Jansen JLC, Arvin E: Wastewater Treatment: Biological and Chemical Processes. 1995, Heidelberg: Springer

    Google Scholar 

  24. Roels JA: Macroscopic Thermodynamics and the Description of Growth and Product Formation in Microorganisms. Acs Symposium Series. 1983, 207: 295-322.

    Article  CAS  Google Scholar 

  25. Lee JH, Lim HC, Hong J: Application of nonsingular transformation to online optimal control of poly-beta-hydroxybutyrate fermentation. Journal of Biotechnology. 1997, 55 (3): 135-150. 10.1016/S0168-1656(97)00064-3.

    Article  CAS  Google Scholar 

  26. Lee JH, Hong J, Lim HC: Experimental optimization of fed-batch culture for poly-beta-hydroxybutyric acid production. Biotechnol Bioeng. 1997, 56 (6): 697-705. 10.1002/(SICI)1097-0290(19971220)56:6<697::AID-BIT13>3.0.CO;2-5.

    Article  CAS  PubMed  Google Scholar 

  27. Tohyama M, Patarinska T, Qiang ZW, Shimizu K: Modeling of the mixed culture and periodic control for PHB production. Biochemical Engineering Journal. 2002, 10 (3): 157-173. 10.1016/S1369-703X(01)00184-X.

    Article  CAS  Google Scholar 

  28. , : Standard methods for the examination of water and wastewater. 1989, Washington DC: American Public Health Association

    Google Scholar 

  29. Stephanopoulos GN, Aristidou AA, Nielsen J: Metabolic Engineering: Principles and Methodologies. 1998, San Diego: Academic Press

    Google Scholar 

  30. Edwards JS, Ibarra RU, Palsson BO: In silico predictions of Escherichia coli metabolic capabilities are consistent with experimental data. Nat Biotechnol. 2001, 19 (2): 125-130. 10.1038/84379

    Article  CAS  PubMed  Google Scholar 

  31. Beun JJ, Paletta F, Van Loosdrecht MCM, Heijnen JJ: Stoichiometry and kinetics of poly-beta-hydroxybutyrate metabolism in aerobic, slow growing, activated sludge cultures. Biotechnol Bioeng. 2000, 67 (4): 379-389. 10.1002/(SICI)1097-0290(20000220)67:4<379::AID-BIT1>3.0.CO;2-2

    Article  CAS  PubMed  Google Scholar 

  32. Sin G, Guisasola A, De Pauw DJW, Baeza JA, Carrera J, Vanrolleghem PA: A new approach for modelling simultaneous storage and growth processes for activated sludge systems under aerobic conditions. Biotechnol Bioeng. 2005, 92 (5): 600-613. 10.1002/bit.20741

    Article  CAS  PubMed  Google Scholar 

  33. Frigon D, Muyzer G, van Loosdrecht M, Raskin L: rRNA and poly-beta-hydroxybutyrate dynamics in bioreactors subjected to feast and famine cycles. Appl Environ Microbiol. 2006, 72 (4): 2322-2330. 10.1128/AEM.72.4.2322-2330.2006

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Stouthamer AH: Theoretical Study on Amount of ATP Required for Synthesis of Microbial Cell Material. Antonie Van Leeuwenhoek Journal of Microbiology. 1973, 39 (3): 545-565. 10.1007/BF02578899.

    Article  CAS  Google Scholar 

  35. Smolders GJF, Vandermeij J, Vanloosdrecht MCM, Heijnen JJ: Stoichiometric Model of the Aerobic Metabolism of the Biological Phosphorus Removal Process. Biotechnol Bioeng. 1994, 44 (7): 837-848. 10.1002/bit.260440709.

    Article  CAS  PubMed  Google Scholar 

  36. Oehmen A, Zeng RJ, Keller J, Yuan Z: Modeling the Aerobic Metabolism of Polyphosphate-Accumulating Organisms Enriched with Propionate as a Carbon Source. Water Environment Research. 2007, 79 (13): 2477-2486.

    CAS  PubMed  Google Scholar 

  37. Shi HD, Shiraishi M, Shimizu K: Metabolic flux analysis for biosynthesis of poly(beta-hydroxybutyric acid) in Alcaligenes eutrophus from various carbon sources. J Ferment Bioeng. 1997, 84 (6): 579-587. 10.1016/S0922-338X(97)81915-0.

    Article  CAS  Google Scholar 

  38. Aulenta F, Dionisi D, Majone M, Parisi A, Ramadori R, Tandoi V: Effect of periodic feeding in sequencing batch reactor on substrate uptake and storage rates by a pure culture of Amaricoccus kaplicensis. Water Research. 2003, 37 (11): 2764-2772. 10.1016/S0043-1354(03)00059-9

    Article  CAS  PubMed  Google Scholar 

  39. Beun JJ, Dircks K, Van Loosdrecht MCM, Heijnen JJ: Poly-beta-hydroxybutyrate metabolism in dynamically fed mixed microbial cultures. Water Research. 2002, 36 (5): 1167-1180. 10.1016/S0043-1354(01)00317-7

    Article  CAS  PubMed  Google Scholar 

  40. Gujer W, Henze M, Mino T, van Loosdrecht M: Activated Sludge Model No. 3. Water Science and Technology. 1999, 39 (1): 183-193. 10.1016/S0273-1223(98)00785-9.

    Article  CAS  Google Scholar 

Download references


This work was supported by Fundação para a Ciência e a Tecnologia (FCT) through the project POCI/BIO/55789/2004 and by the Integrated Project no. 026515-2: Bioproduction – Sustainable Microbial and Biocatalytic Production of Advanced Functional Materials. J. Dias, A. Oehmen and L. Serafim acknowledge FCT for grants SFRH/BD/13714/2003, SFRH/BPD/41486/2007 and SFRH/BPD/14663/2003.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Rui Oliveira.

Additional information

Authors' contributions

All authors read and approved the final manuscript. JMLD and RO developed the software. JMLD, AO and RO participated in the model implementation. LSS and PCL performed experimental work. MAMR and RO designed and coordinated the study. JMLD, AO and RO drafted the manuscript.

João ML Dias, Adrian Oehmen, Luísa S Serafim, Paulo C Lemos, Maria AM Reis contributed equally to this work.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Dias, J.M., Oehmen, A., Serafim, L.S. et al. Metabolic modelling of polyhydroxyalkanoate copolymers production by mixed microbial cultures. BMC Syst Biol 2, 59 (2008).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: