本研究由来自麻省理工学院（Massachusetts Institute of Technology）化学系的 Henry R. Kilgore 和 Ronald T. Raines*（通讯作者）共同完成。该研究论文于2020年4月15日在线发表于《Journal of Physical Chemistry B》期刊，并于同年刊出（DOI: 10.1021/acs.jpcb.0c02272）。
该研究属于生物物理与物理化学领域，具体聚焦于生物大分子（尤其是蛋白质）中二硫键的光物理性质。二硫键广泛存在于天然产物、多肽和蛋白质中，对其结构稳定性、折叠和功能至关重要。长期以来，科研人员观察到不同二硫键的最大吸收波长（λmax）存在显著差异（从250纳米到500纳米不等），但其背后的量子化学起源尚不完全明确，且缺乏一个预测模型。因此，本研究旨在探究二硫键发色团（chromophore）的根本起源，并建立其光吸收特性与分子几何构型之间的定量关系。研究的目标是阐明立体电子效应（stereoelectronic effects）如何调控二硫键的激发能，进而发展出一种能够根据吸收光谱预测二硫键构象，或根据构象设计特定光物理性质二硫键的方法。
详细工作流程： 本研究采用了结合实验测量与理论计算的综合方法，其流程系统而严谨，可分为以下几个主要步骤：
1. 样品准备与紫外-可见光谱（UV-Vis Spectroscopy）实验： 研究首先选用了两种典型的含二硫键分子作为实验对象：硫辛酸（lipoic acid）和二硫苏糖醇（dithiothreitol, DTT）的氧化形式。样品购自Sigma-Aldrich。实验方法是将这些二硫化物配制成10 mM的溶液（溶剂为20 mM、pH 7.5的Bis-Tris缓冲液），然后使用安捷伦Cary 60紫外-可见分光光度计记录其吸收光谱，获取实验最大吸收波长（λmax_exp）。此外，为了探究溶剂环境（介电常数）对二硫键吸收强度（摩尔吸光系数）的影响，研究还设计了混合溶剂实验：分别配制硫辛酸和二硫苏糖醇二硫化物在水和二氧六环中的1 mM溶液，然后以不同比例混合，形成具有不同介电常数的溶液体系，并测量其吸光度变化。这部分实验数据为后续理论模型的验证和补充提供了关键依据。
2. 量子化学计算与模型构建： 这是本研究的核心部分，主要依赖高斯（Gaussian）软件进行计算。所有计算均在M06-2X/6-311+G(d,p)的理论水平上进行。 * 体系选择与构象优化： 研究选择了一系列环状二硫化物作为主要计算模型，这些环状分子因其环大小的不同（包含1到6个亚甲基），其C-S-S-C二面角（θ）被限制在特定的、各不相同的值。首先对这些环状二硫化物的基态（S0）几何结构进行优化，并确认无虚频，从而得到能量最低的稳定构象及其对应的θ计算值。 * 激发态计算： 使用含时密度泛函理论（Time-Dependent Density Functional Theory， TD-DFT）计算上述优化后分子的激发性质。计算了前20个单重态激发，其中能量最低的激发所对应的波长被记录为计算最大吸收波长（λmax_comp）。 * 自然跃迁轨道分析： 为了直观理解电子从基态（S0）到第一激发态（S1）跃迁时电子密度的重组情况，研究采用了Martin的自然跃迁轨道（Natural Transition Orbitals）方法。该方法是高斯软件中实现的一种分析工具，能够提供比传统分子轨道更清晰的电子跃迁图像。使用GaussView 6软件对轨道进行可视化（等值面设为0.045）。 * 构象约束计算： 为了在更广泛的θ值范围内（而不仅仅是环状分子所限定的几个离散值）研究构象与吸收波长的关系，研究选取了无环的二乙基二硫醚（diethyl disulfide）作为模型。通过在高斯计算中约束其C-S-S-C二面角θ为一系列固定值（从-90°到90°），分别进行几何优化和TD-DFT计算，从而获得一系列（θ, λmax_comp, 构象约束能E-E0）的数据点。 * 端基异构效应分析： 为了进一步探究超共轭（hyperconjugation）等立体电子效应的影响，研究构建了一系列X取代的二甲基二硫醚（XH2C-S-S-CH2X， X = Me, F, Cl, OMe）模型。对每种分子，分别计算了有利于端基异构效应（anomeric effect，即硫原子孤对电子向相邻C-X键的σ*反键轨道捐赠电子）的构象和不利于该效应的“无关”构象的自由能，并计算两者之差（ΔG_anomeric）。同时，计算了处于端基异构效应构象时二硫化物的λmax_comp。
3. 数据分析与函数拟合： 所有获得的数据（包括实验的λmax_exp与计算的θ， 计算的λmax_comp与θ， 以及构象约束能E-E0与θ等）均使用MATLAB R2018a软件进行绘图和分析。关键的一步是将λmax与θ的关系数据拟合为一个正弦函数方程，该过程类比了核磁共振中著名的Karplus方程。
主要研究结果： 研究通过上述流程，获得了以下层层递进的重要结果：
1. 确立了C-S-S-C二面角θ与最大吸收波长λmax之间的定量正弦关系： 通过分析环状二硫化物的实验和计算数据，研究者发现，无论是实验测得的λmax_exp还是理论计算的λmax_comp，都与计算的C-S-S-C二面角θ呈现出完美的正弦曲线关系（图1a， 1b）。二者高度吻合（图1c， R²=0.98）。这直接证明了二硫键的发色团性质主要由其三维几何构型决定。研究者成功地将该关系拟合为以下经验方程（公式1）： λmax(θ) = A – B sin|θ| + C sin²|θ| （其中-90° ≤ θ ≤ 90°） 对于实验数据，拟合参数为A=460 nm， B=207 nm， C=8.42 nm。该方程类似于Karplus方程，实现了通过光学吸收光谱（λmax）来预测二硫键关键构象参数（θ）的可能性，反之亦然。例如，利用该方程预测蛋白质中常见的二硫键构象（θ ≈ ±90°）对应的λmax约为257 nm（实验）或262 nm（计算），这与已知的胱氨酸在约260 nm处的吸收峰非常一致。
2. 揭示了该关系的量子化学起源——立体电子效应： 自然跃迁轨道的可视化分析（图2， 图3a）为上述经验关系提供了直观的机理解释。分析表明，基态（S0）的最高占据分子轨道（HOMO）主要由硫原子的孤对电子组成，并受到α位（及较弱程度的β位）原子的贡献。当二硫键构象趋于重叠（θ→0°）时，硫原子孤对电子与α位取代基之间的超共轭相互作用减弱，同时硫原子孤对电子之间的电子-电子排斥增强。这两种立体电子效应的共同作用，导致基态（S0）的能量被抬高（去稳定化）。由于激发态（S1）的能量受构象影响相对较小，因此基态能量的升高直接导致了S0→S1跃迁所需能量的增加，即吸收波长λmax发生蓝移（向短波方向移动）。这一解释在二乙基二硫醚的构象约束计算中得到进一步支持：构象约束能（E-E0， 代表偏离最稳定构象的能量代价）的变化趋势与激发能（即1/λmax）的变化趋势完全同步（图3b）。这表明，任何通过构象或取代基扰动二硫键局部电子环境、从而改变其基态稳定性的因素，都会直接影响其光吸收性质。
3. 验证了超共轭（端基异构效应）对吸收波长的调控作用： 在取代的二甲基二硫醚模型中，研究者发现，当分子采取有利于端基异构效应的构象时，其基态因孤对电子向σ*C-X轨道的捐赠而得到稳定，从而需要更高的能量才能激发。计算结果显示，端基异构构象与非端基异构构象的自由能差（ΔG_anomeric）与处于端基异构构象时分子的λmax_comp存在明确的负相关（图4b）：ΔG_anomeric越负（即端基异构稳定化效应越强），λmax_comp越短（蓝移越大）。自然跃迁轨道图像（图4c）清晰地展示了在端基异构构象下，电子密度从硫原子向C-X反键轨道的重新分布。这直接证明了除二面角外，通过引入特定的取代基（如吸电子的F、Cl）来增强超共轭效应，是调控二硫键光物理性质的另一有效手段。
4. 发现了吸光度与溶剂极性的关系： 补充实验表明，二硫键的摩尔吸光系数（而非最大吸收波长）受溶剂环境影响。在极性溶剂中，硫辛酸和DTT二硫化物的吸光度比在非极性溶剂中降低约一倍（图S2c）。研究者推测这可能源于溶剂对跃迁偶极矩的微扰。这一发现提示，在分析蛋白质疏水核心中的二硫键时，其在非极性环境中的较高吸光度可能更有利于检测。
结论与意义： 本研究的结论是：二硫键的光物理性质是可预测且可调控的。其最大吸收波长源于电子-电子排斥和超共轭等立体电子效应，这些效应通过C-S-S-C二面角以及α位取代基来调节基态（S0）的能量。研究者建立的类Karplus正弦方程（公式1）为实现这一预测和调控提供了简单而有力的工具。
该研究的科学价值在于，首次从量子化学层面清晰阐明了二硫键发色团巨大波长变化范围的结构起源，将宏观的光谱现象与微观的分子轨道和立体电子效应直接联系起来。其应用价值广泛：在基础研究方面，为生物化学家提供了一种通过常规紫外光谱便捷表征蛋白质或肽链中二硫键构象（特别是C-S-S-C二面角）的新方法；在合成与材料化学方面，为设计具有特定光吸收性质（例如，用于光裂解、光传感或作为光谱探针）的二硫键分子提供了理论指导原则。
研究亮点： 1. 重要发现： 建立了二硫键C-S-S-C二面角与其最大吸收波长之间精确的正弦函数定量关系，填补了该领域长期存在的知识空白。 2. 方法新颖性： 创造性地将核磁共振构象分析中的Karplus方程思想引入到紫外光谱对二硫键构象的表征中，实现了跨学科方法的迁移应用。 3. 机理阐释深刻： 综合运用TD-DFT计算和自然跃迁轨道可视化，清晰揭示了立体电子效应（电子排斥和超共轭）调控激发能的物理图像，使经验关系具备了坚实的理论根基。 4. 研究系统性： 从环状模型到无环模型，从几何效应到电子效应（端基异构），从气相计算到溶液实验，研究设计层次分明，证据链完整，逻辑严密。
其他有价值的内容： 该研究在引言部分回顾了自20世纪50年代Melvin Calvin开创性工作以来的相关研究历史，将当前工作置于经典研究的脉络之中，体现了科学的传承与发展。此外，研究将论文题献给Martin Karplus教授90岁生日，并在文中多次与Karplus方程类比，也体现了对学术前辈的致敬和对科学方法普适性的思考。支持信息中包含了详尽的实验数据表、补充光谱图和计算所得的原子坐标，确保了研究的可重复性和可验证性。