The fully automatic processing of crystals of macromolecules has presented a unique opportunity to gather information on the samples that is not usually recorded. This has proved invaluable in improving sample-location, characterization and data-collection algorithms. After operating for four years, MASSIF-1 has now processed over 56 000 samples, gathering information at each stage, from the volume of the crystal to the unit-cell dimensions, the space group, the quality of the data collected and the reasoning behind the decisions made in data collection. This provides an unprecedented opportunity to analyse these data together, providing a detailed landscape of macromolecular crystals, intimate details of their contents and, importantly, how the two are related. The data show that mosaic spread is unrelated to the size or shape of crystals and demonstrate experimentally that diffraction intensities scale in proportion to crystal volume and molecular weight. It is also shown that crystal volume scales inversely with molecular weight. The results set the scene for the development of X-ray crystallography in a changing environment for structural biology.

Copper-containing nitrite reductases (CuNiRs) that convert NO2− to NO via a CuCAT–His–Cys–CuET proton-coupled redox system are of central importance in nitrogen-based energy metabolism. These metalloenzymes, like all redox enzymes, are very susceptible to radiation damage from the intense synchrotron-radiation X-rays that are used to obtain structures at high resolution. Understanding the chemistry that underpins the enzyme mechanisms in these systems requires resolutions of better than 2 Å. Here, for the first time, the damage-free structure of the resting state of one of the most studied CuNiRs was obtained by combining X-ray free-electron laser (XFEL) and neutron crystallography. This represents the first direct comparison of neutron and XFEL structural data for any protein. In addition, damage-free structures of the reduced and nitrite-bound forms have been obtained to high resolution from cryogenically maintained crystals by XFEL crystallography. It is demonstrated that AspCAT and HisCAT are deprotonated in the resting state of CuNiRs at pH values close to the optimum for activity. A bridging neutral water (D2O) is positioned with one deuteron directed towards AspCAT Oδ1 and one towards HisCAT N∊2. The catalytic T2Cu-ligated water (W1) can clearly be modelled as a neutral D2O molecule as opposed to D3O+ or OD−, which have previously been suggested as possible alternatives. The bridging water restricts the movement of the unprotonated AspCAT and is too distant to form a hydrogen bond to the O atom of the bound nitrite that interacts with AspCAT. Upon the binding of NO2− a proton is transferred from the bridging water to the Oδ2 atom of AspCAT, prompting electron transfer from T1Cu to T2Cu and reducing the catalytic redox centre. This triggers the transfer of a proton from AspCAT to the bound nitrite, enabling the reaction to proceed.

Meta-magnetic shape-memory alloys combine ferroelastic order with ferromagnetic order and exhibit attractive multifunctional properties, but they are extremely brittle, showing hardly any tensile deformability, which impedes their practical application. Here, for the first time, an Ni–Cu–Co–Mn–In microwire has been developed that simultaneously exhibits a magnetic field-induced first-order meta-magnetic phase transition and huge tensile superelasticity. A temperature-dependent in situ synchrotron high-energy X-ray diffraction investigation reveals that the martensite of this Ni43.7Cu1.5Co5.1Mn36.7In13 microwire shows a monoclinic six-layered modulated structure and the austenite shows a cubic structure. This microwire exhibits an oligocrystalline structure with bamboo grains, which remarkably reduces the strain incompatibility during deformation and martensitic transformation. As a result, huge tensile superelasticity with a recoverable strain of 13% is achieved in the microwire. This huge tensile superelasticity is in agreement with our theoretical calculations based on the crystal structure and lattice correspondence of austenite and martensite and the crystallographic orientation of the grains. Owing to the large magnetization difference between austenite and martensite, a pronounced magnetic field-induced magnetostructural transition is achieved in the microwire, which could give rise to a variety of magnetically driven functional properties. For example, a large magnetocaloric effect with an isothermal entropy change of 12.7 J kg−1 K−1 (under 5 T) is obtained. The realization of magnetic-field- and tensile-stress-induced structural transformations in the microwire may pave the way for exploiting the multifunctional properties under the coupling of magnetic field and stress for applications in miniature multifunctional devices.

Reliable sample delivery and efficient use of limited beam time have remained bottlenecks for serial crystallography (SX). Using a high-intensity polychromatic X-ray beam in combination with a newly developed charge-integrating JUNGFRAU detector, we have applied the method of fixed-target SX to collect data at a rate of 1 kHz at a synchrotron-radiation facility. According to our data analysis for the given experimental conditions, only about 3 000 diffraction patterns are required for a high-quality diffraction dataset. With indexing rates of up to 25%, recording of such a dataset takes less than 30 s.

The present work reports on the charge and spin density modelling of YTiO3 in its ferromagnetic state (TC = 27 K). Accurate polarized neutron diffraction and high-resolution X-ray diffraction (XRD) experiments were carried out on a single crystal at the ORPHÉE reactor (LLB) and SPRING8 synchrotron source. The experimental data are modelled by the spin resolved pseudo-atomic multipolar model (Deutsch et al., 2012). The refinement strategy is discussed and the result of this electron density modelling is compared with that from XRD measured at 100 K and with density functional theory calculations. The results show that the spin and charge densities around the Ti atom have lobes directed away from the O atoms, confirming the filling of the t2g orbitals of the Ti atom. The dxy orbital is less populated than dxz and dyz, which is a sign of a partial lift of degeneracy of the t2g orbitals. This study confirms the orbital ordering at low temperature (20 K), which is already present in the paramagnetic state above the ferromagnetic transition (100 K).

For Heusler-type Ni–Mn–Ga ferromagnetic shape-memory alloys, the configuration of the martensite variants is a decisive factor in achieving a large magnetic shape-memory effect through field-induced variant reorientation. Based upon the spatially resolved electron backscatter diffraction technique, the microstructural evolution associated with the martensitic transformation from austenite to seven-layered modulated (7M) martensite was investigated on a polycrystalline Ni53Mn22Ga25 alloy. It was clearly shown that grain interior nucleation led to the formation of diamond-shaped 7M martensite within the parent austenite matrix. This diamond microstructure underwent further growth through an isotropic expansion with the coordinated outward movement of four side habit planes, followed by an anisotropic elongation with the forward extension of a type-I twin pair. A two-step growth model is proposed to describe the specific morphology and crystallography of 7M martensite. In addition, the habit planes were revealed to possess a stepped structure, with the {1 0 1}A plane as the terrace and the {0 1 0}A plane as the step. The characteristic combination of martensite variants and the underlying mechanism of self-accommodation in the martensitic transformation have been analysed in terms of the minimum total transformation strain, where the deformation gradient matrix was constructed according to the experimentally determined orientation relationship between the two phases. The present results may deepen the understanding of special martensite microstructures during the martensitic transformation in ferromagnetic shape-memory alloys.

The interoperability of chemical and biological crystallographic data is a key challenge to research and its application to pharmaceutical design. Research attempting to combine data from the two disciplines, small-molecule or chemical crystallography (CX) and macromolecular crystallography (MX), will face unique challenges including variations in terminology, software development, file format and databases which differ significantly from CX to MX. This perspective overview spans the two disciplines and originated from the investigation of protein binding to model radiopharmaceuticals. The opportunities of interlinked research while utilizing the two databases of the CSD (Cambridge Structural Database) and the PDB (Protein Data Bank) will be highlighted. The advantages of software that can handle multiple file formats and the circuitous route to convert organometallic small-molecule structural data for use in protein refinement software will be discussed. In addition some pointers to avoid being shipwrecked will be shared, such as the care which must be taken when interpreting data precision involving small molecules versus proteins.

As is well known, polymers commonly form lamellar crystals, and these assemble further into lamellar stacks and spherulites during quiescent crystallization. Fifty years ago, Vonk and Kortleve constructed the classical small-angle X-ray scattering theory (SAXS) for a lamellar system, in which it was assumed that the lamellar stack had an infinite lateral size [Vonk & Kortleve (1967), Kolloid Z. Z. Polym. 220, 19–24]. Under this assumption, only crystal planes satisfying the Bragg condition can form strong scattering, and the scattering from the lamellar stack arises from the difference between the scattering intensities in the amorphous and crystalline layers, induced by the incident X-ray beam. This assumption is now deemed unreasonable. In a real polymer spherulite, the lamellar crystal commonly has dimensions of only a few hundred nanometres. At such a limited lateral size, lamellar stacks in a broad orientation have similar scattering, so interference between these lamellar stacks must be considered. Scattering from lamellar stacks parallel to the incident X-ray beam also needs to be considered when total reflection occurs. In this study, various scattering contributions from lamellar stacks in a spherulite are determined. It is found that, for a limited lateral size, the scattering induced by the incident X-ray beam is not the main origin of SAXS. It forms double peaks, which are not observed in real scattering because of destructive interference between the lamellar stacks. The scattering induced by the evanescent wave is the main origin. It can form a similar interference pattern to that observed in a real SAXS measurement: a Guinier region in the small-q range, a signal region in the intermediate-q range and a Porod region in the high-q range. It is estimated that, to avoid destructive interference, the lateral size needs to be greater than 11 µm, which cannot be satisfied in a real lamellar system. Therefore, SAXS in a real polymer system arises largely from the scattering induced by the evanescent wave. Evidence for the existence of the evanescent wave was identified in the scattering of isotactic polypropyl­ene. This study corrects a long-term misunderstanding of SAXS in a polymer lamellar system.

High-throughput X-ray crystal structures of protein–ligand complexes are critical to pharmaceutical drug development. However, cryocooling of crystals and X-ray radiation damage may distort the observed ligand binding. Serial femtosecond crystallography (SFX) using X-ray free-electron lasers (XFELs) can produce radiation-damage-free room-temperature structures. Ligand-binding studies using SFX have received only modest attention, partly owing to limited beamtime availability and the large quantity of sample that is required per structure determination. Here, a high-throughput approach to determine room-temperature damage-free structures with excellent sample and time efficiency is demonstrated, allowing complexes to be characterized rapidly and without prohibitive sample requirements. This yields high-quality difference density maps allowing unambiguous ligand placement. Crucially, it is demonstrated that ligands similar in size or smaller than those used in fragment-based drug design may be clearly identified in data sets obtained from <1000 diffraction images. This efficiency in both sample and XFEL beamtime opens the door to true high-throughput screening of protein–ligand complexes using SFX.

The disposition of functional groups can induce variations in the nature and type of interactions and hence affect the molecular recognition and self-assembly mechanism in cocrystals. To better understand the formation of cocrystals on a molecular level, the effects of disposition of functional groups on the formation of cocrystals were systematically and comprehensively investigated using cresol isomers (o-, m-, p-cresol) as model compounds. Consistency and variability in these cocrystals containing positional isomers were found and analyzed. The structures, molecular recognition and self-assembly mechanism of supramolecular synthons in solution and in their corresponding cocrystals were verified by a combined experimental and theoretical calculation approach. It was found that the heterosynthons (heterotrimer or heterodimer) combined with O—H⋯N hydrogen bonding played a significant role. Hirshfeld surface analysis and computed interaction energy values were used to determine the hierarchical ordering of the weak interactions. The quantitative analyses of charge transfers and molecular electrostatic potential were also applied to reveal and verify the reasons for consistency and variability. Finally, the molecular recognition, self-assembly and evolution process of the supramolecular synthons in solution were investigated. The results confirm that the supramolecular synthon structures formed initially in solution would be carried over to the final cocrystals, and the supramolecular synthon structures are the precursors of cocrystals and the information memory of the cocrystallization process, which is evidence for classical nucleation theory.

Indanomycin is biosynthesized by a hybrid nonribosomal peptide synthase/polyketide synthase (NRPS/PKS) followed by a number of `tailoring' steps to form the two ring systems that are present in the mature product. It had previously been hypothesized that the indane ring of indanomycin was formed by the action of IdmH using a Diels–Alder reaction. Here, the crystal structure of a selenomethionine-labelled truncated form of IdmH (IdmH-Δ99–107) was solved using single-wavelength anomalous dispersion (SAD) phasing. This truncated variant allows consistent and easy crystallization, but importantly the structure was used as a search model in molecular replacement, allowing the full-length IdmH structure to be determined to 2.7 Å resolution. IdmH is a homodimer, with the individual protomers consisting of an α+β barrel. Each protomer contains a deep hydrophobic pocket which is proposed to constitute the active site of the enzyme. To investigate the reaction catalysed by IdmH, 88% of the backbone NMR resonances were assigned, and using chemical shift perturbation of [15N]-labelled IdmH it was demonstrated that indanomycin binds in the active-site pocket. Finally, combined quantum mechanical/molecular mechanical (QM/MM) modelling of the IdmH reaction shows that the active site of the enzyme provides an appropriate environment to promote indane-ring formation, supporting the assignment of IdmH as the key Diels–Alderase catalysing the final step in the biosynthesis of indanomycin through a similar mechanism to other recently characterized Diels–Alderases involved in polyketide-tailoring reactions. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at https://proteopedia.org/w/Journal:IUCrJ:S2052252519012399.

In the past three years, Dirac half-metals (DHMs) have attracted considerable attention and become a high-profile topic in spintronics becuase of their excellent physical properties such as 100% spin polarization and massless Dirac fermions. Two-dimensional DHMs proposed recently have not yet been experimentally synthesized and thus remain theoretical. As a result, their characteristics cannot be experimentally confirmed. In addition, many theoretically predicted Dirac materials have only a single cone, resulting in a nonlinear electromagnetic response with insufficient intensity and inadequate transport carrier efficiency near the Fermi level. Therefore, after several attempts, we have focused on a novel class of DHMs with multiple Dirac crossings to address the above limitations. In particular, we direct our attention to three-dimensional bulk materials. In this study, the discovery via first principles of an experimentally synthesized DHM LaNiO3 with many Dirac cones and complete spin polarization near the Fermi level is reported. It is also shown that the crystal structures of these materials are strongly correlated with their physical properties. The results indicate that many rhombohedral materials with the general formula LnNiO3 (Ln = La, Ce, Nd, Pm, Gd, Tb, Dy, Ho, Er, Lu) in the space group R3c are potential DHMs with multiple Dirac cones.

The class B family of G-protein-coupled receptors (GPCRs) has long been a paradigm for peptide hormone recognition and signal transduction. One class B GPCR, the glucagon-like peptide-1 receptor (GLP-1R), has been considered as an anti-diabetes drug target and there are several peptidic drugs available for the treatment of this overwhelming disease. The previously determined structures of inactive GLP-1R in complex with two negative allosteric modulators include ten thermal-stabilizing mutations that were selected from a total of 98 designed mutations. Here we systematically summarize all 98 mutations we have tested and the results suggest that the mutagenesis strategy that strengthens inter-helical hydro­phobic interactions shows the highest success rate. We further investigate four back mutations by thermal-shift assay, crystallization and molecular dynamic simulations, and conclude that mutation I1962.66bF increases thermal stability intrinsically and that mutation S2714.47bA decreases crystal packing entropy extrinsically, while mutations S1932.63bC and M2333.36bC may be dispensable since these two cysteines are not di­sulfide-linked. Our results indicate intrinsic connections between different regions of GPCR transmembrane helices and the current data suggest a general mutagenesis principle for structural determination of GPCRs and other membrane proteins.

Tannerella forsythia is an oral dysbiotic periodontopathogen involved in severe human periodontal disease. As part of its virulence factor armamentarium, at the site of colonization it secretes mirolysin, a metallopeptidase of the unicellular pappalysin family, as a zymogen that is proteolytically auto-activated extracellularly at the Ser54–Arg55 bond. Crystal structures of the catalytically impaired promirolysin point mutant E225A at 1.4 and 1.6 Å revealed that latency is exerted by an N-terminal 34-residue pro-segment that shields the front surface of the 274-residue catalytic domain, thus preventing substrate access. The catalytic domain conforms to the metzincin clan of metallopeptidases and contains a double calcium site, which acts as a calcium switch for activity. The pro-segment traverses the active-site cleft in the opposite direction to the substrate, which precludes its cleavage. It is anchored to the mature enzyme through residue Arg21, which intrudes into the specificity pocket in cleft sub-site S1'. Moreover, residue Cys23 within a conserved cysteine–glycine motif blocks the catalytic zinc ion by a cysteine-switch mechanism, first described for mammalian matrix metallopeptidases. In addition, a 1.5 Å structure was obtained for a complex of mature mirolysin and a tetradecapeptide, which filled the cleft from sub-site S1' to S6'. A citrate molecule in S1 completed a product-complex mimic that unveiled the mechanism of substrate binding and cleavage by mirolysin, the catalytic domain of which was already preformed in the zymogen. These results, including a preference for cleavage before basic residues, are likely to be valid for other unicellular pappalysins derived from archaea, bacteria, cyanobacteria, algae and fungi, including archetypal ulilysin from Methanosarcina acetivorans. They may further apply, at least in part, to the multi-domain orthologues of higher organisms.

For serial femtosecond crystallography at X-ray free-electron lasers, which entails collection of single-pulse diffraction patterns from a constantly refreshed supply of microcrystalline sample, delivery of the sample into the X-ray beam path while maintaining low background remains a technical challenge for some experiments, especially where this methodology is applied to relatively low-ordered samples or those difficult to purify and crystallize in large quantities. This work demonstrates a scheme to encapsulate biological samples using polymer thin films and graphene to maintain sample hydration in vacuum conditions. The encapsulated sample is delivered into the X-ray beam on fixed targets for rapid scanning using the Roadrunner fixed-target system towards a long-term goal of low-background measurements on weakly diffracting samples. As a proof of principle, we used microcrystals of the 24 kDa rapid encystment protein (REP24) to provide a benchmark for polymer/graphene sandwich performance. The REP24 microcrystal unit cell obtained from our sandwiched in-vacuum sample was consistent with previously established unit-cell parameters and with those measured by us without encapsulation in humidified helium, indicating that the platform is robust against evaporative losses. While significant scattering from water was observed because of the sample-deposition method, the polymer/graphene sandwich itself was shown to contribute minimally to background scattering.

Although a plethora of metal complexes have been characterized, those having multifunctional properties are very rare. This article reports three isotypical complexes, namely [Cu(benzoate)L2], where L = 4-styryl­pyridine (4spy) (1), 2'-fluoro-4-styryl­pyridine (2F-4spy) (2) and 3'-fluoro-4-styryl­pyridine (3F-4spy) (3), which show photosalient behavior (photoinduced crystal mobility) while they undergo [2+2] cyclo­addition. These crystals also exhibit anisotropic thermal expansion when heated from room temperature to 200°C. The overall thermal expansion of the crystals is impressive, with the largest volumetric thermal expansion coefficients for 1, 2 and 3 of 241.8, 233.1 and 285.7 × 10−6 K−1, respectively, values that are comparable to only a handful of other reported materials known to undergo colossal thermal expansion. As a result of the expansion, their single crystals occasionally move by rolling. Altogether, these materials exhibit unusual and hitherto untapped solid-state properties.

The preferred orientation growth characteristics and surface roughness of polycrystalline bis­muth (Bi) thin films fabricated on glass substrates using the molecular beam epitaxy method were investigated at temperatures ranging from 18 to 150°C. The crystallization and morphology were analyzed in detail and the polycrystalline metal film structure-zone model (SZM) was modified to fit the polycrystalline Bi thin film. The boundary temperature between Zone T and Zone II in the SZM shifted to higher temperatures with the increase in film thickness or the decrease of growth rate. Furthermore, the effect of the thickness and surface roughness on the transport properties was investigated, especially for Bi thin films in Zone II. A two-transport channels model was adopted to reveal the influence of the film thickness on the competition between the metallic surface states and the semiconducting bulk states, which is consistent with the results of Bi single-crystal films. Therefore, the polycrystalline Bi thin films are expected to replace the single-crystal films in the application of spintronic devices.

Using the typical WC–Co cemented carbide as an example, the interactions of dislocations within the ceramic matrix and the binder metal, as well as the possible cooperation and competition between the matrix and binder during deformation of the nanocrystalline cermets, were studied by molecular dynamics simulations. It was found that at the same level of strain, the dislocations in Co have more complex configurations in the cermet with higher Co content. With loading, the ratio between mobile and sessile dislocations in Co becomes stable earlier in the high-Co cermet. The strain threshold for the nucleation of dislocations in WC increases with Co content. At the later stage of deformation, the growth rate of WC dislocation density increases more rapidly in the cermet with lower Co content, which exhibits an opposite tendency compared with Co dislocation density. The relative contribution of Co and WC to the plasticity of the cermet varies in the deformation process. With a low Co content, the density of WC dislocations becomes higher than that of Co dislocations at larger strains, indicating that WC may contribute more than Co to the plasticity of the nanocrystalline cermet at the final deformation stage. The findings in the present work will be applicable to a large variety of ceramic–metal composite materials.

A series of Co2−xRuxMnSi (x = 0, 0.25, 0.5, 0.75, 1) Heusler compounds were successfully synthesized. The heat-treatment conditions were crucial to make the materials form a single phase with a Heusler structure. With increasing Ru content, the half-metallic gap, lattice parameters and magnetization are continuously adjustable in a wide range. The Co2−xRuxMnSi (x = 0, 0.25) compounds are rigorous half-metals and show a T3 dependence of resistance at low temperature. The Co2−xRuxMnSi (x = 0.5, 0.75, 1) Heusler compounds are the nearly half-metallic materials and show a semiconductive dependence of resistance at low temperature. The experimental magnetization is consistent with that in theory and follows the Slater–Pauling rule. The Curie temperature is higher than 750 K for all Co2−xRuxMnSi Heusler compounds.

The first ab initio aspherical structure refinement against experimental X-ray structure factors for polypeptides and proteins using a fragmentation approach to break up the protein into residues and solvent, thereby speeding up quantum-crystallographic Hirshfeld atom refinement (HAR) calculations, is described. It it found that the geometric and atomic displacement parameters from the new fragHAR method are essentially unchanged from a HAR on the complete unfragmented system when tested on dipeptides, tripeptides and hexapeptides. The largest changes are for the parameters describing H atoms involved in hydrogen-bond interactions, but it is shown that these discrepancies can be removed by including the interacting fragments as a single larger fragment in the fragmentation scheme. Significant speed-ups are observed for the larger systems. Using this approach, it is possible to perform a highly parallelized HAR in reasonable times for large systems. The method has been implemented in the TONTO software.

The non-steroidal anti-inflammatory drugs mefenamic acid (MFA) and tolfenamic acid (TFA) have a close resemblance in their molecular scaffold, whereby a methyl group in MFA is substituted by a chloro group in TFA. The present study demonstrates the isomorphous nature of these compounds in a series of their multicomponent solids. Furthermore, the unique nature of MFA and TFA has been demonstrated while excavating their alternate solid forms in that, by varying the drug (MFA or TFA) to coformer [4-di­methyl­amino­pyridine (DMAP)] stoichiometric ratio, both drugs have produced three different types of multicomponent crystals, viz. salt (1:1; API to coformer ratio), salt hydrate (1:1:1) and cocrystal salt (2:1). Interestingly, as anticipated from the close similarity of TFA and MFA structures, these multicomponent solids have shown an isomorphous relation. A thorough characterization and structural investigation of the new multicomponent forms of MFA and TFA revealed their similarity in terms of space group and structural packing with isomorphic nature among the pairs. Herein, the experimental results are generalized in a broader perspective for predictably identifying any possible new forms of comparable compounds by mapping their crystal structure landscapes. The utility of such an approach is evident from the identification of polymorph VI of TFA from hetero-seeding with isomorphous MFA form I from acetone–methanol (1:1) solution. That aside, a pseudopolymorph of TFA with di­methyl­formamide (DMF) was obtained, which also has some structural similarity to that of the solvate MFA:DMF. These new isostructural pairs are discussed in the context of solid form screening using structural landscape similarity.

3D electron diffraction (3DED) has been used to follow polymorph evolution in the crystallization of glycine from aqueous solution. The three polymorphs of glycine which exist under ambient conditions follow the stability order β < α < γ. The least stable β polymorph forms within the first 3 min, but this begins to yield the α-form after only 1 min more. Both structures could be determined from continuous rotation electron diffraction data collected in less than 20 s on crystals of thickness ∼100 nm. Even though the γ-form is thermodynamically the most stable polymorph, kinetics favour the α-form, which dominates after prolonged standing. In the same sample, some β and one crystallite of the γ polymorph were also observed.

Co-crystallization is a phenomenon widely employed to enhance the physio-chemical and biological properties of active pharmaceutical ingredients (APIs). Exemestane, or 6-methyl­ideneandrosta-1,4-diene-3,17-dione, is an anabolic steroid used as an irreversible steroidal aromatase inhibitor, which is in clinical use to treat breast cancer. The present study deals with the synthesis of co-crystals of exemestane with thio­urea by liquid-assisted grinding. The purity and homogeneity of the exemestane–thio­urea (1:1) co-crystal were confirmed by single-crystal X-ray diffraction followed by thermal stability analysis on the basis of differential scanning calorimetry and thermogravimetric analysis. Detailed geometric analysis of the co-crystal demonstrated that a 1:1 co-crystal stoichiometry is sustained by N—H⋯O hydrogen bonding between the amine (NH2) groups of thio­urea and the carbonyl group of exemestane. The synthesized co-crystal exhibited potent urease inhibition activity in vitro (IC50 = 3.86 ± 0.31 µg ml−1) compared with the API (exemestane), which was found to be inactive, and the co-former (thio­urea) (IC50 = 21.0 ± 1.25 µg ml−1), which is also an established tested standard for urease inhibition assays in vitro. The promising results of the present study highlight the significance of co-crystallization as a crystal engineering tool to improve the efficacy of pharmaceutical ingredients. Furthermore, the role of various hydrogen bonds in the crystal stability is successfully analysed quantitatively using Hirshfeld surface analysis.

With the increasing trend of using microcrystals and intense microbeams at synchrotron X-ray beamlines, radiation damage becomes a more pressing problem. Theoretical calculations show that the photoelectrons that primarily cause damage can escape microcrystals. This effect would become more pronounced with decreasing crystal size as well as at higher energies. To prove this effect, data from cryocooled lysozyme crystals of dimensions 5 × 3 × 3 and 20 × 8 × 8 µm mounted on cryo-transmission electron microscopy (cryo-TEM) grids were collected at 13.5 and 20.1 keV using a PILATUS CdTe 2M detector, which has a similar quantum efficiency at both energies. Accurate absorbed doses were calculated through the direct measurement of individual crystal sizes using scanning electron microscopy after the experiment and characterization of the X-ray microbeam. The crystal lifetime was then quantified based on the D1/2 metric. In this first systematic study, a longer crystal lifetime for smaller crystals was observed and crystal lifetime increased at higher X-ray energies, supporting the theoretical predictions of photoelectron escape. The use of detector technologies specifically optimized for data collection at energies above 20 keV allows the theoretically predicted photoelectron escape to be quantified and exploited, guiding future beamline-design choices.

Serial crystallography has enabled the study of complex biological questions through the determination of biomolecular structures at room temperature using low X-ray doses. Furthermore, it has enabled the study of protein dynamics by the capture of atomically resolved and time-resolved molecular movies. However, the study of many biologically relevant targets is still severely hindered by high sample consumption and lengthy data-collection times. By combining serial synchrotron crystallography (SSX) with 3D printing, a new experimental platform has been created that tackles these challenges. An affordable 3D-printed, X-ray-compatible microfluidic device (3D-MiXD) is reported that allows data to be collected from protein microcrystals in a 3D flow with very high hit and indexing rates, while keeping the sample consumption low. The miniaturized 3D-MiXD can be rapidly installed into virtually any synchrotron beamline with only minimal adjustments. This efficient collection scheme in combination with its mixing geometry paves the way for recording molecular movies at synchrotrons by mixing-triggered millisecond time-resolved SSX.

This paper recounts the first successful cryo-cooling of protein crystals that demonstrated the reduction in X-ray damage to macromolecular crystals. The project was suggested by David C. Phillips in 1965 at the Royal Institution of Great Britain and continued in 1967 at the Weizmann Institute of Science, where the first cryo-cooling experiments were performed on lysozyme crystals, and was completed in 1969 at Purdue University on lactate dehydrogenase crystals. A 1970 publication in Acta Crystallographica described the cryo-procedures, the use of cryo-protectants to prevent ice formation, the importance of fast, isotropic cryo-cooling and the collection of analytical data showing more than a tenfold decrease in radiation damage in cryo-cooled lactate dehydrogenase crystals. This was the first demonstration of any method that reduced radiation damage in protein crystals, which provided crystallographers with suitable means to employ synchrotron X-ray sources for protein-crystal analysis. Today, fifty years later, more than 90% of the crystal structures deposited in the Protein Data Bank have been cryo-cooled.

Bicontinuous cubic structures in soft matter consist of two intertwining labyrinths separated by a partitioning layer. Combining experiments, numerical modelling and techniques in differential geometry, we investigate twinning defects in bicontinuous cubic structures. We first demonstrate that a twin boundary is most likely to occur at a plane that cuts the partitioning layer almost perpendicularly, so that the perturbation caused by twinning remains minimal. This principle can be used as a criterion to identify potential twin boundaries, as demonstrated through detailed investigations of mesoporous silica crystals characterized by diamond and gyroid surfaces. We then discuss that a twin boundary can result from a stacking fault in the arrangement of inter-lamellar attachments at an early stage of structure formation. It is further shown that enhanced curvature fluctuations near the twin boundary would cost energy because of geometrical frustration, which would be eased by a crystal distortion that is experimentally observed.

Innovative new crystallographic methods are facilitating structural studies from ever smaller crystals of biological macromolecules. In particular, serial X-ray crystallography and microcrystal electron diffraction (MicroED) have emerged as useful methods for obtaining structural information from crystals on the nanometre to micrometre scale. Despite the utility of these methods, their implementation can often be difficult, as they present many challenges that are not encountered in traditional macromolecular crystallography experiments. Here, XFEL serial crystallography experiments and MicroED experiments using batch-grown microcrystals of the enzyme cyclophilin A are described. The results provide a roadmap for researchers hoping to design macromolecular microcrystallography experiments, and they highlight the strengths and weaknesses of the two methods. Specifically, we focus on how the different physical conditions imposed by the sample-preparation and delivery methods required for each type of experiment affect the crystal structure of the enzyme.

During single-crystal-to-single-crystal (SCSC) phase transitions, a polymorph of a compound can transform to a more stable form while remaining in the solid state. By understanding the mechanism of these transitions, strategies can be developed to control this phenomenon. This is particularly important in the pharmaceutical industry, but also relevant for other industries such as the food and agrochemical industries. Although extensive literature exists on SCSC phase transitions in inorganic crystals, it is unclear whether their classications and mechanisms translate to molecular crystals, with weaker interactions and more steric hindrance. A comparitive study of SCSC phase transitions in aliphatic linear-chain amino acid crystals, both racemates and quasi-racemates, is presented. A total of 34 transitions are considered and most are classified according to their structural change during the transition. Transitions without torsional changes show very different characteristics, such as transition temperature, enthalpy and free energy, compared with transitions that involve torsional changes. These differences can be rationalized using classical nucleation theory and in terms of a difference in mechanism; torsional changes occur in a molecule-by-molecule fashion, whereas transitions without torsional changes involve cooperative motion with multiple molecules at the same time.

Characterizing and controlling the uniformity of nanoparticles is crucial for their application in science and technology because crystalline defects in the nanoparticles strongly affect their unique properties. Recently, ultra-short and ultra-bright X-ray pulses provided by X-ray free-electron lasers (XFELs) opened up the possibility of structure determination of nanometre-scale matter with Å spatial resolution. However, it is often difficult to reconstruct the 3D structural information from single-shot X-ray diffraction patterns owing to the random orientation of the particles. This report proposes an analysis approach for characterizing defects in nanoparticles using wide-angle X-ray scattering (WAXS) data from free-flying single nanoparticles. The analysis method is based on the concept of correlated X-ray scattering, in which correlations of scattered X-ray are used to recover detailed structural information. WAXS experiments of xenon nanoparticles, or clusters, were conducted at an XFEL facility in Japan by using the SPring-8 Ångstrom compact free-electron laser (SACLA). Bragg spots in the recorded single-shot X-ray diffraction patterns showed clear angular correlations, which offered significant structural information on the nanoparticles. The experimental angular correlations were reproduced by numerical simulation in which kinematical theory of diffraction was combined with geometric calculations. We also explain the diffuse scattering intensity as being due to the stacking faults in the xenon clusters.

X-ray diffraction studies of crystals under pressure and quantitative experimental charge density analysis are among the most demanding types of crystallographic research. A successful feasibility study of the electron density in the mineral grossular under 1 GPa pressure conducted at the CRISTAL beamline at the SOLEIL synchrotron is presented in this work. A single crystal was placed in a diamond anvil cell, but owing to its special design (wide opening angle), short synchrotron wavelength and the high symmetry of the crystal, data with high completeness and high resolution were collected. This allowed refinement of a full multipole model of experimental electron distribution. Results are consistent with the benchmark measurement conducted without a diamond-anvil cell and also with the literature describing investigations of similar structures. Results of theoretical calculations of electron density distribution on the basis of dynamic structure factors mimic experimental findings very well. Such studies allow for laboratory simulations of processes which take place in the Earth's mantle.

Human septins 3, 9 and 12 are the only members of a specific subgroup of septins that display several unusual features, including the absence of a C-terminal coiled coil. This particular subgroup (the SEPT3 septins) are present in rod-like octameric protofilaments but are lacking in similar hexameric assemblies, which only contain representatives of the three remaining subgroups. Both hexamers and octamers can self-assemble into mixed filaments by end-to-end association, implying that the SEPT3 septins may facilitate polymerization but not necessarily function. These filaments frequently associate into higher order complexes which associate with biological membranes, triggering a wide range of cellular events. In the present work, a complete compendium of crystal structures for the GTP-binding domains of all of the SEPT3 subgroup members when bound to either GDP or to a GTP analogue is provided. The structures reveal a unique degree of plasticity at one of the filamentous interfaces (dubbed NC). Specifically, structures of the GDP and GTPγS complexes of SEPT9 reveal a squeezing mechanism at the NC interface which would expel a polybasic region from its binding site and render it free to interact with negatively charged membranes. On the other hand, a polyacidic region associated with helix α5', the orientation of which is particular to this subgroup, provides a safe haven for the polybasic region when retracted within the interface. Together, these results suggest a mechanism which couples GTP binding and hydrolysis to membrane association and implies a unique role for the SEPT3 subgroup in this process. These observations can be accounted for by constellations of specific amino-acid residues that are found only in this subgroup and by the absence of the C-terminal coiled coil. Such conclusions can only be reached owing to the completeness of the structural studies presented here.

Alongshan virus (ALSV) is an emerging human pathogen that was identified in China and rapidly spread to the European continent in 2019, raising concerns about public health. ALSV belongs to the distinct Jingmenvirus group within the Flaviviridae family with segmented RNA genomes. While segments 2 and 4 of the ALSV genome encode the VP1–VP3 proteins of unknown origin, segments 1 and 3 encode the NS2b–NS3 and NS5 proteins, which are related to Flavivirus nonstructural proteins, suggesting an evolutionary link between segmented and unsegmented viruses within the Flaviviridae family. Here, the enzymatic activity of the ALSV NS3-like helicase (NS3-Hel) was characterized and its crystal structure was determined to 2.9 Å resolution. ALSV NS3-Hel exhibits an ATPase activity that is comparable to those measured for Flavivirus NS3 helicases. The structure of ALSV NS3-Hel exhibits an overall fold similar to those of Flavivirus NS3 helicases. Despite the limited amino-acid sequence identity between ALSV NS3-Hel and Flavivirus NS3 helicases, structural features at the ATPase active site and the RNA-binding groove remain conserved in ALSV NS3-Hel. These findings provide a structural framework for drug design and suggest the possibility of developing a broad-spectrum antiviral drug against both Flavivirus and Jingmenvirus.

Developing methods to determine high-resolution structures from micrometre- or even submicrometre-sized protein crystals has become increasingly important in recent years. This applies to both large protein complexes and membrane proteins, where protein production and the subsequent growth of large homogeneous crystals is often challenging, and to samples which yield only micro- or nanocrystals such as amyloid or viral polyhedrin proteins. The versatile macromolecular crystallography microfocus (VMXm) beamline at Diamond Light Source specializes in X-ray diffraction measurements from micro- and nanocrystals. Because of the possibility of measuring data from crystalline samples that approach the resolution limit of visible-light microscopy, the beamline design includes a scanning electron microscope (SEM) to visualize, locate and accurately centre crystals for X-ray diffraction experiments. To ensure that scanning electron microscopy is an appropriate method for sample visualization, tests were carried out to assess the effect of SEM radiation on diffraction quality. Cytoplasmic polyhedrosis virus polyhedrin protein crystals cryocooled on electron-microscopy grids were exposed to SEM radiation before X-ray diffraction data were collected. After processing the data with DIALS, no statistically significant difference in data quality was found between datasets collected from crystals exposed and not exposed to SEM radiation. This study supports the use of an SEM as a tool for the visualization of protein crystals and as an integrated visualization tool on the VMXm beamline.

Nanocrystallography has transformed our ability to interrogate the atomic structures of proteins, peptides, organic molecules and materials. By probing atomic level details in ordered sub-10 nm regions of nanocrystals, scanning nanobeam electron diffraction extends the reach of nanocrystallography and in principle obviates the need for diffraction from large portions of one or more crystals. Scanning nanobeam electron diffraction is now applied to determine atomic structures from digitally defined regions of beam-sensitive peptide nanocrystals. Using a direct electron detector, thousands of sparse diffraction patterns over multiple orientations of a given crystal are recorded. Each pattern is assigned to a specific location on a single nanocrystal with axial, lateral and angular coordinates. This approach yields a collection of patterns that represent a tilt series across an angular wedge of reciprocal space: a scanning nanobeam diffraction tomogram. Using this diffraction tomogram, intensities can be digitally extracted from any desired region of a scan in real or diffraction space, exclusive of all other scanned points. Intensities from multiple regions of a crystal or from multiple crystals can be merged to increase data completeness and mitigate missing wedges. It is demonstrated that merged intensities from digitally defined regions of two crystals of a segment from the OsPYL/RCAR5 protein produce fragment-based ab initio solutions that can be refined to atomic resolution, analogous to structures determined by selected-area electron diffraction. In allowing atomic structures to now be determined from digitally outlined regions of a nanocrystal, scanning nanobeam diffraction tomography breaks new ground in nanocrystallography.

The structure of a decagonal quasicrystal in the Zn58Mg40Y2 (at.%) alloy was studied using electron diffraction and atomic resolution Z-contrast imaging techniques. This stable Frank–Kasper Zn–Mg–Y decagonal quasicrystal has an atomic structure which can be modeled with a rhombic/hexagonal tiling decorated with icosahedral units at each vertex. No perfect decagonal clusters were observed in the Zn–Mg–Y decagonal quasicrystal, which differs from the Zn–Mg–Dy decagonal crystal with the same space group P10/mmm. Y atoms occupy the center of `dented decagon' motifs consisting of three fat rhombic and two flattened hexagonal tiles. About 75% of fat rhombic tiles are arranged in groups of five forming star motifs, while the others connect with each other in a `zigzag' configuration. This decagonal quasicrystal has a composition of Zn68.3Mg29.1Y2.6 (at.%) with a valence electron concentration (e/a) of about 2.03, which is in accord with the Hume–Rothery criterion for the formation of the Zn-based quasicrystal phase (e/a = 2.0–2.15).

Nucleation of crystals from solution is fundamental to many natural and industrial processes. In this work, the molecular mechanism of conformational polymorphism nucleation and the links between the molecular conformation in solutions and in crystals were investigated in detail by using 5-nitro­furazone as the model compound. Different polymorphs were prepared, and the conformations in solutions obtained by dissolving different polymorphs were analysed and compared. The solutions of 5-nitro­furazone were proven to contain multiple conformers through quantum chemical computation, Raman spectra analysis, 2D nuclear Overhauser effect spectroscopy spectra analysis and molecular dynamics simulation. The conformational evolution and desolvation path was illustrated according to the 1H NMR spectra of solutions with different concentrations. Finally, based on all the above analysis, the molecular conformational evolution path during nucleation of 5-nitro­furazone was illustrated. The results presented in this work shed a new light on the molecular mechanism of conformational polymorphism nucleation in solution.

The performance of automated model building in crystal structure determination usually decreases with the resolution of the experimental data, and may result in fragmented models and incorrect side-chain assignment. Presented here are new methods for machine-learning-based docking of main-chain fragments to the sequence and for their sequence-independent connection using a dedicated library of protein fragments. The combined use of these new methods noticeably increases sequence coverage and reduces fragmentation of the protein models automatically built with ARP/wARP.

Olfactomedins are a family of modular proteins found in multicellular organisms that all contain five-bladed β-propeller olfactomedin (OLF) domains. In support of differential functions for the OLF propeller, the available crystal structures reveal that only some OLF domains harbor an internal calcium-binding site with ligands derived from a triad of residues. For the myocilin OLF domain (myoc-OLF), ablation of the ion-binding site (triad Asp, Asn, Asp) by altering the coordinating residues affects the stability and overall structure, in one case leading to misfolding and glaucoma. Bioinformatics analysis reveals a variety of triads with possible ion-binding characteristics lurking in OLF domains in invertebrate chordates such as Arthropoda (Asp–Glu–Ser), Nematoda (Asp–Asp–His) and Echinodermata (Asp–Glu–Lys). To test ion binding and to extend the observed connection between ion binding and distal structural rearrangements, consensus triads from these phyla were installed in the myoc-OLF. All three protein variants exhibit wild-type-like or better stability, but their calcium-binding properties differ, concomitant with new structural deviations from wild-type myoc-OLF. Taken together, the results indicate that calcium binding is not intrinsically destabilizing to myoc-OLF or required to observe a well ordered side helix, and that ion binding is a differential feature that may underlie the largely elusive biological function of OLF propellers.

The bacterial flavin adenine dinucleotide (FAD)-dependent glucose dehydrogenase complex derived from Burkholderia cepacia (BcGDH) is a representative molecule of direct electron transfer-type FAD-dependent dehydrogenase complexes. In this study, the X-ray structure of BcGDHγα, the catalytic subunit (α-subunit) of BcGDH complexed with a hitchhiker protein (γ-subunit), was determined. The most prominent feature of this enzyme is the presence of the 3Fe–4S cluster, which is located at the surface of the catalytic subunit and functions in intramolecular and intermolecular electron transfer from FAD to the electron-transfer subunit. The structure of the complex revealed that these two molecules are connected through disulfide bonds and hydrophobic interactions, and that the formation of disulfide bonds is required to stabilize the catalytic subunit. The structure of the complex revealed the putative position of the electron-transfer subunit. A comparison of the structures of BcGDHγα and membrane-bound fumarate reductases suggested that the whole BcGDH complex, which also includes the membrane-bound β-subunit containing three heme c moieties, may form a similar overall structure to fumarate reductases, thus accomplishing effective electron transfer.

Theoretically, crystals with supercells exist at a unique crossroads where they can be considered as either a large unit cell with closely spaced reflections in reciprocal space or a higher dimensional superspace with a modulation that is commensurate with the supercell. In the latter case, the structure would be defined as an average structure with functions representing a modulation to determine the atomic location in 3D space. Here, a model protein structure and simulated diffraction data were used to investigate the possibility of solving a real incommensurately modulated protein crystal using a supercell approximation. In this way, the answer was known and the refinement method could be tested. Firstly, an average structure was solved by using the `main' reflections, which represent the subset of the reflections that belong to the subcell and in general are more intense than the `satellite' reflections. The average structure was then expanded to create a supercell and refined using all of the reflections. Surprisingly, the refined solution did not match the expected solution, even though the statistics were excellent. Interestingly, the corresponding superspace group had multiple 3D daughter supercell space groups as possibilities, and it was one of the alternate daughter space groups that the refinement locked in on. The lessons learned here will be applied to a real incommensurately modulated profilin–actin crystal that has the same superspace group.

For the extraction of the best possible X-ray diffraction data from macromolecular crystals, accurate positioning of the crystals with respect to the X-ray beam is crucial. In addition, information about the shape and internal defects of crystals allows the optimization of data-collection strategies. Here, it is demonstrated that the X-ray beam available on the macromolecular crystallo­graphy beamline P14 at the high-brilliance synchrotron-radiation source PETRA III at DESY, Hamburg, Germany can be used for high-energy phase-contrast microtomography of protein crystals mounted in an optically opaque lipidic cubic phase matrix. Three-dimensional tomograms have been obtained at X-ray doses that are substantially smaller and on time scales that are substantially shorter than those used for diffraction-scanning approaches that display protein crystals at micrometre resolution. Adding a compound refractive lens as an objective to the imaging setup, two-dimensional imaging at sub-micrometre resolution has been achieved. All experiments were performed on a standard macromolecular crystallography beamline and are compatible with standard diffraction data-collection workflows and apparatus. Phase-contrast X-ray imaging of macromolecular crystals could find wide application at existing and upcoming low-emittance synchrotron-radiation sources.

Serial crystallography is having an increasing impact on structural biology. This emerging technique opens up new possibilities for studying protein structures at room temperature and investigating structural dynamics using time-resolved X-ray diffraction. A limitation of the method is the intrinsic need for large quantities of well ordered micrometre-sized crystals. Here, a method is presented to screen for conditions that produce microcrystals of membrane proteins in the lipidic cubic phase using a well-based crystallization approach. A key advantage over earlier approaches is that the progress of crystal formation can be easily monitored without interrupting the crystallization process. In addition, the protocol can be scaled up to efficiently produce large quantities of crystals for serial crystallography experiments. Using the well-based crystallization methodology, novel conditions for the growth of showers of microcrystals of three different membrane proteins have been developed. Diffraction data are also presented from the first user serial crystallography experiment performed at MAX IV Laboratory.

As part of the Virus-X Consortium that aims to identify and characterize novel proteins and enzymes from bacteriophages and archaeal viruses, the genes of the putative lytic proteins XepA from Bacillus subtilis prophage PBSX and YomS from prophage SPβ were cloned and the proteins were subsequently produced and functionally characterized. In order to elucidate the role and the molecular mechanism of XepA and YomS, the crystal structures of these proteins were solved at resolutions of 1.9 and 1.3 Å, respectively. XepA consists of two antiparallel β-sandwich domains connected by a 30-amino-acid linker region. A pentamer of this protein adopts a unique dumbbell-shaped architecture consisting of two discs and a central tunnel. YomS (12.9 kDa per monomer), which is less than half the size of XepA (30.3 kDa), shows homology to the C-terminal part of XepA and exhibits a similar pentameric disc arrangement. Each β-sandwich entity resembles the fold of typical cytoplasmic membrane-binding C2 domains. Only XepA exhibits distinct cytotoxic activity in vivo, suggesting that the N-terminal pentameric domain is essential for this biological activity. The biological and structural data presented here suggest that XepA disrupts the proton motive force of the cytoplasmatic membrane, thus supporting cell lysis.

The refinement of biomolecular crystallographic models relies on geometric restraints to help to address the paucity of experimental data typical in these experiments. Limitations in these restraints can degrade the quality of the resulting atomic models. Here, an integration of the full all-atom Amber molecular-dynamics force field into Phenix crystallographic refinement is presented, which enables more complete modeling of biomolecular chemistry. The advantages of the force field include a carefully derived set of torsion-angle potentials, an extensive and flexible set of atom types, Lennard–Jones treatment of nonbonded interactions and a full treatment of crystalline electrostatics. The new combined method was tested against conventional geometry restraints for over 22 000 protein structures. Structures refined with the new method show substantially improved model quality. On average, Ramachandran and rotamer scores are somewhat better, clashscores and MolProbity scores are significantly improved, and the modeling of electrostatics leads to structures that exhibit more, and more correct, hydrogen bonds than those refined using traditional geometry restraints. In general it is found that model improvements are greatest at lower resolutions, prompting plans to add the Amber target function to real-space refinement for use in electron cryo-microscopy. This work opens the door to the future development of more advanced applications such as Amber-based ensemble refinement, quantum-mechanical representation of active sites and improved geometric restraints for simulated annealing.

Although the success of molecular-replacement techniques requires the solution of a six-dimensional problem, this is often subdivided into two three-dimensional problems. REMO09 is one of the programs which have adopted this approach. It has been revisited in the light of a new probabilistic approach which is able to directly derive conditional distribution functions without passing through a previous calculation of the joint probability distributions. The conditional distributions take into account various types of prior information: in the rotation step the prior information may concern a non-oriented model molecule alone or together with one or more located model molecules. The formulae thus obtained are used to derive figures of merit for recognizing the correct orientation in the rotation step and the correct location in the translation step. The phases obtained by this new version of REMO09 are used as a starting point for a pipeline which in its first step extends and refines the molecular-replacement phases, and in its second step creates the final electron-density map which is automatically interpreted by CAB, an automatic model-building program for proteins and DNA/RNA structures.

Corrections are published for the article by Caldararu et al. [(2019), Acta Cryst. D75, 368–380].

Noncrystallographic symmetry (NCS) averaging following molecular-replacement phasing is generally the major technique used to solve a structure with several molecules in one asymmetric unit, such as a spherical icosahedral viral particle. As an alternative method to NCS averaging, a new approach to optimize or to refine the electron density directly under NCS constraints is proposed. This method has the same effect as the conventional NCS-averaging method but does not include the process of Fourier synthesis to generate the electron density from amplitudes and the corresponding phases. It has great merit for the solution of structures with limited data that are either twinned or incomplete at low resolution. This method was applied to the case of the T = 1 shell-domain subviral particle of Penaeus vannamei nodavirus with data affected by twinning using the REFMAC5 refinement software.