Atomistic-Level Description of the Covalent Inhibition of SARS-CoV-2 Papain-like Protease

Abstract

Inhibition of the papain-like protease (PLpro) of SARS-CoV-2 has been demonstrated to be a successful target to prevent the spreading of the coronavirus in the infected body. In this regard, covalent inhibitors, such as the recently proposed VIR251 ligand, can irreversibly inactivate PLpro by forming a covalent bond with a specific residue of the catalytic site (Cys¹¹¹), through a Michael addition reaction. An inhibition mechanism can therefore be proposed, including four steps: (i) ligand entry into the protease pocket; (ii) Cys¹¹¹ deprotonation of the thiol group by a Brønsted-Lowry base; (iii) Cys¹¹¹-S^- addition to the ligand; and (iv) proton transfer from the protonated base to the covalently bound ligand. Evaluating the energetics and PLpro conformational changes at each of these steps could aid the design of more efficient and selective covalent inhibitors. For this aim, we have studied by means of MD simulations and QM/MM calculations the whole mechanism. Regarding the first step, we show that the inhibitor entry in the PLpro pocket is thermodynamically favorable only when considering the neutral Cys¹¹¹, that is, prior to the Cys¹¹¹ deprotonation. For the second step, MD simulations revealed that His²⁷² would deprotonate Cys¹¹¹ after overcoming an energy barrier of ca. 32 kcal/mol (at the QM/MM level), but implying a decrease of the inhibitor stability inside the protease pocket. This information points to a reversible Cys¹¹¹ deprotonation, whose equilibrium is largely shifted toward the neutral Cys¹¹¹ form. Although thermodynamically disfavored, if Cys¹¹¹ is deprotonated in close proximity to the vinylic carbon of the ligand, then covalent binding takes place in an irreversible way (third step) to form the enolate intermediate. Finally, due to Cys¹¹¹-S^- negative charge redistribution over the bound ligand, proton transfer from the initially protonated His²⁷² is favored, finally leading to an irreversibly modified Cys¹¹¹ and a restored His²⁷². These results elucidate the selectivity of Cys¹¹¹ to enable formation of a covalent bond, even if a weak proton acceptor is available, as His²⁷².

Keywords: SARS-CoV-2; covalent inhibitor; free energy calculations; molecular dynamics; papain-like protease.