解决问题
蛋白模拟
小小编程
机器学习
要看文献
成长思考
自娱自乐
vina批量对接
type
status
date
slug
summary
tags
category
icon
password
vina --help
prepare_receptor
成功的workflow
直接用rdkit将smiles转换为pdbqt
或者用ADFR将mol2转换为pdbqt。(但是用rdkit将smiles转为mol2的格式有问题的。可以用ds导出mol2.)下载ADFR,设置环境变量。
prepare_receptor -r dpp4.pdb
prepare_ligand -l CDSS.mol2
直接用autodock打开受体文件找到对接盒子大小Grid-Grid box,记得长宽高乘那个放大比例。
写好config文件
vina --config config.txt
config文件
遇到的问题
用openbabel将fasta直接转为mol2,
obabel -ifasta *.fasta -omol2 -m -h -gen3d
能转换,但是全都重复生成两个。
先转换为sdf
obabel -ifasta *.fasta -osdf -O *.sdf --gen3d
再转为mol2
obabel -isdf *.sdf -omol2 -m -h -gen3d
报错,变成了一个原点,所有的原子坐标在一个位置。
dock.bat文件
验证了一个问题。
说明写的bat脚本有问题
参考资料
反思:不要解决不存在的问题。直接用sequence转换出问题就改为smiles,不要折磨自己,能解决就行,管他呢。需要一定的思考钻研,但不要把自己卡进去。
Loading...
Last update: 2025-03-28
🎉终于搭建完成🎉
-- 感谢您的支持 ---
👏欢迎阅读👏