Skip to main content

Continuous symmetry and chirality measures: approximate algorithms for large molecular structures


Quantifying imperfect symmetry of molecules can help explore the sources, roles and extent of structural distortion. Based on the established methodology of continuous symmetry and chirality measures, we develop a set of three-dimensional molecular descriptors to estimate distortion of large structures. These three-dimensional geometrical descriptors quantify the gap between the desirable symmetry (or chirality) and the actual one. They are global parameters of the molecular geometry, intuitively defined, and have the ability to detect even minute structural changes of a given molecule across chemistry, including organic, inorganic, and biochemical systems. Application of these methods to large structures is challenging due to countless permutations that are involved in the symmetry operations and have to be accounted for. Our approach focuses on iteratively finding the approximate direction of the symmetry element in the three-dimensional space, and the relevant permutation. Major algorithmic improvements over previous versions are described, showing increased accuracy, reliability and structure preservation. The new algorithms are tested for three sets of molecular structures including pillar[5]arene complexes with Li+, C100 fullerenes, and large unit cells of metal organic frameworks. These developments complement our recent algorithms for calculating continuous symmetry and chirality measures for small molecules as well as protein homomers, and simplify the usage of the full set of measures for various research goals, in molecular modeling, QSAR and cheminformatics.


Symmetry is an eye-catching phenomenon that expresses the beauty and mystery of nature. In chemistry, it is frequently perceived as a driving force that controls the shape of molecular structures, defines selection rules for the interaction of light and matter, and determines the terms and mechanisms of chemical reactions [1]. While perfect symmetry is conceptually appealing, numerous experimental and computational studies show that actual structures are only approximately symmetric [2,3,4,5,6]. Various reasons lead to distortion of otherwise symmetric molecules, including conformational flexibility, dynamics, chemical processes, physical conditions, crystallization conditions and the chemical environment. In many of these cases, describing the molecules at hand from the perspective of their deviation from the original symmetric geometry provides deeper understanding of the molecular systems. Such a treatment can highlight anomalous cases, and shed light on mechanisms of symmetry breaking. Mathematically, the description is based on the treatment of symmetry as a continuous quantitative parameter of the molecular structure, rather than a binary yes/no property. This approach was developed by Avnir and coworkers in the early 90's of the twentieth century, in the form of the continuous symmetry measure (CSM) [7,8,9], the continuous chirality measure (CCM) [10], and the continuous shape measure (CShM) [11]. In recent years we had considerably improved the algorithms that calculates the CSM and CCM for small-to-medium sized molecules [12] as well as for protein homomers [13]. Here we extend our methodology to large molecular systems which are beyond the reach of previous algorithms.

Given a distorted molecular structure and a symmetry point group G, the CSM algorithm searches for the nearest structure that belongs to G and maintains the same connectivity as the original molecule. The distance between the original structure and the nearest symmetric structure defines the symmetry measure. The continuous chirality measure follows by calculating the minimum CSM with respect to all the achiral point groups, Sn. In the last three decades, this set of symmetry and chirality measures were extensively applied for describing various chemical phenomena in studies related to the crystal structure of inorganic compounds [3, 14, 15], reaction paths and reactivity [16,17,18], dynamics and temperature [19,20,21], protein structure and activity [22,23,24], quantitative structure–activity relationship (QSAR) [25, 26], and many more. Beyond measuring symmetry and chirality, the CSMs and CCM can be used to estimate structural elongation, planarity and conformational flexibility, for both symmetric and asymmetric structures [20, 24]. In addition, the method was applied to other fields such as image processing [27] and archeology [28, 29].

Finding the nearest symmetric structure is the main challenge of the CSM method, as it is unknown a priory, and may change even between conformers of the same molecule. This property turns CSMs into powerful global descriptors of the three-dimensional structure capable of distinguishing between various conformers of the same molecule. This ability stands at the heart of the concept of symmetry maps, an analysis tool for distorted structures [30, 31]. From a mathematical perspective, a structure that belongs to G is represented by a three dimensional vector (which represents the direction of the symmetry operation) and a permutation of the set of atoms (which represents the action of the symmetry operation on the molecule's atoms). As the number of atoms increases, the number of possible permutations increases as well, and calculating the CSMs becomes computationally intensive. To overcome this obstacle, we improved the algorithm for small-to medium-sized structures, by utilizing the connectivity map of the molecule to scan only structure-preserving permutations [12]. This improvement increased the accuracy and speed of the calculation, making the method applicable to much larger molecules compared to the original algorithm [8].

For very large molecular systems, with branched structures and complex connectivity maps, scanning all the permutations, even only the ones that preserve the connectivity map, becomes computationally intensive. For this purpose, Dryzun et al. [9] developed an approximate algorithm to calculate CSMs. Instead of searching for all possible permutations, their algorithm iteratively searches for an approximate direction of the symmetry element and its related permutation, until convergence is reached. While this approach is relatively fast, it may result with a permutation that breaks the connectivity map of the structure. To overcome this obstacle, we recently modified the method for protein homomers, utilizing the amino acids sequence to reduce the size of symmetry-equivalent groups of atoms, and force the code to preserve both the sequence and the chains structure. The Hungarian algorithm [32] was applied to efficiently solve the assignment problem and find the best permutation [13]. This approach is based on our prior knowledge of the protein sequence, and may become less effective when such information is absent, for example in the case of supramolecular structures and nanostructures. Finding the nearest symmetric structure in such cases without losing information on their connectivity maps requires a different methodology. In this study we present several algorithms for approximate calculation of symmetry and chirality measures that differ by their level of structure preservation and efficiency. Our set of algorithms provides a comprehensive toolkit for structural analysis that can be used to explore different sources of distortion, including conformational and topological distortion. In what follows we review the CSM methodology, present the details of the new algorithms, and use them to analyze three sets of molecules with various levels of approximate symmetry: pillar[5]arenes, C100 fullerenes and large unit cells of metal organic frameworks (MOFs).


Review of the CSM method

Let us briefly review the fundamentals of the CSM methodology [7, 8, 12]. We consider a given molecule A of N atoms, that belongs to the symmetry point group G, where G is either Cn (n = 2, 3, 4, 5,…) or Sn (n = 1, 2, 4, 6,…). Recall that by definition, S1 = Cs and S2 = Ci. Let \({{\mathbf{Q}}} = \left\{ {{{\mathbf{Q}}}_k :\,\,1 \le k \le N} \right\}\) be the set of coordinate vectors of the molecule's atoms, and let \({{\mathbf{Q}}}_0 = \frac{1}{N}\sum_{k = 1}^N {{{\mathbf{Q}}}_k }\) be its geometric center of mass. We are looking for a symmetry operation \(T\), which generates a cyclic point group of type G. Note that T is a rotation (either proper or improper) by an angle of \(360^{\circ} / n\). In both cases, \(T\) is determined by a 3-dimensional direction vector, which we denote by \(\nu_{sym}\). The continuous symmetry measure (CSM) is defined by:

$$S\left( G \right) = 100 \cdot M\left( G \right)/D$$


$$M(G) = \min \left[ {\sum_{k = 1}^N {\left| {{\mathbf{Q}}_k - {\mathbf{P}}_k} \right|^2 } } \right];\quad D = \sum_{k = 1}^N {\left| {{\mathbf{Q}}_k- {\mathbf{Q}}_0 } \right|^2 }$$

and the minimum is over all the symmetric (i.e. \(T\)-invariant) structures \(\left\{ {{{\mathbf{P}}}_k :\,\,1 \le k \le N} \right\}\) and all possible direction vectors \(\nu_{sym}\). Equivalently,

$$M(G) = \frac{1}{2n}\min \sum_{i = 1}^n {\sum_{k = 1}^N {\left| {T^i {{\mathbf{Q}}}_k - {{\mathbf{Q}}}_{\pi^i (k)} } \right|^2 } }$$

where \(\pi\) is a permutation of the set of atoms \(\left\{ {1,2,\ldots,N} \right\}\) which preserves the atom types and the molecule's connectivity map, and the cycles of \(\pi\) are of size 1, 2, or n. Note that the value of 2 is only allowed when G = Sn or G = C2.

In our previous work [12, 13] we developed two algorithms for evaluating the CSM: The first one finds the exact value of M(G) in Eq. (3) by enumerating the structure preserving permutations, that is, permutations that satisfy the condition: π(i) ↔ π( j) if and only if i ↔ j, for all pairs of atoms (i, j). Here we denote by \(i \leftrightarrow j\) the existence of a bond between atoms i and j. This algorithm, which we will call here the exact algorithm, shows excellent performance for small and medium sized molecules, and for larger molecules with a small number of structural symmetries [12]. For example, for fullerene, C60, all the atoms are in the same equivalence class. There are 2.73 × 1043 permutations that define a C2 operation, but only 32 of them preserve the structure of the molecule [12].

The second method focuses on protein homomers. Here the number of atoms does not allow scanning all structure preserving permutations. We therefore find an approximation for the CSM value by performing permutation-direction iterations (as explained below) [13]. A partial reduction of the number of permutations is achieved by exploiting the amino-acids sequence to define equivalence classes based on the types of the atom as well as its residue's designation and sequence number. Thus, a Cα of alanine with sequence number 3 on chain A can only be permuted with Cα of alanine with sequence number 3 on another chain, but not with Cα of different residues, or other alanine residues with different sequence numbers. The rest of the permutation is found using the Hungarian algorithm [32]. Further improvements are achieved by exploiting the protein polypeptide chain structure, making sure that the permutation does not break the chains, and carry each chain in its entirety to another chain: for protein homomers with more than two chains, the Hungarian method is used at the chain level as well, in order to find bijections between the different chains and piece them together to a complete permutation [13].

Methods for general large molecules

As molecules become larger, the size of the equivalence classes increases, and the feasibility of a CSM calculation that performs an exhaustive search over all possible permutations reduces considerably. Unlike proteins, the atoms of large molecular structures do not have sequence identifiers that can help reduce the number of possible permutations. The challenge with such molecules is to find a good approximation to the CSM within a reasonable calculation time. We have developed several algorithms that successfully face this challenge, as described below.

Permutation-direction iterations

Our algorithms are based on the process, first described by Dryzun et al. [9], of going back and forth between estimating the direction vector \(\nu_{sym}\) and estimating the permutation \(\pi\). An initial guess for a direction vector is based on the best line or plane that fits the geometric centers of all the equivalence groups of the molecule at hand, after outliers are removed. Two additional perpendicular direction vectors are taken into account in order to equally span the three-dimensional space [9].

Next, we perform an iterative process in which we update each time the permutation \(\pi\) and the vector \(\nu_{sym}\):

Given the current vector \(\nu_{sym}\), the permutation is estimated by using a greedy algorithm. One first calculates the distance matrix \(A = \left( {A_{ij} } \right)\), where

$$A_{ij} = \left| {TQ_i - Q_j } \right|^2 \quad {\mathrm{for}}\quad 1 \le i \le N,\quad 1 \le j \le N$$

and \(T\) is the symmetry operation corresponding to the vector \(\nu_{sym}\).

Let the smallest entry in the matrix \(A\) be \(A_{i_0 j_0 }\). The permutation value \(\pi (i_0 )\) is set to \(j_0\), and the row and column of this entry are greyed out. Then, the smallest entry from the remaining values in the matrix is chosen, to set another permutation value, and so on.

Given the permutation, the optimal vector for this permutation \(\nu_{sym}\) is found by an exact analytical method [9, 12], which later serves to find a new permutation and so on, until convergence is achieved.

The Hungarian algorithm

A more sophisticated approach to finding the permutation, taken in our previous work [13] applies the Hungarian algorithm [32] to the matrix \(A\). The algorithm finds a permutation \(\pi\) for which the sum \(\sum_{i = 1}^N {A_{i\pi (i)} }\) is minimal. This approach takes into account the interaction between different choices of the permutation values.

We should note that both approaches—the greedy algorithm and the Hungarian algorithm [32]—yield permutations that are not guaranteed to preserve the structure of the molecule. In practice, this depends on the symmetry level of the molecule. For molecules of high symmetry level, the approximate algorithms will tend to find the right permutation, and consequently, the bonding structure will be preserved, either fully or with a high preservation rate. For highly asymmetric molecules, this will not always hold, as permutations which do not preserve the bonding structure could attain lower values of the target function, thereby creating imprecisions in the CSM calculation. We elaborate more on this issue in the Results section.

The Fibonacci lattice algorithm

The original algorithm of permutation-vector iterations [9] starts with a direction vector which is an educated guess, hoping that through the iterations, the direction vector will converge to the optimal one. However, there is a risk that the iterative process will converge to a local minimum and not to the global one. In order to increase the reliability of the method, we devised the following strategy: A set of \(k\) vectors on the unit sphere, \(S\), is chosen such that the points are as equally spaced as possible on the sphere. A sequence of permutation-vector iterations is performed, starting at each of the vectors in \(S\).

Finding a set of \(M\) unit vectors which are equally spaced on the sphere is, in general, an unsolved mathematical problem. However, there is a well known way to produce \(M\) unit vectors which are almost equally spaced. This set is called The Fibonacci lattice [33, 34]. It is defined as follows:

For any \(0 \le k \le M - 1\),

$$\theta = \pi \left( {\sqrt {5} + 1} \right)k;\quad x = 1 - \frac{2k}{{M - 1}};\quad r = \sqrt {1 - x^2 }$$
$$u_k = (x,r\cos \theta ,r\sin \theta )$$

The set \({{\mathcal{F}}}_M = \left\{ {u_k |0 \le k \le M - 1} \right\}\) is called the Fibonacci lattice of size \(M\). Figure 1 presents a Fibonacci lattice of size 500.Footnote 1

Fig. 1
figure 1

Fibonacci lattice of size 500

We use all of the vectors in \({{\mathcal{F}}}_M\) as starting points for the permutation-direction iterations. The implementation of this algorithm can be sped up by parallel computation: The set \({{\mathcal{F}}}_M\) is divided into subsets, and each processor performs the iterations with the vectors in its assigned subset as initial vectors.

The approximate structure preserving algorithm

As discussed above, scanning only structure preserving permutations is considerably advantageous over the original CSM algorithm [7, 8], making it faster and more accurate [12]. However, for large molecules its running time can be too long for a feasible or efficient calculation. On the other hand, the approach of the approximate permutation-direction iterations algorithm is operative on any molecular size, but can often yield permutations that do not preserve the structure. In these cases, the nearest symmetric structure, which serves as the reference structure for the CSM calculation, may lose its chemical essence making the CSM result less informative. We now present a new algorithm, which attempts to bridge these two approaches and have the benefits of both. The algorithm performs permutation-direction iterations as in the approximate algorithm, but instead of using the Hungarian algorithm [32] for estimating the permutation, it searches for the best structure preserving permutation with respect to the distances matrix (4). As detailed below, this algorithm uses information from the distances matrix to prioritize the search for the permutation.

Recall that the exact (structure preserving) algorithm for small molecules [12] performs a recursive enumeration of the permutations in the following way: at each step, we define a partial assignment of permutation values; these assignments determine restrictions on the permutation values of the other atoms, according to the connectivity map of the molecule and the expected cycle structure of the permutation. An atom with a minimal number of permutation options is chosen, and for each option, this permutation value is assigned and the algorithm continues recursively.

In our new algorithm, which we call the approximate structure-preserving algorithm, in each permutation-vector iteration, the permutation is chosen as follows: We build the distances matrix (4), and perform a search for a structure preserving permutation as in the exact algorithm, but in the recursive step we prioritize assignment values \(\pi (i) = j\) for which \(A_{ij}\) is small: Given the options \(j_1 ,\ldots,j_k\) for values of \(\pi (i)\), we sort them according to the value \(A_{ij}\) and perform the recursive calls from the minimum value and up. Whenever the algorithm reaches a full permutation \(\pi\), the sum \(S(\pi ) = \sum_{i = 1}^N {A_{i\pi (i)} }\) is calculated, and in the end, the permutation for which \(S(\pi )\) is minimal is returned. As soon as the first permutation is obtained, the corresponding \(S(\pi )\) is set as a threshold for the next partial permutations. For a partial permutation \(\mu\), defined for a set \(T\) of atoms, if the sum \(\sum_{i \in T}^{\,} {A_{i\mu (i)} }\) is already larger than the minimal value of \(S(\pi )\) so far obtained, there is no point in carrying on the search from \(\mu\). By prioritizing assignment values with small corresponding entries in the distance matrix, we increase the likelihood of finding the best permutation fast. The algorithm contains a time limit, after which the search is terminated.

Obtaining upper and lower bounds on the CSM

For large molecules, when the exact algorithm is unfeasible, it is advantageous to obtain upper and lower bounds for the CSM. Instead of bounding the expression \(M(G)\) in (2), we bound the related expression

$$\widehat{M}(G) = \frac{1}{2}\min \sum_{k = 1}^N {\left| {T{{\mathbf{Q}}}_k - {{\mathbf{Q}}}_{\pi (k)} } \right|^2 }$$

In this expression the average in (2) over \(i = 1,2,\ldots,n\) is replaced by the value for \(i = 1\). Note that for \(n \le 2\) (i.e. when G is Cs or Ci), \(\widehat{M}(G) = M(G)\), and for \(n > 2\), \(\widehat{M}(G)\) is a reasonable approximation of \(M(G)\). Indeed, the approximate algorithms settle for minimization of \(\widehat{M}(G)\) [9, 13].

Consider the approximate algorithm, paired with the Hungarian method [32] for finding the permutation, and the Fibonacci lattice method for the initial direction vectors, with a large number of initial points. We claim that this method produces a lower bound for \(\widehat{M}(G)\) (which is a lower bound to the CSM when G is Cn or Sn with \(n \le 2\)). Indeed, if the Fibonacci lattice is dense enough, one of its vectors will be close enough to the optimal direction vector; and the permutation found by the Hungarian algorithm [32] for this vector, will satisfy \(\frac{1}{2}\sum_{k = 1}^N {\left| {T{\mathbf{Q}}_k - {\mathbf{Q}}_{\pi (k)} } \right|^2 } \le \widehat{M}(G)\), as the Hungarian method finds the minimum over all permutations, not just the structure preserving ones.

Consider, on the other hand, the approximate structure-preserving algorithm. This algorithm calculates the minimum over some of the direction vectors \(v_0\), and some of the structure preserving permutations \(\pi\), of the expression \(\frac{1}{2}\min \sum_{k = 1}^N {\left| {T{{\mathbf{Q}}}_k - {{\mathbf{Q}}}_{\pi (k)} } \right|^2 }\). Its value is, therefore, bigger or equal to \(\widehat{M}(G)\).

Results and discussion

In order to assess the efficiency and accuracy of the algorithms we tested them with three sets of molecules with various levels and sources of distortion. The first two sets were chosen so that both the exact and approximate CSM calculations can be performed, in order to compare the results and evaluate the performance of the various approximate algorithms. The first consisted of a set of highly flexible conformers of pillar[5]arene complexed with a Li+ ion in the gas phase. The second consisted of isomers of C100 fullerene. Different combinations of pentagons and hexagons that characterize the topology of each isomer, create diverse ellipsoidal-like structures with various symmetries. As a third set, we analyzed the crystallographic unit-cells of several MOFs, each with thousands of atoms, to test the applicability of the approximate algorithms on very large structures. For each set we calculated symmetry and chirality measures, the level of structure preservation, and the time of calculation. Details of each analysis are described next.

Pillar[5]arene-Li+ complexes

Pillar [n]arene represents a class of supramolecular systems, originally synthesized by Ogoshi and coworkers [35, 36], with pillar-shaped architecture and double rim structure. Recent studies show that their superior host–guest abilities, planar chirality, and the ability to undergo self-assembly processes with highly symmetrical structures stem from their unique shape [37,38,39]. Pillar[5]arene, presented in Fig. 2, is the most common member of this family, for which five-fold rotational symmetry has been reported [36]. The high number of rotatable bonds gives rise to hundreds of conformers, which are generally close in energy. As will be shown below, the actual structure of these conformers can often be quite far from symmetry. In order to correctly describe the manifold of conformers of this molecular system, and distinguish between them, tools that provide a global description of the geometrical parameters are needed. Such description can contribute to better understanding of host–guest interactions, as was recently shown for crown ethers [40]. The size of the pillar[5]arene molecule, with 75 atoms, allows for CSM calculation with the exact algorithm. On the other hand, the flexibility leads to high distortion levels, making the evaluation of the CSM with the approximate permutation-direction iterations approach, quite challenging. Therefore, this system is suitable for benchmark purposes. Although the host itself can be distorted, for the purpose of testing our algorithms we enforced distortion even further and inserted a Li+ ion to the center of the host. Previously we showed that the symmetry of the related host, 18-crown-6, is strongly influenced by the presence of alkali ions, particularly Li+, and that this distortion is related to the host–guest binding energy [40]. The small size of the ion prevents direct interaction with all the oxygen atoms of the host, forcing the last to fold in order to increase its interaction with the guest. The cavity of pillar[5]arene is larger than 18-crown-6, and it is therefore reasonable to assume that its tendency to distort with respect to either the fivefold or twofold symmetries will be significant when a Li+ ion is used as a guest.

Fig. 2
figure 2

Pillar[5]arene with Li+ ion in a perfectly symmetric D5 conformation

Conformation analysis

A pillar[5]arene molecule with perfect D5 symmetry with Li+ ion in its center (Fig. 2), was subjected to conformation analysis calculation in the gas phase, using the LowModeMD algorithm [41] implemented in MOE [42]. The Amber10:EHT basis set was used with an energy window of 20 kcal/mol. This calculation created 647 conformers which were filtered to remove duplicates based on their energy, radius of gyration and the CSM with respect to C5, C2 and Cs point groups (were the last serves as chirality measure as explained below). The exact algorithm of the CSM was used for this purpose. We used thresholds of 0.0001 kcal/mol for the energy, 0.05 for the radius of gyration, and a relative threshold of 5% for the CSMs. If the differences between two conformers, for all parameters, were equal or below these thresholds, they were considered equal, and one of them was deleted. The lowest energy conformers within an energy window of 8 kcal/mol were then optimized with Gaussian [43], at the M06/6-31g(d) level including D3-dispersion corrections [44], followed by a second stage of duplicates filtering, as described above. This process resulted with 159 conformers within an energy range of 20.4 kcal/mol. These define our pillar[5]arene data set.

Symmetry analysis of pillar[5]arenes

For each conformer of the pillar[5]arene set we calculated the distortion with respect to C5, C2 and Cs point groups for the host molecule without the Li+ ion and excluding the hydrogen atoms. Table 1 summarizes results based on the exact CSM algorithm, teaching that none of the conformers in our dataset had perfect C5 symmetry. Two conformers had perfect C2 symmetry. We note that S(Cs) is a measure of chirality in this case since the distance of the pillar[5]arene structure from higher order achiral structures (belonging to e.g., Ci, S4, S6,… point groups) is larger. As seen in Table 1, all the conformers are chiral to some degree. Although our focus is benchmark of symmetry algorithms, a note about the energy is in place. Generally, direct correlation between energy and symmetry or chirality was not detected for this set. The complex with the minimal energy is not the most symmetric. This is not surprising given the size of the Li+ ion compared to the pillar[5]arene cavity and the expected prevalence of host–guest electronic interactions over the symmetry of the host.

Table 1 Descriptive statistics for 159 conformers of pillar[5]arene

Figure 3 presents correlation plots between the exact algorithm and two approximate algorithms: The Hungarian algorithm and the approximate structure-preservation algorithm for S(C5). Figure 4 presents similar correlations for S(C2) with four algorithms (Hungarian, approximate structure preservation, greedy and Fibonacci lattice with 100 directions coupled with the Hungarian algorithm). Additional figures are provided in Additional file 1 : Results for the other approximate approaches for S(C5) which were similar to the Hungarian algorithm (Additional file 1: Fig. S1), and results for S(Cs) which were similar to S(C2) (Additional file 1: Fig. S2). It is striking to see how the approximate structure-preserving algorithm outperforms all other algorithms, with perfect correlation for S(C5) and very good correlations for S(C2) and S(Cs). For a D5-symmetric pillar[5]arene, there is a single C5 rotation axis, and five different C2 rotation axes. Naturally, when the number of possible permutations is small, it is more likely that an algorithm that is forced to preserve the structure will converge to the correct direction of the rotation axis in space. When the algorithm is not forced to preserve the structure, the calculation can converge with a different permutation, which can either be related to a direction of true rotation axis in the molecule with higher CSM, or a direction that is related to a permutation that does not preserve the structure. Therefore, even a Fibonacci lattice with 100 directions in space provides only an approximate value for the CSM. The Fibonacci method can potentially be improved by increasing the number of directions on the expanse of the computational effort. However, calculations with 500 directions for pillar[5]arene had minor effect on S(C5) and negligible effect on S(C2) and S(Cs). When more directions are added, the calculation becomes much slower, and is not justified here. Nevertheless, as discussed above, the method does provide a lower bound to the CSM for C2 and Cs symmetries. Comparing the correlations of the Hungarian and greedy algorithms, both of them are less accurate as compared with the other algorithms, with the first showing somewhat better correlation with the exact algorithm.

Fig. 3
figure 3

Approximate versus exact values of S(C5) for the pillar[5]arene dataset. A Hungarian algorithm B Approximate structure preservation algorithm. Magenta line represents the y = x curve

Fig. 4
figure 4

Approximate versus exact values of S(C2) for the pillar[5]arene dataset. A Hungarian algorithm B Approximate structure preservation algorithm. C Greedy algorithm. D Fibonacci lattice with 100 directions. Magenta line represents the y = x curve

Another striking finding from Fig. 4 and Additional file 1: Fig. S2 is the high agreement between all the algorithms when the distortion level is low. As seen in Fig. 3B and Additional file 1: Fig. S1, when the CSM is smaller than ~ 2, all the approximate algorithms reach the same CSM as the exact algorithm. A value of 2 for pillar[5]arene means that the distortion is not very high. As an example, Additional file 1: Fig. S3 shows a distorted structure for which S(C2) = 1.8918 superimposed on its nearest C2-symmetric structure. For other molecules the specific threshold of 2 may vary, but the principle remains: when the distortion is relatively low, all algorithms reach the same CSM. This finding provides us with a very basic guideline for interpreting the values of approximate CSMs: as structural descriptors, they can distinguish between conformers regardless of their symmetry, but their interpretation as structure preserving symmetry measures is most likely guaranteed when the level of distortion is low. When the distortion is high, the nearest symmetric structure may not preserve the chemical essence of the structure in terms of its chemical bonds, yet still maintains its symmetric shape.

Another measure of accuracy of the approximate CSMs is the level of structure preservation with respect to the connectivity map of the molecules [12], displayed in Table 2. The approximate structure preservation algorithm reaches 100% by definition, while the other algorithms preserve 63–90% of the structure, depending on the algorithm and point group. Comparing Tables 1 and 2 we note that structure preservation is better when the CSM is small. Thus, structure preservation with respect to the C2 and Cs point groups is better than with respect to C5. In addition, as one may expect, the structure preservation with the Hungarian algorithm is somewhat better as compared with the greedy algorithm, and further improved upon using the Fibonacci lattice approach, particularly for C2 and Cs, where more than one direction is possible. Finally we recall that true Cs symmetry is unlikely for pillar[5]arene and it is used here due to its interpretation as a continuous measure of chirality.

Table 2 Structure preservation during approximate CSM calculations for 159 conformers of pillar[5]arene

Last but not least, let us discuss the time of calculation. We ran the code using 1 core on our Intel(R) Xeon(R) Gold 6130 CPU@2.10 GHz server. Table 3 presents the total time for calculating the CSM for 20 conformers (one after the other) taken from the pillar[5]arene data set. Most of the calculations were completed within a few seconds and even less than that, with an average time per molecule of 0.05–6.51 s. The size of the molecule makes the exact algorithm faster than all others, since it scans a smaller number of permutations. Among the approximate approaches, the greedy algorithm is the fastest, with the Hungarian and the approximate structure preservation algorithms coming next with negligible differences between them. The Fibonacci lattice algorithm naturally requires longer time for computation, up to ~ 100 times longer when 500 directions are taken into account.

Table 3 User time (in seconds) for consecutive CSM calculation of 20 conformers of Pillar[5]arene

C100 fullerenes

The Fullerene family of molecules is an allotrope form of carbon, characterized with hollow structures. Fullerene molecules generally display high symmetry, with the most abundant fullerene, C60, showing icosahedral symmetry [45]. Chiral fullerenes have also been documented [46, 47]. Fullerenes have numerous applications that exploit their unique symmetry, in host–guest chemistry, solar cells, catalysis, drug design, and cancer treatments [45,46,47,48,49,50]. Understanding the symmetry and chirality, of fullerenes, particularly when it is only approximate, can shed light on different distortive processes. Here we focus on fullerenes with 100 carbon atoms. Theoretical studies show that this system have 450 topological isomers that follow the isolated-pentagon-rule [51, 52]. C100 ions were recently detected experimentally as original constituent of aerosol samples [53].

While applicable for an exact CSM calculation, C100 fullerene provides an algorithmic challenge for approximate CSM estimation due to variability in topology and the large number of permutations. Coordinates of the 450 isomers of C100 fullerene were downloaded from the fullerene library [52] without modifications. These structures are based on the Yoshida’s fullerene library and were further optimized by Tománek using the fast Dreiding-like force field [52, 54]. The data set is divided to 336 topologically asymmetric isomers (that belong to the C1 point group), 62 isomers with C2 symmetry, and 31 isomers with Cs symmetry. The rest 21 isomers display higher symmetry. Examples of three isomers with different symmetries are presented in Fig. 5. Our first goal was to test the set of algorithms on symmetric structures in order to see whether the approximate approaches can identify the correct symmetry. All methods were able to correctly estimate S(C2) for the C2-symmetric isomers and S(Cs) for the Cs-symmetric isomers. The CSM values varied between 0.0000 and 0.0004. The small non-zero CSM values result from negligible numerical inaccuracies of the original coordinates which are reasonable to ignore in this case. No significant differences were found between the different algorithms for these molecules.

Fig. 5
figure 5

Examples of three isomers of C100 from the Yoshida fullerene library that belong to different point groups. A Isomer 134 with C1 symmetry; B Isomer 221 with C3 symmetry shown with blue line; C Isomer 126 with Cs symmetry shown with a red reflection plane

The set of 336 asymmetric fullerenes was more challenging. For this set, we calculated the CCM as the minimum of S(Cs) and S(Ci). The CSM with respect to the S4 point group was also calculated but found higher than S(Cs) for all isomers and will not be discussed here. Since the asymmetry stems from the connectivity map of the molecule and not from conformational distortion, the nearest symmetric or achiral structure, regardless of the algorithm or the point group at hand, can not be a real molecule and atoms' overlap is to be expected for this structure. Nevertheless, the CSM or CCM still have a mathematical meaning of a symmetry or chirality measure, since the reference structure found by the algorithm is the closest structure that belongs to the desired point group, as originally defined by Avnir and coworkers [7, 9, 10].

To elaborate on this point, let us look at the distortion with respect to inversion symmetry of the asymmetric C100 isomers. The exact algorithm returns the value of 100 for S(Ci) for each one of the 336 isomers, since the nearest structure with inversion symmetry, that maintains the connectivity map of the molecule, collapses to the center of mass, and all atoms overlap each other. In this sense the exact CSM cannot distinguish between the different isomers. The Hungarian and greedy algorithms on the other hand, find Ci-symmetric structures with lower percentages of atoms' overlap, that retain the hollow shape of the fullerene. Figure 6 displays the nearest Ci-symmetric structures of the C1 isomer presented in Fig. 5A, calculated with the Hungarian and greedy algorithms. The resulting structures differ from each other and so are the values of S(Ci). The ability to distinguish between the original isomers makes the approximate CSMs better structural descriptors for this set of isomers. Therefore, relaxing the requirement of structure preservation is advantageous in this case.

Fig. 6
figure 6

Nearest structure with inversion symmetry for the C1-isomer #134 of the C100 fullerene set calculated by two approximate algorithms. A Hungarian algorithm, S(Ci) = 0.7296. B Greedy algorithm, S(Ci) = 1.2325. Original structure is presented in Fig. 5A

Which of the approximate CSM algorithms provides a more accurate answer to S(Ci)? Our main goal is to find the minimal CSM, but we also attempt to reach a reference structure that maintains, as much as possible, the chemical essence of the original molecule. For isomers of C100 with C1 symmetry, calculation of S(Ci) with the Hungarian algorithm outperforms the greedy algorithm in both criteria: It finds smaller CSM values (Table 4) and reaches higher percentages of structure preservation as discussed below. The Fibonacci lattice algorithm is unnecessary here since the inversion point must be at the center of mass of the molecule and there is no need to test further directions in space. The approximate structure preserving algorithm reaches the value of 100 as the exact algorithm, with 100% of structure preservation, as expected.

Table 4 Descriptive statistics of CSM for 336 asymmetric isomers of C100

Repeating this analysis for S(Cs) we found smaller CSM values as compared with S(Ci), for all algorithms and all the isomers. We can thus claim that CCM = S(Cs) in this case, that is, the chirality of the C100 isomers stems from the lack of reflection symmetry. Figure 7 presents a box and whisker plot of the CCM calculated with the different algorithms. Using the exact algorithm (as well as the approximate structure preserving algorithm) the nearest structure with reflection symmetry collapses to a planar surface, with a mean CSM of 27. In other words, the exact CSM can be interpreted as a planarity measure in this case. With the approximate algorithms, the nearest symmetric structures maintain the ellipsoid shape by overlapping fewer atoms, leading to much smaller CSM values. We note that unlike the exact algorithm, the approximate algorithms predict that different isomers will have different chirality levels. Additional file 1: Fig. S4 displays the nearest symmetric structure with reflection symmetry that was obtained for isomer #134 presented in Fig. 5A. In this sense the approximate CSM is a good 3D-shape descriptor. The Fibonacci lattice algorithm improves the Hungarian algorithm and finds reference structures with smaller CSM values. Increasing the number of directions to 500 improves the results even further, but with minor effect, therefore not shown here. Figure 8 shows that as compared with S(Ci), structure preservation for reflection symmetry is better with all algorithms, and further improves with the Fibonacci lattice algorithm. Table 4 summarizes descriptive statistics of these calculations and Table 5 presents data on the time of calculation for a subset of 20 C100 fullerenes on our server. Unlike the pillar[5]arene data set, here the exact algorithm was slower than some of the approximate algorithms. Particularly for S(Ci), an approximate algorithm that does not attempt to preserve the structure is much faster than an exact calculation. The greedy algorithm was the fastest algorithm for both S(Ci) and S(Cs).

Fig. 7
figure 7

Box and whisker plot of CCM values for 336 asymmetric C100 isomers, calculated with different methods. Box boundaries represent 25–75% of the data, horizontal line within the box is the median and the white star is the mean value in each box

Fig. 8
figure 8

Box and whisker plot of the percentages of structure preservation for approximate CSM calculations for 336 asymmetric C100 isomers. Left: S(Ci). Right: S(Cs). Box boundaries represent 25–75% of the data, horizontal line within the box is the median and the white star is the mean value in each box

Table 5 User time (in seconds) for consecutive CSM calculation of 20 isomers of C100 fullerenes

Metal organic frameworks

As a third challenge we focused on crystal structures for which the exact algorithm cannot scan the permutation space in a reasonable time frame. Our purpose was to test whether perfect symmetry can be identified by approximate CSM algorithms, to estimate the time of calculation, and to evaluate the extent by which the structures of the initial sets of coordinates are preserved in the final sets. Coordinates of three crystals of MOFs with very large unit cells were retrieved from the Crystallographic Open Database (COD) [55]. Structures that contain all the molecules whose centroids fit inside a single unit cell were constructed using Mercury [56], and are presented in Fig. 9. Continuous symmetry measures with respect to several point groups were calculated with the Hungarian algorithm and are summarized in Table 6. Perfect symmetry, as expected from the space group, is clearly identified with 100% structure preservation. Our method correctly predicts that the unit cell of UHM-25-Ala-Boc [57] (Fig. 9A) that belongs to the P432 space group is chiral, while the other two structures are achiral. We note that the non-zero CSM values of UHM-25-Ala-Boc are relatively small. This stems from the large number of atoms that significantly increases the normalization factor, \(D\) in (1). Results for the greedy algorithm were similar to the Hungarian algorithm, with a bit higher CSM in the lack of perfect symmetry in most of the cases (Additional file 1: Table S1), on the expense of lower structure preservation (Additional file 1: Table S2). With regards to time of calculation, the greedy algorithm was significantly faster, and ranged between ca. 4–35 s per structure. The Hungarian algorithm required longer time in most cases, ranging between ca. 4–221 s per structure (Additional file 1: Table S3). The Fibonacci lattice algorithm required much longer calculation time and is not justified here.

Fig. 9
figure 9

A set of three MOFs, viewed along their unit cell c axis, without the hydrogens atoms taken from the COD IDs: A 4002650, B 4003149, C 4002646

Table 6 Approximate CSMs of unit cells of MOFs, calculated with the Hungarian algorithm

Guidelines for effective usage

The CSM code offers several algorithms and many parameters that affect the accuracy and speed of the calculation, as well as the chemical and mathematical interpretation of the results. Choosing between all possible options is not always trivial, and may benefit from a trial and error strategy. Nevertheless, based on our current findings and previous studies [12, 13], we provide here several guidelines for the choice of algorithms for CSM and CCM analysis of molecular structures.

Small-to-medium sized molecules are suitable for an exact calculation, based on the structure preserving permutation algorithm [12]. The number of atoms is the first indicator of size. A molecule with up to several hundreds of atoms is generally considered medium sized here, but the bonding structure and the size of the symmetry equivalent atoms' groups is also important, as they affect the number of possible permutations. Ignoring the hydrogen atoms is a common workaround that reduces the number of permutations and can speed up the calculation. While the hydrogen atoms naturally affect the CSM value of a given molecule, they often have minor effect on the overall distortion trends of a set of related molecules. Another option for size reduction is to analyze the core structure of large molecules, particularly when this structure has a symmetric topology. Sets of molecular derivatives with a common skeleton, can especially benefit from this type of analysis [15].

Large structures for which the exact algorithm is not practical, are handled with approximate CSM calculations. Several algorithms can be used, and are selected as an interplay between calculation time and accuracy:

  1. 1.

    The Hungarian algorithm is our default choice for approximate CSM calculation, presenting a reasonable compromise between speed and accuracy. This algorithm also excels in distinguishing between different molecules, and as such can be used with medium-sized molecules as well, like the fullerenes discussed above.

  2. 2.

    The approximate structure-preserving algorithm is more accurate than the Hungarian algorithm in its ability to find permutations that preserve the chemical essence of the structure. The resulting CSM values correlate well with the corresponding values of the exact algorithm, as seen for the sets of pillar[5]arene complexes and the fullerenes presented here. The drawback is that these calculations take longer time.

  3. 3.

    The greedy algorithm is often the fastest choice among all approximate algorithms. It is a relatively crude approximation, showing the lowest accuracy in terms of structure preservation. Nevertheless, like the Hungarian algorithm it excels in terms of distinguishing between molecules.

  4. 4.

    The Fibonacci sphere algorithm with up to 100 directions can improve the accuracy of the Hungarian algorithm and find permutations that lead to smaller CSM values. The drawback is a longer calculation time. It is particularly useful for molecules with large equivalence groups of atoms, (e.g., highly symmetric nanoparticles) for which a slight change of the initial direction can help the code to find a better permutation. Increasing the number of directions beyond 100 can increase the accuracy along with the time of calculation.

  5. 5.

    In the special case of protein oligomers, one first need to make sure that chains' length and sequence are equalized. Then the Hungarian algorithm is recommended together with special features of sequence preservation discussed in our previous publication [13].

Calculation of the CCM follows a similar strategy as the CSM with one additional guideline. In many cases of chiral or asymmetric topology, the chirality stems from the lack of reflection symmetry. That is, S(Cs) is expected to be smaller than S(Ci) or S(Sn) with n > 2. Therefore, it is often enough to calculate S(Cs) in order to speed up the calculation of the CCM, as was done here for pillar[5]arenes.

The CSM and the alignment problem

The problem of calculating the continuous symmetry measure is somewhat related to the problem of molecular alignment and RMSD (root mean square deviation) calculations, where one is interested in finding the best fit between two structures.

More concretely, given two sets of atom positions—\({{\mathbf{P}}} = \left\{ {{{\mathbf{P}}}_k :\,\,1 \le k \le N} \right\}\) and \({\mathbf{Q}} = \left\{{\mathbf{Q}}_k :\,\,1 \le k \le N \right\}\), alignment algorithms look for a permutation \(\pi\) of \(\left\{ {1,2,\ldots,N} \right\}\) and a spatial transformation \(T\) such that the RMSD is minimized:

$$RMSD = \sqrt {{\frac{1}{N}\sum_{k = 1}^N {\left| {TP_k - Q_{\pi (k)} } \right|^2 } }}$$

This problem has been widely discussed in the literature [60,61,62].

If we take \({{\mathbf{P}}} = {{\mathbf{Q}}}\), then the alignment problem becomes very similar to the problem of finding the CSM. However, when we look for symmetry, we have to exclude the trivial solution where \(T\) is the identity transformation. We also impose restrictions on the operation \(T\). In the context of point symmetry, \(T\) must generate a finite group. This means, for example, that if \(T\) is a rotation around an axis vector, then the angle of rotation must be a multiple of \(360^{\circ} /n\) where \(n\) is the order of the symmetry operation. In contrast, when we look for mere alignment, \(T\) can be a rotation at any angle.

The continuous shape measure, CShM, originally developed by Pinsky and Avnir [11] as an extension of the CSM method, follows a similar strategy. It becomes a symmetry measure when the reference shape is symmetric. Various applications of this approach for small molecules were published through the years for the analysis of e.g., coordination complexes and their distortion pathways, particularly with respect to platonic solids [30, 31, 63]. The techniques described here could be applied, with some adjustments, to calculate shape measures as well as solving the alignment problems for large molecules. Progress in this direction has already been documented, in the form of alignment algorithms that calculate the RMSD while taking symmetry into considerations, utilizing, for example, the Hungarian algorithm [64, 65] and the greedy algorithm with some level of structure preservation [66].


The continuous symmetry and chirality measures estimate the distance between a given molecule and its nearest symmetric (or achiral) structure. Preserving the connectivity map of the original molecule, although desired, was not a requirement in the original CSM algorithm [7, 10]. Consequently, the approach was limited by the size of molecules that could be handled, an obstacle that was partially overcame by applying an approximate approach based on a permutation-direction search using the greedy algorithm [9]. In recent years, our group considerably improved the method. We introduced the concept of structure preserving permutations in order to reduce the number of permutations the code needs to scan, while preserving the connectivity map of the original molecule [12]. We further improved the approximate algorithm for protein homomers by applying the Hungarian algorithm instead of the greedy algorithm in conjunction with utilizing the sequence structure to reduce the number of permutations [13]. In this work, we took the method one step forward and developed a set of approximate algorithms that are suitable to any type of molecular structure. Our main goal is to provide a quantitative estimation to the symmetry and chirality content of molecules of any size. Our set of approximate CSM algorithms should thus be viewed as a collection of 3D geometrical descriptors for global structural analysis that fulfills the main criteria for good descriptors, in terms of interpretability, ability to differentiate between isomers and conformers, applicability for local structures, continuity, usage simplicity and efficiency [67]. As such, these descriptors are suitable for characterizing structural changes and play part in QSAR/QSPR modeling.

The accuracy of all the approximate methods described in this paper depends on the distortion level of the molecule with respect to the desired symmetry point group G. For molecules which are only slightly distorted, all the algorithms tend to find a permutation that preserve the structure (in terms of its connectivity map) either fully or with a high preservation rate. For highly distorted molecules, this does not always hold, as permutations which do not preserve the structure could attain lower values of the target function, thereby creating a bias between the exact and approximate CSMs.

Based on our calculations of three different sets of molecules (Pillar[5]arenes complexes, C100 fullerenes, and MOFs), guidelines for effective usage of the different algorithms were specified. These are based on the ability of the code to find structure preserving permutations (wherever relevant), the speed of the calculation and the ability of the CSM to distinguish between the studied molecules. In summary, we emphasize the ease and efficiency of using the CSM and CCM approach for structural analysis, which make them applicable as robust three-dimensional geometrical descriptors that can be used to follow dynamical processes and statistical studies on quantitative structure–activity and structure-properties relationships.

Availability of data and materials

The datasets supporting the conclusions of this article are available at: Free online calculators of the CSM and CCM for small molecules and protein homomers are available at The CSM software is freely available. Project name: CSM; Project homepage:; Archived version: 1.3.7.b1; Operating system(s): Linux, Windows; Programming language: Python, c++; Other requirements: OpenBabel 3.1.1, Cython 0.29.14, scipy 1.7.3, c++ compiler, conda, numpy; License: GNU-GPL version 2; Any restrictions to use by non-academics: Not applicable.


  1. The concise presentation of the Fibonacci lattice in Eqs. (5) and (6) is based on a code that was published under the Stack Overflow platform by the user Fnord. After testing this code and confirming its reliability, we applied it in the CSM software. The original code is available at: Accessed 2023, July 28.



Continuous chirality measure


Crystallographic open database


Continuous shape measure


Continuous symmetry measure


Metal organic framework


Quantitative structure activity relationship


Quantitative structure property relationship


Root mean square deviation


  1. Hargittai M, Hargittai I (2009) Symmetry through the Eyes of a Chemist, 3rd edn. Springer, Dordrecht

    Book  Google Scholar 

  2. Goodsell DS, Olson AJ (2000) Structural symmetry and protein function. Annu Rev Biophys Biomol Struct 29:105–153.

    Article  CAS  PubMed  Google Scholar 

  3. Alemany P, Casanova D, Alvarez S, Dryzun C, Avnir D (2017) Continuous symmetry measures: a new tool in quantum chemistry. In: Parrill AL, Lipkowitz KB (eds) Reviews in computational chemistry. Wiley, London, pp 289–352

    Chapter  Google Scholar 

  4. Zwart PH, Grosse-Kunstleve RW, Lebedev AA, Murshudov GN, Adams PD (2008) Surprises and pitfalls arising from (pseudo)symmetry. ACTA Crystallogr Sect D STRUCTURAL Biol 64:99–107.

    Article  CAS  Google Scholar 

  5. Pratt Brock C (2023) Prevalence and significance of approximate symmetry in organic Pc Structures. Helv Chim Acta.

    Article  Google Scholar 

  6. Nespolo M, Benahsene AH (2021) Symmetry and chirality in crystals. J Appl Crystallogr 54:1594–1599.

    Article  CAS  Google Scholar 

  7. Zabrodsky H, Peleg S, Avnir D (1992) Continuous symmetry measures. J Am Chem Soc 114:7843–7851.

    Article  CAS  Google Scholar 

  8. Pinsky M, Dryzun C, Casanova D, Alemany P, Avnir D (2008) Analytical methods for calculating continuous symmetry measures and the chirality measure. J Comput Chem 29:2712–2721.

    Article  CAS  PubMed  Google Scholar 

  9. Dryzun C, Zait A, Avnir D (2011) Quantitative symmetry and chirality—a fast computational algorithm for large structures: proteins, macromolecules, nanotubes, and unit cells. J Comput Chem 32:2526–2538.

    Article  CAS  PubMed  Google Scholar 

  10. Zabrodsky H, Avnir D (1995) Continuous symmetry measures. 4. Chirality. J Am Chem Soc 117:462–473.

    Article  CAS  Google Scholar 

  11. Pinsky M, Avnir D (1998) Continuous symmetry measures. 5. The classical Polyhedra. Inorg Chem 37:5575–5582.

    Article  CAS  PubMed  Google Scholar 

  12. Alon G, Tuvi-Arad I (2018) Improved algorithms for symmetry analysis: structure preserving permutations. J Math Chem 56:193–212.

    Article  CAS  Google Scholar 

  13. Tuvi-Arad I, Alon G (2019) Improved algorithms for quantifying the near symmetry of proteins: complete side chains analysis. J Cheminform 11:39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Alvarez S, Alemany P, Avnir D (2005) Continuous chirality measures in transition metal chemistry. Chem Soc Rev 34:313–326.

    Article  CAS  PubMed  Google Scholar 

  15. Kaspi-Kaneti AWWA, Tuvi-Arad I (2018) Twisted and bent out of shape: symmetry and chirality analysis of substituted ferrocenes. Organometallics 37:3314–3321.

    Article  CAS  Google Scholar 

  16. Tuvi-Arad I, Avnir D (2010) Determining symmetry changes during a chemical reaction: the case of diazene isomerization. J Math Chem 47:1274–1286.

    Article  CAS  Google Scholar 

  17. Tuvi-Arad I, Avnir D (2012) Symmetry-enthalpy correlations in diels-alder reactions. Chem A Eur J 18:10014–10020.

    Article  CAS  Google Scholar 

  18. Kaspi-Kaneti AWAW, Barroso J, Merino G, Avnir D, Garzón IL, Tuvi-Arad I (2020) Head to tail distortion wave characterizes the enantiomerization of helicenes. J Org Chem 85:15415–15421.

    Article  CAS  PubMed  Google Scholar 

  19. Tuvi-Arad I, Stirling A (2016) The distortive nature of temperature—a symmetry analysis. Isr J Chem 56:1067–1075.

    Article  CAS  Google Scholar 

  20. Fossépré M, Tuvi-Arad I, Beljonne D, Richeter S, Clément S, Surin M (2020) Binding mode multiplicity and multiscale chirality in the supramolecular assembly of DNA and a π-conjugated polymer. ChemPhysChem.

    Article  PubMed  Google Scholar 

  21. Tuvi-Arad I, Rozgonyi T, Stirling A (2013) Effect of temperature and substitution on cope rearrangement: a symmetry perspective. J Phys Chem A 117:12726–12733.

    Article  CAS  PubMed  Google Scholar 

  22. Bonjack-Shterengartz M, Avnir D (2015) The near-symmetry of proteins. Proteins Struct Funct Bioinforma 83:722–734.

    Article  CAS  Google Scholar 

  23. Keinan S, Avnir D (2000) Quantitative symmetry in structure-activity correlations: the near C-2 symmetry of inhibitor/HIV protease complexes. J Am Chem Soc 122:4378–4384

    Article  CAS  Google Scholar 

  24. Shalit Y, Tuvi-Arad I (2020) Side chain flexibility and the symmetry of protein homodimers. PLoS ONE 15:e0235863.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Crippen GM (2008) Chirality descriptors in QSAR. Curr Comput Aided Drug Des 4:259–264

    Article  CAS  Google Scholar 

  26. Jamróz MH, Rode JE, Ostrowski S, Lipiński PFJ, Dobrowolski JC (2012) Chirality measures of α-amino acids. J Chem Inf Model 52:1462–1479.

    Article  CAS  PubMed  Google Scholar 

  27. Milner D, Raz S, Hel-Or H, Keren D, Nevo E (2007) A new measure of symmetry and its application to classification of bifurcating structures. Pattern Recognit 40:2237–2250.

    Article  Google Scholar 

  28. Saragusti I, Sharon I, Katzenelson O, Avnir D (1998) Quantitative analysis of the symmetry of artefacts: Lower paleolithic handaxes. J Archaeol Sci 25:817–825.

    Article  Google Scholar 

  29. Iovita R, Tuvi-Arad I, Moncel M-HM-HMH, Despriee J, Voinchet P, Bahain JJJ-JJ-J, Despriée J, Voinchet P, Bahain JJJ-JJ-J, Desprieae J, Voinchet P, Bahain JJJ-JJ-J (2017) High handaxe symmetry at the beginning of the European Acheulian: the data from la Noira (France) in context. PLoS ONE 12:e0177063.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Keinan S, Avnir D (2001) Continuous symmetry analysis of tetrahedral/planar distortions. Copper chlorides and other AB(4) species. Inorg Chem 40:318–323.

    Article  CAS  PubMed  Google Scholar 

  31. Casanova D, Cirera J, Llunell M, Alemany P, Avnir D, Alvarez S (2004) Minimal distortion pathways in polyhedral rearrangements. J Am Chem Soc 126:1755–1763.

    Article  CAS  PubMed  Google Scholar 

  32. Munkres J (1957) Algorithms for the assignment and transportation problems. J Soc Ind Appl Math 5:32–38.

    Article  Google Scholar 

  33. González Á (2009) Measurement of areas on a sphere using Fibonacci and latitude–longitude lattices. Math Geosci 42:49.

    Article  Google Scholar 

  34. Kildgaard JV, Mikkelsen KV, Bilde M, Elm J (2018) Hydration of atmospheric molecular clusters: a new method for systematic configurational sampling. J Phys Chem A 122:5026–5036.

    Article  CAS  PubMed  Google Scholar 

  35. Ogoshi T, Kanai S, Fujinami S, Yamagishi TA, Nakamoto Y (2008) Para-bridged symmetrical pillar[5]arenes: their Lewis acid catalyzed synthesis and host-guest property. J Am Chem Soc 130:5022–5023.

    Article  CAS  PubMed  Google Scholar 

  36. Ogoshi T, Yamagishi TA, Nakamoto Y (2016) Pillar-shaped macrocyclic hosts pillar[n]arenes: new key players for supramolecular chemistry. Chem Rev 116:7937–8002.

    Article  CAS  PubMed  Google Scholar 

  37. Wu X, Gao L, Hu XY, Wang LY (2016) Supramolecular drug delivery systems based on water-soluble pillar[n]arenes. Chem Rec 16:1216–1227.

    Article  CAS  PubMed  Google Scholar 

  38. Chen L, Cai Y, Feng W, Yuan L (2019) Pillararenes as macrocyclic hosts: a rising star in metal ion separation. Chem Commun 55:7883–7898.

    Article  CAS  Google Scholar 

  39. Zhang HC, Liu ZN, Xin FF, Zhao YL (2020) Metal-ligated pillararene materials: from chemosensors to multidimensional self-assembled architectures. Coord Chem Rev.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Shalit Y, Tuvi-Arad I (2021) Symmetry-binding correlations of crown ethers complexes with Li+ and Na+. ACS Omega 6:1–5.

    Article  CAS  Google Scholar 

  41. Labute P (2010) LowModeMD—implicit low-mode velocity filtering applied to conformational search of macrocycles and protein loops. J Chem Inf Model 50:792–800.

    Article  CAS  PubMed  Google Scholar 

  42. Molecular Operating Environment (MOE) (2019). 2019.12.

  43. Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, Scalmani G, Barone V, Petersson GA, Nakatsuji H, Li X, Caricato M, Marenich AV, Bloino J, Janesko BG, Gomperts R, Mennucci B, Hratchian HP, Ortiz JV, Izmaylov AF, Sonnenberg JL, Williams, Ding F, Lipparini F, Egidi F, Goings J, Peng B, Petrone A, Henderson T, Ranasinghe D, Zakrzewski VG, Gao J, Rega N, Zheng G, Liang W, Hada M, Ehara M, Toyota K, Fukuda R, Hasegawa J, Ishida M, Nakajima T, Honda Y, Kitao O, Nakai H, Vreven T, Throssell K, Montgomery JA Jr, Peralta JE, Ogliaro F, Bearpark MJ, Heyd JJ, Brothers EN, Kudin KN, Staroverov VN, Keith TA, Kobayashi R, Normand J, Raghavachari K, Rendell AP, Burant JC, Iyengar SS, Tomasi J, Cossi M, Millam JM, Klene M, Adamo C, Cammi R, Ochterski JW, Martin RL, Morokuma K, Farkas O, Foresman JB, Fox DJ (2016) Gaussian 16 Rev. C.01

  44. Grimme S, Antony J, Ehrlich S, Krieg H (2010) A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J Chem Phys 132:154104.

    Article  CAS  PubMed  Google Scholar 

  45. Glukhova OE, Druzhinin AA, Zhbanov AI, Rezkov AG (2005) Structure of high-symmetry fullerenes. J Struct Chem 46:501–507.

    Article  CAS  Google Scholar 

  46. Thilgen C, Diederich F (2006) Structural aspects of fullerene chemistry—a journey through fullerene chirality. Chem Rev 106:5049–5135.

    Article  CAS  PubMed  Google Scholar 

  47. Fernández-García JM, Evans PJ, Filippone S, Herranz MÁ, Martín N (2019) Chiral molecular carbon nanostructures. Acc Chem Res 52:1565–1574.

    Article  CAS  PubMed  Google Scholar 

  48. Popov AA, Yang S, Dunsch L (2013) Endohedral fullerenes. Chem Rev 113:5989–6113.

    Article  CAS  PubMed  Google Scholar 

  49. Fernandes NB, Shenoy RUK, Kajampady MK, DCruz CEM, Shirodkar RK, Kumar L, Verma R (2022) Fullerenes for the treatment of cancer: an emerging tool. Environ Sci Pollut Res 29:58607–58627.

    Article  CAS  Google Scholar 

  50. Ganesamoorthy R, Sathiyan G, Sakthivel P (2017) Review: fullerene based acceptors for efficient bulk heterojunction organic solar cell applications. Sol Energy Mater Sol Cells 161:102–148.

    Article  CAS  Google Scholar 

  51. Zhao X, Goto H, Slanina Z (2004) C100 IPR fullerenes: temperature-dependent relative stabilities based on the Gibbs function. Chem Phys 306:93–104.

    Article  CAS  Google Scholar 

  52. Tománek D, Frederick N (2013) Cn fullerenes. Accessed 21 Feb 2023

  53. dos Santos FN, Nascimento MM, da Rocha GO, de Andrade JB (2023) The occurrence of pristine and functionalized fullerenes as constituents of airborne aerosols. Sci Rep 13:1–9.

    Article  CAS  Google Scholar 

  54. Tománek D (2014) Guide through the nanocarbon jungle. Morgan & Claypool Publishers, San Rafael

    Book  Google Scholar 

  55. Gražulis S, Daškevič A, Merkys A, Chateigner D, Lutterotti L, Quirós M, Serebryanaya NR, Moeck P, Downs RT, Le Bail A (2012) Crystallography open database (COD): an open-access collection of crystal structures and platform for world-wide collaboration. Nucl Acids Res 40:420–427.

    Article  CAS  Google Scholar 

  56. Macrae CF, Sovago I, Cottrell SJ, Galek PTA, McCabe P, Pidcock E, Platings M, Shields GP, Stevens JS, Towler M, Wood PA (2020) Mercury 4.0: from visualization to analysis, design and prediction. J Appl Crystallogr 53:226–235.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Sartor M, Stein T, Hoffmann F, Fröba M (2016) A new set of isoreticular, homochiral metal–organic frameworks with UCP topology. Chem Mater 28:519–528.

    Article  CAS  Google Scholar 

  58. Gao W-Y, Cai R, Pham T, Forrest KA, Hogan A, Nugent P, Williams K, Wojtas L, Luebke R, Weseliński ŁJ, Zaworotko MJ, Space B, Chen Y-S, Eddaoudi M, Shi X, Ma S (2015) Remote stabilization of copper paddlewheel based molecular building blocks in metal–organic frameworks. Chem Mater 27:2144–2151.

    Article  CAS  Google Scholar 

  59. Lee SJ, Doussot C, Baux A, Liu L, Jameson GB, Richardson C, Pak JJ, Trousselet F, Coudert F-X, Telfer SG (2016) Multicomponent metal–organic frameworks as defect-tolerant materials. Chem Mater 28:368–375.

    Article  CAS  Google Scholar 

  60. Kabsch W (1978) A discussion of the solution for the best rotation to relate two sets of vectors. Acta Crystallogr Sect A 34:827–828.

    Article  Google Scholar 

  61. Labute P, Williams C, Feher M, Sourial E, Schmidt JM (2001) Flexible alignment of small molecules. J Med Chem 44:1483–1490.

    Article  CAS  PubMed  Google Scholar 

  62. Zhang Y, Skolnick J (2005) TM-align: a protein structure alignment algorithm based on the TM-score. Nucleic Acids Res 33:2302–2309.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Alvarez S (2005) Polyhedra in (inorganic) chemistry. Dalt Trans 13:2209–2233.

    Article  CAS  Google Scholar 

  64. Allen WJ, Rizzo RC (2014) Implementation of the hungarian algorithm to account for ligand symmetry and similarity in structure-based design. J Chem Inf Model 54:518–529.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Temelso B, Mabey JM, Kubota T, Appiah-Padi N, Shields GC (2017) ArbAlign: a tool for optimal alignment of arbitrarily ordered isomers using the Kuhn–Munkres algorithm. J Chem Inf Model 57:1045–1054.

    Article  CAS  PubMed  Google Scholar 

  66. Bell EW, Zhang Y (2019) DockRMSD: an open-source tool for atom mapping and RMSD calculation of symmetric molecules through graph isomorphism. J Cheminform.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Randić M (1991) Generalized molecular descriptors. J Math Chem 7:155–168.

    Article  Google Scholar 

Download references


We are sincerely grateful for fruitful discussions with Prof. David Avnir (The Hebrew University of Jerusalem). The programming of the new code was done by the Research Software Company, Israel. We are thankful to Mr. Sagiv Barhoom (The Open University) for his help in programming and Ms. Yaffa Shalit (The Open University) for her help in testing the code.


Supported by the Israel Science Foundation (Grant No. 411/15), and the Open University Research Fund (Grants Nos. 102128, 102558, 511711).

Author information

Authors and Affiliations



GA and IT-A mutually contributed to conceptualization and theoretical development. GA developed the mathematical formulation and algorithmic framework, IT-A contributed to the chemical application, validation, formal analysis and overall supervision of the research. YB-H performed the calculations on pillar[5]arenes and the C100 data sets. GA and IT-A wrote the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Gil Alon or Inbal Tuvi-Arad.

Ethics declarations

Competing interests

The authors declare no competing financial interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Additional Tables and Figures that support the research findings.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Alon, G., Ben-Haim, Y. & Tuvi-Arad, I. Continuous symmetry and chirality measures: approximate algorithms for large molecular structures. J Cheminform 15, 106 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: