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

david.lagorce
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]
O=C1=C2N=NNC2C(=O)=C2NN=NC12

how to get around the problem?

thanks a lot

David