Difference between revisions of "Bioinfromatics HW 1"
imported>CHUNG YEON LEE (Created page with "<p><strong>1. Learn <a href="http://biolecture.org/index.php/Perl" title="Perl">Perl</a> programming language and learn bioinformatics using it.</strong></p> <p>1) <a ...") |
imported>CHUNG YEON LEE |
||
| Line 4: | Line 4: | ||
<p>2) <a href="http://biolecture.org/index.php?title=Pick_5_protein_sequences_and_predict_their_secondary_structures_using_available_prediction_programs&action=edit&redlink=1" title="Pick 5 protein sequences and predict their secondary structures using available prediction programs (page does not exist)">Pick 5 protein sequences and predict their secondary structures using available prediction programs</a></p> | <p>2) <a href="http://biolecture.org/index.php?title=Pick_5_protein_sequences_and_predict_their_secondary_structures_using_available_prediction_programs&action=edit&redlink=1" title="Pick 5 protein sequences and predict their secondary structures using available prediction programs (page does not exist)">Pick 5 protein sequences and predict their secondary structures using available prediction programs</a></p> | ||
| − | |||
| − | |||
<p>#!/bin/perl</p> | <p>#!/bin/perl</p> | ||
| Line 95: | Line 93: | ||
<p> </p> | <p> </p> | ||
| − | <p> | + | <p><img alt="" src="/ckfinder/userfiles/images/1(4).png" style="height:651px; width:1031px" /></p> |
| + | |||
| + | <p><img alt="" src="/ckfinder/userfiles/images/2(1).png" style="height:625px; width:883px" /></p> | ||
<p style="margin-left:5.9pt"><strong>2. Learn various bioinformatics program and tasks</strong></p> | <p style="margin-left:5.9pt"><strong>2. Learn various bioinformatics program and tasks</strong></p> | ||
Latest revision as of 03:42, 16 June 2016
1. Learn Perl programming language and learn bioinformatics using it.
2) Pick 5 protein sequences and predict their secondary structures using available prediction programs
#!/bin/perl
my $dna = 'ATGGCAATGCAGATGCAGCTTGAAGCAAATGCAGATACTTCAGTGGAAGAAGAAAGCTTTGGCCCACAACCCATTTCACGGTTAGAGCAGTGTGGCATAAATGCCAACGATGTGAAGAAATTGGAAGAAGCTGGATTCCATACTGTGGAGGCTGTTGCCTATGCGCCAAAGAAGGAGCTAATAAATATTAAGGGAATTAGTGAAGCCAAAGCTGATAAAATTCTGGCTGAGGCAGCTAAATTAGTTCCAATGGGTTTCACCACTGCAACTGAATTCCACCAAAGGCGGTCAGAGATCATACAGATTACTACTGGCTCCAAAGAGCTTGACAAACTACTTCAAGGTGGAATTGAGACTGGATCTATCACAGAAATGTTTGGAGAATTCCGAACTGGGAAGACCCAGATCTGTCATACGCTAGCTGTCACCTGCCAGCTTCCCATTGACCGGGGTGGAGGTGAAGGAAAGGCCATGTACATTGACACTGAGGGTACCTTTAGGCCAGAACGGCTGCTGGCAGTGGCTGAGAGGTATGGTCTCTCTGGCAGTGATGTCCTGGATAATGTAGCATATGCTCGAGCGTTCAACACAGACCACCAGACCCAGCTCCTTTATCAAGCATCAGCCATGATGGTAGAATCTAGGTATGCACTGCTTATTGTAGACAGTGCCACCGCCCTTTACAGAACAGACTACTCGGGTCGAGGTGAGCTTTCAGCCAGGCAGATGCACTTGGCCAGGTTTCTGCGGATGCTTCTGCGACTCGCTGATGAGTTTGGTGTAGCAGTGGTAATCACTAATCAGGTGGTAGCTCAAGTGGATGGAGCAGCGATGTTTGCTGCTGATCCCAAAAAACCTATTGGAGGAAATATCATCGCCCATGCATCAACAACCAGATTGTATCTGAGGAAAGGAAGAGGGGAAACCAGAATCTGCAAAATCTACGACTCTCCCTGTCTTCCTGAAGCTGAAGCTATGTTCGCCATTAATGCAGATGGAGTGGGAGATGCCAAAGACTGA';
my $protein = '';
my $codon;
sub Codon_table {
my($codon) = @_;
if( $codon =~ /TCA/i ){ return 'S' }
elsif( $codon =~ /TCC/i ){ return 'S' }
elsif( $codon =~ /TCG/i ){ return 'S' }
elsif( $codon =~ /TCT/i ){ return 'S' }
elsif( $codon =~ /TTC/i ){ return 'F' }
elsif( $codon =~ /TTT/i ){ return 'F' }
elsif( $codon =~ /TTA/i ){ return 'L' }
elsif( $codon =~ /TTG/i ){ return 'L' }
elsif( $codon =~ /TAC/i ){ return 'Y' }
elsif( $codon =~ /TAT/i ){ return 'Y' }
elsif( $codon =~ /TAA/i ){ return '*' }
elsif( $codon =~ /TAG/i ){ return '*' }
elsif( $codon =~ /TGC/i ){ return 'C' }
elsif( $codon =~ /TGT/i ){ return 'C' }
elsif( $codon =~ /TGA/i ){ return '*' }
elsif( $codon =~ /TGG/i ){ return 'W' }
elsif( $codon =~ /CTA/i ){ return 'L' }
elsif( $codon =~ /CTC/i ){ return 'L' }
elsif( $codon =~ /CTG/i ){ return 'L' }
elsif( $codon =~ /CTT/i ){ return 'L' }
elsif( $codon =~ /CCA/i ){ return 'P' }
elsif( $codon =~ /CCC/i ){ return 'P' }
elsif( $codon =~ /CCG/i ){ return 'P' }
elsif( $codon =~ /CCT/i ){ return 'P' }
elsif( $codon =~ /CAC/i ){ return 'H' }
elsif( $codon =~ /CAT/i ){ return 'H' }
elsif( $codon =~ /CAA/i ){ return 'Q' }
elsif( $codon =~ /CAG/i ){ return 'Q' }
elsif( $codon =~ /CGA/i ){ return 'R' }
elsif( $codon =~ /CGC/i ){ return 'R' }
elsif( $codon =~ /CGG/i ){ return 'R' }
elsif( $codon =~ /CGT/i ){ return 'R' }
elsif( $codon =~ /ATA/i ){ return 'I' }
elsif( $codon =~ /ATC/i ){ return 'I' }
elsif( $codon =~ /ATT/i ){ return 'I' }
elsif( $codon =~ /ATG/i ){ return 'M' }
elsif( $codon =~ /ACA/i ){ return 'T' }
elsif( $codon =~ /ACC/i ){ return 'T' }
elsif( $codon =~ /ACG/i ){ return 'T' }
elsif( $codon =~ /ACT/i ){ return 'T' }
elsif( $codon =~ /AAC/i ){ return 'N' }
elsif( $codon =~ /AAT/i ){ return 'N' }
elsif( $codon =~ /AAA/i ){ return 'K' }
elsif( $codon =~ /AAG/i ){ return 'K' }
elsif( $codon =~ /AGC/i ){ return 'S' }
elsif( $codon =~ /AGT/i ){ return 'S' }
elsif( $codon =~ /AGA/i ){ return 'R' }
elsif( $codon =~ /AGG/i ){ return 'R' }
elsif( $codon =~ /GTA/i ){ return 'V' }
elsif( $codon =~ /GTC/i ){ return 'V' }
elsif( $codon =~ /GTG/i ){ return 'V' }
elsif( $codon =~ /GTT/i ){ return 'V' }
elsif( $codon =~ /GCA/i ){ return 'A' }
elsif( $codon =~ /GCC/i ){ return 'A' }
elsif( $codon =~ /GCG/i ){ return 'A' }
elsif( $codon =~ /GCT/i ){ return 'A' }
elsif( $codon =~ /GAC/i ){ return 'D' }
elsif( $codon =~ /GAT/i ){ return 'D' }
elsif( $codon =~ /GAA/i ){ return 'E' }
elsif( $codon =~ /GAG/i ){ return 'E' }
elsif( $codon =~ /GGA/i ){ return 'G' }
elsif( $codon =~ /GGC/i ){ return 'G' }
elsif( $codon =~ /GGG/i ){ return 'G' }
elsif( $codon =~ /GGT/i ){ return 'G' }
else{print STDERR "Wrong!!\nPlease check your sequence\n\n";
exit;}
}
for(my $i=0; $i < (length($dna) - 2) ; $i += 3) {
$codon = substr($dna,$i,3);
$protein .= Codon_table($codon);
}
print "\n <DNA>\n\n$dna\n\n <Protein>\n\n$protein\n\n";
exit;
2. Learn various bioinformatics program and tasks
1) Access http://koreanreference.org/home/index.html and search for TERT_HUMAN DNA region
2) Install a protein 3D structure viewer program on your computer
: I chose Pymol
