1
1
# !/usr/bin/perl
2
2
# # Pombert Lab, Illinois Tech, 2021
3
3
my $name = " prepare_visualizations.pl" ;
4
- my $version = " 0.6b " ;
5
- my $updated = " 2022-05-01 " ;
4
+ my $version = " 0.7 " ;
5
+ my $updated = " 2022-05-13 " ;
6
6
7
7
use strict;
8
8
use warnings;
50
50
' m|match=s' => \$match_file ,
51
51
' p|pdb=s' => \$pdb ,
52
52
' r|rcsb=s@{1,}' => \@rcsb ,
53
- ' rlist=s' => \$rcsb_list ,
54
53
' k|keep' => \$keep ,
55
54
' o|out=s' => \$outdir ,
56
55
' l|log=s' => \$log_file ,
75
74
close LOG;
76
75
}
77
76
78
- # # Loading data from RCSB PDB list file
79
- my %rcsb_titles ;
80
- if ($rcsb_list ){
81
- # # Creating a database of RSCB stuctures and their descriptions; PDB 4-letter code => description
82
- open DB, " <" , " $rcsb_list " or die " Can't open tab-delimited file $rcsb_list : $! \n " ;
83
- while (my $line = <DB>){
84
- chomp $line ;
85
- my @columns = split (" \t " , $line );
86
- my $pdb_locus = $columns [0];
87
- my $chain_or_title = $columns [1];
88
- my $description = $columns [2];
89
- $rcsb_titles {$pdb_locus }{$chain_or_title } = $description ;
90
- }
91
- }
92
-
93
77
# # Load predicted pdb filenames into database
94
78
if ($log_file ){
95
79
open LOG," >>" ," $log_file " or die " \n [WARNING]\t Unable to access $log_file : $! \n " ;
149
133
# # Check the PDB headers for proteins that are in the selection provided
150
134
if ($line =~ / ^###/ ){
151
135
my $filename = (fileparse($line ))[0];
152
- ($model_tag ) = $filename =~ / ### (\S +)\; ?/ ;
136
+ if ($aligner eq " gesamt" ){
137
+ ($model_tag ) = $filename =~ / ### (\S +)\; / ;
138
+ }
139
+ elsif ($aligner eq " foldseek" ){
140
+ ($model_tag ) = $filename =~ / ### (\S +)/ ;
141
+ }
142
+
153
143
}
154
144
# # Store the matching RCSB .ent.gz filepath under the locus tag
155
145
else {
@@ -226,4 +216,4 @@ sub Check_Mand_Args{
226
216
die " \n [ERROR]\t GESAMT descriptive match file not provided\n\n $usage \n\n " unless $match_file ;
227
217
die " \n [WARNING]\t RCSB PDB directory(s) not provided, no visualizations will be made" unless @rcsb ;
228
218
die " \n [ERROR]\t Predicted PDB directory not provided\n\n $usage \n\n " unless $pdb ;
229
- }
219
+ }
0 commit comments