Structure-based Discovery Research Articles

Structure-based discovery of first inhibitors targeting the helicase activity of human PIF1.

PIF1 is a conserved helicase and G4 DNA binding and unwinding enzyme, with roles in genome stability. Human PIF1 (hPIF1) is poorly understood, but its functions can become critical for tumour cell survival during oncogene-driven replication stress. Here we report the discovery, via an X-ray crystallographic fragment screen (XChem), of hPIF1 DNA binding and unwinding inhibitors. A structure was obtained with a 4-phenylthiazol-2-amine fragment bound in a pocket between helicase domains 2A and 2B, with additional contacts to Valine 258 from domain 1A. The compound makes specific interactions, notably through Leucine 548 and Alanine 551, that constrain conformational adjustments between domains 2A and 2B, previously linked to ATP hydrolysis and DNA unwinding. We next synthesized a range of related compounds and characterized their effects on hPIF1 DNA-binding and helicase activity in vitro, expanding the structure activity relationship (SAR) around the initial hit. A systematic analysis of clinical cancer databases is also presented here, supporting the notion that hPIF1 upregulation may represent a specific cancer cell vulnerability. The research demonstrates that hPIF1 is a tractable target through 4-phenylthiazol-2-amine derivatives as inhibitors of its helicase action, setting a foundation for creation of a novel class of anti-cancer therapeutics.

Design of a Cereblon construct for crystallographic and biophysical studies of protein degraders

The ubiquitin E3 ligase cereblon (CRBN) is the target of therapeutic drugs thalidomide and lenalidomide and is recruited by most targeted protein degraders (PROTACs and molecular glues) in clinical development. Biophysical and structural investigation of CRBN has been limited by current constructs that either require co-expression with the adaptor DDB1 or inadequately represent full-length protein, with high-resolution structures of degrader ternary complexes remaining rare. We present the design of CRBNmidi, a construct that readily expresses from E. coli with high yields as soluble, stable protein without DDB1. We benchmark CRBNmidi for wild-type functionality through a suite of biophysical techniques and solve high-resolution co-crystal structures of its binary and ternary complexes with degraders. We qualify CRBNmidi as an enabling tool to accelerate structure-based discovery of the next generation of CRBN based therapeutics.

Structure-Based Discovery of Phytocompounds from Azadirachta indica as Potential Inhibitors of Thioredoxin Glutathione Reductase in Schistosoma mansoni.

Schistosomiasis, a parasitic disease caused by Schistosoma species such as S. haematobium, S. mansoni, and S. japonicum, poses a significant global health burden. The thioredoxin glutathione reductase (TGR) enzyme, crucial for maintaining the parasite's redox balance and preventing oxidative stress, has been identified as a promising target for anti-schistosomal drug development. This study aims to identify potential TGR inhibitors from Azadirachta indica phytochemicals using molecular modeling approaches. We screened 60 compounds derived from A. indica bark and leaves through molecular docking to assess their binding affinity, followed by the evaluation of binding-free energies for the most promising candidates. Drug-likeness and pharmacokinetic properties were assessed, and molecular dynamics simulations were conducted to explore the conformational stability of the protein-ligand complexes. Our findings revealed that several A. indica compounds exhibited significantly lower docking scores (up to -9.669 kcal/mol) compared to the standard drug praziquantel (-4.349 kcal/mol). Notably, Isorhamnetin, Isomargolonone, Nimbaflavone, Quercetin, and Nimbionol demonstrated strong interactions with TGR, although Isorhamnetin showed potential mutagenicity. Further binding free energy calculations and molecular dynamics simulations confirmed the stability of Isomargolonone, Nimbionol, and Quercetin as potential TGR inhibitors. In conclusion, these findings suggest that Isomargolonone, Nimbionol, and Quercetin warrant further experimental validation as promising candidates for anti-schistosomal therapy.

Discovery and optimization of 4-(imidazo[1,2-a]pyrimidin-3-yl)thiazol-2-amine derivatives as novel phosphodiesterase 4 inhibitors.

Phosphodiesterases (PDEs) are important intracellular enzymes that hydrolyze the second messengers cAMP and/or cGMP. Now several studies have shown that PDE4 received particular attention due to which it represents the most prominent cAMP-metabolizing enzyme involved in many diseases. In this study, we performed prescreening of our internal compound library and discovered the compound (PTC-209) with moderate PDE4 inhibitory activity (IC50 of 4.78 ± 0.08μM). And a series of 4-(imidazo[1,2-a]pyrimidin-3-yl)thiazol-2-amine derivatives as novel PDE4 inhibitors starting from PTC-209 were successfully designed and synthesized using a structure-based discovery strategy. L19, the most potent inhibitor, exhibited good inhibitory activity (IC50 of 0.48 ± 0.02μM) and remarkable metabolic stability in rat liver microsomes. Our study presents an example of discovery novel PDE4 inhibitors, which would be helpful for design and optimization of novel inhibitors in future.

Structure-based discovery of a 4,5-Dihydropyrazole-cored PET ligand for imaging of receptor-interacting serine/threonine-protein kinase 1 (RIPK1) in the brain

Receptor-interacting serine/threonine-protein kinase 1 (RIPK1) regulates programmed cell death and inflammation, contributing to a wide range of human pathologies, including inflammatory disorders, neurodegenerative conditions, and cancer. Despite this, no RIPK1 positron emission tomography (PET) ligand with significant in vivo specificity has been reported to date. In this work, we designed and synthesized a new family of dihydropyrazole-cored ligands suitable for 18F-labeling at the late stage. Among these, WL8 showed a strong binding affinity to RIPK1 (EC50 = 19.9 nM, Kd = 25 nM) and was successfully labeled with 18F in the 6-position of pyridine ring, yielding a high radiochemistry yield of 27.9 % (decay-corrected) and a high molar activity of 18.8–31.2 GBq/μmol. In in vitro autoradiography, [18F]WL8 showed some specific binding in the brain sections of rats and lipopolysaccharide (LPS) model mice. Preliminary PET studies in rat brains revealed that [18F]WL8 could efficiently penetrate the blood-brain barrier and was rapidly washed out. As anticipated, [18F]WL8 exhibited a high initial uptake (brain2min = 4.80 % ID/g) in mouse brains, followed by a rapid washout (brain60min = 0.14 % ID/g), although no clear specific binding to RIPK1 was observed. Moderate in vivo stability was noted for [18F]WL8 in mouse brains with 35.2 % of the parent fraction remaining after 30 min post-administration. Altogether, our work broadens the landscape and offers a new chemotype for RIPK1 PET ligand development.

Structure-Based Design and Discovery of a Potent and Cell-Active LC3A/B Covalent Inhibitor.

Autophagy is a highly conserved cellular homeostasis maintenance mechanism in eukaryotes. Microtubule-associated protein light chain 3 (LC3) plays a crucial role in autophagy. It has multiple pairs of protein-protein interactions (PPIs) with other proteins, and these PPIs have an effect on the regulation of autophagosome formation and the recruitment of autophagic substrates. In our previous work, a small molecule covalent inhibitor DC-LC3in-D5 which could inhibit LC3A/B PPIs was identified, but a detailed study of structure-activity relationships (SARs) was lacking. Herein, a new molecule LC3in-C42 was discovered utilizing the hybridization of advantageous fragments, whose potency (IC50 = 7.6 nM) had been greatly improved compared with that of DC-LC3in-D5. LC3in-C42 inhibits autophagy at the cellular level and its efficacy far exceeds that of DC-LC3in-D5. Thus far, LC3in-C42 stands as the most potent LC3A/B small molecule inhibitor. LC3in-C42 could serve as a powerful tool for LC3A/B protein and autophagy research.

Structure-based discovery of novel piperidine-biphenyl-DAPY derivatives as non-nucleoside reverse transcriptase inhibitors featuring improved potency, safety, and selectivity: From piperazine-biphenyl-DAPYs to piperidine-biphenyl-DAPYs

Starting from our previously reported nonnucleoside reverse transcriptase inhibitor (NNRTI, 3), continuous efforts were made to enhance its potency and safety through a structure-based drug design strategy. This led to the discovery of a series of novel piperidine-biphenyl-diarylpyrimidines (DAPYs). Compound 10p, the most active compound in this series, exhibited an EC50 value of 6 nM against wide-type HIV-1 strain, which was approximately 560-fold more potent than the initial compound 3 (EC50 = 3.36 μM). Furthermore, significant improvements were observed in cytotoxicity and selectivity (CC50 > 202.17 μM, SI > 33144) compared to compound 3 (CC50 = 14.84 μM, SI = 4). Additionally, compound 10p demonstrated increased inhibitory activity against clinically mutant virus strains (EC50 = 7–63 nM). Further toxicity evaluation revealed that compound 10p exhibited minimal CYP enzyme and hERG inhibition. Importantly, single-dose acute toxicity testing did not result in any fatalities or noticeable pathological damage in mice. Therefore, compound 10p can be regarded as a lead candidate for guiding further development of biphenyl-diarylpyrimidine NNRTIs with favorable druggability for HIV therapy.

Structure-based discovery of F. religiosa phytochemicals as potential inhibitors against Monkeypox (mpox) viral protein

Outbreaks of Monkeypox (mpox) in over 100 non-endemic countries in 2022 represented a serious global health concern. Once a neglected disease, mpox has become a global public health issue. A42R profilin-like protein from mpox (PDB ID: 4QWO) represents a potential new lead for drug development and may interact with various synthetic and natural compounds. In this report, the interaction of A42R profilin-like protein with six phytochemicals found in the medicinal plant Ficus religiosa (abundant in India) was examined. Based on the predicted and compared protein–ligand binding energies, biological properties, IC50 values and toxicity, two compounds, kaempferol (C-1) and piperine (C-4), were selected. ADMET characteristics and quantitative structure–activity relationship (QSAR) of these two compounds were determined, and molecular dynamics (MD) simulations were performed. In silico examination of the kaempferol (C-1) and piperine (C-4) interactions with A42R profilin-like protein gave best-pose ligand-binding energies of –6.98 and –5.57 kcal/mol, respectively. The predicted IC50 of C-1 was 7.63 μM and 82 μM for C-4. Toxicity data indicated that kaempferol and piperine are non-mutagenic, and the QSAR data revealed that piperlongumine (5.92) and piperine (5.25) had higher log P values than the other compounds examined. MD simulations of A42R profilin-like protein in complex with C-1 and C-4 were performed to examine the stability of the ligand–protein interactions. As/C and C-4 showed the highest affinity and activities, they may be suitable lead candidates for developing mpox therapeutic drugs. This study should facilitate discovering and synthesizing innovative therapeutics to address other infectious diseases.

AlphaFold2 structures guide prospective ligand discovery.

AlphaFold2 (AF2) models have had wide impact but mixed success in retrospective ligand recognition. We prospectively docked large libraries against unrefined AF2 models of the σ2 and serotonin 2A (5-HT2A) receptors, testing hundreds of new molecules and comparing results with those obtained from docking against the experimental structures. Hit rates were high and similar for the experimental and AF2 structures, as were affinities. Success in docking against the AF2 models was achieved despite differences between orthosteric residue conformations in the AF2 models and the experimental structures. Determination of the cryo-electron microscopy structure for one of the more potent 5-HT2A ligands from the AF2 docking revealed residue accommodations that resembled the AF2 prediction. AF2 models may sample conformations that differ from experimental structures but remain low energy and relevant for ligand discovery, extending the domain of structure-based drug design.

Structure-Based Design, Virtual Screening, and Discovery of Novel Patulin Derivatives as Biogenic Photosystem II Inhibiting Herbicides.

Computer-aided design usually gives inspirations and has become a vital strategy to develop novel pesticides through reconstructing natural lead compounds. Patulin, an unsaturated heterocyclic lactone mycotoxin, is a new natural PSII inhibitor and shows significant herbicidal activity to various weeds. However, some evidence, especially the health concern, prevents it from developing as a bioherbicide. In this work, molecular docking and toxicity risk prediction are combined to construct interaction models between the ligand and acceptor, and design and screen novel derivatives. Based on the analysis of a constructed patulin-Arabidopsis D1 protein docking model, in total, 81 derivatives are designed and ranked according to quantitative estimates of drug-likeness (QED) values and free energies. Among the newly designed derivatives, forty-five derivatives with better affinities than patulin are screened to further evaluate their toxicology. Finally, it is indicated that four patulin derivatives, D3, D6, D34, and D67, with higher binding affinity but lower toxicity than patulin have a great potential to develop as new herbicides with improved potency.

Macrocyclic Azapeptide Nitriles: Structure-Based Discovery of Potent SARS-CoV-2 Main Protease Inhibitors as Antiviral Drugs.

Given the crucial role of the main protease (Mpro) in the replication cycle of SARS-CoV-2, this viral cysteine protease constitutes a high-profile drug target. We investigated peptidomimetic azapeptide nitriles as auspicious, irreversibly acting inhibitors of Mpro. Our systematic approach combined an Mpro active-site scanning by combinatorially assembled azanitriles with structure-based design. Encouraged by the bioactive conformation of open-chain inhibitors, we conceptualized the novel chemotype of macrocyclic azanitriles whose binding mode was elucidated by cocrystallization. This strategy provided a favorable entropic contribution to target binding and resulted in the development of the extraordinarily potent Mpro inhibitor 84 with an IC50 value of 3.23 nM and a second-order rate constant of inactivation, kinac/Ki, of 448,000 M-1s-1. The open-chain Mpro inhibitor 58, along with the macrocyclic compounds 83 and 84, a broad-spectrum anticoronaviral agent, demonstrated the highest antiviral activity with EC50 values in the single-digit micromolar range. Our findings are expected to promote the future development of peptidomimetic Mpro inhibitors as anti-SARS-CoV-2 agents.

Structure-Based Discovery of Potential HPV E6 and EBNA1 Inhibitors: Implications for Cervical Cancer Treatment

Cervical cancer is the fourth most diagnosed cancer and the fourth leading cause of cancer death in women globally. Its onset and progression have been attributed to high-risk human papillomavirus (HPV) types, especially 16 and 18, while the Epstein–Barr virus (EBV) is believed to also significantly contribute to cervical cancer growth. The E6 protein associated with high-risk HPV strains, such as HPV16 and HPV18, is known for its role in promoting cervical cancer and other anogenital cancers. E6 proteins contribute to the malignant transformation of infected cells by targeting and degrading tumor suppressor proteins, especially p53. On the other hand, EBV nuclear antigen 1 (EBNA1) plays a crucial role in the maintenance and replication of the EBV genome in infected cells. EBNA1 is believed to increase HPV E6 and E7 levels, as well as c-MYC, and BIRC5 cellular genes in the HeLa cell line, implying that HPV/EBV co-infection accelerates cervical cancer onset and growth. Thus, the E6 and EBNA1 antigens of HPV and EBV, respectively, are attractive targets for cervical cancer immunotherapy. This study, therefore, virtually screened for potential drug candidates with good binding affinity to all three oncoviral proteins, HPV16 E6, HPV18 E6, and EBNA1. The compounds were further subjected to ADMET profiling, biological activity predictions, molecular dynamics (MD) simulations, and molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) calculations. A total of six compounds comprising ZINC000013380012, ZINC000070454124, ZINC000014588133, ZINC000085568136, ZINC000095909247, and ZINC000085597263 demonstrated very strong affinity (≤−60 kJ/mol) to the three oncoviral proteins (EBNA1, HPV16 E6, and HPV18 E6) after being subjected to docking, MD, and MM/PBSA. These compounds demonstrated relatively stronger binding than the controls used, inhibitors of EBNA1 (VK-1727) and HPV E6 (baicalein and gossypetin). Biological activity predictions also corroborated their antineoplastic, p53-enhancing, Pin1 inhibitory, and JAK2 inhibitory activities. Further experimental testing is required to validate the ability of the shortlisted compounds to silence the insidious effects of HPV E6 and EBNA1 proteins in cervical cancers.

Structure-based discovery of CFTR potentiators and inhibitors

The cystic fibrosis transmembrane conductance regulator (CFTR) is a crucial ion channel whose loss of function leads to cystic fibrosis, whereas its hyperactivation leads to secretory diarrhea. Small molecules that improve CFTR folding (correctors) or function (potentiators) are clinically available. However, the only potentiator, ivacaftor, has suboptimal pharmacokinetics and inhibitors have yet to be clinically developed. Here, we combine molecular docking, electrophysiology, cryo-EM, and medicinal chemistry to identify CFTR modulators. We docked ∼155 million molecules into the potentiator site on CFTR, synthesized 53 test ligands, and used structure-based optimization to identify candidate modulators. This approach uncovered mid-nanomolar potentiators, as well as inhibitors, that bind to the same allosteric site. These molecules represent potential leads for the development of more effective drugs for cystic fibrosis and secretory diarrhea, demonstrating the feasibility of large-scale docking for ion channel drug discovery.

Repurposing DrugBank compounds as potential Plasmodium falciparum class 1a aminoacyl tRNA synthetase multi-stage pan-inhibitors with a specific focus on mitomycin

Plasmodium falciparum aminoacyl tRNA synthetases (PfaaRSs) are potent antimalarial targets essential for proteome fidelity and overall parasite survival in every stage of the parasite's life cycle. So far, some of these proteins have been singly targeted yielding inhibitor compounds that have been limited by incidences of resistance which can be overcome via pan-inhibition strategies. Hence, herein, for the first time, we report the identification and in vitro antiplasmodial validation of Mitomycin (MMC) as a probable pan-inhibitor of class 1a (arginyl(A)-, cysteinyl(C), isoleucyl(I)-, leucyl(L), methionyl(M), and valyl(V)-) PfaaRSs which hypothetically may underlie its previously reported activity on the ribosomal RNA to inhibit protein translation and biosynthesis. We combined multiple in silico structure-based discovery strategies that first helped identify functional and druggable sites that were preferentially targeted by the compound in each of the plasmodial proteins: Ins1-Ins2 domain in Pf-ARS; anticodon binding domain in Pf-CRS; CP1-editing domain in Pf-IRS and Pf-MRS; C-terminal domain in Pf-LRS; and CP-core region in Pf-VRS. Molecular dynamics studies further revealed that MMC allosterically induced changes in the global structures of each protein. Likewise, prominent structural perturbations were caused by the compound across the functional domains of the proteins. More so, MMC induced systematic alterations in the binding of the catalytic nucleotide and amino acid substrates which culminated in the loss of key interactions with key active site residues and ultimate reduction in the nucleotide-binding affinities across all proteins, as deduced from the binding energy calculations. These altogether confirmed that MMC uniformly disrupted the structure of the target proteins and essential substrates. Further, MMC demonstrated IC50 < 5 μM against the Dd2 and 3D7 strains of parasite making it a good starting point for malarial drug development. We believe that findings from our study will be important in the current search for highly effective multi-stage antimalarial drugs.

Structure-Based Discovery of Allosteric Inhibitors Targeting a New Druggable Site in the Respiratory Syncytial Virus Polymerase.

Respiratory syncytial virus (RSV) is a major cause of severe lower respiratory infections for which effective treatment options remain limited. Herein, we employed a computational structure-based design strategy aimed at identifying potential targets for a new class of allosteric inhibitors. Our investigation led to the discovery of a previously undisclosed allosteric binding site within the RSV polymerase, the large (L) protein. This discovery was achieved through a combination of virtual screening and molecular dynamics simulations. Subsequently, we identified two inhibitors, 6a and 10b, which both exhibited promising antiviral activity in the low micromolar range. Resistance profiling revealed a distinctive pattern in how RSV evaded treatment with this class of inhibitors. This pattern strongly suggested that this class of small molecules was targeting a new binding site in the RSV L protein, aligning with the computational predictions made in our study. This study paves the way for the development of more potent inhibitors for combating RSV infections by targeting a new druggable pocket within the RdRp which does not overlap with previously known resistance sites.

Structure-based discovery of potent myosin inhibitors to guide antiparasite drug development

Monogenea, a prevalent parasite in aquaculture, poses significant threats to the industry, leading to substantial losses. Current preventive measures have proven insufficient, necessitating the development of novel and effective anti-parasitic drugs. In this investigation, we obtained the full-length myosin cDNA sequence by analyzing three-generation transcriptome data, revealing a 5817-base sequence encoding 1938 amino acids. Subsequently, we modeled and analyzed the characteristics of the secondary and tertiary of myosin, pinpointing the crucial functional region within the motor domain (amino acids 1–768). The prokaryotic expression of this domain yielded a protein of 87.44 kDa, confirmed as myosin by Western Blotting. Molecular docking identified ASN439 as the key amino acid residue involved in arctigenin and myosin binding, a result corroborated by site-directed mutagenesis, affirming the active cavity of this interaction. Chalcone and shikonin were chosen from a virtual sieve of molecular library of natural drugs based on the active cavity. Chalcone and shikonin exhibited EC50 values of 1.085 mg/L and 0.371 mg/L, respectively, with corresponding IC50 values for myosin of 0.44 mM and 0.14 mM. Given its superior activity and structure, shikonin was selected for further optimization of drug molecule design, culminating in the discovery of 1,4-naphthoquinone as a potent antiparasitic agent. This compound demonstrated an EC50 of 0.047 mg/L, LC50 of 0.23 mg/L, and a TI index of 4.893. These findings collectively highlight the potential of shikonin and 1,4-naphthoquinone as alternative compounds to control Gyrodactylus infections. Further optimization of medicinal chemistry holds promise for the development of more potent 1,4-naphthoquinone analogues, offering prospects for future anthelmintic control through combinatorial or replacement strategies.

Structure-based discovery of potent CARM1 inhibitors for colorectal cancer therapy

Coactivator-associated arginine methyltransferase 1 (CARM1) plays an important role in cell proliferation and gene expression, and is highly expressed in a variety of tumor tissues. Guided by our previous reported structure of DCPR049_12, we focused on designing and evaluating selective CARM1 inhibitors, resulting in the identification of compound 11f as a promising lead candidate. Compound 11f displayed potent inhibition of CARM1 (IC50 = 9 nM). Comprehensive evaluations, including in vitro metabolic stability assessments, molecular modelling, cellular studies, and in vivo anti-tumor studies, confirmed that it induced cancer cell apoptosis and specifically inhibited CARM1's methylation function. Notably, compound 11f displayed significant anti-proliferative effects on colorectal cancer cell lines, showcasing its potential for targeted therapies against CARM1-related diseases. This study provides valuable insights for the future development of specific and effective CARM1 inhibitors.

Structure-Based Discovery of High-Affinity Small Molecule Ligands and Development of Tool Probes to Study the Role of Chitinase-3-Like Protein 1.

Chitinase-3-like-1 (CHI3L1), also known as YKL-40, is a glycoprotein linked to inflammation, fibrosis, and cancer. This study explored CHI3L1's interactions with various oligosaccharides using microscale thermophoresis (MST) and AlphaScreen (AS). These investigations guided the development of high-throughput screening assays to assess interference of small molecules in binding between CHI3L1 and biotinylated small molecules or heparan sulfate-based probes. Small molecule binders of YKL-40 were identified in our chitotriosidase inhibitors library with MST and confirmed through X-ray crystallography. Based on cocrystal structures of potent hit compounds with CHI3L1, small molecule probes 19 and 20 were designed for an AS assay. Structure-based optimization led to compounds 30 and 31 with nanomolar activities and drug-like properties. Additionally, an orthogonal AS assay using biotinylated heparan sulfate as a probe was developed. The compounds' affinity showed a significant correlation in both assays. These screening tools and compounds offer novel avenues for investigating the role of CHI3L1.

Ligand and structure-based discovery of phosphorus-containing compounds as potential metalloproteinase inhibitors

ABSTRACT In this study, a methodology is proposed, combining ligand- and structure-based virtual screening tools, for the identification of phosphorus-containing compounds as inhibitors of zinc metalloproteases. First, we use Dragon molecular descriptors to develop a Linear Discriminant Analysis classification model, which is widely validated according to the OECD principles. This model is simple, robust, stable and has good discriminating power. Furthermore, it has a defined applicability domain and it is used for virtual screening of the DrugBank database. Second, docking experiments are carried out on the identified compounds that showed good binding energies to the enzyme thermolysin. Considering the potential toxicity of phosphorus-containing compounds, their toxicological profile is evaluated according to Protox II. Of the five molecules evaluated, two show carcinogenic and mutagenic potential at small LD50, not recommended as drugs, while three of them are classified as non-toxic, and could constitute a starting point for the development of new vasoactive metalloprotease inhibitor drugs. According to molecular dynamics simulation, two of them show stable interactions with the active site maintaining coordination with the metal. A high agreement is evident between QSAR, docking and molecular dynamics results, demonstrating the potentialities of the combination of these tools.

Structure-Based Discovery of the SARS-CoV-2 Main Protease Noncovalent Inhibitors from Traditional Chinese Medicine.

Traditional Chinese medicine (TCM) has been extensively employed for the treatment of coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, there is demand for discovering more SARS-CoV-2 Mpro inhibitors with diverse scaffolds to optimize anti-SARS-CoV-2 lead compounds. In this study, comprehensive in silico and in vitro assays were utilized to determine the potential inhibitors from TCM compounds against SARS-CoV-2 Mpro, which is an important therapeutic target for SARS-CoV-2. The ensemble docking analysis of 18263 TCM compounds against 15 SARS-CoV-2 Mpro conformations identified 19 TCM compounds as promising candidates. Further in vitro testing validated three compounds as inhibitors of SARS-CoV-2 Mpro and showed IC50 values of 4.64 ± 0.11, 7.56 ± 0.78, and 11.16 ± 0.26 μM, with EC50 values of 12.25 ± 1.68, 15.58 ± 0.77, and 29.32 ± 1.25 μM, respectively. Molecular dynamics (MD) simulations indicated that the three complexes remained stable over the last 100 ns of production run. An analysis of the binding mode revealed that the active compounds occupy different subsites (S1, S2, S3, and S4) of the active site of SARS-CoV-2 Mpro via specific poses through noncovalent interactions with key amino acids (e.g., HIS 41, ASN 142, GLY 143, MET 165, GLU 166, or GLN 189). Overall, this study provides evidence indicating that the three natural products obtained from TCM could be further used for anti-COVID-19 research, justifying the investigation of Chinese herbal medicinal ingredients as bioactive constituents for therapeutic targets.