Potential allosteric pockets identification of glucagon receptor based on molecular dynamics simulations

Glucagon receptor (GCGR) is an important target for the treatment of type 2 diabetes mellitus. Although several small molecules with antagonistic activity have been discovered, so far, only one small molecule binding site has been resolved. To discover more novel allosteric pockets and allosteric molecules, we started with the unique full-length inactive conformation of GCGR and applied all-atom molecular dynamics (MD) simulations to obtain extensive dynamic conformations of the GCGR/glucagon complex. For the first time, MDpocket, FTMove and FTMap were used to detect allosteric pockets in simulation trajectories, selecting 4 stable pockets with a total of 14 structures as templates for virtual screening. From the results of virtual screening, 14 compounds were ultimately selected after a series of filtering steps. The cAMP accumulation assay indicated that compound gs6 has antagonistic activity, and MD simulations further revealed the allosteric mechanism of gs6. We are the first to identify new allosteric pockets and allosteric molecules in simulation trajectories of the GCGR/glucagon complex, providing a reference for research on other G-protein-coupled receptors (GPCR). However, there is still considerable room for improvement, such as using more simulation methods to obtain a richer set of dynamic conformations.