Re: [Rdkit-discuss] mass replacement of External R-groups with many substituents

2017-03-16 Thread Greg Landrum
Combining Steve's and Chris' answer gets to how I would do it with reactions: In [17]: core = Chem.MolFromSmiles('*c1c(C)1(O)') In [18]: chain = Chem.MolFromSmiles('CN*') In [19]: rxn = AllChem.ReactionFromSmarts('[c:1][#0].[#0][*:2]>>[c:1]-[*:2]') In [20]: ps =

Re: [Rdkit-discuss] mass replacement of External R-groups with many substituents

2017-03-16 Thread Stephen Roughley
You can match a dummy atom (*) with the SMARTS [#0] Steve On 16 Mar 2017 16:43, "Chris Earnshaw" wrote: > Hi Brian > > I'm by no means an expert in RDKit with Python, but until someone else > comes along, here are a few thoughts. > > Your reaction SMARTS specifically

Re: [Rdkit-discuss] mass replacement of External R-groups with many substituents

2017-03-16 Thread Chris Earnshaw
Hi Brian I'm by no means an expert in RDKit with Python, but until someone else comes along, here are a few thoughts. Your reaction SMARTS specifically defines aromatic carbons joined by single bonds which won't match an incoming benzene ring, and it's a bit redundant to specify that aromatic