Membrane Protein Simulation with Desmond
来源:互联网 发布:宏投网络 编辑:程序博客网 时间:2024/04/29 00:12
Desmond supports algorithms typically used to perform fast and accurate molecular dynamics. Long-range electrostatic energy and forces can be calculated using particle-mesh-based Ewald techniques.Constraints can be enforced using the M-SHAKE algorithm. These approaches can be used in combination with time-scale splitting (RESPA-based) integration schemes. Here are the generally steps for membrane protein simulation in Desmond.
1. System Building
This can be done very easily with Maestro GUI from Schrodinger.
menu-->application-->desmond-->system builder
this will need 7 steps to be finished as shown in the snapshot. The membrane in Desmond are all pre-equilibrated which is quite important for initial membrane system. step 2 in the system building step is also quite important, for we can place the membrane in pre-aligned structures so that our protein can be embedded in correct origination in membrane.
The default size for membrane system is 10X10X10 A^3. It is necessary to adjust this according to the according to our membrane protein size. Sometimes, membrane system may be broken if the system is too small.
step6 can add salt to the system, namely 0.15 M NaCl usually,to mimic physiology environment.
Output files:
desmond_setup-out.cms
2. Assign Force Field
We can only assign OPLS_FF in the initial step. However, most of us would like to use other force field such as CHARMM or AMBER. In the latest version of Desmond, the following FFs are available:
$SCHRODINGER/run -FROM desmond viparr.py [options] input-file output-file
for example:
$SCHRODINGER/run -FROM desmond viparr.py -f charmm27 -f tip3p desmond_setup-out.cms charmm.cms
Viparr uses atomic numbers and bond structure to match residues to templates. Thus if there are non-standard atom or residue PDB names in the .cms file, there is no need to modify them to match the names used in the force field.
You may also use multiple force fields when parameters assigned to a particular residue by one force field override parameters assigned by another force field. This situation is similar to combining force fields; residues in the chemical system match templates in more than one of the specified force fields and all matching force fields are applied. However, in this case two or more force fields provide parameters for the same term; for example, two force fields provide parameters for the angle between atoms 1, 2, and 3 causing a conflict. The conflict is resolved by allowing the first force field(that matches the residue) to take precedence. The order is the order in which the force fields were
specified on the command line.
3. Add Constraints
While the System Builder automatically adds all the necessary constraints to the .cms file, Viparr does not. The build_constraints.py script is all but mandatory; otherwise Desmond cannot make use of the SHAKE algorithm to constrain bonds between heavy atoms and hydrogens. The build_constraints.py script should be run on your .cms file after every use of viparr. Viparr removes the constraints each time.
If you process a .mae or .cms file with Viparr you will have to add the constraints in a separate step by running the following command:
$SCHRODINGER/run -FROM desmond build_constraints.py input.cms constraints.cms
4. Relax System
It can be difficult to relax freshly built protein-membrane systems. In particular, penetration of the water between the protein than the lipids can be problematic and require very lengthy simulations to correct. A relaxation protocol that should reduce or eliminate such problems is available by running the script relax_membrane.py from the command line.
a. Run the script to prepare the necessary input files:
$SCHRODINGER/run relax_membrane.py -i constraints.cms -t temperature -j protein-membrane
b. Run the membrane relaxation protocol using the command:
$SCHRODINGER/utilities/multisim -JOBNAME protein-membrane -HOST myhost -mode umbrella -cpu cpus -i protein-membrane-in.cms -m protein-membrane.msj -o protein-membrane-out.cms
5. MD Production
Now we get the relaxed membrane system protein-membrane-out.cms at hand, and MD can be produced either by importing it to Maestro or by command line.
$SCHRODINGER/utilities/multisim -JOBNAME myMD -HOST my-16cpu -maxjob 0 -cpu "2 2 4" -m desmond_md_job.msj -c desmond_md_job.cfg -description "Molecular dynamics"-mode umbrella -iprotein-membrane-out.cms -o desmond_md_job-out.cms -DISP append > log &
Tips:
We can setup host, CPU, tmp directory from /home/user/.schrodigner/schrodinger.hosts. Command line based jobs can be monitored from the tmp directory which we set in schrodinger.hosts file.
Reference
Desmond User Manual
Desmond Tutorial
Desmond User's Guide
download:
http://www.deshawresearch.com/downloads/download_desmond.cgi/
resource:
http://www.deshawresearch.com/resources_desmond.html
user group:
http://groups.google.com/group/desmond-md-users?hl=en
helps:
Desmond@DEShawResearch.com
http://blog.sciencenet.cn/blog-355217-458994.html
1. System Building
This can be done very easily with Maestro GUI from Schrodinger.
menu-->application-->desmond-->system builder
this will need 7 steps to be finished as shown in the snapshot. The membrane in Desmond are all pre-equilibrated which is quite important for initial membrane system. step 2 in the system building step is also quite important, for we can place the membrane in pre-aligned structures so that our protein can be embedded in correct origination in membrane.
The default size for membrane system is 10X10X10 A^3. It is necessary to adjust this according to the according to our membrane protein size. Sometimes, membrane system may be broken if the system is too small.
step6 can add salt to the system, namely 0.15 M NaCl usually,to mimic physiology environment.
Output files:
desmond_setup-out.cms
2. Assign Force Field
We can only assign OPLS_FF in the initial step. However, most of us would like to use other force field such as CHARMM or AMBER. In the latest version of Desmond, the following FFs are available:
charmm22nocmap
charmm27
charmm32
charmm36_lipids
amber94
amber96
amber99
amber99SB
amber99SB-ILDN
amber03
oplsaa_impact_2001
oplsaa_impact_2005
oplsaa_ions_Jensen_2005
tip3p
spce
tip3p_charmm
tip4p
tip4pew
tip5p
FF can be assigned by command:charmm27
charmm32
charmm36_lipids
amber94
amber96
amber99
amber99SB
amber99SB-ILDN
amber03
oplsaa_impact_2001
oplsaa_impact_2005
oplsaa_ions_Jensen_2005
tip3p
spce
tip3p_charmm
tip4p
tip4pew
tip5p
$SCHRODINGER/run -FROM desmond viparr.py [options] input-file output-file
for example:
$SCHRODINGER/run -FROM desmond viparr.py -f charmm27 -f tip3p desmond_setup-out.cms charmm.cms
Viparr uses atomic numbers and bond structure to match residues to templates. Thus if there are non-standard atom or residue PDB names in the .cms file, there is no need to modify them to match the names used in the force field.
You may also use multiple force fields when parameters assigned to a particular residue by one force field override parameters assigned by another force field. This situation is similar to combining force fields; residues in the chemical system match templates in more than one of the specified force fields and all matching force fields are applied. However, in this case two or more force fields provide parameters for the same term; for example, two force fields provide parameters for the angle between atoms 1, 2, and 3 causing a conflict. The conflict is resolved by allowing the first force field(that matches the residue) to take precedence. The order is the order in which the force fields were
specified on the command line.
3. Add Constraints
While the System Builder automatically adds all the necessary constraints to the .cms file, Viparr does not. The build_constraints.py script is all but mandatory; otherwise Desmond cannot make use of the SHAKE algorithm to constrain bonds between heavy atoms and hydrogens. The build_constraints.py script should be run on your .cms file after every use of viparr. Viparr removes the constraints each time.
If you process a .mae or .cms file with Viparr you will have to add the constraints in a separate step by running the following command:
$SCHRODINGER/run -FROM desmond build_constraints.py input.cms constraints.cms
4. Relax System
It can be difficult to relax freshly built protein-membrane systems. In particular, penetration of the water between the protein than the lipids can be problematic and require very lengthy simulations to correct. A relaxation protocol that should reduce or eliminate such problems is available by running the script relax_membrane.py from the command line.
a. Run the script to prepare the necessary input files:
$SCHRODINGER/run relax_membrane.py -i constraints.cms -t temperature -j protein-membrane
b. Run the membrane relaxation protocol using the command:
$SCHRODINGER/utilities/multisim -JOBNAME protein-membrane -HOST myhost -mode umbrella -cpu cpus -i protein-membrane-in.cms -m protein-membrane.msj -o protein-membrane-out.cms
5. MD Production
Now we get the relaxed membrane system protein-membrane-out.cms at hand, and MD can be produced either by importing it to Maestro or by command line.
$SCHRODINGER/utilities/multisim -JOBNAME myMD -HOST my-16cpu -maxjob 0 -cpu "2 2 4" -m desmond_md_job.msj -c desmond_md_job.cfg -description "Molecular dynamics"-mode umbrella -iprotein-membrane-out.cms -o desmond_md_job-out.cms -DISP append > log &
Tips:
We can setup host, CPU, tmp directory from /home/user/.schrodigner/schrodinger.hosts. Command line based jobs can be monitored from the tmp directory which we set in schrodinger.hosts file.
Reference
Desmond User Manual
Desmond Tutorial
Desmond User's Guide
download:
http://www.deshawresearch.com/downloads/download_desmond.cgi/
resource:
http://www.deshawresearch.com/resources_desmond.html
user group:
http://groups.google.com/group/desmond-md-users?hl=en
helps:
Desmond@DEShawResearch.com
http://blog.sciencenet.cn/blog-355217-458994.html
- Membrane Protein Simulation with Desmond
- Business Process Modeling and Simulation with UML
- Board Simulation with SO-DIMM slot
- learning R with swirl-simulation(模拟)
- Modelling the Wireless Propagation Channel: A simulation approach with Matlab
- omnet找不到libinet.dll+Simulation terminated with exit code : -1073741515
- omnet++ Simulation terminated with exit code: -1073741819 的一个解决方案
- Co-Simulation of C with Verilog in VCS
- RH-AS5下安装Desmond
- NCBI推出blastp加速服务(Accelerated protein-protein BLAST)
- In Silico: 3D Animation and Simulation of Cell Biology with Maya and MEL
- Behavioral Simulation,Structural Simulation,Timing Simulation
- Rosetta Protein Prediction Tools-----clean_pdb.py
- Cloth Simulation
- Simulation Robot
- Game Simulation
- Rigidbody Simulation
- game simulation
- 2012年度最佳 Web 前端开发工具和框架——《下篇》
- vxworks下MUX接口的使用
- 计算广告学习笔记3.5 受众定向-数据加工和交易
- 拯救大学生
- 安卓手机怎么连接电脑上网,如何设置图文教程详解
- Membrane Protein Simulation with Desmond
- iOS系统GCD学习(8):Dispatch Queue挂起恢复和目标指定
- junit走起~~~
- Datashee常用名词
- VxWorks网络协议栈的MUX接口
- MyEclipse8.6 安装aptana插件
- java通过poi解析Excel文件,直接读取客户端文件
- 利用jdbc实现批处理
- C++ 和 汇编 混合编程幼儿园经验总结**持续更新** 学习笔记贴