Classification of Neisseria meningitidis genomes with a bag-of-words approach and machine learning

Whole genome sequencing of bacteria is important to enable strain classification. Using entire genomes as an input to machine learning (ML) models would allow rapid classification of strains while using information from multiple genetic elements. We developed a "bag-of-words" approach to encode, using SentencePiece or k-mer tokenization, entire bacterial genomes and analyze these with ML. Initial model selection identified SentencePiece with 8,000 and 32,000 words as the best approach for genome tokenization. We then classified in Neisseria meningitidis genomes the capsule B group genotype with 99.6% accuracy and the multifactor invasive phenotype with 90.2% accuracy, in an independent test set. Subsequently, in silico knockouts of 2,808 genes confirmed that the ML model predictions aligned with our current understanding of the underlying biology. To our knowledge, this is the first ML method using entire bacterial genomes to classify strains and identify genes considered relevant by the classifier.

Controlling astrocyte-mediated synaptic pruning signals for schizophrenia drug repurposing with deep graph networks

Schizophrenia is a debilitating psychiatric disorder, leading to both physical and social morbidity. Worldwide 1% of the population is struggling with the disease, with 100,000 new cases annually only in the United States. Despite its importance, the goal of finding effective treatments for schizophrenia remains a challenging task, and previous work conducted expensive large-scale phenotypic screens. This work investigates the benefits of Machine Learning for graphs to optimize drug phenotypic screens and predict compounds that mitigate abnormal brain reduction induced by excessive glial phagocytic activity in schizophrenia subjects. Given a compound and its concentration as input, we propose a method that predicts a score associated with three possible compound effects, i.e., reduce, increase, or not influence phagocytosis. We leverage a high-throughput screening to prove experimentally that our method achieves good generalization capabilities. The screening involves 2218 compounds at five different concentrations. Then, we analyze the usability of our approach in a practical setting, i.e., prioritizing the selection of compounds in the SWEETLEAD library. We provide a list of 64 compounds from the library that have the most potential clinical utility for glial phagocytosis mitigation. Lastly, we propose a novel approach to computationally validate their utility as possible therapies for schizophrenia.

Literature Mining and Mechanistic Graphical Modelling to Improve mRNA Vaccine Platforms

RNA vaccines represent a milestone in the history of vaccinology. They provide several advantages over more traditional approaches to vaccine development, showing strong immunogenicity and an overall favorable safety profile. While preclinical testing has provided some key insights on how RNA vaccines interact with the innate immune system, their mechanism of action appears to be fragmented amid the literature, making it difficult to formulate new hypotheses to be tested in clinical settings and ultimately improve this technology platform. Here, we propose a systems biology approach, based on the combination of literature mining and mechanistic graphical modeling, to consolidate existing knowledge around mRNA vaccines mode of action and enhance the translatability of preclinical hypotheses into clinical evidence. A Natural Language Processing (NLP) pipeline for automated knowledge extraction retrieved key biological evidences that were joined into an interactive mechanistic graphical model representing the chain of immune events induced by mRNA vaccines administration. The achieved mechanistic graphical model will help the design of future experiments, foster the generation of new hypotheses and set the basis for the development of mathematical models capable of simulating and predicting the immune response to mRNA vaccines.

Early outcome detection for COVID-19 patients

With the outbreak of COVID-19 exerting a strong pressure on hospitals and health facilities, clinical decision support systems based on predictive models can help to effectively improve the management of the pandemic. We present a method for predicting mortality for COVID-19 patients. Starting from a large number of clinical variables, we select six of them with largest predictive power, using a feature selection method based on genetic algorithms and starting from a set of COVID-19 patients from the first wave. The algorithm is designed to reduce the impact of missing values in the set of variables measured, and consider only variables that show good accuracy on validation data. The final predictive model provides accuracy larger than 85% on test data, including a new patient cohort from the second COVID-19 wave, and on patients with imputed missing values. The selected clinical variables are confirmed to be relevant by recent literature on COVID-19.

QSPcc reduces bottlenecks in computational model simulations

Mathematical models have grown in size and complexity becoming often computationally intractable. In sensitivity analysis and optimization phases, critical for tuning, validation and qualification, these models may be run thousands of times. Scientific programming languages popular for prototyping, such as MATLAB and R, can be a bottleneck in terms of performance. Here we show a compiler-based approach, designed to be universal at handling engineering and life sciences modeling styles, that automatically translates models into fast C code. At first QSPcc is demonstrated to be crucial in enabling the research on otherwise intractable Quantitative Systems Pharmacology models, such as in rare Lysosomal Storage Disorders. To demonstrate the full value in seamlessly accelerating, or enabling, the R&D efforts in natural sciences, we then benchmark QSPcc against 8 solutions on 24 real-world projects from different scientific fields. With speed-ups of 22000x peak, and 1605x arithmetic mean, our results show consistent superior performances.

Lombardo and colleagues present QSPcc, a computational code compiler designed to convert code from popular scientific programming languages, such as MATLAB or R, into fast-running C code. This reduces the computational load required for complex modelling approaches and reduces user investment learning additional complex languages.

A robust computational pipeline for model-based and data-driven phenotype clustering

MOTIVATION

Precision medicine is a promising field that proposes, in contrast to a one-size-fits-all approach, the tailoring of medical decisions, treatments or products. In this context, it is crucial to introduce innovative methods to stratify a population of patients on the basis of an accurate system-level knowledge of the disease. This is particularly important in very challenging conditions, where the use of standard statistical methods can be prevented by poor data availability or by the need of oversimplifying the processes regulating a complex disease.

RESULTS

We define an innovative method for phenotype classification that combines experimental data and a mathematical description of the disease biology. The methodology exploits the mathematical model for inferring additional subject features relevant for the classification. Finally, the algorithm identifies the optimal number of clusters and classifies the samples on the basis of a subset of the features estimated during the model fit. We tested the algorithm in two test cases: an in silico case in the context of dyslipidemia, a complex disease for which a large population of patients has been generated, and a clinical test case, in the context of a lysosomal rare disorder, for which the amount of available data was limited. In both the scenarios, our methodology proved to be accurate and robust, and allowed the inference of an additional phenotype division that the experimental data did not show.

AVAILABILITY AND IMPLEMENTATION

The code to reproduce the in silico results has been implemented in MATLAB v.2017b and it is available in the Supplementary Material.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

DNA methylation during human adipogenesis and the impact of fructose

BACKGROUND

Increased adipogenesis and altered adipocyte function contribute to the development of obesity and associated comorbidities. Fructose modified adipocyte metabolism compared to glucose, but the regulatory mechanisms and consequences for obesity are unknown. Genome-wide methylation and global transcriptomics in SGBS pre-adipocytes exposed to 0, 2.5, 5, and 10 mM fructose, added to a 5-mM glucose-containing medium, were analyzed at 0, 24, 48, 96, 192, and 384 h following the induction of adipogenesis.

RESULTS

Time-dependent changes in DNA methylation compared to baseline (0 h) occurred during the final maturation of adipocytes, between 192 and 384 h. Larger percentages (0.1% at 192 h, 3.2% at 384 h) of differentially methylated regions (DMRs) were found in adipocytes differentiated in the glucose-containing control media compared to adipocytes differentiated in fructose-supplemented media (0.0006% for 10 mM, 0.001% for 5 mM, and 0.005% for 2.5 mM at 384 h). A total of 1437 DMRs were identified in 5237 differentially expressed genes at 384 h post-induction in glucose-containing (5 mM) control media. The majority of them inversely correlated with the gene expression, but 666 regions were positively correlated to the gene expression.

CONCLUSIONS

Our studies demonstrate that DNA methylation regulates or marks the transformation of morphologically differentiating adipocytes (seen at 192 h), to the more mature and metabolically robust adipocytes (as seen at 384 h) in a genome-wide manner. Lower (2.5 mM) concentrations of fructose have the most robust effects on methylation compared to higher concentrations (5 and 10 mM), suggesting that fructose may be playing a signaling/regulatory role at lower concentrations of fructose and as a substrate at higher concentrations.

Computational approaches to understanding nutrient metabolism and metabolic disorders

Computational methods are becoming more and more essential to elucidate biological systems. Many different approaches exist with pros and cons. This paper reviews the most useful technologies focusing on nutrient metabolism and metabolic disorders. Space limitation prevents from exploring the examples in details, but pointers to the relevant papers are reported.

Multi-omics integration-a comparison of unsupervised clustering methodologies

With the recent developments in the field of multi-omics integration, the interest in factors such as data preprocessing, choice of the integration method and the number of different omics considered had increased. In this work, the impact of these factors is explored when solving the problem of sample classification, by comparing the performances of five unsupervised algorithms: Multiple Canonical Correlation Analysis, Multiple Co-Inertia Analysis, Multiple Factor Analysis, Joint and Individual Variation Explained and Similarity Network Fusion. These methods were applied to three real data sets taken from literature and several ad hoc simulated scenarios to discuss classification performance in different conditions of noise and signal strength across the data types. The impact of experimental design, feature selection and parameter training has been also evaluated to unravel important conditions that can affect the accuracy of the result.

A network-based approach to identify deregulated pathways and drug effects in metabolic syndrome

Metabolic syndrome is a pathological condition characterized by obesity, hyperglycemia, hypertension, elevated levels of triglycerides and low levels of high-density lipoprotein cholesterol that increase cardiovascular disease risk and type 2 diabetes. Although numerous predisposing genetic risk factors have been identified, the biological mechanisms underlying this complex phenotype are not fully elucidated. Here we introduce a systems biology approach based on network analysis to investigate deregulated biological processes and subsequently identify drug repurposing candidates. A proximity score describing the interaction between drugs and pathways is defined by combining topological and functional similarities. The results of this computational framework highlight a prominent role of the immune system in metabolic syndrome and suggest a potential use of the BTK inhibitor ibrutinib as a novel pharmacological treatment. An experimental validation using a high fat diet-induced obesity model in zebrafish larvae shows the effectiveness of ibrutinib in lowering the inflammatory load due to macrophage accumulation.

Regulatory Crosstalk of Doxorubicin, Estradiol and TNFα Combined Treatment in Breast Cancer-derived Cell Lines

We present a new model of ESR1 network regulation based on analysis of Doxorubicin, Estradiol, and TNFα combination treatment in MCF-7. We used Doxorubicin as a therapeutic agent, TNFα as marker and mediator of an inflammatory microenvironment and 17β-Estradiol (E2) as an agonist of Estrogen Receptors, known predisposing factor for hormone-driven breast cancer, whose pharmacological inhibition reduces the risk of breast cancer recurrence. Based on the results of transcriptomics analysis, we found 71 differentially expressed genes that are specific for the combination treatment with Doxorubicin + Estradiol + TNFα in comparison with single or double treatments. The responsiveness to the triple treatment was examined for seven genes by qPCR, of which six were validated, and then extended to four additional cell lines differing for p53 and/or ER status. The results of differential regulation enrichment analysis highlight the role of the ESR1 network that included 36 of 71 specific differentially expressed genes. We propose that the combined activation of p53 and NF-kB transcription factors significantly influences ligand-dependent, ER-driven transcriptional responses, also of the ESR1 gene itself. These results provide a model of coordinated interaction of TFs to explain the Doxorubicin, E2 and TNFα induced repression mechanisms.

Stochastic simulation algorithms for computational systems biology: Exact, approximate, and hybrid methods

Nowadays, mathematical modeling is playing a key role in many different research fields. In the context of system biology, mathematical models and their associated computer simulations constitute essential tools of investigation. Among the others, they provide a way to systematically analyze systems perturbations, develop hypotheses to guide the design of new experimental tests, and ultimately assess the suitability of specific molecules as novel therapeutic targets. To these purposes, stochastic simulation algorithms (SSAs) have been introduced for numerically simulating the time evolution of a well-stirred chemically reacting system by taking proper account of the randomness inherent in such a system. In this work, we review the main SSAs that have been introduced in the context of exact, approximate, and hybrid stochastic simulation. Specifically, we will introduce the direct method (DM), the first reaction method (FRM), the next reaction method (NRM) and the rejection-based SSA (RSSA) in the area of exact stochastic simulation. We will then present the τ-leaping method and the chemical Langevin method in the area of approximate stochastic simulation and an implementation of the hybrid RSSA (HRSSA) in the context of hybrid stochastic-deterministic simulation. Finally, we will consider the model of the sphingolipid metabolism to provide an example of application of SSA to computational system biology by exemplifying how different simulation strategies may unveil different insights into the investigated biological phenomenon. This article is categorized under: Models of Systems Properties and Processes > Mechanistic Models Analytical and Computational Methods > Computational Methods.

Cross-disease analysis of Alzheimer's disease and type-2 Diabetes highlights the role of autophagy in the pathophysiology of two highly comorbid diseases

Evidence is accumulating that the main chronic diseases of aging Alzheimer's disease (AD) and type-2 diabetes mellitus (T2DM) share common pathophysiological mechanisms. This study aimed at applying systems biology approaches to increase the knowledge of the shared molecular pathways underpinnings of AD and T2DM. We analysed transcriptomic data of post-mortem AD and T2DM human brains to obtain disease signatures of AD and T2DM and combined them with protein-protein interaction information to construct two disease-specific networks. The overlapping AD/T2DM network proteins were then used to extract the most representative Gene Ontology biological process terms. The expression of genes identified as relevant was studied in two AD models, 3xTg-AD and ApoE3/ApoE4 targeted replacement mice. The present transcriptomic data analysis revealed a principal role for autophagy in the molecular basis of both AD and T2DM. Our experimental validation in mouse AD models confirmed the role of autophagy-related genes. Among modulated genes, Cyclin-Dependent Kinase Inhibitor 1B, Autophagy Related 16-Like 2, and insulin were highlighted. In conclusion, the present investigation revealed autophagy as the central dys-regulated pathway in highly co-morbid diseases such as AD and T2DM allowing the identification of specific genes potentially involved in disease pathophysiology which could become novel targets for therapeutic intervention.

A comparison of deterministic and stochastic approaches for sensitivity analysis in computational systems biology

With the recent rising application of mathematical models in the field of computational systems biology, the interest in sensitivity analysis methods had increased. The stochastic approach, based on chemical master equations, and the deterministic approach, based on ordinary differential equations (ODEs), are the two main approaches for analyzing mathematical models of biochemical systems. In this work, the performance of these approaches to compute sensitivity coefficients is explored in situations where stochastic and deterministic simulation can potentially provide different results (systems with unstable steady states, oscillators with population extinction and bistable systems). We consider two methods in the deterministic approach, namely the direct differential method and the finite difference method, and five methods in the stochastic approach, namely the Girsanov transformation, the independent random number method, the common random number method, the coupled finite difference method and the rejection-based finite difference method. The reviewed methods are compared in terms of sensitivity values and computational time to identify differences in outcome that can highlight conditions in which one approach performs better than the other.

Correction: Identification of an early transcriptomic signature of insulin resistance and related diseases in lymphomonocytes of healthy subjects

[This corrects the article DOI: 10.1371/journal.pone.0182559.].

Efficient finite-difference method for computing sensitivities of biochemical reactions

Sensitivity analysis of biochemical reactions aims at quantifying the dependence of the reaction dynamics on the reaction rates. The computation of the parameter sensitivities, however, poses many computational challenges when taking stochastic noise into account. This paper proposes a new finite-difference method for efficiently computing sensitivities of biochemical reactions. We employ propensity bounds of reactions to couple the simulation of the nominal and perturbed processes. The exactness of the simulation is preserved by applying the rejection-based mechanism. For each simulation step, the nominal and perturbed processes under our coupling strategy are synchronized and often jump together, increasing their positive correlation and hence reducing the variance of the estimator. The distinctive feature of our approach in comparison with existing coupling approaches is that it only needs to maintain a single data structure storing propensity bounds of reactions during the simulation of the nominal and perturbed processes. Our approach allows to compute sensitivities of many reaction rates simultaneously. Moreover, the data structure does not require to be updated frequently, hence improving the computational cost. This feature is especially useful when applied to large reaction networks. We benchmark our method on biological reaction models to prove its applicability and efficiency.

A reverse metabolic approach to weaning: in silico identification of immune-beneficial infant gut bacteria, mining their metabolism for prebiotic feeds and sourcing these feeds in the natural product space

BACKGROUND

Weaning is a period of marked physiological change. The introduction of solid foods and the changes in milk consumption are accompanied by significant gastrointestinal, immune, developmental, and microbial adaptations. Defining a reduced number of infections as the desired health benefit for infants around weaning, we identified in silico (i.e., by advanced public domain mining) infant gut microbes as potential deliverers of this benefit. We then investigated the requirements of these bacteria for exogenous metabolites as potential prebiotic feeds that were subsequently searched for in the natural product space.

RESULTS

Using public domain literature mining and an in silico reverse metabolic approach, we constructed probiotic-prebiotic-food associations, which can guide targeted feeding of immune health-beneficial microbes by weaning food; analyzed competition and synergy for (prebiotic) nutrients between selected microbes; and translated this information into designing an experimental complementary feed for infants enrolled in a pilot clinical trial ( http://www.nourishtoflourish.auckland.ac.nz/ ).

CONCLUSIONS

In this study, we applied a benefit-oriented microbiome research strategy for enhanced early-life immune health. We extended from "classical" to molecular nutrition aiming to identify nutrients, bacteria, and mechanisms that point towards targeted feeding to improve immune health in infants around weaning. Here, we present the systems biology-based approach we used to inform us on the most promising prebiotic combinations known to support growth of beneficial gut bacteria ("probiotics") in the infant gut, thereby favorably promoting development of the immune system.

Incorporating extrinsic noise into the stochastic simulation of biochemical reactions: A comparison of approaches

The stochastic simulation algorithm (SSA) has been widely used for simulating biochemical reaction networks. SSA is able to capture the inherently intrinsic noise of the biological system, which is due to the discreteness of species population and to the randomness of their reciprocal interactions. However, SSA does not consider other sources of heterogeneity in biochemical reaction systems, which are referred to as extrinsic noise. Here, we extend two simulation approaches, namely, the integration-based method and the rejection-based method, to take extrinsic noise into account by allowing the reaction propensities to vary in time and state dependent manner. For both methods, new efficient implementations are introduced and their efficiency and applicability to biological models are investigated. Our numerical results suggest that the rejection-based method performs better than the integration-based method when the extrinsic noise is considered.

Combined use of protein biomarkers and network analysis unveils deregulated regulatory circuits in Duchenne muscular dystrophy

Although the genetic basis of Duchenne muscular dystrophy has been known for almost thirty years, the cellular and molecular mechanisms characterizing the disease are not completely understood and an efficacious treatment remains to be developed. In this study we analyzed proteomics data obtained with the SomaLogic technology from blood serum of a cohort of patients and matched healthy subjects. We developed a workflow based on biomarker identification and network-based pathway analysis that allowed us to describe different deregulated pathways. In addition to muscle-related functions, we identified other biological processes such as apoptosis, signaling in the immune system and neurotrophin signaling as significantly modulated in patients compared with controls. Moreover, our network-based analysis identified the involvement of FoxO transcription factors as putative regulators of different pathways. On the whole, this study provided a global view of the molecular processes involved in Duchenne muscular dystrophy that are decipherable from serum proteome.

Clinical and Vitamin Response to a Short-Term Multi-Micronutrient Intervention in Brazilian Children and Teens: From Population Data to Interindividual Responses

SCOPE

Micronutrients are in small amounts in foods, act in concert, and require variable amounts of time to see changes in health and risk for disease. These first principles are incorporated into an intervention study designed to develop new experimental strategies for setting target recommendations for food bioactives for populations and individuals.

METHODS AND RESULTS

A 6-week multivitamin/mineral intervention is conducted in 9-13 year olds. Participants (136) are (i) their own control (n-of-1); (ii) monitored for compliance; (iii) measured for 36 circulating vitamin forms, 30 clinical, anthropometric, and food intake parameters at baseline, post intervention, and following a 6-week washout; and (iv) had their ancestry accounted for as modifier of vitamin baseline or response. The same intervention is repeated the following year (135 participants). Most vitamins respond positively and many clinical parameters change in directions consistent with improved metabolic health to the intervention. Baseline levels of any metabolite predict its own response to the intervention. Elastic net penalized regression models are identified, and significantly predict response to intervention on the basis of multiple vitamin/clinical baseline measures.

CONCLUSIONS

The study design, computational methods, and results are a step toward developing recommendations for optimizing vitamin levels and health parameters for individuals.

A computationally driven analysis of the polyphenol-protein interactome

Polyphenol-rich foods are part of many nutritional interventions aimed at improving health and preventing cardiometabolic diseases (CMDs). Polyphenols have oxidative, inflammatory, and/or metabolic effects. Research into the chemistry and biology of polyphenol bioactives is prolific but knowledge of their molecular interactions with proteins is limited. We mined public data to (i) identify proteins that interact with or metabolize polyphenols, (ii) mapped these proteins to pathways and networks, and (iii) annotated functions enriched within the resulting polyphenol-protein interactome. A total of 1,395 polyphenols and their metabolites were retrieved (using Phenol-Explorer and Dictionary of Natural Products) of which 369 polyphenols interacted with 5,699 unique proteins in 11,987 interactions as annotated in STITCH, Pathway Commons, and BindingDB. Pathway enrichment analysis using the KEGG repository identified a broad coverage of significant pathways of low specificity to particular polyphenol (sub)classes. When compared to drugs or micronutrients, polyphenols have pleiotropic effects across many biological processes related to metabolism and CMDs. These systems-wide effects were also found in the protein interactome of the polyphenol-rich citrus fruits, used as a case study. In sum, these findings provide a knowledgebase for identifying polyphenol classes (and polyphenol-rich foods) that individually or in combination influence metabolism.

Identification of an early transcriptomic signature of insulin resistance and related diseases in lymphomonocytes of healthy subjects

Insulin resistance is considered to be a pathogenetic mechanism in several and diverse diseases (e.g. type 2 diabetes, atherosclerosis) often antedating them in apparently healthy subjects. The aim of this study is to investigate with a microarray based approach whether IR per se is characterized by a specific pattern of gene expression. For this purpose we analyzed the transcriptomic profile of peripheral blood mononuclear cells in two groups (10 subjects each) of healthy individuals, with extreme insulin resistance or sensitivity, matched for BMI, age and gender, selected within the MultiKnowledge Study cohort (n = 148). Data were analyzed with an ad-hoc rank-based classification method. 321 genes composed the gene set distinguishing the insulin resistant and sensitive groups, within which the "Adrenergic signaling in cardiomyocytes" KEGG pathway was significantly represented, suggesting a pattern of increased intracellular cAMP and Ca2+, and apoptosis in the IR group. The same pathway allowed to discriminate between insulin resistance and insulin sensitive subjects with BMI >25, supporting his role as a biomarker of IR. Moreover, ASCM pathway harbored biomarkers able to distinguish healthy and diseased subjects (from publicly available data sets) in IR-related diseases involving excitable cells: type 2 diabetes, chronic heart failure, and Alzheimer's disease. The altered gene expression profile of the ASCM pathway is an early molecular signature of IR and could provide a common molecular pathogenetic platform for IR-related disorders, possibly representing an important aid in the efforts aiming at preventing, early detecting and optimally treating IR-related diseases.

Efficient Constant-Time Complexity Algorithm for Stochastic Simulation of Large Reaction Networks

Exact stochastic simulation is an indispensable tool for a quantitative study of biochemical reaction networks. The simulation realizes the time evolution of the model by randomly choosing a reaction to fire and update the system state according to a probability that is proportional to the reaction propensity. Two computationally expensive tasks in simulating large biochemical networks are the selection of next reaction firings and the update of reaction propensities due to state changes. We present in this work a new exact algorithm to optimize both of these simulation bottlenecks. Our algorithm employs the composition-rejection on the propensity bounds of reactions to select the next reaction firing. The selection of next reaction firings is independent of the number reactions while the update of propensities is skipped and performed only when necessary. It therefore provides a favorable scaling for the computational complexity in simulating large reaction networks. We benchmark our new algorithm with the state of the art algorithms available in literature to demonstrate its applicability and efficiency.

A hybrid stochastic model of folate-mediated one-carbon metabolism: Effect of the common C677T MTHFR variant on de novo thymidylate biosynthesis

Folate-mediated one-carbon metabolism (FOCM) is an interconnected network of metabolic pathways, including those required for the de novo synthesis of dTMP and purine nucleotides and for remethylation of homocysteine to methionine. Mouse models of folate-responsive neural tube defects (NTDs) indicate that impaired de novo thymidylate (dTMP) synthesis through changes in SHMT expression is causative in folate-responsive NTDs. We have created a hybrid computational model comprised of ordinary differential equations and stochastic simulation. We investigated whether the de novo dTMP synthesis pathway was sensitive to perturbations in FOCM that are known to be associated with human NTDs. This computational model shows that de novo dTMP synthesis is highly sensitive to the common MTHFR C677T polymorphism and that the effect of the polymorphism on FOCM is greater in folate deficiency. Computational simulations indicate that the MTHFR C677T polymorphism and folate deficiency interact to increase the stochastic behavior of the FOCM network, with the greatest instability observed for reactions catalyzed by serine hydroxymethyltransferase (SHMT). Furthermore, we show that de novo dTMP synthesis does not occur in the cytosol at rates sufficient for DNA replication, supporting empirical data indicating that impaired nuclear de novo dTMP synthesis results in uracil misincorporation into DNA.

Graphical Modeling Meets Systems Pharmacology

A main source of failures in systems projects (including systems pharmacology) is poor communication level and different expectations among the stakeholders. A common and not ambiguous language that is naturally comprehensible by all the involved players is a boost to success. We present bStyle, a modeling tool that adopts a graphical language close enough to cartoons to be a common media to exchange ideas and data and that it is at the same time formal enough to enable modeling, analysis, and dynamic simulations of a system. Data analysis and simulation integrated in the same application are fundamental to understand the mechanisms of actions of drugs: a core aspect of systems pharmacology.

A Novel Insulin/Glucose Model after a Mixed-Meal Test in Patients with Type 1 Diabetes on Insulin Pump Therapy

Current closed-loop insulin delivery methods stem from sophisticated models of the glucose-insulin (G/I) system, mostly based on complex studies employing glucose tracer technology. We tested the performance of a new minimal model (GLUKINSLOOP 2.0) of the G/I system to characterize the glucose and insulin dynamics during multiple mixed meal tests (MMT) of different sizes in patients with type 1 diabetes (T1D) on insulin pump therapy (continuous subcutaneous insulin infusion, CSII). The GLUKINSLOOP 2.0 identified the G/I system, provided a close fit of the G/I time-courses and showed acceptable reproducibility of the G/I system parameters in repeated studies of identical and double-sized MMTs. This model can provide a fairly good and reproducible description of the G/I system in T1D patients on CSII, and it may be applied to create a bank of “virtual” patients. Our results might be relevant at improving the architecture of upcoming closed-loop CSII systems.

Simulation of biochemical reactions with time-dependent rates by the rejection-based algorithm

We address the problem of simulating biochemical reaction networks with time-dependent rates and propose a new algorithm based on our rejection-based stochastic simulation algorithm (RSSA) [Thanh et al., J. Chem. Phys. 141(13), 134116 (2014)]. The computation for selecting next reaction firings by our time-dependent RSSA (tRSSA) is computationally efficient. Furthermore, the generated trajectory is exact by exploiting the rejection-based mechanism. We benchmark tRSSA on different biological systems with varying forms of reaction rates to demonstrate its applicability and efficiency. We reveal that for nontrivial cases, the selection of reaction firings in existing algorithms introduces approximations because the integration of reaction rates is very computationally demanding and simplifying assumptions are introduced. The selection of the next reaction firing by our approach is easier while preserving the exactness.

On the rejection-based algorithm for simulation and analysis of large-scale reaction networks

Stochastic simulation for in silico studies of large biochemical networks requires a great amount of computational time. We recently proposed a new exact simulation algorithm, called the rejection-based stochastic simulation algorithm (RSSA) [Thanh et al., J. Chem. Phys. 141(13), 134116 (2014)], to improve simulation performance by postponing and collapsing as much as possible the propensity updates. In this paper, we analyze the performance of this algorithm in detail, and improve it for simulating large-scale biochemical reaction networks. We also present a new algorithm, called simultaneous RSSA (SRSSA), which generates many independent trajectories simultaneously for the analysis of the biochemical behavior. SRSSA improves simulation performance by utilizing a single data structure across simulations to select reaction firings and forming trajectories. The memory requirement for building and storing the data structure is thus independent of the number of trajectories. The updating of the data structure when needed is performed collectively in a single operation across the simulations. The trajectories generated by SRSSA are exact and independent of each other by exploiting the rejection-based mechanism. We test our new improvement on real biological systems with a wide range of reaction networks to demonstrate its applicability and efficiency.

Crowdsourced estimation of cognitive decline and resilience in Alzheimer's disease

Identifying accurate biomarkers of cognitive decline is essential for advancing early diagnosis and prevention therapies in Alzheimer's disease. The Alzheimer's disease DREAM Challenge was designed as a computational crowdsourced project to benchmark the current state-of-the-art in predicting cognitive outcomes in Alzheimer's disease based on high dimensional, publicly available genetic and structural imaging data. This meta-analysis failed to identify a meaningful predictor developed from either data modality, suggesting that alternate approaches should be considered for prediction of cognitive performance.

A network analysis of cofactor-protein interactions for analyzing associations between human nutrition and diseases

The involvement of vitamins and other micronutrients in intermediary metabolism was elucidated in the mid 1900's at the level of individual biochemical reactions. Biochemical pathways remain the foundational knowledgebase for understanding how micronutrient adequacy modulates health in all life stages. Current daily recommended intakes were usually established on the basis of the association of a single nutrient to a single, most sensitive adverse effect and thus neglect interdependent and pleiotropic effects of micronutrients on biological systems. Hence, the understanding of the impact of overt or sub-clinical nutrient deficiencies on biological processes remains incomplete. Developing a more complete view of the role of micronutrients and their metabolic products in protein-mediated reactions is of importance. We thus integrated and represented cofactor-protein interaction data from multiple and diverse sources into a multi-layer network representation that links cofactors, cofactor-interacting proteins, biological processes, and diseases. Network representation of this information is a key feature of the present analysis and enables the integration of data from individual biochemical reactions and protein-protein interactions into a systems view, which may guide strategies for targeted nutritional interventions aimed at improving health and preventing diseases.

Systems biology approaches to study the molecular effects of caloric restriction and polyphenols on aging processes

Worldwide population is aging, and a large part of the growing burden associated with age-related conditions can be prevented or delayed by promoting healthy lifestyle and normalizing metabolic risk factors. However, a better understanding of the pleiotropic effects of available nutritional interventions and their influence on the multiple processes affected by aging is needed to select and implement the most promising actions. New methods of analysis are required to tackle the complexity of the interplay between nutritional interventions and aging, and to make sense of a growing amount of -omics data being produced for this purpose. In this paper, we review how various systems biology-inspired methods of analysis can be applied to the study of the molecular basis of nutritional interventions promoting healthy aging, notably caloric restriction and polyphenol supplementation. We specifically focus on the role that different versions of network analysis, molecular signature identification and multi-omics data integration are playing in elucidating the complex mechanisms underlying nutrition, and provide some examples on how to extend the application of these methods using available microarray data.

Diversity of key players in the microbial ecosystems of the human body

Coexisting bacteria form various microbial communities in human body parts. In these ecosystems they interact in various ways and the properties of the interaction network can be related to the stability and functional diversity of the local bacterial community. In this study, we analyze the interaction network among bacterial OTUs in 11 locations of the human body. These belong to two major groups. One is the digestive system and the other is the female genital tract. In each local ecosystem we determine the key species, both the ones being in key positions in the interaction network and the ones that dominate by frequency. Beyond identifying the key players and discussing their biological relevance, we also quantify and compare the properties of the 11 networks. The interaction networks of the female genital system and the digestive system show totally different architecture. Both the topological properties and the identity of the key groups differ. Key groups represent four phyla of prokaryotes. Some groups appear in key positions in several locations, while others are assigned only to a single body part. The key groups of the digestive and the genital tracts are totally different.

Qualitative and Quantitative Protein Complex Prediction Through Proteome-Wide Simulations

Despite recent progress in proteomics most protein complexes are still unknown. Identification of these complexes will help us understand cellular regulatory mechanisms and support development of new drugs. Therefore it is really important to establish detailed information about the composition and the abundance of protein complexes but existing algorithms can only give qualitative predictions. Herein, we propose a new approach based on stochastic simulations of protein complex formation that integrates multi-source data--such as protein abundances, domain-domain interactions and functional annotations--to predict alternative forms of protein complexes together with their abundances. This method, called SiComPre (Simulation based Complex Prediction), achieves better qualitative prediction of yeast and human protein complexes than existing methods and is the first to predict protein complex abundances. Furthermore, we show that SiComPre can be used to predict complexome changes upon drug treatment with the example of bortezomib. SiComPre is the first method to produce quantitative predictions on the abundance of molecular complexes while performing the best qualitative predictions. With new data on tissue specific protein complexes becoming available SiComPre will be able to predict qualitative and quantitative differences in the complexome in various tissue types and under various conditions.

Efficient rejection-based simulation of biochemical reactions with stochastic noise and delays

We propose a new exact stochastic rejection-based simulation algorithm for biochemical reactions and extend it to systems with delays. Our algorithm accelerates the simulation by pre-computing reaction propensity bounds to select the next reaction to perform. Exploiting such bounds, we are able to avoid recomputing propensities every time a (delayed) reaction is initiated or finished, as is typically necessary in standard approaches. Propensity updates in our approach are still performed, but only infrequently and limited for a small number of reactions, saving computation time and without sacrificing exactness. We evaluate the performance improvement of our algorithm by experimenting with concrete biological models.

The role of breast-feeding in infant immune system: a systems perspective on the intestinal microbiome

BACKGROUND

The human intestinal microbiota changes from being sparsely populated and variable to possessing a mature, adult-like stable microbiome during the first 2 years of life. This assembly process of the microbiota can lead to either negative or positive effects on health, depending on the colonization sequence and diet. An integrative study on the diet, the microbiota, and genomic activity at the transcriptomic level may give an insight into the role of diet in shaping the human/microbiome relationship. This study aims at better understanding the effects of microbial community and feeding mode (breast-fed and formula-fed) on the immune system, by comparing intestinal metagenomic and transcriptomic data from breast-fed and formula-fed babies.

RESULTS

We re-analyzed a published metagenomics and host gene expression dataset from a systems biology perspective. Our results show that breast-fed samples co-express genes associated with immunological, metabolic, and biosynthetic activities. The diversity of the microbiota is higher in formula-fed than breast-fed infants, potentially reflecting the weaker dependence of infants on maternal microbiome. We mapped the microbial composition and the expression patterns for host systems and studied their relationship from a systems biology perspective, focusing on the differences.

CONCLUSIONS

Our findings revealed that there is co-expression of more genes in breast-fed samples but lower microbial diversity compared to formula-fed. Applying network-based systems biology approach via enrichment of microbial species with host genes revealed the novel key relationships of the microbiota with immune and metabolic activity. This was supported statistically by data and literature.

Novel drug target identification for the treatment of dementia using multi-relational association mining

Dementia is a neurodegenerative condition of the brain in which there is a progressive and permanent loss of cognitive and mental performance. Despite the fact that the number of people with dementia worldwide is steadily increasing and regardless of the advances in the molecular characterization of the disease, current medical treatments for dementia are purely symptomatic and hardly effective. We present a novel multi-relational association mining method that integrates the huge amount of scientific data accumulated in recent years to predict potential novel targets for innovative therapeutic treatment of dementia. Owing to the ability of processing large volumes of heterogeneous data, our method achieves a high performance and predicts numerous drug targets including several serine threonine kinase and a G-protein coupled receptor. The predicted drug targets are mainly functionally related to metabolism, cell surface receptor signaling pathways, immune response, apoptosis, and long-term memory. Among the highly represented kinase family and among the G-protein coupled receptors, DLG4 (PSD-95), and the bradikynin receptor 2 are highlighted also for their proposed role in memory and cognition, as described in previous studies. These novel putative targets hold promises for the development of novel therapeutic approaches for the treatment of dementia.

SCUDO: a tool for signature-based clustering of expression profiles

SCUDO (Signature-based ClUstering for DiagnOstic purposes) is an online tool for the analysis of gene expression profiles for diagnostic and classification purposes. The tool is based on a new method for the clustering of profiles based on a subject-specific, as opposed to disease-specific, signature. Our approach relies on construction of a reference map of transcriptional signatures, from both healthy and affected subjects, derived from their respective mRNA or miRNA profiles. A diagnosis for a new individual can then be performed by determining the position of the individual's transcriptional signature on the map. The diagnostic power of our method has been convincingly demonstrated in an open scientific competition (SBV Improver Diagnostic Signature Challenge), scoring second place overall and first place in one of the sub-challenges.

Genetic associations with micronutrient levels identified in immune and gastrointestinal networks

The discovery of vitamins and clarification of their role in preventing frank essential nutrient deficiencies occurred in the early 1900s. Much vitamin research has understandably focused on public health and the effects of single nutrients to alleviate acute conditions. The physiological processes for maintaining health, however, are complex systems that depend upon interactions between multiple nutrients, environmental factors, and genetic makeup. To analyze the relationship between these factors and nutritional health, data were obtained from an observational, community-based participatory research program of children and teens (age 6–14) enrolled in a summer day camp in the Delta region of Arkansas. Assessments of erythrocyte S-adenosylmethionine (SAM) and S-adenosylhomocysteine (SAH), plasma homocysteine (Hcy) and 6 organic micronutrients (retinol, 25-hydroxy vitamin D3, pyridoxal, thiamin, riboflavin, and vitamin E), and 1,129 plasma proteins were performed at 3 time points in each of 2 years. Genetic makeup was analyzed with 1 M SNP genotyping arrays, and nutrient status was assessed with 24-h dietary intake questionnaires. A pattern of metabolites (met_PC1) that included the ratio of erythrocyte SAM/SAH, Hcy, and 5 vitamins were identified by principal component analysis. Met_PC1 levels were significantly associated with (1) single-nucleotide polymorphisms, (2) levels of plasma proteins, and (3) multilocus genotypes coding for gastrointestinal and immune functions, as identified in a global network of metabolic/protein–protein interactions. Subsequent mining of data from curated pathway, network, and genome-wide association studies identified genetic and functional relationships that may be explained by gene–nutrient interactions. The systems nutrition strategy described here has thus associated a multivariate metabolite pattern in blood with genes involved in immune and gastrointestinal functions.

Methylation potential associated with diet, genotype, protein, and metabolite levels in the Delta Obesity Vitamin Study

Micronutrient research typically focuses on analyzing the effects of single or a few nutrients on health by analyzing a limited number of biomarkers. The observational study described here analyzed micronutrients, plasma proteins, dietary intakes, and genotype using a systems approach. Participants attended a community-based summer day program for 6-14 year old in 2 years. Genetic makeup, blood metabolite and protein levels, and dietary differences were measured in each individual. Twenty-four-hour dietary intakes, eight micronutrients (vitamins A, D, E, thiamin, folic acid, riboflavin, pyridoxal, and pyridoxine) and 3 one-carbon metabolites [homocysteine (Hcy), S-adenosylmethionine (SAM), and S-adenosylhomocysteine (SAH)], and 1,129 plasma proteins were analyzed as a function of diet at metabolite level, plasma protein level, age, and sex. Cluster analysis identified two groups differing in SAM/SAH and differing in dietary intake patterns indicating that SAM/SAH was a potential marker of nutritional status. The approach used to analyze genetic association with the SAM/SAH metabolites is called middle-out: SNPs in 275 genes involved in the one-carbon pathway (folate, pyridoxal/pyridoxine, thiamin) or were correlated with SAM/SAH (vitamin A, E, Hcy) were analyzed instead of the entire 1M SNP data set. This procedure identified 46 SNPs in 25 genes associated with SAM/SAH demonstrating a genetic contribution to the methylation potential. Individual plasma metabolites correlated with 99 plasma proteins. Fourteen proteins correlated with body mass index, 49 with group age, and 30 with sex. The analytical strategy described here identified subgroups for targeted nutritional interventions.

Consensus statement understanding health and malnutrition through a systems approach: the ENOUGH program for early life

Nutrition research, like most biomedical disciplines, adopted and often uses experimental approaches based on Beadle and Tatum's one gene-one polypeptide hypothesis, thereby reducing biological processes to single reactions or pathways. Systems thinking is needed to understand the complexity of health and disease processes requiring measurements of physiological processes, as well as environmental and social factors, which may alter the expression of genetic information. Analysis of physiological processes with omics technologies to assess systems' responses has only become available over the past decade and remains costly. Studies of environmental and social conditions known to alter health are often not connected to biomedical research. While these facts are widely accepted, developing and conducting comprehensive research programs for health are often beyond financial and human resources of single research groups. We propose a new research program on essential nutrients for optimal underpinning of growth and health (ENOUGH) that will use systems approaches with more comprehensive measurements and biostatistical analysis of the many biological and environmental factors that influence undernutrition. Creating a knowledge base for nutrition and health is a necessary first step toward developing solutions targeted to different populations in diverse social and physical environments for the two billion undernourished people in developed and developing economies.

Algorithmic modeling quantifies the complementary contribution of metabolic inhibitions to gemcitabine efficacy

Gemcitabine (2,2-difluorodeoxycytidine, dFdC) is a prodrug widely used for treating various carcinomas. Gemcitabine exerts its clinical effect by depleting the deoxyribonucleotide pools, and incorporating its triphosphate metabolite (dFdC-TP) into DNA, thereby inhibiting DNA synthesis. This process blocks the cell cycle in the early S phase, eventually resulting in apoptosis. The incorporation of gemcitabine into DNA takes place in competition with the natural nucleoside dCTP. The mechanisms of indirect competition between these cascades for common resources are given with the race for DNA incorporation; in clinical studies dedicated to singling out mechanisms of resistance, ribonucleotide reductase (RR) and deoxycytidine kinase (dCK) and human equilibrative nucleoside transporter1 (hENT1) have been associated to efficacy of gemcitabine with respect to their roles in the synthesis cascades of dFdC-TP and dCTP. However, the direct competition, which manifests itself in terms of inhibitions between these cascades, remains to be quantified. We propose an algorithmic model of gemcitabine mechanism of action, verified with respect to independent experimental data. We performed in silico experiments in different virtual conditions, otherwise difficult in vivo, to evaluate the contribution of the inhibitory mechanisms to gemcitabine efficacy. In agreement with the experimental data, our model indicates that the inhibitions due to the association of dCTP with dCK and the association of gemcitabine diphosphate metabolite (dFdC-DP) with RR play a key role in adjusting the efficacy. While the former tunes the catalysis of the rate-limiting first phosphorylation of dFdC, the latter is responsible for depletion of dCTP pools, thereby contributing to gemcitabine efficacy with a dependency on nucleoside transport efficiency. Our simulations predict the existence of a continuum of non-efficacy to high-efficacy regimes, where the levels of dFdC-TP and dCTP are coupled in a complementary manner, which can explain the resistance to this drug in some patients.

Biological network inference for drug discovery

A better understanding of the pathophysiology should help deliver drugs whose targets are involved in the causative processes underlying a disease. Biological network inference uses computational methods for deducing from high-throughput experimental data, the topology and the causal structure of the interactions among the drugs and their targets. Therefore, biological network inference can support and contribute to the experimental identification of both gene and protein networks causing a disease as well as the biochemical networks of drugs metabolism and mechanisms of action. The resulting high-level networks serve as a foundational basis for more detailed mechanistic models and are increasingly used in drug discovery by pharmaceutical and biotechnology companies. We review and compare recent computational technologies for network inference applied to drug discovery.

Inferring biochemical reaction pathways: the case of the gemcitabine pharmacokinetics

BACKGROUND

The representation of a biochemical system as a network is the precursor of any mathematical model of the processes driving the dynamics of that system. Pharmacokinetics uses mathematical models to describe the interactions between drug, and drug metabolites and targets and through the simulation of these models predicts drug levels and/or dynamic behaviors of drug entities in the body. Therefore, the development of computational techniques for inferring the interaction network of the drug entities and its kinetic parameters from observational data is raising great interest in the scientific community of pharmacologists. In fact, the network inference is a set of mathematical procedures deducing the structure of a model from the experimental data associated to the nodes of the network of interactions. In this paper, we deal with the inference of a pharmacokinetic network from the concentrations of the drug and its metabolites observed at discrete time points.

RESULTS

The method of network inference presented in this paper is inspired by the theory of time-lagged correlation inference with regard to the deduction of the interaction network, and on a maximum likelihood approach with regard to the estimation of the kinetic parameters of the network. Both network inference and parameter estimation have been designed specifically to identify systems of biotransformations, at the biochemical level, from noisy time-resolved experimental data. We use our inference method to deduce the metabolic pathway of the gemcitabine. The inputs to our inference algorithm are the experimental time series of the concentration of gemcitabine and its metabolites. The output is the set of reactions of the metabolic network of the gemcitabine.

CONCLUSIONS

Time-lagged correlation based inference pairs up to a probabilistic model of parameter inference from metabolites time series allows the identification of the microscopic pharmacokinetics and pharmacodynamics of a drug with a minimal a priori knowledge. In fact, the inference model presented in this paper is completely unsupervised. It takes as input the time series of the concetrations of the parent drug and its metabolites. The method, applied to the case study of the gemcitabine pharmacokinetics, shows good accuracy and sensitivity.

Model-based clustering reveals vitamin D dependent multi-centrality hubs in a network of vitamin-related proteins

Background

Nutritional systems biology offers the potential for comprehensive predictions that account for all metabolic changes with the intricate biological organization and the multitudinous interactions between the cellular proteins. Protein-protein interaction (PPI) networks can be used for an integrative description of molecular processes. Although widely adopted in nutritional systems biology, these networks typically encompass a single category of functional interaction (i.e., metabolic, regulatory or signaling) or nutrient. Incorporating multiple nutrients and functional interaction categories under an integrated framework represents an informative approach for gaining system level insight on nutrient metabolism.

Results

We constructed a multi-level PPI network starting from the interactions of 200 vitamin-related proteins. Its final size was 1,657 proteins, with 2,700 interactions. To characterize the role of the proteins we computed 6 centrality indices and applied model-based clustering. We detected a subgroup of 22 proteins that were highly central and significantly related to vitamin D. Immune system and cancer-related processes were strongly represented among these proteins. Clustering of the centralities revealed a degree of redundancy among the indices; a repeated analysis using subsets of the centralities performed well in identifying the original set of 22 most central proteins.

Conclusions

Hierarchical and model-based clustering revealed multi-centrality hubs in a vitamin PPI network and redundancies among the centrality indices. Vitamin D-related proteins were strongly represented among network hubs, highlighting the pervasive effects of this nutrient. Our integrated approach to network construction identified promiscuous transcription factors, cytokines and enzymes - primarily related to immune system and cancer processes - representing potential gatekeepers linking vitamin intake to disease.

Stochastic simulation of the spatio-temporal dynamics of reaction-diffusion systems: the case for the bicoid gradient

Reaction-diffusion systems are mathematical models that describe how the concentrations of substances distributed in space change under the influence of local chemical reactions, and diffusion which causes the substances to spread out in space. The classical representation of a reaction-diffusion system is given by semi-linear parabolic partial differential equations, whose solution predicts how diffusion causes the concentration field to change with time. This change is proportional to the diffusion coefficient. If the solute moves in a homogeneous system in thermal equilibrium, the diffusion coefficients are constants that do not depend on the local concentration of solvent and solute. However, in nonhomogeneous and structured media the assumption of constant intracellular diffusion coefficient is not necessarily valid, and, consequently, the diffusion coefficient is a function of the local concentration of solvent and solutes. In this paper we propose a stochastic model of reaction-diffusion systems, in which the diffusion coefficients are function of the local concentration, viscosity and frictional forces. We then describe the software tool Redi (REaction-DIffusion simulator) which we have developed in order to implement this model into a Gillespie-like stochastic simulation algorithm. Finally, we show the ability of our model implemented in the Redi tool to reproduce the observed gradient of the bicoid protein in the Drosophila Melanogaster embryo. With Redi, we were able to simulate with an accuracy of 1% the experimental spatio-temporal dynamics of the bicoid protein, as recorded in time-lapse experiments obtained by direct measurements of transgenic bicoidenhanced green fluorescent protein.

Snazer: the simulations and networks analyzer

BACKGROUND

Networks are widely recognized as key determinants of structure and function in systems that span the biological, physical, and social sciences. They are static pictures of the interactions among the components of complex systems. Often, much effort is required to identify networks as part of particular patterns as well as to visualize and interpret them.From a pure dynamical perspective, simulation represents a relevant way-out. Many simulator tools capitalized on the "noisy" behavior of some systems and used formal models to represent cellular activities as temporal trajectories. Statistical methods have been applied to a fairly large number of replicated trajectories in order to infer knowledge.A tool which both graphically manipulates reactive models and deals with sets of simulation time-course data by aggregation, interpretation and statistical analysis is missing and could add value to simulators.

RESULTS

We designed and implemented Snazer, the simulations and networks analyzer. Its goal is to aid the processes of visualizing and manipulating reactive models, as well as to share and interpret time-course data produced by stochastic simulators or by any other means.

CONCLUSIONS

Snazer is a solid prototype that integrates biological network and simulation time-course data analysis techniques.

Elucidation of functional consequences of signalling pathway interactions

BACKGROUND

A great deal of data has accumulated on signalling pathways. These large datasets are thought to contain much implicit information on their molecular structure, interaction and activity information, which provides a picture of intricate molecular networks believed to underlie biological functions. While tremendous advances have been made in trying to understand these systems, how information is transmitted within them is still poorly understood. This ever growing amount of data demands we adopt powerful computational techniques that will play a pivotal role in the conversion of mined data to knowledge, and in elucidating the topological and functional properties of protein - protein interactions.

RESULTS

A computational framework is presented which allows for the description of embedded networks, and identification of common shared components thought to assist in the transmission of information within the systems studied. By employing the graph theories of network biology - such as degree distribution, clustering coefficient, vertex betweenness and shortest path measures - topological features of protein-protein interactions for published datasets of the p53, nuclear factor kappa B (NF-kappaB) and G1/S phase of the cell cycle systems were ascertained. Highly ranked nodes which in some cases were identified as connecting proteins most likely responsible for propagation of transduction signals across the networks were determined. The functional consequences of these nodes in the context of their network environment were also determined. These findings highlight the usefulness of the framework in identifying possible combination or links as targets for therapeutic responses; and put forward the idea of using retrieved knowledge on the shared components in constructing better organised and structured models of signalling networks.

CONCLUSION

It is hoped that through the data mined reconstructed signal transduction networks, well developed models of the published data can be built which in the end would guide the prediction of new targets based on the pathway's environment for further analysis. Source code is available upon request.