[Rdkit-discuss] sampling of ring conformation for docking

2024-05-13 Thread Pavel Polishchuk
regards, Pavel ___ Rdkit-discuss mailing list Rdkit-discuss@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/rdkit-discuss

[Rdkit-discuss] 6th Advanced in silico Drug Design workshop in Olomouc

2022-12-09 Thread Pavel Polishchuk
solve real chemoinformatic tasks and win prizes.   https://www.kfc.upol.cz/6add   Please feel free to share this information with those who can be interested in participation in such event. Thank you! Kind regards, Pavel ___ Rdkit-discuss mailing

Re: [Rdkit-discuss] Comparison of the docking conformations with the X-ray pose

2022-03-30 Thread Pavel Polishchuk
suitable for your needs. Kind regards, Pavel. On 30/03/2022 18:16, Enrico Martinez wrote: Dear RDKIT users! I am dealing with the analysis of the results of the docking poses calculated via VINA and saved into the multi-model pdb. I need to find a possibility to compare each docking pose

Re: [Rdkit-discuss] 2dpharmacophores

2021-08-12 Thread Pavel Polishchuk
Hi Chris, Error reported a bin (0,2,0) which does not satisfy the triangle rule. I cannot check it, but I suppose that the error may disappear if you set trianglePruneBins=False in SigFactory, but this would be a workaround not a proper solution I think. Pavel. On 12/08/2021 10:25, Chris

Re: [Rdkit-discuss] rejoining pairs of fragments after fragmenting a molecule

2021-04-01 Thread Pavel Polishchuk
Hi Ling,   this can be a workaround if RDKit does not have a built-in function to extract a submolecule by atom ids. You may assign atom property labels to these atoms and then looping over atoms in EditableMol remove those ones which do not have this property assigned. Kind regards, Pavel

Re: [Rdkit-discuss] ConstrainedEmbed issue

2020-07-07 Thread Pavel Polishchuk
ils embed template mol in rdkit, no AddHs 0.1012757033705761 Pavel. On 07/07/2020 21:41, Sunhwan Jo wrote: Makes sense :) On Jul 7, 2020, at 12:35 PM, Sereina Riniker mailto:sereina.rini...@gmail.com>> wrote: Dear Pavel and Sunhwan, Please note that hydrogens should always

[Rdkit-discuss] ConstrainedEmbed issue

2020-07-07 Thread Pavel Polishchuk
to solve this properly? Kind regards, Pavel. 1.mol Description: MOL mdl chemical test ___ Rdkit-discuss mailing list Rdkit-discuss@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/rdkit-discuss

[Rdkit-discuss] conformations for cycles with multiple stereocenters

2020-02-21 Thread Pavel Polishchuk
, Chem.rdDistGeom.ETKDGv2())   mol.GetConformers() Kind regards, Pavel. ___ Rdkit-discuss mailing list Rdkit-discuss@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/rdkit-discuss

Re: [Rdkit-discuss] 5th in silico drug design workshop in Olomouc

2019-11-21 Thread Pavel Polishchuk
Sorry for the mistake, the dates are 3-7 of *February* 2020. Pavel. On 20/11/2019 09:39, Pavel Polishchuk wrote: Dear colleagues,   we would like to invite you on the 5th Drug Design workshop which will be held 3-7 January 2020 in Olomouc (Czech Republic). It is focused on practical

[Rdkit-discuss] 5th in silico drug design workshop in Olomouc

2019-11-20 Thread Pavel Polishchuk
, Pavel. -- Dr. Pavel Polishchuk senior researcher Institute of Molecular and Translational Medicine Faculty of Medicine and Dentistry Palacky University Hněvotínská 1333/5 779 00 Olomouc Czech Republic +420 585632298 ___ Rdkit-discuss mailing list Rdkit

Re: [Rdkit-discuss] add hydrogen after its removal

2019-03-22 Thread Pavel
Thank you, Paolo! That was not obvious for me how RDKit manage hydrogens in this case. Pavel. On 22/03/2019 12:15, Paolo Tosco wrote: Hi Pavel, After you have first called AddHs(), all hydrogens in your molecule are now in the molecule graph as real atoms, and there are no more implicit

[Rdkit-discuss] add hydrogen after its removal

2019-03-22 Thread Pavel
]([H])[c]([H])[c]1[H] remove H [H][O][c]1[c][c]([H])[c]([H])[c]([H])[c]1[H] add Hs [H][O][c]1[c][c]([H])[c]([H])[c]([H])[c]1[H] update property cache [H][O][c]1[c][c]([H])[c]([H])[c]([H])[c]1[H] Pavel. ___ Rdkit-discuss

Re: [Rdkit-discuss] 3D pharmacophore-based screening

2019-02-20 Thread Pavel
Hi Pawan,   I can suggest to use our developed module for 3D pharmacophore representation: https://github.com/DrrDom/pmapper   I added an example to answer your particular question: https://github.com/DrrDom/pmapper/blob/master/examples/screen_example.ipynb   Hope this will help, Pavel

[Rdkit-discuss] 4th advanced in silico drug design workshop in Olomouc

2018-11-06 Thread Pavel
tasks and win prizes.   http://fch.upol.cz/en/research/conferences-workshops/4add/   Please feel free to share this information to those who can be interested in participation in such event. Thank you. Kind regards, Pavel. -- Dr. Pavel Polishchuk senior researcher Institute of Molecular

Re: [Rdkit-discuss] Fingerprint collision and machine learning

2018-10-10 Thread Pavel
Hi Michal,   I think if you can provide several examples of structures having identical bitstrings this will help a lot to better understand the issue. Pavel. On 10/10/18 14:15, Michal Krompiec wrote: Hi Thomas, Radius 2, 2048 bits, 5200 data points. On Wed, 10 Oct 2018 at 13:13, Thomas

[Rdkit-discuss] similarity maps for regression models

2018-01-24 Thread Pavel
Hello,   I'm curious, would it be possible to use regression models for similarity maps drawing or only classification ones?   If someone has a piece of code, I would appreciate this. Because I tried to modify the example given in the cookbook, but without success. Pavel

[Rdkit-discuss] access C++ code from Python

2017-12-21 Thread Pavel
= Chem.MolFromSmiles('F[C@@H](Cl)Br') m = Chem.AddHs(m) AllChem.EmbedMolecule(m) AllChem.MMFFOptimizeMolecule(m, confId=0) Descriptors3D.PMI1(m, 0) Pavel. -- Check out the vibrant tech community on one of the world's most engaging tech

Re: [Rdkit-discuss] canonical SMILES of a fragment

2017-08-02 Thread Pavel Polishchuk
, 1, Chem.rdchem.BondType.SINGLE) m.AddBond(1, 2, Chem.rdchem.BondType.SINGLE) m.AddBond(1, 3, Chem.rdchem.BondType.SINGLE) s = Chem.MolToSmiles(m, isomericSmiles=True) re.sub('\[[0-9]+([a-z]+)H?[0-9]?\]', '\\1', s) # remove isotope in output SMILES OUTPUT: 'CC(c)[*:1]' Pavel. On 08/02

[Rdkit-discuss] canonical SMILES of a fragment

2017-08-01 Thread Pavel Polishchuk
as aromatic m2.GetAtomWithIdx(3).SetAtomMapNum(1) m2.AddBond(0, 1, Chem.rdchem.BondType.SINGLE) m2.AddBond(1, 2, Chem.rdchem.BondType.SINGLE) m2.AddBond(1, 3, Chem.rdchem.BondType.SINGLE) Chem.MolToSmiles(m2) OUTPUT: 'CC(c)[*:1]' Pavel

Re: [Rdkit-discuss] Using RDKit in PyCharm and Anaconda on Windows

2017-06-01 Thread Pavel Polishchuk
I had some issues to run rdkit from Python console in PyCharm (4.5.5) on Linux. After recent installation of PyCharm 2017.1.3 it started to work. Maybe updating PyCharm will help on Win as well. Pavel. On 05/30/2017 10:10 PM, West, Richard wrote: We're having trouble getting RDKit to work

[Rdkit-discuss] RemoveHs and [H][*:1] mol

2017-05-30 Thread Pavel Polishchuk
1ccc(C[*:1])cc1 [H][*:1] [*:1] If this is a bug I can create an issue on github Pavel. -- Check out the vibrant tech community on one of the world's most engaging tech sites, Slashdot.org! http://sdm.lin

Re: [Rdkit-discuss] Redundant hydrogen atoms

2017-05-28 Thread Pavel Polishchuk
print(Chem.MolToSmiles(ps[0][0])) print(Chem.MolToSmiles(ps[0][1])) for a in ps[0][0].GetAtoms(): if a.GetAtomicNum() == 16 and a.GetFormalCharge() == 1: a.SetNoImplicit(True) print(Chem.MolToSmiles(ps[0][0])) print(Chem.MolToSmiles(ps[0][1])) Pavel. On 05/28/2017 08:26 AM, Nitzan Tzanani

Re: [Rdkit-discuss] canonical smiles for fragments with map numbers

2017-05-27 Thread Pavel Polishchuk
smiles it would be easy to relabel attachment points back. Thank you again! Pavel. On 05/27/2017 03:03 PM, Brian Kelley wrote: Pavel, this isn't exactly trivial so I went ahead and made an example. The basics are that atomMaps are canonicalized, i.e. their value is used in the generation

[Rdkit-discuss] custom fingerprints in PostgreSQL

2017-04-26 Thread Pavel Polishchuk
= Generate.Gen2DFingerprint(mol, sigFactory) Pavel. -- Check out the vibrant tech community on one of the world's most engaging tech sites, Slashdot.org! http://sdm.link/slashdot ___ Rdkit-discuss

[Rdkit-discuss] Pharm2D.SigFactory and psql issues

2017-04-24 Thread Pavel Polishchuk
1, 1, 2), (1, 2, 0), (1, 2, 1), (1, 2, 2), (2, 0, 1), (2, 0, 2), (2, 1, 0), (2, 1, 1), (2, 1, 2), (2, 2, 0), (2, 2, 1), (2, 2, 2)], 0] Kind regards, Pavel. -- Check out the vibrant tech community on one of the world's

Re: [Rdkit-discuss] DeleteSubstructs vs ReplaceSubstructs

2017-03-31 Thread Pavel Polishchuk
Hi Maxim, if you change your query to SMARTS it would be possible to delete what you want m=Chem.MolFromSmiles('C1(C2=NC=CC=C2)=CC=CC(C)=C1') ss = Chem.MolFromSmarts('c1c1C') frag = AllChem.DeleteSubstructs(m, ss) print(Chem.MolToSmiles(frag)) Pavel. On 03/31/2017 07:41 AM

[Rdkit-discuss] CanonicalRankAtoms

2017-03-28 Thread Pavel Polishchuk
[*])c.CN([*])C [0, 0, 3, 2] Kind regards, Pavel. -- Check out the vibrant tech community on one of the world's most engaging tech sites, Slashdot.org! http://sdm.link/slashdot _

Re: [Rdkit-discuss] FindAtomEnvironmentOfRadiusN

2017-03-27 Thread Pavel Polishchuk
? Pavel. On 03/27/2017 03:53 PM, Peter Gedeck wrote: Hello, The atom numbers start with 0. From the middle atom, there are no environments with radius 2. You will get a result if you use the first (=0) or the last (=2) atom. Try this: m = Chem.MolFromSmiles("NC

[Rdkit-discuss] FindAtomEnvironmentOfRadiusN

2017-03-27 Thread Pavel Polishchuk
1) Chem.MolToSmiles(Chem.PathToSubmol(m, i)) returns "" In the latter case I expected the same output "NCO". Were my expectations mistaken? Kind regards, Pavel. -- Check out the vibrant tech community on one

Re: [Rdkit-discuss] delete a substructure

2017-03-08 Thread Pavel Polishchuk
be reasonable to revise SMARTS to avoid such overlapping or create a list of rules (maybe hierarchical) which will define valid and not valid overlappings. Pavel. On 03/08/2017 06:32 PM, Chenyang Shi wrote: Dear Hongbin, I tried your method on a molecule, 4-Methylsalicylic acid (CC1=CC(=C(C=C1)C

Re: [Rdkit-discuss] Question about generating configurational isomerism

2017-01-27 Thread Pavel Polishchuk
Hi Jacob, you need to call AssignStereochemistry with force=True parameter Chem.AssignStereochemistry(mol, force=True) Pavel. On 01/28/2017 05:43 AM, Jacob Durrant wrote: I'm trying to set the configuration of a molecule with a double bond, but it doesn't seem to be working. Here's my

Re: [Rdkit-discuss] multiline legend in MolsToGridImage

2016-12-20 Thread Pavel Polishchuk
Thank you, that helped! I noticed that I need to substantially increase the height of the image to make the second line visible. If I'll want to draw more lined I should increase the height further. Pavel. On 12/20/2016 03:14 PM, Peter Gedeck wrote: Hello, I thought we had removed all

Re: [Rdkit-discuss] Generating all stereochem possibilities from smile

2016-12-09 Thread Pavel Polishchuk
risk :) If you will find mistakes let me know. Pavel. On 12/10/2016 03:44 AM, Brian Cole wrote: So I may need a little guidance on this one from someone with a little more historical knowledge of RDKit. I found the findPotentialStereoBonds function inside Chirality.cpp that appears to do what we

Re: [Rdkit-discuss] Extracting SMILES from text

2016-12-02 Thread Pavel
. Pavel. On 12/02/2016 11:11 AM, Greg Landrum wrote: An initial start on some regexps that match SMILES is here: https://gist.github.com/lsauer/1312860/264ae813c2bd2c27a769d261c8c6b38da34e22fb that may also be useful On Fri, Dec 2, 2016 at 11:07 AM, Alexis Parenty <alexis.parenty.h...@gmail.