4.9 md-write

Initial structures for molecular dynamics simulation

比较活性位点对接结果和D3DOCK对接结果
  • 75.2 /home/databank/zyzhou/project/hsbd/known-pkt
  • (base) dddc@SR665-V3:/home/databank/zyzhou/project/hsbd/known-pkt$ grep -A 1 min */sys-pre/ligand/dock/*
  • 3cl/sys-pre/ligand/dock/smina_56083-03-5-dock-3CL.sdf:> <minimizedAffinity>
  • 3cl/sys-pre/ligand/dock/smina_56083-03-5-dock-3CL.sdf--7.67074
  • --
  • 3cl/sys-pre/ligand/dock/smina_59870-65-4-dock-3CL.sdf:> <minimizedAffinity>
  • 3cl/sys-pre/ligand/dock/smina_59870-65-4-dock-3CL.sdf--7.91815
  • --
  • pl/sys-pre/ligand/dock/smina_1486-70-0-dock-PL.sdf:> <minimizedAffinity>
  • pl/sys-pre/ligand/dock/smina_1486-70-0-dock-PL.sdf--7.90081
  • --
  • pl/sys-pre/ligand/dock/smina_56083-03-5-dock-PL.sdf:> <minimizedAffinity>
  • pl/sys-pre/ligand/dock/smina_56083-03-5-dock-PL.sdf--9.61422
  • --
  • pl/sys-pre/ligand/dock/smina_91433-17-9-dock-PL.sdf:> <minimizedAffinity>
  • pl/sys-pre/ligand/dock/smina_91433-17-9-dock-PL.sdf--7.90676
绘制初始结构比较图,展示对接打分
D3DOCK的对接中,要仔细想想 二聚体/三聚体的氨基酸序号
先全部按照Tleap-MD的编号来,后面看情况把第二条链上的序号减去第一条链(经检查,是完全一致的)
PL_dimer 314*2 314ILE

Dynamic behavior of protein-ligand complex

rmsd,rmsf都在这儿展示
75.3 的数据基本都转移到75.2了
  • 75.2 /home/databank/zyzhou/project/hsbd/4bioactive-d3dock/analysis/rmsd/complex
  • # 之前的图片都是CAS号命名的,要转成TCMID
图已有,要看怎么描述
google scholar"rmsd source:Journal source:of source:the source:American source:Chemical source:Society"
看看高水平杂志怎么描述的
rmsf同理
mmgbsa,rmsd.rmsf已整理好
mmgbsa t检验的单双尾,type要在verify一下
还差残基能量分解和残基接触频率
  • 23 /home/databank_70t/zzy/project/hsbd/D3Targets-dock/workdir/5bioactive-dock-all-3cl_pl/analysis
就单纯的给出5个体系MD模拟较稳定的这个结论

Predicted highly frequent residues of protein-ligand binding

TCMID-27476 and TCMID-27300 MD中结合位点的证据:
ARG298 这个残基有点意思,可展开讨论
针对AT7519-Mpro 复合物做的MD模拟,可参考:
Mpro 可药位点的一些讨论:
PLpro也是类似的思路:
找该活性位点的关键氨基酸,是否在高接触频率的残基中,以此展开讨论。
  • 75.2 /home/databank/zyzhou/project/hsbd/4bioactive-d3dock/analysis/cluster/ligand_rmsd
  • # 做了根据配体RMSD的聚类