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 =
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
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
3 matches
Mail list logo