<div dir="ltr">Thank you Elaine !<div>Sorry I didn't cc the list before.</div><div><br></div><div>I think I need to be a little bit more specific. I am doing the following:</div><div><br></div><div>open 1ycr</div><div>name frozen pocket1 (#1/B:19 & sidechain) @<6 & (#1/A & sidechain)</div><div><br></div><div>and I want to list the distances between the atoms of sel1 and sel2 that are lower than 6A</div><div>I could use contacts, but it includes things that are over 6A.</div><div><br></div><div>Thank you!</div><div><br></div><div>Juliana</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El mar, 15 de jun. de 2021 a la(s) 14:59, Elaine Meng (<a href="mailto:meng@cgl.ucsf.edu">meng@cgl.ucsf.edu</a>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Of course, if you only want to use alpha-carbons instead of all atoms of all residues, you could define the atom sets that way, e.g.<br>
<br>
name frozen blah #1/A:100-200@ca<br>
contacts #1/A:55@ca restrict blah overlap -1000 log true<br>
<br>
Also, I was using named selection because you asked about it. You do not need to use a named selection. Instead you can just specify both atom sets in "contacts," e.g.<br>
<br>
contacts #1/A:55@ca restrict #1/A:100-200@ca overlap -1000 log true<br>
<br>
Elaine<br>
<br>
> On Jun 15, 2021, at 10:48 AM, Elaine Meng via ChimeraX-users <<a href="mailto:chimerax-users@cgl.ucsf.edu" target="_blank">chimerax-users@cgl.ucsf.edu</a>> wrote:<br>
> <br>
> <br>
>> On Jun 15, 2021, at 9:33 AM, Juliana Glavina <<a href="mailto:jglavina@gmail.com" target="_blank">jglavina@gmail.com</a>> wrote:<br>
>> <br>
>> Hi Elaine,<br>
>> I am trying to calculate the distance between one residue and a set of residues (within a name selection). I couldn't do this in Chimera either and I couldn't find the answer online.<br>
>> Thanks!<br>
>> Juliana<br>
> <br>
> Hi Juliana,<br>
> (We recommend using the chimerax-users address CC'd here for questions, so that the Q & A will be archived and searchable.)<br>
> <br>
> One way to do it in ChimeraX is with the "contacts" command or Contacts tool. In this case you aren't necessarily using it to find contacts, but instead to measure and list atom-atom distances.<br>
> <br>
> You say "residues" but really all the distances are measured from atom to atom. You will get all by all distances between two sets of atoms, where the first set is your first residue and the second set is all the residues in your named selection.<br>
> <br>
> In ChimeraX, you can name the current selection or some other specification with the "name" command. To name the current selection "blah" it would be command:<br>
> <br>
> name frozen blah sel<br>
> <br>
> See "name" help:<br>
> <<a href="https://rbvi.ucsf.edu/chimerax/docs/user/commands/name.html" rel="noreferrer" target="_blank">https://rbvi.ucsf.edu/chimerax/docs/user/commands/name.html</a>><br>
> <br>
> Then if you want to measure all atom-atom distances between (1) atoms in residue 55 of chain A of model 1 and (2) your named selection "blah" and show them in the Log, you could use a command something like<br>
> <br>
> contacts #1/A:55 restrict blah overlap -1000 log true<br>
> <br>
> The large negative overlap value is used to force measuring all the distances even if the atoms are 1000 Angstroms apart, instead of the usual behavior of finding only the contacts. The Log will have lines like the following, where the last column is the distance between atomic centers:<br>
> <br>
> 1728 contacts<br>
> atom1 atom2 overlap distance<br>
> /A BGC 310 C6 /A TYR 10 CE2 0.018 3.622<br>
> /A BGC 310 O6 /A HIS 152 NE2 -0.134 2.834<br>
> /A BGC 310 O4 /A TRP 183 CG -0.202 3.272<br>
> /A BGC 310 C3 /A TRP 183 CE3 -0.253 3.893<br>
> /A BGC 310 C2 /A PHE 16 CD1 -0.323 3.963<br>
> /A BGC 310 C5 /A HIS 152 NE2 -0.360 3.880<br>
> /A BGC 310 C6 /A HIS 152 NE2 -0.366 3.886<br>
> /A BGC 310 O4 /A TRP 183 CD2 -0.424 3.494<br>
> /A BGC 310 C5 /A TRP 183 CE2 -0.429 3.919<br>
> /A BGC 310 O6 /A HIS 152 CD2 -0.490 3.710<br>
> /A BGC 310 C3 /A TRP 183 CZ3 -0.495 4.135<br>
> /A BGC 310 O4 /A TRP 183 CD1 -0.499 3.719<br>
> /A BGC 310 O6 /A HIS 152 CE1 -0.556 3.776<br>
> /A BGC 310 C5 /A TRP 183 NE1 -0.596 4.116<br>
> [...]<br>
> <br>
> The lines are in order of largest VDW overlap which is similar to shortest distance but not exactly the same. If you save to text file you can then sort by column, for example, by the second residue number to group those together, or by the distance.<br>
> <br>
> The "contacts" command has a lot of options including saving results to text file, as detailed here:<br>
> <<a href="https://rbvi.ucsf.edu/chimerax/docs/user/commands/clashes.html" rel="noreferrer" target="_blank">https://rbvi.ucsf.edu/chimerax/docs/user/commands/clashes.html</a>><br>
> <br>
> You can do the same thing in Chimera, name selection and use "findclash" but the exact commands and their syntax are slightly different.<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>
> <br>
> <br>
> <br>
> <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://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users" rel="noreferrer" target="_blank">https://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users</a><br>
<br>
</blockquote></div>