Enzyme kinetics model for the coronavirus main protease including dimerization and ligand binding

The coronavirus main protease (MPro) plays a pivotal role in viral replication and is the target of several antivirals against SARS-CoV-2. In some species, CRCs of MPro enzymatic activity can exhibit biphasic behavior in which low ligand concentrations activate the enzyme whereas higher ones inhibit it. While this behavior has been attributed to ligand-induced dimerization, quantitative enzyme kinetics models have not been fit to it. Here, we develop a kinetic model integrating dimerization and ligand binding. We perform a Bayesian regression to globally fit the model to multiple types of biochemical and biophysical data. The reversible covalent inhibitor GC376 strongly induces dimerization and binds to the dimer with no cooperativity. In contrast, the fluorescent peptide substrate has a minor effect on dimerization but binds to the dimer with positive cooperativity. The biphasic concentration response curve occurs because compared to substrate, the inhibitor accelerates turnover in the opposite catalytic site.