Hi
I'm interested in the number of pi electrons. I have a molecule which has
smiles C=Cc1ccccc1. And
mol = Chem.MolFromSmiles('C=Cc1ccccc1')
for atom in mol.GetAtoms():
npibond = Chem.AtomPairs.Utils.NumPiElectrons(atom)
print npibond
gives
1
1
1
1
1
1
1
1
as expected.
However the molecule comes from a reaction and is coded as
[CH2:3]=[CH:1][c:2]1[cH:5][cH:7][cH:9][cH:8][cH:6]1
Using this in the above gives
3
2
1
1
1
1
1
1
Surely the atom-mapping should not make a difference?
Thanks
Eleanor
--
Dr Eleanor Gardiner
Information School
University of Sheffield
------------------------------------------------------------------------------
Comprehensive Server Monitoring with Site24x7.
Monitor 10 servers for $9/Month.
Get alerted through email, SMS, voice calls or mobile push notifications.
Take corrective actions from your mobile device.
http://p.sf.net/sfu/Zoho
_______________________________________________
Rdkit-discuss mailing list
Rdkit-discuss@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/rdkit-discuss
