Re: [Rdkit-discuss] [Rdkit-announce] RDKit 2018.03.2 release

2018-06-09 Thread Lukas Pravda
Hi Drew, are you sure that you are installing 64bit version? This https://anaconda.org/rdkit/rdkit says that the version you have installed is present in win-32 build if I am not mistaken. Best, Lukas From: Drew Gibson via Rdkit-discuss Reply-To: Drew Gibson Date: Saturday, 9 June

[Rdkit-discuss] Generate depiction matching 2D structure

2018-01-11 Thread Lukas Pravda
Dear all, I’ve just recently started using rdkit in python. Btw. A very nice piece of work. First, I’d like to generate 2D depiction of molecules I’m constructing from 3D coordinate data. I’m aware of methods such as GenerateDepictionMatching2DStructure(…) and

[Rdkit-discuss] Programatic access to the mol sanitation process results

2018-02-22 Thread Lukas Pravda
Dear rdkiters, I’m constructing molecules from scratch using python 3.5.4 and RDKit 2017.09.2 and due to the variety of reasons some of them are violating general principles of chemistry in a way implemented in rdkit, so I’m getting information like: Explicit valence for atom # 14 N, 4,

[Rdkit-discuss] FW: Protein Data Bank in Europe is looking for bioinformaticians

2018-08-07 Thread Lukas Pravda
here: https://www.ebi.ac.uk/pdbe/about/jobs Cheers, Lukas Pravda, Ph.D. Bioinformatician/Scientific Programmer Protein Data Bank in Europe (PDBe) European Bioinformatics Institute (EMBL-EBI) Wellcome Trust Genome Campus, Hinxton, Cambridge CB10 1SD United Kingdom

Re: [Rdkit-discuss] Programatic access to the mol sanitation process results

2018-03-09 Thread Lukas Pravda
: Greg Landrum <greg.land...@gmail.com> Date: Thursday, 22 February 2018 at 13:32 To: Lukas Pravda <lpra...@ebi.ac.uk> Cc: RDKit Discuss <Rdkit-discuss@lists.sourceforge.net> Subject: Re: [Rdkit-discuss] Programatic access to the mol sanitation process results Hi Lukas, On

Re: [Rdkit-discuss] Programatic access to the mol sanitation process results

2018-03-16 Thread Lukas Pravda
)     error_msg = out_stream.getvalue().decode('utf-8')     print(error_msg) Lukas From: Peter Gedeck <peter.ged...@gmail.com> Date: Friday, 9 March 2018 at 15:02 To: Lukas Pravda <lpra...@ebi.ac.uk> Cc: Greg Landrum <greg.land...@gmail.com>, <Rdkit-discuss@lists.sourceforge.net

[Rdkit-discuss] How to set bond width with use of MolDraw2DSVG

2018-10-30 Thread Lukas Pravda
Hi all, First of all, my configuration is following: macOS: 10.14 conda: 4.5.11 python: 3.6.6 rdkit: 2018.09.01 I just tried to set bond width for 2D SVG images and I run into situation I don’t understand much (nothing surprising ☺). I’m aware of 2 ways how to generate SVG images

[Rdkit-discuss] Coordgen library questions

2018-10-09 Thread Lukas Pravda
Hi all, I’m playing with the Coordgen library inside rdkit and I have a couple of questions I could not figure out by myself. Hopefully someone more experienced will know. [comment] The way one can pass a scaling factor to the bond size is very unintuitive. If I don’t provide any

Re: [Rdkit-discuss] Warning as error

2019-01-21 Thread Lukas Pravda
Hi Jean-Marc, Just a thought, but SDMolSupplier has a lazy eval, if I am not mistaken. Technically you should get all the rdkit warnings and errors at the time of processing that bit of the sdf file. You can always read the stderror output, parse it and throw exception every time a 'funny'

Re: [Rdkit-discuss] Bond tags in SVGs

2018-12-04 Thread Lukas Pravda
: Greg Landrum Date: Tuesday, 4 December 2018 at 17:24 To: Lukas Pravda Cc: RDKIT mailing list Subject: Re: [Rdkit-discuss] Bond tags in SVGs Hi Lukas, There's not currently a way to do this at the moment. The closest you can get is by calling AddMoleculeMetadata(): In [6]: d

[Rdkit-discuss] Bond tags in SVGs

2018-12-03 Thread Lukas Pravda
Hi all, I was wondering if there is a way how you can tag elements (bonds) in the svg created by rdkit. i.e. transform something like this: Into: Or similar. I’ve found possibility of tagging atoms in the SVG using Draw.rdMolDraw2D.MolDraw2DSVG.drawOptions() method that

[Rdkit-discuss] PDBe-KB ligand centric pages survey

2019-04-04 Thread Lukas Pravda
, including structural and functional annotations, domains, ligand-binding sites and interfaces. In the next step we would like to present a similar aggregated view from a small molecule / ligand perspective. Thank you for your time, Lukas -- Lukas Pravda, Ph.D. Bioinformatician/Scientific

Re: [Rdkit-discuss] Atom coordinates from PDB-file

2019-02-25 Thread Lukas Pravda
Hi Illimar, If you need to access coordinates without creating conformer object do you really need to use rdkit I the first place? PDB file is column based format, so extracting coordinates for atoms for example with python is very straightforward. Lukas -- Lukas

Re: [Rdkit-discuss] Bond tags in SVGs

2019-03-12 Thread Lukas Pravda
Hi Greg, I was wondering If you managed to create the Mac build you were talking about some time ago. Also I wonder If this functionality is going to be part of the next RDKit release? Best, Lukas From: Greg Landrum Date: Tuesday, 5 February 2019 at 14:45 To: Lukas Pravda Cc

Re: [Rdkit-discuss] Bond tags in SVGs

2019-02-05 Thread Lukas Pravda
at 17:49 To: Lukas Pravda Cc: RDKIT mailing list Subject: Re: [Rdkit-discuss] Bond tags in SVGs Hi Lukas, I had a chance to do a bit of work on this recently and I'd be interested to hear your feedback. Bonds are now tagged with their bond IDs (using classes) and the "Tag

Re: [Rdkit-discuss] Bond tags in SVGs

2019-02-05 Thread Lukas Pravda
If it is not too much trouble to ask, please build it for mac os (10.14.3) python 3.6.x. Thanks! Lukas From: Greg Landrum Date: Tuesday, 5 February 2019 at 13:40 To: Lukas Pravda Cc: RDKIT mailing list Subject: Re: [Rdkit-discuss] Bond tags in SVGs On Tue, Feb 5, 2019 at 12

Re: [Rdkit-discuss] GenerateDepictionMatching2DStructure question

2019-05-23 Thread Lukas Pravda
Hi Pat, >From my experience rdkit uses more or less 1.5 units bond length for 2D >depictions. So it makes sense if you rescale your template so that the bond >length is 1.5. This is the code snippet I use for the same thing to upscale template with bond lengths 1.0 to 1.5 import

[Rdkit-discuss] Get num of heavy atoms returns incorrect value

2019-05-01 Thread Lukas Pravda
Dear all, I construct my own rdkit.Mol objects from mmcif files. I wanted to use mol.GetNumAtoms(onlyExplicit=True) to get the number of heavy atoms for that molecule, however, I have noticed that the function returns all the time number of all atoms in the molecule including hydrogens (47

Re: [Rdkit-discuss] Get num of heavy atoms returns incorrect value

2019-05-01 Thread Lukas Pravda
(now explicitOnly) returns heavy atoms. Thanks Lukas From: Paolo Tosco Date: Wednesday, 1 May 2019 at 15:32 To: Lukas Pravda , RDKIT mailing list Subject: Re: [Rdkit-discuss] Get num of heavy atoms returns incorrect value mol.GetNumHeavyAtoms

Re: [Rdkit-discuss] Which method to prefer for computing 2D coordinates

2019-04-10 Thread Lukas Pravda
Evangelidis Date: Tuesday, 9 April 2019 at 15:43 To: Lukas Pravda Cc: Jose Manuel Gally , RDKIT mailing list Subject: Re: [Rdkit-discuss] Which method to prefer for computing 2D coordinates Hello Lukas, I am also struggling with 2D coordinate generation quite a long time as well as what

Re: [Rdkit-discuss] Which method to prefer for computing 2D coordinates

2019-04-09 Thread Lukas Pravda
Hi Jose, As you have shown there is no single method which would be perfect for everything. If you don’t care that much about speed, the possible solution could be to compute coordinates with all three approaches and then simply select the best conformer based on some criteria. The solution I

Re: [Rdkit-discuss] Fwd: sending warnings to stderr/stdout to files

2019-09-25 Thread Lukas Pravda
') Chem.WrapLogs() mol = Chem.MolFromSmiles('c1c1(C)(C)') log.close() sys.stderr = saved_std_err works Lukas From: Brian Lee Date: Wednesday, 25 September 2019 at 20:55 To: Cc: RDKIT mailing list , Lukas Pravda Subject: Re: [Rdkit-discuss] Fwd: sending warnings to stderr

Re: [Rdkit-discuss] Stereochemistry in rdkit

2019-10-30 Thread Lukas Pravda
atom tags (which randomly fixed those couple of issues, but apparently broke everything else down. I was wondering how does rdkit work out R/S from inchi string? Lukas From: Greg Landrum Date: Wednesday, 30 October 2019 at 04:28 To: Lukas Pravda Cc: RDKIT mailing list Subject: Re

Re: [Rdkit-discuss] Constructing a mol object from a PDB ligand

2019-12-17 Thread Lukas Pravda
Hi IIllimar, I don’t really know what your use case is, so it may be completely useless. However, just to add my two cents, we've created a package that builds on the top of rdkit and parses PDB ligand definitions from cif files. You can find the package here:

Re: [Rdkit-discuss] Fwd: sending warnings to stderr/stdout to files

2019-09-24 Thread Lukas Pravda
Hi Mike, The following code works for me: import sys from io import StringIO from rdkit import Chem saved_std_err = sys.stderr log = sys.stderr = StringIO() Chem.WrapLogs() # do whatever you want with rdkit whatever_used_to_be_printed_by_rdkit_in_console_as_str =

Re: [Rdkit-discuss] RDKit installation problem

2020-08-02 Thread Lukas Pravda
Hi Sebastian, Quickly looking at the available builds in the rdkit conda channel (https://anaconda.org/rdkit/rdkit) it appears that you are pulling windows 32-bit version of rdkit. Perhaps this is caused by the fact that you use 32bit version of conda? Try installing 64-bit version of conda

[Rdkit-discuss] Manganese ion as a radical?

2020-07-27 Thread Lukas Pravda
Dear rdkit community, I’m not quite sure if this is more of an rdkit or a chemistry related question. I’d like to understand why a manganese ion has 3 radical electrons when interpreted by rdkit. I have not seen radicals in any other metal ion so far. The code to get the depiction looks