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.

The Inferential Binding Sites of GCGR for Small Molecules Using Protein Dynamic Conformations and Crystal Structures.

Glucagon receptor (GCGR) is a class B1 G-protein-coupled receptor that plays a crucial role in maintaining human blood glucose homeostasis and is a significant target for the treatment of type 2 diabetes mellitus (T2DM). Currently, six small molecules (Bay 27-9955, MK-0893, MK-3577, LY2409021, PF-06291874, and LGD-6972) have been tested or are undergoing clinical trials, but only the binding site of MK-0893 has been resolved. To predict binding sites for other small molecules, we utilized both the crystal structure of the GCGR and MK-0893 complex and dynamic conformations. We docked five small molecules and selected the best conformation based on binding mode, docking score, and binding free energy. We performed MD simulations to verify the binding mode of the selected small molecules. Moreover, when selecting conformations, results of competitive binding were referred to. MD simulation indicated that Bay 27-9955 exhibits moderate binding stability in Pocket 3. MK-3577, LY2409021, and PF-06291874 exhibited highly stable binding to Pocket 2, consistent with experimental results. However, LY2409021 may also bind to Pocket 5. Additionally, LGD-6972 exhibited relatively stable binding in Pocket 5. We also conducted structural modifications of LGD-6972 based on the results of MD simulations and predicted its analogues' bioavailability, providing a reference for the study of GCGR small molecules.