In Silico Elucidation of Key Drivers of Staphyloccocus aureus-Staphyloccocus epidermidis-Induced Skin Damage in Atopic Dermatitis Lesions.

Staphylococcus aureus (SA) colonizes and can damage skin in atopic dermatitis lesions, despite being commonly found with Staphylococcus epidermidis (SE), a commensal that can inhibit SA's virulence and kill SA. In this study, we developed an in silico model, termed a virtual skin site, describing the dynamic interplay between SA, SE, and the skin barrier in atopic dermatitis lesions to investigate the mechanisms driving skin damage by SA and SE. We generated 106 virtual skin sites by varying model parameters to represent different skin physiologies and bacterial properties. In silico analysis revealed that virtual skin sites with no skin damage in the model were characterized by parameters representing stronger SA and SE growth attenuation than those with skin damage. This inspired an in silico treatment strategy combining SA-killing with an enhanced SA-SE growth attenuation, which was found through simulations to recover many more damaged virtual skin sites to a non-damaged state, compared with SA-killing alone. This study demonstrates that in silico modelling can help elucidate the key factors driving skin damage caused by SA-SE colonization in atopic dermatitis lesions and help propose strategies to control it, which we envision will contribute to the design of promising treatments for clinical studies.

gcFront: a tool for determining a Pareto front of growth-coupled cell factory designs.

MOTIVATION: A widely applicable strategy to create cell factories is to knockout (KO) genes or reactions to redirect cell metabolism so that chemical synthesis is made obligatory when the cell grows at its maximum rate. Synthesis is thus growth-coupled, and the stronger the coupling the more deleterious any impediments in synthesis are to cell growth, making high producer phenotypes evolutionarily robust. Additionally, we desire that these strains grow and synthesize at high rates. Genome-scale metabolic models can be used to explore and identify KOs that growth-couple synthesis, but these are rare in an immense design space, making the search difficult and slow. RESULTS: To address this multi-objective optimization problem, we developed a software tool named gcFront-using a genetic algorithm it explores KOs that maximize cell growth, product synthesis and coupling strength. Moreover, our measure of coupling strength facilitates the search so that gcFront not only finds a growth-coupled design in minutes but also outputs many alternative Pareto optimal designs from a single run-granting users flexibility in selecting designs to take to the lab. AVAILABILITY AND IMPLEMENTATION: gcFront, with documentation and a workable tutorial, is freely available at GitHub: https://github.com/lLegon/gcFront and archived at Zenodo, DOI: 10.5281/zenodo.5557755. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Trade-Offs in Biosensor Optimization for Dynamic Pathway Engineering.

Recent progress in synthetic biology allows the construction of dynamic control circuits for metabolic engineering. This technology promises to overcome many challenges encountered in traditional pathway engineering, thanks to its ability to self-regulate gene expression in response to bioreactor perturbations. The central components in these control circuits are metabolite biosensors that read out pathway signals and actuate enzyme expression. However, the construction of metabolite biosensors is a major bottleneck for strain design, and a key challenge is to understand the relation between biosensor dose-response curves and pathway performance. Here we employ multiobjective optimization to quantify performance trade-offs that arise in the design of metabolite biosensors. Our approach reveals strategies for tuning dose-response curves along an optimal trade-off between production flux and the cost of an increased expression burden on the host. We explore properties of control architectures built in the literature and identify their advantages and caveats in terms of performance and robustness to growth conditions and leaky promoters. We demonstrate the optimality of a control circuit for glucaric acid production in Escherichia coli, which has been shown to increase the titer by 2.5-fold as compared to static designs. Our results lay the groundwork for the automated design of control circuits for pathway engineering, with applications in the food, energy, and pharmaceutical sectors.

Designing an irreversible metabolic switch for scalable induction of microbial chemical production.

Bacteria can be harnessed to synthesise high-value chemicals. A promising strategy for increasing productivity uses inducible control systems to switch metabolism from growth to chemical synthesis once a large population of cell factories are generated. However, use of expensive chemical inducers limits scalability of this approach for biotechnological applications. Switching using cheap nutrients is an appealing alternative, but their tightly regulated uptake and consumption again limits scalability. Here, using mathematical models of fatty acid uptake in E. coli as an exemplary case study, we unravel how the cell's native regulation and program of induction can be engineered to minimise inducer usage. We show that integrating positive feedback loops into the circuitry creates an irreversible metabolic switch, which, requiring only temporary induction, drastically reduces inducer usage. Our proposed switch should be widely applicable, irrespective of the product of interest, and brings closer the realization of scalable and sustainable microbial chemical production.

Metabolite Sequestration Enables Rapid Recovery from Fatty Acid Depletion in Escherichia coli.

Microbes adapt their metabolism to take advantage of nutrients in their environment. Such adaptations control specific metabolic pathways to match energetic demands with nutrient availability. Upon depletion of nutrients, rapid pathway recovery is key to release cellular resources required for survival under the new nutritional conditions. Yet, little is known about the regulatory strategies that microbes employ to accelerate pathway recovery in response to nutrient depletion. Using the fatty acid catabolic pathway in Escherichia coli, here, we show that fast recovery can be achieved by rapid release of a transcriptional regulator from a metabolite-sequestered complex. With a combination of mathematical modeling and experiments, we show that recovery dynamics depend critically on the rate of metabolite consumption and the exposure time to nutrients. We constructed strains with rewired transcriptional regulatory architectures that highlight the metabolic benefits of negative autoregulation over constitutive and positive autoregulation. Our results have wide-ranging implications for our understanding of metabolic adaptations, as well as for guiding the design of gene circuitry for synthetic biology and metabolic engineering.IMPORTANCE Rapid metabolic recovery during nutrient shift is critical to microbial survival, cell fitness, and competition among microbiota, yet little is known about the regulatory mechanisms of rapid metabolic recovery. This work demonstrates a previously unknown mechanism where rapid release of a transcriptional regulator from a metabolite-sequestered complex enables fast recovery to nutrient depletion. The work identified key regulatory architectures and parameters that control the speed of recovery, with wide-ranging implications for the understanding of metabolic adaptations as well as synthetic biology and metabolic engineering.

Dynamic metabolic control: towards precision engineering of metabolism.

Advances in metabolic engineering have led to the synthesis of a wide variety of valuable chemicals in microorganisms. The key to commercializing these processes is the improvement of titer, productivity, yield, and robustness. Traditional approaches to enhancing production use the "push-pull-block" strategy that modulates enzyme expression under static control. However, strains are often optimized for specific laboratory set-up and are sensitive to environmental fluctuations. Exposure to sub-optimal growth conditions during large-scale fermentation often reduces their production capacity. Moreover, static control of engineered pathways may imbalance cofactors or cause the accumulation of toxic intermediates, which imposes burden on the host and results in decreased production. To overcome these problems, the last decade has witnessed the emergence of a new technology that uses synthetic regulation to control heterologous pathways dynamically, in ways akin to regulatory networks found in nature. Here, we review natural metabolic control strategies and recent developments in how they inspire the engineering of dynamically regulated pathways. We further discuss the challenges of designing and engineering dynamic control and highlight how model-based design can provide a powerful formalism to engineer dynamic control circuits, which together with the tools of synthetic biology, can work to enhance microbial production.

Fundamental Design Principles for Transcription-Factor-Based Metabolite Biosensors.

Metabolite biosensors are central to current efforts toward precision engineering of metabolism. Although most research has focused on building new biosensors, their tunability remains poorly understood and is fundamental for their broad applicability. Here we asked how genetic modifications shape the dose-response curve of biosensors based on metabolite-responsive transcription factors. Using the lac system in Escherichia coli as a model system, we built promoter libraries with variable operator sites that reveal interdependencies between biosensor dynamic range and response threshold. We developed a phenomenological theory to quantify such design constraints in biosensors with various architectures and tunable parameters. Our theory reveals a maximal achievable dynamic range and exposes tunable parameters for orthogonal control of dynamic range and response threshold. Our work sheds light on fundamental limits of synthetic biology designs and provides quantitative guidelines for biosensor design in applications such as dynamic pathway control, strain optimization, and real-time monitoring of metabolism.

Signaling Tug-of-War Delivers the Whole Message.

How do cells transmit biochemical signals accurately? It turns out, pushing and pulling can go a long way.

MUFINS: multi-formalism interaction network simulator.

Systems Biology has established numerous approaches for mechanistic modeling of molecular networks in the cell and a legacy of models. The current frontier is the integration of models expressed in different formalisms to address the multi-scale biological system organization challenge. We present MUFINS (MUlti-Formalism Interaction Network Simulator) software, implementing a unique set of approaches for multi-formalism simulation of interaction networks. We extend the constraint-based modeling (CBM) framework by incorporation of linear inhibition constraints, enabling for the first time linear modeling of networks simultaneously describing gene regulation, signaling and whole-cell metabolism at steady state. We present a use case where a logical hypergraph model of a regulatory network is expressed by linear constraints and integrated with a Genome-Scale Metabolic Network (GSMN) of mouse macrophage. We experimentally validate predictions, demonstrating application of our software in an iterative cycle of hypothesis generation, validation and model refinement. MUFINS incorporates an extended version of our Quasi-Steady State Petri Net approach to integrate dynamic models with CBM, which we demonstrate through a dynamic model of cortisol signaling integrated with the human Recon2 GSMN and a model of nutrient dynamics in physiological compartments. Finally, we implement a number of methods for deriving metabolic states from ~omics data, including our new variant of the iMAT congruency approach. We compare our approach with iMAT through the analysis of 262 individual tumor transcriptomes, recovering features of metabolic reprogramming in cancer. The software provides graphics user interface with network visualization, which facilitates use by researchers who are not experienced in coding and mathematical modeling environments.

Integrating Kinetic Model of E. coli with Genome Scale Metabolic Fluxes Overcomes Its Open System Problem and Reveals Bistability in Central Metabolism.

An understanding of the dynamics of the metabolic profile of a bacterial cell is sought from a dynamical systems analysis of kinetic models. This modelling formalism relies on a deterministic mathematical description of enzyme kinetics and their metabolite regulation. However, it is severely impeded by the lack of available kinetic information, limiting the size of the system that can be modelled. Furthermore, the subsystem of the metabolic network whose dynamics can be modelled is faced with three problems: how to parameterize the model with mostly incomplete steady state data, how to close what is now an inherently open system, and how to account for the impact on growth. In this study we address these challenges of kinetic modelling by capitalizing on multi-'omics' steady state data and a genome-scale metabolic network model. We use these to generate parameters that integrate knowledge embedded in the genome-scale metabolic network model, into the most comprehensive kinetic model of the central carbon metabolism of E. coli realized to date. As an application, we performed a dynamical systems analysis of the resulting enriched model. This revealed bistability of the central carbon metabolism and thus its potential to express two distinct metabolic states. Furthermore, since our model-informing technique ensures both stable states are constrained by the same thermodynamically feasible steady state growth rate, the ensuing bistability represents a temporal coexistence of the two states, and by extension, reveals the emergence of a phenotypically heterogeneous population.

Interrogation of global mutagenesis data with a genome scale model of Neisseria meningitidis to assess gene fitness in vitro and in sera.

BACKGROUND: Neisseria meningitidis is an important human commensal and pathogen that causes several thousand deaths each year, mostly in young children. How the pathogen replicates and causes disease in the host is largely unknown, particularly the role of metabolism in colonization and disease. Completed genome sequences are available for several strains but our understanding of how these data relate to phenotype remains limited. RESULTS: To investigate the metabolism of N. meningitidis we generated and then selected a representative Tn5 library on rich medium, a minimal defined medium and in human serum to identify genes essential for growth under these conditions. To relate these data to a systems-wide understanding of the pathogen's biology we constructed a genome-scale metabolic network: Nmb_iTM560. This model was able to distinguish essential and non-essential genes as predicted by the global mutagenesis. These essentiality data, the library and the Nmb_iTM560 model are powerful and widely applicable resources for the study of meningococcal metabolism and physiology. We demonstrate the utility of these resources by predicting and demonstrating metabolic requirements on minimal medium, such as a requirement for phosphoenolpyruvate carboxylase, and by describing the nutritional and biochemical status of N. meningitidis when grown in serum, including a requirement for both the synthesis and transport of amino acids. CONCLUSIONS: This study describes the application of a genome scale transposon library combined with an experimentally validated genome-scale metabolic network of N. meningitidis to identify essential genes and provide novel insight into the pathogen's metabolism both in vitro and during infection.

Modeling and simulation of the main metabolism in Escherichia coli and its several single-gene knockout mutants with experimental verification.

BACKGROUND: It is quite important to simulate the metabolic changes of a cell in response to the change in culture environment and/or specific gene knockouts particularly for the purpose of application in industry. If this could be done, the cell design can be made without conducting exhaustive experiments, and one can screen out the promising candidates, proceeded by experimental verification of a select few of particular interest. Although several models have so far been proposed, most of them focus on the specific metabolic pathways. It is preferred to model the whole of the main metabolic pathways in Escherichia coli, allowing for the estimation of energy generation and cell synthesis, based on intracellular fluxes and that may be used to characterize phenotypic growth. RESULTS: In the present study, we considered the simulation of the main metabolic pathways such as glycolysis, TCA cycle, pentose phosphate (PP) pathway, and the anapleorotic pathways using enzymatic reaction models of E. coli. Once intracellular fluxes were computed by this model, the specific ATP production rate, the specific CO2 production rate, and the specific NADPH production rate could be estimated. The specific ATP production rate thus computed was used for the estimation of the specific growth rate. The CO2 production rate could be used to estimate cell yield, and the specific NADPH production rate could be used to determine the flux of the oxidative PP pathway. The batch and continuous cultivations were simulated where the changing patterns of extracellular and intra-cellular metabolite concentrations were compared with experimental data. Moreover, the effects of the knockout of such pathways as Ppc, Pck and Pyk on the metabolism were simulated. It was shown to be difficult for the cell to grow in Ppc mutant due to low concentration of OAA, while Pck mutant does not necessarily show this phenomenon. The slower growth rate of the Ppc mutant was properly estimated by taking into account the lower specific ATP production rate. In the case of Pyk mutant, the enzyme level regulation was made clear such that Pyk knockout caused PEP concentration to be up-regulated and activated Ppc, which caused the increase in MAL concentration and backed up reduced PYR through Mez, resulting in the phenotypic growth characteristics similar to the wild type. CONCLUSIONS: It was shown to be useful to simulate the main metabolism of E. coli for understanding metabolic changes inside the cell in response to specific pathway gene knockouts, considering the whole main metabolic pathways. The comparison of the simulation result with the experimental data indicates that the present model could simulate the effect of the specific gene knockouts to the changes in the metabolisms to some extent.