Hydrodynamic Modeling and Its Application in AUC

Strong non-ideality effects at low protein concentrations: considerations for elongated proteins

A recent investigation was aimed at obtaining structural information on a highly extended protein via SEC-MALS-SAXS. Significantly broadened elution peaks were observed, reminiscent of a phenomenon known as viscous fingering. This phenomenon is usually observed above 50 mg/mL for proteins like bovine serum albumin (BSA). Interestingly, the highly extended protein (Brpt5.5) showed viscous fingering at concentrations lower than 5 mg/mL. The current study explores this and other non-ideal behavior, emphasizing the presence of these effects at relatively low concentrations for extended proteins. BSA, Brpt5.5, and a truncated form of Brpt5.5 referred to as Brpt1.5 are studied systematically using size-exclusion chromatography (SEC), sedimentation velocity analytical ultracentrifugation (AUC), and viscosity. The viscous fingering effect is quantified using two approaches and is found to correlate well with the intrinsic viscosity of the proteins-Brpt5.5 exhibits the most severe effect and is the most extended protein tested in the study. By AUC, the hydrodynamic non-ideality was measured for each protein via global analysis of a concentration series. Compared to BSA, both Brpt1.5 and Brpt5.5 showed significant non-ideality that could be easily visualized at concentrations at or below 5 mg/mL and 1 mg/mL, respectively. A variety of relationships were examined for their ability to differentiate the proteins by shape using information from AUC and/or viscosity. Furthermore, these relationships were also tested in the context of hydrodynamic modeling. The importance of considering non-ideality when investigating the structure of extended macromolecules is discussed.

Revisiting macromolecular hydration with HullRadSAS

Hydration of biological macromolecules is important for their stability and function. Historically, attempts have been made to describe the degree of macromolecular hydration using a single parameter over a narrow range of values. Here, we describe a method to calculate two types of hydration: surface shell water and entrained water. A consideration of these two types of hydration helps to explain the "hydration problem" in hydrodynamics. The combination of these two types of hydration allows accurate calculation of hydrodynamic volume and related macromolecular properties such as sedimentation and diffusion coefficients, intrinsic viscosities, and the concentration-dependent non-ideality identified with sedimentation velocity experiments.

A database of calculated solution parameters for the AlphaFold predicted protein structures

Recent spectacular advances by AI programs in 3D structure predictions from protein sequences have revolutionized the field in terms of accuracy and speed. The resulting “folding frenzy” has already produced predicted protein structure databases for the entire human and other organisms’ proteomes. However, rapidly ascertaining a predicted structure’s reliability based on measured properties in solution should be considered. Shape-sensitive hydrodynamic parameters such as the diffusion and sedimentation coefficients (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D_{t(20,w)}^{0}}$$\end{document}Dt(20,w)0, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${s_{{\left( {{20},w} \right)}}^{{0}} }$$\end{document}s20,w0) and the intrinsic viscosity ([η]) can provide a rapid assessment of the overall structure likeliness, and SAXS would yield the structure-related pair-wise distance distribution function p(r) vs. r. Using the extensively validated UltraScan SOlution MOdeler (US-SOMO) suite, a database was implemented calculating from AlphaFold structures the corresponding \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D_{t(20,w)}^{0}}$$\end{document}Dt(20,w)0, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${s_{{\left( {{20},w} \right)}}^{{0}} }$$\end{document}s20,w0, [η], p(r) vs. r, and other parameters. Circular dichroism spectra were computed using the SESCA program. Some of AlphaFold’s drawbacks were mitigated, such as generating whenever possible a protein’s mature form. Others, like the AlphaFold direct applicability to single-chain structures only, the absence of prosthetic groups, or flexibility issues, are discussed. Overall, this implementation of the US-SOMO-AF database should already aid in rapidly evaluating the consistency in solution of a relevant portion of AlphaFold predicted protein structures.

Recent advances in the UltraScan SOlution MOdeller (US-SOMO) hydrodynamic and small-angle scattering data analysis and simulation suite

The UltraScan SOlution MOdeller (US-SOMO) is a comprehensive, public domain, open-source suite of computer programs centred on hydrodynamic modelling and small-angle scattering (SAS) data analysis and simulation. We describe here the advances that have been implemented since its last official release (#3087, 2017), which are available from release #3141 for Windows, Linux and Mac operating systems. A major effort has been the transition from the legacy Qt3 cross platform software development and user interface library to the modern Qt5 release. Apart from improved graphical support, this has allowed the direct implementation of the newest, almost two-orders of magnitude faster version of the ZENO hydrodynamic computation algorithm for all operating systems. Coupled with the SoMo-generated bead models with overlaps, ZENO provides the most accurate translational friction computations from atomic-level structures available (Rocco and Byron Eur Biophys J 44:417-431, 2015a), with computational times comparable with or faster than those of other methods. In addition, it has allowed us to introduce the direct representation of each atom in a structure as a (hydrated) bead, opening interesting new modelling possibilities. In the small-angle scattering (SAS) part of the suite, an indirect Fourier transform Bayesian algorithm has been implemented for the computation of the pairwise distance distribution function from SAS data. Finally, the SAS HPLC module, recently upgraded with improved baseline correction and Gaussian decomposition of not baseline-resolved peaks and with advanced statistical evaluation tools (Brookes et al. J Appl Cryst 49:1827-1841, 2016), now allows automatic top-peak frame selection and averaging.

HullRad: Fast Calculations of Folded and Disordered Protein and Nucleic Acid Hydrodynamic Properties

Hydrodynamic properties are useful parameters for estimating the size and shape of proteins and nucleic acids in solution. The calculation of such properties from structural models informs on the solution properties of these molecules and complements corresponding structural studies. Here we report, to our knowledge, a new method to accurately predict the hydrodynamic properties of molecular structures. This method uses a convex hull model to estimate the hydrodynamic volume of the molecule and is orders of magnitude faster than common methods. It works well for both folded proteins and ensembles of conformationally heterogeneous proteins and for nucleic acids. Because of its simplicity and speed, the method should be useful for the modification of computer-generated, intrinsically disordered protein ensembles and ensembles of flexible, but folded, molecules in which rapid calculation of experimental parameters is needed. The convex hull method is implemented in a Python script called HullRad. The use of the method is facilitated by a web server and the code is freely available for batch applications.