Hi JP: On Tue, Apr 12, 2011 at 5:35 AM, JP <jeanpaul.ebe...@inhibox.com> wrote: > > Copying the example verbatim from the documentation (section 3.5 working > with 3d molecules) doesn't work: > >>>> from rdkit import Chem >>>> from rdkit.Chem import AllChem >>>> m = Chem.MolFromSmiles('C1CCC1OC') >>>> m2 = Chem.AddHs(m) >>>> m > <rdkit.Chem.rdchem.Mol object at 0x1f02bb0> >>>> m2 > <rdkit.Chem.rdchem.Mol object at 0x1f021a0> >>>> AllChem.UFFOptimizeMolecule(m2) > Traceback (most recent call last): > File "<stdin>", line 1, in <module> > ValueError: Bad Conformer Id > > Any ideas why?
You forgot one step, before optimizing the molecule you need to add 3D coordinates: >>> AllChem.EmbedMolecule(m2) 0 -greg ------------------------------------------------------------------------------ Forrester Wave Report - Recovery time is now measured in hours and minutes not days. Key insights are discussed in the 2010 Forrester Wave Report as part of an in-depth evaluation of disaster recovery service providers. Forrester found the best-in-class provider in terms of services and vision. Read this report now! http://p.sf.net/sfu/ibm-webcastpromo _______________________________________________ Rdkit-discuss mailing list Rdkit-discuss@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/rdkit-discuss