Seleno-Functionalization of Quercetin Improves the Non-Covalent Inhibition of Mpro and Its Antiviral Activity in Cells against SARS-CoV-2

The development of new antiviral drugs against SARS-CoV-2 is a valuable long-term strategy to protect the global population from the COVID-19 pandemic complementary to the vaccination. Considering this, the viral main protease (M^pro) is among the most promising molecular targets in light of its importance during the viral replication cycle. The natural flavonoid quercetin 1 has been recently reported to be a potent M^pro inhibitor in vitro, and we explored the effect produced by the introduction of organoselenium functionalities in this scaffold. In particular, we report here a new synthetic method to prepare previously inaccessible C-8 seleno-quercetin derivatives. By screening a small library of flavonols and flavone derivatives, we observed that some compounds inhibit the protease activity in vitro. For the first time, we demonstrate that quercetin (1) and 8-(p-tolylselenyl)quercetin (2d) block SARS-CoV-2 replication in infected cells at non-toxic concentrations, with an IC₅₀ of 192 μM and 8 μM, respectively. Based on docking experiments driven by experimental evidence, we propose a non-covalent mechanism for M^pro inhibition in which a hydrogen bond between the selenium atom and Gln189 residue in the catalytic pocket could explain the higher M^pro activity of 2d and, as a result, its better antiviral profile.     

A Structural Comparison of SARS-CoV-2 Main Protease and Animal Coronaviral Main Protease Reveals Species-Specific Ligand Binding and Dimerization Mechanism

Animal coronaviruses (CoVs) have been identified to be the origin of Severe Acute Respiratory Syndrome (SARS)-CoV, Middle East respiratory syndrome (MERS)-CoV, and probably SARS-CoV-2 that cause severe to fatal diseases in humans. Variations of zoonotic coronaviruses pose potential threats to global human beings. To overcome this problem, we focused on the main protease (M^pro), which is an evolutionary conserved viral protein among different coronaviruses. The broad-spectrum anti-coronaviral drug, GC376, was repurposed to target canine coronavirus (CCoV), which causes gastrointestinal infections in dogs. We found that GC376 can efficiently block the protease activity of CCoV M^pro and can thermodynamically stabilize its folding. The structure of CCoV M^pro in complex with GC376 was subsequently determined at 2.75 Å. GC376 reacts with the catalytic residue C144 of CCoV M^pro and forms an (R)- or (S)-configuration of hemithioacetal. A structural comparison of CCoV M^pro and other animal CoV M^pros with SARS-CoV-2 M^pro revealed three important structural determinants in a substrate-binding pocket that dictate entry and release of substrates. As compared with the conserved A141 of the S1 site and P188 of the S4 site in animal coronaviral M^pros, SARS-CoV-2 M^pro contains N142 and Q189 at equivalent positions which are considered to be more catalytically compatible. Furthermore, the conserved loop with residues 46–49 in animal coronaviral M^pros has been replaced by a stable α-helix in SARS-CoV-2 M^pro. In addition, the species-specific dimerization interface also influences the catalytic efficiency of CoV M^pros. Conclusively, the structural information of this study provides mechanistic insights into the ligand binding and dimerization of CoV M^pros among different species.     

Identification of Novel Small Molecules as Inhibitors of Hepatitis C Virus by Structure-Based Virtual Screening

Hepatitis C virus (HCV) NS3/NS4A serine protease is essential for viral replication, which is regarded as a promising drug target for developing direct-acting anti-HCV agents. In this study, sixteen novel compounds with cell-based HCV replicon activity ranging from 3.0 to 28.2 μM (IC₅₀) were successfully identified by means of structure-based virtual screening. Compound 5 and compound 11, with an IC₅₀ of 3.0 μM and 5.1 μM, respectively, are the two most potent molecules with low cytotoxicity.     

Iterated Virtual Screening-Assisted Antiviral and Enzyme Inhibition Assays Reveal the Discovery of Novel Promising Anti-SARS-CoV-2 with Dual Activity

Unfortunately, COVID-19 is still a threat to humankind and has a dramatic impact on human health, social life, the world economy, and food security. With the limited number of suggested therapies under clinical trials, the discovery of novel therapeutic agents is essential. Here, a previously identified anti-SARS-CoV-2 compound named Compound 13 (1,2,5-Oxadiazole-3-carboximidic acid, 4,4′-(methylenediimino) bis,bis[[(2-hydroxyphenyl)methylene]hydrazide) was subjected to an iterated virtual screening against SARS-CoV-2 M^pro using a combination of Ligand Designer and PathFinder. PathFinder, a computational reaction enumeration tool, was used for the rapid generation of enumerated structures via default reaction library. Ligand designer was employed for the computerized lead optimization and selection of the best structural modification that resulted in a favorable ligand–protein complex. The obtained compounds that showed the best binding to M^pro were re-screened against TMPRSS2, leading to the identification of 20 shared compounds. The compounds were further visually inspected, which resulted in the identification of five shared compounds M1–5 with dual binding affinity. In vitro evaluation and enzyme inhibition assay indicated that M3, an analogue of Compound 13 afforded by replacing the phenolic moiety with pyridinyl, possesses an improved antiviral activity and safety. M3 displayed in vitro antiviral activity with IC₅₀ 0.016 µM and M^pro inhibition activity with IC₅₀ 0.013 µM, 7-fold more potent than the parent Compound 13 and potent than the antivirals drugs that are currently under clinical trials. Moreover, M3 showed potent activity against human TMPRSS2 and furin enzymes with IC₅₀ 0.05, and 0.08 µM, respectively. Molecular docking, WaterMap analysis, molecular dynamics simulation, and R-group analysis confirmed the superiority of the binding fit to M3 with the target enzymes. WaterMap analysis calculated the thermodynamic properties of the hydration site in the binding pocket that significantly affects the biological activity. Loading M3 on zinc oxide nanoparticles (ZnO NPs) increased the antiviral activity of the compound 1.5-fold, while maintaining a higher safety profile. In conclusion, lead optimized discovery following an iterated virtual screening in association with molecular docking and biological evaluation revealed a novel compound named M3 with promising dual activity against SARS-CoV-2. The compound deserves further investigation for potential clinical-based studies.     

The Se–S Bond Formation in the Covalent Inhibition Mechanism of SARS-CoV-2 Main Protease by Ebselen-like Inhibitors: A Computational Study

The inhibition mechanism of the main protease (M^pro) of SARS-CoV-2 by ebselen (EBS) and its analog with a hydroxyl group at position 2 of the benzisoselenazol-3(2H)-one ring (EBS-OH) was studied by using a density functional level of theory. Preliminary molecular dynamics simulations on the apo form of M^pro were performed taking into account both the hydrogen donor and acceptor natures of the Nδ and Nε of His41, a member of the catalytic dyad. The potential energy surfaces for the formation of the Se–S covalent bond mediated by EBS and EBS-OH on M^pro are discussed in detail. The EBS-OH shows a distinctive behavior with respect to EBS in the formation of the noncovalent complex. Due to the presence of canonical H-bonds and noncanonical ones involving less electronegative atoms, such as sulfur and selenium, the influence on the energy barriers and reaction energy of the Minnesota hybrid meta-GGA functionals M06, M06-2X and M08HX, and the more recent range-separated hybrid functional wB97X were also considered. The knowledge of the inhibition mechanism of M^pro by the small protease inhibitors EBS or EBS-OH can enlarge the possibilities for designing more potent and selective inhibitor-based drugs to be used in combination with other antiviral therapies.     

In silico allicin induced S-thioallylation of SARS-CoV-2 main protease

Coronavirus disease 2019 (COVID-19) is an ongoing pandemic caused due to new coronavirus infection with ³716075 deaths across the world as reported by the World Health Organization (WHO). SARS-CoV-2 main protease (M^pro) plays a vital role in the replication of coronavirus and thus an attractive target for the screening of inhibitors for the therapy of COVID-19. The preclinical drugs ebselen and PX-12 are potent inhibitors of SARS-CoV-2 M^pro and covalently modifies the active site Cys-145 residue of M^pro through selenosulfide/disulfide. In the current report, using virtual screening methods, reactive sulfur species allicin is subjecting for covalent docking at the active site of SARS-CoV-2 M^pro using PX-12 as a benchmark reference compound. The results indicate that allicin induces dual S-thioallylation of Cys-145 and Cys-85/ Cys-156 residues of SARS-CoV-2 M^pro. Using density functional theory (DFT), Gibbs free energy change (DG) is calculated for the putative reactions between N-acetylcysteine amide thiol and allicin/allyl sulfenic acid. The overall reaction is exergonic and allyl disulfide of Cys-145 residue of M^pro is involved in a sulfur mediated hydrogen bond. The results indicate that allicin causes dual S-thioallylation of SARS-CoV-2 M^pro which may be of interest for treatment and attenuation of ongoing coronavirus infection.     

De novo ssRNA Aptamers against the SARS-CoV-2 Main Protease: In Silico Design and Molecular Dynamics Simulation

Herein, we have generated ssRNA aptamers to inhibit SARS-CoV-2 M^pro, a protease necessary for the SARS-CoV-2 coronavirus replication. Because there is no aptamer 3D structure currently available in the databanks for this protein, first, we modeled an ssRNA aptamer using an entropic fragment-based strategy. We refined the initial sequence and 3D structure by using two sequential approaches, consisting of an elitist genetic algorithm and an RNA inverse process. We identified three specific aptamers against SARS-CoV-2 M^pro, called MApta^pro, MApta^pro-IR1, and MApta^pro-IR2, with similar 3D conformations and that fall in the dimerization region of the SARS-CoV-2 M^pro necessary for the enzymatic activity. Through the molecular dynamic simulation and binding free energy calculation, the interaction between the MApta^pro-IR1 aptamer and the SARS-CoV-2 M^pro enzyme resulted in the strongest and the highest stable complex; therefore, the ssRNA MApta^pro-IR1 aptamer was selected as the best potential candidate for the inhibition of SARS-CoV-2 M^pro and a perspective therapeutic drug for the COVID-19 disease.     

In Vitro Inhibition of Renal OCT2 and MATE1 Secretion by Antiemetic Drugs

The organic cation transporter 2 (OCT2) and multidrug and toxin extrusion protein 1 (MATE1) mediate the renal secretion of drugs. Recent studies suggest that ondansetron, a 5-HT₃ antagonist drug used to prevent nausea and vomiting, can inhibit OCT2- and MATE1-mediated transport. The purpose of this study was to test the ability of five 5-HT₃ antagonist drugs to inhibit the OCT2 and MATE1 transporters. The transport of the OCT2/MATE1 probe substrate ASP⁺ was assessed using two models: (1) HEK293 kidney cells overexpressing human OCT2 or MATE1, and (2) MDCK cells transfected with human OCT2 and MATE1. In HEK293 cells, the inhibition of ASP⁺ uptake by OCT2 listed in order of potency was palonosetron (IC₅₀: 2.6 μM) > ondansetron > granisetron > tropisetron > dolasetron (IC₅₀: 85.4 μM) and the inhibition of ASP⁺ uptake by MATE1 in order of potency was ondansetron (IC₅₀: 0.1 μM) > palonosetron = tropisetron > granisetron > dolasetron (IC₅₀: 27.4 μM). Ondansetron (0.5–20 μM) inhibited the basolateral-to-apical transcellular transport of ASP⁺ up to 64%. Higher concentrations (10 and 20 μM) of palonosetron, tropisetron, and dolasetron similarly reduced the transcellular transport of ASP⁺. In double-transfected OCT2-MATE1 MDCK cells, ondansetron at concentrations of 0.5 and 2.5 μM caused significant intracellular accumulation of ASP⁺. Taken together, these data suggest that 5-HT₃ antagonist drugs may inhibit the renal secretion of cationic drugs by interfering with OCT2 and/or MATE1 function.     

Design,Synthesis and Cytotoxic Evaluation of o-Carboxamido Stilbene Analogues

Resveratrol, a natural stilbene found in grapes and wines exhibits a wide range of pharmacological properties. Resveratrol is also known as a good chemopreventive agent for inhibiting carcinogenesis processes that target kinases, cyclooxygenases, ribonucleotide reductase and DNA polymerases. A total of 19 analogues with an amide moiety were synthesized and the cytotoxic effects of the analogues on a series of human cancer cell lines are reported. Three compounds 6d, 6i and 6n showed potent cytotoxicity against prostate cancer DU-145 (IC₅₀ = 16.68 μM), colon cancer HT-29 (IC₅₀ = 7.51 μM) and breast cancer MCF-7 (IC₅₀ = 21.24 μM), respectively, which are comparable with vinblastine. The resveratrol analogues were synthesized using the Heck method.     

Human Superantibodies to 3CLpro Inhibit Replication of SARS-CoV-2 across Variants

Broadly effective and safe anti-coronavirus agent is existentially needed. Major protease (3CL^pro) is a highly conserved enzyme of betacoronaviruses. The enzyme plays pivotal role in the virus replication cycle. Thus, it is a good target of a broadly effective anti-Betacoronavirus agent. In this study, human single-chain antibodies (HuscFvs) of the SARS-CoV-2 3CL^pro were generated using phage display technology. The 3CL^pro-bound phages were used to infect Escherichia coli host for the production the 3CL^pro-bound HuscFvs. Computerized simulation was used to guide the selection of the phage infected-E. coli clones that produced HuscFvs with the 3CL^pro inhibitory potential. HuscFvs of three phage infected-E. coli clones were predicted to form contact interface with residues for 3CL^pro catalytic activity, substrate binding, and homodimerization. These HuscFvs were linked to a cell-penetrating peptide to make them cell-penetrable, i.e., became superantibodies. The superantibodies blocked the 3CL^pro activity in vitro, were not toxic to human cells, traversed across membrane of 3CL^pro-expressing cells to co-localize with the intracellular 3CL^pro and most of all, they inhibited replication of authentic SARS-CoV-2 Wuhan wild type and α, β, δ, and Omicron variants that were tested. The superantibodies should be investigated further towards clinical application as a safe and broadly effective anti-Betacoronavirus agent.     

Synthesis and Biological Evaluation of Novel Dehydroabietic Acid Derivatives Conjugated with Acyl-Thiourea Peptide Moiety as Antitumor Agents

A series of dehydroabietic acid (DHAA) acyl-thiourea derivatives were designed and synthesized as potent antitumor agents. The in vitro pharmacological screening results revealed that the target compounds exhibited potent cytotoxicity against HeLa, SK-OV-3 and MGC-803 tumor cell lines, while they showed lower cytotoxicity against HL-7702 normal human river cells. Compound 9n (IC₅₀ = 6.58 ± 1.11 μM) exhibited the best antitumor activity against the HeLa cell line and even displayed more potent inhibitory activity than commercial antitumor drug 5-FU (IC₅₀ = 36.58 ± 1.55 μM). The mechanism of representative compound 9n was then studied by acridine orange/ethidium bromide staining, Hoechst 33,258 staining, JC-1 mitochondrial membrane potential staining, TUNEL assay and flow cytometry, which illustrated that this compound could induce apoptosis in HeLa cells. Cell cycle analysis indicated that compound 9n mainly arrested HeLa cells in the S phase stage. Further investigation demonstrated that compound 9n induced apoptosis of HeLa cells through a mitochondrial pathway.     

Access to New Cytotoxic Triterpene and Steroidal Acid-TEMPO Conjugates by Ugi Multicomponent-Reactions

Multicomponent reactions, especially the Ugi-four component reaction (U-4CR), provide powerful protocols to efficiently access compounds having potent biological and pharmacological effects. Thus, a diverse library of betulinic acid (BA), fusidic acid (FA), cholic acid (CA) conjugates with TEMPO (nitroxide) have been prepared using this approach, which also makes them applicable in electron paramagnetic resonance (EPR) spectroscopy. Moreover, convertible amide modified spin-labelled fusidic acid derivatives were selected for post-Ugi modification utilizing a wide range of reaction conditions which kept the paramagnetic center intact. The nitroxide labelled betulinic acid analogue 6 possesses cytotoxic effects towards two investigated cell lines: prostate cancer PC3 (IC₅₀ 7.4 ± 0.7 μM) and colon cancer HT29 (IC₅₀ 9.0 ± 0.4 μM). Notably, spin-labelled fusidic acid derivative 8 acts strongly against these two cancer cell lines (PC3: IC₅₀ 6.0 ± 1.1 μM; HT29: IC₅₀ 7.4 ± 0.6 μM). Additionally, another fusidic acid analogue 9 was also found to be active towards HT29 with IC₅₀ 7.0 ± 0.3 μM (CV). Studies on the mode of action revealed that compound 8 increased the level of caspase-3 significantly which clearly indicates induction of apoptosis by activation of the caspase pathway. Furthermore, the exclusive mitochondria targeting of compound 18 was successfully achieved, since mitochondria are the major source of ROS generation.     

The Pharmacological Profile of Second Generation Pyrovalerone Cathinones and Related Cathinone Derivative

Pyrovalerone cathinones are potent psychoactive substances that possess a pyrrolidine moiety. Pyrovalerone-type novel psychoactive substances (NPS) are continuously detected but their pharmacology and toxicology are largely unknown. We assessed several pyrovalerone and related cathinone derivatives at the human norepinephrine (NET), dopamine (DAT), and serotonin (SERT) uptake transporters using HEK293 cells overexpressing each respective transporter. We examined the transporter-mediated monoamine efflux in preloaded cells. The receptor binding and activation potency was also assessed at the 5-HT_1A, 5-HT_2A, 5-HT_2B, and 5-HT_2C receptors. All pyrovalerone cathinones were potent DAT (IC₅₀ = 0.02–8.7 μM) and NET inhibitors (IC₅₀ = 0.03–4.6 μM), and exhibited no SERT activity at concentrations < 10 μM. None of the compounds induced monoamine efflux. NEH was a potent DAT/NET inhibitor (IC₅₀ = 0.17–0.18 μM). 4F-PBP and NEH exhibited a high selectivity for the DAT (DAT/SERT ratio = 264–356). Extension of the alkyl chain enhanced NET and DAT inhibition potency, while presence of a 3,4-methylenedioxy moiety increased SERT inhibition potency. Most compounds did not exhibit any relevant activity at other monoamine receptors. In conclusion, 4F-PBP and NEH were selective DAT/NET inhibitors indicating that these substances likely produce strong psychostimulant effects and have a high abuse liability.     

Synthesis and Anti-Breast Cancer Evaluation of Novel N-(Guanidinyl)benzenesulfonamides

A series of 4-(substituted)-N-(guanidinyl)benzenesulfonamides bearing biologically active pyrazole, pyrimidine and pyridine moieties were prepared and evaluated for their anticancer activity against human tumor breast cell line (MCF7). These sulfonamides showed promising activity with IC₅₀ values ranging from 49.5 to 70.2 μM. The structure-activity relationship of the synthesized compounds was studied. Interestingly, it was found that the most potent compounds in this study were the corresponding 2-cyanoacrylate 3, 3-oxobutanoate 4, pyrazole 6, pyridine 9 and pyrazole 13. Compounds 7 and 8 are nearly as active as Doxorubicin as reference drug with (IC₅₀ values = 70.2, 68.1 μM), while compounds 5, 10 and 11 exhibited a moderate activity.     

Computational Selectivity Assessment of Protease Inhibitors against SARS-CoV-2

The pandemic of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) poses a serious global health threat. Since no specific therapeutics are available, researchers around the world screened compounds to inhibit various molecular targets of SARS-CoV-2 including its main protease (M^pro) essential for viral replication. Due to the high urgency of these discovery efforts, off-target binding, which is one of the major reasons for drug-induced toxicity and safety-related drug attrition, was neglected. Here, we used molecular docking, toxicity profiling, and multiple molecular dynamics (MD) protocols to assess the selectivity of 33 reported non-covalent inhibitors of SARS-CoV-2 M^pro against eight proteases and 16 anti-targets. The panel of proteases included SARS-CoV M^pro, cathepsin G, caspase-3, ubiquitin carboxy-terminal hydrolase L1 (UCHL1), thrombin, factor Xa, chymase, and prostasin. Several of the assessed compounds presented considerable off-target binding towards the panel of proteases, as well as the selected anti-targets. Our results further suggest a high risk of off-target binding to chymase and cathepsin G. Thus, in future discovery projects, experimental selectivity assessment should be directed toward these proteases. A systematic selectivity assessment of SARS-CoV-2 M^pro inhibitors, as we report it, was not previously conducted.     

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 (3CL^pro) by means of docking and molecular dynamic (MD) simulations. The analysis of MD trajectories of 3CL^pro with PF-07321332, α-ketoamide, lopinavir, and ritonavir revealed that 3CL^pro–PF-07321332 and 3CL^pro–α-ketoamide complexes remained stable compared with 3CL^pro–ritonavir and 3CL^pro–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 3CL^pro. 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 3CL^pro 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.     

Mass Spectrometric Assays Reveal Discrepancies in Inhibition Profiles for the SARS-CoV-2 Papain-Like Protease

The two SARS-CoV-2 proteases, i. e. the main protease (M^pro) and the papain-like protease (PL^pro), which hydrolyze the viral polypeptide chain giving functional non-structural proteins, are essential for viral replication and are medicinal chemistry targets. We report a high-throughput mass spectrometry (MS)-based assay which directly monitors PL^pro catalysis in vitro. The assay was applied to investigate the effect of reported small-molecule PL^pro inhibitors and selected M^pro inhibitors on PL^pro catalysis. The results reveal that some, but not all, PL^pro inhibitor potencies differ substantially from those obtained using fluorescence-based assays. Some substrate-competing M^pro inhibitors, notably PF-07321332 (nirmatrelvir) which is in clinical development, do not inhibit PL^pro. Less selective M^pro inhibitors, e. g. auranofin, inhibit PL^pro, highlighting the potential for dual PL^pro/M^pro inhibition. MS-based PL^pro assays, which are orthogonal to widely employed fluorescence-based assays, are of utility in validating inhibitor potencies, especially for inhibitors operating by non-covalent mechanisms.     

Luminescent Assay for the Screening of SARS-CoV-2 MPro Inhibitors

Since the outbreak of SARS-CoV-2 in December 2019 millions of infections have been reported globally. The viral chymotrypsin-like main protease (M^Pro) exhibits a crucial role in viral replication and represents a relevant target for antiviral drug development. In order to screen potential M^Pro inhibitors we developed a luminescent assay using a peptide based probe containing a cleavage site specific for M^Pro. This assay was validated showing IC₅₀ values similar to those reported in the literature for known M^Pro inhibitors and can be used to screen new inhibitors.