#!/usr/bin/perl

%amino_acids= (
 'A' => 'Alanine',
 'C' => 'Cysteine',
 'D' => 'Aspartate',
 'E' => 'Glutamate',
 'F' => 'Phenylalanine',
 'G' => 'Glycine',
 'H' => 'Histidine',
 'I' => 'Isoleucine',
 'K' => 'Lysine',
 'L' => 'Leucine',
 'M' => 'Methionine',
 'N' => 'Asparagine',
 'P' => 'Proline',
 'Q' => 'Glutamine',
 'R' => 'Arginine',
 'S' => 'Serine',
 'T' => 'Threonine',
 'V' => 'Valine',
 'W' => 'Tryptophane',
 'Y' => 'Tyrosine'
);

%aas=(
ALA => 'A',
CYS => 'C',
ASP => 'D',
GLU => 'E',
PHE => 'F',
GLY => 'G',
HIS => 'H',
ILE => 'I',
LYS => 'K',
LEU => 'L',
MET => 'M',
ASN => 'N',
PRO => 'P',
GLN => 'Q',
ARG => 'R',
SER => 'S',
THR => 'T',
VAL => 'V',
TRP => 'W',
TYR => 'Y'
);


if (scalar @ARGV == 0) {
  print "il manque un nom de fichier PDB\n";
  exit(1);
}

$file=$ARGV[0];
print "ouverture de $file\n";

open(FICHIER,"<".$file) or die("impossible d'ouvrir le fichier $file");