Interaction of the prototypical α-ketoamide inhibitor with the SARS-CoV-2 main protease active site in silico: Molecular dynamic simulations highlight the stability of the ligand-protein complex

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causes an illness known as COVID-19, which has been declared a global pandemic with over 2 million confirmed cases and 137,000 deaths in 185 countries and regions at the time of writing (16 April 2020), over a quarter of these cases being in the United States. In the absence of a vaccine, or an approved effective therapeutic, there is an intense interest in repositioning available drugs or designing small molecule antivirals. In this context, in silico modelling has proven to be an invaluable tool. An important target is the SARS-CoV-2 main protease (M^pro), involved in processing translated viral proteins. Peptidomimetic α-ketoamides represent prototypical inhibitors of M^pro. A recent attempt at designing a compound with enhanced pharmacokinetic properties has resulted in the synthesis and evaluation of the α-ketoamide 13b analogue. Here, we performed molecular docking and molecular dynamics simulations to further characterize the interaction of α-ketoamide 13b with the active site of the SARS-CoV-2 M^pro. We included the widely used antibiotic, amoxicillin, for comparison. Our findings indicate that α-ketoamide 13b binds more tightly (predicted GlideScore = -8.7 and -9.2 kcal/mol for protomers A and B, respectively), to the protease active site compared to amoxicillin (-5.0 and -4.8 kcal/mol). Further, molecular dynamics simulations highlight the stability of the interaction of the α-ketoamide 13b ligand with the SARS-CoV-2 M^pro (ΔG = -25.2 and -22.3 kcal/mol for protomers A and B). In contrast, amoxicillin interacts unfavourably with the protease (ΔG = +32.8 kcal/mol for protomer A), with unbinding events observed in several independent simulations. Overall, our findings are consistent with those previously observed, and highlight the need to further explore the α-ketoamides as potential antivirals for this ongoing COVID-19 pandemic.     

Interaction of small molecules with the SARS-CoV-2 main protease in silico and in vitro validation of potential lead compounds using an enzyme-linked immunosorbent assay

Caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the COVID-19 pandemic is ongoing, with no proven safe and effective vaccine to date. Further, effective therapeutic agents for COVID-19 are limited, and as a result, the identification of potential small molecule antiviral drugs is of particular importance. A critical antiviral target is the SARS-CoV-2 main protease (M^pro), and our aim was to identify lead compounds with potential inhibitory effects. We performed an initial molecular docking screen of 300 small molecules, which included phenolic compounds and fatty acids from our OliveNet™ library (224), and an additional group of curated pharmacological and dietary compounds. The prototypical α-ketoamide 13b inhibitor was used as a control to guide selection of the top 30 compounds with respect to binding affinity to the M^pro active site. Further studies and analyses including blind docking were performed to identify hypericin, cyanidin-3-O-glucoside and SRT2104 as potential leads. Molecular dynamics simulations demonstrated that hypericin (ΔG = -18.6 and -19.3 kcal/mol), cyanidin-3-O-glucoside (ΔG = -50.8 and -42.1 kcal/mol), and SRT2104 (ΔG = -8.7 and -20.6 kcal/mol), formed stable interactions with the M^pro active site. An enzyme-linked immunosorbent assay indicated that, albeit, not as potent as the covalent positive control (GC376), our leads inhibited the M^pro with activity in the micromolar range, and an order of effectiveness of hypericin and cyanidin-3-O-glucoside > SRT2104 > SRT1720. Overall, our findings, and those highlighted by others indicate that hypericin and cyanidin-3-O-glucoside are suitable candidates for progress to in vitro and in vivo antiviral studies.     

Fullerenes against COVID-19: Repurposing C60 and C70 to Clog the Active Site of SARS-CoV-2 Protease

The persistency of COVID-19 in the world and the continuous rise of its variants demand new treatments to complement vaccines. Computational chemistry can assist in the identification of moieties able to lead to new drugs to fight the disease. Fullerenes and carbon nanomaterials can interact with proteins and are considered promising antiviral agents. Here, we propose the possibility to repurpose fullerenes to clog the active site of the SARS-CoV-2 protease, M^pro. Through the use of docking, molecular dynamics, and energy decomposition techniques, it is shown that C₆₀ has a substantial binding energy to the main protease of the SARS-CoV-2 virus, M^pro, higher than masitinib, a known inhibitor of the protein. Furthermore, we suggest the use of C₇₀ as an innovative scaffold for the inhibition of SARS-CoV-2 M^pro. At odds with masitinib, both C₆₀ and C₇₀ interact more strongly with SARS-CoV-2 M^pro when different protonation states of the catalytic dyad are considered. The binding of fullerenes to M^pro is due to shape complementarity, i.e., vdW interactions, and is aspecific. As such, it is not sensitive to mutations that can eliminate or invert the charges of the amino acids composing the binding pocket. Fullerenic cages should therefore be more effective against the SARS-CoV-2 virus than the available inhibitors such as masinitib, where the electrostatic term plays a crucial role in the binding.     

Computational Screening of Phenylamino-Phenoxy-Quinoline Derivatives against the Main Protease of SARS-CoV-2 Using Molecular Docking and the ONIOM Method

In the search for new anti-HIV-1 agents, two forms of phenylamino-phenoxy-quinoline derivatives have been synthesized, namely, 2-phenylamino-4-phenoxy-quinoline and 6-phenylamino-4-phenoxy-quinoline. In this study, the binding interactions of phenylamino-phenoxy-quinoline derivatives and six commercially available drugs (hydroxychloroquine, ritonavir, remdesivir, S-217622, N3, and PF-07321332) with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease (M^pro) were investigated using molecular docking and the ONIOM method. The molecular docking showed the hydrogen bonding and hydrophobic interactions of all the compounds in the pocket of SARS-CoV-2 main protease (M^pro), which plays an important role for the division and proliferation of the virus into the cell. The binding free energy values between the ligands and M^pro ranged from −7.06 to −10.61 kcal/mol. The molecular docking and ONIOM results suggested that 4-(2′,6′-dimethyl-4′-cyanophenoxy)-2-(4″-cyanophenyl)-aminoquinoline and 4-(4′-cyanophenoxy)-2-(4″-cyanophenyl)-aminoquinoline have low binding energy values and appropriate molecular properties; moreover, both compounds could bind to M^pro via hydrogen bonding and Pi-Pi stacking interactions with amino acid residues, namely, HIS41, GLU166, and GLN192. These amino acids are related to the proteolytic cleavage process of the catalytic triad mechanisms. Therefore, this study provides important information for further studies on synthetic quinoline derivatives as antiviral candidates in the treatment of SARS-CoV-2.     

Design of Novel Enantiopure Dispirooxindolopyrrolidine-Piperidones as Promising Candidates toward COVID-19: Asymmetric Synthesis,Crystal Structure and In Silico Studies

Despite the effectiveness of COVID-19 vaccines, there is still an urgent need for discovering new anti-viral drugs to address the awful spread and transmission of the rapidly modifiable virus. In this study, the ability of a small library of enantiomerically pure spirooxindolopyrrolidine-grafted piperidones to inhibit the main protease of SARS-CoV-2 (M^pro) is evaluated. These spiroheterocycles were synthesized by 1,3-dipolar cycloaddition of various stabilized azomethine ylides with chiral dipolarophiles derived from N-[(S)-(-)-methylbenzyl]-4-piperidone. The absolute configuration of contiguous carbons was confirmed by a single crystal X-ray diffraction analysis. The binding of these compounds to SARS-CoV-2 M^pro was investigated using molecular docking and molecular dynamics simulation. Three compounds 4a, 4b and 4e exhibited stable binding modes interacting with the key subsites of the substrate-binding pocket of SARS-CoV-2 M^pro. The synthesized compounds represent potential leads for the development of novel inhibitors of SARS-CoV-2 main protease protein for COVID-19 treatment.     

Inhibitor binding influences the protonation states of histidines in SARS-CoV-2 main protease

The main protease (M^pro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics. Recently, many high-resolution apo and inhibitor-bound structures of M^pro, a cysteine protease, have been determined, facilitating structure-based drug design. M^pro plays a central role in the viral life cycle by catalyzing the cleavage of SARS-CoV-2 polyproteins. In addition to the catalytic dyad His41–Cys145, M^pro contains multiple histidines including His163, His164, and His172. The protonation states of these histidines and the catalytic nucleophile Cys145 have been debated in previous studies of SARS-CoV M^pro, but have yet to be investigated for SARS-CoV-2. In this work we have used molecular dynamics simulations to determine the structural stability of SARS-CoV-2 M^pro as a function of the protonation assignments for these residues. We simulated both the apo and inhibitor-bound enzyme and found that the conformational stability of the binding site, bound inhibitors, and the hydrogen bond networks of M^pro are highly sensitive to these assignments. Additionally, the two inhibitors studied, the peptidomimetic N3 and an α-ketoamide, display distinct His41/His164 protonation-state-dependent stabilities. While the apo and the N3-bound systems favored N_δ (HD) and N_ϵ (HE) protonation of His41 and His164, respectively, the α-ketoamide was not stably bound in this state. Our results illustrate the importance of using appropriate histidine protonation states to accurately model the structure and dynamics of SARS-CoV-2 M^pro in both the apo and inhibitor-bound states, a necessary prerequisite for drug-design efforts.

The main protease (M^pro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics.     

Novel Small-Molecule Scaffolds as Candidates against the SARS Coronavirus 2 Main Protease: A Fragment-Guided in Silico Approach

The ongoing pandemic caused by the novel coronavirus has been the greatest global health crisis since the Spanish flu pandemic of 1918. Thus far, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in over 1 million deaths, and there is no cure or vaccine to date. The recently solved crystal structure of the SARS-CoV-2 main protease has been a major focus for drug-discovery efforts. Here, we present a fragment-guided approach using ZINCPharmer, where 17 active fragments known to bind to the catalytic centre of the SARS-CoV-2 main protease (SARS-CoV-2 M^pro) were used as pharmacophore queries to search the ZINC databases of natural compounds and natural derivatives. This search yielded 134 hits that were then subjected to multiple rounds of in silico analyses, including blind and focused docking against the 3D structure of the main protease. We scrutinised the poses, scores, and protein–ligand interactions of 15 hits and selected 7. The scaffolds of the seven hits were structurally distinct from known inhibitor scaffolds, thus indicating scaffold novelty. Our work presents several novel scaffolds as potential candidates for experimental validation against SARS-CoV-2 M^pro.     

Discovery of SARS-CoV-2 Mpro peptide inhibitors from modelling substrate and ligand binding

The main protease (M^pro) of SARS-CoV-2 is central to viral maturation and is a promising drug target, but little is known about structural aspects of how it binds to its 11 natural cleavage sites. We used biophysical and crystallographic data and an array of biomolecular simulation techniques, including automated docking, molecular dynamics (MD) and interactive MD in virtual reality, QM/MM, and linear-scaling DFT, to investigate the molecular features underlying recognition of the natural M^pro substrates. We extensively analysed the subsite interactions of modelled 11-residue cleavage site peptides, crystallographic ligands, and docked COVID Moonshot-designed covalent inhibitors. Our modelling studies reveal remarkable consistency in the hydrogen bonding patterns of the natural M^pro substrates, particularly on the N-terminal side of the scissile bond. They highlight the critical role of interactions beyond the immediate active site in recognition and catalysis, in particular plasticity at the S2 site. Building on our initial M^pro-substrate models, we used predictive saturation variation scanning (PreSaVS) to design peptides with improved affinity. Non-denaturing mass spectrometry and other biophysical analyses confirm these new and effective ‘peptibitors’ inhibit M^pro competitively. Our combined results provide new insights and highlight opportunities for the development of M^pro inhibitors as anti-COVID-19 drugs.

The main protease (M^pro) of SARS-CoV-2 is central to viral maturation and is a promising drug target. In silico methods reveal structural aspects of how it binds to its 11 natural cleavage sites, the design of novel peptide inhibitors, and insights into drug design.     

The S1′–S3′ Pocket of the SARS-CoV-2 Main Protease Is Critical for Substrate Selectivity and Can Be Targeted with Covalent Inhibitors

The main protease (M^pro) of SARS-CoV-2 is a well-characterized target for antiviral drug discovery. To date, most antiviral drug discovery efforts have focused on the S4–S1′ pocket of M^pro; however, it is still unclear whether the S1′–S3′ pocket per se can serve as a new site for drug discovery. In this study, the S1′–S3′ pocket of M^pro was found to differentially recognize viral peptidyl substrates. For instance, S3′ in M^pro strongly favors Phe or Trp, and S1′ favors Ala. The peptidyl inhibitor D-4–77, which possesses an α-bromoacetamide warhead, was discovered to be a promising inhibitor of M^pro, with an IC₅₀ of 0.95 μM and an antiviral EC₅₀ of 0.49 μM. The M^pro/inhibitor co-crystal structure confirmed the binding mode of the inhibitor to the S1′–S3′ pocket and revealed a covalent mechanism. In addition, D-4–77 functions as an immune protectant and suppresses SARS-CoV-2 M^pro-induced antagonism of the host NF-κB innate immune response. These findings indicate that the S1′–S3′ pocket of SARS-CoV-2 M^pro is druggable, and that inhibiting SARS-CoV-2 M^pro can simultaneously protect human innate immunity and inhibit virion assembly.     

In Silico Identification and Validation of Organic Triazole Based Ligands as Potential Inhibitory Drug Compounds of SARS-CoV-2 Main Protease

The SARS-CoV-2 virus is highly contagious to humans and has caused a pandemic of global proportions. Despite worldwide research efforts, efficient targeted therapies against the virus are still lacking. With the ready availability of the macromolecular structures of coronavirus and its known variants, the search for anti-SARS-CoV-2 therapeutics through in silico analysis has become a highly promising field of research. In this study, we investigate the inhibiting potentialities of triazole-based compounds against the SARS-CoV-2 main protease (M^pro). The SARS-CoV-2 main protease (M^pro) is known to play a prominent role in the processing of polyproteins that are translated from the viral RNA. Compounds were pre-screened from 171 candidates (collected from the DrugBank database). The results showed that four candidates (Bemcentinib, Bisoctrizole, PYIITM, and NIPFC) had high binding affinity values and had the potential to interrupt the main protease (M^pro) activities of the SARS-CoV-2 virus. The pharmacokinetic parameters of these candidates were assessed and through molecular dynamic (MD) simulation their stability, interaction, and conformation were analyzed. In summary, this study identified the most suitable compounds for targeting Mpro, and we recommend using these compounds as potential drug molecules against SARS-CoV-2 after follow up studies.     

Establishing an Analogue Based In Silico Pipeline in the Pursuit of Novel Inhibitory Scaffolds against the SARS Coronavirus 2 Papain-Like Protease

The ongoing coronavirus pandemic has been a burden on the worldwide population, with mass fatalities and devastating socioeconomic consequences. It has particularly drawn attention to the lack of approved small-molecule drugs to inhibit SARS coronaviruses. Importantly, lessons learned from the SARS outbreak of 2002–2004, caused by severe acute respiratory syndrome coronavirus 1 (SARS-CoV-1), can be applied to current drug discovery ventures. SARS-CoV-1 and SARS-CoV-2 both possess two cysteine proteases, the main protease (M^pro) and the papain-like protease (PL^pro), which play a significant role in facilitating viral replication, and are important drug targets. The non-covalent inhibitor, GRL-0617, which was found to inhibit replication of SARS-CoV-1, and more recently SARS-CoV-2, is the only PL^pro inhibitor co-crystallised with the recently solved SARS-CoV-2 PL^pro crystal structure. Therefore, the GRL-0617 structural template and pharmacophore features are instrumental in the design and development of more potent PL^pro inhibitors. In this work, we conducted scaffold hopping using GRL-0617 as a reference to screen over 339,000 ligands in the chemical space using the ChemDiv, MayBridge, and Enamine screening libraries. Twenty-four distinct scaffolds with structural and electrostatic similarity to GRL-0617 were obtained. These proceeded to molecular docking against PL^pro using the AutoDock tools. Of two compounds that showed the most favourable predicted binding affinities to the target site, as well as comparable protein-ligand interactions to GRL-0617, one was chosen for further analogue-based work. Twenty-seven analogues of this compound were further docked against the PL^pro, which resulted in two additional hits with promising docking profiles. Our in silico pipeline consisted of an integrative four-step approach: (1) ligand-based virtual screening (scaffold-hopping), (2) molecular docking, (3) an analogue search, and, (4) evaluation of scaffold drug-likeness, to identify promising scaffolds and eliminate those with undesirable properties. Overall, we present four novel, and lipophilic, scaffolds obtained from an exhaustive search of diverse and uncharted regions of chemical space, which may be further explored in vitro through structure-activity relationship (SAR) studies in the search for more potent inhibitors. Furthermore, these scaffolds were predicted to have fewer off-target interactions than GRL-0617. Lastly, to our knowledge, this work contains the largest ligand-based virtual screen performed against GRL-0617.     

Allosteric Inhibition of the SARS-CoV-2 Main Protease: Insights from Mass Spectrometry Based Assays**

The SARS-CoV-2 main protease (M^pro) cleaves along the two viral polypeptides to release non-structural proteins required for viral replication. M^Pro is an attractive target for antiviral therapies to combat the coronavirus-2019 disease. Here, we used native mass spectrometry to characterize the functional unit of M^pro. Analysis of the monomer/dimer equilibria reveals a dissociation constant of K_d=0.14±0.03 μM, indicating M^Pro has a strong preference to dimerize in solution. We characterized substrate turnover rates by following temporal changes in the enzyme-substrate complexes, and screened small molecules, that bind distant from the active site, for their ability to modulate activity. These compounds, including one proposed to disrupt the dimer, slow the rate of substrate processing by ≈35 %. This information, together with analysis of the x-ray crystal structures, provides a starting point for the development of more potent molecules that allosterically regulate M^Pro activity.     

From Repurposing to Redesign: Optimization of Boceprevir to Highly Potent Inhibitors of the SARS-CoV-2 Main Protease

The main protease (M^pro) of the betacoronavirus SARS-CoV-2 is an attractive target for the development of treatments for COVID-19. Structure-based design is a successful approach to discovering new inhibitors of the M^pro. Starting from crystal structures of the M^pro in complexes with the Hepatitis C virus NS3/4A protease inhibitors boceprevir and telaprevir, we optimized the potency of the alpha-ketoamide boceprevir against the M^pro by replacing its P1 cyclobutyl moiety by a γ-lactam as a glutamine surrogate. The resulting compound, MG-78, exhibited an IC₅₀ of 13 nM versus the recombinant M^pro, and similar potency was observed for its P1′ N-methyl derivative MG-131. Crystal structures confirmed the validity of our design concept. In addition to SARS-CoV-2 M^pro inhibition, we also explored the activity of MG-78 against the M^pro of the alphacoronavirus HCoV NL63 and against enterovirus 3C proteases. The activities were good (0.33 µM, HCoV-NL63 M^pro), moderate (1.45 µM, Coxsackievirus 3C^pro), and relatively poor (6.7 µM, enterovirus A71 3C^pro), respectively. The structural basis for the differences in activities was revealed by X-ray crystallo-graphy. We conclude that the modified boceprevir scaffold is suitable for obtaining high-potency inhibitors of the coronavirus M^pros but further optimization would be needed to target enterovirus 3C^pros efficiently.     

Inhibition of Cysteine Proteases by 6,6′-Dihydroxythiobinupharidine (DTBN) from Nuphar lutea

The specificity of inhibition by 6,6′-dihydroxythiobinupharidine (DTBN) on cysteine proteases was demonstrated in this work. There were differences in the extent of inhibition, reflecting active site structural-steric and biochemical differences. Cathepsin S (IC₅₀ = 3.2 μM) was most sensitive to inhibition by DTBN compared to Cathepsin B, L and papain (IC₅₀ = 1359.4, 13.2 and 70.4 μM respectively). DTBN is inactive for the inhibition of M^pro of SARS-CoV-2. Docking simulations suggested a mechanism of interaction that was further supported by the biochemical results. In the docking results, it was shown that the cysteine sulphur of Cathepsin S, L and B was in close proximity to the DTBN thiaspirane ring, potentially forming the necessary conditions for a nucleophilic attack to form a disulfide bond. Covalent docking and molecular dynamic simulations were performed to validate disulfide bond formation and to determine the stability of Cathepsins-DTBN complexes, respectively. The lack of reactivity of DTBN against SARS-CoV-2 M^pro was attributed to a mismatch of the binding conformation of DTBN to the catalytic binding site of M^pro. Thus, gradations in reactivity among the tested Cathepsins may be conducive for a mechanism-based search for derivatives of nupharidine against COVID-19. This could be an alternative strategy to the large-scale screening of electrophilic inhibitors.     

In silico exploration of binding potentials of anti SARS-CoV-1 phytochemicals against main protease of SARS-CoV-2

The phytochemicals can play complementary medicine compared to synthetic drugs considering their natural origin, safety, and low cost. Phytochemicals hold a key position for the expansion of drug development against corona viruses and need better consideration to the agents that have already been shown to display effective activity against various strains of corona viruses. In this study, we performed molecular docking studies on potential forty seven phytochemicals which are SARS-CoV-1 M^pro inhibitors to identify potential candidate against the main proteins of SARS-CoV-2. In Silico Molecular docking studies revealed that phytochemicals 16 (Broussoflavan A), 22 (Dieckol), 31 (Hygromycin B), 45 (Sinigrin) and 46 (Theaflavin-3,3′-digallate) exhibited excellent SARS-CoV-2 M^pro inhibitors. Furthermore, supported by Molecular dynamics (MD) simulation analysis such as Root Mean Square Deviation (RMSD), Root Mean Square Fluctuation (RMSF), Radius of gyration (Rg) and H-bond interaction analysis. We expect that our findings will provide designing principles for new corona virus strains and establish important frameworks for the future development of antiviral drugs.     

Ensemble Docking Coupled to Linear Interaction Energy Calculations for Identification of Coronavirus Main Protease (3CLpro) Non-Covalent Small-Molecule Inhibitors

SARS-CoV-2, or severe acute respiratory syndrome coronavirus 2, represents a new strain of Coronaviridae. In the closing 2019 to early 2020 months, the virus caused a global pandemic of COVID-19 disease. We performed a virtual screening study in order to identify potential inhibitors of the SARS-CoV-2 main viral protease (3CL^pro or M^pro). For this purpose, we developed a novel approach using ensemble docking high-throughput virtual screening directly coupled with subsequent Linear Interaction Energy (LIE) calculations to maximize the conformational space sampling and to assess the binding affinity of identified inhibitors. A large database of small commercial compounds was prepared, and top-scoring hits were identified with two compounds singled out, namely 1-[(R)-2-(1,3-benzimidazol-2-yl)-1-pyrrolidinyl]-2-(4-methyl-1,4-diazepan-1-yl)-1-ethanone and [({(S)-1-[(1H-indol-2-yl)methyl]-3-pyrrolidinyl}methyl)amino](5-methyl-2H-pyrazol-3-yl)formaldehyde. Moreover, we obtained a favorable binding free energy of the identified compounds, and using contact analysis we confirmed their stable binding modes in the 3CL^pro active site. These compounds will facilitate further 3CL^pro inhibitor design.     

The Inhibitory Effects of Plant Derivate Polyphenols on the Main Protease of SARS Coronavirus 2 and Their Structure–Activity Relationship

The main protease (M^pro) is a major protease having an important role in viral replication of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the novel coronavirus that caused the pandemic of 2020. Here, active M^pro was obtained as a 34.5 kDa protein by overexpression in E. coli BL21 (DE3). The optimal pH and temperature of M^pro were 7.5 and 37 °C, respectively. M^pro displayed a K_m value of 16 μM with Dabcyl-KTSAVLQ↓SGFRKME-Edans. Black garlic extract and 49 polyphenols were studied for their inhibitory effects on purified M^pro. The IC₅₀ values were 137 μg/mL for black garlic extract and 9–197 μM for 15 polyphenols. The mixtures of tannic acid with puerarin, daidzein, and/or myricetin enhanced the inhibitory effects on M^pro. The structure–activity relationship of these polyphenols revealed that the hydroxyl group in C3′, C4′, C5′ in the B-ring, C3 in the C-ring, C7 in A-ring, the double bond between C2 and C3 in the C-ring, and glycosylation at C8 in the A-ring contributed to inhibitory effects of flavonoids on M^pro.     

The Discovery of Small Allosteric and Active Site Inhibitors of the SARS-CoV-2 Main Protease via Structure-Based Virtual Screening and Biological Evaluation

The main protease enzyme (M^pro) of SARS-CoV-2 is one of the most promising targets for COVID-19 treatment. Accordingly, in this work, a structure-based virtual screening of 3.8 million ligand libraries was carried out. After rigorous filtering, docking, and post screening assessments, 78 compounds were selected for biological evaluation, 3 of which showed promising inhibition of the M^pro enzyme. The obtained hits (CB03, GR04, and GR20) had reasonable potencies with K_i values in the medium to high micromolar range. Interestingly, while our most potent hit, GR20, was suggested to act via a reversible covalent mechanism, GR04 was confirmed as a noncompetitive inhibitor that seems to be one of a kind when compared to the other allosteric inhibitors discovered so far. Moreover, all three compounds have small sizes (~300 Da) with interesting fittings in their relevant binding sites, and they possess lead-like characteristics that can introduce them as very attractive candidates for the future development of COVID-19 treatments.     

Non-Toxic Dimeric Peptides Derived from the Bothropstoxin-I Are Potent SARS-CoV-2 and Papain-like Protease Inhibitors

The COVID-19 outbreak has rapidly spread on a global scale, affecting the economy and public health systems throughout the world. In recent years, peptide-based therapeutics have been widely studied and developed to treat infectious diseases, including viral infections. Herein, the antiviral effects of the lysine linked dimer des-Cys¹¹, Lys¹²,Lys¹³-(pBthTX-I)₂K ((pBthTX-I)₂K)) and derivatives against SARS-CoV-2 are reported. The lead peptide (pBthTX-I)₂K and derivatives showed attractive inhibitory activities against SARS-CoV-2 (EC₅₀ = 28–65 µM) and mostly low cytotoxic effect (CC₅₀ > 100 µM). To shed light on the mechanism of action underlying the peptides’ antiviral activity, the Main Protease (M^pro) and Papain-Like protease (PL^pro) inhibitory activities of the peptides were assessed. The synthetic peptides showed PL^pro inhibition potencies (IC₅₀s = 1.0–3.5 µM) and binding affinities (K_d = 0.9–7 µM) at the low micromolar range but poor inhibitory activity against M^pro (IC₅₀ > 10 µM). The modeled binding mode of a representative peptide of the series indicated that the compound blocked the entry of the PL^pro substrate toward the protease catalytic cleft. Our findings indicated that non-toxic dimeric peptides derived from the Bothropstoxin-I have attractive cellular and enzymatic inhibitory activities, thereby suggesting that they are promising prototypes for the discovery and development of new drugs against SARS-CoV-2 infection.     

Identification of SARS-CoV-2 Main Protease Inhibitors from a Library of Minor Cannabinoids by Biochemical Inhibition Assay and Surface Plasmon Resonance Characterized Binding Affinity

The replication of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is mediated by its main protease (M^pro), which is a plausible therapeutic target for coronavirus disease 2019 (COVID-19). Although numerous in silico studies reported the potential inhibitory effects of natural products including cannabis and cannabinoids on SARS-CoV-2 M^pro, their anti-M^pro activities are not well validated by biological experimental data. Herein, a library of minor cannabinoids belonging to several chemotypes including tetrahydrocannabinols, cannabidiols, cannabigerols, cannabichromenes, cannabinodiols, cannabicyclols, cannabinols, and cannabitriols was evaluated for their anti-M^pro activity using a biochemical assay. Additionally, the binding affinities and molecular interactions between the active cannabinoids and the M^pro protein were studied by a biophysical technique (surface plasmon resonance; SPR) and molecular docking, respectively. Cannabinoids tetrahydrocannabutol and cannabigerolic acid were the most active M^pro inhibitors (IC₅₀ = 3.62 and 14.40 μM, respectively) and cannabigerolic acid had a binding affinity KD=2.16×10−4 M). A preliminary structure and activity relationship study revealed that the anti-Mpro effects of cannabinoids were influenced by the decarboxylation of cannabinoids and the length of cannabinoids’ alkyl side chain. Findings from the biochemical, biophysical, and computational assays support the growing evidence of cannabinoids’ inhibitory effects on SARS-CoV-2 M^pro.