Hi,

I have a molecule that RDKit claims to have 9 rotatable bonds but the
chemists in my office claim to have 6. It's smiles is
CCCCNC(=O)C1CC[NH+](CC(=O)NC(N)=O)CC1
and my colleagues are concerned that the amide and urea groups are being
counted as rotatable incorrectly. Is this a bug, a design decision or a
feature?

The below code snippet shows the problem:

>>> import rdkit
>>> from rdkit import Chem
>>> from rdkit.Chem import AllChem
>>> mol = Chem.MolFromSmiles('CCCCNC(=O)C1CC[NH+](CC(=O)NC(N)=O)CC1')
>>> print AllChem.CalcNumRotatableBonds(mol)
9


Thanks,

Yours,

Toby Wright

--
InhibOx Ltd
Oxford
------------------------------------------------------------------------------
Rapidly troubleshoot problems before they affect your business. Most IT 
organizations don't have a clear picture of how application performance 
affects their revenue. With AppDynamics, you get 100% visibility into your 
Java,.NET, & PHP application. Start your 15-day FREE TRIAL of AppDynamics Pro!
http://pubads.g.doubleclick.net/gampad/clk?id=84349831&iu=/4140/ostg.clktrk
_______________________________________________
Rdkit-discuss mailing list
Rdkit-discuss@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/rdkit-discuss

Reply via email to