Gaucher Task07 Protocol
From Bioinformatikpedia
SCWRL
scwrl -i 2nt0_A.pdb -s 2nt0_A_1.seq -o 2nt0_A_1.pdb > 2nt0_A_1.out
FoldX
list.txt
2nt0_A.pdb
mutant_list.txt:
ARPCIPKSFGYSSVVCVCNATYCDSFDPPTFPALGTFSRYESTRSGRRMELSMGPIQANHTGTGLLLTLQPEQKFQKVKG FGGAMTDAAALNILALSPPAQNLLLKSYFSEEGIGYNIIRVPMASCDFSIRTYTYADTPDDFQLHNFSLPEEDTKLKIPL IHRALQLAQRPVSLLASPWTSPTWLKTNGAVNGKGSLKGQPGDIYHQTWARYFVKFLDAYAEHKLQFWAVTAENEPSAGL LSGYPFQCLGFTPEHQRDFIARDLGPTLANSTHHNVRLLMLDDQRLLLPHWAKVVLTDPEAAKYVHGIAVHWYLDFLAPA KATLGETHRLFPNTMLFASEACVGSKFWEQSVRLGSWDRGMQYSHSIITNLLYHVVGWTDWNLALNPEGGPNWVRNFVDS PIIVDITKDTFYKQPMFYHLGHFSKFIPEGSQRVGLVASQKNDLDAVALMHPDGSAVVVVLNRSSKDVPLTIKDPAVGFL ETISPGYSIHTYLWHRQ ARPCIPKSFGYSSVVCVCNATYCDSFDPPTFPALGTFSRYESTRSGRRMELSMGPIQANRTGTGLLLTLQPEQKFQKVKG FGGAMTDAAALNILALSPPAQNLLLKSYFSEEGIGYNIIRVPMASCDFSIRTYTYADTPDDFQLHNFSLPEEDTKLKIPL IHRALQLAQRPVSLLASPWTSPTWLKTNGAVNGKGSLKGQPGDIYHQTWARYFVKFLDAYAEHKLQFWAVTAENEPSAGL LSGYPFQCLGFTPEHQRDFIARDLGPTLANSTHHNVRLLMLDDQRLLLPHWAKVVLTDPEAAKYVHGIAVHWYLDFLAPA KATLGETHRLFPNTMLFASEACVGSKFWEQSVRLGSWDRGMQYSHSIITNLLYHVVGWTDWNLALNPEGGPNWVRNFVDS PIIVDITKDTFYKQPMFYHLGHFSKFIPEGSQRVGLVASQKNDLDAVALMHPDGSAVVVVLNRSSKDVPLTIKDPAVGFL ETISPGYSIHTYLWHRQ ...
run.txt
2nt0_A.pdb <TITLE>FOLDX_runscript; <JOBSTART>#; <PDBS>#; <BATCH>list.txt; <COMMANDS>FOLDX_commandfile; <BuildModel>#,mutant_file.txt; <END>#; <OPTIONS>FOLDX_optionfile; <Temperature>298; <R>#; <pH>4.5; <IonStrength>0.050; <water>-CRYSTAL; <metal>-CRYSTAL; <VdWDesign>2; <OutPDB>true; <pdb_hydrogens>false; <complex_with_DNA> true; <END>#; <JOBEND>#; <ENDFILE>#;
foldx -runfile run.txt
Gromacs
repairPDB:
/opt/SS12-Practical/scripts/repairPDB 2nt0_A_1.pdb -jprot > output/2nt0_A_1_re.pdb
SCWRL:
/opt/SS12-Practical/scwrl4/Scwrl4 -i output/2nt0_A_1_re.pdb -o output/2nt0_A_1_re_sc.pdb
repairPDB again:
/opt/SS12-Practical/scripts/repairPDB output/2nt0_A_1_re_sc.pdb -jprot > output/2nt0_A_1.pdb
pdb2gmx:
/opt/SS12-Practical/gromacs/bin/pdb2gmx -f output/2nt0_A_1.pdb -o output/2nt0_A_1.gro -p output/2nt0_A_1.top -water tip3p -ff amber03
Create a MDP file with the following content:
title = PBSA minimization in vacuum cpp = /usr/bin/cpp # Preprocessing # -DFLEXIBLE: include flexible water in stead of rigid water into topology # -DPOSRES: includes posre.itp into topology, used for position restraints. define = -DFLEXIBLE -DPOSRES # Implicit solvent - simulation with implicit solvent using the Generalized Born formalism implicit_solvent = GBSA # Run control - steepest descent algorithm for energy minimization integrator = steep # Run control - tolerance emtol = 1.0 # Run control - maximum number of steps to integrate or minimize nsteps = 500 # Output - frequency to write energies to energy file nstenergy = 1 # Tables - group(s) to write to energy file energygrps = System # Neighbour searching - makes grid in the box and only check atoms in neighboring grid cells when constructing a new neighbor list every nstlist steps ns_type = grid # Electrostatics - twin range cut-off's with neighborlist cut-off rlist and Coulomb cut-off rcoulomb coulombtype = cut-off rcoulomb = 1.0 # VDW - distance for the LJ or Buckingham cut-off rvdw = 1.0 # Bonds - no constraints except for those defined explicitly in the topology constraints = none # Neighbour searching - Use no periodic boundary conditions, ignore the box. pbc = no