Efficient Analysis of Small-Angle Scattering Curves for Large Biomolecular Assemblies Using Monte Carlo Methods

Structure elicitation from small angle scattering (i.e. SAXS, SANS) curves of large biomolecular assemblies is notoriously challenging. This is because the simulation of high-resolution features in the structure of large macromolecular assemblies, such as de-novo protein assemblies, is computationally demanding when it needs to cover a broad range of length scales. Conventional methods such as the numerical approximation to the Debye equation or the use of spherical harmonics do not scale well as the size of the assembly increases, which limits their application to small structures (e.g. individual proteins). This work explores the effectiveness of a Monte Carlo method to simulate and fit scattering curves for large biomolecular assemblies spanning over ranges covering atomic and molecular detail, as well as large scale (100’s of nm) features. Due to its excellent speed and scalability, it can be combined with a fitting algorithm to extract structural features from experimental small-angle scattering curves in biomolecular assemblies that are otherwise intractable for interpretation. This work first demonstrates the effectiveness of the tool using experimental small angle X-ray scattering (SAXS) data from tile-like proteins that assemble into 1D tube-like macromolecular structures. The diameter distribution of tubes is extracted from SAXS fits, and this quantitatively compared with distributions from electron microscopy. SAXS data is also obtained from 2D sheet-like protein assemblies and the proposed method is used to quantify structural features such as the separation distance between protein building blocks and the flexing of the sheet. An open-source implementation of the methodology is provided for use in a broad range of biological systems involving multi-scale scattering analysis.