The directive of the protein: how does cytochrome P450 select the mechanism of dopamine formation?

Origins of Catalysis in Non-Heme Fe(II)/2-Oxoglutarate-Dependent Histone Lysine Demethylase KDM4A with Differently Methylated Histone H3 Peptides

Histone lysine demethylase 4 A (KDM4A), a non-heme Fe(II)/2-oxoglutarate (2OG) dependent oxygenase that catalyzes the demethylation of tri-methylated lysine residues at the 9, 27, and 36 positions of histone H3 (H3 K9me3, H3 K27me3, and H3 K36me3). These methylated residues show contrasting transcriptional roles; therefore, understanding KDM4A's catalytic mechanisms with these substrates is essential to explain the factors that control the different sequence-dependent demethylations. In this study, we use molecular dynamics (MD)-based combined quantum mechanics/molecular mechanics (QM/MM) methods to investigate determinants of KDM4A catalysis with H3 K9me3, H3 K27me3 and H3 K36me3 substrates. In KDM4A-H3(5-14)K9me3 and KDM4A-H3(23-32)K27me3 ferryl complexes, the O-H distance positively correlates with the activation barrier of the rate-limiting step, however in the KDM4A-H3(32-41)K36me3, no direct one-to-one relationship was found implying that the synergistic effects between the geometric parameters, second sphere interactions and the intrinsic electric field contribute for the effective catalysis for this substrate. The intrinsic electric field along the Fe-O bond changes between the three complexes and shows a positive correlation with the HAT activation barrier, suggesting that modulating electric field can be used for fine engineering KDM catalysis with a specific substrate. The results reveal how KDM4A uses a combination of strategies to enable near equally efficient demethylation of different H3Kme3 residues.

Structural diversity, evolutionary origin, and metabolic engineering of plant specialized benzylisoquinoline alkaloids

Covering: up to June 2024Benzylisoquinoline alkaloids (BIAs) represent a diverse class of plant specialized metabolites derived from L-tyrosine, exhibiting significant pharmacological properties such as anti-microbial, anti-spasmodic, anti-cancer, cardiovascular protection, and analgesic effects. The industrial production of valuable BIAs relies on extraction from plants; however, challenges concerning their low concentration and efficiency hinder drug development. Hence, alternative approaches, including biosynthesis and chemoenzymatic synthesis, have been explored. Model species like Papaver somniferum and Coptis japonica have played a key role in unraveling the biosynthetic pathways of BIAs; however, many aspects, particularly modified steps like oxidation and methylation, remain unclear. Critical enzymes, e.g., CYP450s and methyltransferases, play a substantial role in BIA backbone formation and modification, which is essential for understanding the origin and adaptive evolution of these plant specialized metabolites. This review comprehensively analyzes the structural diversity of reported BIAs and their distribution in plant lineages. In addition, the progress in understanding biosynthesis, evolution, and catalytic mechanisms underlying BIA biosynthesis is summarized. Finally, we discuss the progress and challenges in metabolic engineering, providing valuable insights into BIA drug development and the sustainable utilization of BIA-producing plants.

QM/MM Study of the Reaction Mechanism of L-Tyrosine Hydroxylation Catalyzed by the Enzyme CYP76AD1

We have studied the hydroxylation mechanism of l-Tyr by the heme-dependent enzyme CYP76AD1 from the sugar beet (Beta vulgaris). This enzyme has a promising biotechnological application in modified yeast strains to produce medicinal alkaloids, an alternative to the traditional opium poppy harvest. A generative machine learning software based on AlphaFold was used to build the structure of CYP76AD1 since there are no structural data for this specific enzyme. After model validation, l-Tyr was docked in the active site of CYP76AD1 to assemble the reactive complex, whose catalytic distances remained stable throughout the 100 ns of MD simulation. Subsequent QM/MM calculations elucidated that l-Tyr hydroxylation occurs in two steps: hydrogen abstraction from l-Tyr by CpdI, forming an l-Tyr radical, and subsequent radical rebound, corresponding to a rate-limiting step of 16.0 kcal·mol-1. Our calculations suggest that the hydrogen abstraction step should occur in the doublet state, while the radical rebound should happen in the quartet state. The clarification of the reaction mechanism of CYP76AD1 provides insights into the rational optimization of the biosynthesis of alkaloids to eliminate the use of opium poppy.

Unusual catalytic strategy by non-heme Fe(ii)/2-oxoglutarate-dependent aspartyl hydroxylase AspH

Biocatalytic C-H oxidation reactions are of important synthetic utility, provide a sustainable route for selective synthesis of important organic molecules, and are an integral part of fundamental cell processes. The multidomain non-heme Fe(ii)/2-oxoglutarate (2OG) dependent oxygenase AspH catalyzes stereoselective (3R)-hydroxylation of aspartyl- and asparaginyl-residues. Unusually, compared to other 2OG hydroxylases, crystallography has shown that AspH lacks the carboxylate residue of the characteristic two-His-one-Asp/Glu Fe-binding triad. Instead, AspH has a water molecule that coordinates Fe(ii) in the coordination position usually occupied by the Asp/Glu carboxylate. Molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) studies reveal that the iron coordinating water is stabilized by hydrogen bonding with a second coordination sphere (SCS) carboxylate residue Asp721, an arrangement that helps maintain the six coordinated Fe(ii) distorted octahedral coordination geometry and enable catalysis. AspH catalysis follows a dioxygen activation-hydrogen atom transfer (HAT)-rebound hydroxylation mechanism, unusually exhibiting higher activation energy for rebound hydroxylation than for HAT, indicating that the rebound step may be rate-limiting. The HAT step, along with substrate positioning modulated by the non-covalent interactions with SCS residues (Arg688, Arg686, Lys666, Asp721, and Gln664), are essential in determining stereoselectivity, which likely proceeds with retention of configuration. The tetratricopeptide repeat (TPR) domain of AspH influences substrate binding and manifests dynamic motions during catalysis, an observation of interest with respect to other 2OG oxygenases with TPR domains. The results provide unique insights into how non-heme Fe(ii) oxygenases can effectively catalyze stereoselective hydroxylation using only two enzyme-derived Fe-ligating residues, potentially guiding enzyme engineering for stereoselective biocatalysis, thus advancing the development of non-heme Fe(ii) based biomimetic C-H oxidation catalysts, and supporting the proposal that the 2OG oxygenase superfamily may be larger than once perceived.

Free Energy Profile Decomposition Analysis for QM/MM Simulations of Enzymatic Reactions

In enzyme mechanistic studies and mutant design, it is highly desirable to know the individual residue contributions to the reaction free energy and barrier. In this work, we show that such free energy contributions from each residue can be readily obtained by postprocessing ab initio quantum mechanical molecular mechanical (ai-QM/MM) free energy simulation trajectories. Specifically, through a mean force integration along the minimum free energy pathway, one can obtain the electrostatic, polarization, and van der Waals contributions from each residue to the free energy barrier. Separately, a similar analysis procedure allows us to assess the contribution from different collective variables along the reaction coordinate. The chorismate mutase reaction is used to demonstrate the utilization of these two trajectory analysis tools.

QM/MM Modeling Aided Enzyme Engineering in Natural Products Biosynthesis

Natural products and their derivatives are widely used across various industries, particularly pharmaceuticals. Modern engineered biosynthesis provides an alternative way of producing and meeting the growing need for diverse natural products. Natural enzymes, on the other hand, often exhibit unsatisfactory catalytic characteristics and necessitate further enzyme engineering modifications. QM/MM, as a powerful and extensively used computational tool in the field of enzyme catalysis, has been increasingly applied in rational enzyme engineering over the past decade. In this review, we summarize recent advances in QM/MM computational investigation on enzyme catalysis and enzyme engineering for natural product biosynthesis. The challenges and perspectives for future QM/MM applications aided enzyme engineering in natural product biosynthesis will also be discussed.

Applications and Potential of In Silico Approaches for Psychedelic Chemistry

Molecular-level investigations of the Central Nervous System have been revolutionized by the development of computational methods, computing power, and capacity advances. These techniques have enabled researchers to analyze large amounts of data from various sources, including genomics, in vivo, and in vitro drug tests. In this review, we explore how computational methods and informatics have contributed to our understanding of mental health disorders and the development of novel drugs for neurological diseases, with a special focus on the emerging field of psychedelics. In addition, the use of state-of-the-art computational methods to predict the potential of drug compounds and bioinformatic tools to integrate disparate data sources to create predictive models is also discussed. Furthermore, the challenges associated with these methods, such as the need for large datasets and the diversity of in vitro data, are explored. Overall, this review highlights the immense potential of computational methods and informatics in Central Nervous System research and underscores the need for continued development and refinement of these techniques and more inclusion of Quantitative Structure-Activity Relationships (QSARs).

Unraveling the Major Differences between the Trinuclear Cyclopentadienylmetal Carbonyl Chemistry of Cobalt and That of Nickel-A Theoretical Study

The geometries and energetics of the trinuclear cyclopentadienylmetal carbonyls Cp₃M₃(CO)_n (Cp = η⁵-C₅H₅); M = Co, Ni; n = 3, 2, 1, 0) have been investigated by density functional theory. The cobalt and nickel systems are found to be rather different owing to the different electronic configurations of the metal atoms. For cobalt, the small calculated energy separation of 5.0 kcal/mol between the two lowest-energy singlet Cp₃Co₃(μ₃-CO)(μ-CO)₂ and Cp₃Co₃(μ-CO)₃ tricarbonyl structures accounts for the experimental results of both isomers as stable species that can be isolated and structurally characterized by X-ray crystallography. The corresponding Cp₃Ni₃(CO)₃ species in the nickel system are predicted not to be viable owing to exothermic CO dissociation to give the experimentally observed very stable Cp₃Ni₃(μ-CO)₂, which is found to be the lowest-energy isomer by a substantial margin of ∼25 kcal/mol. In all of the low-energy Cp₃M₃(CO)_n (n = 2, 1) structures, including that of the experimentally known triplet spin state Cp₃Co₃(μ₃-CO)₂, all of the carbonyl groups are face-bridging or face-semi-bridging μ₃-CO groups bonded to all three metal atoms of the M₃ triangle. In the lowest-energy carbonyl-free Cp₃M₃ (M = Co, Ni) structures, agostic C-H-M interactions are found using hydrogens of the Cp rings. In addition, the lowest-energy Cp₃Ni₃ is the only structure among all of the low-energy Cp₃M₃(CO)_n (M = Co, Ni; n = 3, 2, 1, 0) structures in which each Cp ring is a bridging rather than terminal ligand.

Evaluating the Transition State Stabilization/Destabilization Effects of the Electric Fields from Scaffold Residues by a QM/MM Approach

The protein scaffolds of enzymes not only provide structural support for the catalytic center but also exert preorganized electric fields for electrostatic catalysis. In recent years, uniform oriented external electric fields (OEEFs) have been widely applied to enzymatic reactions to mimic the electrostatic effects of the environment. However, the electric fields exerted by individual residues in proteins may be quite heterogeneous across the active site, with varying directions and strengths at different positions of the active site. Here, we propose a QM/MM-based approach to evaluate the effects of the electric fields exerted by individual residues in the protein scaffold. In particular, the heterogeneity of the residue electric fields and the effect of the native protein environment can be properly accounted for by this QM/MM approach. A case study of the O-O heterolysis reaction in the catalytic cycle of TyrH shows that (1) for scaffold residues that are relatively far from the active site, the heterogeneity of the residue electric field in the active site is not very significant and the electrostatic stabilization/destabilization due to each residue can be well approximated with the interaction energy between a uniform electric field and the QM region dipole; (2) for scaffold residues near the active site, the residue electric fields can be highly heterogeneous along the breaking O-O bond. In such a case, approximating the residue electric fields as uniform fields may misrepresent the overall electrostatic effect of the residue. The present QM/MM approach can be applied to evaluate the residues' electrostatic impact on enzymatic reactions, which also can be useful in computational optimization of electric fields to boost the enzyme catalysis.

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.

Local Electric Fields: From Enzyme Catalysis to Synthetic Catalyst Design

This Mini-Review Article outlines recent advances in the study of local electric field (LEF) governed enzyme catalysis and the application of the LEF principle in synthetic catalyst design. We start by discussing the electrostatics principles that drive enzyme catalysis, and its experimental verifications through vibrational Stark spectroscopy. Subsequently, we describe aspects of LEFs other than catalysis, i.e., induction of mechanistic crossovers, among others. Here, we focus on the early work done using computational tools, along with some recent contributions. Following an in-depth discussion of the role of LEFs in enzyme catalysis, we then highlight some recent works on designed local electric fields (D-LEF) and their applications in organic synthesis. Subsequently, we turn to D-LEFs in synthetic enzymes and supramolecular systems (cf. the work by the Head-Gordon group). We end by discussing some of the software packages that have been developed to analyze local electric fields computationally. Overall, the present Mini-Review Article paints an insightful picture of the current state of the art using LEF in enzyme catalysis and its application for further bioengineering and synthetic organic frameworks in a broad perspective.

Oriented internal electrostatic fields: an emerging design element in coordination chemistry and catalysis

The power of oriented electrostatic fields (ESFs) to influence chemical bonding and reactivity is a phenomenon of rapidly growing interest. The presence of strong ESFs has recently been implicated as one of the most significant contributors to the activity of select enzymes, wherein alignment of a substrate's changing dipole moment with a strong, local electrostatic field has been shown to be responsible for the majority of the enzymatic rate enhancement. Outside of enzymology, researchers have studied the impacts of "internal" electrostatic fields via the addition of ionic salts to reactions and the incorporation of charged functional groups into organic molecules (both experimentally and computationally), and "externally" via the implementation of bulk fields between electrode plates. Incorporation of charged moieties into homogeneous inorganic complexes to generate internal ESFs represents an area of high potential for novel catalyst design. This field has only begun to materialize within the past 10 years but could be an area of significant impact moving forward, since it provides a means for tuning the properties of molecular complexes via a method that is orthogonal to traditional strategies, thereby providing possibilities for improved catalytic conditions and novel reactivity. In this perspective, we highlight recent developments in this area and offer insights, obtained from our own research, on the challenges and future directions of this emerging field of research.

Electrostatic vs. inductive effects in phosphine ligand donor properties and reactivity

Enhanced rates and selectivity in enzymes are enabled in part by precisely tuned electric fields within active sites. Analogously, the use of charged groups to leverage electrostatics in molecular systems is a promising strategy to tune reactivity. However, separation of the through space and through bond effects of charged functional groups is a long standing challenge that limits the rational application of electric fields in molecular systems. To address this challenge we developed a method using the phosphorus selenium coupling value (J _P-Se) of anionic phosphine selenides to quantify the electrostatic contribution of the borate moiety to donor strength. In this analysis we report the synthesis of a novel anionic phosphine, PPh₂CH₂BF₃K, the corresponding tetraphenyl phosphonium and tetraethyl ammonium selenides [PPh₄][SePPh₂CH₂BF₃] and [TEA][SePPh₂CH₂BF₃], and the Rh carbonyl complex [PPh₄][Rh(acac)(CO)(PPh₂(CH₂BF₃))]. Solvent-dependent changes in J _P-Se were fit using Coulomb's law and support up to an 80% electrostatic contribution to the increase in donor strength of [PPh₄][SePPh₂CH₂BF₃] relative to SePPh₂Et, while controls with [TEA][SePPh₂CH₂BF₃] exclude convoluting ion pairing effects. Calculations using explicit solvation or point charges effectively replicate the experimental data. This J _P-Se method was extended to [PPh₄][SePPh₂(2-BF₃Ph)] and likewise estimates up to a 70% electrostatic contribution to the increase in donor strength relative to SePPh₃. The use of PPh₂CH₂BF₃K also accelerates C-F oxidative addition reactivity with Ni(COD)₂ by an order of magnitude in comparison to the comparatively donating neutral phosphines PEt₃ and PCy₃. This enhanced reactivity prompted the investigation of catalytic fluoroarene C-F borylation, with improved yields observed for less fluorinated arenes. These results demonstrate that covalently bound charged functionalities can exert a significant electrostatic influence under common solution phase reaction conditions and experimentally validate theoretical predictions regarding electrostatic effects in reactivity.

Local Electric Fields Dictate Function: The Different Product Selectivities Observed for Fatty Acid Oxidation by Two Deceptively Very Similar P450-Peroxygenases OleT and BSβ

Cytochrome P450 peroxygenases use hydrogen peroxide to hydroxylate long-chain fatty acids by bypassing the use of O₂ and a redox partner. Among the peroxygenases, P450_OleT uniquely performs decarboxylation of fatty acids and production of terminal olefins. This route taken by P450_OleT is intriguing, and its importance is augmented by the practical importance of olefin production. As such, this mechanistic choice merits elucidation. To address this puzzle, we use hybrid QM/MM calculations and MD simulations for the OleT enzyme as well as for the structurally analogous enzyme, P450_BSβ. The study of P450_OleT reveals that the protonated His85 in the wild-type P450_OleT plays a crucial role in steering decarboxylation activity by stabilizing the corresponding hydroxoiron(IV) intermediate (Cpd II). In contrast, for P450_BSβ in which Q85 replaces H85, the respective Cpd II species is unstable and it reacts readily with the substrate radical by rebound, producing hydroxylation products. As shown, this single-site difference creates in P450_OleT a local electric field (LEF), which is significantly higher than that in P450_BSβ. In turn, these LEF differences are responsible for the different stabilities of the respective Cpd II/radical intermediates and hence for different functions of the two enzymes. P450_BSβ uses the common rebound mechanism and leads to hydroxylation, whereas P450_OleT proceeds via decarboxylation and generates terminal olefins. Olefin production projects the power of a single residue to alter the LEF and the enzyme's function.

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.

Effect of External Electric Field on Tetrel Bonding Interactions in (FTF3···FH) Complexes (T = C, Si, Ge, and Sn)

A quantum chemical study was accomplished on the σ-hole interactions of the barely explored group IV elements, for the first time, in the absence and presence of the positively and negatively directed external electric field (EEF). The analyses of molecular electrostatic potential addressed the occurrence of the σ-hole on all the inspected tetrel atoms, confirming their salient versatility to engage in σ-hole interactions. MP2 energetic findings disclosed the occurrence of favorable σ-hole interactions within the tetrel bonding complexes. The tetrel bonding interactions became stronger in the order of C < Si < Ge < Sn for F-T-F3···FH complexes with the largest interaction energy amounting to -19.43 kcal/mol for the optimized F-Sn-F3···FH complex under the influence of +0.020 au EEF. The interaction energy conspicuously evolved by boosting the magnitude of the positively directed EEF value and declining the negatively directed EEF one. The decomposition analysis for the interaction energies was also executed in terms of symmetry-adapted perturbation theory, illuminating the dominant electrostatic contribution to all the studied complexes' interactions except carbon-based interactions controlled by dispersion forces. The outcomes that emerged from the current work reported significantly how the direction and strength of the EEF affect the tetrel-bonding interactions, leading to further improvements in the forthcoming studies of supramolecular chemistry and materials science.

Tris(Butadiene) Compounds versus Butadiene Oligomerization in Second-Row Transition Metal Chemistry: Effects of Increased Ligand Fields

The geometries, energetics, and preferred spin states of the second-row transition metal tris(butadiene) complexes (C₄H₆)₃M (M = Zr–Pd) and their isomers, including the experimentally known very stable molybdenum derivative (C₄H₆)₃Mo, have been examined by density functional theory. Such low-energy structures are found to have low-spin singlet and doublet spin states in contrast to the corresponding derivatives of the first-row transition metals. The three butadiene ligands in the lowest-energy (C₄H₆)₃M structures of the late second-row transition metals couple to form a C₁₂H₁₈ ligand that binds to the central metal atom as a hexahapto ligand for M = Pd but as an octahapto ligand for M = Rh and Ru. However, the lowest-energy (C₄H₆)₃M structures of the early transition metals have three separate tetrahapto butadiene ligands for M = Zr, Nb, and Mo or two tetrahapto butadiene ligands and one dihapto butadiene ligand for M = Tc. The low energy of the experimentally known singlet (C₄H₆)₃Mo structure contrasts with the very high energy of its experimentally unknown singlet chromium (C₄H₆)₃Cr analog relative to quintet (C₁₂H₁₈)Cr isomers with an open-chain C₁₂H₁₈ ligand.

External electric field effects on the σ-hole and lone-pair hole interactions of group V elements: a comparative investigation

σ-hole and lone-pair (lp) hole interactions of trivalent pnicogen-bearing (ZF3) compounds were comparatively scrutinized, for the first time, under field-free and external electric field (EEF) conditions. Conspicuously, the sizes of the σ-hole and lp-hole were increased by applying an EEF along the positive direction, while the sizes of both holes decreased through the reverse EEF direction. The MP2 energetic calculations of ZF3⋯FH/NCH complexes revealed that σ-holes exhibited more impressive interaction energies compared to the lp-holes. Remarkably, the strengths of σ-hole and lp-hole interactions evolved with the increment of the positive value of the considered EEF; i.e., the interaction energy increased as the utilized EEF value increased. Unexpectedly, under field-free conditions, nitrogen-bearing complexes showed superior strength for their lp-hole interactions than phosphorus-bearing complexes. However, the reverse picture was exhibited for the interaction energies of nitrogen- and phosphorus-bearing complexes interacting within lp-holes by applying the high values of a positively directed EEF. These results significantly demonstrate the crucial influence of EEF on the strength of σ-hole and lp-hole interactions, which in turn leads to an omnipresent enhancement for variable fields, including biological simulations and material science.

The molecular mechanism of P450-catalyzed amination of the pyrrolidine derivative of lidocaine: insights from multiscale simulations

Nitrogen heterocycles are key and prevalent motifs in drugs. Evolved variants of cytochrome P450_BM3 (CYP102A1) from Bacillus megaterium employ high-valent oxo-iron(iv) species to catalyze the synthesis of imidazolidine-4-ones via an intramolecular C–H amination. Herein, we use multi-scale simulations, including classical molecular dynamics (MD) simulations, quantum mechanical/molecular mechanical (QM/MM) calculations and QM calculations, to reveal the molecular mechanism of the intramolecular C–H amination of the pyrrolidine derivative of lidocaine bearing cyclic amino moieties catalyzed by the variant RP/FV/EV of P450_BM3, which bears five mutations compared to wild type. Our calculations show that overall catalysis includes both the enzymatic transformation in P450 and non-enzymatic transformation in water solution. The enzymatic transformation involves the exclusive hydroxylation of the C–H bond of the pyrrolidine derivative of lidocaine, leading to the hydroxylated intermediate, during which the substrate radical would be bypassed. The following dehydration and C–N coupling reactions are found to be much favored in aqueous situation compared to that in the non-polar protein environment. The present findings expand our understanding of the P450-catalyzed C(sp₃)–H amination reaction.

Nitrogen heterocycles are key and prevalent motifs in drugs.

Deciphering Reaction Determinants of Altered-Activity CYP2D6 Variants by Well-Tempered Metadynamics Simulation and QM/MM Calculations

The xenobiotic metabolizing enzyme CYP2D6 is the P450 cytochrome family member with the highest rate of polymorphism. This causes changes in the enzyme activity and specificity, which can ultimately lead to adverse reactions during drug treatment. To avoid or lower CYP-related toxicity risks, prediction of the most likely positions within a molecule where a metabolic reaction might occur is paramount. In order to obtain accurate predictions, it is crucial to understand all phenomena within the active site of the enzyme that contribute to an efficient substrate recognition and the subsequent catalytic reaction together with their relative weight within the overall thermodynamic context. This study aims to define the weight of the driving forces upon the C-H bond activation within CYP2D6 wild-type and a clinically relevant allelic variant with increased activity (CYP2D6*53) featuring two amino acid mutations in close vicinity of the heme. First, we investigated the steric and electrostatic complementarity of the substrate bufuralol using well-tempered metadynamics simulations with the aim to obtain the free energy profiles for each site of metabolism (SoM) within the different active sites. Second, the stereoelectronic complementarity was determined for each SoM within the two different active-site environments. Relying on the well-tempered metadynamics simulation energy profiles of each SoM, we identified the binding mode that was closest to the preferred transition-state geometry for efficient C-H bond activation. The binding modes were then used as starting structures for the quantum mechanics/molecular mechanics calculations performed to quantify the corresponding activation barriers. Our results show the relevance of the steric component in orienting the SoM in an energetically accessible position toward the heme. However, the corresponding intrinsic reactivity and electronic complementarity within the active site must be accurately evaluated in order to obtain a meaningful reaction prediction, from which the predominant SoM can be determined. The F120I mutation lowered the activation barrier for the major site and one of the minor SoMs. However, it had an impact neither on the CYP2D6 enantioselectivity preference of the oxidation reaction nor on the stereoselectivity from the substrate point of view.

Computational Studies on the Mechanism and Origin of the Different Regioselectivities of Manganese Porphyrin-Catalyzed C-H Bond Hydroxylation and Amidation of Equilenin Acetate

The manganese porphyrin-catalyzed C-H bond hydroxylation and amidation of equilenin acetate developed by Breslow and his co-worker have been investigated with density functional theory (DFT) calculations. The hydroxylation of C(sp²)-H bond of equilenin acetate leading to the 6-hydroxylated product is more favorable than the hydroxylation of C(sp³)-H bond of equilenin acetate, leading to the 11β-hydroxylation product. The computational results suggest that the C(sp²)-H bond hydroxylation of equilenin acetate undergoes an oxygen-atom-transfer mechanism, which is more favorable than the C(sp³)-H bond hydroxylation undergoing the hydrogen-atom-abstraction/oxygen-rebound (HAA/OR) mechanism by 1.6 kcal/mol. That is why, the 6-hydroxylated product is the major product and the 11β-hydroxylated product is the minor product. In contrast, the 11β-amidated product is the only observed product in manganese porphyrin-catalyzed amidation reaction. The benzylic amidation undergoes a hydrogen-atom-abstraction/nitrogen-rebound (HAA/NR) mechanism, in which hydrogen atom abstraction is followed by nitrogen rebound, leading to the 11β-amidated product. The benzylic C(sp³)-H bond amidation at the C-11 position is more favorable than aromatic amidation at the C-6 position by 4.9 kcal/mol. Therefore, the DFT computational results are consistent with the experiments that manganese porphyrin-catalyzed C-H bond hydroxylation and amidation of equilenin acetate have different regioselectivities.

Engineering Escherichia coli for High-Yielding Hydroxytyrosol Synthesis from Biobased l-Tyrosine

Hydroxytyrosol (HT) is a natural antioxidant with many associated health benefits. In this study, we established efficient enzymatic cascades for the synthesis of HT from biobased l-tyrosine. First, a dopamine-mediated route for HT production was investigated. The combination of native hydroxylase (HpaBC) from Escherichia coli and l-DOPA decarboxylase (DODC) from Pseudomonas putida could efficiently convert 5 mM l-tyrosine into dopamine with conversion above 90%. However, further incorporation of monoamine oxidase (MAO) from Micrococcus luteus and phenylacetaldehyde reductase (PAR) from Solanum lycopersicum only resulted in 3.47 mM HT with 69.4% conversion. Therefore, a second enzyme cascade that comprises HpaBC from E. coli, l-amino acid deaminase (LAAD) from Proteus mirabilis, α-keto acid decarboxylase (ARO10) from Saccharomyces cerevisiae, and PAR from S. lycopersicum was designed. This enzymatic route showed higher catalytic activity and efficiently synthesized HT. The 24.27 mM HT was obtained from 25 mM l-tyrosine with a high conversion of 97.1%, and 32.35 mM HT was produced using 50 mM l-tyrosine, which represents the highest HT titer using l-tyrosine as a substrate reported to date. In summary, we have developed a green and sustainable platform for efficient HT enzymatic synthesis.

Solvent Organization and Rate Regulation of a Menshutkin Reaction by Oriented External Electric Fields are Revealed by Combined MD and QM/MM Calculations

When and how do external electric fields (EEFs) lead to catalysis in the presence of a (polar or nonpolar) solvent? This is the question that is addressed here using a combination of molecular dynamics (MD) simulations, quantum mechanical/molecular mechanical calculations with EEF, and quantum mechanical/(local) electric field calculations. The paper focuses on a model reaction, the Menshutkin reaction between CH3I and pyridine in three solvents of varying polarity. Using MD simulations, we find that the EEF causes the solvent to undergo organization; the solvent molecules gradually align with the applied field as the field strength increases. The collective orientation of the solvent molecules modifies the electrostatic environment around the Menshutkin species and induces a global electric field pointing in the opposite direction of the applied EEF. The combination of these two entangled effects leads to partial or complete screening of the EEF, with the extent of screening being proportional to the polarity/polarizability of the solvent. Nevertheless, we find that catalysis of the Menshutkin reaction inevitably emerges once the EEF exceeds the opposing field of the organizing solvent, i.e., once polarization of the Menshutkin complex is observed to set in. Overall, our analysis provides a lucid and pictorial interpretation of the behavior of solutions in the presence of EEFs and indicates that EEF-mediated catalysis should, in principle, be feasible in bulk setups, especially for nonpolar and mildly polar solvents. By application of the charge-transfer paradigm, it is shown that the emergence of OEEF catalysis in solution can be generalized to other reactions as well.

Perfluoroolefin complexes versus perfluorometallacycles and perfluorocarbene complexes in cyclopentadienylcobalt chemistry

Fluorocarbons have been shown experimentally by Baker and coworkers to combine with the cyclopentadienylcobalt (CpCo) moiety to form fluoroolefin and fluorocarbene complexes as well as fluorinated cobaltacyclic rings. In this connection density functional theory (DFT) studies on the cyclopentadienylcobalt fluorocarbon complexes CpCo(L)(C_nF_2n) (L = CO, PMe₃; n = 3 and 4) indicate structures with perfluoroolefin ligands to be the lowest energy structures followed by perfluorometallacycle structures and finally by structures with perfluorocarbene ligands. Thus, for the CpCo(L)(C₃F₆) (L = CO, PMe₃) complexes, the perfluoropropene structure has the lowest energy, followed by the perfluorocobaltacyclobutane structure and the perfluoroisopropylidene structure less stable by 8 to 11 kcal mol^-1, and the highest energy perfluoropropylidene structure less stable by more than 12 kcal mol^-1. For the two metal carbene structures Cp(L)Co[double bond, length as m-dash]C(CF₃)₂ and Cp(L)Co[double bond, length as m-dash]CF(C₂F₅), the former is more stable than the latter, even though the latter has Fischer carbene character. For the CpCo(L)(C₄F₈) (L = CO, PMe₃) complexes, the perfluoroolefin complex structures have the lowest energies, followed by the perfluorometallacycle structures at 10 to 20 kcal mol^-1, and the structures with perfluorocarbene ligands at yet higher energies more than 20 kcal mol^-1 above the lowest energy structure. This is consistent with the experimentally observed isomerization of the perfluorinated cobaltacyclobutane complexes CpCo(PPh₂Me)(-CFR-CF₂-CF₂-) (R = F, CF₃) to the perfluoroolefin complexes CpCo(PPh₂Me)(RCF[double bond, length as m-dash]CF₂) in the presence of catalytic quantities of HN(SO₂CF₃)₂. Further refinement of the relative energies by the state-of-the-art DLPNO-CCSD(T) method gives results essentially consistent with the DFT results summarized above.

Computational Biotransformation Profile of Emerging Phenolic Pollutants by Cytochromes P450: Phenol-Coupling Mechanism

Phenols are ubiquitous environmental pollutants, whose biotransformation involving phenol coupling catalyzed by cytochromes P450 may produce more lipophilic and toxic metabolites. Density functional theory (DFT) computations were performed to explore the debated phenol-coupling mechanisms, taking triclosan as a model substrate. We find that a diradical pathway facilitated by compound I and protonated compound II of P450 is favored vs alternative radical addition or electron-transfer mechanisms. The identified diradical coupling resembles a "two-state reactivity" from compound I characterized by significantly high rebound barriers of the phenoxy radicals, which can be formulated into three equations for calculating the ratio [coupling]/[hydroxylation]. A higher barrier for rebound than for H-abstraction in high-spin triclosan can facilitate the phenoxy radical dissociation and thus enable phenol coupling, while H-abstraction/radical rebound causing phenol hydroxylation via minor rebound barriers mostly occurs via the low-spin state. Therefore, oxidation of triclosan by P450 fits the first equation with a ratio [coupling]/[hydroxylation] of 1:4, consistent with experimental data indicating different extents of triclosan coupling (6-40%). The high rebound barrier of phenoxy radicals, as a key for the mechanistic identification of phenol coupling vs hydroxylation, originates from their weak electron donor ability due to spin aromatic delocalization. We envision that the revealed mechanism can be extended to the cross-coupling reactions between different phenolic pollutants, and the coupling reactions of several other aromatic pollutants, to infer unknown metabolites.

Molecular Docking as a Promising Predictive Model for Silver Nanoparticle-Mediated Inhibition of Cytochrome P450 Enzymes

Cytochrome P450 (CYP) enzymes are responsible for oxidative metabolisms of a large number of xenobiotics. In this study, we investigated interactions of silver nanoparticles (AgNPs) and silver ions (Ag⁺) with six CYP isoforms, namely, CYP1A2, CYP2C9, CYP2C19, CYP2D6, CYP2E1, and CYP3A4, within CYP-specific inhibitor-binding pockets by molecular docking and quantum mechanical (QM) calculations. The docking results revealed that the Ag₃ cluster, not Ag⁺, interacted with key amino acids of CYP2C9, CYP2C19, and CYP2D6 within a distance of about 3 Å. Moreover, the QM analysis confirmed that the amino acid residues of these CYP enzymes strongly interacted with the Ag₃ cluster, providing more insight into the mechanism of the potential inhibition of CYP enzyme activities. Interestingly, these results are consistent with previous in vitro data indicating that AgNPs inhibited activities of CYP2C and CYP2D in rat liver microsomes. It is suggested that the Ag₃ cluster is a minimal unit of AgNPs for in silico modeling. In summary, we demonstrated that molecular docking, together with QM analysis, is a promising tool to predict AgNP-mediated CYP inhibition. These methods are useful for deeper understanding of reaction mechanisms and could be used for other nanomaterials.

Conformational turn triggers regio-selectivity in the bioactivation of thiophene-contained compounds mediated by cytochrome P450

In the present work, we performed Density Functional Theory calculations to explore the bioactivation mechanism of thiophene-containing molecules mediated by P450s. For this purpose, relatively large size compounds, 2,5-diaminothiophene derivatives were selected particularly for this investigation. Here we found the successive regio-selectivity triggered by conformational turn played a significant role in the occurrence of bioactivation. 2,5-Diaminothiophene was oxidized to a 2,5-diimine thiophene-reactive intermediate by Compound I (Cpd I) through successive activations of two N-H bonds (H3-N11 and H1-N6). This reaction exhibited three special characteristics: (1) self-controlled regio-selectivity during the oxidation process. There was a large scale of conformational turn in the abstraction of the first H atom which triggers the selection of the second H for abstraction. (2) Proton-shuttle mechanism. In high spin (HS) state, proton-shuttle mechanism was observed for the abstraction of the second H atom. (3) Spin-selective manner. In protein environment, the energy barrier in HS state was much lower than that in low spin state. The novel proposed bioactivation mechanism of 2,5-diaminothiophene compounds can help us in rational design of thiophene-contained drugs avoiding the occurrence of bioactivation.

Metabolites of Vinca Alkaloid Vinblastine: Tubulin Binding and Activation of Nausea-Associated Receptors

Vinblastine (VLB) is an antimitotic drug that binds to the vinca site of tubulin. The molecule possesses a high molecular weight and a complex chemical structure with many possibilities of metabolization. Despite advances in drug discovery research in reducing drug toxicity, the cause and mechanism of VLB-induced adverse drug reactions (ADRs) remains poorly understood. VLB is metabolized to at least 35 known metabolites, which have been identified and collected in this present work. This study also explores how VLB metabolites affect nausea-associated receptors such as muscarinic, dopaminergic, and histaminic. The metabolites have stronger binding interactions than acetylcholine (ACh) for muscarinic M₁, M₄, and M₅ receptors and demonstrate similar binding profiles to that of the natural substrate, ACh. The affinities of VLB metabolites to dopaminergic and histaminic receptors, their absorption, distribution, metabolism, excretion, toxicity properties, and the superiority of VLB to ACh for binding to M₅R, indicate their potential to trigger activation of nausea-associated receptors during chemotherapy with VLB. It has been shown that metabolite 20-hydroxy-VLB (metabolite 10) demonstrates a stronger binding affinity to the vinca site of tubulin than VLB; however, they have similar modes of action. VLB and metabolite 10 have similar gastric solubility (FaSSGF), intestinal solubility (FeSSIF), and log P values. Metabolite 10 has a more acceptable pharmacokinetic profile than VLB, a better gastric and intestinal solubility. Furthermore, metabolite 10 was found to be less bound to plasma proteins than VLB. These are desired and essential features for effective drug bioavailability. Metabolite 10 is not a substrate of CYP2D6 and thus is less likely to cause drug-drug interactions and ADRs compared to its parent drug. The hydroxyl group added upon metabolism of VLB suggests that it can also be a reasonable starting compound for designing the next generation of antimitotic drugs to overcome P-glycoprotein-mediated multidrug resistance, which is often observed with vinca alkaloids.

Tris(butadiene) Metal Complexes of the First-Row Transition Metals versus Coupling of Butadiene to Eight- and Twelve-Carbon Hydrocarbon Chains

The role that zerovalent nickel plays in catalyzing the trimerization of butadiene to 1,5,9-cyclododecatriene conveys interest in the properties of the tris(butadiene)metal complexes (C₄H₆)₃M. In this connection the complexes (C₄H₆)₃M (M = Ti-Ni) of the first-row transition metals have been investigated by density functional theory. The intermediate C₁₂H₁₈Ni which has been isolated in the nickel-catalyzed trimerization of butadiene but is too unstable for X-ray crystallography is suggested here to have an open-chain hexahapto η^3,3-C₁₂H₁₈ ligand rather than the octahapto such ligand suggested by some investigators. The lowest energy (C₄H₆)₃M structures of the other first-row transition metals from vanadium to cobalt are found to have related structures with open-chain C₁₂H₁₈ ligands having hapticities ranging from four to eight with hexahapto structures being most common. The nickel and cobalt (C₁₂H₁₈)M derivatives favor low-spin singlet and doublet spin states, respectively, whereas the manganese derivative (C₁₂H₁₈)Mn favors the high-spin sextet state corresponding to the half-filled d⁵ shell of Mn(II). A (C₄H₆)₃Cr structure with three separate tetrahapto butadiene ligands analogous to the very stable (η⁴-C₄H₆)₃M (M = Mo, W) with the favored 18-electron metal configuration is found to be a very high energy structure relative to isomers containing an open-chain C₁₂H₁₈ ligand.

Mechanistic insights into a non-heme 2-oxoglutarate-dependent ethylene-forming enzyme: selectivity of ethylene-formation versusl-Arg hydroxylation

The ethylene-forming enzyme (EFE) is a unique member of the Fe(ii)- and 2-oxoglutarate-dependent (Fe/2OG) oxygenases. It converts 2OG into ethylene plus three CO2 molecules (ethylene-forming reaction) and also catalyzes the C5 hydroxylation of l-arginine coupled to the oxidative decarboxylation of 2OG (l-Arg hydroxylation reaction). To uncover the mechanisms of the dual transformations by EFE, quantum mechanical/molecular mechanical (QM/MM) calculations were carried out. Based on the results, a branched mechanism was proposed. An FeII-peroxysuccinate complex with a dissociated CO2 generated through the nucleophilic attack of the superoxo moiety of the Fe-O2 species on the keto carbon of 2OG is the key common intermediate in both reactions. A competition between the subsequent CO2 insertion (a key step in the ethylene-forming pathway) and the O-O bond cleavage (leading to the formation of succinate) governs the product selectivity. The calculated reaction barriers suggested that the CO2 insertion is favored over the O-O bond cleavage. This is consistent with the product preference observed in experiments. By comparison with the results of AsqJ (an Fe/2OG oxygenase that leads to substrate oxidation exclusively), the protein environment was found to be crucial for the selectivity. Further calculations demonstrated that the local electric field of the protein environment in EFE promotes ethylene formation by acting as a charge template, exemplifying the importance of the electrostatic interaction in enzyme catalysis. These findings offer mechanistic insights into the EFE catalysis and provide important clues for better understanding the unique ethylene-forming capability of EFE compared with other Fe/2OG oxygenases.