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 [email protected]<mailto:[email protected]> |
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:[email protected]]
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 [email protected] https://lists.sourceforge.net/lists/listinfo/rdkit-discuss
