ABSTRACT
There are several ways in which mathematical modeling is used in fermentation control, but mechanistic mathematical genome-scale models of metabolism within the cell have not been applied or implemented so far. As part of the metabolic engineering task setting, we propose that metabolite fluxes and/or biomass growth rate be used to search for a fermentation steady state marker rule. During fermentation, the bioreactor control system can automatically detect the desired steady state using a logical marker rule. The marker rule identification can be also integrated with the production growth coupling approach, as presented in this study. A design of strain with marker rule is demonstrated on genome scale metabolic model iML1515 of Escherichia coli MG1655 proposing two gene deletions enabling a measurable marker rule for succinate production using glucose as a substrate. The marker rule example at glucose consumption 10.0 is: IF (specific growth rate µ is above 0.060 h-1, AND CO2 production under 1.0, AND ethanol production above 5.5), THEN succinate production is within the range 8.2-10, where all metabolic fluxes units are mmol ∗ gDW-1 ∗ h-1. An objective function for application in metabolic engineering, including productivity features and rule detecting sensor set characterizing parameters, is proposed. Two-phase approach to implementing marker rules in the cultivation control system is presented to avoid the need for a modeler during production.
ABSTRACT
Finding the best knockout strategy for coupling biomass growth and production of a target metabolite using a mathematic model of metabolism is a challenge in biotechnology. In this research, a three-step method named OptEnvelope is presented based on finding minimal set of active reactions for a target point in the feasible solution space (envelope) using a mixed-integer linear programming formula. The method initially finds the reduced desirable solution space envelope in the product versus biomass plot by removing all inactive reactions. Then, with reinsertion of the deleted reactions, OptEnvelope attempts to reduce the number of knockouts so that the desirable production envelope is preserved. Additionally, OptEnvelope searches for envelopes with higher minimum production rates or fewer knockouts by evaluating different target points within the desired solution space. It is possible to limit the maximal number of knockouts. The method was implemented on metabolic models of E. coli and S. cerevisiae to test the method benchmarking the capability of these industrial microbes for overproduction of acetate and glycerol under aerobic conditions and succinate and ethanol under anaerobic conditions. The results illustrate that OptEnvelope is capable to find multiple strong coupled envelopes located in the desired solution space because of its novel target point oriented strategy of envelope search. The results indicate that E. coli is more appropriate to produce acetate and succinate while S. cerevisiae is a better host for glycerol production. Gene deletions for some of the proposed reaction knockouts have been previously reported to increase the production of these metabolites in experiments. Both organisms are suitable for ethanol production, however, more knockouts for the adaptation of E. coli are required. OptEnvelope is available at https://github.com/lv-csbg/optEnvelope.
Subject(s)
Escherichia coli , Saccharomyces cerevisiae , Escherichia coli/metabolism , Saccharomyces cerevisiae/genetics , Saccharomyces cerevisiae/metabolism , Glycerol/metabolism , Acetates/metabolism , Succinates/metabolism , Ethanol/metabolism , Metabolic Networks and PathwaysABSTRACT
Docosahexaenoic acid (DHA) is one of the most important long-chain polyunsaturated fatty acids (LC-PUFAs), with numerous health benefits. Crypthecodinium cohnii, a marine heterotrophic dinoflagellate, is successfully used for the industrial production of DHA because it can accumulate DHA at high concentrations within the cells. Glycerol is an interesting renewable substrate for DHA production since it is a by-product of biodiesel production and other industries, and is globally generated in large quantities. The DHA production potential from glycerol, ethanol and glucose is compared by combining fermentation experiments with the pathway-scale kinetic modeling and constraint-based stoichiometric modeling of C. cohnii metabolism. Glycerol has the slowest biomass growth rate among the tested substrates. This is partially compensated by the highest PUFAs fraction, where DHA is dominant. Mathematical modeling reveals that glycerol has the best experimentally observed carbon transformation rate into biomass, reaching the closest values to the theoretical upper limit. In addition to our observations, the published experimental evidence indicates that crude glycerol is readily consumed by C. cohnii, making glycerol an attractive substrate for DHA production.
Subject(s)
Dinoflagellida/metabolism , Docosahexaenoic Acids/metabolism , Models, Theoretical , Biomass , Ethanol/metabolism , Fermentation , Glucose/metabolism , Glycerol/metabolismABSTRACT
BACKGROUND: Metformin has been used for the treatment of type 2 diabetes for over 60 years; however, its mechanism of pharmacological action is not fully clear. Different hypotheses exist regarding metformin distribution and redistribution mechanisms between plasma and erythrocytes/red blood cells (RBCs). OBJECTIVE: We aimed to test the hypothesis that the metformin distribution between plasma and RBC occurs via concentration difference-driven passive transport and estimated transport rate coefficient values based on metformin concentration time series in plasma and RBCs from in vivo studies. METHODS: An ordinary differential equation (ODE) system with two compartments was used to describe diffusion-based passive transport between plasma and RBCs. Metformin concentration time series in plasma and RBCs of 35 individuals were used for metformin transport parametrization. Plasma concentration has been approximated by biexponential decline. RESULTS: A single passive transport coefficient, k = 0.044 ± 0.014 (h-1), can be applied, describing the uptake and release transport rate versus the linear equation v = k × (Mpl - MRBC), where Mpl is the metformin concentration in plasma and MRBC is the metformin concentration in RBCs. CONCLUSIONS: Our research suggests that passive transport can explain metformin distribution dynamics between plasma and RBCs because transport speed is proportional to the metformin concentration difference and independent of the transport direction. Concentration difference-driven passive transport can explain the mechanism of faster metformin distribution to RBCs the first few hours after administration, and faster release and domination of the redistribution transport rate after metformin concentration in plasma becomes smaller than in RBCs.