<div dir="ltr">Dear Elaine,<div><br></div><div>Thank you for your help! The individual balls/ beads of the genome models turn out to be atoms.</div><div><br></div><div>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:</div><div><font face="monospace" style="background-color:rgb(238,238,238)" color="#cc0000">align #1 to #2 #1@@serial_number>=1 & #1@@serial_number<=49 #2@@serial_number>=1 & #2@@serial_number<=49</font></div><div><font face="monospace" color="#cc0000" style="background-color:rgb(238,238,238)">align #1 to #2 #1@@serial_number>=1 & #1@@serial_number<=49 #3@@serial_number>=1 & #3@@serial_number<=49 </font><font face="monospace" style="background-color:rgb(238,238,238)" color="#cc0000"><br></font></div><br>Combining the sequential commands together results in "Unequal number of atoms to pair, 49 and 0":<div><font color="#cc0000"><font face="monospace" style="background-color:rgb(238,238,238)">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
</font><br></font></div><div><div><br></div><div>How should I align four models to each other? Thank you.</div><div><br></div><div>Best regards,</div><div>Zichen "Cardiff" Jiang<br><blockquote style="margin:0 0 0 40px;border:none;padding:0px"><div><br></div></blockquote><div><div><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Mar 27, 2021 at 6:35 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>
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>
</blockquote></div></div></div></div></div></div>