A network-based model for drug repurposing in Rheumatoid Arthritis

Background The identification of drug repositioning targets through computational methods has the potential to provide a fast, inexpensive alternative to traditional drug discovery process. Diseases where complicated pathophysiology evolves over time are excellent targets for such methods. One such disease is Rheumatoid Arthritis (RA), a chronic inflammatory autoimmune disease, where the drug survival at 5-years is less than 50% for patients treated with disease-modifying anti-rheumatic drugs (DMARDs). Methods and Findings We have developed a network-based approach for drug repositioning that takes into account the human interactome network, proximity measures between drug targets and disease-associated genes, potential side-effects, genome-wide gene expression and disease modules that emerge through pertinent analysis. We found that all DMARDs, except for hydroxychloroquine (HCQ), were found to be significantly proximal to RA-related genes. Application of the method on anti-diabetic agents, statins and H2 receptor blockers identified anti-diabetic agents – gliclazide, sitagliptin and metformin – that have similar network signatures with the DMARDs. Subsequent in-vitro experiments on mouse fibroblast NIH-3T3 cells validated the findings and the down regulation of six key RA-related inflammatory genes. Our analysis further argues that the combination of HCQ and/or sulfasalazine with methotrexate (MTX) is predicted to have an additive synergistic effect in treatment based on network complementarity. Similarly, leflunomide and tofacitinib were found to be suitable alternatives upon chemoresistance to MTX-based double/triple therapy, given the complementary network signatures and overlapping critical target hubs. Conclusions Our results corroborate that computational methods that are based on network proximity, among other contextual information can help narrow down the drug candidates for drug repositioning, as well as support decisions for combinatorial drug treatment that is tailored to patient’s needs. Author summary The network-based proximity between drug targets and disease genes can provide novel insights on the repercussion, interplay, and reposition of drugs in the context of disease. Disease-modifying anti-rheumatic drugs (DMARDs) located significantly close to rheumatoid arthritis (RA)-associated genes and RA-relevant pathways. Three anti-diabetic agents were identified to have an anti-inflammatory effect like DMARDs. We built RA disease module encompassing the emerging small-molecule targets and functional neighbors, which better explained the RA pathophysiology. By proximity and network robustness, tofacitinib and tocilizumab were the most potent for RA disease module as a single agent, and this is consistent with clinical observation. Side effects of clinical importance are predictable by measuring network-based proximity between drug targets and side effect protein. Network-based drug-disease proximity offers a novel and clinically actionable information about drugs and opens the new possibility to drug combination and reposition.

[1]  M. Scherer,et al.  Attenuated Wnt/β-catenin signalling mediates methotrexate chemotherapy-induced bone loss and marrow adiposity in rats. , 2012, Bone.

[2]  P. Casey,et al.  Targeting Ras signaling through inhibition of carboxyl methylation: An unexpected property of methotrexate , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[3]  C. Werning [Rheumatoid arthritis]. , 1983, Medizinische Monatsschrift fur Pharmazeuten.

[4]  A. Barabasi,et al.  Uncovering disease-disease relationships through the incomplete interactome , 2015, Science.

[5]  R. Altman,et al.  Pharmacogenomics Knowledge for Personalized Medicine , 2012, Clinical pharmacology and therapeutics.

[6]  J. Mesirov,et al.  The Molecular Signatures Database (MSigDB) hallmark gene set collection. , 2015, Cell systems.

[7]  Raveendhara R. Bannuru,et al.  American College of Rheumatology Guideline for the Treatment of Rheumatoid Arthritis , 2015 .

[8]  I. Tagkopoulos,et al.  Compendium of Synovial Signatures Identifies Pathologic Characteristics for Predicting Treatment Response in Rheumatoid Arthritis , 2018, bioRxiv.

[9]  O. Serri,et al.  Benefits of gliclazide in the atherosclerotic process: decrease in monocyte adhesion to endothelial cells. , 2003, Metabolism: clinical and experimental.

[10]  Albert-László Barabási,et al.  A DIseAse MOdule Detection (DIAMOnD) Algorithm Derived from a Systematic Analysis of Connectivity Patterns of Disease Proteins in the Human Interactome , 2015, PLoS Comput. Biol..

[11]  A. Barabasi,et al.  A disease module in the interactome explains disease heterogeneity, drug response and captures novel pathways and genes in asthma. , 2015, Human molecular genetics.

[12]  Núria Queralt-Rosinach,et al.  DisGeNET: a discovery platform for the dynamical exploration of human diseases and their genes , 2015, Database J. Biol. Databases Curation.

[13]  J. Isaacs,et al.  Mechanism of action of methotrexate in rheumatoid arthritis, and the search for biomarkers , 2016, Nature Reviews Rheumatology.

[14]  Lon Phan,et al.  Phenotype–Genotype Integrator (PheGenI): synthesizing genome-wide association study (GWAS) data with existing genomic resources , 2013, European Journal of Human Genetics.

[15]  S. Momohara,et al.  Molecular aspects of rheumatoid arthritis: chemokines in the joints of patients , 2008, The FEBS journal.

[16]  Monica L. Mo,et al.  Global reconstruction of the human metabolic network based on genomic and bibliomic data , 2007, Proceedings of the National Academy of Sciences.

[17]  T. Witte Methotrexate as combination partner of TNF inhibitors and tocilizumab. What is reasonable from an immunological viewpoint? , 2015, Clinical Rheumatology.

[18]  P. Bork,et al.  Systematic identification of proteins that elicit drug side effects , 2013, Molecular systems biology.

[19]  E. Schadt Molecular networks as sensors and drivers of common human diseases , 2009, Nature.

[20]  G. Mills,et al.  Functional variomics and network perturbation: connecting genotype to phenotype in cancer , 2017, Nature Reviews Genetics.

[21]  Teresa M. Przytycka,et al.  Identifying Causal Genes and Dysregulated Pathways in Complex Diseases , 2011, PLoS Comput. Biol..

[22]  Rafael C. Jimenez,et al.  The IntAct molecular interaction database in 2012 , 2011, Nucleic Acids Res..

[23]  Henning Hermjakob,et al.  The Reactome pathway Knowledgebase , 2015, Nucleic acids research.

[24]  Lin Yang,et al.  Emerging roles of dipeptidyl peptidase 4 inhibitors: anti-inflammatory and immunomodulatory effect and its application in diabetes mellitus. , 2014, Canadian journal of diabetes.

[25]  Kakajan Komurov,et al.  An integrated network platform for contextual prioritization of drugs and pathways. , 2015, Molecular bioSystems.

[26]  Sandhya Rani,et al.  Human Protein Reference Database—2009 update , 2008, Nucleic Acids Res..

[27]  Albert-László Barabási,et al.  Error and attack tolerance of complex networks , 2000, Nature.

[28]  Ben S. Wittner,et al.  Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1 , 2009, Nature.

[29]  Bai Zhang,et al.  Network Biology in Medicine and Beyond , 2014, Circulation. Cardiovascular genetics.

[30]  Ruth Nussinov,et al.  Structure and dynamics of molecular networks: A novel paradigm of drug discovery. A comprehensive review , 2012, Pharmacology & therapeutics.

[31]  Damian Szklarczyk,et al.  STITCH 5: augmenting protein–chemical interaction networks with tissue and affinity data , 2015, Nucleic Acids Res..

[32]  Jacob Köhler,et al.  Addressing the problems with life-science databases for traditional uses and systems biology , 2006, Nature Reviews Genetics.

[33]  Deng Pan,et al.  DGIdb 2.0: mining clinically relevant drug–gene interactions , 2015, Nucleic Acids Res..

[34]  Min-jung Park,et al.  Metformin Attenuates Experimental Autoimmune Arthritis through Reciprocal Regulation of Th17/Treg Balance and Osteoclastogenesis , 2014, Mediators of inflammation.

[35]  J. Eckel,et al.  Soluble DPP4 induces inflammation and proliferation of human smooth muscle cells via protease-activated receptor 2. , 2014, Biochimica et biophysica acta.

[36]  A. Vinayagam,et al.  A Directed Protein Interaction Network for Investigating Intracellular Signal Transduction , 2011, Science Signaling.

[37]  Livia Perfetto,et al.  MINT, the molecular interaction database: 2012 update , 2011, Nucleic Acids Res..

[38]  P. V. van Riel,et al.  Ups and downs in the treatment strategies of rheumatoid arthritis. , 2011, Rheumatology.

[39]  Gary D Bader,et al.  Enrichment Map: A Network-Based Method for Gene-Set Enrichment Visualization and Interpretation , 2010, PloS one.

[40]  J. Mesirov,et al.  The Molecular Signatures Database Hallmark Gene Set Collection , 2015 .

[41]  Bin Zhang,et al.  PhosphoSitePlus: a comprehensive resource for investigating the structure and function of experimentally determined post-translational modifications in man and mouse , 2011, Nucleic Acids Res..

[42]  Kara Dolinski,et al.  The BioGRID interaction database: 2017 update , 2016, Nucleic Acids Res..

[43]  Alexander E. Kel,et al.  TRANSFAC®: transcriptional regulation, from patterns to profiles , 2003, Nucleic Acids Res..

[44]  Damian Szklarczyk,et al.  The STRING database in 2017: quality-controlled protein–protein association networks, made broadly accessible , 2016, Nucleic Acids Res..

[45]  R. Westhovens,et al.  The efficacy and safety of abatacept in rheumatoid arthritis , 2010, Therapeutic advances in musculoskeletal disease.

[46]  I. Komuro,et al.  A DPP-4 Inhibitor Suppresses Fibrosis and Inflammation on Experimental Autoimmune Myocarditis in Mice , 2015, PloS one.

[47]  J. Smolen,et al.  Rheumatoid arthritis therapy reappraisal: strategies, opportunities and challenges , 2015, Nature Reviews Rheumatology.

[48]  Minoru Kanehisa,et al.  KEGG: new perspectives on genomes, pathways, diseases and drugs , 2016, Nucleic Acids Res..

[49]  F. Houssiau,et al.  Global Molecular Effects of Tocilizumab Therapy in Rheumatoid Arthritis Synovium , 2014, Arthritis & rheumatology.

[50]  John Wong,et al.  EULAR recommendations for the management of rheumatoid arthritis with synthetic and biological disease-modifying antirheumatic drugs , 2010, Annals of the rheumatic diseases.

[51]  David S. Wishart,et al.  DrugBank: a knowledgebase for drugs, drug actions and drug targets , 2007, Nucleic Acids Res..

[52]  Brad T. Sherman,et al.  Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources , 2008, Nature Protocols.

[53]  Pei-Ming Yang,et al.  Methotrexate induces apoptosis through p53/p21-dependent pathway and increases E-cadherin expression through downregulation of HDAC/EZH2. , 2011, Biochemical pharmacology.

[54]  angesichts der Corona-Pandemie,et al.  UPDATE , 1973, The Lancet.

[55]  A. Barabasi,et al.  Network-based in silico drug efficacy screening , 2016, Nature Communications.

[56]  Hans-Werner Mewes,et al.  CORUM: the comprehensive resource of mammalian protein complexes , 2007, Nucleic Acids Res..

[57]  Tsutomu Takeuchi,et al.  EULAR recommendations for the management of rheumatoid arthritis with synthetic and biological disease-modifying antirheumatic drugs: 2016 update , 2010, Annals of the rheumatic diseases.

[58]  Susumu Goto,et al.  KEGG for integration and interpretation of large-scale molecular data sets , 2011, Nucleic Acids Res..

[59]  Amy R. Cameron,et al.  Anti-Inflammatory Effects of Metformin Irrespective of Diabetes Status , 2016, Circulation research.

[60]  P. Bork,et al.  A side effect resource to capture phenotypic effects of drugs , 2010, Molecular systems biology.

[61]  Massimo Marchiori,et al.  Error and attacktolerance of complex network s , 2004 .