SER A
ARG A
GLU A
LYS A
ASN A
GLN A
PRO A
LYS A
PRO A
SER A
PRO A
LYS A
ARG A
####
ATOM 2 CA PRO A 6 2.437 5.295 -23.029 1.00 0.00 C
ATOM 3 C PRO A 6 2.307 5.113 -21.515 1.00 0.00 C
ATOM 4 O PRO A 6 3.276 4.869 -20.823 1.00 0.00 O
ATOM 15 N TRP A 7 1.114 5.225 -20.998 1.00 0.00 N
ATOM 16 CA TRP A 7 0.916 5.054 -19.531 1.00 0.00 C
ATOM 17 C TRP A 7 0.772 3.567 -19.197 1.00 0.00 C <>
####
sub parse_pdb() {
my $file = shift;
my %hash_pdb = ();
# Opening of PDB file after dowloading
open(FH, $file);
@atomrecord = ;
foreach my $record (@atomrecord)
{
chomp($record);
if($record =~/^ATOM/) {
@temp = split(/\s+/, $record);
$con_key = "";
$con_key = "$temp[3]"." "."$temp[4]";
$count = 0;
if( $hash_pdb{$con_key} ) {
if($val_seen{$temp[5]}) { }
else {
$hash_pdb{$con_key} = $hash_pdb{$con_key}."\t"."$temp[5]";
}
}
else {
$hash_pdb{$con_key} = $temp[5];
}
$val_seen{$temp[5]} = 1;
}
if($record =~ /^ENDMDL/){
print "i am going out of loop\n";
last;
}
}
return \%hash_pdb;
}