通过分子建模方法探索半胱天冬酶突变和翻译后修饰

本协议使用生物分子模拟包，并描述了用于模拟野生型半胱天冬酶及其突变形式的分子动力学（MD）方法。MD方法允许评估半胱天冬酶结构的动态演变以及突变或翻译后修饰的潜在影响。

在使用分子建模的启动和执行中发挥关键作用。我们研究了翻译后修饰和突变对半胱天冬酶结构和功能的影响。我们描述了分子动力学方法，该方法在原子水平上引入结构修饰后，给出了蛋白质进化的观点。

这种方法对于理解调节程序细胞死亡和相关疾病的机制以及开发新的有效疗法具有特殊意义。我们使用最流行的工具之一 Amber 20 演示分子建模功能。但是所提出的协议也可以与该软件的正式版本一起使用。

以下视频描述了如何制备半胱天冬酶结构以进行模拟以及对这种蛋白质野生型和修饰形式进行纤毛研究的分步说明。要检索选定的蛋白质数据库或P D B结构，请使用下载文件下拉列表并单击P D B格式。删除备注和连接数据，并在P D B文件中的单独蛋白质链之间插入T E R卡。

要准备起始模型，请从琥珀工具包启动T-LEAP程序。然后加载F F 14 S B力场，通过在命令行中输入指示的命令，用分子力学以及水分子和原子离子（如钠和氯化物）的参数来描述蛋白质。接下来，加载 P D B 文件并为氢构建坐标，从而创建一个名为 mol 的对象。

使用 check mol 命令检查可能导致问题的内部不一致。在蛋白质周围创建一个溶剂盒。然后通过输入电荷摩尔来检查总电荷并添加反离子以中和系统。

通过输入指示的命令创建拓扑 P R M 顶部文件和坐标 I N P C R D 文件。完成后，退出T-lEAP计划。进行能量最小化的第一阶段，以优化添加的氢原子和水分子的位置，同时通过重原子的位置约束保持蛋白质坐标固定。

通过输入指示的命令运行 P M E M D 程序。在我控制数据时遵循所需的参数。P分子拓扑结构、力场参数和原子名称。

C 初始坐标。O 用户可读日志输出 R 最终坐标。REF，位置约束的参考坐标。

接下来，无限制地执行能量最小化的第二阶段，以使用指示命令的输入优化整个系统。该阶段旨在将系统从零加热到 300 开尔文。通过使用给定的命令作为输入，在恒定体积下以50皮秒的速度对蛋白质原子进行位置约束的加热过程。

遵循所需的参数，因为 X 坐标集保存在分子动力学轨迹上。下一阶段是调整水的密度并获得蛋白质的平衡状态所必需的。在成功达到平衡后，使用指示的命令在恒定压力下以 300 开尔文的速度执行 500 皮秒的平衡，没有任何约束。

在恒定压力下进行 10 纳秒或更长时间的生产分子动力学模拟，并生成轨迹文件，以便使用指示的命令对蛋白质结构进行后续分析。该程序的并行版本，P M E M D M P I 或 G P U 加速版 P M E M D cuda 可用于计算机集群和超级计算机。长分子动力学模拟可以分成几个部分并按顺序执行。

本分析采用分子动力学建模工作流程研究了野生型半胱天冬酶二号及其丝氨酸-384丙氨酸突变体。丝氨酸-384丙氨酸取代诱导活性位点精氨酸-378的重要构象变化。此外，还表明丝氨酸-384丙氨酸取代影响了活性位点中精氨酸残基的底物识别，从而损害了铸基活性。

该位点定向诱变和生化测试表明，Serine-384丙氨酸突变阻断了半胱天冬酶二加工中的酶活性并抑制了癌细胞的凋亡。所描述的方法可用于评估半胱天冬酶二以及涉及各种癌症疾病的其他半胱天冬酶中的氨基酸突变和翻译后修饰的影响。