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.     

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.     

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.     

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.     

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.     

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.     

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.     

Synthesis and Biochemical Evaluation of 8H-Indeno[1,2-d]thiazole Derivatives as Novel SARS-CoV-2 3CL Protease Inhibitors

The COVID-19 pandemic caused by SARS-CoV-2 is a global burden on human health and economy. The 3-Chymotrypsin-like cysteine protease (3CL^pro) becomes an attractive target for SARS-CoV-2 due to its important role in viral replication. We synthesized a series of 8H-indeno[1,2-d]thiazole derivatives and evaluated their biochemical activities against SARS-CoV-2 3CL^pro. Among them, the representative compound 7a displayed inhibitory activity with an IC₅₀ of 1.28 ± 0.17 μM against SARS-CoV-2 3CL^pro. Molecular docking of 7a against 3CL^pro was performed and the binding mode was rationalized. These preliminary results provide a unique prototype for the development of novel inhibitors against SARS-CoV-2 3CL^pro.     

In Silico Mining of Terpenes from Red-Sea Invertebrates for SARS-CoV-2 Main Protease (Mpro) Inhibitors

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent for the COVID-19 pandemic, which generated more than 1.82 million deaths in 2020 alone, in addition to 83.8 million infections. Currently, there is no antiviral medication to treat COVID-19. In the search for drug leads, marine-derived metabolites are reported here as prospective SARS-CoV-2 inhibitors. Two hundred and twenty-seven terpene natural products isolated from the biodiverse Red-Sea ecosystem were screened for inhibitor activity against the SARS-CoV-2 main protease (M^pro) using molecular docking and molecular dynamics (MD) simulations combined with molecular mechanics/generalized Born surface area binding energy calculations. On the basis of in silico analyses, six terpenes demonstrated high potency as M^pro inhibitors with ΔG_binding ≤ −40.0 kcal/mol. The stability and binding affinity of the most potent metabolite, erylosides B, were compared to the human immunodeficiency virus protease inhibitor, lopinavir. Erylosides B showed greater binding affinity towards SARS-CoV-2 M^pro than lopinavir over 100 ns with ΔG_binding values of −51.9 vs. −33.6 kcal/mol, respectively. Protein–protein interactions indicate that erylosides B biochemical signaling shares gene components that mediate severe acute respiratory syndrome diseases, including the cytokine- and immune-signaling components BCL2L1, IL2, and PRKC. Pathway enrichment analysis and Boolean network modeling were performed towards a deep dissection and mining of the erylosides B target–function interactions. The current study identifies erylosides B as a promising anti-COVID-19 drug lead that warrants further in vitro and in vivo testing.     

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.     

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.     

Selective covalent targeting of SARS-CoV-2 main protease by enantiopure chlorofluoroacetamide

The coronavirus disease 2019 (COVID-19) pandemic has necessitated the development of antiviral agents against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The main protease (M^pro) is a promising target for COVID-19 treatment. Here, we report an irreversible SARS-CoV-2 M^pro inhibitor possessing chlorofluoroacetamide (CFA) as a warhead for the covalent modification of M^pro. Ugi multicomponent reaction using chlorofluoroacetic acid enabled the rapid synthesis of dipeptidic CFA derivatives that identified 18 as a potent inhibitor of SARS-CoV-2 M^pro. Among the four stereoisomers, (R,R)-18 exhibited a markedly higher inhibitory activity against M^pro than the other isomers. Reaction kinetics and computational docking studies suggest that the R configuration of the CFA warhead is crucial for the rapid covalent inhibition of M^pro. Our findings highlight the prominent influence of the CFA chirality on the covalent modification of proteinous cysteines and provide the basis for improving the potency and selectivity of CFA-based covalent inhibitors.

Chlorofluoroacetamide (CFA) was used as the warhead for covalent targeting of SARS-CoV-2 M^pro. The chirality at CFA showed marked influence on inhibitory activity, suggesting stereospecific activation of CFA for cysteine modification in the protein.     

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.     

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.     

Site mapping and small molecule blind docking reveal a possible target site on the SARS-CoV-2 main protease dimer interface

The SARS-CoV-2 virus is causing COVID-19 resulting in an ongoing pandemic with serious health, social, and economic implications. Much research is focused in repurposing or identifying new small molecules which may interact with viral or host-cell molecular targets. An important SARS-CoV-2 target is the main protease (M^pro), and the peptidomimetic α-ketoamides represent prototypical experimental inhibitors. The protease is characterised by the dimerization of two monomers each which contains the catalytic dyad defined by Cys¹⁴⁵ and His⁴¹ residues (active site). Dimerization yields the functional homodimer. Here, our aim was to investigate small molecules, including lopinavir and ritonavir, α-ketoamide 13b, and ebselen, for their ability to interact with the M^pro. The sirtuin 1 agonist SRT1720 was also used in our analyses. Blind docking to each monomer individually indicated preferential binding of the ligands in the active site. Site-mapping of the dimeric protease indicated a highly reactive pocket in the dimerization region at the domain III apex. Blind docking consistently indicated a strong preference of ligand binding in domain III, away from the active site. Molecular dynamics simulations indicated that ligands docked both to the active site and in the dimerization region at the apex, formed relatively stable interactions. Overall, our findings do not obviate the superior potency with respect to inhibition of protease activity of covalently-linked inhibitors such as α-ketoamide 13b in the M^pro active site. Nevertheless, along with those from others, our findings highlight the importance of further characterisation of the M^pro active site and any potential allosteric sites.     

Structural Insight into the Binding of Cyanovirin-N with the Spike Glycoprotein,Mpro and PLpro of SARS-CoV-2: Protein–Protein Interactions,Dynamics Simulations and Free Energy Calculations

The emergence of COVID-19 continues to pose severe threats to global public health. The pandemic has infected over 171 million people and claimed more than 3.5 million lives to date. We investigated the binding potential of antiviral cyanobacterial proteins including cyanovirin-N, scytovirin and phycocyanin with fundamental proteins involved in attachment and replication of SARS-CoV-2. Cyanovirin-N displayed the highest binding energy scores (−16.8 ± 0.02 kcal/mol, −12.3 ± 0.03 kcal/mol and −13.4 ± 0.02 kcal/mol, respectively) with the spike protein, the main protease (M^pro) and the papainlike protease (PL^pro) of SARS-CoV-2. Cyanovirin-N was observed to interact with the crucial residues involved in the attachment of the human ACE2 receptor. Analysis of the binding affinities calculated employing the molecular mechanics-Poisson–Boltzmann surface area (MM-PBSA) approach revealed that all forms of energy, except the polar solvation energy, favourably contributed to the interactions of cyanovirin-N with the viral proteins. With particular emphasis on cyanovirin-N, the current work presents evidence for the potential inhibition of SARS-CoV-2 by cyanobacterial proteins, and offers the opportunity for in vitro and in vivo experiments to deploy the cyanobacterial proteins as valuable therapeutics against COVID-19.     

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.     

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.