Theoretical study of the inhibition mechanism of human 20S proteasome by dihydroeponemycin

Electrostatic Preorganization in Three Distinct Heterogeneous Proteasome β-Subunits

The origin of the enzyme's powerful role in accelerating chemical reactions is one of the most critical and still widely discussed questions. It is already accepted that enzymes impose an electrostatic field onto their substrates by adopting complex three-dimensional structures; therefore, the preorganization of electric fields inside protein active sites has been proposed as a crucial contributor to catalytic mechanisms and rate constant enhancement. In this work, we focus on three catalytically active β-subunits of 20S proteasomes with low sequence identity (∼30%) whose active sites, although situated in an electrostatically miscellaneous environment, catalyze the same chemical reaction with similar catalytic efficiency. Our in silico experiments reproduce the experimentally observed equivalent reactivity of the three sites and show that obliteration of the electrostatic potential in all active sites would deprive the enzymes of their catalytic power by slowing down the chemical process by a factor of 1035. To regain enzymatic efficiency, besides catalytic Thr1 and Lys33 residues, the presence of aspartic acid in position 17 and an aqueous solvent is required, proving that the electrostatic potential generated by the remaining residues is insignificant for catalysis. Moreover, it was found that the gradual decay of atomic charges on Asp17 strongly correlates with the enzyme's catalytic rate deterioration as well as with a change in the charge distributions due to introduced mutations. The computational procedure used and described here may help identify key residues for catalysis in other biomolecular systems and consequently may contribute to the process of designing enzyme-like synthetic catalysts.

Discovery of Novel Nonpeptidic Proteasome Inhibitors Using Covalent Virtual Screening and Biological Evaluation

Many reported proteasome inhibitors, including the three clinically approved inhibitors, bortezomib, carfilzomib, and ixazomib, have peptidic structures. In this study, using a hybrid and versatile strategy for covalent virtual screening by combining warhead screening and preprocessing with GOLD and CovDock software that were applied to the ZINC virtual library, we identified multiple proteasome inhibitors with new nonpeptidic structural scaffolds. Proteasome inhibition assays confirmed the inhibitory activities of these new compounds. These results demonstrate the effectiveness of our computational strategy for large-scale covalent virtual screening. Furthermore, these identified proteasome inhibitors may serve as starting points for the development of a new class of nonpeptidic therapeutic agents.

Elucidation of the α-Ketoamide Inhibition Mechanism: Revealing the Critical Role of the Electrostatic Reorganization Effect of Asp17 in the Active Site of the 20S Proteasome

The 20S proteasome is an attractive drug target for the development of anticancer agents because it plays an important role in cellular protein degradation. It has a threonine residue that can act as a nucleophile to attack inhibitors with an electrophilic warhead, forming a covalent adduct. Fundamental understanding of the reaction mechanism between covalent inhibitors and the proteasome may assist the design and refinement of compounds with the desired activity. In this study, we investigated the covalent inhibition mechanism of an α-keto phenylamide inhibitor of the proteasome. We calculated the noncovalent binding free energy using the PDLD/S-LRA/β method and the reaction free energy through the empirical valence bond method (EVB). Several possible reaction pathways were explored. Subsequently, we validated the calculated activation and reaction free energies of the most plausible pathways by performing kinetic experiments. Furthermore, the effects of different ionization states of Asp17 on the activation energy at each step were also discussed. The results revealed that the ionization states of Asp17 remarkably affect the activation energies and there is an electrostatic reorganization of Asp17 during the course of the reaction. Our results demonstrate the critical electrostatic effect of Asp17 in the active site of the 20S proteasome.

Covalent Inhibition of the Human 20S Proteasome with Homobelactosin C Inquired by QM/MM Studies

20S proteasome is a main player in the protein degradation pathway in the cytosol, thus intervening in multiple pivotal cellular processes. Over the years the proteasome has emerged as a crucial target for the treatment of many diseases such as neurodegenerative diseases, cancer, autoimmune diseases, developmental disorders, cystic fibrosis, diabetes, cardiac diseases, atherosclerosis, and aging. In this work, the mechanism of proteasome covalent inhibition with bisbenzyl-protected homobelactosin C (hBelC) was explored using quantum mechanics/molecular mechanics (QM/MM) methods. Molecular dynamic simulations were used to describe key interactions established between the hBelC and its unique binding mode in the primed site of the β5 subunit. The free energy surfaces were computed to characterize the kinetics and thermodynamics of the inhibition process. This study revealed that although the final inhibition product for hBelC is formed according to the same molecular mechanism as one described for hSalA, the free energy profile of the reaction pathway differs significantly from the one previously reported for γ-lactam-β-lactone containing inhibitors in terms of the height of the activation barrier as well as the stabilization of the final product. Moreover, it was proved that high stabilization of the covalent adduct formed between β5-subunit and hBelC, together with the presence of aminocarbonyl side chain in the structure of the inhibitor which prevents the hydrolysis of the ester bond from taking place, determines its irreversible character.

Mechanistic and thermodynamic characterization of oxathiazolones as potent and selective covalent immunoproteasome inhibitors

The ubiquitin–proteasome system is responsible for the degradation of proteins and plays a critical role in key cellular processes. While the constitutive proteasome (cPS) is expressed in all eukaryotic cells, the immunoproteasome (iPS) is primarily induced during disease processes, and its inhibition is beneficial in the treatment of cancer, autoimmune disorders and neurodegenerative diseases. Oxathiazolones were reported to selectively inhibit iPS over cPS, and the inhibitory activity of several oxathiazolones against iPS was experimentally determined. However, the detailed mechanism of the chemical reaction leading to irreversible iPS inhibition and the key selectivity drivers are unknown, and separate characterization of the noncovalent and covalent inhibition steps is not available for several compounds. Here, we investigate the chemical reaction between oxathiazolones and the Thr1 residue of iPS by quantum mechanics/molecular mechanics (QM/MM) simulations to establish a plausible reaction mechanism and to determine the rate-determining step of covalent complex formation. The modelled binding mode and reaction mechanism are in line with the selective inhibition of iPS versus cPS by oxathiazolones. The k_inact value of several ligands was estimated by constructing the potential of mean force of the rate-determining step by QM/MM simulations coupled with umbrella sampling. The equilibrium constant K_i of the noncovalent complex formation was evaluated by classical force field-based thermodynamic integration. The calculated K_i and k_inact values made it possible to analyse the contribution of the noncovalent and covalent steps to the overall inhibitory activity. Compounds with similar intrinsic reactivities exhibit varying selectivities for iPS versus cPS owing to subtle differences in the binding modes that slightly affect K_i, the noncovalent affinity, and importantly alter k_inact, the covalent reactivity of the bound compounds. A detailed understanding of the inhibitory mechanism of oxathiazolones is useful in designing iPS selective inhibitors with improved drug-like properties.

Simulating the directional translocation of a substrate by the AAA+ motor in the 26S proteasome

This work explored the molecular origin of substrate translocation by the AAA+ motor of the 26S proteasome. This exploration was performed by combining different simulation approaches including calculations of binding free energies, coarse-grained simulations, and considerations of the ATP hydrolysis energy. The simulations were used to construct the free energy landscape for the translocation process. This included the evaluation of the conformational barriers in different translocation steps. Our simulation reveals that the substrate translocation by the AAA+ motor is guided in part by electrostatic interactions. We also validated the experimental observation that bulkier residues in pore loop 1 are responsible for substrate translocation. However, our calculation also reveals that the lysine residues prior to the bulkier residues (conserved along pore loop 1) are also important for the translocation process. We believe that this computational study can help in guiding the ongoing research of the proteasome.

Hydrogen bonding catalysis by water in epoxide ring opening reaction

Water can act as catalyst is perhaps the most intriguing property reported of this molecule in the last decade. However, despite being an integral part of many enzyme structures, the role of water in catalyzing enzymatic reactions remains sparsely studied. In a recent study, we have shown that the epoxide ring opening in aspartate proteases follows a two-step process involving water. In this work, we attempt to unravel the electronic basis of the co-catalytic role of water in the epoxide ring opening reaction by employing high-level quantum mechanical calculations at M06-2X/6-31+G(d,p) level of accuracy. Our computed electron density and its reduced gradient show that water anchor the reactant molecules through strong H-bond bridges. In addition, the strong ionizing power of water allows better charge delocalization to stabilize the transition states and oxyanion intermediate. Electrostatic analyses suggest greater charge transfer from the aspartates to the epoxide in the transition state, which is found to be exergonic in nature rendering a low-barrier reaction compared to a control system where water was omitted in the reaction field. This elucidated mechanism at electronic level could promote further research to search for the co-catalytic role of water in other enzymes.

Theoretical Studies of the Acid-Base Equilibria in a Model Active Site of the Human 20S Proteasome

The 20S proteasome is a macromolecule responsible for the chemical step in the ubiquitin-proteasome system of degrading unnecessary and unused proteins of the cell. It plays a central role both in the rapid growth of cancer cells and in viral infection cycles. Herein, we present a computational study of the acid-base equilibria in an active site of the human proteasome (caspase-like), an aspect which is often neglected despite the crucial role protons play in the catalysis. As example substrates, we take the inhibition by epoxy- and boronic acid-containing warheads. We have combined cluster quantum mechanical calculations, replica exchange molecular dynamics, and Bayesian optimization of nonbonded potential terms in the inhibitors. In relation to the latter, we propose an easily scalable approach for the reevaluation of nonbonded potentials making use of the hybrid quantum mechanics molecular mechanics dynamics information. Our results show that coupled acid-base equilibria need to be considered when modeling the inhibition mechanism. The coupling between a neighboring lysine and the reacting threonine is not affected by the presence of the studied inhibitors.

Mechanism of inhibition of SARS-CoV-2 Mpro by peptidyl Michael acceptor explained by QM/MM simulations and design of new derivatives with tunable chemical reactivity

The SARS-CoV-2 main protease (Mpro) is essential for replication of the virus responsible for the COVID-19 pandemic, and one of the main targets for drug design. Here, we simulate the inhibition process of SARS-CoV-2 Mpro with a known Michael acceptor (peptidyl) inhibitor, N3. The free energy landscape for the mechanism of the formation of the covalent enzyme-inhibitor product is computed with QM/MM molecular dynamics methods. The simulations show a two-step mechanism, and give structures and calculated barriers in good agreement with experiment. Using these results and information from our previous investigation on the proteolysis reaction of SARS-CoV-2 Mpro, we design two new, synthetically accessible N3-analogues as potential inhibitors, in which the recognition and warhead motifs are modified. QM/MM modelling of the mechanism of inhibition of Mpro by these novel compounds indicates that both may be promising candidates as drug leads against COVID-19, one as an irreversible inhibitor and one as a potential reversible inhibitor.

Exploring the Proteolysis Mechanism of the Proteasomes

The proteasome is a key protease in the eukaryotic cells which is responsible for various important cellular processes such as the control of the cell cycle, immune responses, protein homeostasis, inflammation, apoptosis, and the response to proteotoxic stress. Acting as a major molecular machine for protein degradation, proteasome first identifies damaged or obsolete regulatory proteins by attaching ubiquitin chains and subsequently utilizes conserved pore loops of the heterohexameric ring of AAA+ (ATPases associated with diverse cellular activities) to pull and mechanically unfold and translocate the misfolded protein to the active site for proteolysis. A detailed knowledge of the reaction mechanism for this proteasomal proteolysis is of central importance, both for fundamental understanding and for drug discovery. The present study investigates the mechanism of the proteolysis by the proteasome with full consideration of the protein's flexibility and its impact on the reaction free energy. Major attention is paid to the role of the protein electrostatics in determining the activation barriers. The reaction mechanism is studied by considering a small artificial fluorogenic peptide substrate (Suc-LLVY-AMC) and evaluating the activation barriers and reaction free energies for the acylation and deacylation steps, by using the empirical valence bond method. Our results shed light on the proteolysis mechanism and thus should be important for further studies of the proteasome action.

Water Plays a Cocatalytic Role in Epoxide Ring Opening Reaction in Aspartate Proteases: A QM/MM Study

Aspartate proteases are potential targets for various diseases, and many of their inhibitors are FDA-approved drugs. However, these peptidomimetic and reversibly bound drugs become ineffective upon prolonged use. Attempts have been made to design and synthesize various nonpeptidic epoxide-based irreversible inhibitors to combat the drug-resistance enigma. Here, we study the mechanism of epoxide ring opening in two widely studied aspartate proteases, HIV-1 protease and pepsin. Our results from QM/MM molecular dynamics show that the epoxide ring opening in aspartate proteases follow a two-step mechanism with the formation of an oxyanion intermediate, stabilized by a set of water molecules in the protein active site. These water molecules by virtue of "low-barrier hydrogen bonds" with the epoxide ring reduce the intrinsic reaction barrier while remaining structurally unperturbed, thus playing a cocatalytic role. We validated our results by reproducing the experimentally observed protease/pepsin-epoxide covalent complexes as end products. The observed stability of our oxyanion intermediate in a four-water-coordinated state is also consistent with the reported stable state of the hydroxide ion in water as OH^-(H₂O)₄. Our study could pave the way for the design of new class "HIV protease irreversible inhibitors" from the acquired knowledge of the structures of intermediate and transition states traced during the explored reaction mechanism.