rdchiral icon indicating copy to clipboard operation
rdchiral copied to clipboard

Explicit hydrogens

Open kjelljorner opened this issue 3 years ago • 1 comments

We've encountered a case where we want to run a reaction with explicit hydrogens. From what I can see, the default behavior of rdchiral is to strip them at the stage of creating the reactants. Would it be possible to do just change that line to allow explicit hydrogens, or would that break things further along in the code?

https://github.com/connorcoley/rdchiral/blob/01cca8c7f5b0946187f47928738730040de06b16/rdchiral/initialization.py#L180

A possible change would be:

params = Chem.SmilesParserParams()
params.removeHs = False
reactants = Chem.MolFromSmiles(reactant_smiles, params)

kjelljorner avatar Dec 07 '21 22:12 kjelljorner

I made a fork with this change that we will test out: https://github.com/kjelljorner/rdchiral. Can open a pull request if needed.

kjelljorner avatar Dec 07 '21 22:12 kjelljorner