<div dir="ltr">Dear Elaine,<div><br></div><div>Just to check if I understood you correctly. I can only superimpose model #1 to #2 in one pairwise alignment job. Then with a separate pairwise alignment job, I can superimpose #1 to #3. I cannot superimpose more than two models simultaneously. Is my understanding correct? Thank you.</div><div><br></div><div>Best regards,</div><div>Zichen "Cardiff" Jiang</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Mar 28, 2021 at 6:24 PM Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu">meng@cgl.ucsf.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi Cardiff,<br>
You can only do pairwise alignments. There is no way to best-fit one model to multiple others at the same time.<br>
<br>
However, you can still superimpose all of the models together by picking one as the reference (the one that stays in place) and then matching each of the other three to that one. E.g. if #1 is the reference, you could match #2 to #1, then #3 to #1, then #4 to #1. <br>
<br>
I hope this helps,<br>
Elaine<br>
-----<br>
Elaine C. Meng, Ph.D. <br>
UCSF Chimera(X) team<br>
Department of Pharmaceutical Chemistry<br>
University of California, San Francisco<br>
<br>
> On Mar 28, 2021, at 3:17 PM, Cardiff Jiang <<a href="mailto:z8jiang@ucsd.edu" target="_blank">z8jiang@ucsd.edu</a>> wrote:<br>
> <br>
> Dear Elaine,<br>
> <br>
> Thank you for your help! The individual balls/ beads of the genome models turn out to be atoms.<br>
> <br>
> Is it possible to align four genome models, each with a different number of "atoms"? The smallest atom count is 49 and using the following command sequentially didn't work. The commands aligned model #1 to #2, then separated #1 from #2, finally aligned #1 to #3, and so on:<br>
> align #1 to #2 #1@@serial_number>=1 & #1@@serial_number<=49 #2@@serial_number>=1 & #2@@serial_number<=49<br>
> align #1 to #2 #1@@serial_number>=1 & #1@@serial_number<=49 #3@@serial_number>=1 & #3@@serial_number<=49 <br>
> <br>
> Combining the sequential commands together results in "Unequal number of atoms to pair, 49 and 0":<br>
> align #1 to #2-4 #1@@serial_number>=1 & #1@@serial_number<=49 #2@@serial_number>=1 & #2@@serial_number<= 49 #3@@serial_number>=1 & #3@@serial_number<=49 #4@@serial_number>=1 & #4@@serial_number<= 49 <br>
> <br>
> How should I align four models to each other? Thank you.<br>
> <br>
> Best regards,<br>
> Zichen "Cardiff" Jiang<br>
> <br>
> On Sat, Mar 27, 2021 at 6:35 PM Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu" target="_blank">meng@cgl.ucsf.edu</a>> wrote:<br>
> Hi Cardiff,<br>
> It is impossible to tell from a picture what the atoms are named or how the residues are numbered. You have to figure out how to specify an equal number of "atoms" from the two files, based on how the atoms are named and how the residues are numbered. Matchmaker is for biopolymers (proteins or nucleic acids made of real atoms, not your chromosome models), so you would need to use "align" after figuring out how to specify equal numbers of "atoms." <br>
> <br>
> You might need to look at your PDB files in a text editor to understand its naming/numbering. Do they really contain 40 residues? Or is each one really a single residue with 40 atoms? <br>
> <br>
> "#1:1-40" specifies residues numbered 1-40 in model #1 because the colon symbol ":" means residues. Command "select #1:1-40" will report how many "atoms" that specifies. It might be 0 atoms if the specification is wrong. If you really meant "atoms" numbered 1-40 (not residues) the atom specification should instead be "#1@@serial_number>=1 & #1@@serial_number<=40"<br>
> <br>
> How to specify atoms in the command line, where # : @ are symbols for hierarchical levels model residue atom, and @@ refers to atom attributes:<br>
> <<a href="https://urldefense.com/v3/__https://rbvi.ucsf.edu/chimerax/docs/user/commands/atomspec.html__;!!Mih3wA!RYB58EhB-zYqHVp_yrqydhbVxbphw5nBtvKXW9lzRfrV6VCXRpMW0mfZVh0Ft8b3$" rel="noreferrer" target="_blank">https://urldefense.com/v3/__https://rbvi.ucsf.edu/chimerax/docs/user/commands/atomspec.html__;!!Mih3wA!RYB58EhB-zYqHVp_yrqydhbVxbphw5nBtvKXW9lzRfrV6VCXRpMW0mfZVh0Ft8b3$</a> ><br>
> <br>
> I hope this helps,<br>
> Elaine<br>
> -----<br>
> Elaine C. Meng, Ph.D. <br>
> UCSF Chimera(X) team<br>
> Department of Pharmaceutical Chemistry<br>
> University of California, San Francisco<br>
> <br>
> > On Mar 27, 2021, at 11:45 AM, Cardiff Jiang <<a href="mailto:z8jiang@ucsd.edu" target="_blank">z8jiang@ucsd.edu</a>> wrote:<br>
> > <br>
> > Dear ChimeraX User List,<br>
> > <br>
> > How can I compare two chromosome 3D models (PDB files) that have different numbers of "atoms"?<br>
> > <br>
> > Matchmaker tool reports, "Reference and/or match model contains no nucleic or amino acid chains. Use the command-line 'align' command to superimpose small molecules/ligands." align #1 toAtoms #2 reports, "Unequal number of atoms to pair, 49 and 63." align #1 toAtoms #2 matchNumbering true reports, "Pairing dropped 49 atoms and 63 reference atoms. No atoms paired for alignment." align #1:1-40 to #2:1-40 reports, "No atoms paired for alignment."<br>
> > <br>
> > How should I compare two chromosome 3D models below? Thank you.<br>
> > <image.png><br>
> > <br>
> > Best regards,<br>
> > Zichen "Cardiff" Jiang<br>
> > _______________________________________________<br>
> > ChimeraX-users mailing list<br>
> > <a href="mailto:ChimeraX-users@cgl.ucsf.edu" target="_blank">ChimeraX-users@cgl.ucsf.edu</a><br>
> > Manage subscription:<br>
> > <a href="https://urldefense.com/v3/__https://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users__;!!Mih3wA!RYB58EhB-zYqHVp_yrqydhbVxbphw5nBtvKXW9lzRfrV6VCXRpMW0mfZVs5P__U-$" rel="noreferrer" target="_blank">https://urldefense.com/v3/__https://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users__;!!Mih3wA!RYB58EhB-zYqHVp_yrqydhbVxbphw5nBtvKXW9lzRfrV6VCXRpMW0mfZVs5P__U-$</a> <br>
> <br>
<br>
</blockquote></div>