- Define merging opportunities by creating tags of LHS hits in Fragalysis
- Download target from Fragalysis and place the .zip archive in the repo
- Setup target in BulkDock
cp -v ZIKA_NS3_helicase.zip $BULK/TARGETS
cd $BULK
python -m bulkdock extract ZIKA_NS3_helicase
python -m bulkdock setup ZIKA_NS3_helicase
- Copy the
aligned_files
directory from$BULK/TARGETS/ZIKA_NS3_helicase/aligned_files
into this repository
cd -
cp -rv $BULK/TARGETS/ZIKA_NS3_helicase/aligned_files .
- run the notebook
hippo/1_merge_prep.ipynb
For each merging hypothesis (i.e. RNA-Cleft)
- go to the fragmenstein subdirectory
cd fragmenstein
- queue fragmenstein job
sb.sh --job-name "ZIKA_NS3_helicase_RNA-Cleft_fragmenstein" --mem 16000 $HOME2/slurm/run_bash_with_conda.sh run_fragmenstein.sh RNA-Cleft
This will create outputs in the chosen RNA-Cleft subdirectory:
-
RNA-Cleft_fstein_bulkdock_input.csv
: use this for BulkDock placement -
output
: fragmenstein merge poses in subdirectories -
output.sdf
: fragmenstein merge ligand conformers -
output.csv
: fragmenstein merge metadata -
placement with bulkdock
cp -v RNA-Cleft_fstein_bulkdock_input.csv $BULK/INPUTS/TARGET_RNA-Cleft_fstein.csv
cd $BULK
python -m bulkdock place ZIKA_NS3_helicase INPUTS/TARGET_RNA-Cleft_fstein.csv
- monitor placements (until complete)
python -m bulkdock status
- export Fragalysis SDF
sb.sh --job-name "bulkdock_out" $HOME2/slurm/run_python.sh -m bulkdock to-fragalysis TARGET OUTPUTS/SDF_FILE METHOD_NAME
Running Fragment Knitting currently requires access to a specific VM known as graph-sw-2
. If you don't have access, skip this section
-
git add
,commit
andpush
the contents ofaligned_files
andknitwork
to the repository -
git clone
the repository ongraph-sw-2
- navigate to the
knitwork
subdirectory
Then, for each merging hypothesis:
- Run the "fragment" step of FragmentKnitwork:
./run_fragment.sh RNA-Cleft
- Run the pure "knitting" step of FragmentKnitwork:
./run_knitwork_pure.sh RNA-Cleft
- Run the impure "knitting" step of FragmentKnitwork:
./run_knitwork_impure.sh RNA-Cleft
- Create the BulkDock inputs:
python to_bulkdock.py RNA-Cleft
-
git add
,commit
andpush
the CSVs created by the previous step - back on
cepheus-slurm
pull the latest changes - Run BulkDock placement as for Fragmenstein above
cp -v RNA-Cleft_knitwork_pure.csv $BULK/INPUTS/TARGET_RNA-Cleft_knitwork_pure.csv
cp -v RNA-Cleft_knitwork_impure.csv $BULK/INPUTS/TARGET_RNA-Cleft_knitwork_impure.csv
cd $BULK
python -m bulkdock place ZIKA_NS3_helicase INPUTS/TARGET_RNA-Cleft_knitwork_pure.csv
python -m bulkdock place ZIKA_NS3_helicase INPUTS/TARGET_RNA-Cleft_knitwork_impure.csv
- Export Fragalysis SDF as for Fragmenstein