Search

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.

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.

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. [ABSTRACT FROM AUTHOR]

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 C3v-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. [ABSTRACT FROM AUTHOR]

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-9 s) with the largest ones (1016 s). [ABSTRACT FROM AUTHOR]

Two-dimensional (2D) Raman-THz spectroscopy in the frequency of up to 7 THz has been applied to study the crystalline β-phase of bromoform (CHBr3). As for liquid CHBr3, 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 CHBr3 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. [ABSTRACT FROM AUTHOR]

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 (CHBr3) 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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 T1 ∼ 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Vibrational energy transfer (VET) between two isotopologues of [Re(dcb)(CO)3Br] immobilized on a TiO2 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/r6-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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

A three-point velocity correlation function ⟨ v (t1 + t2) v 2 (t1) v (0)⟩ is introduced for a better understanding of the recent 2D-Raman-THz spectroscopy of the intermolecular degrees of freedoms of water and aqueous salt solutions. This correlation function reveals echoes in the presence of inhomogeneous broadening, which are coined "velocity echoes." In analogy to the well-known two-point velocity correlation function ⟨ v (t) v (0)⟩, it reflects the density of states (DOS) of the system under study without having to amend them with transition dipoles and transition polarizabilities. The correlation function can be calculated from equilibrium trajectories and converges extremely quickly. After deriving the theory, the information content of the three-point velocity correlation function is first tested based on a simple harmonic oscillator model with Langevin dynamics. Subsequently, velocity echoes of TIP4P/2005 water are calculated as a function of temperature, covering ambient conditions, the supercooled regime and amorphous ice, as well as upon addition of various salts. The experimentally observed trends can be reproduced qualitatively with the help of computationally very inexpensive molecular dynamics simulations. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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 S1 back into the S0 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. [ABSTRACT FROM AUTHOR]

We present experimental and simulated 2D IR spectra of some high-pressure forms of isotope-pure D2O ice and compare the results to those of ice Ih published previously [F. Perakis and P. Hamm, Phys. Chem. Chem. Phys. 14, 6250 (2012); L. Shi et al., ibid. 18, 3772 (2016)]. Ice II, ice V, and ice XIII have been chosen for this study, since this selection covers many aspects of the polymorphism of ice. That is, ice II is a hydrogen-ordered phase of ice, in contrast to ice Ih, while ice V and ice XIII are a hydrogen-disordered/ordered couple that shares essentially the same oxygen structure and hydrogenbonded network. For the transmission 2D IR spectroscopy, a novel method had to be developed for the preparation of ultrathin films (1-2 μm) of high-pressure ices with good optical quality. We also simulated 2D IR spectra based on molecular dynamics simulations connected to a vibrational exciton picture. These simulations agree with the experimental results in a semi-quantitative manner for ice II, while the same approach failed for ice V and ice XIII. From the perspective of 2D IR spectroscopy, ice II appears to be more inhomogeneously broadened than ice Ih, despite its hydrogen-order, which we attribute to the fact that ice II is structurally more complex with four distinguishable hydrogen bonds that mix due to exciton coupling. Ice V and ice XIII, on the other hand, behave as expected with the hydrogen-disordered case (ice V) being more inhomogenously broadened. Furthermore, in all hydrogen-ordered forms (ice II and ice XIII), cross peaks could be identified in the anisotropic 2D IR spectrum, whose signs reveal the relative direction of the corresponding excitonic states. [ABSTRACT FROM AUTHOR]

We have investigated the site-specific folding kinetics of a photo-switchable cross-linked [alpha]-helical peptide by using single [sup.13]C = [sup.18]O isotope labeling together with time-resolved IR spectroscopy. We observe that the folding times differ from site to site by a factor of eight at low temperatures (6[degrees]C), whereas at high temperatures (45[degrees]C), the spread is considerably smaller. The trivial sum of the site signals coincides with the overall folding signal of the unlabeled peptide, and different sites fold in a noncooperative manner. Moreover, one of the sites exhibits a decrease of hydrogen bonding upon folding, implying that the unfolded state at low temperature is not unstructured. Molecular dynamics simulations at low temperature reveal a stretched-exponential behavior which originates from parallel folding routes that start from a kinetically partitioned unfolded ensemble. Different metastable structures (i.e., traps) in the unfolded ensemble have a different ratio of loop and helical content. Control simulations of the peptide at high temperature, as well as without the cross-linker at low temperature, show faster and simpler (i.e., single-exponential) folding kinetics. The experimental and simulation results together provide strong evidence that the rate-limiting step in formation of a structurally constrained [alpha]-helix is the escape from heterogeneous traps rather than the nucleation rate. This conclusion has important implications for an [alpha]-helical segment within a protein, rather than an isolated [alpha]-helix, because the cross-linker is a structural constraint similar to those present during the folding of a globular protein. cooperativity | infrared spectroscopy | molecular dynamics simulation | peptide folding

Two photocleavable disulfide-bridged tetrapeptides [cyclo(Boc-Cys-Pro-Aib-Cys-OMe) and cyclo(Boc-Cys-Pro-Phe-Cys-OMe)] which provide insight into cyclic peptide dynamics are reported.

Author(s): Christoph Kolano (corresponding author) [1]; Jan Helbing [1]; Mariusz Kozinski [1]; Wolfram Sander [2]; Peter Hamm (corresponding author) [1] X-ray crystallography and nuclear magnetic resonance measurements provide us with [...]

Nuclear quantum effects in liquid water have profound implications for several of its macroscopic properties related to the structure, dynamics, spectroscopy, and transport. Although several of water's macroscopic properties can be reproduced by classical descriptions of the nuclei using interaction potentials effectively parameterized for a narrow range of its phase diagram, a proper account of the nuclear quantum effects is required to ensure that the underlying molecular interactions are transferable across a wide temperature range covering different regions of that diagram. When performing an analysis of the hydrogen-bonded structural networks in liquid water resulting from the classical (class) and quantum (qm) descriptions of the nuclei with two interaction potentials that are at the two opposite ends of the range in describing quantum effects, namely the flexible, pair-wise additive q-TIP4P/F, and the flexible, polarizable TTM3-F, we found that the (class) and (qm) results can be superimposed over the temperature range T = 250-350 K using a surprisingly simple, linear scaling of the two temperatures according to T(qm) = α T(class) + ΔT, where α = 0.99 and ΔT = -6 K for q-TIP4P/F and α = 1.24 and ΔT = -64 K for TTM3-F. This simple relationship suggests that the structural networks resulting from the quantum and classical treatment of the nuclei with those two very different interaction potentials are essentially similar to each other over this extended temperature range once a model-dependent linear temperature scaling lawis applied. [ABSTRACT FROM AUTHOR]