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

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
EEEKKTWGTVFKTLKSLYKTHACYEYNHIFPLLEKYCGFHEDNIPQLEDVSQFLQTCTGF
RLRPVAGLLSSRDFLGGLAFRVFHCTQYIRHGSKPMYTPEPDICHELLGHVPLFSDRSFA
QFSQEIGLASLGAPDEYIEKLATIYWFTVEFGLCKQGDSIKAYGAGLLSSFGELQYCLSE
KPKLLPLELEKTAIQNYTVTEFQPLYYVAESFNDAKEKVRNFAATIPRPFSVRYDPYTQR

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,

knowns   = '1j8u',
sequence = '2pah',
assess_methods=(assess.DOPE, assess.GA341))

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

Also done for 2phm and 3luy.

I-TASSER

We used the I-TASSER server for the protein structure and function prediction.

First, we set the fasta sequence of our protein PAH (P00439) as query sequence into the big field on the beginning. Then, we have three different templates (1J8U, 2PHm and 3LUY) given and so we have to run three different jobs with I-TASSER. Therefore, we type the PDB ID into the field "Specify template without alignment" in Option I. Since the protein 1J8U has a 100% identity to our protein (2PAH), we did not set any cutoff value for this job. However, we set the cutoff for the PDB ID of 2PHM to 90% and for 3LUY to 35%. For this purpose, one has to set the cutoff in percent in Option II into the field "Exclude homologous templates". Subsequent, one can run the job by clicking onto the button "Run I-TASSER. Since the server is overloaded with jobs, one run takes about 60 hours for completion and everyone can only submit one job.