《Protein Science》：Structural insights into signal transduction in photoactivated adenylyl cyclase through molecular dynamics simulations
编辑推荐：
来自Oscillatoria acuminata的OaPAC是一种同源二聚体蛋白，由蓝光利用黄素结构域（BLUF domain）和腺苷酸环化酶结构域（AC domain）组成，两者通过长卷曲螺旋α3-螺旋连接。在蓝光吸收后，涉及黄素发色团、Y6和Q48的氢键网
来自Oscillatoria acuminata的OaPAC是一种同源二聚体蛋白，由蓝光利用黄素结构域（BLUF domain）和腺苷酸环化酶结构域（AC domain）组成，两者通过长卷曲螺旋α3-螺旋连接。在蓝光吸收后，涉及黄素发色团、Y6和Q48的氢键网络重排触发信号通过α3-螺旋向AC结构域传播。尽管近期对光态下W90–M92开关构象的晶体学结构提供了重要的结构信息，但信号转导的原子层面机制仍不清楚。在本研究中，研究人员对四种OaPAC模型进行了分子动力学模拟：野生型（WT）暗态、W90–M92未开关的WT光态、W90–M92已开关的WT光态以及W90–M92未开关的L111A/L115A光态。结果表明，已开关的WT光态模型的结构波动约为未开关的WT光态的两倍。此外，研究人员在已开关的WT光态中识别出两种不同构象状态之间的显著转变——“开放”和“闭合”——其特征由α3-螺旋与“手柄”区域之间的角度定义。值得注意的是，闭合构象在已开关的WT光态中占有率最高。相比之下，无功能的L111A/L115A突变体则被锁定在刚性的开放构象中，反映了其活性的丧失，并将结构柔性与酶活性联系起来。这些发现表明，W90–M92构象开关是OaPAC信号转导的一个关键的多步触发器，诱导瞬时的非对称（开放–闭合）构象和可能促进酶激活的屈曲运动。该研究将为先进光遗传学执行器的理性设计提供关键的原子层面见解。
论文解读文章
研究背景与目的
环磷酸腺苷（cyclic adenosine monophosphate, cAMP）是最通用的第二信使之一，在将细胞外信号转化为细胞内响应中发挥关键作用。腺苷酸环化酶（adenylyl cyclase, AC）是催化三磷酸腺苷（adenosine triphosphate, ATP）生成cAMP的关键酶。光激活腺苷酸环化酶（photoactivated adenylyl cyclase, PAC）能够在光照刺激下产生cAMP，无需G蛋白参与，因此成为精准时空控制细胞功能的强大工具，广泛应用于神经科学和心脏病学等领域。
来自Oscillatoria acuminata的OaPAC是一个同源二聚体多结构域蛋白，每个单体由蓝光利用黄素结构域（BLUF domain，残基1–102）和腺苷酸环化酶结构域（AC domain，残基127–366）组成，两者通过卷曲螺旋α3-螺旋（残基103–125）和称为“手柄”（handle）的α4-螺旋（残基127–143）连接。BLUF光感受器区域与AC催化区域相距约45–48 ?，排除了直接相互作用。因此，信号传递被认为是通过α3-螺旋进行的。
尽管近年来在OaPAC的结构解析方面取得了重要进展，包括暗态和光态的晶体结构，以及通过飞秒X射线晶体学（serial femtosecond x-ray crystallography, SFX）捕获的W90–M92开关构象，但光吸收如何在BLUF结构域触发AC结构域的催化激活这一分子机制仍不完全清楚。特别是，BLUF结构域的结构变化如何通过45–48 ?的距离动态传递到AC活性位点，这一原子层面的机制尚属空白。本研究旨在通过分子动力学（molecular dynamics, MD）模拟填补这一认知缺口。
研究方法概述
研究人员采用分子动力学模拟作为核心技术手段，构建了四种OaPAC模拟模型：野生型（wild-type, WT）暗态、W90–M92未开关的WT光态、W90–M92已开关的WT光态以及W90–M92未开关的L111A/L115A光态。初始结构分别基于暗态晶体结构（PDB ID: 4YUS）和W90–M92开关光态晶体结构（PDB ID: 8QFG）。发色团区域的原子电荷和坐标取自前期混合量子力学/分子力学（quantum mechanics/molecular mechanics, QM/MM）自由能几何优化结果。模拟采用AMBER ff14SB力场和TIP3P水模型，每个模型进行五个独立3 μs的生产运行。轨迹分析使用CPPTRAJ模块和MDTraj Python库进行。
研究结果
2.1 OaPAC的整体结构动力学
通过计算主链均方根涨落（root-mean-square fluctuation, RMSF），研究人员发现已开关WT光态的RMSF值约为未开关WT光态的两倍，表明W90–M92开关构象对于结构传播至关重要。二级结构分析显示，各模型在全长蛋白中的二级结构谱几乎相同，但已开关WT光态的C末端无序区域增加。这表明OaPAC激活不涉及大规模二级结构转变，而是依赖于BLUF结构域特定区域的局部结构变化。
2.2 α3-螺旋和手柄区域的构象转变
在已开关WT光态中，α3-螺旋（残基103–125）与手柄区域（残基127–143）之间的角度表现出显著的构象转变，分为开放构象（角度>40°）和闭合构象（角度<40°）。已开关WT光态的角度分布广泛分散，而未开关WT光态集中在两个晶体结构值附近。定量数据显示，已开关WT光态中不对称开放–闭合构象的占有率最高（63.7%），而L111A/L115A光态则完全锁定在开放构象（90.0%），表明W90–M92开关诱导了OaPAC信号转导所需的大规模结构重排。
2.3 氢键重排驱动β-片层移位和结构域间波动
氢键分析确认了β5-片层向C末端的两位点移位在模拟中保持稳定。在已开关WT光态中，β4-片层也相对于其他模型向C末端移位。跨单体氢键分析揭示了K77(A)–E117(B)盐桥在已开关WT光态中被显著破坏，而在WT暗态和未开关光态中保持稳定。这种“系绳”作用的丧失可能触发了BLUF核心与α3-螺旋之间的结构波动增加，从而促进信号转导。
2.4 单体内和单体间相互作用的接触图谱分析
接触图谱分析显示，在已开关与未开关WT光态的对比中，β4-片层与α4-螺旋之间的相互作用发生显著变化，反映手柄区域的构象转变。同时，α3-螺旋与β2/β3-片层的接触高度局部化。在未开关WT光态与暗态的对比中，研究人员识别出α3-螺旋根部与β5-片层之间的特异性残基对，这些可作为W90–M92开关的结构前体。
2.5 信号转导的残基层面动力学
β5-α3连接子（残基97–102）在已开关WT光态中展现出更高的螺旋形成倾向，这可能作为容纳β5-片层下游两位点移位的机械触发器。M92侧链动力学在暗态和未开关光态之间存在明显差异，光态下M92向I22靠近，这可能是W90–M92开关的前体行为。W90与P107之间的距离分析显示，未开关WT光态中6–7 ?范围的特定群体消失，进一步支持了W90动力学的调制。
2.6 主成分分析揭示的优势动力学与功能构象模式
主成分分析（principal component analysis, PCA）识别出三种主要模式：模式1为屈曲运动（BLUF与AC结构域间距减小），模式2为扭转运动，模式3为AC结构域的打开运动。在已开关WT光态中，模式1的贡献显著增加至54.7%，成为优势模式。这种屈曲运动被认为是驱动AC结构域构象变化的关键因素。相比之下，L111A/L115A变异体的所有模式缩放特征值均极小，表明该突变无法诱导达到瞬时状态所需的功能性运动。
讨论与结论
本研究通过MD模拟强有力地支持了W90–M92开关是驱动OaPAC信号转导的关键事件的假说。RMSF分析表明，该开关诱导了最显著的结构波动，而仅靠黄素周围氢键网络的重排不足以驱动大规模波动。L111A/L115A突变抑制了BLUF结构域中的这些波动，突显了结构柔性与酶活性之间的强相关性。然而，结构柔性本身不应被视为催化能力的直接证据，更可能是适当的柔性有助于访问与信号转导相关的特定集体运动和构象转变。
基于研究发现，研究人员提出了OaPAC的精细化信号转导机制：首先，蓝光照射诱导Q48侧链的酰胺-亚氨酸互变异构化和旋转，重排涉及黄素发色团、Y6和Q48的氢键网络；其次，这一重排增加了M92侧链的构象自由度，使其能够破坏BLUF结构域的疏水核心；第三，结构波动的增加使W90发生构象变化，通过W90–M92开关实现β5-片层向新构象的转变；最后，α3-螺旋与手柄区域之间的角度从开放转变为闭合构象，伴随屈曲运动使BLUF结构域与AC结构域靠近，最终驱动AC结构域发生催化反应所需的构象变化。
本研究的结构见解为理解优化PAC作为光遗传学工具所需机制奠定了基础，为先进光遗传学执行器的理性设计提供了关键的原子层面信息。
