Pomotrelvir and Nirmatrelvir Binding and Reactivity with SARS-CoV-2 Main Protease: Implications for Resistance Mechanisms from Computations.

We investigate the inhibition mechanism between pomotrelvir and the SARS-CoV-2 main protease using molecular mechanics and quantum mechanics/molecular mechanics simulations. Alchemical transformations where each Pi group of pomotrelvir was transformed into its counterpart in nirmatrelvir were performed to unravel the individual contribution of each group to the binding and reaction processes. We have shown that while a γ-lactam ring is preferred at position P1, a δ-lactam ring is a good alternative for the design of inhibitors for variants presenting mutations at position 166.

Understanding the Interactions of Ubiquitin-Specific Protease 7 with Its Substrates through Molecular Dynamics Simulations: Insights into the Role of Its C-Terminal Domains in Substrate Recognition.

Ubiquitin-specific protease 7 (USP7) is a deubiquitinase enzyme that plays a critical role in regulating various cellular processes by cleaving ubiquitin molecules from target proteins. The C-terminal loop (CTL) motif is a specific region at the C-terminal end of the USP7 enzyme. Recent experiments suggest that the CTL motif plays a role in USP7's catalytic activity by contributing to the enzyme's structural stability, substrate recognition, and catalytic efficiency.

Activation and friction in enzymatic loop opening and closing dynamics.

Protein loop dynamics have recently been recognized as central to enzymatic activity, specificity and stability. However, the factors controlling loop opening and closing kinetics have remained elusive. Here, we combine molecular dynamics simulations with string-method determination of complex reaction coordinates to elucidate the molecular mechanism and rate-limiting step for WPD-loop dynamics in the PTP1B enzyme.

Electrostatics as a Guiding Principle in Understanding and Designing Enzymes.

Enzyme design faces challenges related to the implementation of the basic principles that govern the catalytic activity in natural enzymes. In this work, we revisit basic electrostatic concepts that have been shown to explain the origin of enzymatic efficiency like preorganization and reorganization. Using magnitudes such as the electrostatic potential and the electric field generated by the protein, we explain how these concepts work in different enzymes and how they can be used to rationalize the consequences of point mutations.

Insights into the Enhancement of the Poly(ethylene terephthalate) Degradation by FAST-PETase from Computational Modeling.

Polyethylene terephthalate (PET) is the most abundant polyester plastic, widely used in textiles and packaging, but, unfortunately, it is also one of the most discarded plastics after one use. In the last years, the enzymatic biodegradation of PET has sparked great interest owing to the discovery and subsequent mutation of PETase-like enzymes, able to depolymerize PET. FAST-PETase is one of the best enzymes hitherto proposed to efficiently degrade PET, although the origin of its efficiency is not completely clear.

The impact of SARS-CoV-2 3CL protease mutations on nirmatrelvir inhibitory efficiency. Computational insights into potential resistance mechanisms.

The use of antiviral drugs can promote the appearance of mutations in the target protein that increase the resistance of the virus to the treatment. This is also the case of nirmatrelvir, a covalent inhibitor of the 3CL protease, or main protease, of SARS-CoV-2. In this work we show how the by-residue decomposition of noncovalent interactions established between the drug and the enzyme, in combination with an analysis of naturally occurring mutations, can be used to detect potential mutations in the 3CL protease conferring resistance to nirmatrelvir.

Testing Affordable Strategies for the Computational Study of Reactivity in Cysteine Proteases: The Case of SARS-CoV-2 3CL Protease Inhibition.

Cysteine proteases are an important target for the development of inhibitors that could be used as drugs to regulate the activity of these kinds of enzymes involved in many diseases, including COVID-19. For this reason, it is important to have methodological tools that allow a detailed study of their activity and inhibition, combining computational efficiency and accuracy. We here explore the performance of different quantum mechanics/molecular mechanics methods to explore the inhibition reaction mechanism of the SARS-CoV-2 3CL protease with a hydroxymethyl ketone derivative.

Inhibition Mechanism of SARS-CoV-2 Main Protease with Ketone-Based Inhibitors Unveiled by Multiscale Simulations: Insights for Improved Designs*.

We present the results of classical and QM/MM simulations for the inhibition of SARS-CoV-2 3CL protease by a hydroxymethylketone inhibitor, PF-00835231. In the noncovalent complex the carbonyl oxygen atom of the warhead is placed in the oxyanion hole formed by residues 143 to 145, while P1-P3 groups are accommodated in the active site with interactions similar to those observed for the peptide substrate. According to alchemical free energy calculations, the P1' hydroxymethyl group also contributes to the binding free energy.

Computational simulations on the binding and reactivity of a nitrile inhibitor of the SARS-CoV-2 main protease.

We present a detailed computational analysis of the binding mode and reactivity of the novel oral inhibitor PF-07321332 developed against the SARS-CoV-2 3CL protease. Alchemical free energy calculations suggest that positions P3 and P4 could be susceptible to improvement in order to get a larger binding strength. QM/MM simulations unveil the reaction mechanism for covalent inhibition, showing that the nitrile warhead facilitates the recruitment of a water molecule for the proton transfer step.

Multiscale Simulations of SARS-CoV-2 3CL Protease Inhibition with Aldehyde Derivatives. Role of Protein and Inhibitor Conformational Changes in the Reaction Mechanism.

We here investigate the mechanism of SARS-CoV-2 3CL protease inhibition by one of the most promising families of inhibitors, those containing an aldehyde group as a warhead. These compounds are covalent inhibitors that inactivate the protease, forming a stable hemithioacetal complex. Inhibitor is a potent inhibitor that has been already tested in vitro and in animals.

A microscopic description of SARS-CoV-2 main protease inhibition with Michael acceptors. Strategies for improving inhibitor design.

The irreversible inhibition of the main protease of SARS-CoV-2 by a Michael acceptor known as N3 has been investigated using multiscale methods. The noncovalent enzyme-inhibitor complex was simulated using classical molecular dynamics techniques and the pose of the inhibitor in the active site was compared to that of the natural substrate, a peptide containing the Gln-Ser scissile bond. The formation of the covalent enzyme-inhibitor complex was then simulated using hybrid QM/MM free energy methods.

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

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.

Why are some Enzymes Dimers? Flexibility and Catalysis in Dihydrofolate Reductase.

Dihydrofolate Reductase from (TmDFHFR) is a dimeric thermophilic enzyme that catalyzes the hydride transfer from the cofactor NADPH to dihydrofolate less efficiently than other DHFR enzymes, such as the mesophilic analogue DHFR (EcDHFR). Using QM/MM potentials we show that the reduced catalytic efficiency of TmDHFR is most likely due to differences in the amino acid sequence that stabilize the M20 loop in an open conformation, which prevents the formation of some interactions in the transition state and increases the number of water molecules in the active site. However, dimerization provides two advantages to the thermophilic enzyme; it protects its structure against denaturation by reducing thermal fluctuations and it provides a less negative activation entropy, toning down the increase of the activation free energy with temperature.

Translocation of enzymes into a mesoporous MOF for enhanced catalytic activity under extreme conditions.

Enzymatic catalysis is of great importance to the chemical industry. However, we are still scratching the surface of the potential of biocatalysis due to the limited operating range of enzymes in harsh environments or their low recyclability. The role of Metal-Organic Frameworks (MOFs) as active supports to help overcome these limitations, mainly by immobilization and stabilization of enzymes, is rapidly expanding.

Temperature dependence of dynamic, tunnelling and kinetic isotope effects in formate dehydrogenase.

The origin of the catalytic power of enzymes has been a question of debate for a long time. In this regard, the possible contribution of protein dynamics in enzymatic catalysis has become one of the most controversial topics. In the present work, the hydride transfer step in the formate dehydrogenase (FDH EC 1.

Isotope Substitution of Promiscuous Alcohol Dehydrogenase Reveals the Origin of Substrate Preference in the Transition State.

The origin of substrate preference in promiscuous enzymes was investigated by enzyme isotope labelling of the alcohol dehydrogenase from Geobacillus stearothermophilus (BsADH). At physiological temperature, protein dynamic coupling to the reaction coordinate was insignificant. However, the extent of dynamic coupling was highly substrate-dependent at lower temperatures.

Minimization of dynamic effects in the evolution of dihydrofolate reductase.

Protein isotope labeling is a powerful technique to probe functionally important motions in enzyme catalysis and can be applied to investigate the conformational dynamics of proteins. Previous investigations have indicated that dynamic coupling is detrimental to catalysis by dihydrofolate reductase (DHFR) from the mesophile (EcDHFR). Comparison of DHFRs from organisms adapted to survive at a wide range of temperatures suggests that dynamic coupling in DHFR catalysis has been minimized during evolution; it arises from reorganizational motions needed to facilitate charge transfer events.

Chemical Ligation and Isotope Labeling to Locate Dynamic Effects during Catalysis by Dihydrofolate Reductase.

Chemical ligation has been used to alter motions in specific regions of dihydrofolate reductase from E. coli and to investigate the effects of localized motional changes on enzyme catalysis. Two isotopic hybrids were prepared; one with the mobile N-terminal segment containing heavy isotopes ((2) H, (13) C, (15) N) and the remainder of the protein with natural isotopic abundance, and the other one with only the C-terminal segment isotopically labeled.

Protein isotope effects in dihydrofolate reductase from Geobacillus stearothermophilus show entropic-enthalpic compensatory effects on the rate constant.

Catalysis by dihydrofolate reductase from the moderately thermophilic bacterium Geobacillus stearothermophilus (BsDHFR) was investigated by isotope substitution of the enzyme. The enzyme kinetic isotope effect for hydride transfer was close to unity at physiological temperatures but increased with decreasing temperatures to a value of 1.65 at 5 °C.

Linking electrostatic effects and protein motions in enzymatic catalysis. A theoretical analysis of catechol o-methyltransferase.

The role of protein motions in enzymatic catalysis is the subject of a hot scientific debate. We here propose the use of an explicit solvent coordinate to analyze the impact of environmental motions during the reaction process. The example analyzed here is the reaction catalyzed by catechol O-methyltransferase, a methyl transfer reaction from S-adenosylmethionine (SAM) to the nucleophilic oxygen atom of catecholate.

Heavy enzymes--experimental and computational insights in enzyme dynamics.

The role of protein motions in the chemical step of enzyme-catalyzed reactions is the subject of an open debate in the scientific literature. The systematic use of isotopically substituted enzymes has been revealed as a useful tool to quantify the role of these motions. According to the Born-Oppenheimer approximation, changing the mass of the protein does not change the forces acting on the system but alters the frequencies of the protein motions, which in turn can affect the rate constant.

Increased dynamic effects in a catalytically compromised variant of Escherichia coli dihydrofolate reductase.

Isotopic substitution ((15)N, (13)C, (2)H) of a catalytically compromised variant of Escherichia coli dihydrofolate reductase, EcDHFR-N23PP/S148A, has been used to investigate the effect of these mutations on catalysis. The reduction of the rate constant of the chemical step in the EcDHFR-N23PP/S148A catalyzed reaction is essentially a consequence of an increase of the quasi-classical free energy barrier and to a minor extent of an increased number of recrossing trajectories on the transition state dividing surface. Since the variant enzyme is less well set up to catalyze the reaction, a higher degree of active site reorganization is needed to reach the TS.

Unraveling the role of protein dynamics in dihydrofolate reductase catalysis.

Protein dynamics have controversially been proposed to be at the heart of enzyme catalysis, but identification and analysis of dynamical effects in enzyme-catalyzed reactions have proved very challenging. Here, we tackle this question by comparing an enzyme with its heavy ((15)N, (13)C, (2)H substituted) counterpart, providing a subtle probe of dynamics. The crucial hydride transfer step of the reaction (the chemical step) occurs more slowly in the heavy enzyme.