Measuring electric fields and noncovalent interactions using the vibrational stark effect

Probing the structure of water in individual living cells

Water regulates or even governs a wide range of biological processes. Despite its fundamental importance, surprisingly little is known about the structure of intracellular water. Herein we employ a Raman micro-spectroscopy technique to uncover the composition, abundance and vibrational spectra of intracellular water in individual living cells. In three different cell types, we show a small but consistent population (~3%) of non-bulk-like water. It exhibits a weakened hydrogen-bonded network and a more disordered tetrahedral structure. We attribute this population to biointerfacial water located in the vicinity of biomolecules. Moreover, our whole-cell modeling suggests that all soluble (globular) proteins inside cells are surrounded by, on average, one full molecular layer (about 2.6 Angstrom) of biointerfacial water. Furthermore, relative invariance of biointerfacial water is observed among different single cells. Overall, our study not only opens up experimental possibilities of interrogating water structure in vivo but also provides insights into water in life.

Observation of phonon Stark effect

Stark effect, the electric-field analogue of magnetic Zeeman effect, is one of the celebrated phenomena in modern physics and appealing for emergent applications in electronics, optoelectronics, as well as quantum technologies. While in condensed matter it has prospered only for excitons, whether other collective excitations can display Stark effect remains elusive. Here, we report the observation of phonon Stark effect in a two-dimensional quantum system of bilayer 2H-MoS2. The longitudinal acoustic phonon red-shifts linearly with applied electric fields and can be tuned over ~1 THz, evidencing giant Stark effect of phonons. Together with many-body ab initio calculations, we uncover that the observed phonon Stark effect originates fundamentally from the strong coupling between phonons and interlayer excitons (IXs). In addition, IX-mediated electro-phonon intensity modulation up to ~1200% is discovered for infrared-active phonon A2u. Our results unveil the exotic phonon Stark effect and effective phonon engineering by IX-mediated mechanism, promising for a plethora of exciting many-body physics and potential technological innovations.

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.

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.

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.

Cation Dependence of Enniatin B/Membrane-Interactions Assessed Using Surface-Enhanced Infrared Absorption (SEIRA) Spectroscopy

Enniatins are mycotoxins with well-known antibacterial, antifungal, antihelmintic and antiviral activity, which have recently come to attention as potential mitochondriotoxic anticancer agents. The cytotoxicity of enniatins is traced back to ionophoric properties, in which the cyclodepsipeptidic structure results in enniatin:cation-complexes of various stoichiometries proposed as membrane-active species. In this work, we employed a combination of surface-enhanced infrared absorption (SEIRA) spectroscopy, tethered bilayer lipid membranes (tBLMs) and density functional theory (DFT)-based computational spectroscopy to monitor the cation-dependence (Mz+=Na+, K+, Cs+, Li+, Mg2+, Ca2+) on the mechanism of enniatin B (EB) incorporation into membranes and identify the functionally relevant EBn : Mz+ complexes formed. We find that Na+ promotes a cooperative incorporation, modelled via an autocatalytic mechanism and mediated by a distorted 2 : 1-EB2 : Na+ complex. K+ (and Cs+) leads to a direct but less efficient insertion into membranes due to the adoption of "ideal" EB2 : K+ sandwich complexes. In contrast, the presence of Li+, Mg2+, and Ca2+ causes a (partial) extraction of EB from the membrane via the formation of "belted" 1 : 1-EB : Mz+ complexes, which screen the cationic charge less efficiently. Our results point to a relevance of the cation dependence for the transport into the malignant cells where the mitochondriotoxic anticancer activity is exerted.

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.

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.

Electrostatics as a Guiding Principle in Understanding and Designing Enzymes

Enzyme design faces challenges related to the implementation of the basic principles that govern the catalytic activity in natural enzymes. In this work, we revisit basic electrostatic concepts that have been shown to explain the origin of enzymatic efficiency like preorganization and reorganization. Using magnitudes such as the electrostatic potential and the electric field generated by the protein, we explain how these concepts work in different enzymes and how they can be used to rationalize the consequences of point mutations. We also discuss examples of protein design in which electrostatic effects have been implemented. For the near future, molecular simulations, coupled with the use of machine learning methods, can be used to implement electrostatics as a guiding principle for enzyme designs.

Influence of electric double layer rigidity on CO adsorption and electroreduction rate

Understanding the structure of the electric double layer (EDL) is critical for designing efficient electrocatalytic processes. However, the interplay between reactant adsorbates and the concentrated ionic species within the EDL remains an aspect that has yet to be fully explored. In the present study, we employ electrochemical CO reduction on Cu as a model reaction to reveal the significant impact of EDL structure on CO adsorption. By altering the sequence of applying negative potential and elevating CO pressure, we discern two distinct EDL structures with varying cation density and CO coverage. Our findings demonstrate that the EDL comprising densely packed cations substantially hinders CO adsorption on the Cu as opposed to the EDL containing less compact cations. These two different EDL structures remained stable over the course of our experiments, despite their identical initial and final conditions, suggesting an insurmountable kinetic barrier present in between. Moreover, we show that the size and identity of cations play decisive roles in determining the properties of the EDL in CO electroreduction on Cu. This study presents a refined adaptation of the classical Gouy-Chapman-Stern model and highlights its catalytic importance, which bridges the mechanistic gap between the EDL structure and cathodic reactions.

On the Role of a Conserved Tryptophan in the Chromophore Pocket of Cyanobacteriochrome

The cyanobacteriochrome Slr1393 can be photoconverted between a red (Pr) and green absorbing form (Pg). The recently determined crystal structures of both states suggest a major movement of Trp496 from a stacking interaction with ring D of the phycocyanobilin (PCB) chromophore in Pr to a position outside the chromophore pocket in Pg. Here, we investigated the role of this amino acid during photoconversion in solution using engineered protein variants in which Trp496 was substituted by natural and non-natural amino acids. These variants and the native protein were studied by various spectroscopic techniques (UV-vis absorption, fluorescence, IR, NIR and UV resonance Raman) complemented by theoretical approaches. Trp496 is shown to affect the electronic transition of PCB and to be essential for the thermal equilibrium between Pr and an intermediate state O600. However, Trp496 is not required to stabilize the tilted orientation of ring D in Pr, and does not play a role in the secondary structure changes of Slr1393 during the Pr/Pg transition. The present results confirm the re-orientation of Trp496 upon Pr → Pg conversion, but do not provide evidence of a major change in the microenvironment of this residue. Structural models indicate the penetration of water molecules into the chromophore pocket in both Pr and Pg states and thus water-Trp contacts, which can readily account for the subtle spectral changes between Pr and Pg. Thus, we conclude that reorientation of Trp496 during the Pr-to-Pg photoconversion in solution is not associated with a major change in the dielectric environment in the two states.

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.

Studying the cation dependence of CO2 reduction intermediates at Cu by in situ VSFG spectroscopy

The nature of the electrolyte cation is known to have a significant impact on electrochemical reduction of CO2 at catalyst|electrolyte interfaces. An understanding of the underlying mechanism responsible for catalytic enhancement as the alkali metal cation group is descended is key to guide catalyst development. Here, we use in situ vibrational sum frequency generation (VSFG) spectroscopy to monitor changes in the binding modes of the CO intermediate at the electrochemical interface of a polycrystalline Cu electrode during CO2 reduction as the electrolyte cation is varied. A CObridge mode is observed only when using Cs+, a cation that is known to facilitate CO2 reduction on Cu, supporting the proposed involvement of CObridge sites in CO coupling mechanisms during CO2 reduction. Ex situ measurements show that the cation dependent CObridge modes correlate with morphological changes of the Cu surface.

Voltage-Induced Inversion of Band Bending and Photovoltages at Semiconductor/Liquid Interfaces

At semiconductor/liquid interfaces, the surface potential and photovoltages are produced by a combination of band bending and quasi-Fermi-level splitting at the semiconductor surface, which are usually treated in a qualitative fashion. As such, it is important to develop quantitative metrics for the band energies and photovoltaics at these interfaces. Here, we present a spectroscopic method for monitoring the photovoltages produced at semiconductor/liquid junctions. The surface reporter molecule mercaptobenzonitrile (MBN) is functionalized on the photoelectrode surface (p-type silicon) and is measured using in situ surface-enhanced Raman scattering (SERS) spectroscopy with a water immersion lens under electrochemical working conditions. In particular, the vibrational frequency of the C≡N stretch mode (ωCN) around 2225 cm-1 is sensitive to the local electric field in solution at the electrode/electrolyte interface via the vibrational Stark effect. Over the applied potential range from -0.8 to 0.6 V vs Ag/AgCl, we observe ωCN to increase from 2220 to 2229 cm-1 (at low laser power). As the incident laser power is increased from 83.5 μW to 13.3 mW, we observe additional shifts of ΔωCN = ±1 cm-1, corresponding to photovoltages produced at the semiconductor/liquid interface ΔV = ±0.2 V. Based on Mott-Schottky measurements, the flat band potential (FBP) occurs at -0.39 V vs Ag/AgCl. For applied potentials above the FBP, we observe ΔωCN > 0 (i.e., blue-shifts ∼1 cm-1) corresponding to positive photovoltages, whereas for applied potentials below the flat band potential, we observe ΔωCN < 0 (i.e., red-shifts ∼1 cm-1) corresponding to negative photovoltages. These spectroscopic observations reveal voltage-induced changes in the band bending at the semiconductor/liquid junction that, thus far, have been difficult to measure.

Vibrational Coupling Infrared Nanocrystallography

Coupling between molecular vibrations leads to collective vibrational states with spectral features sensitive to local molecular order. This provides spectroscopic access to the low-frequency intermolecular energy landscape. In its nanospectroscopic implementation, this technique of vibrational coupling nanocrystallography (VCNC) offers information on molecular disorder and domain formation with nanometer spatial resolution. However, deriving local molecular order relies on prior knowledge of the transition dipole magnitude and crystal structure of the underlying ordered phase. Here we develop a quantitative model for VCNC by relating nano-FTIR collective vibrational spectra to the molecular crystal structure from X-ray crystallography. We experimentally validate our approach at the example of a metal organic porphyrin complex with a carbonyl ligand as the probe vibration. This framework establishes VCNC as a powerful tool for measuring low-energy molecular interactions, wave function delocalization, nanoscale disorder, and domain formation in a wide range of molecular systems.

Nitrile Substituents at the Conjugated Dipyridophenazine Moiety as Infrared Redox Markers in Electrochemically Reduced Heteroleptic Ru(II) Polypyridyl Complexes

Ruthenium(II) complexes [Ru(tap)2(NN)]2+ (tap = 1,4,5,8-tetraazaphenanthrene, NN = 11-cyano-dipyrido[3,2-a:2',3'-c]phenazine (11-CN-dppz) and 11,12-dicyano-dipyrido[3,2-a:2',3'-c]phenazine (11,12-CN-dppz)) feature the C≡N groups as infrared (IR)-active redox markers. They were studied by cyclic voltammetry, UV-vis, and IR spectroelectrochemistry (SEC), and density functional theory calculations to assign the four 1e- reduction waves R1-R4 observed in dichloromethane. Generally, the NN ligands are reduced first (R1). For [Ru(tap)2(11,12-CN-dppz)]2+, R1 is sufficiently separated from R2 and delocalized over both tap ligands. Accordingly, IR SEC conducted at R1 shows a large red shift of the νs,as(C≡N) modes by -18/-28 cm-1, accompanied by a 4-fold enhancement of the νs(C≡N) intensity, comparably with reference data for free 11,12-CN-dppz. The first tap-based reduction of spin-doublet [Ru(tap)2(11,12-CN-dppz)]+ to spin-triplet [Ru(tap)2(11,12-CN-dppz)] at R2 decreased ν(C≡N) by merely -2 cm-1, while the intensity enhancement reached an overall factor of 8. Comparably, a red shift of ν(C≡N) by -27 cm-1 resulted from the 1e- reduction of [Ru(tap)2(11-CN-dppz)]2+ at R1 (poorly resolved from R2), and the intensity enhancement was roughly 3-fold. Concomitant 1e- reductions of the tap ligands (R2 and R3) caused only minor ν(C≡N) shifts of -3 cm-1 and increased the absorbance by overall factors of 6.5 and 8, respectively.

Enticing a Proton using Single Ammonia Molecule as Bait

In microhydrated acid-solvent clusters, deprotonation of an acid is assisted by a critical number of solvent molecules and a solvent electric field. Born-Oppenheimer molecular dynamics simulations reveal that trifluoroacetic acid undergoes spontaneous proton transfer in water clusters, with the critical number being five. Acetic acid and phenol, on the other hand, do not dissociate even in the presence of a large number of water molecules (in excess of 40). The addition of a single ammonia molecule to the water cluster, which interacts directly with the protic group, lowers the critical number of solvent water molecules required for proton transfer to three and seven in the case of acetic acid and phenol, respectively. The population of the undissociated and the proton-transferred structures get dispersed to form separate islands on the electric field versus the O-H distance representation with the cusp representing the critical values. The critical electric fields for the spontaneous proton transfer are around 254, 237, and 318 MV cm-1 for trifluoroacetic acid, acetic acid, and phenol, respectively. In the case of phenol, the free energy profiles suggest that proton transfer to the ammonia moiety embedded in water promotes proton transfer efficiently due to the higher basicity of ammonia and enhanced hydrogen bonding network of solvent water, vis-à-vis phenol-ammonia clusters.

Mapping enzyme activity in living systems by real-time mid-infrared photothermal imaging of nitrile chameleons

Simultaneous spatial mapping of the activity of multiple enzymes in a living system can elucidate their functions in health and disease. However, methods based on monitoring fluorescent substrates are limited. Here, we report the development of nitrile (C≡N)-tagged enzyme activity reporters, named nitrile chameleons, for the peak shift between substrate and product. To image these reporters in real time, we developed a laser-scanning mid-infrared photothermal imaging system capable of imaging the enzymatic substrates and products at a resolution of 300 nm. We show that when combined, these tools can map the activity distribution of different enzymes and measure their relative catalytic efficiency in living systems such as cancer cells, Caenorhabditis elegans, and brain tissues, and can be used to directly visualize caspase-phosphatase interactions during apoptosis. Our method is generally applicable to a broad category of enzymes and will enable new analyses of enzymes in their native context.

In Situ Spectroscopic Elucidation of the Electrochemical Potential Drop at Polyelectrolytes/Au Interfaces

Polyelectrolytes have been widely applied in electrochemical devices. Understanding the polyelectrolyte/electrode interfaces is pivotal for polyelectrolyte-based applications. Here, we measured the electrochemical potential drop and the local activity of the mobile ion of H+ or OH- at the polyelectrolytes/Au interfaces by in situ electrochemical surface-enhanced Raman spectroscopy and voltammetry in three-electrode cells. We found that the potential dependences of the electrochemical potential drop in polyelectrolytes were smaller than that in conventional electrolyte solutions. The interfacial activity of H+ or OH- was much lower than that of bulk polyelectrolytes. The potential-dependent molecular dynamics simulations showed that the mobility of ionomers of polyelectrolytes in an electrostatic field was limited by a polymer matrix. These results suggested a characteristically thicker compact layer in the electrical double layer of a polyelectrolyte/electrode interface due to the accumulation of mobile H+ or OH- with a thicker hydration layer and immobile ionomers.

Triple-Bond Vibrations: Emerging Applications in Energy and Biological Sciences

Triple bonds, such as that formed between two carbon atoms (i.e., C≡C) or that formed between one carbon atom and one nitrogen atom (i.e., C≡N), afford unique chemical bonding and hence vibrational characteristics. As such, they are not only frequently used to construct molecules with tailored chemical and/or physical properties but also employed as vibrational probes to provide site-specific chemical and/or physical information at the molecular level. Herein, we offer our perspective on the emerging applications of various triple-bond vibrations in energy and biological sciences with a focus on C≡C and C≡N triple bonds.

Electrostatic Impact of Brefeldin A on Thiocyanate Probes Surrounding the Interface of Arf1-BFA-ARNO4M, a Protein-Drug-Protein Complex

Protein-protein interactions regulate many cellular processes, making them ideal drug candidates. Design of such drugs, however, is hindered by a lack of understanding of the factors that contribute to the interaction specificity. Specific protein-protein complexes possess both structural and electrostatic complementarity, and while structural complementarity of protein complexes has been extensively investigated, fundamental understanding of the complicated networks of electrostatic interactions at these interfaces is lacking, thus hindering the rational design of orthosterically binding small molecules. To better understand the electrostatic interactions at protein interfaces and how a small molecule could contribute to and fit within that environment, we used a model protein-drug-protein system, Arf1-BFA-ARNO4M, to investigate how small molecule brefeldin A (BFA) perturbs the Arf1-ARNO4M interface. By using nitrile probe labeled Arf1 sites and measuring vibrational Stark effects as well as temperature dependent infrared shifts, we measured changes in the electric field and hydrogen bonding at this interface upon BFA binding. At all five probe locations of Arf1, we found that the vibrational shifts resulting from BFA binding corroborate trends found in Poisson-Boltzmann calculations of surface potentials of Arf1-ARNO4M and Arf1-BFA-ARNO4M, where BFA contributes negative electrostatic potential to the protein interface. The data also corroborate previous hypotheses about the mechanism of interfacial binding and confirm that alternating patches of hydrophobic and polar interactions lead to BFA binding specificity. These findings demonstrate the impact of BFA on this protein-protein interface and have implications for the design of other interfacial drug candidates.

Substrate Positioning Dynamics Involves a Non-Electrostatic Component to Mediate Catalysis

Substrate positioning dynamics (SPD) orients the substrate in the active site, thereby influencing catalytic efficiency. However, it remains unknown whether SPD effects originate primarily from electrostatic perturbation inside the enzyme or can independently mediate catalysis with a significant non-electrostatic component. In this work, we investigated how the non-electrostatic component of SPD affects transition state (TS) stabilization. Using high-throughput enzyme modeling, we selected Kemp eliminase variants with similar electrostatics inside the enzyme but significantly different SPD. The kinetic parameters of these mutants were experimentally characterized. We observed a valley-shaped, two-segment linear correlation between the TS stabilization free energy (converted from kinetic parameters) and substrate positioning index (a metric to quantify SPD). The energy varies by approximately 2 kcal/mol. Favorable SPD was observed for the distal mutant R154W, increasing the proportion of reactive conformations and leading to the lowest activation free energy. These results indicate the substantial contribution of the non-electrostatic component of SPD to enzyme catalytic efficiency.

Vibrational Probe at the Electrochemical Interface: Dependence on Plasmon Coupling and Potential of the Lineshape in Two-Dimensional Infrared Spectroscopy

Two-dimensional infrared spectroscopy of vibrational probes at an electrode surface shows promise for studying the structural dynamics at an active electrochemical interface. This interface is a complex environment where the solution structures in response to the applied potential. A strategy for achieving the necessary monolayer sensitivity is to use a plasmonically active electrode, which enhances the electromagnetic fields that produce the spectroscopic response. Here, we show how the coupling between the plasmon and the vibrations of the molecular monolayer impacts the FTIR and 2D IR spectroscopy, with an emphasis on the electrochemical potential difference spectra. We show how mixing between the vibrational and plasmonic states gives rise to the distortions that are observed in these measurements. This provides an important step toward 2D IR measurements of vibrational probes at the electrochemical interface as a tool for probing the structural dynamics in the double layer.

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.

The catalytic reaction of cytochrome c oxidase probed by in situ gas titrations and FTIR difference spectroscopy

Cytochrome c oxidase (CcO) is a transmembrane heme‑copper metalloenzyme that catalyzes the reduction of O2 to H2O at the reducing end of the respiratory electron transport chain. To understand this reaction, we followed the conversion of CcO from Rhodobacter sphaeroides between several active-ready and carbon monoxide-inhibited states via attenuated total reflection Fourier-transform infrared (ATR FTIR) difference spectroscopy. Utilizing a novel gas titration setup, we prepared the mixed-valence, CO-inhibited R2CO state as well as the fully-reduced R4 and R4CO states and induced the "active ready" oxidized state OH. These experiments are performed in the dark yielding FTIR difference spectra exclusively triggered by exposure to O2, the natural substrate of CcO. Our data demonstrate that the presence of CO at heme a3 does not impair the catalytic oxidation of CcO when the cycle starts from the fully-reduced states. Interestingly, when starting from the R2CO state, the release of the CO ligand upon purging with inert gas yield a product that is indistinguishable from photolysis-induced states. The observed changes at heme a3 in the catalytic binuclear center (BNC) result from the loss of CO and are unrelated to electronic excitation upon illumination. Based on our experiments, we re-evaluate the assignment of marker bands that appear in time-resolved photolysis and perfusion-induced experiments on CcO.

Electro-inductive Effect Dominates Vibrational Frequency Shifts of Conjugated Probes on Gold Electrodes

Interfacial electric fields play a critical role in electrocatalysis and are often characterized by using vibrational probes attached to an electrode surface. Understanding the physical principles dictating the impact of the applied electrode potential on the vibrational probe frequency is important. Herein, a comparative study is performed for two molecular probes attached to a gold electrode. Both probes contain a nitrile (CN) group, but 4-mercaptobenzonitrile (4-MBN) exhibits continuous conjugation from the electrode through the nitrile group, whereas this conjugation is interrupted for 2-(4-mercaptophenyl)acetonitrile (4-MPCN). Periodic density functional theory calculations predict that the CN vibrational frequency shift of the 4-MBN system is dominated by induction, which is a through-bond polarization effect, leading to a strong potential dependence that does not depend significantly on the orientation of the CN bond relative to the surface. In contrast, the CN vibrational frequency shift of the 4-MPCN system is influenced less by induction and more by through-space electric field effects, leading to a weaker potential dependence and a greater orientation dependence. These theoretical predictions were confirmed by surface-enhanced Raman spectroscopy experiments. Balancing through-bond and through-space electrostatic effects may assist in the fundamental understanding and design of electrocatalytic systems.

Comment on "Applicability of perturbed matrix method for charge transfer studies at bio/metallic interfaces: a case of azurin" by O. Kontkanen, D. Biriukov and Z. Futera, Phys. Chem. Chem. Phys., 2023, , 12479

Polarizability is a fundamental property of all molecular systems describing the deformation of the molecular electronic density in response to an applied electric field. The question of whether polarizability of the active site needs to be included in theories of enzymatic activity remains open. Hybrid quantum mechanical/molecular mechanical calculations are hampered by difficulties faced by many quantum-chemistry algorithms to provide sufficiently accurate estimates of the anisotropic second-rank tensor of molecular polarizability. In this Comment, we provide general theoretical arguments for the values of polarizability of the quantum region or a molecule which have to be reproduced by electronic structure calculations.

Good Vibrations Report on the DNA Quadruplex Binding of an Excited State Amplified Ruthenium Polypyridyl IR Probe

The nitrile containing Ru(II)polypyridyl complex [Ru(phen)2(11,12-dCN-dppz)]2+ (1) is shown to act as a sensitive infrared probe of G-quadruplex (G4) structures. UV-visible absorption spectroscopy reveals enantiomer sensitive binding for the hybrid htel(K) and antiparallel htel(Na) G4s formed by the human telomer sequence d[AG3(TTAG3)3]. Time-resolved infrared (TRIR) of 1 upon 400 nm excitation indicates dominant interactions with the guanine bases in the case of Λ-1/htel(K), Δ-1/htel(K), and Λ-1/htel(Na) binding, whereas Δ-1/htel(Na) binding is associated with interactions with thymine and adenine bases in the loop. The intense nitrile transient at 2232 cm-1 undergoes a linear shift to lower frequency as the solution hydrogen bonding environment decreases in DMSO/water mixtures. This shift is used as a sensitive reporter of the nitrile environment within the binding pocket. The lifetime of 1 in D2O (ca. 100 ps) is found to increase upon DNA binding, and monitoring of the nitrile and ligand transients as well as the diagnostic DNA bleach bands shows that this increase is related to greater protection from the solvent environment. Molecular dynamics simulations together with binding energy calculations identify the most favorable binding site for each system, which are in excellent agreement with the observed TRIR solution study. This study shows the power of combining the environmental sensitivity of an infrared (IR) probe in its excited state with the TRIR DNA "site effect" to gain important information about the binding site of photoactive agents and points to the potential of such amplified IR probes as sensitive reporters of biological environments.

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.

Phenol as a Tethering Group to Gold Surfaces: Stark Response and Comparison to Benzenethiol

Understanding the adsorption of organic molecules on metals is important in numerous areas of surface science, including electrocatalysis, electrosynthesis, and biosensing. While thiols are commonly used to tether organic molecules on metals, it is desirable to broaden the range of anchoring groups. In this study, we use a combined spectroelectrochemical and computational approach to demonstrate the adsorption of 4-cyanophenols (CPs) on polycrystalline gold. Using the nitrile stretching vibration as a marker, we confirm the adsorption of CP on the gold electrode and compare our results with those obtained for the thiol counterpart, 4-mercaptobenzonitirle (MBN). Our results reveal that CP adsorbs on the gold electrode via the OH linker, as evidenced by the similarity in the direction and magnitude of the nitrite Stark shifts for CP and MBN. This finding paves the way for exploring new approaches to modify electrode surfaces for controlled reactivity. Furthermore, it highlights adsorption on metals as an important step in the electroreactivity of phenols.

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.

EnzyHTP Computational Directed Evolution with Adaptive Resource Allocation

Directed evolution facilitates enzyme engineering via iterative rounds of mutagenesis. Despite the wide applications of high-throughput screening, building "smart libraries" to effectively identify beneficial variants remains a major challenge in the community. Here, we developed a new computational directed evolution protocol based on EnzyHTP, a software that we have previously reported to automate enzyme modeling. To enhance the throughput efficiency, we implemented an adaptive resource allocation strategy that dynamically allocates different types of computing resources (e.g., GPU/CPU) based on the specific need of an enzyme modeling subtask in the workflow. We implemented the strategy as a Python library and tested the library using fluoroacetate dehalogenase as a model enzyme. The results show that compared to fixed resource allocation where both CPU and GPU are on-call for use during the entire workflow, applying adaptive resource allocation can save 87% CPU hours and 14% GPU hours. Furthermore, we constructed a computational directed evolution protocol under the framework of adaptive resource allocation. The workflow was tested against two rounds of mutational screening in the directed evolution experiments of Kemp eliminase (KE07) with a total of 184 mutants. Using folding stability and electrostatic stabilization energy as computational readout, we identified all four experimentally observed target variants. Enabled by the workflow, the entire computation task (i.e., 18.4 μs MD and 18,400 QM single-point calculations) completes in 3 days of wall-clock time using ∼30 GPUs and ∼1000 CPUs.

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

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

C≡N Stretching Frequency as a Convenient Reporter of Charge Separation in Molecular Systems

Previously, several studies have shown that, for a set of structurally related nitrile compounds, there could be a linear relationship between the total charge on the nitrile group (qCN) and its stretching frequency (νCN). However, it is unclear whether the corresponding frequency and charge properties of structurally different nitrile compounds can be described by a single linear νCN-qCN relationship. Herein, we compute the qCN magnitudes of a large number of nitrile-containing molecules whose νCN values cover a spectral range of ca. 200 cm-1 and are measured under different experimental conditions. Our results reveal that νCN indeed exhibits a linear dependence on qCN, with a slope of 637 ± 30 cm-1/charge. Because the nitrile moiety is a commonly used building block in electronic donor-acceptor (D-A) molecular systems, we believe that this linear relationship will find utility in a wide range of applications where such D-A constructs are used, such as in organic photovoltaic assemblies. In addition, we apply this linear relationship to characterize the degree of charge transfer upon photoexcitation of two indole derivatives, 5-cyanoindole and 6-cyanoindole, and are able to show that in both cases, the fluorescence emission arises from a charge-transfer or La state.

Structural Characterization of Fluorescent Proteins Using Tunable Femtosecond Stimulated Raman Spectroscopy

The versatile functions of fluorescent proteins (FPs) as fluorescence biomarkers depend on their intrinsic chromophores interacting with the protein environment. Besides X-ray crystallography, vibrational spectroscopy represents a highly valuable tool for characterizing the chromophore structure and revealing the roles of chromophore-environment interactions. In this work, we aim to benchmark the ground-state vibrational signatures of a series of FPs with emission colors spanning from green, yellow, orange, to red, as well as the solvated model chromophores for some of these FPs, using wavelength-tunable femtosecond stimulated Raman spectroscopy (FSRS) in conjunction with quantum calculations. We systematically analyzed and discussed four factors underlying the vibrational properties of FP chromophores: sidechain structure, conjugation structure, chromophore conformation, and the protein environment. A prominent bond-stretching mode characteristic of the quinoidal resonance structure is found to be conserved in most FPs and model chromophores investigated, which can be used as a vibrational marker to interpret chromophore-environment interactions and structural effects on the electronic properties of the chromophore. The fundamental insights gained for these light-sensing units (e.g., protein active sites) substantiate the unique and powerful capability of wavelength-tunable FSRS in delineating FP chromophore properties with high sensitivity and resolution in solution and protein matrices. The comprehensive characterization for various FPs across a colorful palette could also serve as a solid foundation for future spectroscopic studies and the rational engineering of FPs with diverse and improved functions.

AMOEBA Force Field Predicts Accurate Hydrogen Bond Counts of Nitriles in SNase by Revealing Water-Protein Interaction in Vibrational Absorption Frequencies

Precisely quantifying the magnitude and direction of electric fields in proteins has long been an outstanding challenge in understanding biological functions. Nitrile vibrational Stark effect probes have been shown to be minimally disruptive to the protein structure and can be better direct reporters of local electrostatic field in the native state of a protein than other measures such as pKa shifts of titratable residues. However, interpretations of the connection between measured vibrational energy and electric field rely on the accurate molecular understanding of interactions of the nitrile group and its environment, particularly from hydrogen bonding. In this work, we compared the extent of hydrogen bonding calculated in two common force fields, the fixed charge force field Amber03 and polarizable force field AMOEBA, at 10 locations of cyanocysteine (CNC) in staphylococcal nuclease (SNase) against the experimental nitrile absorption frequency in terms of full width at half-maximum (FWHM) and frequency temperature line slope (FTLS). We observed that the number of hydrogen bonds correlated well in AMOEBA trajectories with respect to both the FWHM (r = 0.88) and the FTLS (r = -0.85), whereas the correlation of Amber03 trajectories was less reliable because the Amber03 force field predicted more hydrogen bonds in some mutants. Moreover, we demonstrated that contributions from the interactions between CNC and nearby water molecules were significant in AMOEBA trajectories but were not predicted by Amber03. We conclude that although the nitrile absorption peak shape could be qualitatively predicted by the fixed charge Amber03 force field, the detailed electrostatic environment measured by the nitrile probe in terms of the extent of hydrogen bonding could only be accurately observed in the AMOEBA trajectories, where the permanent dipole, quadrupole, and dipole-induced-dipole polarizable interactions were all taken into account. The significance of this finding to the goal of accurately predicting electric fields in complex biomolecular environments is discussed.

Vibrational Spectroscopy of Phytochromes

Phytochromes are biological photoswitches that translate light into physiological functions. Spectroscopic techniques are essential tools for molecular research into these photoreceptors. This review is directed at summarizing how resonance Raman and IR spectroscopy contributed to an understanding of the structure, dynamics, and reaction mechanism of phytochromes, outlining the substantial experimental and theoretical challenges and describing the strategies to master them. It is shown that the potential of the various vibrational spectroscopic techniques can be most efficiently exploited using integral approaches via a combination of theoretical methods as well as other experimental techniques.

Charge and Solvent Effects on the Redox Behavior of Vanadyl Salen-Crown Complexes

The incorporation of charged groups proximal to a redox active transition metal center can impact the local electric field, altering redox behavior and enhancing catalysis. Vanadyl salen (salen = N,N'-ethylenebis(salicylideneaminato)) complexes functionalized with a crown ether containing a nonredox active metal cation (V-Na, V-K, V-Ba, V-La, V-Ce, and V-Nd) were synthesized. The electrochemical behavior of this series of complexes was investigated by cyclic voltammetry in solvents with varying polarity and dielectric constant (ε) (acetonitrile, ε = 37.5; N,N-dimethylformamide, ε = 36.7; and dichloromethane, ε = 8.93). The vanadium(V/IV) reduction potential shifted anodically with increasing cation charge compared to a complex lacking a proximal cation (ΔE_1/2 > 900 mV in acetonitrile and >700 mV in dichloromethane). In contrast, the reduction potential for all vanadyl salen-crown complexes measured in N,N-dimethylformamide was insensitive to the magnitude of the cationic charge, regardless of the electrolyte or counteranion used. Titration studies of N,N-dimethylformamide into acetonitrile resulted in cathodic shifting of the vanadium(V/IV) reduction potential with increasing concentration of N,N-dimethylformamide. Binding constants of N,N-dimethylformamide (log(K_DMF)) for the series of crown complexes show increased binding affinity in the order of V-La > V-Ba > V-K > (salen)V(O), indicating an enhancement of Lewis acid/base interaction with increasing cationic charge. The redox behavior of (salen)V(O) and (salen-OMe)V(O) (salen-OMe = N,N'-ethylenebis(3-methoxysalicylideneamine) was also investigated and compared to the crown-containing complexes. For (salen-OMe)V(O), a weak association of triflate salt at the vanadium(IV) oxidation state was observed through cyclic voltammetry titration experiments, and cation dissociation upon oxidation to vanadium(V) was identified. These studies demonstrate the noninnocent role of solvent coordination and cation/anion effects on redox behavior and, by extension, the local electric field.

Electric-field induced entropic effects in liquid water

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

Evaluating the Transition State Stabilization/Destabilization Effects of the Electric Fields from Scaffold Residues by a QM/MM Approach

The protein scaffolds of enzymes not only provide structural support for the catalytic center but also exert preorganized electric fields for electrostatic catalysis. In recent years, uniform oriented external electric fields (OEEFs) have been widely applied to enzymatic reactions to mimic the electrostatic effects of the environment. However, the electric fields exerted by individual residues in proteins may be quite heterogeneous across the active site, with varying directions and strengths at different positions of the active site. Here, we propose a QM/MM-based approach to evaluate the effects of the electric fields exerted by individual residues in the protein scaffold. In particular, the heterogeneity of the residue electric fields and the effect of the native protein environment can be properly accounted for by this QM/MM approach. A case study of the O-O heterolysis reaction in the catalytic cycle of TyrH shows that (1) for scaffold residues that are relatively far from the active site, the heterogeneity of the residue electric field in the active site is not very significant and the electrostatic stabilization/destabilization due to each residue can be well approximated with the interaction energy between a uniform electric field and the QM region dipole; (2) for scaffold residues near the active site, the residue electric fields can be highly heterogeneous along the breaking O-O bond. In such a case, approximating the residue electric fields as uniform fields may misrepresent the overall electrostatic effect of the residue. The present QM/MM approach can be applied to evaluate the residues' electrostatic impact on enzymatic reactions, which also can be useful in computational optimization of electric fields to boost the enzyme catalysis.

Vibrational Solvatochromism Study of the C-H···O Improper Hydrogen Bond

The improper C-H···O hydrogen bond is an important weak interaction, with broad implications for protein and nucleic acid structure, molecular recognition, enzyme catalysis, and drug interaction. Despite its wide identification in crystal structures, the general existence of C-H···O hydrogen bonds remains elusive especially for natural C-H groups in bulk aqueous solutions at room temperature. Vibrational spectroscopy is a promising methodology to tackle this challenge, as formation of C-H···O hydrogen bonds usually causes shifts of the C-H stretch frequency. Yet, prior observations are inconclusive, as they are all based on a simple blue-shift in aqueous solution and cannot distinguish if it is an effect caused by solvent reorganization or a specific hydrogen-bonding interaction. In this work, we used vibrational solvatochromism as a calibration of the solvent reorganization effect and identified a specific H-bonding interaction. We performed vibrational solvatochromism study of C-H(D) of multiple alcohol molecules including the CH mode of CD₃CH(OH)CD₃ and the CD₃ modes of CD₃OH, CD₃CH₂OH, and CD₃CH(OH)CD₃ in a series of solvents. We found an abnormal blue-shift of the Raman frequency of the C-H and C-D bonds at both the C_α and C_β positions of alcohols in water, which lies in an opposite direction to the expected trend due to vibrational solvatochromism. This experimental evidence supports that the improper C-H···O hydrogen bonds might generally exist between nonpolarized C-H and water in liquid solutions at room temperature.

Tools for analyzing protonation states and for tracing proton transfer pathways with examples from the Rb. sphaeroides photosynthetic reaction centers

Protons participate in many reactions. In proteins, protons need paths to move in and out of buried active sites. The vectorial movement of protons coupled to electron transfer reactions establishes the transmembrane electrochemical gradient used for many reactions, including ATP synthesis. Protons move through hydrogen bonded chains of waters and hydroxy side chains via the Grotthuss mechanism and by proton binding and release from acidic and basic residues. MCCE analysis shows that proteins exist in a large number of protonation states. Knowledge of the equilibrium ensemble can provide a rational basis for setting protonation states in simulations that fix them, such as molecular dynamics (MD). The proton path into the Q_B site in the bacterial reaction centers (RCs) of Rb. sphaeroides is analyzed by MD to provide an example of the benefits of using protonation states found by the MCCE program. A tangled web of side chains and waters link the cytoplasm to Q_B. MCCE analysis of snapshots from multiple trajectories shows that changing the input protonation state of a residue in MD biases the trajectory shifting the proton affinity of that residue. However, the proton affinity of some residues is more sensitive to the input structure. The proton transfer networks derived from different trajectories are quite robust. There are some changes in connectivity that are largely restricted to the specific residues whose protonation state is changed. Trajectories with Q_B^•- are compared with earlier results obtained with Q_B [Wei et. al Photosynthesis Research volume 152, pages153-165 (2022)] showing only modest changes. While introducing new methods the study highlights the difficulty of establishing the connections between protein conformation.

Distinguishing between the Electrostatic Effects and Explicit Ion Interactions in a Stark Probe

Vibrational Stark probes are incisive tools for measuring local electric fields in a wide range of chemical environments. The interpretation of the frequency shift often gets complicated due to the specific interactions of the probe, such as hydrogen bonding and Lewis bonding. Therefore, it is important to distinguish between the pure electrostatic response and the response due to such specific interactions. Here we report a molecular system that is sensitive to both the Stark effect from a single ion and the explicit Lewis bonding of ions with the probe. The molecule consists of a crown ether with an appended benzonitrile. The crown captures cations of various charges, and the electric field from the ions is sensed by the benzonitrile probe. Additionally, the lone pair of the benzonitrile can engage in Lewis interactions with some of the ions by donating partial charge density to the ions. Our system exhibits both of these effects and therefore is a suitable test bed for distinguishing between the pure electrostatic and the Lewis interactions. Our computational results show that the electrostatic influence of the ion is operative at large distances, while the Lewis interaction becomes important only within distances that permit orbital overlap. Our results may be useful for using the nitrile probe for measuring electrostatic and coordination effects in complex ionic environments such as the electrode-electrolyte interfaces.

Forced Interactions: Ionic Polymers at Charged Surfactant Interfaces

Characterizing electrostatic interactions at heterogeneous interfaces is critical for developing a fundamental description of the dynamic processes at charged interfaces. Water-in-oil reverse micelles (RMs) offer a high degree of tunability across composition, polarity, and temperature, making them ideal systems for studying interactions at heterogeneous liquid-liquid interfaces. In the present study, we use a combination of ultrafast two-dimensional infrared spectroscopy and molecular dynamics (MD) simulations to determine the picosecond interfacial dynamics in RMs containing binary compositions of sorbitan monostearate and anionic or cationic cosurfactants, which are used to tune the ratio of charged to nonionic surfactants at the interface. The positively charged polyethylenimine (PEI) polymer is encapsulated within the RMs, and the carbonyl stretching mode of sorbitan monostearate reports on the interfacial hydrogen-bond populations and dynamics. The results show that hydrogen-bond populations are altered through the inclusion of both negatively and positively charged cosurfactants. Charged surfactants increase interfacial water penetration into the surfactant layer, and the surface localization of polymers decreases water penetration. Local hydrogen-bond dynamics undergo a slowdown with the inclusion of charged surfactants, and the encapsulation of polymers results in similar effects, irrespective of the charge.

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.

Quantification of the Real Plasmonic Field Transverse Distribution in a Nanocavity Using the Vibrational Stark Effect

Quantifying the real plasmonic field strength experimentally has been long pursued in expanding the applications related to plasmonic enhancement. However, it is still an enormous challenge to determine the inhomogeneous plasmonic field distribution. Here, self-assembled monolayers (SAMs) of 4-mercaptobenzonitrile (MBN) are sandwiched as a gap spacer in a nanoparticle-on-mirror (NPoM) structure, effectively forming ultrahigh field enhancement to observe Stark shifts of the chemical bond. Transverse position-dependent Stark shifts of ν(C═C) and ν(C≡N) in the individual nanocavity measured by surface-enhanced Raman scattering (SERS) experiment combined with the Stark tuning rate by density functional theory (DFT) simulation accurately revealed the inhomogeneous plasmonic field transverse distribution and quantified the transverse plasmonic field strength up to ∼1.9 × 10⁹ V/m, which matches the value predicted by finite element method (FEM) simulation. This work deepens the insight into plasmon-based technologies and will coordinate high-resolution techniques such as tip-enhanced Raman spectroscopy (TESR) to reveal the real plasmonic field distribution.

Force Decomposition Analysis: A Method to Decompose Intermolecular Forces into Physically Relevant Component Contributions

Computational quantum chemistry can be more than just numerical experiments when methods are specifically adapted to investigate chemical concepts. One important example is the development of energy decomposition analysis (EDA) to reveal the physical driving forces behind intermolecular interactions. In EDA, typically the interaction energy from a good-quality density functional theory (DFT) calculation is decomposed into multiple additive components that unveil permanent and induced electrostatics, Pauli repulsion, dispersion, and charge-transfer contributions to noncovalent interactions. Herein, we formulate, implement, and investigate decomposing the forces associated with intermolecular interactions into the same components. The resulting force decomposition analysis (FDA) is potentially useful as a complement to the EDA to understand chemistry, while also providing far more information than an EDA for data analysis purposes such as training physics-based force fields. We apply the FDA based on absolutely localized molecular orbitals (ALMOs) to analyze interactions of water with sodium and chloride ions as well as in the water dimer. We also analyze the forces responsible for geometric changes in carbon dioxide upon adsorption onto (and activation by) gold and silver anions. We also investigate how the force components of an EDA-based force field for water clusters, namely MB-UCB, compare to those from force decomposition analysis.

Electroabsorption and Stark Fluorescence Spectroscopies of Thioflavin T

Thioflavin T (ThT) is a typical fluorescent marker for detecting the formation of amyloid fibrils, because its fluorescence intensity increases by more than 2 orders of magnitude upon complexation with the fibrils. Strong electrostatic fields on protein surfaces are known to be a significant factor in chemical reactions and biological functions. Therefore, ThT bound to amyloid fibrils must experience strong electric fields. This study employed electroabsorption and Stark fluorescence spectroscopies to clarify the effects of external electric fields on the photophysics of ThT. The absorption spectrum shows two bands ascribed to locally excited (LE) and charge transfer (CT) states. Coupling between the LE and CT states is enhanced in the presence of an external electric field, resulting in fluorescence quenching. The electric field strength of the amyloid fibril surface was inferred from the fluorescence quenching efficiency of ThT.

TFRegNCI: Interpretable Noncovalent Interaction Correction Multimodal Based on Transformer Encoder Fusion

The interpretability is an important issue for end-to-end learning models. Motivated by computer vision algorithms, an interpretable noncovalent interaction (NCI) correction multimodal (TFRegNCI) is proposed for NCI prediction. TFRegNCI is based on RegNet feature extraction and a transformer encoder fusion strategy. RegNet is a network design paradigm that mainly focuses on local features. Meanwhile, the Vision Transformer is also leveraged for feature extraction, because it can capture global features better than RegNet while lowering the computational cost. Using a transformer encoder as the fusion strategy rather than multilayer perceptron can enhance model performance, due to its emphasis on important features with less parameters. Therefore, the proposed TFRegNCI achieved high accurate prediction (mean absolute error of ∼0.1 kcal/mol) comparing with the coupled cluster single double (triple) (CCSD(T)) benchmark. To further improve the model efficiency, TFRegNCI applies two-dimensional (2D) inputs transformed from three-dimensional (3D) electron density cubes, which saves time (30%), while the model accuracy remains. To improve model interpretability, a visualization module, Gradient-weighted Regression Activation Mapping (Grad-RAM) has been embedded. Grad-RAM is promoted from the classification algorithm, Gradient-weighted Class Activation Mapping, to perform feature visualization for the regression task. With Grad-RAM, the visual location map for features in deep learning models can be displayed. The feature map visualizations suggest that the 2D model has the similar performance as the 3D model, because of equally effective feature extractions from electron density. Moreover, the valid feature region on the location map by the 3D model is consistent with the NCIPLOT NCI isosurface. It is confirmed that the model does extract significant features related to the NCI interaction. The interpretable analyses are carried out through molecular orbital contribution on effective features. Thereby, the proposed model is likely to be a promising tool to reveal some essential information on NCIs, with regard to the level of electronic theory.

Rethinking Vibrational Stark Spectroscopy: Peak Shifts, Line Widths, and the Role of Non-Stark Solvent Coupling

A vibration's transition frequency is partly determined by the first-order Stark effect, which accounts for the electric field experienced by the mode. Using ultrafast infrared pump-probe and FT-IR spectroscopies, we characterized both the 0 → 1 and 1 → 2 vibrational transitions' field-dependent peak positions and line widths of the CN stretching mode of benzonitrile (BZN) and phenyl selenocyanate (PhSeCN) in ten solvents. We present a theoretical model that decomposes the observed line width into a field-dependent Stark contribution and a field-independent non-Stark solvent coupling contribution (NSC). The model demonstrates that the field-dependent peak position is independent of the line width, even when the NSC dominates the latter. Experiments show that when the Stark tuning rate is large compared to the NSC (PhSeCN), the line width has a field dependence, albeit with major NSC-induced excursions from linearity. When the Stark tuning rate is small relative to the NSC (BZN), the line width is field-independent. BZN's line widths are substantially larger for the 1 → 2 transition, indicating a 1 → 2 transition enhancement of the NSC. Additionally, we examine, theoretically and experimentally, the difference in the 0 → 1 and 1 → 2 transitions' Stark tuning rates. Second-order perturbation theory combined with density functional theory explain the difference and show that the 1 → 2 transition's Stark tuning rate is ∼10% larger. The Stark tuning rate of PhSeCN is larger than BZN's for both transitions, consistent with the theoretical calculations. This study provides new insights into vibrational line shape components and a more general understanding of the vibrational response to external electric fields.