How fast monoamine oxidases decompose adrenaline? Kinetics of isoenzymes A and B evaluated by empirical valence bond simulation

Brunner syndrome caused by point mutation explained by multiscale simulation of enzyme reaction

Brunner syndrome is a disorder characterized by intellectual disability and impulsive, aggressive behavior associated with deficient function of the monoamine oxidase A (MAO-A) enzyme. These symptoms (along with particularly high serotonin levels) have been reported in patients with two missense variants in MAO-A (p.R45W and p.E446K). Herein, we report molecular simulations of the rate-limiting step of MAO-A-catalyzed serotonin degradation for these variants. We found that the R45W mutation causes a 6000-fold slowdown of enzymatic function, whereas the E446K mutation causes a 450-fold reduction of serotonin degradation rate, both of which are practically equivalent to a gene knockout. In addition, we thoroughly compared the influence of enzyme electrostatics on the catalytic function of both the wild type MAO-A and the p.R45W variant relative to the wild type enzyme, revealing that the mutation represents a significant electrostatic perturbation that contributes to the barrier increase. Understanding genetic disorders is closely linked to understanding the associated chemical mechanisms, and our research represents a novel attempt to bridge the gap between clinical genetics and the underlying chemical physics.

The catalytic mechanism of the mitochondrial methylenetetrahydrofolate dehydrogenase/cyclohydrolase (MTHFD2)

Methylenetetrahydrofolate dehydrogenase/cyclohydrolase (MTHFD2) is a new drug target that is expressed in cancer cells but not in normal adult cells, which provides an Achilles heel to selectively kill cancer cells. Despite the availability of crystal structures of MTHFD2 in the inhibitor- and cofactor-bound forms, key information is missing due to technical limitations, including (a) the location of absolutely required Mg2+ ion, and (b) the substrate-bound form of MTHFD2. Using computational modeling and simulations, we propose that two magnesium ions are present at the active site whereby (i) Arg233, Asp225, and two water molecules coordinate MgA2+, while MgA2+ together with Arg233 stabilize the inorganic phosphate (Pi); (ii) Asp168 and three water molecules coordinate MgB2+, and MgB2+ further stabilizes Pi by forming a hydrogen bond with two oxygens of Pi; (iii) Arg201 directly coordinates the Pi; and (iv) through three water-mediated interactions, Asp168 contributes to the positioning and stabilization of MgA2+, MgB2+ and Pi. Our computational study at the empirical valence bond level allowed us also to elucidate the detailed reaction mechanisms. We found that the dehydrogenase activity features a proton-coupled electron transfer with charge redistribution connected to the reorganization of the surrounding water molecules which further facilitates the subsequent cyclohydrolase activity. The cyclohydrolase activity then drives the hydration of the imidazoline ring and the ring opening in a concerted way. Furthermore, we have uncovered that two key residues, Ser197/Arg233, are important factors in determining the cofactor (NADP+/NAD+) preference of the dehydrogenase activity. Our work sheds new light on the structural and kinetic framework of MTHFD2, which will be helpful to design small molecule inhibitors that can be used for cancer treatment.

Pairing structural reconstruction with catalytic competence to evaluate the mechanisms of key enzymes in the folate-mediated one-carbon pathway

Mammalian metabolism comprises a series of interlinking pathways that include two major cycles: the folate and methionine cycles. The folate-mediated metabolic cycle uses several oxidation states of tetrahydrofolate to carry activated one-carbon units to be readily used and interconverted within the cell. They are required for nucleotide synthesis, methylation and metabolism, and particularly for proliferation of cancer cells. Based on the latest progress in genome-wide CRISPR loss-of-function viability screening of 789 cell lines, we focus on the most cancer-dependent enzymes in this pathway, especially those that are hyperactivated in cancer, to provide new insight into the chemical basis for cancer drug development. Since the complete 3D structure of several of these enzymes of the one-carbon pathway in their active form are not available, we used homology modelling integrated with the interpretation of the reaction mechanism. In addition, have reconstructed the most likely scenario for the reactions taking place paired with their catalytic competence that provides a testable framework for this pathway.

Why Monoamine Oxidase B Preferably Metabolizes N-Methylhistamine over Histamine: Evidence from the Multiscale Simulation of the Rate-Limiting Step

Histamine levels in the human brain are controlled by rather peculiar metabolic pathways. In the first step, histamine is enzymatically methylated at its imidazole N^τ atom, and the produced N-methylhistamine undergoes an oxidative deamination catalyzed by monoamine oxidase B (MAO-B), as is common with other monoaminergic neurotransmitters and neuromodulators of the central nervous system. The fact that histamine requires such a conversion prior to oxidative deamination is intriguing since MAO-B is known to be relatively promiscuous towards monoaminergic substrates; its in-vitro oxidation of N-methylhistamine is about 10 times faster than that for histamine, yet this rather subtle difference appears to be governing the decomposition pathway. This work clarifies the MAO-B selectivity toward histamine and N-methylhistamine by multiscale simulations of the rate-limiting hydride abstraction step for both compounds in the gas phase, in aqueous solution, and in the enzyme, using the established empirical valence bond methodology, assisted by gas-phase density functional theory (DFT) calculations. The computed barriers are in very good agreement with experimental kinetic data, especially for relative trends among systems, thereby reproducing the observed MAO-B selectivity. Simulations clearly demonstrate that solvation effects govern the reactivity, both in aqueous solution as well as in the enzyme although with an opposing effect on the free energy barrier. In the aqueous solution, the transition-state structure involving histamine is better solvated than its methylated analog, leading to a lower barrier for histamine oxidation. In the enzyme, the higher hydrophobicity of N-methylhistamine results in a decreased number of water molecules at the active side, leading to decreased dielectric shielding of the preorganized catalytic electrostatic environment provided by the enzyme. This renders the catalytic environment more efficient for N-methylhistamine, giving rise to a lower barrier relative to histamine. In addition, the transition state involving N-methylhistamine appears to be stabilized by the surrounding nonpolar residues to a larger extent than with unsubstituted histamine, contributing to a lower barrier with the former.

Quantum chemical (QM:MM) investigation of the mechanism of enzymatic reaction of tryptamine and N,N-dimethyltryptamine with monoamine oxidase A

The endogenous psychedelic (mind-altering) N,N-dimethyltryptamine (DMT) molecule has an important role in tissue protection, regeneration, and immunity via sigma-1 receptor activation as its natural ligand. The immunologic properties of DMT suggest this biogenic compound should be investigated thoroughly in other aspects as well. In our in silico project, we examined the metabolism of DMT and its primary analogue, the tryptamine (T), by the monoamine oxidase (MAO) flavoenzyme. MAO has two isoforms, MAO-A and MAO-B. MAOs perform the oxidation of various monoamines by their flavin adenine dinucleotide (FAD) cofactor. Two-layer QM:MM calculations at the ONIOM(M06-2X/6-31++G(d,p):UFF=QEq) level were performed including the whole enzyme to explore the potential energy surface (PES) of the reactions. Our findings reinforced that a hybrid mechanism, a mixture of pure H⁺ and H^- transfer pathways, describes precisely the rate-determining step of amine oxidation as suggested by earlier works. Additionally, our results show that the oxidation of tertiary amine DMT requires a lower activation barrier than the primary amine T. This may reflect a general rule, thus we recommend further investigations. Furthermore, we demonstrated that at pH 7.4 the protonated form of these substrates enter the enzyme. As the deprotonation of substrates is crucial, we presumed protonated cofactor, FADH⁺, may form. Surprisingly, the activation barriers are much lower compared to FAD with both substrates. Therefore, we suggest further investigations in this direction.

Adrenaline induces calcium signal in astrocytes and vasoconstriction via activation of monoamine oxidase

Adrenaline or epinephrine is a hormone playing an important role in physiology. It is produced de-novo in the brain in very small amounts compared to other catecholamines, including noradrenaline. Although the effects of adrenaline on neurons have been extensively studied, much less is known about the action of this hormone on astrocytes. Here, we studied the effects of adrenaline on astrocytes in primary co-culture of neurons and astrocytes. Application of adrenaline induced calcium signal in both neurons and astrocytes, but only in neurons this effect was dependent on α- and β-receptor antagonists. The effects of adrenaline on astrocytes were less dependent on adrenoreceptors: the antagonist carvedilol had only moderate effect on the calcium signal and the agonist of adrenoreceptors methoxamine induced a signal only in small proportion of the cells. We found that adrenaline in astrocytes activates phospholipase C and subsequent release of calcium from the endoplasmic reticulum. Calcium signal in astrocytes is initiated by the metabolism of adrenaline by the monoamine oxidase (MAO), which activates reactive oxygen species production and induces lipid peroxidation. Inhibitor of MAO selegiline inhibited both adrenaline-induced calcium signal in astrocytes and the vasoconstriction that indicates an important role for monoamine oxidase in adrenaline-induced signalling and function.

How Monoamine Oxidase A Decomposes Serotonin: An Empirical Valence Bond Simulation of the Reactive Step

The enzyme-catalyzed degradation of the biogenic amine serotonin is an essential regulatory mechanism of its level in the human organism. In particular, monoamine oxidase A (MAO A) is an important flavoenzyme involved in the metabolism of monoamine neurotransmitters. Despite extensive research efforts, neither the catalytic nor the inhibition mechanisms of MAO enzymes are currently fully understood. In this article, we present the quantum mechanics/molecular mechanics simulation of the rate-limiting step for the serotonin decomposition, which consists of hydride transfer from the serotonin methylene group to the N5 atom of the flavin moiety. Free-energy profiles of the reaction were computed by the empirical valence bond method. Apart from the enzymatic environment, the reference reaction in the gas phase was also simulated, facilitating the estimation of the catalytic effect of the enzyme. The calculated barrier for the enzyme-catalyzed reaction of 14.82 ± 0.81 kcal mol^-1 is in good agreement with the experimental value of 16.0 kcal mol^-1, which provides strong evidence for the validity of the proposed hydride-transfer mechanism. Together with additional experimental and computational work, the results presented herein contribute to a deeper understanding of the catalytic mechanism of MAO A and flavoenzymes in general, and in the long run, they should pave the way toward applications in neuropsychiatry.

Hydride Abstraction as the Rate-Limiting Step of the Irreversible Inhibition of Monoamine Oxidase B by Rasagiline and Selegiline: A Computational Empirical Valence Bond Study

Monoamine oxidases (MAOs) catalyze the degradation of a very broad range of biogenic and dietary amines including many neurotransmitters in the brain, whose imbalance is extensively linked with the biochemical pathology of various neurological disorders, and are, accordingly, used as primary pharmacological targets to treat these debilitating cognitive diseases. Still, despite this practical significance, the precise molecular mechanism underlying the irreversible MAO inhibition with clinically used propargylamine inhibitors rasagiline and selegiline is still not unambiguously determined, which hinders the rational design of improved inhibitors devoid of side effects current drugs are experiencing. To address this challenge, we present empirical valence bond QM/MM simulations of the rate-limiting step of the MAO inhibition involving the hydride anion transfer from the inhibitor α-carbon onto the N5 atom of the flavin adenin dinucleotide (FAD) cofactor. The proposed mechanism is strongly supported by the obtained free energy profiles, which confirm a higher reactivity of selegiline over rasagiline, while the calculated difference in the activation Gibbs energies of ΔΔG^‡ = 3.1 kcal mol^-1 is found to be in very good agreement with that from the measured literature k_inact values that predict a 1.7 kcal mol^-1 higher selegiline reactivity. Given the similarity with the hydride transfer mechanism during the MAO catalytic activity, these results verify that both rasagiline and selegiline are mechanism-based irreversible inhibitors and offer guidelines in designing new and improved inhibitors, which are all clinically employed in treating a variety of neuropsychiatric and neurodegenerative conditions.

Homology modeling of human GABA-AT and devise some novel and potent inhibitors via computer-aided drug design techniques

γ-aminobutyric acid aminotransferase (GABA-AT) is a pyridoxal 5'-phosphate (PLP)-dependent enzyme which degrades γ-aminobutyric (GABA) in the brain. GABA is an important inhibitory neurotransmitter that plays important neurological roles in the brain. Therefore, GABA-AT is an important drug target which regulates the GABA level. Novel and potent drug development to inhibit GABA-AT is still very challenging task. In this study, we aimed to devise novel and potent inhibitors against GABA-AT using computer-aided drug design (CADD) tools. However, the human GABA-AT crystal structure is not available yet, and we built the 3D structure of human GABA-AT based on the crystal structure of pig's liver (Sus Scrofa) enzyme as a template. The generated model was validated with numerous tools such as ProSA and PROCHECK. A set of selected well-known inhibitors have been tested against the modeled GABA-AT. Molecular docking studies have been accomplished via application of Genetic Optimization for Ligand Docking (GOLD), Vina and Autodock 4.2 software to search for potent inhibitors. The best two candidate inhibitors have been computationally examined for absorption, distribution, metabolism, elimination and toxicity descriptors (ADMET) and Lipinski's rule of 5. Lastly, molecular dynamics (MD) simulations were carried out to inspect the ligands' binding mode and stability of the active site of human GABA-AT over time. The top ranked ligands exhibited reliable stability throughout the MD simulation. The selected compounds are promising candidates and might be tested experimentally for the inhibition of human GABA-AT enzyme. Communicated by Ramaswamy H. Sarma.

Combined computational and experimental studies on cysteine-sulfadiazine adduct formation

The electrochemical characterization of sulfadiazine-cysteine (SD-CYS) adduct formation was performed in phosphate buffer (pH 7) on the basis of voltammetric current and peak potential measurements. Due to the association of cysteine with sulfadiazine, the reduction peak currents of mercuric and mercurous cysteine thiolates decreased and their peak potentials simultaneously shifted to less negative potentials. By using the current changes of mercurous cysteine thiolate, it was determined that cysteine and sulfadiazine are associated with a 1:1 stoichiometry with a conditional association constant of 1.99 ×104 M-1 . In addition to experimental studies, a computational approach was carried out to study the geometrical parameters, electron densities, and UV-Vis absorption spectra of sulfadiazine and SDCYS adduct in water. Calculated (B3LYP/6-311++G(d,p) level) and experimental UV-Vis absorption spectra of the compounds were found to be in good agreement in water. The computational study suggests that cysteine bound to the C(5) on the pyrimidine ring via SH-group nucleophilic attack. Computational results reveal that sulfadiazine and its derivatives effectively bind cysteine and may lead to new molecules/drugs to target cysteine.

Nuclear quantum effects in enzymatic reactions: simulation of the kinetic isotope effect of phenylethylamine oxidation catalyzed by monoamine oxidase A

By using computational techniques for quantizing nuclear motion one can accurately reproduce kinetic isotope effect of enzymatic reactions, as demonstrated for phenylethylamine oxidation catalyzed by the monoamine oxidase A enzyme.

Monoamine oxidase-B inhibitors as potential neurotherapeutic agents: An overview and update

Monoamine oxidase (MAO) inhibitors have made significant contributions and remain an indispensable approach of molecular and mechanistic diversity for the discovery of antineurodegenerative drugs. However, their usage has been hampered by nonselective and/or irreversible action which resulted in drawbacks like liver toxicity, cheese effect, and so forth. Hence, the search for selective MAO inhibitors (MAOIs) has become a substantial focus in current drug discovery. This review summarizes our current understanding on MAO-A/MAO-B including their structure, catalytic mechanism, and biological functions with emphases on the role of MAO-B as a potential therapeutic target for the development of medications treating neurodegenerative disorders. It also highlights the recent developments in the discovery of potential MAO-B inhibitors (MAO-BIs) belonging to diverse chemical scaffolds, arising from intensive chemical-mechanistic and computational studies documented during past 3 years (2015-2018), with emphases on their potency and selectivity. Importantly, readers will gain knowledge of various newly established MAO-BI scaffolds and their development potentials. The comprehensive information provided herein will hopefully accelerate ideas for designing novel selective MAO-BIs with superior activity profiles and critical discussions will inflict more caution in the decision-making process in the MAOIs discovery.

Dopamine Autoxidation Is Controlled by Acidic pH

We studied the reaction mechanism of dopamine autoxidation using quantum chemical methods. Unlike other biogenic amines important in the central nervous system, dopamine and noradrenaline are capable of undergoing a non-enzymatic autoxidative reaction giving rise to a superoxide anion that further decomposes to reactive oxygen species. The reaction in question, which takes place in an aqueous solution, is as such not limited to the mitochondrial membrane where scavenging enzymes such as catalase and superoxide dismutase are located. With the experimental rate constant of 0.147 s⁻¹, the dopamine autoxidation reaction is comparably as fast as the monoamine oxidase B catalyzed dopamine decomposition with a rate constant of 1 s⁻¹. By using quantum chemical calculations, we demonstrated that the rate-limiting step is the formation of a hydroxide ion from a water molecule, which attacks the amino group that enters intramolecular Michael addition, giving rise to a pharmacologically inert aminochrome. We have shown that for dopamine stability on a time scale of days, it is essential that the pH value of the synaptic vesicle interior is acidic. The pathophysiologic correlates of the results are discussed in the context of Parkinson's disease as well as the pathology caused by long-term amphetamine and cocaine administration.

Why does the Y326I mutant of monoamine oxidase B decompose an endogenous amphetamine at a slower rate than the wild type enzyme? Reaction step elucidated by multiscale molecular simulations

This work investigates the Y326I point mutation effect on the kinetics of oxidative deamination of phenylethylamine (PEA) catalyzed by the monoamine oxidase B (MAO B) enzyme. PEA is a neuromodulator capable of affecting the plasticity of the brain and is responsible for the mood enhancing effect caused by physical exercise. Due to a similar functionality, PEA is often regarded as an endogenous amphetamine. The rate limiting step of the deamination was simulated at the multiscale level, employing the Empirical Valence Bond approach for the quantum treatment of the involved valence states, whereas the environment (solvated protein) was represented with a classical force field. A comparison of the reaction free energy profiles delivered by simulation of the reaction in the wild type MAO B and its Y326I mutant yields an increase in the barrier by 1.06 kcal mol^-1 upon mutation, corresponding to a roughly 6-fold decrease in the reaction rate. This is in excellent agreement with the experimental kinetic studies. Inspection of simulation trajectories reveals possible sources of the point mutation effect, namely vanishing favorable electrostatic interactions between PEA and a Tyr326 side chain and an increased amount of water molecules at the active site due to the replacement of tyrosine by a less spacious isoleucine residue, thereby increasing the dielectric shielding of the catalytic environment provided by the enzyme.