Why is the Omicron main protease of SARS-CoV-2 less stable than its wild-type counterpart? A crystallographic, biophysical, and theoretical study of the free enzyme and its complex with inhibitor 13b-K.

During the continuing evolution of SARS-CoV-2, the Omicron variant of concern emerged in the second half of 2021 and has been dominant since November that year. Along with its sublineages, it has maintained a prominent role ever since. The Nsp5 main protease (Mpro) of the Omicron virus is characterized by a single dominant mutation, P132H. Here we determined the X-ray crystal structures of the P132H mutant (or O-Mpro) as free enzyme and in complex with the Mpro inhibitor, the alpha-ketoamide 13b-K, and we conducted enzymology, biophysical as well as theoretical studies to characterize the O-Mpro. We found that O-Mpro has a similar overall structure and binding with 13b-K; however, it displays lower enzymatic activity and lower thermal stability compared to the WT-Mpro (with "WT" referring to the original Wuhan-1 strain). Intriguingly, the imidazole ring of His132 and the carboxylate plane of Glu240 are in a stacked configuration in the X-ray structures determined here. The empirical folding free energy calculations suggest that the O-Mpro dimer is destabilized relative to the WT-Mpro due to the less favorable van der Waals interactions and backbone conformation in the individual protomers. The all-atom continuous constant pH molecular dynamics (MD) simulations reveal that His132 and Glu240 display coupled titration. At pH 7, His132 is predominantly neutral and in a stacked configuration with respect to Glu240 which is charged. In order to examine whether the Omicron mutation eases the emergence of further Mpro mutations, we also determined crystal structures of the relatively frequent P132H+T169S double mutant but found little evidence for a correlation between the two sites.

Integrative Approach to Dissect the Drug Resistance Mechanism of the H172Y Mutation of SARS-CoV-2 Main Protease.

Nirmatrelvir is an orally available inhibitor of SARS-CoV-2 main protease (Mpro) and the main ingredient of Paxlovid, a drug approved by the U.S. Food and Drug Administration for high-risk COVID-19 patients. Recently, a rare natural mutation, H172Y, was found to significantly reduce nirmatrelvir's inhibitory activity. As the COVID-19 cases skyrocket in China and the selective pressure of antiviral therapy builds in the US, there is an urgent need to characterize and understand how the H172Y mutation confers drug resistance. Here, we investigated the H172Y Mpro's conformational dynamics, folding stability, catalytic efficiency, and inhibitory activity using all-atom constant pH and fixed-charge molecular dynamics simulations, alchemical and empirical free energy calculations, artificial neural networks, and biochemical experiments. Our data suggest that the mutation significantly weakens the S1 pocket interactions with the N-terminus and perturbs the conformation of the oxyanion loop, leading to a decrease in the thermal stability and catalytic efficiency. Importantly, the perturbed S1 pocket dynamics weaken the nirmatrelvir binding in the P1 position, which explains the decreased inhibitory activity of nirmatrelvir. Our work demonstrates the predictive power of the combined simulation and artificial intelligence approaches, and together with biochemical experiments, they can be used to actively surveil continually emerging mutations of SARS-CoV-2 Mpro and assist the optimization of antiviral drugs. The presented approach, in general, can be applied to characterize mutation effects on any protein drug targets.

An Integrative Approach to Dissect the Drug Resistance Mechanism of the H172Y Mutation of SARS-CoV-2 Main Protease.

Nirmatrelvir is an orally available inhibitor of SARS-CoV-2 main protease (Mpro) and the main ingredient of PAXLOVID, a drug approved by FDA for high-risk COVID-19 patients. Recently, a rare natural mutation, H172Y, was found to significantly reduce nirmatrelvir's inhibitory activity. As the COVID-19 cases skyrocket in China and the selective pressure of antiviral therapy builds up in the US, there is an urgent need to characterize and understand how the H172Y mutation confers drug resistance. Here we investigated the H172Y Mpro's conformational dynamics, folding stability, catalytic efficiency, and inhibitory activity using all-atom constant pH and fixed-charge molecular dynamics simulations, alchemical and empirical free energy calculations, artificial neural networks, and biochemical experiments. Our data suggests that the mutation significantly weakens the S1 pocket interactions with the N-terminus and perturbs the conformation of the oxyanion loop, leading to a decrease in the thermal stability and catalytic efficiency. Importantly, the perturbed S1 pocket dynamics weakens the nirma-trelvir binding in the P1 position, which explains the decreased inhibitory activity of nirmatrelvir. Our work demonstrates the predictive power of the combined simulation and artificial intel-ligence approaches, and together with biochemical experiments they can be used to actively surveil continually emerging mutations of SARS-CoV-2 Mpro and assist the discovery of new antiviral drugs. The presented workflow can be applicable to characterize mutation effects on any protein drug targets.

H172Y mutation perturbs the S1 pocket and nirmatrelvir binding of SARS-CoV-2 main protease through a nonnative hydrogen bond.

Nirmatrelvir is an orally available inhibitor of SARS-CoV-2 main protease (Mpro) and the main ingredient of PAXLOVID, a drug approved by FDA for high-risk COVID-19 patients. Although the prevalent Mpro mutants in the SARS-CoV-2 Variants of Concern (e.g., Omicron) are still susceptible to nirmatrelvir, a rare natural mutation, H172Y, was found to significantly reduce nirmatrelvir's inhibitory activity. As the selective pressure of antiviral therapy may favor resistance mutations, there is an urgent need to understand the effect of the H172Y mutation on Mpro's structure, function, and drug resistance. Here we report the molecular dynamics (MD) simulations as well as the measurements of stability, enzyme kinetics of H172Y Mpro, and IC50 value of nirmatrelvir. Simulations showed that mutation disrupts the interactions between the S1 pocket and N terminus of the opposite protomer. Intriguingly, a native hydrogen bond (H-bond) between Phe140 and the N terminus is replaced by a transient H-bond between Phe140 and Tyr172. In the ligand-free simulations, strengthening of this nonnative H-bond is correlated with disruption of the conserved aromatic stacking between Phe140 and His163, leading to a partial collapse of the oxyanion loop. In the nirmatrelvir-bound simulations, the nonnative H-bond is correlated with the loss of an important H-bond between Glu166 and nirmatrelvir's lactam nitrogen at P1 position. These results are consistent with the newly reported X-ray structures of H172Y Mpro and suggest a mechanism by which the H172Y substitution perturbs the S1 pocket, leading to the decreased structural stability and binding affinity, which in turn explains the drastic reduction in catalytic activity and antiviral susceptibility.

A Guide to the Continuous Constant pH Molecular Dynamics Methods in Amber and CHARMM [Article v1.0].

Like temperature and pressure, solution pH is an important environmental variable in biomolecular simulations. Virtually all proteins depend on pH to maintain their structure and function. In conventional molecular dynamics (MD) simulations of proteins, pH is implicitly accounted for by assigning and fixing protonation states of titratable sidechains. This is a significant limitation, as the assigned protonation states may be wrong and they may change during dynamics. In this tutorial, we guide the reader in learning and using the various continuous constant pH MD methods in Amber and CHARMM packages, which have been applied to predict pK a values and elucidate proton-coupled conformational dynamics of a variety of proteins including enzymes and membrane transporters.

Electrostatic interaction optimization improves catalytic rates and thermotolerance on xylanases.

Understanding the aspects that contribute to improving proteins' biochemical properties is of high relevance for protein engineering. Properties such as the catalytic rate, thermal stability, and thermal resistance are crucial for applying enzymes in the industry. Different interactions can influence those biochemical properties of an enzyme. Among them, the surface charge-charge interactions have been a target of particular attention. In this study, we employ the Tanford-Kirkwood solvent accessibility model using the Monte Carlo algorithm (TKSA-MC) to predict possible interactions that could improve stability and catalytic rate of a WT xylanase (XynAWT) and its M6 xylanase (XynAM6) mutant. The modeling prediction indicates that mutating from a lysine in position 99 to a glutamic acid (K99E) favors the native state stabilization in both xylanases. Our lab results showed that mutated xylanases had their thermotolerance and catalytic rate increased, which conferred higher processivity of delignified sugarcane bagasse. The TKSA-MC approach employed here is presented as an efficient computational-based design strategy that can be applied to improve the thermal resistance of enzymes with industrial and biotechnological applications.