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.

Observation of a single protein by ultrafast X-ray diffraction.

The idea of using ultrashort X-ray pulses to obtain images of single proteins frozen in time has fascinated and inspired many. It was one of the arguments for building X-ray free-electron lasers. According to theory, the extremely intense pulses provide sufficient signal to dispense with using crystals as an amplifier, and the ultrashort pulse duration permits capturing the diffraction data before the sample inevitably explodes. This was first demonstrated on biological samples a decade ago on the giant mimivirus. Since then, a large collaboration has been pushing the limit of the smallest sample that can be imaged. The ability to capture snapshots on the timescale of atomic vibrations, while keeping the sample at room temperature, may allow probing the entire conformational phase space of macromolecules. Here we show the first observation of an X-ray diffraction pattern from a single protein, that of Escherichia coli GroEL which at 14 nm in diameter is the smallest biological sample ever imaged by X-rays, and demonstrate that the concept of diffraction before destruction extends to single proteins. From the pattern, it is possible to determine the approximate orientation of the protein. Our experiment demonstrates the feasibility of ultrafast imaging of single proteins, opening the way to single-molecule time-resolved studies on the femtosecond timescale.

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.

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.

Decreasing ultrafast x-ray pulse durations with saturable absorption and resonant transitions.

Saturable absorption is a nonlinear effect where a material's ability to absorb light is frustrated due to a high influx of photons and the creation of electron vacancies. Experimentally induced saturable absorption in copper revealed a reduction in the temporal duration of transmitted x-ray laser pulses, but a detailed account of changes in opacity and emergence of resonances is still missing. In this computational work, we employ nonlocal thermodynamic equilibrium plasma simulations to study the interaction of femtosecond x rays and copper. Following the onset of frustrated absorption, we find that a K-M resonant transition occurring at highly charged states turns copper opaque again. The changes in absorption generate a transient transparent window responsible for the shortened transmission signal. We also propose using fluorescence induced by the incident beam as an alternative source to achieve shorter x-ray pulses. Intense femtosecond x rays are valuable to probe the structure and dynamics of biological samples or to reach extreme states of matter. Shortened pulses could be relevant for emerging imaging techniques.

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 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.

Distinguishing between Similar Miniproteins with Single-Molecule Nanopore Sensing: A Computational Study.

A nanopore is a tool in single-molecule sensing biotechnology that offers label-free identification with high throughput. Nanopores have been successfully applied to sequence DNA and show potential in the study of proteins. Nevertheless, the task remains challenging due to the large variability in size, charges, and folds of proteins. Miniproteins have a small number of residues, limited secondary structure, and stable tertiary structure, which can offer a systematic way to reduce complexity. In this computational work, we theoretically evaluated sensing two miniproteins found in the human body using a silicon nitride nanopore. We employed molecular dynamics methods to compute occupied-pore ionic current magnitudes and electronic structure calculations to obtain interaction strengths between pore wall and miniprotein. From the interaction strength, we derived dwell times using a mix of combinatorics and numerical solutions. This latter approach circumvents typical computational demands needed to simulate translocation events using molecular dynamics. We focused on two miniproteins potentially difficult to distinguish owing to their isotropic geometry, similar number of residues, and overall comparable structure. We found that the occupied-pore current magnitudes not to vary significantly, but their dwell times differ by 1 order of magnitude. Together, these results suggest a successful identification protocol for similar miniproteins.

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.

Structural Heterogeneity in Single Particle Imaging Using X-ray Lasers.

One of the challenges facing single particle imaging with ultrafast X-ray pulses is the structural heterogeneity of the sample to be imaged. For the method to succeed with weakly scattering samples, the diffracted images from a large number of individual proteins need to be averaged. The more the individual proteins differ in structure, the lower the achievable resolution in the final reconstructed image. We use molecular dynamics to simulate two globular proteins in vacuum, fully desolvated as well as with two different solvation layers, at various temperatures. We calculate the diffraction patterns based on the simulations and evaluate the noise in the averaged patterns arising from the structural differences and the surrounding water. Our simulations show that the presence of a minimal water coverage with an average 3 Å thickness will stabilize the protein, reducing the noise associated with structural heterogeneity, whereas additional water will generate more background noise.

The role of transient resonances for ultra-fast imaging of single sucrose nanoclusters.

Intense x-ray free-electron laser (XFEL) pulses hold great promise for imaging function in nanoscale and biological systems with atomic resolution. So far, however, the spatial resolution obtained from single shot experiments lags averaging static experiments. Here we report on a combined computational and experimental study about ultrafast diffractive imaging of sucrose clusters which are benchmark organic samples. Our theoretical model matches the experimental data from the water window to the keV x-ray regime. The large-scale dynamic scattering calculations reveal that transient phenomena driven by non-linear x-ray interaction are decisive for ultrafast imaging applications. Our study illuminates the complex interplay of the imaging process with the rapidly changing transient electronic structures in XFEL experiments and shows how computational models allow optimization of the parameters for ultrafast imaging experiments.

Femtosecond bond breaking and charge dynamics in ultracharged amino acids.

Historically, structure determination of nanocrystals, proteins, and macromolecules required the growth of high-quality crystals sufficiently large to diffract X-rays efficiently while withstanding radiation damage. The development of the X-ray free-electron laser has opened the path toward high resolution single particle imaging, and the extreme intensity of the X-rays ensures that enough diffraction statistics are collected before the sample is destroyed by radiation damage. Still, recovery of the structure is a challenge, in part due to the partial fragmentation of the sample during the diffraction event. In this study, we use first-principles based methods to study the impact of radiation induced ionization of six amino acids on the reconstruction process. In particular, we study the fragmentation and charge rearrangement to elucidate the time scales involved and the characteristic fragments occurring.

Electrospray sample injection for single-particle imaging with x-ray lasers.

The possibility of imaging single proteins constitutes an exciting challenge for x-ray lasers. Despite encouraging results on large particles, imaging small particles has proven to be difficult for two reasons: not quite high enough pulse intensity from currently available x-ray lasers and, as we demonstrate here, contamination of the aerosolized molecules by nonvolatile contaminants in the solution. The amount of contamination on the sample depends on the initial droplet size during aerosolization. Here, we show that, with our electrospray injector, we can decrease the size of aerosol droplets and demonstrate virtually contaminant-free sample delivery of organelles, small virions, and proteins. The results presented here, together with the increased performance of next-generation x-ray lasers, constitute an important stepping stone toward the ultimate goal of protein structure determination from imaging at room temperature and high temporal resolution.

Erratum: Experimental strategies for imaging bioparticles with femtosecond hard X-ray pulses. Corrigendum.

[This corrects the article DOI: 10.1107/S2052252517003591.].

Considerations for three-dimensional image reconstruction from experimental data in coherent diffractive imaging.

Diffraction before destruction using X-ray free-electron lasers (XFELs) has the potential to determine radiation-damage-free structures without the need for crystallization. This article presents the three-dimensional reconstruction of the Melbournevirus from single-particle X-ray diffraction patterns collected at the LINAC Coherent Light Source (LCLS) as well as reconstructions from simulated data exploring the consequences of different kinds of experimental sources of noise. The reconstruction from experimental data suffers from a strong artifact in the center of the particle. This could be reproduced with simulated data by adding experimental background to the diffraction patterns. In those simulations, the relative density of the artifact increases linearly with background strength. This suggests that the artifact originates from the Fourier transform of the relatively flat background, concentrating all power in a central feature of limited extent. We support these findings by significantly reducing the artifact through background removal before the phase-retrieval step. Large amounts of blurring in the diffraction patterns were also found to introduce diffuse artifacts, which could easily be mistaken as biologically relevant features. Other sources of noise such as sample heterogeneity and variation of pulse energy did not significantly degrade the quality of the reconstructions. Larger data volumes, made possible by the recent inauguration of high repetition-rate XFELs, allow for increased signal-to-background ratio and provide a way to minimize these artifacts. The anticipated development of three-dimensional Fourier-volume-assembly algorithms which are background aware is an alternative and complementary solution, which maximizes the use of data.

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.

Reproducibility of single protein explosions induced by X-ray lasers.

Single particle imaging (SPI) using X-ray pulses has become increasingly attainable with the advent of high-intensity free electron lasers. Eliminating the need for crystallized samples enables structural studies of molecules previously inaccessible by conventional crystallography. While this emerging technique already demonstrates substantial promise, some obstacles need to be overcome before SPI can reach its full potential. One such problem is determining the spatial orientation of the sample at the time of X-ray interaction. Existing solutions rely on diffraction data and are computationally demanding and sensitive to noise. In this in silico study, we explore the possibility of aiding these methods by mapping the ion distribution as the sample undergoes a Coulomb explosion following the intense ionization. By detecting the ions ejected from the fragmented sample, the orientation of the original sample should be possible to determine. Knowledge of the orientation has been shown earlier to be of substantial advantage in the reconstruction of the original structure. 150 explosions of each of twelve separate systems - four polypeptides with different amounts of surface bound water - were simulated with molecular dynamics (MD) and the average angular distribution of carbon and sulfur ions was investigated independently. The results show that the explosion maps are reproducible in both cases, supporting the idea that orientation information is preserved. Additional water seems to restrict the carbon ion trajectories further through a shielding mechanism, making the maps more distinct. For sulfurs, water has no significant impact on the trajectories, likely due to their higher mass and greater ionization cross section, indicating that they could be of particular interest. Based on these findings, we conclude that explosion data can aid spatial orientation in SPI experiments and could substantially improve the capabilities of the novel technique.

Hit detection in serial femtosecond crystallography using X-ray spectroscopy of plasma emission.

Serial femtosecond crystallography is an emerging and promising method for determining protein structures, making use of the ultrafast and bright X-ray pulses from X-ray free-electron lasers. The upcoming X-ray laser sources will produce well above 1000 pulses per second and will pose a new challenge: how to quickly determine successful crystal hits and avoid a high-rate data deluge. Proposed here is a hit-finding scheme based on detecting photons from plasma emission after the sample has been intercepted by the X-ray laser. Plasma emission spectra are simulated for systems exposed to high-intensity femtosecond pulses, for both protein crystals and the liquid carrier systems that are used for sample delivery. The thermal radiation from the glowing plasma gives a strong background in the XUV region that depends on the intensity of the pulse, around the emission lines from light elements (carbon, nitrogen, oxygen). Sample hits can be reliably distinguished from the carrier liquid based on the characteristic emission lines from heavier elements present only in the sample, such as sulfur. For buffer systems with sulfur present, selenomethionine substitution is suggested, where the selenium emission lines could be used both as an indication of a hit and as an aid in phasing and structural reconstruction of the protein.

Experimental strategies for imaging bioparticles with femtosecond hard X-ray pulses.

This study explores the capabilities of the Coherent X-ray Imaging Instrument at the Linac Coherent Light Source to image small biological samples. The weak signal from small samples puts a significant demand on the experiment. Aerosolized Omono River virus particles of ∼40 nm in diameter were injected into the submicrometre X-ray focus at a reduced pressure. Diffraction patterns were recorded on two area detectors. The statistical nature of the measurements from many individual particles provided information about the intensity profile of the X-ray beam, phase variations in the wavefront and the size distribution of the injected particles. The results point to a wider than expected size distribution (from ∼35 to ∼300 nm in diameter). This is likely to be owing to nonvolatile contaminants from larger droplets during aerosolization and droplet evaporation. The results suggest that the concentration of nonvolatile contaminants and the ratio between the volumes of the initial droplet and the sample particles is critical in such studies. The maximum beam intensity in the focus was found to be 1.9 × 1012 photons per µm2 per pulse. The full-width of the focus at half-maximum was estimated to be 500 nm (assuming 20% beamline transmission), and this width is larger than expected. Under these conditions, the diffraction signal from a sample-sized particle remained above the average background to a resolution of 4.25 nm. The results suggest that reducing the size of the initial droplets during aerosolization is necessary to bring small particles into the scope of detailed structural studies with X-ray lasers.