Autodock问题集锦课件

Autodock问题集锦1.autodock给出的bindingenergy可靠性如何？Autodock给出的结合自由能仅仅是一个打分函数！很粗略的这个仅仅用来评判小分子与受体结合得是否好的一个标准。就如同GOLD或者Surflex等对接软件给出的chemscoreGOLDsocreCscore等等这样的打分函数！是由半经验公式得到的，并不能准确的描述结合的自由能。除了autodock还有以下软件可以计算自由能或打分a.原来的Insight2005（现在的DS）中的docking，可以给出两者相互作用能。包括静电能和范德华能。b.dock也能算，MM/PBSA打分，另外schrodinger旗下的macromodel模块，prime模块都能算bindingenergy2.用autodock做对接，请问里面格子大小及其中心的选择依据格子的中心和你的活性位点有关~一般设在你的活性位点上即可，大小只要能包括你觉得可能发生作用的区域就可以了~3.AUTODOCK产生的氢键原理在AUTODOCK4.0做对接，对接完成后看蛋白与配体相互作用的时候，有一个氢键形成，是蛋白中的一个蛋氨酸上的-S-与配体中的-O-形成氢键该O上没有氢，两边连的是C查看TXT文本，该S原子的类型为SD，也就是S是氢键给体AUTODOCK里定义的是S是氢键给体，O是氢键受体。对于这个H键，显而易见，O是氢键受体，但是S是如何成为氢键给体的呢？它的两个键也是接的是C求助解释下，S与O是如何形成氢键的？如果两个杂原子间距（小于两个氢键键长）和角度在一定范围内，很多软件会认为它们能通过H原子介导形成氢键相互作用，即共用一个氢原子，各软件定义的范围不同，此氢原子亦可来自于水分子，当然你对接过程中可能没有添加水分子，实际上分子间空隙是会填充水的另外建议你用Sybyl看氢键会好一些4.使用autodock之前需要进行能量最小化优化吗？对接时，计算机会自动寻找小分子不同的构象进行对接，但寻找的方式有随机性，在你设置对接次数大时可能影响不大，但如果设置的次数少就可能没找到最低能构象，如果对接之前能量最小化了，不管设置的参数怎样，最低能构象一定在里面5.关于AutoDock输出结果中的能量值下面是autodock运行打分值最好的输出结果Rank:1_1BindingEnergy:-15.1kI:8.58pMIntermolecularEnergy:-14.39InternalEnergy:-5.38TorsionalEnergy:3.84UnboundExtendedEnergy:-0.83ClusterRMS:0.0RefRMS:6.79tutorial中说DockEnergy=intermolecularenergy+internalenergy那么在此例中DockEnergy=-14.39+(-5.38)=-19.77,而BindingEnergy=-15.1,两者并不相等,但两者是平行的,按这两个能量值进行排序都得到相同的排序结果我很困惑的是如果BindingEnergy就是我们关心的BindingFreeEnergy(结合自由能△G),那么DockEnergy又是什么含义?DockEnergyandBindingEnergy有什么不同,他们都是怎么定义的?还有这里BindingEnergy的单位是kcal/mol还是kJ/mol或是其它?关于dockedEnergy:对接能量＝分子间能量+分子内能量关于freeEnergy:自由能＝分子间能量项+扭转熵因子罚分项（这一项的值是一个常数乘以配体内可旋转键数目，此处我怀疑这一项的值要么是那个unboundextendedenergy,是这一项值的可能性比较大，因为这两项加起来很接近freeEnergy的值，但又不是绝对相等，小数点后总差点。）分子间的结合自由能计算不应包括分子内能,这点我非常的认同.我的理解是文献中引用的能量值△G(结合自由能,bindingfreeenergy)就是输出中的bindingenergy,即△G=bindingfreeenergy=bindingenergy,AutoDocK结果中的△G应该直接计算得出,不可能再经过计算得到.freeEnergy与Bindingenergy虽然很接近,但不绝对相等.很可能是它们的定义略有不同,而不能将他们等同起来,因为我觉得AutoDocK那么广泛使用的程序,不可能在这点上没有做到精确.BindingEnergy=IntermolecularEnergy+InternalEnergy+TorsionalEnergy-UnboundExtendedEnergy,单位是kcal/mol。internalenergy可能是指配体结合前后内能的差值，不是内能的绝对值。AutoDocK中评价函数中不应该包含内能(差值),设想若对接的起始构象不同,理论上如果对接好的化,它们会得到同一构象和取向,那么这里internalenergy(差值)是不同的,显然不能包含在AutoDocK中评价函数中我查阅了徐莜杰的书及相关资料以及tutarial中的解释,得出一下结论△G(bind)(评价函数)=△G(vdw)+△G(H-bond)+△G(ele)+△G(del-sol)+△G(tor)经过计算证实前四项包含在intermolecularenergy中,另外,freeenergy=intermolecularenergy+torsionentropypenalty正如2楼兄弟所说freeEnergy:自由能＝分子间能量项+扭转熵因子罚分项（这一项的值是一个常数乘以配体内可旋转键数目，这一项的值是那个unboundextendedenergy,可能性比较大，因为这两项加起来很接近freeEnergy的值，但又不是绝对相等，小数点后总差点。）对接结果中TorsionalEnergy:3.84UnboundExtendedEnergy:-0.83这两个值是固定的,即不论你排名如何,构象如何,他们的值都是这个.这两个能量表示的是同一个东西.所以推测UnboundExtendedEnergy=△G(tor)这样△G(bind)(评价函数)=△G(vdw)+△G(H-bond)+△G(ele)+△G(del-sol)+△G(tor)=intemolecarenergy+△G(tor)=freeenergy约=bindingenergy而dockedenergy在这里对我们来说没有多少价值6.Analyze方法(八)/*Analyze*/A：1.Analyze-Docking-(檔名.dlg)2.Analyze-Conformations-Play(可讓ligand變更位置，ps.已對結好的位置)3.Macromolecule-(檔名.pdbqt)(protein)(八)/*Analyze*/B：1.Analyze-Grids-OpenOther[選擇步驟(五)中的九個由autogrid產生的文件，每一個文件都代表不同的力場]2.Analyze-Conformations-Play(ligand可配合任一力場，做autodock所產生出的九個docking位置，並判斷它們的關係)(八)/*Analyze*/C：1.Analyze-Conformations-ExtractHistogram計⋯算(分子對接，ligand與proteinbindsite的分數高低並由高到