Re: [Rdkit-discuss] UpdatePropertyCache() after RunReactants

2017-01-18 Thread Greg Landrum
Coming back to an old one that was mostly answered already: On Thu, Jan 12, 2017 at 3:11 AM, Curt Fischer wrote: > > I recently wanted to use RDKit to model the famous copper-catalyzed > cycloaddition of alkynes and azides. > > I eventually got things working, kind of, but had two questions. Fi

Re: [Rdkit-discuss] UpdatePropertyCache() after RunReactants

2017-01-12 Thread Curt Fischer
Oof, thanks for the pointer MIchal. I'm sorry I didn't read the docs carefully enough! ~CF On Thu, Jan 12, 2017 at 9:32 AM, Michal Krompiec wrote: > You need to sanitize the products, just run Chem.SanitizeMol on each > molecule. See http://www.rdkit.org/docs/GettingStartedInPython.html# > chem

Re: [Rdkit-discuss] UpdatePropertyCache() after RunReactants

2017-01-12 Thread Michal Krompiec
You need to sanitize the products, just run Chem.SanitizeMol on each molecule. See http://www.rdkit.org/docs/GettingStartedInPython.html#chemical-reactions : "the molecules that are produced by the chemical reaction processing code are not sanitized". Best, Michal On 12 January 2017 at 17:22, Cur

Re: [Rdkit-discuss] UpdatePropertyCache() after RunReactants

2017-01-12 Thread Curt Fischer
What makes you think the molecules are nonsensical? They look OK to me. Converting to SMILES before doing any UpdatePropertyCache() stuff *products_tuples = copper_click.RunReactants((diyne, azide))products = list(chain(*products_tuples))print [Chem.MolToSmiles(prod) for prod in products]*

Re: [Rdkit-discuss] UpdatePropertyCache() after RunReactants

2017-01-12 Thread Brian Kelley
The outputs of reaction are a bit confusing. Reactions can have multiple product templates so the output of RunReactants is a list of list of molecules. For products in result: For molecule in products: Molecule.UpdatePropertyCache() However, it looks like your reaction is generating non