Computing translational diffusion and sedimentation coefficients: an evaluation of experimental data and programs

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.

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.

Sedimentation Velocity Methods for the Characterization of Protein Heterogeneity and Protein Affinity Interactions

Sedimentation velocity analytical ultracentrifugation is a powerful and versatile tool for the characterization of proteins and macromolecular complexes in solution. The direct modeling of the sedimentation process using modern computational strategies allows among others to assess the homogeneity/heterogeneity state of protein samples and to characterize protein associations. In this chapter, we will provide theoretical backgrounds and protocols to analyze the size distribution of protein samples and to determine the affinity of protein-protein hetero-associations.

Analytical ultracentrifugation (AUC): a seminal tool offering multiple solutions

Analytical ultracentrifugation (AUC) remains a highly versatile and widely applicable tool for the analysis of macromolecules and their interactions. The current state-of-the-art was demonstrated at a recent international meeting held in Glasgow, Scotland, in July 2017, the 23rd International Analytical Ultracentrifugation Workshop and Symposium. This special issue showcases the reports made at the meeting, which concerned the application of AUC to a wide range of topics in biochemical and polymer science including antibody and membrane protein characterisation, and protein-carbohydrate interactions. Presentations on development and testing of new instrumentation and methods of analysis were a particular feature of the meeting, including the optimisation of experimental protocols, and the latest optimised computational approaches to experimental simulation and the modelling of macromolecular structures.

Structural and functional characterisation of the entry point to pyocyanin biosynthesis in Pseudomonas aeruginosa defines a new 3-deoxy-d-arabino-heptulosonate 7-phosphate synthase subclass

In Pseudomonas aeruginosa (Pae), the shikimate pathway end product, chorismate, serves as the last common precursor for the biosynthesis of both primary aromatic metabolites, including phenylalanine, tyrosine and tryptophan, and secondary aromatic metabolites, including phenazine-1-carboxylic acid (PCA) and pyocyanin (PYO). The enzyme 3-deoxy-d-arabino-heptulosonate 7-phosphate synthase (DAH7PS) catalyses the first committed step of the shikimate pathway, en route to chorismate. P. aeruginosa expresses multiple, distinct DAH7PSs that are associated with either primary or secondary aromatic compound biosynthesis. Here we report the structure of a type II DAH7PS, encoded by phzC as part of the duplicated phenazine biosynthetic cluster, from P. aeruginosa (PAO1) revealing for the first time the structure of a type II DAH7PS involved in secondary metabolism. The omission of the structural elements α_2a and α_2b, relative to other characterised type II DAH7PSs, leads to the formation of an alternative, dimeric, solution-state structure for this type II DAH7PS with an oligomeric interface that has not previously been characterised and that does not facilitate the formation of aromatic amino acid allosteric binding sites. The sequence similarity and, in particular, the common N-terminal extension suggest a common origin for the type II DAH7PSs from P. aeruginosa. The results described in the present study support an expanded classification of the type II DAH7PSs as type II_A and type II_B based on sequence characteristics, structure and function of the resultant proteins, and on defined physiological roles within primary or secondary metabolism.

GRPY: An Accurate Bead Method for Calculation of Hydrodynamic Properties of Rigid Biomacromolecules

Two main problems that arise in the context of hydrodynamic bead modeling are an inaccurate treatment of bead overlaps and the necessity of using volume corrections when calculating intrinsic viscosity. We present a formalism based on the generalized Rotne-Prager-Yamakawa approximation that successfully addresses both of these issues. The generalized Rotne-Prager-Yamakawa method is shown to be highly effective for the calculation of transport properties of rigid biomolecules represented as assemblies of spherical beads of different sizes, both overlapping and nonoverlapping. We test the method on simple molecular shapes as well as real protein structures and compare its performance with other computational approaches.

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.

Physicochemical improvement of rabbit derived single-domain antibodies by substitutions with amino acids conserved in camelid antibodies

Recently, we showed that immunized rabbit heavy chain variable regions (rVHs) can have strong antigen binding activity comparable to that of the camelid variable domain of the heavy chain of heavy chain antibody (VHH). These rVHs lack the light chain variable regions (rVLs), which exist in the authentic Fab format; thus, molecular surfaces at the interface region of rVHs are exposed to solvent. This physical feature may change physicochemical properties, such as causing reduced stability. By overcoming potential physicochemical issues through engineering the interface region, rVHs could become more useful as single-domain antibodies. In this study, we substituted amino acid residues conserved at the interface region of rVHs with those of VHHs. These substitutions included V37F, involving substitution of a residue in the hydrophobic core with a bulkier hydrophobic amino acid, and G44E/L45R, involving double substitutions of highly exposed residues with more hydrophilic ones. As expected, biophysical and structural characterizations showed that the V37F substitution markedly enhanced the thermal stability through increased hydrophobic packing, while G44E/L45R substitutions greatly reduced hydrophobicity of the interface. The quadruple substitutions of V37F/G44E/L45R/F91Y resulted in not only enhancements of thermal stability and reduction in hydrophobicity, both in an additive manner, but also synergistic improvement of purification yield. This quadruple mutant exhibited greatly reduced non-specific binding with improved colloidal stability owing to the reduced hydrophobicity. The approach used in this study should further enhance the utility of rVHs and promote research and development of single-domain antibodies.

Insight into adaptive remodeling of the rotor ring complex of the bacterial flagellar motor

The bacterial flagellar motor rotates in both counterclockwise (CCW) and clockwise (CW) directions. FliG, FliM and FliN form the C ring on the cytoplasmic face of the MS ring made of a transmembrane protein, FliF. The C ring acts not only as a rotor but also as a switch of the direction of motor rotation. FliG consists of three domains: FliG_N, FliG_M and FliG_C. FliG_N directly binds to FliF. Intermolecular interactions between FliG_M and FliG_C drive FliG ring formation. FliG_M is responsible for the interaction with FliM. FliG_C is involved in the interaction with the stator protein MotA. Adaptive remodeling of the C ring occurs when the motor switches between the CCW and CW states. However, it remained unknown how. Here, we report the effects of a CW-locked deletion mutation (ΔPEV) in FliG of Thermotaoga maritia (Tm-FliG) on FliG-FliG and FliG-FliM interactions. The PEV deletion stabilized the intramolecular interaction between FliG_M and FliG_C, thereby suppressing the oligomerization of Tm-FliG_MC in solution. This deletion also induced a conformational change of Helix_MC connecting FliG_M and FliG_C to reduce the binding affinity of Tm-FliG_MC for FliM. We will discuss adaptive remodeling of the C ring responsible for flagellar motor switching.

Disassembly of the self-assembled, double-ring structure of proteasome α7 homo-tetradecamer by α6

The 20S core particle of the eukaryotic proteasome is composed of two α- and two β-rings, each of which is a hetero-heptamer composed of seven homologous but distinct subunits. Although formation of the eukaryotic proteasome is a highly ordered process assisted by assembly chaperones, α7, an α-ring component, has the unique property of self-assembling into a homo-tetradecamer. We used biophysical methods to characterize the oligomeric states of this proteasome subunit and its interaction with α6, which makes direct contacts with α7 in the proteasome α-ring. We determined a crystal structure of the α7 tetradecamer, which has a double-ring structure. Sedimentation velocity analytical ultracentrifugation and mass spectrometric analysis under non-denaturing conditions revealed that α7 exclusively exists as homo-tetradecamer in solution and that its double-ring structure is disassembled upon the addition of α6, resulting in a 1:7 hetero-octameric α6–α7 complex. Our findings suggest that proteasome formation involves the disassembly of non-native oligomers, which are assembly intermediates.

Protein-protein interactions: a supra-structural phenomenon demanding trans-disciplinary biophysical approaches

Responsive formation of protein:protein interaction (PPI) upon diverse stimuli is a fundament of cellular function. As a consequence, PPIs are complex, adaptive entities, and exist in structurally heterogeneous interplays defined by the energetic states of the free and complexed protomers. The biophysical and structural investigations of PPIs consequently demand hybrid approaches, implementing orthogonal methods and strategies for global data analysis. Currently, impressive developments in hardware and software within several methodologies define a new era for the biostructural community. Data can be obtained at increasing resolution, at relevant time-scales and under increasingly relevant experimental conditions, intricate data are interpreted reliably, and the questions posed and answered grow in complexity. With this review, highlights from the study of PPIs using a multitude of biophysical methods, are reported. The aim is to depict how the elucidation of the interplay of structures requires the interplay of methods.