Computational Screening of T-Muurolol for an Alternative Antibacterial Solution against Staphylococcus aureus Infections: An In Silico Approach for Phytochemical-Based Drug Discovery

Staphylococcus aureus infections present a significant threat to the global healthcare system. The increasing resistance to existing antibiotics and their limited efficacy underscores the urgent need to identify new antibacterial agents with low toxicity to effectively combat various S. aureus infections. Hence, in this study, we have screened T-muurolol for possible interactions with several S. aureus-specific bacterial proteins to establish its potential as an alternative antibacterial agent. Based on its binding affinity and interactions with amino acids, T-muurolol was identified as a potential inhibitor of S. aureus lipase, dihydrofolate reductase, penicillin-binding protein 2a, D-Ala:D-Ala ligase, and ribosome protection proteins tetracycline resistance determinant (RPP TetM), which indicates its potentiality against S. aureus and its multi-drug-resistant strains. Also, T-muurolol exhibited good antioxidant and anti-inflammatory activity by showing strong binding interactions with flavin adenine dinucleotide (FAD)-dependent nicotinamide adenine dinucleotide phosphate (NAD(P)H) oxidase, and cyclooxygenase-2. Consequently, molecular dynamics (MD) simulation and recalculating binding free energies elucidated its binding interaction stability with targeted proteins. Furthermore, quantum chemical structure analysis based on density functional theory (DFT) depicted a higher energy gap between the highest occupied molecular orbital and lowest unoccupied molecular orbital (EHOMO-LUMO) with a lower chemical potential index, and moderate electrophilicity suggests its chemical hardness and stability and less polarizability and reactivity. Additionally, pharmacological parameters based on ADMET, Lipinski’s rules, and bioactivity score validated it as a promising drug candidate with high activity toward ion channel modulators, nuclear receptor ligands, and enzyme inhibitors. In conclusion, the current findings suggest T-muurolol as a promising alternative antibacterial agent that might be a potential phytochemical-based drug against S. aureus. This study also suggests further clinical research before human application.


Introduction
The population growth and invasion dynamics of S. aureus are undoubtedly opportunistic, as they subvert the human immune surveillance machinery by evoking the metastatic colony dissemination from superficial to life-threatening levels [1].The source of bacterial dissemination mainly evolved from the mucosal contacts, as the mucosal layer mainly favours harbouring of bacterial strains.Therefore, the mucosal surfaces of the nose, throat, vaginal wall, open wound, and gastrointestinal tract impart favourable colonization carriages of S. aureus bacteraemia (SAB) [2].Colonization is also frequent in human juveniles and patients who have already been vulnerable to Human immunodeficiency virus (HIV) and diabetes [3].Despite the backdrop of prevalence and prognostic outcomes of SAB in the last decade, it has been witnessed the emergence of multi-drug-resistant (MDR) strains to empiric therapies, antibiotics, and large-scale treatment options [4].Although, the present array of antibiotics, including methicillin, vancomycin, daptomycin, and linezolid, multi-drug-resistant S. aureus has emerged as the most formidable bacterial strain, causing life-threatening infections and illnesses in humans [5].The estimated numbers of invasive disease cases and morbidity-mortality in the United States due to methicillin-resistant S. aureus (MRSA) pathogenesis were 94,360 and 18,650 in 2005, respectively [2].Similarly, MRSA infection is associated with annual healthcare costs of about $3 billion a year and is set to increase over time [6].Additionally, SAB-associated infections are also reported in the context of immunotherapies and oncological treatments [7].The rampant and escalating severity of diseases, coupled with the failure of rational therapeutic approaches, due to the evolution of MDR S. aureus strains, might lead one to speculate that humans could become highly vulnerable and unmatched in facing the uncontrolled infectious threat posed by S. aureus [3].Therefore, considering the systemic toxicity and limited efficacy of current antibacterial medications, identifying the underlying mechanisms of this epidemic and developing alternative antibacterial drugs should be the prime goal of the current research.These new drugs should aim to achieve high success rates against S. aureus and prevent obsolescence, thereby mitigating the uncertainty of public health and improving patient care for SAB-associated infections.
Plant-derived phytochemicals can be adopted as sustainable alternative antibacterial agents.Over the last three decades, phytochemicals have undergone extensive experimental validation, establishing themselves as emerging antibacterial agents with significant potential [8].In the treatment of penicillin G-resistant strains of S. aureus, many phytochemicals have successfully passed experimental tests against various MDR strains [9].Phytochemicals-based antibacterials not only possess antibacterial activity but also act as an anti-inflammatory and antioxidant agent that can dramatically improve the targeting of the infection site while minimizing systemic exposure and associated toxicity [10].Therefore, different classes of phytochemicals, such as phenols, alkaloids, terpenoids, flavonoids, carotenoids, organosulfur compounds, and coumarins, etc., are actively available in the market as food and drug administration (FDA)-approved drugs and used in herbal medicinal products [11].Phytochemicals have been shown to impede major bacterial MDR factors, including bacterial virulence factors, ion pumps, bacterial cell wall toxins, cellular microstructure, replication machinery, and membrane permeability [12].Therefore, this study focuses on the in silico repurposing of a new generation of phytochemicals against S. aureus, aiming to explore alternative antibacterial therapeutic approaches.
T-muurolol is a cadinene sesquiterpenoid with versatile characteristics, as several studies have demonstrated its diverse biological activities, including antimicrobial, antioxidant, anti-inflammatory, and antitermitic effects [13][14][15].Although T-muurolol is a potential candidate with various reported antimicrobial activities and could serve as an alternative solution to combat S. aureus infections, the phytochemical-based drug-reprofiling approach is still in its early stages.Currently, there are no such research articles detailing the re-profiling of T-muurolol as an alternative treatment for suppressing SAB-associated infections, and its molecular interaction with S. aureus-associated proteins remains largely unexplored.Therefore, this study investigates the potential of T-muurolol as an anti-staphylococcal agent through molecular docking, pharmacokinetics, molecular dynamics, and chemical structure analysis.

Molecular Docking and Interaction Analysis
The molecular interactions between T-muurolol and tested bacteria-specific proteins for S. aureus, along with various proteins related to antioxidant activity and inflammation in humans, are summarized in Tables 1 and 2. T-muurolol exhibited a good binding affinity with all the proteins used in the study between −7.5 kcal/mol and −4.3 kcal/mol.T-muurolol demonstrated a strong binding affinity for S. aureus lipase (6KSI), showing a binding energy of −7.1 kcal/mol, which is notably higher than its affinity for other tested pathogenic proteins.A hydrogen bond of bond length 3.78 Ǻ with A chains of amino acid residue Asp236, along with non-weak interactions with A chains of Leu250, Lys249, and Tyr240, were found in this interaction (Figure 1C,D).
Table 1.Binding free-energy values (kcal/mol) of T-muurolol as ligand along with ascorbic acid, ibuprofen, and 2-oxazolidinone as references for good antioxidant, anti-inflammatory, and antibacterial agents, respectively.T-muurolol also demonstrated significant binding interactions (with a binding energy of −7.5 kcal/mol) with dihydrofolate reductase (3SRW), a key enzyme responsible for bacterial DNA synthesis.Notably, T-muurolol formed a hydrogen bond of 5.03 Å with amino acid Phe93, along with one group of electrostatic and two groups of hydrophobic interactions (Pi-Sigma, Alkyl, and Pi-Alkyl) with the amino acid residues of Ala8, Leu6, Leu21, Phe93, and Val32, which are mainly responsible for strong interaction (Figure 1A,B).Similarly, MDR S. aureus penicillin-binding protein 2a (1MWU) is a key drug target responsible for broad-spectrum beta-lactam resistance.T-muurolol exhibited a binding affinity of −6.1 kcal/mol with 1MWU, forming a 4.42 Å hydrogen bond with Glu294, along with electrostatic and hydrophobic interactions involving Ala276, Lys289, and Tyr272, contributing to strong binding interactions (Figure 2A,B).The robust interactions between PBPs and T-muurolol can effectively inhibit MRSA by disrupting peptidoglycan cross-linking and impeding bacterial cell wall synthesis.Likewise, D-Ala: D-Ala ligase (3N8D), a potential target for vancomycin-resistant Staphylococcus aureus (VRSA) due to its role in ligase impairment, is affected by T-muurolol, evidenced by significant binding energy with 3N8D (−5.9 kcal/mol).Furthermore, T-muurolol demonstrates substantial binding affinities with RPP TetM (3J25) with a binding energy of −6.6 kcal/mol.These proteins serve as key targets in the context of tetracycline-resistant Staphylococcus aureus (TetRSA), emphasizing the potential of T-muurolol as a therapeutic agent against resistant bacterial strains.Moreover, T-muurolol displayed favourable binding interactions with FAD-dependent NAD(P)H oxidase(2CDU) (−6.4 kcal/mol) among all the tested antioxidant proteins.Additionally, it demonstrated strong binding affinity with cyclooxygenase-2 (1CX2) (−6.3 kcal/mol), which was higher compared to the reference anti-inflammatory drug ibuprofen (Figure 3) (Table 1

Absorption
The water solubility of tested compounds, cell permeability utilizing the colon carcinoma (Caco-2) cell line, human intestinal absorption, skin permeability, and whether the molecule is a P-glycoprotein substrate or inhibitor are the primary parameters for evaluating drug absorption criteria.The current results suggested that T-muurolol had a lower solubility than water, whereas the reference antibiotic drug 2-oxazolidinone seemed to be highly soluble (Table 3).Caco-2 cell permeability influences the ultimate bioavailability; according to Chandra et al. [16], a drug with a value of > 0.90 is deemed highly permeable.The obtained results exhibited the high permeability of T-muurolol in the Caco-2 cell line with a value of 1.479.Similarly, all the tested compounds showed a higher gastrointestinal (GI) absorption percentage, which was more than 90% (Table 3).A study by Saha et al. [17] explained that the human intestine is primarily an active location for drug absorption, with more than 30% being considered rapidly absorbed.Our results also indicated that T-muurolol, along with reference drugs, is neither a P-glycoprotein substrate nor an inhibitor (Table 3).

Absorption
The water solubility of tested compounds, cell permeability utilizing the colon carcinoma (Caco-2) cell line, human intestinal absorption, skin permeability, and whether the molecule is a P-glycoprotein substrate or inhibitor are the primary parameters for evaluating drug absorption criteria.The current results suggested that T-muurolol had a lower solubility than water, whereas the reference antibiotic drug 2-oxazolidinone seemed to be highly soluble (Table 3).Caco-2 cell permeability influences the ultimate bioavailability; according to Chandra et al. [16], a drug with a value of > 0.90 is deemed highly permeable.The obtained results exhibited the high permeability of T-muurolol in the Caco-2 cell line with a value of 1.479.Similarly, all the tested compounds showed a higher gastrointestinal (GI) absorption percentage, which was more than 90% (Table 3).A study by Saha et al. [17] explained that the human intestine is primarily an active location for drug absorption, with more than 30% being considered rapidly absorbed.Our results also indicated that T-muurolol, along with reference drugs, is neither a P-glycoprotein substrate nor an inhibitor (Table 3).

Distribution
Here, drug distribution is generally analyzed based on the steady-state volume of distribution (VDss), and fraction of unbound value to human plasma.The steady-state volume of distribution (VDss) is a key pharmacokinetic parameter for determining drug dosages, indicating the theoretical volume needed to achieve the same blood plasma concentration as the administered dose.Higher VDss values suggest greater distribution into tissues rather than plasma.For antibiotics and antivirals, extensive tissue distribution is often desired.VDss is considered low if log(VDss) < −0.15 and high if >0.45 [17].T-muurolol showed a moderate VDss value of 0.42 log L/kg, which suggests its great distribution over tissue (Table 3).A drug's effectiveness is influenced by its binding to blood proteins; more binding facilitates membrane crossing [18].The fraction unbound to human plasma should be between 0.02 and 1.0 [18].T-muurolol showed good fraction unbound values with 0.28 (Table 3).

Metabolism
The cytochrome P450 (CYP) enzyme group includes isozymes that metabolize drugs, fatty acids, steroids, bile acids, and carcinogens.Drug metabolism depends on whether a substance is a CYP substrate or inhibitor.A study found that T-muurolol and the reference drugs tested are neither substrates nor inhibitors of CYP enzymes (Table 3), indicating they will be metabolized effectively and not hindered by the body's biological processes.

Excretion
Organic cation transporter 2 (OCT2) is a renal uptake transporter crucial for drug deposition and clearance from the kidneys [16].Excretion is assessed by total clearance and whether the compound is a substrate of OCT2 or not.The predicted results showed that T-muurolol is not a substrate of OCT2 (Table 3), suggesting it may be eliminated through a different route.Also, T-muurolol showed a moderate excretion clearance efficacy with a log(CLtot) value of 1.085 mL/min/kg.

Toxicity
The toxicity of the drug compounds was assessed using AMES results, human maximum tolerated dose, oral rat-acute toxicity, hepatotoxicity, skin sensitization, minnow toxicity, and human ether-a-go-go gene (hERG) inhibition.The AMES result indicated that T-muurolol is non-mutagenic and non-carcinogenic (Table 3).The maximum recommended tolerance dose (MRTD) estimates human toxic doses and is considered low if below log 0.477 (mg/kg/day) [17].So, T-muurolol showed less toxicity, with an MRTD value of log 0.343 mg/kg/day, whereas the reference drugs seemed to be highly toxic.The results also showed that none of the tested compounds are non-hERG inhibitors.However, T-muurolol was predicted to be skin-sensitive, although no hepatotoxic prediction was established.All of the tested compounds showed higher oral rat-acute toxicity (LD 50 ) (Table 3) which indicated less lethality compared to one with a lower LD 50 value [19].If a molecule's log LC50 (concentration causing 50% fathead minnow mortality) is below 0.5 mM (log LC 50 < −0.3), it is deemed highly acutely toxic.Current findings suggest that T-muurolol and reference drugs are less toxic, with significantly higher scores than the mentioned LC 50 threshold (Table 3).

Drug-Likeness and Bio-Activity Analysis
The drug-likeness property of the compounds is based on Lipinski's rule of five, which encompasses crucial molecular properties that impact a drug's pharmacokinetics in the human body.Current results showed that none of the tested compounds violated Lipinski's parameters, such as molecular weight, number of hydrogen bond donors and acceptors, and octanol-water partition coefficient (Table 4).The topological polar surface area (TPSA) calculation indicates a drug's bioavailability and its hydrogen bonding potential.All tested compounds have a TPSA range of 20.23-38.33Å (Table 4), significantly lower than the typical range of 160 Å [20].Our results also indicated an ideal bioavailability score of 0.55 for T-muurolol which will be well-absorbed by the human body [21].Bioactivity scores for the drug compound can be calculated based on parameters including kinase inhibition, protease inhibition, enzyme activity inhibition, binding to G protein-coupled receptors (GPCRs) and nuclear receptors, and ion channel modulation.The bioactivity score indicated that T-muurolol is highly active toward ion channel modulators, nuclear receptor ligands, and enzyme inhibitors with moderate activity toward GPCR ligands, but at the same time, it is inactive toward kinase inhibitors and protease inhibitors (Table 5).* GPCRs-G protein-coupled receptors, ICM-ion channel modulator, KI-kinase inhibitor, NRL-nuclear receptor ligand, PI-protease inhibitor, EI-enzyme inhibitor.A molecule with a bioactivity score above 0.00 is likely to exhibit significant biological activities; scores ranging from −0.50 to 0.00 are considered moderately active, while those below −0.50 are presumed inactive.

MD Simulation Analysis
Structural integrity, stability, compactness protein folding of apo-protein (ligand-free protein), and holo-proteins (proteins with ligands) had been critically analyzed based on some variable characteristics, namely RMSD (root mean square deviation), RMSF (root mean square fluctuation), R g (radius of gyration), hydrogen bonds (HB) formation, SASA (solvent accessible surface area), and binding free energy decomposition (MM-PBSA), after 100 ns of MD run.The superimposition of T-muurolol on the targeted proteins was shown in Figure 4 to observe the time-dependent protein-drug interaction.
nuclear receptor ligands, and enzyme inhibitors with moderate activity toward GPCR ligands, but at the same time, it is inactive toward kinase inhibitors and protease inhibitors (Table 5).* GPCRs-G protein-coupled receptors, ICM-ion channel modulator, KI-kinase inhibitor, NRLnuclear receptor ligand, PI-protease inhibitor, EI-enzyme inhibitor.A molecule with a bioactivity score above 0.00 is likely to exhibit significant biological activities; scores ranging from −0.50 to 0.00 are considered moderately active, while those below −0.50 are presumed inactive.

MD Simulation Analysis
Structural integrity, stability, compactness protein folding of apo-protein (ligand-free protein), and holo-proteins (proteins with ligands) had been critically analyzed based on some variable characteristics, namely RMSD (root mean square deviation), RMSF (root mean square fluctuation), Rg (radius of gyration), hydrogen bonds (HB) formation, SASA (solvent accessible surface area), and binding free energy decomposition (MM-PBSA), after 100 ns of MD run.The superimposition of T-muurolol on the targeted proteins was shown in Figure 4 to observe the time-dependent protein-drug interaction.

Monitoring Conformational Stability of the Backbone
The root mean square deviation (RMSD) of the six apo-protein and holo-protein system was illustrated in Figure 5, where the apo-protein and holo-protein system of 1CX2, 1MWU,

Residual Flexibility Analysis
The thermodynamic stability of protein-ligand complexes is determined by residual fluctuations (RMSF) of proteins in the presence of ligands.Comparing the RMSF of apoproteins and holo-proteins characterizes the conformational stability and flexibility of amino acid Cα atoms.Larger RMSF peaks in holo-protein residues indicate greater flexibility, while smaller peaks indicate rigidity, reflecting maximal and minimal interactions with their ligands and surrounding polar molecules, respectively [22]. Figure 6 gives us the rationale insights about the comparative analysis of apo-protein and holo-protein specific average RMSF values which ranges as 0.17 ± 0.02 nm for 1CX2, 0.26 ± 0.001 nm for 1MWU, 0.21 ± 0.04 nm for 3N8D, 0.13 ± 0.01 nm for 3SRW, 0.13 ± 0.005 nm for 6KSI, and 0.32 ± 0.24 nm for 3J25.From the average RMSF values of combined apo-protein and holo-protein systems, 3J25 holo-protein showed the largest increase in fluctuation, while 6KSI holo-protein

Residual Flexibility Analysis
The thermodynamic stability of protein-ligand complexes is determined by residual fluctuations (RMSF) of proteins in the presence of ligands.Comparing the RMSF of apoproteins and holo-proteins characterizes the conformational stability and flexibility of amino acid Cα atoms.Larger RMSF peaks in holo-protein residues indicate greater flexibility, while smaller peaks indicate rigidity, reflecting maximal and minimal interactions with their ligands and surrounding polar molecules, respectively [22]. Figure 6 gives us the rationale insights about the comparative analysis of apo-protein and holo-protein specific average RMSF values which ranges as 0.17 ± 0.02 nm for 1CX2, 0.26 ± 0.001 nm for 1MWU, 0.21 ± 0.04 nm for 3N8D, 0.13 ± 0.01 nm for 3SRW, 0.13 ± 0.005 nm for 6KSI, and 0.32 ± 0.24 nm for 3J25.From the average RMSF values of combined apo-protein and holoprotein systems, 3J25 holo-protein showed the largest increase in fluctuation, while 6KSI holo-protein showed the smallest change compared to their apo-proteins.This indicates that 3J25 maintained residual flexibility, whereas 6KSI retained rigidity in the presence of specific ligands.Among the six holo-proteins, 6KSI demonstrated the best interactions with ligands and polar molecules, while 3J25 showed the least interaction, as supported by binding free energy calculations in Table 6 for T-muurolol and protein interactions.Figure 6 shows that higher RMSF values for specific holo-protein residues 80-100 for 1CX2-ibuprofen, 390-410 for 1MWU-T-muurolol, 75-78 for 3N8D-2-Oxazolidinone, 175-240 for 6KSI-2-Oxazolidinone, 400-550 for 3J25-2-Oxazolidinone, and 550-620 for 3J25-T-muurolol indicate that proteins associated with T-muurolol maintained greater residual rigidity compared to their control drug complexes.This suggests that T-muurolol exhibited more stable interactions with all proteins over the 100 ns MD timespan compared to the control drugs.
showed the smallest change compared to their apo-proteins.This indicates that 3J25 maintained residual flexibility, whereas 6KSI retained rigidity in the presence of specific ligands.Among the six holo-proteins, 6KSI demonstrated the best interactions with ligands and polar molecules, while 3J25 showed the least interaction, as supported by binding free energy calculations in Table 6 for T-muurolol and protein interactions.Figure 6 shows that higher RMSF values for specific holo-protein residues 80-100 for 1CX2-ibuprofen, 390-410 for 1MWU-T-muurolol, 75-78 for 3N8D-2-Oxazolidinone, 175-240 for 6KSI-2-Oxazolidinone, 400-550 for 3J25-2-Oxazolidinone, and 550-620 for 3J25-T-muurolol indicate that proteins associated with T-muurolol maintained greater residual rigidity compared to their control drug complexes.This suggests that T-muurolol exhibited more stable interactions with all proteins over the 100 ns MD timespan compared to the control drugs.Structural compactness is analyzed by R g (radius of gyration) values.A steady and less fluctuating R g value indicates a more compact protein structure.R g , the distance between the centre of mass of all protein atoms, is crucial for assessing stable protein folding [23].As shown in Figure 7, owing to the volumetric size of proteins, the R g values for the apo-protein and holo-protein systems varied as follows: 2.44-2.54nm for 1CX2, 3.55-3.80nm for 1MWU, 2.10-2.40nm for 3N8D, 1.54-1.65 nm for 3SRW, 2.05-2.18nm for 6KSI, and 2.80-3.70nm for 3J25.The comparative R g trajectory of apo-protein and holo-protein systems in Figure 7A, B, D shows that 1CX2, 1MWU, and 3SRW systems consistently maintained steady R g values with moderate fluctuations.The holo-proteins closely overlapped with their corresponding apo-proteins, indicating that holo-proteins retained compactness similar to apo-proteins despite ligand interactions.Conversely, Figure 7C,E show that 3N8D-T-muurolol and 6KSI-T-muurolol exhibited greater compactness and integrity compared to their apo-protein and protein-control drug counterparts.The R g trajectories for 3N8D-T-muurolol and 6KSI-T-muurolol were steadier and declined gradually after 40 and 60 ns, respectively, over the 100 ns MD run.Additionally, in Figure 7F, the 3J25 apo-protein showed steady behaviour with moderate R g fluctuations compared to its holo-proteins.This indicates that 3J25 holo-proteins progressively lost compactness and stability in protein folding after successive ligand interactions.
fluctuating Rg value indicates a more compact protein structure.Rg, the distance between the centre of mass of all protein atoms, is crucial for assessing stable protein folding [23].As shown in Figure 7, owing to the volumetric size of proteins, the Rg values for the apo-protein and holo-protein systems varied as follows: 2.44-2.54nm for 1CX2, 3.55-3.80nm for 1MWU, 2.10-2.40nm for 3N8D, 1.54-1.65 nm for 3SRW, 2.05-2.18nm for 6KSI, and 2.80-3.70nm for 3J25.The comparative Rg trajectory of apo-protein and holo-protein systems in Figure 7A, B, D shows that 1CX2, 1MWU, and 3SRW systems consistently maintained steady Rg values with moderate fluctuations.The holo-proteins closely overlapped with their corresponding apo-proteins, indicating that holo-proteins retained compactness similar to apo-proteins despite ligand interactions.Conversely, Figure 7C,E show that 3N8D-T-muurolol and 6KSI-Tmuurolol exhibited greater compactness and integrity compared to their apo-protein and protein-control drug counterparts.The Rg trajectories for 3N8D-T-muurolol and 6KSI-Tmuurolol were steadier and declined gradually after 40 and 60 ns, respectively, over the 100 ns MD run.Additionally, in Figure 7F, the 3J25 apo-protein showed steady behaviour with moderate Rg fluctuations compared to its holo-proteins.This indicates that 3J25 holo-proteins progressively lost compactness and stability in protein folding after successive ligand interactions.

SASA Landscape Analysis
SASA measures the protein surface area accessible to solvent molecules.With TIP3P water, hydrophilic apo-proteins showed increased SASA, while hydrophobic ones showed decreased SASA due to unfolding and folding.Ligand binding can shield surface residues, altering hydrophilicity and hydrophobicity, thus impacting the SASA and relative size of holo-proteins compared to apo-proteins [24].As illustrated in Figure 8, SASA varied among apo-protein and holo-protein systems according to protein size: ~(240-280) nm 2 for 1CX2, ~(308-337) nm 2 for 1MWU, ~(160-185) nm 2 for 3N8D, ~(90-100) nm 2 for 3SRW, ~(165-195) nm 2 for 6KSI, and ~(335-420) nm 2 for 3J25.In the 1CX2, 3N8D, and 6KSI systems, SASA of T-muurolol-bound holo-proteins significantly decreased during 50-100 ns, 40-100 ns, and 10-80 ns MD intervals, respectively, compared to their apo-proteins and protein-control drug complexes.This suggests substantial T-muurolol binding to surface residues during these intervals, leading to increased hydrophobicity, size shrinkage, and protein folding for 1CX2-T-muurolol, 3N8D-T-muurolol, and 6KSI-T-muurolol.Conversely, in Figure 8F, 3J25-T-muurolol showed increased SASA, indicating higher hydrophilicity than its apo-protein and protein-control drug complex during 10-85 ns.This suggests that ligand interactions led to greater hydrophilicity, protein unfolding, and relative size expansion for 3J25-T-muurolol compared to 3J25 apo-protein and 3J25-2-Oxazolidinone. Furthermore, substantial overlapping SASA trajectories and moderate alterations in timedependent SASA values between apo-protein and holo-protein systems of 1MWU and 3SRW suggest comparatively less significant ligand-protein interactions in these systems.Except for 1MWU holo-proteins, the differing SASA profiles of T-muurolol-protein interaction complexes indicate that hydrophobic and Van Der Waals interactions predominantly contributed to binding affinity in T-muurolol holo-proteins, as supported by the MM-PBSA andscape of T-muurolol-protein interactions in Table 6.
drug complexes.This suggests substantial T-muurolol binding to surface residues during these intervals, leading to increased hydrophobicity, size shrinkage, and protein folding for 1CX2-T-muurolol, 3N8D-T-muurolol, and 6KSI-T-muurolol.Conversely, in Figure 8F, 3J25-T-muurolol showed increased SASA, indicating higher hydrophilicity than its apo-protein and protein-control drug complex during 10-85 ns.This suggests that ligand interactions led to greater hydrophilicity, protein unfolding, and relative size expansion for 3J25-Tmuurolol compared to 3J25 apo-protein and 3J25-2-Oxazolidinone. Furthermore, substantial overlapping SASA trajectories and moderate alterations in time-dependent SASA values between apo-protein and holo-protein systems of 1MWU and 3SRW suggest comparatively less significant ligand-protein interactions in these systems.Except for 1MWU holo-proteins, the differing SASA profiles of T-muurolol-protein interaction complexes indicate that hydrophobic and Van Der Waals interactions predominantly contributed to binding affinity in T-muurolol holo-proteins, as supported by the MM-PBSA andscape of T-muurolol-protein interactions in Table 6.

H Bonds Formation Analysis
GROningen MOlecular Simulation (GROMACS) analysis predicts ligand-protein complex stability by assessing hydrogen bond formation.The occupancy of H bonds in the interaction complexes of protein-ligands aids in heightening electrostatic energy flow [23].In Figure 9

H Bonds Formation Analysis
GROningen MOlecular Simulation (GROMACS) analysis predicts ligand-protein complex stability by assessing hydrogen bond formation.The occupancy of H bonds in the interaction complexes of protein-ligands aids in heightening electrostatic energy flow [23].In Figure 9

Evaluation of MM-PBSA Binding Free Energy Decomposition
Table 6 presents the binding free energy decomposition for protein-ligand interactions involving T-muurolol with six macromolecules (1CX2, 1MWU, 3N8D, 3SRW, 6KSI, and 3J25).The MM-PBSA binding free energies (ΔGBA) are approximately ~(−39.48)KJ/mol (1CX2), ~(−55.58)KJ/mol (1MWU), ~(−21.39)KJ/mol (3N8D), ~(−11.03)KJ/mol (3SRW), ~(−50.64)(6KSI), and ~(−5.01)KJ/mol (3J25).Van der Waals interactions contributed the most to ΔGBA, indicating that short-range dipole-dipole, permanent, and transient hydrophobic interactions are crucial in the protein-ligand binding stability.In addition to van der Waals energy, electrostatic and polar solvation energies reflect hydrogen bond formation, desolvation, and solvation effects among charged dipoles and polar and non-polar molecules in the MD environment.The protein-ligand complexes are surrounded by the polar solvent TIP3P (water) molecules and ions (K ⁺ and Cl − ).Thus, the conformational stability and folding of the proteins are heavily influenced by polar solvation and electrostatic energies.Polar solvation impacts the interactions by shielding the proteins' polar hydrophilic residues, decreasing the hydrophobic binding affinity of ligands.This indicates that an increase in polar solvation energy negatively correlates with ligand-protein interactions, as it reduces the effectiveness of other binding energies (SASA, van der Waals, hydrophobic, and electrostatic).Consequently, higher polar solvation energy negatively impacts the overall stability of strong protein-ligand interactions [25].Strong protein-ligand interactions influenced by polar solvation release positive energy, while SASA and electrostatic energies favour negative energy flow, indicating exothermic binding.Polar solvents like TIP3P enhance SASA interactions for hydrophilic surfaces, promoting protein folding and increasing hydrophobic interactions.Lower SASA in protein-ligand complexes correlates with better binding affinity, strengthening hydrophobic interactions and stabilizing the complex [24,26].As shown in Table 6, the 1MWU-T-muurolol complex had the highest ΔGBA due to increased polar solvation energy and decreased electrostatic, van der Waals, and hydrophobic energies.This indicates significant polar solvation with TIP3P molecules interacting with surface hydrophilic residues, reducing SASA and other binding energies.Conversely, 6KSI maintained the strongest interaction with a ΔGBA of ~(−50.64)KJ/mol through maximized hydrophobic and electrostatic interactions and minimized polar solvation effects over a 100  Van der Waals interactions contributed the most to ∆G BA , indicating that short-range dipole-dipole, permanent, and transient hydrophobic interactions are crucial in the protein-ligand binding stability.In addition to van der Waals energy, electrostatic and polar solvation energies reflect hydrogen bond formation, desolvation, and solvation effects among charged dipoles and polar and non-polar molecules in the MD environment.The protein-ligand complexes are surrounded by the polar solvent TIP3P (water) molecules and ions (K + and Cl − ).Thus, the conformational stability and folding of the proteins are heavily influenced by polar solvation and electrostatic energies.Polar solvation impacts the interactions by shielding the proteins' polar hydrophilic residues, decreasing the hydrophobic binding affinity of ligands.This indicates that an increase in polar solvation energy negatively correlates with ligand-protein interactions, as it reduces the effectiveness of other binding energies (SASA, van der Waals, hydrophobic, and electrostatic).Consequently, higher polar solvation energy negatively impacts the overall stability of strong protein-ligand interactions [25].Strong protein-ligand interactions influenced by polar solvation release positive energy, while SASA and electrostatic energies favour negative energy flow, indicating exothermic binding.Polar solvents like TIP3P enhance SASA interactions for hydrophilic surfaces, promoting protein folding and increasing hydrophobic interactions.Lower SASA in protein-ligand complexes correlates with better binding affinity, strengthening hydrophobic interactions and stabilizing the complex [24,26].As shown in Table 6, the 1MWU-T-muurolol complex had the highest ∆G BA due to increased polar solvation energy and decreased electrostatic, van der Waals, and hydrophobic energies.This indicates significant polar solvation with TIP3P molecules interacting with surface hydrophilic residues, reducing SASA and other binding energies.Conversely, 6KSI maintained the strongest interaction with a ∆G BA of ~(−50.64)KJ/mol through maximized hydrophobic and electrostatic interactions and minimized polar solvation effects over a 100 ns MD interval.This suggests that higher polar solvation energy and ∆G BA are not definitive indicators of strong protein-ligand interactions [25].Thus, it is suggested that the combined contributions of van der Waals, electrostatic, SASA, and hydrophobic interactions in enhancing ∆G BA support robust ligand-protein interactions.

Evaluation of MM-PBSA Binding Free Energy Decomposition
2.6.Quantum Chemical Computational Analysis 2.6.1.Frontier Molecular Orbital (FMO) Analysis Molecular orbital (MO) analysis, especially the examination of frontier molecular orbitals (FMOs) such as the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), is essential for understanding a molecule's electronic structure and reactivity.FMOs offer valuable insights into molecular interactions and help predict the chemical reactivity and stability of compounds.Figure 10A displays the HOMO and LUMO molecular orbitals for the optimized geometry of T-muurolol.The HOMO and LUMO energies for T-muurolol are −6.37 eV and 0.52 eV, respectively, with an energy gap (∆E Gap ) of 6.9 eV.In the FMO plots, the green and red lobes represent the orbital phases, with one colour indicating a positive phase and the other a negative phase.Hyperconjugation stabilizes molecular orbitals through overlap between an occupied orbital and a neighbouring electron-deficient orbital.As shown in Figure 10A, the HOMO is localized over the double bond region, typical of high electron density in alkenes, with some electron density around the hydroxyl group, indicating potential nucleophilic attack sites.The LUMO, on the other hand, shows electron density around the carbon atoms adjacent to the double bond and the hydroxyl group.The LUMO lobes, like those of the HOMO, also display phases (green and red) but are located in different regions of the molecule.The significant electron density around the double bond in both orbitals highlights this region as crucial for the molecule's reactivity.
Int. J. Mol.Sci.2024, 25, x FOR PEER REVIEW 16 of 24 ns MD interval.This suggests that higher polar solvation energy and ΔGBA are not definitive indicators of strong protein-ligand interactions [25].Thus, it is suggested that the combined contributions of van der Waals, electrostatic, SASA, and hydrophobic interactions in enhancing ΔGBA support robust ligand-protein interactions.

Frontier Molecular Orbital (FMO) Analysis
Molecular orbital (MO) analysis, especially the examination of frontier molecular orbitals (FMOs) such as the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), is essential for understanding a molecule's electronic structure and reactivity.FMOs offer valuable insights into molecular interactions and help predict the chemical reactivity and stability of compounds.Figure 10A displays the HOMO and LUMO molecular orbitals for the optimized geometry of T-muurolol.The HOMO and LUMO energies for T-muurolol are −6.37 eV and 0.52 eV, respectively, with an energy gap (ΔEGap) of 6.9 eV.In the FMO plots, the green and red lobes represent the orbital phases, with one colour indicating a positive phase and the other a negative phase.Hyperconjugation stabilizes molecular orbitals through overlap between an occupied orbital and a neighbouring electron-deficient orbital.As shown in Figure 10A, the HOMO is localized over the double bond region, typical of high electron density in alkenes, with some electron density around the hydroxyl group, indicating potential nucleophilic attack sites.The LUMO, on the other hand, shows electron density around the carbon atoms adjacent to the double bond and the hydroxyl group.The LUMO lobes, like those of the HOMO, also display phases (green and red) but are located in different regions of the molecule.The significant electron density around the double bond in both orbitals highlights this region as crucial for the molecule's reactivity.A Molecular electrostatic potential (MEP) map visually represents the charge distribution across a molecule, highlighting positive and negative potential areas.This is useful for predicting molecular interactions, particularly electrostatic interactions.The MEP map of T-muurolol (Figure 10B) emphasizes the hydroxyl group as a key functional site due to its distinct polar regions.The map shows regions of positive potential (blue, 4.528 × 10 2 −e) that are likely to attract nucleophiles and areas of negative potential (red, −4.528 × 10 −2 e) that are likely to attract electrophiles, with the rest of the molecule being relatively neutral (green).The pronounced blue and red regions around the hydroxyl group indicate significant polarity, making this group in T-muurolol a reactive site for hydrogen bonding and polar interactions.The neutral green regions, covering most of the hydrocarbon skeleton, suggest that these parts of the molecule are less reactive but contribute to its overall shape and steric properties (Figure 10B).

Global Chemical Reactivity Parameters
In computational chemistry, Koopmans' theorem [27] serves as a powerful tool for estimating key electronic properties of molecules, enabling the derivation of essential global reactivity parameters.By applying this theorem, we can derive several crucial parameters that provide quantitative insights into the behaviour of the T-muurolol molecule, including electronegativity (χ = −1/2 (E HOMO + E LUMO )), chemical potential (µ = −χ), global hardness index (η = 1/2 (E LUMO − E HOMO )), global softness index (σ = 1/η), and global electrophilicity index (ω = µ 2 /2η).All parameters listed above (Table 7) provide a quantitative framework for predicting and rationalizing the chemical behaviour of T-muurolol molecules, aiding in the design and synthesis of new compounds with desired reactivity profiles.

Discussion
T-muurolol is a bioactive sesquiterpene predominantly found as a major compound in the essential oils of several plants, including Alpinia zerumbet, Pulicaria somalensis, Calocedrus formosana, Calendula officinalis, and Chamaecyparis obtusa with reported antifungal, antibacterial, and antioxidant activities [13][14][15]28,29].Sesquiterpenes, due to their lipophilicity, integrate into bacterial lipid membranes and passively diffuse through the bilayer, disrupting membrane integrity.This increases permeability and destabilizes membrane proteins, compromising bacterial function.In Gram-positive bacteria, such as S. aureus, the porous peptidoglycan layer allows sesquiterpenes like T-muurolol to diffuse through and penetrate the cytoplasmic membrane, where they inhibit specific bacterial targets [30].With a documented history of bioactivity, T-muurolol could be an effective agent in inhibiting S. aureus-associated infections.For that, it was screened against various drug targets of S. aureus to evaluate its efficacy as an alternative antibacterial agent.
T-muurolol exhibited a strong binding interaction with S. aureus lipase (SAL) which is a triacylglycerol esterase, a key virulence factor that degrades immune-responsive lipids, inhibits innate immune cell activation, and disrupts host immune recognition [31].So, SAL is a major drug target for suppressing S. aureus pathogenesis.Strong interactions with SAL indicate its potency as an inhibitor.However, few SAL inhibitors have been introduced until now, suggesting that T-muurolol could be an effective agent for targeting SAL in further drug development.Similarly, dihydrofolate reductase, a key protein in the thymidine synthesis pathway essential for DNA synthesis, is crucial for bacterial survival [32].The strong binding interaction of T-muurolol with this enzyme suggests its efficacy in inhibiting bacterial replication.
MDR S. aureus strains are major contributors to S. aureus-associated infections, posing significant management challenges and life-threatening risks to humans.T-muurolol has demonstrated strong molecular interactions with various drug targets specific to MDR S. aureus strains.For MRSA, T-muurolol targets penicillin-binding protein 2a, thereby inhibiting cell wall biosynthesis.It also exhibits a high binding affinity for D-Ala:D-Ala ligase, a critical drug target in VRSA.Additionally, T-muurolol interacts strongly with the RPP TetM protein in complex with the 70S ribosome, which is implicated in tetracycline resistance in TetRSA by mutating the interaction point of domain IV with the tetracycline binding site [33].These interactions suggest that T-muurolol can inhibit TetRSA susceptibility by effectively targeting this protein.Furthermore, T-muurolol exhibits effective binding interactions with FAD-dependent NAD(P)H oxidase, which is involved in the ROS mechanism in the human body [34], and cyclooxygenase-2, which catalyzes the biosynthesis of pro-inflammatory prostanoids.These interactions establish its antioxidant and antiinflammatory activities.These results consequently indicate the potential of T-muurolol as an effective agent for mitigating S. aureus-associated inflammation and cellular oxidative stress.Moreover, MD simulation analysis elucidated the stability of T-muurolol interactions with selected proteins.RMSD analysis revealed enhanced backbone stability across all T-muurolol interactions, compared to protein-control drug interactions.Similarly, RMSF analysis indicated stable interaction of T-muurolol with all proteins over the 100 ns MD timespan, surpassing control drugs.Except for the penicillin-binding protein 2a complex, differing SASA profiles of T-muurolol-protein interaction complexes suggested that hydrophobic and van der Waals interactions predominantly contributed to T-muurolol's binding affinity.R g values of proteins and protein-ligand complexes highlighted increased compactness and structural integrity in the presence of T-muurolol relative to control drugs.Hydrogen bonds remained consistent with docking scores, exhibiting minimal alterations in protein-ligand interaction dynamics.Our simulation study, consistent with the docking results, reveals that the 1MWU and 3SRW complexes exhibit less stable interactions compared to other T-muurolol-protein complexes, with a strong binding affinity driven primarily by hydrophobic and electrostatic interactions, rather than hydrogen bonding.Furthermore, the recalculated binding energies of selected drugs using the MM-PBSA method after a 100 ns MD simulation demonstrated T-muurolol's superior binding energy toward targeted proteins.
Additionally, quantum chemical structure analysis revealed a ∆E Gap of 6.9 eV for T-muurolol, indicating greater stability and lower polarizability, resulting in higher chemical hardness [35].A higher hardness value suggests lower reactivity [20], indicating the less reactive nature of T-muurolol.However, molecules with smaller HOMO-LUMO energy gaps are more interactive and inclined to participate in chemical events involving bond breaking or formation [36].T-muurolol, with higher HOMO-LUMO values, signifies less anticipation in bond formations like covalent, hydrogen, and ionic bonds, correlating with molecular interaction studies showing fewer bond formations during protein-ligand interactions.A lower global electrophilicity index suggests a strong nucleophilic molecule, while a higher value indicates a strong electrophilic characteristic [35,37].T-muurolol exhibits a moderate electrophilic nature with a value of 1.24 eV, indicating moderate facilitation of bond formation with biomolecules.Chemical potential reflects the electron escaping tendency from the molecule, with a lower value implying greater stability, indicating the greater structural stability of T-muurolol.Additionally, MEP analysis highlights the hydroxyl group in T-muurolol as the functional active site for hydrogen bonding and polar interactions.
Based on ADMET pharmacokinetics analysis, T-muurolol exhibits lower water solubility but possesses good absorption and distribution properties within the human body.Toxicity parameter analysis indicates that T-muurolol can be a promising drug candidate for inhibiting S. aureus bacterial-specific proteins, with very low toxicity observed.These findings suggest T-muurolol as an antibacterial agent against S. aureus, which has the benefit of reducing infection-related inflammation and oxidative stress.Several studies have demonstrated the antibacterial potential of T-muurolol, as a key bioactive compound in essential oils.For instance, Su et al. [38] identified T-muurolol as a major bioactive compound in Diospyros discolor essential oil, exhibited strong antibacterial activity against pathogens, including S. aureus, with MICs ranging from 31.25 to 500 µg/mL.In Fraxinus angustifolia essential oil, T-muurolol was found as an abundant compound and effective against both Gram-positive and Gram-negative bacteria, with MICs between 20 and 300 µg/mL [39].Additionally, Hirtellina lobelia essential oil, rich in T-muurolol, showed activity against MDR bacteria and fungi, with MICs ranging from 4 to 512 µg/mL [40].These findings underscore T-muurolol as a significant contributor to the antibacterial efficacy of these essential oils, highlighting its potential as a potent antibacterial agent.
Phytochemical-based drug discovery using molecular docking and ADMET analysis provides a cost-effective method for identifying novel compounds and addressing multidrug resistance.However, in silico predictions are limited by structural data and may not fully capture accurate, effective drug concentration, pharmacokinetics, or toxicity.Future research should incorporate deep learning and advanced computational models to improve accuracy.Additionally, further in vitro and in vivo studies will be performed to validate T-muurolol's efficacy against S. aureus and its MDR strains and compare them with in silico findings.This study lays the groundwork for subsequent investigations, and the follow-up studies can provide insight into diverse prevention strategies for S. aureus infections.
The crystal structures of all proteins were retrieved from the RCSB PDB (accessed on March 12, 2024) and were prepared for docking studies by removing unwanted chains and residues using UCSF Chimera 1.16.
For ligand preparation, the 3D SDF structures of T-muurolol, ascorbic acid, ibuprofen, and 2-oxazolidinone were downloaded from PubChem (accessed on March 12, 2024).In this study, 2-oxazolidinone was selected as the reference compound due to its role as the fundamental structural backbone of the oxazolidinone class of antibiotics, which is primarily responsible for the antibacterial efficacy of all derivatives within this group and effective against all MDR S. aureus.These compounds were imported into Avogadro 1.2.0 and subjected to energy minimization using the MMFF94 force field with the steepest descent algorithm.The structures were then converted to PDBQT format for docking.

Molecular Docking and Binding Analysis
AutoDock Tools 1.5.7 (Scripps Research Institute, La Jolla, CA, USA) was used to predict protein-ligand interactions.Water molecules were removed, and polar hydrogens and Kollman charges were added to the protein structures.The protein and ligand files were converted to PDBQT format.Active sites were selected based on blind-docking and enclosed within a 3D affinity grid centred on the active site residues.Docking was performed via command prompt, as described by Gupta et al. [42].Binding energies were recorded, and initial visualizations were performed with BIOVIA Discovery Studio Visualizer (BIOVIA, San Diego, CA, USA).The follow-up detailed analysis of amino acid and ligand interactions was also performed using the same.Compounds with the best binding affinities were selected for further molecular dynamics simulation.
Furthermore, for cross-checking the binding affinities of the ligands with proteins, docking was performed in the CB-Dock2 online server (https://cadd.labshare.cn/cb-dock2/index.php,accessed on 12 March 2024).

Pharmacokinetics Prediction and Bio-Activity Analysis
Pharmacokinetics parameters related to absorption, distribution, metabolism, excretion, and toxicity (ADMET) play a substantial role in the detection of novel drug candidates.To predict candidate molecules using in silico methods, the pkCSM (https: //biosig.lab.uq.edu.au/pkcsm/prediction,accessed on 12 March 2024) online tool was used.Parameters such as AMES toxicity, maximum tolerated dose (human), hERG I and hERG II inhibitory effects, oral rat acute toxicities, hepatotoxicity, skin sensitization, and fathead minnow toxicity were explored.In addition to these, molecular weight, hydrogen bond acceptor, hydrogen bond donor, number of rotatable bonds, topological polar surface area, octanol/water partition coefficient, and number of violations of Lipinski's rule of five were also surveyed using SwissADME (http://www.swissadme.ch/,accessed on 18 March 2024).The molecular properties and bioactivity score of the ligands were obtained from the molinspiration web server (https://www.molinspiration.com/cgi/properties,accessed on 18 March 2024).

Molecular Dynamics and Simulation
The six best protein-ligand complexes from the molecular docking study, selected based on the lowest binding energy and optimal docked pose, were chosen for MD simulation.The macromolecular structures (PDB ID_1CX2, 1MWU, 6KSI, and 3N8D) consisted of multiple subunits that were segregated; their additional homodomains had been removed via Discovery Studio Visualizer v20.1.0.19195 (https://www.3ds.com/products/biovia/discovery-studio, accessed on 18 April 2024), followed by the removal of adjacent heteroatoms via similar fashion.Then, the modified structures were further modified to become compatible with in silico virtual screening and molecular dynamic simulation (MDS).A comparison map of the dynamic characteristics of targeted proteins and their protein-ligand complexes using GROMACS 2020.1 according to the described methods by Kandasamy et al. [40].In this study, the PDB structures of the proteins were transformed to gmx (gromacs) format using the CHARMM36 force field (http://www.charmm-gui.org/,accessed on 18 April 2024).The parameter files for the docked complexes were created using the guidelines outlined in the GROMACS course.The topology and parameter files for protein and ligand were generated on the CHARMM-GUI server, and selected protein-ligand complex files were prepared using the CHARMM-GUI Multicomponent Assembler (MCA).The cube enclosing the system was sized based on 'Calculated Volume Fraction', 'Minimum Recommended Size Length', and 'Maximum Volume Fraction', with a 20 Å buffer to prevent periodic image overlap during MD simulation.Following the final system size determination, the solvent composition was calculated, and the system was solvated with TIP3P water and 150 mM K + and Cl − ions per the Solution Builder protocol.GROMACS files for proteins and protein-ligand complexes were extracted from CHRMM-GUI files for 100 ns (nanoseconds) with 2 femtoseconds (fs) steps.
The analysis of root means square deviation (RMSD), root means square fluctuations (RMSF) of backbone atoms, solvent accessible surface area (SASA), radius of gyration (R g ), and hydrogen bond (HB) formation were analyzed after the successful completion of MD simulation of specific protein-ligand interaction complexes with a time interval of 20 ns throughout the whole 100 ns of MD trajectory.All the graphs (RMSD, RMSF, SASA, R g , and HB) were plotted by using Qtgrace (https://qtgrace.sourceforge.io/,accessed on 18 April 2024).

Binding Free Energy Analysis
Apart from in silico docking-based binding affinity analysis between proteins and specific ligands, the strength of protein-ligand interactions was further assessed using MD simulation-specific binding free energy analysis.This was accomplished by calculating the binding free energies with the MM-PBSA (Molecular Mechanics Poisson-Boltzmann Surface Area) method.Using this method, the GROMACS function "g_mmpbsa" captured MD trajectory files of protein-ligand interaction complexes throughout 100 ns at 5 ns intervals by using the following equation: where ∆G BA , ∆E MM , ∆G PBSA , and T∆S MM are denoted as average free energy, average molecular mechanics energy, solvation energy, and solute configuration entropy, respectively [43].

Computational Chemical Structure Analysis
The chemical structure of the T-muurolol molecule was optimized using the DFT/B3LYP method with the 6-311G(d,p) basis set in the Gaussian-09 programme (Gaussian 09, Revision A.02) [44].This optimized geometry was then used to generate the frontier molecular orbitals (FMOs) and molecular electrostatic potential (MEP) plots for the T-muurolol molecule.Additionally, the chemical potential of the structure was examined using DFT.

Conclusions
The present study explored the potential of T-muurolol as an alternative antibacterial agent against S. aureus.Several in silico approaches, such as molecular docking, MD simulation, quantum chemical structure analysis, and pharmacokinetic prediction, have been utilized to screen the efficacy of T-muurolol as a drug.As a result, T-muurolol was identified as a potential antibacterial agent that can inhibit several S. aureus-associated bacterial proteins responsible for pathogenesis, bacterial survival, and multi-drug resistance.Along with it, T-muurolol is also identified as an antioxidant and anti-inflammatory agent that signifies its role in mitigating bacterial infection-associated inflammation and cellular oxidative stress in the human body.Further, quantum chemical structure analysis revealed its structural stability and less reactive nature.The results from the ADMET and bioactivity analysis strongly indicate that T-muurolol has the potential to act as a drug and demonstrate significant biological activities.Moreover, this study suggests that T-muurolol has the potential to counteract S. aureus infections as an alternative antibacterial solution.However, further detailed in vitro and in vivo experimental validation is required to confirm its activity.
, 6KSI-T-muurolol consistently formed one intermolecular hydrogen bond throughout the MD interval ~(26-93) ns, indicating strong bond stability.Conversely, 1MWU-T-muurolol displayed intermittent bond formation at ~85 ns, suggesting less stable interactions.In this case, the H bond is not sufficient to serve as significant evidence of a stable binding interaction between T-muurolol and 1MWU.Additionally, 1CX2-Tmuurolol showed regular and irregular bond formation patterns at specific intervals, while 3N8D-T-muurolol exhibited intermittent bond formation across multiple MD time spans.3SRW-T-muurolol maintained regular bond formation, whereas 3J25-T-muurolol formed bonds infrequently at various MD intervals.muurololmaintained regular bond formation, whereas 3J25-T-muurolol formed bonds infrequently at various MD intervals.

Figure 9 .
Figure 9. Illustration of the hydrogen bonds' intensity up to MD time span of 100,000 ps (picoseconds) formed between T-muurolol and respective proteins depicted in the figure itself.

Figure 9 .
Figure 9. Illustration of the hydrogen bonds' intensity up to MD time span of 100,000 ps (picoseconds) formed between T-muurolol and respective proteins depicted in the figure itself.

Figure 10 .
Figure 10.(A) HOMO and LUMO molecular orbitals of the T-muurolol molecule, along with their energy gap reported in electron volts (eV) using an iso-surface value of 0.02 electron.(B) Molecular

Figure 10 .
Figure 10.(A) HOMO and LUMO molecular orbitals of the T-muurolol molecule, along with their energy gap reported in electron volts (eV) using an iso-surface value of 0.02 electron.(B) Molecular electrostatic potential map of T-muurolol molecule (red atom: O; blue atom: N; grey atom: C; white atom: H).These plots are derived from optimizing the T-muurolol molecule using the DFT/B3LYP/6 −311G(d,p) basis set level.Molecular electrostatic potential analysis.

Table 2 .
Top hit interactions of T-muurolol with proteins and their hydrogen bonds and several other interactions.

Table 5 .
Prediction of bioactivity score of compounds.

Table 5 .
Prediction of bioactivity score of compounds.

Table 6 .
Binding free energy (KJ/mol) calculation of the following proteins based on MD interaction with T-muurolol.

Table 6 .
Binding free energy (KJ/mol) calculation of the following proteins based on MD interaction with T-muurolol.