Structural basis for the in vitro efficacy of nirmatrelvir against SARS-CoV-2 variants.

The COVID-19 pandemic continues to be a public health threat with emerging variants of SARS-CoV-2. Nirmatrelvir (PF-07321332) is a reversible, covalent inhibitor targeting the main protease (Mpro) of SARS-CoV-2 and the active protease inhibitor in PAXLOVID (nirmatrelvir tablets and ritonavir tablets). However, the efficacy of nirmatrelvir is underdetermined against evolving SARS-CoV-2 variants. Here, we evaluated the in vitro catalytic activity and potency of nirmatrelvir against the Mpro of prevalent variants of concern (VOCs) or variants of interest (VOIs): Alpha (α, B.1.1.7), Beta (ß, B.1.351), Delta (δ, B1.617.2), Gamma (γ, P.1), Lambda (λ, B.1.1.1.37/C37), Omicron (ο, B.1.1.529), as well as the original Washington or wildtype strain. These VOCs/VOIs carry prevalent mutations at varying frequencies in the Mpro specifically for α, ß, γ (K90R), λ (G15S), and ο (P132H). In vitro biochemical enzymatic assay characterization of the enzyme kinetics of the mutant Mpros demonstrates that they are catalytically comparable to wildtype. We found that nirmatrelvir has similar potency against each mutant Mpro including P132H that is observed in the Omicron variant with a Ki of 0.635 nM as compared to a Ki of 0.933 nM for wildtype. The molecular basis for these observations were provided by solution-phase structural dynamics and structural determination of nirmatrelvir bound to the ο, λ, and ß Mpro at 1.63 to 2.09 Å resolution. These in vitro data suggest that PAXLOVID has the potential to maintain plasma concentrations of nirmatrelvir many-fold times higher than the amount required to stop the SARS-CoV-2 VOC/VOI, including Omicron, from replicating in cells.

Mechanism of Caspase-1 Inhibition by Four Anti-inflammatory Drugs Used in COVID-19 Treatment.

The inflammatory protease caspase-1 is associated with the release of cytokines. An excessive number of cytokines (a "cytokine storm") is a dangerous consequence of COVID-19 infection and has been indicated as being among the causes of death by COVID-19. The anti-inflammatory drug colchicine (which is reported in the literature to be a caspase-1 inhibitor) and the corticosteroid drugs, dexamethasone and methylprednisolone, are among the most effective active compounds for COVID-19 treatment. The SERM raloxifene has also been used as a repurposed drug in COVID-19 therapy. In this study, inhibition of caspase-1 by these four compounds was analyzed using computational methods. Our aim was to see if the inhibition of caspase-1, an important biomolecule in the inflammatory response that triggers cytokine release, could shed light on how these drugs help to alleviate excessive cytokine production. We also measured the antioxidant activities of dexamethasone and colchicine when scavenging the superoxide radical using cyclic voltammetry methods. The experimental findings are associated with caspase-1 active site affinity towards these compounds. In evaluating our computational and experimental results, we here formulate a mechanism for caspase-1 inhibition by these drugs, which involves the active site amino acid Cys285 residue and is mediated by a transfer of protons, involving His237 and Ser339. It is proposed that the molecular moiety targeted by all of these drugs is a carbonyl group which establishes a S(Cys285)-C(carbonyl) covalent bond.

Binding mechanism of inhibitors to SARS-CoV-2 main protease deciphered by multiple replica molecular dynamics simulations.

The outbreak caused by SARS-CoV-2 has received extensive worldwide attention. As the main protease (Mpro) in SARS-CoV-2 has no human homologues, it is feasible to reduce the possibility of targeting the host protein by accidental drugs. Thus, Mpro has been an attractive target of efficient drug design for anti-SARS-CoV-2 treatment. In this work, multiple replica molecular dynamics (MRMD) simulations, principal component analysis (PCA), free energy landscapes (FELs), and the molecular mechanics-generalized Born surface area (MM-GBSA) method were integrated together to decipher the binding mechanism of four inhibitors masitinib, O6K, FJC and GQU to Mpro. The results indicate that the binding of four inhibitors clearly affects the structural flexibility and internal dynamics of Mpro along with dihedral angle changes of key residues. The analysis of FELs unveils that the stability in the relative orientation and geometric position of inhibitors to Mpro is favorable for inhibitor binding. Residue-based free energy decomposition reveals that the inhibitor-Mpro interaction networks involving hydrogen bonding interactions and hydrophobic interactions provide significant information for the design of potent inhibitors against Mpro. The hot spot residues including H41, M49, F140, N142, G143, C145, H163, H164, M165, E166 and Q189 identified by computational alanine scanning are considered as reliable targets of clinically available inhibitors inhibiting the activities of Mpro.

Isolation and characterization of ACE-I inhibitory peptides from ribbonfish for a potential inhibitor of the main protease of SARS-CoV-2: An in silico analysis.

Recently, multifunctional fish peptides (FWPs) have gained a lot of attention because of their different biological activities. In the present study, three angiotensin-I converting enzyme (ACE-I) inhibitory peptides [Ala-Pro-Asp-Gly (APDG), Pro-Thr-Arg (PTR), and Ala-Asp (AD)] were isolated and characterized from ribbonfish protein hydrolysate (RFPH) and described their mechanism of action on ACE activity. As per the results, peptide PTR showed ≈ 2 and 2.5-fold higher enzyme inhibitory activity (IC50 = 0.643 ± 0.0011 µM) than APDG (IC50 = 1.061 ± 0.0127 µM) and AD (IC50 = 2.046 ± 0.0130 µM). Based on experimental evidence, peptides were used for in silico analysis to check the inhibitory activity of the main protease (PDB: 7BQY) of SARS-CoV-2. The results of the study reveal that PTR (-46.16 kcal/mol) showed higher binding affinity than APDG (-36.80 kcal/mol) and AD (-30.24 kcal/mol) compared with remdesivir (-30.64 kcal/mol). Additionally, physicochemical characteristics of all the isolated peptides exhibited appropriate pharmacological properties and were found to be nontoxic. Besides, 20 ns molecular dynamic simulation study confirms the rigid nature, fewer confirmation variations, and binding stiffness of the peptide PTR with the main protease of SARS-CoV-2. Therefore, the present study strongly suggested that PTR is the perfect substrate for inhibiting the main protease of SARS-CoV-2 through the in silico study, and this potential drug candidate may promote the researcher for future wet lab experiments.

In Silico Evaluation of Iranian Medicinal Plant Phytoconstituents as Inhibitors against Main Protease and the Receptor-Binding Domain of SARS-CoV-2.

The novel coronavirus disease 2019 (COVID-19) is caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which initially appeared in Wuhan, China, in December 2019. Elderly individuals and those with comorbid conditions may be more vulnerable to this disease. Consequently, several research laboratories continue to focus on developing drugs to treat this infection because this disease has developed into a global pandemic with an extremely limited number of specific treatments available. Natural herbal remedies have long been used to treat illnesses in a variety of cultures. Modern medicine has achieved success due to the effectiveness of traditional medicines, which are derived from medicinal plants. The objective of this study was to determine whether components of natural origin from Iranian medicinal plants have an antiviral effect that can prevent humans from this coronavirus infection using the most reliable molecular docking method; in our case, we focused on the main protease (Mpro) and a receptor-binding domain (RBD). The results of molecular docking showed that among 169 molecules of natural origin from common Iranian medicinal plants, 20 molecules (chelidimerine, rutin, fumariline, catechin gallate, adlumidine, astragalin, somniferine, etc.) can be proposed as inhibitors against this coronavirus based on the binding free energy and type of interactions between these molecules and the studied proteins. Moreover, a molecular dynamics simulation study revealed that the chelidimerine-Mpro and somniferine-RBD complexes were stable for up to 50 ns below 0.5 nm. Our results provide valuable insights into this mechanism, which sheds light on future structure-based designs of high-potency inhibitors for SARS-CoV-2.

A Novel Class of Norovirus Inhibitors Targeting the Viral Protease with Potent Antiviral Activity In Vitro and In Vivo.

Human noroviruses (HuNoVs) are the most common cause of viral gastroenteritis resulting annually in ~219,000 deaths and a societal cost of ~USD 60 billion, and no antivirals or vaccines are available. Here, we assess the anti-norovirus activity of new peptidomimetic aldehydes related to the protease inhibitor rupintrivir. The early hit compound 4 inhibited the replication of murine norovirus (MNV) and the HuNoV GI.1 replicon in vitro (EC50 ~1 µM) and swiftly cleared the HuNoV GI.1 replicon from the cells. Compound 4 still inhibits the proteolytic activity. We selected a resistant GI.1 replicon, with a mutation (I109V) in a highly conserved region of the viral protease, conferring a low yield of resistance against compound 4 and rupintrivir. After testing new derivatives, compound 10d was the most potent (EC50 nanomolar range). Molecular docking indicated that the aldehyde group of compounds 4 and 10d bind with Cys139 in the HuNoV 3CL protease by a covalent linkage. Finally, compound 10d inhibited the replication of HuNoV GII.4 in infected zebrafish larvae, and PK studies in mice showed an adequate profile.

Fragment-based in silico design of SARS-CoV-2 main protease inhibitors.

3CLpro is essential for SARS-CoV-2 replication and infection; its inhibition using small molecules is a potential therapeutic strategy. In this study, a comprehensive crystallography-guided fragment-based drug discovery approach was employed to design new inhibitors for SARS-CoV-2 3CLpro. All small molecules co-crystallized with SARS-CoV-2 3CLpro with structures deposited in the Protein Data Bank were used as inputs. Fragments sitting in the binding pocket (87) were grouped into eight geographical types. They were interactively coupled using various synthetically reasonable linkers to generate larger molecules with divalent binding modes taking advantage of two different fragments' interactions. In total, 1,251 compounds were proposed, and 7,158 stereoisomers were screened using Glide (standard precision and extra precision), AutoDock Vina, and Prime MMGBSA. The top 22 hits having conformations approaching the linear combination of their constituent fragments were selected for MD simulation on Desmond. MD simulation suggested 15 of these did adopt conformations very close to their constituent pieces with far higher binding affinity than either constituent domain alone. These structures could provide a starting point for the further design of SARS-CoV-2 3CLpro inhibitors with improved binding, and structures are provided.

Screening potential FDA-approved inhibitors of the SARS-CoV-2 major protease 3CLpro through high-throughput virtual screening and molecular dynamics simulation.

It has been confirmed that the new coronavirus SARS-CoV-2 caused the global pandemic of coronavirus disease 2019 (COVID-19). Studies have found that 3-chymotrypsin-like protease (3CLpro) is an essential enzyme for virus replication, and could be used as a potential target to inhibit SARS-CoV-2. In this work, 3CLpro was used as the target to complete the high-throughput virtual screening of the FDA-approved drugs, and Indinavir and other 10 drugs with high docking scores for 3CLpro were obtained. Studies on the binding pattern of 3CLpro and Indinavir found that Indinavir could form the stable hydrogen bond (H-bond) interactions with the catalytic dyad residues His41-Cys145. Binding free energy study found that Indinavir had high binding affinity with 3CLpro. Subsequently, molecular dynamics simulations were performed on the 3CLpro and 3CLpro-Indinavir systems, respectively. The post-dynamic analyses showed that the conformational state of the 3CLpro-Indinavir system transformed significantly and the system tended to be more stable. Moreover, analyses of the residue interaction network (RIN) and H-bond occupancy revealed that the residue-residue interaction at the catalytic site of 3CLpro was significantly enhanced after binding with Indinavir, which in turn inactivated the protein. In short, through this research, we hope to provide more valuable clues against COVID-19.

Hispidin and Lepidine E: Two Natural Compounds and Folic Acid as Potential Inhibitors of 2019-novel Coronavirus Main Protease (2019- nCoVMpro), Molecular Docking and SAR Study.

BACKGROUND: 2019-nCoVis, a novel coronavirus was isolated and identified in 2019 in the city of Wuhan, China. On February 17, 2020 and according to the World Health Organization, 71, 429 confirmed cases worldwide were identified, among them 2162 new cases were recorded in the last 24 hours. One month later, the confirmed cases jumped to 179111, with 11525 new cases in the last 24 hours, with 7426 total deaths. No drug or vaccine is present at the moment for human and animal coronavirus. METHODS: The inhibition of 3CL hydrolase enzyme provides a promising therapeutic principle for developing treatments against CoViD-19. The 3CLpro (Mpro) is known for involving in counteracting the host innate immune response. RESULTS: This work presents the inhibitory effect of some natural compounds against 3CL hydrolase enzyme, and explains the main interactions in inhibitor-enzyme complex. Molecular docking study was carried out using Autodock Vina. By screening several molecules, we identified three candidate agents that inhibit the main protease of coronavirus. Hispidin, lepidine E, and folic acid are bound tightly in the enzyme, therefore strong hydrogen bonds have been formed (1.69-1.80Å) with the active site residues. CONCLUSION: This study provides a possible therapeutic strategy for CoViD-19.

PRIMoRDiA: A Software to Explore Reactivity and Electronic Structure in Large Biomolecules.

Plenty of enzymes with structural data do not have their mechanism of catalysis elucidated. Reactivity descriptors, theoretical quantities generated from resolved electronic structure, provide a way to predict and rationalize chemical processes of such systems. In this Application Note, we present PRIMoRDiA (PRIMoRDiA Macromolecular Reactivity Descriptors Access), a software built to calculate the reactivity descriptors of large biosystems by employing an efficient and accurate treatment of the large output files produced by quantum chemistry packages. Here, we show the general implementation details and the software main features. Calculated descriptors were applied for a set of enzymatic systems in order to show their relevance for biological studies and the software potential for use in large scale. Also, we test PRIMoRDiA to aid in the interaction depiction between the SARS-CoV-2 main protease and a potential inhibitor.

Interactive Molecular Dynamics in Virtual Reality Is an Effective Tool for Flexible Substrate and Inhibitor Docking to the SARS-CoV-2 Main Protease.

The main protease (Mpro) of the SARS-CoV-2 virus is one focus of drug development efforts for COVID-19. Here, we show that interactive molecular dynamics in virtual reality (iMD-VR) is a useful and effective tool for creating Mpro complexes. We make these tools and models freely available. iMD-VR provides an immersive environment in which users can interact with MD simulations and so build protein complexes in a physically rigorous and flexible way. Recently, we have demonstrated that iMD-VR is an effective method for interactive, flexible docking of small molecule drugs into their protein targets (Deeks et al. PLoS One 2020, 15, e0228461). Here, we apply this approach to both an Mpro inhibitor and an oligopeptide substrate, using experimentally determined crystal structures. For the oligopeptide, we test against a crystallographic structure of the original SARS Mpro. Docking with iMD-VR gives models in agreement with experimentally observed (crystal) structures. The docked structures are also tested in MD simulations and found to be stable. Different protocols for iMD-VR docking are explored, e.g., with and without restraints on protein backbone, and we provide recommendations for its use. We find that it is important for the user to focus on forming binding interactions, such as hydrogen bonds, and not to rely on using simple metrics (such as RMSD), in order to create realistic, stable complexes. We also test the use of apo (uncomplexed) crystal structures for docking and find that they can give good results. This is because of the flexibility and dynamic response allowed by the physically rigorous, atomically detailed simulation approach of iMD-VR. We make our models (and interactive simulations) freely available. The software framework that we use, Narupa, is open source, and uses commodity VR hardware, so these tools are readily accessible to the wider research community working on Mpro (and other COVID-19 targets). These should be widely useful in drug development, in education applications, e.g., on viral enzyme structure and function, and in scientific communication more generally.

Novel cyclohexanone compound as a potential ligand against SARS-CoV-2 main-protease.

No commercially available drug candidate has yet been devised which is unique to and not repurposed against SARS-CoV-2 and has high efficacy or safe toxicity profile or both. Taking curcumin as a reference compound, we identified a new commercially available cyclohexanone compound, ZINC07333416 with binding energy (-8.72 kcal/mol) better than that of popularly devised anti-Covid-19 drugs like viral protease inhibitor Lopinavir, nucleoside analogue Remdesivir and the repurposed drug hydroxychloroquine when targeted to the active-site of SARS-CoV-2 Main protease (Mpro) through docking studies. The ligand ZINC07333416 exhibits crucial interactions with major active site residues of SARS-CoV-2 Mpro viz. Cys145 and His41 involving in the protease activity; as well as GLU-166 and ASN-142 which plays the pivotal role in the protein-dimerization. The protein-ligand stable interaction was further confirmed with molecular dynamics simulation (MDS) studies. Based on virtual assessment, ZINC07333416 also have significant values in terms of medicinal chemistry, pharmacokinetics, synthetic accessibility and anti-viral activity that encourage its experimental applications against COVID-19.

Comparative Antiviral Efficacy of Viral Protease Inhibitors against the Novel SARS-CoV-2 In Vitro.

The recent outbreak of novel coronavirus pneumonia (COVID-19) caused by a new coronavirus has posed a great threat to public health. Identifying safe and effective antivirals is of urgent demand to cure the huge number of patients. Virus-encoded proteases are considered potential drug targets. The human immunodeficiency virus protease inhibitors (lopinavir/ritonavir) has been recommended in the global Solidarity Trial in March launched by World Health Organization. However, there is currently no experimental evidence to support or against its clinical use. We evaluated the antiviral efficacy of lopinavir/ritonavir along with other two viral protease inhibitors in vitro, and discussed the possible inhibitory mechanism in silico. The in vitro to in vivo extrapolation was carried out to assess whether lopinavir/ritonavir could be effective in clinical. Among the four tested compounds, lopinavir showed the best inhibitory effect against the novel coronavirus infection. However, further in vitro to in vivo extrapolation of pharmacokinetics suggested that lopinavir/ritonavir could not reach effective concentration under standard dosing regimen [marketed as Kaletra®, contained lopinavir/ritonavir (200 mg/50 mg) tablets, recommended dosage is 400 mg/10 mg (2 tablets) twice daily]. This research concluded that lopinavir/ritonavir should be stopped for clinical use due to the huge gap between in vitro IC50 and free plasma concentration. Nevertheless, the structure-activity relationship analysis of the four inhibitors provided further information for de novel design of future viral protease inhibitors of SARS-CoV-2.

Structure-Based Virtual Screening to Discover Potential Lead Molecules for the SARS-CoV-2 Main Protease.

The COVID-19 disease is caused by a new strain of the coronavirus family (SARS-CoV-2), and it has affected at present millions of people all over the world. The indispensable role of the main protease (Mpro) in viral replication and gene expression makes this enzyme an attractive drug target. Therefore, inhibition of SARS-CoV-2 Mpro as a proposition to halt virus ingression is being pursued by scientists globally. Here we carried out a study with two objectives: the first being to perform comparative protein sequence and 3D structural analysis to understand the effect of 12 point mutations on the active site. Among these, two mutations, viz., Ser46 and Phe134, were found to cause a significant change at the active sites of SARS-CoV-2. The Ser46 mutation present at the entrance of the S5 subpocket of SARS-CoV-2 increases the contribution of other two hydrophilic residues, while the Phe134 mutation, present in the catalytic cysteine loop, can cause an increase in catalytic efficiency of Mpro by facilitating fast proton transfer from the Cys145 to His41 residue. It was observed that active site remained conserved among Mpro of both SARS-CoVs, except at the entrance of the S5 subpocket, suggesting sustenance of substrate specificity. The second objective was to screen the inhibitory effects of three different data sets (natural products, coronaviruses main protease inhibitors, and FDA-approved drugs) using a structure-based virtual screening approach. A total of 73 hits had a combo score >2.0. Eight different structural scaffold classes were identified, such as one/two tetrahydropyran ring(s), dipeptide/tripeptide/oligopeptide, large (approximately 20 atoms) cyclic peptide, and miscellaneous. The screened hits showed key interactions with subpockets of the active site. Further, molecular dynamics studies of selected screened compounds confirmed their perfect fitting into the subpockets of the active site. This study suggests promising structures that can fit into the SARS-CoV-2 Mpro active site and also offers direction for further lead optimization and rational drug design.