IUCrJ

Time-resolved crystallography can be inaccessible for membrane proteins and other difficult-to-crystallize targets for which suitable crystals are scarce. Here we present an acquisition strategy that records consecutive diffraction images during a small crystal rotation synchronized to a single activation pulse. Matching image segments are merged across crystals to reconstruct successive time windows, enabling a single crystal ensemble to capture an extended structural trajectory. Using the Mesh&Collect workflow at a synchrotron, we obtained ten consecutive 7.5-ms time windows at 2 [A] resolution from 1,255 microcrystals of a bacteriorhodopsin triple mutant at room temperature. An additional 790 crystals extended the observation window to 500 ms. Rotation-based processing, retrospective time binning and dose control maximize information extracted from each crystal, while compatible measurements can be pooled across acquisitions. This approach makes millisecond time-resolved crystallography practical for crystal-limited targets at synchrotrons, with further development potentially extending the temporal resolution down to 1 ms.

X-ray fiber diffraction patterns exhibit four-fold symmetry that can be exploited, through folding and averaging, to improve signal-to-noise ratio. Accurate folding requires a precise sub-pixel estimate of the symmetry center and precise orientation of the meridional pattern axis to the fiber axis: small center or angular errors blur diffraction features, reduce layer-line sharpness, and introduce errors in spacing measurements. A pixel-level estimate is often too imprecise for this purpose, and detector gaps further complicate the alignment objective. We formulate the masked quadrant-folding problem, define a four-quadrant symmetry loss that consistently excludes invalid pixels, and evaluate several refinement strategies: hierarchical coarse-to-fine grid search; ECC-based rigid registration with global center/orientation correction fitting; ECC registration followed by local gradient refinement; and a hybrid that appends a local grid search on a cropped pattern. Direct gradient optimization from the rough QF alignment was found to be unreliable. Grid search provides a robust, interpretable baseline that directly minimizes the folding objective but is substantially slower than registration; ECC gives a fast near-correct alignment, and the hybrid closes the accuracy gap to brute-force search at a fraction of its runtime. On real datasets with calibration data, applying a calibration center with optimized rotation is effectively optimal. The hybrid center-refinement method has been integrated into the MuscleX package.

Damage to biological specimens by the electron beam is the fundamental resolution-limiting factor in cryoelectron microscopy (cryo-EM) single particle analysis. There is, however, currently no method to accurately infer fluence-dependent changes to the specimen structure during electron irradiation. We develop a Bayesian framework to fit a sequence of atomic models to a series of cryo-EM reconstructions produced at increasing fluence. In particular, our algorithm is able to infer the ensemble average position and atomic displacement parameter of every atom in the macromolecule as a function of fluence. Application of the algorithm to cryo-EM datasets shows that the molecule expands during imaging and identifies environment-dependent variations in beam-induced damage. We use our results to propose a stochastic process model of this phenomenon. We envisage that our method will lead to a better mechanistic understanding of radiation damage to biological specimens and may contribute to efforts to mitigate its effects.

Three-dimensional electron diffraction (3DED) is an emerging and useful technique for solving molecular structures of small and biological macro-molecules from nanometre-sized crystals. We present our automated data processing workflow for 3DED datasets collected at Diamond Light Sources electron Bio-Imaging Centre (eBIC). For this purpose, we developed a package called AutoED. The processing pipeline includes data collection, analysis of the beam position, metadata gathering, file conversion, and finally data processing using xia2 (which supports both DIALS and XDS). The processing results are captured in a summary report produced by AutoED. Our main goal is to reduce the workload of electron diffraction scientists, but also to enforce good standards already used in macromolecular crystallography (MX). All the collected 3DED datasets are automatically converted into NeXus data format which is considered a Gold Standard for MX. This standardized data format allows for all the relevant metadata about the experiment to be kept together with diffraction images. We also discuss the methods used in AutoED to determine the electron beam position on diffraction images.

Visualizing macromolecules within their native cellular context by cryo-electron tomography (cryo-ET) is fundamentally limited by the trade-off between field of view and resolution: Capturing high-resolution information about biomolecules requires high magnification, which restricts the field of view and obscures the cellular context in which those biomolecules function. Collecting montage data by tiling the electron beam over the region of interest offers one solution, although traditional round electron beams cause excessive radiation damage across overlapping regions. We previously made electron beams square in shape, enabling montage collection with minimal overlap and thereby reducing excessive exposure and loss of high-resolution information. Here, we create a pipeline for collecting and processing montage cryo-ET data with square electron beams. We show that square beam montages retain high-resolution information by reconstructing virus-like particles to 3.5 [A] resolution using sub-tomogram averaging, and apply the workflow to imaging a glial cell and malaria parasite lamellae over fields of view up to 65 m2. We also provide a comprehensive protocol to make square beams accessible to the community.

The visualisation of large biological samples in cryo-electron tomography requires both few-nanometre resolution and a wide, ideally >10 micrometrre, field of view. For a single image, increasing magnification trades increased resolution for a smaller field of view. Montaging, where high-magnification images are stitched together to increase field of view, addresses the trade-off but is dose-inefficient with round beams that are standard in transmission electron microscopy (TEM), since this illumination shape does not tessellate. Our workflow greatly simplifies montaging schemes for dose-sensitive specimens, allowing large montages of any shape - from individual cells or other features of interest up to whole lamellae - to be collected at high resolution without compromising electron dose efficiency. We demonstrate the workflow by acquiring montage tomograms of whole yeast (Saccharomyces cerevisiae) cells in cryo-FIB lamellae and of human neurons.

Time-resolved serial crystallography (TR-SX) has emerged as a powerful method for capturing ultrafast structural dynamics in proteins. TR-SX continues to produce remarkable studies, revealing previously unobserved transient states and providing deeper insights into processes such as drug targeting, DNA repair, and photosynthesis. However, extracting weak structural signals from noisy time-resolved datasets remains a major challenge. Robust computational methods are therefore required to isolate the signals associated with the underlying transient states. Importantly, this should be performed in reciprocal space to preserve compatibility with established downstream structure refinement workflows. Here, we introduce a framework for kinetic decomposition directly in reciprocal space that enables separation of kinetically distinct structural states. The method decomposes crystallographic data according to a predefined kinetic model, improving the recovery of weak transient signals and enhancing mechanistic interpretation from limited time-resolved datasets. We validate the framework using simulated data based on a previously published time-resolved crystallography study and demonstrate its application to a new TR-SX dataset comprising 17 time points. We show that the method separates the reciprocal space signatures of four intermediates by incorporating kinetic information from a predefined reaction model. This establishes a workflow for extracting kinetic states directly from time-resolved X-ray diffraction data that can be seamlessly integrated into existing crystallographic structure-determination pipelines.

Microtubules are cytoskeletal filaments that are involved in intracellular transport, cell division, and motility. Despite their biological importance, determining their high-resolution structures via cryo-electron microscopy remains a significant technical challenge due to their polymorphisms and pseudo-helical assembly. Current processing workflows are complex, often requiring the integration of multiple software packages and custom scripts, which creates a steep learning curve for many research groups. To address these limitations, we introduce CsMT, a streamlined workflow implemented entirely within the CryoSPARC environment and using synthetic references. CsMT simplifies microtubule reconstruction by utilizing a novel protofilament-pair classification approach, which effectively handles the inherent pseudo-symmetry and structural heterogeneity of microtubules with minimal manual intervention. Our workflow is versatile, capable of processing both undecorated and decorated microtubules while accurately determining seams and performing high-resolution refinement. We demonstrate the efficacy of this workflow by achieving a 2.3 and 2.7 [A] resolution reconstruction of homotypic and heterotypic maps of undecorated microtubules, matching the best-resolved microtubule structures in the field. By unifying the pipeline into a single and portable workflow, CsMT enhances reproducibility and accessibility, empowering more laboratories to explore the structural biology of microtubules and associated proteins, yielding new insights into their function.

Solution X-ray scattering provides unique structural information on biomolecules under biological conditions, resolving conformational heterogeneity and time-resolved structural changes. The scattering profiles contain limited information, and interpretation largely relies on fitting candidate structures guided by priors. Direct reconstruction of electron density maps is desirable, but so far has been prevented by the difficulty of incorporating such prior knowledge. Here we propose XSSDense, a framework that couples a variational autoencoder trained on electron densities from predicted or simulated protein ensembles with a genetic algorithm to refine densities against scattering data. We validate XSSDense on synthetic data for crambin, recover the conformational heterogeneity of the unfolded state of Avena sativa light-oxygen-voltage sensing domain 2, resolve a de-novo density for the pre-unfolding state of the same protein, and provide a new structural description of the signalling-state ensemble of photoactive yellow protein. XSSDense enables structurally grounded electron density reconstructions that intrinsically capture conformational heterogeneity.

Proteins naturally adopt multiple conformations in mediating cellular processes, and ensemble models are used to fit X-ray crystallography data that captures this heterogeneity. In practice, ensemble refinement produces only minor improvements in agreement with experimental data (R-free) over single-conformation models. It has recently been shown that ensemble models are universally trapped, or "tangled"; refinement algorithms strain each individual conformation in the model to fit the electron density in its immediate vicinity, missing more harmonious ways to arrange the collection of protein conformations to fit the electron density. Here, we demonstrate that this type of trap may be escaped by formulating the construction of low-energy conformations from individual conformer coordinates as an integer linear programming problem. The method successfully recovered the two original protein conformations from a previously published synthetic dataset that traps current refinement methods. Inclusion of the method in an automated refinement procedure with real data is shown to improve R-free and reduce geometric strain in a four-conformation model by comparison with controls. Applying this method in combination with human input and fitting low-occupancy waters to density features in the bulk solvent, we produce models for deposited datasets of three separate 14-19 kDa proteins with greatly improved geometry and R-factors. This includes a 0.77 [A] six-conformation model of the SARS-CoV-2 macrodomain Mac1 (PDB ID: 44PS) with an R-work of 4.7% and an R-free of 6.4%.

Proteins are molecular machines whose dynamics govern their function. Using time-resolved crystallography, the basis of protein mechanics can be uncovered by physical perturbation on protein crystals. Perturbations are often vector fields that, in crystals, break symmetry. Lacking existing tools for analyzing symmetry breaking, we built open-source software to do so. We applied this method to three time-resolved crystallography experiments and analyzed the effect of vectorial perturbations on macromolecular crystal structure. We demonstrate that our method can handle symmetry reduction in all 65 Sohncke space groups.

Cryo-electron tomography of 100-200 nm lamellae prepared by focused ion beam milling enables structural analysis of macromolecules within cells, tissues, and small organisms, but lamella preparation remains time-consuming and limits throughput. Although automated workflows enable unsupervised fabrication, they commonly produce variable results because fixed, pre-calibrated parameters are applied without accounting for lamella-to-lamella variability. To overcome these limitations, we developed an Adaptive Milling approach that adapts milling parameters based on imaging feedback. We implemented this approach for the critical polishing stage as Adaptive Polishing, which uses machine-learning-based segmentation of scanning electron microscopy images to identify lamella features predictive of imminent collapse. These are used to periodically assess whether to continue or stop thinning during polishing. We performed a systematic comparison of the quality of lamellae prepared with Adaptive Polishing and conventional automation, and demonstrate that Adaptive Polishing produces thinner lamellae with greater consistency in thickness, and performs well for in situ cryo-electron tomography. Adaptive Polishing was tested on three different microscopes with minimal model retraining required, demonstrating transferability of the approach. We implemented Adaptive Polishing as an open-source plugin for fibsemOS together with machine-learning model training software and pre-trained models.

Atomistic simulations often struggle to capture the full conformational landscape of biomolecular complexes on experimentally relevant timescales, leaving gaps in our understanding of functionally relevant thermodynamics and kinetics. Cryo-electron microscopy (cryo-EM) provides high-resolution structural snapshots spanning multiple conformational states, but translating these static images into dynamic mechanisms remains challenging. Existing approaches for integrating cryo-EM data with molecular simulations either perturb the underlying dynamics through artificial biasing potentials or remain confined to conformations already present in the prior ensemble. We introduce cryoWEight, an automated and extensible Python framework that couples iterative weighted ensemble (WE) simulations with ensemble reweighting against cryo-EM particle images to steer simulation trajectories towards conformational states consistent with cryo-EM populations. We have found that reweighting in the image space introduces a directional bias towards the target distribution even when the prior ensemble has no overlap with the target distribution in the collective variable space. The reweighting requires only the state likelihoods that set the ensemble weights, and any observable that defines these likelihoods can serve as the target distribution, therefore generalizing beyond the application of cryo-EM. cryoWEight is validated on three proteins: chignolin, the N-terminal domain of ribosomal protein L9, and adenylate kinase. The framework converges towards the target conformational states across different sampling strategies and noise levels, providing a scalable route to integrate cryo-EM data with enhanced sampling approaches.

Accurate assessment of ligand coordinate-density consistency across different resolutions remains challenging in macromolecular crystallography. We introduce the atomic Box Correlation Coefficient (aBCC), an atom-level metric for evaluating the consistency between ligand atomic coordinates and electron density in a resolution-standardized framework. To predict aBCC values from electron-density maps, we developed QAEmap, a machine-learning model based on three-dimensional convolutional neural networks (3D-CNNs). The model was trained using Fourier-truncated electron-density maps and corresponding ligand coordinates generated from high-resolution structures in the Protein Data Bank. It was evaluated using both Fourier-truncated electron-density maps and experimentally determined PDB structures. was evaluated using both Fourier-truncated electron-density maps and experimentally determined PDB structures.The prediction accuracy gradually decreased with decreasing resolution, but remained reliable up to [~]3.5 [A]. These results demonstrate that aBCC enables resolution-standardized atom-wise evaluation of coordinate-density consistency across different resolutions and provide a foundation for further development and refinement of machine learning-based coordinate validation. SynopsisWe introduce the atomic box correlation coefficient (aBCC), a machine learning-based metric for the resolution-standardized atom-level evaluation of ligand coordinate-density consistency in crystallographic structures. aBCC provides a common framework for assessing and communicating the local coordinate reliability between structural biologists and researchers in structure-based drug discovery.

Single-particle cryo-electron microscopy (cryo-EM) pose estimation is traditionally solved anew for each dataset, where iterative refinement is done from scratch while the estimator learns to store the molecule in its weights. In this work, we show that pose inference is a generalizable, specimen-agnostic operation when conditioned explicitly on a reference volume. We introduce ARCHER, an amortized contrastive classifier that models the pose posterior over a discrete rotation grid. Trained across a variety of protein structures, it operates zero-shot without retraining per structure. This transferability is grounded in Fourier-space information mechanics, where all specimen dependence is captured by the reference structure's power spectrum and spatial extent. ARCHER achieves a median angular error of 5.0{degrees} on 100 held-out test structures and 2.5{degrees} on experimental particles, matching dedicated estimators within 0.16[A] in 3D reconstruction. Crucially, downstream conformational signal is preserved. The leading conformational coordinate correlates at 0.97 with deposited benchmarks, faithfully reconstructing free-energy basins and mobile domains. These results overall demonstrate that cryo-EM pose estimation can be generalized across different structures.

The ability of cryo-electron microscopy (cryo-EM) to interrogate the atomic structure and motion of biomolecules has motivated the development of a wide range of algorithms for heterogeneity analysis. However, evaluating and comparing these methods remains challenging because ground-truth structures are generally unknown for experimental samples. Here, we introduce the 2026 Community-Wide Assessment of Heterogeneous Reconstruction Algorithms (CAHRA), a community-wide challenge centered on three benchmark datasets that probe distinct tasks in heterogeneity analysis. These include (1) compositional heterogeneity arising from mixtures of distinct protein complexes, (2) continuous conformational variability and atomic modeling, and (3) entanglement between molecular conformation and particle pose. We describe the design, construction, and validation of these datasets, as well as their use in the recently completed CAHRA Challenge.

Cryo-electron ptychography is an emerging technique for studying radiation-sensitive biological specimens, developed from four-dimensional transmission electron microscopy (4D-STEM). Although ptychography has achieved ultrahigh resolution beyond conventional transmission electron microscopy limits for radiation-resistant samples, its application to frozen hydrated biological specimens currently remains at sub-nanometer resolution. Here we overcome this limitation by implementing electron counting with a hybrid-pixel detector, establishing key technical foundations for near-atomic-resolution cryo-ptychography. This counting approach significantly improves weak signal detection in convergent-beam electron diffraction, enabling ptychographic reconstruction at doses below 1 e-/[A]{superscript 2}. Additionally, we found beam-induced motion is effectively eliminated within single scans, suggesting conventional cryoEMs dose-fractionation approach may need reevaluation. Demonstrating high contrast under both low-dose and tilted conditions, along with achieving 3.59 [A] resolution for the [~]700 kDa T20S proteasome, we validate cryo-electron ptychography as a viable general imaging modality that could complement or surpass conventional phase-contrast cryoEM methods.

Directional anisotropy in electron density provides key information about chemical bonding that is not readily accessible from conventional electron-density maps. Here, a model-independent framework is presented for decomposing experimental structure factors into angular components using spherical harmonics. Reciprocal-space projection onto spherical harmonics followed by standard Fourier synthesis yields angularly filtered density maps. The{ell} = 0 component captures the isotropic part of the density, while the{ell} = 1 components resemble px, py and pz-like dipolar functions that highlight directional electronic structure. Applications to high-resolution datasets, including urea, the Gly-Ala dipeptide and a 0.97 [A]{beta}-lactamase structure, reveal chemically interpretable dipolar features associated with carbonyl and amide bonds, N-H interactions and aromatic{pi} systems. Quantitative analysis using bond-centred sampling demonstrates stable dipolar signatures that remain detectable under moderate resolution truncation. These results establish spherical-harmonic angular decomposition as a practical framework for extracting directional electronic information from crystallographic electron-density maps. SynopsisAngular decomposition of experimental structure factors reveals dipolar anisotropy and directional electron-density features that are directly meaningful for chemical interpretation.

Here, we present MC-Bayes, a Python-based script for processing cryo-electron microscopy EER movies on one or more GPUs using MotionCor3 in a user-friendly manner. Further, it generates the .star files necessary for RELION to perform Bayesian polishing (a.k.a.: reference-based motion correction) with EER movies. Until now, Bayesian polishing of EER data was only possible if the CPU-based "RELIONCor" implementation of MotionCor2 was used, which is sub-optimal on GPU-heavy cryo-EM processing systems. This wrapper was created for those facilities and/or users who may have (many) powerful GPUs, but for whatever reason have few CPU cores or less system RAM. Leveraging MotionCor3, MC-Bayes allows motion correction of EER data 2 or more times faster (depending on system) than the RELION CPU implementation, except in circumstances where dozens or hundreds of CPU cores with high quantities of system RAM can be utilised.

Preferred particle orientation, which results from the adsorption of proteins to the air-water interface, is a frequently encountered issue in single-particle cryo-electron microscopy (cryo-EM). Particles then present only a limited number of views, which are oversampled during data acquisition at the expense of throughput, while other views may be missing altogether. In severe cases, a high-resolution reconstruction becomes unattainable. We have recently established that ultrasonic excitation of the sample during vitrification improves the angular distribution by detaching particles from the air-water interface and reshuffling their orientations. Here, we describe an implementation of this principle for a plunge freezing apparatus, in which the ultrasonic transducer is integrated into the head of the pair of tweezers holding the sample, and show that this design is able to reduce preferred orientation. We demonstrate this for hemagglutinin, the canonical example of preferred orientation, and show that the fraction of side and tilted views increases with the excitation amplitude. This allows us to obtain a complete reconstruction of the trimer at 2.5 [A].