Name: 通过分子建模方法探索半胱天冬酶突变和翻译后修饰 (视频) | JoVE Journal
Uploaded: 2022-10-13T14:00:00.000+00:00
Duration: 5 min 56 s
Description: 探索分子动力学模拟方法研究半胱天冬酶突变和翻译后修饰对蛋白质结构与功能的影响。该方法使用Amber 20工具，可评估氨基酸突变在癌症相关半胱天冬酶中的作用机制。

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

本研究得到了俄罗斯科学基金会（项目编号：17-75-20102，用于实验方案开发）的资助。代表性结果部分所述实验（磷酸化分析）得到了斯德哥尔摩癌症协会（181301）和瑞典癌症协会（190345）的支持。

1. 系统准备
注意：天然蛋白和突变体蛋白的分子模型是基于从蛋白质数据库（Protein Data Bank）获得的合适晶体结构构建的15,16。
<ol>
	<li>为获取选定的 PDB 结构，使用 Download Files 下拉菜单并点击 PDB Format。从 PDB 文件中移除注释和连接性数据，并在不同蛋白质链之间插入 TER 卡片。可选地，通过将组氨酸残基名称 HIS 替换为 HIE、HID 或 HIP 来设定其质子化状态（补充文件 1）。 
	注意：保留晶体结构中解析出的水分子（如果 PDB 文件中存在）是合理的。</li>
	<li>为准备起始模型，启动 tleap 程序（AmberTools 软件包），然后输入用于操作 tleap 对象的命令。
	<ol>
		<li>启动程序：tleap。加载 ff14SB 力场以用分子力学方法描述蛋白质：&#62; source leaprc.protein.ff14SB。</li>
		<li>加载水分子及原子离子（Na+、Cl-）的力场：&#62; source leaprc.water.tip3p。</li>
		<li>加载 PDB 文件 并构建氢原子坐标，创建一个名为 mol 的对象：&#62; mol = loadpdb protein.pdb。</li>
		<li>检查可能引起问题的内部不一致性：&#62; check mol。 
		注意：如果存在二硫键，必须手动指定。将共价连接的半胱氨酸残基名称 CYS 替换为 CYX，并在 tleap 中输入以下命令以在 SG 原子间建立化学键：bond mol.X.SG mol.Y.SG（X 和 Y 为半胱氨酸残基编号）。</li>
		<li>在蛋白质周围创建溶剂盒（TIP3P 水模型，12 Å 距离）：&#62; solvatebox mol TIP3PBOX 12.0。</li>
		<li>检查总电荷：&#62; charge mol，并添加抗衡离子（Na+ 或 Cl-）以中和体系： 
		&#62; addions mol Na+ 0。</li>
		<li>生成拓扑 prmtop 文件 和坐标 inpcrd 文件（用于运行模拟的输入文件）：&#62; saveamberparm mol protein.prmtop protein.inpcrd。退出 tleap 程序：&#62; quit。 
		&#8203;注意：可使用 ambpdb 工具轻松将 Amber 格式的坐标文件转换为 PDB 格式（参见《用户手册》，见 材料表）。</li>
	</ol>
	</li>
</ol>
2. 能量最小化
注意：能量最小化是必要的，可消除初始体系中原子间的不良接触和重叠，这些因素在进行分子动力学模拟时会导致体系不稳定。
<ol>
	<li>执行能量最小化的第一阶段（最陡下降法 2,500 步 + 共轭梯度法 2,500 步），在对重原子施加位置限制以固定蛋白质坐标的同时，优化所添加的氢原子和水分子的位置。
	<ol>
		<li>按如下方式运行 pmemd 程序： 
		pmemd -O -i min1.in -p protein.prmtop -c protein.inpcrd -o protein_min1.out 
		-r protein_min1.rst -ref protein.inpcrd。</li>
		<li>遵循所需的参数：-i FILE 控制数据；-p FILE 分子拓扑、力场参数、原子名称；-c FILE 初始坐标；-o FILE 用户可读的日志输出；-r FILE 最终坐标；-ref FILE 用于位置限制的参考坐标。</li>
		<li>在输入文件 min1.in 中指定选项： 
		&amp;#38;cntrl 
		imin=1, maxcyc=5000, ncyc=2500, 
		cut=10.0, ntb=1, 
		ntc=1, ntf=1, 
		ntpr=10, 
		ntr=1, 
		restraintmask=':1-517 &amp;#38; !@H=', 
		restraint_wt=2.0 
		/ 
		注意：输入选项说明：imin=1 执行最小化；maxcyc=5000 最小化循环的最大次数；ncyc=2500 在 ncyc 次循环后，最小化方法从最陡下降法切换为共轭梯度法；cut=10.0 非键截断距离（Å）；ntb=1 施加周期性边界条件，恒定体积；ntc=1 对键长不施加约束（以获得更好的能量收敛性）；ntf=1 计算所有相互作用；ntpr=10 每 ntpr 步将能量信息输出至 out 文件；ntr=1 使用谐波势能限制指定原子的标志；restraintmask=':1-517 &amp;#38; !@H=' 指定受限制的原子；restraint_wt=2.0 位置限制的权重（kcal/mol∙Å2）。</li>
	</ol>
	</li>
	<li>执行能量最小化的第二阶段（最陡下降法 5,000 步 + 共轭梯度法 5,000 步），在无限制条件下优化整个体系。
	<ol>
		<li>使用 min2.in 和 protein_min1.rst 作为输入文件：pmemd -O -i min2.in -p protein.prmtop 
		-c protein_min1.rst -o protein_min2.out -r protein_min2.rst。</li>
		<li>在输入文件 min2.in 中指定选项： 
		&amp;#38;cntrl 
		imin=1, maxcyc=10000, ncyc=5000, 
		cut=10.0, ntb=1, 
		ntc=1, ntf=1, 
		ntpr=10 
		/</li>
	</ol>
	</li>
</ol>
3. 加热
注意：此阶段旨在将体系从 0 K 加热至 300 K。由于基于 PDB 文件的初始模型不包含原子速度信息，因此需要为原子分配初始速度。
<ol>
	<li>对蛋白质原子施加位置限制进行加热过程（50 ps，恒定体积）。
	<ol>
		<li>使用 heat.in 和 protein_min2.rst 作为输入文件：pmemd -O -i heat.in -p protein.prmtop 
		-c protein_min2.rst -o protein_heat.out 
		-r protein_heat.rst -x protein_heat.mdcrd -ref protein_min2.rst</li>
		<li>遵循所需参数：-x FILE 表示在分子动力学轨迹中保存的坐标集合。</li>
		<li>在输入文件 heat.in 中指定选项： 
		&amp;#38;cntrl 
		imin=0, irest=0, ntx=1, 
		nstlim=25000, dt=0.002, 
		ntc=2, ntf=2, 
		cut=10.0, ntb=1, 
		ntpr=500, ntwx=500, 
		ntt=3, gamma_ln=2.0, 
		tempi=0.0, temp0=300.0, 
		ntr=1, restraintmask=':1-517', 
		restraint_wt=1.0, 
		nmropt=1 
		/ 
		&amp;#38;wt TYPE='TEMP0', 
		istep1=0, istep2=25000, 
		value1=0.1, value2=300.0 / 
		&amp;#38;wt TYPE='END' / 
		注：输入参数说明：imin=0 执行分子动力学模拟；irest=0 开始新的模拟；ntx=1 不从 rst 文件中读取初始速度；nstlim=25000 分子动力学步数；dt=0.002 时间步长（ps）；ntc=2 使用 SHAKE 算法约束含氢键的键长；ntf=2 不计算被约束键的力；ntwx=500 每 ntwx 步将坐标写入 mdcrd 文件；ntt=3 使用 Langevin 恒温器；gamma_ln=2.0 Langevin 动力学的碰撞频率（ps−1）；tempi=0.0 初始温度；temp0=300.0 参考温度；nmropt=1 读取 &amp;#38;wt 名单中指定的参数；TYPE='TEMP0' 目标温度逐步变化。</li>
	</ol>
	</li>
</ol>
4. 平衡
注意：此步骤对于调节水的密度并获得蛋白质的平衡状态是必要的。
<ol>
	<li>在300 K下进行平衡，不施加任何限制（500 ps，恒压）。
	<ol>
		<li>使用 平衡液 和 protein_heat.rst 作为输入： pmemd -O -i 平衡.in -p 蛋白质.prmtop 
		-c protein_heat.rst -o protein_equil.out -r protein_equil.rst -x protein_equil.mdcrd.</li>
		<li>在输入文件中指定选项 平衡中: 
		&amp;#38;cntrl 
		imin=0，irest=1，ntx=5， 
		nstlim=250000, 
		dt=0.002, 
		ntc=2，ntf=2， 
		cut=10.0，ntb=2，ntp=1，taup=2.0 
		ntpr=1000，ntwx=1000，ntwr=50000， 
		ntt=3，gamma_ln=2.0， 
		temp0=300.0 
		/ 
		注意：输入选项： irest=1 重新启动模拟； ntx=5 坐标和速度从先前生成的文件中读取 rst 文件； ntb=2 周期性边界&#160;施加恒定压力； ntp=1 各向同性压力缩放； τau_p = 2.0 压力弛豫时间（ps）； ntwr=50000 每个 ntwr 步骤，这些 rst 文件已写入，确保从崩溃中恢复。</li>
	</ol>
	</li>
</ol>
5. 产量动态
<ol>
	<li>达到平衡后，进行恒压条件下的生产阶段分子动力学模拟（10 ns 或更长时间），并生成轨迹文件，用于后续分析蛋白质结构。
	<ol>
		<li>使用 prod.in 和 protein_equil.rst 作为输入文件：pmemd -O -i prod.in -p protein.prmtop 
		-c protein_equil.rst -o protein_prod.out -r protein_prod.rst -x protein_prod.mdcrd。</li>
		<li>在输入文件 prod.in 中指定选项： 
		&amp;#38;cntrl 
		imin=0, irest=1, ntx=5, 
		nstlim=5000000, 
		dt=0.002, 
		ntc=2, ntf=2, 
		cut=10.0, ntb=2, ntp=1, taup=2.0, 
		ntpr=1000, ntwx=1000, ntwr=50000, 
		ntt=3, gamma_ln=2.0, 
		temp0=300.0, ig=-1 
		/ 
		注意：可在计算机集群和超级计算机上使用 pmemd 的并行版本（pmemd.MPI）或 GPU 加速版本（pmemd.cuda）。长时间的分子动力学模拟可拆分为多个阶段依次执行。强烈建议每次重启模拟时设置 ig=-1（随机种子选项）。可使用 ptraj 或 cpptraj 程序对坐标轨迹进行分析和处理。更多详细信息，请参阅软件用户手册。</li>
	</ol>
	</li>
</ol>

<ol><li>Olsson, M., Zhivotovsky, B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Caspases+and+cancer%5BTitle%5D)+AND+%22Cell+Death+and+Differentiation%22%5BJournal%5D)">Caspases and cancer.</a> Cell Death and Differentiation. 18 (9), 1441-1449 (2011).</li>
<li>Lavrik, I. N., Golks, A., Krammer, P. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Caspase%3A+Pharmacological+manipulation+of+cell+death%5BTitle%5D)+AND+%22Journal+of+Clinical+Investigation%22%5BJournal%5D)">Caspase: Pharmacological manipulation of cell death.</a> Journal of Clinical Investigation. 115 (10), 2665-2672 (2005).</li>
<li>Degterev, A., Boyce, M., Yuan, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+decade+of+caspases%5BTitle%5D)+AND+%22Oncogene%22%5BJournal%5D)">A decade of caspases.</a> Oncogene. 22 (53), 8543-8567 (2003).</li>
<li>Pop, C., Salvesen, G. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Human+caspases%3A+Activation%2C+specificity%2C+and+regulation%5BTitle%5D)+AND+%22Journal+of+Biological+Chemistry%22%5BJournal%5D)">Human caspases: Activation, specificity, and regulation.</a> Journal of Biological Chemistry. 284 (33), 21777-21781 (2009).</li>
<li>Zamaraev, A. V., Kopeina, G. S., Prokhorova, E. A., Zhivotovsky, B., Lavrik, I. N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Post-translational+modification+of+caspases%3A+The+other+side+of+apoptosis+regulation%5BTitle%5D)+AND+%22Trends+in+Cell+Biology%22%5BJournal%5D)">Post-translational modification of caspases: The other side of apoptosis regulation.</a> Trends in Cell Biology. 27 (5), 322-339 (2017).</li>
<li>Pearlman, S. M., Serber, Z., Ferrell, J. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+mechanism+for+the+evolution+of+phosphorylation+sites%5BTitle%5D)+AND+%22Cell%22%5BJournal%5D)">A mechanism for the evolution of phosphorylation sites.</a> Cell. 147 (4), 934-946 (2011).</li>
<li>Zamaraev, A. V., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Requirement+for+Serine-384+in+Caspase-2+processing+and+activity%5BTitle%5D)+AND+%22Cell+Death+and+Disease%22%5BJournal%5D)">Requirement for Serine-384 in Caspase-2 processing and activity.</a> Cell Death and Disease. 11 (10), 825(2020).</li>
<li>Dagbay, K. B., Hill, M. E., Barrett, E., Hardy, J. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Tumor-associated+mutations+in+caspase-6+negatively+impact+catalytic+efficiency%5BTitle%5D)+AND+%22Biochemistry%22%5BJournal%5D)">Tumor-associated mutations in caspase-6 negatively impact catalytic efficiency.</a> Biochemistry. 56 (34), 4568-4577 (2017).</li>
<li>Ando, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Cancer-associated+missense+mutations+of+caspase-8+activate+nuclear+factor-%CE%BAB+signaling%5BTitle%5D)+AND+%22Cancer+Science%22%5BJournal%5D)">Cancer-associated missense mutations of caspase-8 activate nuclear factor-κB signaling.</a> Cancer Science. 104 (8), 1002-1008 (2013).</li>
<li>Case, D. A., et al. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aCase%2C+DA+%22AMBER+2020%22">AMBER 2020</a>. , University of California, San Francisco. (2020).</li>
<li>Salomon-Ferrer, R., Case, D. A., Walker, R. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((An+overview+of+the+Amber+biomolecular+simulation+package%5BTitle%5D)+AND+%22WIREs+Computational+Molecular+Science%22%5BJournal%5D)">An overview of the Amber biomolecular simulation package.</a> WIREs Computational Molecular Science. 3 (2), 198-210 (2013).</li>
<li>Roe, D. R., Cheatham, T. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((PTRAJ+and+CPPTRAJ%3A+Software+for+processing+and+analysis+of+molecular+dynamics+trajectory+data%5BTitle%5D)+AND+%22Journal+of+Chemical+Theory+and+Computation%22%5BJournal%5D)">PTRAJ and CPPTRAJ: Software for processing and analysis of molecular dynamics trajectory data.</a> Journal of Chemical Theory and Computation. 9 (7), 3084-3095 (2013).</li>
<li>Maier, J. A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((ff14SB%3A+Improving+the+accuracy+of+protein+side+chain+and+backbone+parameters+from+ff99SB%5BTitle%5D)+AND+%22Journal+of+Chemical+Theory+and+Computation%22%5BJournal%5D)">ff14SB: Improving the accuracy of protein side chain and backbone parameters from ff99SB.</a> Journal of Chemical Theory and Computation. 11 (8), 3696-3713 (2015).</li>
<li>Salomon-Ferrer, R., Götz, A. W., Poole, D., Le Grand, S., Walker, R. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Routine+microsecond+molecular+dynamics+simulations+with+AMBER+on+GPUs.+2.+Explicit+solvent+particle+mesh+ewald%5BTitle%5D)+AND+%22Journal+of+Chemical+Theory+and+Computation%22%5BJournal%5D)">Routine microsecond molecular dynamics simulations with AMBER on GPUs. 2. Explicit solvent particle mesh ewald.</a> Journal of Chemical Theory and Computation. 9 (9), 3878-3888 (2013).</li>
<li>Berman, H. M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+Protein+Data+Bank%5BTitle%5D)+AND+%22Nucleic+Acids+Research%22%5BJournal%5D)">The Protein Data Bank.</a> Nucleic Acids Research. 28 (1), 235-242 (2000).</li>
<li>Burley, S. K., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((RCSB+Protein+Data+Bank%3A+Powerful+new+tools+for+exploring+3D+structures+of+biological+macromolecules+for+basic+and+applied+research+and+education+in+fundamental+biology%2C+biomedicine%2C+biotechnology%2C+bioengineering+and+energy+sciences%5BTitle%5D)+AND+%22Nucleic+Acids+Research%22%5BJournal%5D)">RCSB Protein Data Bank: Powerful new tools for exploring 3D structures of biological macromolecules for basic and applied research and education in fundamental biology, biomedicine, biotechnology, bioengineering and energy sciences.</a> Nucleic Acids Research. 49, 437-451 (2021).</li>
<li>Martinez, X., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Molecular+graphics%3A+Bridging+structural+biologists+and+computer+scientists%5BTitle%5D)+AND+%22Structure%22%5BJournal%5D)">Molecular graphics: Bridging structural biologists and computer scientists.</a> Structure. 27 (11), 1617-1623 (2019).</li>
<li>Lee, T. S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((GPU-accelerated+molecular+dynamics+and+free+energy+methods+in+Amber18%3A+Performance+enhancements+and+new+features%5BTitle%5D)+AND+%22Journal+of+Chemical+Information+and+Modeling%22%5BJournal%5D)">GPU-accelerated molecular dynamics and free energy methods in Amber18: Performance enhancements and new features.</a> Journal of Chemical Information and Modeling. 58 (10), 2043-2050 (2018).</li>
<li>Jäger, R., Zwacka, R. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+enigmatic+roles+of+caspases+in+tumor+development%5BTitle%5D)+AND+%22Cancers%22%5BJournal%5D)">The enigmatic roles of caspases in tumor development.</a> Cancers. 2 (4), 1952-1979 (2010).</li>
</ol>

作者声明无利益冲突。

所述的分子动力学方法可用于对胱天蛋白酶（caspase）的野生型和突变型进行建模，借助生物分子模拟软件包实现。本文讨论了该方法中若干重要的技术问题。首先，需从蛋白质数据库（Protein Data Bank）中选择一个具有代表性的胱天蛋白酶晶体结构。值得注意的是，胱天蛋白酶的单体和二聚体形式均可接受。建议选择分辨率较高且缺失残基数量最少的结构。某些残基的质子化状态可在输入的PDB文件中手动设定。例如，在图1中，催化残基His277的Nε2原子上连接了一个氢原子（HIE形式）。其次，必须在蛋白质周围构建溶剂盒子，以进行后续的显式溶剂模拟。需要进行能量最小化，以优化所添加的氢原子和水分子的坐标。第三，在开始生产性动力学模拟之前，必须完成加热和平衡过程。在加热阶段，应确保蛋白质表面的空腔和结合口袋被水分子充分填充（如有必要，可增加分子动力学步数）。在平衡阶段，需通过分析蛋白质主链原子相对于初始位置的均方根偏差来确认其达到平衡构象12。第四，在生产性动力学模拟过程中，必须监测蛋白质结构的构象变化。为此，应通过拟合主链原子将轨迹帧与起始结构叠加，然后使用分子可视化工具（如VMD、PyMOL等）分析构象变化17。如有需要，可增加分子动力学步数，或将其拆分为多个连续执行的模拟段。第五，建议使用GPU加速版本的分子动力学模拟程序（pmemd.cuda），其性能显著优于传统的CPU实现方式14,18。
所述方法的一个可能局限性在于蛋白质数据库（Protein Data Bank）中某些caspase结构的可获得性。然而，在PDB中已有超过900条条目属于该类别 &#34;caspase&#34; 请求表明，caspase 的结构已被深入研究。当晶体结构的分辨率较低（即细节层次粗糙）或缺少对 caspase 功能至关重要的某些片段时，可能会出现一些困难。
总之，分子动力学（MD）建模方法已被证明在预测基因突变、翻译后修饰或分子间相互作用后蛋白质结构变化方面具有有效性。在细胞死亡领域应用MD模拟，为研究翻译后修饰调控半胱天冬酶（caspase）的分子机制提供了重要机遇。在此背景下，长时间的MD模拟可用于评估功能区域的动态行为以及氨基酸替换的潜在影响，从而阐明与突变或修饰的caspase-2及其他caspase相关的分子机制。例如，在胃肠道癌症以及中枢和周围神经系统肿瘤中已检测到起始型caspase基因（caspase-2/caspase-8/caspase-9/caspase-10）的错义突变19。综上所述，对caspase进行MD建模是一种强大的in silico方法，有助于理解调控程序性细胞死亡及其相关疾病发生的机制，并推动新型有效疗法的开发。

细胞凋亡是调控多细胞生物形态发生和组织稳态研究最为广泛的一种细胞过程。细胞凋亡可由多种外部或内部刺激引发，例如死亡受体的激活、细胞周期信号的紊乱、DNA损伤、内质网（ER）应激以及多种细菌和病毒感染1。半胱天冬酶（caspases）是细胞凋亡中的关键执行分子，通常根据其结构域特征及在caspase级联反应中的位置分为两类：起始型（initiator caspases，包括caspase-2、caspase-8、caspase-9和caspase-10）和效应型（effector caspases，包括caspase-3、caspase-6和caspase-7）2,3。当接收到细胞死亡信号时，起始型caspase会与衔接分子相互作用，促进其近端诱导的二聚化和自剪切，从而形成有活性的酶。效应型caspase则通过起始型caspase的切割而被激活，并通过裂解多种细胞底物来执行下游的凋亡程序4。
起始型和效应型caspase的成熟与功能受到多种不同细胞内机制的调控，其中翻译后修饰在细胞死亡调控中发挥着不可或缺的作用5。添加修饰基团（如磷酸化、亚硝基化、甲基化或乙酰化）或蛋白质（如泛素化或SUMO化）可改变caspase的酶活性，或影响其蛋白构象与稳定性，从而调控凋亡过程。定点突变技术被广泛用于研究潜在的翻译后修饰位点并阐明其功能。通常将假定的修饰位点突变为不能被进一步修饰的其他氨基酸，例如将可能被磷酸化的丝氨酸和苏氨酸突变为丙氨酸，或将赖氨酸的泛素化位点替换为精氨酸。另一种策略是引入特定氨基酸以模拟翻译后修饰状态（例如，谷氨酸和天冬氨酸已被用于模拟磷酸化的丝氨酸或苏氨酸）6。然而，某些位于活性位点附近或关键位置的氨基酸替换可能改变caspase的结构，干扰其催化活性，并抑制凋亡性细胞死亡7。在caspase基因的肿瘤相关错义突变中也可观察到类似效应。例如，caspase-6的肿瘤相关突变R259H导致底物结合口袋中环状结构的构象改变，降低了底物的有效催化转换率8。在头颈部鳞状细胞癌中发现的caspase-8的G325A氨基酸替换可能损害caspase-8的活性，进而调控核因子-kB（NF-kB）信号通路，促进肿瘤发生9。
为了评估氨基酸替换对半胱天冬酶（caspase）结构和功能的潜在影响，可采用分子建模方法。本文描述了使用生物分子模拟软件包（Amber）对野生型半胱天冬酶及其突变体进行建模的分子动力学（MD）方法。分子动力学方法能够展示引入突变后蛋白质结构的动态演化过程。Amber软件包最初由Peter Kollman团队开发，现已成为生物分子模拟中最常用的软件工具之一10,11,12,13。该软件分为两个部分：（1）AmberTools，是一组常用于系统准备（如原子类型分配、添加氢原子和显式水分子等）和轨迹分析的程序；（2）Amber，其核心是pmemd模拟程序。AmberTools为免费软件包（同时也是安装Amber的前提条件），而Amber则需通过独立的许可协议并支付费用获取。在超级计算机上或利用图形处理器（GPU）进行并行模拟，可显著提升蛋白质结构动力学科学研究的计算性能14。目前最新的可用软件版本为AmberTools21和Amber20，但所描述的实验方案同样适用于早期版本。

<table><tbody><tr><td>Amber20</td><td>加利福尼亚大学旧金山分校</td><td></td><td>用于分子动力学模拟的软件&lt;br/&gt; http://ambermd.org</td></tr><tr><td>AmberTools21</td><td>加利福尼亚大学旧金山分校</td><td></td><td>用于分子建模与分析的软件&lt;br/&gt; http://ambermd.org</td></tr></tbody></table>

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

Exploring Caspase Mutations and Post-Translational Modification by Molecular Modeling Approaches

细胞凋亡是一种程序性细胞死亡，可清除受损细胞，并调控多细胞生物的发育过程和组织稳态。半胱天冬酶（caspase）是一类半胱氨酸蛋白酶家族，在细胞凋亡的启动和执行过程中发挥关键作用。半胱天冬酶的成熟及其活性通过翻译后修饰以高度动态的方式进行精细调控。为了评估翻译后修饰的影响，通常将潜在修饰位点突变为不受任何修饰影响的氨基酸残基。例如，将丝氨酸残基替换为丙氨酸或天冬氨酸。然而，此类替换可能改变半胱天冬酶活性位点的构象，从而干扰其催化活性和细胞功能。此外，位于关键位置的其他氨基酸残基的突变也可能破坏半胱天冬酶的结构与功能，导致细胞凋亡过程紊乱。为避免使用突变残基所带来的困难，可直接采用分子建模方法来估算氨基酸替换对半胱天冬酶结构的潜在影响。本实验方案利用生物分子模拟软件包（Amber）和超级计算机资源，对野生型半胱天冬酶及其突变体进行建模，以检测突变对蛋白质结构和功能的影响。

本方案可直接应用于半胱天冬酶（caspase）的翻译后修饰或致病性突变的研究。本节展示了分子动力学（MD）建模的工作流程（图1），该方法已成功用于caspase-2的研究7。通过体外（in vitro）定点突变技术将潜在的磷酸化位点（Ser/Thr突变为Ala）并结合生物化学方法，研究发现Ser384Ala突变可阻止caspase-2的加工过程，并抑制其酶活性以及凋亡性细胞死亡的诱导（补充图1）。然而，质谱分析和Phos-Tag技术均未证实Ser384位点发生磷酸化7。为了阐明Ser384如何调控caspase-2的活性，研究人员对该蛋白的三维结构进行了分子动力学（MD）建模分析（图1）。
使用1pyo晶体结构构建了野生型和突变型caspase-2的分子模型（可用的caspase-2结构列于补充表1中）。Ser384Ala突变体通过去除Ser384残基中的Oγ原子而构建（在PDB中，Ser和Ala残基的区别在于是否含有Oγ原子）。晶体结构中通常不包含氢原子坐标，因此向蛋白质结构中添加了氢原子，随后用厚度为12 Å的TIP3P水层进行溶剂化，以进行后续的显式溶剂模拟。添加钠离子以中和体系电荷。所得的caspase-2起始模型根据分子动力学（MD）建模流程，使用min1.in、min2.in、heat.in、equil.in和prod.in控制数据文件，依次进行能量最小化、平衡和随后的10 ns MD模拟。
在caspase-2的晶体结构中，Ser384与Arg219和Arg378共同参与构成活性位点腔体的表面。Arg219和Arg378与底物的羧基形成氢键，而Ser384不直接与底物发生相互作用。分子动力学（MD）模拟证实，Ser384Ala取代未影响催化残基Cys320（亲核体）和His277（广义碱），但导致Arg378发生显著的构象变化。其胍基转向体相溶剂，使得Nε原子无法与底物的羧基形成关键的氢键（图1）。在天然酶中，Ser384的Oγ原子（部分负电荷）与Arg378的胍基（正电荷）之间存在静电相互作用。然而，丝氨酸的取代显然破坏了这一相互作用。因此，研究表明Ser384Ala取代影响了caspase-2活性位点中精氨酸残基对底物的识别，从而削弱了酶活性以及诱导凋亡性细胞死亡的能力。该发现的机制在进化上似乎保守，并在其他caspase家族成员中普遍存在。因此，分子动力学模拟的应用不仅验证了生化实验结果，还进一步揭示了caspase活性中心分子结构的新见解。
<img alt="蛋白质动力学示意图；野生型及Ser384Ala突变体的PDB结构；溶剂化、能量最小化。" class="xfigimg" src="/files/ftp_upload/64206/64206fig1v2.jpg" data-altupdatedat="1762425657000" data-alt="图1"> 
图1：用于研究caspase结构的分子动力学建模流程。 依据方案部分的说明，对野生型caspase-2及其Ser384Ala突变体进行了研究。结果显示，Ser384Ala替换在活性位点残基Arg378处诱导了重要的构象变化。 <a href="https://www.jove.com/files/ftp_upload/64206/64206fig01large.jpg" target="_blank">请点击此处查看该图的放大版本。</a>
补充图1：caspase-2在细胞死亡中的功能。 在外源性和内源性刺激作用下，野生型caspase-2可通过自蛋白水解机制被激活。活化的caspase-2切割Bid蛋白，进而促进线粒体外膜通透化，诱导细胞凋亡。Ser384Ala突变可阻止caspase-2的活化及其介导的细胞死亡。 <a href="https://www.jove.com/files/ftp_upload/64206/Supplementary%20Figure%201%20fin.jpg" target="_blank">请点击此处下载该文件。</a>
补充表 1：蛋白质数据库中可获取的 Caspase-2 结构。 结构按发布日期排序。<a href="https://www.jove.com/files/ftp_upload/64206/Supplementary%20Table%201.xlsx" target="_blank">请点击此处下载该表格。</a>
补充文件1：编辑PDB文件的不同步骤。 <a href="https://www.jove.com/files/ftp_upload/64206/Supplementary%20File%201.docx" target="_blank">请点击此处下载该文件。</a>

本实验方案采用一种生物分子模拟软件，描述了用于模拟野生型半胱天冬酶（caspase）及其突变体的分子动力学（MD）方法。分子动力学方法可用于评估半胱天冬酶结构的动态演化过程，以及突变或翻译后修饰可能产生的影响。

细胞凋亡是一种程序性细胞死亡，可清除受损细胞，并调控多细胞生物的发育过程和组织稳态。

在启动和执行过程中，利用分子建模发挥关键作用。我们研究翻译后修饰和突变对Caspase结构与功能的影响。我们描述了分子动力学方法，该方法可在原子水平上引入结构修饰后，揭示蛋白质演变过程的动态视图。这种方法对于理解调控程序性细胞死亡及其相关疾病的机制并开发新的有效疗法具有重要意义。我们展示了使用最流行的工具之一 Amber 20 进行分子建模的能力。但本研究所呈现的方案也可用于该软件的正式版本。以下视频逐步介绍了如何准备Caspase结构以进行模拟，并对该蛋白质的野生型和修饰型进行ascilica分析。要获取选定的蛋白质数据银行（PDB）结构，请使用下载文件下拉列表并点击PDB格式。从PDB文件中删除备注和连接性数据，并在独立的蛋白质链之间插入TER卡片。要准备起始模型，请从 Amber Tools 软件包中启动 T-LEAP 程序。然后在命令行中输入指定命令，加载 F F 14 S B 力场，以用于通过分子力学描述蛋白质，以及水分子和原子离子（如钠离子和氯离子）的参数。接下来，加载 P D B 文件，并为氢原子构建坐标，创建一个名为 mol 的对象。使用 check mol 命令检查可能导致问题的内部不一致性。在蛋白质周围创建溶剂盒。然后通过输入 charge mol 检查总电荷，并添加反离子以中和体系。输入指定命令以创建拓扑文件 P R M top 和坐标文件 I N P C R D。完成后，退出 T-lEAP 程序。进行第一阶段的能量最小化，以优化所添加的氢原子和水分子的位置，同时通过对重原子施加位置限制来固定蛋白质坐标。通过输入指定命令来运行 P M E M D 程序。按照我控制的数据提供所需的参数。P 分子拓扑结构、力场参数和原子名称。C 初始坐标。O 用户可读的日志输出。R 最终坐标。REF，位置限制的参考坐标。接下来，使用指定命令的输入参数，在无约束条件下进行能量最小化的第二阶段，以优化整个系统。该阶段的目标是将体系从零开尔文加热至300开尔文。通过使用给定命令作为输入，在恒定体积下对蛋白质原子施加位置限制，用50皮秒完成加热过程。按照分子动力学轨迹中保存的X坐标集合提供所需的参数。下一步需要调节水的密度，并获得蛋白质的平衡状态。在达到平衡后，使用指定命令在300开尔文、恒压条件下进行500皮秒的平衡，且不施加任何限制。在恒定压力下进行一次长达10纳秒或更长时间的生产阶段分子动力学模拟，并使用指定命令生成轨迹文件，用于后续对蛋白质结构的分析。可在计算机集群和超级计算机上使用程序的并行版本 PMEMD.MPI 或 GPU 加速版本 PMEMD.cuda。长时间的分子动力学模拟可拆分为多个阶段，依次连续执行。在本项分析中，通过分子动力学模拟工作流程研究了野生型Caspase-2及其丝氨酸-384丙氨酸突变体。丝氨酸-384到丙氨酸的替换在活性位点残基精氨酸-378处诱导了重要的构象变化。此外，研究显示，丝氨酸-384到丙氨酸的替换影响了活性位点中精氨酸残基对底物的识别，从而损害了caspase的催化活性。位点定向诱变和生化实验证明，丝氨酸-384到丙氨酸的突变阻断了Caspase-2的酶活性，抑制了癌细胞的凋亡。该方法可用于评估氨基酸突变及翻译后修饰对Caspase-2以及其他参与多种癌症疾病的Caspase家族蛋白功能的影响。

exploring-caspase-mutations-post-translational-modification-molecular

Research

科研

JoVE Journal

Biology

生物学

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

心肌梗死（MI）不仅会导致心脏功能受损，还会引发杏仁核的细胞凋亡，杏仁核是参与心肌梗死后行为学改变的脑区。本方案描述了如何诱导心肌梗死、采集杏仁核组织，并检测其中的caspase-3活性——一种细胞凋亡的标志物。...

通过光谱荧光法测定心肌梗死后大鼠杏仁核中的 Caspase-3 活性

为了在整体动物中检测含有低水平半胱天冬酶（caspase）活性的健康细胞，研究人员为黑腹果蝇（Drosophila）开发了高灵敏度的生物传感器，命名为CaspaseTracker。在无死亡刺激条件下、经优化培养的成年动物体内各内脏器官的长寿健康细胞中，均可检测到依赖半胱天冬酶的生物传感器活性。...

体内 生物传感器追踪非凋亡性半胱天冬酶活性 果蝇

本方案描述了半胱天冬酶双分子荧光互补（BiFC）技术；这是一种基于成像的方法，可用于可视化诱导起始半胱天冬酶的邻近效应，这是其激活的第一步。...

利用双分子荧光互补技术揭示半胱天冬酶激活通路

本方案描述了从人血液样本中获取单核细胞来源的巨噬细胞（MDM）的工作流程，一种在不损害细胞活力和行为的前提下，将炎症性半胱天冬酶双分子荧光互补（BiFC）报告系统高效导入人MDM的简便方法，以及一种基于成像的技术，用于在活细胞中检测炎症性半胱天冬酶的激活。...

人单核细胞来源巨噬细胞中炎症性半胱天冬酶诱导接近现象的可视化

甲型流感病毒（IAV）感染会激活半胱天冬酶（caspases），从而切割宿主和病毒蛋白，这些切割产物具有促病毒和抗病毒功能。通过使用抑制剂、RNA干扰、定点突变以及蛋白质印迹和逆转录定量PCR技术，研究人员鉴定了在感染的哺乳动物细胞中切割宿主蛋白cortactin和组蛋白去乙酰化酶的半胱天冬酶。...

鉴定在甲型流感病毒感染过程中切割蛋白质的半胱天冬酶及其基序

本文介绍了一种详细的免疫组织化学实验方案，用于在复杂组织中识别、验证和靶向具有功能相关性的半胱天冬酶（caspases）。...

视网膜组织中免疫染色Caspase-9的定量分析

本方案描述了一种全面的方法，用于检测在两种刺激条件下caspase激活（caspase-1、caspase-3、caspase-7、caspase-8、caspase-9 和 caspase-11）的反应 体外 和 体内...

评估半胱天冬酶激活以检测天然免疫细胞死亡

本文介绍了一种表达和纯化重组蛋白的实验方案 果蝇 caspases Dronc 和 Drice，及其应用 体外...

体外 使用纯化重组蛋白的切割实验 果蝇 用于底物筛选的Caspase酶

本实验方案描述了两种通过荧光底物结合流式细胞术或荧光分光光度计来检测caspase活性的方法。...

使用荧光定量法或流式细胞术检测Caspase活性

本文介绍了一种针对经活性羰基物种修饰的氨基酸残基进行优化与参数化的实验方案，该方案可适用于蛋白质体系。本方案的步骤包括结构设计与优化、电荷分配、参数构建以及蛋白质体系的准备。...

利用分子动力学模拟在计算机中合成带有反应性羰基的修饰氨基酸以评估结构效应

探索分子动力学模拟方法研究半胱天冬酶突变和翻译后修饰对蛋白质结构与功能的影响。该方法使用Amber 20工具，可评估氨基酸突变在癌症相关半胱天冬酶中的作用机制。

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

探索更多视频