Exploring the anticancer and antibacterial potential of naphthoquinone derivatives: a comprehensive computational investigation

This study investigates the potential of 2-(4-butylbenzyl)-3-hydroxynaphthalene-1,4-dione (11) and its 12 derivatives as anticancer and biofilm formation inhibitors for methicillin-resistant staphylococcus aureus using in silico methods. The study employed various computational methods, including molecular dynamics simulation molecular docking, density functional theory, and global chemical descriptors, to evaluate the interactions between the compounds and the target proteins. The docking results revealed that compounds 9, 11, 13, and ofloxacin exhibited binding affinities of −7.6, −7.9, −7.5, and −7.8 kcal mol−1, respectively, against peptide methionine sulfoxide reductase msrA/msrB (PDB: 3E0M). Ligand (11) showed better inhibition for methicillin-resistant staphylococcus aureus msrA/msrB enzyme. The complex of the 3E0M-ligand 11 remained highly stable across all tested temperatures (300, 305, 310, and 320 K). Principal Component Analysis (PCA) was employed to evaluate the behavior of the complex at various temperatures (300, 305, 310, and 320 K), demonstrating a total variance of 85%. Convergence was confirmed by the eigenvector’s cosine content value of 0.43, consistently displaying low RMSD values, with the minimum observed at 310 K. Furthermore, ligand 11 emerges as the most promising candidate among the compounds examined, showcasing notable potential when considering a combination of in vitro, in vivo, and now in silico data. While the naphthoquinone derivative (11) remains the primary candidate based on comprehensive in silico studies, further analysis using Frontier molecular orbital (FMO) suggests while the Egap value of compound 11 (2.980 eV) and compound 13 (2.975 eV) is lower than ofloxacin (4.369 eV), indicating their potential, so it can be a statement that compound 13 can also be investigated in further research.


GRAPHICAL ABSTRACT 1 Introduction
Naphthoquinones have emerged as a pivotal compound class in drug development, celebrated for their versatility in synthesizing bioactive compounds.These compounds demonstrate a range of biological activities, including antimicrobial, antitumor, antioxidant, antimalarial, and neuroprotective effects, primarily attributed to their distinctive redox properties (Fiorito et al., 2016;Malik et al., 2021).Their electron-accepting capabilities form highly reactive radicals that interact with crucial biological molecules like DNA, enzymes, and proteins, a central aspect contributing to their efficacy in various therapeutic applications.The escalating global health crisis posed by multidrug-resistant bacteria, particularly Methicillinresistant staphylococcus aureus (MRSA), has spurred research into alternative therapeutic agents.Notably, naphthoquinones, especially 1,4-naphthoquinone and its derivatives, exhibit promise in this context.Their antibacterial potential against resistant strains underscores their significance in current medical research, providing an innovative approach to address this urgent health concern (Stefani et al., 2012;Coenen et al., 2013).
Our research concentrates explicitly on the synthesis and assessment of derivatives of lawsone, a naturally occurring hydroxynaphthoquinone within the naphthoquinone family.These derivatives manifest a spectrum of biological activities, encompassing antifungal, antitumor, and antiviral properties (Deurenberg and Stobberingh, 2008;Stefani and Goglio, 2010).Their mechanisms of action as antibacterial agents are diverse, involving the disruption of protein and nucleic acid synthesis, modulation of redox processes, and induction of apoptosis in bacterial cells (Loomba et al., 2010;Rossi et al., 2014;Wang et al., 2022).The wide range of activities underscores the potential of lawsone derivatives as versatile therapeutic agents.In the field of cancer treatment, naphthoquinones have exhibited significant potential.They have been identified to target multiple molecular pathways in cancerous cells, leading to apoptosis and inhibition of cell growth.This mode of action, combined with their capacity to inhibit crucial cellular processes selectively, positions them as promising candidates in oncology (Harada et al., 2014;Brixius-Anderko and Schott, 2019).The versatility of naphthoquinones is further emphasized by their capability to interfere with the electron transport chain in both eukaryotic and prokaryotic cells, highlighting their potential to address a broad range of diseases.Figure 1 illustrates the chemical structures of the 13 naphthoquinone derivatives which has been utilized, obtained from Song R et al. (2020).
Alongside lab experiments, computer-based studies have been vital in exploring the potential of naphthoquinone derivatives.These virtual investigations give us insights into how naphthoquinones interact with different biological targets, helping us identify promising leads for further experiments (Goldberg et al., 2015;Jahanbad and Esfahlan, 2017;Khosla et al., 2022).Our study uses these computer techniques to evaluate new naphthoquinone derivatives.We analyze their physical and chemical properties, how they move through the body, and any potentially harmful effects, aiming to find candidates with the best antibacterial and anticancer properties (Sahoo et al., 2022;Rahman et al., 2022).Uddin et al.'s research highlights certain chemical derivatives' potential cancer-fighting abilities (Esha et al., 2023;Uddin et al., 2023).The first study (Uddin et al., 2023) looked into benzylidene malononitrile and ethyl 2-Cyano-3-phenylacrylate derivatives, while the next one (Esha et al., 2023) focused on fluoro flavonoid derivatives.Both studies provide valuable insights into the anticancer properties of these compounds.Using computer simulations, the second study (Esha et al., 2023) gives a computational perspective that can guide more experiments and potentially have applications in medicine and biology.
In this study, we performed in silico investigations to assess the potential of recently developed naphthoquinone derivatives (1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(13) in comparison to established pharmaceutical compounds-specifically, ciprofloxacin (D1), ofloxacin (D2), and vancomycin (D3) (Figures 1,  2) as potential agents with antibacterial and anticancer properties.Various computational techniques, including molecular docking, PASS predictions, molecular dynamics simulations, and DFT calculations, were utilized to scrutinize their physicochemical properties, pharmacokinetic and pharmacodynamics profiles, and potential toxicological effects.Additionally, molecular dynamics simulation analyses were conducted on the active site of proteinligand interactions to evaluate the stability of the protein-ligand complexes.The outcomes of these studies suggest that these derivatives represent a promising novel class of anticancer and antibacterial agents, emphasizing the need for further exploration in drug development.In summary, naphthoquinones constitute a significant and versatile class of compounds in the pharmaceutical landscape, showing promise in combating multidrug-resistant bacteria and various types of cancer.Our research contributes to this expanding field by synthesizing and evaluating new derivatives, aiming to develop effective therapeutic agents against these challenging medical conditions.

Molecular docking 2.3.1 Preparation of ligands
The reference pharmaceutical compounds, specifically ciprofloxacin, ofloxacin, and vancomycin, were obtained from the PubChem database in SDF (structure-data file) format (Kim et al., 2016).The naphthoquinone derivatives' three-dimensional (3D) molecular structures were also acquired in SDF format from PubChem.Subsequently, the structural optimization of all naphthoquinone derivatives was conducted using the B3LYP/6-31G (d,p) computational method with the Gaussian 16 software package.The visualization and representation of these naphthoquinone structures were achieved using GaussView 6 software.To ensure the geometric stability and conformational accuracy of the ligands chosen for molecular docking, energy minimization (EM) procedures were applied.Following this, the conversion of these optimized structures into the PDBQT (Protein Data Bank, Partial Charge, and Atom Type) format was carried out using the OpenBabel plugin integrated within PyRx 0.8 software, available at https://pyrx.sourceforge.io/(Dallakyan and Olson, 2015).This involved retrieving reference pharmaceutical compounds and naphthoquinone derivatives in specific file formats, optimizing their structures, and converting them for subsequent molecular docking studies.
(A) Residues in interaction within 3E0M, (B) the distribution of non-covalent interactions, and (C) a map illustrating the interaction between residues in 3E0M-naphthoquinone derivatives 9, 11, 13, and the reference drug of ofloxacin complex.
edu/ (Laskowski et al., 2013).Additionally, ProSA was utilized to calculate the Z-scores for these proteins (Rathod et al., 2022).To ensure optimal performance in molecular docking studies, structural optimization of the selected proteins was carried out using Chimera version 1.16 (Goddard et al., 2007).The Chimera 1.16 Dock Prep tool was used to prepare the proteins for docking, involving the addition of a +2 charge to the heme group, assignment of charges to standard residues using AMBER ff14SB, utilization of Gasteiger charges for other residues, and addition of hydrogen atoms.The preparation of peptide methionine sulfoxide reductase msrA/msrB (PDB: 3E0M) followed a similar process, except not removing any chains.

Protein-ligand docking
In our study, we employed AutoDock Vina software for proteinligand docking.We focused our docking investigations on the target protein and various ligands.To encompass the entire protein, we established a grid box centered on it.We must note that we ensured the stability of both proteins and ligands throughout the entire docking process.Additionally, we validated the protein-ligand complexes by conducting re-docking experiments.To identify the amino acids interacting with the ligands, we employed UCSF Chimera version 1.16 (Goddard et al., 2007).Three-dimensional structures for molecular docking images were generated using Pymol version 2.5 and Chimera 1.16.Finally, the BIOVIA Discovery Studio (Baroroh et al., 2023) visualizer facilitated the visualization of binding modes in protein-ligand interactions, enabled the examination of 2D protein interactions with the ligands, and assisted in determining hydrogen density around residues engaged in interactions with the protein.

Molecular dynamics simulation
Molecular dynamics simulations were carried out using Galaxy Europe to investigate the interactions between ligands 9, 11, 13, and ofloxacin as reference drugs and the chain of 3E0M.These simulations were conducted with GROMACS version 2021.6 (Van Der Spoel et al., 2005) and the AMBER99SB force field (Hansson et al., 2002;Showalter and Brüschweiler, 2007), offering a detailed account of atomic interactions.Molecular dynamics simulations are a robust computational technique, surpassing less computationally intensive docking methods.They provide high precision and offer insights into complex system behavior, often inaccessible through experimental means (Hansson et al., 2002;Showalter and Brüschweiler, 2007).In this study, we performed additional MD simulations for compounds 9, 11, 13, and ofloxacin with the A chain of 3E0M to validate previous findings.The protein topology parameters were generated using the Galaxy European Server (Bray et al., 2020), a widely used molecular modeling and simulation platform.The simulation parameters included the TIP3P water model, and hydrogen atoms were initially excluded from the GROMACS setup.Hydrogens were added to the ligand at pH 7.4, and during MD topology generation, the molecule's charge was maintained at 0 with a multiplicity of 1.The gaff force field was applied for parameterization.The ligand and protein files were combined, and a structural configuration was established within a 1 nm triclinic box.The system was solved with SPC water molecules in a triclinic box, and sodium and chloride ions were introduced to attain standard salt concentrations and neutralize the system (Presti et al., 2016).
An equilibration process was conducted to ensure system stability, utilizing position-restrained dynamics (NVT) at 300 K for 3,000 ps, implemented with the leapfrog algorithm (Cuendet and van Gunsteren, 2007;Esha et al., 2023;Uddin et al., 2023).Following this equilibration phase, the system underwent a production run for an additional 3,000 ps under constant temperature and pressure conditions.Subsequently, the system was further simulated for 20 ns at the same temperature and pressure.Various GROMACS utility tools, including gmx rmsd, gmx gyrate, gmx rmsf, and 'gmx hbond', were employed to generate graphs for RMSD, radius of gyration, RMSF, and hydrogen bonds of the ligands at 300 K over the 20 ns period.Furthermore, principal component analysis (PCA) was utilized to assess the stability of protein complexes, proving highly advantageous (Biswal et al., 2023;Panda et al., 2023).PCA was carried out for the ligand-protein complexes (9, 11,13, and ofloxacin) at 300 K, as well as for compound 11 at temperature settings of 300 K, 305 K, 310 K, and 320 K, using the Bio3D package of GALAXY Europe server (Grant et al., 2006;Afgan et al., 2018;Kumar et al., 2022).The cosine content of the ligands was also evaluated.These MD simulations have provided valuable insights into the behavior of the protein-ligand complexes, aiding in a better comprehension of the underlying physical principles governing the structure and function of biological macromolecules.

Analysis of frontier molecular orbitals (FMO)
Evaluating naphthoquinone and its derivatives (1-13) through fragment molecular analysis (FMO) has provided valuable insights into their chemical reactivity and stability.A broader energy gap (Egap) indicates heightened stability and reduced reactivity, whereas a narrower gap signifies softness, corresponding to increased reactivity and diminished stability.In addition to Egap, the energy levels of the Highest Occupied Molecular Orbital (HOMO) and Lowest Unoccupied Molecular Orbital (LUMO) play a significant role in understanding the electron-donating and -accepting characteristics of these compounds.To elaborate further, compounds with a narrow  3 and Supplementary Figures S1-S15.Table 1 indicates that compounds 9, 10, 11, 12, 13, and D2 (ofloxacin) exhibited energy gaps (Egap) of 3.202 eV, 2.998 eV, 2.980 eV, 2.977 eV, 2.975 eV, and 4.369 eV, respectively.Among these, compounds 9 to 13 displayed low energy gaps, suggesting good reactivity and stability, requiring less energy for electron promotion from HOMO to LUMO orbital.Compared to D2, their energy gaps were lower, indicating greater stability and a willingness to participate in reactions.Compound 9 had the highest Egap of 3.202 eV, signifying lower reactivity but increased hardness and stability.Furthermore, compound 13 showed the lowest HOMO-LUMO gap (2.975 eV) with a dipole moment of 2.504 D, implying efficient biological activities.The energy gap values followed this order: compound 9 (3.202 eV) > compound 10 (2.988 eV) > compound 11 (2.98 eV) > compound 12 (2.977eV) > compound 13 (2.975eV).Compounds 9, 11, and 13 were chosen for further evaluation through additional computational analysis and molecular dynamics.

Molecular docking calculation
Molecular docking is a vital computational approach for evaluating the binding affinity between a ligand and a protein or receptor.This method offers valuable insights into how compounds interact with a protein's active site residues and influence cellular processes.In our study, we performed molecular docking experiments involving naphthoquinone derivatives (1-13) and three reference drugs-ciprofloxacin, ofloxacin, and vancomycin-against four specific proteins: MSR (3E0M), MERS-CoV (4RNA), EFGR (4UV7_L), and MDD (4DPT), as detailed in the table.To ensure the quality of the protein structure (PDB: 3E0M), we conducted various assessments, including PROCHECK, The progression of the root mean square deviation (RMSD) (A) and the root mean square fluctuation (RMSF) (B) for the combined docked complex involving the protein (PDB: 3E0M) and ligands 9 (black line), 11 (blue line), 13 (red line), and ofloxacin (green line) during the 20 ns MD simulation.
ERRAT, Verify3D, ProSA-web server, and ProSA-web energy server, as illustrated in Figure 4.These evaluations provided valuable insights into the protein's structure.The structure met high-quality standards, as evidenced by the Ramachandran plot, overall quality factor, Verify3D score, and ProSA-web energy plot.In summary, the results confirm the validation and excellent quality of the protein structure, providing a solid foundation for our molecular docking Studies.
Table 2 illustrates the docking results, indicating that compounds 1-13 demonstrate notably higher binding affinities with the 3E0M protein compared to the other proteins in the study.A recent study reported that compound 11 exhibited a lower minimum inhibitory concentration in vitro.Our findings further highlight that compound 11 displays the highest binding affinity among all tested proteins, mainly showing strong affinities for Peptide methionine sulfoxide reductase msrA/msrB (PDB: 3E0M) −7.9 kcal mol −1 , aligning with our desired outcome.Compound 11 also exhibits binding affinity, Mevalonate diphosphate decarboxylase (PDB: 4DPT) −7 kcal mol −1 , The middle east respiratory syndrome coronavirus (MERS-CoV) papain like protease (PLpro) (PDB: 4RNA) −6.2 kcal mol −1 , and Epidermal Growth Factor Receptor (PDB: 4UV7) −6.5 kcal mol −1 .Table 2 in the SI shows that in the evaluation of MSR (PDB: 3E0M), all reference drugs demonstrated lower docking values ranging from −7.3 to −8 kcal mol −1 , in contrast to compounds 9 (−7.6 kcal mol −1 ), 13 (−7.5 kcal mol −1 ), and 11 (−7.9 kcal mol −1 ).Additionally, the results indicate that most of the examined docking values for our compounds surpass those of the reference drugs used in antibacterial and cancer treatment.To summarize, compounds 11 and ofloxacin (−7.8 kcal mol −1 ) exhibit strong binding affinities for all targeted proteins, while some compounds, such as 7 (−5.5 kcal mol −1 ) and 8 (−5.8 kcal mol −1 ), 12 (−6.2kcal mol −1 ), exhibit lower binding affinities, possibly attributed to differences in their chemical structures.
When treated with naphthoquinone derivatives, the functions of key enzymes like MERS-CoV papain-like protease (PLpro), Epidermal Growth Factor Receptor (EGFR), and Mevalonate Diphosphate Decarboxylase are potentially altered, impacting their roles in disease processes.Naphthoquinone derivatives, known for their broad biological activities, could inhibit MERS-CoV PLpro, disrupting the virus's replication and ability to evade the immune system.This inhibition could be crucial in developing new antiviral strategies against MERS-CoV.In the case of EGFR, which plays a significant role in cancer cell proliferation, naphthoquinone derivatives might act as inhibitors, blocking EGFR's signaling pathways that lead to tumor growth and metastasis.This suggests a promising avenue for cancer therapy, particularly in EGFR-driven cancers.Lastly, targeting Mevalonate Diphosphate Decarboxylase, a key enzyme in cholesterol biosynthesis, with these derivatives could reduce cholesterol synthesis.This action holds potential in treating disorders related to cholesterol metabolism, such as cardiovascular diseases.Thus, applying naphthoquinone derivatives to these enzymes offers a multifaceted therapeutic approach, spanning antiviral, anticancer, and cholesterol-lowering effects, demonstrating the versatility of these compounds in disease treatment and management.
In compound 11 and its interaction with the 3E0M protein, four distinct interactions have been identified, as outlined in Table 3 and   FIGURE 9 The progression of the radius of gyration (Rg) (A) and the hydrogen bond count (HBs) (B) for the combined docked complex involving the protein (PDB: 3E0M) and ligand 11 during the 20 ns MD simulation.
Frontiers in Chemistry frontiersin.org13  -S18).These interactions include four hydrophobic interactions, specifically Pi-Alkyl interactions with MET A:116 and ARG A:73 on the residue and Pi-sigma interactions with LEU A:69.Additionally, one interaction is classified as a conventional hydrogen bond at the ALA A: 165 residue.Compound 13's protein-ligand interaction analysis with 3E0M revealed one conventional hydrogen interaction at the VAL A:465 residue, one halogen interaction with GLY A:367, and four Pi-Alkyl interactions involving residues ALA A:556, ARG A: 415, ALA A:366, and VAL A:606.Compound 9 exhibited a binding affinity for 3E0M (−7.4 kcal/mol) with seven interactions in total, as shown in Figure 5 and Supplementary Figure S16.These interactions included two conventional hydrogen bonds (ARG A:261, ILE A: 162), a Carbon Hydrogen Bond with ASP A:163, and Pi-Alkyl interactions involving MET A:116, ALA A:164, and ARG A:261.
The pivotal finding of this study highlights compound 11 as an exceptionally promising candidate for drug development due to its profound interaction with the 3E0M protein (Figure 7).Compound 11 demonstrates superior binding affinity to the 3E0M protein compared to other compounds and reference drugs like ofloxacin Temperature (K) versus Time (ps) plots for (A) the combined system at four different temperature conditions (300, 305, 310, and 320 K), and (B) RMSD for the combined docked complex involving the protein (PDB: 3E0M) and ligand 11 (300 K: black line, 305 K: red line, 310 K: blue line, 320 K: green line) throughout the 20 ns MD simulation.and ciprofloxacin, engaging four distinct amino acid residues.This comprehensive and specific binding profile underscores the remarkable affinity of compound 11 for the target protein.Based on the molecular docking and interaction analyses, compound 11 exhibits heightened binding affinity and forms a stronger interaction with the target protein compared to other compounds.This indicates its potential as a more effective inhibitor of proteins 3E0M and 4DPT, which are relevant in anti-cancer treatment and antibacterial activities.Consequently, this in silico study suggests that compound 11 holds substantial promise as an inhibitor of bacterial cell membrane proteins.
Our findings reveal that, with the exceptions of compounds 8 and 10, most compounds exhibited moderate blood-brain barrier (BBB) penetration values, falling within the range of 0.625-1.Ciprofloxacin and ofloxacin displayed medium absorption values of 0.714 and 1, respectively, while vancomycin exhibited low absorption at −0.914.Moreover, the majority of naphthoquinone derivatives demonstrated favorable synthetic accessibility, ranging from 1 to 3.68, comparable to standard drugs.Compound 11, with a synthetic accessibility score of 3.97, and ofloxacin, with a score of 2.44, exemplified these characteristics.Our predictions regarding the potential risk of hERG activity inhibitors, with values ranging from −0.4 to −0.82 (Shadrack and Ndesendo, 2017), indicated that compound 11 achieved the desired hERG pIC50 value of −0.6441.In contrast, all other compounds remained within the reference range.These results affirm that all compounds (1-13) adhere to ADMET standards, signifying drug-likeness (bioavailability), and underscore their potential for novel drug development.Specifically, compound 11 exhibits favorable properties across all ADMET parameters, including moderate BBB penetration, high intestinal absorption, and a low risk of hERG inhibition.These characteristics align with the desired profile for a potential therapeutic agent.

Molecular dynamics simulation
We conducted molecular dynamics (MD) simulations spanning 20 ns for compounds 9, 11, 13, and ofloxacin in conjunction with the MSR protein (PDB: 3E0M).Previous molecular docking investigations revealed robust binding affinities of −7.6, −7.9, −7.5, and −7.8 kcal/mol for compounds 9, 11, 13, and ofloxacin, respectively, with the 3E0M protein, indicating their potential as promising drug candidates.The MD simulations aimed to assess the stability and interactions within the protein-ligand complex throughout the 20 ns trajectory.Various parameters, including RMSD (root-mean-square deviation), RMSF (root-mean-square fluctuation), Rg (radius of gyration) values, potential energies, temperature, hydrogen bonding, and principal component analysis (PCA), were employed for the simulation analysis.Additionally, the complex underwent simulations at four distinct temperatures (300, 305, 310, and 320 K) to explore its structural alterations under varying thermal conditions.This approach provided insights into the stability and adaptability of the protein-ligand complex under different physiological conditions.MD simulation evolution of RMSD for graphs (A) merged protein-ligand complex of ligands 9,11,13 and ofloxacin between PDB protein: 3E0M and modeled ligands 9 (light red line), 11 (grey line), 13 (yellow line) and reference drug, ofloxacin (purple line), (B) Protein-ligand complex of ligand 11 during 200 ns MD simulation.
The analysis of RMSD values obtained from the MD simulations provides insights into the stability and conformational dynamics of ligand-protein complexes (refer to Figure 8A and Supplementary Figures S19-S21).Notably, the complex involving compound 11 and protein 3E0M displayed the highest conformational stability, with RMSD values ranging from 0.2 to 0.55 nm (blue line), indicating minimal deviation from the initial structure.The complex of compound 13 and 3E0M exhibited the secondhighest stability, showing RMSD values from 0.3 to 0.6 nm, although it demonstrated larger fluctuations, reaching up to 0.7 nm.In contrast, the complexes involving ofloxacin and 3E0M, as well as compound 9 with 3E0M, displayed lower stability, evident from increased scattering in their RMSD plots, with values spanning from 0.2 to 0.7 nm and 0.1 to 0.8 nm, respectively.
In addition to RMSD, root-mean-square fluctuation (RMSF) assessment was utilized to gauge the flexibility of individual amino acid residues within the protein-ligand complexes throughout the MD simulations.RMSF values offer insights into the dynamic behavior of specific amino acids, indicating their level of flexibility or rigidity.Analyzing the RMSF values for each compound interacting with protein 3E0M revealed variations in residual flexibility among the complexes.The observed fluctuations ranged approximately from 0.04 to 0.13 nm for ligand 11, 0.03-0.9nm for the target protein 3E0M, and 0.05-0.9nm for the protein-ligand complex (see Figure 8B and Supplementary Figures S22-S24).These values suggest that a majority of amino acid residues in the complex exhibited fluctuations.Notably, regions directly interacting with ligand 11 exhibited low fluctuations, while other regions containing amino acids without direct interactions displayed more significant fluctuations.The complex involving protein 3E0M and ligand 11 demonstrated the highest RMSF values, indicating greater flexibility compared to other complexes (refer to Figure 8B and Supplementary Figures S22-S24).This implies that ligand 11 may induce a more dynamic conformation in the protein compared to the other ligands.The complex of protein 3E0M with ligand 13 also displayed moderate flexibility, while the complexes of protein 3E0M with ligand 9 and ofloxacin-3E0M exhibited relatively higher rigidity.
The variations in RMSF values suggest diverse flexibility within the protein-ligand complexes, potentially influencing their biological functions.The analysis of RMSF provides valuable insights into the structural dynamics of protein-ligand interactions, complementing the information derived from RMSD analysis.These findings indicate the stability of both the ligand and the protein throughout the simulation, with minor deviations in atomic positions.Notably, the observed fluctuations in specific ligand atoms may carry significance, potentially related to particular functional or interactional aspects of the protein.However, further investigation is necessary to ascertain these fluctuations' precise implications and potential impacts on the biological activities associated with 3E0M.In this computational study, it has been assessed that the radius of gyration (Rg) is a parameter used to characterize the size and shape of macromolecules like protein-ligand complexes.The Rg value determined for the rotein-ligand complex ranged from 2.12 to 2.45 nm, as depicted in Figure 9A (see Supplementary Figures S25-S28).Notably, all findings indicate that including the ligand did not induce substantial alterations in the overall size or shape of the protein.This outcome aligns with the notion that the binding of the ligand to the protein occurred in a manner that preserved the protein's tertiary structure without causing significant disruptions.Furthermore, the comparison in Rg values between the protein and the protein-ligand complex suggests that the binding event did not bring about notable structural modifications within the protein.Nevertheless, it is crucial to acknowledge that while Rg values provide insights into a protein's general size and shape, they may not unveil the precise molecularlevel structural alterations that have transpired.Further experiments, such as X-ray crystallography or NMR spectroscopy, may be indispensable for a more comprehensive understanding of the molecular interactions between the protein and the ligand.The results of the MD simulations emphasize the significance of the number of hydrogen bonds (HB) formed between the protein and the ligand-a critical factor in sustaining the stability of the proteinligand system.These findings are depicted in the time-dependent assessment of intermolecular hydrogen bonds during the 20 ns simulation, as shown in Figure 9B and Supplementary Figure S29.The simulation conducted at the standard body temperature exhibited the most consistent HB profile, with the number of hydrogen bonds fluctuating between zero and four throughout the entire 20 ns simulation period.Notably, between the 5 ns and 15 ns timeframe (Figure 9B), ligand 11 demonstrated the highest occurrence of hydrogen bond formations with the residues of the target protein (PDB: 3E0M).Initially, the analysis of interactions within the docked complex did not clearly indicate the presence of hydrogen bonds, as seen in Figure 8B.However, the subsequent 20 ns MD simulation revealed the emergence of hydrogen bond interactions for ligands 9, 11, 13, and the reference ligand (ofloxacin).
Upon comparing the outcomes of simulations involving compounds 9, 11, 13, and ofloxacin with the MSR protein (PDB: 3E0M), it is evident that compound 11 establishes a more robust interaction with 3E0M.This suggests its potential as an inhibitor of the MSR protein, which holds significance due to its crucial role in cancer and bacterial inflammation.Previous in vitro and in vivo experiments have underscored the promise of compound 11 as a candidate for an anti-cancer and antibacterial drug.This is attributed to its potent anti-mycobacterial activity, low cytotoxicity, and unique lipophilic vehicle feature.The complex comprising compound 9 and MSR showed HOMO-LUMO combined results slightly lower than ofloxacin, which exhibited the most pronounced negative binding affinity and met all criteria for consideration as a drug candidate, as substantiated by the molecular simulation results at 310 K. To further corroborate our identification of compound 11 as a promising anticancer and antibacterial drug candidate, additional molecular dynamics simulations were performed at three different temperatures: 300, 305, 310, and 320 K for the complex of ligand 11 with 3E0M, as depicted in Figure 10A.
The RMSD analysis of the ligand 11-MSR protein complex revealed its optimum stability at 310 K, with minimal deviation and a peak value of 0.8 nm.At 305 K, the complex exhibited the secondhighest stability, with RMSD values fluctuating between 0.1 and 0.4 nm and peaking at 0.5 nm.In contrast, the complex demonstrated reduced stability at 320 K and 300 K, manifesting more pronounced deviations.At these temperatures, the RMSD values ranged from 0.2 to 0.4 nm and 0.1 to 0.7 nm, respectively (see Figure 10B and Supplementary Figures S30-S44).Furthermore, the system exhibited sustained stability at different temperatures (300, 305, 310, and 320 K) over the 20 ns simulation, as depicted in Figure 10A.Additionally, conformational principal component analysis (PCA) based on Cα atoms was applied to the molecular dynamics (MD) trajectories of the protein (PDB: 3E0M) and ligand 11 at various temperatures (refer to Figure 11 and Supplementary Figures S45-S48).This study utilized PCA to evaluate the protein's variance, collective motions, and conformational changes during the MD simulations, as shown in Figure 11 (Table 7 and Supplementary Figure S17).The PCA analysis was performed on the MD trajectories of the target protein (ID: 3E0M) and the ligand 11 complex at 300, 305, 310, and 320 K, using the Bio3D package.The resulting plots illustrated eigenvalues versus eigenvectors, highlighting the primary motions extracted from the trajectories, with a focus on the first three eigenvectors (PC1, PC2, and PC3) and color dots representing the variance captured by the eigenvectors.
Table 7 provides a concise overview of the primary motions observed in the (ligand 11-3E0M) complex at various temperatures, specifically focusing on a subset and comparing the first three eigenvectors (PC1, PC2, and PC3).Notably, the complex simulated at 310 K (physiological body temperature) exhibited the most significant variation in PC1 (24.39) and the least in PC2 (19.19) and PC3 (14.14) concerning internal motion within the MD trajectory.This implies a robust interaction between compound 11 and 3E0M, a finding supported by both in vivo and in vitro investigations.However, further clinical applications are necessary to validate the efficacy and safety of compound 11 as an MSR-3E0M inhibitor.
In previous in vitro investigations, ligand 11 demonstrated significant potency and specificity as an inhibitor of biofilm formation and various bacterial diseases.Moreover, this ligand shows promise as a potential cancer treatment.Docking results revealed that ligand 11 exhibited the highest binding affinity (−7.9 kcal mol −1 ) with MSR-3E0M, while its minimum inhibitory concentration (MIC) values against MSSA, MSRA-ATCCBAA1717, and MSRA-ATCCBAA-44 (Table 8) ranged from 1.25 to 2.5 (Song et al., 2020).
In comparison, ofloxacin, with a binding affinity of −7.3 kcal mol −1 , displayed MIC values ranging from 0.25 to 16. Ofloxacin is a promising candidate for treating cancer and bacterial diseases, given its efficacy and selectivity against methicillin-resistant staphylococcus aureus (MRSA).It targets the enzyme responsible for biofilm formation, causes cell membrane damage, chelates intracellular iron ions, and generates intracellular reactive oxygen species.Ligand 11 also demonstrated stable binding to methicillin-resistant staphylococcus aureus, with consistent RMSD and RMSF results.Additionally, it maintained a stable conformation at human body temperature.Ligands 9 and 13 also exhibited potential as inhibitors of bacterial diseases, boasting lower Egap values and higher docking affinity than ofloxacin.These ligands formed interactions with amino acid residues that had positive results in previous studies.Particularly, ligand 11 showed the highest binding affinity for the A chain of MSR-3E0M.
To validate our findings and determine whether 2-(4butylbenzyl)-3-hydroxynaphthalene-1,4-dione ( 11) is optimal, we extended our investigation to include two additional naphthoquinone derivatives (9 and 13) and the standard drug, ofloxacin, which exhibited the most favorable docking results.Subsequently, we conducted 200 ns MD simulation revealed molecular dynamics simulations, and the outcomes are presented in Figure 12, comprising three RMSD panels: a) Merge proteinligand complex of ligands 9, 11, 13 and ofloxacin, b) protein-ligand complex 11 and c) Protein-ligand complex of ofloxacin.Figure 12 shows the root mean square deviation (RMSD) analysis for the three naphthoquinones and the standard drug.In Figure 12A) the yellow line corresponds to ligand 11, while the grey, red, and purple lines represent ligands 13, 9, and the standard drug, ofloxacin, respectively.Naphthoquinone derivative 11 is the only one exhibiting an RMSD value close to 0.35 nm, whereas the others collectively yield a value over 0.60 nm except the reference drug ofloxacin (0.40 nm).This notably higher combined value indicates that these ligands possess dissimilar shape orientations and exhibit less structural similarity.Notably, the naphthoquinone derivative (11) displays a lower RMSD value than the other compounds, as evident in Figures 12A, B. This underscores its unique structural characteristics and suggests its potential as the most promising candidate.For more detailed information, including individual graphs for each naphthoquinone with respect to the reference molecule (PDB: 3E0M), please refer to the supplementary information (Supplementary Figure S45).
Figure 13 illustrates the merged radius of gyration for three naphthoquinones and the standard drug.In Figure 13A, we observe a radius of gyration range of 0.32-0.35nm for naphthoquinone derivative 11, 0.37-0.41nm for the standard drug (ofloxacin), and 0.38-0.42nm for ligand 9 and 0.5-0.6 nm for ligand 13.It is important to note that higher values in the radius of gyration indicate lower structural similarity and fewer interactions among molecules.Remarkably, the ligand 11 exhibits the highest structural similarity and interactions among its molecules due to its lower radius of gyration.Detailed radius of gyration values and graphs for each naphthoquinone derivative can be found in the supplementary information (see Supplementary Figure S47).In conclusion, simulation results for ligand 11 provided valuable insights, indicating a robust interaction with 3E0M.These findings suggest favorable interactions with the protein, supported by both computational and laboratory experiments.Nevertheless, it is crucial to emphasize that further in vivo studies are necessary to validate the efficacy and safety of ligand 11 as an inhibitor of 3E0M.In conclusion, simulation results for ligand 11 provided valuable insights, indicating a robust interaction with 3E0M.These findings suggest favorable interactions with the protein, supported by both computational and laboratory experiments.Nevertheless, it is crucial to emphasize that further in vivo studies are necessary to validate the efficacy and safety of ligand 11 as an inhibitor of 3E0M.

Conclusion
Our in-depth computational analysis explored the medical applications of 13 recently synthesized naphthoquinone derivatives (1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(13)(14), focusing on their interactions with four essential target proteins.Using advanced techniques like FMO calculations, molecular docking, and dynamic MD simulations, our study revealed the exceptional binding affinity of compound 11 with protein 3E0M, surpassing in vitro results (Song et al., 2020).This underscores compound 11's considerable therapeutic potential, supported by substantial effects on 3E0M cell lines comparable to ofloxacin.All synthesized compounds met strict drug-likeness criteria, showcasing notable antibacterial and anticancer properties in computational assessments.Validating their pharmaceutical potential necessitates essential in vivo experiments, marking a significant stride toward innovative therapeutic solutions in the medical field.
) data: energy gap (Egap), ionization potential (IP), electron affinity (EA), binding energy (H), electronegativity (χ), chemical potential (μ), global hardness (η), softness (σ), electrophilicity (ω), and dipole moment (Debye) for compounds (1-13) at 298.15 K. a Ligand E LUMO (eV) E HOMO (eV) E gap (eV) IP (eV) EA (eV) M (eV) χ (eV) μ (eV) η (eV) σ (eV) ω (eV) (d,p) level of theory, yielded optimized structural parameters, as elaborated in Supplementary TablesS1-S16.In a prior investigation,Uddin et al. (2023) comprehensively explored the structural parameters, including bond lengths and angles, for seven compounds.The calculations were conducted at the B3LYP level of theory, utilizing various basis sets, namely, 6-31G(d,p),6- 311G(d,p), and 6-311++G(d,p).Their findings indicated that B3LYP/6-31G(d,p) exhibited good agreement with experimental results.Furthermore, Frontier Molecular Orbital (FMO) analysis was undertaken to elucidate the energy levels and spatial distributions of the HOMO (highest occupied molecular orbital) and LUMO (lowest unoccupied molecular orbital).These Frontier orbitals, positioned at the outermost boundary of the molecule, play a pivotal role in determining its electronic behavior.The HOMO energy reflects the molecule's electron-donating ability, while the LUMO energy characterizes its electron-accepting potential.An ideal scenario is characterized by a high HOMO value and a low LUMO value, indicating favorable electron transfer properties.Gaussview 6 software(Dennington et al., 2016) was employed to visualize the electrostatic potential distribution across the molecules.Molecular electrostatic potential (MEP), represented through a spectrum of colors, reveals the net electrostatic potential arising from the total charge distribution at nucleophilic and electrophilic sites within a compound.Electrophilic regions, indicating the strongest attraction, are depicted in blue, while nucleophilic areas, signifying the highest repulsion, are colored red.Green zones represent a neutral electrostatic potential.Additionally, the dipole moments of the compounds were calculated using Gaussview 6.The FMO energy gap, representing molecular stability, is determined by the difference between HOMO and LUMO FIGURE 3 (A) Molecular orbitals of isodensity surfaces (0.02 electrons Bohr −3 surface) (red = electron-rich, blue = electron-deficient) of HOMO and LUMO; (B) DOS plot and HOMO-LUMO energy gap; (C) Maps of electrostatic potential (0.02 electrons Bohr −3 surface) (red = electron-rich, blue = electrondeficient) for the compound 11.

FIGURE 4
FIGURE 4 Validation and quality evaluation of protein (PDB: 3E0M) by using (A) Rama-chandran plot analysis using the Procheck tool, (B) Z-Score prediction from the ProSA-web Server, (C) ERRAT, (D) Local model quality and (E) Varify3D plot showing amino acids in favored regions.

FIGURE 7
FIGURE 7 Molecular docking poses: (A) Protein ligand binding cavity; (B) Ligand placement in protein pocket; (C) Active site positioning; and (D) Visualization of hydrogen bonding in solid for the modeled 3E0M protein with ligand 11.

Figure 3
Figure 3 illustrates Molecular Electrostatic Potential (MEP) maps for compound 11.Electrophilic sites (depicted in blue) function as electron acceptors, while nucleophilic sites (shown in red) serve as electron donors.Partially nucleophilic regions are represented in yellow.The color-coded electrostatic potential highlights negative potential (in red, orange, and yellow shades) around electronegative atoms such as oxygen, indicating electrophilic reactivity.Conversely, positive potential (in blue) is observed over hydrogen atoms, signaling nucleophilic reactivity.Green areas denote the neutral potential.MEP maps provide insights into molecular stability by indicating consistent potential and are also helpful in studying molecular recognition and understanding interactions based on electrostatic complementarity.

Figures
Figures 5, 6 (Supplementary FiguresS16-S18).These interactions include four hydrophobic interactions, specifically Pi-Alkyl interactions with MET A:116 and ARG A:73 on the residue and Pi-sigma interactions with LEU A:69.Additionally, one interaction is classified as a conventional hydrogen bond at the ALA A: 165 residue.Compound 13's protein-ligand interaction analysis with 3E0M revealed one conventional hydrogen interaction at the VAL A:465 residue, one halogen interaction with GLY A:367, and four Pi-Alkyl interactions involving residues ALA A:556, ARG A: 415, ALA A:366, and VAL A:606.Compound 9 exhibited a binding

FIGURE 10
FIGURE 10 conform to Lipinski and Veber guidelines.Lipinski's criteria, crucial for preparing a compound for oral use, include five key conditions: a.Molecular weight (MW) of 500 g/mol or less b.Partition coefficient (log P) of 5 or less c.No more than 5 hydrogen bond donors (HBD) d.A maximum of

FIGURE 11
FIGURE 11Principal component analysis (PCA) of MD trajectories for the target protein (PDB: 3E0M) and ligand 11 complex at (A) 300 K, (B) 305 K, (C) 310 K, and (D) 320 K (Intermediate states are marked by white dots, energetically unstable conformations are represented by blue dots with scattering, and stable conformation states are denoted by red dots).

FIGURE 13 MD
FIGURE 13MD simulation evolution of Rg for merged graphs (A) ligands 9, 11, 13 and ofloxacin and (B) merged docked protein-ligand complex between the PDB protein: 3E0M and modeled ligands 11 and reference drug (ofloxacin) during 200 ns MD simulation.

TABLE 1 Molecular
orbital (MO

TABLE 6
In silico prediction of selected ADMET parameters for all compounds (1-13).

TABLE 8
Comparative analysis of the selectivity and potency of compounds of naphthoquinone derivatives (9, 11, 13, and ofloxacin), utilizing docking results and minimum inhibitory values.

TABLE 7
Variability in principal components revealed via PCA for the target protein (3E0M) and ligand 11 complex at various temperatures.