Department of Bioinformatics, Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran.
School of Biological Science, Institute for Research in Fundamental Sciences (IPM), Tehran, Iran.
J Mol Graph Model. 2023 Jul;122:108495. doi: 10.1016/j.jmgm.2023.108495. Epub 2023 Apr 18.
Exploring allosteric inhibition and the discovery of new inhibitor binding sites are important studies in protein regulation mechanisms and drug discovery. Structural and network-based analyses of trajectories resulting from molecular dynamics (MD) simulations have been developed to discover protein dynamics, landscape, functions, and allosteric regions. Here, an experimentally suggested non-competitive inhibitor, xanthene-11v, was considered to explore its allosteric inhibition mechanism in α-glucosidase MAL12. Comparative structural and network analyses were applied to eight 250 ns independent MD simulations, four of which were performed in the free state and four of which were performed in ligand-bound forms. Projected two-dimensional free energy landscapes (FEL) were constructed from the probabilistic distribution of conformations along the first two principal components. The post-simulation analyses of the coordinates, side-chain torsion angles, non-covalent interaction networks, network communities, and their centralities were performed on α-glucosidase conformations and the intermediate sub-states. Important communities of residues have been found that connect the allosteric site to the active site. Some of these residues like Thr307, Arg312, TYR344, ILE345, Phe357, Asp406, Val407, Asp408, and Leu436 are the key messengers in the transition pathway between allosteric and active sites. Evaluating the probability distribution of distances between gate residues including Val407 in one community and Phe158, and Pro65 in another community depicted the closure of this gate due to the inhibitor binding. Six macro states of protein were deduced from the topology of FEL and analysis of conformational preference of free and ligand-bound systems to these macro states shows a combination of lock-and-key, conformational selection, and induced fit mechanisms are effective in ligand binding. All these results reveal structural states, allosteric mechanisms, and key players in the inhibition pathway of α-glucosidase by xanthene-11v.
探索别构抑制和发现新的抑制剂结合位点是蛋白质调节机制和药物发现的重要研究内容。基于结构和网络的轨迹分析方法已经被开发出来,用于发现蛋白质的动力学、景观、功能和别构区域。在这里,考虑了一种实验建议的非竞争性抑制剂,即香豆素-11v,以探索其在 MAL12 葡萄糖苷酶中的别构抑制机制。比较结构和网络分析应用于 8 个 250ns 独立的 MD 模拟,其中 4 个是在自由状态下进行的,4 个是在配体结合形式下进行的。从构象沿前两个主成分的概率分布构建二维自由能景观(FEL)。对α-葡萄糖苷酶构象和中间亚态的坐标、侧链扭转角、非共价相互作用网络、网络社区及其中心度进行了后模拟分析。已经发现了一些重要的残基社区,它们将别构位点与活性位点连接起来。这些残基中的一些,如 Thr307、Arg312、TYR344、ILE345、Phe357、Asp406、Val407、Asp408 和 Leu436,是别构和活性位点之间过渡途径中的关键信使。评估包括一个社区中的 Val407 在内的门控残基之间的距离概率分布,以及另一个社区中的 Phe158 和 Pro65,描绘了由于抑制剂结合导致该门关闭。从 FEL 的拓扑结构和自由和配体结合系统对这些宏观状态的构象偏好分析中推断出蛋白质的六个宏观状态,表明锁钥、构象选择和诱导契合机制的组合在配体结合中是有效的。所有这些结果揭示了香豆素-11v 抑制α-葡萄糖苷酶的结构状态、别构机制和抑制途径中的关键参与者。