Computational and Theoretical Chemistry 1210 (2022) 113658

Contents lists available at ScienceDirect

Computational and Theoretical Chemistry

journal homepage: www.elsevier.com/locate/comptc

A DFT study of the chemical reactivity properties, spectroscopy and

bioactivity scores of bioactive flavonols
Eric O. Akintemi a, Krishna K. Govender b, c, Thishana Singh a, *
a
School of Chemistry and Physics, University of KwaZulu-Natal, P.M.B. X54001, Durban 4000, South Africa
b
Department of Chemical Sciences, University of Johannesburg, P.O. Box. 17011, Doornfontein Campus, 2028 Johannesburg, South Africa
c
Council for Scientific and Industrial Research, National Integrated Cyber Infrastructure, Centre for High Performance Computing, 15 Lower Hope Road, Rosebank, Cape
Town 7700, South Africa

A R T I C L E I N F O A B S T R A C T

Keywords: Density function theory calculations was used to determine the molecular parameters, electronic and chemical
DFT reactivity descriptors, spectroscopy, and non-linear optical properties, electronic dipole moment, polarizability
Flavonols and hyperpolarizability of fifteen (15) flavonol aglycones (no sugar moiety) from plant sources to investigate
Chemical reactivity properties
their possible application as drug candidates. Geometry optimisations was carried out using the hybrid functional
Spectroscopy
Bioactivity scores
and basis set: M06-2X/6-31+G(d,p). Our calculations show that all the flavonols investigated are chemically
reactive. Their reactivity is greatest in water hence making them suitable drug candidates since this is the ideal
medium for drug delivery. The highest negative charge on the oxygen atom of the hydroxyl and high positive
charge on the H atom of C5-OH are vital for antioxidant activity. The most reactive species, from reactivity
descriptor calculations, is Gossypetin. All the flavonols are active as enzyme inhibitors and moderately active as
G-protein-coupled receptors, ion channel modulators and protease inhibitors. The physicochemical properties
show the flavonols have good bioavailability. All the compounds agree with Lipinski’s rule of 5, signifying
potential use as oral active drugs.

1. Introduction leaves, flowers and outer part of the plant namely the skin and peel, with
reduction in concentration towards the core. Only minute quantities of
Flavonoids are a class of multi-phenolic secondary metabolites found flavonols are present in parts of the plant below the soil surface, with the
in fruits and edible vegetables [1]. Over 9000 flavonoids have been significant exception of onions [11–14].
isolated from different plant sources [2]. The health benefits and dietary Theoretically calculated parameters to depict and explain the radical
importance of polyphenols has drawn the attention of researchers [3,4] scavenging activity of antioxidant phenols were first proposed, dis­
due to their high anti-oxidant, anti-inflammatory, anti-bacterial, anti- cussed and tested in computational studies performed in the 1990s
viral and anti-tumor properties [1,5,6]. Research reports confirm that [15–17]. The correlation of experimentally determined phenolic O–H
these properties arise from their ability to quench free radicals [7,8]. As bond dissociation enthalpies (BDE) and half-peak oxidation potentials
a result, there is an emerging study in the area of food, chemical and (Ep/2) [18,19] with the rate constants for reactions of phenols with
pharmaceuticals towards unraveling sources of anti-oxidants to inhibit different free radicals prompted extensive studies of theoretical pa­
radical induced ailments. Furthermore, the clarification of polyphenol rameters associated with physicochemical properties, which are actually
anti-oxidant activities is of vital importance in gaining knowledge about the computed BDE and ionization potentials, IP.
other biological features [9]. Mendes et al., [20] carried out a computational investigation on the
The flavonoid family comprises of organic compounds characterised antioxidant potential of myrcetin 3,4-di-O-α-L-rhamnopyranoside, a
by phenol structural units. Based on structural differences, the poly­ double glycosylated derivative of Myricetin. The bond dissociation en­
phenols are classified as flavones, flavonols, chalcones, lignans and ergy obtained suggests that the derivative has an antioxidant potential
hydroxycoumarins. Flavonols have diverse chemical structure and that is as good as the parent molecule. The bond dissociation enthalpy of
characteristics [10]. Flavonol glycosides are found primarily in the the O–H group was computationally determined by Wright and co-

* Corresponding author.
E-mail address: [email protected] (T. Singh).

https://doi.org/10.1016/j.comptc.2022.113658
Received 8 January 2022; Received in revised form 21 February 2022; Accepted 26 February 2022
Available online 1 March 2022
2210-271X/© 2022 Elsevier B.V. All rights reserved.
E.O. Akintemi et al. Computational and Theoretical Chemistry 1210 (2022) 113658

workers [7] and was the theoretical parameter successfully used to Geometry optimisations was carried out at 298.15 K in vacuum using the
measure the H-atom donor ability of flavonoids. An attempt was made hybrid functional and basis set: M06-2X/6-31+G(d,p) [30]. This was
by Woldu and Mai [21] to investigate the thermodynamics and kinetics chosen since it showed good results in published reports [20,21]. Using
of hydrogen atom withdrawal from flavonoids and related phenolic the gas-phase optimized structures, a single-point frequency calculation
compounds using theoretical computations. This study shed light on the in water (ε = 78.3553) and n-octanol (ε = 9.8629) was done to obtain
fact that some common elements such as hydrogen, influence the ther­ the free energy for evaluating the octanol–water partition co-efficient.
modynamics and kinetics of the chain breaking ability of flavonoids. All frequency calculations in solvent was performed using the SCRF
From their theoretical findings they further predicted that for flavonoids keyword together with the SMD [31] solvent model. A single-point gas-
and related polyphenols to serve as radical scavengers, they must donate phase optimisation was carried out using the time-dependent density
a pair of hydrogen atoms where the total enthalpy of reaction does not functional theory (TD-DFT) [32] keyword and solving for 20 excited
exceed the enthalpy of the homolytic dissociation of the hydrogen states. The resulting molecules were submitted for optimisation in water
molecule (H2), 104.206 kcal mol− 1. and ethanol (ε = 24.852) which were the experimental solvents for UV.
Weber and co-workers [22] assembled a dataset of flavonoid com­ Drug-likeness scores and bioactivity prediction was carried out using the
pounds with pre-determined Trolox Equivalent Antioxidant Capacity online program Molinspiration cheminformatics [33]. Also, provided in
(TEAC) in 2,2′ -azino-bis(3-ethylbenzothiazoline-6-sulfonic acid) radical the Supporting Information is a list of keywords used for the different
scavenging assay and developed different classifications and analytical calculations described above as well as the XYZ coordinates of all the gas
models using calculated descriptors, including bond lengths, bond or­ phase optimized structures.
ders, dipole moment, enthalpy of formation, polarizability, atom and
substituent charges, electronic and total energies, and frontier molecular 3. Results and discussion
orbital energies of the flavonoid compounds. The extrapolated Partial
Least Square (PLS) model listed the molecular polarizability and the 3.1. Chemical reactivity descriptors
total charge of the substituent at position 5 in A-ring (Fig. 1) as the most
important descriptors for TEAC. Similarly, Principal Component Anal­ In this study, DFT was used to give a concise structure-chemical
ysis (PCA), hierarchical cluster analysis and k-nearest neighbours’ reactivity report of the flavonols to provide an understanding of the
techniques by Werner and his team38 listed molecular polarizability, chemical reactivity properties. All the chemical reactivity descriptors
charge on the carbon at the position 3, and total charge at positions 5 are presented in Table 1. Energies from the frontier molecular orbitals
and 3′ as the chief descriptors for classifying flavonoids as having low (FMOs) were used to calculate chemical reactivity descriptors arising
and high radical capacities. from conceptual DFT (the contour surfaces of the FMOs are in the
Over the past few years, there has been increased interest in the use Supporting Information). The energy gap, Egap (Eq. (1)) is defined as the
of bioactive organic molecules in the area of solid state research energy difference between EL (the lowest unoccupied molecular orbital,
[23–26]. The reason for this interest is their potential as nonlinear optics LUMO) and EH (energy of the highest occupied molecular orbital,
in electro-optics and conductors. Organic compounds have a high degree HOMO). This determines the chemical reactivity and kinetic stability of
of delocalization since they are formed by weak van der Waal’s and a molecule. The lower the energy gap, the more reactive a molecule is,
hydrogen bonds. The result is that they are optically more nonlinear which indicates that, it is more susceptible to react with other molecular
than inorganic materials. species, such as for example proteins or enzymes.
In this study, density function theory (DFT), will be used to present
Egap = EL − EH (1)
detailed chemical reactivities for fifteen (15) flavonol aglycones (no
sugar moiety) of plant sources to investigate their possible application as It can be seen (Fig. 3) that the energy gap is slightly lower in water
drug candidates (Fig. 2). Quantum mechanical calculations were used to than in ethanol. Since the preferred medium for drug delivery is water,
determine the molecular parameters, electronic and chemical reactivity this means that all the flavonols are reactive, with Gossypetin having the
descriptors, spectroscopy, and non-linear optical properties. Electronic lowest Egap. The lower energy gap corresponds to a higher concentration
dipole moment, polarizability and hyperpolarizability, which are non- of the molecules, and hence greater the activity of the polyphenol in that
linear optical properties (NLO), were also calculated to demonstrate medium. Although Gossypetin and Myricetin have the same molecular
the application of bioactive compounds in materials science. formula (Table S2), their energy gaps are not the same (Table 1). This
difference can be attributed to the position of the hydroxyl substituent.
2. Computational methodology In Myricetin, the sixth hydroxyl substituent is situated at C5′ on the
phenyl chain while in Gossypetin, it is attached at C8 in the main
All the structures were obtained from PubChem [27]. With the aid of chromone skeleton. The lower energy gap for Gossypetin may be
the GaussView 6.0.16 [28], we verified that all the functional groups attributed to the more nucleophilic benzyl group in the chromone
and substituents at the various positions were correct. All the calcula­ backbone, and allows electronic interactions which are pivotal to the
tions were done using the Gaussian 16 Rev. B.01 [29] program. flavonols’ reactivity. Thus, more hydroxyl groups on the chromone
structure is favourable for the radical scavenging activity of flavonols.
Although 3-hydroxyflavone and Natsudaidain have high inhibitory
activity, their energy gaps indicate no reactivity in the aqueous layer,
but imply probable reactivity in the lipid layers.
The electronic chemical potential, µ [34], which characterizes the
escaping tendency of electrons, is expressed by Eq. (2). There was no
significant difference in the chemical potential for the 15 flavonols.
Chemical potential values in water range from − 4.50 eV (3-hydroxy­
flavone) to − 4.15 eV (Gossypetin) and in ethanol from − 4.47 eV (3-
hydroxyflavone) to − 4.12 eV (Gossypetin). The low chemical potential
of Gossypetin in water further supports this flavonols’ good chemical
reactivity.
EL + EH
μ= (2)
2
Fig. 1. Flavonol substituents at positions 5 and 3 .
′

2
E.O. Akintemi et al. Computational and Theoretical Chemistry 1210 (2022) 113658

Fig. 2. 2D Structures of the 15 flavonols investigated.

Global hardness, η [35], is the resistance to charge transfer and has Electronegativity, χ (Eq. (5)), is a measure of electron distribution in
been successful in the validation of chemical processes [36], is calcu­ the molecule [37]. The high electronegativity obtained for these mole­
lated using Eq. (3). Global softness, s (Eq. (4)), is the inverse of global cules indicate they are good electron acceptors. In water, 3-hydroxyfla­
hardness, η and demonstrates the molecule’s susceptibility to charge vone had the highest χ value (4.50 eV) while Gossypetin had the lowest
transfer. Gossypetin has the smallest η value (2.862 eV) and invariably (4.15 eV). A possible explanation for high electronegativity for 3-
the highest s value of 0.349 eV in water. On the other hand, Natsudai­ hydroxyflavone is the lack of the –OH group on the phenyl moiety
dain has η value of 3.107 eV and s value of 0.321 eV. Hence, Gossypetin which leads to the molecule having a partially positive side. Gossypetin,
is the softest molecule, most susceptible to charge transfer and the most on the other hand is almost surrounded by –OH groups, hence an even
chemically reactive. Natsudaidain is the hardest molecule of all the electron distribution which accounts for the reduced electronegativity.
flavonols, most resistant to charge transfer and the least reactive.
χ= − μ (5)
E − EH
η= L (3) The electrophilicity index (ω), describes the stabilization energy of
2
the molecule when saturated by electrons approaching from the sur­
1 roundings; and is expressed by Eq. (6) [37]. The electron donating
s= (4)
η power and the electron accepting power are denoted as ω− and ω+, (Eqs.

3
E.O. Akintemi et al. Computational and Theoretical Chemistry 1210 (2022) 113658

Table 1
Chemical reactivity descriptors for the 15 flavonols.a
EL EH Egap η S μ χ ω ω− ω+ Δω± μD I

3Hf(g) − 1.36 − 7.56 6.200 3.100 0.322 − 4.46 4.46 3.222 10.443 2.749 13.193 5.236 7.962
(w) − 1.42 − 7.58 6.161 3.080 0.324 − 4.50 4.50 3.291 11.856 2.849 14.706 8.012 6.153
(eth) − 1.39 − 7.56 6.168 3.084 0.324 − 4.47 4.47 3.249 11.747 2.793 14.540 7.714 6.184

Aza(g) − 0.98 − 7.33 6.346 3.173 0.315 − 4.15 4.15 2.725 9.561 2.085 11.646 4.899 7.626
(w) − 1.08 − 7.28 6.196 3.098 0.322 − 4.18 4.18 2.827 10.615 2.244 12.859 7.907 5.916
(eth) − 1.04 − 7.26 6.218 3.109 0.321 − 4.15 4.15 2.778 10.490 2.177 12.667 7.583 5.944

Fis(g) − 1.14 − 7.30 6.155 3.077 0.324 − 4.22 4.22 2.899 9.782 2.343 12.125 6.663 7.649
(w) − 1.25 − 7.23 5.979 2.989 0.334 − 4.24 4.24 3.017 11.029 2.534 13.564 10.188 5.880
(eth) − 1.22 − 7.22 5.999 2.999 0.333 − 4.22 4.22 2.968 10.907 2.466 13.374 9.830 5.894

Gal(g) − 1.42 − 7.47 6.049 3.024 0.330 − 4.44 4.44 3.271 10.434 2.850 13.284 4.617 7.925
(w) − 1.36 − 7.39 6.035 3.017 0.331 − 4.37 4.37 3.177 11.488 2.729 14.218 6.566 6.102
(eth) − 1.34 − 7.38 6.041 3.020 0.331 − 4.36 4.36 3.151 11.421 2.694 14.116 6.387 6.130

Gos(g) − 1.28 − 7.08 5.796 2.898 0.345 − 4.18 4.18 3.018 9.775 2.578 12.354 5.654 7.481
(w) − 1.28 − 7.01 5.724 2.862 0.349 − 4.15 4.15 3.009 10.886 2.584 13.470 8.201 5.687
(eth) − 1.26 − 6.99 5.732 2.866 0.348 − 4.12 4.12 2.975 10.796 2.536 13.333 7.974 5.707

Kfd(g) − 1.07 − 7.34 6.262 3.131 0.319 − 4.20 4.20 2.829 9.715 2.232 11.947 5.236 7.706
(w) − 1.12 − 7.30 6.179 3.089 0.323 − 4.21 4.21 2.874 10.737 2.307 13.044 6.726 5.979
(eth) − 1.09 − 7.28 6.189 3.094 0.323 − 4.19 4.19 2.841 10.651 2.263 12.914 6.622 6.007

Kfr(g) − 1.27 − 7.33 6.063 3.031 0.329 − 4.30 4.30 3.058 10.033 2.568 12.602 5.675 7.731
(w) − 1.25 − 7.22 5.969 2.984 0.335 − 4.23 4.23 3.005 10.993 2.521 13.514 8.272 5.899
(eth) − 1.22 − 7.20 5.980 2.990 0.334 − 4.21 4.21 2.973 10.910 2.476 13.387 8.040 5.922

Iso(g) − 1.11 − 7.32 6.210 3.105 0.322 − 4.21 4.21 2.864 9.749 2.287 12.036 6.086 7.681
(w) − 1.13 − 7.23 6.091 3.045 0.328 − 4.18 4.18 2.875 10.697 2.326 13.024 8.368 5.901
(eth) − 1.11 − 7.21 6.104 3.052 0.327 − 4.16 4.16 2.845 10.620 2.286 12.906 8.184 5.950

Mor(g) − 1.35 − 7.34 5.991 2.995 0.333 − 4.35 4.35 3.162 10.185 2.721 12.907 3.921 7.688
(w) − 1.25 − 7.15 5.893 2.946 0.339 − 4.20 4.20 3.000 10.943 2.532 13.476 6.106 5.812
(eth) − 1.24 − 7.14 5.902 2.951 0.338 − 4.19 4.19 2.978 10.886 2.501 13.388 5.875 5.834

Myr(g) − 1.32 − 7.34 6.014 3.007 0.332 − 4.33 4.33 3.122 10.123 2.663 12.786 6.476 7.717
(w) − 1.29 − 7.18 5.888 2.944 0.339 − 4.23 4.23 3.047 11.067 2.595 13.662 9.548 5.880
(eth) − 1.27 − 7.17 5.902 2.951 0.338 − 4.22 4.22 3.019 10.997 2.554 13.552 9.265 7.151

Nat(g) − 0.99 − 7.30 6.305 3.152 0.317 − 4.14 4.14 2.729 9.542 2.097 11.640 4.852 7.420
(w) − 1.20 − 7.41 6.215 3.107 0.321 − 4.31 4.31 2.990 11.068 2.446 13.514 7.127 5.892
(eth) − 1.15 − 7.38 6.230 3.115 0.320 − 4.26 4.26 2.925 10.899 2.360 13.259 6.888 5.905

Pac(g) − 1.14 − 7.22 6.084 3.042 0.328 − 4.18 4.18 2.880 9.699 2.335 12.034 3.342 7.386
(w) − 1.24 − 7.26 6.016 3.008 0.332 − 4.25 4.25 3.005 11.015 2.510 13.526 5.179 5.932
(eth) − 1.21 − 7.23 6.019 3.009 0.332 − 4.22 4.22 2.967 10.913 2.460 13.373 4.984 5.931

Que(g) − 1.33 − 7.32 5.998 2.999 0.333 − 4.33 4.33 3.126 10.121 2.672 12.793 5.173 7.714
(w) − 1.27 − 7.16 5.888 2.944 0.339 − 4.22 4.22 3.024 11.005 2.564 13.570 7.643 5.895
(eth) − 1.39 − 7.56 5.899 3.084 0.324 − 4.47 4.47 3.249 11.747 2.793 14.540 7.415 5.911

Rhz(g) − 1.05 − 7.32 6.274 3.137 0.318 − 4.19 4.19 2.800 9.665 2.193 11.859 8.619 7.675
(w) − 1.11 − 7.25 6.141 3.070 0.325 − 4.18 4.18 2.850 10.653 2.284 12.938 11.992 5.911
(eth) − 1.08 − 7.24 6.156 3.078 0.324 − 4.16 4.16 2.817 10.568 2.239 12.807 11.721 5.946

Rht(g) − 1.34 − 7.17 5.829 2.914 0.343 − 4.26 4.26 3.118 9.991 2.702 12.693 5.895 7.525
(w) − 1.30 − 7.10 5.806 2.903 0.344 − 4.20 4.20 3.044 11.018 2.610 13.692 8.630 5.814
(eth) − 1.28 − 7.09 5.806 2.903 0.344 − 4.19 4.19 3.027 10.972 2.587 13.560 8.419 5.824
a
Units for the chemical reactivity descriptors are in electron volts (eV) EL: LUMO energy, EH: HOMO energy, Egap: HOMO-LUMO energy gap, η: Chemical hardness,
S: Chemical softness, μ: Chemical potential, χ: Electronegativity, ω: Electrophilicity index, ω-: Electron donating power, ω+: Electron accepting power, Δω±: Net
electrophilicity, μD: Dipole moment, I: Ionization potential, g: gas phase, w: water and eth: ethanol.

4
E.O. Akintemi et al. Computational and Theoretical Chemistry 1210 (2022) 113658

Fig. 3. Comparison of the energy gaps in gas-phase, water and ethanol.

(7) and (8)) respectively [38]. The multiphilic descriptor [39] or net hyperpolarizability (β) are indicators for the application of organic
electrophilicty [1] Δω±, (Eq. (9)) gives a better description of intermo­ compounds as non-linear optical (NLO) materials are reported [43–47].
lecular reactivity [40]. The results show that on average the flavonols The polarizablities and hyperpolarizabilities illustrate the reaction of a
have an electron donating power, ω− of approximately 10 eV which is system in an applied electric field. They are essential features of atomic
much higher than the electron accepting power, ω+ of about 3 eV. This and molecular systems. The donor and acceptor substituents offer the
means that the flavonols are strongly nucleophilic and weakly required ground-state charge asymmetry, while the π-conjugated system
electrophilic. offers a pathway for the reallocation of electric charges as a response to
the control of electric fields [46].
(EL + EH )2 Since there were no experimental NLO values for these flavonols,
ω= (6)
4(EL − EH ) urea and pure tetramethylurea were chosen as the reference compounds
[48,49]. Urea and pure tetramethylurea have first hyperpolarizabilities
(3EH − EL )2 β, 0.32 × 10− 30 esu and 0.41 × 10− 30 esu respectively [50]. Using these
ω− = (7)
16η values as the comparative threshold for the study, it can be inferred that
all these natural flavonols have the capacity for non-linear optical ap­
(EH + 3EL )2 plications, since their β values are greater than those for both urea and
ω+ = (8)
16η pure tetramethylurea. Interestingly, the comparison shows that

Δω± = ω+ + ω− (9)
The ionization potential, I (Eq. (10)) is the difference between the Table 2
positively charged aromatic hydroxyl species (ArOH+) and neutral ar­ Non-linear optical properties.
omatic hydroxyl species (ArOH) [41,42]. Ionization potential is also α (esu) Δα (esu) β (esu)
inversely proportional to reactivity so the lower the ionization potential, 3-Hydroxyflavone − 1.47 × 10− 23
4.05 × 10− 24
0.71 × 10− 30

the higher the chemical reactivity. Generally, the ionization potentials Azaleatin − 1.89 × 10− 23
2.75 × 10− 24
1.41 × 10− 30

in water were the lowest when compared to the other solvents. This Fisetin − 1.65 × 10− 23
5.58 × 10− 24
1.63 × 10− 30
23 24 30
reiterates the fact that the chemical reactivity of the flavonols is best in Galagin − 1.62 × 10− 5.85 × 10− 0.59 × 10−
23 24 30
Gossypetin 1.90 × 10− 3.54 × 10− 2.00 × 10−
water since this medium provides favourable interaction via hydrogen
−
23 24 30
Isorhamnetin − 1.93 × 10− 1.44 × 10− 1.61 × 10−
bonding. Gossypetin has the lowest I value (5.687 eV) in water relative Kaempferide − 1.78 × 10− 23
6.24 × 10− 24
1.13 × 10− 30

to other flavonols thus making it the most chemically reactive molecule. Kaempferol − 1.72 × 10− 23
3.14 × 10− 24
1.30 × 10− 30
23 24 30
Morin − 1.82 × 10− 2.38 × 10− 1.25 × 10−
I = EArOH + − EArOH (10) Myricetin − 1.89 × 10− 23
1.88 × 10− 24
1.60 × 10− 30
23 24 30
Natsudaidain − 2.48 × 10− 3.19 × 10− 0.56 × 10−
The relationship between energy gap and ionization potential is Pachypodol − 2.03 × 10− 23
2.00 × 10− 24
1.16 × 10− 30

presented in the Supporting Information (Table S3). Quercetin − 1.85 × 10− 23

2.17 × 10− 24
1.34 × 10− 30
23 24 30
Rhamnetin − 1.97 × 10− 4.62 × 10− 1.31 × 10−
23 24 30
3.2. Polarizabilties and hyperpolarizabilities Rhamnazin − 2.07 × 10− 2.12 × 10− 2.10 × 10−

α is linear or mean polarizability, Δα is anisotropy of polarizability, and β is first

Polarizabilities (α), anisotropy of polarizability Δα, and first order order hyperpolarizability. See Supporting Information for equations.

5
E.O. Akintemi et al. Computational and Theoretical Chemistry 1210 (2022) 113658

Gossypetin and Rhamnazin (β values of 2.00 × 10− 30 esu and 2.10 × the scavenging activity of flavonols reveal that water gives higher con­
10− 30 esu, respectively, Table 2) are each around 6 times more active, centrates than ethanol. Rhaman et al. [55] reported that a water extract
Isorhamnetin and Myricetin are 5 times as active, Azaleatin, Kaemp­ solution of C. asiatica had higher polyphenols and flavonol content than
feride, Kaempferol, Morin, Quercetin and Rhamnetin are 4 times as those of 100% ethanol solution extract. The work of Fang Tian and co-
active, Pachypodol is 3 times as active while 3-hydroxyflavone, Galagin workers [56] also showed that the total polyphenol (TPP) content ob­
and Natsudaidain are 2 times as active as urea for non-linear optical tained from a water extract of Galla chinensis is much more than that
applications. obtained for the ethanol extract. These properties also reflected in their
Sankar and Philip [51] described NLO as the study of the interaction half maximum effective concentration (EC50) measurements.
of intense light fields with matter. They state that molecules which are
anisotropic and display a linear anisotropic polarizability tend to show 3.4. Molecular electrostatic potential surface
an optically isotropic behavior in the bulk when their orientation is
randomly distributed. In our study we found that all compounds exhibit The 3-D molecular electrostatic potential (MEP) surfaces (Fig. 5) aid
some measure of linear polarizability, α and anisotropy of polarizability, in determining the reactivity of molecules [57–60]. The blue region
Δα. Hence, all should have optical isotropic behavior if their orienta­ indicates the positive potential of the molecule and hence is the area
tions are randomly distributed. where there is proton repulsion. The red region has negative potential
It should be noted that the moderate NLO properties displayed by the and represents proton attraction [61]. The MEP surfaces for all the fla­
15 flavonols can be improved by derivatisation if materials scientists vonols show an overall region of positive potential around protons of the
consider bioactive compounds for applications in optoelectronics and as C7, C4′ and C3′ . Conversely, areas of negative potential are around the
conductors [23,24]. oxygen atoms of the hydroxyl groups C3 and C5 as well as the C4
carbonyl group. This further verifies the necessity of the hydroxyl and
carbonyl groups for effective antioxidant properties of flavonols.
3.3. Atomic charge and electrostatic potential Kumar and Pandey [62] report that the B ring hydroxyl configuration
is important for the reactive oxygen species (ROS) and reactive
Electrical charge in the molecule are the driving force of electrostatic hydrogen species (RHS) scavenging ability of flavonols. The reason for
interactions (Fig. 4, see the Supporting Information for all flavonols). It this is that the ring offers a hydrogen and an electron to the hydroxyl as
is proven that electron density or charge are important in chemical re­ well as the peroxyl radical, and still maintains a stable flavonoid radical.
actions and physicochemical properties of compounds [52].Charge- Since C4′ and C3′ are on the B ring and have an area of positive potential
based descriptors are widely employed as chemical reactivity indices or it will bind to the negatively charged side groups of amino acids such as
as a measure of weak intermolecular interactions. Many quantum- aspartic acid and glutamic acid. Conversely, the areas of negative po­
chemical descriptors are derived from the partial charge distribution tential will have an affinity for the positively charged side groups of
in a molecule or from the electron density on specific atoms [52]. The amino acids (lysine, arginine and histidine) and sugar moieties to
calculated net atomic charges, Q [53] describe the molecular polarity of become β-glycosides.
the molecules.
A uniform distribution of charge was found for all molecules. For 3.5. UV–vis spectroscopy
Gossypetin, which had the lowest energy gap and was the most reactive,
it had the highest negatively charged atom (− 0.714 a.u.) which is the UV–visible spectroscopy has emerged as the most favoured tech­
oxygen atom of the hydroxyl group at C3′ position (Fig. 4). This charge is nique to evaluate the phenolic properties and content of bioactive
significant for the antioxidant property of the molecules due to the role compounds. The phenolic ring in bioactive molecules have the ability to
the –OH group plays in scavenging radicals. On the other hand, the absorb UV light. Some phenolic compounds also show absorption in the
highest positively charged atom (0.518 a.u.) is the hydrogen atom of the visible region since they are coloured. Hence, these features, make
hydroxyl (–OH) group at the C5 position of the ring A which is also UV–visible spectroscopy the appropriate technique to examine bioactive
important for the reactivity of the compound since the presence of –OH phenolic compounds. [63]. To determine the total antioxidant activity,
at position C5 is important for antioxidant activities. Secondly, this in bioactive compounds, the absorbance is measured against water and
hydrogen atom is involved in free radical inhibition activity because used to calculate the inhibition percent. The absorbance value is also
phenolic compounds inhibit free radicals by transfer of the hydrogen used to calculate total phenol content and total flavonoid [64]. In this
atom from its hydroxyl group [54]. study a single-point (SP) energy calculation provides the HOMO and
Water is a polar solvent whereas ethanol is less polar than water due LUMO energies (Table 1). All the flavonols have two or three charac­
to its non-polar hydrocarbon. Studies on the effect of solvent polarity on teristic absorption peaks in each medium (gas-phase, ethanol or water).
Each absorption peak is associated with an electron transition from the
ground state to its excited state and has a corresponding excitation
energy.
It has been reported that ethanol and water or a mixture of these
solvents is suitable since the hydroxyl groups are soluble in the solvents
[65–67]. We chose water as the solvent for the UV–vis calculation since
this is preferred medium for drug delivery. The maximum absorption
wavelength, λmax, of the flavonols in water and the excitation energies
are presented in the Supporting Information (Table S4). This study
shows that the maximum wavelength of absorption is inversely pro­
portional to the excitation energy (Fig. 6). Hence, the promotion of an
electron from the ground state to the excited state at longer wavelengths
require low energy, and vice versa. Gossypetin with the lowest excita­
tion energy and the highest wavelength, implies that the compound is
the most active flavonol. Azaleatin, Natsudaidain and 3-hydroxyflavone
have high excitation energies indicating low activity.
Time-dependent density functional theory (TD-DFT) calculations
Fig. 4. Atomic charge distribution in Gossypetin. gave information on electronic transitions such as which orbital was

6
E.O. Akintemi et al. Computational and Theoretical Chemistry 1210 (2022) 113658

Fig. 5. Molecular electrostatic potential (MEP) surface for the 15 flavonols.

involved in the transitions, extinction coefficient, oscillator strength,

excitation state energies and the electronic states of the transitions. The AB → AB*, AB** → A* + B*
strongest transitions for each molecule in the aqueous medium is pre­
sented (Table 3). where A* and B* are the resulting radicals. For flavonols, the resulting
Molar extinction coefficient, also known as molar absorptivity, is a radical is the flavonol radical (Ar*) and the hydroxyl radical (OH*)
measure of how strongly a chemical species attenuates light at a given [71,72]. The flavonol radical is very reactive. It binds to appropriate
wavelength [68]. All molecules have high molar absorptivity hence, sites to decrease lipid peroxidation and suppresses other scavenging
they strongly attenuate light. species.
Oscillator strength is a measure of the probability of excitation [69].
The strongest transitions (which are presented in Table 3) were chosen
based on the highest oscillator strengths compared to other transitions. 3.6. Molecular and physicochemical properties
Oscillator strength values >0.1 were obtained for all the flavonols
indicating the high probability of excited states. Ferry et al. [70] re­ The protocol of adsorption, distribution, metabolism and excretion
ported that electronic excitation produces molecules which may relax to (ADME) are very important considerations in the development of new
give radicals through homolytic bond splitting as follows: drug molecules. These are the parameters that play a critical part in
outlining the dose amounts, dosing intervals and overall safety of oral

7
E.O. Akintemi et al. Computational and Theoretical Chemistry 1210 (2022) 113658

Fig. 6. Relationship between excitation energy and maximum absorption wavelength.

Table 3
Strongest transitions of electronic absorptions in water.
Orbitals Extinction coefficient Oscillator strength 1st Excitation state energy (eV) Electronic state

3-Hydroxyflavone 62 → 63 0.69593 0.5853 4.0629 1– A′

Azaleatin 82 → 83 0.67399 0.6304 4.0893 1– A
Fisetin 74 → 75 0.68233 0.7263 3.9648 1– A′
Galagin 70 → 71 0.69162 0.5167 3.8835 1– A′
Gossypetin 82 → 83 0.69094 0.4697 3.6362 1– A′
Isorhamnetin 82 → 83 0.67397 0.5832 3.9387 1– A
Kaempferide 78 → 79 0.68557 0.5917 3.9739 1– A
Kaempferol 74 → 75 0.68985 0.6727 3.8482 1– A′
Morin 78 → 79 0.69057 0.7088 3.7732 1– A′
Myricetin 82 → 83 0.68062 0.6660 3.8173 1– A′
Natsudaidain 110 → 111 0.67673 0.5571 4.0803 1– A
Pachypodol 90 → 91 0.66378 0.6832 3.9156 1– A
Quercetin 78 → 79 0.68405 0.6886 3.8101 1– A′
Rhamnetin 82 → 83 0.68520 0.7309 3.7401 1– A
Rhamnazin 86 → 87 0.67131 0.5721 3.9819 1– A

drug administration [73]. Hence determining physicochemical proper­ and toxicology, as long-established by a large numeral of literature data
ties are vital [74] because these properties have an impact on various [79,80].
ADME attributes and can be easily controlled by medicinal chemists. In a study from Li et al. [81] they report that for active flavonols the
The primary objective is to find the desirable balance of physicochem­ log P value is principally 1.6 upward. So, values below this imply a
ical properties for safe drug administration [75]. reduced inhibitory activity. Gossypetin and Myricetin, which are
The logarithm of octanol–water partition coefficient (log POW) is a structural isomers with the molecular formula C15H10O8, have low
broadly acknowledged measure of lipophilicity and is determined for miLogP values (miLogP refers to the octanol–water partition coefficient
numerous compounds [76]. POW is the concentration ratio of the com­ as reported in Molinspiration) (1.42 and 1.39, respectively) (Table 4).
pound distributed between n-octanol and water [77,78]. Note: in this Despite the value implying reduced inhibitory activity they display
report, the subscript OW for the octanol–water partition coefficient, is strong chemical reactivity as already discussed.
omitted. The partition coefficient is an exceedingly important physico­ Vlahovic et al. [82] as well as Keen [83] and Shargel et al. [84] report
chemical parameter in medicinal chemistry, beneficial in pharmacology that hydrophobic substances have a high octanol–water partition

8
E.O. Akintemi et al. Computational and Theoretical Chemistry 1210 (2022) 113658

Table 4
Physicochemical parameters for flavonols.
miLogPa TPSAb Natomsc MWd nONe nOHNHf nviolationsg nrotbh Volumei

3-Hydroxyflavone 3.45 50.44 18 238.24 3 1 0 1 208.01

Azaleatin 1.96 120.36 23 316.26 7 4 0 2 257.61
Fisetin 1.97 111.12 21 286.24 6 4 0 1 232.07
Galagin 2.65 90.89 20 270.24 5 3 0 1 224.05
Gossypetin 1.42 151.58 23 318.24 8 6 1 1 248.10
Isorhamnetin 2.71 100.13 22 300.27 6 3 0 2 249.59
Kaempferide 2.17 111.12 21 286.24 6 4 0 1 232.07
Kaempferol 1.99 120.36 23 316.26 7 4 0 2 257.61
Morin 1.88 131.35 22 302.24 7 5 0 1 240.08
Myricetin 1.39 151.58 23 318.24 8 6 1 1 248.10
Natsudaidain 3.08 105.84 30 418.40 9 1 0 7 361.29
Pachypodol 2.80 98.37 25 344.32 7 2 0 4 292.67
Quercetin 1.68 131.35 22 302.24 7 5 0 1 240.08
Rhamnetin 2.53 109.36 24 330.29 7 3 0 3 275.14
Rhamnazin 2.22 120.36 23 316.26 7 4 0 2 257.61
a
Octanol–water partition coefficient as reported in Molinspiration.
b
Total polar surface area in Å2.
c
Total number of atoms.
d
Molecular Weight in g mol− 1.
e
Number of hydrogen bond acceptors.
f
Number of hydrogen bond donors.
g
Number of violations to Lipinski’s rule of 5.
h
Number of rotatable bonds.
i
Volume in in Å3. Formulae and total energy are presented in Table S2.

coefficient which is distributed primarily in the lipid bilayers of cells. active it should obey the following rules: number of hydrogen bond
Hydrophilic substances on the other hand have a low octanol–water donors (nOHNH) ≤ 5, number of hydrogen bond acceptors (nON) ≤ 10,
partition co-efficient and are distributed in blood plasma. They also molecular weight (MW) ≤ 500, partition coefficient (miLogP) ≤ 5, and
showed the linear relationship of the octanol–water partition coefficient number of violations (nviolations) ≤ 2 [91]. All the flavonols comply
with log P. This explains the high miLogP values obtained for Natsu­ with the Lipinski rule of 5 (Table 4), confirming that all are orally active
daidain and 3-hydroxyflavone (3.08 and 3.45, respectively) (Table 3). drug molecules. We found that the molecular weight of the flavonols is
Our deduction is that these two compounds are hydrophobic and are directly proportional to the volume of the compounds (Fig. 7). The
easily distributed in the lipid bilayers. On the other hand, Gossypetin Pearson correlation co-efficient (r) of 0.97 and R2 value of 0.94 verify
and Myricetin with low miLogP values (1.42 and 1.39 respectively) this.
(Table 4) are hydrophilic and will be readily distributed in blood For an average organic molecule, if the bioactivity score is more than
plasma. This also conforms to the high reactivity observed for Gossy­ 0 then it is active, a score ranging from − 0.50 to 0.0 indicates moder­
petin and Myricetin in water. ately active and a score of less than − 0.50 denotes inactivity [92]. The
results are presented in Table S5 (Supporting Information) but are
3.7. Drug likeness/bioactivity score depicted in Fig. 8.
All 15 compounds are excellently active (Fig. 8) as enzyme inhibitors
In drug discovery the critical phase is the hit-to-lead optimization. (EI). Although 3-hydroxyflavone (1) is a synthesized compound, it is a
Medicinal chemists use drug-likeness as the essential guide for drug strong enzyme inhibitor (0.31). The two natural flavonols that are
discovery [85,86]. The ADME parameters referred to as the pharmaco­ closest to this is Gossypetin (5) and Myricetin (10) which have enzyme
kinetic properties are also strongly affected by the physicochemical inhibitor values of 0.30. Galagin (4), Morin (9), and Quercetin (13) with
properties of a drug [87,88]. As already stated a drug-like molecule must an enzyme inhibitor values of 0.28 is also good. Flavonol (1) is inactive
be soluble in water as well as fat. This is necessary since after ingestion, as a nuclear receptor ligand (NRL). Morin (9) and Natsudaidain (11) are
an orally administered drug passes through the intestinal lining into the moderately active, while the remaining 9 are strongly active as kinase
bloodstream and is carried to the target. The molecule then penetrates inhibitors (KI). In addition, all demonstrate moderate activity as G-
the lipid-based cell membrane to get into the target cell [85]. Hydrogen protein-coupled receptors (GPCR), ion channel modulators (ICM) and
bonds have an important role in water solubility since to enable the
molecule to permeate into and through the lipid bilayer membrane, this
bond must break. Additionally, increased hydrogen bonds results in
increased aqueous solubility. The molecule is hence less hydrophobic,
which reduces partitioning [89].
Bioavailability refers to the rate and extent to which the flavonol is
absorbed from the fruit or vegetable and is available at the site of action.
Veber and co-workers [90] carried out oral bioavailability measure­
ments on over a thousand possible drug candidates and report that
compounds will have a high probability of good oral bioavailability if
they meet following criteria: (1) a total polar surface area (TPSA) ≤ 140
Å2 and, (2) if the number of rotatable bonds (nrotb) ≤ 10. All the fla­
vonols fulfil this criteria with the exception of Gossypetin and Myricetin
that have high total polar surface area (151.58 Å2) (Table 4) implying
they have reduced permeability in lipid membranes.
The Lipinski’s rule of 5 states that for a drug molecule to be orally Fig. 7. Relationship between volume and molecular weight.

9
E.O. Akintemi et al. Computational and Theoretical Chemistry 1210 (2022) 113658

Fig. 8. Profile depicting bioactivity score.

protease inhibitors (PI). References

4. Conclusion [1] E. Middleton, C. Kandaswami, T.C. Theoharides, The effects of plant flavonoids on
mammalian cells: implications for inflammation, heart diseases, and cancer,
Pharm. Rev. 52 (2000) 673–751.
All the flavonols investigated are chemically reactive. We have [2] O.M. Andersen, K.R. Markham, Flavonoids: Chemistry, Biochemistry and
shown that the highest negative charge for the oxygen atom on C-3′ Applications, CRC Press, 2005.
[3] D. Guajardo-Flores, S.O. Serna-Saldívar, J.A. Gutiérrez-Uribe, Evaluation of the
hydroxyl and highest positive charge on the H atom of C5-OH to be antioxidant and antiproliferative activities of extracted saponins and flavonols
crucial for antioxidant activity. The energy gap is lower in water and this from germinated black beans (Phaseolus vulgaris L.), Food Chem. 141 (2013)
is important because water is the medium for drug delivery. All the 1497–1503.
[4] M. Plaza, J. Kariuki, C. Turner, Quantification of individual phenolic compounds’
flavonols are reactive with Gossypetin having the lowest Egap. This also contribution to antioxidant capacity in apple: a novel analytical tool based on
highlights the hydrophilic property of flavonols. The potential distri­ liquid chromatography with diode array, electrochemical, and charged aerosol
butions on the MEP surfaces of the flavonols showed the hydroxyl and detection, J. Agric. Food Chem. 62 (2014) 409–418.
[5] X. Zhang, F. Chen, M. Wang, Antioxidant and antiglycation activity of selected
carbonyl groups are required for effective antioxidant properties.
dietary polyphenols in a cookie model, J. Agric. Food Chem. 62 (2014) 1643–1648.
Although the flavonols investigated exhibit moderate NLO properties, [6] T.C. Malig, M.R. Ashkin, A.L. Burman, M. Barday, B.J. Heyne, T.G. Back,
this value can be improved by derivatisation for application in materials Comparison of free-radical inhibiting antioxidant properties of carvedilol and its
science. phenolic metabolites, MedChemComm 8 (2017) 606–615.
[7] J.S. Wright, E.R. Johnson, G.A. DiLabio, Predicting the activity of phenolic
The quantum chemical calculations of their physicochemical prop­ antioxidants: theoretical method, analysis of substituent effects, and application to
erties showed they have good bioavailability. All the compounds agree major families of antioxidants, J. Am. Chem. Soc. 123 (2001) 1173–1183.
with the Lipinski rule of 5, indicating potential use as oral active drugs [8] M. Leopoldini, I.P. Pitarch, N. Russo, M. Toscano, Structure, conformation, and
electronic properties of apigenin, luteolin, and taxifolin antioxidants. A first
and promising drug likeness. All the flavonols are active as enzyme in­ principle theoretical study, J. Phys. Chem. A 108 (2004) 92–96.
hibitors and moderately active as G-protein-coupled receptors, ion [9] R. Praveena, K. Sadasivam, R. Kumaresan, V. Deepha, R. Sivakumar, Experimental
channel modulators and protease inhibitors. The results presented here and DFT studies on the antioxidant activity of a C-glycoside from Rhynchosia
capitata, Spectrochim. Acta A Mol. Biomol. Spectrosc. 103 (2013) 442–452.
can be used to further evaluate and test these flavonols in molecular [10] S.A. Aherne, N.M. O’Brien, Dietary flavonols: chemistry, food content, and
docking studies against infections and diseases. metabolism, Nutr. J. 18 (2002) 75–81.
[11] A. Crozier, M.E. Lean, M.S. McDonald, C. Black, Quantitative analysis of the
flavonoid content of commercial tomatoes, onions, lettuce, and celery, J. Agric.
Declaration of Competing Interest Food Chem. 45 (1997) 590–595.
[12] M.G. Hertog, P.C. Hollman, B. Van de Putte, Content of potentially
The authors declare that they have no known competing financial anticarcinogenic flavonoids of tea infusions, wines, and fruit juices, J. Agric. Food
Chem. 41 (1993) 1242–1246.
interests or personal relationships that could have appeared to influence
[13] M.G. Hertog, D. Kromhout, C. Aravanis, H. Blackburn, R. Buzina, F. Fidanza,
the work reported in this paper. S. Giampaoli, A. Jansen, A. Menotti, S. Nedeljkovic, Flavonoid intake and long-
term risk of coronary heart disease and cancer in the seven countries study, Arch.
Acknowledgement Intern. Med. 155 (1995) 381–386.
[14] M.S. McDonald, M. Hughes, J. Burns, M.E. Lean, D. Matthews, A. Crozier, Survey of
the free and conjugated myricetin and quercetin content of red wines of different
Thishana Singh acknowledges the National Research Foundation geographical origins, J. Agric. Food Chem. 46 (1998) 368–375.
(NRF) South Africa for the Post-PhD Thuthuka grant, the UKZN RO Fund [15] W. Bors, W. Heller, C. Michel, M. Saran, Flavonoids as antioxidants: determination
of radical-scavenging efficiencies, Meth. Enzymol. 186 (1990) 343–355.
and RSC Research Fund Grant (R20-2009) for financial support. All [16] A. von Gadow, E. Joubert, C.F. Hansmann, Comparison of the antioxidant activity
computations were carried out using the computational cluster re­ of aspalathin with that of other plant phenols of rooibos tea (Aspalathus linearis),
sources at the Centre for High Performance Computing (CHPC), Cape α-tocopherol, BHT, and BHA, J. Agric. Food Chem. 45 (1997) 632–638.
[17] M. Sato, N. Ramarathnam, Y. Suzuki, T. Ohkubo, M. Takeuchi, H. Ochi, Varietal
Town, South Africa. differences in the phenolic content and superoxide radical scavenging potential of
wines from different sources, J. Agric. Food Chem. 44 (1996) 37–41.
Appendix A. Supplementary material [18] S.A. Van Acker, L.M. Koymans, A. Bast, Molecular pharmacology of vitamin E:
structural aspects of antioxidant activity, Free Radic. Biol. Med. 15 (1993)
311–328.
Supplementary data to this article can be found online at https://doi. [19] E. Migliavacca, P.A. Carrupt, B. Testa, Theoretical parameters to characterize
org/10.1016/j.comptc.2022.113658. antioxidants. Part 1. The case of vitamin E and analogs, Helv. Chim. Acta 80 (1997)
1613–1626.
[20] R.A. Mendes, S.K. Almeida, I.N. Soares, C.A. Barboza, R.G. Freitas, A. Brown, G.
L. de Souza, A computational investigation on the antioxidant potential of
myricetin 3, 4′ -di-O-α-L-rhamnopyranoside, J. Mol. Model. 24 (2018) 1–8.

10
E.O. Akintemi et al. Computational and Theoretical Chemistry 1210 (2022) 113658

[21] A.S. Woldu, J. Mai, A novel relationship between the radical-scavenging activity of [46] N. Günay, H.A. Pir, D. Avci, V. Atalay, NLO and NBO analysis of sarcosine-maleic
flavonoids and enthalpy of formation revealed with Hartree-Fock computations acid by using HF and B3LYP calculations, J. Chem. 2013 (2013).
and thermochemical deduction, Redox Rep. 17 (2012) 115–130. [47] D. Avcı, Second and third-order nonlinear optical properties and molecular
[22] K.C. Weber, K.M. Honório, A.T. Bruni, A.D. Andricopulo, A.B. da Silva, A partial parameters of azo chromophores: Semiempirical analysis, Spectrochim. Acta Part A
least squares regression study with antioxidant flavonoid compounds, J. Struct. Mol. Biomol. Spectrosc. 82 (2011) 37–43.
Chem. 17 (2006) 307–313. [48] K. Bhat, K.J. Chang, M.D. Aggarwal, W.S. Wang, B.G. Penn, D.O. Frazier, Synthesis
[23] D. Mishurov, A. Voronkin, A.D. Roshal, O. Brovko, Relaxation behavior and and characterization of various Schiff bases for non-linear optical applications,
nonlinear properties of thermally stable polymers based on glycidyl derivatives of Mater. Chem. Phys. 44 (1996) 261–266.
quercetin, Opt. Mater. 57 (2016) 179–184. [49] A. Houlton, J.R. Miller, J. Silver, N. Jassim, M.J. Ahmet, T.L. Axon, D. Bloor, G.
[24] D. Mishurov, A. Voronkin, A.D. Roshal, Synthesis, molecular structure and optical H. Cross, Molecular materials for non-linear optics. Second harmonic generation
properties of glycidyl derivatives of quercetin, Struct. Chem. 27 (1) (2016) and the crystal and molecular structure of the 4-nitrophenylimine of
285–294. ferrocenecarboxaldehyde, Inorg. Chim. Acta 205 (1993) 67–70.
[25] D. Mishurov, A. Voronkin, O. Nedilko, I. Zykina, The influence of different factors [50] I. Ledoux, J. Zyss, Influence of the molecular environment in solution
on exploitation properties of nonlinear optical polymeric materials based on an measurements of the second-order optical susceptibility for urea and derivatives,
epoxy matrix doped with flavonoids, Polym. Test. 87 (2020), 106535. J. Chem. Phys. 73 (1982) 203–213.
[26] A. Kumar, A.K. Srivastava, R. Kumar, G. Brahmachari, N. Misra, Nonlinear Optical [51] P. Sankar, R. Philip, Nonlinear optical properties of nanomaterials, Mater. Charact.
Property Of 6’-Amino-5-Fluoro -2-Oxo-3’-Propyl-2’H-Spiro[Indoline-3,4’-Pyrano (2018) 301–334.
[2,3-C]Pyrazole]-5’-Carbonitrile- A Theoretical Approach, Organic and Medicinal [52] M. Karelson, V.S. Lobanov, A.R. Katritzky, Quantum-chemical descriptors in
Chem, IJ 3 (4) (2017), 555617. QSAR/QSPR studies, Chem. Rev. 96 (1996) 1027–1044.
[27] National Center for Biotechnology Information, https://pubchem.ncbi.nlm.nih. [53] B.W. Clare, C.T. Supuran, Carbonic Anhydrase Activators. 3: Structure-Activity
gov/#query (accessed October 2019). Correlations for a Series of Isozyme I1 Activators, J. Pharm. Sci. 83 (1994)
[28] R. Dennington, T. Keith, J. Millam, K. Eppinnett, W. Hovell, R.J.I. Gilliland, 768–773.
Shawnee Mission, KS, USA, GaussView. Version 6. Semichem, (2016). [54] N.F Santos-Sánchez, R. Salas-Coronado, R. Villanueva-Canongo, B. Hernandez-
[29] M.J. Frisch, G.W. Trucks, H.B. Schlegel, G.E. Scuseria, M.A. Robb, J.R. Cheeseman, Carlos, Antioxidant compounds and their antioxidant mechanism, in: Antioxidants.
G. Scalmani, V. Barone, G.A. Petersson, H. Nakatsuji, X. Li, M. Caricato, A.V. IntechOpen, London, UK, 2019, pp. 1–28.
Marenich, J. Bloino, B.G. Janesko, R. Gomperts, B. Mennucci, H.P. Hratchian, J.V. [55] M. Rahman, S. Hossain, A. Rahaman, N. Fatima, T. Nahar, B. Uddin, M.A. Basunia,
Ortiz, A.F. Izmaylov, J.L. Sonnenberg, Williams, F. Ding, F. Lipparini, F. Egidi, J. Antioxidant activity of Centella asiatica (Linn.) Urban: Impact of extraction solvent
Goings, B. Peng, A. Petrone, T. Henderson, D. Ranasinghe, V.G. Zakrzewski, J. Gao, polarity, J. Pharmacogn. Phytochem. 1 (2013).
N. Rega, G. Zheng, W. Liang, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. [56] F. Tian, B. Li, J. Yang, G. Zhang, Y. Chen, Y. Luo, Antioxidant and antimicrobial
Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, T. Vreven, K. activities of consecutive extracts from Galla chinensis: The polarity affects the
Throssell, J.A. Montgomery Jr., J.E. Peralta, F. Ogliaro, M.J. Bearpark, J.J. Heyd, bioactivities, Food Chem. 113 (2009) 173–179.
E.N. Brothers, K.N. Kudin, V.N. Staroverov, T.A. Keith, R. Kobayashi, J. Normand, [57] I. Alkorta, J.J. Perez, Molecular polarization potential maps of the nucleic acid
K. Raghavachari, A.P. Rendell, J.C. Burant, S.S. Iyengar, J. Tomasi, M. Cossi, J.M. bases, Int. J. Quantum Chem. 57 (1996) 123–135.
Millam, M. Klene, C. Adamo, R. Cammi, J.W. Ochterski, R.L. Martin, K. Morokuma, [58] C.U. Ibeji, G.F. Tolufashe, T. Ntombela, T. Govender, G.E. Maguire, G. Lamichhane,
O. Farkas, J.B. Foresman, D.J. Fox, Gaussian 16 Rev. B.01, Wallingford, CT, 2016. H.G. Kruger, B. Honarparvar, The catalytic role of water in the binding site of l, d-
[30] (a) C. Lee , W. Yang , R.G. Parr , Phys. Rev. B 37 (1988) 785–789 ; (b) A.D. Becke, transpeptidase 2 within acylation mechanism: a QM/MM (ONIOM) modelling,
Phys. Rev. A Mol. Opt. Phys. 38 (1988) 3098–3100 ; (c) R. Krishnan , J.S. Binkley , Tuberculosis 113 (2018) 222–230.
R. Seeger , J. Pople , J. Chem. Phys 72 (1980) 650–654 ; (d) T. Clark , J. [59] G.F. Tolufashe, V.T. Sabe, C.U. Ibeji, M.M. Lawal, T. Govender, G.E. Maguire,
Chandrasekhar , G.W. Spitznagel , P.V.R. Schleyer , J.Comput. Chem. 4 (1983) Inhibition mechanism of L, D-transpeptidase 5 in presence of the β-lactams using
294–301. ONIOM method, J. Mol. Graph. 87 (2019) 204–210.
[31] A.V. Marenich, C.J. Cramer, D.G. Truhlar, Universal solvation model based on [60] E.J. Braga, B.T. Corpe, M.M. Marinho, S.S. Marinho, Molecular electrostatic
solute electron density and on a continuum model of the solvent defined by the potential surface, HOMO–LUMO, and computational analysis of synthetic drug
bulk dielectric constant and atomic surface tensions, J. Phys. Chem. B 113 (2009) Rilpivirine, Int. J. Sci. Eng. Res. 7 (2016) 315–319.
6378–6396. [61] I.F. Afonso, Molecular modeling and evaluation of the structure-activity
[32] F. Trani, G. Scalmani, G.S. Zheng, I. Carnimeo, M.J. Frisch, V. Barone, Time- relationship coupled with in silico pharmacokinetic and toxicological studies of
dependent density functional tight binding: new formulation and benchmark of heterocyclic derivatives with antimicrobial activity, Dissertation, Federal
excited states, J. Chem. Theory Comput. 7 (2011) 3304–3313. University of Rio (2008).
[33] Molinspiration.com, Molinspiration Cheminformatics, https://www.molinspirati [62] S. Kumar, A.K. Pandey, Chemistry and biological activities of flavonoids: an
on.com/cgi-bin/properties, 2021 (Accessed 15 June 2021). overview, Sci. World J. (2013).
[34] R.G. Parr, R.A. Donnelly, M. Levy, W.E. Palke, Electronegativity: the density [63] M. Pisani, P. Astolfi, S. Sabbatini, P. Carloni, Antioxidant activity level, bioactive
functional viewpoint, J. Chem. Phys. 68 (1978) 3801–3807. compounds, colour and spectroscopic analysis (UV-Vis and FT-IR) of flavoured
[35] A.C. Gaudio, A. Korolkovas, Y. Takahata, Quantitative structure-activity drinks made with wine and sour cherries (Prunus cerasus Var. austera), Foods 10
relationships for 1, 4-dihydropyridine calcium channel antagonists (nifedipine (8) (2021) 1953.
analogues): A quantum chemical/classical approach, J. Pharm. Sci. 83 (1994) [64] S. Donkor, C. Larbie, G. Komlaga, B.O. Emikpe, Phytochemical, antimicrobial, and
1110–1115. antioxidant profiles of duranta erecta L. parts, Biochem. Res. Int. 2019 (2019)
[36] P.K. Chattaraj, P. Perez, J. Zevallos, A. Toro-Labbe, Theoretical study of the trans- 8731595.
N2H2→ cis-N2H2 and F2S2→ FSSF reactions in gas and solution phases, J. Mol. [65] R.T. Ramos, I.C. Bezerra, M.R. Ferreira, L.A. Lira, Dual Anti-cholinesterase activity
Struct.-THEOCHEM 580 (2002) 171–182. of Ajoene by in silico and In Vitro studies, Res. J. Pharmacogn. 9 (2017) 253.
[37] S.K. Nasiri, A. Reisi-Vanani, M. Hamadanian, Molecular Structure, Spectroscopic, [66] Y. Chen, J. Wang, D. Wan, Determination of total flavonoids in three Sedum crude
Local and Global Reactivity Descriptors and NBO Analysis of C32H12: A New drugs by UV–Vis spectrophotometry, Pharmacogn. Mag. 6 (2010) 259.
Buckybowl and Sub-Fullerene Structure, Polycyclic Aromat. Compd. 1 (2018) [67] H. Zhu, Y. Wang, Y. Liu, Y. Xia, T. Tang, Analysis of flavonoids in Portulaca
1–12. oleracea L. by UV–vis spectrophotometry with comparative study on different
[38] P.K. Chattaraj, S. Giri, Electrophilicity index within a conceptual DFT framework, extraction technologies, Food Anal. Methods 3 (2010) 90–97.
Annu. Rep. Prog. Chem. Sect. C: Phys. Chem. 105 (2009) 13–39. [68] A.D. McNaught, A. Wilkinson, Compendium of Chemical Terminology, 1669,
[39] J. Padmanabhan, R. Parthasarathi, M. Elango, V. Subramanian, B. Blackwell Science Oxford, 1997.
S. Krishnamoorthy, S. Gutierrez-Oliva, A. Toro-Labbe, D.R. Roy, P.K. Chattaraj, [69] R.L. Martin, Natural transition orbitals, J. Chem. Phys. 118 (2003) 4775–4777.
Multiphilic descriptor for chemical reactivity and selectivity, J. Phys. Chem. A 111 [70] M. Ferry, Y. Ngono-Ravache, C. Aymes-Chodur, M.C. Clochard, X. Coqueret,
(2007) 9130–9138. L. Cortella, E. Pellizzi, S. Rouif, S. Esnouf, Ionizing Radiation Effects in Polymers,
[40] Y.L. Zhang, Nucleophilic reactivity analysis: atom-bond electronegativity Elsevier, 2016.
equalization method Yang, J. Mol. Struct.-THEOCHEM 496 (2000) 139–144. [71] J. Treml, K. Šmejkal, Flavonoids as potent scavengers of hydroxyl radicals, Compr.
[41] R. Shankar, K. Senthilkumar, P. Kolandaivel, Calculation of ionization potential Rev. Food Sci. Food Saf. 15 (2016) 720–738.
and chemical hardness: a comparative study of different methods, Int. J. Quantum [72] A. Pérez-González, E. García-Hernández, E. Chigo-Anota, The antioxidant capacity
Chem. 109 (2009) 764–771. of an imidazole alkaloids family through single-electron transfer reactions, J. Mol.
[42] P. Geerlings, F. De Proft, P. Ayers, Chemical reactivity and the shape function, Model. 26 (11) (2020) 321.
Comput. Theor. Chem. 1 (2007) 1–17. [73] D.F. McGinnity, J. Collington, R.P. Austin, R.J. Riley, Evaluation of human
[43] O.E. Oyeneyin, I.A. Adejoro, B.T. Ogunyemi, O.T. Esan, Structural and solvent pharmacokinetics, therapeutic dose and exposure predictions using marketed oral
dependence on the molecular and nonlinear optical properties of 10-octyl drugs, Curr. Drug Metab. 8 (2007) 463–479.
thiophene-based phenothiazine and substituted derivatives–a theoretical approach, [74] S. Nithiyanantham, L. Palaniappan, M. Lenin, Physico-chemical and bio-chemical
J. Taibah Univ. Sci. 12 (2018) 483–493. nature of pectin with amylase S, J. Bionanosci. 9 (2015), 359362.
[44] O.E. Oyeneyin, A. Ajibade, O.T. Esan, Substituent Effects on the Structural and [75] S. Belaidi, H. Belaidi, A. Kerassa, M. Saoula, D. Bouzidi, Predictive qualitative
Nonlinear Optical Properties of 1-[4-({(E)-[4-(methylsulfanyl) phenyl] structure-property/activity relationships for drug design in some of
methylidene} amino) phenyl] ethanone and Some of its Substituted Derivatives-a antimycobacterial pyrrole derivatives, Quantum Matter 5 (2016) 798–805.
Theoretical Method, Phys. Chem. Res. 6 (2018) 667–683. [76] S.K. Poole, C.F. Poole, Separation methods for estimating octanol–water partition
[45] O. Oyeneyin, Structural and solvent dependence of the electronic properties and coefficients, J. Chromatogr. B Biomed. Appl. 797 (2003) 3–19.
corrosion inhibitive potentials of 1, 3, 4-thiadiazole and Its substituted derivatives-
A theoretical investigation, Phys. Sci. Int. J. (2017) 1–8.

11
E.O. Akintemi et al. Computational and Theoretical Chemistry 1210 (2022) 113658

[77] R.F. Ribeiro, A.V. Marenich, C.J. Cramer, D.G. Truhlar, The solvation, partitioning, [86] H. Zheng, J. Jie, Y. Zheng, Multi-Swarm Chaotic Particle Swarm Optimization for
hydrogen bonding, and dimerization of nucleotide bases: a multifaceted challenge Protein Folding, J. Bionanosci. 7 (2013) 643–648.
for quantum chemistry, Phys. Chem. Chem. Phys. 13 (2011) 10908–10922. [87] B. Wei, Z. Shi, J. Xiao, Y. Xu, L. Lv, In vivo and in vitro antibacterial effect of nano-
[78] A. Jalan, R.W. Ashcraft, R.H. West, W.H. Green, Predicting solvation energies for structured titanium coating incorporated with silver oxide nanoparticles,
kinetic modelling, Annu. Rep. Prog. Chem. Sect. C: Phys. Chem. 106 (2010) J. Biomater. Tissue Eng. 7 (5) (2017) 418–425.
211–258. [88] H.G. Guo, H.S. Guo, F.L. Yao, S.G. Yang, Z. Chen, T. Wang, In vitro study on
[79] S. Balaz, Modeling kinetics of subcellular disposition of chemicals, Chem. Rev. 109 biological potential of tissue-engineered cage by using surface-modified technique:
(2009) 1793–1899. A preliminary evaluation, J. Biomater. Tissue Eng. 6 (2) (2016) 114–121.
[80] J.M. Sangster, Octanol-water Partition Coefficients: Fundamentals and Physical [89] L. Di, E.H. Kerns, Drug-Like Properties: Concepts, Structure Design and Methods
Chemistry, vol. 1, John Wiley & Sons, 1997. from ADME to Toxicity Optimization, Academic Press, USA, 2015.
[81] B.H. Li, W.X. Tian, Inhibitory effects of flavonoids on animal fatty acid synthase, [90] D.F. Veber, S.R. Johnson, H.Y. Cheng, B.R. Smith, K.W. Ward, K.D. Kopple,
J. Biochem. 135 (2004) 85–91. Molecular properties that influence the oral bioavailability of drug candidates,
[82] F. Vlahović, S. Ivahović, M. Zlatar, M. Gruden, Density functional theory J. Med. Chem. 45 (12) (2002) 2615–2623.
calculation of lipophilicity for organophosphate type pesticides, J. Serbian Chem. [91] V.K. Maurya, S. Kumar, A.K. Prasad, M.L. Bhatt, S.K. Saxena, Structure-based drug
Soc. 82 (2017) 1369–1378. designing for potential antiviral activity of selected natural products from
[83] P. Keen, Effect of binding to plasma proteins on the distribution, activity and Ayurveda against SARS-CoV-2 spike glycoprotein and its cellular receptor,
elimination of drugs, Biochem, Pharmacol. Springer (1971) 213–233. Virusdisease 31 (2020) 179–193.
[84] L. Shargel, B.C. Andrew, S. Wu-Pong, Applied Biopharmaceutics & [92] E.A. Alodeani, M. Arshad, M.A. Izhari, Anti-uropathogenic activity, drug likeness,
Pharmacokinetics, vol. 246, Appleton & Lange, Stamford, 1999. physicochemical and molecular docking assessment of (E-)-N′ -(substituted-
[85] J.M. Beale, J. Block, R. Hill, Organic Medicinal and Pharmaceutical Chemistry, benzylidene)-2-(quinolin-8-yloxy) acetohydrazide, Asian Pac. J. Trop. Biomed. 5
Lippincott Williams and Wilkins, New York, Philadelphia, 2010. (2015) 676–683.

12