External Electric Field Will Control the Selectivity of Enzymatic-Like Bond Activations

Quantitative Insight into the Electric Field Effect on CO2 Electrocatalysis via Machine Learning Spectroscopy

During chemical reactions, especially for electrocatalysis and electrosynthesis, the electric field is the most central driving force to regulate the reaction process. However, due to the difficulty of quantitatively measuring the electric field effects caused at the microscopic level, the regulation of electrocatalytic reactions by electric fields has not been well digitally understood yet. Herein, we took the infrared/Raman spectral signals of CO2 molecules as descriptors to quantitatively predict the effects of different electric fields on the catalytic properties. Taking the metal-doped graphitic C3N4 (g-C3N4) catalyst as an example, we theoretically investigated the adsorption mode and energy of CO2 molecules adsorbed on 27 distinct metal single-atom catalysts under different directions and intensities of electric field. Through a machine learning approach, a spectroscopy-property model between infrared/Raman spectral descriptors and adsorption energy/charge transfer was established, which quantified the facilitation of electric field effects on the CO2 catalytic conversion. Meanwhile, based on the attention mechanism, the catalytic insight of the relationship between spectra and adsorption modes was mined, and the inverse prediction of electric field strength from spectra was realized. This work opens a new quantitative pathway for monitoring and regulating electrocatalytic reactions using machine learning spectroscopy.

Toward the Electrostatic Catalysis of Nucleophilic Substitutions: A Surface Chemistry Study of the Menshutkin Reaction

The catalysis of nonredox reactions by external electric fields is one of the most rapidly expanding areas of chemistry. The Menshutkin reaction, a classic example of bimolecular nucleophilic substitution (SN2), involves the conversion of a tertiary amine to a quaternary ammonium salt by coupling it with an alkyl halide. The reaction barrier of the Menshutkin reaction is theoretically predicted to be highly sensitive to the magnitude and direction of an external electric field experienced by the transition state. In this study, we investigate how near-surface electric fields can drive this prototypical nucleophilic substitution by examining the coupling of a diffusive redox-tagged tertiary amine with an electrode-tethered alkyl bromide under a variable external bias. Our findings reveal a competition between electrostatically assisted reactions, solvent effects, and electrochemically triggered side reactions involving radical intermediates. We estimate that only about 5% of the coupling events are attributable to the external field, while the majority of the reaction products originate from electrochemically generated radical intermediates.

What is the Origin of the Regioselective C3-Hydroxylation of L-Arg by the Nonheme Iron Enzyme Capreomycin C?

The nonheme iron dioxygenase capreomycin C (CmnC) hydroxylates a free L-arginine amino acid regio- and stereospecifically at the C3-position as part of the capreomycin antibiotics biosynthesis. Little is known on its structure, catalytic cycle and substrate specificity and, therefore, a comprehensive computational study was performed. A large QM cluster model of CmnC was created of 297 atoms and the mechanisms for C3-H, C4-H and C5-H hydroxylation and C3-C4 desaturation were investigated. All low-energy pathways correspond to radical reaction mechanisms with an initial hydrogen atom abstraction followed by OH rebound to form alcohol product complexes. The work is compared to alternative L-Arg hydroxylating nonheme iron dioxygenases and the differences in active site polarity are compared. We show that a tight hydrogen bonding network in the substrate binding pocket positions the substrate in an ideal orientation for C3-H activation, whereby the polar groups in the substrate binding pocket induce an electric field effect that guides the selectivity.

Electric fields imbue enzyme reactivity by aligning active site fragment orbitals

It is broadly recognized that intramolecular electric fields, produced by the protein scaffold and acting on the active site, facilitate enzymatic catalysis. This field effect can be described by several theoretical models, each of which is intuitive to varying degrees. In this contribution, we show that a fundamental effect of electric fields is to generate electrostatic potentials that facilitate the energetic alignment of reactant frontier orbitals. We apply this model to demystify the impact of electric fields on high-valent iron-oxo heme proteins: catalases, peroxidases, and peroxygenases/monooxygenases. Specifically, we show that this model easily accounts for the observed field-induced changes to the spin distribution within peroxidase active sites and explains the transition between epoxidation and hydroxylation pathways seen in Cytochrome P450 active site models. Thus, for the intuitive interpretation of the chemical effect of the field, the strategy involves analyzing the response of the orbitals of active site fragments, and their energetic alignment. We note that the energy difference between fragment orbitals involved in charge redistribution acts as a measure for the chemical hardness/softness of the reactive complex. This measure, and its sensitivity to electric fields, offers a single parameter model from which to quantitatively assess the effects of electric fields on reactivity and selectivity. Thus, the model provides an additional perspective to describe electrostatic preorganization and offers ways for its manipulation.

Electrostatic Preorganization in Three Distinct Heterogeneous Proteasome β-Subunits

The origin of the enzyme's powerful role in accelerating chemical reactions is one of the most critical and still widely discussed questions. It is already accepted that enzymes impose an electrostatic field onto their substrates by adopting complex three-dimensional structures; therefore, the preorganization of electric fields inside protein active sites has been proposed as a crucial contributor to catalytic mechanisms and rate constant enhancement. In this work, we focus on three catalytically active β-subunits of 20S proteasomes with low sequence identity (∼30%) whose active sites, although situated in an electrostatically miscellaneous environment, catalyze the same chemical reaction with similar catalytic efficiency. Our in silico experiments reproduce the experimentally observed equivalent reactivity of the three sites and show that obliteration of the electrostatic potential in all active sites would deprive the enzymes of their catalytic power by slowing down the chemical process by a factor of 1035. To regain enzymatic efficiency, besides catalytic Thr1 and Lys33 residues, the presence of aspartic acid in position 17 and an aqueous solvent is required, proving that the electrostatic potential generated by the remaining residues is insignificant for catalysis. Moreover, it was found that the gradual decay of atomic charges on Asp17 strongly correlates with the enzyme's catalytic rate deterioration as well as with a change in the charge distributions due to introduced mutations. The computational procedure used and described here may help identify key residues for catalysis in other biomolecular systems and consequently may contribute to the process of designing enzyme-like synthetic catalysts.

Evaluating the Ability of External Electric Fields to Accelerate Reactions in Solution

This study investigates the catalytic effects of external electric fields (EEFs) on two reactions in solution: the Menshutkin reaction and the Chapman rearrangement. Utilizing a scanning tunneling microscope-based break-junction (STM-BJ) setup and monitoring reaction rates through high-performance liquid chromatography connected to a UV detector (HPLC-UV) and ultraperformance liquid chromatography coupled with quadrupole time-of-flight mass spectrometry (UPLC-q-ToF-MS), we observed no rate enhancement for either reaction under ambient conditions. Density functional theory (DFT) calculations indicate that electric field-induced changes in reactant orientation and the minimization of activation energy are crucial factors in determining the efficacy of EEF-driven catalysis. Our findings suggest that the current experimental setups and field strengths are insufficient to catalyze these reactions, underscoring the importance of these criteria in assessing the reaction candidates.

Measuring the Electric Fields of Ions Captured in Crown Ethers

Controlling reactivity with electric fields is a persistent challenge in chemistry. One approach is to tether ions at well-defined locations near a reactive center. To quantify fields arising from ions, we report crown ethers that capture metal cations as field sources and a covalently bound vibrational Stark shift probe as a field sensor. We use experiments and computations in both the gas and liquid phases to quantify the vibrational frequencies of the probe and estimate the electric fields from the captured ions. Cations, in general, blue shift the probe frequency, with effective fields estimated to vary in the range of ∼0.2-3 V/nm in the liquid phase. Comparison of the gas and liquid phase data provides insight into the effects of mutual polarization of the molecule and solvent and screening of the ion's field. These findings reveal the roles of charge, local screening, and geometry in the design of tailored electric fields.

Water Dimer under Electric Fields: An Ab Initio Investigation up to Quantum Accuracy

It is well-established that strong electric fields (EFs) can align water dipoles, partially order the H-bond network of liquid water, and induce water splitting and proton transfers. To illuminate the fundamental behavior of water under external EFs, we present the first benchmark, to the best of our knowledge, of DFT calculations of the water dimer exposed to intense EFs against coupled cluster calculations. The analyses of the vibrational Stark effect and electron density provide a consistent picture of the intermolecular charge transfer effects driven along the H-bond by the increasing applied field at all theory levels. However, our findings prove that at extreme field regimes (∼1-2 V/Å) DFT calculations significantly exaggerate by ∼10-30% the field-induced strengthening of the H-bond, both within the GGA, hybrid GGA, and hybrid meta-GGA approximations. Notably, a linear correlation emerges between the vibrational Stark effect on OH stretching and H-bond strengthening: a 1 kcal mol-1 increase corresponds to an 80 cm-1 red-shift in OH stretching frequency.

My Vision of Electric-Field-Aided Chemistry in 2050

This manuscript outlines my outlook on the development of electric-field (EF)-mediated-chemistry and the vision of its state by 2050. I discuss applications of oriented-external electric-fields (OEEFs) on chemical reactions and proceed with relevant experimental verifications. Subsequently, the Perspective outlines other ways of generating EFs, e.g., by use of pH-switchable charges, ionic additives, water droplets, and so on. A special section summarizes conceptual principles for understanding and predicting OEEF effects, e.g., the "reaction-axis rule", the capability of OEEFs to act as tweezers that orient reactants and accelerate their reaction, etc. Finally, I discuss applications of OEEFs in continuous-flow setups, which may, in principle, scale-up to molar concentrations. The Perspective ends with the vision that by 2050, OEEF usage will change chemical education, if not also the art of making new molecules.

Optimal Oriented External Electric Fields to Trigger a Barrierless Oxaphosphetane Ring Opening Step of the Wittig Reaction

The Wittig reaction is one of the most important processes in organic chemistry for the asymmetric synthesis of olefinic compounds. In view of the increasingly acknowledged potentiality of the electric fields in promoting reactions, here we will consider the effect of the oriented external electric field (OEEF) on the second step of Wittig reaction (i. e. the ring opening oxaphosphetane) in a model system for non-stabilized ylides. In particular, we have determined the optimal direction and strength of the electric field that should be applied to annihilate the reaction barrier of the ring opening through the polarizable molecular electric dipole (PMED) model that we have recently developed. We conclude that the application of the optimal external electric field for the oxaphosphetane ring opening favours a Bestmann-like mechanism.

An Active Site Tyr Residue Guides the Regioselectivity of Lysine Hydroxylation by Nonheme Iron Lysine-4-hydroxylase Enzymes through Proton-Coupled Electron Transfer

Lysine dioxygenase (KDO) is an important enzyme in human physiology involved in bioprocesses that trigger collagen cross-linking and blood pressure control. There are several KDOs in nature; however, little is known about the factors that govern the regio- and stereoselectivity of these enzymes. To understand how KDOs can selectively hydroxylate their substrate, we did a comprehensive computational study into the mechanisms and features of 4-lysine dioxygenase. In particular, we selected a snapshot from the MD simulation on KDO5 and created large QM cluster models (A, B, and C) containing 297, 312, and 407 atoms, respectively. The largest model predicts regioselectivity that matches experimental observation with rate-determining hydrogen atom abstraction from the C4-H position, followed by fast OH rebound to form 4-hydroxylysine products. The calculations show that in model C, the dipole moment is positioned along the C4-H bond of the substrate and, therefore, the electrostatic and electric field perturbations of the protein assist the enzyme in creating C4-H hydroxylation selectivity. Furthermore, an active site Tyr233 residue is identified that reacts through proton-coupled electron transfer akin to the axial Trp residue in cytochrome c peroxidase. Thus, upon formation of the iron(IV)-oxo species in the catalytic cycle, the Tyr233 phenol loses a proton to the nearby Asp179 residue, while at the same time, an electron is transferred to the iron to create an iron(III)-oxo active species. This charged tyrosyl residue directs the dipole moment along the C4-H bond of the substrate and guides the selectivity to the C4-hydroxylation of the substrate.

QM/MM Study Into the Mechanism of Oxidative C=C Double Bond Cleavage by Lignostilbene-α,β-Dioxygenase

The enzymatic biosynthesis of fragrance molecules from lignin fragments is an important reaction in biotechnology for the sustainable production of fine chemicals. In this work we investigated the biosynthesis of vanillin from lignostilbene by a nonheme iron dioxygenase using QM/MM and tested several suggested proposals via either an epoxide or dioxetane intermediate. Binding of dioxygen to the active site of the protein results in the formation of an iron(II)-superoxo species with lignostilbene cation radical. The dioxygenase mechanism starts with electrophilic attack of the terminal oxygen atom of the superoxo group on the central C=C bond of lignostilbene, and the second-coordination sphere effects in the substrate binding pocket guide the reaction towards dioxetane formation. The computed mechanism is rationalized with thermochemical cycles and valence bond schemes that explain the electron transfer processes during the reaction mechanism. Particularly, the polarity of the protein and the local electric field and dipole moments enable a facile electron transfer and an exergonic dioxetane formation pathway.

Computational Study Into the Oxidative Ring-Closure Mechanism During the Biosynthesis of Deoxypodophyllotoxin

The nonheme iron dioxygenase deoxypodophyllotoxin synthase performs an oxidative ring-closure reaction as part of natural product synthesis in plants. How the enzyme enables the oxidative ring-closure reaction of (-)-yatein and avoids substrate hydroxylation remains unknown. To gain insight into the reaction mechanism and understand the details of the pathways leading to products and by-products we performed a comprehensive computational study. The work shows that substrate is bound tightly into the substrate binding pocket with the C7'-H bond closest to the iron(IV)-oxo species. The reaction proceeds through a radical mechanism starting with hydrogen atom abstraction from the C7'-H position followed by ring-closure and a final hydrogen transfer to form iron(II)-water and deoxypodophyllotoxin. Alternative mechanisms including substrate hydroxylation and an electron transfer pathway were explored but found to be higher in energy. The mechanism is guided by electrostatic perturbations of charged residues in the second-coordination sphere that prevent alternative pathways.

Biosystem-Inspired Engineering of Nanozymes for Biomedical Applications

Nanozymes with intrinsic enzyme-mimicking activities have shown great potential to become surrogates of natural enzymes in many fields by virtue of their advantages of high catalytic stability, ease of functionalization, and low cost. However, due to the lack of predictable descriptors, most of the nanozymes reported in the past have been obtained mainly through trial-and-error strategies, and the catalytic efficacy, substrate specificity, as well as practical application effect under physiological conditions, are far inferior to that of natural enzymes. To optimize the catalytic efficacies and functions of nanozymes in biomedical settings, recent studies have introduced biosystem-inspired strategies into nanozyme design. In this review, recent advances in the engineering of biosystem-inspired nanozymes by leveraging the refined catalytic structure of natural enzymes, simulating the behavior changes of natural enzymes in the catalytic process, and mimicking the specific biological processes or living organisms, are introduced. Furthermore, the currently involved biomedical applications of biosystem-inspired nanozymes are summarized. More importantly, the current opportunities and challenges of the design and application of biosystem-inspired nanozymes are discussed. It is hoped that the studies of nanozymes based on bioinspired strategies will be beneficial for constructing the new generation of nanozymes and broadening their biomedical applications.

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.

Designed Local Electric Fields-Promising Tools for Enzyme Engineering

Designing efficient catalysts is one of the ultimate goals of chemists. In this Perspective, we discuss how local electric fields (LEFs) can be exploited to improve the catalytic performance of supramolecular catalysts, such as enzymes. More specifically, this Perspective starts by laying out the fundamentals of how local electric fields affect chemical reactivity and review the computational tools available to study electric fields in various settings. Subsequently, the advances made so far in optimizing enzymatic electric fields through targeted mutations are discussed critically and concisely. The Perspective ends with an outlook on some anticipated evolutions of the field in the near future. Among others, we offer some pointers on how the recent data science/machine learning revolution, engulfing all science disciplines, could potentially provide robust and principled tools to facilitate rapid inference of electric field effects, as well as the translation between optimal electrostatic environments and corresponding chemical modifications.

External electric field, a potential catalyst for C-N cross-coupling reaction

The present work investigates the role of the external electric field (EEF) in boosting the C-N cross-coupling reaction between 2-chlorobenzoic acid and propylamine, by computing the reaction rates and energy barrier. The investigation suggests that the reaction can become barrierless by choosing an electric field in the appropriate direction, resulting in a quadrillionfold increase in the reaction rate in the presence of an EEF. We also found that the efficiency of the electric field depends on the dipole moment of the reactants, and hence, the results of the present work are general in nature and should be applicable to a variety of C-N cross-coupling reactions.

Methane Activation by [AlFeO3 ]+ : the Hidden Spin Selectivity

The performance of heteronuclear cluster [AlFeO3 ]+ in activating methane has been explored by a combination of high-level quantum chemical calculations with gas-phase experiments. At room temperature, [AlFeO3 ]+ is a mixture of 7 [AlFeO3 ]+ and 5 [AlFeO3 ]+ , in which two states lead to different reactivity and chemoselectivity for methane activation. While hydrogen extracted from methane is the only product channel for the 7 [AlFeO3 ]+ /CH4 couple, 5 [AlFeO3 ]+ is able to convert this substrate to formaldehyde. In addition, the introduction of an external electric field may regulate the reactivity and product selectivity. The interesting doping effect of Fe and the associated electronic origins are discussed, which may guide one on the design of Fe-involved catalyst for methane conversion.

High Electric Fields on Water Microdroplets Catalyze Spontaneous and Fast Reactions in Halogen-Bond Complexes

The use of external electric fields as green and efficient catalysts in synthetic chemistry has recently received significant attention for their ability to deliver remarkable control of reaction selectivity and acceleration of reaction rates. Technically, methods of generating high electric fields in the range of 1-10 V/nm are limited, as in-vacuo techniques have obvious scalability issues. The spontaneous high fields at various interfaces promise to solve this problem. In this study, we take advantage of the spontaneous high electric field at the air-water interface of sprayed water microdroplets in the reactions of several halogen bond systems: Nu:--X-X, where Nu: is pyridine or quinuclidine and X is bromine or iodine. The field facilitates ultrafast electron transfer from Nu:, yielding a Nu-X covalent bond and causing the X-X bond to cleave. This reaction occurs in microseconds in microdroplets but takes days to weeks in bulk solution. Density functional theory calculations predict that the reaction becomes barrier-free in the presence of oriented external electric fields, supporting the notion that the electric fields in the water droplets are responsible for the catalysis. We anticipate that microdroplet chemistry will be an avenue rich in opportunities in the reactions facilitated by high electric fields and provides an alternative way to tackle the scalability problem.

An algorithm to find the optimal oriented external electrostatic field for annihilating a reaction barrier in a polarizable molecular system

The use of oriented external electric fields (OEEFs) to promote and control chemical reactivity has motivated many theoretical and computational studies in the last decade to model the action of OEEFs on a molecular system and its effects on chemical processes. Given a reaction, a central goal in this research area is to predict the optimal OEEF (oOEEF) required to annihilate the reaction energy barrier with the smallest possible field strength. Here, we present a model rooted in catastrophe and optimum control theories that allows us to find the oOEEF for a given reaction valley in the potential energy surface (PES). In this model, the effective (or perturbed) PES of a polarizable molecular system is constructed by adding to the original, non-perturbed, PES a term accounting for the interaction of the OEEF with the intrinsic electric dipole and polarizability of the molecular system, so called the polarizable molecular electric dipole (PMED) model. We demonstrate that the oOEEF can be established by locating a point in the original PES with unique topological properties: the optimal barrier breakdown or bond-breaking point (oBBP). The essential feature of the oBBP structure is the fact that this point maintains its topological properties for all the applied OEEFs, also for the unperturbed PES, thus becoming much more relevant than the commonly used minima and transition state structures. The PMED model proposed here has been implemented in an open access package and is shown to successfully predict the oOEEF for two processes: an isomerization reaction of a cumulene derivative and the Huisgen cycloaddition reaction.

Impact of Static-Oriented Electric Fields on the Kinetics of Some Representative Suzuki-Miyaura and Metal-Cluster Mediated Reactions

In order to examine the effect of oriented (static) electric fields (OEF) on the kinetics of some representative Suzuki-Miyaura and metal-cluster mediated reactions at ambient temperatures, density functional theory-based calculations are reported herein. Results indicate that, in general, OEF can facilitate the kinetics of the concerned reactions when applied along the suitable direction (parallel or anti-parallel with respect to the reaction axis). The reverse effect happens if the direction of the OEF is flipped. OEF (when applied along the 'right' direction) helps to polarize the transition states in the desired direction, thereby facilitating favorable bonding interactions. Given the growing need for finding appropriate catalysts among the scientific community, OEF can prove to be a vital route for the same.

Two-State Hydrogen Atom Transfer Reactivity of Unsymmetric [Cu2(O)(NO)]2+ Complexes

We report the temperature-dependent spin switching of dicopper oxo nitrosyl [Cu2(O)(NO)]2+ complexes and their influence on hydrogen atom transfer (HAT) reactivity. Electron paramagnetic resonance (EPR) and Evans method analysis suggest that [Cu2(O)(NO)]2+ complexes transition from the S = 1/2 to the S = 3/2 state around ca. 202 K. At low temperatures (198 K) where S = 3/2 dominates, a strong correlation between the rate of HAT (kHAT) and the population of the S = 1/2 state was identified (R2 = 0.988), suggesting that the HAT by [Cu2(O)(NO)]2+ complexes proceeds by the S = 1/2 isomer. Installation of functional groups that introduce an unsymmetric secondary coordination environment accelerates the HAT rates through perturbation of the spin equilibria. Given the often unsymmetric coordination sphere of bimetallic active sites in natural proteins, we anticipate that similar strategies could be employed by metalloenzymes to control HAT reactions.

How Is Substrate Halogenation Triggered by the Vanadium Haloperoxidase from Curvularia inaequalis?

Vanadium haloperoxidases (VHPOs) are unique enzymes in biology that catalyze a challenging halogen transfer reaction and convert a strong aromatic C-H bond into C-X (X = Cl, Br, I) with the use of a vanadium cofactor and H₂O₂. The VHPO catalytic cycle starts with the conversion of hydrogen peroxide and halide (X = Cl, Br, I) into hypohalide on the vanadate cofactor, and the hypohalide subsequently reacts with a substrate. However, it is unclear whether the hypohalide is released from the enzyme or otherwise trapped within the enzyme structure for the halogenation of organic substrates. A substrate-binding pocket has never been identified for the VHPO enzyme, which questions the role of the protein in the overall reaction mechanism. Probing its role in the halogenation of small molecules will enable further engineering of the enzyme and expand its substrate scope and selectivity further for use in biotechnological applications as an environmentally benign alternative to current organic chemistry synthesis. Using a combined experimental and computational approach, we elucidate the role of the vanadium haloperoxidase protein in substrate halogenation. Activity studies show that binding of the substrate to the enzyme is essential for the reaction of the hypohalide with substrate. Stopped-flow measurements demonstrate that the rate-determining step is not dependent on substrate binding but partially on hypohalide formation. Using a combination of molecular mechanics (MM) and molecular dynamics (MD) simulations, the substrate binding area in the protein is identified and even though the selected substrates (methylphenylindole and 2-phenylindole) have limited hydrogen-bonding abilities, they are found to bind relatively strongly and remain stable in a binding tunnel. A subsequent analysis of the MD snapshots characterizes two small tunnels leading from the vanadate active site to the surface that could fit small molecules such as hypohalide, halide, and hydrogen peroxide. Density functional theory studies using electric field effects show that a polarized environment in a specific direction can substantially lower barriers for halogen transfer. A further analysis of the protein structure indeed shows a large dipole orientation in the substrate-binding pocket that could enable halogen transfer through an applied local electric field. These findings highlight the importance of the enzyme in catalyzing substrate halogenation by providing an optimal environment to lower the energy barrier for this challenging aromatic halide insertion reaction.

Electric-field induced entropic effects in liquid water

Externally applied electric fields in liquid water can induce a plethora of effects with wide implications in electrochemistry and hydrogen-based technologies. Although some effort has been made to elucidate the thermodynamics associated with the application of electric fields in aqueous systems, to the best of our knowledge, field-induced effects on the total and local entropy of bulk water have never been presented so far. Here, we report on classical TIP4P/2005 and ab initio molecular dynamics simulations measuring entropic contributions carried by diverse field intensities in liquid water at room temperature. We find that strong fields are capable of aligning large fractions of molecular dipoles. Nevertheless, the order-maker action of the field leads to quite modest entropy reductions in classical simulations. Albeit more significant variations are recorded during first-principles simulations, the associated entropy modifications are small compared to the entropy change involved in the freezing phenomenon, even at intense fields slightly beneath the molecular dissociation threshold. This finding further corroborates the idea that electrofreezing (i.e., the electric-field-induced crystallization) cannot take place in bulk water at room temperature. In addition, here, we propose a molecular-dynamics-based analysis (3D-2PT) that spatially resolves the local entropy and the number density of bulk water under an electric field, which enables us to map their field-induced changes in the environment of reference H2O molecules. By returning detailed spatial maps of the local order, the proposed approach is capable of establishing a link between entropic and structural modifications with atomistic resolution.

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.

Strong External Electric Fields Reduce Explosive Sensitivity: A Theoretical Investigation into the Reaction Selectivity in NH2NO2∙∙∙NH3

Controlling the selectivity of a detonation initiation reaction of explosive is essential to reduce sensitivity, and it seems impossible to reduce it by strengthening the external electric field. To verify this, the effects of external electric fields on the initiation reactions in NH2NO2∙∙∙NH3, a model system of the nitroamine explosive with alkaline additive, were investigated at the MP2/6-311++G(2d,p) and CCSD(T)/6-311++G(2d,p) levels. The concerted effect in the intermolecular hydrogen exchange is characterized by an index of the imaginary vibrations. Due to the weakened concerted effects by the electric field along the −x-direction opposite to the “reaction axis”, the dominant reaction changes from the intermolecular hydrogen exchange to 1,3-intramolecular hydrogen transference with the increase in the field strengths. Furthermore, the stronger the field strengths, the higher the barrier heights become, indicating the lower sensitivities. Therefore, by increasing the field strength and adjusting the orientation between the field and “reaction axis”, not only can the reaction selectivity be controlled, but the sensitivity can also be reduced, in particular under a super-strong field. Thus, a traditional concept, in which the explosive is dangerous under the super-strong external electric field, is theoretically broken. Compared to the neutral medium, a low sensitivity of the explosive with alkaline can be achieved under the stronger field. Employing atoms in molecules, reduced density gradient, and surface electrostatic potentials, the origin of the reaction selectivity and sensitivity change is revealed. This work provides a new idea for the technical improvement regarding adding the external electric field into the explosive system.

Theoretical Investigation of the X-ray Stark Effect in Small Molecules

We have studied the Stark effect in the soft x-ray region for various small molecules by calculating the field-dependent x-ray absorption spectra. This effect is explained in terms of the response of molecular orbitals (core and valence), the molecular dipole moment, and the molecular geometry to the applied electric field. A number of consistent trends are observed linking the computed shifts in absorption energies and intensities with specific features of the molecular electronic structure. We find that both the virtual molecular orbitals (valence and/or Rydberg) as well as the core orbitals contribute to observed trends in a complementary fashion. This initial study highlights the potential impact of x-ray Stark spectroscopy as a tool to study electronic structure and environmental perturbations at a submolecular scale.

How Can Static and Oscillating Electric Fields Serve in Decomposing Alzheimer's and Other Senile Plaques?

Alzheimer's disease is one of the most common neurodegenerative conditions, which are ascribed to extracellular accumulation of β-amyloid peptides into plaques. This phenomenon seems to typify other related neurodegenerative diseases. The present study uses classical molecular-dynamics simulations to decipher the aggregation-disintegration behavior of β-amyloid peptide plaques in the presence of static and oscillating oriented external electric fields (OEEFs). A long-term disintegration of such plaques is highly desirable since this may improve the prospects of therapeutic treatments of Alzheimer's disease and of other neurodegenerative diseases typified by senile plaques. Our study illustrates the spontaneous aggregation of the β-amyloid, its prevention and breakdown when OEEF is applied, and the fate of the broken aggregate when the OEEF is removed. Notably, we demonstrate that the usage of an oscillating OEEF on β-amyloid aggregates appears to lead to an irreversible disintegration. Insight is provided into the root causes of the various modes of aggregation, as well as into the different fates of OEEF-induced disintegration in oscillating vs static fields. Finally, our simulation results are compared to the well-established TTFields and the Deep Brain Stimulation (DBS) therapies, which are currently used options for treatments of Alzheimer's disease and other related neurodegenerative diseases.

Electric fields drive bond homolysis

Electric fields have been used to control and direct chemical reactions in biochemistry and enzymatic catalysis, yet directly applying external electric fields to activate reactions in bulk solution and to characterize them ex situ remains a challenge. Here we utilize the scanning tunneling microscope-based break-junction technique to investigate the electric field driven homolytic cleavage of the radical initiator 4-(methylthio)benzoic peroxyanhydride at ambient temperatures in bulk solution, without the use of co-initiators or photochemical activators. Through time-dependent ex situ quantification by high performance liquid chromatography using a UV-vis detector, we find that the electric field catalyzes the reaction. Importantly, we demonstrate that the reaction rate in a field increases linearly with the solvent dielectric constant. Using density functional theory calculations, we show that the applied electric field decreases the dissociation energy of the O-O bond and stabilizes the product relative to the reactant due to their different dipole moments.

Protonation state control of electric field induced molecular switching mechanisms

Scanning tunneling microscopy tip-induced deprotonation has been demonstrated experimentally and can be used as an additional control mechanism in electric-field induced molecular switching. The goal of the current work is to establish whether (de)protonation can be used to inhibit or enhance the electric field controlled thermal and photoisomerization processes. Dihydroxyazobenzene is used as a model system, where protonation/deprotonation of the free hydroxyl moiety changes the azo bond order, and so modifies the rate of electric field induced isomerization. Through the combined action of deprotonation and applied field, it was found that the cis-to-trans thermal isomerization barrier could be completely removed, changing the isomerization half-life from the order of several months. In addition, due to the presence of multiple isomerization mechanisms, electric fields could modify the isomerization kinetics by increasing the number of energetically viable isomerization pathways, rather than reducing the activation barrier of the lowest energy pathway. Excited state calculations indicated that the protonation state and electric field could be used together to control the presence of electronic degeneracies along the rotation pathway between S₀/S₁, and along all three pathways between S₁/S₂. This work provides insight into the mechanisms that enable the use of protonation state, light, and electric fields in concert to control molecular switches.

Does the Spin State and Oriented External Electric Field Boost the Efficiency of Fe(II) Pincer Catalyst toward CO2 Hydrogenation Reaction?

In this study, we have explored the catalytic reactivities of four PNP-pincer supported Fe(II) complexes, namely, [(^iPrPN^MeP)FeH₂(CO)] (1), [(^iPrPN^MeP)FeH(CO)(BH₄)] (2), [(^iPrPN^HP)FeH₂(CO)] (3), and [(^iPrPN^MeP)FeH(BH₄)] (4) (^iPrPN^MeP = MeN{CH₂CH₂(PⁱPr₂)}₂ and ^iPrPN^HP = HN{CH₂CH₂(PⁱPr₂)}₂) toward reductive CO₂ hydrogenation for formate production. Our density functional theory and ab initio complete active space self-consistent field study have identified three fundamental steps in this catalytic transformation: (i) anchoring of the CO₂ molecule in the vicinity of the metal using noncovalent interactions, (ii) catalyst regeneration via H₂ cleavage, and (iii) formate rebound step leading to catalytic poisoning. The variations in the catalytic efficiency observed among these catalysts were attributed to either easing of steps (i) and (ii) or the hampering step (iii). This can be achieved in various chemical/non-chemical ways, for instance, (a) incorporation of strong-field ligands such as CO facilitating single-state reactivity and eliminating two-state reactivity that generally enhances the rate and (b) inclusion of Lewis acids such as LiOTf and strong bases found to either avoid catalytic poisoning or ease the H-H cleavages, to enhance the rate of reaction (c) evading mixing of excited open-shell singlet states to the ground closed-shell singlet state that hampers the catalytic regeneration. We have probed the role of oriented external electric fields (OEEFs) in the entire mechanistic profile for the best and worst catalyst, and our study suggests that imposing OEEFs opposite to the reaction axis (z-axis) fastens the catalytic regeneration step and, at the same time, hampers catalytic poisoning. The application of OEEFs is found to regulate the energetics of various spin states and can hamper two-state reactivity, therefore increasing the efficiency. Thus, this study provides insights into the CO₂ hydrogenation mechanism where the role of bases/Lewis acid, ligand design, spin states, and electric field in a particular direction has been established and is, therefore, likely to pave the way forward for a new generation of catalysts.

A catastrophe theory-based model for optimal control of chemical reactions by means of oriented electric fields

The effect of oriented external electric fields (OEEF) on chemical reactivity has been studied theoretically and computationally in the last decades. A central goal in this research area is to predict the orientation and the smallest amplitude electric field that renders a barrierless chemical process with the smallest possible strength. Recently, a model to find the optimal electric field has been proposed and described (Bofill JM et al., J. Chem. Theory Comput. 18:935, 2022). We here proof that this model is based on catastrophe and optimum control theories. Based on both theories a technical treatment of the model is given and applied to a two-dimensional generic example that provides insight into its nature and capability. Finally, the model is applied to determine the optimal OEEF for the trans-to-cis isomerization of a [3]cumulene derivative.

Efficient Production of Epoxy-Norbornane from Norbornene by an Engineered P450 Peroxygenase

Epoxy-norbornane (EPO-NBE) is a crucial building block for the synthesis of various biologically active heterocyclic systems. To develop an efficient protocol for producing EPO-NBE using norbornene (NBE) as a substrate, cytochrome P450 enzyme from Pseudomonas putida (CYP238A1) was examined and its crystal structure (PDB code: 7X53) was resolved. Molecular mechanism analysis showed a high energy barrier related to iron-alkoxy radical complex formation. Therefore, a protein engineering strategy was developed and an optimal CYP238A1^NPV variant containing a local hydrophobic "fence" at the active site was obtained, which increased the H₂ O₂ -dependent epoxidation activity by 7.5-fold compared with that of CYP238A1^WT . Among the "fence", Glu255 participates in an efficient proton transfer system. Whole-cell transformation using CYP238A1^NPV achieved an EPO-NBE yield of 77.6 g ⋅ L^-1 in a 30-L reactor with 66.3 % conversion. These results demonstrate the potential of this system for industrial production of EPO-NBE and provides a new biocatalytic platform for epoxidation chemistry.

Biotransformation of Bisphenol by Human Cytochrome P450 2C9 Enzymes: A Density Functional Theory Study

Bisphenol A (BPA, 2,2-bis-(4-hydroxyphenyl)propane) is used as a precursor in the synthesis of polycarbonate and epoxy plastics; however, its availability in the environment is causing toxicity as an endocrine-disrupting chemical. Metabolism of BPA and their analogues (substitutes) is generally performed by liver cytochrome P450 enzymes and often leads to a mixture of products, and some of those are toxic. To understand the product distributions of P450 activation of BPA, we have performed a computational study into the mechanisms and reactivities using large model structures of a human P450 isozyme (P450 2C9) with BPA bound. Density functional theory (DFT) calculations on mechanisms of BPA activation by a P450 compound I model were investigated, leading to a number of possible products. The substrate-binding pocket is tight, and as a consequence, aliphatic hydroxylation is not feasible as the methyl substituents of BPA cannot reach compound I well due to constraints of the substrate-binding pocket. Instead, we find low-energy pathways that are initiated with phenol hydrogen atom abstraction followed by OH rebound to the phenolic ortho- or para-position. The barriers of para-rebound are well lower in energy than those for ortho-rebound, and consequently, our P450 2C9 model predicts dominant hydroxycumyl alcohol products. The reactions proceed through two-state reactivity on competing doublet and quartet spin state surfaces. The calculations show fast and efficient substrate activation on a doublet spin state surface with a rate-determining electrophilic addition step, while the quartet spin state surface has multiple high-energy barriers that can also lead to various side products including C⁴-aromatic hydroxylation. This work shows that product formation is more feasible on the low spin state, while the physicochemical properties of the substrate govern barrier heights of the rate-determining step of the reaction. Finally, the importance of the second-coordination sphere is highlighted that determines the product distributions and guides the bifurcation pathways.

Effects and Influence of External Electric Fields on the Equilibrium Properties of Tautomeric Molecules

In this review, we have attempted to briefly summarize the influence of an external electric field on an assembly of tautomeric molecules and to what experimentally observable effects this interaction can lead to. We have focused more extensively on the influence of an oriented external electric field (OEEF) on excited-state intramolecular proton transfer (ESIPT) from the studies available to date. The possibilities provided by OEEF for regulating several processes and studying physicochemical processes in tautomers have turned this direction into an attractive area of research due to its numerous applications.

Automated Variable Electric-Field DFT Application for Evaluation of Optimally Oriented Electric Fields on Chemical Reactivity

Recent theoretical work and experiments at molecular junctions have provided a strong conceptualization for the effects of oriented electric fields (OEFs) on organic reactions. Depending on the axis of application, OEFs can increase (or decrease) the reaction rate or distinguish between isomeric pathways. Despite the conceptual elegance of OEFs, which may be applied externally or induced locally, as tools for catalyzing organic reactions, implementation in synthetically relevant systems has been hampered by inefficiencies in evaluating reaction sensitivity to field effects. Herein, we describe the development of the Automated Variable Electric-Field DFT Application (A.V.E.D.A.) for streamlined evaluation of a reaction's susceptibility to OEFs. This open-source software was designed to be accessible for nonexpert users of computational and programming tools. Following initiation by a single command (and with no subsequent intervention) the Linux workflow manages a series of density functional theory calculations and mathematical manipulations to optimize local-minimum and transition-state structures in oriented electric fields of increasing magnitude. The resulting molecular and reaction dipole moments, field-perturbed geometries, and net effective activation energies are compiled for user interpretation. Ten representative pericyclic reactions that showcase the development and evaluation of A.V.E.D.A. are described.

The Performance of Different Water Models on the Structure and Function of Cytochrome P450 Enzymes

Modeling approaches and modern simulations to investigate the biomolecular structure and function rely on various methods. Since water molecules play a crucial role in all sorts of chemistry, the accurate modeling of water molecules is vital for such simulations. In cytochrome P450 (CYP450), in particular, water molecules play a key role in forming active oxidant that ultimately performs oxidation and metabolism. In the present study, we have highlighted the behavior of the three most widely used water models─TIP3P, SPC/E, and OPC─for three different CYP450 enzymes─CYP450_BM3, CYP450_OleT, and CYP450_BSβ─during MD simulations and QM/MM calculations. We studied the various properties, such as RMSD, RMSF, H-bond, water occupancy, and hydrogen atom transfer (HAT), using QM/MM calculations and compared them for all three water models. Our study shows that the stabilities of the enzyme complexes are well maintained in all three water models. However, the OPC water model performs well for the polar active sites, that is, in CYP450_OleT and CYP450_BSβ, while the TIP3P water model is superior for the hydrophobic site, such as CYP450_BM3.

β-Lactamases Evolve against Antibiotics by Acquiring Large Active-Site Electric Fields

A compound bound covalently to an enzyme active site can act either as a substrate if the covalent linkage is readily broken up by the enzyme or as an inhibitor if the bond dissociates slowly. We tracked the reactivity of such bonds associated with the rise of the resistance to penicillin G (PenG) in protein evolution from penicillin-binding proteins (PBPs) to TEM β-lactamases and with the development of avibactam (Avb) to overcome the resistance. We found that the ester linkage in PBP-PenG is resistant to hydrolysis mainly due to the small electric fields present in the protein active site. Conversely, the same linkage in the descendant TEM-PenG experiences large electric fields that stabilize the more charge-separated transition state and thus lower the free energy barrier to hydrolysis. Specifically, the electric fields were improved from -59 to -140 MV/cm in an ancient evolution dating back billions of years, contributing 5 orders of magnitude rate acceleration. This trend continues today in the nullification of newly developed antibiotic drugs. The fast linkage hydrolysis acquired from evolution is counteracted by the upgrade of PenG to Avb whose linkage escapes from the hydrolysis by returning to a low-field environment. Using the framework of electrostatic catalysis, the electric field, an observable from vibrational spectroscopy, provides a unifying physical metric to understand protein evolution and to guide the design of covalent drugs.

Single-Molecule Chemical Reactions Unveiled in Molecular Junctions

Understanding chemical processes at the single-molecule scale represents the ultimate limit of analytical chemistry. Single-molecule detection techniques allow one to reveal the detailed dynamics and kinetics of a chemical reaction with unprecedented accuracy. It has also enabled the discoveries of new reaction pathways or intermediates/transition states that are inaccessible in conventional ensemble experiments, which is critical to elucidating their intrinsic mechanisms. Thanks to the rapid development of single-molecule junction (SMJ) techniques, detecting chemical reactions via monitoring the electrical current through single molecules has received an increasing amount of attention and has witnessed tremendous advances in recent years. Research efforts in this direction have opened a new route for probing chemical and physical processes with single-molecule precision. This review presents detailed advancements in probing single-molecule chemical reactions using SMJ techniques. We specifically highlight recent progress in investigating electric-field-driven reactions, reaction dynamics and kinetics, host–guest interactions, and redox reactions of different molecular systems. Finally, we discuss the potential of single-molecule detection using SMJs across various future applications.

Geometric Structure, Electronic, and Spectral Properties of Metal-free Phthalocyanine under the External Electric Fields

Here, the ground-state structures, electronic structures, polarizability, and spectral properties of metal-free phthalocyanine (H₂Pc) under different external electric fields (EEFs) are investigated. The results show that EEF has an ultrastrong regulation effect on various aspects of H₂Pc; the geometric structures, electronic properties, polarizability, and spectral properties are strongly sensitive to the EEF. In particular, an EEF of 0.025 a.u. is an important control point: an EEF of 0.025 a.u. will bend the benzene ring subunits to the positive and negative x directions of the planar molecule. Flipping the EEF from positive (0.025 a.u.) to negative (-0.025 a.u.) flips also the bending direction of benzene ring subunits. The H₂Pc shows different dipole moments projecting an opposite direction along the x direction (-84 and 84 Debye for EEFs of -0.025 and 0.025 a.u., respectively) under negative and positive EEF, revealing a significant dipole moment transformation. Furthermore, when the EEF is removed, the molecule can be restored to the planar structure. The transformation of the H₂Pc structure can be induced by the EEF, which has potential applications in the molecular devices such as molecular switches or molecular forceps. EEF lowers total energy and reduces highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) gap; especially, an EEF of 0.025 a.u. can reduce the HOMO-LUMO gap from 2.1 eV (in the absence of EEF) to 0.37 eV, and thus, it can enhance the molecular conductivity. The first hyperpolarizability of H₂Pc is 0 in the absence of EEF; remarkably, an EEF of 0.025 a.u. can enhance the first hyperpolarizability up to 15,578 a.u. Therefore, H₂Pc under the EEF could be introduced as a promising innovative nonlinear optical (NLO) nanomaterial such as NLO switches. The strong EEF (0.025 a.u.) causes a large number of new absorption peaks in IR and Raman spectra and causes the redshift of electronic absorption spectra. The changes of EEF can be used to regulate the structure transformation and properties of H₂Pc, which can promote the application of H₂Pc in nanometer fields such as molecular devices.

Enzymatic Control over Reactive Intermediates Enables Direct Oxidation of Alkenes to Carbonyls by a P450 Iron-Oxo Species

The aerobic oxidation of alkenes to carbonyls is an important and challenging transformation in synthesis. Recently, a new P450-based enzyme (aMOx) has been evolved in the laboratory to directly oxidize styrenes to their corresponding aldehydes with high activity and selectivity. The enzyme utilizes a heme-based, high-valent iron-oxo species as a catalytic oxidant that normally epoxidizes alkenes, similar to other catalysts. How the evolved aMOx enzyme suppresses the commonly preferred epoxidation and catalyzes direct carbonyl formation is currently not well understood. Here, we combine computational modelling together with mechanistic experiments to study the reaction mechanism and unravel the molecular basis behind the selectivity achieved by aMOx. Our results describe that although both pathways are energetically accessible diverging from a common covalent radical intermediate, intrinsic dynamic effects determine the strong preference for epoxidation. We discovered that aMOx overrides these intrinsic preferences by controlling the accessible conformations of the covalent radical intermediate. This disfavors epoxidation and facilitates the formation of a carbocation intermediate that generates the aldehyde product through a fast 1,2-hydride migration. Electrostatic preorganization of the enzyme active site also contributes to the stabilization of the carbocation intermediate. Computations predicted that the hydride migration is stereoselective due to the enzymatic conformational control over the intermediate species. These predictions were corroborated by experiments using deuterated styrene substrates, which proved that the hydride migration is cis- and enantioselective. Our results demonstrate that directed evolution tailored a highly specific active site that imposes strong steric control over key fleeting biocatalytic intermediates, which is essential for accessing the carbonyl forming pathway and preventing competing epoxidation.

2,7- and 4,9-Dialkynyldihydropyrene Molecular Switches: Syntheses, Properties, and Charge Transport in Single-Molecule Junctions

This paper describes the syntheses of several functionalized dihydropyrene (DHP) molecular switches with different substitution patterns. Regioselective nucleophilic alkylation of a 5-substituted dimethyl isophthalate allowed the development of a workable synthetic protocol for the preparation of 2,7-alkyne-functionalized DHPs. Synthesis of DHPs with surface-anchoring groups in the 2,7- and 4,9-positions is described. The molecular structures of several intermediates and DHPs were elucidated by X-ray single-crystal diffraction. Molecular properties and switching capabilities of both types of DHPs were assessed by light irradiation experiments, spectroelectrochemistry, and cyclic voltammetry. Spectroelectrochemistry, in combination with density functional theory (DFT) calculations, shows reversible electrochemical switching from the DHP forms to the cyclophanediene (CPD) forms. Charge-transport behavior was assessed in single-molecule scanning tunneling microscope (STM) break junctions, combined with density functional theory-based quantum transport calculations. All DHPs with surface-contacting groups form stable molecular junctions. Experiments show that the molecular conductance depends on the substitution pattern of the DHP motif. The conductance was found to decrease with increasing applied bias.

Oriented External Electric Fields Regurating the Reaction Mechanism of CH4 Oxidation Catalyzed by Fe(IV)-Oxo-Corrolazine: Insight from Density Functional Calculations

Methane is the simplest alkane and can be used as an alternative energy source for oil and coal, but the greenhouse effect caused by its leakage into the air is not negligible, and its conversion into liquid methanol not only facilitates transportation, but also contributes to carbon neutrality. In order to find an efficient method for converting methane to methanol, CH4 oxidation catalyzed by Fe(IV)-Oxo-corrolazine (Fe(IV)-Oxo-Cz) and its reaction mechanism regulation by oriented external electric fields (OEEFs) are systematically studied by density functional calculations. The calculations show that Fe(IV)-Oxo-Cz can abstract one H atom from CH4 to form the intermediate with OH group connecting on the corrolazine ring, with the energy barrier of 25.44 kcal mol-1. And then the product methanol is formed through the following rebound reaction. Moreover, the energy barrier can be reduced to 20.72 kcal mol-1 through a two-state reaction pathway. Furthermore, the effect of OEEFs on the reaction is investigated. We found that OEEFs can effectively regulate the reaction by adjusting the stability of the reactant and the transition state through the interaction of electric field-molecular dipole moment. When the electric field is negative, the energy barrier of the reaction decreases with the increase of electric intensity. Moreover, the OEEF aligned along the intrinsic Fe‒O reaction axis can effectively regulate the ability of forming the OH on the corrolazine ring by adjusting the charges of O and H atoms. When the electric field intensity is -0.010 a.u., the OH can be directly rebounded to the CH3· before it is connecting on the corrolazine ring, thus forming the product directly from the transition state without passing through the intermediate with only an energy barrier of 17.34 kcal mol-1, which greatly improves the selectivity of the reaction.

Boosting the performance of single-atom catalysts via external electric field polarization

Single-atom catalysts represent a unique catalytic system with high atomic utilization and tunable reaction pathway. Despite current successes in their optimization and tailoring through structural and synthetic innovations, there is a lack of dynamic modulation approach for the single-atom catalysis. Inspired by the electrostatic interaction within specific natural enzymes, here we show the performance of model single-atom catalysts anchored on two-dimensional atomic crystals can be systematically and efficiently tuned by oriented external electric fields. Superior electrocatalytic performance have been achieved in single-atom catalysts under electrostatic modulations. Theoretical investigations suggest a universal "onsite electrostatic polarization" mechanism, in which electrostatic fields significantly polarize charge distributions at the single-atom sites and alter the kinetics of the rate determining steps, leading to boosted reaction performances. Such field-induced on-site polarization offers a unique strategy for simulating the catalytic processes in natural enzyme systems with quantitative, precise and dynamic external electric fields.

Double proton transfer in hydrated formic acid dimer: Interplay of spatial symmetry and solvent-generated force on reactivity

The double proton transfer (DPT) reaction in the hydrated formic acid dimer (FAD) is investigated at molecular-level detail. For this, a global and reactive machine learned (ML) potential energy surface (PES) is developed to run extensive (more than 100 ns) mixed ML/MM molecular dynamics (MD) simulations in explicit molecular mechanics (MM) solvent at MP2-quality for the solute. Simulations with fixed – as in a conventional empirical force field – and conformationally fluctuating – as available from the ML-based PES – charge models for FAD show a significant impact on the competition between DPT and dissociation of FAD into two formic acid monomers. With increasing temperature the barrier height for DPT in solution changes by about 10% (∼1 kcal mol⁻¹) between 300 K and 600 K. The rate for DPT is largest, ∼1 ns⁻¹, at 350 K and decreases for higher temperatures due to destabilisation and increased probability for dissociation of FAD. The water solvent is found to promote the first proton transfer by exerting a favourable solvent-induced Coulomb force along the O–H⋯O hydrogen bond whereas the second proton transfer is significantly controlled by the O–O separation and other conformational degrees of freedom. Double proton transfer in hydrated FAD is found to involve a subtle interplay and balance between structural and electrostatic factors.

Simulation of double proton transfer in formic acid dimer by reactive ML potential in explicit molecular mechanics water solvent.

Influence of an Oriented External Electric Field on the Mechanism of Double Proton Transfer between Pyrazole and Guanidine: from an Asynchronous Plateau Transition State to a Synchronous or Stepwise Mechanism

The double proton transfer (DPT) reaction between pyrazole and guanidine, a concerted reaction but strongly asynchronous and presenting a "plateau transition region", has been theoretically reinvestigated in the presence of an external uniform electric field. First, we computed the reaction path by DFT and proposed a very detailed description of the constitutive electronic events, based on the ELF topology and the bond evolution theory. Then, we studied the effect of an oriented external electric field (OEEF) on the reaction mechanism, for an OEEF oriented along the proton transfer axis. We observe that in one direction, the DPT reaction can be transformed into a stepwise reaction, going through a stabilized single proton transferred intermediate. Contrarily, the two proton transfers occur simultaneously when the electric field is applied in the opposite direction. In the latter case, the order in which the two protons are transferred in the same elementary step can even be reversed if the OEEF is intense enough. Finally, it has been shown that the evolution of the double proton transfer reaction in the presence of an electric field could be quantitatively anticipated by analyzing the ELF value at the bifurcation point between V(A, H) proton donor and V(B) proton acceptor of the double hydrogen bonded complex in the entrance channel.