[phenixbb] refine : peculiar TLS tensors

Alexandre OURJOUMTSEV sacha at igbmc.fr
Wed Oct 20 00:32:14 PDT 2010


Hi, Bryan,

Do not going into details of TLS, the very first remark is that the center of reaction that characterizes the MOVEMENT of the group is NOT AT ALL the center of mass that characterizes the GROUP as it is independently of any its movement. They should not be confused.

Sacha

-----Message d'origine-----
De : phenixbb-bounces at phenix-online.org [mailto:phenixbb-bounces at phenix-online.org] De la part de Bryan Lepore
Envoyé : mercredi 20 octobre 2010 00:01
À : phenixbb at phenix-online.org
Objet : [phenixbb] refine : peculiar TLS tensors

[ phenix 1.6.4-486 ]
[ ALSO : tlsanl, refmac, etc. from ccp4 6.1.3 ]
[ molscript/povscript ]

the one-liner : i have some TLS tensors that sometimes do not compare
at all well with refmac/anisoanl tensors.

sometimes, the principal axes of a TLS tensor refines *very* far away
from the center of mass(=center of reaction), or the axes from tlsanl
are exceptionally small as if not scaled properly. the eigenvalues are
really small too. see [*] for troubleshooting and strategy used.

i suspect i am missing something. i also understand this might not be
a phenix question.

-bryan

[*] i compared another .def that gave sensible TLS results with a
different dataset/model, and my {} are all good AFAIK.
anisoanl/refmac/parvati results are reasonable. the phenix refinement
appears well behaved. these are all using tlsextract to get the
eigenvalues/origins, and molscript to check the axes/centers of mass.
input B factors are Beq unless i made an error. the tls group
selections can vary - however, simple groups also can give these
results.

in general, using strategy (excerpt) :

  refine {
    strategy = *individual_sites individual_sites_real_space rigid_body \
               *individual_adp group_adp *tls *occupancies *group_anomalous
    sites {
      individual = None
      torsion_angles = None
      rigid_body = None
    }
    adp {
      individual {
        isotropic = None
        anisotropic = None
      }
      group_adp_refinement_mode = *one_adp_group_per_residue \
                                  two_adp_groups_per_residue group_selection
      group = None
      tls = "chain A and (resid 1:41 or resid 42:159 or resid 195:237
or resid 263:420)"
      tls = "chain A and resid 160:194"
      tls = "chain A and resid 245:252"
    }
[...cut off ... ]

example output of an ill-behaved S tensor from tlsanl :

 SCREW ROTATION AXIS        ABSOLUTE POSITION (A)     PITCH (A)
                               X       Y       Z
   1                         70.667 154.197 308.627    -144.435
   2                        977.334 819.588********    6222.947
   3                        640.7191431.084 471.341     269.382
_______________________________________________
phenixbb mailing list
phenixbb at phenix-online.org
http://phenix-online.org/mailman/listinfo/phenixbb



More information about the phenixbb mailing list