Start-to-end simulation of single-particle imaging using ultra-short pulses at the European X-ray Free-Electron Laser

Computational study of diffraction image formation from XFEL irradiated single ribosome molecule

Single particle imaging at atomic resolution is perhaps one of the most desired goals for ultrafast X-ray science with X-ray free-electron lasers. Such a capability would create great opportunity within the biological sciences, as high-resolution structural information of biosamples that may not crystallize is essential for many research areas therein. In this paper, we report on a comprehensive computational study of diffraction image formation during single particle imaging of a macromolecule, containing over one hundred thousand non-hydrogen atoms. For this study, we use a dedicated simulation framework, SIMEX, available at the European XFEL facility. Our results demonstrate the full feasibility of computational single-particle imaging studies for biological samples of realistic size. This finding is important as it shows that the SIMEX platform can be used for simulations to inform relevant single-particle-imaging experiments and help to establish optimal parameters for these experiments. This will enable more focused and more efficient single-particle-imaging experiments at XFEL facilities, making the best use of the resource-intensive XFEL operation.

Characterization of Biological Samples Using Ultra-Short and Ultra-Bright XFEL Pulses

The advent of X-ray Free Electron Lasers (XFELs) has ushered in a transformative era in the field of structural biology, materials science, and ultrafast physics. These state-of-the-art facilities generate ultra-bright, femtosecond-long X-ray pulses, allowing researchers to delve into the structure and dynamics of molecular systems with unprecedented temporal and spatial resolutions. The unique properties of XFEL pulses have opened new avenues for scientific exploration that were previously considered unattainable. One of the most notable applications of XFELs is in structural biology. Traditional X-ray crystallography, while instrumental in determining the structures of countless biomolecules, often requires large, high-quality crystals and may not capture highly transient states of proteins. XFELs, with their ability to produce diffraction patterns from nanocrystals or even single particles, have provided solutions to these challenges. XFEL has expanded the toolbox of structural biologists by enabling structural determination approaches such as Single Particle Imaging (SPI) and Serial X-ray Crystallography (SFX). Despite their remarkable capabilities, the journey of XFELs is still in its nascent stages, with ongoing advancements aimed at improving their coherence, pulse duration, and wavelength tunability.

Water layer and radiation damage effects on the orientation recovery of proteins in single-particle imaging at an X-ray free-electron laser

The noise caused by sample heterogeneity (including sample solvent) has been identified as one of the determinant factors for a successful X-ray single-particle imaging experiment. It influences both the radiation damage process that occurs during illumination as well as the scattering patterns captured by the detector. Here, we investigate the impact of water layer thickness and radiation damage on orientation recovery from diffraction patterns of the nitrogenase iron protein. Orientation recovery is a critical step for single-particle imaging. It enables to sort a set of diffraction patterns scattered by identical particles placed at unknown orientations and assemble them into a 3D reciprocal space volume. The recovery quality is characterized by a "disconcurrence" metric. Our results show that while a water layer mitigates protein damage, the noise generated by the scattering from it can introduce challenges for orientation recovery and is anticipated to cause problems in the phase retrieval process to extract the desired protein structure. Compared to these disadvantageous effects due to the thick water layer, the effects of radiation damage on the orientation recovery are relatively small. Therefore, minimizing the amount of residual sample solvent should be considered a crucial step in improving the fidelity and resolution of X-ray single-particle imaging experiments.

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.

Expected resolution limits of x-ray free-electron laser single-particle imaging for realistic source and detector properties

The unprecedented intensity of x-ray free-electron laser sources has enabled single-particle x-ray diffraction imaging (SPI) of various biological specimens in both two-dimensional projection and three dimensions (3D). The potential of studying protein dynamics in their native conditions, without crystallization or chemical staining, has encouraged researchers to aim for increasingly higher resolutions with this technique. The currently achievable resolution of SPI is limited to the sub-10 nanometer range, mainly due to background effects, such as instrumental noise and parasitic scattering from the carrier gas used for sample delivery. Recent theoretical studies have quantified the effects of x-ray pulse parameters, as well as the required number of diffraction patterns to achieve a certain resolution, in a 3D reconstruction, although the effects of detector noise and the random particle orientation in each diffraction snapshot were not taken into account. In this work, we show these shortcomings and address limitations on achievable image resolution imposed by the adaptive gain integrating pixel detector noise.

High-fluence and high-gain multilayer focusing optics to enhance spatial resolution in femtosecond X-ray laser imaging

With the emergence of X-ray free-electron lasers (XFELs), coherent diffractive imaging (CDI) has acquired a capability for single-particle imaging (SPI) of non-crystalline objects under non-cryogenic conditions. However, the single-shot spatial resolution is limited to ~5 nanometres primarily because of insufficient fluence. Here, we present a CDI technique whereby high resolution is achieved with very-high-fluence X-ray focusing using multilayer mirrors with nanometre precision. The optics can focus 4-keV XFEL down to 60 nm × 110 nm and realize a fluence of >3 × 10⁵ J cm⁻² pulse⁻¹ or >4 × 10¹² photons μm⁻² pulse⁻¹ with a tenfold increase in the total gain compared to conventional optics due to the high demagnification. Further, the imaging of fixed-target metallic nanoparticles in solution attained an unprecedented 2-nm resolution in single-XFEL-pulse exposure. These findings can further expand the capabilities of SPI to explore the relationships between dynamic structures and functions of native biomolecular complexes.

Here, the authors realize an ultra-high fluence X-ray laser by high-gain multilayer focusing optics. This enables in-solution imaging with 2-nm resolution in a single-pulse exposure, making strides toward biomolecular imaging under physiological conditions.

Skopi: a simulation package for diffractive imaging of noncrystalline biomolecules

X-ray free-electron lasers (XFELs) have the ability to produce ultra-bright femtosecond X-ray pulses for coherent diffraction imaging of biomolecules. While the development of methods and algorithms for macromolecular crystallography is now mature, XFEL experiments involving aerosolized or solvated biomolecular samples offer new challenges in terms of both experimental design and data processing. Skopi is a simulation package that can generate single-hit diffraction images for reconstruction algorithms, multi-hit diffraction images of aggregated particles for training machine learning classifiers using labeled data, diffraction images of randomly distributed particles for fluctuation X-ray scattering algorithms, and diffraction images of reference and target particles for holographic reconstruction algorithms. Skopi is a resource to aid feasibility studies and advance the development of algorithms for noncrystalline experiments at XFEL facilities.

Single-Image Super-Resolution Improvement of X-ray Single-Particle Diffraction Images Using a Convolutional Neural Network

Femtosecond X-ray pulse lasers are promising probes for the elucidation of the multiconformational states of biomolecules because they enable snapshots of single biomolecules to be observed as coherent diffraction images. Multi-image processing using an X-ray free-electron laser has proven to be a successful structural analysis method for viruses. However, the performance of single-particle analysis (SPA) for flexible biomolecules with sizes ≤100 nm remains difficult. Owing to the multiconformational states of biomolecules and noisy character of diffraction images, diffraction image improvement by multi-image processing is often ineffective for such molecules. Herein, a single-image super-resolution (SR) model was constructed using an SR convolutional neural network (SRCNN). Data preparation was performed in silico to consider the actual observation situation with unknown molecular orientations and the fluctuation of molecular structure and incident X-ray intensity. It was demonstrated that the trained SRCNN model improved the single-particle diffraction image quality, corresponding to an observed image with an incident X-ray intensity (approximately three to seven times higher than the original X-ray intensity), while retaining the individuality of the diffraction images. The feasibility of SPA for flexible biomolecules with sizes ≤100 nm was dramatically increased by introducing the SRCNN improvement at the beginning of the various structural analysis schemes.

Tree-Code Based Improvement of Computational Performance of the X-ray-Matter-Interaction Simulation Tool XMDYN

In this work, we report on incorporating for the first time tree-algorithm based solvers into the molecular dynamics code, XMDYN. XMDYN was developed to describe the interaction of ultrafast X-ray pulses with atomic assemblies. It is also a part of the simulation platform, SIMEX, developed for computational single-particle imaging studies at the SPB/SFX instrument of the European XFEL facility. In order to improve the XMDYN performance, we incorporated the existing tree-algorithm based Coulomb solver, PEPC, into the code, and developed a dedicated tree-algorithm based secondary ionization solver, now also included in the XMDYN code. These extensions enable computationally efficient simulations of X-ray irradiated large atomic assemblies, e.g., large protein systems or viruses that are of strong interest for ultrafast X-ray science. The XMDYN-based preparatory simulations can now guide future single-particle-imaging experiments at the free-electron-laser facility, EuXFEL.

In a flash of light: X-ray free electron lasers meet native mass spectrometry

During the last years, X-ray free electron lasers (XFELs) have emerged as X-ray sources of unparalleled brightness, delivering extreme amounts of photons in femtosecond pulses. As such, they have opened up completely new possibilities in drug discovery and structural biology, including studying high resolution biomolecular structures and their functioning in a time resolved manner, and diffractive imaging of single particles without the need for their crystallization. In this perspective, we briefly review the operation of XFELs, their immediate uses for drug discovery and focus on the potentially revolutionary single particle diffractive imaging technique and the challenges which remain to be overcome to fully realize its potential to provide high resolution structures without the need for crystallization, freezing or the need to keep proteins stable at extreme concentrations for long periods of time. As the issues have been to a large extent sample delivery related, we outline a way for native mass spectrometry to overcome these and enable so far impossible research with a potentially huge impact on structural biology and drug discovery, such as studying structures of transient intermediate species in viral life cycles or during functioning of molecular machines.

Effects of radiation damage and inelastic scattering on single-particle imaging of hydrated proteins with an X-ray Free-Electron Laser

We present a computational case study of X-ray single-particle imaging of hydrated proteins on an example of 2-Nitrogenase-Iron protein covered with water layers of various thickness, using a start-to-end simulation platform and experimental parameters of the SPB/SFX instrument at the European X-ray Free-Electron Laser facility. The simulations identify an optimal thickness of the water layer at which the effective resolution for imaging the hydrated sample becomes significantly higher than for the non-hydrated sample. This effect is lost when the water layer becomes too thick. Even though the detailed results presented pertain to the specific sample studied, the trends which we identify should also hold in a general case. We expect these findings will guide future single-particle imaging experiments using hydrated proteins.

Ultraintense, ultrashort pulse X-ray scattering in small molecules

We examine X-ray scattering from an isolated organic molecule from the linear to nonlinear absorptive regime. In the nonlinear regime, we explore the importance of both the coherent and incoherent channels and observe the onset of nonlinear behavior as a function of pulse duration and energy. In the linear regime, we test the sensitivity of the scattering signal to molecular bonding and electronic correlation via calculations using the independent atom model (IAM), Hartree-Fock (HF) and density functional theory (DFT). Finally, we describe how coherent X-ray scattering can be used to directly visualize femtosecond charge transfer and dissociation within a single molecule undergoing X-ray multiphoton absorption.

Transient ionization potential depression in nonthermal dense plasmas at high x-ray intensity

The advent of x-ray free-electron lasers (XFELs), which provide intense ultrashort x-ray pulses, has brought a new way of creating and analyzing hot and warm dense plasmas in the laboratory. Because of the ultrashort pulse duration, the XFEL-produced plasma will be out of equilibrium at the beginning, and even the electronic subsystem may not reach thermal equilibrium while interacting with a femtosecond timescale pulse. In the dense plasma, the ionization potential depression (IPD) induced by the plasma environment plays a crucial role for understanding and modeling microscopic dynamical processes. However, all theoretical approaches for IPD have been based on local thermal equilibrium (LTE), and it has been controversial to use LTE IPD models for the nonthermal situation. In this work, we propose a non-LTE (NLTE) approach to calculate the IPD effect by combining a quantum-mechanical electronic-structure calculation and a classical molecular dynamics simulation. This hybrid approach enables us to investigate the time evolution of ionization potentials and IPDs during and after the interaction with XFEL pulses, without the limitation of the LTE assumption. In our NLTE approach, the transient IPD values are presented as distributions evolving with time, which cannot be captured by conventional LTE-based models. The time-integrated ionization potential values are in good agreement with benchmark experimental data on solid-density aluminum plasma and other theoretical predictions based on LTE. The present work is promising to provide critical insights into nonequilibrium dynamics of dense plasma formation and thermalization induced by XFEL pulses.

An encryption-decryption framework to validating single-particle imaging

We propose an encryption-decryption framework for validating diffraction intensity volumes reconstructed using single-particle imaging (SPI) with X-ray free-electron lasers (XFELs) when the ground truth volume is absent. This conceptual framework exploits each reconstructed volumes' ability to decipher latent variables (e.g. orientations) of unseen sentinel diffraction patterns. Using this framework, we quantify novel measures of orientation disconcurrence, inconsistency, and disagreement between the decryptions by two independently reconstructed volumes. We also study how these measures can be used to define data sufficiency and its relation to spatial resolution, and the practical consequences of focusing XFEL pulses to smaller foci. This conceptual framework overcomes critical ambiguities in using Fourier Shell Correlation (FSC) as a validation measure for SPI. Finally, we show how this encryption-decryption framework naturally leads to an information-theoretic reformulation of the resolving power of XFEL-SPI, which we hope will lead to principled frameworks for experiment and instrument design.

Impact of real mirror profiles inside a split-and-delay unit on the spatial intensity profile in pump/probe experiments at the European XFEL

For the High-Energy-Density (HED) beamline at the SASE2 undulator of the European XFEL, a hard X-ray split-and-delay unit (SDU) has been built enabling time-resolved pump/probe experiments with photon energies between 5 keV and 24 keV. The optical layout of the SDU is based on geometrical wavefront splitting and multilayer Bragg mirrors. Maximum delays between Δτ = ±1 ps at 24 keV and Δτ = ±23 ps at 5 keV will be possible. Time-dependent wavefront propagation simulations were performed by means of the Synchrotron Radiation Workshop (SRW) software in order to investigate the impact of the optical layout, including diffraction on the beam splitter and recombiner edges and the three-dimensional topography of all eight mirrors, on the spatio-temporal properties of the XFEL pulses. The radiation is generated from noise by the code FAST which simulates the self-amplified spontaneous emission (SASE) process. A fast Fourier transformation evaluation of the disturbed interference pattern yields for ideal mirror surfaces a coherence time of τc = 0.23 fs and deduces one of τc = 0.21 fs for the real mirrors, thus with an error of Δτ = 0.02 fs which is smaller than the deviation resulting from shot-to-shot fluctuations of SASE2 pulses. The wavefronts are focused by means of compound refractive lenses in order to achieve fluences of a few hundred mJ mm-2 within a spot width of 20 µm (FWHM) diameter. Coherence effects and optics imperfections increase the peak intensity between 200 and 400% for pulse delays within the coherence time. Additionally, the influence of two off-set mirrors in the HED beamline are discussed. Further, we show the fluence distribution for Δz = ±3 mm around the focal spot along the optical axis. The simulations show that the topographies of the mirrors of the SDU are good enough to support X-ray pump/X-ray probe experiments.

Effect of radiation damage and illumination variability on signal-to-noise ratio in X-ray free-electron laser single-particle imaging

The deterioration of both the signal-to-noise ratio and the spatial resolution in the electron-density distribution reconstructed from diffraction intensities collected at different orientations of a sample is analysed theoretically with respect to the radiation damage to the sample and the variations in the X-ray intensities illuminating different copies of the sample. The simple analytical expressions and numerical estimates obtained for models of radiation damage and incident X-ray pulses may be helpful in planning X-ray free-electron laser (XFEL) imaging experiments and in analysis of experimental data. This approach to the analysis of partially coherent X-ray imaging configurations can potentially be used for analysis of other forms of imaging where the temporal behaviour of the sample and the incident intensity during exposure may affect the inverse problem of sample reconstruction.

Recovery of undamaged electron-density maps in the presence of damage-induced partial coherence in single-particle imaging

Resolving the electronic structure of single biological molecules in their native state was among the primary motivations behind X-ray free-electron lasers. The ultra-short pulses they produce can outrun the atomic motion induced by radiation damage, but the electronic structure of the sample is still significantly modified from its original state. This paper explores the decoherence of the scattered signal induced by temporal evolution of the electronic structure in the sample molecule. It is shown that the undamaged electron density of a single-molecule sample can often be retrieved using only the two most occupied modes from the coherent mode decomposition of the partially coherent diffraction fluence.

Coarse-Grained Diffraction Template Matching Model to Retrieve Multiconformational Models for Biomolecule Structures from Noisy Diffraction Patterns

Biomolecular imaging using X-ray free-electron lasers (XFELs) has been successfully applied to serial femtosecond crystallography. However, the application of single-particle analysis for structure determination using XFELs with 100 nm or smaller biomolecules has two practical problems: the incomplete diffraction data sets for reconstructing 3D assembled structures and the heterogeneous conformational states of samples. A new diffraction template matching method is thus presented here to retrieve a plausible 3D structural model based on single noisy target diffraction patterns, assuming candidate structures. Two concepts are introduced here: prompt candidate diffraction, generated by enhanced sampled coarse-grain (CG) candidate structures, and efficient molecular orientation searching for matching based on Bayesian optimization. A CG model-based diffraction-matching protocol is proposed that achieves a 100-fold speed increase compared to exhaustive diffraction matching using an all-atom model. The conditions that enable multiconformational analysis were also investigated by simulated diffraction data for various conformational states of chromatin and ribosomes. The proposed method can enable multiconformational analysis, with a structural resolution of at least 20 Å for 270-800 Å flexible biomolecules, in experimental single-particle structure analyses that employ XFELs.

Structural dynamics in proteins induced by and probed with X-ray free-electron laser pulses

X-ray free-electron lasers (XFELs) enable crystallographic structure determination beyond the limitations imposed upon synchrotron measurements by radiation damage. The need for very short XFEL pulses is relieved through gating of Bragg diffraction by loss of crystalline order as damage progresses, but not if ionization events are spatially non-uniform due to underlying elemental distributions, as in biological samples. Indeed, correlated movements of iron and sulfur ions were observed in XFEL-irradiated ferredoxin microcrystals using unusually long pulses of 80 fs. Here, we report a femtosecond time-resolved X-ray pump/X-ray probe experiment on protein nanocrystals. We observe changes in the protein backbone and aromatic residues as well as disulfide bridges. Simulations show that the latter’s correlated structural dynamics are much slower than expected for the predicted high atomic charge states due to significant impact of ion caging and plasma electron screening. This indicates that dense-environment effects can strongly affect local radiation damage-induced structural dynamics.

The local X-ray-induced dynamics that occur in protein crystals during serial femtosecond crystallography (SFX) measurements at XFELs are not well understood. Here the authors performed a time-resolved X-ray pump X-ray probe SFX experiment, and they observe distinct structural changes in the disulfide bridges and peptide backbone of proteins; complementing theoretical approaches allow them to further characterize the details of the X-ray induced ionization and local structural dynamics.

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.

X-ray free electron lasers provide high photon flux to explore single particle diffraction imaging of biological samples. Here the authors present dynamic electronic structure calculations and benchmark them to single-particle XFEL diffraction data of sucrose clusters to predict optimal single-shot imaging conditions.

Is radiation damage the limiting factor in high-resolution single particle imaging with X-ray free-electron lasers?

The prospect of single particle imaging with atomic resolution is one of the scientific drivers for the development of X-ray free-electron lasers. The assumption since the beginning has been that damage to the sample caused by intense X-ray pulses is one of the limiting factors for achieving subnanometer X-ray imaging of single particles and that X-ray pulses need to be as short as possible. Based on the molecular dynamics simulations of proteins in X-ray fields of various durations (5 fs, 25 fs, and 50 fs), we show that the noise in the diffracted signal caused by radiation damage is less than what can be expected from other sources, such as sample inhomogeneity and X-ray shot-to-shot variations. These findings show a different aspect of the feasibility of high-resolution single particle imaging using free-electron lasers, where employing X-ray pulses of longer durations could still provide a useful diffraction signal above the noise due to the Coulomb explosion.

Experimental 3D coherent diffractive imaging from photon-sparse random projections

The routine atomic resolution structure determination of single particles is expected to have profound implications for probing structure-function relationships in systems ranging from energy-storage materials to biological molecules. Extremely bright ultrashort-pulse X-ray sources - X-ray free-electron lasers (XFELs) - provide X-rays that can be used to probe ensembles of nearly identical nanoscale particles. When combined with coherent diffractive imaging, these objects can be imaged; however, as the resolution of the images approaches the atomic scale, the measured data are increasingly difficult to obtain and, during an X-ray pulse, the number of photons incident on the 2D detector is much smaller than the number of pixels. This latter concern, the signal 'sparsity', materially impedes the application of the method. An experimental analog using a conventional X-ray source is demonstrated and yields signal levels comparable with those expected from single biomolecules illuminated by focused XFEL pulses. The analog experiment provides an invaluable cross check on the fidelity of the reconstructed data that is not available during XFEL experiments. Using these experimental data, it is established that a sparsity of order 1.3 × 10-3 photons per pixel per frame can be overcome, lending vital insight to the solution of the atomic resolution XFEL single-particle imaging problem by experimentally demonstrating 3D coherent diffractive imaging from photon-sparse random projections.

Incorporating particle symmetry into orientation determination in single-particle imaging

In coherent-diffraction-imaging experiments X-ray diffraction patterns of identical particles are recorded. The particles are injected into the X-ray free-electron laser (XFEL) beam in random orientations. If the particle has symmetry, finding the orientation of a pattern can be ambiguous. With some modifications, the correlation-maximization method can find the relative orientations of the diffraction patterns for the case of symmetric particles as well. After convergence, the correlation maps show the symmetry of the particle and can be used to determine the symmetry elements and their orientations. The C factor, slightly modified for the symmetric case, can indicate the consistency of the assembled three-dimensional intensity distribution.