| 3 | 1/1 | 返回列表 |
| 查看: 1050 | 回復(fù): 2 | ||
cfmzxf84木蟲 (正式寫手)
|
[求助]
求助AMBER的能量分解和氫鍵分析in文件
|
|
求助,做小分子和蛋白的能量分解和氫鍵分析,看了官網(wǎng)的教程,但是一直做不對。求知道的同學(xué)發(fā)兩個in文件學(xué)習(xí)下參數(shù)。 主要目的:小分子與殘基結(jié)合能分解;小分子與殘基氫鍵分析 謝謝 |
金蟲 (小有名氣)
|
MM-PBSA的 input file &general endframe=500, interval=1, verbose=1 / #&gb # igb=5, saltcon=0.100 #/ &pb istrng=0 / &decomp idecomp=2, dec_verbose=3, print_res='1-15,36-44,64,88-95,158-161,176-199' (這個是你要分解的殘基) ptraj氫鍵掃描的腳本: ############################################################################# #--This H-bonding analysis script is created based on one script available #--on the wed targeting h-bond formation of protein and DNA, first, whose #--original writer should be appreciated ! #-- #--This enhanced script modified by me is able to capture all the donors & #--acceptors of standard residues from protein-protein, protein-DNA as well #--as protein and RNA in principle. #-- #--Usage: #--For protein-RNA: Actually, this script is exactly applicable to pro-RNA! #--For protein-DNA: Please comment corresponding lines related to 2'OH in # RNA, also,remember to substitute DG for DG (for instance) #--For protein-protein: Just need to comment all the line regarding DNA/RNA. #-- #--Some issues are inevitable even though I have checked it many times.Please #--feel free to contact me in either following ways for useful discussion: #--E-mail: vip_zzc@emails.bjut.edu.cn #-- QQ: 350 2826 89 #-- Written by Zhicheng Zuo, 5/3/2011 ############################################################################# ######### laod trajectory ############ ###################################### trajin 1abc.mdcrd/mdcrd.gz/binpos trajin etc. trajin etc. ###################################### #-- Donors from amino acids sidechains ###################################### donor mask :GLN@OE1 donor mask :GLN@NE2 donor mask :ASN@OD1 donor mask :ASN@ND2 donor mask :TYR@OH donor mask :ASP@OD1 donor mask :ASP@OD2 donor mask :GLU@OE1 donor mask :GLU@OE2 donor mask :SER@OG donor mask :THR@OG1 donor mask :HIS@ND1 donor mask :HIE@ND1 donor mask :HID@NE2 ###################################### #-- Donors from nucleotides: DNA/RNA ###################################### ###### DC/DC ####### donor mask :DC@P donor mask :DC@O1P donor mask :DC@O2P donor mask :DC@O5' donor mask :DC@O4' donor mask :DC@N1 donor mask :DC@N4 donor mask :DC@N3 donor mask :DC@O2 donor mask :DC@O2' donor mask :DC@O3' ###### DG/DG ####### donor mask :DG@P donor mask :DG@O1P donor mask :DG@O2P donor mask :DG@O5' donor mask :DG@O4' donor mask :DG@N9 donor mask :DG@N7 donor mask :DG@O6 donor mask :DG@N1 donor mask :DG@N2 donor mask :DG@N3 donor mask :DG@O2' donor mask :DG@O3' ###### DU/DT ####### donor mask :DU@P donor mask :DU@O1P donor mask :DU@O2P donor mask :DU@O5' donor mask :DU@O4' donor mask :DU@N1 donor mask :DU@O4 donor mask :DU@N3 donor mask :DU@O2 donor mask :DU@O2' donor mask :DU@O3' ###### DA/DA ####### donor mask :DA@P donor mask :DA@O1P donor mask :DA@O2P donor mask :DA@O5' donor mask :DA@O4' donor mask :DA@N9 donor mask :DA@N7 donor mask :DA@N6 donor mask :DA@N1 donor mask :DA@N3 donor mask :DA@O2' donor mask :DA@O3' ####################################### #--Acceptors from animo acid sidechains ####################################### acceptor mask :ASN@ND2 :ASN@HD21 acceptor mask :ASN@ND2 :ASN@HD22 acceptor mask :TYR@OH :TYR@HH acceptor mask :GLN@NE2 :GLN@HE21 acceptor mask :GLN@NE2 :GLN@HE22 acceptor mask :TRP@NE1 :TRP@HE1 acceptor mask :LYS@NZ :LYS@HZ1 acceptor mask :LYS@NZ :LYS@HZ2 acceptor mask :LYS@NZ :LYS@HZ3 acceptor mask :SER@OG :SER@HG acceptor mask :THR@OG1 :THR@HG1 acceptor mask :ADG@NH2 :ADG@HH21 acceptor mask :ADG@NH2 :ADG@HH22 acceptor mask :ADG@NH1 :ADG@HH11 acceptor mask :ADG@NH1 :ADG@HH12 acceptor mask :ADG@NE :ADG@HE acceptor mask :HIS@NE2 :HIS@HE2 acceptor mask :HIE@NE2 :HIE@HE2 acceptor mask :HID@ND1 :HID@HD1 acceptor mask :HIP@ND1,NE2 :HIP@HE2,HD1 ####################################### #--Acceptors from nucleotides ####################################### acceptor mask :DC@N4 :DC@H41 acceptor mask :DC@N4 :DC@H42 acceptor mask :DC@O2' :DC@HO'2 acceptor mask :DG@N1 :DG@H1 acceptor mask :DG@N2 :DG@H21 acceptor mask :DG@N2 :DG@H22 acceptor mask :DG@O2' :DG@HO'2 acceptor mask :DU@N3 :DU@H3 acceptor mask :DU@O2' :DU@HO'2 acceptor mask :DA@N6 :DA@H61 acceptor mask :DA@N6 :DA@H62 acceptor mask :DA@O2' :DA@HO'2 ############################################################################################# ##-- NOTE: Pay particular attention to hydrogen associoted with 2'-OH in RNA, the "mask name" ##-- must always be HO'2 other than HO2' or others , otherwise, it cannot be recognized. ############################################################################################# ########################################################### #-- Backbone donors and acceptors for amino acid sidechains # N-H for prolines do not exist so are not in the mask ########################################################### donor mask @O ############################################# #-Terminal residues have different atom names ############################################# donor mask @OXT acceptor mask :1@N :1@H1 acceptor mask :1@N :1@H2 acceptor mask :1@N :1@H3 ############################################################################################ #--Specify distance and angle cutoff of H-bond, in this case, set to 3.5 and 120 repectively. #--To print out information related to time series, keyword "series" is to be included. ############################################################################################ hbond distance 3.5 angle 120 print .05 series hbt ############################################################################################ #--In addition, water or/and ions interacting with a particular donor or acceptor can be also #--keeped track of. Refer to documentation of hydrogen bonding facility in AmberTools manual! ############################################################################################ hbond distance 3.5 angle 120 solventneighbor 6 solventdonor WAT O solventacceptor WAT O H1 solventacceptor WAT O H2 hbond distance 3.5 angle 120 solventneighbor 6 solventdonor Na+ Na+ solventacceptor Na+ Na+ Na+ hbond distance 3.5 angle 120 solventneighbor 6 solventdonor Cl- Cl- solventacceptor Cl- Cl- Cl- |
金蟲 (小有名氣)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +6 | PolarBear11 2026-03-26 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
[考研] 0856材料化工調(diào)劑 總分330 +8 | zhubinhao 2026-03-27 | 8/400 |
|
|
[考研] 322求調(diào)劑 +3 | 舊吢 2026-03-24 | 3/150 |
|
|
[考研] 求調(diào)劑323材料與化工 +7 | 1124361 2026-03-24 | 7/350 |
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap- 2026-03-22 | 10/500 |
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
[考研] 318求調(diào)劑 +3 | plum李子 2026-03-23 | 3/150 |
|
|
[考研] 化學(xué)調(diào)劑 +6 | yzysaa 2026-03-21 | 6/300 |
|
|
[考研] 300分,材料,求調(diào)劑,英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
[考研] 292求調(diào)劑 +4 | 鵝鵝鵝額額額額?/a> 2026-03-24 | 4/200 |
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-21 | 6/300 |
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 6/300 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境,總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|
|
[考研] 考研調(diào)劑 +3 | 呼呼?~+123456 2026-03-21 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
[考研] A區(qū)線材料學(xué)調(diào)劑 +5 | 周周無極 2026-03-20 | 5/250 |
|