Hi RDKitters, This isn't strictly an RDKit based issue, but I guess someone here will be able to help me. I'm trying to identify the index of atoms in a molecule which have one or more hydrogen bonded to them but are not bonded to a dummy atom ([*]).
>>> mol = Chem.MolFromSmiles('[*]CCC[*]') >>> patt = Chem.MolFromSmarts('[*H1,H2,H3,H4]!~[#0]') >>> print mol.GetSubstructMatches(patt) () I'd expect this to return (1). But it would appear that I'm in the wrong, as Pipeline Pilot give the same result. Is there a SMARTS query that will give me the matches that I want or do I have to go back to iterating over atoms? Many thanks in advance. Best, Nick Nicholas C. Firth | PhD Student | Cancer Therapeutics The Institute of Cancer Research | 15 Cotswold Road | Belmont | Sutton | Surrey | SM2 5NG T 020 8722 4033 | E nicholas.fi...@icr.ac.uk<mailto:nicholas.fi...@icr.ac.uk> | W www.icr.ac.uk<http://www.icr.ac.uk/> | Twitter @ICRnews<https://twitter.com/ICRnews> Facebook www.facebook.com/theinstituteofcancerresearch<http://www.facebook.com/theinstituteofcancerresearch> Making the discoveries that defeat cancer [cid:image001.gif@01CE053D.51D3C4E0] The Institute of Cancer Research: Royal Cancer Hospital, a charitable Company Limited by Guarantee, Registered in England under Company No. 534147 with its Registered Office at 123 Old Brompton Road, London SW7 3RP. This e-mail message is confidential and for use by the addressee only. If the message is received by anyone other than the addressee, please return the message to the sender by replying to it and then delete the message from your computer and network.
<<inline: image001.gif>>
------------------------------------------------------------------------------
_______________________________________________ Rdkit-discuss mailing list Rdkit-discuss@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/rdkit-discuss