Implications of short time scale dynamics on long time processes

Solvent Dynamics of Aqueous Halides before and after Photoionization

Electron transfer reactions can be strongly influenced by solvent dynamics. We study the photoionization of halides in water as a model system for such reactions. There are no internal nuclear degrees of freedom in the solute, allowing the dynamics of the solvent to be uniquely identified. We simulate the equilibrium solvent dynamics for Cl^-, Br^-, I^-, and their respective neutral atoms in water, comparing quantum mechanical/molecular mechanical (QM/MM) and classical molecular dynamics (MD) methods. On the basis of the obtained configurations, we calculate the extended X-ray absorption fine structure (EXAFS) spectra rigorously based on the MD snapshots and compare them in detail with other theoretical and experimental results available in the literature. We find our EXAFS spectra based on QM/MM MD simulations in good agreement with their experimental counterparts for the ions. Classical MD simulations for the ions lead to EXAFS spectra that agree equally well with the experiment when it comes to the oscillatory period of the signal, even though they differ from the QM/MM radial distribution functions extracted from the MD. The amplitude is, however, considerably overestimated. This suggests that to judge the reliability of theoretical simulation methods or to elucidate fine details of the atomistic dynamics of the solvent based on EXAFS spectra, the amplitude as well as the oscillatory period need to be considered. If simulations fail qualitatively, as does the classical MD for the aqueous neutral halogen atoms, the resulting EXAFS will also be strongly affected in both oscillatory period and amplitude. The good reliability of QM/MM-based EXAFS simulations, together with clear qualitative differences in the EXAFS spectra found between halides and their atomic counterparts, suggests that a combined theory and experimental EXAFS approach is suitable for elucidating the nonequilibrium solvent dynamics in the photoionization of halides and possibly also for electron transfer reactions in more complex systems.

Conformational Heterogeneity and Frustration of the Tumor Suppressor p53 as Tuned by Punctual Mutations

The conformational heterogeneity of the p53 tumor suppressor, the wild-type (p53wt) and mutated forms, was investigated by a computational approach, including the modeling and all atoms of the molecular dynamics (MD) simulations. Four different punctual mutations (p53R175H, p53R248Q, p53R273H, and p53R282W) which are known to affect the DNA binding and belong to the most frequent hot-spot mutations in human cancers, were taken into consideration. The MD trajectories of the wild-type and mutated p53 forms were analyzed by essential dynamics to extract the relevant collective motions and by the frustration method to evaluate the degeneracy of the energy landscape. We found that p53 is characterized by wide collective motions and its energy landscape exhibits a rather high frustration level, especially in the regions involved in the binding to physiological ligands. Punctual mutations give rise to a modulation of both the collective motions and the frustration of p53, with different effects depending on the mutation. The regions of p53wt and of the mutated forms characterized by a high frustration level are also largely involved in the collective motions. Such a correlation is discussed also in connection with the intrinsic disordered character of p53 and with its central functional role.

Structure and diffusive dynamics of aspartate α-decarboxylase (ADC) liganded with D-serine in aqueous solution

Incoherent neutron spectroscopy, in combination with dynamic light scattering, was used to investigate the effect of ligand binding on the center-of-mass self-diffusion and internal diffusive dynamics of Escherichia coli aspartate α-decarboxylase (ADC). The X-ray crystal structure of ADC in complex with the d-serine inhibitor was also determined, and molecular dynamics simulations were used to further probe the structural rearrangements that occur as a result of ligand binding. These experiments reveal that d-serine forms hydrogen bonds with some of the active site residues, that higher order oligomers of the ADC tetramer exist on ns–ms time-scales, and also show that ligand binding both affects the ADC internal diffusive dynamics and appears to further increase the size of the higher order oligomers.

Neutron spectroscopy, dynamic light scattering, X-ray diffraction, and MD-simulations were used to investigate the effect of ligand binding on the structure and diffusive dynamics of Escherichia coli aspartate alpha-decarboxylase.

Double proton transfer in hydrated formic acid dimer: Interplay of spatial symmetry and solvent-generated force on reactivity

The double proton transfer (DPT) reaction in the hydrated formic acid dimer (FAD) is investigated at molecular-level detail. For this, a global and reactive machine learned (ML) potential energy surface (PES) is developed to run extensive (more than 100 ns) mixed ML/MM molecular dynamics (MD) simulations in explicit molecular mechanics (MM) solvent at MP2-quality for the solute. Simulations with fixed – as in a conventional empirical force field – and conformationally fluctuating – as available from the ML-based PES – charge models for FAD show a significant impact on the competition between DPT and dissociation of FAD into two formic acid monomers. With increasing temperature the barrier height for DPT in solution changes by about 10% (∼1 kcal mol⁻¹) between 300 K and 600 K. The rate for DPT is largest, ∼1 ns⁻¹, at 350 K and decreases for higher temperatures due to destabilisation and increased probability for dissociation of FAD. The water solvent is found to promote the first proton transfer by exerting a favourable solvent-induced Coulomb force along the O–H⋯O hydrogen bond whereas the second proton transfer is significantly controlled by the O–O separation and other conformational degrees of freedom. Double proton transfer in hydrated FAD is found to involve a subtle interplay and balance between structural and electrostatic factors.

Simulation of double proton transfer in formic acid dimer by reactive ML potential in explicit molecular mechanics water solvent.

Effect of dopamine-functionalization, charge and pH on protein corona formation around TiO2 nanoparticles

Inorganic nanoparticles (NPs) are gaining increasing attention in nanomedicine because of their stimuli responsiveness, which allows combining therapy with diagnosis. However, little information is known about their interaction with intracellular or plasma proteins when they are introduced in a biological environment. Here we present atomistic molecular dynamics (MD) simulations investigating the case study of dopamine-functionalized TiO₂ nanoparticles and two proteins that are overexpressed in cancer cells, i.e. PARP1 and HSP90, since experiments proved them to be the main components of the corona in cell cultures. The mechanism and the nature of the interaction (electrostatic, van der Waals, H-bonds, etc.) is unravelled by defining the protein residues that are more frequently in contact with the NPs, the extent of contact surface area and the variations in the protein secondary structures, at different pH and ionic strength conditions of the solution where they are immersed to simulate a realistic biological environment. The effects of the NP surface functionalization and charge are also considered. Our MD results suggest that less acidic intracellular pH conditions in the presence of cytosolic ionic strength enhance PARP1 interaction with the nanoparticle, whereas the HSP90 contribution is partly weakened, providing a rational explanation to existing experimental observations.

Transfer learned potential energy surfaces: accurate anharmonic vibrational dynamics and dissociation energies for the formic acid monomer and dimer

The vibrational dynamics of the formic acid monomer (FAM) and dimer (FAD) is investigated from machine-learned potential energy surfaces at the MP2 (PES_MP2) and transfer-learned (PES_TL) to the CCSD(T) levels of theory. The normal mode (MAEs of 17.6 and 25.1 cm⁻¹) and second order vibrational perturbation theory (VPT2, MAEs of 6.7 and 17.1 cm⁻¹) frequencies from PES_TL for all modes below 2000 cm⁻¹ for FAM and FAD agree favourably with experiment. For the OH stretch mode the experimental frequencies are overestimated by more than 150 cm⁻¹ for both FAM and FAD from normal mode calculations. Conversely, VPT2 calculations on PES_TL for FAM reproduce the experimental OH frequency to within 22 cm⁻¹. For FAD the VPT2 calculations find the high-frequency OH stretch at 3011 cm⁻¹, compared with an experimentally reported, broad (∼100 cm⁻¹) absorption band with center frequency estimated at ∼3050 cm⁻¹. In agreement with earlier reports, MD simulations at higher temperature shift the position of the OH-stretch in FAM to the red, consistent with improved sampling of the anharmonic regions of the PES. However, for FAD the OH-stretch shifts to the blue and for temperatures higher than 1000 K the dimer partly or fully dissociates using PES_TL. Including zero-point energy corrections from diffusion Monte Carlo simulations for FAM and FAD and corrections due to basis set superposition and completeness errors yields a dissociation energy of D₀ = −14.23 ± 0.08 kcal mol⁻¹ compared with an experimentally determined value of −14.22 ± 0.12 kcal mol⁻¹.

Neural network based PESs are constructed for formic acid monomer and dimer at the MP2 and transfer learned to the CCSD(T) level of theory. The PESs are used to study the vibrational dynamics and dissociation energy of the molecules.

Femtosecond X-ray spectroscopy of haem proteins

We discuss our recently reported femtosecond (fs) X-ray emission spectroscopy results on the ligand dissociation and recombination in nitrosylmyoglobin (MbNO) in the context of previous studies on ferrous haem proteins. We also present a preliminary account of femtosecond X-ray absorption studies on MbNO, pointing to the presence of more than one species formed upon photolysis.

Spectroscopy, Dynamics, and Hydration of S-Nitrosylated Myoglobin

S-Nitrosylation, the covalent addition of NO to the thiol side chain of cysteine, is an important post-transitional modification that can alter the function of various proteins. The structural dynamics and vibrational spectroscopy of S-nitrosylation in the condensed phase are investigated for the methyl-capped cysteine model system and for myoglobin. Using conventional point charge and physically more realistic multipolar force fields for the -SNO group, it is found that the SN- and NO-stretch and the SNO-bend vibrations can be located and distinguished from the other protein modes for simulations of MbSNO at 50 K. The finding of stable cis- and trans-MbSNO agrees with experimental findings on other proteins as is the observation of buried -SNO. For MbSNO the observed relocation of the EF loop in the simulations by ∼3 Å is consistent with the available X-ray structure, and the conformations adopted by the -SNO label are in good overall agreement with the X-ray structure. Despite the larger size of the -SNO group compared with -SH, MbSNO recruits more water molecules in the first two hydration shells due to stronger electrostatic interactions. Similarly, when comparing the hydration between the A- and H-helices, they differ by up to 30% between WT and MbSNO. This suggests that local hydration can also be significantly modulated through nitrosylation.

Dopamine-Decorated TiO2 Nanoparticles in Water: A QM/MM vs an MM Description

Nanoparticle functionalization is a modern strategy in nanotechnology to build up devices for several applications. Modeling fully decorated metal oxide nanoparticles of realistic size (few nanometers) in an aqueous environment is a challenging task. In this work, we present a case study relevant for solar-light exploitation and for biomedical applications, i.e., a dopamine-functionalized TiO₂ nanoparticle (1700 atoms) in bulk water, for which we have performed an extensive comparative investigation with both MM and QM/MM approaches of the structural properties and of the conformational dynamics. We have used a combined multiscale protocol for a more efficient exploration of the complex conformational space. On the basis of the results of this study and of some QM and experimental data, we have defined strengths and limitations of the existing force field parameters. Our findings will be useful for an improved modeling and simulation of many other similar hybrid bioinorganic nanosystems in an aqueous environment that are pivotal in a broad range of nanotechnological applications.

Deciphering the Curly Arrow Representation and Electron Flow for the 1,3-Dipolar Rearrangement between Acetonitrile Oxide and (1S,2R,4S)-2-Cyano-7-oxabicyclo[2.2.1]hept-5-en-2-yl Acetate Derivatives

This study is focused on describing the molecular mechanism beyond the molecular picture provided by the evolution of molecular orbitals, valence bond structures along the reaction progress, or conceptual density functional theory. Using bonding evolution theory (BET) analysis, we have deciphered the mechanism of the 1,3-dipolar rearrangement between acetonitrile oxide and (1S,2R,4S)-2-cyano-7-oxabicyclo[2.2.1]hept-5-en-2-yl acetate derivatives. The BET study revealed that the formation of the C-C bond takes place via a usual sharing model before the O-C one that is also formed in the halogenated species through a not very usual sharing model. The mechanism includes depopulation of the electron density at the N-C triple bond and creation of the V(N) and V(C) monosynaptic basins, depopulation of the former C-C double bond with the creation of V(C,C) basins, and final formation of the V(O,C) basin associated with the O-C bond. The topological changes along the reaction pathway take place in a highly synchronous way. BET provides a convenient quantitative method for deriving curly arrows and electron flow representation to unravel molecular mechanisms.

Single-molecule FRET methods to study the dynamics of proteins at work

Feynman commented that "Everything that living things do can be understood in terms of the jiggling and wiggling of atoms". Proteins can jiggle and wiggle large structural elements such as domains and subunits as part of their functional cycles. Single-molecule fluorescence resonance energy transfer (smFRET) is an excellent tool to study conformational dynamics and decipher coordinated large-scale motions within proteins. smFRET methods introduced in recent years are geared toward understanding the time scales and amplitudes of function-related motions. This review discusses the methodology for obtaining and analyzing smFRET temporal trajectories that provide direct dynamic information on transitions between conformational states. It also introduces correlation methods that are useful for characterizing intramolecular motions. This arsenal of techniques has been used to study multiple molecular systems, from membrane proteins through molecular chaperones, and we examine some of these studies here. Recent exciting methodological novelties permit revealing very fast, submillisecond dynamics, whose relevance to protein function is yet to be fully grasped.

Probing the Differential Dynamics of the Monomeric and Dimeric Insulin from Amide-I IR Spectroscopy

The monomer-dimer equilibrium for insulin is one of the essential steps in forming the receptor-binding competent monomeric form of the hormone. Despite this importance, the thermodynamic stability, in particular for modified insulins, is quite poorly understood, in part, due to experimental difficulties. This work explores one- and two-dimensional infrared spectroscopy in the range of the amide-I band for the hydrated monomeric and dimeric wild-type hormone. It is found that for the monomer the frequency fluctuation correlation function (FFCF) and the one-dimensional infrared spectra are position sensitive. The spectra of the -CO probes at the dimerization interface (residues Phe24, Phe25, and Tyr26) split and indicate an asymmetry despite the overall (formal) point symmetry of the dimer structure. Also, the decay times of the FFCF for the same -CO probe in the monomer and the dimer can differ by up to 1 order of magnitude, for example, for residue PheB24, which is solvent exposed for the monomer but at the interface for the dimer. The spectroscopic shifts correlate approximately with the average number of hydration waters and the magnitude of the FFCF at time zero. However, this correlation is only qualitative due to the heterogeneous and highly dynamical environment. Based on density functional theory calculations, the dominant contribution for solvent-exposed -CO is found to arise from the surrounding water (∼75%), whereas the protein environment contributes considerably less. The results suggest that infrared spectroscopy is a positionally sensitive probe of insulin dimerization, in particular in conjunction with isotopic labeling of the probe.

Vibrational Spectroscopy of N3- in the Gas and Condensed Phase

Azido-derivatized amino acids are potentially useful, positionally resolved spectroscopic probes for studying the structural dynamics of proteins and macromolecules in solution. To this end, a computational model for the vibrational modes of N₃^- based on accurate electronic structure calculations and a reproducing kernel Hilbert space representation of the potential energy surface for the internal degrees of freedom is developed. Fully dimensional quantum bound state calculations yield the antisymmetric stretch vibration at 1974 cm^-1 compared with 1986 cm^-1 from experiment. This mode shifts by 64 cm^-1 (from the frequency distribution) and 74 cm^-1 (from the IR line shape) to the blue, respectively, compared with 61 cm^-1 from experiment for N₃^- in water. The decay time of the frequency fluctuation correlation function is 1.1 ps, which is in good agreement with experiment (1.2-1.3 ps) and the full width at half maximum of the asymmetric stretch in solution is 18.5 cm^-1 compared with 25.2 cm^-1 from experiment. A computationally more efficient analysis based on instantaneous normal modes is shown to provide comparable, albeit somewhat less quantitative results compared to solving the three-dimensional Schrödinger equation for the fundamental vibrations.

Vibrational Energy in Proteins Correlates with Topology

The exchange of vibrational energy in proteins is crucial for their function. Here, we establish a connection between quantities related to it with geometry-based properties such as the proteins' residues coordination number. This relation is proven by molecular simulation in a neuro-pharmacologically relevant transmembrane receptor. The connection demonstrated here paves the way to studies of protein allostery and conformational changes based solely on protein structure.

Perspective: Opportunities for ultrafast science at SwissFEL

We present the main specifications of the newly constructed Swiss Free Electron Laser, SwissFEL, and explore its potential impact on ultrafast science. In light of recent achievements at current X-ray free electron lasers, we discuss the potential territory for new scientific breakthroughs offered by SwissFEL in Chemistry, Biology, and Materials Science, as well as nonlinear X-ray science.