Difference between revisions of "Workflow homology modelling glucocerebrosidase"
From Bioinformatikpedia
(→MODELLER) |
|||
| Line 3: | Line 3: | ||
== MODELLER == |
== MODELLER == |
||
| − | === |
+ | === Pairwise Sequence Alignments === |
| + | |||
| − | # Save target protein sequence in PIR-format: TARGET.pir |
||
| + | '''1. Preparation of the Alignment File''' |
||
| − | # Save PDB-file of template sequence: TEMPLATE.pdb |
||
| + | * Save target protein sequence in PIR-format: TARGET.pir |
||
| − | #: If PDB-file consists of several chains: split pdb file with the help of [http://structure.usc.edu/splitpdb/ splitpdb] (note that minor changes are needed, so that ATOM coordinates get listed in the resulting PDB-file instead of HETATOMS). |
||
| + | * Save PDB-file of template sequence: TEMPLATE.pdb |
||
| − | # Run the following Python script with command '<code>mod9.9 align.py</code>' to create a target-template alignment in PIR-format: |
||
| + | *: If PDB-file consists of several chains: split pdb file with the help of [http://structure.usc.edu/splitpdb/ splitpdb] (note that minor changes are needed, so that ATOM coordinates get listed in the resulting PDB-file instead of HETATOMS). |
||
| + | * Run the following Python script with command '<code>mod9.9 align.py</code>' to create a target-template alignment in PIR-format: |
||
<code> |
<code> |
||
log.verbose() |
log.verbose() |
||
| Line 19: | Line 21: | ||
aln.write(file='TARGET_TEMPLATE.pap', alignment_format='PAP')</code> |
aln.write(file='TARGET_TEMPLATE.pap', alignment_format='PAP')</code> |
||
| − | + | '''2. Modelling of the Target Structure''' |
|
| − | + | * Run the following Python script with command '<code>mod9.9 model.py</code>' to model the structure of the target sequence: |
|
| − | + | *: Note that all files (alignment- and structure file) must be in the same folder |
|
<code> |
<code> |
||
from modeller.automodel import * |
from modeller.automodel import * |
||
| Line 32: | Line 34: | ||
a.make() |
a.make() |
||
</code> |
</code> |
||
| + | |||
| + | |||
| + | === Multiple Sequence Alignments === |
||
| + | |||
| + | '''1. Preparation of the Alignment File ''' |
||
| + | * Save target protein sequence in PIR-format: TARGET.pir |
||
| + | * Save PDB-files of template sequences: TEMPLATE_x.pdb |
||
| + | * Run the following Python script with command '<code>mod9.9 align_templates.py</code>' to create an alignment of the templates in PIR-format: |
||
| + | <code> |
||
| + | log.verbose() |
||
| + | env = environ() |
||
| + | aln = alignment(env) |
||
| + | for (code, chain) in ((TEMPLATE_1, CHAIN), (TEMPLATE2, CHAIN), ...): |
||
| + | mdl = model(env, file=code, model_segment=('FIRST:'+chain, 'LAST:'+chain)) |
||
| + | aln.append_model(mdl, atom_files=code, align_codes=code+chain) |
||
| + | aln.salign() |
||
| + | aln.write(file='msa.ali', alignment_format='PIR') |
||
| + | aln.write(file='msa.pap', alignment_format='PAP') |
||
| + | </code> |
||
| + | * Run the following Python script with command '<code>mod9.9 align_target.py</code>' to add the TARGET sequence to the multiple sequence alignment: |
||
| + | <code> |
||
| + | from modeller import * |
||
| + | log.verbose() |
||
| + | env = environ() |
||
| + | aln = alignment(env) |
||
| + | aln.append(file='msa.ali', align_codes='all') |
||
| + | aln_block = len(aln) |
||
| + | aln.append(file='TARGET.pir', align_codes='TARGET') |
||
| + | aln.salign() |
||
| + | aln.write(file='msa.ali', alignment_format='PIR') |
||
| + | aln.write(file='msa.pap', alignment_format='PAP') |
||
| + | </code> |
||
| + | |||
| + | '''2. Modelling of the Target Strucutre''' |
||
| + | * Modify the following line of the python script given for the pairwise sequence alignments: |
||
| + | : Original: <code>a = automodel (env, alnfile = 'TARGET_TEMPLATE.ali', knowns = 'TEMPLATE', sequence = 'TARGET')</code> |
||
| + | : Modification: <code>a = automodel(env, alnfile='msa.ali', knowns=('TEMPLATE_1', 'TEMPLATE_2', ...), sequence='TARGET') |
||
== I-TASSER == |
== I-TASSER == |
||
Revision as of 18:00, 12 June 2011
Detailed workflow of the different homology modelling approaches for glucocerebrosidase. Return to overview.
Contents
MODELLER
Pairwise Sequence Alignments
1. Preparation of the Alignment File
- Save target protein sequence in PIR-format: TARGET.pir
- Save PDB-file of template sequence: TEMPLATE.pdb
- If PDB-file consists of several chains: split pdb file with the help of splitpdb (note that minor changes are needed, so that ATOM coordinates get listed in the resulting PDB-file instead of HETATOMS).
- Run the following Python script with command '
mod9.9 align.py' to create a target-template alignment in PIR-format:
log.verbose()
env = environ()
aln = alignment(env)
mdl= model(env, file='TEMPLATE')
aln.append_model(mdl, align_codes='TEMPLATE')
aln.append(file='TARGET.pir', align_codes=('TARGET'))
aln.align(gap_penalties_1d=(-600,-400))
aln.write(file='TARGET_TEMPLATE.ali', alignment_format='PIR')
aln.write(file='TARGET_TEMPLATE.pap', alignment_format='PAP')
2. Modelling of the Target Structure
- Run the following Python script with command '
mod9.9 model.py' to model the structure of the target sequence:- Note that all files (alignment- and structure file) must be in the same folder
from modeller.automodel import *
log.verbose()
env = environ()
env.io.atom_files_directory =
a = automodel (env, alnfile = 'TARGET_TEMPLATE.ali', knowns = 'TEMPLATE', sequence = 'TARGET')
a.starting_model = 1
a.ending_model = 1
a.make()
Multiple Sequence Alignments
1. Preparation of the Alignment File
- Save target protein sequence in PIR-format: TARGET.pir
- Save PDB-files of template sequences: TEMPLATE_x.pdb
- Run the following Python script with command '
mod9.9 align_templates.py' to create an alignment of the templates in PIR-format:
log.verbose()
env = environ()
aln = alignment(env)
for (code, chain) in ((TEMPLATE_1, CHAIN), (TEMPLATE2, CHAIN), ...):
mdl = model(env, file=code, model_segment=('FIRST:'+chain, 'LAST:'+chain))
aln.append_model(mdl, atom_files=code, align_codes=code+chain)
aln.salign()
aln.write(file='msa.ali', alignment_format='PIR')
aln.write(file='msa.pap', alignment_format='PAP')
- Run the following Python script with command '
mod9.9 align_target.py' to add the TARGET sequence to the multiple sequence alignment:
from modeller import *
log.verbose()
env = environ()
aln = alignment(env)
aln.append(file='msa.ali', align_codes='all')
aln_block = len(aln)
aln.append(file='TARGET.pir', align_codes='TARGET')
aln.salign()
aln.write(file='msa.ali', alignment_format='PIR')
aln.write(file='msa.pap', alignment_format='PAP')
2. Modelling of the Target Strucutre
- Modify the following line of the python script given for the pairwise sequence alignments:
- Original:
a = automodel (env, alnfile = 'TARGET_TEMPLATE.ali', knowns = 'TEMPLATE', sequence = 'TARGET') - Modification:
a = automodel(env, alnfile='msa.ali', knowns=('TEMPLATE_1', 'TEMPLATE_2', ...), sequence='TARGET')
I-TASSER
- Webserver: http://zhanglab.ccmb.med.umich.edu/I-TASSER/
- Input: protein sequence in FASTA-format.
I-TASSER provides the possibility to exclude homologous structures with a certain sequence identity cut-off, which was used in this analysis as well.
SWISS-MODEL
Automated Mode
- Webserver: http://swissmodel.expasy.org/workspace/index.php?func=modelling_simple1
- Input: Sequence in fasta format and template PDB-ID (optional)
The automated mode should only be used, if target and template share more than 50% of sequence identity.
Alignment Mode
- Webserver: http://swissmodel.expasy.org/workspace/index.php?func=modelling_align1
- Input: Target-Template Alignment in different formats (FASTA, CLUSTALW, ...)
To create the Alignments needed as input, the tool ClustalW2 was used with standard settings. Additionally the Alignment created with MODELLER was used for 2WNW
