Mercurial > repos > cpt > cpt_psm_recombine

comparison lib/CPT/Bio/SAR.pm @ 1:97ef96676b48 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload commit 94b0cd1fff0826c6db3e7dc0c91c0c5a8be8bb0c
author cpt
date Mon, 05 Jun 2023 02:51:26 +0000
parents
children
comparison
equal deleted inserted replaced
0:b18e8268bf4e 1:97ef96676b48
1 package CPT::Bio::SAR;
2 use strict;
3 use warnings;
4 use autodie;
5 use Moose;
6
7 sub filter_sar {
8 my ($self, @seqs) = @_;
9 my @good;
10 foreach(@seqs){
11 if(has_sar_motif($_)){
12 push(@good, $_);
13 }
14 }
15 return @good;
16 }
17
18 sub has_sar_motif {
19 my ( $self, $seq ) = @_;
20
21 return 0 if(length $seq < 40);
22
23 my $reg_a = qr/([^DEKR]{3}K[^DEKR]{8,}[^DER]{1}[^DEKR]{3})/;
24 my $reg_b = qr/([KR]{1,}[^DEKR]{12,}[^DER]{1}[^DEKR]{3})/;
25
26 my $first40 = substr( $seq, 0, 40 );
27
28 # there is a transmembrane domain in the first 40 AAs
29 # there is at least one positive charged AAs in front of the TMD
30 if ( $first40 =~ $reg_a || $first40 =~ $reg_b ) {
31 my $modi1st40 = $first40;
32 # Cut out the match, and then add the whole thing to the end.
33 my $t4homology =
34 #substr($seq,0, $-[0] ), # Before the match
35 #substr($seq, $-[0], ($+[0] - $-[0])), # the match
36 substr($seq, $+[0]). # After the match
37 $first40;
38 $t4homology = substr($t4homology, 0 , 40);
39
40 if ( $t4homology =~ qr/E[A-Z]{8}[DC][A-Z]{4,5}T/ ) {
41 return 1;
42 }
43 }
44 return 0;
45 }
46
47 no Moose;
48 1;
49
50 __END__
51
52 =pod
53
54 =encoding UTF-8
55
56 =head1 NAME
57
58 CPT::Bio::SAR
59
60 =head1 VERSION
61
62 version 1.99.4
63
64 =head1 AUTHOR
65
66 Eric Rasche <rasche.eric@yandex.ru>
67
68 =head1 COPYRIGHT AND LICENSE
69
70 This software is Copyright (c) 2014 by Eric Rasche.
71
72 This is free software, licensed under:
73
74 The GNU General Public License, Version 3, June 2007
75
76 =cut