23a0d3bfa8246b8f93fc515667024a1270f61a00
embuild.md
... | ... | @@ -1,75 +0,0 @@ |
1 | -#!/usr/bin/env bash |
|
2 | - |
|
3 | -## SBGrid EMBuild example - EMD-0346 |
|
4 | -## http://huanglab.phys.hust.edu.cn/EMBuild/ |
|
5 | -## Args: none |
|
6 | -## |
|
7 | -## James Vincent - biogrids.org |
|
8 | -## vincent@hkl.hms.harvard.edu |
|
9 | -## Mar 08, 2023 |
|
10 | - |
|
11 | -:<<ACKNOWLEDGE |
|
12 | - |
|
13 | -If your use of SBGrid compiled software was an important element in |
|
14 | -your publication, please include the following statement in your work: |
|
15 | -"Software used in the project was installed and configured by SBGrid |
|
16 | -(cite: eLife 2013;2:e01456, Collaboration gets the most out of software.) |
|
17 | - |
|
18 | -ACKNOWLEDGE |
|
19 | - |
|
20 | -# BCH E2 slurm |
|
21 | -#SBATCH --partition=bch-gpu |
|
22 | -#SBATCH --gres=gpu:Tesla_K:1 |
|
23 | -#SBATCH --mem-per-cpu=24G |
|
24 | -#SBATCH --time=01:00:00 |
|
25 | -#SBATCH --job-name=JV-embuild |
|
26 | - |
|
27 | - |
|
28 | -# STRIDE dependency must be downloaded/installed by user |
|
29 | -curl -kLO http://webclu.bio.wzw.tum.de/stride/stride.tar.gz |
|
30 | -# Untars in current directory - make a subdir first |
|
31 | -mkdir stride |
|
32 | -mv stride.tar.gz stride |
|
33 | -cd stride |
|
34 | -tar zxf stride.tar.gz |
|
35 | -make |
|
36 | -# Absolute stride path required in call to EMBuild |
|
37 | -stride_path=`pwd`/stride |
|
38 | -echo "STRIDE path: ${stride_path} " |
|
39 | -cd .. |
|
40 | - |
|
41 | - |
|
42 | -# Load cuda |
|
43 | -module load cuda/11.2 |
|
44 | - |
|
45 | -# Start SBGrid environment |
|
46 | -source /programs/sbgrid.shrc |
|
47 | - |
|
48 | -# MCP dir needed in calls |
|
49 | -mcp_dir=/programs/x86_64-linux/embuild/1.0/bin/mcp |
|
50 | - |
|
51 | -# embuild requires large ulimit - segfault otherwise |
|
52 | -ulimit -s 1048576 |
|
53 | - |
|
54 | - |
|
55 | -# Download example data |
|
56 | -curl -kLO http://huanglab.phys.hust.edu.cn/EMBuild/6N52.tgz |
|
57 | -tar zxf 6N52.tgz |
|
58 | -cd ./6N52/ |
|
59 | - |
|
60 | -python.embuild ${mcp_dir}/mcp-predict.py -i emd_0346.map -o 0346_MC.mrc -m ${mcp_dir}/model_state_dict |
|
61 | - |
|
62 | -SWORD -i 6N52_A_AF2.pdb -m 15 -v > 6N52_A_AF2_SWORD.out |
|
63 | -SWORD -i 6N52_B_AF2.pdb -m 15 -v > 6N52_B_AF2_SWORD.out |
|
64 | - |
|
65 | -asgdom 6N52_A_AF2.pdb 6N52_A_AF2_SWORD.out 6N52_A_AF2_doms.pdb |
|
66 | -asgdom 6N52_B_AF2.pdb 6N52_B_AF2_SWORD.out 6N52_B_AF2_doms.pdb |
|
67 | - |
|
68 | -cat *_doms.pdb > init_chains.pdb |
|
69 | - |
|
70 | -EMBuild 0346_MC.mrc init_chains.pdb 4.0 6N52_EMBuild.pdb |
|
71 | - |
|
72 | -rearrangepdb 6N52_EMBuild.pdb 6N52_EMBuild_rearranged.pdb |
|
73 | - |
|
74 | -EMBuild 0346_MC.mrc 6N52_EMBuild_rearranged.pdb 4.0 6N52_EMBuild_scored.pdb \ |
|
75 | - --score -stride ${stride_path} |
|
... | ... | \ No newline at end of file |