Search

Background Identifying protein complexes from protein–protein interaction (PPI) networks is a crucial task, and many related algorithms have been developed to solve this issue. These algorithms usually consider a node’s direct neighbors and ignore resource allocation and second-order neighbors. The effective use of such information is crucial to protein complex detection.Results To overcome this deficiency, this paper proposes a new protein complex identification method based on node-local topological properties and gene expression information on a new weighted PPI network, named NLPGE-WPN (joint node-local topological properties and gene expression information on weighted PPI network). First, based on the resource allocation of the PPI network and gene expression, a new weight metric is designed to describe the interaction between proteins. Second, our method constructs a series of dense complex cores based on density and network diameter constraints; the final complexes are recognized by expanding the second-order neighbor nodes of core complexes. Experimental results demonstrate that this algorithm has improved the performances of precision and f-measure, which is more valid in identifying protein complexes．Conclusions This identification method is simple and can accurately identify more complexes by integrating node-local properties and gene expression on PPI weighted networks.

Sepsis-induced acute respiratory distress syndrome (ARDS) remains a major threat to human health without effective therapeutic drugs. Previous studies demonstrated the power of gene expression profiling to reveal pathological changes associated with sepsis-induced ARDS. However, there is still a lack of systematic data mining framework for identifying potential targets for treatment. In this study, we demonstrated the feasibility of druggable targets prediction based on gene expression data. Through the functional enrichment analysis of microarray-based expression profiles between sepsis-induced ARDS and non-sepsis ARDS samples, we revealed genes involved in anti-microbial infection immunity were significantly altered in sepsis-induced ARDS. Protein–protein interaction (PPI) network analysis highlighted TOP2A gene as the key regulator in the dysregulated gene network of sepsis-induced ARDS. We were also able to predict several therapeutic drug candidates for sepsis-induced ARDS using Connectivity Map (Cmap) database, among which doxorubicin was identified to interact with TOP2A with a high affinity similar to its endogenous ligand. Overall, our findings suggest that doxorubicin could be a potential therapeutic for sepsis-induced ARDS by targeting TOP2A, which requires further investigation and validation. The whole study relies on publicly available dataset and publicly accessible database or bioinformatic tools for data mining. Therefore, our study benchmarks a workflow for druggable target prediction which can be widely applicable in the search of targets in other pathological conditions., Graphical Abstract

Molecular glues are an intriguing therapeutic modality that harness small-molecules to induce interactions between proteins that typically do not interact. However, such molecules are rare and have been discovered fortuitously, thus limiting their potential as a general strategy for therapeutic intervention. We postulated that natural products bearing one or more electrophilic sites may be an unexplored source of new molecular glues, potentially acting through multi-covalent attachment. Using chemoproteomic platforms, we show that members of the manumycin family of polyketides, which bear multiple potentially reactive sites, target C374 of the putative E3 ligase UBR7 in breast cancer cells and engage in molecular glue interactions with the neo-substrate tumor-suppressor TP53, leading to p53 transcriptional activation and cell death. Our results reveal a novel anti-cancer mechanism of this natural product family and highlight the potential for combining chemoproteomics and multi-covalent natural products for the discovery of new molecular glues.

The genome of the SARS-CoV-2 Omicron variant (B.1.1.529) was released on November 22, 2021, which has caused a flurry of media attention due the large number of mutations it contains. These raw data have spurred questions around vaccine efficacy. Given that neither the structural information nor the experimentally-derived antibody interaction of this variant are available, we have turned to predictive computational methods to model the mutated structure of the spike protein's receptor binding domain and posit potential changes to vaccine efficacy. In this study, we predict some structural changes in the receptor-binding domain that may reduce antibody interaction without completely evading existing neutralizing antibodies (and therefore current vaccines).

Alzheimer’s Disease (AD) is a progressive neurodegenerative disorder and the most common cause of dementia; however, early diagnosis of the disease is challenging. Research suggests that biomarkers found in blood, such as microRNAs (miRNA), may be promising for AD diagnostics. Experimental data on miRNA–target interactions (MTI) associated with AD are scattered across databases and publications, thus making the identification of promising miRNA biomarkers for AD difficult. In response to this, a list of experimentally validated AD-associated MTIs was obtained from miRTarBase. Cytoscape was used to create a visual MTI network. STRING software was used for protein–protein interaction analysis and mirPath was used for pathway enrichment analysis. Several targets regulated by multiple miRNAs were identified, including: BACE1, APP, NCSTN, SP1, SIRT1, and PTEN. The miRNA with the highest numbers of interactions in the network were: miR-9, miR-16, miR-34a, miR-106a, miR-107, miR-125b, miR-146, and miR-181c. The analysis revealed seven subnetworks, representing disease modules which have a potential for further biomarker development. The obtained MTI network is not yet complete, and additional studies are needed for the comprehensive understanding of the AD-associated miRNA targetome.

Drug-likeness quantification is useful for screening drug candidates. Quantitative estimates of drug-likeness (QED) are commonly used to assess quantitative drug efficacy but are not suitable for screening compounds targeting protein-protein interactions (PPIs), which have recently gained attention. Therefore, we developed a quantitative estimate index for compounds targeting PPIs (QEPPI), specifically for early-stage screening of PPI-targeting compounds. QEPPI is an extension of the QED method for PPI-targeting drugs that models physicochemical properties based on the information available for drugs/compounds, specifically those reported to act on PPIs. FDA-approved drugs and compounds in iPPI-DB, which comprise PPI inhibitors and stabilizers, were evaluated using QEPPI. The results showed that QEPPI is more suitable than QED for early screening of PPI-targeting compounds. QEPPI was also considered an extended concept of the “Rule-of-Four” (RO4), a PPI inhibitor index. We evaluated the discriminatory performance of QEPPI and RO4 for datasets of PPI-target compounds and FDA-approved drugs using F-score and other indices. The F-scores of RO4 and QEPPI were 0.451 and 0.501, respectively. QEPPI showed better performance and enabled quantification of drug-likeness for early-stage PPI drug discovery. Hence, it can be used as an initial filter to efficiently screen PPI-targeting compounds.

Proteins rarely carry out their cellular functions in isolation. Instead, eukaryotic proteins engage in about six interactions with other proteins on average. The aggregated protein interactome of an organism forms a “hairy ball”-type protein-protein interaction (PPI) network. Yet, in a typical human cell, only about half of all proteins are expressed at a particular time. Hence, it has become common practice to prune the full PPI network to the subset of expressed proteins. If RNAseq data is available, one can further resolve the specific protein isoforms present in a cell or tissue. Here, we review various approaches, software tools and webservices that enable users to construct context-specific or tissue-specific PPI networks and how these are rewired between two cellular conditions. We illustrate their different functionalities on the example of the interactions involving the human TNR6 protein. In an outlook, we describe how PPI networks may be integrated with epigenetic data or with data on the activity of splicing factors.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), better known as COVID-19, has become a current threat to humanity. The second wave of the SARS-CoV-2 virus has hit many countries, and the confirmed COVID-19 cases are quickly spreading. Therefore, the epidemic is still passing the terrible stage. Having idiopathic pulmonary fibrosis (IPF) and chronic obstructive pulmonary disease (COPD) are the risk factors of the COVID-19, but the molecular mechanisms that underlie IPF, COPD, and CVOID-19 are not well understood. Therefore, we implemented transcriptomic analysis to detect common pathways and molecular biomarkers in IPF, COPD, and COVID-19 that help understand the linkage of SARS-CoV-2 to the IPF and COPD patients. Here, three RNA-seq datasets (GSE147507, GSE52463, and GSE57148) from Gene Expression Omnibus (GEO) is employed to detect mutual differentially expressed genes (DEGs) for IPF, and COPD patients with the COVID-19 infection for finding shared pathways and candidate drugs. A total of 65 common DEGs among these three datasets were identified. Various combinatorial statistical methods and bioinformatics tools were used to build the protein–protein interaction (PPI) and then identified Hub genes and essential modules from this PPI network. Moreover, we performed functional analysis under ontologies terms and pathway analysis and found that IPF and COPD have some shared links to the progression of COVID-19 infection. Transcription factors–genes interaction, protein–drug interactions, and DEGs-miRNAs coregulatory network with common DEGs also identified on the datasets. We think that the candidate drugs obtained by this study might be helpful for effective therapeutic in COVID-19.

14-3-3σ is an acidic homodimer protein with more than one hundred different protein partners associated with oncogenic signaling and cell cycle regulation. This review aims to highlight the crucial role of 14-3-3σ in controlling tumor growth and apoptosis and provide a detailed discussion on the structure–activity relationship and binding interactions of the most recent 14-3-3σ protein-protein interaction (PPI) modulators reported to date, which has not been reviewed previously. This includes the new fusicoccanes stabilizers (FC-NAc, DP-005), fragment stabilizers (TCF521-123, TCF521-129, AZ-003, AZ-008), phosphate-based inhibitors (IMP, PLP), peptide inhibitors (2a–d), as well as inhibitors from natural sources (85531185, 95911592). Additionally, this review will also include the discussions of the recent efforts by a different group of researchers for understanding the binding mechanisms of existing 14-3-3σ PPI modulators. The strategies and state-of-the-art techniques applied by various group of researchers in the discovery of a different chemical class of 14-3-3σ modulators for cancer are also briefly discussed in this review, which can be used as a guide in the development of new 14-3-3σ modulators in the near future.

Synovial sarcoma (SS) is a highly aggressive soft tissue tumor with high risk of local recurrence and metastasis. However, the mechanisms underlying SS metastasis are still largely unclear. The purpose of this study is to screen metastasis-associated biomarkers in SS by integrated bioinformatics analysis. Two mRNA datasets (GSE40018 and GSE40021) were selected to analyze the differentially expressed genes (DEGs). Using the Database for Annotation, Visualization and Integrated Discovery (DAVID) and gene set enrichment analysis (GSEA), functional and pathway enrichment analyses were performed for DEGs. Then, the protein-protein interaction (PPI) network was constructed via the Search Tool for the Retrieval of Interacting Genes (STRING) database. The module analysis of the PPI network and hub genes validation were performed using Cytoscape software. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of the hub genes were performed using WEB-based GEne SeT AnaLysis Toolkit (WebGestalt). The expression levels and survival analysis of hub genes were further assessed through Gene Expression Profiling Interactive Analysis (GEPIA) and the Kaplan-Meier plotter database. In total, 213 overlapping DEGs were identified, of which 109 were upregulated and 104 were downregulated. GO analysis revealed that the DEGs were predominantly involved in mitosis and cell division. KEGG pathways analysis demonstrated that most DEGs were significantly enriched in cell cycle pathway. GSEA revealed that the DEGs were mainly enriched in oocyte meiosis, cell cycle and DNA replication pathways. A key module was identified and 10 hub genes (CENPF, KIF11, KIF23, TTK, MKI67, TOP2A, CDC45, MELK, AURKB, and BUB1) were screened out. The expression and survival analysis disclosed that the 10 hub genes were upregulated in SS patients and could result in significantly reduced survival. Our study identified a series of metastasis-associated biomarkers involved in the progression of SS, and may provide novel therapeutic targets for SS metastasis.

Ganoderma lucidum has salutary effects on tumor treatment, including pancreatic cancer and hepatocellular carcinoma. However, the molecular mechanisms underlying Ganoderma lucidum therapy is obscure. In this study, the Hepa1-6-bearing C57 BL/6 mouse model was utilized to explore the therapeutic efficacy of Ganoderma lucidum extract (GLE), documenting that it could effectively inhibit tumor growth. The microRNA (miRNA) profiles of GLE-treated and untreated mice were detected, and 25 differentially expressed (DE) miRNAs were determined, including 24 up-expressed and one down-expressed miRNAs. Using the ClusterOne algorithm, 8 hub miRNAs were isolated from the established miRNA-target network. The qRT-PCR assay demonstrated that these 8 miRNAs were up-expressed in the GLE treated tumor mice. Furthermore, the mRNA profiles showed that there are 76 DE mRNAs between GLE treated and model groups. The protein-protein interaction (PPI) network shows that Cntn1, Irs1, Nfkbia, Rybp and Ywhaz playing important roles, and qRT-PCR further revealed they were down-expressed in GLE treated Hepa1-6-bearing C57 BL/6 mice. The rebuilt miRNA-target network was shown that these 5 mRNAs were regulated by mmu-mir-23a-5p, -3102-3p, -337-3p, and -467a-3p, respectively. This study suggested that these 4 interesting miRNAs were potential biomarkers for evaluation of GLE efficacy, which may down-regulate the expression of Cntn1, Irs1, Nfkbia, Rybp and Ywhaz, and mediate many signaling pathways occurring in tumor treatment.

Background The number of applications of deep learning algorithms in bioinformatics is increasing as they usually achieve superior performance over classical approaches, especially, when bigger training datasets are available. In deep learning applications, discrete data, e.g. words or n-grams in language, or amino acids or nucleotides in bioinformatics, are generally represented as a continuous vector through an embedding matrix. Recently, learning this embedding matrix directly from the data as part of the continuous iteration of the model to optimize the target prediction – a process called ‘end-to-end learning’ – has led to state-of-the-art results in many fields. Although usage of embeddings is well described in the bioinformatics literature, the potential of end-to-end learning for single amino acids, as compared to more classical manually-curated encoding strategies, has not been systematically addressed. To this end, we compared classical encoding matrices, namely one-hot, VHSE8 and BLOSUM62, to end-to-end learning of amino acid embeddings for two different prediction tasks using three widely used architectures, namely recurrent neural networks (RNN), convolutional neural networks (CNN), and the hybrid CNN-RNN. Results By using different deep learning architectures, we show that end-to-end learning is on par with classical encodings for embeddings of the same dimension even when limited training data is available, and might allow for a reduction in the embedding dimension without performance loss, which is critical when deploying the models to devices with limited computational capacities. We found that the embedding dimension is a major factor in controlling the model performance. Surprisingly, we observed that deep learning models are capable of learning from random vectors of appropriate dimension. Conclusion Our study shows that end-to-end learning is a flexible and powerful method for amino acid encoding. Further, due to the flexibility of deep learning systems, amino acid encoding schemes should be benchmarked against random vectors of the same dimension to disentangle the information content provided by the encoding scheme from the distinguishability effect provided by the scheme.

In silico binding site location and pose prediction for a molecule targeted at a large protein surface is a challenging task. We report a blind test with two peptidomimetic molecules that bind the flu virus hemagglutinin (HA) surface antigen, JNJ7918 and JNJ4796 (recently disclosed in van Dongen et al., Science, 2019, 363). Tests with a series of conventional approaches such as rigid (receptor) docking against available X-ray crystal structures or against an ensemble of structures generated by quick methodologies (NMA, homology modeling) gave mixed results, due to the shallowness and flexibility of the binding site and the sheer size of the target. However, tests with our Monte Carlo platform PELE in two protocols involving either exploration of the whole protein surface (global exploration), or the latter followed by refinement of best solutions (local exploration) yielded remarkably good results by locating the actual binding site and generating binding modes that recovered all native contacts found in the X-ray structures. Thus, the Monte Carlo scheme of PELE seems promising as a quick methodology to overcome the challenge of identifying entirely unknown binding sites and modes for protein–protein disruptors.

Microarray data of non-functioning pituitary adenomas (NFPAs) were analyzed to disclose novel genes and pathways involved in NFPA tumorigenesis. Raw microarray data were downloaded from Gene Expression Omnibus. Data pre-treatment and differential analysis were conducted using packages in R. Functional and pathway enrichment analyses were performed using package GOs-tats. A protein-protein interaction (PPI) network was constructed using server STRING and Cytoscape. Known genes involved in pituitary adenomas (PAs), were obtained from the Comparative Toxicogenomics Database. A total of 604 differentially expressed genes (DEGs) were identifed between NFPAs and controls, including 177 up- and 427 down-regulated genes. Jak-STAT and p53 signaling pathways were significantly enriched by DEGs. The PPI network of DEGs was constructed, containing 99 up- and 288 down-regulated known disease genes (e.g. EGFR and ESR1) as well as 16 up- and 17 down-regulated potential novel NFPAs-related genes (e.g. COL4A5, LHX3, MSN, and GHSR). Genes like COL4A5, LHX3, MSN, and GHSR and pathways such as p53 signaling and Jak-STAT signaling, might participate in NFPA development. Although further validations are required, these findings might provide guidance for future basic and therapy researches.

Our study's goal was to screen novel biomarkers that could accurately predict the progression and prognosis of bladder cancer (BC). Firstly, we used the Gene Expression Omnibus (GEO) dataset GSE37815 to screen differentially expressed genes (DEGs). Secondly, we used the DEGs to construct a co-expression network by weighted gene co-expression network analysis (WGCNA) in GSE71576. We then screened the brown module, which was significantly correlated with the histologic grade (r = 0.85, p = 1e-12) of BC. We conducted functional annotation on all genes of the brown module and found that the genes of the brown module were mainly significantly enriched in "cell cycle" correlation pathways. Next, we screened out two real hub genes (ANLN, HMMR) by combining WGCNA, protein-protein interaction (PPI) network and survival analysis. Finally, we combined the GEO datasets (GSE13507, GSE37815, GSE31684, GSE71576). Oncomine, Human Protein Atlas (HPA), and The Cancer Genome Atlas (TCGA) dataset to confirm the predict value of the real hub genes for BC progression and prognosis. A gene-set enrichment analysis (GSEA) revealed that the real hub genes were mainly enriched in "bladder cancer" and "cell cycle" pathways. A survival analysis showed that they were of great significance in predicting the prognosis of BC. In summary, our study screened and confirmed that two biomarkers could accurately predict the progression and prognosis of BC, which is of great significance for both stratification therapy and the mechanism study of BC.

High-throughput technology has generated large-scale protein interaction data, which is crucial in our understanding of biological organisms. Many complex identification algorithms have been developed to determine protein complexes. However, these methods are only suitable for dense protein interaction networks, because their capabilities decrease rapidly when applied to sparse protein&ndash, protein interaction (PPI) networks. In this study, based on penalized matrix decomposition (PMD), a novel method of penalized matrix decomposition for the identification of protein complexes (i.e., PMDpc) was developed to detect protein complexes in the human protein interaction network. This method mainly consists of three steps. First, the adjacent matrix of the protein interaction network is normalized. Second, the normalized matrix is decomposed into three factor matrices. The PMDpc method can detect protein complexes in sparse PPI networks by imposing appropriate constraints on factor matrices. Finally, the results of our method are compared with those of other methods in human PPI network. Experimental results show that our method can not only outperform classical algorithms, such as CFinder, ClusterONE, RRW, HC-PIN, and PCE-FR, but can also achieve an ideal overall performance in terms of a composite score consisting of F-measure, accuracy (ACC), and the maximum matching ratio (MMR).

International audience; MADS family transcription factors are crucial during plant reproductive development, and have evolved a complex protein–protein interaction (PPI) network. Proteins of the SEPALLATA (SEP) clade are required for tetramer formation and can act as critical ‘hubs’ in the network. Rümpler et al. (2018) have now provided quantitative measures of the contribution of individual amino acids to cooperative DNA binding, laying a foundation for predicting MADS tetramer formation based on primary sequence. It is an important step forward in understanding how cooperativity affects processes from flowering time to floral organ identity

Background Budding yeast, S. cerevisiae, has been used extensively as a model organism for studying cellular processes in evolutionarily distant species, including humans. However, different human tissues, while inheriting a similar genetic code, exhibit distinct anatomical and physiological properties. Specific biochemical processes and associated biomolecules that differentiate various tissues are not completely understood, neither is the extent to which a unicellular organism, such as yeast, can be used to model these processes within each tissue. Results We present a novel framework to systematically quantify the suitability of yeast as a model organism for different human tissues. To this end, we develop a computational method for dissecting the global human interactome into tissue-specific cellular networks. By individually aligning these networks with the yeast interactome, we simultaneously partition the functional space of human genes, and their corresponding pathways, based on their conservation both across species and among different tissues. Finally, we couple our framework with a novel statistical model to assess the conservation of tissue-specific pathways and infer the overall similarity of each tissue with yeast. We further study each of these subspaces in detail, and shed light on their unique biological roles in the human tissues. Conclusions Our framework provides a novel tool that can be used to assess the suitability of the yeast model for studying tissue-specific physiology and pathophysiology in humans. Many complex disorders are driven by a coupling of housekeeping (universally expressed in all tissues) and tissue-selective (expressed only in specific tissues) dysregulated pathways. While tissue-selective genes are significantly associated with the onset and development of a number of tissue-specific pathologies, we show that the human-specific subset has even higher association. Consequently, they provide excellent candidates as drug targets for therapeutic interventions. Electronic supplementary material The online version of this article (doi:10.1186/s12918-015-0253-0) contains supplementary material, which is available to authorized users.

Type VI protein secretion systems (T6SSs) are specialized transport apparatus which can target both eukaryotic and prokaryotic cells and play key roles in host–pathogen–microbiota interactions. Therefore, T6SSs have attracted much attention as a research topic during the past ten years. In this review, we particularly summarized the T6SS antibacterial function, which involves an interesting offensive and defensive mechanism of the effector–immunity (E–I) pairs. The three main categories of effectors that target the cell wall, membranes, and nucleic acids during bacterial interaction, along with their corresponding immunity proteins are presented. We also discuss structural analyses of several effectors and E–I pairs, which explain the offensive and defensive mechanisms underpinning T6SS function during bacterial competition for niche-space, as well as the bioinformatics, proteomics, and protein–protein interaction (PPI) methods used to identify and characterize T6SS mediated E–I pairs. Additionally, we described PPI methods for verifying E–I pairs.

Background: H. sapiens-M. tuberculosis H37Rv protein-protein interaction (PPI) data are very important information to illuminate the infection mechanism of M. tuberculosis H37Rv. But current H. sapiens-M. tuberculosis H37Rv PPI data are very scarce. This seriously limits the study of the interaction between this important pathogen and its host H. sapiens. Computational prediction of H. sapiens-M. tuberculosis H37Rv PPIs is an important strategy to fill in the gap. Domain-domain interaction (DDI) based prediction is one of the frequently used computational approaches in predicting both intra-species and inter-species PPIs. However, the performance of DDI-based host-pathogen PPI prediction has been rather limited. Results: We develop a stringent DDI-based prediction approach with emphasis on (i) differences between the specific domain sequences on annotated regions of proteins under the same domain ID and (ii) calculation of the interaction strength of predicted PPIs based on the interacting residues in their interaction interfaces. We compare our stringent DDI-based approach to a conventional DDI-based approach for predicting PPIs based on gold standard intra-species PPIs and coherent informative Gene Ontology terms assessment. The assessment results show that our stringent DDI-based approach achieves much better performance in predicting PPIs than the conventional approach. Using our stringent DDI-based approach, we have predicted a small set of reliable H. sapiens-M. tuberculosis H37Rv PPIs which could be very useful for a variety of related studies. We also analyze the H. sapiens-M. tuberculosis H37Rv PPIs predicted by our stringent DDI-based approach using cellular compartment distribution analysis, functional category enrichment analysis and pathway enrichment analysis. The analyses support the validity of our prediction result. Also, based on an analysis of the H. sapiens-M. tuberculosis H37Rv PPI network predicted by our stringent DDI-based approach, we have discovered some important properties of domains involved in host-pathogen PPIs. We find that both host and pathogen proteins involved in host-pathogen PPIs tend to have more domains than proteins involved in intra-species PPIs, and these domains have more interaction partners than domains on proteins involved in intra-species PPI. Conclusions: The stringent DDI-based prediction approach reported in this work provides a stringent strategy for predicting host-pathogen PPIs. It also performs better than a conventional DDI-based approach in predicting PPIs. We have predicted a small set of accurate H. sapiens-M. tuberculosis H37Rv PPIs which could be very useful for a variety of related studies.

Smoking is the primary cause of lung cancer and is linked to 85% of lung cancer cases. However, how lung cancer develops in patients with smoking history remains unclear. Systems approaches that combine human protein-protein interaction (PPI) networks and gene expression data are superior to traditional methods. We performed these systems to determine the role that smoking plays in lung cancer development and used the support vector machine (SVM) model to predict PPIs. By defining expression variance (EV), we found 520 dynamic proteins (EV>0.4) using data from the Human Protein Reference Database and Gene Expression Omnibus Database, and built 7 dynamic PPI subnetworks of lung cancer in patients with smoking history. We also determined the primary functions of each subnetwork: signal transduction, apoptosis, and cell migration and adhesion for subnetwork A; cell-sustained angiogenesis for subnetwork B; apoptosis for subnetwork C; and, finally, signal transduction and cell replication and proliferation for subnetworks D–G. The probability distribution of the degree of dynamic protein and static protein differed, clearly showing that the dynamic proteins were not the core proteins which widely connected with their neighbor proteins. There were high correlations among the dynamic proteins, suggesting that the dynamic proteins tend to form specific dynamic modules. We also found that the dynamic proteins were only correlated with the expression of selected proteins but not all neighbor proteins when cancer occurred.

Background Essential hypertension (EH) is a complex disease as a consequence of interaction between environmental factors and genetic background, but the pathogenesis of EH remains elusive. The emerging tools of network medicine offer a platform to explore a complex disease at system level. In this study, we aimed to identify the key proteins and the biological regulatory pathways involving in EH and further to explore the molecular connectivities between these pathways by the topological analysis of the Protein-protein interaction (PPI) network. Result The extended network including one giant network consisted of 535 nodes connected via 2572 edges and two separated small networks. 27 proteins with high BC and 28 proteins with large degree have been identified. NOS3 with highest BC and Closeness centrality located in the centre of the network. The backbone network derived from high BC proteins presents a clear and visual overview which shows all important regulatory pathways for blood pressure (BP) and the crosstalk between them. Finally, the robustness of NOS3 as central protein and accuracy of backbone were validated by 287 test networks. Conclusion Our finding suggests that blood pressure variation is orchestrated by an integrated PPI network centered on NOS3.