Model-guided identification of a therapeutic strategy to reduce hyperammonemia in liver diseases.

BACKGROUND & AIMS: Recently, spatial-temporal/metabolic mathematical models have been established that allow the simulation of metabolic processes in tissues. We applied these models to decipher ammonia detoxification mechanisms in the liver. METHODS: An integrated metabolic-spatial-temporal model was used to generate hypotheses of ammonia metabolism. Predicted mechanisms were validated using time-resolved analyses of nitrogen metabolism, activity analyses, immunostaining and gene expression after induction of liver damage in mice. Moreover, blood from the portal vein, liver vein and mixed venous blood was analyzed in a time dependent manner. RESULTS: Modeling revealed an underestimation of ammonia consumption after liver damage when only the currently established mechanisms of ammonia detoxification were simulated. By iterative cycles of modeling and experiments, the reductive amidation of alpha-ketoglutarate (α-KG) via glutamate dehydrogenase (GDH) was identified as the lacking component. GDH is released from damaged hepatocytes into the blood where it consumes ammonia to generate glutamate, thereby providing systemic protection against hyperammonemia. This mechanism was exploited therapeutically in a mouse model of hyperammonemia by injecting GDH together with optimized doses of cofactors. Intravenous injection of GDH (720 U/kg), α-KG (280 mg/kg) and NADPH (180 mg/kg) reduced the elevated blood ammonia concentrations (>200 µM) to levels close to normal within only 15 min. CONCLUSION: If successfully translated to patients the GDH-based therapy might provide a less aggressive therapeutic alternative for patients with severe hyperammonemia.

Integrated metabolic spatial-temporal model for the prediction of ammonia detoxification during liver damage and regeneration.

UNLABELLED: The impairment of hepatic metabolism due to liver injury has high systemic relevance. However, it is difficult to calculate the impairment of metabolic capacity from a specific pattern of liver damage with conventional techniques. We established an integrated metabolic spatial-temporal model (IM) using hepatic ammonia detoxification as a paradigm. First, a metabolic model (MM) based on mass balancing and mouse liver perfusion data was established to describe ammonia detoxification and its zonation. Next, the MM was combined with a spatial-temporal model simulating liver tissue damage and regeneration after CCl4 intoxication. The resulting IM simulated and visualized whether, where, and to what extent liver damage compromised ammonia detoxification. It allowed us to enter the extent and spatial patterns of liver damage and then calculate the outflow concentrations of ammonia, glutamine, and urea in the hepatic vein. The model was validated through comparisons with (1) published data for isolated, perfused livers with and without CCl4 intoxication and (2) a set of in vivo experiments. Using the experimentally determined portal concentrations of ammonia, the model adequately predicted metabolite concentrations over time in the hepatic vein during toxin-induced liver damage and regeneration in rodents. Further simulations, especially in combination with a simplified model of blood circulation with three ammonia-detoxifying compartments, indicated a yet unidentified process of ammonia consumption during liver regeneration and revealed unexpected concomitant changes in amino acid metabolism in the liver and at extrahepatic sites. CONCLUSION: The IM of hepatic ammonia detoxification considerably improves our understanding of the metabolic impact of liver disease and highlights the importance of integrated modeling approaches on the way toward virtual organisms.

Optimality in the zonation of ammonia detoxification in rodent liver.

The rodent liver eliminates toxic ammonia. In mammals, three enzymes (or enzyme systems) are involved in this process: glutaminase, glutamine synthetase and the urea cycle enzymes, represented by carbamoyl phosphate synthetase. The distribution of these enzymes for optimal ammonia detoxification was determined by numerical optimization. This in silico approach predicted that the enzymes have to be zonated in order to achieve maximal removal of toxic ammonia and minimal changes in glutamine concentration. Using 13 compartments, representing hepatocytes, the following predictions were generated: glutamine synthetase is active only within a narrow pericentral zone. Glutaminase and carbamoyl phosphate synthetase are located in the periportal zone in a non-homogeneous distribution. This correlates well with the paradoxical observation that in a first step glutamine-bound ammonia is released (by glutaminase) although one of the functions of the liver is detoxification by ammonia fixation. The in silico approach correctly predicted the in vivo enzyme distributions also for non-physiological conditions (e.g. starvation) and during regeneration after tetrachloromethane (CCl4) intoxication. Metabolite concentrations of glutamine, ammonia and urea in each compartment, representing individual hepatocytes, were predicted. Finally, a sensitivity analysis showed a striking robustness of the results. These bioinformatics predictions were validated experimentally by immunohistochemistry and are supported by the literature. In summary, optimization approaches like the one applied can provide valuable explanations and high-quality predictions for in vivo enzyme and metabolite distributions in tissues and can reveal unknown metabolic functions.

Data- and knowledge-based modeling of gene regulatory networks: an update.

Gene regulatory network inference is a systems biology approach which predicts interactions between genes with the help of high-throughput data. In this review, we present current and updated network inference methods focusing on novel techniques for data acquisition, network inference assessment, network inference for interacting species and the integration of prior knowledge. After the advance of Next-Generation-Sequencing of cDNAs derived from RNA samples (RNA-Seq) we discuss in detail its application to network inference. Furthermore, we present progress for large-scale or even full-genomic network inference as well as for small-scale condensed network inference and review advances in the evaluation of network inference methods by crowdsourcing. Finally, we reflect the current availability of data and prior knowledge sources and give an outlook for the inference of gene regulatory networks that reflect interacting species, in particular pathogen-host interactions.

Computational prediction of molecular pathogen-host interactions based on dual transcriptome data.

Inference of inter-species gene regulatory networks based on gene expression data is an important computational method to predict pathogen-host interactions (PHIs). Both the experimental setup and the nature of PHIs exhibit certain characteristics. First, besides an environmental change, the battle between pathogen and host leads to a constantly changing environment and thus complex gene expression patterns. Second, there might be a delay until one of the organisms reacts. Third, toward later time points only one organism may survive leading to missing gene expression data of the other organism. Here, we account for PHI characteristics by extending NetGenerator, a network inference tool that predicts gene regulatory networks from gene expression time series data. We tested multiple modeling scenarios regarding the stimuli functions of the interaction network based on a benchmark example. We show that modeling perturbation of a PHI network by multiple stimuli better represents the underlying biological phenomena. Furthermore, we utilized the benchmark example to test the influence of missing data points on the inference performance. Our results suggest that PHI network inference with missing data is possible, but we recommend to provide complete time series data. Finally, we extended the NetGenerator tool to incorporate gene- and time point specific variances, because complex PHIs may lead to high variance in expression data. Sample variances are directly considered in the objective function of NetGenerator and indirectly by testing the robustness of interactions based on variance dependent disturbance of gene expression values. We evaluated the method of variance incorporation on dual RNA sequencing (RNA-Seq) data of Mus musculus dendritic cells incubated with Candida albicans and proofed our method by predicting previously verified PHIs as robust interactions.

Basic regulatory principles of Escherichia coli's electron transport chain for varying oxygen conditions.

For adaptation between anaerobic, micro-aerobic and aerobic conditions Escherichia coli's metabolism and in particular its electron transport chain (ETC) is highly regulated. Although it is known that the global transcriptional regulators FNR and ArcA are involved in oxygen response it is unclear how they interplay in the regulation of ETC enzymes under micro-aerobic chemostat conditions. Also, there are diverse results which and how quinones (oxidised/reduced, ubiquinone/other quinones) are controlling the ArcBA two-component system. In the following a mathematical model of the E. coli ETC linked to basic modules for substrate uptake, fermentation product excretion and biomass formation is introduced. The kinetic modelling focusses on regulatory principles of the ETC for varying oxygen conditions in glucose-limited continuous cultures. The model is based on the balance of electron donation (glucose) and acceptance (oxygen or other acceptors). Also, it is able to account for different chemostat conditions due to changed substrate concentrations and dilution rates. The parameter identification process is divided into an estimation and a validation step based on previously published and new experimental data. The model shows that experimentally observed, qualitatively different behaviour of the ubiquinone redox state and the ArcA activity profile in the micro-aerobic range for different experimental conditions can emerge from a single network structure. The network structure features a strong feed-forward effect from the FNR regulatory system to the ArcBA regulatory system via a common control of the dehydrogenases of the ETC. The model supports the hypothesis that ubiquinone but not ubiquinol plays a key role in determining the activity of ArcBA in a glucose-limited chemostat at micro-aerobic conditions.

Inference of dynamical gene-regulatory networks based on time-resolved multi-stimuli multi-experiment data applying NetGenerator V2.0.

BACKGROUND: Inference of gene-regulatory networks (GRNs) is important for understanding behaviour and potential treatment of biological systems. Knowledge about GRNs gained from transcriptome analysis can be increased by multiple experiments and/or multiple stimuli. Since GRNs are complex and dynamical, appropriate methods and algorithms are needed for constructing models describing these dynamics. Algorithms based on heuristic approaches reduce the effort in parameter identification and computation time. RESULTS: The NetGenerator V2.0 algorithm, a heuristic for network inference, is proposed and described. It automatically generates a system of differential equations modelling structure and dynamics of the network based on time-resolved gene expression data. In contrast to a previous version, the inference considers multi-stimuli multi-experiment data and contains different methods for integrating prior knowledge. The resulting significant changes in the algorithmic procedures are explained in detail. NetGenerator is applied to relevant benchmark examples evaluating the inference for data from experiments with different stimuli. Also, the underlying GRN of chondrogenic differentiation, a real-world multi-stimulus problem, is inferred and analysed. CONCLUSIONS: NetGenerator is able to determine the structure and parameters of GRNs and their dynamics. The new features of the algorithm extend the range of possible experimental set-ups, results and biological interpretations. Based upon benchmarks, the algorithm provides good results in terms of specificity, sensitivity, efficiency and model fit.