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 
<anthony.n...@ndcn.ox.ac.uk<mailto:anthony.n...@ndcn.ox.ac.uk>> 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
Rdkit-discuss@lists.sourceforge.net<mailto:Rdkit-discuss@lists.sourceforge.net>
https://eur02.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.sourceforge.net%2Flists%2Flistinfo%2Frdkit-discuss&amp;data=04%7C01%7Cjhjensen%40chem.ku.dk%7C57815f9f2e754e262a5408d895fa14d0%7Ca3927f91cda14696af898c9f1ceffa91%7C0%7C1%7C637424248612593024%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C2000&amp;sdata=vCo931FmtaS2hJ%2FCHVrfzFaJO32h4kCaA252sUHmxrs%3D&amp;reserved=0

_______________________________________________
Rdkit-discuss mailing list
Rdkit-discuss@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/rdkit-discuss

Reply via email to