IUCrJ

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.

Cryo-electron microscopy (cryo-EM) of biological specimens is limited by radiation damage and a low signal-to-noise ratio (SNR). Here, we show that reducing the illuminated area substantially slows the observed diffraction decay in protein microcrystals. We further show that narrow parallel-beam electron diffraction from thin non-crystalline biological specimens provides substantially higher reciprocal-space SNR than conventional cryo-EM imaging. We developed a multimodal scanning workflow, 4D-para-STEM, that records narrow-beam diffraction patterns together with corresponding images. Using viruses, peptide assemblies, and microtubules, we demonstrate interpretable diffraction signals from both crystalline and non-crystalline biological specimens. Together, these results show that narrow parallel-beam scanning reduces observed radiation damage and improves the SNR in cryo-EM.

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.

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.

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.

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.

Cryo-electron microscopy (cryo-EM) has become a central tool for structure-based drug discovery, yet ligand-binding sites often remain substantially less well resolved than the surrounding protein, limiting reliable atomic interpretation. Although deep-learning methods have substantially improved overall cryo-EM map quality, their predominantly protein-focused training limits their ability to recover ligand density. Here we present CryoLigATE, a deep learning framework specifically designed to enhance densities associated with protein-bound ligands in cryo-EM maps. We curated a chemically and structurally diverse dataset of more than 6,000 protein-ligand complexes from the EMDB and PDB, encompassing drug-like molecules, lipids, steroids, carbohydrates and other ligand classes, and trained a hybrid convolutional-transformer network to enhance local density around binding pockets. During inference, CryoLigATE automatically extracts the target region from a preliminary atomic model, requiring no manual map preparation and completing localized refinement in seconds on a desktop GPU. Evaluation on an independent test set of 649 complexes demonstrates substantial improvements in ligand resolvability, particularly for maps with poorly resolved binding sites, while preserving high-quality experimental densities. The enhanced maps recover chemically meaningful features, including ligand functional groups and topological continuity, enabling more confident atomic modeling. By learning the structural diversity of ligand features, CryoLigATE addresses a longstanding limitation of cryo-EM map enhancement and provides a useful framework for improving structural interpretation and structure-guided drug discovery.

Accurate estimation of the contrast transfer function (CTF) of tilt images is a critical first step in cryo electron tomography (cryoET), enabling reliable recovery of high-resolution structural information from thick, heterogeneous specimens. This challenge is especially acute in in situ cryoET, where macromolecules are imaged in their native cellular environment, often at high tilt and through substantial specimen thickness, with correspondingly low signal-to-noise ratios. Although CTF parameters can be later refined using reference-based approaches, accurate initial estimates are critical for downstream processing and the interpretability of tomographic reconstructions, yet they remain difficult to automate. Here, we present Quinoa, a software package designed to address these challenges. Quinoa first validates the tilt geometry and assesses data quality to generate robust initial estimates of defocus and phase shift. These estimates are then refined through optimization of a single global model, enabling precise fitting of the per-image defoci, tilt-dependent astigmatisms, time-dependent phase shifts, the specimen orientation (rotation, tilt and pitch) and the specimen thickness. Notably, and as a key distinguishing feature of this approach is that Quinoa fits equiphase-binned polar power spectra. This substantially reduces the computational cost of optimization without sacrificing accuracy, enabling more progressive and exhaustive refinement passes that further improve robustness. We validated Quinoa using both simulated and experimental data and benchmarked its performance against Warp, Ctfplotter, CTFMeasure, and AreTomo. Our results show that Quinoa is the most robust approach across all simulated cases, maintaining high accuracy even in the simultaneous presence of severe astigmatism, high specimen inclination and variable phase shift. Integrated recovery mechanisms further allow Quinoa to adapt automatically to a wide range of pixel sizes, defoci, astigmatisms and specimen thicknesses. Despite fitting a more complex and dynamic model, Quinoa remains extremely efficient due to extensive GPU acceleration, making it well suited for real-time monitoring during data collection as well as high-throughput offline batch processing. By improving automated CTF estimation in challenging tomographic data, Quinoa supports more accurate structural analysis of cells and tissues in situ.

Images in the electron microscope are formed by electron scattering and focusing. The spherical geometry of these processes gives rise to two coherent, conjugate spherical wave fronts, known as Ewald spheres. These spheres are associated with the two halves of the contrast transfer function (CTF), and their widths are determined by the focal gradient through the specimen. To properly correct for the CTF, each half of the CTF must be applied to an image individually and integrated into the reconstruction into the corresponding Ewald sphere. Theory indicates that this dual Ewald sphere reconstruction method should recover the maximal amount of information possible. This method was compared to the other reconstruction methods commonly used: the projection approximation (ignoring the Ewald sphere), the simple insertion and the single sideband methods. In simulated reconstructions the dual Ewald sphere method recovered the most information when the correct half of the CTF is matched to the corresponding Ewald sphere. If the wrong half is matched, the result worse than the projection approximation method. Examining reconstructions from real data indicated that the dual Ewald sphere method performs at least as well as the simple insertion method, but not as good as in simulations. The likely reason is the two-fold ambiguity in the assigned orientations of the particle images, which remains an issue to pursue in further studies. In conclusion, the dual Ewald sphere reconstruction method may offer the best way to calculate very high resolution reconstructions when the micrograph quality warrants it. HighlightsO_LIThe dual Ewald sphere reconstruction corrects for the two halves of the CTF. C_LIO_LIThe signs of the two halves of the CTF must correspond to the focal gradient. C_LIO_LIDetermining the focal gradient for individual particle images remains unresolved. C_LIO_LIComplex reconstructions indicate any real space phases are artifacts. C_LI

Cryogenic electron microscopy (cryoEM) is now routinely used for high-resolution structure determination of biological macromolecules. However, many biological specimens exhibit varying degrees of preferred orientation on cryoEM grids, resulting in uneven sampling of three-dimensional Fourier space. This orientation bias produces anisotropic reconstruction artifacts and, in severe cases, can exacerbate particle misalignment during iterative refinement, thereby limiting the success rate of near-atomic resolution cryoEM structure determination. This protocol provides a practical guide for applying spIsoNet, a self-supervised deep-learning method, to mitigate preferred-orientation issues in cryoEM reconstructions. We describe two complementary workflows: (1) map Anisotropy Correction to correct anisotropic artifacts of cryoEM maps and (2) particle Misalignment Correction, which integrates spIsoNet with RELION external reconstruction to improve particle-pose estimation. We demonstrate these workflows using two influenza hemagglutinin (HA) trimer datasets representing moderate and severe degrees of preferred-orientation bias. The protocol includes installation instructions, parameter-selection guidance, quality-control checkpoints and troubleshooting advice, and can typically be completed in ~7 hours on a workstation equipped with four NVIDIA A100 GPUs. Together, these workflows provide step-by-step guidance for using the open-source spIsoNet software to mitigate the preferred-orientation problem directly from experimental data.

The quality of single particle cryo-EM reconstructions can be severely degraded when an insufficient variety of 3D particle orientations is present in the image data, limiting downstream model building and interpretation. However, it is often difficult to ascertain whether or not a particular dataset suffers from such preferred orientation since the required orientation coverage depends on target geometry, alignment accuracy, and particle quality. To simplify diagnosis of preferred orientation, we present two complementary methods. First, the conical Fourier Shell Correlation Area Ratio (cFAR) compares the worst- and best-correlating conical regions of 3D Fourier space to quantify half-map anisotropy into a single, easily interpretable score ranging from zero to one. Second, Relative Signal, a companion to cFAR, directly relates signal content to viewing direction so that under-sampled views can be identified. We characterize our methods and compare them to existing anisotropy detection approaches on synthetic data and on 14 real datasets that span sundry molecular weights and structure types. Implementations of both cFAR and Relative Signal are included in CryoSPARC v4.5 and later versions.

Four-dimensional scanning transmission electron microscopy (4D-STEM) enables the acquisition of diffraction patterns at every probe position in a dense array. For imaging applications this approach offers significant benefits in terms of spatial resolution and contrast enhancement. In this work, we present the development of a synchronous scan generator integrated with SerialEM software to enable automation of complex experimental protocols such as tomography. The proposed hardware functions as an interface between SerialEM, the scan controls of the microscope, a fast annular dark-field detector, and a synchronized trigger for a pixelated detector. Our previous implementation, named SavvyScan, relied on a dedicated computer equipped with a multichannel acquisition and signal-generation cards, as well as a separate microcontroller for synchronization. Here, we report a low-cost implementation based on a Red Pitaya board, utilizing direct programming of its embedded FPGA and Linux server components. We provide detailed instructions for system installation and operation, along with practical guidance for modifying the source code. System performance is validated through oscilloscope measurements and imaging of a replica grating sample. The utility of the approach is further demonstrated by generating a 3D electron tomogram of a cryogenic sample of mitochondria from a tilt series of shadow montage projections.

Ordered water molecules mediate many protein functions, including stability, ligand binding, and catalysis. Predicting their positions with sub-angstrom accuracy would support protein design, binding affinity prediction, and automated model building in X-ray crystallography and cryo-EM. However, water molecule prediction lags behind protein and other molecule structure predictions. Here, we introduce WaterFlow, a flow-matching-based generator model and confidence model for predicting the positions of ordered water molecules in protein structures. WaterFlow outperforms the existing state of the art at every precision level. We demonstrate that WaterFlow can accurately predict ground truth modeled water molecules, including those around protein-ligand interactions and on predicted structures. We also show that WaterFlow predictions fit well directly to experimental data, and therefore propose that it may be used for both prediction and modeling water molecules. This includes novel predictions that are often associated with positive electron difference density, meaning the model places water molecules at sites the original structure depositions omitted. We use this improved model to address the data constraint. By mapping the Pareto front of achievable accuracy of water molecule prediction, alongside analysis of different training data schemas, we quantified the trade-off between data quantity and data quality, demonstrating that the diversity of high-quality structures is limiting the possible results. Overall, WaterFlow predicts ordered water to serve as a solvent module for structure-based drug design and for water molecule placement during crystallographic refinement.