research papers
Wholepattern fitting technique in serial femtosecond nanocrystallography
^{a}ARC Centre of Excellence for Advanced Molecular Imaging, School of Physics, University of Melbourne, Parkville, Victoria 3010, Australia, and ^{b}CSIRO Manufacturing Flagship, Parkville, Victoria 3052, Australia
^{*}Correspondence email: s.williams3@student.unimelb.edu.au
Serial femtosecond Xray crystallography (SFX) has created new opportunities in the field of structural analysis of protein nanocrystals. The intensity and timescale characteristics of the Xray freeelectron laser sources used in SFX experiments necessitate the analysis of a large collection of individual crystals of variable shape and quality to ultimately solve a single, average via indexing and merged into a threedimensional data set, as is done for conventional crystallography data. In this way, serial femtosecond crystallography data lie in between conventional crystallography data and powder diffraction data, sharing features of both. The extremely small sizes of nanocrystals, as well as the possible imperfections of their crystallite structure, significantly affect the diffraction pattern and raise the question of how best to extract accurate structurefactor moduli from serial crystallography data. Here it is demonstrated that wholepattern fitting techniques established for onedimensional powder can be feasibly extended to higher dimensions for the analysis of merged SFX diffraction data. It is shown that for very small crystals, wholepattern fitting methods are more accurate than Monte Carlo integration methods that are currently used.
Ensembles of crystals are commonly encountered in powder diffraction, but serial crystallography is different because each crystal is measured individually and can be orientedKeywords: protein nanocrystallography; peakshape analysis; wholepattern fitting; Xray freeelectron lasers; nanocrystals; SFX; protein structure; XFEL.
1. Introduction
Serial femtosecond Xray crystallography (SFX) (Chapman et al., 2011; Spence et al., 2012), which emerged with the commissioning of hard Xray freeelectron laser (XFEL) sources, provides a unique opportunity for modern biology to conduct structural analysis of proteins which have previously been inaccessible to study because of the extremely small size of crystals that they form (e.g. submicron). The capabilities of the SFX approach have been successfully tested in recent studies of submicron and nanoscale protein crystals (Chapman et al., 2011; Boutet et al., 2012; Johansson et al., 2012, 2013; Koopmann et al., 2012; Redecke et al., 2013; Aquila et al., 2012; Liu et al., 2013; Demirci et al., 2013; Kupitz et al., 2014).
The SFX technique involves illuminating a stream of randomly oriented protein crystals of various sizes and orientations by an extremely bright and ultrashort (tens or hundreds of femtoseconds) XFEL source and properly merging the obtained diffraction data. SFX experiments are performed in this manner due to the destructive nature of the XFEL source for which a single exposure can be expected to cause the disintegration of a nanocrystal (Neutze et al., 2000). Diffraction information is obtained from different orientations of the of the target structure through continual replenishment of individual crystal samples in the XFEL beam. Whilst some characteristics of single crystals can be obtained from processing individual SFX diffraction patterns, the solution of a threedimensional requires the processing of large numbers of SFX diffraction patterns. The solution obtained is then an average crystal unitcell structure found from data that are intrinsically based on distributions of both crystal sizes and qualities.
The SFX approach has a number of problems that have not been so critical to protein crystallography until now. In an SFX experiment, a single crystal effectively stands still during Xray illumination due to the femtosecond timescale of an XFEL pulse. A collected diffraction pattern may then be thought of as a diffraction `snapshot' that represents the diffraction of a single XFEL pulse from a single particle of finite size and a unique orientation. Furthermore, the small size of the protein crystals illuminated by the XFEL source can create broad intensity distributions around Bragg reflections (Yefanov et al., 2014). Since there is no time for any rotations of the sample during the XFEL pulse, only partial information about the crystal shape transform is recorded on the diffraction pattern. Consequently, the shapes of the observed Bragg reflections may vary significantly from shot to shot and within a single image (White et al., 2012). The SFX data set consists of twodimensional diffraction patterns comprising partially recorded information from different crystals. To resolve these issues and to obtain the structurefactor moduli of Bragg reflections, which contain encoded molecular structural information, the SFX approach has relied on the Monte Carlo integration method (Kirian et al., 2010, 2011; White et al., 2012; White, 2014), in which a large number of diffraction patterns are separately analysed. For each diffraction pattern, the recorded intensities are summed within a fixed integration volume around each Bragg reflection and the obtained integrated intensities are averaged over all `snapshots'. The choice of the integration area is critical for the integration method to work and the accuracy of this procedure determines the accuracy of further structural analysis. Current crystallographic programs, such as MOSFLM (Leslie & Powell, 2007), CrystFEL (White et al., 2012) or Cheetah (Barty et al., 2014), use either circular or rectangular integration areas for the analysis. It has also been suggested that proper integration of diffracted intensities within the Wigner–Seitz cell around each Bragg reflection could be used (Kirian et al., 2010, 2011).
The approaches currently used in protein crystallography to extract structurefactor moduli from the diffraction pattern (Leslie & Powell, 2007; Kabsch, 2010; Kirian et al., 2010; White et al., 2012; Barty et al., 2014) rely on the segregation method, i.e. the diffraction pattern is considered as a discrete set of completely isolated Bragg reflections. In this case, the intensity distribution within a predefined region around a given Bragg reflection is used to extract structurefactor moduli. Here we present a wholepattern fitting technique that uses a continuous description of the merged diffraction data. This holds some similarities with profilefitting methods used for singlecrystal data, such as MOSFLM (Leslie & Powell, 2007) and XDS (Kabsch, 2010), and those recently introduced for SFX data analysis, such as nXDS (Kabsch, 2014). A key difference is that the latter approaches rely on the segregation method for the extraction of structurefactor moduli with the assumption that Bragg reflections can be separated and integrated. These approaches are based on the scaling of learned peak shapes found from strong peaks that are assumed to be isolated, rather than using analytical expressions to model a continuous intensity distribution. Moreover, the integration approach relies on the statement that the resulting shapes of Bragg reflections on the SFX diffraction pattern are governed by the averaged crystal shape transform only (Kirian et al., 2010; White et al., 2012), which is the same for all Bragg reflections. Thus, assuming a defectfree structure of protein nanocrystals, the structurefactor moduli can be extracted from SFX data by the integration of diffracted intensities within a predefined region around all Bragg reflections. We will later refer to this as the `integration approach'.
We have recently shown (Dilanian et al., 2013), however, that the size and the quality of individual protein nanocrystals illuminated by an Xray source during the SFX experiment significantly affect the resulting SFX diffraction pattern. Given the asymptotic behaviours of the Bragg reflections due to the size distribution of the nanocrystals and the scattering in interBragg regions due to the size of individual crystals, the assumption that the Bragg reflections are isolated cannot always be satisfied. The smaller the nanocrystals and the bigger the unitcell parameters of the protein crystal, the stronger is the influence of asymptotes of the nearest Bragg reflections on the intensity distribution of a given Bragg reflection. Moreover, the large surfacetovolume ratio of protein nanocrystals leads to significant contributions from surface effects, such as structural disorder, impurities and lattice distortions near the surface of the crystal etc. (Feher & Kam, 1985; Grant & Saville, 1994; Caylor et al., 1999; Malkin & Thorne, 2004), to the diffraction pattern. In such cases the shape of Bragg reflections will not be exclusively governed by the averaged crystal shape transform and, therefore, will not be identical for all reflections, varying with the scattering vector (Dilanian et al., 2013). Consequently, restriction of the integration areas may lead to an incorrect estimation of the corresponding structurefactor moduli.
In this paper, we present an approach to the extraction of the structurefactor moduli of Bragg reflections from SFX data which resolves the issue mentioned above. A key idea in this approach is the treatment of the merged SFX diffraction data set as a continuous function of the scattering vector, q, and not as a discrete set of Bragg reflections. Such a treatment is demonstrated in the threedimensional merging of whole twodimensional SFX diffraction patterns by Yefanov et al. (2014). This involves the mapping of whole twodimensional diffraction patterns (according to the orientations of individual particles) to threedimensional q space for further analysis. The result is a diffraction data volume comprised of an ensemble of particles of different sizes and structural qualities that retains the dimensionality of the reciprocal As indicated by Yefanov et al., it is possible to take into account fluctuations in the incident pulse intensity and beam convergence during an SFX experiment by weighting intensities during merging based on singleshot spectra. From this point of view, there is a similarity in the formation of a merged SFX diffraction data set and a powder diffraction pattern, where the shapes of Bragg reflections are similarly formed by a collection of independent scatterers of varying characteristics. Differences exist in that individual SFX diffraction patterns may be first processed before merging and that the merged data need not be collapsed into a function of the diffracted intensity with respect to the scattering vector magnitude (i.e. in one dimension) due to the ability to estimate the orientation of individual crystals. Instead, twodimensional or threedimensional merged SFX diffraction data sets can be formed. For a sufficient number of independently scattered particles, the shapes of the resulting Bragg reflections in the merged SFX diffraction data volume, as well as the asymptotic behaviours of their tails, will be governed by the statistical properties of distribution functions, which are characterized by the variations of particle size and structural imperfections (Suortti & Jennings, 1977; Suortti et al., 1979; Young & Wiles, 1982), and not by properties of individual crystallites. Consequently, the structurefactor moduli may be extracted from the merged SFX data by fitting the whole SFX diffraction pattern using analytical peakshape functions, which are defined over the entire range of q space, rather than only within restricted areas around individual Bragg reflections. Moreover, the fitting approach allows us to incorporate characteristics of Bragg reflections and the scattering in interBragg regions into the analysis via the adjustable parameters of the peakshape function, such as unitcell parameters, width of the Bragg reflection or asymmetry of the reflection. Here, we demonstrate that the fitting approach may provide more accurate and robust results in extraction of structurefactor moduli from protein nanocrystals compared to the integration approaches developed so far.
2. General considerations
We base our analysis on the consideration of finite crystals and on the assumption that the merged SFX diffraction pattern collected from a stream of nanoscale protein crystals of various sizes and qualities is a continuous function of the scattering vector, q. We start with the assumption that the scattering factor for a finite crystal can be expressed as
where is the kth Bragg reflection and is a series of functions centred at Bragg positions, which arises from the finite We will refer to as the structurefactor modulus of the kth Bragg reflection here. Equation (1) is appropriate for crystals with whole unit cells (Ino & Minami, 1979) and can account for lattice imperfections like strain or lattice defects. For a perfect finite lattice, is identical for all values of k, but in general it can vary for different Bragg peaks. We note that in the nanocrystallography literature (Kirian et al., 2010, 2014; Spence et al., 2011), the scattering factor for a finite crystal is commonly expressed as a product of continuous unitcell and lattice scattering factors, which for an ideal crystal with whole unit cells can be written in the form of equation (1) (see Appendix A). Equation (1) cannot describe models of the crystal surface that prevent the continuous scattering factor of the , from being completely characterized by discrete samples of (Ino & Minami, 1979), such as the presence of incomplete unit cells on the crystal surface.
is the Bragg position of theAssuming the validity of equation (1), the diffracted distribution is given by
The second summation term in equation (2c) is typically neglected (James, 1954; Hosemann & Bagchi, 1962; Guinier, 1963; Ino & Minami, 1979). For sufficiently large crystals, equation (2c) is well approximated by the first summation term solely (James, 1954; Guinier, 1963). If the crystal surface is modelled by a random shift of the lattice with respect to the crystal's centre, Ino & Minami (1979) show that the second summation can also disappear for small crystals. We note, however, that this model also requires a continuous dependence on the scattering factor, and is not consistent with the underlying assumptions made within equation (1), such as the presence of whole unit cells. In fact, the second summation contains phase information that is targeted by new methods for directly phasing SFX data (Spence et al., 2011). In light of the progress in the direct phasing of finite crystals (Elser, 2013; Liu et al., 2014; Kirian et al., 2014, 2015), it is unlikely that the second summation is entirely absent for SFX data. Here, however, our goal is to obtain an improved estimate of the structurefactor moduli by accurately modelling the first summation of equation (2c) and not to pursue the additional phase information contained in the second summation term. We do so without invoking the physical assumptions of Ino & Minami (1979) of random shifts. We test the validity of ignoring the second term via simulation in §4. Further work may consider the contribution of this term more rigorously.
The model used here for the threedimensional diffracted intensity distribution, merged from a collection of N finite crystals, can then be written as
This holds similarities to the onedimensional intensity distribution used in the profilefitting analysis of powder diffraction data (Rietveld, 1967), which can be given by
The above expression represents a wholepattern fitting scheme in which a set of sampling points within a diffraction pattern, , and an associated set of peak profile functions, , are chosen to evaluate the continuous function, . Based on the aforementioned similarities of SFX and powder diffraction data, we propose that the average shapetransform function contained in our model, , may also be modelled by continuous peakshape functions, . For the proposed wholepattern fitting analysis, we define the modelled intensity distribution of the merged SFX diffraction as
where the summation is performed over all Bragg reflections, contributing to the intensity distribution, , for a given scattering vector, q. is the structurefactor modulus of the kth Bragg reflection, which we seek to accurately extract. It should be noted that in the modelled continuous intensity distributions [equations (3), (5)], all contributions from individual crystals are given equal weight.
With this model, each point of the resulting SFX diffraction pattern is affected by all Bragg reflections via the corresponding peakshape functions, . The intensity variation between neighbouring Bragg reflections will be governed by the collective contribution of both reflections. The degree of this `overlapping' of the peakshape functions depends on averaged size of the crystallites illuminated by Xrays, the size of the and the structural quality of protein crystals. From this point of view, Bragg reflections can be considered as isolated only if the corresponding peakshape functions tend to zero within the Wigner–Seitz cell region around each Bragg reflection, which allows one to disregard the scattering in the interBragg regions. In this case, the aforementioned segregation method provides correct estimations of the structurefactor moduli. In general, however, the collective contribution of Bragg reflections to the intensity distribution cannot be neglected and the scattering in interBragg regions should be included in the analysis of the diffraction pattern (Dilanian et al., 2013). Thus, consideration of the diffraction pattern as a continuous function of the scattering vector [equations (3), (5)], which forms by a collective contribution of all Bragg reflections via the corresponding peakshape functions, may provide a more comprehensive approach for the analysis of diffraction data from a stream of nanoscale protein crystals.
The presence of whole unit cells has been assumed in the simulations presented here. The effects of partial unit cells present on the surface of protein nanocrystals for structural analysis have been considered by others (e.g. Liu et al., 2014; Kirian et al., 2014). Because of the finite size of nanocrystals, partial unit cells can influence intensity distributions between Bragg positions and introduce further ambiguities regarding the definition of the for a finite crystal. While the results of simulations presented here have only contained whole unit cells, we will consider the presence of partial unit cells in future work for which modifications may be introduced to the model used for the scattering factor for a finite crystal.
Two different schemes for extracting structurefactor moduli are considered here – an integration approach, based upon the consideration of discrete sets of Bragg reflections, and a wholepattern fitting approach, based upon a continuous treatment of the diffraction data. The choice between these approaches depends on the average size of protein crystals, their size distribution and imperfections in et al., 2013).
According to our previous analysis, the consideration of the diffraction pattern as a continuous function of the scattering vector is crucial when the fraction of unit cells adjacent to the surface of the protein crystal exceeds 10% of the total number of unit cells. This suggests an effective size limit of 1 µm for a protein crystal with an average unitcell parameter of 100 Å (DilanianIntegration and wholepattern fitting approaches are used here to extract structurefactor moduli from simulated diffraction patterns in order to evaluate the effectiveness and accuracy of both approaches in the analysis of SFX data from protein nanocrystals of various sizes. In the case of the fitting approach, equation (5), we extend wholepattern fitting procedures developed for the analysis of powder diffraction data (Le Bail et al., 1988) from fitting of onedimensional intensity profiles to higher dimensions. Comparative analysis is presented in §4.
3. Fitting procedure
The fitting procedure outlined here is conducted in a similar manner to that described by Le Bail et al. (1988). To begin, knowledge of the (i.e. unitcell parameters and spacegroup symmetry) and experimental geometry are used to predict Bragg peak locations in the modelled intensity distribution [equation (5)]. These positions can later be refined. Initial structurefactor moduli, , can be estimated by Wilson's statistics (Giacovazzo, 2011). The initial modelled distribution [equation (5)] is calculated continuously using the estimates of structurefactor moduli and peakshape distributions, based upon the chosen form of the peakshape function and peakshape parameters that are input by the user.
The peakshape parameters are refined in each iteration of the wholepattern fitting procedure by minimizing the errorcost function,
where is the measured intensity and is the calculated continuous intensity distribution [equation (5)]. The summation is performed over all measured points of the diffraction pattern, rather than being assessed only at Bragg positions.
As in Le Bail analysis (Le Bail et al., 1988), the structurefactor moduli are estimated at the conclusion of each parameter cycle and are fed in to the next cycle as the input structurefactor moduli. We propose the estimation of the structurefactor moduli based upon the form of modelled intensity distribution [equation (5)]. This can be evaluated at Bragg locations to give
where I is a column matrix with elements, I_{j}, given by the observed peak intensity in the merged diffraction data set for the jth Bragg reflection and F is a column matrix with elements, F_{j}, given by the squared structurefactor moduli, . The square matrix, P, has elements, P_{ij}, given by the contribution from the peakshape distribution for the ith reflection relative to the position of the jth reflection, i.e. . The P_{ij} elements are evaluated based on the current peakshape parameters, . After each cycle, the structurefactor moduli can then be estimated on the basis of the current peak shapes through inversion of the P matrix, using in our case the Gaussian elimination method, to find the F column matrix elements that satisfy equation (7) for the merged diffraction intensity distribution. In this manner, structurefactor moduli are iteratively estimated at the end of each fitting cycle.
The accuracy of the fitting process is monitored by two R factors, R_{wp } and R_{B}. The standard definitions of these factors are used; see the comprehensive review by Hill & Fischer (1990). After a good fit between the modelled intensity distribution and merged intensity distribution is achieved, the current structurefactor moduli are extracted. Wilson statistics can then be employed to include scaling and thermal correction factors to the extracted structurefactor moduli (Giacovazzo, 2011).
3.1. Peakshape function
In the analysis presented here, we have chosen to model the peak shapes formed from the merged diffraction data with normalized pseudoVoigt functions. This choice was motivated primarily by the pseudoVoigt function's widespread use in the analysis of powder diffraction data collected from crystal samples of varying size, shape and quality (e.g. Young & Wiles, 1982; Langford, 1999) and by the function's considerable flexibility in form. The modelling of peak shapes with analytical functions such as the pseudoVoigt function also improves computational speed and efficiency. It should be noted, however, that the general scheme presented here is based on the broader form of the modelled intensity distribution expressed in equation (5), for which alternative peakshape distributions could be readily incorporated. Improvements to this approach might be made with greater theoretical consideration of the expected form of peakshape distributions.
In this case, the calculated intensity distribution of the kth Bragg reflection from equation (5) can be written as
where represents normalized pseudoVoigt functions, are adjustable peakshape parameters, and M is the total number of adjustable parameters. In general, each of the adjustable parameters, , can be represented by an analytical function of the scattering vector, q, which takes into account possible variations of shapes of individual Bragg reflections with respect to the scattering vector.
We define an extended pseudoVoigt peakshape function in the following way:
where . and are normalized Gaussian and Lorentzian functions, respectively, with widths defined by and . The parameters determine the weighting of both functions in the peak shape.
4. Results and discussion
In this section we show the results of the extraction of the structurefactor moduli of the Bragg reflections from merged twodimensional diffraction patterns using the wholepattern fitting approach. The following is presented to demonstrate the feasibility of the extension of established powder
techniques to higherdimensional SFX diffraction data. While the simulated diffraction data correspond to twodimensional crystallographic planes, we expect this technique to be readily extendable to threedimensional diffraction data.4.1. Simulations
The sugarbinding domain of langerin protein with the F241L mutation (Chabrol et al., 2015; PDB entry 4AK8 ) was used as a test case for simulation studies. This protein has P4_{2} spacegroup symmetry and unitcell parameters of a = b = 79.959, c = 90.419 Å. Simulations were performed to test the applicability of the wholepattern fitting method to data merged from crystals of varying mean size. This was investigated by simulating several sets of needlelike nanocrystals of langerin with the numbers of unit cells in each crystal randomly generated from lognormal distributions. Both the shape and size of individual crystals were allowed to vary by sampling from independent size distributions that were taken to correspond to the orthogonal dimensions of a crystal. Finitelattice transforms [see equations (16) and (17) in the Appendix and Dilanian et al. (2013), for details] were calculated along the [100] and [010] directions of the reciprocal crystal lattice.
The mean numbers of unit cells were varied between sets of simulated crystals as follows. The mean number of unit cells in the Y direction, , was varied from 60 unit cells progressively down to 30 unit cells in increments of 5, while the mean number of unit cells in the X direction, , was fixed at 10 unit cells. Each generated set involved a total of 2000 nanocrystals. The orientations of simulated crystals were constrained to diffract into a selected crystallographic plane, the (hk0) plane. A significant increase in the total number of crystals would be required to fill a complete threedimensional diffraction volume.
Within each simulated set of crystals, individual diffraction patterns were calculated for each crystal on a 1024 × 1024 array [see Dilanian et al. (2013) for more details] using the published atomic positions and Debye–Waller factors (PDB entry 4AK8 ). The crystal set corresponding to the smallest average size was also calculated on a 340 × 340 pixel array. Knowledge of the such as the unitcell parameters and symmetry operations, was employed during the calculation of individual diffraction patterns and the subsequent merging. All diffraction patterns were calculated in the (hk0) crystallographic plane to 5.0 Å resolution. Additional crystal disorder effects were not included in this simulation study.
Merging of the individual diffraction patterns was performed for each collection of crystals to create seven twodimensional diffraction patterns defined by varying mean crystal dimensions. Since we are only considering the twodimensional diffraction patterns in the (hk0) crystallographic plane, the dimension of the nanocrystals along the [001] direction is not relevant in this case. The proposed approach, however, can be applied to any arbitrarily oriented reciprocallattice plane and extended to the fitting of threedimensional diffraction data.
4.2. Analysis
The accurate and robust extraction of structurefactor moduli is the primary objective of the presented wholepattern fitting technique. Molecular structural information is encoded within these values and the structural analysis of protein crystals for diminishing crystal size requires suitable methods with which to read out such structural information. Structurefactor moduli estimated from wholepattern fitting, , are compared here both to the values, , calculated from published structural data (PDB entry 4AK8 ) and also to the values, , found from the integration approach. The respective relative errors, = and = , of the extracted structurefactor moduli are assessed for both the wholepattern fitting and the integration methods. It is shown here that the accuracy of extracted structurefactor moduli via the integration of merged diffraction data is sensitive to the chosen integration area, whereas the values found from the presented wholepattern fitting approach do not require integration and exhibit improved accuracy for diffraction data merged from small average crystal sizes.
In this investigation, a variety of integration areas were trialled for the extraction of structurefactor moduli via integration for the test diffraction data from the protein langerin. Scaling and thermal correction factors were estimated and applied for all sets of structurefactor moduli found from each integration area. This was performed using Wilson statistics (Giacovazzo, 2011), as with the structurefactor moduli extracted via wholepattern fitting. The average relative error, , was assessed, where N_{p} is the number of extracted structurefactor moduli. Figs. 1(a)–1(g) show the average relative errors from integration for the seven merged diffraction patterns for a variety of integration areas. The integration area was defined by a circle with the radius, , varied from zero (the Bragg position) to half of the distance between the nearest Bragg reflection, , where is the reciprocallattice parameter, which corresponds, in our case, to the boundary of the Wigner–Seitz cell along the [100] (or [010]) direction of the The average error found from wholepattern fitting is also indicated in Figs. 1(a)–1(g). These values are represented by solid lines given that the extracted structurefactor moduli are independent of integration area with this technique. It should be noted that in each instance and in both the wholepattern fitting and the integration methods, a fixed set of several very weak reflections has been excluded from the assessment of average relative errors. These were excluded due to the weighted errors being considerably higher in these cases and thus having a significant impact on the overall averages. The problem of weak reflections will be discussed later in this section.
The wholepattern fitting approach is shown in Fig. 1 to improve the accuracy of extracted structurefactor moduli for the smallest average crystal sizes. This is shown in Figs. 1(a)–1(b) and, most significantly, in Fig. 1(h), for which the smallest average crystal size is represented with fewer pixels simulated. Similar accuracy is achieved by both approaches for larger average crystal sizes, shown in Figs. 1(c)–1(g). The relative success of the integration method shown in Figs. 1(a)–1(g) can be understood from the size of the integration areas: the smallest average errors were found from the smallest integration lengths. It follows that, as the integration length is decreased, the extracted values approach the values of the heights of the Bragg reflections. Given that additional sources of structural disorder were not included in the simulated diffraction data, it is to be expected that the peak heights of well resolved, highsignal Bragg reflections would provide accurate structurefactor moduli. However, both the number of pixels between adjacent Bragg positions and the signaltonoise ratio of the diffraction data would often be lower for experimentally collected SFX data. The average errors found via integration indicate that intensity contributions from neighbouring peaks can quickly begin to contribute as the integration area is expanded. The accuracy of the wholepattern fitting approach does not exhibit this sensitivity.
The wholepattern fitting approach allows peak contributions to extend continuously throughout the modelled diffracted intensity distribution. Limited integration or calculation areas are not required for the extraction of structurefactor moduli, as outlined earlier. As in powder diffraction Le Bail analysis, the wholepattern fitting approach allows for individual reflections to be isolated and untangled from the contributions of nearby Bragg reflections. Average relative errors found [Figs. 1(a)–1(g)] indicate the stability of structurefactor moduli extracted by the wholepattern fitting formulation. This approach is shown to produce an accuracy comparable to the integration method using optimal integration areas, while the results from integration decrease in quality for varying integration area. Improved accuracy from wholepattern fitting for diffraction data merged from the smallest average crystal sizes is also evident. This favours the use of the wholepattern fitting approach particularly for the analysis of data from small average crystal sizes in SFX experiments, with stable results found that are independent of calculation areas.
An average of slightly over 30 pixels was held between the positions of neighbouring Bragg reflections due to the large number of pixels used in this simulation study. This provided the capacity for detailed distributions of peak shapes to test the accuracy of the extended pseudoVoigt function and of the wholepattern fitting approach. However, it is acknowledged that these pixel numbers are not necessarily realistic in the case of the analysis of highresolution SFX experimental data. To test this, simulations were performed of the diffraction pattern on a 340 × 340 pixel array to 5.0 Å resolution for a collection of 2000 nanoscale crystals of the test protein with mean dimensions of and unit cells along orthogonal realspace dimensions. In this case, as few as 9–10 pixels exist from Bragg peak to Bragg peak. Both the wholepattern fitting and the integration approaches were used to extract structurefactor moduli. A selection of extracted values and corresponding relative errors, , is provided in Table S1 (in the supporting information). Within Table S1, several weak reflections have been highlighted that have been excluded from the calculated average error, . The structurefactor moduli values shown in Table S1 are provided with the optimal integration area for the integration method. The dependence of the average relative error upon these areas was also tested and is shown in Fig. 1(h). It is apparent that the wholepattern fitting approach was able to extract structurefactor moduli to an accuracy similar to the earlier and more pixeldense case. Accuracy was significantly diminished for the structurefactor moduli extracted via integration, however. The dependence on integration area was also further exacerbated. This can be expected due to the fineness of integration area choices being limited by the fewer number of pixels present. It is worth noting that the ability of the integration method to approach the peakheight values of the Bragg reflections is hindered by the lower density of pixels and that the accuracy of extracted structurefactor moduli is poorly affected as a result. The wholepattern fitting approach does not appear to exhibit the same sensitivity. In cases where the peak heights of Bragg reflections cannot be expected to be reliably extracted for structurefactor moduli estimation, wholepattern fitting techniques may be valuable.
4.3. Peakshape parameters
Peakshape parameters found from the fitting procedures are given in Table 1 for each of the simulated diffraction patterns. No prior knowledge of the protein structure was used in fitting procedures, excluding the unitcell parameters for the estimation of accurate Bragg peak positions. Fitting procedures were restricted to a single quadrant of the (hk0) crystallographic plane given that, in this case, the tetragonal symmetry of the test protein structure can be exploited. Selected sections in the (hk0) crystallographic plane of the simulated and fitted diffraction patterns are shown in Fig. 2 for the sets of generated crystallites with the largest average unitcell dimension, . An example of cross sections of the simulated and fitted diffracted intensity distributions is shown in Fig. 3. Fig. 3 contains cross sections of the simulated and fitted intensity distributions for the smallest average unitcell case, , through Bragg reflections along the (h, 14, 0) crystallographic direction. Fitting of the simulated diffraction patterns was achieved with R factors of R_{wp} = 0.17, R_{B} = 0.18 and R_{wp} = 0.15, R_{B} = 0.16 for the cases with the smallest, , and the largest, , average unitcell dimensions, respectively.

The set of mixing parameters, , obtained from the fitting procedures, Table 1, demonstrates the dependence of the character of the peakshape description upon the average size of the collection of crystallites. Put differently, diffraction patterns simulated with smaller average crystal sizes obtained smaller average mixing parameters during fitting procedures. This indicates a trend of greater contribution from the Lorentzian component in the fitted pseudoVoigt functions for diffraction data merged from crystals of smaller average size. This might be interpreted as indicating more diffuse scattering in these cases, with less localized distributions of intensities found to be fitted around Bragg reflections.
Throughout the fitting process, inversion of the P matrix was employed to estimate structurefactor moduli according to current peakshape parameters. The condition number can be estimated by the ratio of the largest to the smallest singular value of the P matrix to test whether the system of linear equations, equation (7), is well conditioned. Condition numbers calculated from final peakshape parameters were found to be in the range of 1.1–1.2, indicating that the formulation used was well conditioned in the cases considered here, allowing the simple Gaussian elimination method to be used for matrix inversion.
4.4. `Overlapping' Bragg peak contributions
Each individual term in the summation of equation (5) represents the intensity distribution for a given Bragg reflection untangled from the nearby Bragg reflections. The intensity distribution around a given Bragg reflection in the merged diffraction pattern can be represented by
where is the modelled intensity distribution for the kth peak and is the combined contribution of other Bragg peaks around the Bragg position, . This second term defines the contribution of surrounding Bragg reflections to a given Bragg reflection in our modelled intensity distribution. The magnitude of this term determines the extent to which the Bragg reflections can be considered as isolated. The combined intensity distribution, , is used to model the total intensity distribution from the kth peak during wholepattern fitting.
In order to evaluate the influence of the nearestneighbour reflections on the intensity distribution of the Bragg reflection, we calculated the first and the second terms of equation (10) using the peakshape parameters obtained via wholepattern fitting. The contribution of two terms was calculated for various radial distances, , from the Bragg position of a selected reflection. Similarly to the previous analysis (Fig. 1), the radial distance varied in the range . Only the intensities exceeding 2% of the peak intensity of the corresponding Bragg reflection were considered in the calculations. The analysis was performed on two merged diffraction patterns generated from a collection of 2000 nanoscale crystals of the test protein with mean dimensions of , , and , , respectively. As one can see from Fig. 4, the contribution of the second term increases with the distance from the Bragg position of the selected reflection, quickly reaching 10% (or more) of the total intensity. The obtained results are in a good agreement with results shown in Fig. 1, in that the wider the integration area, the stronger the influence of the nearest Bragg reflections on the intensity distribution of the selected reflection and, consequently, the bigger the error in the determination of the corresponding integrated intensity. The contribution of the second term of equation (10) to the total intensity of the Bragg reflection can differ for different reflections. In particular, the influence is negligible when a strong Bragg reflection is surrounded by weak reflections. Conversely, the intensity distribution of a weak Bragg reflection surrounded by very strong reflections will largely be determined by the contributions of surrounding reflections.
This may be the cause of some difficulties present in both of the extraction methods considered here, which make the accurate estimation of structurefactor moduli for several weak reflections problematic. Further development of the wholepattern fitting method may be required for the treatment of weak reflections, together with the use of more robust and numerically stable algorithms for less well conditioned P matrices. In the current approach, all Bragg reflections are fitted simultaneously with either common peakshape parameters or common dependences of the peakshape parameters on the scattering vector. The weak reflections that were poorly estimated were often surrounded by much stronger reflections; this may have affected the accuracy with which these sets of reflections could be modelled with this approach. Further extensions could be made by performing the individual fitting of small numbers of selected weak reflections, following wholepattern fitting procedures. A similar approach is present in powder known as partial profile relaxation (Izumi, 2003). This would be expected to increase the accuracy of structurefactor moduli for weak reflections by allowing those with nearby strong reflections to have independent peakshape parameters.
5. Conclusion
Presented here is an approach that builds upon an established analysis technique in powder diffraction, the wholepattern fitting method (Le Bail et al., 1988). The impetus for this can be seen in the similarities between powder diffraction and merged SFX data sets. In both instances, peakshape distributions are formed by the shape, size and disorder characteristics of a large set of independent scatterers. It is shown here that appropriate extension of the wholepattern fitting technique can be used to closely model SFX diffraction patterns and to extract integrated intensity information of Bragg reflections. This follows the work of Dilanian et al. (2013) in the use of a continuous description of the diffraction pattern obtained from a distribution of protein nanocrystals.
This analysis indicates that the wholepattern fitting method is a feasible approach for the extraction of intensity information from SFX data. Parameters obtained from fitting procedures show some dependence on the average size of contributing crystals. Flexibility is provided in the form of peakshape distributions to fit merged diffraction patterns from crystals of varying mean sizes using smooth analytical functions adopted from the analysis of powder diffraction data. Further improvements in this approach might be made with closer consideration of the theoretical basis of the average shapetransform distributions formed by collections of nanocrystals.
The strength of the wholepattern fitting approach in isolating Bragg reflections is particularly desirable in cases where substantial interBragg scattering occurs – such as when unitcell parameters are large yet the dimensions of the contributing crystals are small. This supports the application of the wholepattern fitting approach in the analysis of highresolution experimental SFX data from large collections of protein nanocrystals.
Several important aspects of the SFX approach were outside the scope of the analysis presented in this article. We have analysed the defectfree structure of the protein crystal. This allowed all peaks to be fitted with the same , parameters. In general, each of the adjustable parameters and can be represented by an analytical function of the scattering vector, q, to take into account effects of the structural disorder on the diffraction pattern of protein nanocrystals. The unitcell parameters may also vary from crystal to crystal, or even within one crystal. This will lead to a shift of Bragg reflections from ideal positions. All of these problems are well known in singlecrystal and powder diffraction crystallography (e.g. Ungár & Gubicza, 2007; Palosz et al., 2003) and can further be incorporated into the analysis of SFX data. In our analysis we considered all the molecules which form incomplete unit cells on the surface of the protein crystals as independently scattered objects. In this case, molecules from incomplete unit cells contribute only to the background scattering (Welberry, 2004). This assumption is incorrect, however, for extremely small crystals, comprised of several molecular clusters, when the influence of such molecules on the diffraction pattern is strong (Chen & Millane, 2013). It should also be noted that the definitions of `crystal lattice' and the `unit cell' are not entirely clear in this case and a separate analysis of this situation is required.
APPENDIX A
The scattering factor for a finite crystal
The scattering factor for a nanocrystal is often expressed as the product of a finite lattice transform and unitcell scattering factor. Here we describe the alternative expression used throughout this work [equations (1) and (15)] and show its equivalence to the former formulation [equation (16)] for the ideal case of a finite crystal composed solely of whole, identically ordered unit cells. This should be viewed simply as an example presented for an idealized case. Discrepancies in the scattering factor formulations can be expected to arise for ensembles of real crystals for which the diffracted intensity distributions will be shaped by variations in factors such as size, shape, disorder and strain.
The formulation of the diffracted ) and Ewald (1940), the electron density of a finite crystal, , can be defined as
from a finite crystal is briefly reviewed here. Following representations of Patterson (1939where r is a vector in real space and is the electron density of an effectively infinite such as in macroscopic crystals used in conventional crystallography. The shape of the crystal is denoted and can be written as
Here, V_{C} is the volume of the crystal and determines its boundary. We adopt the further definition of Guinier (1963) that is bounded such that whole unit cells are contained within the crystal volume. The scattering amplitude, F_{C}(q), can then be represented as
where is the scattering amplitude from an effectively infinite crystal and is the Fourier transform of . Ignoring nuclear motion, the scattering amplitude from an infinite crystal can be expressed as
where is the scattering factor of the
sampled at the Bragg positions, , with the Dirac delta function distribution, . This provides the scattering factor for the finite crystal aswhich is identical to the form of the scattering factor provided in equation (1). Recent work (e.g. Kirian et al., 2010, 2014; Spence et al., 2011) has used a formulation of the scattering factor for a nanocrystal that can be expressed as
where is the continuous scattering factor of the
and is the Fourier transform of the lattice function of the nanocrystal, . The lattice function in real space can be written as a finite series of Dirac delta functions located at the corner of each occupiedwhere N denotes the number of unit cells in the nanocrystal. We seek to demonstrate here that both formulations [i.e. equations (15) and (16)] are equivalent for the ideal case considered in this appendix.
The scattering factor of the
, is a bandlimited function and thus can be reproduced with a minimal sampling rate given by the Shannon–Nyquist sampling, which is given here by the Bragg locations. The scattering factor for the can consequently be expressed aswhere is the Fourier transform for the shape function of the
centred on the Bragg locations, . The unitcell shape function can be defined aswhere V_{U} is the volume of the and defines its boundary. Knowledge of at the Shannon–Nyquist sampling points, , is sufficient to interpolate to find at other positions using the unitcell shape transform, .
The scattering factor for the nanocrystal given in equation (16) can then be expressed as
The product, , can be written as
In our approach, the following relationship is used:
Equations (21) and (22) are not equal for general expressions due to convolution not holding associativity properties with multiplication. These equations are equivalent in this work, however, due to the characteristics of the particular functions considered here. This can be shown by considering the inverse Fourier transform of each expression. The inverse Fourier transform of equation (21) is given by
where indicates the inverse Fourier transform operation and K is an integer given by the dot product of a lattice point vector, , and a reciprocallattice point vector, . The reciprocal relationship of the points and arises from the bandlimited nature of the scattering factor of the . Given that K is an integer, the phase factor is equal to one. This provides the inverse Fourier transform of equation (21) as
This is equal to the inverse Fourier transform of equation (22), given the following:
The equality of equations (24) and (25) arises from the characteristics of the lattice function, , and from the sampling points in . These characteristics allow the convolution operator to be applied equivalently in equations (21) and (22). The order of the application of convolution operations between shape functions and unitcell electrondensity distributions can create differences in models for the scattering factor for a nanocrystal when the presence of incomplete unit cells is considered (Ino & Minami, 1979; Beyerlein, 2011). Alternative scattering factor models have been evaluated in powder diffraction studies (Ino & Minami, 1984); further consideration of the influence of crystal boundary effects in ensembles of finite crystals may be required in future work. The product, , can then be shown to be equivalent to the Fourier transform of the shape function of the crystal centred at the Bragg positions, , given that
where indicates the Fourier transform operation and the presence of whole unit cells has been assumed. Substitution into equation (20) then provides the result,
showing the equivalence of the models for the scattering factor of a nanocrystal for the cases we consider here.
A consequence of this is that it is implicitly assumed that the sampling rate for the scattering factor for the i.e. interpolation between the Bragg locations.
, is the Shannon–Nyquist sampling rate when considering `overlapping' Bragg peak distributions in this paper. The spreading of Bragg peak distributions can also be understood as incorporating the interpolation of to locations in between these sampling points,Supporting information
Supplementary section containing structurefactor moduli extracted via wholepattern fitting and integration. DOI: 10.1107/S2052252516001238/cw5009sup1.pdf
Acknowledgements
The authors acknowledge the support of the Australian Research Council through its Centre of Excellence for Advanced Molecular Imaging.
References
Aquila, A. et al. (2012). Opt. Express, 20, 2706–2716. Web of Science CrossRef CAS PubMed Google Scholar
Barty, A., Kirian, R. A., Maia, F. R. N. C., Hantke, M., Yoon, C. H., White, T. A. & Chapman, H. (2014). J. Appl. Cryst. 47, 1118–1131. Web of Science CrossRef CAS IUCr Journals Google Scholar
Beyerlein, K. R. (2011). Simulation and modeling of the powder diffraction pattern from nanoparticles: studying the influence of surface strain. PhD thesis, Georgia Institute of Technology, USA. Retrieved from https://smartech.gatech.edu/handle/1853/41211 . Google Scholar
Boutet, S. et al. (2012). Science, 337, 362–364. CrossRef CAS PubMed Google Scholar
Caylor, C. L., Dobrianov, I., Lemay, S. G., Kimmer, C., Kriminski, S., Finkelstein, K. D., Zipfel, W., Webb, W. W., Thomas, B. R., Chernov, A. A. & Thorne, R. E. (1999). Proteins, 36, 270–281. CrossRef PubMed CAS Google Scholar
Chabrol, E., Thépaut, M., DezutterDambuyant, C., Vivès, C., Marcoux, J., Kahn, R., ValladeauGuilemond, J., Vachette, P., Durand, D. & Fieschi, F. (2015). Biophys. J. 108, 666–677. CrossRef CAS PubMed Google Scholar
Chapman, H. N. et al. (2011). Nature (London), 470, 73–77. Web of Science CrossRef CAS PubMed Google Scholar
Chen, J. P. J. & Millane, R. P. (2013). J. Opt. Soc. Am. A, 30, 2627–2634. CrossRef Google Scholar
Demirci, H. et al. (2013). Acta Cryst. F69, 1066–1069. Web of Science CrossRef IUCr Journals Google Scholar
Dilanian, R. A., Streltsov, V. A., Quiney, H. M. & Nugent, K. A. (2013). Acta Cryst. A69, 108–118. Web of Science CrossRef CAS IUCr Journals Google Scholar
Elser, V. (2013). Acta Cryst. A69, 559–569. Web of Science CrossRef CAS IUCr Journals Google Scholar
Ewald, P. P. (1940). Proc. Phys. Soc. 52, 167–174. CrossRef CAS Google Scholar
Feher, G. & Kam, Z. (1985). Methods Enzymol. 114, 77–112. CrossRef CAS PubMed Google Scholar
Giacovazzo, C. (2011). Editor. Fundamentals of Crystallography, 3rd ed. Oxford University Press. Google Scholar
Grant, M. L. & Saville, D. A. (1994). J. Phys. Chem. 98, 10358–10367. CrossRef CAS Web of Science Google Scholar
Guinier, A. (1963). Xray Diffraction in Crystals, Imperfect Crystals, and Amorphous Bodies. San Francisco: W. H. Freeman and Company. Google Scholar
Hill, R. J. & Fischer, R. X. (1990). J. Appl. Cryst. 23, 462–468. CrossRef CAS Web of Science IUCr Journals Google Scholar
Hosemann, R. & Bagchi, S. N. (1962). Direct Analysis of Diffraction by Matter. Amsterdam: NorthHolland Publishing Company. Google Scholar
Ino, T. & Minami, N. (1979). Acta Cryst. A35, 163–170. CrossRef CAS IUCr Journals Web of Science Google Scholar
Ino, T. & Minami, N. (1984). Acta Cryst. A40, 538–544. CrossRef CAS Web of Science IUCr Journals Google Scholar
Izumi, F. (2003). J. Ceram. Soc. Jpn, 111, 617–623. Web of Science CrossRef CAS Google Scholar
James, R. W. (1954). The Crystalline State, Vol. II, The Optical Principles of the Diffraction of Xrays. London: G. Bell and Sons. Google Scholar
Johansson, L. C. et al. (2012). Nat. Methods, 9, 263–265. Web of Science CrossRef CAS PubMed Google Scholar
Johansson, L. C. et al. (2013). Nat. Commun. 4, 2911. Web of Science PubMed Google Scholar
Kabsch, W. (2010). Acta Cryst. D66, 125–132. Web of Science CrossRef CAS IUCr Journals Google Scholar
Kabsch, W. (2014). Acta Cryst. D70, 2204–2216. Web of Science CrossRef IUCr Journals Google Scholar
Kirian, R. A., Bean, R. J., Beyerlein, K. R., Barthelmess, M., Yoon, C. H., Wang, F., Capotondi, F., Pedersoli, E., Barty, A. & Chapman, H. N. (2015). Phys. Rev. X, 5, 011015. Google Scholar
Kirian, R. A., Bean, R. J., Beyerlein, K. R., Yefanov, O. M., White, T. A., Barty, A. & Chapman, H. N. (2014). Philos. Trans. R. Soc. Lond. B Biol. Sci. 369, 20130331. CrossRef PubMed Google Scholar
Kirian, R. A., Wang, X., Weierstall, U., Schmidt, K. E., Spence, J. C. H., Hunter, M., Fromme, P., White, T., Chapman, H. N. & Holton, J. (2010). Opt. Express, 18, 5713–5723. Web of Science CrossRef PubMed Google Scholar
Kirian, R. A., White, T. A., Holton, J. M., Chapman, H. N., Fromme, P., Barty, A., Lomb, L., Aquila, A., Maia, F. R. N. C., Martin, A. V., Fromme, R., Wang, X., Hunter, M. S., Schmidt, K. E. & Spence, J. C. H. (2011). Acta Cryst. A67, 131–140. Web of Science CrossRef CAS IUCr Journals Google Scholar
Koopmann, R. et al. (2012). Nat. Methods, 9, 259–262. Web of Science CrossRef CAS PubMed Google Scholar
Kupitz, C. et al. (2014). Nature (London), 513, 261–265. Web of Science CrossRef CAS PubMed Google Scholar
Langford, J. I. (1999). Defect and Microstructure Analysis by Diffraction, edited by R. L. Snyder, J. Fiala & H.J. Bunge, pp. 59–81. Oxford University Press. Google Scholar
Le Bail, A., Duroy, H. & Fourquet, J. L. (1988). Mater. Res. Bull. 23, 447–452. CrossRef CAS Web of Science Google Scholar
Leslie, A. G. W. & Powell, H. R. (2007). Evolving Methods for Macromolecular Crystallography. NATO Science Series, Vol. 245, pp. 45–51. Google Scholar
Liu, W. et al. (2013). Science, 342, 1521–1524. Web of Science CrossRef CAS PubMed Google Scholar
Liu, H., Zatsepin, N. A. & Spence, J. C. H. (2014). IUCrJ, 1, 19–27. Web of Science CrossRef CAS PubMed IUCr Journals Google Scholar
Malkin, A. J. & Thorne, R. E. (2004). Methods, 34, 273–299. Web of Science CrossRef PubMed CAS Google Scholar
Neutze, R., Wouts, R., van der Spoel, D., Weckert, E. & Hajdu, J. (2000). Nature (London), 406, 752–757. Web of Science CrossRef PubMed CAS Google Scholar
Palosz, B., Grzanka, E., Gierlotka, S., Stel'makh, S., Pielaszek, R., Lojkowski, W., Bismayer, U., Neuefeind, J., Weber, H.P. & Palosz, W. (2003). Phase Transitions, 76, 171–185. CrossRef CAS Google Scholar
Patterson, A. L. (1939). Phys. Rev. 56, 972–977. CrossRef CAS Google Scholar
Redecke, L. et al. (2013). Science, 339, 227–230. Web of Science CrossRef CAS PubMed Google Scholar
Rietveld, H. M. (1967). Acta Cryst. 22, 151–152. CrossRef CAS IUCr Journals Web of Science Google Scholar
Spence, J. C. H., Kirian, R. A., Wang, X., Weierstall, U., Schmidt, K. E., White, T., Barty, A., Chapman, H. N., Marchesini, S. & Holton, J. (2011). Opt. Express, 19, 2866–2873. Web of Science CrossRef CAS PubMed Google Scholar
Spence, J. C. H., Weierstall, U. & Chapman, H. M. (2012). Rep. Prog. Phys. 75, 102601. CrossRef PubMed Google Scholar
Suortti, P., Ahtee, M. & Unonius, L. (1979). J. Appl. Cryst. 12, 365–369. CrossRef CAS IUCr Journals Web of Science Google Scholar
Suortti, P. & Jennings, L. D. (1977). Acta Cryst. A33, 1012–1027. CrossRef CAS IUCr Journals Google Scholar
Ungár, T. & Gubicza, J. (2007). Z. Kristallogr. 222, 114–128. Google Scholar
Welberry, T. R. (2004). Diffuse Xray Scattering and Models of Disorder. Oxford University Press. Google Scholar
White, A. (2014). Philos. Trans. R. Soc. Lond. B Biol. Sci. 369, 20130330. CrossRef PubMed Google Scholar
White, T. A., Kirian, R. A., Martin, A. V., Aquila, A., Nass, K., Barty, A. & Chapman, H. N. (2012). J. Appl. Cryst. 45, 335–341. Web of Science CrossRef CAS IUCr Journals Google Scholar
Yefanov, O., Gati, C., Bourenkov, G., Kirian, R. A., White, T. A., Spence, J. C. H., Chapman, H. N. & Barty, A. (2014). Philos. Trans. R. Soc. Lond. B Biol. Sci. 369, 20130333. CrossRef PubMed Google Scholar
Young, R. A. & Wiles, D. B. (1982). J. Appl. Cryst. 15, 430–438. CrossRef CAS Web of Science IUCr Journals Google Scholar
This is an openaccess article distributed under the terms of the Creative Commons Attribution (CCBY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.