# Development of an in silico model of DPPH• free radical scavenging capacity: prediction of antioxidant activity of a coumarin type compounds

Elizabeth Goya Jorge, Anita Maria Rayar, Stephen Jones Barigye, María Elisa Jorge Rodríguez, Maité Sylla-Iyarreta Veitía

Elisa Jorge Rodríguez, Maité Sylla-Iyarreta Veitía

► To cite this version:

Elizabeth Goya Jorge, Anita Maria Rayar, Stephen Jones Barigye, María Elisa Jorge Rodríguez, Maité Sylla-Iyarreta Veitía. Development of an in silico model of DPPH• free radical scavenging capacity: prediction of antioxidant activity of a coumarin type compounds. International Journal of Molecular Sciences, 2016, 17 (6), pp.1-10. 10.3390/ijms17060881. hal-01487390

HAL Id: hal-01487390

https://hal.science/hal-01487390

Submitted on 4 Sep 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Development of an *in Silico* Model of DPPH• Free Radical Scavenging Capacity: Prediction of Antioxidant Activity of Coumarin Type Compounds

Elizabeth Goya Jorge<fnref n="1" />, Anita Maria Rayar<fnref n="2" />, Stephen J. Barigye<fnref n="3" />, María Elisa Jorge Rodríguez<fnref n="1" /> and Maité Sylla-Iyarreta Veitía<fnref n="2" />*

Academic Editor: Humberto González-Díaz

Received: 24 March 2016; Accepted: 25 May 2016; Published: 7 June 2016

**Abstract:** A quantitative structure-activity relationship (QSAR) study of the 2,2-diphenyl-1-picrylhydrazyl (DPPH•) radical scavenging ability of 1373 chemical compounds, using DRAGON molecular descriptors (MD) and the neural network technique, a technique based on the multilayer multilayer perceptron (MLP), was developed. The built model demonstrated a satisfactory performance for the training (R<fnref n="2" /> = 0.713) and test set (Q²ext = 0.654), respectively. To gain greater insight on the relevance of the MD contained in the MLP model, sensitivity and principal component analyses were performed. Moreover, structural and mechanistic interpretation was carried out to comprehend the relationship of the variables in the model with the modeled property. The constructed MLP model was employed to predict the radical scavenging ability for a group of coumarin-type compounds. Finally, in order to validate the model's predictions, an *in vitro* assay for one of the compounds (4-hydroxycoumarin) was performed, showing a satisfactory proximity between the experimental and predicted pIC₅₀ values.

**Keywords:** artificial neural networks; MLP; antioxidant; QSAR; DPPH•; free radical scavenger; coumarin

## 1. Introduction

The discovery of antioxidant agents has attracted much attention in recent years, because oxidative damage is related to many pathological conditions [1]. Compounds whose antioxidant activity is based on scavenge free radicals are among the most important and studied antioxidants. Numerous *in vivo* and *in vitro* methods have been developed to measure the antioxidant capacity and effectiveness, but there is no universal technique by which this property can be estimated precisely and quantitatively [2–5]. The *in vitro* scavenging capacity can be evaluating using biological oxidants such as nitric oxide radical (NO•) and hydroxyl radical (OH•), or with nonbiological oxidants such as 2,2'-azino-bis-(3-ethylbenzothiazoline-6-sulphonate) radical cation (ABTS) and 2,2-diphenyl-1-picrylhydrazyl radical (DPPH•). The capturing of the DPPH• radical is one of the best-known methods [2,6,7].

A great number of coumarin derivatives have been studied for their biochemical and pharmacological profiles. Some studies suggest that these compounds may significantly affect the function of various mammalian cellular systems. Specifically, their antioxidant effect has been

<footnote n="1">¹ Pharmacy Department, Faculty of Chemistry and Pharmacy, Central University “Marta Abreu” of Las Villas, C-54830 Santa Clara, Cuba; egoyaj@gmail.com or egoyaj@yandex.com (E.G.J.); elisajorge@yahoo.es (M.E.J.R.)</footnote>

<footnote n="2">² Equipe de Chimie Moléculaire du Laboratoire CMGPCE, EA 7341, Conservatoire National des Arts et Métiers, 2 rue Conté, 75003 Paris, France; anitarayar@hotmail.fr</footnote>

<footnote n="3">³ Department of Chemistry, Federal University of Lavras, P.O. Box 3037, 37200-000 Lavras, Brazil; sjbarigye@gmail.com</footnote>

<footnote>* Correspondence: maite.sylla@lecnam.net; Tel.: +33-1-58-80-84-82</footnote>
explored, since the structural features of this group of compounds suggest that they can exhibit this pharmacologic property [8–10].

Chemoinformatics tools are used in the modeling of the antiradical activity, as well as other biological properties, given their advantages in saving time and resources [11,12]. Complex associations between biological activities and chemical features of compounds have been approached with mathematical models of quantitative structure-activity relationships (QSAR). [13]. Several statistical and machine learning methods have been widely used in the literature to build models for the QSAR study of chemical compounds. For the last two decades, artificial neural network (ANN) techniques have increasingly found applicability in QSAR studies, thanks to their ability to map non-linear relations between structural characteristics of chemical compounds and their chemical/biological behavior [14].

The objective of this study was to develop an ANN model in order to relate the chemical compounds' scavenging ability of the DPPH• radical with the corresponding structural features, also known as molecular descriptors (MD), in a set of 1373 molecules. Then, the built model was used to predict the antioxidant activity of a group of coumarin derivatives and an *in vitro* study of 4-hydroxycoumarin was performed to corroborate the result predicted by the network.

## 2. Results

### 2.1. Modeling

**Compounds:** The whole dataset of 1373 chemical compounds was optimized for three-dimensional (3D) coordinates with CORINA software (Molecular Networks GmbH: Erlangen, Germany) [15], a process that displayed 13 molecules as structural outliers. For building the model, the dataset was divided into training and test sets, comprised of 1017 and 339 compounds, respectively. Four compounds were excluded as they showed atypical behavior in the clusters division. From the amalgamation schedule of the hierarchical cluster analysis, nine clusters were determined and posteriorly used in the K-means clustering experiment where 25% of the compounds from each cluster were used as the test set.

**Molecular descriptors selection:** The x/x correlation filter (based on a correlation coefficient threshold of 0.90) and the standardized entropy measure (less than 0.30) available in the MobyDigs software (TALETE srl., Milano, Italy) [16] allowed for the reduction of the original data matrix to 478 MD. Posteriorly, supervised feature selection based on multiple linear regression coupled with the genetic algorithm was employed to select 14 MD, considered as the most significant for the modeling of DPPH• free radical scavenging capacity, i.e., MATS2e, BELe6, HATS3u, H2v, R7v, nN-N, nImidazoles, C-005, C-020, O-057, O-060, GVWAI-50, B02 (O-S) and B07 (O-S).

**QSAR method:** The multilayer perceptron (MLP) neural model was constructed using the DPPH• scavenging capacity of 1356 molecules. Outlier diagnosis was performed using the William's plot where data points with residuals greater than two standard deviation units were deemed to be outliers. Additionally, an analysis of the leverage values was performed to avoid the model's dependence on particular data points. Therefore, 31 compounds were separated from the model and the network was finally constructed with 1325 molecules. Several options were explored, giving, as the best final configuration of the network, the one based on the Quasi-Newton method as the training algorithm: Broyden-Fletcher-Goldfarb-Shanno [17], with 90 interaction cycles (BFGS 90); the error function was the sum-of-squares (SOS), and Tanh and Logistic were employed as activation functions used for the hidden and output layers, respectively. The 14 MD selected as the most significant variables were used as input neurons in the neural network model and the remaining architecture was comprised of nine neurons in the hidden layer and one output (MLP 14-9-1).
## 2.2. Performance and Predictive Capacity of the Model

The neural MLP model was trained and posteriorly validated over the test set. The correlation coefficient for the training set was R² = 0.713, and for test set it was Q²ₐxt = 0.654. The relation between the targets, output and standard residuals pIC₅₀ values on the MLP model for the studied compounds is shown in Figure 1 below; the correlation between the experimental and predicted values for the training and test samples is shown in Figure 2. The whole pIC₅₀ target and predicted values are available in Supporting Information (Tables S1 and S2).

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Regression plane of relation between targets, output and standard residuals values of the analyzed variable (pIC₅₀).</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Correlation between experimental and predicted pIC₅₀.</figcaption></figure>
## 2.3. Relative Importance of the Variables in the Model

Sensitivity analysis (SA): In order to determine the relative importance of the MD used as variables in the ANN model, SA was performed. Sensitivity is a statistical parameter measured as the difference between standard deviation (SDE) values when each MD is considered as an input (SDE(n)), and when the same MD is excluded (SDE(n-1)), with both values computed over the same dataset. Greater differences are associated with higher relevance for the excluded MD [17]. Figure 3 shows the SA for the 14 MD employed in the construction of the MLP 14-9-1 network model.

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Sensitivity analysis of the MD for the MLP 14-9-1 model.</figcaption></figure>

Principal component analysis (PCA): The variable importance was also explored with a PCA using STATISTICA 8.0 software (StatSoft Inc., Tulsa, OK, USA) [17]. This method offers information about the significance of the variables and their correlation which can be used as a diagnostic tool for model interpretation. The results are shown in Figure 4 below.

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Diagram of variable importance according to the PCA.</figcaption></figure>

## 2.4. Prediction of Coumarin Derivatives Scavenging

The MLP model was used to predict the DPPH• scavenging capacity of coumarin derivatives. Several studies have suggested that these chemical compounds possess the appropriate characteristics for potential antioxidant activity [8–11]. However, bearing in mind that any model may only appropriately predict the behavior of data points in its applicability domain (AD), the inclusion of these compounds in the model's chemical space was assessed using the Ambit Discovery software (Nina
Jeliazkova, Sofia, Bulgaria) [18]. With this analysis, the coumarin-type compounds were confirmed to lie in the network's AD. The predicted results (expressed as pIC₅₀ values) for each molecule are depicted in Table 1.

<table><caption>**Table 1.** Predictions of the pIC₅₀ values.</caption><thead><tr><th>N°</th><th>Structures</th><th>IUPAC Name</th><th>Predicted pIC<sub>50</sub></th></tr></thead><tbody><tr><td>1</td><td><img src="image_1.png" alt="Chemical structure of 4-(4-(trifluoromethyl)phenyl)-3,4-dihydro-2-methoxy-2-methylpyrano[3,2-c]chromen-5(2H)-one"></td><td>4-(4-(trifluoromethyl)phenyl)-3,4-dihydro-2-methoxy-2-methylpyrano[3,2-c]chromen-5(2H)-one</td><td>3.881</td></tr><tr><td>2</td><td><img src="image_2.png" alt="Chemical structure of 3,4-dihydro-2-methoxy-2-methyl-4-(4-nitrophenyl)pyrano[3,2-c]chromen-5(2H)-one"></td><td>3,4-dihydro-2-methoxy-2-methyl-4-(4-nitrophenyl)pyrano[3,2-c]chromen-5(2H)-one</td><td>3.951</td></tr><tr><td>3</td><td><img src="image_3.png" alt="Chemical structure of 3,4-dihydro-2-methoxy-4-(4-methoxyphenyl)-2-methylpyrano[3,2-c]chromen-5(2H)-one"></td><td>3,4-dihydro-2-methoxy-4-(4-methoxyphenyl)-2-methylpyrano[3,2-c]chromen-5(2H)-one</td><td>3.829</td></tr><tr><td>4</td><td><img src="image_4.png" alt="Chemical structure of 3,4-dihydro-2-methoxy-2-methyl-4-phenylpyrano[3,2-c]chromen-5(2H)-one"></td><td>3,4-dihydro-2-methoxy-2-methyl-4-phenylpyrano[3,2-c]chromen-5(2H)-one</td><td>3.944</td></tr><tr><td>5</td><td><img src="image_5.png" alt="Chemical structure of 3,4-dihydro-2-methoxy-2-methyl-4-p-tolylpyrano[3,2-c]chromen-5(2H)-one"></td><td>3,4-dihydro-2-methoxy-2-methyl-4-p-tolylpyrano[3,2-c]chromen-5(2H)-one</td><td>3.946</td></tr><tr><td>6</td><td><img src="image_6.png" alt="Chemical structure of 4-(4-fluorophenyl)-3,4-dihydro-2-methoxy-2-methylpyrano[3,2-c]chromen-5(2H)-one"></td><td>4-(4-fluorophenyl)-3,4-dihydro-2-methoxy-2-methylpyrano[3,2-c]chromen-5(2H)-one</td><td>3.959</td></tr><tr><td>7</td><td><img src="image_7.png" alt="Chemical structure of 4-(4-tert-butylphenyl)-3,4-dihydro-2-methoxy-2-methylpyrano[3,2-c]chromen-5(2H)-one"></td><td>4-(4-tert-butylphenyl)-3,4-dihydro-2-methoxy-2-methylpyrano[3,2-c]chromen-5(2H)-one</td><td>3.903</td></tr><tr><td>8</td><td><img src="image_8.png" alt="Chemical structure of 4-hydroxy-3-(3-oxo-1-phenylbutyl)-2H-chromen-2-one"></td><td>4-hydroxy-3-(3-oxo-1-phenylbutyl)-2H-chromen-2-one</td><td>2.326</td></tr><tr><td>9</td><td><img src="image_9.png" alt="Chemical structure of 4-hydroxy-3-(1-(4-methoxyphenyl)-3-oxobutyl)-2H-chromen-2-one"></td><td>4-hydroxy-3-(1-(4-methoxyphenyl)-3-oxobutyl)-2H-chromen-2-one</td><td>2.291</td></tr><tr><td>10</td><td><img src="image_10.png" alt="Chemical structure of 3-(1-(4-(trifluoromethyl)phenyl)-3-oxobutyl)-4-hydroxy-2H-chromen-2-one"></td><td>3-(1-(4-(trifluoromethyl)phenyl)-3-oxobutyl)-4-hydroxy-2H-chromen-2-one</td><td>2.217</td></tr><tr><td>11</td><td><img src="image_11.png" alt="Chemical structure of 4-hydroxy-3-(1-(4-nitrophenyl)-3-oxobutyl)-2H-chromen-2-one"></td><td>4-hydroxy-3-(1-(4-nitrophenyl)-3-oxobutyl)-2H-chromen-2-one</td><td>2.288</td></tr><tr><td>12</td><td><img src="image_12.png" alt="Chemical structure of 4-hydroxy-3-(3-oxo-1-p-tolylbutyl)-2H-chromen-2-one"></td><td>4-hydroxy-3-(3-oxo-1-p-tolylbutyl)-2H-chromen-2-one</td><td>2.318</td></tr><tr><td>13</td><td><img src="image_13.png" alt="Chemical structure of 3-(1-(4-fluorophenyl)-3-oxobutyl)-4-hydroxy-2H-chromen-2-one"></td><td>3-(1-(4-fluorophenyl)-3-oxobutyl)-4-hydroxy-2H-chromen-2-one</td><td>2.341</td></tr><tr><td>14</td><td><img src="image_14.png" alt="Chemical structure of 3-(1-(4-tert-butylphenyl)-3-oxobutyl)-4-hydroxy-2H-chromen-2-one"></td><td>3-(1-(4-tert-butylphenyl)-3-oxobutyl)-4-hydroxy-2H-chromen-2-one</td><td>2.179</td></tr><tr><td>15</td><td><img src="image_15.png" alt="Chemical structure of 4-hydroxy-2H-chromen-2-one"></td><td>4-hydroxy-2H-chromen-2-one</td><td>3.421</td></tr></tbody></table>
## 2.5. *In Vitro* Assay

The *in vitro* study of 4-hydroxycoumarin (compound 15) was developed as an experimental corroboration of the predicted value. The pIC₅₀ experimental result attained, according to the method described in Section 4.3, was 3.443.

# 3. Discussion

## 3.1. *Database and Neural Network*

The database of 1373 compounds with their corresponding DPPH• free radical activity values is, to the best of the authors' knowledge, the most diverse and largest that has been reported until this moment, and it will allow for deeper study of the structure-antiradical activity relationships of chemical compounds.

The calibration (R²) and external validation (Q²ext) values (i.e., 0.713 and 0.654, respectively) of the MLP model are above the limits established for model acceptance [19], and thus indicate the satisfactory fitness and predictive capacity of the obtained model.

## 3.2. *Analysis of the Molecular Descriptors*

The analysis of the relative importance and a mechanistic interpretation of the 14 variables included in the MLP model provide a deeper understanding of the chemical information codified and its relationship with the modeled property. As may be observed in Figures 3 and 4, the results obtained in the analysis of the relative importance of each MD based on the SA and PCA methods, respectively, exhibit great similarity. A more detailed analysis of these variables revealed the following:

- All the one-dimensional (1D) molecular descriptors included in the 14 variables, *O*-060, *O*-057, *C*-005, *C*-020, *nN-N*, *n*-imidazoles, are "Atom-centered fragments" or "Functional group counts" descriptors based in functional groups with the presence of some electronegative atom (O, N, S, Se, halogens) [20].

- On the other hand, 2D molecular descriptors are related with the presence or absence of the O-S bond (in the case of B02 (O-S) and B07 (O-S)) or, such as MATS2e y BELe6, described properties correlated with Sanderson atomic electronegativity [20].

- The 3D selected MD (HATS3u, H2v y R7v) are GATEWAY descriptors weighted by atomic van der Waals volumes [20].

- The GVWAI-50 is a drug-like molecular properties descriptor, and its values are provided only for compounds having C, H, O, N, S, Se, P, B, Si, and halogens [20].

The chemical interpretation on the importance of these MD may be retrieved from the analysis of the two currently accepted mechanisms for antioxidant activity [21], specifically to deactivate a free radical. One of them is the H-atom transfer, in which the radical reaction chain is interrupted. The intermediate molecule in this sequence of reactions has to be stable, a condition that may be obtained in molecules with electron-donating heteroatoms, or high electronic delocalization. The other mechanism is the electron transfer, in which the radical cation is first formed, followed by rapid and reversible deprotonation in solution; however, if the radical cation formed and has a sufficient lifetime, it can attack biologic molecules, even causing mutagenic effects. In both mechanisms, the presence of high electronegative atoms contributes positively to the neutralization of intermediate species, thus preventing cellular damage [22–29]. Additionally, weighting the atoms according to their chemical environment allows us to evaluate the capacity of the molecules to interact with DPPH•, and thus their scavenging capacity.

## 3.3. *Prediction of Antioxidant Activities of a Group of Coumarins*

Recent advances in drug discovery have resulted in an increase in the number of synthetic and naturally occurring molecules available for testing using *in vitro* assays for the scavenging ability of the
DPPH• radical. Virtual screening allows for prior assessment of the potential bioactivity of chemical compounds, and thus providing key guidelines in posterior experimental work [30,31].

Coumarins form a large class of phenolic compounds occurring in plants [32]. There are numerous research initiatives aimed at studying the effects of coumarins with several positions of the hydroxyl groups and other substitutions on the scavenging activity of different radicals, including DPPH• [22–26]. The series of coumarin-type compounds used in this study may be divided for analysis into two groups, according to the structural analogy: Cy-analog(Compounds 1–7): *Cyclocoumarol analogous* and Wf-analog(Compounds 8–14): *Warfarine analogous*.

The DPPH• scavenging capacity predictions for the group of coumarin derivatives using the MLP model constitute the first approximation on the degree of possible antioxidant activity for this group of compounds.

The two groups of compounds have significantly different values of pIC₅₀, as can be observed in Table 1. Wf-analogs clearly seem to be less effective in DPPH• radical capturing because their values of pIC₅₀ were much lower (below 2.4). The non-substituted compound 4-hydroxycoumarin (compound 15) showed an intermediate value of pIC₅₀ (3.421). On the other hand, the pIC₅₀ Cy-analog values were much higher (over 3.8). These results indicated a superior ability of cyclocoumarol derivatives for scavenging the DPPH• radical. Nevertheless, a more detailed analysis of the structural features is needed.

### 3.4. *In Vitro Assay*

The *in vitro* corroboration of the MLP model prediction showed satisfactory proximity between the experimental and predicted pIC₅₀ values (i.e., 3.443 and 3.421, respectively). Therefore, the built MLP model may reliably be applied in the search for new antioxidant compounds.

## 4. Materials and Methods

### 4.1. *Preparation of Cases and Variables*

The results of the experimentally determined scavenging ability of the DPPH• radical, (expressed in terms of the IC₅₀ values) for 1373 molecules extracted from 181 scientific reports in the literature (See Supplementary Information, Table S3) produced a comprehensive and diverse dataset of compounds for the posterior analysis. All the structures were optimized using CORINA software (Molecular Networks GmbH: Erlangen, Germany) [15], and STANDARDIZER software available in the ChemAxon package (ChemAxon Ltd., Budapest, Hungary) [33]. Outlier diagnosis of the modeled compounds was performed in different moments of the study.

Furthermore, the dataset was rationally divided into training and test sets using hierarchical and K-Means clustering methods available in the STATISTICA 8.0 software (StatSoft Inc., Tulsa, OK, USA) [17]. Complete linkage was employed as the linkage rule and squared Euclidean distance as distance measure in the case of the former, while the optimal number of clusters for the latter were determined from the amalgamation schedule of obtained joining tree.

The parameterization of the structures was performed using 3224 molecular descriptors implemented in the DRAGON 5.5 software (TALETE srl, Milano, Italy) [34]. The correlation filter of Dragon software was applied to reduce the number of variables. Additionally, multiple linear regression coupled with the genetic algorithm in MobyDigs software (TALETE srl, Milano, Italy) [34], was employed to select the final subset of variables used in the ANN building.

The coumarin derivatives were analyzed with the Ambit Discovery software (Nina Jeliazkova, Sofia, Bulgaria) [18] to assess if they were included in the applicability domain of the MLP model, and were later optimized, and parameterized with the molecular descriptors contained in the built model.
## 4.2. *Development of ANN Model*

The QSAR model was developed as a Multilayer Perceptron Neural Network using STATISTICA 8.0 software (StatSoft Inc., Tulsa, OK, USA) [18]. The response variable values from the scavenging ability of the molecules (IC₅₀) were transformed to their corresponding pIC₅₀ values (–log IC₅₀).

## 4.3. *In Vitro DPPH• Assay*

The free radical scavenging activity of the 4-hydroxycoumarin was measured using the stable DPPH• radical, according to Blois's method [35]. Briefly, 3 mL of each sample solution was prepared in methanol at different concentrations (150–750 µg/mL) and was added, individually, 1 mL of DPPH• solution (0.1 mM). The mixture was shaken vigorously and left in the dark for 30 min. Then, the absorbance was measured in a Spectrophotometer (Thermo Scientific™ GENESYS 10S UV-Vis, Waltham, MA, USA) at 517 nm. BHT was used as reference in the experimental assay. This procedure was repeated three times for reproducibility. The capability to scavenge the DPPH• radical was expressed as IC₅₀ (concentration of antioxidant that produces 50% of absorbance inhibition).

# 5. Conclusions

The scavenging capacity of the DPPH• radical is one of the most common methods for evaluating *in vitro* antiradical activity. An MLP neural network model was constructed to relate the structure of 1373 molecules with their scavenging activity. This model was validated using both internal and external validation techniques, showing a good predictive ability. The constructed network was used to predict the antioxidant activity of a set of coumarin-type compounds. An *in vitro* assay to further validate the predictive capacity of the built model demonstrated satisfactory proximity between experimental and predicted values, thus corroborating the performance of the model.

**Supplementary Materials:** Supplementary materials can be found at http://www.mdpi.com/1422-0067/17/6/881/s1.

**Acknowledgments:** Ministry of National Education, Research and Technology is gratefully acknowledged for the graduate fellowship awarded to Anita Maria Rayar. The authors gratefully acknowledge support for part of this work through the program PYTHAGORAS II of EPEAEK II (MIS: 97436/073). We would like to thank Biobyte Corp. 201 West 4th Street, Suite 204, Claremont, CA 91711, USA, for free access to the C-QSAR program. The authors also thank the "Service de Coopération et d'Action Culturelle" of Embassy of France in Cuba for the financial support.

**Author Contributions:** Elizabeth Goya Jorge and Stephen J. Barigye conceived and designed the chemometric model and are responsible for its evaluation. María Elisa Jorge Rodríguez performed the *in vitro* antioxidant assay. The French team, Anita Maria Rayar and Maité Sylla-Iyarreta Veitía, provided the coumarin derivative structures to develop the chemical library and contributed to the analysis of results. All authors contributed to the drafting and revision of the article and approved the final version.

**Conflicts of Interest:** The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the article, and in the decision to publish the results.

## References

1. Valko, M.; Leibfritz, D.; Moncol, J.; Cronin, M.; Mazur, M.; Telser, J. Free radicals and antioxidants in normal physiological functions and human disease. *Int. J. Biochem. Cell Biol.* **2007**, *39*, 44–84. [CrossRef] [PubMed]

2. Magalhaes, L.M.; Segundo, M.A.; Reis, S.; Lima, J.L.F.C. Methodological aspects about *in vitro* evaluation of antioxidant properties. *Anal. Chim. Acta* **2008**, *613*, 1–19. [CrossRef] [PubMed]

3. Huang, D.; Ou, B.; Prior, R.L. The chemistry behind antioxidant capacity assays. *J. Agric. Food Chem.* **2005**, *53*, 1841–1856. [CrossRef] [PubMed]

4. Scalzo, J.; Mezzetti, B.; Battino, M. Total antioxidant capacity evaluation: Critical steps for assaying berry antioxidant features. *BioFactors* **2005**, *23*, 221–227. [CrossRef] [PubMed]

5. Niki, E. Assessment of Antioxidant Capacity *in vitro* and *in vivo*. *Free Rad. Biol. Med.* **2010**, *49*, 503–515. [CrossRef] [PubMed]
6. Gunars, T.; Grzegorz, B. Determination of antiradical and antioxidant activity: Basic principles and new insights. *Acta Biochim. Pol.* **2010**, *57*, 139–142.

7. Kedare, S.; Singh, R. Genesis and development of DPPH method of antioxidant assay. *J. Food Sci. Technol.* **2011**, *48*, 412–422. [CrossRef] [PubMed]

8. Naceur, H.; Fischmeister, C.; Puerta, M.; Valerga, P. A rapid access to new coumarinyl chalcone and substituted chromeno[4,3-c]pyrazol-4(1H)-ones and their antibacterial and DPPH radical scavenging activities. *Med. Chem. Res.* **2011**, *20*, 522–530.

9. Puerta, M.; Naceur, H.; Valerga, P. Synthesis, structure, antimicrobial and antioxidant investigations of dicoumarol and related compounds. *Eur. J. Med. Chem.* **2008**, *43*, 2541–2548.

10. Fylaktakidou, K.; Hadjipavlou-Litina, D.; Litinas, K.; Nicolaides, D. Natural and synthetic coumarin derivatives with antiinflammatory/antioxidant activities. *Curr. Pharm.* **2004**, *10*, 3813–3833. [CrossRef]

11. Razo-Hernández, R.; Pineda-Urbina, K.; Velazco-Medel, M.; Villanueva-García, M.; Sumaya-Martínez, T.; Martínez-Martínez, F.; Gómez-Sandoval, Z. QSAR study of the DPPH• radical scavenging activity of coumarin derivatives and xanthine oxidase inhibition by molecular docking. *Cent. Eur. J. Chem.* **2014**, *12*, 1067–1080. [CrossRef]

12. Todeschini, R.; Consonni, V.; Gramatica, P. Chemometrics in QSAR. In *Comprehensive Chemometrics Chemical and Biochemical Data Analysis*; Brown, S., Tauler, R., Walczak, B., Eds.; Elsevier: Ámsterdam, The Netherlands, 2009; Volume 4, pp. 129–172.

13. Eriksson, L.; Jaworska, J.; Worth, A.; Cronin, M.; McDowell, R.; Gramatica, P. Methods for reliability and uncertainty assessment and for applicability evaluations of classification- and regression-based QSARs. *Environ. Health Perspect.* **2003**, *111*, 1361–1375. [CrossRef] [PubMed]

14. Aoyama, T.; Suzuki, Y.; Ichikawa, H. Neural networks applied to structure-activity relationships. *J. Med. Chem.* **1990**, *33*, 905–908. [CrossRef] [PubMed]

15. CORINA Classic (software); Molecular Networks GmbH: Erlangen, Germany, 2015. Available online: http://www.molecular-networks.com (accessed on 9 March 2015).

16. TALETE srl. *MobyDigs* (software), v1.0; Milano, Italy, 2004. Available online: http://www.talete.mi.it/mobydigs.htm (accessed on 2 June 2015).

17. StatSoft, Inc. *STATISTICA* (software), v8.0; Tulsa, OK, USA, 2007. Available online: http://www.statsoft.com (accessed on 7 June 2015).

18. Jeliazkova, N. *Ambit Discovery* (software), v0.04; Sofia, Bulgaria, 2005. Available online: http://www.ambit.acad.bg (accessed on 5 June 2015).

19. Tropsha, A. Best Practices for QSAR model development, validation, and exploitation. *Mol. Inf.* **2010**, *29*, 476–488. [CrossRef]

20. Todeschini, R.; Consonni, V. *Handbook of Molecular Descriptors*, 1st ed.; Wiley-VCH: Weinheim, Germany, 2000.

21. Wright, J.; Johnson, E.; DiLabio, G. Predicting the activity of phenolic antioxidants: Theoretical method, analysis of substituent effects, and application to major families of antioxidants. *J. Am. Chem. Soc.* **2001**, *123*, 1173–1183. [CrossRef] [PubMed]

22. Chen, Y.; Xiao, H.; Zheng, J.; Liang, G. Structure-thermodynamics-antioxidant activity relationships of selected natural phenolic acids and derivatives: An experimental and theoretical evaluation. *PLoS ONE* **2015**, *10*. [CrossRef] [PubMed]

23. Gupta, S.; Matthew, S.; Abreu, P.M.; Aires-de-Sousa, J. QSAR analysis of phenolic antioxidants using MOLMAP descriptors of local properties. *Bioorg. Med. Chem.* **2006**, *14*, 1199–1206. [CrossRef] [PubMed]

24. Jing, P.; Zhao, S.-J.; Jian, W.-J.; Qian, B.-J.; Dong, Y.; Pang, J. Quantitative studies on structure-DPPH scavenging activity relationships of food phenolic acids. *Molecules* **2012**, *17*, 12910–12924. [CrossRef] [PubMed]

25. Khlebnikov, A.I.; Schepetkin, I.A.; Domina, N.G.; Kirpotina, L.N.; Quinn, M.T. Improved quantitative structure-activity relationship models to predict antioxidant activity of flavonoids in chemical, enzymatic, and cellular systems. *Bioorg. Med. Chem.* **2007**, *15*, 1749–1770. [CrossRef] [PubMed]

26. Mitra, I.; Saha, A.; Roy, K. Quantitative structure-activity relationship modeling of antioxidant activities of hidroxybenzalacetones using quantum chemical, physicochemical and spatial descriptors. *Chem. Biol. Drug. Des.* **2009**, *73*, 526–536. [CrossRef] [PubMed]

27. Mitra, I.; Saha, A.; Roy, K. Chemometric modeling of free radical scavenging activities of flavone derivatives. *Eur. J. Med. Chem.* **2010**, *45*, 5071–5079. [CrossRef] [PubMed]
28. Worachartcheewan, A.; Nantasenamat, C.; Isarankura-Na-Ayudhya, C.; Prachaiasittikul, S.; Prachaiasittikul, V. Predicting the free radical scavenging activity of curcumin derivatives. *Chemometr. Intell. Lab. Syst.* **2011**, *109*, 207–216. [CrossRef]

29. Yamagami, C.; Motohashi, N.; Emoto, T.; Hamasaki, A.; Tanahashi, T.; Nagakura, N.; Takeuchi, Y. Quantitative structure-activity relationship analyses of antioxidant and free radical scavenging activities for hydroxybenzalacetones. *Bioorg. Med. Chem. Lett.* **2004**, *14*, 5629–5633. [CrossRef] [PubMed]

30. Kadhum, A.; Al-Amiery, A.; Musa, A.; Mohamad, A. The antioxidant activity of new coumarin derivatives. *Int. J. Mol. Sci.* **2011**, *12*, 5747–5761. [CrossRef] [PubMed]

31. McInnes, C. Virtual screening strategies in drug discovery. *Curr. Opin. Chem. Biol.* **2007**, *11*, 494–502. [CrossRef] [PubMed]

32. Cseke, L.J.; Kirakosyan, A.; Kaufman, P.B.; Warber, S.; Duke, J.A.; Brielmann, H.L. *Natural Products from Plants*, 2nd ed.; CRC Press: New York, NY, USA, 2006.

33. ChemAxon Ltd. *ChemAxon, 6.1.0 (software)*; Budapest, Hungary, 2013. Available online: http://www.chemaxon.com (accessed on 10 March 2015).

34. TALETE srl. *DRAGON for Windows (Software for Molecular Descriptor Calculations)*, V 5.5; Milano, Italy, 2007. Available online: http://www.talete.mi.it/ (accessed on 15 March 2015).

35. Blois, M. Antioxidant determinations by the use of a stable free radical. *Nature* **1958**, *181*, 1199–1200. [CrossRef]

<figure><img src="image_7.png" /></figure>

© 2016 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC-BY) license (http://creativecommons.org/licenses/by/4.0/).