The R-factor gap in macromolecular crystallography: an untapped potential for insights on accurate structures

Uncovering Protein Ensembles: Automated Multiconformer Model Building for X-ray Crystallography and Cryo-EM

In their folded state, biomolecules exchange between multiple conformational states that are crucial for their function. Traditional structural biology methods, such as X-ray crystallography and cryogenic electron microscopy (cryo-EM), produce density maps that are ensemble averages, reflecting molecules in various conformations. Yet, most models derived from these maps explicitly represent only a single conformation, overlooking the complexity of biomolecular structures. To accurately reflect the diversity of biomolecular forms, there is a pressing need to shift towards modeling structural ensembles that mirror the experimental data. However, the challenge of distinguishing signal from noise complicates manual efforts to create these models. In response, we introduce the latest enhancements to qFit, an automated computational strategy designed to incorporate protein conformational heterogeneity into models built into density maps. These algorithmic improvements in qFit are substantiated by superior R f r e e and geometry metrics across a wide range of proteins. Importantly, unlike more complex multicopy ensemble models, the multiconformer models produced by qFit can be manually modified in most major model building software (e.g. Coot) and fit can be further improved by refinement using standard pipelines (e.g. Phenix, Refmac, Buster). By reducing the barrier of creating multiconformer models, qFit can foster the development of new hypotheses about the relationship between macromolecular conformational dynamics and function.

Application of density matrix Wigner transforms for ultrafast macromolecular and chemical x-ray crystallography

Time-Resolved Serial Femtosecond Crystallography (TR-SFX) conducted at X-ray Free Electron Lasers (XFELs) has become a powerful tool for capturing macromolecular structural movies of light-initiated processes. As the capabilities of XFELs advance, we anticipate that a new range of coherent control and structural Raman measurements will become achievable. Shorter optical and x-ray pulse durations and increasingly more exotic pulse regimes are becoming available at free electron lasers. Moreover, with high repetition enabled by the superconducting technology of European XFEL (EuXFEL) and Linac Coherent Light Source (LCLS-II) , it will be possible to improve the signal-to-noise ratio of the light-induced differences, allowing for the observation of vibronic motion on the sub-Angstrom level. To predict and assign this coherent motion, which is measurable with a structural technique, new theoretical approaches must be developed. In this paper, we present a theoretical density matrix approach to model the various population and coherent dynamics of a system, which considers molecular system parameters and excitation conditions. We emphasize the use of the Wigner transform of the time-dependent density matrix, which provides a phase space representation that can be directly compared to the experimental positional displacements measured in a TR-SFX experiment. Here, we extend the results from simple models to include more realistic schemes that include large relaxation terms. We explore a variety of pulse schemes using multiple model systems using realistic parameters. An open-source software package is provided to perform the density matrix simulation and Wigner transformations. The open-source software allows us to define any arbitrary level schemes as well as any arbitrary electric field in the interaction Hamiltonian.

Structure is beauty, but not always truth

Structural biology, as powerful as it is, can be misleading. We highlight four fundamental challenges: interpreting raw experimental data; accounting for motion; addressing the misleading nature of in vitro structures; and unraveling interactions between drugs and "anti-targets." Overcoming these challenges will amplify the impact of structural biology on drug discovery.

Deep residual networks for crystallography trained on synthetic data

The use of artificial intelligence to process diffraction images is challenged by the need to assemble large and precisely designed training data sets. To address this, a codebase called Resonet was developed for synthesizing diffraction data and training residual neural networks on these data. Here, two per-pattern capabilities of Resonet are demonstrated: (i) interpretation of crystal resolution and (ii) identification of overlapping lattices. Resonet was tested across a compilation of diffraction images from synchrotron experiments and X-ray free-electron laser experiments. Crucially, these models readily execute on graphics processing units and can thus significantly outperform conventional algorithms. While Resonet is currently utilized to provide real-time feedback for macromolecular crystallography users at the Stanford Synchrotron Radiation Lightsource, its simple Python-based interface makes it easy to embed in other processing frameworks. This work highlights the utility of physics-based simulation for training deep neural networks and lays the groundwork for the development of additional models to enhance diffraction collection and analysis.

Interpreting macromolecular diffraction through simulation

This chapter discusses the use of diffraction simulators to improve experimental outcomes in macromolecular crystallography, in particular for future experiments aimed at diffuse scattering. Consequential decisions for upcoming data collection include the selection of either a synchrotron or free electron laser X-ray source, rotation geometry or serial crystallography, and fiber-coupled area detector technology vs. pixel-array detectors. The hope is that simulators will provide insights to make these choices with greater confidence. Simulation software, especially those packages focused on physics-based calculation of the diffraction, can help to predict the location, size, shape, and profile of Bragg spots and diffuse patterns in terms of an underlying physical model, including assumptions about the crystal's mosaic structure, and therefore can point to potential issues with data analysis in the early planning stages. Also, once the data are collected, simulation may offer a pathway to improve the measurement of diffraction, especially with weak data, and might help to treat problematic cases such as overlapping patterns.

MD simulations of macromolecular crystals: Implications for the analysis of Bragg and diffuse scattering

Some of our most detailed information about structure and dynamics of macromolecules comes from X-ray-diffraction studies in crystalline environments. More than 170,000 atomic models have been deposited in the Protein Data Bank, and the number of observations (typically of intensities of Bragg diffraction peaks) is generally quite large, when compared to other experimental methods. Nevertheless, the general agreement between calculated and observed intensities is far outside the experimental precision, and the majority of scattered photons fall between the sharp Bragg peaks, and are rarely taken into account. This chapter considers how molecular dynamics simulations can be used to explore the connections between microscopic behavior in a crystalline lattice and observed scattering intensities, and point the way to new atomic models that could more faithfully recapitulate Bragg intensities and extract useful information from the diffuse scattering that lies between those peaks.

Progress in detection of and correction for low-energy contamination

Contamination with low-energy radiation leads to an increased number of weighted residuals being larger in absolute terms than three standard uncertainties. For a Gaussian distribution, these rare events occur only in 0.27% of all cases, which is a small number for small- to medium-sized data sets. The correct detection of rare events - and an adequate correction procedure - thus relies crucially on correct standard uncertainties, which are often not available [Henn (2019), Crystallogr. Rev. 25, 83-156]. It is therefore advisable to use additional, more robust, metrics to complement the established ones. These metrics are developed here and applied to reference data sets from two different publications about low-energy contamination. Other systematic errors were found in the reference data sets. These errors compromise the correction procedures and may lead to under- or overcompensation. This can be demonstrated clearly with the new metrics. Empirical correction procedures generally may be compromised or bound to fail in the presence of other systematic errors. The following systematic errors, which were found in the reference data sets, need to be corrected for prior to application of the low-energy contamination correction procedure: signals of 2λ contamination, extinction, disorder, twinning, and too-large or too-low standard uncertainties (this list may not be complete). All five reference data sets of one publication show a common resolution-dependent systematic error of unknown origin. How this affects the correction procedure can be stated only after elimination of this error. The methodological improvements are verified with data published by other authors.

Water Networks in Photosystem II Using Crystalline Molecular Dynamics Simulations and Room-Temperature XFEL Serial Crystallography

Structural dynamics of water and its hydrogen-bonding networks play an important role in enzyme function via the transport of protons, ions, and substrates. To gain insights into these mechanisms in the water oxidation reaction in Photosystem II (PS II), we have performed crystalline molecular dynamics (MD) simulations of the dark-stable S1 state. Our MD model consists of a full unit cell with 8 PS II monomers in explicit solvent (861 894 atoms), enabling us to compute the simulated crystalline electron density and to compare it directly with the experimental density from serial femtosecond X-ray crystallography under physiological temperature collected at X-ray free electron lasers (XFELs). The MD density reproduced the experimental density and water positions with high fidelity. The detailed dynamics in the simulations provided insights into the mobility of water molecules in the channels beyond what can be interpreted from experimental B-factors and electron densities alone. In particular, the simulations revealed fast, coordinated exchange of waters at sites where the density is strong, and water transport across the bottleneck region of the channels where the density is weak. By computing MD hydrogen and oxygen maps separately, we developed a novel Map-based Acceptor-Donor Identification (MADI) technique that yields information which helps to infer hydrogen-bond directionality and strength. The MADI analysis revealed a series of hydrogen-bond wires emanating from the Mn cluster through the Cl1 and O4 channels; such wires might provide pathways for proton transfer during the reaction cycle of PS II. Our simulations provide an atomistic picture of the dynamics of water and hydrogen-bonding networks in PS II, with implications for the specific role of each channel in the water oxidation reaction.

Self-Supervised Deep Learning for Model Correction in the Computational Crystallography Toolbox

The Computational Crystallography Toolbox (cctbx) is open-source software that allows for processing of crystallographic data, including from serial femtosecond crystallography (SFX), for macromolecular structure determination. We aim to use the modules in cctbx to determine the oxidation state of individual metal atoms in a macromolecule. Changes in oxidation state are reflected in small shifts of the atom's X-ray absorption edge. These energy shifts can be extracted from the diffraction images recorded in serial femtosecond crystallography, given knowledge of a forward physics model. However, as the diffraction changes only slightly due to the absorption edge shift, inaccuracies in the forward physics model make it extremely challenging to observe the oxidation state. In this work, we describe the potential impact of using self-supervised deep learning to correct the scientific model in cctbx and provide uncertainty quantification. We provide code for forward model simulation and data analysis, built from cctbx modules, at https://github.com/gigantocypris/SPREAD, which can be integrated with machine learning. We describe open questions in algorithm development to help spur advances through dialog between crystallographers and machine learning researchers. New methods could help elucidate charge transfer processes in many reactions, including key events in photosynthesis.

Emerging Time-Resolved X-Ray Diffraction Approaches for Protein Dynamics

Proteins guide the flows of information, energy, and matter that make life possible by accelerating transport and chemical reactions, by allosterically modulating these reactions, and by forming dynamic supramolecular assemblies. In these roles, conformational change underlies functional transitions. Time-resolved X-ray diffraction methods characterize these transitions either by directly triggering sequences of functionally important motions or, more broadly, by capturing the motions of which proteins are capable. To date, most successful have been experiments in which conformational change is triggered in light-dependent proteins. In this review, I emphasize emerging techniques that probe the dynamic basis of function in proteins lacking natively light-dependent transitions and speculate about extensions and further possibilities. In addition, I review how the weaker and more distributed signals in these data push the limits of the capabilities of analytical methods. Taken together, these new methods are beginning to establish a powerful paradigm for the study of the physics of protein function.

Robust total X-ray scattering workflow to study correlated motion of proteins in crystals

The breathing motions of proteins are thought to play a critical role in function. However, current techniques to study key collective motions are limited to spectroscopy and computation. We present a high-resolution experimental approach based on the total scattering from protein crystals at room temperature (TS/RT-MX) that captures both structure and collective motions. To reveal the scattering signal from protein motions, we present a general workflow that enables robust subtraction of lattice disorder. The workflow introduces two methods: GOODVIBES, a detailed and refinable lattice disorder model based on the rigid-body vibrations of a crystalline elastic network; and DISCOBALL, an independent method of validation that estimates the displacement covariance between proteins in the lattice in real space. Here, we demonstrate the robustness of this workflow and further demonstrate how it can be interfaced with MD simulations towards obtaining high-resolution insight into functionally important protein motions.

Crystal diffraction prediction and partiality estimation using Gaussian basis functions

The recent diversification of macromolecular crystallographic experiments including the use of pink beams, convergent electron diffraction and serial snapshot crystallography has shown the limitations of using the Laue equations for diffraction prediction. This article gives a computationally efficient way of calculating approximate crystal diffraction patterns given varying distributions of the incoming beam, crystal shapes and other potentially hidden parameters. This approach models each pixel of a diffraction pattern and improves data processing of integrated peak intensities by enabling the correction of partially recorded reflections. The fundamental idea is to express the distributions as weighted sums of Gaussian functions. The approach is demonstrated on serial femtosecond crystallography data sets, showing a significant decrease in the required number of patterns to refine a structure to a given error.

Molecular-dynamics simulation methods for macromolecular crystallography

It is investigated whether molecular-dynamics (MD) simulations can be used to enhance macromolecular crystallography (MX) studies. Historically, protein crystal structures have been described using a single set of atomic coordinates. Because conformational variation is important for protein function, researchers now often build models that contain multiple structures. Methods for building such models can fail, however, in regions where the crystallographic density is difficult to interpret, for example at the protein-solvent interface. To address this limitation, a set of MD-MX methods that combine MD simulations of protein crystals with conventional modeling and refinement tools have been developed. In an application to a cyclic adenosine monophosphate-dependent protein kinase at room temperature, the procedure improved the interpretation of ambiguous density, yielding an alternative water model and a revised protein model including multiple conformations. The revised model provides mechanistic insights into the catalytic and regulatory interactions of the enzyme. The same methods may be used in other MX studies to seek mechanistic insights.

Simulating diffraction photographs based on molecular dynamics trajectories of a protein crystal: a new option to examine structure-solving strategies in protein crystallography

A molecular dynamics (MD)-based pipeline has been designed and implemented to emulate the entire process of collecting diffraction photographs and calculating crystallographic structures of proteins from them. Using a structure of lysozyme solved in-house, a supercell comprising 125 (5 × 5 × 5) crystal unit cells containing a total of 1000 protein molecules and explicit interstitial solvent was constructed. For this system, two 300 ns MD trajectories at 298 and 250 K were recorded. A series of snapshots from these trajectories were then used to simulate a fully realistic set of diffraction photographs, which were further fed into the standard pipeline for structure determination. The resulting structures show very good agreement with the underlying MD model not only in terms of coordinates but also in terms of B factors; they are also consistent with the original experimental structure. The developed methodology should find a range of applications, such as optimizing refinement protocols to solve crystal structures and extracting dynamics information from diffraction data or diffuse scattering.

Combined Structural Analysis and Molecular Dynamics Reveal Penicillin-Binding Protein Inhibition Mode with β-Lactones

β-Lactam antibiotics comprise one of the most widely used therapeutic classes to combat bacterial infections. This general scaffold has long been known to inhibit bacterial cell wall biosynthesis by inactivating penicillin-binding proteins (PBPs); however, bacterial resistance to β-lactams is now widespread, and new strategies are urgently needed to target PBPs and other proteins involved in bacterial cell wall formation. A key requirement in the identification of strategies to overcome resistance is a deeper understanding of the roles of the PBPs and their associated proteins during cell growth and division, such as can be obtained with the use of selective chemical probes. Probe development has typically depended upon known PBP inhibitors, which have historically been thought to require a negatively charged moiety that mimics the C-terminus of the PBP natural peptidoglycan substrate, d-Ala-d-Ala. However, we have identified a new class of β-lactone-containing molecules that interact with PBPs, often in an isoform-specific manner, and do not incorporate this C-terminal mimetic. Here, we report a series of structural biology experiments and molecular dynamics simulations that we utilized to evaluate specific binding modes of this novel PBP inhibitor class. In this work, we obtained <2 Å resolution X-ray structures of four β-lactone probes bound to PBP1b from Streptococcus pneumoniae. Despite their diverging recognition modes beyond the site of covalent modification, these four probes all efficiently labeled PBP1b, as well as other PBPs from S. pneumoniae. From these structures, we analyzed protein-ligand interactions and characterized the β-lactone-bound active sites using in silico mutagenesis and molecular dynamics. Our approach has clarified the dynamic interaction profile in this series of ligands, expanding the understanding of PBP inhibitor binding.

Medical contrast agents as promising tools for biomacromolecular SAXS experiments

Lanthanide-based complexes are presented as a promising class of molecules for efficient SAXS contrast-variation experiments. Their interactions and contrast properties are analyzed for an oligomeric protein and a protein–RNA complex.

Small-angle X-ray scattering (SAXS) has become an indispensable tool in structural biology, complementing atomic-resolution techniques. It is sensitive to the electron-density difference between solubilized biomacromolecules and the buffer, and provides information on molecular masses, particle dimensions and interactions, low-resolution conformations and pair distance-distribution functions. When SAXS data are recorded at multiple contrasts, i.e. at different solvent electron densities, it is possible to probe, in addition to their overall shape, the internal electron-density profile of biomacromolecular assemblies. Unfortunately, contrast-variation SAXS has been limited by the range of solvent electron densities attainable using conventional co-solutes (for example sugars, glycerol and salt) and by the fact that some biological systems are destabilized in their presence. Here, SAXS contrast data from an oligomeric protein and a protein–RNA complex are presented in the presence of iohexol and Gd-HPDO3A, two electron-rich molecules that are used in biomedical imaging and that belong to the families of iodinated and lanthanide-based complexes, respectively. Moderate concentrations of both molecules allowed solvent electron densities matching those of proteins to be attained. While iohexol yielded higher solvent electron densities (per mole), it interacted specifically with the oligomeric protein and precipitated the protein–RNA complex. Gd-HPDO3A, while less efficient (per mole), did not disrupt the structural integrity of either system, and atomic models could be compared with the SAXS data. Due to their elevated solubility and electron density, their chemical inertness, as well as the possibility of altering their physico-chemical properties, lanthanide-based complexes represent a class of molecules with promising potential for contrast-variation SAXS experiments on diverse biomacromolecular systems.

Integral equation models for solvent in macromolecular crystals

The solvent can occupy up to ∼70% of macromolecular crystals, and hence, having models that predict solvent distributions in periodic systems could improve the interpretation of crystallographic data. Yet, there are few implicit solvent models applicable to periodic solutes, and crystallographic structures are commonly solved assuming a flat solvent model. Here, we present a newly developed periodic version of the 3D-reference interaction site model (RISM) integral equation method that is able to solve efficiently and describe accurately water and ion distributions in periodic systems; the code can compute accurate gradients that can be used in minimizations or molecular dynamics simulations. The new method includes an extension of the Ornstein–Zernike equation needed to yield charge neutrality for charged solutes, which requires an additional contribution to the excess chemical potential that has not been previously identified; this is an important consideration for nucleic acids or any other charged system where most or all the counter- and co-ions are part of the “disordered” solvent. We present several calculations of proteins, RNAs, and small molecule crystals to show that x-ray scattering intensities and the solvent structure predicted by the periodic 3D-RISM solvent model are in closer agreement with the experiment than are intensities computed using the default flat solvent model in the refmac5 or phenix refinement programs, with the greatest improvement in the 2 to 4 Å range. Prospects for incorporating integral equation models into crystallographic refinement are discussed.

Universally Accessible Structural Data on Macromolecular Conformation, Assembly, and Dynamics by Small Angle X-Ray Scattering for DNA Repair Insights

Structures provide a critical breakthrough step for biological analyses, and small angle X-ray scattering (SAXS) is a powerful structural technique to study dynamic DNA repair proteins. As toxic and mutagenic repair intermediates need to be prevented from inadvertently harming the cell, DNA repair proteins often chaperone these intermediates through dynamic conformations, coordinated assemblies, and allosteric regulation. By measuring structural conformations in solution for both proteins, DNA, RNA, and their complexes, SAXS provides insight into initial DNA damage recognition, mechanisms for validation of their substrate, and pathway regulation. Here, we describe exemplary SAXS analyses of a DNA damage response protein spanning from what can be derived directly from the data to obtaining super resolution through the use of SAXS selection of atomic models. We outline strategies and tactics for practical SAXS data collection and analysis. Making these structural experiments in reach of any basic and clinical researchers who have protein, SAXS data can readily be collected at government-funded synchrotrons, typically at no cost for academic researchers. In addition to discussing how SAXS complements and enhances cryo-electron microscopy, X-ray crystallography, NMR, and computational modeling, we furthermore discuss taking advantage of recent advances in protein structure prediction in combination with SAXS analysis.

Target highlights in CASP14: Analysis of models by structure providers

The biological and functional significance of selected Critical Assessment of Techniques for Protein Structure Prediction 14 (CASP14) targets are described by the authors of the structures. The authors highlight the most relevant features of the target proteins and discuss how well these features were reproduced in the respective submitted predictions. The overall ability to predict three-dimensional structures of proteins has improved remarkably in CASP14, and many difficult targets were modeled with impressive accuracy. For the first time in the history of CASP, the experimentalists not only highlighted that computational models can accurately reproduce the most critical structural features observed in their targets, but also envisaged that models could serve as a guidance for further studies of biologically-relevant properties of proteins.

Detecting anomalies in X-ray diffraction images using convolutional neural networks

Our understanding of life is based upon the interpretation of macromolecular structures and their dynamics. Almost 90% of currently known macromolecular models originated from electron density maps constructed using X-ray diffraction images. Even though diffraction images are critical for structure determination, due to their vast amounts and noisy, non-intuitive nature, their quality is rarely inspected. In this paper, we use recent advances in machine learning to automatically detect seven types of anomalies in X-ray diffraction images. For this purpose, we utilize a novel X-ray beam center detection algorithm, propose three different image representations, and compare the predictive performance of general-purpose classifiers and deep convolutional neural networks (CNNs). In benchmark tests on a set of 6,311 X-ray diffraction images, the proposed CNN achieved between 87% and 99% accuracy depending on the type of anomaly. Experimental results show that the proposed anomaly detection system can be considered suitable for early detection of sub-optimal data collection conditions and malfunctions at X-ray experimental stations.

MagC is a NplC/P60-like member of the α-2-macroglobulin Mag complex of Pseudomonas aeruginosa that interacts with peptidoglycan

Bacterial α-2 macroglobulins (A2Ms) structurally resemble the large spectrum protease inhibitors of the eukaryotic immune system. In Pseudomonas aeruginosa, MagD acts as an A2M and is expressed within a six-gene operon encoding the MagA-F proteins. In this work, we employ isothermal calorimetry (ITC), analytical ultracentrifugation (AUC), and X-ray crystallography to investigate the function of MagC and show that MagC associates with the macroglobulin complex and with the peptidoglycan (PG). However, the catalytic residues of MagC display an inactive conformation that could suggest that it binds to PG but does not degrade it. We hypothesize that MagC could serve as an anchor between the MagD macroglobulin and the PG and could provide stabilization and/or regulation for the entire complex.

X-ray scattering reveals disordered linkers and dynamic interfaces in complexes and mechanisms for DNA double-strand break repair impacting cell and cancer biology

Evolutionary selection ensures specificity and efficiency in dynamic metastable macromolecular machines that repair DNA damage without releasing toxic and mutagenic intermediates. Here we examine non‐homologous end joining (NHEJ) as the primary conserved DNA double‐strand break (DSB) repair process in human cells. NHEJ has exemplary key roles in networks determining the development, outcome of cancer treatments by DSB‐inducing agents, generation of antibody and T‐cell receptor diversity, and innate immune response for RNA viruses. We determine mechanistic insights into NHEJ structural biochemistry focusing upon advanced small angle X‐ray scattering (SAXS) results combined with X‐ray crystallography (MX) and cryo‐electron microscopy (cryo‐EM). SAXS coupled to atomic structures enables integrated structural biology for objective quantitative assessment of conformational ensembles and assemblies in solution, intra‐molecular distances, structural similarity, functional disorder, conformational switching, and flexibility. Importantly, NHEJ complexes in solution undergo larger allosteric transitions than seen in their cryo‐EM or MX structures. In the long‐range synaptic complex, X‐ray repair cross‐complementing 4 (XRCC4) plus XRCC4‐like‐factor (XLF) form a flexible bridge and linchpin for DNA ends bound to KU heterodimer (Ku70/80) and DNA‐PKcs (DNA‐dependent protein kinase catalytic subunit). Upon binding two DNA ends, auto‐phosphorylation opens DNA‐PKcs dimer licensing NHEJ via concerted conformational transformations of XLF‐XRCC4, XLF–Ku80, and LigIV^BRCT–Ku70 interfaces. Integrated structures reveal multifunctional roles for disordered linkers and modular dynamic interfaces promoting DSB end processing and alignment into the short‐range complex for ligation by LigIV. Integrated findings define dynamic assemblies fundamental to designing separation‐of‐function mutants and allosteric inhibitors targeting conformational transitions in multifunctional complexes.

Vagabond: bond-based parametrization reduces overfitting for refinement of proteins

Structural biology methods have delivered over 150 000 high-resolution structures of macromolecules, which have fundamentally altered our understanding of biology and our approach to developing new medicines. However, the description of molecular flexibility is instrinsically flawed and in almost all cases, regardless of the experimental method used for structure determination, there remains a strong overfitting bias during molecular model building and refinement. In the worst case this can lead to wholly incorrect structures and thus incorrect biological interpretations. Here, by reparametrizing the description of these complex structures in terms of bonds rather than atomic positions, and by modelling flexibility using a deterministic ensemble of structures, it is demonstrated that structures can be described using fewer parameters than in conventional refinement. The current implementation, applied to X-ray diffraction data, significantly reduces the extent of overfitting, allowing the experimental data to reveal more biological information in electron-density maps.

The PDB and the ribosome

This essay, which was written to commemorate the 50th anniversary of the Protein Data Bank, opens with some comments about the intentions of the scientists who pressed for its establishment and the nature of services it provides. It includes a brief account of the events that resulted in the determination of the crystal structure of the large ribosomal subunit from Haloarcula marismortui. The magnitude of the challenge the first ribosome crystal structures posed for the PDB is commented upon, and in the description of subsequent developments in the ribosome structure field that follows, it is pointed out that cryo-EM has replaced X-ray crystallography as the method of choice for investigating ribosome structure.

FMX - the Frontier Microfocusing Macromolecular Crystallography Beamline at the National Synchrotron Light Source II

Two new macromolecular crystallography (MX) beamlines at the National Synchrotron Light Source II, FMX and AMX, opened for general user operation in February 2017 [Schneider et al. (2013). J. Phys. Conf. Ser. 425, 012003; Fuchs et al. (2014). J. Phys. Conf. Ser. 493, 012021; Fuchs et al. (2016). AIP Conf. Proc. SRI2015, 1741, 030006]. FMX, the micro-focusing Frontier MX beamline in sector 17-ID-2 at NSLS-II, covers a 5-30 keV photon energy range and delivers a flux of 4.0 × 1012 photons s-1 at 1 Å into a 1 µm × 1.5 µm to 10 µm × 10 µm (V × H) variable focus, expected to reach 5 × 1012 photons s-1 at final storage-ring current. This flux density surpasses most MX beamlines by nearly two orders of magnitude. The high brightness and microbeam capability of FMX are focused on solving difficult crystallographic challenges. The beamline's flexible design supports a wide range of structure determination methods - serial crystallography on micrometre-sized crystals, raster optimization of diffraction from inhomogeneous crystals, high-resolution data collection from large-unit-cell crystals, room-temperature data collection for crystals that are difficult to freeze and for studying conformational dynamics, and fully automated data collection for sample-screening and ligand-binding studies. FMX's high dose rate reduces data collection times for applications like serial crystallography to minutes rather than hours. With associated sample lifetimes as short as a few milliseconds, new rapid sample-delivery methods have been implemented, such as an ultra-high-speed high-precision piezo scanner goniometer [Gao et al. (2018). J. Synchrotron Rad. 25, 1362-1370], new microcrystal-optimized micromesh well sample holders [Guo et al. (2018). IUCrJ, 5, 238-246] and highly viscous media injectors [Weierstall et al. (2014). Nat. Commun. 5, 3309]. The new beamline pushes the frontier of synchrotron crystallography and enables users to determine structures from difficult-to-crystallize targets like membrane proteins, using previously intractable crystals of a few micrometres in size, and to obtain quality structures from irregular larger crystals.

qFit 3: Protein and ligand multiconformer modeling for X-ray crystallographic and single-particle cryo-EM density maps

New X-ray crystallography and cryo-electron microscopy (cryo-EM) approaches yield vast amounts of structural data from dynamic proteins and their complexes. Modeling the full conformational ensemble can provide important biological insights, but identifying and modeling an internally consistent set of alternate conformations remains a formidable challenge. qFit efficiently automates this process by generating a parsimonious multiconformer model. We refactored qFit from a distributed application into software that runs efficiently on a small server, desktop, or laptop. We describe the new qFit 3 software and provide some examples. qFit 3 is open-source under the MIT license, and is available at https://github.com/ExcitedStates/qfit-3.0.

Beyond integration: modeling every pixel to obtain better structure factors from stills

Most crystallographic data processing methods use pixel integration. In serial femtosecond crystallography (SFX), the intricate interaction between the reciprocal lattice point and the Ewald sphere is integrated out by averaging symmetrically equivalent observations recorded across a large number (104-106) of exposures. Although sufficient for generating biological insights, this approach converges slowly, and using it to accurately measure anomalous differences has proved difficult. This report presents a novel approach for increasing the accuracy of structure factors obtained from SFX data. A physical model describing all observed pixels is defined to a degree of complexity such that it can decouple the various contributions to the pixel intensities. Model dependencies include lattice orientation, unit-cell dimensions, mosaic structure, incident photon spectra and structure factor amplitudes. Maximum likelihood estimation is used to optimize all model parameters. The application of prior knowledge that structure factor amplitudes are positive quantities is included in the form of a reparameterization. The method is tested using a synthesized SFX dataset of ytterbium(III) lysozyme, where each X-ray laser pulse energy is centered at 9034 eV. This energy is 100 eV above the Yb3+ L-III absorption edge, so the anomalous difference signal is stable at 10 electrons despite the inherent energy jitter of each femtosecond X-ray laser pulse. This work demonstrates that this approach allows the determination of anomalous structure factors with very high accuracy while requiring an order-of-magnitude fewer shots than conventional integration-based methods would require to achieve similar results.

Paired refinement under the control of PAIREF

Crystallographic resolution is a key characteristic of diffraction data and represents one of the first decisions an experimenter has to make in data evaluation. Conservative approaches to the high-resolution cutoff determination are based on a number of criteria applied to the processed X-ray diffraction data only. However, high-resolution data that are weaker than arbitrary cutoffs can still result in the improvement of electron-density maps and refined structure models. Therefore, the impact of reflections from resolution shells higher than those previously used in conservative structure refinement should be analysed by the paired refinement protocol. For this purpose, a tool called PAIREF was developed to provide automation of this protocol. As a new feature, a complete cross-validation procedure has also been implemented. Here, the design, usage and control of the program are described, and its application is demonstrated on six data sets. The results prove that the inclusion of high-resolution data beyond the conventional criteria can lead to more accurate structure models.

Making a difference in multi-data-set crystallography: simple and deterministic data-scaling/selection methods

Phasing by single-wavelength anomalous diffraction (SAD) from multiple crystallographic data sets can be particularly demanding because of the weak anomalous signal and possible non-isomorphism. The identification and exclusion of non-isomorphous data sets by suitable indicators is therefore indispensable. Here, simple and robust data-selection methods are described. A multi-dimensional scaling procedure is first used to identify data sets with large non-isomorphism relative to clusters of other data sets. Within each cluster that it identifies, further selection is based on the weighted ΔCC1/2, a quantity representing the influence of a set of reflections on the overall CC1/2 of the merged data. The anomalous signal is further improved by optimizing the scaling protocol. The success of iterating the selection and scaling steps was verified by substructure determination and subsequent structure solution. Three serial synchrotron crystallography (SSX) SAD test cases with hundreds of partial data sets and one test case with 62 complete data sets were analyzed. Structure solution was dramatically simplified with this procedure, and enabled solution of the structures after a few selection/scaling iterations. To explore the limits, the procedure was tested with much fewer data than originally required and could still solve the structure in several cases. In addition, an SSX data challenge, minimizing the number of (simulated) data sets necessary to solve the structure, was significantly underbid.

Towards the spatial resolution of metalloprotein charge states by detailed modeling of XFEL crystallographic diffraction

Oxidation states of individual metal atoms within a metalloprotein can be assigned by examining X-ray absorption edges, which shift to higher energy for progressively more positive valence numbers. Indeed, X-ray crystallography is well suited for such a measurement, owing to its ability to spatially resolve the scattering contributions of individual metal atoms that have distinct electronic environments contributing to protein function. However, as the magnitude of the shift is quite small, about +2 eV per valence state for iron, it has only been possible to measure the effect when performed with monochromated X-ray sources at synchrotron facilities with energy resolutions in the range 2-3 × 10-4 (ΔE/E). This paper tests whether X-ray free-electron laser (XFEL) pulses, which have a broader bandpass (ΔE/E = 3 × 10-3) when used without a monochromator, might also be useful for such studies. The program nanoBragg is used to simulate serial femtosecond crystallography (SFX) diffraction images with sufficient granularity to model the XFEL spectrum, the crystal mosaicity and the wavelength-dependent anomalous scattering factors contributed by two differently charged iron centers in the 110-amino-acid protein, ferredoxin. Bayesian methods are then used to deduce, from the simulated data, the most likely X-ray absorption curves for each metal atom in the protein, which agree well with the curves chosen for the simulation. The data analysis relies critically on the ability to measure the incident spectrum for each pulse, and also on the nanoBragg simulator to predict the size, shape and intensity profile of Bragg spots based on an underlying physical model that includes the absorption curves, which are then modified to produce the best agreement with the simulated data. This inference methodology potentially enables the use of SFX diffraction for the study of metalloenzyme mechanisms and, in general, offers a more detailed approach to Bragg spot data reduction.

Factors influencing estimates of coordinate error for molecular replacement

Good prior estimates of the effective root-mean-square deviation (r.m.s.d.) between the atomic coordinates of the model and the target optimize the signal in molecular replacement, thereby increasing the success rate in difficult cases. Previous studies using protein structures solved by X-ray crystallography as models showed that optimal error estimates (refined after structure solution) were correlated with the sequence identity between the model and target, and with the number of residues in the model. Here, this work has been extended to find additional correlations between parameters of the model and the target and hence improved prior estimates of the coordinate error. Using a graph database, a curated set of 6030 molecular-replacement calculations using models that had been solved by X-ray crystallography was analysed to consider about 120 model and target parameters. Improved estimates were achieved by replacing the sequence identity with the Gonnet score for sequence similarity, as well as by considering the resolution of the target structure and the MolProbity score of the model. This approach was extended by analysing 12 610 additional molecular-replacement calculations where the model was determined by NMR. The median r.m.s.d. between pairs of models in an ensemble was found to be correlated with the estimated r.m.s.d. to the target. For models solved by NMR, the overall coordinate error estimates were larger than for structures determined by X-ray crystallography, and were more highly correlated with the number of residues.

JUNGFRAU detector for brighter x-ray sources: Solutions for IT and data science challenges in macromolecular crystallography

In this paper, we present a data workflow developed to operate the adJUstiNg Gain detector FoR the Aramis User station (JUNGFRAU) adaptive gain charge integrating pixel-array detectors at macromolecular crystallography beamlines. We summarize current achievements for operating at 9 GB/s data-rate a JUNGFRAU with 4 Mpixel at 1.1 kHz frame-rate and preparations to operate at 46 GB/s data-rate a JUNGFRAU with 10 Mpixel at 2.2 kHz in the future. In this context, we highlight the challenges for computer architecture and how these challenges can be addressed with innovative hardware including IBM POWER9 servers and field-programmable gate arrays. We discuss also data science challenges, showing the effect of rounding and lossy compression schemes on the MX JUNGFRAU detector images.

Controlled dehydration, structural flexibility and gadolinium MRI contrast compound binding in the human plasma glycoprotein afamin

Afamin, which is a human blood plasma glycoprotein, a putative multifunctional transporter of hydrophobic molecules and a marker for metabolic syndrome, poses multiple challenges for crystallographic structure determination, both practically and in analysis of the models. Several hundred crystals were analysed, and an unusual variability in cell volume and difficulty in solving the structure despite an ∼34% sequence identity with nonglycosylated human serum albumin indicated that the molecule exhibits variable and context-sensitive packing, despite the simplified glycosylation in insect cell-expressed recombinant afamin. Controlled dehydration of the crystals was able to stabilize the orthorhombic crystal form, reducing the number of molecules in the asymmetric unit from the monoclinic form and changing the conformational state of the protein. An iterative strategy using fully automatic experiments available on MASSIF-1 was used to quickly determine the optimal protocol to achieve the phase transition, which should be readily applicable to many types of sample. The study also highlights the drawback of using a single crystallographic structure model for computational modelling purposes given that the conformational state of the binding sites and the electron density in the binding site, which is likely to result from PEGs, greatly varies between models. This also holds for the analysis of nonspecific low-affinity ligands, where often a variety of fragments with similar uncertainty can be modelled, inviting interpretative bias. As a promiscuous transporter, afamin also seems to bind gadoteridol, a magnetic resonance imaging contrast compound, in at least two sites. One pair of gadoteridol molecules is located near the human albumin Sudlow site, and a second gadoteridol molecule is located at an intermolecular site in proximity to domain IA. The data from the co-crystals support modern metrics of data quality in the context of the information that can be gleaned from data sets that would be abandoned on classical measures.

The binding of precipitant ions in the tetragonal crystals of hen egg white lysozyme

The bonds between lysozyme molecules and precipitant ions in single crystals grown with chlorides of several metals are analysed on the basis of crystal structure data. Crystals of tetragonal hen egg lysozyme (HEWL) were grown with chlorides of several alkali and transition metals (LiCl, NaCl, KCl, NiCl₂ and CuCl₂) as precipitants and the three-dimensional structures were determined at 1.35 Å resolution by X-ray diffraction method. The positions of metal and chloride ions attached to the protein were located, divided into three groups and analysed. Some of them, in accordance with the recently proposed and experimentally confirmed crystal growth model, provide connections in protein dimers and octamers that are precursor clusters in the crystallization lysozyme solution. The first group, including Cu⁺², Ni⁺² and Na⁺¹ cations, binds specifically to the protein molecule. The second group consists of metal and chloride ions bound inside the dimers and octamers. The third group of ions can participate in connections between the octamers that are suggested as building units during the crystal growth. The arrangement of chloride and metal ions associated with lysozyme molecule at all stages of the crystallization solution formation and crystal growth is discussed.Communicated by Ramaswamy H. Sarma.

Liquid-like and rigid-body motions in molecular-dynamics simulations of a crystalline protein

To gain insight into crystalline protein dynamics, we performed molecular-dynamics (MD) simulations of a periodic 2 × 2 × 2 supercell of staphylococcal nuclease. We used the resulting MD trajectories to simulate X-ray diffraction and to study collective motions. The agreement of simulated X-ray diffraction with the data is comparable to previous MD simulation studies. We studied collective motions by analyzing statistically the covariance of alpha-carbon position displacements. The covariance decreases exponentially with the distance between atoms, which is consistent with a liquidlike motions (LLM) model, in which the protein behaves like a soft material. To gain finer insight into the collective motions, we examined the covariance behavior within a protein molecule (intraprotein) and between different protein molecules (interprotein). The interprotein atom pairs, which dominate the overall statistics, exhibit LLM behavior; however, the intraprotein pairs exhibit behavior that is consistent with a superposition of LLM and rigid-body motions (RBM). Our results indicate that LLM behavior of global dynamics is present in MD simulations of a protein crystal. They also show that RBM behavior is detectable in the simulations but that it is subsumed by the LLM behavior. Finally, the results provide clues about how correlated motions of atom pairs both within and across proteins might manifest in diffraction data. Overall, our findings increase our understanding of the connection between molecular motions and diffraction data and therefore advance efforts to extract information about functionally important motions from crystallography experiments.

Resolution and dose dependence of radiation damage in biomolecular systems

The local Fourier-space relation between diffracted intensity I, diffraction wavevector q and dose D, , is key to probing and understanding radiation damage by X-rays and energetic particles in both diffraction and imaging experiments. The models used in protein crystallography for the last 50 years provide good fits to experimental I(q) versus nominal dose data, but have unclear physical significance. More recently, a fit to diffraction and imaging experiments suggested that the maximum tolerable dose varies as q ^-1 or linearly with resolution. Here, it is shown that crystallographic data have been strongly perturbed by the effects of spatially nonuniform crystal irradiation and diffraction during data collection. Reanalysis shows that these data are consistent with a purely exponential local dose dependence, = I ₀(q)exp[-D/D _e(q)], where D _e(q) ∝ q ^α with α ≃ 1.7. A physics-based model for radiation damage, in which damage events occurring at random locations within a sample each cause energy deposition and blurring of the electron density within a small volume, predicts this exponential variation with dose for all q values and a decay exponent α ≃ 2 in two and three dimensions, roughly consistent with both diffraction and imaging experiments over more than two orders of magnitude in resolution. The B-factor model used to account for radiation damage in crystallographic scaling programs is consistent with α = 2, but may not accurately capture the dose dependencies of structure factors under typical nonuniform illumination conditions. The strong q dependence of radiation-induced diffraction decays implies that the previously proposed 20-30 MGy dose limit for protein crystallography should be replaced by a resolution-dependent dose limit that, for atomic resolution data sets, will be much smaller. The results suggest that the physics underlying basic experimental trends in radiation damage at T ≃ 100 K is straightforward and universal. Deviations of the local I(q, D) from strictly exponential behavior may provide mechanistic insights, especially into the radiation-damage processes responsible for the greatly increased radiation sensitivity observed at T ≃ 300 K.

Evolving SAXS versatility: solution X-ray scattering for macromolecular architecture, functional landscapes, and integrative structural biology

Small-angle X-ray scattering (SAXS) has emerged as an enabling integrative technique for comprehensive analyses of macromolecular structures and interactions in solution. Over the past two decades, SAXS has become a mainstay of the structural biologist’s toolbox, supplying multiplexed measurements of molecular shape and dynamics that unveil biological function. Here, we discuss evolving SAXS theory, methods, and applications that extend the field of small-angle scattering beyond simple shape characterization. SAXS, coupled with size-exclusion chromatography (SEC-SAXS) and time-resolved (TR-SAXS) methods, is now providing high-resolution insight into macromolecular flexibility and ensembles, delineating biophysical landscapes, and facilitating high-throughput library screening to assess macromolecular properties and to create opportunities for drug discovery. Looking forward, we consider SAXS in the integrative era of hybrid structural biology methods, its potential for illuminating cellular supramolecular and mesoscale structures, and its capacity to complement high-throughput bioinformatics sequencing data. As advances in the field continue, we look forward to proliferating uses of SAXS based upon its abilities to robustly produce mechanistic insights for biology and medicine.

Challenge data set for macromolecular multi-microcrystallography

A synthetic data set demonstrating a particularly challenging case of indexing ambiguity in the context of radiation damage was generated. This set shall serve as a standard benchmark and reference point for the ongoing development of new methods and new approaches to robust structure solution when single-crystal methods are insufficient. Of the 100 short wedges of data, only the first 36 are currently necessary to solve the structure by `cheating', or using the correct reference structure as a guide. The total wall-clock time and number of crystals required to solve the structure without cheating is proposed as a metric for the efficacy and efficiency of a given multi-crystal automation pipeline.

Structural Control of Nonnative Ligand Binding in Engineered Mutants of Phosphoenolpyruvate Carboxykinase

Protein engineering to alter recognition underlying ligand binding and activity has enormous potential. Here, ligand binding for Escherichia coli phosphoenolpyruvate carboxykinase (PEPCK), which converts oxaloacetate into CO₂ and phosphoenolpyruvate as the first committed step in gluconeogenesis, was engineered to accommodate alternative ligands as an exemplary system with structural information. From our identification of bicarbonate binding in the PEPCK active site at the supposed CO₂ binding site, we probed binding of nonnative ligands with three oxygen atoms arranged to resemble the bicarbonate geometry. Crystal structures of PEPCK and point mutants with bound nonnative ligands thiosulfate and methanesulfonate along with strained ATP and reoriented oxaloacetate intermediates and unexpected bicarbonate were determined and analyzed. The mutations successfully altered the bound ligand position and orientation and its specificity: mutated PEPCKs bound either thiosulfate or methanesulfonate but never both. Computational calculations predicted a methanesulfonate binding mutant and revealed that release of the active site ordered solvent exerts a strong influence on ligand binding. Besides nonnative ligand binding, one mutant altered the Mn²⁺ coordination sphere: instead of the canonical octahedral ligand arrangement, the mutant in question had an only five-coordinate arrangement. From this work, critical features of ligand binding, position, and metal ion cofactor geometry required for all downstream events can be engineered with small numbers of mutations to provide insights into fundamental underpinnings of protein-ligand recognition. Through structural and computational knowledge, the combination of designed and random mutations aids in the robust design of predetermined changes to ligand binding and activity to engineer protein function.

Bringing diffuse X-ray scattering into focus

X-ray crystallography is experiencing a renaissance as a method for probing the protein conformational ensemble. The inherent limitations of Bragg analysis, however, which only reveals the mean structure, have given way to a surge in interest in diffuse scattering, which is caused by structure variations. Diffuse scattering is present in all macromolecular crystallography experiments. Recent studies are shedding light on the origins of diffuse scattering in protein crystallography, and provide clues for leveraging diffuse scattering to model protein motions with atomic detail.

Small angle X-ray scattering and cross-linking for data assisted protein structure prediction in CASP 12 with prospects for improved accuracy

Experimental data offers empowering constraints for structure prediction. These constraints can be used to filter equivalently scored models or more powerfully within optimization functions toward prediction. In CASP12, Small Angle X-ray Scattering (SAXS) and Cross-Linking Mass Spectrometry (CLMS) data, measured on an exemplary set of novel fold targets, were provided to the CASP community of protein structure predictors. As solution-based techniques, SAXS and CLMS can efficiently measure states of the full-length sequence in its native solution conformation and assembly. However, this experimental data did not substantially improve prediction accuracy judged by fits to crystallographic models. One issue, beyond intrinsic limitations of the algorithms, was a disconnect between crystal structures and solution-based measurements. Our analyses show that many targets had substantial percentages of disordered regions (up to 40%) or were multimeric or both. Thus, solution measurements of flexibility and assembly support variations that may confound prediction algorithms trained on crystallographic data and expecting globular fully-folded monomeric proteins. Here, we consider the CLMS and SAXS data collected, the information in these solution measurements, and the challenges in incorporating them into computational prediction. As improvement opportunities were only partly realized in CASP12, we provide guidance on how data from the full-length biological unit and the solution state can better aid prediction of the folded monomer or subunit. We furthermore describe strategic integrations of solution measurements with computational prediction programs with the aim of substantially improving foundational knowledge and the accuracy of computational algorithms for biologically-relevant structure predictions for proteins in solution.

Robust Production, Crystallization, Structure Determination, and Analysis of [Fe-S] Proteins: Uncovering Control of Electron Shuttling and Gating in the Respiratory Metabolism of Molybdopterin Guanine Dinucleotide Enzymes

[Fe-S] clusters are essential cofactors in all domains of life. They play many biological roles due to their unique abilities for electron transfer and conformational control. Yet, producing and analyzing Fe-S proteins can be difficult and even misleading if not done anaerobically. Due to unique redox properties of [Fe-S] clusters and their oxygen sensitivity, they pose multiple challenges and can lose enzymatic activity or cause their component proteins to be structurally disordered due to [Fe-S] cluster oxidation and loss in air. Here we highlight tested protocols and strategies enabling efficient and stable [Fe-S] protein production, purification, crystallization, X-ray diffraction data collection, and structure determination. From multiple high-resolution anaerobic crystal structures, we furthermore analyze exemplary data defining [Fe-S] clusters, substrate entry, and product exit for the functional oxidation states of type II molybdo-bis(molybdopterin guanine dinucleotide) (Mo-bisMGD) enzymes. Notably, these enzymes perform electron shuttling between quinone pools and specific substrates to catalyze respiratory metabolism. The identified structure-activity relationships for this enzyme class have broad implications germane to perchlorate environments on Earth and Mars extending to an alternative mechanism underlying metabolic origins for the evolution of the oxygen atmosphere. Integrated structural analyses of type II Mo-bisMGD enzymes unveil novel distinctive shared molecular mechanisms for dynamic control of substrate entry and product release gated by hydrophobic residues. Collective findings support a prototypic model for type II Mo-bisMGD enzymes including insights for a fundamental molecular mechanistic understanding of selectivity and regulation by a conformationally gated channel with general implications for [Fe-S] cluster respiratory enzymes.

Lifetimes and spatio-temporal response of protein crystals in intense X-ray microbeams

Serial synchrotron-based crystallography using intense microfocused X-ray beams, fast-framing detectors and protein microcrystals held at 300 K promises to expand the range of accessible structural targets and to increase overall structure-pipeline throughputs. To explore the nature and consequences of X-ray radiation damage under microbeam illumination, the time-, dose- and temperature-dependent evolution of crystal diffraction have been measured with maximum dose rates of 50 MGy s-1. At all temperatures and dose rates, the integrated diffraction intensity for a fixed crystal orientation shows non-exponential decays with dose. Non-exponential decays are a consequence of non-uniform illumination and the resulting spatial evolution of diffracted intensity within the illuminated crystal volume. To quantify radiation-damage lifetimes and the damage state of diffracting crystal regions, a revised diffraction-weighted dose (DWD) is defined and it is shown that for Gaussian beams the DWD becomes nearly independent of actual dose at large doses. An apparent delayed onset of radiation damage seen in some intensity-dose curves is in fact a consequence of damage. Intensity fluctuations at high dose rates may arise from the impulsive release of gaseous damage products. Accounting for these effects, data collection at the highest dose rates increases crystal radiation lifetimes near 300 K (but not at 100 K) by a factor of ∼1.5-2 compared with those observed at conventional dose rates. Improved quantification and modeling of the complex spatio-temporal evolution of protein microcrystal diffraction in intense microbeams will enable more efficient data collection, and will be essential in improving the accuracy of structure factors and structural models.

A multi-crystal method for extracting obscured crystallographic states from conventionally uninterpretable electron density

In macromolecular crystallography, the rigorous detection of changed states (for example, ligand binding) is difficult unless signal is strong. Ambiguous (‘weak' or ‘noisy') density is experimentally common, since molecular states are generally only fractionally present in the crystal. Existing methodologies focus on generating maximally accurate maps whereby minor states become discernible; in practice, such map interpretation is disappointingly subjective, time-consuming and methodologically unsound. Here we report the PanDDA method, which automatically reveals clear electron density for the changed state—even from inaccurate maps—by subtracting a proportion of the confounding ‘ground state'; changed states are objectively identified from statistical analysis of density distributions. The method is completely general, implying new best practice for all changed-state studies, including the routine collection of multiple ground-state crystals. More generally, these results demonstrate: the incompleteness of atomic models; that single data sets contain insufficient information to model them fully; and that accuracy requires further map-deconvolution approaches.

Building a ligand into a weak region of an electron density map of a protein is a subjective process. Here, the authors present a new method to obtain a clear electron density for a bound ligand based on multi-crystal experiments and 3D background correction.

The N14 anti-afamin antibody Fab: a rare VL1 CDR glycosylation, crystallographic re-sequencing, molecular plasticity and conservative versus enthusiastic modelling

The monoclonal antibody N14 is used as a detection antibody in ELISA kits for the human glycoprotein afamin, a member of the albumin family, which has recently gained interest in the capture and stabilization of Wnt signalling proteins, and for its role in metabolic syndrome and papillary thyroid carcinoma. As a rare occurrence, the N14 Fab is N-glycosylated at Asn26L at the onset of the VL1 antigen-binding loop, with the α-1-6 core fucosylated complex glycan facing out of the L1 complementarity-determining region. The crystal structures of two non-apparent (pseudo) isomorphous crystals of the N14 Fab were analyzed, which differ significantly in the elbow angles, thereby cautioning against the overinterpretation of domain movements upon antigen binding. In addition, the map quality at 1.9 Å resolution was sufficient to crystallographically re-sequence the variable VL and VH domains and to detect discrepancies in the hybridoma-derived sequence. Finally, a conservatively refined parsimonious model is presented and its statistics are compared with those from a less conservatively built model that has been modelled more enthusiastically. Improvements to the PDB validation reports affecting ligands, clashscore and buried surface calculations are suggested.