《Biochimica et Biophysica Acta (BBA) - Biomembranes》：Pore-level quantitative structure-activity relationship (QSAR) for water permeation rate in aquaporins
编辑推荐：
水通道蛋白(AQ(G)Ps)是水跨膜运输的关键蛋白，其渗透机制的理解对设计高性能仿生膜至关重要。本研究针对AQ(G)Ps水渗透率的结构决定因素不明确这一挑战，创新性地构建了孔道水平的定量构效关系(QSAR)模型。通过整合多种孔道描述符与实验水渗透率数据，研究不仅验证了单因素分析中已知的决定因素，更在多变量背景下揭示了新的孔道描述符，为理性设计基于AQ(G)Ps的过滤装置提供了新见解，有助于加深对这些重要大分子在健康和疾病中功能的理解。
在生命体中，水分子高效而选择性地穿越细胞膜，是维持细胞和机体稳态的基础生命过程。承担这一关键任务的主角，是一类名为水通道蛋白(Aquaporins, AQPs)和水甘油通道蛋白(aquaglyceroporins, AQGPs)的膜蛋白。它们不仅是维持植物耐盐、耐旱、动物细胞体积调节、神经兴奋等生理功能的核心，其异常表达还与多种人类非传染性和传染性疾病相关。除了生物学上的重要意义，AQ(G)Ps以其近乎完美的渗透效率和选择透过性，成为了设计下一代仿生过滤膜的绝佳蓝图，尤其是在海水淡化、废水处理等领域展现出替代传统反渗透(RO)或热脱盐技术的巨大潜力，有望实现显著的节能和成本节约。然而，尽管AQ(G)Ps的研究已持续数十年，其渗透特性和选择性机制的精确、全面理解仍然不完整。以往的研究大多聚焦于单个结构特征，如NPA、SF（Selectivity Filter, 选择性滤器）基序的氨基酸组成、孔道形状大小、孔道衬里氢键数量等。如何将这些分散的、单因素的研究成果整合起来，在一个多变量的框架下系统地揭示决定AQ(G)Ps水渗透效率的关键结构特征，是当前领域面临的一个重要挑战。
为了回答上述问题，来自意大利国家研究委员会化学科学与技术研究所“Giulio Natta”的研究团队，包括Juan José Galano-Frutos、Luca Bergamasco、Paolo Vigo、Matteo Morciano、Matteo Fasano、Davide Pirolli、Eliodoro Chiavazzo和Maria Cristina de Rosa，在《Biochimica et Biophysica Acta (BBA) - Biomembranes》期刊上发表了题为“Pore-level quantitative structure-activity relationship (QSAR) for water permeation rate in aquaporins”的研究。该研究首次将广泛应用于药物化学的定量构效关系(QSAR)方法，创新性地应用于AQ(G)Ps的孔道水平分析。研究人员构建了一个概念验证性的、数据驱动的方法学框架，旨在建立一组具有高分辨率三维结构的AQ(G)Ps的孔道结构与物理化学描述符，与其实验报告的水渗透率之间的定量关系。
研究人员为开展此研究，主要采用了以下几个关键技术方法：(1) 结构数据收集与处理：从PROSITE和InterPro数据库检索主要内在蛋白(Major Intrinsic Proteins, MIPs)家族成员，并从RCSB PDB档案中获取了22个具有实验三维结构的代表性AQ(G)Ps蛋白，为每个蛋白选择最高分辨率的、无突变的结构进行分析。(2) 渗透率数据汇编：系统回顾文献，从实验和分子动力学模拟(Molecular Dynamics Simulations, MDS)研究中，汇编了这些AQ(G)Ps的单通道水渗透率(pf)数据。特别地，他们还整理了由Wachlmayr等人近期提出的、经过方法学修正的校正渗透率(pfcorr)数据集。(3) 孔道描述符计算：为每个AQ(G)Ps结构计算了13个孔道描述符，涵盖几何特征（平均直径、最小直径、三维形状规则性）、亲/疏水性、基于晶体学B因子的局部柔性指标、氢键位点数量以及孔道前庭区域的静电/溶剂化特征。计算使用了PoreWalker服务器、Discovery Studio软件以及基于Horner等人方法开发的自定义Python脚本。(4) QSAR建模与验证：针对包含实验pf值（n=10）和校正pfcorr值（n=8）的两个训练集，采用LASSO回归和普通最小二乘法进行多变量线性回归建模。通过内部交叉验证、方差膨胀因子(VIF)分析、y-随机化检验来评估模型的稳健性和预测能力，并利用四个缺乏实验结构但拥有pfcorr值的AQPs的AlphaFold2预测结构作为外部测试集，对优选模型进行验证。
结果
3.1. 关于AQ(G)Ps的渗透率
研究首先系统汇编了水渗透率数据。分析发现，实验测得的pf值与分子动力学模拟(MDS)估算的值一致性较差。此外，基于Wachlmayr等人对传统停流(stopped-flow)方法学偏误的修正，许多历史报告的pf值可能被高估。因此，研究人员分别针对未经修正的实验pf数据集和经过方法学修正的pfcorr数据集进行了独立的建模，以评估数据校正对模型性能和物理解释性的影响。
3.2. 孔道描述符分析
研究人员计算了广泛的孔道描述符。例如，通过PoreWalker分析获得了孔道直径沿轴向(z-coordinate)的分布图，并通过将所有AQ(G)Ps结构的孔道以其第二个NPA基序中天冬酰胺的β碳为参考点(z=0)进行对齐，实现了孔道坐标的均一化处理。计算了孔道衬里残基的平均疏水性/亲水性指数、选择性滤器(SF)区域（z坐标2-6 ?）内残基的平均B因子（作为局部柔性的间接指标）、单列水线区域潜在的氢键位点数量(NH)，以及孔道前庭内带正/负电荷的残基数(Res+/Res-)和带电原子位点数(Sites+/Sites-)。
3.3. 针对实验pf的孔道水平QSAR建模
对包含10个AQ(G)Ps的实验pf数据集进行建模。经过变量筛选和多重共线性检验，最终得到一个包含三个孔道描述符的稳健双变量模型。该模型具有良好的拟合优度(R2= 0.88)和内部预测能力(loo-Q2= 0.76)。模型表明，更高的水渗透率与更大的平均孔道直径和孔道衬里更高的平均疏水性正相关。这一发现与之前的单因素研究相符，即更宽的孔道有利于水传输，而疏水环境可能减少了水分子与孔道壁之间过强的相互作用，从而降低了传输能垒。
3.4. 针对校正pfcorr的孔道水平QSAR建模
对包含8个AQ(G)Ps的校正pfcorr数据集进行建模。经过更严格的筛选，得到一个同样包含三个描述符但构成不同的最优三变量模型。该模型具有优异的拟合优度(R2= 0.99)和内部预测能力(loo-Q2= 0.95)。此模型揭示的决定因素与前一模型有显著不同。它表明，更高的校正渗透率与更小的平均孔道直径、选择性滤器(SF)区域更高的平均B因子（暗示更高的局部柔性）以及孔道前庭区域更少的带正电荷位点(Sites+) 相关联。这个模型不仅成功复现了先前单因素研究中关于孔道尺寸和静电环境对渗透率影响的结论，更重要的是，它首次在QSAR框架下，将SF区域的局部构象柔性与水渗透效率定量地联系了起来。研究人员推测，SF区域适当的柔性可能有利于水分子的重新定向和快速通过。
3.5. 使用缺乏实验结构的AQPs进行外部验证
为了测试模型的泛化能力，研究人员选取了四个拥有pfcorr实验值但缺乏实验三维结构的AQPs（来自不同物种的AQP3, AQP9, PsPIP2;1, AtPIP2;1），使用AlphaFold2预测了它们的结构，并计算了优选三变量模型所需的描述符。将描述符代入模型后，预测的ln(pfcorr)值与实验报告的ln(pfcorr)值表现出良好的一致性，皮尔逊相关系数达到0.94。这一结果强有力地支持了该QSAR模型的预测能力和可靠性，表明该方法即使对于缺乏高分辨率晶体结构的同源蛋白也具有一定的应用潜力。
结论与讨论
本研究成功开发并验证了一个针对AQ(G)Ps水渗透率的孔道水平QSAR方法学框架。作为概念验证，该研究展示了如何将计算结构描述符与实验渗透率数据相结合，以揭示驱动水传输效率的关键、且可能是新颖的孔道特征。
研究的主要结论和重要意义在于：
- 1.
方法论创新：首次将QSAR建模系统性地应用于AQ(G)Ps功能研究，为理解这类复杂膜蛋白的渗透性质提供了一个全新的、数据驱动的、多变量的分析视角。
- 2.
关键决定因素的再发现与新发现：通过建模，不仅确认了已知的几何（孔径）和静电（前庭电荷）因素的重要性，更重要的是，首次在QSAR模型中识别出SF区域的局部构象柔性是水渗透率的一个关键决定因素。这为理解AQ(G)Ps的功能机制增添了新的维度。
- 3.
数据质量的重要性：研究分别基于原始实验数据和经过方法学修正的数据构建模型，得到了物理解释迥异的结论。这突显了精确、可靠的实验测量数据对于构建有意义的预测模型和获得可靠生物学见解至关重要。使用校正数据集得到的模型具有更优的统计性能和更合理的物理解释，这为未来相关研究中的数据标准化和处理提供了参考。
- 4.
模型具有良好的预测与泛化能力：优选模型在内部验证中表现优异，并且能够较为准确地预测一组外部测试蛋白（基于AI预测结构）的渗透率，证明了模型的稳健性和潜在应用价值。
- 5.
为理性设计提供指导：该研究所建立的多特征框架和识别出的关键孔道描述符，为理性设计和优化基于AQ(G)Ps的仿生过滤膜(ABMs) 以及设计人工水通道(AWCs) 提供了分子水平的见解和潜在的设计靶点。例如，可以通过工程化手段微调SF区域的氨基酸组成以调节其柔性，或优化孔道前庭的电荷分布，从而可能定制出具有特定渗透性能的过滤材料。
- 6.
深化分子机制理解：从更基础的生物学角度看，这项研究有助于深化对AQ(G)Ps在健康和疾病状态下分子功能的理解，可能为针对这些通道的疾病干预策略提供新思路。
总之，这项研究不仅为AQ(G)Ps的研究领域引入了一种有力的新分析方法，其得出的具体结论也为未来高性能仿生膜的设计和AQ(G)Ps的生物学功能调控开辟了新的可能性。
