Abstract
Coherent diffraction imaging (CDI) allows the retrieval of an isolated object’s structure, such as a macromolecule, from its diffraction pattern. CDI requires the fulfillment of two conditions: the imaging radiation must be coherent and the object must be isolated. We discuss that it is possible to directly retrieve the molecular structure from its diffraction pattern, which was acquired neither with coherent radiation nor from an individual molecule. This is provided that the molecule exhibits periodicity in one direction, as in the case of fiber diffraction. We demonstrate that, when we apply iterative phase retrieval methods to a fiber diffraction pattern, the repeating unit; that is, the molecule structure, can directly be reconstructed without any prior modeling. For example, we recover the the DNA double helix’s structure in three-dimensions from its two-dimensional X-ray fiber diffraction pattern, Photograph (Photo) 51, which was acquired in Raymond Gosling and Rosalind Franklin’s famous experiment at a resolution of 3.4 Å.
© 2018 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. Introduction
In 1953 James Watson and Francis Crick deduced the structure of B-DNA from its X-ray fiber diffraction pattern which was acquired by Raymond Gosling, a PhD student of Rosalind Franklin [1–3]. Fiber diffraction [4–6] is employed for studying the molecular structure of elongated assemblies of repeating units. In the case of a subunit consisting of a few atoms, such as the DNA nucleotide, the atomic positions can be determined by combining X-ray diffraction data with the known chemical composition of the unit. In the case of larger subunits, such as proteins, the fiber diffraction pattern can provide the information about the shape of the subunit. However, to determine its detailed structure X-ray crystallography data are required. Coherent diffractive imaging (CDI) [7] is a relatively modern technique which unlike other techniques allows to directly visualize the structure of an isolated object, a macromolecule for example [7–17]. In CDI, the structure of a sample is reconstructed from its diffraction pattern by applying numerical iterative phase retrieval algorithms. For realizing CDI, two requirements must to be fulfilled: the object under study must be isolated and the imaging radiation must be coherent, although some attempts to employ partially coherent waves have also been reported [18]. The second requirement of a high brilliance coherent X-ray source is fulfilled by using a synchrotron or one of the recently build X-ray free-electron lasers. The very short X-ray pulses of the latter are to be used for imaging individual biological macromolecules [19].
In this study we show that a fiber diffraction pattern can be treated as if it was acquired in a CDI experiment. As an example, we reconstruct the three-dimensional (3D) DNA double helical structure from its two-dimensional (2D) fiber diffraction pattern by applying a phase retrieval algorithm. In this way, no modeling is thus required. The oversampling condition normally needed in the case of CDI can be relaxed when imaging a periodical structure like a DNA. Moreover, the molecule under study does not need to be isolated; the molecular structure can be directly recovered from its fiber diffraction pattern, a diffraction pattern created by an ensemble of molecules acquired with a partially coherent radiation.
2. Fiber diffraction and coherent diffraction imaging (CDI)
A fiber diffraction experiment, which typically employs a partially coherent X-ray source, is sketched in Fig. 1(a). A coherence length which ensures coherent addition of the waves scattered by adjacent atoms to create a peak in the diffraction pattern is sufficient. The molecules in a fiber are approximately aligned along the fiber axis, and, unlike in a crystal, they are randomly rotated along their axes. Fibers show helical symmetry [4–6]. The molecular structure is recovered in a similar to crystallographic experiment approach, starting with an initial model which is then refined to achieve a good agreement between atomic positions and intensity maxima in the fiber diffraction pattern.
A CDI experiment is sketched is in Fig. 1(b). Here, a coherent wave is employed, and the scattered wavefront in the far-field distribution is given by the Fourier transform of the distribution of atoms within the molecule. If the complex-valued distribution of the scattered wavefront could be detected, a simple inverse Fourier transform would give the molecular structure. However, since only the intensity can be recorded, the phase information is missing, but can be recovered provided the intensity distribution is sampled at a frequency finer than at least twice the Nyquist frequency (oversampled) [7]. A finer sampling of the intensity in the far-field (Fourier) domain is equivalent to zero-padding of the sample in the object domain. In reality it means that the sample must be isolated and surrounded by a non-scattering medium. The size of the sample given by S0 × S0 should not exceed in each direction, where is the extent of the reconstructed field of view and is the pixel size in the Fourier domain. Provided this oversampling condition is fulfilled, the number of equations (that is equal to pixels of the sampled diffraction pattern) exceeds the number of unknowns (that is equal to for a complex-valued object) and the sample structure can then uniquely be recovered from its diffraction pattern [22]. This sample structure recovery is conventionally done by applying numerical iterative methods [23,24]. Examples of various objects, like bacteria or viruses, which were recovered by CDI from their diffraction patterns can be found in the literature [7–17].
3. Fiber diffraction imaging (FDI)
The method considered in the present study is the reconstruction of the repeating unit structure obtained by iterative phase retrieval of a fiber diffraction pattern. In order to distinguish the approach considered in this work from conventional fiber diffraction and CDI we name it as fiber diffraction imaging (FDI).
Given the density of an object its Fourier transform is given by
where is the coordinate in the object domain and is the coordinate in the Fourier domain. For digitally sampled signals Eq. (1) can be re-written as:where and are the pixel sizes in the object, respectively Fourier domain. When a fast Fourier transform (FFT) is applied to calculate the Fourier transform given by Eq. (2):and the pixel size in the object and Fourier domains are related by the following formula:The measured intensity distribution is given by:which constitutes a set of equations.Now, we consider the formation of a diffraction pattern of a sample which is an infinite one-dimensional (1D) chain of repeating units. The distribution of such a sample can be written as:
where is the transmission function of the repeating unit, is the Shah (Dirac comb) function describing the periodical positions of the repeating units. The intensity of the diffracted wave is given by:where is the Fourier transform of the repeating unit, FT denotes the Fourier transform, and is the Shah-function describing the positions of the diffraction peaks. According to Eq. (7) the intensity distribution of the resulting diffraction pattern is given by a product of the diffraction pattern of the repeating unit and the diffraction peaks distribution.For a perfect infinite crystal, the signal between the diffraction peaks is zero and the intensity of the diffraction peaks is modulated by the diffraction pattern of a repeating unit [25]. It has recently been shown that when the crystal is not perfect (as a typical realistic crystal is), then the signal between the diffraction peaks is non-zero and creates a so-called “incoherent” background which allows mapping of and consequent reconstruction of the repeating unit [26]. In the present study we do not consider an imperfect crystal and our method is thus not based on using the “incoherent” background between the diffraction peaks. We assume that the crystal is close to being perfect and that the intensity between the diffraction peaks is almost zero.
The principle of our approach is illustrated in Fig. 2. Figure 2(a) shows the repeating unit distribution - a cartoon of a “cat”. The distribution of the “cat” occupies almost the entire region of the repeating unit cell and thus it is not oversampled within the repeating unit. The corresponding simulated pattern is shown in Fig. 2(b). Here, the diffraction pattern was simulated by solving the analytical integral for far-field diffraction as given by Eq. (5), without applying fast Fourier transform (FFT) to avoid wrapping effects [27]. The diffraction pattern shown in Fig. 2(b) is not oversampled, it is sampled exactly at the Nyquist frequency, so that the inverse Fourier transform of the complex-valued distribution gives the object distribution, as shown in Fig. 2(c).
Next, a diffraction pattern was simulated by solving the analytical integral for far-field diffraction as given by Eq. (5), where the pixel size is four times smaller and the number of pixels N is four times larger than in the previously simulated diffraction pattern, shown in Fig. 2(d). The simulated diffraction pattern is sampled at 4 times the Nyquist frequency, and thus oversampled, with an oversampling ratio σ = 4. The inverse Fourier transform of the simulated complex-valued distribution delivers the object distribution which is zero-padded, as shown in Fig. 2(e).
Next, we create the diffraction pattern of an infinite 1D chain of repeating units, as described by Eq. (7), by multiplying the oversampled diffraction pattern of the repeating unit with a function The function is given by the Fourier transform of the 1D Dirac comb-function and in two-dimensions it is a set of equidistant lines as shown in Fig. 2(f). The period of is determined as follows. The size of the repeating unit is given by and the size of the field of view is given by with The oversampling ratio is given by where is the size of the repeating unit in pixels. The peaks due to the structure of period are expected at k-values of Or, in pixels: , that is at each σth pixel. In the simulated example considered here σ = 4 and the non-zero horizontal lines appear at every 4th line, as shown in Fig. 2(f). The product of and that is is shown in Fig. 2(g). Note that this diffraction pattern is not simply the FFT of 4 repeating units, this diffraction pattern is the diffraction pattern of an infinite 1D chain of repeating units. The inverse Fourier transform of results in the repeating units arranged into vertical line, as illustrated in Fig. 2(h).
Now, we discuss the possibility of reconstructing the repeating unit distribution from a diffraction pattern of a 1D chain of the repeating units. The intensity distribution of such diffraction pattern is given by Eq. (7). For a 2D diffraction pattern the number of equations is N2. For a 2D real-valued signal, according to Friedel's law [28], the measured diffraction pattern exhibits central symmetry. Thus, the total number of equation is N2/2. In the corresponding diffraction pattern, only every σth line has non-zero values which results in a total number of lines with non-zero values: Therefore, the total number of equations for a such a sample is given by The number of unknowns corresponds to the number of pixels of the repeating unit: In order for the system of equations to have a solution, the number of equations must exceed the number of unknowns which translates into: or When the oversampling condition is fulfilled the repeating unit structure can be recovered from the corresponding diffraction pattern. A typical fiber diffraction pattern has a four-quadrant symmetry, and thus the number of equations is reduced to which leads to the oversampling condition We would like to emphasize that only real-valued objects (repeating unit) are considered here because it follows from the four-quadrant symmetry of a typical fiber diffraction pattern.
Similar considerations can be applied for a diffraction pattern of a 2D crystal. In this case, the number of equations corresponds to the number of measured peaks and equals to for a crystal with repeating unit of identical size in both dimensions. Since the number of unknowns is the system of equations Eq. (7), that is, the structure of the repeating unit, can thus not be solved.
The 2D diffraction pattern shown in Fig. 2(g) was reconstructed by applying iterative phase retrieval procedure employing the hybrid input-output (HIO) algorithm with a constant cylindrical object support mask (more details are provided in Appendix A). The resulting 2D reconstruction is shown in Fig. 2(i). It exhibits the repeating units arranged into a vertical line, where the distribution of the repeating unit (“cat”) is correctly reconstructed. Phase retrieval transfer function (PRTF) [9,29] was calculated by averaging 10 reconstructions with the least error (details on calculation of error are provided in Appendix A), shown in Fig. 2(j); PRTF exhibits relatively high values, decreasing to about 0.7 at the highest frequency.
Since the oversampling ratio is sufficient to recover the repeating unit, and thus the problem can have a solution, there are various ways of finding the solution. For example, the phase retrieval algorithm can be modified to reconstruct only one single repeating unit as follows: (1) In the object plane the object support mask is set to the area of (2) In the Fourier domain, the values of the pixels where the measured intensity was zero are updated with the values obtained after each iteration. This way, during the reconstruction, the zero-intensity values in the Fourier domain are recovered to non-zero values and the diffraction pattern turns into a diffraction pattern of the repeating unit, as shown in Fig. 2(k). A single reconstruction with the least error is shown in Fig. 2(l) and the repeating unit distribution obtained by averaging 10 reconstructions with the least error is shown in Fig. 2(m) exhibiting the correct recovered distribution of the “cat”. PRTF was calculated by averaging 10 reconstructions with the least error (details on calculation of error are provided in Appendix A), shown in Fig. 2(n); PRTF exhibits relatively high values, decreasing to about 0.8 at the highest frequency.
An example of recovery non-binary repeating unit from the corresponding diffraction pattern is shown in Fig. 3. This example demonstrates that the method in principle should allow recovery of scatterers (such as chemical elements, for example) with different scattering amplitudes. Some artifact noise can be observed in the reconstructions shown in Fig. 2(m) and 3(d) which can be greatly reduced if a tight support mask is applied in the object domain. This can be done for example, by applying the shrinkwrap algorithm [24] which automatically adjusts the support mask tightly to the object shape.
In all reconstructions obtained from the simulated examples and shown here, no multiple solutions were observed, the reconstructions obtained with the lowest errors all looked very similar. We can explain the uniqueness of the solution by the fact (already mentioned above) that diffraction pattern of a sample that is a 1D chain of the repeating units can be represented as a product of the oversampled diffraction pattern of the repeating unit and a function This means that the resulting diffraction pattern can be considered as the oversampled diffraction pattern of the repeating unit where some of measured amplitudes are missing. It has been previously demonstrated that correct sample distribution can be reconstructed from its sufficiently oversampled diffraction pattern even if only limited amount of the measured amplitudes is available (even if only 10% of the measure amplitudes are available) [30]. This situation is similar to the example shown in Figs. 2(k) – 2(m), where the zero-intensity values in the Fourier domain are recovered to non-zero values and simultaneously the repeating unit is reconstructed during the iterative procedure.
Previously, Zuo et al reported the reconstruction of a carbon nanotube where the sample was not constrained along the direction of the nanotube, which implies that the diffraction pattern was not oversampled in that direction [31]. Although Zuo et al did not discuss the uniqueness of their reconstruction, their result confirms our hypothesis that the repeating unit distribution can be reconstructed from a diffraction pattern of a 1D chain of the repeating units, provided the diffraction pattern is sufficiently oversampled to recover the repeating unit.
4. Simulations
Before dealing with the experimental data, we show that it is possible to reconstruct a fragment of a 3D helical molecule structure from a 3D diffraction pattern of an infinite helical molecular structure. Thereafter, we demonstrate that the 3D helical molecule structure can be reconstructed from a 2D slice of its 3D angular-averaged diffraction pattern.
4.1 Diffraction pattern of a single helix
We considered a sample in the form of a 3D helical curve, the helix axis is oriented along the -axis. The far-field distribution of the scattered intensity is provided by the Fourier transform of the distribution of the scatterers. For a helical arrangement it has an analytical solution expressed through Bessel functions [32]:
where , are the n-th order Bessel functions of the first kind, is the radius of the helix, , is the imaginary number, and are the coordinates in the Fourier domain. Amplitude of and thus the intensity in the far-fieldhave non-zero values only at where is the period of the helix [32].For a discontinuous helix the total intensity is given by
Here, we accounted for the zero-th (the term ) and the first order (the terms and ) of diffraction on a discontinuous lattice, being the period of the lattice (the -distance between the scatterers).Figure 4 shows the diffraction patterns of the repeating unit and the diffraction patterns of an infinite 1D chain of the repeating units, for both continuous and discontinuous single helices. The parameters of the simulations shown in Figs. 4 – 6 were chosen to match those of an X-ray diffraction experiment on B-DNA [20]: Å is the distance between the two helices along the -axis, Å is the period of the helical turn, Å is the distance between the scatterers along the -axis with 10 scatterers per helical turn. The distributions were sampled with 128 × 128 × 128 pixels.
From Figs. 4(a) – 4(e), it is evident that the repeating unit in form of one continuous helical turn provides the “X”-shaped distribution of the intensity in the diffraction pattern. When repeating units are arranged periodically, the diffraction lines (analogous to diffraction peaks for crystals) appear in the diffraction pattern, see Figs. 4(f) – 4(j). It should be noted that 3D diffraction pattern of single continuous helix according to Eq. (9) exhibits no angular dependency, which is also illustrated in the simulations shown in Fig. 4(f) – 4(j). Therefore, a single 2D diffraction pattern is sufficient for reproducing the entire 3D diffraction pattern by a simple rotation.
A discontinuous helix (made up of individual scatterers), with 10 scatterers per helical turn, produces a diffraction pattern which can be interpreted as a superposition of continuous helix diffraction patterns shifted along the -axis, see Figs. 4(k) – 4(t). Here, the 3D diffraction pattern of single discontinuous helix according to Eq. (10) exhibits angular dependency, which is also illustrated in the simulations shown in Fig. 4(k) – 4(t).
4.2 Diffraction pattern of a discontinuous double helix
For a discontinuous double helix the intensity distribution in the diffraction pattern is given by:
where the modulating term accounts for the second helix shifted by relatively to the first helix.The 3D diffraction pattern of a discontinuous double helix simulated from the analytical solution given by Eq. (11) is shown in Fig. 5(a). The diffraction pattern exhibits a 10-fold angular dependency in the -plane (shown in Fig. 5(b)), as provided by the 10-fold symmetry of the arrangement of the scatterers in the -plane. The 2D intensity distributions in the , and -planes are shown in Figs. 5(b) – 5(d), correspondingly. Since the simulated 3D object considered here is a real-valued object, according to Friedel's law, the intensity of its diffraction pattern is centrosymmetric [28]. It should be noted that the diffraction pattern simulated according to Eq. (11) corresponds to a helix not finite (but infinite) in z-direction.
A 3D object reconstructed from its 3D diffraction pattern by applying an iterative phase retrieval algorithm (the reconstruction algorithm is described in Appendix A) is presented in Fig. 5(e), showing scatterers arranged into a double helix.
Projected distributions were calculated from the reconstructed 3D distribution as follows:
The projected distribution exhibits ten scatterers arranged into a circle, as shown in Fig. 5(f). The projected distributions and exhibit a discontinuous double helix, as shown in Figs. 5(g) and 5(h), respectively.4.3 Angular-averaged diffraction pattern
As demonstrated in the previous subsection, a 3D diffraction pattern of a discontinuous double helix exhibits an angular dependency. However, in a fiber diffraction experiment, the individual molecules in the fiber are randomly rotated around their axes and the resulting diffraction pattern is thus an angular-averaged diffraction pattern. The intensity distribution of an angular-averaged 3D diffraction pattern is obtained by averaging the intensity distribution given by Eq. (11) over angle :
where , and are integer numbers which satisfy the following equations:, , and . An angular-averaged 3D diffraction pattern, simulated according to Eq. (13), and its slices are shown in Figs. 6(a) – 6(d). In Eq. (13), the terms which provided the angular ten-fold symmetry observed in the -planes and the atomic arrangement information associated with them are gone (compare Figs. 5(b) and 6(b)). However, the terms related to the first-order diffraction ( and ) which are also associated with the angular atomic arrangement remain. Thus, an angular-averaged 3D diffraction pattern of a helical structure described by Eq. (13) still contains information about angular arrangement of the scatterers.A 3D object reconstructed from its angular-averaged diffraction pattern (the reconstruction algorithm is described in Appendix A) is shown in Fig. 6(e). Surprisingly, the projected distributions exhibit again a discontinuous double helix (Figs. 6(f) – 6(h)) just as the reconstruction obtained from non-angular-averaged diffraction pattern (Figs. 5(f) – 5(h)). Here, during the reconstruction the following observation was made: while the overall double helical shape was stable and the distances between scatterers remained constant, the position of the scatterers were “oscillating” around their averaged positions.
The successful reconstruction of double helix from its angular-averaged diffraction pattern can be explained as follows. The angular-averaged 3D diffraction pattern (Fig. 6(a) – 6(d), described by Eq. (13) does not exhibit angular dependency. This is the same situation as for not angular-averaged diffraction pattern of single continuous helix (Fig. 4(f) – 4(j) described by Eq. (9)), which also does not exhibit angular dependency but 3D helical structure can be successfully recovered. The following consideration can be applied. According to the projection-slice theorem, the projected object can be reconstructed from the 2D slice through the origin of the 3D reciprocal space (for more details, see Appendix B). By rotating the slice, the angle of the projection of the object is accordingly rotated. For a helical structure, a rotation of helix around its axis is equivalent to shifting the helix along that axis. Since a lateral shift of the object only causes an additional linear phase in its Fourier spectrum, the intensity of the diffraction pattern remains unchanged when the object is laterally shifted. Thus, when the helical structure is rotated around its axis, the corresponding 2D diffraction pattern remains almost unchanged (compare Figs. 4(i) and 4(j) which are diffraction patterns of discontinuous double helix rotated by 90°).
The successful reconstruction of individual scatterers arranged into a double helix from its angular-averaged diffraction pattern can be explained as follows. An angular-averaged 3D diffraction pattern of a helical structure (as described by Eq. (13)) still contains information about angular arrangement of the scatterers, which can be in principle extracted in the reconstruction. When the helical structure is rotated around its axis, the corresponding 2D diffraction pattern remains almost unchanged (compare Figs. 4(i) and 4(j) for single continuous helix and Fig. 5(c) and 5(d) for discontinuous double helix). The only changes are associated with the precise positions of the scatterers. This uncertainty is within the distance between the scatterers.
A single 2D slice of the angular-averaged diffraction pattern is therefore sufficient for reproducing the entire angular averaged 3D diffraction pattern by a simple rotation. Thus, a single 2D diffraction pattern obtained from a fiber should be sufficient for the reconstruction of the 3D structure of the molecule (repeating unit).
4.4 Reconstruction of other helical structures
Figure 7 shows further examples of helical structures reconstructed from their 2D angular-averaged diffraction patterns.
5. Reconstruction of the experimental X-ray DNA diffraction pattern, Photo 51
To proof the concept, we selected the fiber diffraction pattern of DNA acquired by Gosling and Franklin, which was called by the BBC as “the most important photo ever taken“. The well-known dramatic story behind the Photograph 51 demonstrates how important the knowledge of the DNA molecule structure, not just its chemical composition, was. In this section we reconstruct the 3D DNA double helical shape from its 2D fiber diffraction pattern without using any prior information about the structure or any modeling.
The experimental setup, as it was employed in the X-ray fiber diffraction experiment by Gosling and Franklin, is sketched in Fig. 1(a). A fibre of calf thymus DNA, 50 μm in diameter, was placed at a distance of 15 mm in front of a photographic film which was mounted onto a cylindrical holder and kept under continuous X-ray radiation of wavelength Å for about 62 hours [20]. Parameters of Photo 51 are provided in Appendix C. A high-resolution image of the DNA diffraction pattern Photo 51 was downloaded as a scan of Raymond Gosling PhD thesis [20], it was optimized (as described in Appendix D), and is shown in Fig. 8(a). Interpretation of Photo 51 is provided in Appendix E. While the double helical structure of DNA and the periodicity of the bases can be predicted from the positions of the peaks (layer lines) in the diffraction pattern, applying a Fourier transform to the diffraction pattern does not reveal a helical structure (demonstrated in Appendix F).
By comparing the simulated and the experimental diffraction patterns, shown in Figs. 6(c) and 8(a), respectively, one immediately notices that the experimental diffraction pattern appears to be rotationally blurred. The rotational blurring can be caused by DNA molecules in the fiber being slightly rotated relatively to the axis of the fiber, and/or by bending segments of DNA where the bending can reach 20° when molecules are packed into a crystal [21], both situations are illustrated in Fig. 8(b). A good match was achieved between the experimental diffraction pattern and a sum of the simulated diffraction patterns when each simulated diffraction pattern was randomly rotated around its center with the rotations normally distributed with the standard deviation of 9° (more details are provided in Appendix G). This led to the conclusion that the DNA molecules in the fiber or their segments were randomly rotated around their axes with the angle of rotations normally distributed with the standard deviation of 9°. Next, the experimental diffraction pattern was rotationally deblurred and the result is shown in Fig. 8(c) (the deblurring procedure is described in Appendix H). It should be noted that the deblurring procedure is optional, but without deblurring only a finite fraction of the sample can be reconstructed (an example is provided in Appendix I).
Figure 8(d) shows the 3D diffraction pattern which was subject to the reconstruction. The low-resolution reconstructions, obtained from the 3D diffraction pattern where the signal corresponding to higher than 5-th order peaks (layer lines) was set to zero, clearly show a double helix, Figs. 8(e) and 8(g). The high-resolution reconstructions, obtained from the entire diffraction pattern, are shown in Figs. 8(h) - 8(j). The projected distributions calculated from the reconstructed 3D distribution, depicted in Figs. 8(f) and 8(i) show that the reconstructed double helix exhibits not a hollow but a dense center, which agrees well with the molecular B-DNA model [33]. High-resolution reconstructions, shown in Fig. 8(h) and 8(j), demonstrate that the helices consist not from point-like scatterers aligned along the double helix curves, but of stick-like features connecting the opposite sides of the helix. These features can be attributed to the reconstructed bases. All of the atoms that constitute a base contribute to the X-ray scattering diffraction pattern with the maxima of the scattered amplitudes for 8 keV X-rays: 1.00 (H), 6.00 (C), 7.03 (N), 8.05 (O), 15.31 (P) [34]. It is very fortunate that atoms such as P and O, that scatter the strongest due to their relatively high electron density, are positioned almost exactly on the curves that form the double helix. A closer look at the high-resolution reconstruction shows that the individual bases are not orthogonal to the helix axis but rather slightly tilted, Fig. 8(j).
6. Summary and conclusions
We have shown that the reconstruction of the repeating unit can be obtained from the diffraction pattern of a 1D chain of the repeating units by conventional phase retrieval methods, provided the diffraction pattern is sufficiently oversampled to recover the repeating unit. PRTFs calculated for reconstructions obtained from simulated diffraction patterns exhibit relatively high values reaching 0.7 – 0.8 at the highest frequency, thus indicating high reproducibility of the reconstructions obtained by our method.
We demonstrated our method by applying it to fiber diffraction data where typically long assemblies of repeating units are imaged. Due to the partial coherence of the sources typically applied in fiber diffraction experiments, the interference of wavefronts scattered by different molecules in the fiber can be neglected and thus the resulting diffraction pattern can consequently well be approximated by a diffraction pattern of an individual molecule (repeating unit). A fiber diffraction pattern can thus be treated as if it was recorded in a CDI experiment. For a periodical object, the CDI requirement that the sample must be isolated (zero-padded in a mathematical sense) is replaced by the much more relaxed condition that the size of the reconstructed area should exceed at least twice the size of the repeating unit. When applying our method to an angular-averaged diffraction pattern given by a fiber diffraction experiment, the 3D structure of the molecule can be reconstructed from the 2D diffraction pattern without any need for modeling or prior information about the molecule as it is required by conventional reconstruction procedure in fiber diffraction methods [4–6]. The retrieved structure is not a uniquely-defined structure of a particular isolated molecule, but is a molecular structure averaged over an ensemble, that is - structure of the repeating unit. The achievable resolution by this approach is defined by the highest resolvable peak in the diffraction pattern. For example, for the X-ray DNA diffraction pattern Photo 51, the individual bases are resolved in the reconstruction which corresponds to the resolution of 3.4 Å. While our approach recovers the geometrical arrangement of the scatterers, the contrast of a reconstructed scatterer is proportional to its scattering amplitude and thus the chemical nature of individual scatterers composing the molecule can in principle be derived.
Fiber diffraction allows imaging of some of the molecules that cannot be imaged by X-ray crystallography or CDI, where our technique can simplify the structure determination. The vast amount of already acquired fiber diffraction patterns can in principle be post-processed by using our technique and molecular structures can directly be viewed, as we have demonstrated with the historical DNA diffraction pattern.
Recently, Seuring et al demonstrated an application of the X-ray free electron laser radiation for imaging individual tobacco mosaic viruses (a helical structure) deposited onto graphene [35] where the resulting diffraction patterns are analogous to that obtained in a fiber diffraction experiment [36]. This opens up a new application for the fiber diffraction technique where the reconstruction method proposed here can serve as a simple tool for the reconstruction of the 3D molecular structure from a 2D experimental fiber diffraction pattern.
Appendix A Phase retrieval algorithms
- (a) Hybrid input-output algorithm
The employed hybrid input–output (HIO) algorithm [23] with the feedback parameter includes the following steps. In the first iteration, random phases are combined with the measured amplitudes. The real part of the inverse Fourier transform of the resulting complex-valued distribution provides the object distribution ( for the first iteration). Then, the iterative procedure is applied, which includes the following steps:
- (i) A Fourier transform of provides the complex-valued distribution , where are the coordinates in the detector plane.
- (ii) The iterated amplitudes are replaced with the measured amplitudes and thus an updated distribution is formed: . At the points, where the measured intensities are missing (central part of diffraction pattern), the iterated amplitudes are adapted for the next iteration: .
- (iii) The real part of the inverse Fourier transform of provides .
where denotes the set of the points where satisfies the support and the reality and positivity constraints, so that. The resultant is plugged into the next iteration starting at (i).
The object support mask is selected in form of a cylindrical mask with a fixed radius of 9 pixels.
The quality of the reconstruction was evaluated by calculating the error function in the object plane as it is more appropriate for the HIO algorithm [22,23]:
For the 3D reconstructions, 1 out of 10 reconstructions with the least error was selected. Although one reconstruction with the least error is in principle sufficient, averaging over several reconstructions with the least error provides a final reconstruction with superior contrast. While we did such an averaging in the 2D case, it was however not feasible for the 3D case due to limited computation power.
- (b) Reconstruction of the experimental diffraction pattern
The reconstruction of the experimental diffraction pattern was done by applying the error-reduction [23] algorithm. The missing intensity values due to the beam stop were replaced with the iterated values after each iteration. The radius of the cylindrical object support mask in the object was varied and the reconstructions with the least error were obtained when the radius of the cylindrical object support mask was 9 pixels (corresponding to approximately 12 Å). The quality of the retrieved complex-valued amplitudes in the detector plane was evaluated by calculating the mismatch between the measured and the iterated amplitudes, or the error [23]:
where are the experimentally measured amplitudes at the detector, are the iterated amplitudes, i and j are the pixel numbers , and the missing pixels were excluded from the summation. A total of 200 iterations were done for each reconstruction, 10 reconstructions were obtained and the reconstruction with the least error was selected.
Appendix B Reconstruction of slices of diffraction pattern
The intensity of a 3D diffraction pattern of an object is given by:
where and are the coordinates in the real and the Fourier space, respectively. According to the projection-slice theorem, each slice of a 3D diffraction pattern through its origin is a 2D diffraction pattern of the projected object. For example, the 2D intensity distribution at is given bywhere is the projection defined by Eq. (12).2D reconstructions obtained from the slices of the 3D diffraction pattern in the -plane are shown in Fig. 9. The parameters of the simulation of the 3D diffraction pattern are the same mentioned in the main text: Å is the distance between the two helices along the -axis, Å is the period of the helical turn, Å is the distance between the scatterers along the -axis, 10 scatterers per helical turn. The simulated 3D diffraction pattern was sampled with 128 × 128 × 128 pixels. The shown reconstructions were obtained from the slices of the original and the angular-averaged 3D diffraction patterns. The reconstructions were obtained by applying the HIO algorithm [23] with the feedback parameter and a constant cylindrical object support mask. The final reconstruction was obtained by alignment [37] and averaging of 10 out of 100 reconstructions with the least error. During the alignment procedure the individual reconstructions were cyclically shifted so that the replacing values are not zeros but the values of the object distribution obtained by periodical shift of the object distribution.
In both cases, the reconstructions exhibit a projected double helix where individual scatterers can be resolved. In the case of the reconstruction from the slice of the 3D diffraction pattern, every reconstruction converges to an exact solution of the projected double helix, as shown in Fig. 9(b). In the case of the reconstruction from the slice of the angular-averaged 3D diffraction pattern the projected double helix is slightly distorted. It appears as if it is viewed under a slightly different angle, which is caused by the averaging over all the angular information. Such distortions were observed for some, but not for all reconstructions. However, such distortions were not observed when reconstructing the 3D structure from the entire 3D angular-averaged diffraction pattern. Although slightly distorted, the reconstruction in Fig. 9(d) shows a double helix consisting of individual scatterers.
Appendix C Parameters of Photo 51
- (a) DNA sample
The sample was a fiber of calf thymus DNA, 50 μm in diameter, which gives about 6.25 × 106 DNA molecules for a single DNA molecule of 20 Å in diameter.
- (b) Experiment
In this historical experiment, the photographic film was mounted onto a cylindrical holder with a hole in the center to allow the direct beam to go through, which also helped to the optical alignment of the sample orientation in the beam by looking through the hole with a microscope. The distance between the sample and the photographic film was 15 mm [20] (page 66). The X-ray radiation was provided by an Ehrenberg-Spear Fine Focus tube [20] (page 52) with a copper cathode, which resulted in a wavelength of the obtained X-rays of about Å [20] (page 82).
- (c) Image of Photo 51
The scanned page with the Photo 51 was sampled at 2475 × 3182 pixel with 96 dpi resolution and the effective size of the Photo 51 was about 1200 × 1200 pixel.
- (d) Real size of the diffraction pattern
When a periodical structure with a period is illuminated with a coherent radiation of wavelength , a peak at the diffraction angle is observed in the diffraction pattern where satisfies the following equation
For Å and wavelength Å [20], which gives and therefore
With a distance between the sample and the photographic film of 15 mm this peak (layer line) is found at a distance from the center of the diffraction pattern. Thus, the real-world size of the X-ray fiber diffraction pattern between the top and the bottom intense peaks (layer lines) is about 16 mm × 16 mm.
Appendix D Optimization of Photo 51
The diffraction pattern obtained as a scan of Photo 51, was sampled with 1200 × 1200 pixels. Then it was optimized in the following way:
- 1. The image was centered based on the positions of the second-order diffraction peaks (layer lines), thus giving .
- 2. In general, an ideal fiber diffraction pattern exhibits a four quadrant symmetry. A real experimental fiber diffraction pattern can exhibit mirror-symmetry only with respect to the meridian due non-perfect perpendicular orientation of the fiber and the incident beam. The corresponding geometrical distortion can be corrected [38]. In the Photograph 51, the positions of the peaks (layer lines) relatively to the center were found almost symmetrical, so that no corrections for fiber tilt needed to be implemented.
- 2. Individual photographic film defects were numerically filtered out by setting their value to the average value of the surrounding pixels, thus giving .
- 3. The contrast of the image was inverted as .
- 4. A constant background was subtracted: .
where .
- 6. Photo 51 was acquired on photographic film placed over a cylindrical surface, and therefore it is given in the coordinates , where is the length of the arc along the cylindrical surface, as illustrated in Fig. 10.
For each coordinate, there is a corresponding coordinate on a cylindrical surface, which satisfies
where is the radius of the cylinder. The coordinates are given by
where we used the fact that the measured intensities are on cylindrical surface so that and where Note that which corresponds to the intensities measured on the Ewald sphere. The image was transformed to the -coordinates as follows (similar to the procedure described in Ref [39].).
- (i) 2D array of -coordinates was created.
- (iii) The intensity at this coordinate was extracted and assigned to the intensity at the coordinate.
The resulting image in -coordinates is
- 7. The image was zero-padded to 1290 × 1290 pixels, thus giving .
- 8. The image was down-sampled to 129 × 129 pixels by binning with the bin size of 10, thus giving . The odd number of pixels ensures that the zero-order diffraction is represented by a single pixel position in the center of diffraction pattern.
- 9. The image was cropped to 128 × 128 pixels, thus giving .
- 10. After the binning and cropping, the image was double-checked to be centrosymmetric by applying the following transformation:
where , to ensure that the zero-order diffraction intensity is given at , and the intensities at pixels (for example, at pixels and ) have identical values.
Appendix E Interpretation of Photo 51
Figure 11 illustrates interpretation of Photo 51 B-form DNA diffraction pattern. The range of -vectors, where and is the diffraction angle, reaches approximately Å−1. “X”-form distribution of the diffraction peaks (layer lines) is an indication of a helical structure, as shown in Fig. 11(a). As indicated in Fig. 11(b), the intensity exhibits non-zero values at , where is the period of the helical turn and is the order of diffraction . The distance between the lines Å−1 corresponds to the period of the helical turn, Å. Figure 11(c) demonstrates that the broad extended peaks (layer lines) on the top and the bottom at Å−1 are formed by diffraction on small periodical features - the base pairs. The distance between the base pairs is therefore Å.
The missing diffraction spots, as shown in Fig. 11(d), indicate that there is a modulation of the interference pattern. The modulation is caused by interference between the waves diffracted at two helical structures that are shifted relatively to each other along their long axis by – a double helix. The diffraction pattern of a double helix can be represented as a diffraction pattern of a single helix whose intensity is modulated by function, which turns into zero at Thus, . In Photo 51, the missing order of diffraction is which gives Å. The shift between the two helices is approximately Å.
The position of the maxima of the diffraction spots is related to the radius of the DNA helix, as indicated in Fig. 11(e). The intensity distribution in the individual spots in lines 0, 1, 2... are described by the Bessel functions of the first kind: where is the radius of the DNA. From the positions of the maxima of the diffraction spots, the radius can be recovered. For example, for the second line diffraction spots, as indicated in the figure, the intensity distribution is described by which has the first maximum at . The position of the maximum evaluated from the diffraction pattern is about Å, which gives Å, which agrees well with the DNA radius obtained by Gosling [20].
Appendix F Fourier transform of Photo 51 and of simulated 2D diffraction patterns
The Fourier transform of a diffraction pattern of an object delivers the autocorrelation of the object distribution function. Figure 12 shows the amplitudes of the Fourier transform of the optimized Photo 51 and of the simulated diffraction pattern of a discontinuous double helix. From Fig. 12 it is evident that only a limited amount of structural information can be obtained from the Fourier spectra.
For the simulated discontinuous double helix, the amplitude of the Fourier transform of its diffraction pattern, see Figs. 12(a) and 12(b), exhibit a structure which resembles two out-of phase cosines with twenty fringes per cosine period. A double helical structure cannot be established here.
The situation becomes even worse for the diffraction pattern which mimics the experimental diffraction pattern by implementing rotations of the molecules in the fiber, shown in Figs. 12(c) and 12(d). Here, the amplitude of the Fourier transform does not show any features of a discontinuous double helix, not even a periodical cosine-like structure. Similar situation is observed for Photo 51, as shown in Figs. 12(e) and 12(f).
Appendix G Effect of rotation of individual DNA molecules within fiber
The effect of rotated individual DNA molecules in the fiber is shown in Fig. 13. The diffraction pattern of a DNA fiber was simulated by incoherently adding 106 diffraction patterns of individual discontinuous double helices. Each individual diffraction pattern was rotated around the direction of the incident beam and the rotation angles were Gaussian distributed with a certain the standard deviation In this manner, diffraction patterns for were simulated. Each simulated diffraction pattern was compared to the experimental diffraction pattern (optimized Photo 51) by calculating The least Error was found for the simulated rotationally blurred diffraction pattern with which we name RBDP0. The corresponding sample arrangement and the simulated diffraction pattern RBDP0 are shown in Fig. 13(c) – 13(d). Thus, the DNA molecules in the fiber imaged in Photo 51 were rotated around their axis, where the rotation angles were Gaussian distributed with a standard deviation of .
Appendix H Rotational motion deblurring
It is assumed that the experimental diffraction pattern of DNA fiber is an incoherent sum of the diffraction patterns of rotated individual DNA. The rotation angles, , are Gaussian distributed as with the standard deviation of 9°. There are currently no numerical routines that would allow for compensation of this kind of “rotational shaking” blur, although some algorithms for rotational deblurring have been presented in the literature [40]. To carry out such deblurring we notice that rotation of an image by an angle is equivalent to a linear shift by along angle coordinate when the image is presented in the polar -coordinates. Thus, as a first step, the rotationally blurred experimental diffraction pattern is transformed from the spherical -coordinates to the polar -coordinates, as shown in Figs. 14(a) – 14(b). Then, at each -value, 1D intensity distribution is extracted, and it is deconvoluted with the known point spread function , where The deconvolution is performed as an iterative deconvolution according to Gold’s algorithm [41], as described below. Finally, the deblurred diffraction pattern is transformed from the polar -coordinates back to the -coordinates, as shown in Figs. 14(d) – 14(e).
Gold’s iterative deconvolution algorithm
The iterative loop includes the following steps [41]:
Here is the blurred distribution, is the deblurred distribution at the i-th iteration, and is the point-spread function. The division in step (ii) is calculated aswhere is a small addend to avoid division by zero,Figure 15 shows the results of the iterative deblurring applied to the simulated rotationally blurred diffraction pattern RBDP0. It demonstrates that the many diffraction peaks (layer lines) can be recovered although they still exhibit a somewhat blurry appearance.
Appendix I Reconstruction of the rotationally blurred angular-averaged diffraction pattern
Figure 16 demonstrates that a 3D angular-averaged diffraction pattern, obtained from a 2D rotationally blurred diffraction pattern, can in principle be reconstructed, but only a fraction of the sample can be retrieved. This can be explained as follows: a finite section of the molecule which is close to the center of rotation occupies the same position as in the same but not rotated sample.
Funding
University of Zurich.
Acknowledgments
Financial support of the University of Zurich is acknowledged.
References
1. J. D. Watson and F. H. C. Crick, “Molecular structure of nucleic acids; a structure for deoxyribose nucleic acid,” Nature 171(4356), 737–738 (1953). [CrossRef] [PubMed]
2. M. H. F. Wilkins, A. R. Stokes, and H. R. Wilson, “Molecular structure of deoxypentose nucleic acids,” Nature 171(4356), 738–740 (1953). [CrossRef] [PubMed]
3. R. E. Franklin and R. G. Gosling, “Molecular configuration in sodium thymonucleate,” Nature 171(4356), 740–741 (1953). [CrossRef] [PubMed]
4. D. A. Marvin and C. Nave, “X-ray fiber diffraction,” in Structural Molecular Biology: Methods and Applications, D. B. Davies, W. Saenger, and S. S. Danyluk, eds. (Springer New York, Boston, MA, 1982), pp. 3–44.
5. K. C. Holmes, “50 years of fiber diffraction,” J. Struct. Biol. 170(2), 184–191 (2010). [CrossRef] [PubMed]
6. D. A. Marvin, “Fibre diffraction studies of biological macromolecules,” Prog. Biophys. Mol. Biol. 127, 43–87 (2017). [CrossRef] [PubMed]
7. J. W. Miao, P. Charalambous, J. Kirz, and D. Sayre, “Extending the methodology of X-ray crystallography to allow imaging of micrometre-sized non-crystalline specimens,” Nature 400(6742), 342–344 (1999). [CrossRef]
8. J. Miao, K. O. Hodgson, T. Ishikawa, C. A. Larabell, M. A. LeGros, and Y. Nishino, “Imaging whole Escherichia coli bacteria by using single-particle x-ray diffraction,” Proc. Natl. Acad. Sci. U.S.A. 100(1), 110–112 (2003). [CrossRef] [PubMed]
9. D. Shapiro, P. Thibault, T. Beetz, V. Elser, M. Howells, C. Jacobsen, J. Kirz, E. Lima, H. Miao, A. M. Neiman, and D. Sayre, “Biological imaging by soft X-ray diffraction microscopy,” Proc. Natl. Acad. Sci. U.S.A. 102(43), 15343–15346 (2005). [CrossRef] [PubMed]
10. C. Song, H. Jiang, A. Mancuso, B. Amirbekian, L. Peng, R. Sun, S. S. Shah, Z. H. Zhou, T. Ishikawa, and J. Miao, “Quantitative imaging of single, unstained viruses with coherent x rays,” Phys. Rev. Lett. 101(15), 158101 (2008). [CrossRef] [PubMed]
11. G. J. Williams, E. Hanssen, A. G. Peele, M. A. Pfeifer, J. Clark, B. Abbey, G. Cadenazzi, M. D. de Jonge, S. Vogt, L. Tilley, and K. A. Nugent, “High-resolution X-ray imaging of Plasmodium falciparum-infected red blood cells,” Cytometry A 73A(10), 949–957 (2008). [CrossRef] [PubMed]
12. X. Huang, J. Nelson, J. Kirz, E. Lima, S. Marchesini, H. Miao, A. M. Neiman, D. Shapiro, J. Steinbrener, A. Stewart, J. J. Turner, and C. Jacobsen, “Soft X-ray diffraction microscopy of a frozen hydrated yeast cell,” Phys. Rev. Lett. 103(19), 198101 (2009). [CrossRef] [PubMed]
13. Y. Nishino, Y. Takahashi, N. Imamoto, T. Ishikawa, and K. Maeshima, “Three-dimensional visualization of a human chromosome using coherent X-ray diffraction,” Phys. Rev. Lett. 102(1), 018101 (2009). [CrossRef] [PubMed]
14. J. Nelson, X. Huang, J. Steinbrener, D. Shapiro, J. Kirz, S. Marchesini, A. M. Neiman, J. J. Turner, and C. Jacobsen, “High-resolution x-ray diffraction microscopy of specifically labeled yeast cells,” Proc. Natl. Acad. Sci. U.S.A. 107(16), 7235–7239 (2010). [CrossRef] [PubMed]
15. R. N. Wilke, M. Priebe, M. Bartels, K. Giewekemeyer, A. Diaz, P. Karvinen, and T. Salditt, “Hard X-ray imaging of bacterial cells: nano-diffraction and ptychographic reconstruction,” Opt. Express 20(17), 19232–19254 (2012). [CrossRef] [PubMed]
16. M. M. Seibert, T. Ekeberg, F. R. N. C. Maia, M. Svenda, J. Andreasson, O. Jönsson, D. Odić, B. Iwan, A. Rocker, D. Westphal, M. Hantke, D. P. DePonte, A. Barty, J. Schulz, L. Gumprecht, N. Coppola, A. Aquila, M. Liang, T. A. White, A. Martin, C. Caleman, S. Stern, C. Abergel, V. Seltzer, J.-M. Claverie, C. Bostedt, J. D. Bozek, S. Boutet, A. A. Miahnahri, M. Messerschmidt, J. Krzywinski, G. Williams, K. O. Hodgson, M. J. Bogan, C. Y. Hampton, R. G. Sierra, D. Starodub, I. Andersson, S. Bajt, M. Barthelmess, J. C. H. Spence, P. Fromme, U. Weierstall, R. Kirian, M. Hunter, R. B. Doak, S. Marchesini, S. P. Hau-Riege, M. Frank, R. L. Shoeman, L. Lomb, S. W. Epp, R. Hartmann, D. Rolles, A. Rudenko, C. Schmidt, L. Foucar, N. Kimmel, P. Holl, B. Rudek, B. Erk, A. Hömke, C. Reich, D. Pietschner, G. Weidenspointner, L. Strüder, G. Hauser, H. Gorke, J. Ullrich, I. Schlichting, S. Herrmann, G. Schaller, F. Schopper, H. Soltau, K.-U. Kühnel, R. Andritschke, C.-D. Schröter, F. Krasniqi, M. Bott, S. Schorb, D. Rupp, M. Adolph, T. Gorkhover, H. Hirsemann, G. Potdevin, H. Graafsma, B. Nilsson, H. N. Chapman, and J. Hajdu, “Single mimivirus particles intercepted and imaged with an x-ray laser,” Nature 470(7332), 78–81 (2011). [CrossRef] [PubMed]
17. T. Ekeberg, M. Svenda, C. Abergel, F. R. N. C. Maia, V. Seltzer, J.-M. Claverie, M. Hantke, O. Jönsson, C. Nettelblad, G. van der Schot, M. Liang, D. P. DePonte, A. Barty, M. M. Seibert, B. Iwan, I. Andersson, N. D. Loh, A. V. Martin, H. Chapman, C. Bostedt, J. D. Bozek, K. R. Ferguson, J. Krzywinski, S. W. Epp, D. Rolles, A. Rudenko, R. Hartmann, N. Kimmel, and J. Hajdu, “Three-dimensional reconstruction of the giant mimivirus particle with an X-ray free-electron laser,” Phys. Rev. Lett. 114(9), 098102 (2015). [CrossRef] [PubMed]
18. L. W. Whitehead, G. J. Williams, H. M. Quiney, D. J. Vine, R. A. Dilanian, S. Flewett, K. A. Nugent, A. G. Peele, E. Balaur, and I. McNulty, “Diffractive imaging using partially coherent x rays,” Phys. Rev. Lett. 103(24), 243902 (2009). [CrossRef] [PubMed]
19. J. C. H. Spence, U. Weierstall, and H. N. Chapman, “X-ray lasers for structural and dynamic biology,” Rep. Prog. Phys. 75(10), 102601 (2012). [CrossRef] [PubMed]
20. R. G. Gosling, “X-ray diffraction studies of desoxyribose nucleic acid,” PhD Thesis, King's College London (1954).
21. H. R. Drew, R. M. Wing, T. Takano, C. Broka, S. Tanaka, K. Itakura, and R. E. Dickerson, “Structure of a B-DNA dodecamer: conformation and dynamics,” Proc. Natl. Acad. Sci. U.S.A. 78(4), 2179–2183 (1981). [CrossRef] [PubMed]
22. J. W. Miao, D. Sayre, and H. N. Chapman, “Phase retrieval from the magnitude of the Fourier transforms of nonperiodic objects,” J. Opt. Soc. Am. A 15(6), 1662–1669 (1998). [CrossRef]
23. J. R. Fienup, “Phase retrieval algorithms: a comparison,” Appl. Opt. 21(15), 2758–2769 (1982). [CrossRef] [PubMed]
24. S. Marchesini, H. He, H. N. Chapman, S. P. Hau-Riege, A. Noy, M. R. Howells, U. Weierstall, and J. C. H. Spence, “X-ray image reconstruction from a diffraction pattern alone,” Phys. Rev. B Condens. Matter Mater. Phys. 68(14), 140101 (2003). [CrossRef] [PubMed]
25. J. C. H. Spence, R. A. Kirian, X. Wang, U. Weierstall, K. E. Schmidt, T. White, A. Barty, H. N. Chapman, S. Marchesini, and J. Holton, “Phasing of coherent femtosecond X-ray diffraction from size-varying nanocrystals,” Opt. Express 19(4), 2866–2873 (2011). [CrossRef] [PubMed]
26. K. Ayyer, O. M. Yefanov, D. Oberthür, S. Roy-Chowdhury, L. Galli, V. Mariani, S. Basu, J. Coe, C. E. Conrad, R. Fromme, A. Schaffer, K. Dörner, D. James, C. Kupitz, M. Metz, G. Nelson, P. L. Xavier, K. R. Beyerlein, M. Schmidt, I. Sarrou, J. C. H. Spence, U. Weierstall, T. A. White, J. H. Yang, Y. Zhao, M. Liang, A. Aquila, M. S. Hunter, J. S. Robinson, J. E. Koglin, S. Boutet, P. Fromme, A. Barty, and H. N. Chapman, “Macromolecular diffractive imaging using imperfect crystals,” Nature 530(7589), 202–206 (2016). [CrossRef] [PubMed]
27. T. Latychevskaia, “Iterative phase retrieval in coherent diffractive imaging: practical issues,” Appl. Opt. 57(25), 7187–7197 (2018). [CrossRef] [PubMed]
28. G. Friedel, “Sur les symétries cristallines que peut révéler la diffraction des rayons Rontgen,” Cr Acad Sci C Chim 157, 1533–1536 (1913).
29. H. N. Chapman, A. Barty, S. Marchesini, A. Noy, S. P. Hau-Riege, C. Cui, M. R. Howells, R. Rosen, H. He, J. C. H. Spence, U. Weierstall, T. Beetz, C. Jacobsen, and D. Shapiro, “High-resolution ab initio three-dimensional x-ray diffraction microscopy,” J. Opt. Soc. Am. A 23(5), 1179–1200 (2006). [CrossRef] [PubMed]
30. H. B. Chang, Y. F. Lou, M. K. Ng, and T. Y. Zeng, “Phase retrieval from incomplete magnitude information via total variation regularization,” SIAM J. Sci. Comput. 38(6), A3672–A3695 (2016). [CrossRef]
31. J. M. Zuo, I. Vartanyants, M. Gao, R. Zhang, and L. A. Nagahara, “Atomic resolution imaging of a carbon nanotube from diffraction intensities,” Science 300(5624), 1419–1421 (2003). [CrossRef] [PubMed]
32. W. Cochran, F. H. C. Crick, and V. Vand, “The structure of synthetic polypeptides. I. The transform of atoms on a helix,” Acta Crystallogr. A 5(5), 581–586 (1952). [CrossRef]
33. H.-L. Ng, M. L. Kopka, and R. E. Dickerson, “The structure of a stable intermediate in the A ←> B DNA helix transition,” Proc. Natl. Acad. Sci. U.S.A. 97(5), 2035–2039 (2000). [CrossRef] [PubMed]
34. NIST, “NIST X-ray form factors table ” https://physics.nist.gov/PhysRefData/FFast/html/form.html.
35. R. R. Nair, P. Blake, J. R. Blake, R. Zan, S. Anissimova, U. Bangert, A. P. Golovanov, S. V. Morozov, A. K. Geim, K. S. Novoselov, and T. Latychevskaia, “Graphene as a transparent conductive support for studying biological molecules by transmission electron microscopy,” Appl. Phys. Lett. 97(15), 153102 (2010). [CrossRef]
36. C. Seuring, K. Ayyer, E. Filippaki, M. Barthelmess, J. N. Longchamp, P. Ringler, T. Pardini, D. H. Wojtas, M. A. Coleman, K. Dörner, S. Fuglerud, G. Hammarin, B. Habenstein, A. E. Langkilde, A. Loquet, A. Meents, R. Riek, H. Stahlberg, S. Boutet, M. S. Hunter, J. Koglin, M. Liang, H. M. Ginn, R. P. Millane, M. Frank, A. Barty, and H. N. Chapman, “Femtosecond X-ray coherent diffraction of aligned amyloid fibrils on low background graphene,” Nat. Commun. 9(1), 1836 (2018). [CrossRef] [PubMed]
37. M. Guizar-Sicairos, S. T. Thurman, and J. R. Fienup, “Efficient subpixel image registration algorithms,” Opt. Lett. 33(2), 156–158 (2008). [CrossRef] [PubMed]
38. R. E. Franklin and R. G. Gosling, “The structure of Sodium Thymonucleate Fibres 2. The cylindrically symmetrical Patterson function,” Acta Crystallogr. 6(8), 678–685 (1953). [CrossRef]
39. T. Latychevskaia and H.-W. Fink, “Practical algorithms for simulation and reconstruction of digital in-line holograms,” Appl. Opt. 54(9), 2424–2434 (2015). [CrossRef] [PubMed]
40. Q. Shan, W. Xiong, and J. Y. Tia, “Rotational Motion deblurring of a rigid object from a single image,” IEEE 11th International Conference on Computer Vision. 16, 738–745 (2007). [CrossRef]
41. R. Gold, “An iterative unfolding method for matrices,” in Tech. Rep. USAEC ANL-6984 Argonne National Laboratory(Argonne, 1964).