Abstract
The high conservation of the three subtypes of glycogen phosphorylase (GP) presents significant challenges for specific inhibitor studies targeting GP. Our prior screening revealed that compound 1 exhibited unequal inhibitory activity against the three GP subtypes, with a noticeable effect against brain GP (PYGB). The commercially available ingliforib demonstrated potent inhibitory activity specifically against liver GP (PYGL). To guide the further design and screening of high-specificity inhibitors, the possible reasons for the differential inhibitory activity of two compounds against different GP subtypes were analyzed, with ingliforib as a reference, through molecular docking and molecular dynamics simulations. Initially, the study predicted the binding modes of ligands with the three GP receptor subtypes using molecular docking. Subsequently, this was validated by molecular dynamics experiments, and possible amino acid residues that had important interactions were explored. The strong correlation between the calculated interaction free energies and experimental inhibitory activity implied the reasonable binding conformations of the compounds. These findings offer insight into the different inhibitory activity of compound 1 and ingliforib against all three GP subtypes and provide guidance for the design of specific target molecules that regulate subtype selectivity.