Journal of Biomolecular Structure and Dynamics, cilt.42, sa.2, ss.819-833, 2024 (SCI-Expanded)
A complete investigation to understand the pathways that could be affected by glycyrrhizin (licorice), as anti-breast cancer (BC) agent, has not been performed to date. This study aims to investigate the pathways involved in the anti-cancer activity of glycyrrhizin against BC. For this purpose, the target genes of glycyrrhizin were obtained from the ChEMBL database. The BC-associated genes for three types of BC (breast carcinoma, malignant neoplasm of breast, and triple-negative breast neoplasms) were retrieved from DisGeNET. The target genes of glycyrrhizin and the BC-associated genes were compared, and the genes with disease specificity index (DSI) > 0.6 were selected for further evaluation using in silico methods. The protein-protein interaction (PPI) network was constructed, and the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were analyzed. The potential complexes were further evaluated using molecular dynamics (MD) simulation. The results revealed that among 80 common genes, ten genes had DSI greater than 0.6, which included POLK, TACR2, MC3R, TBXAS1, HH1R, SLCO4A1, NPY2R, ADRA2C, ADRA1A, and SLCO2B1. The binding affinity of glycyrrhizin to the cognate proteins and binding characteristics were assessed using molecular docking and binding free energy calculations (MM/GBSA). POLK, TBXAS1, and ADRA1A showed the highest binding affinity with −8.9, −9.3, and −9.6 kcal/mol, respectively. The final targets had an association with BC at several stages of tumor growth. By affecting these targets, glycyrrhizin could influence and control BC efficiently. MD simulation suggested the pathways triggered by the complex glycyrrhizin-ADRA1A were more likely to happen. Communicated by Ramaswamy H. Sarma.