linux集群上的SCHRODINGER 分子对接工作流程

已于 2023-03-27 16:38:12 修改
阅读量4.8k 收藏 31
点赞数 13
分类专栏： 药学软件使用教程 文章标签： linux 药学 glide
于 2022-08-26 10:52:00 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/recher_he1107/article/details/126537558
版权

文章目录

前言

glide作为广为人知的分子对接软件，提供了非常方便的各类型对接工作，本篇文章旨在介绍如何方便的在无图形化界面的liinux集群中运行整个分子对接流程

一、总体工作流简介

在这里插入图片描述
总体流程大致可以输入，预处理和分子对接三个部分
• 输入为蛋白文件和小分子2D结构数据库。
• 输入文件经过ligprep和grid_gen_from_pv.py进行预处理， 其中ligprep是将小分子2D结构转3D。grid_gen_from_pv.py是从蛋白复合物中的已知配体生成蛋白结结合位点网格。
• 分子对接依赖于GSVrun程序进行，GSVrun包括多种常用的对接模式，以及对接过程中的可调参数。详见第4节。
• 输出文件包括每个对接阶段的小分子3D构象输出以及得分。

二、ligprep参数介绍

Bash
usage: ligprep [options] (-ismi|-icsv|-imae|-isd) infile (-osd|-omae) outfile

General:
  -h, -help             Print brief help message.
  -long_help            Print exhaustive help message.
  -inp <filename>       Read arguments from the supplied input file.
  -sif_docs             List supported simplified input file format (SIF)
                        keywords.

Ionization:
  -epik                 Use Epik for ionization and tautomerization
                        (Recommended, overrides -i).
  -emb, -epik_metal_binding
                        Run Epik with the metal_binding option so that states
                        appropriate for interactions with metal ions in
                        protein binding pockets are also generated.
  -i {0,1,2}            Ionization treatment: 0 - do not neutralize or ionize,
                        1 - neutralize only, 2 - neutralize and ionize
                        (Default: 1). Note that -epik option overrides -i and
                        always implies ionization and tautomerization.
  -ph <number>          Effective/target pH.
  -pht <number>         pH tolerance for generated structures.

Stereoisomers:
  -ac                   Do not respect existing chirality properties and do
                        not respect chiralities from the input geometry.
                        Generate stereoisomers for all chiral centers up to
                        the number permitted (specified using the -s option).
                        This is equivalent to "Generate all combinations" in
                        the Ligand Preparation user interface. Default
                        behavior is to respect only explicitly indicated
                        chiralities.
  -g                    Respect chiralities from input geometry when
                        generating stereoisomers.
  -s #                  Generate up to this many (#) stereoisomers per input
                        structure. (Default: 32).

Force-field based geometry optimization:
  -bff {14,16}          Force-field to be used for the final geometry
                        optimization. Default: 14 (OPLS_2005). For S-OPLS
                        specify 16.

Standard and job control options:
  -NJOBS NJOBS          Divide the overall job into NJOBS subjobs.
  -NSTRUCTS NSTRUCTS    Divide the overall job into subjobs with no more than
                        NSTRUCTS structures.
  -HOST <hostname>      Run job remotely on the indicated host entry.
  -WAIT                 Do not return a prompt until the job completes.
  -LOCAL                Do not use temporary directory for job files. Keep
                        files in the current directory.

Typical usage:
  Sample protonation states and stereoisomers, optimize geometry
  using OPLS force-field, and store results in Maestro format:

    ligprep <input> -omae outfile.maegz -epik

  where <input> is one of:
    -ismi infile.smi
    -icsv infile.csv
    -isd  infile.sd
    -imae infile.mae

一般情况下默认参数即可，在slurm集群上采用如下代码提交任务(ligprep的cpu利用率不高， 一般建议Njobs设置为cpu核数的2~3倍)：

srun -c 4 -J ${name}_ligprep /home/xiaxichen/hyd/Schrodinger-soft/ligprep -ismi ${name}.smi -omae ${name}.maegz -epik -NJOBS 12 -WAIT

三、grid_gen_from_pv.py 参数介绍

SCHRODINGER本身其实就提供了很多python脚本以实现各种功能
脚本路径：path/to/schrodinger/mmshare-v5.4/python/script
脚本执行方式 $SCHRODINGER/run XXX.py [Options]
脚本查询：https://www.schrodinger.com/scriptcenter

Bash
Options:
  -v, -version          Show the program version and exit.
  -h, -help             Show this help message and exit.
  -o OUT_BASE, -outbase OUT_BASE
                        Base name for gridgen input to be produced.
  -a, -allprint         Report all levels of status message.

可以看到参数非常少，连基本的网格大小都不可以设置。所以建议还是在本地创建好grid文件在传到集群上为宜，首先是所需时间其实也非常短，其次就是能够加深对蛋白结构的理解和更多的网格设置选项。

四、GVSrun参数介绍

Bash
Warning: The GVSrun need a $compound_library environmental variable, or you can use -d to specify a library.

Perform Virtual Screening Workflow.

Usage: GVSrun [OPTION] <parameter>

Input parameter:
  -i        Gird file input.
        Use a file name (Multiple files are wrapped in "", and split by ' ') or regular expression to represent your input Grid file, default is *.zip.
  -D    Which databases do you want to screen? The Database basic path is <>.
        These database are Supported:
            Custom_DB-ab-site-Fast  GVSrun  GVXrun_help.txt  plmd  README.md  Scripts  XDock  
  -d    Provide your own database path, the compounds files are recommended as maegz or SDF file format.
  -R    Optional, reference Ligands correlated with the grid or use for RMSD calculation. 
            As *.mae, *.maegz (for query structure) or *.phypo (for pharmacophore hypothesis).
  -m    Running Mode: a serial combination of different computing tasks.
            Show all mode and task in this option using "-O".
    @ Available Running Mode: <Fast>
        Fast: Fast Virtual Screening.
        Normal: Filter The Drug-like compounds and dock screening.
        Reference: Virtual Screening with reference ligand restrain.
        Prep_Normal: Virtual Screening for un-prepared compounds database.
        *Normal_MMGBSA: Virtual Screening and MMGBSA re-scoring.
        *Cov_Screening: Virtual Screening to discover covalent durg.
        *Induce_Fit_Screening: Induce fit screening.
        *QM_Screening: Virtual Screening and QMMM re-docking.
        Local: Local docking and screening were carried out using the input ligand structure.
        Shape_Screening: perform screening based on reference ligand shape.
        Advance: dock screening and clustering.
        Advance_Shape_Screening: perform screening based on reference ligand shape.
    @ or Custom Running Mode, e.g. -m "EDL+R+HTVS+CD"

Control parameters:
  -F    Force Fields for docking stage (SP/XP/MMGBSA...), OPLS_2005, OPLS3e or OPLS4. <OPLS4>
  -f    Force Fields for other stage (HTVS/LigPrep...), OPLS_2005, OPLS3e or OPLS4. <OPLS_2005>
  -W    Define a Max/Min Molecular weight for MW module, such as "100:400" <100:400>
  -T    Set a Job Name. Default is "Grid_name-Database_name-Run_Mode".
  -C    Aattach residue number on receptor, required in Covalent Docking.
            e.g. "cys:A:1425", the A is chain name and 145 is the atom number(Heavy atom). 
            The cys is residue name for A:1425, Supported residues: cys, ser, lys.
  -q    Define a "DFT:Basis_Set" to QM/QMMM, default is "B3LYP-D3(BJ):6-311G**".
            Other QM setting: "B3LYP-D3M(BJ):6-311G+**","M06‑2X:def2-tzvpp(-g)","wB97M‑V:cc-pVTZ-pp"
  -p    set a PH:PHT for ligPrep, default is 7.0:2.0, means 7.0±2.0. <7.0:2.0>
  -s    Set a SMARTs Expression for compounds filter at first step. such as [B]([O])[O].
  -E    Shape Screening Options. <10000:mmod:rapid:1200>
            The paramter means: "keep_num:atomtypes:shape_sample_method:max_confs";
            NOTE: element and qsar are alternative atomtypes; thorough are alternative shape_sample_method.

OUTPUT parameters:
  -a    The number of Output compounds per Screening Task. <5%>
            e.g. 10% means reatin top 10% compounds.
                10000 means reatin top 10000 compounds.
  -b    The number of Output compounds after Standard docking Task. <4000>
  -c    The number of conformations generated by each ligand in the docking task. <1>
  -e    The number of candidates to IFT/CD/MMGBSA/QMMM. <500>

Job control:
  -H        Host Name of your Queue, defult is CPU.
  -N    The max number of subjobs. <100>
  -G    The number of Glide subjobs. <90>
  -P    The number of Prime subjobs. <10>
  -L    The number of LigPrep subjobs. <30>
  -A    The number of phase subjobs. <10>
  -Q    The number of Qsite subjobs. <10>
  -K    The number of QIKPROP subjobs. <10>
  -M    The number of MACROMODEL subjobs. <10>
  -S        Your Schrodinger path. </home/xiaxichen/hyd/Schrodinger-soft>

GVSrun通过环境变量compound_library来获取系统中化合物数据库的位置；如果不存在该变量需要利用-d参数给定化合物数据库路径（纯SDF或MAEGZ格式的多个文件，不可套叠文件夹）除此以外，脚本需要知道安装好的薛定谔路径：

  • ~/.bashrc中指定好**$SCHRODINGER**环境变量，这也是正常安装薛定谔已经完成的操作；
  • 也可以修改脚本，直接更改默认的SCHRODINGER路径（脚本第22行修改);
  • 可以使用各种参数进行组合进行所需分子对接的定制，在这里重点介绍Running mode参数（其余参数如上所示就不加赘述了），GVSrun一共拥有12种Running mode分别对应这常见的12种分子对接方案，他们分别是：

在这里插入图片描述

  • 以及如何制定新的Running_model（以Advance为例），每个Running_model其实就是个各项任务的组合，对于新的Running_model制定所需要添加的内容有三个部分：方法注释，运行逻辑， 任务运行代码：
#===========================
#添加方法注释
#===========================
@Available Running Mode:
     Advance: dock screening.
    "HTVS_Normal+SP_ExtensionA+SP_Enhanced" #任务组合，代码93行

@Screening: Reduce the number of compounds rough docking. Related to -a option.
    HTVS_Normal: The standard HTVS screening.

@Standard_Docking: Standard docking was performed to rank the compounds. Related to -b and -c option.
    SP_ExtensionA: Reward intramolecular ligand hydrogen bonds, and accept halogens as H-bond acceptors.
    SP_Enhanced: Increase the depth of sampling for larger grid or More accurate poses. 

#===========================
#添加运行逻辑
#===========================
Parse_Running_Mode(){
    # Parse_Running_Mode Running_Mode
    Running_Mode=$1
    if [ $Running_Mode == "Normal" ]; then
        Pipeline="RDL+HTVS_Normal+QIKPROP+R5R+SP_ExtensionA"
    elif [ $Running_Mode == "Prep_Normal" ]; then
        Pipeline="No_Dup+RDL+IONIZE+HTVS_Normal+QIKPROP+R5R+SP_ExtensionA"
    elif [ $Running_Mode == "Normal_MMGBSA" ]; then
        Pipeline="RDL+HTVS_Normal+QIKPROP+R5R+SP_ExtensionA+MMGBSA_EN"
    elif [ $Running_Mode == "Reference" ]; then
        Pipeline="HTVS_REF+SP_REF+QIKPROP+R5R"
    elif [ $Running_Mode == "Induce_Fit_Screening" ]; then
        Pipeline="IFT_pre+IFT"
    elif [ $Running_Mode == "Cov_Screening" ]; then
        Pipeline="R+RDL+HTVS_Normal+SP_ExtensionA+CD"
    elif [ $Running_Mode == "Fast" ]; then
        Pipeline="HTVS_Normal+SP_Normal"
    elif [ $Running_Mode == "QM_Screening" ]; then
        Pipeline="RDL+HTVS_Normal+QIKPROP+R5R+SP_ExtensionA+QM_redock+RMSD"
    elif [ $Running_Mode == "Shape_Screening" ];then
        Pipeline="PhaseShape+HTVS_local+SP_local"
    elif [ $Running_Mode == "Local" ];then
        Pipeline="HTVS_local+SP_local+MMGBSA_OPT"
    elif [ $Running_Mode == "Advance" ];then
        Pipeline="HTVS_Normal+SP_ExtensionA+SP_Enhanced"
    elif [ $Running_Mode == "Advance_Shape_Screening" ];then
        Pipeline="HTVS_Shape+SP_Shape"
    else 
        Pipeline=$Running_Mode
        Running_Mode="User_Defined"
    fi
}

#===========================
#添加任务运行代码
#===========================
HTVS_Normal(){
    KEEP_NUM=`parse_num_or_percentage ${HTVS_out_num}`
cat<<HTVS >> ${Inp_Name}
[STAGE:PRE_DOCK_HTVS]
    STAGECLASS   gencodes.RecombineStage
    INPUTS   ${Ligand_name},
    OUTPUTS   HTVS_RECOMBINE_OUT,
    NUMOUT   njobs
    OUTFORMAT   maegz
    MIN_SUBJOB_STS   4000
    MAX_SUBJOB_STS   40000
    GENCODES   YES
    OUTCOMPOUNDFIELD   s_vsw_compound_code
    OUTVARIANTFIELD   s_vsw_variant
    UNIQUEFIELD   NONE
[STAGE:DOCK_HTVS]
    STAGECLASS   glide.DockingStage
    INPUTS   HTVS_RECOMBINE_OUT, GRID
    OUTPUTS   HTVS_OUT,
    RECOMBINE   NO
    PRECISION   HTVS
    UNIQUEFIELD   s_vsw_compound_code
    ${KEEP_NUM}

五、相关文件

grid_gen_from_pv：
链接：https://pan.baidu.com/s/1ymh7e2eesWmU5yMFgxK6mA?pwd=6666
GVSrun：
https://pan.baidu.com/s/1wTLhDBruY5F0iA-ndx944A?pwd=6666

六、时间预估

流程名速度1万个分子所需时长
ligprep1.5 秒/分子4.1 小时/1万分子
glide_HTVS0.5 秒/分子1.5 小时/1万分子
glide_SP9 秒/分子25 小时/1万分子
glide_XP100 秒/分子280 小时/1万分子

六、 特别感谢

感谢上海科技大学的王林写的GVSrun脚本并将其开源，用起来非常方便。附上他的github主页：
https://github.com/Wang-Lin-boop

贺俊宏
关注
专栏目录
Schrodinger软件介绍（以及Glide的工作流）
01-24
1、计算机辅助药物设计薛定谔软件； 2、虚拟筛选之分子对接； 3、Glide。
Linux的常见问题解决
SchrodingerY-博客
02-09 422
这些都是实际使用过程中遇到的问题，记录一下，同时提供给大家参考 删除 /mnt/sdb/backups文件夹下 15天前的 文件名包含 “gz” 的文件： find /mnt/sdb/backups -mtime +15 -name "*gz*" -exec rm -rf {} \; scp文件传输命令: scp local_file remote_username@remote_ip...
分子对接软件（薛定谔Schrodinger) 使用
pang0416的博客
03-16 7213
蛋白-配体相互作用、柔性对接及基于结构的虚拟筛选、基于药效团虚拟筛选、构效关系模型预测分子活性
schrodinger 薛定谔安装与卸载
发呆的比目鱼的博客
04-04 7445
schrodinger 薛定谔安装
schrodinger 薛定谔
周迪新的博客
09-24 7239
链接: https://pan.baidu.com/s/1Duh4GsoytIklvZp7TdHpWg 密码: wmr1 Linux系统版 进入 Schrodinger_Suites_2021-2_Linux-x86_64 sudo ./INSTALL 中间按默认的来即可，也可以自己调整安装目录。安装目录一定要记住。 安装完，把patcher 拷贝到安装目录 sudo ./patcher 点击Patch即可。 启动程序，进入安装目录 ./maestro ...
使用Schrödinger Python API系列教程 -- 介绍 (一)
发呆的比目鱼的博客
04-12 1703
使用Schrödinger Python API系列教程(一)
Linux系统安装AutoDockTools、AutoGrid和AutoDock并实现分子对接（详细讲解）
热门推荐
carlos_zhao的博客
02-09 1万+
linux系统为centos7.0 按照下面操作安装AutoDockTools（MGLTools） 1. 下载安装包，链接（https://ccsb.scripps.edu/mgltools/downloads/），我的linux系统是64位，下载红色标记源代码，拷贝入linux目标文件夹下。 2. shell终端进行安装： 1. 终端进入安装包所在文件夹，并获取root权限；赋予代码安装权限； cd /home/... #目标文件夹，根据目标程序文件位置进行更改 su #获
linux医学软件,Schrodinger Suites 2015-2 Linux版简单安装经验分享
weixin_29234239的博客
05-03 2182
首先感谢网友蔡剑雄药物制剂提供的安装文件。我是在ubuntu16.04桌面系统里安装的。1. 挂载schrodinger2015-2.iso2. 打开终端，进入刚挂载的安装光盘所在目录(cd /nmt/schrodinger2015-2 [这是我的挂载位置])3. 按照安装文件说明，先检查系统要求(./platform -s)，检查机器配置(./platform -l)。如果没有问题就进行下一步...
2023年linux镜像下载地址，包括linux官网、国内镜像下载，超详细教学，小白也能学会
04-03
Linux操作系统是世界上最流行的开源操作系统之一，它被广泛用于服务器、超级计算机、嵌入式设备以及个人电脑等场景。本文将详细介绍如何在2023年下载Linux镜像文件，包括从Linux官网以及国内的一些镜像站点获取镜像...
SSFT.rar_Schrodinger matlab_ssft_薛定谔方程
07-15
在“SSFT.m”文件中，我们可以期待看到如何利用MATLAB的内置函数或者自定义算法来实现时间或空间域上的数值解法。 短时傅里叶变换是STFT的核心，它通过将信号分解成一系列重叠的小窗口，并对每个窗口进行傅里叶变换...
SchrodingerWellPython:二维任意时间形状的Schrodinger方程求解器
03-25
“Schrodinger井Python”是一个强大的工具，能够处理任意形状的二维势阱，对于研究具有复杂几何结构的量子系统非常有用。此程序基于Python编程语言，利用其丰富的科学计算库，如NumPy、SciPy等，进行数值求解。...
Split Step Fourier Method.rar_method_schrodinger_split_split ste
07-14
This M-File solves Nonlinear schrodinger equation and display the results in 3d graphics along with it the Pulse broadening ratio and phase shift has also been calculated and displayed
薛定谔maestro 安装
weixin_63016274的博客
06-13 2570
maestro 安装
Maestro 薛定谔软件简单分子对接案例
weixin_42357472的博客
01-15 6985
#参考：Maestro 薛定谔软件使用：Maestro 薛定谔对接：sid=287811。
Linux系统安装AutoDockTools、AutoGrid和AutoDock并实现分子对接 分子对接软件 CentOs Ubuntu系统适用（详细讲解）
GHZ2443063059的博客
09-10 2269
linux系统Centos7 | Ubuntu系统安装分子对接软件AutoDockTools并实现分子对接
redis集群搭建
003号猫
08-04 122
1、安装redis（每台实例都是相同操作） redis安装在Linux上，redis是c语言开发的，需要依赖以下环境。 yum -y install cpp binutils glibc glibc-kernheaders glibc-common glibc-devel gcc make gcc-c++ libstdc++-devel tcl #如果安装过则不会执行（没有规则可以创建目标“gcc-c++”。 停止） make gcc-c++ libstdc++-devel tcl #进入root目录，
Linux连接脚本详解
Olivertien的专栏
10-15 663
转载自：http://blog.csdn.net/jiangyuping_fyl/article/details/7293420 一、 概论 每一个链接过程都由链接脚本(linker script, 一般以lds作为文件的后缀名)控制. 链接脚本主要用于规定如何把输入文件内的section放入输出文件内, 并控制输出文件内各部分在程序地址空间内的布局. 但你也可以用连接命令做一
AutoDock Vina批量分子对接|Linux系统+Windows系统|R语言+Bash+Python|22.6.17更新
Squirrelity的博客
06-17 1861
本文是关于一个受体与多个配体的分子对接批处理。 由于作者本人也是赶鸭子上架，写出来的脚本比较简单，请见谅。 此外，由于对某个语言某个包的喜爱与使用，作者用到的编程工具比较杂，但是可以作为思路参考。当然，如果您具备条件，也可以选择复刻，本文章的代码都被跑过，请根据自己需要修改参数（内含中文不能直接跑）与编程。.........
【Linux】：信号与信号产生
最新发布
Yikefore的博客
09-12 934
信号的基本概念、对信号的理解、信号的多种产生方式以及核心转储的详细介绍！
薛定谔软件如何进行分子对接
06-10
薛定谔软件(Schrodinger)提供了多个分子对接工具，其中较为常用的是Glide和Docking等。以下是Glide和Docking进行分子对接的基本步骤： 1. Glide分子对接： - 打开Schrodinger软件并加载配体和受体分子结构。 - 在主菜单栏中选择“Glide”>“Ligand Docking”。 - 在Glide对接向导中，设置对接相关参数，例如搜索方式、打分函数等。 - 点击“Run”按钮开始对接计算。 - 对接结果可以在Glide工作区中查看和分析。 2. Docking分子对接： - 打开Schrodinger软件并加载配体和受体分子结构。 - 在主菜单栏中选择“Docking”>“Ligand Docking”。 - 在Docking对接向导中，设置对接相关参数，例如搜索方式、打分函数等。 - 点击“Run”按钮开始对接计算。 - 对接结果可以在Docking工作区中查看和分析。 以上仅是分子对接的基本步骤，具体操作细节和参数设置需要根据具体情况进行调整。另外，分子对接是一个复杂的计算过程，需要充分理解原理和使用方法，以获得准确的结果。
评论 12
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值