examples/embuild
... ...
@@ -0,0 +1,75 @@
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