Identification of 14 Known Drugs as Inhibitors of the Main Protease of SARS-CoV-2

A strategy to optimize the peptide-based inhibitors against different mutants of the spike protein of SARS-CoV-2

The SARS-CoV-2 virus has caused high-priority health concerns at a global level. Vaccines have stalled the proliferation of viruses to some extent. Yet, the emergence of newer, potentially more infectious, and dangerous mutants such as Delta and Omicron are among the major challenges in finding a more permanent solution for this pandemic. The effectiveness of antivirals Molnupiravir and Paxlovid, authorized for emergency use by the FDA, are yet to be assessed on a larger population. Patients with a high risk of disease progression have received treatment with antibody-cocktail. Most of the mutations leading to the new lineage of SARS-CoV-2 are found in the spike protein of this virus that plays a key role in facilitating host entry. The current study has investigated how to modify a promising peptide-based inhibitor of spike protein, LCB3, against common mutations, N501Y and K417N in the target protein so that it retains its efficacy against the spike protein. LCB3 being a prototype for protein-based inhibitors is an ideal testing system to learn about protein-based inhibitors. This study proposes the substitutions of amino acid residues of LCB3 inhibitor using a structure-based approach that considers free energy decomposition of residues, the distance between atoms, and the interaction among amino acids. The binding free energy calculations suggest a possible improvement in the binding affinity of existing inhibitor LCB3 to the mutant forms of the S-protein using simple substitutions at specific positions of the inhibitor. This approach, being general, can be used in different inhibitors and other mutations and help in fighting against SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

Potential COVID-19 Therapies from Computational Repurposing of Drugs and Natural Products against the SARS-CoV-2 Helicase

Repurposing of existing drugs is a rapid way to find potential new treatments for SARS-CoV-2. Here, we applied a virtual screening approach using Autodock Vina and molecular dynamic simulation in tandem to screen and calculate binding energies of repurposed drugs against the SARS-CoV-2 helicase protein (non-structural protein nsp13). Amongst the top hits from our study were antivirals, antihistamines, and antipsychotics, plus a range of other drugs. Approximately 30% of our top 87 hits had published evidence indicating in vivo or in vitro SARS-CoV-2 activity. Top hits not previously reported to have SARS-CoV-2 activity included the antiviral agents, cabotegravir and RSV-604; the NK1 antagonist, aprepitant; the trypanocidal drug, aminoquinuride; the analgesic, antrafenine; the anticancer intercalator, epirubicin; the antihistamine, fexofenadine; and the anticoagulant, dicoumarol. These hits from our in silico SARS-CoV-2 helicase screen warrant further testing as potential COVID-19 treatments.

Synthesis, crystal structure, DFT and molecular docking studies of N-acetyl-2,4-[diaryl-3-azabicyclo[3.3.1]nonan-9-yl]-9-spiro-4'-acetyl-2'-(acetylamino)-4',9-dihydro-[1',3',4']-thiadiazoles: A potential SARS-nCoV-2 Mpro (COVID-19) inhibitor

In this paper, we describe the synthesis and crystal structure analysis of N-acetyl-2,4-[diphenyl-3-azabicyclo[3.3.1]nonan-9-yl]-9-spiro-4′-acetyl-2′-(acetylamino)-4′,9-dihydro-[1′,3′,4′]-thiadiazole (3a) and N-acetyl- 2,4-[bis(p-methoxyphenyl)-3-azabicyclo[3.3.1]nonan-9-yl]-9-spiro-4′-acetyl-2′-(acetylamino)-4′,9-dihydro-[1′,3′,4′]-thiadiazole (3b). The title compounds 3a and 3b are characterized by 1D NMR and single crystal x-ray diffraction analysis. Non-covalent interactions in a molecule were identified by Hirshfeld surface (d_norm contacts and 2D fingerprint plot) analysis. In addition, the existence of chalcogen bond (S•••O bond) in the molecular structures (3a and 3b) are described by NCI-RDG and QTAIM analysis. NBO analysis is employed to describe the orbital interactions and electron transfer between sulfur and oxygen atoms. Molecular docking is carried out for compounds 3a and 3b with COVID-19 viral protein SARS-nCoV-2 M^pro (PDB ID: 6LU7).

Primer for Designing Main Protease (Mpro) Inhibitors of SARS-CoV-2

The COVID-19 outbreak has been devastating, with hundreds of millions of infections and millions of deaths reported worldwide. In response, the application of structure-activity relationships (SAR) upon experimentally validated inhibitors of SARS-CoV-2 main protease (Mpro) may provide an avenue for the identification of new lead compounds active against COVID-19. Upon the basis of information gleaned from a combination of reported crystal structures and the docking of experimentally validated inhibitors, four "rules" for designing potent Mpro inhibitors have been proposed. The aim here is to guide medicinal chemists toward the most probable hits and to provide guidance on repurposing available structures as Mpro inhibitors. Experimental examination of our own previously reported inhibitors using the four "rules" identified a potential lead compound, the cathepsin inhibitor GB111-NH2, that was 2.3 times more potent than SARS-CoV-2 Mpro inhibitor N3.

Discovery of 2-Thiobenzimidazoles as Noncovalent Inhibitors of SARS-CoV-2 Main Protease

The discovery of antiviral agents against SARS-CoV-2 is an important step toward ending the COVID-19 pandemic and to tackle future outbreaks. In this context, the main protease (M^pro) represents an ideal target for developing coronavirus antivirals, being conserved among different strains and essential for survival. In this work, using in silico tools, we created and validated a docking protocol able to predict binders to the catalytic site of M^pro. The following structure-based virtual screening of a subset of the ZINC library (over 4.3 million unique structures), led to the identification of a hit compound having a 2-thiobenzimidazole scaffold. The inhibitory activity was confirmed using a FRET-based proteolytic assay against recombinant M^pro. Structure-activity relationships were obtained with the synthesis of a small library of analogs, guided by the analysis of the docking pose. Our efforts led to the identification of a micromolar M^pro inhibitor (IC₅₀ = 14.9 µM) with an original scaffold possessing ideal drug-like properties (predicted using the QikProp function) and representing a promising lead for the development of a novel class of coronavirus antivirals.

Quercetin and luteolin are single-digit micromolar inhibitors of the SARS-CoV-2 RNA-dependent RNA polymerase

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has rapidly become a global health pandemic. Among the viral proteins, RNA-dependent RNA polymerase (RdRp) is responsible for viral genome replication and has emerged as one of the most promising targets for pharmacological intervention against SARS-CoV-2. To this end, we experimentally tested luteolin and quercetin for their ability to inhibit the RdRp enzyme. These two compounds are ancestors of flavonoid natural compounds known for a variety of basal pharmacological activities. Luteolin and quercetin returned a single-digit IC₅₀ of 4.6 µM and 6.9 µM, respectively. Then, through dynamic docking simulations, we identified possible binding modes of these compounds to a recently published cryo-EM structure of RdRp. Collectively, these data indicate that these two compounds are a valid starting point for further optimization and development of a new class of RdRp inhibitors to treat SARS-CoV-2 and potentially other viral infections.

Exploiting protease activation for therapy

Proteases have crucial roles in homeostasis and disease; and protease inhibitors and recombinant proteases in enzyme replacement therapy have become key therapeutic applications of protease biology across several indications. This review briefly summarises therapeutic approaches based on protease activation and focuses on how recent insights into the spatial and temporal control of the proteolytic activation of growth factors and interleukins are leading to unique strategies for the discovery of new medicines. In particular, two emerging areas are covered: the first is based on antibody therapies that target the process of proteolytic activation of the pro-form of proteins rather than their mature form; the second covers a potentially new class of biopharmaceuticals using engineered, proteolytically activable and initially inactive pro-forms of antibodies or effector proteins to increase specificity and improve the therapeutic window.

[Nirmatrelvir plus ritonavir (Paxlovid) a potent SARS-CoV-2 3CLpro protease inhibitor combination]

All coronavirus, including SARS-CoV-2, encode two proteases needed for the processing of PP1A and PP1AB polyproteins. The main protease 3CL (chemotripsine-like) gives rise to the formation of NSP11/16 proteins. The 3CL protease has been constituted as one of the possible therapeutic targets for the development of antiviral drugs against SARS-COV-2 due to its highly conserved sequence and structure among all coronaviruses. During the SARS-COV-1 pandemic, a hydroxymethyl ketone derivative (PF-00835231) was identified with an intense inhibitory activity against the 3CL protease. Subsequent chemical modifications gave rise to derivative PF-07321332 (nirmatrelvir) which has shown a high antiviral efficacy against SARS-COV-2. The company's data indicate that it is capable of reducing 89% the risk of hospitalization and death of patients infected with hardly adverse effects. Its effectiveness improves if it is administered orally in the first 24-48 hours and the duration of treatment has been established between 3-5 days. The commercial form has been associated with the antiviral ritonavir that has shown the metabolism of nirmatrelvir, lengthening its average life. This antiviral would be effective against current and future viral variants, since 3CL is not modified in them. The FDA approved this antiviral in November 2021 and EMA is in the final evaluation phase.

The effect of various compounds on the COVID mechanisms, from chemical to molecular aspects

The novel coronavirus that caused COVID-19 pandemic is SARS-CoV-2. Although various vaccines are currently being used to prevent the disease's severe consequences, there is still a need for medications for those who become infected. The SARS-CoV-2 has a variety of proteins that have been studied extensively since the virus's advent. In this review article, we looked at chemical to molecular aspects of the various structures studied that have pharmaceutical activity and attempted to find a link between drug activity and compound structure. For example, designing of the compounds which bind to the allosteric site and modify hydrogen bonds or the salt bridges can disrupt SARS-CoV2 RBD–ACE2 complex. It seems that quaternary ammonium moiety and quinolin-1-ium structure could act as a negative allosteric modulator to reduce the tendency between spike-ACE2. Pharmaceutical structures with amino heads and hydrophobic tails can block envelope protein to prevent making mature SARS-CoV-2. Also, structures based on naphthalene pharmacophores or isosteres can form a strong bond with the PLpro and form a π-π and the Mpro's active site can be occupied by octapeptide compounds or linear compounds with a similar fitting ability to octapeptide compounds. And for protein RdRp, it is critical to consider pH and pKa so that pKa regulation of compounds to comply with patients is very effective, thus, the presence of tetrazole, phenylpyrazole groups, and analogs of pyrophosphate in the designed drugs increase the likelihood of the RdRp active site inhibition. Finally, it can be deduced that designing hybrid drug molecules along with considering the aforementioned characteristics would be a suitable approach for developing medicines in order to accurate targeting and complete inhibition this virus.

Nucleic Acids as Biotools at the Interface between Chemistry and Nanomedicine in the COVID-19 Era

The recent development of mRNA vaccines against the SARS-CoV-2 infection has turned the spotlight on the potential of nucleic acids as innovative prophylactic agents and as diagnostic and therapeutic tools. Until now, their use has been severely limited by their reduced half-life in the biological environment and the difficulties related to their transport to target cells. These limiting aspects can now be overcome by resorting to chemical modifications in the drug and using appropriate nanocarriers, respectively. Oligonucleotides can interact with complementary sequences of nucleic acid targets, forming stable complexes and determining their loss of function. An alternative strategy uses nucleic acid aptamers that, like the antibodies, bind to specific proteins to modulate their activity. In this review, the authors will examine the recent literature on nucleic acids-based strategies in the COVID-19 era, focusing the attention on their applications for the prophylaxis of COVID-19, but also on antisense- and aptamer-based strategies directed to the diagnosis and therapy of the coronavirus pandemic.

How Ethanolic Disinfectants Disintegrate Coronavirus Model Membranes: A Dissipative Particle Dynamics Simulation Study

We have developed dissipative particle dynamics models for pure dipalmitoylphosphatidylcholine (DPPC), dioleoylphosphatidylcholine (DOPC), and dimyristoylphosphatidylcholine (DMPC) as well as their binary and ternary mixed membranes, as coronavirus model membranes. The stabilities of pure and mixed membranes, surrounded by aqueous solutions containing up to 70 mol % ethanol (alcoholic disinfectants), have been investigated at room temperature. We found that aqueous solutions containing 5-10 mol % ethanol already have a significant weakening effect on the pure and mixed membranes. The magnitude of the effect depends on the membrane composition and the ethanol concentration. Ethanol permeabilizes the membrane, causing its lateral swelling and thickness shrinking and reducing the orientational order of the hydrocarbon tail of the bilayer. The free energy barrier for the permeation of ethanol in the bilayers is considerably reduced by the ethanol uptake. The rupture-critical ethanol concentrations causing the membrane failure are 20.7, 27.5, and 31.7 mol % in the aqueous phase surrounding pure DMPC, DOPC, and DPPC membranes, respectively. Characterizing the failure of lipid membranes by a machine-learning neural network framework, we found that all mixed binary and/or ternary membranes disrupt when immersed in an aqueous solution containing a rupture-critical ethanol concentration, ranging from 20.7 to 31.7 mol %, depending on the composition of the membrane; the DPPC-rich membranes are more intact, while the DMPC-rich membranes are least intact. Due to the tight packing of long, saturated hydrocarbon tails in DPPC, increasing the DPPC content of the mixed membrane increases its stability against the disinfectant. At high DPPC concentrations, where the DOPC and DMPC molecules are confined between the DPPC lipids, the ordered hydrocarbon tails of DPPC also induce order in the DOPC and DMPC molecules and, hence, stabilize the membrane more. Our simulations on pure and mixed membranes of a diversity of compositions reveal that a maximum ethanol concentration of 32 mol % (55 wt %) in the alcohol-based disinfectants is enough to disintegrate any membrane composed of these three lipids.

Anti-Fungal Drug Anidulafungin Inhibits SARS-CoV-2 Spike-Induced Syncytia Formation by Targeting ACE2-Spike Protein Interaction

Drug repositioning continues to be the most effective, practicable possibility to treat COVID-19 patients. The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus enters target cells by binding to the ACE2 receptor via its spike (S) glycoprotein. We used molecular docking-based virtual screening approaches to categorize potential antagonists, halting ACE2-spike interactions by utilizing 450 FDA-approved chemical compounds. Three drug candidates (i.e., anidulafungin, lopinavir, and indinavir) were selected, which show high binding affinity toward the ACE2 receptor. The conformational stability of selected docked complexes was analyzed through molecular dynamics (MD) simulations. The MD simulation trajectories were assessed and monitored for ACE2 deviation, residue fluctuation, the radius of gyration, solvent accessible surface area, and free energy landscapes. The inhibitory activities of the selected compounds were eventually tested in-vitro using Vero and HEK-ACE2 cells. Interestingly, besides inhibiting SARS-CoV-2 S glycoprotein induced syncytia formation, anidulafungin and lopinavir also blocked S-pseudotyped particle entry into target cells. Altogether, anidulafungin and lopinavir are ranked the most effective among all the tested drugs against ACE2 receptor-S glycoprotein interaction. Based on these findings, we propose that anidulafungin is a novel potential drug targeting ACE2, which warrants further investigation for COVID-19 treatment.

Molecular insights to the binding interactions of APNS containing HIV-protease inhibitors against SARS-CoV-2 Mpro: an in silico approach towards drug repurposing

SARS-CoV-2 M^pro is one of the most vital enzymes of the new coronavirus-2 (SARS-CoV-2) and is a crucial target for drug discovery. Unfortunately, there is not any potential drugs available to combat the action of SARS-CoV-2 M^pro. Based on the reports HIV-protease inhibitors can be applied against the SARS by targeting the SARS-CoV-1 M^pro, we have chosen few clinically trialed experimental and allophenylnorstatine (APNS) containing HIV-protease inhibitors (JE-2147, JE-533, KNI-227, KNI-272 & KNI-1931), to examine their binding affinities with SARS-CoV-2 M^pro and to assess their potential to check for a possible drug candidate against the protease. Here, we have chosen a methodology to understand the binding mechanism of these five inhibitors to SARS-CoV-2 M^pro by merging molecular docking, molecular dynamics (MD) simulation and MM-PBSA based free energy calculations. Our estimations disclose that JE-2147 is highly effective (ΔG_Bind = -28.31 kcal/mol) due to an increased favorable van der Waals (ΔE_vdw) interactions and decreased solvation (ΔG_solv) energies between the inhibitor and viral protease. JE-2147 shows a higher level of interactions as compared to JE-533 (-6.85 kcal/mol), KNI-227 (-18.36 kcal/mol), KNI-272 (-15.69 kcal/mol) and KNI-1931 (-21.59 kcal/mol) against SARS-CoV-2 M^pro. Binding contributions of important residues (His41, Met49, Cys145, His164, Met165, Glu166, Pro168, Gln189, etc.) from the active site or near the active site regions with ≥1.0 kcal/mol suggest a potent binding of the inhibitors. It is anticipated that the current study of binding interactions of these APNS containing inhibitors can pitch some valuable insights to design the significantly effective anti-SARS-CoV-2 M^pro drugs.Communicated by Ramaswamy H. Sarma.

Validation and invalidation of SARS-CoV-2 main protease inhibitors using the Flip-GFP and Protease-Glo luciferase assays

SARS-CoV-2 main protease (M^pro) is one of the most extensively exploited drug targets for COVID-19. Structurally disparate compounds have been reported as M^pro inhibitors, raising the question of their target specificity. To elucidate the target specificity and the cellular target engagement of the claimed M^pro inhibitors, we systematically characterize their mechanism of action using the cell-free FRET assay, the thermal shift-binding assay, the cell lysate Protease-Glo luciferase assay, and the cell-based FlipGFP assay. Collectively, our results have shown that majority of the M^pro inhibitors identified from drug repurposing including ebselen, carmofur, disulfiram, and shikonin are promiscuous cysteine inhibitors that are not specific to M^pro, while chloroquine, oxytetracycline, montelukast, candesartan, and dipyridamole do not inhibit M^pro in any of the assays tested. Overall, our study highlights the need of stringent hit validation at the early stage of drug discovery.

Computationally driven discovery of SARS-CoV-2 Mpro inhibitors: from design to experimental validation

We report a fast-track computationally driven discovery of new SARS-CoV-2 main protease (M^pro) inhibitors whose potency ranges from mM for the initial non-covalent ligands to sub-μM for the final covalent compound (IC₅₀ = 830 ± 50 nM). The project extensively relied on high-resolution all-atom molecular dynamics simulations and absolute binding free energy calculations performed using the polarizable AMOEBA force field. The study is complemented by extensive adaptive sampling simulations that are used to rationalize the different ligand binding poses through the explicit reconstruction of the ligand–protein conformation space. Machine learning predictions are also performed to predict selected compound properties. While simulations extensively use high performance computing to strongly reduce the time-to-solution, they were systematically coupled to nuclear magnetic resonance experiments to drive synthesis and for in vitro characterization of compounds. Such a study highlights the power of in silico strategies that rely on structure-based approaches for drug design and allows the protein conformational multiplicity problem to be addressed. The proposed fluorinated tetrahydroquinolines open routes for further optimization of M^pro inhibitors towards low nM affinities.

The dominant binding mode of the QUB-00006-Int-07 main protease inhibitor during absolute binding free energy simulations.

Association of Amlodipine with the Risk of In-Hospital Death in Patients with COVID-19 and Hypertension: A Reanalysis on 184 COVID-19 Patients with Hypertension

Association between calcium channel blockers (CCBs) or functional inhibitors of acid sphingomyelinase (FIASMAs) use and decreased mortality in people with COVID-19 has been reported in recent studies. Since amlodipine is both a CCB and a FIASMA, the aim of this study was to investigate the association between chronic amlodipine use and the survival of people with hypertension infected with COVID-19. This retrospective cohort study used data extracted from the medical records of adult inpatients with hypertension and laboratory-confirmed COVID-19 between 1 March 2020 and 31 August 2020 with definite outcomes (discharged from hospital or deceased) from Erasme Hospital (Brussels, Belgium). We re-analyzed the data of the retrospective cohort study using only the 184 patients (103 males, 81 females) with a mean age of 69.54 years (SD = 14.6) with hypertension. The fifty-five participants (29.9%) receiving a chronic prescription of amlodipine were compared with the 129 patients who did not receive a chronic prescription of amlodipine. Univariate and multivariate logistic regressions were used to explore the relationships between mortality and sex, age, comorbidities, smoking, and amlodipine status. Out of the 184 participants, 132 (71.7%) survived and 52 (28.3%) died. The mortality rates were, respectively, 12.73% (n = 7) and 34.88% (n = 45) for the amlodipine and non-amlodipine groups. Multivariate logistic regression was significant (Chi square (5) = 29.11; p < 0.0001). Chronic kidney disease and malignant neoplasm were significant predictors as well as amlodipine status. For chronic kidney disease and malignant neoplasm, the odds ratio with 95% confidence interval (95% CI) were, respectively, 2.16 (95% CI: 1.04−4.5; p = 0.039) and 2.46 (95% CI: 1.01−6.01; p = 0.047). For amlodipine status the odds ratio was 0.29 (95% CI: 0.11−0.74; p = 0.009). The result of the present study suggests that amlodipine may be associated with reduced mortality in people with hypertension infected with COVID-19. Further research and randomized clinical trials are needed to confirm the potential protective effect of amlodipine in people with hypertension infected with COVID-19.

Inhibition of the main protease of SARS-CoV-2 (Mpro) by repurposing/designing drug-like substances and utilizing nature's toolbox of bioactive compounds

The emergence of the Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) has resulted in a long pandemic, with numerous cases and victims worldwide and enormous consequences on social and economic life. Although vaccinations have proceeded and provide a valuable shield against the virus, the approved drugs are limited and it is crucial that further ways to combat infection are developed, that can also act against potential mutations. The main protease (Mpro) of the virus is an appealing target for the development of inhibitors, due to its importance in the viral life cycle and its high conservation among different coronaviruses. Several compounds have shown inhibitory potential against Mpro, both in silico and in vitro, with few of them also having entered clinical trials. These candidates include: known drugs that have been repurposed, molecules specifically designed based on the natural substrate of the protease or on structural moieties that have shown high binding affinity to the protease active site, as well as naturally derived compounds, either isolated or in plant extracts. The aim of this work is to collectively present the results of research regarding Mpro inhibitors to date, focusing on the function of the compounds founded by in silico simulations and further explored by in vitro and in vivo assays. Creating an extended portfolio of promising compounds that may block viral replication by inhibiting Mpro and by understanding involved structure-activity relationships, could provide a basis for the development of effective solutions against SARS-CoV-2 and future related outbreaks.

Evaluation of COVID-19 protease and HIV inhibitors interactions

The epidemic of the novel coronavirus disease (COVID-19) that started in 2019 has evoked an urgent demand for finding new potential therapeutic agents. In this study, we performed a molecular docking of anti-HIV drugs to refine HIV protease inhibitors and nucleotide analogues to target COVID-19. The evaluation was based on docking scores calculated by AutoDock Vina and top binding poses were analyzed. Our results suggested that lopinavir, darunavir, atazanavir, remdesivir, and tipranavir have the best binding affinity for the 3-chymotrypsin-like protease of COVID-19. The comparison of the binding sites of three drugs, namely, darunavir, atazanavir and remdesivir, showed an overlap region of the protein pocket. Our study showed a strong affinity between lopinavir, darunavir, atazanavir, tipranavir and COVID-19 protease. However, their efficacy should be confirmed by in vitro studies since there are concerns related to interference with their active sites.

Molecular simulations of proteins: From simplified physical interactions to complex biological phenomena

Molecular dynamics simulation is the most popular computational technique for investigating the structural and dynamical behaviour of proteins, in search of the molecular basis of their function. Far from being a completely settled field of research, simulations are still evolving to best capture the essential features of the atomic interactions that govern a protein's inner motions. Modern force fields are becoming increasingly accurate in providing a physical description adequate to this purpose, and allow us to model complex biological systems under fairly realistic conditions. Furthermore, the use of accelerated sampling techniques is improving our access to the observation of progressively larger molecular structures, longer time scales, and more hidden functional events. In this review, the basic principles of molecular dynamics simulations and a number of key applications in the area of protein science are summarized, and some of the most important results are discussed. Examples include the study of the structure, dynamics and binding properties of 'difficult' targets, such as intrinsically disordered proteins and membrane receptors, and the investigation of challenging phenomena like hydration-driven processes and protein aggregation. The findings described provide an overall picture of the current state of this research field, and indicate new perspectives on the road ahead to the upcoming future of molecular simulations.

A new inactive conformation of SARS-CoV-2 main protease

The SARS‐CoV‐2 main protease (M^pro) has a pivotal role in mediating viral genome replication and transcription of the coronavirus, making it a promising target for drugs against the COVID‐19 pandemic. Here, a crystal structure is presented in which M^pro adopts an inactive state that has never been observed before, called new‐inactive. It is shown that the oxyanion loop, which is involved in substrate recognition and enzymatic activity, adopts a new catalytically incompetent conformation and that many of the key interactions of the active conformation of the enzyme around the active site are lost. Solvation/desolvation energetic contributions play an important role in the transition from the inactive to the active state, with Phe140 moving from an exposed to a buried environment and Asn142 moving from a buried environment to an exposed environment. In new‐inactive M^pro a new cavity is present near the S2′ subsite, and the N‐terminal and C‐terminal tails, as well as the dimeric interface, are perturbed, with partial destabilization of the dimeric assembly. This novel conformation is relevant both for comprehension of the mechanism of action of M^pro within the catalytic cycle and for the successful structure‐based drug design of antiviral drugs.

Molecular Docking and Molecular Dynamics Simulations Discover Curcumin Analogue as a Plausible Dual Inhibitor for SARS-CoV-2

Recently, the world has been witnessing a global pandemic with no effective therapeutics yet, while cancer continues to be a major disease claiming many lives. The natural compound curcumin is bestowed with multiple medicinal applications in addition to demonstrating antiviral and anticancer activities. In order to elucidate the impact of curcumin on COVID-19 and cancer, the current investigation has adapted several computational techniques to unfold its possible inhibitory activity. Accordingly, curcumin and similar compounds and analogues were retrieved and assessed for their binding affinities at the binding pocket of SARS-CoV-2 main protease and DDX3. The best binding pose was escalated to molecular dynamics simulation (MDS) studies to assess the time dependent stability. Our findings have rendered one compound that has demonstrated good molecular dock score complemented by key residue interactions and have shown stable MDS results inferred by root mean square deviation (RMSD), radius of gyration (Rg), binding mode, hydrogen bond interactions, and interaction energy. Essential dynamics results have shown that the systemadapts minimum energy conformation to attain a stable state. The discovered compound (curA) could act as plausible inhibitor against SARS-CoV-2 and DDX3. Furthermore, curA could serve as a chemical scaffold for designing and developing new compounds.

Phytochemicals of Euphorbia hirta L. and Their Inhibitory Potential Against SARS-CoV-2 Main Protease

Euphorbia hirta L. is a medicinal plant widely used in the Philippines and across tropical Asia against various diseases, including respiratory disorders. In this study, the phytochemical components of E. hirta were investigated in silico for their potential to inhibit the severe acute respiratory syndrome-coronavirus-2 main protease (SARS-CoV-2 Mpro), a coronavirus disease 2019 (COVID-19) drug target that plays a critical role in the infection process of SARS-CoV-2. Phytochemical mining in tandem with virtual screening (PM-VS) was the strategy implemented in this study, which allows efficient preliminary in silico assessment of the COVID-19 therapeutic potential of the reported phytochemicals from the plant. The main rationale for considering E. hirta in the investigation was its reported efficacy against respiratory disorders. It is very promising to investigate the phytochemicals of E. hirta for their potential efficacy against diseases, such as COVID-19, that also target the respiratory system. A total of 298 E. hirta phytochemicals were comprehensively collected from the scientific literature. One hundred seventy of these phytochemicals were computed through molecular docking and were shown to have comparable or better binding properties (promising inhibitors) toward SARS-CoV-2 Mpro than known in vitro inhibitors. In connection to our previous work considering different medicinal plants, antiviral compounds were also rediscovered from the phytochemical composition of E. hirta. This finding provides additional basis for the potential of the plant (or its phytochemicals) as a COVID-19 therapeutic directly targeting drug targets such as SARS-CoV-2 Mpro and/or addressing respiratory-system-related symptoms. The study also highlights the utility of PM-VS, which can be efficiently implemented in the preliminary steps of drug discovery and development.

Computational Identification of Possible Allosteric Sites and Modulators of the SARS-CoV-2 Main Protease

In this study, we target the main protease (M^pro) of the SARS-CoV-2 virus as it is a crucial enzyme for viral replication. Herein, we report three plausible allosteric sites on M^pro that can expand structure-based drug discovery efforts for new M^pro inhibitors. To find these sites, we used mixed-solvent molecular dynamics (MixMD) simulations, an efficient computational protocol that finds binding hotspots through mapping the surface of unbound proteins with 5% cosolvents in water. We have used normal mode analysis to support our claim of allosteric control for these sites. Further, we have performed virtual screening against the sites with 361 hits from M^pro screenings available through the National Center for Advancing Translational Sciences (NCATS). We have identified the NCATS inhibitors that bind to the remote sites better than the active site of M^pro, and we propose these molecules may be allosteric regulators of the system. After identifying our sites, new X-ray crystal structures were released that show fragment molecules in the sites we found, supporting the notion that these sites are accurate and druggable.

Medicinal chemistry strategies towards the development of effective SARS-CoV-2 inhibitors

Novel therapies are urgently needed to improve global treatment of SARS-CoV-2 infection. Herein, we briefly provide a concise report on the medicinal chemistry strategies towards the development of effective SARS-CoV-2 inhibitors with representative examples in different strategies from the medicinal chemistry perspective.

A Computer-Aided Approach for the Discovery of D-Peptides as Inhibitors of SARS-CoV-2 Main Protease

The SARS-CoV-2 main protease, also known as 3-chymotrypsin-like protease (3CLpro), is a cysteine protease responsible for the cleavage of viral polyproteins pp1a and pp1ab, at least, at eleven conserved sites, which leads to the formation of mature nonstructural proteins essential for the replication of the virus. Due to its essential role, numerous studies have been conducted so far, which have confirmed 3CLpro as an attractive drug target to combat Covid-19 and have reported a vast number of inhibitors and their co-crystal structures. Despite all the ongoing efforts, D-peptides, which possess key advantages over L-peptides as therapeutic agents, have not been explored as potential drug candidates against 3CLpro. The current work fills this gap by reporting an in silico approach for the discovery of D-peptides capable of inhibiting 3CLpro that involves structure-based virtual screening (SBVS) of an in-house library of D-tripeptides and D-tetrapeptides into the protease active site and subsequent rescoring steps, including Molecular Mechanics Generalized-Born Surface Area (MM-GBSA) free energy calculations and molecular dynamics (MD) simulations. In vitro enzymatic assays conducted for the four top-scoring D-tetrapeptides at 20 μM showed that all of them caused 55–85% inhibition of 3CLpro activity, thus highlighting the suitability of the devised approach. Overall, our results present a promising computational strategy to identify D-peptides capable of inhibiting 3CLpro, with broader application in problems involving protein inhibition.

Quantum-Chemical Quasi-Docking for Molecular Dynamics Calculations

The quantum quasi-docking procedure is used to compare the docking accuracies of two quantum-chemical semiempirical methods, namely, PM6-D3H4X and PM7. Quantum quasi-docking is an approximation to quantum docking. In quantum docking, it is necessary to search directly for the global minimum of the energy of the protein-ligand complex calculated by the quantum-chemical method. In quantum quasi-docking, firstly, we look for a wide spectrum of low-energy minima, calculated using the MMFF94 force field, and secondly, we recalculate the energies of all these minima using the quantum-chemical method, and among these recalculated energies we determine the lowest energy and the corresponding ligand position. Both PM6-D3H4X and PM7 are novel methods that describe well-dispersion interactions, hydrogen and halogen bonds. The PM6-D3H4X and PM7 methods are used with the COSMO implicit solvent model as it is implemented in the MOPAC program. The comparison is made for 25 high quality protein-ligand complexes. Firstly, the docking positioning accuracies have been compared, and we demonstrated that PM7+COSMO provides better positioning accuracy than PM6-D3H4X. Secondly, we found that PM7+COSMO demonstrates a much higher correlation between the calculated and measured protein-ligand binding enthalpies than PM6-D3H4X. For future quantum docking PM7+COSMO is preferable, but the COSMO model must be improved.

Structure-based inhibitor design and repurposing clinical drugs to target SARS-CoV-2 proteases

SARS-CoV-2, the coronavirus responsible for the current COVID-19 pandemic, encodes two proteases, 3CLpro and PLpro, two of the main antiviral research targets. Here we provide an overview of the structures and functions of 3CLpro and PLpro and examine strategies of structure-based drug designing and drug repurposing against these proteases. Rational structure-based drug design enables the generation of potent and target-specific antivirals. Drug repurposing offers an attractive prospect with an accelerated turnaround. Thus far, several protease inhibitors have been identified, and some candidates are undergoing trials that may well prove to be effective antivirals against SARS-CoV-2.

Antiviral cyclic peptides targeting the main protease of SARS-CoV-2

Antivirals that specifically target SARS-CoV-2 are needed to control the COVID-19 pandemic. The main protease (M^pro) is essential for SARS-CoV-2 replication and is an attractive target for antiviral development. Here we report the use of the Random nonstandard Peptide Integrated Discovery (RaPID) mRNA display on a chemically cross-linked SARS-CoV-2 M^pro dimer, which yielded several high-affinity thioether-linked cyclic peptide inhibitors of the protease. Structural analysis of M^pro complexed with a selenoether analogue of the highest-affinity peptide revealed key binding interactions, including glutamine and leucine residues in sites S1 and S2, respectively, and a binding epitope straddling both protein chains in the physiological dimer. Several of these M^pro peptide inhibitors possessed antiviral activity against SARS-CoV-2 in vitro with EC₅₀ values in the low micromolar range. These cyclic peptides serve as a foundation for the development of much needed antivirals that specifically target SARS-CoV-2.

RaPID mRNA display was used for the discovery of antiviral cyclic peptides that potently and selectively inhibit SARS-CoV-2 M^pro. The most potent inhibitor exhibited a novel binding mode, interacting with residues across the homodimer interface.

Pharmacologic profiling reveals lapatinib as a novel antiviral against SARS-CoV-2 in vitro

The emergence of SARS-CoV-2 virus has resulted in a worldwide pandemic, but effective antiviral therapies are not widely available. To improve treatment options, we conducted a high-throughput screen to uncover compounds that block SARS-CoV-2 infection. A minimally pathogenic human betacoronavirus (OC43) was used to infect physiologically-relevant human pulmonary fibroblasts (MRC5) to facilitate rapid antiviral discovery in a preclinical model. Comprehensive profiling was conducted on more than 600 compounds, with each compound arrayed across 10 dose points. Our screening revealed several FDA-approved agents that can attenuate both OC43 and SARS-CoV-2 viral replication, including lapatinib, doramapimod, and 17-AAG. Importantly, lapatinib inhibited SARS-CoV-2 RNA replication by over 50,000-fold. Further, both lapatinib and doramapimod could be combined with remdesivir to improve antiviral activity in cells. These findings reveal novel therapeutic avenues that could limit SARS-CoV-2 infection.

A Review of the Current Landscape of SARS-CoV-2 Main Protease Inhibitors: Have We Hit the Bullseye Yet?

In this review, we collected 1765 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) M-pro inhibitors from the bibliography and other sources, such as the COVID Moonshot project and the ChEMBL database. This set of inhibitors includes only those compounds whose inhibitory capacity, mainly expressed as the half-maximal inhibitory concentration (IC50) value, against M-pro from SARS-CoV-2 has been determined. Several covalent warheads are used to treat covalent and non-covalent inhibitors separately. Chemical space, the variation of the IC50 inhibitory activity when measured by different methods or laboratories, and the influence of 1,4-dithiothreitol (DTT) are discussed. When available, we have collected the values of inhibition of viral replication measured with a cellular antiviral assay and expressed as half maximal effective concentration (EC50) values, and their possible relationship to inhibitory potency against M-pro is analyzed. Finally, the most potent covalent and non-covalent inhibitors that simultaneously inhibit the SARS-CoV-2 M-pro and the virus replication in vitro are discussed.

Pre-Steady-State Kinetics of the SARS-CoV-2 Main Protease as a Powerful Tool for Antiviral Drug Discovery

The design of effective target-specific drugs for COVID-19 treatment has become an intriguing challenge for modern science. The SARS-CoV-2 main protease, M^pro, responsible for the processing of SARS-CoV-2 polyproteins and production of individual components of viral replication machinery, is an attractive candidate target for drug discovery. Specific M^pro inhibitors have turned out to be promising anticoronaviral agents. Thus, an effective platform for quantitative screening of M^pro-targeting molecules is urgently needed. Here, we propose a pre-steady-state kinetic analysis of the interaction of M^pro with inhibitors as a basis for such a platform. We examined the kinetic mechanism of peptide substrate binding and cleavage by wild-type M^pro and by its catalytically inactive mutant C145A. The enzyme induces conformational changes of the peptide during the reaction. The inhibition of M^pro by boceprevir, telaprevir, GC-376, PF-00835231, or thimerosal was investigated. Detailed pre-steady-state kinetics of the interaction of the wild-type enzyme with the most potent inhibitor, PF-00835231, revealed a two-step binding mechanism, followed by covalent complex formation. The C145A M^pro mutant interacts with PF-00835231 approximately 100-fold less effectively. Nevertheless, the binding constant of PF-00835231 toward C145A M^pro is still good enough to inhibit the enzyme. Therefore, our results suggest that even noncovalent inhibitor binding due to a fine conformational fit into the active site is sufficient for efficient inhibition. A structure-based virtual screening and a subsequent detailed assessment of inhibition efficacy allowed us to select two compounds as promising noncovalent inhibitor leads of SARS-CoV-2 M^pro.

A Genetic Trap in Yeast for Inhibitors of SARS-CoV-2 Main Protease

The ongoing COVID-19 pandemic urges searches for antiviral agents that can block infection or ameliorate its symptoms. Using dissimilar search strategies for new antivirals will improve our overall chances of finding effective treatments. Here, we have established an experimental platform for screening of small molecule inhibitors of the SARS-CoV-2 main protease in Saccharomyces cerevisiae cells, genetically engineered to enhance cellular uptake of small molecules in the environment. The system consists of a fusion of the Escherichia coli toxin MazF and its antitoxin MazE, with insertion of a protease cleavage site in the linker peptide connecting the MazE and MazF moieties. Expression of the viral protease confers cleavage of the MazEF fusion, releasing the MazF toxin from its antitoxin, resulting in growth inhibition. In the presence of a small molecule inhibiting the protease, cleavage is blocked and the MazF toxin remains inhibited, promoting growth. The system thus allows positive selection for inhibitors. The engineered yeast strain is tagged with a fluorescent marker protein, allowing precise monitoring of its growth in the presence or absence of inhibitor. We detect an established main protease inhibitor by a robust growth increase, discernible down to 1 μM. The system is suitable for robotized large-scale screens. It allows in vivo evaluation of drug candidates and is rapidly adaptable for new variants of the protease with deviant site specificities.

IMPORTANCE The COVID-19 pandemic may continue for several years before vaccination campaigns can put an end to it globally. Thus, the need for discovery of new antiviral drug candidates will remain. We have engineered a system in yeast cells for the detection of small molecule inhibitors of one attractive drug target of SARS-CoV-2, its main protease, which is required for viral replication. The ability to detect inhibitors in live cells brings the advantage that only compounds capable of entering the cell and remain stable there will score in the system. Moreover, because of its design in yeast cells, the system is rapidly adaptable for tuning the detection level and eventual modification of the protease cleavage site in the case of future mutant variants of the SARS-CoV-2 main protease or even for other proteases.

The P3 O - Tert -Butyl-Threonine is Key to High Cellular and Antiviral Potency for Aldehyde-Based SARS-CoV-2 Main Protease Inhibitors

As an essential enzyme to SARS-CoV-2, main protease (M ^Pro ) is a viable target to develop antivirals for the treatment of COVID-19. By varying chemical compositions at both P2 and P3 sites and the N -terminal protection group, we synthesized a series of M ^Pro inhibitors that contain β -(S-2-oxopyrrolidin-3-yl)-alaninal at the P1 site. These inhibitors have a large variation of determined IC ₅₀ values that range from 4.8 to 650 nM. The determined IC ₅₀ values reveal that relatively small side chains at both P2 and P3 sites are favorable for achieving high in vitro M ^Pro inhibition potency, the P3 site is tolerable toward unnatural amino acids with two alkyl substituents on the α -carbon, and the inhibition potency is sensitive toward the N -terminal protection group. X-ray crystal structures of M ^Pro bound with 16 inhibitors were determined. All structures show similar binding patterns of inhibitors at the M ^Pro active site. A covalent interaction between the active site cysteine and a bound inhibitor was observed in all structures. In M ^Pro , large structural variations were observed on residues N142 and Q189. All inhibitors were also characterized on their inhibition of M ^Pro in 293T cells, which revealed their in cellulo potency that is drastically different from their in vitro enzyme inhibition potency. Inhibitors that showed high in cellulo potency all contain O - tert -butyl-threonine at the P3 site. Based on the current and a previous study, we conclude that O - tert -butyl-threonine at the P3 site is a key component to achieve high cellular and antiviral potency for peptidyl aldehyde inhibitors of M ^Pro . This finding will be critical to the development of novel antivirals to address the current global emergency of concerning the COVID-19 pandemic.

The N -Terminal Carbamate is Key to High Cellular and Antiviral Potency for Boceprevir-Based SARS-CoV-2 Main Protease Inhibitors

Boceprevir is an HCV NSP3 inhibitor that has been explored as a repurposed drug for COVID-19. It inhibits the SARS-CoV-2 main protease (M ^Pro ) and contains an α-ketoamide warhead, a P1 β-cyclobutylalanyl moiety, a P2 dimethylcyclopropylproline, a P3 tert -butyl-glycine, and a P4 N -terminal tert -butylcarbamide. By introducing modifications at all four positions, we synthesized 20 boceprevir-based M ^Pro inhibitors including PF-07321332 and characterized their M ^Pro inhibition potency in test tubes ( in vitro ) and human host cells ( in cellulo ). Crystal structures of M ^Pro bound with 10 inhibitors and antiviral potency of 4 inhibitors were characterized as well. Replacing the P1 site with a β-(S-2-oxopyrrolidin-3-yl)-alanyl (opal) residue and the warhead with an aldehyde leads to high in vitro potency. The original moieties at P2, P3 and the P4 N -terminal cap positions in boceprevir are better than other tested chemical moieties for high in vitro potency. In crystal structures, all inhibitors form a covalent adduct with the M ^Pro active site cysteine. The P1 opal residue, P2 dimethylcyclopropylproline and P4 N -terminal tert -butylcarbamide make strong hydrophobic interactions with M ^Pro , explaining high in vitro potency of inhibitors that contain these moieties. A unique observation was made with an inhibitor that contains an P4 N -terminal isovaleramide. In its M ^Pro complex structure, the P4 N -terminal isovaleramide is tucked deep in a small pocket of M ^Pro that originally recognizes a P4 alanine side chain in a substrate. Although all inhibitors show high in vitro potency, they have drastically different in cellulo potency in inhibiting ectopically expressed M ^Pro in human 293T cells. All inhibitors including PF-07321332 with a P4 N -terminal carbamide or amide have low in cellulo potency. This trend is reversed when the P4 N -terminal cap is changed to a carbamate. The installation of a P3 O-tert -butyl-threonine improves in cellulo potency. Three molecules that contain a P4 N -terminal carbamate were advanced to antiviral tests on three SARS-CoV-2 variants. They all have high potency with EC ₅₀ values around 1 μM. A control compound with a nitrile warhead and a P4 N -terminal amide has undetectable antiviral potency. Based on all observations, we conclude that a P4 N -terminal carbamate in a boceprevir derivative is key for high antiviral potency against SARS-CoV-2.

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

Mass spectrometry 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 SARS-CoV-2 3CL protease inhibitors. Primary screening of 32,033 chemical samples was completed in 12 hours. Among the hits showing selective, dose-dependent 3CL-inhibitory activity, 8 compounds showed antiviral activity in cell-based assay.

SARS-CoV-2 Antiviral Therapy

The development of effective antiviral therapy for COVID-19 is critical for those awaiting vaccination, as well as for those who do not respond robustly to vaccination. This review summarizes 1 year of progress in the race to develop antiviral therapies for COVID-19, including research spanning preclinical and clinical drug development efforts, with an emphasis on antiviral compounds that are in clinical development or that are high priorities for clinical development. The review is divided into sections on compounds that inhibit SARS-CoV-2 enzymes, including its polymerase and proteases; compounds that inhibit virus entry, including monoclonal antibodies; interferons; and repurposed drugs that inhibit host processes required for SARS-CoV-2 replication. The review concludes with a summary of the lessons to be learned from SARS-CoV-2 drug development efforts and the challenges to continued progress.

Targeting SARS-CoV-2 main protease by teicoplanin: A mechanistic insight by docking, MM/GBSA and molecular dynamics simulation

First emerged in late December 2019, the outbreak of novel severe acute respiratory syndrome corona virus-2 (SARS-CoV-2) pandemic has instigated public-health emergency around the globe. Till date there is no specific therapeutic agent for this disease and hence, the world is craving to identify potential antiviral agents against SARS-CoV-2. The main protease (M^Pro) is considered as an attractive drug target for rational drug design against SARS-CoV-2 as it is known to play a crucial role in the viral replication and transcription. Teicoplanin is a glycopeptide class of antibiotic which is regularly used for treating Gram-positive bacterial infections, has shown potential therapeutic efficacy against SARS-CoV-2 in vitro. Therefore, in this study, a mechanistic insight of intermolecular interactions between teicoplanin and SARS-CoV-2 M^Pro has been scrutinized by molecular docking. Both monomeric and dimeric forms of M^Pro was used in docking involving blind as well as defined binding site based on the known inhibitor. Binding energies of teicoplanin-M^Pro complexes were estimated by Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) computations from docking and simulated trajectories. The dynamic and thermodynamics constraints of docked drug in complex with target proteins under specific physiological conditions was ascertained by all-atom molecular dynamics simulation of 100 ns trajectory. Root mean square deviation and fluctuation of carbon α chain justified the stability of the bound complex in biological environments. The outcomes of current study are supposed to be fruitful in rational design of antiviral drugs against SARS-CoV-2.

Automated discovery of noncovalent inhibitors of SARS-CoV-2 main protease by consensus Deep Docking of 40 billion small molecules

Recent explosive growth of 'make-on-demand' chemical libraries brought unprecedented opportunities but also significant challenges to the field of computer-aided drug discovery. To address this expansion of the accessible chemical universe, molecular docking needs to accurately rank billions of chemical structures, calling for the development of automated hit-selecting protocols to minimize human intervention and error. Herein, we report the development of an artificial intelligence-driven virtual screening pipeline that utilizes Deep Docking with Autodock GPU, Glide SP, FRED, ICM and QuickVina2 programs to screen 40 billion molecules against SARS-CoV-2 main protease (Mpro). This campaign returned a significant number of experimentally confirmed inhibitors of Mpro enzyme, and also enabled to benchmark the performance of twenty-eight hit-selecting strategies of various degrees of stringency and automation. These findings provide new starting scaffolds for hit-to-lead optimization campaigns against Mpro and encourage the development of fully automated end-to-end drug discovery protocols integrating machine learning and human expertise.

Structural, Electronic, and Electrostatic Determinants for Inhibitor Binding to Subsites S1 and S2 in SARS-CoV-2 Main Protease

Creating small-molecule antivirals specific for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) proteins is crucial to battle coronavirus disease 2019 (COVID-19). SARS-CoV-2 main protease (Mpro) is an established drug target for the design of protease inhibitors. We performed a structure-activity relationship (SAR) study of noncovalent compounds that bind in the enzyme's substrate-binding subsites S1 and S2, revealing structural, electronic, and electrostatic determinants of these sites. The study was guided by the X-ray/neutron structure of Mpro complexed with Mcule-5948770040 (compound 1), in which protonation states were directly visualized. Virtual reality-assisted structure analysis and small-molecule building were employed to generate analogues of 1. In vitro enzyme inhibition assays and room-temperature X-ray structures demonstrated the effect of chemical modifications on Mpro inhibition, showing that (1) maintaining correct geometry of an inhibitor's P1 group is essential to preserve the hydrogen bond with the protonated His163; (2) a positively charged linker is preferred; and (3) subsite S2 prefers nonbulky modestly electronegative groups.

Colloidal Aggregators in Biochemical SARS-CoV-2 Repurposing Screens

To fight COVID-19, much effort has been directed toward in vitro drug repurposing. Here, we investigate the impact of colloidal aggregation, a common screening artifact, in these repurposing campaigns. We tested 56 drugs reported as active in biochemical assays for aggregation by dynamic light scattering and by detergent-based enzyme counter screening; 19 formed colloids at concentrations similar to their literature IC₅₀'s, and another 14 were problematic. From a common repurposing library, we further selected another 15 drugs that had physical properties resembling known aggregators, finding that six aggregated at micromolar concentrations. This study suggests not only that many of the drugs repurposed for SARS-CoV-2 in biochemical assays are artifacts but that, more generally, at screening-relevant concentrations, even drugs can act artifactually via colloidal aggregation. Rapid detection of these artifacts will allow the community to focus on those molecules that genuinely have potential for treating COVID-19.

Target identification for repurposed drugs active against SARS-CoV-2 via high-throughput inverse docking

Screening already approved drugs for activity against a novel pathogen can be an important part of global rapid-response strategies in pandemics. Such high-throughput repurposing screens have already identified several existing drugs with potential to combat SARS-CoV-2. However, moving these hits forward for possible development into drugs specifically against this pathogen requires unambiguous identification of their corresponding targets, something the high-throughput screens are not typically designed to reveal. We present here a new computational inverse-docking protocol that uses all-atom protein structures and a combination of docking methods to rank-order targets for each of several existing drugs for which a plurality of recent high-throughput screens detected anti-SARS-CoV-2 activity. We demonstrate validation of this method with known drug-target pairs, including both non-antiviral and antiviral compounds. We subjected 152 distinct drugs potentially suitable for repurposing to the inverse docking procedure. The most common preferential targets were the human enzymes TMPRSS2 and PIKfyve, followed by the viral enzymes Helicase and PLpro. All compounds that selected TMPRSS2 are known serine protease inhibitors, and those that selected PIKfyve are known tyrosine kinase inhibitors. Detailed structural analysis of the docking poses revealed important insights into why these selections arose, and could potentially lead to more rational design of new drugs against these targets.

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

The main protease (Mpro) 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 Mpro 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 Mpro 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 Mpro-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 Mpro competitively. Our combined results provide new insights and highlight opportunities for the development of Mpro inhibitors as anti-COVID-19 drugs.

Haste makes waste: A critical review of docking-based virtual screening in drug repurposing for SARS-CoV-2 main protease (M-pro) inhibition

This review makes a critical evaluation of 61 peer‐reviewed manuscripts that use a docking step in a virtual screening (VS) protocol to predict SARS‐CoV‐2 M‐pro (M‐pro) inhibitors in approved or investigational drugs. Various manuscripts predict different compounds, even when they use a similar initial dataset and methodology, and most of them do not validate their methodology or results. In addition, a set of known 150 SARS‐CoV‐2 M‐pro inhibitors extracted from the literature and a second set of 81 M‐pro inhibitors and 113 inactive compounds obtained from the COVID Moonshot project were used to evaluate the reliability of using docking scores as feasible predictors of the potency of a SARS‐CoV‐2 M‐pro inhibitor. Using two SARS‐CoV‐2 M‐pro structures and five protein‐ligand docking programs, we proved that the correlation between the pIC₅₀ and docking scores is not good. Neither was any correlation found between the pIC₅₀ and the ∆G calculated with an MM‐GBSA method. When a group of experimentally known inactive compounds was added, neither the docking scores or the ∆G were able to distinguish between compounds with or without M‐pro experimental inhibitory activity. Performances improved when covalent and noncovalent inhibitors were treated separately, but were not good enough to fully support using a docking score as a cutoff value for selecting new putative M‐pro inhibitors or predicting the relative bioactivity of a compound by comparison with a reference compound. The two sets of known SARS‐CoV‐2 M‐pro inhibitors presented here could be used for validating future VS protocols which aim to predict M‐pro inhibitors.

Development of a Novel Pharmacophore Model Guided by the Ensemble of Waters and Small Molecule Fragments Bound to SARS-CoV-2 Main Protease

Recent fragment‐based drug design efforts have generated huge amounts of information on water and small molecule fragment binding sites on SARS‐CoV‐2 M^pro and preference of the sites for various types of chemical moieties. However, this information has not been effectively utilized to develop automated tools for in silico drug discovery which are routinely used for screening large compound libraries. Utilization of this information in the development of pharmacophore models can help in bridging this gap. In this study, information on water and small molecule fragments bound to M^pro has been utilized to develop a novel Water Pharmacophore (Waterphore) model. The Waterphore model can also implicitly represent the conformational flexibilities of binding pockets in terms of pharmacophore features. The Waterphore model derived from 173 apo‐ or small molecule fragment‐bound structures of M^pro has been validated by using a dataset of 68 known bioactive inhibitors and 78 crystal structure bound inhibitors of SARS‐CoV‐2 M^pro. It is encouraging to note that, even though no inhibitor data has been used in developing the Waterphore model, it could successfully identify the known inhibitors from a library of decoys with a ROC‐AUC of 0.81 and active hit rate (AHR) of 70 %. The Waterphore model is also general enough for potential applications for other drug targets.

Pan-3C Protease Inhibitor Rupintrivir Binds SARS-CoV-2 Main Protease in a Unique Binding Mode

Rupintrivir targets the 3C cysteine proteases of the picornaviridae family, which includes rhinoviruses and enteroviruses that cause a range of human diseases. Despite being a pan-3C protease inhibitor, rupintrivir activity is extremely weak against the homologous 3C-like protease of SARS-CoV-2. In this study, the crystal structures of rupintrivir were determined bound to enterovirus 68 (EV68) 3C protease and the 3C-like main protease (Mpro) from SARS-CoV-2. While the EV68 3C protease-rupintrivir structure was similar to previously determined complexes with other picornavirus 3C proteases, rupintrivir bound in a unique conformation to the active site of SARS-CoV-2 Mpro splitting the catalytic cysteine and histidine residues. This bifurcation of the catalytic dyad may provide a novel approach for inhibiting cysteine proteases.

Asunaprevir, a Potent Hepatitis C Virus Protease Inhibitor, Blocks SARS-CoV-2 Propagation

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has become a global health concern. Various SARS-CoV-2 vaccines have been developed and are being used for vaccination worldwide. However, no therapeutic agents against coronavirus disease 2019 (COVID-19) have been developed so far; therefore, new therapeutic agents are urgently needed. In the present study, we evaluated several hepatitis C virus direct-acting antivirals as potential candidates for drug repurposing against COVID-19. Theses include asunaprevir (a protease inhibitor), daclatasvir (an NS5A inhibitor), and sofosbuvir (an RNA polymerase inhibitor). We found that asunaprevir, but not sofosbuvir and daclatasvir, markedly inhibited SARS-CoV-2-induced cytopathic effects in Vero E6 cells. Both RNA and protein levels of SARS-CoV-2 were significantly decreased by treatment with asunaprevir. Moreover, asunaprevir profoundly decreased virion release from SARS-CoV-2-infected cells. A pseudoparticle entry assay revealed that asunaprevir blocked SARS-CoV-2 infection at the binding step of the viral life cycle. Furthermore, asunaprevir inhibited SARS-CoV-2 propagation in human lung Calu-3 cells. Collectively, we found that asunaprevir displays broad-spectrum antiviral activity and therefore might be worth developing as a new drug repurposing candidate for COVID-19.