Jin-Mei Lai

PhosphoPOINT: a comprehensive human kinase interactome and phospho-protein database

MOTIVATION To fully understand how a protein kinase regulates biological processes, it is imperative to first identify its substrate(s) and interacting protein(s). However, of the 518 known human serine/threonine/tyrosine kinases, 35% of these have known substrates, while 14% of the kinases have identified substrate recognition motifs. In contrast, 85% of the kinases have protein-protein interaction (PPI) datasets, raising the possibility that we might reveal potential kinase-substrate pairs from these PPIs. RESULTS PhosphoPOINT, a comprehensive human kinase interactome and phospho-protein database, is a collection of 4195 phospho-proteins with a total of 15 738 phosphorylation sites. PhosphoPOINT annotates the interactions among kinases, with their down-stream substrates and with interacting (phospho)-proteins to modulate the kinase-substrate pairs. PhosphoPOINT implements various gene expression profiles and Gene Ontology cellular component information to evaluate each kinase and their interacting (phospho)-proteins/substrates. Integration of cSNPs that cause amino acids change with the proteins with the phosphoprotein dataset reveals that 64 phosphorylation sites result in a disease phenotypes when changed; the linked phenotypes include schizophrenia and hypertension. PhosphoPOINT also provides a search function for all phospho-peptides using about 300 known kinase/phosphatase substrate/binding motifs. Altogether, PhosphoPOINT provides robust annotation for kinases, their downstream substrates and their interaction (phospho)-proteins and this should accelerate the functional characterization of kinomemediated signaling. AVAILABILITY PhosphoPOINT can be freely accessed in http://kinase. bioinformatics.tw/. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

Ortholog-based protein-protein interaction prediction and its application to inter-species interactions

BackgroundThe rapid growth of protein-protein interaction (PPI) data has led to the emergence of PPI network analysis. Despite advances in high-throughput techniques, the interactomes of several model organisms are still far from complete. Therefore, it is desirable to expand these interactomes with ortholog-based and other methods.ResultsOrthologous pairs of 18 eukaryotic species were expanded and merged with experimental PPI datasets. The contributions of interologs from each species were evaluated. The expanded orthologous pairs enable the inference of interologs for various species. For example, more than 32,000 human interactions can be predicted. The same dataset has also been applied to the prediction of host-pathogen interactions. PPIs between P. falciparum calmodulin and several H. sapiens proteins are predicted, and these interactions may contribute to the maintenance of host cell Ca2+ concentration. Using comparisons with Bayesian and structure-based approaches, interactions between putative HSP40 homologs of P. falciparum and the H. sapiens TNF receptor associated factor family are revealed, suggesting a role for these interactions in the interference of the human immune response to P. falciparum.ConclusionThe PPI datasets are available from POINT http://point.bioinformatics.tw/ and POINeT http://poinet.bioinformatics.tw/. Further development of methods to predict host-pathogen interactions should incorporate multiple approaches in order to improve sensitivity, and should facilitate the identification of targets for drug discovery and design.

Detection of the inferred interaction network in hepatocellular carcinoma from EHCO (Encyclopedia of Hepatocellular Carcinoma genes Online)

BackgroundThe significant advances in microarray and proteomics analyses have resulted in an exponential increase in potential new targets and have promised to shed light on the identification of disease markers and cellular pathways. We aim to collect and decipher the HCC-related genes at the systems level.ResultsHere, we build an integrative platform, the E ncyclopedia of H epatocellular C arcinoma genes O nline, dubbed EHCO http://ehco.iis.sinica.edu.tw, to systematically collect, organize and compare the pileup of unsorted HCC-related studies by using natural language processing and softbots. Among the eight gene set collections, ranging across PubMed, SAGE, microarray, and proteomics data, there are 2,906 genes in total; however, more than 77% genes are only included once, suggesting that tremendous efforts need to be exerted to characterize the relationship between HCC and these genes. Of these HCC inventories, protein binding represents the largest proportion (~25%) from Gene Ontology analysis. In fact, many differentially expressed gene sets in EHCO could form interaction networks (e.g. HBV-associated HCC network) by using available human protein-protein interaction datasets. To further highlight the potential new targets in the inferred network from EHCO, we combine comparative genomics and interactomics approaches to analyze 120 evolutionary conserved and overexpressed genes in HCC. 47 out of 120 queries can form a highly interactive network with 18 queries serving as hubs.ConclusionThis architectural map may represent the first step toward the attempt to decipher the hepatocarcinogenesis at the systems level. Targeting hubs and/or disruption of the network formation might reveal novel strategy for HCC treatment.

POINeT: protein interactome with sub-network analysis and hub prioritization

BackgroundProtein-protein interactions (PPIs) are critical to every aspect of biological processes. Expansion of all PPIs from a set of given queries often results in a complex PPI network lacking spatiotemporal consideration. Moreover, the reliability of available PPI resources, which consist of low- and high-throughput data, for network construction remains a significant challenge. Even though a number of software tools are available to facilitate PPI network analysis, an integrated tool is crucial to alleviate the burden on querying across multiple web servers and software tools.ResultsWe have constructed an integrated web service, POINeT, to simplify the process of PPI searching, analysis, and visualization. POINeT merges PPI and tissue-specific expression data from multiple resources. The tissue-specific PPIs and the numbers of research papers supporting the PPIs can be filtered with user-adjustable threshold values and are dynamically updated in the viewer. The network constructed in POINeT can be readily analyzed with, for example, the built-in centrality calculation module and an integrated network viewer. Nodes in global networks can also be ranked and filtered using various network analysis formulas, i.e., centralities. To prioritize the sub-network, we developed a ranking filtered method (S3) to uncover potential novel mediators in the midbody network. Several examples are provided to illustrate the functionality of POINeT. The network constructed from four schizophrenia risk markers suggests that EXOC4 might be a novel marker for this disease. Finally, a liver-specific PPI network has been filtered with adult and fetal liver expression profiles.ConclusionThe functionalities provided by POINeT are highly improved compared to previous version of POINT. POINeT enables the identification and ranking of potential novel genes involved in a sub-network. Combining with tissue-specific gene expression profiles, PPIs specific to selected tissues can be revealed. The straightforward interface of POINeT makes PPI search and analysis just a few clicks away. The modular design permits further functional enhancement without hampering the simplicity. POINeT is available at http://poinet.bioinformatics.tw/.

Gene expression profiles of the aurora family kinases.

The evolutionarily conserved Aurora family kinases, a family of mitotic serine/threonine kinases, has three members in humans (Aurora-A, -B and -C). Overexpression of Aurora family members, particularly Aurora-A, has been reported in many human cancers and cell lines. In this study, we present evidence based on comparative gene expression analysis via quantitative RT-PCR to delineate the relative contributions of these kinases in 60 cell lines and statistical analysis in five different human cancer microarray datasets. The analysis demonstrated the selective upregulation of these Aurora members in various cancers. In general, Aurora-A exhibited the highest expression levels, with substantially decreased quantities of the Aurora-C transcript detected relative to Aurora-A and -B. Moreover, to characterize the roles of each Aurora member, which share many similarities, we investigated the expression profiles of the family in normal tissues and a panel of different phases of the HeLa cell cycle. Finally, both Aurora-A and -B were overexpressed in a majority of esophageal tumor tissues in comparison to the normal variants. Taken together, the results show that each Aurora member exhibits distinct expression patterns, implying that they are engaged in different biological processes to accomplish more elaborate cell physiological functions in higher organisms.

Gene Expression-Based Chemical Genomics Identifies Potential Therapeutic Drugs in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is an aggressive tumor with a poor prognosis. Currently, only sorafenib is approved by the FDA for advanced HCC treatment; therefore, there is an urgent need to discover candidate therapeutic drugs for HCC. We hypothesized that if a drug signature could reverse, at least in part, the gene expression signature of HCC, it might have the potential to inhibit HCC-related pathways and thereby treat HCC. To test this hypothesis, we first built an integrative platform, the “Encyclopedia of Hepatocellular Carcinoma genes Online 2”, dubbed EHCO2, to systematically collect, organize and compare the publicly available data from HCC studies. The resulting collection includes a total of 4,020 genes. To systematically query the Connectivity Map (CMap), which includes 6,100 drug-mediated expression profiles, we further designed various gene signature selection and enrichment methods, including a randomization technique, majority vote, and clique analysis. Subsequently, 28 out of 50 prioritized drugs, including tanespimycin, trichostatin A, thioguanosine, and several anti-psychotic drugs with anti-tumor activities, were validated via MTT cell viability assays and clonogenic assays in HCC cell lines. To accelerate their future clinical use, possibly through drug-repurposing, we selected two well-established drugs to test in mice, chlorpromazine and trifluoperazine. Both drugs inhibited orthotopic liver tumor growth. In conclusion, we successfully discovered and validated existing drugs for potential HCC therapeutic use with the pipeline of Connectivity Map analysis and lab verification, thereby suggesting the usefulness of this procedure to accelerate drug repurposing for HCC treatment.

Graptopetalum Paraguayense Ameliorates Chemical-Induced Rat Hepatic Fibrosis In Vivo and Inactivates Stellate Cells and Kupffer Cells In Vitro

Background Graptopetalum paraguayense (GP) is a folk herbal medicine with hepatoprotective effects that is used in Taiwan. The aim of this study was to evaluate the hepatoprotective and antifibrotic effects of GP on experimental hepatic fibrosis in both dimethylnitrosamine (DMN)- and carbon tetrachloride (CCl4)-induced liver injury rats. Methods Hepatic fibrosis-induced rats were fed with the methanolic extract of GP (MGP) by oral administration every day. Immunohistochemistry, biochemical assays, and Western blot analysis were performed. The effects of MGP on the expression of fibrotic markers and cytokines in the primary cultured hepatic stellate cells (HSCs) and Kupffer cells, respectively, were evaluated. Results Oral administration of MGP significantly alleviated DMN- or CCl4-induced liver inflammation and fibrosis. High levels of alanine transaminase, aspartate transaminase, bilirubin, prothrombin activity and mortality rates also decreased in rats treated with MGP. There were significantly decreased hydroxyproline levels in therapeutic rats compared with those of the liver-damaged rats. Collagen I and alpha smooth muscle actin (α-SMA) expression were all reduced by incubation with MGP in primary cultured rat HSCs. Furthermore, MGP induced apoptotic cell death in activated HSCs. MGP also suppressed lipopolysaccharide-stimulated rat Kupffer cell activation by decreasing nitric oxide, tumor necrosis factor-α and interleukin-6 production, and increasing interleukin-10 expression. Conclusions The results show that the administration of MGP attenuated toxin-induced hepatic damage and fibrosis in vivo and inhibited HSC and Kupffer cell activation in vitro, suggesting that MGP might be a promising complementary or alternative therapeutic agent for liver inflammation and fibrosis.

Cliques in mitotic spindle network bring kinetochore-associated complexes to form dependence pathway.

The mitotic spindle is an essential molecular machine for chromosome segregation during mitosis. Achieving a better understanding of its organization at the topological level remains a daunting task. To determine the functional connections among 137 mitotic spindle proteins, a protein–protein interaction network among queries was constructed. Many hub proteins, which connect more than one query and serve as highly plausible candidates for expanding the mitotic spindle proteome, are ranked by conventional degree centrality and a new subnetwork specificity score. Evaluation of the ranking results by literature reviews and empirical verification of SEPT6, a novel top‐ranked hub, suggests that the subnetwork specificity score could enrich for putative spindle‐related proteins. Topological analysis of this expanded network shows the presence of 30 3‐cliques and six 4‐cliques (fully connected subgraphs) that, respectively, reside in eight kinetochore‐associated complexes, of which seven are evolution conserved. Notably, these complexes strikingly form dependence pathways for the assembly of the kinetochore complex. These analyses indicate the feasibility of using network topology, i.e. cliques, to uncover novel pathways to accelerate our understanding of potential biological processes.

Using an in Situ Proximity Ligation Assay to Systematically Profile Endogenous Protein–Protein Interactions in a Pathway Network

Signal transduction pathways in the cell require protein-protein interactions (PPIs) to respond to environmental cues. Diverse experimental techniques for detecting PPIs have been developed. However, the huge amount of PPI data accumulated from various sources poses a challenge with respect to data reliability. Herein, we collected ∼ 700 primary antibodies and employed a highly sensitive and specific technique, an in situ proximity ligation assay, to investigate 1204 endogenous PPIs in HeLa cells, and 557 PPIs of them tested positive. To overview the tested PPIs, we mapped them into 13 PPI public databases, which showed 72% of them were annotated in the Human Protein Reference Database (HPRD) and 8 PPIs were new PPIs not in the PubMed database. Moreover, TP53, CTNNB1, AKT1, CDKN1A, and CASP3 were the top 5 proteins prioritized by topology analyses of the 557 PPI network. Integration of the PPI-pathway interaction revealed that 90 PPIs were cross-talk PPIs linking 17 signaling pathways based on Reactome annotations. The top 2 connected cross-talk PPIs are MAPK3-DAPK1 and FAS-PRKCA interactions, which link 9 and 8 pathways, respectively. In summary, we established an open resource for biological modules and signaling pathway profiles, providing a foundation for comprehensive analysis of the human interactome.

Inferring a transcriptional regulatory network of the cytokinesis-related genes by network component analysis

BackgroundNetwork Component Analysis (NCA) is a network structure-driven framework for deducing regulatory signal dynamics. In contrast to principal component analysis, which can be employed to select the high-variance genes, NCA makes use of the connectivity structure from transcriptional regulatory networks to infer dynamics of transcription factor activities. Using the budding yeast Saccharomyces cerevisiae as a model system, we aim to deduce regulatory actions of cytokinesis-related genes, using precise spatial proximity (midbody) and/or temporal synchronicity (cytokinesis) to avoid full-scale computation from genome-wide databases.ResultsNCA was applied to infer regulatory actions of transcription factor activity from microarray data and partial transcription factor-gene connectivity information for cytokinesis-related genes, which were a subset of genome-wide datasets. No literature has so far discussed the inferred results through NCA are independent of the scale of the gene expression dataset. To avoid full-scale computation from genome-wide databases, four cytokinesis-related gene cases were selected for NCA by running computational analysis over the transcription factor database to confirm the approach being scale-free. The inferred dynamics of transcription factor activity through NCA were independent of the scale of the data matrix selected from the four cytokinesis-related gene sets. Moreover, the inferred regulatory actions were nearly identical to published observations for the selected cytokinesis-related genes in the budding yeast; namely, Mcm1, Ndd1, and Fkh2, which form a transcription factor complex to control expression of the CLB2 cluster (i.e. BUD4, CHS2, IQG1, and CDC5).ConclusionIn this study, using S. cerevisiae as a model system, NCA was successfully applied to infer similar regulatory actions of transcription factor activities from two various microarray databases and several partial transcription factor-gene connectivity datasets for selected cytokinesis-related genes independent of data sizes. The regulated action for four selected cytokinesis-related genes (BUD4, CHS2, IQG1, and CDC5) belongs to the M-phase or M/G1 phase, consistent with the empirical observations that in S. cerevisiae, the Mcm1-Ndd1-Fkh2 transcription factor complex can regulate expression of the cytokinesis-related genes BUD4, CHS2, IQG1, and CDC5. Since Bud4, Iqg1, and Cdc5 are highly conserved between human and yeast, results obtained from NCA for cytokinesis in the budding yeast can lead to a suggestion that human cells should have the transcription regulator(s) as the budding yeast Mcm1-Ndd1-Fkh2 transcription factor complex in controlling occurrence of cytokinesis.