Deacylation Stage Research Articles

A rationale for the unlike potency of avibactam and ETX2514 against OXA-24 β-lactamase

Diazabicyclooctanone inhibitors such as ETX2514 and avibactam have shown enhanced inhibitory performance to fight the antibiotic resistance developed by pathogens. However, avibactam is ineffective against Acinetobacter baumannii infections, unlike ETX2514. The molecular basis for this difference has not been tackled from a molecular approach, precluding the knowledge of relevant information. In this article, the mechanisms involved in the inhibition of OXA-24 by ETX2514 and avibactam are studied theoretically by hybrid QM/MM calculations. The results show that both inhibitors share the same inhibition mechanisms, comprising acylation a deacylation stages. The involved mechanisms include the same number of steps, transition states and intermediates; although they differ in the involved activation barriers. This difference accounts for the dissimilar inhibitory ability of both inhibitors. The molecular reason for this is the endocyclic double bond in the piperidine ring of ETX2514 increasing the ring strain and chemical reactivity on the N6 and C7 atoms, besides the methyl substituent, which enhance the hydrophobic character of the ring. Furthermore, Lys218 and the carboxylated Lys84 of ETX2514, play a crucial role in the mechanism by coordinating their protonation states in an on/off (protonated/deprotonated) manner, favoring the proton transference between the residues and the inhibitor.

Unrevealing the Proteolytic Activity of RgpB Gingipain from Computational Simulations.

Alzheimer's disease represents one of the greatest medical concerns for today's population and health services. Its multifactorial inherent nature represents a challenge for its treatment and requires the development of a broad spectrum of drugs. Recently, the cysteine protease gingipain RgpB has been related to neurodegenerative diseases, including Alzheimer's disease, and its inhibition appears to be a promising neuroprotective strategy. Given these features, a computational study that integrates molecular dynamics (MD) simulations with classical and hybrid quantum mechanics/molecular mechanics (QM/MM) potentials was carried out to unravel the atomistic details of RgpB activity. First, a preliminary study based on principal component analysis (PCA), determined the protonation state of the Cys/His catalytic dyad, as well as the crucial role of a flexible loop that favors reactive interactions of the catalytic residues and the peptide in the precatalytic state in its closed conformation. Then, different mechanisms were explored by means of QM/MM MD simulations. The most favorable mechanism consists of two stages. First is an acylation stage that takes place in two steps where, initially, the sulfur atom of the C244 residue attacks the carbonylic carbon of the peptide and the proton of the C244 residue is transferred to the amino group of the peptide in a concerted manner. Subsequently, the peptide bond is broken, and a fragment of the peptide is released. After that, the deacylation stage takes place in a single step where a water molecule attacks the carbonylic carbon of the peptide and a proton of the water is transferred to the C244 residue. The free energy barrier of the rate limiting step is in very good agreement with available experimental data. The mechanism exhibits an unusual role of H211 residue compared with other cysteine proteases but a crucial role of the peptide in triggering the catalysis. Notably, the atomic and energetic particularities found represent a significant contribution to the comprehension of the reaction mechanism and a great opportunity for the design of efficient inhibitors of gingipain RgpB.

A theoretical approach for the acylation/deacylation mechanisms of avibactam in the reversible inhibition of KPC-2.

Klebsiella pneumoniae carbapenemase (KPC-2) is the most commonly encountered class A β-lactamase variant worldwide, which confer high-level resistance to most available antibiotics. In this article we address the issue by a combined approach involving molecular dynamics simulations and hybrid quantum mechanics/molecular mechanics calculations. The study contributes to improve the understanding, at molecular level, of the acylation and deacylation stages of avibactam involved in the inhibition of KPC-2. The results show that both mechanisms, acylation and deacylation, the reaction occur via the formation of a tetrahedral intermediate. The formation of this intermediate corresponds to the rate limiting stage. The activation barriers are 19.5 kcal/mol and 23.0 kcal/mol for the acylation and deacylation stages, respectively. The associated rate constants calculated, using the Eyring equation, are 1.2 × 10−1 and 3.9 × 10−4 (s−1). These values allow estimating a value of 3.3 × 10−3 for the inhibition constant, in good agreement with the experimental value.Supplementary InformationThe online version contains supplementary material available at 10.1007/s10822-021-00408-3.

Molecular mechanism of chromogenic substrate hydrolysis in the active site of human carboxylesterase-1

Human carboxylesterases are involved in the protective processes of detoxification during the hydrolytic metabolism of xenobiotics. Knowledge of the molecular mechanisms of substrates hydrolysis in the enzymes active site is necessary for the rational drug design. In this work, the molecular mechanism of the hydrolysis reaction of para-nitrophenyl acetate in the active site of human carboxylesterase was determined using modern methods of molecular modeling. According to the combined method of quantum mechanics/molecular mechanics calculations, the chemical reaction occurs within four elementary steps, including two steps of the acylation stage, and two steps of the deacylation stage. All elementary steps have low energy barriers, with the gradual lowering of the intermediate energies that stimulates reaction in the forward direction. The molecular docking was used to estimate the binding constants of the enzyme-substrate complex and the dissociation constant of enzyme-product complexes. The effective kinetic parameters of the enzymatic hydrolysis in the active site of carboxylesterase are determined by numerical solution of the differential kinetic equations.

Catalytic Hydrolysis Mechanism of Cocaine by Human Carboxylesterase 1: An Orthoester Intermediate Slows Down the Reaction.

Human carboxylesterase 1 (hCES1) is a major carboxylesterase in the human body and plays important roles in the metabolism of a wide variety of substances, including lipids and drugs, and therefore is attracting more and more attention from areas including lipid metabolism, pharmacokinetics, drug–drug interactions, and prodrug activation. In this work, we studied the catalytic hydrolysis mechanism of hCES1 by the quantum mechanics computation method, using cocaine as a model substrate. Our results support the four-step theory of the esterase catalytic hydrolysis mechanism, in which both the acylation stage and the deacylation stage include two transition states and a tetrahedral intermediate. The roles and cooperation of the catalytic triad, S221, H468, and E354, were also analyzed in this study. Moreover, orthoester intermediates were found in hCES1-catalyzed cocaine hydrolysis reaction, which significantly elevate the free energy barrier and slow down the reaction. Based on this finding, we propose that hCES1 substrates with β-aminocarboxylester structure might form orthoester intermediates in hCES1-catalyzed hydrolysis, and therefore prolong their in vivo half-life. Thus, this study helps to clarify the catalytic mechanism of hCES1 and elucidates important details of its catalytic process, and furthermore, provides important insights into the metabolism of hCES1 substrates and drug designing.

A water-assisted nucleophilic mechanism utilized by BphD, the meta-cleavage product hydrolase in biphenyl degradation

As members of the α/β-hydrolase superfamily, Meta-cleavage product (MCP) hydrolases generally utilize a Ser-His-Asp catalytic triad to hydrolyze the cleavage of CC bond during the aerobic catabolism of aromatic compounds by bacteria. BphD is one kind of MCP hydrolase that catalyzes the hydrolysis of 2-hydroxy-6-oxo-6-phenylhexa-2,4-dienoic acid (HOPDA) to 2-hydroxypenta-2,4-dienoic acid (HPD) and benzoate. In this article, a combined quantum mechanics and molecule mechanics (QM/MM) approach has been employed to explore the reaction mechanism of BphD from Burkholderia xenovorans LB400. On the basis of the recently resolved crystal structures, three computational models have been constructed. Our calculation results reveal that BphD utilizes a water-assisted nucleophilic mechanism, which contains acylation and deacylation stages. In acylation reaction, an active site water molecule assists the proton transfer from Ser112 to the carbanion intermediate (substrate) by forming hydrogen bonds with Ser112 and His265, and this proton transfer is in concert with the nucleophilic attack of deprotonated Ser112 on the C6-carbonyl of substrate to form the acylated intermediate. In deacylation, the Asp237-His265 dyad acts as a general base to activate the hydrolytic water, whose nucleophilic attack leads to the collapses of acyl-enzyme intermediate. The acylation and deacylation process correspond to the highest energy barriers of 21.0 and 23.9kcal/mol, respectively. During the catalytic reaction, the active site water and Asp237-His265 dyad play an important role for each elementary steps.

Computational Study of the Catalytic Mechanism of the Cruzain Cysteine Protease

Cysteine proteases of the papain family are involved in many diseases, making them attractive targets for developing drugs. In this paper, different catalytic mechanisms of cruzain cysteine protease have been studied on the basis of molecular dynamics simulations within hybrid quantum mechanics/molecular mechanics potentials. The obtained free energy surfaces have allowed characterizing every single step along the mechanisms. The results confirm that the full process can be divided into an acylation and a deacylation stage, but important differences with respect to previous studies can be deduced from our calculations. Thus, our calculations suggest that the acylation stage takes place in a stepwise mechanism where a proton from a conserved His159 is transferred first to the N1 atom of the peptide and, after a transient intermediate is located, the Cys25 attacks the carbonyl carbon atom. The stabilization of the activated Cys25 is achieved by an effect of the local environment through interactions with re...

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.

Reaction Pathway and Free Energy Profiles for Butyrylcholinesterase-Catalyzed Hydrolysis of Acetylthiocholine

The catalytic mechanism for butyrylcholineserase (BChE)-catalyzed hydrolysis of acetylthiocholine (ATCh) has been studied by performing pseudobond first-principles quantum mechanical/molecular mechanical-free energy (QM/MM-FE) calculations on both acylation and deacylation of BChE. Additional quantum mechanical (QM) calculations have been carried out, along with the QM/MM-FE calculations, to understand the known substrate activation effect on the enzymatic hydrolysis of ATCh. It has been shown that the acylation of BChE with ATCh consists of two reaction steps including the nucleophilic attack on the carbonyl carbon of ATCh and the dissociation of thiocholine ester. The deacylation stage includes nucleophilic attack of a water molecule on the carboxyl carbon of substrate and dissociation between the carboxyl carbon of substrate and hydroxyl oxygen of Ser198 side chain. QM/MM-FE calculation results reveal that the acylation of BChE is rate-determining. It has also been demonstrated that an additional substrate molecule binding to the peripheral anionic site (PAS) of BChE is responsible for the substrate activation effect. In the presence of this additional substrate molecule at PAS, the calculated free energy barrier for the acylation stage (rate-determining step) is decreased by ~1.7 kcal/mol. All of our computational predictions are consistent with available experimental kinetic data. The overall free energy barriers calculated for BChE-catalyzed hydrolysis of ATCh at regular hydrolysis phase and substrate activation phase are ~13.6 and ~11.9 kcal/mol, respectively, which are in reasonable agreement with the corresponding experimentally derived activation free energies of 14.0 kcal/mol (for regular hydrolysis phase) and 13.5 kcal/mol (for substrate activation phase).

Reaction Pathway and Free Energy Profile for Butyrylcholinesterase-Catalyzed Hydrolysis of Acetylcholine

A catalytic mechanism for the butyrylcholinesterase (BChE)-catalyzed hydrolysis of acetylcholine (ACh) has been studied by performing pseudobond first-principles quantum mechanical/molecular mechanical-free energy calculations on both acylation and deacylation of BChE. It has been shown that the acylation with ACh includes two reaction steps, including nucleophilic attack on the carbonyl carbon of ACh and dissociation of choline ester. The deacylation stage includes nucleophilic attack of a water molecule on the carboxyl carbon of the substrate and dissociation between the carboxyl carbon of the substrate and the hydroxyl oxygen of the Ser198 side chain. Notably, despite the fact that acetylcholinesterase (AChE) and BChE are very similar enzymes, the acylation of BChE with ACh is rate-determining, which is remarkably different from the AChE-catalyzed hydrolysis of ACh, in which the deacylation is rate-determining. The computational prediction is consistent with available experimental kinetic data. The overall free energy barrier calculated for BChE-catalyzed hydrolysis of ACh is 13.8 kcal/mol, which is in good agreement with the experimentally derived activation free energy of 13.3 kcal/mol.

Characterization of a complete cycle of acetylcholinesterase catalysis by ab initio QM/MM modeling

The reaction mechanism of acetylcholine hydrolysis by acetylcholinesterase, including both acylation and deacylation stages from the enzyme-substrate (ES) to the enzyme-product (EP) molecular complexes, is examined by using an ab initio type quantum mechanical - molecular mechanical (QM/MM) approach. The density functional theory PBE0/aug-6-31+G* method for a fairly large quantum part trapped inside the native protein environment, and the AMBER force field parameters in the molecular mechanical part are employed in computations. All reaction steps, including the formation of the first tetrahedral intermediate (TI1), the acylenzyme (EA) complex, the second tetrahedral intermediate (TI2), and the EP complex, are modeled at the same theoretical level. In agreement with the experimental rate constants, the estimated activation energy barrier of the deacylation stage is slightly higher than that for the acylation phase. The critical role of the non-triad Glu202 amino acid residue in orienting lytic water molecule and in stabilizing the second tetrahedral intermediate at the deacylation stage of the enzymatic process is demonstrated.

A rate-limiting step of enteropeptidase hydrolysis

It has been shown for the first time that deacylation is the rate-limiting stage in the enteropeptidase-catalyzed hydrolysis of highly efficient oligopeptide substrates containing four Asp residues in positions P2-P5. On the other hand, the rate-limiting stage in the hydrolysis of low-efficiency peptide substrates containing less than four Asp or Glu residues in positions P2-P5 is acylation, as has previously been suggested for all amide and peptide substrates of serine proteases on the basis of the classic studies by Bender et al. The method of introduction of an additional nucleophile or another effector that selectively affects the deacylation stage was used to determine the rate-limiting stage in the enteropeptidase hydrolysis of Nalpha-benzyloxycarbonyl-L-lysine thiobenzyl ester, the highly efficient amide substrate GlyAsp4-Lys beta-naphthyl amide, and the low-efficiency peptide substrate VLSAADK-GNVKAAWG (where a hyphen denotes the hydrolysis site). The English version of the paper: Russian Journal of Bioorganic Chemistry, 2008, vol. 34, no. 2; see also http://www.maik.ru.

Fundamental reaction mechanism for cocaine hydrolysis in human butyrylcholinesterase.

Butyrylcholinesterase (BChE)-cocaine binding and the fundamental pathway for BChE-catalyzed hydrolysis of cocaine have been studied by molecular modeling, molecular dynamics (MD) simulations, and ab initio calculations. Modeling and simulations indicate that the structures of the prereactive BChE/substrate complexes for (-)-cocaine and (+)-cocaine are all similar to that of the corresponding prereactive BChE/butyrylcholine (BCh) complex. The overall binding of BChE with (-)-cocaine and (+)-cocaine is also similar to that proposed with butyrylthiocholine and succinyldithiocholine, i.e., (-)- or (+)-cocaine first slides down the substrate-binding gorge to bind to Trp-82 and stands vertically in the gorge between Asp-70 and Trp-82 (nonprereactive complex) and then rotates to a position in the catalytic site within a favorable distance for nucleophilic attack and hydrolysis by Ser-198 (prereactive complex). In the prereactive complex, cocaine lies horizontally at the bottom of the gorge. The fundamental catalytic hydrolysis pathway, consisting of acylation and deacylation stages similar to those for ester hydrolysis by other serine hydrolases, was proposed on the basis of the simulated prereactive complex and confirmed theoretically by ab initio reaction coordinate calculations. Both the acylation and deacylation follow a double-proton-transfer mechanism. The calculated energetic results show that within the chemical reaction process the highest energy barrier and Gibbs free energy barrier are all associated with the first step of deacylation. The calculated ratio of the rate constant (k(cat)) for the catalytic hydrolysis to that (k(0)) for the spontaneous hydrolysis is approximately 9.0 x 10(7). The estimated k(cat)/k(0) value of approximately 9.0 x 10(7) is in excellent agreement with the experimentally derived k(cat)/k(0) value of approximately 7.2 x 10(7) for (+)-cocaine, whereas it is approximately 2000 times larger than the experimentally derived k(cat)/k(0) value of approximately 4.4 x 10(4) for (-)-cocaine. All of the results suggest that the rate-determining step of the BChE-catalyzed hydrolysis of (+)-cocaine is the first step of deacylation, whereas for (-)-cocaine the change from the nonprereactive complex to the prereactive complex is rate-determining and has a Gibbs free energy barrier higher than that for the first step of deacylation by approximately 4 kcal/mol. A further analysis of the structural changes from the nonprereactive complex to the prereactive complex reveals specific amino acid residues hindering the structural changes, providing initial clues for the rational design of BChE mutants with improved catalytic activity for (-)-cocaine.

Active site dynamics of acyl‐chymotrypsin

The motions of water molecules, the acyl moiety, the catalytic triad, and the oxyanion binding site of acyl-chymotrypsin were studied by means of a stochastic boundary molecular dynamics simulation. A water molecule that could provide the nucleophilic OH- for the deacylation stage of the catalysis was found to be trapped between the imidazole ring of His-57 and the carbonyl carbon of the acyl group. It makes a hydrogen bond with the N epsilon 2 of His-57 and is held in place through a network of hydrogen-bonded water molecules in the active site. The water molecule was found as close as 2.8 A to the carbonyl carbon. This appears to be due to the constraints imposed by nonbonded interaction in the active site. Configurations were found in which one hydrogen of the trapped water shared a bifurcated hydrogen bond with His-57-N epsilon 2 and Ser-195-O gamma, with the water oxygen very close to the carbonyl carbon. The existence of such a water molecule suggests that large movement of the His-57 imidazole ring between positions suitable for providing general-base catalyzed assistance and for providing general-acid catalyzed assistance may not be required during the reaction. The simulation indicates that the side chains of residues involved in catalysis (i.e., His-57, Ser-195, and Asp-102) are significantly less flexible than other side chains in the protein. The 40% reduction in rms fluctuations is consistent with a comparable reduction calculated from the temperature factors obtained in the X-ray crystallographic data of gamma-chymotrypsin. The greater rigidity of active site residues seems to result from interconnected hydrogen bonding networks among the residues and between the residues and the solvent water in the active site.

Comparison of Deacylation Rates of Chymotryptic Catalysis within an Enantiomeric Pair of p-Nitrophenyl Esters.

p-Nitrophenyl esters carrying a chiral acyl group were synthesized. These compounds were shown to meet the requirements of chymotrypsin for the specific binding and the acylation. Therefore, the behavior of p-nitrophenyl esters with chymotrypsin is considered to be virtually identical to that of p-amidinophenyl esters with trypsin which were proposed as inverse substrates for the enzyme. These esters, derived from each pair of enantiomers, have been successfully used for the analysis of enantiomeric preference of chymotrypsin at the deacylation stage. The chiral requirement of the enzyme active site for the catalytic efficiency was discussed.

Induced activation in the deacylation step of tryptic hydrolysis. An application of "inverse substrates" to mechanistic studies of the enzyme.

Trypsin [EC 3.4.21.4]-catalyzed hydrolysis of "inverse substrates" was investigated kinetically. "Inverse substrates" for trypsin are specific substrates in which the arrangement of the site-specific group is reversed compared to that of the normal substrate, e.g., a cationic center is included in the leaving group instead of being in the acyl moiety (Tanizawa, K., Kasaba, Y., & Kanaoka, Y. (1977) J. Am. Chem. Soc. 99, 4485-4488). Acyl enzyme intermediates formed specifically from these substrates are advantageous for the mechanistic analysis of trypsin action, since the cationic group liberated from the acyl moiety can no longer exhibit specific interaction with the enzyme binding site in the subsequent deacylation stage. Remarkable rate acceleration at the deacylation step was observed on adding amidinium or ammonium compounds. The effects of the size of the acyl moiety and the charged molecule on the acceleration were examined. Latent properties of 1-butylamine as an activator were found in the present study. Based on these observations, it is suggested that a cationic molecule which can be well accommodated together with the acyl group within the active center cleft causes rate enhancement, with associated conformational changes.

The effect of temperature on the individual stages of the hydrolysis of non-specific-p-nitrophenol esters by alpha-chymotrypsin.

Precise studies were performed on the effect of temperature on the rate and equilibrium parameters characterizing the individual stages of the alpha-chymotrypsin-catalysed hydrolysis of non-specific p-nitrophenol esters at pH 7.40 and 8.50. At both pH values the results indicate that a sharp kinetic anomaly is observed in Arrhenius plots of these parameters for the binding and acylation stages of the process, but not for the deacylation stage. Detailed comparison with other kinetic studies was made, and a comparison with thermal transitions observed in alpha-chymotrypsin by using physical techniques was attempted. A detailed discussion of possible causes of the anomalies is given.