Re: [Rdkit-discuss] RMSD value between two non-covalent molecular complexes

2017-06-22 Thread Francois BERENGER
On 06/22/2017 11:20 PM, gosia olejniczak wrote: Hi again, i found where the problem was (it seems): as i was reading in molecules from sdf file through "SDMolSupplier" by doing: suppl = Chem.SDMolSupplier(filename) the hydrogen atoms were removed (what was not obvious since e.g. GetNumAtom

Re: [Rdkit-discuss] RMSD value between two non-covalent molecular complexes

2017-06-22 Thread gosia olejniczak
Hi again, i found where the problem was (it seems): as i was reading in molecules from sdf file through "SDMolSupplier" by doing: suppl = Chem.SDMolSupplier(filename) the hydrogen atoms were removed (what was not obvious since e.g. GetNumAtoms with "heavy" flag was returning the total number of

Re: [Rdkit-discuss] RMSD value between two non-covalent molecular complexes

2017-06-21 Thread Paolo Tosco
Hi Gosia, you may want to look at the CombineMols() function in the rdmolops module. If you make a single molecule out of the two fragments, then you should manage to get what you are looking for. Cheers, p. On 06/21/17 14:47, gosia olejniczak wrote: dear RDKit experts, is it possible to c