Skip to content

Conversation

@sbliven
Copy link
Member

@sbliven sbliven commented Aug 30, 2024

This fixes rcsb/symmetry#118.

It switches the determination of the symmetry group and rotation axes to use the CESymm alignment directly. Before it was recomputed using QuatSymmDetector, which failed with short repeats.

A new CE-Symm version will be released after the next biojava release.

Most are in docstrings. One misspelled public method has been deprecated.
[rcsb/symmetry#118](rcsb/symmetry#118)

The problem was actually downstream from CE-Symm in the QuatSymmDetector
code used for assigning a symmetry group. The code was re-aligning
subunits using CE. This fails for short alignments of <15 residues
(two AFP with default settings). Thus no atoms get aligned and an
exception was thrown whenever an output format required information
about the symmetry group.

This fix re-uses the original CE-Symm alignment for the symmetry group
detection. This might change CE-Symm and QuatSymm results for proteins
with small subunits or borderline scores. 

Changes:
- New `SubunitClustererMethod.MANUAL` enum value
- Improve `SymmetryTools.getQuaternarySymmetry`
Better refinement now finds slightly larger helical rotation,
moving it from R to H.
Copy link
Contributor

@josemduarte josemduarte left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

LGTM thanks!

@sbliven sbliven merged commit b50e990 into biojava:master Sep 4, 2024
@heuermh
Copy link
Member

heuermh commented Sep 4, 2024

Thank you, @sbliven!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

Random error detecting symmetry

3 participants