Current state and challenges for dynamic metabolic modeling

Emerging methods for genome-scale metabolic modeling of microbial communities

Genome-scale metabolic models (GEMs) are consolidating as platforms for studying mixed microbial populations, by combining biological data and knowledge with mathematical rigor. However, deploying these models to answer research questions can be challenging due to the increasing number of available computational tools, the lack of universal standards, and their inherent limitations. Here, we present a comprehensive overview of foundational concepts for building and evaluating genome-scale models of microbial communities. We then compare tools in terms of requirements, capabilities, and applications. Next, we highlight the current pitfalls and open challenges to consider when adopting existing tools and developing new ones. Our compendium can be relevant for the expanding community of modelers, both at the entry and experienced levels.

Review on computer-assisted biosynthetic capacities elucidation to assess metabolic interactions and communication within microbial communities

Microbial communities thrive through interactions and communication, which are challenging to study as most microorganisms are not cultivable. To address this challenge, researchers focus on the extracellular space where communication events occur. Exometabolomics and interactome analysis provide insights into the molecules involved in communication and the dynamics of their interactions. Advances in sequencing technologies and computational methods enable the reconstruction of taxonomic and functional profiles of microbial communities using high-throughput multi-omics data. Network-based approaches, including community flux balance analysis, aim to model molecular interactions within and between communities. Despite these advances, challenges remain in computer-assisted biosynthetic capacities elucidation, requiring continued innovation and collaboration among diverse scientists. This review provides insights into the current state and future directions of computer-assisted biosynthetic capacities elucidation in studying microbial communities.

Control Theory and Systems Biology: Potential Applications in Neurodegeneration and Search for Therapeutic Targets

Control theory, a well-established discipline in engineering and mathematics, has found novel applications in systems biology. This interdisciplinary approach leverages the principles of feedback control and regulation to gain insights into the complex dynamics of cellular and molecular networks underlying chronic diseases, including neurodegeneration. By modeling and analyzing these intricate systems, control theory provides a framework to understand the pathophysiology and identify potential therapeutic targets. Therefore, this review examines the most widely used control methods in conjunction with genomic-scale metabolic models in the steady state of the multi-omics type. According to our research, this approach involves integrating experimental data, mathematical modeling, and computational analyses to simulate and control complex biological systems. In this review, we find that the most significant application of this methodology is associated with cancer, leaving a lack of knowledge in neurodegenerative models. However, this methodology, mainly associated with the Minimal Dominant Set (MDS), has provided a starting point for identifying therapeutic targets for drug development and personalized treatment strategies, paving the way for more effective therapies.

Monitoring and modelling the glutamine metabolic pathway: a review and future perspectives

BACKGROUND

Analysis of the glutamine metabolic pathway has taken a special place in metabolomics research in recent years, given its important role in cell biosynthesis and bioenergetics across several disorders, especially in cancer cell survival. The science of metabolomics addresses the intricate intracellular metabolic network by exploring and understanding how cells function and respond to external or internal perturbations to identify potential therapeutic targets. However, despite recent advances in metabolomics, monitoring the kinetics of a metabolic pathway in a living cell in situ, real-time and holistically remains a significant challenge.

AIM

This review paper explores the range of analytical approaches for monitoring metabolic pathways, as well as physicochemical modeling techniques, with a focus on glutamine metabolism. We discuss the advantages and disadvantages of each method and explore the potential of label-free Raman microspectroscopy, in conjunction with kinetic modeling, to enable real-time and in situ monitoring of the cellular kinetics of the glutamine metabolic pathway.

KEY SCIENTIFIC CONCEPTS

Given its important role in cell metabolism, the ability to monitor and model the glutamine metabolic pathways are highlighted. Novel, label free approaches have the potential to revolutionise metabolic biosensing, laying the foundation for a new paradigm in metabolomics research and addressing the challenges in monitoring metabolic pathways in living cells.

Exogenous metabolite feeding on altering antibiotic susceptibility in Gram-negative bacteria through metabolic modulation: a review

BACKGROUND

The rise of antimicrobial resistance at an alarming rate is outpacing the development of new antibiotics. The worrisome trends of multidrug-resistant Gram-negative bacteria have enormously diminished existing antibiotic activity. Antibiotic treatments may inhibit bacterial growth or lead to induce bacterial cell death through disruption of bacterial metabolism directly or indirectly. In light of this, it is imperative to have a thorough understanding of the relationship of bacterial metabolism with antimicrobial activity and leverage the underlying principle towards development of novel and effective antimicrobial therapies.

OBJECTIVE

Herein, we explore studies on metabolic analyses of Gram-negative pathogens upon antibiotic treatment. Metabolomic studies revealed that antibiotic therapy caused changes of metabolites abundance and perturbed the bacterial metabolism. Following this line of thought, addition of exogenous metabolite has been employed in in vitro, in vivo and in silico studies to activate the bacterial metabolism and thus potentiate the antibiotic activity.

KEY SCIENTIFIC CONCEPTS OF REVIEW

Exogenous metabolites were discovered to cause metabolic modulation through activation of central carbon metabolism and cellular respiration, stimulation of proton motive force, increase of membrane potential, improvement of host immune protection, alteration of gut microbiome, and eventually facilitating antibiotic killing. The use of metabolites as antimicrobial adjuvants may be a promising approach in the fight against multidrug-resistant pathogens.

Kinetic Modeling of Saccharomyces cerevisiae Central Carbon Metabolism: Achievements, Limitations, and Opportunities

Central carbon metabolism comprises the metabolic pathways in the cell that process nutrients into energy, building blocks and byproducts. To unravel the regulation of this network upon glucose perturbation, several metabolic models have been developed for the microorganism Saccharomyces cerevisiae. These dynamic representations have focused on glycolysis and answered multiple research questions, but no commonly applicable model has been presented. This review systematically evaluates the literature to describe the current advances, limitations, and opportunities. Different kinetic models have unraveled key kinetic glycolytic mechanisms. Nevertheless, some uncertainties regarding model topology and parameter values still limit the application to specific cases. Progressive improvements in experimental measurement technologies as well as advances in computational tools create new opportunities to further extend the model scale. Notably, models need to be made more complex to consider the multiple layers of glycolytic regulation and external physiological variables regulating the bioprocess, opening new possibilities for extrapolation and validation. Finally, the onset of new data representative of individual cells will cause these models to evolve from depicting an average cell in an industrial fermenter, to characterizing the heterogeneity of the population, opening new and unseen possibilities for industrial fermentation improvement.

Glycolysis/gluconeogenesis specialization in microbes is driven by biochemical constraints of flux sensing

Central carbon metabolism is highly conserved across microbial species, but can catalyze very different pathways depending on the organism and their ecological niche. Here, we study the dynamic reorganization of central metabolism after switches between the two major opposing pathway configurations of central carbon metabolism, glycolysis, and gluconeogenesis in Escherichia coli, Pseudomonas aeruginosa, and Pseudomonas putida. We combined growth dynamics and dynamic changes in intracellular metabolite levels with a coarse-grained model that integrates fluxes, regulation, protein synthesis, and growth and uncovered fundamental limitations of the regulatory network: After nutrient shifts, metabolite concentrations collapse to their equilibrium, rendering the cell unable to sense which direction the flux is supposed to flow through the metabolic network. The cell can partially alleviate this by picking a preferred direction of regulation at the expense of increasing lag times in the opposite direction. Moreover, decreasing both lag times simultaneously comes at the cost of reduced growth rate or higher futile cycling between metabolic enzymes. These three trade-offs can explain why microorganisms specialize for either glycolytic or gluconeogenic substrates and can help elucidate the complex growth patterns exhibited by different microbial species.

Simultaneous Quantification of the Concentration and Carbon Isotopologue Distribution of Polar Metabolites in a Single Analysis by Gas Chromatography and Mass Spectrometry

¹³C-isotope tracing is a frequently employed approach to study metabolic pathway activity. When combined with the subsequent quantification of absolute metabolite concentrations, this enables detailed characterization of the metabolome in biological specimens and facilitates computational time-resolved flux quantification. Classically, a ¹³C-isotopically labeled sample is required to quantify ¹³C-isotope enrichments and a second unlabeled sample for the quantification of metabolite concentrations. The rationale for a second unlabeled sample is that the current methods for metabolite quantification rely mostly on isotope dilution mass spectrometry (IDMS) and thus isotopically labeled internal standards are added to the unlabeled sample. This excludes the absolute quantification of metabolite concentrations in ¹³C-isotopically labeled samples. To address this issue, we have developed and validated a new strategy using an unlabeled internal standard to simultaneously quantify metabolite concentrations and ¹³C-isotope enrichments in a single ¹³C-labeled sample based on gas chromatography-mass spectrometry (GC/MS). The method was optimized for amino acids and citric acid cycle intermediates and was shown to have high analytical precision and accuracy. Metabolite concentrations could be quantified in small tissue samples (≥20 mg). Also, we applied the method on ¹³C-isotopically labeled mammalian cells treated with and without a metabolic inhibitor. We proved that we can quantify absolute metabolite concentrations and ¹³C-isotope enrichments in a single ¹³C-isotopically labeled sample.

Computational/in silico methods in drug target and lead prediction

Drug-like compounds are most of the time denied approval and use owing to the unexpected clinical side effects and cross-reactivity observed during clinical trials. These unexpected outcomes resulting in significant increase in attrition rate centralizes on the selected drug targets. These targets may be disease candidate proteins or genes, biological pathways, disease-associated microRNAs, disease-related biomarkers, abnormal molecular phenotypes, crucial nodes of biological network or molecular functions. This is generally linked to several factors, including incomplete knowledge on the drug targets and unpredicted pharmacokinetic expressions upon target interaction or off-target effects. A method used to identify targets, especially for polygenic diseases, is essential and constitutes a major bottleneck in drug development with the fundamental stage being the identification and validation of drug targets of interest for further downstream processes. Thus, various computational methods have been developed to complement experimental approaches in drug discovery. Here, we present an overview of various computational methods and tools applied in predicting or validating drug targets and drug-like molecules. We provide an overview on their advantages and compare these methods to identify effective methods which likely lead to optimal results. We also explore major sources of drug failure considering the challenges and opportunities involved. This review might guide researchers on selecting the most efficient approach or technique during the computational drug discovery process.

Application of Systems Engineering Principles and Techniques in Biological Big Data Analytics: A Review

In the past few decades, we have witnessed tremendous advancements in biology, life sciences and healthcare. These advancements are due in no small part to the big data made available by various high-throughput technologies, the ever-advancing computing power, and the algorithmic advancements in machine learning. Specifically, big data analytics such as statistical and machine learning has become an essential tool in these rapidly developing fields. As a result, the subject has drawn increased attention and many review papers have been published in just the past few years on the subject. Different from all existing reviews, this work focuses on the application of systems, engineering principles and techniques in addressing some of the common challenges in big data analytics for biological, biomedical and healthcare applications. Specifically, this review focuses on the following three key areas in biological big data analytics where systems engineering principles and techniques have been playing important roles: the principle of parsimony in addressing overfitting, the dynamic analysis of biological data, and the role of domain knowledge in biological data analytics.

Metabolic Modelling as a Framework for Metabolomics Data Integration and Analysis

Metabolic networks are regulated to ensure the dynamic adaptation of biochemical reaction fluxes to maintain cell homeostasis and optimal metabolic fitness in response to endogenous and exogenous perturbations. To this end, metabolism is tightly controlled by dynamic and intricate regulatory mechanisms involving allostery, enzyme abundance and post-translational modifications. The study of the molecular entities involved in these complex mechanisms has been boosted by the advent of high-throughput technologies. The so-called omics enable the quantification of the different molecular entities at different system layers, connecting the genotype with the phenotype. Therefore, the study of the overall behavior of a metabolic network and the omics data integration and analysis must be approached from a holistic perspective. Due to the close relationship between metabolism and cellular phenotype, metabolic modelling has emerged as a valuable tool to decipher the underlying mechanisms governing cell phenotype. Constraint-based modelling and kinetic modelling are among the most widely used methods to study cell metabolism at different scales, ranging from cells to tissues and organisms. These approaches enable integrating metabolomic data, among others, to enhance model predictive capabilities. In this review, we describe the current state of the art in metabolic modelling and discuss future perspectives and current challenges in the field.

Strategies for the Biosynthesis of Pharmaceuticals and Nutraceuticals in Microbes from Renewable Feedstock

BACKGROUNDS

Abundant and renewable biomaterials serve as ideal substrates for the sustainable production of various chemicals, including natural products (e.g., pharmaceuticals and nutraceuticals). For decades, researchers have been focusing on how to engineer microorganisms and developing effective fermentation processes to overproduce these molecules from biomaterials. Despite many laboratory achievements, it remains a challenge to transform some of these into successful industrial applications.

RESULTS

Here, we review recent progress in strategies and applications in metabolic engineering for the production of natural products. Modular engineering methods, such as a multidimensional heuristic process markedly improve efficiencies in the optimization of long and complex biosynthetic pathways. Dynamic pathway regulation realizes autonomous adjustment and can redirect metabolic carbon fluxes to avoid the accumulation of toxic intermediate metabolites. Microbial co-cultivation bolsters the identification and overproduction of natural products by introducing competition or cooperation of different species. Efflux engineering is applied to reduce product toxicity or to overcome storage limitation and thus improves product titers and productivities.

CONCLUSION

Without dispute, many of the innovative methods and strategies developed are gradually catalyzing this transformation from the laboratory into the industry in the biosynthesis of natural products. Sometimes, it is necessary to combine two or more strategies to acquire additive or synergistic benefits. As such, we foresee a bright future of the biosynthesis of pharmaceuticals and nutraceuticals in microbes from renewable biomaterials.

Hyperpolarized product selective saturating-excitations for determination of changes in metabolic reaction rates in real-time

Investigation of hyperpolarized substrate metabolism has been showing utility in real-time determination of in-cell and in vivo enzymatic activities. Intracellular reaction rates may vary during the course of a measurement, even on the very short time scales of visibility on hyperpolarized MR, due to many factors such as the availability of the substrate and co-factors in the intracellular space. Despite this potential variation, the kinetic analysis of hyperpolarized signals typically assumes that the same rate constant (and in many cases, the same rate) applies throughout the course of the reaction as observed via the build-up and decay of the hyperpolarized signals. We demonstrate here an acquisition approach that can null the need for such an assumption and enable the detection of instantaneous changes in the rate of the reaction during an ex vivo hyperpolarized investigation, (i.e. in the course of the decay of one hyperpolarized substrate dose administered to a viable tissue sample ex vivo). This approach utilizes hyperpolarized product selective saturating-excitation pulses. Similar pulses have been previously utilized in vivo for spectroscopic imaging. However, we show here favorable consequences to kinetic rate determinations in the preparations used. We implement this acquisition strategy for studies on perfused tissue slices and develop a theory that explains why this particular approach enables the determination of changes in enzymatic rates that are monitored via the chemical conversions of hyperpolarized substrates. Real-time changes in intracellular reaction rates are demonstrated in perfused brain, liver, and xenograft breast cancer tissue slices and provide another potential differentiation parameter for tissue characterization.

Flux prediction using artificial neural network (ANN) for the upper part of glycolysis

The selection of optimal enzyme concentration in multienzyme cascade reactions for the highest product yield in practice is very expensive and time-consuming process. The modelling of biological pathways is a difficult process because of the complexity of the system. The mathematical modelling of the system using an analytical approach depends on the many parameters of enzymes which rely on tedious and expensive experiments. The artificial neural network (ANN) method has been successively applied in different fields of science to perform complex functions. In this study, ANN models were trained to predict the flux for the upper part of glycolysis as inferred by NADH consumption, using four enzyme concentrations i.e., phosphoglucoisomerase, phosphofructokinase, fructose-bisphosphate-aldolase, triose-phosphate-isomerase. Out of three ANN algorithms, the neuralnet package with two activation functions, “logistic” and “tanh” were implemented. The prediction of the flux was very efficient: RMSE and R² were 0.847, 0.93 and 0.804, 0.94 respectively for logistic and tanh functions using a cross validation procedure. This study showed that a systemic approach such as ANN could be used for accurate prediction of the flux through the metabolic pathway. This could help to save a lot of time and costs, particularly from an industrial perspective. The R-code is available at: https://github.com/DSIMB/ANN-Glycolysis-Flux-Prediction.

Approaches to Computational Strain Design in the Multiomics Era

Modern omics analyses are able to effectively characterize the genetic, regulatory, and metabolic phenotypes of engineered microbes, yet designing genetic interventions to achieve a desired phenotype remains challenging. With recent developments in genetic engineering techniques, timelines associated with building and testing strain designs have been greatly reduced, allowing for the first time an efficient closed loop iteration between experiment and analysis. However, the scale and complexity associated with multi-omics datasets complicates manual biological reasoning about the mechanisms driving phenotypic changes. Computational techniques therefore form a critical part of the Design-Build-Test-Learn (DBTL) cycle in metabolic engineering. Traditional statistical approaches can reduce the dimensionality of these datasets and identify common motifs among high-performing strains. While successful in many studies, these methods do not take full advantage of known connections between genes, proteins, and metabolic networks. There is therefore a growing interest in model-aided design, in which modeling frameworks from systems biology are used to integrate experimental data and generate effective and non-intuitive design predictions. In this mini-review, we discuss recent progress and challenges in this field. In particular, we compare methods augmenting flux balance analysis with additional constraints from fluxomic, genomic, and metabolomic datasets and methods employing kinetic representations of individual metabolic reactions, and machine learning. We conclude with a discussion of potential future directions for improving strain design predictions in the omics era and remaining experimental and computational hurdles.

Combined In-Cell NMR and Simulation Approach to Probe Redox-Dependent Pathway Control

Dynamic response of intracellular reaction cascades to changing environments is a hallmark of living systems. As metabolism is complex, mechanistic models have gained popularity for describing the dynamic response of cellular metabolism and for identifying target genes for engineering. At the same time, the detailed tracking of transient metabolism in living cells on the subminute time scale has become amenable using dynamic nuclear polarization-enhanced ¹³C NMR. Here, we suggest an approach combining in-cell NMR spectroscopy with perturbation experiments and modeling to obtain evidence that the bottlenecks of yeast glycolysis depend on intracellular redox state. In pre-steady-state glycolysis, pathway bottlenecks shift from downstream to upstream reactions within a few seconds, consistent with a rapid decline in the NAD⁺/NADH ratio. Simulations using mechanistic models reproduce the experimentally observed response and help identify unforeseen biochemical events. Remaining inaccuracies in the computational models can be identified experimentally. The combined use of rapid injection NMR spectroscopy and in silico simulations provides a promising method for characterizing cellular reactions with increasing mechanistic detail.

Metabolic Labeling of Cultured Mammalian Cells for Stable Isotope-Resolved Metabolomics: Practical Aspects of Tissue Culture and Sample Extraction

Stable isotope-resolved metabolomics (SIRM) methods are used increasingly by cancer researchers to probe metabolic pathways and identify vulnerabilities in cancer cells. Analytical and computational advances are being made constantly, but tissue culture and sample extraction procedures are often variable and not elaborated in the literature. This chapter discusses basic aspects of tissue culture practices as they relate to the use of stable isotope tracers and provides a detailed metabolic labeling and metabolite extraction procedure designed to maximize the amount of information that can be obtained from a single tracer experiment.

Integrating -omics data into genome-scale metabolic network models: principles and challenges

At genome scale, it is not yet possible to devise detailed kinetic models for metabolism because data on the in vivo biochemistry are too sparse. Predictive large-scale models for metabolism most commonly use the constraint-based framework, in which network structures constrain possible metabolic phenotypes at steady state. However, these models commonly leave many possibilities open, making them less predictive than desired. With increasingly available -omics data, it is appealing to increase the predictive power of constraint-based models (CBMs) through data integration. Many corresponding methods have been developed, but data integration is still a challenge and existing methods perform less well than expected. Here, we review main approaches for the integration of different types of -omics data into CBMs focussing on the methods' assumptions and limitations. We argue that key assumptions - often derived from single-enzyme kinetics - do not generally apply in the context of networks, thereby explaining current limitations. Emerging methods bridging CBMs and biochemical kinetics may allow for -omics data integration in a common framework to provide more accurate predictions.

Multidimensional heuristic process for high-yield production of astaxanthin and fragrance molecules in Escherichia coli

Optimization of metabolic pathways consisting of large number of genes is challenging. Multivariate modular methods (MMMs) are currently available solutions, in which reduced regulatory complexities are achieved by grouping multiple genes into modules. However, these methods work well for balancing the inter-modules but not intra-modules. In addition, application of MMMs to the 15-step heterologous route of astaxanthin biosynthesis has met with limited success. Here, we expand the solution space of MMMs and develop a multidimensional heuristic process (MHP). MHP can simultaneously balance different modules by varying promoter strength and coordinating intra-module activities by using ribosome binding sites (RBSs) and enzyme variants. Consequently, MHP increases enantiopure 3S,3′S-astaxanthin production to 184 mg l⁻¹ day⁻¹ or 320 mg l⁻¹. Similarly, MHP improves the yields of nerolidol and linalool. MHP may be useful for optimizing other complex biochemical pathways.

Achieving high titer yield and productivity of target chemicals in industrial organism depends on multidimensional pathway optimization. Here, the authors use a refined modular method called multidimensional heuristic process to improve production of astaxanthin, nerolidol and linalool in E. coli.

Flux analysis of inborn errors of metabolism

Patients with an inborn error of metabolism (IEM) are deficient of an enzyme involved in metabolism, and as a consequence metabolism reprograms itself to reach a new steady state. This new steady state underlies the clinical phenotype associated with the deficiency. Hence, we need to know the flux of metabolites through the different metabolic pathways in this new steady state of the reprogrammed metabolism. Stable isotope technology is best suited to study this. In this review the progress made in characterizing the altered metabolism will be presented. Studies done in patients to estimate the residual flux through the metabolic pathway affected by enzyme deficiencies will be discussed. After this, studies done in model systems will be reviewed. The focus will be on glycogen storage disease type I, medium-chain acyl-CoA dehydrogenase deficiency, propionic and methylmalonic aciduria, urea cycle defects, phenylketonuria, and combined D,L-2-hydroxyglutaric aciduria. Finally, new developments are discussed, which allow the tracing of metabolic reprogramming in IEM on a genome-wide scale. In conclusion, the outlook for flux analysis of metabolic derangement in IEMs looks promising.

Metabolic engineering of Pichia pastoris

Besides its use for efficient production of recombinant proteins the methylotrophic yeast Pichia pastoris (syn. Komagataella spp.) has been increasingly employed as a platform to produce metabolites of varying origin. We summarize here the impressive methodological developments of the last years to model and analyze the metabolism of P. pastoris, and to engineer its genome and metabolic pathways. Efficient methods to insert, modify or delete genes via homologous recombination and CRISPR/Cas9, supported by modular cloning techniques, have been reported. An outstanding early example of metabolic engineering in P. pastoris was the humanization of protein glycosylation. More recently the cell metabolism was engineered also to enhance the productivity of heterologous proteins. The last few years have seen an increased number of metabolic pathway design and engineering in P. pastoris, mainly towards the production of complex (secondary) metabolites. In this review, we discuss the potential role of P. pastoris as a platform for metabolic engineering, its strengths, and major requirements for future developments of chassis strains based on synthetic biology principles.

Corynebacterium glutamicum Chassis C1*: Building and Testing a Novel Platform Host for Synthetic Biology and Industrial Biotechnology

Targeted top-down strategies for genome reduction are considered to have a high potential for providing robust basic strains for synthetic biology and industrial biotechnology. Recently, we created a library of 26 genome-reduced strains of Corynebacterium glutamicum carrying broad deletions in single gene clusters and showing wild-type-like biological fitness. Here, we proceeded with combinatorial deletions of these irrelevant gene clusters in two parallel orders, and the resulting library of 28 strains was characterized under various environmental conditions. The final chassis strain C1* carries a genome reduction of 13.4% (412 deleted genes) and shows wild-type-like growth behavior in defined medium with d-glucose as carbon and energy source. Moreover, C1* proves to be robust against several stresses (including oxygen limitation) and shows long-term growth stability under defined and complex medium conditions. In addition to providing a novel prokaryotic chassis strain, our results comprise a large strain library and a revised genome annotation list, which will be valuable sources for future systemic studies of C. glutamicum.

Stage-specific metabolic features of differentiating neurons: Implications for toxicant sensitivity

Developmental neurotoxicity (DNT) may be induced when chemicals disturb a key neurodevelopmental process, and many tests focus on this type of toxicity. Alternatively, DNT may occur when chemicals are cytotoxic only during a specific neurodevelopmental stage. The toxicant sensitivity is affected by the expression of toxicant targets and by resilience factors. Although cellular metabolism plays an important role, little is known how it changes during human neurogenesis, and how potential alterations affect toxicant sensitivity of mature vs. immature neurons. We used immature (d0) and mature (d6) LUHMES cells (dopaminergic human neurons) to provide initial answers to these questions. Transcriptome profiling and characterization of energy metabolism suggested a switch from predominantly glycolytic energy generation to a more pronounced contribution of the tricarboxylic acid cycle (TCA) during neuronal maturation. Therefore, we used pulsed stable isotope-resolved metabolomics (pSIRM) to determine intracellular metabolite pool sizes (concentrations), and isotopically non-stationary ¹³C-metabolic flux analysis (INST ¹³C-MFA) to calculate metabolic fluxes. We found that d0 cells mainly use glutamine to fuel the TCA. Furthermore, they rely on extracellular pyruvate to allow continuous growth. This metabolic situation does not allow for mitochondrial or glycolytic spare capacity, i.e. the ability to adapt energy generation to altered needs. Accordingly, neuronal precursor cells displayed a higher sensitivity to several mitochondrial toxicants than mature neurons differentiated from them. In summary, this study shows that precursor cells lose their glutamine dependency during differentiation while they gain flexibility of energy generation and thereby increase their resistance to low concentrations of mitochondrial toxicants.

Regulation-Structured Dynamic Metabolic Model Provides a Potential Mechanism for Delayed Enzyme Response in Denitrification Process

In a recent study of denitrification dynamics in hyporheic zone sediments, we observed a significant time lag (up to several days) in enzymatic response to the changes in substrate concentration. To explore an underlying mechanism and understand the interactive dynamics between enzymes and nutrients, we developed a trait-based model that associates a community's traits with functional enzymes, instead of typically used species guilds (or functional guilds). This enzyme-based formulation allows to collectively describe biogeochemical functions of microbial communities without directly parameterizing the dynamics of species guilds, therefore being scalable to complex communities. As a key component of modeling, we accounted for microbial regulation occurring through transcriptional and translational processes, the dynamics of which was parameterized based on the temporal profiles of enzyme concentrations measured using a new signature peptide-based method. The simulation results using the resulting model showed several days of a time lag in enzymatic responses as observed in experiments. Further, the model showed that the delayed enzymatic reactions could be primarily controlled by transcriptional responses and that the dynamics of transcripts and enzymes are closely correlated. The developed model can serve as a useful tool for predicting biogeochemical processes in natural environments, either independently or through integration with hydrologic flow simulators.

Formulation, construction and analysis of kinetic models of metabolism: A review of modelling frameworks

Kinetic models are critical to predict the dynamic behaviour of metabolic networks. Mechanistic kinetic models for large networks remain uncommon due to the difficulty of fitting their parameters. Recent modelling frameworks promise new ways to overcome this obstacle while retaining predictive capabilities. In this review, we present an overview of the relevant mathematical frameworks for kinetic formulation, construction and analysis. Starting with kinetic formalisms, we next review statistical methods for parameter inference, as well as recent computational frameworks applied to the construction and analysis of kinetic models. Finally, we discuss opportunities and limitations hindering the development of larger kinetic reconstructions.

Genome-scale modeling of yeast: chronology, applications and critical perspectives

Over the last 15 years, several genome-scale metabolic models (GSMMs) were developed for different yeast species, aiding both the elucidation of new biological processes and the shift toward a bio-based economy, through the design of in silico inspired cell factories. Here, an historical perspective of the GSMMs built over time for several yeast species is presented and the main inheritance patterns among the metabolic reconstructions are highlighted. We additionally provide a critical perspective on the overall genome-scale modeling procedure, underlining incomplete model validation and evaluation approaches and the quest for the integration of regulatory and kinetic information into yeast GSMMs. A summary of experimentally validated model-based metabolic engineering applications of yeast species is further emphasized, while the main challenges and future perspectives for the field are finally addressed.

Parameter identifiability analysis and visualization in large-scale kinetic models of biosystems

Background

Kinetic models of biochemical systems usually consist of ordinary differential equations that have many unknown parameters. Some of these parameters are often practically unidentifiable, that is, their values cannot be uniquely determined from the available data. Possible causes are lack of influence on the measured outputs, interdependence among parameters, and poor data quality. Uncorrelated parameters can be seen as the key tuning knobs of a predictive model. Therefore, before attempting to perform parameter estimation (model calibration) it is important to characterize the subset(s) of identifiable parameters and their interplay. Once this is achieved, it is still necessary to perform parameter estimation, which poses additional challenges.

Methods

We present a methodology that (i) detects high-order relationships among parameters, and (ii) visualizes the results to facilitate further analysis. We use a collinearity index to quantify the correlation between parameters in a group in a computationally efficient way. Then we apply integer optimization to find the largest groups of uncorrelated parameters. We also use the collinearity index to identify small groups of highly correlated parameters. The results files can be visualized using Cytoscape, showing the identifiable and non-identifiable groups of parameters together with the model structure in the same graph.

Results

Our contributions alleviate the difficulties that appear at different stages of the identifiability analysis and parameter estimation process. We show how to combine global optimization and regularization techniques for calibrating medium and large scale biological models with moderate computation times. Then we evaluate the practical identifiability of the estimated parameters using the proposed methodology. The identifiability analysis techniques are implemented as a MATLAB toolbox called VisId, which is freely available as open source from GitHub (https://github.com/gabora/visid).

Conclusions

Our approach is geared towards scalability. It enables the practical identifiability analysis of dynamic models of large size, and accelerates their calibration. The visualization tool allows modellers to detect parts that are problematic and need refinement or reformulation, and provides experimentalists with information that can be helpful in the design of new experiments.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-017-0428-y) contains supplementary material, which is available to authorized users.