Re: [Rdkit-discuss] Help mapping atoms between two files

2020-04-04 Thread Paolo Tosco
Hi Gustavo, that's because a2.pdb has formal charges whereas a1.pdb has none. If you neutralize a2 then you'll find a match. You could use |Chem.MolStandardize.rdMolStandardize.||Uncharger| or something as below: from rdkit import Chem a1 = Chem.MolFromPDBFile("a1.pdb") a2 = Chem.MolFromPDBF

[Rdkit-discuss] Help mapping atoms between two files

2020-04-04 Thread Gustavo Seabra
HI all, I'm trying to use get the substructure matches between two different PDB files with the same molecule, but different atom order and naming. However, GetSubstructMatches Just returns nothing, i.e. no matches (files attached): For example: >>> ref_mol = Chem.MolFromPDBFile(str("a1.pdb")) >