Fundamental reaction pathway for peptide metabolism by proteasome: insights from first-principles quantum mechanical/molecular mechanical free energy calculations

Landscape of In Silico Tools for Modeling Covalent Modification of Proteins: A Review on Computational Covalent Drug Discovery

Covalent drug discovery has been a challenging research area given the struggle of finding a sweet balance between selectivity and reactivity for these drugs, the lack of which often leads to off-target activities and hence undesirable side effects. However, there has been a resurgence in covalent drug design following the success of several covalent drugs such as boceprevir (2011), ibrutinib (2013), neratinib (2017), dacomitinib (2018), zanubrutinib (2019), and many others. Design of covalent drugs includes many crucial factors, where "evaluation of the binding affinity" and "a detailed mechanistic understanding on covalent inhibition" are at the top of the list. Well-defined experimental techniques are available to elucidate these factors; however, often they are expensive and/or time-consuming and hence not suitable for high throughput screens. Recent developments in in silico methods provide promise in this direction. In this report, we review a set of recent publications that focused on developing and/or implementing novel in silico techniques in "Computational Covalent Drug Discovery (CCDD)". We also discuss the advantages and disadvantages of these approaches along with what improvements are required to make it a great tool in medicinal chemistry in the near future.

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.

The chiral pyridoxal-catalyzed biomimetic Mannich reaction: the mechanism and origin of stereoselectivity

A biomimetic organocatalyst with a pyridoxal-like structure is one of the most successful examples of catalyzing organic reactions under mild conditions in an asymmetric synthesis field.

Organocatalytic insertion into C–B bonds by in situ generated carbene: mechanism, role of the catalyst, and origin of stereoselectivity

The calculated results showed that C–H⋯X (X = I, Br) hydrogen bond interactions should be the determining factors of stereoselectivity, and FMO overlap mode and ELF analyses along IRC results were performed to explore the nature of carbene insertion.

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.

N-Heterocyclic Carbene-Catalyzed Asymmetric C-O Bond Construction Between Benzoic Acid and o-Phthalaldehyde: Mechanism and Origin of Stereoselectivity

A computational study was contributed to explore the origin of stereoselectivity of NHC-mediated cyclization reaction between benzoic acid and o-phthalaldehyde for asymmetric construction of phthalidyl ester. The most energetically favorable pathway mainly includes the following steps: (1) nucleophilic attack on carbonyl carbon of o-phthalaldehyde by catalyst NHC, (2) formation of Breslow intermediate, (3) oxidation by DQ, (4) asymmetric formation of dual C-O bonds, and (5) dissociation of catalyst with the product. The C-O bond formation was testified as the stereoselectivity-determining step, the R-configurational pathway is more energetically favorable than the S-configurational one. The non-covalent interaction (NCI) and atom-in-molecule (AIM) analyses were performed to reveal that the O-H ⋅⋅⋅ O and C-H ⋅⋅⋅ O hydrogen-bond interactions are the key factors for controlling the stereoselectivity. The detailed mechanism and origin of stereoselectivity give useful insights for understanding organocatalytic reactions for asymmetric construction of C-O bond.

Mechanisms of Proteolytic Enzymes and Their Inhibition in QM/MM Studies

Experimental evidence for enzymatic mechanisms is often scarce, and in many cases inadvertently biased by the employed methods. Thus, apparently contradictory model mechanisms can result in decade long discussions about the correct interpretation of data and the true theory behind it. However, often such opposing views turn out to be special cases of a more comprehensive and superior concept. Molecular dynamics (MD) and the more advanced molecular mechanical and quantum mechanical approach (QM/MM) provide a relatively consistent framework to treat enzymatic mechanisms, in particular, the activity of proteolytic enzymes. In line with this, computational chemistry based on experimental structures came up with studies on all major protease classes in recent years; examples of aspartic, metallo-, cysteine, serine, and threonine protease mechanisms are well founded on corresponding standards. In addition, experimental evidence from enzyme kinetics, structural research, and various other methods supports the described calculated mechanisms. One step beyond is the application of this information to the design of new and powerful inhibitors of disease-related enzymes, such as the HIV protease. In this overview, a few examples demonstrate the high potential of the QM/MM approach for sophisticated pharmaceutical compound design and supporting functions in the analysis of biomolecular structures.

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.

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.

Insights into Lewis base-catalyzed chemoselective [3 + 2] and [3 + 4] annulation reactions of MBH carbonates

The chemoselectivity of Lewis base-catalyzed [3 + 2] and [3 + 4] annulation reactions of MBH carbonates can be predicted using FMO overlap analysis.

Origin and stabilization of axial chirality in the construction of naphthyl-C2-indoles: a DFT study

A novel method for predicting the stabilization of axial chirality in the construction of naphthyl-C2-indoles has been suggested for the first time.

Insights into NHC-catalyzed oxidative α-C(sp3)–H activation of aliphatic aldehydes and cascade [2 + 3] cycloaddition with azomethine imines

The NHC catalyst is identified to promote [2 + 3] cycloaddition by avoiding the poor FMO overlap mode in theory.

Insights into Ag(i)-catalyzed addition reactions of amino alcohols to electron-deficient olefins: competing mechanisms, role of catalyst, and origin of chemoselectivity

The competing mechanisms of Ag(i)-catalyzed chemoselective addition reactions of amino alcohols and electron-deficient olefins leading to the O-adduct or N-adduct products were systematically studied with density functional theory methods. Calculations indicate that the AgHMDS/dppe versus AgOAc/dppe catalytic systems can play different roles and thereby generate two different products. The AgHMDS/dppe system works as a Brønsted base to deprotonate the amino alcohol OH to form the Ag–O bond, which leads to formation of the O-adduct. In contrast, the AgOAc/dppe system mainly acts as a Lewis acid to coordinate with O and N atoms of the amino alcohol, but it cannot act as the Brønsted base to further activate the OH group because of its weaker basicity. Therefore, the AgOAc/dppe catalyzed reaction has a mechanism that is similar to the non-catalyzed reaction, and generates the same N-adduct. The obtained insights will be important for rational design of the various kinds of cooperatively catalyzed chemoselective addition reactions, including the use of the less nucleophilic hydroxyl groups of unprotected amino alcohols.

The origin of the chemoselectivities of Ag(i)-catalyzed addition reactions of amino alcohols to olefin has been predicted for the first time.

Theoretical Study on DBU-Catalyzed Insertion of Isatins into Aryl Difluoronitromethyl Ketones: A Case for Predicting Chemoselectivity Using Electrophilic Parr Function

The possible mechanisms of 1,8-diazabicyclo[5.4.0]undec-7-ene (DBU)-catalyzed chemoselective insertion of N-methyl isatin into aryl difluoronitromethyl ketone to synthesize 3,3-disubstituted and 2,2-disubstituted oxindoles have been studied in this work. As revealed by calculated results, the reaction occurs via two competing paths, including α and β carbonyl paths, and each path contains five steps, that is, nucleophilic addition of DBU to ketone, C-C bond cleavage affording difluoromethylnitrate anion and phenylcarbonyl-DBU cation, nucleophilic addition of difluoromethylnitrate anion to carbonyl carbon of N-methyl isatin, acyl transfer process, and dissociation of DBU and product. The computational results suggest that nucleophilic additions on different carbonyl carbons of N-methyl isatin via α and β carbonyl paths would lead to different products in the third step, and β carbonyl path associated with the main product 3,3-disubstituted oxindole is more energetically favorable, which is consistent with the experimental observations. Noteworthy, electrophilic Parr function can be successfully applied for exactly predicting the activity of reaction site and reasonably explaining the chemoselectivity. In addition, the distortion/interaction and noncovalent interaction analyses show that much more hydrogen bond interactions should be responsible for the lower energy of the transition state associated with β carbonyl path. The obtained insights would be valuable for the rational design of efficient organocatalysts for this kind of reactions with high selectivities.

Catalytic Mechanisms for Cofactor-Free Oxidase-Catalyzed Reactions: Reaction Pathways of Uricase-Catalyzed Oxidation and Hydration of Uric Acid

First-principles quantum mechanical/molecular mechanical (QM/MM)-free energy calculations have been performed to uncover how uricase catalyzes metabolic reactions of uric acid (UA), demonstrating that the entire reaction process of UA in uricase consists of two stages-oxidation followed by hydration. The oxidation consists of four steps: (1) chemical transformation from 8-hydroxyxythine to an anionic radical via a proton transfer along with an electron transfer, which is different from the previously proposed electron-transfer mechanism that involves a dianion intermediate (UA^2-) during the catalytic reaction process; (2) proton transfer to the O₂^- anion (radical); (3) diradical recombination to form a peroxo intermediate; (4) dissociation of H₂O₂ to generate the dehydrourate. Hydration, for the most favorable pathway, is initiated by the nucleophilic attack of a water molecule on dehydrourate, along with a concerted proton transfer through residue Thr69 in the catalytic site. According to the calculated free energy profile, the hydration is the rate-determining step, and the corresponding free energy barrier of 16.2 kcal/mol is consistent with that derived from experimental kinetic data, suggesting that the computational insights into the catalytic mechanisms are reasonable. The mechanistic insights not only provide a mechanistic base for future rational design of uricase mutants with improved catalytic activity against uric acid as an improved enzyme therapy, but also are valuable for understanding a variety of other cofactor-free oxidase-catalyzed reactions involving an oxygen molecule.

Fundamental reaction pathway and free energy profile of proteasome inhibition by syringolin A (SylA)

In this study, molecular dynamics (MD) simulations and first-principles quantum mechanical/molecular mechanical free energy (QM/MM-FE) calculations have been performed to uncover the fundamental reaction pathway of proteasome with a representative inhibitor syringolin A (SylA). The calculated results reveal that the reaction process consists of three steps. The first step is a proton transfer process, activating Thr1-O(γ) directly by Thr1-N(z) to form a zwitterionic intermediate. The next step is a nucleophilic attack on the olefin carbon of SylA by the negatively charged Thr1-O(γ) atom. The last step is a proton transfer from Thr1-N(z) to another olefin carbon of SylA to complete the inhibition reaction process. The calculated free energy profile demonstrates that the second step should be the rate-determining step and has the highest free energy barrier of 24.6 kcal mol(-1), which is reasonably close to the activation free energy (∼22.4-23.0 kcal mol(-1)) derived from the available experimental kinetic data. In addition, our computational results indicate that no water molecule can assist the rate-determining step, since the second step is not involved in a proton transfer process. The obtained mechanistic insights should be valuable for understanding the inhibition process of proteasome by SylA and structurally related inhibitors at a molecular level, and thus provide a solid mechanistic base and valuable clues for future rational design of novel, more potent inhibitors of proteasome.

Theoretical investigations toward TMEDA-catalyzed [2 + 4] annulation of allenoate with 1-aza-1,3-diene: mechanism, regioselectivity, and role of the catalyst

A computational study on the reaction between allenoate and 1-aza-1,3-diene catalyzed by TMEDA has been performed using the DFT method.

Reaction Pathway for Cocaine Hydrolase-Catalyzed Hydrolysis of (+)-Cocaine

A recently designed and discovered cocaine hydrolase (CocH), engineered from human butyrylcholinesterase (BChE), has been proven promising as a novel enzyme therapy for treatment of cocaine overdose and addiction because it is highly efficient in catalyzing hydrolysis of naturally occurring (-)-cocaine. It has been known that the CocH also has a high catalytic efficiency against (+)-cocaine, a synthetic enantiomer of cocaine. Reaction pathway and the corresponding free energy profile for the CocH-catalyzed hydrolysis of (+)-cocaine have been determined, in the present study, by performing first-principles pseudobond quantum mechanical/molecular mechanical (QM/MM)-free energy (FE) calculations. Acordingt to the QM/MM-FE results, the catalytic hydrolysis process is initiated by the nucleophilic attack on carbonyl carbon of (-)-cocaine benzoyl ester via hydroxyl oxygen of S198 side chain, and the second reaction step (i.e. dissociation of benzoyl ester) is rate-determining. This finding for CocH-catalyzed hydrolysis of (+)-cocaine is remarkably different from that for the (+)-cocaine hydrolysis catalyzed by bacterial cocaine esterase in which the first reaction step of the deacylation is associated with the highest free energy barrier (~17.9 kcal/mol). The overall free energy barrier (~16.0 kcal/mol) calculated for the acylation stage of CocH-catalyzed hydrolysis of (+)-cocaine is in good agreement with the experimental free energy barrier of ~14.5 kcal/mol derivated from the experimental kinetic data.

Organocatalyzed asymmetric Michael addition by an efficient bifunctional carbohydrate–thiourea hybrid with mechanistic DFT analysis

A series of thiourea based bifunctional organocatalysts having d-glucose as a core scaffold were synthesized and examined as catalysts for the asymmetric Michael addition reaction of aryl/alkyl trans-β-nitrostyrenes over cyclohexanone and other Michael donors having active methylene.

DFT study on the reaction mechanisms and stereoselectivities of NHC-catalyzed [2 + 2] cycloaddition between arylalkylketenes and electron-deficient benzaldehydes

In this paper, two possible mechanisms (mechanisms A and B) on the stereoselective [2 + 2] cycloaddition of aryl(alkyl)ketenes and electron-deficient benzaldehydes catalyzed by N-heterocyclic carbenes (NHCs) have been investigated using density functional theory (DFT). Our calculated results indicate that the favorable mechanism (mechanism A) includes three processes: the first step is the nucleophilic attack on the arylalkylketene by the NHC catalyst to form an intermediate, the second step is the [2 + 2] cycloaddition of the intermediate and benzaldehyde for the formation of a β-lactone, and the last step is the dissociation of the NHC catalyst and the β-lactone. Notably, the [2 + 2] cycloaddition, in which two chiral centers associated with four configurations (SS, RR, SR and RS) are formed, is demonstrated to be both the rate- and stereoselectivity-determining step. Moreover, the reaction pathway associated with the SR configuration is the most favorable pathway and leads to the main product, which is in good agreement with the experimental results. Furthermore, the analysis of global and local reactivity indexes has been performed to explain the role of the NHC catalyst in the [2 + 2] cycloaddition reaction. Therefore, this study will be of great use for the rational design of more efficient catalysts for this kind of cycloaddition.

Mechanistic and stereoselectivity study for the reaction of trifluoropyruvates with arylpropenes catalyzed by a cationic Lewis acid rhodium complex

The mechanism and stereoselectivity of a Lewis acid catalyzed carbonyl–ene reaction of trifluoropyruvates with arylpropenes have been investigated using a DFT method.

Reaction pathways and free energy profiles for cholinesterase-catalyzed hydrolysis of 6-monoacetylmorphine

As the most active metabolite of heroin, 6-monoacetylmorphine (6-MAM) can penetrate into the brain for the rapid onset of heroin effects. The primary enzymes responsible for the metabolism of 6-MAM to the less potent morphine in humans are acetylcholinesterase (AChE) and butyrylcholinesterase (BChE). The detailed reaction pathways for AChE- and BChE-catalyzed hydrolysis of 6-MAM to morphine have been explored, for the first time, in the present study by performing first-principles quantum mechanical/molecular mechanical free energy calculations. It has been demonstrated that the two enzymatic reaction processes follow similar catalytic reaction mechanisms, and the whole catalytic reaction pathway for each enzyme consists of four reaction steps. According to the calculated results, the second reaction step associated with the transition state TS2(a)/TS2(b) should be rate-determining for the AChE/BChE-catalyzed hydrolysis, and the free energy barrier calculated for the AChE-catalyzed hydrolysis (18.3 kcal mol(-1)) is 2.5 kcal mol(-1) lower than that for the BChE-catalyzed hydrolysis (20.8 kcal mol(-1)). The free energy barriers calculated for the AChE- and BChE-catalyzed reactions are in good agreement with the experimentally derived activation free energies (17.5 and 20.7 kcal mol(-1) for the AChE- and BChE-catalyzed reactions, respectively). Further structural analysis reveals that the aromatic residues Phe295 and Phe297 in the acyl pocket of AChE (corresponding to Leu286 and Val288 in BChE) contribute to the lower energy of TS2(a) relative to TS2(b). The obtained structural and mechanistic insights could be valuable for use in future rational design of a novel therapeutic treatment of heroin abuse.

DFT study on the mechanisms and diastereoselectivities of Lewis acid-promoted ketene-alkene [2 + 2] cycloadditions: what is the role of Lewis acid in the ketene and C = X (X = O, CH₂, and NH) [2 + 2] cycloaddition reactions?

The detailed mechanisms and diastereoselectivities of Lewis acid-promoted ketene-alkene [2 + 2] cycloaddition reactions have been studied by density functional theory (DFT). Four possible reaction channels, including two noncatalyzed diastereomeric reaction channels (channels A and B) and two Lewis acid (LA) ethylaluminum dichloride (EtAlCl2) catalyzed diastereomeric reaction channels (channels C and D), have been investigated in this work. The calculated results indicate that channel A (associated with product R-configurational cycloputanone) is more energy favorable than channel B (associated with the other product S-configurational cyclobutanone) under noncatalyzed condition, but channel D leading to S-configurational cyclobutanone is more energy-favorable than channel C, leading to R-configurational cycloputanone under a LA-promoted condition, which is consistent with the experimental results. And Lewis acid can make the energy barrier of ketene-alkene [2 + 2] cycloaddition much lower. In order to explore the role of LA in ketene and C = X (X = O, CH2, and NH) [2 + 2] cycloadditions, we have tracked and compared the interaction modes of frontier molecular orbitals (FMOs) along the intrinsic reaction coordinate (IRC) under the two different conditions. Besides by reducing the energy gap between the FMOs of the reactants, our computational results demonstrate that Lewis acid lowers the energy barrier of the ketene and C = X [2 + 2] cycloadditions by changing the overlap modes of the FMOs, which is remarkably different from the traditional FMO theory. Furthermore, analysis of global reactivity indexes has also been performed to explain the role of LA catalyst in the ketene-alkene [2 + 2] cycloaddition reaction.