±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 1814  |  »Ø¸´: 2

MolAICal

гæ (³õÈëÎÄ̳)

[½»Á÷] ʹÓÃMolAICal»ùÓÚNAMDÄ£Äâ½á¹û¼ÆËãС·Ö×Ӻ͵°°×MM/GBSAµÄ½Ì³Ì ÒÑÓÐ1È˲ÎÓë

ʹÓÃMolAICal»ùÓÚNAMDÄ£Äâ½á¹û¼ÆËãС·Ö×Ӻ͵°°×MM/GBSAµÄ½Ì³Ì


¸ü¶à½Ì³Ì£¨º¬Ó¢ÎĽ̳̣©Çë¼ûÈçÏ£º
MolAICal¹Ù·½Ö÷Ò³£ºhttps://molaical.github.io
MolAICal ÎÄÕ½éÉÜ£ºhttps://doi.org/10.1093/bib/bbaa161
MolAICalÖÐÎIJ©¿Í£ºhttps://molaical.github.io/cntutorial.html
MolAICal blogspot£ºhttps://qblab.blogspot.com


1. ¼ò½é
ÔÚ±¾½Ì³ÌÖнéÉÜÁË»ùÓÚNAMDµÄ·Ö×Ó¶¯Á¦Ñ§Ä£Äâ½á¹û£¬Ê¹ÓÃMolAICal¼ÆËãС·Ö×ÓºÍMproµ°°×ÊÜÌåMM/GBSAµÄ·½·¨¡£±¾½Ì³ÌÖ»ÊÇÒ»¸ö¼òµ¥ÑÝʾ¡£ÎªÁ˽ÚÊ¡ÔËÐм°´æ´¢¿Õ¼ä£¬±¾½Ì³Ì½öÑ¡ÔñÁËMpro¸´ºÏÎï·Ö×Ó¶¯Á¦Ñ§Ä£ÄâµÄ25Ö¡ÓÃÓÚ¼ÆËã¡£

2.¹¤¾ß
2.1. ËùÐèÈí¼þÏÂÔØµØÖ·
1) MolAICal : https://molaical.github.io
2) NAMD: https://www.ks.uiuc.edu/Research/namd/

2.2. ²Ù×÷ʾÀýÎļþ
ËùÓÐÓõ½µÄ²Ù×÷½Ì³ÌÎļþ¾ù¿ÉÔÚÏÂÃæµÄÍøÕ¾ÏÂÔØ£º
https://github.com/MolAICal/tutorials/tree/master/004-MMGBSA

3. ²Ù×÷Á÷³Ì
תµ½ÒÔÏÂĿ¼£º
#> cd 004-MMGBSA

3.1. ÌáÈ¡µ°°×ÓëÅäÌ帴ºÏÎïµÄ¹ì¼£Îļþ
CODE:
#> vmd -dispdev text -psf "mpro.psf" -e stripDCD.vmd -args protein,or,resname,LIG "mpro.dcd"  "complex" mpro.psf mpro.pdb

-args: ÆäÓ÷¨ÀàËÆVMDÈí¼þÖеġ°atomselect¡±ÃüÁ±ÈÈç"atomselect top protein or resname LIG"£¬´Ë´¦¶ººÅ"," ´ú±í¿Õ¸ñ¡£ÆäÖнű¾ÎļþstripDCD.vmd¿ÉÒÔÔÚ±¾½Ì³Ì²ÄÁÏ»òMolAICalÈí¼þµÄ¡°scripts¡±Ä¿Â¼ÀïÃæÕÒµ½¡£

Ö´ÐÐÉÏÊöÃüÁîºóÉú³Écomplex.psf, complex.pdb ºÍ complex.dcdÎļþ¡£½«¡°GBIS¡± ºÍ¡°sasa¡±²ÎÊýÉèÖÃΪon¡£´ò¿ª²¢°´ÕÕÏÂÎÄÄÚÈÝÐ޸ġ°complex.conf¡±Îļþ£º
--------------------------------------------------------------------------
structure          complex.psf
coordinates        complex.pdb
outputName         complex

paraTypeCharmm      on
parameters          par_all36_prot.prm
parameters          par_all36_cgenff.prm
parameters          ligand.str
parameters          toppar_water_ions.str

coorfile open dcd complex.dcd
--------------------------------------------------------------------------
±¾½Ì³ÌÖÐÃüÁîÔÚCPUÉÏÔËÐС£Äã¿ÉÒÔÑ¡ÔñGPU½øÐÐÔËËã¡£ÔÚLinuxϵͳÏÂÔËÐÐNAMDÃüÁÈçÏ£º
CODE:
#> namd2 +p3 complex.conf >& complex.log &

ÆäÖзûºÅ¡°&¡±´ú±í³ÌÐòÔÚLinuxϵͳÖнøÐкǫ́ÔËÐУ¬Èç¹ûÄãʹÓõÄÊÇWindows²Ù×÷ϵͳ£¬Çë²»ÒªÓá°&¡±£¬ÀýÈ磬ÃüÁî»»³ÉÕâÑù£º
CODE:
#> namd2 +p3 complex.conf > complex.log

3.2. ½öÌáÈ¡µ°°×µÄ¹ì¼£Îļþ£º
CODE:
#> vmd -dispdev text -psf "mpro.psf" -e stripDCD.vmd -args protein "mpro.dcd" "protein" mpro.psf mpro.pdb

ÉÏÊöÃüÁî»áÉú³É protein.psf, protein.pdbºÍprotein.dcd¡£´ò¿ª¡°protein.conf¡±£¬²Î¿¼ ¡°complex.conf¡±ÐÞ¸ÄÏà¹Ø²ÎÊý¡£

±¾½Ì³ÌÖÐÃüÁîÔÚCPUÉÏÔËÐС£Äã¿ÉÒÔÑ¡ÔñGPU½øÐÐÔËËã¡£ÔÚLinuxϵͳÏÂÔËÐÐNAMDÃüÁÈçÏ£º
CODE:
#> namd2 +p3 protein.conf >& protein.log &

3.3. ½öÌáÈ¡ÅäÌåµÄ¹ì¼£Îļþ£º
CODE:
#> vmd -dispdev text -psf "mpro.psf" -e stripDCD.vmd -args resname,LIG "mpro.dcd" "ligand" mpro.psf mpro.pdb

ÉÏÊöÃüÁî»áÉú³É ligand.psf, ligand.pdbºÍligand.dcd¡£´ò¿ª¡°ligand.conf¡±£¬²Î¿¼ ¡°complex.conf¡±ÐÞ¸ÄÏà¹Ø²ÎÊý¡£

±¾½Ì³ÌÖÐÃüÁîÔÚCPUÉÏÔËÐС£Äã¿ÉÒÔÑ¡ÔñGPU½øÐÐÔËËã¡£ÔÚLinuxϵͳÏÂÔËÐÐNAMDÃüÁÈçÏ£º
CODE:
#> namd2 +p3 ligand.conf >& ligand.log &   

4. ÓÃMolAICal¼ÆËãMM/GBSA
CODE:
#> molaical.exe -mmgbsa -c complex.log -r protein.log -l ligand.log

Êä³ö½á¹ûÖиø³öÏÂÎÄËùʾµÄ½áºÏ×ÔÓÉÄÜ¡÷G£º
--------------------------------------------------------------------------
delta E(internal): -4.0000007572871255E-6
delta E(electrostatic) + deltaG(sol): 7.702936000001536
delta E(VDW) + deltaG(sol): -44.43611599999989
delta G binding: -36.73318399999911
--------------------------------------------------------------------------
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹ØÉ̼ÒÍÆ¼ö: (ÎÒÒ²ÒªÔÚÕâÀïÍÆ¹ã)

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

zxhuizhang

ľ³æ (ÕýʽдÊÖ)

¡ï
Сľ³æ: ½ð±Ò+0.5, ¸ø¸öºì°ü£¬Ð»Ð»»ØÌû
ÎÒÏëÎÊһϣ¬¹ØÓÚlinuxϵͳµÄµçÄÔÅäÖÃÎÊÌâ¡£´ó¸ÅÐèҪʲôÑù×ÓµÄ
¼û´Ë
2Â¥2021-02-02 12:12:09
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

MolAICal

гæ (³õÈëÎÄ̳)

ÒýÓûØÌû:
2Â¥: originally posted by zxhuizhang at 2021-02-02 12:12:09
ÎÒÏëÎÊһϣ¬¹ØÓÚlinuxϵͳµÄµçÄÔÅäÖÃÎÊÌâ¡£´ó¸ÅÐèҪʲôÑù×ÓµÄ

Èç¹ûÒѾ­ÅÜÍêmd, ½ö½öÓÃmolaicalËã, Ò»°ãÆÕͨ¼ÒÓõÄlinuxϵͳ¾ÍÐС£¡£µ±È»Ó÷þÎñÆ÷¸üºÃÁË
3Â¥2021-04-27 12:26:15
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ MolAICal µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] ÉúÎïѧ308·ÖÇóµ÷¼Á£¨Ò»Ö¾Ô¸»ª¶«Ê¦´ó£© +6 ÏàÐűػá¹ââÍòÕ 2026-04-06 7/350 2026-04-06 11:51 by lijunpoly
[Óлú½»Á÷] 3-äåßÁà¤-4-¼×È©ºÏ³É 10+3 Ñ©¸âÒ²Äܱ¬Õ¨ 2026-03-31 9/450 2026-04-06 00:24 by nBuï®
[¿¼ÑÐ] 304Çóµ÷¼Á +7 c297914 2026-04-05 8/400 2026-04-05 22:13 by hemengdong
[¿¼ÑÐ] 296Çóµ÷¼Á +3 Íô£¡£¿£¡ 2026-04-05 4/200 2026-04-05 20:13 by à£à£à£0119
[¿¼ÑÐ] 288Çóµ÷¼Á +7 ûÓдð°¸_ 2026-04-05 7/350 2026-04-05 17:01 by yulian1987
[¿¼ÑÐ] 304Çóµ÷¼Á£¨085602£¬¹ýËļ¶£¬Ò»Ö¾Ô¸985£© +15 »¯¹¤ÈË999 2026-04-04 15/750 2026-04-05 16:25 by ÎÒÊǵç·çÉÈr
[¿¼ÑÐ] 344²ÄÁÏÓ뻯¹¤µ÷¼Á +9 µ÷¼ÁÉϰ¶«^ 2026-04-03 9/450 2026-04-04 23:10 by happyddm
[¿¼ÑÐ] 321Çóµ÷¼Á +6 ÈÏÕæÇóÉÏѧ 2026-04-03 6/300 2026-04-04 19:51 by dongzh2009
[¿¼ÑÐ] 26µ÷¼Á 086003 +6 ʧ»îµÄϸ°û 2026-04-04 6/300 2026-04-04 09:50 by zhangdingwa
[¿¼ÑÐ] 311Çóµ÷¼Á +20 zchqwer 2026-04-01 22/1100 2026-04-03 22:09 by lglzsd
[¿¼ÑÐ] 288Çóµ÷¼Á Ò»Ö¾Ô¸¹þ¹¤´ó ²ÄÁÏÓ뻯¹¤ +39 ÂåÉñ¸ç¸ç 2026-03-31 41/2050 2026-04-03 21:51 by qlm5820
[¿¼ÑÐ] 285Çóµ÷¼Á +5 AZMK 2026-04-03 8/400 2026-04-03 18:17 by AZMK
[¿¼ÑÐ] ÍÁľˮÀû328·ÖÇóµ÷¼Á +6 ¼²·çÖª¾¢²Ý666 2026-04-02 6/300 2026-04-03 11:38 by znian
[¿¼ÑÐ] 325·Ö»¯Ñ§µ÷¼Á +5 15771691647 2026-04-02 5/250 2026-04-03 09:58 by ChemPharm
[¿¼ÑÐ] Ò»Ö¾Ô¸ÏÃÃÅ´óѧ²ÄÁϹ¤³Ìר˶354ÕÒµ÷¼Á£¡£¡£¡ +8 ±´ß±µ±µ 2026-03-30 8/400 2026-04-03 09:41 by hypershenger
[¿¼ÑÐ] 346Çóµ÷¼Á +5 Ö£³ÏÀÖ 2026-04-02 5/250 2026-04-02 16:38 by SZW_UJN
[¿¼ÑÐ] 298ÇóBÇøµ÷¼Á +4 zzz£¬£¬r 2026-04-02 5/250 2026-04-02 12:17 by ÍÁľ˶ʿÕÐÉú
[¿¼ÑÐ] Ò»Ö¾Ô¸±±¾©¿Æ¼¼´óѧ085601²ÄÁϹ¤³ÌÓ¢Ò»Êý¶þ³õÊÔ×Ü·Ö335Çóµ÷¼Á +5 Ë«ÂíβƦÀϰå2 2026-03-31 5/250 2026-04-01 09:04 by oooqiao
[¿¼ÑÐ] ºÏ·ÊÇøÓòÐÔÖØµãÒ»±¾ÕÐÊÕµ÷¼Á +4 6266jl 2026-03-30 8/400 2026-03-31 18:43 by 6266jl
[¿¼ÑÐ] Ò»Ö¾Ô¸Î÷µç085401ÊýÒ»Ó¢Ò»299Çóµ÷¼Á Áù¼¶521 +4 °®³Ô´óѼÀæ 2026-03-31 4/200 2026-03-31 11:51 by ²«»÷518
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û