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

Noncanonical Amino Acids in Biocatalysis

In recent years, powerful genetic code reprogramming methods have emerged that allow new functional components to be embedded into proteins as noncanonical amino acid (ncAA) side chains. In this review, we will illustrate how the availability of an expanded set of amino acid building blocks has opened a wealth of new opportunities in enzymology and biocatalysis research. Genetic code reprogramming has provided new insights into enzyme mechanisms by allowing introduction of new spectroscopic probes and the targeted replacement of individual atoms or functional groups. NcAAs have also been used to develop engineered biocatalysts with improved activity, selectivity, and stability, as well as enzymes with artificial regulatory elements that are responsive to external stimuli. Perhaps most ambitiously, the combination of genetic code reprogramming and laboratory evolution has given rise to new classes of enzymes that use ncAAs as key catalytic elements. With the framework for developing ncAA-containing biocatalysts now firmly established, we are optimistic that genetic code reprogramming will become a progressively more powerful tool in the armory of enzyme designers and engineers in the coming years.

Local electric field in nanocavities dictates the vibrational relaxation dynamics of interfacial molecules

Plasmonic nanocavities enable the generation of strong light-matter coupling and exhibit great potential in plasmon-mediated chemical reactions (PMCRs). Although an electric field generated by nanocavities (E n) has recently been reported, its effect on the vibrational energy relaxation (VER) of the molecules in the nanocavities has not been explored. In this study, we reveal the impact of an electric field sensed by molecules (para-substituted thiophenol derivatives) in a nanocavity (E f) on VER processes by employing advanced time-resolved femtosecond sum frequency generation vibrational spectroscopy (SFG-VS) supplemented by electrochemical measurements. The magnitude of E n is almost identical (1.0 ± 0.2 V nm-1) beyond the experimental deviation while E f varies from 0.3 V nm-1 to 1.7 V nm-1 depending on the substituent. An exponential correlation between E f and the complete recovery time of the ground vibrational C[double bond, length as m-dash]C state (T 2) of the phenyl ring is observed. Substances with a smaller T 2 are strongly correlated with the reported macroscopic chemical reactivity. This finding may aid in enriching the current understanding of PMCRs and highlights the possibility of regulating vibrational energy flow into desired reaction coordinates by using a local electric field.

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.

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.

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.

Development of a universal method for vibrational analysis of the terminal alkyne C≡C stretch

The terminal alkyne C≡C stretch has a large Raman scattering cross section in the "silent" region for biomolecules. This has led to many Raman tag and probe studies using this moiety to study biomolecular systems. A computational investigation of these systems is vital to aid in the interpretation of these results. In this work, we develop a method for computing terminal alkyne vibrational frequencies and isotropic transition polarizabilities that can easily and accurately be applied to any terminal alkyne molecule. We apply the discrete variable representation method to a localized version of the C≡C stretch normal mode. The errors of (1) vibrational localization to the terminal alkyne moiety, (2) anharmonic normal mode isolation, and (3) discretization of the Born-Oppenheimer potential energy surface are quantified and found to be generally small and cancel each other. This results in a method with low error compared to other anharmonic vibrational methods like second-order vibrational perturbation theory and to experiments. Several density functionals are tested using the method, and TPSS-D3, an inexpensive nonempirical density functional with dispersion corrections, is found to perform surprisingly well. Diffuse basis functions are found to be important for the accuracy of computed frequencies. Finally, the computation of vibrational properties like isotropic transition polarizabilities and the universality of the localized normal mode for terminal alkynes are demonstrated.

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.

The hydrogen-bonding dynamics of water to a nitrile-functionalized electrode is modulated by voltage according to ultrafast 2D IR spectroscopy

We report the hydrogen-bonding dynamics of water to a nitrile-functionalized and plasmonic electrode surface as a function of applied voltage. The surface-enhanced two-dimensional infrared spectra exhibit hydrogen-bonded and non-hydrogen-bonded nitrile features in similar proportions, plus cross peaks between the two. Isotopic dilution experiments show that the cross peaks arise predominantly from chemical exchange between hydrogen-bonded and non-hydrogen-bonded nitriles. The chemical exchange rate depends upon voltage, with the hydrogen bond of the water to the nitriles breaking 2 to 3 times slower (>63 vs. 25 ps) under a positive as compared to a negative potential. Spectral diffusion created by hydrogen-bond fluctuations occurs on a ~1 ps timescale and is moderately potential-dependent. Timescales from molecular dynamics simulations agree qualitatively with the experiment and show that a negative voltage causes a small net displacement of water away from the surface. These results show that the voltage applied to an electrode can alter the timescales of solvent motion at its interface, which has implications for electrochemically driven reactions.

Molecular overlay-guided design of new CDK2 inhibitor thiazepinopurines: Synthesis, anticancer, and mechanistic investigations

Adopting the molecular overlay approach, three novel sets of thiazepinopurines with expected cytotoxicity and CDK2 inhibition potential were designed and synthesized. This was accomplished through the heteroannelation of purines, for the first time, with thiazepine. The obtained thiazepinopurines derivatives were assessed for their cytotoxicity toward tumor cells of three different types, HepG2, MCF-7, and PC-3 as well as one normal cell (WI38). Among the studied compounds, 3b and 3c exhibited significant antiproliferative activity against tumor cells presenting IC50 range of 5.52-17.09 µM in comparison with Roscovitine (9.32-13.82 µM). Additionally, both compounds displayed superior selectivity indices (SI = 3.00-7.15) toward tested cancer cells. The 4-chlorophenyl analog 3b has shown the best selectivity index, and hence it has been subjected to additional investigations to determine its proper mechanistic effect. Accordingly, the CDK2 inhibition potential, apoptosis induction, and cell cycle analysis of MCF-7 were evaluated. Results revealed that this analog displayed a potent CDK2 inhibition potential with an IC50 value of 0.219 µM. Findings also showed that 3b was thought to arrest MCF-7 cell cycle at S phase together with apoptosis induction by the increased expression of Bax, Caspase-8, and -9 markers with a concomitant decrease in Bcl-2 expression. Besides, the probable interaction of 3b with CDK2 binding pocket was investigated by molecular docking.

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.

Solid state solvation: a fresh view

The design of efficient organic electronic devices, including OLEDs, OPVs, luminescent solar concentrators, etc., relies on the optimization of relevant materials, often constituted by an active (functional) dye embedded in a matrix. Understanding solid state solvation (SSS), i.e. how the properties of the active dye are affected by the matrix, is therefore an issue of fundamental and technological relevance. Here an extensive experimental and theoretical investigation is presented shedding light on this, somewhat controversial, topic. The spectral properties of the dye at equilibrium, i.e. absorption and Raman spectra, are not affected by the matrix dynamics. Reliable estimates of the matrix polarity are then obtained from an analysis of the micro-Raman spectra of polar dyes. Specifically, to establish a reliable polarity scale, the spectra of DCM or NR dispersed in amorphous matrices are compared with the spectra of the same dyes in liquid solvents with known polarity. On the other hand, steady-state emission spectra obtained in solid matrices depend in a highly non-trivial way on the matrix polarity and its dynamics. An extensive experimental and theoretical analysis of the time-resolved emission spectra of NR in a very large time window (15 fs-15 ns) allows us to validate this dye as a good probe of the dielectric dynamics of the surrounding medium. We provide a first assessment of the relaxation dynamics of two matrices (mCBPCN and DPEPO) of interest for OLED application, unambiguously demonstrating that the matrix readjusts for at least 15 ns after the dye photoexcitation.

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.

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.

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.

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.

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

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

Spectroscopic Characterization of the Divalent Metal Docking Motif to Isolated Cyanobenzoate: Direct Observation of Tridentate Binding to ortho-Cyanobenzoate and Implications for the CN Response

Cryogenic ion vibrational spectra of D₂-tagged cyanobenzoate (CBA) derivatives are obtained and analyzed to characterize the intrinsic spectroscopic responses of the -CO₂^- headgroup to its location on the ring in both the isolated anions and the cationic complexes with divalent metal ions, M²⁺ (M = Mg, Ca, Sr). The benzonitrile functionality establishes the different ring isomers (para, meta, ortho) according to the location of the carboxylate and provides an additional reporter on the molecular response to the proximal charge center. The aromatic carboxylates display shifts slightly smaller than those observed for a related aliphatic system upon metal ion complexation. Although the CBA anions display very similar band patterns for all three ring positions, upon complexation with metal ions, the ortho isomer yields dramatically different spectral responses in both the -CO₂^- moiety and the CN group. This behavior is traced to the emergence of a tridentate binding motif unique to the ortho isomer in which the metal ions bind to both the oxygen atoms of the carboxylate group and the N atom of the cyano group. In that configuration, the -CO₂^- moiety is oriented perpendicular to the phenyl ring, and the CN stretching fundamental is both strong and red-shifted relative to its behavior in the isolated neutral. The behaviors of the metal-bound ortho complexes occur in contrast to the usual blue shifts associated with "Lewis" type binding of metal ions end-on to -CN. The origins of these spectroscopic features are analyzed with the aid of electronic structure calculations, which also explore differences expected for complexation of monovalent cations to the ortho carboxylate. The resulting insights have implications for understanding the balance between electrostatic and steric interactions at metal binding sites in chemical and biological systems.

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

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

Genetically encoded non-canonical amino acids reveal asynchronous dark reversion of chromophore, backbone and side-chains in EL222

Photoreceptors containing the light-oxygen-voltage (LOV) domain elicit biological responses upon excitation of their flavin mononucleotide (FMN) chromophore by blue light. The mechanism and kinetics of dark-state recovery are not well understood. Here we incorporated the non-canonical amino acid p-cyanophenylalanine (CNF) by genetic code expansion technology at forty-five positions of the bacterial transcription factor EL222. Screening of light-induced changes in infrared (IR) absorption frequency, electric field and hydration of the nitrile groups identified residues CNF31 and CNF35 as reporters of monomer/oligomer and caged/decaged equilibria, respectively. Time-resolved multi-probe UV/Visible and IR spectroscopy experiments of the lit-to-dark transition revealed four dynamical events. Predominantly, rearrangements around the A'α helix interface (CNF31 and CNF35) precede FMN-cysteinyl adduct scission, folding of α-helices (amide bands), and relaxation of residue CNF151. This study illustrates the importance of characterizing all parts of a protein and suggests a key role for the N-terminal A'α extension of the LOV domain in controlling EL222 photocycle length. This article is protected by copyright. All rights reserved.

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

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

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.

Potential Distribution across Model Membranes

Membrane models assembled on electrodes are widely used tools to study potential-dependent molecular processes at or in membranes. However, the relationship between the electrode potential and the potential across the membrane is not known. Here we studied lipid bilayers immobilized on mixed self-assembled monolayers (SAM) on Au electrodes. The mixed SAM was composed of thiol derivatives of different chain lengths such that between the islands of the short one, mercaptobenzonitrile (MBN), and the tethered lipid bilayer an aqueous compartment was formed. The nitrile function of MBN, which served as a reporter group for the vibrational Stark effect (VSE), was probed by surface-enhanced infrared absorption spectroscopy to determine the local electric field as a function of the electrode potential for pure MBN, mixed SAM, and the bilayer system. In parallel, we calculated electric fields at the VSE probe by molecular dynamics (MD) simulations for different charge densities on the metal, thereby mimicking electrode potential changes. The agreement with the experiments was very good for the calculations of the pure MBN SAM and only slightly worse for the mixed SAM. The comparison with the experiments also guided the design of the bilayer system in the MD setups, which were selected to calculate the electrode potential dependence of the transmembrane potential, a quantity that is not directly accessible by the experiments. The results agree very well with estimates in previous studies and thus demonstrate that the present combined experimental-theoretical approach is a promising tool for describing potential-dependent processes at biomimetic interfaces.

Electric fields control water-gated proton transfer in cytochrome c oxidase

Aerobic life is powered by membrane-bound enzymes that catalyze the transfer of electrons to oxygen and protons across a biological membrane. Cytochrome c oxidase (CcO) functions as a terminal electron acceptor in mitochondrial and bacterial respiratory chains, driving cellular respiration and transducing the free energy from O₂ reduction into proton pumping. Here we show that CcO creates orientated electric fields around a nonpolar cavity next to the active site, establishing a molecular switch that directs the protons along distinct pathways. By combining large-scale quantum chemical density functional theory (DFT) calculations with hybrid quantum mechanics/molecular mechanics (QM/MM) simulations and atomistic molecular dynamics (MD) explorations, we find that reduction of the electron donor, heme a, leads to dissociation of an arginine (Arg438)-heme a₃ D-propionate ion-pair. This ion-pair dissociation creates a strong electric field of up to 1 V Å^-1 along a water-mediated proton array leading to a transient proton loading site (PLS) near the active site. Protonation of the PLS triggers the reduction of the active site, which in turn aligns the electric field vectors along a second, "chemical," proton pathway. We find a linear energy relationship of the proton transfer barrier with the electric field strength that explains the effectivity of the gating process. Our mechanism shows distinct similarities to principles also found in other energy-converting enzymes, suggesting that orientated electric fields generally control enzyme catalysis.