Re: [Rdkit-discuss] Exhaustive Library Enumeration

2018-01-17 Thread Christos Kannas
Hi Andy, A better option is to sanitize the products of a reaction enumeration before using them as reactants. Look at this example from RDKit "Getting Started" documentation. Note that the molecules that are produced by the chemical reaction processing code are not sanitized, as this artificial

Re: [Rdkit-discuss] Exhaustive Library Enumeration

2018-01-17 Thread Andy Jennings
Hi Christos, Many thanks for the reply. I hadn't appreciated that the presence of a single invalid reagent would bring the entire thing crashing down, rather than issuing a warning/error and moving onto other molecules in the set. Good to know, and I'll have to be less lazy in my code ;-) Best,

Re: [Rdkit-discuss] Exhaustive Library Enumeration

2018-01-17 Thread Christos Kannas
Hi Andy, The reason that your code breaks is that the second product of the third iteration ( 'NN(Cc1c1)(Cc1c1)Cc1c1') is not a valid molecule. And when calling Chem.MolFromSmiles( 'NN(Cc1c1)(Cc1c1)Cc1c1') it creates a None object. So you have to filter out the

[Rdkit-discuss] Exhaustive Library Enumeration

2018-01-17 Thread Andy Jennings
Hi RDKitters, I have a question and an observation on the topic of library enumeration. First, the question: is there a call within RDKit to trigger the exhaustive reaction of reagents? For example, if I have two reagents - a primary amine and an akyl chloride - can I tell RDKit to enumerate the