Abstract:We present a detailed theoretical analysis of the reaction mechanism of proteolysis catalyzed by the main protease of SARS-CoV-2. Using multiscale simulation methods, we have characterized the interactions established by a peptidic substrate in the active site, and then we have explored the free energy landscape associated with the acylation and deacylation steps of the proteolysis reaction, characterizing the transition states of the process. Our mechanistic proposals can explain most of the experimental observations made on the highly similar ortholog protease of SARS-CoV. We point to some key interactions that may facilitate the acylation process and thus can be crucial in the design of more specific and efficient inhibitors of the main protease activity. In particular, from our results, the P1′ residue can be a key factor to improve the thermodynamics and kinetics of the inhibition process.The Supporting Information is available free of charge at <a class="ext-link" href="/doi/10.1021/acscatal.0c03420?goto=supporting-info">https://pubs.acs.org/doi/10.1021/acscatal.0c03420</a>.Plot of the string method convergence; RMSD plots of MD simulations; comparison between the active sites of the two protomers; distribution of S–C distances in the dimer with one and two substrates; a free energy profile corresponding to the acyl–enzyme complex; study of the impact of QM level on the acylation profile; alternative mechanisms for the acylation; deacylation mechanism when His41 is the general base; details of free energy calculations; and a gas-phase evaluation of the performance of the selected functional and basis set (<a class="ext-link" href="/doi/suppl/10.1021/acscatal.0c03420/suppl_file/cs0c03420_si_001.pdf">PDF</a>)Formation of a catalytic dyad ion pair from which the acylation step can proceed (<a class="ext-link" href="/doi/suppl/10.1021/acscatal.0c03420/suppl_file/cs0c03420_si_002.mp4">MP4</a>)Novel mechanism where the N-terminal group of the first peptide fragment acts as the general base catalyzing the hydrolysis of the acyl–enzyme complex (<a class="ext-link" href="/doi/suppl/10.1021/acscatal.0c03420/suppl_file/cs0c03420_si_003.mp4">MP4</a>)This article has not yet been cited by other publications.

Molecular Dynamics Simulations Investigate the Pathway of Substrate Entry Active Site of Rhomboid Protease

An in vitro platform for the enzymatic characterization of the rhomboid protease RHBDL4

A Combination of Solid-State NMR and MD Simulations Reveals the Binding Mode of a Rhomboid Protease Inhibitor

Enhanced activity of Rhizomucor miehei lipase by directed evolution with simultaneous evolution of the propeptide

The rhomboid protease family: a decade of progress on function and mechanism

Structural Analysis of a Rhomboid Family Intramembrane Protease Reveals a Gating Mechanism for Substrate Entry

Structure and Dynamics of the Rhomboid Protease GlpG in Liposomes Studied by Solid-State NMR.

Insight Derived from Molecular Dynamics Simulation into Substrate-Induced Changes in Protein Motions of Proteinase K.

Rhomboid distorts lipids to break the viscosity-imposed speed limit of membrane diffusion

Molecular Docking and Dynamics Simulation Improving Thermophilic Protease Activity of Phpi

Broad substrate specificity and catalytic mechanism of Pseudomonas stutzeri L-rhamnose isomerase: insights from QM/MM molecular dynamics simulations.

In Silico Substrate-Binding Profiling for SARS-CoV-2 Main Protease (Mpro) Using Hexapeptide Substrates

Modification and Simulation of Rhizomucor Miehei Lipase: the Influence of Surficial Electrostatic Interaction on Enantioselectivity.

Unraveling the SARS-CoV-2 Main Protease Mechanism Using Multiscale Methods

Substrate Entering and Product Leaving Trajectories Predict an Engulfing Dynamic for the Major Conformational Change of the &946;-Lactam Acylase

Insight into the activity of SARS main protease: Molecular dynamics study of dimeric and monomeric form of enzyme.

The Conformational Space of the SARS-CoV-2 Main Protease Active Site Loops is Determined by Ligand Binding and Interprotomer Allostery

Probing the interaction mechanism of small molecule inhibitors with matriptase based on molecular dynamics simulation and free energy calculations.

Insight Derived from Molecular Dynamics Simulation into Dynamics and Molecular Motions of Cuticle-Degrading Serine Protease Ver112

Insights Derived from Molecular Dynamics Simulation into the Molecular Motions of Serine Protease Proteinase K

Membrane cholesterol as regulator of human rhomboid protease RHBDL4