Repurposing of phytomedicine-derived bioactive compounds with promising anti-SARS-CoV-2 potential: Molecular docking,MD simulation and drug-likeness/ADMET studies

In view of the potential of traditional plant-based remedies (or phytomedicines) in the management of COVID-19, the present investigation was aimed at finding novel anti-SARS-CoV-2 molecules by in silico screening of bioactive phytochemicals (database) using computational methods and drug repurposing approach. A total of 160 compounds belonging to various phytochemical classes (flavonoids, limonoids, saponins, triterpenoids, steroids etc.) were selected (as initial hits) and screened against three specific therapeutic targets (Mpro/3CLpro, PLpro and RdRp) of SARS-CoV-2 by docking, molecular dynamics simulation and drug-likeness/ADMET studies. From our studies, six phytochemicals were identified as notable ant-SARS-CoV-2 agents (best hit molecules) with promising inhibitory effects effective against protease (Mpro and PLpro) and polymerase (RdRp) enzymes. These compounds are namely, ginsenoside Rg2, saikosaponin A, somniferine, betulinic acid, soyasapogenol C and azadirachtin A. On the basis of binding modes and dynamics studies of protein–ligand intercations, ginsenoside Rg2, saikosaponin A, somniferine were found to be the most potent (in silico) inhibitors potentially active against Mpro, PLpro and RdRp, respectively. The present investigation can be directed towards further experimental studies in order to confirm the anti-SARS-CoV-2 efficacy along with toxicities of identified phytomolecules.     

Potential targets of severe acute respiratory syndrome coronavirus 2 of clinical drug fluvoxamine: Docking and molecular dynamics studies to elucidate viral action

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has continued evolving for survival and adaptation by mutating itself into different variants of concern, including omicron. Several studies and clinical trials found fluvoxamine, an Food and Drug Administration-approved antidepressant drug, to be effective at preventing mild coronavirus disease 2019 (COVID-19) from progressing to severe diseases. However, the mechanism of fluvoxamine's direct antiviral action against COVID-19 is still unknown. Fluvoxamine was docked with 11 SARS-CoV-2 targets and subjected to stability, conformational changes, and binding free energy analyses to explore its mode of action. Of the targets, nonstructural protein 14 (NSP14), main protease (Mpro), and papain-like protease (PLpro) had the best docking scores with fluvoxamine. Consistent with the docking results, it was confirmed by molecular dynamics simulations that the NSP14 N7-MTase ((N7-guanine)-methyltransferase)–fluvoxamine, Mpro–fluvoxamine, and PLpro–fluvoxamine complexes are stable, with the lowest binding free energies of −105.1, −82.7, and − 38.5 kJ/mol, respectively. A number of hotspot residues involved in the interaction were also identified. These include Glu166, Asp187, His41, and Cys145 in Mpro, Gly163 and Arg166 in PLpro, and Glu302, Gly333, and Phe426 in NSP14, which could aid in the development of better antivirals against SARS-CoV-2.     

High-throughput screening identifies established drugs as SARS-CoV-2 PLpro inhibitors

A new coronavirus(SARS-CoV-2)has been identified as the etiologic agent for the COVID-19 outbreak.Currently,effective treatment options remain very limited for this disease;therefore,there is an urgent need to identify new anti-COVID-19 agents.In this study,we screened over 6,000 compounds that included approved drugs,drug candidates in clinical trials,and pharmacologically active compounds to identify leads that target the SARS-CoV-2 papain-like protease(PLpro).Together with main protease(Mpro),PLpro is responsible for processing the viral replicase polyprotein into functional units.There-fore,it is an attractive target for antiviral drug develop-ment.Here we discovered four compounds,YM155,cryptotanshinone,tanshinone I and GRL0617 that inhibit SARS-CoV-2 PLpro with IC50 values ranging from 1.39 to 5.63 pmol/L.These compounds also exhibit strong antiviral activities in cell-based assays.YM155,an anti-cancer drug candidate in clinical trials,has the most potent antiviral activity with an EC50 value of 170 nmol/L.In addition,we have determined the crystal structures of this enzyme and its complex with YM155,revealing a unique binding mode.YM155 simultaneously targets three"hot"spots on PLpro,including the substrate-binding pocket,the interferon stimulating gene product 15(ISG15)binding site and zinc finger motif.Our results demonstrate the efficacy of this screening and repur-posing strategy,which has led to the discovery of new drug leads with clinical potential for COVID-19 treatments.     

Multi-targeted molecular docking,pharmacokinetics, and drug-likeness evaluation of coumarin based compounds targeting proteins involved in development of COVID-19

COVID-19 is a progressing pandemic of coronavirus disease-2019, which had drowned the whole world in a deep sorrow sea. Uncountable deaths were extending the list of deaths every single day. The present research was aimed to study the multi-target interaction of coumarins against COVID-19 using molecular docking analysis. The structure of coumarin compounds was checked for ADME and Lipinski rule of five by using SwissADME, an online tool. SARS-CoV-2 proteins such as RdRp, PLpro, Mpro and spike protein were collected from the Protein Data Bank. The molecular docking study was performed in the PyRx tool, and the molecular interactions were visualised by Discovery Studio Visualizer. All the coumarin compounds used in the study were obeyed Lipinski’s rule of 5 without any violations. All the three designed derivatives of phenprocoumon, hymecromone, and psoralen were showed high binding affinity and prominent interactions with the drug target. The presence of –OH groups in the compound, His41, a catalytic dyad in Mpro, number of and the distance of hydrogen bond interactions with SARS-CoV-2 targets was accountable for the high binding attractions. The modified drug structures possess better binding efficacy towards at least three targets compared to their parent compounds. Further, molecular dynamic studies can be suggested to find the ligand–protein complex stability. The present study outcome reveals that the designed coumarins can be synthesised and examined as a potent inhibitory drug of SARS-CoV-2.     

新型冠状病毒主蛋白酶抑制剂的筛选方法研究进展

新型冠状病毒肺炎(coronavirus disease 2019,COVID-19)席卷全球,具有较高的传染性和死亡率,但目前尚缺乏安全有效的COVID-19疫苗与治疗药物.新型冠状病毒主蛋白酶(main protease,Mpro)的进化高度保守,在调控新冠病毒RNA复制中具有重要的生物学功能,已成为新型广谱抗冠状...     

Targeting papain-like protease for broad-spectrum coronavirus inhibition

The emergence of SARS-CoV-2 variants of concern and repeated outbreaks of coronavirus epidemics in the past two decades emphasize the need for next-generation pan-coronaviral therapeutics. Drugging the multi-functional papain-like protease (PLpro) domain of the viral nsp3 holds promise. However, none of the known coronavirus PLpro inhibitors has been shown to be in vivo active. Herein, we screened a structurally diverse library of 50,080 compounds for potential coronavirus PLpro inhibitors and identified a noncovalent lead inhibitor F0213 that has broad-spectrum anti-coronaviral activity, including against the Sarbecoviruses (SARS-CoV-1 and SARS-CoV-2), Merbecovirus (MERS-CoV), as well as the Alphacoronavirus (hCoV-229E and hCoV-OC43). Importantly, F0213 confers protection in both SARS-CoV-2-infected hamsters and MERS-CoV-infected human DPP4-knockin mice. F0213 possesses a dual therapeutic functionality that suppresses coronavirus replication via blocking viral polyprotein cleavage, as well as promoting antiviral immunity by antagonizing the PLpro deubiquitinase activity. Despite the significant difference of substrate recognition, mode of inhibition studies suggest that F0213 is a competitive inhibitor against SARS2-PLpro via binding with the 157K amino acid residue, whereas an allosteric inhibitor of MERS-PLpro interacting with its 271E position. Our proof-of-concept findings demonstrated that PLpro is a valid target for the development of broad-spectrum anti-coronavirus agents. The orally administered F0213 may serve as a promising lead compound for combating the ongoing COVID-19 pandemic and future coronavirus outbreaks.Supplementary InformationThe online version contains supplementary material available at 10.1007/s13238-022-00909-3.     

新型冠状病毒PLpro蛋白酶结构与功能的生物信息学分析

2019年12月,由新型冠状病毒(SARS-CoV-2)引起的新型冠状病毒肺炎(COVID-19)在中国武汉暴发。SARS-CoV-2的基因组编码2种病毒蛋白酶,即木瓜样蛋白酶(Papain-like protease,PLpro)和3C样蛋白酶(3C-like protease)。其中PLpro是SARS-CoV-2复制酶复合体(RC)形成的重要调节蛋白分子,对于病毒基因组转录和复制至关重要。因此,将SARS-CoV-2 PLpro作为药物的靶点对COVID-19的治疗具有积极意义。本研究应用生物信息学工具分析新型冠状病毒的木瓜样蛋白酶的结构和功能,首先利用BLAST和BioEdit获取SARS-CoV-2 PLpro蛋白酶(SC2-PLpro)及其同源蛋白的氨基酸序列,并利用BLAST和MEGA 6.0进行同源性分析。之后,利用ProParam和Proscale分别对SC2-PLpro蛋白酶的理化性质、亲水性和疏水性进行分析。然后,通过SOMPA、ScanProsite和InterPro分别预测SC2-PLpro蛋白酶的二级结构和功能区域,进一步利用SignalP 4.0和TMHMM对SC2-PLpro蛋白酶的信号肽和跨膜区进行分析。最后,通过SWISS-MODEL对SARS-CoV-2 PLpro蛋白酶进行三级结构同源建模。结果显示,对SARS-CoV-2 PLpro蛋白酶与已报道的PLpro蛋白酶进行多序列比对后,发现SARS-CoV-2 nsp3的746~1063段氨基酸与多种冠状病毒PLpro蛋白酶氨基酸序列高度相似。同时,同源性分析发现SARS-CoV-2与蝙蝠冠状病毒的PLpro蛋白酶具有同源性,其中与QHR63299、AVP78030相似性最高。对SC2-PLpro进行理化性质预测结果显示,其由318个氨基酸所组成,为稳定亲水性蛋白。二级结构预测结果显示SC2-PLpro主要含有α-螺旋、延伸链、β-转角、无规卷曲,四种结构贯穿整条氨基酸链。进一步进行功能分析,发现其具有完整的催化三联体、锌结合域、泛素样N末端结构域,故推测该蛋白具有去泛素化的功能。然后,信号肽假说和跨膜结构域分析结果表明,SC2-PLpro既不是分泌蛋白,也不属于跨膜蛋白。本研究提示,生物信息学分析SC2-PLpro为稳定性亲水蛋白,属于非跨膜蛋白,比较保守,具有去泛素化的功能,利用此功能可以进一步规避宿主的固有免疫反应。通过制备PLpro蛋白酶小分子抑制剂,可能有助于治疗新型冠状病毒肺炎。     

In silico investigation of ACE2 and the main protease of SARS-CoV-2 with phytochemicals from Myristica fragrans (Houtt.) for the discovery of a novel COVID-19 drug

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes coronavirus disease 2019 (COVID-19), is a new coronavirus strain that was first reported in December 2019 in Wuhan, China. A specific treatment for COVID-19 has yet to be identified. Potential therapeutic targets include SARS-CoV-2 main protease (Mpro) and the SARS-CoV-2 spike-ACE2 interaction. Molecular docking, molecular dynamics (MD), solvent screening for the extraction of the specified compounds, and prediction of the drug properties of certain molecules were the methods used in this study to investigate compounds from the medicinal plant Myristica fragrans, which is one of twelve herbs in Prasachandaeng remedy (PSD). ArgusLab, AutoDock Vina, and AutoDock were used to perform docking tasks. The examined ligands were compared with panduratin A as a standard (Kanjanasirirat et al., 2020), which is a promising medicinal plant molecule for the treatment of COVID-19. Molecular docking revealed that malabaricones B and C and licarins A, B and C bound to SARS-CoV-2/ACE2 and SARS-CoV-2 Mpro with low binding energies compared to that of the standard ligand. Furthermore, appropriate solvent usage is important. Acetone was selected by COSMOquick software for compound extraction in this investigation because it can extract large amounts of all five of the abovementioned M. fragrans compounds. Furthermore, the drug-like properties of these compounds were studied utilizing the Lipinski, Veber, and Ghose criteria. The results revealed that these M. fragrans compounds have potential as effective medicines to combat the COVID-19 pandemic. However, to assess the therapeutic potential of these ligands, additional research is needed, which will use our findings as a foundation.     

Drug repurposing for SARS-CoV-2 main protease: Molecular docking and molecular dynamics investigations

The current novel corona virus illness (COVID-19) is a developing viral disease that was discovered in 2019. There is currently no viable therapeutic strategy for this illness management. Because traditional medication development and discovery has lagged behind the threat of emerging and re-emerging illnesses like Ebola, MERS-CoV, and, more recently, SARS-CoV-2. Drug developers began to consider drug repurposing (or repositioning) as a viable option to the more traditional drug development method. The goal of drug repurposing is to uncover new uses for an approved or investigational medicine that aren't related to its original use. The main benefits of this strategy are that there is less developmental risk and that it takes less time because the safety and pharmacologic requirements are met. The main protease (Mpro) of corona viruses is one of the well-studied and appealing therapeutic targets. As a result, the current research examines the molecular docking of Mpro (PDB ID: 5R81) conjugated repurposed drugs. 12,432 approved drugs were collected from ChEMBL and drugbank libraries, and docked separately into the receptor grid created on 5R81, using the three phases of molecular docking including high throughput virtual screening (HTVS), standard precision (SP), and extra precision (XP). Based on docking scores and MM-GBSA binding free energy calculation, top three drugs (kanamycin, sulfinalol and carvedilol) were chosen for further analyses for molecular dynamic simulations.     

新型三明治样荧光偏振筛选模型在新型冠状病毒主蛋白酶小分子抑制剂筛选中的应用北大核心CSCD

新型冠状病毒主蛋白酶(main protease, Mpro)通过水解多聚蛋白质体(polyprotein)调控病毒基因组RNA复制,且人体不存在其同源蛋白酶,这使Mpro成为抗新型冠状病毒药物开发的理想靶标之一。本研究基于荧光偏振技术(fluorescence polarization,FP)和生物素-亲和素反应(biotin-avidin system, BAS)原理,成功地建立了三明治样荧光偏振筛选模型用于Mpro小分子抑制剂的快速筛选。通过对天然产物化合物库进行高通量筛选,发现了漆树酸(anacardic acid,AA)是Mpro的竞争型抑制剂,1,2,3,4,6-O-五没食子酰葡萄糖(1,2,3,4,6-O-pentagalloylglucose,PGG)是Mpro的混合型抑制剂,且已报道的部分抑制剂是非特异性Mpro小分子抑制剂。文中建立的三明治样荧光偏振筛选模型具有良好的简便性、灵敏性和稳定性,初步证实了漆树酸和PGG是一类新型苗头化合物,建立科学严谨的活性评价体系对于抗新型冠状病毒药物的筛选与发现是至关重要的。     

Assessing Activity and Inhibition of Middle East Respiratory Syndrome Coronavirus Papain-Like and 3C-Like Proteases Using Luciferase-Based Biosensors

Middle East respiratory syndrome coronavirus (MERS-CoV) is associated with an outbreak of more than 90 cases of severe pneumonia with high mortality (greater than 50%). To date, there are no antiviral drugs or specific therapies to treat MERS-CoV. To rapidly identify potential inhibitors of MERS-CoV replication, we expressed the papain-like protease (PLpro) and the 3-chymotrypsin-like protease (3CLpro) from MERS-CoV and developed luciferase-based biosensors to monitor protease activity in cells. We show that the expressed MERS-CoV PLpro recognizes and processes the canonical CoV-PLpro cleavage site RLKGG in the biosensor. However, existing CoV PLpro inhibitors were unable to block MERS-CoV PLpro activity, likely due to the divergence of the amino acid sequence in the drug binding site. To investigate MERS-CoV 3CLpro activity, we expressed the protease in context with flanking nonstructural protein 4 (nsp4) and the amino-terminal portion of nsp6 and detected processing of the luciferase-based biosensors containing the canonical 3CLpro cleavage site VRLQS. Importantly, we found that a small-molecule inhibitor that blocks replication of severe acute respiratory syndrome (SARS) CoV and murine CoV also inhibits the activity of MERS-CoV 3CLpro. Overall, the protease expression and biosensor assays developed here allow for rapid evaluation of viral protease activity and the identification of protease inhibitors. These biosensor assays can now be used to screen for MERS-CoV-specific or broad-spectrum coronavirus PLpro and 3CLpro inhibitors.     

SARS-COV-2 Coronavirus Papain-like Protease PLpro as an Antiviral Target for Inhibitors of Active Site and Protein–Protein Interactions

Biophysics - The papain-like protease PLpro of the SARS-CoV-2 coronavirus is a multifunctional enzyme that catalyzes the proteolytic processing of two viral polyproteins, pp1a and pp1ab. PLpro also...     

Comparison of SARS and NL63 papain-like protease binding sites and binding site dynamics: inhibitor design implications

The human severe acute respiratory syndrome coronavirus (SARS-CoV) and the NL63 coronaviruses are human respiratory pathogens for which no effective antiviral treatment exists. The papain-like cysteine proteases encoded by the coronavirus (SARS-CoV: PLpro; NL63: PLP1 and PLP2) represent potential targets for antiviral drug development. Three recent inhibitor-bound PLpro structures highlight the role of an extremely flexible six-residue loop in inhibitor binding. The high binding site plasticity is a major challenge in computational drug discovery/design efforts. From conventional molecular dynamics and accelerated molecular dynamics (aMD) simulations, we find that with conventional molecular dynamics simulation, PLpro translationally samples the open and closed conformation of BL2 loop on a picosecond-nanosecond timescale but does not reproduce the peptide bond inversion between loop residues Tyr269 and Gln270 that is observed on inhibitor GRL0617 binding. Only aMD simulation, starting from the closed loop conformation, reproduced the 180° ?-ψ dihedral rotation back to the open loop state. The Tyr-Gln peptide bond inversion appears to involve a progressive conformational change of the full loop, starting at one side, and progressing to the other. We used the SARS-CoV apo X-ray structure to develop a model of the NL63-PLP2 catalytic site. Superimposition of the PLP2 model on the PLpro X-ray structure identifies binding site residues in PLP2 that contribute to the distinct substrate cleavage site specificities between the two proteases. The topological and electrostatic differences between the two protease binding sites also help explain the selectivity of non-covalent PLpro inhibitors.     

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.     

The papain-like protease of severe acute respiratory syndrome coronavirus has deubiquitinating activity

Replication of the genomic RNA of severe acute respiratory syndrome coronavirus (SARS-CoV) is mediated by replicase polyproteins that are processed by two viral proteases, papain-like protease (PLpro) and 3C-like protease (3CLpro). Previously, we showed that SARS-CoV PLpro processes the replicase polyprotein at three conserved cleavage sites. Here, we report the identification and characterization of a 316-amino-acid catalytic core domain of PLpro that can efficiently cleave replicase substrates in trans-cleavage assays and peptide substrates in fluorescent resonance energy transfer-based protease assays. We performed bioinformatics analysis on 16 papain-like protease domains from nine different coronaviruses and identified a putative catalytic triad (Cys1651-His1812-Asp1826) and zinc-binding site. Mutagenesis studies revealed that Asp1826 and the four cysteine residues involved in zinc binding are essential for SARS-CoV PLpro activity. Molecular modeling of SARS-CoV PLpro suggested that this catalytic core may also have deubiquitinating activity. We tested this hypothesis by measuring the deubiquitinating activity of PLpro by two independent assays. SARS CoV-PLpro hydrolyzed both diubiquitin and ubiquitin-7-amino-4-methylcoumarin (AMC) substrates, and hydrolysis of ubiquitin-AMC is approximately 180-fold more efficient than hydrolysis of a peptide substrate that mimics the PLpro replicase recognition sequence. To investigate the critical determinants recognized by PLpro, we performed site-directed mutagenesis on the P6 to P2' residues at each of the three PLpro cleavage sites. We found that PLpro recognizes the consensus cleavage sequence LXGG, which is also the consensus sequence recognized by cellular deubiquitinating enzymes. This similarity in the substrate recognition sites should be considered during the development of SARS-CoV PLpro inhibitors.     

Screening of drug databank against WT and mutant main protease of SARS-CoV-2: Towards finding potential compound for repurposing against COVID-19

Although several pharmacological agents are under investigation to be repurposed as therapeutic against COVID-19, not much success has been achieved yet. So, the search for an effective and active option for the treatment of COVID-19 is still a big challenge. The Spike protein (S), RNA-dependent RNA polymerase (RdRp), and Main protease (Mpro) are considered to be the primary therapeutic drug target for COVID-19. In this study we have screened the drugbank compound library against the Main Protease. But our search was not limited to just Mpro. Like other viruses, SARS-CoV-2, have also acquired unique mutations. These mutations within the active site of these target proteins may be an important factor hindering effective drug candidate development. In the present study we identified important active site mutations within the SARS-CoV-2 Mpro (Y54C, N142S, T190I and A191V). Further the drugbank database was computationally screened against Mpro and the selected mutants. Finally, we came up with the common molecules effective against the wild type (WT) and all the selected Mpro. The study found Imiglitazar, was found to be the most active compound against the wild type of Mpro. While PF-03715455 (Y54C), Salvianolic acid A (N142S and T190I), and Montelukast (A191V) were found to be most active against the other selected mutants. It was also found that some other compounds such as Acteoside, 4-Amino-N- {4-[2-(2,6-Dimethyl-Phenoxy)-Acetylamino]-3-Hydroxy-1-Isobutyl-5-Phenyl-Pentyl}-Benzamide, PF-00610355, 4-Amino-N-4-[2-(2,6-Dimethyl-Phenoxy)-Acetylamino]-3-Hydroxy-1-Isobutyl-5-Phenyl-Pentyl}-Benzamide and Atorvastatin were showing high efficacy against the WT as well as other selected mutants. We believe that these molecules will provide a better and effective option for the treatment of COVID-19 clinical manifestations.     

Ligand-based design,synthesis, computational insights,and in vitro studies of novel N-(5-Nitrothiazol-2-yl)-carboxamido derivatives as potent inhibitors of SARS-CoV-2 main protease

The global outbreak of the COVID-19 pandemic provokes scientists to make a prompt development of new effective therapeutic interventions for the battle against SARS-CoV-2. A new series of N-(5-nitrothiazol-2-yl)-carboxamido derivatives were designed and synthesised based on the structural optimisation principle of the SARS-CoV Mpro co-crystallized WR1 inhibitor. Notably, compound 3b achieved the most promising anti-SARS-CoV-2 activity with an IC₅₀ value of 174.7 µg/mL. On the other hand, compounds 3a, 3b, and 3c showed very promising SARS-CoV-2 Mpro inhibitory effects with IC₅₀ values of 4.67, 5.12, and 11.90 µg/mL, respectively. Compound 3b docking score was very promising (−6.94 kcal/mol) and its binding mode was nearly similar to that of WR1. Besides, the molecular dynamics (MD) simulations of compound 3b showed its great stability inside the binding pocket until around 40 ns. Finally, a very promising SAR was concluded to help to design more powerful SARS-CoV-2 Mpro inhibitors shortly.