Rational strain design with minimal phenotype perturbation.

Devising genetic interventions for desired cellular phenotypes remains challenging regarding time and resources. Kinetic models can accelerate this task by simulating metabolic responses to genetic perturbations. However, exhaustive design evaluations with kinetic models are computationally impractical, especially when targeting multiple enzymes. Here, we introduce a framework for efficiently scouting the design space while respecting cellular physiological requirements. The framework employs mixed-integer linear programming and nonlinear simulations with large-scale nonlinear kinetic models to devise genetic interventions while accounting for the network effects of these perturbations. Importantly, it ensures the engineered strain's robustness by maintaining its phenotype close to that of the reference strain. The framework, applied to improve the anthranilate production in E. coli, devises designs for experimental implementation, including eight previously experimentally validated targets. We expect this framework to play a crucial role in future design-build-test-learn cycles, significantly expediting the strain design compared to exhaustive design enumeration.

From microbiome composition to functional engineering, one step at a time.

SUMMARYCommunities of microorganisms (microbiota) are present in all habitats on Earth and are relevant for agriculture, health, and climate. Deciphering the mechanisms that determine microbiota dynamics and functioning within the context of their respective environments or hosts (the microbiomes) is crucially important. However, the sheer taxonomic, metabolic, functional, and spatial complexity of most microbiomes poses substantial challenges to advancing our knowledge of these mechanisms. While nucleic acid sequencing technologies can chart microbiota composition with high precision, we mostly lack information about the functional roles and interactions of each strain present in a given microbiome. This limits our ability to predict microbiome function in natural habitats and, in the case of dysfunction or dysbiosis, to redirect microbiomes onto stable paths. Here, we will discuss a systematic approach (dubbed the N+1/N-1 concept) to enable step-by-step dissection of microbiome assembly and functioning, as well as intervention procedures to introduce or eliminate one particular microbial strain at a time. The N+1/N-1 concept is informed by natural invasion events and selects culturable, genetically accessible microbes with well-annotated genomes to chart their proliferation or decline within defined synthetic and/or complex natural microbiota. This approach enables harnessing classical microbiological and diversity approaches, as well as omics tools and mathematical modeling to decipher the mechanisms underlying N+1/N-1 microbiota outcomes. Application of this concept further provides stepping stones and benchmarks for microbiome structure and function analyses and more complex microbiome intervention strategies.

Metabolic interaction models recapitulate leaf microbiota ecology.

Resource allocation affects the structure of microbiomes, including those associated with living hosts. Understanding the degree to which this dependency determines interspecies interactions may advance efforts to control host-microbiome relationships. We combined synthetic community experiments with computational models to predict interaction outcomes between plant-associated bacteria. We mapped the metabolic capabilities of 224 leaf isolates from Arabidopsis thaliana by assessing the growth of each strain on 45 environmentally relevant carbon sources in vitro. We used these data to build curated genome-scale metabolic models for all strains, which we combined to simulate >17,500 interactions. The models recapitulated outcomes observed in planta with >89% accuracy, highlighting the role of carbon utilization and the contributions of niche partitioning and cross-feeding in the assembly of leaf microbiomes.

Optimal enzyme utilization suggests that concentrations and thermodynamics determine binding mechanisms and enzyme saturations.

Deciphering the metabolic functions of organisms requires understanding the dynamic responses of living cells upon genetic and environmental perturbations, which in turn can be inferred from enzymatic activity. In this work, we investigate the optimal modes of operation for enzymes in terms of the evolutionary pressure driving them toward increased catalytic efficiency. We develop a framework using a mixed-integer formulation to assess the distribution of thermodynamic forces and enzyme states, providing detailed insights into the enzymatic mode of operation. We use this framework to explore Michaelis-Menten and random-ordered multi-substrate mechanisms. We show that optimal enzyme utilization is achieved by unique or alternative operating modes dependent on reactant concentrations. We find that in a bimolecular enzyme reaction, the random mechanism is optimal over any other ordered mechanism under physiological conditions. Our framework can investigate the optimal catalytic properties of complex enzyme mechanisms. It can further guide the directed evolution of enzymes and fill in the knowledge gaps in enzyme kinetics.

Dynamics of CLIMP-63 S-acylation control ER morphology.

The complex architecture of the endoplasmic reticulum (ER) comprises distinct dynamic features, many at the nanoscale, that enable the coexistence of the nuclear envelope, regions of dense sheets and a branched tubular network that spans the cytoplasm. A key player in the formation of ER sheets is cytoskeleton-linking membrane protein 63 (CLIMP-63). The mechanisms by which CLIMP-63 coordinates ER structure remain elusive. Here, we address the impact of S-acylation, a reversible post-translational lipid modification, on CLIMP-63 cellular distribution and function. Combining native mass-spectrometry, with kinetic analysis of acylation and deacylation, and data-driven mathematical modelling, we obtain in-depth understanding of the CLIMP-63 life cycle. In the ER, it assembles into trimeric units. These occasionally exit the ER to reach the plasma membrane. However, the majority undergoes S-acylation by ZDHHC6 in the ER where they further assemble into highly stable super-complexes. Using super-resolution microscopy and focused ion beam electron microscopy, we show that CLIMP-63 acylation-deacylation controls the abundance and fenestration of ER sheets. Overall, this study uncovers a dynamic lipid post-translational regulation of ER architecture.

Symbolic kinetic models in python (SKiMpy): intuitive modeling of large-scale biological kinetic models.

MOTIVATION: Large-scale kinetic models are an invaluable tool to understand the dynamic and adaptive responses of biological systems. The development and application of these models have been limited by the availability of computational tools to build and analyze large-scale models efficiently. The toolbox presented here provides the means to implement, parameterize and analyze large-scale kinetic models intuitively and efficiently. RESULTS: We present a Python package (SKiMpy) bridging this gap by implementing an efficient kinetic modeling toolbox for the semiautomatic generation and analysis of large-scale kinetic models for various biological domains such as signaling, gene expression and metabolism. Furthermore, we demonstrate how this toolbox is used to parameterize kinetic models around a steady-state reference efficiently. Finally, we show how SKiMpy can implement multispecies bioreactor simulations to assess biotechnological processes. AVAILABILITY AND IMPLEMENTATION: The software is available as a Python 3 package on GitHub: https://github.com/EPFL-LCSB/SKiMpy, along with adequate documentation. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Dynamic partitioning of branched-chain amino acids-derived nitrogen supports renal cancer progression.

Metabolic reprogramming is critical for tumor initiation and progression. However, the exact impact of specific metabolic changes on cancer progression is poorly understood. Here, we integrate multimodal analyses of primary and metastatic clonally-related clear cell renal cancer cells (ccRCC) grown in physiological media to identify key stage-specific metabolic vulnerabilities. We show that a VHL loss-dependent reprogramming of branched-chain amino acid catabolism sustains the de novo biosynthesis of aspartate and arginine enabling tumor cells with the flexibility of partitioning the nitrogen of the amino acids depending on their needs. Importantly, we identify the epigenetic reactivation of argininosuccinate synthase (ASS1), a urea cycle enzyme suppressed in primary ccRCC, as a crucial event for metastatic renal cancer cells to acquire the capability to generate arginine, invade in vitro and metastasize in vivo. Overall, our study uncovers a mechanism of metabolic flexibility occurring during ccRCC progression, paving the way for the development of novel stage-specific therapies.

A workflow for annotating the knowledge gaps in metabolic reconstructions using known and hypothetical reactions.

Advances in medicine and biotechnology rely on a deep understanding of biological processes. Despite the increasingly available types and amounts of omics data, significant knowledge gaps remain, with current approaches to identify and curate missing annotations being limited to a set of already known reactions. Here, we introduce Network Integrated Computational Explorer for Gap Annotation of Metabolism (NICEgame), a workflow to identify and curate nonannotated metabolic functions in genomes using the ATLAS of Biochemistry and genome-scale metabolic models (GEMs). To resolve gaps in GEMs, NICEgame provides alternative sets of known and hypothetical reactions, assesses their thermodynamic feasibility, and suggests candidate genes to catalyze these reactions. We identified metabolic gaps and applied NICEgame in the latest GEM of Escherichia coli, iML1515, and enhanced the E. coli genome annotation by resolving 47% of these gaps. NICEgame, applicable to any GEM and functioning from open-source software, should thus enhance all GEM-based predictions and subsequent biotechnological and biomedical applications.

ARBRE: Computational resource to predict pathways towards industrially important aromatic compounds.

Synthetic biology and metabolic engineering rely on computational search tools for predictions of novel biosynthetic pathways to industrially important compounds, many of which are derived from aromatic amino acids. Pathway search tools vary in their scope of covered reactions and compounds, as well as in metrics for ranking and evaluation. In this work, we present a new computational resource called ARBRE: Aromatic compounds RetroBiosynthesis Repository and Explorer. It consists of a comprehensive biochemical reaction network centered around aromatic amino acid biosynthesis and a computational toolbox for navigating this network. ARBRE encompasses over 33'000 known and 390'000 novel reactions predicted with generalized enzymatic reactions rules and over 74'000 compounds, of which 19'000 are known to biochemical databases and 55'000 only to PubChem. Over 1'000 molecules that were solely part of the PubChem database before and were previously impossible to integrate into a biochemical network are included into the ARBRE reaction network by assigning enzymatic reactions. ARBRE can be applied for pathway search, enzyme annotation, pathway ranking, visualization, and network expansion around known biochemical pathways and products of lignin degradation to predict valuable compound derivations. In line with the standards of open science, we have made the toolbox freely available to the scientific community on git (https://github.com/EPFL-LCSB/ARBRE) and we provide the web-version at http://lcsb-databases.epfl.ch/arbre/. We envision that ARBRE will provide the community with a new computational resource and comprehensive search tool to predict and rank pathways towards industrially important aromatic compounds.

Computational tools and resources for designing new pathways to small molecules.

The metabolic engineering community relies on computational methods for pathway design to produce important small molecules in microbial hosts. Metabolic network databases are continuously curated and updated with known and novel reactions that expand the known biochemistry based on different sets of enzymatic reaction rules. To address the complexity of the metabolic networks, elaborate methods were developed to transform them into computable graphs, navigate them, and construct the best possible pathways. However, the recent experimental research points to the new challenges and opportunities for the computational pathway design. Here, we review the most recent advances, especially in the last two years, in computational discovery of new pathways and their prospects for expanding metabolic capabilities. We draw attention to the potential ways of improvement for pathway design algorithms, including the expansion of Design-Build-Test-Learn cycle to novel compounds and reactions and the standardization for the reaction rules and metabolic reaction databases.

Expanding biochemical knowledge and illuminating metabolic dark matter with ATLASx.

Metabolic "dark matter" describes currently unknown metabolic processes, which form a blind spot in our general understanding of metabolism and slow down the development of biosynthetic cell factories and naturally derived pharmaceuticals. Mapping the dark matter of metabolism remains an open challenge that can be addressed globally and systematically by existing computational solutions. In this work, we use 489 generalized enzymatic reaction rules to map both known and unknown metabolic processes around a biochemical database of 1.5 million biological compounds. We predict over 5 million reactions and integrate nearly 2 million naturally and synthetically-derived compounds into the global network of biochemical knowledge, named ATLASx. ATLASx is available to researchers as a powerful online platform that supports the prediction and analysis of biochemical pathways and evaluates the biochemical vicinity of molecule classes ( https://lcsb-databases.epfl.ch/Atlas2 ).

Reconstructing Kinetic Models for Dynamical Studies of Metabolism using Generative Adversarial Networks.

Kinetic models of metabolism relate metabolic fluxes, metabolite concentrations and enzyme levels through mechanistic relations, rendering them essential for understanding, predicting and optimizing the behaviour of living organisms. However, due to the lack of kinetic data, traditional kinetic modelling often yields only a few or no kinetic models with desirable dynamical properties, making the analysis unreliable and computationally inefficient. We present REKINDLE (Reconstruction of Kinetic Models using Deep Learning), a deep-learning-based framework for efficiently generating kinetic models with dynamic properties matching the ones observed in cells. We showcase REKINDLE's capabilities to navigate through the physiological states of metabolism using small numbers of data with significantly lower computational requirements. The results show that data-driven neural networks assimilate implicit kinetic knowledge and structure of metabolic networks and generate kinetic models with tailored properties and statistical diversity. We anticipate that our framework will advance our understanding of metabolism and accelerate future research in biotechnology and health.

NICEdrug.ch, a workflow for rational drug design and systems-level analysis of drug metabolism.

The discovery of a drug requires over a decade of intensive research and financial investments - and still has a high risk of failure. To reduce this burden, we developed the NICEdrug.ch resource, which incorporates 250,000 bioactive molecules, and studied their enzymatic metabolic targets, fate, and toxicity. NICEdrug.ch includes a unique fingerprint that identifies reactive similarities between drug-drug and drug-metabolite pairs. We validated the application, scope, and performance of NICEdrug.ch over similar methods in the field on golden standard datasets describing drugs and metabolites sharing reactivity, drug toxicities, and drug targets. We use NICEdrug.ch to evaluate inhibition and toxicity by the anticancer drug 5-fluorouracil, and suggest avenues to alleviate its side effects. We propose shikimate 3-phosphate for targeting liver-stage malaria with minimal impact on the human host cell. Finally, NICEdrug.ch suggests over 1300 candidate drugs and food molecules to target COVID-19 and explains their inhibitory mechanism for further experimental screening. The NICEdrug.ch database is accessible online to systematically identify the reactivity of small molecules and druggable enzymes with practical applications in lead discovery and drug repurposing.

A genome-scale metabolic model of Saccharomyces cerevisiae that integrates expression constraints and reaction thermodynamics.

Eukaryotic organisms play an important role in industrial biotechnology, from the production of fuels and commodity chemicals to therapeutic proteins. To optimize these industrial systems, a mathematical approach can be used to integrate the description of multiple biological networks into a single model for cell analysis and engineering. One of the most accurate models of biological systems include Expression and Thermodynamics FLux (ETFL), which efficiently integrates RNA and protein synthesis with traditional genome-scale metabolic models. However, ETFL is so far only applicable for E. coli. To adapt this model for Saccharomyces cerevisiae, we developed yETFL, in which we augmented the original formulation with additional considerations for biomass composition, the compartmentalized cellular expression system, and the energetic costs of biological processes. We demonstrated the ability of yETFL to predict maximum growth rate, essential genes, and the phenotype of overflow metabolism. We envision that the presented formulation can be extended to a wide range of eukaryotic organisms to the benefit of academic and industrial research.

Spatio-temporal modeling of the crowding conditions and metabolic variability in microbial communities.

The metabolic capabilities of the species and the local environment shape the microbial interactions in a community either through the exchange of metabolic products or the competition for the resources. Cells are often arranged in close proximity to each other, creating a crowded environment that unevenly reduce the diffusion of nutrients. Herein, we investigated how the crowding conditions and metabolic variability among cells shape the dynamics of microbial communities. For this, we developed CROMICS, a spatio-temporal framework that combines techniques such as individual-based modeling, scaled particle theory, and thermodynamic flux analysis to explicitly incorporate the cell metabolism and the impact of the presence of macromolecular components on the nutrients diffusion. This framework was used to study two archetypical microbial communities (i) Escherichia coli and Salmonella enterica that cooperate with each other by exchanging metabolites, and (ii) two E. coli with different production level of extracellular polymeric substances (EPS) that compete for the same nutrients. In the mutualistic community, our results demonstrate that crowding enhanced the fitness of cooperative mutants by reducing the leakage of metabolites from the region where they are produced, avoiding the resource competition with non-cooperative cells. Moreover, we also show that E. coli EPS-secreting mutants won the competition against the non-secreting cells by creating less dense structures (i.e. increasing the spacing among the cells) that allow mutants to expand and reach regions closer to the nutrient supply point. A modest enhancement of the relative fitness of EPS-secreting cells over the non-secreting ones were found when the crowding effect was taken into account in the simulations. The emergence of cell-cell interactions and the intracellular conflicts arising from the trade-off between growth and the secretion of metabolites or EPS could provide a local competitive advantage to one species, either by supplying more cross-feeding metabolites or by creating a less dense neighborhood.

The influence of the crowding assumptions in biofilm simulations.

Microorganisms are frequently organized into crowded structures that affect the nutrients diffusion. This reduction in metabolite diffusion could modify the microbial dynamics, meaning that computational methods for studying microbial systems need accurate ways to model the crowding conditions. We previously developed a computational framework, termed CROMICS, that incorporates the effect of the (time-dependent) crowding conditions on the spatio-temporal modeling of microbial communities, and we used it to demonstrate the crowding influence on the community dynamics. To further identify scenarios where crowding should be considered in microbial modeling, we herein applied and extended CROMICS to simulate several environmental conditions that could potentially boost or dampen the crowding influence in biofilms. We explore whether the nutrient supply (rich- or low-nutrient media), the cell-packing configuration (square or hexagonal spherical cell arrangement), or the cell growing conditions (planktonic state or biofilm) modify the crowding influence on the growth of Escherichia coli. Our results indicate that the growth rate, the abundance and appearance time of different cell phenotypes as well as the amount of by-products secreted to the medium are sensitive to some extent to the local crowding conditions in all scenarios tested, except in rich-nutrient media. Crowding conditions enhance the formation of nutrient gradient in biofilms, but its effect is only appreciated when cell metabolism is controlled by the nutrient limitation. Thus, as soon as biomass (and/or any other extracellular macromolecule) accumulates in a region, and cells occupy more than 14% of the volume fraction, the crowding effect must not be underestimated, as the microbial dynamics start to deviate from the ideal/expected behaviour that assumes volumeless cells or when a homogeneous (reduced) diffusion is applied in the simulation. The modeling and simulation of the interplay between the species diversity (cell shape and metabolism) and the environmental conditions (nutrient quality, crowding conditions) can help to design effective strategies for the optimization and control of microbial systems.

NICEpath: Finding metabolic pathways in large networks through atom-conserving substrate-product pairs.

MOTIVATION: Finding biosynthetic pathways is essential for metabolic engineering of organisms to produce chemicals, biodegradation prediction of pollutants and drugs, and for the elucidation of bioproduction pathways of secondary metabolites. A key step in biosynthetic pathway design is the extraction of novel metabolic pathways from big networks that integrate known biological, as well as novel, predicted biotransformations. However, the efficient analysis and the navigation of big biochemical networks remain a challenge. RESULTS: Here, we propose the construction of searchable graph representations of metabolic networks. Each reaction is decomposed into pairs of reactants and products, and each pair is assigned a weight, which is calculated from the number of conserved atoms between the reactant and the product molecule. We test our method on a biochemical network that spans 6546 known enzymatic reactions to show how our approach elegantly extracts biologically relevant metabolic pathways from biochemical networks, and how the proposed network structure enables the application of efficient graph search algorithms that improve navigation and pathway identification in big metabolic networks. The weighted reactant-product pairs of an example network and the corresponding graph search algorithm are available online. The proposed method extracts metabolic pathways fast and reliably from big biochemical networks, which is inherently important for all applications involving the engineering of metabolic networks. AVAILABILITY AND IMPLEMENTATION: https://github.com/EPFL-LCSB/nicepath. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Constraint-based metabolic control analysis for rational strain engineering.

The advancements in genome editing techniques over the past years have rekindled interest in rational metabolic engineering strategies. While Metabolic Control Analysis (MCA) is a well-established method for quantifying the effects of metabolic engineering interventions on flows in metabolic networks and metabolite concentrations, it does not consider the physiological limitations of the cellular environment and metabolic engineering design constraints. We report here a constraint-based framework, Network Response Analysis (NRA), for rational genetic strain design. NRA is cast as a Mixed-Integer Linear Programming problem that integrates MCA, Thermodynamically-based Flux Analysis (TFA), biologically relevant constraints, as well as genome editing restrictions into a comprehensive platform for identifying metabolic engineering targets. We show that the NRA formulation and its core constraints are equivalent to the ones of Flux Balance Analysis (FBA) and TFA, which allows it to be used for a wide range of optimization criteria and with various physiological constraints. We also show how the parametrization and introduction of biological constraints enhance the NRA formulation compared to the classical MCA approach, and we demonstrate its features and its ability to generate multiple alternative optimal strategies given several user-defined boundaries and objectives. In summary, NRA is a sophisticated alternative to classical MCA for rational metabolic engineering that accommodates the incorporation of physiological data at metabolic flux, metabolite concentration, and enzyme expression levels.

The effects of model complexity and size on metabolic flux distribution and control: case study in Escherichia coli.

BACKGROUND: Significant efforts have been made in building large-scale kinetic models of cellular metabolism in the past two decades. However, most kinetic models published to date, remain focused around central carbon pathways or are built around ad hoc reduced models without clear justification on their derivation and usage. Systematic algorithms exist for reducing genome-scale metabolic reconstructions to build thermodynamically feasible and consistently reduced stoichiometric models. However, it is important to study how network complexity affects conclusions derived from large-scale kinetic models built around consistently reduced models before we can apply them to study biological systems. RESULTS: We reduced the iJO1366 Escherichia Coli genome-scale metabolic reconstruction systematically to build three stoichiometric models of different size. Since the reduced models are expansions around the core subsystems for which the reduction was performed, the models are nested. We present a method for scaling up the flux profile and the concentration vector reference steady-states from the smallest model to the larger ones, whilst preserving maximum equivalency. Populations of kinetic models, preserving similarity in kinetic parameters, were built around the reference steady-states and their metabolic sensitivity coefficients (MSCs) were computed. The MSCs were sensitive to the model complexity. We proposed a metric for measuring the sensitivity of MSCs to these structural changes. CONCLUSIONS: We proposed for the first time a workflow for scaling up the size of kinetic models while preserving equivalency between the kinetic models. Using this workflow, we demonstrate that model complexity in terms of networks size has significant impact on sensitivity characteristics of kinetic models. Therefore, it is essential to account for the effects of network complexity when constructing kinetic models. The presented metric for measuring MSC sensitivity to structural changes can guide modelers and experimentalists in improving model quality and guide synthetic biology and metabolic engineering. Our proposed workflow enables the testing of the suitability of a kinetic model for answering certain study-specific questions. We argue that the model-based metabolic design targets that are common across models of different size are of higher confidence, while those that are different could be the objective of investigations for model improvement.

A computational workflow for the expansion of heterologous biosynthetic pathways to natural product derivatives.

Plant natural products (PNPs) and their derivatives are important but underexplored sources of pharmaceutical molecules. To access this untapped potential, the reconstitution of heterologous PNP biosynthesis pathways in engineered microbes provides a valuable starting point to explore and produce novel PNP derivatives. Here, we introduce a computational workflow to systematically screen the biochemical vicinity of a biosynthetic pathway for pharmaceutical compounds that could be produced by derivatizing pathway intermediates. We apply our workflow to the biosynthetic pathway of noscapine, a benzylisoquinoline alkaloid (BIA) with a long history of medicinal use. Our workflow identifies pathways and enzyme candidates for the production of (S)-tetrahydropalmatine, a known analgesic and anxiolytic, and three additional derivatives. We then construct pathways for these compounds in yeast, resulting in platforms for de novo biosynthesis of BIA derivatives and demonstrating the value of cheminformatic tools to predict reactions, pathways, and enzymes in synthetic biology and metabolic engineering.