Find Cavities detects cavities or pockets in an atomic structure using KVFinder, a method based on surface shape, as described in:
parKVFinder: A thread-level parallel approach in biomolecular cavity detection. Guerra JVDS, Ribeiro-Filho HV, Bortot LO, Honorato RV, Pereira JGC, Lopes-de-Oliveira PS. SoftwareX. 2020 Jul-Dec;100606.
pyKVFinder: an efficient and integrable Python package for biomolecular cavity detection and characterization in data science. Guerra JVDS, Ribeiro-Filho HV, Jara GE, Bortot LO, Pereira JGC, Lopes-de-Oliveira PS. BMC Bioinformatics. 2021 Dec 20;22(1):607.
Find Cavities installs and runs pyKVFinder locally. It is also implemented as the kvfinder command. The authors of KVFinder (at the Brazilian Biosciences National Lab) also provide a web service, but it is not used in the ChimeraX implementation. See also: measurements, Protein-Ligand Binding Sites tutorial, MOLE Channel Coloring tutorial, recipe for reading CASTp pockets
Find Cavities can be started from the Binding Analysis or Structure Analysis section of the Tools menu.
One or more atomic models should be chosen from the list that appears. Clicking OK runs the calculation, and tabulates the results for each atomic model in a separate cavity list window. These results are also reported in the Log. PDB HET residues (solvent, ligands, and ions) are ignored, and default parameters are used. For details on the method and/or to use nondefault values, see the kvfinder command.
In the graphics window, the cavities are shown as clouds of points on a grid. These collections of dots are opened as (pseudo)atomic models in ChimeraX, each cavity as a separate submodel assigned a unique color.
Running Find Cavities (or the kvfinder command with showTool true) on an atomic model opens a separate window listing the cavities that were found in the model, with sortable columns:
If a cavity list is closed but the corresponding cavity models still exist, the list can be reshown by choosing it from the bottom section of the Tools menu.
One or more cavities can be chosen from the list by clicking and dragging with the left mouse button; Ctrl-click (or command-click if using a Mac) toggles whether a row is chosen. Chosen rows are highlighted in the dialog, and the following options (if turned on) are applied to the corresponding cavity models:
"Nearby" atoms/residues are within [distance] angstroms of cavity points – distance criterion for the nearby checkbox options above
The Find Cavities tool includes the option to show surfaces around the clouds of grid points. It might be interesting to color these surfaces by properties from the surrounding protein. However, since only the cavity points (and not the protein atoms) were used to calculate the surfaces, the coloring command must specify using the protein atoms to calculate the properties for coloring these surfaces.
For example, if the option to show cavity surfaces is turned on, and the protein structure is model #1 and the computed cavity models are #1.2.1, #1.2.2, (etc.), the following commands could be used to color the surfaces by Coulombic electrostatic potential (ESP) and and molecular lipophilic potential (MLP), respectively:
coulombic (#1 & protein) surfaces #1.2 offset -1
mlp (#1 & protein) surfaces #1.2
A negative offset is used to show the ESP inside the cavity surfaces, where binding might occur. For the usual case of coloring a protein surface by its own potential, a positive offset (default 1.4 Å) is used to show the potential outside the surface, where interactions with other molecules would occur.
These types of coloring can be combined with surface transparency, for example:
transparency #1.2 30