Hi Ruben,

I ran into this issue a few months ago. Take a look at this thread :
https://www.mail-archive.com/rdkit-discuss@lists.sourceforge.net/msg07311.html

And also Connor Coley's fork of rdkit:


https://github.com/connorcoley/Rdkit
<https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_connorcoley_Rdkit&d=DwMFaQ&c=WO-RGvefibhHBZq3fL85hQ&r=UPzYrSHLXjnX3tYn90C8Ljjzb-yfrb1UtMOxOFh-tKk&m=84fKplcxw5Qax5fYEDsUlzLtw3JZn_wawz-VmLyJl7c&s=hRPMSu46pUw2GC3cj95hZINhb3t2HcQm3b7wFHqQ6qg&e=>

Example of use:

https://github.com/connorcoley/ochem_predict_nn/blob/master/data/generate_candidates_edits_fullgrants.py
<https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_connorcoley_ochem-5Fpredict-5Fnn_blob_master_data_generate-5Fcandidates-5Fedits-5Ffullgrants.py&d=DwMFaQ&c=WO-RGvefibhHBZq3fL85hQ&r=UPzYrSHLXjnX3tYn90C8Ljjzb-yfrb1UtMOxOFh-tKk&m=84fKplcxw5Qax5fYEDsUlzLtw3JZn_wawz-VmLyJl7c&s=XXN9Jsptw1Z3CQbZSTmp_QoIv04C1lAKujRMc7iYS8Q&e=>


Cheers,
Jennifer


On Thu, Dec 21, 2017 at 10:02 AM Ruben Van de Vijver <
ruben.vandevij...@ugent.be> wrote:

> Dear all,
>
>
>
> I have the following lines to create the products of a hydrogen
> abstraction reaction by a radical atom. The length of the product tuple is
> 6 as expected from the 6 hydrogen atoms of ethane that can be abstracted. A
> first question that I have is, is there a way to know for each set of
> products which hydrogen got abstracted from which carbon atom? I would need
> some sort of mapping of the atom indices of the reactants to the indices
> given by the reaction smarts.
>
> Secondly, if there is a reaction with two reactants, and both reactants
> match the two reactant templates, is it possible to know which reactant
> template was used on which reactant?
>
>
>
> rxn =
> AllChem.ReactionFromSmarts("[C;v4:0][H:1].[C;v3:2]>>[C;v3:0].[H:1][C;v4:2]")
>
> react1 = Chem.AddHs(Chem.MolFromSmiles('CC'))
>
> react2 = Chem.AddHs(Chem.MolFromSmiles('[CH3]'))
>
> prods = rxn.RunReactants((react1,react2))
>
> print len(prods)
>
>
>
> Thanks,
>
> Ruben
>
>
>
>
>
> [image:
> https://styleguide.ugent.be/files/uploads/logo_UGent_EN_RGB_2400_kleur_witbg.png]
>
> *Ruben Van de Vijver*
>
> *Postdoctoral Researcher*
>
> *Laboratory for Chemical Technology (LCT) *
>
> Department of Materials, Textiles and Chemical Engineering
>
> Technologiepark 914, 9052 Ghent, Belgium
>
> T: +32 9 264 56 75 <+32%209%20264%2056%2075>
>
> http://www.lct.ugent.be/
> <https://urldefense.proofpoint.com/v2/url?u=http-3A__www.lct.ugent.be_&d=DwMFAg&c=WO-RGvefibhHBZq3fL85hQ&r=UPzYrSHLXjnX3tYn90C8Ljjzb-yfrb1UtMOxOFh-tKk&m=48dydnCN1kbtEacvTIbqjovuTdygOqW5nQcY737lpJY&s=xJNXx2tQPNA3yX85GbgCnML9xyfGlMbgUjEa58IOX2o&e=>
> http://helpdesk.UGent.be/e-maildisclaimer.php
> <https://urldefense.proofpoint.com/v2/url?u=http-3A__helpdesk.ugent.be_e-2Dmaildisclaimer.php&d=DwMFAg&c=WO-RGvefibhHBZq3fL85hQ&r=UPzYrSHLXjnX3tYn90C8Ljjzb-yfrb1UtMOxOFh-tKk&m=48dydnCN1kbtEacvTIbqjovuTdygOqW5nQcY737lpJY&s=0hES28XZd8FZIDeS0u6ciknsFOYlTlqFHzPH8M_w70E&e=>
>
>
>
>
>
> ------------------------------------------------------------------------------
> Check out the vibrant tech community on one of the world's most
> engaging tech sites, Slashdot.org! http://sdm.link/slashdot
> _______________________________________________
> Rdkit-discuss mailing list
> Rdkit-discuss@lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/rdkit-discuss
>
------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
Rdkit-discuss mailing list
Rdkit-discuss@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/rdkit-discuss

Reply via email to