Comprehensive Insights into the Catalytic Mechanism of Middle East Respiratory Syndrome 3C-Like Protease and Severe Acute Respiratory Syndrome 3C-Like Protease

A SARS-CoV-2 Mpro fluorescent sensor for exploring pharmacodynamic substances from traditional Chinese medicine

The main protease of SARS-CoV-2 (SARS-CoV-2 Mpro) plays a critical role in the replication and life cycle of the virus. Currently, how to screen SARS-CoV-2 Mpro inhibitors from complex traditional Chinese medicine (TCM) is the bottleneck for exploring the pharmacodynamic substances of TCM against SARS-CoV-2. In this study, a simple, cost-effective, rapid, and selective fluorescent sensor (TPE-S-TLG sensor) was designed with an AIE (aggregation-induced emission) probe (TPE-Ph-In) and the SARS-CoV-2 Mpro substrate (S-TLG). The TPE-S-TLG sensor was characterized using UV-Vis absorption spectroscopy, fluorescence spectroscopy, dynamic light scattering (DLS), transmission electron microscopy (TEM), zeta potential, and Fourier transform infrared (FTIR) spectroscopy techniques. The limit of detection of this method to detect SARS-CoV-2 Mpro was measured to be 5 ng mL-1. Furthermore, the TPE-S-TLG sensor was also successfully applied to screen Mpro inhibitors from Xuebijing injection using the separation and collection of the HPLC-fully automatic partial fraction collector (HPLC-FC). Six active compounds, including protocatechualdehyde, chlorogenic acid, hydroxysafflower yellow A, caffeic acid, isoquercetin, and pentagalloylglucose, were identified using UHPLC-Q-TOF/MS that could achieve 90% of the Mpro inhibition rate for the Xuebijing injection. Accordingly, the strategy can be broadly applied in the detection of disease-related proteases as well as screening active substances from TCM.

Insights into the Activation of Unfolded Protein Response Mechanism during Coronavirus Infection

Coronaviruses represent a significant class of viruses that affect both animals and humans. Their replication cycle is strongly associated with the endoplasmic reticulum (ER), which, upon virus invasion, triggers ER stress responses. The activation of the unfolded protein response (UPR) within infected cells is performed from three transmembrane receptors, IRE1, PERK, and ATF6, and results in a reduction in protein production, a boost in the ER's ability to fold proteins properly, and the initiation of ER-associated degradation (ERAD) to remove misfolded or unfolded proteins. However, in cases of prolonged and severe ER stress, the UPR can also instigate apoptotic cell death and inflammation. Herein, we discuss the ER-triggered host responses after coronavirus infection, as well as the pharmaceutical targeting of the UPR as a potential antiviral strategy.

An orally bioavailable SARS-CoV-2 main protease inhibitor exhibits improved affinity and reduced sensitivity to mutations

Inhibitors of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease (Mpro) such as nirmatrelvir (NTV) and ensitrelvir (ETV) have proven effective in reducing the severity of COVID-19, but the presence of resistance-conferring mutations in sequenced viral genomes raises concerns about future drug resistance. Second-generation oral drugs that retain function against these mutants are thus urgently needed. We hypothesized that the covalent hepatitis C virus protease inhibitor boceprevir (BPV) could serve as the basis for orally bioavailable drugs that inhibit SARS-CoV-2 Mpro more efficiently than existing drugs. Performing structure-guided modifications of BPV, we developed a picomolar-affinity inhibitor, ML2006a4, with antiviral activity, oral pharmacokinetics, and therapeutic efficacy similar or superior to those of NTV. A crucial feature of ML2006a4 is a derivatization of the ketoamide reactive group that improves cell permeability and oral bioavailability. Last, ML2006a4 was found to be less sensitive to several mutations that cause resistance to NTV or ETV and occur in the natural SARS-CoV-2 population. Thus, anticipatory design can preemptively address potential resistance mechanisms to expand future treatment options against coronavirus variants.

Reprogramming the Transition States to Enhance C-N Cleavage Efficiency of Rhodococcus opacusl-Amino Acid Oxidase

l-Amino acid oxidase (LAAO) is an important biocatalyst used for synthesizing α-keto acids. LAAO from Rhodococcus opacus (RoLAAO) has a broad substrate spectrum; however, its low total turnover number limits its industrial use. To overcome this, we aimed to employ crystal structure-guided density functional theory calculations and molecular dynamic simulations to investigate the catalytic mechanism. Two key steps were identified: S → [TS1] in step 1 and Int1 → [TS2] in step 2. We reprogrammed the transition states [TS1] and [TS2] to reduce the identified energy barrier and obtain a RoLAAO variant capable of catalyzing 19 kinds of l-amino acids to the corresponding high-value α-keto acids with a high total turnover number, yield (≥95.1 g/L), conversion rate (≥95%), and space-time yields ≥142.7 g/L/d in 12-24 h, in a 5 L reactor. Our results indicated the promising potential of the developed RoLAAO variant for use in the industrial production of α-keto acids while providing a potential catalytic-mechanism-guided protein design strategy to achieve the desired physical and catalytic properties of enzymes.

AI-driven covalent drug design strategies targeting main protease (mpro) against SARS-CoV-2: structural insights and molecular mechanisms

The emergence of new SARS-CoV-2 variants has raised concerns about the effectiveness of COVID-19 vaccines. To address this challenge, small-molecule antivirals have been proposed as a crucial therapeutic option. Among potential targets for anti-COVID-19 therapy, the main protease (Mpro) of SARS-CoV-2 is important due to its essential role in the virus's life cycle and high conservation. The substrate-binding region of the core proteases of various coronaviruses, including SARS-CoV-2, SARS-CoV, and Middle East respiratory syndrome coronavirus (MERS-CoV), could be used for the generation of new protease inhibitors. Various drug discovery methods have employed a diverse range of strategies, targeting both monomeric and dimeric forms, including drug repurposing, integrating virtual screening with high-throughput screening (HTS), and structure-based drug design, each demonstrating varying levels of efficiency. Covalent inhibitors, such as Nirmatrelvir and MG-101, showcase robust and high-affinity binding to Mpro, exhibiting stable interactions confirmed by molecular docking studies. Development of effective antiviral drugs is imperative to address potential pandemic situations. This review explores recent advances in the search for Mpro inhibitors and the application of artificial intelligence (AI) in drug design. AI leverages vast datasets and advanced algorithms to streamline the design and identification of promising Mpro inhibitors. AI-driven drug discovery methods, including molecular docking, predictive modeling, and structure-based drug repurposing, are at the forefront of identifying potential candidates for effective antiviral therapy. In a time when COVID-19 potentially threat global health, the quest for potent antiviral solutions targeting Mpro could be critical for inhibiting the virus.Communicated by Ramaswamy H. Sarma.

Identification of potential targets against SARS-CoV-2 of antiviral drugs based on photoaffinity probes

Facing the sudden outbreak of coronavirus disease 2019 (COVID-19), it is extremely urgent to develop effective antiviral drugs against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Drug repurposing is a promising strategy for the treatment of COVID-19. To identify the precise target protein of marketed medicines, we initiate a chemical biological program to identify precise target of potential antivirus drugs. In this study, two types of recombinant human coronavirus SARS-CoV-2 RdRp protein capturing probes with various photoaffinity labeling units were designed and synthesized based on the structure of FDA-approved drugs stavudine, remdesivir, acyclovir, and aladenosine. Fortunately, it was found that one novel photoaffinity probe, RD-1, could diaplayed good affinity with SARS-CoV-2 RdRp around the residue ARG_553. In addition, RD-1 probe also exhibited potent inhibitory activity against 3CLpro protease. Taken together, our findings will elucidate the structural basis for the efficacy of marketed drugs, and explore a rapid and efficient strategy of drug repurposing based on the identification of new targets. Moreover, these results could also provide a scientific basis for the clinical application of marketed drugs.

Large library docking for novel SARS-CoV-2 main protease non-covalent and covalent inhibitors

Antiviral therapeutics to treat SARS-CoV-2 are needed to diminish the morbidity of the ongoing COVID-19 pandemic. A well-precedented drug target is the main viral protease (MPro ), which is targeted by an approved drug and by several investigational drugs. Emerging viral resistance has made new inhibitor chemotypes more pressing. Adopting a structure-based approach, we docked 1.2 billion non-covalent lead-like molecules and a new library of 6.5 million electrophiles against the enzyme structure. From these, 29 non-covalent and 11 covalent inhibitors were identified in 37 series, the most potent having an IC50 of 29 and 20 μM, respectively. Several series were optimized, resulting in low micromolar inhibitors. Subsequent crystallography confirmed the docking predicted binding modes and may template further optimization. While the new chemotypes may aid further optimization of MPro inhibitors for SARS-CoV-2, the modest success rate also reveals weaknesses in our approach for challenging targets like MPro versus other targets where it has been more successful, and versus other structure-based techniques against MPro itself.

Identification of a Putative SARS-CoV-2 Main Protease Inhibitor through In Silico Screening of Self-Designed Molecular Library

There have been outbreaks of SARS-CoV-2 around the world for over three years, and its variants continue to evolve. This has become a major global health threat. The main protease (M^pro, also called 3CL^pro) plays a key role in viral replication and proliferation, making it an attractive drug target. Here, we have identified a novel potential inhibitor of M^pro, by applying the virtual screening of hundreds of nilotinib-structure-like compounds that we designed and synthesized. The screened compounds were assessed using SP docking, XP docking, MM-GBSA analysis, IFD docking, MD simulation, ADME/T prediction, and then an enzymatic assay in vitro. We finally identified the compound V291 as a potential SARS-CoV-2 M^pro inhibitor, with a high docking affinity and enzyme inhibitory activity. Moreover, the docking results indicate that His41 is a favorable amino acid for pi-pi interactions, while Glu166 can participate in salt-bridge formation with the protonated primary or secondary amines in the screened molecules. Thus, the compounds reported here are capable of engaging the key amino acids His41 and Glu166 in ligand-receptor interactions. A pharmacophore analysis further validates this assertion.

Comprehensive Understanding of the Kinetic Behaviors of Main Protease from SARS-CoV-2 and SARS-CoV: New Data and Comparison to Published Parameters

The main protease (M^pro) is a promising drug target for inhibiting the coronavirus due to its conserved properties and lack of homologous genes in humans. However, previous studies on M^pro's kinetic parameters have been confusing, hindering the selection of accurate inhibitors. Therefore, obtaining a clear view of M^pro's kinetic parameters is necessary. In our study, we investigated the kinetic behaviors of M^pro from SARS-CoV-2 and SARS-CoV using both FRET-based cleavage assay and the LC-MS method, respectively. Our findings indicate that the FRET-based cleavage assay could be used for preliminary screening of M^pro inhibitors, while the LC-MS method should be applied to select the effective inhibitors with higher reliability. Furthermore, we constructed the active site mutants (H41A and C145A) and measured the kinetic parameters to gain a deeper understanding of the atomic-level enzyme efficiency reduction compared to the wild type. Overall, our study provides valuable insights for inhibitor screening and design by offering a comprehensive understanding of M^pro's kinetic behaviors.

Unraveling the catalytic mechanism of SARS-CoV-2 papain-like protease with allosteric modulation of C270 mutation using multiscale computational approaches

Papain-like protease (PL^pro) is a promising therapeutic target against SARS-CoV-2, but its restricted S1/S2 subsites pose an obstacle in developing active site-directed inhibitors. We have recently identified C270 as a novel covalent allosteric site for SARS-CoV-2 PL^pro inhibitors. Here we present a theoretical investigation of the proteolysis reaction catalyzed by the wild-type SARS-CoV-2 PL^pro as well as the C270R mutant. Enhanced sampling MD simulations were first performed to explore the influence of C270R mutation on the protease dynamics, and sampled thermodynamically favorable conformations were then submitted to MM/PBSA and QM/MM MD simulations for thorough characterization of the protease-substrate binding and covalent reactions. The disclosed proteolysis mechanism of PL^pro, as characterized by the occurrence of proton transfer from the catalytic C111 to H272 prior to the substrate binding and with deacylation being the rate-determining step of the whole proteolysis process, is not completely identical to that of the 3C-like protease, another key cysteine protease of coronaviruses. The C270R mutation alters the structural dynamics of the BL2 loop that indirectly impairs the catalytic function of H272 and reduces the binding of the substrate with the protease, ultimately showing an inhibitory effect on PL^pro. Together, these results provide a comprehensive understanding at the atomic level of the key aspects of SARS-CoV-2 PL^pro proteolysis, including the catalytic activity allosterically regulated by C270 modification, which is crucial to the follow-up inhibitor design and development.

Novel dithiocarbamates selectively inhibit 3CL protease of SARS-CoV-2 and other coronaviruses

Since end of 2019, the global and unprecedented outbreak caused by the coronavirus SARS-CoV-2 led to dramatic numbers of infections and deaths worldwide. SARS-CoV-2 produces two large viral polyproteins which are cleaved by two cysteine proteases encoded by the virus, the 3CL protease (3CL^pro) and the papain-like protease, to generate non-structural proteins essential for the virus life cycle. Both proteases are recognized as promising drug targets for the development of anti-coronavirus chemotherapy. Aiming at identifying broad spectrum agents for the treatment of COVID-19 but also to fight emergent coronaviruses, we focused on 3CL^pro that is well conserved within this viral family. Here we present a high-throughput screening of more than 89,000 small molecules that led to the identification of a new chemotype, potent inhibitor of the SARS-CoV-2 3CL^pro. The mechanism of inhibition, the interaction with the protease using NMR and X-Ray, the specificity against host cysteine proteases and promising antiviral properties in cells are reported.

Effectiveness of ruminal xylanase with an extra proline-rich C-terminus on lignocellulosic biomass degradation

The efficient degradation of plant polysaccharides in agricultural waste requires xylanases with high catalytic activity. In this study, the C-terminal proline-rich GH10 xylanase XynA from sheep rumen was investigated using product analysis, structural characterization, truncated and site-directed mutagenesis, molecular dynamics simulation, and application evaluation, revealing that the proline-rich C-terminus contributes to the interaction at the substrate-binding pocket to reduce the binding free energy. Compared to the C-terminally truncated enzyme XynA-Tr, XynA has a more favorable conformation for proton transfer and affinity attack, facilitating the degradation of oligomeric and beechwood xylan without altering the hydrolysis pattern. Moreover, both the reduced sugar yield and weight loss of the pretreated wheat bran, corn cob, and corn stalk hydrolyzed by XynA for 12 h increased by more than 30 %. These findings are important to better understand the relationship between enzyme activities and their terminal regions and suggest candidate materials for lignocellulosic biomass utilization.

Identifying the natural compound Catechin from tropical mangrove plants as a potential lead candidate against 3CLpro from SARS-CoV-2: An integrated in silico approach

SARS-CoV-2, a member of beta coronaviruses, is a single-stranded, positive-sense RNA virus responsible for the COVID-19 pandemic. With global fatalities of the pandemic exceeding 4.57 million, it becomes crucial to identify effective therapeutics against the virus. A protease, 3CL^pro, is responsible for the proteolysis of viral polypeptides into functional proteins, which is essential for viral pathogenesis. This indispensable activity of 3CL^pro makes it an attractive target for inhibition studies. The current study aimed to identify potential lead molecules against 3CL^pro of SARS-CoV-2 using a manually curated in-house library of antiviral compounds from mangrove plants. This study employed the structure-based virtual screening technique to evaluate an in-house library of antiviral compounds against 3CL^pro of SARS-CoV-2. The library was comprised of thirty-three experimentally proven antiviral molecules extracted from different species of tropical mangrove plants. The molecules in the library were virtually screened using AutoDock Vina, and subsequently, the top five promising 3CL^pro-ligand complexes along with 3CL^pro-N3 (control molecule) complex were subjected to MD simulations to comprehend their dynamic behaviour and structural stabilities. Finally, the MM/PBSA approach was used to calculate the binding free energies of 3CLpro complexes. Among all the studied compounds, Catechin achieved the most significant binding free energy (-40.3 ± 3.1 kcal/mol), and was closest to the control molecule (-42.8 ± 5.1 kcal/mol), and its complex with 3CL^pro exhibited the highest structural stability. Through extensive computational investigations, we propose Catechin as a potential therapeutic agent against SARS-CoV-2. Communicated by Ramaswamy H. Sarma.

Design and in silico investigation of novel Maraviroc analogues as dual inhibition of CCR-5/SARS-CoV-2 Mpro

A sudden increase in life-threatening COVID-19 infections around the world inflicts global crisis and emotional trauma. In current study two druggable targets, namely SARS-COV-2 M^pro and CCR-5 were selected due to their significant nature in the viral life cycle and cytokine molecular storm respectively. The systematic drug repurposing strategy has been utilized to recognize inhibitory mechanism through extensive in silico investigation of novel Maraviroc analogues as promising inhibitors against SARS-CoV-2 M^pro and CCR-5. The dual inhibition specificity approach implemented in present study using molecular docking, molecular dynamics (MD), principal component analysis (PCA), free energy landscape (FEL) and MM/PBSA binding energy studies. The proposed Maraviroc analogues obtained from in silico investigation could be easily synthesized and constructive in developing significant drug against COVID-19 pandemic, with essentiality of their in vivo/in vitro evaluation to affirm the conclusions of this study. This will further fortify the concept of single drug targeting dual inhibition mechanism for treatment of COVID-19 infection and complications.

Evaluation of the anti-diabetic drug sitagliptin as a novel attenuate to SARS-CoV-2 evidence-based in silico: molecular docking and molecular dynamics

The current outbreak of COVID-19 cases worldwide has been responsible for a significant number of deaths, especially in hospitalized patients suffering from comorbidities, such as obesity, diabetes, hypertension. The disease not only has prompted an interest in the pathophysiology, but also it has propelled a massive race to find new anti-SARS-CoV-2 drugs. In this scenario, known drugs commonly used to treat other diseases have been suggested as alternative or complementary therapeutics. Herein we propose the use of sitagliptin, an inhibitor of dipeptidyl peptidase-4 (DPP₄) used to treat type-II diabetes, as an agent to block and inhibit the activity of two proteases, 3CL^pro and PL^pro, related to the processing of SARS-CoV-2 structural proteins. Inhibition of these proteases may possibly reduce the viral load and infection on the host by hampering the synthesis of new viruses, thus promoting a better outcome. In silico assays consisting in the modeling of the ligand sitagliptin and evaluation of its capacity to interact with 3CL^pro and PL^pro through the prediction of the ligand bioactivity, molecular docking, overlapping of crystal structures, and molecular dynamic simulations were conducted. The experiments indicate that sitagliptin can interact and bind to both targets. However, this interaction seems to be stronger and more stable to 3CL^pro (ΔG = −7.8 kcal mol⁻¹), when compared to PL^pro (ΔG = −7.5 kcal mol⁻¹). This study suggests that sitagliptin may be suitable to treat COVID-19 patients, beyond its common use as an anti-diabetic medication. In vivo studies may further support this hypothesis.

Prediction of Recurrent Mutations in SARS-CoV-2 Using Artificial Neural Networks

Predicting SARS-CoV-2 mutations is difficult, but predicting recurrent mutations driven by the host, such as those caused by host deaminases, is feasible. We used machine learning to predict which positions from the SARS-CoV-2 genome will hold a recurrent mutation and which mutations will be the most recurrent. We used data from April 2021 that we separated into three sets: a training set, a validation set, and an independent test set. For the test set, we obtained a specificity value of 0.69, a sensitivity value of 0.79, and an Area Under the Curve (AUC) of 0.8, showing that the prediction of recurrent SARS-CoV-2 mutations is feasible. Subsequently, we compared our predictions with updated data from January 2022, showing that some of the false positives in our prediction model become true positives later on. The most important variables detected by the model's Shapley Additive exPlanation (SHAP) are the nucleotide that mutates and RNA reactivity. This is consistent with the SARS-CoV-2 mutational bias pattern and the preference of some host deaminases for specific sequences and RNA secondary structures. We extend our investigation by analyzing the mutations from the variants of concern Alpha, Beta, Delta, Gamma, and Omicron. Finally, we analyzed amino acid changes by looking at the predicted recurrent mutations in the M-pro and spike proteins.

Ultrafast one-minute electronic detection of SARS-CoV-2 infection by 3CLpro enzymatic activity in untreated saliva samples

Since its onset in December 2019, severe acute respiratory syndrome coronavirus 2, SARS-CoV-2, has caused over 6.5 million deaths worldwide as of October 2022. Attempts to curb viral transmission rely heavily on reliable testing to detect infections since a large number of transmissions are carried through asymptomatic individuals. Many available detection methods fall short in terms of reliability or point-of-care applicability. Here, we report an electrochemical approach targeting a viral proteolytic enzyme, 3CLpro, as a marker of active infection. We detect proteolytic activity directly from untreated saliva within one minute of sample incubation using a reduction-oxidation pH indicator. Importantly, clinical tests of saliva samples from 50 subjects show accurate detection of SARS-CoV-2, with high sensitivity and specificity, validated by PCR testing. These, coupled with our platform's ultrafast detection, simplicity, low cost and point-of-care compatibility, make it a promising method for the real-world SARS-CoV-2 mass-screening.

Recent Drug Development and Medicinal Chemistry Approaches for the Treatment of SARS-CoV-2 Infection and COVID-19

COVID-19, caused by SARS-CoV-2 infection, continues to be a major public health crisis around the globe. Development of vaccines and the first cluster of antiviral drugs has brought promise and hope for prevention and treatment of severe coronavirus disease. However, continued development of newer, safer, and more effective antiviral drugs are critically important to combat COVID-19 and counter the looming pathogenic variants. Studies of the coronavirus life cycle revealed several important biochemical targets for drug development. In the present review, we focus on recent drug design and medicinal chemistry efforts in small molecule drug discovery, including the development of nirmatrelvir that targets viral protein synthesis and remdesivir and molnupiravir that target viral RdRp. These are recent FDA approved drugs for the treatment of COVID-19.

Autoprocessing and oxyanion loop reorganization upon GC373 and nirmatrelvir binding of monomeric SARS-CoV-2 main protease catalytic domain

The monomeric catalytic domain (residues 1–199) of SARS-CoV-2 main protease (MPro^1-199) fused to 25 amino acids of its flanking nsp4 region mediates its autoprocessing at the nsp4-MPro^1-199 junction. We report the catalytic activity and the dissociation constants of MPro^1-199 and its analogs with the covalent inhibitors GC373 and nirmatrelvir (NMV), and the estimated monomer-dimer equilibrium constants of these complexes. Mass spectrometry indicates the presence of the accumulated adduct of NMV bound to MPro^WT and MPro^1-199 and not of GC373. A room temperature crystal structure reveals a native-like fold of the catalytic domain with an unwound oxyanion loop (E state). In contrast, the structure of a covalent complex of the catalytic domain-GC373 or NMV shows an oxyanion loop conformation (E* state) resembling the full-length mature dimer. These results suggest that the E-E* equilibrium modulates autoprocessing of the main protease when converting from a monomeric polyprotein precursor to the mature dimer.

Structural characterization and catalytic activity of SARS-CoV-2 main protease reveal minimal interface regions enabling dimer formation driven by inhibitor-induced conformational changes of the oxyanion loop.

Novel covalent and non-covalent complex-based pharmacophore models of SARS-CoV-2 main protease (Mpro) elucidated by microsecond MD simulations

As the world enters its second year of the pandemic caused by SARS-CoV-2, intense efforts have been directed to develop an effective diagnosis, prevention, and treatment strategies. One promising drug target to design COVID-19 treatments is the SARS-CoV-2 M^pro. To date, a comparative understanding of M^pro dynamic stereoelectronic interactions with either covalent or non-covalent inhibitors (depending on their interaction with a pocket called S1' or oxyanion hole) has not been still achieved. In this study, we seek to fill this knowledge gap using a cascade in silico protocol of docking, molecular dynamics simulations, and MM/PBSA in order to elucidate pharmacophore models for both types of inhibitors. After docking and MD analysis, a set of complex-based pharmacophore models was elucidated for covalent and non-covalent categories making use of the residue bonding point feature. The highest ranked models exhibited ROC-AUC values of 0.93 and 0.73, respectively for each category. Interestingly, we observed that the active site region of M^pro protein-ligand complex undergoes large conformational changes, especially within the S2 and S4 subsites. The results reported in this article may be helpful in virtual screening (VS) campaigns to guide the design and discovery of novel small-molecule therapeutic agents against SARS-CoV-2 M^pro protein.

The structure-based design of peptidomimetic inhibitors against SARS-CoV-2 3C like protease as Potent anti-viral drug candidate

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), as the pathogen of coronavirus disease 2019 (COVID-19), has infected millions of people and took hundreds of thousands of lives. Unfortunately, there is deficiency of effective medicines to prevent or treat COVID-19. 3C like protease (3CL^Pro) of SARS-CoV-2 is essential to the viral replication and transcription, and is an attractive target to develop anti-SARS-CoV-2 agents. Targeting on the 3CL^Pro, we screened our protease inhibitor library and obtained compound 10a as hit to weakly inhibit the SARS-CoV-2 3CL^Pro, and determined the co-crystal structure of 10a and the protease. Based on the deep understanding on the protein-ligand complexes between the hit and SARS-CoV-2 3CL^Pro, we designed a series of peptidomimetic inhibitors, with outstanding inhibitory activity against SARS-CoV-2 3CL^Pro and excellent anti-viral potency against SARS-CoV-2. The protein-ligand complexes of the other key inhibitors with SARS-CoV-2 3CL^Pro were explicitly described by the X-ray co-crystal study. All such results suggest these peptidomimetic inhibitors could be further applied as encouraging drug candidates.

Specific anion effect on properties of HRV 3C protease

Specific salts effect is intensively studied from the prospective of modification of different physico-chemical properties of biomacromolecules. Limited knowledge of the specific salts effect on enzymes led us to address the influence of five sodium anions: sulfate, phosphate, chloride, bromide, and perchlorate, on catalytic and conformational properties of human rhinovirus-14 (HRV) 3C protease. The enzyme conformation was monitored by circular dichroism spectrum (CD) and by tyrosines fluorescence. Stability and flexibility of the enzyme have been analyzed by CD in the far-UV region, differential scanning calorimetry and molecular dynamics simulations, respectively. We showed significant influence of the anions on the enzyme properties in accordance with the Hofmeister effect. The HRV 3C protease in the presence of kosmotropic anions, in contrast with chaotropic anions, exhibits increased stability, rigidity. Correlations of stabilization effect of anions on the enzyme with their charge density and the rate constant of the enzyme with the viscosity B-coefficients of anions suggest direct interaction of the anions with HRV 3C protease. The role of stabilization and decreased fluctuation of the polypeptide chain of HRV 3C protease on its activation in the presence of kosmotropic anions is discussed within the frame of the macromolecular rate theory.

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.

Maturation of the SARS-CoV-2 virus is regulated by dimerization of its main protease

SARS-CoV-2 main protease (M^pro) involved in COVID-19 is required for maturation of the virus and infection of host cells. The key question is how to block the activity of M^pro. By combining atomistic simulations with machine learning, we found that the enzyme regulates its own activity by a collective allosteric mechanism that involves dimerization and binding of a single substrate. At the core of the collective mechanism is the coupling between the catalytic site residues, H41 and C145, which direct the activity of M^pro dimer, and two salt bridges formed between R4 and E290 at the dimer interface. If these salt bridges are mutated, the activity of M^pro is blocked. The results suggest that dimerization of main proteases is a general mechanism to foster coronavirus proliferation, and propose a robust drug-based strategy that does not depend on the frequently mutating spike proteins at the viral envelope used to develop vaccines.

Multiple protonation states in ligand-free SARS-CoV-2 main protease revealed by large-scale quantum molecular dynamics simulations

The main protease (M^pro) in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) catalyzes the cleavage of polyproteins for viral replication. Here, large-scale quantum molecular dynamics and metadynamics simulations for ligand-free M^pro were performed, where all the atoms were treated quantum-mechanically, focusing on elucidation of the controversial active-site protonation state. The simulations clarified that the interconverting multiple protonation states exist in unliganded M^pro, and the catalytically relevant ion-pair state is more stable than the neutral state, which is consistent with neutron crystallography. The results highlight the importance of the ion-pair state for repurposing or discovering antiviral drugs that target M^pro.

Covalent narlaprevir- and boceprevir-derived hybrid inhibitors of SARS-CoV-2 main protease

Emerging SARS-CoV-2 variants continue to threaten the effectiveness of COVID-19 vaccines, and small-molecule antivirals can provide an important therapeutic treatment option. The viral main protease (M^pro) is critical for virus replication and thus is considered an attractive drug target. We performed the design and characterization of three covalent hybrid inhibitors BBH-1, BBH-2 and NBH-2 created by splicing components of hepatitis C protease inhibitors boceprevir and narlaprevir, and known SARS-CoV-1 protease inhibitors. A joint X-ray/neutron structure of the M^pro/BBH-1 complex demonstrates that a Cys145 thiolate reaction with the inhibitor’s keto-warhead creates a negatively charged oxyanion. Protonation states of the ionizable residues in the M^pro active site adapt to the inhibitor, which appears to be an intrinsic property of M^pro. Structural comparisons of the hybrid inhibitors with PF-07321332 reveal unconventional F···O interactions of PF-07321332 with M^pro which may explain its more favorable enthalpy of binding. BBH-1, BBH-2 and NBH-2 exhibit comparable antiviral properties in vitro relative to PF-07321332, making them good candidates for further design of improved antivirals.

Three covalent hybrid inhibitors of SARS-CoV-2 main protease (Mpro) have been designed and compared to Pfizer’s nirmatrelvir (PF-07321332), providing atomic and thermodynamic details of their binding to the enzyme, and antiviral potency.

Discovery of 4'-O-methylscutellarein as a potent SARS-CoV-2 main protease inhibitor

The coronavirus disease 2019 (COVID-19) pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has resulted in millions of deaths and seriously threatened public health and safety. Despite COVID-19 vaccines being readily popularized worldwide, targeted therapeutic agents for the treatment of this disease remain very limited. Here, we studied the inhibitory activity of the scutellarein and its methylated derivatives against SARS-CoV-2 main protease (M^pro) by the fluorescence resonance energy transfer (FRET) assay. Among all the methylated derivatives we studied, 4'-O-methylscutellarein exhibited the most promising enzyme inhibitory activity in vitro, with the half-maximal inhibitory concentration value (IC₅₀) of 0.40 ± 0.03 μM. Additionally, the mechanism of action of the hits was further characterized through enzyme kinetic studies and molecular docking. Overall, our results implied that 4'-O-methylscutellarein could be a primary lead compound with clinical potential for the development of inhibitors against the SARS-CoV-2 M^pro.

Peptidomimetic α-Acyloxymethylketone Warheads with Six-Membered Lactam P1 Glutamine Mimic: SARS-CoV-2 3CL Protease Inhibition, Coronavirus Antiviral Activity, and in Vitro Biological Stability

Recurring coronavirus outbreaks, such as the current COVID-19 pandemic, establish a necessity to develop direct-acting antivirals that can be readily administered and are active against a broad spectrum of coronaviruses. Described in this Article are novel α-acyloxymethylketone warhead peptidomimetic compounds with a six-membered lactam glutamine mimic in P1. Compounds with potent SARS-CoV-2 3CL protease and in vitro viral replication inhibition were identified with low cytotoxicity and good plasma and glutathione stability. Compounds 15e, 15h, and 15l displayed selectivity for SARS-CoV-2 3CL protease over CatB and CatS and superior in vitro SARS-CoV-2 antiviral replication inhibition compared with the reported peptidomimetic inhibitors with other warheads. The cocrystallization of 15l with SARS-CoV-2 3CL protease confirmed the formation of a covalent adduct. α-Acyloxymethylketone compounds also exhibited antiviral activity against an alphacoronavirus and non-SARS betacoronavirus strains with similar potency and a better selectivity index than remdesivir. These findings demonstrate the potential of the substituted heteroaromatic and aliphatic α-acyloxymethylketone warheads as coronavirus inhibitors, and the described results provide a basis for further optimization.

Discovery of SARS-CoV-2 3CLPro Peptidomimetic Inhibitors through the Catalytic Dyad Histidine-Specific Protein-Ligand Interactions

As the etiological agent for the coronavirus disease 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) challenges the ongoing efforts of vaccine development and drug design. Due to the accumulating cases of breakthrough infections, there are urgent needs for broad-spectrum antiviral medicines. Here, we designed and examined five new tetrapeptidomimetic anti-SARS-CoV-2 inhibitors targeting the 3C-Like protease (3CL^Pro), which is highly conserved among coronaviruses and essential for viral replications. We significantly improved the efficacy of a ketoamide lead compound based on high-resolution co-crystal structures, all-atom simulations, and binding energy calculations. The inhibitors successfully engaged the catalytic dyad histidine residue (H41) of 3CLPro as designed, and they exhibited nanomolar inhibitory capacity as well as mitigated the viral loads of SARS-CoV-2 in cellular assays. As a widely applicable design principle, our results revealed that the potencies of 3CL^Pro-specific drug candidates were determined by the interplay between 3CL^Pro H41 residue and the peptidomimetic inhibitors.

In silico screening-based discovery of novel covalent inhibitors of the SARS-CoV-2 3CL protease

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) 3CL protease (3CL^pro) has been regarded as an extremely promising antiviral target for the treatment of coronavirus disease 2019 (COVID-19). Here, we carried out a virtual screening based on commercial compounds database to find novel covalent non-peptidomimetic inhibitors of this protease. It allowed us to identify 3 hit compounds with potential covalent binding modes, which were evaluated through an enzymatic activity assay of the SARS-CoV-2 3CL^pro. Moreover, an X-ray crystal structure of the SARS-CoV-2 3CL^pro in complex with compound 8, the most potent hit with an IC₅₀ value of 8.50 μM, confirmed the covalent binding of the predicted warhead to the catalytic residue C145, as well as portrayed interactions of the compound with S1’ and S2 subsites at the ligand binding pocket. Overall, the present work not merely provided an experiment-validated covalent hit targeting the SARS-CoV-2 3CL^pro, but also displayed a prime example to seeking new covalent small molecules by a feasible and effective computational approach.

Perspectives on SARS-CoV-2 Main Protease Inhibitors

The main protease (Mpro) plays a crucial role in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) replication and is highly conserved, rendering it one of the most attractive therapeutic targets for SARS-CoV-2 inhibition. Currently, although two drug candidates targeting SARS-CoV-2 Mpro designed by Pfizer are under clinical trials, no SARS-CoV-2 medication is approved due to the long period of drug development. Here, we collect a comprehensive list of 817 available SARS-CoV-2 and SARS-CoV Mpro inhibitors from the literature or databases and analyze their molecular mechanisms of action. The structure-activity relationships (SARs) among each series of inhibitors are discussed. Additionally, we broadly examine available antiviral activity, ADMET (absorption, distribution, metabolism, excretion, and toxicity), and animal tests of these inhibitors. We comment on their druggability or drawbacks that prevent them from becoming drugs. This Perspective sheds light on the future development of Mpro inhibitors for SARS-CoV-2 and future coronavirus diseases.

Structural and functional characterization of NEMO cleavage by SARS-CoV-2 3CLpro.. [PMID: 34816264 PMCID: PMC8609902 DOI: 10.1101/2021.11.11.468228] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

In addition to its essential role in viral polyprotein processing, the SARS-CoV-2 3C-like (3CLpro) protease can cleave human immune signaling proteins, like NF-κB Essential Modulator (NEMO) and deregulate the host immune response. Here, in vitro assays show that SARS-CoV-2 3CLpro cleaves NEMO with fine-tuned efficiency. Analysis of the 2.14 Å resolution crystal structure of 3CLpro C145S bound to NEMO_226–235 reveals subsites that tolerate a range of viral and host substrates through main chain hydrogen bonds while also enforcing specificity using side chain hydrogen bonds and hydrophobic contacts. Machine learning- and physics-based computational methods predict that variation in key binding residues of 3CLpro-NEMO helps explain the high fitness of SARS-CoV-2 in humans. We posit that cleavage of NEMO is an important piece of information to be accounted for in the pathology of COVID-19.

Michaelis-like complex of SARS-CoV-2 main protease visualized by room-temperature X-ray crystallography

SARS-CoV-2 emerged at the end of 2019 to cause an unprecedented pandemic of the deadly respiratory disease COVID-19 that continues to date. The viral main protease (Mpro) is essential for SARS-CoV-2 replication and is therefore an important drug target. Understanding the catalytic mechanism of Mpro, a cysteine protease with a catalytic site comprising the noncanonical Cys145-His41 dyad, can help in guiding drug design. Here, a 2.0 Å resolution room-temperature X-ray crystal structure is reported of a Michaelis-like complex of Mpro harboring a single inactivating mutation C145A bound to the octapeptide Ac-SAVLQSGF-CONH2 corresponding to the nsp4/nsp5 autocleavage site. The peptide substrate is unambiguously defined in subsites S5 to S3' by strong electron density. Superposition of the Michaelis-like complex with the neutron structure of substrate-free Mpro demonstrates that the catalytic site is inherently pre-organized for catalysis prior to substrate binding. Induced fit to the substrate is driven by P1 Gln binding in the predetermined subsite S1 and rearrangement of subsite S2 to accommodate P2 Leu. The Michaelis-like complex structure is ideal for in silico modeling of the SARS-CoV-2 Mpro catalytic mechanism.

Potential SARS-CoV-2 3CLpro inhibitors from chromene, flavonoid and hydroxamic acid compound based on FRET assay, docking and pharmacophore studies

This present study reports some natural products and one hydroxamic acid synthetic compound which were previously reported as matrix metalloproteinase-9 (MMP-9) inhibitors to be evaluated for their inhibition toward severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) 3-chymotrypsin-like protease (3CLpro). This enzyme is one of the proteins responsible for this coronaviral replication. Two herbal methanolic extracts i.e., Averrhoa carambola leaves and Ageratum conyzoides aerial part demonstrate >50% inhibition at 1000 µg/mL. Interestingly, apigenin, one of flavonoids, demonstrates 92% inhibition at 250 µg/mL (925 µM) as well as hydroxamic acid compound, N-isobutyl-N-(4-methoxyphenylsulfonyl)glycyl hydroxamic acid (NNGH), which shows 69% inhibition at 100 µM. The in vitro results are supported by the docking studies revealing that the binding mode of both compounds is mainly by interacting with GLU166 residue in the hydrophobic pocket of the 3CLpro. Pharmacophore mapping further supported the results by confirming that the in vitro activities of both compounds are due to their pharmacophore features employing hydrogen bond acceptor (HBA), hydrogen bond donor (HBD) and hydrophobic. Gas Chromatography-Mass Spectrometry (GC–MS) analysis reported chromene compounds in Ageratum conyzoides aerial part methanolic extract are potential to be this enzyme inhibitor candidate. These all results reflect their potencies to be SARS-CoV-2 inhibitors through 3CLpro inhibition mechanism.

Near-physiological-temperature serial crystallography reveals conformations of SARS-CoV-2 main protease active site for improved drug repurposing

The COVID-19 pandemic has resulted in 198 million reported infections and more than 4 million deaths as of July 2021 (covid19.who.int). Research to identify effective therapies for COVID-19 includes: (1) designing a vaccine as future protection; (2) de novo drug discovery; and (3) identifying existing drugs to repurpose them as effective and immediate treatments. To assist in drug repurposing and design, we determine two apo structures of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease at ambient temperature by serial femtosecond X-ray crystallography. We employ detailed molecular simulations of selected known main protease inhibitors with the structures and compare binding modes and energies. The combined structural and molecular modeling studies not only reveal the dynamics of small molecules targeting the main protease but also provide invaluable opportunities for drug repurposing and structure-based drug design strategies against SARS-CoV-2.

Neutron crystallography for the elucidation of enzyme catalysis

Hydrogen atoms and hydration water molecules in proteins are indispensable for many biochemical processes, especially enzymatic catalysis. The locations of hydrogen atoms in proteins are usually predicted based on X-ray structures, but it is still very difficult to know the ionization states of the catalytic residues, the hydration structure of the protein, and the characteristics of hydrogen-bonding interactions. Neutron crystallography allows the direct observation of hydrogen atoms that play crucial roles in molecular recognition and the catalytic reactions of enzymes. In this review, we present the current status of neutron crystallography in structural biology and recent neutron structural analyses of three enzymes: ascorbate peroxidase, the main protease of severe acute respiratory syndrome coronavirus 2, and copper-containing nitrite reductase.

Benchmarking the Ability of Common Docking Programs to Correctly Reproduce and Score Binding Modes in SARS-CoV-2 Protease Mpro

The coronavirus SARS-CoV-2 main protease, Mpro, is conserved among coronaviruses with no human homolog and has therefore attracted significant attention as an enzyme drug target for COVID-19. The number of studies targeting Mpro for in silico screening has grown rapidly, and it would be of great interest to know in advance how well docking methods can reproduce the correct ligand binding modes and rank these correctly. Clearly, current attempts at designing drugs targeting Mpro with the aid of computational docking would benefit from a priori knowledge of the ability of docking programs to predict correct binding modes and score these correctly. In the current work, we tested the ability of several leading docking programs, namely, Glide, DOCK, AutoDock, AutoDock Vina, FRED, and EnzyDock, to correctly identify and score the binding mode of Mpro ligands in 193 crystal structures. None of the codes were able to correctly identify the crystal structure binding mode (lowest energy pose with root-mean-square deviation < 2 Å) in more than 26% of the cases for noncovalently bound ligands (Glide: top performer), whereas for covalently bound ligands the top score was 45% (EnzyDock). These results suggest that one should perform in silico campaigns of Mpro with care and that more comprehensive strategies including ligand free energy perturbation might be necessary in conjunction with virtual screening and docking.

New insights into the catalytic mechanism of the SARS-CoV-2 main protease: an ONIOM QM/MM approach

SARS-CoV-2 M^pro, also known as the main protease or 3C-like protease, is a key enzyme involved in the replication process of the virus that is causing the COVID-19 pandemic. It is also the most promising antiviral drug target targeting SARS-CoV-2 virus. In this work, the catalytic mechanism of M^pro was studied using the full model of the enzyme and a computational QM/MM methodology with a 69/72-atoms QM region treated at DLPNO-CCSD(T)/CBS//B3LYP/6-31G(d,p):AMBER level and including the catalytic important oxyanion-hole residues. The transition state of each step was fully characterized and described together with the related reactants and products. The rate-limiting step of the catalytic process is the hydrolysis of the thioester-enzyme adduct, and the calculated barrier closely agrees with the available kinetic data. The calculated Gibbs free energy profile, together with the full atomistic detail of the structures involved in catalysis, can now serve as valuable models for the rational drug design of transition state analogs as new inhibitors targeting the SARS-CoV-2 virus.

Tuning Proton Transfer Thermodynamics in SARS-CoV-2 Main Protease: Implications for Catalysis and Inhibitor Design

The catalytic reaction in SARS-CoV-2 main protease is activated by a proton transfer (PT) from Cys145 to His41. The same PT is likely also required for the covalent binding of some inhibitors. Here we use a multiscale computational approach to investigate the PT thermodynamics in the apo enzyme and in complex with two potent inhibitors, N3 and the α-ketoamide 13b. We show that with the inhibitors the free energy cost to reach the charge-separated state of the active-site dyad is lower, with N3 inducing the most significant reduction. We also show that a few key sites (including specific water molecules) significantly enhance or reduce the thermodynamic feasibility of the PT reaction, with selective desolvation of the active site playing a crucial role. The approach presented is a cost-effective procedure to identify the enzyme regions that control the activation of the catalytic reaction and is thus also useful to guide the design of inhibitors.

Potential 3-chymotrypsin-like cysteine protease cleavage sites in the coronavirus polyproteins pp1a and pp1ab and their possible relevance to COVID-19 vaccine and drug development

Coronavirus (CoV) 3-chymotrypsin (C)-like cysteine protease (3CLpro ) is a target for anti-CoV drug development and drug repurposing because along with papain-like protease, it cleaves CoV-encoded polyproteins (pp1a and pp1ab) into nonstructural proteins (nsps) for viral replication. However, the cleavage sites of 3CLpro and their relevant nsps remain unclear, which is the subject of this perspective. Here, we address the subject from three standpoints. First, we explore the inconsistency in the cleavage sites and relevant nsps across CoVs, and investigate the function of nsp11. Second, we consider the nsp16 mRNA overlapping of the spike protein mRNA, and analyze the effect of this overlapping on mRNA vaccines. Finally, we study nsp12, whose existence depends on ribosomal frameshifting, and investigate whether 3CLpro requires a large number of inhibitors to achieve full inhibition. This perspective helps us to clarify viral replication and is useful for developing anti-CoV drugs with 3CLpro as a target in the current coronavirus disease 2019 (COVID-19) pandemic.

Insights into SARS-CoV-2: Medicinal Chemistry Approaches to Combat Its Structural and Functional Biology

Coronavirus disease 2019, caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is still a pandemic around the world. Currently, specific antiviral drugs to control the epidemic remain deficient. Understanding the details of SARS-CoV-2 structural biology is extremely important for development of antiviral agents that will enable regulation of its life cycle. This review focuses on the structural biology and medicinal chemistry of various key proteins (Spike, ACE2, TMPRSS2, RdRp and M^pro) in the life cycle of SARS-CoV-2, as well as their inhibitors/drug candidates. Representative broad-spectrum antiviral drugs, especially those against the homologous virus SARS-CoV, are summarized with the expectation they will drive the development of effective, broad-spectrum inhibitors against coronaviruses. We are hopeful that this review will be a useful aid for discovery of novel, potent anti-SARS-CoV-2 drugs with excellent therapeutic results in the near future.

Nano-size dependence in the adsorption by the SARS-CoV-2 spike protein over gold colloid

Gold nano-particles were coated with the spike protein (S protein) of SARS-CoV-2 and exposed to increasingly acidic conditions. Their responses were investigated by monitoring the surface plasmon resonance (SPR) band shift. As the external pH was gradually changed from neutral pH to pH ∼2 the peak of the SPR band showed a significant red-shift, with a sigmoidal feature implying the formation of the gold-protein aggregates. The coating of S protein changed the surface property of the gold enough to extract the coverage fraction of protein over nano particles, Θ, which did not exhibit clear nano-size dependence. The geometrical simulation to explain Θ showed the average axial length to be a = 7. 25 nm and b =8.00 nm when the S-protein was hypothesized as a prolate shape with spiking-out orientation. As the pH value externally hopped between pH∼3 and pH∼10, a behavior of reversible protein folding was observed for particles with diameters >30 nm. It was concluded that S protein adsorption conformation was impacted by the size (diameter, d) of a core nano-gold, where head-to-head dimerized S protein was estimated for d ≤ 80 nm and a parallel in opposite directions formation for d = 100 nm.