A Comprehensive Understanding of Enzymatic Catalysis by Hydroxynitrile Lyases with S Stereoselectivity from the α/β-Hydrolase Superfamily: Revised Role of the Active-Site Lysine and Kinetic Behavior of Substrate Delivery and Sequential Product Release

Protective effect of Eucommia ulmoides Oliver male flowers on ethanol‐induced DNA damage in mouse cerebellum and cerebral cortex

Ethanol is a principal ingredient of alcoholic beverages with potential neurotoxicity and genotoxicity, and the ethanol‐associated oxidative DNA damage in the central nervous system is well documented. Natural product may offer new options to protect the brain against ethanol‐induced neurotoxicity. The male flower of Eucommia ulmoides (EUF) Oliver has been extensively utilized as the tea, the healthy hot drink on the market. In this study, 19 constituents in the effective fraction of EUF were identified by ultra‐performance liquid chromatography–tandem mass spectrometry (UPLC–MS/MS). In the single‐cell gel electrophoresis assay, EUF was observed to ameliorate DNA damage in mouse cerebellum and cerebral cortex caused by acute ethanol administration, which was further confirmed by the morphological observation. The protective effects of EUF were associated with increasing total superoxide dismutase (T‐SOD) and glutathione peroxidase (GSH‐PX) activities, and a decrease in nitric oxide (NO), malondialdehyde (MDA), 8‐hydroxy‐2′‐deoxyguanosine (8‐OHdG), and kelch‐like ECH‐associated protein‐1 (Keap1) levels. Molecular docking results demonstrated that compounds 4, 7, 9, and 16 from EUF have a strong affinity to the Keap1 Kelch domain to hinder the interaction of nuclear factor‐erythroid 2‐related factor 2 (Nrf2) with Keap1. These findings suggest that EUF is a potent inhibitor of ethanol‐induced brain injury possibly via the inhibition of oxidative stress.

Degradation of pesticides diazinon and diazoxon by phosphotriesterase: insight into divergent mechanisms from QM/MM and MD simulations

Enzymatic hydrolysis by phosphotriesterase (PTE) is one of the most effective ways of degrading organophosphorus pesticides, but the catalytic efficiency depends on the structural features of substrates. Here the enzymatic degradation of diazinon (DIN) and diazoxon (DON), characterized by PS and PO, respectively, have been investigated by QM/MM calculations and MM MD simulations. Our calculations demonstrate that the hydrolysis of DON (with PO) is inevitably initiated by the nucleophilic attack of the bridging-OH^- on the phosphorus center, while for DIN (with PS), we proposed a new degradation mechanism, initiated by the nucleophilic attack of the Zn_α-bound water molecule, for its low-energy pathway. For both DIN and DON, the hydrolytic reaction is predicted to be the rate-limiting step, with energy barriers of 18.5 and 17.7 kcal mol^-1, respectively. The transportation of substrates to the active site, the release of the leaving group and the degraded product are generally verified to be favorable by MD simulations via umbrella sampling, both thermodynamically and dynamically. The side-chain residues Phe132, Leu271 and Tyr309 play the gate-switching role to manipulate substrate delivery and product release. In comparison with the DON-enzyme system, the degraded product of DIN is more easily released from the active site. These new findings will contribute to the comprehensive understanding of the enzymatic degradation of toxic organophosphorus compounds by PTE.

QM/MM and MM MD simulations on decontamination of the V-type nerve agent VX by phosphotriesterase: Toward a comprehensive understanding of steroselectivity and activity

Due to deadly toxicity and high environmental stability of the nerve agent VX, an efficient decontamination approach is desperately needed in tackling its severe threat to human secu-rity. The enzymatic...

Computer-aided understanding and engineering of enzymatic selectivity

Enzymes offering chemo-, regio-, and stereoselectivity enable the asymmetric synthesis of high-value chiral molecules. Unfortunately, the drawback that naturally occurring enzymes are often inefficient or have undesired selectivity toward non-native substrates hinders the broadening of biocatalytic applications. To match the demands of specific selectivity in asymmetric synthesis, biochemists have implemented various computer-aided strategies in understanding and engineering enzymatic selectivity, diversifying the available repository of artificial enzymes. Here, given that the entire asymmetric catalytic cycle, involving precise interactions within the active pocket and substrate transport in the enzyme channel, could affect the enzymatic efficiency and selectivity, we presented a comprehensive overview of the computer-aided workflow for enzymatic selectivity. This review includes a mechanistic understanding of enzymatic selectivity based on quantum mechanical calculations, rational design of enzymatic selectivity guided by enzyme-substrate interactions, and enzymatic selectivity regulation via enzyme channel engineering. Finally, we discussed the computational paradigm for designing enzyme selectivity in silico to facilitate the advancement of asymmetric biosynthesis.

Combination of multiple computational methods revealing specific sub-sectional recognition and hydrogen-bond dependent transportation of CKII peptide fragment in O-GlcNAc transferase

•

Mechanism of CKII peptide recognition, transportation and binding in OGT is obtained.

•

Peptide delivery is strong exothermic, highly dependent on hydrogen bond network.

•

Typical ‘spread’ & ‘V’ conformation change noticed for peptide accompanies stable OGT.

•

Specific subsection of peptide has diverse performance in its recognition and delivery.

•

Multiple methods combination may be used in other bio-system with flexible substrate.

O-linked β-N-acetyl-D-glucosamine (O-GlcNAc) transferase (OGT) is an essential enzyme in many cellular physiological catalytic reactions that regulates protein O-GlcNAcylation. Aberrant O-GlcNAcylation is related to insulin resistance, diabetic complications, cancer and neurodegenerative diseases. Understanding the peptide delivery in OGT is significant in comprehending enzymatic catalytic process, target-protein recognition and pathogenic mechanism. Herein extensive molecular dynamics (MD) simulations combined with various techniques are utilized to study the recognizing and binding mechanism of peptide fragment extracted from casein kinase II by OGT from atomic level. The residues of His496, His558, Thr633, Lys634, and Pro897 are demonstrated to play a dominant role in the peptide stabilization via hydrogen bonds and σ-π interaction, whose van der Waals and non-polar solvent effects provide the main driving force. In addition, two channels are identified. The delivery mode, mechanism together with thermodynamic and dynamic characterizations for the most favorable channel are determined. The peptide is more inclined to be recognized by OGT through the cavity comprised of residues 799–812, 893–899, and 865–871, and Tyr13-terminal is prior recognized to Met26-terminal. The transportation process is accompanied with conformation changes between the “spread” and “V” shapes. The whole process is strong exothermic that is highly dependent on the variation of hydrogen bond interactions between peptide and OGT as well as the performance of different subsections of peptide. Besides that, multiple computational methods combinations may contribute meaningfully to calculation of similar bio-systems with long and flexible substrate.

Chemical compounds with a neuroprotective effect from the seeds of Celosia argentea L

Oxidative stress plays a central role in the common pathophysiology of neurodegenerative diseases such as Alzheimer's disease, amyotrophic lateral sclerosis, and Parkinson's disease. Antioxidant therapy has been suggested for the prevention and treatment of neurodegenerative diseases. Compounds derived from natural sources may offer the potential for new treatment options. Semen Celosiae is a traditional Chinese edible herbal medicine with a long history in China and exhibits wide-reaching biological activities such as hepatoprotective, anti-tumor, anti-diarrheal, anti-diabetic, anti-oxidant, etc. In this study, nine saponins and two phenylacetonitrile glycosides were isolated from Semen Celosiae and their structures were identified using ESI-MS and NMR techniques. Among them, compounds 1 and 2 have not been previously reported. The total concentrations of the five triterpenoid saponins and the two phenylacetonitrile glycosides were 3.348 mg g^-1 and 0.187 mg g^-1, respectively, suggesting that Semen Celosiae is a novel viable source of the two kinds of compounds. These compounds were observed to significantly attenuate t-BHP-induced neuronal damage by effectively enhancing cell viability and decreasing reactive oxygen species generation and cell apoptosis rate in NSC-34 cells. Furthermore, compounds 1 and 7 reduced the ratios of cleaved caspase-3: caspase-3 and cleaved caspase-7: caspase-7 and the level of cytochrome C, while they increased the levels of SOD1 and Beclin 1. These findings suggest that compounds 1-11 are potent inhibitors of neuron injury elicited by t-BHP, possibly via inhibition of oxidative stress and apoptosis, and activation of autophagy; therefore they may be valuable leads for future therapeutic development.

Molecular Dynamics Revealing a Detour-Forward Release Mechanism of Tacrine: Implication for the Specific Binding Characteristics in Butyrylcholinesterase

Butyrylcholinesterase (BChE) is a non-specific enzyme with clinical pharmacological and toxicological significance, which was a renewed interest as therapeutic target in Alzheimer's disease (AD) nowadays. Here, all-atom molecular dynamics simulations of butyrylcholinesterase with tacrine complex were designed to characterize inhibitor binding modes, strengths, and the hydrogen-bond dependent non-covalent release mechanism. Four possible release channels were identified, and the most favorable channel was determined by random acceleration molecular dynamics molecular dynamics (RAMD MD) simulations. The thermodynamic and dynamic properties as well as the corresponding Detour-forward delivery mechanism were determined according to the classical molecular dynamics (MD) simulations accompanied with umbrella sampling. The free energy barrier of the tacrine release process for the most beneficial pathway is about 10.95 kcal/mol, which is related to the non-covalent interactions from the surrounding residues, revealing the specific binding characteristics in the active site. The residues including Asp70, Ser79, Trp82, Gly116, Thr120, Tyr332, and His438 were identified to play major roles in the stabilization of tacrine in the pocket of BChE, where hydrogen bonding and π-π interactions are significant factors. Tyr332 and Asp70, which act as gate keepers, play crucial roles in the substrate delivery. The present results provide a basic understanding for the ligand transport mechanism depending on the BChE enzymatic environment, which is useful for the design of BChE inhibitors in the future.

Molecular dynamics simulations of an engineered T4 lysozyme exclude helix to sheet transition, and provide insights into long distance, intra-protein switchable motion

An engineered variant of T4 lysozyme serves as a model for studying induced remote conformational changes in a full protein context. The design involves a duplicated surface helix, flanked by two loops, that switches between two different conformations spanning about 20 Å. Molecular dynamics simulations of the engineered protein, up to 1 μs, rule out α-helix to β-sheet transitions within the duplicated helix as suggested by others. These simulations highlight how the use of different force fields can lead to radical differences in the structure of the protein. In addition, Markov state modeling and transition path theory were employed to map a 6.6 μs simulation for possible early intermediate states and to provide insights into the onset of the switching motion. The putative intermediates involve the folding of one helical turn in the C-terminal loop through energy driven, sequential rearrangement of nearby salt bridges around the key residue Arg63. These results provide a first step towards understanding the energetics and dynamics of a rather complicated intra-protein motion.

Linking inhibitor motions to proteolytic stability of sunflower trypsin inhibitor-1

The remarkable capability of an enzyme isn't only determined by its active site but also controlled by the environment. To unravel the environment role in catalysis, the dynamic motions as well as the static mechanism need to be studied. In this work, QM/MM MD simulations were employed to study the proteolysis process of SFTI-1 and BiKF, which revealed that a combination of static non-bonded interactions and dynamic motions along the reaction coordinate can account for the different hydrolysis rates between them. A comparison among SFTI-1 and three analogs with similar non-bonded interactions further revealed a positive correlation between the mobility of inhibitors and the hydrolysis rates. Apart from the cyclic backbone and disulfide bond, intramolecular hydrogen bonds also increase the rigidity of the backbone of inhibitors, and therefore hinder inhibitor motions to resist proteolysis. These new detailed mechanistic insights suggest the need to consider inhibitor motions in the rational design of peptide inhibitors.

Besides the non-bonded interactions, inhibitor motions especially rotation of the scissile bond also influence proteolytic stability.

Insight into the delivery channel and selectivity of multiple binding sites in bovine serum albumin towards naphthalimide-polyamine derivatives

Naphthalimide derivatives are types of small-molecule anticancer drug candidates; however, their negative factors and potential side effects make their application limited. The pharmacophores select a direct access into the tumor cells as the first choice; this can reduce the side effect of the anti-cancer drugs on the normal cells. Herein, the delivery and binding of the naphthalimide-polyamine complex assisted by the bovine serum albumin (BSA) protein have been studied by combining several molecular dynamic simulations. The plausible transportation channels and the most favorable pathways for the delivery of the naphthalimide-polyamine complex to two drug sites (DSI and DSII), their thermodynamic and dynamic properties and the mechanisms have been discussed in detail. The residues His287 and Phe394 acted as guards in the DSI and DSII, respectively, which played a gating-switch role by flipping the ring from open to close during the compound delivery. The binding mode, binding energy and substituent effects have been also identified. The two drug sites have different preferences towards the compound with the electron-withdrawing and electron-donating substituents, and their strong interactions are more sensitive to the number of the substituent groups. The naphthalimide-polyamine complexes are more likely to choose DSI, both thermodynamically and dynamically, as compared to DSII. This selective specificity of these two drug sites manipulated by the electron-withdrawing and electron-donating substituents is quite promising for the design of new naphthalimide drugs.

Molecular insight into chymotrypsin inhibitor 2 resisting proteolytic degradation

Chymotrypsin inhibitor 2 (CI2) is a special serine protease inhibitor which can resist hydrolysis for several days with a rapid equilibrium between the Michaelis complex and acyl-enzyme intermediate. The energies and conformational changes for subtilisin-catalyzed proteolysis of CI2 were examined in this paper for the first time by employing pseudo bond ab initio QM/MM MD simulations. In the acylation reaction, a low-barrier hydrogen bond between His64 and Asp32 in the transition state together with the lack of covalent backbone constraints makes the peptide bonds of CI2 break more easily than in other serine protease inhibitors. After acyl-enzyme formation, molecular dynamics simulations showed that the access of hydrolytic water to the active site requires partial dissociation of the leaving group. However, retention of the leaving group mainly by the intra- and inter-molecular H-bonding networks hinders the access of water and retards the deacylation reaction. Instead of the dissociation constant of inhibitors, we suggest employing the free energy at the acyl-enzyme state to predict the relative hydrolysis rates of CI2 mutants, which are testified by the experimental relative hydrolysis rates.

Uridine diphosphate release mechanism in O-N-acetylglucosamine (O-GlcNAc) transferase catalysis

O-linked N-acetylglucosamine transferase (OGT) is an essential enzyme that catalyzes the covalent bonding of N-acetylglucosamine (GlcNAc) to the hydroxyl group of a serine or threonine in the target protein. It plays an important role in many important cellular physiological catalytic reactions. Here, we determine the binding mode and the binding free energy of the OGT product (uridine diphosphate, UDP) as well as the hydrogen-bond-dependent release mechanism using extensive molecular dynamic simulations. The Lys634, Asn838, Gln839, Lys842, His901, and Asp925 residues were identified to play a major role in the UDP stabilization in the active site of OGT, where hydrogen bonding and π-π interactions mainly occur. The calculations on the mutant forms support our results. Sixteen possible release channels were identified while the two most favorable channels were determined using random acceleration molecular dynamics (RAMD) simulations combined with the constant velocity pulling (PCV) method. The thermodynamic and dynamic properties as along with the corresponding mechanism were determined and discussed according to the umbrella sampling technique. For the most optimal channel, the main free energy barrier is 13 kcal/mol, which probably originates from the hydrogen bonds between UDP and the Ala896 and Asp925 residues. Moreover, the unstable hydrogen bonds and the rollback of the ligand likely cause the other two small obstacles. This work clarifies the ligand transport mechanism in the OGT enzymatic process and is a great resource for designing inhibitors based on UDP or UDP-GlcNAc.

Catalytic promiscuity of the non-native FPP substrate in the TEAS enzyme: non-negligible flexibility of the carbocation intermediate

By QM(DFT)/MM MD simulations, it has been revealed that the non-native substrate catalytic promiscuity of TEAS (one of the sesquiterpene cyclases) is mostly attributable to its notable conformational flexibility of the branching intermediate bisabolyl cation.

Product release mechanism and the complete enzyme catalysis cycle in yeast cytosine deaminase (yCD): A computational study

Yeast cytosine deaminase (yCD) is critical in gene-directed enzyme prodrug therapy as it catalyzes the hydrolytic deamination of cytosine. The product (uracil) release process is considered as rate-limiting in the whole enzymatic catalysis and includes the cleavage of the uracil-metal bond and the delivery of free uracil out of the reactive site. Herein extensive combined random acceleration molecular dynamics (RAMD) and molecular dynamics (MD) simulations coupled with the umbrella sampling technique have been performed to study the product transport mechanism. Five channels have been identified, and the thermodynamic and dynamic characterizations for the two most favorable channels have been determined and analyzed. The free energy barrier for the most beneficial pathway is about 13kcal/mol and mainly results from the cleavage of hydrogen bonds between the ligand uracil and surrounding residues Asn51, Glu64, and Asp155. The conjugated rings of Phe114 and Trp152 play gating and guiding roles in the product delivery via π⋯π van der Waals interactions with the product. Finally, the full cycle of the enzymatic catalysis has been determined, making the whole process computationally more precise.