Curcumin Analogues as a Potential Drug against Antibiotic Resistant Protein, -Lactamases and L, D-Transpeptidases Involved in Toxin Secretion in Salmonella typhi: A Computational Approach

Typhoid fever caused by the bacteria Salmonella typhi gained resistance through multidrugresistant S. typhi strains. One of the reasons behind β-lactam antibiotic resistance is -lactamase. L, D-Transpeptidases is responsible for typhoid fever as it is involved in toxin release that results in typhoid fever in humans. A molecular modeling study of these targeted proteins was carried out by various methods, such as homology modeling, active site prediction, prediction of disease-causing regions, and by analyzing the potential inhibitory activities of curcumin analogs by targeting these proteins to overcome the antibiotic resistance. The five potent drug candidate compounds were identified to be natural ligands that can inhibit those enzymes compared to controls in our research. The binding affinity of both the Go-Y032 and NSC-43319 were found against β-lactamase was −7.8 Kcal/mol in AutoDock, whereas, in SwissDock, the binding energy was −8.15 and −8.04 Kcal/mol, respectively. On the other hand, the Cyclovalone and NSC-43319 had an equal energy of −7.60 Kcal/mol in AutoDock, whereas −7.90 and −8.01 Kcal/mol in SwissDock against L, D-Transpeptidases. After the identification of proteins, the determination of primary and secondary structures, as well as the gene producing area and homology modeling, was accomplished. The screened drug candidates were further evaluated in ADMET, and pharmacological properties along with positive drug-likeness properties were observed for these ligand molecules. However, further in vitro and in vivo experiments are required to validate these in silico data to develop novel therapeutics against antibiotic resistance.


Introduction
Typhoid is a usual illness in economically handicapped countries where public health settings are very poor. A globally estimated 12-27 million people get stricken with typhoid fever each year, whereas the overall yearly estimated incidence lies between 292 and 395 cases per 100,000 people in Bangladesh. This infection-causing agent is an anaerobic Gram-negative rod, namely Salmonella enterica serotype Typhi (S. typhi), a highly conserved serovar subspecies of S. enterica, which is transmitted by the fecal-oral route and can infect connected by small peptides. Here, the enzyme L, D-Transpeptidases plays its role in introducing cross-links within L-and D-amino acids that comprise the peptides (Figure 1). This PG remodeling by L, D-Transpeptidases is necessary for TtsA to position the typhoid toxin for its proper release [32]. Here, L, D-Transpeptidases can be a major target for in silico studies as it plays a vital role in the secretion of typhoid toxin, and there is no effective drug available to inhibit it without exhibiting any side effects. For example, drug carbapenem and copper can inactivate L, D-Transpeptidases, yet they are associated with diarrhea, nausea, vomiting, skin rash, low blood pressure, anemia, heart problems, etc. [36]. Moreover, several antibiotics are working alone or coupled with β-lactamase inhibitors (Avibactam, Clavulanic acid (clavulanate), Relebactam, Sulbactum, Tozobactum, etc.), which have many adverse effects such as gastrointestinal complications, impairment of nervous system, hematological effects, and dermatological abnormalities, including Stevens-Johnson syndrome, toxic epidermal necrolysis, and drug-induced eosinophilia, etc. [37][38][39][40][41][42][43]. Curcumin, the main bioactive component of turmeric (Curcuma longa L.), has been shown to be a powerful antioxidant, anti-inflammatory, antibacterial, antifungal, and antiviral agent in many studies [44]. Curcumin has been shown to be antibacterial against Staphylococcus aureus (S. aureus). Curcumin has significantly more effective antibacterial properties when combined with other antibacterial drugs, as revealed by in vitro experiments [45]. Curcumin inhibits bacterial growth due to its structural properties and the production of anti-oxidative chemicals. Through the bacterial quorum sensing regulatory system, curcumin can decrease bacterial virulence factors, reduce bacterial biofilm formation, and restrict bacterial adherence to host receptors [46]. Curcumin's potential antibacterial action makes it a viable option for enhancing the inhibitory impact of current antimicrobial drugs through synergism [47]. It decreased Salmonella enterica serovar Typhimurium's motility by reducing the length of the flagellar filament (from 8 m to 5 m) and lowering its density (4 or 5 flagella/bacterium instead of 8 or 9 flagella/bacterium). Curcumin therapy reduced the proportion of flagellated bacteria from 84 percent to 59 percent [48].
As curcumin has antibacterial properties, including antioxidant and anti-inflammatory properties, we selected 70 curcumin analogues as ligands in this study. The schematic representation of the methodology applied in the present study is displayed in Figure 2.

Drug-Likeness Features Interpretation
About 70 curcumin molecules were selected as ligands from literature to target the proteins β-lactamase and L, D-Transpeptidases of S. typhi. Each of the molecules was assessed for Lipinski's rule of five or not [49,50]. Molinspiration Cheminformatics server (https://www.molinspi-ration.com/cgi-bin/properties accessed on 8 October 2020) was applied to experiment various drug-like parameters of the ligand molecules [51,52] (Table S2). Molinspiration allows for the prediction of significant molecular parameters (logP, polar surface area, number of hydrogen bond donors and acceptors, and so on), as well as the prediction of bioactivity scores for the most relevant therapeutic targets (GPCR ligands, kinase inhibitors, ion channel modulators, nuclear receptors). Compounds that did not comply with the rule were excluded from further study.

Protein Preparation for Docking Sequence Retrieval
The complete protein sequences of β-lactamase and L, D-Transpeptidases were retrieved from NCBI (https://www.ncbi.nlm.nih.gov/ accessed on 12 October 2020) in the standard FASTA format.

Physiochemical Property Identification
The physical and chemical parameters of the proteins, including molecular weight (MW), theoretical pI, amino acid composition, estimated half-life, instability index, aliphatic index, etc., were computed using the ProtParam tool of ExPasy server as it evaluates physicochemical data (molecular weight, theoretical pI, amino acid composition, atomic composition, extinction coefficient, estimated half-life, instability index, aliphatic index, and grand average of hydropathicity (GRAVY)) from a protein sequence. (http://web. expasy.org/program/ accessed on 12 October 2020) [53].

Secondary and Tertiary Structure Prediction
The online tool SOPMA [54] (https://npsa-prabi.ibcp.fr/NPSA/npsa_sopma.html accessed on 15 October 2020) was applied for the secondary structure prediction of the proteins. Homology Modeling was performed in Swiss-Model, as its workspace is a Webbased modeling expert system that is integrated. We search a library of experimental protein architectures for acceptable templates for a specified target protein. A three-dimensional model of the target protein is constructed based on a sequence alignment between the target protein and the template structure [55][56][57][58][59] (https://swissmodel.expasy.org/, accessed on 18 October 2020) to construct the tertiary structures of the target proteins using the threedimensional structure of a related protein as a template. Homology modeling was also done by Phyre2. It is a web-based collection of tools for predicting and analyzing protein structure, function, and mutations. Phyre2's goal is to provide biologists with a simple and intuitive interface to cutting-edge protein bioinformatics tools. Phyre2, which builds 3D models, predicts ligand binding sites and analyzes the influence of amino acid changes (e.g., no synonymous SNPs (nsSNPs)) for a user's protein sequence using advanced distant homology detection algorithms [60].

Disordered Regions Prediction
Disordered regions present in protein molecules remain unstable in the native state. To find out the disordered regions in proteins for which they lack a fixed tertiary structure, the Protein Disorder prediction System (PrDOS) server [61] (http://prdos.hgc.jp/cgi-bin/ top.cgi, accessed on 18 October 2020) was exploited. This server predicts the disordered regions based on both local amino acid sequence and the template or homologous proteins through the SVM algorithm and PSI-BLAST, respectively. The prediction method consists of two predictors: one based on local amino acid sequence information and the other on template proteins. For each residue, the server aggregates the findings of the two predictors and delivers a two-state prediction (order/disorder) and a disorder probability.

Validation of Tertiary Protein Model
Model validation was carried out in PROCHECK [62,63] (http://www.ebi.ac.uk/ thornton-srv/software/PROCHECK/, accessed on 20 December 2021), which assesses the stereo-chemical quality of a protein structure, i.e., how normal or unusual the pattern of the protein residues is, compared with a fined-tuned, high-resolution structure of a protein. PROCHECK evaluates the stereochemical quality of a protein structure by generating a series of PostScript graphs that analyze its overall and residue-by-residue geometry. It contains PROCHECK-NMR, which is used to check the quality of structures solved by NMR.

Active Site Prediction
For the prediction of active sites in the proteins where ligands will likely bind, two servers CAST P (Computed Atlas of Surface Topography of proteins) [64] (http://sts.bioe. uic.edu/castp/ accessed on 20 October 2020) and COACH [65,66] (http://zhanglab.ccmb. med.umich.edu/COACH/ accessed on 22 October 2020) were utilized. CAST P implements the theoretical and algorithmic results of computational geometry to predict the ligandbinding sites. It has several advantages: (1) pockets and cavities are recognized analytically, (2) the boundary between the bulk solvent and the pocket is accurately specified, and (3) all derived parameters are rotationally invariant, do not need discretization, and do not make use of dot surface or grid points. On the contrary, the COACH server applies two comparative methods, TM-SITE and S-SITE, to identify active sites in the protein.

Ligand and Protein Preparation for Docking
The ligand molecules were minimized in the Avogadro software using the mmff94 force field. Then, the protein structure was minimized in YASARA software using the AMBER14 force field. The docking program was carried out in the AutoDock Vina program. The ligands that will give the best results will be docked again using the SwissDock server. We used two servers to check the validity and to build a strong hypothesis. AutoDock Vina, a novel molecular docking and virtual screening application, has been introduced. Vina's local optimization process employs a powerful gradient optimization algorithm. The gradient computation essentially provides the optimization algorithm with a "feeling of direction" from a single evaluation. Vina may speed up processing by making use of many CPUs or CPU cores by employing multithreading. SwissDock is a webserver dedicated to doing protein-ligand docking simulations in an easy and beautiful manner. SwissDock is protein-ligand docking software with a simple and integrated interface that is based on EADock DSS.
AutoDock Vina PDB file of protein was converted into a pdbqt file using the Auto Dock Vina tool. Then, every selected ligand file was converted to a PDB file using pymol because the AutoDock tool can only recognize the PDB format. Then, the PDB file of the ligand was converted into a pdbqt file, which is the criterion for the AutoDock run. After docking, we visualized the docked file using pymol.
SwissDock: At first, we set up the protein PDB file manipulating PDB code, which was retrieved from Uniprot or RCSB PDB. We specified chain A in β-lactamase. Then, ligands were selected from the ZINC database, but those that were not present in this database were uploaded in Mol2 format. The server provides very fast, fast, and accurate results.

ADME/T Prediction
ADME/T describes the Absorption, Distribution, Metabolism, Excretion and Toxicity of a drug-like substance. These properties account for the success of a drug in clinical trials. Therefore, in silico ADME/T profile examination of the candidate drugs is a prerequisite for the fruitful measure of drug designing expenditure [67,68]. The best 8 ligands (based on the docking score) were utilized to speculate their drug-like potential by observing pharmacokinetic and pharmacodynamics features. ADME/T profile of all the chosen ligands were calculated utilizing admetSAR 2.0 (http://lmmd.ecust.edu.cn/admetsar2/ accessed on 23 October 2020) and pkCSM (http://biosig.unimelb.edu.au/pkcsm/ accessed on 23 October 2020) server [69,70]. admetSAR presents an easy-to-use interface for searching for ADME/T (Absorption, Distribution, Metabolism, Excretion, and Toxicity) attributes profiling by name, CASRN, and similarity search. With QSAR models, admetSAR can predict around 50 ADMET endpoints. The pkCSM signatures were effectively employed to create predictive regression and classification models across five different pharmacokinetic property classes.

Validation of Tertiary Protein Model
ERRAT software was used to validate the protein model. ERRAT is a software that verifies crystallographically determined protein structures.

Pharmacological and Biological Activity Prediction
The elected ligands were employed to determine their pharmacological and biological activities accurately by using Prediction of Activity Spectra of Substances (PASS) Online (http://www.pharmaexpert.ru/passonline/ accessed on 26 October 2020) and Molinspiration Cheminformatics server [71]. These methods are used in conjunction with recognized compounds present in the database, depending on the structure-activity relationship (SAR). PASS Online predicts about 4000 different types of biological activity, such as pharmacological effects, mechanisms of action, toxic and unfavorable effects, interactions with metabolic enzymes and transporters, gene expression influence, and so on. Molinspiration Cheminformatics is also useful software to predict pharmacological and biological activities.

Pred. P450 Site of Metabolism Iction
In silico methods can contribute significantly to the prediction of drug metabolism sites focusing on the experimental view of the drug designing process. By conducting bioassay, these sites impart the knowledge of the molecules' metabolic susceptibility and their fate inside the body [72]. RS-WebPredictor (http://reccr.chem.rpi.edu/Software/ RS-WebPredictor/ accessed on 1 November 2020), an online server, was used to predict the best sites of drug metabolism mediated by CYP2C9, CYP2D6, and CYP3A4, three promiscuous isoforms of Cytochrome P450 (CYP) family. Predictions may be made for the promiscuous CYP isozymes 2C9, 2D6, and 3A4, as well as CYPs 1A2, 2A6, 2B6, 2C8, 2C19, and 2E1. The RS-WebPredictor service is the first publicly available server that predicts the regioselectivity of the last six isozymes.

Drug-Likeness Features Interpretation
Among 70 types of ligands, all the compounds except DM1 were compatible with the drug-likeness features. Most of the selected compounds did not display any violations to Lipinski's rules, which indicated the pharmacokinetic conformity of these compounds. Hence, all these compounds were taken into account for the next phases of the study (Tables S1 and S2).

Physiochemical Property Identification
Using the ProtParam tool of the ExPasy server, the physical and chemical parameters of the proteins were analyzed. β-lactamases is 286 amino acids long, weighted at 31515.20 grams, and has an instability index of 40.74. L, D-Transpeptidases is 615 amino acids long, 67812.49 (gm) weighted, with an instability index value of 43.79. The Instability Index is a metric for determining whether a protein will remain stable in a test tube. Table 1 shows the physicochemical properties of these two proteins.

Secondary and Tertiary Structure Prediction
The SOPMA tool was used to predict the secondary structures of these two proteins. The values of alpha helix were 49.30% and 39.35% for β-lactamase and L, D-Transpeptidases, respectfully. The value of the extended strand was also greater for β-lactamase (12.94%) compared with the L, D-Transpeptidases (11.54%). L, D-Transpeptidases (43.58%) exceeds the β-lactamase (29.37%) in the case of a random coil. The values of 3 10 helix, pi helix, beta bridge, bend region, and ambiguous status were 0.00% for both proteins (Table 2). The three-dimensional structures of β-lactamase and L, D-Transpeptidases were predicted in Swiss-Model web tools. The biounit oligo state of both proteins was a monomer. The template displayed 0.61 sequence similarities, coverage score of 1.0, and 24-286 range for β-lactamase, whereas L, D-transpeptidases showed 0.62 sequence similarity, coverage score of 0.95, and 37-615 range (Table 3). In Swiss-Model, two models for β-lactamase and three models for L, D-Transpeptidases were predicted based on the top 31 and 50 templates, respectively. On the contrary, the top 20 models were predicted by Phyre 2 for both proteins each. The best fit built by the two servers for β-lactamase had a confidence score of 100 when modeling 263 amino acid residues at positions 24-286, and for L, D-transpeptidases, it modelled 505 residues at positions 37-615 with a confidence score of 100.

Validation of Tertiary Protein Model
In ERRAT, the overall quality factor of β-lactamase and L, D-Transpeptidases was is 97.2549 and 84.1141, respectively. The two proteins also passed the verified 3D in their respective prediction results. The β-lactamase (EC 3.5.2.6) has an 11% disease-causing region, and its active sites are acyl ester intermediate (position 70) and proton acceptor (position 168). The Ramachandran plot analysis showed that both proteins delineated more than 90% of the amino acid residues in the most favored regions. The number of nonglycine and non-proline residues is 228 among 263 residues. Twenty-one glycine residues and 12 proline residues are present (Ramachandran plot), and 93.4% of residues are in the favored region. The L, D-Transpeptidases is 615 amino acids long, and its molecular weight is 67812.49. The active site of this protein stays in the 528 position. Here, 91.8% residues are in the favored region with 427 residues excluding glycine and proline, 8 terminal residues other than Gly and Pro, 31 glycine (represented as a triangle), and 3 proline residues (Table 4 and Figures 3 and 4).

Active Site Prediction
Using CAST P (Computed Atlas of Surface Topography of proteins) and COACH servers, the active sites of β-lactamase and L, D-Transpeptidases were predicted ( Figures 3F and 4F).

Molecular Docking
The best five ligands were selected based on docking experiments among 70 ligands. Go-Y032, NSC-43319, Cyclovalone, Salsalate, and Cyclocurcumin showed the best docking results against the β-lactamase (1ZG4) enzyme (Tables 5 and 6 and Figure 5), and NSC-43319, Cyclovalone, Cyclocurcumin, Difluorinated curcumin, and Go-Y032 showed the best docking results against the L, D-Transpeptidases (6NTW) enzyme (Table 5, Table 6 and Figure 5). Two controls were selected for each of the enzymes. Clavulanic acid and Tazobactum were docked against β-lactamase as they are existing drugs. These existing drugs showed low affinity to other ligands. For L, D-Transpeptidases, Carbapenem and Cephalosporin were selected as controls. These controls were also showed lower affinity.        3.2.7. ADME/T Prognosis ADME/T profiling was carried out for the ligands that gave the best docking scores in the molecular docking study (Table S6) and control group (Tables S3 and S4). Intestinal absorption and oral bioavailability were high for all of the ligands. Caco-2 permeability was high for the ligands Go-Y032 and NSC-43319, while the remaining ligand showed low Caco-2 permeability. Furthermore, NSC-43319 and Cyclovalone were predicted as substrates of P-glycoproteins, while the remaining was non-substrates of membrane Pglycoproteins. All of the ligands were able to enter the blood-brain barrier, excluding Cyclocurcumin. All of them exhibited substrate specificity of CYP3A4 except Cyclovalone. The NSC-43319, Salsalate, and Cyclovalone showed no specificity as substrates of CYP3C9. All of them were negative for AMES toxicity, and compounds demonstrated the inhibitions of the Human ether-a-go-go related gene (hERG) channel and also a less toxic profile.

Prediction of Pharmacological and Biological Activity
The screened ligands were inspected for the pharmacology-related study (Table S5) involving their accordance with Antibacterial, Bacterial Efflux Pump Inhibitor, Antibiotic Anthracycline-like activity, β-lactamase Inhibitor, Anti-mycobacterial and antibiotic activities. Here, Cyclovalone showed all of these pharmacological activities, whereas Go-Y032 was observed to have all of these properties except Antibiotic activity. The NSC-43319 showed antibacterial, bacterial efflux pump inhibitor, β-lactamase inhibitor, and anti-mycobacterial activities. In this study, Go-Y032, NSC-43319, and Cyclovalone were found as the best-performing ligands (Table 7). Thereafter, these five ligands were analyzed to observe whether they function against G protein-coupled receptor (GPCR) ligand, protein kinase, ion channels, enzyme protease, nuclear receptor ligand, etc., (Table S6).

Prediction of P450-Mediated Sites of Metabolism (SOMs)
The most potential ligands were investigated for the prediction of their possible sites of metabolism (SOMs) by three main Cytochrome P450 (CYP) isozymes, i.e., 3A4, 2C9, and 2D6, respectively (Table 8). For CYP2D6 and CYP2C9 isoforms, Go-Y032, and for CYP3A4 and CYP2C9 isoforms, NSC-43319 have shown identical metabolism sites. However, for all of the isozymes, cyclovalone exhibited different metabolism sites.

Discussion
Molecular docking is the most significant approach in in silico drug designing. It assesses the binding affinity of a protein-ligand complex in the form of binding energy using computer algorithms. The lower the binding energy, the higher the affinity of the ligand bound to the target [73]. Five ligands gave the best free binding energies in the docking experiment conducted by AutoDock Vina and SwissDock, which included Go-Y032, NSC-43319 with β-lactamase enzyme, and NSC-43319, Cyclovalone with L, D-Transpeptidases (Tables 5 and 6). The docking results were compared with controls. For the β-lactamase enzyme, the docking results with other ligands were compared with Clavulanic acid and Tazobactum (Table 9). Clavulanic acid and tazobactam are all plasmidmediated β-lactamase inhibitors. Several studies have concluded that Clavulanic acid inhibits extended-spectrum TEM and SHV β-lactamases. They expressed lower affinity in our study compared to other ligands. Several studies reported that β-lactam antibiotics could work against L, D-Transpeptidases. Carbapenem and Cephalosporin are antibiotics in the beta-lactam class that kill bacteria by attaching to penicillin-binding proteins and blocking bacterial cell wall formation. As these are existing drugs, we compared their activities with our selected ligands. These controls were showed very low binding affinities in docking. Therefore, the ligands we selected have a better chance of working against those enzymes or as antibacterial drugs. After docking, the ligands were introduced to ADME/T prediction tools admetSAR 2.0 and pkCSM software for in silico prediction of ADME/T properties. This time-saving and cost-effective approach helps establish a candidate molecule as a promising drug through in vitro experiments [74,75]. A leading concern for drugs that mainly target the central nervous system (CNS) is that they need to be able to permeate across the blood-brain barrier. The most common route for drug delivery is the oral delivery system, through which the delivered drug enters into the intestine; hence, an investigated drug needs to be highly absorbed in the intestine. Cell membrane proteins, such as P-glycoproteins, facilitate the movement of many drugs through the cell membrane. Caco2 permeability of a drug reflects the permeability across the intestinal lining of humans, as this tissue is widely used in in vitro drug permeability studies due to its small intestinal mucosa-like behavior when cultured [75][76][77]. The Cytochrome P450 enzyme family focuses on the regulation of drug interaction, biotransformation, and their elimination outside the body. Acute toxicity, delayed removal, and eventual drug compound failure within the human body result from these enzymes' inhibitory activity of the drugs [78][79][80]. The purpose of in silico AMES toxicity test is to determine the toxicity and mutagenicity of the chemicals [81,82]. Voltagegated potassium ion channels, namely hERG channels, are involved in the transport of potassium ions through the cell membrane, which may be subjected to off-target drug interaction resulting in inhibition. Therefore, proper screening is necessary to see whether the investigated drugs have inhibitory activity on these transporters [83]. Renal organic cation transporter 2 (OCT2) plays a major role in the removal of drugs and xenobiotics via the kidney. It is considered that the substrates of this transporter protein are quickly excreted by urine [84]. In the ADME/T test, all the selected ligands showed almost similar properties (Table S6).
Afterwards, the pharmacological and biological activities of the ligands were carried out in the PASS online server and Molinspiration Cheminformatics server, respectively. Pharmacological activity (PASS prediction) is determined in terms of the likelihood of activity (Pa) and the likelihood of inactivity (Pi) of a drug, and the result of the prediction ranges between 0 and 1. The pharmacological activity of the drug is deemed possible if Pa > Pi [85]. The probability of anti-mycobacterial activity (Pa) for NSC-43319 and Cyclovalone was between 0.5-0.7, while for all ligands, Pa of all activities was <0.5, implying the unlikeliness of their activities [86]. However, Cyclovalone, Go-Y032, and NSC-43319 showed more satisfactory outcomes than others (Tables S5 and S6). Assessing biological activities against the most influential drug targets in the human body, such as G protein-coupled receptors (GPCRs), ion channels, enzymes, nuclear receptors, etc., is crucial because, when coupled with them, a drug mediates its therapeutic activity inside the body [87]. Probability scores for Petasiphenol representing activity against the targets were comparatively significant (Table S5).
Finally, the ligands were assessed in the RS-WebPredictor server to predict the probable sites where their metabolism will be likely to occur. Almost similar metabolism sites were reported for Go-Y032, NSC-43319, Cyclovalone, Difluorinated curcumin, which exhibited multiple sites of metabolism except for Cyclocurcumin and Salsalate. These two compounds showed few sites of metabolism compared to others (Table 8). After a successive screening of all potent ligands, this study suggests Go-Y032, Cyclovalone, NSC-43319, Difluorinated curcumin, Salsalate, and Cyclocurcumin as the best inhibitors of the enzymes β-lactamase and L, D-Transpeptidases, respectively, while NSC-43319, Cyclovalone, Cyclocurcumin, and Go-Y032 are recommended for both. These compounds may not express better performances at all screening tests. As our study was based on a computational approach, further in vitro analysis is needed to prove their activities. These compounds may give better performance in an in vitro approach than the controls we chose.

Conclusions
This study aimed to develop curcumin derivatives as potential drug candidates against S. typhi causing typhoid fever. The target proteins included two crucial proteins of this bacterium, β-lactamase and L, D-Transpeptidases, which are involved in acquiring antibioticresistance properties and toxin secretion, respectively. Although several antibiotics and β-lactam inhibitors are available to treat this fatal illness, none is effective enough to avoid side effects that eventually result in adverse complications. After surpassing multiple stages in assessing drug-like properties, among 70 ligands, three, including Go-Y032, Cyclovalone, and NSC-43319, were reported as the best-performing ligands. Further investigations based on in vivo and in vitro experiments are needed to ascertain the use of these ligands as drugs to treat S. typhi infection. Additionally, the other four ligands providing a satisfactory docking performance are also recommended for further wet laboratory investigation.