Resource allocation modeling for autonomous prediction of plant cell phenotypes.

Predicting the plant cell response in complex environmental conditions is a challenge in plant biology. Here we developed a resource allocation model of cellular and molecular scale for the leaf photosynthetic cell of Arabidopsis thaliana, based on the Resource Balance Analysis (RBA) constraint-based modeling framework. The RBA model contains the metabolic network and the major macromolecular processes involved in the plant cell growth and survival and localized in cellular compartments. We simulated the model for varying environmental conditions of temperature, irradiance, partial pressure of CO2 and O2, and compared RBA predictions to known resource distributions and quantitative phenotypic traits such as the relative growth rate, the C:N ratio, and finally to the empirical characteristics of CO2 fixation given by the well-established Farquhar model. In comparison to other standard constraint-based modeling methods like Flux Balance Analysis, the RBA model makes accurate quantitative predictions without the need for empirical constraints. Altogether, we show that RBA significantly improves the autonomous prediction of plant cell phenotypes in complex environmental conditions, and provides mechanistic links between the genotype and the phenotype of the plant cell.

BioSimulators: a central registry of simulation engines and services for recommending specific tools.

Computational models have great potential to accelerate bioscience, bioengineering, and medicine. However, it remains challenging to reproduce and reuse simulations, in part, because the numerous formats and methods for simulating various subsystems and scales remain siloed by different software tools. For example, each tool must be executed through a distinct interface. To help investigators find and use simulation tools, we developed BioSimulators (https://biosimulators.org), a central registry of the capabilities of simulation tools and consistent Python, command-line and containerized interfaces to each version of each tool. The foundation of BioSimulators is standards, such as CellML, SBML, SED-ML and the COMBINE archive format, and validation tools for simulation projects and simulation tools that ensure these standards are used consistently. To help modelers find tools for particular projects, we have also used the registry to develop recommendation services. We anticipate that BioSimulators will help modelers exchange, reproduce, and combine simulations.

BiPOm: a rule-based ontology to represent and infer molecule knowledge from a biological process-centered viewpoint.

BACKGROUND: Managing and organizing biological knowledge remains a major challenge, due to the complexity of living systems. Recently, systemic representations have been promising in tackling such a challenge at the whole-cell scale. In such representations, the cell is considered as a system composed of interlocked subsystems. The need is now to define a relevant formalization of the systemic description of cellular processes. RESULTS: We introduce BiPOm (Biological interlocked Process Ontology for metabolism) an ontology to represent metabolic processes as interlocked subsystems using a limited number of classes and properties. We explicitly formalized the relations between the enzyme, its activity, the substrates and the products of the reaction, as well as the active state of all involved molecules. We further showed that the information of molecules such as molecular types or molecular properties can be deduced by automatic reasoning using logical rules. The information necessary to populate BiPOm can be extracted from existing databases or existing bio-ontologies. CONCLUSION: BiPOm provides a formal rule-based knowledge representation to relate all cellular components together by considering the cellular system as a whole. It relies on a paradigm shift where the anchorage of knowledge is rerouted from the molecule to the biological process. AVAILABILITY: BiPOm can be downloaded at https://github.com/SysBioInra/SysOnto.

Automated generation of bacterial resource allocation models.

Resource Balance Analysis (RBA) is a computational method based on resource allocation, which performs accurate quantitative predictions of whole-cell states (i.e. growth rate, metabolic fluxes, abundances of molecular machines including enzymes) across growth conditions. We present an integrated workflow of RBA together with the Python package RBApy. RBApy builds bacterial RBA models from annotated genome-scale metabolic models by adding descriptions of cellular processes relevant for growth and maintenance. The package includes functions for model simulation and calibration and for interfacing to Escher maps and Proteomaps for visualization. We demonstrate that RBApy faithfully reproduces results obtained by a hand-curated and experimentally validated RBA model for Bacillus subtilis. We also present a calibrated RBA model of Escherichia coli generated from scratch, which obtained excellent fits to measured flux values and enzyme abundances. RBApy makes whole-cell modelling accessible for a wide range of bacterial wild-type and engineered strains, as illustrated with a CO2-fixing Escherichia coli strain. AVAILABILITY: RBApy is available at /https://github.com/SysBioInra/RBApy, under the licence GNU GPL version 3, and runs on Linux, Mac and Windows distributions.

Resource allocation in living organisms.

Quantitative prediction of resource allocation for living systems has been an intensive area of research in the field of biology. Resource allocation was initially investigated in higher organisms by using empirical mathematical models based on mass distribution. A challenge is now to go a step further by reconciling the cellular scale to the individual scale. In the present paper, we review the foundations of modelling of resource allocation, particularly at the cellular scale: from small macro-molecular models to genome-scale cellular models. We enlighten how the combination of omic measurements and computational advances together with systems biology has contributed to dramatic progresses in the current understanding and prediction of cellular resource allocation. Accurate genome-wide predictive methods of resource allocation based on the resource balance analysis (RBA) framework have been developed and ensure a good trade-off between the complexity/tractability and the prediction capability of the model. The RBA framework shows promise for a wide range of applications in metabolic engineering and synthetic biology, and for pursuing investigations of the design principles of cellular and multi-cellular organisms.

Translation elicits a growth rate-dependent, genome-wide, differential protein production in Bacillus subtilis.

Complex regulatory programs control cell adaptation to environmental changes by setting condition-specific proteomes. In balanced growth, bacterial protein abundances depend on the dilution rate, transcript abundances and transcript-specific translation efficiencies. We revisited the current theory claiming the invariance of bacterial translation efficiency. By integrating genome-wide transcriptome datasets and datasets from a library of synthetic gfp-reporter fusions, we demonstrated that translation efficiencies in Bacillus subtilis decreased up to fourfold from slow to fast growth. The translation initiation regions elicited a growth rate-dependent, differential production of proteins without regulators, hence revealing a unique, hard-coded, growth rate-dependent mode of regulation. We combined model-based data analyses of transcript and protein abundances genome-wide and revealed that this global regulation is extensively used in B. subtilis We eventually developed a knowledge-based, three-step translation initiation model, experimentally challenged the model predictions and proposed that a growth rate-dependent drop in free ribosome abundance accounted for the differential protein production.

Optimal resource allocation enables mathematical exploration of microbial metabolic configurations.

Central to the functioning of any living cell, the metabolic network is a complex network of biochemical reactions. It may also be viewed as an elaborate production system, integrating a diversity of internal and external signals in order to efficiently produce the energy and the biochemical precursors to ensure all cellular functions. Even in simple organisms like bacteria, it shows a striking level of coordination, adapting to very different growth media. Constraint-based models constitute an efficient mathematical framework to compute optimal metabolic configurations, at the scale of a whole genome. Combining the constraint-based approach "Resource Balance Analysis" with combinatorial optimization techniques, we propose a general method to explore these configurations, based on the inference of logical rules governing the activation of metabolic fluxes in response to diverse extracellular media. Using the concept of partial Boolean functions, we notably introduce a novel tractable algorithm to infer monotone Boolean functions on a minimal support. Monotonicity seems particularly relevant in this context, since the orderliness exhibited by the metabolic network's dynamical behavior is expected to give rise to relatively simple rules. First results are promising, as the application of the method on Bacillus subtilis central carbon metabolism allows to recover known regulations as well as to investigate lesser known parts of the global regulatory network.

Comprehensive absolute quantification of the cytosolic proteome of Bacillus subtilis by data independent, parallel fragmentation in liquid chromatography/mass spectrometry (LC/MS(E)).

In the growing field of systems biology, the knowledge of protein concentrations is highly required to truly understand metabolic and adaptational networks within the cells. Therefore we established a workflow relying on long chromatographic separation and mass spectrometric analysis by data independent, parallel fragmentation of all precursor ions at the same time (LC/MS(E)). By prevention of discrimination of co-eluting low and high abundant peptides a high average sequence coverage of 40% could be achieved, resulting in identification of almost half of the predicted cytosolic proteome of the Gram-positive model organism Bacillus subtilis (>1,050 proteins). Absolute quantification was achieved by correlation of average MS signal intensities of the three most intense peptides of a protein to the signal intensity of a spiked standard protein digest. Comparative analysis with heavily labeled peptides (AQUA approach) showed the use of only one standard digest is sufficient for global quantification. The quantification results covered almost four orders of magnitude, ranging roughly from 10 to 150,000 copies per cell. To prove this method for its biological relevance selected physiological aspects of B. subtilis cells grown under conditions requiring either amino acid synthesis or alternatively amino acid degradation were analyzed. This allowed both in particular the validation of the adjustment of protein levels by known regulatory events and in general a perspective of new insights into bacterial physiology. Within new findings the analysis of "protein costs" of cellular processes is extremely important. Such a comprehensive and detailed characterization of cellular protein concentrations based on data independent, parallel fragmentation in liquid chromatography/mass spectrometry (LC/MS(E)) data has been performed for the first time and should pave the way for future comprehensive quantitative characterization of microorganisms as physiological entities.

Quantitative prediction of genome-wide resource allocation in bacteria.

Predicting resource allocation between cell processes is the primary step towards decoding the evolutionary constraints governing bacterial growth under various conditions. Quantitative prediction at genome-scale remains a computational challenge as current methods are limited by the tractability of the problem or by simplifying hypotheses. Here, we show that the constraint-based modeling method Resource Balance Analysis (RBA), calibrated using genome-wide absolute protein quantification data, accurately predicts resource allocation in the model bacterium Bacillus subtilis for a wide range of growth conditions. The regulation of most cellular processes is consistent with the objective of growth rate maximization except for a few suboptimal processes which likely integrate more complex objectives such as coping with stressful conditions and survival. As a proof of principle by using simulations, we illustrated how calibrated RBA could aid rational design of strains for maximizing protein production, offering new opportunities to investigate design principles in prokaryotes and to exploit them for biotechnological applications.

BasyLiCA: a tool for automatic processing of a Bacterial Live Cell Array.

UNLABELLED: Live Cell Array (LCA) technology allows the acquisition of high-resolution time-course profiles of bacterial gene expression by the systematic assessment of fluorescence in living cells carrying either transcriptional or translational fluorescent protein fusion. However, the direct estimation of promoter activities by time-dependent derivation of the fluorescence datasets generates high levels of noise. Here, we present BasyLiCA, a user-friendly open-source interface and database dedicated to the automatic storage and standardized treatment of LCA data. Data quality reports are generated automatically. Growth rates and promoter activities are calculated by tunable discrete Kalman filters that can be set to incorporate data from biological replicates, significantly reducing the impact of noise measurement in activity estimations. AVAILABILITY: The BasyLiCA software and the related documentation are available at http://genome.jouy.inra.fr/basylica.

RBAtools: a programming interface for Resource Balance Analysis models.

Motivation: Efficient resource allocation can contribute to an organism's fitness and can improve evolutionary success. Resource Balance Analysis (RBA) is a computational framework that models an organism's growth-optimal proteome configurations in various environments. RBA software enables the construction of RBA models on genome scale and the calculation of medium-specific, growth-optimal cell states including metabolic fluxes and the abundance of macromolecular machines. However, existing software lacks a simple programming interface for non-expert users, easy to use and interoperable with other software. Results: The python package RBAtools provides convenient access to RBA models. As a flexible programming interface, it enables the implementation of custom workflows and the modification of existing genome-scale RBA models. Its high-level functions comprise simulation, model fitting, parameter screens, sensitivity analysis, variability analysis and the construction of Pareto fronts. Models and data are represented as structured tables and can be exported to common data formats for fluxomics and proteomics visualization. Availability and implementation: RBAtools documentation, installation instructions and tutorials are available at https://sysbioinra.github.io/rbatools/. General information about RBA and related software can be found at rba.inrae.fr.

A comparative transcriptomic, fluxomic and metabolomic analysis of the response of Saccharomyces cerevisiae to increases in NADPH oxidation.

BACKGROUND: Redox homeostasis is essential to sustain metabolism and growth. We recently reported that yeast cells meet a gradual increase in imposed NADPH demand by progressively increasing flux through the pentose phosphate (PP) and acetate pathways and by exchanging NADH for NADPH in the cytosol, via a transhydrogenase-like cycle. Here, we studied the mechanisms underlying this metabolic response, through a combination of gene expression profiling and analyses of extracellular and intracellular metabolites and 13 C-flux analysis. RESULTS: NADPH oxidation was increased by reducing acetoin to 2,3-butanediol in a strain overexpressing an engineered NADPH-dependent butanediol dehydrogenase cultured in the presence of acetoin. An increase in NADPH demand to 22 times the anabolic requirement for NADPH was accompanied by the intracellular accumulation of PP pathway metabolites consistent with an increase in flux through this pathway. Increases in NADPH demand were accompanied by the successive induction of several genes of the PP pathway. NADPH-consuming pathways, such as amino-acid biosynthesis, were upregulated as an indirect effect of the decrease in NADPH availability. Metabolomic analysis showed that the most extreme modification of NADPH demand resulted in an energetic problem. Our results also highlight the influence of redox status on aroma production. CONCLUSIONS: Combined 13 C-flux, intracellular metabolite levels and microarrays analyses revealed that NADPH homeostasis, in response to a progressive increase in NADPH demand, was achieved by the regulation, at several levels, of the PP pathway. This pathway is principally under metabolic control, but regulation of the transcription of PP pathway genes can exert a stronger effect, by redirecting larger amounts of carbon to this pathway to satisfy the demand for NADPH. No coordinated response of genes involved in NADPH metabolism was observed, suggesting that yeast has no system for sensing NADPH/NADP+ ratio. Instead, the induction of NADPH-consuming amino-acid pathways in conditions of NADPH limitation may indirectly trigger the transcription of a set of PP pathway genes.

A constraint-based model analysis of the metabolic consequences of increased NADPH oxidation in Saccharomyces cerevisiae.

Controlling the amounts of redox cofactors to manipulate metabolic fluxes is emerging as a useful approach to optimizing byproduct yields in yeast biotechnological processes. Redox cofactors are extensively interconnected metabolites, so predicting metabolite patterns is challenging and requires in-depth knowledge of how the metabolic network responds to a redox perturbation. Our aim was to analyze comprehensively the metabolic consequences of increased cytosolic NADPH oxidation during yeast fermentation. Using a genetic device based on the overexpression of a modified 2,3-butanediol dehydrogenase catalyzing the NADPH-dependent reduction of acetoin into 2,3-butanediol, we increased the NADPH demand to between 8 and 40-fold the anabolic demand. We developed (i) a dedicated constraint-based model of yeast fermentation and (ii) a constraint-based modeling method based on the dynamical analysis of mass distribution to quantify the in vivo contribution of pathways producing NADPH to the maintenance of redox homeostasis. We report that yeast responds to NADPH oxidation through a gradual increase in the flux through the PP and acetate pathways, providing 80% and 20% of the NADPH demand, respectively. However, for the highest NADPH demand, the model reveals a saturation of the PP pathway and predicts an exchange between NADH and NADPH in the cytosol that may be mediated by the glycerol-DHA futile cycle. We also reveal the contribution of mitochondrial shuttles, resulting in a net production of NADH in the cytosol, to fine-tune the NADH/NAD(+) balance. This systems level study helps elucidate the physiological adaptation of yeast to NADPH perturbation. Our findings emphasize the robustness of yeast to alterations in NADPH metabolism and highlight the role of the glycerol-DHA cycle as a redox valve, providing additional NADPH from NADH under conditions of very high demand.

BiPSim: a flexible and generic stochastic simulator for polymerization processes.

Detailed whole-cell modeling requires an integration of heterogeneous cell processes having different modeling formalisms, for which whole-cell simulation could remain tractable. Here, we introduce BiPSim, an open-source stochastic simulator of template-based polymerization processes, such as replication, transcription and translation. BiPSim combines an efficient abstract representation of reactions and a constant-time implementation of the Gillespie's Stochastic Simulation Algorithm (SSA) with respect to reactions, which makes it highly efficient to simulate large-scale polymerization processes stochastically. Moreover, multi-level descriptions of polymerization processes can be handled simultaneously, allowing the user to tune a trade-off between simulation speed and model granularity. We evaluated the performance of BiPSim by simulating genome-wide gene expression in bacteria for multiple levels of granularity. Finally, since no cell-type specific information is hard-coded in the simulator, models can easily be adapted to other organismal species. We expect that BiPSim should open new perspectives for the genome-wide simulation of stochastic phenomena in biology.

The bacterial interlocked process ONtology (BiPON): a systemic multi-scale unified representation of biological processes in prokaryotes.

BACKGROUND: High-throughput technologies produce huge amounts of heterogeneous biological data at all cellular levels. Structuring these data together with biological knowledge is a critical issue in biology and requires integrative tools and methods such as bio-ontologies to extract and share valuable information. In parallel, the development of recent whole-cell models using a systemic cell description opened alternatives for data integration. Integrating a systemic cell description within a bio-ontology would help to progress in whole-cell data integration and modeling synergistically. RESULTS: We present BiPON, an ontology integrating a multi-scale systemic representation of bacterial cellular processes. BiPON consists in of two sub-ontologies, bioBiPON and modelBiPON. bioBiPON organizes the systemic description of biological information while modelBiPON describes the mathematical models (including parameters) associated with biological processes. bioBiPON and modelBiPON are related using bridge rules on classes during automatic reasoning. Biological processes are thus automatically related to mathematical models. 37% of BiPON classes stem from different well-established bio-ontologies, while the others have been manually defined and curated. Currently, BiPON integrates the main processes involved in bacterial gene expression processes. CONCLUSIONS: BiPON is a proof of concept of the way to combine formally systems biology and bio-ontology. The knowledge formalization is highly flexible and generic. Most of the known cellular processes, new participants or new mathematical models could be inserted in BiPON. Altogether, BiPON opens up promising perspectives for knowledge integration and sharing and can be used by biologists, systems and computational biologists, and the emerging community of whole-cell modeling.

Conversion of Glycerol to 3-Hydroxypropanoic Acid by Genetically Engineered Bacillus subtilis.

3-Hydroxypropanoic acid (3-HP) is an important biomass-derivable platform chemical that can be converted into a number of industrially relevant compounds. There have been several attempts to produce 3-HP from renewable sources in cell factories, focusing mainly on Escherichia coli, Klebsiella pneumoniae, and Saccharomyces cerevisiae. Despite the significant progress made in this field, commercially exploitable large-scale production of 3-HP in microbial strains has still not been achieved. In this study, we investigated the potential of Bacillus subtilis as a microbial platform for bioconversion of glycerol into 3-HP. Our recombinant B. subtilis strains overexpress the two-step heterologous pathway containing glycerol dehydratase and aldehyde dehydrogenase from K. pneumoniae. Genetic engineering, driven by in silico optimization, and optimization of cultivation conditions resulted in a 3-HP titer of 10 g/L, in a standard batch cultivation. Our findings provide the first report of successful introduction of the biosynthetic pathway for conversion of glycerol into 3-HP in B. subtilis. With this relatively high titer in batch, and the robustness of B. subtilis in high density fermentation conditions, we expect that our production strains may constitute a solid basis for commercial production of 3-HP.

Reconstruction of the Regulatory Network for Bacillus subtilis and Reconciliation with Gene Expression Data.

We introduce a manually constructed and curated regulatory network model that describes the current state of knowledge of transcriptional regulation of Bacillus subtilis. The model corresponds to an updated and enlarged version of the regulatory model of central metabolism originally proposed in 2008. We extended the original network to the whole genome by integration of information from DBTBS, a compendium of regulatory data that includes promoters, transcription factors (TFs), binding sites, motifs, and regulated operons. Additionally, we consolidated our network with all the information on regulation included in the SporeWeb and Subtiwiki community-curated resources on B. subtilis. Finally, we reconciled our network with data from RegPrecise, which recently released their own less comprehensive reconstruction of the regulatory network for B. subtilis. Our model describes 275 regulators and their target genes, representing 30 different mechanisms of regulation such as TFs, RNA switches, Riboswitches, and small regulatory RNAs. Overall, regulatory information is included in the model for ∼2500 of the ∼4200 genes in B. subtilis 168. In an effort to further expand our knowledge of B. subtilis regulation, we reconciled our model with expression data. For this process, we reconstructed the Atomic Regulons (ARs) for B. subtilis, which are the sets of genes that share the same "ON" and "OFF" gene expression profiles across multiple samples of experimental data. We show how ARs for B. subtilis are able to capture many sets of genes corresponding to regulated operons in our manually curated network. Additionally, we demonstrate how ARs can be used to help expand or validate the knowledge of the regulatory networks by looking at highly correlated genes in the ARs for which regulatory information is lacking. During this process, we were also able to infer novel stimuli for hypothetical genes by exploring the genome expression metadata relating to experimental conditions, gaining insights into novel biology.

Global network reorganization during dynamic adaptations of Bacillus subtilis metabolism.

Adaptation of cells to environmental changes requires dynamic interactions between metabolic and regulatory networks, but studies typically address only one or a few layers of regulation. For nutritional shifts between two preferred carbon sources of Bacillus subtilis, we combined statistical and model-based data analyses of dynamic transcript, protein, and metabolite abundances and promoter activities. Adaptation to malate was rapid and primarily controlled posttranscriptionally compared with the slow, mainly transcriptionally controlled adaptation to glucose that entailed nearly half of the known transcription regulation network. Interactions across multiple levels of regulation were involved in adaptive changes that could also be achieved by controlling single genes. Our analysis suggests that global trade-offs and evolutionary constraints provide incentives to favor complex control programs.

Condition-dependent transcriptome reveals high-level regulatory architecture in Bacillus subtilis.

Bacteria adapt to environmental stimuli by adjusting their transcriptomes in a complex manner, the full potential of which has yet to be established for any individual bacterial species. Here, we report the transcriptomes of Bacillus subtilis exposed to a wide range of environmental and nutritional conditions that the organism might encounter in nature. We comprehensively mapped transcription units (TUs) and grouped 2935 promoters into regulons controlled by various RNA polymerase sigma factors, accounting for ~66% of the observed variance in transcriptional activity. This global classification of promoters and detailed description of TUs revealed that a large proportion of the detected antisense RNAs arose from potentially spurious transcription initiation by alternative sigma factors and from imperfect control of transcription termination.

Reconstruction and analysis of the genetic and metabolic regulatory networks of the central metabolism of Bacillus subtilis.

BACKGROUND: Few genome-scale models of organisms focus on the regulatory networks and none of them integrates all known levels of regulation. In particular, the regulations involving metabolite pools are often neglected. However, metabolite pools link the metabolic to the genetic network through genetic regulations, including those involving effectors of transcription factors or riboswitches. Consequently, they play pivotal roles in the global organization of the genetic and metabolic regulatory networks. RESULTS: We report the manually curated reconstruction of the genetic and metabolic regulatory networks of the central metabolism of Bacillus subtilis (transcriptional, translational and post-translational regulations and modulation of enzymatic activities). We provide a systematic graphic representation of regulations of each metabolic pathway based on the central role of metabolites in regulation. We show that the complex regulatory network of B. subtilis can be decomposed as sets of locally regulated modules, which are coordinated by global regulators. CONCLUSION: This work reveals the strong involvement of metabolite pools in the general regulation of the metabolic network. Breaking the metabolic network down into modules based on the control of metabolite pools reveals the functional organization of the genetic and metabolic regulatory networks of B. subtilis.