Investigating the metabolic capabilities of Mycobacterium tuberculosis H37Rv using the in silico strain iNJ661 and proposing alternative drug targets

<p>Abstract</p> <p>Background:</p> <p><it>Mycobacterium tuberculosis </it>continues to be a major pathogen in the third world, killing almost 2 million people a year by the most recent estimates. Even in industrialized countries, the emergence of multi-drug resistant (MDR) strains of tuberculosis hails the need to develop additional medications for treatment. Many of the drugs used for treatment of tuberculosis target metabolic enzymes. Genome-scale models can be used for analysis, discovery, and as hypothesis generating tools, which will hopefully assist the rational drug development process. These models need to be able to assimilate data from large datasets and analyze them.</p> <p>Results:</p> <p>We completed a bottom up reconstruction of the metabolic network of <it>Mycobacterium tuberculosis </it>H37Rv. This functional <it>in silico </it>bacterium, <it>iNJ</it>661, contains 661 genes and 939 reactions and can produce many of the complex compounds characteristic to tuberculosis, such as mycolic acids and mycocerosates. We grew this bacterium <it>in silico </it>on various media, analyzed the model in the context of multiple high-throughput data sets, and finally we analyzed the network in an 'unbiased' manner by calculating the Hard Coupled Reaction (HCR) sets, groups of reactions that are forced to operate in unison due to mass conservation and connectivity constraints.</p> <p>Conclusion:</p> <p>Although we observed growth rates comparable to experimental observations (doubling times ranging from about 12 to 24 hours) in different media, comparisons of gene essentiality with experimental data were less encouraging (generally about 55%). The reasons for the often conflicting results were multi-fold, including gene expression variability under different conditions and lack of complete biological knowledge. Some of the inconsistencies between <it>in vitro </it>and <it>in silico </it>or <it>in vivo </it>and <it>in silico </it>results highlight specific loci that are worth further experimental investigations. Finally, by considering the HCR sets in the context of known drug targets for tuberculosis treatment we proposed new alternative, but equivalent drug targets.</p

Spatiotemporal modeling of microbial metabolism

Background Microbial systems in which the extracellular environment varies both spatially and temporally are very common in nature and in engineering applications. While the use of genome-scale metabolic reconstructions for steady-state flux balance analysis (FBA) and extensions for dynamic FBA are common, the development of spatiotemporal metabolic models has received little attention. Results We present a general methodology for spatiotemporal metabolic modeling based on combining genome-scale reconstructions with fundamental transport equations that govern the relevant convective and/or diffusional processes in time and spatially varying environments. Our solution procedure involves spatial discretization of the partial differential equation model followed by numerical integration of the resulting system of ordinary differential equations with embedded linear programs using DFBAlab, a MATLAB code that performs reliable and efficient dynamic FBA simulations. We demonstrate our methodology by solving spatiotemporal metabolic models for two systems of considerable practical interest: (1) a bubble column reactor with the syngas fermenting bacterium Clostridium ljungdahlii; and (2) a chronic wound biofilm with the human pathogen Pseudomonas aeruginosa. Despite the complexity of the discretized models which consist of 900 ODEs/600 LPs and 250 ODEs/250 LPs, respectively, we show that the proposed computational framework allows efficient and robust model solution. Conclusions Our study establishes a new paradigm for formulating and solving genome-scale metabolic models with both time and spatial variations and has wide applicability to natural and engineered microbial systems

Integrating Flux Balance Analysis into Kinetic Models to Decipher the Dynamic Metabolism of Shewanella oneidensis MR-1

Shewanella oneidensis MR-1 sequentially utilizes lactate and its waste products (pyruvate and acetate) during batch culture. To decipher MR-1 metabolism, we integrated genome-scale flux balance analysis (FBA) into a multiple-substrate Monod model to perform the dynamic flux balance analysis (dFBA). The dFBA employed a static optimization approach (SOA) by dividing the batch time into small intervals (i.e., ∼400 mini-FBAs), then the Monod model provided time-dependent inflow/outflow fluxes to constrain the mini-FBAs to profile the pseudo-steady-state fluxes in each time interval. The mini-FBAs used a dual-objective function (a weighted combination of “maximizing growth rate” and “minimizing overall flux”) to capture trade-offs between optimal growth and minimal enzyme usage. By fitting the experimental data, a bi-level optimization of dFBA revealed that the optimal weight in the dual-objective function was time-dependent: the objective function was constant in the early growth stage, while the functional weight of minimal enzyme usage increased significantly when lactate became scarce. The dFBA profiled biologically meaningful dynamic MR-1 metabolisms: 1. the oxidative TCA cycle fluxes increased initially and then decreased in the late growth stage; 2. fluxes in the pentose phosphate pathway and gluconeogenesis were stable in the exponential growth period; and 3. the glyoxylate shunt was up-regulated when acetate became the main carbon source for MR-1 growth

Optimization of fed-batch Saccharomyces cerevisiae fermentation using dynamic flux balance models

We developed a dynamic flux balance model for fed-batch Saccharomyces cereVisiae fermentation that couples a detailed steady-state description of primary carbon metabolism with dynamic mass balances on key extracellular species. Model-based dynamic optimization is performed to determine fed-batch operating policies that maximize ethanol productivity and/or ethanol yield on glucose. The initial volume and glucose concentrations, the feed flow rate and dissolved oxygen concentration profiles, and the final batch time are treated as decision variables in the dynamic optimization problem. Optimal solutions are generated to analyze the tradeoff between maximal productivity and yield objectives. We find that for both cases the prediction of a microaerobic region is significant. The optimization results are sensitive to network model parameters for the growth associated maintenance and P/O ratio. The results of our computational study motivate continued development of dynamic flux balance models and further exploration of their application to productivity optimization in biochemical reactors