Hi Elaine,

Thanks a lot for this. The models that I am analyzing are separated by distances in the range of 20 to  >100 angstroms, for which the contacts command/interfaces becomes rather slow. Sounds like I might need a solution outside the Chimera universe.

Thanks for your help,
Best, Matthias

Holen Sie sich Outlook für Android

From: Elaine Meng <meng@cgl.ucsf.edu>
Sent: Tuesday, September 28, 2021 5:19:17 PM
To: Vorländer,Matthias Kopano <matthias.vorlaender@imp.ac.at>
Cc: chimerax-users@cgl.ucsf.edu <chimerax-users@cgl.ucsf.edu>
Subject: Re: [chimerax-users] Possible to measure the shortest distance between two models?
 
Hi Matthias,
The "contacts" or "clashes" command gives all-by-all distances (and VDW overlaps) between two sets of atoms.  You can write all the info to a file or to a log.  It does not automatically only give you the shortest distance, but it sorts the output so that the shortest is first.

For example (these structures are nonphysically on top of each other, but work for illustrating the commands):

open 1zik
open 1gcn
contacts #1 restrict #2 distanceOnly 5 log true

...shows in the log all the pairwise distances between the two sets of atoms that are within 5 angstroms, with the shortest distance first.
  
See measurements:
<https://rbvi.ucsf.edu/chimerax/docs/user/measurements.html>

Clashes or contacts command:
<https://rbvi.ucsf.edu/chimerax/docs/user/commands/clashes.html>

For scripting, instead of "log" you could use "saveFile" to get an output file for each structure pair.  What you get in the Log or file looks like this:

Ignore distances between atoms separated by 4 bonds or less
Detect intra-residue distances: False
Detect intra-molecule distances: True

3030 distances
       atom1                 atom2          overlap  distance
1gcn #1/A ARG 18 CZ   1zik #2/A LYS 15 O     4.767    0.233
1gcn #1/A LYS 12 O    1zik #2/A VAL 9 CA     4.568    0.432
1gcn #1/A ARG 18 CD   1zik #2/A LYS 16 CA    4.540    0.460
1gcn #1/A ALA 19 O    1zik #2/B LYS 15 CB    4.449    0.551
1gcn #1/A TRP 25 CE3  1zik #2/B GLU 22 CB    4.415    0.585
1gcn #1/A SER 8 OG    1zik #2/A LYS 3 CE     4.407    0.593
1gcn #1/A PHE 22 CE2  1zik #2/B HIS 18 N     4.392    0.608
1gcn #1/A LYS 12 NZ   1zik #2/A LEU 5 CB     4.356    0.644
1gcn #1/A LYS 12 CD   1zik #2/A LEU 5 O      4.335    0.665
1gcn #1/A GLN 20 CA   1zik #2/B LYS 15 CD    4.309    0.691
[... etc. ...]

I hope this helps,
Elaine
-----
Elaine C. Meng, Ph.D.                       
UCSF Chimera(X) team
Department of Pharmaceutical Chemistry
University of California, San Francisco

On Sep 28, 2021, at 2:15 AM, Vorländer,Matthias Kopano via ChimeraX-users <chimerax-users@cgl.ucsf.edu> wrote:

Dear all, 
I am interested in measuring the shortest distance between two models (PDBs) without a priori knowledge of which atoms will be the closest ones. I am planning to do this for a large number of model pairs, therefore it would be great if this would be scriptable. Is this possible using ChimeraX?
Thanks a lot in advance, 
Best,
Matthias