Identification of natural compounds as potent inhibitors of SARS-CoV-2 main protease using combined docking and molecular dynamics simulations

Coronavirus disease 2019 (COVID-19) has emerged from China and globally affected the entire population through the human-to-human transmission of a newly emerged virus called severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The genome of SARS-CoV-2 encodes several proteins that are essential for multiplication and pathogenesis. The main protease (M^pro or 3CL^pro) of SARS-CoV-2 plays a central role in its pathogenesis and thus is considered as an attractive drug target for the drug design and development of small-molecule inhibitors. We have employed an extensive structure-based high-throughput virtual screening to discover potential natural compounds from the ZINC database which could inhibit the M^pro of SARS-CoV-2. Initially, the hits were selected on the basis of their physicochemical and drug-like properties. Subsequently, the PAINS filter, estimation of binding affinities using molecular docking, and interaction analyses were performed to find safe and potential inhibitors of SARS-CoV-2 M^pro. We have identified ZINC02123811 (1-(3-(2,5,9-trimethyl-7-oxo-3-phenyl-7H-furo[3,2-g]chromen-6-yl)propanoyl)piperidine-4-carboxamide), a natural compound bearing appreciable affinity, efficiency, and specificity towards the binding pocket of SARS-CoV-2 M^pro. The identified compound showed a set of drug-like properties and preferentially binds to the active site of SARS-CoV-2 M^pro. All-atom molecular dynamics (MD) simulations were performed to evaluate the conformational dynamics, stability and interaction mechanism of M^pro with ZINC02123811. MD simulation results indicated that M^pro with ZINC02123811 forms a stable complex throughout the trajectory of 100 ns. These findings suggest that ZINC02123811 may be further exploited as a promising scaffold for the development of potential inhibitors of SARS-CoV-2 M^pro to address COVID-19.     

An in silico approach unveils the potential of antiviral compounds in preclinical and clinical trials as SARS-CoV-2 omicron inhibitors

The increased transmissibility and highly infectious nature of the new variant of concern (VOC) that is severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron and lack of effective therapy need the rapid discovery of therapeutic antivirals against it. The present investigation aimed to identify antiviral compounds that would be effective against SARS-CoV-2 Omicron. In this study, molecular docking experiments were carried out using the recently reported experimental structure of omicron spike protein in complex with human angiotensin-converting enzyme 2 (ACE2) and various antivirals in preclinical and clinical trial studies. Out of 36 tested compounds, Abemaciclib, Dasatinib and Spiperone are the three top-ranked molecules which scored binding energies of ?10.08 kcal/mol, ?10.06 kcal/mol and ?9.54 kcal/mol respectively. Phe338, Asp339, and Asp364 are crucial omicron receptor residues involved in hydrogen bond interactions, while other residues were mostly involved in hydrophobic interactions with the lead molecules. The identified lead compounds also scored well in terms of drug-likeness. Molecular dynamics (MD) simulation, essential dynamics (ED) and entropic analysis indicate the ability of these molecules to modulate the activity of omicron spike protein. Therefore, Abemaciclib, Dasatinib and Spiperone are likely to be viable drug-candidate molecules that can block the interaction between the omicron spike protein and the host cellular receptor ACE2. Though our findings are compelling, more research into these molecules is needed before they can be employed as drugs to treat SARS-CoV-2 omicron infections.     

In silico structure-based discovery of a SARS-CoV-2 main protease inhibitor

The Coronavirus Disease 2019 (COVID-19) pandemic caused by the novel lineage B betacoroanvirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in significant mortality, morbidity, and socioeconomic disruptions worldwide. Effective antivirals are urgently needed for COVID-19. The main protease (M^pro) of SARS-CoV-2 is an attractive antiviral target because of its essential role in the cleavage of the viral polypeptide. In this study, we performed an in silico structure-based screening of a large chemical library to identify potential SARS-CoV-2 M^pro inhibitors. Among 8,820 compounds in the library, our screening identified trichostatin A, a histone deacetylase inhibitor and an antifungal compound, as an inhibitor of SARS-CoV-2 M^pro activity and replication. The half maximal effective concentration of trichostatin A against SARS-CoV-2 replication was 1.5 to 2.7µM, which was markedly below its 50% effective cytotoxic concentration (75.7µM) and peak serum concentration (132µM). Further drug compound optimization to develop more stable analogues with longer half-lives should be performed. This structure-based drug discovery platform should facilitate the identification of additional enzyme inhibitors of SARS-CoV-2.     

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 (M^pro) 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 M^pro 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 M^pro specifically for α, β, γ (K90R), λ (G15S), and ο (P132H). In vitro biochemical enzymatic assay characterization of the enzyme kinetics of the mutant M^pros demonstrates that they are catalytically comparable to wildtype. We found that nirmatrelvir has similar potency against each mutant M^pro 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 β M^pro 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.     

Studies of Cytotoxicity Effects,SARS-CoV-2 Main Protease Inhibition,and in Silico Interactions of Synthetic Chalcones

SARS-CoV-2 main protease (M^pro) plays an essential role in proteolysis cleavage that promotes coronavirus replication. Thus, attenuating the activity of this enzyme represents a strategy to develop antiviral agents. We report inhibitory effects against M^pro of 40 synthetic chalcones, and cytotoxicity activities, hemolysis, and in silico interactions of active compounds. Seven of them bearing a (E)-3-(furan-2-yl)-1-arylprop-2-en-1-one skeleton ( 10 , 28 , and 35–39 ) showed enzyme inhibition with IC₅₀ ranging from 13.76 and 36.13 μM. Except for 35 and 36 , other active compounds were not cytotoxic up to 150 μM against THP-1 and Vero cell lines. Compounds 10 , and 35–39 showed no hemolysis while 28 was weakly hemotoxic at 150 μM. Moreover, molecular docking showed interactions between compound 10 and M^pro (PDBID 5RG2 and 5RG3) with proximity to cys145 and His41, suggesting a covalent binding. Products of the reaction between chalcones and cyclohexanethiol indicated that this binding could be a Michael addition type.     

Scutellaria baicalensis extract and baicalein inhibit replication of SARS-CoV-2 and its 3C-like protease in vitro

COVID-19 has become a global pandemic and there is an urgent call for developing drugs against the virus (SARS-CoV-2). The 3C-like protease (3CL^pro) of SARS-CoV-2 is a preferred target for broad spectrum anti-coronavirus drug discovery. We studied the anti-SARS-CoV-2 activity of S. baicalensis and its ingredients. We found that the ethanol extract of S. baicalensis and its major component, baicalein, inhibit SARS-CoV-2 3CL^pro activity in vitro with IC₅₀’s of 8.52 µg/ml and 0.39 µM, respectively. Both of them inhibit the replication of SARS-CoV-2 in Vero cells with EC₅₀’s of 0.74 µg/ml and 2.9 µM, respectively. While baicalein is mainly active at the viral post-entry stage, the ethanol extract also inhibits viral entry. We further identified four baicalein analogues from other herbs that inhibit SARS-CoV-2 3CL^pro activity at µM concentration. All the active compounds and the S. baicalensis extract also inhibit the SARS-CoV 3CL^pro, demonstrating their potential as broad-spectrum anti-coronavirus drugs.     

Structural Basis of Main Proteases of Coronavirus Bound to Drug Candidate PF-07304814

New variants of the severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2) emerged and spread rapidly all over the world, which strongly supports the need for pharmacological options to complement vaccine strategies. Main protease (M^pro or 3CL^pro) is a critical enzyme in the life cycle of SARS-CoV-2 and appears to be highly conserved among different genera of coronaviruses, making it an ideal target for the development of drugs with broad-spectrum property. PF-07304814 developed by Pfizer is an intravenously administered inhibitor targeting SARS-CoV-2 M^pro. Here we showed that PF-07304814 displays broad-spectrum inhibitory activity against M^pros from multiple coronaviruses. Crystal structures of M^pros of SARS-CoV-2, SARS-CoV, MERS-CoV, and HCoV-NL63 bound to the inhibitor PF-07304814 revealed a conserved ligand-binding site, providing new insights into the mechanism of inhibition of viral replication. A detailed analysis of these crystal structures complemented by comprehensive comparison defined the key structural determinants essential for inhibition and illustrated the binding mode of action of M^pros from different coronaviruses. In view of the importance of M^pro for the medications of SARS-CoV-2 infection, insights derived from the present study should accelerate the design of pan-coronaviral main protease inhibitors that are safer and more effective.     

Caffeic acid derivatives (CAFDs) as inhibitors of SARS-CoV-2: CAFDs-based functional foods as a potential alternative approach to combat COVID-19

BackgroundSARS-CoV-2, an emerging strain of coronavirus, has affected millions of people from all the continents of world and received worldwide attention. This emerging health crisis calls for the urgent development of specific therapeutics against COVID-19 to potentially reduce the burden of this emerging pandemic.PurposeThis study aims to evaluate the anti-viral efficacy of natural bioactive entities against COVID-19 via molecular docking and molecular dynamics simulation.MethodsA library of 27 caffeic-acid derivatives was screened against 5 proteins of SARS-CoV-2 by using Molegro Virtual Docker 7 to obtain the binding energies and interactions between compounds and SARS-CoV-2 proteins. ADME properties and toxicity profiles were investigated via www.swissadme.ch web tools and Toxtree respectively. Molecular dynamics simulation was performed to determine the stability of the lead-protein interactions.ResultsOur obtained results has uncovered khainaoside C, 6-O-Caffeoylarbutin, khainaoside B, khainaoside C and vitexfolin A as potent modulators of COVID-19 possessing more binding energies than nelfinavir against COVID-19 M^pro, Nsp15, SARS-CoV-2 spike S2 subunit, spike open state and closed state structure respectively. While Calceolarioside B was identified as pan inhibitor, showing strong molecular interactions with all proteins except SARS-CoV-2 spike glycoprotein closed state. The results are supported by 20 ns molecular dynamics simulations of the best complexes.ConclusionThis study will hopefully pave a way for development of phytonutrients-based antiviral therapeutic for treatment or prevention of COVID-19 and further studies are recommended to evaluate the antiviral effects of these phytochemicals against SARS-CoV-2 in in vitro and in vivo models.     

A COVID moonshot: assessment of ligand binding to the SARS-CoV-2 main protease by saturation transfer difference NMR spectroscopy

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the etiological cause of the coronavirus disease 2019, for which no effective antiviral therapeutics are available. The SARS-CoV-2 main protease (M^pro) is essential for viral replication and constitutes a promising therapeutic target. Many efforts aimed at deriving effective M^pro inhibitors are currently underway, including an international open-science discovery project, codenamed COVID Moonshot. As part of COVID Moonshot, we used saturation transfer difference nuclear magnetic resonance (STD-NMR) spectroscopy to assess the binding of putative M^pro ligands to the viral protease, including molecules identified by crystallographic fragment screening and novel compounds designed as M^pro inhibitors. In this manner, we aimed to complement enzymatic activity assays of M^pro performed by other groups with information on ligand affinity. We have made the M^pro STD-NMR data publicly available. Here, we provide detailed information on the NMR protocols used and challenges faced, thereby placing these data into context. Our goal is to assist the interpretation of M^pro STD-NMR data, thereby accelerating ongoing drug design efforts.

     

Improved SARS-CoV-2 main protease high-throughput screening assay using a 5-carboxyfluorescein substrate

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) as a global threat to human health has highlighted the need for the development of novel therapies targeting current and emerging coronaviruses with pandemic potential. The coronavirus main protease (M^pro, also called 3CL^pro) is a validated drug target against coronaviruses and has been heavily studied since the emergence of SARS-CoV-2 in late 2019. Here, we report the biophysical and enzymatic characterization of native M^pro, then characterize the steady-state kinetics of several commonly used FRET substrates, fluorogenic substrates, and six of the 11 reported SARS-CoV-2 polyprotein cleavage sequences. We then assessed the suitability of these substrates for high-throughput screening. Guided by our assessment of these substrates, we developed an improved 5-carboxyfluorescein-based FRET substrate, which is better suited for high-throughput screening and is less susceptible to interference and false positives than existing substrates. This study provides a useful framework for the design of coronavirus M^pro enzyme assays to facilitate the discovery and development of therapies targeting M^pro.     

The Molecular Mechanism of Domain Swapping of the C-Terminal Domain of the SARS-Coronavirus Main Protease

In three-dimensional domain swapping, two protein monomers exchange a part of their structures to form an intertwined homodimer, whose subunits resemble the monomer. Several viral proteins domain swap to increase their structural complexity or functional avidity. The main protease (M^pro) of the severe acute respiratory syndrome (SARS) coronavirus proteolyzes viral polyproteins and has been a target for anti-SARS drug design. Domain swapping in the α-helical C-terminal domain of M^pro (M^proC) locks M^pro into a hyperactive octameric form that is hypothesized to promote the early stages of viral replication. However, in the absence of a complete molecular understanding of the mechanism of domain swapping, investigations into the biological relevance of this octameric M^pro have stalled. Isolated M^proC can exist as a monomer or a domain-swapped dimer. Here, we investigate the mechanism of domain swapping of M^proC using coarse-grained structure-based models and molecular dynamics simulations. Our simulations recapitulate several experimental features of M^proC folding. Further, we find that a contact between a tryptophan in the M^proC domain-swapping hinge and an arginine elsewhere forms early during folding, modulates the folding route, and promotes domain swapping to the native structure. An examination of the sequence and the structure of the tryptophan containing hinge loop shows that it has a propensity to form multiple secondary structures and contacts, indicating that it could be stabilized into either the monomer- or dimer-promoting conformations by mutations or ligand binding. Finally, because all residues in the tryptophan loop are identical in SARS-CoV and SARS-CoV-2, mutations that modulate domain swapping may provide insights into the role of octameric M^pro in the early-stage viral replication of both viruses.     

N-Terminal Finger Stabilizes the S1 Pocket for the Reversible Feline Drug GC376 in the SARS-CoV-2 Mpro Dimer

The main protease (M^pro, also known as 3CL protease) of SARS-CoV-2 is a high priority drug target in the development of antivirals to combat COVID-19 infections. A feline coronavirus antiviral drug, GC376, has been shown to be effective in inhibiting the SARS-CoV-2 main protease and live virus growth. As this drug moves into clinical trials, further characterization of GC376 with the main protease of coronaviruses is required to gain insight into the drug’s properties, such as reversibility and broad specificity. Reversibility is an important factor for therapeutic proteolytic inhibitors to prevent toxicity due to off-target effects. Here we demonstrate that GC376 has nanomolar K_i values with the M^pro from both SARS-CoV-2 and SARS-CoV strains. Restoring enzymatic activity after inhibition by GC376 demonstrates reversible binding with both proteases. In addition, the stability and thermodynamic parameters of both proteases were studied to shed light on physical chemical properties of these viral enzymes, revealing higher stability for SARS-CoV-2 M^pro. The comparison of a new X-ray crystal structure of M^pro from SARS-CoV complexed with GC376 reveals similar molecular mechanism of inhibition compared to SARS-CoV-2 M^pro, and gives insight into the broad specificity properties of this drug. In both structures, we observe domain swapping of the N-termini in the dimer of the M^pro, which facilitates coordination of the drug’s P1 position. These results validate that GC376 is a drug with an off-rate suitable for clinical trials.     

Molecular identification and biological control of Ralstonia solanacearum from wilt of papaya by natural compounds and Bacillus subtilis: An integrated experimental and computational study

Ralstonia solanacearum is a harmful pathogen that causes severe wilt disease in several vegetables. In the present study, we identified R. solanacearum from wilt of papaya by 16S rRNA PCR amplification. Virulence ability of R. solanacearum was determined by amplification of approximately 1500 bp clear band of hrpB gene. Further, in-vitro seed germination assay showed that R. solanacearum reduced the germination rate up to 26.21%, 34% and 33.63% of cucumber, bottle guard and pumpkin seeds, respectively whereas shoot and root growth were also significantly decreased. Moreover, growth inhibition of R. solanacearum was recorded using antibacterial compound from medicinal plant and antagonistic B. subtilis. Petroleum ether root extract of Rauvolfia serpentina showed highest 22 ± 0.04 mm diameter of zone of inhibition where methanolic extract of Cymbopogon citratus and ethanolic extract of Lantana camara exhibited 20 ± 0.06 mm and 20 ± 0.01 mm zone of inhibition against R. solanacearum, respectively. In addition, bioactive compounds of B. subtilis inhibited R. solanacearum growth by generating 17 ± 0.09 mm zone of inhibition. To unveil the inhibition mechanism, we adopted chemical-protein interaction network and molecular docking approaches where we found that, rutin from C. citratus interacts with citrate (Si)-synthase and dihydrolipoyl dehydrogenase of R. solanacearum with binding affinity of −9.7 kcal/mol and −9.5 kcal/mol while quercetin from B. subtillis interacts with the essential protein F0F1 ATP synthase subunit alpha of the R. solancearum with binding affinity of −6.9 kcal/mol and inhibit the growth of R. solanacearum. Our study will give shed light on the development of eco-friendly biological control of wilt disease of papaya.     

A Crystallographic Snapshot of SARS-CoV-2 Main Protease Maturation Process

SARS-CoV-2 is the causative agent of COVID-19. The dimeric form of the viral M^pro is responsible for the cleavage of the viral polyprotein in 11 sites, including its own N and C-terminus. The lack of structural information for intermediary forms of M^pro is a setback for the understanding its self-maturation process. Herein, we used X-ray crystallography combined with biochemical data to characterize multiple forms of SARS-CoV-2 M^pro. For the immature form, we show that extra N-terminal residues caused conformational changes in the positioning of domain-three over the active site, hampering the dimerization and diminishing its activity. We propose that this form preludes the cis and trans-cleavage of N-terminal residues. Using fragment screening, we probe new cavities in this form which can be used to guide therapeutic development. Furthermore, we characterized a serine site-directed mutant of the M^pro bound to its endogenous N and C-terminal residues during dimeric association stage of the maturation process. We suggest this form is a transitional state during the C-terminal trans-cleavage. This data sheds light in the structural modifications of the SARS-CoV-2 main protease during its self-maturation process.     

Extending the Record of Terpenes from Soft Coral Sarcophyton mililatensis

In the present study, a new polyoxygenated cembranoid named sarcomililatol H ( 1 ) as well as six known terpenes 2 – 7 with different skeletons were isolated from South China Sea soft coral Sarcophyton mililatensis. Based on the comprehensive analyses of 1D and 2D NMR spectroscopic data, the structure of the new compound 1 was established. This new cembranoid was characterized by the presence of the rarely encountered tetrahydropyran ring with the ether linkage across C-2 and C-12. By applying the time-dependent density functional theory electronic circular dichroism (TDDFT ECD) approach, the absolute configuration of sarcomililatol H ( 1 ) was determined. All of the isolates were subjected to the anti-inflammatory and anti-tumor bioassays. However, none of them was active in these evaluations. Additionally, the preliminary virtual screening of inhibitory against SARS-CoV-2 by molecular docking showed that diterpene 1 could be regarded as a SARS-CoV-2 main protease (M^pro) inhibitor (binding energy: −7.63 kcal/mol). The discovery of these terpenes has expanded the chemical diversity and complexity of terpenes from the species S. mililatensis.     

Investigation of changes in protein stability and substrate affinity of 3CL-protease of SARS-CoV-2 caused by mutations

3CL^pro of SARS-CoV-2 is one of the enzymes required for the replication process of the virus responsible for the COVID-19 pandemic. In this study, changes in protein stability and substrate affinity caused by mutations were investigated to stir the development of potent inhibitors. Sequence data of samples were obtained from the NCBI Virus database. Mutation analyses were performed with RDP4 and MegaX. 3CL^pro tertiary models were created using Robetta. Molecular docking for peptidomimetic substrate and inhibitor ligand was done with Autodock v4.2 and Haddock v2.4. Protein stability analysis was performed using mCSM stability and DynaMut2. Twenty-four missense mutations in 3CL^pro were identified in this study. Changes in the 3CL^pro structure induced by the mutations Met49Thr, Leu167Ser, and Val202Ala resulted in significant levels of instability (-2.029,-2.612,-2.177 kcal.mol^-1, respectively). The lowest interaction energy for substrate was -58.7 kcal.mol^-1 and -62.6 kcal.mol^-1 in wild-type and mutant, respectively. The lowest docking energy for ligand was -6.19 and -9.52 kcal.mol^-1 for wild-type and mutant, respectively. This study reports for the first time that mutations cause increased substrate affinity of 3CL^pro from SARS-CoV-2. This research provides important data for the development of potent peptidomimetic inhibitors for the treatment of COVID-19. Keywords: 3CL-protease, mutation analysis, protein stability, SARS-CoV-2 genome, substrate affinity     

The structural basis of conserved residue variant effect on enzyme activity of UGT2B15

UDP-glucuronosyltransferase 2B15 (UGT2B15) is a crucial phase II drug-metabolizing enzyme, which glucuronidates various compounds, including clinical drugs and hormones. Mutants might affect glucuronidation, leading to a disruption of drug metabolism in vivo and decrease of therapeutic effect. Here, we mainly analyzed two representative mutants, H401P and L446S, on UGT2B15 activity using glucuronidation assays, molecular dynamic (MD) simulation and X-ray diffraction methods. The enzyme activity of L446S obviously increased six-fold than the wild type, although the enzyme activities of P191L, T374A, and H401P were lost apparently. Furthermore, we used MD simulations to calculate the energy change in the catalytic process of H401P and L446S, and the results indicated the free binding energies of H401P mutant to oxazepam and UDPGA were ?30.98 ± 1.00 kcal/mol and ?36.42 ± 1.04 kcal/mol, respectively, increased obviously compared to wild type, suggesting the mutation on position 401 had a crucial effect on the catalysis. Moreover, the three-dimensional structure of UGT2B15 C-terminal domain L446S was determined through protein crystallography and X-ray diffraction technology and the results suggested that one more hydrogen bonding between S446 and K410 was formed in the S446 crystal structure, compared to the wild type. Isothermal titration calorimetry assay further revealed the K_d values of C-terminal domain of UGT2B15 harbored L446S towards the cofactor UDPGA was similar to the value of wild type. Above all, our results pointed out that H401P and L446S affected the enzyme activity by different mechanism. Our work provided a helpful mechanism for variance explained in the UGTs catalyzation process.