Molecular dynamics simulation of tryptophan hydroxylase-1: binding modes and free energy analysis to phenylalanine derivative inhibitors

Lasiodiplodiapyrones A and B, Pyrone-Preussomerin Adducts with Highly Strained Polycyclic Ring Systems from Lasiodiplodia pseudotheobromae

Lasiodiplodiapyrones A and B (1 and 2), two new preussomerin derivatives, possessing an unexpected 6-methyl-4H-furo[3,2-c]pyran-4-one moiety and a highly functionalized conjoint and complicated polycyclic ring system, along with two known congeners (3 and 4), were isolated from the fungus Lasiodiplodia pseudotheobromae. Their structures including absolute configurations were determined by spectroscopic analyses, Mosher's method, and ECD calculations. A biosynthetic pathway was proposed to explain the origin of lasiodiplodiapyrones A and B as well as their relationship with preussomerins. Compounds 1-4 showed suppressive effects on the production of NO with IC₅₀ values of 4.8 ± 0.3, 8.5 ± 1.1, 5.9 ± 0.8, and 12.8 ± 1.3 μM, respectively.

Four undescribed ergostane-type steroids from Lasiodiplodia pseudotheobromae and their neuroprotective activity

Four undescribed ergostane-type steroids, (22E,24R)-4α,5α-epoxyergosta-9α,14β-dihydroxy-7,22-diene-3,6-dione, (22E,24R)-4α,5α-epoxyergosta-9α,14α-dihydroxy-7,22-diene-3,6-dione, 12α-hydroxyergosta-7,22,24(28)-triene-3-one, and 3β,12α-dihydroxyergosta-7,24(28)-diene, along with a known congener (22E,24R)-9α,14β-dihydroxyergosta-4,7,22-triene-3,6-dione, were isolated from the fungus Lasiodiplodia pseudotheobromae. Their structures were elucidated using NMR, HRESIMS, ECD calculation, and X-ray diffraction analyses. (22E,24R)-4α,5α-epoxyergosta-9α,14β-dihydroxy-7,22-diene-3,6-dione and (22E,24R)-4α,5α-epoxyergosta-9α,14α-dihydroxy-7,22-diene-3,6-dione are a pair of C-14 epimers possessing an unusual epoxy group between C-4 and C-5, which was demonstrated using single-crystal X-ray diffraction analyses. The absolute configurations of 12α-hydroxyergosta-7,22,24(28)-triene-3-one and 3β,12α-dihydroxyergosta-7,24(28)-diene were determined by ECD calculations. Moreover, 3β,12α-dihydroxyergosta-7,24(28)-diene exhibited neuroprotective activity in vitro in glutamate-treated SH-SY5Y cell lines.

Infection Dynamics of ATG8 in Leishmania: Balancing Autophagy for Therapeutics

In many regions of the world, Leishmaniasis is a cause of substantial mortality and ailment. Due to impediment in available treatment, development of novel and effective treatments is indispensable. Significance of autophagy has been accentuated in infectious disease as well as in Leishmaniasis, and it is having capability to be manifested as a therapeutic target. By evincing autophagy as a novel therapeutic regime, this study emphasized on the critical role of ATG4.1-ATG8 and ATG5-ATG12 complexes in Leishmania species. The objective here was to identify ATG8 as a potential therapeutic target in Leishmania. R71T, P56E, R18P are the significant mutations which shows detrimental effect on ATG8 while Arg276, Arg73, Cys75 of ATG4.1 and Val88, Pro89, Glu116, Asn117, and Gly120 are interacting residues of ATG8. Along with this, we also bring into spotlight an enticing role of Thiabendazole derivatives that interferes with the survival mechanisms by targeting ATG8. Further, the study claims that thiabendazole can be a potential drug candidate to target autophagy process in the infectious disease Leishmaniasis.

Structure based identification of first-in-class fragment inhibitors that target the NMN pocket of M. tuberculosis NAD+-dependent DNA ligase A

NAD⁺-dependent DNA ligase (LigA) is the essential replicative ligase in bacteria and differs from ATP-dependent counterparts like the human DNA ligase I (HligI) in several aspects. LigA uses NAD⁺ as the co-factor while the latter uses ATP. Further, the LigA carries out enzymatic activity with a single divalent metal ion in the active site while ATP-dependent ligases use two metal ions. Instead of the second metal ion, LigA have a unique NMN binding subdomain that facilitates the orientation of the β-phosphate and NMN leaving group. LigA are therefore attractive targets for new anti-bacterial therapeutic development. Others and our group have earlier identified several LigA inhibitors that mainly bind to AMP binding site of LigA. However, no inhibitor is known to bind to the unique NMN binding subdomain. We initiated a fragment inhibitor discovery campaign against the M. tuberculosis LigA based on our co-crystal structure of adenylation domain with AMP and NMN. The study identified two fragments, 4-(4-fluorophenyl)-4,5,6,7-tetrahydro-3H imidazo[4,5-c] pyridine and N-(4-methylbenzyl)-1H-pyrrole-2-carboxamide, that bind to the NMN site. The fragments inhibit LigA with IC₅₀ of 16.9 and 28.7 µM respectively and exhibit MIC of ~20 and 60 µg/ml against a temperature sensitive E. coli GR501 ligA^ts strain, rescued by MtbLigA. Co-crystal structures of the fragments with the adenylation domain of LigA show that they mimic the interactions of NMN. Overall, our results suggest that the NMN binding-site is a druggable target site for developing anti-LigA therapeutic strategies.

Lignans from Schisandra chinensis ameliorate alcohol and CCl4-induced long-term liver injury and reduce hepatocellular degeneration via blocking ETBR

ETHNOPHARMACOLOGICAL RELEVANCE

Chemical hepatotoxicity, especially alcoholic liver injury (ALI), commonly occurs in young and middle-aged people who drink heavily. ALI is extremely harmful and can induce severe disease states, such as hepatitis, liver fibrosis, cirrhosis, or liver cancer, which are similar to CCl₄-induced liver disease states in animals. In recent studies, the pathological changes of hepatocytes and the hepatic stellate cell have shown a significant connection between endoplasmic reticulum (ER) stress and the development of liver pathology in patients. However, the detailed pathological mechanism needs to be further studied. Schisandra chinensis, (S. chinensis), a fruit-bearing vine used in Traditional Chinese Medicine (TCM), has been used to treat chronic or acute diseases, including liver disease. S. chinensis-derived lignans (SCDLs) in particular have been shown to alleviate liver pathological changes.

AIM OF THE STUDY

This study sought to elucidate the mechanisms underlying SCDL-mediated hepatoprotection.

MATERIALS AND METHODS

We first used in silico target prediction and computational simulation methods to identify putative lignan-binding targets relative to the hepatoprotective effect. A gene microarray analysis was performed to identify differently expressed genes that might have significance in the disease pathological process. We then used histological analyses in a mice hepatotoxicity model to test the effectiveness of SCDLs in vivo, and a hepatocellular toxicity model to analyze the candidate-compound-mediated hepatoprotection and expression states of the key targets in vitro.

RESULTS

The in silico analysis results indicated that endothelin receptor B (ETBR/EDNRB) is likely a significant node during the liver pathological change process and a promising key target for the SCDL compound schisantherin D on the hepatoprotective effect; experimental studies showed that schisantherin D alleviated the EtOH- and ET-1-induced HL-7702 cell (belongs to liver parenchymal cell lines) injury ratio, decreased the expression of ETBR, and inhibited ECMs and ET-1 secretion in LX-2 cells (one form of hepatic stellate cells). SCDLs ameliorated EtOH- and CCl₄-induced fibrosis formation in mice liver tissue. Liver tissue western blots of SCDL-treated mice showed downregulated α-SMA, ETBR, PLCβ, CHOP, Bax, and the apoptotic factors of cleaved-caspase 12, cleaved-caspase 9, and cleaved-caspase 3 hinted at an anti-apoptosis and hepatoprotective effect. The SCDL treatment also elevated serum glutathione (GSH) and reduced the serum-transforming growth factor-β1 (TGF-β1) level.

CONCLUSION

The findings indicated that SCDLs prevent hepatotoxicity via their anti-fibrotic, anti-oxidant, and anti-apoptosis properties. ETBR may be the key factor in promoting chemical hepatotoxicity.

Identification of potential leukocyte antigen-related protein (PTP-LAR) inhibitors through 3D QSAR pharmacophore-based virtual screening and molecular dynamics simulation

Owing to its negative regulatory role in insulin signaling, protein tyrosine phosphatase of leukocyte antigen-related protein (PTP-LAR) was widely thought as a potential drug target for diabetes. Now, it was urgent to search for potential LAR inhibitors targeting diabetes. Initially, the pharmacophore models of LAR inhibitors were established with the application of the HypoGen module. The cost analysis, test set validation, as well as Fischer's test was used to verify the efficiency of pharmacophore model. Then, the best pharmacophore model (Hypo-1-LAR) was applied for the virtual screening of the ZINC database. And 30 compounds met the Lipinski's rule of five. Among them, 10 compounds with better binding affinity than the known LAR inhibitor (BDBM50296375) were discovered by docking studies. Finally, molecular dynamics simulations and post-analysis experiments (RMSD, RMSF, PCA, DCCM and RIN) were conducted to explore the effect of ligands (ZINC97018474 and Compound 1) on LAR and preliminary understand why ZINC97018474 had better inhibitory activity than Compound 1 (BDBM50296375). Communicated by Ramaswamy H. Sarma.

Design, synthesis and biological evaluation of pyrimidine derivatives as novel CDK2 inhibitors that induce apoptosis and cell cycle arrest in breast cancer cells

Cyclin-dependent kinase 2 (CDK2) plays a key role in eukaryotic cell cycle progression which could facilitate the transition from G1 to S phase. The dysregulation of CDK2 is closely related to many cancers. CDK2 is utilized as one of the most studied kinase targets in oncology. In this article, 24 benzamide derivatives were designed, synthesized and investigated for the inhibition activity against CDK2. Our results revealed that the compound 25 is a potent CDK2 inhibitor exhibiting a broad spectrum anti-proliferative activity against several human breast cancer cells. Additionally, compound 25 could block cell cycle at G0 or G1 and induce significant apoptosis in MDA-MB-468 cells. These findings highlight a rationale for further development of CDK2 inhibitors to treat human breast cancer.

Design, synthesis and structure-activity relationship of a focused library of β-phenylalanine derivatives as novel eEF2K inhibitors with apoptosis-inducing mechanisms in breast cancer

Eukaryotic elongation factor 2 kinase (eEF2K) is a Ca²⁺/calmudulin-dependent protein kinase, belonging to a small family of an atypical Ser/Thr-protein kinase. eEF2K has been recently reported to be highly activated or overexpressed in many types of cancer; therefore, eEF2K would be regarded as a promising therapeutic target. In this study, we discovered a β-phenylalanine scaffold by virtual high-throughput screening, as well as designed and synthesized 46 derivatives with assessment of inhibition activity against eEF2K and cytotoxicity. After several rounds of kinase and anti-proliferative activity screening, we discovered an eEF2K inhibitor (21l) with best eEF2K enzymatic activity (IC₅₀ of 5.5 μM) and anti-proliferative activity (MDA-MB-231 cells, IC₅₀ of 12.6 μM, MDA-MB-436 cells, IC₅₀ of 19.8 μM). Moreover, we found that 21l could induce cell death via the apoptotic pathways in MDA-MB-231 and MDA-MB-436 cells. Subsequently, we evaluated its anti-tumor activity and apoptosis-inducing mechanisms in vivo. These results suggested that 21l inhibited tumor growth by apoptosis in the xenograft mouse model of breast cancer (MDA-MB-231 and MDA-MB-436). Collectively, our results demonstrate a novel small-molecule inhibitor targeting eEF2K with mechanism of apoptosis and a therapeutic potential in breast cancer.

Synthesis and biological evaluation of benzimidazole derivatives as the G9a Histone Methyltransferase inhibitors that induce autophagy and apoptosis of breast cancer cells

G9a (also known as KMT1C or EHMT2) is initially identified as a H3K9 methyltransferase that specifically mono- and dimethylates 'Lys-9' of histone H3 (H3K9me1 and H3K9me2, respectively) in euchromatin. It is overexpressed in various human cancers and employed as a promising target in cancer therapy. We discovered a benzoxazole scaffold through virtual high-throughput screening, and designed, synthesized 24 derivatives and investigated for inhibition of G9a. After several rounds of kinase and anti-proliferative activity screening, we discovered a potent G9a antagonist (GA001) with an IC₅₀ value of 1.32μM that could induce autophagy via AMPK in MCF7 cells. In addition, we found high concentration of GA001 could induce apoptosis via p21-Bim signal cascades in MCF7 cells. Our results highlight a new approach for the development of a novel drug targeting G9a with a potential to induce autophagy and apoptosis for future breast cancer therapy.

Ecdysterones from Rhaponticum carthamoides (Willd.) Iljin reduce hippocampal excitotoxic cell loss and upregulate mTOR signaling in rats

Glutamate-induced excitotoxicity is a key pathological mechanism in many neurological disease states. Ecdysterones derived from Rhaponticum carthamoides (Willd.) Iljin (RCI) have been shown to alleviate glutamate-induced neuronal damage; although their mechanism of action is unclear, some data suggest that they enhance signaling in the mechanistic target of rapamycin (mTOR) signaling pathway. This study sought to elucidate the mechanisms underlying ecdysterone-mediated neuroprotection. We used in silico target prediction and simulation methods to identify putative ecdysterone binding targets, and to specifically identify those that represent nodes where several neurodegenerative diseases converge. We then used histological analyses in a rat hippocampal excitotoxicity model to test the effectiveness of ecdysterones in vivo. We found that RCI-derived ecdysterones should bind to glutamatergic NMDA-type receptors (NMDARs); specifically, in vivo modeling showed binding to the GRIN2B subunit of NMDARs, which was found also to be a node of convergence in several neurodegenerative disease pathways. Computerized network construction by using pathway information from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database showed putative links between GRIN2B and mTOR pathway elements including phosphoinositide-3kinase (PI3K), mTOR, and protein kinase C (PKC); these elements are associated with neuronal survival. Brain tissue western blots of ecdysterone-treated rats showed upregulated PI3K, Akt, mTOR, and phosphorylated Akt and mTOR, and down regulated GRIN2B and the apoptotic enzyme cleaved caspase-3. Ecdysterone treatment also prevented glutamate-induced rat hippocampal cell loss. In summary, RCI-derived ecdysterones appear to prevent glutamatergic excitotoxicity by increasing mTOR/Akt/PI3K signaling activity.

L1198F Mutation Resensitizes Crizotinib to ALK by Altering the Conformation of Inhibitor and ATP Binding Sites

The efficacy of anaplastic lymphoma kinase (ALK) positive non-small-cell lung cancer (NSCLC) treatment with small molecule inhibitors is greatly challenged by acquired resistance. A recent study reported the newest generation inhibitor resistant mutation L1198F led to the resensitization to crizotinib, which is the first Food and Drug Administration (FDA) approved drug for the treatment of ALK-positive NSCLC. It is of great importance to understand how this extremely rare event occurred for the purpose of overcoming the acquired resistance of such inhibitors. In this study, we exploited molecular dynamics (MD) simulation to dissect the molecular mechanisms. Our MD results revealed that L1198F mutation of ALK resulted in the conformational change at the inhibitor site and altered the binding affinity of ALK to crizotinib and lorlatinib. L1198F mutation also affected the autoactivation of ALK as supported by the identification of His1124 and Tyr1278 as critical amino acids involved in ATP binding and phosphorylation. Our findings are valuable for designing more specific and potent inhibitors for the treatment of ALK-positive NSCLC and other types of cancer.

Crystal structure-based discovery of a novel synthesized PARP1 inhibitor (OL-1) with apoptosis-inducing mechanisms in triple-negative breast cancer

Poly (ADP-ribose) polymerase-1 (PARP1) is a highly conserved enzyme focused on the self-repair of cellular DNA damage. Until now, numbers of PARP inhibitors have been reported and used for breast cancer therapy in recent years, especially in TNBC. However, developing a new type PARP inhibitor with distinctive skeleton is alternatively promising strategy for TNBC therapy. In this study, based on co-crystallization studies and pharmacophore-docking-based virtual screening, we discovered a series of dihydrodibenzo[b,e]-oxepin compounds as PARP1 inhibitors. Lead optimization result in the identification of compound OL-1 (2-(11-(3-(dimethylamino)propylidene)-6,11- dihydrodibenzo[b,e]oxepin )-2-yl)acetohydrazide), which has a novel chemical scaffold and unique binding interaction with PARP1 protein. OL-1 demonstrated excellent potency (inhibiting PARP1 enzyme activity with IC₅₀ = 0.079 μM), as well as inhibiting PARP-modulated PARylation and cell proliferation in MDA-MB-436 cells (BRAC1 mutation). In addition, OL-1 also inhibited cell migration that closely related to cancer metastasis and displayed remarkable anti-tumor efficacy in MDA-MB-436 xenograft model without apparent toxicities. These findings highlight a new small-molecule PAPR1 inhibitor (OL-1) that has the potential to impact future TNBC therapy.

Ursolic acid derivatives as bone anabolic agents targeted to tryptophan hydroxylase 1 (Tph-1)

Tryptophan hydroxylase 1 (Tph-1) initiates the biosynthesis of peripheral serotonin. As peripheral serotonin suppresses bone formation, inhibitor of Tph-1 provides a useful tool to discover anabolic agents for osteoporosis. In the present study, series of ursolic acid (UA) derivatives were synthesized, and their inhibitory activity on serotonin biosynthesis and cytotoxicity were evaluated. Among the derivatives, 8d with potent inhibitory activity on serotonin was applied for further research. The data revealed that 8d significantly inhibited protein and mRNA expressions of Tph-1, and an SPR study indicated that 8d directly interacted to Tph-1 with a binding affinity of KD=15.09μM. Oral administration of 8d significantly prevented bone loss via suppressing serotonin biosynthesis without estrogenic side-effects in ovariectomized (OVX) rats.

Systems biology network-based discovery of a small molecule activator BL-AD008 targeting AMPK/ZIPK and inducing apoptosis in cervical cancer

The aim of this study was to discover a small molecule activator BL-AD008 targeting AMPK/ZIPK and inducing apoptosis in cervical cancer. In this study, we systematically constructed the global protein-protein interaction (PPI) network and predicted apoptosis-related protein connections by the Naïve Bayesian model. Then, we identified some classical apoptotic PPIs and other previously unrecognized PPIs between apoptotic kinases, such as AMPK and ZIPK. Subsequently, we screened a series of candidate compounds targeting AMPK/ZIPK, synthesized some compounds and eventually discovered a novel dual-target activator (BL-AD008). Moreover, we found BL-AD008 bear remarkable anti-proliferative activities toward cervical cancer cells and could induce apoptosis by death-receptor and mitochondrial pathways. Additionally, we found that BL-AD008-induced apoptosis was affected by the combination of AMPK and ZIPK. Then, we found that BL-AD008 bear its anti-tumor activities and induced apoptosis by targeting AMPK/ZIPK in vivo. In conclusion, these results demonstrate the ability of systems biology network to identify some key apoptotic kinase targets AMPK and ZIPK; thus providing a dual-target small molecule activator (BL-AD008) as a potential new apoptosis-modulating drug in future cervical cancer therapy.

Patent Highlights

A snapshot of noteworthy recent developments in the patent literature of relevance to pharmaceutical and medical research and development.

In silico discovery of potential VEGFR-2 inhibitors from natural derivatives for anti-angiogenesis therapy

Angiogenesis is the growth of new capillaries from existing blood vessels that supply oxygen and nutrients and provide gateways for immune surveillance. Abnormal vessel growth in term of excessive angiogenesis is a hallmark of cancer, inflammatory and eye diseases. VEGFR-2 (vascular endothelial growth factor receptor 2) dominating the process of angiogenesis has led to approval of therapeutic inhibitors and is becoming a promising target for anti-angiogenic drugs. Notwithstanding these successes, the clinical use of current VEGFR-2 blockers is more challenging than anticipated. Taking axitinib as a reference drug, in our study we found three potent VEGFR-2 inhibitors (ZINC08254217, ZINC08254138, and ZINC03838680) from natural derivatives. Each of the three inhibitors acquired a better grid score than axitinib (−62.11) when docked to VEGFR-2. Molecular dynamics simulations demonstrated that ZINC08254217– and ZINC08254138–VEGFR-2 complexes were more stable than axitinib. Similar to bind free energy for axitinib (−54.68 kcal/mol), such for ZINC03838680, ZINC08254217, and ZINC08254138 was −49.37, −43.32, and −32.73 kcal/mol respectively. These results suggested these three compounds could be candidate drugs against angiogenesis, with comparable VEGFR-2 binding affinity of axitinib. Hence findings in our study are able to provide valuable information on discovery of effective anti-angiogenesis therapy.

Combining structure-based pharmacophore modeling, virtual screening, and in silico ADMET analysis to discover novel tetrahydro-quinoline based pyruvate kinase isozyme M2 activators with antitumor activity

Compared with normal differentiated cells, cancer cells upregulate the expression of pyruvate kinase isozyme M2 (PKM2) to support glycolytic intermediates for anabolic processes, including the synthesis of nucleic acids, amino acids, and lipids. In this study, a combination of the structure-based pharmacophore modeling and a hybrid protocol of virtual screening methods comprised of pharmacophore model-based virtual screening, docking-based virtual screening, and in silico ADMET (absorption, distribution, metabolism, excretion and toxicity) analysis were used to retrieve novel PKM2 activators from commercially available chemical databases. Tetrahydroquinoline derivatives were identified as potential scaffolds of PKM2 activators. Thus, the hybrid virtual screening approach was applied to screen the focused tetrahydroquinoline derivatives embedded in the ZINC database. Six hit compounds were selected from the final hits and experimental studies were then performed. Compound 8 displayed a potent inhibitory effect on human lung cancer cells. Following treatment with Compound 8, cell viability, apoptosis, and reactive oxygen species (ROS) production were examined in A549 cells. Finally, we evaluated the effects of Compound 8 on mice xenograft tumor models in vivo. These results may provide important information for further research on novel PKM2 activators as antitumor agents.

Tryptophan hydroxylase 1 (Tph-1)-targeted bone anabolic agents for osteoporosis

Tryptophan hydroxylase 1 (Tph-1), the principal enzyme for peripheral serotonin biosynthesis, provides a novel target to design anabolic agents for osteoporosis. Here, we present a design, synthesis of a novel series of ursolic acid derivatives under the guidance of docking technique, and bioevaluation of the derivatives using RBL2H3 cells and ovariectomized (OVX) rats. Of the compounds, 9a showed a potent inhibitory activity on serotonin biosynthesis. Further investigations revealed that 9a, as an efficient Tph-1 binder identified by SPR (estimated KD: 6.82 μM), suppressed the protein and mRNA expressions of Tph-1 and lowered serotonin contents in serum and gut without influence on brain serotonin. Moreover, oral administration of 9a elevated serum level of N-terminal propeptide of procollagen type 1 (P1NP), a bone formation marker, and improved bone microarchitecture without estrogenic side effects in ovariectomized rats. Collectively, 9a may serve as a new candidate for bone anabolic drug discovery.

Fast, metadynamics-based method for prediction of the stereochemistry-dependent relative free energies of ligand-receptor interactions

The computational approach applicable for the molecular dynamics (MD)-based techniques is proposed to predict the ligand-protein binding affinities dependent on the ligand stereochemistry. All possible stereoconfigurations are expressed in terms of one set of force-field parameters [stereoconfiguration-independent potential (SIP)], which allows for calculating all relative free energies by only single simulation. SIP can be used for studying diverse, stereoconfiguration-dependent phenomena by means of various computational techniques of enhanced sampling. The method has been successfully tested on the β2-adrenergic receptor (β2-AR) binding the four fenoterol stereoisomers by both metadynamics simulations and replica-exchange MD. Both the methods gave very similar results, fully confirming the presence of stereoselective effects in the fenoterol-β2-AR interactions. However, the metadynamics-based approach offered much better efficiency of sampling which allows for significant reduction of the unphysical region in SIP.

Novel hybrid virtual screening protocol based on molecular docking and structure-based pharmacophore for discovery of methionyl-tRNA synthetase inhibitors as antibacterial agents

Methione tRNA synthetase (MetRS) is an essential enzyme involved in protein biosynthesis in all living organisms and is a potential antibacterial target. In the current study, the structure-based pharmacophore (SBP)-guided method has been suggested to generate a comprehensive pharmacophore of MetRS based on fourteen crystal structures of MetRS-inhibitor complexes. In this investigation, a hybrid protocol of a virtual screening method, comprised of pharmacophore model-based virtual screening (PBVS), rigid and flexible docking-based virtual screenings (DBVS), is used for retrieving new MetRS inhibitors from commercially available chemical databases. This hybrid virtual screening approach was then applied to screen the Specs (202,408 compounds) database, a structurally diverse chemical database. Fifteen hit compounds were selected from the final hits and shifted to experimental studies. These results may provide important information for further research of novel MetRS inhibitors as antibacterial agents.