Empirical Valence Bond Research Articles

L-DOPA Autoxidation: An Empirical Valence Bond Simulation of the Reactive Step.

L-DOPA, or levodopa, plays an important role in the treatment of Parkinson's disease, a debilitating neurological disorder. It acts as a precursor to dopamine, a neurotransmitter crucial for the regulation of motor functions. Administered orally, L-DOPA easily crosses the blood-brain barrier and converts into dopamine in the brain, relieving symptoms such as tremors and rigidity. However, its prolonged use can lead to complications. A significant concern with L-DOPA is its conversion to dopaquinone, a quinone metabolite that enters the redox cycle and continuously produces hydrogen peroxide. In addition, L-DOPA, which resembles tyrosine with an additional hydroxyl group, can randomly incorporate into the proteins of dopaminergic neurons and thus become an additional source of oxidative stress in Parkinson's patients. In this study, we scrutinized the rate-limiting step of L-DOPA autoxidation in aqueous solution. The reaction we studied is an intramolecular Michael addition concerted with a proton transfer from the amino group. Using the Empirical Valence Bond (EVB) method, we computed the free energy profiles of the reaction in water. The calculated barrier of 30.93 ± 1.12 kcal/mol is in reasonable agreement with the experimental barrier of 27.55 kcal/mol. This agreement confirms the validity of the studied mechanism and demonstrates the applicability of our simulation methodology for studying the autoxidation kinetics of L-DOPA within proteins.

Toward transferable empirical valence bonds: Making classical force fields reactive.

The empirical valence bond technique allows classical force fields to model reactive processes. However, parametrization from experimental data or quantum mechanical calculations is required for each reaction present in the simulation. We show that the parameters present in the empirical valence bond method can be predicted using a neural network model and the SMILES strings describing a reaction. This removes the need for quantum calculations in the parametrization of the empirical valence bond technique. In doing so, we have taken the first steps toward defining a new procedure for enabling reactive atomistic simulations. This procedure would allow researchers to use existing classical force fields for reactive simulations, without performing additional quantum mechanical calculations.

Assessing the Catalytic Role of Native Glucagon Amyloid Fibrils.

Glucagon stands out as a pivotal peptide hormone, instrumental in controlling blood glucose levels and lipid metabolism. While the formation of glucagon amyloid fibrils has been documented, their biological functions remain enigmatic. Recently, we demonstrated experimentally that glucagon amyloid fibrils can act as catalysts in several biological reactions including esterolysis, lipid hydrolysis, and dephosphorylation. Herein, we present a multiscale quantum mechanics/molecular mechanics (QM/MM) simulation of the acylation step in the esterolysis of para-nitrophenyl acetate (p-NPA), catalyzed by native glucagon amyloid fibrils, serving as a model system to elucidate their catalytic function. This step entails a concerted mechanism, involving proton transfer from serine to histidine, followed by the nucleophilic attack of the serine oxy anion on the carbonyl carbon of p-NPA. We computed the binding energy and free-energy profiles of this reaction using the protein-dipole Langevin-dipole (PDLD) within the linear response approximation (LRA) framework (PDLD/S-LRA-2000) and the empirical valence bond (EVB) methods. This included simulations of the reaction in an aqueous environment and in the fibril, enabling us to estimate the catalytic effect of the fibril. Our EVB calculations obtained a barrier of 23.4 kcal mol-1 for the enzyme-catalyzed reaction compared to the experimental value of 21.9 kcal mol-1 (and a calculated catalytic effect of 3.2 kcal mol-1 compared to the observed effect of 4.7 kcal mol-1). This close agreement together with the barrier reduction when transitioning from the reference solution reaction to the amyloid fibril provides supporting evidence to the catalytic role of glucagon amyloid fibrils. Moreover, employing the PDLD/S-LRA-2000 approach further reinforced exclusively the enzyme's catalytic role. The results presented in this study contribute significantly to our understanding of the catalytic role of glucagon amyloid fibrils, marking, to the best of our knowledge, the first-principles mechanistic investigation of fibrils using QM/MM methods. Therefore, our findings offer fruitful insights for future research into the mechanisms of related amyloid catalysis.

Why Do Empirical Valence Bond Simulations Yield Accurate Arrhenius Plots?

Computer simulations of the temperature dependence of enzyme reactions using the empirical valence bond (EVB) method have proven to give very accurate results in terms of the thermodynamic activation parameters. Here, we analyze the reasons for why such simulations are able to correctly capture activation enthalpies and entropies and how sensitive these quantities are to parametrization of the reactive potential energy function. We examine first the solution reference reaction for the enzyme ketosteroid isomerase, which corresponds to the acetate catalyzed deprotonation of the steroid in water. The experimentally determined activation parameters for this reaction turn out to be remarkably well reproduced by the calculations. By modifying the EVB potential so that the activation and reaction free energies become significantly shifted, we show that the activation entropy is basically invariant to such changes and that ΔS⧧ is instead determined by the specific mixture of the underlying force fields in the transition state region. The coefficients of this mixture do not change appreciably when the EVB potential is modified within reasonable limits, and hence, the estimate of ΔS⧧ becomes very robust. This is further verified by examining a more complex concerted hydride and proton transfer reaction in the enzyme hydroxybutyrate dehydrogenase.

Accurate Computation of Thermodynamic Activation Parameters in the Chorismate Mutase Reaction from Empirical Valence Bond Simulations.

Chorismate mutase (CM) enzymes have long served as model systems for benchmarking new methods and tools in computational chemistry. Despite the enzymes' prominence in the literature, the extent of the roles that activation enthalpy and entropy play in catalyzing the conversion of chorismate to prephenate is still subject to debate. Knowledge of these parameters is a key piece in fully understanding the mechanism of chorismate mutases. Within this study, we utilize EVB/MD free energy perturbation calculations at a range of temperatures, allowing us to extract activation enthalpies and entropies from an Arrhenius plot of activation free energies of the reaction catalyzed by a monofunctional Bacillus subtilis CM and the promiscuous enzyme isochorismate pyruvate lyase of Pseudomonas aeruginosa. In comparison to the uncatalyzed reaction, our results show that both enzyme-catalyzed reactions exhibit a substantial reduction in activation enthalpy, while the effect on activation entropy is relatively minor, demonstrating that enzyme-catalyzed CM reactions are enthalpically driven. Furthermore, we observe that the monofunctional CM from B. subtilis more efficiently catalyzes this reaction than its promiscuous counterpart. This is supported by a structural analysis of the reaction pathway at the transition state, from which we identified key residues explaining the enthalpically driven nature of the reactions and also the difference in efficiencies between the two enzymes.

Elucidation of the α-Ketoamide Inhibition Mechanism: Revealing the Critical Role of the Electrostatic Reorganization Effect of Asp17 in the Active Site of the 20S Proteasome.

The 20S proteasome is an attractive drug target for the development of anticancer agents because it plays an important role in cellular protein degradation. It has a threonine residue that can act as a nucleophile to attack inhibitors with an electrophilic warhead, forming a covalent adduct. Fundamental understanding of the reaction mechanism between covalent inhibitors and the proteasome may assist the design and refinement of compounds with the desired activity. In this study, we investigated the covalent inhibition mechanism of an α-keto phenylamide inhibitor of the proteasome. We calculated the noncovalent binding free energy using the PDLD/S-LRA/β method and the reaction free energy through the empirical valence bond method (EVB). Several possible reaction pathways were explored. Subsequently, we validated the calculated activation and reaction free energies of the most plausible pathways by performing kinetic experiments. Furthermore, the effects of different ionization states of Asp17 on the activation energy at each step were also discussed. The results revealed that the ionization states of Asp17 remarkably affect the activation energies and there is an electrostatic reorganization of Asp17 during the course of the reaction. Our results demonstrate the critical electrostatic effect of Asp17 in the active site of the 20S proteasome.

Efficient Empirical Valence Bond Simulations with GROMACS.

We describe a protocol to perform empirical valence bond (EVB) simulations using GROMACS software. EVB is a fast and reliable method that allows one to determine the reaction free-energy profiles in complex systems, such as enzymes, by employing classical force fields to represent a chemical reaction. Therefore, running EVB simulations is basically as fast as any classical molecular dynamics simulation, and the method uses standard free-energy calculations to map the free-energy change along a given reaction path. To exemplify and validate our EVB implementation, we replicated two cases of our earlier enzyme simulations. One of these addresses the decomposition of the activation free energy into its enthalpic and entropic components, and the other is focused on calculating the overall catalytic effect of the enzyme compared to the same reaction in water. These two examples give virtually identical results to those obtained with programs that were specifically designed for EVB simulations and show that the GROMACS implementation is robust and can be used for very large systems.

Caracal: A Versatile Ring Polymer Molecular Dynamics Simulation Package.

A new open-source program package named Caracal covering simulations of molecular systems with ring polymer molecular dynamics (RPMD) is presented. It combines a powerful RPMD implementation including chemical reaction rate calculations and biased periodic and nonperiodic samplings with a collection of easy to set up potential energy surface (PES) methodologies, thus delivering an all-inclusive approach. Most implemented PESs are based on the QMDFF and EVB-QMDFF methods. Where the quantum mechanically derived force field (QMDFF) can be set up for an arbitrary molecular system in a black-box fashion, the empirical valence bond (EVB)-QMDFF connects two QMDFFs and is able to represent the PES of a chemical reaction. With our previously published flavors of this composite method, PESs for almost arbitrary gas phase thermal ground state reactions can be set up. Given an optimized reaction path, the mechanism of the reaction can be classified and RPMD rate constants can be obtained via umbrella sampling and recrossing calculations on an EVB-QMDFF PES. Further, QMDFFs can be polymerized for the description of liquid systems. In this paper, the internal structure as well as the handling philosophy of Caracal are outlined. Further, examples of the different possible kinds of calculations are given.

Examination of the performance of degraded Nafion membrane with graphene oxide

The degradation of Nafion membrane easily occurs after a long run of proton exchange membrane fuel cell. Graphene-oxide (GO) is regarded as an ideal material to form composite membrane. In this work, a multistate empirical valence bond (MS-EVB) method and molecular dynamics method are employed to investigate the effect of GO on transport and solvation characteristics of protons in the degraded hydration membrane. The hydrophilic cluster and its connectivity are analyzed. The interaction between GO, Nafion membrane and water molecules is discussed. The influence of GO on the thermal conduction of Nafion membrane is also evaluated. The results demonstrate that with the addition of GO, the water cluster connectivity of Nafion membrane becomes better. The distribution of hydronium ions and membranes in the water is more uniform, which increases the proton conductivity. The larger GO may not interact well with water molecules, but promotes the thermal conductivity of degraded Nafion membrane more obviously.

Acidic Conditions Impact Hydrophobe Transfer across the Oil-Water Interface in Unusual Ways.

Molecular dynamics simulation and enhanced free energy sampling are used to study hydrophobic solute transfer across the water-oil interface with explicit consideration of the effect of different electrolytes: hydronium cation (hydrated excess proton) and sodium cation, both with chloride counterions (i.e., dissociated acid and salt, HCl and NaCl). With the Multistate Empirical Valence Bond (MS-EVB) methodology, we find that, surprisingly, hydronium can to a certain degree stabilize the hydrophobic solute, neopentane, in the aqueous phase and including at the oil-water interface. At the same time, the sodium cation tends to "salt out" the hydrophobic solute in the expected fashion. When it comes to the solvation structure of the hydrophobic solute in the acidic conditions, hydronium shows an affinity to the hydrophobic solute, as suggested by the radial distribution functions (RDFs). Upon consideration of this interfacial effect, we find that the solvation structure of the hydrophobic solute varies at different distances from the oil-liquid interface due to a competition between the bulk oil phase and the hydrophobic solute phase. Together with an observed orientational preference of the hydroniums and the lifetime of water molecules in the first solvation shell of neopentane, we conclude that hydronium stabilizes to a certain degree the dispersal of neopentane in the aqueous phase and eliminates any salting out effect in the acid solution; i.e., the hydronium acts like a surfactant. The present molecular dynamics study provides new insight into the hydrophobic solute transfer across the water-oil interface process, including for acid and salt solutions.

Loop dynamics and evolvability in the TIM-barrel proteins of histidine and tryptophan biosynthesis.

Recent years have witnessed increasing interest in the role of conformational dynamics in enzyme evolution and engineering. In this context, TIM-barrel proteins are important model systems, as the TIM-barrel is ubiquitous, highly evolvable, supports a diverse range of chemistry, and the active site is typically decorated by one or more mobile loops that play an important role in catalysis and/or regulation of these enzymes. Here, the TIM-barrel proteins of histidine and tryptophan biosynthesis (HisA, TrpF and PriA) stand out, as their active sites are decorated by up to three long mobile loops, that move independently of each other, and the conformations of which are critically important for catalysis. These motions appear to be ligand-gated and occur on sufficiently slow timescales (ms or slower) to be at least partially rate-limiting. In recent work, we combined conventional and enhanced molecular dynamics simulations with empirical valence bond simulations to investigate the conformational behavior of the catalytic loops in these enzymes, and the link between their dynamical behavior and substrate selectivity and catalysis. We have explored the role of loop dynamics in facilitating the bifunctionality of PriA towards substrates ProFAR and PRA (precursors for histidine and tryptophan biosynthesis, respectively), as well as during the real-time evolution of HisA to yield HisA-specific, TrpF-specific and bifunctional variants. Finally, the insights this provides emphasizes the potential of loop engineering as a powerful tool to manipulate the selectivity and substrate scope of the TIM-barrel scaffold, which has historically been one of the most important scaffolds for de novo enzyme design.

Operando Condition Reaction Modeling Shows Highly Dynamic Attachment of Oligomeric Ruthenium Catalysts

To increase the stability and current density of molecular-catalyst-based electroanodes for water oxidation, immobilization of the catalysts at the electrode surface is a common strategy. A prominent example is the oligomerized Ru(tda) molecular catalyst, which showed outstanding current densities even at neutral pH values. One of the most challenging aspects of immobilized catalysts is to understand the interaction between the catalyst and the surface under operando conditions. Experiments are often performed under model conditions, and computational methods to study reaction steps are typically limited to a few hundred atoms. In this study, we combined three computational methods, density functional theory electronic structure computations, molecular dynamics for large-scale simulations of the catalyst–solid interaction, and empirical valence bond for reaction modeling the catalyst at the interface of a large carbon support and a phosphate water buffer. These techniques allowed us to explore the combined effects of solvent, hydrophobic directionality, and electric field on the attachment and reactivity of a Ru(tda) pentamer at a graphene surface. Our simulations have a perfect agreement with the experimental characterization under model conditions. However, we find that under operando conditions, where the catalyst is oxidized to the active RuV state, with a phosphate-containing electrolyte and an applied electric field, the attachment is completely reversed compared to the model conditions with RuII and organic solvents. This reversed attachment leads to a water-excluded region close to the active RuV═O center. The EVB reaction modeling showed that the reaction could still proceed to form an O–O bond via an oxide relay mechanism, where a dangling carboxylate reacts with the oxo via nucleophilic attack. We find that the activation energies are identical in water solution and at the electrode surface, showing how this mechanism is key to highly active molecular water oxidation catalysts immobilized on surfaces. Since attachment to surfaces could have a strong, and often negative, influence on the reactions, this study provides a guideline on how to model reactions without compromising the complexity of the electrode environment.

OntoPESScan: An Ontology for Potential Energy Surface Scans.

In this work, a new OntoPESScan ontology is developed for the semantic representation of one-dimensional potential energy surface (PES) scans, a central concept in computational chemistry. This ontology is developed in line with knowledge graph principles and The World Avatar (TWA) project. OntoPESScan is linked to other ontologies for chemistry in TWA, including OntoSpecies, which helps uniquely identify species along the PES and access their properties, and OntoCompChem, which allows the association of potential energy surfaces with quantum chemical calculations and the concepts used to derive them. A force-field fitting agent is also developed that makes use of the information in the OntoPESScan ontology to fit force fields to reactive surfaces of interest on the fly by making use of the empirical valence bond methodology. This agent is demonstrated to successfully parametrize two cases, namely, a PES scan on ethanol and a PES scan on a localized π-radical PAH hypothesized to play a role in soot formation during combustion. OntoPESScan is an extension to the capabilities of TWA and, in conjunction with potential further ontological support for molecular dynamics and reactions, will further progress toward an open, continuous, and self-growing knowledge graph for chemistry.

Q-RepEx: A Python pipeline to increase the sampling of empirical valence bond simulations

The exploration of chemical systems occurs on complex energy landscapes. Comprehensively sampling rugged energy landscapes with many local minima is a common problem for molecular dynamics simulations. These multiple local minima trap the dynamic system, preventing efficient sampling. This is a particular challenge for large biochemical systems with many degrees of freedom. Replica exchange molecular dynamics (REMD) is an approach that accelerates the exploration of the conformational space of a system, and thus can be used to enhance the sampling of complex biomolecular processes. In parallel, the empirical valence bond (EVB) approach is a powerful approach for modeling chemical reactivity in biomolecular systems. Here, we present an open-source Python-based tool that interfaces with the Q simulation package, and increases the sampling efficiency of the EVB free energy perturbation/umbrella sampling approach by means of REMD. This approach, Q-RepEx, both decreases the computational cost of the associated REMD-EVB simulations, and opens the door to more efficient studies of biochemical reactivity in systems with significant conformational fluctuations along the chemical reaction coordinate.

Analyzing the Reaction of Orotidine 5'-Phosphate Decarboxylase as a Way to Examine Some Key Catalytic Proposals.

This study analyzes the origin of enzyme catalysis by focusing on the reaction of orotidine 5'-phosphate decarboxylase (ODCase). This reaction involves an enormous catalytic effect of 23 kcal/mol that has been attributed to reactant state destabilization associated with the use of binding energy through the so-called Circe effect. However, our early studies and subsequent key experiments have shown that the presumed effect of the binding energy (namely, the strain exerted by a bond to a phosphate group) does not contribute to the catalysis. In this study, we perform quantitative empirical valence bond calculations that reproduce the catalytic effect of ODCase and the effect of removing the phosphate side chain. The calculations demonstrate that the effect of the phosphate is due to a change in reorganization energy and should not be described as an induced fit effect. Similarly, we show that the overall catalytic effect is due to electrostatic transition state stabilization, which again reflects the smaller reorganization energy in the enzyme than in water. We also elaborate on the problems with the induced fit proposal, including the fact that it does not serve to tell us what the actual origin of the action of the catalytic effect is. In addition to the above points, we use this paper to discuss misconceptions about the meaning of the preorganization effect, as well as other misunderstandings of what is being done in consistent calculations of enzyme catalysis.

Exploring the Role of Chemical Reactions in the Selectivity of Tyrosine Kinase Inhibitors.

A variety of diseases are associated with tyrosine kinase enzymes that activate many proteins via signal transduction cascades. The similar ATP-binding pockets of these tyrosine kinases make it extremely difficult to design selective covalent inhibitors. The present study explores the contribution of the chemical reaction steps to the selectivity of the commercialized inhibitor acalabrutinib over the Bruton's tyrosine kinase (BTK) and the interleukin-2-inducible T-cell kinase (ITK). Ab initio and empirical valence bond (EVB) simulations of the two kinases indicate that the most favorable reaction path involves a water-assisted mechanism of the 2-butynamide reactive group of acalabrutinib. BTK reacts with acalabrutinib with a substantially lower barrier than ITK, according to our calculated free-energy profile and kinetic simulations. Such a difference is due to the microenvironment of the active site, as further supported by a sequence-based analysis of specificity determinants for several commercialized inhibitors. Our study involves a new approach of simulating directly the IC50 and inactivation efficiency keff, instead of using the standard formulas. This new strategy is particularly important in studies of covalent inhibitors with a very exothermic bonding step. Overall, our results demonstrate the importance of understanding the chemical reaction steps in designing selective covalent inhibitors for tyrosine kinases.

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 , while together with Arg233 stabilize the inorganic phosphate (Pi); (ii) Asp168 and three water molecules coordinate , and 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 , 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.

Complex Loop Dynamics Underpin Activity, Specificity, and Evolvability in the (βα)8 Barrel Enzymes of Histidine and Tryptophan Biosynthesis.

Enzymes are conformationally dynamic, and their dynamical properties play an important role in regulating their specificity and evolvability. In this context, substantial attention has been paid to the role of ligand-gated conformational changes in enzyme catalysis; however, such studies have focused on tremendously proficient enzymes such as triosephosphate isomerase and orotidine 5′-monophosphate decarboxylase, where the rapid (μs timescale) motion of a single loop dominates the transition between catalytically inactive and active conformations. In contrast, the (βα)8-barrels of tryptophan and histidine biosynthesis, such as the specialist isomerase enzymes HisA and TrpF, and the bifunctional isomerase PriA, are decorated by multiple long loops that undergo conformational transitions on the ms (or slower) timescale. Studying the interdependent motions of multiple slow loops, and their role in catalysis, poses a significant computational challenge. This work combines conventional and enhanced molecular dynamics simulations with empirical valence bond simulations to provide rich details of the conformational behavior of the catalytic loops in HisA, PriA, and TrpF, and the role of their plasticity in facilitating bifunctionality in PriA and evolved HisA variants. In addition, we demonstrate that, similar to other enzymes activated by ligand-gated conformational changes, loops 3 and 4 of HisA and PriA act as gripper loops, facilitating the isomerization of the large bulky substrate ProFAR, albeit now on much slower timescales. This hints at convergent evolution on these different (βα)8-barrel scaffolds. Finally, our work reemphasizes the potential of engineering loop dynamics as a tool to artificially manipulate the catalytic repertoire of TIM-barrel proteins.

Essential Functional Interplay of the Catalytic Groups in Acid Phosphatase.

The cooperative interplay between the functional devices of a preorganized active site is fundamental to enzyme catalysis. An in-depth understanding of this phenomenon is central to elucidating the remarkable efficiency of natural enzymes and provides an essential benchmark for enzyme design and engineering. Here, we study the functional interconnectedness of the catalytic nucleophile (His18) in an acid phosphatase by analyzing the consequences of its replacement with aspartate. We present crystallographic, biochemical, and computational evidence for a conserved mechanistic pathway via a phospho-enzyme intermediate on Asp18. Linear free-energy relationships for phosphoryl transfer from phosphomonoester substrates to His18/Asp18 provide evidence for the cooperative interplay between the nucleophilic and general-acid catalytic groups in the wild-type enzyme, and its substantial loss in the H18D variant. As an isolated factor of phosphatase efficiency, the advantage of a histidine compared to an aspartate nucleophile is ∼104-fold. Cooperativity with the catalytic acid adds ≥102-fold to that advantage. Empirical valence bond simulations of phosphoryl transfer from glucose 1-phosphate to His and Asp in the enzyme explain the loss of activity of the Asp18 enzyme through a combination of impaired substrate positioning in the Michaelis complex, as well as a shift from early to late protonation of the leaving group in the H18D variant. The evidence presented furthermore suggests that the cooperative nature of catalysis distinguishes the enzymatic reaction from the corresponding reaction in solution and is enabled by the electrostatic preorganization of the active site. Our results reveal sophisticated discrimination in multifunctional catalysis of a highly proficient phosphatase active site.

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.