Skip to main content

Using the CDK's group module

There is a new module and package for the CDK that is currently under review. This is a short guide to how to use it - to help both reviewers and users.

The basic idea is that molecules can be considered as a kind of graph, and that one useful thing to calculate about such graphs is the automorphism group that preserves element labels and/or bond labels. To put it another way, calculating the symmetries of the molecule - although I should point out that it's not quite the same as the crystallographic symmetry groups.

As a simple example, consider these two molecules (1,4-cylohexadiene and 4h-pyran) :


They are numbered from 0-5 for programming convenience; on the right each molecule has a table of automorphisms written as permutations in cycle notation. It should be fairly obvious that - for example - the H-Flip sends atom 0 to atom 4, 1 to 3, and fixes 2 and 5. Only the H-Flip is an automorphism for 4h-pyran, due to the oxygen atom.

The code to do this is fairly short and really just involves creating an AtomDiscretePartitionRefiner and then calling the getAutomorphismGroup(IAtomContainer) method. This returns a PermutationGroup which stores the automorphisms. What you do with them then is up to you...

There is a corresponding class to find automorphisms of the bonds of an atom container. This may be less useful, but here is napthalene as an example:


Note that the bonds are numbered, not the atoms; also the two different double-bond arrangements are called a and b for reference. The a form has only the V-Flip automorphism that swaps bonds (1, 2), (3, 10) and so on.

Finally, what are the actual uses in chemistry for this? Well, one possibility is external symmetry numbers (interesting reference, actually) - as also mentioned in this post. Another is molecule generation; it's used heavily in AMG. A future possibility might also be using it in CIP or other chirality code.

Comments

Unknown said…
There is already some related code in the CDK which does symmetry of atoms. You might be interested in this, it is the EquivalentClassPartitioner and the function getTopoEquivClassbyHuXu
gilleain said…
This comment has been removed by the author.
gilleain said…
Hmm. I thought it double-posted so I deleted the duplicate comment - and now it's gone..

Anyway, it was :

"Good point, Patrik - one of the strengths of the CDK is that it has multiple solutions. It's also one of the weaknesses!

There was also an ancient branch that had a class to find the symmetries from the 3D structure, that could be integrated somehow. It's a little difficult to make packages written by different authors to work neatly together without making large changes. Some sort of interface, perhaps... "
Anonymous said…
Appreciate the recommendation. Will try it out.

Popular posts from this blog

Adamantane, Diamantane, Twistane

After cubane, the thought occurred to look at other regular hydrocarbons. If only there was some sort of classification of chemicals that I could use look up similar structures. Oh wate, there is . Anyway, adamantane is not as regular as cubane, but it is highly symmetrical, looking like three cyclohexanes fused together. The vertices fall into two different types when colored by signature: The carbons with three carbon neighbours (degree-3, in the simple graph) have signature (a) and the degree-2 carbons have signature (b). Atoms of one type are only connected to atoms of another - the graph is bipartite . Adamantane connects together to form diamondoids (or, rather, this class have adamantane as a repeating subunit). One such is diamantane , which is no longer bipartite when colored by signature: It has three classes of vertex in the simple graph (a and b), as the set with degree-3 has been split in two. The tree for signature (c) is not shown. The graph is still bipartite accordin

Király's Method for Generating All Graphs from a Degree Sequence

After posting about the Hakimi-Havel  theorem, I received a nice email suggesting various relevant papers. One of these was by Zoltán Király  called " Recognizing Graphic Degree Sequences and Generating All Realizations ". I have now implemented a sketch of the main idea of the paper, which seems to work reasonably well, so I thought I would describe it. See the paper for details, of course. One focus of Király's method is to generate graphs efficiently , by which I mean that it has polynomial delay. In turn, an algorithm with 'polynomial delay' takes a polynomial amount of time between outputs (and to produce the first output). So - roughly - it doesn't take 1s to produce the first graph, 10s for the second, 2s for the third, 300s for the fourth, and so on. Central to the method is the tree that is traversed during the search for graphs that satisfy the input degree sequence. It's a little tricky to draw, but looks something like this: At the top

1,2-dichlorocyclopropane and a spiran

As I am reading a book called "Symmetry in Chemistry" (H. H. Jaffé and M. Orchin) I thought I would try out a couple of examples that they use. One is 1,2-dichlorocylopropane : which is, apparently, dissymmetric because it has a symmetry element (a C2 axis) but is optically active. Incidentally, wedges can look horrible in small structures - this is why: The box around the hydrogen is shaded in grey, to show the effect of overlap. A possible fix might be to shorten the wedge, but sadly this would require working out the bounds of the text when calculating the wedge, which has to be done at render time. Oh well. Another interesting example is this 'spiran', which I can't find on ChEBI or ChemSpider: Image again courtesy of JChempaint . I guess the problem marker (the red line) on the N suggests that it is not a real compound? In any case, some simple code to determine potential chiral centres (using signatures) finds 2 in the cyclopropane structure, and 4 in the