Study1 – extracellular insoluble light (EX NSL LT) products
For the EX NSL LT category, we choose an EX NSL LT LQD CMD product as the representative base case. Note that EX NSL LT products float on the top and are thus naturally separated from the dense microbial cells. The key parameters used are shown in Table 2. The base case parameters have been chosen to represent typical cases reported in the literature for products belonging to EX NSL LT class such as terpenes [58, 59, 96]. The analysis range for the product titer has been chosen up to 250 g/L to reflect the possibilities of achieving higher product concentrations due to advances in metabolic engineering and bioreactor designs [97,98,99].
Superstructure and optimal solution
By simplifying the separation scheme of an EX NSL product (provided in Additional file 1: Figure S. A1), we obtain the superstructure for EX NSL LT LQD CMD product shown in Fig. 2. In Stage I – phase isolation, Sdm1 (sedimentation), and Cnt1 (centrifugation) separate the product as a top phase (isolated from the cells at the bottom) while removing water at the same time. Ftt (filtration) functions to only remove the cells. In Stage II – concentration, Sdm2, Cnt2 and Mbr (membranes) can separate the product from water, which can be bypassed if enough concentration has been achieved in Stage I. In Stage III – refinement, Ads (adsorption) further purifies the product, which can also be bypassed if the optimization model decides that the final product specifications have been satisfied. The optimization model comprises 601 equations, 486 continuous variables and 33 discrete variables and takes about a minute or two to solve for one set of parameter values.
After solving the superstructure optimization model, we obtain the cost-minimal separation network (Cnt1 in Stage I followed by Cnt2 in Stage II), as represented by the highlighted parts in Fig. 2, as well as the corresponding cost distribution. The total minimum cost is 0.65 $/kg, with the feed accounting for 60% of the total cost, Stage I accounting for 38%, and Stages II for 2% of the total costs. Cnt2 is also selected because Cnt1 alone is not able to concentrate the product stream from 5 g/L (~ 0.5 wt%) to the required final purity of 95 wt% due to concentrating factor limitations (see centrifugation (Cnt) parameters in Additional File).
Analysis for EX NSL LT LQD CMD product
Performance of phase isolation technologies
Since the major cost component in the optimal network is Stage I, we vary parameters related to the phase isolation technologies: sedimentation (Sdm), centrifugation (Cnt) and filtration (Ftt). The performance for Sdm1 and Cnt1 is defined in terms of “efficiency” of the separation of product from the aqueous phase. For Ftt, it is defined as the retention factor of cells on the retentate side of the filter. We vary the Cnt1 efficiency and Ftt retention factor between 70 and 95% simultaneously [81, 85, 100, 101] and run the optimization model to obtain the optimal separation network and cost for each combination of these two parameters. This analysis requires the solution of about 150–200 optimization problems to generate the complete heatmap shown in Fig. 3. Note that for the efficiency of Sdm1, there is limited scope for performance enhancement [81], therefore, we fix it to 70%.
As a result, regardless of the Ftt retention factor variation, Ftt is not selected in the optimal solution (see Fig. 3a). This is because Ftt can only separate the cells from the product and water, while Cnt1 or Sdm1 functions to remove water at the same time. The optimal technologies selected are noted in the corresponding regions in Fig. 3a. The readers can identify the optimal separation network in Fig. 2 accordingly. The same notation is used in the subsequent figures. It can be seen, when the Cnt1 efficiency is below 76.5%, Sdm1 is selected in Stage I, and Cnt2 is selected in Stage II due to lower cost; otherwise, Cnt1 and Cnt2 is selected. Since the efficiency of Cnt2 is still fixed to 80%, Cnt2 is always selected as the optimal technology in Stage II.
Product titer
Sdm1 and Cnt1 equipment sizes are the major phase isolation cost drivers, and they depend on the product titer in the feed entering the separation network (affecting the total input stream flow rate). Product titer depends on the microbial strain, substrate utilization, microbial-conversion pathways, and bioreactor design. It has the potential to be altered by metabolic engineering tools [17, 102,103,104,105]. We vary the product titer from 1 to 250 g/L and obtain the costs and the corresponding optimal technology selection in Fig. 3b.
It can be seen that at high titer (≥ 32 g/L), Cnt1 alone is able to achieve the required product purity. When the titer is 14–32 g/L, another concentration technology is required in Stage II, and Mbr is preferred to Sdm2 and Cnt2 because the product loss is lower. When the titer is 4–14 g/L, Cnt2 becomes a better option than Mbr in Stage II because the low titer requires large Mbr equipment size and more costs associated with Mbr replacements. When the titer is 1.57–4 g/L, Sdm1 is a cheaper option than Cnt1 in Stage I because the major cost of Sdm is equipment cost, which is scaled with the equipment size based on the power scaling rule; however, in addition to equipment cost, Cnt also has electricity cost, which is scaled linearly with the equipment size. Therefore, when the titer is low (leading to large equipment size), Cnt1 electricity cost is high, and thus Sdm1 is a cheaper option. Finally, when the titer is < 1.57 g/L, Sdm1 in Stage I is limited by its maximum concentrating factor and thus is not able to concentrate the product enough for Stage II to reach the required purity. Therefore, Cnt1, with a higher concentrating factor, is selected.
Extension to other classes of EX NSL LT products
The other properties determining the superstructure for EX NSL LT products are the physical state (LQD/SLD) and intended use (CMD/SPC) (see Fig. 1b). The SPC (specialty) products have a high purity requirement, thus require additional treatment technologies which make their recovery more expensive in comparison to commodity products.
For a SLD product, Ftt in Stage I cannot separate the product from the cells, and thus it should be removed from the superstructure. However, since Ftt is not selected in the optimal network, our base case analysis is still able to reflect the potential technologies and the overall separation network for SLD product.
For an SPC product, separation technologies in Stage III need to be more stringent to meet the purity requirements. Thus, technologies such as Chr (chromatography) and Blc (bleaching - to remove pigments) can be included in the superstructure, and a similar analysis can be performed to identify the impact of variation in technology parameters in Stage III. However, the Stage III parameters are not selected for further analysis because, for high-value specialty chemicals, quality is a major concern and cost minimization becomes secondary [106].
Study2 – extracellular insoluble heavy (EX NSL HV) products
For this category, we choose EX NSL HV SLD CMD product as a representative base case. The key parameters used are the same with those in the EX NSL LT base case (refer Table 2). Other parameters can be found in the Additional file. Some representative products belonging to the EX NSL HV class are hyaluronic acid [54,55,56], and xanthan gum [62, 63].
Superstructure and optimal solution
By simplifying the separation scheme of an EX NSL product (see Additional file 1: Figure S.A1), we obtain the superstructure for EX NSL HV SLD CMD product, as shown in Fig. 4. Since the product is extracellular (EX), Stage I consists of phase isolation and cell removal. Phase isolation (i.e., the separation of the product-containing phase from other components in the stream), can be achieved using Slb (solubilization) or Ddg (differential digestion). Slb is used to dissolve the product in a suitable solvent to separate it from cells and other solid impurities. Ddg is used to dissolve the Non-Product Containing Materials (NPCM). Cell removal technologies include Sdm, Cnt, Ftt, and Mbr. Multiple technologies may be required in series depending on the initial product and cell concentration in the inlet stream. Flc (flocculation), as an optional pretreatment technology, can enhance the separation efficiency of subsequent tasks by enabling the formation of flocs of cells, which are then easier to isolate from the aqueous phase.
Stage II is required if the process stream undergoes solubilization (Slb) in Stage I. Stage II includes precipitation (Prc) followed by the phase separation technologies which are similar to the ones used for cell removal. If Ddg is selected in Stage I, then Stage II is not required. Stage III involves bleaching (Blc) and drying (Dry) options for final product purification.
The optimization model comprises 1618 equations, 1137continuous variables and 30 discrete variables and takes about a minute or two to solve to optimality for one set of parameter values. The solution of the superstructure optimization model yields the cost-minimal separation network, as represented by the highlighted parts in Fig. 4, as well as the corresponding cost distribution shown on the left bar alongside the superstructure. The technologies selected in Stage I are Ddg (differential digestion) for phase isolation and Cnt3 (centrifugation) and Cnt4 for cell removal. The final product refinement involves Blc (bleaching) to remove undesired color imparting impurities and Dry (drying) to retrieve product in the solid form. The overall process cost is 4.20 $/kg, where the separation cost contribution is $3.81/kg (91%). Stage I is the highest cost contributor (51%).
Analysis for EX NSL HV SLD CMD product
Performance of phase isolation and cell removal technologies
Since Ddg is used for phase isolation, Stage I is the major cost contributor in the optimal configuration, 41% of the overall cost). Its competing technology is Slb, hence, we vary the required amount of the NPCM digestion agent and the solubilizing solvent, as shown in Fig. 5a. Ddg is the preferred option even if the digestion agent is required in higher amounts as compared to Slb solvent. This is because Slb selection adds additional cost in Stage II, where Prc (precipitation) followed by phase separation is required.
The second major cost component is Blc (28% of the overall cost). However, there is limited room for performance improvements since the Blc efficiency is typically high ~ 99%. This is reasonable as Blc is a Stage III technology which already has the input stream with the product in concentrated form with traces of color-imparting impurities.
The next major cost component is Cnt3 (8.4% of the overall cost) for cell removal, and its performance affects product loss and Stage III cost. Therefore, we also vary the Cnt3 efficiency and Ftt3 retention factor between 70 and 95% simultaneously, while fixing Sdm3 efficiency to 70%, for the same reason discussed in the previous case study for EX NSL LT products. The analysis is shown in Fig. 5b. Cnt3 is the preferred option in most cases because its capital cost is lower than that of Sdm3, and Ftt3 filter replacement costs are high.
Product titer
We vary the product titer from 1 to 250 g/L. The costs and the corresponding optimal technology selection are obtained and presented in Fig. 5c. We observe that when the titer is greater than 52 g/L, Ddg-Cnt3-Dry (differential digestion-centrifugation-drying) is the optimal selection; from 15 to 52 g/L, further concentration by Cnt4 (centrifuge) is needed to achieve the final purity requirement; from 2.2 to 15 g/L, the amount of color imparting impurities is substantial compared with the amount of product present, and thus Blc (bleaching) is required; when the titer is less than 2.2 g/L, MF4 (microfiltration) replaces Cnt4 due to limitations of Cnt4 concentrating factor. Note that we have assumed that color imparting impurities exist and thus Blc is required in some cases. To account for cases where no such impurities exist (and thus Blc is not necessary), we also present the corresponding cost-titer curve in Additional file 1: Figure S.A2. Sometimes larger product titers for NSL products at high temperatures can lead to viscous process streams. However, in our technology models, we have accounted for heat dissipation effects in a centrifuge and considered appropriate cooling utility to maintain the normal operating temperatures.
Extension to other classes of EX NSL HV products
The other properties determining the superstructure for EX NSL HV products are the physical state (LQD/SLD) and intended use (CMD/SPC) (see Fig. 1b). For a LQD product, instead of using Ddg or Slb (see Fig. 4), Ftt or Mbr (depending on the size of solid impurities) can be used directly to separate the product from cells, followed by concentration using Sdm, Cnt or Mbr. Also, Dry in Stage III is only applicable to SLD products and thus should be removed. Therefore, for LQD product, the separation cost will be lower. We modify the base case superstructure to account for LQD product, and the optimal technology selections are Flc-Ftt-MF-Blc (flocculation-filtration-microfiltration-bleaching), and the cost is decreased from 4.2 $/kg (for SLD) to 3.52 $/kg (for LQD), as shown in Additional file 1: Figure S.A3. For a SPC product, separation technologies in Stage III need to be more stringent for high purity requirements. Similar considerations mentioned for EX NSL LT products will be valid here.
Study3 – extracellular soluble (EX SOL) products
For this category, we choose EX SOL LQD NVL CMD product as the base case. The key parameters used are the same as those in the EX NSL LT base case (refer Table 2). The remaining parameters can be found in the Additional file 1: Table S.C.B1. The typical examples for EX SOL products from microbial sources include propanediol [107, 108] and butanediol [64,65,66,67].
Superstructure and optimal solution
By simplifying the separation scheme of an EX SOL product (see Additional file 1: Figure S.B1), we obtain the superstructure for EX SOL LQD NVL CMD product, as shown in Fig. 6. In Stage I, cells are first removed, and thus the aqueous phase containing the product is isolated. Then, in Stage II, Dst (distillation), Atpe (aqueous two-phase extraction), and Ext (extraction) are considered as concentrating technologies. Finally, in Stage III, Pvp (pervaporation) can be used to remove small amount of remaining impurities, if necessary. The optimization model comprises 969 equations, 754 continuous variables and 25 discrete variables and takes about 30 min to 5 h to solve the problem to optimality for one set of parameter values. This is due to the presence of the distillation and extraction model equations which include logarithmic and exponential terms involving mole fractions of components.
After solving the superstructure model, we obtain the cost-minimal separation network represented by the highlighted parts in Fig. 6, as well as the corresponding cost distribution. The technologies selected in stage I are Flc (flocculation), Cnt1 (centrifuge), and Cnt2. Atpe, followed by MF3 (microfiltration) and MF4, is selected for concentration in Stage II. Stage III is bypassed. The total minimum cost is 3.98 $/kg, with the feed accounting for 7%, Stage I 16% and Stage II 77%. Stage II is the major cost component because the separation for soluble product requires concentration of the product present in water-rich phase.
Alternative concentration options in stage II
Atpe (aqueous two-phase extraction) is selected in the optimal network for the base case. However, we also analyze the separation networks when Ext (extraction) or Dst1 (distillation) is selected in Stage II, thus accounting for cases where effective Atpe using alcohol-salt [109], polymer-salt [110] or polymer-polymer [111, 112] phases for the product may not exist.
When Ext (as well as its auxiliary Dst 2) is selected (by setting the binary variable for Ext to 1 to ensure selection), the technologies selected (see Fig. 7a) in Stage I include Flc, Cnt1, and Cnt2. Stage III is bypassed. The overall process cost is 4.06 $/kg, and the separation cost contribution is 3.77 $/kg (94%). Stage II is still the major cost contributor (80%).
When Dst1 is selected, the technologies selected (see Fig. 7b) in Stage I include Flc, Cnt1, and Cnt2. Stage III is bypassed. The overall process cost is 8.19 $/kg, and the separation cost contribution is 7.93 $/kg (97%). Stage II is still the major cost contributor (90%). Direct Dst is costly because a large amount of water needs to be vaporized in Dst1, which increases the heating utility costs in the system.
We can see a major difference in the cost distribution for soluble products based on their localization. From the previous study on intracellular products [51], the major cost driver was Stage I irrespective of the other properties of the product, however, for extracellular products, Stage II costs take precedence and have a major impact on the overall downstream separation network selection.
Analysis for EX SOL LQD NVL CMD product
Since Stage II is the major cost driver, we perform the analysis with varying parameters related to Atpe, Ext and Dst1.
Partition coefficient for Atpe
The major cost driver of Atpe is the equipment size, which is a function of the flowrates of feed and added separating agents (alcohol/polymer for the top phase and salt for the bottom phase). The amount of agents required is affected by the partition coefficient for the top phase (KpT), which we vary from 3 to 10 while fixing Atpe selection. These are valid choices for the KpT values based on the literature for typical Atpe systems [66, 109]. In Fig. 8a, we observe that the overall cost decreases from 4.25 $/kg to 3.63 $/kg with the increase of partition coefficient from 3 to 10, which enables more product to be extracted into the top phase and thus reduces the amount of separating agents required. When the Atpe partition coefficient is below 3.7 (thus cost> 4.06 $/kg), Atpe becomes a more expensive option than Ext, assuming the base case parameters for Ext (see Fig. 7a).
Ext parameters –solvent solubility, cost and partition coefficient
For Ext, the solvent is the major cost contributor, and thus we vary solubility of solvent in water from 0.0002 to 0.03 kg/kg and vary the cost of solvent from 0.2 to 1.5 $/kg while fixing Ext selection, and the analysis result is shown in Fig. 8b. Compared with the base case, if the solubility decreases from 0.02 to 0.0002 kg/kg and the solvent cost decreases from 1.2 to 0.2 $/kg, then the overall cost will be reduced from 4.06 to 1.06 $/kg, a 74% reduction. Thus, if both parameters can be improved such that the overall cost is lower than 3.98 $/kg, then Ext becomes a cheaper option than Atpe, assuming the base case parameters for Atpe (see Fig. 6). Also, note that the change in partition coefficient usually does not have a significant impact on the Ext cost when compared with the solvent cost and solubility (see Additional file 1: Figure S.B2). Therefore, targeting solvents with low water solubility and cost, even if the partition coefficients are low, can help reduce cost because less solvent will be lost.
Dst parameter- relative volatility
For Dst1, we vary the volatility of water relative to the product from 1.1 to 3, and the result is shown in Fig. 8c, where the cost ranges from 4.34 to 20.35 $/kg. Thus, even if relative volatility of 3 can be achieved, direct Dst is still more expensive than Atpe or Ext, assuming their base case parameters.
Product titer
We vary the product titer from 1 to 250 g/L. As a result, the costs and the corresponding optimal technology selection are obtained and presented in Fig. 9. We observe that when the titer is greater than 7.5 g/L, Ext selection is optimal; otherwise, Atpe selection is optimal.
Extension to other classes of EX SOL products
The other properties determining the superstructure for EX SOL products are physical state (LQD/SLD), volatility (VOL/NVL), and intended use (CMD/SPC) (see Fig. 1b).
For a SLD product (such as a soluble salt), Mbr and Prc can be used for product concentration as an alternative to Dst, Atpe, and Ext in Stage II. Also, in Stage III, Ads (adsorption), Crs (crystallization) and Dry (drying) can be considered. However, the cost will not likely be influenced because the desired product purity is already achieved without Stage III in the base case.
For a volatile (VOL) product, the product will be obtained at the top instead of at the bottom in Dst. Also, direct Dst is typically cheaper than Ext or Atpe when the relative volatility is greater than 1.05 [51, 113]. For a specialty (SPC) product, the same argument about additional purification and stringent requirements in Stage III technologies will be valid due to high purity requirements.