生物分子建模中的一项关键技能是在蛋白质中显示和注释活性位点。本技术演示使用了四种流行的免费大分子可视化程序:iCn3D、Jmol、PyMOL 和 UCSF ChimeraX。
生物分子建模中的一项关键技能是在蛋白质中显示和注释活性位点。本技术演示使用了四种流行的免费大分子可视化程序:iCn3D、Jmol、PyMOL 和 UCSF ChimeraX。
生物分子可视化技能对于理解生命科学中的关键概念至关重要,例如结构与功能的关系以及分子间的相互作用。多种程序允许学习者操作三维结构,生物分子建模能够促进主动学习、培养计算能力,并弥合二维教科书图像与生命真实三维结构之间的差距。该领域的一项关键技能是构建蛋白质活性位点模型,以展示大分子中能够与小分子(即配体)发生相互作用的部分,并呈现其结合相互作用方式。在本实验方案中,我们使用四种免费的大分子建模程序来描述这一过程:iCn3D、Jmol/JSmol、PyMOL 和 UCSF ChimeraX。本指南适用于希望学习特定程序基础操作的学生,也适用于将生物分子建模纳入课程教学的教师。该方案使用户能够使用某一特定可视化程序构建活性位点模型,或尝试多种可用的免费程序。本方案选用的人类葡萄糖激酶作为建模对象,它是己糖激酶的一种同工酶,可催化糖酵解的第一步反应。该酶结合了其一种底物以及一种非反应性底物类似物,使用户能够分析催化复合物中的相互作用。
理解分子世界的表征对于成为生物分子科学领域的专家至关重要1,因为对这些图像的解读是理解生物功能的关键2。学习者通常首先通过教科书中以二维形式呈现的细胞膜、细胞器、大分子等图像来接触大分子,但生物学的真实情况是,这些结构本质上是三维的,理解其性质需要能够可视化并从三维模型中提取信息的方法。
因此,高年级分子生命科学课程中生物分子可视化素养的培养已受到关注,多篇文章报道了教授和评估可视化技能的重要性及困难1,3,4,5,6,7,8,9。针对这些文章的回应是,课堂干预措施的数量有所增加,通常是在单一机构的一个学期内实施,利用分子可视化程序和模型来攻克难点概念2,10,11,12,13,14,15。此外,研究人员还致力于表征学生如何使用生物分子可视化程序和/或模型来学习特定主题16,17,18,19。我们所在的 BioMolViz 团队已提出一个框架,将可视化素养中的总体主题细分为学习目标和具体目标,以指导此类教学干预20,21,并组织培训工作坊,指导教师使用该框架进行逆向设计评估工具,以测量可视化素养技能22。
所有这些工作的核心都依赖于一项关键技能:利用生物分子可视化软件来操控大分子结构的能力。这些工具在不同的平台上独立开发,因此在操作和使用方式上可能各具特点。这就需要针对不同软件提供特定的操作说明,同时选择一款用户熟悉的软件对于持续应用至关重要。
除了掌握三维结构操作的基本技能(如旋转、选择和修改模型)之外,一个主要目标是构建蛋白质活性位点的模型。该过程有助于学习者深入理解生物分子可视化框架(BioMolViz Framework)所提出的三个核心主题:分子相互作用、配体/修饰以及结构与功能之间的关系20,21。
用于生物分子可视化的四种常用程序包括:Jmol/JSmol23、iCn3D24、PyMOL25 和 UCSF Chimera26,27。我们建议初次使用 Chimera 的用户采用 UCSF ChimeraX,这是 Chimera 分子可视化程序的下一代版本,也是目前官方支持的程序版本。
在本实验方案中,我们演示如何使用这四种程序中的每一种,对人葡萄糖激酶的活性位点进行建模,其中包含一个结合的底物类似物复合物(PDB 编号:3FGU),并展示测量结果以说明特定的结合相互作用28。该模型代表了酶的催化复合物。为了捕捉催化前状态的活性位点,将一种不可水解的 ATP 类似物结合到葡萄糖激酶的活性位点。这种磷氨酰磷酸腺苷酯(ANP)在此位置含有一个磷-氮键,取代了通常的磷-氧连接键。活性位点中还包含葡萄糖(模型中记为 BCG)和镁离子(记为 MG)。此外,结构中存在一个钾离子(K),来源于结晶过程中使用的氯化钾溶剂。该离子对生物学功能不具关键作用,且位于活性位点之外。

图1:ATP/ANP 结构。 三磷酸腺苷(ATP)结构与磷氨酰基磷酸-腺苷酯(ANP)的比较。请点击此处查看该图的放大版本。
该方案展示了如何筛选底物类似物复合物中结合的配体,并鉴定与结合复合物距离在5 Å以内的活性位点残基,从而识别出能够形成相关分子相互作用(包括疏水作用和范德华作用)的氨基酸和水分子。
最初调整显示方式,以卡通形式展示蛋白质的大部分结构,同时将活性位点的氨基酸残基以棍状形式显示,以呈现蛋白质的相关原子并突出分子相互作用。在每个程序的操作步骤3完成后,这些显示方式已被应用,蛋白质的视图在不同程序间保持相似(图2)。在操作流程结束时,隐藏蛋白质的卡通表示,以简化视图,突出活性位点。

图2:不同程序间的结构比较。 在执行“调整显示方式”步骤后(各方案的步骤2或步骤3),对程序中3FGU结构的比较。请点击此处查看此图的放大版本。
活性位点的氨基酸和结合配体采用CPK着色法29,30。该着色方案可区分在线条、棍状、球棍和空间填充模型中显示的不同化学元素的原子。在CPK着色方案中,氢为白色,氮为蓝色,氧为红色,硫为黄色,磷为橙色。传统上,碳用黑色表示,但在现代应用中,碳的颜色可能有所不同。
在晶体结构中,氢原子不可见,尽管上述每种程序均能预测其位置。在大型大分子结构中添加氢原子可能会遮挡视图,因此本方案中不显示氢原子。因此,这些结构中的氢键将通过测量两个杂原子中心之间的距离(例如,氧到氧、氧到氮)来表示。
程序概述
可下载的图形用户界面(GUI):PyMOL(版本 2.4.1)、ChimeraX(版本 1.2.5)和 Jmol(版本 1.8.0_301)均为基于图形用户界面的分子建模工具。这三种界面均提供命令行以输入代码,同时许多相同的功能也可通过 GUI 中的菜单和按钮实现。这些程序命令行的一个共同特点是,用户可使用键盘上的上下箭头键加载并重新执行先前输入的命令。
基于网页的图形用户界面:iCn3D(I-see-in-3D)是一种基于 WebGL 的浏览器工具,用于在网页上交互式查看三维大分子结构和化学物质,无需安装独立应用程序。它不依赖命令行操作,但完整的网页版本提供可编辑的命令日志。JSmol 是 Jmol 的 JavaScript 或 HTML5 版本,可在网站或网页浏览器窗口中使用,其操作方式与 Jmol 非常相似。JSmol 可用于创建包含动画的在线教程。
Proteopedia31,32、Jmol 中的 FirstGlance33 以及位于密尔沃基工程学院生物分子建模中心的 JSmol 网络界面(JUDE)是基于 Jmol 的在线设计环境的典型示例34。Proteopedia 维基是一个教学工具,允许用户对大分子结构进行建模,并在网站内创建包含这些模型的页面35。该 Proteopedia 场景编写工具基于 JSmol 构建,集成了图形用户界面(GUI),并提供了 Jmol GUI 中所不具备的附加功能。
Jmol 和 iCn3D 基于 Java 编程语言;JSmol 使用 Java 或 HTML5,而 PyMOL 和 ChimeraX 基于 Python 编程语言。这些程序均可加载蛋白质数据库文件,此类文件可从 RCSB 蛋白质数据库根据一个四位字母数字组成的 PDB ID 下载36,37。最常见的文件类型是含有 .pdb 扩展名的蛋白质数据库(PDB)文件,以及含有 .cif 扩展名的晶体学信息文件(CIF 或 mmCIF)。CIF 已取代 PDB 成为蛋白质数据库的默认文件类型,但两种文件格式均可在这些程序中使用。使用 CIF 文件与 PDB 文件时,序列/结构的显示方式可能存在细微差异;然而,两种文件的功能相似,此处不再详细讨论这些差异。分子建模数据库(MMDB)是美国国家生物技术信息中心(NCBI)开发的产品,为蛋白质数据库结构的一个子集,其中关联了分类信息(例如生物学特征、保守的蛋白质结构域)38。iCn3D 是 NCBI 开发的工具,能够加载包含 MMDB 数据的 PDB 文件。
要查看模型,用户可从该结构对应的 Protein Data Bank 专用页面下载所需文件(例如, https://www.rcsb.org/structure/3FGU),然后使用程序中的下拉 文件菜单打开该结构。所有程序也可通过界面直接加载结构文件,该方法在实验方案中有详细说明。
ChimeraX、Jmol 和 PyMOL 的图形用户界面各自包含一个或多个控制台窗口,可通过拖动角落调整其大小。iCn3D 和 JSmol 完全运行于网络浏览器中。使用 iCn3D 时,由于屏幕尺寸和分辨率的不同,用户可能需要在弹出窗口内滚动,以查看全部菜单项。
本文详述的实验方案提供了一种利用各程序展示酶活性位点的简便方法。需要注意的是,每个程序中的操作步骤均可通过多种方式执行。例如,在ChimeraX中,同一任务可通过下拉菜单、顶部工具栏或命令行来完成。建议希望深入了解某一特定程序的用户查阅这些程序提供的在线教程、使用手册和维基文档39,40,41,42,43,44,45,46。
现有这些程序的手册和教程将本方案中的各项内容作为独立的任务进行介绍。为了显示一个活性位点,用户必须从各种手册和教程中自行整合所需的各项操作。本文通过提供一个用于建模带有分子相互作用的标记活性位点的线性操作流程,补充了现有的教程资源,为用户提供了可用于活性位点建模的逻辑思路,该思路可进一步应用于其他模型和程序中。

图 3:ChimeraX 图形用户界面。 ChimeraX 图形用户界面,标注了下拉菜单、工具栏、结构查看器和命令行。 请点击此处查看此图的放大版本。

图 4:iCn3D 图形用户界面。 iCn3D 图形用户界面,标注了下拉菜单、工具栏、结构查看器、命令日志、选择集弹出窗口以及序列和注释弹出菜单。请点击此处查看此图的放大版本。

图5:Jmol 图形用户界面。 Jmol 图形用户界面,标注了下拉菜单、工具栏、结构查看器、弹出菜单以及控制台/命令行。 请点击此处查看此图的放大版本。

图6:PyMOL 图形用户界面。 PyMOL 图形用户界面,标注了下拉菜单、结构查看窗口、名称/对象面板、鼠标控制菜单以及命令行。 请点击此处查看此图的放大版本。
注意:每个程序的流程均分为十个总体步骤:(1)将结构载入程序;(2)识别活性位点中的配体;(3)调整结构显示方式;(4)选择距活性位点5 Å范围内的残基以定义活性位点;(5)显示酶与活性位点配体之间的相互作用;(6)以棍状形式显示侧链,并显示/调整活性位点中的水分子;(7)简化结构;(8)标记配体及形成氢键的侧链;(9)在任意阶段保存渲染结果,以便后续继续编辑或与他人共享;(10)保存用于嵌入或打印的图像。其中,步骤1、4以及7至10在各程序中均相同;然而,由于各程序操作方式不同,部分流程在交换步骤2/3与5/6的顺序后执行效率更高。
1. UCSF ChimeraX 操作流程
注意:触控板与鼠标控制。旋转操作:点击并拖动,或使用双指拖动(鼠标操作:左键点击并拖动)。缩放操作:捏合或张开手指(Mac)或按住 Ctrl 键并使用双指移动(PC)(鼠标操作:滚动滚轮)。平移操作(即移动整个结构):按住 Option 键并点击拖动(Mac),或按住 Shift 键并点击拖动(PC)(鼠标操作:右键点击并拖动)。重新居中:使用界面顶部的下拉菜单,点击 操作 > 视图。
2. iCn3D 操作流程
注意:触控板和鼠标控制:旋转操作,请点击并拖动(鼠标:左键点击并拖动)。缩放操作,请捏合或张开手指(鼠标:滚动滚轮)。平移(即移动整个结构),请用两个手指点击并拖动(鼠标:右键点击并拖动)。重新居中,请将光标悬停在顶部下拉菜单中的“View”上,然后点击 Center Selection。
3. Jmol 实验方案
注意:触控板和鼠标操作说明:旋转操作为点击并拖动(鼠标操作:左键点击并拖动)。缩放操作为使用两个手指垂直滚动(鼠标操作:按住 Shift 键 + 左键点击并垂直拖动)。平移操作(即移动整个结构)为 Control + Alt + 点击并拖动(PC),Control + Option + 点击并拖动(Mac)。重新居中操作:在结构查看器窗口的空白区域按下 Shift 键并双击。
4. PyMOL 操作流程
注意:触控板和鼠标操作说明:旋转时,请点击并拖动(鼠标操作:左键点击并拖动)。缩放时,请用双指捏合或张开(鼠标操作:右键点击并拖动)。平移(即移动整个结构)时,请按住 Control 键 + 点击并拖动(鼠标操作:按住 Command 键 + 左键点击并拖动)。如需重新居中,请前往右侧对象面板,点击 A > 定向或居中。
每项程序成功执行后,将得到一个聚焦于活性位点的分子模型,其中活性位点的残基和配体以棍状形式显示,蛋白的卡通表示被隐藏,配体采用对比色方案呈现。相互作用的氨基酸残基应标示其残基编号,并用线条显示氢键和离子相互作用。这些特征的存在可通过模型的目视检查来确定。
为了便于进行此项检查,并帮助用户判断是否正确执行了实验方案的各个步骤,我们提供了动画图示,以展示每一步骤完成后结构的形态。对于 ChimeraX、iCn3D、Jmol 和 PyMOL,分别如图7-10所示。
图7:ChimeraX 实验方案输出结果。 该动画图示展示了 ChimeraX 实验方案的步骤1.1–1.8。请点击此处下载该图。
图8:iCn3D 协议输出结果。 动态图示展示 iCn3D 协议的步骤 2.1–2.8。请点击此处下载该图。
图9:Jmol 实验方案输出结果。 动态图示展示 Jmol 实验方案的步骤 3.1–3.8。请点击此处下载该图。
图10:PyMOL 实验方案输出结果。 动态图示展示 PyMOL 实验方案的步骤 4.1–4.8。请点击此处下载该图。
影响这些操作结果最常见的错误是选择不当,导致结构的一部分以不期望的渲染方式显示。这通常是由于在结构本身或显示菜单按钮上的误点击所致。一个次优结果的示例是,模型中包含活性位点之外的残基以棍状形式显示。用户可通过目视检查以棍状显示的残基,确认它们是否位于活性位点配体附近,从而初步判断是否发生了此类错误。评估显示的残基是否在活性位点配体5Å范围内的进阶方法是,使用各程序内置的测量工具,测量邻近配体与活性位点残基之间的距离。测量工具的使用超出了本文的范围;然而,我们鼓励感兴趣的用户查阅大量在线教程,深入了解此类分析的具体操作。
我们展示了一个由于在 PyMOL 中错误点击名称/对象面板而导致该实验方案执行不理想的特定示例。此错误会以棍状模型显示整个蛋白质,而不是仅用该表示方式显示活性位点,如图11所示。

图11:阴性结果。 阴性结果的示例。在 PyMOL 中错误地选中了完整的卡通图并显示为棍状模型。请点击此处查看该图的放大版本。
为排除故障,用户需要先隐藏整个模型的棍状表示(在名称/对象面板中标记为 3FGU),然后仅对名为“active”的选定部分使用 PyMOL 中的隐藏和显示按钮/命令来显示其棍状表示。一旦用户能够为模型的不同部分创建合适的选定区域,并有效地显示或隐藏它们,从此类错误中恢复模型就相对简单了。尽管重新开始实验流程并再次逐步操作看似诱人,但我们鼓励用户不要拘泥于既定步骤,大胆对模型进行探索和尝试。根据我们的经验,通过解决显示错误的过程,有助于加深对建模软件的理解和掌握。
每个程序成功执行方案后的最终输出结果并列显示在图12中。视图采用相似的取向,以便用户比较不同程序所创建模型的外观。

图 12:不同程序的最终结构对比。 对本方案结束时各程序生成的活性位点结构进行比较。A:ChimeraX,B:iCn3D,C:Jmol,D:PyMOL。PyMOL 的活性位点标签包含所有活性位点残基及配体,而其他程序的输出仅标注了形成氢键的侧链。请点击此处查看该图的高清版本。
本方案概述了酶活性位点建模的十个步骤,并将其应用于四种常用的生物分子建模软件。该方案的关键步骤包括:识别活性位点中的配体,选择距活性位点5 Å以内的残基以界定活性位点区域,以及展示酶与活性位点配体之间的相互作用。区分与生物学功能相关的配体至关重要,因为这使用户能够确定5 Å范围内可能参与配体结合的氨基酸残基。最后,利用程序展示分子相互作用,使用户能够掌握理解促进结合的分子相互作用所必需的技能。
基于计算机的分子建模方法的一个局限性在于其对特定命令和语法的依赖性。尽管生物化学实验方案可能允许操作步骤中有微小变动,但基于计算机的研究若未严格遵循既定流程,则可能产生截然不同的最终结果。这一点在使用命令行界面时尤为重要,因为必须使用程序特定的语法才能获得预期输出,而标点符号或大小写的看似微不足道的更改都可能导致命令执行失败。每种程序通常都有相应的各种 Wiki 页面和用户手册,用户可在其中查找并排查命令行输入的问题;因此,用户应仔细关注命令语法的每一个细节。尽管大多数分子可视化程序都包含“撤销”命令,但由于界面复杂,“撤销”功能并不总能准确地逆转上一步操作。因此,建议经常保存当前工作状态,特别是对于新用户而言更为重要。
模型本身所使用的数据也可能带来进一步的局限性。尽管蛋白质数据库(Protein Data Bank)中的固有标准确保了一定程度的一致性,但分子可视化程序的使用者常常会在蛋白质的渲染图像中遇到意想不到的效果。首先,大多数结构是通过X射线晶体学确定的,该方法仅提供蛋白质的一个单一模型;然而,核磁共振(NMR)结构通常由多个模型组成,可逐个进行可视化。其次,通过晶体学或冷冻电子显微镜实验确定的结构可能包含某些无法明确其位置的原子,这些原子在蛋白质的某些表示形式中会呈现为缺口。蛋白质结构中的侧链可能存在多种构象,在以棍状模式显示时,会表现为从同一氨基酸主链上伸出的两组不同结构。甚至短段的主链也可能存在此类替代构象,有时配体在活性位点中也会以多种结合构象叠加存在。
对于晶体结构,已提交的三维坐标包含非对称单元的所有组分,这些信息足以重构蛋白质晶体的重复单元。有时,该结构所包含的蛋白质链可能多于蛋白质的生物学活性形式(例如,胎儿血红蛋白突变体,PDB ID: 4MQK)。相反,某些程序可能不会自动加载生物学活性单元中的全部蛋白链。例如,在ChimeraX、PyMOL和Jmol中使用本实验方案所述命令获取SARS-CoV-2主蛋白酶(PDB ID: 6Y2E)结构时,仅加载了生物学活性二聚体的一半(由两条蛋白链组成)。尽管对命令稍作修改即可加载完整的生物学活性二聚体,但对于建模软件的新手用户而言,这一操作可能并不直观。另一个可能出现的问题是活性位点或底物本身的识别。晶体学实验中使用了多种分子,这些分子可能被建模到最终结构中。例如,硫酸根分子可能结合在活性位点的磷酸结合位点,也可能结合在与催化机制无关的其他区域。这些分子可能妨碍对活性位点的正确识别,甚至可能使学生误认为它们参与了催化机制。
用户可能希望将此方法应用于其他活性位点或结合位点。在今后涉及新蛋白质活性位点分析的工作中,若要应用本实验方案,用户需要确定哪些结合的配体与蛋白质功能相关。某些配体并不参与蛋白质功能,而是实验所用溶剂或结晶条件的产物(例如 3FGU 模型中存在的钾离子)。关键配体应通过查阅原始文献来确认。通过实践,以及在适用的情况下掌握命令行语法,用户将能够将该方案应用于任意酶活性位点的建模,并对其选择的其他大分子进行建模。
鉴定和分析结合的底物与配体对于阐明分子机制以及基于结构的药物设计至关重要,这些研究已直接推动了包括获得性免疫缺陷综合征(艾滋病)和 COVID-19 在内的疾病治疗手段的进步47,48,49,50,51,52。尽管不同的分子可视化程序提供各具特色的界面和用户体验,但大多数程序具备相似的功能。在生物化学高年级学生中培养生物分子可视化能力至关重要,他们需要熟悉结构可视化及其图像生成工具4,20,53。这有助于学生超越教科书和期刊文章中二维图像的解读,更轻松地基于结构数据提出自己的科学假设54,从而为未来科研人员应对公共卫生问题、深化对生物化学过程的理解做好准备。
综上所述,本实验方案详细介绍了使用四种主流的免费大分子建模程序进行活性位点建模的方法。我们的社区 BioMolViz 采用与具体软件无关的生物分子建模教学方法。本文并未对各程序的功能进行批评或比较,尽管用户在试用不同程序后,可能会发现其中某一程序在某些大分子建模方面更符合个人偏好。我们邀请读者使用 BioMolViz 教学框架,该框架明确了本方案所涉及的基于生物分子可视化学习的学习目标与教学目的,并通过 BioMolViz 社区网站 http://biomolviz.org 获取有关生物分子可视化教学与学习的更多资源。
作者声明,他们对本文所述研究不存在任何相关或重大的经济利益。
本工作的经费由美国国家科学基金会提供:
提升本科 STEM 教育资助项目(资助编号 #1712268)
本科生物学教育中的研究协调网络(资助编号 # 1920270)
我们感谢西菲尔德大学的 Karsten Theis 博士就 Jmol 问题进行的有益讨论。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| ChimeraX(版本 1.2.5)https://www.rbvi.ucsf.edu/chimerax/ | |||
| 计算机 | 任意 | ||
| iCn3D(仅限网络版:https://www.ncbi.nlm.nih.gov/Structure/icn3d/full.html) | |||
| Java(用于 Jmol)https://java.com/en/download/ | |||
| Jmol(版本 1.8.0_301)http://jmol.sourceforge.net/ | |||
| 鼠标(可选) | 任意 | ||
| PyMOL(版本 2.4.1 - 教育版):https://pymol.org/2 仅限教育用途版本:https://pymol.org/edu/?q=educational |