Targeting omicron and other reported SARS-CoV-2 lineages by potent inhibitors of main protease 3CL Mpro: Molecular simulation analysis

Abstract

Besides vaccines, anti-SARS-CoV-2 drugs are urgently needed to cease the pandemic. Clinical use of its protease inhibitor therapy for COVID-19 has been reported in this Journal by Deng et al.1Deng L. et al.Arbidol combined with LPV/r versus LPV/r alone against corona virus disease 2019: a retrospective cohort study.J Infect. 2020; 81: e1-e5Abstract Full Text Full Text PDF PubMed Scopus (322) Google Scholar The SARS-CoV-2 3CL Mpro protease plays essential roles in viral replication, transcription, inducing cytokines IL-1β, IL-6, TNF-α, in Calu-3 and THP1 cells.2Li W. et al.SARS-CoV-2 Nsp5 activates NF-kappaB pathway by upregulating SUMOylation of MAVS.Front Immunol. 2021; 12750969Google Scholar To comment on the utilization of protease inhibitors, we here identify and probe other potential candidates besides Arbidol and LPV/r.1Deng L. et al.Arbidol combined with LPV/r versus LPV/r alone against corona virus disease 2019: a retrospective cohort study.J Infect. 2020; 81: e1-e5Abstract Full Text Full Text PDF PubMed Scopus (322) Google Scholar We picked the Mpro structure in full-length that retains a cysteine containing catalytic pocket, unlike chymotrypsin-like enzymes and performed structure-based virtual screening and molecular dynamic simulation on 25 FDA approved antiviral and 4000 unknown lead molecules. Among the reported inhibitors silymarin, ferulic acid and p-coumaric acid identified as best binders with overall energies of -15.0411 Kcal/mol, -10.6320 Kcal/mol and -10.7634 Kcal/mol, respectively. Contrary, from the list of the unknown compounds, ZINC64219754, ZINC64219742 and ZINC64219736 showed binding affinities of -11.6969 Kcal/mol, -11.1067 Kcal/mol and -10.7840 Kcal/mol (Table 1). In line, silymarin, ferulic acid and p-coumaric acid have already shown high potency as antiviral drugs against flaviviruses, togaviruses, adenoviruses and dengue virus.3Liu C.H. et al.Antiviral activities of silymarin and derivatives.Molecules. 2019; 24: 1-15Crossref Scopus (40) Google Scholar,4Musarra-Pizzo M. et al.Antiviral activity exerted by natural products against human viruses.Viruses. 2021; 13: 1-30Crossref Scopus (18) Google Scholar Silymarin engaged Thr26, Asn142, Ser144, Cys145 and His163 and contributed to six hydrogen bonds and van der Waal (vdW) interactions. In the case of Ferulic Acid and p-coumaric acid, besides vdW interactions, six and four hydrogen bonds were observed, respectively. Interestingly, all the three known compounds stayed consistent in capturing the catalytic residue Cys145, which was recently a reported case with MG-132 against 3CL Mpro.5Costanzi E. et al.Structural and biochemical analysis of the dual inhibition of MG-132 against SARS-CoV-2 main protease (Mpro/3CLpro) and human cathepsin-L.Int J Mol Sci. 2021; 22Crossref Scopus (6) Google ScholarTable 1List of the scores and interactions retrieved during docking analyses.Name/Zinc IDS-scoreHydrogen bondsVan der Waals and Hydrophobic interactionsSilymarin-15.0411Thr26, Asn142, Ser144, Cys145, His163Thr24, Thr25, Phe140, Leu141, Gly143, Glu166Ferulic Acid-10.6320Leu141, Gly143, Ser144, Cys145, His163, Glu166Phe140, Asn142, Met165p-coumaric acid-10.7634Phe140, Cys145, GLU166, HIS172Leu141, Asn142, Gly143, Ser144, His163ZINC64219754-11.6969His41, Gly143Leu27, Met49, Phe140, Leu141, Asn142, Cys145, His164, Met165, Glu166, Arg188, Gln189, Thr190, Gln192ZINC64219742-11.1067His41, Cys145Thr25, Thr26, Leu27, Cys44, Thr45, Ser46, Met49, Phe140, Leu141, Asn142, Gly143, Ser144, His163, His164, Glu166ZINC64219736-10.7840His41, Asn142, Gly143, His164Leu27, Phe140, Leu141, Gly143, Ser144, Cys145, His163, Met165, Glu166 Open table in a new tab The docked ZINC64219754 depicted two hydrogen bonds coupled with vdW interactions, where one of the interactions was achieved with catalytic His41 compared to the Cys145, observed in known drugs. The presence of multiple carboxylate groups makes the environment of the compound more electrophilic and facilitates the formation of electrostatic interactions with active site residues. The imidazole Nε1 atom of His41 forms a hydrogen contact with conserved H2Ocat molecule whereas the Nε2 atom engaged carbonyl oxygen of the inhibitor. To our interest, the compound ZINC64219742 captured both the catalytic residues, His41 and Cys145 at distances of 2.5 Å and 3.1 Å. Finally, the interaction pattern of ZINC64219736 was observed to constitute three hydrogen bonds and vdW interaction with amino acids residues of His41, Asn142, Gly143 and His164 at the range of 2.1 – 2.9 Å (Fig. 1A). Overall, our docking analyses suggest that all the three unknown binders retain the capability to arrest the catalytic residues and in particular ZINC64219742 indicated the ability of capturing the dyad. Accessing the stability of protein–ligand interactions and ligand-induced changes in the protein structure, all the six 3CL Mpro docked complexes were tested in 50ns each of MD simulation. All the complexes successfully attained the equilibration state and stayed in the RMSD range of 1 – 4 Å, suggesting a stable behavior of the Mpro with all the six compounds (Fig. 1B). Analyzing the interaction panel for ferulic acid, we found that only two H-bond interactions with Glu166 and Gly143 stayed stable in simulation, while interaction with Cys145 and Ser144 were lost in the first 2ns range. Contrary, two new and relatively stable interactions were attained with Phe140 and His163. On the other hand, Silymarin showed the best fitness and achieved few more H-bond interactions but interestingly all the interactions remained highly stable. Apart from Glu166, p-coumaric acid replaced all the interactions achieved in docking with Gln189, His41, Ser46 and Tyr54. The ZINC compounds, ZINC64219736 formed two new stable contacts with Asn277 and Asn238 along with other part time interactions but unfortunately have lost all the contacts observed in docked complexes. Similar to its docked complex, ZINC64219742 and ZINC64219754 stayed in the active pocket region through hydrophobic interactions (Fig. 2A). To conclude, status of the tested compounds during simulation, silymarin appeared the most native compound in the active pocket of Mpro forming H-interactions, followed by the p-coumaric and ferulic acid. Contrary, all the three compounds from the ZINC database stayed in the active-pocket relying on hydrophobic interactions. It is noteworthy that all the complexes in simulation have attained equilibrated state whereas the RMSF profile did not indicate any spontaneous jump, suggesting the stable nature of these compounds in the active pocket region. Many variants of SARS-Cov-2 have been reported since pandemic including alpha to omicron. Based on the observed number of mutations in these variants, we found that only beta and omicron variants harbor single mutations each, i.e., K90R and P132H, respectively (Fig. 2B). Both the mutations are distant and not involved in any interaction. Recently, two existing drugs talampicillin and lurasidone have shown binding energies of -11.17 Kcal/mol each and 2 drug-like compounds ZINC000015988935 and ZINC000103558522 exhibited binding energies of -12.39 and -12.36 Kcal/mol which suggests be 3CL Mpro protease inhibitors.6Elmezayen A.D. et al.Drug repurposing for coronavirus (COVID-19): in silico screening of known drugs against coronavirus 3CL hydrolase and protease enzymes.J Biomol Struct Dyn. 2021; 39: 2980-2992Crossref PubMed Scopus (224) Google Scholar We observed that in the presence of H2Ocat the His41 and Asp187 refrained to involve the direct H-interaction. In a study, candidates with binding energies -6.1 to -7.75 kcal/mol were suggested better inhibition than nelfinavir based on the target interactions7Sargolzaei M. Effect of nelfinavir stereoisomers on coronavirus main protease: molecular docking, molecular dynamics simulation and MM/GBSA study.J Mol Gr Model. 2021; 103107803PubMed Google Scholar and among the 36 tested derivatives 20c, 24c, 30c, 34c, 35c, and 36c were identified best inhibitors based on the target receptor interactions with Gln189, Cys145, and His41 found significant. Derivatives of 2.5-diaminobenzophenone chosen molecules that achieved hydrophobic and vdW interactions with Thr25, Thr26, Leu27, Cys44, Thr45, Ser46, Met49, Phe140, Leu141, Asn142, Gly143, Ser144, His163, residues and H-bond interactions with residues His41, Leu141, Gly143, Ser144, Cys145, His163 and His164 were suggested potential candidates.8Ouassaf M. et al.Combined docking methods and molecular dynamics to identify effective antiviral 2, 5-diaminobenzophenonederivatives against SARS-CoV-2.J King Saud Univ Sci. 2021; 33101352Crossref PubMed Scopus (8) Google Scholar Silymarin could also inhibit the dengue virus in vitro and were well tolerated in Vero cells with a half-maximum cytotoxic concentration (CC50) of 749.70 µg/mL against DENV-3, and prevented viral entry (72.46%) into the cells with binding to the viral envelope.9Low Z.X. et al.Antiviral activity of silymarin and baicalein against dengue virus.Sci Rep. 2021; 11: 21221Crossref PubMed Scopus (4) Google Scholar Overall, we further suggest six possible candidates and particularly silymarin that attained 15.0411Kcal/mol in docking and gained further fitness in terms of H-bonds in the pocket as compared to the other five candidates, suggesting the potent clinically targeting omicron through main protease 3CL Mpro besides targeting entry.10Kazybay B. Ahmad A. Mu C. Mengdesh D. Xie Y. Omicron N501Y mutation among SARS-CoV-2 lineages: insilico analysis of potent binding to tyrosine kinase and hypothetical repurposed medicine.Travel Med Infect Dis. 2022; 45 (Jan-FebEpub 2021 Dec 17. PMID: 34929375; PMCID: PMC8677628)102242https://doi.org/10.1016/j.tmaid.2021.102242Crossref PubMed Scopus (8) Google Scholar N/A. N/A. All the data materials mentioned are available.