如何从 pdb 文件中取出某些元素

Question

我正在尝试从 pdb 文件中取出某些列。我已经在我的代码中删除了所有以 ATOM 开头的行。出于某种原因，我的子功能不起作用，我不知道在哪里或如何调用它们。 我的代码是：

open (FILE, $ARGV[0])
    or die "Could not open file\n";

my @newlines;
 while ( my $line = <FILE> ) {
    if ($line =~ m/^ATOM.*/) {
    push @newlines, $line;
    }
}

my $atomcount = @newlines;
#print "@newlines\n";
#print "$atomcount\n";

##############################################################
#This function will take out the element from each line
#The element is from column 77 and contains one or two letters

sub atomfreq {
    foreach my $record1(@newlines) {
      my $element = substr($record1, 76, 2);
      print "$element\n";
      return;
    }
}

################################################################
#This function will take out the residue name from each line
#The element is from column 18 and contains 3 letters

sub resfreq {
    foreach my $record2(@newlines) {
      my $residue = substr($record2, 17, 3);
      print "$residue\n";
      return;
    }
}

Answer 1

正如@Ossip 在this answer 中所说，您只需要调用您的函数：

sub atomfreq {
    ...
}

sub resfreq {
    ...
}

atomfreq();
resfreq();

但是我不确定这些函数是否符合您的预期，因为 cmets 暗示他们应该从 @987654325 打印 每个 $residue 和 $element @ 大批。您已经在for 循环中放置了一个return 语句，它将立即从整个函数（及其for 循环）返回，因此它只会打印第一个$residue 或$element。因为函数不应该返回任何内容，所以您可以删除该语句：

sub atomfreq {
    foreach my $record1(@newlines) {
        my $element = substr($record1, 76, 2);
        print "$element\n";
    }
}

sub resfreq {
    foreach my $record2(@newlines) {
        my $residue = substr($record2, 17, 3);
        print "$residue\n";
    }
}

atomfreq();
resfreq();

Answer 2

您可以直接在其他代码下调用它们，如下所示：

atomfreq();
resfreq();