ca

BGS Open Day 2024 - British Geological Survey

BGS Open Day 2024  British Geological Survey











ca

2-Amino-5-oxo-4-(thio­phen-2-yl)-5,6,7,8-tetra­hydro-4H-chromene-3-carbo­nitrile

In the crystal structure of the title compound, inter­molecular N—H⋯N and N—H⋯O hydrogen bonds form a two-dimensional supra­molecular network along the ac plane, contributing to the cohesion of the crystal.




ca

2,4-Di­chloro-6-{N-[2-(tri­fluoro­meth­yl)phen­yl]carboximido­yl}phenol

The title compound was synthesized by the condensation between tri­fluoro­methyl­aniline and di­chloro­salicyl­aldehyde by nucleophilic addition, forming a hemiaminal, followed by a dehydration to generate an imine. The compound crystallizes in an ortho­rhom­bic Pbca (Z = 8) space group with a dihedral angle of 44.70 (5)° between the two aromatic rings. In the crystal, the mol­ecules pack together to form a zigzag pattern along the c axis.




ca

Ethyl (2RS,3SR,4RS)-1-ethyl-2-(furan-2-yl)-4-hy­droxy-5-oxopyrrolidine-3-carboxyl­ate

The crystal structure of a pyrrolidine analogue obtained from the stereoselective reduction of the enolic form of 4-hy­droxy-2-furyl-pyrrole­carboxyl­ate is described.




ca

From `crystallographic accuracy' to `thermodynamic accuracy': a redetermination of the crystal structure of calcium atorvastatin trihydrate (Lipitor®)

The crystal structure of calcium atorvastatin trihydrate was redetermined from previously published synchrotron powder diffraction data to give a much-improved agreement with two independent density-functional theory calculations.




ca

Atypical phase transition, twinning and ferroelastic domain structure in bis(ethylenediammonium) tetrabromozincate(II) bromide, [NH3(CH2)2NH3]2[ZnBr4]Br2

A unique phase transition, twinning and ferroelastic domain structure in [NH3(CH2)2NH3]2[ZnBr4]Br2 is found. The new additional domain structure is observed at the phase transition on heating, which is preserved after cooling to room temperature.




ca

Morphological control for hollow rod crystals of a photochromic di­aryl­ethene on spherulites by surface properties of substrates

1,2-Bis(3,5-di­methyl-2-thienyl)perfluoro­cyclo­pentene formed its own spherulites by sublimation onto the hydro­philic surfaces of the (0001) planes of α-quartz and sapphire substrates. The formation of different morphologies of these spherulites was attributed to the surface properties of each substrate. Depending on the morphology of the spherulites, hollow rod crystals with cross sections of different sizes and shapes and branching structures were generated on the surfaces of the spherulites.




ca

Selective Acceleration and Inhibition of Crystal Growth of Glass Carbamazepine by Low-Concentration Poly(ethylene oxide):Effects of Drug Polymorph

Low-concentration poly(ethylene oxide) exhibit the polymorph-dependent effects on both the surface and bulk crystal growth of carbamazepine polymorphs. These polymorph-dependent effects of PEO were mainly attributed to the polymer enrichment at the interface and different crystal surface-polymer interactions.




ca

Atypical phase transition, twinning and ferroelastic domain structure in bis(ethylenediammonium) tetrabromozincate(II) bromide, [NH3(CH2)2NH3]2[ZnBr4]Br2

Single-crystal growth, differential thermal analysis (DTA), derivative thermogravimetry (DTG), differential scanning calorimetry (DSC), X-ray structural studies and polarized microscopy observations of bis(ethylenediammonium) tetrabromozincate(II) bromide [NH3(CH2)2NH3]2[ZnBr4]Br2 are presented. A reversible phase transition is described. At room temperature, the complex crystallizes in the monoclinic system. In some cases, the single crystals are twinned into two or more large domains of ferroelastic type with domain walls in the (100) crystallographic plane. DTA and DTG measurements show chemical stability of the crystal up to ∼538 K. In the DSC studies, a reversible isostructural phase transition was revealed at ∼526/522 K on heating/cooling run, respectively. Optical observation on the heating run reveals that at the phase transition the plane of twinning (domain wall) does not disappear and additionally the appearance of a new domain structure of ferroelastic type with domain walls in the planes (101), (101), (100) and (001) is observed. The domain structure pattern is preserved after cooling to the room-temperature phase and the symmetry of this phase is unchanged.




ca

Morphological control for hollow rod crystals of a photochromic diarylethene on spherulites by surface properties of substrates

Sublimation methods utilizing the surface properties of substrates can address the challenge of controlling hollow morphologies in rod crystals. Spherulites were formed on the hydrophilic surface of the (0001) planes of α-quartz and sapphire substrates by sublimation of 1,2-bis(3,5-dimethyl-2-thienyl)perfluorocyclopentene (1a). Various types of hollow morphologies, distinguished by the size and shape of their cross sections and by the presence or absence of branching structures, were formed separately on α-quartz and sapphire substrates. Such precise control of the hollow morphologies was attributed to the wettability of each substrate, leading to the formation of spherulites of 1a. In addition, it was indicated that the formation process of the surface morphologies of spherulites was associated with the hollow morphologies of rod crystals of 1a.




ca

From `crystallographic accuracy' to `thermodynamic accuracy': a redetermination of the crystal structure of calcium atorvastatin trihydrate (Lipitor®)

With ever-improving quantum-mechanical computational methods, the accuracy requirements for experimental crystal structures increase. The crystal structure of calcium atorvastatin trihydrate, which has 56 degrees of freedom when determined with a real-space algorithm, was determined from powder diffraction data by Hodge et al. [Powder Diffr. (2020), 35, 136–143]. The crystal structure was a good fit to the experimental data, indicating that the electron density had been captured essentially correctly, but two independent quantum-mechanical calculations disagreed with the experimental structure and with each other. Using the same experimental data, the crystal structure was redetermined from scratch and it was shown that it can be reproduced within a root-mean-square Cartesian displacement of 0.1 Å by two independent quantum-mechanical calculations. The consequences for the calculated energies and solubilities are described.




ca

Polymorph sampling with coupling to extended variables: enhanced sampling of polymorph energy landscapes and free energy perturbation of polymorph ensembles

A novel approach to computationally enhance the sampling of molecular crystal structures is proposed and tested. This method is based on the use of extended variables coupled to a Monte Carlo based crystal polymorph generator. Inspired by the established technique of quasi-random sampling of polymorphs using the rigid molecule constraint, this approach represents molecular clusters as extended variables within a thermal reservoir. Polymorph unit-cell variables are generated using pseudo-random sampling. Within this framework, a harmonic coupling between the extended variables and polymorph configurations is established. The extended variables remain fixed during the inner loop dedicated to polymorph sampling, enforcing a stepwise propagation of the extended variables to maintain system exploration. The final processing step results in a polymorph energy landscape, where the raw structures sampled to create the extended variable trajectory are re-optimized without the thermal coupling term. The foundational principles of this approach are described and its effectiveness using both a Metropolis Monte Carlo type algorithm and modifications that incorporate replica exchange is demonstrated. A comparison is provided with pseudo-random sampling of polymorphs for the molecule coumarin. The choice to test a design of this algorithm as relevant for enhanced sampling of crystal structures was due to the obvious relation between molecular structure variables and corresponding crystal polymorphs as representative of the inherent vapor to crystal transitions that exist in nature. Additionally, it is shown that the trajectories of extended variables can be harnessed to extract fluctuation properties that can lead to valuable insights. A novel thermodynamic variable is introduced: the free energy difference between ensembles of Z' = 1 and Z' = 2 crystal polymorphs.





ca

X-ray crystallographic structure of a novel enantiopure chiral isothiourea with potential applications in enantioselective synthesis

The synthesis of a chiral isothiourea, namely, (4aR,8aR)-3-phenyl-4a,5,6,7,8,8a-hexahydrobenzo[4,5]imidazo[2,1-b]thiazol-9-ium bromide, C15H17N2S+·Br−, with potential organocatalytic and anti-inflammatory activity is reported. The preparation of the heterocycle of interest was carried out in two high-yielding steps. The hydrobromide salt of the isothiourea of interest provided suitable crystals for X-ray diffraction analysis, the results of which are reported. Salient observations from this analysis are the near perpendicular arrangement of the phenyl ring and the mean plane of the heterocycle. This conformational characteristic may be relevant with regard the stereoselectivity induced by the chiral isothiourea in asymmetric reactions. Furthermore, evidence was found for the existence of an S...Br− halogen bond.





ca

Analytical models representing X-ray form factors of ions

Parameters in analytical models for X-ray form factors of ions f0(s), based on the inverse Mott–Bethe formula involving a variable number of Gaussians, are determined for a wide range of published data sets {s, f0(s)}. The models reproduce the calculated form-factor values close to what is expected from a uniform statistical distribution with limits determined by their precision. For different ions associated with the same atom, the number of Gaussians in the models decreases with increasing net positive charge.




ca

Modelling dynamical 3D electron diffraction intensities. I. A scattering cluster algorithm

Three-dimensional electron diffraction (3D-ED) is a powerful technique for crystallographic characterization of nanometre-sized crystals that are too small for X-ray diffraction. For accurate crystal structure refinement, however, it is important that the Bragg diffracted intensities are treated dynamically. Bloch wave simulations are often used in 3D-ED, but can be computationally expensive for large unit cell crystals due to the large number of diffracted beams. Proposed here is an alternative method, the `scattering cluster algorithm' (SCA), that replaces the eigen-decomposition operation in Bloch waves with a simpler matrix multiplication. The underlying principle of SCA is that the intensity of a given Bragg reflection is largely determined by intensity transfer (i.e. `scattering') from a cluster of neighbouring diffracted beams. However, the penalty for using matrix multiplication is that the sample must be divided into a series of thin slices and the diffracted beams calculated iteratively, similar to the multislice approach. Therefore, SCA is more suitable for thin specimens. The accuracy and speed of SCA are demonstrated on tri-iso­propyl silane (TIPS) pentacene and rubrene, two exemplar organic materials with large unit cells.




ca

Modelling dynamical 3D electron diffraction intensities. II. The role of inelastic scattering

The strong interaction of high-energy electrons with a crystal results in both dynamical elastic scattering and inelastic events, particularly phonon and plasmon excitation, which have relatively large cross sections. For accurate crystal structure refinement it is therefore important to uncover the impact of inelastic scattering on the Bragg beam intensities. Here a combined Bloch wave–Monte Carlo method is used to simulate phonon and plasmon scattering in crystals. The simulated thermal and plasmon diffuse scattering are consistent with experimental results. The simulations also confirm the empirical observation of a weaker unscattered beam intensity with increasing energy loss in the low-loss regime, while the Bragg-diffracted beam intensities do not change significantly. The beam intensities include the diffuse scattered background and have been normalized to adjust for the inelastic scattering cross section. It is speculated that the random azimuthal scattering angle during inelastic events transfers part of the unscattered beam intensity to the inner Bragg reflections. Inelastic scattering should not significantly influence crystal structure refinement, provided there are no artefacts from any background subtraction, since the relative intensity of the diffracted beams (which includes the diffuse scattering) remains approximately constant in the low energy loss regime.




ca

Parameterized absorptive electron scattering factors

In electron diffraction, thermal atomic motion produces incoherent scattering over a relatively wide angular range, which appears as a diffuse background that is usually subtracted from measurements of Bragg spot intensities in structure solution methods. The transfer of electron flux from Bragg spots to diffuse scatter is modelled using complex scattering factors f + if' in the Bloch wave methodology. In a two-beam Einstein model the imaginary `absorptive' scattering factor f' can be obtained by the evaluation of an integral containing f over all possible scattering angles. While more sophisticated models of diffuse scatter are widely used in the electron microscopy community, it is argued in this paper that this simple model is appropriate for current structure solution and refinement methods. The two-beam model is a straightforward numerical calculation, but even this simplistic approach can become time consuming for simulations of materials with large numbers of atoms in the unit cell and/or many incident beam orientations. Here, a parameterized form of f' is provided for 103 elements as neutral, spherical atoms that reduces calculation time considerably.




ca

Automated selection of nanoparticle models for small-angle X-ray scattering data analysis using machine learning

Small-angle X-ray scattering (SAXS) is widely used to analyze the shape and size of nanoparticles in solution. A multitude of models, describing the SAXS intensity resulting from nanoparticles of various shapes, have been developed by the scientific community and are used for data analysis. Choosing the optimal model is a crucial step in data analysis, which can be difficult and time-consuming, especially for non-expert users. An algorithm is proposed, based on machine learning, representation learning and SAXS-specific preprocessing methods, which instantly selects the nanoparticle model best suited to describe SAXS data. The different algorithms compared are trained and evaluated on a simulated database. This database includes 75 000 scattering spectra from nine nanoparticle models, and realistically simulates two distinct device configurations. It will be made freely available to serve as a basis of comparison for future work. Deploying a universal solution for automatic nanoparticle model selection is a challenge made more difficult by the diversity of SAXS instruments and their flexible settings. The poor transferability of classification rules learned on one device configuration to another is highlighted. It is shown that training on several device configurations enables the algorithm to be generalized, without degrading performance compared with configuration-specific training. Finally, the classification algorithm is evaluated on a real data set obtained by performing SAXS experiments on nanoparticles for each of the instrumental configurations, which have been characterized by transmission electron microscopy. This data set, although very limited, allows estimation of the transferability of the classification rules learned on simulated data to real data.




ca

Permissible domain walls in monoclinic ferroelectrics. Part II. The case of MC phases

Monoclinic ferroelectric phases are prevalent in various functional materials, most notably mixed-ion perovskite oxides. These phases can manifest as regularly ordered long-range crystallographic structures or as macroscopic averages of the self-assembled tetragonal/rhombohedral nanodomains. The structural and physical properties of monoclinic ferroelectric phases play a pivotal role when exploring the interplay between ferroelectricity, ferroelasticity, giant piezoelectricity and multiferroicity in crystals, ceramics and epitaxial thin films. However, the complex nature of this subject presents challenges, particularly in deciphering the microstructures of monoclinic domains. In Paper I [Biran & Gorfman (2024). Acta Cryst. A80, 112–128] the geometrical principles governing the connection of domain microstructures formed by pairing MAB type monoclinic domains were elucidated. Specifically, a catalog was established of `permissible domain walls', where `permissible', as originally introduced by Fousek & Janovec [J. Appl. Phys. (1969), 40, 135–142], denotes a mismatch-free connection between two monoclinic domains along the corresponding domain wall. The present article continues the prior work by elaborating on the formalisms of permissible domain walls to describe domain microstructures formed by pairing the MC type monoclinic domains. Similarly to Paper I, 84 permissible domain walls are presented for MC type domains. Each permissible domain wall is characterized by Miller indices, the transformation matrix between the crystallographic basis vectors of the domains and, crucially, the expected separation of Bragg peaks diffracted from the matched pair of domains. All these parameters are provided in an analytical form for easy and intuitive interpretation of the results. Additionally, 2D illustrations are provided for selected instances of permissible domain walls. The findings can prove valuable for various domain-related calculations, investigations involving X-ray diffraction for domain analysis and the description of domain-related physical properties.




ca

Bond topology of chain, ribbon and tube silicates. Part II. Geometrical analysis of infinite 1D arrangements of (TO4)n− tetrahedra

In Part I of this series, all topologically possible 1-periodic infinite graphs (chain graphs) representing chains of tetrahedra with up to 6–8 vertices (tetrahedra) per repeat unit were generated. This paper examines possible restraints on embedding these chain graphs into Euclidean space such that they are compatible with the metrics of chains of tetrahedra in observed crystal structures. Chain-silicate minerals with T = Si4+ (plus P5+, V5+, As5+, Al3+, Fe3+, B3+, Be2+, Zn2+ and Mg2+) have a grand nearest-neighbour 〈T–T〉 distance of 3.06±0.15 Å and a minimum T⋯T separation of 3.71 Å between non-nearest-neighbour tetrahedra, and in order for embedded chain graphs (called unit-distance graphs) to be possible atomic arrangements in crystals, they must conform to these metrics, a process termed equalization. It is shown that equalization of all acyclic chain graphs is possible in 2D and 3D, and that equalization of most cyclic chain graphs is possible in 3D but not necessarily in 2D. All unique ways in which non-isomorphic vertices may be moved are designated modes of geometric modification. If a mode (m) is applied to an equalized unit-distance graph such that a new geometrically distinct unit-distance graph is produced without changing the lengths of any edges, the mode is designated as valid (mv); if a new geometrically distinct unit-distance graph cannot be produced, the mode is invalid (mi). The parameters mv and mi are used to define ranges of rigidity of the unit-distance graphs, and are related to the edge-to-vertex ratio, e/n, of the parent chain graph. The program GraphT–T was developed to embed any chain graph into Euclidean space subject to the metric restraints on T–T and T⋯T. Embedding a selection of chain graphs with differing e/n ratios shows that the principal reason why many topologically possible chains cannot occur in crystal structures is due to violation of the requirement that T⋯T > 3.71 Å. Such a restraint becomes increasingly restrictive as e/n increases and indicates why chains with stoichiometry TO<2.5 do not occur in crystal structures.




ca

Online carbohydrate 3D structure validation with the Privateer web app

Owing to the difficulties associated with working with carbohydrates, validating glycan 3D structures prior to deposition into the Protein Data Bank has become a staple of the structure-solution pipeline. The Privateer software provides integrative methods for the validation, analysis, refinement and graphical representation of 3D atomic structures of glycans, both as ligands and as protein modifiers. While Privateer is free software, it requires users to install any of the structural biology software suites that support it or to build it from source code. Here, the Privateer web app is presented, which is always up to date and available to be used online (https://privateer.york.ac.uk) without installation. This self-updating tool, which runs locally on the user's machine, will allow structural biologists to simply and quickly analyse carbohydrate ligands and protein glycosylation from a web browser whilst retaining all confidential information on their devices.




ca

Expression, purification and crystallization of the photosensory module of phytochrome B (phyB) from Sorghum bicolor

Sorghum, a short-day tropical plant, has been adapted for temperate grain production, in particular through the selection of variants at the MATURITY loci (Ma1–Ma6) that reduce photoperiod sensitivity. Ma3 encodes phytochrome B (phyB), a red/far-red photochromic biliprotein photoreceptor. The multi-domain gene product, comprising 1178 amino acids, autocatalytically binds the phytochromobilin chromophore to form the photoactive holophytochrome (Sb.phyB). This study describes the development of an efficient heterologous overproduction system which allows the production of large quantities of various holoprotein constructs, along with purification and crystallization procedures. Crystals of the Pr (red-light-absorbing) forms of NPGP, PGP and PG (residues 1–655, 114–655 and 114–458, respectively), each C-terminally tagged with His6, were successfully produced. While NPGP crystals did not diffract, those of PGP and PG diffracted to 6 and 2.1 Å resolution, respectively. Moving the tag to the N-terminus and replacing phytochromobilin with phycocyanobilin as the ligand produced PG crystals that diffracted to 1.8 Å resolution. These results demonstrate that the diffraction quality of challenging protein crystals can be improved by removing flexible regions, shifting fusion tags and altering small-molecule ligands.




ca

Structural and biochemical characterization of the M405S variant of Desulfovibrio vulgaris formate dehydrogenase

Molybdenum- or tungsten-dependent formate dehydrogenases have emerged as significant catalysts for the chemical reduction of CO2 to formate, with biotechnological applications envisaged in climate-change mitigation. The role of Met405 in the active site of Desulfovibrio vulgaris formate dehydrogenase AB (DvFdhAB) has remained elusive. However, its proximity to the metal site and the conformational change that it undergoes between the resting and active forms suggests a functional role. In this work, the M405S variant was engineered, which allowed the active-site geometry in the absence of methionine Sδ interactions with the metal site to be revealed and the role of Met405 in catalysis to be probed. This variant displayed reduced activity in both formate oxidation and CO2 reduction, together with an increased sensitivity to oxygen inactivation.




ca

X-ray crystal structure of a designed rigidified imaging scaffold in the ligand-free conformation

Imaging scaffolds composed of designed protein cages fused to designed ankyrin repeat proteins (DARPins) have enabled the structure determination of small proteins by cryogenic electron microscopy (cryo-EM). One particularly well characterized scaffold type is a symmetric tetrahedral assembly composed of 24 subunits, 12 A and 12 B, which has three cargo-binding DARPins positioned on each vertex. Here, the X-ray crystal structure of a representative tetrahedral scaffold in the apo state is reported at 3.8 Å resolution. The X-ray crystal structure complements recent cryo-EM findings on a closely related scaffold, while also suggesting potential utility for crystallographic investigations. As observed in this crystal structure, one of the three DARPins, which serve as modular adaptors for binding diverse `cargo' proteins, present on each of the vertices is oriented towards a large solvent channel. The crystal lattice is unusually porous, suggesting that it may be possible to soak crystals of the scaffold with small (≤30 kDa) protein cargo ligands and subsequently determine cage–cargo structures via X-ray crystallography. The results suggest the possibility that cryo-EM scaffolds may be repurposed for structure determination by X-ray crystallography, thus extending the utility of electron-microscopy scaffold designs for alternative structural biology applications.




ca

Preliminary X-ray diffraction and ligand-binding analyses of the N-terminal domain of hypothetical protein Rv1421 from Mycobacterium tuberculosis H37Rv

Mycobacterium tuberculosis can reside and persist in deep tissues; latent tuberculosis can evade immune detection and has a unique mechanism to convert it into active disease through reactivation. M. tuberculosis Rv1421 (MtRv1421) is a hypothetical protein that has been proposed to be involved in nucleotide binding-related metabolism in cell-growth and cell-division processes. However, due to a lack of structural information, the detailed function of MtRv1421 remains unclear. In this study, a truncated N-terminal domain (NTD) of MtRv1421, which contains a Walker A/B-like motif, was purified and crystallized using PEG 400 as a precipitant. The crystal of MtRv1421-NTD diffracted to a resolution of 1.7 Å and was considered to belong to either the C-centered monoclinic space group C2 or the I-centered orthorhombic space group I222, with unit-cell parameters a = 124.01, b = 58.55, c = 84.87 Å, β = 133.12° or a = 58.53, b = 84.86, c = 90.52 Å, respectively. The asymmetric units of the C2 or I222 crystals contained two or one monomers, respectively. In terms of the binding ability of MtRv1421-NTD to various ligands, uridine diphosphate (UDP) and UDP-N-acetylglucosamine significantly increased the melting temperature of MtRv1421-NTD, which indicates structural stabilization through the binding of these ligands. Altogether, the results reveal that a UDP moiety may be required for the interaction of MtRv1421-NTD as a nucleotide-binding protein with its ligand.