Protein Interactome Research Articles

Mapping the Human Cell Surface Interactome: A Key to Decode Cell-to-Cell Communication.

Proteins on the surfaces of cells serve as physical connection points to bridge one cell with another, enabling direct communication between cells and cohesive structure. As biomedical research makes the leap from characterizing individual cells toward understanding the multicellular organization of the human body, the binding interactions between molecules on the surfaces of cells are foundational both for computational models and for clinical efforts to exploit these influential receptor pathways. To achieve this grander vision, we must assemble the full interactome of ways surface proteins can link together. This review investigates how close we are to knowing the human cell surface protein interactome. We summarize the current state of databases and systematic technologies to assemble surface protein interactomes, while highlighting substantial gaps that remain. We aim for this to serve as a road map for eventually building a more robust picture of the human cell surface protein interactome.

Drug repositioning for rosacea disease: Biological TARGET identification, molecular docking, pharmacophore mapping, and molecular dynamics analysis

Rosacea is a chronic dermatological condition that currently lacks a clear treatment approach due to an uncomprehensive knowledge of its pathogenesis. The main obstacle lies in understanding its etiology and the mode of action of the different drugs used. This study aims to clarify these aspects by employing drug repositioning. Using an in silico approach, we performed a transcriptomic analysis comparing samples from individuals with diverse types of rosacea to those from healthy controls to identify genes deregulated in this disease. Subsequently, we realized molecular docking and molecular dynamics studies to assess the binding affinity of drugs currently used to treat rosacea and drugs that target proteins interacting with, and thus affecting, proteins deregulated in rosacea. Our findings revealed that the downregulation of SKAP2 and upregulation of S100A7A in rosacea, could be involved in the pathogenesis of the disease. Furthermore, considering the drugs currently used for rosacea management, we demonstrated stable interactions between isotretinoin and BFH772 with SKAP2, and permethrin and PAC-14028 with S100A7A. Similarly, considering drugs targeting SKAP2 and S100A7A interactome proteins, we found that pitavastatin and dasatinib exert stable interactions with SKAP2, and lovastatin and tirbanibulin with S100A7A. In addition, we determine that the types of bonds involved in the interactions were different in SKAP2 from S100A7A. The drug-SKAP2 interactions are hydrogen bonds, whereas the drug-S100A7A interactions are of the hydrophobic type. In conclusion, our study provides evidence for the possible contribution of SKAP2 and S100A7A to rosacea pathology. Furthermore, it provides significant information on the molecular interactions between drugs and these proteins, highlighting the importance of considering structural features and binding interactions in the design of targeted therapies for skin disorders such as rosacea.

When less is more - a fast TurboID knock-in approach for high-sensitivity endogenous interactome mapping.

In recent years, proximity labeling has established itself as an unbiased and powerful approach to map the interactome of specific proteins. Although physiological expression of labeling enzymes is beneficial for the mapping of interactors, generation of the desired cell lines remains time-consuming and challenging. Using our established pipeline for rapid generation of C- and N-terminal CRISPR-Cas9 knock-ins (KIs) based on antibiotic selection, we were able to compare the performance of commonly used labeling enzymes when endogenously expressed. Endogenous tagging of the µ subunit of the adaptor protein (AP)-1 complex with TurboID allowed identification of known interactors and cargo proteins that simple overexpression of a labeling enzyme fusion protein could not reveal. We used the KI strategy to compare the interactome of the different AP complexes and clathrin and were able to assemble lists of potential interactors and cargo proteins that are specific for each sorting pathway. Our approach greatly simplifies the execution of proximity labeling experiments for proteins in their native cellular environment and allows going from CRISPR transfection to mass spectrometry analysis and interactome data in just over a month.

Proximity labelling of pro-interleukin-1α reveals evolutionary conserved nuclear interactions.

Interleukin-1α is a suggested dual-function cytokine that diverged from interleukin-1β in mammals potentially by acquiring additional biological roles that relate to highly conserved regions in the pro-domain of interleukin-1α, including a nuclear localisation sequence and histone acetyltransferase-binding domains. Why evolution modified pro-interleukin-1α's subcellular location and protein interactome, and how this shaped interleukin-1α's intracellular role, is unknown. Here we show that TurboID proximity labelling with pro-interleukin-1α suggests a nuclear role for pro-interleukin-1α that involves interaction with histone acetyltransferases, including EP300. We also identify and validate inactivating mutations in the pro-interleukin-1α nuclear localisation sequence of multiple mammalian species, including toothed whales, castorimorpha and marsupials. However, histone acetyltransferase-binding domains are conserved in those species that have lost pro-interleukin-1α nuclear localisation. Together, these data suggest that histone acetyltransferase binding and nuclear localisation occurred together, and that while some species lost the nuclear localisation sequence in their pro-interleukin-1α, histone acetyltransferase binding ability was maintained. The nuclear localisation sequence was lost from several distinct species at different evolutionary times, suggesting convergent evolution, and that the loss of the nuclear localisation sequence confers some important biological outcome.

Transcriptomics and weighted protein network analyses of the LRRK2 protein interactome reveal distinct molecular signatures for sporadic and LRRK2 Parkinson’s Disease

Mutations in the LRRK2 gene are the most common genetic cause of familial Parkinson’s Disease (LRRK2-PD) and an important risk factor for sporadic PD (sPD). Multiple clinical trials are ongoing to evaluate the benefits associated with the therapeutical reduction of LRRK2 kinase activity. In this study, we described the changes of transcriptomic profiles (whole blood mRNA levels) of LRRK2 protein interactors in sPD and LRRK2-PD cases as compared to healthy controls with the aim of comparing the two PD conditions. We went on to model the protein-protein interaction (PPI) network centred on LRRK2, which was weighted to reflect the transcriptomic changes on expression and co-expression levels of LRRK2 protein interactors. Our results showed that LRRK2 interactors present both similar and distinct alterations in expression levels and co-expression behaviours in the sPD and LRRK2-PD cases; suggesting that, albeit being classified as the same disease based on clinical features, LRRK2-PD and sPD display significant differences from a molecular perspective. Interestingly, the similar changes across the two PD conditions result in decreased connectivity within a topological cluster of the LRRK2 PPI network associated with protein metabolism/biosynthesis and ribosomal metabolism suggesting protein homoeostasis and ribosomal dynamics might be affected in both sporadic and familial PD in comparison with controls.

Trans-Interacting Plasma Membrane Proteins and Binding Partner Identification.

Plasma membrane proteins (PMPs) play critical roles in a myriad of physiological and disease conditions. A unique subset of PMPs functions through interacting with each other in trans at the interface between two contacting cells. These trans-interacting PMPs (tiPMPs) include adhesion molecules and ligands/receptors that facilitate cell-cell contact and direct communication between cells. Among the tiPMPs, a significant number have apparent extracellular binding domains but remain orphans with no known binding partners. Identification of their potential binding partners is therefore important for the understanding of processes such as organismal development and immune cell activation. While a number of methods have been developed for the identification of protein binding partners in general, very few are applicable to tiPMPs, which interact in a two-dimensional fashion with low intrinsic binding affinities. In this review, we present the significance of tiPMP interactions, the challenges of identifying binding partners for tiPMPs, and the landscape of method development. We describe current avidity-based screening approaches for identifying novel tiPMP binding partners and discuss their advantages and limitations. We conclude by highlighting the importance of developing novel methods of identifying new tiPMP interactions for deciphering the complex protein interactome and developing targeted therapeutics for diseases.

Identification of Host Proteins Involved in Hepatitis B Virus Genome Packaging.

A critical part of the hepatitis B virus (HBV) life cycle is the packaging of the pregenomic RNA (pgRNA) into nucleocapsids. While this process is known to involve several viral elements, much less is known about the identities and roles of host proteins in this process. To better understand the role of host proteins, we isolated pgRNA and characterized its protein interactome in cells expressing either packaging-competent or packaging-incompetent HBV genomes. We identified over 250 host proteins preferentially associated with pgRNA from the packaging-competent version of the virus. These included proteins already known to support capsid formation, enhance viral gene expression, catalyze nucleocapsid dephosphorylation, and bind to the viral genome, demonstrating the ability of the approach to effectively reveal functionally significant host-virus interactors. Three of these host proteins, AURKA, YTHDF2, and ATR, were selected for follow-up analysis. RNA immunoprecipitation qPCR (RIP-qPCR) confirmed pgRNA-protein association in cells, and siRNA knockdown of the proteins showed decreased encapsidation efficiency. This study provides a template for the use of comparative RNA-protein interactome analysis in conjunction with virus engineering to reveal functionally significant host-virus interactions.

Glucocorticoid receptor controls atopic dermatitis inflammation via functional interactions with P63 and autocrine signaling in epidermal keratinocytes

Atopic dermatitis (AD), a prevalent chronic inflammatory disease with multifactorial etiology, features epidermal barrier defects and immune overactivation. Synthetic glucocorticoids (GCs) are widely prescribed for treating AD due to their anti-inflammatory actions; however, mechanisms are incompletely understood. Defective local GC signaling due to decreased production of endogenous ligand and/or GC receptor (GR) levels was reported in prevalent inflammatory skin disorders; whether this is a consequence or contributing factor to AD pathology is unclear. To identify the chromatin-bound cell-type-specific GR protein interactome in keratinocytes, we used rapid immunoprecipitation of endogenous proteins and mass spectrometry identifying 145 interactors that increased upon dexamethasone treatment. GR-interacting proteins were enriched in p53/p63 signaling, including epidermal transcription factors with critical roles in AD pathology. Previous analyses indicating mirrored AD-like phenotypes between P63 overexpression and GR loss in epidermis, and our data show an intricate relationship between these transcription factors in human keratinocytes, identifying TP63 as a direct GR target. Dexamethasone treatment counteracted transcriptional up-regulation of inflammatory markers by IL4/IL13, known to mimic AD, causing opposite shifts in GR and P63 genomic binding. Indeed, IL4/IL13 decreased GR and increased P63 levels in cultured keratinocytes and human epidermal equivalents (HEE), consistent with GR down-regulation and increased P63 expression in AD lesions vs normal skin. Moreover, GR knockdown (GRKD) resulted in constitutive increases in P63, phospho-P38 and S100A9, IL6, and IL33. Also, GRKD culture supernatants showed increased autocrine production of TH2-/TH1-/TH17-TH22-associated factors including IL4, CXCL10, CXCL11, and CXCL8. GRKD HEEs showed AD-like features including hyperplasia and abnormal differentiation, resembling phenotypes observed with GR antagonist or IL4/IL13 treatment. The simultaneous GR/P63 knockdown partially reversed constitutive up-regulation of inflammatory genes in GRKD. In summary, our data support a causative role for GR loss in AD pathogenesis via functional interactions with P63 and autocrine signaling in epidermal keratinocytes.

Nematode Extracellular Protein Interactome Expands Connections between Signaling Pathways.

Multicellularity was accompanied by the emergence of new classes of cell surface and secreted proteins. The nematode C. elegans is a favorable model to study cell surface interactomes, given its well-defined and stereotyped cell types and intercellular contacts. Here we report our C. elegans extracellular interactome dataset, the largest yet for an invertebrate. Most of these interactions were unknown, despite recent datasets for flies and humans, as our collection contains a larger selection of protein families. We uncover new interactions for all four major axon guidance pathways, including ectodomain interactions between three of the pathways. We demonstrate that a protein family known to maintain axon locations are secreted receptors for insulins. We reveal novel interactions of cystine-knot proteins with putative signaling receptors, which may extend the study of neurotrophins and growth-factor-mediated functions to nematodes. Finally, our dataset provides insights into human disease mechanisms and how extracellular interactions may help establish connectomes.

Comparative Proteomics and Interactome Analysis of the SARS-CoV-2 Nucleocapsid Protein in Human and Bat Cell Lines.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of COVID-19 and responsible for the global coronavirus pandemic which started in 2019. Despite exhaustive efforts to trace its origins, including potential links with pangolins and bats, the precise origins of the virus remain unclear. Bats have been recognized as natural hosts for various coronaviruses, including the Middle East respiratory coronavirus (MERS-CoV) and the SARS-CoV. This study presents a comparative analysis of the SARS-CoV-2 nucleocapsid protein (N) interactome in human and bat cell lines. We identified approximately 168 cellular proteins as interacting partners of SARS-CoV-2 N in human cells and 196 cellular proteins as interacting partners with this protein in bat cells. The results highlight pathways and events that are both common and unique to either bat or human cells. Understanding these interactions is crucial to comprehend the reasons behind the remarkable resilience of bats to viral infections. This study provides a foundation for a deeper understanding of host-virus interactions in different reservoirs.

Evidence of mitochondria origin of SARS-CoV-2 double-membrane vesicles: a review.

Coronavirus Disease-19 (COVID-19) pandemic is caused by SARS-CoV-2 that has infected more than 600 million people and killed more than 6 million people worldwide. This infection affects mainly certain groups of people that have high susceptibility to present severe COVID-19 due to comorbidities. Moreover, the long-COVID-19 comprises a series of symptoms that may remain in some patients for months after infection that further compromises their health. Thus, since this pandemic is profoundly affecting health, economy, and social life of societies, a deeper understanding of viral replication cycle could help to envisage novel therapeutic alternatives that limit or stop COVID-19. Several findings have unexpectedly discovered that mitochondria play a critical role in SARS-CoV-2 cell infection. Indeed, it has been suggested that this organelle could be the origin of its replication niches, the double membrane vesicles (DMV). In this regard, mitochondria derived vesicles (MDV), involved in mitochondria quality control, discovered almost 15 years ago, comprise a subpopulation characterized by a double membrane. MDV shedding is induced by mitochondrial stress, and it has a fast assembly dynamic, reason that perhaps has precluded their identification in electron microscopy or tomography studies. These and other features of MDV together with recent SARS-CoV-2 protein interactome and other findings link SARS-CoV-2 to mitochondria and support that these vesicles are the precursors of SARS-CoV-2 induced DMV. In this work, the morphological, biochemical, molecular, and cellular evidence that supports this hypothesis is reviewed and integrated into the current model of SARS-CoV-2 cell infection. In this scheme, some relevant questions are raised as pending topics for research that would help in the near future to test this hypothesis. The intention of this work is to provide a novel framework that could open new possibilities to tackle SARS-CoV-2 pandemic through mitochondria and DMV targeted therapies.

Residues in the fructose-binding pocket are required for ketohexokinase-A activity

Excessive fructose consumption is a primary contributor to the global surges in obesity, cancer, and metabolic syndrome. Fructolysis is not robustly regulated and is initiated by ketohexokinase (KHK). In this study, we determined the crystal structure of KHK-A, one of two human isozymes of KHK, in the apo-state at 1.85 Å resolution, and we investigated the roles of residues in the fructose-binding pocket by mutational analysis. Introducing alanine at D15, N42, or N45 inactivated KHK-A, whereas mutating R141 or K174 reduced activity and thermodynamic stability. Kinetic studies revealed that the R141A and K174A mutations reduced fructose affinity by 2- to 4-fold compared to WT KHK-A, without affecting ATP affinity. Molecular dynamics simulations provided mechanistic insights into the potential roles of the mutated residues in ligand coordination and the maintenance of an open state in one monomer and a closed state in the other. Protein–protein interactome analysis indicated distinct expression patterns and downregulation of partner proteins in different tumor tissues, warranting a re-evaluation of KHK’s role in cancer development and progression. The connections between different cancer genes and the KHK signaling pathway suggest that KHK is a potential target for preventing cancer metastasis. This study enhances our understanding of KHK-A's structure and function and offers valuable insights into potential targets for developing treatments for obesity, cancer, and metabolic syndrome.

Comparing the BAD Protein Interactomes in 2D and 3D Cell Culture Using Proximity Labeling.

Protein-protein interaction studies using proximity labeling techniques, such as biotin ligase-based BioID, have become integral in understanding cellular processes. Most studies utilize conventional 2D cell culture systems, potentially missing important differences in protein behavior found in 3D tissues. In this study, we investigated the protein-protein interactions of a protein, Bcl-2 Agonist of cell death (BAD), and compared conventional 2D culture conditions to a 3D system, wherein cells were embedded within a 3D extracellular matrix (ECM) mimic. Using BAD fused to the engineered biotin ligase miniTurbo (BirA*), we identified both overlapping and distinct BAD interactomes under 2D and 3D conditions. The known BAD binding proteins 14-3-3 isoforms and Bcl-XL interacted with BAD in both 2D and 3D. Of the 131 BAD-interactors identified, 56% were specific to 2D, 14% were specific to 3D, and 30% were common to both conditions. Interaction network analysis demonstrated differential associations between 2D and 3D interactomes, emphasizing the impact of the culture conditions on protein interactions. The 2D-3D overlap interactome encapsulated the apoptotic program, which is a well-known role of BAD. The 3D unique pathways were enriched in ECM signaling, suggestive of hitherto unknown functions for BAD. Thus, exploring protein-protein interactions in 3D provides novel clues into cell behavior. This exciting approach has the potential to bridge the knowledge gap between tractable 2D cell culture and organoid-like 3D systems.

Next-generation Drosophila protein interactome map and its functional implications

We describe a next-generation Drosophila protein interaction map-"DPIM2"-established from affinity purification-mass spectrometry of 5,805 baits, covering the largest fraction of the Drosophila proteome. The network contains 32,668 interactions among 3,644 proteins, organized into 632 clusters representing putative functional modules. Our analysis expands the pool of known protein interactions in Drosophila, provides annotation for poorly studied genes, and postulates previously undescribed protein interaction relationships. The predictive power and functional relevance of this network are probed through the lens of the Notch signaling pathway, and we find that newly identified members of complexes that include known Notch modifiers can also modulate Notch signaling. DPIM2 allows direct comparisons with a recently published human protein interaction network, defining the existence of functional interactions conserved across species. Thus, DPIM2 defines a valuable resource for predicting protein co-complex memberships and functional associations as well as generates functional hypotheses regarding specific protein interactions.

Bidirectional linkage of DNA barcodes for the multiplexed mapping of higher-order protein interactions in cells.

Capturing the full complexity of the diverse hierarchical interactions in the protein interactome is challenging. Here we report a DNA-barcoding method for the multiplexed mapping of pairwise and higher-order protein interactions and their dynamics within cells. The method leverages antibodies conjugated with barcoded DNA strands that can bidirectionally hybridize and covalently link to linearize closely spaced interactions within individual 3D protein complexes, encoding and decoding the protein constituents and the interactions among them. By mapping protein interactions in cancer cells and normal cells, we found that tumour cells exhibit a larger diversity and abundance of protein complexes with higher-order interactions. In biopsies of human breast-cancer tissue, the method accurately identified the cancer subtype and revealed that higher-order protein interactions are associated with cancer aggressiveness.

Understanding Longevity: SIN-3 and DAF-16 Revealed as Independent Players in Lifespan Regulation.

Aging is the process of gradual physio-biochemical deterioration. Although aging is inevitable, healthy aging is the key to individual and communal well-being. Therefore, it is essential to understand the regulation of aging. SIN-3/Sin-3 is a unique regulatory protein that regulates aging without DNA-binding activity. It functions by establishing multiple protein interactions. To understand the functional mechanism of this transcriptional regulator, the Caenorhabditis elegans protein interactome was assessed for SIN-3 interactions. DAF-16/FOXO emerged as one of the leading contenders for SIN-3-mediated regulation of aging. This study looks at the concerted role of SIN-3 and DAF-16 proteins in lifespan regulation. Phenotypic profiling for the mutants of these genes shows the functional accord between these 2 proteins with similar functions in stress response and vital biological processes. However, there were no significant physical interactions when checked for protein-protein interaction between SIN-3 and DAF-16 proteins. C. elegans genomics and transcriptomics data also indicated the possibilities of concerted gene regulation. This genetic regulation is more likely related to SIN-3 dominance on DAF-16 function. Overall, SIN-3 and DAF-16 proteins have strong functional interactions that ensure healthy aging. The influence of SIN-3 on DAF-16-mediated stress response is one of their convergence points in longevity regulation.

The Human Ganglioside Interactome in Live Cells Revealed Using Clickable Photoaffinity Ganglioside Probes.

Gangliosides, sialic acid bearing glycosphingolipids, are components of the outer leaflet of plasma membranes of all vertebrate cells. They contribute to cell regulation by interacting with proteins in their own membranes (cis) or their extracellular milieu (trans). As amphipathic membrane constituents, gangliosides present challenges for identifying their ganglioside protein interactome. To meet these challenges, we synthesized bifunctional clickable photoaffinity gangliosides, delivered them to plasma membranes of cultured cells, then captured and identified their interactomes using proteomic mass spectrometry. Installing probes on ganglioside lipid and glycan moieties, we captured cis and trans ganglioside-protein interactions. Ganglioside interactomes varied with the ganglioside structure, cell type, and site of the probe (lipid or glycan). Gene ontology revealed that gangliosides engage with transmembrane transporters and cell adhesion proteins including integrins, cadherins, and laminins. The approach developed is applicable to other gangliosides and cell types, promising to provide insights into molecular and cellular regulation by gangliosides.

DIPG-64. UNRAVELING H3K27M INTERACTIONS IN DIFFUSE MIDLINE GLIOMA

Abstract BACKGROUND Diffuse midline glioma (DMG) harboring the histone H3 K27M mutation, also known as diffuse intrinsic pontine glioma when it occurs in the pons, is a uniformly fatal pediatric brain tumor. We sought to better understand protein-protein interactions of the mutant histone in DMG. METHODS To address this, we carried out endogenous co-immunoprecipitation of histone H3 K27M from patient-derived DMG neurosphere cell lines, followed by LC-MS/MS identification of binding partners, and validation of binding by Western Blot. Moving beyond identification, we investigated the functional consequences of these interactions. RESULTS We identified putative novel interactors of H3 K27M with notable binding partners including histone chaperones and RNA helicases of the DDX family. We then explored the disruption of downstream effects resulting from this interaction by pharmacologically inhibiting DDX3 using the small molecule inhibitor RK-33, and assessed impact on viability and radiation sensitivity. CONCLUSIONS Understanding the protein interactome of the mutant histone can inform alternative ways for targeting DMG, thus enriching the therapeutic toolbox against this aggressive tumor. Functional characterization of binding partners holds the potential to identify novel therapeutic targets that cooperate with H3 K27M. The use of inhibitors targeting these partners could specifically destabilize oncohistone-dependent homeostasis and survival pathways and increase sensitivity to radiation and other therapies.

MDB-06. ALTERNATIVE SPLICING MEDIATED BY OTX2 DRIVES THE DEVELOPMENT OF GROUP 3 MEDULLOBLASTOMA

Abstract BACKGROUND OTX2 is a transcription factor and known molecular hallmark of Group 3 medulloblastoma (MB). While previous studies have focused on its association with active enhancer elements to modulate the chromatin landscape, here, we demonstrate an unexpected role for OTX2 in regulating Group 3 MB alternative splicing (AS). METHODS We used TurboID to broadly interrogate the Group 3 MB OTX2 protein interactome and identify associations with splicing regulators. RNA-sequencing was performed to identify OTX2-driven AS events, a subset of which were validated in Group 3 and Group 4 MB patient tumors and the developing human cerebellar rhombic lip (RL). Splice blocking morpholinos were used to mimic OTX2-regulated splicing events and to determine their role in MB cell fate decisions both in vitro and in vivo. RESULTS OTX2 associates with the Large Assembly of Splicing Regulators (LASR) complex through protein-protein interactions and modulates the localization of a LASR protein to help stabilize the complex and sustain a Group 3 MB splicing program. Motif enrichment analysis demonstrates that members of the LASR complex are also significantly enriched near OTX2-regulated alternatively spliced exons. OTX2 is found on nascent mRNA, and its role in AS is independent of DNA-regulatory functions. OTX2 controls a pro-tumorigenic splicing program that is mirrored in the human cerebellar RL origins. Among the OTX2-regulated spliced genes, PPHLN1 is expressed in the most primitive rhombic lip stem cells, and targeting PPHLN1 splicing reduces tumor growth, while enhancing differentiation and survival in vivo. CONCLUSIONS While splicing defects are commonly observed in adult cancers, the functional role of AS in the most aggressive MBs has never been explored. These findings identify OTX2-mediated AS as a major determinant of cell fate decisions and highlight AS events as novel therapeutic targets to drive Group 3 MB differentiation.

Predicting herb-disease associations using network-based measures in human protein interactome

BackgroundNatural herbs are frequently used to treat diseases or to relieve symptoms in many countries. Moreover, as their safety has been proven for a long time, they are considered as main sources of new drug development. However, in many cases, the herbs are still prescribed relying on ancient records and/or traditional practices without scientific evidences. More importantly, the medicinal efficacy of the herbs has to be evaluated in the perspective of MCMT (multi-compound multi-target) effects, but most efforts focus on identifying and analyzing a single compound experimentally. To overcome these hurdles, computational approaches which are based on the scientific evidences and are able to handle the MCMT effects are needed to predict the herb-disease associations.ResultsIn this study, we proposed a network-based in silico method to predict the herb-disease associations. To this end, we devised a new network-based measure, WACP (weighted average closest path length), which not only quantifies proximity between herb-related genes and disease-related genes but also considers compound compositions of each herb. As a result, we confirmed that our method successfully predicts the herb-disease associations in the human protein interactome (AUROC = 0.777). In addition, we observed that our method is superior than the other simple network-based proximity measures (e.g. average shortest and closest path length). Additionally, we analyzed the associations between Brassica oleracea var. italica and its known associated diseases more specifically as case studies. Finally, based on the prediction results of the WACP, we suggested novel herb-disease pairs which are expected to have potential relations and their literature evidences.ConclusionsThis method could be a promising solution to modernize the use of the natural herbs by providing the scientific evidences about the molecular associations between the herb-related genes targeted by multiple compounds and the disease-related genes in the human protein interactome.