Exploring serial crystallography for drug discovery

In this work, serial crystallography is applied to a drug discovery target and the room-temperature ligand-bound complexes are used to explore temperature-dependent structural differences.


Introduction
Structure-based drug design is a very well established concept and is integral to drug discovery since it can improve the speed and quality of pharmaceutical candidate compound discovery.Structural methods are commonly used throughout a drugdesign project from early hit-finding to optimization of potency and properties for in vivo delivery (Maveyraud & Mourey, 2020).Thus, all major pharmaceutical companies and many biotech companies have access to structural biology capabilities.Successful structure-based drug design requires structure determination to high quality on a timescale matching that of the design cycle.Conventional macromolecular X-ray crystallography, performed by collecting X-ray diffraction data while continuously rotating a single protein crystal kept at cryogenic temperature, is highly streamlined and the work horse of structure-based drug design.It does, however, require a very robust crystallization protocol and includes the labour-intense step of harvesting crystals.Single-particle cryo-electron microscopy (cryo-EM) has gained much attention during the recent years as an attractive alternative to X-ray crystallography as it does not require that the protein of interest is prone to form well ordered crystals (Guaita et al., 2022).However, both traditional cryo-crystallography and single-particle cryo-EM are limited in that they only provide structural information about the target protein in the frozen state.
Serial crystallography (SX) is a relatively new concept where X-ray diffraction data are collected on micrometresized crystals at room temperature (RT).The method was first developed for use at X-ray free-electron laser facilities, where the extreme intensity of the X-ray pulses required a new way of collecting diffraction data.To circumvent this problem, injection devices were developed that supply a stream of very small crystals across the X-ray beam such that each crystal that intercepts the X-rays gives rise to one diffraction image (Chapman et al., 2011).Pioneering experiments established that, although the data recorded is built up of partial reflections from randomly oriented crystals, it is possible to obtain high-quality structural information using SX (Boutet et al., 2012).Since then, the concept has been further developed at synchrotron facilities (Gati et al., 2014;Martin-Garcia, 2021;Caramello & Royant, 2024), giving rise to the term serial synchrotron crystallography (SSX).There are now a number of purpose-built synchrotron beamlines that allow many more users to explore the method than what XFEL facilities can cater for, including T-REXX at PETRA III, Germany; ID29 at the ESRF, France; and MicroMAX at MAX IV, Sweden.The growth of the SX method has also led to a wide variety of sample delivery methods being developed.This includes different variants of jets, including the high-viscosity extruder injector, suitable for viscous samples such as membrane proteins crystallized in lipidic cubic phase (Weierstall et al., 2014;Vakili et al., 2023;Ghosh et al., 2023).There is also a multitude of different fixed-target devices available where the crystals are dispensed onto a surface that is translated across the X-ray beam in a grid-like fashion.This includes highly specialized chips where the surrounding atmosphere of the crystal can be controlled (Roedig et al., 2017), as well as simple X-ray-transparent membranes that are sandwiched around the crystals and mounted onto a standard crystallography pin (Roedig et al., 2017(Roedig et al., , 2016;;Owen et al., 2017;Carrillo et al., 2023;Mehrabi et al., 2020), enabling data collection at any high-focus synchrotron beamline equipped with a fast translating goniometer.A large advantage of fixed-target devices is that sample consumption can typically be minimized compared with jets, a prerequisite for SX to be a realistic alternative for the majority of proteins.
A scientific area where SX has had great impact is within time-resolved structural studies, where the method has opened up unique possibilities previously out of reach for traditional cryo-crystallography (Bra ¨nde ´n & Neutze, 2021).Of more general interest, it has been suggested that it may be important to study protein structures at RT as opposed to in their frozen state to reduce the risk of temperature artefacts.Several studies have made interesting observations of structural differences related to temperature (Fenwick et al., 2014;Ebrahim et al., 2022;Keedy et al., 2018), including the loss of information about transient binding sites when working at cryogenic temperature (Fischer et al., 2015).Moreover, RT structural information provided an increased understanding of the link between inhibitor potency and binding conformation in the study of a cancer target (Milano et al., 2021).Recently, a crystallographic screen of 143 small-compound fragment binders targeting protein tyrosine phosphatase PTP1B was performed using rotational crystallography at RT (Skaist Mehlman et al., 2023).By comparison with earlier cryocrystallography data, it was shown that fewer ligands bound at RT, but also that unique binding poses were identified in the RT experiment, including interactions at new binding sites that had previously not been observed.Finally, SX has a potential advantage in that it may be possible to avoid the manual handling of crystals required in traditional cryo-crystallography, where the protein crystals are typically harvested and manipulated by hand.Future automation of crystal handling would be highly interesting, not the least for use in high-throughput drug discovery campaigns.
To explore the utility of SX in a drug discovery setting, we selected a pharmaceutically relevant target -soluble epoxide hydrolase (sEH) -as a model system.sEH is an enzyme that binds epoxides and converts them to their corresponding diols (Newman et al., 2003).The human protein is 62.5 kDa in size and consists of a C-terminal and an N-terminal domain connected by a short proline-rich linker (Fig. 1).The Cterminal domain exhibits hydrolase activity and the N-terminal domain is suggested to have phosphatase activity.Three residues within the C-terminal domain active site, ASP335, TYR383 and TYR466, make up the catalytic triad (Fig. 1).The inhibition of sEH can effectively maintain endogenous epoxyeicosatrienoic acid (EET) levels and reduce dihydroxyeicosatrienoic acid (DHET) levels, thus providing therapeutic potential for cardiovascular, central nervous system and metabolic diseases (Imig & Hammock, 2009;Chiamvimonvat et al., 2007;Sun et al., 2022;Shen & Hammock, 2012;Hashimoto, 2019;Va ´zquez et al., 2023).Structures of sEH in complex with a large number of inhibitors have been solved previously using conventional cryocrystallography (O ¨ster et al., 2015).Thus, we set out to generate a set of RT protein-inhibitor complexes to compare with those previously solved and explore whether the use of SX would provide additional structural information of interest in a drug discovery setting.

Protein production and purification
As full length sEH has a disordered C-terminal, a C-terminally truncated construct of sEH(1-548) was expressed in Spodoptera frugiperda (SF9) insect cells as previously described (O ¨ster et al., 2015).100 ml of cells were used for each purification following the published protocol with the exception that two interconnected 5 ml HisTrap FF Crude (Cytivia) columns were used instead of one to avoid saturation of the column.The purified protein was concentrated to 15-20 mg ml À 1 and stored in a buffer containing 20 mM Tris-HCl pH 8, 150 mM NaCl, 1 mM TCEP and 10% glycerol.

Crystallization
A concentrated seed solution was generated by growing larger crystals of sEH, 100-1000 mm in size, following a previously established protocol (O ¨ster et al., 2015).In brief, vapour diffusion crystallizations were set up using a protein concentration of 15-20 mg ml À 1 and a 1:1 ratio of protein: precipitant solution in 24-well sitting drop plates (Cryschem M Plate, Hampton research) with a well solution containing 32-38% PEG3350, 0.1 M LiSO 4 , 0.1 M Tris pH 8.5 and a drop volume of 5-10 ml.The drops were streak-seeded using crystal seeds previously obtained using this protocol.Drops containing well formed crystals were transferred into a 2 ml Eppendorf tube, typically resulting in a volume of �100 ml.Two microseed beads (Molecular Dimensions) were added to the tube and the crystals were crushed by vortexing the sample for 3 � 5 min, or until no large fragments of crystals were visible in an optical microscope.This seed production protocol was inspired by a previously published method to produce seeds for microcrystallization (Dods et al., 2017).
Optimization of batch crystallization conditions was carried out in a micro-batch setting using sitting-drop plates, referred to as hybrid crystallization.This was achieved by composing the well solution such that the concentration of the individual components corresponds to the concentration in the crystallization drop after the addition of protein.A crystallization screen was carried out in this fashion by varying the seed concentration and seed volume.In addition, the concentrations of protein, PEG3350, Li 2 SO 4 and Tris-HCl (pH 8.5) were optimized, as well as the ratio of protein to precipitant solution (see Table S1 of the supporting information for details).
Optimized microcrystals for SX data collection were produced as follows (Fig. 2).A crystallization solution was prepared by mixing 10%(v/v) seed solution into a precipitant solution containing 34% PEG3350, 0.1 M Tris-HCl (pH 8.5) and 0.1 M Li 2 SO 4 .Protein at 14 mg ml À 1 was then mixed with the crystallization solution at a 1:4 protein:precipitant/seed solution ratio in a 0.5 ml Eppendorf tube and the solution was vortexed.The total microcrystallization volume was typically 50-100 ml.The sample was kept upright to allow the crystals to sediment to the bottom of the tube and was then stored at 20 � C.
Each compound was dissolved in 99.5% DMSO to make up a stock solution with a concentration of 500 mM used to prepare the soaking solution.The soaking solution was prepared by diluting the dissolved compound in a buffer [43.5%(w/v)PEG3350, 0.1 M Tris-HCl (pH 8.5) and 0.1 M Li 2 SO 4 ], spun down to remove any precipitated ligand and mixed into the crystal slurry at a final concentration of 10-50 mM.Microcrystals were incubated with the compound for 2-72 h before data collection.

Data collection and processing
All data collection was performed at the MX beamline BioMAX of MAX IV Laboratory using a fixed-target SX setup (Shilova et al., 2020).To increase the crystal hit-rate, the crystals were concentrated by centrifugation for 5-30 s at 1500g.Next, 1.8-2 ml of the crystal solution was sandwiched between two Silson membranes (SiRN-5.0-200-2.5-1000)and manually mounted onto the goniometer using a reusable pin.Room-temperature diffraction data were collected on a Dectris Eiger 16M hybrid pixel detector by raster grid scanning using a 20 � 20 mm X-ray beam with 100% transmission and an exposure time of 11 ms.
The datasets varied in size between 50 000 and 100 000 images.Hit finding, indexing and integration were achieved using indexamajig in the CrystFEL suite (White et al., 2012;White, 2019).Partialator was used for merging and scaling of the data.The high-resolution limit of each dataset was determined based on a CC 1/2 value of at least 30 and an I/�(I) value of at least 0.7, unless the CC 1/2 parameter indicated that the data in a lower-resolution shell were problematic, in which case a lower resolution cut-off was selected.The structures were solved by molecular replacement in the ccp4i suite (Agirre et al., 2023) using a structure of sEH solved by conventional crystallography at cryo-temperature (cryo-T) as a search model (PDB entry 9exm).The model was subsequently refined using BUSTER (Bricogne et al., 2017), while Coot (Emsley & Cowtan, 2004) was used for real-space rebuilding and refinement.Ligand restraints were generated with either writedict or Grade (Smart et al., 2011), and AFITT (Open Eye Scientific Software; Wlodek et al., 2006).The software PyMOL (Schrodinger, 2015) was used to produce the figures.For data processing and refinement statistics, see Table 1.

Development of a new crystallization work-flow for SX.
A key factor for a successful SX study is the ability to generate a sample of microcrystals homogeneous in size, of high quality and in sufficient quantity (Dods et al., 2017).Thus, a robust crystallization protocol is required.For the purpose of generating the quantities needed, batch crystallization has typically been used for growing microcrystals (Beale et al., 2019).A severe limitation when using batch crystallization, however, is that it typically requires very large volumes of protein already at the screening stage.Another disadvantage is that visual inspection of crystallization progress is difficult.Our so-called hybrid crystallization approach, described below, presents a solution to both these problems.
Initial microcrystallization trials of sEH using vapour diffusion in sitting-drop plates with the addition of seeds showed that reducing the protein concentration or increasing the concentrations of Li 2 SO 4 and PEG3350 decreased the size of the crystals.In this study, we aimed for a crystal size of 20-40 mm in length as this was expected to give a good balance between high diffraction and a low number of overlapping crystals on the grids.
The main issue with the initial sitting-drop experiments was the low number of crystals generated in each drop.This is problematic as a low crystal density gives a low hit-rate (i.e. a low number of recorded diffraction images).It was possible to increase the number of crystals by several orders of magnitude by adding seed solution to the drop; however, this also induced the formation of precipitate.The problem was in part overcome by screening each new seed stock and testing dilutions between 1:16 and 1:256 of the seed stock in 43%(w/v) PEG3350 0.1 M Tris-HCl (pH 8.5) and Li 2 SO 4 , where 10%(v/v) seed stock was added to the precipitant solution in each case, to ensure an optimal number of crystallization nuclei.If available, using a device such as a hemocytometer for crystal counting (Boudes et al., 2017) would likely have improved the reproducibility of this step, as has been shown by others (Shoeman et al., 2023).Crystals could be detected after 1 h at 20 � C. To reduce the speed of crystal formation and thereby potentially improve the crystal packing, crystallization experiments were also performed at 12 and 4 � C. The lower temperatures reduced the speed of crystal growth but also induced the formation of aggregates, thus all further setups were carried out at 20 � C.
Batch crystallization is typically the preferred method to produce large amounts of microcrystals.However, in the case of sEH, optimal crystallization conditions for vapour diffusion experiments with drop sizes of 1-3 ml did not translate well into batch conditions with final volumes of 50-100 ml.Since it is not feasible to screen a large variety of conditions in batch due to the high protein consumption, we developed an alternative method that we refer to as hybrid crystallization.
Hybrid crystallization is done in regular sitting-drop plates where the solution in the reservoir well is a mixture of the precipitant solution and the protein buffer in the same ratio as that used between the precipitant solution and protein in the drop.Therefore, the drop size will not change over the course of the experiment as in a regular vapour-diffusion setup.Instead, the experiment is very similar to a batch setup, but without risking that the small-volume drop dries out.The main advantage of this micro-batch setup is that sample consumption can be kept very low while screening for optimal crystallization conditions.Moreover, the hybrid crystallization approach allows the screening to be performed in plates which facilitates visual inspection during the course of the experiment.Finally, the hybrid crystallization method can take advantage of automated dispensing to make it less labour intensive.Using a drop volume of 4-10 ml and a reservoir volume of 500 ml, a large number of conditions were tested.It was found that the protein to precipitant/seed solution ratio influenced the thickness of the crystals (Fig. 3) and thereby the resolution obtained.If the width of a single crystal was below 10 mm then the resolution in most cases would be 3 A ˚or worse.The optimal conditions identified were a precipitant solution of 34% PEG3350, 0.1 M Tris-HCl (pH 8.5) and 0.1 M Li 2 SO 4 with the addition of 10%(v/v) seeds in combination with a protein:precipitant/seed solution ratio of 1:4.This condition could be successfully translated to larger-scale batch setups of 50-100 ml in tubes for microcrystallization of sEH.Microcrystals of sEH were sensitive to transportation, which resulted in crystals either melting or aggregating in transit.This was possibly due to temperature variations during shipping, and the best diffraction was obtained from crystals produced on-site at the synchrotron facility.However, to make SX viable for drug discovery applications, where the optimal working model relies on crystals being produced in-house and then sent to a synchrotron for data collection, it is crucial to develop a logistic chain that can maintain an environment where crystals can be shipped in solution or dispensed onto a chip.

Ligand introduction.
For this study, we selected a subset of ligands that had previously been structurally characterized in complex with sEH using conventional cryo-crystallography (O ¨ster et al.) in order to be able to make a direct comparison of the resulting structures.The selected ligands range from small weakly binding compounds (i.e.fragments) to larger potent inhibitors of sEH (Table S2 and Fig. S1 of the supporting information).Crystals of ligand-bound protein were prepared by soaking the crystals in a solution of highconcentration compound (10-50 mM final concentration) containing 43.5%(w/v) PEG3350, 0.1 M Tris-HCl (pH 8.5) and 0.1 M Li 2 SO 4 for a duration of 2 h or longer.We noted that several of the ligands were not fully solubilized under these conditions.Although it is quite common in soaking experiments that the ligand precipitates, it is seldom detrimental in conventional cryo-crystallography where the crystal is fished out of the drop prior to data collection.However, in fixedtarget SX experiments, the crystal solution is added to a chip and thus ligand precipitate may result in high background and powder diffraction which negatively affect the data quality.Thus, more care must be taken when working with poorly soluble ligands in SX experiments.To reduce the problem caused by ligand precipitation, the ligand solution was centrifuged before addition to the crystal slurry.
3.1.3.Data collection and processing.RT diffraction data were collected at the BioMAX beamline of MAX IV from crystals soaked with seven different compounds.A number of different fixed-target chips were explored for the SX data collection where the best results were achieved using the Silson membranes.We found that an X-ray flux of at least 1 � 10 12 photons s À 1 with 100% transmission were required for good diffraction.Each chip gave rise to 15 000-20 000 images on average and required 15-20 min to assemble and collect data from.Depending on the hit-rate and the number of diffraction images aimed at, a complete dataset was compiled out of data collected from 4-6 chips.Thus, each dataset was acquired in about 1-2 h.Interestingly, the fraction of indexed images varies between 6 and 69%, with compound 2 being the outlier on the low side.The main reason for this large variation lies in how the data collection grids are drawn on the chips (i.e.how carefully areas without crystals or of low crystal density were avoided).Other reasons may be variations in crystal density of the crystal batch, and that some compounds cause precipitate to form in the crystalline sample, thus lowering the indexing rate.As further discussed in Section 3.3, the high redundancy of the resulting data (up to 1700, see Table 1) suggests that the number of images could be drastically reduced.Combined optimization of the crystal density and data collection strategy could thus significantly reduce the data collection time and thereby make SX a more feasible option from a synchrotron usage perspective.Data collection, processing and refinement statistics are presented in Table 1.

Room-temperature structures of sEH and comparison with cryo-temperature structures
In this study, we present seven structures solved by SX at resolutions ranging from 2.1 to 2.5 A ˚resulting from sEH microcrystals soaked with different ligands.In addition, we solved an RT structure in the absence of a ligand (the socalled apo structure) to be used for comparison with the ligated complexes.Structures in complex with the same ligands previously solved at cryo-T were reported to be of similar resolutions [1.9-2.5A ˚ (O ¨ster et al., 2015)].The microcrystals belonged to the space group P6 5 22 which is the same as previously published sEH data.
Overall, the structures solved at RT agree well with the previously reported cryo-crystallography structures with a root-mean-square deviation of the C� positions of 0.3 A between the apo structures solved at cryo-T and RT [PDB Microcrystals of sEH.The effect of using different protein:precipitant/seed solution ratios during screening for optimal conditions using the hybrid crystallization method is shown.(a) A protein:precipitant/seed solution ratio of 1 is used.(b) A protein:precipitant/seed solution ratio of 2 is used.A ratio of 2 gives a lower number of crystals but a more homogeneous sample with fewer overlapping crystals.
5ahx (O ¨ster et al., 2015) and 8qvm, respectively].Interestingly, the apo structure determined by SX at RT harbours an extended-differences density in the active site, which is not present in the apo structure determined by cryocrystallography.The shape of the density agrees well with a bound PEG fragment, and was modelled as such (Fig. 5).

Ligand binding
All the ligands with the exception of compound 7 could be identified in the active site binding pocket (Fig. 4).The large, potent ligands (compounds 1-4) have distinct chemical features, such as the adamatyl moiety of compounds 1 and 2, and the tetrahedral sulfur atoms of compounds 2 and 3 (Fig. S1) which are readily recognizable and could be unambiguously modelled in the electron-density maps.
Among the ligands, three compounds have a molecular weight below 300 Da and can be classified as fragments.Typically, fragments are weak binders and are therefore often associated with ambiguous electron density features due to partial occupancy (Schiebel et al., 2016).In addition, it is not uncommon to find several copies of the same fragment bound in a structure when performing fragment screening due to the high compound concentrations used.Thus, maps from fragment-soaking experiments can be more difficult to interpret than those generated with more potent compounds.
The electron density corresponding to compound 5 has a shape which clearly fits the bi-aryl ring with the adjacent sulfur atom, displaying a stronger peak than the surrounding atoms (Fig. S2).An alternative effort to explain the density as a bound PEG fragment, based on the apo structure, gave an unconvincing result.There is a diffuse electron density feature extending from the hydroxyl group of compound 5 which is modelled as a water molecule.
The electron density map calculated from data collected on crystals soaked with compound 6 is more complicated to interpret and different options were explored (Fig. S3).The omit-map electron density extends over a large part of the binding pocket and cannot be accounted for by a single PEG fragment bound in the central position as in the apo structure.Also, attempts to fit the density with two PEG fragments do not produce a convincing result.Instead, as compound 6 contains a bromide atom, the presence of a strong peak in the omit-map electron density ('left-hand side' of the pocket) clearly identifies the position of one compound 6 molecule bound (although not at full occupancy).This binding pose is identical to that found in the cryo-T structure.Moreover, the omit-map density also suggests the binding of a second copy of compound 6 in the active site ('right-hand side' of the pocket), again similar to what is found in the previously solved cryo-T structure.The density in this part of the pocket is of lower quality and it is more difficult to obtain a completely satisfactory result.Overall, the best fit of the electron density is Figure 4 RT structures of sEH in complex with compounds.The outer images show the different compounds that form complexes numbered according to Table S2.The F o F c omit difference electron density map (green) is contoured at +3.0� in each case.The central panel displays an overlay of the seven RT sEH complex structures zoomed-in on the active site.
obtained with a model that includes 50% partially occupied PEG fragment in the same position as in the apo structure, in combination with two copies of compound 6 at 50% occupancy (Fig. S3).
Finally, the structure generated from microcrystals soaked with compound 7 shows a large positive difference electrondensity feature in the active site.However, the compound cannot be reasonably modelled and the density is instead well fitted with a PEG fragment (Fig. 5).To rule out the possibility that compound 7 binds at lower occupancy in the RT structure, a model with the PEG fragment at 50% occupancy and compound 7 in the adjacent position at 50% occupancy was tested but results in a very poor fit of the electron density (Fig. S4).
With the exception of compound 7, the structures agree well with the corresponding cryo-T structures.Interestingly, the RT structure solved from crystals soaked with compound 7 is distinctly different from that solved at cryo-T, where compound 7 is well fitted in the electron density at a position adjacent to the PEG position and not coordinated with the catalytic triad [Fig.5(c)].The fact that differences in ligandbinding features can be observed upon comparing structures solved at cryo-T and RT has been observed previously (Fischer et al., 2015).
3.3.1.Water molecules, B factors and loops.A significant difference between the sEH structures solved at RT and at cryo-T is the number of water molecules that can be modelled.
The water content is on average more than two times higher in the structures solved using conventional cryo-crystallography (Fig. 6).This is likely due to the fact that the elevated temperature increases disorder of water molecules that are not highly coordinated, as has been observed by others (Nakasako, 1999(Nakasako, , 2001)).Despite this difference in water content, we find that many well ordered water molecules in the interior of the catalytic domain are also conserved in the SX structures.
As expected, we observe higher B factors in the RT structures compared with the cryo-T structures.This difference can be attributed to several factors, including increased thermal motion and enhanced flexibility when the structure is not cryotrapped.The crystal lattice of microcrystals may be more resilient to conformational changes, thus allowing the dynamics of the protein to be observed.Radiation damage can also be manifested in the form of high B factors, and care must be taken to stay below the dose limit of 5 � 10 5 Gy when X-ray radiation induced damage is typically observed in data collected at RT (Schneps et al., 2022).In our experiments, we calculate the radiation dose to �4.5 � 10 4 Gy at each grid position of the chip (Bury et al., 2018).In a recent study of a cryptochrome using fixed-target SX, it was noted that the refined B factors correlated well with those predicted from molecular-dynamics simulations (Schneps et al., 2022), pointing towards the protein being in a fully solvated, functional state in the microcrystals.The structures of sEH display in some loop-regions at both temperatures, including those surrounding the active site.Interestingly, the B factors of the residues around the active site are significantly elevated in the RT structures relative to the cryo-T structures, suggesting a dynamic behaviour of the active site under ambient conditions that is masked in the frozen state (Fig. 7).This type of enhanced loop mobility of functionally important regions was reported recently in a lower-resolution SX study (Schneps et al., 2022), where it was suggested that information gained from analysing dynamic behaviour could be useful to identify key structural elements in a protein.
One of the most pronounced differences related to temperature in the structures of sEH is that several of the loops are more well defined in the RT models.Thus, it seems that it is possible to obtain more information on structural regions that are typically less ordered by use of SX.In our SX structures, we can confidently model the loop composed of residues 65 to 95 in the N-terminal domain, whereas this loop is either not possible to model at all in the cryo-T structures or only associated with partial electron density (Fig. 8).Interestingly, this loop is close to the putative phosphatase region of sEH (Cronin et al., 2003;Matsumoto et al., 2019;Va ´zquez et al., 2023).The fact that RT diffraction data can provide superior electron density maps of flexible regions has been noted previously in a study of a bacterial reaction centre using XFEL radiation (Dods et al., 2017).The reason for this temperature effect may be related to the fact that flash-cooling induces an artificial thermal heterogeneity in the structure (Halle, 2004).In contrast, it was recently shown that lipid tails and detergent molecules were not as well resolved at higher temperature (Ba ˚th et al., 2022).

Potential for optimization of the SX workflow
A serious limitation of how SX is used today is the timeconsuming nature of the method.In the study presented here, the collection of each dataset typically took 1-2 h.We see several possible avenues to explore in order to achieve a more efficient workflow.In the case of our model system sEH, the  speed of data collection could be enhanced by collecting sufficient -but not excessive -amounts of data.To investigate how the quality of the electron density map is affected by reducing the number of images, we solved the structure of sEH in complex with compound 5 using a subset of the data containing �5000 indexed diffraction images (to be compared with the �65 000 indexed images used for the structure presented in Table 1).The structure resulting from the subset of the data is of similar resolution (2.2A ˚based on the subset of data compared with 2.1 A ˚for the full dataset) and the ligand can be unambiguously placed in the active site electron density.Thus, it is clear that the very high redundancy, and thereby data collection time, of several of our datasets is not necessary and there is potential to reduce the data collection time by up to 5� using the current setup.Moreover, there are other sample delivery methods to be explored for efficient SX data collection (Zielinski et al., 2022).
A future development that would significantly reduce the data collection time is to make the fixed targets compatible with the robotics systems at the beamlines so that mounting and dismounting of the chips could be achieved without manual handling.
Further, over the last couple of years, automatic data processing pipelines have been successfully established at many SX beamlines.This development offers rapid feedback regarding, for example, resolution during the experiment, and also reduces the time that has to be invested into each structure after the data are collected.
Another limitation that is currently associated with the method is the lack of reliable ways to ship or transport the microcrystals at RT without damaging them.On-chip crystallization has the potential to overcome this problem and may be a future path towards enabling remote collection also of SX data as it would circumvent the need for manual crystal handling during the experiment.

Conclusions
We have developed a robust workflow for efficient production of microcrystals for SX experiments.This involves the use of a hybrid crystallization approach to screen for suitable conditions using low amounts of protein.The method allows easy evaluation of crystallization trials and is readily scalable to larger volumes.Hybrid crystallization was successfully used to find microcrystallization conditions for the pharmaceutical target sEH, and based on this, larger amounts of microcrystals were produced by batch crystallization.We were able to solve RT high-quality structures of sEH in complex with six different ligands by SSX using a fixed-target setup.In comparison with previously published cryo-T structures of the same protein-ligand complexes, we found that the ligands generally bind in a similar fashion in the RT and the cryo-T structures.However, we observed one example of a temperature-dependent difference in the binding of a lowmolecular-weight compound, information that may influence the chemical design in a drug discovery project.In addition, we observe differences in protein conformations and solvent molecules both in the vicinity of the active site and in the more remote regions.Thus, it may be beneficial to consider structural information at RT in a structure-based drug design project, in particular when working with weakly binding compounds.
Future studies will need to address challenges in working with poorly soluble ligands that cause precipitation, as well as issues with transporting and shipping crystals at RT.The most pressing challenge for SX from the perspective of pharmaceutical drug discovery is to tackle the slow turnover and labour-intense nature of this method.Indeed, there is potential to optimize the SX workflow by developing a more automated crystal handling process.If successful, SX could also become an attractive alternative for high-throughput drug discovery campaigns.

Figure 1
Figure 1 Structure of sEH.(Left) Structure of inhibitor-bound sEH solved at cryogenic temperature (PDB entry 5ake; O ¨ster et al., 2015).The active site with the bound compound 4 (purple) is located in the C-terminal domain.(Right) Zoom-in of the active site where the compound (purple) interacts with the catalytic triad (Asp335, Tyr383 and Tyr466).

Figure 2
Figure 2Microcrystallization and data collection workflow.Panel 1: production of crystal seeds from macro crystals using seed beads and vortexing.Panel 2: hybrid crystallization used for screening to find optimal batch crystallization conditions.Panel 3: batch crystallization with seeding and vortexing to induce nucleation.Panel 4: soaking of compound into the microcrystals.Panel 5: dispensing of soaked crystals onto the chip and sealing.Panel 6: fixedtarget raster grid-scan data collection on ligand-soaked crystals at a synchrotron source.

Figure 5
Figure 5Temperature-dependent difference in compound binding mode.The RT structure from crystals soaked with compound 7 is compared with the RT apo structure and the cryo-T structure in complex with compound 7. (a) In the RT structure of sEH obtained from crystals soaked with compound 7, the active-site density is best fitted with a PEG fragment (PDB entry 8qwg).(b) RT apo structure with a PEG fragment bound in the active site (PDB entry 8qvm).(c) In the cryo-T structure from crystals soaked with compound 7 (PDB entry 5ai8; O ¨ster et al., 2015), the compound is bound as well as a sulfate ion and a water molecule.For each structure, the 2F o F c electron density map is contoured at 1� (blue) and the F o F c electron density map at +3.5� (green).

Figure 6
Figure 6Comparison of water molecules modelled at RT and at cryo-T.(Left) RT sEH structure (in complex with compound 3, PDB entry 8qvk) with water molecules in red.(Right) Corresponding cryo-T structure (PDB entry 5ake) with water molecules in black.The central panel displays an overlay of the water molecules.

Figure 7 B
Figure 7 B-factor comparison.A B-factor comparison is shown for one RT structure (PDB entry 8qwg) and its cryo-T counterpart (PDB entry 5ai8).At the top, the RT (left-hand side) and cryo-T (right-hand side) sEH structures are displayed colored according to B factors.At the bottom, the average B factors of main chain atoms are plotted according to residue number (RT structure -orange, cryo-T structure -blue).

Figure 8
Figure 8Ordering of loops.An example of a loop (Pro65 to Ala95) with a higher degree of order in the RT structures compared with in the cryo-T structures.(a) 2F o F c electron density map associated with crystals soaked with compound 7 at RT (PDB entry 8qwg).(b) 2F o F c electron density map associated with the corresponding cryo-T structure (PDB entry 5ai8).The 2F o F c electron density maps are contoured at 1� (blue).

Table 1
Data collection, processing and refinement statistics.