Jerkwin分享 http://blog.sciencenet.cn/u/Jerkwin

博文

吴伟:auto-martini的安装和使用简介

已有 351 次阅读 2019-8-6 07:07 |系统分类:科研笔记

  • 2019-08-05 17:50:22

在文献Tristan Bereau, Kurt Kremer; Automated Parametrization of the Coarse-Grained Martini Force Field for Small Organic Molecules; J. Chem. Theory Comput. 11(6):2783-2791, 2015; 10.1021/acs.jctc.5b00056中提出了自动将小分子划分为Martini珠子的方法, 作者还提供了一个python脚本. 这里简单介绍下它的使用方法.

安装

  1. 首先从github下载auto-martini程序包放到合适目录

  2. 下载Windows版的anaconda2. 注意, 这里要选择python2.7版本的, 因为auto-martini是用python2.7写的, 不兼容python3之后的语法

  3. 安装完anaconda之后, 以管理员身份运行anaconda, 进入auto-martini包所在位置

  4. 安装auto-martini所依赖的四个程序包

pip install numpypip install beautifulsouppip install requestsconda create -c rdkit -n my-rdkit-env rdkit

安装完之后, 就可以使用了.

使用方法

命令

python auto-martini [-h] (--sdf SDF | --smi SMI) --mol MOLNAME [--xyz XYZ] [--gro GRO] [--verbose] [--fpred]

说明

  • --sdf--smi: 输入文件, 指定其中一个就可以. 可以使用openbabel将pdb或其他格式的文件转化为sdf或者smi文件.

  • --mol: 必须选项, 输出文件中残基的名称

  • --xyz--gro: 可选的输出文件

  • --verbose--fpred: 无法找到符合的参数时, 使用按原子或按片段判别珠子的方法, 准确度较差

示例

使用Gaussview画出想要的小分子

可以直接保存为sdf文件, 但好像Gaussview生成的sdf文件在auto-martini中使用会出错(未详细考察), 所以我们保存为pdb文件, 再用openbabel转换为sdf文件.

将该文件放到auto-martini目录下. 执行

python auto_martini --sdf TEG.sdf --mol TEG --gro TEG_CG.gro

会得到如下结果

同时生成gro文件

注意事项

  • 如果提示smi2alogps这一步出错, 那你得先检查一下 http://vcclab.org/web/alogps 网站是否能够打开, 因为auto-martini会将你的分子上传到alogps网站上, 生成所需的logP等信息

◆本文地址: https://jerkwin.github.io/2019/08/05/吴伟-auto-martini的安装和使用简介/, 转载请注明◆
◆评论问题: https://jerkwin.herokuapp.com/category/3/博客, 欢迎留言◆



http://blog.sciencenet.cn/blog-548663-1192582.html

上一篇:gmx_mmpbsa使用说明
下一篇:物美价廉:GROMACS 2018在GPU节点上的使用

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备14006957 )

GMT+8, 2019-9-15 16:12

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部