Kinetic comparison of all eleven viral polyprotein cleavage site processing events by SARS-CoV-2 main protease using a linked protein FRET platform.

The main protease (Mpro) remains an essential therapeutic target for COVID-19 post infection intervention given its critical role in processing the majority of viral proteins encoded by the genome of severe acute respiratory syndrome related coronavirus 2 (SARS-CoV-2). Upon viral entry, the +ssRNA genome is translated into two long polyproteins (pp1a or the frameshift-dependent pp1ab) containing all the nonstructural proteins (nsps) required by the virus for immune modulation, replication, and ultimately, virion assembly. Included among these nsps is the cysteine protease Mpro (nsp5) which self-excises from the polyprotein, dimerizes, then sequentially cleaves 11 of the 15 cut-site junctions found between each nsp within the polyprotein. Many structures of Mpro (often bound to various small molecule inhibitors or peptides) have been detailed recently, including structures of Mpro bound to each of the polyprotein cleavage sequences, showing that Mpro can accommodate a wide range of targets within its active site. However, to date, kinetic characterization of the interaction of Mpro with each of its native cleavage sequences remains incomplete. Here, we present a robust and cost-effective FRET based system that benefits from a more consistent presentation of the substrate that is also closer in organization to the native polyprotein environment compared to previously reported FRET systems that use chemically modified peptides. Using this system, we were able to show that while each site maintains a similar Michaelis constant, the catalytic efficiency of Mpro varies greatly between cut-site sequences, suggesting a clear preference for the order of nsp processing.

Structural basis of rosmarinic acid inhibitory mechanism on SARS-CoV-2 main protease.

The SARS-CoV-2 coronavirus is characterized by high mutation rates and significant infectivity, posing ongoing challenges for therapeutic intervention. To address potential challenges in the future, the continued development of effective drugs targeting SARS-CoV-2 remains an important task for the scientific as well as the pharmaceutical community. The main protease (Mpro) of SARS-CoV-2 is an ideal therapeutic target for COVID-19 drug development, leading to the introduction of various inhibitors, both covalent and non-covalent, each characterized by unique mechanisms of action and possessing inherent strengths and limitations. Natural products, being compounds naturally present in the environment, offer advantages such as low toxicity and diverse activities, presenting a viable source for antiviral drug development. Here, we identified a natural compound, rosmarinic acid, which exhibits significant inhibitory effects on the Mpro of the SARS-CoV-2. Through detailed structural biology analysis, we elucidated the precise crystal structure of the complex formed between rosmarinic acid and SARS-CoV-2 Mpro, revealing the molecular basis of its inhibitory mechanism. These findings not only enhance our understanding of the antiviral action of rosmarinic acid, but also provide valuable structural information and mechanistic insights for the further development of therapeutic strategies against SARS-CoV-2.

Improved SARS-CoV-2 main protease high-throughput screening assay using a 5-carboxyfluorescein substrate.

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) as a global threat to human health has highlighted the need for the development of novel therapies targeting current and emerging coronaviruses with pandemic potential. The coronavirus main protease (Mpro, also called 3CLpro) is a validated drug target against coronaviruses and has been heavily studied since the emergence of SARS-CoV-2 in late 2019. Here, we report the biophysical and enzymatic characterization of native Mpro, then characterize the steady-state kinetics of several commonly used FRET substrates, fluorogenic substrates, and six of the 11 reported SARS-CoV-2 polyprotein cleavage sequences. We then assessed the suitability of these substrates for high-throughput screening. Guided by our assessment of these substrates, we developed an improved 5-carboxyfluorescein-based FRET substrate, which is better suited for high-throughput screening and is less susceptible to interference and false positives than existing substrates. This study provides a useful framework for the design of coronavirus Mpro enzyme assays to facilitate the discovery and development of therapies targeting Mpro.

Engineering and Characterization of Avian Coronavirus Mutants Expressing Fluorescent Reporter Proteins from the Replicase Gene.

Infectious bronchitis virus (IBV) is an avian coronavirus that causes infectious bronchitis, an acute and highly contagious respiratory disease of chickens. IBV evolution under the pressure of comprehensive and widespread vaccination requires surveillance for vaccine resistance, as well as periodic vaccine updates. Reverse genetics systems are very valuable tools in virology, as they facilitate rapid genetic manipulation of viral genomes, thereby advancing basic and applied research. We report here the construction of an infectious clone of IBV strain Beaudette as a bacterial artificial chromosome (BAC). The engineered full-length IBV clone allowed the rescue of an infectious virus that was phenotypically indistinguishable from the parental virus. We used the infectious IBV clone and examined whether an enhanced green fluorescent protein (EGFP) can be produced by the replicase gene ORF1 and autocatalytically released from the replicase polyprotein through cleavage by the main coronavirus protease. We show that IBV tolerates insertion of the EGFP ORF at the 3' end of the replicase gene, between the sequences encoding nsp13 and nsp16 (helicase, RNA exonuclease, RNA endonuclease, and RNA methyltransferase). We further show that EGFP is efficiently cleaved from the replicase polyprotein and can be localized in double-membrane vesicles along with viral RNA polymerase and double-stranded RNA, an intermediate of IBV genome replication. One of the engineered reporter EGFP viruses were genetically stable during passage in cultured cells. We demonstrate that the reporter EGFP viruses can be used to study virus replication in host cells and for antiviral drug discovery and development of diagnostic assays. IMPORTANCE Reverse genetics systems based on bacterial artificial chromosomes (BACs) are the most valuable systems in coronavirus research. Here, we describe the establishment of a reverse genetics system for the avian coronavirus strain Beaudette, the most intensively studied strain. We cloned a copy of the avian coronavirus genome into a BAC vector and recovered infectious virus in permissive cells. We used the new system to construct reporter viruses that produce enhanced green fluorescent protein (EGFP). The EGFP coding sequence was inserted into 11 known cleavage sites of the major coronavirus protease in the replicase gene ORF1. Avian coronavirus tolerated the insertion of the EGFP coding sequence at three sites. The engineered reporter viruses replicated with parental efficiency in cultured cells and were sufficiently genetically stable. The new system facilitates functional genomics of the avian coronavirus genome but can also be used for the development of novel vaccines and anticoronaviral drugs.

Human coronavirus 3CL proteases cleave septins and disrupt Hedgehog signaling, causing ciliary dysfunction.

Coronaviruses target ciliate cells causing the loss of cilia, acute rhinorrheas, and other ciliopathies. The loss of ciliary function may help the virus infect, replicate, and spread. However, the molecular mechanisms by which coronaviruses cause ciliary defects are still unclear. Herein we demonstrate how coronavirus infection and severe acute respiratory syndrome coronavirus2 3CL protease induce cilia dysfunction by targeting a host protein septin that is required for the structure and function of cilia. Further, we demonstrate that coronaviruses and 3CL protease lead to the cleavage of several septin proteins (SEPT2, -6, and -9), producing cleaved obstructive fragments. Furthermore, ectopic expression of cleaved SEPT2 fragments shows defective ciliogenesis, disoriented septin filaments, and ablated Sonic Hedgehog (SHH) signaling in a protease activity-dependent manner. We present that the 3CLpro inhibitors are potent and prevent abnormal ciliary structures and SHH signaling. These results provide useful insights into the general mechanisms underlying ciliary defects caused by coronaviruses, which, in turn, facilitate virus spread and prove that preclinical and clinical 3CL protease inhibitors may prove useful as therapeutics for treating ciliary defects of coronaviruses.

Design, synthesis and biological evaluation of covalent peptidomimetic 3CL protease inhibitors containing nitrile moiety.

In this paper, a series of peptidomimetic SARS-CoV-2 3CL protease inhibitors with new P2 and P4 positions were synthesized and evaluated. Among these compounds, 1a and 2b exhibited obvious 3CLpro inhibitory activities with IC50 of 18.06 nM and 22.42 nM, respectively. 1a and 2b also showed excellent antiviral activities against SARS-CoV-2 in vitro with EC50 of 313.0 nM and 170.2 nM, respectively, the antiviral activities of 1a and 2b were 2- and 4-fold better than that of nirmatrelvir, respectively. In vitro studies revealed that these two compounds had no significant cytotoxicity. Further metabolic stability tests and pharmacokinetic studies showed that the metabolic stability of 1a and 2b in liver microsomes was significantly improved, and 2b had similar pharmacokinetic parameters to that of nirmatrelvir in mice.

Computational Screening of Inhibitory Compounds for SARS-Cov-2 3CL Protease with a Database Consisting of Approved and Investigational Chemicals.

Computational screening is one of the fundamental techniques in drug discovery. Each compound in a chemical database is bound to the target protein in virtual, and candidate compounds are selected from the binding scores. In this work, we carried out combinational computation of docking simulation to generate binding poses and molecular mechanics calculation to estimate binding scores. The coronavirus infectious disease has spread worldwide, and effective chemotherapy is strongly required. The viral 3-chymotrypsin-like (3CL) protease is a good target of low molecular-weight inhibitors. Hence, computational screening was performed to search for inhibitory compounds acting on the 3CL protease. As a preliminary assessment of the performance of this approach, we used 51 compounds for which inhibitory activity had already been confirmed. Docking simulations and molecular mechanics calculations were performed to evaluate binding scores. The preliminary evaluation suggested that our approach successfully selected the inhibitory compounds identified by the experiments. The same approach was applied to 8820 compounds in a database consisting of approved and investigational chemicals. Hence, docking simulations, molecular mechanics calculations, and re-evaluation of binding scores including solvation effects were performed, and the top 200 poses were selected as candidates for experimental assays. Consequently, 25 compounds were chosen for in vitro measurement of the enzymatic inhibitory activity. From the enzymatic assay, 5 compounds were identified to have inhibitory activities against the 3CL protease. The present work demonstrated the feasibility of a combination of docking simulation and molecular mechanics calculation for practical use in computational virtual screening.

Eupatin, a Flavonoid, Inhibits Coronavirus 3CL Protease and Replication.

The coronavirus disease 2019 (COVID-19) pandemic has caused more than six million deaths worldwide since 2019. Although vaccines are available, novel variants of coronavirus are expected to appear continuously, and there is a need for a more effective remedy for coronavirus disease. In this report, we isolated eupatin from Inula japonica flowers and showed that it inhibits the coronavirus 3 chymotrypsin-like (3CL) protease as well as viral replication. We showed that eupatin treatment inhibits SARS-CoV-2 3CL-protease, and computational modeling demonstrated that it interacts with key residues of 3CL-protease. Further, the treatment decreased the number of plaques formed by human coronavirus OC43 (HCoV-OC43) infection and decreased viral protein and RNA levels in the media. These results indicate that eupatin inhibits coronavirus replication.

The Ethanol Extract of Evodiae Fructus and Its Ingredient, Rutaecarpine, Inhibit Infection of SARS-CoV-2 and Inflammatory Responses.

COVID-19, derived from SARS-CoV-2, has resulted in millions of deaths and caused unprecedented socioeconomic damage since its outbreak in 2019. Although the vaccines developed against SARS-CoV-2 provide some protection, they have unexpected side effects in some people. Furthermore, new viral mutations reduce the effectiveness of the current vaccines. Thus, there is still an urgent need to develop potent non-vaccine therapeutics against this infectious disease. We recently established a series of detecting platforms to screen a large library of Chinese medicinal herbs and phytochemicals. Here, we reveal that the ethanolic extract of Evodiae Fructus and one of its components, rutaecarpine, showed promising potency in inhibiting the activity of 3C-like (3CL) protease, blocking the entry of the pseudo-typed SARS-CoV-2 (including wild-type and omicron) into cultured cells. In addition, inflammatory responses induced by pseudo-typed SARS-CoV-2 were markedly reduced by Evodiae Fructus extract and rutaecarpine. Together our data indicate that the herbal extract of Evodiae Fructus and rutaecarpine are potent anti-SARS-CoV-2 agents, which might be considered as a treatment against COVID-19 in clinical applications.

Discovery of Novel Chinese Medicine Compounds Targeting 3CL Protease by Virtual Screening and Molecular Dynamics Simulation.

The transmission and infectivity of COVID-19 have caused a pandemic that has lasted for several years. This is due to the constantly changing variants and subvariants that have evolved rapidly from SARS-CoV-2. To discover drugs with therapeutic potential for COVID-19, we focused on the 3CL protease (3CLpro) of SARS-CoV-2, which has been proven to be an important target for COVID-19 infection. Computational prediction techniques are quick and accurate enough to facilitate the discovery of drugs against the 3CLpro of SARS-CoV-2. In this paper, we used both ligand-based virtual screening and structure-based virtual screening to screen the traditional Chinese medicine small molecules that have the potential to target the 3CLpro of SARS-CoV-2. MD simulations were used to confirm these results for future in vitro testing. MCCS was then used to calculate the normalized free energy of each ligand and the residue energy contribution. As a result, we found ZINC15676170, ZINC09033700, and ZINC12530139 to be the most promising antiviral therapies against the 3CLpro of SARS-CoV-2.

The inhibitory effects of PGG and EGCG against the SARS-CoV-2 3C-like protease.

The coronavirus disease (COVID-19) pandemic, resulting from human-to-human transmission of a novel severe acute respiratory syndrome coronavirus (SARS-CoV-2), has led to a global health crisis. Given that the 3 chymotrypsin-like protease (3CLpro) of SARS-CoV-2 plays an indispensable role in viral polyprotein processing, its successful inhibition halts viral replication and thus constrains virus spread. Therefore, developing an effective SARS-CoV-2 3CLpro inhibitor to treat COVID-19 is imperative. A fluorescence resonance energy transfer (FRET)-based method was used to assess the proteolytic activity of SARS-CoV-2 3CLpro using intramolecularly quenched fluorogenic peptide substrates corresponding to the cleavage sequence of SARS-CoV-2 3CLpro. Molecular modeling with GEMDOCK was used to simulate the molecular interactions between drugs and the binding pocket of SARS-CoV-2 3CLpro. This study revealed that the Vmax of SARS-CoV-2 3CLpro was about 2-fold higher than that of SARS-CoV 3CLpro. Interestingly, the proteolytic activity of SARS-CoV-2 3CLpro is slightly more efficient than that of SARS-CoV 3CLpro. Meanwhile, natural compounds PGG and EGCG showed remarkable inhibitory activity against SARS-CoV-2 3CLpro than against SARS-CoV 3CLpro. In molecular docking, PGG and EGCG strongly interacted with the substrate binding pocket of SARS-CoV-2 3CLpro, forming hydrogen bonds with multiple residues, including the catalytic residues C145 and H41. The activities of PGG and EGCG against SARS-CoV-2 3CLpro demonstrate their inhibition of viral protease activity and highlight their therapeutic potentials for treating SARS-CoV-2 infection.

Inhibitors of Coronavirus 3CL Proteases Protect Cells from Protease-Mediated Cytotoxicity.

We describe a mammalian cell-based assay to identify coronavirus 3CL protease (3CLpro) inhibitors. This assay is based on rescuing protease-mediated cytotoxicity and does not require live virus. By enabling the facile testing of compounds across a range of 15 distantly related coronavirus 3CLpro enzymes, we identified compounds with broad 3CLpro-inhibitory activity. We also adapted the assay for use in compound screening and in doing so uncovered additional severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) 3CLpro inhibitors. We observed strong concordance between data emerging from this assay and those obtained from live-virus testing. The reported approach democratizes the testing of 3CLpro inhibitors by developing a simplified method for identifying coronavirus 3CLpro inhibitors that can be used by the majority of laboratories, rather than the few with extensive biosafety infrastructure. We identified two lead compounds, GC376 and compound 4, with broad activity against all 3CL proteases tested, including 3CLpro enzymes from understudied zoonotic coronaviruses. IMPORTANCE Multiple coronavirus pandemics have occurred over the last 2 decades. This has highlighted a need to be proactive in the development of therapeutics that can be readily deployed in the case of future coronavirus pandemics. We developed and validated a simplified cell-based assay for the identification of chemical inhibitors of 3CL proteases encoded by a wide range of coronaviruses. This assay is reporter free, does not require specialized biocontainment, and is optimized for performance in high-throughput screening. By testing reported 3CL protease inhibitors against a large collection of 3CL proteases with variable sequence similarity, we identified compounds with broad activity against 3CL proteases and uncovered structural insights into features that contribute to their broad activity. Furthermore, we demonstrated that this assay is suitable for identifying chemical inhibitors of proteases from families other than 3CL proteases.

Application of Acoustic Ejection MS System to High-Throughput Screening for SARS-CoV-2 3CL Protease Inhibitors.

MS is a powerful methodology for chemical screening to directly quantify substrates and products of enzymes, but its low throughput has been an issue. Recently, an acoustic liquid-handling apparatus (Echo®) used for rapid nano-dispensing has been coupled to a high-sensitivity mass spectrometer to create the Echo® MS system, and we applied this system to screening of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) 3CL protease inhibitors. Primary screening of 32033 chemical samples was completed in 12 h. Among the hits showing selective, dose-dependent 3CL-inhibitory activity, 8 compounds showed antiviral activity in cell-based assay.

Dimerization Tendency of 3CLpros of Human Coronaviruses Based on the X-ray Crystal Structure of the Catalytic Domain of SARS-CoV-2 3CLpro.

3CLpro of SARS-CoV-2 is a promising target for developing anti-COVID19 agents. In order to evaluate the catalytic activity of 3CLpros according to the presence or absence of the dimerization domain, two forms had been purified and tested. Enzyme kinetic studies with a FRET method revealed that the catalytic domain alone presents enzymatic activity, despite it being approximately 8.6 times less than that in the full domain. The catalytic domain was crystallized and its X-ray crystal structure has been determined to 2.3 Å resolution. There are four protomers in the asymmetric unit. Intriguingly, they were packed as a dimer though the dimerization domain was absent. The RMSD of superimposed two catalytic domains was 0.190 for 182 Cα atoms. A part of the long hinge loop (LH-loop) from Gln189 to Asp197 was not built in the model due to its flexibility. The crystal structure indicates that the decreased proteolytic activity of the catalytic domain was due to the incomplete construction of the substrate binding part built by the LH-loop. A structural survey with other 3CLpros showed that SARS-CoV families do not have interactions between DM-loop due to the conformational difference at the last turn of helix α7 compared with others. Therefore, we can conclude that the monomeric form contains nascent enzyme activity and that its efficiency increases by dimerization. This new insight may contribute to understanding the behavior of SARS-CoV-2 3CLpro and thus be useful in developing anti-COVID-19 agents.

A Study of Drug Repurposing to Identify SARS-CoV-2 Main Protease (3CLpro) Inhibitors.

The outbreak of coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) wreaked havoc all over the world. Although vaccines for the disease have recently become available and started to be administered to the population in various countries, there is still a strong and urgent need for treatments to cure COVID-19. One of the safest and fastest strategies is represented by drug repurposing (DRPx). In this study, thirty compounds with known safety profiles were identified from a chemical library of Phase II-and-up compounds through a combination of SOM Biotech's Artificial Intelligence (AI) technology, SOMAIPRO, and in silico docking calculations with third-party software. The selected compounds were then tested in vitro for inhibitory activity against SARS-CoV-2 main protease (3CLpro or Mpro). Of the thirty compounds, three (cynarine, eravacycline, and prexasertib) displayed strong inhibitory activity against SARS-CoV-2 3CLpro. VeroE6 cells infected with SARS-CoV-2 were used to find the cell protection capability of each candidate. Among the three compounds, only eravacycline showed potential antiviral activities with no significant cytotoxicity. A further study is planned for pre-clinical trials.

EGCG, a green tea polyphenol, inhibits human coronavirus replication in vitro.

COVID-19 pandemic results in record high deaths in many countries. Although a vaccine for SARS-CoV-2 is now available, effective antiviral drugs to treat coronavirus diseases are not available yet. Recently, EGCG, a green tea polyphenol, was reported to inhibit SARS-CoV-2 3CL-protease, however the effect of EGCG on coronavirus replication is unknown. In this report, human coronavirus HCoV-OC43 (beta coronavirus) and HCoV-229E (alpha coronavirus) were used to examine the effect of EGCG on coronavirus. EGCG treatment decreases 3CL-protease activity of HCoV-OC43 and HCoV-229E. Moreover, EGCG treatment decreased HCoV-OC43-induced cytotoxicity. Finally, we found that EGCG treatment decreased the levels of coronavirus RNA and protein in infected cell media. These results indicate that EGCG inhibits coronavirus replication.

Repurposing existing drugs: identification of SARS-CoV-2 3C-like protease inhibitors.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for coronavirus disease 2019 (COVID-19). Since its emergence, the COVID-19 pandemic has not only distressed medical services but also caused economic upheavals, marking urgent the need for effective therapeutics. The experience of combating SARS-CoV and MERS-CoV has shown that inhibiting the 3-chymotrypsin-like protease (3CLpro) blocks the replication of the virus. Given the well-studied properties of FDA-approved drugs, identification of SARS-CoV-2 3CLpro inhibitors in an FDA-approved drug library would be of great therapeutic value. Here, we screened a library consisting of 774 FDA-approved drugs for potent SARS-CoV-2 3CLpro inhibitors, using an intramolecularly quenched fluorescence (IQF) peptide substrate. Ethacrynic acid, naproxen, allopurinol, butenafine hydrochloride, raloxifene hydrochloride, tranylcypromine hydrochloride, and saquinavir mesylate have been found to block the proteolytic activity of SARS-CoV-2 3CLpro. The inhibitory activity of these repurposing drugs against SARS-CoV-2 3CLpro highlights their therapeutic potential for treating COVID-19 and other Betacoronavirus infections.

Rational design of flavonoid based potential inhibitors targeting SARS-CoV 3CL protease for the treatment of COVID-19.

The current outbreak of Coronavirus Disease 2019 (COVID-19) pandemic has reported thousands of deaths worldwide due to the rapid transmission rate and the lack of antiviral drugs and vaccinations. There is an urgent need to develop potential antiviral drug candidates for the prevention of COVID-19 infection. In the present study, a series of potential inhibitors targeting SARS-CoV 3CL protease were rationally designed by incorporating gamma lactam ring, and various fluoro substituted heterocyclic ring systems to the flavonoid scaffold. The prediction of drug-likeness, oral bioavailability, toxicity, synthetic accessibility, and ADMET properties was made by computational means. Quercetin was used as standard. The binding affinity of the ligands towards the 3CL protease target was examined using docking simulations. The designed ligands possess favourable pharmacokinetic and pharmacodynamic properties. Ligand L4, L8, and L14 appeared to be the lead compounds in the series and can be considered for further in-vivo and in-vitro validation.

Exploring the Binding Mechanism of PF-07321332 SARS-CoV-2 Protease Inhibitor through Molecular Dynamics and Binding Free Energy Simulations.

The novel coronavirus disease, caused by severe acute respiratory coronavirus 2 (SARS-CoV-2), rapidly spreading around the world, poses a major threat to the global public health. Herein, we demonstrated the binding mechanism of PF-07321332, α-ketoamide, lopinavir, and ritonavir to the coronavirus 3-chymotrypsin-like-protease (3CLpro) by means of docking and molecular dynamic (MD) simulations. The analysis of MD trajectories of 3CLpro with PF-07321332, α-ketoamide, lopinavir, and ritonavir revealed that 3CLpro-PF-07321332 and 3CLpro-α-ketoamide complexes remained stable compared with 3CLpro-ritonavir and 3CLpro-lopinavir. Investigating the dynamic behavior of ligand-protein interaction, ligands PF-07321332 and α-ketoamide showed stronger bonding via making interactions with catalytic dyad residues His41-Cys145 of 3CLpro. Lopinavir and ritonavir were unable to disrupt the catalytic dyad, as illustrated by increased bond length during the MD simulation. To decipher the ligand binding mode and affinity, ligand interactions with SARS-CoV-2 proteases and binding energy were calculated. The binding energy of the bespoke antiviral PF-07321332 clinical candidate was two times higher than that of α-ketoamide and three times than that of lopinavir and ritonavir. Our study elucidated in detail the binding mechanism of the potent PF-07321332 to 3CLpro along with the low potency of lopinavir and ritonavir due to weak binding affinity demonstrated by the binding energy data. This study will be helpful for the development and optimization of more specific compounds to combat coronavirus disease.

SARS-CoV-2 Main Protease Active Site Ligands in the Human Metabolome.

In late 2019, a global pandemic occurred. The causative agent was identified as a member of the Coronaviridae family, called severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). In this study, we present an analysis on the substances identified in the human metabolome capable of binding the active site of the SARS-CoV-2 main protease (Mpro). The substances present in the human metabolome have both endogenous and exogenous origins. The aim of this research was to find molecules whose biochemical and toxicological profile was known that could be the starting point for the development of antiviral therapies. Our analysis revealed numerous metabolites-including xenobiotics-that bind this protease, which are essential to the lifecycle of the virus. Among these substances, silybin, a flavolignan compound and the main active component of silymarin, is particularly noteworthy. Silymarin is a standardized extract of milk thistle, Silybum marianum, and has been shown to exhibit antioxidant, hepatoprotective, antineoplastic, and antiviral activities. Our results-obtained in silico and in vitro-prove that silybin and silymarin, respectively, are able to inhibit Mpro, representing a possible food-derived natural compound that is useful as a therapeutic strategy against COVID-19.