Integrative analysis of ocular complications in atherosclerosis unveils pathway convergence and crosstalk

Abstract

Atherosclerosis is a life-threatening disease and a major cause of mortalities worldwide. While many of the atherosclerotic sequelae are reflected as microvascular effects in the eye, the molecular mechanisms of their development is not yet known. In this study, we employed a systems biology approach to unveil the most significant events and key molecular mediators of ophthalmic sequelae caused by atherosclerosis. Literature mining was used to identify the proteins involved in both atherosclerosis and ophthalmic diseases. A protein–protein interaction (PPI) network was prepared using the literature-mined seed nodes. Network topological analysis was carried out using Cytoscape, while network nodes were annotated using database for annotation, visualization and integrated discovery in order to identify the most enriched pathways and processes. Network analysis revealed that mitogen-activated protein kinase 1 (MAPK1) and protein kinase C occur with highest betweenness centrality, degree and closeness centrality, thus reflecting their functional importance to the network. Our analysis shows that atherosclerosis-associated ophthalmic complications are caused by the convergence of neurotrophin signaling pathways, multiple immune response pathways and focal adhesion pathway on the MAPK signaling pathway. The PPI network shares features with vasoregression, a process underlying multiple vascular eye diseases. Our study presents a first clear and composite picture of the components and crosstalk of the main pathways of atherosclerosis-induced ocular diseases. The hub bottleneck nodes highlight the presence of molecules important for mediating the ophthalmic complications of atherosclerosis and contain five established drug targets for future therapeutic modulation efforts.     

Atopic dermatitis-associated protein interaction network lead to new insights in chronic sulfur mustard skin lesion mechanisms

Chronic sulfur mustard skin lesions (CSMSLs) are the most common complications of sulfur mustard exposure; however, its mechanism is not completely understood.According to clinical signs, there are similarities between CSMSL and atopic dermatitis (AD). In this study, proteomic results of AD were reviewed and the AD-associated protein–protein interaction network (PIN) was analyzed. According to centrality measurements, 16 proteins were designated as pivotal elements in AD mechanisms. Interestingly, most of these proteins had been reported in some sulfur mustard-related studies in late and acute phases separately. Based on the gene enrichment analysis, aging, cell response to stress, cancer, Toll- and NOD-like receptor and apoptosis signaling pathways have the greatest impact on the disease. By the analysis of directed protein interaction networks, it is concluded that TNF, IL-6, AKT1, NOS3 and CDKN1A are the most important proteins. It is possible that these proteins play role in the shared complications of AD and CSMSL including xerosis and itching.     

Correlation between Oncogenic Mutations and Parameter Sensitivity of the Apoptosis Pathway Model

One of the major breakthroughs in oncogenesis research in recent years is the discovery that, in most patients, oncogenic mutations are concentrated in a few core biological functional pathways. This discovery indicates that oncogenic mechanisms are highly related to the dynamics of biologic regulatory networks, which govern the behaviour of functional pathways. Here, we propose that oncogenic mutations found in different biological functional pathways are closely related to parameter sensitivity of the corresponding networks. To test this hypothesis, we focus on the DNA damage-induced apoptotic pathway—the most important safeguard against oncogenesis. We first built the regulatory network that governs the apoptosis pathway, and then translated the network into dynamics equations. Using sensitivity analysis of the network parameters and comparing the results with cancer gene mutation spectra, we found that parameters that significantly affect the bifurcation point correspond to high-frequency oncogenic mutations. This result shows that the position of the bifurcation point is a better measure of the functionality of a biological network than gene expression levels of certain key proteins. It further demonstrates the suitability of applying systems-level analysis to biological networks as opposed to studying genes or proteins in isolation.     

Dominating biological networks

Proteins are essential macromolecules of life that carry out most cellular processes. Since proteins aggregate to perform function, and since protein-protein interaction (PPI) networks model these aggregations, one would expect to uncover new biology from PPI network topology. Hence, using PPI networks to predict protein function and role of protein pathways in disease has received attention. A debate remains open about whether network properties of "biologically central (BC)" genes (i.e., their protein products), such as those involved in aging, cancer, infectious diseases, or signaling and drug-targeted pathways, exhibit some topological centrality compared to the rest of the proteins in the human PPI network.To help resolve this debate, we design new network-based approaches and apply them to get new insight into biological function and disease. We hypothesize that BC genes have a topologically central (TC) role in the human PPI network. We propose two different concepts of topological centrality. We design a new centrality measure to capture complex wirings of proteins in the network that identifies as TC those proteins that reside in dense extended network neighborhoods. Also, we use the notion of domination and find dominating sets (DSs) in the PPI network, i.e., sets of proteins such that every protein is either in the DS or is a neighbor of the DS. Clearly, a DS has a TC role, as it enables efficient communication between different network parts. We find statistically significant enrichment in BC genes of TC nodes and outperform the existing methods indicating that genes involved in key biological processes occupy topologically complex and dense regions of the network and correspond to its "spine" that connects all other network parts and can thus pass cellular signals efficiently throughout the network. To our knowledge, this is the first study that explores domination in the context of PPI networks.     

Computational Modeling of Allosteric Regulation in the Hsp90 Chaperones: A Statistical Ensemble Analysis of Protein Structure Networks and Allosteric Communications

A fundamental role of the Hsp90 chaperone in regulating functional activity of diverse protein clients is essential for the integrity of signaling networks. In this work we have combined biophysical simulations of the Hsp90 crystal structures with the protein structure network analysis to characterize the statistical ensemble of allosteric interaction networks and communication pathways in the Hsp90 chaperones. We have found that principal structurally stable communities could be preserved during dynamic changes in the conformational ensemble. The dominant contribution of the inter-domain rigidity to the interaction networks has emerged as a common factor responsible for the thermodynamic stability of the active chaperone form during the ATPase cycle. Structural stability analysis using force constant profiling of the inter-residue fluctuation distances has identified a network of conserved structurally rigid residues that could serve as global mediating sites of allosteric communication. Mapping of the conformational landscape with the network centrality parameters has demonstrated that stable communities and mediating residues may act concertedly with the shifts in the conformational equilibrium and could describe the majority of functionally significant chaperone residues. The network analysis has revealed a relationship between structural stability, global centrality and functional significance of hotspot residues involved in chaperone regulation. We have found that allosteric interactions in the Hsp90 chaperone may be mediated by modules of structurally stable residues that display high betweenness in the global interaction network. The results of this study have suggested that allosteric interactions in the Hsp90 chaperone may operate via a mechanism that combines rapid and efficient communication by a single optimal pathway of structurally rigid residues and more robust signal transmission using an ensemble of suboptimal multiple communication routes. This may be a universal requirement encoded in protein structures to balance the inherent tension between resilience and efficiency of the residue interaction networks.     

A systematic survey of centrality measures for protein-protein interaction networks

Background

Numerous centrality measures have been introduced to identify “central” nodes in large networks. The availability of a wide range of measures for ranking influential nodes leaves the user to decide which measure may best suit the analysis of a given network. The choice of a suitable measure is furthermore complicated by the impact of the network topology on ranking influential nodes by centrality measures. To approach this problem systematically, we examined the centrality profile of nodes of yeast protein-protein interaction networks (PPINs) in order to detect which centrality measure is succeeding in predicting influential proteins. We studied how different topological network features are reflected in a large set of commonly used centrality measures.

Results

We used yeast PPINs to compare 27 common of centrality measures. The measures characterize and assort influential nodes of the networks. We applied principal component analysis (PCA) and hierarchical clustering and found that the most informative measures depend on the network’s topology. Interestingly, some measures had a high level of contribution in comparison to others in all PPINs, namely Latora closeness, Decay, Lin, Freeman closeness, Diffusion, Residual closeness and Average distance centralities.

Conclusions

The choice of a suitable set of centrality measures is crucial for inferring important functional properties of a network. We concluded that undertaking data reduction using unsupervised machine learning methods helps to choose appropriate variables (centrality measures). Hence, we proposed identifying the contribution proportions of the centrality measures with PCA as a prerequisite step of network analysis before inferring functional consequences, e.g., essentiality of a node.

     

Bioinformatics role of the WGCNA analysis and co-expression network identifies of prognostic marker in lung cancer

Lung cancer is the most talked about cancer in the world. It is also one of the cancers that currently has a high mortality rate. The aim of our research is to find more effective therapeutic targets and prognostic markers for human lung cancer. First, we download gene expression data from the GEO database. We performed weighted co-expression network analysis on the selected genes, we then constructed scale-free networks and topological overlap matrices, and performed correlation modular analysis with the cancer group. We screened the 200 genes with the highest correlation in the cyan module for functional enrichment analysis and protein interaction network construction, found that most of them focused on cell division, tumor necrosis factor-mediated signaling pathways, cellular redox homeostasis, reactive oxygen species biosynthesis, and other processes, and were related to the cell cycle, apoptosis, HIF-1 signaling pathway, p53 signaling pathway, NF-κB signaling pathway, and several cancer disease pathways are involved. Finally, we used the GEPIA website data to perform survival analysis on some of the genes with GS > 0.6 in the cyan module. CBX3, AHCY, MRPL12, TPGB, TUBG1, KIF11, LRRC59, MRPL17, TMEM106B, ZWINT, TRIP13, and HMMR was identified as an important prognostic factor for lung cancer patients. In summary, we identified 12 mRNAs associated with lung cancer prognosis. Our study contributes to a deeper understanding of the molecular mechanisms of lung cancer and provides new insights into drug use and prognosis.     

Localized network centrality and essentiality in the yeast–protein interaction network

It has been suggested that a close relationship exists between gene essentiality and network centrality in protein–protein interaction networks. However, recent studies have reported somewhat conflicting results on this relationship. In this study, we investigated whether essential proteins could be inferred from network centrality alone. In addition, we determined which centrality measures describe the essentiality well. For this analysis, we devised new local centrality measures based on several well‐known centrality measures to more precisely describe the connection between network topology and essentiality. We examined two recent yeast protein–protein interaction networks using 40 different centrality measures. We discovered a close relationship between the path‐based localized information centrality and gene essentiality, which suggested underlying topological features that represent essentiality. We propose that two important features of the localized information centrality (proper representation of environmental complexity and the consideration of local sub‐networks) are the key factors that reveal essentiality. In addition, a random forest classifier showed reasonable performance at classifying essential proteins. Finally, the results of clustering analysis using centrality measures indicate that some network clusters are closely related with both particular biological processes and essentiality, suggesting that functionally related proteins tend to share similar network properties.     

Preeclampsia: a bioinformatics approach through protein-protein interaction networks analysis

ABSTRACT: BACKGROUND: In this study we explored preeclampsia through a bioinformatics approach. We create a comprehensive genes/proteins dataset by the analysis of both public proteomic data and text mining of public scientific literature. From this dataset the associated protein-protein interaction network has been obtained. Several indexes of centrality have been explored for hubs detection as well as the enrichment statistical analysis of metabolic pathway and disease. RESULTS: We confirmed the well known relationship between preeclampsia and cardiovascular diseases but also identified statistically significant relationships with respect to cancer and aging. Moreover, significant metabolic pathways such as apoptosis, cancer and cytokine-cytokine receptor interaction have also been identified by enrichment analysis. We obtained FLT1, VEGFA, FN1, F2 and PGF genes with the highest scores by hubs analysis; however, we also found other genes as PDIA3, LYN, SH2B2 and NDRG1 with high scores. CONCLUSIONS: The applied methodology not only led to the identification of well known genes related to preeclampsia but also to propose new candidates poorly explored or completely unknown in the pathogenesis of preeclampsia, which eventually need to be validated experimentally. Moreover, new possible connections were detected between preeclampsia and other diseases that could open new areas of research. More must be done in this area to resolve the identification of unknown interactions of proteins/genes and also for a better integration of metabolic pathways and diseases.     

Probabilistic model of the human protein-protein interaction network

A catalog of all human protein-protein interactions would provide scientists with a framework to study protein deregulation in complex diseases such as cancer. Here we demonstrate that a probabilistic analysis integrating model organism interactome data, protein domain data, genome-wide gene expression data and functional annotation data predicts nearly 40,000 protein-protein interactions in humans-a result comparable to those obtained with experimental and computational approaches in model organisms. We validated the accuracy of the predictive model on an independent test set of known interactions and also experimentally confirmed two predicted interactions relevant to human cancer, implicating uncharacterized proteins into definitive pathways. We also applied the human interactome network to cancer genomics data and identified several interaction subnetworks activated in cancer. This integrative analysis provides a comprehensive framework for exploring the human protein interaction network.     

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.     

Network analysis and the impact of Aflibercept on specific mediators of angiogenesis in HUVEC cells

Angiogenesis, inflammation and endothelial cells’ migration and proliferation exert fundamental roles in different diseases. However, more studies are needed to identify key proteins and pathways involved in these processes. Aflibercept has received the approval of the US Food and Drug Administration (FDA) for the treatment of wet AMD and colorectal cancer. Moreover, the effect of Aflibercept on VEGFR2 downstream signalling pathways has not been investigated yet. Here, we integrated text mining data, protein-protein interaction networks and multi-experiment microarray data to specify candidate genes that are involved in VEGFA/VEGFR2 signalling pathways. Network analysis of candidate genes determined the importance of the nominated genes via different centrality parameters. Thereupon, several genes—with the highest centrality indexes—were recruited to investigate the impact of Aflibercept on their expression pattern in HUVEC cells. Real-time PCR was performed, and relative expression of the specific genes revealed that Aflibercept modulated angiogenic process by VEGF/PI3KA/AKT/mTOR axis, invasion by MMP14/MMP9 axis and inflammation-related angiogenesis by IL-6-STAT3 axis. Data showed Aflibercept simultaneously affected these processes and determined the nominated axes that had been affected by the drug. Furthermore, integrating the results of Aflibercept on expression of candidate genes with the current network analysis suggested that resistance against the Aflibercept effect is a plausible process in HUVEC cells.     

Human papillomavirus elevated genetic biomarker signature by statistical algorithm

Head and neck squamous cell carcinoma (HNSCC) is the one of the most frequently found cancers in the world. The aim of the study was to find the genes responsible and enriched pathways associated with HNSCC using bioinformatics and survival analysis methods. A total of 646 patients with HNSCC based on clinical information were considered for the study. HNSCC samples were grouped according to the parameters (RFS, DFS, PFS, or OS). The probe ID of these 11 genes was retrieved by Affymetrix using the NetAffx Query algorithm. The protein–protein interaction (PPI) network and Kaplan–Meier curve were used to find associations among the genes' expression data. We found that among these 11 genes, nine genes, CCNA1, MMP3, FLRT3, GJB6, ZFR2, PITX2, SYCP2, MEI1, and UGT8 were significant (p < .05). A survival plot was drawn between the p value and gene expression. This study helped us find the nine significant genes which play vital roles in HNSCC along with their key pathways and their interaction with other genes in the PPI network. Finally, we found the biomarker index for relapse time and risk factors for HNSCC in cancer patients.     

The organisational structure of protein networks: revisiting the centrality–lethality hypothesis

Protein networks, describing physical interactions as well as functional associations between proteins, have been unravelled for many organisms in the recent past. Databases such as the STRING provide excellent resources for the analysis of such networks. In this contribution, we revisit the organisation of protein networks, particularly the centrality–lethality hypothesis, which hypothesises that nodes with higher centrality in a network are more likely to produce lethal phenotypes on removal, compared to nodes with lower centrality. We consider the protein networks of a diverse set of 20 organisms, with essentiality information available in the Database of Essential Genes and assess the relationship between centrality measures and lethality. For each of these organisms, we obtained networks of high-confidence interactions from the STRING database, and computed network parameters such as degree, betweenness centrality, closeness centrality and pairwise disconnectivity indices. We observe that the networks considered here are predominantly disassortative. Further, we observe that essential nodes in a network have a significantly higher average degree and betweenness centrality, compared to the network average. Most previous studies have evaluated the centrality–lethality hypothesis for Saccharomyces cerevisiae and Escherichia coli; we here observe that the centrality–lethality hypothesis hold goods for a large number of organisms, with certain limitations. Betweenness centrality may also be a useful measure to identify essential nodes, but measures like closeness centrality and pairwise disconnectivity are not significantly higher for essential nodes.     

Molecular Architecture of Spinal Cord Injury Protein Interaction Network

Spinal cord injury (SCI) is associated with complex pathophysiological processes that follow the primary traumatic event and determine the extent of secondary damage and functional recovery. Numerous reports have used global and hypothesis-driven approaches to identify protein changes that contribute to the overall pathology of SCI in an effort to identify potential therapeutic interventions. In this study, we use a semi-automatic annotation approach to detect terms referring to genes or proteins dysregulated in the SCI literature and develop a curated SCI interactome. Network analysis of the SCI interactome revealed the presence of a rich-club organization corresponding to a “powerhouse” of highly interacting hub-proteins. Studying the modular organization of the network have shown that rich-club proteins cluster into modules that are specifically enriched for biological processes that fall under the categories of cell death, inflammation, injury recognition and systems development. Pathway analysis of the interactome and the rich-club revealed high similarity indicating the role of the rich-club proteins as hubs of the most prominent pathways in disease pathophysiology and illustrating the centrality of pro-and anti-survival signal competition in the pathology of SCI. In addition, evaluation of centrality measures of single nodes within the rich-club have revealed that neuronal growth factor (NGF), caspase 3, and H-Ras are the most central nodes and potentially an interesting targets for therapy. Our integrative approach uncovers the molecular architecture of SCI interactome, and provide an essential resource for evaluating significant therapeutic candidates.     

Structure-Based Network Analysis of Activation Mechanisms in the ErbB Family of Receptor Tyrosine Kinases: The Regulatory Spine Residues Are Global Mediators of Structural Stability and Allosteric Interactions

The ErbB protein tyrosine kinases are among the most important cell signaling families and mutation-induced modulation of their activity is associated with diverse functions in biological networks and human disease. We have combined molecular dynamics simulations of the ErbB kinases with the protein structure network modeling to characterize the reorganization of the residue interaction networks during conformational equilibrium changes in the normal and oncogenic forms. Structural stability and network analyses have identified local communities integrated around high centrality sites that correspond to the regulatory spine residues. This analysis has provided a quantitative insight to the mechanism of mutation-induced “superacceptor” activity in oncogenic EGFR dimers. We have found that kinase activation may be determined by allosteric interactions between modules of structurally stable residues that synchronize the dynamics in the nucleotide binding site and the αC-helix with the collective motions of the integrating αF-helix and the substrate binding site. The results of this study have pointed to a central role of the conserved His-Arg-Asp (HRD) motif in the catalytic loop and the Asp-Phe-Gly (DFG) motif as key mediators of structural stability and allosteric communications in the ErbB kinases. We have determined that residues that are indispensable for kinase regulation and catalysis often corresponded to the high centrality nodes within the protein structure network and could be distinguished by their unique network signatures. The optimal communication pathways are also controlled by these nodes and may ensure efficient allosteric signaling in the functional kinase state. Structure-based network analysis has quantified subtle effects of ATP binding on conformational dynamics and stability of the EGFR structures. Consistent with the NMR studies, we have found that nucleotide-induced modulation of the residue interaction networks is not limited to the ATP site, and may enhance allosteric cooperativity with the substrate binding region by increasing communication capabilities of mediating residues.     

Identification of significant pathways in gastric cancer based on protein-protein interaction networks and cluster analysis

Gastric cancer is one of the most common and lethal cancers worldwide. However, despite its clinical importance, the regulatory mechanisms involved in the aggressiveness of this cancer are still poorly understood. A better understanding of the biology, genetics and molecular mechanisms of gastric cancer would be useful in developing novel targeted approaches for treating this disease. In this study we used protein-protein interaction networks and cluster analysis to comprehensively investigate the cellular pathways involved in gastric cancer. A primary immunodeficiency pathway, focal adhesion, ECM-receptor interactions and the metabolism of xenobiotics by cytochrome P450 were identified as four important pathways associated with the progression of gastric cancer. The genes in these pathways, e.g., ZAP70, IGLL1, CD79A, COL6A3, COL3A1, COL1A1, CYP2C18 and CYP2C9, may be considered as potential therapeutic targets for gastric cancer.