Potential Mechanisms of Biejiajian Pill in the Treatment of Diabetic Atherosclerosis Based on Network Pharmacology, Molecular Docking, and Molecular Dynamics Simulation

BACKGROUND: Biejiajian pill (BJJP), a classical traditional Chinese formula, has been reported that it has an effective treatment for diabetic atherosclerosis in recent years, but its underlying mechanisms remain elusive. The study aimed to explore the potential mechanisms of BJJP on diabetic atherosclerosis by integrating network pharmacology, molecular docking, and molecular dynamics simulation. METHODS: The active components of BJJP were collected by TCMSP and TCMID, and then the potential targets were obtained from the SwissTargetPrediction database. The targets related to diabetic atherosclerosis were identified from the GeneCards and OMIM databases. The intersection of the potential targets regulated by active components of BJJP and the targets of diabetic atherosclerosis were common targets, which were visualized by the Venn diagram. The common targets were imported into the STRING database to construct a protein-protein interaction (PPI) network. The network of “Medicine-Compound-Target” was constructed with Cytoscape 3.7.1 software. GO functional enrichment analysis and KEGG pathway enrichment analysis were performed using the DAVID database and visualized through bioinformatics. The intersecting targets were input into Cytoscape 3.7.1 software, and the Network Analyzer tool was employed to screen out the key targets. Then molecular docking was used to verify the binding affinity between the active compounds and the key targets, and molecular dynamics simulation was used to investigate the stability of the binding models. RESULTS: A total of 81 active components, 186 targets of BJJP, and 4041 targets of diabetic atherosclerosis were obtained. Furthermore, 121 overlapping targets were identified. GO functional enrichment analysis revealed that these targets were correlated with the oxidation-reduction process, negative regulation of apoptotic process, inflammatory response, and other biological processes. The results of the KEGG pathway enrichment analysis showed that the common targets mainly participated in proteoglycans in cancer, PPAR signaling pathway, adherens junction, insulin resistance, HIF-1 signaling pathway, PI3K-Akt signaling pathway, etc. The results of molecular docking confirmed that the core active components in BJJP could bind well to the key targets. Results from molecular dynamics simulation showed that the binding energies of AKT1-Luteolin, MMP9-quercetin, and MMP9-luteolin complexes were −28.93 kJ·mol(−1), −37.12 kJ·mol(−1), and −62.91 kJ·mol(−1), respectively. CONCLUSION: The study revealed that BJJP is characterized as multicomponent, multitarget, and multipathway to treat diabetic atherosclerosis, which is helpful to provide ideas and a basis for pharmacological research and clinical application in the future.