Extreme electric fields power catalysis in the active site of ketosteroid isomerase

Evaluating aliphatic CF, CF2, and CF3 groups as vibrational Stark effect reporters

Given the extensive use of fluorination in molecular design, it is imperative to understand the solvation properties of fluorinated compounds and the impact of the C-F bond on electrostatic interactions. Vibrational spectroscopy can provide direct insights into these interactions by using the C-F bond stretching [v(C-F)] as an electric field probe through the vibrational Stark effect (VSE). In this work, we explore the VSE of the three basic patterns of aliphatic fluorination, i.e., mono-, di-, and trifluorination in CF, CF2, and CF3 groups, respectively, and compare their response to the well-studied aromatic v(C-F). Magnitudes (i.e., Stark tuning rates) and orientations of the difference dipole vectors of the v(C-F)-containing normal modes were determined using density functional theory and a molecular dynamics (MD)-assisted solvatochromic analysis of model compounds in solvents of varying polarity. We obtain Stark tuning rates of 0.2-0.8 cm-1/(MV/cm), with smallest and largest electric field sensitivities for CFaliphatic and CF3,aliphatic, respectively. While average electric fields of solvation were oriented along the main symmetry axis of the CFn, and thus along its static dipole, the Stark tuning rate vectors were tilted by up to 87° potentially enabling to map electrostatics in multiple dimensions. We discuss the influence of conformational heterogeneity on spectral shifts and point out the importance of multipolar and/or polarizable MD force fields to describe the electrostatics of fluorinated molecules. The implications of this work are of direct relevance for studies of fluorinated molecules as found in pharmaceuticals, fluorinated peptides, and proteins.

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.

Practical Synthesis of Chiral α-Aminophosphonates with Weak Bonding Organocatalysis at ppm Loading

α-Aminophosphonic acids as an important class of bioisosteres of α-amino acids demonstrate various biologically important activities. We report here the development of a highly enantioselective isomerization of α-iminophosphonates enabled by an extraordinarily efficient organocatalyst. This organocatalyst afforded a total turnover number (TON) of 20,000-1,000,000 for a wide range of α-alkyl iminophosphonates. Even at a parts-per-million (ppm) loading, this catalyst achieved a complete reaction in greater than 93% enantiomeric excess (ee). Computational studies revealed that this small-molecule catalyst achieved enzyme-like efficiency via a network of weak bonding interactions that effectively preorganized the substrate and catalyst toward a transition-state-like complex. Considering the substrate tolerance, catalytic efficiency, and mechanism, this organocatalyst could be regarded as a small-molecule isomerase.

Long-range electrostatic effects from intramolecular Lewis acid binding influence the redox properties of cobalt-porphyrin complexes

A CoII-porphyrin complex (1) with an appended aza-crown ether for Lewis acid (LA) binding was synthesized and characterized. NMR spectroscopy and electrochemistry show that cationic group I and II LAs (i.e., Li+, Na+, K+, Ca2+, Sr2+, and Ba2+) bind to the aza-crown ether group of 1. The binding constant for Li+ is comparable to that observed for a free aza-crown ether. LA binding causes an anodic shift in the CoII/CoI couple of between 10 and 40 mV and also impacts the CoIII/CoII couple. The magnitude of the anodic shift of the CoII/CoI couple varies linearly with the strength of the LA as determined by the pKa of the corresponding metal-aqua complex, with dications giving larger shifts than monocations. The extent of the anodic shift of the CoII/CoI couple also increases as the ionic strength of the solution decreases. This is consistent with electric field effects being responsible for the changes in the redox properties of 1 upon LA binding and provides a novel method to tune the reduction potential. Density functional theory calculations indicate that the bound LA is 5.6 to 6.8 Å away from the CoII ion, demonstrating that long-range electrostatic effects, which do not involve changes to the primary coordination sphere, are responsible for the variations in redox chemistry. Compound 1 was investigated as a CO2 reduction electrocatalyst and shows high activity but rapid decomposition.

Perturbation of Fermi Resonance on Hydrogen-Bonded > C═O: 2D IR Studies of Small Ester Probes

We utilized linear and 2D infrared spectroscopy to analyze the carbonyl stretching modes of small esters in different solvents. Particularly noteworthy were the distinct carbonyl spectral line shapes in aqueous solutions, prompting our investigation of the underlying factors responsible for these differences. Through our experimental and theoretical calculations, we identified the presence of the hydrogen-bond-induced Fermi resonance as the primary contributor to the varied line shapes of small esters in aqueous solutions. Furthermore, our findings revealed that the skeletal deformation mode plays a crucial role in the Fermi resonance for all small esters. Specifically, the first overtone band of the skeletal deformation mode intensifies when hydrogen bonds form with the carbonyl group of esters, whereas such coupling is rare in aprotic organic solvents. These spectral insights carry significant implications for the utilization of esters as infrared probes in both biological and chemical systems.

The impact of solvation on the structure and electric field strength in Li+GlyGly complexes

To scrutinise the impact of electric fields on the structure and vibrations of biomolecules in the presence of water, we study the sequential solvation of lithium diglycine up to three water molecules with cryogenic infrared action spectroscopy. Conformer-specific IR-IR spectroscopy and H2O/D2O isotopic substitution experiments provide most of the information required to decipher the structure of the observed conformers. Additional confirmation is provided by scaled harmonic vibrational frequency calculations using MP2 and DFT. The first water molecule is shown to bind to the Li+ ion, which weakens the electric field experienced by the peptide and as a consequence, also the strength of an internal NH⋯NH2 hydrogen bond in the diglycine backbone. The strength of this hydrogen bond decreases approximately linearly with the number of water molecules as a result of the decreasing electric field strength and coincides with an increase in the number of conformers observed in our spectra. The addition of two water molecules is already sufficient to change the preferred conformation of the peptide backbone, allowing for Li+ coordination to the lone pair of the terminal amine group.

Electric-field-assisted proton coupling enhanced oxygen evolution reaction

The discovery of Mn-Ca complex in photosystem II stimulates research of manganese-based catalysts for oxygen evolution reaction (OER). However, conventional chemical strategies face challenges in regulating the four electron-proton processes of OER. Herein, we investigate alpha-manganese dioxide (α-MnO2) with typical MnIV-O-MnIII-HxO motifs as a model for adjusting proton coupling. We reveal that pre-equilibrium proton-coupled redox transition provides an adjustable energy profile for OER, paving the way for in-situ enhancing proton coupling through a new "reagent"- external electric field. Based on the α-MnO2 single-nanowire device, gate voltage induces a 4-fold increase in OER current density at 1.7 V versus reversible hydrogen electrode. Moreover, the proof-of-principle external electric field-assisted flow cell for water splitting demonstrates a 34% increase in current density and a 44.7 mW/cm² increase in net output power. These findings indicate an in-depth understanding of the role of proton-incorporated redox transition and develop practical approach for high-efficiency electrocatalysis.

Deciphering the Microdroplet Acceleration Factors of Aza-Michael Addition Reactions

Microdroplet chemistry is emerging as a great tool for accelerating reactions by several orders of magnitude. Several unique properties such as extreme pHs, interfacial electric fields (IEFs), and partial solvation have been reported to be responsible for the acceleration; however, which factor plays the key role remains elusive. Here, we performed quantum chemical calculations to explore the underlying mechanisms of an aza-Michael addition reaction between methylamine and acrylamide. We showed that the acceleration in methanol microdroplets results from the cumulative effects of several factors. The acidic surface of the microdroplet plays a dominating role, leading to a decrease of ∼9 kcal/mol in the activation barrier. We speculated that the dissociation of both methanol and trace water contributes to the surface acidity. An IEF of 0.1 V/Å can further decrease the barrier by ∼2 kcal/mol. Partial solvation has a negligible effect on lowering the activation barrier in microdroplets but can increase the collision frequency between reactants. With acidity revealed to be the major accelerating factor for methanol droplets, reactions on water microdroplets should have even higher rates because water is more acidic. Both theoretically and experimentally, we confirmed that water microdroplets significantly accelerate the aza-Michael reaction, achieving an acceleration factor that exceeds 107. This work elucidates the multifactorial influences on the microdroplet acceleration mechanism, and with such detailed mechanistic investigations, we anticipate that microdroplet chemistry will be an avenue rich in opportunities in the realm of green synthesis.

The Scope of the Applicability of Non-relativistic DFT Calculations of NMR Chemical Shifts in Pyridine-Metal Complexes for Applied Applications

Heavy metals are toxic, but it is impossible to stop using them. Considering the variety of molecular systems in which they can be present, the multicomponent nature and disorder of the structure of such systems, one of the most effective methods for studying them is NMR spectroscopy. This determines the need to calculate NMR chemical shifts for expected model systems. For elements beyond the third row of the periodic table, corrections for relativistic effects are necessary when calculating NMR parameters. Such corrections may be necessary even for light atoms due to the shielding effect of a neighboring heavy atom. This work examines the extent to which non-relativistic DFT calculations are able to reproduce experimental 15N and 113Cd NMR chemical shift tensors in pyridine-metal coordination complexes. It is shown that while for the calculation of 15N NMR chemical shift tensors there is no real need to consider relativistic corrections, for 113Cd, on the contrary, none of the tested calculation methods could reproduce the experimentally obtained tensor to any extent correctly.

Potential Dependent Reorientation Controlling Activity of a Molecular Electrocatalyst

The activity of molecular electrocatalysts depends on the interplay of electrolyte composition near the electrode surface, the composition and morphology of the electrode surface, and the electric field at the electrode-electrolyte interface. This interplay is challenging to study and often overlooked when assessing molecular catalyst activity. Here, we use surface specific vibrational sum frequency generation (VSFG) spectroscopy to study the solvent and potential dependent activation of Mo(bpy)(CO)4, a CO2 reduction catalyst, at a polycrystalline Au electrode. We find that the parent complex undergoes potential dependent reorientation at the electrode surface when a small amount of N-methyl-2-pyrrolidone (NMP) is present. This preactivates the complex, resulting in greater yields at less negative potentials, of the active electrocatalyst for CO2 reduction.

Critical Evaluation of Polarizable and Nonpolarizable Force Fields for Proteins Using Experimentally Derived Nitrile Electric Fields

Molecular dynamics (MD) simulations are frequently carried out for proteins to investigate the role of electrostatics in their biological function. The choice of force field (FF) can significantly alter the MD results, as the simulated local electrostatic interactions lack benchmarking in the absence of appropriate experimental methods. We recently reported that the transition dipole moment (TDM) of the popular nitrile vibrational probe varies linearly with the environmental electric field, overcoming well-known hydrogen bonding (H-bonding) issues for the nitrile frequency and, thus, enabling the unambiguous measurement of electric fields in proteins (J. Am. Chem. Soc. 2022, 144 (17), 7562-7567). Herein, we utilize this new strategy to enable comparisons of experimental and simulated electric fields in protein environments. Specifically, previously determined TDM electric fields exerted onto nitrile-containing o-cyanophenylalanine residues in photoactive yellow protein are compared with MD electric fields from the fixed-charge AMBER FF and the polarizable AMOEBA FF. We observe that the electric field distributions for H-bonding nitriles are substantially affected by the choice of FF. As such, AMBER underestimates electric fields for nitriles experiencing moderate field strengths; in contrast, AMOEBA robustly recapitulates the TDM electric fields. The FF dependence of the electric fields can be partly explained by the presence of additional negative charge density along the nitrile bond axis in AMOEBA, which is due to the inclusion of higher-order multipole parameters; this, in turn, begets more head-on nitrile H-bonds. We conclude by discussing the implications of the FF dependence for the simulation of nitriles and proteins in general.

Electrically driven proton transfer promotes Brønsted acid catalysis by orders of magnitude

Electric fields play a key role in enzymatic catalysis and can enhance reaction rates by 100,000-fold, but the same rate enhancements have yet to be achieved in thermochemical heterogeneous catalysis. In this work, we probe the influence of catalyst potential and interfacial electric fields on heterogeneous Brønsted acid catalysis. We observed that variations in applied potential of ~380 mV led to a 100,000-fold rate enhancement for 1-methylcyclopentanol dehydration, which was catalyzed by carbon-supported phosphotungstic acid. Mechanistic studies support a model in which the interfacial electrostatic potential drop drives quasi-equilibrated proton transfer to the adsorbed substrate prior to rate-limiting C-O bond cleavage. Large increases in rate with potential were also observed for the same reaction catalyzed by Ti/TiOyHx and for the Friedel Crafts acylation of anisole with acetic anhydride by carbon-supported phosphotungstic acid.

Theoretical Improvements in Enzyme Efficiency Associated with Noisy Rate Constants and Increased Dissipation

Previous studies have revealed the extraordinarily large catalytic efficiency of some enzymes. High catalytic proficiency is an essential accomplishment of biological evolution. Natural selection led to the increased turnover number, kcat, and enzyme efficiency, kcat/KM, of uni-uni enzymes, which convert a single substrate into a single product. We added or multiplied random noise with chosen rate constants to explore the correlation between dissipation and catalytic efficiency for ten enzymes: beta-galactosidase, glucose isomerase, β-lactamases from three bacterial strains, ketosteroid isomerase, triosephosphate isomerase, and carbonic anhydrase I, II, and T200H. Our results highlight the role of biological evolution in accelerating thermodynamic evolution. The catalytic performance of these enzymes is proportional to overall entropy production-the main parameter from irreversible thermodynamics. That parameter is also proportional to the evolutionary distance of β-lactamases PC1, RTEM, and Lac-1 when natural or artificial evolution produces the optimal or maximal possible catalytic efficiency. De novo enzyme design and attempts to speed up the rate-limiting catalytic steps may profit from the described connection between kinetics and thermodynamics.

Molecular mechanisms and evolutionary robustness of a color switch in proteorhodopsins

Proteorhodopsins are widely distributed photoreceptors from marine bacteria. Their discovery revealed a high degree of evolutionary adaptation to ambient light, resulting in blue- and green-absorbing variants that correlate with a conserved glutamine/leucine at position 105. On the basis of an integrated approach combining sensitivity-enhanced solid-state nuclear magnetic resonance (ssNMR) spectroscopy and linear-scaling quantum mechanics/molecular mechanics (QM/MM) methods, this single residue is shown to be responsible for a variety of synergistically coupled structural and electrostatic changes along the retinal polyene chain, ionone ring, and within the binding pocket. They collectively explain the observed color shift. Furthermore, analysis of the differences in chemical shift between nuclei within the same residues in green and blue proteorhodopsins also reveals a correlation with the respective degree of conservation. Our data show that the highly conserved color change mainly affects other highly conserved residues, illustrating a high degree of robustness of the color phenotype to sequence variation.

Electrostatic catalysis of a click reaction in a microfluidic cell

Electric fields have been highlighted as a smart reagent in nature's enzymatic machinery, as they can directly trigger or accelerate chemical processes with stereo- and regio-specificity. In enzymatic catalysis, controlled mass transport of chemical species is also key in facilitating the availability of reactants in the active reaction site. However, recent progress in developing a clean catalysis that profits from oriented electric fields is limited to theoretical and experimental studies at the single molecule level, where both the control over mass transport and scalability cannot be tested. Here, we quantify the electrostatic catalysis of a prototypical Huisgen cycloaddition in a large-area electrode surface and directly compare its performance to the conventional Cu(I) catalysis. Our custom-built microfluidic cell enhances reagent transport towards the electrified reactive interface. This continuous-flow microfluidic electrostatic reactor is an example of an electric-field driven platform where clean large-scale electrostatic catalytic processes can be efficiently implemented and regulated.

Studying Local Electrostatics by Terahertz Spectroscopy Using Amines as a Probe

In a previous study[1] we could show that a large amplitude mode of the zwitterion glycine can serve as a sensitive probe for protonation and allows to deduce local pKa values. Here we show that the underlying concept is more general: We present the results of a pH dependent measurement of Terahertz-FTIR (THz-FTIR) spectra of solvated amines, i. e. Diethylamine (DEA), Triethylamine (TEA), and Diisopropylamine (DiPA). We show that amines serve as a sensitive, label free probe for local protonation. Protonation of the amines yield intensity changes which can be quantified by precise THz spectroscopy (30 cm-1 -450 cm-1 ). A detailed analysis allows us to correlate the titration spectra of solvated amines in the THz range with pKa values. This demonstrates the potential of THz spectroscopy to probe the charge state of biomolecules in water in a label free manner.

Computational Characterization of the Reactivity of Compound I in Unspecific Peroxygenases

Unspecific peroxygenases (UPOs) are emerging as promising biocatalysts for selective oxyfunctionalization of unactivated C-H bonds. However, their potential in large-scale synthesis is currently constrained by suboptimal chemical selectivity. Improving the selectivity of UPOs requires a deep understanding of the molecular basis of their catalysis. Recent molecular simulations have sought to unravel UPO's selectivity and inform their design principles. However, most of these studies focused on substrate-binding poses. Few researchers have investigated how the reactivity of CpdI, the principal oxidizing intermediate in the catalytic cycle, influences selectivity in a realistic protein environment. Moreover, the influence of protein electrostatics on the reaction kinetics of CpdI has also been largely overlooked. To bridge this gap, we used multiscale simulations to interpret the regio- and enantioselective hydroxylation of the n-heptane substrate catalyzed by Agrocybe aegerita UPO (AaeUPO). We comprehensively characterized the energetics and kinetics of the hydrogen atom-transfer (HAT) step, initiated by CpdI, and the subsequent oxygen rebound step forming the product. Notably, our approach involved both free energy and potential energy evaluations in a quantum mechanics/molecular mechanics (QM/MM) setting, mitigating the dependence of results on the choice of initial conditions. These calculations illuminate the thermodynamics and kinetics of the HAT and oxygen rebound steps. Our findings highlight that both the conformational selection and the distinct chemical reactivity of different substrate hydrogen atoms together dictate the regio- and enantio-selectivity. Building on our previous study of CpdI's formation in AaeUPO, our results indicate that the HAT step is the rate-limiting step in the overall catalytic cycle. The subsequent oxygen rebound step is swift and retains the selectivity determined by the HAT step. We also pinpointed several polar and charged amino acid residues whose electrostatic potentials considerably influence the reaction barrier of the HAT step. Notably, the Glu196 residue is pivotal for both the CpdI's formation and participation in the HAT step. Our research offers in-depth insights into the catalytic cycle of AaeUPO, which will be instrumental in the rational design of UPOs with enhanced properties.

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.

Enhanced active-site electric field accelerates enzyme catalysis

The design and improvement of enzymes based on physical principles remain challenging. Here we demonstrate that the principle of electrostatic catalysis can be leveraged to substantially improve a natural enzyme's activity. We enhanced the active-site electric field in horse liver alcohol dehydrogenase by replacing the serine hydrogen-bond donor with threonine and replacing the catalytic Zn2+ with Co2+. Based on the electric field enhancement, we make a quantitative prediction of rate acceleration-50-fold faster than the wild-type enzyme-which was in close agreement with experimental measurements. The effects of the hydrogen bonding and metal coordination, two distinct chemical forces, are described by a unified physical quantity-electric field, which is quantitative, and shown here to be additive and predictive. These results suggest a new design paradigm for both biological and non-biological catalysts.

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

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

Understanding Cysteine Reactivity in Protein Environments with Electric Fields

The role cysteine residues play in proteins is mediated by their protonation state, whereby the thiolate form of the side chain is highly reactive while the thiol form is more inert. However, the pKa of cysteine residues is hard to predict as it can differ widely from its reference value in solution, an effect that is accentuated by local effects in the heterogeneous protein environment. Here, we present a new approach to the prediction of cysteine reactivity based on electric field calculations at the thiol/thiolate group. We validated our approach by predicting the protonation state of cysteine residues in different protein environments (in the active site, at the protein surface, and buried within the protein interior), including Cys-25 in papaya protease omega, which was proven problematic for the more traditional constant pH molecular dynamics (MD) technique. We predict pKa shifts consistent with experimental observations, and the decomposition of the electric fields into contributions from molecular fragments provides a direct handle to rationalize local pH and pKa effects in proteins without introducing parameters other than those of the force field used for MD simulations.

Differences in Conformational Sampling and Intrinsic Electric Fields Drive Ion Binding in Telomeric and TERRA G-Quadruplexes

The formation of G-quadruplexes (GQs) occurs in guanine-rich sequences of DNA and RNA, producing highly stable and structurally diverse noncanonical nucleic acid structures. GQs play crucial roles in regulating transcription, translation, and replication and maintaining the genome, among others; thus, changes to their structures can lead to diseases such as cancer. Previous studies using polarizable molecular dynamics simulations have shown differences in ion binding properties between telomeric and telomeric repeat-containing RNA GQs despite architectural similarities. Here, we used volume-based metadynamics and repulsive potential simulations in conjunction with polarizable force fields to quantify the impact of ion binding on the GQ dynamics and ion binding free energies. Furthermore, we describe how GQs exert electric fields on their surroundings to link dynamics with variations in the electronic structure. Our findings provide new insights into the energetic, physical, and conformational properties of GQs and expose subtle but important differences between DNA and RNA GQs with the same fold.

Internal charge-transfer in a metal-catalyzed oxidative addition reaction turns an inhibitive electric field stimulus to catalytic

In a metal-catalyzed oxidative addition, an oriented external electric field (EEF) catalyzes the reaction along one direction and inhibits it when applied in the opposite direction. Beyond a threshold value, the inhibitory direction becomes catalyzing by swapping the metal-to-ligand charge transfer (MLCT) to ligand-to-metal charge-transfer (LMCT) or vice versa. The change in direction of the charge-transfer mechanism triggers the inversion of the dipole moment along the reaction axis, that results in the resurgence of catalysis. The charge-transfer mechanism in metal-catalyzed oxidative addition is tunable by EEF.

From random to rational: improving enzyme design through electric fields, second coordination sphere interactions, and conformational dynamics

Enzymes are versatile and efficient biological catalysts that drive numerous cellular processes, motivating the development of enzyme design approaches to tailor catalysts for diverse applications. In this perspective, we investigate the unique properties of natural, evolved, and designed enzymes, recognizing their strengths and shortcomings. We highlight the challenges and limitations of current enzyme design protocols, with a particular focus on their limited consideration of long-range electrostatic and dynamic effects. We then delve deeper into the impact of the protein environment on enzyme catalysis and explore the roles of preorganized electric fields, second coordination sphere interactions, and protein dynamics for enzyme function. Furthermore, we present several case studies illustrating successful enzyme-design efforts incorporating enzyme strategies mentioned above to achieve improved catalytic properties. Finally, we envision the future of enzyme design research, spotlighting the challenges yet to be overcome and the synergy of intrinsic electric fields, second coordination sphere interactions, and conformational dynamics to push the state-of-the-art boundaries.

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.

Direct measurements of biomolecular electrostatics through experiments

Biomolecular electrostatics has been a subject of computational investigations based on 3D structures. This situation is changing because emerging experimental tools allow us to quantitatively investigate biomolecular electrostatics without any use of structure information. Now, electrostatic potentials around biomolecules can directly be measured for many residues simultaneously by nuclear magnetic resonance (NMR) spectroscopy. This NMR method can be used to study electrostatic aspects of various processes, including macromolecular association and liquid-liquid phase separation. Applications to structurally flexible biomolecules such as intrinsically disordered proteins are particularly useful. The new tools also facilitate examination of theoretical models and methods for biomolecular electrostatics.

Assessment of an Electrostatic Energy-Based Charge Model for Modeling the Electrostatic Interactions in Water Solvent

The protein force field based on the restrained electrostatic potential (RESP) charges has limitations in accurately describing hydrogen bonding interactions in proteins. To address this issue, we propose an alternative approach called the electrostatic energy-based charges (EEC) model, which shows improved performance in describing electrostatic interactions (EIs) of hydrogen bonds in proteins. In this study, we further investigate the performance of the EEC model in modeling EIs in water solvent. Our findings demonstrate that the fixed EEC model can effectively reproduce the quantum mechanics/molecular mechanics (QM/MM)-calculated EIs between a water molecule and various water solvent environments. However, to achieve the same level of computational accuracy, the electrostatic potential (ESP) charge model needs to fluctuate according to the electrostatic environment. Our analysis indicates that the requirement for charge adjustments depends on the specific mathematical and physical representation of EIs as a function of the environment for deriving charges. By comparing with widely used empirical water models calibrated to reproduce experimental properties, we confirm that the performance of the EEC model in reproducing QM/MM EIs is similar to that of general purpose TIP4P-like water models such as TIP4P-Ew and TIP4P/2005. When comparing the computed 10,000 distinct EI values within the range of -40 to 0 kcal/mol with the QM/MM results calculated at the MP2/aug-cc-pVQZ/TIP3P level, we noticed that the mean unsigned error (MUE) for the EEC model is merely 0.487 kcal/mol, which is remarkably similar to the MUE values of the TIP4P-Ew (0.63 kcal/mol) and TIP4P/2005 (0.579 kcal/mol) models. However, both the RESP method and the TIP3P model exhibit a tendency to overestimate the EIs, as evidenced by their higher MUE values of 1.761 and 1.293 kcal/mol, respectively. EEC-based molecular dynamics simulations have demonstrated that, when combined with appropriate van der Waals parameters, the EEC model can closely reproduce oxygen-oxygen radial distribution function and density of water, showing a remarkable similarity to the well-established TIP4P-like empirical water models. Our results demonstrate that the EEC model has the potential to build force fields with comparable accuracy to more sophisticated empirical TIP4P-like water models.

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.

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.

MOFs as Versatile Catalysts: Synthesis Strategies and Applications in Value-Added Compound Production

Catalysts played a crucial role in advancing modern human civilization, from ancient times to the industrial revolution. Due to high cost and limited availability of traditional catalysts, there is a need to develop cost-effective, high-activity, and nonprecious metal-based electrocatalysts. Metal-organic frameworks (MOFs) have emerged as an ideal candidate for heterogeneous catalysis due to their physicochemical properties, hybrid inorganic/organic structures, uncoordinated metal sites, and accessible organic sections. MOFs are high nanoporous crystalline materials that can be used as catalysts to facilitate polymerization reactions. Their chemical and structural diversity make them effective for various reactions compared to traditional catalysts. MOFs have been applied in gas storage and separation, ion-exchange, drug delivery, luminescence, sensing, nanofilters, water purification, and catalysis. The review focuses on MOF-enabled heterogeneous catalysis for value-added compound production, including alcohol oxidation, olefin oligomerization, and polymerization reactions. MOFs offer tunable porosity, high spatial density, and single-crystal XRD control over catalyst properties. In this review, MOFs were focused on reactions of CO2 fixation, CO2 reduction, and photoelectrochemical water splitting. Overall, MOFs have great potential as versatile catalysts for diverse applications in the future.

Differences in Conformational Sampling and Intrinsic Electric Fields Drive Ion Binding in Telomeric and TERRA G-Quadruplexes

The formation of G-quadruplexes (GQs) occurs in guanine-rich sequences of DNA and RNA, producing highly stable and structurally diverse noncanonical nucleic acid structures. GQs play crucial roles in regulating transcription, translation, and replication; and maintaining the genome, among others, thus changes to their structures can lead to diseases such as cancer. Previous studies using polarizable molecular dynamics simulations have shown differences in ion binding properties between telomeric and TERRA GQs despite architectural similarities. Here, we used volume-based metadynamics and repulsive potential simulations in conjunction with polarizable force fields to quantify the impact of ion binding on GQ dynamics and ion binding free energies. Furthermore, we describe how GQs exert electric fields on their surroundings to link dynamics with variations in electronic structure. Our findings provide new insights into the energetic, physical, and conformational properties of GQs and expose subtle, but important, differences between DNA and RNA GQs with the same fold.

A simple geometrical model of the electrostatic environment around the catalytic center of the ribosome and its significance for the elongation cycle kinetics

The central function of the large subunit of the ribosome is to catalyze peptide bond formation. This biochemical reaction is conducted at the peptidyl transferase center (PTC). Experimental evidence shows that the catalytic activity is affected by the electrostatic environment around the peptidyl transferase center. Here, we set up a minimal geometrical model fitting the available x-ray solved structures of the ribonucleic cavity around the catalytic center of the large subunit of the ribosome. The purpose of this phenomenological model is to estimate quantitatively the electrostatic potential and electric field that are experienced during the peptidyl transfer reaction. At least two reasons motivate the need for developing this quantification. First, we inquire whether the electric field in this particular catalytic environment, made only of nucleic acids, is of the same order of magnitude as the one prevailing in catalytic centers of the proteic enzymes counterparts. Second, the protein synthesis rate is dependent on the nature of the amino acid sequentially incorporated in the nascent chain. The activation energy of the catalytic reaction and its detailed kinetics are shown to be dependent on the mechanical work exerted on the amino acids by the electric field, especially when one of the four charged amino acid residues (R, K, E, D) has previously been incorporated at the carboxy-terminal end of the peptidyl-tRNA. Physical values of the electric field provide quantitative knowledge of mechanical work, activation energy and rate of the peptide bond formation catalyzed by the ribosome. We show that our theoretical calculations are consistent with two independent sets of previously published experimental results. Experimental results for E.coli in the minimal case of the dipeptide bond formation when puromycin is used as the final amino acid acceptor strongly support our theoretically derived reaction time courses. Experimental Ribo-Seq results on E. coli and S. cerevisiae comparing the residence time distribution of ribosomes upon specific codons are also well accounted for by our theoretical calculations. The statistical queueing time theory was used to model the ribosome residence time per codon during nascent protein elongation and applied for the interpretation of the Ribo-Seq data. The hypo-exponential distribution fits the residence time observed distribution of the ribosome on a codon. An educated deconvolution of this distribution is used to estimate the rates of each elongation step in a codon specific manner. Our interpretation of all these results sheds light on the functional role of the electrostatic profile around the PTC and its impact on the ribosome elongation cycle.

Research on Electric Field-Induced Catalysis Using Single-Molecule Electrical Measurement

The role of catalysis in controlling chemical reactions is crucial. As an important external stimulus regulatory tool, electric field (EF) catalysis enables further possibilities for chemical reaction regulation. To date, the regulation mechanism of electric fields and electrons on chemical reactions has been modeled. The electric field at the single-molecule electronic scale provides a powerful theoretical weapon to explore the dynamics of individual chemical reactions. The combination of electric fields and single-molecule electronic techniques not only uncovers new principles but also results in the regulation of chemical reactions at the single-molecule scale. This perspective focuses on the recent electric field-catalyzed, single-molecule chemical reactions and assembly, and highlights promising outlooks for future work in single-molecule catalysis.

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.

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.

Protein Dynamics and Enzymatic Catalysis

This Perspective presents a review of our work and that of others in the highly controversial topic of the coupling of protein dynamics to reaction in enzymes. We have been involved in studying this topic for many years. Thus, this perspective will naturally present our own views, but it also is designed to present an overview of the variety of viewpoints of this topic, both experimental and theoretical. This is obviously a large and contentious topic.

Underlying Role of Hydrophobic Environments in Tuning Metal Elements for Efficient Enzyme Catalysis

The catalytic functions of metalloenzymes are often strongly correlated with metal elements in the active sites. However, dioxygen-activating nonheme quercetin dioxygenases (QueD) are found with various first-row transition-metal ions when metal swapping inactivates their innate catalytic activity. To unveil the molecular basis of this seemingly promiscuous yet metal-specific enzyme, we transformed manganese-dependent QueD into a nickel-dependent enzyme by sequence- and structure-based directed evolution. Although the net effect of acquired mutations was primarily to rearrange hydrophobic residues in the active site pocket, biochemical, kinetic, X-ray crystallographic, spectroscopic, and computational studies suggest that these modifications in the secondary coordination spheres can adjust the electronic structure of the enzyme-substrate complex to counteract the effects induced by the metal substitution. These results explicitly demonstrate that such noncovalent interactions encrypt metal specificity in a finely modulated manner, revealing the underestimated chemical power of the hydrophobic sequence network in enzyme catalysis.

Voltage-Dependent FTIR and 2D Infrared Spectroscopies within the Electric Double Layer Using a Plasmonic and Conductive Electrode

Strong electric fields exist between the electric double layer and charged surfaces. These fields impact molecular structures and chemistry at interfaces. We have developed a transparent electrode with infrared plasmonic enhancement sufficient to measure FTIR and two-dimensional infrared spectra at submonolayer coverages on the surface to which a voltage can be applied. Our device consists of an infrared transparent substrate, a 10-20 nm layer of conductive indium tin oxide (ITO), an electrically resistive layer of 3-5 nm Al₂O₃, and a 3 nm layer of nonconductive plasmonic gold. The materials and thicknesses are set to maximize the surface number density of the monolayer molecules, electrical conductivity, and plasmonic enhancement while minimizing background signals and avoiding Fano line shape distortions. The design was optimized by iteratively characterizing the material roughness and thickness with atomic force microscopy and electron microscopy and by monitoring the plasmon resonance enhancement with spectroscopy. The design is robust to repeated fabrication. This new electrode is tested on nitrile functional groups using a monolayer of 4-mercaptobenzonitrile as well as on CO and CC stretching modes using 4-mercaptobenzoic acid methyl ester. A voltage-dependent Stark shift is observed on both monolayers. We also observe that the transition dipole strength of the CN mode scales linearly with the applied voltage, providing a second way of measuring the surface electric field strength. We anticipate that this cell will enable many new voltage-dependent infrared experiments under applied voltages.

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.

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.

Geometry of Charge Density as a Reporter on the Role of the Protein Scaffold in Enzymatic Catalysis: Electrostatic Preorganization and Beyond

Enzymes host active sites inside protein macromolecules, which have diverse, often incredibly complex, and atom-expensive structures. It is an outstanding question what the role of these expensive scaffolds might be in enzymatic catalysis. Answering this question is essential to both enzymology and the design of artificial enzymes with proficiencies that will match those of the best natural enzymes. Protein rigidifying the active site, contrasted with the dynamics and vibrational motion promoting the reaction, as well as long-range electrostatics (also known as electrostatic preorganization) were all proposed as central contributions of the scaffold to the catalysis. Here, we show that all these effects inevitably produce changes in the quantum mechanical electron density in the active site, which in turn defines the reactivity. The phenomena are therefore fundamentally inseparable. The geometry of the electron density-a scalar field characterized by a number of mathematical features such as critical points-is a rigorous and convenient descriptor of enzymatic catalysis and a reporter on the role of the protein. We show how this geometry can be analyzed, linked to the reaction barriers, and report in particular on intramolecular electric fields in enzymes. We illustrate these tools on the studies of electrostatic preorganization in several representative enzyme classes, both natural and artificial. We highlight the forward-looking aspects of the approach.

Quantum Tunneling in Reactions Modulated by External Electric Fields: Reactivity and Selectivity

Quantum tunneling and external electric fields (EEFs) can promote some reactions. However, the synergetic effect of an EEF on a tunneling-involving reaction and its temperature-dependence is not very clear. In this study, we extensively investigated how EEFs affect three reactions that involve hydrogen- or (ground- and excited-state) carbon-tunneling using reliable DFT, DLPNO-CCSD(T1), and variational transition-state theory methods. Our study revealed that oriented EEFs can significantly reduce the barrier and corresponding barrier width (and vice versa) through more electrostatic stabilization in transition states. These EEF effects enhance the nontunneling and tunneling-involving rates. Such EEF effects also decrease the crossover temperatures and quantum tunneling contribution, albeit with lower and thinner barriers. Moreover, EEFs can modulate and switch on/off the tunneling-driven 1,2-H migration of hydroxycarbenes under cryogenic conditions. Furthermore, our study predicts for the first time that EEF/tunneling synergy can control the chemo- or site-selectivity of one molecule bearing two similar/same reactive sites.

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.

External fields in conceptual density functional theory

The necessity of the recent incorporation of new external variables in the context of conceptual DFT (CDFT) is discussed based on the ever-increasing portfolio of experimental reaction conditions in the endeavor of experimentalists to synthesize new molecules with unprecedented properties. Electric and magnetic fields (ε and B), mechanical forces (F), and confinement are proposed as valuable new variables, extending conventional CDFT and its associated response functions. A finite field approach is used to calculate the evolution of both global and local descriptors in a selected series of atomic and molecular applications, and from it derive new response function involving, with one exception, the first derivative to the field considered. The electric field results, displaying, for example, a case of a field-induced enantioselectivity in the Fukui function, may be instrumental in the recent upsurge of chemistry in oriented external electric fields. The study of atomic electronegativity and hardness in magnetic fields displays a piecewise behavior, associated to configurational jumps upon increasing field strength and reveals an overall compression of their ranges for stronger fields, which may be guiding upon investigating chemistry in extremely high fields like in white dwarfs. The evolution of the electronegativity and hardness of diatomics under mechanical force can elegantly be traced back to differences in their equilibrium distance in the neutral, cationic, and anionic state. The well-known reduction of the polarizability under confinement can be seen as a fore-runner of the increasing hardness of atoms under pressure, presently under investigation. Periodicity showing up in a spontaneous way in the variety of properties is a leitmotiv in this study, as well as the interconnections/analogies between the different response functions.

Electro-detachment of kinesin motor domain from microtubule in silico

Kinesin is a motor protein essential in cellular functions, such as intracellular transport and cell-division, as well as for enabling nanoscopic transport in bio-nanotechnology. Therefore, for effective control of function for nanotechnological applications, it is important to be able to modify the function of kinesin. To circumvent the limitations of chemical modifications, here we identify another potential approach for kinesin control: the use of electric forces. Using full-atom molecular dynamics simulations (247,358 atoms, total time ∼ 4.4 μs), we demonstrate, for the first time, that the kinesin-1 motor domain can be detached from a microtubule by an intense electric field within the nanosecond timescale. We show that this effect is field-direction dependent and field-strength dependent. A detailed analysis of the electric forces and the work carried out by electric field acting on the microtubule-kinesin system shows that it is the combined action of the electric field pulling on the β-tubulin C-terminus and the electric-field-induced torque on the kinesin dipole moment that causes kinesin detachment from the microtubule. It is shown, for the first time in a mechanistic manner, that an electric field can dramatically affect molecular interactions in a heterologous functional protein assembly. Our results contribute to understanding of electromagnetic field-biomatter interactions on a molecular level, with potential biomedical and bio-nanotechnological applications for harnessing control of protein nanomotors.

4-Cyanotryptophan as a Sensitive Fluorescence Probe of Local Electric Field of Proteins

Electrostatic interactions are key determinants of protein structure, dynamics, and function. Since protein electrostatics are nonuniform, assessment of the internal electric fields (EFs) of proteins requires spatial resolution at the amino acid residue level. In this regard, vibrational Stark spectroscopy, in conjunction with various unnatural amino acid-based vibrational probes, has become a common method for site-specific interrogation of protein EFs. However, application of this method is often limited to proteins with relatively high solubility, due to the intrinsically low oscillator strength of vibrational transitions. Therefore, it would be useful to develop an alternative method that can overcome this limitation. To this end, we show that, using solvatochromic study and molecular dynamics simulations, the frequency of maximum emission intensity of the fluorophore of 4-cyanotryptophan (4CN-Trp), 3-methyl-1H-indole-4-carbonitrile, exhibits a linear dependence on the local EF. Since the absorption and emission spectra of 4CN-Trp are easily distinguishable from those of naturally occurring aromatic amino acids, we believe that this linear relationship provides an easier and more sensitive means to determine the local EF of proteins.

Revisiting the Burden Borne by Fumarase: Enzymatic Hydration of an Olefin

Fumarate hydratase (FH) is a remarkable catalyst that decreases the free energy of the catalyzed reaction by 30 kcal mol^-1, much larger than most exceptional enzymes with extraordinary catalytic rates. Two classes of FH are observed in nature: class-I and class-II, which have different folds, yet catalyze the same reversible hydration/dehydration reaction of the dicarboxylic acids fumarate/malate, with equal efficiencies. Using class-I FH from the hyperthermophilic archaeon Methanocaldococcus jannaschii (Mj) as a model along with comparative analysis with the only other available class-I FH structure from Leishmania major (Lm), we provide insights into the molecular mechanism of catalysis in this class of enzymes. The structure of MjFH apo-protein has been determined, revealing that large intersubunit rearrangements occur across apo- and holo-protein forms, with a largely preorganized active site for substrate binding. Site-directed mutagenesis of active site residues, kinetic analysis, and computational studies, including density functional theory (DFT) and natural population analysis, together show that residues interacting with the carboxylate group of the substrate play a pivotal role in catalysis. Our study establishes that an electrostatic network at the active site of class-I FH polarizes the substrate fumarate through interactions with its carboxylate groups, thereby permitting an easier addition of a water molecule across the olefinic bond. We propose a mechanism of catalysis in FH that occurs through transition-state stabilization involving the distortion of the electronic structure of the substrate olefinic bond mediated by the charge polarization of the bound substrate at the enzyme active site.

Connecting Conformational Motions to Rapid Dynamics in Human Purine Nucleoside Phosphorylase

The influence of protein motions on enzyme catalysis remains a topic of active discussion. Protein motions occur across a variety of time scales, from vibrational fluctuations in femtoseconds, to collective motions in milliseconds. There have been numerous studies that show conformational motions may assist in catalysis, protein folding, and substrate specificity. It is also known through transition path sampling studies that rapid promoting vibrations contribute to enzyme catalysis. Human purine nucleoside phosphorylase (PNP) is one enzyme that contains both an important conformational motion and a rapid promoting vibration. The slower motion in this enzyme is associated with a loop motion, that when open allows substrate entry and product release but closes over the active site during catalysis. We examine the differences between an unconstrained PNP structure and a PNP structure with constraints on the loop motion. To investigate possible coupling between the slow and fast protein dynamics, we employed transition path sampling, reaction coordinate identification, electric field calculations, and free energy calculations reported here.

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.