research papers
Shannon sampling based approach for the structural solution of nano-objects by laboratory and synchrotron SAXS data
aInstitute of Crystallography, National Research Council, Bari, 70126, Italy, and bExcelsus Structural Solutions (Swiss) AG, Park Innovaare, Villigen, 5234, Switzerland
*Correspondence e-mail: [email protected], [email protected]
We propose a Shannon sampling based approach for the analysis of laboratory and synchrotron small-angle X-ray scattering (SAXS) data. This method avoids calculation of the pair distance distribution function P(r), which is typically used as an alternative to the model-dependent fitting of SAXS data for determining the size and shape of the scattering object. Using this new approach it is also possible to skip any model-dependent fit of the SAXS data. The method has been tested on SAXS data collected from Polysorbate 20, D-α-tocopheryl polyethylene glycol 1000 succinate spheroidal core–shell micelles and relevant well known proteins (lysozyme, apoferritin, carbonic anhydrase 2) of different molecular weights. The results confirm the potential of the approach, even for laboratory scattering data and in the presence of polydispersity. The method allows the reliable derivation of maximum size, mass, volume, polar/equatorial asymmetry, hydration fraction, shell thickness and core/shell electron-density differences for core–shell structures, directly in the q domain, working on SAXS data measured from dilute samples for which interparticle structure factors are negligible.
Keywords: SAXS; small-angle X-ray scattering; PS20 micelles; VitE-TPGS micelles; Nyquist–Shannon sampling; proteins.
1. Introduction
Small-angle X-ray scattering (SAXS) is a non-destructive, indispensable tool for quantitative nanoscale characterization of soft matter and complex fluids, including micelles, colloids and macromolecules. By meticulously controlling sample preparation – typically using thin-walled capillaries to minimize background and ensure optimal transmission – and employing radiation-controlled exposure strategies, SAXS yields robust structural/morphological information on the scattering nano-objects under investigation, dispersed in a buffer. Data can be collected either on a relative scale, i.e. , or on an absolute scale, as
cm−1, where the normalization to water scattering is based on the isothermal compressibility of water. The value 0.0164 cm−1 corresponds to the theoretical absolute scattering intensity of pure water at 298 K and atmospheric pressure (Orthaber et al., 2000
; Zeng et al., 2017
; Roe, 2000
).
Data interpretation in SAXS relies on the form factor, which describes the scattering from an individual particle, combined with the structure factor – when interparticle interactions are present – to account for those interactions. When interparticle effects are minimized for sufficiently dilute sample solutions, the scattering curve primarily reflects the form factor and, therefore, the shape of the particle. For the form factor, analytical expressions are available to compute the scattering intensity and describe the object shape (sphere, disc, rod, dumbbell, core–shell etc.), often incorporating polydispersity parameters such as size, shape and shell thickness.
The form factor modeling approach has successfully allowed the morphological characterization of spheroidal core–shell micelles, enabling the determination of micellar core size, shell thickness and electron-density contrast. This required the definition of a priori assumptions concerning particle shape, polydispersity, uniformity, layer sharpness and interparticle interactions (Pedersen, 1997
; Glatter & Kratky, 1982
; Castelletto & Hamley, 2002
). In practice, some of these a priori assumptions may not hold, which can invalidate the determination of key structural parameters of unknown scattering nano-objects from SAXS data.
Alternatively, a second class of methods can be used, which builds on the pair distance distribution function computed through indirect Fourier transform of the SAXS data. These techniques avoid shape assumptions but still yield insights into the maximum particle size Dmax, radius of gyration
and internal inhomogeneities (Svergun, 1992
; Hansen, 2000
; Vestergaard & Hansen, 2006
). Within this second class, a theoretical calculation of , compared with the evaluation of experimental
moments, was recently proposed as a possible alternative (De Caro et al., 2025
).
More advanced frameworks employ hybrid or simulation-based models, combining SAXS data with molecular dynamics or Monte Carlo simulations. These methods generate structural ensembles consistent with scattering patterns and capture effects such as dynamic fluctuations, anisotropy in the shell region, or non-ideal solvation (Stuhrmann, 2008
; Chen et al., 2023
; Cezar et al., 2025
). Moreover, contrast variation and global fitting approaches integrate SAXS with complementary scattering techniques (e.g. small-angle neutron scattering, SANS) or with multiple experimental datasets at once, thereby improving robustness and reducing the degeneracy of solutions (König et al., 2020
). In fact, in SAXS model fitting, unconstrained least-squares minimization often yields degenerate solutions. To stabilize and obtain physically meaningful results, one needs to add constraints to the χ2 minimization (Svergun, 1992
; Pedersen, 1997
; Stuhrmann, 2008
; Breßler et al., 2015
; König et al., 2020
).
In the present work, we propose a new method based on the Nyquist–Shannon (NS) sampling theorem. The new approach works directly on the SAXS data, as collected in the q domain, without any need for calculating distributions of distances in real space, namely the pair distance distribution function . The method proved to be effective for deriving many key parameters from SAXS data measured on dilute samples, for which interparticle structure factors are negligible: maximum size of the scattering particle, Dmax, polar/equatorial asymmetry (ellipticity), volume, mass and number of monomers when needed. We successfully tested this new approach on Polysorbate 20 (PS20) and vitamin E D-α-tocopheryl polyethylene glycol succinate (VitE-TPGS) micelles, using both laboratory and synchrotron SAXS data, and on relevant well known proteins (lysozyme, apoferritin, carbonic anhydrase 2) from synchrotron SAXS data. Moreover, we have checked the reliability of the structural parameters derived by the new approach by also determining them by means of a model-dependent least-squares χ2 minimization of the difference profile between experimental and calculated SAXS data. Using the new proposed approach, it is straightforward to derive several key structural parameters of unknown scattering nano-objects from SAXS data; it could also be thought of as a preliminary and quick way to obtain a first estimate of these parameters, before performing time-consuming least-squares χ2 minimization of the difference profile between experimental and calculated SAXS data. Additionally, it could be considered a useful independent and indirect check of the a priori assumptions concerning particle shape, polydispersity, uniformity, layer sharpness etc., which usually characterize SAXS modeling approaches.
2. Theoretical framework
In the following sections we will present the new approach, based on the NS sampling of SAXS data, reporting the main formulae and leaving further mathematical details to the supporting information.
2.1. Determination of the maximum size of the scattering particle in q space
The aim of this section is to demonstrate that it is possible to obtain a reliable estimate of Dmax, the maximum size of the scattering particle from SAXS data, without calculating , through a new approach based on the NS sampling of SAXS intensities.
Moore (1980
) used the Shannon sampling theorem (Shannon, 1948
) to sample the SAXS profile at a selection of q values, called the Shannon channels, here denoted as NS channels (Moore, 1980
; Grant, 2022
), defined by
with n a positive integer and Dmax the maximum particle dimension.
The SAXS profile can be described as an infinite sum of the intensities sampled at qn, i.e.
(Grant, 2022
):
Application of the NS sampling to SAXS data has been shown to be advantageous because it reduces a densely measured and noisy scattering curve to a small set of independent points that fully capture the real structural information, avoiding risks of overfitting. In fact, the structural information content in the SAXS profiles is affected by rotational averaging and other physical factors, thus giving experimental access to only a limited number of NS channels (up to nmax). However, the number of experimentally available NS channels, nmax, is much lower than the typical number of measured intensities () in a SAXS experiment. A measurement step Δq ∼ 0.005 Å−1 of
yields an oversampling ratio OR =
of the information (signal) contained in the experimental SAXS intensity. This redundancy is crucial for robustly constraining the determination of NS channel intensities In, even in the presence of low signal-to-noise ratios. With typical oversampling factors of 50−100, the information content within the SAXS profile is significantly redundant; this compensates for the low signal-to-noise ratio inherent in the determination of the NS channel intensities, in particular for noisy laboratory measurements.
Moreover, considering that the scattering intensity decreases rapidly with increasing magnitude of the scattering vector q (typically following a high negative power of q), often also very few NS channels are sufficient to calculate reliably average quantities determined by the SAXS intensity. In fact, the pair distance distribution function
, calculated as a Fourier transform of a truncated version of the infinite series in equation (2)
, containing only a finite number of terms In, i.e. those experimentally available, suffers very often from severe ripples caused by Fourier termination effects. These ripples significantly affect the Dmax determination, i.e. where goes to zero (Grant, 2022
).
Nevertheless, the NS sampling approach offers the possibility to determine the forward scattering intensity , the average vector length r1 in the nanoparticle and the gyration ratio
, evaluated as sums over the NS channel intensities (Grant, 2022
), up to the In experimentally available. These can be defined by
These quantities are, respectively, moments of order 0, 1 and 2 of . The
moments are mean values, related to integrals of the function
multiplied by some power of r. Being integral quantities, they are much less affected by truncation effects due to the finite nmax experimentally available.
Equations (3)–(5) allow structural quantities to be derived directly from the SAXS profiles measured on dilute samples, for which interparticle structure factors are negligible, if Dmax has been previously estimated. The (Guinier & Fournet, 1955
) permits evaluation of and Rg directly from the experimental SAXS profile at very low q (in the asymptotic limit q →0). These values are related to a low-spatial-resolution approximation of the shape of the scattering particle, therefore requiring a small number of terms in the sum expressions of equations (3)–(5).
Thus, in that case, starting from equations (3) and (5), we can write the following two equations:
where is the normalized gyration radius, and
Here, `mis' refers to `missing' intensities, i.e. related to SAXS intensity values which are unobserved due to the noise affecting measurements, and the subscript gui refers to the However, for the SAXS intensities are orders of magnitude smaller than the measured values for
. For scattering nano-objects we usually have a
scaling of intensities, as a function of q, with n larger than 1, up to 4 for ideally well defined surfaces of globular objects (Porod law). Moreover, there are `alternating signs' in the numerical series. Therefore,
and
are expected to be very small.
For this reason, they can also be estimated assuming a power-law q scaling of the unobserved intensities:
with , or different suitable powers that could be chosen in the case of additional information about the shape or other structural properties of the scattering nano-objects, or after having calculated the Kratky plot (Glatter & Kratky, 1982
). The upper limit in the sums can be chosen arbitrarily. The value of 1000 in the above equations is equivalent to having an infinite upper limit in the sums. Our experimental tests have demonstrated that the small contributions and
respectively, in equations (6) and (7), can also be very often neglected, i.e. the truncation at nmax of equations (6) and (7) already gives reliable results of the
moments (see Table 1
).
‡The SAXS intensity for SASDFN8, reported in the SASBDB repository, has been multiplied by 0.0001 to express it in cm−1. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
By varying Dmax, e.g. in small steps of 0.1 Å, a value usually smaller than the experimental error with which Dmax can be determined, we obtain from equations (6) and (7) two dimensionless functions, and
=
, whose minima give two possible values, Dmax,I and Dmax,r, of the unknown maximum size. We can search for the two minima by changing Dmax in the range from 2.5
to 3.5
, to obtain Dmax,I and Dmax,r which minimize equation (6)
and equation (7)
, respectively, or in other suitable intervals, if needed, when the minima fall out of the above interval. To find these minima, smoothed intensity values are used to reduce the effect of noise, an operation allowed due to the large oversampling of the SAXS signal. For this purpose, q-adaptive filtering sigmoidal kernels, described in the supporting information (Section S2), can be used. This process leaves lower-q SAXS intensities unchanged, averaging only the noisy region of the experimental data at higher q. Utilizing mean-filtered SAXS profiles ensures statistically robust evaluations of NS channel intensities. The large oversampling allows for local averaging around each qn point, thereby incorporating neighboring intensities and enabling the extraction of reliable In values from the noisy data.
The experimental profile is measured on a discrete set of scattering vector magnitudes qi, with i ∈ (1, imax). For any considered Dmax, the smoothed SAXS profile will usually not contain an intensity value I( qi) located exactly at the qi value of the nth NS channel, i.e. corresponding to . We can therefore use the intensity value from the smoothed SAXS profile whose qi is closest to
, i.e. the one minimizing
. The step size of 0.1 Å may be too fine for identifying the minima of the functions
and
=
. Therefore, both functions also usually require smoothing, e.g. with a mean filter. Different levels of smoothing of
and
=
may yield slightly different solutions, Dmax,I and Dmax,r. Among all the solutions obtained by varying the width of the mean filter, we can select only the values of Dmax,I and Dmax,r that satisfy the following condition:
i.e. the values of Dmax that minimize simultaneously both and
=
.
As an example, Fig. 1
shows the functions (black curve) and
=
(red curve) evaluated for laboratory SAXS data collected from PS20 micelles. The experimental conditions used for measuring the SAXS data are described in Section 3
.
| Figure 1 |
By changing the level of smoothing of the two functions, we could find several solutions satisfying both equations (6)
and (7) and the constraint given by equation (12)
. The average values and
of all minima of equations (6)
and (7), satisfying also the condition given by equation (12)
, will give the maximum size of the nano-objects under study, without the need to calculate the pair distance distribution function:
When the Guinier plot is not accurately evaluated, e.g. for inaccurate experimental measurements at low q, due to incorrect sample preparation or experimental settings, we have no solution satisfying equation (12)
, because the two minima are not aligned in terms of abscissa coordinates, i.e. Dmax,I is too different from Dmax,r. Therefore, the result of having the two minima perfectly aligned, as shown in the example of Fig. 1
, is an indirect confirmation that the Guinier plot has been correctly computed, i.e. sample and experiment preparation have been correctly performed. When this condition is not satisfied, the should be recalculated, usually eliminating some SAXS intensities at very low q, leading to solutions satisfying the constraint given by equation (12)
, or alternatively repeating sample preparation or data collection. This finding provides an internal check that Dmax, and
have been consistently determined with respect to each other. Every time the SAXS data allow a meaningful to be obtained, giving
and
values truly representative of the scattering nano-objects, the approach described above gives the possibility to determine a Dmax value consistent with these values.
Table 1
shows the results obtained for Dmax for several proteins and core–shell micelles, using laboratory and synchrotron SAXS data. For comparison, we report also the gyration radius and the I(0) value determined by the described approach based on the NS sampling theorem, as well as the values reported in the Small Angle Scattering Biological Data Bank (SASBDB) repository (https://www.sasbdb.org/; Valentini et al., 2015
), obtained by other software (GNOM) (Svergun, 1992
). Although different approaches give different Dmax, the differences are within the experimental error bars. Thus, we can conclude that there is a very good agreement.
Moreover, the estimation of the truncation errors due to the finite maximum q value experimentally available, reported in Table 1
, shows that they usually are less than 1%. In any case, we have used these small corrections only for the determination of Dmax. Conversely, the values of I(0) and Rg reported in Table 1
have been obtained without considering these small corrections due to finite nmax. Even without these corrections, the results shown in Table 1
demonstrate that the predicted values for these quantities obtained by the NS sampling formulae agree with the values obtained by the within the experimental errors. These results show that truncation of the series does not significantly affect the P(r) moments I(0) and , as these quantities are related to integrals of P(r).
2.2. Determination of the mass and polar/equatorial asymmetry of core–shell structures
The aim of this section is to demonstrate that it is possible to obtain relevant structural information on core–shell structures (mass, volume, ellipticity) through the new approach, based on the NS sampling of SAXS data.
The volume V of a protein of a nearly globular shape, proportional to the ratio between the scattering intensity at q = 0, I(0), and the Porod invariant Q = , can, in principle, be determined directly from SAXS data, even collected in relative units (Fischer et al., 2010
):
The protein molecular weight can be obtained from the volume V (Å3), given by equation (14)
, multiplied by the average density 0.83 × 10−3 kDa Å−3 for proteins in solution, equivalent to a ρp = 1.37 g cm−3 (Fischer et al., 2010
). The calculation of the Porod invariant Q can be performed up to the qmax accessible in the SAXS data, giving an apparent protein volume , instead of the correct volume V, under the hypothesis of a protein globular shape. However, a linear relationship has been derived by Fischer et al. (2010
), to obtain the correct volume from the apparent value :
where V, A and are expressed in Å3, and B is dimensionless. The coefficients of the linear regression, A and B, depend on the maximum scattering vector magnitude
used for the calculation of the invariant Q (i.e. the integral of the Kratky plot) (Fischer et al., 2010
). We note that qmax corresponding to the maximum NS channel available in the SAXS data is usually different from , used for the calculation of the Porod invariant Q.
can be estimated by comparing
at any q and the mean value
, as shown, for example, in the supporting information for apoferritin (SASDFN8) (see Section S4). The scattering vector value where
allows one to estimate
, i.e. where the level of noise prevents access to further structural information from the
profile.
Applying the NS formalism, the following expressions for , V and the molecular weight Mp can be derived (Grant, 2022
):
Here we apply a linear fit of A and B values reported by Fischer et al. (2010
), where has the dimension of a volume (Å3) and
is dimensionless. The linear fit gives
Corrections like , related to the finite nmax of the experimental SAXS data, are not needed in equations (17)–(18), because the corrections due to the finite q range are already considered through the correction
, by means of equations (18)–(20).
The determination of the maximum size, Dmax, and of the Shannon channel intensities obtained with the approach described in the previous section allows us to determine the molecular weight for biological macromolecules and proteins.
For estimating the molecular weight of macromolecules another invariant can be introduced, named the correlation volume but with a surface (length2) unit,
(Rambo & Tainer, 2013
), which is useful when the Q invariant diverges at high scattering vectors in the The use of this invariant leads to the following mass formula, which is concentration independent, holds also for SAXS data in relative units and gives the nano-object volume () per self-correlation length (
) (Rambo & Tainer, 2013
):
Here, is the effective protein volume, which could be different from V obtained by the Porod invariant [equation (17)
], due to its flexibility. Equation (22)
has been derived by Rambo & Tainer (2013
) through a comparison of the new invariant QR with the mass of many proteins taken from the Protein Data Bank (https://www.rcsb.org/). as obtained by equation (22)
is given in kDa if is expressed in Å2 and
in Å. For mixed protein–nucleic acid complexes or RNA alone, other formulae can be derived, with different fitting power-law relationships (Rambo & Tainer, 2013
).
The NS formalism also allows us to express the above quantities in terms of sums on the NS channels (Grant, 2022
):
) is the sine integral function calculated at
. In equations (23)
and (24) we have approximated the infinite sum up to nmax, because equation (22)
already takes into account the effects of the finite maximum value of q, since it has been obtained by the fit over thousands of macromolecules always characterized by a finite qmax (Rambo & Tainer, 2013
).
For spherical macromolecules we have . For spheroidal shapes with sufficiently large polar/equatorial asymmetry and for flexible macromolecules, we have
. In principle, given
and
,
could give indications about the internal spatial correlation, i.e. the effective extension and the coherence of the scattering of any nanoparticle. Instead, equation (17)
gives the effective volume, corrected by the effect of a finite qmax. This effective volume is that of an equivalent sphere with the same scattering as the nanoparticle. Furthermore, given V, if the particle is not a protein, we could obtain its mass from equation (18)
by inserting the correct value of the particle-mass density, substituting 1.37 g cm−3 for the of the molecules constituting the particle, i.e. generalizing equation (18)
as follows:
These considerations allow us to extend the range of application of the above formulae to predict the mass of core–shell spheroidal micelles. In the supporting information we have shown how to derive the polar/equatorial asymmetry (ellipticity) directly by the comparison of and
. The ratio of
and
for spheroidal nano-objects, where
can be determined from equation (25)
, can be related just to the polar/equatorial asymmetry of its shape. When applied to core–shell structures [equation (S37) in Section S5 of the supporting information], in the case of prolate shapes, this ratio gives simply the particle ellipticity; for oblate shapes this ratio gives :
We have tested this approach on:
(i) Well known proteins (lysozyme, apoferritin, carbonic anhydrase 2) of different molecular weights, whose synchrotron SAXS data were downloaded from the SASBDB repository.
(ii) The core–shell nano-objects studied in this work, i.e. apoferritin (SASDFN8), PS20 and VitE-TPGS micelles (laboratory and synchrotron SAXS data).
Results are summarized in Table 2
. We observe a good agreement between and M (experimental mass value reported in the SASBDB repository) for the proteins, and between
and
for core–shell structures (apoferritin and micelles only), with the values derived from a model-dependent fit of the experimental SAXS data, described in detail in the following.
|
Conversely, for the proteins that do not have a core–shell structure, the ratio of the two masses is an indication of the actual internal scattering correlation compared with the value corresponding to an ideal compact scattering object, assumed for the calculation of .
2.3. Determination of the volume, hydration fraction and aggregation number of core–shell structures
The aim of this section is to demonstrate that it is possible to determine the number of monomers or sub-units constituting the core–shell structure through the new approach based on the Shannon sampling formalism of SAXS data, combined with the volume of the scattering particle and its hydration fraction. All these quantities can be derived without the need to perform model-dependent simulations.
The knowledge of the polar/equatorial asymmetry, obtained by equation (26)
(penultimate column of Table 2
), and of the maximum size (third column of Table 1
) for spheroidal core–shell structures allows the calculation of its volume. Moreover, Rambo & Tainer (2013
) have shown that, for mixed protein–nucleic acid complexes or single RNA molecules, other formulae than equation (22)
can be derived, with different fitting power-law relationships of . We could verify if, for core–shell structures, there is a power-law relation of
suitable to derive the non-hydrated mass of the nanostructure and, consequently, the number of monomers constituting the core–shell structures. We would expect that the internal scattering coherence (related to
) should be quantitatively related to the number of monomers
packed in the structure, because in spheroidal micelles the monomers do not pack necessarily with an isotropic order. We have already shown that for core–shell nanostructures the ratio
is related to the ellipticity [equation (26)
]. Therefore, the non-hydrated mass of core–shell structures and, consequently, the number of monomers
could be derived from the following equation:
where is the molar weight mass of the monomer constituting the core–shell structure, and
is given by equation (25)
. The above formula is related to quantities and
, which can be derived also by SAXS data on a relative scale. Therefore, equation (27)
could be useful to estimate the aggregation number from SAXS data even collected on a relative scale.
For SAXS intensity on an absolute scale an independent evaluation of the aggregation number can be determined as
where is the q→0 value of the SAXS intensity of one monomer and can be estimated as follows (De Caro et al., 2024
):
c is the concentration of monomers (g cm−3), the Avogadro number,
the classical electron radius, SLD the scattering length density for the monomer and the water molecule (9.44 × 1010 cm−2),
the monomer volume,
the total number of electrons in a monomer,
the total number of electrons in a water molecule, and
the volume of a water molecule. The comparison between the aggregation numbers obtained by the two independent formulae, i.e. equation (27)
versus equations (28)–(29), is useful for validating the prediction from the new formula [equation (27)
].
Moreover, we note that the aggregation numbers derived from the above formulae are independent of the particle volume. By using equation (27)
, the hydrated fraction of the core–shell volume can be calculated from the ratio between the volume of monomers and the volume of the core–shell structure, where
can be readily derived from Dmax and ɛM [equation (26)
]:
The results obtained by the above formulae are summarized in Table 3
.
‡For the VitE-TPGS monomer we have assumed the following chemical formula: hydrophobic chain + succinate linker + PEG 1000 (n = 23): C29H50O2 + C4H4O3 + C46H92O23 = C79H146O28. | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Note that we obtain hydrated volume fractions for the micelles that are too high for the oblate shape compared with the prolate values. This finding implies that the prolate shape is the correct one. For the apoferritin, we obtain a spherical shape, as expected. The new method proposed here for determining prolate/oblate shapes of spheroidal micelles can be considered an alternative approach to existing ones (e.g. Vass et al., 2008
).
Furthermore, the results of Table 3
show that the estimation of the aggregation number obtained from the new formula proposed here [equation (27)
] seems to give reliable numbers of monomers constituting the core–shell nanostructures, independently of those estimated through equations (28)–(29). Within the errors, equation (27)
gives the correct number of sub-units, N = 24, constituting apoferritin.
Moreover, the comparison with Nagg = 127 ± 2 obtained for a sample of VitE-TPGS of the same concentration and the same Dmax, but derived by a different approach based on the P(r) moments (De Caro et al., 2025
), is an important confirmation of the new estimations obtained by equation (27)
and equations (28)–(29) through the NS sampling approach proposed here. Also, the comparison of Nagg = 45 for PS20 with the results obtained in other work (De Caro et al., 2024
; Ehrit et al., 2023
) confirms the same conclusion about the reliability of the Nagg estimations here proposed. The Dmax of the PS20 micelles characterized here is about 1 nm larger than that of the PS20 micelles studied by De Caro et al. (2024
) and Ehrit et al. (2023
), implying a volume about 33% greater and, consequently, 33% more monomers constituting the micelles. The value of Nagg = 34 ± 1 reported by De Caro et al. (2024
) and Ehrit et al. (2023
), incremented by one-third, gives Nagg ∼ 45, i.e. the same as the value obtained here, but by completely different approaches. All these findings confirm the validity of the new proposed formula for estimating the aggregation number of core–shell nanostructures from SAXS data collected on a relative scale [equation (27)
]. Future use of this formula will help to expand and further validate its applicability.
2.4. Determination of the shell thickness of core–shell structures
The aim of this section is to demonstrate that it is possible to determine the shell thickness of core–shell structures, by solving either numerically or graphically an equation based on the balance of the volume of the monomers and the volume of the hydrated molecules inside the core–shell structures.
In the supporting information we have derived the following equation, which gives the possibility to derive the fraction [equation (S10) of the supporting information]:
The numerical solution of equation (32)
enables the determination of the shell thickness since both the shell and core volumes are a function of X [equations (S24)–(S25) and (S27) of the supporting information] and Dmax has already been determined. Fig. 2
shows as a function of X, for the PS20 laboratory data, for the prolate (blue) and oblate (red) shapes. The results for the considered core–shell structures, including apoferritin, have been summarized in Table 4
.
|
| Figure 2 ΔV/ |
The prolate solutions for the micelles shown in Table 4
are in good agreement with the results obtained by the model-dependent simulations, reported in the supporting information (Section S4).
Also, the shell size of the apoferritin, in the core–shell approximation, is correctly derived. For oblate shapes we always obtain values for shell thickness/maximum distance ratios that are too large: for the PS20 laboratory SAXS data X = 0.245 ± 0.005; for the PS20 synchrotron SAXS data X = 0.245 ± 0.005; for the VitE-TPGS laboratory SAXS data X = 0.205 ± 0.005. In all cases, the oblate solution can be discarded because it leads to shell thickness values that are too large, a check that can be readily verified if one has some additional information about the monomers constituting the micelles (e.g. estimated hydrophobic chain length etc.).
2.5. Determination of shell and core electron-density differences
The aim of this section is to demonstrate that, when SAXS data are available in absolute units, it is possible to determine the core and shell electron-density differences of core–shell structures, by solving either numerically or graphically a suitable equation based on the balance of the total number of electrons inside the volume of the scattering particles.
For a core–shell system, we have derived an equation [see equations (S46)–(S48) of Section S6 of the supporting information] that can be solved numerically to determine the ratio of the two electron-density differences =
/
.
This equation has been rewritten here in terms of :
The ratio is the solution of the above equation, once
and X have been determined, as previously described.
can also be obtained for SAXS data on a relative scale. Conversely, for determining the shell and core electron-density differences, we need SAXS data on an absolute scale. For this purpose, we can use the following equations [equations (S48)–(S49) of Section S6 of the supporting information]:
c is the monomer concentration (g cm−3), × 1023 mol−1 is the Avogadro number,
is the classical electron radius (
cm e−1),
is the molar mass of one monomer and
is the number of aggregated monomers within the core–shell structure;
is the volume of the core–shell structure and
is the core volume.
Table 5
reports the electron-density differences obtained for the core–shell structures under study. We have a good agreement between derived here and the values determined by the model-dependent fits of SAXS data, reported in Section S4 of the supporting information, although we note that the solutions of equation (33)
are systematically about 5% higher than the values obtained by the fits of SAXS data. This difference could be caused by the presence of polydispersity, as discussed in Section S4.
|
Nevertheless, we obtain a perfect agreement between the electron-density differences derived here for the VitE-TPGS micelles and those determined by De Caro et al. (2025
) for micelles of the same size (Δρshell-solution = 0.032 ± 0.001 e Å−3, Δρcore-solution = −0.026 ± 14;0.001 e Å−3). Conversely, for the PS20 micelles, if we compare the results derived here with those published elsewhere (De Caro et al., 2024
; Ehrit et al., 2023
), we note a perfect agreement for the values. Indeed, in their work we find Δρcore-solution = −0.030 ± 0.002 e Å−3. Conversely, the value determined by De Caro et al. (2024
) of Δρshell-solution = 0.050 ± 0.005 e Å−3 is quite different from the values determined here from both laboratory and synchrotron SAXS data. However, we have already stressed that the PS20 micelles characterized here are about one-third bigger than those characterized by De Caro et al. (2024
) and Ehrit et al. (2023
). Very recently, Evers et al. (2025
) have shown the great variability of the polyoxyethylene (POE) chain with respect to the nominal 20 units, for commercial PS20. For the monolaurate core we should expect a more stable chemical composition and, consequently, a stable value of the core electron-density difference, as obtained by our analysis. The variability of the length and the structure of the POE chains in commercial PS20 products (Evers et al., 2025
) could easily cause different maximum sizes, shell thickness and shell electron density, as verified in the present analysis.
Finally, the application of our approach to apoferritin leads to a very low core electron-density difference, as expected, for the iron cage function of this protein. From the FASTA sequence of the SASDFN8 apoferritin the nominal number of electrons inside one of the 24 units constituting this protein can be estimated as about 10 × 103. This value can be compared with the value obtained by the electron-density differences derived here, by means of the formula
where the last term gives the number of electrons of an equivalent number of water molecules contained in the un-hydrated volume of the protein. By inserting the electron-density differences determined here (Table 5
) and the other quantities reported in the previous tables for apoferritin, we obtain .
3. Experimental section
In this section we summarize the main experimental data of the SAXS experiments and the results of model-dependent fits of the SAXS experimental data. Further details are reported in the supporting information. The model-dependent fits described here have been performed as a further check of the reliability of the results discussed in Section 2
, which are completely independent from the results reported in the following.
3.1. Laboratory SAXS data
SAXS data were collected by means of a XEUSS3.0 HR system equipped with a motorized dual-X-ray source (Excillum Ga metal jet and Cu Genix3D micro-sources). All measurements reported here were carried out using Ga K radiation (λ = 1.34 Å). Scattering was recorded with an Eiger 2R 1M detector under standard vacuum conditions.
SAXS measurements were performed using 2 mm Biocube capillaries. For PS20 (5 mg ml−1), measurements were conducted at a sample-to-detector distance of 600 mm using C1 collimation (1.3 × 1.3 mm beam) with an exposure time of 60 s per frame for 60 frames (1 h total). Additional measurements were performed at 1800 mm using C4 collimation (0.5 × 0.5 mm beam) with an exposure time of 600 s per frame for 60 frames (10 h total).
For VitE-TPGS (0.4 wt%), measurements at 600 mm with C1 collimation (1.3 × 1.3 mm beam) were performed for 60 s per frame for 60 frames (1 h total). Measurements at 1800 mm with C4 collimation (0.5 × 0.5 mm beam) were performed for 300 s per frame for 60 frames (5 h total).
3.2. Synchrotron SAXS data
SAXS data were collected at the P12 BioSAXS beamline of the PETRA III synchrotron on the DESY campus in Hamburg, Germany. The synchrotron source employs the PETRA III U29 undulator with a gap of 20.95 mm, corresponding to an energy of 10 keV and a wavelength λ = 1.24 Å. Scattering intensities were recorded using a Pilatus 6M hybrid photon detector.
For PS20, measurements were performed at a sample-to-detector distance of 3 m, covering a q range of 0.03–4.43 nm−1, with a beam size at the sample of 250 × 100 µm. Due to the low sample viscosity, an autosampler was used and measurements were carried out in batch mode. For each sample, 40 µl was loaded into a flow-through quartz capillary and measured at 20°C, using an exposure time of 0.1 s per frame for 40 frames (total exposure of 4 s). For VitE-TPGS (0.4 wt%), samples were measured only with the laboratory SAXS instrument.
3.3. Data reduction
To extend the accessible q range, particularly at low scattering vector values in the SAXS regime, laboratory measurements collected at two different sample–detector distances (600 and 1800 mm) were merged.
To this end, the multiple acquisitions collected at each distance were first averaged to improve the signal-to-noise ratio. The averaging procedure was carried out using PRIMUS (ATSAS; Franke et al., 2025
; Manalastas-Cantos et al., 2021
), which performs a preliminary correlation-based comparison of all profiles to identify and exclude datasets affected by spurious signals or acquisition artifacts.
After averaging, the datasets acquired at 600 and 1800 mm were manually merged, with the scattering curves automatically rescaled by the XSACT software (https://www.xenocs.com/), so that they were quantitatively consistent across distances. A q range with a reasonably good signal-to-noise ratio in the low-q region – corresponding to the dataset collected at the longer distance – was selected and combined with the higher-q portion of the dataset collected at the shorter distance. This approach aims to produce the smoothest and most continuous profile possible across the full q range.
Raw SAXS detector frames collected at the synchrotron facility were reduced using the standard beamline pipeline based on pyFAI (https://github.com/silx-kit/pyFAI). 2D detector images were corrected for dark current, environment, detector geometry and masked pixels, and the q scale was calibrated using a silver behenate standard. Corrected 2D patterns were azimuthally integrated to obtain 1D intensity profiles with polarization and solid-angle normalization. Intensities were normalized by incident beam and sample transmission, and solvent/background scattering measured under identical conditions was subtracted to isolate the sample signal. Data were placed on an absolute intensity scale using a calibrated reference standard (typically water). The final reduced 1D scattering curves were used for structural analysis.
3.4. Model-dependent fits of SAXS data
The detailed results obtained from the model-dependent fits of SAXS data have been reported in Section S4 of the supporting information. Here, we summarize the main findings.
Fig. 3
compares the laboratory (red) and synchrotron (black) SAXS data for PS20 micelles (monomer concentration of 5 mg ml−1). The two datasets differ only in their noise level. The NS channels (blue lines), derived as previously described, are also shown. The results of the model-dependent fits reported in the supporting information, Section S4, show that also laboratory data are sufficiently robust to characterize the structure of core–shell spheroidal micelles. Despite the higher level of noise compared with synchrotron radiation, the laboratory measurements yield structural parameters consistent with those obtained at the synchrotron.
| Figure 3 Experimental SAXS profiles for PS20 micelles (monomer concentration of 5 mg ml−1) measured in the laboratory (red) and at the synchrotron (black). The blue lines are the NS channels. |
This is primarily due to the high oversampling ratio of the SAXS data, which provides enough mathematical redundancy to compensate for a lower signal-to-noise ratio. By applying appropriate smoothing techniques – leveraging the fact that the signal is sampled well beyond the Nyquist sampling rate – the underlying scattering profiles become nearly identical, allowing for a reliable retrieval of the morphology of nano-objects without significant loss of accuracy.
Fig. 4
shows the experimental laboratory SAXS data (black) compared with the fit (red) given by equations (S1)–(S21). Details of the fit have been reported in Section S4.
| Figure 4 Experimental SAXS profiles for PS20 micelles (monomer concentration of 5 mg ml−1) measured in the laboratory (black) versus the profile calculated minimizing the χ2 (red), as described in detail in Section S4 of the supporting information. Vertical lines are the NS channels. |
Fig. 5
shows the experimental laboratory SAXS data (black) compared with the model-dependent fit (red) obtained for VitE-TPGS micelles. Details of the fit have been reported in Section S4.
| Figure 5 Experimental SAXS profiles for VitE-TPGS micelles measured in the laboratory (black) versus the profile obtained by the model-dependent fit (red), as described in detail in Section S4 of the supporting information. Vertical lines are the NS channels. |
For a VitE-TPGS solution of the same concentration and for laboratory SAXS data, De Caro et al. (2025
) found Dmax = 122.0 ± 0.1 Å, = 16.7 ± 0.1 Å, X = 0.137 ± 0.002, ɛ = 1.51 ± 0.01, ɛM = 1.32 ± 0.02, Xρ = −1.81 ± 0.05, Rg = 47.3 ± 0.5 Å and
= 0.056 ± 0.001 cm−1, in agreement with the results obtained here from the laboratory data (see Tables 2–4 reported in the previous sections, and Table S3 of Section S4). Therefore, we can conclude that the SAXS profile determined in the laboratory contains all the scattering information needed for a fast and reliable structural solution of the micelles under study by means of the new proposed approach, under the assumptions of the simplified model discussed in the supporting information.
Finally, in Fig. 6
we show the results of the model-dependent fit (red) for the SASDFN8 apoferritin SAXS data (black) derived from horse spleen, taken from the SASBDB repository. From Table 4
, previously reported, and Table S3 of the supporting information, it is evident that there is overall good agreement between the results obtained using our approach and those reported in the SASBDB repository for the SASDFN8 apoferritin, particularly in the `.out file' generated by GNOM (Svergun, 1992
).
| Figure 6 Experimental SAXS profiles for apoferritin (SASDFN8) (black) versus the profile (red) calculated by means of the model-dependent fit, described in Section S4 of the supporting information. Vertical lines are the NS channels. |
4. Related literature
The following references are cited in the supporting information: Doucet et al. (2025
), Lipfert et al. (2007
).
5. Conclusions
In this work, we have presented a method to reliably estimate Dmax directly from the SAXS profile, using a new approach based on the Shannon sampling theorem. The proposed new approach enables the fast determination of several structural quantities characterizing the unknown structure of proteins as well as core–shell spheroidal micelles – such as mass, volume, ellipticity, maximum size, number of monomers, hydrated fraction and shell thickness – without the need to calculate the pair distribution function . The application of the same formalism to relevant known proteins, taken from the SASBDB repository, has shown the reliability of the new approach in deriving the maximum size and mass of proteins. A new formula has been derived for determining the aggregation number of micelles starting from SAXS data on a relative scale. All these quantities can be derived from experimental SAXS data, through the Shannon sampling formalism. For SAXS data on an absolute scale also the shell and core electron-density difference can be derived without any model-dependent fit of SAXS data.
We have further demonstrated that the new proposed approach, even using laboratory SAXS data, can yield reliable morphological results, as validated using experimental measurements on PS20 and VitE-TPGS micelles, despite the higher noise level compared with synchrotron data. A good overall agreement has been obtained for the structural parameters derived for core–shell structures by the new approach both with the model-dependent fit of experimental SAXS data and with the results reported in the literature, under the assumptions of the simplified model discussed in the supporting information.
The approach discussed here, based on the Nyquist–Shannon sampling of the experimental SAXS profile, can be applied to nano-objects of any shape. The new proposed approach involves calculating integral quantities of the SAXS data, through the Nyquist–Shannon sampling formula. Consequently, it does not suffer from limiting effects due to the q-truncation of experimental data, as happens in the Fourier transform of SAXS data used for calculating the pair distance distribution function. Our tests have shown that the new proposed approach is able to give correct structural information about the scattering particles, without the need to calculate the pair distribution function or model-dependent fits of SAXS data. Therefore, this approach is quite useful for high-throughput data analysis, typical of massive SAXS collection (in situ and operando experiments). For core–shell micelles, the structural results obtained by the new approach could also be considered a useful preliminary and quantitative analysis of the structural features that could be refined/checked by a model-dependent fit of the experimental SAXS data.
More articulated modeling of micelle SAXS data, which could take into account complexity related to excluded volume and fluctuation effects in the corona (Svaneborg & Pedersen, 2000
) or coupling of size dispersity and eccentricity (Lebecque et al., 2017
), could benefit from the preliminary quantitative results extracted by the Nyquist–Shannon-based analysis proposed here.
In particular, modeling of surfactant and polymeric micelles builds upon earlier work originally developed for the analysis of relatively noisy SANS data, and more recent approaches continue to evolve, improving information extraction even from sparse SANS datasets (Tung et al., 2025
). Nevertheless, a fast evaluation of some structural properties of micelles and proteins (maximum size, aggregation number of monomers, mass, volume), working directly on SAXS laboratory experimental data without any need of conversion of data in real space, is a useful tool. It can be combined with different techniques (such as SANS) when further characterization is required to define in more detail the structure of the scattering objects.
The new approach will be implemented in the software SUNBIM (Scattarella et al., 2025
). The future use of this software by users will help us to further understand the actual reliability of structural information obtained by this approach in (i) retrieving the structure of core–shell micelles without any model-dependent fitting of the SAXS data and (ii) determining useful preliminary structural information on proteins from SAXS data, such as molecular mass.
Supporting information
Supporting information. DOI: https://doi.org/10.1107/S1600576726004358/fc5088sup1.pdf
Acknowledgements
Open access publishing facilitated by Consiglio Nazionale delle Ricerche, as part of the Wiley–CRUI–CARE agreement.
Conflict of interest
There are no conflicts of interest to declare.
Data availability
Data are available from the SASBDB repository: https://www.sasbdb.org.
Funding information
The authors acknowledge support from the project Potentiating the Italian Capacity for Structural Biology Services in Instruct-ERIC (ITACA.SB) (project No. IR0000009, CUP B53C22001790006), funded by the European Union's NextGenerationEU under the MUR call 3264/2021 PNRR M4/C2/L3.1.1. This work benefitted from the use of the SasView application, originally developed under NSF award DMR-0520547. SasView contains code developed with funding from the European Union's Horizon 2020 research and innovation programme under the SINE2020 project, grant agreement No. 654000.
References
Breßler, I., Kohlbrecher, J. & Thünemann, A. F. (2015). J. Appl. Cryst. 48, 1587–1598. Web of Science CrossRef IUCr Journals Google Scholar
Castelletto, V. & Hamley, I. W. (2002). Curr. Opin. Colloid Interface Sci. 7, 167–172. CrossRef CAS Google Scholar
Cezar, H. M., Bjørnestad, V. A., Prévost, S., Lund, R. & Cascella, M. (2025). Small Struct. 6, 2400553. CrossRef Google Scholar
Chen, R., Song, Y., Wang, Z., Ji, H., Du, Z., Ma, Q., Yang, Y., Liu, X., Li, N. & Sun, Y. (2023). Int. J. Biol. Macromol. 251, 126288. CrossRef PubMed Google Scholar
De Caro, L., Stoll, T., Grandeury, A., Gozzo, F. & Giannini, C. (2024). Pharmaceutics 16, 604. CrossRef PubMed Google Scholar
De Caro, L., Stoll, T., Grandeury, A., Gozzo, F. & Giannini, C. (2025). Pharmaceutics 17, 431. CrossRef PubMed Google Scholar
Doucet, M., Adams, M., Agouzal, N., Alina, G., Anuchitanukul, A., Attala, Z., Bakker, J., Beaucage, P., Berger, J., Bourne, R., Bouwman, W., Bressler, I., Butler, P., Cadwallader-Jones, I., Campbell, K., Cho, J.-H., Cooper-Benun, T., Cortes Hernandez, R., Doutch, J., et al. (2025). SasView. Version 6.0.1. https://doi.org/10.5281/zenodo.14940319. Google Scholar
Ehrit, G., Gräwert, T. W., Göddeke, H., Konarev, P. V., Svergun, D. I. & Nagel, N. (2023). Biophys. J. 122, 3078–3088. CrossRef CAS PubMed Google Scholar
Evers, D., Giebel, J., Nienau, F., Carle, S., Gorissen, S., Buske, J., Herbig, M. E., Garidel, P. & Hagelskamp, E. (2025). Eur. J. Pharm. Biopharm. 216, 114854. CrossRef PubMed Google Scholar
Fischer, H., de Oliveira Neto, M., Napolitano, H. B., Polikarpov, I. & Craievich, A. F. (2010). J. Appl. Cryst. 43, 101–109. Web of Science CrossRef CAS IUCr Journals Google Scholar
Franke, D., Gräwert, T. & Svergun, D. I. (2025). J. Appl. Cryst. 58, 1027–1033. Web of Science CrossRef CAS IUCr Journals Google Scholar
Glatter, O. & Kratky, O. (1982). Small-Angle X-ray Scattering. Academic Press. Google Scholar
Grant, T. D. (2022). J. Appl. Cryst. 55, 1116–1124. Web of Science CrossRef CAS IUCr Journals Google Scholar
Guinier, A. & Fournet, G. (1955). Small-Angle Scattering of X-rays. Wiley. Google Scholar
Hansen, S. (2000). J. Appl. Cryst. 33, 1415–1421. Web of Science CrossRef CAS IUCr Journals Google Scholar
König, N., Willner, L., Carlström, G., Zinn, T., Knudsen, K. D., Rise, F., Topgaard, D. & Lund, R. (2020). Macromolecules 53, 10686–10698. PubMed Google Scholar
Lebecque, S., Crowet, J. M., Nasir, M. N., Deleu, M. & Lins, L. (2017). J. Mol. Graphics Modell. 72, 6–15. CrossRef CAS Google Scholar
le Maire, M., Champeil, P. & Møller, J. V. (2000). Biochim. Biophys. Acta 1508, 86–111. CrossRef PubMed CAS Google Scholar
Lipfert, J., Columbus, L., Chu, V. B., Lesley, S. A. & Doniach, S. (2007). J. Phys. Chem. B 111, 12427–12438. CrossRef PubMed CAS Google Scholar
Manalastas-Cantos, K., Konarev, P. V., Hajizadeh, N. R., Kikhney, A. G., Petoukhov, M. V., Molodenskiy, D. S., Panjkovich, A., Mertens, H. D. T., Gruzinov, A., Borges, C., Jeffries, C. M., Svergun, D. I. & Franke, D. (2021). J. Appl. Cryst. 54, 343–355. Web of Science CrossRef CAS IUCr Journals Google Scholar
Moore, P. B. (1980). J. Appl. Cryst. 13, 168–175. CrossRef CAS IUCr Journals Web of Science Google Scholar
Orthaber, D., Bergmann, A. & Glatter, O. (2000). J. Appl. Cryst. 33, 218–225. Web of Science CrossRef CAS IUCr Journals Google Scholar
Patriati, A., Suparno, N., Sulungbudi, G. T., Mujamilah, M. & Putra, E. G. R. (2020). Indones. J. Chem. 20, 1178–1183. CrossRef CAS Google Scholar
Pedersen, J. S. (1997). Adv. Colloid Interface Sci. 70, 171–210. CrossRef CAS Web of Science Google Scholar
Rambo, R. P. & Tainer, J. A. (2013). Nature 496, 477–481. CrossRef CAS PubMed Google Scholar
Roe, R.-J. (2000). Methods of X-ray and Neutron Scattering in Polymer Science. Oxford University Press. Google Scholar
Scattarella, F., Altamura, D., Sibillano, T., De Caro, L., Siliqi, D. & Giannini, C. (2025). J. Appl. Cryst. 58, 1817–1826. CrossRef CAS IUCr Journals Google Scholar
Shannon, C. E. (1948). Bell Syst. Tech. J. 27, 379–423. CrossRef Web of Science Google Scholar
Stuhrmann, H. B. (2008). Acta Cryst. A64, 181–191. Web of Science CrossRef CAS IUCr Journals Google Scholar
Svaneborg, C. & Pedersen, J. S. (2000). J. Chem. Phys. 112, 9661–9670. CrossRef CAS Google Scholar
Svergun, D. I. (1992). J. Appl. Cryst. 25, 495–503. CrossRef CAS Web of Science IUCr Journals Google Scholar
Tung, C. H., Yip, S., Huang, G. R., Porcar, L., Shinohara, Y., Sumpter, B. G., Ding, L., Do, C. & Chen, W. R. (2025). J. Colloid Interface Sci. 692, 137554. CrossRef PubMed Google Scholar
Valentini, E., Kikhney, A. G., Previtali, G., Jeffries, C. M. & Svergun, D. I. (2015). Nucleic Acids Res. 43, D357–D363. Web of Science CrossRef CAS PubMed Google Scholar
Vass, S., Pedersen, J. S., Pleštil, J., Laggner, P., Rétfalvi, E., Varga, I. & Gilányi, T. (2008). Langmuir 24, 408–417. CrossRef PubMed CAS Google Scholar
Vestergaard, B. & Hansen, S. (2006). J. Appl. Cryst. 39, 797–804. Web of Science CrossRef CAS IUCr Journals Google Scholar
Zeng, J., Bian, F., Wang, J., Li, X., Wang, Y., Tian, F. & Zhou, P. (2017). J. Synchrotron Rad. 24, 509–520. Web of Science CrossRef CAS IUCr Journals Google Scholar
This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.

journal menu
access



