Difference between revisions of "Lab Journal - Task 5 (PAH)"

From Bioinformatikpedia
(Modeller)
Line 1: Line 1:
 
== Model calculation ==
 
== Model calculation ==
 
=== Modeller ===
 
=== Modeller ===
To use Modeller we follow the [https://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php/Using_Modeller_for_TASK_4 tutorial] written by the students 2011.
+
To use Modeller we followed the [https://i12r-studfilesrv.informatik.tu-muenchen.de/wiki/index.php/Using_Modeller_for_TASK_4 tutorial] written by the students 2011.
 
*First the target sequence must be provided into PIR-format (2pah has one Lysine at the end and therefore is one aa longer than P00439???):
 
*First the target sequence must be provided into PIR-format (2pah has one Lysine at the end and therefore is one aa longer than P00439???):
 
<code>
 
<code>

Revision as of 16:21, 5 June 2013

Model calculation

Modeller

To use Modeller we followed the tutorial written by the students 2011.

  • First the target sequence must be provided into PIR-format (2pah has one Lysine at the end and therefore is one aa longer than P00439???):

>P1;2pah
sequence:2pah:::::::0.00: 0.00
MSTAVLENPGLGRKLSDFGQETSYIEDNCNQNGAISLIFSLKEEVGALAKVLRLFEENDV
NLTHIESRPSRLKKDEYEFFTHLDKRSLPALTNIIKILRHDIGATVHELSRDKKKDTVPW
FPRTIQELDRFANQILSYGAELDADHPGFKDPVYRARRKQFADIAYNYRHGQPIPRVEYM
EEEKKTWGTVFKTLKSLYKTHACYEYNHIFPLLEKYCGFHEDNIPQLEDVSQFLQTCTGF
RLRPVAGLLSSRDFLGGLAFRVFHCTQYIRHGSKPMYTPEPDICHELLGHVPLFSDRSFA
QFSQEIGLASLGAPDEYIEKLATIYWFTVEFGLCKQGDSIKAYGAGLLSSFGELQYCLSE
KPKLLPLELEKTAIQNYTVTEFQPLYYVAESFNDAKEKVRNFAATIPRPFSVRYDPYTQR
IEVLDNTQQLKILADSINSEIGILCSALQKIK*

An example python script for template 1j8u:
Creating the alignment <source lang=python> from modeller import * env = environ() aln = alignment(env) mdl = model(env, file='1j8u.pdb', model_segment=('FIRST:@', 'END:')) aln.append_model(mdl, align_codes='1j8u', atom_files='1j8u') aln.append(file='target.pir', align_codes='2pah') aln.align2d() aln.check() aln.write(file='/mnt/home/student/waldraffs/masterpracitcal/Task5/alignments/1j8u-2pah-2d.ali', alignment_format='PIR') aln.malign() aln.check() aln.write(file='/mnt/home/student/waldraffs/masterpracitcal/Task5/alignments/1j8u-2pah.ali', alignment_format='PIR') </source> Single template modeling <source lang=python> from modeller import * from modeller.automodel import * log.verbose() env = environ() a = automodel(env,

           alnfile  = '/mnt/home/student/waldraffs/masterpracitcal/Task5/alignments/1j8u-2pah.ali',   
           knowns   = '1j8u',              
           sequence = '2pah',
           assess_methods=(assess.DOPE, assess.GA341))

a.starting_model= 1 a.ending_model = 1 a.make() </source>

Swissmodel

iTasser

Model evaluation