Name: カスパーゼ変異とダイナミの研究 - JoVEジャーナル
Uploaded: 2022-10-13T14:00:00.000+00:00
Duration: 5 min 56 s
Description: 本研究は、分子動力学シミュレーションを用いて突然変異や翻訳後修飾がカスパーゼの構造と機能にどのように影響するかを調べます。

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. システムの準備
注:天然および変異タンパク質形態の分子モデルは、タンパク質データバンク15,16から得られた適切な結晶構造に基づいて構築されています。
<ol><li>選択したPDB構造を取得するには、[ ファイルのダウンロード ]ドロップダウンリストを使用して、[ PDB形式]をクリックします。注釈と接続データを削除し、PDB ファイル内の個別のタンパク質鎖の間に TER カードを挿入します。必要に応じて、残基名HISをHIE、HID、またはHIP(補足ファイル1)に置き換えて、ヒスチジン残基のプロトン化状態を設定します。 注 : 結晶学的に分解された水分子 (PDB ファイルに存在する場合) を除去しないのが妥当です。</li>
	<li>開始モデルを準備するには、tleap プログラム (AmberTools パッケージ) を起動し、 tleap オブジェクトを操作するコマンドを入力します。<ol><li>プログラムを起動します: tleap。分子力学でタンパク質を記述するためにff14SB力場をロードします: ソースleaprc.protein.ff14SB>。</li>
		<li>水分子と原子イオン(Na+、Cl-)の荷重力場:>源leaprc.water.tip3p。</li>
		<li>PDB ファイルをロードし、水素の座標を構築して、mol: > mol = loadpdb protein.pdb という名前のオブジェクトを作成します。</li>
		<li>問題を引き起こす可能性のある内部の不整合を確認します。 > mol. 注: ジスルフィド ブリッジが存在する場合は、手動で指定する必要があります。共有結合したシステインCYSの名前をCYXに置き換え、 tleap で次のコマンドを入力して、SG原子間に結合を作成します: 結合 mol.X.SG mol.Y.SG (XとYはシステイン残基数です)。</li>
		<li>タンパク質の周囲に溶媒ボックスを作成します(TIP3P水、12 Å距離): >溶媒和ボックスモルTIP3PBOX 12.0。</li>
		<li>総電荷を確認します: >モルを充電し、対イオン(Na + またはCl-)を追加してシステムを中和します。 > 添加モルNa+ 0。</li>
		<li>トポロジ prmtop ファイルと座標 inpcrd ファイル (シミュレーションを実行するための入力) を作成します。 > saveamberparm mol protein.prmtop protein.inpcrd.tleap プログラムを終了する: >終了します。 注:アンバー形式の座標ファイルは、 ambpdb ツールを使用してPDB形式に簡単に変換できます(ユーザーズマニュアル、 材料表を参照)。</li>
	</ol></li>
</ol>2. エネルギー最小化
注意: エネルギー最小化は、MDの実行時に不安定になる始動システム内の原子間の悪い接触やオーバーラップを取り除くために必要です。
<ol><li>エネルギー最小化の第1段階(最急降下アルゴリズムの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 ファイルの初期座標;-o FILE ユーザーが読み取り可能なログ出力。-r ファイルの最終座標;-ref FILE 位置拘束の参照座標。</li>
		<li>入力ファイル min1.in でオプションを指定します。 &cntrl imin=1, maxcyc = 5000, ncyc=2500, カット= 10.0、NTB = 1、 NTC = 1、NTF = 1、 ntpr=10、 ntr=1、 restraintmask=':1-517 & !@H=', restraint_wt=2.0 / 注:入力オプション:imin = 1は最小化を実行します。maxcyc=5000 最小化サイクルの最大数。NCYC = 2500最小化法は、NCYCサイクル後に最急降下から共役勾配に切り替えられます。カット= 10.0非結合カットオフ(Å);ntb = 1周期的な境界が課され、一定のボリュームがあります。ntc=1 結合長に制約は適用されません(エネルギー収束を改善するため)。NTF = 1すべての相互作用が計算されます。NTPR = 10 エネルギー情報は、NTPR ステップごとに OUT ファイルに出力されます。調和ポテンシャルを使用して特定の原子を抑制するためのntr = 1フラグ。restraintmask=':1-517 & !@H=' は拘束された原子を指定する。位置拘束のためのrestraint_wt = 2.0重量(kcal / mol∙Å2)。</li>
	</ol></li>
	<li>エネルギー最小化の第2段階(最も急な降下ステップ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 でオプションを指定します。 &cntrl imin=1, maxcyc=10000, ncyc=5000, カット= 10.0、NTB = 1、 NTC = 1、NTF = 1、 ntpr=10 /</li>
	</ol></li>
</ol>3.暖房
注:この段階は、システムを0 Kから300 Kに加熱することを目的としています。 PDBファイルに基づく開始モデルには速度情報が含まれていないため、初期速度は原子に割り当てられます。
<ol><li>タンパク質原子の位置拘束(50ps、定容)で加熱処理を行います。<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 座標セットは、MD 軌道上に保存されます。</li>
		<li>入力ファイル heat.in でオプションを指定します。 &cntrl imin=0, irest=0, ntx=1, nstlim=25000, dt=0.002, NTC = 2、NTF = 2、 カット= 10.0、NTB = 1、 NTPR = 500、NTWX = 500、 NTT = 3, gamma_ln = 2.0, tempi=0.0, temp0=300.0, NTR=1, 拘束マスク=':1-517', restraint_wt=1.0、 nmropt=1 / &wt 型='TEMP0', istep1=0, istep2=25000, 値 1=0.1、値 2=300.0 / &wt タイプ='終了' / メモ: 入力オプション: imin=0 は MD を実行します。irest=0 新しいシミュレーションを実行します。ntx=1 RST ファイルから初期速度は読み取られません。nstlim=25000 MD ステップ数。dt=0.002 時間ステップ (ps);水素を含むntc = 2結合は、SHAKEアルゴリズムを使用して拘束されます。拘束された結合のNTF = 2力は計算されません。NTWX = 500座標は、NTWXステップごとにMDCRDファイルに書き込まれます。NTT=3ランジュバン温度調節;gamma_ln=2.0 ランジュバンダイナミクスの衝突頻度 (ps−1);温度= 0.0初期温度;温度0 = 300.0基準温度;&wt 名前リストに指定された nmropt=1 パラメーターが読み取られます。TYPE='TEMP0' は目標温度を変化させます。</li>
	</ol></li>
</ol>4.平衡化
注:この段階は、水の密度を調整し、タンパク質の平衡状態を取得するために必要です。
<ol><li>拘束なしで300 Kで平衡化を実行します(500 ps、定圧)。<ol><li>equil.in と protein_heat.rst を入力として使用します。pmemd -O -i equil.in -p protein.prmtop -c protein_heat.rst -o protein_equil.out -r protein_equil.rst -x protein_equil.mdcrd.</li>
		<li>入力ファイル equil.in でオプションを指定します。 &cntrl imin=0, irest=1, ntx=5, nstlim=250000、 dt=0.002、 NTC = 2、NTF = 2、 カット= 10.0、NTB = 2、NTP = 1、タウプ= 2.0、 NTPR = 1000、NTWX = 1000、NTWR = 50000、 NTT = 3, gamma_ln = 2.0, 温度0=300.0 / 注: 入力オプション: irest=1 シミュレーションを再起動します。 NTX = 5 の座標と速度は、以前に生成された RST ファイルから読み取られます。 NTB = 2 周期的な境界が課され、一定の圧力がかかります。 NTP = 1 等方性圧力スケーリング; タウプ=2.0 圧力緩和時間(ps); ntwr=50000 NTWR ステップごとに RST ファイルが書き込まれ、クラッシュから確実に回復します。</li>
	</ol></li>
</ol>5.生産ダイナミクス
<ol><li>平衡が達成されたら、一定の圧力で生産MDシミュレーション(10ns以上)を実行し、その後のタンパク質構造の解析のために軌道ファイルを生成します。<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 でオプションを指定します。 &cntrl imin=0, irest=1, ntx=5, nstlim=5000000, dt=0.002、 NTC = 2、NTF = 2、 カット= 10.0、NTB = 2、NTP = 1、タウプ= 2.0、 NTPR = 1000、NTWX = 1000、NTWR = 50000、 NTT = 3, gamma_ln = 2.0, 温度0=300.0, ig=-1 / 注: pmemd の並列バージョン (pmemd.MPI) または GPU アクセラレーション バージョン (pmemd.cuda) は、コンピューター クラスターおよびスーパーコンピューターで使用できます。長いMDシミュレーションは、いくつかのセグメントに分割され、順番に実行される場合があります。シミュレーションを再起動するたびにig=-1(ランダムシードオプション)を設定することを強くお勧めします。ptraj または cppptraj プログラムは、座標軌道の解析と処理に使用できます。詳細については、ソフトウェアのユーザーズマニュアルを参照してください。</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>

著者は開示する利益相反を持っていません。

記載されたMDアプローチは、生体分子シミュレーションパッケージを使用して、野生型および変異型の両方のカスパーゼをモデル化することを可能にする。ここでは、方法論のいくつかの重要な問題について説明します。まず、カスパーゼの代表的な結晶構造をタンパク質データバンクから選択する必要があります。重要なことに、カスパーゼの単量体および二量体形態の両方が許容される。欠落残基の数を最小限に抑えた高解像度の構造を選択することをお勧めします。一部の残基のプロトン化状態は、入力PDBファイルで手動で設定できます。例えば、図1において、触媒残基His277のNε2原子(HIE型)に水素原子が結合していた。次に、さらなる陽解法溶媒シミュレーションのために、タンパク質の周囲に溶媒ボックスを作成する必要があります。追加された水素原子と水分子の座標を最適化するには、エネルギーの最小化が必要です。第三に、加熱と平衡化は、生産ダイナミクスシミュレーションの前に実行する必要があります。加熱段階では、タンパク質表面の空洞と結合ポケットが水分子で満たされていることを確認する必要があります(必要に応じて、MDステップの数を増やすことができます)。平衡化段階では、タンパク質の平衡立体配座の獲得は、初期位置からのその骨格原子の二乗平均平方根偏差を分析することによって確認される必要がある12。第四に、生産動態シミュレーション中にタンパク質構造の立体構造変化を監視する必要があります。そのためには、骨格原子をフィッティングして軌道系を出発構造に重ね合わせ、分子可視化ツール(VMD、PyMOLなど)を用いて立体構造変化を解析する必要があります。17. 必要に応じて、MDステップの数を増やしたり、シミュレーションを複数のセグメントに分割して順番に実行したりできます。第五に、MDシミュレーションプログラムのGPUアクセラレーションバージョン(pmemd.cuda)を使用することが推奨されており、そのパフォーマンスは従来のCPU実装で達成可能なパフォーマンスを大幅に上回っています14,18。
記載された方法の考えられる制限は、タンパク質データバンクにおける特定のカスパーゼ構造の可用性である。しかし、「カスパーゼ」要求の下で900を超えるPDBエントリが見つかり、カスパーゼ構造が詳細に調査されたことを示しています。結晶構造が低分解能(すなわち、粗い詳細レベル)で得られる場合、またはカスパーゼ機能にとって重要ないくつかの断片を欠く場合、いくつかの困難が生じる可能性があります。
要約すると、MDモデリングアプローチは、遺伝子の突然変異、修飾、または分子間相互作用後のタンパク質の構造変化を予測するのに効果的であることが証明されています。細胞死の分野におけるMDシミュレーションの適用は、翻訳後修飾によるカスパーゼ制御の分子メカニズムを研究する大きな機会を開きます。この点に関して、長期MDシミュレーションは、機能領域の動的挙動およびアミノ酸置換の潜在的な影響を評価し、変異または修飾されたカスパーゼ-2および他のカスパーゼに関連する分子的原因を解明することができる。例えば、イニシエーターカスパーゼ遺伝子(カスパーゼ-2/カスパーゼ-8/カスパーゼ-9/カスパーゼ-10)のミスセンス変異は、消化器癌および中枢神経系および末梢神経系の腫瘍において検出された19。同時に、カスパーゼのMDモデリングは、プログラムされた細胞死および関連する障害を制御するメカニズムを理解し、新しい効果的な治療法を開発するための強力な インシリコ アプローチです。

アポトーシスは、多細胞生物の形態形成と組織恒常性を調節する最も広く研究されている細胞プロセスの1つです。アポトーシスは、死受容体の活性化、細胞周期シグナルの乱れ、DNA損傷、小胞体(ER)ストレス、さまざまな細菌およびウイルス感染など、広範囲の外部または内部刺激によって開始される可能性があります1。カスパーゼ(主要なアポトーシスプレーヤー)は、ドメイン構造とカスパーゼカスケード2,3の場所に応じて、従来、イニシエーター(カスパーゼ-2、カスパーゼ-8、カスパーゼ-9、およびカスパーゼ-10)とエフェクター(カスパーゼ-3、カスパーゼ-6、およびカスパーゼ-7)の2つのグループに分類されます。細胞死シグナルにより、開始剤カスパーゼはアダプター分子と相互作用し、近接誘導二量体化および自動処理を促進して活性酵素を形成します。エフェクターカスパーゼは、開始剤カスパーゼによる切断を介して活性化され、複数の細胞基質4を切断することによって下流の実行ステップを実行する。
開始剤およびエフェクターカスパーゼの成熟および機能は、多数の異なる細胞内メカニズムによって調節されており、その中で翻訳後修飾は細胞死調節に不可欠な役割を果たす5。修飾基(リン酸化、ニトロシル化、メチル化、アセチル化)またはタンパク質(ユビキチン化またはSUMO化)の付加は、カスパーゼの酵素活性またはアポトーシスを調節するタンパク質の立体構造および安定性を変化させます。部位特異的突然変異誘発は、潜在的な翻訳後修飾部位を調査し、それらの役割を識別するために広く適用されています。推定修飾部位は通常、別のアミノ酸に置き換えられますが、これはさらに修飾することはできません。したがって、潜在的にリン酸化されたセリンおよびスレオニンはアラニンに変異し、リジンユビキチン化部位はアルギニンに置換される。別の戦略は、特に翻訳後修飾を模倣するアミノ酸を置換することを含む(例えば、リン酸化セリンまたはスレオニンを模倣するために使用されているグルタミン酸およびアスパラギン酸)6。しかし、活性部位の高い近傍または臨界位置に位置するこれらの置換のいくつかは、カスパーゼ構造を変化させ、触媒活性を乱し、アポトーシス細胞死を抑制する可能性があります7。同様の効果は、カスパーゼ遺伝子における腫瘍関連ミスセンス変異の場合にも観察され得る。例えば、カスパーゼ−6−R259H−の腫瘍関連変異は、基質結合ポケット内のループの立体構造変化をもたらし、基質8の効率的な触媒代謝回転を減少させた。頭頸部扁平上皮癌で同定されたカスパーゼ-8のG325Aアミノ酸置換は、カスパーゼ-8活性を阻害し、核因子-kB(NF-kB)シグナル伝達の調節につながり、腫瘍形成を促進しました9。
カスパーゼの構造と機能に対するアミノ酸置換の潜在的な影響を評価するために、分子モデリングを適用することができます。分子動力学(MD)アプローチは、生体分子シミュレーションパッケージ(Amber)を使用して野生型カスパーゼとその変異型をモデル化するためのこの研究で説明されています。MD法は、突然変異の導入後のタンパク質構造の動的進化のビューを提供します。もともとPeter Kollmanのグループによって開発されたAmberパッケージは、生体分子シミュレーション用の最も人気のあるソフトウェアツールの1つになりました10,11,12,13。このソフトウェアは2つの部分に分かれています:(1)AmberTools、システムの準備(原子タイプの割り当て、水素と明示的な水分子の追加など)と軌道分析に日常的に使用されるプログラムのコレクション。(2)pmemdシミュレーションプログラムを中心とするアンバー。AmberTools は無料のパッケージ (および Amber 自体をインストールするための前提条件) ですが、Amber は別のライセンスと料金体系で配布されます。スーパーコンピュータ上での並列シミュレーションやグラフィックスプロセッシングユニット(GPU)の使用により、タンパク質構造ダイナミクスの科学的研究のパフォーマンスを大幅に向上させることができます14。利用可能な最新のソフトウェアバージョンはAmberTools21とAmber20ですが、説明されているプロトコルは以前のバージョンでも使用できます。

<table><tbody><tr><td>アンバー20</td><td>カリフォルニア大学サンフランシスコ校</td><td></td><td>分子動力学シミュレーション用ソフトウェア http://ambermd.org</td></tr><tr><td>アンバーツール21</td><td>カリフォルニア大学サンフランシスコ校</td><td></td><td>分子モデリングおよび分析用ソフトウェア http://ambermd.org</td></tr></tbody></table>

分子モデリングアプローチによるカスパーゼ変異と翻訳後修飾の探索

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

アポトーシスは、損傷した細胞を排除し、多細胞生物の発生と組織の恒常性を制御するプログラムされた細胞死の一種です。システインプロテアーゼのファミリーであるカスパーゼは、アポトーシスの開始と実行において重要な役割を果たします。カスパーゼの成熟とその活性は、非常に動的な方法で翻訳後修飾によって微調整されます。翻訳後変化の影響を評価するために、潜在的な部位は、修飾に永続的な残基で日常的に変異されます。例えば、セリン残基はアラニンまたはアスパラギン酸で置換される。しかし、そのような置換はカスパーゼ活性部位の立体構造を変化させ、触媒活性および細胞機能の障害をもたらす可能性がある。さらに、重要な位置に位置する他のアミノ酸残基の突然変異もカスパーゼの構造と機能を破壊し、アポトーシス摂動を引き起こす可能性があります。変異残基を採用することの難しさを回避するために、分子モデリングアプローチを容易に適用して、カスパーゼ構造に対するアミノ酸置換の潜在的な影響を推定することができます。本プロトコルは、生物分子シミュレーションパッケージ(Amber)とスーパーコンピューター機能を使用して野生型カスパーゼとその変異型の両方をモデル化し、タンパク質の構造と機能に対する突然変異の影響をテストすることを可能にします。

本プロトコルは、カスパーゼまたは病原性変異の翻訳後修飾の研究に容易に適用することができる。このセクションでは、カスパーゼ-2の研究で成功裏に使用されているMDモデリングワークフローが示されています(図1)。潜在的なリン酸化部位(Ser/ThrからAla)のin vitro部位特異的突然変異誘発と生化学的アプローチを使用して、Ser384Ala変異がカスパーゼ-2プロセシングを妨げ、酵素活性とアポトーシス細胞死の誘導をブロックすることが実証されました(補足図1)。しかし、質量分析もPhos-Tag技術も、Ser384がリン酸化を受けることを確認しませんでした7。Ser384がカスパーゼ-2活性をどのように調節するかを理解するために、3次元タンパク質構造のMDモデリングを実施しました(図1)。
1pyo結晶構造を用いてカスパーゼ-2の野生型および変異型の分子モデルを構築した(利用可能なカスパーゼ-2構造は 別表1に記載)。Ser384Ala変異体は、Ser384残基のO γ原子を除去して作成された(PDBにおいて、 Ser残基およびAla残基はOγ 原子を有することで区別される)。水素座標は通常、結晶構造には存在しません。そこで、タンパク質の構造に水素原子を付加し、さらに陽解な溶媒シミュレーションを可能にするために、12 Åの厚さのTIP3P水によって溶媒和しました。ナトリウムイオンを添加して系を中和した。得られたカスパーゼ-2の出発モデルは、 min1.in、 min2.in、 heat.in、 equil.in、および prod.in 制御データファイルを使用して、MDモデリングプロトコルに従ってエネルギー最小化、平衡化、およびその後の10nsMDシミュレーションにかけられました。
カスパーゼ-2の結晶構造では、Ser384がArg219およびArg378とともに、活性サイトでの空洞表面の形成に関与しています。Arg219とArg378は基質のカルボキシル基と水素結合を形成しますが、Ser384は基質との直接的な相互作用には介入しません。MDシミュレーションでは、Ser384Ala置換は触媒残基Cys320(求核剤)およびHis277(一般塩基)に影響を与えなかったが、Arg378に大きな立体構造変化を誘導することが確認された。そのグアニジニウム基はバルク溶媒に向かったため、Nε 原子は基板のカルボキシル基と必須の水素結合を形成できませんでした(図1)。天然酵素では、Ser384のOγ 原子(部分的な負電荷)とArg378グアニジニウム基(正電荷)との間に静電相互作用が起こる。しかしながら、セリン置換は明らかにこの相互作用を混乱させる。したがって、Ser384Ala置換は、カスパーゼ-2活性部位のアルギニン残基による基質認識に影響を与え、酵素活性とアポトーシス細胞死を引き起こす能力を損なうことが示されました。発見されたメカニズムは進化的に保存されており、他のカスパーゼファミリーのメンバーにも共通しているようです。このように、MDシミュレーションの実施により、生化学的結果を確認し、カスパーゼの活性中心の分子構造に関する新たな洞察を得ることができました。
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/64206/64206fig1v2.jpg"> 図1:カスパーゼ構造の研究に使用されるMDモデリングワークフロー。 野生型カスパーゼ-2およびそのSer384Ala変異体を、プロトコルセクションの指示に従って調査しました。Ser384Ala置換は、活性部位残基Arg378の重要な立体構造変化を誘導することが明らかになった。 <a href="https://www.jove.com/files/ftp_upload/64206/64206fig01large.jpg" target="_blank">この図の拡大版を表示するには、ここをクリックしてください。</a>
補足図1:細胞死におけるカスパーゼ-2の機能。 外因性および内因性の刺激に応答して、野生型カスパーゼ-2は自己タンパク質分解機構によって活性化され得る。活性カスパーゼ-2はBidを切断し、ミトコンドリア外膜透過処理とアポトーシス細胞死を促進します。Ser384Ala変異は、カスパーゼ-2の活性化と細胞死の誘導を防ぎます。 <a href="https://www.jove.com/files/ftp_upload/64206/Supplementary%20Figure%201%20fin.jpg" target="_blank">このファイルをダウンロードするには、ここをクリックしてください。</a>
補足表1:タンパク質データバンクで利用可能なカスパーゼ-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>

本プロトコルは、生体分子シミュレーションパッケージを使用し、野生型カスパーゼおよびその変異形態をモデル化するための分子動力学(MD)アプローチを記載する。MD法は、カスパーゼ構造の動的進化および突然変異または翻訳後修飾の潜在的な影響を評価することを可能にする。

Apoptosis is a type of programmed cell death that eliminates damaged cells and controls the development and tissue homeostasis of multicellular organisms. ...

分子モデリングを使用した開始と実行において重要な役割を果たします。翻訳後修飾と突然変異がカスパーゼの構造と機能に及ぼす影響を研究しています。原子レベルでの構造修飾の導入に続くタンパク質の進化を視野に入れる分子動力学アプローチについて説明します。このようなアプローチは、プログラム細胞死および関連する障害を調節するメカニズムを理解し、新しい効果的な治療法を開発する上で特に重要です。最も一般的なツールの1つであるAmber 20を使用して、分子モデリング機能を実証します。ただし、提示されたプロトコルは、このソフトウェアの正式なバージョンでも使用できます。次のビデオでは、シミュレーション用のカスパーゼ構造を準備し、このタンパク質の野生型および修飾型のアシリカ調査を実行する手順を段階的に説明します。選択したタンパク質データバンクまたはPDB構造を取得するには、ダウンロードファイルのドロップダウンリストを使用して、P D B形式をクリックします。注釈と接続データを削除し、PDBファイル内の別々のタンパク質鎖の間にT E Rカードを挿入します。開始モデルを準備するには、アンバーツールパッケージからT-LEAPプログラムを起動します。次に、F 14 S B力場をロードして、コマンドラインに示されたコマンドを入力して、分子力学と水分子とナトリウムや塩化物などの原子イオンのパラメータでタンパク質を記述します。次に、P D Bファイルをロードし、水素の座標を構築してmolという名前のオブジェクトを作成します。check molコマンドを使用して、問題を引き起こす可能性のある内部不整合を確認します。タンパク質の周りに溶媒ボックスを作成します。次に、電荷molを入力して総電荷を確認し、対イオンを追加して系を中和します。指定されたコマンドを入力して、トポロジ P R M トップ ファイルと座標 I N P C R D ファイルを作成します。完了したら、T-lEAPプログラムを終了します。エネルギー最小化の第1段階を実行して、重原子の位置拘束によってタンパク質座標を固定しながら、追加された水素原子と水分子の位置を最適化します。示されたコマンドを入力して、P M E M D プログラムを実行します。データを制御するときに必要な引数に従ってください。P分子トポロジー、力場パラメータ、原子名。C の初期座標。O ユーザーが読み取り可能なログ出力 R 最終座標。REF、位置拘束の参照座標。次に、エネルギー最小化の第2段階を制約なしで実行し、指定されたコマンドの入力を使用してシステム全体を最適化します。この段階は、システムを0ケルビンから300ケルビンに加熱することを目的としています。与えられたコマンドを入力として、一定体積で50ピコ秒でタンパク質原子の位置拘束で加熱プロセスを実行します。分子動力学軌道上に保存されたX座標セットとして、必要な引数に従います。次の段階は、水の密度を調整し、タンパク質の平衡状態を得るために必要である。平衡が正常に達成された後、示されたコマンドを使用して、一定の圧力で300ケルビンで500ピコ秒の平衡化を拘束なしで実行します。一定の圧力で10ナノ秒以上の生産分子動力学シミュレーションを実行し、示されたコマンドを使用してタンパク質構造の後続の分析のための軌道ファイルを生成します。プログラムの並列バージョンであるP M E M D M P IまたはG P U 加速版P M E M D cudaは、コンピュータクラスタおよびスーパーコンピュータで使用できます。長い分子動力学シミュレーションは、いくつかのセグメントに分割され、順番に実行される場合があります。この解析では、野生型カスパーゼ2とそのセリン-384アラニン変異体を分子動力学モデリングワークフローに従って調査しました。セリン-384アラニン置換は、活性部位残基アルギニン-378に重要な立体配座変化を誘導した。さらに、セリン-384アラニン置換は、キャスト塩基活性を損なう活性部位のアルギニン残基による基質認識に影響を与えることが示された。部位指向突然変異誘発および生化学的試験は、セリン-384アラニン変異がカスパーゼ2のプロセシングにおける酵素活性をブロックし、癌細胞のアポトーシスを抑制することを示した。記載されたアプローチは、カスパーゼ2ならびに様々な癌疾患に関与する他のカスパーゼにおけるアミノ酸突然変異および翻訳後修飾の影響を評価するために使用することができる。

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

Research

JoVE Journal

Biology

生物学

心筋梗塞(MI)は、心機能の障害だけでなく、心筋梗塞の行動結果に関与する脳領域である扁桃体のアポトーシスも引き起こします。このプロトコルは、MIを誘発する方法、扁桃体組織を採取する方法、およびアポトーシスのマーカーであるカスパーゼ-3活性を測定する方法を説明しています。...

心筋梗塞後の分光分析により測定したラットの扁桃体におけるカスパーゼ-3活性

カスパーゼ活性の低いレベルを含む全動物に健康な細胞を検出するために、CaspaseTracker指定された高感度バイオセンサーは、 ショウジョウバエのために生成されました。カスパーゼ依存性のバイオセンサーの活性は、死刺激の非存在下で最適化された条件下で飼育した成体動物の臓器全体長寿命健康な細胞で検出されます。...

インビボでのバイオセンサーは、 ショウジョウバエの非アポトーシスカスパーゼ活性を追跡します

このプロトコルを記述するカスパーゼ分子蛍光補完 (BiFC);自分の活性化の最初のステップであるイニシエーター カスパーゼの誘起近接を視覚化に使用することができますイメージング...

カスパーゼ活性化分子蛍光補完への道をライトアップ

このプロトコルは、ヒト血液サンプルから単球由来マクロファージ(MDM)を取得するワークフロー、細胞の生存率と挙動を損なうことなく炎症性カスパーゼ二分子蛍光相補(BiFC)レポーターをヒトMDMに効率的に導入する簡単な方法、および生細胞における炎症性カスパーゼ活性化を測定するためのイメージングベースのアプローチを記述しています。...

ヒト単球由来マクロファージにおける炎症性カスパーゼ誘導近接の可視化

インフルエンザAウイルス(IAV)感染は、宿主およびウイルスタンパク質を切断するカスパーゼを活性化し、次に、プロウイルスおよび抗ウイルス機能を有する。阻害剤、RNA干渉、部位特異的突然変異誘発、およびウェスタンブロッティングおよびRT-qPCR技術を採用することにより、宿主コルタクチンおよびヒストンデアセチラーゼを切断する感染哺乳類細胞におけるカスパーゼが同定された。...

A型インフルエンザウイルス感染時にタンパク質を切断するカスパーゼとそのモチーフの同定

ここでは、複雑な組織における機能的に関連するカスパーゼを同定、検証、標的化するための詳細な免疫組織化学プロトコルを示します。...

網膜組織における免疫染色カスパーゼ-9の定量

このプロトコルは、パイロトーシス、アポトーシス、ネクロプトーシス、およびパンノプトーシスなどの細胞死経路の開始を決定するための感染、無菌の侮辱、および癌のin vitro および in vivo...

自然免疫細胞死を評価するためのカスパーゼ活性化の評価

ここでは、組換え ショウジョウバエ カスパーゼDroncおよびDriceを発現および精製するためのプロトコル、および in vitro...

インビトロ 基質スクリーニングのための精製組換え ショウジョウバエ カスパーゼを用いた切断アッセイ

本プロトコルは、フローサイトメトリーまたは分光蛍光光度計を使用して蛍光発生基質を介してカスパーゼ活性を測定する2つの方法を記載する。...

蛍光アッセイまたはフローサイトメトリーを使用したカスパーゼ活性の測定

ここでは、タンパク質システムに適応可能な反応性カルボニル種で修飾されたアミノ酸残基の最適化とパラメータ化のためのプロトコールについて説明します。プロトコールのステップには、構造の設計と最適化、電荷の割り当て、パラメータの構築、およびタンパク質システムの調製が含まれます。...

反応性カルボニルをin silicoで修飾したアミノ酸を合成し、分子動力学シミュレーションを用いた構造効果の評価 

本研究は、分子動力学シミュレーションを用いて突然変異や翻訳後修飾がカスパーゼの構造と機能にどのように影響するかを調べます。

分子モデリングアプローチによるカスパーゼ変異と翻訳後修飾の探索

さらに動画を探す