Testing mixing rules for structural and dynamical quantities in multi-component crowded protein solutions.

Crowding effects significantly influence the phase behavior and the structural and dynamic properties of the concentrated protein mixtures present in the cytoplasm of cells or in the blood serum. This poses enormous difficulties for our theoretical understanding and our ability to predict the behavior of these systems. While the use of course grained colloid-inspired models allows us to reproduce the key physical solution properties of concentrated monodisperse solutions of individual proteins, we lack corresponding theories for complex polydisperse mixtures. Here, we test the applicability of simple mixing rules in order to predict solution properties of protein mixtures. We use binary mixtures of the well-characterized bovine eye lens proteins α and γB crystallin as model systems. Combining microrheology with static and dynamic scattering techniques and observations of the phase diagram for liquid-liquid phase separation, we show that reasonably accurate descriptions are possible for macroscopic and mesoscopic signatures, while information on the length scale of the individual protein size requires more information on cross-component interaction.

Lipid exchange of apolipoprotein A-I amyloidogenic variants in reconstituted high-density lipoprotein with artificial membranes.

High-density lipoproteins (HDLs) are responsible for removing cholesterol from arterial walls, through a process known as reverse cholesterol transport. The main protein in HDL, apolipoprotein A-I (ApoA-I), is essential to this process, and changes in its sequence significantly alter HDL structure and functions. ApoA-I amyloidogenic variants, associated with a particular hereditary degenerative disease, are particularly effective at facilitating cholesterol removal, thus protecting carriers from cardiovascular disease. Thus, it is conceivable that reconstituted HDL (rHDL) formulations containing ApoA-I proteins with functional/structural features similar to those of amyloidogenic variants hold potential as a promising therapeutic approach. Here we explored the effect of protein cargo and lipid composition on the function of rHDL containing one of the ApoA-I amyloidogenic variants G26R or L174S by Fourier transformed infrared spectroscopy and neutron reflectometry. Moreover, small-angle x-ray scattering uncovered the structural and functional differences between rHDL particles, which could help to comprehend higher cholesterol efflux activity and apparent lower phospholipid (PL) affinity. Our findings indicate distinct trends in lipid exchange (removal vs. deposition) capacities of various rHDL particles, with the rHDL containing the ApoA-I amyloidogenic variants showing a markedly lower ability to remove lipids from artificial membranes compared to the rHDL containing the native protein. This effect strongly depends on the level of PL unsaturation and on the particles' ultrastructure. The study highlights the importance of the protein cargo, along with lipid composition, in shaping rHDL structure, contributing to our understanding of lipid-protein interactions and their behavior.

Structural characterisation of α-synuclein-membrane interactions and the resulting aggregation using small angle scattering.

The presence of amyloid fibrils is a hallmark of several neurodegenerative diseases. Some amyloidogenic proteins, such as α-synuclein and amyloid ß, interact with lipids, and this interaction can strongly favour the formation of amyloid fibrils. In particular the primary nucleation step, i.e. the de novo formation of amyloid fibrils, has been shown to be accelerated by lipids. However, the exact mechanism of this acceleration is still mostly unclear. Here we use a range of scattering methods, such as dynamic light scattering (DLS) and small angle X-ray and neutron scattering (SAXS and SANS) to obtain structural information on the binding of α-synuclein to model membranes formed from negatively charged lipids and their co-assembly into amyloid fibrils. We find that the model membranes take an active role in the reaction. The binding of α synuclein to the model membranes immediately induces a major structural change in the lipid assembly, which leads to a break-up into small and mostly disc- or rod-like lipid-protein particles. This transition can be reversed by temperature changes or proteolytic protein removal. Incubation of the small lipid-α-synuclein particles for several hours, however, leads to amyloid fibril formation, whereby the lipids are incorporated into the amyloid fibrils.

The Onset of Molecule-Spanning Dynamics in Heat Shock Protein Hsp90.

Protein dynamics have been investigated on a wide range of time scales. Nano- and picosecond dynamics have been assigned to local fluctuations, while slower dynamics have been attributed to larger conformational changes. However, it is largely unknown how fast (local) fluctuations can lead to slow global (allosteric) changes. Here, fast molecule-spanning dynamics on the 100 to 200 ns time scale in the heat shock protein 90 (Hsp90) are shown. Global real-space movements are assigned to dynamic modes on this time scale, which is possible by a combination of single-molecule fluorescence, quasi-elastic neutron scattering and all-atom molecular dynamics (MD) simulations. The time scale of these dynamic modes depends on the conformational state of the Hsp90 dimer. In addition, the dynamic modes are affected to various degrees by Sba1, a co-chaperone of Hsp90, depending on the location within Hsp90, which is in very good agreement with MD simulations. Altogether, this data is best described by fast molecule-spanning dynamics, which precede larger conformational changes in Hsp90 and might be the molecular basis for allostery. This integrative approach provides comprehensive insights into molecule-spanning dynamics on the nanosecond time scale for a multi-domain protein.

High-Density Lipoprotein function is modulated by the SARS-CoV-2 spike protein in a lipid-type dependent manner.

There is a close relationship between the SARS-CoV-2 virus and lipoproteins, in particular high-density lipoprotein (HDL). The severity of the coronavirus disease 2019 (COVID-19) is inversely correlated with HDL plasma levels. It is known that the SARS-CoV-2 spike (S) protein binds the HDL particle, probably depleting it of lipids and altering HDL function. Based on neutron reflectometry (NR) and the ability of HDL to efflux cholesterol from macrophages, we confirm these observations and further identify the preference of the S protein for specific lipids and the consequent effects on HDL function on lipid exchange ability. Moreover, the effect of the S protein on HDL function differs depending on the individuals lipid serum profile. Contrasting trends were observed for individuals presenting low triglycerides/high cholesterol serum levels (LTHC) compared to high triglycerides/high cholesterol (HTHC) or low triglycerides/low cholesterol serum levels (LTLC). Collectively, these results suggest that the S protein interacts with the HDL particle and, depending on the lipid profile of the infected individual, it impairs its function during COVID-19 infection, causing an imbalance in lipid metabolism.

Effects of flexibility in coarse-grained models for bovine serum albumin and immunoglobulin G.

We construct a coarse-grained, structure-based, low-resolution, 6-bead flexible model of bovine serum albumin (BSA, PDB: 4F5S), which is a popular example of a globular protein in biophysical research. The model is obtained via direct Boltzmann inversion using all-atom simulations of a single molecule, and its particular form is selected from a large pool of 6-bead coarse-grained models using two suitable metrics that quantify the agreement in the distribution of collective coordinates between all-atom and coarse-grained Brownian dynamics simulations of solutions in the dilute limit. For immunoglobulin G (IgG), a similar structure-based 12-bead model has been introduced in the literature [Chaudhri et al., J. Phys. Chem. B 116, 8045 (2012)] and is employed here to compare findings for the compact BSA molecule and the more anisotropic IgG molecule. We define several modified coarse-grained models of BSA and IgG, which differ in their internal constraints and thus account for a variation of flexibility. We study denser solutions of the coarse-grained models with purely repulsive molecules (achievable by suitable salt conditions) and address the effect of packing and flexibility on dynamic and static behavior. Translational and rotational self-diffusivity is enhanced for more elastic models. Finally, we discuss a number of effective sphere sizes for the BSA molecule, which can be defined from its static and dynamic properties. Here, it is found that the effective sphere diameters lie between 4.9 and 6.1 nm, corresponding to a relative spread of about ±10% around a mean of 5.5 nm.

Structural and Dynamical Properties of Elastin-Like Peptides near Their Lower Critical Solution Temperature.

Elastin-like peptides (ELPs) are artificially derived intrinsically disordered proteins (IDPs) mimicking the hydrophobic repeat unit in the protein elastin. ELPs are characterized by a lower critical solution temperature (LCST) in aqueous media. Here, we investigate the sequence GVG(VPGVG)3 over a wide range of temperatures (below, around, and above the LCST) and peptide concentrations employing all-atom molecular dynamics simulations, where we focus on the role of intra- and interpeptide interactions. We begin by investigating the structural properties of a single peptide that demonstrates a hydrophobic collapse with temperature, albeit moderate, because the sequence length is short. We observe a change in the interaction between two peptides from repulsive to attractive with temperature by evaluating the potential of mean force, indicating an LCST-like behavior. Next, we explore dynamical and structural properties of peptides in multichain systems. We report the formation of dynamical aggregates with coil-like conformation, in which valine central residues play an important role. Moreover, the lifetime of contacts between chains strongly depends on the temperature and can be described by a power-law decay that is consistent with the LCST-like behavior. Finally, the peptide translational and internal motion are slowed by an increase in the peptide concentration and temperature.

Probing Cage Relaxation in Concentrated Protein Solutions by X-Ray Photon Correlation Spectroscopy.

Diffusion of proteins on length scales of their size is crucial for understanding the machinery of living cells. X-ray photon correlation spectroscopy (XPCS) is currently the only way to access long-time collective diffusion on these length scales, but radiation damage so far limits the use in biological systems. We apply a new approach to use XPCS to measure cage relaxation in crowded α-crystallin solutions. This allows us to correct for radiation effects, obtain missing information on long time diffusion, and support the fundamental analogy between protein and colloid dynamical arrest.

Expanding the Toolbox for Bicelle-Forming Surfactant-Lipid Mixtures.

Bicelles are disk-shaped models of cellular membranes used to study lipid-protein interactions, as well as for structural and functional studies on transmembrane proteins. One challenge for the incorporation of transmembrane proteins in bicelles is the limited range of detergent and lipid combinations available for the successful reconstitution of proteins in model membranes. This is important, as the function and stability of transmembrane proteins are very closely linked to the detergents used for their purification and to the lipids that the proteins are embedded in. Here, we expand the toolkit of lipid and detergent combinations that allow the formation of stable bicelles. We use a combination of dynamic light scattering, small-angle X-ray scattering and cryogenic electron microscopy to perform a systematic sample characterization, thus providing a set of conditions under which bicelles can be successfully formed.

Dynamics of post-occlusion water diffusion in stratum corneum.

Diffusion of water through membranes presents a considerable challenge, as the diffusivity often depends on the local concentration of water. One particular example with strong biological relevance is the stratum corneum (SC) as the primary permeability barrier for the skin. A simple alternative for the constant diffusivity model is provided by the Fujita's two-parameter rational approximation, which captures the experimentally observed fact that the SC diffusion constant for water increases with increasing the water concentration. Based on Fick's law of diffusion, a one-dimensional concentration-dependent diffusion model is developed and applied for the analysis of both the steady-state transepidermal water loss (TEWL) and the non-steady-state so-called skin surface water loss (SSWL) occurred after removal of an occlusion patch from the SC surface. It is shown that some of the age-related changes in the SSWL can be qualitatively explained by the variation of the dimensionless Fujita concentration-dependence parameter.

Short-Time Transport Properties of Bidisperse Suspensions of Immunoglobulins and Serum Albumins Consistent with a Colloid Physics Picture.

The crowded environment of biological systems such as the interior of living cells is occupied by macromolecules with a broad size distribution. This situation of polydispersity might influence the dependence of the diffusive dynamics of a given tracer macromolecule in a monodisperse solution on its hydrodynamic size and on the volume fraction. The resulting size dependence of diffusive transport crucially influences the function of a living cell. Here, we investigate a simplified model system consisting of two constituents in aqueous solution, namely, of the proteins bovine serum albumin (BSA) and bovine polyclonal gamma-globulin (Ig), systematically depending on the total volume fraction and ratio of these constituents. From high-resolution quasi-elastic neutron spectroscopy, the separate apparent short-time diffusion coefficients for BSA and Ig in the mixture are extracted, which show substantial deviations from the diffusion coefficients measured in monodisperse solutions at the same total volume fraction. These deviations can be modeled quantitatively using results from the short-time rotational and translational diffusion in a two-component hard sphere system with two distinct, effective hydrodynamic radii. Thus, we find that a simple colloid picture well describes short-time diffusion in binary mixtures as a function of the mixing ratio and the total volume fraction. Notably, the self-diffusion of the smaller protein BSA in the mixture is faster than the diffusion in a pure BSA solution, whereas the self-diffusion of Ig in the mixture is slower than in the pure Ig solution.

Molecular Flexibility of Antibodies Preserved Even in the Dense Phase after Macroscopic Phase Separation.

Antibody therapies are typically based on high-concentration formulations that need to be administered subcutaneously. These conditions induce several challenges, inter alia a viscosity suitable for injection, sufficient solution stability, and preservation of molecular function. To obtain systematic insights into the molecular factors, we study the dynamics on the molecular level under strongly varying solution conditions. In particular, we use solutions of antibodies with poly(ethylene glycol), in which simple cooling from room temperature to freezing temperatures induces a transition from a well-dispersed solution into a phase-separated and macroscopically arrested system. Using quasi-elastic neutron scattering during in situ cooling ramps and in prethermalized measurements, we observe a strong decrease in antibody diffusion, while internal flexibility persists to a significant degree, thus ensuring the movement necessary for the preservation of molecular function. These results are relevant for a more dynamic understanding of antibodies in high-concentration formulations, which affects the formation of transient clusters governing the solution viscosity.

Temperature and salt controlled tuning of protein clusters.

The formation of molecular assemblies in protein solutions is of strong interest both from a fundamental viewpoint and for biomedical applications. While ordered and desired protein assemblies are indispensable for some biological functions, undesired protein condensation can induce serious diseases. As a common cofactor, the presence of salt ions is essential for some biological processes involving proteins, and in aqueous suspensions of proteins can also give rise to complex phase diagrams including homogeneous solutions, large aggregates, and dissolution regimes. Here, we systematically study the cluster formation approaching the phase separation in aqueous solutions of the globular protein BSA as a function of temperature (T), the protein concentration (cp) and the concentrations of the trivalent salts YCl3 and LaCl3 (cs). As an important complement to structural, i.e. time-averaged, techniques we employ a dynamical technique that can detect clusters even when they are transient on the order of a few nanoseconds. By employing incoherent neutron spectroscopy, we unambiguously determine the short-time self-diffusion of the protein clusters depending on cp, cs and T. We determine the cluster size in terms of effective hydrodynamic radii as manifested by the cluster center-of-mass diffusion coefficients D. For both salts, we find a simple functional form D(cp, cs, T) in the parameter range explored. The calculated inter-particle attraction strength, determined from the microscopic and short-time diffusive properties of the samples, increases with salt concentration and temperature in the regime investigated and can be linked to the macroscopic behavior of the samples.

ApoE and ApoE Nascent-Like HDL Particles at Model Cellular Membranes: Effect of Protein Isoform and Membrane Composition.

Apolipoprotein E (ApoE), an important mediator of lipid transportation in plasma and the nervous system, plays a large role in diseases such as atherosclerosis and Alzheimer's. The major allele variants ApoE3 and ApoE4 differ only by one amino acid. However, this difference has major consequences for the physiological behaviour of each variant. In this paper, we follow (i) the initial interaction of lipid-free ApoE variants with model membranes as a function of lipid saturation, (ii) the formation of reconstituted High-Density Lipoprotein-like particles (rHDL) and their structural characterisation, and (iii) the rHDL ability to exchange lipids with model membranes made of saturated lipids in the presence and absence of cholesterol [1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) or 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) with and without 20 mol% cholesterol]. Our neutron reflection results demonstrate that the protein variants interact differently with the model membranes, adopting different protein conformations. Moreover, the ApoE3 structure at the model membrane is sensitive to the level of lipid unsaturation. Small-angle neutron scattering shows that the ApoE containing lipid particles form elliptical disc-like structures, similar in shape but larger than nascent or discoidal HDL based on Apolipoprotein A1 (ApoA1). Neutron reflection shows that ApoE-rHDL do not remove cholesterol but rather exchange saturated lipids, as occurs in the brain. In contrast, ApoA1-containing particles remove and exchange lipids to a greater extent as occurs elsewhere in the body.

Self-diffusion of nonspherical particles fundamentally conflicts with effective sphere models.

Modeling diffusion of nonspherical particles presents an unsolved and considerable challenge, despite its importance for the understanding of crowding effects in biology, food technology and formulation science. A common approach in experiment and simulation is to map nonspherical objects on effective spheres to subsequently use the established predictions for spheres to approximate phenomena for nonspherical particles. Using numerical evaluation of the hydrodynamic mobility tensor, we show that this so-called effective sphere model fundamentally fails to represent the self-diffusion in solutions of ellipsoids as well as rod-like assemblies of spherical beads. The effective sphere model drastically overestimates the slowing down of self-diffusion down to volume fractions below 0.01. Furthermore, even the linear term relevant at lower volume fraction is inaccurate, linked to a fundamental misconception of effective sphere models. To overcome the severe problems related with the use of effective sphere models, we suggest a protocol to predict the short-time self-diffusion of rod-like systems, based on simulations with hydrodynamic interactions that become feasible even for more complex molecules as the essential observable shows a negligible system-size effect.

Crowding in the Eye Lens: Modeling the Multisubunit Protein ß-Crystallin with a Colloidal Approach.

We present a multiscale characterization of aqueous solutions of the bovine eye lens protein ßH crystallin from dilute conditions up to dynamical arrest, combining dynamic light scattering, small-angle x-ray scattering, tracer-based microrheology, and neutron spin echo spectroscopy. We obtain a comprehensive explanation of the observed experimental signatures from a model of polydisperse hard spheres with additional weak attraction. In particular, the model predictions quantitatively describe the multiscale dynamical results from microscopic nanometer cage diffusion over mesoscopic micrometer gradient diffusion up to macroscopic viscosity. Based on a comparative discussion with results from other crystallin proteins, we suggest an interesting common pathway for dynamical arrest in all crystallin proteins, with potential implications for the understanding of crowding effects in the eye lens.

Evolution of the structure and dynamics of bovine serum albumin induced by thermal denaturation.

Protein denaturation in concentrated solutions consists of the unfolding of the native protein structure, and subsequent cross-linking into clusters or gel networks. While the kinetic evolution of structure has been studied for some cases, the underlying microscopic dynamics of proteins has so far been neglected. However, protein dynamics is essential to understand the specific nature of assembly processes, such as diffusion-limited growth, or vitrification of dense liquids. Here, we present a study on thermal denaturation of concentrated solutions of bovine serum albumin (BSA) in D2O with and without NaCl. Using small-angle scattering, we provide information on structure before, during and after denaturation. Using quasi-elastic neutron scattering, we monitor in real-time the microscopic dynamics and dynamical confinement throughout the entire denaturation process covering protein unfolding and cross-linking. After denaturation, the protein dynamics is slowed down in salty solutions compared to those in pure water, while the stability and dynamics of the native solution appears unaffected by salt. The approach presented here opens opportunities to link microscopic dynamics to emerging structural properties, with implications for assembly processes in soft and biological matter.

Multivalent ions and biomolecules: Attempting a comprehensive perspective.

Ions are ubiquitous in nature. They play a key role for many biological processes on the molecular scale, from molecular interactions, to mechanical properties, to folding, to self-organisation and assembly, to reaction equilibria, to signalling, to energy and material transport, to recognition etc. Going beyond monovalent ions to multivalent ions, the effects of the ions are frequently not only stronger (due to the obviously higher charge), but qualitatively different. A typical example is the process of binding of multivalent ions, such as Ca2+ , to a macromolecule and the consequences of this ion binding such as compaction, collapse, potential charge inversion and precipitation of the macromolecule. Here we review these effects and phenomena induced by multivalent ions for biological (macro)molecules, from the "atomistic/molecular" local picture of (potentially specific) interactions to the more global picture of phase behaviour including, e. g., crystallisation, phase separation, oligomerisation etc. Rather than attempting an encyclopedic list of systems, we rather aim for an embracing discussion using typical case studies. We try to cover predominantly three main classes: proteins, nucleic acids, and amphiphilic molecules including interface effects. We do not cover in detail, but make some comparisons to, ion channels, colloidal systems, and synthetic polymers. While there are obvious differences in the behaviour of, and the relevance of multivalent ions for, the three main classes of systems, we also point out analogies. Our attempt of a comprehensive discussion is guided by the idea that there are not only important differences and specific phenomena with regard to the effects of multivalent ions on the main systems, but also important similarities. We hope to bridge physico-chemical mechanisms, concepts of soft matter, and biological observations and connect the different communities further.

Protein Short-Time Diffusion in a Naturally Crowded Environment.

The interior of living cells is a dense and polydisperse suspension of macromolecules. Such a complex system challenges an understanding in terms of colloidal suspensions. As a fundamental test we employ neutron spectroscopy to measure the diffusion of tracer proteins (immunoglobulins) in a cell-like environment (cell lysate) with explicit control over crowding conditions. In combination with Stokesian dynamics simulation, we address protein diffusion on nanosecond time scales where hydrodynamic interactions dominate over negligible protein collisions. We successfully link the experimental results on these complex, flexible molecules with coarse-grained simulations providing a consistent understanding by colloid theories. Both experiments and simulations show that tracers in polydisperse solutions close to the effective particle radius Reff = ⟨ Ri3⟩1/3 diffuse approximately as if the suspension was monodisperse. The simulations further show that macromolecules of sizes R > Reff ( R < Reff) are slowed more (less) effectively even at nanosecond time scales, which is highly relevant for a quantitative understanding of cellular processes.

Weak Shape Anisotropy Leads to a Nonmonotonic Contribution to Crowding, Impacting Protein Dynamics under Physiologically Relevant Conditions.

The effect of a nonspherical particle shape on the dynamics in crowded solutions presents a significant challenge for a comprehensive understanding of interaction and structural relaxation in biological and soft matter. We report that small deviations from a spherical shape induce a nonmonotonic contribution to the crowding effect on the short-time cage diffusion compared with spherical systems, using molecular dynamics simulations with mesoscale hydrodynamics of a multiparticle collision dynamics fluid in semidilute systems with volume fractions smaller than 0.35. We show that the nonmonotonic effect due to anisotropy is caused by the combination of a reduced relative mobility over the entire concentration range and a looser and less homogeneous cage packing of nonspherical particles. Our finding stresses that nonsphericity induces new complexity, which cannot be accounted for in effective sphere models, and is of great interest in applications such as formulations as well as for the fundamental understanding of soft matter in general and crowding effects in living cells in particular.