PathVisio-MIM: PathVisio plugin for creating and editing Molecular Interaction Maps (MIMs)

PG-path: Modeling and personalizing pharmacogenomics-based pathways

A pharmacogenomics-based pathway represents a series of reactions that occur between drugs and genes in the human body after drug administration. PG-path is a pharmacogenomics-based pathway that standardizes and visualizes the components (nodes) and actions (edges) involved in pharmacokinetic and pharmacodynamic processes. It provides an intuitive understanding of the drug response in the human body. A pharmacokinetic pathway visualizes the absorption, distribution, metabolism, and excretion (ADME) at the systemic level, and a pharmacodynamic pathway shows the action of the drug in the target cell at the cellular-molecular level. The genes in the pathway are displayed in locations similar to those inside the body. PG-path allows personalized pathways to be created by annotating each gene with the overall impact degree of deleterious variants in the gene. These personalized pathways play a role in assisting tailored individual prescriptions by predicting changes in the drug concentration in the plasma. PG-path also supports counseling for personalized drug therapy by providing visualization and documentation.

Using the Semantic Web for Rapid Integration of WikiPathways with Other Biological Online Data Resources

The diversity of online resources storing biological data in different formats provides a challenge for bioinformaticians to integrate and analyse their biological data. The semantic web provides a standard to facilitate knowledge integration using statements built as triples describing a relation between two objects. WikiPathways, an online collaborative pathway resource, is now available in the semantic web through a SPARQL endpoint at http://sparql.wikipathways.org. Having biological pathways in the semantic web allows rapid integration with data from other resources that contain information about elements present in pathways using SPARQL queries. In order to convert WikiPathways content into meaningful triples we developed two new vocabularies that capture the graphical representation and the pathway logic, respectively. Each gene, protein, and metabolite in a given pathway is defined with a standard set of identifiers to support linking to several other biological resources in the semantic web. WikiPathways triples were loaded into the Open PHACTS discovery platform and are available through its Web API (https://dev.openphacts.org/docs) to be used in various tools for drug development. We combined various semantic web resources with the newly converted WikiPathways content using a variety of SPARQL query types and third-party resources, such as the Open PHACTS API. The ability to use pathway information to form new links across diverse biological data highlights the utility of integrating WikiPathways in the semantic web.

Automatically visualise and analyse data on pathways using PathVisioRPC from any programming environment

BACKGROUND

Biological pathways are descriptive diagrams of biological processes widely used for functional analysis of differentially expressed genes or proteins. Primary data analysis, such as quality control, normalisation, and statistical analysis, is often performed in scripting languages like R, Perl, and Python. Subsequent pathway analysis is usually performed using dedicated external applications. Workflows involving manual use of multiple environments are time consuming and error prone. Therefore, tools are needed that enable pathway analysis directly within the same scripting languages used for primary data analyses. Existing tools have limited capability in terms of available pathway content, pathway editing and visualisation options, and export file formats. Consequently, making the full-fledged pathway analysis tool PathVisio available from various scripting languages will benefit researchers.

RESULTS

We developed PathVisioRPC, an XMLRPC interface for the pathway analysis software PathVisio. PathVisioRPC enables creating and editing biological pathways, visualising data on pathways, performing pathway statistics, and exporting results in several image formats in multiple programming environments. We demonstrate PathVisioRPC functionalities using examples in Python. Subsequently, we analyse a publicly available NCBI GEO gene expression dataset studying tumour bearing mice treated with cyclophosphamide in R. The R scripts demonstrate how calls to existing R packages for data processing and calls to PathVisioRPC can directly work together. To further support R users, we have created RPathVisio simplifying the use of PathVisioRPC in this environment. We have also created a pathway module for the microarray data analysis portal ArrayAnalysis.org that calls the PathVisioRPC interface to perform pathway analysis. This module allows users to use PathVisio functionality online without having to download and install the software and exemplifies how the PathVisioRPC interface can be used by data analysis pipelines for functional analysis of processed genomics data.

CONCLUSIONS

PathVisioRPC enables data visualisation and pathway analysis directly from within various analytical environments used for preliminary analyses. It supports the use of existing pathways from WikiPathways or pathways created using the RPC itself. It also enables automation of tasks performed using PathVisio, making it useful to PathVisio users performing repeated visualisation and analysis tasks. PathVisioRPC is freely available for academic and commercial use at http://projects.bigcat.unimaas.nl/pathvisiorpc.

Predicted Role of NAD Utilization in the Control of Circadian Rhythms during DNA Damage Response

The circadian clock is a set of regulatory steps that oscillate with a period of approximately 24 hours influencing many biological processes. These oscillations are robust to external stresses, and in the case of genotoxic stress (i.e. DNA damage), the circadian clock responds through phase shifting with primarily phase advancements. The effect of DNA damage on the circadian clock and the mechanism through which this effect operates remains to be thoroughly investigated. Here we build an in silico model to examine damage-induced circadian phase shifts by investigating a possible mechanism linking circadian rhythms to metabolism. The proposed model involves two DNA damage response proteins, SIRT1 and PARP1, that are each consumers of nicotinamide adenine dinucleotide (NAD), a metabolite involved in oxidation-reduction reactions and in ATP synthesis. This model builds on two key findings: 1) that SIRT1 (a protein deacetylase) is involved in both the positive (i.e. transcriptional activation) and negative (i.e. transcriptional repression) arms of the circadian regulation and 2) that PARP1 is a major consumer of NAD during the DNA damage response. In our simulations, we observe that increased PARP1 activity may be able to trigger SIRT1-induced circadian phase advancements by decreasing SIRT1 activity through competition for NAD supplies. We show how this competitive inhibition may operate through protein acetylation in conjunction with phosphorylation, consistent with reported observations. These findings suggest a possible mechanism through which multiple perturbations, each dominant during different points of the circadian cycle, may result in the phase advancement of the circadian clock seen during DNA damage.

Many physiological processes are regulated by the circadian clock, and we are continuing to learn about the role of the circadian clock in disease. Research in recent years has begun to shed light on the feedback mechanisms that exist between circadian regulation and other processes, including metabolism and the response to DNA damage. A challenge has been to understand the dynamic nature of the protein interactions of these processes, which often involve protein modification as a means of communicating cellular states, such as damaged DNA. Here we have devised a model that simulates an alteration of the circadian clock that is observed during DNA damage response. A novel aspect of this model is the inclusion of SIRT1, a protein that regulates core circadian proteins through modification and helps to repress gene expression. SIRT1 is dependent on a metabolite regulated by the circadian clock and is depleted during DNA damage. In conjunction with a second form of protein modification, our results suggest that multiple forms of protein modification may contribute to the experimentally observed alterations to circadian function.

PathVisio 3: an extendable pathway analysis toolbox

PathVisio is a commonly used pathway editor, visualization and analysis software. Biological pathways have been used by biologists for many years to describe the detailed steps in biological processes. Those powerful, visual representations help researchers to better understand, share and discuss knowledge. Since the first publication of PathVisio in 2008, the original paper was cited more than 170 times and PathVisio was used in many different biological studies. As an online editor PathVisio is also integrated in the community curated pathway database WikiPathways. Here we present the third version of PathVisio with the newest additions and improvements of the application. The core features of PathVisio are pathway drawing, advanced data visualization and pathway statistics. Additionally, PathVisio 3 introduces a new powerful extension systems that allows other developers to contribute additional functionality in form of plugins without changing the core application. PathVisio can be downloaded from http://www.pathvisio.org and in 2014 PathVisio 3 has been downloaded over 5,500 times. There are already more than 15 plugins available in the central plugin repository. PathVisio is a freely available, open-source tool published under the Apache 2.0 license (http://www.apache.org/licenses/LICENSE-2.0). It is implemented in Java and thus runs on all major operating systems. The code repository is available at http://svn.bigcat.unimaas.nl/pathvisio. The support mailing list for users is available on https://groups.google.com/forum/#!forum/wikipathways-discuss and for developers on https://groups.google.com/forum/#!forum/wikipathways-devel.

Gene expression correlations in human cancer cell lines define molecular interaction networks for epithelial phenotype

Using gene expression data to enhance our knowledge of control networks relevant to cancer biology and therapy is a challenging but urgent task. Based on the premise that genes that are expressed together in a variety of cell types are likely to functions together, we derived mutually correlated genes that function together in various processes in epithelial-like tumor cells. Expression-correlated genes were derived from data for the NCI-60 human tumor cell lines, as well as data from the Broad Institute's CCLE cell lines. NCI-60 cell lines that selectively expressed a mutually correlated subset of tight junction genes served as a signature for epithelial-like cancer cells. Those signature cell lines served as a seed to derive other correlated genes, many of which had various other epithelial-related functions. Literature survey yielded molecular interaction and function information about those genes, from which molecular interaction maps were assembled. Many of the genes had epithelial functions unrelated to tight junctions, demonstrating that new function categories were elicited. The most highly correlated genes were implicated in the following epithelial functions: interactions at tight junctions (CLDN7, CLDN4, CLDN3, MARVELD3, MARVELD2, TJP3, CGN, CRB3, LLGL2, EPCAM, LNX1); interactions at adherens junctions (CDH1, ADAP1, CAMSAP3); interactions at desmosomes (PPL, PKP3, JUP); transcription regulation of cell-cell junction complexes (GRHL1 and 2); epithelial RNA splicing regulators (ESRP1 and 2); epithelial vesicle traffic (RAB25, EPN3, GRHL2, EHF, ADAP1, MYO5B); epithelial Ca(+2) signaling (ATP2C2, S100A14, BSPRY); terminal differentiation of epithelial cells (OVOL1 and 2, ST14, PRSS8, SPINT1 and 2); maintenance of apico-basal polarity (RAB25, LLGL2, EPN3). The findings provide a foundation for future studies to elucidate the functions of regulatory networks specific to epithelial-like cancer cells and to probe for anti-cancer drug targets.

SIRT1/PARP1 crosstalk: connecting DNA damage and metabolism

An intricate network regulates the activities of SIRT1 and PARP1 proteins and continues to be uncovered. Both SIRT1 and PARP1 share a common co-factor nicotinamide adenine dinucleotide (NAD+) and several common substrates, including regulators of DNA damage response and circadian rhythms. We review this complex network using an interactive Molecular Interaction Map (MIM) to explore the interplay between these two proteins. Here we discuss how NAD + competition and post-transcriptional/translational feedback mechanisms create a regulatory network sensitive to environmental cues, such as genotoxic stress and metabolic states, and examine the role of those interactions in DNA repair and ultimately, cell fate decisions.

Gene expression profiles of the NCI-60 human tumor cell lines define molecular interaction networks governing cell migration processes

Although there is extensive information on gene expression and molecular interactions in various cell types, integrating those data in a functionally coherent manner remains challenging. This study explores the premise that genes whose expression at the mRNA level is correlated over diverse cell lines are likely to function together in a network of molecular interactions. We previously derived expression-correlated gene clusters from the database of the NCI-60 human tumor cell lines and associated each cluster with function categories of the Gene Ontology (GO) database. From a cluster rich in genes associated with GO categories related to cell migration, we extracted 15 genes that were highly cross-correlated; prominent among them were RRAS, AXL, ADAM9, FN14, and integrin-beta1. We then used those 15 genes as bait to identify other correlated genes in the NCI-60 database. A survey of current literature disclosed, not only that many of the expression-correlated genes engaged in molecular interactions related to migration, invasion, and metastasis, but that highly cross-correlated subsets of those genes engaged in specific cell migration processes. We assembled this information in molecular interaction maps (MIMs) that depict networks governing 3 cell migration processes: degradation of extracellular matrix, production of transient focal complexes at the leading edge of the cell, and retraction of the rear part of the cell. Also depicted are interactions controlling the release and effects of calcium ions, which may regulate migration in a spaciotemporal manner in the cell. The MIMs and associated text comprise a detailed and integrated summary of what is currently known or surmised about the role of the expression cross-correlated genes in molecular networks governing those processes.

PathVisio-Validator: a rule-based validation plugin for graphical pathway notations

PURPOSE

The PathVisio-Validator plugin aims to simplify the task of producing biological pathway diagrams that follow graphical standardized notations, such as Molecular Interaction Maps or the Systems Biology Graphical Notation. This plugin assists in the creation of pathway diagrams by ensuring correct usage of a notation, and thereby reducing ambiguity when diagrams are shared among biologists. Rulesets, needed in the validation process, can be generated for any graphical notation that a developer desires, using either Schematron or Groovy. The plugin also provides support for filtering validation results, validating on a subset of rules, and distinguishing errors and warnings.

A formal MIM specification and tools for the common exchange of MIM diagrams: an XML-Based format, an API, and a validation method

Background

The Molecular Interaction Map (MIM) notation offers a standard set of symbols and rules on their usage for the depiction of cellular signaling network diagrams. Such diagrams are essential for disseminating biological information in a concise manner. A lack of software tools for the notation restricts wider usage of the notation. Development of software is facilitated by a more detailed specification regarding software requirements than has previously existed for the MIM notation.

Results

A formal implementation of the MIM notation was developed based on a core set of previously defined glyphs. This implementation provides a detailed specification of the properties of the elements of the MIM notation. Building upon this specification, a machine-readable format is provided as a standardized mechanism for the storage and exchange of MIM diagrams. This new format is accompanied by a Java-based application programming interface to help software developers to integrate MIM support into software projects. A validation mechanism is also provided to determine whether MIM datasets are in accordance with syntax rules provided by the new specification.

Conclusions

The work presented here provides key foundational components to promote software development for the MIM notation. These components will speed up the development of interoperable tools supporting the MIM notation and will aid in the translation of data stored in MIM diagrams to other standardized formats. Several projects utilizing this implementation of the notation are outlined herein. The MIM specification is available as an additional file to this publication. Source code, libraries, documentation, and examples are available at http://discover.nci.nih.gov/mim.