Many thanks for your examples, Wim. But I was checking the option to save the labels D and T in the molfile for the hydrogen isotopes, as other tools can do.
Regards Santiago [http://www.mestrelab.com/mestrelab/wp-content/uploads/signs/line3.jpg] [http://www.mestrelab.com/mestrelab/wp-content/uploads/signs/M-red-200pxb.jpg]<http://www.mestrelab.com> SANTIAGO FRAGA Software Developer santi...@mestrelab.com<mailto:%20santi...@mestrelab.com> MESTRELAB RESEARCH S.L. PHONE +34881976775 FAX +34981941079 Feliciano Barrera, 9B-Bajo 15706 Santiago de Compostela (SPAIN) Follow us: [Mestrelab Twitter]<https://twitter.com/mestrelab> [Mestrelab Linkedin] <https://www.linkedin.com/company/mestrelab-research> [Canal de YouTube Mestrelab] <https://www.youtube.com/channel/UCf3MVnd3XZflv0acvTv14ww> [MestreBlog] <http://mestrelab.com/blog/> ________________________________ De: Wim Dehaen <wimdeh...@gmail.com> Enviado: lunes, 10 de abril de 2023 18:07 Para: Santiago Fraga <santi...@mestrelab.com> Cc: rdkit-discuss@lists.sourceforge.net <rdkit-discuss@lists.sourceforge.net> Asunto: Re: [Rdkit-discuss] Deuterium/Tritium labels in Molfile rdkit outputs a molfile with correct isotope labels for me using just: mol=Chem.MolFromSmiles("[3H]c1ccccc1[2H]") Chem.MolToMolFile(mol,"test.mol") or labelling the atoms post hoc: mol=Chem.MolFromSmiles("c1ccccc1") mol=Chem.AddHs(mol) mol.GetAtomWithIdx(6).SetIsotope(3) mol.GetAtomWithIdx(7).SetIsotope(2) mol=Chem.RemoveHs(mol) Chem.MolToMolFile(mol,"test2.mol") I hope this helps best wishes wim On Mon, Apr 10, 2023 at 4:43 PM Santiago Fraga <santi...@mestrelab.com<mailto:santi...@mestrelab.com>> wrote: Good afternoon! I am a relatively new user of RDKit, and mainly the C++ API. I am trying to save in a molfile the labels D and T for the hydrogen isotopes. Like in the following molfile: MJ230401 8 8 0 0 0 0 0 0 0 0999 V2000 -0.3572 0.4125 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 -1.0716 0.0000 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 -1.0716 -0.8250 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 -0.3572 -1.2375 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 0.3572 -0.8250 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 0.3572 0.0000 0.0000 C 0 0 0 0 0 0 0 0 0 0 0 0 -0.3572 1.2375 0.0000 T 0 0 0 0 0 0 0 0 0 0 0 0 1.0717 0.4125 0.0000 D 0 0 0 0 0 0 0 0 0 0 0 0 3 4 2 0 0 0 0 4 5 1 0 0 0 0 5 6 2 0 0 0 0 6 1 1 0 0 0 0 1 2 2 0 0 0 0 2 3 1 0 0 0 0 6 8 1 0 0 0 0 1 7 1 0 0 0 0 M END I am trying to set directly the labels in the hydrogen atoms: atom->setProp<string>("atomLabel", "D"); or     atom->setProp<string>("_displayLabel", "D"); But when the molfile is generated the labels are not transferred. It seems also that when reading a mofile including the labels, they are discarded. Many thanks in advance Santiago Fraga _______________________________________________ Rdkit-discuss mailing list Rdkit-discuss@lists.sourceforge.net<mailto:Rdkit-discuss@lists.sourceforge.net> https://lists.sourceforge.net/lists/listinfo/rdkit-discuss
_______________________________________________ Rdkit-discuss mailing list Rdkit-discuss@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/rdkit-discuss
