Fundamental reaction pathway and free energy profile for inhibition of proteasome by Epoxomicin

The Architecture of the Anbu Complex Reflects an Evolutionary Intermediate at the Origin of the Proteasome System

Proteasomes are self-compartmentalizing proteases that function at the core of the cellular protein degradation machinery in eukaryotes, archaea, and some bacteria. Although their evolutionary history is under debate, it is thought to be linked to that of the bacterial protease HslV and the hypothetical bacterial protease Anbu (ancestral beta subunit). Here, together with an extensive bioinformatic analysis, we present the first biophysical characterization of Anbu. Anbu forms a dodecameric complex with a unique architecture that was only accessible through the combination of X-ray crystallography and small-angle X-ray scattering. While forming continuous helices in crystals and electron microscopy preparations, refinement of sections from the crystal structure against the scattering data revealed a helical open-ring structure in solution, contrasting the ring-shaped structures of proteasome and HslV. Based on this primordial architecture and exhaustive sequence comparisons, we propose that Anbu represents an ancestral precursor at the origin of self-compartmentalization.

•

The crystal structure of the bacterial proteasome homolog Anbu has been solved

•

The dodecameric architecture reveals unique features compared with classical proteasomes

•

Bioinformatic analysis places Anbu at the root of the proteasome family

Computational study on NHC-catalyzed enantioselective and chemoselective fluorination of aliphatic aldehydes

A mechanistic study on NHC-catalyzed enantioselective and chemoselective fluorination of aliphatic aldehydes has been performed for the first time.

Insights into the Competing Mechanisms and Origin of Enantioselectivity for N-Heterocyclic Carbene-Catalyzed Reaction of Aldehyde with Enamide

Hydroacylation reactions and aza-benzoin reactions have attracted considerable attention from experimental chemists. Recently, Wang et al. reported an interesting reaction of N-heterocyclic carbene (NHC)-catalyzed addition of aldehyde to enamide, in which both hydroacylation and aza-benzoin reactions may be involved. Thus, understanding the competing relationship between them is of great interest. Now, density functional theory (DFT) investigation was performed to elucidate this issue. Our results reveal that enamide can tautomerize to its imine isomer with the assistance of HCO3-. The addition of NHC to aldehydes formed Breslow intermediate, which can go through cross-coupling with enamide via hydroacylation reaction or its imine isomer via aza-benzoin reaction. The aza-benzoin reaction requires relatively lower free energy barrier than the hydroacylation reaction. The more polar characteristic of C=N group in the imine isomers, and the more advantageous stereoelectronic effect in the carbon-carbon bond forming transition states in aza-benzoin pathway were identified to determine that the imine isomer can react with the Breslow intermediate more easily. Furthermore, the origin of enantioselectivities for the reaction was explored and reasonably explained by structural analyses on key transition states. The work should provide valuable insights for rational design of switchable NHC-catalyzed hydroacylation and aza-benzoin reactions with high stereoselectivity.

Selective depletion of tumour suppressors Deleted in Colorectal Cancer (DCC) and neogenin by environmental and endogenous serine proteases: linking diet and cancer

BACKGROUND

The related tumour suppressor proteins Deleted in Colorectal Cancer (DCC) and neogenin are absent or weakly expressed in many cancers, whereas their insertion into cells suppresses oncogenic behaviour. Serine proteases influence the initiation and progression of cancers although the mechanisms are unknown.

METHODS

The effects of environmental (bacterial subtilisin) and endogenous mammalian (chymotrypsin) serine proteases were examined on protein expression in fresh, normal tissue and human neuroblastoma and mammary adenocarcinoma lines. Cell proliferation and migration assays (chemoattraction and wound closure) were used to examine cell function. Cells lacking DCC were transfected with an ectopic dcc plasmid.

RESULTS

Subtilisin and chymotrypsin selectively depleted DCC and neogenin from cells at nanomolar concentrations without affecting related proteins. Cells showed reduced adherence and increased migration, but after washing they re-attached within 24 h, with recovery of protein expression. These effects are induced by chymotryptic activity as they are prevented by chymostatin and the soybean Bowman-Birk inhibitor typical of many plant protease inhibitors.

CONCLUSIONS

Bacillus subtilis, which secretes subtilisin is widely present in soil, the environment and the intestinal contents, while subtilisin itself is used in meat processing, animal feed probiotics and many household cleaning agents. With chymotrypsin present in chyme, blood and tissues, these proteases may contribute to cancer development by depleting DCC and neogenin. Blocking their activity by Bowman-Birk inhibitors may explain the protective effects of a plant diet. Our findings identify a potential non-genetic contribution to cancer cell behaviour which may explain both the association of processed meats and other factors with cancer incidence and the protection afforded by plant-rich diets, with significant implications for cancer prevention.

Computational Approaches for the Discovery of Human Proteasome Inhibitors: An Overview

Proteasome emerged as an important target in recent pharmacological research due to its pivotal role in degrading proteins in the cytoplasm and nucleus of eukaryotic cells, regulating a wide variety of cellular pathways, including cell growth and proliferation, apoptosis, DNA repair, transcription, immune response, and signaling processes. The last two decades witnessed intensive efforts to discover 20S proteasome inhibitors with significant chemical diversity and efficacy. To date, the US FDA approved to market three proteasome inhibitors: bortezomib, carfilzomib, and ixazomib. However new, safer and more efficient drugs are still required. Computer-aided drug discovery has long being used in drug discovery campaigns targeting the human proteasome. The aim of this review is to illustrate selected in silico methods like homology modeling, molecular docking, pharmacophore modeling, virtual screening, and combined methods that have been used in proteasome inhibitors discovery. Applications of these methods to proteasome inhibitors discovery will also be presented and discussed to raise improvements in this particular field.

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.

DFT perspective toward [3 + 2] annulation reaction of enals with α-ketoamides through NHC and Brønsted acid cooperative catalysis: mechanism, stereoselectivity, and role of NHC

The DFT perspective toward the [3 + 2] annulation reaction through NHC and Brønsted acid cooperative catalysis has been investigated.

A DFT study on PBu3-catalyzed intramolecular cyclizations of N-allylic substituted α-amino nitriles for the formation of functionalized pyrrolidines: mechanisms, selectivities, and the role of catalysts

The detailed mechanisms and stereoselectivities of PBu₃-catalyzed intramolecular cyclizations for the formation of functionalized pyrrolidines have been investigated using a DFT method.

A DFT study on NHC-catalyzed intramolecular aldehyde–ketone crossed-benzoin reaction: mechanism, regioselectivity, stereoselectivity, and role of NHC

A computational study on NHC-catalyzed intramolecular aldehyde–ketone crossed-benzoin reaction has been performed using a DFT method.

A computational study on the N-heterocyclic carbene-catalyzed Csp2–Csp3 bond activation/[4+2] cycloaddition cascade reaction of cyclobutenones with imines: a new application of the conservation principle of molecular orbital symmetry

A new application of the conservation principle of molecular orbital symmetry has been suggested for disclosing the role of NHC.

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.

Water sandwiched by a pair of aromatic rings in a proton-conducting metal–organic framework

Water–benzene interactions are investigated in a metal–organic framework incorporating strong H-bond nets which contribute to distinct proton conduction.

Mechanistic insights into cobalt(ii/iii)-catalyzed C-H oxidation: a combined theoretical and experimental study

Cobalt-mediated C-H functionalization has been the subject of extensive interest in synthetic chemistry, but the mechanisms of many of these reactions (such as the cobalt-catalyzed C-H oxidation) are poorly understood. In this paper, possible mechanisms including single electron transfer (SET) and the concerted metalation-deprotonation (CMD) pathways of the CoII/CoIII-catalyzed alkoxylation of C(sp2)-H bonds have been investigated for the first time using the DFT method. CoII(OAc)2 has been employed as an efficient catalyst in our previous experimental study, but the calculated results unexpectedly indicated that the intermolecular SET pathway with CoIII as the actual catalyst might be the most favorable pathway. To support this theoretical prediction, we have explored a series of Cp*CoIII(CO)I2 catalyzed C(sp2)-H bond alkoxylations, extending the application of cobalt-catalyzed functionalization of C-H bonds. Furthermore, kinetic isotope effect (KIE) data, electron paramagnetic resonance (EPR) data, and TEMPO inhibition experiments also support the SET mechanism in both the Co-catalyzed alkoxylation reactions. Thus, this work should support an understanding of the possible mechanisms of the CoII/CoIII-catalyzed C(sp2)-H functionalization, and also provide an example of the rational design of novel catalytic reactions guided by theoretical calculations.

David and Goliath: chemical perturbation of eukaryotes by bacteria

Environmental microbes produce biologically active small molecules that have been mined extensively as antibiotics and a smaller number of drugs that act on eukaryotic cells. It is known that there are additional bioactives to be discovered from this source. While the discovery of new antibiotics is challenged by the frequent discovery of known compounds, we contend that the eukaryote-active compounds may be less saturated. Indeed, despite there being far fewer eukaryotic-active natural products these molecules interact with a far richer diversity of molecular and cellular targets.

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.

Proaporphine and aporphine alkaloids with acetylcholinesterase inhibitory activity from Stephania epigaea

An unusual proaporphine alkaloid bearing an isopropanenitrile group at isoquinoline nitrogen, named epiganine A (1) and a new aporphine alkaloid, epiganine B (2), together with eight known alkaloids, pronuciferine (3), dehydrodicentrine (4), romerine (5), romeline (6), N-methylcalycinine (7), phanostenine (8), dicentrine (9), and N-methyllaurotetanine (10), were isolated from the roots of Stephania epigaea. The absolute configuration of 1 was determined by calculating electronic circular dichroism (ECD) and comparing with experimental data. Compounds 2 and 4 showed strong acetylcholinesterase (AChE) inhibitory effects with the IC50 values of 4.36 and 2.98μM, respectively. Compounds 5-9 also exhibited potent AChE inhibitory activities.

Mechanisms and stereoselectivities of the Rh(i)-catalyzed carbenoid carbon insertion reaction of benzocyclobutenol with diazoester

Mechanisms and stereoselectivities of a Rh(i)-catalyzed carbenoid carbon insertion reaction of benzocyclobutenol with diazoester have been investigated using the DFT method.

A quantum mechanical study of the mechanism and stereoselectivity of the N-heterocyclic carbene catalyzed [4 + 2] annulation reaction of enals with azodicarboxylates

The mechanism and stereoselectivity of the NHC-catalyzed [4 + 2] annulation reaction of enals with azodicarboxylates have been investigated using the 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.

A new insight into the 5-carboxycytosine and 5-formylcytosine under typical bisulfite conditions: a deamination mechanism study

5-Methylcytosine (5-MeCyt) can be converted to 5-hydroxymethylcytosine (5-hmCyt) in mammalian DNA by the ten-eleven translocation enzymes. The conventional bisulfite sequencing cannot discriminate 5-hmCyt from 5-MeCyt, whereas the oxidation products of 5-hmCyt, 5-carboxycytosine (5-caCyt) and 5-formylcytosine (5-fCyt) enable them to be identified in bisulfite sequencing. This mechanism likely involves the decarboxylation of 5-caCyt and deformylation of 5-fCyt to cytosine (Cyt) before deamination. Another possibility could be a direct bisulfite-induced deamination reaction followed by decarboxylation and deformylation. Here the HSO3(-)-induced direct hydrolytic deamination of 5-caCytN3(+)-SO3(-) (paths A and B) and 5-O(+)fCytN3(+)-SO3(-) (paths C and D) has been explored at the MP2/6-311++G(3df,3pd)//B3LYP/6-311++G(d,p) level. The activation free energy (ΔG(s≠) = 54.16 kJ mol(-1)) of the direct hydrolytic deamination of 5-caCytN3(+)-SO3(-) path A is much lower than the ΔG(s≠) of CytN3(+)-SO3(-) (100.91 kJ mol(-1)) under bisulfite conditions, implying that 5-caCyt may firstly involve a process of deamination. Meanwhile, the ΔG(s≠) (103.84 kJ mol(-1)) of the HSO3(-)-induced direct hydrolytic deamination of 5-O(+)fCytN3(+)-SO3(-) path C is in close proximity to our previous theoretical data for CytN3(+)-SO3(-), indicating that the deamination of 5-fCyt is also likely to occur in the presence of bisulfite. Meanwhile, the HSO3(-)-induced direct hydrolytic deamination of 5-caCytN3(+)-SO3(-) path A and 5-O(+)fCytN3(+)-SO3(-) path C is represented and has been further explored in the presence of one and two water molecules. The results show that both in the gas and aqueous phases, the participation of one and two water molecules makes the HSO3(-)-induced direct hydrolytic deamination of 5-caCytN3(+)-SO3(-) path A unfavorable, whereas the contribution of one and two water molecules facilitates the HSO3(-)-induced direct hydrolytic deamination of 5-O(+)fCytN3(+)-SO3(-) path C.

Is the contribution of cis and trans protonated 5-methylcytosine-SO3(-) isomers equal in the conversion to thymine-SO3(-) under bisulfite conditions? A theoretical perspective

Cytosine (Cyt) can be converted to 5-methylcytosine (5-MeCyt) in CpG sequences of DNA. Conventional bisulfite sequencing can discriminate Cyt from 5-MeCyt, however inappropriate conversion of 5-MeCyt to thymine and a failure to convert Cyt to uracil always occur when Cyt and 5-MeCyt are treated with bisulfite, which would lead to erroneous estimates of DNA methylation densities. Here, the direct hydrolytic deamination of cis (paths A-C) and trans (paths A'-C') 5-MeCytN3(+)-SO3(-) isomers with bisulfite have been explored at the MP2/6-311++G(3df,3pd)//B3LYP/6-311++G(d,p) level. The activation free energies (ΔG(s-a≠)) of the cis and trans 5-MeCytN3(+)-SO3(-) isomers' paths exhibit no obvious differences, implying both isomers may make an equal contribution to the hydrolytic deamination of 5-MeCyt under bisulfite conditions. It is greatly expected that these results could aid experimental scientists to explore new methods to avoid the formation of the deaminated reactants (5-MeCytN3(+)-SO3(-)). Meanwhile, the HSO3(-)-induced direct hydrolytic deamination of cis and trans 5-MeCytN3(+)-SO3(-) isomers is represented by paths A and A', respectively, and has been further explored in the presence of two water molecules. It was found that the contribution of two water molecules renders the HSO3(-)-induced direct hydrolytic deamination of cis and trans 5-MeCytN3(+)-SO3(-) isomers by paths A and A' favourable. In addition, the ΔG(s-a≠) values (85.74-85.34 kJ mol(-1)) of the rate-limiting steps of the two water-mediated paths A and A' are very close to that of the theoretical value for CytN3(+)-SO3(-) (88.18 kJ mol(-1)), implying that the free barrier gap between Cyt and 5-MeCyt is very small under bisulfite conditions. This further suggests that bisulfite sequencing technology may be easily influenced by the external environment.

Enzyme inhibition by hydroamination: design and mechanism of a hybrid carmaphycin-syringolin enone proteasome inhibitor

Hydroamination reactions involving the addition of an amine to an inactivated alkene are entropically prohibited and require strong chemical catalysts. While this synthetic process is efficient at generating substituted amines, there is no equivalent in small molecule-mediated enzyme inhibition. We report an unusual mechanism of proteasome inhibition that involves a hydroamination reaction of alkene derivatives of the epoxyketone natural product carmaphycin. We show that the carmaphycin enone first forms a hemiketal intermediate with the catalytic Thr1 residue of the proteasome before cyclization by an unanticipated intramolecular alkene hydroamination reaction, resulting in a stable six-membered morpholine ring. The carmaphycin enone electrophile, which does not undergo a 1,4-Michael addition as previously observed with vinyl sulfone and α,β-unsaturated amide-based inhibitors, is partially reversible and gives insight into the design of proteasome inhibitors for cancer chemotherapy.

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.

Selective immunoproteasome inhibitors with non-peptide scaffolds identified from structure-based virtual screening

As a major component of the crucial nonlysosomal protein degradation pathway in the cells, the proteasome has been implicated in many diseases such as Alzheimer's disease, Huntington's disease, inflammatory bowel diseases, autoimmune diseases, multiple myeloma (MM) and other cancers. There are two main proteasome subtypes: the constitutive proteasome which is expressed in all eukaryotic cells and the immunoproteasome which is expressed in immune cells and can be induced in other cell types. Majority of currently available proteasome inhibitors are peptide backbone-based, having short half-lives in the body. It is highly desirable to identify novel, immunoproteasome-selective inhibitors with non-peptide scaffolds for development of novel therapeutics. Through combined virtual screening and experimental studies targeting the immunoproteasome, we have identified a set of novel immunoproteasome inhibitors with diverse non-peptide scaffolds. Some of the identified inhibitors have significant selectivity for the immunoproteasome over the constitutive proteasome. Unlike most of the currently available proteasome inhibitors, these new inhibitors lacking electrophilic pharmacophores are not expected to form a covalent bond with proteasome after the binding. These non-peptide scaffolds may provide a new platform for future rational drug design and discovery targeting the immunoproteasome.

A DFT study on the reaction mechanism of dimerization of methyl methacrylate catalyzed by N-heterocyclic carbene

The novel proton transfers assisted by methyl methacrylate lower the energy barrier of the dimerization.

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

Proteasome is the major component of the crucial non-lysosomal protein degradation pathway in the cells, but the detailed reaction pathway is unclear. In this study, first-principles quantum mechanical/molecular mechanical free energy calculations have been performed to explore, for the first time, possible reaction pathways for proteasomal proteolysis/hydrolysis of a representative peptide, succinyl-leucyl-leucyl-valyl-tyrosyl-7-amino-4-methylcoumarin (Suc-LLVY-AMC). The computational results reveal that the most favorable reaction pathway consists of six steps. The first is a water-assisted proton transfer within proteasome, activating Thr1-O(γ). The second is a nucleophilic attack on the carbonyl carbon of a Tyr residue of substrate by the negatively charged Thr1-O(γ), followed by the dissociation of the amine AMC (third step). The fourth step is a nucleophilic attack on the carbonyl carbon of the Tyr residue of substrate by a water molecule, accompanied by a proton transfer from the water molecule to Thr1-N(z). Then, Suc-LLVY is dissociated (fifth step), and Thr1 is regenerated via a direct proton transfer from Thr1-N(z) to Thr1-O(γ). According to the calculated energetic results, the overall reaction energy barrier of the proteasomal hydrolysis is associated with the transition state (TS3(b)) for the third step involving a water-assisted proton transfer. The determined most favorable reaction pathway and the rate-determining step have provided a reasonable interpretation of the reported experimental observations concerning the substituent and isotopic effects on the kinetics. The calculated overall free energy barrier of 18.2 kcal/mol is close to the experimentally derived activation free energy of ∼18.3-19.4 kcal/mol, suggesting that the computational results are reasonable.

Fundamental reaction pathway and free energy profile for butyrylcholinesterase-catalyzed hydrolysis of heroin

The pharmacological function of heroin requires an activation process that transforms heroin into 6-monoacetylmorphine (6-MAM), which is the most active form. The primary enzyme responsible for this activation process in human plasma is butyrylcholinesterase (BChE). The detailed reaction pathway of the activation process via BChE-catalyzed hydrolysis has been explored computationally, for the first time, in this study via molecular dynamics simulation and first-principles quantum mechanical/molecular mechanical free energy calculations. It has been demonstrated that the whole reaction process includes acylation and deacylation stages. The acylation consists of two reaction steps, i.e., the nucleophilic attack on the carbonyl carbon of the 3-acetyl group of heroin by the hydroxyl oxygen of the Ser198 side chain and the dissociation of 6-MAM. The deacylation also consists of two reaction steps, i.e., the nucleophilic attack on the carbonyl carbon of the acyl-enzyme intermediate by a water molecule and the dissociation of the acetic acid from Ser198. The calculated free energy profile reveals that the second transition state (TS2) should be rate-determining. The structural analysis reveals that the oxyanion hole of BChE plays an important role in the stabilization of rate-determining TS2. The free energy barrier (15.9 ± 0.2 or 16.1 ± 0.2 kcal/mol) calculated for the rate-determining step is in good agreement with the experimentally derived activation free energy (~16.2 kcal/mol), suggesting that the mechanistic insights obtained from this computational study are reliable. The obtained structural and mechanistic insights could be valuable for use in the future rational design of a novel therapeutic treatment of heroin abuse.

Reaction pathway and free energy profile for papain-catalyzed hydrolysis of N-acetyl-Phe-Gly 4-nitroanilide

Possible reaction pathways for papain-catalyzed hydrolysis of N-acetyl-Phe-Gly 4-nitroanilide (APGNA) have been studied by performing pseudobond first-principles quantum mechanical/molecular mechanical-free energy (QM/MM-FE) calculations. The whole hydrolysis process includes two stages: acylation and deacylation. For the acylation stage of the catalytic reaction, we have explored three possible paths (A, B, and C) and the corresponding free energy profiles along the reaction coordinates. It has been demonstrated that the most favorable reaction path in this stage is path B consisting of two reaction steps: the first step is a proton transfer to form a zwitterionic form (i.e., Cys-S⁻/His-H⁺ ion-pair), and the second step is the nucleophilic attack on the carboxyl carbon of the substrate accompanied by the dissociation of 4-nitroanilide. The deacylation stage includes the nucleophilic attack of a water molecule on the carboxyl carbon of the substrate and dissociation between the carboxyl carbon of the substrate and the sulfhydryl sulfur of Cys25 side chain. The free energy barriers calculated for the acylation and deacylation stages are 20.0 and 10.7 kcal/mol, respectively. Thus, the acylation is rate-limiting. The overall free energy barrier calculated for papain-catalyzed hydrolysis of APGNA is 20.0 kcal/mol, which is reasonably close to the experimentally derived activation free energy of 17.9 kcal/mol.

Synthesis and pharmacology of proteasome inhibitors

Shortly after the discovery of the proteasome it was proposed that inhibitors could stabilize proteins which ultimately would trigger apoptosis in tumor cells. The essential questions were whether small molecules would be able to inhibit the proteasome without generating prohibitive side effects and how one would derive these compounds. Fortunately, "Mother Nature" has generated a wide variety of natural products that provide distinct selectivities and specificities. The chemical synthesis of these natural products finally provided access to analogues and optimized drugs of which two different classes have been approved for the treatment of malignancies. Despite these achievements, additional lead structures derived from nature are under investigation and will be discussed with regard to their biological potential and chemical challenges.

Catalytic mechanism of cytochrome P450 for N-methylhydroxylation of nicotine: reaction pathways and regioselectivity of the enzymatic nicotine oxidation

The fundamental reaction mechanism of cytochrome P450 2A6 (CYP2A6)-catalyzed N-methylhydroxylation of (S)-(-)-nicotine and the free energy profile have been studied by performing pseudobond first-principles quantum mechanical/molecular mechanical (QM/MM) reaction-coordinate calculations. In the CYP2A6-(S)-(-)-nicotine binding structures that allow for 5'-hydroxylation, the N-methyl group is also sufficiently close to the oxygen of Cpd I for the N-methylhydroxylation reaction to occur. It has been demonstrated that the CYP2A6-catalyzed N-methylhydroxylation reaction is a concerted process involving a hydrogen-transfer transition state on both the quartet and the doublet states. The N-methylhydroxylation reaction proceeds mainly in the doublet state, since the free energy barriers on the doublet state are lower than the corresponding ones on the quartet state. The calculated free energy barriers indicate that (S)-(-)-nicotine oxidation catalyzed by CYP2A6 proceeds with a high regioselective abstraction of the hydrogen at the 5'-position, rather than the hydrogen at the N-methyl group. The predicted regioselectivity of 93% is in agreement with the most recent experimentally reported regioselectivity of 95%. The binding mode of (S)-(-)-nicotine in the active site of CYP2A6 is an important determinant for the stereoselectivity of nicotine (S)-(-)-oxidation, whereas the regioselectivity of (S)-(-)-nicotine oxidation is determined mainly by the free energy barrier difference between the 5'-hydroxylation and N-methylhydroxylation reactions.