シミュレーションの計算において、昨年実施したHER2ダイマーの初期構造決定の結果を解析した結果、病原性バリアントの解析を実施していくにあたって、初期構造としていくつかの出発構造を設定する事を考えた。特に、HER2ホモダイマーとHER2-EGFRヘテロダイマーはそのダイマー状態の構造の報告がなく、構造における信頼性を高めるために温度レプリカ交換分子動力学法(T-REMD)を実施して、それらのダイマーがとりうる構造パターンを広く探索した。
その結果、通常MDで行った結果から得られた複数の構造パターンによる空間内で、T-REMDによって得られる構造のパターンをほぼカバーできている事がわかった。したがって、T-REMDのシミュレーションの構造は熱で構造が壊れている可能性も考慮して、通常MDシミュレーションの構造を使用していく事を判断した。
HERタンパクにおける病的バリアントを対象としたシミュレーションを行っていくための初期構造として、通常MDで得られた構造から代表構造としてHER2ホモダイマーを5種類、HER2-EGFRヘテロダイマーを5種類、HER2-HER3ヘテロダイマーから3種類を選択した。それに対して野生型を含む13のバリアントに対して、169通りの組み合わせで、5 equilで1マイクロ秒のシミュレーションをそれぞて実施して計算を完了させた。 また、同様にHER2モノマーに関しても、野生型を含む13のバリアントに対しても、13種類のバリアントに対してのMDシミュレーションの計算を実施して完了した。
|