The Surface of Ionic Liquids in Water: From an Ionic Tug of War to a Quasi-Ordered Two-Dimensional Layer.

The sustainable development encompasses the search for new materials for energy storage, gas capture, separation, and solvents in industrial processes that can substitute conventional ones in an efficient and clean manner. Ionic liquids (ILs) emerged and have been advanced as alternative materials for such applications, but an obstacle is their hygroscopicity and the effects on their physical properties in the presence of humidity. Several industrial processes depend on the aqueous interfacial properties, and the main focus of this work is the water/IL interface. The behavior of the aqueous ionic liquids at the water-vacuum interface is representative for their water interfacial properties. Using X-ray photoelectron spectroscopy in combination with molecular dynamics simulations we investigate four aqueous IL systems, and provide molecular level insight on the interfacial behaviour of the ionic liquids, such as ion-pair formation, orientation and surface concentration. We find that ionic liquids containing a chloride anion have a lowered surface enrichment due to the low surface propensity of chloride. In contrast, the ionic liquids containing a bistriflimide anion are extremely surface-enriched due to cooperative surface propensity between the cations and anions, forming a two-dimensional ionic liquid on the water surface at low concentrations.

Collision induced unfolding and molecular dynamics simulations of norovirus capsid dimers reveal strain-specific stability profiles.

Collision induced unfolding (CIU) is a method used with ion mobility mass spectrometry to examine protein structures and their stability. Such experiments yield information about higher order protein structures, yet are unable to provide details about the underlying processes. That information can however be provided using molecular dynamics simulations. Here, we investigate the gas-phase unfolding of norovirus capsid dimers from the Norwalk and Kawasaki strains by employing molecular dynamics simulations over a range of temperatures, representing different levels of activation, together with CIU experiments. The dimers have highly similar structures, but their CIU reveals different stability that can be explained by the different dynamics that arises in response to the activation seen in the simulations, including a part of the sequence with previously observed strain-specific dynamics in solution. Our findings show how similar protein variants can be examined using mass spectrometric techniques in conjunction with atomistic molecular dynamics simulations to reveal differences in stability as well as differences in how and where unfolding takes place upon activation.

Heavy element incorporation in nitroimidazole radiosensitizers: molecular-level insights into fragmentation dynamics.

The present study investigates the photofragmentation behavior of iodine-enhanced nitroimidazole-based radiosensitizer model compounds in their protonated form using near-edge X-ray absorption mass spectrometry and quantum mechanical calculations. These molecules possess dual functionality: improved photoabsorption capabilities and the ability to generate species that are relevant to cancer sensitization upon photofragmentation. Four samples were investigated by scanning the generated fragments in the energy regions around C 1s, N 1s, O 1s, and I 3d-edges with a particular focus on NO2+ production. The experimental summed ion yield spectra are explained using the theoretical near-edge X-ray absorption fine structure spectrum based on density functional theory. Born-Oppenheimer-based molecular dynamics simulations were performed to investigate the fragmentation processes.

Shell-dependent photofragmentation dynamics of a heavy-atom-containing bifunctional nitroimidazole radiosensitizer.

Radiation therapy uses ionizing radiation to break chemical bonds in cancer cells, thereby causing DNA damage and leading to cell death. The therapeutic effectiveness can be further increased by making the tumor cells more sensitive to radiation. Here, we investigate the role of the initial halogen atom core hole on the photofragmentation dynamics of 2-bromo-5-iodo-4-nitroimidazole, a potential bifunctional radiosensitizer. Bromine and iodine atoms were included in the molecule to increase the photoionization cross-section of the radiosensitizer at higher photon energies. The fragmentation dynamics of the molecule was studied experimentally in the gas phase using photoelectron-photoion-photoion coincidence spectroscopy and computationally using Born-Oppenheimer molecular dynamics. We observed significant changes between shallow core (I 4d, Br 3d) and deep core (I 3d) ionization in fragment formation and their kinetic energies. Despite the fact, that the ions ejected after deep core ionization have higher kinetic energies, we show that in a cellular environment, the ion spread is not much larger, keeping the damage well-localized.

Enhanced EMC-Advantages of partially known orientations in x-ray single particle imaging.

Single particle imaging of proteins in the gas phase with x-ray free-electron lasers holds great potential to study fast protein dynamics, but is currently limited by weak and noisy data. A further challenge is to discover the proteins' orientation as each protein is randomly oriented when exposed to x-rays. Algorithms such as the expand, maximize, and compress (EMC) exist that can solve the orientation problem and reconstruct the three-dimensional diffraction intensity space, given sufficient measurements. If information about orientation were known, for example, by using an electric field to orient the particles, the reconstruction would benefit and potentially reach better results. We used simulated diffraction experiments to test how the reconstructions from EMC improve with particles' orientation to a preferred axis. Our reconstructions converged to correct maps of the three-dimensional diffraction space with fewer measurements if biased orientation information was considered. Even for a moderate bias, there was still significant improvement. Biased orientations also substantially improved the results in the case of missing central information, in particular in the case of small datasets. The effects were even more significant when adding a background with 50% the strength of the averaged diffraction signal photons to the diffraction patterns, sometimes reducing the data requirement for convergence by a factor of 10. This demonstrates the usefulness of having biased orientation information in single particle imaging experiments, even for a weaker bias than what was previously known. This could be a key component in overcoming the problems with background noise that currently plague these experiments.

MolDStruct: Modeling the dynamics and structure of matter exposed to ultrafast x-ray lasers with hybrid collisional-radiative/molecular dynamics.

We describe a method to compute photon-matter interaction and atomic dynamics with x-ray lasers using a hybrid code based on classical molecular dynamics and collisional-radiative calculations. The forces between the atoms are dynamically determined based on changes to their electronic occupations and the formation of a free electron cloud created from the irradiation of photons in the x-ray spectrum. The rapid transition from neutral solid matter to dense plasma phase allows the use of screened potentials, reducing the number of non-bonded interactions. In combination with parallelization through domain decomposition, the hybrid code handles large-scale molecular dynamics and ionization. This method is applicable for large enough samples (solids, liquids, proteins, viruses, atomic clusters, and crystals) that, when exposed to an x-ray laser pulse, turn into a plasma in the first few femtoseconds of the interaction. We present four examples demonstrating the applicability of the method. We investigate the non-thermal heating and scattering of bulk water and damage-induced dynamics of a protein crystal using an x-ray pump-probe scheme. In both cases, we compare to the experimental data. For single particle imaging, we simulate the ultrafast dynamics of a methane cluster exposed to a femtosecond x-ray laser. In the context of coherent diffractive imaging, we study the fragmentation as given by an x-ray pump-probe setup to understand the evolution of radiation damage in the time range of hundreds of femtoseconds.

Speckle contrast of interfering fluorescence X-rays.

With the development of X-ray free-electron lasers (XFELs), producing pulses of femtosecond durations comparable with the coherence times of X-ray fluorescence, it has become possible to observe intensity-intensity correlations due to the interference of emission from independent atoms. This has been used to compare durations of X-ray pulses and to measure the size of a focusedX-ray beam, for example. Here it is shown that it is also possible to observe the interference of fluorescence photons through the measurement of the speckle contrast of angle-resolved fluorescence patterns. Speckle contrast is often used as a measure of the degree of coherence of the incident beam or the fluctuations of the illuminated sample as determined from X-ray diffraction patterns formed by elastic scattering, rather than from fluorescence patterns as addressed here. Commonly used approaches to estimate speckle contrast were found to suffer when applied to XFEL-generated fluorescence patterns due to low photon counts and a significant variation of the excitation pulse energy from shot to shot. A new method to reliably estimate speckle contrast under such conditions, using a weighting scheme, is introduced. The method is demonstrated by comparing the speckle contrast of fluorescence observed with pulses of 3 fs to 15 fs duration.

Imaging via Correlation of X-Ray Fluorescence Photons.

We demonstrate that x-ray fluorescence emission, which cannot maintain a stationary interference pattern, can be used to obtain images of structures by recording photon-photon correlations in the manner of the stellar intensity interferometry of Hanbury Brown and Twiss. This is achieved utilizing femtosecond-duration pulses of a hard x-ray free-electron laser to generate the emission in exposures comparable to the coherence time of the fluorescence. Iterative phasing of the photon correlation map generated a model-free real-space image of the structure of the emitters. Since fluorescence can dominate coherent scattering, this may enable imaging uncrystallised macromolecules.

Coherent diffractive imaging of proteins and viral capsids: simulating MS SPIDOC.

MS SPIDOC is a novel sample delivery system designed for single (isolated) particle imaging at X-ray Free-Electron Lasers that is adaptable towards most large-scale facility beamlines. Biological samples can range from small proteins to MDa particles. Following nano-electrospray ionization, ionic samples can be m/z-filtered and structurally separated before being oriented at the interaction zone. Here, we present the simulation package developed alongside this prototype. The first part describes how the front-to-end ion trajectory simulations have been conducted. Highlighted is a quadrant lens; a simple but efficient device that steers the ion beam within the vicinity of the strong DC orientation field in the interaction zone to ensure spatial overlap with the X-rays. The second part focuses on protein orientation and discusses its potential with respect to diffractive imaging methods. Last, coherent diffractive imaging of prototypical T = 1 and T = 3 norovirus capsids is shown. We use realistic experimental parameters from the SPB/SFX instrument at the European XFEL to demonstrate that low-resolution diffractive imaging data (q < 0.3 nm-1) can be collected with only a few X-ray pulses. Such low-resolution data are sufficient to distinguish between both symmetries of the capsids, allowing to probe low abundant species in a beam if MS SPIDOC is used as sample delivery.

Observations of phase changes in monoolein during high viscous injection.

Serial crystallography of membrane proteins often employs high-viscosity injectors (HVIs) to deliver micrometre-sized crystals to the X-ray beam. Typically, the carrier medium is a lipidic cubic phase (LCP) media, which can also be used to nucleate and grow the crystals. However, despite the fact that the LCP is widely used with HVIs, the potential impact of the injection process on the LCP structure has not been reported and hence is not yet well understood. The self-assembled structure of the LCP can be affected by pressure, dehydration and temperature changes, all of which occur during continuous flow injection. These changes to the LCP structure may in turn impact the results of X-ray diffraction measurements from membrane protein crystals. To investigate the influence of HVIs on the structure of the LCP we conducted a study of the phase changes in monoolein/water and monoolein/buffer mixtures during continuous flow injection, at both atmospheric pressure and under vacuum. The reservoir pressure in the HVI was tracked to determine if there is any correlation with the phase behaviour of the LCP. The results indicated that, even though the reservoir pressure underwent (at times) significant variation, this did not appear to correlate with observed phase changes in the sample stream or correspond to shifts in the LCP lattice parameter. During vacuum injection, there was a three-way coexistence of the gyroid cubic phase, diamond cubic phase and lamellar phase. During injection at atmospheric pressure, the coexistence of a cubic phase and lamellar phase in the monoolein/water mixtures was also observed. The degree to which the lamellar phase is formed was found to be strongly dependent on the co-flowing gas conditions used to stabilize the LCP stream. A combination of laboratory-based optical polarization microscopy and simulation studies was used to investigate these observations.

Imaging of femtosecond bond breaking and charge dynamics in ultracharged peptides.

X-ray free-electrons lasers have revolutionized the method of imaging biological macromolecules such as proteins, viruses and cells by opening the door to structural determination of both single particles and crystals at room temperature. By utilizing high intensity X-ray pulses on femtosecond timescales, the effects of radiation damage can be reduced. Achieving high resolution structures will likely require knowledge of how radiation damage affects the structure on an atomic scale, since the experimentally obtained electron densities will be reconstructed in the presence of radiation damage. Detailed understanding of the expected damage scenarios provides further information, in addition to guiding possible corrections that may need to be made to obtain a damage free reconstruction. In this work, we have quantified the effects of ionizing photon-matter interactions using first principles molecular dynamics. We utilize density functional theory to calculate bond breaking and charge dynamics in three ultracharged molecules and two different structural conformations that are important to the structural integrity of biological macromolecules, comparing to our previous studies on amino acids. The effects of the ultracharged states and subsequent bond breaking in real space are studied in reciprocal space using coherent diffractive imaging of an ensemble of aligned biomolecules in the gas phase.

Probing aqueous ions with non-local Auger relaxation.

Non-local analogues of Auger decay are increasingly recognized as important relaxation processes in the condensed phase. Here, we explore non-local autoionization, specifically Intermolecular Coulombic Decay (ICD), of a series of aqueous-phase isoelectronic cations following 1s core-level ionization. In particular, we focus on Na+, Mg2+, and Al3+ ions. We unambiguously identify the ICD contribution to the K-edge Auger spectrum. The different strength of the ion-water interactions is manifested by varying intensities of the respective signals: the ICD signal intensity is greatest for the Al3+ case, weaker for Mg2+, and absent for weakly-solvent-bound Na+. With the assistance of ab initio calculations and molecular dynamics simulations, we provide a microscopic understanding of the non-local decay processes. We assign the ICD signals to decay processes ending in two-hole states, delocalized between the central ion and neighbouring water. Importantly, these processes are shown to be highly selective with respect to the promoted water solvent ionization channels. Furthermore, using a core-hole-clock analysis, the associated ICD timescales are estimated to be around 76 fs for Mg2+ and 34 fs for Al3+. Building on these results, we argue that Auger and ICD spectroscopy represents a unique tool for the exploration of intra- and inter-molecular structure in the liquid phase, simultaneously providing both structural and electronic information.

X-ray Induced Fragmentation of Protonated Cystine.

We demonstrate site-specific X-ray induced fragmentation across the sulfur L-edge of protonated cystine, the dimer of the amino acid cysteine. Ion yield NEXAFS were performed in the gas phase using electrospray ionization (ESI) in combination with an ion trap. The interpretation of the sulfur L-edge NEXAFS spectrum is supported by Restricted Open-Shell Configuration Interaction (ROCIS) calculations. The fragmentation pathway of triply charged cystine ions was modeled by Molecular Dynamics (MD) simulations. We have deduced a possible pathway of fragmentation upon excitation and ionization of S 2p electrons. The disulfide bridge breaks for resonant excitation at lower photon energies but remains intact upon higher energy resonant excitation and upon ionization of S 2p. The larger fragments initially formed subsequently break into smaller fragments.

Protein orientation in time-dependent electric fields: orientation before destruction.

Proteins often have nonzero electric dipole moments, making them interact with external electric fields and offering a means for controlling their orientation. One application that is known to benefit from orientation control is single-particle imaging with x-ray free-electron lasers, in which diffraction is recorded from proteins in the gas phase to determine their structures. To this point, theoretical investigations into this phenomenon have assumed that the field experienced by the proteins is constant or a perfect step function, whereas any real-world pulse will be smooth. Here, we explore the possibility of orienting gas-phase proteins using time-dependent electric fields. We performed ab initio simulations to estimate the field strength required to break protein bonds, with 45 V/nm as a breaking point value. We then simulated ubiquitin in time-dependent electric fields using classical molecular dynamics. The minimal field strength required for orientation within 10 ns was on the order of 0.5 V/nm. Although high fields can be destructive for the structure, the structures in our simulations were preserved until orientation was achieved regardless of field strength, a principle we denote "orientation before destruction."

Ultrafast dynamics and scattering of protic ionic liquids induced by XFEL pulses.

X-rays are routinely used for structural studies through scattering, and femtosecond X-ray lasers can probe ultrafast dynamics. We aim to capture the femtosecond dynamics of liquid samples using simulations and deconstruct the interplay of ionization and atomic motion within the X-ray laser pulse. This deconstruction is resolution dependent, as ionization influences the low momentum transfers through changes in scattering form factors, while atomic motion has a greater effect at high momentum transfers through loss of coherence. Our methodology uses a combination of classical molecular dynamics and plasma simulation on a protic ionic liquid to quantify the contributions to the scattering signal and how these evolve with time during the X-ray laser pulse. Our method is relevant for studies of organic liquids, biomolecules in solution or any low-Z materials at liquid densities that quickly turn into a plasma while probed with X-rays.

Propensity, free energy contributions and conformation of primary n-alcohols at a water surface.

Atmospheric aerosols contain organic molecules that serve as cloud condensation nucleation sites and affect the climate. Several experimental and simulation studies have been dedicated to investigate their surface propensity, but the mechanisms that drive them to the water surface are still not fully understood. In this molecular dynamics (MD) simulation study, primary alcohols are considered as a model system representing polar organic molecules. We find that the surface affinity of n-alcohols increases linearly with the length of the hydrophobic tail. By decomposing the adsorption free energy into enthalpy and entropy contributions, we find that the transition from bulk to surface is entropically driven, compatible with the fact that the hydrophobic effect of small solutes is of entropic origin. The enthalpy of surface adsorption is nearly invariant among different n-alcohols because the loss of solvent-alcohol interactions is balanced by a gain in solvent-solvent interactions. Structural analysis shows that, at the surface, the linear alcohols prefer an orientation with the hydrophobic tail pointing out from the surface, whereas the hydroxyl group remains buried in the water. This general behaviour is likely transferable to other small molecules with similar structures but other functional groups that are present in the atmosphere. Therefore, the present study is a step forward toward a general description of organic molecules in aerosols.

Charge State Dependence of Amino Acid Propensity at Water Surface: Mechanisms Elucidated by Molecular Dynamics Simulations.

Atmospheric aerosols contain a variety of compounds, among them free amino acids and salt ions. The pH of the aerosol droplets depends on their origin and environment. Consequently, compounds like free amino acids found in the droplets will be at different charge states, since these states to a great extent depend on the surrounding pH condition. In droplets of marine origin, amino acids are believed to drive salt ions to the water surface and a pH-dependent amino acid surface propensity will, therefore, indirectly affect many processes in atmospheric chemistry and physics such as for instance cloud condensation. To understand the surface propensity of glycine, valine, and phenylalanine at acidic, neutral, and basic pH, we used molecular dynamics (MD) simulations to investigate them at three different charge states in water. Their respective surface propensities were obtained by the means of a potential of mean force (PMF) in an umbrella sampling approach. Glycine was found to have no preference for the surface, while both valine and phenylalanine showed high propensities. Among the charge states of the surface-enriched ones, the cation, representing the amino acids at low pH, was found to have the highest affinity. Free energy decomposition revealed that the driving forces depend strongly on the nature of the amino acid and its charge state. In phenylalanine, the main factor was found to be a substantial entropy gain, likely related to the side chain, whereas in valine, hydrogen bonding to the functional groups leads to favorable energies and, in turn, affects the surface propensity. A significant gain in water-water enthalpy was seen for both valine and phenylalanine.

Ultrafast nonthermal heating of water initiated by an X-ray Free-Electron Laser.

The bright ultrafast pulses of X-ray Free-Electron Lasers allow investigation into the structure of matter under extreme conditions. We have used single pulses to ionize and probe water as it undergoes a phase transition from liquid to plasma. We report changes in the structure of liquid water on a femtosecond time scale when irradiated by single 6.86 keV X-ray pulses of more than 106 J/cm2 These observations are supported by simulations based on molecular dynamics and plasma dynamics of a water system that is rapidly ionized and driven out of equilibrium. This exotic ionic and disordered state with the density of a liquid is suggested to be structurally different from a neutral thermally disordered state.

Glycan-Induced Protein Dynamics in Human Norovirus P Dimers Depend on Virus Strain and Deamidation Status.

Noroviruses are the major cause of viral gastroenteritis and re-emerge worldwide every year, with GII.4 currently being the most frequent human genotype. The norovirus capsid protein VP1 is essential for host immune response. The P domain mediates cell attachment via histo blood-group antigens (HBGAs) in a strain-dependent manner but how these glycan-interactions actually relate to cell entry remains unclear. Here, hydrogen/deuterium exchange mass spectrometry (HDX-MS) is used to investigate glycan-induced protein dynamics in P dimers of different strains, which exhibit high structural similarity but different prevalence in humans. While the almost identical strains GII.4 Saga and GII.4 MI001 share glycan-induced dynamics, the dynamics differ in the emerging GII.17 Kawasaki 308 and rare GII.10 Vietnam 026 strain. The structural aspects of glycan binding to fully deamidated GII.4 P dimers have been investigated before. However, considering the high specificity and half-life of N373D under physiological conditions, large fractions of partially deamidated virions with potentially altered dynamics in their P domains are likely to occur. Therefore, we also examined glycan binding to partially deamidated GII.4 Saga and GII.4 MI001 P dimers. Such mixed species exhibit increased exposure to solvent in the P dimer upon glycan binding as opposed to pure wildtype. Furthermore, deamidated P dimers display increased flexibility and a monomeric subpopulation. Our results indicate that glycan binding induces strain-dependent structural dynamics, which are further altered by N373 deamidation, and hence hint at a complex role of deamidation in modulating glycan-mediated cell attachment in GII.4 strains.

Proton dynamics in molecular solvent clusters as an indicator for hydrogen bond network strength in confined geometries.

Hydrogen bonding leads to the formation of strong, extended intermolecular networks in molecular liquids such as water. However, it is less well-known how robust the network is to environments in which surface formation or confinement effects become prominent, such as in clusters or droplets. Such systems provide a useful way to probe the robustness of the network, since the degree of confinement can be tuned by altering the cluster size, changing both the surface-to-volume ratio and the radius of curvature. To explore the formation of hydrogen bond networks in confined geometries, here we present O 1s Auger spectra of small and large clusters of water, methanol, and dimethyl ether, as well as their deuterated equivalents. The Auger spectra of the clusters and the corresponding macroscopic liquids are compared and evaluated for an isotope effect, which is due to proton dynamics within the lifetime of the core hole (proton-transfer-mediated charge-separation, PTM-CS), and can be linked to the formation of a hydrogen bond network in the system. An isotope effect is observed in water and methanol but not for dimethyl ether, which cannot donate a hydrogen bond at its oxygen site. The isotope effect, and therefore the strength of the hydrogen bond network, is more pronounced in water than in methanol. Its value depends on the average size of the cluster, indicating that confinement effects change proton dynamics in the core ionised excited state.