Fractal and Mathematical Morphology in Intricate Comparison
between Tertiary Protein Structures
Ranjeet Kumar Rout1, Pabitra Pal Choudhury1, B. S. Daya Sagar2, Sk. Sarif Hassan3
1
2

Applied Statistics Unit, Indian Statistical Institute, Kolkata, India,

System Science and Informatics Unit, Indian Statistical Institute, Bangalore, India.

3

Institute of Mathematics and Application, Andharua, Bhubaneswar-751003, India,
Email: ranjeetkumarrout@gmail.com, pabitrapalchoudhury@gmail.com ,
bsdasagar@isibang.ac.in , sarimif@gmail.com

Abstract
Intricate comparison between two given tertiary structures of proteins is as important as the comparison
of their functions. Several algorithms have been devised to compute the similarity and dissimilarity
among protein structures. But, these algorithms compare protein structures by structural alignment of the
protein backbones which are usually unable to determine precise differences. In this paper, an attempt has
been made to compute the similarities and dissimilarities among 3D protein structures using the
fundamental mathematical morphology operations and fractal geometry which can resolve the problem of
real differences. In doing so, two techniques are being used here in determining the superficial structural
(global similarity) and local similarity in atomic level of the protein molecules. This intricate structural
difference would provide insight to Biologists to understand the protein structures and their functions
more precisely.
Keywords: 3D-Protein Structure, Similarities, Mathematical Morphology, Geodesic Dilation, Skeleton,
Fractal Dimension

1. Introduction:
Proteins are made of amino acids chain with its length ranging from 50 to more than 3000. A carbon
atom 𝐶𝛼 is connected to a carboxyl (-COOH) group, an amine (-NH2) group, a hydrogen atom and a
residue (which depends on the specific amino acid) to formulate a single amino acid. The amine group of
an amino acid is covalently bonded by polypeptide bond with the carboxyl group of another amino acid
to form a protein. The sequence of 𝐶𝛼 carbon atoms forms the backbone of the protein. Whenever the
protein is left in its natural environment, it folds to a specific 3D structure. This is due to the forces
between the amino acids such that the total free energy is minimized [1]. This renders a stable 3D protein
structure. Thus, a protein can either be considered as polypeptides sequence of 20 amino acids occurring
1

naturally or as a tertiary structure into which a particular protein folds [2]. The two highly similar amino
acid sequences (primary structures) would not necessarily mean protein functional similarity [2] [3].
Therefore it is evident that comparison of protein functionality requires the intricate comparison of
tertiary structures. The search for an effective solution for tertiary protein structural similarity is justified
because such tools can be of aid to scientists for prediction of the functions of a newly found protein, in
development of procedures for drug design, in the identification of new types of protein architecture, in
the organization of the known database of protein structures by classifying them according to their
structures and can help to discover unexpected evolutionary and functional inter-relations between
proteins [4] [5]. Several algorithms have been devised to compute the similarity between protein
structures but they land up in a difficult computational problem as well as accuracy problem [6]. As, in
many cases there is not even a single superposition that reveals all regions of similarity between
compared proteins (RMSD, DALI, ProSup) [7]. Also, there are many conceptual difficulties associated
with various methods (RMSD, ad hoc scores based on local secondary structure, hydrogen bonding
pattern, burial status, or interaction environment) which have not been resolved [8]. Classical criteria
such as the Root Mean Square Deviation (RMSD) fail to identify similar shapes in a consistent way [9].
To add on various systems have been proposed for structural classification, such as Structural
Classification of Proteins (SCOP), Class Architecture Topology Homology (CATH), Families of
Structurally Similar Proteins (FSSP), and others. The similarity in their cases is computed using
structural alignment algorithms such as DALI, CE, VAST, SSAP and others. Most of these methods are
computationally intensive and time-consuming, especially when searching large databases due to
intrinsic complexity of structural alignment [10]. Also, the prevailing practice in the protein
crystallographic community for computing structural differences is highly inappropriate, in particular
when medium- and low-resolution structures are involved [11]. Geometrical feature like Fractal
dimension of 𝐶𝛼 of the backbone structure of one peptide chain proteins are considered in [12].
Obviously, a more objective method is highly desirable. In this paper these problems have been tried to
resolve by introducing two different methods using Mathematical Morphology and Fractals which would
yield desired output.
The organization of the paper is as follows: in section 2, the basic review of Mathematical
Morphology operations is presented; in section 3, the result and analysis based on the two methods are
discussed; conclusion has been made in section 4.

2. Basics of Mathematical Morphology and Fractal Dimension.
Mathematical Morphology is a widely used paradigm in the field of image processing. Morphological
tools are already very popular for image segmentation, image decomposition etc. Morphological

2

operations like erosion, dilation, opening, closing are used for processing images frequently and produce
results with high accuracy. The definitions of these basic morphological operators are as follows [13].
Erotion: A ⊖ S = a − s: a ∈ A , s ∈ S =

Ms
s∈S

Dilation: A ⊕ S = a + s: a ∈ A , s ∈ S =

Ms
s∈S

Opening ∶ A ∘ S = (A ⊖ S) ⊕ S
Closing ∶ A ⦁ S = A ⊕ S ⊖ S
Where A denotes the shape that is to be transformed and S denotes the structuring element that is used for
the transformation.

2.1. Morphological Skeleton.
Morphological skeleton of every geometrical structure is a subset of the original structure which has the
same connectivity as the original structure from which inference can be drawn. From each point of the
skeleton the distance to the boundary of the original set is the radius of a maximal circle (whose center is
at a point of the skeleton) which touches the boundary at least two different points. The skeleton of an
object gives a clear idea about the shape of the object. For the shape A, and the structuring element S, the
skeleton can be constructed through the operation [14] [15]:
Sk n = (A ⊖ nS)\(A ⊖ nS) ∘ S,

for n = 1, 2, … , N

And the reverse process is as follows, Where, N is the number of performed iterations. Dilating the
skeleton N times iteratively using the multi-scale structuring elements S a shape that is almost same to the
original shape can be achieved.
N
′

A =

Sk n ⊕ nS
n=0

Where nS = S ⊕ S ⊕ … ⊕ S (n times).

2.2. Fractal Dimension.
A fractal dimension is an index for characterizing fractal patterns or sets. The patterns illustrate selfsimilarity and the fractal dimension indicates the extent to which the fractal objects fills a particular
3

Euclidean space in which it is embedded. These dimensions are usually non-integers. The fractal
dimensions can be computed through Box Counting Method which is briefly stated as follows.
Box-Counting Method: This method computes the number of cells required to entirely cover an object,
with grids of cells of varying size. Practically, this is performed by superimposing regular grids over an
object and by counting the number of occupied cells. The logarithm of N(r), the number of occupied cells,
versus the logarithm of

𝟏
𝒓

, where r is the size of one cell, gives a line whose gradient corresponds to the

box dimension [16][17]. To calculate the dimension for a fractal S, the Box-Counting dimension is
defined as,
Dim box(S) = lim𝑛→0

log 𝑁(𝑟)
log

1
r

3. Methods and Results
In this section two different methods are proposed to compute the similarity between tertiary protein
structures in intricate level on the basis of mathematical morphology and fractal dimension.

3.0.1. Tertiary Structure Skeleton and its Fractal Dimension.
The Protein Data Bank (PDB) is the largest and most commonly used repository for any kind of
information regarding proteins. Information like 3D structure, family, function of every protein found till
date is available in PDB. Mainly the X-Ray crystallography and Nuclear Magnetic Resonance is used for
determining the 3D structure of the protein. The 3D structure is represented in (x, y, z) coordinates (with
respect to an arbitrary origin) of the atoms presented in the protein. The '.pdb' files available in the PDB
database contain all the structural information of a protein. Any molecule structure viewer like PyMol,
JMol is able to simulate the 3D protein structure available in the .pdb file. The proteins have an intrinsic
self-similarity as they are hetero-polymers with a variable composition of twenty different amino acids.
Thus, this protein backbone space curve consisting of 𝐶𝛼 atoms motivates us to compare 3D protein
structure on the basis of their fractal features [16]. Following are the steps for calculating the skeletons
and their corresponding fractal dimension.

1. Since the atoms of the tertiary protein structure has three co-ordinates (x, y and z). The
idea of slice representation is to decompose a tertiary structure into a sequence of nonoverlapping pieces, namely slices, by cutting the tertiary structure on its 𝑚 vertices with
𝑚′ planes 𝑝1 , 𝑝2 , 𝑝3 , … , 𝑝𝑚 ′ that are perpendicular to the z-axis and the union of planes
𝑚′ contains all m atoms from the tertiary protein structure. Each slice contains the protein
4

atoms that share the same z-coordinates. Therefore, each slice can be attached with a
unique z-interval.
2. Connected component 𝐶𝑖 is obtained for all slices by the multi-scale opening of the
atoms presented in that slice. 𝐶𝑖 = 𝑃𝑖 ∘ 𝑛𝑆, 𝑤𝑕𝑒𝑟𝑒 𝑖 = 1, 2, … , 𝑚, m is the number
of slices and S is the structuring element.
3. Skeleton of the all connected component
Sk n (𝐶𝑖 ) = (𝐶𝑖 ⊖ nS)\(𝐶𝑖 ⊖ nS) ∘ S for n = 1, 2, … , N and 𝑖 = 1, 2, 3, … , 𝑚.
4. A =

N
n=0 Sk n

⊕ nS

Where, nS = S ⊕ S ⊕ … ⊕ S (n times),

where

A

is

the

connected component of a tertiary protein structure.
5. Compute the fractal Dimension 𝐷𝑝 of the Skeleton A of the corresponding tertiary
protein structure.

Illustration:
Local similarity of the tertiary structure is obtained through the above steps. In Figure 1, a slice is shown
for the protein 2LEP. Each “.” represents an atom. The slice contains all the atoms whose z coordinate is
within -0.600 to -0.700.

After getting a slice, a connected component has been obtained by the multi-

scale opening of the atoms presented in that slice. For each iteration of multi-scale opening the size of
structuring element increased by one. And for the opening, a primitive structuring element of size 𝑛 × 𝑛
where 𝑛 = 1, 2, … 𝑁 is used. By doing this, the tertiary structural atom information has been transformed
into two dimensional planes without losing any information at all. The iterations for the slice shown in
Figure 1 are shown below.

Figure 1 slice for protein 2LEP for z = −0.600 to − 0.700
In Figure 2, a to j shows the multi-scale opening using primitive structuring elements. Starting from size
one, each figure shows the iteration with structuring element larger by ten units from the previous one.
The iterative opening may take a large number of iterations to contain all the atoms in a particular slice.
And there may be more than one plane for a slice. So we dilate the plane with a primitive structuring
element. This reduces the number of planes for each slice. The example of the plane after dilating with a
5

disk shape structuring element of size 25, the resulting plane becomes a single connected component as
shown in k of Figure 2. After acquiring all the planes for a particular protein structure, our next aim is to
find the skeletons for each of the plane shapes. Skeleton of the plane- k is shown in 𝑙 of Figure 2.
If we stack the skeletons for all the planes over each other, then the resulting image gives us an idea of
how the atoms form the overall protein structure in terms of the planes, which are formed by the
coordinates of the atoms. This stacked skeleton is basically a projection of the skeleton of the tertiary
structure of the protein. For the protein 2LEP the skeleton structure is shown Figure 3 as given below,
From the skeleton we have an idea of fractal-like distribution of protein atoms of the 3D protein structure
in the form of plane. Now we can compute the fractal dimension 𝐷𝑝 of the skeleton of the corresponding
tertiary protein structure and use it as the feature of tertiary protein structure. Fractal Dimension 𝐷𝑝 of a
group of protein molecules are given in Table 1 irrespective of their residue.

(𝑎)

(𝑏)

(𝑑)

(𝑔)

(𝑐)

(𝑒)

(𝑕)
6

(𝑓)

(𝑖)

(𝑘)

(j)

(l)

Figure 2. Representation of different iteration of multi-scale opening on One Slice of 2LEP
protein (a to k) and l is the Skeleton of one Slice of 2LEP protein.

Figure 3: The Stacked Skeleton of 2LEP protein.

Table 1: Fractal Dimension of Protein Molecules.
Protein ID
3v2j
3smk
3t0o
4ecs
3v2m
3sy1
4ag2
1cah
1cai
4bij
2lep
1cgi
4eym
2cbc

3.0.2.

Residue
Number
260
236
238
435
260
190
226
259
259
476
69
245
371
260

𝐷𝑝
1.661190e+000
1.620140e+000
1.646469e+000
1.649489e+000
1.656160e+000
1.605381e+000
1.695859e+000
1.661085e+000
1.660481e+000
1.680213e+000
1.549605e+000
1.635992e+000
1.649456e+000
1.661399e+000

Geodesic Dilation and Its Quantification.

From the PDB database the protein structures are viewed by using JMol and the protein structures are
rotated depending on the 3-axis, from which we have collected the 6- faces or views (front, left, right,
7

top, bottom, and back) of each 3D protein structure respectively. To find out self similarity between two
2D images we use geodesic dilation which is a morphological transformation to operate only some part of
the image (as marker) to grow until the boundary of the image and the advantages of this transformation
is that the structuring element can grow at each pixel, according to the image.
The Geodesic Dilation 𝜹𝑿 of an image Y inside X is defined as the intersection of the dilation of Y (with
respect to a structuring element S) with the image X
𝛿 𝑛 𝑋 = 𝑌⨁𝑛𝑆 ∩ 𝑋 𝑤𝑕𝑒𝑟𝑒 𝑛 = 1, 2, … , 𝑁
So Geodesic dilation terminates when all the connected components of X are constructed i.e. idem
potency is reached ∀ 𝑛 > 𝑛0 , 𝛿 (𝑛)𝑋 (𝑌) = 𝛿 (𝑛 0 )𝑋 (𝑌) . The following step are used for calculating
Geodesic Dilation 𝛿𝑝 of the corresponding faces of protein molecules.
Let 𝑓1 , 𝑓2 , … , 𝑓6

and

𝑔1 , 𝑔2 , … , 𝑔6

are six 2D faces of two 3D protein structures

𝑝𝑠 (𝑠𝑜𝑢𝑟𝑐𝑒 𝑝𝑟𝑜𝑡𝑒𝑖𝑛) and 𝑝𝑡 (𝑡𝑎𝑟𝑔𝑒𝑡 𝑝𝑟𝑜𝑡𝑒𝑖𝑛) .
1. The marker
2. 𝛿

𝑛

𝑝𝑠

𝐼𝑖 = 𝑓𝑖 ∩ 𝑔𝑖 , for all 𝑖 = 1, 2, … , 6.

= 𝐼𝑖 ⊕ 𝑛𝑆 ! = 𝑓𝑖 , for all 𝑛 = 1, 2, … , 𝑁, where S is a structuring element.

3. 𝛿 𝑛 𝑝𝑡 = 𝐼𝑖 ⊕ 𝑛𝑆 ! = 𝑔𝑖 , for all 𝑛 = 1, 2, … , 𝑁.
4. The number of dilation

𝛿𝑝 =

6
𝑛
𝑖=1 | 𝛿 𝑝 𝑠

− 𝛿 𝑛 𝑝𝑡 | .

Illustration:
Now we take the front view of two different proteins to compute the global similarity between them. For,
this purpose we consider the front view of 3V2J, 3V2M and the common structural part of both the
protein molecules, which are given in Figure 3. Here we use Geodesic dilation as a parameter for
determining the structure similarity between those protein molecules. Now we can determine the number
of dilation required from the intersection part to both the protein molecules towards constructing the self
similar structure, i.e. the number of geodesic dilation from the marker 𝐼 = 3𝑉2𝐽 ∩ 3𝑉2𝑀 towards 3V2J
and 3V2M is four for each. Similarly for other faces are given in Table 2. So the number of geodesic
dilation 𝛿𝑝 = 2 with respect to each faces, which shows that 3V2J and 3V2M have more structural
similar, whereas for 2LE8 and 2LLS have less structural similarity as 𝛿𝑝 = 64. The Table 2 given below
shows the different geodesic dilation of different protein molecules.

8

Front view of 3V2J

Front view of 3V2M

intersection of 3V2J and 3V2M

Figure 3: The backbone and intersection part of the protein molecules.

Table 2: Geodesic Dilation δ of different faces
Protein
ID
3V2J
3V2M
2LE8
2LLS
1CAH
2CBC
1CAH
4EYM
3V2J
3T0O

4.

Marker I

3V2J ∩ 3V2M
2LE8 ∩ 2LLS
1CAH ∩ 2CBC
1CAH ∩ 4EYM
3V2J ∩ 3T0O

Geodesic Dilation δ of different faces
Front
4
4
10
19
1
1
10
20
12
4

Left
4
4
6
19
1
1
10
6
8
5

Right
5
4
6
21
1
1
9
6
9
5

Top
6
6
7
18
1
2
18
9
12
5

Bottom
5
4
8
15
1
1
18
11
12
11

Back
3
3
10
19
1
1
11
9
11
5

δp
2
64
1
35
35

Result and Discussion.

In our experiment, we downloaded protein molecules from PDB and the result shows that our
methodology performance quite well for comparing tertiary protein structure in intricate level. The
similarity between two protein structures i and j can be computed by using the following equation:
𝜌 = |𝐷𝑝 𝑖 − 𝐷𝑝 𝑗 |.
and Geodesic Dilation 𝛿𝑝 , where 𝜌 the difference between the fractal is dimensions of any two protein
molecules and some experimental results are shown in Table 3. The difference between the fractal
dimensions is essentially measure the difference between the structural complexities. As 𝜌 approaches to
zero, the structures closed to be similar. The experimental result shows that if 𝜌 ≤ 0.008 and 𝛿𝑝 ≤ 12 ,
two protein molecules are similar in structures and functions. Thus, lower difference between fractal
dimensions and Geodesic dilation will ensure high similarity between the proteins which are being
compared. This would become clearer with few examples. For the same from the Table 3 we conclude

9

that the proteins molecule 1cah is more similar to 1cai and 2cbc as their fractal dimension difference 𝜌
=0.000604, 0.000314 and geodesic dilation 𝛿𝑝 = 5 𝑎𝑛𝑑 1 respectively.
Table 3: Difference between Fractal Dimensions of Compared Protein Pairs
Protein ID 1

3v2j

1cah

FD 𝐷𝑝𝑖

Protein-ID 2
3smk
3t0o
4ecs
1.661190e+000
3v2m
3sv1
4ag2
1cai
4bij
2lep
1.661085e+000
1cgi
4eym
2cbc

FD 𝐷𝑝𝑗
1.620140e+000
1.646469e+000
1.649489e+000
1.656160e+000
1.605381e+000
1.695859e+000
1.660481e+000
1.680213e+000
1.549605e+000
1.635992e+000
1.649456e+000
1.661399e+000

ρ
0.04105
0.014721
0.011701
0.00503
0.055809
0.034669
0.000604
0.019128
0.11148
0.025093
0.011629
0.000314

𝛿𝑝
24
35
31
2
28
39
5
21
43
35
23
1

PDB result
12%
18%
2%
100%
28%
31%
100%
41%
57%
50%
39%
100%

4. Conclusion:
In this work, we presented a novel technique to compute the structural similarity of 3D protein structure
using fractal dimension and geodesic dilation in atom levels and proteins backbone structure level
respectively. Compared with the existing methods, fractal dimension and geodesic dilation is easy to
compute and efficient enough to eliminate the limitations encountered in the existing algorithms. In our
experiments, atoms of all the protein structures are divided into slices by fixing the z co-ordinate value.
So only the analysis of the x-y planes is done. This work can be further extended by fixing the x or y coordinate values, i.e. analysis of the x-z and y-z planes of the protein structure. Thus, this article is allowing
us to see the intricate similarity and dissimilarity between tertiary protein structures with enhanced
results.

References:
[1]

G. Lancia, S. Istrail, “Protein Structure Comparison: Algorithms and Applications”. In: Guerra,
C., Istrail, S. (eds.) Mathematical Methods for Protein Structure Analysis and Design. LNCS
(LNBI), Vol. 2666, pp. 1–33. Springer, Heidelberg (2003).

[2]

L. P. Chew; “Exact Computation of Protein Structure Similarity”, Proceeding SCG '06
Proceedings of the twenty-second annual symposium on Computational geometry pp.468474(2006).
10

[3]

J. Galgonek, D. Hoksza, T. Skopal; “SProt: sphere-based protein structure similarity algorithm”.
Galgonek et al. Proteome Science, Vol.-9(Suppl 1):S20 (2011).

[4]

N. Krasnogor, D. A. Pelta; “Measuring the similarity of protein structures by means of the
universal similarity metric”; Bioinformatics, Vol.- 20,pp.1015-1021(2003).

[5]

P. Koehl, “Protein structure similarities”; Curr Opin Struct Biol, Vol.-11, pp.348-353(2001).

[6]

L. Holm, C. Ouzounis, C. Sander, G. Tuparev, G. Vriend; “A database of protein structure
families with common folding motifs”; Vol.-1, Issue-12, pp.1691-1698(1992).

[7]

A. Zemla, “LGA: a method for finding 3D similarities in protein structures “; Nucleic Acids
Research, Vol.-31, Issue-13, pp. 3370-3374 (2003).

[8]

D. Goldman, C. H. Papadimitriouz, S. Istraily; “Algorithmic Aspects of Protein Structure
Similarity". FOCS’ 99 Proceedings of the 40th Annual Symposium on Foundations of Computer
Science, pp. 512.(1999).

[9]

F Guyon, P. Tufféry, “Assessing 3D scores for protein structure fragment mining “, Open Access
Bioinformatics, Vol.-2, pp-67–77(2010).

[10]

I.G. Choi, J. Kwon, S. H. Kim; Local feature frequency profile: A method to measure structural
similarity in proteins; PNAS, Vol. 101; pp.-11(2004).

[11]

G. J. Kleywegt, “Experimental assessment of differences between related protein crystal
structures”; Acta Crystallogr. D Biol. Crystallogr, Vol.-55, pp. 1878-1884(1999).

[12]

C Cui, D. Wang, X. Yuan,“3D protein structures similarity matching based on fractal features”,
SPIE Vol. 5637, pp-567-572(2005).

[13]
[14]
[15]
[16]

L. L. Teo, B. S. Daya Sagar; “Modeling, description, and characterization of fractal pore via
mathematical morphology”; Discrete Dynamics in Nature and Society, Vol. 2006, Article ID
89280; Pages 1–24 (2006).
G. Klette, “Skeletons in digital image processing”, July, (2002).
Petros A. Maragos, Ronald W. Schafer; “Morphological skeleton representation and coding of
binary images”; IEEE Transactions on Acoustics, Speech and Signal processing, Vol. assp-34;
pp. 5(1986).
D. Avnir, O. Biham, D. Lidar, O. Malcai,"Is the geometry of Nature fractal", Science Vol.-279, pp.
39-40(1998).

[17]

K Develi, T Babadagli, "Quantification of natural fracture surfaces using fractal geometry", Math.
Geology, Vol.- 30 , pp. 971-998(1998).

11

