The crystal structures of 16 boron sub­phthalo­cy­an­ines (BsubPcs) with structurally diverse axial groups were analyzed and com­pared to elucidate the impact of the axial group on the inter­molecular π–π inter­actions, axial-group inter­actions, axial bond length and BsubPc bowl depth. π–π inter­actions between the iso­indole units of adjacent BsubPc mol­ecules most often involve concave–concave packing, whereas axial-group inter­actions with adjacent BsubPc mol­ecules tend to favour the convex side of the BsubPc bowl. Furthermore, axial groups that contain O and/or F atoms tend to have significant hy­dro­gen-bonding inter­actions, while axial groups containing arene site(s) can participate in π–π inter­actions with the BsubPc bowl, both of which can strongly influence the crystal packing. Bulky axial groups did tend to disrupt the π–π inter­actions and/or axial-group inter­actions, preventing some of the close packing that is seen in BsubPcs with less bulky axial groups. The atomic radius of the heteroatom bonded to boron directly influences the axial bond length, whereas the axial group has minimal impact on the BsubPc bowl depth. Finally, the crystal growth method did not generally appear to have a significant impact on the solid-state arrangement, with the exception of water occasionally being incorporated into crystal structures when hygroscopic solvents were used. These insights can help with the design and fine-tuning of the solid-state structures of BsubPcs as they continue to be developed as functional materials in organic electronics.

The incommensurately modulated structure of (2S,3S)-2-amino-3-hy­droxy-3-methyl-4-phen­oxy­butanoic acid dihydrate (C11H15NO4·2H2O or I·2H2O) is described in the (3+1)-dimensional superspace group P212121(0β0)000 (β = 0.357). The loss of the three-dimensional periodicity is ascribed to the occupational modulation of one positionally disordered solvent water mol­ecule, where the two positions are related by a small translation [ca 0.666 (9) Å] and ∼168 (5)° rotation about one of its O—H bonds, with an average 0.624 (3):0.376 (3) occupancy ratio. The occupational modulation of this mol­ecule arises due to the com­petition between the different hy­dro­gen-bonding motifs associated with each position. The structure can be very well refined in the average approximation (all satellite reflections disregarded) in the space group P212121, with the water mol­ecule refined as disordered over two positions in a 0.625 (16):0.375 (16) ratio. The refinement in the commensurate threefold supercell approximation in the space group P1121 is also of high quality, with the six corresponding water mol­ecules exhibiting three different occupancy ratios averaging 0.635:0.365.

It is well known that Hirshfeld surfaces provide an easy and straightforward way of analysing inter­molecular inter­actions in the crystal environment. The use of atomic Hirshfeld surfaces has also demonstrated that such surfaces carry information related to chemical bonds which allow a deeper evaluation of the structures. Here we briefly summarize the approach of atomic Hirshfeld surfaces while further evaluating the kind of information that can be retrieved from them. We show that the analysis of the metal-centre Hirshfeld surfaces from structures refined via Hirshfeld Atom Refinement (HAR) allow accurate evaluation of contacts of type M⋯H, and that such contacts can be related to the overall shape of the surfaces. The com­pounds analysed were tetra­aqua­bis­(3-carb­oxy­propionato)metal(II), [M(C4H3O4)2(H2O)4], for metal(II)/M = manganese/Mn, cobalt/Co, nickel/Ni and zinc/Zn. We also evaluate the sensitivity of the surfaces by an investigation of seemingly flat surfaces through analysis of the curvature functions in the direction of C—C bonds. The obtained values not only demonstrate variations in curvature but also show a correlation with the hybridization of the C atoms involved in the bond.

The activation of C—C bonds by transition-metal com­plexes is of continuing inter­est and aceto­nitrile (MeCN) has attracted attention as a cyanide source with com­paratively low toxicity for organic cyanation reactions. A di­iron end-on μ-η1:η1-CN-bridged com­plex was obtained from a crystallization experiment of an open-chain iron–NHC com­plex, namely, μ-cyanido-κ2C:N-bis­{[(aceto­nitrile-κN)[3,3'-bis­(pyridin-2-yl)-1,1'-(methyl­idene)bis­(benzimidazol-2-yl­idene)]iron(II)} tris­(hexa­fluoro­phos­phate), [Fe2(CN)(C2H3N)2(C25H18N6)2](PF6)3. The cyanide appears to originate from the MeCN solvent by C—C bond cleavage or through carbon–hy­dro­gen oxidation.

Herein we report the crystal structures of two ben­zo­di­az­e­pines obtained by reacting N,N'-(4,5-di­amino-1,2-phenyl­ene)bis­(4-methyl­ben­zene­sul­fon­am­ide) (1) or 4,5-(4-methyl­ben­zene­sul­fon­am­ido)­ben­zene-1,2-diaminium dichloride (1·2HCl) with acetone, giving 2,2,4-trimethyl-8,9-bis­(4-methyl­ben­zene­sul­fon­am­ido)-2,3-di­hydro-5H-1,5-ben­zo­di­az­e­pine, C26H30N4O4S2 (2), and 2,2,4-tri­methyl-8,9-bis­(4-methyl­ben­zene­sul­fon­am­ido)-2,3-di­hydro-5H-1,5-ben­zo­di­az­e­pin-1-ium chloride 0.3-hydrate, C26H31N4O4S2+·Cl−·0.3H2O (3). Compounds 2 and 3 were first obtained in attempts to recrystallize 1 and 1·2HCl using acetone as solvent. This solvent reacted with the vicinal di­amines present in the mol­ecular structures, forming a 5H-1,5-ben­zo­di­az­e­pine ring. In the crystal structure of 2, the seven-membered ring of ben­zo­di­az­e­pine adopts a boat-like conformation, while upon protonation, observed in the crystal structure of 3, it adopts an envelope-like conformation. In both crystalline com­pounds, the tosyl­amide N atoms are not in resonance with the arene ring, mainly due to hy­dro­gen bonds and steric hindrance caused by the large vicinal groups in the aromatic ring. At a supra­molecular level, the crystal structure is maintained by a combination of hy­dro­gen bonds and hydro­phobic inter­actions. In 2, amine-to-tosyl N—H⋯O and amide-to-imine N—H⋯N hy­dro­gen bonds can be observed. In contrast, in 3, the chloride counter-ion and water mol­ecule result in most of the hy­dro­gen bonds being of the amide-to-chloride and ammonium-to-chloride N—H⋯Cl types, while the amine inter­acts with the tosyl group, as seen in 2. In conclusion, we report the synthesis of 1, 1·2HCl and 2, as well as their chemical characterization. For 2, two synthetic methods are described, i.e. solvent-mediated crystallization and synthesis via a more efficient and cleaner route as a polycrystalline material. Salt 3 was only obtained as presented, with only a few crystals being formed.

Treating the amide acetanilide (N-phenyl­acetamide, C8H9NO) with aqueous strong acids allowed the structures of five hemi-protonated salt forms of acetanilide to be elucidated. N-(1-Hy­droxy­ethyl­idene)anilinium chloride–N-phenyl­acetamide (1/1), [(C8H9NO)2H][Cl], and the bromide, [(C8H9NO)2H][Br], triiodide, [(C8H9NO)2H][I3], tetra­fluoro­borate, [(C8H9NO)2H][BF4], and di­iodo­bromide hemi(diiodine), [(C8H9NO)2H][I2Br]·0.5I2, analogues all feature centrosymmetric dimeric units linked by O—H⋯O hy­dro­gen bonds that extend into one-dimensional hy­dro­gen-bonded chains through N—H⋯X inter­actions, where X is the halide atom of the anion. Protonation occurs at the amide O atom and results in systematic lengthening of the C=O bond and a corresponding shortening of the C—N bond. The size of these geometric changes is similar to those found for hemi-protonated paracetamol structures, but less than those in fully protonated paracetamol structures. The bond angles of the amide fragments are also found to change on protonation, but these angular changes are also influenced by conformation, namely, whether the amide group is coplanar with the phenyl ring or twisted out of plane.

Three structurally diverse 5-phenyl­tetra­zolato (Tz) Ti, Zr, and Ta com­plexes, namely, (C2H8N)[Ti2(C7H5N4)5(C2H6N)4]·1.45C6H6 or (Me2NH2)[Ti2(NMe2)4(2,3-μ-Tz)3(2-η1-Tz)2]·1.45C6H6, (1·1.45C6H6), [Zr2(C7H5N4)6(C2H6N)2(C2H7N)2]·1.12C6H6·0.382CH2Cl2 or [Zr2(Me2NH)2(NMe2)2(2,3-μ-Tz)3(2-η1-Tz)2(1,2-η2-Tz)]·1.12C6H6·0.38CH2Cl2 (2·1.12C6H6·0.38CH2Cl2), and (C2H8N)2[Ta2(C7H5N4)8(C2H6N)2O]·0.25C7H8 or (Me2NH2)2[Ta2(NMe2)2(2,3-μ-Tz)2(2-η1-Tz)6O]·0.25C7H8 (3·0.25C7H8), where TzH is 5-phenyl-1H-tetra­zole, have been synthesized and structurally characterized. All three com­plexes are dinuclear; the Ti center in 1 is six-coordinate, whereas the Zr and Ta atoms in 2 and 3 are seven-coordinate. The coordination environments of the Ti centers in 1 are similar, and so are the ligations of the Ta centers in 3. In contrast, the two Zr centers in 2 bear a different number of ligands, one of which is a bidentate η2-5-phenyl­tetra­zolato ligand that has not been observed previously for d-block elements. The di­methyl­amido ligand, present in the starting materials, remained un­changed, or was converted to di­methyl­amine and di­methyl­ammonium during the synthesis. Di­methyl­amine coordinates as a neutral ligand, whereas di­methyl­ammonium is retained as a hy­dro­gen-bonded entity bridging Tz ligands.

In this study, we report the results of continuous rotation electron diffraction studies of single DyPO4·nH2O (rhabdophane) nanocrystals. The diffraction patterns can be fit to a trigonal lattice (P3121) with lattice parameters a = 7.019 (5) and c = 6.417 (5) Å. However, there is also a set of diffuse background scattering features present that are associated with a disordered superstructure that is double these lattice parameters and fits with an arrangement of water mol­ecules present in the structure pore. Pair distribution function (PDF) maps based on the diffuse background allowed the extent of the water correlation to be estimated, with 2–3 nm correlation along the c axis and ∼5 nm along the a/b axis.

Historically, knowledge of the mol­ecular packing within the crystal structures of organic semi­con­duc­tors has been instrumental in understanding their solid-state electronic properties. Nowadays, crystal structures are thus becoming increasingly important for enabling engineering properties, understanding poly­mor­phism in bulk and in thin films, exploring dynamics and elucidating phase-transition mech­a­nisms. This review article introduces the most salient and recent results of the field.

Starting from (p-tolyl­sulfin­yl)fer­ro­cene (1), a mixture of the complete series [CpFe{C5H5–n(SOTol-p)n}] (n = 2–4) (2–4) in all regioisomers was obtained. After chromatographic separation, crystals of 1,2-bis­[(4-methyl­benzene)­sulfin­yl]fer­ro­cene, 2a, and 1,3-bis­[(4-methyl­benzene)­sulfin­yl]fer­ro­cene, 2b, both [Fe(C5H5)(C19H17O2S2)], as well as of 1,2,3-tris­[(4-methyl­benzene)­sulfin­yl]fer­ro­cene, [Fe(C5H5)(C26H23O3S3)], 3a, and 1,2,3,4-tetra­kis­[(4-methyl­benzene)­sul­fin­yl]fer­ro­cene ethyl acetate 0.75-solvate, [Fe(C5H5)(C33H29O4S4)]·0.75C4H8O2, 4, could be isolated. Their mol­ecular and crystal structures are compared with each other and also with the so far un­reported structures of related 1,2-bis­(phenyl­sulfan­yl)fer­ro­cene, [Fe(C5H5)(C17H13S2)], 5, and 1,2,3,4-tetra­kis­(phenyl­sulfan­yl)fer­ro­cene, [Fe(C5H5)(C29H21S4)], 6. In all the sulfinyl structures, the O atoms of the S=O groups are in equatorial positions, except for that in tetrasubstituted 4. All the arene rings of these com­pounds (except for one ring in 4) are in axial positions directed away from the Fe atom, mostly in a near perpendicular orientation with respect to the plane of the cyclo­penta­di­en­yl ring. The main inter­molecular inter­actions in the crystals are C—H⋯H—C, C—H⋯π and C—H⋯O, while C—H⋯S inter­actions are much less important, except for tetra­sul­fan­yl com­pound 6. π–π inter­actions (intra­molecular) are only important in com­pound 3a. Hirshfeld analysis shows that dispersion terms are dominant for the inter­action energies of all six com­pounds. In general, the calculated total inter­action energies increase with increasing number of substituents and are higher for the sulfinyl than for the sul­fan­yl groups.

Catalase is an antioxidant enzyme that breaks down hydrogen peroxide (H2O2) into molecular oxygen and water. In all monofunctional catalases the pathway that H2O2 takes to the catalytic centre is via the `main channel'. However, the structure of this channel differs in large-subunit and small-subunit catalases. In large-subunit catalases the channel is 15 Å longer and consists of two distinct parts, including a hydrophobic lower region near the heme and a hydrophilic upper region where multiple H2O2 routes are possible. Conserved glutamic acid and threonine residues are located near the intersection of these two regions. Mutations of these two residues in the Scytalidium thermophilum catalase had no significant effect on catalase activity. However, the secondary phenol oxidase activity was markedly altered, with kcat and kcat/Km values that were significantly increased in the five variants E484A, E484I, T188D, T188I and T188F. These variants also showed a lower affinity for inhibitors of oxidase activity than the wild-type enzyme and a higher affinity for phenolic substrates. Oxidation of heme b to heme d did not occur in most of the studied variants. Structural changes in solvent-chain integrity and channel architecture were also observed. In summary, modification of the main-channel gate glutamic acid and threonine residues has a greater influence on the secondary activity of the catalase enzyme, and the oxidation of heme b to heme d is predominantly inhibited by their conversion to aliphatic and aromatic residues.

Phenuiviridae nucleoprotein is the main structural and functional component of the viral cycle, protecting the viral RNA and mediating the essential replication/transcription processes. The nucleoprotein (N) binds the RNA using its globular core and polymerizes through the N-terminus, which is presented as a highly flexible arm, as demonstrated in this article. The nucleoprotein exists in an `open' or a `closed' conformation. In the case of the closed conformation the flexible N-terminal arm folds over the RNA-binding cleft, preventing RNA adsorption. In the open conformation the arm is extended in such a way that both RNA adsorption and N polymerization are possible. In this article, single-crystal X-ray diffraction and small-angle X-ray scattering were used to study the N protein of Toscana virus complexed with a single-chain camelid antibody (VHH) and it is shown that in the presence of the antibody the nucleoprotein is unable to achieve a functional assembly to form a ribonucleoprotein complex.

To identify starting points for therapeutics targeting SARS-CoV-2, the Paul Scherrer Institute and Idorsia decided to collaboratively perform an X-ray crystallographic fragment screen against its main protease. Fragment-based screening was carried out using crystals with a pronounced open conformation of the substrate-binding pocket. Of 631 soaked fragments, a total of 29 hits bound either in the active site (24 hits), a remote binding pocket (three hits) or at crystal-packing interfaces (two hits). Notably, two fragments with a pose that was sterically incompatible with a more occluded crystal form were identified. Two isatin-based electrophilic fragments bound covalently to the catalytic cysteine residue. The structures also revealed a surprisingly strong influence of the crystal form on the binding pose of three published fragments used as positive controls, with implications for fragment screening by crystallography.

Macromolecular crystallography generally requires the recovery of missing phase information from diffraction data to reconstruct an electron-density map of the crystallized molecule. Most recent structures have been solved using molecular replacement as a phasing method, requiring an a priori structure that is closely related to the target protein to serve as a search model; when no such search model exists, molecular replacement is not possible. New advances in computational machine-learning methods, however, have resulted in major advances in protein structure predictions from sequence information. Methods that generate predicted structural models of sufficient accuracy provide a powerful approach to molecular replacement. Taking advantage of these advances, AlphaFold predictions were applied to enable structure determination of a bacterial protein of unknown function (UniProtKB Q63NT7, NCBI locus BPSS0212) based on diffraction data that had evaded phasing attempts using MIR and anomalous scattering methods. Using both X-ray and micro-electron (microED) diffraction data, it was possible to solve the structure of the main fragment of the protein using a predicted model of that domain as a starting point. The use of predicted structural models importantly expands the promise of electron diffraction, where structure determination relies critically on molecular replacement.

Complex I (proton-pumping NADH:ubiquinone oxidoreductase) is the first component of the mitochondrial respiratory chain. In recent years, high-resolution cryo-EM studies of complex I from various species have greatly enhanced the understanding of the structure and function of this important membrane-protein complex. Less well studied is the structural basis of complex I biogenesis. The assembly of this complex of more than 40 subunits, encoded by nuclear or mitochondrial DNA, is an intricate process that requires at least 20 different assembly factors in humans. These are proteins that are transiently associated with building blocks of the complex and are involved in the assembly process, but are not part of mature complex I. Although the assembly pathways have been studied extensively, there is limited information on the structure and molecular function of the assembly factors. Here, the insights that have been gained into the assembly process using cryo-EM are reviewed.

Low-molecular-weight (LMW) thiols are involved in many processes in all organisms, playing a protective role against reactive species, heavy metals, toxins and antibiotics. Actinobacteria, such as Mycobacterium tuberculosis, use the LMW thiol mycothiol (MSH) to buffer the intracellular redox environment. The NADPH-dependent FAD-containing oxidoreductase mycothiol disulfide reductase (Mtr) is known to reduce oxidized mycothiol disulfide (MSSM) to MSH, which is crucial to maintain the cellular redox balance. In this work, the first crystal structures of Mtr are presented, expanding the structural knowledge and understanding of LMW thiol reductases. The structural analyses and docking calculations provide insight into the nature of Mtrs, with regard to the binding and reduction of the MSSM substrate, in the context of related oxidoreductases. The putative binding site for MSSM suggests a similar binding to that described for the homologous glutathione reductase and its respective substrate glutathione disulfide, but with distinct structural differences shaped to fit the bulkier MSSM substrate, assigning Mtrs as uniquely functioning reductases. As MSH has been acknowledged as an attractive antitubercular target, the structural findings presented in this work may contribute towards future antituberculosis drug development.

Mevalonate kinase is central to the isoprenoid biosynthesis pathway. Here, high-resolution X-ray crystal structures of two mevalonate kinases are presented: a eukaryotic protein from Ramazzottius varieornatus and an archaeal protein from Methanococcoides burtonii. Both enzymes possess the highly conserved motifs of the GHMP enzyme superfamily, with notable differences between the two enzymes in the N-terminal part of the structures. Biochemical characterization of the two enzymes revealed major differences in their sensitivity to geranyl pyrophosphate and farnesyl pyrophosphate, and in their thermal stabilities. This work adds to the understanding of the structural basis of enzyme inhibition and thermostability in mevalonate kinases.

The validation of structural models obtained by macromolecular X-ray crystallography against experimental diffraction data, whether before deposition into the PDB or after, is typically carried out exclusively against the merged data that are eventually archived along with the atomic coordinates. It is shown here that the availability of unmerged reflection data enables valuable additional analyses to be performed that yield improvements in the final models, and tools are presented to implement them, together with examples of the results to which they give access. The first example is the automatic identification and removal of image ranges affected by loss of crystal centering or by excessive decay of the diffraction pattern as a result of radiation damage. The second example is the `reflection-auditing' process, whereby individual merged data items showing especially poor agreement with model predictions during refinement are investigated thanks to the specific metadata (such as image number and detector position) that are available for the corresponding unmerged data, potentially revealing previously undiagnosed instrumental, experimental or processing problems. The third example is the calculation of so-called F(early) − F(late) maps from carefully selected subsets of unmerged amplitude data, which can not only highlight the location and extent of radiation damage but can also provide guidance towards suitable fine-grained parametrizations to model the localized effects of such damage.

The axoneme, a microtubule-based array at the center of every cilium, has been the subject of structural investigations for decades, but only recent advances in cryo-EM and cryo-ET have allowed a molecular-level interpretation of the entire complex to be achieved. The unique properties of the nine doublet microtubules and central pair of singlet microtubules that form the axoneme, including the highly decorated tubulin lattice and the docking of massive axonemal complexes, provide opportunities and challenges for sample preparation, 3D reconstruction and atomic modeling. Here, the approaches used for cryo-EM and cryo-ET of axonemes are reviewed, while highlighting the unique opportunities provided by the latest generation of AI-guided tools that are transforming structural biology.

A considerable bottleneck in serial crystallography at XFEL and synchrotron sources is the efficient production of large quantities of homogenous, well diffracting microcrystals. Efficient high-throughput screening of batch-grown microcrystals and the determination of ground-state structures from different conditions is thus of considerable value in the early stages of a project. Here, a highly sample-efficient methodology to measure serial crystallography data from microcrystals by raster scanning within standard in situ 96-well crystallization plates is described. Structures were determined from very small quantities of microcrystal suspension and the results were compared with those from other sample-delivery methods. The analysis of a two-dimensional batch crystallization screen using this method is also described as a useful guide for further optimization and the selection of appropriate conditions for scaling up microcrystallization.

Lanthanide ions have ideal chemical properties for catalysis, such as hard Lewis acidity, fast ligand-exchange kinetics, high coordination-number preferences and low geometric requirements for coordination. As a result, many small-molecule lanthanide catalysts have been described in the literature. Yet, despite the ability of enzymes to catalyse highly stereoselective reactions under gentle conditions, very few lanthanoenzymes have been investigated. In this work, the mononuclear binding of europium(III) and gadolinium(III) to the active site of a mutant of the model enzyme phosphotriesterase are described using X-ray crystallography at 1.78 and 1.61 Å resolution, respectively. It is also shown that despite coordinating a single non-natural metal cation, the PTE-R18 mutant is still able to maintain esterase activity.

Metalloproteins are ubiquitous in all living organisms and take part in a very wide range of biological processes. For this reason, their experimental characterization is crucial to obtain improved knowledge of their structure and biological functions. The three-dimensional structure represents highly relevant information since it provides insight into the interaction between the metal ion(s) and the protein fold. Such interactions determine the chemical reactivity of the bound metal. The available PDB structures can contain errors due to experimental factors such as poor resolution and radiation damage. A lack of use of distance restraints during the refinement and validation process also impacts the structure quality. Here, the aim was to obtain a thorough overview of the distribution of the distances between metal ions and their donor atoms through the statistical analysis of a data set based on more than 115 000 metal-binding sites in proteins. This analysis not only produced reference data that can be used by experimentalists to support the structure-determination process, for example as refinement restraints, but also resulted in an improved insight into how protein coordination occurs for different metals and the nature of their binding interactions. In particular, the features of carboxylate coordination were inspected, which is the only type of interaction that is commonly present for nearly all metals.

CdaA is the most widespread diadenylate cyclase in many bacterial species, including several multidrug-resistant human pathogens. The enzymatic product of CdaA, cyclic di-AMP, is a secondary messenger that is essential for the viability of many bacteria. Its absence in humans makes CdaA a very promising and attractive target for the development of new antibiotics. Here, the structural results are presented of a crystallographic fragment screen against CdaA from Listeria monocytogenes, a saprophytic Gram-positive bacterium and an opportunistic food-borne pathogen that can cause listeriosis in humans and animals. Two of the eight fragment molecules reported here were localized in the highly conserved ATP-binding site. These fragments could serve as potential starting points for the development of antibiotics against several CdaA-dependent bacterial species.

Type VII secretion (T7S) systems, also referred to as ESAT-6 secretion (ESX) systems, are molecular machines that have gained great attention due to their implications in cell homeostasis and in host–pathogen interactions in mycobacteria. The latter include important human pathogens such as Mycobacterium tuberculosis (Mtb), the etiological cause of human tuberculosis, which constitutes a pandemic accounting for more than one million deaths every year. The ESX-5 system is exclusively found in slow-growing pathogenic mycobacteria, where it mediates the secretion of a large family of virulence factors: the PE and PPE proteins. The secretion driving force is provided by EccC5, a multidomain ATPase that operates using four globular cytosolic domains: an N-terminal domain of unknown function (EccC5DUF) and three FtsK/SpoIIIE ATPase domains. Recent structural and functional studies of ESX-3 and ESX-5 systems have revealed EccCDUF to be an ATPase-like fold domain with potential ATPase activity, the functionality of which is essential for secretion. Here, the crystal structure of the MtbEccC5DUF domain is reported at 2.05 Å resolution, which reveals a nucleotide-free structure with degenerated cis-acting and trans-acting elements involved in ATP binding and hydrolysis. This crystallographic study, together with a biophysical assessment of the interaction of MtbEccC5DUF with ATP/Mg2+, supports the absence of ATPase activity proposed for this domain. It is shown that this degeneration is also present in DUF domains from other ESX and ESX-like systems, which are likely to exhibit poor or null ATPase activity. Moreover, based on an in silico model of the N-terminal region of MtbEccC5DUF, it is hypothesized that MtbEccC5DUF is a degenerated ATPase domain that may have retained the ability to hexamerize. These observations draw attention to DUF domains as structural elements with potential implications in the opening and closure of the membrane pore during the secretion process via their involvement in inter-protomer interactions.

Fragment-based drug design using X-ray crystallography is a powerful technique to enable the development of new lead compounds, or probe molecules, against biological targets. This study addresses the need to determine fragment binding orientations for low-occupancy fragments with incomplete electron density, an essential step before further development of the molecule. Halogen atoms play multiple roles in drug discovery due to their unique combination of electronegativity, steric effects and hydrophobic properties. Fragments incorporating halogen atoms serve as promising starting points in hit-to-lead development as they often establish halogen bonds with target proteins, potentially enhancing binding affinity and selectivity, as well as counteracting drug resistance. Here, the aim was to unambiguously identify the binding orientations of fragment hits for SARS-CoV-2 nonstructural protein 1 (nsp1) which contain a combination of sulfur and/or chlorine, bromine and iodine substituents. The binding orientations of carefully selected nsp1 analogue hits were focused on by employing their anomalous scattering combined with Pan-Dataset Density Analysis (PanDDA). Anomalous difference Fourier maps derived from the diffraction data collected at both standard and long-wavelength X-rays were compared. The discrepancies observed in the maps of iodine-containing fragments collected at different energies were attributed to site-specific radiation-damage stemming from the strong X-ray absorption of I atoms, which is likely to cause cleavage of the C—I bond. A reliable and effective data-collection strategy to unambiguously determine the binding orientations of low-occupancy fragments containing sulfur and/or halogen atoms while mitigating radiation damage is presented.

For cryo-electron tomography (cryo-ET) of beam-sensitive biological specimens, a planar sample geometry is typically used. As the sample is tilted, the effective thickness of the sample along the direction of the electron beam increases and the signal-to-noise ratio concomitantly decreases, limiting the transfer of information at high tilt angles. In addition, the tilt range where data can be collected is limited by a combination of various sample-environment constraints, including the limited space in the objective lens pole piece and the possible use of fixed conductive braids to cool the specimen. Consequently, most tilt series are limited to a maximum of ±70°, leading to the presence of a missing wedge in Fourier space. The acquisition of cryo-ET data without a missing wedge, for example using a cylindrical sample geometry, is hence attractive for volumetric analysis of low-symmetry structures such as organelles or vesicles, lysis events, pore formation or filaments for which the missing information cannot be compensated by averaging techniques. Irrespective of the geometry, electron-beam damage to the specimen is an issue and the first images acquired will transfer more high-resolution information than those acquired last. There is also an inherent trade-off between higher sampling in Fourier space and avoiding beam damage to the sample. Finally, the necessity of using a sufficient electron fluence to align the tilt images means that this fluence needs to be fractionated across a small number of images; therefore, the order of data acquisition is also a factor to consider. Here, an n-helix tilt scheme is described and simulated which uses overlapping and interleaved tilt series to maximize the use of a pillar geometry, allowing the entire pillar volume to be reconstructed as a single unit. Three related tilt schemes are also evaluated that extend the continuous and classic dose-symmetric tilt schemes for cryo-ET to pillar samples to enable the collection of isotropic information across all spatial frequencies. A fourfold dose-symmetric scheme is proposed which provides a practical compromise between uniform information transfer and complexity of data acquisition.

The expansive scientific software ecosystem, characterized by millions of titles across various platforms and formats, poses significant challenges in maintaining reproducibility and provenance in scientific research. The diversity of independently developed applications, evolving versions and heterogeneous components highlights the need for rigorous methodologies to navigate these complexities. In response to these challenges, the SBGrid team builds, installs and configures over 530 specialized software applications for use in the on-premises and cloud-based computing environments of SBGrid Consortium members. To address the intricacies of supporting this diverse application collection, the team has developed the Capsule Software Execution Environment, generally referred to as Capsules. Capsules rely on a collection of programmatically generated bash scripts that work together to isolate the runtime environment of one application from all other applications, thereby providing a transparent cross-platform solution without requiring specialized tools or elevated account privileges for researchers. Capsules facilitate modular, secure software distribution while maintaining a centralized, conflict-free environment. The SBGrid platform, which combines Capsules with the SBGrid collection of structural biology applications, aligns with FAIR goals by enhancing the findability, accessibility, interoperability and reusability of scientific software, ensuring seamless functionality across diverse computing environments. Its adaptability enables application beyond structural biology into other scientific fields.

The determination of the atomic resolution structure of biomacromolecules is essential for understanding details of their function. Traditionally, such a structure determination has been performed with crystallographic or nuclear resonance methods, but during the last decade, cryogenic transmission electron microscopy (cryo-TEM) has become an equally important tool. As the blotting and flash-freezing of the samples can induce conformational changes, external validation tools are required to ensure that the vitrified samples are representative of the solution. Although many validation tools have already been developed, most of them rely on fully resolved atomic models, which prevents early screening of the cryo-TEM maps. Here, a novel and automated method for performing such a validation utilizing small-angle X-ray scattering measurements, publicly available through the new software package AUSAXS, is introduced and implemented. The method has been tested on both simulated and experimental data, where it was shown to work remarkably well as a validation tool. The method provides a dummy atomic model derived from the EM map which best represents the solution structure.

The Azotobacter vinelandii FeSII protein forms an oxygen-resistant complex with the nitrogenase MoFe and Fe proteins. FeSII is an adrenodoxin-type ferredoxin that forms a dimer in solution. Previously, the crystal structure was solved [Schlesier et al. (2016), J. Am. Chem. Soc. 138, 239–247] with five copies in the asymmetric unit. One copy is a normal adrenodoxin domain that forms a dimer with its crystallographic symmetry mate. The other four copies are in an `open' conformation with a loop flipped out exposing the 2Fe–2S cluster. The open and closed conformations were interpreted as oxidized and reduced, respectively, and the large conformational change in the open configuration allowed binding to nitrogenase. Here, the structure of FeSII was independently solved in the same crystal form. The positioning of the atoms in the unit cell is similar to the earlier report. However, the interpretation of the structure is different. The `open' conformation is interpreted as the product of a crystallization-induced domain swap. The 2Fe–2S cluster is not exposed to solvent, but in the crystal its interacting helix is replaced by the same helix residues from a crystal symmetry mate. The domain swap is complicated, as it is unusual in being in the middle of the protein rather than at a terminus, and it creates arrangements of molecules that can be interpreted in multiple ways. It is also cautioned that crystal structures should be interpreted in terms of the contents of the entire crystal rather than of one asymmetric unit.

Histidine can be protonated on either or both of the two N atoms of the imidazole moiety. Each of the three possible forms occurs as a result of the stereochemical environment of the histidine side chain. In an atomic model, comparing the possible protonation states in situ, looking at possible hydrogen bonding and metal coordination, it is possible to predict which is most likely to be correct. A more direct method is described that uses quantum-mechanical methods to calculate, also in situ, the minimum geometry and energy for comparison, and therefore to more accurately identify the most likely proton­ation state.

The Mycobacterium tuberculosis trifunctional enzyme (MtTFE) is an α2β2 tetrameric enzyme in which the α-chain harbors the 2E-enoyl-CoA hydratase (ECH) and 3S-hydroxyacyl-CoA dehydrogenase (HAD) active sites, and the β-chain provides the 3-ketoacyl-CoA thiolase (KAT) active site. Linear, medium-chain and long-chain 2E-enoyl-CoA molecules are the preferred substrates of MtTFE. Previous crystallographic binding and modeling studies identified binding sites for the acyl-CoA substrates at the three active sites, as well as the NAD binding pocket at the HAD active site. These studies also identified three additional CoA binding sites on the surface of MtTFE that are different from the active sites. It has been proposed that one of these additional sites could be of functional relevance for the substrate channeling (by surface crawling) of reaction intermediates between the three active sites. Here, 226 fragments were screened in a crystallographic fragment-binding study of MtTFE crystals, resulting in the structures of 16 MtTFE–fragment complexes. Analysis of the 121 fragment-binding events shows that the ECH active site is the `binding hotspot' for the tested fragments, with 41 binding events. The mode of binding of the fragments bound at the active sites provides additional insight into how the long-chain acyl moiety of the substrates can be accommodated at their proposed binding pockets. In addition, the 20 fragment-binding events between the active sites identify potential transient binding sites of reaction intermediates relevant to the possible channeling of substrates between these active sites. These results provide a basis for further studies to understand the functional relevance of the latter binding sites and to identify substrates for which channeling is crucial.

Several proteins from plant pathogenesis-related family 10 (PR10) are highly abundant in the latex of opium poppy and have recently been shown to play diverse and important roles in the biosynthesis of benzylisoquinoline alkaloids (BIAs). The recent determination of the first crystal structures of PR10-10 showed how large conformational changes in a surface loop and adjacent β-strand are coupled to the binding of BIA compounds to the central hydrophobic binding pocket. A more detailed analysis of these conformational changes is now reported to further clarify how ligand binding is coupled to the formation and cleavage of an intermolecular disulfide bond that is only sterically allowed when the BIA binding pocket is empty. To decouple ligand binding from disulfide-bond formation, each of the two highly conserved cysteine residues (Cys59 and Cys155) in PR10-10 was replaced with serine using site-directed mutagenesis. Crystal structures of the Cys59Ser mutant were determined in the presence of papaverine and in the absence of exogenous BIA compounds. A crystal structure of the Cys155Ser mutant was also determined in the absence of exogenous BIA compounds. All three of these crystal structures reveal conformations similar to that of wild-type PR10-10 with bound BIA compounds. In the absence of exogenous BIA compounds, the Cys59Ser and Cys155Ser mutants appear to bind an unidentified ligand or mixture of ligands that was presumably introduced during expression of the proteins in Escherichia coli. The analysis of conformational changes triggered by the binding of BIA compounds suggests a molecular mechanism coupling ligand binding to the disruption of an intermolecular disulfide bond. This mechanism may be involved in the regulation of biosynthetic reactions in plants and possibly other organisms.

Eukaryotic TIR (Toll/interleukin-1 receptor protein) domains signal via TIR–TIR interactions, either by self-association or by interaction with other TIR domains. In mammals, TIR domains are found in Toll-like receptors (TLRs) and cytoplasmic adaptor proteins involved in pro-inflammatory signaling. Previous work revealed that the MAL TIR domain (MALTIR) nucleates the assembly of MyD88TIR into crystalline arrays in vitro. A microcrystal electron diffraction (MicroED) structure of the MyD88TIR assembly has previously been solved, revealing a two-stranded higher-order assembly of TIR domains. In this work, it is demonstrated that the TIR domain of TLR2, which is reported to signal as a heterodimer with either TLR1 or TLR6, induces the formation of crystalline higher-order assemblies of MyD88TIR in vitro, whereas TLR1TIR and TLR6TIR do not. Using an improved data-collection protocol, the MicroED structure of TLR2TIR-induced MyD88TIR microcrystals was determined at a higher resolution (2.85 Å) and with higher completeness (89%) compared with the previous structure of the MALTIR-induced MyD88TIR assembly. Both assemblies exhibit conformational differences in several areas that are important for signaling (for example the BB loop and CD loop) compared with their monomeric structures. These data suggest that TLR2TIR and MALTIR interact with MyD88 in an analogous manner during signaling, nucleating MyD88TIR assemblies uni­directionally.

Crystal polymorphism serves as a strategy to study the conformational flexibility of proteins. However, the relationship between protein crystal packing and protein conformation often remains elusive. In this study, two distinct crystal forms of a green fluorescent protein variant, NowGFP, are compared: a previously identified monoclinic form (space group C2) and a newly discovered ortho­rhombic form (space group P212121). Comparative analysis reveals that both crystal forms exhibit nearly identical linear assemblies of NowGFP molecules interconnected through similar crystal contacts. However, a notable difference lies in the stacking of these assemblies: parallel in the monoclinic form and perpendicular in the orthorhombic form. This distinct mode of stacking leads to different crystal contacts and induces structural alteration in one of the two molecules within the asymmetric unit of the orthorhombic crystal form. This new conformational state captured by orthorhombic crystal packing exhibits two unique features: a conformational shift of the β-barrel scaffold and a restriction of pH-dependent shifts of the key residue Lys61, which is crucial for the pH-dependent spectral shift of this protein. These findings demonstrate a clear connection between crystal packing and alternative conformational states of proteins, providing insights into how structural variations influence the function of fluorescent proteins.

β-Glucosidase from the thermophilic bacterium Caldicellulosiruptor saccharo­lyticus (Bgl1) has been denoted as having an attractive catalytic profile for various industrial applications. Bgl1 catalyses the final step of in the decomposition of cellulose, an unbranched glucose polymer that has attracted the attention of researchers in recent years as it is the most abundant renewable source of reduced carbon in the biosphere. With the aim of enhancing the thermostability of Bgl1 for a broad spectrum of biotechnological processes, it has been subjected to structural studies. Crystal structures of Bgl1 and its complex with glucose were determined at 1.47 and 1.95 Å resolution, respectively. Bgl1 is a member of glycosyl hydrolase family 1 (GH1 superfamily, EC 3.2.1.21) and the results showed that the 3D structure of Bgl1 follows the overall architecture of the GH1 family, with a classical (β/α)8 TIM-barrel fold. Comparisons of Bgl1 with sequence or structural homologues of β-glucosidase reveal quite similar structures but also unique structural features in Bgl1 with plausible functional roles.

A group of three deep-learning tools, referred to collectively as CHiMP (Crystal Hits in My Plate), were created for analysis of micrographs of protein crystallization experiments at the Diamond Light Source (DLS) synchrotron, UK. The first tool, a classification network, assigns images into categories relating to experimental outcomes. The other two tools are networks that perform both object detection and instance segmentation, resulting in masks of individual crystals in the first case and masks of crystallization droplets in addition to crystals in the second case, allowing the positions and sizes of these entities to be recorded. The creation of these tools used transfer learning, where weights from a pre-trained deep-learning network were used as a starting point and repurposed by further training on a relatively small set of data. Two of the tools are now integrated at the VMXi macromolecular crystallography beamline at DLS, where they have the potential to absolve the need for any user input, both for monitoring crystallization experiments and for triggering in situ data collections. The third is being integrated into the XChem fragment-based drug-discovery screening platform, also at DLS, to allow the automatic targeting of acoustic compound dispensing into crystallization droplets.

The availability of highly accurate protein structure predictions from AlphaFold2 (AF2) and similar tools has hugely expanded the applicability of molecular replacement (MR) for crystal structure solution. Many structures can be solved routinely using raw models, structures processed to remove unreliable parts or models split into distinct structural units. There is therefore an open question around how many and which cases still require experimental phasing methods such as single-wavelength anomalous diffraction (SAD). Here, this question is addressed using a large set of PDB depositions that were solved by SAD. A large majority (87%) could be solved using unedited or minimally edited AF2 predictions. A further 18 (4%) yield straightforwardly to MR after splitting of the AF2 prediction using Slice'N'Dice, although different splitting methods succeeded on slightly different sets of cases. It is also found that further unique targets can be solved by alternative modelling approaches such as ESMFold (four cases), alternative MR approaches such as ARCIMBOLDO and AMPLE (two cases each), and multimeric model building with AlphaFold-Multimer or UniFold (three cases). Ultimately, only 12 cases, or 3% of the SAD-phased set, did not yield to any form of MR tested here, offering valuable hints as to the number and the characteristics of cases where experimental phasing remains essential for macromolecular structure solution.

Pollution from plastics is a global problem that threatens the biosphere for a host of reasons, including the time scale that it takes for most plastics to degrade. Biodegradation is an ideal solution for remediating bioplastic waste as it does not require the high temperatures necessary for thermal degradation and does not introduce additional pollutants into the environment. Numerous organisms can scavenge for bioplastics, such as polylactic acid (PLA) or poly-(R)-hydroxybutyrate (PHB), which they can use as an energy source. Recently, a promiscuous PHBase from the thermophilic soil bacterium Lihuaxuella thermophila (LtPHBase) was identified. LtPHBase can accommodate many substrates, including PHB granules and films and PHB block copolymers, as well as the unrelated polymers polylactic acid (PLA) and polycaprolactone (PCL). LtPHBase uses the expected Ser–His–Asp catalytic triad for hydrolysis at an optimal enzyme activity near 70°C. Here, the 1.75 Å resolution crystal structure of apo LtPHBase is presented and its chemical stability is profiled. Knowledge of its substrate preferences was extended to different-sized PHB granules. It is shown that LtPHBase is highly resistant to unfolding, with barriers typical for thermophilic enzymes, and shows a preference for low-molecular-mass PHB granules. These insights have implications for the long-term potential of LtPHBase as an industrial PHB hydrolase and shed light on the evolutionary role that this enzyme plays in bacterial metabolism.

For protein crystals in which more than two thirds of the volume is occupied by solvent, the featureless nature of the solvent region often generates a constraint that is powerful enough to allow direct phasing of X-ray diffraction data. Practical implementation relies on the use of iterative projection algorithms with good global convergence properties to solve the difficult nonconvex phase-retrieval problem. In this paper, some aspects of phase retrieval using iterative projection algorithms are systematically explored, where the diffraction data and density-value distributions in the protein and solvent regions provide the sole constraints. The analysis is based on the addition of random error to the phases of previously determined protein crystal structures, followed by evaluation of the ability to recover the correct phase set as the distance from the solution increases. The properties of the difference-map (DM), relaxed–reflect–reflect (RRR) and relaxed averaged alternating reflectors (RAAR) algorithms are compared. All of these algorithms prove to be effective for crystallographic phase retrieval, and the useful ranges of the adjustable parameter which controls their behavior are established. When these algorithms converge to the solution, the algorithm trajectory becomes stationary; however, the density function continues to fluctuate significantly around its mean position. It is shown that averaging over the algorithm trajectory in the stationary region, following convergence, improves the density estimate, with this procedure outperforming previous approaches for phase or density refinement.

Fatty acid-derivative prodrugs have been utilized extensively to improve the physicochemical, biopharmaceutical and pharmacokinetic properties of active pharmaceutical ingredients. However, to our knowledge, the crystallization behavior of prodrugs modified with different fatty acids has not been explored. In the present work, a series of paliperidone aliphatic prodrugs with alkyl chain lengths ranging from C4 to C16 was investigated with respect to crystal structure, crystal morphology and crystallization kinetics. The paliperidone derivatives exhibited isostructural crystal packing, despite the different alkyl chain lengths, and crystallized with the dominant (100) face in both melt and solution. The rate of crystallization for paliperidone derivatives in the melt increases with alkyl chain length owing to greater molecular mobility. In contrast, the longer chains prolong the nucleation induction time and reduce the crystal growth kinetics in solution. The results show a correlation between difficulty of nucleation in solution and the interfacial energy. This work provides insight into the crystallization behavior of paliperidone aliphatic prodrugs and reveals that the role of alkyl chain length in the crystallization behavior has a strong dependence on the crystallization method.

Serial electron diffraction (SerialED), which applies a snapshot data acquisition strategy for each crystal, was introduced to tackle the problem of radiation damage in the structure determination of beam-sensitive materials by three-dimensional electron diffraction (3DED). The snapshot data acquisition in SerialED can be realized using both transmission and scanning transmission electron microscopes (TEM/STEM). However, the current SerialED workflow based on STEM setups requires special external devices and software, which limits broader adoption. Here, we present a simplified experimental implementation of STEM-based SerialED on Thermo Fisher Scientific STEMs using common proprietary software interfaced through Python scripts to automate data collection. Specifically, we utilize TEM Imaging and Analysis (TIA) scripting and TEM scripting to access the STEM functionalities of the microscope, and DigitalMicrograph scripting to control the camera for snapshot data acquisition. Data analysis adapts the existing workflow using the software CrystFEL, which was developed for serial X-ray crystallography. Our workflow for STEM SerialED can be used on any Gatan or Thermo Fisher Scientific camera. We apply this workflow to collect high-resolution STEM SerialED data from two aluminosilicate zeolites, zeolite Y and ZSM-25. We demonstrate, for the first time, ab initio structure determination through direct methods using STEM SerialED data. Zeolite Y is relatively stable under the electron beam, and STEM SerialED data extend to 0.60 Å. We show that the structural model obtained using STEM SerialED data merged from 358 crystals is nearly identical to that using continuous rotation electron diffraction data from one crystal. This demonstrates that accurate structures can be obtained from STEM SerialED. Zeolite ZSM-25 is very beam-sensitive and has a complex structure. We show that STEM SerialED greatly improves the data resolution of ZSM-25, compared with serial rotation electron diffraction (SerialRED), from 1.50 to 0.90 Å. This allows, for the first time, the use of standard phasing methods, such as direct methods, for the ab initio structure determination of ZSM-25.

The orientation ordering and assembly behavior of silica–nickel Janus particles in a static external magnetic field were probed by ultra small-angle X-ray scattering (USAXS). Even in a weak applied field, the net magnetic moments of the individual particles aligned in the direction of the field, as indicated by the anisotropy in the recorded USAXS patterns. X-ray photon correlation spectroscopy (XPCS) measurements on these suspensions revealed that the corresponding particle dynamics are primarily Brownian diffusion [Zinn, Sharpnack & Narayanan (2023). Soft Matter, 19, 2311–2318]. At higher fields, the magnetic forces led to chain-like configurations of particles, as indicated by an additional feature in the USAXS pattern. A theoretical framework is provided for the quantitative interpretation of the observed anisotropic scattering diagrams and the corresponding degree of orientation. No anisotropy was detected when the magnetic field was applied along the beam direction, which is also replicated by the model. The method presented here could be useful for the interpretation of oriented scattering patterns from a wide variety of particulate systems. The combination of USAXS and XPCS is a powerful approach for investigating asymmetric colloidal particles in external fields.

The application of grazing-incidence total X-ray scattering (GITXS) for pair distribution function (PDF) analysis using >50 keV X-rays from synchrotron light sources has created new opportunities for structural characterization of supported thin films with high resolution. Compared with grazing-incidence wide-angle X-ray scattering, which is only useful for highly ordered materials, GITXS/PDFs expand such analysis to largely disordered or nanostructured materials by examining the atomic pair correlations dependent on the direction relative to the surface of the supporting substrate. A characterization of nanocrystalline In2O3-derived thin films is presented here with in-plane-isotropic and out-of-plane-anisotropic orientational ordering of the atomic structure, each synthesized using different techniques. The atomic orientations of such films are known to vary based on the synthetic conditions. Here, an azimuthal orientational analysis of these films using GITXS with a single incident angle is shown to resolve the markedly different orientations of the atomic structures with respect to the planar support and the different degrees of long-range order, and hence, the terminal surface chemistries. It is anticipated that orientational analysis of GITXS/PDF data will offer opportunities to extend structural analyses of thin films by providing a means to qualitatively determine the major atomic orientation within nanocrystalline and, eventually, non-crystalline films.

Interoperability of crystallographic data with other disciplines is essential for the smooth and rapid progress of structure-based science in the computer age. Within crystallography and closely related subject areas, there is already a high level of conformance to the generally accepted FAIR principles (that data be findable, accessible, interoperable and reusable) through the adoption of common information exchange protocols by databases, publishers, instrument vendors, experimental facilities and software authors. Driven by the success within these domains, the IUCr has worked closely with CODATA (the Committee on Data of the International Science Council) to help develop the latter's commitment to cross-domain integration of discipline-specific data. The IUCr has, in particular, emphasized the need for standards relating to data quality and completeness as an adjunct to the FAIR data landscape. This can ensure definitive reusable data, which in turn can aid interoperability across domains. A microsymposium at the IUCr 2023 Congress provided an up-to-date survey of data interoperability within and outside of crystallography, expounded using a broad range of examples.

Vancomycin is a glycopeptide antibiotic that for decades has been a mainstay of treatment for persistent bacterial infections. However, the spread of antibiotic resistance threatens its continued utility. In particular, vancomycin-resistant enterococci (VRE) have become a pressing clinical challenge. Vancomycin acts by binding and sequestering the intermediate Lipid II in cell-wall biosynthesis, specifically recognizing a d-alanine-d-alanine dipeptide motif within the Lipid II molecule. VRE achieve resistance by remodeling this motif to either d-alanine-d-lactate or d-alanine-d-serine; the former substitution essentially abolishes recognition by vancomycin of Lipid II, whereas the latter reduces the affinity of the antibiotic by roughly one order of magnitude. The complex of vancomycin bound to d-alanine-d-serine has been crystallized, and its 1.20 Å X-ray crystal structure is presented here. This structure reveals that the d-alanine-d-serine ligand is bound in essentially the same position and same pose as the native d-alanine-d-alanine ligand. The serine-containing ligand appears to be slightly too large to be comfortably accommodated in this way, suggesting one possible contribution to the reduced binding affinity. In addition, two flexible hydroxyl groups – one from the serine side chain of the ligand, and the other from a glucose sugar on the antibiotic – are locked into single conformations in the complex, which is likely to contribute an unfavorable entropic component to the recognition of the serine-containing ligand.

In the search for new active pharmaceutical ingredients, the precise control of the chemistry of cocrystals becomes essential. One crucial step within this chemistry is proton migration between cocrystal coformers to form a salt, usually anticipated by the empirical ΔpKa rule. Due to the effective role it plays in modifying intermolecular distances and interactions, pressure adds a new dimension to the ΔpKa rule. Still, this variable has been scarcely applied to induce proton-transfer reactions within these systems. In our study, high-pressure X-ray diffraction and Raman spectroscopy experiments, supported by DFT calculations, reveal modifications to the protonation states of the 4,4'-bi­pyridine (BIPY) and malonic acid (MA) cocrystal (BIPYMA) that allow the conversion of the cocrystal phase into ionic salt polymorphs. On compression, neutral BIPYMA and monoprotonated (BIPYH+MA−) species coexist up to 3.1 GPa, where a phase transition to a structure of P21/c symmetry occurs, induced by a double proton-transfer reaction forming BIPYH22+MA2−. The low-pressure C2/c phase is recovered at 2.4 GPa on decompression, leading to a 0.7 GPa hysteresis pressure range. This is one of a few studies on proton transfer in multicomponent crystals that shows how susceptible the interconversion between differently charged species is to even slight pressure changes, and how the proton transfer can be a triggering factor leading to changes in the crystal symmetry. These new data, coupled with information from previous reports on proton-transfer reactions between coformers, extend the applicability of the ΔpKa rule incorporating the pressure required to induce salt formation.

This paper uses deep learning to present a proof-of-concept for data-driven chemistry in single-molecule magnets (SMMs). Previous discussions within SMM research have proposed links between molecular structures (crystal structures) and single-molecule magnetic properties; however, these have only interpreted the results. Therefore, this study introduces a data-driven approach to predict the properties of SMM structures using deep learning. The deep-learning model learns the structural features of the SMM molecules by extracting the single-molecule magnetic properties from the 3D coordinates presented in this paper. The model accurately determined whether a molecule was a single-molecule magnet, with an accuracy rate of approximately 70% in predicting the SMM properties. The deep-learning model found SMMs from 20 000 metal complexes extracted from the Cambridge Structural Database. Using deep-learning models for predicting SMM properties and guiding the design of novel molecules is promising.

Evidence that the electronic structure of atoms persists in molecules to a much greater extent than has been usually admitted is presented. This is achieved by resorting to N-electron real-space descriptors instead of one- or at most two-particle projections like the electron or exchange-correlation densities. Here, the 3N-dimensional maxima of the square of the wavefunction, the so-called Born maxima, are used. Since this technique is relatively unknown to the crystallographic community, a case-based approach is taken, revisiting first the Born maxima of atoms in their ground state and then some of their excited states. It is shown how they survive in molecules and that, beyond any doubt, the distribution of electrons around an atom in a molecule can be recognized as that of its isolated, in many cases excited, counterpart, relating this fact with the concept of energetic promotion. Several other cases that exemplify the applicability of the technique to solve chemical bonding conflicts and to introduce predictability in real-space analyses are also examined.