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.
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 Sw
Hi Enrico,
you may look at this script
https://github.com/DrrDom/rdkit-scripts/blob/master/rmsd_rdkit.py
It takes PDBQT as input and reference files and calc rmsd for the
largest common fragment if structures do not match each other. If there
are multiple models in the input file it will c
Dear colleagues,
after covid "holidays" we are cordially invite you on the 6th Drug
Design workshop which will be held 30 January - 3 February 2023 in
Olomouc (Czech Republic). It is focused on practical applications of
different chemoinformatic tools and approaches for drug development.
Th
Hello,
I use RDKit to embed initial conformations for docking. The issue is
with saturated rings. I can use a single random conformer but its
geometry may be unsuitable and the whole molecule will fail to dock. I
can use several starting conformers for docking and to avoid docking of
very s
Hi, Guillaume,
Idon't fully understand your goal.
You may fragment a molecule using certain rule(s). Rules may be
represented as positive and negative SMARTS patterns for bonds. Then you
can manually cleave only those bonds which correspond to positive SMARTS
patterns but not to negative o
I just want to share my script, which I use for enumeration of
stereoisomers. Enumeration of double bonds was added quite recently and
thus I didn't test it extensively.
I put it on github: https://github.com/DrrDom/rdk
It seems to work well on quite complex queries like
CC1CCC(CC1)C1CC=C(C\C(=
Hi,
I try to print multiline legends for molecules in a grid and use the
following code from notebook
from rdkit import Chem
from rdkit.Chem import Draw
# from rdkit.Chem.Draw import IPythonConsole
mols = [Chem.MolFromSmiles(s) for s in ["CCC", "C"]]
legends = ["1\nCCC", "2\nCCC
s = s.replace('\r\n', '\n')
s = s.replace('\n\r', '\n')
s = s.replace('\r', '\n')
lines = s.split('\n')
On Tue, Dec 20, 2016 at 2:45 AM, Pavel Polishchuk
mailto:pavel_polishc...@ukr.net>
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 cod
You might find this link useful -
http://www.rdkit.org/docs/GettingStartedInPython.html#chemical-transformations
However, the issue in your case is SMARTS definitions. If one SMARTS
completely covers another one it would be difficult to understand is it
artifact or not.I think it might be reas
Dear RDKitters,
I found the issue with FindAtomEnvironmentOfRadiusN but this can be a
feature. However, I did not findthis information in help and did not
expect such behavior.
If I apply FindAtomEnvironmentOfRadiusN function to a small molecule
and specify the radius greater than the siz
;
''
Is this more intuitive to you?
Best,
Peter
On Mon, Mar 27, 2017 at 9:35 AM Pavel Polishchuk
mailto:pavel_polishc...@ukr.net>> wrote:
Dear RDKitters,
I found the issue with FindAtomEnvironmentOfRadiusN but this
can be a
feature. However, I did not findt
SSERT(pth.size() == 3);
pth = findAtomEnvironmentOfRadiusN(*mol, 2, 1);
TEST_ASSERT(pth.size() == 4);
pth = findAtomEnvironmentOfRadiusN(*mol, 3, 1);
*TEST_ASSERT(pth.size() == 0);*
delete mol;
}
On Mon, Mar 27, 2017 at 10:00 AM Pavel Polishchuk
mailto:pavel_polishc...@ukr.net>> wrote
Hi,
I experimented a little bit with CanonicalRankAtoms and observed some
unexpected results.
I have two mols (actually sets of fragments):
C[*].n[*].C[*].N[*]
CC[*].CC[*].cn([*])c.CN([*])C
In the first case, pairs of carbons and nitrogens are recognized as
symmetrical
[0, 2, 0, 2]
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, Po
Hello,
I found two issues which are reproducible in 2016.09 and 2017.3 conda
builds.
0. The first issue is related to PostgreSQL. If I try to install it
from conda it asks me to downgrade rdkit to 2016.3. After downgrading of
rdkit psql works well.
1. Generation of 2D pharmacophore f
Hello,
is it possible to store custom fingerprints in psql DB and use them
for similarity search? And how to do this?
I foundtwo commands bfp_to_binary_text(bfp) and
bfp_from_binary_text(bytea)in RDKit cartridge but cannot understand how
to use them.
I want to store pharmacophore finge
indices speed up similarity search or what is a reason to add them?
Pavel.
On 04/26/2017 02:36 PM, Greg Landrum wrote:
Hi Pavel,
On Wed, Apr 26, 2017 at 8:41 AM, Pavel Polishchuk
mailto:pavel_polishc...@ukr.net>> wrote:
Hello,
is it possible to store custom fingerprints i
Hi,
I cannot solve an issue and would like to ask for an advice.
If there are different map numbers for attachment points for the same
fragment different canonical smiles are generated.
I observed such behavior only for fragments with 3 attachment points.
Below is an example.
I'm lookin
e)
Output:
S=c1c([*:1])c(Cl)[nH]c([*:3])c1[*:2]
S=c1c([*:1])c(Cl)[nH]c([*:3])c1[*:2]
S=c1c([*:1])c(Cl)[nH]c([*:3])c1[*:2]
S=c1c([*:1])c(Cl)[nH]c([*:3])c1[*:2]
S=c1c([*:1])c(Cl)[nH]c([*:3])c1[*:2]
S=c1c([*:1])c(Cl)[nH]c([*:3])c1[*:2]
Now, if you want the atomMaps in 1...2...3 output order, we
I do not know the valence model of RDkit. Therefore below there is an
ugly but working solution: set to the S+ atom NoImplicit.
mol = Chem.MolFromSmiles('CC(C)(C)SC')
rxn =
AllChem.ReactionFromSmarts('[CH0:1][S:2][CH3:3]>>[C:1][SH0+:2].[CH3-:3]')
ps = rxn.RunReactants((mol,))
print(Chem.MolT
Hello,
probably this is a message mainly to developers.
I discovered some strange behavior of removeHs function applied for
'[H][*:1]' molecule.
If I create that mol from smiles, RemoveHs does not remove a single H
from the mol:
mmm = Chem.MolFromSmiles('[H][*:1]')
Chem.MolToSmiles(Ch
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 in
Hi all,
canonicalization of fragment SMILES does not work properly. Below
there are two examples of identical fragments. The only difference is
the order of atoms (indices). However, it seems that RDKit
canonicalization does not take into account atom types.
Does someone have an idea how
ts. The
canonicalization code was not really designed to be used in situations
like this.
-greg
[1] The details of the canonicalization algorithm, including the
contents of the atom invariants, are described here:
http://dx.doi.org/10.1021/acs.jcim.5b00543
On Tue, Aug 1, 2017 at 2:
,
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
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
Hello,
I want to generate conformers for a stereoisomeric sugar moiety. The
code below works (loads proccesor) but returns none of them.
But if I remove all stereoconfiguration info in input SMILES the code
generates conformers. Playing with this issue I noticed that conformers
are generat
Hi all,
I have an issue with ConstrainedEmbed and I cannot figure out what
exactly causes this.
I have a molecule C[C@@H]1C1=O with 3D coordinates in 1.mol file
(attached). And I want to generate coordinates for another structure
with this core -
C[C@@H]1CC[C@H](O)CC1=O.
This is us
up and couldn’t finish
embedding without some help. Hope
someone with more in-depth insight can help here.
Anyway, for a walk around, adding H seems to do the trick:
In [39]: mol = AllChem.AddHs(mol_child)
In [40]: AllChem.EmbedMolecule(mol)
Out[40]: 0 # worked
In [41]: AllChem.ConstrainedE
31 matches
Mail list logo