EasyHybrid是一款专为混合量子化学/分子力学（QC/MM）模拟而设计的免费开源图形界面工具，其构建基于pDynamo3库。该软件为研究人员提供了一个直观的环境，用于准备、检查和编辑分子系统，同时支持多种模拟类型，包括反应坐标扫描、分子动力学（MD）、正常模式分析、受约束弹性带（NEB）和伞形采样（umbrella sampling）。EasyHybrid的主要特点是先进的3D可视化功能，支持对大规模生物分子系统的交互式编辑、灵活的原子选择、系统剪枝以提高QC/MM模拟的效率，以及轨道和静电势面的可视化。此外，它还具备自动日志解析和轨迹分析等功能。EasyHybrid将这些工具整合到一个平台上，为量子化学和混合QC/MM模拟提供了一个既熟悉又专业的环境。
在计算化学和分子建模领域，交互式图形界面扮演着至关重要的角色。这些工具不仅仅是简单的可视化工具，它们还提供了分子绘制与编辑、文件格式转换以及生成和提交模拟输入文件等核心功能。目前，已经开发出多种图形工具来满足这一领域不同的需求。例如，wXMacMolPlt、ECCE和GaussView是专门为GAMESS-US、NWChem和Gaussian等量子化学软件设计的界面。而像Avogadro和Gabedit这样的通用分子编辑器则支持多个计算化学软件。Molden也是一个广泛使用的图形界面，具备良好的兼容性。尽管这些工具在处理大规模分子系统时可能面临一定的性能挑战，但它们在结构生物学和分子动力学数据可视化方面表现出色。然而，它们并不特别适合分子编辑任务，也缺乏对量子化学数据的可视化支持。
在处理混合QC/MM系统时，这些工具的局限性变得更加明显。混合QC/MM系统通常由大量使用经典力场（如AMBER或CHARMM）描述的原子组成，而其中一小部分则采用量子化学方法进行计算。目前，还没有一款免费且开源的交互式环境专门用于编辑、操作、执行和分析这类混合系统。虽然商业工具如Schrödinger的Maestro配合QSite模块提供了集成支持，但其高昂的成本限制了其在学术界的普及。另外，一些基于网络的界面，如CHARMM-GUI，虽然在QC/MM输入文件的准备和系统设置期间提供了结构可视化支持，但它们的功能范围有限。因此，EasyHybrid的出现填补了这一空白，为学术界提供了一个免费、开源且全面的平台，支持pDynamo3生态系统中的先进方法。
EasyHybrid的界面采用Python3语言实现，并使用GTK3工具包生成图形窗口。其交互式3D可视化区域作为一个GTK3小部件，由一个名为VISMOL的Python3模块开发而成，该模块与EasyHybrid一同发布，但由同一开发团队维护作为独立项目。这种模块化设计使得VISMOL可以轻松集成到任何基于GTK3的容器应用程序中，为开发者提供了一个灵活的解决方案，用于将分子3D可视化功能嵌入到自己的工具中。VISMOL利用现代OpenGL（版本3.6）进行渲染，除了常用的片段和顶点着色器外，还引入了几何着色器，从而在特定的渲染模式下显著提升性能，特别是在“线”和“棍棒”表示方式中。虽然本文未深入探讨VISMOL的渲染策略，但更多相关信息和源代码可在其官方GitHub仓库中找到。
EasyHybrid的主窗口集成了六个关键组件：一个用于所有界面功能的菜单栏、一个包含常用操作快捷方式的工具栏、一个列出已加载系统的侧边树形视图、一个位于底部的包含操作日志和残基查看器的面板、一个用于总结系统属性的状态栏，以及一个位于中央的交互式3D画布。如图1所示，该界面提供了对混合QC/MM系统的全面管理能力。用户可以在同一个会话中管理并操作多个系统，每个新加载的系统都会自动添加到主树形视图中，并分配唯一的标题和颜色代码。这些颜色代码默认应用于可视对象中的碳原子，有助于用户快速识别哪些原子属于哪个系统。用户可以通过树形视图中的单选按钮在不同系统之间切换，实现针对性的编辑和模拟。可视对象可以来源于模拟输出或从外部坐标文件导入，这些对象是动态的，新导入的坐标集可以更新现有的可视对象或用于生成新的对象，从而允许用户在一个对象中聚合不同的轨迹。轨迹导航可以通过轨迹面板或方向键实现，多个对象可以同时显示或独立切换，为用户提供了比较模拟前后结构变化的灵活性。交互式3D查看器设计用于高效渲染包含数万个原子的系统，确保软件在处理复杂分子系统时具有良好的性能表现。鼠标交互方式借鉴了PyMOL和Coot等工具的常用操作模式，为用户提供了熟悉的操作体验。此外，EasyHybrid还提供了基于滚动的缩放控制，这对于在大型分子系统中导航至关重要，并且在3D画布中支持上下文相关的右键菜单，为用户提供了访问高级选项的便捷途径，如定义量子区域以进行QC计算或调整特定部分的可视化表示。
EasyHybrid提供了两种不同的原子选择机制：视图选择和拾取选择，均通过鼠标点击感兴趣原子的方式来实现。选择机制遵循“点击选中”/“点击取消选中”的逻辑，即同一操作既用于选中原子，也用于取消选中。视图选择以青色圆点表示（用户可以根据需要更改颜色），选中原子的顺序并不重要，且可以选中任意数量的原子。这种选择方式适用于多种任务，包括冻结原子位置、更改原子颜色、定义量子区域、剪枝或删除系统的一部分等。选择是动态的，支持加法和减法修改，允许用户逐步更新选择。EasyHybrid还允许用户在后续模拟步骤中保存和重用这些选择。视图选择模式下的选择标准包括：
- 仅选择点击的单个原子；
- 选择点击原子所属的整个残基；
- 选择与点击原子同一条链的所有原子；
- 仅选择点击的Cα原子；
- 选择包含点击原子的整个分子；
- 选择与溶剂分子相关的原子（例如水分子）。
相比之下，拾取选择模式遵循不同的逻辑。在这种模式下，最多可以选择四个原子，并且选择的顺序是严格保留的。该模式特别适用于需要原子顺序的操作，如计算距离、键角和二面角，或定义反应坐标。在3D可视化中，被选中的原子以彩色球体表示，并带有标签指示其选择顺序，如图S1所示。
为了增强用户体验，EasyHybrid针对包含大量原子的系统，尤其是生物分子系统，提供了多种分子结构表示方式。目前支持的表示类型包括：线框表示、棍棒表示、动态棍棒表示、原子球表示、范德华球表示、丝带或Cα轨迹表示，以及非键合原子的线框表示（如图S2所示）。一旦通过之前的选中操作激活了某种表示方式，它将应用于与该可视对象相关的所有轨迹帧。表示方式采用加法逻辑，用户可以在连续的事件中添加多个表示段，或根据需要移除它们。颜色分配同样遵循加法逻辑，碳原子默认使用与系统相关的参考颜色，但用户可以自由选择模型中的特定部分并为其分配不同的颜色。由于许多分子轨迹涉及原子转移（即键的形成和断裂），EasyHybrid还提供了动态键表示。这种表示方式在轨迹帧之间动态更新，与棍棒表示相似。与其他表示方式不同的是，动态键表示要求用户一次性提供所有需要渲染的原子选择。默认情况下，量子区域内的键使用动态键表示。在pDynamo3中，键的形成和断裂判断标准基于原子共价半径之和乘以一个校正因子，该因子可以调整。此外，EasyHybrid还支持分子轨道、电子密度和静电势面的表面表示。这些表面可以与轨迹链接，允许用户可视化诸如反应坐标扫描过程中前线轨道的变化等信息。图S3展示了在色氨酸合成酶反应坐标扫描中的HOMO表示。
EasyHybrid能够读取和导出pDynamo3的序列化文件，格式包括.pkl和.yaml。这为用户在图形界面之外进行模拟设置和执行提供了灵活性。这些文件包含了系统的全部信息，包括坐标和QC/MM参数。加载后，EasyHybrid会将MM原子以线条形式显示，QC原子以动态的棍棒形式呈现，固定原子则以灰色显示，便于系统检查。坐标还可以导出为计算化学中常用的格式，如.xyz、.pdb、.mol和.mol2。整个会话可以通过EasyHybrid的原生项目文件格式.easy保存，该格式将所有加载的系统和配置设置整合到一个文件中。最直接的加载系统方法是通过坐标文件，但要进行任何模拟，系统必须与一个能量模型相关联。最简单的模型是纯量子化学模型，因为坐标通常已足够。然而，QC-only模拟仅适用于小系统，因为其计算成本较高。
EasyHybrid提供了一个专门的窗口用于设置QC计算（见补充图）。用户可以选择使用pDynamo3原生方法或结合外部软件（如ORCA、xTB和DFTB+）进行计算，所有这些选项都可以通过专用的界面面板访问。每种选项都包含一个辅助窗口，用于设置所需的参数。将系统与分子力学（MM）模型关联则更为复杂，因为它需要拓扑信息以及原子类型和坐标。MM系统可以使用pDynamo3原生支持的力场（如OPLS、CHARMM、DYFF或pDynamo3的通用力场）构建。在这种情况下，用户必须提供一个包含拓扑信息的文件（如.mol2）和一个兼容的参数集。界面会建议默认的参数文件，但用户也可以替换为自己的文件。
另一种方法是使用外部软件（如Amber或CHARMM）准备的MM系统。这些系统可以导入并操作于EasyHybrid中。许多操作依赖于视图选择，例如原子固定和系统剪枝，这两种操作都能显著降低计算成本。原子固定功能允许用户在模拟过程中冻结选定的原子（如几何优化或分子动力学），防止其移动。系统剪枝功能则允许用户基于选定的原子创建一个更小的系统。原子选择的第三种重要用途是在混合QC/MM系统中定义量子区域。在这种情况下，用户需要从一个MM定义的系统开始，然后选择需要作为QC区域的原子。EasyHybrid会自动调整部分包含在QC区域中的残基的剩余MM部分的电荷，以确保总MM电荷为整数，这是大多数QC/MM方法所必需的。该程序还会存储原始电荷，以便在定义新的量子区域时，首先恢复原始电荷，从而最小化可能的误差积累。图2展示了EasyHybrid中定义QC区域的步骤。
EasyHybrid支持多种模拟类型，涵盖了计算化学协议的广泛范围，适用于量子力学（QC）、分子力学（MM）以及混合QC/MM势能模型。其支持的模拟类型包括：
- 单点能量计算：可以使用MM、QC或混合QC/MM势能模型进行能量评估。这些计算对于基准测试、比较系统构型或为后续模拟准备结构非常有用；
- 几何优化：支持使用pDynamo3库中的最速下降法和共轭梯度算法进行结构最小化。用户可以指定优化周期数、收敛标准以及是否在优化过程中保存中间结构的轨迹；
- 分子动力学（MD）模拟：EasyHybrid支持使用三种积分器（速度Verlet、leapfrog和随机Langevin动力学）进行经典MD模拟。界面允许用户定义标准的模拟参数，如时间步长、总模拟时间、温度、压力、轨迹保存频率和随机种子；
- 受约束弹性带（NEB）反应路径计算：NEB用于确定两个端点结构（通常是反应物和产物）之间的最低能量反应路径。界面允许用户定义初始和最终几何构型、中间图像（或副本）的数量以及其他NEB相关参数。这些图像随后会被优化，以在势能面上找到最低能量路径；
- 势能面（PES）扫描：可以通过应用几何约束（通常为距离限制或其线性组合）进行一维（1D）和二维（2D）PES扫描。在每个扫描点上进行几何优化。用户可以定义坐标范围、每个维度的扫描点数以及与约束相关的力常数；
- 通过伞形采样进行势能平均力（PMF）模拟：PMF计算提供了一个自由能曲线，作为所选反应坐标的函数。EasyHybrid中的实现与PES扫描相似，但使用短MD模拟代替几何优化。每个窗口的反应坐标轨迹可以通过pDynamo3中实现的加权直方图分析方法（WHAM）进行后处理，以重建整体自由能面。
所有模拟类型均通过pDynamo3的后端执行，并受益于EasyHybrid的集成可视化、选择和配置工具。对于QC和QC/MM模拟，用户可以选择使用pDynamo3原生方法或与外部引擎（如ORCA、xTB和DFTB+）结合的方法，所有这些方法均可通过专用的界面面板访问。
在EasyHybrid中，由pDynamo3库执行的模拟会生成多种格式的结果。所有pDynamo3过程均设计为输出一个日志文件，其中包含特定模拟的结果。EasyHybrid能够自动读取和解释这些日志文件，并以图形形式显示相关数据。用户可以保存和操作这些图表，从而方便地生成图形和结构表示。日志文件的处理在任何pDynamo3通过EasyHybrid执行的模拟结束后会自动触发，但也可以手动执行，以处理之前生成的EasyHybrid/pDynamo3日志文件。绘图功能由一个名为EasyPlot的自建工具处理，该工具使用Pycairo图形库开发。图3展示了通过EasyPlot对两个反应坐标同时进行的势能面扫描的可视化分析。另一类由pDynamo3生成的输出文件是轨迹文件。这些轨迹文件可能采用原生格式（如pkl）或外部格式（如CRD、NetCDF和DCD），并可能包含原子坐标、能量、反应坐标值、速度等数据。EasyHybrid支持多种pDynamo3轨迹类型，允许用户同时加载多个轨迹，并指定要处理的数据对象。界面还包含一组结构分析工具，包括在轨迹过程中监控多个距离、角度或二面角，以及进行如均方根偏差（RMSD）、结构对齐、重映射等计算。
EasyHybrid是GTKDynamo项目的重大演进，GTKDynamo最初是为了支持pDynamo版本1.7到1.9（基于Python 2）而设计的。在GTKDynamo中，用户界面结合了GTK2窗口和PyMOL 1.8作为3D查看器，以及Matplotlib用于绘图。相比之下，EasyHybrid已更新为使用现代版本的Python和GTK。值得注意的是，PyMOL已被一个基于现代OpenGL/GLSL开发的自定义3D渲染引擎取代，而Matplotlib则被内部开发的EasyPlot模块完全取代，用于图形可视化。EasyHybrid的用户界面设计受到多个领先的分子可视化程序的影响，包括PyMOL、VMD、Avogadro、wXMacMolPlt和Gabedit。作者们承认这些项目对社区的显著贡献，并对其长期影响表示感谢。然而，为了实现本文所述的功能，EasyHybrid是独一无二的，它专门为管理基于pDynamo3的QC/MM系统提供了优化的用户体验。表S1展示了EasyHybrid的一些功能与该领域其他知名免费软件的对比。
总之，本文介绍了EasyHybrid，一款现代、免费且开源的图形环境，专门用于支持量子/经典混合模拟。EasyHybrid旨在与pDynamo3模拟库相辅相成，为分子系统准备、模拟控制、轨迹分析和高级可视化提供了一个统一的平台。通过将Python脚本的灵活性与直观、用户友好的界面相结合，EasyHybrid为刚开始探索计算化学和混合方法的学生，以及寻求加速复杂模拟准备和分析的资深研究人员提供了便捷的途径。其处理大规模系统、直观定义量子区域以及支持NEB和伞形采样等高级协议的能力，使其成为计算化学家和结构生物学家的重要工具。EasyHybrid在Linux系统上免费提供，采用GNU通用公共许可证，并可通过Ubuntu终端（WSL）在Windows上运行，或通过虚拟机在Windows和macOS上运行。更多详细信息可在项目主页上找到。
