Computational analysis of the interactions between Ebselen and derivatives with the active site of the main protease from SARS-CoV-2

The main protease (M^pro) of the novel coronavirus SARS-CoV-2 is a key target for developing antiviral drugs. Ebselen (EbSe) is a selenium-containing compound that has been shown to inhibit Mpro in vitro by forming a covalent bond with the cysteine (Cys) residue in the active site of the enzyme. However, EbSe can also bind to other proteins, like albumin, and low molecular weight compounds that have free thiol groups, such as Cys and glutathione (GSH), which may affect its availability and activity. In this study, we analyzed the M^pro interaction with EbSe, its analogues, and its metabolites with Cys, GSH, and albumin by molecular docking. We also simulated the electronic structure of the generated molecules by density functional theory (DFT) and explored the stability of EbSe and one of its best derivatives, EbSe-2,5-MeClPh, in the catalytic pocket of Mpro through covalent docking and molecular dynamics. Our results show that EbSe and its analogues bound to GSH/albumin have larger distance between the selenium atom of the ligands and the sulfur atom of Cys145 of M^pro than the other compounds. This suggests that EbSe and its GSH/albumin-analogues may have less affinity for the active site of Mpro. EbSe-2,5-MeClPh was found one of the best molecules, and in molecular dynamics simulations, it showed to undergo more conformational changes in the active site of Mpro, in relation to EbSe, which remained stable in the catalytic pocket. Moreover, this study also reveals that all compounds have the potential to interact closely with the active site of Mpro, providing us with a concept of which derivatives may be promising for in vitro analysis in the future. We propose that these compounds are potential covalent inhibitors of M^pro and that organoselenium compounds are molecules that should be studied for their antiviral properties.