溶剂力与对称性破缺对设计蛋白在液-固界面组装的影响

研究人员旨在解决当前蛋白质设计平台无法准确预测蛋白在无机表面组装结果的问题。他们通过高速原子力显微镜（HS-AFM）结合机器学习，结合蒙特卡洛模拟，研究了溶剂力及基底对称性如何影响设计蛋白纳米棒在云母表面的组装。研究发现，界面溶剂结构（如水化层）的对称性破缺，为原本非相互作用的棒状系统引入方向性偏倚，从而稳定了意想不到的二维近晶相。这项成果强调了在设计蛋白-无机界面杂化材料时，必须考虑溶剂力及基底调制的水化结构，为改进蛋白质设计平台提供了关键见解。

在生物材料的世界里，蛋白质扮演着“建筑大师”的角色，能够在纳米尺度上进行精妙的自我组装，构建出具有复杂功能的杂化（hybrid）材料。受此启发，科学家们希望像设计积木一样，从头（de novo）设计蛋白质，让它们按预设的蓝图在无机晶体表面有序排列。这被称为蛋白质设计时代，有望在能源收集、催化、光电子和生物医学等领域带来革命性的应用。然而，现实往往比蓝图复杂。即便研究人员精心设计了蛋白质之间的化学键、静电力或疏水作用，当这些“设计蛋白”真正在固体表面（如云母）上安家落户时，它们常常不按“剧本”出牌，组装出许多意想不到的图案。例如，一种名为RhuA的酶，在溶液里能完美形成预设的二维（2D）晶格，一旦放在云母上，却产生了三种完全没预料到的单层和双层结构。这就像一个建筑队在图纸上画好了完美的房屋结构，一到实地却盖出了风格迥异的建筑。问题出在哪里？显然，当前的设计平台漏掉了一些关键的“幕后推手”。

研究人员敏锐地意识到，这些“不听话”的组装结果，暗示着除了预先设计的蛋白质-蛋白质、蛋白质-基底相互作用之外，另一些力在暗中左右大局。蛋白质虽然结构精确，但在界面上更像一个个带静电“补丁”的胶体颗粒，与周围的溶液和基底都在互动。特别是界面处的溶剂（水）结构，其排列方式受到基底晶体对称性的深刻影响，从而产生所谓的溶剂力。这些力是否就是导致“脱靶”组装的神秘力量？为了解答这个疑问，一个研究团队在《Nature Communications》上发表了一项研究，他们深入探究了溶剂力及基底对称性破缺对设计蛋白质在液-固界面组装的决定性影响。

为了开展这项研究，作者们运用了几个关键的技术方法。他们利用高速原子力显微镜（HS-AFM）在液相中原位、实时地观测蛋白质纳米棒在云母表面的动态组装过程。通过机器学习算法对海量的HS-AFM图像进行自动分析，实现了对单个蛋白纳米棒的位置、取向的精准识别与追踪，从而量化组装动力学和有序度。此外，他们建立了二维巨正则蒙特卡洛（Monte Carlo）模拟，将蛋白质简化为硬棒模型，通过系统地改变化学势、棒状物迁移率和取向偏置，来预测和解释实验中观察到的不同有序相。研究所用的核心蛋白是名为DHR10-mica18的矩形纳米棒状设计蛋白，其通过Rosetta平台设计，旨在与云母的K⁺亚晶格实现静电匹配。

研究结果

组装了DHR蛋白在云母上

研究人员使用了一种长宽比为1:5.6的矩形棒状从头设计螺旋重复（DHR）蛋白DHR10-mica18。该蛋白被设计为通过其54个谷氨酸残基的羧酸根基团与云母（001）解理面上的K< />位点产生静电相互作用。他们选择了两种云母：白云母（m-mica）和氟金云母（f-mica）。两者具有完全相同的K⁺亚晶格，但晶体结构细节不同，这导致其上方的水化层结构存在差异：f-mica上的水化层保持六方对称，而m-mica由于对称性破缺，其第二水化层呈现出条纹状图案。

实验发现，在100 mM KCl溶液中，无论在f-mica还是m-mica上，蛋白纳米棒都组装成一种无序相，由沿着云母晶格三个主轴方向排列的小区域组成。然而，当K⁺浓度升至3 M时，情况发生了分化：在f-mica上，蛋白质仍形成三方向无序相；但在m-mica上，所有纳米棒都沿着m-mica唯一轴的方向完全对齐，并排列成平行的行，在整个基底上形成了高度有序的二维近晶相（smectic phase）。近晶相的出现令人惊讶，因为根据长期以来的胶体系统理论，对于非相互作用的二维矩形棒，仅有向列相（nematic phase）是稳定的，近晶相不应出现。

原位高速AFM结果和机器学习分析以追踪蛋白纳米棒在m-mica上的组装

为了量化组装动力学，研究人员利用HS-AFM原位观察了m-mica在3 M KCl下的组装过程。初始阶段，蛋白质形成高迁移性的二维液相，随后逐渐出现短暂存在的小型对齐纳米棒域，最终这些域扩大、稳定，形成近晶相。通过基于机器学习的图像分析工作流程，他们追踪了每个蛋白质的位置和取向，计算了向列序参数和近晶序参数。分析显示，向列序参数迅速上升并饱和，表明几乎所有纳米棒都沿单一方向排列。而近晶序参数则缓慢增加，在表面覆盖度接近饱和时快速跃升，表明当域紧密排列时达到了渗流阈值，从而实现了高度的平移有序。

原位高速AFM结果和机器学习分析以追踪蛋白纳米棒在f-mica上的组装

相比之下，在f-mica的3 M KCl条件下，初始阶段就能观察到沿三个K⁺亚晶格方向随机取向的单个纳米棒和小域。这些域会波动性生长，但最终形成的是高密度无序相，其向列序参数始终很低，近晶序参数几乎为零。这表明在f-mica上，蛋白质的迁移率较低，使其被“动力学捕获”在无序相中，无法达到平衡态的向列相。

巨正则蒙特卡洛模拟有无偏置的硬棒

为了探究观察到的现象是否纯粹由胶体力引起，研究人员进行了二维巨正则蒙特卡洛模拟。他们模拟了硬棒在两种不同势能景观下的沉积：一种是三方向对称的势（模拟f-mica），另一种是其中一个取向能量上比其他两个有利两倍的准二重对称势（模拟m-mica的对称性破缺）。模拟系统地改变了化学势和棒迁移率。

结果显示，在三重对称势下，模拟结果与f-mica实验一致：在零迁移率时得到无序相，在非零迁移率和高化学势下会出现向列相，但近晶序参数始终近乎为零。然而，当引入哪怕只是两倍的取向偏置（模拟m-mica），只要棒有足够的迁移率和化学势，清晰的近晶序就会出现。这一关键模拟结果证明，实验中观察到的近晶相并非由棒-棒相互作用（如电荷排斥）稳定，而是由基底对称性破缺（通过界面水结构）在势能景观中引入的二重偏置所导致。

结论与讨论

这项研究清晰地表明，设计蛋白质在表面的二维组装并非仅由预设的蛋白质-基底相互作用所主导。与之同等重要的，是源自溶剂-胶体相互作用的熵力。更重要的是，这些胶体力不同于体相溶液中的力，因为基底的对称性被“烙印”在了界面附近的溶剂结构上。具体而言，DHR10-mica18蛋白纳米棒在3 M K⁺条件下的m-mica上组装出近晶相，正是精心设计的与K⁺亚晶格三方向匹配的相互作用，与由m-mica晶格对称性破缺所导致的界面溶剂结构引入的二重偏置，两者共同作用的结果。这种组合效应稳定了原本不可能存在的近晶相。

这项发现具有重要的方法论意义。它指出了一个改进从头蛋白质设计平台（如Rosetta）的分子级补充方案。未来的设计流程需要整合能考虑界面水化结构非均匀性的粗粒度模拟。例如，可以针对每种新的晶体表面，先通过分子动力学或三维AFM确定其水化结构。然后将蛋白质设计粗粒化，捕捉其关键分子特征，并使用能够结合溶剂排阻效应和界面水化结构特异性的计算代码来预测真实的组装结果。如果预测结果与设计目标不符，这一偏差可作为神经网络的损失函数，反过来指导修改蛋白质设计，经过迭代最终收敛到能够实现目标组装的设计。这种将物理模拟与蛋白质设计算法耦合的思路，有望引领开发出更强大的、用于设计生物启发杂化材料的平台，使得在纳米尺度上精确“编程”蛋白质在无机界面的建筑成为可能。