aromatic compounds

Previous Topic Next Topic
classic Classic list List threaded Threaded
1 message Options
Reply | Threaded
Open this post in threaded view

aromatic compounds

I am facing a problem when manipulating an aromatic molecule: [nH]1nnc2c1c(=O)c1[nH]nnc1c2=O
Cycles are not in kekulize form.
When I create a molecule with pybel:
>>> mol = pybel.readstring("smi","O=c1c2nn[nH]c2c(=O)c2[nH]nnc12").next()

 pybel create a connectivity error with 2 pentavalent carbon atom (carbonyls) in the cycle:
>>> mol.write().split()[0]

how to get around the problem?

thanks a lot