Development of thermodynamic optimum searching (TOS) to improve the prediction accuracy of flux balance analysis

Deciphering and designing microbial communities by genome-scale metabolic modelling

Microbial communities are shaped by the complex interactions among organisms and the environment. Genome-scale metabolic models (GEMs) can provide deeper insights into the complexity and ecological properties of various microbial communities, revealing their intricate interactions. Many researchers have modified GEMs for the microbial communities based on specific needs. Thus, GEMs need to be comprehensively summarized to better understand the trends in their development. In this review, we summarized the key developments in deciphering and designing microbial communities using different GEMs. A timeline of selected highlights in GEMs indicated that this area is evolving from the single-strain level to the microbial community level. Then, we outlined a framework for constructing GEMs of microbial communities. We also summarized the models and resources of static and dynamic community-level GEMs. We focused on the role of external environmental and intracellular resources in shaping the assembly of microbial communities. Finally, we discussed the key challenges and future directions of GEMs, focusing on the integration of GEMs with quorum sensing mechanisms, microbial ecology interactions, machine learning algorithms, and automatic modeling, all of which contribute to consortia-based applications in different fields.

Flux Balance Analysis of Mammalian Cell Systems

Flux balance analysis (FBA) is a computational methodology to model and analyze the metabolic behavior of cells. In this chapter, we break down the key steps for formulating an FBA model and other FBA-derived methodologies in the context of mammalian cell biology, including strain design, developing cell line-specific models, and conducting flux sampling. We provide annotated COBRApy code for each step to show how it would work in practice.

Model-driven insights into the effects of temperature on metabolism

Temperature affects cellular processes at different spatiotemporal scales, and identifying the genetic and molecular mechanisms underlying temperature responses paves the way to develop approaches for mitigating the effects of future climate scenarios. A systems view of the effects of temperature on cellular physiology can be obtained by focusing on metabolism since: (i) its functions depend on transcription and translation and (ii) its outcomes support organisms' development, growth, and reproduction. Here we provide a systematic review of modelling efforts directed at investigating temperature effects on properties of single biochemical reactions, system-level traits, metabolic subsystems, and whole-cell metabolism across different prokaryotes and eukaryotes. We compare and contrast computational approaches and theories that facilitate modelling of temperature effects on key properties of enzymes and their consideration in constraint-based as well as kinetic models of metabolism. In addition, we provide a summary of insights from computational approaches, facilitating integration of omics data from temperature-modulated experiments with models of metabolic networks, and review the resulting biotechnological applications. Lastly, we provide a perspective on how different types of metabolic modelling can profit from developments in machine learning and models of different cellular layers to improve model-driven insights into the effects of temperature relevant for biotechnological applications.

The view of microbes as energy converters illustrates the trade-off between growth rate and yield

The application of thermodynamics to microbial growth has a long tradition that originated in the middle of the 20th century. This approach reflects the view that self-replication is a thermodynamic process that is not fundamentally different from mechanical thermodynamics. The key distinction is that a free energy gradient is not converted into mechanical (or any other form of) energy but rather into new biomass. As such, microbes can be viewed as energy converters that convert a part of the energy contained in environmental nutrients into chemical energy that drives self-replication. Before the advent of high-throughput sequencing technologies, only the most central metabolic pathways were known. However, precise measurement techniques allowed for the quantification of exchanged extracellular nutrients and heat of growing microbes with their environment. These data, together with the absence of knowledge of metabolic details, drove the development of so-called black-box models, which only consider the observable interactions of a cell with its environment and neglect all details of how exactly inputs are converted into outputs. Now, genome sequencing and genome-scale metabolic models (GEMs) provide us with unprecedented detail about metabolic processes inside the cell. However, mostly due to computational complexity issues, the derived modelling approaches make surprisingly little use of thermodynamic concepts. Here, we review classical black-box models and modern approaches that integrate thermodynamics into GEMs. We also illustrate how the description of microbial growth as an energy converter can help to understand and quantify the trade-off between microbial growth rate and yield.

Designing microbial communities to maximize the thermodynamic driving force for the production of chemicals

Microbial communities have become a major research focus due to their importance for biogeochemical cycles, biomedicine and biotechnological applications. While some biotechnological applications, such as anaerobic digestion, make use of naturally arising microbial communities, the rational design of microbial consortia for bio-based production processes has recently gained much interest. One class of synthetic microbial consortia is based on specifically designed strains of one species. A common design principle for these consortia is based on division of labor, where the entire production pathway is divided between the different strains to reduce the metabolic burden caused by product synthesis. We first show that classical division of labor does not automatically reduce the metabolic burden when metabolic flux per biomass is analyzed. We then present ASTHERISC (Algorithmic Search of THERmodynamic advantages in Single-species Communities), a new computational approach for designing multi-strain communities of a single-species with the aim to divide a production pathway between different strains such that the thermodynamic driving force for product synthesis is maximized. ASTHERISC exploits the fact that compartmentalization of segments of a product pathway in different strains can circumvent thermodynamic bottlenecks arising when operation of one reaction requires a metabolite with high and operation of another reaction the same metabolite with low concentration. We implemented the ASTHERISC algorithm in a dedicated program package and applied it on E. coli core and genome-scale models with different settings, for example, regarding number of strains or demanded product yield. These calculations showed that, for each scenario, many target metabolites (products) exist where a multi-strain community can provide a thermodynamic advantage compared to a single strain solution. In some cases, a production with sufficiently high yield is thermodynamically only feasible with a community. In summary, the developed ASTHERISC approach provides a promising new principle for designing microbial communities for the bio-based production of chemicals.

Communities of microbes are ubiquitous in nature and also of high relevance for industrial applications, e.g. for the production of biogas. The development and use of non-natural communities for biotechnological applications has become an important subject of research. In this work, we present a new computational method to design synthetic communities with improved capabilities for the synthesis of desired target metabolites. Our method takes a constraint-based metabolic model of an organism as input and searches for a suitable partitioning of the product pathway via different strains of the organism such that the thermodynamic driving force for product synthesis is maximized. Essentially, this approach exploits the fact that having multiple strains allows adjustment of different metabolite concentrations in the different strains by which the thermodynamic driving force for product synthesis can often be increased. We tested this approach with a core and with a genome-scale metabolic network model of Escherichia coli. We found that, for dozens of metabolites, there exist communities with specifically designed strains of E. coli where the maximal thermodynamic driving force can be increased compared to a single E. coli strain. In summary, our presented method provides a new approach, together with a new design principle, for the computational design of microbial communities.

Low Entropy Sub-Networks Prevent the Integration of Metabolomic and Transcriptomic Data

The constantly and rapidly increasing amount of the biological data gained from many different high-throughput experiments opens up new possibilities for data- and model-driven inference. Yet, alongside, emerges a problem of risks related to data integration techniques. The latter are not so widely taken account of. Especially, the approaches based on the flux balance analysis (FBA) are sensitive to the structure of a metabolic network for which the low-entropy clusters can prevent the inference from the activity of the metabolic reactions. In the following article, we set forth problems that may arise during the integration of metabolomic data with gene expression datasets. We analyze common pitfalls, provide their possible solutions, and exemplify them by a case study of the renal cell carcinoma (RCC). Using the proposed approach we provide a metabolic description of the known morphological RCC subtypes and suggest a possible existence of the poor-prognosis cluster of patients, which are commonly characterized by the low activity of the drug transporting enzymes crucial in the chemotherapy. This discovery suits and extends the already known poor-prognosis characteristics of RCC. Finally, the goal of this work is also to point out the problem that arises from the integration of high-throughput data with the inherently nonuniform, manually curated low-throughput data. In such cases, the over-represented information may potentially overshadow the non-trivial discoveries.

Handling Complexity in Animal and Plant Science Research-From Single to Functional Traits: Are We There Yet?

The current knowledge of the main factors governing livestock, crop and plant quality as well as yield in different species is incomplete. For example, this can be evidenced by the persistence of benchmark crop varieties for many decades in spite of the gains achieved over the same period. In recent years, it has been demonstrated that molecular breeding based on DNA markers has led to advances in breeding (animal and crops). However, these advances are not in the way that it was anticipated initially by the researcher in the field. According to several scientists, one of the main reasons for this was related to the evidence that complex target traits such as grain yield, composition or nutritional quality depend on multiple factors in addition to genetics. Therefore, some questions need to be asked: are the current approaches in molecular genetics the most appropriate to deal with complex traits such as yield or quality? Are the current tools for phenotyping complex traits enough to differentiate among genotypes? Do we need to change the way that data is collected and analysed?

Heading in the right direction: thermodynamics-based network analysis and pathway engineering

Thermodynamics-based network analysis through the introduction of thermodynamic constraints in metabolic models allows a deeper analysis of metabolism and guides pathway engineering. The number and the areas of applications of thermodynamics-based network analysis methods have been increasing in the last ten years. We review recent applications of these methods and we identify the areas that such analysis can contribute significantly, and the needs for future developments. We find that organisms with multiple compartments and extremophiles present challenges for modeling and thermodynamics-based flux analysis. The evolution of current and new methods must also address the issues of the multiple alternatives in flux directionalities and the uncertainties and partial information from analytical methods.

Elementary flux modes in a nutshell: properties, calculation and applications

Elementary flux mode (EFM) analysis allows the unbiased decomposition of a metabolic network into minimal functional units, making it a powerful tool for metabolic engineering. While the use of EFM analysis (EFMA) is still limited by the size of the models it can handle, EFMA has been successfully applied to solve real-world metabolic engineering problems. Here we provide a user-oriented introduction to EFMA, provide examples of recent applications, analyze current research strategies to overcome the computational restrictions and give an overview over current approaches, which aim to identify and calculate only biologically relevant EFMs.

Construction and analysis of the model of energy metabolism in E. coli

Genome-scale models of metabolism have only been analyzed with the constraint-based modelling philosophy and there have been several genome-scale gene-protein-reaction models. But research on the modelling for energy metabolism of organisms just began in recent years and research on metabolic weighted complex network are rare in literature. We have made three research based on the complete model of E. coli’s energy metabolism. We first constructed a metabolic weighted network using the rates of free energy consumption within metabolic reactions as the weights. We then analyzed some structural characters of the metabolic weighted network that we constructed. We found that the distribution of the weight values was uneven, that most of the weight values were zero while reactions with abstract large weight values were rare and that the relationship between w (weight values) and v (flux values) was not of linear correlation. At last, we have done some research on the equilibrium of free energy for the energy metabolism system of E. coli. We found that (free energy rate input from the environment) can meet the demand of (free energy rate dissipated by chemical process) and that chemical process plays a great role in the dissipation of free energy in cells. By these research and to a certain extend, we can understand more about the energy metabolism of E. coli.