[Chimera-users] Match command question
meng at cgl.ucsf.edu
Fri Nov 17 10:39:27 PST 2017
I don’t understand the situation well. So each state (unliganded, and in complex with each different ligand) is a tetramer? And there is one pivot point in each monomer of the tetramer? I.e. something like chains A-D, with :430-433 occurring in each chain?
Your match command is trying to match all the pivot points of all four tetramers at once with a single rigid-body movement. It won’t move the monomers separately. To do that, you’d have to split each tetramer into separate models for each monomer and then match the monomers individually.
However, if you do intend to match all the pivot points of all four tetramers at once with a single rigid-body movement, my only idea is that maybe the chain IDs or ordering are different for some of the liganded tetramers than the others. In that case you may have to give a more elaborate atom specification in the commandor edit the files (e.g. change chain IDs) to get the correct pairing. You can mouse over the structures to see if the chain IDs are for the corresponding chains (using the balloon information that appears when you pause over the structure). This is really just a guess but all I can think of based on the description.
I hope this helps,
Elaine C. Meng, Ph.D.
UCSF Chimera(X) team
Department of Pharmaceutical Chemistry
University of California, San Francisco
> On Nov 17, 2017, at 2:31 AM, Yann STERCKX <Yann.Sterckx at vub.be> wrote:
> Dear Chimera users,
> I've been struggling with the match command and was hoping you could help me out.
> I'm working a tetrameric enzyme. When the enzyme binds ligands, its domains rotate over a pivot point in the monomer. What I would like to do is superpose different liganded structures of the enzyme on the four monomer pivot points. Here's what I've done so far:
> * Each liganded state of the tetrameric enzyme is a single model
> * I use the match command to superpose the structures on the monomer pivot points (residues 430-433 in the structure). I used the following command
> match #model1:430-433 #model2:430-433
> What is strange is that I obtain exactly what I want for 2 out of 4 models. Using this command I get an rmsd of 0.77A. But for the other 2 models, match does not seem to work as I obtain an RMSD of 27A and the superposition is clearly off. Including the iterate option does not alleviate the problem. I don't want to use matchmaker because this feature will, at least in my understanding, try to optimally fit two chains (or selections within two chains) on top of each other. If I do that, I do not get to see how the enzyme rotates around the monomer pivot points.
> How come match works really well for 2 out of 4 structures, but fails to perform the same task for the rest? Any pointers on what I might be doing wrong?
> Many thanks and kind regards,
> PS: I've browsed through the archive, but could not really find a thread related to this problem (unless I overlooked it).
> Chimera-users mailing list: Chimera-users at cgl.ucsf.edu
> Manage subscription: http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users
More information about the Chimera-users