Difference between revisions of "MSUP-task7-journal"

From Bioinformatikpedia
(repairPDB)
(FoldX)
 
(2 intermediate revisions by one other user not shown)
Line 31: Line 31:
   
 
== FoldX ==
 
== FoldX ==
  +
{|class="wikitable" border="1" style="border-spacing:0;text-align: center;"
For FoldX we used the following Files
 
  +
|-
  +
! Mutation
  +
! total energy
  +
! Backbone Hbond
  +
! Sidechain Hbond
  +
! Van der Waals
  +
! Electrostatics
  +
! Solvation Polar
  +
! Solvation Hydrophobic
  +
! Van der Waals clashes
  +
! entropy sidechain
  +
! entropy mainchain
  +
! sloop_entropy
  +
! mloop_entropy
  +
! cis_bond
  +
! torsional clash
  +
! backbone clash
  +
! helix dipole
  +
! water bridge
  +
! disulfide
  +
! electrostatic kon
  +
! partial covalent bonds
  +
! energy Ionisation
  +
! Entropy Complex
  +
|-
  +
| N71S
  +
| 0.47
  +
| 0.13
  +
| 0.13
  +
| 0.12
  +
| 0.00
  +
| -0.16
  +
| 0.21
  +
| 0.00
  +
| -0.16
  +
| 0.20
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| -0.00
  +
| -0.02
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
|-
  +
| M82I
  +
| -0.12
  +
| -0.00
  +
| 0.00
  +
| 0.33
  +
| -0.03
  +
| -0.15
  +
| 0.24
  +
| 0.06
  +
| -0.31
  +
| -0.39
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.12
  +
| 0.04
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
|-
  +
| Q125E
  +
| 1.50
  +
| -0.01
  +
| 3.18
  +
| -0.05
  +
| -0.70
  +
| -0.04
  +
| -0.05
  +
| 0.15
  +
| -1.58
  +
| 0.07
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| -0.03
  +
| -0.03
  +
| 0.55
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
|-
  +
| I213T
  +
| 2.21
  +
| -0.38
  +
| -0.49
  +
| 1.04
  +
| 0.00
  +
| 0.44
  +
| 2.54
  +
| -0.08
  +
| -0.49
  +
| -0.38
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.01
  +
| -0.44
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
|-
  +
| C258Y
  +
| 4.55
  +
| 0.01
  +
| 0.17
  +
| -1.78
  +
| 0.07
  +
| 1.11
  +
| -2.84
  +
| 7.22
  +
| 0.31
  +
| 0.01
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.27
  +
| 0.90
  +
| -0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
|-
  +
| T310R
  +
| 8.51
  +
| 0.02
  +
| 0.28
  +
| -1.67
  +
| 0.13
  +
| 3.01
  +
| -1.29
  +
| 4.89
  +
| 1.30
  +
| 0.04
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 2.22
  +
| 0.11
  +
| -0.42
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
|-
  +
| A328T
  +
| 2.77
  +
| -0.40
  +
| -0.41
  +
| -0.63
  +
| -0.01
  +
| 1.14
  +
| -0.63
  +
| 2.76
  +
| 0.51
  +
| -0.92
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 1.36
  +
| -0.01
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
|-
  +
| I361V
  +
| 0.50
  +
| 0.01
  +
| -0.01
  +
| 0.48
  +
| 0.00
  +
| -0.28
  +
| 0.86
  +
| -0.04
  +
| -0.55
  +
| 0.02
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.01
  +
| -0.07
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| -0.01
  +
| 0.00
  +
|-
  +
| N404S
  +
| 0.20
  +
| 0.01
  +
| 0.20
  +
| 0.37
  +
| 0.00
  +
| -0.48
  +
| 0.50
  +
| -0.02
  +
| -0.25
  +
| -0.22
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.09
  +
| -0.11
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
|-
  +
| R429H
  +
| -0.08
  +
| -0.14
  +
| -0.19
  +
| 0.37
  +
| 0.14
  +
| -0.50
  +
| 0.40
  +
| 0.03
  +
| -0.43
  +
| 0.26
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| -0.19
  +
| -0.20
  +
| -0.07
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.00
  +
| 0.25
  +
| 0.00
  +
|-
  +
|}
   
 
== Minimise ==
 
== Minimise ==
Line 38: Line 304:
 
=== repairPDB ===
 
=== repairPDB ===
   
USAGE: repairPDB <PDB FILE> [<RANGE>] [<CHAIN:] [<CHAIN:RANGE>] [OPTIONS]
+
USAGE: repairPDB <PDB FILE> [<RANGE>] [<CHAIN:] [<CHAIN:RANGE>] [OPTIONS]
OPTIONS: [-offset value] offset the residue numbering
+
OPTIONS: [-offset value] offset the residue numbering
 
[-chain char] change Chain ID
 
[-chain char] change Chain ID
 
[-ratom] renumber Atoms
 
[-ratom] renumber Atoms
Line 47: Line 313:
 
[-seq] protein sequence from AA
 
[-seq] protein sequence from AA
 
[-seqrs] protein sequence from SEQRES entries
 
[-seqrs] protein sequence from SEQRES entries
  +
 
 
[-jprot] just Protein OR
 
[-jprot] just Protein OR
 
[-nohoh] no water OR
 
[-nohoh] no water OR
 
[-ssw cutoff] print only waters with B-value below cutoff OR
 
[-ssw cutoff] print only waters with B-value below cutoff OR
 
[-cleansol] remove overlapping solvent for GROMACS
 
[-cleansol] remove overlapping solvent for GROMACS
  +
 
 
[-dna] print DNA only OR
 
[-dna] print DNA only OR
 
[-nodna] do not print DNA
 
[-nodna] do not print DNA
  +
 
=== SCWRL ===
 
=== SCWRL ===
 
resulting sequence:
 
resulting sequence:

Latest revision as of 20:45, 25 June 2012

Scwrl

Sequence extraction from pdb-file:

/opt/SS12-Practical/scripts/repairPDB 2bfd_A.ent -chain A -seq
sequence mutated:
kpqfpgasaefidklefiqpSvisgipiyrvIdrqgqiinpsedphlpkekvlklyksmtll
ntmdrilyesqrEgrisfymtnygegtvgsaaaldntdlvfgareagvlmyrdyplelfmaq
cygnisdlgkgrqpvygckerhfvtissplatqTpqavgaayaakrananrvvicyfgegaa
segdahagfnfaatleYpiiffcrnngyaistptseqyrgdgiaargpgygimsirvdgndv
favynaRkearrravaenqpflieTmtyrigstdhpVsrlrhylqgwwdeeqekawrkqsrr
kvmeafeqaerkpkpSpnllfsdvyqempaqlrkqqeslaHhlqtygehypldhfdk

posmap:

snp-pospdb-pos-mapping
N71SN21S
M82IM32I
Q125EQ75E
I213TI158T
C258YC203Y
T310RT255R
A328TA273T
I361VI285V
N404SN326S
R429HR351H

Commad for SCWRL4:

/opt/SS12-Practical/scwrl4/Scwrl4 -i 2bfd_A.ent -s BCKDHA_2bfd_a.sequence_final -o model_OUT.pdb

FoldX

Mutation total energy Backbone Hbond Sidechain Hbond Van der Waals Electrostatics Solvation Polar Solvation Hydrophobic Van der Waals clashes entropy sidechain entropy mainchain sloop_entropy mloop_entropy cis_bond torsional clash backbone clash helix dipole water bridge disulfide electrostatic kon partial covalent bonds energy Ionisation Entropy Complex
N71S 0.47 0.13 0.13 0.12 0.00 -0.16 0.21 0.00 -0.16 0.20 0.00 0.00 0.00 -0.00 -0.02 0.00 0.00 0.00 0.00 0.00 0.00 0.00
M82I -0.12 -0.00 0.00 0.33 -0.03 -0.15 0.24 0.06 -0.31 -0.39 0.00 0.00 0.00 0.12 0.04 0.00 0.00 0.00 0.00 0.00 0.00 0.00
Q125E 1.50 -0.01 3.18 -0.05 -0.70 -0.04 -0.05 0.15 -1.58 0.07 0.00 0.00 0.00 -0.03 -0.03 0.55 0.00 0.00 0.00 0.00 0.00 0.00
I213T 2.21 -0.38 -0.49 1.04 0.00 0.44 2.54 -0.08 -0.49 -0.38 0.00 0.00 0.00 0.01 -0.44 0.00 0.00 0.00 0.00 0.00 0.00 0.00
C258Y 4.55 0.01 0.17 -1.78 0.07 1.11 -2.84 7.22 0.31 0.01 0.00 0.00 0.00 0.27 0.90 -0.00 0.00 0.00 0.00 0.00 0.00 0.00
T310R 8.51 0.02 0.28 -1.67 0.13 3.01 -1.29 4.89 1.30 0.04 0.00 0.00 0.00 2.22 0.11 -0.42 0.00 0.00 0.00 0.00 0.00 0.00
A328T 2.77 -0.40 -0.41 -0.63 -0.01 1.14 -0.63 2.76 0.51 -0.92 0.00 0.00 0.00 1.36 -0.01 0.00 0.00 0.00 0.00 0.00 0.00 0.00
I361V 0.50 0.01 -0.01 0.48 0.00 -0.28 0.86 -0.04 -0.55 0.02 0.00 0.00 0.00 0.01 -0.07 0.00 0.00 0.00 0.00 0.00 -0.01 0.00
N404S 0.20 0.01 0.20 0.37 0.00 -0.48 0.50 -0.02 -0.25 -0.22 0.00 0.00 0.00 0.09 -0.11 0.00 0.00 0.00 0.00 0.00 0.00 0.00
R429H -0.08 -0.14 -0.19 0.37 0.14 -0.50 0.40 0.03 -0.43 0.26 0.00 0.00 0.00 -0.19 -0.20 -0.07 0.00 0.00 0.00 0.00 0.25 0.00

Minimise

For minimise we used a slightly adapted version of the Task7_Hemochromatosis_Protocol#Minimise shell script from the heterochromatisis group.

Gromacs

repairPDB

USAGE:   repairPDB <PDB FILE> [<RANGE>] [<CHAIN:] [<CHAIN:RANGE>] [OPTIONS]
OPTIONS: [-offset value] offset the residue numbering
         [-chain char]   change Chain ID
         [-ratom]        renumber Atoms
         [-rres]         renumber Residues
         [-noh]          remove hydrogens
         [-het]          do not change HETATM to ATOM for AA
         [-seq]          protein sequence from AA
         [-seqrs]        protein sequence from SEQRES entries

         [-jprot]        just Protein OR
         [-nohoh]        no water OR
         [-ssw cutoff]   print only waters with B-value below cutoff OR
         [-cleansol]     remove overlapping solvent for GROMACS

         [-dna]          print DNA only OR
         [-nodna]        do not print DNA

SCWRL

resulting sequence:

sslddkpqfpgasaefidklefiqpnvisgipiyrvmdrqgqiinpsedphlpkekvlklyksmtllntmdril
yesqrqgrisfymtnygeegthvgsaaaldntdlvfgqyreagvlmyrdyplelfmaqcygnisdlgkgrqmpv
hygckerhfvtissplatqipqavgaayaakrananrvvicyfgegaasegdahagfnfaatlecpiiffcrnn
gyaistptseqyrgdgiaargpgygimsirvdgndvfavynatkearrravaenqpflieamtyrighhstsdd
ssayrsvdevnywdkqdhpisrlrhyllsqgwwdeeqekawrkqsrrkvmeafeqaerkpkpnpnllfsdvyqe
mpaqlrkqqeslarhlqtygehypldhfdk