Functional and structural changes due to a serine to alanine mutation in the active-site flap of enolase

Multifunctional roles of γ-enolase in the central nervous system: more than a neuronal marker

Enolase, a multifunctional protein with diverse isoforms, has generally been recognized for its primary roles in glycolysis and gluconeogenesis. The shift in isoform expression from α-enolase to neuron-specific γ-enolase extends beyond its enzymatic role. Enolase is essential for neuronal survival, differentiation, and the maturation of neurons and glial cells in the central nervous system. Neuron-specific γ-enolase is a critical biomarker for neurodegenerative pathologies and neurological conditions, not only indicating disease but also participating in nerve cell formation and neuroprotection and exhibiting neurotrophic-like properties. These properties are precisely regulated by cysteine peptidase cathepsin X and scaffold protein γ1-syntrophin. Our findings suggest that γ-enolase, specifically its C-terminal part, may offer neuroprotective benefits against neurotoxicity seen in Alzheimer's and Parkinson's disease. Furthermore, although the therapeutic potential of γ-enolase seems promising, the effectiveness of enolase inhibitors is under debate. This paper reviews the research on the roles of γ-enolase in the central nervous system, especially in pathophysiological events and the regulation of neurodegenerative diseases.

Dopamine modification of glycolytic enzymes impairs glycolysis: possible implications for Parkinson's disease

BACKGROUND

Parkinson's disease (PD), a chronic and severe neurodegenerative disease, is pathologically characterized by the selective loss of nigrostriatal dopaminergic neurons. Dopamine (DA), the neurotransmitter produced by dopaminergic neurons, and its metabolites can covalently modify proteins, and dysregulation of this process has been implicated in neuronal loss in PD. However, much remains unknown about the protein targets.

METHODS

In the present work, we designed and synthesized a dopamine probe (DA-P) to screen and identify the potential protein targets of DA using activity-based protein profiling (ABPP) technology in combination with liquid chromatography-tandem mass spectrometry (LC-MS/MS). In situ pull-down assays, cellular thermal shift assays (CETSAs) and immunofluorescence were performed to confirm the DA modifications on these hits. To investigate the effects of DA modifications, we measured the enzymatic activities of these target proteins, evaluated glycolytic stress and mitochondrial respiration by Seahorse tests, and systematically analyzed the changes in metabolites with unbiased LC-MS/MS-based non-targeted metabolomics profiling.

RESULTS

We successfully identified three glycolytic proteins, aldolase A, α-enolase and pyruvate kinase M2 (PKM2), as the binding partners of DA. DA bound to Glu166 of α-enolase, Cys49 and Cys424 of PKM2, and Lys230 of aldolase A, inhibiting the enzymatic activities of α-enolase and PKM2 and thereby impairing ATP synthesis, resulting in mitochondrial dysfunction.

CONCLUSIONS

Recent research has revealed that enhancing glycolysis can offer protection against PD. The present study identified that the glycolytic pathway is vulnerable to disruption by DA, suggesting a promising avenue for potential therapeutic interventions. Safeguarding glycolysis against DA-related disruption could be a potential therapeutic intervention for PD.

Serine 106 preserves the tertiary structure, function, and stability of a cyclophilin from Staphylococcus aureus

SaCyp, a staphylococcal cyclophilin involved in both protein folding and pathogenesis, has a Ser residue at position 106 and a Trp residue at position 136. While Ser 106 of SaCyp aligned with a cyclosporin A (CsA) binding Ala residue, its Trp 136 aligned with a Trp or a Phe residue of most other cyclophilins. To demonstrate the exact roles of Ser 106 and Trp 136 in SaCyp, we have elaborately studied rCyp[S106A] and rCyp[W136A], two-point mutants of a recombinant SaCyp (rCyp) harboring an Ala substitution at positions 106 and 136, respectively. Of the mutants, rCyp[W136A] showed the rCyp-like CsA binding affinity and peptidyl-prolyl cis-trans isomerase (PPIase) activity. Conversely, the PPIase activity, CsA binding affinity, stability, tertiary structure, surface hydrophobicity, and Trp accessibility of rCyp[S106A] notably differed from those of rCyp. The computational experiments also reveal that the structure, dimension, and fluctuation of SaCyp are not identical to those of SaCyp[S106A]. Furthermore, Ser at position 106 of SaCyp, compared to Ala at the same position, formed a higher number of non-covalent bonds with CsA. Collectively, Ser 106 is an indispensable residue for SaCyp that keeps its tertiary structure, function, and stability intact.Communicated by Ramaswamy H. Sarma.

Rational Metabolic Engineering Combined with Biosensor-Mediated Adaptive Laboratory Evolution for l-Cysteine Overproduction from Glycerol in Escherichia coli

l-Cysteine is an important sulfur-containing amino acid with numerous applications in the pharmaceutical and cosmetic industries. The microbial production of l-cysteine has received substantial attention, and the supply of the precursor l-serine is important in l-cysteine biosynthesis. In this study, to achieve l-cysteine overproduction, we first increased l-serine production by deleting genes involved in the pathway of l-serine degradation to glycine (serine hydroxymethyl transferase, SHMT, encoded by glyA genes) in strain 4W (with l-serine titer of 1.1 g/L), thus resulting in strain 4WG with l-serine titer of 2.01 g/L. Second, the serine-biosensor based on the transcriptional regulator NCgl0581 of C. glutamicum was constructed in E. coli, and the validity and sensitivity of the biosensor were demonstrated in E. coli. Then 4WG was further evolved through adaptive laboratory evolution (ALE) combined with serine-biosensor, thus yielding the strain 4WGX with 4.13 g/L l-serine production. Moreover, the whole genome of the evolved strain 4WGX was sequenced, and ten non-synonymous mutations were found in the genome of strain 4WGX compared with strain 4W. Finally, 4WGX was used as the starting strain, and deletion of the l-cysteine desulfhydrases (encoded by tnaA), overexpression of serine acetyltransferase (encoded by cysE) and the key enzyme of transport pathway (encoded by ydeD) were performed in strain 4WGX. The recombinant strain 4WGX-∆tnaA-cysE-ydeD can produce 313.4 mg/L of l-cysteine using glycerol as the carbon source. This work provides an efficient method for the biosynthesis of value-added commodity products associated with glycerol conversion.

Substrate-to-Product Conversion Facilitates Active Site Loop Opening in Yeast Enolase: A Molecular Dynamics Study

Yeast enolase serves as a prototype for metalloenzymes with labile, catalytic active site metal ions and is important for glycolysis and fermentation processes. Herein, microsecond molecular dynamics simulations of the protein-substrate and protein-product complexes are conducted to elucidate the mechanism of the opening of catalytically important active site loops. These simulations indicate that conversion of substrate to product is accompanied by diminished metal coordination and hydrogen-bonding interactions, as well as enhanced loop flexibility. Moreover, free energy simulations show that the loop opening is endergonic when substrate is bound but exergonic when product is bound. Thus, the conversion to product weakens the association of the loop with the ligand and binding site, thereby facilitating the loop opening after catalysis and enabling product release. These insights about active site loop motions in enzyme catalysis may be useful in guiding enzyme design efforts.

Genetic variants in cardiac calcification in Northern Sweden

Extensive coronary calcification without significant stenosis, described as calcific coronary artery disease (CCAD) may cause abnormal myocardial perfusion and hence generalized ischemia. There is a discrepancy in the expression pattern of CCAD compared to the well-known atherosclerotic disease which raises questions about the exact pathophysiology of coronary calcification and whether there is a genetic etiology for it.In this pilot study we studied 3 candidate genes, ectonucleotide pyrophosphatase/phosphodiesterase (ENPP1), ATP Binding Cassette Subfamily C Member 6 (ABCC6), and 5'-Nucleotidase Ecto (NT5E) involved in pyrophosphate (PPi) and inorganic phosphate (Pi) metabolism, which may predispose to coronary arterial or valvular calcification. We studied 70 patients with calcific cardiac disease; 65 with CCAD (age 43-83 years) and 5 with calcific aortic valve disease (CAVD) (age 76-82 years).Five DNA variants potentially affecting protein function were found in 6 patients. One variant is a known disease-causing mutation in the ABCC6 gene. Our findings support that disturbances in the PPi and Pi metabolism might influence the development of CCAD and CAVD. However, segregation in the families must first be performed to ascertain any damaging effect of these variants we have found.We report 4 new genetic variants potentially related to coronary calcification, through the disturbed Pi and PPi metabolism. The search for direct causative genetic variants in coronary artery and aortic valve calcification must be broadened with other genes particularly those involved with Pi and PPi metabolism.

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.

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.

Why does the G117H mutation considerably improve the activity of human butyrylcholinesterase against sarin? Insights from quantum mechanical/molecular mechanical free energy calculations

Human butyrylcholinesterase (BChE) is recognized as the most promising bioscavenger for organophosphorus (OP) warfare nerve agents. The G117H mutant of human BChE has been identified as a potential catalytic bioscavenger with a remarkably improved activity against OP nerve agents such as sarin, but it still does not satisfy the clinical use. For further design of the higher-activity mutants against OP nerve agents, it is essential to understand how the G117H mutation improves the activity. The reaction mechanisms and the free energy profiles for spontaneous reactivation of wild-type BChE and its G117H mutant phosphorylated by sarin have been explored, in this study, by performing first-principles quantum mechanical/molecular mechanical free energy calculations, and the remarkable role of the G117H mutation on the activity has been elucidated. For both the wild-type and G117H mutant enzymes, H438 acts as a general base to initiate the spontaneous reactivation that consists of two reaction steps: the nucleophilic attack at the phosphorus by a water molecule and decomposition of the pentacoordinated phosphorus intermediate. The calculated overall free energy barriers, i.e., 30.2 and 23.9 kcal/mol for the wild type and G117H mutant, respectively, are in good agreement with available experimental kinetic data. On the basis of the calculated results, the mutated residue (H117 in the G117H mutant) cannot initiate the spontaneous reactivation as a general base. Instead, it skews the oxyanion hole and makes the phosphorus more open to the nucleophilic water molecule, resulting in a remarkable change in the rate-determining step and significantly improved catalytic activity of human BChE.

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).

Fundamental reaction pathway and free energy profile for hydrolysis of intracellular second messenger adenosine 3',5'-cyclic monophosphate (cAMP) catalyzed by phosphodiesterase-4

As important drug targets for a variety of human diseases, cyclic nucleotide phosphodiesterases (PDEs) are a superfamily of enzymes sharing a similar catalytic site. We have performed pseudobond first-principles quantum mechanical/molecular mechanical-free energy perturbation (QM/MM-FE) and QM/MM-Poisson-Boltzmann surface area (PBSA) calculations to uncover the detailed reaction mechanism for PDE4-catalyzed hydrolysis of adenosine 3',5'-cyclic monophosphate (cAMP). This is the first report on QM/MM reaction-coordinate calculations including the protein environment of any PDE-catalyzed reaction system, demonstrating a unique catalytic reaction mechanism. The QM/MM-FE and QM/MM-PBSA calculations revealed that the PDE4-catalyzed hydrolysis of cAMP consists of two reaction stages: cAMP hydrolysis (stage 1) and bridging hydroxide ion regeneration (stage 2). The stage 1 includes the binding of cAMP in the active site, nucleophilic attack of the bridging hydroxide ion on the phosphorus atom of cAMP, cleavage of O3'-P phosphoesteric bond of cAMP, protonation of the departing O3' atom, and dissociation of hydrolysis product (AMP). The stage 2 includes the binding of solvent water molecules with the metal ions in the active site and regeneration of the bridging hydroxide ion. The dissociation of the hydrolysis product is found to be rate-determining for the enzymatic reaction process. The calculated activation Gibbs free energy of ≥16.0 and reaction free energy of -11.1 kcal/mol are in good agreement with the experimentally derived activation free energy of 16.6 kcal/mol and reaction free energy of -11.5 kcal/mol, suggesting that the catalytic mechanism obtained from this study is reliable and provides a solid base for future rational drug design.

Structure analysis of Entamoeba histolytica enolase

Entamoeba histolytica enolase (EhENO) reversibly interconverts 2-phosphoglyceric acid (2-PGA) and phosphoenolpyruvic acid (PEP). The crystal structure of the homodimeric EhENO is presented at a resolution of 1.9 Å. In the crystal structure EhENO presents as an asymmetric dimer with one active site in the open conformation and the other active site in the closed conformation. Interestingly, both active sites contain a copurified 2-PGA molecule. While the 2-PGA molecule in the closed active site closely resembles the conformation known from other enolase-2-PGA complexes, the conformation in the open active site is different. Here, 2-PGA is shifted approximately 1.6 Å away from metal ion I, most likely representing a precatalytic situation.

Catalytic mechanism of cytochrome P450 for 5'-hydroxylation of nicotine: fundamental reaction pathways and stereoselectivity

A series of computational methods were used to study how cytochrome P450 2A6 (CYP2A6) interacts with (S)-(-)-nicotine, demonstrating that the dominant molecular species of (S)-(-)-nicotine in CYP2A6 active site exists in the free base state (with two conformations, SR(t) and SR(c)), despite the fact that the protonated state is dominant for the free ligand in solution. The computational results reveal that the dominant pathway of nicotine metabolism in CYP2A6 is through nicotine free base oxidation. Further, first-principles quantum mechanical/molecular mechanical free energy (QM/MM-FE) calculations were carried out to uncover the detailed reaction pathways for the CYP2A6-catalyzed nicotine 5'-hydroxylation reaction. In the determined CYP2A6-(S)-(-)-nicotine binding structures, the oxygen of Compound I (Cpd I) can abstract a hydrogen from either the trans-5'- or the cis-5'-position of (S)-(-)-nicotine. CYP2A6-catalyzed (S)-(-)-nicotine 5'-hydroxylation consists of two reaction steps, that is, the hydrogen transfer from the 5'-position of (S)-(-)-nicotine to the oxygen of Cpd I (the H-transfer step), followed by the recombination of the (S)-(-)-nicotine moiety with the iron-bound hydroxyl group to generate the 5'-hydroxynicotine product (the O-rebound step). The H-transfer step is rate-determining. The 5'-hydroxylation proceeds mainly with the stereoselective loss of the trans-5'-hydrogen, that is, the 5'-hydrogen trans to the pyridine ring. The calculated overall stereoselectivity of ∼97% favoring the trans-5'-hydroxylation is close to the observed stereoselectivity of 89-94%. This is the first time it has been demonstrated that a CYP substrate exists dominantly in one protonation state (cationic species) in solution, but uses its less-favorable protonation state (neutral free base) to perform the enzymatic reaction.

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.

Engineering the enolase magnesium II binding site: implications for its evolution

The glycolytic enzyme enolase catalyzes the reversible elimination of water from 2-phosphoglycerate (2-PGA) to form phosphoenolpyruvate (PEP). Two magnesium ions in the active site are thought to facilitate the reaction by activation of the C2 proton of 2-PGA and charge stabilization of the intermediate. The initial abstraction of a proton from a carboxylic acid is common to all members of the enolase superfamily, yet in all other known members of this superfamily, only one magnesium ion (MgI) per active site is sufficient to promote catalysis. We wanted to further investigate the importance of the second magnesium ion (MgII) for the catalytic mechanism of yeast enolase 1. Toward this end, we removed all MgII coordinating residues and replaced substrate-MgII interactions by introducing positively charged side chains. High-resolution crystal structures and activity assays show that the introduced positively charged side chains effectively prohibit MgII binding but fail to promote catalysis. We conclude that enolase is inactive without MgII, yet control mutants without additional positively charged side chains retain basal enolase activity through binding of magnesium to 2-PGA in an open active site without the help of MgII coordinating residues. Thus, we believe that ancestral enolase activity might have evolved in a member of the enolase superfamily that provides only the necessary catalytic residues and the binding site for MgI. Additionally, precatalytic binding of 2-PGA to the apo state of enolase was observed.

Reaction pathway and free-energy barrier for reactivation of dimethylphosphoryl-inhibited human acetylcholinesterase

The dephosphorylation/reactivation mechanism and the corresponding free-energy profile of the dimethylphosphoryl-inhibited conjugate of human acetylcholinesterase (AChE) has been studied by performing first-principles quantum mechanical/molecular mechanical free-energy (QM/MM-FE) calculations. On the basis of the QM/MM-FE results, for the favorable reaction pathway, the entire dephosphorylation/reactivation process consists of three reaction steps, including the nucleophilic water attack on the P atom, the spatial reorganization of the dimethylphosphoryl group, and the dissociation between the dimethylphosphoryl group and Ser203 of AChE. The overall free-energy barrier for the entire dephosphorylation/reactivation reaction is found to be the free-energy change from the initial reactant to the transition state associated with the spatial reorganization step, and the calculated overall free-energy barrier (20.1 to 23.5 kcal/mol) is reasonably close to the experimentally derived activation free energy of 22.3 kcal/mol. In addition, key amino acid residues and their specific roles in the reaction process have been identified.

Fundamental reaction mechanism and free energy profile for (-)-cocaine hydrolysis catalyzed by cocaine esterase

The fundamental reaction mechanism of cocaine esterase (CocE)-catalyzed hydrolysis of (-)-cocaine and the corresponding free energy profile have been studied by performing pseudobond first-principles quantum mechanical/molecular mechanical free energy (QM/MM-FE) calculations. On the basis of the QM/MM-FE results, the entire hydrolysis reaction consists of four reaction steps, including the nucleophilic attack on the carbonyl carbon of (-)-cocaine benzoyl ester by the hydroxyl group of Ser117, dissociation of (-)-cocaine benzoyl ester, nucleophilic attack on the carbonyl carbon of (-)-cocaine benzoyl ester by water, and finally dissociation between the (-)-cocaine benzoyl group and Ser117 of CocE. The third reaction step involving the nucleophilic attack of a water molecule was found to be rate-determining, which is remarkably different from (-)-cocaine hydrolysis catalyzed by wild-type butyrylcholinesterase (BChE; where the formation of the prereactive BChE-(-)-cocaine complex is rate-determining) or its mutants containing Tyr332Gly or Tyr332Ala mutation (where the first chemical reaction step is rate-determining). Besides, the role of Asp259 in the catalytic triad of CocE does not follow the general concept of the "charge-relay system" for all serine esterases. The free energy barrier calculated for the rate-determining step of CocE-catalyzed hydrolysis of (-)-cocaine is 17.9 kcal/mol, which is in good agreement with the experimentally derived activation free energy of 16.2 kcal/mol. In the present study, where many sodium ions are present, the effects of counterions are found to be significant in determining the free energy barrier. The finding of the significant effects of counterions on the free energy barrier may also be valuable in guiding future mechanistic studies on other charged enzymes.

Highly dissociative and concerted mechanism for the nicotinamide cleavage reaction in Sir2Tm enzyme suggested by ab initio QM/MM molecular dynamics simulations

Sir2 enzymes catalyze the NAD+-dependent protein deacetylation and play critical roles in epigenetics, cell death, and lifespan regulation. In spite of a current flurry of experimental studies, the catalytic mechanism for this unique and important class of enzymes remains elusive. Employing on-the-fly Born-Oppenheimer molecular dynamics simulations with the B3LYP/6-31G(d) QM/MM potential and the umbrella sampling method, we have characterized the initial step of the Sir2Tm-catalyzed reaction, which is also the most controversial portion of its mechanism. Our results indicate that the nicotinamide cleavage reaction employs a highly dissociative and concerted displacement mechanism: the cleavage of the glycosidic bond is facilitated by the nucleophilic participation of the acetyl-lysine, and the dissociative transition state has a significant oxocarbenium ion character. During this step of the reaction, the Sir2Tm enzyme strongly stabilizes the covalent O-alkylamidate intermediate whereas its effect on the transition state is quite minimal. In addition, functional roles of key residues and motifs have been elucidated. This work further demonstrates the feasibility and applicability of the state-of-the-art ab initio QM/MM molecular dynamics approach in simulating enzyme reactions.

Ab initio quantum mechanical/molecular mechanical molecular dynamics simulation of enzyme catalysis: the case of histone lysine methyltransferase SET7/9

To elucidate enzyme catalysis through computer simulation, a prerequisite is to reliably compute free energy barriers for both enzyme and solution reactions. By employing on-the-fly Born-Oppenheimer molecular dynamics simulations with the ab initio quantum mechanical/molecular mechanical approach and the umbrella sampling method, we have determined free energy profiles for the methyl-transfer reaction catalyzed by the histone lysine methyltransferase SET7/9 and its corresponding uncatalyzed reaction in aqueous solution, respectively. Our calculated activation free energy barrier for the enzyme catalyzed reaction is 22.5 kcal/mol, which agrees very well with the experimental value of 20.9 kcal/mol. The difference in potential of mean force between a corresponding prereaction state and the transition state for the solution reaction is computed to be 30.9 kcal/mol. Thus, our simulations indicate that the enzyme SET7/9 plays an essential catalytic role in significantly lowering the barrier for the methyl-transfer reaction step. For the reaction in solution, it is found that the hydrogen bond network near the reaction center undergoes a significant change, and there is a strong shift in electrostatic field from the prereaction state to the transition state, whereas for the enzyme reaction, such an effect is much smaller and the enzyme SET7/9 is found to provide a preorganized electrostatic environment to facilitate the methyl-transfer reaction. Meanwhile, we find that the transition state in the enzyme reaction is a little more dissociative than that in solution.

Catalytic mechanism and product specificity of the histone lysine methyltransferase SET7/9: an ab initio QM/MM-FE study with multiple initial structures

Histone lysine methylation is emerging as an important mechanism to regulate chromatin structure and gene activity. To provide theoretical understanding of its reaction mechanism and product specificity, ab initio quantum mechanical/molecular mechanical free energy (QM/MM-FE) calculations and molecular dynamics simulations have been carried out to investigate the histone lysine methyltransferase SET7/9. It is found that the methyl-transfer reaction catalyzed by SET7/9 is a typical in-line S(N)2 nucleophilic substitution reaction with a transition state of 70% dissociative character. The calculated average free energy barrier at the MP2(6-31+G) QM/MM level is 20.4 +/- 1.1 kcal/mol, consistent with the activation barrier of 20.9 kcal/mol estimated from the experimental reaction rate. The barrier fluctuation has a strong correlation with the nucleophilic attack distance and angle in the reactant complex. The calculation results show that the product specificity of SET7/9 as a monomethyltransferase is achieved by disrupting the formation of near-attack conformations for the dimethylation reaction.

Structure and catalytic properties of an engineered heterodimer of enolase composed of one active and one inactive subunit

Enolase is a dimeric enzyme that catalyzes the interconversion of 2-phospho-D-glycerate and phosphoenolpyruvate. This reversible dehydration is effected by general acid-base catalysis that involves, principally, Lys345 and Glu211 (numbering system of enolase 1 from yeast). The crystal structure of the inactive E211Q enolase shows that the protein is properly folded. However, K345 variants have, thus far, failed to crystallize. This problem was solved by crystallization of an engineered heterodimer of enolase. The heterodimer was composed of an inactive subunit that has a K345A mutation and an active subunit that has N80D and N126D surface mutations to facilitate ion-exchange chromatographic separation of the three dimeric species. The structure of this heterodimeric variant, in complex with substrate/product, was obtained at 1.85 A resolution. The structure was compared to a new structure of wild-type enolase obtained from crystals belonging to the same space group. Asymmetric dimers having one subunit exhibiting two of the three active site loops in an open conformation and the other in a conformation having all three loops closed appear in both structures. The K345A subunit of the heterodimer is in the loop-closed conformation; its Calpha carbon atoms closely match those of the corresponding subunit of wild-type enolase (root-mean-squared deviation of 0.23 A). The kcat and kcat/Km values of the heterodimer are approximately half those of the N80D/N126D homodimer, which suggests that the subunits in solution are kinetically independent. A comparison of enolase structures obtained from crystals belonging to different space groups suggests that asymmetric dimers can be a consequence of the asymmetric positioning of the subunits within the crystal lattice.

Changing the metal ion selectivity of rabbit muscle enolase by mutagenesis: effects of the G37A and G41A mutations

During the reaction catalyzed by enolase, a mobile loop, residues 36-45, closes over the active site. In order to probe the role of this loop movement in catalysis, the glycines at positions 37 and 41 of rabbit muscle enolase (beta beta) have been mutated to alanines. The mutant forms-G37A and G41A-of enolase are both active, but have altered selectivity for divalent cations. G37A, when assayed with Mg(2+), has 12% the activity of the wild type. However, it is twice as active as wild type when assayed with Mn(2+), Zn(2+), or Co(2+). G41A has 4% the activity of the wild type with Mg(2+), is more active than wild type with Co(2+), and slightly less active than wild type with Mn(2+) and Zn(2+). The kinetic isotope effect for both mutants is greater than that of the wild type with all 4 divalent cations. These results indicate that the flexibility of this loop has subtle effects on catalytic activity.

How does the cAMP-dependent protein kinase catalyze the phosphorylation reaction: an ab initio QM/MM study

We have carried out density functional theory QM/MM calculations on the catalytic subunit of cAMP-dependent protein kinase (PKA). The QM/MM calculations indicate that the phosphorylation reaction catalyzed by PKA is mainly dissociative, and Asp166 serves as the catalytic base to accept the proton delivered by the substrate peptide. Among the key interactions in the active site, the Mg(2+) ions, glycine rich loop, and Lys72 are found to stabilize the transition state through electrostatic interactions. On the other hand, Lys168, Asn171, Asp184, and the conserved waters bound to Mg(2+) ions do not directly contribute to lower the energy barrier of the phosphorylation reaction, and possible roles for these residues are proposed. The QM/MM calculations with different QM/MM partition schemes or different initial structures yield consistent results. In addition, we have carried out 12 ns molecular dynamics simulations on both wild type and K168A mutated PKA, respectively, to demonstrate that the catalytic role of Lys168 is to keep ATP and substrate peptide in the near-attack reactive conformation.

Enzymatic function of loop movement in enolase: preparation and some properties of H159N, H159A, H159F, and N207A enolases

The hypothesis that His159 in yeast enolase moves on a polypeptide loop to protonate the phosphoryl of 2-phosphoglycerate to initiate its conversion to phosphoenolpyruvate was tested by preparing H159N, H159A, and H159F enolases. These have 0.07%-0.25% of the native activity under standard assay conditions and the pH dependence of maximum velocities of H159A and H159N mutants is markedly altered. Activation by Mg2+ is biphasic, with the smaller Mg2+ activation constant closer to that of the "catalytic" Mg2+ binding site of native enolase and the larger in the mM range in which native enolase is inhibited. A third Mg2+ may bind to the phosphoryl, functionally replacing proton donation by His159. N207A enolase lacks an intersubunit interaction that stabilizes the closed loop(s) conformation when 2-phosphoglycerate binds. It has 21% of the native activity, also exhibits biphasic Mg2+ activation, and its reaction with the aldehyde analogue of the substrate is more strongly inhibited than is its normal enzymatic reaction. Polypeptide loop(s) closure may keep a proton from His159 interacting with the substrate phosphoryl oxygen long enough to stabilize a carbanion intermediate.