the chemokine proteins are a class of small molecules that play a significant role in leukocyte trafficking during immune response ( 1 ) .
ccr1 , as one of the chemokine receptors , is expressed on a number of human cells such as monocytes , macrophages , dendritic cells , and t cells ( 2).a large number of studies have provided strong evidences for a significant role of the chemokines , rantes ( regulated upon activation , normal t cell expressed and presumably secreted ) , and mip-1a ( macrophage inflammatory protein-1 ) in chronic inflammatory diseases . because mip-1a and rantes are agonists for ccr1 , antagonists for this protein may be helpful in treatment of these diseases ( 3 ) .
qsar has become a very well known discipline in the drug discovery researches ( 46 ) . the basis of such relationships is the assumption that the variation of bioactivity of molecules , as expressed by pic50 , can be regressed with changes in molecular descriptors .
development of qsar involves selection of most informative independent variables to describe different sets of molecules and the application of various algorithms , such as multiple linear regression or ann to construct the qsar model .
the advantage of ann is it inherent power to reveal non - linear relationships between idependent and dependent variables in the derivation of the qsar models .
a quantitative structure activity relationship ( qsar ) model based on artificial neural networks ( ann ) was developed to study the activities of 29 derivatives of 3-amino-4-(2-(2-(4-benzylpiperazin-1-yl)-2-oxoethoxy ) phenylamino ) cyclobutenedione as c - c chemokine receptor type 1(ccr1 ) inhibitors ( 7 ) .
biological and chemical data from 29 derivatives of 4 - 3-amino-4-(2-(2-(4-benzylpiperazin-1-yl)-2-oxoethoxy)phenylamino)cyclobutenedione were used in this study ( table 1 ) ( 7 ) .
all the structures were drawn and optimized using the semiempirical quantum - chemical routine of am1 implemented in hyperchem ( 8) . 
 
the structures and biological activities of compounds 
  the structures and biological activities of compounds 913 
  the structures and biological activities of compounds 1427 
  the structures and biological activities of compounds 2829 
  compounds selected as test set in order to model the biological activities of the studied compounds , four classes of descriptors were calculated : constitutional , geometrical , topological , and functional group using dragon ( 9 ) on the minimal energy conformations .
the data set was devided into the training and the testing sets based on kennard and stone algorithm ( 10 ) .
principal component analysis was used to compress a pool of descriptors into principal components ( pcs ) as new variables .
after that , a model using a nonlinear regression model , artificial neural network , was constructed to make a relationship between pcs and the pic50 .
a feed - forward ann with error back - propagation learning algorithm was applied for model building .
pca was performed on the data set that gives 14 significant pcs ( % variance explaind>1 ) .
therefore , the next steps of study were restricted to these 14 pcs . plotting of first pc vs. second pc showed none of the compounds is outlier ( fig .
clarification of the theory of the artificial neural networks in details has been adequately described elsewhere ( 11 ) and some relevant remarks is presented . 
 
the first two components ( pc1 , and pc2 ) from the principal component analysis of the 29 studied molecules .
the first layer namely input layer has ni neurons , and function of this layer is reception of information ( i.e. inputs ) which transfers them to all neurons in the next layer called the hidden layer that their number are indicated by nh .
the neurons in the hidden layer calculate a weighted sum of the inputs that is subsequently transformed by a linear or non - linear function .
the last layer is the output layer and its neurons handle the output from the network and it calculate response vector .
the function of synapses is connection of input layer to hidden layer and hidden layer to output layer .
the manner in which each node transforms its input depends on the  weights  and bias of the node , which are modifiable .
ba - ann network was trained with the training set of molecules using a back propagation algorithm followed by conjugate gradient descent in the second stage ( 47 ) .
rmsecv was then employed as tool to select optimum value of various parameters ( 6 ) .
figure 2 shows the effect of the different number of pcs on predictability of developed model . on the based of this figure , the ann has the highest degree of predictability when the number of pcs is 4 .
the parameters of network which should be optimized are learning rate , number of neuron in hidden layer , momentum.the optimal values for these parameters as it is shown in figures 3a - c are 1.2 , 9 , and 0.9 , respectively .
optimization of number of neurons in hidden layer ( a ) , momentum ( b ) , and , learning rate ( c ) .
the predicted activity of the ann calculated values of pic50 versus the experimental values are shown in figure 4 and reported in table 1 .
as it was expected , the calculated values are in good agreement with the experimental values . calculated vs. experimental activity of the investigated compounds in training and test sets 
various statistical criteria for ann model were calculated and reported in table 3 ( 58 ) .
the external predictability of a proposed model was generally tested using test sets and rcv2 
 .
the satisfactory prediction of the values of the inhibitory activity of the test set compounds demonstrates the efficacy of the qsar in predicting the activities of external molecules . 
 
various statistical parameter for developed pc- ann model 
 r 
 = square regression coefficient 
 rmse = root mean square error 
 press = predicted error sum of square for training set 
 
 k=yiyiyi2fk'=yiyiyi2 

pca was performed on the data set that gives 14 significant pcs ( % variance explaind>1 ) .
therefore , the next steps of study were restricted to these 14 pcs . plotting of first pc vs. second pc showed none of the compounds is outlier ( fig .
clarification of the theory of the artificial neural networks in details has been adequately described elsewhere ( 11 ) and some relevant remarks is presented . 
 
the first two components ( pc1 , and pc2 ) from the principal component analysis of the 29 studied molecules .
the first layer namely input layer has ni neurons , and function of this layer is reception of information ( i.e. inputs ) which transfers them to all neurons in the next layer called the hidden layer that their number are indicated by nh .
the neurons in the hidden layer calculate a weighted sum of the inputs that is subsequently transformed by a linear or non - linear function .
the last layer is the output layer and its neurons handle the output from the network and it calculate response vector .
the function of synapses is connection of input layer to hidden layer and hidden layer to output layer .
the manner in which each node transforms its input depends on the  weights  and bias of the node , which are modifiable .
ba - ann network was trained with the training set of molecules using a back propagation algorithm followed by conjugate gradient descent in the second stage ( 47 ) .
rmsecv was then employed as tool to select optimum value of various parameters ( 6 ) .
figure 2 shows the effect of the different number of pcs on predictability of developed model . on the based of this figure , the ann has the highest degree of predictability when the number of pcs is 4 .
the parameters of network which should be optimized are learning rate , number of neuron in hidden layer , momentum.the optimal values for these parameters as it is shown in figures 3a - c are 1.2 , 9 , and 0.9 , respectively .
optimization of number of neurons in hidden layer ( a ) , momentum ( b ) , and , learning rate ( c ) .
the predicted activity of the ann calculated values of pic50 versus the experimental values are shown in figure 4 and reported in table 1 .
as it was expected , the calculated values are in good agreement with the experimental values . calculated vs. experimental activity of the investigated compounds in training and test sets 
various statistical criteria for ann model were calculated and reported in table 3 ( 58 ) .
the external predictability of a proposed model was generally tested using test sets and rcv2 
 .
the satisfactory prediction of the values of the inhibitory activity of the test set compounds demonstrates the efficacy of the qsar in predicting the activities of external molecules . 
 
various statistical parameter for developed pc- ann model 
 r 
 = square regression coefficient 
 rmse = root mean square error 
 press = predicted error sum of square for training set 
 
 k=yiyiyi2fk'=yiyiyi2 

pca was performed on the data set that gives 14 significant pcs ( % variance explaind>1 ) .
therefore , the next steps of study were restricted to these 14 pcs . plotting of first pc vs. second pc showed none of the compounds is outlier ( fig .
clarification of the theory of the artificial neural networks in details has been adequately described elsewhere ( 11 ) and some relevant remarks is presented . 
 
the first two components ( pc1 , and pc2 ) from the principal component analysis of the 29 studied molecules .
the first layer namely input layer has ni neurons , and function of this layer is reception of information ( i.e. inputs ) which transfers them to all neurons in the next layer called the hidden layer that their number are indicated by nh .
the neurons in the hidden layer calculate a weighted sum of the inputs that is subsequently transformed by a linear or non - linear function .
the last layer is the output layer and its neurons handle the output from the network and it calculate response vector .
the function of synapses is connection of input layer to hidden layer and hidden layer to output layer .
the manner in which each node transforms its input depends on the  weights  and bias of the node , which are modifiable .
ba - ann network was trained with the training set of molecules using a back propagation algorithm followed by conjugate gradient descent in the second stage ( 47 ) .
rmsecv was then employed as tool to select optimum value of various parameters ( 6 ) .
figure 2 shows the effect of the different number of pcs on predictability of developed model . on the based of this figure , the ann has the highest degree of predictability when the number of pcs is 4 .
the parameters of network which should be optimized are learning rate , number of neuron in hidden layer , momentum.the optimal values for these parameters as it is shown in figures 3a - c are 1.2 , 9 , and 0.9 , respectively .
optimization of number of neurons in hidden layer ( a ) , momentum ( b ) , and , learning rate ( c ) .
the predicted activity of the ann calculated values of pic50 versus the experimental values are shown in figure 4 and reported in table 1 .
as it was expected , the calculated values are in good agreement with the experimental values . calculated vs. experimental activity of the investigated compounds in training and test sets 
various statistical criteria for ann model were calculated and reported in table 3 ( 58 ) .
the external predictability of a proposed model was generally tested using test sets and rcv2 
 .
the satisfactory prediction of the values of the inhibitory activity of the test set compounds demonstrates the efficacy of the qsar in predicting the activities of external molecules . 
 
various statistical parameter for developed pc- ann model 
 r 
 = square regression coefficient 
 rmse = root mean square error 
 press = predicted error sum of square for training set 
 
 k=yiyiyi2fk'=yiyiyi2 

in the developed model , a network including a fully connected three layer , feed - forward ann model trained with a back - propagation learning algorithm was used .
the input of the network was the eigenvalue ranked pcs , the number of them which entered neural network varied from 1 to 14 , of which 4 pcs of them were selected as input of networks . by using this number of pcs the best results on the basis of lowest rmsecv in the output of network were obtained . since there are no exact theoretical principles for choosing the appropriate network topology , before the training of the network , the adjustable parameters as number of nodes in the hidden layer , transfer function , learning rate and etc .
the values resulting from hidden layer are transferred to the last layer , which contains a single neuron representing the predicted activity .
various ann architectures were run with the four selected pcs as input . in each run , the neuron architecture and parameters were optimized to reach the lowest rmsecv as the performances of the resulted models . according to the criteria proposed by tropsha and roy ( 46 ) , for testing the reliability and the robustness of qsar models ,
as discussed above , the calculated pcs have meaning physicochemically , but they may be employed for building statistical models which help the medicinal chemist limit the number of compounds to be synthesized .
for instance , medicinal chemist can propose a training set comprised of molecules which have the characters of two or more chemical classes with the smallest amount of similarity .
therefore , the qsar model was used to estimate inhibitory activities of a few suggested compounds .
the general structures of four suggested compounds and also their calculated activities are reported in table 4 .
the relative high predicted activity of the tested compounds suggest further study such as synthesis of other compounds with such chemical structures . 
 
the main objective of this study was to define and establish a qsar model to predict bioactivity of a series of 3-amino-4-(2-(2-(4-benzylpiperazin-1-yl)-2-oxoethoxy ) phenylamino ) cyclobutenedione derivatives as novel ccr1 antagonists without any knowledge of the under study system .
calculated pcs were used to make model of the relationship between the molecule structures of the studied compounds and the corresponding bioactivities .
the study showed that the calculated pcs as input variable to network can improve the predictive ability of the neural networks .
moreover , the suggested qsar model was based on nonlinear ann approach , which can be employed to simulate any kinds of complex correlation or function relationship in a given multivariable system .
i.e. , ann approach is more appropriate for modeling where no clearly defined mathematical model for a system is available .
therefore , accurate , well - organized and intelligent qsar model for the bioactivity will be influential for drug design and development .
 background and purpose of the studya quantitative structure activity relationship ( qsar ) model based on artificial neural networks ( ann ) was developed to study the activities of 29 derivatives of 3-amino-4-(2-(2-(4-benzylpiperazin-1-yl)-2-oxoethoxy ) phenylamino ) cyclobutenedione as c - c chemokine receptor type 1(ccr1 ) inhibitors.methodsa feed - forward ann with error back - propagation learning algorithm was used for model building which was achieved by optimizing initial learning rate , learning momentum , epoch and the number of hidden neurons.resultsgood results were obtained with a root mean square error ( rmse ) and correlation coefficients ( r 
 2 ) of 0.189 and 0.906 for the training and 0.103 and 0.932 prediction sets , respectively.conclusionthe results reflect a nonlinear relationship between the principal components obtained from calculated molecular descriptors and the inhibitory activities of the investigated molecules . 
INTRODUCTION
METHODS
RESULT AND DISCUSSION
None
Results
Discussion
CONCLUSION