Enhanced Sensitivity of Nonsmall Cell Lung Cancer with Acquired Resistance to Epidermal Growth Factor Receptor-Tyrosine Kinase Inhibitors to Phenformin: The Roles of a Metabolic Shift to Oxidative Phosphorylation and Redox Balance.

BACKGROUND: Although the efficacy of epidermal growth factor receptor-tyrosine kinase inhibitor (EGFR- TKI) therapy has been proven in non-small cell lung cancer (NSCLC) patients, acquired resistance to EGFR-TKIs presents a serious clinical problem. Hence, the identification of new therapeutic strategy is needed to treat EGFR-TKI-resistant NSCLC. METHODS: Acquired EGFR-TKI-resistant lung cancer cell lines (HCC827, H1993, and H292 cells with acquired resistance to gefitinib or erlotinib) were used for cell-based studies. IncuCyte live cell analysis system and XFp analyzer were used for the determination of cell proliferation and energy metabolism, respectively. In vivo anticancer effect of phenformin was assessed in xenografts implanting HCC827 and gefitinib-resistant HCC827 (HCC827 GR) cells. RESULTS: HCC827 GR and erlotinib-resistant H1993 (H1993 ER) cells exhibited different metabolic properties compared with their respective parental cells, HCC827, and H1993. In EGFR-TKI-resistant NSCLC cells, glycolysis markers including the glucose consumption rate, intracellular lactate level, and extracellular acidification rate were decreased; however, mitochondrial oxidative phosphorylation (OXPHOS) markers including mitochondria-driven ATP production, mitochondrial membrane potential, and maximal OXPHOS capacity were increased. Cell proliferation and tumor growth were strongly inhibited by biguanide phenformin via targeting of mitochondrial OXPHOS complex 1 in EGFR-TKI-resistant NSCLC cells. Inhibition of OXPHOS resulted in a reduced NAD+/NADH ratio and intracellular aspartate levels. Recovery of glycolysis by hexokinase 2 overexpression in erlotinib-resistant H292 (H292 ER) cells significantly reduced the anticancer effects of phenformin. CONCLUSION: Long-term treatment with EGFR-TKIs causes reactivation of mitochondrial metabolism, resulting in vulnerability to OXPHOS inhibitor such as phenformin. We propose a new therapeutic option for NSCLC with acquired EGFR-TKI resistance that focuses on cancer metabolism.

Erratum: Modulation of mRNA and lncRNA expression dynamics by the Set2-Rpd3S pathway.

This corrects the article DOI: 10.1038/ncomms13534.

Modulation of mRNA and lncRNA expression dynamics by the Set2-Rpd3S pathway.

H3K36 methylation by Set2 targets Rpd3S histone deacetylase to transcribed regions of mRNA genes, repressing internal cryptic promoters and slowing elongation. Here we explore the function of this pathway by analysing transcription in yeast undergoing a series of carbon source shifts. Approximately 80 mRNA genes show increased induction upon SET2 deletion. A majority of these promoters have overlapping lncRNA transcription that targets H3K36me3 and deacetylation by Rpd3S to the mRNA promoter. We previously reported a similar mechanism for H3K4me2-mediated repression via recruitment of the Set3C histone deacetylase. Here we show that the distance between an mRNA and overlapping lncRNA promoter determines whether Set2-Rpd3S or Set3C represses. This analysis also reveals many previously unreported cryptic ncRNAs induced by specific carbon sources, showing that cryptic promoters can be environmentally regulated. Therefore, in addition to repression of cryptic transcription and modulation of elongation, H3K36 methylation maintains optimal expression dynamics of many mRNAs and ncRNAs.

Transcriptomic analysis of insulin-sensitive tissues from anti-diabetic drug treated ZDF rats, a T2DM animal model.

Gene expression changes have been associated with type 2 diabetes mellitus (T2DM); however, the alterations are not fully understood. We investigated the effects of anti-diabetic drugs on gene expression in Zucker diabetic fatty (ZDF) rats using oligonucleotide microarray technology to identify gene expression changes occurring in T2DM. Global gene expression in the pancreas, adipose tissue, skeletal muscle, and liver was profiled from Zucker lean control (ZLC) and anti-diabetic drug treated ZDF rats compared with those in ZDF rats. We showed that anti-diabetic drugs regulate the expression of a large number of genes. We provided a more integrated view of the diabetic changes by examining the gene expression networks. The resulting sub-networks allowed us to identify several biological processes that were significantly enriched by the anti-diabetic drug treatment, including oxidative phosphorylation (OXPHOS), systemic lupus erythematous, and the chemokine signaling pathway. Among them, we found that white adipose tissue from ZDF rats showed decreased expression of a set of OXPHOS genes that were normalized by rosiglitazone treatment accompanied by rescued blood glucose levels. In conclusion, we suggest that alterations in OXPHOS gene expression in white adipose tissue may play a role in the pathogenesis and drug mediated recovery of T2DM through a comprehensive gene expression network study after multi-drug treatment of ZDF rats.

Ogura-CMS in Chinese cabbage (Brassica rapa ssp. pekinensis) causes delayed expression of many nuclear genes.

We investigated the mechanism regulating cytoplasmic male sterility (CMS) in Brassica rapa ssp. pekinensis using floral bud transcriptome analyses of Ogura-CMS Chinese cabbage and its maintainer line in B. rapa 300-K oligomeric probe (Br300K) microarrays. Ogura-CMS Chinese cabbage produced few and infertile pollen grains on indehiscent anthers. Compared to the maintainer line, CMS plants had shorter filaments and plant growth, and delayed flowering and pollen development. In microarray analysis, 4646 genes showed different expression, depending on floral bud size, between Ogura-CMS and its maintainer line. We found 108 and 62 genes specifically expressed in Ogura-CMS and its maintainer line, respectively. Ogura-CMS line-specific genes included stress-related, redox-related, and B. rapa novel genes. In the maintainer line, genes related to pollen coat and germination were specifically expressed in floral buds longer than 3mm, suggesting insufficient expression of these genes in Ogura-CMS is directly related to dysfunctional pollen. In addition, many nuclear genes associated with auxin response, ATP synthesis, pollen development and stress response had delayed expression in Ogura-CMS plants compared to the maintainer line, which is consistent with the delay in growth and development of Ogura-CMS plants. Delayed expression may reduce pollen grain production and/or cause sterility, implying that mitochondrial, retrograde signaling delays nuclear gene expression.

hiPathDB: a human-integrated pathway database with facile visualization.

One of the biggest challenges in the study of biological regulatory networks is the systematic organization and integration of complex interactions taking place within various biological pathways. Currently, the information of the biological pathways is dispersed in multiple databases in various formats. hiPathDB is an integrated pathway database that combines the curated human pathway data of NCI-Nature PID, Reactome, BioCarta and KEGG. In total, it includes 1661 pathways consisting of 8976 distinct physical entities. hiPathDB provides two different types of integration. The pathway-level integration, conceptually a simple collection of individual pathways, was achieved by devising an elaborate model that takes distinct features of four databases into account and subsequently reformatting all pathways in accordance with our model. The entity-level integration creates a single unified pathway that encompasses all pathways by merging common components. Even though the detailed molecular-level information such as complex formation or post-translational modifications tends to be lost, such integration makes it possible to investigate signaling network over the entire pathways and allows identification of pathway cross-talks. Another strong merit of hiPathDB is the built-in pathway visualization module that supports explorative studies of complex networks in an interactive fashion. The layout algorithm is optimized for virtually automatic visualization of the pathways. hiPathDB is available at http://hiPathDB.kobic.re.kr.

GARNET--gene set analysis with exploration of annotation relations.

BACKGROUND: Gene set analysis is a powerful method of deducing biological meaning for an a priori defined set of genes. Numerous tools have been developed to test statistical enrichment or depletion in specific pathways or gene ontology (GO) terms. Major difficulties towards biological interpretation are integrating diverse types of annotation categories and exploring the relationships between annotation terms of similar information. RESULTS: GARNET (Gene Annotation Relationship NEtwork Tools) is an integrative platform for gene set analysis with many novel features. It includes tools for retrieval of genes from annotation database, statistical analysis & visualization of annotation relationships, and managing gene sets. In an effort to allow access to a full spectrum of amassed biological knowledge, we have integrated a variety of annotation data that include the GO, domain, disease, drug, chromosomal location, and custom-defined annotations. Diverse types of molecular networks (pathways, transcription and microRNA regulations, protein-protein interaction) are also included. The pair-wise relationship between annotation gene sets was calculated using kappa statistics. GARNET consists of three modules--gene set manager, gene set analysis and gene set retrieval, which are tightly integrated to provide virtually automatic analysis for gene sets. A dedicated viewer for annotation network has been developed to facilitate exploration of the related annotations. CONCLUSIONS: GARNET (gene annotation relationship network tools) is an integrative platform for diverse types of gene set analysis, where complex relationships among gene annotations can be easily explored with an intuitive network visualization tool (http://garnet.isysbio.org/ or http://ercsb.ewha.ac.kr/garnet/).

Clinical validation of colorectal cancer biomarkers identified from bioinformatics analysis of public expression data.

PURPOSE: Identification of novel biomarkers of cancer is important for improved diagnosis, prognosis, and therapeutic intervention. This study aimed to identify marker genes of colorectal cancer (CRC) by combining bioinformatics analysis of gene expression data and validation experiments using patient samples and to examine the potential connection between validated markers and the established oncogenes such as c-Myc and K-ras. EXPERIMENTAL DESIGN: Publicly available data from GenBank and Oncomine were meta-analyzed leading to 34 candidate marker genes of CRC. Multiple case-matched normal and tumor tissues were examined by RT-PCR for differential expression, and 9 genes were validated as CRC biomarkers. Statistical analyses for correlation with major clinical parameters were carried out, and RNA interference was used to examine connection with major oncogenes. RESULTS: We show with high confidence that 9 (ECT2, ETV4, DDX21, RAN, S100A11, RPS4X, HSPD1, CKS2, and C9orf140) of the 34 candidate genes are expressed at significantly elevated levels in CRC tissues compared to normal tissues. Furthermore, high-level expression of RPS4X was associated with nonmucinous cancer cell type and that of ECT2 with lack of lymphatic invasion while upregulation of CKS2 was correlated with early tumor stage and lack of family history of CRC. We also demonstrate that RPS4X and DDX21 are regulatory targets of c-Myc and ETV4 is downstream to K-ras signaling. CONCLUSIONS: We have identified multiple novel biomarkers of CRC. Further analyses of their function and connection to signaling pathways may reveal potential value of these biomarkers in diagnosis, prognosis, and treatment of CRC.

Age-dependent evolution of the yeast protein interaction network suggests a limited role of gene duplication and divergence.

Proteins interact in complex protein-protein interaction (PPI) networks whose topological properties-such as scale-free topology, hierarchical modularity, and dissortativity-have suggested models of network evolution. Currently preferred models invoke preferential attachment or gene duplication and divergence to produce networks whose topology matches that observed for real PPIs, thus supporting these as likely models for network evolution. Here, we show that the interaction density and homodimeric frequency are highly protein age-dependent in real PPI networks in a manner which does not agree with these canonical models. In light of these results, we propose an alternative stochastic model, which adds each protein sequentially to a growing network in a manner analogous to protein crystal growth (CG) in solution. The key ideas are (1) interaction probability increases with availability of unoccupied interaction surface, thus following an anti-preferential attachment rule, (2) as a network grows, highly connected sub-networks emerge into protein modules or complexes, and (3) once a new protein is committed to a module, further connections tend to be localized within that module. The CG model produces PPI networks consistent in both topology and age distributions with real PPI networks and is well supported by the spatial arrangement of protein complexes of known 3-D structure, suggesting a plausible physical mechanism for network evolution.

A critical assessment of Mus musculus gene function prediction using integrated genomic evidence.

BACKGROUND: Several years after sequencing the human genome and the mouse genome, much remains to be discovered about the functions of most human and mouse genes. Computational prediction of gene function promises to help focus limited experimental resources on the most likely hypotheses. Several algorithms using diverse genomic data have been applied to this task in model organisms; however, the performance of such approaches in mammals has not yet been evaluated. RESULTS: In this study, a standardized collection of mouse functional genomic data was assembled; nine bioinformatics teams used this data set to independently train classifiers and generate predictions of function, as defined by Gene Ontology (GO) terms, for 21,603 mouse genes; and the best performing submissions were combined in a single set of predictions. We identified strengths and weaknesses of current functional genomic data sets and compared the performance of function prediction algorithms. This analysis inferred functions for 76% of mouse genes, including 5,000 currently uncharacterized genes. At a recall rate of 20%, a unified set of predictions averaged 41% precision, with 26% of GO terms achieving a precision better than 90%. CONCLUSION: We performed a systematic evaluation of diverse, independently developed computational approaches for predicting gene function from heterogeneous data sources in mammals. The results show that currently available data for mammals allows predictions with both breadth and accuracy. Importantly, many highly novel predictions emerge for the 38% of mouse genes that remain uncharacterized.

Inferring mouse gene functions from genomic-scale data using a combined functional network/classification strategy.

The complete set of mouse genes, as with the set of human genes, is still largely uncharacterized, with many pieces of experimental evidence accumulating regarding the activities and expression of the genes, but the majority of genes as yet still of unknown function. Within the context of the MouseFunc competition, we developed and applied two distinct large-scale data mining approaches to infer the functions (Gene Ontology annotations) of mouse genes from experimental observations from available functional genomics, proteomics, comparative genomics, and phenotypic data. The two strategies - the first using classifiers to map features to annotations, the second propagating annotations from characterized genes to uncharacterized genes along edges in a network constructed from the features - offer alternative and possibly complementary approaches to providing functional annotations. Here, we re-implement and evaluate these approaches and their combination for their ability to predict the proper functional annotations of genes in the MouseFunc data set. We show that, when controlling for the same set of input features, the network approach generally outperformed a naive Bayesian classifier approach, while their combination offers some improvement over either independently. We make our observations of predictive performance on the MouseFunc competition hold-out set, as well as on a ten-fold cross-validation of the MouseFunc data. Across all 1,339 annotated genes in the MouseFunc test set, the median predictive power was quite strong (median area under a receiver operating characteristic plot of 0.865 and average precision of 0.195), indicating that a mining-based strategy with existing data is a promising path towards discovering mammalian gene functions. As one product of this work, a high-confidence subset of the functional mouse gene network was produced - spanning >70% of mouse genes with >1.6 million associations - that is predictive of mouse (and therefore often human) gene function and functional associations. The network should be generally useful for mammalian gene functional analyses, such as for predicting interactions, inferring functional connections between genes and pathways, and prioritizing candidate genes. The network and all predictions are available on the worldwide web.

Using structural motif descriptors for sequence-based binding site prediction.

BACKGROUND: Many protein sequences are still poorly annotated. Functional characterization of a protein is often improved by the identification of its interaction partners. Here, we aim to predict protein-protein interactions (PPI) and protein-ligand interactions (PLI) on sequence level using 3D information. To this end, we use machine learning to compile sequential segments that constitute structural features of an interaction site into one profile Hidden Markov Model descriptor. The resulting collection of descriptors can be used to screen sequence databases in order to predict functional sites. RESULTS: We generate descriptors for 740 classified types of protein-protein binding sites and for more than 3,000 protein-ligand binding sites. Cross validation reveals that two thirds of the PPI descriptors are sufficiently conserved and significant enough to be used for binding site recognition. We further validate 230 PPIs that were extracted from the literature, where we additionally identify the interface residues. Finally we test ligand-binding descriptors for the case of ATP. From sequences with Swiss-Prot annotation "ATP-binding", we achieve a recall of 25% with a precision of 89%, whereas Prosite's P-loop motif recognizes an equal amount of hits at the expense of a much higher number of false positives (precision: 57%). Our method yields 771 hits with a precision of 96% that were not previously picked up by any Prosite-pattern. CONCLUSION: The automatically generated descriptors are a useful complement to known Prosite/InterPro motifs. They serve to predict protein-protein as well as protein-ligand interactions along with their binding site residues for proteins where merely sequence information is available.

The many faces of protein-protein interactions: A compendium of interface geometry.

A systematic classification of protein-protein interfaces is a valuable resource for understanding the principles of molecular recognition and for modelling protein complexes. Here, we present a classification of domain interfaces according to their geometry. Our new algorithm uses a hybrid approach of both sequential and structural features. The accuracy is evaluated on a hand-curated dataset of 416 interfaces. Our hybrid procedure achieves 83% precision and 95% recall, which improves the earlier sequence-based method by 5% on both terms. We classify virtually all domain interfaces of known structure, which results in nearly 6,000 distinct types of interfaces. In 40% of the cases, the interacting domain families associate in multiple orientations, suggesting that all the possible binding orientations need to be explored for modelling multidomain proteins and protein complexes. In general, hub proteins are shown to use distinct surface regions (multiple faces) for interactions with different partners. Our classification provides a convenient framework to query genuine gene fusion, which conserves binding orientation in both fused and separate forms. The result suggests that the binding orientations are not conserved in at least one-third of the gene fusion cases detected by a conventional sequence similarity search. We show that any evolutionary analysis on interfaces can be skewed by multiple binding orientations and multiple interaction partners. The taxonomic distribution of interface types suggests that ancient interfaces common to the three major kingdoms of life are enriched by symmetric homodimers. The classification results are online at http://www.scoppi.org.

Hispidin derivatives from the mushroom Inonotus xeranticus and their antioxidant activity.

In an effort to identify antioxidants from edible and medicinal mushrooms, three new hispidin derivatives, methylinoscavin A (2), inoscavin B (4), and methylinoscavin B (5), together with the known compounds inoscavin A and phelligridin F, were isolated from the methanolic extract of the fruiting bodies of Inonotus xeranticus. Their structures were determined on the basis of spectroscopic analyses.

SCOPPI: a structural classification of protein-protein interfaces.

SCOPPI, the structural classification of protein-protein interfaces, is a comprehensive database that classifies and annotates domain interactions derived from all known protein structures. SCOPPI applies SCOP domain definitions and a distance criterion to determine inter-domain interfaces. Using a novel method based on multiple sequence and structural alignments of SCOP families, SCOPPI presents a comprehensive geometrical classification of domain interfaces. Various interface characteristics such as number, type and position of interacting amino acids, conservation, interface size, and permanent or transient nature of the interaction are further provided. Proteins in SCOPPI are annotated with Gene Ontology terms, and the ontology can be used to quickly browse SCOPPI. Screenshots are available for every interface and its participating domains. Here, we describe contents and features of the web-based user interface as well as the underlying methods used to generate SCOPPI's data. In addition, we present a number of examples where SCOPPI becomes a useful tool to analyze viral mimicry of human interface binding sites, gene fusion events, conservation of interface residues and diversity of interface localizations. SCOPPI is available at http://www.scoppi.org.

Equivalent binding sites reveal convergently evolved interaction motifs.

MOTIVATION: Much research has been devoted to the characterization of interaction interfaces found in complexes with known structure. In this context, the interactions of non-homologous domains at equivalent binding sites are of particular interest, as they can reveal convergently evolved interface motifs. Such motifs are an important source of information to formulate rules for interaction specificity and to design ligands based on the common features shared among diverse partners. RESULTS: We develop a novel method to identify non-homologous structural domains which bind at equivalent sites when interacting with a common partner. We systematically apply this method to all pairs of interactions with known structure and derive a comprehensive database for these interactions. Of all non-homologous domains, which bind with a common interaction partner, 4.2% use the same interface of the common interaction partner (excluding immunoglobulins and proteases). This rises to 16% if immunoglobulin and proteases are included. We demonstrate two applications of our database: first, the systematic screening for viral protein interfaces, which can mimic native interfaces and thus interfere; and second, structural motifs in enzymes and its inhibitors. We highlight several cases of virus protein mimicry: viral M3 protein interferes with a chemokine dimer interface. The virus has evolved the motif SVSPLP, which mimics the native SSDTTP motif. A second example is the regulatory factor Nef in HIV which can mimic a kinase when interacting with SH3. Among others the virus has evolved the kinase's PxxP motif. Further, we elucidate motif resemblances in Baculovirus p35 and HIV capsid proteins. Finally, chymotrypsin is subject to scrutiny wrt. its structural similarity to subtilisin and wrt. its inhibitor's similar recognition sites. SUPPLEMENTARY INFORMATION: A database is online at scoppi.biotec.tu-dresden.de/abac/.

Survey of the geometric association of domain-domain interfaces.

Considering the limited success of the most sophisticated docking methods available and the amount of computation required for systematic docking, cataloging all the known interfaces may be an alternative basis for the prediction of protein tertiary and quaternary structures. We classify domain interfaces according to the geometry of domain-domain association. By applying a simple and efficient method called "interface tag clustering," more than 4,000 distinct types of domain interfaces are collected from Protein Quaternary Structure Server and Protein Data Bank. Given a pair of interacting domains, we define "face" as the set of interacting residues in each single domain and the pair of interacting faces as an "interface." We investigate how the geometry of interfaces relates to a network of interacting protein families, such as how many different binding orientations are possible between two families or whether a family uses distinct surfaces or the same surface when the family has diverse interaction partners from various families. We show there are, on average, 1.2-1.9 different types of interfaces between interacting domains and a significant number of family pairs associate in multiple orientations. In general, a family tends to use distinct faces for each partner when the family has diverse interaction partners. Each face is highly specific to its interaction partner and the binding orientation. The relative positions of interface residues are generally well conserved within the same type of interface even between remote homologs. The classification result is available at http://www.biotec.tu-dresden.de/~wkim/supplement.

Large scale statistical prediction of protein-protein interaction by potentially interacting domain (PID) pair.

Protein-protein interaction plays a critical role in biological processes. The identification of interacting proteins by computational methods can provide new leads in functional studies of uncharacterized proteins without performing extensive experiments. We developed a database for the potentially interacting domain pairs (PID) extracted from a dataset of experimentally identified interacting protein pairs (DIP: database of interacting proteins) with InterPro, an integrated database of protein families, domains and functional sites. In developing protein interaction databases and predictive methods, sensitive statistical scoring systems is critical to provide a reliability index for accurate functional analysis of interaction networks. We present a statistical scoring system, named "PID matrix score" as a measure of the interaction probability (interactability) between domains. This system provided a valuable tool for functional prediction of unknown proteins. For the evaluation of PID matrix, cross validation was performed with subsets of DIP data. The prediction system gives about 50% sensitivity and more than 98% specificity, which implies that the information for interacting proteins pairs could be enriched about 30 fold with the PID matrix. It is demonstrated that mapping of the genome-wide interaction network can be achieved by using the PID matrix.