分子シミュレーション関連:amber:解析方法

差分

このページの2つのバージョン間の差分を表示します。

この比較画面へのリンク

両方とも前のリビジョン 前のリビジョン
次のリビジョン		 前のリビジョン
分子シミュレーション関連:amber:解析方法 [2021/10/15 14:53] – [距離解析] mumeiyamibito分子シミュレーション関連:amber:解析方法 [2022/08/18 13:48] (現在) – mumeiyamibito
行 1: 行 1:
 ====== 解析方法 ====== ====== 解析方法 ======
 ===== RMSD 解析 ===== ===== RMSD 解析 =====
-  * 指定された対象の RMSD の経時変化を解析する +  * 指定された対象の RMSD の経時変化を解析する。 
-  * $RMSD = \frac {1} {N} \sqrt{\sum_{i = 1}^{N} (a_i - b_i)^2}$ +  * $RMSD = \sqrt{\frac {\sum_{i = 1}^{N} (a_i - b_i)^2}{N}}$ 
-  * ''cpptraj'' に同梱のパッケージで解析ができる+  * ''cpptraj'' に同梱のパッケージで解析ができる。
   * パラメータファイル (''test.prmtop'') と座標ファイル (''test.inpcrd'' / ''*.nc'' / ''*.mdcrd'') が必要   * パラメータファイル (''test.prmtop'') と座標ファイル (''test.inpcrd'' / ''*.nc'' / ''*.mdcrd'') が必要
-  * 対象 ''TARGET_MASK'' の RMSD を測定する\\ <code bash>+  * 対象 ''TARGET_MASK'' の RMSD を測定する。\\ <code bash>
 $ cpptraj $ cpptraj
 > parm test.prmtop > parm test.prmtop
行 13: 行 13:
 </code> </code>
     * ''TARGET_MASK'': 対象のマスク     * ''TARGET_MASK'': 対象のマスク
-    * ''first'': 最初のフレームをリファレンスにする+    * ''first'': 最初のフレームをリファレンスにする。
     * ''OUTPUT'': 出力ファイル     * ''OUTPUT'': 出力ファイル
     * ''time'': 各時間の刻み幅     * ''time'': 各時間の刻み幅
行 19: 行 19:
  
 ===== RMSF解析 ===== ===== RMSF解析 =====
-  * 指定された対象の RMSF を解析する+  * 指定された対象の RMSF を解析する。
   * RMSF はゆらぎのパラメータ   * RMSF はゆらぎのパラメータ
-  * ''cpptraj'' に同梱のパッケージで解析ができる+  * ''cpptraj'' に同梱のパッケージで解析ができる。
   * パラメータファイル (''test.prmtop'') と座標ファイル (''test.inpcrd'' / ''*.nc'' / ''*.mdcrd'') が必要\\ <code bash>   * パラメータファイル (''test.prmtop'') と座標ファイル (''test.inpcrd'' / ''*.nc'' / ''*.mdcrd'') が必要\\ <code bash>
 $ cpptraj $ cpptraj
行 33: 行 33:
  
 ===== 水素結合解析 ===== ===== 水素結合解析 =====
-  * 指定された対象の水素結合数の経時変化を解析する +  * 指定された対象の水素結合数の経時変化を解析する。 
-  * ''cpptraj'' に同梱のパッケージで解析ができる+  * ''cpptraj'' に同梱のパッケージで解析ができる。
   * パラメータファイル (''test.prmtop'') と座標ファイル (''test.inpcrd'' / ''*.nc'' / ''*.mdcrd'') が必要\\ <code bash>   * パラメータファイル (''test.prmtop'') と座標ファイル (''test.inpcrd'' / ''*.nc'' / ''*.mdcrd'') が必要\\ <code bash>
 $ cpptraj $ cpptraj
行 88: 行 88:
  
 ===== 水和水分析 ===== ===== 水和水分析 =====
-  * 指定された対象の周囲の水分子数をカウントする +  * 指定された対象の周囲の水分子数をカウントする。 
-  * ''cpptraj'' に同梱のパッケージで解析ができる+  * ''cpptraj'' に同梱のパッケージで解析ができる。
   * パラメータファイル (''test.prmtop'') と座標ファイル (''test.inpcrd'' / ''*.nc'' / ''*.mdcrd'') が必要\\ <code bash>   * パラメータファイル (''test.prmtop'') と座標ファイル (''test.inpcrd'' / ''*.nc'' / ''*.mdcrd'') が必要\\ <code bash>
 $ cpptraj $ cpptraj
行 104: 行 104:
  
 ===== 距離解析 ===== ===== 距離解析 =====
-  * 指定された対象間の距離の経時変化を解析する +  * 指定された対象間の距離の経時変化を解析する。 
-  * ''cpptraj'' に同梱のパッケージで解析ができる+  * ''cpptraj'' に同梱のパッケージで解析ができる。
   * パラメータファイル (''test.prmtop'') と座標ファイル (''test.inpcrd'' / ''*.nc'' / ''*.mdcrd'') が必要   * パラメータファイル (''test.prmtop'') と座標ファイル (''test.inpcrd'' / ''*.nc'' / ''*.mdcrd'') が必要
-  * 対象1 ''TARGET1_MASK'' と対象2 ''TARGET2_MASK'' の重心 (中心) 間の距離を測定する\\ <code bash>+  * 対象1 ''TARGET1_MASK'' と対象2 ''TARGET2_MASK'' の重心 (中心) 間の距離を測定する。\\ <code bash>
 $ cpptraj $ cpptraj
 > parm test.prmtop > parm test.prmtop
行 118: 行 118:
     * ''OUTPUT'': 出力ファイル     * ''OUTPUT'': 出力ファイル
     * ''geom'': 幾何学的中心で測定 (指定しなければ、重心で測定)     * ''geom'': 幾何学的中心で測定 (指定しなければ、重心で測定)
 +
 +===== 水分子等を一定距離で切り出した構造の抽出 =====
 +  * 生体分子から一定距離にある水分子のみを残した構造をトラジェクトリから抽出する。\\ <code bash>
 +$ cpptraj
 +> parm test.prmtop
 +> trajin test.inpcrd
 +> mask MASK_BIOMOL<:DIST maskpdb OUTPUT.pdb
 +> go
 +</code>
 +    * ''MASK_BIOMOL'': 生体分子のマスク
 +    * ''DIST'': 生体分子からの距離 (例: '':5.0'' = 残基単位で切り出す; ''@5.0'' = 原子単位で切り出す)
 +    * ''OUTPUT.pdb'': 出力するファイル (''OUTPUT.pdb.1'' のように、ファイル名の後にフレーム番号が付けられ、フレーム毎に出力される)
 +    * カウンターイオンが含まれており、水分子のみを切り出したい場合は、''(:1-20)|((:1-20<:5.0)&WAT)'' のように、論理演算子を使う。
 +    * 距離の算出には trajin のそれぞれのフレームが用いられるが、別のリファレンス構造を用いたい場合は、''reference'' と ''activeref REFERENCE_INDEX'' を指定する。\\ <code bash>
 +$ cpptraj
 +> parm test.prmtop
 +> trajin test.inpcrd
 +> reference test2.inpcrd
 +> activeref 0
 +> mask MASK_BIOMOL<:DIST maskpdb OUTPUT.pdb
 +> go
 +</code>
  
 ===== カーブフィッティング ===== ===== カーブフィッティング =====
  • 分子シミュレーション関連/amber/解析方法.1634277187.txt.gz
  • 最終更新: 2021/10/15 14:53
  • by mumeiyamibito