Hello,
When a molecule is fragmented with Chem.FragmentOnBonds, how can I check if
there still is a path between two given atoms?
For example, this snippet breaks C1-C2 and C3-C4 bonds in cyclohexane,
producing two fragments:

from rdkit import Chem
cx=Chem.MolFromSmiles('C1CCCCC1')
bonds=(cx.GetBondBetweenAtoms(1,2).GetIdx(),
cx.GetBondBetweenAtoms(3,4).GetIdx())
fragmented=Chem.FragmentOnBonds(cx,bonds)

Is there a way in RDKit to check if there is a path between, say, C1 and
C2, without having to manually traverse the whole molecular graph? I guess
the alternative is to load everything into a graph library like networkx...

Best wishes,
Michal Krompiec
Merck KGaA
_______________________________________________
Rdkit-discuss mailing list
Rdkit-discuss@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/rdkit-discuss

Reply via email to