Computational and Preclinical Prediction of the Antimicrobial Properties of an Agent Isolated from Monodora myristica: A Novel DNA Gyrase Inhibitor

The African nutmeg (Monodora myristica) is a medically useful plant. We, herein, aimed to critically examine whether bioactive compounds identified in the extracted oil of Monodora myristica could act as antimicrobial agents. To this end, we employed the Schrödinger platform as the computational tool to screen bioactive compounds identified in the oil of Monodora myristica. Our lead compound displayed the highest potency when compared with levofloxacin based on its binding affinity. The hit molecule was further subjected to an Absorption, Distribution, Metabolism, Excretion (ADME) prediction, and a Molecular Dynamics (MD) simulation was carried out on molecules with PubChem IDs 529885 and 175002 and on three standards (levofloxacin, cephalexin, and novobiocin). The MD analysis results demonstrated that two molecules are highly compact when compared to the native protein; thereby, this suggests that they could affect the protein on a structural and a functional level. The employed computational approach demonstrates that conformational changes occur in DNA gyrase after the binding of inhibitors; thereby, this resulted in structural and functional changes. These findings expand our knowledge on the inhibition of bacterial DNA gyrase and could pave the way for the discovery of new drugs for the treatment of multi-resistant bacterial infections.


Introduction
Bacterial infections cause more deaths in the developing nations than the developed nations [1]. A report commissioned by the UK government anticipates an increase in deaths caused by antimicrobial-resistance by 2050, with a prediction that 10 million people will fall victim to antibiotic-resistant infections annually: roughly 4.73 million in Asia and 4.15 million in Africa, in contrast to 0.39 and 0.32 million in Europe and the US, respectively [2]. Several countries identified as lower-middle income countries are more likely to lack access to clean water, be faced with sanitation and hygiene problems, and have greater numbers of their people be affected by bacterial pathogens. Moreover, the malnourished, immunocompromised, and HIV-positive citizens of these countries are particularly affected [2].
DNA gyrase, a member of the family of type IIA topoisomerases, is an essential protein in bacteria, as it introduces negative supercoils in their DNA via an ATP-dependent mechanism. The enzyme is an A2B2 tetramer consisting of GyrA and GyrB, and its inhibition has catastrophic effects on bacterial cell proliferation and survival; thus can serve as a major target for the development of antibacterial agents [3][4][5][6]. The DNA gyrase complex requires the transition between several conformational states to cleave duplex DNA, pass the DNA strand through the break, and relegate the DNA fragments [4][5][6][7].
Recently, several antibiotics such as levofloxacin have gained prominent importance in the market. However, bacterial resistance and their adverse effects (including dizziness, blunt vision, and headache) represent a serious concern. The re-evaluation of DNA gyrase as a target for the development of new antibiotics received considerable applause from researchers aiming to find a lasting solution to this menace. DNA gyrase inhibitors could act as effective agents against a broad range of pathogens, bind to new sites in these proteins, and inhibit their enzymes via a mechanism that is different from that of quinolones [8][9][10][11].
In recent years, a tremendous amount of research was conducted in order to explore the pharmacological utility of medicinal or plant-derived natural products [12]. Over the years, plants have been used for medicinal purposes by diverse people and cultures throughout the world, and many of them have been effective in treatment of various illness and diseases. The use of plants for medicinal purposes continues to this day, mostly in the form of traditional medicine, which is now recognized by the World Health Organization as a building block for primary health care [13][14][15].
African nutmeg (Monodora myristica), belongs to the family of Ananacea, and it is an important species of the evergreen forest of West Africa (mostly common in the Southern part of Nigeria) [16,17]. The seeds that are embedded in the white sweet-smelling pulp of the sub-spherical fruit of Monodora myristica are of particular economic importance, whereas the kernel when ground to powder is a popular condiment used for the preparation of pepper soup that serves as a stimulant for the relief from constipation and is also used to control the passive uterine hemorrhage that occurs immediately after birth [17,18].
Many bioactive compounds were identified in the oil extract of the plant as reported by Adewole et al. [12], including 1,3,3-trimethyl-2-oxabicyclo[2.2.2]octan-6-ol, called eucalyptol; it is a natural organic compound and a colorless liquid that is used as an ingredient in many brands of mouthwash and cough suppressants. Eucalyptol is known to control the hypersecretion of airway mucus and asthma, as it can suppress arachidonic acid metabolism and cytokine production in human monocytes [19,20]. It was also found to possess antibiotic activities [21].
The therapeutic potentials of fatty acids (FAs) were found to be particularly promising for the treatment of various microbial infections. The antibacterial activity of FAs is wellknown in the literature and represents a promising option for the development of the next generation of antibacterial agents [22]. Moreover, FAs are highly involved in living organisms' defense systems against numerous pathogens, including multidrug-resistant bacteria. Previous studies established a computational approach to the investigation of protein-ligand interactions and their stability in relation to the structure and the function of a receptor [22]. We previously identified hit molecules from plant extracts that are theoretically capable of protein inhibition [23][24][25][26][27]. Herein, we employed computational prediction tools in order to study the binding energy and interaction of these ligands against DNA gyrase. Finally, a Molecular Dynamics (MD) simulation was undertaken in order to investigate the stability of the interaction over a period of 100 ns. Therefore, this study utilized the characterized secondary metabolites deriving from Monodora myristicaextracted oil as potential antimicrobial agents that could suppress the bacterial cell growth, and it employed a computational approach to predict alternative compounds with drug-like properties and better inhibitory activities.

Results
The estimated number of hydrogen bonds that would be donated by the solute to water molecules in an aqueous solution is proportional to the number of hydrogen bonds that will be accepted by the solute when they interact with water molecules. Chemical structures of eucalyptol when compared with standard ligands are shown in Figure 1. The eucalyptol showed a successfully processed molecule on predicted central nervous system activity of +1 as shown Table 1. the next generation of antibacterial agents [22]. Moreover, FAs are highly involved ing organisms' defense systems against numerous pathogens, including multidr sistant bacteria. Previous studies established a computational approach to the inv tion of protein-ligand interactions and their stability in relation to the structure an function of a receptor [22]. We previously identified hit molecules from plant extrac are theoretically capable of protein inhibition [23][24][25][26][27]. Herein, we employed com tional prediction tools in order to study the binding energy and interaction of the ands against DNA gyrase. Finally, a Molecular Dynamics (MD) simulation was u taken in order to investigate the stability of the interaction over a period of 100 ns. fore, this study utilized the characterized secondary metabolites deriving from Mo myristica-extracted oil as potential antimicrobial agents that could suppress the ba cell growth, and it employed a computational approach to predict alternative comp with drug-like properties and better inhibitory activities.

Results
The estimated number of hydrogen bonds that would be donated by the so water molecules in an aqueous solution is proportional to the number of hydrogen that will be accepted by the solute when they interact with water molecules. Ch structures of eucalyptol when compared with standard ligands are shown in Figure  eucalyptol showed a successfully processed molecule on predicted central nervou tem activity of +1 as shown Table 1.   Therefore, with the crucial analysis carried out on the efficacy of the molecules, all of the examined parameters were within Lipinski's rule of five (ROF) cut-off range for the test compound 1,3,3-trimethyl-2-oxabicyclo[2.2.2]octan-6-ol and showed no violation of Lipinski's ROF and Veber rules as shown in Table 2. Furthermore, Table 3 shows that the ligand of interest (529885) exerted no inhibition upon the Pgp substrate when compared with cephalexin and levofloxacin, which were predicted to be inhibitors of the Pgp substrate. This showed that the ligand revealed an ATPdependent drug efflux pump for xenobiotic compounds with broad substrate specificity. Table 3. Comparative pharmacokinetics and toxicological analyses of studied molecules. Pharmacokinetics analyses were performed with SwissADME, and the toxicological analyses were performed with Pred-hERG and AdmetSAR servers. produced excellent results when compared with those of cephalexin and levofloxacin (Table 4).

PubChem GI Abs
The docking analysis approach already revealed several suitable drug molecules for the target of interest [28][29][30]. This was also the case after undertaking a similar approach with the DNA gyrase subunit A N-terminal domain and the E. coli DNA gyrase-DNA binding and cleavage domain in State 1 receptor complex with 1,3,3-trimethyl-2oxabicyclo[2.2.2]octan-6-ol, producing a competitive docking score (−4.115 kcal/mol). At the same time, the docking scores for cephalexin and levofloxacin were −6.13 kcal/mol and −4.658 kcal/mol, respectively (Table 5). This same inhibitory effect was demonstrated with the highest level of glide ligand efficiency, and it is indicative of a more suitable binding target of interest against the control ligands. Furthermore, based on specific interactions, the ARD98 in the DNA gyrase subunit A N-terminal domain receptor was involved in hydrogen bonding with the atoms showing no π-π interaction or salt bridge formed with 1,3,3-trimethyl-2-oxabicyclo[2.2.2]octan-6-ol, whereas GLY120, ARG98, SER118, and THR230 were involved in hydrogen bonding in the DNA gyrase subunit A N-terminal domain receptor with atoms in cephalexin and levofloxacin. Cephalexin had ARG98 at the salt bridge on interaction, whereas levofloxacin had no atoms at the salt bridge interaction. Moreover, GLU381 at the salt bridge interacted with levofloxacin but had no interaction at the salt bridge of cephalexin. On the other hand, VAL420 and GLU381 were involved in E. coli DNA gyrase-DNA binding and cleavage domain in State 1 receptor hydrogen bonding with the atoms in cephalexin and levofloxacin, whereas ALA421 and VAL420 were involved in the hydrogen bonding with atoms with no π-π interaction or salt bridge formed in 1,3,3-trimethyl-2-oxabicyclo[2.2.2]octan-6-ol, as  Table 5. This is an indication of good measurable binding affinities for receptor residues with ligand contribution to the flexibility of the target.
The molecular analysis undertaken on the 3D and 2D structures of the complexes formed by the top three compounds (i.e., 1,3,3-trimethyl-2-oxabicyclo[2.2.2]octan-6-ol acetate, cephalexin, and levofloxacin) with the target proteins revealed that the ligand of interest can occupy the active site of the enzyme (Figures 2-4).
in the DNA gyrase subunit A N-terminal domain receptor with atoms in cephalexin and levofloxacin. Cephalexin had ARG98 at the salt bridge on interaction, whereas levofloxacin had no atoms at the salt bridge interaction. Moreover, GLU381 at the salt bridge interacted with levofloxacin but had no interaction at the salt bridge of cephalexin. On the other hand, VAL420 and GLU381 were involved in E. coli DNA gyrase-DNA binding and cleavage domain in State 1 receptor hydrogen bonding with the atoms in cephalexin and levofloxacin, whereas ALA421 and VAL420 were involved in the hydrogen bonding with atoms with no π-π interaction or salt bridge formed in 1,3,3-trimethyl-2-oxabicyclo[2.2.2]octan-6-ol, as shown in Table 5. This is an indication of good measurable binding affinities for receptor residues with ligand contribution to the flexibility of the target.
The molecular analysis undertaken on the 3D and 2D structures of the complexes formed by the top three compounds (i.e., 1,3,3-trimethyl-2-oxabicyclo[2.2.2]octan-6-ol acetate, cephalexin, and levofloxacin) with the target proteins revealed that the ligand of interest can occupy the active site of the enzyme (Figures 2-4).  The three compounds, like the standard ligands, interacted with the DNA gyrase subunit A N-terminal domain at the active site amino acid residues SER118, THR230, and ARG98, respectively, and they also interacted with the E. coli DNA gyrase-DNA binding and cleavage domain in State 1 at the active site amino acid residues ALA421, VAL420, and HIS471, respectively ( Figures 5-7).
The dynamic behavior of the protein residues was further studied by examining the RMSF pattern and by calculating the mean fluctuations of the Apo protein in complex with the standard drugs and the studied molecules (Table 6). In agreement with the RMSD results, the fluctuation patterns of the complex-forming molecules were relatively the same as that of the Apo protein. The highest mean fluctuation was observed when the Apo protein was complexed with levofloxacin (1.57 Å), followed by 1KIJ-529885 (1.35 Å). The Apo protein had the smallest mean fluctuation of 1.      The three compounds, like the standard ligands, interacted with the DNA gyrase subunit A N-terminal domain at the active site amino acid residues SER118, THR230, and ARG98, respectively, and they also interacted with the E. coli DNA gyrase-DNA binding and cleavage domain in State 1 at the active site amino acid residues ALA421, VAL420, and HIS471, respectively (Figures 5-7).   The three compounds, like the standard ligands, interacted with the DNA gyrase subunit A N-terminal domain at the active site amino acid residues SER118, THR230, and ARG98, respectively, and they also interacted with the E. coli DNA gyrase-DNA binding and cleavage domain in State 1 at the active site amino acid residues ALA421, VAL420, and HIS471, respectively (Figures 5-7).    The dynamic behavior of the protein residues was further studied by examining the RMSF pattern and by calculating the mean fluctuations of the Apo protein in complex with the standard drugs and the studied molecules (Table 6). In agreement with the RMSD results, the fluctuation patterns of the complex-forming molecules were relatively the same as that of the Apo protein.

Discussion
Despite the important role of bacteria in the environment in which we live, bacterial infections exert a large impact on public health across the globe [31]. Although bacterial infections are considered to be easier to treat than viral infections, the resistance of the bacteria to several antimicrobial agents is becoming a growing concern with potentially devastating consequences [31]. To date, attention is mainly directed toward the development of small molecules for the targeting of bacterial infections; however, the current clinical pipeline is insufficient in tackling the sporadic emergence and spread of antimicrobial resistance [32]. Therefore, we herein employed computational screening methods in order to assess the therapeutic potential of secondary metabolites deriving from Monodora myristica-extracted oil as antimicrobial agents [12]. The undertaken computational analysis revealed that 1,3,3-trimethyl-2-oxabicyclo[2.2.2]octan-6-ol can display a better binding affinity when compared with approved antimicrobial agents, and it is characterized by a competitive docking score and glide ligand efficiency against levofloxacin (as shown in Table 6). The inhibitory potential of the protein-ligand complex was revealed by the docking score (Table 6) [33]. It is known that ligands are predicted to have a hydrogen bond donor score ≤5 and a hydrogen bond acceptor score ≤10 with successful processed molecular structures on a predicted central nervous system activity of +2. Therefore, the ligand of interest showed a successfully processed molecule on a predicted central nervous system activity of +1. An attempt to change any of the single molecular properties of the aforementioned could result in absorption and bioavailability [34]. However, the main test ligand (1,3,3-trimethyl-2-oxabicyclo[2.2.2]octan-6-ol) exhibited the most desirable pharmacological potential, with the parameters fitting within the ROF cut-off range and not violating the rule.
Furthermore, the ligand 1,3,3-trimethyl-2-oxabicyclo[2.2.2]octan-6-ol exhibited a good metabolic profile with the Pgp substrate, and thereby, it showed no inhibitory effects on the enzyme's activity when compared with those of standard antimicrobial agents cephalexin and levofloxacin (Table 2). Moreover, the blood-brain barrier's (BBB) permeant inhibitory property of the target ligand 1,3,3-trimethyl-2-oxabicyclo[2.2.2]octan-6-ol was increased, whereas it was not inhibited in the case of cephalexin. The potentiality of any drug candidate is traceable to oral bioavailability with a uniqueness in its functions, and it would enable the drug to pass through a cellular membrane. Otherwise, the drug could easily be trapped within these barriers, thereby posing a serious health challenge. Among these functions are desolvation, diffusion, resolvation, and physicochemical properties such as lipophilicity [35]. The context of the absorption, distribution, and excretion of drugs provides several ways to identify pharmacokinetic models, such as human intestinal absorption, aqueous solubility (absorption), p-glycoprotein inhibition (distribution), and renal organic cationic transporters inhibition (excretion). Therefore, it is important to identify the aforementioned pharmacokinetic parameters, as they provide additional knowledge for drug acceptability or rejection in any drug development process [36]. The hit molecule exhibited its permeability effect in the human intestinal membrane and the BBB.
Furthermore, in an attempt to validate the pharmacological potential of the hit compound, the native ligand 1,3,3-trimethyl-2-oxabicyclo[2.2.2]octan-6-ol was docked with the drug targets, and the intermolecular interactions with most of the crucial amino acid residues were revealed (Table 6); thereby, it suggests an interaction of the amino acids with the ligand at the binding site, which could infer that the selected hit served as an antimicrobial agent.
Moreover, MD simulation studies were performed for 100 ns to analyze the atomic level changes in DNA gyrase with respect to the timescale. In order to gain insight into the effect of the selected compounds as efficient inhibitors of DNA gyrase, their stability in the protein's active site was investigated by using Schrödinger version 2022-1 to conduct the MD simulation on the native protein and in complex with the compounds. MD properties including the root mean square deviation (RMSD), the root mean square fluctuation (RMSF), the radius of gyration (rGyr), and the solvent-accessibility surface area (SASA) were analyzed.
Furthermore, the Analysis of the best-predicted docking pose interaction as shown in (Figure 8) and (Figure 9) revealed E. coli DNA gyrase-DNA binding domain with 1,3,3trimethyl-2-oxabicyclo[2.2.2]octan-6-ol complexed with 1KIJ-cephalexin, as superimposed with the standard. Also, the RMSD and the RMSF were plotted in order to examine the stability ( Figure 10A) and the flexibility ( Figure 10B) of the Apo protein compared with its interaction with inhibitory molecules. The rGyr and the SASA were plotted in order to analyze the compactness ( Figure 10C) and to compute the secondary structure of the complexes ( Figure 10D). in (Figure 8) and (Figure 9) revealed E. coli DNA gyrase-DNA binding domain with 1,3,3trimethyl-2-oxabicyclo[2.2.2]octan-6-ol complexed with 1KIJ-cephalexin, as superimposed with the standard. Also, the RMSD and the RMSF were plotted in order to examine the stability ( Figure 10A) and the flexibility ( Figure 10B) of the Apo protein compared with its interaction with inhibitory molecules. The rGyr and the SASA were plotted in order to analyze the compactness ( Figure 10C) and to compute the secondary structure of the complexes ( Figure 10D).   The structural stability of the protein alone and in complex with the molecules w interpreted by the deviation of the plot and the mean ± SEM (Table 6. The deviation p terns of the each of these complexes and of the Apo protein were relatively the sam throughout the simulation time of 100 ns. The RMSD alteration was observed in all of t RMSD patterns at 30-60 ns, after which they all attained an equilibrium until the end the simulation period. By comparing the mean ± SEM of the plots, the 1KIJ-529885 co plex exhibited the lowest RMSD value (2.87 Å ), whereas the 1KIJ-27447 complex (3.89 appeared to be higher when compared to the Apo protein (3.16 Å ) and the 1KIJ-levoflo acin complex (3.38 Å ). The differences in mean deviation explain the impact of the int action of the molecules in the binding pocket of the receptor, and thus, provide a basis further investigation.
Subsequently, the rGyr for the complexes were analyzed in order to deduce th changes in compactness over the 100-ns simulation period ( Figure 10C and Table 5). T plot shows that the 1KIJ-529885 complex had the smallest rGyr mean value of 2.13 Å , a it was stable from the onset and throughout the 100-ns simulation period. The 1KIJ-274 complex exhibited a mean rGyr value of 3.85 Å , and its pattern fluctuation was observ The structural stability of the protein alone and in complex with the molecules was interpreted by the deviation of the plot and the mean ± SEM ( Table 6). The deviation patterns of the each of these complexes and of the Apo protein were relatively the same throughout the simulation time of 100 ns. The RMSD alteration was observed in all of the RMSD patterns at 30-60 ns, after which they all attained an equilibrium until the end of the simulation period. By comparing the mean ± SEM of the plots, the 1KIJ-529885 complex exhibited the lowest RMSD value (2.87 Å), whereas the 1KIJ-27447 complex (3.89 Å) appeared to be higher when compared to the Apo protein (3.16 Å) and the 1KIJlevofloxacin complex (3.38 Å). The differences in mean deviation explain the impact of the interaction of the molecules in the binding pocket of the receptor, and thus, provide a basis for further investigation.
Subsequently, the rGyr for the complexes were analyzed in order to deduce their changes in compactness over the 100-ns simulation period ( Figure 10C and Table 5). This plot shows that the 1KIJ-529885 complex had the smallest rGyr mean value of 2.13 Å, and it was stable from the onset and throughout the 100-ns simulation period. The 1KIJ-27447 complex exhibited a mean rGyr value of 3.85 Å, and its pattern fluctuation was observed at the beginning of the simulation time until 30 ns, after which there were no observable changes. The 1KIJ-levofloxacin complex had the highest mean rGyr of 3.97 Å; however, a noticeable drop in its pattern was observed at 60 ns, but the pattern was stable until 90 ns, after which it returned to its original position until the end of the simulation. Taken together, the two molecules displayed a slight pattern of deviation and fluctuation when compared to the Apo protein and the standard complex. Furthermore, these complexes showed even greater compactness than that exhibited by the standard (levofloxacin). All of these results demonstrate that the studied molecules formed highly stable complexes with DNA gyrase.
In addition, the SASA examined the interaction between the protein surface and surrounding solvent molecules. The interactions during the protein folding may be crucial to its stability and rearrangement with respect to its structure. As such, the SASA values of the complexes with respect to the standard were computed ( Figure 10D and Table 5). The standard complex had a mean SASA value of 72.78 Å, whereas both the 1KIJ-27447 and the 1KIJ-529885 had higher mean SASA values of 127.30 Å and 150.99 Å, respectively. The observed SASA fluctuation is associated with the rearrangement of the amino acid residues from either the accessible region or the buried region. This may contribute to alterations in the protein-ligand overall structure.

Materials and Methods
The computational tools used for this study were part of the Schrödinger Platform (version 2020-3 for Windows 10).

Preparation of Ligands
A library of ten small molecular weight compounds was used in this in silico study, which was based on the bioactive compounds we previously identified and isolated from the oil of Monodora myristica [12]. The preparation of the ligands was carried out by using the Ligprep module of Maestro 20.3 with an OPLS3 force field at pH 7.0 ± 2.0 [37]. The "desalt" and "generate tautomers" options were selected, and the stereoisomer computation was left to generate at most 32 binding sites per ligand. The output format was left as that of Maestro's defaults.

Preparation of Proteins
The proteins of interest (Protein Data Bank IDs: 3ILW, 1KIJ, and 6RKU) were retrieved from the Research Collaboratory for Structural Bioinformatics database (http: //www.rcsb.org/pdb accessed on 2 November 2022) and were uploaded to the workspace of Maestro 20.3. The proteins were prepared via the protein preparation wizard of the Schrödinger platform. During the pre-processing of the proteins, bond orders were assigned, water molecules were deleted within 2.8 Å from het groups, and het states were set at pH 7.0 ± 2.0 [38]. Subsequently, hydrogen bonds were added, and ions were removed. In the "Refine" tab, the H-bond network was optimized by using PROPKA, and water molecules with less than three H-bonds to non-water molecules were removed [39]. Energy minimization was carried out by using the OPLS3 force field with an RMSD at 0.30 Å.

Receptor Grid Generation
The grid file receptor was generated by using the receptor grid generation panel that represents the active sites of the receptor for the simulation of glide ligand docking. The ligand-binding site was defined by the co-crystalized ligand within the workspace or the sitemap, a module of the Schrödinger platform. The van der Waals radii of the receptor atoms with partial atomic charge were set at a scaling factor of 1.0, and a partial cut-off of 0.25 was applied in order to soften the potential for the nonpolar parts of the receptor.

Molecular Docking Analysis
The compounds prepared from PubChem were docked into the active sites of the protein by using extra precision with the ligand sampling set as "non-refined". Prior to the docking of the prepared compounds with said protein, the co-crystallized ligand was docked into the active site of the protein to predict binding affinity and molecular interaction.

Predictions of ADME Properties and Toxicological Potential
The Absorption, Distribution, Metabolism, Excretion (ADME), the toxicity of the compounds, and the molecular properties of the lead compounds were predicted by using the Qikprop, the Maestro 20.3, and the AdmetSAR 2.0 web-based tools, respectively [40].

Binding Free Energy Calculation
The free energy binding of the protein-ligand complexes was accessed in order to determine the stability of these complexes via the Prime MM-GBSA program (Schrödinger suite Maestro 20.3). The Prime MM-GBSA panel was used in order to calculate the binding free energy for the ligands' protein complexes by using the MM-GBSA technology available with Prime [41,42]. The OPLS3 force field was selected, and VSGB was used as the continuum solvent model. All other options were set as default.

MD Simulations and Trajectory Analysis
The MD simulation for the native protein (1KIJ) and the three complexes (1KIJlevofloxacin, 1KIJ-27447, and 1KIJ-529885) was carried by using Schrödinger 2022 version 1 with Maestro version 13.1.137, MM share version 5.7.137, and Windows-x64 Platform. The system setup and the MD preparation and trajectory analysis methods were similar to those of a previous method of ours, with slight modification [28,29]. Briefly, the docked complexes were individually subjected to MD by using the Desmond module of the Schrödinger software with an OPLS 2005 force field. The protein-ligand complex was bounded with a predefined transferable intermolecular potential with a 3-point water model in an orthorhombic box. The volume of the box was minimized, and the overall charge of the system was neutralized by adding sodium and chloride ions to mimic physiological conditions. The temperature and the pressure were kept constant at 310 Kelvin and 1.01325 bar, respectively, by using a Nose-Hoover thermostat and a Martyna-Tobias-Klein barostat made from United State. The simulation relaxation was undertaken by using an NPT ensemble after considering the number of atoms, the pressure, and the timescale. During the MD simulation, the long-range electrostatic interactions were calculated by using the particle mesh Ewald method. In addition, the MD simulation was carried out for 100 ns, and the trajectory sampling was set at an interval of 100 ps with 1000 frame numbers. Simulation outputs were analyzed and visualized by a simulation interaction diagram and an MS-MD trajectory analysis. The MD analysis was done in replicate to avoid variation. Data were plotted by using OriginPro version 9.

Conclusions
The resistance of bacteria to several antimicrobial agents has become an issue of major concern with potentially devastating consequences if not approached seriously. The herein identified hit candidate displayed satisfactory pharmacokinetic properties, exhibited interesting binding affinities, and did not violate Lipinski's ROF; thereby, the results indicate a safe treatment option for bacterial infections. Moreover, our computational analyses provide the theoretical inhibitory activities of the identified molecules against DNA gyrase.