Bacterial Quorum-Sensing Molecules as Promising Natural Inhibitors of Candida albicans Virulence Dimorphism: An In Silico and In Vitro Study

Farnesol, a self-secreted quorum-sensing molecule (QSM) of Candida albicans, has been known to limit yeast-to-hyphal transition by blocking the RAS1–cAMP–PKA pathway. In a similar fashion, certain bacterial QSMs have also been reported to be successful in attenuating C. albicans biofilm and hyphal formation at relatively high cell density. This prompted us to investigate the antihyphal efficacy of certain bacterial QSMs through virtual docking against seminal drug targets, viz., CYCc and RAS1, that have been reported to be the hallmark players in C. albicans dimorphic virulence cascade. Against this backdrop, 64 QSMs belonging to five different bacterial QS signaling systems were subjected to initial virtual screening with farnesol as reference. Data of the virtual screening unveiled QSMs belonging to diketopiperazines (DKPs), i.e., 3-benzyl-6-isobutylidene-2,5-piperazinedione (QSSM 1157) and cyclo(l-Pro-l-Leu) (QSSM 1112), as potential inhibitors of CYCc and RAS1 with binding energies of −8.2 and −7.3 kcal mol−1, respectively. Further, the molecular dynamics simulations (for 50 ns) of CYCc-QSSM 1157 and RAS1-QSSM 1112 complexes revealed the mean ligand root mean square deviation (RMSD) values of 0.35 and 0.27 Å, respectively, which endorsed the rigid nature, less fluctuation in binding stiffness, and conformation of binding complexes. Furthermore, the identified two QSMs were found to be good in solubility, absorption, and permeation and less toxic in nature, as revealed by pharmacokinetics and toxicity analyses. In addition, the in vitro antihyphal assays using liquid and solid media, germ-tube experiment, and microscopic analysis strongly validated DKP-QSSM 1112 as a promising inhibitor of hyphal transition. Taken together, the present study unequivocally proves that DKPs can be used as potent inhibitors of C. albicans virulence dimorphism.


INTRODUCTION
Despite the advances in modern medicine, the management of infectious diseases has become more challenging, as microbial pathogens consistently break down every antimicrobial wall through a phenomenon called "antimicrobial resistance" (Llor and Bjerrum, 2014;Dadgostar, 2019). At times, the situation becomes worse with fungus-associated infections owing to the limitation of therapeutic options (Roemer and Krysan, 2014). Among the various fungal species, Candida albicans is the most representative pathogen at clinical setup, which causes severe contagious infections in humans (Jha and Kumar, 2018). C. albicans is a diploid polymorphic fungus that asymptomatically colonizes various niches of healthy humans as commensal. However, given the opportunity, C. albicans turns into a pathogen and causes clinically diverse infections, especially in individuals with disturbed immune surveillance or other debilitating conditions (Garcia et al., 2021). Various endogenous and exogenous signals facilitate C. albicans to manifest its devastating pathogenesis. Nevertheless, a plethora of studies on C. albicans pathogenesis at the molecular level have reinforced the morphological transition from yeast to hyphae as the hallmark event, which in turn triggers all other virulence traits (at any level of pathogenic circumstances) in a cascade manner (Calderone and Fonzi, 2001;Banerjee et al., 2019). Therefore, identifying molecules with potency to target dimorphic switching has been considered as one of the promising alternatives to effectively combat the infections associated with antifungal-resistant C. albicans by nullifying the phenomenon of selection pressure.
To date, numerous natural/synthetic compounds have been reported for their remarkable antihyphal efficacy against C. albicans (Manoharan et al., 2017;Priya and Pandian, 2020). However, farnesol, being a self-produced quorum-sensing molecule (QSM) of C. albicans, creates a high level of curiosity due to its unique mode of action. The presence of QSM in C. albicans was first reported by Hornby et al. (Hornby et al., 2001), and they made it clear that besides external signal, internally secreted farnesol also modulates the regulation of yeast-tohyphal transition and subsequent biofilm formation in accordance with cell density. After the discovery of farnesol, few other research groups have tried to find out its accurate mode of action in inhibiting yeast-to-hyphal transitions (Nickerson and Atkin, 2017;Chen et al., 2018), but a clear understanding is still lacking. It has been well studied and reported that farnesol negatively regulates the RAS1-cAMP-PKA pathway by hindering the functions of CYR1 by binding with one of its active domains called CYCc (Davis-Hanna et al., 2008;Hall et al., 2011). The GTP-bound RAS1 (activated) directly interacts with CYR1 to catalyze the synthesis of cAMP that eventually activates PKA leading to the expression of several virulence genes in response to different environmental cues (Wang, 2013). Thus, the RAS1-cAMP-PKA pathway seems to be responsible for the expression of several genes associated with yeast-to-hyphal transition as well as other virulence traits and establishes infection (Wang, 2013). Furthermore, farnesol produced by in situ planktonic cells also limits the C. albicans biofilm formation owing to the strong association between hyphal morphology and biofilm development (Ramage et al., 2002). Apart from the fungal world, farnesol has also been shown to thwart the virulence such as biofilm and lipase productions in bacteria, viz., Staphylococcus aureus, and sensitizes the drug-resistant S. aureus to gentamicin (Jabra-Rizk et al., 2006).
Given the prominence of farnesol in effectively interfering with the yeast-to-hyphal switch, several natural and synthetic farnesol analogs have been investigated for their improved inhibitory efficacy against C. albicans dimorphic switching (Shchepin et al., 2003;Lee et al., 2020). Following the same paradigm, certain QSMs of bacteria have been envisaged and well demonstrated to greatly influence the morphogenesis of C. albicans. For instance, the QSMs N-(3-oxododecanoyl)-Lhomoserine lactone (3OC12HSL), cis-2-dodecenoic acid (BDSF), and trans-2-decenoic acid (SDSF) produced by Pseudomonas aeruginosa, Burkholderia cenocepacia, and Streptococcus mutans, respectively, have been reported to substantially limit C. albicans hyphae without altering its basic metabolism (Grainha et al., 2020). Besides, it has also been deciphered that 3OC12HSL interferes with the C. albicans yeast-to-hyphal transition by mimicking the function of farnesol. The other two QSMs, namely, BDSF and SDSF, have been reported to follow a distinct mode of action unlike 3OC12HSL (Hogan et al., 2004;Grainha et al., 2020). Naturally, such antagonistic interactions between bacteria and fungi occur in human microbiota wherein the bacterial QSMs play an indispensable role in controlling fungal growth to maintain the ecological balance. Therefore, understanding the mighty role of bacterial QSMs, it is envisaged that these bacterial QSMs could be an effective agent to probe against fungal filamentation by mimicking the action mechanism of the farnesol molecule.
Against this backdrop, the present study was focused to screen and envisage bacterial QSMs with potency to inhibit the seminal protein targets of C. albicans filamentation through in silico and in vitro analyses. Owing to the essentiality of CYR1 and RAS1 in farnesol-mediated hyphal inhibition, they were chosen as potential drug targets in the current study (Wang, 2013;Boutet et al., 2016). Overall, the current study is the first of its kind that provides potential insights into bacterial QSMs mediated inter-microbial cross-talk. Further, the outcome of the study will prompt the researchers working in the arena of alternative medicine to refocus/revisit the bacterial QSMs as a novel hyphal inhibitor against pathogenic fungi.

Homology Modeling
The experimentally solved crystal structure of CYCc and RAS1 was not publicly available in the Protein Data Bank (PDB). Therefore, three-dimensional (3D) structures of these proteins were modeled on the basis of homology modeling using SWISS-Model (https://swissmodel.expasy.org/). Initially, the amino acid sequences of these target proteins, i.e., CYCc (UniProt ID. P0CY32-1) and RAS1 (UniProt ID. A0A1D8PR83-1; aa length. 1,329-1,466), were retrieved from UniProt KnowledgeBase (UniProtKB) in FASTA format. The resulting sequence was submitted to PDB-BLAST in order to identify the appropriate templates for structural modeling. Finally, the 3D structures of CYCs and RAS1 were generated using SWISS-Model (Boutet et al., 2016).

Validation of Modeled Protein
The stereochemical qualities and reliability of modeled structures were verified by PROCHECK (Program to check the stereochemical quality of protein structures), a program that relies on the Ramachandran plot for structure validation. The number of factors such as overall G-factor, the total number of amino acid residues in the core, and allowed, generously allowed, and disallowed regions were considered for choosing the best model (Laskowski et al., 2005).

Protein Preparation
Further, the validated modeled structures of CYCs and RAS1 were prepared for docking using AutoDock. In order to minimize the energy, the solvated structures were refined through the computation of Gasteiger charges with the addition of the polar hydrogens. Finally, refined protein structures were subjected to further analysis.

Ligand Selection and Preparation
A total of 64 bacterial QSMs screened in this study were retrieved from the SigMol database created by Rajput et al. (2016) (http:// bioinfo.imtech.res.in/manojk/sigmol) (Rajput et al., 2016). Basic information about the ligands such as molecular weight and 3D structure (.SDF) was retrieved from PubChem and ChemSpider databases, respectively. Furthermore, the AutoDock formats (.pdbqt) were generated through AutoDock Vina in PyRx 0.8 software for the molecular docking simulation (Rolta et al., 2020).

Docking-Based Virtual Screening
In order to identify the potential inhibitors against putative active sites of target proteins CYCc and RAS1, virtual screening was performed using PyRx, 0.8. The minimized proteins and ligand molecules were imported to PyRx, and blind docking was performed (Schuttelkopf and van Aalten, 2004). The full grid map was constructed around the complete target proteins. The resulting docking conformations were ranked based on their binding affinities in comparison with the reference molecule (farnesol). Finally, the best-docked complexes were visualized using Discovery studio visualize and ligplot+.

Molecular Dynamics Simulation
The best interactive ligand molecules with their target proteins were subjected to 50 ns of molecular dynamics (MD) simulation by Gromacs 5.1.4 simulation package (http://www.gromacs.org/) using all-atom optimized potentials for liquid simulations (OPLS-AA) force field (Abdulazeez, 2019). The topological parameter files of proteins and ligand molecules were prepared from GROMOS96 force field and PRODRG2 web server, respectively. Initially, docked complexes were kept in a cubic box surrounded by SPC, a water molecule (water density 1.0) using SPC/E water models, and further, the whole system was neutralized by the addition of Na + ions. The solvated structures were energy minimized using GROMOS54a7 force field and equilibrated by running simulations for 100 ps under constant NVT (300 K temperature) and NPT (1 atm pressure) ensemble. The equilibrated complexes were further extended to run production MD simulation for 50 ns at constant temperature and pressure. The GROMACS package was used to calculate the structural properties such as root mean square deviation (RMSD), root mean square fluctuation (RMSF), hydrogen bond radius of gyration, and solvent accessible surface area. Finally, the resulting values were plotted as graphs using Origin Pro.

Estimation of Pharmacokinetics and Toxicity Profile
The antihyphal inhibitors screened via molecular docking and molecular simulations were checked for various pharmacokinetic parameters using SwissADME online tool (http://www. swissadme.ch) (Rodrigues et al., 2020). The smileys of both ligands were collected from SigMol database and added into SwissADME as an input to analyze the principles of ADMET (absorption, distribution, metabolism, and elimination), pharmacokinetics, drug likeliness, and medical chemistry. In order to predict the toxicity dosage and toxicity class of the ligands, ProTox-II server (https://tox-new.charite.de/protox_II/) was utilized (Banerjee et al., 2018).

Strain and Culture Condition
The test organism C. albicans (ATCC 90028) used in the present study was procured from HiMedia, India. The fungus was maintained in Sabouraud Dextrose Agar (SDA) plates and routinely cultured in yeast extract peptone dextrose (YEPD) broth at 37°C until used for tests. The culture with 0.1 optical density (OD) (1 × 10 6 CFU/ml) was used as the inoculum to perform antihyphal assays. The hyphal assay was performed in spider medium (consisting of mannitol 1%, K 2 HPO 4 0.2%, and nutrient broth 1%) to allow hyphal elongation.

In Vitro Filamentation Assay
In order to substantiate the in silico finding, in vitro antihyphal assay was performed using farnesol as the positive control (Priya and Pandian, 2020). For this experiment, one of the predicted ligands, i.e., QSSM 1112, was employed to assess its inhibitory propensity on C. albicans dimorphic switching. The QSSM 1112 (purity >98% based on high-performance liquid chromatograph) used in the current study was reported in our previous study (Gowrishankar et al., 2014) and deployed for assays by dissolving in methanol at the final concentration of 100 mg/ml. Briefly, C. albicans cells (1 × 10 6 CFU/ml) were dispensed into 24-well microtiter plates (MTPs) containing 1 ml of spider medium supplemented with QSSM 1112 and farnesol at various concentrations (0-1,024 µg/ml) along with 10% of fetal bovine serum (FBS). After incubation at 37°C for 24 h, the pellets were resuspended in 10 µl of phosphate-buffered saline (PBS), and the ratio of yeast to hyphal cells was further confirmed by a light microscope at ×400 magnification (Nikon Eclipse 80i, Japan).

Germ-Tube Inhibition Assay
The efficiency of QSSM 1112 in inhibiting the C. albicans germtube formation was evaluated through germ-tube assay as described by Zhang et al. (2011) (Zhang et al., 2011). Here again, farnesol was used as a positive control for the comparative analysis. In brief, C. albicans cells (1 OD) were used to inoculate into the 200 µl of FBS along with QSSM 1112 at its hyphal inhibitory concentration (HIC), i.e., 64 µg/ml. The tubes were incubated at 37°C for 5 h to allow germ-tube induction. At different time intervals, i.e., 0, 2, 4, and 6 h, the cells were removed aseptically and examined under a light microscope at ×400 magnification in order to assess the germ-tube induction (Nikon Eclipse 80i, Japan).
Additionally, the inhibitory propensity of QSSM 1112 on preformed germ-tube formation was evaluated. Initially, the C. albicans cells were allowed for germ-tube induction following the protocol mentioned above. Then, the QSSM 1112 (at HIC) was manifested on germinated C. albicans cells and incubated at 37°C for 3 h. The cells without the active agent were considered as untreated control. Subsequent to incubation, the morphological changes between QSSM 1112-treated and untreated control cells were observed under a light microscope at ×400 magnification (Nikon Eclipse 80i, Japan).

Fluorescence Microscopic Analysis
To further ascertain the QSSM 1112 mediated hyphal and germtube inhibition, fluorescence microscopic analysis was performed. For hyphal assay, the C. albicans cells were allowed to form hyphal formation on a 1-cm 2 glass slide in 24-well MTPs containing spider medium supplemented with farnesol and QSSM 1112 at their HIC. After incubation, non-adherent planktonic cells were removed, and adhered hyphal cells on the slides were stained with 1% acridine orange for 20 min.
For germ-tube inhibition assay, C. albicans cells were grown in test tubes containing spider medium supplemented with farnesol and QSSM 1112 at their HIC for 24 h. Subsequent to incubation, cell pellets were collected by centrifugation at 8,000 rpm for 10 min. Then, pellets were resuspended in 100 µl of PBS and followed by staining with 20 µl of 1% acridine orange for 20 min. After incubation, the difference in cell morphology was visualized under fluorescence microscopy (Ts2R, ECLIPSE, and Japan).

RESULTS AND DISCUSSION
Clinical complications associated with the conventional antifungals immensely demand a new treatment strategy against the dimorphic fungus C. albicans. One of the proven strategies that propel in the arena of alternative medicine to get rid of drug resistance is "disarming the virulence without disturbing the metabolism" of any fungus (Charro and Mota, 2015). As hyphal transition is a hallmark event that lionizes commensal yeasts as opportunistic pathogens, blocking this morphological switch is imperative to control infections associated with C. albicans (Roemer and Krysan, 2014;Banerjee et al., 2019). So far, many molecules with antihyphal proficiency have been reported from natural and organic sources against C. albicans (Manoharan et al., 2017;Priya and Pandian, 2020). However, the QSMs of bacteria remain unexplored for antihyphal propensity. Few reports have been focused on the competitive antagonistic interactions between bacteria and C. albicans, which signify that bacterial QSMs could play a key role in targeting and blocking yeast-to-hyphal transition (Hogan et al., 2004;Nickerson and Atkin, 2017). This intertaxon chemical communication between bacteria and fungi prompted us to investigate the effect of bacterial QSMs as natural hyphal inhibitors against C. albicans through in silico approach. Therefore, in the current study, 64 QSMs belonging to five different bacterial QS signaling systems such as acylated homoserine lactones (AHLs), diketopiperazines (DKPs), diffusible signal factors (DSFs), autoinducer-2 (AI-2), and 4hydroxy-2-alkylquinolines (HAQs) were selected and subjected to virtual docking against the two identified seminal drug targets (CYCc and RAS1) of C. albicans. It was anticipated that the outcome of the present study would render a basic understanding of the interaction of bacterial QSMs against important hyphal proteins of C. albicans.

Homology Modeling
Primarily, homology modeling of target protein structures was performed. Protein crystal structure is important for understanding the mechanisms of biological systems and protein complexes to manipulate or inhibit protein function (Vyas et al., 2012). Owing to the unavailability of resolved crystal structure for the target protein in PDB, computational-based techniques were used for 3D structure prediction. Homology modeling offers 3D models of protein structures when only sequence data are available (Satyanarayana et al., 2018).
Since the experimentally solved 3D structure of selected target proteins, i.e., CYCc and RAS1, were not available in freely accessible public databases, homology modeling was done using SWISS-MODEL to obtain the good models. The primary sequences of target proteins CYCc (Q9P977_CANAX) and RAS1 (P0CY32_CANAX) were retrieved from UNIPROTKB database. Based on the results of SWISS-MODEL's Critical Assessment of Structure Prediction (CASP), adenylate cyclase from Trypanosoma brucei (PDB: 1fx2) and GTPase KRas from Homo sapiens (PDB: 6god) were identified to have high sequence similarity, i.e., 64.71% and 34.78% for CYCc and RAS1, respectively. In general, a protein sequence with a minimum 30% identity with target protein is adequate to be considered as an appropriate template for accurate homology modeling (Fiser, 2010). Correspondingly, 1fx2 and 6god were selected as templates for building CYCc and RAS1 models. The 3D structures of modeled proteins are shown in Figures 1A, B.
Subsequently, the qualities of generated protein models were validated through the Ramachandran plot constructed using PROCHECK module from SAVES Server (https://saves.mbi. ucla.edu/). As can be seen from Figure 1C, the predicted 3D structure of CYCc has 95.0%, 12.5%, 0.0%, and 0.0% of residues in the most favored, additional allowed, generously allowed, and disallowed regions, respectively. Similarly, the RAS1 structure has 94.8%, 5.2%, 0.0%, and 0.0% of residues in the most favored, additional allowed, generously allowed, and disallowed regions, respectively ( Figure 1D). According to the Ramachandran plot, a model structure upholding more than 90% of residues in the most favored region is deemed to be as accurate as like the experimentally solved crystal structure (2-Å resolution) (Sobolev et al., 2020). Interestingly, both the generated models showed ≥90% residues [i.e., CYCc (95.0%) and RAS1 (94.8)] in the most favored region. Besides, both the models did not cover any residue in the disallowed region, which ascertained that the SWISS-MODEL predicted models as reasonable and reliable for further molecular docking studies.

Protein Preparation
Generally, the 3D structures built through homology modeling used to have unfavorable bond lengths, bond angles, torsion angles, and contacts. Therefore, energy minimization was done using UCSF Chimera to normalize local bond length and angle and to remove close contacts in the geometric chain for the predicted proteins. Initially, the refined proteins were energy minimized using AMBER ff14SB and Gasteiger charges for standard and non-standard amino acid residues, respectively. Finally, the minimized protein structures were saved in ".pdb" extension and used for further molecular docking.
The 3D structures of 64 bacterial QSMs were retrieved from SigMoL database (Rajput et al., 2016). Open Babel platform was used to equalize the ligand format to ".sdf" extension, which was collected from different primary databases integrated with SigMoL database. The reformatted ligands were energy minimized with uff force field using Open Babel wizard on PYRX 0.8. The ligands were saved in ".pdbqt" extension for further molecular docking and protein-ligand simulation analyses.

Structure-Based Virtual Screening
As the virtual screening workflow offers selective filtration of ligands (lead drug) that are capable of interacting with protein macromolecules (drug target) with good binding energy, it has been an imperative tool in drug discovery (Maia et al., 2020). In the current study, virtual screening was done for 64 QSMs belonging to five different bacterial QS signaling systems, and the best hits are listed (Tables 1 and 2) based on their binding energy toward the selected two seminal drug targets, viz., CYCc and RAS1. The obtained binding energies and hydrogen bond interactions of docked ligands with target proteins were summarized in Tables 1 and 2.
In this study, a well-known antihyphal QSM of C. albicans, viz., farnesol, was used as the positive control. It showed a binding affinity of −6.4 and −7.0 kcal mol −1 toward CYCc and RAS1, respectively. Two filtration processes were done for 64 QSMs to select one best hit against each of the two target proteins. Initially, the top three lead molecules were selected. However, on comparing the binding affinity and binding profile with reference molecule farnesol, the best hit was identified.

Quorum-Sensing Molecules' Interaction Mode Against CYCc
As shown in Table 1, the binding affinity of docked QSMs against the target protein CYCc ranged from −4.9 to −8.2 kcal mol −1 . The maximum binding affinity of −8.2 kcal mol −1 was displayed by a QSM named cyclo(L-Phe-L-4-OH-Pro), which belongs to DKP signaling system. The initial selection was based on the binding affinity toward the target protein in which out of 31 AHLs, 15 DKPs, 14 DSFs, 2 HAQs and 2 AI-2, 7 AHLs, 9 DKPs, and 2 HAQs were considered for selecting the top three hits, as they exhibited higher binding affinity than the farnesol. On the other hand, the QSMs belonging to AI-2 and DSFs displayed considerably lower binding energy than the positive control farnesol, and hence, these were excluded for the selection of the top three hits.
First, the ligand QSSM 1157 isolated from bacterium Streptomyces albulus has shown two hydrogen bonds with the residues of Asp6 and Asp76 at 3.02 and 3.12 distances, respectively. It also showed nine hydrophobic bonds with the residues of Phe4, Glu47, Ala50, Met52, Thr74, Ala75, Arg105, Ser135, and Ser138. Notably, S. albulus has already been reported to be a good producer of the fungicide cycloheximide ( Figure 2Aa) (Yin et al., 2014). The antifungal effect of S. albulus synthesized dipeptides; i.e., e-polylysine was demonstrated against C. albicans (Shima et al., 1984). In addition, the antibiofilm activity of a metabolite produced by other Streptomyces species, i.e., Streptomyces chrestomyceticus, has been reported against C. albicans (Srivastava and Dubey, 2016). Taken together, the present docking result with other existing investigations suggests the plausible existence of intercommunication cross-talk between Streptomyces species and C. albicans, which in a deeper sense certainly requires further investigation.
QSSM 1153, also belonging to DKP signaling system and isolated from soil-dwelling Gram-positive bacterium Saccharothrix espanaensis, has the required structure to interact through two hydrogen bonds with Asp76 and Ser138 at 2.96 and 2.73 distances, respectively. Similar to QSSM 1141, it can also form 10 hydrophobic bonds with the amino acid residues of Phe4, Thr 5, Asp6, Lys45, Glu47, Ala50, Met52, Thr74, Arg105, and Ser135. It is noteworthy to mention here that a new class of antibiotics, namely, saccharomicins produced by S. espanaensis, has already been reported to have potent antimicrobial efficacy against certain pathogenic bacteria and yeast ( Figure 2Ac) (Singh et al., 2000).
Based on the farnesol pattern of hydrogen bond formation and non-hydrophobic interactions with the target proteins CYCc and RAS1, the top three QSMs were examined, and one best ligand molecule was selected. Accordingly, QSSM 1157 was selected as the best ligand molecule among the three DKPs. It can establish two hydrogen bonds (at Asp6 and Asp76) and five nonhydrophobic bonds (at Phe4, Glu47, Ala50, Met52, and Thr74) at the same position as like farnesol with a greater binding affinity score. Hence, this ligand molecule was taken forward for MD analysis to understand its interaction with CYCc protein.

Quorum-Sensing Molecules' Interaction Mode Against RAS1
As tabulated (Table 2), the 64 docked QS ligand molecules displayed a range of binding affinity with a maximum of −8.5 kcal mol −1 and a minimum of −5.4 kcal mol −1 against the modeled crystal structure of RAS1 protein. In particular, a QSM named p-coumaroyl-HSL belonging to the AHL system demonstrated the highest binding affinity (−8.5 kcal mol −1 ). In parallel, the reference molecule (farnesol) exhibited strong binding interactions toward RAS1 with a binding energy of −7 kcal mol −1 ( Figure 2B).
For the identification of the top three hits against RAS1, the same selection criterion used for the CYC1 target protein was employed. Accordingly, the QSMs belonging to DKP (9), AHL (18), and HAQ (2) QS signaling systems were selected and considered for the next level of screening, as they displayed better binding energy than the reference molecule. On the other hand, the QSMs of AI-2 and DSFs showcased marginal binding energy than the biomolecules of other QS systems and reference molecules, and hence they were excluded from further consideration.    Figure 2Bd). The positive control (farnesol) interacts with the RAS1 protein by establishing four hydrogen bonds with the amino acid residues of Asp120, Ser147, Ala148, and Lys149 at the distances of 2.91, 2.86, 2.97, and 2.89, respectively.
The QS molecule QSSM 1151 has been identified from various bacterial sources, namely, Pseudomonas putida, Burkholderia cepacia, Shewanella baltica, and Streptomyces spp. It stabilized the interaction with RAS1 through three hydrogen bonds with Ala148, Lys149, and Asp120 (at binding distances of 3.10, 3.10, and 3.05, respectively) and eight hydrophobic bonds with amino acid residues of Gly16, Ala19, Phe29, Asp31, Asn117, Lys118, Leu121, and Ser147 ( Figure 2Bb). A study conducted by Li et al. (2007) revealed the antifungal efficacy of a novel compound from B. cepacia. Further, B. cepacia has been reported as a source of production of antibiotics such as pyrrolnitrin, altericidins, and cepacin. Furthermore, another study signified the antifungal efficacy of extract from Streptomyces spp. against C. albicans (Human et al., 2016). Notably, the antagonistic mechanism of B. cepacia through its signaling molecules, i.e., BDSF, toward C. albicans has already been well investigated by Boon et al. (2008), since the QSM of B. cepacia showed a better binding interaction with RAS1, which further signifies the possible interspecies communication-mediated regulation of C. albicans virulence dimorphism by B. cepacia QSM.
Out of the top three hits, the best one QSM against RAS1 was selected based on the same interaction criteria followed for CYCc. As a result, the QS ligand molecule QSSM 1112 was selected as the best ligand, as it established four strong hydrogen Q18 bonds (at Asp120, Ser147, Ala148 and Lys 149) and six nonhydrophilic bonds (Gly14, Gly16, Ala19, Phe29, Tyr33 and Asn117 ) with RAS1 and taken forward to investigate the QSSM 1112-RAS1 complex stability.
Overall, out of 64 QSMs belonging to five different QS signaling systems used in the virtual docking analysis, two QSMs, viz., QSSM 1157 and QSSM 1112, belonging to DKPs were selected as the best hits against CYCc and RAS1, respectively. The 3D interaction of QSSM 1157 and QSSM 1112 against CYCc and RAS1 is shown in Figure 3. Moreover, the data of virtual screening obviously demonstrated the phenomenal binding efficacy of QSMs from DKP QS system than other docked signaling systems. In contrast, a couple of earlier studies have reported the inhibitory efficacy of AHLs as well as DSFs toward C. albicans hyphal dimorphism in a dosedependent manner (Zhang et al., 2011;Grainha et al., 2020).
Despite having better binding energy, AHL could not occupy the top hit list owing to its weak binding profile when compared with farnesol. On the other hand, the well-reported structural and functional farnesol analogs of DSFs have displayed a very weak binding energy among the used signaling systems. The weak binding profile of DSFs toward target proteins (CYCc and RAS1) could be because of a distinct antihyphal mode of action different from that of the action mechanism of farnesol. In line with this, a previous finding has also demonstrated that, unlike AHLs and farnesol, the antihyphal mechanism of BDSFs (a B. cenocepacia QSM) follows a different mode of action (Hall et al., 2011).
The data of the current study unveiled the excellent binding mode of DKPs over essential hyphal proteins. DKPs are the simplest cyclic dipeptides that have been utilized for myriads of pharmaceutical propensities over the past few decades (de Carvalho and Abraham, 2012). The drug likeliness properties of DKPs attracted the attention of researchers to use them for a wide range of biological activities including antitumor, antiviral, antifungal, antibacterial, and antihyperglycemic efficacies (Rhee, 2004). Besides, the QS intervening efficacy of DKPs by acting as a QS signaling molecule during bacterial cross-talk has been reviewed critically by Bellezza et al. (2014).

Molecular Dynamics and Simulation
In order to understand the stability and interatomic motions of CYCc-QSSM 1157 and RAS1-QSSM 1112 complexes, MD   simulation analysis was performed using GROMACS, which would offer essential information regarding the physical movement of atoms and molecules in a biophysical system using computational approaches (Barnes and Williams, 1996). Therefore, the predicted protein-ligand complexes, i.e., CYCc-QSSM 1157 and RAS1-QSSM 1112 complexes, were subjected to 50-ns MD simulation. The subsequent MD run trajectories were examined for protein and ligand RMSDs, protein RMSF, and hydrogen bond interactions.

Root Mean Square Deviation and Root Mean Square Fluctuation Analyses
The stabilized interaction of protein-ligand complexes during MD simulation was evaluated by estimating their respective RMSD. Generally, the higher RMSD value signifies the instability of inspected molecules due to the high level of conformational changes within the complex. Additionally, the ligand molecule with a higher RMSD value against its respective protein-ligand complexes would suggest the insufficient binding of the ligand to the target protein (Opo et al., 2021). Thereby, RMSD analysis would reveal how often a protein-ligand interaction is integrated and stabilized. As shown in Figures 4C, D   QSSM 1157 and RAS1-QSSM 1112l are stable in their conformations.
To gain more insights into the stability of modeled proteins, the RMSF was estimated during the MD simulation. The RMSF of the Ca-atoms demonstrated the residual flexibility of each protein. Further, the results also showed the amino acid oscillations in catalytic and non-catalytic sites. For CYCc-QSSM 1157 and RAS1-QSSM 1112, the average protein RMSF was in the range of 0.1-0.73 and 0.06-0.47 Å, respectively ( Figures 4A, B). Further, negligible fluctuations were observed in disallowed and looping regions of CYCc and RAS1 proteins, which obviously signifies the stability of both predicted proteins throughout the simulation reactions.

H-Bond Interaction Analysis
As the effective binding of ligand to the target protein is mainly attributed to H-bond formation, the QSSM 1157 and QSSM 1112 were monitored for H-bond formation toward CYCc and RAS1 proteins throughout the simulation run. To identify the system stability during the simulation procedure, the number of hydrogen bonds formed vs. simulation duration was calculated. The hydrogen bonds stability was evaluated by means of 50-ns molecular simulation ( Figures 4E, F). The maximum number of hydrogen bonds formed for CYCc-QSSM 1157 and RAS1-QSSM 1112 was 9 and 7, respectively, indicating the greater affinity of ligand attachment to the protein binding pockets. The MD calculated hydrogen bond formation was in agreement with the result of molecular docking analysis. The number of hydrogen bonds formed by the QSSM 1157 and QSSM 1112 over CYCc and RAS1 throughout the simulation showed enough ligand binding strength to the proteins. Furthermore, the hydrogen bond ensures greater affinity and stability of ligand to the proteins.

Pharmacokinetic Prediction Using SwissADME Server
In the de novo drug development process, most of the drugs have failed in the final clinical phase due to their inappropriate pharmacokinetics and toxicity profile. To overcome this issue, pharmacokinetic and ADME prediction of the new drug candidate has to be conducted at the initial stages (Lee et al., 2003). To attain a desirable in vivo response, there should be an equilibrium between pharmacokinetic and pharmacodynamic qualities (Tuntland et al., 2014). Presently, in silico-based ADME screening is being utilized for the identification of the most promising compounds devoid of drug attrition.
In the present study, numerous parameters such as cell permeability, drug solubility, partition coefficient, gastrointestinal (GI) absorption, polar surface area (PSA), and drug homology were evaluated through virtual screening approaches provided by SwissADME interface. The results obtained in SwissADME are listed in Table 3, Supplementary Table 1, which contain ADMET, pharmacokinetics, drug likeliness, and medical chemistry of predicted ligands QSSM 1157 and QSSM 1112.

Drug Likeliness Prediction
Drug likeness of a ligand molecule is predicted by comparing the structure or other merit parameters with known drugs. In SwissADME, the drugability of the ligands was evaluated by comparing various properties with pre-developed drug likeliness modules such as Lipinski's rule of five, Ghose's rule, Veber's rule, Egan's rule, and Muegge's rule. Each of the mentioned modules contains different parameters that estimate ligand drug likeliness. Molecular weight, hydrogen donors, hydrogen acceptors, rotatable bonds, and PSA are various ligand properties that have been taken into account to evaluate a ligand drug's transport characteristics.
As can be seen in Table 3, the structures of QSSM 1157 and QSSM 1112 possess a better drug likeliness without any deviation in Lipinski's rule of five, Ghose's rule, Veber's rule, Egan's rule, and Muegge's rule as their values fall in the desirable range. To be precise, all the parameters, viz., molecular weight, rotatable bonds, number of hydrogen acceptors, number of hydrogen donors, PSA, molar refractivity, log P, and clogP, have fallen within the acceptable range. Furthermore, the QSSM 1157 and QSSM 1112 showed better lipophilicity (XLOGP3 between −0.7 and +5.0) and good water solubility scores. This result confirms that the chemical structure and orientation of both ligands will not cause any irrelevant metabolic interactions during drug formulation.

Pharmacokinetics Prediction
The SwissADME pharmacokinetics evaluation revealed that both the ligands have no blood-brain barrier (BBB) permeation and no substrate of permeability glycoprotein (P-gp) ( Table 3). Furthermore, the QSSM 1157 and QSSM 1112 showed high and low GI absorption, respectively.
The cytochrome P450 (CYP) is a hemeprotein that is responsible for the metabolism of drugs at various stages. Usually, CYP will have exhibited as five major isomers such as CYP1A2, CYP3A4, CYP2C9, CYP2C19, and CYP2D6. Inhibition of CYP pathway may lead to unwanted toxic effects due to reduced drug metabolism (Lynch and Price, 2007). Hence, the property of QSSM 1157 and QSSM 1112 to inhibit CYP was evaluated. The CYP interaction results showed that the QSSM 1112 is not a CYP inhibitor, whereas QSSM 1157 is an inhibitor of one of the isomers of CYP1A2. The obtained result clearly evidences that both QSSM 1157 and QSSM 1112 have a lesser possibility to interfere with CYP, which guarantees the drug selectivity and bioavailability. The Pan-Assay Interference compounds (PAINS) and Brenk technique were used to identify potentially toxic and unstable fragments present in the intact structure. In many biochemical and pharmacological trials, the PAINS evaluate the misleading hit compounds showing false-positive biological yields (Baell and Nissink, 2018). As expected, both QSSM 1157 and QSSM 1112 had "zero alerts" in PAINS prediction.
Although the QSSM 1112 does not have any violation in the medicinal chemistry profile, the Brenk list predicted QSSM 1157 as a Michael acceptor, indicating the possibility for a hazardous interaction of its fragments with other biological molecules. Since the QSSM 1157 showed inadequate toxicity in the Brenks-Michael acceptor alert, it is a necessity to evaluate the complete toxicity profile of the screened ligands. Hence, it is necessary to evaluate the toxicity profile of the screened ligands.

Toxicity Prediction via ProTox-II Server
The prediction of chemical toxicity is a crucial step in the drug development process. Computational-based toxicity prediction plays a significant role in fairly reducing the time and the number of animals to be used in in vivo studies for toxicity prediction (Raies and Bajic, 2016). ProTox-II uses 33 various models for predicting molecular similarity, fragment propensities, and machine learning to predict various toxicity endpoints including acute toxicity, hepatotoxicity, cytotoxicity, carcinogenicity, mutagenicity, and immunotoxicity.
The molecular weight distribution profile of both ligands was estimated by the ProTox-II server, and it was in the range of 258.32 to 319.67 for QSSM 1157 and 222.24 to 319.67 for QSSM 1112 (Figures 5A, 6A).
The toxicity was estimated by recognizing the ligand fragments and structure. The various toxicity profiles of both ligands estimated through ProTox-II server are given in Tables 4, 5 and Figures 5, 6. The toxicity radar chart in Figures 5C, 6C was intended to quickly illustrate the confidence of positive toxicity results compared with the average of its class. The LD50 values of QSSM 1112 and QSSM 1157 were predicted to be 1,600 and 2,287 mg/kg, respectively. The distribution of dose values for QSSM 1112 and QSSM 1157 with its mean dose value is represented in Figures 5A, 6A, respectively. The radar chart shows the bioavailability of the selected hit ligands. The pink area of the radar indicates the ideal area for properties such as flexibility (FLEX), insolubility (INSOLU), unsaturation (INSATU), lipophilicity (LIPO), polarity (POLAR), and SIZE. Both the selected ligands match with the recommended size of 500 g mol −1 as indicated by Lipinski for successful medication candidates. The POLAR was assessed for the selected ligands using their total area of polarity (TPSA). The TPSA value of a polar molecule should be 20 and 130 Å. The TPSA of QSSM 1157 and 1112 falls within the allowed range. Through INSOLU of the target ligands depicted with their respective ESOL and ESOL Class, it was found that QSSM 1157 and QSSM 1112 displayed high and moderate watersoluble properties, respectively.
The carbon fraction Sp3 (CSP3) and the rotatable bond count should not be greater than 9 and should be within the range of 0.5-1 (CSP3). It has been used to assess the INSATU and FLEX of the identified hit ligands. Both the selected hits QSSM 1157 (CSP3:0.33) and QSSM 1112 (CSP3: 0.45) showed CSP3 values fairly near csp3 range. In general, the numbers of rotatable bonds should not exceed nine; data showed that QSSM 1157 and QSSM 1112 possess 3 and 2 rotatable bonds, respectively. XLogP3 and ESOL (Log S) with the suggested range of − 0.7 to + 5.0 and 0 to 6, respectively, were used to access the LIPO and INSOLU profile of the selected ligands. Interestingly, both the selected ligands fall within the recommended range. Taken together, the data of pharmacokinetic and pharmacodynamic analyses of the selected hit molecules underline their highest plausibility for drug development.
Both the active agents were analyzed for various toxicity profiles and are listed in Tables 4, 5. QSSM 1157 showed mild toxicity possibility rates in hepatotoxicity, carcinogenicity, and cytotoxicity profiles, whereas QSSM 1112 demonstrated mild toxicity rates in carcinogenicity and cytotoxicity profiles.
As listed in Tables 4, 5, both ligands showed no toxicity in Tox21 receptor signaling pathway profiles. The overall toxicity class was predicted to be 4 out of 6 for QSSM 1112 and 5 out of 6 for QSSM 1157, which indicates the lesser toxicity of ligands. Overall, the analyses unambiguously predict that the identified QSM-based hyphal inhibitors would be safe for biological use.

In Vitro Antihyphal Assay Revealed Dose-Dependent Inhibitory Potency of QSSM 1112 Against Candida albicans Dimorphism
At the end of the virtual screening, the obtained result clearly showed that QSSM 1157 and QSSM 1112 have a similar kind of interaction like farnesol against CYCc and RAS1, respectively. Among the two predicted hits, biomolecule QSSM 1112 (active against RAS1) has been well studied and documented by our research group for its remarkable quorum quelling, antibiofilm, and antivirulence efficacies toward various multidrug-resistant pathogens, viz., S. mutans (Gowrishankar et al., 2014), S. aureus (Gowrishankar et al., 2016a), Listeria monocytogenes (Gowrishankar et al., 2016b), Staphylococcus epidermidis, and Serratia marcescens (Gowrishankar and Pandian, 2017;Gowrishankar et al., 2019). This information led us to be more curious to validate the in vitro efficiency of QSSM 1112 in inhibiting yeast-to-hyphal transition due to its in-house availability. Quite a few earlier studies on QSSM 1112 have already been reported for its synthesis (from Achromobacter xylosoxidans and Streptomyces sp.) and exhibition of antagonistic activity against aflatoxins of Aspergillus parasiticus and rice blast fungus Pyricularia oryzae (Yan et al., 2004;Hall et al., 2011). Except for the anticandidal effect of QSSM 1112 (Rhee, 2004), so far, none of the studies has been focused on its inhibitory potency toward C. albicans virulence dimorphism. Hence, in order to accomplish our curiosity in validating the in silico docking result, in vitro antihyphal assay was performed using QSSM 1112, with farnesol as the positive control.
For this, C. albicans cells were cultured with and without the active agent (QSSM 1112) in liquid spider medium supplemented with 10% of hyphal inducer (FBS). As expected, the result of the liquid antihyphal assay signified that QSSM 1112 was potent enough to inhibit the yeast-to-hyphal transition in a concentration-dependent manner ( Figure 7A). Of the various concentrations (0-1,024 µg/ml) used, it was found that QSSM 1112 successfully thwarts the hyphal formation from the maximu m (1,024 µg/m l) to minimu m (64 µg/ml) concentrations. Therefore, the lowest concentration that brings a substantial hyphal inhibitory condition was considered to be the minimal HIC. Consequently, a concentration of 64 µg/ml was considered as HIC of QSSM 1112, as it demonstrated a phenomenal antihyphal efficacy in a similar fashion to higher concentration (1,024 µg/ml) as well. As depicted in Figure 7A, the cells devoid of active compound (untreated control) have displayed a highly dense criss-cross architecture of hyphal cell population along with few yeast cells. On the other hand, the cells exposed with QSSM 1112 have showcased only yeast cell population even under strong hyphae-induced conditions, which clearly proves that the existence of QSSM 1112 disarmed the C. albicans cells to undergo hyphae induction. Moreover, the consistency of QSSM 1112 in expressing the sustained hyphal inhibitory potency even in the presence of serum would make the QSSM 1112 an interesting drug candidate in the treatment against invasive C. albicans infection. Farnesol was used as the positive control in all the in vitro assays. In a similar fashion to QSSM 1112, farnesol also exhibited a profound concentration-dependent inhibition on filamentous growth of C. albicans. At the lowest concentration of 32 µg/ml, the farnesol has shown complete hyphal inhibition (0% hyphal) ( Figure 7C). Although the influence of farnesol on the C. albicans dimorphic switching has already been demonstrated, the exact concentration required for yeast-to-hyphal inhibition varied from study to study owing to multiple associated factors. Mosel et al. (2005) reported that the concentration of farnesol required for the inhibition of yeast-to-hyphal transition has varied depending on the medium (Mosel et al., 2005). They found that 2 µM (0.44 µg/ml) and 250 µM (55.29 µg/ml) of farnesol were required to inhibit 50% of C. albicans hyphal formation in defined media and serum-containing media, respectively (Mosel et al., 2005). However, our study demonstrated that 32 µg/ml of farnesol was potent enough to inhibit the filamentation of C. albicans in both spider medium and serum. To further confirm the propensity of QSSM 1112 in inhibiting the C. albicans dimorphic behavior, solid spider agar assay was done. For that, overnight C. albicans culture was spotted on the spider agar medium supplemented with and without QSSM 1112 at its HIC ( Figure 7B). After 2 days of incubation, the agar plates were directly visualized under the light microscope to observe the change in colony morphology upon QSSM 1112 treatment. The micrograph of untreated control has shown rough colony morphology with densely elongated hyphal cells. At the same time, the micrographs of QSSM 1112 and farnesol-treated cells have portrayed smooth colony morphology with abridged hyphal formation. In addition, after 7 days of incubation, the agar plates were documented using a gel documentation system ( Figure 7C). The results signified the remarkable reduction of hyphal protrusion manifested with QSSM 1112 compared with the untreated control.
The hyphal inhibitory effect of farnesol was better than that of QSSM 1112, as it exhibited a substantial impact even at a very low concentration, i.e., even 2 µg/ml with 25% of hyphal inhibition. In contrast, QSSM 1112 at 2 µg/ml displayed a compactly elongated hyphal formation similar to that of control.

QSSM 1112 (at Hyphal Inhibitory Concentration) Averts Candida albicans Germ-Tube Formation
The dimorphic growth of C. albicans is closely linked with the thin filamentous outgrowth extending from the blastospore known as "germ-tube." The extended germ-tube later develops into true hyphae. Hence, the efficacy of QSSM 1112 (HIC) on C. albicans germ-tube formation was evaluated by culturing the C. albicans cells in FBS. In parallel, the effect of farnesol on germtube formation was evaluated for comparative purposes. The C. albicans cells were imaged in different time intervals (from 0 to 6 h) for comparison of germ-tube formation with and without QSSM 1112. As can be seen in Figure 8A, both QSSM 1112treated and untreated control cells appeared as yeast form at the 0-h time interval. With the increase in time, few numbers of cells in untreated control samples showcased budding yeast morphology (a thin filamentous tube extending from yeast), and the true germ-tube germination was observed in the control sample at the time interval of 6 h. However, the cells treated with QSSM 1112 and farnesol (at HIC) portrayed complete spread of yeast cells, lacking the classic characteristics of hyphal elongation morphologies, viz., budding yeast followed by badminton shaped germ-tubes. This inevitably reinforces the data of in vitro liquid and solid antihyphal assays, which in turn signifies that QSSM A B FIGURE 8 | (A) The influence of QSSM 1112 (at HIC) and farnesol on Candida albicans germ-tube formation at various time intervals (0-6 h). The QSSM 1112 manifested into the FBS along with the C. albicans cells. The micrograph evidently demonstrates the blockage of germ-tube formation in C. albicans cells treated with QSSM 1112 compared with that of untreated control cells (which encompasses classic badminton structures, i.e., true germ-tubes of C. albicans). (B) The impact of QSSM 1112 (at HIC) and farnesol on preformed germ-tube formation. There is no significant difference between untreated control and QSSM 1112-and farnesol-treated cells. FBS, fetal bovine serum; HIC, hyphal inhibitory concentration. 1112 followed the same mode of hyphal inhibitory action as that of C. albicans QS signaling molecule farnesol. In line with the present study, a similar kind of bacterial QSMs mediated germtube inhibition was observed by Zhang et al. (2011), wherein the two DSFs, i.e., cis-and trans-2-dodecenoic acid, have arrested the C. albicans cells from undergoing germ-tube induction (Zhang et al., 2011).
In most of the cases, the antihyphal inhibitors not only halt the germ-tube germination but also have the capability to revert the preformed germ-tube into yeast form. In order to know whether the QSSM 1112 exposure has any effect on preformed germ-tube, QSSM 1112 was added to wells that had preexisting C. albicans germ-tube. After 5-h time exposure, the preformed germ-tube cells were visualized under a light microscope. As can be seen from Figure 8B, no considerable change was observed in the germ-tube formation of both QSSM 1112-treated and untreated control cells. It was also observed that the positive control farnesol also had no effect; it neither reverted nor destroyed the germ-tube cells. This further confirms that QSSM 1112 and farnesol are ably averting C. albicans dimorphism before germination, whereas they could not have any effect on C. albicans cells that had germinated already. This result is in agreement with the observation of Mosel et al. where the farnesol appears to be ineffective on preexisting hyphae even at a very high concentration (Mosel et al., 2005). Although both farnesol and QSSM 1112 did not effectively inhibit preformed germ-tube formation, it has been reported that exposure with farnesol could sensitize the drug-resistant C. albicans to antifungal treatment. With the use of this same principle, several studies also demonstrated the synergistic combinatorial effect of farnesol over the conventional antifungals (amphotericin B, fluconazole, and micafungin) against C. albicans biofilm formation. In this way, QSSM also could be used in synergistic therapy to enhance the efficacy of a single dose.

Fluorescence Microscopic Analysis Reinforces the QSSM 1112-Mediated Hyphae and Germ-Tube Inhibition
In order to further reinforce the data obtained from a light microscopic analysis, fluorescence microscopic analysis was carried out. As seen in Figures 9A, B, QSSM 1112 is sufficient enough to inhibit the germ-tube and hyphal formation in C. albicans. Through in vitro validation of in silico results, the present investigation stands as a "proof of concept" for the reliability of in silico-based drug prediction in the drug development process.

CONCLUSION
The current study investigates the antihyphal efficacy of bacterial QSMs through in silico and in vitro analyses against C. albicans dimorphic switching. Results of initial virtual screening analysis revealed hyphal inhibitory potential of the two DKPs, viz., QSSM 1112 and QSSM 1157, against CYCc and RAS1, respectively. Further, 50-ns MD simulation demonstrated the stable binding of ligands to target proteins with less conformational fluctuation and greater affinity. In addition, the pharmacokinetics and toxicity prediction revealed that the identified QSM hyphal inhibitors are safe for biological use. These QSMs displayed better lipophilicity and water solubility scores, lesser toxicity, and no BBB permeation. Further, in vitro investigations such as liquid and solid antihyphal assays and germ-tube inhibition assays unequivocally established QSSM 1112 as a promising hyphal inhibitor. On the whole, through in silico and in vitro analyses, the present study opens up new avenues for considering the bacterial QSMs as potential inhibitors against fungi especially C. albicans. This study also highlights that the antagonistic interaction between bacteria and fungi can be utilized as a tool to overcome certain fungal infections.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

AUTHOR CONTRIBUTIONS
RJ designed the research idea, analyzed the data, carried out in vitro analysis and result interpretation, and wrote the original draft. NH carried out virtual screening, molecular docking and simulation analysis, and in silico result interpretation and prepared the figures. SG designed the research idea; performed the formal analysis, funding acquisition, investigation, administration, and validation; and critically revised the original draft. SP contributed to the conceptualization, supervision, funding acquisition, and critical revision of the original draft. All authors contributed to the article and approved the submitted version.