Dear Dr. Pettersen and Dr. Meng,
Thank you for your last email. I'm sorry that I am not so familiar with
Python, so could you help me to use the Chimera to produce the distance
matrix of the listed 45 pdbs below? Otherwise, is there any easy way to do this
just by running the program instead of using the Match->Align by hand? Thank
you so much!
The list of the 45 pdbs:
1il2a,1asza,1kmmc,1adyb,1adjb,1b76b,1efwa,11asb,1j5wb,1bbwa,
1lylb,1e1oa,1pysa,1b8aa,1httb,1kmna,1qe0b,1g51a,1evlc,1evka,
1nj8d,1kmmb,1c0aa,1j5wa,1serb,1h4qa,1eqrc,1sera,1nj1a,1h4sb,
1atib,1h4vb,1b70a,1eqra,1kmnb,1adya,1adjc,1hc7a,1httc,1eqrb,
1g51b,1qe0a,1b76a,1atia,1efwb.
Best regards!
sincerely,
Minming Li
--------------------------
Minming Li
PhD student,
Interdisciplinary Life Science Program (PULSe Program),
Purdue University,
West Lafayette, IN. USA.
Email:
li362@purdue.eduQuoting Eric Pettersen <
pett@cgl.ucsf.edu>:
On Mar 30, 2010, at 9:04 PM, Elaine Meng wrote:
However, there are some issues that may affect whether you decide to
use Chimera to calculate these values for your structures:
(1) You would need to run Match->Align for each pair. There is no
option to do all pairwise comparisons on multiple structures. You
would need to script looping through all pairs (with python or shell
scripting).
(2) There is no Chimera command for Match->Align. Thus Match->Align
would also need to be run via python. I don't have scripting
expertise, so someone else would have to provide more details if you
decide to try this.
Given that you said you wanted a complete pairwise matrix of ~50
structures, that would be ~1200 unique pairs, which is obviously way
too many to do by hand. Currently Match->Align is awkward to script,
even in Python, since all the important code is internal to the Match-
Align dialog and basically impossible to call externally. If you
are sure that you would like to use Chimera for this computation I can
reorganize the Match->Align code so that it could be called
externally. I don't think the reorganization will be difficult. You
would still need to resort to a Python script to make the Match->Align
calls and loop through your models, so you would need some rudimentary
familiarity with Python. Let me know.
--Eric
Eric Pettersen
UCSF Computer Graphics Lab
http://www.cgl.ucsf.edu