Difference between revisions of "Structure-based mutation analysis"
m (→FoldX) |
m (→SCWRL) |
||
Line 33: | Line 33: | ||
SCWRL is able to predict side-chain conformations by using a backbone given from an experimental solved structure. |
SCWRL is able to predict side-chain conformations by using a backbone given from an experimental solved structure. |
||
− | Usage: |
+ | '''Usage:''' |
* use only chain A of backbone pdb: <code>1A6ZA.pdb</code> |
* use only chain A of backbone pdb: <code>1A6ZA.pdb</code> |
||
* extract amino acid sequence and change it to lowercase: <code>aa.txt</code> |
* extract amino acid sequence and change it to lowercase: <code>aa.txt</code> |
||
Line 39: | Line 39: | ||
* cmd: <code>scwrl4 -i 1A6ZA.pdb -s aa_x.txt -o 1A6ZA_mutant.pdb</code> |
* cmd: <code>scwrl4 -i 1A6ZA.pdb -s aa_x.txt -o 1A6ZA_mutant.pdb</code> |
||
− | Results: |
+ | '''Results:''' |
to be done |
to be done |
||
Revision as of 14:38, 30 June 2011
General
According to the UniProt entry about HFE_HUMAN are three 3D-structures of the HFE_HUMAN available, which are listed below. We have chosen the '1A6Z' because it has the best resolution, a very good R-Value (it measures the quality of the model obtained from the crystallographic data), a pH near the physiological optimum and is as good as complete. '1DE4' has a slightly better R-Value and pH, but this PDB also includes the transferrin receptor, which we do not need and do not want in our structure. '1C42' is only a hypothetical model, so we exclude it from further research.
All stereochemistrical properties of the structure are shown in the figure to the right<ref>Lebrón JA, Bennett MJ, Vaughn DE, Chirino AJ, Snow PM, Mintier GA, Feder JN, Bjorkman PJ.: Crystal structure of the hemochromatosis protein HFE and characterization of its interaction with transferrin receptor.</ref>.
PDB | Method | Resolution (Å) | Chain | R-Value | R-Free | pH |
---|---|---|---|---|---|---|
1A6Z | X-ray | 2.60 | A/C | 0.233 | 0.277 | 6.5 |
1C42 | model | - | A | - | - | - |
1DE4 | X-ray | 2.80 | A/D/G | 0.231 | 0.265 | 8.0 |
Mapping
Because we have no annotation about the active site so we just visualized the mutations at the '1A6Z' structure. Secondly we are using the same mutations as in Task 6 and have therefore the same problemes with their visualization (only 7 of 10 visualizeable).
All mutations are scattered accross the protein with no affinity to some special region or secondary structure element. Mutations shown in red are near glycosylation positions and mutations in yellow are near disulfidbonds.
SCWRL
SCWRL is able to predict side-chain conformations by using a backbone given from an experimental solved structure.
Usage:
- use only chain A of backbone pdb:
1A6ZA.pdb
- extract amino acid sequence and change it to lowercase:
aa.txt
- introduce each mutation into on seperated
aa_x.txt
file as capital - cmd:
scwrl4 -i 1A6ZA.pdb -s aa_x.txt -o 1A6ZA_mutant.pdb
Results: to be done
Energy comparison
FoldX
FoldX scores the importance of amino acid interactions according to the overall stability of the protein and calculates the energy.
Usage:
- create a runfile tutorial and adjust all default parameters to known (if possible):
runfile.txt
- create a listfile of all pdb files that should be included in energy calculation:
listfile.txt
- cmd:
sudo ./foldx -runfile runfile.txt > output.txt
Results: to be done.
References
<references />