Search

Protein dynamics form a critical bridge between protein structure and function, yet the impact of evolution on ultrafast processes inside proteins remains enigmatic. This study delves deep into nanosecond-scale protein dynamics of a structurally and functionally conserved protein across species separated by almost a billion years, investigating ten homologs in complex with their ligand. By inducing a photo-triggered destabilization of the ligand inside the binding pocket, we resolved distinct kinetic footprints for each homolog via transient infrared spectroscopy. Strikingly, we found a cascade of rearrangements within the protein complex which manifest in time points of increased dynamic activity conserved over hundreds of millions of years within a narrow window. Among these processes, one displays a subtle temporal shift correlating with evolutionary divergence, suggesting reduced selective pressure in the past. Our study not only uncovers the impact of evolution on molecular processes in a specific case, but has also the potential to initiate a field of scientific inquiry within molecular paleontology, where species are compared and classified based on the rapid pace of protein dynamic processes; a field which connects the shortest conceivable time scale in living matter (10[Formula: see text] s) with the largest ones (10[Formula: see text] s)., Competing Interests: Competing interests statement:The authors declare no competing interest.

Several ways to electronically synchronize different types of amplified femtosecond laser systems are presented based on a single freely programmable electronics hardware: arbitrary-detuning asynchronous optical sampling (ADASOPS), as well as actively locking two femtosecond laser oscillators, albeit not necessarily to the same round-trip frequency. They allow us to rapidly probe a very wide range of timescales, from picoseconds to potentially seconds, in a single transient absorption experiment without the need to move any delay stage. Experiments become possible that address a largely unexplored aspect of many photochemical reactions, in particular in the context of photo-catalysis as well as photoactive proteins, where an initial femtosecond trigger very often initiates a long-lasting cascade of follow-up processes. The approach is very versatile and allows us to synchronize very different lasers, such as a Ti:Sa amplifier and a 100 kHz Yb-laser system. The jitter of the synchronization, and therewith the time-resolution in the transient experiment, lies in the range from 1 to 3 ps, depending on the method. For illustration, transient IR measurements of the excited state solvation and decay of a metal carbonyl complex as well as the full reaction cycle of bacteriorhodopsin are shown. The pros and cons of the various methods are discussed, with regard to the scientific question one might want to address, and also with regard to the laser systems that might be already existent in a laser lab.

ConspectusUltrafast spectroscopy and imaging have become tools utilized by a broad range of scientists involved in materials, energy, biological, and chemical sciences. Commercialization of ultrafast spectrometers including transient absorption spectrometers, vibrational sum frequency generation spectrometers, and even multidimensional spectrometers have put these advanced spectroscopy measurements into the hands of practitioners originally outside the field of ultrafast spectroscopy. There is now a technology shift occurring in ultrafast spectroscopy, made possible by new Yb-based lasers, that is opening exciting new experiments in the chemical and physical sciences. Amplified Yb-based lasers are not only more compact and efficient than their predecessors but also, most importantly, operate at many times the repetition rate with improved noise characteristics in comparison to the previous generation of Ti:sapphire amplifier technologies. Taken together, these attributes are enabling new experiments, generating improvements to long-standing techniques, and affording the transformation of spectroscopies to microscopies. This Account aims to show that the shift to 100 kHz lasers is a transformative step in nonlinear spectroscopy and imaging, much like the dramatic expansion that occurred with the commercialization of Ti:sapphire laser systems in the 1990s. The impact of this technology will be felt across a great swath of scientific communities. We first describe the technology landscape of amplified Yb-based laser systems used in conjunction with 100 kHz spectrometers operating with shot-to-shot pulse shaping and detection. We also identify the range of different parametric conversion and supercontinuum techniques which now provide a path to making pulses of light optimal for ultrafast spectroscopy. Second, we describe specific instances from our laboratories of how the amplified Yb-based light sources and spectrometers are transformative. For multiple probe time-resolved infrared and transient 2D IR spectroscopy, the gain in temporal span and signal-to-noise enables dynamical spectroscopy measurements from femtoseconds to seconds. These gains widen the applicability of time-resolved infrared techniques across a range of topics in photochemistry, photocatalysis, and photobiology as well as lower the technical barriers to implementation in a laboratory. For 2D visible spectroscopy and microscopy with white light, as well as 2D IR imaging, the high repetition rates of these new Yb-based light sources allow one to spatially map 2D spectra while maintaining high signal-to-noise in the data. To illustrate the gains, we provide examples of imaging applications in the study of photovoltaic materials and spectroelectrochemistry.

Cyanobacteriochromes are bistable photoreceptor proteins with desirable photochemical properties for biotechnological applications, such as optogenetics or fluorescence microscopy. Here, we investigate Slr1393-g3, a cyanobacteriochrome that reversibly photoswitches between a red-absorbing (Pr) and green-absorbing (Pg) form. We applied advanced IR spectroscopic methods to track the sequence of intermediates during the photocycle over many orders of magnitude in time. In the conversion from Pg to Pr, we have revealed a new intermediate with distinct spectroscopic features in the IR, which precedes Pr formation using transient IR spectroscopy. In addition, stationary and transient 2D IR experiments measured the vibrational couplings between different groups of the chromophore and the protein in these intermediate states, as well as their structural disorder. Anharmonic QM/MM calculations predict spectra in good agreement with experimental 2D IR spectra of the initial and final states of the photocycle. They facilitate the assignment of the IR spectra that serve as a basis for the interpretation of the spectroscopic results and suggest structural changes of the intermediates along the photocycle.

The allosteric protein MCL-1 and its natural inhibitors, the BH3-only proteins PUMA, BIM, and NOXA regulate apoptosis by interacting promiscuously within an entangled binding network. Little is known about the transient processes and dynamic conformational fluctuations that are the basis for the formation and stability of the MCL-1/BH3-only complex. In this study, we designed photoswitchable versions of MCL-1/PUMA and MCL-1/NOXA, and investigated the protein response after an ultrafast photo-perturbation with transient infrared spectroscopy. We observed partial α-helical unfolding in all cases, albeit on strongly varying timescales (1.6 ns for PUMA, 9.7 ns for the previously studied BIM, and 85 ns for NOXA). These differences are interpreted as a BH3-only-specific "structural resilience" to defy the perturbation while remaining in MCL-1's binding pocket. Thus, the presented insights could help to better understand the differences between PUMA, BIM, and NOXA, the promiscuity of MCL-1, in general, and the role of the proteins in the apoptotic network.

Theoretical calculations of the low-frequency anharmonic couplings of the β-phase of crystalline bromoform are presented based on density functional theory quantum chemistry calculations. The electrical and mechanical anharmonicities between intra- and intermolecular modes are calculated, revealing that the electrical anharmonicity dominates the cross-peak intensities in the 2D Raman-THz response and crystalline, as well as liquid, bromoform. Furthermore, the experimentally observed difference in relative cross-peak intensities between the two intramolecular modes of bromoform and the intermolecular modes can be explained by the C 3v -symmetry of bromoform in combination with orientational averaging. The good agreement with the experimental results provides further evidence for our interpretation that the 2D Raman-THz response of bromoform is, indeed, related to the anharmonic coupling between the intra- and intermolecular modes.

Halogenated groups are relevant in pharmaceutical applications and potentially useful spectroscopic probes for infrared spectroscopy. In this work, the structural dynamics and infrared spectroscopy of para -fluorophenol (F-PhOH) and phenol (PhOH) is investigated in the gas phase and in water using a combination of experiment and molecular dynamics (MD) simulations. The gas phase and solvent dynamics around F-PhOH and PhOH is characterized from atomistic simulations using empirical energy functions with point charges or multipoles for the electrostatics, Machine Learning (ML) based parametrizations and with full ab initio (QM) and mixed Quantum Mechanical/Molecular Mechanics (QM/MM) simulations with a particular focus on the CF- and OH-stretch region. The CF-stretch band is heavily mixed with other modes whereas the OH-stretch in solution displays a characteristic high-frequency peak around 3600 cm -1 most likely associated with the -OH group of PhOH and F-PhOH together with a characteristic progression below 3000 cm -1 due to coupling with water modes which is also reproduced by several of the simulations. Solvent and radial distribution functions indicate that the CF-site is largely hydrophobic except for simulations using point charges which renders them unsuited for correctly describing hydration and dynamics around fluorinated sites. The hydrophobic character of the CF-group is particularly relevant for applications in pharmaceutical chemistry with a focus on local hydration and interaction with the surrounding protein.

The protein MCL-1 is a crucial factor in regulating apoptosis, the programmed cell death, and thus plays a major role in numerous cancer types. The allosteric protein MCL-1 is naturally moderated by the BH3-only peptide BIM, which binds at its canonical binding groove. In its isolated form, BIM is disordered but assumes an α-helical shape when bound by MCL-1. The underlying binding mechanism (i.e., induced fit vs conformational selection), as well as the time scales of the signal cascade subsequent to binding, are not understood. Here, an artificially photoswitchable variant of the MCL-1/BIM complex was designed and investigated by transient infrared spectroscopy. By destabilizing the α-helix of BIM with a covalently linked azobenzene photoswitch, the dynamical response of the whole complex upon an ultrafast photo-perturbation was characterized. While the destabilized and partially unfolded BIM still binds to MCL-1, a step-like cascade of structural rearrangements of both, MCL-1 and BIM was detected, spanning a wide range of time scales from pico- to microseconds. The results indicate that BIM binds according to an induced fit mechanism, while the structural adaptations of MCL-1 may constitute an allosteric signal., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

The spectroscopy and structural dynamics of a deep eutectic mixture (KSCN/acetamide) with varying water content is investigated from 2D IR (with the C-N stretch vibration of the SCN - anions as the reporter) and THz spectroscopy. Molecular dynamics simulations correctly describe the nontrivial dependence of both spectroscopic signatures depending on water content. For the 2D IR spectra, the MD simulations relate the steep increase in the cross-relaxation rate at high water content to the parallel alignment of packed SCN - anions. Conversely, the nonlinear increase of the THz absorption with increasing water content is mainly attributed to the formation of larger water clusters. The results demonstrate that a combination of structure-sensitive spectroscopies and molecular dynamics simulations provides molecular-level insights into the emergence of heterogeneity of such mixtures by modulating their composition.

Phytochromes are ubiquitous photoreceptor proteins that undergo a significant refolding of secondary structure in response to initial photoisomerization of the chromophoric group. This process is important for the signal transduction through the protein and thus its regulatory function in different organisms. Here, we employ two-dimensional infrared absorption (2D-IR) spectroscopy, an ultrafast spectroscopic technique that is sensitive to vibrational couplings, to study the photoreaction of bacterial phytochrome Agp1. By calculating difference spectra with respect to the photoactivation, we are able to isolate sharp difference cross-peaks that report on local changes in vibrational couplings between different sites of the chromophore and the protein. These results indicate inter alia that a dipole coupling between the chromophore and the so-called tongue region plays a role in stabilizing the protein in the light-activated state.

Two-dimensional (2D) Raman-THz spectroscopy in the frequency of up to 7 THz has been applied to study the crystalline β-phase of bromoform (CHBr 3 ). As for liquid CHBr 3 , cross peaks are observed, which, however, sharpen up in the crystalline sample and split into assignable sub-contributions. In the Raman dimension, the frequency positions of these cross peaks coincide with the intramolecular bending modes of the CHBr 3 molecules and in the THz dimension with the IR-active lattice modes of the crystal. This work expands the applicability of this new 2D spectroscopic technique to solid samples at cryogenic temperatures. Furthermore, it provides new experimental evidence that the cross peaks, indeed, originate from the coupling between intra- and intermolecular vibrational modes.

Controlling the activity of proteins with azobenzene photoswitches is a potent tool for manipulating their biological function. With the help of light, it is possible to change binding affinities, control allostery or manipulate complex biological processes, for example. Additionally, owing to their intrinsically fast photoisomerization, azobenzene photoswitches can serve as triggers that initiate out-of-equilibrium processes. Such switching of the activity initiates a cascade of conformational events that can be accessed with time-resolved methods. In this Review, we show how the potency of azobenzene photoswitching can be combined with transient spectroscopic techniques to disclose the order of events and experimentally observe biomolecular interactions in real time. This strategy will further our understanding of how a protein can accommodate, adapt and readjust its structure to answer an incoming signal, revealing more of the dynamical character of proteins., (© 2021. Springer Nature Limited.)

A stop-flow sample delivery system for transient spectroscopy is presented, which is, in particular, suited for laser-based instruments (quantum-cascade lasers or amplified femtosecond lasers) with excitation pulse repetition rates in the range 10-100 Hz. Two pulsing micro-valves are mounted onto a flow cuvette designed for transient IR spectroscopy, which is integrated into a flow cycle driven by a peristaltic pump. The performance of the system is demonstrated with transient IR experiments of the trans-to-cis photoisomerization of a water-soluble azobenzene derivative. The sample stands still when the micro-valves are closed and is pushed out from the probe beam focus on a 1 ms timescale when opening the micro-valves. The setup is extremely sample efficient. It needs only small sample volumes, and at the same time, it enables excitation of a large fraction of molecules in solution.

We present a 2D-Raman-terahertz (THz) setup with multichannel (single-shot) THz detection, utilizing two crossed echelons, in order to reduce the acquisition time of typical 2D-Raman-THz experiments from days to a few hours. This speed-up is obtained in combination with a high repetition rate (100 kHz) Yb-based femtosecond laser system and a correspondingly fast array detector. The wavelength of the Yb-laser (1030 nm) is advantageous, since it assures almost perfect phase matching in GaP for THz generation and detection and since the dispersion in the transmissive echelons is minimal. 2D-Raman-THz test measurements on liquid bromoform (CHBr 3 ) are reported. An enhancement of a factor ∼5.8 in signal-to-noise ratio is obtained for single-shot detection when compared to conventional step-scanning measurements in the THz time domain, corresponding to a speed-up of acquisition time of ∼34.

The dynamics of peptide-protein binding and unbinding of a variant of the RNase S system has been investigated. To initiate the process, a photoswitchable azobenzene moiety has been covalently linked to the S-peptide, thereby switching its binding affinity to the S-protein. Transient fluorescence quenching was measured with the help of a time-resolved fluorometer, which has been specifically designed for these experiments and is based on inexpensive light-emitting diodes and laser diodes only. One mutant shows on-off behavior with no specific binding detectable in one of the states of the photoswitch. Unbinding is faster by at least 2 orders of magnitude, compared to that of other variants of the RNase S system. We conclude that unbinding is essentially barrier-less in that case, revealing the intrinsic dynamics of the unbinding event, which occurs on a time scale of a few hundred microseconds in a strongly stretched-exponential manner.

The phototriggered unbinding of the intrinsically disordered S-peptide from the RNase S complex is studied with the help of transient IR spectroscopy, covering a wide range of time scales from 100 ps to 10 ms. To that end, an azobenzene moiety has been linked to the S-peptide in a way that its helicity is disrupted by light, thereby initiating its complete unbinding. The full sequence of events is observed, starting from unfolding of the helical structure of the S-peptide on a 20 ns time scale while still being in the binding pocket of the S-protein, S-peptide unbinding after 300 μs, and the structural response of the S-protein after 3 ms. With regard to the S-peptide dynamics, the binding mechanism can be classified as an induced fit, while the structural response of the S-protein is better described as conformational selection.

While much is known about different allosteric regulation mechanisms, the nature of the allosteric signal and the time scale on which it propagates remains elusive. The PDZ3 domain from postsynaptic density-95 protein is a small protein domain with a terminal third α-helix, i.e., the α3-helix, which is known to be allosterically active. By cross-linking the allosteric helix with an azobenzene moiety, we obtained a photocontrollable PDZ3 variant. Photoswitching triggers its allosteric transition, resulting in a change in binding affinity of a peptide to the remote binding pocket. Using time-resolved infrared and UV/vis spectroscopy, we follow the allosteric signal transduction and reconstruct the timeline in which the allosteric signal propagates through the protein within 200 ns.

Recently, re-purposing of cyanobacterial photoreceptors as optogentic actuators enabled light-regulated protein expression in different host systems. These new bi-stable optogenetic tools enable interesting new applications, but their light-driven working mechanism remains largely elusive on a molecular level. Here, we study the optogenetic cyanobacteriochrome Am1-c0023g2 with isotope labeling and two dimensional infrared (2D-IR) spectroscopy. Isotope labeling allows us to isolate two site-specific carbonyl marker modes from the overwhelming mid-IR signal of the peptide backbone vibrations. Unlike conventional difference-FTIR spectroscopy, 2D-IR is sensitive to homogeneous and inhomogeneous broadening mechanisms of these two vibrational probes in the different photostates of the protein. We analyse the 2D-IR line shapes in the context of available structural models and find that they reflect the hydrogen-bonding environment of these two marker groups.

This review aims to promote the role of transient IR spectroscopy to investigate molecular-based photocatalytic water reduction. Examples are discussed in which this method has been successfully applied to elucidate reaction mechanisms. Focus is given to kinetic changes and their consequences when a photochemical water reduction system, which is functional and well understood in solution, is brought onto a metal oxide surface.

A new application of high-repetition rate, femtosecond Yb-laser/amplifier systems is introduced: transient 2D IR spectroscopy covering the time range from micro- to milliseconds. This approach intertwines the measurement of 2D IR spectra with the time separation from an actinic pump pulse and utilizes the high repetition rate of these lasers systems in two ways: by offering a high time resolution (10 µs) and by enabling the measurement of many 2D IR spectra. The well-studied photocycle of bacteriorhodopsin is used as a demonstration object in this proof-of-principle experiment.

Isotopic substitutions largely affect the dielectric relaxation dynamics of hydrogen-bonded liquid water; yet, the role of the altered molecular masses and nuclear quantum effects has not been fully established. To disentangle these two effects we study the dielectric relaxation of light (H216O), heavy (D216O) and heavy-oxygen (H218O) water at temperatures ranging from 278 to 338 K. Upon 16O/18O exchange, we find that the relaxation time of the collective orientational relaxation mode of water increases by 4-5%, in quantitative agreement with the enhancement of viscosity. Despite the rotational character of dielectric relaxation, the increase is consistent with a translational mass factor. For H/D substitution, the slow-down of the relaxation time is more pronounced and also shows a strong temperature dependence. In addition to the classical mass factor, the enhancement of the relaxation time for D216O can be described by an apparent temperature shift of 7.2 K relative to H216O, which is higher than the 6.5 K shift reported for viscosity. As this shift accounts for altered zero-point energies, the comparison suggests that the underlying thermally populated states relevant to the activation of viscous flow and dielectric relaxation differ.

Resonant and off-resonant mid-infrared pump-probe spectroscopy is used to measure the vibrational dynamics of CO adsorbed to thin (0.2 nm, 2 nm, and 10 nm) heterogeneous Pt layers in an aqueous solution. The transient signals observed with resonant pumping are dominated by vibrational relaxation of the CO internal stretch vibration with a lifetime of T 1 ∼ 3 ps in all cases. Off-resonant pumping suppresses that contribution to the signal and singles out a signal, which is attributed to heating of the metal layer as well as transient desorption of the CO molecules. Due to the small photon energy (0.2 eV) used as pump pulses, the mechanism of desorption must be thermal, in which case the desorption yield depends exclusively on the fluence of absorbed light and not its wavelength. The thin Pt layers facilitate CO desorption, despite a relatively low pump pulse fluence, as they concentrate the absorbed energy in a small volume.

This work explores the solid-liquid interface of a rhenium-tricarbonyl complex embedded in a layer of zirconium oxide deposited by atomic layer deposition (ALD). Time-resolved and steady state infrared spectroscopy were applied to reveal the correlations between the thickness of the ALD layer and the spectroscopic response of the system. We observed a transition of the molecular environment from flexible to rigid, as well as limitations to ligand exchange and excited state quenching on the embedded complexes, when the ALD layer is roughly of the same height as the molecules.

Over the last decades, photoreceptive proteins were extensively studied with biophysical methods to gain a fundamental understanding of their working mechanisms and further guide the development of optogenetic tools. Time-resolved infrared (IR) spectroscopy is one of the key methods to access their functional non-equilibrium processes with high temporal resolution but has the major drawback that experimental data are usually highly complex. Linking the spectral response to specific molecular events is a major obstacle. Here, we investigate a cyanobacteriochrome photoreceptor with a combined approach of transient absorption spectroscopy in the visible and IR spectral regions. We obtain kinetic information in both spectral regions by analysis with two different fitting methods: global multiexponential fitting and lifetime analysis. We investigate the ground state dynamics that follow photoexcitation in both directions of the bi-stable photocycle (Pr * and Pg * ) in the nanosecond and microsecond time regimes. We find two ground state intermediates associated with the decay of Pr * and four with Pg * and report the macroscopic time constants of their interconversions. One of these processes is assigned to a structural change in the protein backbone.

Allosteric regulation is an innate control in most metabolic and signalling cascades that enables living organisms to adapt to the changing environment by tuning the affinity and regulating the activity of target proteins. For a microscopic understanding of this process, a protein system has been designed in such a way that allosteric communication between the binding and allosteric site can be observed in both directions. To that end, an azobenzene-derived photoswitch has been linked to the α3-helix of the PDZ3 domain, arguably the smallest allosteric protein with a clearly identifiable binding and allosteric site. Photo-induced trans-to-cis isomerisation of the photoswitch increases the binding affinity of a small peptide ligand to the protein up to 120-fold, depending on temperature. At the same time, ligand binding speeds up the thermal cis-to-trans back-isomerisation rate of the photoswitch. Based on the energetics of the four states of the system (cis vs trans and ligand-bound vs free), the concept of an allosteric force is introduced, which can be used to drive chemical reactions.

The majority of 2D IR spectrometers operate at 1-10 kHz using Ti:Sapphire laser technology. We report a 2D IR spectrometer designed around Yb:KGW laser technology that operates shot-to-shot at 100 kHz. It includes a home-built OPA, a mid-IR pulse shaper, and custom-designed electronics with optional on-chip processing. We report a direct comparison between Yb:KGW and Ti:Sapphire based 2D IR spectrometers. Even though the mid-IR pulse energy is much lower for the Yb:KGW driven system, there is an 8x improvement in signal-to-noise over the 1 kHz Ti:Sapphire driven spectrometer to which it is compared. Experimental data is shown for sub-millimolar concentrations of amides. Advantages and disadvantages of the design are discussed, including thermal background that arises at high repetition rates. This fundamental spectrometer design takes advantage of newly available Yb laser technology in a new way, providing a straightforward means of enhancing sensitivity.

Vibrational energy transfer (VET) between two isotopologues of [Re(dcb)(CO) 3 Br] immobilized on a TiO 2 surface is studied with the help of 2D IR spectroscopy in dependence of surface coverage. To dilute the molecules on the surface, and thereby control the intermolecular distances, two different diluents have been used: a third isotopologue of the same molecule and 4-cyanobenzoic acid. As expected, the VET rate decreases with dilution. For a quantitative investigation of the distance dependence of the VET rate, we analyze the data based on an excitonic model. This model reveals the typical 1/r 6 -distance dependence for a dimer of a donor and acceptor, similar to the nuclear Overhauser effect in NMR spectroscopy or Förster resonant energy transfer in electronic spectroscopy. However, VET becomes a collective phenomenon on the surface, with the existence of a network of coupled molecules and its disappearance below a percolation threshold, dominating the concentration dependence of the VET rate.

While allostery is of paramount importance for protein regulation, the underlying dynamical process of ligand (un)binding at one site, resulting time evolution of the protein structure, and change of the binding affinity at a remote site are not well understood. Here the ligand-induced conformational transition in a widely studied model system of allostery, the PDZ2 domain, is investigated by transient infrared spectroscopy accompanied by molecular dynamics simulations. To this end, an azobenzene-derived photoswitch is linked to a peptide ligand in a way that its binding affinity to the PDZ2 domain changes upon switching, thus initiating an allosteric transition in the PDZ2 domain protein. The subsequent response of the protein, covering four decades of time, ranging from ∼1 ns to ∼μs, can be rationalized by a remodeling of its rugged free-energy landscape, with very subtle shifts in the populations of a small number of structurally well-defined states. It is proposed that structurally and dynamically driven allostery, often discussed as limiting scenarios of allosteric communication, actually go hand-in-hand, allowing the protein to adapt its free-energy landscape to incoming signals., Competing Interests: The authors declare no competing interest.

The asymmetry of Fano line shapes observed for metal-adsorbate systems is reflected in two-dimensional infrared (2D IR) spectroscopy as a distorted spectrum. A phenomenological correction scheme is proposed that transforms distorted 2D IR spectra into conventional spectra. To that end, a phase correction factor is first derived from the IR absorption spectrum of the sample by symmetrizing the asymmetric line shape and subsequently applied to the distorted 2D IR spectra. The concept is illustrated for a model system consisting of an organic molecule ( p -mercaptobenzonitrile) adsorbed on a sputter-coated metal layer (Au). The correction scheme reveals conventional, easily interpretable 2D IR spectra.

The 2D-Raman-THz response in all possible time-orderings (Raman-THz-THz, THz-Raman-THz, and THz-THz-Raman) of amorphous water ice is calculated in two ways: from atomistic molecular dynamics simulations and with the help of a Feynman diagram model, the latter of which power-expands the potential energy surface and the dipole and polarizability surfaces up to leading order. Comparing both results allows one to dissect the 2D-Raman-THz response into contributions from mechanical anharmonicity, as well as electrical dipole and polarizability anharmonicities. Mechanical anharmonicity dominates the 2D-Raman-THz response of the hydrogen-bond stretching and hydrogen-bond bending bands of water, and dipole anharmonicity dominates that of the librational band, while the contribution of polarizability anharmonicity is comparably weak. A distinct echo of the hydrogen-bond stretching band is observed for the THz-Raman-THz pulse sequence, again dominated by mechanical anharmonicity. A peculiar mechanism is discussed, which is based on the coupling between the many normal modes within the hydrogen-bond stretching band and which will inevitably generate such an echo for an amorphous structure.

Hydridocarbonyl complexes, a class of industrially relevant catalysts, contain both the M-H and M-CO moieties. Here, using two-dimensional infrared spectroscopy, we examine the coupling of the typically weak M-H stretching mode and the intense M(C≡O) mode. By studying a series of Ir(I)- and Ir(III)-based hydridocarbonyl complexes, we show that the arrangement of the H and CO ligands in a trans configuration leads to strong vibrational coupling and mode delocalization. In contrast, a cis arrangement leads to no coupling, with the localized M-H mode having a much larger anharmonicity. These results highlight a promising strategy for enhancing the M-H vibration by intensity borrowing from the strong C≡O modes in a trans configuration, allowing for direct determination by infrared spectroscopy of both the oxidation (by frequency shifts) and the protonation state (via vibrational coupling) of the complex, in mechanistic studies of proton-coupled electron transfer reactions.

Hybrid two-dimensional (2D) Raman-terahertz spectroscopy with the Raman-terahertz-terahertz (RTT) pulse sequence is used to explore the ultrafast intra- and intermolecular degrees of freedom of liquid bromoform (CHBr 3 ) in the frequency range of 1-8 THz. Cross peaks observed in these 2D spectra are assigned to the coupling between the narrow intramolecular modes of the molecules and the much broader intermolecular degrees of freedom of the liquid. This assignment is based on the frequency position of the cross peaks; however, it is shown that these frequency positions can be deduced accurately only when properly taking into account the convolution of the molecular response with the instrument response function of the experimental setup, the latter of which distorts the 2D spectra considerably. The assignment is supported by additional experiments on diiodomethane (CH 2 I 2 ), which has only one intramolecular mode in the frequency range of the experiment, and hence excludes the possibility of intramolecular couplings.

An azobenzene-derived photoswitch has been covalently cross-linked to two sites of the S-peptide in the RNase S complex in a manner that the α-helical content of the S-peptide reduces upon cis -to- trans isomerization of the photoswitch. Three complementary experimental techniques have been employed, isothermal titration calorimetry, circular dichroism spectroscopy and intrinsic tyrosine fluorescence quenching, to determine the binding affinity of the S-peptide to the S-protein in the two states of the photoswitch. Five mutants with the photoswitch attached to different sites of the S-peptide have been explored, with the goal to maximize the change in binding affinity upon photoswitching, and to identify the mechanisms that determine the binding affinity. With regard to the first goal, one mutant has been identified, which binds with reasonable affinity in the one state of the photoswitch, while specific binding is completely switched off in the other state. With regard to the second goal, accompanying molecular dynamics simulations combined with a quantitative structure activity relationship revealed that the α-helicity of the S-peptide in the binding pocket correlates surprisingly well with measured dissociation constants. Moreover, the simulations show that both configurations of all S-peptides exhibit quite well-defined structures, even in apparently disordered states.

Allosteric regulation plays an important role in many biological processes, such as signal transduction, transcriptional regulation, and metabolism. Allostery is rooted in the fundamental physical properties of macromolecular systems, but its underlying mechanisms are still poorly understood. A collection of contributions to a recent interdisciplinary CECAM (Center Européen de Calcul Atomique et Moléculaire) workshop is used here to provide an overview of the progress and remaining limitations in the understanding of the mechanistic foundations of allostery gained from computational and experimental analyses of real protein systems and model systems. The main conceptual frameworks instrumental in driving the field are discussed. We illustrate the role of these frameworks in illuminating molecular mechanisms and explaining cellular processes, and describe some of their promising practical applications in engineering molecular sensors and informing drug design efforts., (Copyright © 2019. Published by Elsevier Ltd.)

The 2D Raman-terahertz (THz) response of liquid water is studied in dependence of temperature and isotope substitution ([Formula: see text]O, [Formula: see text]O, and [Formula: see text]O). In either case, a very short-lived (i.e., between 75 and 95 fs) echo is observed that reports on the inhomogeneity of the low-frequency intermolecular modes and hence, on the heterogeneity of the hydrogen bond networks of water. The echo lifetime slows down by about 20% when cooling the liquid from room temperature to the freezing point. Furthermore, the echo lifetime of [Formula: see text]O is [Formula: see text] slower than that of [Formula: see text]O, and both can be mapped on each other by introducing an effective temperature shift of [Formula: see text] K. In contrast, the temperature-dependent echo lifetimes of [Formula: see text]O and [Formula: see text]O are the same within error. [Formula: see text]O and [Formula: see text]O have identical masses, yet [Formula: see text]O is much closer to [Formula: see text]O in terms of nuclear quantum effects. It is, therefore, concluded that the echo is a measure of the structural inhomogeneity of liquid water induced by nuclear quantum effects., Competing Interests: The authors declare no conflict of interest.

2D-Raman-THz spectroscopy of liquid water, which has been presented recently [J. Savolainen et al., Proc. Natl. Acad. Sci. U. S. A. 110, 20402 (2013)], directly probes the intermolecular degrees of freedom of the hydrogen-bond network. However, being a relatively new technique, its information content is not fully explored to date. While the spectroscopic signal can be simulated based on molecular dynamics simulation in connection with a water force field, it is difficult to relate spectroscopic signatures to the underlying microscopic features of the force field. Here, a completely different approach is taken that starts from an as simple as possible model, i.e., a single vibrational mode with electrical and mechanical anharmonicity augmented with homogeneous and inhomogeneous broadening. An intuitive Feynman diagram picture is developed for all possible pulse sequences of hybrid 2D-Raman-THz spectroscopy. It is shown that the model can explain the experimental data essentially quantitatively with a very small set of parameters, and it is tentatively concluded that the experimental signal originates from the hydrogen-bond stretching vibration around 170 cm -1 . Furthermore, the echo observed in the experimental data can be quantified by fitting the model. A dominant fraction of its linewidth is attributed to quasi-inhomogeneous broadening in the slow-modulation limit with a correlation time of 370 fs, reflecting the lifetime of the hydrogen-bond networks giving rise to the absorption band.

Three structurally very different antifreeze proteins (AFPs) are studied, addressing the question as to what extent the hypothesized preordering-binding mechanism is still relevant in the second solvation layer of the protein and beyond. Assuming a two-state model of water, the solvation layers are analyzed with the help of molecular dynamics simulations together with a Markov state model, which investigates the local tedrahedrality of the water hydrogen-bond network around a given water molecule. It has been shown previously that this analysis can discriminate the high-entropy, high-density state of the liquid (HDL) from its more structured low-density state (LDL). All investigated proteins, regardless of whether they are an AFP or not, have a tendency to increase the amount of HDL in their second solvation layer. The ice binding site (IBS) of the antifreeze proteins counteracts that trend, with either a hole in the HDL layer or a true excess of LDL. The results correlate to a certain extent with recent experiments, which have observed ice-like vibrational (VSFG) spectra for the water atop the IBS of only a subset of antifreeze proteins. It is concluded that the preordering-binding mechanism indeed seems to play a role but is only part of the overall picture.

The noncanonical amino acid azidohomoalanine (Aha) is known to be an environment-sensitive infrared probe for the site-specific investigation of protein structure and dynamics. Here, the capability of that label is explored to detect protein-ligand interactions by incorporating it in the vicinity of the binding groove of a PDZ2 domain. Circular dichroism and isothermal titration calorimetry measurements reveal that the perturbation of the protein system by mutation is negligible, with minimal influence on protein stability and binding affinity. Two-dimensional infrared spectra exhibit small (1-3 cm -1 ) but clearly measurable red shifts of the Aha vibrational frequency upon binding of two different peptide ligands, while accompanying molecular dynamics simulations suggest that these red shifts are induced by polar contacts with side chains of the peptide ligands. Hence, Aha is a versatile and minimally invasive vibrational label that is not only able to report on large structural changes during, e.g., protein folding, but also on very subtle changes of the electrostatic environment upon ligand binding.

A polarizable water model is presented which has been calibrated against experimental THz and Raman spectra of bulk water. These low-frequency spectra directly probe the dynamics, and thereby intermolecular interactions, on time scales relevant to molecular motions. The model is based on the TL4P force field developed recently by Tavan and co-workers [J. Phys. Chem. B 117, 9486 (2013)], which has been designed to be transferable between different environments; in particular, to correctly describe the electrostatic properties of both the isolated water molecule in the gas-phase and the liquid water at ambient conditions. Following this design philosophy, TL4P was amended with charge transfer across hydrogen-bonded dimers as well as an anisotropic polarizability in order to correctly reproduce the THz and Raman spectra. The thermodynamic and structural properties of the new model are of equal quality as those of TL4P, and at the same time, an almost quantitative agreement with the spectroscopic data could be achieved. Since TL4P is a rigid model with a single polarizable site, it is computationally very efficient, while the numerical overhead for the addition of charge transfer and the anisotropic polarizability is minor. Overall, the model is expected to be well suited for, e.g., large scale simulations of 2D-Raman-THz spectra or biomolecular simulations.

The response of water re-solvating a charge-transfer dye (deprotonated Coumarin 343) after photoexcitation has been measured by means of transient THz spectroscopy. Two steps of increasing THz absorption are observed, a first ∼10 ps step on the time scale of Debye relaxation of bulk water and a much slower step on a 3.9 ns time scale, the latter of which reflecting heating of the bulk solution upon electronic relaxation of the dye molecules from the S 1 back into the S 0 state. As an additional reference experiment, the hydroxyl vibration of water has been excited directly by a short IR pulse, establishing that the THz signal measures an elevated temperature within ∼1 ps. This result shows that the first step upon dye excitation (10 ps) is not limited by the response time of the THz signal; it rather reflects the reorientation of water molecules in the solvation layer. The apparent discrepancy between the relatively slow reorientation time and the general notion that water is among the fastest solvents with a solvation time in the sub-picosecond regime is discussed. Furthermore, non-equilibrium molecular dynamics simulations have been performed, revealing a close-to-quantitative agreement with experiment, which allows one to disentangle the contribution of heating to the overall THz response from that of water orientation.

While the theory of protein folding is well developed, including concepts such as rugged energy landscape, folding funnel, etc., the same degree of understanding has not been reached for the description of the dynamics of allosteric transitions in proteins. This is not only due to the small size of the structural change upon ligand binding to an allosteric site, but also due to challenges in designing experiments that directly observe such an allosteric transition. On the basis of recent pump-probe-type experiments (Buchli et al. 2013 Proc. Natl Acad. Sci. USA 110 , 11 725-11 730. (doi:10.1073/pnas.1306323110)) and non-equilibrium molecular dynamics simulations (Buchenberg et al. 2017 Proc. Natl Acad. Sci. USA 114 , E6804-E6811. (doi:10.1073/pnas.1707694114)) studying an photoswitchable PDZ2 domain as model for an allosteric transition, we outline in this perspective how such a description of allosteric communication might look. That is, calculating the dynamical content of both experiment and simulation (which agree remarkably well with each other), we find that allosteric communication shares some properties with downhill folding, except that it is an 'order-order' transition. Discussing the multiscale and hierarchical features of the dynamics, the validity of linear response theory as well as the meaning of 'allosteric pathways', we conclude that non-equilibrium experiments and simulations are a promising way to study dynamical aspects of allostery.This article is part of a discussion meeting issue 'Allostery and molecular machines'., (© 2018 The Author(s).)

The vibrational dynamics of Pt-H on a nanostructured platinum surface has been examined by ultrafast infrared spectroscopy. Three bands are observed at 1800, 2000, and 2090 cm -1 , which are assigned to Pt-CO in a bridged and linear configuration and Pt-H, respectively. Lifetime analysis revealed a time constant of (0.8 ± 0.1) ps for the Pt-H mode, considerably shorter than that of Pt-CO because of its stronger coupling to the metal substrate. Two-dimensional attenuated total reflection infrared spectroscopy provided additional evidence for the assignment based on the anharmonic shift, which is large in the case of Pt-H (90 cm -1 ), in agreement with the density functional theory calculations. The absorption cross section of Pt-H is smaller than that of the very strong Pt-CO vibration by only a modest factor of ∼1.5-3. Because Pt-H is transiently involved in catalytic water splitting on Pt, the present spectroscopic characterization paves the way for in-operando kinetic studies of such reactions.

[This corrects the article DOI: 10.1063/1.4996448.].

Due to their very nature, ultrafast phenomena are often accompanied by the occurrence of nonadiabatic effects. From a theoretical perspective, the treatment of nonadiabatic processes makes it necessary to go beyond the (quasi) static picture provided by the time-independent Schrödinger equation within the Born-Oppenheimer approximation and to find ways to tackle instead the full time-dependent electronic and nuclear quantum problem. In this review, we give an overview of different nonadiabatic processes that manifest themselves in electronic and nuclear dynamics ranging from the nonadiabatic phenomena taking place during tunnel ionization of atoms in strong laser fields to the radiationless relaxation through conical intersections and the nonadiabatic coupling of vibrational modes and discuss the computational approaches that have been developed to describe such phenomena. These methods range from the full solution of the combined nuclear-electronic quantum problem to a hierarchy of semiclassical approaches and even purely classical frameworks. The power of these simulation tools is illustrated by representative applications and the direct confrontation with experimental measurements performed in the National Centre of Competence for Molecular Ultrafast Science and Technology.

Intermolecular vibrational energy transfer in monolayers of isotopically mixed rhenium carbonyl complexes at solid-liquid interfaces is investigated with the help of ultrafast 2D Attenuated Total Reflectance Infrared (2D ATR IR) spectroscopy in dependence of plasmonic surface enhancement effects. Dielectric and plasmonic materials are used to demonstrate that plasmonic effects have no impact on the vibrational energy transfer rate in a regime of moderate IR surface enhancement (enhancement factors up to ca. 30). This result can be explained with the common image-dipole picture. The vibrational energy transfer rate thus can be used as a direct observable to determine intermolecular distances on surfaces, regardless of their plasmonic properties.

The transfer of charge at the molecular level plays a fundamental role in many areas of chemistry, physics, biology and materials science. Today, more than 60 years after the seminal work of R. A. Marcus, charge transfer is still a very active field of research. An important recent impetus comes from the ability to resolve ever faster temporal events, down to the attosecond time scale. Such a high temporal resolution now offers the possibility to unravel the most elementary quantum dynamics of both electrons and nuclei that participate in the complex process of charge transfer. This review covers recent research that addresses the following questions. Can we reconstruct the migration of charge across a molecule on the atomic length and electronic time scales? Can we use strong laser fields to control charge migration? Can we temporally resolve and understand intramolecular charge transfer in dissociative ionization of small molecules, in transition-metal complexes and in conjugated polymers? Can we tailor molecular systems towards specific charge-transfer processes? What are the time scales of the elementary steps of charge transfer in liquids and nanoparticles? Important new insights into each of these topics, obtained from state-of-the-art ultrafast spectroscopy and/or theoretical methods, are summarized in this review.

This review provides a comprehensive overview of the structural dynamics in topical gas- and condensed-phase systems on multiple length and time scales. Starting from vibrationally induced dissociation of small molecules in the gas phase, the question of vibrational and internal energy redistribution through conformational dynamics is further developed by considering coupled electron/proton transfer in a model peptide over many orders of magnitude. The influence of the surrounding solvent is probed for electron transfer to the solvent in hydrated I - . Next, the dynamics of a modified PDZ domain over many time scales is analyzed following activation of a photoswitch. The hydration dynamics around halogenated amino acid side chains and their structural dynamics in proteins are relevant for iodinated TyrB26 insulin. Binding of nitric oxide to myoglobin is a process for which experimental and computational analyses have converged to a common view which connects rebinding time scales and the underlying dynamics. Finally, rhodopsin is a paradigmatic system for multiple length- and time-scale processes for which experimental and computational methods provide valuable insights into the functional dynamics. The systems discussed here highlight that for a comprehensive understanding of how structure, flexibility, energetics, and dynamics contribute to functional dynamics, experimental studies in multiple wavelength regions and computational studies including quantum, classical, and more coarse grained levels are required.

Recent years have seen dramatic developments in the technology of intense pulsed light sources in the THz frequency range. Since many dipole-active excitations in solids and molecules also lie in this range, there is now a tremendous potential to use these light sources to study linear and nonlinear dynamics in such systems. While several experimental investigations of THz-driven dynamics in solid-state systems have demonstrated a variety of interesting linear and nonlinear phenomena, comparatively few efforts have been made to drive analogous dynamics in molecular systems. In the present Perspective article, we discuss the similarities and differences between THz-driven dynamics in solid-state and molecular systems on both conceptual and practical levels. We also discuss the experimental parameters needed for these types of experiments and thereby provide design criteria for a further development of this new research branch. Finally, we present a few recent examples to illustrate the rich physics that may be learned from nonlinear THz excitations of phonons in solids as well as inter-molecular vibrations in liquid and gas-phase systems.