Revealing Monoamine Oxidase B Catalytic Mechanisms by Means of the Quantum Chemical Cluster Approach

Why Does Monoamine Oxidase (MAO) Catalyze the Oxidation of Some Tetrahydropyridines?

Results pertaining to the mechanism of the oxidation of the tertiary amine 1-methyl-4-(1-methyl-1-H-pyrrol-2-yl)-1,2,3,6-tetrahydropyridine (MMTP, a close analog of the Parkinsonism inducing compound MPTP) by 3-methyllumiflavin (3MLF), a chemical model for the FAD cofactor of monoamine oxidase, are reported. MMTP and related compounds are among the few tertiary amines that are monoamine oxidase B (MAO-B) substrates. The MMTP/3MLF reaction is catalytic in the presence of O2 and the results under anaerobic conditions strongly suggest the involvement of radical intermediates, consistent with a single electron transfer mechanism. These observations support a new hypothesis to explain the MAO-catalyzed oxidations of amines. In general, electron transfer is thermodynamically unfavorable, and as a result, most 1° and 2° amines react via one of the currently accepted polar pathways. Steric constraints prevent 3° amines from reacting via a polar pathway. Those select 3° amines that are MAO substrates possess certain structural features (e. g., a C-H bond that is α- both to nitrogen and a C=C) that dramatically lower the pKa of the corresponding radical cation. Consequently, the thermodynamically unfavorable electron transfer equilibrium is driven towards products by an extremely favorable deprotonation step in the context of Le Chatelier's principle.

Cluster energy prediction based on multiple strategy fusion whale optimization algorithm and light gradient boosting machine

BACKGROUND

Clusters, a novel hierarchical material structure that emerges from atoms or molecules, possess unique reactivity and catalytic properties, crucial in catalysis, biomedicine, and optoelectronics. Predicting cluster energy provides insights into electronic structure, magnetism, and stability. However, the structure of clusters and their potential energy surface is exceptionally intricate. Searching for the global optimal structure (the lowest energy) among these isomers poses a significant challenge. Currently, modelling cluster energy predictions with traditional machine learning methods has several issues, including reliance on manual expertise, slow computation, heavy computational resource demands, and less efficient parameter tuning.

RESULTS

This paper introduces a predictive model for the energy of a gold cluster comprising twenty atoms (referred to as Au20 cluster). The model integrates the Multiple Strategy Fusion Whale Optimization Algorithm (MSFWOA) with the Light Gradient Boosting Machine (LightGBM), resulting in the MSFWOA-LightGBM model. This model employs the Coulomb matrix representation and eigenvalue solution methods for feature extraction. Additionally, it incorporates the Tent chaotic mapping, cosine convergence factor, and inertia weight updating strategy to optimize the Whale Optimization Algorithm (WOA), leading to the development of MSFWOA. Subsequently, MSFWOA is employed to optimize the parameters of LightGBM for supporting the energy prediction of Au20 cluster.

CONCLUSIONS

The experimental results show that the most stable Au20 cluster structure is a regular tetrahedron with the lowest energy, displaying tight and uniform atom distribution, high geometric symmetry. Compared to other models, the MSFWOA-LightGBM model excels in accuracy and correlation, with MSE, RMSE, and R2 values of 0.897, 0.947, and 0.879, respectively. Additionally, the MSFWOA-LightGBM model possesses outstanding scalability, offering valuable insights for material design, energy storage, sensing technology, and biomedical imaging, with the potential to drive research and development in these areas.

Computational Mechanistic Study of l-Aspartate Oxidase by ONIOM Method

l-Aspartate oxidase (Laspo) is responsible for the oxidation of l-aspartate into iminoaspartate using flavin as a cofactor. During this process flavin is reduced, and it can be reoxidized by either molecular oxygen or fumarate. The overall fold and the catalytic residues of Laspo are similar to succinate dehydrogenase and fumarate reductase. On the basis of deuterium kinetic isotope effects as well as other kinetic and structural data, it is proposed that the enzyme can catalyze the oxidation of l-aspartate through a mechanism similar to amino acid oxidases. It is suggested that a proton is removed from the α-amino group, while a hydride is transferred from C2 to flavin. It is also suggested that the hydride transfer is a rate-limiting step. However, there is still an ambiguity about the stepwise or concerted mechanism of hydride- and proton-transfer steps. In this study, we formulated some computational models to study the hydride-transfer mechanism using the crystal structure of Escherichia colil-aspartate oxidase in complexes with succinate. The calculations involved our own N-layered integrated molecular orbital and molecular mechanics method, and we evaluated the geometry and energetics of the hydride/proton-transfer processes while probing the roles of active site residues. Based on the calculations, it is concluded that proton- and hydride-transfer steps are decoupled, and a stepwise mechanism might be operative as opposed to the concerted one.

Computational Insights into β-Carboline Inhibition of Monoamine Oxidase A

Monoamine oxidases (MAOs) are an important group of enzymes involved in the degradation of neurotransmitters and their imbalanced mode of action may lead to the development of various neuropsychiatric or neurodegenerative disorders. In this work, we report the results of an in-depth computational study in which we performed a static and a dynamic analysis of a series of substituted β-carboline natural products, found mainly in roasted coffee and tobacco smoke, that bind to the active site of the MAO-A isoform. By applying molecular docking in conjunction with structure-based pharmacophores and molecular dynamics simulations coupled with dynamic pharmacophores, we extensively investigated the geometric aspects of MAO-A binding. To gain insight into the energetics of binding, we used the linear interaction energy (LIE) method and determined the key anchors that allow productive β-carboline binding to MAO-A. The results presented herein could be applied in the rational structure-based design and optimization of β-carbolines towards preclinical candidates that would target the MAO-A enzyme and would be applicable especially in the treatment of mental disorders such as depression.

Roles of selected non-P450 human oxidoreductase enzymes in protective and toxic effects of chemicals: review and compilation of reactions

This is an overview of the metabolic reactions of drugs, natural products, physiological compounds, and other (general) chemicals catalyzed by flavin monooxygenase (FMO), monoamine oxidase (MAO), NAD(P)H quinone oxidoreductase (NQO), and molybdenum hydroxylase enzymes (aldehyde oxidase (AOX) and xanthine oxidoreductase (XOR)), including roles as substrates, inducers, and inhibitors of the enzymes. The metabolism and bioactivation of selected examples of each group (i.e., drugs, “general chemicals,” natural products, and physiological compounds) are discussed. We identified a higher fraction of bioactivation reactions for FMO enzymes compared to other enzymes, predominately involving drugs and general chemicals. With MAO enzymes, physiological compounds predominate as substrates, and some products lead to unwanted side effects or illness. AOX and XOR enzymes are molybdenum hydroxylases that catalyze the oxidation of various heteroaromatic rings and aldehydes and the reduction of a number of different functional groups. While neither of these two enzymes contributes substantially to the metabolism of currently marketed drugs, AOX has become a frequently encountered route of metabolism among drug discovery programs in the past 10–15 years. XOR has even less of a role in the metabolism of clinical drugs and preclinical drug candidates than AOX, likely due to narrower substrate specificity.

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.

Predicting Effects of Site-Directed Mutagenesis on Enzyme Kinetics by QM/MM and QM Calculations: A Case of Glutamate Carboxypeptidase II

Quantum and molecular mechanics (QM/MM) and QM-only (cluster model) modeling techniques represent the two workhorses in mechanistic understanding of enzyme catalysis. One of the stringent tests for QM/MM and/or QM approaches is to provide quantitative answers to real-world biochemical questions, such as the effect of single-point mutations on enzyme kinetics. This translates into predicting the relative activation energies to 1-2 kcal·mol^-1 accuracy; such predictions can be used for the rational design of novel enzyme variants with desired/improved characteristics. Herein, we employ glutamate carboxypeptidase II (GCPII), a dizinc metallopeptidase, also known as the prostate specific membrane antigen, as a model system. The structure and activity of this major cancer antigen have been thoroughly studied, both experimentally and computationally, which makes it an ideal model system for method development. Its reaction mechanism is quite well understood: the reaction coordinate comprises a "tetrahedral intermediate" and two transition states and experimental activation Gibbs free energy of ∼17.5 kcal·mol^-1 can be inferred for the known k_cat ≈ 1 s^-1. We correlate experimental kinetic data (including the E424H variant, newly characterized in this work) for various GCPII mutants (k_cat = 8.6 × 10^-5 s^-1 to 2.7 s^-1) with the energy profiles calculated by QM/MM and QM-only (cluster model) approaches. We show that the near-quantitative agreement between the experimental values and the calculated activation energies (ΔH^⧧) can be obtained and recommend the combination of the two protocols: QM/MM optimized structures and cluster model (QM) energetics. The trend in relative activation energies is mostly independent of the QM method (DFT functional) used. Last but not least, a satisfactory correlation between experimental and theoretical data allows us to provide qualitative and fairly simple explanations of the observed kinetic effects which are thus based on a rigorous footing.

An electrostatic duel: subtle differences in the catalytic performance of monoamine oxidase A and B isoenzymes elucidated at the residue level using quantum computations

The origin of the immense catalytic power of enzymes remains one of the biggest unresolved questions in biochemistry, with electrostatics being one of the main contenders. Herein, we report results that not only confirm that electrostatics is the driving force behind enzyme catalysis, but also that it is capable of tuning subtle differences in the catalytic performance between structurally similar enzymes, as demonstrated using the example of isoenzymes, monoamine oxidases A and B. Using our own computationally efficient multiscale model [A. Prah, et al., ACS Catal., 2019, 9, 1231] we analyzed the rate-limiting step of the reaction between phenylethylamine and both isoenzymes and deduced that the electrostatic environment provided by isoenzyme B has a perceivably higher catalytic influence on all the considered parameters of the reaction (energy barrier, charge transfer, dipole moment, and HOMO-LUMO gap). This is in full agreement with the available experimental kinetic data and with our own simulations of the reaction in question. In-depth analysis of individual amino acid contributions of both isoenzymes to the barrier (based on the interaction between the electric field provided by the enzyme and the dipole moment of the reacting moiety) shows that the majority of the difference between the isoenzymes can be attributed to a small number of sizable differences between the aligned amino acid pairs, whereas in most of the pairs the difference in contribution to the barrier is vanishingly small. These results suggest that electrostatics largely controls the substrate selectivity of enzymes and validates our approach as being capable of discerning fine nuances in the selectivity of structurally related isoenzymes.

Computational Analysis of the Nicotine Oxidoreductase Mechanism by the ONIOM Method

Nicotine oxidoreductase (NicA2) is a monoamine oxidase (MAO)-based flavoenzyme that catalyzes the oxidation of S-nicotine into N-methylmyosmine. Due to its nanomolar binding affinity toward nicotine, it is seen as an ideal candidate for the treatment of nicotine addiction. Based on the crystal structure of the substrate-bound enzyme, hydrophobic interactions mainly govern the binding of the substrate in the active site through Trp108, Trp364, Trp427, and Leu217 residues. In addition, Tyr308 forms H-bonding with the pyridyl nitrogen of the substrate. Experimental and computational studies support the hydride transfer mechanism for MAO-based enzymes. In this mechanism, a hydride ion transfers from the substrate to the flavin cofactor. In this study, computational models involving the ONIOM method were formulated to study the hydride transfer mechanism based on the crystal structure of the enzyme-substrate complex. The geometry and energetics of the hydride transfer mechanism were analyzed, and the roles of active site residues were highlighted.

Alkylammonium Cation Affinities of Nitrogenated Organobases: The Roles of Hydrogen Bonding and Proton Transfer

Alkylammonium cation affinities of 64 nitrogen-containing organobases, as well as the respective proton transfer processes from the alkylammonium cations to the base, have been computed in the gas phase by using DFT methods. The guanidine bases show the highest proton transfer values (191.9-233 kJ mol^-1 ) whereas the cis-2,2'-biimidazole presents the largest affinity towards the alkylammonium cations (>200 kJ mol^-1 ) values. The resulting data have been compared with the experimentally reported proton affinities of the studied nitrogen-containing organobases revealing that the propensity of an organobase for the proton transfer process increases linearly with its proton affinity. This work can provide a tool for designing senors for bioactive compounds containing amino groups that are protonated at physiological pH.

Theoretical Study on the Mechanism of the Acylate Reaction of β-Lactamase

Using density functional theory and a cluster approach, we study the reaction potential surface and compute Gibbs free energies for the acylate reaction of β-lactamase with penicillin G, where the solvent effect is important and taken into consideration. Two reaction paths are investigated: one is a multi-step process with a rate-limit energy barrier of 19.1 kcal/mol, which is relatively small, and the reaction can easily occur; the other is a one-step process with a barrier of 45.0 kcal/mol, which is large and thus makes the reaction hard to occur. The reason why the two paths have different barriers is explained.

Binding mechanism of naringenin with monoamine oxidase - B enzyme: QM/MM and molecular dynamics perspective

The reduced level of dopamine at midbrain (substantia nigra) leads to Parkinson disease by the influence of monoamine oxidation process of monoamine oxidase B (MAO-B) enzyme. This disease mostly affects the aged people. Reports outline that the naringenin molecule acts as an inhibitor of MAO-B enzyme and it potentially prevents the development of PD. To elucidate the binding mechanism of naringenin with MAO-B, we performed the molecular docking, QM/MM and molecular dynamics (MD) simulations. The molecular docking results confirm that the naringenin strongly binds with the substrate binding site of MAO-B enzyme (-12.0 kcal/mol). The low values of RMSD, RMSF and Rg indicate that the naringenin - MAO-B complex is stable over the entire period of MD simulation. Naringenin forms strong interaction with the orient keeper residue Tyr326 and other binding site residues Leu171, Glu206 and these interactions were maintained throughout the MD simulation. It is also important to block the function of MAO-B enzyme. The QM/MM study coupled with the charge density analysis reveals the charge density distribution and the strength of intermolecular interactions of naringenin-MAO-B complex. The above results suggest that this molecule is a potential inhibitor of MAO-B enzyme.

Mechanistic study of L-6-hydroxynicotine oxidase by DFT and ONIOM methods

L-6-Hydroxynicotine oxidase (LHNO) is a member of monoamine oxidase (MAO) family and catalyzes conversion of (S)-6-hydroxynicotine to 6-hydroxypseudooxynicotine during bacterial degradation of nicotine. Recent studies indicated that the enzyme catalyzes oxidation of carbon-nitrogen bond instead of previously proposed carbon-carbon bond. Based on kinetics and mutagenesis studies, Asn166, Tyr311, and Lys287 as well as an active site water molecule have roles in the catalysis of the enzyme. A number of studies including experimental and computational methods support hydride transfer mechanism in MAO family as a common mechanism in which a hydride ion transfer from amine substrate to flavin cofactor is the rate-limiting step. In this study, we formulated computational models to study the hydride transfer mechanism using crystal structure of enzyme-substrate complex. The calculations involved ONIOM and DFT methods, and we evaluated the geometry and energetics of the hydride transfer process while probing the roles of active site residues. Based on the calculations involving hydride, radical, and polar mechanisms, it was concluded that hydride transfer mechanism is the only viable mechanism for LHNO.

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.

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.

Advances in the Treatment of Explicit Water Molecules in Docking and Binding Free Energy Calculations

Background:

The inclusion of direct effects mediated by water during the ligandreceptor recognition is a hot-topic of modern computational chemistry applied to drug discovery and development. Docking or virtual screening with explicit hydration is still debatable, despite the successful cases that have been presented in the last years. Indeed, how to select the water molecules that will be included in the docking process or how the included waters should be treated remain open questions.

Objective:

In this review, we will discuss some of the most recent methods that can be used in computational drug discovery and drug development when the effect of a single water, or of a small network of interacting waters, needs to be explicitly considered.

Results:

Here, we analyse the software to aid the selection, or to predict the position, of water molecules that are going to be explicitly considered in later docking studies. We also present software and protocols able to efficiently treat flexible water molecules during docking, including examples of applications. Finally, we discuss methods based on molecular dynamics simulations that can be used to integrate docking studies or to reliably and efficiently compute binding energies of ligands in presence of interfacial or bridging water molecules.

Conclusions:

Software applications aiding the design of new drugs that exploit water molecules, either as displaceable residues or as bridges to the receptor, are constantly being developed. Although further validation is needed, workflows that explicitly consider water will probably become a standard for computational drug discovery soon.

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.

Computational Insight into the Mechanism of the Irreversible Inhibition of Monoamine Oxidase Enzymes by the Antiparkinsonian Propargylamine Inhibitors Rasagiline and Selegiline

Monoamine oxidases (MAOs) are flavin adenine dinucleotide containing flavoenzymes that catalyze the degradation of a range of brain neurotransmitters, whose imbalance is extensively linked with the pathology of various neurological disorders. This is why MAOs have been the central pharmacological targets in treating neurodegeneration for more than 60 years. Still, despite this practical importance, the precise chemical mechanisms underlying the irreversible inhibition of the MAO B isoform with clinical drugs rasagiline (RAS) and selegiline (SEL) remained unknown. Here we employed a combination of MD simulations, MM-GBSA binding free energy evaluations, and QM cluster calculations to show the MAO inactivation proceeds in three steps, where, in the rate-limiting first step, FAD utilizes its N5 atom to abstracts a hydride anion from the inhibitor α-CH₂ group to ultimately give the final inhibitor-FAD adduct matching crystallographic data. The obtained free energy profiles reveal a lower activation energy for SEL by 1.2 kcal mol^-1 and a higher reaction exergonicity by 0.8 kcal mol^-1, with the former being in excellent agreement with experimental ΔΔG^‡_EXP = 1.7 kcal mol^-1, thus rationalizing its higher in vivo reactivity over RAS. The calculated ΔG_BIND energies confirm SEL binds better due to its bigger size and flexibility allowing it to optimize hydrophobic C-H···π and π···π interactions with residues throughout both of enzyme's cavities, particularly with FAD, Gln206 and four active site tyrosines, thus overcoming a larger ability of RAS to form hydrogen bonds that only position it in less reactive orientations for the hydride abstraction. Offered results elucidate structural determinants affecting the affinity and rates of the inhibition reaction that should be considered to cooperate when designing more effective compounds devoid of untoward effects, which are of utmost significance and urgency with the growing prevalence of brain diseases.

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.

On the practical aspects of characterising monoamine oxidase inhibition in vitro

The development of novel inhibitors of human monoamine oxidase enzymes with improved pharmacodynamic and pharmacokinetic profiles has, in the past, been hampered by limited access to enzyme, by assay protocols offering limited throughput, and by inappropriate analyses of kinetic data. More recently, high-level expression of human enzymes in yeast has facilitated thorough examinations of steady-state enzyme behaviour that have led to improvements in our understanding of the mathematical underpinnings of kinetic analyses of monoamine oxidases. However, with these improvements have come a realisation that to be useful, more data points across wider concentration ranges are required. In turn, many discontinuous assay approaches, such as those involving radiolabelled substrates or chromatographic separation of product from substrate, have been rendered somewhat obsolete. Justification for the use of a platereader-based approach to assess the effects of novel inhibitors on monamine oxidases is provided, along with details of experimental design optimised to address the unexpectedly complex kinetics followed by these enzymes. Potential sources of error are discussed, and comments provided on techniques that may enhance the quality of experimental data.

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.

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

This work scrutinizes kinetics of decomposition of adrenaline catalyzed by monoamine oxidase (MAO) A and B enzymes, a process controlling the levels of adrenaline in the central nervous system and other tissues. Experimental kinetic data for MAO A and B catalyzed decomposition of adrenaline are reported only in the form of the maximum reaction rate. Therefore, we estimated the experimental free energy barriers form the kinetic data of closely related systems using regression method, as was done in our previous study. By using multiscale simulation on the Empirical Valence Bond (EVB) level, we studied the chemical reactivity of the MAO A catalyzed decomposition of adrenaline and we obtained a value of activation free energy of 17.3 ± 0.4 kcal/mol. The corresponding value for MAO B is 15.7 ± 0.7 kcal/mol. Both values are in good agreement with the estimated experimental barriers of 16.6 and 16.0 kcal/mol for MAO A and MAO B, respectively. The fact that we reproduced the kinetic data and preferential catalytic effect of MAO B over MAO A gives additional support to the validity of the proposed hydride transfer mechanism. Furthermore, we demonstrate that adrenaline is preferably involved in the reaction in a neutral rather than in a protonated form due to considerably higher barriers computed for the protonated adrenaline substrate. The results are discussed in the context of chemical mechanism of MAO enzymes and possible applications of multiscale simulation to rationalize the effects of MAO activity on adrenaline level.

Multiscale simulation of monoamine oxidase catalyzed decomposition of phenylethylamine analogs

Phenylethylamine (PEA) is an endogenous amphetamine and its levels are increased by physical activity. As other biogenic monoamines, it is decomposed by monoamine oxidase (MAO) enzymes. The chemical mechanism of MAO, and flavoenzymes in general, is a subject of heated debate. We have previously shown that the rate-limiting step of MAO catalysis involves a hydride transfer from the substrate methylene group vicinal to the amino group to the N5 atom of the lumiflavin co-factor moiety. By using multiscale simulation on the Empirical Valence Bond (EVB) level, we studied the chemical reactivity of the monoamine oxidase B catalyzed decomposition of PEA and its two derivatives: p-chloro-β-methylphenylamine (p-CMP) and p-methoxy-β-methylphenethylamine (p-MMP). We calculated activation free energies of 17.1kcal/mol (PEA), 18.4kcal/mol (p-MMP) and 20.0kcal/mol (p-CMP), which are in excellent agreement with the experimental values of 16.7kcal/mol for PEA and 18.3kcal/mol for p-MMP, while the experimental value for p-CMP is not available. This gives strong support to the validity of our hydride transfer mechanism for both MAO A and B isoforms. The results are discussed in the context of the interplay between MAO point mutations and neuropsychiatric disorders.

Recent Trends in Quantum Chemical Modeling of Enzymatic Reactions

The quantum chemical cluster approach is a powerful method for investigating enzymatic reactions. Over the past two decades, a large number of highly diverse systems have been studied and a great wealth of mechanistic insight has been developed using this technique. This Perspective reviews the current status of the methodology. The latest technical developments are highlighted, and challenges are discussed. Some recent applications are presented to illustrate the capabilities and progress of this approach, and likely future directions are outlined.

Mechanistic Studies of an Amine Oxidase Derived from d-Amino Acid Oxidase

The flavoprotein d-amino acid oxidase has long served as a paradigm for understanding the mechanism of oxidation of amino acids by flavoproteins. Recently, a mutant d-amino acid oxidase (Y228L/R283G) that catalyzed the oxidation of amines rather than amino acids was described [Yasukawa, K., et al. (2014) Angew. Chem., Int. Ed. 53, 4428-4431]. We describe here the use of pH and kinetic isotope effects with (R)-α-methylbenzylamine as a substrate to determine whether the mutant enzyme utilizes the same catalytic mechanism as the wild-type enzyme. The effects of pH on the steady-state and rapid-reaction kinetics establish that the neutral amine is the substrate, while an active-site residue, likely Tyr224, must be uncharged for productive binding. There is no solvent isotope effect on the k_cat/K_m value for the amine, consistent with the neutral amine being the substrate. The deuterium isotope effect on the k_cat/K_m value is pH-independent, with an average value of 5.3, similar to values found with amino acids as substrates for the wild-type enzyme and establishing that there is no commitment to catalysis with this substrate. The k_cat/K_O2 value is similar to that seen with amino acids as the substrate, consistent with the oxidative half-reaction being unperturbed by the mutation and with flavin oxidation preceding product release. All of the data are consistent with the mutant enzyme utilizing the same mechanism as the wild-type enzyme, transfer of hydride from the neutral amine to the flavin.

What a Difference a Methyl Group Makes: The Selectivity of Monoamine Oxidase B Towards Histamine and N-Methylhistamine

Monoamine oxidase (MAO) enzymes 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. Although sharing around 70 % sequence identity, both MAO A and B isoforms differ in substrate affinities and inhibitor sensitivities. Inhibitors that act on MAO A are used to treat depression, due to their ability to raise serotonin concentrations, whereas MAO B inhibitors decrease dopamine degradation and improve motor control in patients with Parkinson disease. Despite this functional importance, the factors affecting MAO selectivity are poorly understood. Here, we used a combination of molecular dynamics (MD) simulations, molecular mechanics with Poisson-Boltzmann and surface area solvation (MM-PBSA) binding free energy evaluations, and quantum mechanical (QM) cluster calculations to address the unexpected, yet challenging MAO B selectivity for N-methylhistamine (NMH) over histamine (HIS), differing only in a single methyl group distant from the reactive ethylamino center. This study shows that a dominant selectivity contribution is offered by a lower activation free energy for NMH by 2.6 kcal mol^-1 , in excellent agreement with the experimental ΔΔG^≠_EXP =1.4 kcal mol^-1 , together with a more favorable reaction exergonicity and active-site binding. This study also confirms the hydrophobic nature of the MAO B active site and underlines the important role of Ile199, Leu171, and Leu328 in properly orienting substrates for the reaction.

The Structure of Murine N1-Acetylspermine Oxidase Reveals Molecular Details of Vertebrate Polyamine Catabolism

N¹-Acetylspermine oxidase (APAO) catalyzes the conversion of N¹-acetylspermine or N¹-acetylspermidine to spermidine or putrescine, respectively, with concomitant formation of N-acetyl-3-aminopropanal and hydrogen peroxide. Here we present the structure of murine APAO in its oxidized holo form and in complex with substrate. The structures provide a basis for understanding molecular details of substrate interaction in vertebrate APAO, highlighting a key role for an asparagine residue in coordinating the N¹-acetyl group of the substrate. We applied computational methods to the crystal structures to rationalize previous observations with regard to the substrate charge state. The analysis suggests that APAO features an active site ideally suited for binding of charged polyamines. We also reveal the structure of APAO in complex with the irreversible inhibitor MDL72527. In addition to the covalent adduct, a second MDL72527 molecule is bound in the active site. Binding of MDL72527 is accompanied by altered conformations in the APAO backbone. On the basis of structures of APAO, we discuss the potential for development of specific inhibitors.

Complete sampling of an enzyme reaction pathway: a lesson from gas phase simulations

With proper sampling strategy, convergence of free energy profiles of biomolecular reactions in the gas phase can be achieved in microseconds of simulation.

Empirical Valence Bond Simulations of the Hydride-Transfer Step in the Monoamine Oxidase A Catalyzed Metabolism of Noradrenaline

Monoamine oxidases (MAOs) A and B are flavoenzymes responsible for the metabolism of biogenic amines, such as dopamine, serotonin, and noradrenaline (NA), which is why they have been extensively implicated in the etiology and course of various neurodegenerative disorders and, accordingly, used as primary pharmacological targets to treat these debilitating cognitive diseases. The precise chemical mechanism through which MAOs regulate the amine concentration, which is vital for the development of novel inhibitors, is still not unambiguously determined in the literature. In this work, we present atomistic empirical valence bond simulations of the rate-limiting step of the MAO-A-catalyzed NA (norepinephrine) degradation, involving hydride transfer from the substrate α-methylene group to the flavin moiety of the flavin adenine dinucleotide prosthetic group, employing the full dimensionality and thermal fluctuations of the hydrated enzyme, with extensive configurational sampling. We show that MAO-A lowers the free energy of activation by 14.3 kcal mol^-1 relative to that of the same reaction in aqueous solution, whereas the calculated activation free energy of ΔG^‡ = 20.3 ± 1.6 kcal mol^-1 is found to be in reasonable agreement with the correlated experimental value of 16.5 kcal mol^-1. The results presented here strongly support the fact that both MAO-A and MAO-B isoforms function by the same hydride-transfer mechanism. We also considered a few point mutations of the "aromatic cage" tyrosine residue (Tyr444Phe, Tyr444Leu, Tyr444Trp, Tyr444His, and Tyr444Glu), and the calculated changes in the reaction barriers are in agreement with the experimental values, thus providing further support to the proposed mechanism.

ONIOM calculations on serotonin degradation by monoamine oxidase B: insight into the oxidation mechanism and covalent reversible inhibition

Monoamine oxidase (MAO) is an enzyme which catalyzes the oxidation of neurotransmitter amines and regulates their level. There are two forms of the enzyme with 70% similarity, known as MAO-A and MAO-B. MAO inhibitors are used in the treatment of neurological disorders such as depression, Parkinson's and Alzheimer's diseases. Therefore, understanding the chemical steps of MAO catalyzed amine oxidation is crucial for rational drug design. However, despite many experimental studies and recent computational efforts in the literature, the amine oxidation mechanism by MAO enzymes is still controversial. The polar nucleophilic mechanism and hydride transfer mechanisms are under debate in recent QM/MM studies. In this study, the serotonin oxidation mechanism by MAO was explored via the ONIOM (QM : QM) methodology at the M06-2X/6-31+G(d,p):PM6 level. A modified MAO mechanism involving a covalent reversible inhibition step via formation of flavin N5 ylide was proposed. This mechanism can be used to modulate the potency and reversibility of novel mechanism-based covalent inhibitors by intelligent modifications of the structure of the inhibitors. NBO donor-acceptor analysis confirms that the rate-determining αC-H cleavage step is a hybrid of hydride and proton transfer where hydride transfer dominates over the proton transfer. The functional role of covalent FAD was also investigated by calculating the activation energy of noncovalent FAD models where a 22 fold decrease in the rate of catalysis was predicted. Geometrical features imply that the function of the covalent bond in FAD might be to maintain the correct geometry and conformation for a more efficient catalysis.

Red emission fluorescent probes for visualization of monoamine oxidase in living cells

Here we report two novel red emission fluorescent probes for the highly sensitive and selective detection of monoamine oxidase (MAO) with large Stokes shift (227 nm). Both of the probes possess solid state fluorescence and can accomplish the identification of MAO on test papers. The probe MAO-Red-1 exhibited a detection limit down to 1.2 μg mL⁻¹ towards MAO-B. Moreover, the cleavage product was unequivocally conformedby HPLC and LCMS and the result was in accordance with the proposed oxidative deamination mechanism. The excellent photostability of MAO-Red-1 was proved both in vitro and in vivo through fluorescent kinetic experiment and laser exposure experiment of confocal microscopy, respectively. Intracellular experiments also confirmed the low cytotoxity and exceptional cell imaging abilities of MAO-Red-1. It was validated both in HeLa and HepG2 cells that MAO-Red-1 was capable of reporting MAO activity through the variation of fluorescence intensity.

The Use of Multiscale Molecular Simulations in Understanding a Relationship between the Structure and Function of Biological Systems of the Brain: The Application to Monoamine Oxidase Enzymes

HIGHLIGHTS Computational techniques provide accurate descriptions of the structure and dynamics of biological systems, contributing to their understanding at an atomic level.Classical MD simulations are a precious computational tool for the processes where no chemical reactions take place.QM calculations provide valuable information about the enzyme activity, being able to distinguish among several mechanistic pathways, provided a carefully selected cluster model of the enzyme is considered.Multiscale QM/MM simulation is the method of choice for the computational treatment of enzyme reactions offering quantitative agreement with experimentally determined reaction parameters.Molecular simulation provide insight into the mechanism of both the catalytic activity and inhibition of monoamine oxidases, thus aiding in the rational design of their inhibitors that are all employed and antidepressants and antiparkinsonian drugs. Aging society and therewith associated neurodegenerative and neuropsychiatric diseases, including depression, Alzheimer's disease, obsessive disorders, and Parkinson's disease, urgently require novel drug candidates. Targets include monoamine oxidases A and B (MAOs), acetylcholinesterase (AChE), butyrylcholinesterase (BChE), and various receptors and transporters. For rational drug design it is particularly important to combine experimental synthetic, kinetic, toxicological, and pharmacological information with structural and computational work. This paper describes the application of various modern computational biochemistry methods in order to improve the understanding of a relationship between the structure and function of large biological systems including ion channels, transporters, receptors, and metabolic enzymes. The methods covered stem from classical molecular dynamics simulations to understand the physical basis and the time evolution of the structures, to combined QM, and QM/MM approaches to probe the chemical mechanisms of enzymatic activities and their inhibition. As an illustrative example, the later will focus on the monoamine oxidase family of enzymes, which catalyze the degradation of amine neurotransmitters in various parts of the brain, the imbalance of which is associated with the development and progression of a range of neurodegenerative disorders. Inhibitors that act mainly on MAO A are used in the treatment of depression, due to their ability to raise serotonin concentrations, while MAO B inhibitors decrease dopamine degradation and improve motor control in patients with Parkinson disease. Our results give strong support that both MAO isoforms, A and B, operate through the hydride transfer mechanism. Relevance of MAO catalyzed reactions and MAO inhibition in the context of neurodegeneration will be discussed.

Why p-OMe- and p-Cl-β-Methylphenethylamines Display Distinct Activities upon MAO-B Binding

Despite their structural and chemical commonalities, p-chloro-β-methylphenethylamine and p-methoxy-β-methylphenethylamine display distinct inhibitory and substrate activities upon MAO-B binding. Density Functional Theory (DFT) quantum chemical calculations reveal that β-methylation and para-substitution underpin the observed activities sustained by calculated transition state energy barriers, attained conformations and key differences in their interactions in the enzyme’s substrate binding site. Although both compounds meet substrate requirements, it is clear that β-methylation along with the physicochemical features of the para-substituents on the aromatic ring determine the activity of these compounds upon binding to the MAO B-isoform. While data for a larger set of compounds might lend generality to our conclusions, our experimental and theoretical results strongly suggest that the contrasting activities displayed depend on the conformations adopted by these compounds when they bind to the enzyme.

Insights into enzyme point mutation effect by molecular simulation: phenylethylamine oxidation catalyzed by monoamine oxidase A

The I335Y point mutation effect on the kinetics of phenylethylamine decomposition catalyzed by monoamine oxidase A was elucidated by means of molecular simulation. The established empirical valence bond methodology was used in conjunction with the free energy perturbation sampling technique and a classical force field representing the state of reactants and products. The methodology allows for the simulation of chemical reactions, in the present case the breaking of the α-C-H bond in a phenylethylamine substrate and the subsequent hydrogen transfer to the flavin cofactor, resulting in the formation of the N-H bond on flavin. The empirical parameters were calibrated against the experimental data for the simulated reaction in a wild type protein and then used for the calculation of the reaction free energy profile in the I335Y mutant. In very good agreement with the measured kinetic data, mutation increases the free energy barrier for the rate limiting step by slightly more than 1 kcal mol(-1) and consequently decreases the rate constant by about an order of magnitude. The magnitude of the computed effect slightly varies with simulation settings, but always remains in reasonable agreement with the experiment. Analysis of trajectories reveals a major change in the interaction between phenyl rings of the substrate and the neighboring Phe352 residue upon the I335Y mutation due to the increased local polarity, leading to an attenuated quadrupole interaction between the rings and destabilization of the transition state. Additionally, the increased local polarity in the mutant allows for a larger number of water molecules to be present near the active site, effectively shielding the catalytic effect of the enzyme and contributing to the increased barrier.