Investigation of the catalytic mechanism of Sir2 enzyme with QM/MM approach: SN1 vs SN2?

Exponential averaging versus umbrella sampling for computing the QM/MM free energy barrier of the initial step of the desuccinylation reaction catalyzed by sirtuin 5

The computational characterization of enzymatic reactions poses a great challenge which arises from the high dimensional and often rough potential energy surfaces commonly explored by static QM/MM methods such as adiabatic mapping (AM). The present study highlights the difficulties in estimating free energy barriers via exponential averaging over AM pathways. Based on our previous study [von der Esch et al., J. Chem. Theory Comput., 2019, 15, 6660-6667], where we analyzed the first reaction step of the desuccinylation reaction catalyzed by human sirtuin 5 (SIRT5) by means of QM/MM adiabatic mapping and machine learning, we use, here, umbrella sampling to compute the free energy profile of the initial reaction step. The computational investigations show that the initial step of the desuccinylation reaction proceeds via an S_N2-type reaction mechanism in SIRT5, suggesting that the first step of the deacylation reactions catalyzed by sirtuins is highly conserved. In addition, the direct comparison of the extrapolated free energy barrier from minimal energy paths and the computed free energy path from umbrella sampling further underlines the importance of extensive sampling.

Finding Reactive Configurations: A Machine Learning Approach for Estimating Energy Barriers Applied to Sirtuin 5

Sirtuin 5 is a class III histone deacetylase that, unlike its classification, mainly catalyzes desuccinylation and demanoylation reactions. It is an interesting drug target that we use here to test new ideas for calculating reaction pathways of large molecular systems such as enzymes. A major issue with most schemes (e.g., adiabatic mapping) is that the resulting activation barrier height heavily depends on the chosen educt conformation. This makes the selection of the initial structure decisive for the success of the characterization. Here, we apply machine learning to a large number of molecular dynamics frames and potential energy barriers obtained by quantum mechanics/molecular mechanics calculations in order to identify (1) suitable start-conformations for reaction path calculations and (2) structural features relevant for the first step of the desuccinylation reaction catalyzed by Sirtuin 5. The latter generally aids the understanding of reaction mechanisms and important interactions in active centers. Using our novel approach, we found eleven key features that govern the reactivity. We were able to estimate reaction barriers with a mean absolute error of 3.6 kcal/mol and identified reactive configurations.

Unraveling the Unbinding Pathways of Products Formed in Catalytic Reactions Involved in SIRT1-3: A Random Acceleration Molecular Dynamics Simulation Study

Sirtuins are a family of nicotinamide adenine dinucleotide (NAD⁺)-dependent enzymes, which undergo robust deacetylase activity, resulting in the production of nicotinamide. It is well known that nicotinamide, which is one of the products, can also act as an inhibitor for further deacetylation process by forming NAD⁺ again. Hence, the removal of nicotinamide from sirtuins is a demanding process, and the mechanistic understanding of the process remains elusive. In this investigation, we have made an attempt to unravel the unbinding pathways of nicotinamide from SIRT1, SIRT2, and SIRT3 (SIRT1-3) using Random Acceleration Molecular Dynamics (RAMD) Simulations, and we have successfully identified various unbinding channels. The selectivity of the egression channel is determined by using a thorough analysis of the frequency of egression trajectories. Similarly, various inhibitors have been docked with the active sites of SIRT1-3, and their egression pathways have been investigated to understand whether they follow the same egression pathway as that of nicotinamide. The residues that are responsible for the unbinding pathways have been determined from the analysis of RAMD trajectories. From these results, it is clear that phenylalanine and histidine residues play major roles in the egression of inhibitors. Additionally, the key residues Leu, Pro, Met, Phe, Tyr, and Ile are found to control the release by acting as gateway residues. The role of these residues from different egression channels has been studied by carrying out mutations with alanine residue. This is the first report on sirtuins, which demonstrates the novel unbinding pathways for nicotinamide/inhibitors. This work provides new insights for developing more promising SIRT1-3 inhibitors.

Biophysical characterization of hit compounds for mechanism-based enzyme activation

Across all families of enzymes, only a dozen or so distinct classes of non-natural small molecule activators have been characterized, with only four known modes of activation among them. All of these modes of activation rely on naturally evolved binding sites that trigger global conformational changes. Among the enzymes that are of greatest interest for small molecule activation are the seven sirtuin enzymes, nicotinamide adenine dinucleotide (NAD⁺)-dependent protein deacylases that play a central role in the regulation of healthspan and lifespan in organisms ranging from yeast to mammals. However, there is currently no understanding of how to design sirtuin-activating compounds beyond allosteric activators of SIRT1-catalyzed reactions that are limited to particular substrates. Here, we introduce a general mode of sirtuin activation that is distinct from the known modes of enzyme activation. Based on the conserved mechanism of sirtuin-catalyzed deacylation reactions, we establish biophysical properties of small molecule modulators that can in principle result in enzyme activation for diverse sirtuins and substrates. Building upon this framework, we propose strategies for the identification, characterization and evolution of hits for mechanism-based enzyme activating compounds.

The chemical biology of sirtuins

This article reviews the tremendous accomplishments achieved during the past few years in the field of chemical biology for the physiologically and therapeutically important sirtuin family of N^ε-acyl-lysine deacylase enzymes.

Mechanism of inhibition of the human sirtuin enzyme SIRT3 by nicotinamide: computational and experimental studies

Sirtuins are key regulators of many cellular functions including cell growth, apoptosis, metabolism, and genetic control of age-related diseases. Sirtuins are themselves regulated by their cofactor nicotinamide adenine dinucleotide (NAD+) as well as their reaction product nicotinamide (NAM), the physiological concentrations of which vary during the process of aging. Nicotinamide inhibits sirtuins through the so-called base exchange pathway, wherein rebinding of the reaction product to the enzyme accelerates the reverse reaction. We investigated the mechanism of nicotinamide inhibition of human SIRT3, the major mitochondrial sirtuin deacetylase, in vitro and in silico using experimental kinetic analysis and Molecular Mechanics-Poisson Boltzmann/Generalized Born Surface Area (MM-PB(GB)SA) binding affinity calculations with molecular dynamics sampling. Through experimental kinetic studies, we demonstrate that NAM inhibition of SIRT3 involves apparent competition between the inhibitor and the enzyme cofactor NAD+, contrary to the traditional characterization of base exchange as noncompetitive inhibition. We report a model for base exchange inhibition that relates such kinetic properties to physicochemical properties, including the free energies of enzyme-ligand binding, and estimate the latter through the first reported computational binding affinity calculations for SIRT3:NAD+, SIRT3:NAM, and analogous complexes for Sir2. The computational results support our kinetic model, establishing foundations for quantitative modeling of NAD+/NAM regulation of mammalian sirtuins during aging and the computational design of sirtuin activators that operate through alleviation of base exchange inhibition.

Theoretical insights into catalytic mechanism of protein arginine methyltransferase 1

Protein arginine methyltransferase 1 (PRMT1), the major arginine asymmetric dimethylation enzyme in mammals, is emerging as a potential drug target for cancer and cardiovascular disease. Understanding the catalytic mechanism of PRMT1 will facilitate inhibitor design. However, detailed mechanisms of the methyl transfer process and substrate deprotonation of PRMT1 remain unclear. In this study, we present a theoretical study on PRMT1 catalyzed arginine dimethylation by employing molecular dynamics (MD) simulation and quantum mechanics/molecular mechanics (QM/MM) calculation. Ternary complex models, composed of PRMT1, peptide substrate, and S-adenosyl-methionine (AdoMet) as cofactor, were constructed and verified by 30-ns MD simulation. The snapshots selected from the MD trajectory were applied for the QM/MM calculation. The typical SN2-favored transition states of the first and second methyl transfers were identified from the potential energy profile. Deprotonation of substrate arginine occurs immediately after methyl transfer, and the carboxylate group of E144 acts as proton acceptor. Furthermore, natural bond orbital analysis and electrostatic potential calculation showed that E144 facilitates the charge redistribution during the reaction and reduces the energy barrier. In this study, we propose the detailed mechanism of PRMT1-catalyzed asymmetric dimethylation, which increases insight on the small-molecule effectors design, and enables further investigations into the physiological function of this family.

Sirtuin Deacetylation Mechanism and Catalytic Role of the Dynamic Cofactor Binding Loop

Sirtuins constitute a novel family of protein deacetylases and play critical roles in epigenetics, cell death, and metabolism. In spite of numerous experimental studies, the key and most complicated stage of its NAD+-dependent catalytic mechanism remains to be elusive. Herein by employing Born-Oppenheimer ab initio QM/MM molecular dynamics simulations, a state-of-the-art computational approach to study enzyme reactions, we have characterized the complete deacetylation mechanism for a sirtuin enzyme, determined its multistep free-energy reaction profile, and elucidated essential catalytic roles of the conserved dynamic cofactor binding loop. These new detailed mechanistic insights would facilitate the design of novel mechanism-based sirtuin modulators.

Modeling methods for studying post-translational and transcriptional modifying enzymes

Biological catalysis is a complex chemical process that involves not only electronic reorganization in the substrate but also the reorganization of the catalyst. This complexity is even larger in the case of post-transcriptional and post-translational modifications which may involve the interaction between two biomacromolecules. However, the development over the past decades of new computational methods and strategies is offering a detailed molecular picture of the catalytic event and a deep understanding of the mechanisms of chemical reactions in biological environments. Here we review the efforts made in the last years to model catalysis in post-transcriptional and post-translational processes. We stress on the advantages and problems of the different computational strategies and their applicability in different cases.

Investigation of the acetylation mechanism by GCN5 histone acetyltransferase

The histone acetylation of post-translational modification can be highly dynamic and play a crucial role in regulating cellular proliferation, survival, differentiation and motility. Of the enzymes that mediate post-translation modifications, the GCN5 of the histone acetyltransferase (HAT) proteins family that add acetyl groups to target lysine residues within histones, has been most extensively studied. According to the mechanism studies of GCN5 related proteins, two key processes, deprotonation and acetylation, must be involved. However, as a fundamental issue, the structure of hGCN5/AcCoA/pH3 remains elusive. Although biological experiments have proved that GCN5 mediates the acetylation process through the sequential mechanism pathway, a dynamic view of the catalytic process and the molecular basis for hGCN5/AcCoA/pH3 are still not available and none of theoretical studies has been reported to other related enzymes in HAT family. To explore the molecular basis for the catalytic mechanism, computational approaches including molecular modeling, molecular dynamic (MD) simulation and quantum mechanics/molecular mechanics (QM/MM) simulation were carried out. The initial hGCN5/AcCoA/pH3 complex structure was modeled and a reasonable snapshot was extracted from the trajectory of a 20 ns MD simulation, with considering post-MD analysis and reported experimental results. Those residues playing crucial roles in binding affinity and acetylation reaction were comprehensively investigated. It demonstrated Glu80 acted as the general base for deprotonation of Lys171 from H3. Furthermore, the two-dimensional QM/MM potential energy surface was employed to study the sequential pathway acetylation mechanism. Energy barriers of addition-elimination reaction in acetylation obtained from QM/MM calculation indicated the point of the intermediate ternary complex. Our study may provide insights into the detailed mechanism for acetylation reaction of GCN5, and has important implications for the discovery of regulators against GCN5 enzymes and related HAT family enzymes.

SIRT2 ablation has no effect on tubulin acetylation in brain, cholesterol biosynthesis or the progression of Huntington's disease phenotypes in vivo

Huntington's disease (HD) is a devastating neurodegenerative disorder for which there are no disease-modifying treatments. The molecular pathogenesis of HD is complex and many mechanisms and cellular processes have been proposed as potential sites of therapeutic intervention. However, prior to embarking on drug development initiatives, it is essential that therapeutic targets can be validated in mammalian models of HD. Previous studies in invertebrate and cell culture HD models have suggested that inhibition of SIRT2 could have beneficial consequences on disease progression. SIRT2 is a NAD⁺-dependent deacetylase that has been proposed to deacetylate α-tubulin, histone H4 K16 and to regulate cholesterol biogenesis – a pathway which is dysregulated in HD patients and HD mouse models. We have utilized mice in which SIRT2 has been reduced or ablated to further explore the function of SIRT2 and to assess whether SIRT2 loss has a beneficial impact on disease progression in the R6/2 mouse model of HD. Surprisingly we found that reduction or loss of SIRT2 had no effect on the acetylation of α-tubulin or H4K16 or on cholesterol biosynthesis in the brains of wild type mice. Equally, genetic reduction or ablation of SIRT2 had no effect on HD progression as assessed by a battery of physiological and behavioural tests. Furthermore, we observed no change in aggregate load or levels of soluble mutant huntingtin transprotein. Intriguingly, neither the constitutive genetic loss nor acute pharmacological inhibition of SIRT2 affected the expression of cholesterol biosynthesis enzymes in the context of HD. Therefore, we conclude that SIRT2 inhibition does not modify disease progression in the R6/2 mouse model of HD and SIRT2 inhibition should not be prioritised as a therapeutic option for HD.

Catalytic mechanism investigation of lysine-specific demethylase 1 (LSD1): a computational study

Lysine-specific demethylase 1 (LSD1), the first identified histone demethylase, is a flavin-dependent amine oxidase which specifically demethylates mono- or dimethylated H3K4 and H3K9 via a redox process. It participates in a broad spectrum of biological processes and is of high importance in cell proliferation, adipogenesis, spermatogenesis, chromosome segregation and embryonic development. To date, as a potential drug target for discovering anti-tumor drugs, the medical significance of LSD1 has been greatly appreciated. However, the catalytic mechanism for the rate-limiting reductive half-reaction in demethylation remains controversial. By employing a combined computational approach including molecular modeling, molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics (QM/MM) calculations, the catalytic mechanism of dimethylated H3K4 demethylation by LSD1 was characterized in details. The three-dimensional (3D) model of the complex was composed of LSD1, CoREST, and histone substrate. A 30-ns MD simulation of the model highlights the pivotal role of the conserved Tyr761 and lysine-water-flavin motif in properly orienting flavin adenine dinucleotide (FAD) with respect to substrate. The synergy of the two factors effectively stabilizes the catalytic environment and facilitated the demethylation reaction. On the basis of the reasonable consistence between simulation results and available mutagenesis data, QM/MM strategy was further employed to probe the catalytic mechanism of the reductive half-reaction in demethylation. The characteristics of the demethylation pathway determined by the potential energy surface and charge distribution analysis indicates that this reaction belongs to the direct hydride transfer mechanism. Our study provides insights into the LSD1 mechanism of reductive half-reaction in demethylation and has important implications for the discovery of regulators against LSD1 enzymes.