Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

restraint subcommand should likely specify "CA" for a second selection #29

Closed
edikedik opened this issue Nov 8, 2024 · 3 comments · Fixed by #30
Closed

restraint subcommand should likely specify "CA" for a second selection #29

edikedik opened this issue Nov 8, 2024 · 3 comments · Fixed by #30
Assignees
Labels
bug Something isn't working

Comments

@edikedik
Copy link

edikedik commented Nov 8, 2024

Using the tool to generate restraints for a two-chain protein (Fab), the output looks as follows:

assign ( resid 4 and segid B and name CA ) ( resid 81 and segid B ) 14.5 0.0 0.0

...

The distance is calculated between the CA atoms of a residue pair, but the atom is only specified for the first residue. Then, using these restraints during the [flexref] stage, chains that were expected to keep their mutual orientations drifted away from each other instead. Including "CA" into the second selection solved this problem, hence this issue.

@amjjbonvin amjjbonvin added the bug Something isn't working label Nov 8, 2024
@rvhonorato
Copy link
Member

Thanks for reporting this @edikedik, could you please tell us which type of restraint are you trying to generate - which subcomand are you using?

@edikedik
Copy link
Author

@rvhonorato the subcommand is restraint as I specified in the title :)
I'm running it on a two-chain protein (antibody) preprocessed to impose a uniform chain ID (B) and continuous residue numbering.

@rvhonorato rvhonorato changed the title restraint should likely specify "CA" for a second selection restraint subcommand should likely specify "CA" for a second selection Nov 11, 2024
@rvhonorato rvhonorato linked a pull request Nov 11, 2024 that will close this issue
@rvhonorato
Copy link
Member

Fixed in v0.6.2, thanks again @edikedik. Please let me know if you run into any other issues.

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 a pull request may close this issue.

3 participants