Study on the Dynamics and Key Sites of the μ-Opioid Receptor Based on the Anisotropic Network Model
Opioids mainly achieve analgesic effects by selectively acting on μ-opioid receptors (MOR) and activating inhibitory G proteins. Since its use is often accompanied by serious adverse reactions such as addiction and respiratory depression, the development of drugs targeting opioid receptors has become a major challenge. To investigate the structural dynamics and functional key sites of MOR, we first analyzed the residue fluctuations under low-frequency and high-frequency motion modes and the cross-correlation fluctuations between residues based on the Anisotropic network model (ANM). Then, a complex network of amino acids is constructed to identify key residues associated with allosteric signaling. The results showed that the lowest mode effectively identified functional regions within the structure, the fastest mode identified hotspot residues important for stabilizing the MOR inactive structure, and motion correlation analysis revealed the transduction mechanism of allosteric signals in the transmembrane helical region. Additionally, the complex network identified key residues that play important roles in ligand binding and signaling. This work helps to understand the dynamics characteristics and key sites of μ-opioid receptors, and provides useful reference information for drug design.
μ-Opioid Receptor MOR
阿片受体(opioid receptors, ORs)负责调控镇痛、奖赏
2022年,Robertson等人通过冷冻电镜技术解析了与拮抗剂Alvimopan结合的非活性MOR (残基65-352)结构
目前,许多计算模拟方法被应用于MOR变构机制的研究。Yuan等使用长时间尺度的分子动力学(molecular dynamics, MD)模拟探究了水分子和钠离子在MOR跨膜信号传导过程中的作用
本工作基于ANM研究MOR的结构动力学与功能间的关系,并结合复杂网络模型识别蛋白质中与变构通信相关的关键位点。
从蛋白质结构数据库 (Protein data bank, PDB)
(1)
其中,γ是弹簧的力常数, ,N为残基数,上标T表示向量的转置,H表示3N × 3N的海森矩阵,其元素为:
(2)
当 时,子矩阵Hij为:
(3)
当 时,子矩阵Hii为:
(4)
海森矩阵分解得到的特征向量表示蛋白质的运动模式,残基的均方涨落(mean square fluctuations, MSF)和残基间的交叉相关涨落(correlated fluctuation, CoF)为:
(5)
(6)
其中,kB是Boltzmann常数,T为绝对温度。归一化的运动交叉相关性为:
(7)
Cij的取值范围在[−1, 1]之间,其中正值表示残基间的运动方向相同,而负值则表示运动方向相反。Cij的绝对值越大,意味着残基间的运动相关性越强。
蛋白质的结构和功能依赖于氨基酸残基之间复杂的相互作用网络。在加权的复杂网络模型
(8)
特征路径长度(Characteristic Path Length, CPL)是网络中所有节点对之间最短路径的平均长度:
(9)
其中N和Np分别是节点数和节点对数,di j表示从节点i到j累计边权重最小的路径长度。节点k对网络内信息通信的贡献是通过从网络中移除节点k后CPL的变化(ΔCPLk)来衡量
(10)
其中,ΔCPLk是移除节点k后特征路径长度的变化, 是依次移除各节点后特征路径长度变化的平均值,σ是标准差。
温度因子(B-factor)可以体现分子的局部柔性。对MOR构建模型时,通过最大化理论B-factor和实验之间的皮尔森相关系数(Pearson correlation coefficient, PCC)
慢运动模式通常代表与蛋白功能相关的大幅度集合运动,
在快运动模式下,表现出较大运动幅度的残基往往被视为动力学上的热点残基,这些残基在稳定蛋白质空间结构中起着至关重要的作用
为研究MOR中残基之间的功能性运动耦合,基于公式(7)计算了残基涨落之间的动态互相关矩阵 (Dynamic cross-correlation matrix, DCCM),
接收和传递变构信号的重要残基通常是蛋白质相互作用网络的中心,位于大多数残基之间的最短路径上。通过依次移除氨基酸网络中的节点及其相连通路,计算特征路径长度变化的Z-score分数,来量化
残基对网络内通信的影响,结果如
本工作研究了μ阿片受体 (MOR) 的动力学特性,以及对蛋白变构信号转导起重要作用的残基。在最佳参数下构建的受体的ANM模型,计算得到的B-factor与实验值的PCC达到0.79,较好地再现了各结构域的柔性,证明ANM可用于研究MOR的结构动力学与功能的关系。在慢运动模式下,ECL2和ICL3有大规模的集合运动,猜测ECL2的移动有助于配体结合口袋的打开,ICL3的构象变化负责调节受体活性,涨落比较小的残基与MOR的二聚化、激活密切相关。快运动模式下识别出了可以稳定受体内外两侧非活性结构的关键残基。运动相关性分析表明变构信号主要是通过TM2、TM3、TM5、TM6和TM7的耦合构象变化传导的。最后,加权的氨基酸网络模型预测了MOR的功能性关键残基,大部分结果与实验数据吻合,个别残基还未被证实,可为未来的实验研究提供参考。本工作有助于理解MOR的变构通讯机制,对相关药物设计具有指导意义。