<br>#!/usr/local/bin/perl<br>use LWP::Simple;<br>use strict;<br>use warnings;<br>use LWP 5.64;<br>use LWP::UserAgent;<br>print "Quale entry di PDB vuoi aprire? ";<br>my $pdb=<STDIN>;<br>chomp ($pdb);<br>my $url1 = "<a href="http://www.ebi.ac.uk/thornton-srv/databases/cgi-bin/pdbsum/GetPage.pl?pdbcode=$pdb" target="_blank">http://www.ebi.ac.uk/thornton-srv/databases/cgi-bin/pdbsum/GetPage.pl?pdbcode=$pdb</a>";<br>
my $content = get $url1;<br>die "Couldn't get $url1" unless defined $content;<br>while ($content=~ m/uniprot_id=(.*?)&/g)<br><br> {<br> my $code=$1;<br> print "ID uniprot: $code\n"; <br>
chomp ($code);<br> my $browser =LWP::UserAgent->new;<br> my $url = URI-> new ("<a href="http://www.uniprot.org/uniprot/$code.fasta" target="_blank">http://www.uniprot.org/uniprot/$code.fasta</a>");<br>
my $response=$browser-> get ($url);<br>
my $fasta=$response->content;<br> print "$fasta" ;<br> };<br> <br><br><br><br><br>this
is my script...first I must digit the pdb code (form my test I use
1xmj, that I'm sure it has a ligand) I want to explore, (that's the
code of a protein structure) and then I go to the web page...to extract
another code, that'uniprot code, so the script thakes that code and put
it into another web site and estract and print the aminoacidic
sequence of the protein, then i need to know something about the
ligand...and I need to access to the file running.pl that brings some
information which I have to extract....and i don't know what to do<br>