Skip to content

Conversation

@mb2055
Copy link
Contributor

@mb2055 mb2055 commented Nov 20, 2025

In their current form the COM restraints used by the ATM implementation always apply distance restraints that pull the bound ligand to within 0Å (+- some tolerance) of the centre of mass of the protein. This restraint works fine when the binding site is within the tolerance distance of the centre of mass, but completely breaks down beyond that point, as can be seen in #470 .

This PR adds proper calculation of the starting distance between both ligands and the centre of mass of the protein (which can either be automatically calculated or passed as a list of indices by the user).

This functionality has been tested on the system mentioned in #470 . This is what the system looks like after solvation, prior to any min/eq:
solvated

and this is what it now looks like after being run through the short minimisation/equilibration pipeline shown in the ATM tutorial:

mineq_final

Closes #470.

@mb2055 mb2055 temporarily deployed to biosimspace-build November 20, 2025 10:30 — with GitHub Actions Inactive
@mb2055 mb2055 temporarily deployed to biosimspace-build November 20, 2025 10:30 — with GitHub Actions Inactive
@mb2055 mb2055 temporarily deployed to biosimspace-build November 20, 2025 10:30 — with GitHub Actions Inactive
@mb2055 mb2055 requested a review from lohedges November 20, 2025 10:31
@mb2055 mb2055 temporarily deployed to biosimspace-build November 20, 2025 12:03 — with GitHub Actions Inactive
@mb2055 mb2055 temporarily deployed to biosimspace-build November 20, 2025 12:03 — with GitHub Actions Inactive
@mb2055 mb2055 temporarily deployed to biosimspace-build November 20, 2025 12:03 — with GitHub Actions Inactive
@mb2055 mb2055 temporarily deployed to biosimspace-build November 20, 2025 12:03 — with GitHub Actions Inactive
Copy link
Contributor

@lohedges lohedges left a comment

Choose a reason for hiding this comment

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

Thanks for the quick fix, @mb2055.

@lohedges lohedges merged commit ead4092 into devel Nov 20, 2025
4 of 5 checks passed
@lohedges lohedges deleted the fix_470 branch November 20, 2025 13:45
lohedges added a commit that referenced this pull request Nov 20, 2025
lohedges added a commit that referenced this pull request Nov 20, 2025
@lohedges lohedges added the bug Something isn't working label Nov 20, 2025
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

bug Something isn't working

Projects

None yet

Development

Successfully merging this pull request may close these issues.

Bound ligand falls out of the pocket during ATMMinimization

3 participants