Density modification is a standard step to provide a route for routine structure solution by any experimental phasing method, with single-wavelength or multi-wavelength anomalous diffraction being the most popular methods, as well as to extend fragments or incomplete models into a full solution. The effect of density modification on the starting maps from either source is illustrated in the case of SHELXE . The different modes in which the program can run are reviewed; these include less well known uses such as reading external phase values and weights or phase distributions encoded in Hendrickson–Lattman coefficients. Typically in SHELXE , initial phases are calculated from experimental data, from a partial model or map, or from a combination of both sources. The initial phase set is improved and extended by density modification and, if the resolution of the data and the type of structure permits, polyalanine tracing. As a feature to systematically eliminate model bias from phases derived from predicted models, the trace can be set to exclude the area occupied by the starting model. The trace now includes an extension into the gamma position or hydrophobic and aromatic side chains if a sequence is provided, which is performed in every tracing cycle. Once a correlation coefficient of over 30% between the structure factors calculated from such a trace and the native data indicates that the structure has been solved, the sequence is docked in all model-building cycles and side chains are fitted if the map supports it. The extensions to the tracing algorithm brought in to provide a complete model are discussed. The improvement in phasing performance is assessed using a set of tests.

The Collaborative Computational Project No. 4 (CCP4) is a UK-led international collective with a mission to develop, test, distribute and promote software for macromolecular crystallography. The CCP4 suite is a multiplatform collection of programs brought together by familiar execution routines, a set of common libraries and graphical interfaces. The CCP4 suite has experienced several considerable changes since its last reference article, involving new infrastructure, original programs and graphical interfaces. This article, which is intended as a general literature citation for the use of the CCP4 software suite in structure determination, will guide the reader through such transformations, offering a general overview of the new features and outlining future developments. As such, it aims to highlight the individual programs that comprise the suite and to provide the latest references to them for perusal by crystallographers around the world.

Genetic CLN5 variants are associated with childhood neurodegeneration and Alzheimer’s disease; however, the molecular function of ceroid lipofuscinosis neuronal protein 5 (Cln5) is unknown. We solved the Cln5 crystal structure and identified a region homologous to the catalytic domain of members of the N1pC/P60 superfamily of papain-like enzymes. However, we observed no protease activity for Cln5; and instead, we discovered that Cln5 and structurally related PPPDE1 and PPPDE2 have efficient cysteine palmitoyl thioesterase (S-depalmitoylation) activity using fluorescent substrates. Mutational analysis revealed that the predicted catalytic residues histidine-166 and cysteine-280 are critical for Cln5 thioesterase activity, uncovering a new cysteine-based catalytic mechanism for S-depalmitoylation enzymes. Last, we found that Cln5-deficient neuronal progenitor cells showed reduced thioesterase activity, confirming live cell function of Cln5 in setting S-depalmitoylation levels. Our results provide new insight into the function of Cln5, emphasize the importance of S-depalmitoylation in neuronal homeostasis, and disclose a new, unexpected enzymatic function for the N1pC/P60 superfamily of proteins.

Density modification is a standard step to provide a route for routine structure solution by any experimental phasing method -with SAD and MAD being the most popular ones- as well as to extend fragments or incomplete models into a full solution. The effect of density modification on the starting maps from either source is illustrated in the case of SHELXE. The different modes in which the program can run are reviewed; these include less well-known uses such as reading external phase values and weights or phase distributions encoded in Hendrickson-Lattman coefficients. Typically in SHELXE, initial phases are calculated from experimental data, from a partial model or map, or from a combination of both sources. The initial phase set is improved and extended by density modification and, if the resolution of the data and the type of structure permits, poly-alanine tracing. The trace now includes an extension into the gamma position or hydrophobic and aromatic side chains if a sequence is provided, which is performed in every tracing cycle. Once a correlation coefficient over 30% between the structure factors calculated from such a trace and the native data indicates that the structure has been solved, in all model building cycles sequence is docked and side chains are fitted if the map supports it. The extensions to the tracing algorithm brought in to provide a complete model are discussed. The improvement in phasing performance is assessed using a set of tests. Synopsis Side chain tracing now completes model building in SHELXE to enhance density modification. All alternative SHELXE modes, using single or combined sources of starting phase information, are described. ![Formula][1] ### Competing Interest Statement The authors have declared no competing interest. * CC : Correlation coefficient expressed as percentage [1]: /embed/graphic-1.gif

A new aspherical scattering factor formalism has been implemented in the crystallographic least-squares refinement program SHELXL. The formalism relies on Gaussian functions and can optionally complement the independent atom model to take into account the deformation of electron-density distribution due to chemical bonding and lone pairs. Asphericity contributions were derived from the electron density obtained from quantum-chemical density functional theory computations of suitable model compounds that contain particular chemical environments, as defined by the invariom formalism. Thanks to a new algorithm, invariom assignment for refinement in SHELXL is automated. A suitable parameterization for each chemical environment within the new model was achieved by metaheuristics. Figures of merit, precision and accuracy of crystallographic least-squares refinements improve significantly upon using the new model.

In contrast to twinning by merohedry, the reciprocal lattices of the different domains of non-merohedral twins do not overlap exactly. This leads to three kinds of reflections: reflections with no overlap, reflections with an exact overlap and reflections with a partial overlap of a reflection from a second domain. This complicates the unit-cell determination, indexing, data integration and scaling of X-ray diffraction data. However, with hindsight it is possible to detwin the data because there are reflections that are not affected by the twinning. In this article, the successful solution and refinement of one mineral, one organometallic and two protein non-merohedral twins using a common strategy are described. The unit-cell constants and the orientation matrices were determined by the program CELL_NOW. The data were then integrated with SAINT. TWINABS was used for scaling, empirical absorption corrections and the generation of two different data files, one with detwinned data for structure solution and refinement and a second one for (usually more accurate) structure refinement against total integrated intensities. The structures were solved by experimental phasing using SHELXT for the first two structures and SHELXC/D/E for the two protein structures; all models were refined with SHELXL.

The open-source Python program PDB2INS is designed to prepare a .ins file for refinement with SHELXL [Sheldrick (2015). Acta Cryst. C71, 3-8], taking atom coordinates and other information from a Protein Data Bank (PDB)-format file. If PDB2INS is provided with a four-character PDB code, both the PDB file and the accompanying mmCIF-format reflection data file (if available) are accessed via the internet from the PDB public archive [Read et al. (2011). Structure, 19, 1395-1412] or optionally from the PDB_REDO server [Joosten, Long, Murshudov & Perrakis (2014). IUCrJ, 1, 213-220]. The SHELX-format .ins (refinement instructions and atomic coordinates) and .hkl (reflection data) files can then be generated without further user intervention, appropriate restraints etc. being added automatically. PDB2INS was tested on the 23 974 X-ray structures deposited in the PDB between 2008 and 2018 that included reflection data to 1.7 Å or better resolution in a recognizable format. After creating the two input files for SHELXL without user intervention, ten cycles of conjugate-gradient least-squares refinement were performed. For 96% of these structures PDB2INS and SHELXL completed successfully without error messages.

For the purpose of this article, experimental phasing is understood to mean the determination of macromolecular structures by exploiting small intensity differences of Friedel opposites and possibly of reflections measured at different wavelengths or for heavy-atom derivatives, without the use of specific structural models. TheSHELXprograms provide a robust and efficient route for routine structure solution by the SAD, MAD and related methods, but involve a number of simplifying assumptions that may limit their applicability in borderline cases. The substructure atoms (i.e.those with significant anomalous scattering) are first located by direct methods, and the experimental data are then used to estimate phase shifts that are added to the substructure phases to obtain starting phases for the native reflections. These are then improved by density modification and, if the resolution of the data and the type of structure permit, polyalanine tracing. A number of extensions to the tracing algorithm are discussed; these are designed to improve its performance at low resolution. Given native data to 2.5 Å resolution or better, a correlation coefficient greater than 25% between the structure factors calculated from such a trace and the native data is usually a good indication that the structure has been solved.

The CCP4 (Collaborative Computational Project, Number 4) software suite for macromolecular structure determination by X-ray crystallography groups brings together many programs and libraries that, by means of well established conventions, interoperate effectively without adhering to strict design guidelines. Because of this inherent flexibility, users are often presented with diverse, even divergent, choices for solving every type of problem. Recently, CCP4 introduced CCP4i2, a modern graphical interface designed to help structural biologists to navigate the process of structure determination, with an emphasis on pipelining and the streamlined presentation of results. In addition, CCP4i2 provides a framework for writing structure-solution scripts that can be built up incrementally to create increasingly automatic procedures.

As early as 1961, Rich, Davies, Crick and Watson proposed a parallel double helix RNA structure for poly(rA) on the basis of diffuse fibre diffraction photographs [1], using similar arguments to those used to deduce the antiparallel double helix structure for DNA.It took over 50 years before this could be confirmed by a crystal structure determination [2].However the solution of this structure by ab initio direct methods required one week on an eight-CPU workstation.We recently obtained synchrotron data from better quality crystals of significantly higher resolution and completeness.With the new data the structure can be solved in a much shorter time and the data quality appears to be sufficient for a charge density analysis using the Invariom approach [3].

Figure 1.The Fragon process.The correct solution can be selected based on the correlation coefficient (CC) after density modification.

The quality of diffraction data obtained using silver and molybdenum microsources has been compared for six model compounds with a wide range of absorption factors. The experiments were performed on two 30 W air-cooled Incoatec IµS microfocus sources with multilayer optics mounted on a Bruker D8 goniometer with a SMART APEX II CCD detector. All data were analysed, processed and refined using standard Bruker software. The results show that Ag Kα radiation can be beneficial when heavy elements are involved. A numerical absorption correction based on the positions and indices of the crystal faces is shown to be of limited use for the highly focused microsource beams, presumably because the assumption that the crystal is completely bathed in a (top-hat profile) beam of uniform intensity is no longer valid. Fortunately the empirical corrections implemented inSADABS, although originally intended as a correction for absorption, also correct rather well for the variations in the effective volume of the crystal irradiated. In three of the cases studied (two Ag and one Mo) the finalSHELXL R1 against all data after application of empirical corrections implemented inSADABSwas below 1%. Since such corrections are designed to optimize the agreement of the intensities of equivalent reflections with different paths through the crystal but the same Bragg 2θ angles, a further correction is required for the 2θ dependence of the absorption. For this,SADABSuses the transmission factor of a spherical crystal with a user-defined value of μr(where μ is the linear absorption coefficient andris the effective radius of the crystal); the best results are obtained whenris biased towards the smallest crystal dimension. The results presented here suggest that the IUCr publication requirement that a numerical absorption correction must be applied for strongly absorbing crystals is in need of revision.

The improvements in the crystal structure refinement program SHELXL have been closely coupled with the development and increasing importance of the CIF (Crystallographic Information Framework) format for validating and archiving crystal structures. An important simplification is that now only one file in CIF format (for convenience, referred to simply as 'a CIF') containing embedded reflection data and SHELXL instructions is needed for a complete structure archive; the program SHREDCIF can be used to extract the .hkl and .ins files required for further refinement with SHELXL. Recent developments in SHELXL facilitate refinement against neutron diffraction data, the treatment of H atoms, the determination of absolute structure, the input of partial structure factors and the refinement of twinned and disordered structures. SHELXL is available free to academics for the Windows, Linux and Mac OS X operating systems, and is particularly suitable for multiple-core processors.

ARCIMBOLDO allows ab initio phasing of macromolecular structures below atomic resolution by exploiting the location of small model fragments combined with density modification in a multisolution frame. The model fragments can be either secondary-structure elements predicted from the sequence or tertiary-structure fragments. The latter can be derived from libraries of typical local folds or from related structures, such as a low-homology model that is unsuccessful in molecular replacement. In all ARCIMBOLDO applications, fragments are searched for sequentially. Correct partial solutions obtained after each fragment-search stage but lacking the necessary phasing power can, if combined, succeed. Here, an analysis is presented of the clustering of partial solutions in reciprocal space and of its application to a set of different cases. In practice, the task of combining model fragments from an ARCIMBOLDO run requires their referral to a common origin and is complicated by the presence of correct and incorrect solutions as well as by their not being independent. The F-weighted mean phase difference has been used as a figure of merit. Clustering perfect, non-overlapping fragments dismembered from test structures in polar and nonpolar space groups shows that density modification before determining the relative origin shift enhances its discrimination. In the case of nonpolar space groups, clustering of ARCIMBOLDO solutions from secondary-structure models is feasible. The use of partially overlapping search fragments provides a more favourable circumstance and was assessed on a test case. Applying the devised strategy, a previously unknown structure was solved from clustered correct partial solutions.

K+ channels exist in virtually all living organisms and play fundamental biological roles by conducting K+ ions in a rapid and selective way. Previous models proposed that K+ ions and water molecules alternately occupy four binding sites of the selectivity filter and permeate collectively when the channel is open and activated. However, our latest computational electrophysiology molecular dynamics simulations [1] indicate that this model can not yield measurable conductance under physiologically relevant transmembrane potential. Instead, the rapid ion permeation only occurs when there are no water molecules in the selectivity filter and K+ ions form direct contacts. This has been further confirmed by Brownian dynamics simulations in a simplified model. A revisit to the previous X-ray crystallography data has validated the direct contacts of K+ ions in the selectivity filter, which is in good agreement with classical and QM/MM molecular dynamics simulations. Therefore, we propose a direct Column knock-on mechanism that can account for the high ion conductance of K+ channels. [1] C. Kutzner et al., Biophys. J. 101, 809-817 (2011).