苯基丙氨酸代谢物（苯乙酸、苯乳酸与苯丙酮酸）在稀水体系中的结构-动力学行为特征：一项分子动力学研究

《Chemical Thermodynamics and Thermal Analysis》：Features of the structural-dynamic behavior of phenylalanine metabolites (phenylacetic, phenyllactic and phenylpyruvic acids) in dilute aqueous systems: a molecular dynamics study

　　

本研究采用分子动力学模拟考察苯乙酸（PA）、苯乳酸（PL）与苯丙酮酸（PP）的结构特征对其在水体系中结构与动力学行为的影响，并建立苯基丙氨酸代谢物的动力学与其在苯酮尿症中致病效应之间的关系。代谢物的亲水性按 PA < PL < PP 顺序递增，其中苯乙酸盐（PA）表现出最显著的疏水特性，而苯丙酮酸盐（PP）表现出亲水特性。研究表明苯乳酸盐（PL）呈现超扩散（super-diffusion）机制，这与其在苯酮尿症中结构损伤作用较弱相关；与之相对，PA 与 PP 以缓慢的亚扩散（sub-diffusion）模式运动，可能促进代谢物的局部累积并聚集成神经毒性纤维。研究人员在 307 K 下对高度稀释至单体互不作用的 PA、PP、PL 水溶液进行模拟，借助径向分布函数（RDF）、Wendt-Abraham 参数 RWA、均方位移（MSD）及自扩散系数等分析，阐明苯环与 –CH2– 决定疏水区，R2 基团（PA 的 –OOH–、PL 的 –C2O(OH)2–、PP 的 –C2O2(OH)–）决定亲水区及氢键（HB）网络差异，进而调控聚集路径：PA 与 PP 因慢速亚扩散与 π–π 堆积倾向成纤维，PL 因超扩散无有序水合结构而不形成明确聚集。该工作为苯酮尿症代谢应激下小分子自组装的物化链路提供了单分子水合视角。

研究背景与立项依据

苯酮尿症（PKU）患者因苯丙氨酸羟化酶不足，苯丙氨酸转而积累为苯丙酮酸（PP）等代谢物，并伴随苯乙酸（PA）、苯乳酸（PL）等衍生物在体液中异常富集。已有实验证据提示此类小分子可自组装为纳米纤维并参与代谢病发病，但单体结构如何决定其在水中的水合壳层、扩散异常与后续聚集，尚缺乏从单分子水平出发的定量描述。研究人员以 PA、PL、PP 三种含苯环芳香酸为对象，在无限稀释条件排除分子间聚集干扰，专注溶质–溶剂相互作用与动力学异质性，试图把分子结构差异、亲水/疏水分配、异常扩散机制与 PKU 神经毒性的结构损伤联系起来。论文发表于《Chemical Thermodynamics and Thermal Analysis》。

主要技术方法

研究人员采用 DL_POLY_4.0 执行经典分子动力学（MD），溶质参数基于 OPLS 力场，水采用 SPC/E 势；在 NPT 与 NVT 系综下分步平衡后生产约 6 ns 轨迹，温度 307 K，体系为无限稀释单分子水溶液。结构分析用径向分布函数 G_XY(r) 与 Wendt-Abraham 参数 R^WA 判定水合壳层与局域序；动力学分析用均方位移（MSD）指数 α、时间平均均方位移（TAMSD）、自散射中间函数 F_s(q,t) 及其导数、速度自相关函数（VAF）经 Green–Kubo 关系求自扩散系数 D。无生物样本队列，纯计算研究。

研究结果

3.1 PA、PP、PL–水体系的局域结构

通过 G_XY(r) 分析，苯环碳 CA 与水氧 O^W 第一峰约 3.8 ?、自由空间延至约 6.0 ?；–CH₂– 区自由空间 3.5–6.0 ?，确认 –C₆H₅– 与 –CH₂– 为共同疏水来源且水合半径与 R2 无关。R2 基团区：PA 的 O^PA/OH^PA、PL 的 O^PL/O^OAL/OH、PP 的 O^PP/OH^PP 均形成 ~1.8 ? 氢键，水合壳层数随基团极性与羰基数增加。R^WA 沿 O^PA(0.53)→O^PP(0.89)→O^PL(0.91) 与 OH^PP(0.345)→OH^PA(0.69)→OH^PL(0.92) 增大，表明 PP、PL 羟基氧侧局域水序更低、亲水性更强；综合得亲水性 PA < PL < PP。

3.2 PA、PP、PL–水体系的动力学异质性

MSD 显示 0–25 ps 三者均处振动–跳跃扩散（Frenkel 框架）；>25 ps 分化：PA 与 PP 的 α<1 为亚扩散，PL 在 >80 ps 转 α>1 超扩散。F_s(q,t) 导数在 3 ps 与 6 ps 出现机制切换。自扩散系数：水的 D 随 PL(5.2×10^–9)→PP(4.9×10^–9)→PA(4.7×10^–9) m²·s^–1 微降（纯水 3.14×10^–9）；溶质 D 为 PP(0.22×10^–9)< />^–9)< />^–9) 错误顺序应修正为 PP(0.22)< />< />^–9 m²·s^–1。PP 出现多段随机行走 t_RW（10–14 ps）与水合液态笼形结构，运动迟滞；PL 跨接区间最长（~125 ps）且超扩散，不易局域累积。

讨论与结论翻译

讨论指出，亲水端经水层氢键“胶连”、疏水苯环经 π–π 堆积，使 PA 与 PP 在较低浓度即可成纤维（PA 阈值 0.6 mg·mL^–1），PP 先成胶束再延展纤维；PL 两亲性弱且无序水合，不形成明确聚集，对应其在 PKU 中结构损伤角色较弱。慢速亚扩散延长 PA/PP 在局域液体环境的滞留时间，促进代谢物局部累积与神经毒性自组装；超扩散的 PL 则持续重排水网、避免有序装配。

结论部分：RDF 与 R^WA 表明亲水性由 R2 基团结构决定，疏水性由 –C₆H₅– 与 –CH₂– 决定；亲水性序 PA < PL < PP。PA 与 PP 以迟滞亚扩散运动，PL 以超扩散运动，前者利于 π–π 驱动纤维生成，后者因无有序水合壳而不聚集。从异常亚扩散到有序分子组装（胶束、纤维）的转变为苯酮尿症代谢应激提供了物化链接。