Modeling Enzymatic Enantioselectivity using Quantum Chemical Methodology

Asymmetric photoenzymatic incorporation of fluorinated motifs into olefins

Enzymes capable of assimilating fluorinated feedstocks are scarce. This situation poses a challenge for the biosynthesis of fluorinated compounds used in pharmaceuticals, agrochemicals, and materials. We developed a photoenzymatic hydrofluoroalkylation that integrates fluorinated motifs into olefins. The photoinduced promiscuity of flavin-dependent ene-reductases enables the generation of carbon-centered radicals from iodinated fluoroalkanes, which are directed by the photoenzyme to engage enantioselectively with olefins. This approach facilitates stereocontrol through interaction between a singular fluorinated unit and the enzyme, securing high enantioselectivity at β, γ, or δ positions of fluorinated groups through enzymatic hydrogen atom transfer-a process that is notably challenging with conventional chemocatalysis. This work advances enzymatic strategies for integrating fluorinated chemical feedstocks and opens avenues for asymmetric synthesis of fluorinated compounds.

The astrochemical evolutionary traits of phospholipid membrane homochirality

Compartmentalization is crucial for the evolution of life. Present-day phospholipid membranes exhibit a high level of complexity and species-dependent homochirality, the so-called lipid divide. It is possible that less stable, yet more dynamic systems, promoting out-of-equilibrium environments, facilitated the evolution of life at its early stages. The composition of the preceding primitive membranes and the evolutionary route towards complexity and homochirality remain unexplained. Organics-rich carbonaceous chondrites are evidence of the ample diversity of interstellar chemistry, which may have enriched the prebiotic milieu on early Earth. This Review evaluates the detections of simple amphiphiles - likely ancestors of membrane phospholipids - in extraterrestrial samples and analogues, along with potential pathways to form primitive compartments on primeval Earth. The chiroptical properties of the chiral backbones of phospholipids provide a guide for future investigations into the origins of phospholipid membrane homochirality. We highlight a plausible common pathway towards homochirality of lipids, amino acids, and sugars starting from enantioenriched monomers. Finally, given their high recalcitrance and resistance to degradation, lipids are among the best candidate biomarkers in exobiology.

Computational Study of the Fries Rearrangement Catalyzed by Acyltransferase from Pseudomonas protegens

The acyltransferase from Pseudomonas protegens (PpATase) catalyzes in nature the reversible transformation of monoacetylphloroglucinol to diacetylphloroglucinol and phloroglucinol. Interestingly, this enzyme has been shown to catalyze the promiscuous transformation of 3-hydroxyphenyl acetate to 2',4'-dihydroxyacetophenone, representing a biological version of the Fries rearrangement. In the present study, we report a mechanistic investigation of this activity of PpATase using quantum chemical calculations. A detailed mechanism is proposed, and the energy profile for the reaction is presented. The calculations show that the acylation of the enzyme is highly exothermic, while the acetyl transfer back to the substrate is only slightly exothermic. The deprotonation of the C6-H of the substrate is rate-limiting, and a remote aspartate residue (Asp137) is proposed to be the general base group in this step. Analysis of the binding energies of various acetyl acceptors shows that PpATase can promote both intramolecular and intermolecular Fries rearrangement towards diverse compounds.

Learning from Protein Engineering by Deconvolution of Multi-Mutational Variants

This review analyzes a development in biochemistry, enzymology and biotechnology that originally came as a surprise. Following the establishment of directed evolution of stereoselective enzymes in organic chemistry, the concept of partial or complete deconvolution of selective multi-mutational variants was introduced. Early deconvolution experiments of stereoselective variants led to the finding that mutations can interact cooperatively or antagonistically with one another, not just additively. During the past decade, this phenomenon was shown to be general. In some studies, molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) computations were performed in order to shed light on the origin of non-additivity at all stages of an evolutionary upward climb. Data of complete deconvolution can be used to construct unique multi-dimensional rugged fitness pathway landscapes, which provide mechanistic insights different from traditional fitness landscapes. Along a related line, biochemists have long tested the result of introducing two point mutations in an enzyme for mechanistic reasons, followed by a comparison of the respective double mutant in so-called double mutant cycles, which originally showed only additive effects, but more recently also uncovered cooperative and antagonistic non-additive effects. We conclude with suggestions for future work, and call for a unified overall picture of non-additivity and epistasis.

Quantum chemical modeling of enantioselective sulfoxidation and epoxidation reactions by indole monooxygenase VpIndA1

Indole monooxygenases (IMOs) are enzymes from the family of Group E monooxygenases, requiring flavin adenine dinucleotide (FAD) for their activities. IMOs play important roles in both sulfoxidation and epoxidation reactions. The broad substrate range and high selectivity of IMOs make them promising biocatalytic tools for synthesizing chiral compounds. In the present study, quantum chemical calculations using the cluster approach were performed to investigate the reaction mechanism and the enantioselectivity of the IMO from Variovorax paradoxus EPS (VpIndA1). The sulfoxidation of methyl phenyl sulfide (MPS) and the epoxidation of indene were chosen as the representative reactions. The calculations confirmed that the FADOOH intermediate is the catalytic species in the VpIndA1 reactions. The oxidation of MPS adopts a one-step mechanism involving the direct oxygen-transfer from FADOOH to the substrate and the proton transfer from the -OH group back to FAD, while the oxidation of indene follows a stepwise mechanism involving a carbocation intermediate. It was computationally predicted that VpIndA1 prefers the formation of (S)-product for the MPS sulfoxidation and (1S,2R)-product for the indene epoxidation, consistent with the experimental observations. Importantly, the factors controlling the stereo-preference of the two reactions are identified. The findings in the present study provide valuable insights into the VpIndA1-catalyzed reactions, which are essential for the rational design of this enzyme and other IMOs for industrial applications. It is also worth emphasizing that the quantum chemical cluster approach is again demonstrated to be powerful in studying the enantioselectivity of enzymatic reactions.

Understanding the Enzyme (S)-Norcoclaurine Synthase Promiscuity to Aldehydes and Ketones

The (S)-norcoclaurine synthase from Thalictrum flavum (TfNCS) stereoselectively catalyzes the Pictet-Spengler reaction between dopamine and 4-hydroxyphenylacetaldehyde to give (S)-norcoclaurine. TfNCS can catalyze the Pictet-Spengler reaction with various aldehydes and ketones, leading to diverse tetrahydroisoquinolines. This substrate promiscuity positions TfNCS as a highly promising enzyme for synthesizing fine chemicals. Understanding carbonyl-containing substrates' structural and electronic signatures that influence TfNCS activity can help expand its applications in the synthesis of different compounds and aid in protein optimization strategies. In this study, we investigated the influence of the molecular properties of aldehydes and ketones on their reactivity in the TfNCS-catalyzed Pictet-Spengler reaction. Initially, we compiled a library of reactive and unreactive compounds from previous publications. We also performed enzymatic assays using nuclear magnetic resonance to identify some reactive and unreactive carbonyl compounds, which were then included in the library. Subsequently, we employed QSAR and DFT calculations to establish correlations between substrate-candidate structures and reactivity. Our findings highlight correlations of structural and stereoelectronic features, including the electrophilicity of the carbonyl group, to the reactivity of aldehydes and ketones toward the TfNCS-catalyzed Pictet-Spengler reaction. Interestingly, experimental data of seven compounds out of fifty-three did not correlate with the electrophilicity of the carbonyl group. For these seven compounds, we identified unfavorable interactions between them and the TfNCS. Our results demonstrate the applications of in silico techniques in understanding enzyme promiscuity and specificity, with a particular emphasis on machine learning methodologies, DFT electronic structure calculations, and molecular dynamic (MD) simulations.

Reduction of molecular oxygen in flavodiiron proteins - Catalytic mechanism and comparison to heme-copper oxidases

The family of flavodiiron proteins (FDPs) plays an important role in the scavenging and detoxification of both molecular oxygen and nitric oxide. Using electrons from a flavin mononucleotide cofactor molecular oxygen is reduced to water and nitric oxide is reduced to nitrous oxide and water. While the mechanism for NO reduction in FDPs has been studied extensively, there is very little information available about O2 reduction. Here we use hybrid density functional theory (DFT) to study the mechanism for O2 reduction in FDPs. An important finding is that a proton coupled reduction is needed after the O2 molecule has bound to the diferrous diiron active site and before the OO bond can be cleaved. This is in contrast to the mechanism for NO reduction, where both NN bond formation and NO bond cleavage occurs from the same starting structure without any further reduction, according to both experimental and computational results. This computational result for the O2 reduction mechanism should be possible to evaluate experimentally. Another difference between the two substrates is that the actual OO bond cleavage barrier is low, and not involved in rate-limiting the reduction process, while the barrier connected with bond cleavage/formation in the NO reduction process is of similar height as the rate-limiting steps. We suggest that these results may be part of the explanation for the generally higher activity for O2 reduction as compared to NO reduction in most FDPs. Comparisons are also made to the O2 reduction reaction in the family of heme‑copper oxidases.

Stereocontrol via Propeller Chirality in FLP-Catalyzed Asymmetric Hydrogenation

Utilization of chiral frustrated Lewis pairs as catalysts in enantioselective hydrogenation of unsaturated molecules represents a promising approach in asymmetric synthesis. In our effort to improve our current understanding of the factors governing the stereoselectivity in these catalytic processes, herein we examined the mechanism of direct hydrogenation of aromatic enamines catalyzed by a binaphthyl-based chiral amino-borane. Our computational analysis reveals that only one particular conformer of the key borohydride reaction intermediate can be regarded as a reactive form of this species. This borohydride conformer has a well-defined chiral propeller shape, which induces facial selectivity in the hydride transfer to pro-chiral iminium intermediates. The propeller chirality of the reactive borohydride conformer is generated by the axially chiral binaphthyl scaffold of the amino-borane catalyst through stabilizing π-π stacking interactions. This new computational insight can be readily used to interpret the high degree of stereoinduction observed for these reactions. We expect that the concept of chirality relay could be further exploited in catalyst design endeavors.

Role of Explicit Solvation in Computational Modeling of Chemical Reactions: Mechanism of Cu(I) Transfer Between Thiolate-Based Chelators in Water

Solvation plays important roles in controlling the thermodynamic and kinetic aspects of chemical reactions. The conventional approaches to treat solvation via electronic structure methods are likely to become inadequate, when the reacting solutes have strong electrostatic and hydrogen bonding interactions with the solvent and undergo significant structural changes during the course of the reaction. In this article, we present evidence of such solvent and structural effects in the computational study of the Cu(I) transfer reaction between thiolate-based chelators dithiobutylamine (DTBA) and dithiotheritol (DTT) in water, inspired from biological copper trafficking phenomena. We propose a general solution to the problem by combining classical molecular dynamics (MD) simulations of the bulk system and static quantum chemistry calculations. The fluctuating solvation shell was estimated from MD, and energetics was assessed by averaging QM energies of a series of molecular clusters constructed from the MD snapshots. Applying this approach, we propose a reaction pathway with estimates of relative intermediate stabilities and barriers, which suggest the overall reaction to be reversible in nature and likely to go through both two and three coordinated intermediates, confirming previous studies of similar protein analogues. An interesting fact that emerged from our study was the strong indication that the rate-determining step is the deprotonation of initial thiol bound Cu(I) complex, without involving any Cu(I)-S bonds. The proposed method will lead to a better treatment of solvations, and these mechanistic insights will aid our understanding of biological copper(I) trafficking.

Enzymatic Synthesis of Noncanonical α-Amino Acids Containing γ-Tertiary Alcohols

Noncanonical amino acids (ncAAs) containing tertiary alcohols are valuable as precursors of natural products and active pharmaceutical ingredients. However, the assembly of such ncAA scaffolds from simple material by C-C bond formation remains a challenging task due to the presence of multiple stereocenters and large steric hindrance. In this study, we present a novel solution to this problem through highly selective enzymatic decarboxylative aldol addition. This method allows for the streamlined assembly of multifunctionalized ncAAs with γ-tertiary alcohols from readily available materials, such as L -aspartatic acid and isatins, vicinal diones and keto esters. The modularity of electrophiles furnished four classes of ncAAs with decent efficiency as well as excellent site and stereocontrol. Computational modeling was employed to gain detailed insight into the catalytic mechanism and to provide a rationale for the observed selectivities. The method offers a single-step approach to producing multifunctionalized ncAAs, which can be directly utilized in peptide synthesis and bioactivity assessment.

Catalytic divergencies in the mechanism of L-arginine hydroxylating nonheme iron enzymes

Many enzymes in nature utilize a free arginine (L-Arg) amino acid to initiate the biosynthesis of natural products. Examples include nitric oxide synthases, which generate NO from L-Arg for blood pressure control, and various arginine hydroxylases involved in antibiotic biosynthesis. Among the groups of arginine hydroxylases, several enzymes utilize a nonheme iron(II) active site and let L-Arg react with dioxygen and α-ketoglutarate to perform either C3-hydroxylation, C4-hydroxylation, C5-hydroxylation, or C4-C5-desaturation. How these seemingly similar enzymes can react with high specificity and selectivity to form different products remains unknown. Over the past few years, our groups have investigated the mechanisms of L-Arg-activating nonheme iron dioxygenases, including the viomycin biosynthesis enzyme VioC, the naphthyridinomycin biosynthesis enzyme NapI, and the streptothricin biosynthesis enzyme OrfP, using computational approaches and applied molecular dynamics, quantum mechanics on cluster models, and quantum mechanics/molecular mechanics (QM/MM) approaches. These studies not only highlight the differences in substrate and oxidant binding and positioning but also emphasize on electronic and electrostatic differences in the substrate-binding pockets of the enzymes. In particular, due to charge differences in the active site structures, there are changes in the local electric field and electric dipole moment orientations that either strengthen or weaken specific substrate C-H bonds. The local field effects, therefore, influence and guide reaction selectivity and specificity and give the enzymes their unique reactivity patterns. Computational work using either QM/MM or density functional theory (DFT) on cluster models can provide valuable insights into catalytic reaction mechanisms and produce accurate and reliable data that can be used to engineer proteins and synthetic catalysts to perform novel reaction pathways.

Remote stereocontrol with azaarenes via enzymatic hydrogen atom transfer

Strategies for achieving asymmetric catalysis with azaarenes have traditionally fallen short of accomplishing remote stereocontrol, which would greatly enhance accessibility to distinct azaarenes with remote chiral centres. The primary obstacle to achieving superior enantioselectivity for remote stereocontrol has been the inherent rigidity of the azaarene ring structure. Here we introduce an ene-reductase system capable of modulating the enantioselectivity of remote carbon-centred radicals on azaarenes through a mechanism of chiral hydrogen atom transfer. This photoenzymatic process effectively directs prochiral radical centres located more than six chemical bonds, or over 6 Å, from the nitrogen atom in azaarenes, thereby enabling the production of a broad array of azaarenes possessing a remote γ-stereocentre. Results from our integrated computational and experimental investigations underscore that the hydrogen bonding and steric effects of key amino acid residues are important for achieving such high stereoselectivities.

Evidence of a Distinctive Enantioselective Binding Mode for the Photoinduced Radical Cyclization of α-Chloroamides in Ene-Reductases

We demonstrate here through molecular simulations and mutational studies the origin of the enantioselectivity in the photoinduced radical cyclization of α-chloroacetamides catalyzed by ene-reductases, in particular the Gluconobacter oxidans ene-reductase and the Old Yellow Enzyme 1, which show opposite enantioselectivity. Our results reveal that neither the π-facial selectivity model nor a protein-induced selective stabilization of the transition states is able to explain the enantioselectivity of the radical cyclization in the studied flavoenzymes. We propose a new enantioinduction scenario according to which enantioselectivity is indeed controlled by transition-state stability; however, the relative stability of the prochiral transition states is not determined by direct interaction with the protein but is rather dependent on an inherent degree of freedom within the substrate itself. This intrinsic degree of freedom, distinct from the traditional π-facial exposure mode, can be controlled by the substrate conformational selection upon binding to the enzyme.

Impacts of QM region sizes and conformation numbers on modelling enzyme reactions: a case study of polyethylene terephthalate hydrolase

A quantum mechanics/molecular mechanics (QM/MM) approach is a broadly used tool in computational enzymology. Treating the QM region with a high-level DFT method is one of the important branches. Here, taking leaf-branch compost cutinase-catalyzed polyethylene terephthalate depolymerization as an example, the convergence behavior of energy barriers as well as key structural and charge features with respect to the size of the QM region (up to 1000 atoms) is systematically investigated. BP86/6-31G(d)//CHARMM and M06-2X/6-311G(d,p)//CHARMM level of theories were applied for geometry optimizations and single-point energy calculations, respectively. Six independent enzyme conformations for all the four catalytic steps (steps (i)-(iv)) were considered. Most of the twenty-four cases show that at least 500 QM atoms are needed while only two rare cases show that ∼100 QM atoms are sufficient for convergence when only a single conformation was considered. This explains why most previous studies showed that 500 or more QM atoms are required while a few others showed that ∼100 QM atoms are sufficient for DFT/MM calculations. More importantly, average energy barriers and key structural/charge features from six conformations show an accelerated convergence than that in a single conformation. For instance, to reach energy barrier convergence (within 2.0 kcal mol-1) for step (ii), only ∼100 QM atoms are required if six conformations are considered while 500 or more QM atoms are needed with a single conformation. The convergence is accelerated to be more rapid if hundreds and thousands of conformations were considered, which aligns with previous findings that only several dozens of QM atoms are required for convergence with semi-empirical QM/MM MD simulations.

Mutexa: A Computational Ecosystem for Intelligent Protein Engineering

Protein engineering holds immense promise in shaping the future of biomedicine and biotechnology. This Review focuses on our ongoing development of Mutexa, a computational ecosystem designed to enable "intelligent protein engineering". In this vision, researchers will seamlessly acquire sequences of protein variants with desired functions as biocatalysts, therapeutic peptides, and diagnostic proteins through a finely-tuned computational machine, akin to Amazon Alexa's role as a versatile virtual assistant. The technical foundation of Mutexa has been established through the development of a database that combines and relates enzyme structures and their respective functions (e.g., IntEnzyDB), workflow software packages that enable high-throughput protein modeling (e.g., EnzyHTP and LassoHTP), and scoring functions that map the sequence-structure-function relationship of proteins (e.g., EnzyKR and DeepLasso). We will showcase the applications of these tools in benchmarking the convergence conditions of enzyme functional descriptors across mutants, investigating protein electrostatics and cavity distributions in SAM-dependent methyltransferases, and understanding the role of nonelectrostatic dynamic effects in enzyme catalysis. Finally, we will conclude by addressing the future steps and fundamental challenges in our endeavor to develop new Mutexa applications that assist the identification of beneficial mutants in protein engineering.

Mechanism of CO2 Reduction to Methanol with H2 on an Iron(II)-scorpionate Catalyst

CO2 utilization is an important process in the chemical industry with great environmental power. In this work we show how CO2 and H2 can be reacted to form methanol on an iron(II) center and highlight the bottlenecks for the reaction and what structural features of the catalyst are essential for efficient turnover. The calculations predict the reactions to proceed through three successive reaction cycles that start with heterolytic cleavage of H2 followed by sequential hydride and proton transfer processes. The H2 splitting process is an endergonic process and hence high pressures will be needed to overcome this step and trigger the hydrogenation reaction. Moreover, H2 cleavage into a hydride and proton requires a metal to bind hydride and a nearby source to bind the proton, such as an amide or pyrazolyl group, which the scorpionate ligand used here facilitates. As such the computations highlight the non-innocence of the ligand scaffold through proton shuttle from H2 to substrate as an important step in the reaction mechanism.

Enantioselective High-Throughput Assay Showcased for the Identification of (R)- as well as (S)-Selective Unspecific Peroxygenases for C-H Oxidation

Identifying (bio)catalysts displaying high enantio-/stereoselectivity is a fundamental prerequisite for the advancement of asymmetric catalysis. Herein, a high-throughput, stereoselective screening assay is reported that gives information on enantioselectivity, stereopreference and activity as showcased for peroxygenase-catalyzed hydroxylation. The assay is based on spectrophotometric analysis of the simultaneous formation of NAD(P)H from the alcohol dehydrogenase catalyzed enantioselective oxidation of the sec-alcohol product formed in the peroxygenase reaction. The assay was applied to investigate a library comprising 44 unspecific peroxygenases (UPOs) containing 25 UPOs not reported yet. Thereby, previously non-described wild-type UPOs displaying (S)- as well as (R)-stereoselectivity for the hydroxylation of representative model substrates were identified, reaching up to 98 % ee for the (R)- and 94 % ee for the (S)-enantiomer. Homology models with concomitant docking studies indicated the structural reason for the observed complementary stereopreference.

Role of mutations in a chemoenzymatic enantiodivergent C(sp3)-H insertion: exploring the mechanism and origin of stereoselectivity

New-to-nature enzymes have emerged as powerful catalysts in recent years for streamlining various stereoselective organic transformations. While synthetic strategies employing engineered enzymes have witnessed proliferating success, there is limited clarity on the mechanistic front and more so when considering molecular-level insights into the role of selected mutations, dramatically escalating catalytic competency and selectivity. We have investigated the mechanism and correlation between mutations and exquisite stereoselectivity of a lactone carbene insertion into the C(sp3)-H bond of substituted aniline, catalyzed by two mutants of a cytochrome P450 variant, "P411" (engineered through directed evolution) in which the axial cysteine has been mutated to serine, utilizing various computational tools. The pivotal role of S264 and L/R328 mutations in the active site has been delineated computationally using two cluster models, thus rationalizing the enantiodivergence. This report provides much-needed insights into the origin of enantiodivergence, furnishing a mechanistic framework for understanding the anchoring effects of H-bond donor residues with the lactone ring. This study is expected to have important implications in the rational design of stereodivergent enzymes and toward successful in silico enzyme designing.

How to Compute Atomistic Insight in DFT Clusters: The REG-IQA Approach

The relative energy gradient (REG) method is paired with the topological energy partitioning method interacting quantum atoms (IQA), as REG-IQA, to provide detailed and unbiased knowledge on the intra- and interatomic interactions. REG operates on a sequence of geometries representing a dynamical change of a system. Its recent application to peptide hydrolysis of the human immunodeficiency virus-1 (HIV-1) protease (PDB code: 4HVP) has demonstrated its full potential in recovering reaction mechanisms and through-space electrostatic and exchange-correlation effects, making it a compelling tool for analyzing enzymatic reactions. In this study, the computational efficiency of the REG-IQA method for the 133-atom HIV-1 protease quantum mechanical system is analyzed in every detail and substantially improved by means of three different approaches. The first approach of smaller integration grids for IQA integrations reduces the computational overhead by about a factor of 3. The second approach uses the line-simplification Ramer-Douglas-Peucker (RDP) algorithm, which outputs the minimal number of geometries necessary for the REG-IQA analysis for a predetermined root mean squared error (RMSE) tolerance. This cuts the computational time of the whole REG analysis by a factor of 2 if an RMSE of 0.5 kJ/mol is considered. The third approach consists of a "biased" or "unbiased" selection of a specific subset of atoms of the whole initial quantum mechanical model wave-function, which results in more than a 10-fold speed-up per geometry for the IQA calculation, without deterioration of the outcome of the REG-IQA analysis. Finally, to show the capability of these approaches, the findings gathered from the HIV-1 protease system are also applied to a different system named haloalcohol dehalogenase (HheC). In summary, this study takes the REG-IQA method to a computationally feasible and highly accurate level, making it viable for the analysis of a multitude of enzymatic systems.

Uncovering the Mechanism of Azepino-Indole Skeleton Formation via Pictet-Spengler Reaction by Strictosidine Synthase: A Quantum Chemical Investigation

Strictosidine synthase (STR) catalyzes the Pictet-Spengler (PS) reaction of tryptamine and secologanin to produce strictosidine. Recent studies demonstrated that the enzyme can also catalyze the reaction of non-natural substrates to form new alkaloid skeletons. For example, the PS condensation of 1H-indole-4-ethanamine with secologanin could be promoted by the STR from Rauvolfia serpentina (RsSTR) to generate a rare class of skeletons with a seven-membered ring, namely azepino-[3,4,5-cd]-indoles, which are precursors for the synthesis of new compounds displaying antimalarial activity. In the present study, the detailed reaction mechanism of RsSTR-catalyzed formation of the rare seven-membered azepino-indole skeleton through the PS reaction was revealed at the atomic level by quantum chemical calculations. The structures of the transition states and intermediates involved in the reaction pathway were optimized, and the energetics of the complete reaction were analyzed. Based on our calculation results, the most likely pathway of the enzyme-catalyzed reaction was determined, and the rate-determining step of the reaction was clarified. The mechanistic details obtained in the present study are important in understanding the promiscuous activity of RsSTR in the formation of the rare azepino-indole skeleton molecule and are also helpful in designing STR enzymes for the synthesis of other new alkaloid skeleton molecules.

Periodicity of the Affinity of Lanmodulin for Trivalent Lanthanides and Actinides: Structural and Electronic Insights from Quantum Chemical Calculations

Lanmodulin (LanM) is the first identified macrochelator that has naturally evolved to sequester ions of rare earth elements (REEs) such as Y and all lanthanides, reversibly. This natural protein showed a 10⁶ times better affinity for lanthanide cations than for Ca, which is a naturally abundant and biologically relevant element. Recent experiments have shown that its metal ion binding activity can be further extended to some actinides, like Np, Pu, and Am. For this reason, it was thought that LanM could be adopted for the separation of REE ions and actinides, thus increasing the interest in its potential use for industry-oriented applications. In this work, a systematic study of the affinity of LanM for lanthanides and actinides has been carried out, taking into account all trivalent ions belonging to the 4f (from La to Lu) and 5f (from Ac to Lr) series, starting from their chemistry in solution. On the basis of a recently published nuclear magnetic resonance structure, a model of the LanM-binding site was built and a detailed structural and electronic description of initial aquo- and LanM-metal ion complexes was provided. The obtained binding energies are in agreement with the available experimental data. A possible reason that could explain the origin of the affinity of LanM for these metal ions is also discussed.

The Quantum Chemical Cluster Approach in Biocatalysis

ConspectusThe quantum chemical cluster approach has been used for modeling enzyme active sites and reaction mechanisms for more than two decades. In this methodology, a relatively small part of the enzyme around the active site is selected as a model, and quantum chemical methods, typically density functional theory, are used to calculate energies and other properties. The surrounding enzyme is modeled using implicit solvation and atom fixing techniques. Over the years, a large number of enzyme mechanisms have been solved using this method. The models have gradually become larger as a result of the faster computers, and new kinds of questions have been addressed. In this Account, we review how the cluster approach can be utilized in the field of biocatalysis. Examples from our recent work are chosen to illustrate various aspects of the methodology. The use of the cluster model to explore substrate binding is discussed first. It is emphasized that a comprehensive search is necessary in order to identify the lowest-energy binding mode(s). It is also argued that the best binding mode might not be the productive one, and the full reactions for a number of enzyme-substrate complexes have therefore to be considered to find the lowest-energy reaction pathway. Next, examples are given of how the cluster approach can help in the elucidation of detailed reaction mechanisms of biocatalytically interesting enzymes, and how this knowledge can be exploited to develop enzymes with new functions or to understand the reasons for lack of activity toward non-natural substrates. The enzymes discussed in this context are phenolic acid decarboxylase and metal-dependent decarboxylases from the amidohydrolase superfamily. Next, the application of the cluster approach in the investigation of enzymatic enantioselectivity is discussed. The reaction of strictosidine synthase is selected as a case study, where the cluster calculations could reproduce and rationalize the selectivities of both the natural and non-natural substrates. Finally, we discuss how the cluster approach can be used to guide the rational design of enzyme variants with improved activity and selectivity. Acyl transferase from Mycobacterium smegmatis serves as an instructive example here, for which the calculations could pinpoint the factors controlling the reaction specificity and enantioselectivity. The cases discussed in this Account highlight thus the value of the cluster approach as a tool in biocatalysis. It complements experiments and other computational techniques in this field and provides insights that can be used to understand existing enzymes and to develop new variants with tailored properties.

A Combined Computational–Experimental Study on the Substrate Binding and Reaction Mechanism of Salicylic Acid Decarboxylase

Salicylic acid decarboxylase (SDC) from the amidohydrolase superfamily (AHS) catalyzes the reversible decarboxylation of salicylic acid to form phenol. In this study, the substrate binding mode and reaction mechanism of SDC were investigated using computational and crystallographic methods. Quantum chemical calculations show that the enzyme follows the general mechanism of AHS decarboxylases. Namely, the reaction begins with proton transfer from a metal-coordinated aspartic acid residue (Asp298 in SDC) to the C1 of salicylic acid, which is followed by the C–C bond cleavage, to generate the phenol product and release CO2. Interestingly, the calculations show that SDC is a Mg-dependent enzyme rather than the previously proposed Zn-dependent, and the substrate is shown to be bidentately coordinated to the metal center in the catalysis, which is also different from the previous proposal. These predictions are corroborated by the crystal structure of SDC solved in complex with the substrate analogue 2-nitrophenol. The mechanistic insights into SDC in the present study provide important information for the rational design of the enzyme.

The Mechanism of Biochemical NO‐Sensing: Insights from Computational Chemistry

The binding of small gas molecules such as NO and CO plays a major role in the signaling routes of the human body. The sole NO‐receptor in humans is soluble guanylyl cyclase (sGC) – a histidine‐ligated heme protein, which, upon NO binding, activates a downstream signaling cascade. Impairment of NO‐signaling is linked, among others, to cardiovascular and inflammatory diseases. In the present work, we use a combination of theoretical tools such as MD simulations, high‐level quantum chemical calculations and hybrid QM/MM methods to address various aspects of NO binding and to elucidate the most likely reaction paths and the potential intermediates of the reaction. As a model system, the H‐NOX protein from Shewanella oneidensis (So H‐NOX) homologous to the NO‐binding domain of sGC is used. The signaling route is predicted to involve NO binding to form a six‐coordinate intermediate heme‐NO complex, followed by relatively facile His decoordination yielding a five‐coordinate adduct with NO on the distal side with possible isomerization to the proximal side through binding of a second NO and release of the first one. MD simulations show that the His sidechain can quite easily rotate outward into solvent, with this motion being accompanied in our simulations by shifts in helix positions that are consistent with this decoordination leading to significant conformational change in the protein.

Recognition in the Domain of Molecular Chirality: From Noncovalent Interactions to Separation of Enantiomers

It is not a coincidence that both chirality and noncovalent interactions are ubiquitous in nature and synthetic molecular systems. Noncovalent interactivity between chiral molecules underlies enantioselective recognition as a fundamental phenomenon regulating life and human activities. Thus, noncovalent interactions represent the narrative thread of a fascinating story which goes across several disciplines of medical, chemical, physical, biological, and other natural sciences. This review has been conceived with the awareness that a modern attitude toward molecular chirality and its consequences needs to be founded on multidisciplinary approaches to disclose the molecular basis of essential enantioselective phenomena in the domain of chemical, physical, and life sciences. With the primary aim of discussing this topic in an integrated way, a comprehensive pool of rational and systematic multidisciplinary information is provided, which concerns the fundamentals of chirality, a description of noncovalent interactions, and their implications in enantioselective processes occurring in different contexts. A specific focus is devoted to enantioselection in chromatography and electromigration techniques because of their unique feature as "multistep" processes. A second motivation for writing this review is to make a clear statement about the state of the art, the tools we have at our disposal, and what is still missing to fully understand the mechanisms underlying enantioselective recognition.

Reaction Mechanism of Human PAICS Elucidated by Quantum Chemical Calculations

Human PAICS is a bifunctional enzyme that is involved in the de novo purine biosynthesis, catalyzing the conversion of aminoimidazole ribonucleotide (AIR) into N-succinylcarboxamide-5-aminoimidazole ribonucleotide (SAICAR). It comprises two distinct active sites, AIR carboxylase (AIRc) where the AIR is initially converted to carboxyaminoimidazole ribonucleotide (CAIR) by reaction with CO₂ and SAICAR synthetase (SAICARs) in which CAIR then reacts with an aspartate to form SAICAR, in an ATP-dependent reaction. Human PAICS is a promising target for the treatment of various types of cancer, and it is therefore of high interest to develop a detailed understanding of its reaction mechanism. In the present work, density functional theory calculations are employed to investigate the PAICS reaction mechanism. Starting from the available crystal structures, two large models of the AIRc and SAICARs active sites are built and different mechanistic proposals for the carboxylation and phosphorylation–condensation mechanisms are examined. For the carboxylation reaction, it is demonstrated that it takes place in a two-step mechanism, involving a C–C bond formation followed by a deprotonation of the formed tetrahedral intermediate (known as isoCAIR) assisted by an active site histidine residue. For the phosphorylation–condensation reaction, it is shown that the phosphorylation of CAIR takes place before the condensation reaction with the aspartate. It is further demonstrated that the three active site magnesium ions are involved in binding the substrates and stabilizing the transition states and intermediates of the reaction. The calculated barriers are in good agreement with available experimental data.

Feasible Cluster Model Method for Simulating the Redox Potentials of Laccase CueO and Its Variant

Laccases are regarded as versatile green biocatalysts, and recent scientific research has focused on improving their redox potential for broader industrial and environmental applications. The density functional theory (DFT) quantum mechanics approach, sufficiently rigorous and efficient for the calculation of electronic structures, is conducted to better comprehend the connection between the redox potential and the atomic structural feature of laccases. According to the crystal structure of wild type laccase CueO and its variant, a truncated miniature cluster model method was established in this research. On the basic of thermodynamic cycle, the overall Gibbs free energy variations before and after the one-electron reduction were calculated. It turned out that the trends of redox potentials to increase after variant predicted by the theoretical calculations correlated well with those obtained by experiments, thereby validating the feasibility of this cluster model method for simulating the redox potentials of laccases.

Engineered P450 Atom-Transfer Radical Cyclases are Bifunctional Biocatalysts: Reaction Mechanism and Origin of Enantioselectivity

New-to-nature radical biocatalysis has recently emerged as a powerful strategy to tame fleeting open-shell intermediates for stereoselective transformations. In 2021, we introduced a novel metalloredox biocatalysis strategy that leverages the innate redox properties of the heme cofactor of P450 enzymes, furnishing new-to-nature atom-transfer radical cyclases (ATRCases) with excellent activity and stereoselectivity. Herein, we report a combined computational and experimental study to shed light on the mechanism and origins of enantioselectivity for this system. Molecular dynamics and quantum mechanics/molecular mechanics (QM/MM) calculations revealed an unexpected role of the key beneficial mutation I263Q. The glutamine residue serves as an essential hydrogen bond donor that engages with the carbonyl moiety of the substrate to promote bromine atom abstraction and enhance the enantioselectivity of radical cyclization. Therefore, the evolved ATRCase is a bifunctional biocatalyst, wherein the heme cofactor enables atom-transfer radical biocatalysis, while the hydrogen bond donor residue further enhances the activity and enantioselectivity. Unlike many enzymatic stereocontrol rationales based on a rigid substrate binding model, our computations demonstrate a high degree of rotational flexibility of the allyl moiety in an enzyme-substrate complex and succeeding intermediates. Therefore, the enantioselectivity is controlled by the radical cyclization transition states rather than the substrate orientation in ground-state complexes in the preceding steps. During radical cyclization, anchoring effects of the Q263 residue and steric interactions with the heme cofactor concurrently control the π-facial selectivity, allowing for highly enantioselective C-C bond formation. Our computational findings are corroborated by experiments with ATRCase mutants generated from site-directed mutagenesis.

Lignin peroxidase-catalyzed direct oxidation of trace organic pollutants through a long-range electron transfer mechanism: Using propranolol as an example

In this work, lignin peroxidase (LiP) was extracted for the in vitro degradation of a persistent compound (propranolol, PPN). The results showed that 94.2% of PPN was degraded at 30 U L^-1 LiP activity and 10 mg L^-1 PPN. The PPN degradation rate increased from 33.5% to 94.2% when the veratryl alcohol (VA) concentration varied from 0 to 180 µM, but decreased to 73.1% with further VA addition. This phenomenon confirmed that VA was indispensable, however, it also acted as a competitive inhibitor of PPN oxidation. Computational analysis revealed that the Trp171…iron porphyrin (TRP-FeP) path was responsible for specific substrate (e.g., VA) transformation, and another long-range electron transfer (LRET) path through His-Asp…FeP (HSP-FeP) was discovered for non-specific substrate (e.g., PPN) degradation. These two electron-transfer routes shared one catalytic center, and VA protected the enzyme from H₂O₂-dependent inactivation. The HSP-FeP path transformed PPN through single electron transfer or H abstraction mechanisms. In addition, hydroxyl radicals generated in the LiP/H₂O₂ system were involved in the hydroxylation of the PPN intermediates. Possible degradation pathways were deduced using these degradation mechanisms and mass-spectrometry analysis. The multipath degradation mechanism endowed LiP with a remarkable capacity for removing various recalcitrant pollutants in environmental remediation.

Second coordination sphere effects on the mechanistic pathways for dioxygen activation by a ferritin: involvement of a Tyr radical and the identification of a cation binding site

Ferritins are ubiquitous diiron enzymes involved in iron(II) detoxification and oxidative stress responses and can act as metabolic iron stores. The overall reaction mechanisms of ferritin enzymes are still unclear, particularly concerning the role of the conserved, near catalytic center Tyr residue. Thus, we carried out a computational study of a ferritin using a large cluster model of well over 300 atoms including its first- and second-coordination sphere. The calculations reveal important insight into the structure and reactivity of ferritins. Specifically, the active site Tyr residue delivers a proton and electron in the catalytic cycle prior to iron(II) oxidation. In addition, the calculations highlight a likely cation binding site at Asp65, which through long-range electrostatic interactions, influences the electronic configuration and charge distributions of the metal center. The results are consistent with experimental observations but reveal novel detail of early mechanistic steps that lead to an unusual mixed-valent iron(III)-iron(II) center.

Successes and challenges in multiscale modelling of artificial metalloenzymes: the case study of POP-Rh2 cyclopropanase

Molecular modelling applications in metalloenzyme design are still scarce due to a series of challenges. On top of that, the simulations of metal-mediated binding and the identification of catalytic competent geometries require both large conformational exploration and simulation of fine electronic properties. Here, we demonstrate how the incorporation of new tools in multiscale strategies, namely substrate diffusion exploration, allows taking a step further. As a showcase, the enantioselective profiles of the most outstanding variants of an artificial Rh₂-based cyclopropanase (GSH, HFF and RFY) developed by Lewis and co-workers (Nat. Commun., 2015, 6, 7789 and Nat. Chem., 2018, 10, 318-324) have been rationalized. DFT calculations on the free-cofactor-mediated process identify the carbene insertion and the cyclopropanoid formation as crucial events, the latter being the enantiodetermining step, which displays up to 8 competitive orientations easily altered by the protein environment. The key intermediates of the reaction were docked into the protein scaffold showing that some mutated residues have direct interaction with the cofactor and/or the co-substrate. These interactions take the form of a direct coordination of Rh in GSH and HFF and a strong hydrophobic patch with the carbene moiety in RFY. Posterior molecular dynamics sustain that the cofactor induces global re-arrangements of the protein. Finally, massive exploration of substrate diffusion, based on the GPathFinder approach, defines this event as the origin of the enantioselectivity in GSH and RFY. For HFF, fine molecular dockings suggest that it is likely related to local interactions upon diffusion. This work shows how modelling of long-range mutations on the catalytic profiles of metalloenzymes may be unavoidable and software simulating substrate diffusion should be applied.

β-Methyltryptamine Provoking the Crucial Role of Strictosidine Synthase Tyr151-OH for Its Stereoselective Pictet-Spengler Reactions to Tryptoline-type Alkaloids

Strictosidine synthase (STR), the gate enzyme for monoterpenoid indole alkaloid biosynthesis, catalyzes the Pictet-Spengler reaction (PSR) of various tryptamine derivatives with secologanin assisted by "indole sandwich" stabilization. Continuous exploration with β-methyltryptamine (IPA) stereoselectively delivered the C6-methylstrictosidines and C6-methylvincosides by enzymatic and nonenzymatic PSR, respectively. Unexpectedly, the first "nonindole sandwich" binding mode was witnessed by the X-ray structures of STR1-ligand complexes. Site-directed mutagenesis revealed the critical cryptic role of the hydroxyl group of Tyr151 in IPA biotransformation. Further computational calculations demonstrated the adjustable IPA position in STR1 upon the binding of secologanin, and Tyr151-OH facilitates the productive PSR binding mode via an advantageous hydrogen-bond network. Further chemo-enzymatic manipulation of C6-methylvincosides successfully resulted in the discovered antimalarial framework (IC50 = 0.92 μM).

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.

Unraveling the mechanism of enantio-controlling switches of an alcohol dehydrogenase toward sterically small ketone

Efficient synthesis of chiral compounds under mild conditions is highly desirable in the chemical and pharmaceutical communities, but it often faces difficulties. Although various enzymes have been harnessed as biocatalysts...

A mono-oxo-bridged binuclear iron(iii) complex with a Fe–O–Fe angle of 180.0° and its catalytic activity for hydrogen evolution

There are potential applications for binuclear oxo-bridged iron complexes in biological processes, catalysts, and magnetic materials, and the design of new bridged-iron complexes is important.

Computational Study of Mechanism and Enantioselectivity of Imine Reductase from Amycolatopsis orientalis

Imine reductases (IREDs) are NADPH-dependent enzymes (NADPH=nicotinamide adenine dinucleotide phosphate) that catalyze the reduction of imines to amines. They exhibit high enantioselectivity for a broad range of substrates, making them of interest for biocatalytic applications. In this work, we have employed density functional theory (DFT) calculations to elucidate the reaction mechanism and the origins of enantioselectivity of IRED from Amycolatopsis orientalis. Two substrates are considered, namely 1-methyl-3,4-dihydroisoquinoline and 2-propyl-piperideine. A model of the active site is built on the basis of the available crystal structure. For both substrates, different binding modes are first evaluated, followed by calculation of the hydride transfer transition states from each complex. We have also investigated the effect of mutations of certain important active site residues (Tyr179Ala and Asn241Ala) on the enantioselectivity. The calculated energies are consistent with the experimental observations and the analysis of transition states geometries provides insights into the origins of enantioselectivity of this enzyme.

Biodegradation of herbicides by a plant nonheme iron dioxygenase: mechanism and selectivity of substrate analogues

Aryloxyalkanoate dioxygenases are unique herbicide biodegrading nonheme iron enzymes found in plants and hence, from environmental and agricultural point of view they are important and valuable. However, they often are substrate specific and little is known on the details of the mechanism and the substrate scope. To this end, we created enzyme models and calculate the mechanism for 2,4-dichlorophenoxyacetic acid biodegradation and 2-methyl substituted analogs by density functional theory. The work shows that the substrate binding is tight and positions the aliphatic group close to the metal center to enable a chemoselective reaction mechanism to form the C 2 -hydroxy products, whereas the aromatic hydroxylation barriers are well higher in energy. Subsequently, we investigated the metabolism of R - and S -methyl substituted inhibitors and show that these do not react as efficiently as 2,4-dichlorophenoxyacetic acid substrate due to stereochemical clashes in the active site and particularly for the R -isomer give high rebound barriers.

Computational Prediction of ω-Transaminase Specificity by a Combination of Docking and Molecular Dynamics Simulations

ω-Transaminases (ω-TAs) catalyze the conversion of ketones to chiral amines, often with high enantioselectivity and specificity, which makes them attractive for industrial production of chiral amines. Tailoring ω-TAs to accept non-natural substrates is necessary because of their limited substrate range. We present a computational protocol for predicting the enantioselectivity and catalytic selectivity of an ω-TA from Vibrio fluvialis with different substrates and benchmark it against 62 compounds gathered from the literature. Rosetta-generated complexes containing an external aldimine intermediate of the transamination reaction are used as starting conformations for multiple short independent molecular dynamics (MD) simulations. The combination of molecular docking and MD simulations ensures sufficient and accurate sampling of the relevant conformational space. Based on the frequency of near-attack conformations observed during the MD trajectories, enantioselectivities can be quantitatively predicted. The predicted enantioselectivities are in agreement with a benchmark dataset of experimentally determined ee% values. The substrate-range predictions can be based on the docking score of the external aldimine intermediate. The low computational cost required to run the presented framework makes it feasible for use in enzyme design to screen thousands of enzyme variants.

The Reactivity of Human and Equine Estrogen Quinones towards Purine Nucleosides

Conjugated estrogen medicines, which are produced from the urine of pregnant mares for the purpose of menopausal hormone replacement therapy (HRT), contain the sulfate conjugates of estrone, equilin, and equilenin in varying proportions. The latter three steroid sex hormones are highly similar in molecular structure as they only differ in the degree of unsaturation of the sterane ring “B”: the cyclohexene ring in estrone (which is naturally present in both humans and horses) is replaced by more symmetrical cyclohexadiene and benzene rings in the horse-specific (“equine”) hormones equilin and equilenin, respectively. Though the structure of ring “B” has only moderate influence on the estrogenic activity desired in HRT, it might still significantly affect the reactivity in potential carcinogenic pathways. In the present theoretical study, we focus on the interaction of estrogen orthoquinones, formed upon metabolic oxidation of estrogens in breast cells with purine nucleosides. This multistep process results in a purine base loss in the DNA chain (depurination) and the formation of a “depurinating adduct” from the quinone and the base. The point mutations induced in this manner are suggested to manifest in breast cancer development in the long run. We examine six reactions between deoxyadenosine and deoxyguanosine as nucleosides and estrone-3,4-quinone, equilin-3,4-quinone, and equilenin-3,4-quinone as mutagens. We performed DFT calculations to determine the reaction mechanisms and establish a structure–reactivity relationship between the degree of unsaturation of ring “B” and the expected rate of DNA depurination. As quinones might be present in the cytosol in various protonated forms, we introduce the concept of “effective barriers” to account for the different reactivity and different concentrations of quinone derivatives. According to our results, both equine estrogens have the potential to facilitate depurination as the activation barrier of one of the elementary steps (the initial Michael addition in the case of equilenin and the rearomatization step in the case of equilin) significantly decreases compared to that of estrone. We conclude that the appearance of exogenous equine estrogen quinones due to HRT might increase the risk of depurination-induced breast cancer development compared to the exposure to endogenous estrone metabolites. Still, further studies are required to identify the rate-limiting step of depurination under intracellular conditions to reveal whether the decrease in the barriers affects the overall rate of carcinogenesis.

Computer-aided understanding and engineering of enzymatic selectivity

Enzymes offering chemo-, regio-, and stereoselectivity enable the asymmetric synthesis of high-value chiral molecules. Unfortunately, the drawback that naturally occurring enzymes are often inefficient or have undesired selectivity toward non-native substrates hinders the broadening of biocatalytic applications. To match the demands of specific selectivity in asymmetric synthesis, biochemists have implemented various computer-aided strategies in understanding and engineering enzymatic selectivity, diversifying the available repository of artificial enzymes. Here, given that the entire asymmetric catalytic cycle, involving precise interactions within the active pocket and substrate transport in the enzyme channel, could affect the enzymatic efficiency and selectivity, we presented a comprehensive overview of the computer-aided workflow for enzymatic selectivity. This review includes a mechanistic understanding of enzymatic selectivity based on quantum mechanical calculations, rational design of enzymatic selectivity guided by enzyme-substrate interactions, and enzymatic selectivity regulation via enzyme channel engineering. Finally, we discussed the computational paradigm for designing enzyme selectivity in silico to facilitate the advancement of asymmetric biosynthesis.

Electrostatic Perturbations from the Protein Affect C-H Bond Strengths of the Substrate and Enable Negative Catalysis in the TmpA Biosynthesis Enzyme

The nonheme iron dioxygenase 2-(trimethylammonio)-ethylphosphonate dioxygenase (TmpA) is an enzyme involved in the regio- and chemoselective hydroxylation at the C¹ -position of the substrate as part of the biosynthesis of glycine betaine in bacteria and carnitine in humans. To understand how the enzyme avoids breaking the weak C² -H bond in favor of C¹ -hydroxylation, we set up a cluster model of 242 atoms representing the first and second coordination sphere of the metal center and substrate binding pocket, and investigated possible reaction mechanisms of substrate activation by an iron(IV)-oxo species by density functional theory methods. In agreement with experimental product distributions, the calculations predict a favorable C¹ -hydroxylation pathway. The calculations show that the selectivity is guided through electrostatic perturbations inside the protein from charged residues, external electric fields and electric dipole moments. In particular, charged residues influence and perturb the homolytic bond strength of the C¹ -H and C² -H bonds of the substrate, and strongly strengthens the C² -H bond in the substrate-bound orientation.