Improving health and quality of life in our society is a key focus of drug development. Methods for drug discovery are being optimized in multiple ways to reduce costs and timelines. Crystallographic fragment screening (CFS) is increasingly being employed as an early screening method in drug discovery projects. Here, we demonstrate that selecting the optimal protein crystal form can significantly impact hit rates. Two CFS campaigns are carried out against the two crystal forms of the SARS-CoV-2 main protease, using the same fragment library and an almost identical experimental setup. Although both crystal forms exhibit similar diffraction properties, the observed hit rates in the two campaigns differ significantly. A hit rate of 3% is determined for the monoclinic crystals, while a hit rate of 16% is observed for the orthorhombic crystals. These findings are consistent with the more open molecular packing in the orthorhombic crystals, where the solvent channels leading to the active sites are approximately twice the size of those in the monoclinic crystal form. Our results highlight the critical importance of the crystal form in a crystallographic screening, identifying it as one of the most important parameters to optimize when preparing a CFS campaign.
Regardless of the massive global efforts to combat the virus causing SARS-CoV-2 syndrome, the infection remains a substantial health challenge worldwide in the years following 2019. To this direction, targeting the main viral protease Mpro has been proposed as a tractable and particularly promising approach toward developing effective and safe COVID-19 antivirals. By applying an integrated workflow combining a previously developed in silico consensus ranking protocol with two orthogonal in vitro methods, the NCI/DTP repository is screened and the discovery of an original naphthol scaffold with Mpro inhibitory properties is reported. The hit is characterized in terms of structure and binding thermodynamics by combining X-ray crystallography and isothermal titration calorimetry where a binding affinity constant of 1.55 μM is determined. The compound is further evaluated against virus-infected cells, where an EC50 value of 7.23 μM and comparable toxicity with nirmatrelvir are measured. Chemical synthesis is additionally employed to facilitate optimal exploration of the structure-activity relationship landscape regarding the new hit. By integrating computational, biophysical, and enzymatic methods, the suggested approach allows the combination of a structural hypothesis with functional evidence and shows its capacity toward identifying and rationally optimizing structurally original noncovalent Mpro inhibitors.
Abstract Viral proteases represent validated targets for direct-acting antivirals and the treatment of associated infections. In co-crystal structures of M pro of SARS-CoV-2 with peptidomimetic inhibitors, we noticed a spatial proximity of sidechains filling the S1’ and S2 pockets, as well as those filling S3 and S1 pockets. To enhance molecular rigidity, the proximal residues were conformationally fixed by macrocyclization. We report the synthesis of two macrocyclic series, i.e. exocyclic nitriles with linked P3 and P1 residues and endocyclic α-ketoamides with linked P1’ and P2 residues, and characterize their binding modes and bioactivities. The 17-membered macrocyclic α-ketoamide 20 f inhibited M pro (IC₅₀ = 370 nM) and exerted anti-SARS-CoV-2 effects (EC₅₀ = 1.9 μM). Leveraging structural similarities between M pro and the 3C pro of enterovirus D68, we describe with two co-crystal structures how α-ketoamide macrocycles bound to and inhibited the enteroviral protease. Notably, 20 f exhibited very potent antiviral activities with EC₅₀‘s of 33, 133, and 146 nM against EV-D68, EV-A71, and CVB3, respectively. The study demonstrates how broad-spectrum activity can be achieved with direct-acting antivirals.
The SARS-CoV-2 main protease (Mpro), an enzyme essential for viral replication and lacking a human homologue, has emerged as a highly attractive target for the development of novel antiviral agents. Although several Mpro inhibitors have been developed — some receiving regulatory approval — their use is sometimes limited by drug-drug interactions. In this study, we designed and synthesized peptidomimetic SARS-CoV-2 Mpro inhibitors incorporating a novel thiazolyl 4-carboxylate ketone warhead, previously employed by our group in the development of cytosolic phospholipase A2 inhibitors. The synthesized compounds were evaluated for their in vitro inhibitory potency against SARS-CoV-2 Mpro, and a highly potent Mpro inhibitor (GK730) was identified (IC50 5.75 nM). The melting temperature of the Mpro-GK730 complex revealed high stability, consistent with the high inhibitory potency. The X-ray crystal structures of inhibitors GK729 and GK730 bound to Mpro were determined, providing insights into the binding interactions and mechanism of action. Studies on the host cell proteases cathepsin B and L showed that GK730 did not inhibit cathepsin B, while exhibited weak inhibition of cathepsin L. Furthermore, GK730 demonstrated an EC50 value of 5.70 μM against a wild-type SARS-CoV-2 strain in Vero E6 cells and minimal cytotoxicity (CC50 value greater than 100 μM).
In more and more drug discovery projects, crystallographic fragment screening (CFS) is employed as an early screening method. Here, we demonstrate that choosing the right crystal form has a profound influence on the hit rates and hence success and speed of downstream lead generation. Two CFS campaigns with the same fragment library and an almost identical experimental setup were carried out against the two crystal forms of the SARS-CoV-2 main protease.While both crystal forms exhibit similar diffraction properties, the observed hit rates in the two campaigns were vastly different. For the monoclinic crystals a hit rate of 3% was determined, while a hit rate of 16% was observed for the orthorhombic crystals. These findings align with the more open molecular packing in the orthorhombic crystals where the solvent channels leading to the active sites are about twice larger than in the monoclinic crystal form. Our results highlight the critical importance of the crystal system in a crystallographic fragment-screening campaign and identify this parameter as one of the most important ones to be optimized during preparation of a campaign. ### Competing Interest Statement The authors have declared no competing interest. BMBF, 05K20CB1 iNext Discovery, 871037 Horizon 2020 program, FE2166/1-1 Helmholtz Gemeinschaft
The coronavirus membrane protein (M) is the main organizer of coronavirus assembly1-3. Here, we report on an M-targeting molecule, CIM-834, that blocks the assembly of SARS-CoV-2. CIM-834 was obtained through high-throughput phenotypic antiviral screening followed by medicinal-chemistry efforts and target elucidation. CIM-834 inhibits the replication of SARS-CoV-2 (including a broad panel of variants) and SARS-CoV. In SCID mice and Syrian hamsters intranasally infected with SARS-CoV-2, oral treatment reduced lung viral titres to nearly undetectable levels, even (as shown in mice) when treatment was delayed until 24 h before the end point. Treatment of infected hamsters prevented transmission to untreated sentinels. Transmission electron microscopy studies show that virion assembly is completely absent in cells treated with CIM-834. Single-particle cryo-electron microscopy reveals that CIM-834 binds and stabilizes the M protein in its short form, thereby preventing the conformational switch to the long form, which is required for successful particle assembly. In conclusion, we have discovered a new druggable target in the replication cycle of coronaviruses and a small molecule that potently inhibits it.
The main protease Mpro is a clinically validated target to treat infections by the coronavirus SARS-CoV-2. Among the first reported Mpro inhibitors was the peptidomimetic α-ketoamide 13b, whose cocrystal structure with Mpro paved the way for multiple lead-finding studies. We established structure-activity relationships for the 13b series by modifying residues at the P1', P3, and P4 sites. Guided by cocrystal structures, we reduced the P1' substituent size to better fill the pocket and added a fluorine substituent to the pyridone ring, enabling a new hydrogen bond with Gln189 in P3. Among 22 novel analogues, 6d and 12d inhibited Mpro with IC50s of 110 nM and 40 nM, improving the potency of 13b by up to 9.5-fold. Compound 6d had pronounced antiviral activity with an EC50 of 1.6 μM and was stable in plasma and microsomes. The study illustrates the potential of structure-based design to systematically improve peptidomimetic α-ketoamides.
The SARS-CoV-2 main protease (Mpro) is a well-established target for antiviral drug development One such inhibitor, nirmatrelvir, in combination with ritonavir as a booster, has already been introduced into the market, under the name Paxlovid. However, being an RNA virus, SARS-CoV-2 is prone to the emergence of resistance mutations. A number of such mutations have been characterized, although they have not yet been shown to play a significant role in clinical settings; these include S144A, E166V, H172Y, and Q189K. We recombinantly produced these mutants and studied the corresponding proteins using X-ray crystallography, enzymology, and biophysical approaches. We discuss the potential of each mutant to lead to a widespread nirmatrelvir resistance scenario. We also demonstrate that one of our own inhibitors (13b-K), while showing some degree of cross-resistance with nirmatrelvir, exhibits much higher inhibitory activity against the Mpro carrying the E166V mutation.
The first Federation of European Biochemical Societies Advanced Course on macromolecular crystallization was launched in the Czech Republic in October 2004. Over the past two decades, the course has developed into a distinguished event, attracting students, early career postdoctoral researchers and lecturers. The course topics include protein purification, characterization and crystallization, covering the latest advances in the field of structural biology. The many hands-on practical exercises enable a close interaction between students and teachers and offer the opportunity for students to crystallize their own proteins. The course has a broad and lasting impact on the scientific community as participants return to their home laboratories and act as nuclei by communicating and implementing their newly acquired knowledge and skills.
Besides vaccines, the development of antiviral drugs targeting SARS-CoV-2 is critical for preventing future COVID outbreaks. The SARS-CoV-2 main protease (M pro ), a cysteine protease with essential functions in viral replication, has been validated as an effective drug target. Here, we show that M pro is subject to redox regulation in vitro and reversibly switches between the enzymatically active dimer and the functionally dormant monomer through redox modifications of cysteine residues. These include a disulfide-dithiol switch between the catalytic cysteine C145 and cysteine C117, and generation of an allosteric cysteine-lysine-cysteine SONOS bridge that is required for structural stability under oxidative stress conditions, such as those exerted by the innate immune system. We identify homo- and heterobifunctional reagents that mimic the redox switching and inhibit M pro activity. The discovered redox switches are conserved in main proteases from other coronaviruses, e.g. MERS-CoV and SARS-CoV, indicating their potential as common druggable sites.
During the continuing evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the Omicron variant of concern emerged in the second half of 2021 and has been dominant since November of that year. Along with its sublineages, it has maintained a prominent role ever since. The Nsp5 main protease (Mpro) of the Omicron virus is characterized by a single dominant mutation, P132H. Here we determined the X-ray crystal structures of the P132H mutant (or O-Mpro) as a free enzyme and in complex with the Mpro inhibitor, the alpha-ketoamide 13b-K, and we conducted enzymological, biophysical as well as theoretical studies to characterize the O-Mpro. We found that O-Mpro has a similar overall structure and binding with 13b-K; however, it displays lower enzymatic activity and lower thermal stability compared to the WT-Mpro (with “WT” referring to the prototype strain). Intriguingly, the imidazole ring of His132 and the carboxylate plane of Glu240 are in a stacked configuration in the X-ray structures determined here. Empirical folding free energy calculations suggest that the O-Mpro dimer is destabilized relative to the WT-Mpro due to less favorable van der Waals interactions and backbone conformations in the individual protomers. All-atom continuous constant-pH molecular dynamics (MD) simulations reveal that His132 and Glu240 display coupled titration. At pH 7, His132 is predominantly neutral and in a stacked configuration with respect to Glu240 which is charged. In order to examine whether the Omicron mutation eases the emergence of further Mpro mutations, we also analyzed the P132H+T169S double mutant, which is characteristic of the BA.1.1.2 lineage. However, we found little evidence of a correlation between the two mutation sites.
The main protease from coronaviruses and the 3C protease from enteroviruses play a crucial role in processing viral polyproteins, making them attractive targets for the development of antiviral agents. In this study, we employed a combinatorial chemistry approach-HyCoSuL-to compare the substrate specificity profiles of the main and 3C proteases from alphacoronaviruses, betacoronaviruses, and enteroviruses. The obtained data demonstrate that coronavirus M(pro)s exhibit overlapping substrate specificity in all binding pockets, whereas the 3C(pro) from enterovirus displays slightly different preferences toward natural and unnatural amino acids at the P4-P2 positions. However, chemical tools such as substrates, inhibitors, and activity-based probes developed for SARS-CoV-2 M-pro can be successfully applied to investigate the activity of the M-pro from other coronaviruses as well as the 3C(pro) from enteroviruses. Our study provides a structural framework for the development of broad-spectrum antiviral compounds.
During the continuing evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the Omicron variant of concern emerged in the second half of 2021 and has been dominant since November of that year. Along with its sublineages, it has maintained a prominent role ever since. The Nsp5 main protease (Mpro) of the Omicron virus is characterized by a single dominant mutation, P132H. Here we determined the X-ray crystal structures of the P132H mutant (or O-Mpro) as a free enzyme and in complex with the Mpro inhibitor, the alpha-ketoamide 13b-K, and we conducted enzymological, biophysical, as well as theoretical studies to characterize the O-Mpro. We found that O-Mpro has a similar overall structure and binding with 13b-K; however, it displays lower enzymatic activity and lower thermal stability compared to the WT-Mpro (with "WT" referring to the prototype strain). Intriguingly, the imidazole ring of His132 and the carboxylate plane of Glu240 are in a stacked configuration in the X-ray structures determined here. Empirical folding free energy calculations suggest that the O-Mpro dimer is destabilized relative to the WT-Mpro due to less favorable van der Waals interactions and backbone conformations in the individual protomers. All-atom continuous constant-pH molecular dynamics (MD) simulations reveal that His132 and Glu240 display coupled titration. At pH 7, His132 is predominantly neutral and in a stacked configuration with respect to Glu240 which is charged. In order to examine whether the Omicron mutation eases the emergence of further Mpro mutations, we also analyzed the P132H+T169S double mutant, which is characteristic of the BA.1.1.2 lineage. However, we found little evidence of a correlation between the two mutation sites.
Proteases encoded by SARS-CoV-2 constitute a promising target for new therapies against COVID-19. SARS-CoV-2 main protease (Mpro, 3CLpro) and papain-like protease (PLpro) are responsible for viral polyprotein cleavage-a process crucial for viral survival and replication. Recently it was shown that 2-phenylbenzisoselenazol-3(2H)-one (ebselen), an organoselenium anti-inflammatory small-molecule drug, is a potent, covalent inhibitor of both the proteases and its potency was evaluated in enzymatic and antiviral assays. In this study, we screened a collection of 34 ebselen and ebselen diselenide derivatives for SARS-CoV-2 PLpro and Mpro inhibitors. Our studies revealed that ebselen derivatives are potent inhibitors of both the proteases. We identified three PLpro and four Mpro inhibitors superior to ebselen. Independently, ebselen was shown to inhibit the N7-methyltransferase activity of SARS-CoV-2 nsp14 protein involved in viral RNA cap modification. Hence, selected compounds were also evaluated as nsp14 inhibitors. In the second part of our work, we employed 11 ebselen analogues-bis(2-carbamoylaryl)phenyl diselenides-in biological assays to evaluate their anti-SARS-CoV-2 activity in Vero E6 cells. We present their antiviral and cytoprotective activity and also low cytotoxicity. Our work shows that ebselen, its derivatives, and diselenide analogues constitute a promising platform for development of new antivirals targeting the SARS-CoV-2 virus.
The coronavirus disease 2019 (COVID-19) pandemic has stimulated tremendous efforts to develop therapeutic strategies that target severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and/or human proteins to control viral infection, encompassing hundreds of potential drugs and thousands of patients in clinical trials. So far, a few small-molecule antiviral drugs (nirmatrelvir-ritonavir, remdesivir and molnupiravir) and 11 monoclonal antibodies have been marketed for the treatment of COVID-19, mostly requiring administration within 10 days of symptom onset. In addition, hospitalized patients with severe or critical COVID-19 may benefit from treatment with previously approved immunomodulatory drugs, including glucocorticoids such as dexamethasone, cytokine antagonists such as tocilizumab and Janus kinase inhibitors such as baricitinib. Here, we summarize progress with COVID-19 drug discovery, based on accumulated findings since the pandemic began and a comprehensive list of clinical and preclinical inhibitors with anti-coronavirus activities. We also discuss the lessons learned from COVID-19 and other infectious diseases with regard to drug repurposing strategies, pan-coronavirus drug targets, in vitro assays and animal models, and platform trial design for the development of therapeutics to tackle COVID-19, long COVID and pathogenic coronaviruses in future outbreaks.
Nirmatrelvir is an orally available inhibitor of SARS-CoV-2 main protease (Mpro) and the main ingredient of Paxlovid, a drug approved by the U.S. Food and Drug Administration for high-risk COVID-19 patients. Recently, a rare natural mutation, H172Y, was found to significantly reduce nirmatrelvir's inhibitory activity. As the COVID-19 cases skyrocket in China and the selective pressure of antiviral therapy builds in the US, there is an urgent need to characterize and understand how the H172Y mutation confers drug resistance. Here, we investigated the H172Y Mpro's conformational dynamics, folding stability, catalytic efficiency, and inhibitory activity using all-atom constant pH and fixed-charge molecular dynamics simulations, alchemical and empirical free energy calculations, artificial neural networks, and biochemical experiments. Our data suggest that the mutation significantly weakens the S1 pocket interactions with the N-terminus and perturbs the conformation of the oxyanion loop, leading to a decrease in the thermal stability and catalytic efficiency. Importantly, the perturbed S1 pocket dynamics weaken the nirmatrelvir binding in the P1 position, which explains the decreased inhibitory activity of nirmatrelvir. Our work demonstrates the predictive power of the combined simulation and artificial intelligence approaches, and together with biochemical experiments, they can be used to actively surveil continually emerging mutations of SARS-CoV-2 Mpro and assist the optimization of antiviral drugs. The presented approach, in general, can be applied to characterize mutation effects on any protein drug targets.
The development of novel agents to combat COVID-19 is of high importance. SARS-CoV-2 main protease (Mpro) is a highly attractive target for the development of novel antivirals and a variety of inhibitors have already been developed. Accumulating evidence on the pathobiology of COVID-19 has shown that lipids and lipid metabolizing enzymes are critically involved in the severity of the infection. The purpose of the present study was to identify an inhibitor able to simultaneously inhibit both SARS-CoV-2 Mpro and phospholipase A2 (PLA2), an enzyme which plays a significant role in inflammatory diseases. Evaluating several PLA2 inhibitors, we demonstrate that the previously known potent inhibitor of Group IIA secretory PLA2, GK241, may also weakly inhibit SARS-CoV-2 Mpro. Molecular mechanics docking and molecular dynamics calculations shed light on the interactions between GK241 and SARS-CoV-2 Mpro. 2-Oxoamide GK241 may represent a lead molecular structure for the development of dual PLA2 and SARS-CoV-2 Mpro inhibitors.
The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19, has led to a pandemic, that continues to be a huge public health burden. Despite the availability of vaccines, there is still a need for small-molecule antiviral drugs. In an effort to identify novel and drug-like hit matter that can be used for subsequent hit-to-lead optimization campaigns, we conducted a high-throughput screening of a 160 K compound library against SARS-CoV-2, yielding a 1-heteroaryl-2-alkoxyphenyl analog as a promising hit. Antiviral profiling revealed this compound was active against various beta-coronaviruses and preliminary mode-of-action experiments demonstrated that it interfered with viral entry. A systematic structure–activity relationship (SAR) study demonstrated that a 3- or 4-pyridyl moiety on the oxadiazole moiety is optimal, whereas the oxadiazole can be replaced by various other heteroaromatic cycles. In addition, the alkoxy group tolerates some structural diversity.
Nirmatrelvir is an orally available inhibitor of SARS-CoV-2 main protease (Mpro) and the main ingredient of PAXLOVID, a drug approved by FDA for high-risk COVID-19 patients. Although the prevalent Mpro mutants in the SARS-CoV-2 Variants of Concern (e.g., Omicron) are still susceptible to nirmatrelvir, a rare natural mutation, H172Y, was found to significantly reduce nirmatrelvir's inhibitory activity. As the selective pressure of antiviral therapy may favor resistance mutations, there is an urgent need to understand the effect of the H172Y mutation on Mpro's structure, function, and drug resistance. Here we report the molecular dynamics (MD) simulations as well as the measurements of stability, enzyme kinetics of H172Y Mpro, and IC50 value of nirmatrelvir. Simulations showed that mutation disrupts the interactions between the S1 pocket and N terminus of the opposite protomer. Intriguingly, a native hydrogen bond (H-bond) between Phe140 and the N terminus is replaced by a transient H-bond between Phe140 and Tyr172. In the ligand-free simulations, strengthening of this nonnative H-bond is correlated with disruption of the conserved aromatic stacking between Phe140 and His163, leading to a partial collapse of the oxyanion loop. In the nirmatrelvir-bound simulations, the nonnative H-bond is correlated with the loss of an important H-bond between Glu166 and nirmatrelvir's lactam nitrogen at P1 position. These results are consistent with the newly reported X-ray structures of H172Y Mpro and suggest a mechanism by which the H172Y substitution perturbs the S1 pocket, leading to the decreased structural stability and binding affinity, which in turn explains the drastic reduction in catalytic activity and antiviral susceptibility.
SARS-CoV-2 is the causative agent behind the COVID-19 pandemic. The main protease (Mpro, 3CLpro) of SARS-CoV-2 is a key enzyme that processes polyproteins translated from the viral RNA. Mpro is therefore an attractive target for the design of inhibitors that block viral replication. We report the diastereomeric resolution of the previously designed SARS-CoV-2 Mpro α-ketoamide inhibitor 13b. The pure (S,S,S)-diastereomer, 13b-K, displays an IC50 of 120 nM against the Mpro and EC50 values of 0.8-3.4 μM for antiviral activity in different cell types. Crystal structures have been elucidated for the Mpro complexes with each of the major diastereomers, the active (S,S,S)-13b (13b-K), and the nearly inactive (R,S,S)-13b (13b-H); results for the latter reveal a novel binding mode. Pharmacokinetic studies show good levels of 13b-K after inhalative as well as after peroral administration. The active inhibitor (13b-K) is a promising candidate for further development as an antiviral treatment for COVID-19.