Comprehensive Network Analysis Reveals the Targets and Potential Multitarget Drugs of Type 2 Diabetes Mellitus

Abstract

Type 2 diabetes mellitus (T2DM) is a metabolic disease with increasing prevalence and mortality year by year. The purpose of this study was to explore new therapeutic targets and candidate drugs for multitargets by single-cell RNA expression profile analysis, network pharmacology, and molecular docking. Single-cell RNA expression profiling of islet β cell samples between T2DM patients and nondiabetic controls was conducted to identify important subpopulations and the marker genes. The potential therapeutic targets of T2DM were identified by the overlap analysis of insulin-related genes and diabetes-related genes, the construction of protein-protein interaction network, and the molecular complex detection (MCODE) algorithm. The network distance method was employed to determine the potential drugs of the target. Molecular docking and molecular dynamic simulations were carried out using AutoDock Vina and Gromacs2019, respectively. Eleven cell clusters were identified by single-cell RNA sequencing (scRNA-seq) data, and three of them (C2, C8, and C10) showed significant differences between T2DM samples and normal samples. Eight genes from differential cell clusters were found from differential cell clusters to be associated with insulin activity and T2DM. The MCODE algorithm built six key subnetworks, with five of them correlating with inflammatory pathways and immune cell infiltration. Importantly, CCR5 was a gene within the key subnetworks and was differentially expressed between normal samples and T2DM samples, with the highest area under the ROC curve (AUC) of 82.5% for the diagnosis model. A total of 49 CCR5-related genes were screened, and DB05494 was identified as the most potential drug with the shortest distance to CCR5-related genes. Molecular docking illustrated that DB05494 stably bound with CCR5 (-8.0 kcal/mol) through multiple hydrogen bonds (LYS26, TYR37, TYR89, CYS178, and GLN280) and hydrophobic bonds (TRP86, PHE112, ILE198, TRP248, and TYR251). This study identified CCR5 as a potential therapeutic target and screened DB05494 as a potential drug for T2DM treatment.