与DNA链相互作用的蛋白质分子

《Sensing and Bio-Sensing Research》：Protein molecules engaging DNA strands

　　

本研究探讨了与DNA链耦合的蛋白质分子中的相干电子输运现象，该系统在紧束缚（TB）框架内建模，并采用朗道尔-布蒂克尔（LB）形式主义进行分析。研究人员分析了三种具有不同氢键拓扑结构的蛋白质构象，以考察其对透射谱和电流-电压（I-V）特性的影响，而DNA链则被建模为鱼骨型（FBM）电极。所得的I-V特性呈现出显著的非线性行为，这源于杂化结的能级依赖透射谱。电极温度的影响通过电子通路的热展宽效应纳入考量，揭示了随着温度升高，输运行为从阶梯状逐渐向准线性转变。此外，蛋白质结构的构象变化被证明可改变可用传输通道的数量和对齐方式，从而导致输运响应出现可测量的差异。这些结果为蛋白质构象和热效应如何共同影响蛋白质-DNA杂化系统中的相干电子输运提供了模型层面的认识。尽管本研究纯粹是理论性的，未包含明确的环境效应或直接实验对比，但其为分子结构和蛋白质-DNA耦合在调控相干输运中的作用提供了定性物理认识，从而为生物启发纳米电子结的未来研究奠定了理论框架。

**研究背景与问题提出**

脱氧核糖核酸（DNA）作为由腺嘌呤（A）与胸腺嘧啶（T）、鸟嘌呤（G）与胞嘧啶（C）通过氢键配对形成的双螺旋生物大分子，因其π电子堆积结构和电荷传输能力，在纳米电子学与生物传感技术领域展现出作为纳米导线的应用潜力。与此同时，蛋白质因其在纳米生物电子学中的重要作用以及参与生物体内电子过程的特性，被视为生物电子应用的关键候选材料。蛋白质与DNA的相互作用是DNA复制、修复、重组、基因表达调控及染色体组织等众多细胞过程的基础，这些相互作用使细胞能够解读和利用储存在DNA中的遗传信息以合成蛋白质并执行基本生物学功能。然而，尽管已有研究分别探索了DNA或孤立蛋白质的紧束缚模型，但将多种蛋白质构象与DNA电极耦合并在有限温度条件下进行分析的工作仍较为缺乏，特别是氢键拓扑结构对传输通道的影响及其与热效应的相互作用尚未得到系统研究。因此，开展这项研究旨在填补这一空白，揭示蛋白质构象和氢键拓扑结构如何与DNA介导的耦合及电极温度共同决定相干电子输运特性。

**研究内容、结论与意义**

研究人员采用粗粒化的三维晶格模型表示蛋白质，以鱼骨模型（FBM）描述DNA电极，运用紧束缚框架结合朗道尔-布蒂克尔形式主义，在非平衡格林函数框架下计算了完整的I-V响应。研究构建了三种不同的蛋白质构象（P1、P2、P3），每种构象由36个氨基酸残基组成，但氢键排列方式不同：P1为基本结构，氨基酸仅通过共价肽键连接；P2和P3则同时包含肽键和非共价的氢键相互作用。将这些蛋白质构象置于两条作为电极的DNA链之间，研究人员分析了电极温度变化（T=0.001、0.004、0.008，对应物理温度3.48 K、13.92 K、27.84 K）对输运特性的影响，并进行了蛋白质-电极耦合强度的敏感性分析。研究得出的主要结论包括：蛋白质-DNA杂化系统中的相干电子输运对分子构象和温度敏感，但显示出对耦合参数中等变化的稳定性；氢键拓扑结构作为结构控制参数，能够有选择地开启或抑制电子通路；I-V特性可作为区分不同蛋白质构象（包括正常和突变结构）的敏感指标。该论文发表于《Sensing and Bio-Sensing Research》，其意义在于为蛋白质-DNA杂化结作为可调谐且可靠的生物电子和纳米尺度器件组分提供了理论基础，建立了用于未来生物启发纳米电子结研究的理论框架。

**主要关键技术方法**

研究采用紧束缚哈密顿量描述系统，其中蛋白质器件由包含在位能和最近邻电子耦合的项组成，温度效应通过相邻格点间扭曲角的高斯分布涨落引入；电极采用鱼骨模型的二次量子化哈密顿量描述，求解薛定谔方程获得波矢；运用朗道尔-布蒂克尔理论结合费希尔-李（Fisher–Lee）关系计算透射概率，通过电极自能矩阵和蛋白质器件的推迟格林函数求解电流；态密度由格林函数迹的虚部提取。样本方面，三种蛋白质构象来源于文献报道的基于立方晶格自避行走设计的36氨基酸残基链，其中P2和P3包含13个氢键以再现疾病相关突变。

**研究结果**

**温度对输运特性的影响**：研究人员通过图2展示了三种蛋白质构象在三个不同电极温度下的I-V特性、透射概率（TP）谱和态密度（DOS）。结果显示，I-V曲线呈现明显的非线性行为，表现为电流的急剧增加后接近平坦的平台。随着电极温度升高，这些阶梯状特征逐渐平滑，趋于更线性的外观。对应地，TP谱中的峰展宽并移位，DOS谱则显示出构象依赖的移位和共振强度变化，与TP峰及非线性电流响应直接相关。这一现象归因于电极中电子能量分布的热展宽：低温时费米-狄拉克分布尖锐，仅允许窄能量窗口内的电子参与输运；温度升高后，更宽能量范围内的电子可同时访问多个邻近传输通道，导致离散共振峰被平均化，阶梯结构逐渐平滑。此外，P2和P3构象中更明显的线性响应可归因于更多导电通道的可用性。

**蛋白质构象对输运特性的影响**：图3展示了固定温度T=0.001时蛋白质构象对输运性质的影响。分子构象差异导致TP谱中峰的高度和能量位置发生显著变化，进而直接影响I-V曲线中阶梯的高度和宽度。P2和P3分子显示出比P1增强的电导率，表现为I-V特性中更陡的斜率和更高的饱和电流；尽管P2和P3含有相同数量的氢键，P3构象因氢键空间排列更有效促进电荷传输而表现出比P2更大的饱和电流。这些差异源于蛋白质结构内有效电子耦合的不同，其决定了可用电荷传输路径的数量和效率。DOS分析表明，尽管整体电子谱仍以有限系统中相干输运的典型离散共振态为特征，但氢键拓扑和有效电子连接性的差异导致了可测量的峰位移和强度变化；传输共振一般与DOS极大值对应，但在谱边缘存在偏差，部分电子态因弱耦合或部分局域化而未对输运做出同等贡献。

**P1构象的敏感性分析**：图4展示了对P1构象在固定温度下改变蛋白质-电极耦合强度t₃（-1.5 eV、-1.75 eV、-2.0 eV）的敏感性分析结果。t₃的变化主要影响I-V曲线的整体电流尺度，耦合越强电流绝对值越大，但特征性的阶梯状结构在所有t₃值下均保持，表明基本电子输运趋势对电极耦合的合理变化具有鲁棒性，不依赖于单一参数选择。

**讨论与结论总结**

研究讨论指出，本工作的建模选择是审慎的：采用粗粒化三维晶格表示蛋白质和鱼骨模型TB描述DNA电极，能够在保持计算可处理性的同时保留基本物理机制，从而隔离结构连通性在决定输运路径中的作用。虽然该表示不能涵盖轨道取向或长程静电效应等详细几何效应，但保留了耦合位点网络介导的相干传输基本物理。与包含多个平行输运通道、链特异性耦合和额外链间相互作用的多通道模型相比，本研究采用的单有效通道表述虽更简化，但报告的核心定性趋势——特别是传输被结构连通性调控——预期将仍然存在。研究同时承认，当前模型聚焦相干电子输运，未明确包含电子-声子耦合、溶剂相互作用或退相干过程等环境效应，这些效应在真实生物环境中会引入非弹性散射和相位破缺机制；当前结果应作为对蛋白质-DNA杂化系统中相干输运机制的定性趋势和物理认识加以解读。

研究结论部分指出，通过紧束缚框架结合朗道尔-布蒂克尔形式主义，对与DNA基电极连接的蛋白质分子中的电子输运进行了研究，分析了三种不同蛋白质构象的I-V特性和TP。在均匀热条件下，电流呈现显著的非线性响应，I-V曲线中的阶梯状特征直接对应于透射谱中峰的位置和幅度。电极温度的变化显著影响这些输运特性，逐步平滑阶梯状行为并驱使系统向更线性响应转变。不同蛋白质构象相关的电子耦合参数变化修改了I-V曲线的阶梯结构和透射峰的能量位置。对P1构象在固定温度下进行的蛋白质-电极耦合强度敏感性分析表明，尽管电流的绝对幅度随耦合强度缩放，但I-V响应的特征阶梯状特征在定性上保持不变，表明基本输运行为对电极耦合的合理变化具有鲁棒性。总体而言，这些发现突显了蛋白质-DNA杂化系统中的相干电子输运对分子构象和温度敏感，但对耦合参数的中等变化表现出稳定性。这种鲁棒性，加上对结构和热因素的已证明敏感性，强调了蛋白质-DNA结作为生物电子和纳米尺度器件应用中可调谐且可靠组分的潜力。