Le Monday 25 June 2012 à 03:41:56PM, Modeller Caretaker a écrit :
> On 6/25/12 10:33 AM, Thomas Gaillard wrote:
> >I have noticed that some HETATM records are read by Modeller and
> >automatically handled as corresponding standard residues (for example
> >MSE -> MET)
>
> Actually MSE -> MET is the only such translation done.
OK, I was saying "some HETATM records" because of the following
paragraph in MODELLER manual, also mentioning MEX and ABU:
http://salilab.org/modeller/9.9/manual/node171.html
Some special handling is done for certain commonly-used modified
residues. The MSE residue type is mapped to the regular MET amino acid,
and the SE atom in this residue is mapped to SD. The MEX and ABU residue
types are both mapped to CYS. In all cases, these residues are read as
regular amino acids (ATOM records) regardless of whether they are marked
as ATOM or HETATM in the PDB file. There is no need to turn on
env.io.hetatm.
> >I would like to know if you could suggest a way to disable this
> >behaviour and really ignore all HETATM records.
>
> Unfortunately there isn't any way to do that - it is hardcoded into
> the PDB reader. Why do you want to treat MSE as HETATMs anyway?
> There may be a better way to achieve your goal.
Because when MODELLER is integrated with other programs, some do not
make the same choice of automatically treating HETATM MSE as ATOM MET
and this can lead to confusion. I think this behaviour should not be
hardcoded.
> Ben Webb, Modeller Caretaker
Additionnally, I think your mapping of MSE SE atom to MET SD atom is slightly
incorrect in the eyes of the PDB standard:
HETATM XXXX SE MSE ... SE
should be mapped to:
ATOM XXXX SD MET ... S
and not to:
ATOM XXXX SD MET ... SE
The atom name format and element name are incorrect.
--
Thomas Gaillard
Maître de conférences
Laboratoire de Biochimie
Ecole Polytechnique
91128 Palaiseau cedex
tel: +33 1 69 33 48 62
fax: +33 1 69 33 49 09
http://bioc.polytechnique.fr/~gaillard