Electric Fields at the Active Site of an Enzyme: Direct Comparison of Experiment with Theory

Methods for Theoretical Treatment of Local Fields in Proteins and Enzymes

Electric fields generated by protein scaffolds are crucial in enzymatic catalysis. This review surveys theoretical approaches for detecting, analyzing, and comparing electric fields, electrostatic potentials, and their effects on the charge density within enzyme active sites. Pioneering methods like the empirical valence bond approach rely on evaluating ionic and covalent resonance forms influenced by the field. Strategies employing polarizable force fields also facilitate field detection. The vibrational Stark effect connects computational simulations to experimental Stark spectroscopy, enabling direct comparisons. We highlight how protein dynamics induce fluctuations in local fields, influencing enzyme activity. Recent techniques assess electric fields throughout the active site volume rather than only at specific bonds, and machine learning helps relate these global fields to reactivity. Quantum theory of atoms in molecules captures the entire electron density landscape, providing a chemically intuitive perspective on field-driven catalysis. Overall, these methodologies show protein-generated fields are highly dynamic and heterogeneous, and understanding both aspects is critical for elucidating enzyme mechanisms. This holistic view empowers rational enzyme engineering by tuning electric fields, promising new avenues in drug design, biocatalysis, and industrial applications. Future directions include incorporating electric fields as explicit design targets to enhance catalytic performance and biochemical functionalities.

Electrostatically tuning radical addition and atom abstraction reactions with distonic radical ions

Although electrostatic catalysis can enhance the kinetics and selectivity of reactions to produce greener synthetic processes, the highly directional nature of electrostatic interactions has limited widespread application. In this study, the influence of oriented electric fields (OEF) on radical addition and atom abstraction reactions are systematically explored with ion-trap mass spectrometry using structurally diverse distonic radical ions that maintain spatially separated charge and radical moieties. When installed on rigid molecular scaffolds, charged functional groups lock the magnitude and orientation of the internal electric field with respect to the radical site, creating an OEF which tunes the reactivity across the set of gas-phase carbon-centred radical reactions. In the first case, OEFs predictably accelerate and decelerate the rate of molecular oxygen addition to substituted phenyl, adamantyl, and cubyl radicals, depending on the polarity of the charged functional group and dipole orientation. In the second case, OEFs modulate competition between chlorine and hydrogen atom abstraction from chloroform based on interactions between charge polarity, dipole orientation, and radical polarizability. Importantly, this means the same charge polarity can induce different changes to reaction selectivity. Quantum chemical calculations of these reactions with DSD-PBEP86-D3(BJ)/aug-cc-pVTZ show correlations between the barrier heights and the experimentally determined reaction kinetics. Field effects are consistent between phenyl and cubyl scaffolds, pointing to through-space rather than through-bond field effects, congruent with computations showing that the same effects can be mimicked by point charges. These results experimentally demonstrate how internal OEFs generated by carefully placed charged functional groups can systematically control radical reactions.

Alkali metal cation adsorption-induced surface polarization in polymeric carbon nitride for enhanced photocatalytic hydrogen peroxide production

Photocatalytic hydrogen peroxide (H2O2) generation on the catalyst surface from oxygen is an electron-demanding process, making the construction of an electron-rich surface highly advantageous. In this study, a localized electric field was observed on the surface of polymeric carbon nitride (g-C3N4) when alkali metal cations were adsorbed onto it. These fields effectively inhibited surface carrier recombination and extended their lifespan, thereby enhancing H2O2 production. As a result, g-C3N4 achieved a superior H2O2 yield of 2.25 mM after 1 h in a 0.25 M K+ solution, which was 2.06 times greater than that (1.09 mM) achieved in a pure solvent. Notably, the increase in photocatalytic efficiency showed a remarkable dependence on ion species. At low concentrations, H2O2 generation efficiency was in the order of Li+ < Na+ < K+ < Rb+ < Cs+. However, after optimizing the ion concentration, the highest H2O2 production was achieved in a solution containing K+ instead of Cs+. Molecular dynamics simulations and temperature-dependent photocatalysis experiments revealed that the synergistic interaction between adsorption energy and adsorption distance was crucial in governing the extent to which alkali metal cation adsorption enhanced g-C3N4 photocatalytic H2O2 production. This study provides theoretical insights for the design of materials for electron-demanding photocatalysis and aids in understanding variations in photocatalytic behavior in natural waters.

Utilization of Electric Fields to Modulate Molecular Activities on the Nanoscale: From Physical Properties to Chemical Reactions

As a primary energy source, electricity drives broad fields from everyday electronic circuits to industrial chemical catalysis. From a chemistry viewpoint, studying electric field effects on chemical reactivity is highly important for revealing the intrinsic mechanisms of molecular behaviors and mastering chemical reactions. Recently, manipulating the molecular activity using electric fields has emerged as a new research field. In addition, because integration of molecules into electronic devices has the natural complementary metal-oxide-semiconductor compatibility, electric field-driven molecular devices meet the requirements for both electronic device miniaturization and precise regulation of chemical reactions. This Review provides a timely and comprehensive overview of recent state-of-the-art advances, including theoretical models and prototype devices for electric field-based manipulation of molecular activities. First, we summarize the main approaches to providing electric fields for molecules. Then, we introduce several methods to measure their strengths in different systems quantitatively. Subsequently, we provide detailed discussions of electric field-regulated photophysics, electron transport, molecular movements, and chemical reactions. This review intends to provide a technical manual for precise molecular control in devices via electric fields. This could lead to development of new optoelectronic functions, more efficient logic processing units, more precise bond-selective control, new catalytic paradigms, and new chemical reactions.

Hydrogen Bond Blueshifts in Nitrile Vibrational Spectra Are Dictated by Hydrogen Bond Geometry and Dynamics

Vibrational Stark effect (VSE) spectroscopy has become one of the most important experimental approaches to determine the strength of noncovalent, electrostatic interactions in chemistry and biology and to quantify their influence on structure and reactivity. Nitriles (C≡N) have been widely used as VSE probes, but their application has been complicated by an anomalous hydrogen bond (HB) blueshift which is not encompassed within the VSE framework. We present an empirical model describing the anomalous HB blueshift in terms of H-bonding geometry, i.e., as a function of HB distance and angle with respect to the C≡N group. This model is obtained by comparing vibrational observables from density functional theory and electrostatics from the polarizable AMOEBA force field, and it provides a physical explanation for the HB blueshift in terms of underlying multipolar and Pauli repulsion contributions. Additionally, we compare predicted blueshifts with experimental results and find our model provides a useful, direct framework to analyze HB geometry for rigid HBs, such as within proteins or chemical frameworks. In contrast, nitriles in highly dynamic H-bonding environments like protic solvents are no longer a function solely of geometry; this is a consequence of motional narrowing, which we demonstrate by simulating IR spectra. Overall, when HB geometry and dynamics are accounted for, an excellent correlation is found between observed and predicted HB blueshifts. This correlation includes different types of nitriles and HB donors, suggesting that our model is general and can aid in understanding HB blueshifts wherever nitriles can be implemented.

Temperature and external fields in conceptual density functional theory

Until quite recently, Conceptual DFT (CDFT) was mainly based on the energy functional, E[N,v], where the number of electrons N and the external potential v are state variables. One of the strengths of CDFT, however, is the ease with which additional and/or different state variables can be incorporated. Here, the incorporation of new variables-namely temperature and external fields-is discussed, outlining the motivation for these extensions, sketching their theoretical/computational context, and presenting some elucidative examples. Using the Grand Canonical Ensemble, finite temperature can be introduced, ameliorating the N-differentiability problem and leading to new well-behaved chemical reactivity concepts. Incorporating temperature as an additional fundamental variable enables us to formulate a novel suite of "thermodynamic" reactivity descriptors, including important concepts like the heat capacity of electronic systems. The mathematical structure underpinning the set of (well-behaved) finite-temperature reactivity indices can guide the formulation of plausible definitions for local analogs of global descriptors. This endeavor is especially significant in the case of local hardness, a concept that has remained elusive since the inception of CDFT. The ever-increasing portfolio of experimental reaction conditions to creatively synthesize new molecules, needs the introduction of various external "fields" like electric and magnetic fields (ε and B), mechanical forces (F) and pressure (P) to describe the state of the chemical system. The conventional energy functional can be expressed in a general form, E[N,v,X], where X denotes the "field". Response functions to changes in the field can then be defined in analogy to classical thermodynamics. The electric field results display a case of a field-induced selectivity in a reaction channel of the Fukui function. Remarkably, atomic electronegativity and hardness in magnetic fields display a piecewise behavior in magnetic fields, associated to configurational jumps upon increasing field strength. The overall compression of their ranges for stronger fields may be insightful when investigating chemistry in extremely high fields. The electronegativity and hardness of diatomics under mechanical force can be traced back to changes in equilibrium distances in the neutral, cationic and anionic state, parallel with the evolution of an intrinsic atomic volume under pressure.

Linear and Two-Dimensional Infrared Spectroscopy of the Multifunctional Vibrational Probe, 3-(4-Azidophenyl) Propiolonitrile. Deperturbing a Fermi Triad by Isotopic Substitution

Infrared probes are chemical moieties whose vibrational modes are used to obtain spectroscopic information about structural dynamics of complex systems; in particular, of biomacromolecules. Here, we explore the vibrational spectroscopy and dynamics of a reagent, 3-(4-azidophenyl)propiolonitrile (AzPPN), for selectively tagging thiols in protein environments with a multifunctional infrared probe containing both, an azide and a nitrile chromophore. The linear infrared spectrum of AzPPN is heavily perturbed in the antisymmetric azide stretching region as a result of accidental Fermi resonances. Isotopically labeling the azide group at the β-position deperturbs the spectrum considerably and reveals two combination tones that mix with the antisymmetric stretching fundamental into a Fermi triad of hybrid vibrational excitations. Moreover, two-dimensional infrared (2DIR) spectra were recorded for 15Nβ-labeled AzPPN, which reveal waiting-time-dependent spectral shifts of diagonal peaks and dynamic buildups of cross peaks. The 2DIR-spectral evolution is indicative of intramolecular distribution of the pump-induced excess vibrational energy into low-frequency modes of the molecule that are coupled to either the azide or the nitrile stretching transition dipoles. Finally, IR-pump/IR-probe spectra with selective narrowband excitation reveal a time constant of 2.3 ps for intramolecular vibrational redistribution (IVR) and 18 ps for the final energy dissipation into the solvent. The cross-peak dynamics corroborate a notion in which IVR within the AzPPN-molecule is an irreversible process.

Investigation on molecular and biomolecular spectroscopy of the novel 2BCA molecule to analyse its biological activities and binding interaction with nipah viral protein

The molecule of 2-Biphenyl Carboxylic Acid (2BCA), which contains peculiar features, was explored making use of density functional theory (DFT) and experimental approaches in the area of quantum computational research. The optimised structure, atomic charges, vibrational frequencies, electrical properties, electrostatic potential surface (ESP), natural bond orbital analysis and potential energy surface (PES) were obtained applying the B3LYP approach with the 6-311++ G (d,p) basis set.. The 2BCA molecule was examined for possible conformers using a PES scan. The methods applied for spectral analyses included FT-IR, FT-RAMAN, NMR, and UV-Vis results. Vibrational frequencies for all typical modes of vibration were found using the Potential Energy Distribution (PED) data. The UV-Vis spectrum was simulated using the TD-DFT technique, which is also seen empirically. The Gauge-Invariant Atomic Orbital (GIAO) approach was employed to model and study the 13C and 1H NMR spectra of the 2BCA molecule in a CDCL3 solution. The spectra were then exploited experimentally to establish their chemical shifts. To predict the donor and acceptor interaction, the NBO analysis was used. The electrostatic potential surface was employed to anticipate the locations of nucleophilic and electrophilic sites. Hirshfeld surfaces and their related fingerprint plots are exploited for the investigation of intermolecular interactions. Reduced Density Gradient (RDG) helps to measure and illustrate electron correlation effects, offering precise insights into chemical bonding, reactivity, and the electronic structure of 2BCA. According to Lipinski and Veber's drug similarity criteria, 2BCA exhibits the typical physicochemical and pharmacokinetic properties that make it a potential oral pharmaceutical candidate. According to the findings of a molecular docking study, the 2BCA molecule has promise as a treatment agent for the Nipah virus (PDB ID: 6 EB9), which causes severe respiratory and neurological symptoms in humans.

The Evolution of the Acylation Mechanism in β-Lactamase and Rapid Protein Dynamics

β-Lactamases are a class of well-studied enzymes that are known to have existed since billions of years ago, starting as a defense mechanism to stave off competitors and are now enzymes responsible for antibiotic resistance. Using ancestral sequence reconstruction, it is possible to study the crystal structure of a laboratory resurrected 2-3 billion year-old β-lactamase. Comparing the ancestral enzyme to its modern counterpart, a TEM-1 β-lactamase, the structural changes are minor, and it is probable that dynamic effects play an important role in the evolution of function. We used molecular dynamics simulations and employed transition path sampling methods to identify the presence of rate-enhancing dynamics at the femtosecond level in both systems, found that these fast motions are more efficiently coordinated in the modern enzyme, and examined how specific dynamics can pinpoint evolutionary effects that are essential for improving enzymatic catalysis.

SCN as a local probe of protein structural dynamics

The dynamics of lysozyme is probed by attaching -SCN to all alanine residues. The one-dimensional infrared spectra exhibit frequency shifts in the position of the maximum absorption of 4 cm-1, which is consistent with experiments in different solvents and indicates moderately strong interactions of the vibrational probe with its environment. Isotopic substitution 12C → 13C leads to a redshift by -47 cm-1, which agrees quantitatively with experiments for CN-substituted copper complexes in solution. The low-frequency, far-infrared part of the protein spectra contains label-specific information in the difference spectra when compared with the wild type protein. Depending on the position of the labels, local structural changes are observed. For example, introducing the -SCN label at Ala129 leads to breaking of the α-helical structure with concomitant change in the far-infrared spectrum. Finally, changes in the local hydration of SCN-labeled alanine residues as a function of time can be related to the reorientation of the label. It is concluded that -SCN is potentially useful for probing protein dynamics, both in the high-frequency part (CN-stretch) and in the far-infrared part of the spectrum.

Unnatural Amino Acids for Biological Spectroscopy and Microscopy

Due to advances in methods for site-specific incorporation of unnatural amino acids (UAAs) into proteins, a large number of UAAs with tailored chemical and/or physical properties have been developed and used in a wide array of biological applications. In particular, UAAs with specific spectroscopic characteristics can be used as external reporters to produce additional signals, hence increasing the information content obtainable in protein spectroscopic and/or imaging measurements. In this Review, we summarize the progress in the past two decades in the development of such UAAs and their applications in biological spectroscopy and microscopy, with a focus on UAAs that can be used as site-specific vibrational, fluorescence, electron paramagnetic resonance (EPR), or nuclear magnetic resonance (NMR) probes. Wherever applicable, we also discuss future directions.

Nitrile Vibrational Lifetimes as Probes of Local Electric Fields

Optical measurements of electric fields have wide-ranging applications in the fields of chemistry and biology. Previously, such measurements focused on shifts in intensity or frequency. Here, we show that nitrile vibrational lifetimes can report local electric fields through ultrasensitive picosecond mid-infrared-near-infrared double-resonance fluorescence spectro-microscopy on Rhodamine 800. Using a robust convolution fitting approach, we observe that the nitrile vibrational lifetimes are strongly linearly correlated (R2 = 0.841) with solvent reaction fields. Supported by density functional theory, we rationalize this trend through a doorway model of intramolecular vibrational energy redistribution. This work provides new fundamental insights into the nature of vibrational energy flow in large polyatomic molecular systems and establishes a theoretical basis for electric field sensing with vibrational lifetimes, offering a new experimental dimension for probing intracellular electrostatics.

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.

Controlling the charge-transfer dynamics of two-level systems around avoided crossings

Two-level quantum systems are fundamental physical models that continue to attract growing interest due to their crucial role as a building block of quantum technologies. The exact analytical solution of the dynamics of these systems is central to control theory and its applications, such as that to quantum computing. In this study, we reconsider the two-state charge transfer problem by extending and using a methodology developed to study (pseudo)spin systems in quantum electrodynamics contexts. This approach allows us to build a time evolution operator for the charge transfer system and to show new opportunities for the coherent control of the system dynamics, with a particular emphasis on the critical dynamic region around the transition state coordinate, where the avoided crossing of the energy levels occurs. We identify and propose possible experimental implementations of a class of rotations of the charge donor (or acceptor) that endow the electronic coupling matrix element with a time-dependent phase that can be employed to realize controllable coherent dynamics of the system across the avoided level crossing. The analogy of these rotations to reference frame rotations in generalized semiclassical Rabi models is discussed. We also show that the physical rotations in the charge-transfer systems can be performed so as to implement quantum gates relevant to quantum computing. From an exquisitely physical-mathematical viewpoint, our approach brings to light situations in which the time-dependent state of the system can be obtained without resorting to the special functions appearing in the Landau-Zener approach.

Electrochemistry in sensing of molecular interactions of proteins and their behavior in an electric field

Electrochemical methods can be used not only for the sensitive analysis of proteins but also for deeper research into their structure, transport functions (transfer of electrons and protons), and sensing their interactions with soft and solid surfaces. Last but not least, electrochemical tools are useful for investigating the effect of an electric field on protein structure, the direct application of electrochemical methods for controlling protein function, or the micromanipulation of supramolecular protein structures. There are many experimental arrangements (modalities), from the classic configuration that works with an electrochemical cell to miniaturized electrochemical sensors and microchip platforms. The support of computational chemistry methods which appropriately complement the interpretation framework of experimental results is also important. This text describes recent directions in electrochemical methods for the determination of proteins and briefly summarizes available methodologies for the selective labeling of proteins using redox-active probes. Attention is also paid to the theoretical aspects of electron transport and the effect of an external electric field on the structure of selected proteins. Instead of providing a comprehensive overview, we aim to highlight areas of interest that have not been summarized recently, but, at the same time, represent current trends in the field.

Bond-selective fluorescence imaging with single-molecule sensitivity

Bioimaging harnessing optical contrasts and chemical specificity is of vital importance in probing complex biology. Vibrational spectroscopy based on mid-infrared (mid-IR) excitation can reveal rich chemical information about molecular distributions. However, its full potential for bioimaging is hindered by the achievable sensitivity. Here, we report bond selective fluorescence-detected infrared-excited (BonFIRE) spectral microscopy. BonFIRE employs two-photon excitation in the mid-IR and near-IR to upconvert vibrational excitations to electronic states for fluorescence detection, thus encoding vibrational information into fluorescence. The system utilizes tuneable narrowband picosecond pulses to ensure high sensitivity, biocompatibility, and robustness for bond-selective biological interrogations over a wide spectrum of reporter molecules. We demonstrate BonFIRE spectral imaging in both fingerprint and cell-silent spectroscopic windows with single-molecule sensitivity for common fluorescent dyes. We then demonstrate BonFIRE imaging on various intracellular targets in fixed and live cells, neurons, and tissues, with promises for further vibrational multiplexing. For dynamic bioanalysis in living systems, we implement a high-frequency modulation scheme and demonstrate time-lapse BonFIRE microscopy of live HeLa cells. We expect BonFIRE to expand the bioimaging toolbox by providing a new level of bond-specific vibrational information and facilitate functional imaging and sensing for biological investigations.

Revisiting the influence of pH on 1JCαH and chemical shifts of glycine and alanine short oligopeptides

The pH dependence of several NMR parameters of glycine and alanine short oligopeptides has been reported previously in different studies. Here we have thoroughly examined, summarized and demonstrated the dependences of 1H, 13C and 15N chemical shifts and protonation states of amino acids using two-dimensional NMR experiments. Nevertheless, 1JCαH one bond spin-spin coupling constants are more informative and convenient for determination of the position and protonation state of glycine and alanine residue in the oligopeptide chain. In particular, for various oligopeptides (up to six residues), it was shown that the pH dependence of 1JCαH of N-terminal glycine and alanine residues is larger than that of C-terminal groups, and in backbone residues, it is not influenced by pH and only slightly depends on the position of the amino acid residue in the chain.

Dioxygen and glucose force motion of the electron-transfer switch in the iron(III) flavohemoglobin-type nitric oxide dioxygenase

Kinetic and structural investigations of the flavohemoglobin-type NO dioxygenase have suggested critical roles for transient Fe(III)O2 complex formation and O2-forced movements affecting hydride transfer to the FAD cofactor and electron-transfer to the Fe(III)O2 complex. Stark-effect theory together with structural models and dipole and internal electrostatic field determinations provided a semi-quantitative spectroscopic method for investigating the proposed Fe(III)O2 complex and O2-forced movements. Deoxygenation of the enzyme causes Stark effects on the ferric heme Soret and charge-transfer bands revealing the Fe(III)O2 complex. Deoxygenation also elicits Stark effects on the FAD that expose forces and motions that create a more restricted NADH access to FAD for hydride transfer and switch electron-transfer off. Glucose also forces the enzyme toward an off state. Amino acid substitutions at the B10, E7, E11, G8, D5, and F7 positions influence the Stark effects of O2 on resting heme spin states and FAD consistent with the proposed roles of the side chains in the enzyme mechanism. Deoxygenation of ferric myoglobin and hemoglobin A also induces Stark effects on the hemes suggesting a common 'oxy-met' state. The ferric myoglobin and hemoglobin heme spectra are also glucose-responsive. A conserved glucose or glucose-6-phosphate binding site is found bridging the BC-corner and G-helix in flavohemoglobin and myoglobin suggesting novel allosteric effector roles for glucose or glucose-6-phosphate in the NO dioxygenase and O2 storage functions. The results support the proposed roles of a ferric O2 intermediate and protein motions in regulating electron-transfer during NO dioxygenase turnover.

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.

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.

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.

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.

Blue Fluorescence of Cyano-tryptophan Predicts Local Electrostatics and Hydrogen Bonding in Biomolecules

Cyano-tryptophan is an unnatural fluorescent amino acid that emits in the visible region. Along with the structural similarity with tryptophan, the unique photophysical properties of this fluorophore make it an ideal probe for biophysical research. Herein, combining fluorescence spectroscopy, infrared spectroscopy, and molecular dynamics simulations, we show that the cyano-tryptophan's emission energy quantifies the underlying bond-specific noncovalent interactions in terms of the electric field. We further report the use of fluorophore's emission energy to predict its hydrogen bond characteristics. We demonstrate that combining experiments with molecular dynamics simulations can provide the hydrogen bonding status of the nitrile moiety. In addition, we report a method to differentiate between aqueous and nonaqueous hydrogen-bonding partners. Using a phenomenological approach, we demonstrate that the presence of the cyano-indole moiety is responsible for the distinct correlations between the fluorophore's emission and the electrostatic forces on the nitrile bond. As indole is a privileged scaffold for both native amino acids and nucleobases, cyano-indoles will have many multifaceted applications.

Setting Boundaries for Statistical Mechanics

Statistical mechanics has grown without bounds in space. Statistical mechanics of noninteracting point particles in an unbounded perfect gas is widely used to describe liquids like concentrated salt solutions of life and electrochemical technology, including batteries. Liquids are filled with interacting molecules. A perfect gas is a poor model of a liquid. Statistical mechanics without spatial bounds is impossible as well as imperfect, if molecules interact as charged particles, as nearly all atoms do. The behavior of charged particles is not defined until boundary structures and values are defined because charges are governed by Maxwell's partial differential equations. Partial differential equations require boundary structures and conditions. Boundary conditions cannot be defined uniquely 'at infinity' because the limiting process that defines 'infinity' includes such a wide variety of structures and behaviors, from elongated ellipses to circles, from light waves that never decay, to dipolar fields that decay steeply, to Coulomb fields that hardly decay at all. Boundaries and boundary conditions needed to describe matter are not prominent in classical statistical mechanics. Statistical mechanics of bounded systems is described in the EnVarA system of variational mechanics developed by Chun Liu, more than anyone else. EnVarA treatment does not yet include Maxwell equations.

Diversity oriented clicking delivers β-substituted alkenyl sulfonyl fluorides as covalent human neutrophil elastase inhibitors

Diversity Oriented Clicking (DOC) is a discovery method geared toward the rapid synthesis of functional libraries. It combines the best attributes of both classical and modern click chemistries. DOC strategies center upon the chemical diversification of core "SuFExable" hubs-exemplified by 2-Substituted-Alkynyl-1-Sulfonyl Fluorides (SASFs)-enabling the modular assembly of compounds through multiple reaction pathways. We report here a range of stereoselective Michael-type addition pathways from SASF hubs including reactions with secondary amines, carboxylates, 1H-1,2,3-triazole, and halides. These high yielding conjugate addition pathways deliver unprecedented β-substituted alkenyl sulfonyl fluorides as single isomers with minimal purification, greatly enriching the repertoire of DOC and holding true to the fundamentals of modular click chemistry. Further, we demonstrate the potential for biological function - a key objective of click chemistry - of this family of SASF-derived molecules as covalent inhibitors of human neutrophil elastase.

High Power Irradiance Dependence of Charge Species Dynamics in Hybrid Perovskites and Kinetic Evidence for Transient Vibrational Stark Effect in Formamidinium

Hybrid halide perovskites materials have the potential for both photovoltaic and light-emitting devices. Relatively little has been reported on the kinetics of charge relaxation upon intense excitation. In order to evaluate the illumination power density dependence on the charge recombination mechanism, we have applied a femtosecond transient mid-IR absorption spectroscopy with strong excitation to directly measure the charge kinetics via electron absorption. The irradiance-dependent relaxation processes of the excited, photo-generated charge pairs were quantified in polycrystalline MAPbI₃, MAPbBr₃, and (FAPbI₃)_0.97(MAPbBr₃)_0.03 thin films that contain either methylamonium (MA) or formamidinium (FA). This report identifies the laser-generated charge species and provides the kinetics of Auger, bimolecular and excitonic decay components. The inter-band electron-hole (bimolecular) recombination was found to dominate over Auger recombination at very high pump irradiances, up to the damage threshold. The kinetic analysis further provides direct evidence for the carrier field origin of the vibrational Stark effect in a formamidinium containing perovskite material. The results suggest that radiative excitonic and bimolecular recombination in MAPbI₃ at high excitation densities could support light-emitting applications.

Nitrile Infrared Intensities Characterize Electric Fields and Hydrogen Bonding in Protic, Aprotic, and Protein Environments

Nitriles are widely used vibrational probes; however, the interpretation of their IR frequencies is complicated by hydrogen bonding (H-bonding) in protic environments. We report a new vibrational Stark effect (VSE) that correlates the electric field projected on the -C≡N bond to the transition dipole moment and, by extension, the nitrile peak area or integrated intensity. This linear VSE applies to both H-bonding and non-H-bonding interactions. It can therefore be generally applied to determine electric fields in all environments. Additionally, it allows for semiempirical extraction of the H-bonding contribution to the blueshift of the nitrile frequency. Nitriles were incorporated at H-bonding and non-H-bonding protein sites using amber suppression, and each nitrile variant was structurally characterized at high resolution. We exploited the combined information available from variations in frequency and integrated intensity and demonstrate that nitriles are a generally useful probe for electric fields.

Atomistic Simulations for Reactions and Vibrational Spectroscopy in the Era of Machine Learning─Quo Vadis?

Atomistic simulations using accurate energy functions can provide molecular-level insight into functional motions of molecules in the gas and in the condensed phase. This Perspective delineates the present status of the field from the efforts of others and some of our own work and discusses open questions and future prospects. The combination of physics-based long-range representations using multipolar charge distributions and kernel representations for the bonded interactions is shown to provide realistic models for the exploration of the infrared spectroscopy of molecules in solution. For reactions, empirical models connecting dedicated energy functions for the reactant and product states allow statistically meaningful sampling of conformational space whereas machine-learned energy functions are superior in accuracy. The future combination of physics-based models with machine-learning techniques and integration into all-purpose molecular simulation software provides a unique opportunity to bring such dynamics simulations closer to reality.

Impact of Ultrafast Electric Field Changes on Photoreceptor Protein Dynamics

Studies on photoreceptors provide a wealth of information on cofactor and protein dynamics on the microsecond to seconds time-scale. Up to now, ultrafast dynamics addresses mainly the cofactor or chromophore, but ultrafast protein dynamics are poorly understood. Increasing evidence show that protein responses can occur even faster than the cofactor dynamics. The causal reason for the ultrafast protein response cannot be explained by the localized cofactor excitation or its excited-state decay, alone. We propose a Coulomb interaction mechanism started by a shock wave and stabilized by a dipole moment change at least partially responsible for coherent oscillations in proteins, protonation changes, water dislocations, and protein changes prior to and beyond chromophore's excited-state decay. Photoexcitation changes the electron density distribution of the chromophore within a few femtoseconds: The Coulomb shock wave affects polar groups, hydrogen bonds, and protein bound water molecules. The process occurs on a time-scale even faster than excited-state decay of the chromophore. We discuss studies on selected photoreceptors in light of this mechanism and its impact on a detailed understanding of protein dynamics.

ATP Synthase: A Moonlighting Enzyme with Unprecedented Functions

ATP synthase’s intrinsic molecular electrostatic potential (MESP) adds constructively to, and hence reinforces, the chemiosmotic voltage. This ATP synthase voltage represents a new free energy term that appears to have...

Quantum vibration perturbation approach with polyatomic probe in simulating infrared spectra

The quantitative prediction of vibrational spectra of chromophore molecules in solution is challenging and numerous methods have been developed. In this work, we present a quantum vibration perturbation (QVP) approach, which is a procedure that combines molecular quantum vibration and molecular dynamics with perturbation theory. In this framework, an initial Newtonian molecular dynamics simulation is performed, followed by a substitution process to embed molecular quantum vibrational wave functions into the trajectory. The instantaneous vibrational frequency shift at each time step is calculated using the Rayleigh-Schrödinger perturbation theory, where the perturbation operator is the difference in the vibrational potential between the reference chromophore and the perturbed chromophore in the environment. Semi-classical statistical mechanics is employed to obtain the spectral lineshape function. We validated our method using HCOOH·nH₂O (n = 1-2) clusters and HCOOH aqueous solution as examples. The QVP method can be employed for rapid prediction of the vibrational spectrum of a specific mode in solution.

Protein-Ligand Binding Molecular Details Revealed by Terahertz Optical Kerr Spectroscopy: A Simulation Study

Picosecond fast motions and their involvement in the biochemical processes such as protein-ligand binding has engaged significant attention. Terahertz optical Kerr spectroscopy (OKE) has the superior potential to probe these fast motions directly. Application of OKE in protein-ligand binding study is, however, limited by the difficulty of quantitative atomistic interpretation, and the calculation of Kerr spectrum for entire solvated protein complex was considered not yet feasible, due to the lack of one consistent polarizable model for both configuration sampling and polarizability calculation. Here, we analyzed the biochemical relevance of OKE to the lysozyme-triacetylchitotriose binding based on the first OKE simulation using one consistent Drude polarizable model. An analytical multipole and induced dipole scheme was employed to calculate the off-diagonal Drude polarizability more efficiently and accurately. Further theoretical analysis revealed how the subtle twisting and stiffening of aromatic protein residues' spatial arrangement as well as the confinement of small water clusters between ligand and protein cavity due to the ligand binding can be examined using Kerr spectroscopy. Comparison between the signals of bound complex and that of uncorrelated protein/ligand demonstrated that binding action alone has reflection in the OKE spectrum. Our study indicated OKE as a powerful terahertz probe for protein-ligand binding chemistry and dynamics.

Electric Fields Influence Intramolecular Vibrational Energy Relaxation and Line Widths

Intramolecular vibrational energy relaxation (IVR) is fundamentally important to chemical dynamics. We show that externally applied electric fields affect IVR and vibrational line widths by changing the anharmonic couplings and frequency detunings between modes. We demonstrate this effect in benzonitrile for which prior experimental results show a decrease in vibrational line width as a function of applied electric field. We identify three major channels for IVR that depend on electric field. In the dominant channel, the electric field affects the frequency detuning, while in the other two channels, variation of anharmonic couplings as a function of field is the underlying mechanism. Consistent with experimental results, we show that the combination of all channels gives rise to reduced line widths with increasing electric field in benzonitrile. Our results are relevant for controlling IVR with external or internal fields and for gaining a more complete interpretation of line widths of vibrational Stark probes.

Modeling of the Response of Hydrogen Bond Properties on an External Electric Field: Geometry, NMR Chemical Shift, Spin-Spin Scalar Coupling

The response of the geometric and NMR properties of molecular systems to an external electric field has been studied theoretically in a wide field range. It has been shown that this adduct under field approach can be used to model the geometric and spectral changes experienced by molecular systems in polar media if the system in question has one and only one bond, the polarizability of which significantly exceeds the polarizability of other bonds. If this requirement is met, then it becomes possible to model even extreme cases, for example, proton dissociation in hydrogen halides. This requirement is fulfilled for many complexes with one hydrogen bond. For such complexes, this approach can be used to facilitate a detailed analysis of spectral changes associated with geometric changes in the hydrogen bond. For example, in hydrogen-bonded complexes of isocyanide C≡¹⁵N-¹H⋯X, ¹J(¹⁵N¹H) depends exclusively on the N-H distance, while δ(¹⁵N) is also slightly influenced by the nature of X.

Quantification of the electric field inside protein active sites and fullerenes

While electrostatic interactions are exceedingly accountable for biological functions, no simple method exists to directly estimate or measure the electrostatic field in protein active sites. The electrostatic field inside the protein is generally inferred from the shift in the vibrational stretching frequencies of nitrile and thionitrile probes at the active sites through several painstaking and time-consuming experiments like vibrational Stark effect spectroscopy (VSS). Here we present a simple, fast, and reliable methodology, which can efficiently predict the vibrational Stark tuning rates (VSRs) of a large variety of probes within 10% error of the reported experimental data. Our methodology is based on geometry optimization and frequency calculations in the presence of an external electric field to predict the accurate VSR of newly designed nitrile/thionitrile probes. A priori information of VSRs is useful for difficult experiments such as catalytic/enzymatic study and in structural biology. We also applied our methodology successfully to estimate the electric field inside fullerenes and nano-onions, which is encouraging for researchers to adopt it for further applications in materials science and supramolecular chemistry.

Vibrational Stark shift spectroscopy of catalysts under the influence of electric fields at electrode-solution interfaces

External control of chemical processes is a subject of widespread interest in chemical research, including control of electrocatalytic processes with significant promise in energy research. The electrochemical double-layer is the nanoscale region next to the electrode/electrolyte interface where chemical reactions typically occur. Understanding the effects of electric fields within the electrochemical double layer requires a combination of synthesis, electrochemistry, spectroscopy, and theory. In particular, vibrational sum frequency generation (VSFG) spectroscopy is a powerful technique to probe the response of molecular catalysts at the electrode interface under bias. Fundamental understanding can be obtained via synthetic tuning of the adsorbed molecular catalysts on the electrode surface and by combining experimental VSFG data with theoretical modelling of the Stark shift response. The resulting insights at the molecular level are particularly valuable for the development of new methodologies to control and characterize catalysts confined to electrode surfaces. This Perspective article is focused on how systematic modifications of molecules anchored to surfaces report information concerning the geometric, energetic, and electronic parameters of catalysts under bias attached to electrode surfaces.

Heterogeneous electrocatalysis: characterization of interfacial electric field within the electrochemical double layer.

Site-selective dynamics of azidolysozyme

The spectroscopic response of and structural dynamics around all azido-modified alanine residues (AlaN₃) in lysozyme are characterized. It is found that AlaN₃ is a positionally sensitive probe for the local dynamics, covering a frequency range of ∼15 cm^-1 for the center frequency of the line shape. This is consistent with findings from selective replacements of amino acids in PDZ2, which reported a frequency span of ∼10 cm^-1 for replacements of Val, Ala, or Glu by azidohomoalanine. For the frequency fluctuation correlation functions, the long-time decay constants τ₂ range from ∼1 to ∼10 ps, which compares with experimentally measured correlation times of 3 ps. Attaching azide to alanine residues can yield dynamics that decays to zero on the few ps time scale (i.e., static component Δ₀ ∼ 0 ps^-1) or to a remaining, static contribution of ∼0.5 ps^-1 (corresponding to 2.5 cm^-1), depending on the local environment on the 10 ps time scale. The magnitude of the static component correlates qualitatively with the degree of hydration of the spectroscopic probe. Although attaching azide to alanine residues is found to be structurally minimally invasive with respect to the overall protein structure, analysis of the local hydrophobicity indicates that the hydration around the modification site differs for modified and unmodified alanine residues, respectively.

Interfacial Polarization and Ionic Structure at the Ionic Liquid-Metal Interface Studied by Vibrational Spectroscopy and Molecular Dynamics Simulations

Ionic liquids (ILs) have both fundamental and practical value in interfacial science and electrochemistry. However, understanding their behavior near a surface is challenging because of strong Coulomb interactions and large and irregular ionic sizes, which affect both their structure and energetics. To understand this problem, we present a combined experimental and computational study using a vibrational probe molecule, 4-mercaptobenzonitrile, inserted at the junction between a metal and a variety of ILs. The vibrational frequency of the nitrile in the probe molecule reports on the local solvation environment and the electrostatic field at this junction. Within the ethylmethyl imidazolium (EMIM⁺) cation family of ILs, we varied the anions over a range of sizes and types. Complementing our surface spectroscopy, we also ran molecular dynamics simulations of these interfaces to better understand the ionic structures that produced the measured fields. The magnitude of the frequency shifts, and thereby fields, shows a general correlation with the size of anions, with larger anions corresponding to smaller fields. We find that the source of this correlation is partial intercalation of smaller anions into the probe monolayer, resulting in tighter packing of ionic layers near the surface. Larger anions reduce the overall lateral ion packing density near the surface, which reduces the net charge per unit area and explains the smaller observed fields. The insight from this work is important for developing a fundamental picture of concentrated electrolytes near interfaces and can help with designing ILs to create tailored electric fields near an electrode.

Modeling of Solute-Solvent Interactions Using an External Electric Field-From Tautomeric Equilibrium in Nonpolar Solvents to the Dissociation of Alkali Metal Halides

An implicit account of the solvent effect can be carried out using traditional static quantum chemistry calculations by applying an external electric field to the studied molecular system. This approach allows one to distinguish between the effects of the macroscopic reaction field of the solvent and specific solute-solvent interactions. In this study, we report on the dependence of the simulation results on the use of the polarizable continuum approximation and on the importance of the solvent effect in nonpolar solvents. The latter was demonstrated using experimental data on tautomeric equilibria between the pyridone and hydroxypyridine forms of 2,6-di-tert-butyl-4-hydroxy-pyridine in cyclohexane and chloroform.

Multiscale modeling shows that dielectric differences make NaV channels faster than KV channels

The generation of action potentials in excitable cells requires different activation kinetics of voltage-gated Na (Na_V) and K (K_V) channels. Na_V channels activate much faster and allow the initial Na⁺ influx that generates the depolarizing phase and propagates the signal. Recent experimental results suggest that the molecular basis for this kinetic difference is an amino acid side chain located in the gating pore of the voltage sensor domain, which is a highly conserved isoleucine in K_V channels but an equally highly conserved threonine in Na_V channels. Mutagenesis suggests that the hydrophobicity of this side chain in Shaker K_V channels regulates the energetic barrier that gating charges cross as they move through the gating pore and control the rate of channel opening. We use a multiscale modeling approach to test this hypothesis. We use high-resolution molecular dynamics to study the effect of the mutation on polarization charge within the gating pore. We then incorporate these results in a lower-resolution model of voltage gating to predict the effect of the mutation on the movement of gating charges. The predictions of our hierarchical model are fully consistent with the tested hypothesis, thus suggesting that the faster activation kinetics of Na_V channels comes from a stronger dielectric polarization by threonine (Na_V channel) produced as the first gating charge enters the gating pore compared with isoleucine (K_V channel).

Extending conceptual DFT to include additional variables: oriented external electric field

The extension of the E = E[N, v] functional for exploring chemical reactivity in a conceptual DFT context to include external electric fields is discussed. Concentrating on the case of a homogeneous field the corresponding response functions are identified and integrated, together with the conventional response functions such as permanent dipole moment and polarizability, in an extended response function tree associated with the E = E[N, v, ε] functional. In a case study on the dihalogens F2, Cl2, Br2, I2 the sensitivity of condensed atomic charges (∂q/∂ε) is linked to the polarizability of the halogen atoms. The non-integrated (∂ρ(r)/∂ε) response function, directly related to the field induced density change, is at the basis of these features. It reveals symmetry breaking for a perpendicular field, not detectable in its atom condensed counterpart, and accounts for the induced dipole moment directly related to the molecular polarizability. The much higher sensitivity of the electronic chemical potential/electronegativity as compared to the chemical hardness is highlighted. The response of the condensed Fukui functions to a parallel electric field increases when going down in the periodic table and is interpreted in terms of the extension of the outer contours in the non-condensed Fukui function. In the case of a perpendicular field the (∂f(r)/∂ε) response function hints at stereoselectivity with a preferential side of attack which is not retrieved in its condensed form. In an application the nucleophilic attack on the carbonyl group in H2CO is discussed. Similar to the dihalogens, stereoselectivity is displayed in the Fukui function for nucleophilic attack (f+) in the case of a perpendicular electric field, and opposite to the one that would arise based on the induced density. Disentangling the expression for the evolution of the Fukui function in the presence of an electric field reveals that this difference can be traced back to local differences in the polarization or induced density between the anionic and the neutral system. This difference may be exploited, e.g. for an appropriately substituted H2CO, to generate enantioselectivity.

Recent Progress in the Manipulation of Molecules with DC Electric Fields

The structure and reactivity of a molecule in the condensed phase are governed by its intermolecular interactions with the surrounding environment. The multipole expansion of each molecule in the condensed phase indicates that the intermolecular interactions are essentially electrostatic (e.g., ion-dipole, dipole-dipole, dipole-quadrupole, dipole-induced dipole). The electrostatic field is a fundamental language of intermolecular communications. Therefore, understanding the influence of the electrostatic field on a molecule, that is, the mechanisms by which an electrostatic field manipulates a molecule, from the perspective of molecular structure, energy states, and dynamics is indispensable for illustrating and, by extension, controlling the chemistry in molecular systems.In this Account, we describe the recent progress made in manipulation of molecular processes using an external DC electrostatic field. An electrostatic field with unprecedentedly high strength (≤4 × 10⁸ V/m) was applied in a controlled manner across a molecular film sample using the ice film nanocapacitor method. This field strength is comparable in magnitude to that of weak intermolecular interactions such as van der Waals interactions in the condensed phases. The samples were prepared using a thin film growing technique in vacuum to obtain the desired chemically tailored molecular systems. The examples of prepared systems included small molecules and molecular clusters isolated in cryogenic Ar matrices, frozen molecular films in amorphous or crystalline phase, and interfaces of multilayered molecular films. The response of the molecules to the external field was monitored by reflection-absorption infrared spectroscopy. This approach allowed us to investigate a variety of molecular systems with various intermolecular strength and environments under the influence of strong electrostatic fields. The range of observed molecular behaviors includes the manipulation of molecular orientation, intramolecular dynamics, and proton transfer reactions as an example of stereodynamic control of chemical reactivity. These observations improve our understanding of molecular behaviors in strong electric fields and broaden our perspective on electrostatic manipulation of molecules. This information is also relevant to a variety of research topics in physical and biological sciences where electric fields play a role in molecular and biological functions.

Dynamics and Infrared Spectrocopy of Monomeric and Dimeric Wild Type and Mutant Insulin

The infrared spectroscopy and dynamics of -CO labels in wild type and mutant insulin monomer and dimer are characterized from molecular dynamics simulations using validated force fields. It is found that the spectroscopy of monomeric and dimeric forms in the region of the amide-I vibration differs for residues B24-B26 and D24-D26, which are involved in dimerization of the hormone. Also, the spectroscopic signatures change for mutations at position B24 from phenylalanine, which is conserved in many organisms and is known to play a central role in insulin aggregation, to alanine or glycine. Using three different methods to determine the frequency trajectories (solving the nuclear Schrödinger equation on an effective 1-dimensional potential energy curve, using instantaneous normal modes, and using parametrized frequency maps) leads to the same overall conclusions. The spectroscopic response of monomeric WT and mutant insulin differs from that of their respective dimers, and the spectroscopy of the two monomers in the dimer is also not identical. For the WT and F24A and F24G monomers, spectroscopic shifts are found to be ∼20 cm^-1 for residues (B24-B26) located at the dimerization interface. Although the crystal structure of the dimer is that of a symmetric homodimer, dynamically the two monomers are not equivalent on the nanosecond time scale. Together with earlier work on the thermodynamic stability of the WT and the same mutants, it is concluded that combining computational and experimental infrared spectroscopy provides a potentially powerful way to characterize the aggregation state and dimerization energy of modified insulins.

Role of Brownian Particle Velocity in Bioelectronic Emissions of DNA

Hypothesis: If double stranded DNA (dsDNA) is a charged biomolecule that moves in Earth's magnetic field at a Brownian velocity, then dsDNA may emit bioelectromagnetic waves at energies that reflect discrete genetic states. Methods: This work leverages the Planck-Einstein-de Broglie relationship and applies this concept to Brownian velocity of dsDNA within a cell, to describe the relationship between dsDNA mass, the average Brownian velocity of dsDNA within a cell, and the theoretical wavelengths at which DNA may emit bioelectromagnetic waves. Results: Theoretical emission wavelengths of dsDNA, derived from first principles, were found to correlate closely with experimentally observed emission wavelengths from spectroscopic measurements across various cellular systems in the literature. Conclusion: This work provides a conceptual basis for the potential for unification of bioelectromagnetism with Brownian motion, to elucidate how electromagnetic information can be generated at a subcellular level in biological systems. The implications of how finite mass changes in dsDNA can result in discrete emission wavelengths on electromagnetic timescales is discussed through the lens of genomics. Future refinements of this fundamental methodology may provide a conceptual basis to address previously unexplained multilevel phenomena in the field of biology and is general enough to be extended to other charged biomolecules at a subcellular level. Further exploration in this area could lead to new biological tool development that may augment current genomics methods.