Skip to content

Commit fe0fd33

Browse files
committed
mods to screen nextflow workflows
1 parent d53c49e commit fe0fd33

File tree

4 files changed

+40
-52
lines changed

4 files changed

+40
-52
lines changed

src/nextflow/docking/rdock.nf

Lines changed: 0 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -65,8 +65,6 @@ process rdock {
6565
*/
6666
process results {
6767

68-
publishDir './', mode: 'copy'
69-
7068
input:
7169
file part from docked_parts.collect()
7270

src/nextflow/rdkit/screen-dataset.nsd.nf

Lines changed: 8 additions & 24 deletions
Original file line numberDiff line numberDiff line change
@@ -14,7 +14,7 @@ target = file(params.input)
1414

1515
process splitter {
1616

17-
container 'informaticsmatters/rdkit_pipelines'
17+
container 'informaticsmatters/rdkit_pipelines:latest'
1818

1919
input:
2020
file target
@@ -45,40 +45,24 @@ process rdkitScreen {
4545

4646
process joiner {
4747

48-
container 'informaticsmatters/rdkit_pipelines'
48+
container 'informaticsmatters/rdkit_pipelines:latest'
4949

50-
publishDir baseDir, pattern: "{output.data.gz,output.metadata}"
50+
publishDir baseDir, mode: 'link'
5151

5252
input:
53+
file 'splitter_metrics.txt' from splitter_metrics
5354
file parts from screened_parts.collect()
5455

5556
output:
56-
file 'output_metrics.txt' into joiner_metrics
57+
file 'output_metrics.txt'
5758
file 'output.data.gz'
5859
file 'output.metadata'
5960

6061
"""
6162
zcat $parts | python -m pipelines_utils_rdkit.filter -if sdf -of json -o output --meta
62-
"""
63-
}
64-
65-
process metrics {
66-
67-
container 'informaticsmatters/rdkit_pipelines'
68-
69-
publishDir baseDir
70-
71-
input:
72-
file 'splitter_metrics.txt' from splitter_metrics
73-
file 'joiner_metrics.txt' from joiner_metrics
74-
75-
output:
76-
file 'output_metrics.txt'
77-
78-
"""
79-
grep '__InputCount__' splitter_metrics.txt | sed s/__InputCount__/RDKitScreen/ > output_metrics.txt
63+
mv output_metrics.txt joiner_metrics.txt
64+
grep '__InputCount__' splitter_metrics.txt | sed s/__InputCount__/RDKitScreen/ > output_metrics.txt
8065
grep '__InputCount__' splitter_metrics.txt >> output_metrics.txt
8166
grep '__OutputCount__' joiner_metrics.txt >> output_metrics.txt
82-
"""
67+
"""
8368
}
84-

src/nextflow/rdkit/screen-multi-dataset.nsd.nf

Lines changed: 9 additions & 25 deletions
Original file line numberDiff line numberDiff line change
@@ -15,7 +15,7 @@ query = file(params.query)
1515

1616
process splitter {
1717

18-
container 'informaticsmatters/rdkit_pipelines'
18+
container 'informaticsmatters/rdkit_pipelines:latest'
1919

2020
input:
2121
file target
@@ -31,7 +31,7 @@ process splitter {
3131

3232
process rdkitScreen {
3333

34-
container 'informaticsmatters/rdkit_pipelines'
34+
container 'informaticsmatters/rdkit_pipelines:latest'
3535

3636
input:
3737
file part from target_parts
@@ -47,40 +47,24 @@ process rdkitScreen {
4747

4848
process joiner {
4949

50-
container 'informaticsmatters/rdkit_pipelines'
50+
container 'informaticsmatters/rdkit_pipelines:latest'
5151

52-
publishDir baseDir, pattern: "{output.data.gz,output.metadata}"
52+
publishDir baseDir, mode: 'link'
5353

5454
input:
55+
file 'splitter_metrics.txt' from splitter_metrics
5556
file parts from screened_parts.collect()
5657

5758
output:
58-
file 'output_metrics.txt' into joiner_metrics
59+
file 'output_metrics.txt'
5960
file 'output.data.gz'
6061
file 'output.metadata'
6162

6263
"""
6364
zcat $parts | python -m pipelines_utils_rdkit.filter -if sdf -of json -o output --meta
64-
"""
65-
}
66-
67-
process metrics {
68-
69-
container 'informaticsmatters/rdkit_pipelines'
70-
71-
publishDir baseDir
72-
73-
input:
74-
file 'splitter_metrics.txt' from splitter_metrics
75-
file 'joiner_metrics.txt' from joiner_metrics
76-
77-
output:
78-
file 'output_metrics.txt'
79-
80-
"""
81-
grep '__InputCount__' splitter_metrics.txt | sed s/__InputCount__/RDKitScreen/ > output_metrics.txt
65+
mv output_metrics.txt joiner_metrics.txt
66+
grep '__InputCount__' splitter_metrics.txt | sed s/__InputCount__/RDKitScreen/ > output_metrics.txt
8267
grep '__InputCount__' splitter_metrics.txt >> output_metrics.txt
8368
grep '__OutputCount__' joiner_metrics.txt >> output_metrics.txt
84-
"""
69+
"""
8570
}
86-

test-nextflow.sh

Lines changed: 23 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -4,7 +4,7 @@
44

55
set -e
66

7-
echo 'Running screen+config'
7+
echo 'Running screen+conformers in basic mode'
88
nextflow run src/nextflow/rdkit/screen+conformers.nf -c src/nextflow/rdkit/screen.config -with-docker
99

1010

@@ -55,4 +55,26 @@ ln ../data/dhfr_3d.metadata ligands.metadata
5555
docker run -it --rm -v $PWD:$PWD:z -w $PWD -v /var/run/docker.sock:/var/run/docker.sock informaticsmatters/nextflow-docker:0.30.2 sh -c 'nextflow run nextflow.nf -c nextflow.config --num_dockings 1 --limit 40 --chunk 5 -with-docker'
5656
cd ..
5757

58+
echo 'Running screen in squonk mode'
59+
sudo rm -rf tmp/*
60+
cd tmp
61+
ln ../src/nextflow/rdkit/screen-dataset.nsd.nf nextflow.nf
62+
ln ../src/nextflow/rdkit/screen-dataset.nsd.config nextflow.config
63+
ln ../data/dhfr_3d.data.gz input.data.gz
64+
ln ../data/dhfr_3d.metadata ligands.metadata
65+
docker run -it --rm -v $PWD:$PWD:z -w $PWD -v /var/run/docker.sock:/var/run/docker.sock informaticsmatters/nextflow-docker:0.30.2\
66+
sh -c 'nextflow run nextflow.nf -c nextflow.config -with-docker --chunk 100 --simmin 0.5 --qsmiles "OC(=O)C1=CC=C(NC2=NC3=C(CN=C(C4=CC(Cl)=CC=C34)C3=C(F)C=CC=C3F)C=N2)C=C1"'
67+
cd ..
68+
69+
echo 'Running screen-multi in squonk mode'
70+
sudo rm -rf tmp/*
71+
cd tmp
72+
ln ../src/nextflow/rdkit/screen-multi-dataset.nsd.nf nextflow.nf
73+
ln ../src/nextflow/rdkit/screen-multi-dataset.nsd.config nextflow.config
74+
ln ../data/dhfr_3d.data.gz target.data.gz
75+
ln ../data/nci100.data.gz query.data.gz
76+
docker run -it --rm -v $PWD:$PWD:z -w $PWD -v /var/run/docker.sock:/var/run/docker.sock informaticsmatters/nextflow-docker:0.30.2\
77+
sh -c 'nextflow run nextflow.nf -c nextflow.config -with-docker --chunk 100 --simmin 0.55'
78+
cd ..
79+
5880
sudo rm -rf tmp/*

0 commit comments

Comments
 (0)