Deciphering Molecular Mechanisms of Cutaneous Leishmaniasis, Pathogenesis and Drug Repurposing through Systems Biology

Abstract

Cutaneous leishmaniasis (CL) is a major health problem caused by an intracellular pathogen of the genus Leishmania. CL results in morphologically distinct skin injuries, ranging from nodules to plaques and ulcers, which persist as a recuperating incessant injury depending on the type of contaminating parasite. There is still no effective treatment to reduce the skin lesions in patients infected with CL. The aim of this study was to develop strategies to treat skin lesions in CL patients. We retrieved the transcriptomic data of skin lesions from patients with CL and normal skin from the gene Expression Omnibus (GEO) database. The protein-protein interaction network (PPIN) was constructed using the STRING database and Cytoscape v3.10.1 software. Critical genes were identified by topological network analysis and cluster detection. Finally, gene ontology and repurposing drugs for critical genes were determined. CD8A, IFNG, IL-6, PTPRC, CCR7, TLR2, GSTA5, CYBB, IL-12RB2, ITGB2, FCGR3A, CTLA4, and IFNG were identified as the critical genes in PPIN and subnetworks. Enrichment analysis revealed that T-cell receptor signaling, toll-like receptor signaling, cytokine-cytokine receptor interaction, graft-versus-host disease, leishmaniasis, chemokine signaling, primary immunodeficiency, and Th17 cell differentiation were the major pathways associated with critical genes. The drug repurposing results identified cyclosporine, rituximab, infliximab, blinatumomab, and methylprednisolone as candidates for treatment of CL. After validating our model with available experimental data, we found that critical molecules and drug candidates play a crucial role in the treatment of skin lesions caused by Leishmania in prospective studies.