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.

Managing the Redox Potential of PCET in Grotthuss-Type Proton Wires

Expanding proton-coupled electron transfer to multiproton translocations (MPCET) provides a bioinspired mechanism to transport protons away from the redox site. This expansion has been accomplished by separating the initial phenolic proton donor from the pyridine-based terminal proton acceptor by a Grotthuss-type proton wire made up of concatenated benzimidazoles that form a hydrogen-bonded network. However, it was found that the midpoint potential of the phenol oxidation that launched the Grotthuss-type proton translocations is a function of the number of benzimidazoles in the hydrogen-bonded network; it becomes less positive (i.e., a weaker oxidant) as the number of bridging benzimidazoles increases. Herein, we report a strategy to maintain the high redox potential necessary for oxidative processes relevant to artificial photosynthesis, e.g., water oxidation and long-range MPCET processes for managing protons. The integrated structural and functional roles of the benzimidazole-based bridge provide sites for substitution of the benzimidazoles with electron-withdrawing groups (e.g., trifluoromethyl groups). Such substitution increases the midpoint potential of the phenoxyl radical/phenol couple so that proton translocations over ∼11 Å become thermodynamically comparable to that of an unsubstituted system where one proton is transferred over ∼2.5 Å. The extended, substituted system maintains the hydrogen-bonded network; infrared spectroelectrochemistry confirms reversible proton translocations from the phenol to the pyridyl terminal proton acceptor upon oxidation and reduction. Theory supports the change in driving force with added electron-withdrawing groups and provides insight into the role of electron density and electrostatic potential in MPCET processes associated with these Grotthuss-type proton translocations.

Artificial Neural Networks as Propagators in Quantum Dynamics

The utilization of artificial neural networks (ANNs) provides strategies for accelerating molecular simulations. Herein, ANNs are implemented as propagators of the time-dependent Schrödinger equation to simulate the quantum dynamics of systems with time-dependent potentials. These ANN propagators are trained to map nonstationary wavepackets from a given time to a future time within the discrete variable representation. Each propagator is trained for a specified time step, and iterative application of the propagator enables the propagation of wavepackets over long time scales. Such ANN propagators are developed and applied to one- and two-dimensional proton transfer systems, which exhibit nuclear quantum effects such as hydrogen tunneling. These ANN propagators are trained for either a specific time-independent potential or general potentials that can be time-dependent. Hierarchical, multiple time step algorithms enable parallelization, and the extension to higher dimensions is straightforward. This strategy is applicable to quantum dynamical simulations of diverse chemical and biological processes.

Multi PCET in symmetrically substituted benzimidazoles

Proton-coupled electron transfer (PCET) reactions depend on the hydrogen-bond connectivity between sites of proton donors and acceptors. The 2-(2′-hydroxyphenyl) benzimidazole (BIP) based systems, which mimic the natural Tyr_Z-His190 pair of Photosystem II, have been useful for understanding the associated PCET process triggered by one-electron oxidation of the phenol. Substitution of the benzimidazole by an appropriate terminal proton acceptor (TPA) group allows for two-proton translocations. However, the prototropic properties of substituted benzimidazole rings and rotation around the bond linking the phenol and the benzimidazole can lead to isomers that interrupt the intramolecular hydrogen-bonded network and thereby prevent a second proton translocation. Herein, a strategic symmetrization of a benzimidazole based system with two identical TPAs yields an uninterrupted network of intramolecular hydrogen bonds regardless of the isomeric form. NMR data confirms the presence of a single isomeric form in the disubstituted system but not in the monosubstituted system in certain solvents. Infrared spectroelectrochemistry demonstrates a two-proton transfer process associated with the oxidation of the phenol occurring at a lower redox potential in the disubstituted system relative to its monosubstituted analogue. Computational studies support these findings and show that the disubstituted system stabilizes the oxidized two-proton transfer product through the formation of a bifurcated hydrogen bond. Considering the prototropic properties of the benzimidazole heterocycle in the context of multiple PCET will improve the next generation of novel, bioinspired constructs built by concatenated units of benzimidazoles, thus allowing proton translocations at nanoscale length.

Proton-coupled electron transfer (PCET) reactions depend on the hydrogen-bond connectivity between sites of proton donors and acceptors.

Artificial Neural Networks as Mappings between Proton Potentials, Wave Functions, Densities, and Energy Levels

Artificial neural networks (ANNs) have become important in quantum chemistry. Herein, applications to nuclear quantum effects, such as zero-point energy, vibrationally excited states, and hydrogen tunneling, are explored. ANNs are used to solve the time-independent Schrödinger equation for single- and double-well potentials representing hydrogen-bonded molecular systems capable of proton transfer. ANN mappings are trained to predict the lowest five proton vibrational energies, wave functions, and densities from the proton potentials and to predict the excited state proton vibrational energies and densities from the proton ground state density. For the inverse problem, ANN mappings are trained to predict the proton potential from the proton vibrational energy levels or the proton ground state density. This latter mapping is theoretically justified by the first Hohenberg-Kohn theorem establishing a one-to-one correspondence between the external potential and the ground state density. ANNs for two- and three-dimensional systems are also presented to illustrate the straightforward extension to higher dimensions.

Inhomogeneity of Interfacial Electric Fields at Vibrational Probes on Electrode Surfaces

Electric fields control chemical reactivity in a wide range of systems, including enzymes and electrochemical interfaces. Characterizing the electric fields at electrode-solution interfaces is critical for understanding heterogeneous catalysis and associated energy conversion processes. To address this challenge, recent experiments have probed the response of the nitrile stretching frequency of 4-mercaptobenzonitrile (4-MBN) attached to a gold electrode to changes in the solvent and applied electrode potential. Herein, this system is modeled with periodic density functional theory using a multilayer dielectric continuum treatment of the solvent and at constant applied potentials. The impact of the solvent dielectric constant and the applied electrode potential on the nitrile stretching frequency computed with a grid-based method is in qualitative agreement with the experimental data. In addition, the interfacial electrostatic potentials and electric fields as a function of applied potential were calculated directly with density functional theory. Substantial spatial inhomogeneity of the interfacial electric fields was observed, including oscillations in the region of the molecular probe attached to the electrode. These simulations highlight the microscopic inhomogeneity of the electric fields and the role of molecular polarizability at electrode-solution interfaces, thereby demonstrating the limitations of mean-field models and providing insights relevant to the interpretation of vibrational Stark effect experiments.