key: cord-0071521-gmdpehy2
authors: Bitew, Mamaru; Desalegn, Tegene; Demissie, Taye B.; Belayneh, Anteneh; Endale, Milkyas; Eswaramoorthy, Rajalakshmanan
title: Pharmacokinetics and drug-likeness of antidiabetic flavonoids: Molecular docking and DFT study
date: 2021-12-10
journal: PLoS One
DOI: 10.1371/journal.pone.0260853
sha: 90349b269af1387fd133fa67e7582aa7bfeba007
doc_id: 71521
cord_uid: gmdpehy2

Computer aided toxicity and pharmacokinetic prediction studies attracted the attention of pharmaceutical industries as an alternative means to predict potential drug candidates. In the present study, in-silico pharmacokinetic properties (ADME), drug-likeness, toxicity profiles of sixteen antidiabetic flavonoids that have ideal bidentate chelating sites for metal ion coordination were examined using SwissADME, Pro Tox II, vNN and ADMETlab web tools. Density functional theory (DFT) calculations were also employed to calculate quantum chemical descriptors of the compounds. Molecular docking studies against human alpha amylase were also conducted. The results were compared with the control drugs, metformin and acarbose. The drug-likeness prediction results showed that all flavonoids, except myricetin, were found to obey Lipinski’s rule of five for their drug like molecular nature. Pharmacokinetically, chrysin, wogonin, genistein, baicalein, and apigenin showed best absorption profile with human intestinal absorption (HIA) value of ≥ 30%, compared to the other flavonoids. Baicalein, butein, ellagic acid, eriodyctiol, Fisetin and quercetin were predicted to show carcinogenicity. The flavonoid derivatives considered in this study are predicted to be suitable molecules for CYP3A probes, except eriodyctiol which interacts with P-glycoprotein (p-gp). The toxicological endpoints prediction analysis showed that the median lethal dose (LD(50)) values range from 159–3919 mg/Kg, of which baicalein and quercetin are found to be mutagenic whereas butein is found to be the only immunotoxin. Molecular docking studies showed that the significant interaction (-7.5 to -8.3 kcal/mol) of the studied molecules in the binding pocket of the α-amylase protein relative to the control metformin with the crucial amino acids Asp 197, Glu 233, Asp 197, Glu 233, Trp 59, Tyr 62, His 101, Leu 162, Arg 195, His 299 and Leu 165. Chrysin was predicted to be a ligand with high absorption and lipophilicity with 84.6% absorption compared to metformin (78.3%). Moreover, quantum chemical, ADMET, drug-likeness and molecular docking profiles predicted that chrysin is a good bidentate ligand.

a1111111111 a1111111111 a1111111111 a1111111111 a1111111111

Computational pharmacology is a fast-growing research field focusing on the development of techniques for employing software and databases to generate and analyze molecular, biological and medical data from diverse sources [1, 2] . The methods have been used since the starting of a drug discovery to screen and identify new lead compounds available in molecular libraries [3] [4] [5] , enabling parallel optimization of compound efficacy, drug-like properties and thereby playing a role in improving the quality of drug candidates [6] . The design and development of drug molecules requires earlier assessment of pharmacokinetic properties, absorption, distribution, metabolism and excretion (ADME), at a stage where multiple compounds are considered, but access to the physical samples is limited [3, 7] . Due to limited pharmacokinetics and toxicity profile of new drug-like molecules, only 11% of the drugs succeed to enter clinical development stage to reach the market [4, 8] . As a result, the pharmaceutical industries remain under immense pressure to counter the high rate of attrition in drug development that triggered an increase in the interest of computer aided toxicity and pharmacokinetic profile predictions [9] . It is evident that computer aided methods save time, human and material resources in the search of new drug-like molecules in contrast to experimental techniques [10, 11] . As a result, there has been marvelous advancement in the area of computer aided drug design and computational chemistry. These methods have been used for screening of new chemical species and their chemical properties [12] . On the other hand, web-based platforms-ADMETlab, Pro Tox-II and SwissADME-have been developed to predict ADMET properties [3, 4, 13] .

Diabetes mellitus (DM) is a group of metabolic conditions in which high blood sugar levels transpire over an extended period [14] [15] [16] [17] . Modern medical treatment methods for diabetes including injection and oral hypoglycemic drugs have some toxic or side effects, economic pressures, and so on [14, 18] . Medicinal plants have different classes of phytochemicals with bioactive compounds of various therapeutic and pharmacological potentials [19, 20] , of which polyphenols are among the most studied compounds [21, 22] . The presence of bioactive phytochemical compounds such as phenols and flavonoids have been used to confirm the therapeutic usefulness of a given plant [23, 24] . Their bioactivity is originated from their structural features that give them strong tendency to scavenge free radicals, reduce oxidized chemical entities, bring metal chelation, and inhibit enzymes [25] .

Metal based organic compounds that are synthesized using phytochemicals as ligands are serving as metalo-therapeutics, and they are emerging strategies to treat type II diabetes mellitus. They are known to significantly improve hyperglycemia, glucose intolerance, insulin resistance, and decrease level of serum triglyceride and plasma glucose [26, 27] . Previous studies indicated that metal-flavonoid complexes have better antioxidant activity than the free ligands alone. This is due to the attainment of extra resonating site and better stability of flavonoid metal complex radical than free flavonoids possess [16, 28] [40, 41] .

In the present study, in-silico pharmacokinetic (ADME) properties, drug-likeness, toxicity profiles and toxicological endpoints of sixteen antidiabetic flavonoids that have ideal bidentate chelating sites for metal ion coordination were examined using SwissADME, Pro Tox II, vNN and ADMETlab web tools, Molecular docking and DFT calculations.

The overall procedures followed for the ligand screening process is shown in Fig 1, whereas the structures of the flavonoids studied, which are of plant origin, are presented in Fig 2. The molecules which are structurally similar with compound 9, having a bidentate coordination cite and showed antidiabetic activity in-vitro, were retrieved from literature [30, 42] . The SwissADME [3] web tool from Swiss Institute of Bioinformatics (SIB) was used to convert the two dimensional structures into their simplified molecular input line entry system (SMILES) [43] .

The molecular docking study of the target flavonoids was made by following the reported method. The 3D structure of human pancreatic α-amylase was obtained from Protein Data Bank [44] 

Pharmaceutically significant descriptors and physically relevant properties of the target flavonoids were predicted [45, 46] . The physicochemical properties (molar refractivity, topological polar surface area, number of hydrogen bond donors/number of hydrogen bond acceptors), lipophilicity (logP O/w ), pharmacokinetics properties (GI absorption, BBB permeation, P-gp substrate, cytochrome-P enzymes inhibition, skin permeation (log Kp)) which are critical parameters for prediction of the absorption and distribution of drugs within the body [12] , and drug likeness (Lipinski's rule of five) were predicted via SwissADME [3] . The percent absorption (% Abs) of the ligands were calculated using the reported formula % Abs = 109-0.345 TPSA [47] . The ADMETlab [4] web server was used to predict physicochemical properties (logS, logD, logP); Absorption: CaCO 3 Permeability (CaCO2), P-gp inhibitor/substrate (Pgp), Human Intestinal Absorption (HIA), 20% bioavailability (F20%), and 30% bioavailability (F30%)); Distribution: Plasma Protein Binding (PPB), Blood Brain Barrier (BBB) and Volume Distribution (VD) and Excretion: Half Life (T 1/2 ) and Clearance (CI). The properties of CaCO 2 , VD, PPB, CI, T 1/2 were expressed numerically, whereas the rest of the pharmacokinetic parameters were expressed categorically [4] . Metabolism related parameter, human liver microsomal (HLM) assay was predicted using variable nearest neighbor (vNN) method [5] . Besides, water solubility, LogP and LogS values of the flavonoid derivatives were retrieved from drug bank [48] , which are predicted using ALOGPS 2.1 and ChemAxon. The molecular descriptors and properties generated from all webtool/webservers were considered to screen the best candidate from the sixteen flavonoids [12] .

A total of 19 parameters were predicted to study the toxicity profile of the sixteen flavonoids and the two controls. The toxicological endpoints (Hepatotoxicity, Carcinogenicity, Immunotoxicity, Mutagenicity) and the level of toxicity (LD 50 , mg/Kg) of the studied flavonoids were determined using ProTox-II server [13, 49] . Toxicity profiles-the hERG potassium channel inhibition (cardiotoxicity), H-HT (Human Hepatotoxicity), AMES (Ames Mutagenicity) and SkinSen (Skin sensitization)-were predicted using ADMETlab [4] . Drug induced liver injury (DILI), mitochondrial membrane potential (MMP) toxicity and cytotoxicity parameters were predicted using the vNN method [5] .

Geometry optimizations and frequency calculations were performed using the Gaussian 09 [50] program package and the results were visualized using Gaussview 06 software. Geometries of the compounds for the density functional theory (DFT) calculations in which the DFT/ B3LYP [51] [52] [53] hybrid functional was used together with 6-311++G(d, p) basis set [54] . Such an approach has been shown to give reliable results for related organic molecules [55] [56] [57] . The optimized geometries were confirmed to be real minima without any imaginary vibrational frequency by performing vibrational frequency calculations at the same level of theory. Quantum chemical descriptor set was chosen to obtain pertinent electronic and structural features [45] . The frontier molecular orbitals (FMOs), energy gap (ΔE = E LUMO -E HOMO ), electronegativity (χ = -½ (E HOMO + E LUMO )), electronic chemical potential (μ = ½ (E HOMO + E LUMO ) = -χ), global chemical hardness (η = ½ ( ELUMO -E HOMO )), global softness (σ = 1/2η), global electrophilicity index (ω = μ 2 /2η), and nucleophilicity index (Nu = 1/ω), and dipole moment were also calculated and analyzed at the same level of theory [15, 58] .

Out of the sixteen flavonoids listed in Fig 2, fifteen of them do not violate Lipinski's rule of five [59] for oral availability, demonstrating that the compounds have drug-like molecular nature (Table 1 ) except myricetin (16) holding six NHBDs (> 5). It is important to note that Lipinski's rule of five is essential for rational drug design and it has been suggested that the low permeability or poor absorption for a given compounds results when it violates one of Lipinski's rule of five [12] . Molecular weight, LogP and number of hydrogen bond acceptors (NHBAs) of all flavonoids are within the recognized values of less than 500, 3 and 10, respectively. It has been predicted that the molecules that have similar molecular mass (Apigenin (1) and Baicalein (2); Fisetin (7) and Kaempferol (8); Isorhamnetin (11), Luteolin (12) , and Naringenin (14) ) showed similarities in TPSA, %Abs and MR ( Table 1 ). The molecules with a TPSA of 140 Å 2 and above would be poorly absorbed (< 10% fractional absorption), while those with a TPSA 60 Å 2 would be well absorbed (> 90%) [60] . The TPSA and %Abn of the flavonoid derivatives showed the preference of apigenin (1), baicalein (3), chrysin (4), genistein (10) , and wogonin (15) with better and comparable percent absorption to the control metformin (18). It is found that the %Abn and TPSA values follow the order chrysin (4) > wogonin (15) > metformin (18) > genistein (10) = apigenin (1) = baicalein (2) .

The results also showed an inverse relationship between TPSA and %Abs, and direct relationship between %Abs and MR. The analysis of the %Abs of the sixteen studied flavonoids and the two controls for their percent absorption (Table 1) showed that chrysin (4) have the highest percent absorption (84.6%) than all the other flavonoids and the two controls, whereas metformin (18) displayed better percent absorption (78.3%) than acarbose. These results are in good agreement with a previous report on pharmacokinetic properties of apigenin (1), eriodyctiol (6), fisetin (7), kaempferol (8) and quercetin (9) via molinspiration and pkCSM servers [61] . The calculated result of LogS showed that chrysin (4) (-3.709) is the least soluble molecule and ellagic acid (5) (-3.076) has the highest solubility. The relatively high solubility of ellagic acid (5) is due to the presence of four hydroxyl and two ketone pharmacophores [46] relative to the other flavonoids studied in this work. The computed results for the distribution coefficient (logD 7.4 ) of the flavonoids (Table 1) showed that they have high solubility with low permeability by passive transcellular diffusion [62] . It is predicted that acarbose (17) has moderate solubility, moderate permeability and low metabolism with logD 7.4 value of 2.543 and very low lipophilicity of -8.565.

The distribution coefficients (logP) of the flavonoids (Tables 1 and 2 ) range between 1.313 and 2.871. The results of the logP values showed that they have optimal lipophilicity (Optimal: 0 < LogP < 3). The logP results for the two controls, acarbose (17) and metformin (18), are less than zero, -1.034 and -8.565, respectively, suggesting that the two compounds have poor membrane permeability. The logP of apigenin (1) (3.07) and chrysin (4) (3.44) were found to be greater than three as predicted in ALOGPS (Table 2 ) which might be due to the small number of hydroxyl groups present in both compounds.

The water solubility (mg/mL) of some of the flavonoids considered in this work was retrieved from the drug bank (predicted using ALOGPS) and tabulated in Table 2 . It is predicted that ellagic acid (5) has high solubility (0.823 mg/mL), whereas chrysin (4) is found to have low solubility (0.105 mg/mL). These results are in agreement with those predicted using Table 2 . Comparison of solubility, partition coefficient, water solubility, and P Ka of flavonoids predicted using different methods. ADMETlab as well as with the experimental findings suggesting that chrysin (4) has low aqueous solubility and high absorption or absorption enhancer role [63] . The pKa results range from 5.54 to 7.91 indicating that the studied flavonoids fall under the category of neutral organic acids, in good agreement with experimental pKa of apigenin (1), ellagic acid (5) and quercetin (9) (7.12, 6.69, and 7.17, respectively). The relatively high pKa of naringenin (14) (7.91), apigenin (1) (6.63), chrysin (4) (6.58), luteolin (12) (6.57) and genistein (10) (6.55) may cause difficulty of deprotonation, suggesting that the ligands have stronger coordination ability with metals [64] .

On the other hand, the lipophilicity (LogP o/w ) study showed that wogonin (15) (2.880), chrysin (4) (2.871), genistein (10) (2.557), apigenin (1) and 2 (2.577) are highly lipophilic than the other molecules, whereas quercetin (9) and morin (13) (1.988) were found to be less lipophilic ligands next to ellagic acid (5) (1.313) ( Table 2 ). The presence of one methoxy group in the benzoyl part of wogonin (15) makes it more lipophilic than chrysin (4) . It has also been found that the structural resemblance in benzoyl or cinnamoyl part of the flavonoid derivatives play a role in their physicochemical properties. The calculated LogP o/w of the flavonoid derivatives considered in this work is found to range between 1.313 and 2.880. Our predicted results are in good agreement with the previously reported P o/w values of 1.92, 2.27, 3.11 and 1.82 for apigenin (1), eriodyctiol (6), fisetin (7) and kaempferol (8), respectively [3] . Although further clinical study is required, the predicted results suggest that oral and intestinal absorptions are possible.

To evaluate the absorption property of the flavonoids, we predicted Caco-2-Permeability (Papp), P-glycoprotein (P-gp), substrate or inhibitor (P-gpinh/P-gpsub), Human intestinal absorption (HIA), 20% bioavailability (F 20 ), and 30% bioavailability (F 30 ) as presented (Table 3 ). It is found that the values of Caco-2-Permeability are in the range from -4.732 to -6.63 cm/s, suggesting they have low to moderate permeability. The high negative Caco-2-Permeability (cm/s) of myricetin (16) (-6.63 cm/s), morin (13) (-6.206 cm/s), kaempferol (8) (-6.163 cm/s) and ellagic acid (5) (-6.001 cm/s) is mainly due to the presence of large number of NHBDs and NHBAs in their structures which have the potential to rendere their permeability. The two controls, acarbose (17) and metformin (18) are predicted to have Caco-2 permeability values of -5.803 cm/s and -7.128 cm/s, respectively. The predicted HIA data of the sixteen flavonoids and the two control drugs are represented categorically as 0 and 1 ( Table 3 ) which showed that apigenin (1), chrysin (4), fisetin (7), genistein (10), wogonin (15) , and metformin (18) have high HIA (� 30%) whereas acarbose (17) and the rest six flavonoids are predicted to have HIA below 30% (Table 3) . Flavonoids baicalein (2), chrysin (4), genistein (10), wogonin (15) and the control, metformin (18) predicted to show bioavailability greater than 30%, whereas the rest of the flavonoids are predicted to show bioavailability between 20 to 30%. Acarbose (17) is known for its extremely low bioavailability and less than 2% of absorption as active drug [48] . We predicted low bioavailability (<20%, Table 3 ) and extremely low percent absorption of acarbose (-1.8%, Table 1 ).

Plasma Protein Binding (PPB), Blood Brain Barrier (BBB) permeability, and Volume of distribution (VD) were parameters considered to study the distribution of drug candidates ( Table 4 ). The PPB of apigenin (1), eriodyctiol (5), fisetin (7), kaempferol (8) and quercetin (9) are predicted to be above 90% to bind to common blood proteins significantly, indicating that the pharmacologic effect of these compounds is low. The control drugs, acarbose (17) and metformin (18) weakly bind to common blood proteins. This could be one of the reasons that make acarbose as a working drug, while it has very low percent absorption (-1.8%).

Baicalein (2), chrysin (4), ellagic acid (5), eriodyctiol (6), wogonin (15) and metformin (18) were predicted to permeate BBB using ADMETlab. However, our prediction using Swis-sADME webtool showed that only chrysin (4) can permeate BBB. On contrary to previous report on ellagic acid (5) studied using pkCSM server [6] , our study showed positive result predicted by ADMETlab for ellagic acid (5) to permeate BBB.

The volume of distribution (VD) in L/Kg were predicted to be all negative, indicating higher distribution of the flavonoid derivatives and the two control drugs in plasma than in tissues, confined to plasma proteins. The hydrophilic character of the flavonoids contributed to low values of VD. The more negative value of VD for quercetin (-1.56) and less negative value for wogonin (-0.32) followed by chrysin (-0.35) among flavonoids were also presented in Table 3 . The high VD values of wogonin (15) and chrysin (4)-relative to the other flavonoids -showed the consistency of absorption and distribution parameters.

The skin permeability, Kp, values obtained for all the flavonoid derivatives considered in this study are in the range of -7.4 to -5.35 cm/s (Table 4 ) inferring low skin permeability [3] . It is found that log Kp increases (more negative) with the increase in molecular size ( Table 4 ). The results indicated that eriodyctiol (6) is the only flavonoid to act as drug transporter of P-glycoprotein substrate from target flavonoids considered (Table 4 ). This may lead to pharmacokinetics-related drug-drug interactions leading to toxic or other unwanted adverse effects due to the lower clearance and accumulation of the drug or its metabolites [65] . In contrary, it is the only molecule that can act as CYP1A2 inhibitor (Table 4 ). Its interaction with p-gp together with its carcinogenicity will lead eriodyctiol (6) to be excluded from the list of bidentate ligands for metal complexation. Except ellagic acid (5), the rest of the flavonoid derivatives considered in this study are predicted to interact with CYP3A, making them very suitable molecules as CYP3A4 probes (Table 4 ). CYP3A enzymes are important determinants in the therapeutic efficacy and toxicity of numerous drugs. The interactions at the CYP3A level are often the cause of pronounced drug-drug interactions [66] . Baicalein (2), chrysin (4), fisetin (7), kaempferol (8), quercetin (9), genistein (10), isorhamnetin (11), luteolin (12), morin (15), naringenin (14) and wogonin (15) are predicted as D6 inhibitors. Previous study by May & Schindler, [67] showed that the interaction of antidiabetic drugs with other drugs in CYP3A4 and CYP2D6 inhibition mechanism has an additive effect for antidiabetic drugs [67] . In the present work, all the flavonoid derivatives except ellagic acid (5) were found to inhibit CYP3A4. These results suggest that the flavonoids studied in this work potentially have antidiabetic effect in agreement with previous reports [42, 68] . In this study, only butein (3) and wogonin (15) are predicted to interact with 2C9. Previous report indicated that the inhibition and induction of the CYP enzymes are probably the common causes for most drug interactions [65] . Hence, the predicted carcinogenicity and immunotoxicity (Table 5 ) together with its interaction with 2C9 would lead butein as a flavonoid derivative with high safety risk. In general, our predicted results could be used as important inputs for further experimental studies.

The human liver microsomal (HLM) stability and mitochondrial membrane potential (MMP) parameters of the flavonoid derivatives predicted via the vNN method (Table 4) showed that most of the flavonoid derivatives showed mitochondrial toxicity. Ellagic acid (5) and acarbose (17) were predicted to be inactive for MMP. Apigenin (1) was active for HLM prediction, and could be metabolized rapidly by the liver that will reduce its effective therapeutic concentration in the body [5] .

Absence of toxicity is one of the major factors for selecting a compound as a therapeutic candidate [69] . Herein, the organ toxicity (hepatotoxicity) and toxicological endpoints (carcinogenicity, immunotoxicity, mutagenicity and cytotoxicity) of the sixteen flavonoid derivatives were predicted. The predicted descriptors are presented in Table 5 . The results show that all of the compounds are not hepatotoxic and cytotoxic. Of the considered flavonoid derivatives, baicalein (2), myricetin (16), luteolin (12) and quercetin (9) and butein (3) show two toxicological endpoints; carcinogenicity and mutagenicity, and carcinogenicity and immunotoxicity, respectively. On the other hand, ellagic acid (5), eriodyctiol (6), fisetin (7) and morin (15) show one toxicological endpoint (carcinogenicity). Apigenin (1), chrysin (4), kaempferol (8), wogonin (15) , naringenin (14) , isorhamnetin (11) , and genistein (10) are free from any of the predicted toxicological endpoints. Acute toxicity predictions indicate that the flavonoids considered are not fetal; however, they can be categorized as harmful toxic classes.

The median lethal dose (LD 50 ) values were found to be in the range from 159-3919 mg/Kg (Table 5 ). According to the globally harmonized system of classification of labeling of chemicals (as described in Pro Tox II), quercetin (9), fisetin (7) , and myricetin (16) are predicted as toxic (Class III), whereas apigenin (1), butein (3), elagic acid (5) and eriodyctiol (6) are harmful (Class IV) and the rest were categorized as "may be harmful" (Class V) toxicity classes. The results show that most of the drug-like compounds have the tendency to have carcinogenic activities than any toxicological endpoints.

Besides the ADMETlab descriptors, vNN descriptors such as hERG Blockers (hERG), human hepatotoxicity (H-HT), Ames Mutagenicity (AMES), Skin sensitization (SkinSen) were used to predict the toxicity of the flavonoid derivatives and the two controls. The results are presented in Table 6 . The numbers 1 and 0 in the category section (Table 6 ) stand for positive and negative, respectively, for the predicted toxicities. Except apigenin (1), baicalein (2) and genistein (10) , the rest of the flavonoid derivatives and the two control drugs were not predicted to be human ether-a-go-go related gene (hERG) blockers. However, most of them were predicted to show H-HT and mutagenicity ( Table 6 ). The vNN predicted drug induced liver injury (DILI) and cytotoxicity results (Table 6 ) showed that most of the flavonoid derivatives are not cytotoxic which is in line with results from Pro Tox II prediction. Apigenin (1) and acarbose (17) were predicted with inactive and active, respectively, to induce liver injury. The results also indicate the tendency of acarbose (17) to cause DILI (Table 6 ), which is in agreement with previous report on the risk of hepatotoxicity of acarbose (17) [70] and the probability of acarbose-induced hepatotoxicity [71] .

The frontier molecular orbital analysis (HOMO-LUMO) was performed to predict the reactivity of the flavonoid derivatives. DFT chemical descriptors: electronegativity (χ), electronic chemical potential (μ), globalchemical hardness (η), global softness (σ), global electrophilicity index (ω) and nucleophilicity index (Nu) were derived from the energy gap (Eg = ELUMO-E HOMO ), are presented in Table 7 . Density functional theory-based quantum chemical descriptors (DFTB-QCDs) help to identify the associations between chemical structures and properties, and play an essential role in chemistry, environment, pharmacology, and health science research [72, 73] . The quantum chemical descriptors such as E HOMO and E LUMO play an important role in predicting the energy gap of molecules and determining the anti-oxidant, Table 6 . Toxicity profile of flavonoids and the two controls predicted using different methods.

ADMETlab vNN hERG H-HT AMES SkinSen DILI Cytotoxicity

The symbols X, p and Ø represents inactive, active and no high confidence prediction available, respectively.

https://doi.org/10.1371/journal.pone.0260853.t006 reactivity, hardness and softness of molecules [73] . The calculated energy gap between HOMO and LUMU for the studied flavonoid derivatives were in the range from 2.428 eV (chrysin (4)) to 7.982eV (genistein (10)), whereas metformin (18) has the largest band gap. The high HOMO energy (-4.734 eV), the small energy gap (2.428 eV), small chemical hardness (0.045 eV), and large softness (11.207 eV) calculated for chrysin (4) molecule makes the molecule a good candidate as a chelating agent. On the other hand, the high dipole moment (8.293 Debyee) value found for chrysin (4) makes the molecule a preferred ligand for biological activity (antimicrobial). Among the 16 flavonoid derivatives, the highest HOMO-LUMO energy gap, hardness, and high nucleophilicity were predicted for genistein (10) ( Table 7) . Based on the quantum chemical reactivity descriptors, the isoflavone class of flavonoids, genistein (10) has the highest chemical hardness (0.147 eV) and the lowest electrophilicity (0.013 eV). Whereas the highest electrophilicity (0.309 eV) and the lowest chemical hardness (0.052 eV), next to chrysin (4) (0.045 eV), was predicted for the chalcone family, butein (3). The electrophilicity, electronegativity and softness quantum chemical descriptors for apigenin (1)/genistein (10) were found to be 0.077/0.147 eV, 6.494/3.409 eV and 0.156/0.013 eV, respectively, indicating good chelating nature of apigenin (1) over genistein (10) .

α-Amylase has been a pivotal target in anti-diabetic drug development. In the present work, we studied the interactions between the Human α-amylase (PDB ID: 4w93) protein against the flavonoids and compared the results against the clinical drug metformin. The 2D and 3D representations of apigenin (1)-, chrysin (4)-and the control metformin (18)-protein interactions are given from Figs 3-5, and the rest of flavonoids are given in S1-S7 Figs in S1 File. The binding affinity, H-bond and residual interaction of flavonoids and metformin (18) (anti-diabetic drug) were summarized in Table 8 

The tested flavonoids showed a binding energy ranging from -7.5 to -8.3 kcal/mol. All the investigated compounds showed strong binding score compared to the clinical drug metformin (18) with -4.5 kcal/mol. The top seven flavonoids with low binding energies were apigenin (1), baicalein (2), chrysin (4), ellagic acid (5), kaempferol (8), quercetin (9) , and luteolin (12) . A previous work reported that compound quercetin (9) has low binding energy (-8.0 kcal/ mol) against α-glucosidase [74] , inferring the inhibitory potential of the compound towards α- The lower binding energies predicted for compounds apigenin (1), chrysin (4), kaempferol (8), quercetin (9), luteolin (12) relative to the control metformin (18) also indicated that they fit well in the binding pocket of the human α-amylase forming a stable inhibitor-protein complex.

In this study we presented the pharmacokinetic, drug-likeness, and toxicity profile of sixteen flavonoid derivatives and two control drugs (metformin (18) and acarbose (17)). The results show appropriate logP, LogS, and high HIA values for chrysin (4), wogonin (15) , baicalein (2), genistein (10) and apigenin (1), suggesting better absorption and less toxicity. This has an implication that the flavonoids can be used as bidentate metal complexing agents, of which chrysin (4) was found to be more ideal for metal complexation with high absorption and lipophilicity with 84.6% absorption compared to metformin (18) (78.3%). The tendency of chrysin (4) to cross BBB and act on CNS was also predicted using SwissADME and ADMETlab predictions. Baicalein (2), butein (3), ellagic acid (5), eriodyctiol (6), fisetin (7) and quercetin (9) were predicted to show carcinogenicity. Most of the molecules are predicted as CYP3A probes, except eriodyctiol (6) which interacts with p-gp. Toxicological endpoints prediction analysis indicated median lethal dose (LD 50 ) values ranging from 159-3919 mg/Kg, of which baicalein (2) and quercetin (9) were found to be mutagenic, whereas butein (3) and morin (15) were found to be immunotoxin among the studied flavonoid derivatives. Based on the findings 

Van-der Walls interactions 1 -8.3 from the predicted drug-likeness, toxicity, pharmacokinetic properties together with quantum chemical descriptors and molecular docking, chrysin (4) is found to be a molecule of potential chelating agent to synthesize coordination compounds for further in vitro and in vivo studies. Moreover, the findings from the DFT calculated quantum chemical reactivity descriptors indicated that chrysin (4) is most likely to interact with metal ions to form stable coordination compounds than the other flavonoids. Overall, the results obtained in this study showed better pharmacokinetics and drug likeness profiles for chrysin (4), wogonin (15) , genistein (10), baicalein (2) and apigenin (1) , and showed toxicity alerts to the toxicological endpoints (above class three toxicity). However, further detailed experimental study is required to confirm the findings reported in this work. 

In silico pharmacology for drug discovery: methods for virtual ligand screening and profiling

Computational approaches in preclinical studies on drug discovery and development

SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules

ADMETlab: a platform for systematic ADMET evaluation based on a comprehensively collected ADMET database

vNN web server for ADMET predictions

Ellagic acid, Kaempferol, and Quercetin from Acacia nilotica: Promising combined drug with multiple mechanisms of action

The influence of lipophilicity in drug discovery and design

Knowledge-based chemoinformatic approaches to drug discovery. Drug discovery today

An analysis of the attrition of drug candidates from four major pharmaceutical companies

VoteDock: consensus docking method for prediction of protein-ligand interactions

In silico and in vitro identification of candidate SIRT1 activators from Indonesian medicinal plants compounds database. Computational biology and chemistry

Design, synthesis and biological evaluation of antimalarial activity of new derivatives of 2, 4, 6-s-triazine

ProTox-II: a webserver for the prediction of toxicity of chemicals

Antidiabetic potential of flavonoids from traditional Chinese medicine: A review. The American journal of Chinese medicine

Synthesis and structure investigation of novel pyrimidine-2, 4, 6-trione derivatives of highly potential biological activity as antidiabetic agent

Estrogenic and antiestrogenic activities of flavonoidphytochemicals through estrogen receptorbinding-dependent and-independent mechanisms

Antidiabetic potential of medicinal plants and their active components

In silico pharmacokinetics and molecular docking studies of lead compounds derived from diospyros mespiliformis

The amylase inhibitor montbretin A reveals a new glycosidase inhibition motif

Quantitative Structure Activity Relationship Study of 2, 4, 6-Trisubstituted-s-triazine Derivatives as Antimalarial Inhibitors of Plasmodium Falciparum Dihydrofolate Reductase. Chemical biology & drug design

Quantitative Structure-Property Relationship (Correlation Analysis) of Phosphonic Acid-Based Chelates in Design of MRI Contrast Agent. Chemical biology & drug design

Theoretical study of molecular structure, pKa, lipophilicity, solubility, absorption, and polar surface area of some hypoglycemic agents

DrugBank 5.0: a major update to the DrugBank database for

ProTox: a web server for the in silico prediction of rodent oral toxicity

Gaussian 09, Revision A.02

Density-functional thermochemistry. III. The role of exact exchange

Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density

Ab initio calculation of vibrational absorption and circular dichroism spectra using density functional force fields. The Journal of physical chemistry

Self-consistent molecular orbital methods. XX. A basis set for correlated wave functions

DFT as a Powerful Predictive Tool in Photoredox Catalysis: Redox Potentials and Mechanistic Analysis

An acylation-Finkelstein approach to radioiodination of bioactives: The role of amide group anchimeric assistance

Theoretical study on chemical fixation of carbon dioxide with aziridine into cyclic carbamate catalysed by purine/HI system

Design and synthesis of three Fe (III) mixedligand complexes: Exploration of their biological and phenoxazinone synthase-like activities

Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings

Rapid calculation of polar molecular surface area and its application to the prediction of transport phenomena. 1. Prediction of intestinal absorption

silico pharmacokinetic and molecular docking studies of natural flavonoids and synthetic indole chalcones against essential proteins of SARS-CoV-2

Methods and principles in medicinal chemistry

Improved In vivo Effect of Chrysin as an Absorption Enhancer Via the Preparation of Ternary Solid Dispersion with Brij® L4 and Aminoclay

What can p K a and NBO charges of the ligands tell us about the water and thermal stability of metal organic frameworks

Characteristics and common properties of inhibitors, inducers, and activators of CYP enzymes. Drug metabolism reviews

A critical analysis of the interplay between cytochrome P450 3A and P-glycoprotein: recent insights from knockout and transgenic mice

Clinically and pharmacologically relevant interactions of antidiabetic drugs. Therapeutic advances in endocrinology and metabolism

Pé rez-Sá nchez H, Apostolides Z. Exploring african medicinal plants for potential anti-diabetic compounds with the DIA-DB inverse virtual screening web server

In silico study of garlic (Allium sativum L.)-derived compounds molecular interactions with α-glucosidase

Alpha-glucosidase inhibitors and hepatotoxicity in type 2 diabetes: a systematic review and meta-analysis

Hepatotoxicity associated with acarbose therapy

Chemistry, quantum mechanics and reductionism: perspectives in theoretical chemistry

Quantum-chemical descriptors in QSAR/QSPR studies

Molecular docking studies for discovery of plant-derived a-glucosidase inhibitors

Adama Science and Technology University, Ethiopia and Debre Markos University (Ethiopia) are duly acknowledged for providing research and administrative facilities, respectively, provided to Mr. Mamaru Bitew. TBD acknowledges Department of Chemistry, University of Botswana for the research facilities. Computer time was provided by CERIT-SC (CERIT Scientific Cloud) of the Czech Republic, and greatly acknowledged.

Conceptualization: Mamaru Bitew, Rajalakshmanan Eswaramoorthy.