The plethora of conformations obtained from biomolecular simulations make it difficult to analyze conformational changes. Geometrically Stable Sustructures (GeoStaS) is a program with graphical interface which, based on molecular conformations, divides a molecule into "dynamic domains", i.e., parts that stay relatively rigid in all conformations. This division simplifies the description of differences between the conformations.
The algorithm implemented in GeoStaS is based on searching for geometrical similarities between atomic motions. The pairwise atomic movement similarity matrix (AMSM; left panel in the figure) is constructed and then clustered with the use of either specifically adapted nearest-neighbor clustering algorithm (which suggests an optimal solution) or a hierarchical merging algorithm (which requires inputting the number of desired domains).
GeoStaS is implemented in Java with user-friendly the graphical interface. It is fast (although it depends on the size of the trajectories the user uploads) and robust - it can analyze conformations of both proteins and nucleic acids. There are two types of output from GeoStaS: (a) the molecule divided into dynamic domains, in the PDB format (visualized in the middle and right panels in the figure), and (b) AMSM, in text format.
NEW! The article describing details of the algorithm and tests is published in Journal of Chemical Theory and Computation! (DOI: 10.1021/ct300206j). Please, cite this article if using GeoStaS.
Downloads and source code available at GeoStaS page in BitBucket
Questions? Comments? e-mail me!