Ultrafast terahertz Stark spectroscopy reveals the excited-state dipole moments of retinal in bacteriorhodopsin.

The photoinduced all-trans to 13-cis isomerization of the retinal Schiff base represents the ultrafast first step in the reaction cycle of bacteriorhodopsin (BR). Extensive experimental and theoretical work has addressed excited-state dynamics and isomerization via a conical intersection with the ground state. In conflicting molecular pictures, the excited state potential energy surface has been modeled as a pure S[Formula: see text] state that intersects with the ground state, or in a 3-state picture involving the S[Formula: see text] and S[Formula: see text] states. Here, the photoexcited system passes two crossing regions to return to the ground state. The electric dipole moment of the Schiff base in the S[Formula: see text] and S[Formula: see text] state differs strongly and, thus, its measurement allows for assessing the character of the excited-state potential. We apply the method of ultrafast terahertz (THz) Stark spectroscopy to measure electric dipole changes of wild-type BR and a BR D85T mutant upon electronic excitation. A fully reversible transient broadening and spectral shift of electronic absorption is induced by a picosecond THz field of several megavolts/cm and mapped by a 120-fs optical probe pulse. For both BR variants, we derive a moderate electric dipole change of 5 [Formula: see text] 1 Debye, which is markedly smaller than predicted for a neat S[Formula: see text]-character of the excited state. In contrast, S[Formula: see text]-admixture and temporal averaging of excited-state dynamics over the probe pulse duration gives a dipole change in line with experiment. Our results support a picture of electronic and nuclear dynamics governed by the interaction of S[Formula: see text] and S[Formula: see text] states in a 3-state model.

Transient Terahertz Stark Effect: A Dynamic Probe of Electric Interactions in Polar Liquids.

Electric forces acting on molecules in liquids at ambient temperature fluctuate at terahertz (THz) frequencies with a direct impact on their electronic and optical properties. We introduce the transient THz Stark effect to modify the electronic absorption spectra of dye molecules and, thus, elucidate and determine the underlying molecular interactions and dynamics. Picosecond electric fields of megavolts/cm induce a nonequilibrium response of the prototypical Betaine-30 in polar solution that is probed via transient absorption changes. The field-induced broadening of the absorption band follows the THz intensity in time, with a minor impact of solvent dynamics. The ground and excited state dipole energies in the THz field govern this response, allowing for a quantification of electric forces in a structurally frozen molecular environment.

Short-Range Cooperative Slow-down of Water Solvation Dynamics Around SO4 2--Mg2+ Ion Pairs.

The presence of ions affects the structure and dynamics of water on a multitude of length and time scales. In this context, pairs of Mg2+ and SO4 2- ions in water constitute a prototypical system for which conflicting pictures of hydration geometries and dynamics have been reported. Key issues are the molecular pair and solvation shell geometries, the spatial range of electric interactions, and their impact on solvation dynamics. Here, we introduce asymmetric SO4 2- stretching vibrations as new and most specific local probes of solvation dynamics that allow to access ion hydration dynamics at the dilute concentration (0.2 M) of a native electrolyte environment. Highly sensitive heterodyne 2D-IR spectroscopy in the fingerprint region of the SO4 2- ions around 1100 cm-1 reveals a specific slow-down of solvation dynamics for hydrated MgSO4 and for Na2SO4 in the presence of Mg2+ ions, which manifests as a retardation of spectral diffusion compared to aqueous Na2SO4 solutions in the absence of Mg2+ ions. Extensive molecular dynamics and density functional theory QM/MM simulations provide a microscopic view of the observed ultrafast dephasing and hydration dynamics. They suggest a molecular picture where the slow-down of hydration dynamics arises from the structural peculiarities of solvent-shared SO4 2--Mg2+ ion pairs.

Field-Induced Electron Generation in Water: Solvation Dynamics and Many-Body Interactions.

The solvated electron represents an elementary quantum system in a liquid environment. Electrons solvated in water have raised strong interest because of their prototypical properties, their role in radiation chemistry, and their relevance for charge separation and transport. Nonequilibrium dynamics of photogenerated electrons in water occur on ultrafast time scales and include charge transfer, localization, and energy dissipation processes. We present new insight into the role of fluctuating electric fields of the liquid for generating electrons in the presence of an external terahertz field and address polaronic many-body properties of solvated electrons. This Perspective combines a review of recent results from experiment and theory with a discussion of basic electric interactions of electrons in water.

Coherent polaron dynamics of electrons solvated in polar liquids.

An electron solvated in a polar liquid is an elementary quantum system with properties governed by electric interactions with a fluctuating molecular environment. In the prevailing single particle picture, the quantum ground and excited states are determined by a self-consistent potential, as defined by the particular local configuration of the solvation shell. This description neglects collective many-body excitations, which arise from the coupling of electronic degrees of freedom and nuclear motions of the environment. While recent experiments have demonstrated collective nonequilbrium electronic-nuclear motion, i.e. polaron excitations in liquid water, their relevance in the broader context of polar liquids has remained unexplored. Here, we study the nonequilibrium dielectric response of the, compared to water, less polar alcohols isopropanol, and ethylene glycol, that also display a different hydrogen bond pattern. We demonstrate that ultrafast relaxation of photogenerated electrons impulsively induces coherent charge oscillations, which persist for some 10 ps. They emit electric waves in a frequency range from 0.1 to 2 THz, depending on electron concentration. Oscillation frequencies and line shapes are reproduced by a unified polaron picture for alcohols and water, which is based on a Clausius-Mossotti local field approach for the THz dielectric function. The analysis suggests a longitudinal character of many-body polaron excitations and a weak coupling to transverse excitations, supported by the underdamped character of charge oscillations. Polaron dynamics are governed by the long-range Coulomb interaction between an excess electron and several thousands of polar solvent molecules, while local electron solvation geometries play a minor role.

The mutual interactions of RNA, counterions and water - quantifying the electrostatics at the phosphate-water interface.

The structure and dynamics of polyanionic biomolecules, like RNA, are decisively determined by their electric interactions with the water molecules and the counterions in the environment. The solvation dynamics of the biomolecules involves a subtle balance of non-covalent and many-body interactions with structural fluctuations due to thermal motion occurring in a femto- to subnanosecond time range. This complex fluctuating many particle scenario is crucial in defining the properties of biological interfaces with far reaching significance for the folding of RNA structures and for facilitating RNA-protein interactions. Given the inherent complexity, suited model systems, carefully calibrated and benchmarked by experiments, are required to quantify the relevant interactions of RNA with the aqueous environment. In this feature article we summarize our recent progress in the understanding of the electrostatics at the biological interface of double stranded RNA (dsRNA) and transfer RNA (tRNA). Dimethyl phosphate (DMP) is introduced as a viable and rigorously accessible model system allowing the interaction strength with water molecules and counterions, their relevant fluctuation timescales and the spatial reach of interactions to be established. We find strong (up to ≈90 MV cm-1) interfacial electric fields with fluctuations extending up to ≈20 THz and demonstrate how the asymmetric stretching vibration νAS(PO2)- of the polarizable phosphate group can serve as the most sensitive probe for interfacial interactions, establishing a rigorous link between simulations and experiment. The approach allows for the direct interfacial observation of interactions of biologically relevant Mg2+ counterions with phosphate groups in contact pair geometries via the rise of a new absorption band imposed by exchange repulsion interactions at short interatomic distances. The systematic extension to RNA provides microscopic insights into the changes of the hydration structure that accompany the temperature induced melting of the dsRNA double helix and quantify the ionic interactions in the folded tRNA. The results show that pairs of negatively charged phosphate groups and Mg2+ ions represent a key structural feature of RNA embedded in water. They highlight the importance of binding motifs made of contact pairs in the electrostatic stabilization of RNA structures that have a strong impact on the surface potential and enable the fine tuning of the local electrostatic properties which are expected to be relevant for mediating the interactions between biomolecules.

Contact pairs of RNA with magnesium ions-electrostatics beyond the Poisson-Boltzmann equation.

The electrostatic interaction of RNA with its aqueous environment is most relevant for defining macromolecular structure and biological function. The attractive interaction of phosphate groups in the RNA backbone with ions in the water environment leads to the accumulation of positively charged ions in the first few hydration layers around RNA. Electrostatics of this ion atmosphere and the resulting ion concentration profiles have been described by solutions of the nonlinear Poisson-Boltzmann equation and atomistic molecular dynamics (MD) simulations. Much less is known on contact pairs of RNA phosphate groups with ions at the RNA surface, regarding their abundance, molecular geometry, and role in defining RNA structure. Here, we present a combined theoretical and experimental study of interactions of a short RNA duplex with magnesium (Mg2+) ions. MD simulations covering a microsecond time range give detailed hydration geometries as well as electrostatics and spatial arrangements of phosphate-Mg2+ pairs, including both pairs in direct contact and separated by a single water layer. The theoretical predictions are benchmarked by linear infrared absorption and nonlinear two-dimensional infrared spectra of the asymmetric phosphate stretch vibration which probes both local interaction geometries and electric fields. Contact pairs of phosphate groups and Mg2+ ions are identified via their impact on the vibrational frequency position and line shape. A quantitative analysis of infrared spectra for a range of Mg2+-excess concentrations and comparison with fluorescence titration measurements shows that on average 20-30% of the Mg2+ ions interacting with the RNA duplex form contact pairs. The experimental and MD results are in good agreement. In contrast, calculations based on the nonlinear Poisson-Boltzmann equation fail in describing the ion arrangement, molecular electrostatic potential, and local electric field strengths correctly. Our results underline the importance of local electric field mapping and molecular-level simulations to correctly account for the electrostatics at the RNA-water interface.

Electron Transfer-Induced Active Site Structural Relaxation in 64-Photolyase of Drosophila melanogaster.

While catalytic electron flow and photoreactivation of CPD-photolyases are increasingly understood, the microscopic details of the 64-photolyase repair mechanism are perpetually debated. Here, we investigate in long-time (µs) molecular dynamics simulations combined with extensive quantum mechanical/molecular mechanical (QM/MM) simulations the primary electron transfer (ET) reactions in 64-photolyase of Drosophila melanogaster (D. melanogaster). The characterization of the relative energetics of locally excited and charge separated states in the (6-4) photoproduct enzyme repair complex reveals a charge-separated state involving the adenine moiety of the FADH- cofactor that facilitates reduction of the photoproduct. Microscopic details of the collective reaction coordinate of ET reactions are identified that involve the reorganization of the hydrogen bond network and structural relaxation of the active site. The simulations reveal complex active site relaxation dynamics involving distinguished amino acids (Lys246, His365, and His369), conformational reorganization of the hydroxyl group of the (6-4) photoproduct, and a strengthening of hydrogen bonds with immobilized water molecules. In particular, rotation of the Lys246 side chain is found to impose a double-well character along the reaction coordinate of the ET reaction. Our findings suggest that the primary ET reactions in the (6-4) photoproduct enzyme repair complex of D. melanogaster are governed by a complex multi-minima active site relaxation dynamics and potentially precede the equilibration of the protein. ET pathways mediated by the adenine moiety and the 5' side of the photoproduct are proposed to be relevant for triggering the catalytic (6-4) photoproduct reactivation.

Phosphate Vibrations Probe Electric Fields in Hydrated Biomolecules: Spectroscopy, Dynamics, and Interactions.

Electric interactions have a strong impact on the structure and dynamics of biomolecules in their native water environment. Given the variety of water arrangements in hydration shells and the femto- to subnanosecond time range of structural fluctuations, there is a strong quest for sensitive noninvasive probes of local electric fields. The stretching vibrations of phosphate groups, in particular the asymmetric (PO2)- stretching vibration νAS(PO2)-, allow for a quantitative mapping of dynamic electric fields in aqueous environments via a field-induced redshift of their transition frequencies and concomitant changes of vibrational line shapes. We present a systematic study of νAS(PO2)- excitations in molecular systems of increasing complexity, including dimethyl phosphate (DMP), short DNA and RNA duplex structures, and transfer RNA (tRNA) in water. A combination of linear infrared absorption, two-dimensional infrared (2D-IR) spectroscopy, and molecular dynamics (MD) simulations gives quantitative insight in electric-field tuning rates of vibrational frequencies, electric field and fluctuation amplitudes, and molecular interaction geometries. Beyond neat water environments, the formation of contact ion pairs of phosphate groups with Mg2+ ions is demonstrated via frequency upshifts of the νAS(PO2)- vibration, resulting in a distinct vibrational band. The frequency positions of contact geometries are determined by an interplay of attractive electric and repulsive exchange interactions.

Infrared and NMR Spectroscopic Fingerprints of the Asymmetric H7 + O3 Complex in Solution.

The front cover artwork is provided by the groups of Prof. Ehud Pines (BGU, Israel) and Dr. Benjamin Fingerhut (MBI, Berlin). The image shows a scientist integrating experiments with theory for resolving the structural diffusion of the aqueous proton in acetonitrile providing a novel view on the Grotthuss mechanism. Read the full text of the Article at 10.1002/cphc.202001046.

Terahertz Polaron Oscillations of Electrons Solvated in Liquid Water.

The terahertz (THz) response of solvated electrons in liquid water is studied in nonlinear ultrafast pump-probe experiments. Free electrons with concentrations from c_{e}=4 to 140×10^{-6} moles/liter are generated by high-field THz or near-infrared multiphoton excitation. The time-resolved change of the dielectric function as mapped by broadband THz pulses exhibits pronounced oscillations persisting up to 30 ps. Their frequency increases with electron concentration from 0.2 to 1.5 THz. The oscillatory response is assigned to impulsively excited coherent polarons involving coupled electron and water shell motions with a frequency set by the local electric field.

Infrared and NMR Spectroscopic Fingerprints of the Asymmetric H7+ O3 Complex in Solution.

Infrared (IR) absorption in the 1000-3700 cm-1 range and 1 H NMR spectroscopy reveal the existence of an asymmetric protonated water trimer, H7+ O3, in acetonitrile. The core H7+ O3 motif persists in larger protonated water clusters in acetonitrile up to at least 8 water molecules. Quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations reveal irreversible proton transport promoted by propagating the asymmetric H7+ O3 structure in solution. The QM/MM calculations allow for the successful simulation of the measured IR absorption spectra of H7+ O3 in the OH stretch region, which reaffirms the assignment of the H7+ O3 spectra to a hybrid-complex structure: a protonated water dimer strongly hydrogen-bonded to a third water molecule with the proton exchanging between the two possible shared-proton Zundel-like centers. The H7+ O3 structure lends itself to promoting irreversible proton transport in presence of even one additional water molecule. We demonstrate how continuously evolving H7+ O3 structures may support proton transport within larger water solvates.

Magnesium Contact Ions Stabilize the Tertiary Structure of Transfer RNA: Electrostatics Mapped by Two-Dimensional Infrared Spectra and Theoretical Simulations.

Ions interacting with hydrated RNA play a central role in defining its secondary and tertiary structure. While spatial arrangements of ions, water molecules, and phosphate groups have been inferred from X-ray studies, the role of electrostatic and other noncovalent interactions in stabilizing compact folded RNA structures is not fully understood at the molecular level. Here, we demonstrate that contact ion pairs of magnesium (Mg2+) and phosphate groups embedded in local water shells stabilize the tertiary equilibrium structure of transfer RNA (tRNA). Employing dialyzed tRNAPhe from yeast and tRNA from Escherichia coli, we follow the population of Mg2+ sites close to phosphate groups of the ribose-phosphodiester backbone step by step, combining linear and nonlinear infrared spectroscopy of phosphate vibrations with molecular dynamics simulations and ab initio vibrational frequency calculations. The formation of up to six Mg2+/phosphate contact pairs per tRNA and local field-induced reorientations of water molecules balance the phosphate-phosphate repulsion in nonhelical parts of tRNA, thus stabilizing the folded structure electrostatically. Such geometries display limited sub-picosecond fluctuations in the arrangement of water molecules and ion residence times longer than 1 µs. At higher Mg2+ excess, the number of contact ion pairs per tRNA saturates around 6 and weakly interacting ions prevail. Our results suggest a predominance of contact ion pairs over long-range coupling of the ion atmosphere and the biomolecule in defining and stabilizing the tertiary structure of tRNA.

On the role of non-diagonal system-environment interactions in bridge-mediated electron transfer.

Bridge-mediated electron transfer (ET) between a donor and an acceptor is prototypical for the description of numerous most important ET scenarios. While multi-step ET and the interplay of sequential and direct superexchange transfer pathways in the donor-bridge-acceptor (D-B-A) model are increasingly understood, the influence of off-diagonal system-bath interactions on the transfer dynamics is less explored. Off-diagonal interactions account for the dependence of the ET coupling elements on nuclear coordinates (non-Condon effects) and are typically neglected. Here, we numerically investigate with quasi-adiabatic propagator path integral simulations the impact of off-diagonal system-environment interactions on the transfer dynamics for a wide range of scenarios in the D-B-A model. We demonstrate that off-diagonal system-environment interactions can have profound impact on the bridge-mediated ET dynamics. In the considered scenarios, the dynamics itself does not allow for a rigorous assignment of the underlying transfer mechanism. Furthermore, we demonstrate how off-diagonal system-environment interaction mediates anomalous localization by preventing long-time depopulation of the bridge B and how coherent transfer dynamics between donor D and acceptor A can be facilitated. The arising non-exponential short-time dynamics and coherent oscillations are interpreted within an equivalent Hamiltonian representation of a primary reaction coordinate model that reveals how the complex vibronic interplay of vibrational and electronic degrees of freedom underlying the non-Condon effects can impose donor-to-acceptor coherence transfer on short timescales.

Regulatory Impact of the C-Terminal Tail on Charge Transfer Pathways in Drosophila Cryptochrome.

Interconnected transcriptional and translational feedback loops are at the core of the molecular mechanism of the circadian clock. Such feedback loops are synchronized to external light entrainment by the blue light photoreceptor cryptochrome (CRY) that undergoes conformational changes upon light absorption by an unknown photoexcitation mechanism. Light-induced charge transfer (CT) reactions in Drosophila CRY (dCRY) are investigated by state-of-the-art simulations that reveal a complex, multi-redox site nature of CT dynamics on the microscopic level. The simulations consider redox-active chromophores of the tryptophan triad (Trp triad) and further account for pathways mediated by W314 and W422 residues proximate to the C-terminal tail (CTT), thus avoiding a pre-bias to specific W-mediated CT pathways. The conducted dissipative quantum dynamics simulations employ microscopically derived model Hamiltonians and display complex and ultrafast CT dynamics on the picosecond timescale, subtly balanced by the electrostatic environment of dCRY. In silicio point mutations provide a microscopic basis for rationalizing particular CT directionality and demonstrate the degree of electrostatic control realized by a discrete set of charged amino acid residues. The predicted participation of CT states in proximity to the CTT relates the directionality of CT reactions to the spatial vicinity of a linear interaction motif. The results stress the importance of CTT directional charge transfer in addition to charge transfer via the Trp triad and call for the use of full-length CRY models including the interactions of photolyase homology region (PHR) and CTT domains.

Field-Induced Tunneling Ionization and Terahertz-Driven Electron Dynamics in Liquid Water.

Liquid water at ambient temperature displays ultrafast molecular motions and concomitant fluctuations of very strong electric fields originating from the dipolar H2O molecules. We show that such random intermolecular fields induce the tunnel ionization of water molecules, which becomes irreversible if an external terahertz (THz) pulse imposes an additional directed electric field on the liquid. Time-resolved nonlinear THz spectroscopy maps charge separation, transport, and localization of the released electrons on a few-picosecond time scale. The highly polarizable localized electrons modify the THz absorption spectrum and refractive index of water, a manifestation of a highly nonlinear response. Our results demonstrate how the interplay of local electric field fluctuations and external electric fields allows for steering charge dynamics and dielectric properties in aqueous systems.

Change of Hydration Patterns upon RNA Melting Probed by Excitations of Phosphate Backbone Vibrations.

The water hydration shell has a decisive impact on the structural and functional properties of RNA. Changes of the RNA structure upon melting and in biochemical processes are accompanied by a change of hydration patterns, a process that is barely characterized. To discern hydration geometries around the backbone phosphate groups of an RNA double helix at the molecular level, we combine two-dimensional infrared spectroscopy of phosphate vibrations with theoretical simulations. There are three distinct coexisting hydration motifs of the RNA A-helix: an ordered chain-like arrangement of water molecules with links between neighboring phosphate groups, separate local hydration shells of up to six water molecules, and hydrated phosphate/counterion contact pairs. RNA disordering upon melting is connected with a transition from predominantly ordered water structures to local hydration shells around phosphate units. Structural fluctuations are dominated by librational water motions occurring on a 300 fs time scale, without exchange between hydration motifs.

Contact Ion Pairs of Phosphate Groups in Water: Two-Dimensional Infrared Spectroscopy of Dimethyl Phosphate and ab Initio Simulations.

The interaction of phosphate groups with ions in an aqueous environment has a strong impact on the structure and folding processes of DNA and RNA. The dynamic variety of ionic arrangements, including both contact pairs and water separated ions, and the molecular coupling mechanisms are far from being understood. In a combined experimental and theoretical approach, we address the properties of contact ion pairs of the prototypical system dimethyl phosphate with Na+, Ca2+, and Mg2+ ions in water. Linear and femtosecond two-dimensional infrared (2D-IR) spectroscopy of the asymmetric (PO2)- stretching vibration separates and characterizes the different species via their blue-shifted vibrational signatures and 2D-IR line shapes. Phosphate-magnesium contact pairs stand out as the most compact geometry while the contact pairs with Ca2+ and Na+ display a wider structural variation. Microscopic density functional theory simulations rationalize the observed frequency shifts and reveal distinct differences between the contact geometries.

Coupled excitation energy and charge transfer dynamics in reaction centre inspired model systems.

Functional operation conditions of reaction centre core complexes require the tight coupling of exciton states to concomitant charge separation. Rigorous theoretical treatment of such integrated excitation energy transfer (EET) and charge transfer (CT) dynamics is particularly challenging due to (i) appreciable system sizes, (ii) inter-site and system-bath couplings of similar magnitude that render the Born-Markov approximation invalid, (iii) substantial reorganization energies of CT states, and (iv) the presence of complex structured spectral densities due to vibrational modes of the surroundings. We present numerical simulations on bacterial reaction centre (bRC) inspired model systems that utilize the recently developed MACGIC-iQUAPI method [Richter et al., J. Chem. Phys., 2017, 146, 214101]. The simulations demonstrate that the method provides a rigorous framework for the investigation of such integrated EET-CT dynamics. First, the applicability of the MACGIC-iQUAPI method is explored for a transition from monotonically decaying to oscillatory system-bath influence coefficients, a behavior inherently imposed by structured bath spectral densities. Tightly coupled EET and CT dynamics is further addressed for an excitonic subsystem that resembles strong coupling of special pair states and serves as donor towards a generic bridge-acceptor system. By solving the dissipative quantum dynamics of such bRC inspired model systems, the quenching of excitonic coherence on the hundreds of femtoseconds timescale is explored via a variation of the bridge state energetics, resembling a continuous transition from sequential to superexchange mediated CT regimes. Further, the simulations explore the influence of resonant vibrational modes on the quenching of excitonic coherence via CT. The results reveal a moderate influence of vibrational mode on charge separation dynamics in regimes of biologically relevant EET and CT dynamics.