> I'm working on an OBFormat for parsing SMILES using my new Smiley parser. The 
> initial tests look good but there is a specific problem regarding aromatic 
> nitrogens. When I read a molecule using the new format and write out the 
> caonical smiles, all aromatic nitrogens are written out as "[nH]". Does 
> anyone know how the original SMILES parser handles this or how I should fix 
> this issue?

I assume you mean that, say, pyridine "c1ccccn1" becomes [nH], correct?

I believe the existing SMILES code handles this so:
>     if (arom)
>       {
>         atom->SetAromatic();
>         atom->SetSpinMultiplicity(2); // CM 18 Sept 2003
>       }
>     else
>       atom->ForceImplH();//ensures atom is never hydrogen deficient
> 
> 
>     // Untrue, but necessary to avoid perception being called in 
> OBAtom::IsAromatic()
>     // on incomplete molecule. Undone at end of function.
>     mol.SetAromaticPerceived();


This is done to handle the "lowercase SMILES might be radical" extension.

Hope that helps,
-Geoff
------------------------------------------------------------------------------
Monitor your physical, virtual and cloud infrastructure from a single
web console. Get in-depth insight into apps, servers, databases, vmware,
SAP, cloud infrastructure, etc. Download 30-day Free Trial.
Pricing starts from $795 for 25 servers or applications!
http://p.sf.net/sfu/zoho_dev2dev_nov
_______________________________________________
OpenBabel-Devel mailing list
OpenBabel-Devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/openbabel-devel

Reply via email to