Visual Dynamics 3.0 中的新增功能

摘要

Visual Dynamics 是一个开源工具，可使用 Gromacs 加速分子动力学仿真的实现和学习。所提出的协议将指导您轻松完成在 ACPYPE 中制备的蛋白质-配体模拟的步骤以及其他模拟模型的一般步骤。

摘要

视觉动力学 （VD） 是一种 Web 工具，旨在促进在 Gromacs 中执行的分子动力学 （MD） 的使用和应用，允许不熟悉计算的用户运行短时模拟以进行验证、演示和教学目的。确实，量子方法是最准确的。但是，目前没有计算可行性来执行 MD 执行的实验。这里描述的工具在过去几年中不断改进。该协议将描述使用先前在 ACPYPE 中制备的蛋白质-配体复合物在 VD 中运行模拟所需的内容，以及其他可用模拟模型的一些一般说明。对于详细模拟，将使用来自间日疟原虫与抑制剂 D5 （PDB ID： 4mgv） 复合的 FK506 结合蛋白，并提供使用的所有文件。请注意，此协议将告诉使用每个选项来实现所呈现的相同结果，但这些选项不一定是唯一可用的选项。

引言

根据 IUPAC 的定义，MD 是一种模拟程序，包括根据牛顿运动定律计算分子中原子的运动或固体、液体和气体中单个原子或分子的运动。作用在原子上的力是模拟原子运动所必需的，通常使用分子力学¹ 中的力场来计算。它可以应用于任何寻求在分子水平（通常是原子水平）上提取信息的现象²。

MD 是纳入生物信息学的技术之一，特别是结构生物信息学。有了它，就可以获得生物分子结构的动力学和热力学特性。例如，大分子稳定性、变构位点的鉴定、酶活性机制的阐明、分子识别和与小分子复合物的性质、蛋白质之间的结合、蛋白质折叠及其水合作用³。此外，MD 支持广泛的研究，包括分子设计（广泛用于药物设计）、确定结构及其细化（X 射线、NMR 和蛋白质建模）³。就非量子模拟而言，在 MD 结束时获得的结果是最丰富和最完整的⁴。由于大量近似的数量，经典 MD 的效率比充分考虑生物分子系统物理学的预期要有效得多。值得注意的是，量子动力学效应通常被忽略³。然而，实施 MD 实验并非易事⁵.它需要计算知识，尤其是 Linux 终端，因为大多数结构生物信息学软件都是为此而设计的。即使有了这些知识，学习 Gromacs 命令和参数化也是另一个陡峭的学习曲线。

自 1977 年⁶ 首次应用于生物学以来，由于计算处理的增加和编码的改进，已经发生了很大变化。二十多年前，第一款用于生物问题的 MD 软件问世，即 Gromacs⁷、AMBER⁸ 和 NAMD⁹。

自第一个版本以来，这些软件仍然是使用和引用最多的软件。然而，它们继续面临同样的常见实现困难，这些困难困扰着非计算机专家的研究人员⁵ （computer specialists）。有些具有复杂的安装和配置步骤，有时需要有关它将运行的硬件的广泛知识才能充分利用它，以及高度以计算机为中心的技术文档。除了命令行和无限参数之外，还需要一种更简单的方法来与它们交互。

接口充当要执行的逻辑过程与人类¹⁰ 之间的中介。随着计算资源的改进，软件执行方式的范式也在发展。第一个数字范式是命令行界面 （CLI），随后演变为已知的图形用户界面 （GUI）¹¹。遵循进化周期，万维网（或简称 WEB）产生的界面被认为是 GUI¹¹ 的演变。这三种范式目前共存，具体取决于开发人员。CLI 应用程序在操作系统控制台上使用文本命令。GUI 应用程序（也称为图形桌面）使用由窗口、按钮和其他组件组成的图形界面。它是特定于操作系统的，并且已针对操作系统进行了预编程。与 CLI 的主要区别在于使用鼠标作为人机交互中的附加元素¹²。WEB 应用程序尽管与 GUI 混淆，但开发起来更复杂，但用途更广，并且是迄今为止最灵活的操作。此外，它们仅依赖于称为浏览器的解释软件，这使得客户端应用程序可以通过独立于操作系统的网络与服务器通信¹³。

结构生物信息学软件最常使用 CLI 和 GUI 范式。使用 CLI 的经典软件的一些示例包括用于相似性建模的 Modeller¹⁴ 、用于分子对接的 Autodock¹⁵ 和用于分子动力学的 Gromacs¹⁶ 。采用 GUI 类型的软件示例包括 SwissPDBviewer¹⁷、Pymol¹⁸、VMD¹⁹、UCSF Chimera²⁰、Autodock tools¹⁵、PyRx²¹、Biovia²²、Maestro²³ 和 Moe²⁴ 等。

随着超文本标记语言版本 5 （HTML5）²⁵、级联样式表 （CSS）²⁶ 和 Javascript²⁷ 技术等的出现，许多结构生物信息学应用程序可以被引入 WEB，从而变得更加易于访问。相似性建模 WEB 服务器的示例包括 MODWEB²⁸（使用 Modeller¹⁴ 作为后端）和 Swissmodel²⁹。用于分子对接的 Web 应用程序服务器的示例包括 Haddock³⁰、Swissdock³¹、Cluspro³²、Dockthor³³ 等。

虽然结构分析、建模和对接方法从 CLI 范式发展到 GUI，最后发展到 WEB，但 MD 仍然主要由命令行执行（CLI 类型）支持。已经出现了一些很好的举措来改善这种状况。这些举措的示例包括在现有软件中实施插件，例如 VMD³⁴ 的 QwikMD 插件、PyMOL 的 GROMACS 插件和 UCSF Chimera²⁰ 中的分子动力学模拟选项，一些新的和更简单的 CLI 应用程序，例如 ASGARD³⁵、Gmx_qk³⁶ 和 CHAPERONg³⁷，以及强大的 Web 平台 BioBB-Wfs³⁸.尽管使用这些插件和应用程序是一项进步，但对于大多数不熟练的研究人员来说，它们的实现仍然是一个挑战。常见的困难包括安装和配置 MD 软件时出现问题，这通常会影响模拟⁵ 的完全执行。

2022 年，位于 Fiocruz Rondônia³⁹ 的 Laboratório de Bioinformática e Química Medicinal 提供了用于基于 Web 的计算模拟的 Visual Dynamics 软件。它的初始版本是用 Python 和 Flask 构建的，允许仅用 2 ns 模拟具有游离蛋白（脱辅基酶）的系统。随后，它被增强为包括使用 PRODRG⁴⁰ 制备配体的自动模拟版本。

VD 旨在帮助结构生物物理学、生物技术和相关领域所有在计算知识方面有局限性的研究人员;该工具允许这些研究人员从任何操作系统测试他们涉及 MD 模拟的假设，而无需访问高性能计算机 （HPC）。这项工作的目的是介绍 Visual Dynamics 3.0 版的新功能。此外，它旨在为该工具引入更新的使用协议，并强调未来需要解决的限制，以及截至目前的使用统计数据（图 1）。

研究方案

1. 访问软件和新用户注册

1. 访问 Visual Dynamics （VD） 网页。单击右上角的 +注册 图标以创建一个帐户。注册以使用该软件。
   注意：仅允许使用机构电子邮件地址。用户注册获得批准后，将收到一封电子邮件通知。
2. 单击右上角的 Login 以访问系统登录屏幕。填写用户名/电子邮件和密码字段，然后单击 Log in。登录后，用户将可以访问其仿真提交区域。他们还可以查看 VD 的教程和使用情况统计信息。

2. 脱辅基酶模拟提交

1. 单击左侧边栏中的 New Simulation（新建模拟 ）。在出现的屏幕上，单击按钮 APO，它指的是脱辅基酶。
2. 上传免费的蛋白质 4mvg.pdb 文件。选择 AMBER94 力场 （或任何其他适当的选项）。
   注：4mvg.pdb 可从 补充文件 1 获得，也可直接从蛋白质数据库 （PDB） 下载。
3. 选择 TIP3P Water 模型。选择 Cubic Box。选择蛋白质和盒边缘之间的 0.5 nm 距离 。
   注意：选择的选项只是建议。所有其他选项在 VD 中都有效。
4. 选中选项 Run in Our Servers 以使用提交的文件和参数执行模拟。单击 Run simulation（运行仿真）后，用户将在屏幕上看到仿真步骤的演变，并在完成后收到有关仿真状态的电子邮件通知。

3. 提交在 ACPYPE 中制备的酶与配体复合的模拟

1. 使用 UCSF Chimera²⁰，打开 蛋白质-配体复合物 4mgv.pdb⁴¹，在“选择”下，单击 “残基”，然后将代码设置为 D5I。然后，在文件下，单击 保存 PDB，选择 仅保存所选原子，将文件名设置为 ligand.pdb，然后单击 保存。
2. 将上一步生成的文件 ligand.pdb 提交给 Bio2Byte ACPYPE 服务器⁴²，从输出文件中，ligand_NEW.itp 和 ligand_NEW.pdb 将是本实验中使用和提供的。
   注意：ligand_NEW.itp 和 ligand_NEW.pdb 文件可以从 补充文件 2 和 补充文件 3 中获取。
3. 单击左侧边栏上的 New Simulation（新建模拟 ）。点击 按钮蛋白质 + 配体 （在 ACPYPE 中制备）。
4. 上传免费的蛋白质 4mvg.pdb 文件。选择在 ACPYPE 中准备的 配体文件 ：ligand_NEW.itp 和 ligand_NEW.pdb。选择 AMBER94 力场。
5. 选择 TIP3P Water 模型。选择 Cubic Box。选择蛋白质和盒边缘之间的 0.5 nm 距离 。选择的选项只是建议。所有其他选项在 VD 中都有效。
6. 选中选项 Run in Our Servers 以使用提交的文件和参数执行模拟。单击 Run simulation（运行仿真）后，用户将在屏幕上看到仿真步骤的演变，并在完成后收到有关仿真状态的电子邮件通知。

4. 访问仿真结果

1. 单击左侧边栏中的 My Simulations（我的模拟 ）。
2. 单击 Download MDP Files（下载 MDP 文件 ），将平台使用的仿真配置文件下载到用户的计算机上。
3. 下载模拟元素，如下所述。
   1. 单击 Commands 将平台执行的命令列表下载到用户的计算机上。
   2. 单击 GROMACS Log 将包含 .gmx 命令输出的日志文件按顺序下载到用户的计算机上。
   3. 单击 “结果 ”可将 .gmx 命令生成的文件（如 _npt.gro、_pr.edr、_pr.tpr、_pr_PBC.gro 和 pr_PBC.xtc）下载到用户的计算机上。
   4. 单击 Figure Graphics 将用于分析每个仿真步骤的图像和 .xvg 格式的图形下载到用户的计算机上。

结果

VD 提供完全自主的仿真执行，不需要用户干预或用户提供的计算资源。将模拟提交到执行后，用户可以离开它，关闭他们的机器，模拟将继续运行。它还允许用户从任何设备访问结果，无论是笔记本电脑还是移动设备。

作为通过 WEB 在自动模式下使用 VD 的一个例子，使用间日疟原虫的 FK506 结合蛋白与抑制剂 D5 （PDB ID： 4mgv）⁴¹ 的结构，对在 ACPYPE 中?...

讨论

自动化流程并不容易，但也比从头开始重新编程系统要容易。Gromacs 是目前最流行的分子模拟软件，并且不断更新。它最初是由格罗宁根大学生物物理化学系开发的，现在由斯德哥尔摩大学的生命科学实验室维护⁴³.

对于任何新用户来说，学习模拟技术都是一个漫长的旅程。VD 作为一种替代方案出现，以促进此学习过程，并在投入资源运行更长时间的模拟之前...

材料

Name	Company	Catalog Number	Comments
ACPYPE Server	Bio2Byte		Available at https://www.bio2byte.be/acpype/
GRACE software	Plasma Laboratory at the Weizmann Institute of Science		Available at https://plasma-gate.weizmann.ac.il/Grace/
GROMACS software	GROMACS Team		Installation instructions at https://manual.gromacs.org/current/install-guide/index.html
The structure of the FK506-binding protein From Plasmodium vivax complexed with the inhibitor D5	RCSB Protein Data Bank		Available at https://www.rcsb.org/structure/4mgv Already contains the ligand complexed to the macromolecule.