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

chiral centers are not respected during generation #1

Open
missuse opened this issue Mar 17, 2022 · 1 comment
Open

chiral centers are not respected during generation #1

missuse opened this issue Mar 17, 2022 · 1 comment

Comments

@missuse
Copy link

missuse commented Mar 17, 2022

Hi iwatobipen,

and thanks for a cool tool.

I started experimenting with it to generate some glycoside conformations and noticed that chirality is not respected during this.

See attached images:

this one is as it should be:
glucose_1

this one is off:
glucose_2

File to reproduce the problem is the 3D sdf file from:

PUBCHEM_COMPOUND_CID: 88708

code I used is

confgen -i Gentiopicroside.sdf -o Gentiopicroside_chiral.sdf -t 0.1 --numconf 50 -r true

I am not an rdkit user to make suggestions.

Thanks,

Milan

@missuse missuse changed the title chiral centers with H are not respected during generation chiral centers are not respected during generation Mar 17, 2022
@missuse
Copy link
Author

missuse commented Mar 17, 2022

It seems to me that after this line:

mol = Chem.AddHs(Chem.MolFromMolFile(input), addCoords=True)

information about chirality should be added:

Chem.rdmolops.AssignStereochemistryFrom3D(mol)
Chem.rdmolops.AssignAtomChiralTagsFromStructure(mol)

It looks like this fixed it.

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

No branches or pull requests

1 participant