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

add Z subcommand #14

Merged
merged 21 commits into from
Aug 22, 2024
Merged

add Z subcommand #14

merged 21 commits into from
Aug 22, 2024

Conversation

rvhonorato
Copy link
Member

@rvhonorato rvhonorato commented Aug 5, 2024

I understand the idea behind this is to, given some residue selections, 1) define a plane between them, 2) draw a grid of shape beads around the geometric center of the selection that is perpendicular to this plane and 3) orient this plane at the global z-axis

Also considering the way haddock works, I assume the shape beads needs to be centered at the origin because If everything is centered at the origin then the z-coordinates of the shape will be +/-, then we just use CNS's attr z gt syntax, right?

test01

Other considerations;

  • The dimensions of the bead grid (5x5, 10x10, etc) is an argument - wouldn't it be better if it was automatically determined based on the radius of gyration? > After discussing with @VGPReys it's best indeed if we keep these as a parameter, since all sorts of dynamic assignment might incur in underestimation
  • Should these beads be a grid or circular? > Same here, best to be a grid to avoid undersampling

@rvhonorato rvhonorato linked an issue Aug 5, 2024 that may be closed by this pull request
@rvhonorato rvhonorato marked this pull request as ready for review August 13, 2024 15:33
Comment on lines +1 to +21
# Generate `z` restraints to keep the molecule aligned to the Z-axis

This subcommand generates restraints to keep the molecule aligned to the Z-axis. This is useful when you want to keep the molecule in a specific orientation during docking.

As input you need to pass at least one selection of residues, this will be used to define a plane perpendicular to the Z-axis. The restraints will be generated to keep the molecule aligned to this plane.

## Usage

To run the `z` subcommand, you need to provide the path to the PDB file, the selection of residues, the output file which will contain the shape beads, the number of beads to generate, and the distance between the beads. For example:

```bash
./haddock-restraints z \
--residues 19,83,145,119,167 \
path/to/the/input.pdb \
path/to/the/output/shape.pdb \
20 \ # spacing between the beads in angstrom - 20A
6 \ # grid size in dimension - 6x6
> z_restraints.tbl
```

Further, follow the HADDOCK documentation **\<pending\>** to use the generated restraints in the docking protocol.
Copy link
Member Author

Choose a reason for hiding this comment

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

@amjjbonvin please check

@rvhonorato rvhonorato merged commit 34db258 into main Aug 22, 2024
2 checks passed
@rvhonorato rvhonorato deleted the 13-z-surface-and-restraints-definition branch August 22, 2024 11:25
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.

Z-surface and restraints definition
2 participants