Gaucher Disease: Task 09 - Lab Journal: run minimise.sh
run_minimise.sh
<source lang='bash'>
- !/bin/sh
- first run repairPDB to remove hydrogens, waters and hetero atoms:
- /opt/SS12-Practical/scripts/repairPDB <input.pdb> -noh -jprot > <output.pdb>
- WT
cd ~/master_practical/Assignment9_Structure-based_mutation_analysis/ /opt/SS12-Practical/scripts/repairPDB scwrl/2V3E_B_WT_backbone_IN.pdb -noh -jprot > minimise/input_noh_nohoh/2V3E_B_WT.pdb
- SCWRL
for mutation in "S38R" "N102S" "G202E" "N370S" "L444P"; do /opt/SS12-Practical/scripts/repairPDB scwrl/2V3E_B_${mutation}_model_OUT.pdb -noh -jprot > minimise/input_noh_nohoh/scwrl/2V3E_B_${mutation}_scwrl.pdb; done
- foldX
for i in 1 2 3 4 5; do /opt/SS12-Practical/scripts/repairPDB foldX/output_files/2V3E_B_${i}.pdb -noh -jprot > minimise/input_noh_nohoh/foldx/2V3E_B_${mutation}_foldx.pdb; done
- then minimise:
- /opt/SS12-Practical/minimise/minimise <input.pdb> <output.pdb> > <log file>
- WT
cd ~/master_practical/Assignment9_Structure-based_mutation_analysis/minimise/output/ minimise="/opt/SS12-Practical/minimise/minimise"
base="2V3E_B_WT" $minimise ../input_noh_nohoh/${base}.pdb wt_minimised/${base}_minimised_round1.pdb > wt_minimised/${base}_minimised_round1_log.txt
- 4 more rounds:
for n in 2 3 4 5; do let p=$n-1; $minimise wt_minimised/${base}_minimised_round${p}.pdb wt_minimised/${base}_minimised_round${n}.pdb > wt_minimised/${base}_minimised_round${n}_log.txt; done
- foldX
for mutation in "S38R" "N102S" "G202E" "N370S" "L444P"; do base="2V3E_B_${mutation}_foldx" $minimise ../input_noh_nohoh/foldx/${base}.pdb foldx_minimised/${base}_minimised_round1.pdb > foldx_minimised/${base}_minimised_round1_log.txt
- #4 more rounds:
for n in 2 3 4 5; do let p=$n-1; $minimise foldx_minimised/${base}_minimised_round${p}.pdb foldx_minimised/${base}_minimised_round${n}.pdb > foldx_minimised/${base}_minimised_round${n}_log.txt; done done
- SCWRL
- ->did not work with SCWRL models after removing Hs and HETATMs repairPDB! different errors,two last:
- Unable to find atom 1H in ALA
- Could not get GATOM for (ALA11H )
- ran openbabel locally like this:
cd ~/Documents/Studium/Master/SS13/MasterPractical/Assignment9_Structure-based_mutation_analysis/minimise/input_noh_nohoh/scwrl/ for mutation in "S38R" "N102S" "G202E" "N370S" "L444P"; do sudo /home/maria/Documents/Studium/Master/SS13/MasterPractical/Assignment9_Structure-based_mutation_analysis/openbabel/build/bin/obabel -ipdb 2V3E_B_${mutation}_scwrl.pdb -opdb 2V3E_B_${mutation}_swrl_obabel-d.pdb -d > 2V3E_B_${mutation}_scwrl_obabel-d.pdb done
- -> sudo and piping in the same pdb file as given after -opdb, because it only wrote in the file like this.
- -> First mutation created 2 models, separated manually in one model per file like this:
- sed -n '2,7874p' 2V3E_B_S38R_scwrl_obabel-d.pdb > 2V3E_B_S38R_scwrl_obabel-d_model1.pdb
- [sed -n '7877,15749p' 2V3E_B_S38R_scwrl_obabel-d.pdb > 2V3E_B_S38R_scwrl_obabel-d_model2.pdb]
- [checked that repairPDB -jprot (withoutl -noh), then obabel -d makes no difference.
- ran minimise with the first model of S38R:
base="2V3E_B_S38R_scwrl_obabel-d_model1" $minimise ../input_noh_nohoh/scwrl/${base}.pdb scwrl_minimised/${base}_minimised_round1.pdb > scwrl_minimised/${base}_minimised_round1_log.txt
- 4 more rounds:
for n in 2 3 4 5; do let p=$n-1; $minimise scwrl_minimised/${base}_minimised_round${p}.pdb scwrl_minimised/${base}_minimised_round${n}.pdb > scwrl_minimised/${base}_minimised_round${n}_log.txt; done
- ran minimise with the models of the remaining four mutations:
for mutation in "N102S" "G202E" "N370S" "L444P"; do base="2V3E_B_${mutation}_scwrl_obabel-d" $minimise ../input_noh_nohoh/scwrl/${base}.pdb scwrl_minimised/${base}_minimised_round1.pdb > scwrl_minimised/${base}_minimised_round1_log.txt #4 more rounds: for n in 2 3 4 5; do let p=$n-1; $minimise scwrl_minimised/${base}_minimised_round${p}.pdb scwrl_minimised/${base}_minimised_round${n}.pdb > scwrl_minimised/${base}_minimised_round${n}_log.txt; done done </source>