Printer Friendly

QSAR Study of flavonoids as platelet-activating factor (PAF) receptor binding antagonists.

INTRODUCTION

Platelet activating factor (PAF) (1-0-hexadecyl/octadecyl-2-acetyl-sn-glycero-3-phosphocholine) [1] is a potent phospholipid mediator. This protein is formed by several different cells such as eosinophils, macrophages. platelets, neutrophils and vascular endothelial cells. Its function was to transmit signals between neighboring cells and acts as a hormone, cytokines, and other signaling molecules [2]. The PAF unregulated signaling system however, can trigger inflammatory and thrombotic cascades, amplify these cascades when acting with other mediators, and mediate molecular and cellular interactions (cross talk) between inflammation and thrombosis.

This biologically active phospholipid involved in several pathophysiological condition such as inflammation, atherosclerosis [3] asthma [4] and thrombosis [5]. PAF is also involved in tissue injury during endotoxemia [1] and also participate in several physiological responses including attraction and activation of leucocytes [6], increased in vascular permeability [7], cardiovascular actions, [8], effects in liver [9] and the central nervous system. Platelets could be activated by stimuli such as thrombin to synthesis and secrete PAF[3]. The attachment of PAF to it receptor form a series of intracellular signals to finally cause pathophysiological and physiological processes. The PAF actions could be prevented by blocking its receptor with an antagonist. Thus, it is important to find a compound(s) that could act as the PAF receptor binding antagonist to prevent a variety of inflammation, respiratory, immunological and cardiovascular disorder [10].

One of the natural PAF antagonistsisginkgolides from Ginkobiloba [11]. Three quantitative structure activity relationship (QSAR) approachesi.eCoMFA, CoMSIA, and HQSAR, were used to investigate the relationship between 117 ginkgolide analogues and their bioactivities against PAF receptor. The models revealed that steric, electrostatic, hydrophobicity, and individual atoms affect molecular bioactivity as antagonists of PAF.

However, ginkgolides has a cage skeleton consisting of six five-membered rings, therefore, are very tough to be synthesized and they have their specific mode of action against selected pathophysiological condition cause by PAF.

Flavonoids, another natural inhibitor against PAF receptor [12] could be found in many plant species and have simpler skeleton, which is easy to synthesize and modify. Hence an investigation of the relationship between the flavonoids analogues structures and their activity as the PAF receptor binding antagonist, embarked in this research, so that the understanding derived from this study can be used to predict and synthesize highly active flavonoids.

MATERIALS AND METHODS

The structure of flavonoids 1-11 and their PAF receptor binding inhibitory activity were obtained from literature [12]. The data were expressed as [IC.sub.50]([micro]M)and converted into log [IC.sub.50] in order to improve the normal distribution of the experimental data points.

The structures (Figure 1) were drawn using Chem3D ultra 7 software. Their initial geometry were optimized using molecular mechanics methods and the resultant structures were subjected to quantum mechanical self-consistent field (SCF) calculations using MOPAC software [13] and RM1 Semi empirical quantum chemical method [14]. Gradient norm limit of 0.0 kcal/A was used as stopping criteria for optimized structures. The optimized geometries along with the molecular orbital and vibrational frequency information were be used to calculate over 400 molecular descriptors containing constitutional, topological, geometrical, charge-related, semi empirical and thermo-dynamical descriptors using CODESSA software [16]. The CODESSA software v2.6 were used to calculate descriptors from the produced electronic structural data and finally to get linear correlations between the calculated and experimental PAF receptor binding inhibitory activity values using either best multi linear (BMLR) or Heuristic method.

Among the several possible QSAR models, the best QSAR could be determined from several statistical criteria such as [R.sup.2] and the F value. Cross validation is one of the model validation techniques to internally check the on the derived model. It gives the estimation of the predictive power of the model to predict the values of the untested compounds. It was performed in this study using the leave-one-out method where one compound is removed from the data set and its activity is predicted using the model derived from the rest of the data set [18]. The corresponding squared cross-validated correlation coefficient ([R.sup.2.sub.CV]) for the selected model is calculated automatically by the validation module implemented in CODESSA software v 2.6 [16].

RESULTS AND DISCUSSION

QSA Rmodeling:

From the Heuristic method, thirteen linear calculations were obtained (Table 1), where as the best multi linear calculation gave five linear correlations (Table 2). Over fitting phenomenon occurs when exceedingly large number of descriptors was used in the generation of the best multi linear model for a given property [15]. The over fitting was avoided by the procedure called break point technique as illustrated in Fig. 2 that shows the change in slope in the plot of [R.sup.2] versus the number of descriptors added [16]. If the addition of a new descriptor leads to only a small (<0.02) increase in the correlation coefficient, then the added descriptor is adding very little information to the QSAR model. The procedure was stopped when the difference between [R.sup.2] of the two consequent regression equations was less than or equal 0.02 [16]. This method managed to avoid over fitting of the regression equations by monitoring the increase of [R.sup.2] in the equations with successive number of descriptors involved. Hence, the three descriptors correlation obtained using best multi linear method was selected as the best QSAR model and the parameters for this equation (Eq. 1) is given in Table 3.

The best QSAR model with the parameter could be expressed by the equation:

Log [IC.sub.50] = 0.18592 {td)-0.15194 (tm)+ 6.9251 (me)-61.354 (1)

The correlations between the calculated and experimental activity according to the above equation is shown in Fig. 3. The experimental and predicted log [IC.sub.50] of the flavonoids 1-11 is presented in table 4. The result shows small differences between the experimental and predicted log [IC.sub.50] values (less than 5%) implying good predictive power of the QSAR model.

This model has the correlation coefficient ([R.SUP.2]) value of 0.9632, cross-validated correlation coefficient ([R.sup.2.sub.CV]) value of 0.9279, squared standard deviation (S2) value of 0.0012 and Fisher criterion (F) value of 61.0291. The close agreement between [R.sub.2] and [R.sup.2.sub.CV] shows the good predictive power of the QSAR model.

The type of descriptors involved in the equation 1 are total dipole of the molecule (td) (electronic descriptor), total molecular 2-center resonance energy divided by the no of atoms (tm) (quantum chemical descriptor) and min e-n attraction for H atom(me) (electronic descriptor) (Table 3). The equation showed that the increasing in total dipole of the molecule and min e-n attraction for a H atom as well as the decreasing in total molecular 2-center resonance energy/no of atoms would increase the [IC.sub.50] value of PAF receptor antagonist activity which showed the reducing in the activity.

The total dipole of the molecule (td) is given by the following formula [16]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2)

Where [empty set] are the molecular orbitals, r are the electron position operator; [Z.sub.a] are the a-th atomic nuclear charge and [[??].sub.a] is the position vector of a-th atomic nucleus. Based on the QSAR model (Eq. 1), the log [IC.sub.50] value is found to be directly proportion to the dipole moment and thus implying that the inhibition activity will drop if the polar functional groups, such as hydroxyl groups, are introduced. This explains why compounds 4 and 7 which contains no hydroxyl groups are more active than other compounds containing hydroxyl groups. The resonance energy between given two atomic species can be expressed as [16]:

[E.sub.R](AB) = [[summation].sub.[mu][member of]] [[summation].sub.Av[member of]B] [P.sub.[mu]v] [[beta].sub.[mu]v] (3)

Where A s the given atomic species; B another atoms; [P.sub.[mu]v] density matrix elements over atomic basis {[micro]v}; [[beta].sub.[micro]v] is resonance integrals on atomic basis{[micro]v}.

This resonance energy is much related to the benzene ring in the flavonoids structure, thus indicates the importance of resonance in the activity as the antagonist of PAF receptor binding.

The min nuclear-electron attraction energy for hydrogen atom is defined as [16]:

[E.sub.ne] (AB) = [[summation].sub.B] [[summation].sub.[micro]v[member of]A] [P.sub.[micro]v]([mu][absolute value of [Z.sub.B]/[R.sub.iB]]v) (4)

The first summation is performed over all atomic nuclei in the flavonoid molecule (B) whereas the second

summation is carried out over all atomic orbitals at a given atom, i.e. hydrogen (A) [16]. The term,([mu][absolute value of [Z.sub.B]/[R.sub.iB]]), represent the nuclear-electron attraction integrals on the given atomic basis. This energy may be related to the conformational such as rotational inversional changes in the molecule. The presence of this term in the QSAR equation emphasizes the importance of conformation in predicting the activity of the molecule.

Method validation:

It is important to determine the stability and predictive power of a predictive model such as QSAR model. If the data set is very large, then the data set can be divided into training set which will be used for obtaining the QSAR and tested against the external (test) set. However when the data set is smaller, the internal validation set such as leave-one-out method could be used to evaluate the predictive power of the developed QSAR.

The leave-one-out validation method implemented in the CODESSA software gave the corresponding cross-validated correlation coefficients ([R.sup.2.sub.CV]) to all models. The cross-validated correlation coefficient ([R.sup.2.sub.CV]) for the best three descriptor QSAR equation (Eq. 1) was 0.9279, which is very close to the [R.sup.2] (0.9632). The small difference between the [R.sup.2] and [R.sup.2.sub.CV] suggest the good predictive ability of the model.

Conclusions:

The quantum chemical quantitative structure activity relationship approach was applied to study the correlation between the PAF receptor binding antagonist activity and the chemical structure of flavonoids 1-11. This approach is very powerful, and can be applied without the necessity of measuring any additional physicchemical properties of the active molecules. The polarity of the substituent functional groups, conformation of the active molecules and the resonance energies are found to play significant roles in determining the activity of the molecules. The obtained QSAR was found to have good predictive power, and hence could be used to determine or construct new structure with better inhibition activity.

ARTICLE INFO

Article history:

Received 14 Feb 2014

Received in revised form 24 February 2014

Accepted 29 March 2014

Available online 14 April 2014

REFERENCES

[1] Fernando, L.D., L.N. Maria, F.G. Sagrario, A.G. Miguel and S.C. Mariano, 1989. J. Clin. Invest., 83: 1733-1740.

[2] Prescott, S.M., G.A. Zimmberman, D.M. Stafforini and T.M. McIntyre, 2000. Annu. Rev. Biochem., 69: 419-45.

[3] John, V.M., P.V. Maria, S. Dominique, N. Ewa and D.T. Alexandras, 2006. Arterioscler. Thromb. Vasc Biol., 20: 80-88.

[4] Zimmerman, G.A., T.M. McIntyre, S.M. Prescott and D.M. Stafforini, 2002. Crit, Care Med., 30: A294-S301.

[5] Zhang, J., C. Li, X. Guo and G.J. Wang, 2001. Tradit. Chin. Med., 21: 299-302.

[6] Camussi, A., P. Landi, M. Bertolini, 1988. nalysis of variety crosses to develop early base population for reciprocal recurrent selection in maize.Maydica, 33: 269-281.

[7] Stimler, N.P., C.M. Bloor, T.E. Hugli, R.L. Wykle, C.E. McCall, J.T. O'Flaherty, 1981. Am. J. Pathol., 105: 64-69.

[8] Montrucchio, G., G. Alloati and G. Camussi, 2000. Physiol. Rev., 80: 1669-1690.

[9] Sakaguchi, T., S. Nakamura, S. Suzuki, T. Oda, A. Lahiyama and T. Okamoto, 1999. Hepatology, 30: 959-967.

[10] Moharam, B.A., I. Jantan, J. Jalil and K. Shaari, 2012. Molecules, 15: 7840-7848.

[11] Zhu, W., G. Chen, L. Hu, X. Luo, C. Gui, C. Luo, C.M. Puah, K. Chen and H. Jiang, 2005. Bioorg Med. Chem., 13(2): 313-22.

[12] Ling, S.K., M.P. Mazura and M. Salbiah, 2007. Biol. Pharm. Bull., 30(6): 1150-1152.

[13] Stewart, J.J.P., MOPAC 6.0, 1990. Frank, J. Seiler Research Laboratory, United States Air Force Academy, Colorado, USA.

[14] Rocha, G.B., R.O. Freire, A.M. Simas and J.J.P. Stewart, 2006. Comput. Chem., 27: 1101-111.

[15] Jurgen, B., 2004. Chemoinformatics: Concepts, methods and tools for Drug Discovery, Springer International publisher science, USA.

[16] Kartrizky, A.R., L.M. Pacureanu, S. Slavov, D.A. Dobchev and M. Karelson, 2006. Biooorganic and medicinal chemistry, 14: 7490-7500.

[17] Karelson, M., 2000. Molecular descriptors in QSAR/QSPR, Wiley-Interscience.

[18] CODESSA PRO user manual retrieved from www.codessa-pro.com/manuals/manual.htmon 12 September 2013.

(1) N.H. Abdullah, (1) L.S. Kiong, (2) I.A. Noorbatcha

(1) Natural Product Division, Forest Research Institute Malaysia, 52109 Kepong, Selangor, Malaysia

(2) Department of Biotechnology Engineering, Faculty of Engineering, International Islamic University Malaysia, Jalan Gombak 53100 Kuala Lumpur, Malaysia

Corresponding Author: Ibrahim Ali Noorbatcha, Department of Biotechnology Engineering, Faculty of Engineering, International Islamic University Malaysia, Jalan Gombak 53100 Kuala Lumpur, Malaysia. E-mail: ibrahiman@iium.edu.my. Tel: 03-6196 5453.

Table 1: The statistical parameters for the QSAR models
selected using Heuristic method for flavonoids 1-11.

                                                        No. of
 No.    [R.sup.2]      F      [S.sup.2]   [R.sup.2]   descriptor

  1      0.9777     43.8808    0.0013      0.9340         5
  2      0.9776     43.6726    0.0013      0.9236         5
  3      0.9771     42.6824    0.0014      0.9068         5
  4      0.9769     42.2275    0.0014      0.9300         5
  5      0.9763     41.1531    0.0014      0.9173         5
  6      0.9762     41.0771    0.0014      0.9398         5
  7      0.9760     40.6914    0.0014      0.9070         5
  8      0.9759     40.4520    0.0014      0.8982         5
  9      0.9756     40.0457    0.0015      0.8912         5
 10      0.9748     38.6960    0.0015      0.9078         5
 11      0.9671     44.0566    0.0016      0.9102         4
 12      0.9482     42.7007    0.0022      0.8853         3
 13      0.9644     40.6386    0.0018      0.9018         4

Table 2: The statistical parameters for the QSAR models selected
using BML method for flavonoids (1-11).

No.   [R.sup.2]      F      [S.sup.2]    R2cv     No. of
                                                 descriptor

1      0.8233     18.6362    0.0015     0.6650       2
2      0.9632     61.0291    0.0012     0.9279       3
3      0.9761     61.3906    0.0012     0.9416       4
4      0.9792     47.0306    0.0016     0.9356       5
5      0.9800     32.7146    0.0066     0.4544       6

Table 3: The coefficients, statistical errors and the t-
test values for thebest nonlinear three descriptorQSAR
models selected using BML method for flavonoids 1-11.

                    [+ or -]
S/N        X         [DELTA]X    t-test           Descriptor

0     -6.1354e+01   7.3171e+00   -8.3850   Intercept

1     1.8592e-01    1.5696e-02   11.8450   Total dipole of
                                             the molecule (td)
2     -1.5194e-01   3.6146e-02   -4.2035   Total molecular 2-center
                                             resonance energy/
                                             no of atoms (tm)
3     6.9251e+00    8.6434e-01   8.0120    Minimum e-n attraction
                                             for a H atom (me)

The best QSAR model with the parameter could be expressed by the
equation:

Log [IC.sub.50] = 0.18592 (td)-0.15194 (tm)+ 6.9251(me) - 61.354    (1)

Table 4: The calculated and experimental values of the best
correlation for flavonoids 1-11.

Compound         R1              R2       R3       R4

1                OH              H        H    C[H.sub.3]
2                 H          C[H.sub.3]   OH   C[H.sub.3]
3                OH          C[H.sub.3]   H    C[H.sub.3]
4                 H              H        H    C[H.sub.3]
5                OH              H        H    C[H.sub.3]
6          O-glu (6 [less        H        OH       H
            than or equal
              to] 1)rham
7          O-glu (6 [less        H        H        H
            than or equal
              to] 1)rham
8                OH              H        OH   C[H.sub.3]
9                OH          C[H.sub.3]   OH       H
10               OH              H        OH   C[H.sub.3]
11               OH          C[H.sub.3]   OH       H

Compound   Experimental    Predicted    Difference
           [IC.sub.50]    [IC.sub.50]

1             1.9000        1.8434      0.0566
2             1.6900        1.6796      0.0104
3             1.6800        1.6724      0.0076
4             1.4400        1.4253      0.0147
5             1.7900        1.8434      0.0534
6             1.6800        1.6854      0.0054

7             1.2900        1.2937      0.0037

8             1.6000        1.6571      0.0571
9             1.5900        1.5733      0.0167
10            1.7200        1.7345      0.0145
11            1.8200        1.7920      0.0280
COPYRIGHT 2014 American-Eurasian Network for Scientific Information
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2014 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Author:Abdullah, N.H.; Kiong, L.S.; Noorbatcha, I.A.
Publication:Advances in Environmental Biology
Article Type:Report
Geographic Code:9MALA
Date:Mar 1, 2014
Words:2774
Previous Article:Characterizations of thermoplastic starch and starch reinforced montmorillonite clay nanocomposite.
Next Article:Screening of novel acidified solvents for maximal antimicrobial peptide extraction from Zophobas morio fabricius.
Topics:

Terms of use | Privacy policy | Copyright © 2020 Farlex, Inc. | Feedback | For webmasters