网中之缕：糖类、腈类与唑类选择性前生物反应自由能的计算研究

《Agronomy》：Strands in a Web: Free Energies of Select Prebiotic Reactions with Sugars, Nitriles and Azoles

　　

研究人员利用计算化学方法，测定了C₂和C₃糖类分别与HCN和NH₂CN反应生成氨基酸和唑类的自由能变化。氨基酸和唑类可进一步与初始糖类反应，生成更复杂的分子，从而加入具有众多物种的前生物化学复杂网络。对于所选择的反应，这些计算提供了一个尚不完整的指南，用以判断哪些途径在热力学和/或动力学上有利，以及从少数简单反应物出发可能观察到哪些物种。所研究的反应是网络中的缕缕支线，将前生物化学的不同领域与生命起源化学这一更大网络连接起来。

生命化学起源的核心问题已从“能否由简单小分子合成生命相关有机物”转变为“为什么生物系统仅利用前生物合成所产生的众多化合物中的一小部分”。评估反应的热力学与动力学可为这一难题提供线索。甲酰糖(formose)体系通过简单的醛醇(aldol)反应由甲醛生成多种糖类，但生命利用的仅是其中少数；而HCN等腈类因具有亲核碳，可与糖类经由Strecker机制(Strecker mechanism)生成氨基酸。研究人员此前已构建甲醛和HCN、甲酰糖循环等局部自由能图，但更大的前生物反应网络中许多关键连接仍缺乏统一的能量数据。为此，本研究选取网络中的四条支路，考察糖类、腈类、唑类和氨基酸之间反应的自由能与能垒，以判断哪些通路在热力学和/或动力学上可行，以及从简单反应物出发可能观察到哪些物种。

研究人员采用与既往工作一致的计算化学方案。量子化学计算使用Q-Chem 6.2，在B3LYP-D3/6-311G**水平优化几何并计算电子能，用SMD隐式溶剂模型(implicit solvent model)计算水相溶剂化能；热力学修正取298.15 K，并参照Deubel与Lau的方法将溶解熵设为气相熵的一半。过渡态(transition state)优化加入显式水分子辅助质子转移。所有物种的相对自由能G_rel以CO₂、H₂O、H₂、H₂S和NH₃为基准。该方案与乙二醇醛和甲基乙二醛的实验结果一致，水相相对自由能误差约0.5 kcal/mol；涉及羰基化学的能垒系统误差为2–3 kcal/mol。

### 3.1 糖类与HCN生成氨基酸
以C₂糖乙醇醛(1a)与HCN/NH₃经Strecker途径生成丝氨酸，其硫代类似物巯基乙醛(1b)对应生成半胱氨酸。计算表明，第一步加NH₃放能4.6 kcal/mol，能垒仅1.2 kcal/mol；脱水成亚胺吸热4.2 kcal/mol，能垒18.6 kcal/mol；HCN加成是主要热力学驱动步，放能15.4 kcal/mol，能垒7.7 kcal/mol；后续腈水合为决速步，能垒34.8 kcal/mol。总反应1a→丝氨酸放能24.0 kcal/mol，1b→半胱氨酸放能24.8 kcal/mol，与甘氨酸合成相近。

### 3.2 糖类与氨腈生成唑类
乙醇醛与NH₂CN经加成、关环、互变异构、脱水生成噁唑；巯基乙醛对应生成噻唑。第一步加成放能1.3 kcal/mol，能垒17.4 kcal/mol；关环放能11.1 kcal/mol，能垒23.8 kcal/mol；后续步骤中最高能垒约24.8 kcal/mol。总体放能分别为噁唑13.6 kcal/mol、噻唑14.6 kcal/mol。硫代类似物三步能垒较氧类似物低3–5 kcal/mol，因此噻唑形成在动力学上可能更容易。

### 3.3 糖–唑反应
D-甘油醛与噁唑经加成、关环生成呋喃型氨基噁唑啉，总放能9.8 kcal/mol；计算所得非对映体能量差约1.7 kcal/mol，与实验产物44:30:13:8的分布一致。竞争性六元环吡喃型产物平均比呋喃型不稳定约2 kcal/mol，平衡比例约为10:1，与实验相符。另一竞争反应是噁唑氨基对醛的亲核加成，动力学上先形成半缩醛胺，但热力学上最终仍以氨基噁唑啉为主。乙醇醛与噻唑的醛醇式加成为吸热(+1.4 kcal/mol)且能垒高(+26.4 kcal/mol)，与实验观察到噻唑可作为动力学捕获剂、不产生额外复杂产物的结果一致。

### 3.4 脱水甘油醛与半胱氨酸偶联生成LacCys
甲基乙二醛与L-半胱氨酸经半硫缩醛(HTA)、烯醇化后形成硫酯，再经环化、开环得到LacCys。初始加成放能约10 kcal/mol，能垒3 kcal/mol；硫酯形成总体放能1 kcal/mol，但能垒约26 kcal/mol，为决速步；随后环化吸热4 kcal/mol，开环放能18 kcal/mol。计算还发现，以氨基作为亲核体的副反应路径能垒高达30 kcal/mol，难以显著进行；多个内酯副产物中，五元环内酯在热力学上接近平衡，可能降低LacCys产率。计算得到的速率与实验速率常数大致在同一数量级。

研究人员在讨论中指出，该方案只计算中性非两性离子物种，未考虑pH、浓度、结晶等效应，熵修正和隐式溶剂模型均为简化处理，密度泛函与基组级别有限；但与现有标准条件下的实验结果相比，误差可接受，这得益于误差抵消。结论认为，四条支路分别展示了C–C键形成的新模式，包括HCN加成和醛醇式糖–唑加成、唑类的动力学捕获，以及硫酯中间体在前生物代谢偶联中的潜在作用。这些结果扩展了此前自由能图的覆盖范围，为实验研究提供基线热力学和动力学参数，有助于从前生物化学复杂网络中筛选出通向生命起源的可行途径。