Dear folks,

I have a couple of questions about the alignment and fitting functions
of pymol.

I've been using the 'align' and 'fit' functions of pymol apparently
quite successfully, but I found it difficult to make the 'rms' or
'rms_cur' functions actually report anything. I can take a structure
with identical sequence in two conformations and 'fit' it properly or I
can use 'align' and get an RMS value, but 'rms' doesn't report anything
directly. Now I did find that I could get the value with:

  print cmd.rms('selection1','selection2')

but I was surprised that I couldn't get it directly.  After looking in
the code, I found I could do:

rms (selection1),(selection2),quiet=0

and then the value was reported.  Should that not be the default?

Second question: Is it possible to get the identity of the the atoms
used by 'align'?  I know from 'help align' that 

  "If object is not None, then align will create an object which
   indicates which atoms were paired between the two structures"

but could we possibly have access to a list of atoms or residues
matched?

Cheers,
Robert
-- 
Robert L. Campbell, Ph.D.                         <r...@post.queensu.ca>
Senior Research Associate                            phone: 613-533-6821
Dept. of Biochemistry, Queen's University,             fax: 613-633-2497
Kingston, ON K7L 3N6  Canada
    PGP Fingerprint: 9B49 3D3F A489 05DC B35C  8E33 F238 A8F5 F635 C0E2
         (out of date web site:http://biophysics.med.jhmi.edu/rlc)

Reply via email to