Investigating the binding affinity, molecular dynamics, and ADMET properties of 2,3-dihydrobenzofuran derivatives as an inhibitor of fungi, bacteria, and virus protein

2,3-Dihydrobenzofurans (DHB) have proposed as advantages structures, and used as chemical entresol to design small compound libraries. The present study illustrates to explore 2,3-dihydrobenzofurans(DHB) in comparison to selected some derivatives drugs by using molecular docking and molecular dynamics, as well as ADMET studies. The online database “Molinspiration online server” was used to detect the physicochemical pharmacokinetics and drug likeness score of DHB drugs. For estimation of molecular docking, six pathogens, such as Aspergillus niger (PDB id: 1kum), Candida albicans (3dra), Escherichia coli (6og7), Salmonella typhi (4k6l), Influenza (1ru7), and Hepatitis C (4tyd), were chosen due to close biological studies. From Molinspiration online server has showed that DHB did not violate the “Lipinski five rule” as drugs, leading compound for molecular docking exhibited the potential interaction to the active residue. The binding affinity of DHB2 (−7.00 kcal/mol) against 3dra was higher than DHB8 (−6.40 kcal/mol) and DHB (5.70 kcal/mol) for compounds. The results of molecular docking show that the compounds mentioned in this study are not equally effective against pathogens, such as fungi, viruses, and bacteria. However, DHB2, DHB3, and DHB 8 compounds can work against almost given pathogens which results are derived from auto dock vina in terms of binding affinity around 6.00 kcal/mol, and Fire Dock has values from about 38.0 to 42.0 kcal/mol. To explore the dynamic nature of the interaction, 50 ns molecular dynamics (MD) simulation was performed on the selected protein-DHB complexes. Thus, DHB 8 has greater potential to interact for further for fungi. Finding from this study can play an effective role as a drug in any biological system. This study as well recommends to researchers to synthesize these DHBs for evaluation of its biological activity.

numbers [5] and reaction mechanism [6]. Due to vast applications, molecular modeling has become an essential part of concurrent drug discovery processes though new molecules [7,8]. In current yeast, the charge of these methods has raised significantly and it takes over a decennia for too much short shred of compounds to pass the drug discovery pipeline from basic screening shocks or leads, chemical optimization, and clinical trials before launching into the market as a drug [9]. The pathway and methodologies used in drug design have changed over time, parasitic, and leading new technological advances to evolve the several bottlenecks found along the path [10].
There is various scope used for docking, including SUR-FLEX, FlexX, FRED, GLIDE, GOLD and DOCK-6, and has been evaluating and these programs verified to produce unfailing poses in numerous docking studies on drug molecules [11]. To 1990, the principal issues were operating to discovery and chemical synthesis of drug-like molecules [12], the emersion of combinatorial chemistry [13], gene technology [14], and high-throughput tests which has shifted the focus on various computational tools due to save the consumption of time, experimental cost and labor [15][16][17][18]. On the other hand, poor absorption, distribution, metabolism, and excretion (ADME) properties of new drugs have captured more attention [19], besides the lipophilicity plays an important role on drug discovery [20]. To predict the binding of protein with strong interaction partners as ligands or drugs, protein docking has performed by computational tools and selected the bind site of protein with binding affinity score. Rarey et al. reported on docking and its applications thought the all procedure that it is closely saying about how a drug molecule relative orientation to protein, interaction with protein, and position after interaction [21,22]. Besides, molecular docking is the best methods in structure-based drug design, and it can foretell the binding-conformation of very small molecule ligands to the appropriate target binding site of protein [23]. Two characterization of the binding behavior plays an important role in rational design of drugs which are protein and the ligand as complementary surfaces [24,25] and ligand-protein pairwise interaction energies [26].
2,3-Dihydrobenzofurans has been recognized as key pharmacophores in the field of many natural and synthetic bioactive molecules [27]. The most important place of dihydrobenzofuran and their derivatives are used as biological significance molecules, such as anti-cancer [28], anti-HIV [29], anti-inflammatory or antibiotic [30], antitubercular evaluations [31], leishmanicidal agents [32], and induce immunogenic cell death (ICD) of specific cancer cells [33]. Moreover, they have been found to induce the secretion of an endogenous anticancer factor, antileishmania drugs, and immunomodulatory activities [34].
Similarly, natural products and many organic compounds were obtained from synthetic methodologies from 2,3dihydrobenzofuran moieties [35]. Due to enormous applications and biological significances, we have designed this study as new leading compounds as drug against some fungi, bacteria, and virus proteins. Following a successful strategy of repeated steps of molecular docking calculations, computational study, binding assessment, and ADMET studies have employed on this work.
1ru7 Intracellular transport of viral protein in host cell, receptormediated endocytosis of virus by host cell, viral genome packaging, viral protein processing [41].
4tyd Transformation of host cell by virus, serine-type peptidase activity, ATP binding, metal ion binding [42].
Then, the crystal structure of the protease was optimized and checked by Swiss-PDB viewer software packages (version 4.1.0) [43] based on their least energy. Some significant factors, such as improper bond order, side chain geometry, and missing hydrogen, were observed in the crystal structure of the protease. PyMol (version 1.1) [44] software package was used to erase all the hetero atoms, water molecules, and inhibitor present. Molecular docking by using AutodockVinawas done for molecular docking, and insights of these results Patch-Dock web server were used for ligand-receptor docking [45]. 1000 solutions with area, six-dimensional transformation space and score were gained from Patch-Dock server, and then all solutions were afterward shifted into FireDock to retouch the ten best solutions associated with global energy. The adopted complexes from Fire-Dock were aligned according to minimum global binding energy. In the end, the Discovery Studio 4.1 Client was employed for the visualization of binding modes of the receptor and ligands [46]. The grid box in Autodock-Vina was generated targeting the active site of the protein, where the center was at X = 34.92, Y = 87.74, Z = 80.16, and the grid box was maximum different dimensions for six proteins (unit of the dimensions, Å). Finally, we docked total seventy-two of twelve ligand (DHB) against six proteins. ClusPro online Docking was performed on the best five selected docking result obtained from FireDock and AutoDockVina study.

Molecular dynamic (MD) simulations of selected DHB
MD simulation was devoted to underpin the docking results gained for the best antiviral drugs. 100 ns MD simulation was performed for holo-form (drug-protein).
YASARA dynamics [48] program was used for all simulation applying AMBER14 force field [49]. In the presence of a water solvent, the total system was equilibrated with 0.9% NaCl at 298 K temperature. By simulation, the atomy Mesh Ewald algorithm was discussed for radial electrostatic interactions. A cubic cell was propagate within 20 Å on every side of process and periodic boundary circumstance was proposed during the simulation. A time step of 1.25 fs was maintained to carry out 10 ns MD simulation and the snapshots were taken at every 100 ps.

Drug-like properties
DHB membrane permeability, molecular properties, for instance, and bioavailability of guidance compounds relay on some basic properties of molecules which are molecular weight (MW), partition coefficient (logP), and number of hydrogen bond acceptors/donors, and those properties are developed by Lipinski et al. known as Lipinski "rule of five" shoes (Table S2) [50]. As a matter of principle, an orally active drug should have no more than one interruption of the following status: (1) no more than five hydrogen bond donors, (2) no more than ten hydrogen bond acceptor, (3) molecular mass of less than 500 daltons (Da), and (4) an octanol-water partition coefficient of not more than five. As all the numbers of these recommendations are multiple of 5, the rules are termed Lipinski's rule of five [50]. Table 1 represents all DHB compounds which followed the rule of five indicating the good bioavailability score of DHB. The druglikeness score of command molecules is ascertained with a multitude of GPCR, ion channel modulator, a kinase inhibitor, nuclear receptor ligands, protease inhibitor, and enzyme inhibitor which has been applied to look into the efficiency of molecules to tether for drug development. Srivastava et al. [51] expatiated that the richer the bioactivity score has a higher likelihood of the specific molecule to be active as a drug. It is well documented that, if bioactivity score is more than 0, the molecules have better biological activity. If the bioactivity score is −0.5 to 0, the molecules have moderate activity and less than −0.5, the molecules have no biological activity [52]. The results of bioactivity data indicated that DHB 1 and DHB 2 were highly active as enzyme inhibitor, nuclear receptor ligand, ion channel modulator, GPCR ligands, and protease inhibitor. Other compound (DHB 1, 6-10, 12) were predicted as moderately active as GPCR ligands and all.

ADMET properties
Absorption, Distribution, Metabolism, Excretion and Toxicity (ADMET) properties of drug candidates or environmental chemicals play a key role in drug discovery and environmental hazard assessment. The admetSAR [47,53] was developed as a comprehensive source and free tool for the prediction of chemical ADMET properties [54]. Total eight ADMET parameters including human intestinal absorption (HIA), blood-brain barrier (BBB), human oral bioavailability, water solubility, acute oral toxicity, and Tetrahymenapyriformisp, IGC50, were tested for the twelve DHB by AdmetSAR prediction [55]. The results are summarized in Table 2. The results provide the information for molecules that the DHB are safer to use, because twelve DHB are perceived to be non cancerogenic, and have good passive gastrointestinal absorption and blood-brain barrier.
The rationale of the blood-brain barrier has to reveal that it protects against circulating toxins or pathogens which are responsible for brain infections. Sometimes, it is accountable for allowing vital nutrients to reach the brain. From Table 2, it is observed that all DHBs show positive activity to their functions. The human oral bioavailability (HOB) represent the therapeutic variability or dose fraction of new chemical. Moreover, HOB indicates the positive when the logK(%F) is greater than zero [56]. Thus, DHB 4, 9, and 11 show the positive value and others are negative value. On the other hand, fish aquatic toxicity prediction illustrates that these compounds may harm aquatic ecology. Besides, DHB may also damage the acute ecology which is noticed from the data of acute oral toxicity, and it is found from 1.514 to 3.041 kg/mol with varying compounds. Finally, all DHBs are highly water soluble which indicates that it has high affinity to water, and easily dissociate with water when it will be taken in human body as drug.

Molecular docking
Molecular docking are usually expressed by the binding affinity of a drug with the protein of pathogen, besides it also mentions about the number of hydrogen bonds, hydrophobic bonds, polar and non-polar bond, and van dar Waal force interaction between drug and protein.
There were some studies for development of anticancer drugs [57], antibiotics, and drug discovery by small molecules [58]. This study consists of two type molecular docking, such as auto dock score and FireDock score. The standard value of auto dock vina is 6.0 or more for  [59] that docking score has been considered as the Gibbs free energy of binding, and the Gibbs free energy of binding is more −6.0 kcal mol −1 indicates the active drug and inactive (Gibbs free energy of binding) is less −6.0 kcal mol −1 . Molecular docking by using Patch-Dock web server and refine 1000 solutions with Fire-Dock web server molecular docking is performed to find out the best candidates among the 12 DHB based on their minimum global binding energy. FireDock were ranked according to minimum global binding energy. Global energy of the DHB is distributed within the range from −28.68 to −47.53 kcal/mol shown in Table 3 and Fig. 1. Selected DHB are screened primarily using Patch-Dock and FireDock web server to find out the best candidates, and then, the AutoDockVina was employed to understand their fitness. Six best molecular docking shown in Fig. 2. The binding affinity of AutoDockVina and FireDock score of all DHB are shown in Table 3.

Molecular dynamics simulation
MD simulation for each complex with three selected DHB derivatives like DHB 2, DHB 3, and DHB 8 has performed against two proteins (1kum and 6og7) for 10 ns and one protein complex (3dra) for 50ns, because the DHB 2, DHB 3, and DHB 8 derivatives are found as the highest binding affinity against these pathogens. The MD simulation of DHB3 against 6og7 protein demonstrates the lowest RMSD value than the other complexes which indicates its greater stability, and keeps in the pocket of proteins. As a result, DHB3 has to be considered as drug. Moreover, the lowest RMSF value of the DHB3 against 6og7 has also found after MD simulation. Secondly, the DHB2 against protein (3dra)  illustrates binding affinity (−7.00 kcal/mol); in this case, it was selected for MD simulation for fitting in the protein pocket. In general, the DHB2 against 1kum, DHB3 against 1kum, DHB8 against 1kum, DHB2 against 6og7, DHB3 against 6og7, and DHB8 against 6og7 schemes are chosen for MD simulation on based of binding affinity from molecular docking studies. Analysis of RMSD, RMSF and Backbone of selected three complexes (DHB2 against 1kum, DHB3 against 1kum, DHB8 against 1kum, DHB2 against 6og7, DHB3 against 6og7, DHB8 against 6og7) with protein at 10 ns MD simulations are presented in Fig. 3.

Discussions
If the oral bioactivity score (%F) of the molecule is greater than 0.00, has the important biological activities and scores between 0.50 and 0.00, will be considered to be fairly active, and if the value is less than 0.50, it is guessed to be inactive. The gained values of druglikeness score stated that DHB 2 followed the good drug-likeness along with other standard drugs, for instance DHB 3 and DHB 8 (Table 1). Molecular docking is a computational tool that virtually seeks to predict a complex of two binding partners, such as drugs and protein or macromolecules. Withal, the molecular docking is equivalent the interaction of ligand and protein as a specific scoring function through the binding free energy which is the sum of intermolecular interaction between them after docking. The binding energy is also expressed by Gibbs free energy of binding from docking. Regarding a standard drugs, the binding affinity is equal to 6.00 kcal mol −1 or more [59,60]. In general, the AutoDock approves a simultaneous sample process to arrangement with side chain flexibility, active site of protein and side chains of the receptor. Thus AutoDock evaluates the binding affinity of ligand and protein interaction which is the sum of the binding constant (K d ) and the Gibbs free energy (ΔG L ) [61], and 6.00 kcal mol −1 is the starting value for good drugs small structural molecules for theoretical value [62,63].
Analysis of RMSD, RMSF, and Backbone of DHB2 against 3dra complex is also listed in Fig. 4 and 5ns snapshot showing in Fig. 5. Finally, it has been found that the stability and protein interact with drug molecules are almost right to the height range of MD study, and it has provided information as a drug.
Analysis of the nonbonding interactions of the best two DHB (2, 3 & 8) with the six proteins reveals that the selected compounds interact with either both or at least one catalytic residue detected by FireDock shown in Table 4 and Fig. 1. For molecular docking for DHB2+4tyd, the eight hydrogen bonds and eight hydrophobic bonds are showed in Table 4. The result of DHB3+6og7 docking is illustrated in Table 4 which are listed by five hydrogen bonds and eight hydrophobic bonds. In case of DHB2+1kum docking, there are found the three hydrogen bonds and seven hydrophobic bonds showing in Table 4. We observed maximum four interaction of TRP563 and three interaction of THR525 in DHB2+1kum. Maximum pi-alkyl interaction is observed in DHB2 + 4tyd. From these studies, the best binding affinity vote is going to DHB2 against 3dra shown in Fig. 6, because it contains the maximum protein residue interactions with drugs or DHB2, as well as active sites. In addition, from Fig.  5, it is illustrated the evaluation of the distribution of drugs within the biological system through the MD simulation that DHB2 remains in protein pocket with wide range of dynamics which is other supporter information for the best drugs (Fig 1).

Conclusions
Summarizing from the above discussion, it is also found that the compounds used in the given work show great interest in biological activity. The first thing is that came  to light when designing the molecules was the replacement of the benzofuran and its derivatives using changing the site chain by the alkyl chain. First of all, the SAR was studied where it was shown that with the change of alkyl the biological properties will be the lightness, talking molecular docking studies and ADME T. Secondly, the twelve compounds used in this chapter are found to have more or less biological activity, of which DHB2, DHB3, and DHB8 may exhibit the most biological activity. The reason molecular docking has been used as the basis for determining their properties is because it makes it easier to understand how much value a drug can attach to a particular protein of pathogens. Since the value of talking is considered to be close to 6 in any of the previous standard cases, the value of DHB2, DHB3 and DHB8 is close to or above 6. On the other hand, to be used as a molecular drug, it is first confirmed by molecular screening, then it goes to the human body and a molecular dynamic study is done to find out what dynamics might be. Molecular dynamics shows that DHB2, DHB3, and DHB8 drugs are located in the pocket of the protein. As a result, it can play an effective role as a drug in any biological system. In addition, AMDET studies have been performed for their side effects, suggesting that not all drugs have anti-cancer activity and that they exhibit aquatic and fish toxicity.