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

RdkitConformerFeaturizer #3378

Merged
merged 17 commits into from
May 10, 2023
Merged

RdkitConformerFeaturizer #3378

merged 17 commits into from
May 10, 2023

Conversation

tonydavis629
Copy link
Collaborator

Description

This PR implements a featurizer which uses RDkit to generate conformers and 3d coordinates

Type of change

Please check the option that is related to your PR.

  • Bug fix (non-breaking change which fixes an issue)
  • New feature (non-breaking change which adds functionality)
  • Breaking change (fix or feature that would cause existing functionality to not work as expected)
    • In this case, we recommend to discuss your modification on GitHub issues before creating the PR
  • Documentations (modification for documents)

Checklist

  • My code follows the style guidelines of this project
  • Run yapf -i <modified file> and check no errors (yapf version must be 0.32.0)
  • Run mypy -p deepchem and check no errors
  • Run flake8 <modified file> --count and check no errors
  • Run python -m doctest <modified file> and check no errors
  • I have performed a self-review of my own code
  • I have commented my code, particularly in hard-to-understand areas
  • I have made corresponding changes to the documentation
  • I have added tests that prove my fix is effective or that my feature works
  • New unit tests pass locally with my changes
  • I have checked my code and corrected any misspellings

@tonydavis629
Copy link
Collaborator Author

The main issue which I think we need to discuss is the handling of conformers in the featurization. We input one smiles, and get get out n number of conformers. How are those packaged? I think what makes the most sense is a list of GraphData objects. This requires a little bit different handling when doing batching, this is handled by flattening the batch and then forming a BatchGraphData, as seen in the doctest and unit test. The other alternative is a BatchGraphData object, but batching multiple BatchGraphData objects will require a change to GraphData, so I opted out of that.

The other question is whether we want to enable limits on the RMSD error for a conformer to be returned, because if you set self.num_conformers too high you will get unrealistic molecules. Not sure how important that is, or if we can leave that up to the user to decide.

(18, 3)
"""

def __init__(self, num_conformers: int = 1):
Copy link
Member

Choose a reason for hiding this comment

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

Can you a docstring?

Copy link
Collaborator Author

Choose a reason for hiding this comment

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

done

Copy link
Member

@rbharath rbharath left a comment

Choose a reason for hiding this comment

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

Documenting offline discussion, @tonydavis629 to proceed with fixed number of conformers, a quality rejection threshold for conformers, and to either round-robin repeat or use Nones in the case that there are insufficient high quality conformers.

Copy link
Member

@rbharath rbharath left a comment

Choose a reason for hiding this comment

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

LGTM, feel free to merge in once CI is clear

@tonydavis629 tonydavis629 merged commit 6655d17 into deepchem:master May 10, 2023
23 of 33 checks passed
@tonydavis629 tonydavis629 deleted the 3dfeat branch May 10, 2023 16: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.

None yet

2 participants