International Journal in Foundations of Computer Science & Technology (IJFCST), Vol.5, No.3, May 2015

AGENT ENABLED MINING OF DISTRIBUTED
PROTEIN DATA BANKS
G. S. Bhamra1, A. K. Verma2 and R. B. Patel3
1

M. M. University, Mullana, Haryana, 133207 - India
2
Thapar University, Patiala, Punjab, 147004- India
3
Chandigarh College of Engineering & Technology, Chandigarh- 160019- India

ABSTRACT
Mining biological data is an emergent area at the intersection between bioinformatics and data mining
(DM). The intelligent agent based model is a popular approach in constructing Distributed Data Mining
(DDM) systems to address scalable mining over large scale distributed data. The nature of associations
between different amino acids in proteins has also been a subject of great anxiety. There is a strong need to
develop new models and exploit and analyze the available distributed biological data sources. In this study,
we have designed and implemented a multi-agent system (MAS) called Agent enriched Quantitative
Association Rules Mining for Amino Acids in distributed Protein Data Banks (AeQARM-AAPDB). Such
globally strong association rules enhance understanding of protein composition and are desirable for
synthesis of artificial proteins. A real protein data bank is used to validate the system.

KEYWORDS
Knowledge Discovery, Association Rules, Intelligent Agents, Multi-Agent System, Bioinformatics.

1. INTRODUCTION
Data Mining (DM) is a process to automatically extract some interesting and valid data patterns
or trends representing knowledge, implicitly stored in large databases [1], [2]. Distributed Data
Mining (DDM) is concerned with application of classical DM procedures in a distributed
computing environment trying to make the best of the available resources. In DDM, DM takes
place both locally at each geographically distributed site and at a global level where the local
knowledge is merged in order to discover global knowledge. A DDM system is a very complex
entity that is comprised of many components; mining algorithms; communication subsystems;
resources management; task scheduling; user interface etc. It should provide efficient access to
both distributed data and computing resources; monitor the entire mining procedure; and present
results to users in appropriate formats. A successful DDM system is also flexible enough to adapt
to various situations [3], [4], [5], [6], [7].
Intelligent software agent technology is an interdisciplinary technology. The motivating idea of
this technology is the development and efficient utilization of autonomous software objects called
agents, which have access to geographically distributed and heterogeneous information resources
to simplify the complexities of distributed computing. They are autonomous, adaptive, reactive,
pro-active, social, cooperative, collaborative and flexible. They also support temporal continuity
and mobility within the network. An intelligent agent with mobility feature is known as Mobile
Agent (MA). MA migrates from node to node in a heterogeneous network without losing its
operability. It can continue to function even if the user is disconnected from the network. On
reaching at a node MA is delivered to an Agent Execution Environment (AEE) where its
executable parts are started running. Upon completion of the desired task, it delivers the results to
DOI:10.5121/ijfcst.2015.5303

25

International Journal in Foundations of Computer Science & Technology (IJFCST), Vol.5, No.3, May 2015

the home node. With MA, a single serialized object is transmitted over the network carrying the
small amount of resultant data only thus reducing the consumption of network bandwidth, latency
(response time delay) and network traffic. An AEE or Mobile Agent Platform (MAP), is server
application that provides the appropriate functionality to MAs to authenticate, execute,
communicate (with other agents, users, and other platforms), migrate to other platform, and use
system resources in a secure way. A Multi Agent System (MAS) is distributed application
comprised of multiple interacting intelligent agent components [8].
Table 1. Single Letter codes for Amino Acids

Sr.
No.

Amino Acid Name

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20

Alanine
Cysteine
Aspartic Acid
Glutamic Acid
Phenylalanine
Glycine
Histidine
Isoleucine
Lysine
Leucine
Methionine
Asparagine
Proline
Glutamine
Arginine
Serine
Threonine
Valine
Tryptophan
Tyrosine

Single
Letter
Code
A
C
D
E
F
G
H
I
K
L
M
N
P
Q
R
S
T
V
W
Y

Three
Letter
Code
Ala
Cys
Asp
Glu
Phe
Gly
His
Ile
Lys
Leu
Met
Asn
Pro
Gln
Arg
Ser
Thr
Val
Trp
Tyr

Bioinformatics or computational molecular biology aims at automated analysis and the
management of high-throughput biological data as well as modelling and simulation of complex
biological systems. Bioinformatics has very much changed since the first sequence alignment
algorithms in 1970s [9]. Today in-silico analysis is a fundamental component of biomedical
research. Bioinformatics has now encompasses a wide range of subject areas from structural
biology, genomics to gene expression studies [10], [11]. The post-genomic era has resulted into
availability of the enormous amount of distributed biological data sets that require suitable tools
and methods for modelling and analyzing biological processes and sequences. The bioinformatics
research community feels a strong need to develop new models and exploit and analyze the
available genomes [12]. The protein sequences are made up of 20 types of Amino Acids (AA).
Each AA is represented by a single letter alphabet as shown in Table 1. Unique 3-dimensional
structure of each protein is decided completely by its amino-acid sequence. Proteins are important
constituents of cellular machinery of any organism and the functioning of proteins heavily
depends upon its AA sequence. A slight change in the sequence might completely change the
functioning of the protein. Therefore, the nature of associations between different amino acids has
been a subject of great anxiety. In [20], authors mine the association rules among amino acids in
protein sequences. As per the literature available, no researcher has applied the agent technology
to mine association rules for amino acids from distributed protein data sets. The amalgamation of
26

International Journal in Foundations of Computer Science & Technology (IJFCST), Vol.5, No.3, May 2015

the DDM and MAS provides rewarding solution in terms of security, scalability, storage cost,
computation cost and communication cost. Mining biological data is an emerging area and
continues to be an extremely important problem, both for DM and for biological sciences [21].
State of the art in use of agent technology in bioinformatics is reviewed in [22].
The rest of the paper is organised as follows. Section 2 discusses the distributed association rule
mining and preliminary notations used in this paper. A running environment for the proposed
system is present in Section 3 along with algorithms used for various components and protein
data banks used in this study. The layout and working of various agents involved in the proposed
system and their algorithms are also discussed. System is implemented in Java and its
performance is studied in Section 4 and finally the article is concluded in Section 5.

2. DISTRIBUTED ASSOCIATION RULE MINING

{

}

Let DB = T j , j = 1K D be a transactional dataset of size D where each transaction T is assigned
an identifier ( TID ) and I = {di ,i = 1K m} , total m data items in DB . A set of items in a particular
transaction T is called itemset or pattern. An itemset, P = {d i ,i = 1K k } , which is a set of k data
items in a particular transaction T and P ⊆ I , is called k-itemset. Support of an itemset,
No_of_T_containing_P
s ( P) =
% is the frequency of occurrence of itemset P in DB , where
D
No_of_T_containing_P is the support count (sup_count) of itemset P . Frequent Itemsets (FIs)

are the itemset that appear in DB frequently, i.e., if s ( P ) ≥ min_th_sup (given minimum
threshold support), then P is a frequent k-itemset. Finding such FIs plays an essential role in
miming the interesting relationships among itemsets. Frequent Itemset Mining (FIM) is the task
of finding the set of all the subsets of FIs in a transactional database. It is CPU and input/output
intensive task, mainly because of the large size of the datasets involved [2].

Association Rules (ARs) first introduced in [13], are used to discover the associations (or cooccurrences) among item in a database. AR is an implication of the form
P ⇒ Q [support,confidence ] where, P ⊂ I , Q ⊂ I and P and Q are disjoint itemsets, i.e.,
P ∩ Q = ∅ . An AR is measured in terms of its support ( s ) and confidence ( c ) factors. An AR
P ⇒ Q is said to be strong if s ( P ⇒ Q ) ≥ min_th_sup (given minimum threshold support) and
c ( P ) ≥ min_th_conf (given minimum threshold confidence). Association Rule Mining (ARM)
today is one of the most important aspects of DM tasks. In ARM all the strong ARs are generated
from the FIs. The ARM can be viewed as two step process [14], [15].
1. Find all the frequent k-itemsets( Lk )
2. Generate Strong ARs from Lk
a. For each frequent itemset, l ∈ Lk , generate all non empty subsets of l .
b. For every non empty subset s of l , output the rule “ s ⇒ ( l − s ) ”, if
sup_count ( l )
sup_count ( s )

≥ min_th_conf

Distributed Association Rule Mining (DARM) generates the globally strong association rules
from the global FIs in a distributed environment. Because of an intrinsic data skew property of the
distributed database, it is desirable to mine the global rules for the global business decisions and

27

International Journal in Foundations of Computer Science & Technology (IJFCST), Vol.5, No.3, May 2015

the local rules for the local business decisions. Comparative analysis of existing agent based
DARM systems can be found in [23].
Few preliminaries notations and definitions required for defining DARM and to make this study
self contained are as follows:

•

S = {Si ,i = 1K n} , n distributed sites.

•
•

SCENTRAL , Central Site.
PDBi = {PR m , m = 1K X i } , Protein Data Bank of X i protein records at site Si , each
PRm has two main parts; first part is the description headers ( PD ) for structural
classification of protein and the second part contains protein sequence ( PS ), i.e.,
the chain of amino acids. Snapshot of PDB1 is shown in Appendix A.1.
FPDBi = {PR k ,k = 1K Di } , Filtered Protein Data Bank of Di protein records at site Si ,
where length of PS in each PR is in the range ≥ 50 and < 400 . Snapshot of FPDB1 is
shown in Appendix A.2.
AAFi , Data bank of amino acids frequency for each PR ∈ FPDB at site Si . AAF1 is
shown in Appendix A.3. AAF2 and AAF3 are not shown due to space constraint.
BDBi , Boolean Data Bank that contains a value ‘1’ if the frequency of an amino acid lies
in the specific range and ‘0’ otherwise at site Si . Snapshot of BDB1 is shown in
Appendix A.5.
IDBi , Itemset Data Bank of frequency partition items to map boolean value ‘1’ with its
frequency partition item number at site Si . Snapshot of IDB1 is shown in Appendix A.6.

•

•
•

•
•

LFI
k ( i ) , Local frequent k-itemsets at site Si .

•

FI
LFISC
k ( i ) , List of support count ∀Itemset ∈ Lk ( i ) .

•

LLSAR
, List of locally strong association rules at site Si .
i

•

LTLSAR =

•

LTFI
k

•

LGFI
k

•

LGSAR
CENTRAL , List of Globally strong association rule.

n
LSAR
i =1 i

U L
=U L
=I L

n
FI
i =1 k ( i )
n
FI
i =1 k ( i )

, List of total locally strong association rules.

, List of total frequent k-itemsets.
, List of global frequent k-itemsets.

FISC
LSAR
Local Knowledge Base (LKB), at site Si , comprises of LFI
which can provide
k ( i ) , Lk ( i ) and Li

reference to the local supervisor for local decisions. Global Knowledge Base (GKB), at SCENTRAL ,
GFI
comprises of LTLSAR , LTFI
and LGSAR
k , Lk
CENTRAL for the global decision making. Like ARM, DARM
task can also be viewed as two-step process [15]:

1. Find the global frequent k-itemset ( LGFI
) from the distributed Local frequent k-itemsets
k
( LFI
k ( i ) ) from the partitioned datasets.
GFI
2. Generate globally strong association rules ( LGSAR
.
CENTRAL ) from Lk

28

International Journal in Foundations of Computer Science & Technology (IJFCST), Vol.5, No.3, May 2015

3. PROPOSED AEQARM-AAPDB SYSTEM
3.1. Environment for the proposed system
Every MAS needs an underlying AEE to provide running infrastructure on which agents can be
deployed and tested. A running environment has been designed in Java to execute all the DM
agents involved in the proposed system. Various attributes of the MA are encapsulated within a
data structure known as AgentProfile . AgentProfile contains the name of MA ( AgentName ),
version number ( AgentVersion ), entire byte code ( BC ), list of nodes to be visited by MA, i.e.,
itinerary plan ( LNODES ) , type of the itinerary ( ItinType ) which can be serial or parallel, a
reference of current execution state ( AObject ) and an additional data structure known as
Briefcase that acts as a result bag of MA to store final resultant knowledge ( Result_Si ) at a
particular site. Computational time ( CPUTime ) taken by a MA at a particular site is also stored
in Result_Si . In addition to results, Briefcase also contains the system time for start of agent
journey ( TripTimestart ), system time for end of journey ( TripTimeend ) and total round trip time of
MA ( TripTime ) calculated using the formula TripTime ← TripTimeend − TripTimestart . This
environment consists of the following three components:

•

•

Data Mining Agent Execution Environment (DM_AEE): It is the key component that
acts as a Server. DM_AEE is deployed on any distributed sites Si and is responsible for
receiving, executing and migrating all the visiting data mining (DM) agents. It receives
the incoming AgentProfile at site Si , retrieves the entire BC of agent and save it with
AgentName.class in the local file system of the site Si after that execution of the agent is
started using AObject . Steps are shown in Algorithm 1.
Agent Launcher (AL): It acts a Client at agent launching station ( SCENTRAL ) and launches
the goal oriented DM agents on behalf of the user through a user interface to the
DM_AEE running at the distributed sites. Agent Pool (or Zone) at SCENTRAL is a repository
of all mobile as well as stationary agents (SAs). AL first reads and stores AgentName in
AgentProfile . The entire BC of the AgentName is loaded from the Agent Pool and
stored in AgentProfile . LNODES and ItinType are retrieved and stored in AgentProfile .
TripTimestart is maintained in Briefcase which is further added to AgentProfile . In case

•

of parallel computing model, i.e., if ItinType = Parallel AL creates LNODES number of
clones of the specific MA and dispatches each clone in parallel to all the sites listed in
LNODES . Here LNODES number of threads are created to dispatch the MAs in parallel. Each
clone has AgentVersion starting from 1 to size of LNODES , which is used to identify each
clone on the network. Before dispatching the clone of a MA to DM_AEE, the current
state of the newly created ith clone object ( AObject[i ] ) is also stored in AgentProfile .
AL also contacts the Result Manager (RM) for processing the Briefcase of an agent.
Detailed steps are given in Algorithm 2.
Result Manager (RM): It manages and processes the Briefcase of all MAs. RM is either
contacted by a MA for submitting its results or by AL for processing the results of the
specific MA. On completion of itinerary, each DM agent submits its results to RM which
computes total round trip time ( TripTime ) of that MA and saves it in the Briefcase of
that agent. If ItinType = Parallel then it saves the AgentProfile of all the clones of the
agent with AgentName in a collection LAllProfiles
AgentName . It is assumed that all the clones report
29

International Journal in Foundations of Computer Science & Technology (IJFCST), Vol.5, No.3, May 2015

their results to RM. RM may be equipped with the feature of non reporting clones by
issuing an alert to AL for that clone with specific AgentVersion . AL then launch a new
clone for the specific AgentVersion for the specific site. When it is contacted by AL for
processing the results of a specific agent it sends back a collection LAllProfiles
AgentName for all the
clones of that agent. Steps are defined in Algorithm 3.
Algorithm 1 DATA MINING AGENT EXECUTION ENVIRONMENT (DM_AEE)
1: procedure DM_AEE( )
2: while TRUE do
3:
AgentPofile ← listen and receive AgentProfile at Si
4:
AgentName ← get AgentName from AgentProfile
5:
BC ← retrieve the BC of agent from AgentProfile
6:
save the BC with AgentName.class in the local file system of Si

> current state

7:
AObject ← get AObject from AgentProfile
8:
AObject .run()
9:
end while
10: end procedure

> start executing mobile agent

Algorithm 2 AGENT LAUNCHER (AL) FOR AEQARM-AAPDB
1: procedure AL( )
2: option ← read option (dispatch / result)
3: switch option do
4:

> dispatch the mobile agent to DM_AEE

case dispatch

5:
6:

AgentName ← read Mobile Agent's name
BC ← load entire byte code of AgentName from AgentPool

7:

add AgentName and BC to AgentProfile

8:
9:

LNODES ← read Itinerary (IP addresses) of mobile agent
ItinType ← read ItinType ( Serial / Parallel)

10:
11:

add ItinType to AgentProfile
if ItinType = "Parallel" then
NODES

> Parallel Itinerary
> Array of Agent Objects for clone references

12:

AObject[L

13:
14:

TripTimestart ← get system time for start of the agent journey
add TripTimestart to Briefcase

15:
16:
17:
18:
19:

add Breifcase to AgentProfile
switch AgentName do
case PDBFA
for i ← 1, LNODES .size do
AgentVersion ← i
add AgentVersion to AgentProfile

20:

.size]

21:

NodeAddress ← LNODES .get(i)

22:

AObject[i] ← new PDBFA(AgentProfile)
Add AObject[i] to AgentProfile

23:
24:
25:

> for each node in the itinerary

> get an IP address
> clone’s state

Transfer AgentProfile to DM_AEE at NodeAddress
end for
30

International Journal in Foundations of Computer Science & Technology (IJFCST), Vol.5, No.3, May 2015
26:
27:
28:
29:

end case
case AAFFA
for i ← 1, LNODES .size do
AgentVersion ← i

30:

add AgentVersion to AgentProfile

31:

NodeAddress ← LNODES .get(i)

32:
33:

AObject[i] ← new AAFFA(AgentProfile)
Add AObject[i] to AgentProfile

34:
35:
36:
37:
38:
39:
40:

> for each node in the itinerary

> get an IP address
> clone’s state

Transfer AgentProfile to DM_AEE at NodeAddress
end for
end case
case FMIDBGA
max_freq ← read maximum frequency range for amino acids
for i ← 1, LNODES .size do
AgentVersion ← i

> for each node in the itinerary

41:

add AgentVersion to AgentProfile

42:

NodeAddress ← LNODES .get(i)

43:
44:

AObject[i] ← new FMIDBGA(AgentProfile, max_freq)
Add AObject[i] to AgentProfile

45:
46:
47:
48:
49:
50:
51:
52:

> get an IP address

min_c ← read minimum threshold confidence
for i ← 1, LNODES .size do
AgentVersion ← i

> for each node in the itinerary

53:

add AgentVersion to AgentProfile

54:

NodeAddress ← LNODES .get(i)
AObject[i] ← new LKGA_P(AgentProfile,min_s,min_c)

55:
56:
57:
58:
59:
60:
61:
62:

for i ← 1, LNODES .size do
AgentVersion ← i
add AgentVersion to AgentProfile

64:
65:

NodeAddress ← LNODES .get(i)

71:

> get an IP address
> clone’s state

Add AObject[i] to AgentProfile
Transfer AgentProfile to DM_AEE at NodeAddress
end for
end case
case LKCA_P

63:

66:
67:
68:
69:
70:

> clone’s state

Transfer AgentProfile to DM_AEE at NodeAddress
end for
end case
case LKGA_P
min_s ← read minimum threshold support

> for each node in the itinerary

AObject[i] ← new LKCA_P(AgentProfile)
Add AObject[i] to AgentProfile

> get an IP address
> clone’s state

Transfer AgentProfile to DM_AEE at NodeAddress
end for
end case
case GKDA_P
CENTRAL
LGSAR
CENTRAL ← load LCANTRAL generated by GKGA at SCENTRAL

31

International Journal in Foundations of Computer Science & Technology (IJFCST), Vol.5, No.3, May 2015
72:
73:

add LGSAR
CENTRAL to Briefcase
add updated Briefcase to AgentProfile

74:
75:
76:

for i ← 1, LNODES .size do
AgentVersion ← i
add AgentVersion to AgentProfile

77:

NodeAddress ← LNODES .get(i)

78:
79:
80:
81:
82:
83:
84:
85:
86:
87:
88:
89:
90:

add AgentName and ItinType to LAgentInfo
if ItinType = " Parallel " then
Pr ofile
AgentInfo
LAll
AgentName ← contact RM for L

92:
93:
94:

switch AgentName do
case PDBFA
for all AgentProfile ∈ LAllProfile
AgentName do

101:
102:
103:
104:
105:
106:
107:
108:
109:
110:
111:
112:
113:
114:
115:
116:
117:
118:

> get an IP address

AObject[i] ← new GKDA_P(AgentProfile)
Add AObject[i] to AgentProfile
> clone’s state
Transfer AgentProfile to DM_AEE at NodeAddress
end for
end case
end switch
end if
end case
case result
> process the results of mobile agent
AgentName ← read mobile agent's name
ItinType ← read mobile agent's itinerary type

91:

95:
96:
97:
98:
99:
100:

> for each node in the itinerary

> for each clone

Briefcase ← retrieve Briefcase from AgentProfile

process the Briefcase of PDBFA clone
end for
end case
case AAFFA
for all AgentProfile ∈ LAllProfile
AgentName do

> for each clone

Briefcase ← retrieve Briefcase from AgentProfile

process the Briefcase of AAFFA clone
end for
end case
case FMIDBGA
for all AgentProfile ∈ LAllProfile
AgentName do

> for each clone

Briefcase ← retrieve Briefcase from AgentProfile

process the Briefcase of FMIDBGA clone
end for
end case
case LKGA_P
for all AgentProfile ∈ LAllProfile
AgentName do

> for each clone

Briefcase ← retrieve Briefcase from AgentProfile

process the Briefcase of LKGA_P clone
end for
end case
case LKCA_P
call RIGKGA(LAllProfile
AgentName )

> stationary agent
32

International Journal in Foundations of Computer Science & Technology (IJFCST), Vol.5, No.3, May 2015
119:
120:
121:

end case
case GKDA_P
for all AgentProfile ∈ LAllProfile
AgentName do

> for each clone

122:
Briefcase ← retrieve Briefcase from AgentProfile
123:
process the Briefcase of GKDA_P clone
124:
end for
125:
end case
126:
end switch
127:
end if
128:
end case
129:
end switch
130: end procedure

Algorithm 3 RESULT MANAGER(RM)
1: procedure RM( )
2: while TRUE do
3:
listen and receive the incomming request
4:
if contacted by a mobile agent for submitting results from site Si then
5:
AgentProfile ← receive the incomming AgentProfile from site Si
6:
ItinType ← retrieve ItinType from AgentProfile

7:
8:
9:
10:
11:
12:
13:
14:
15:
16:

Briefcase ← retrieve mobile agent's Briefcase from AgentProfile

TripTimestart ← retrieve TripTimestart from Briefcase
TripTimeend ← retrieve TripTimeend from Briefcase
TripTime ← TripTimeend − TripTimestart
add TripTime to Briefcase
add updated Briefcase to AgentProfile
if ItinType = "Parallel" then
AgentName ← retrieve AgentName from AgentProfile
AgentVersion ← retrieve AgentVersion from AgentProfile
if AgentVersion = 1 then

17:

add AgentProfile to LAllProfiles
AgentName

18:

save LAllProfiles
AgentName at S CENTRAL

19:
20:

end if
if AgentVersion > 1 then

21:

retrieve LAllProfiles
AgentName from S CENTRAL

22:

add AgentProfile to LAllProfiles
AgentName

23:

save updated LAllProfiles
AgentName at SCENTRAL

24:
25:
26:
27:
28:
29:
30:
31:

end if
end if
end if
if contacted by AgentLauncher for processing the results then
AgentName ← retrieve AgentName from AgentProfile

ItinType ← retrieve ItinType from incomming LAgentInfo
if ItinType = " Parallel " then
AllProfiles
LAllProfiles
AgentName ← load LAgentName from SCENTRAL

33

International Journal in Foundations of Computer Science & Technology (IJFCST), Vol.5, No.3, May 2015
32:

dispatch LAllProfiles
AgentName to AgentLauncher

33:
end if
34:
end if
35: end while
36: end procedure

Figure 1. AeQARM-AAPDB MAS

3.2. Protein Data Bank
AeQARM-AAPDB system is tested on a real dataset of protein sequences, the Astral SCOP [16],
[17 ] version 1.75 genetic domain sequence subsets, based on PDB SEQRES records with less
than 40% identity to each other [18]. There are total of 10569 Protein records in this dataset. This
single PDB is divided into 3 unites ( PDB1 , PDB2 and PDB3 ) of 3523 protein records in each
and are stored at three distributed sites. Each PDBi is further filtered to generate FPDBi for PS
length range ≥ 50 and ≤ 400 . A total of only 9633 (3341 ( FPDB1 ) + 3253 ( FPDB2 ) + 3039
( FPDB3 )) such filtered records are considered for the mining. The frequencies of 20 amino acids
in each protein sequence are retrieved and stored in AAFi which is further mapped into BDBi for

34

International Journal in Foundations of Computer Science & Technology (IJFCST), Vol.5, No.3, May 2015

each amino acid having 15 frequency ranges (partitions) resulting into 300 amino acids
attribute, value pairs as shown in Appendix A.4.

3.3. Layout and working of AeQARM-AAPDB system
AeQARM-AAPDB MAS is shown in Figure 1. This MAS consists of total seven agents, clones
of six MAs in serial number 1 to 6 are dispatched from SCENTRAL with parallel itinerary migration
and one at serial number 7 is an intelligent stationary agent (SA) running at SCENTRAL to perform
the different tasks. The CPU time taken by a MA while processing on each site along with some
other specific information is carried back in the result bag at SCENTRAL . Relationship among these
agents and their working behaviour are given as follows.
1. Protein Data Bank Filtering Agent (PDBFA): Clones of this MA is dispatched in parallel
to each distributed site by AL. It carries the AgentProfile along with it and filters PDBi to
generate FPDBi at each site Si . PDBFA carries back the computational time ( CPUTime )
at each site Si and TripTimeend .
2. Amino Acids Frequency Finder Agent (AAFFA): Every clone of this MA carries the
AgentProfile along with it and finds the frequencies of each amino acids in every protein
sequence record in FPDBi to create AAFi at each site Si . AAFFA carries back the
computational time ( CPUTime ) at each site Si and TripTimeend .
3. Frequency Mapping and Itemset Data Bank Generater Agent (FMIDBGA): Every clone
of this MA carries the AgentProfile and maxfrq (the given maximum frequency range for
amino acids) along with it. It divides the frequencies of each amino acids in AAFi into
intervals and maps the frequencies into Boolean values to create BDBi for frequency
intervals and further maps it to IDBi for frequency interval items.
4. Local Knowledge Generator Agent with Parallel Itinerary (LKGA_P): Every cloned
LKGA_P carries the AgentProfile , min_th_sup and min_th_conf along with it. This
agent is embedded with Apriori algorithm [19] for generating all the frequent k-itemset
lists. This agent may be equipped with decision making capability to select other FIM
algorithms based on the density of the dataset at a particular site. It first performs the FIM
FISC
to generate and store LFI
k (i) and Lk (i) at site Si by scanning the local IDBi at that site with
the constraint of min_th_sup . It then performs the ARM applying the constraint of
FISC
LSAR
by using the LFI
list also
min_th_conf to generate and store LLSAR
k (i) and Lk (i) . Li
i

contains support and confidence for a particular association rule along with site name. It
carries back the computational time ( CPUTime ) at each site Si and TripTimeend .
5. Local Knowledge Collector Agent with Parallel Itinerary (LKCA_P): Every cloned
LKGA_P carries the AgentProfile and collects the list of local frequent k-itemset ( LFI
k (i ) )
and list of locally strong association rules ( LLSAR
) generated by LKGA_P. It carries back
i
these distributed results in the result bag to RM at SCENTRAL where these results are
integrated with the help of a RIA stationary agent. In addition to this resultant knowledge it
also carries back the computational time ( CPUTime ) at each site Si and TripTimeend .
6. Global Knowledge Dispatcher Agent with Parallel Itinerary (GKDA_P): Every cloned
GKDA_P carries AgentProfile containing global knowledge ( LGSAR
CENTRAL ) for further

35

International Journal in Foundations of Computer Science & Technology (IJFCST), Vol.5, No.3, May 2015

decision making and comparing with local knowledge at that site. It carries back the
computational time ( CPUTime ) at each site Si and TripTimeend .
7. Result Integration and Global Knowledge Generater Agent (RIGKGA): It is a stationary
agent at SCENTRAL , mainly used for processing the result bags of all clones of LKCA_P. It
GFI
creates a list of total frequent k-itemset ( LTFI
and a
k ), a list of global frequent itemset, Lk

list of total locally strong association rules ( LTLSAR ). LTFI
and LTLSAR are further processed to
k
generate and store LGSAR
CENTRAL list.

Figure 2. Control Panel of AeQARM-AAPDB
Table 2. Network Configuration

LAN Configuration
IP a
Network
SCENTRAL
Intel b
MS c
192.168.46.5
NW d
S1
Intel b
MS c
192.168.46.212
NW d
b
c
S2
Intel
MS
192.168.46.189
NW d
S3
Intel b
MS c
192.168.46.213
NW d
a. IP address with Mask: 255.255.255.0 and Gateway 192.168.46.1
b. Intel Pentium Dual Core(3.40 GHz, 3.40 GHz) with 512 MB RAM
c. Microsoft Windows XP Professional ver. 2002
d. Network Speed: 100 Mbps and Network Adaptor: 82566DM-2 Gigabit NIC
Site Name

Processor

OS

4. IMPLEMENTATION AND PERFORMANCE STUDY
All the agents as well as Control Panel of the system, as shown in Figure 2, are implemented in
Java. The required configuration for the deployment of the system is shown in Table 2 with
additional deployment of DM_AEE at each distributed site and AL and RM at SCENTRAL .

36

International Journal in Foundations of Computer Science & Technology (IJFCST), Vol.5, No.3, May 2015

Figure 7. Globally strong association rules for globally frequent 2-itemsets

Figure 8. Globally strong association rules for frequent 2- amino acids
FISC
LFI
k (1) and Lk (1) at site S1 generated by LKGA_P agent with 20% min_th_sup are shown in

Appendix B.1. Locally strong association rules ( L1LSAR ) generated by LKGA_P for frequent item
numbers at site S1 are shown in Appendix B.2 and the same for their corresponding amino acids
frequency range are shown in Appendix B.3. Globally strong association rules ( LGSAR
CENTRAL ) for the
globally frequent itemsets generated by RIGKGA are shown in Figure 7. When these item
numbers are mapped with their corresponding amino acids frequency ranges then globally strong
37

International Journal in Foundations of Computer Science & Technology (IJFCST), Vol.5, No.3, May 2015

quantitative association rules for frequent 2 amino acids are obtained and shown in Figure 8. The
results, as shown in Figure 8, reveals that:
•
•

< Cysteine: 0..2 > is strongly associated with < Histidine: 3..5 >, < Methionine:0..2 >, <
Methionine: 3..5 >, < Tryptophan: 0..2 >, and < Tyrosine..3::5 >
< Tryptophan: 0..2 > is strongly associated with < Histidine: 0..2 >, <Histidine: 3..5 >, <
Methionine: 0..2 >, < Methionine: 3..5 >, < Glutamine: 3..5 > and < Tyrosine: 3..5 >

5. CONCLUSION
Mobile agents strongly qualify for designing distributed applications. DDM, when clubbed with
the agent technology, makes a promising alliance that gives favourable results and provides a
rewarding solution in managing Big Data with ever increasing size. In this study, a MAS called
AeQARM-AAPDB to mine the strong quantitative association rules among amino acids present
in primary structure of the proteins from the distributed proteins data sets using intelligent agents
is presented. Such globally strong association rules are used in understanding of protein
composition and are desirable for synthesis of artificial proteins. Agent-based bio-data mining
leaves the technical details of choosing mining algorithms, forming hybrid system, and preparing
specific data format to the intelligent system itself because such requirements are unreasonable
for most biologists. It alleviates the technical difficulty while enhance the reusability of the
mining algorithms and available datasets. By applying multi-agent based distributed bio-data
mining, the computing load can be balanced and the computational effort can be achieved in a
parallel and scalable manner.

REFERENCES
[1]

U. M. Fayyad, G. Piatetsky-Shapiro, P. Smyth & R. Uthurusamy, (1996) Advances in Knowledge
Discovery and Data Mining, AAAI/MIT Press.
[2] J. Han & M. Kamber, (2006) Data Mining: Concepts and Techniques, 2nd ed. Morgan Kaufmann.
[3] B. -H. Park & H. Kargupta, (2002) “Distributed Data Mining: Algorithms, Systems, and
Applications,” Department of Computer Science and Electrical Engineering, University of Maryland
Baltimore
County,
1000
Hilltop
Circle
Baltimore,
MD
21250,
Available:
http://www.csee.umbc.edu/_hillol/PUBS/review.pdf
[4] G. Tsoumakas & I. Vlahavas, (2009) “Distributed Data Mining”, Department of Informatics, Aristotle
University
of
Thessaloniki,
Thessaloniki,
Greece,
Available:
http://talos.csd.auth.gr/tsoumakas/publications/E1.pdf
[5] H. Kargupta & P. Chan, (2000) Advances in Distributed and Parallel Knowledge Discovery.
AAAI/MIT Press.
[6] Y. Fu, (2001) “Distributed Data Mining: An Overview”, Department of Computer Science, University
of Missouri-Rolla, Available: http://academic.csuohio.edu/fuy/Pub/tcdp01.pdf
[7] A. O. Ogunde, O. Folorunso, A. S. Sodiya & G. O. Ogunleye, (2011) “A Review of Some Issues and
Challenges in Current Agent-Based Distributed Association Rule Mining”, Asian Journal of
Information Technology, vol. 10, no. 2, pp. 84–95.
[8] G. S. Bhamra, R. B. Patel & A. K. Verma, (2014) “Intelligent Software Agent Technology: An
Overview”, International Journal of Computer Applications (IJCA), vol. 89, no. 2, pp. 19–31.
[9] A. Gibbs & G. McIntyre, (1970) “The diagram, a method for comparing sequences”, European
Journal of Biochemistry, vol. 16, pp.1-11.
[10] J. C. Setubal & J. Meidanis, (1997) “Introduction to Computational Molecular Biology”, PWS
Publishing Company.
[11] T. Attwood & D. Parry-Smith, (1999) “Introduction to Bioinformatics”, First ed. Prentice Hall.
[12] Z. Ezziane, (2006) “Applications of artificial intelligence in bioinformatics: A review”, Expert Systems with
Applications, vol. 30, pp. 2-10.

38

International Journal in Foundations of Computer Science & Technology (IJFCST), Vol.5, No.3, May 2015
[13] R. Agrawal, T. Imielinski & A. Swami, (1993) “Mining association rules between sets of items in

[14]
[15]
[16]
[17]

[18]
[19]

[20]

[21]
[22]
[23]

large databases”, in Proceedings of the ACM-SIGMOD International Conference of Management of
Data, pp. 207–216.
R. Agrawal & J. C. Shafer, (1996) “Parallel mining of association rules”, IEEE Transaction on
Knowledge and Data Engineering, vol. 8, no. 6, pp. 962–969.
M. J. Zaki, (1999) “Parallel and distributed association mining: a survey”, IEEE Concurrency, vol. 7,
no. 4, pp. 14–25.
Loredana Lo Conte, Steven E. Brenner, et al., (2002) “SCOP database in 2002: refinements
accommodate structural genomics”, Nucleic Acids Research, vol. 30, no. 1, pp. 264-267.
Alexey G. Murzin, Steven E. Brenner, et al., (1995) “SCOP: a structural classification of proteins
database for the investigation of sequences and structures”, Journal of Molecular Biology, vol. 247,
no. 4, pp.536-540.
AstralSCOP, (2009) Astral SCOP version 1.75 genetic domain sequence subsets with less than 40%
identity to each other, http://scop.berkeley.edu/astral/ver=1.75.
R. Agrawal and R. Srikant, (1994) “Fast Algorithms for Mining Association Rules in Large
Databases”, In Proceedings of the 20th International Conference on Very Large Data Bases
(VLDB'94), pp. 487-499.
Nitin Gupta, Nitin Mangal, Kamal Tiwari, & Pabitra Mitra, (2006) “Mining Quantitative Association
Rules in Protein Sequences”, In Graham J. Williams & Simeon J. Simeon, editors, Data Mining,
Lecture Notes in Computer Science, vol. 3755, pp. 273-281.
Qiang Yang & Xindong Wu, (2006) “10 Challenging Problems in Data Mining Research”,
International Journal of Information Technology & Decision Making, vol. 5, no. 4, pp. 597-604.
Gurpreet S. Bhamra, Anil K. Verma & Ram B. Patel, (2014) “Agent Technology in Bioinformatics: A
Review”, International Journal of Next-Generation Computing, vol. 5, no. 2, pp.176-189.
Gurpreet S. Bhamra, Anil K. Verma & Ram B. Patel, (2015) “Agent Based Frameworks for
Distributed Association Rule Mining: An Analysis”, International Journal in Foundations of
Computer Science & Technology (IJFCST), vol. 5, no. 1, pp. 11-22.

AUTHORS
Gurpreet Singh Bhamra is currently working as Assistant Professor at Department
of Computer Science and Engineering, M. M. University, Mullana, Haryana. He
received his B.Sc. (Computer Sc.) and MCA from Kurukshetra University,
Kurukshetra in 1995 and 1998, respectively. He is pursuing Ph.D. from Department
of Computer Science and Engineering, Thapar University, Patiala, Punjab. He is in
teaching since 1998. He has published 12 research papers in International/National
Journals and International Conferences. He is a Life Member of Computer Society of
India. His research interests are in Distributed Computing, Distributed Data Mining,
Mobile Agents and Bio-informatics.
Dr. Anil Kumar Verma is currently working as Associate Professor at Department
of Computer Science & Engineering, Thapar University, Patiala. He received his
B.S., M.S. and Ph.D. in 1991, 2001 and 2008 respectively, majoring in Computer
science and engineering. He has worked as Lecturer at M.M.M. Engineering College,
Gorakhpur from 1991 to 1996. He joined Thapar Institute of Engineering &
Technology in 1996 as a Systems Analyst in the Computer Centre and is presently
associated with the same Institute. He has been a visiting faculty to many
institutions. He has published over 100 papers in referred journals and conferences
(India and Abroad). He is a MISCI (Turkey), LMCSI (Mumbai), GMAIMA (New
Delhi). He is a certified software quality auditor by MoCIT, Govt. of India. His
research interests include wireless networks, routing algorithms and securing ad hoc networks and data
mining.

39

International Journal in Foundations of Computer Science & Technology (IJFCST), Vol.5, No.3, May 2015
Dr. Ram Bahadur Patel is currently working as Professor at Department of
Computer Science & Engineering, Chandigarh College of Engineering &
Technology, Chandigarh. He received PhD from IIT Roorkee in Computer Science &
Engineering, PDF from Highest Institute of Education, Science & Technology
(HIEST), Athens, Greece, MS (Software Systems) from BITS Pilani and B. E. in
Computer Engineering from M. M. M. Engineering College, Gorakhpur, UP. Dr.
Patel is in teaching and research since 1991. He has supervised 36 M. Tech, 7 M.
Phil. and 8 PhD Thesis. He is currently supervising 6 PhD students. He has published
130 research papers in International/National Journals and Conferences. He is
member of ISTE (New Delhi), IEEE (USA). His current research interests are in
Mobile & Distributed Computing, Mobile Agent Security and Fault Tolerance and Sensor Network.

APPENDIX A – DATASETS USED FOR AEQARM-AAPDB SYSTEM
A.1 PDB1 at site S1
This is a real dataset of protein sequences, the Astral Structural Classification of Proteins (SCOP)
[16], [17] version 1.75 genetic domain sequence subsets, based on PDB SEQRES records with
less than 40% identity to each other. SCOP database is comprehensive ordering of all proteins of
known structure, according to their evolutionary and structural relationships. Protein domains in
the SCOP are grouped into species and hierarchically classified into families, superfamilies, folds
and classes. In this dataset each protein record starts with greater than (`>') character followed by
protein description headers according to SCOP. An additional hash (`#') character is inserted as a
separator character between protein description headers and protein sequence of amino acids. A
total of 3523 protein records are stored in this protein data bank ( PDB1 ) at site S1 .

A.2 FPDB1 at site S1
PDB1 shown in Appendix A.1 is further filtered to generate FPDB1 for the protein sequence
length range ≥ 50 and < 400 . A total of 3341 such filtered protein records are obtained.

40

International Journal in Foundations of Computer Science & Technology (IJFCST), Vol.5, No.3, May 2015

A.3 AAF1 generated from FPDB1 at site S1
Each column heading represents single letter code of an amino acid. Each row represents the
frequencies of amino acids in a particular protein sequence record.

A.4 Amino acids frequency ranges (15 ranges for each amino acid)
Item No. column represents Serial No. from 1 to 300. Data in the AA Frequency Range column
represents single letter code of amino acid along with a frequency range. Frequency of each
amino acid is divided into 15 partitions (ranges) resulting into 300 items for 20 amino acids. For
example if frequency of amino acid Alanine (A) in a protein record is 22 it lies at Item No. 8.
41

International Journal in Foundations of Computer Science & Technology (IJFCST), Vol.5, No.3, May 2015

A.5 BDB1 at site S1
This Boolean data bank ( BDB1 ) is created using AAF1 as shown in Appendix A.3 and amino
acids frequency range table as shown in Appendix A.4. Each column heading represents an Item
No. as shown in Appendix A.4. Each amino acid has 15 frequency partitions and for each amino
acid, a boolean value ‘1’ is put under that Item No. in which frequency of amino acid lies in a
particular protein record otherwise a value ‘0’ is considered. For example it clear from AAF1 that
the frequency of amino acid Alanine (A) in the 1st protein record is 22 and this frequency lies at
Item No. 8 in frequency range table as shown in Appendix A.4, so a boolean value ‘1’ is put
under Item No. 8 in the 1st protein record in BDB1 .

42

International Journal in Foundations of Computer Science & Technology (IJFCST), Vol.5, No.3, May 2015

A.6 IDB1 at site S1
This Itemset data bank ( IDB1 ) is created using BDB1 as shown in Appendix A.4. Each record in
this dataset consist of 20 Item Numbers, one for each of amino acids for which a boolean value
‘1’ is stored in BDB1 .

APPENDIX B – RESULTANT KNOWLEDGE OF AEQARM-AAPDB SYSTEM
FISC
B.1 LFI
k (1) and Lk (1) at site S1

List of frequent k-itemset, i.e., LFI
k (1) is represented by column L and column SC shows the
support count of the corresponding frequent k-itemset, i.e., LFISC
at site S1 . These frequent
k (1)
itemsets and their support counts are obtained by processing the Itemset Data Bank ( IDB1 ) as
shown in Appendix A.6.

43

International Journal in Foundations of Computer Science & Technology (IJFCST), Vol.5, No.3, May 2015

B.2 L1LSAR for Item numbers at site S1
Column L represents frequent k-itemset and column AR(support,confidence) shows the list of
locally strong association rules, i.e., L1LSAR at site S1 . Each strong rule has its associated support
and confidence factor. The minimum threshold support is taken as 20% and minimum threshold
confidence as 50% for generating the strong rules by making use of the data as shown in
Appendix B.1.

44

International Journal in Foundations of Computer Science & Technology (IJFCST), Vol.5, No.3, May 2015

B.3 L1LSAR for corresponding Amino acid frequency ranges at site S1
Replace Item No. in Appendix B.2 data with its corresponding Amino acid frequency range as
shown in table Appendix A.4.

45

