One option is to use AllChem.MMFFOptimizeMolecule(mol3D,
ignoreInterfragInteractions=False), but I am not sure MMFF can handle Au.
Another option is to define the Au-S bond ('C(C1C(C(C(C(O1)S[Au])O)O)O)O’) and
use AllChem.UFFOptimizeMolecule(mol3D)
Best regards, Jan
On 1 Dec 2020, at 13.36, Anthony Nash
<[email protected]<mailto:[email protected]>> wrote:
Dear all,
I'm new to RDKit and cheminformatics in general. I'm using the latest RDKit
libraries. Any suggestions you can offer would be kindly received.
Using the canonical SMILES C(C1C(C(C(C(O1)[S-])O)O)O)O.[Au+] of
Aurothioglucose (Pubchem CID: 454937) I've generated a 3D structure using the
python RDKit code:
mol = Chem.MolFromSmiles(self.canonicalSMILES)
mol3D = Chem.AddHs(mol)
AllChem.EmbedMolecule(mol3D)
AllChem.MMFFOptimizeMolecule(mol3D)
Unfortunately, the mol3D representation has Au+ right in the middle of and in
the plane of the benzene rings, too far from the negatively charged sulfur. I'm
new at generating structures from SMILES. Are there any steps I'm missing that
could help correct the placement of Au+?
Thanks
Anthony
_______________________________________________
Rdkit-discuss mailing list
[email protected]<mailto:[email protected]>
https://eur02.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.sourceforge.net%2Flists%2Flistinfo%2Frdkit-discuss&data=04%7C01%7Cjhjensen%40chem.ku.dk%7C57815f9f2e754e262a5408d895fa14d0%7Ca3927f91cda14696af898c9f1ceffa91%7C0%7C1%7C637424248612593024%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000&sdata=vCo931FmtaS2hJ%2FCHVrfzFaJO32h4kCaA252sUHmxrs%3D&reserved=0
_______________________________________________
Rdkit-discuss mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/rdkit-discuss
