Dear Henrique, You can try something like this:
import rdkit from rdkit import Chem from rdkit.Chem import AllChem from rdkit.Chem.rdMolDescriptors import CalcCoulombMat smiles = 'CN1C=NC2=C1C(=O)N(C(=O)N2C)C' caffeine = Chem.MolFromSmiles(smiles, sanitize=True) caffeine = Chem.AddHs(caffeine) AllChem.EmbedMolecule(caffeine) AllChem.UFFOptimizeMolecule(caffeine) CM = CalcCoulombMat(caffeine) list(CM[0]) The CM object contains a list of Coulomb matrices calculated for the molecule that I think it is generated by randomizing the arrangement of atoms, but I am not totally sure. It seems that you need first to optimize the molecule to get its 3D conformation, and then you can pass this object as input to the CalcCoulombMat function. One way to access the values of each Coulomb Matrix is by transforming it into a list or a NumPy array. I hope this can help you. Max Pinheiro Jr --------------------------------------------- Postdoctoral researcher Aix-Marseille Université, France Institut de Chimie Radicalaire Em dom., 20 de set. de 2020 às 14:58, Henrique Castro < henrique...@outlook.com> escreveu: > Dear colleagues, how are you? > I have almost zero experience with RDKit, so forgive me for such a basic > request. > Could anyone provide me an example of the new Coulomb Matrix generator in > RDKit? I intend to iterate through my organometallic molecules stored in a > pandas DataFrame. > > Thank you in advance > > -- > Henrique C. S. Junior > > _______________________________________________ > Rdkit-discuss mailing list > 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
