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.

Dinitrogen Activation within Frustrated Lewis Pairs Is Promoted by Adding External Electric Fields.

This study uses computational means to explore the feasibility of N2 cleavage by frustrated Lewis pair (FLPs) species. The employed FLP systems are phosphane/borane (1) and carbene/borane (2). Previous studies show that 1 and 2 react with H2 and CO2 but do not activate N2. The present study demonstrates that N2 is indeed inert, and its activation requires augmentation of the FLPs by an external tool. As we demonstrate here, FLP-mediated N2 activation can be achieved by an external electric field oriented along the reaction axis of the FLP. Additionally, the study demonstrates that FLP -N2 activation generates useful nitrogen compound, e.g., hydrazine (H2N-NH2). In summary, we conclude that FLP effectively activates N2 in tandem with oriented external electric fields (OEEFs), which play a crucial role. This FLP/OEEF combination may serve as a general activator of inert molecules.

On the engineering of reductase-based-monooxygenase activity in CYP450 peroxygenases.

Recent bioengineering of CYP450OleT shows that peroxide-based CYP450OleT can be converted to a reductase-based self-sufficient enzyme, which is capable of showing efficient hydroxylation and decarboxylation activity for a wide range of substrates. The so-generated enzyme creates several mechanistic puzzles: (A) as CYP450 peroxygenases lack the conventional acid-alcohol pair, what is the source of two protons that are required to create the ultimate oxidant Cpd I? (B) Why is it only CYP450OleT that shows the reductase-based activity but no other CYP members? The present study provides a mechanistic solution to these puzzles using comprehensive MD simulations and hybrid QM/MM calculations. We show that the fusion of the reductase domain to the heme-binding domain triggers significant conformational rearrangement, which is gated by the propionate side chain, which constitutes a new water aqueduct via the carboxylate end of the substrate that ultimately participates in Cpd I formation. Importantly, such well-synchronized choreographies are controlled by remotely located Tyr359, which senses the fusion of reductase and communicates to the heme domain via non-covalent interactions. These findings provide crucial insights and a broader perspective which enables us to make a verifiable prediction: thus, the catalytic activity is not only limited to the first or second catalytic shell of an enzyme. Furthermore, it is predicted that reinstatement of tyrosine at a similar position in other members of CYP450 peroxygenases can convert these enzymes to reductase-based monooxygenases.

Flavin-N5OOH Functions as both a Powerful Nucleophile and a Base in the Superfamily of Flavoenzymes.

Flavoenzymes can mediate a large variety of oxidation reactions through the activation of oxygen. However, the O2 activation chemistry of flavin enzymes is not yet fully exploited. Normally, the O2 activation occurs at the C4a site of the flavin cofactor, yielding the flavin C4a-(hydro)hydroperoxyl species in monooxygenases or oxidases. Using extensive MD simulations, QM/MM calculations and QM calculations, our studies reveal the formation of the common nucleophilic species, Flavin-N5OOH, in two distinct flavoenzymes (RutA and EncM). Our studies show that Flavin-N5OOH acts as a powerful nucleophile that promotes C-N cleavage of uracil in RutA, and a powerful base in the deprotonation of substrates in EncM. We reason that Flavin-N5OOH can be a common reactive species in the superfamily of flavoenzymes, which accomplish generally selective general base catalysis and C-X (X=N, S, Cl, O) cleavage reactions that are otherwise challenging with solvated hydroxide ion base. These results expand our understanding of the chemistry and catalysis of flavoenzymes.

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.

Valence Bond Theory Allows a Generalized Description of Hydrogen Bonding.

This paper describes the nature of the hydrogen bond (HB), B:---H-A, using valence bond theory (VBT). Our analysis shows that the most important HB interactions are polarization and charge transfer, and their corresponding sum displays a pattern that is identical for a variety of energy decomposition analysis (EDA) methods. Furthermore, the sum terms obtained with the different EDA methods correlate linearly with the corresponding VB quantities. The VBT analysis demonstrates that the total covalent-ionic resonance energy (RECS) of the HB portion (B---H in B:---H-A) correlates linearly with the dissociation energy of the HB, ΔEdiss. In principle, therefore, RECS(HB) can be determined by experiment. The VBT wavefunction reveals that the contributions of ionic structures to the HB increase the positive charge on the hydrogen of the corresponding external/free O-H bonds in, for example, the water dimer compared with a free water molecule. This increases the electric field of the external O-H bonds of water clusters and contributes to bringing about catalysis of reactions by water droplets and in water-hydrophobic interfaces.

A porphyrin-based molecular cage guided by designed local-electric field is highly selective and efficient.

The present work outlines a general methodology for designing efficient catalytic machineries that can easily be tweaked to meet the demands of the target reactions. This work utilizes a principle of the designed local electric field (LEF) as the driver for an efficient catalyst. It is demonstrated that by tweaking the LEF, we can catalyze the desired hydroxylation products with enantioselectivity that can be changed at will. Using computation tools, we caged a synthetic analog of heme porphyrin (HM1) and investigated the pharmaceutically relevant conversion of tetralin to tetralol, inside the modified supramolecular cage. The QM/MM calculations demonstrate a resulting catalytic efficiency with virtually absolute R-selectivity for the tetralin hydroxylation. Our calculations show that the LEF of the supramolecular cage and HM1 exert a strong electric field along the Fe-O reaction axis, which is the main driving force for enhanced reactivity. At the same time, the supramolecular cage applies a lateral LEF that regulates the enantioselectivity. We further demonstrate that swapping the charged/polar substitution in the supramolecular cage switches the lateral LEF which changes the enantioselectivity of hydroxylation from R to S.

The Reactivity-Enhancing Role of Water Clusters in Ammonia Aqueous Solutions.

Among the many prototypical acid-base systems, ammonia aqueous solutions hold a privileged place, owing to their omnipresence in various planets and their universal solvent character. Although the theoretical optimal water-ammonia molar ratio to form NH4+ and OH- ion pairs is 50:50, our ab initio molecular dynamics simulations show that the tendency of forming these ionic species is inversely (directly) proportional to the amount of ammonia (water) in ammonia aqueous solutions, up to a water-ammonia molar ratio of ∼75:25. Here we prove that the reactivity of these liquid mixtures is rooted in peculiar microscopic patterns emerging at the H-bonding scale, where the highly orchestrated motion of 5 solvating molecules modulates proton transfer events through local electric fields. This study demonstrates that the reaction of water with NH3 is catalyzed by a small cluster of water molecules, in which an H atom possesses a high local electric field, much like the effect observed in catalysis by water droplets [ PNAS 2023, 120, e2301206120].

Two-Way Catalysis in a Diels-Alder Reaction Limits Inhibition Induced by an External Electric Field.

Oriented external electric fields (EEFs) act as catalysts that can induce selectivity in chemical reactions. The responses of the Diels-Alder (DA) reaction between butadiene and ethylene (BDE-DA) as well as cyclopentadiene and ethylene (CPDE-DA) towards EEF stimuli are investigated here using density functional theory (B3LYP) calculations. EEF is a vector that catalyzes the reaction in one direction while inhibiting it in the opposite direction. Here we report that the inhibitive direction becomes rate-enhancing after some increase in the EEF. The EEF value that brings about the maximum possible inhibition for the reaction is defined as the electrostatic resistance point (ERP). The possibility of both normal and inverse electron-demand DA reactions causes catalytic activity in both directions of the EEF starting at a unique ERP value. The C5 substituents of cyclopentadiene control the ERP values depending upon the resistance power that the functional group provides against the EEF. The endo and exo diastereomeric transition states of the DA reaction have distinct ERP values and the difference (ΔERP) provides the through-space electrostatic contribution to the stereoselectivity on a relative scale. Thus, the ERP values can be used as a gauge for the electrostatic interactions between substituent groups and external stimuli.

Electro-Freezing of Supercooled Water Is Induced by Hydrated Al3+ and Mg2+ Ions: Experimental and Theoretical Studies.

This work reports that the octahedral hydrated Al3+ and Mg2+ ions operate within electrolytic cells as kosmotropic (long-range order-making) "ice makers" of supercooled water (SCW). 10-5 M solutions of hydrated Al3+ and Mg2+ ions each trigger, near the cathode (-20 ± 5 V), electro-freezing of SCW at -4 °C. The hydrated Al3+ ions do so with 100% efficiency, whereas the Mg2+ ions induce icing with 40% efficiency. In contrast, hydrated Na+ ions, under the same experimental conditions, do not induce icing differently than pure water. As such, our study shows that the role played by Al3+ and Mg2+ ions in water electro-freezing is impacted by two synchronous effects: (1) a geometric effect due to the octahedral packing of the coordinated water molecules around the metallic ions, and (2) the degree of polarization which these two ions induce and thereby acidify the coordinated water molecules, which in turn imparts them with an ice-like structure. Long-duration molecular dynamics (MD) simulations of the Al3+ and Mg2+ indeed reveal the formation of "ice-like" hexagons in the vicinity of these ions. Furthermore, the MD shows that these hexagons and the electric fields of the coordinate water molecules give rise to ultimate icing. As such, the MD simulations provide a rational explanation for the order-making properties of these ions during electro-freezing.

Voltage-driven control of single-molecule keto-enol equilibrium in a two-terminal junction system.

Keto-enol tautomerism, describing an equilibrium involving two tautomers with distinctive structures, provides a promising platform for modulating nanoscale charge transport. However, such equilibria are generally dominated by the keto form, while a high isomerization barrier limits the transformation to the enol form, suggesting a considerable challenge to control the tautomerism. Here, we achieve single-molecule control of a keto-enol equilibrium at room temperature by using a strategy that combines redox control and electric field modulation. Based on the control of charge injection in the single-molecule junction, we could access charged potential energy surfaces with opposite thermodynamic driving forces, i.e., exhibiting a preference for the conducting enol form, while the isomerization barrier is also significantly reduced. Thus, we could selectively obtain desired and stable tautomers, which leads to significant modulation of the single-molecule conductance. This work highlights the concept of single-molecule control of chemical reactions on more than one potential energy surface.

Oriented External Electric-Field Effects on the Activation of Aryl CO Bond in Anisole using Rh(PEP) (E=Al, B, Ga) Catalysts.

The present work reports a DFT-based mechanistic investigation of aryl C-O bond activation in anisole catalysed by a Rh-Al pincer-type complex at room temperature. The study is extended to analogues Rh-E complexes based on Group 13 elements (E=B/Ga). Our results show a preference for the heterolytic cleavage pathway over oxidative addition in the C-O bond activation. The calculated barriers range from 16 to 36 kcal/mol, following the order: E=Al

Dichotomy of Delocalization/Localization and Charge-Shift Bonding in Germanazene and its Heavier Group 14 Analogues: a Valence Bond Study.

We present here a valence bond analysis of structure and π-delocalization in Ge3 (NH)3 , which models germanazene that was prepared by Power et al. To get a broader perspective, we explore the entire E3 (NH)3 series (E=C, Si, Ge, Sn, Pb). Thus, while (4n+2)π systems of carbon rings are aromatic with cyclic π-delocalization, the E3 (NH)3 rings are dominated by a nonbonded structure, wherein π-lone pairs are localized on the N atoms. Nevertheless, these molecules enjoy large covalent-ionic resonance energies of 153.0, 86.6, 74.2, 61.2, and 58.9 kcal/mol, respectively, for E=C, Si, Ge, Sn, Pb. The covalent-ionic mixing in E3 (NH)3 creates π-systems, which are stabilized by charge-shift bonding. Thus, unlike in benzene, in Ge3 (NH)3 delocalization of π-electron pairs of the N atoms is primarily confined to the domains of their adjacent Ge atoms. These features carry over to the substituted germanazene, Ge3 (NAr)3 (Ar=Ph).

How the Conformational Movement of the Substrate Drives the Regioselective C-N Bond Formation in P450 TleB: Insights from Molecular Dynamics Simulations and Quantum Mechanical/Molecular Mechanical Calculations.

P450 TleB catalyzes the oxidative cyclization of the dipeptide N-methylvalyl-tryptophanol into indolactam V through selective intramolecular C-H bond amination at the indole C4 position. Understanding its catalytic mechanism is instrumental for the engineering or design of P450-catalyzed C-H amination reactions. Using multiscale computational methods, we show that the reaction proceeds through a diradical pathway, involving a hydrogen atom transfer (HAT) from N1-H to Cpd I, a conformational transformation of the substrate radical species, and a second HAT from N13-H to Cpd II. Intriguingly, the conformational transformation is found to be the key to enabling efficient and selective C-N coupling between N13 and C4 in the subsequent diradical coupling reaction. The underlined conformational transformation is triggered by the first HAT, which proceeds with an energy-demanding indole ring flip and is followed by the facile approach of the N13-H group to Cpd II. Detailed analysis shows that the internal electric field (IEF) from the protein environment plays key roles in the transformation process, which not only provides the driving force but also stabilizes the flipped conformation of the indole radical. Our simulations provide a clear picture of how the P450 enzyme can smartly modulate the selective C-N coupling reaction. The present findings are in line with all available experimental data, highlighting the crucial role of substrate dynamics in controlling this highly valuable reaction.

Decarboxylation and Protonation Enigma in the H85Q Mutant of Cytochrome P450OleT.

Cytochrome P450OleT (CYP450OleT), a member of CYP450 peroxygenases, catalyzes unusual decarboxylation activity. Unlike other members of the peroxygenases family, CYP450OleT possesses a histidine at the 85th position, which was supposed to be the root cause of the decarboxylation activity in CYP450OleT. This work addresses the His85 → Gln mutant paradox, where mutation of His → Gln still shows efficient decarboxylation activity in CYP450OleT. The MD simulation of the H85Q mutant of CYP450OleT shows that in the absence of the histidine at the 85th position, an Asp239 plays a similar role via a well-organized water channel. Our simulation shows that such a water channel is vital for the optimal substrate positioning needed for the decarboxylation activity and is gated by the Q85-N242 residue pair. Interestingly, the MD simulation of the WT CYP450BSß shows a closed channel that blocks access to the Glu236 (analogous residue to Asp239 in CYP450OleT), and therefore, CYP450BSß shows low decarboxylation activity.

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.

Mechanistic Insight into Peptidyl-Cysteine Oxidation by the Copper-Dependent Formylglycine-Generating Enzyme.

The copper-dependent formylglycine-generating enzyme (FGE) catalyzes the oxygen-dependent oxidation of specific peptidyl-cysteine residues to formylglycine. Our QM/MM calculations provide a very likely mechanism for this transformation. The reaction starts with dioxygen binding to the tris-thiolate CuI center to form a triplet CuII -superoxide complex. The rate-determining hydrogen atom abstraction involves a triplet-singlet crossing to form a CuII -OOH species that couples with the substrate radical, leading to a CuI -alkylperoxo intermediate. This is accompanied by proton transfer from the hydroperoxide to the S atom of the substrate via a nearby water molecule. The subsequent O-O bond cleavage is coupled with the C-S bond breaking that generates the formylglycine and a CuII -oxyl complex. Moreover, our results suggest that the aldehyde oxygen of the final product originates from O2 , which will be useful for future experimental work.

Local Electric Fields: From Enzyme Catalysis to Synthetic Catalyst Design.

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

How Do Preorganized Electric Fields Function in Catalytic Cycles? The Case of the Enzyme Tyrosine Hydroxylase.

Nature has devised intrinsic electric fields (IEFs) that are engaged in electrostatic catalysis of enzymes. But, how does the IEF target its function in enzymes that involve several reaction steps in catalytic cycles? To decipher the impact of the IEF on the catalytic cycle of an enzyme system, we have performed molecular dynamics and quantum-mechanical/molecular-mechanical (QM/MM) simulations on tyrosine hydroxylase (TyrH). The catalytic cycle of TyrH involves two reaction stages: the activation of H2O2 to form the active species of compound I (Cpd I), in the first stage, and the Cpd I-mediated hydroxylation of l-tyrosine to l-DOPA, in the second stage. For the first stage, the QM/MM calculations show that a heme-propionate group functions as a base to catalyze the O-O heterolysis reaction. For the second stage, the study reveals that the reaction is initiated by the His88-mediated proton-coupled electron transfer followed by the oxygen atom transfer from compound II (Cpd II) to the l-Tyr substrate. Importantly, our calculations demonstrate that the IEF in TyrH is optimized to promote the O-O bond heterolysis that generates the active species of the enzyme, Cpd I. However, the same IEF slows down the subsequent aromatic hydroxylation. Thus, the IEF in the TyrH enzymes does not catalyze the product formation step, but will selectively boost one or more challenging steps in the catalytic cycle. These findings have general implications on O2/H2O2-dependent metalloenzymes, which can expand our understanding of how nature has used electric fields as "smart reagents" in modulating the catalytic reactivity.