Integrating Genome-Scale Metabolic Models with Patient Plasma Metabolome to Study Endothelial Metabolism In Situ.

Patient blood samples are invaluable in clinical omics databases, yet current methodologies often fail to fully uncover the molecular mechanisms driving patient pathology. While genome-scale metabolic models (GEMs) show promise in systems medicine by integrating various omics data, having only exometabolomic data remains a limiting factor. To address this gap, we introduce a comprehensive pipeline integrating GEMs with patient plasma metabolome. This pipeline constructs case-specific GEMs using literature-based and patient-specific metabolomic data. Novel computational methods, including adaptive sampling and an in-house developed algorithm for the rational exploration of the sampled space of solutions, enhance integration accuracy while improving computational performance. Model characterization involves task analysis in combination with clustering methods to identify critical cellular functions. The new pipeline was applied to a cohort of trauma patients to investigate shock-induced endotheliopathy using patient plasma metabolome data. By analyzing endothelial cell metabolism comprehensively, the pipeline identified critical therapeutic targets and biomarkers that can potentially contribute to the development of therapeutic strategies. Our study demonstrates the efficacy of integrating patient plasma metabolome data into computational models to analyze endothelial cell metabolism in disease contexts. This approach offers a deeper understanding of metabolic dysregulations and provides insights into diseases with metabolic components and potential treatments.

A Protocol for the Automatic Construction of Highly Curated Genome-Scale Models of Human Metabolism.

Genome-scale metabolic models (GEMs) have emerged as a tool to understand human metabolism from a holistic perspective with high relevance in the study of many diseases and in the metabolic engineering of human cell lines. GEM building relies on either automated processes that lack manual refinement and result in inaccurate models or manual curation, which is a time-consuming process that limits the continuous update of reliable GEMs. Here, we present a novel algorithm-aided protocol that overcomes these limitations and facilitates the continuous updating of highly curated GEMs. The algorithm enables the automatic curation and/or expansion of existing GEMs or generates a highly curated metabolic network based on current information retrieved from multiple databases in real time. This tool was applied to the latest reconstruction of human metabolism (Human1), generating a series of the human GEMs that improve and expand the reference model and generating the most extensive and comprehensive general reconstruction of human metabolism to date. The tool presented here goes beyond the current state of the art and paves the way for the automatic reconstruction of a highly curated, up-to-date GEM with high potential in computational biology as well as in multiple fields of biological science where metabolism is relevant.

A Multi-Level Systems Biology Analysis of Aldrin's Metabolic Effects on Prostate Cancer Cells.

Although numerous studies support a dose-effect relationship between Endocrine disruptors (EDs) and the progression and malignancy of tumors, the impact of a chronic exposure to non-lethal concentrations of EDs in cancer remains unknown. More specifically, a number of studies have reported the impact of Aldrin on a variety of cancer types, including prostate cancer. In previous studies, we demonstrated the induction of the malignant phenotype in DU145 prostate cancer (PCa) cells after a chronic exposure to Aldrin (an ED). Proteins are pivotal in the regulation and control of a variety of cellular processes. However, the mechanisms responsible for the impact of ED on PCa and the role of proteins in this process are not yet well understood. Here, two complementary computational approaches have been employed to investigate the molecular processes underlying the acquisition of malignancy in prostate cancer. First, the metabolic reprogramming associated with the chronic exposure to Aldrin in DU145 cells was studied by integrating transcriptomics and metabolomics via constraint-based metabolic modeling. Second, gene set enrichment analysis was applied to determine (i) altered regulatory pathways and (ii) the correlation between changes in the transcriptomic profile of Aldrin-exposed cells and tumor progression in various types of cancer. Experimental validation confirmed predictions revealing a disruption in metabolic and regulatory pathways. This alteration results in the modification of protein levels crucial in regulating triacylglyceride/cholesterol, linked to the malignant phenotype observed in Aldrin-exposed cells.

Endothelial Cell Phenotypes Demonstrate Different Metabolic Patterns and Predict Mortality in Trauma Patients.

In trauma patients, shock-induced endotheliopathy (SHINE) is associated with a poor prognosis. We have previously identified four metabolic phenotypes in a small cohort of trauma patients (N = 20) and displayed the intracellular metabolic profile of the endothelial cell by integrating quantified plasma metabolomic profiles into a genome-scale metabolic model (iEC-GEM). A retrospective observational study of 99 trauma patients admitted to a Level 1 Trauma Center. Mass spectrometry was conducted on admission samples of plasma metabolites. Quantified metabolites were analyzed by computational network analysis of the iEC-GEM. Four plasma metabolic phenotypes (A-D) were identified, of which phenotype D was associated with an increased injury severity score (p < 0.001); 90% (91.6%) of the patients who died within 72 h possessed this phenotype. The inferred EC metabolic patterns were found to be different between phenotype A and D. Phenotype D was unable to maintain adequate redox homeostasis. We confirm that trauma patients presented four metabolic phenotypes at admission. Phenotype D was associated with increased mortality. Different EC metabolic patterns were identified between phenotypes A and D, and the inability to maintain adequate redox balance may be linked to the high mortality.

Enhancing control systems of higher plant culture chambers via multilevel structural mechanistic modelling.

Modelling higher plant growth is of strategic interest for modern agriculture as well as for the development of bioregenerative life support systems for space applications, where crop growth is expected to play an essential role. The capability of constraint-based metabolic models to cope the diel dynamics of plants growth is integrated into a multilevel modelling approach including mass and energy transfer and enzyme kinetics. Lactuca sativa is used as an exemplary crop to validate, with experimental data, the approach presented as well as to design a novel model-based predictive control strategy embedding metabolic information. The proposed modelling strategy predicts with high accuracy the dynamics of gas exchange and the distribution of fluxes in the metabolic network whereas the control architecture presented can be useful to manage higher plants chambers and open new ways of merging metabolome and control algorithms.

Metabolic systems analysis identifies a novel mechanism contributing to shock in patients with endotheliopathy of trauma (EoT) involving thromboxane A2 and LTC4.

Purpose: Endotheliopathy of trauma (EoT), as defined by circulating levels of syndecan-1 ≥ 40 ng/mL, has been reported to be associated with significantly increased transfusion requirements and a doubled 30-day mortality. Increased shedding of the glycocalyx points toward the endothelial cell membrane composition as important for the clinical outcome being the rationale for this study. Results: The plasma metabolome of 95 severely injured trauma patients was investigated by mass spectrometry, and patients with EoT vs. non-EoT were compared by partial least square-discriminant analysis, identifying succinic acid as the top metabolite to differentiate EoT and non-EoT patients (VIP score = 3). EoT and non-EoT patients' metabolic flux profile was inferred by integrating the corresponding plasma metabolome data into a genome-scale metabolic network reconstruction analysis and performing a functional study of the metabolic capabilities of each group. Model predictions showed a decrease in cholesterol metabolism secondary to impaired mevalonate synthesis in EoT compared to non-EoT patients. Intracellular task analysis indicated decreased synthesis of thromboxanA2 and leukotrienes, as well as a lower carnitine palmitoyltransferase I activity in EoT compared to non-EoT patients. Sensitivity analysis also showed a significantly high dependence of eicosanoid-associated metabolic tasks on alpha-linolenic acid as unique to EoT patients. Conclusions: Model-driven analysis of the endothelial cells' metabolism identified potential novel targets as impaired thromboxane A2 and leukotriene synthesis in EoT patients when compared to non-EoT patients. Reduced thromboxane A2 and leukotriene availability in the microvasculature impairs vasoconstriction ability and may thus contribute to shock in EoT patients. These findings are supported by extensive scientific literature; however, further investigations are required on these findings.

Soft-sensors application for automated feeding control in high-throughput mammalian cell cultures.

The ever-increasing demand for biopharmaceuticals has created the need for improving the overall productivity of culture processes. One such operational concept that is considered is fed-batch operations as opposed to batch operations. However, optimal fed-batch operations require complete knowledge of the cell culture to optimize the culture conditions and the nutrients feeding. For example, when using high-throughput small-scale bioreactors to test multiple clones that do not behave the same, depletion or overfeeding of some key components can occur if the feeding strategy is not individually optimized. Over the recent years, various solutions for real-time measuring of the main cell culture metabolites have been proposed. Still, the complexity in the implementation of these techniques has limited their use. Soft-sensors present an opportunity to overcome these limitations by indirectly estimating these variables in real-time. This manuscript details the development of a new soft-sensor-based fed-batch strategy to maintain substrate concentration (glucose and glutamine) at optimal levels in small-scale multiparallel Chinese Hamster Ovary Cells cultures. Two alternatives to the standard feeding strategy were tested: an OUR soft-sensor-based strategy for glucose and glutamine (Strategy 1) and a dual OUR for glutamine and CO2 /alkali addition for glucose soft-sensor strategy (Strategy 2). The results demonstrated the applicability of the OUR soft-sensor-based strategy to optimize glucose and glutamine feedings, which yielded a 21% increase in final viable cell density (VCD) and a 31% in erythropoietin titer compared with the reference one. However, CO2 /alkali addition soft-sensor suffered from insufficient data to relate alkali addition with glucose consumption. As a result, the culture was overfed with glucose resulting in a 4% increase on final VCD, but a 9% decrease in final titer compared with the Reference Strategy.

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.

Stoichiometric gene-to-reaction associations enhance model-driven analysis performance: Metabolic response to chronic exposure to Aldrin in prostate cancer.

BACKGROUND: Genome-scale metabolic models (GSMM) integrating transcriptomics have been widely used to study cancer metabolism. This integration is achieved through logical rules that describe the association between genes, proteins, and reactions (GPRs). However, current gene-to-reaction formulation lacks the stoichiometry describing the transcript copies necessary to generate an active catalytic unit, which limits our understanding of how genes modulate metabolism. The present work introduces a new state-of-the-art GPR formulation that considers the stoichiometry of the transcripts (S-GPR). As case of concept, this novel gene-to-reaction formulation was applied to investigate the metabolic effects of the chronic exposure to Aldrin, an endocrine disruptor, on DU145 prostate cancer cells. To this aim we integrated the transcriptomic data from Aldrin-exposed and non-exposed DU145 cells through S-GPR or GPR into a human GSMM by applying different constraint-based-methods. RESULTS: Our study revealed a significant improvement of metabolite consumption/production predictions when S-GPRs are implemented. Furthermore, our computational analysis unveiled important alterations in carnitine shuttle and prostaglandine biosynthesis in Aldrin-exposed DU145 cells that is supported by bibliographic evidences of enhanced malignant phenotype. CONCLUSIONS: The method developed in this work enables a more accurate integration of gene expression data into model-driven methods. Thus, the presented approach is conceptually new and paves the way for more in-depth studies of aberrant cancer metabolism and other diseases with strong metabolic component with important environmental and clinical implications.

Minimizing Clonal Variation during Mammalian Cell Line Engineering for Improved Systems Biology Data Generation.

Mammalian cells are widely used to express genes for basic biology studies and biopharmaceuticals. Current methods for generation of engineered cell lines introduce high genomic and phenotypic diversity, which hamper studies of gene functions and discovery of novel cellular mechanisms. Here, we minimized clonal variation by integrating a landing pad for recombinase-mediated cassette exchange site-specifically into the genome of CHO cells using CRISPR and generated subclones expressing four different recombinant proteins. The subclones showed low clonal variation with high consistency in growth, transgene transcript levels and global transcriptional response to recombinant protein expression, enabling improved studies of the impact of transgenes on the host transcriptome. Little variation over time in subclone phenotypes and transcriptomes was observed when controlling environmental culture conditions. The platform enables robust comparative studies of genome engineered CHO cell lines and can be applied to other mammalian cells for diverse biological, biomedical and biotechnological applications.

Model-driven discovery of long-chain fatty acid metabolic reprogramming in heterogeneous prostate cancer cells.

Epithelial-mesenchymal-transition promotes intra-tumoral heterogeneity, by enhancing tumor cell invasiveness and promoting drug resistance. We integrated transcriptomic data for two clonal subpopulations from a prostate cancer cell line (PC-3) into a genome-scale metabolic network model to explore their metabolic differences and potential vulnerabilities. In this dual cell model, PC-3/S cells express Epithelial-mesenchymal-transition markers and display high invasiveness and low metastatic potential, while PC-3/M cells present the opposite phenotype and higher proliferative rate. Model-driven analysis and experimental validations unveiled a marked metabolic reprogramming in long-chain fatty acids metabolism. While PC-3/M cells showed an enhanced entry of long-chain fatty acids into the mitochondria, PC-3/S cells used long-chain fatty acids as precursors of eicosanoid metabolism. We suggest that this metabolic reprogramming endows PC-3/M cells with augmented energy metabolism for fast proliferation and PC-3/S cells with increased eicosanoid production impacting angiogenesis, cell adhesion and invasion. PC-3/S metabolism also promotes the accumulation of docosahexaenoic acid, a long-chain fatty acid with antiproliferative effects. The potential therapeutic significance of our model was supported by a differential sensitivity of PC-3/M cells to etomoxir, an inhibitor of long-chain fatty acid transport to the mitochondria.

Unveiling the Metabolic Changes on Muscle Cell Metabolism Underlying p-Phenylenediamine Toxicity.

Rhabdomyolysis is a disorder characterized by acute damage of the sarcolemma of the skeletal muscle leading to release of potentially toxic muscle cell components into the circulation, most notably creatine phosphokinase (CK) and myoglobulin, and is frequently accompanied by myoglobinuria. In the present work, we evaluated the toxicity of p-phenylenediamine (PPD), a main component of hair dyes which is reported to induce rhabdomyolysis. We studied the metabolic effect of this compound in vivo with Wistar rats and in vitro with C2C12 muscle cells. To this aim we have combined multi-omic experimental measurements with computational approaches using model-driven methods. The integrative study presented here has unveiled the metabolic disorders associated to PPD exposure that may underlay the aberrant metabolism observed in rhabdomyolys disease. Animals treated with lower doses of PPD (10 and 20 mg/kg) showed depressed activity and myoglobinuria after 10 h of treatment. We measured the serum levels of aspartate aminotransferase (AST), alanine aminotransferase (ALT), and creatine kinase (CK) in rats after 24, 48, and 72 h of PPD exposure. At all times, treatment with PPD at higher doses (40 and 60 mg/kg) showed an increase of AST and ALT, and also an increase of lactate dehydrogenase (LDH) and CK after 24 h. Blood packed cell volume and hemoglobin levels, as well as organs weight at 48 and 72 h, were also measured. No significant differences were observed in these parameters under any condition. PPD induce cell cycle arrest in S phase and apoptosis (40% or early apoptotic cells) on mus musculus mouse C2C12 cells after 24 h of treatment. Incubation of mus musculus mouse C2C12 cells with [1,2-13C2]-glucose during 24 h, subsequent quantification of 13C isotopologues distribution in key metabolites of glucose metabolic network and a computational fluxomic analysis using in-house developed software (Isodyn) showed that PPD is inhibiting glycolysis, non-oxidative pentose phosphate pathway, glycogen turnover, and ATPAse reaction leading to a reduction in ATP synthesis. These findings unveil the glucose metabolism collapse, which is consistent with a decrease in cell viability observed in PPD-treated C2C12 cells and with the myoglubinuria and other effects observed in Wistar Rats treated with PPD. These findings shed new light on muscle dysfunction associated to PPD exposure, opening new avenues for cost-effective therapies in Rhabdomyolysis disease.

Molecular mechanisms underlying COPD-muscle dysfunction unveiled through a systems medicine approach.

MOTIVATION: Skeletal muscle dysfunction is a systemic effect in one-third of patients with chronic obstructive pulmonary disease (COPD), characterized by high reactive-oxygen-species (ROS) production and abnormal endurance training-induced adaptive changes. However, the role of ROS in COPD remains unclear, not least because of the lack of appropriate tools to study multifactorial diseases. RESULTS: We describe a discrete model-driven method combining mechanistic and probabilistic approaches to decipher the role of ROS on the activity state of skeletal muscle regulatory network, assessed before and after an 8-week endurance training program in COPD patients and healthy subjects. In COPD, our computational analysis indicates abnormal training-induced regulatory responses leading to defective tissue remodeling and abnormal energy metabolism. Moreover, we identified tnf, insr, inha and myc as key regulators of abnormal training-induced adaptations in COPD. The tnf-insr pair was identified as a promising target for therapeutic interventions. Our work sheds new light on skeletal muscle dysfunction in COPD, opening new avenues for cost-effective therapies. It overcomes limitations of previous computational approaches showing high potential for the study of other multi-factorial diseases such as diabetes or cancer. CONTACT: jroca@clinic.ub.es or martacascante@ub.eduSupplementary information: Supplementary data are available at Bioinformatics online.

Metabolic Reprogramming and Dependencies Associated with Epithelial Cancer Stem Cells Independent of the Epithelial-Mesenchymal Transition Program.

In solid tumors, cancer stem cells (CSCs) can arise independently of epithelial-mesenchymal transition (EMT). In spite of recent efforts, the metabolic reprogramming associated with CSC phenotypes uncoupled from EMT is poorly understood. Here, by using metabolomic and fluxomic approaches, we identify major metabolic profiles that differentiate metastatic prostate epithelial CSCs (e-CSCs) from non-CSCs expressing a stable EMT. We have found that the e-CSC program in our cellular model is characterized by a high plasticity in energy substrate metabolism, including an enhanced Warburg effect, a greater carbon and energy source flexibility driven by fatty acids and amino acid metabolism and an essential reliance on the proton buffering capacity conferred by glutamine metabolism. An analysis of transcriptomic data yielded a metabolic gene signature for our e-CSCs consistent with the metabolomics and fluxomics analyses that correlated with tumor progression and metastasis in prostate cancer and in 11 additional cancer types. Interestingly, an integrated metabolomics, fluxomics, and transcriptomics analysis allowed us to identify key metabolic players regulated at the post-transcriptional level, suggesting potential biomarkers and therapeutic targets to effectively forestall metastasis. Stem Cells 2016;34:1163-1176.

Systems Medicine: from molecular features and models to the clinic in COPD.

BACKGROUND AND HYPOTHESIS: Chronic Obstructive Pulmonary Disease (COPD) patients are characterized by heterogeneous clinical manifestations and patterns of disease progression. Two major factors that can be used to identify COPD subtypes are muscle dysfunction/wasting and co-morbidity patterns. We hypothesized that COPD heterogeneity is in part the result of complex interactions between several genes and pathways. We explored the possibility of using a Systems Medicine approach to identify such pathways, as well as to generate predictive computational models that may be used in clinic practice. OBJECTIVE AND METHOD: Our overarching goal is to generate clinically applicable predictive models that characterize COPD heterogeneity through a Systems Medicine approach. To this end we have developed a general framework, consisting of three steps/objectives: (1) feature identification, (2) model generation and statistical validation, and (3) application and validation of the predictive models in the clinical scenario. We used muscle dysfunction and co-morbidity as test cases for this framework. RESULTS: In the study of muscle wasting we identified relevant features (genes) by a network analysis and generated predictive models that integrate mechanistic and probabilistic models. This allowed us to characterize muscle wasting as a general de-regulation of pathway interactions. In the co-morbidity analysis we identified relevant features (genes/pathways) by the integration of gene-disease and disease-disease associations. We further present a detailed characterization of co-morbidities in COPD patients that was implemented into a predictive model. In both use cases we were able to achieve predictive modeling but we also identified several key challenges, the most pressing being the validation and implementation into actual clinical practice. CONCLUSIONS: The results confirm the potential of the Systems Medicine approach to study complex diseases and generate clinically relevant predictive models. Our study also highlights important obstacles and bottlenecks for such approaches (e.g. data availability and normalization of frameworks among others) and suggests specific proposals to overcome them.

Workforce preparation: the Biohealth computing model for Master and PhD students.

The article addresses the strategic role of workforce preparation in the process of adoption of Systems Medicine as a driver of biomedical research in the new health paradigm. It reports on relevant initiatives, like CASyM, fostering Systems Medicine at EU level. The chapter focuses on the BioHealth Computing Program as a reference for multidisciplinary training of future systems-oriented researchers describing the productive interactions with the Synergy-COPD project.

Cancer cell metabolism as new targets for novel designed therapies.

Metabolic processes are altered in cancer cells, which obtain advantages from this metabolic reprogramming in terms of energy production and synthesis of biomolecules that sustain their uncontrolled proliferation. Due to the conceptual progresses in the last decade, metabolic reprogramming was recently included as one of the new hallmarks of cancer. The advent of high-throughput technologies to amass an abundance of omic data, together with the development of new computational methods that allow the integration and analysis of omic data by using genome-scale reconstructions of human metabolism, have increased and accelerated the discovery and development of anticancer drugs and tumor-specific metabolic biomarkers. Here we review and discuss the latest advances in the context of metabolic reprogramming and the future in cancer research.

Relevance of the MEK/ERK signaling pathway in the metabolism of activated macrophages: a metabolomic approach.

The activation of immune cells in response to a pathogen involves a succession of signaling events leading to gene and protein expression, which requires metabolic changes to match the energy demands. The metabolic profile associated with the MAPK cascade (ERK1/2, p38, and JNK) in macrophages was studied, and the effect of its inhibition on the specific metabolic pattern of LPS stimulation was characterized. A [1,2-[(13)C](2)]glucose tracer-based metabolomic approach was used to examine the metabolic flux distribution in these cells after MEK/ERK inhibition. Bioinformatic tools were used to analyze changes in mass isotopomer distribution and changes in glucose and glutamine consumption and lactate production in basal and LPS-stimulated conditions in the presence and absence of the selective inhibitor of the MEK/ERK cascade, PD325901. Results showed that PD325901-mediated ERK1/2 inhibition significantly decreased glucose consumption and lactate production but did not affect glutamine consumption. These changes were accompanied by a decrease in the glycolytic flux, consistent with the observed decrease in fructose-2,6-bisphosphate concentration. The oxidative and nonoxidative pentose phosphate pathways and the ratio between them also decreased. However, tricarboxylic acid cycle flux did not change significantly. LPS activation led to the opposite responses, although all of these were suppressed by PD325901. However, LPS also induced a small decrease in pentose phosphate pathway fluxes and an increase in glutamine consumption that were not affected by PD325901. We concluded that inhibition of the MEK/ERK cascade interferes with central metabolism, and this cross-talk between signal transduction and metabolism also occurs in the presence of LPS.

Knowledge management for systems biology a general and visually driven framework applied to translational medicine.

BACKGROUND: To enhance our understanding of complex biological systems like diseases we need to put all of the available data into context and use this to detect relations, pattern and rules which allow predictive hypotheses to be defined. Life science has become a data rich science with information about the behaviour of millions of entities like genes, chemical compounds, diseases, cell types and organs, which are organised in many different databases and/or spread throughout the literature. Existing knowledge such as genotype-phenotype relations or signal transduction pathways must be semantically integrated and dynamically organised into structured networks that are connected with clinical and experimental data. Different approaches to this challenge exist but so far none has proven entirely satisfactory. RESULTS: To address this challenge we previously developed a generic knowledge management framework, BioXM™, which allows the dynamic, graphic generation of domain specific knowledge representation models based on specific objects and their relations supporting annotations and ontologies. Here we demonstrate the utility of BioXM for knowledge management in systems biology as part of the EU FP6 BioBridge project on translational approaches to chronic diseases. From clinical and experimental data, text-mining results and public databases we generate a chronic obstructive pulmonary disease (COPD) knowledge base and demonstrate its use by mining specific molecular networks together with integrated clinical and experimental data. CONCLUSIONS: We generate the first semantically integrated COPD specific public knowledge base and find that for the integration of clinical and experimental data with pre-existing knowledge the configuration based set-up enabled by BioXM reduced implementation time and effort for the knowledge base compared to similar systems implemented as classical software development projects. The knowledgebase enables the retrieval of sub-networks including protein-protein interaction, pathway, gene--disease and gene--compound data which are used for subsequent data analysis, modelling and simulation. Pre-structured queries and reports enhance usability; establishing their use in everyday clinical settings requires further simplification with a browser based interface which is currently under development.