Protein self-diffusion in crowded solutions

Understanding dynamics in coarse-grained models. IV. Connection of fine-grained and coarse-grained dynamics with the Stokes-Einstein and Stokes-Einstein-Debye relations

Applying an excess entropy scaling formalism to the coarse-grained (CG) dynamics of liquids, we discovered that missing rotational motions during the CG process are responsible for artificially accelerated CG dynamics. In the context of the dynamic representability between the fine-grained (FG) and CG dynamics, this work introduces the well-known Stokes-Einstein and Stokes-Einstein-Debye relations to unravel the rotational dynamics underlying FG trajectories, thereby allowing for an indirect evaluation of the effective rotations based only on the translational information at the reduced CG resolution. Since the representability issue in CG modeling limits a direct evaluation of the shear stress appearing in the Stokes-Einstein and Stokes-Einstein-Debye relations, we introduce a translational relaxation time as a proxy to employ these relations, and we demonstrate that these relations hold for the ambient conditions studied in our series of work. Additional theoretical links to our previous work are also established. First, we demonstrate that the effective hard sphere radius determined by the classical perturbation theory can approximate the complex hydrodynamic radius value reasonably well. Furthermore, we present a simple derivation of an excess entropy scaling relationship for viscosity by estimating the elliptical integral of molecules. In turn, since the translational and rotational motions at the FG level are correlated to each other, we conclude that the "entropy-free" CG diffusion only depends on the shape of the reference molecule. Our results and analyses impart an alternative way of recovering the FG diffusion from the CG description by coupling the translational and rotational motions at the hydrodynamic level.

In Silico View of Crowding: Biomolecular Processes to Nanomaterial Design

It is widely accepted that deciphering biomolecular structure and function requires going beyond the single-molecule or single-complex paradigm. The densely packed macromolecules, cosolutes, and metabolites in the living cell impose crowding effects on the biomolecular structure and dynamics that need to be accounted for. Molecular simulations have proven to be a powerful tool to advance the current molecular-level understanding of such a highly concentrated, complex milieu. This Mini-Review focuses on summarizing the understanding achieved so far for the effects of crowding on biomolecular processes using computational methods, along with highlighting a new direction in employing crowding as a tool for tunable nanomaterial design. The two schools of thought that form the pillars of the current understanding of crowding effects are discussed. The investigation of crowded solutions using physics-based models that encompass different time and length scales to mimic the intracellular environment are described. The limitations and challenges faced by the current models and simulation methods are addressed, highlighting the gaps to be filled for better agreement with experiments. Crowding can also act as an effective tool to modulate the structure-property-function relationships of nanomaterials, leading to the development of novel functional materials. A few recent studies, mostly experimental, have been summarized in this direction. The Mini-Review concludes with an outlook for future developments in this field in order to enable accurate mimicking of the intracellular environment using simulations and to bridge the gap between biological processes and nanomaterial design.

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.

High-pressure X-ray photon correlation spectroscopy at fourth-generation synchrotron sources

A new experimental setup combining X-ray photon correlation spectroscopy (XPCS) in the hard X-ray regime and a high-pressure sample environment has been developed to monitor the pressure dependence of the internal motion of complex systems down to the atomic scale in the multi-gigapascal range, from room temperature to 600 K. The high flux of coherent high-energy X-rays at fourth-generation synchrotron sources solves the problems caused by the absorption of diamond anvil cells used to generate high pressure, enabling the measurement of the intermediate scattering function over six orders of magnitude in time, from 10-3 s to 103 s. The constraints posed by the high-pressure generation such as the preservation of X-ray coherence, as well as the sample, pressure and temperature stability, are discussed, and the feasibility of high-pressure XPCS is demonstrated through results obtained on metallic glasses.

Computational Approaches to Predict Protein-Protein Interactions in Crowded Cellular Environments

Investigating protein-protein interactions is crucial for understanding cellular biological processes because proteins often function within molecular complexes rather than in isolation. While experimental and computational methods have provided valuable insights into these interactions, they often overlook a critical factor: the crowded cellular environment. This environment significantly impacts protein behavior, including structural stability, diffusion, and ultimately the nature of binding. In this review, we discuss theoretical and computational approaches that allow the modeling of biological systems to guide and complement experiments and can thus significantly advance the investigation, and possibly the predictions, of protein-protein interactions in the crowded environment of cell cytoplasm. We explore topics such as statistical mechanics for lattice simulations, hydrodynamic interactions, diffusion processes in high-viscosity environments, and several methods based on molecular dynamics simulations. By synergistically leveraging methods from biophysics and computational biology, we review the state of the art of computational methods to study the impact of molecular crowding on protein-protein interactions and discuss its potential revolutionizing effects on the characterization of the human interactome.

Inelastic neutron scattering and spectroscopy methods to characterize dynamics in colloidal and soft matter systems

Colloidal systems and soft materials are well suited to neutron scattering, and the community has readily adopted elastic scattering techniques to investigate their structure. Due to their unique properties, neutrons may also be used to characterize the dynamics of soft materials over a wide range of length and time scales in situ. Both static structures and an understanding of how molecules move about their equilibrium positions is essential if we are to deliver on the promise of rationally designing soft materials. In this review we introduce the basics of neutron spectroscopy and explore the ways in which inelastic neutron scattering can be used to study colloidal and soft materials. Illustrative examples are chosen that highlight the phenomena suitable for investigation using this suite of techniques.

Decoding the Role of the Global Proteome Dynamics for Cellular Thermal Stability

Molecular mechanisms underlying the thermal response of cells remain elusive. On the basis of the recent result that the short-time diffusive dynamics of the Escherichia coli proteome is an excellent indicator of temperature-dependent bacterial metabolism and death, we used neutron scattering (NS) spectroscopy and molecular dynamics (MD) simulations to investigate the sub-nanosecond proteome mobility in psychro-, meso-, and hyperthermophilic bacteria over a wide temperature range. The magnitude of thermal fluctuations, measured by atomic mean square displacements, is similar among all studied bacteria at their respective thermal cell death. Global roto-translational motions turn out to be the main factor distinguishing the bacterial dynamical properties. We ascribe this behavior to the difference in the average proteome net charge, which becomes less negative for increasing bacterial thermal stability. We propose that the chemical-physical properties of the cytoplasm and the global dynamics of the resulting proteome are fine-tuned by evolution to uphold optimal thermal stability conditions.

Effects of Crowding and Cosolutes on Biomolecular Function at Extreme Environmental Conditions

The extent of the effect of cellular crowding and cosolutes on the functioning of proteins and cells is manifold and includes the stabilization of the biomolecular systems, the excluded volume effect, and the modulation of molecular dynamics. Simultaneously, it is becoming increasingly clear how important it is to take the environment into account if we are to shed light on biological function under various external conditions. Many biosystems thrive under extreme conditions, including the deep sea and subseafloor crust, and can take advantage of some of the effects of crowding. These relationships have been studied in recent years using various biophysical techniques, including neutron and X-ray scattering, calorimetry, FTIR, UV-vis and fluorescence spectroscopies. Combining knowledge of the structure and conformational dynamics of biomolecules under extreme conditions, such as temperature, high hydrostatic pressure, and high salinity, we highlight the importance of considering all results in the context of the environment. Here we discuss crowding and cosolute effects on proteins, nucleic acids, membranes, and live cells and explain how it is possible to experimentally separate crowding-induced effects from other influences. Such findings will contribute to a better understanding of the homeoviscous adaptation of organisms and the limits of life in general.

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.

Bridging attraction of condensed bovine serum albumin solution in the presence of trivalent ions: A SANS study

The bridging attraction of condensed bovine serum albumin (BSA) solution (D2O) in the presence of yttrium chloride (YCl3) was studied by small angle neutron scattering (SANS). With increasing the concentration of YCl3 (cY) from 3 to 15 mM and from 15 to 100 mM, the intensity in low-q region increases and then decreases. Combining the tri-axial ellipsoid (TaE) geometry and the multi-component sticky hard sphere (SHS) potential, a SHS-TaE model was established to quantitatively determine the size and distribution of particles. In this way, the structural mechanism of the aggregation-redissolution process in protein solution was demonstrated and discussed. As cY increases from 3 to 100 mM, the SHS radius rL decreases from ca. 2.97 to 2.50 nm, suggesting that the relatively well dispersed BSAs may form aggregates with various polydispersities. The axis a increases from 1.88 to 2.30 nm, while b and c decrease from 3.53 to 3.23 nm and from 4.12 to 3.55 nm, respectively. (RgTaE decreases from ca. 2.57 to 2.38 nm). Moreover, the scattering length density (SLD) of BSA decreases from 3.67 to 1.56 × 10-6 Å-2. All these results consistently indicate a strengthened attraction and the BSA molecules might shrink and tune out to be more like of oblate ellipsoid with increasing the amount of YCl3.

Heterogeneous anomalous transport in cellular and molecular biology

It is well established that a wide variety of phenomena in cellular and molecular biology involve anomalous transport e.g. the statistics for the motility of cells and molecules are fractional and do not conform to the archetypes of simple diffusion or ballistic transport. Recent research demonstrates that anomalous transport is in many cases heterogeneous in both time and space. Thus single anomalous exponents and single generalised diffusion coefficients are unable to satisfactorily describe many crucial phenomena in cellular and molecular biology. We consider advances in the field ofheterogeneous anomalous transport(HAT) highlighting: experimental techniques (single molecule methods, microscopy, image analysis, fluorescence correlation spectroscopy, inelastic neutron scattering, and nuclear magnetic resonance), theoretical tools for data analysis (robust statistical methods such as first passage probabilities, survival analysis, different varieties of mean square displacements, etc), analytic theory and generative theoretical models based on simulations. Special emphasis is made on high throughput analysis techniques based on machine learning and neural networks. Furthermore, we consider anomalous transport in the context of microrheology and the heterogeneous viscoelasticity of complex fluids. HAT in the wavefronts of reaction-diffusion systems is also considered since it plays an important role in morphogenesis and signalling. In addition, we present specific examples from cellular biology including embryonic cells, leucocytes, cancer cells, bacterial cells, bacterial biofilms, and eukaryotic microorganisms. Case studies from molecular biology include DNA, membranes, endosomal transport, endoplasmic reticula, mucins, globular proteins, and amyloids.

Protein translational diffusion as a way to detect intermolecular interactions

In this work, we analyze the information on the protein intermolecular interactions obtained from macromolecular diffusion. We have shown that the most hopeful results are given by our approach based on analysis of protein translational self-diffusion and collective diffusion obtained by dynamic light scattering and pulsed-field gradient NMR (PFG NMR) spectroscopy with the help of Vink's approach to analyze diffusion motion of particles by frictional formalism of non-equilibrium thermodynamics and the usage of the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory of colloid particles interactions in electrolyte solutions. Early we have shown that integration of Vink's theory with DLVO provides a reliable basis for uniform interpreting of PFG NMR and DLS experiments on concentration dependence of diffusion coefficients. Basic details of theoretical and mathematical procedures and a broad analysis of experimental attestation of proposed conception on proteins of various structural form, size, and shape are presented. In the present review, the main capabilities of our approach obtain the details of intermolecular interactions of proteins with different shapes, internal structures, and mass. The universality of Vink's approach is experimentally shown, which gives the appropriate description of experimental results for proteins of complicated structure and shape.

Biophysical Determinants for the Viscosity of Concentrated Monoclonal Antibody Solutions

Monoclonal antibodies (mAbs) are particularly relevant for therapeutics due to their high specificity and versatility, and mAb-based drugs are hence used to treat numerous diseases. The increased patient compliance of self-administration motivates the formulation of products for subcutaneous (SC) administration. The associated challenge is to formulate highly concentrated antibody solutions to achieve a significant therapeutic effect, while limiting their viscosity and preserving their physicochemical stability. Protein-protein interactions (PPIs) are in fact the root cause of several potential problems concerning the stability, manufacturability, and delivery of a drug product. The understanding of macroscopic viscosity requires an in-depth knowledge on protein diffusion, PPIs, and self-association/aggregation. Here, we study the self-diffusion of different mAbs of the IgG1 subtype in aqueous solution as a function of the concentration and temperature by quasi-elastic neutron scattering (QENS). QENS allows us to probe the short-time self-diffusion of the molecules and therefore to determine the hydrodynamic mAb cluster size and to gain information on the internal mAb dynamics. Small-angle neutron scattering (SANS) is jointly employed to probe structural details and to understand the nature and intensity of PPIs. Complementary information is provided by molecular dynamics (MD) simulations and viscometry, thus obtaining a comprehensive picture of mAb diffusion.

Anomalous Colloidal Motion under Strong Confinement

Diffusion of biological macromolecules in the cytoplasm is a paradigm of colloidal diffusion in an environment characterized by a strong restriction of the accessible volume. This makes of the understanding of the physical rules governing colloidal diffusion under conditions mimicking the reduction in accessible volume occurring in the cell cytoplasm, a problem of a paramount importance. This work aims to study how the thermal motion of spherical colloidal beads in the inner cavity of giant unilamellar vesicles (GUVs) is modified by strong confinement conditions, and the viscoelastic character of the medium. Using single particle tracking, it is found that both the confinement and the environmental viscoelasticity lead to the emergence of anomalous motion pathways for colloidal microbeads encapsulated in the aqueous inner cavity of GUVs. This anomalous diffusion is strongly dependent on the ratio between the volume of the colloidal particle and that of the GUV under consideration as well as on the viscosity of the particle's liquid environment. Therefore, the results evidence that the reduction of the free volume accessible to colloidal motion pushes the diffusion far from a standard Brownian pathway as a result of the change in the hydrodynamic boundary conditions driving the particle motion.

Optimized OPEP Force Field for Simulation of Crowded Protein Solutions

Macromolecular crowding has profound effects on the mobility of proteins, with strong implications on the rates of intracellular processes. To describe the dynamics of crowded environments, detailed molecular models are needed, capturing the structures and interactions arising in the crowded system. In this work, we present OPEPv7, which is a coarse-grained force field at amino-acid resolution, suited for rigid-body simulations of the structure and dynamics of crowded solutions formed by globular proteins. Using the OPEP protein model as a starting point, we have refined the intermolecular interactions to match the experimentally observed dynamical slowdown caused by crowding. The resulting force field successfully reproduces the diffusion slowdown in homogeneous and heterogeneous protein solutions at different crowding conditions. Coupled with the lattice Boltzmann technique, it allows the study of dynamical phenomena in protein assemblies and opens the way for the in silico rheology of protein solutions.

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.

Fundamental advances in hydrogels for the development of the next generation of smart delivery systems as biopharmaceuticals

Recent advances in developing and applying therapeutic peptides for anticancer, antimicrobial and immunomodulatory remedies have opened a new era in therapeutics. This development has resulted in the engineering of new biologics as part of a concerted effort by the pharmaceutical industry. Many alternative routes of administration and delivery vehicles, targeting better patient compliance and optimal therapeutic bioavailability, have emerged. However, the design of drug delivery systems to protect a range of unstable macromolecules, including peptides and proteins, from high temperatures, acidic environments, and enzymatic degradation remains a priority. Herein, we give chronological insights in the development of controlled-release drug delivery systems that occurred in the last 70 years or so. Subsequently, we summarise the key physicochemical characteristics of hydrogels contributing to the development of protective delivery systems concerning drug-targeted delivery in the chronospatial domain for biopharmaceuticals. Furthermore, we shed some light on promising hydrogels that can be utilised for systemic bioactive administration.

Effect of crowding, compartmentalization and nanodomains on protein modification and redox signaling - current state and future challenges

Biological milieus are highly crowded and heterogeneous systems where organization of macromolecules within nanodomains (e.g. membraneless compartments) is vital to the regulation of metabolic processes. There is an increasing interest in understanding the effects that such packed environments have on different biochemical and biological processes. In this context, the redox biochemistry and redox signaling fields are moving towards investigating oxidative processes under conditions that exhibit these key features of biological systems in order to solve existing paradigms including those related to the generation and transmission of specific redox signals within and between cells in both normal physiology and under conditions of oxidative stress. This review outlines the effects that crowding, nanodomain formation and altered local viscosities can have on biochemical processes involving proteins, and then discusses some of the reactions and pathways involving proteins and oxidants that may, or are known to, be modulated by these factors. We postulate that knowledge of protein modification processes (e.g. kinetics, pathways and product formation) under conditions that mimic biological milieus, will provide a better understanding of the response of cells to endogenous and exogenous stressors, and their role in ageing, signaling, health and disease.

Incoherent Neutron Scattering and Terahertz Time-Domain Spectroscopy on Protein and Hydration Water

Incoherent inelastic and quasi-elastic neutron scattering (INS) and terahertz time-domain spectroscopy (THz-TDS) are spectroscopy methods that directly detect molecular dynamics, with an overlap in the measured energy regions of each method. Due to the different characteristics of their probes (i.e., neutron and light), the information obtained and the sample conditions suitable for each method differ. In this review, we introduce the differences in the quantum beam properties of the two methods and their associated advantages and disadvantages in molecular spectroscopy. Neutrons are scattered via interaction with nuclei; one characteristic of neutron scattering is a large incoherent scattering cross-section of a hydrogen atom. INS records the auto-correlation functions of atomic positions. By using the difference in neutron scattering cross-sections of isotopes in multi-component systems, some molecules can be selectively observed. In contrast, THz-TDS observes the cross-correlation function of dipole moments. In water-containing biomolecular samples, the absorption of water molecules is particularly large. While INS requires large-scale experimental facilities, such as accelerators and nuclear reactors, THz-TDS can be performed at the laboratory level. In the analysis of water molecule dynamics, INS is primarily sensitive to translational diffusion motion, while THz-TDS observes rotational motion in the spectrum. The two techniques are complementary in many respects, and a combination of the two is very useful in analyzing the dynamics of biomolecules and hydration water.

Light-Induced Conformational Flexibility of the Orange Carotenoid Protein Studied by Quasielastic Neutron Scattering with In Situ Illumination

The orange carotenoid protein plays a vital role in the photoprotection of cyanobacteria and exhibits a significant structural change upon photoactivation. A rarely considered aspect is the importance of internal protein dynamics in facilitating the structural transition to the active state. In this study, we use quasielastic neutron scattering under (in situ) blue light illumination for the first time to directly probe the protein dynamics of the orange carotenoid protein in the dark-adapted and active states. This shows that the localized internal dynamics of amino acid residues is significantly enhanced upon photoactivation. This is attributed to the photoinduced structural changes exposing larger areas of the protein surface to the solvent, thus resulting in a higher degree of motional freedom. However, the flexibility of the W288A mutant assumed to mimic the active state structure is found to be different, thus highlighting the importance of in situ experiments.

Effects of External Perturbations on Protein Systems: A Microscopic View

Protein folding can be viewed as the origami engineering of biology resulting from the long process of evolution. Even decades after its recognition, research efforts worldwide focus on demystifying molecular factors that underlie protein structure-function relationships; this is particularly relevant in the era of proteopathic disease. A complex co-occurrence of different physicochemical factors such as temperature, pressure, solvent, cosolvent, macromolecular crowding, confinement, and mutations that represent realistic biological environments are known to modulate the folding process and protein stability in unique ways. In the current review, we have contextually summarized the substantial efforts in unveiling individual effects of these perturbative factors, with major attention toward bottom-up approaches. Moreover, we briefly present some of the biotechnological applications of the insights derived from these studies over various applications including pharmaceuticals, biofuels, cryopreservation, and novel materials. Finally, we conclude by summarizing the challenges in studying the combined effects of multifactorial perturbations in protein folding and refer to complementary advances in experiment and computational techniques that lend insights to the emergent challenges.

Crowding and confinement act in concert to slow DNA diffusion within cell-sized droplets

Dynamics of biological macromolecules, such as DNA, in crowded and confined environments are critical to understanding cellular processes such as transcription, infection, and replication. However, the combined effects of cellular confinement and crowding on macromolecular dynamics remain poorly understood. Here, we use differential dynamic microscopy to investigate the diffusion of large DNA molecules confined in cell-sized droplets and crowded by dextran polymers. We show that confined and crowded DNA molecules exhibit universal anomalous subdiffusion with scaling that is insensitive to the degree of confinement and crowding. However, effective DNA diffusion coefficients Deff decrease up to 2 orders of magnitude as droplet size decreases—an effect that is enhanced by increased crowding. We mathematically model the coupling of crowding and confinement by combining polymer scaling theories with confinement-induced depletion effects. The generality and tunability of our system and models render them applicable to elucidating wide-ranging crowded and confined systems.

•

DNA diffusion measured in cell-sized droplets with differential dynamic microscopy

•

Combination of crowding and confinement leads to subdiffusion and slowing

•

Diffusion coefficients of DNA decrease strongly with decreasing droplet size

•

Polymer scaling theories and depletion effects predict observed dynamics

Resolving molecular diffusion and aggregation of antibody proteins with megahertz X-ray free-electron laser pulses

X-ray free-electron lasers (XFELs) with megahertz repetition rate can provide novel insights into structural dynamics of biological macromolecule solutions. However, very high dose rates can lead to beam-induced dynamics and structural changes due to radiation damage. Here, we probe the dynamics of dense antibody protein (Ig-PEG) solutions using megahertz X-ray photon correlation spectroscopy (MHz-XPCS) at the European XFEL. By varying the total dose and dose rate, we identify a regime for measuring the motion of proteins in their first coordination shell, quantify XFEL-induced effects such as driven motion, and map out the extent of agglomeration dynamics. The results indicate that for average dose rates below 1.06 kGy μs^-1 in a time window up to 10 μs, it is possible to capture the protein dynamics before the onset of beam induced aggregation. We refer to this approach as correlation before aggregation and demonstrate that MHz-XPCS bridges an important spatio-temporal gap in measurement techniques for biological samples.

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.

Structure and diffusive dynamics of aspartate α-decarboxylase (ADC) liganded with D-serine in aqueous solution

Incoherent neutron spectroscopy, in combination with dynamic light scattering, was used to investigate the effect of ligand binding on the center-of-mass self-diffusion and internal diffusive dynamics of Escherichia coli aspartate α-decarboxylase (ADC). The X-ray crystal structure of ADC in complex with the d-serine inhibitor was also determined, and molecular dynamics simulations were used to further probe the structural rearrangements that occur as a result of ligand binding. These experiments reveal that d-serine forms hydrogen bonds with some of the active site residues, that higher order oligomers of the ADC tetramer exist on ns–ms time-scales, and also show that ligand binding both affects the ADC internal diffusive dynamics and appears to further increase the size of the higher order oligomers.

Neutron spectroscopy, dynamic light scattering, X-ray diffraction, and MD-simulations were used to investigate the effect of ligand binding on the structure and diffusive dynamics of Escherichia coli aspartate alpha-decarboxylase.

Effect of Protein–Protein Interactions on Translational Diffusion of Spheroidal Proteins

One of the commonly accepted approaches to estimate protein–protein interactions (PPI) in aqueous solutions is the analysis of their translational diffusion. The present review article observes a phenomenological approach to analyze PPI effects via concentration dependencies of self- and collective translational diffusion coefficient for several spheroidal proteins derived from the pulsed field gradient NMR (PFG NMR) and dynamic light scattering (DLS), respectively. These proteins are rigid globular α-chymotrypsin (ChTr) and human serum albumin (HSA), and partly disordered α-casein (α-CN) and β-lactoglobulin (β-Lg). The PPI analysis enabled us to reveal the dominance of intermolecular repulsion at low ionic strength of solution (0.003–0.01 M) for all studied proteins. The increase in the ionic strength to 0.1–1.0 M leads to the screening of protein charges, resulting in the decrease of the protein electrostatic potential. The increase of the van der Waals potential for ChTr and α-CN characterizes their propensity towards unstable weak attractive interactions. The decrease of van der Waals interactions for β-Lg is probably associated with the formation of stable oligomers by this protein. The PPI, estimated with the help of interaction potential and idealized spherical molecular geometry, are in good agreement with experimental data.

Multi-length scale structural investigation of lysozyme self-assembly

Reactive amyloid oligomers are responsible for cytotoxicity in amyloid pathologies and because of their unstable nature characterizing their behavior is a challenge. The physics governing the self-assembly of proteins in crowded conditions is extremely complex and its comprehension, despite its paramount relevance to understanding molecular mechanisms inside cells and optimizing pharmaceutical processes, remains inconclusive. Here, we focus on the amyloid oligomerization process in self-crowded lysozyme aqueous solutions in acidic conditions. We reveal that the amyloid oligomers form at high protein concentration and low pH. Through multi-length scale spectroscopic investigations, we find that amyloid oligomers can further interconnect with each other by weak and non-specific interactions forming an extended network that leads to the percolation of the whole system. Our multi-length scale structural analysis follows the thermal history of amyloid oligomers from different perspectives and highlights the impact of hierarchical self-assembly of biological macromolecules on functional properties.

•

Use of multi-length scale spectroscopies to characterize unstable amyloid oligomers

•

Lysozyme form thermo-labile amyloid oligomers in self-crowded conditions

•

Amyloid oligomers interact and form an extended hydrogel network

•

Amyloid oligomers are responsible for the existence of the hydrogel matrix

Dynamical Behavior of Disordered Regions in Disease-Related Proteins Revealed by Quasielastic Neutron Scattering

Background and Objectives: Intrinsically disordered proteins (IDPs) and proteins containing intrinsically disordered regions (IDRs) are known to be involved in various human diseases. Since the IDPs/IDRs are fluctuating between many structural substrates, the dynamical behavior of the disease-related IDPs/IDRs needs to be characterized to elucidate the mechanisms of the pathogenesis of the diseases. As protein motions have a hierarchy ranging from local side-chain motions, through segmental motions of loops or disordered regions, to diffusive motions of entire molecules, segmental motions, as well as local motions, need to be characterized. Materials and Methods: Combined analysis of quasielastic neutron scattering (QENS) spectra with the structural data provides information on both the segmental motions and the local motions of the IDPs/IDRs. Here, this method is applied to re-analyze the QENS spectra of the troponin core domain (Tn-CD), various mutants of which cause the pathogenesis of familial cardiomyopathy (FCM), and α-synuclein (αSyn), amyloid fibril formation of which is closely related to the pathogenesis of Parkinson's disease, collected in the previous studies. The dynamical behavior of wild-type Tn-CD, FCM-related mutant Tn-CD, and αSyn in the different propensity states for fibril formation is characterized. Results: In the Tn-CD, the behavior of the segmental motions is shown to be different between the wild type and the mutant. This difference is likely to arise from changes in the intramolecular interactions, which are suggested to be related to the functional aberration of the mutant Tn-CD. In αSyn, concerted enhancement of the segmental motions and the local motions is observed with an increased propensity for fibril formation, suggesting the importance of these motions in fibril formation. Conclusions: Characterization of the segmental motions as well as the local motions is thus useful for discussing how the changes in dynamical behavior caused by the disease-related mutations and/or environmental changes could be related to the functional and/or behavioral aberrations of these proteins.

Bioinspired Ice-Binding Materials for Tissue and Organ Cryopreservation

Cryopreservation of tissues and organs can bring transformative changes to medicine and medical science. In the past decades, limited progress has been achieved, although cryopreservation of tissues and organs has long been intensively pursued. One key reason is that the cryoprotective agents (CPAs) currently used for cell cryopreservation cannot effectively preserve tissues and organs because of their cytotoxicity and tissue destructive effect as well as the low efficiency in controlling ice formation. In stark contrast, nature has its unique ways of controlling ice formation, and many living organisms can effectively prevent freezing damage. Ice-binding proteins (IBPs) are regarded as the essential materials identified in these living organisms for regulating ice nucleation and growth. Note that controversial results have been reported on the utilization of IBPs and their mimics for the cryopreservation of tissues and organs, that is, some groups revealed that IBPs and mimics exhibited unique superiorities in tissues cryopreservation, while other groups showed detrimental effects. In this perspective, we analyze possible reasons for the controversy and predict future research directions in the design and construction of IBP inspired ice-binding materials to be used as new CPAs for tissue cryopreservation after briefly introducing the cryo-injuries and the challenges of conventional CPAs in the cryopreservation of tissues and organs.

Self-Diffusive Properties of the Intrinsically Disordered Protein Histatin 5 and the Impact of Crowding Thereon: A Combined Neutron Spectroscopy and Molecular Dynamics Simulation Study

Intrinsically disordered proteins (IDPs) are proteins that, in comparison with globular/structured proteins, lack a distinct tertiary structure. Here, we use the model IDP, Histatin 5, for studying its dynamical properties under self-crowding conditions with quasi-elastic neutron scattering in combination with full atomistic molecular dynamics (MD) simulations. The aim is to determine the effects of crowding on the center-of-mass diffusion as well as the internal diffusive behavior. The diffusion was found to decrease significantly, which we hypothesize can be attributed to some degree of aggregation at higher protein concentrations, (≥100 mg/mL), as indicated by recent small-angle X-ray scattering studies. Temperature effects are also considered and found to, largely, follow Stokes–Einstein behavior. Simple geometric considerations fail to accurately predict the rates of diffusion, while simulations show semiquantitative agreement with experiments, dependent on assumptions of the ratio between translational and rotational diffusion. A scaling law that previously was found to successfully describe the behavior of globular proteins was found to be inadequate for the IDP, Histatin 5. Analysis of the MD simulations show that the width of the distribution with respect to diffusion is not a simplistic mirroring of the distribution of radius of gyration, hence, displaying the particular features of IDPs that need to be accounted for.

Non-Gaussian, transiently anomalous, and ergodic self-diffusion of flexible dumbbells in crowded two-dimensional environments: Coupled translational and rotational motions

We employ Langevin-dynamics simulations to unveil non-Brownian and non-Gaussian center-of-mass self-diffusion of massive flexible dumbbell-shaped particles in crowded two-dimensional solutions. We study the intradumbbell dynamics of the relative motion of the two constituent elastically coupled disks. Our main focus is on effects of the crowding fraction ϕ and of the particle structure on the diffusion characteristics. We evaluate the time-averaged mean-squared displacement (TAMSD), the displacement probability-density function (PDF), and the displacement autocorrelation function (ACF) of the dimers. For the TAMSD at highly crowded conditions of dumbbells, e.g., we observe a transition from the short-time ballistic behavior, via an intermediate subdiffusive regime, to long-time Brownian-like spreading dynamics. The crowded system of dimers exhibits two distinct diffusion regimes distinguished by the scaling exponent of the TAMSD, the dependence of the diffusivity on ϕ, and the features of the displacement-ACF. We attribute these regimes to a crowding-induced transition from viscous to viscoelastic diffusion upon growing ϕ. We also analyze the relative motion in the dimers, finding that larger ϕ suppress their vibrations and yield strongly non-Gaussian PDFs of rotational displacements. For the diffusion coefficients D(ϕ) of translational and rotational motion of the dumbbells an exponential decay with ϕ for weak and a power-law variation D(ϕ)∝(ϕ-ϕ^{★})^{2.4} for strong crowding is found. A comparison of simulation results with theoretical predictions for D(ϕ) is discussed and some relevant experimental systems are overviewed.

Notes on Fitting and Analysis Frameworks for QENS Spectra of (Soft) Colloid Suspensions

With continuously improving signal-to-noise ratios, a statistically sound analysis of quasi-elastic neutron scattering (QENS) spectra requires to fit increasingly complex models which poses several challenges. Simultaneous fits of the spectra for all recorded values of the momentum transfer become a standard approach. Spectrometers at spallation sources can have a complicated non-Gaussian resolution function which has to be described most accurately. At the same time, to speed up the fitting, an analytical convolution with this resolution function is of interest. Here, we discuss basic concepts to efficient approaches for fits of QENS spectra based on standard MATLAB and Python fit algorithms. We illustrate the fits with example data from IN16B, BASIS, and BATS.

Ligand Dynamics in Nanocrystal Solids Studied with Quasi-Elastic Neutron Scattering

Nanocrystal surfaces are commonly populated by organic ligands, which play a determining role in the optical, electronic, thermal, and catalytic properties of the individual nanocrystals and their assemblies. Understanding the bonding of ligands to nanocrystal surfaces and their dynamics is therefore important for the optimization of nanocrystals for different applications. In this study, we use temperature-dependent, quasi-elastic neutron scattering (QENS) to investigate the dynamics of different surface bound alkanethiols in lead sulfide nanocrystal solids. We select alkanethiols with mono- and dithiol terminations, as well as different backbone types and lengths. QENS spectra are collected both on a time-of-flight spectrometer and on a backscattering spectrometer, allowing us to investigate ligand dynamics in a time range from a few picoseconds to nanoseconds. Through model-based analysis of the QENS data, we find that ligands can either (1) precess around a central axis, while simultaneously rotating around their own molecular axis, or (2) only undergo uniaxial rotation with no precession. We establish the percentage of ligands undergoing each type of motion, the average relaxation times, and activation energies for these motions. We determine, for example, that dithiols which link facets of neighboring nanocrystals only exhibit uniaxial rotation and that longer ligands have higher activation energies and show smaller opening angles of precession due to stronger ligand-ligand interactions. Generally, this work provides insight into the arrangement and dynamics of ligands in nanocrystal solids, which is key to understanding their mechanical and thermal properties, and, more generally, highlights the potential of QENS for studying ligand behavior.

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 (c_p) and the concentrations of the trivalent salts YCl₃ and LaCl₃ (c_s). 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 c_p, c_s 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(c_p, c_s, 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.

Conformational dynamics of a multidomain protein by neutron scattering and computational analysis

The flexible conformations of a multidomain protein are responsible for its biological functions. Although MurD, a 47-kDa protein that consists of three domains, sequentially changes its domain conformation from an open form to a closed form through a semiclosed form in its enzymatic reaction, the domain dynamics in each conformation remains unclear. In this study, we verify the conformational dynamics of MurD in the corresponding three states (apo and ATP- and inhibitor-bound states) with a combination of small-angle x-ray and neutron scattering (SAXS and SANS), dynamic light scattering (DLS), neutron backscattering (NBS), neutron spin echo (NSE) spectroscopy, and molecular dynamics (MD) simulations. Applying principal component analysis of the MD trajectories, twisting and open-closed domain modes are identified as the major collective coordinates. The deviations of the experimental SAXS profiles from the theoretical calculations based on the known crystal structures become smaller in the ATP-bound state than in the apo state, and a further decrease is evident upon inhibitor binding. These results suggest that domain motions of the protein are suppressed step by step of each ligand binding. The DLS and NBS data yield collective and self-translational diffusion constants, respectively, and we used them to extract collective domain motions in nanometer and nanosecond scales from the NSE data. In the apo state, MurD shows both twisting and open-closed domain modes, whereas an ATP binding suppresses twisting domain motions, and a further reduction of open-closed mode is seen in the inhibitor-binding state. These observations are consistent with the structure modifications measured by the small-angle scattering as well as the MD simulations. Such changes in the domain dynamics associated with the sequential enzymatic reactions should be related to the affinity and reaction efficiency with a ligand that binds specifically to each reaction state.

Effects of Macromolecular Crowding on Nanoparticle Diffusion: New Insights from Mössbauer Spectroscopy

In this work, we used Mössbauer spectroscopy as a new approach for experimental quantification of the self-diffusion coefficient (D_Mössbauer) and hydrodynamic (HD) size of iron-containing nanoparticles (NPs) in complex crowded solutions, mimicking cell cytoplasm. As a probe, we used 9 nm cobalt ferrite NPs (CFNs) dispersed in solutions of bovine serum albumin (BSA) with a volume fraction (φ_BSA) of 0-0.2. Our results show that the broadening of Mössbauer spectra is highly sensitive to the diffusion of CFNs, while when φ_BSA = 0.2, the CFN-normalized diffusivity is reduced by 86% compared to that of a protein-free solution. CFN colloids were also studied by dynamic light scattering (DLS). Comparison of the experimental data shows that DLS significantly underestimates the diffusion coefficient of CFNs and, consequently, overestimates the HD size of CFNs at φ_BSA > 0, which cannot be attributed to the formation of the BSA monolayer on the surface of CFNs.

Computational Resources and Strategies to Construct Single-Molecule Models of FISH

Currently, the interactions occurring between oligonucleotides and the cellular envelope of bacteria are not fully resolved at the molecular level. Understanding these interactions is essential to gain insights on how to improve the internalization of the tagged oligonucleotides during fluorescence in situ hybridization (FISH). Agent-based modeling (ABM) is a promising in silico tool to dynamically simulate FISH and bring forward new knowledge on this process. Notably, it is important to simulate the whole bacterial cell, including the different layers of the cell envelope, given that the oligonucleotide must cross the envelope to reach its target in the cytosol. In addition, it is also important to characterize other molecules in the cell to best emulate the cell and represent molecular crowding. Here, we review the main information that should be compiled to construct an ABM on FISH and provide a practical example of an oligonucleotide targeting the 23S rRNA of Escherichia coli .

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.

Towards molecular movies with X-ray photon correlation spectroscopy

In this perspective article we highlight research opportunities and challenges in probing structural dynamics of molecular systems using X-ray Photon Correlation Spectroscopy (XPCS). The development of new X-ray sources, such as 4th generation storage rings and X-ray free-electron lasers (XFELs), provides promising new insights into molecular motion. Employing XPCS at these sources allows to capture a very broad range of timescales and lengthscales, spanning from femtoseconds to minutes and atomic scales to the mesoscale. Here, we discuss the scientific questions that can be addressed with these novel tools for two prominent examples: the dynamics of proteins in biomolecular condensates and the dynamics of supercooled water. Finally, we provide practical tips for designing and estimating feasibility of XPCS experiments as well as on detecting and mitigating radiation damage.

Revealing and Attenuating the Electrostatic Properties of Tubulin and Its Polymers

Tubulin is an electrostatically negative protein that forms cylindrical polymers termed microtubules, which are crucial for a variety of intracellular roles. Exploiting the electrostatic behavior of tubulin and microtubules within functional microfluidic and optoelectronic devices is limited due to the lack of understanding of tubulin behavior as a function of solvent composition. This work displays the tunability of tubulin surface charge using dimethyl sulfoxide (DMSO) for the first time. Increasing the DMSO volume fractions leads to the lowering of tubulin's negative surface charge, eventually causing it to become positive in solutions >80% DMSO. As determined by electrophoretic mobility measurements, this change in surface charge is directionally reversible, i.e., permitting control between -1.5 and + 0.2 cm²  (V s)^-1 . When usually negative microtubules are exposed to these conditions, the positively charged tubulin forms tubulin sheets and aggregates, as revealed by an electrophoretic transport assay. Fluorescence-based experiments also indicate that tubulin sheets and aggregates colocalize with negatively charged g-C₃ N₄ sheets while microtubules do not, further verifying the presence of a positive surface charge. This study illustrates that tubulin and its polymers, in addition to being mechanically robust, are also electrically tunable.

Single-Nanoparticle Orientation Sensing by Deep Learning

This paper describes a computational imaging platform to determine the orientation of anisotropic optical probes under differential interference contrast (DIC) microscopy. We established a deep-learning model based on data sets of DIC images collected from metal nanoparticle optical probes at different orientations. This model predicted the in-plane angle of gold nanorods with an error below 20°, the inherent limit of the DIC method. Using low-symmetry gold nanostars as optical probes, we demonstrated the detection of in-plane particle orientation in the full 0-360° range. We also showed that orientation predictions of the same particle were consistent even with variations in the imaging background. Finally, the deep-learning model was extended to enable simultaneous prediction of in-plane and out-of-plane rotation angles for a multibranched nanostar by concurrent analysis of DIC images measured at multiple wavelengths.

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.

Effect of red blood cell shape changes on haemoglobin interactions and dynamics: a neutron scattering study

By using a combination of experimental neutron scattering techniques, it is possible to obtain a statistical perspective on red blood cell (RBC) shape in suspensions, and the inter-relationship with protein interactions and dynamics inside the confinement of the cell membrane. In this study, we examined the ultrastructure of RBC and protein-protein interactions of haemoglobin (Hb) in them using ultra-small-angle neutron scattering and small-angle neutron scattering (SANS). In addition, we used the neutron backscattering method to access Hb motion on the ns time scale and Å length scale. Quasi-elastic neutron scattering (QENS) experiments were performed to measure diffusive motion of Hb in RBCs and in an RBC lysate. By using QENS, we probed both internal Hb dynamics and global protein diffusion, on the accessible time scale and length scale by QENS. Shape changes of RBCs and variation of intracellular Hb concentration were induced by addition of the Na⁺-selective ionophore monensin and the K⁺-selective one, valinomycin. The experimental SANS and QENS results are discussed within the framework of crowded protein solutions, where free motion of Hb is obstructed by mutual interactions.

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 D₂O 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.

Packing and dynamics of a protein solution approaching the jammed state

The packing of proteins and their collective behavior in crowded media is crucial for the understanding of biological processes. Here we study the structural and dynamical evolution of solutions of the globular protein bovine serum albumin with increasing concentration via drying using small angle X-ray scattering and dynamic light scattering. We probe an evolving correlation peak on the scattering profile, corresponding to the inter-protein distance, ξ, which decreases following a power law of the protein volume fraction, φ. The rate of decrease in ξ becomes faster above a protein concentration of ∼200 mg ml-1 (φ = 0.15). The power law exponent changes from 0.33, which is typical of colloidal or protein solutions, to 0.41. During the entire drying process, we observe the development and the growth of two-step relaxation dynamics with increasing φ as revealed by dynamic light scattering. We find three different regimes of the dependence of ξ as a function of φ. In the dilute regime (φ < 0.22), protein molecules are far apart from each other compared to their size. In this case, the dynamics mainly corresponds to Brownian motion. At an intermediate concentration (0.22 < φ < 0.47), inter-protein distances become comparable to the size of protein molecules, leading to a preferential orientation of the ellipsoidal protein molecules along with a possible deformation. In this regime, the dynamics shows two distinct relaxation times. At a very high concentration (φ > 0.47), the system reaches a jammed state. Subsequently, the secondary relaxation time in this state becomes extremely slow. In this state, the protein molecules have approximately one hydration layer. This study contributes to the understanding of protein molecular packing in crowded environments and the phenomenon of density-driven jamming for soft matter systems.

Tracking Internal and Global Diffusive Dynamics During Protein Aggregation by High-Resolution Neutron Spectroscopy

Proteins can misfold and form either amorphous or organized aggregates with different morphologies and features. Aggregates of amyloid nature are pathological hallmarks in so-called protein conformational diseases, including Alzheimer's and Parkinson's. Evidence prevails that the transient early phases of the reaction determine the aggregate morphology and toxicity. As a consequence, real-time monitoring of protein aggregation is of utmost importance. Here, we employed time-resolved neutron backscattering spectroscopy to follow center-of-mass self-diffusion and nano- to picosecond internal dynamics of lysozyme during aggregation into a specific β-sheet rich superstructure, called particulates, formed at the isoelectric point of the protein. Particulate formation is found to be a one-step process, and protein internal dynamics, to remain unchanged during the entire aggregation process. The time-resolved neutron backscattering spectroscopy approach developed here, in combination with standard kinetics assays, provides a unifying framework in which dynamics and conformational transitions can be related to the different aggregation pathways.

Multiscale Coarse-Grained Approach to Investigate Self-Association of Antibodies

Self-association of therapeutic monoclonal antibodies (mabs) are thought to modulate the undesirably high viscosity observed in their concentrated solutions. Computational prediction of such a self-association behavior is advantageous early during mab drug candidate selection when material availability is limited. Here, we present a coarse-grained (CG) simulation method that enables microsecond molecular dynamics simulations of full-length antibodies at high concentrations. The proposed approach differs from others in two ways: first, charges are assigned to CG beads in an effort to reproduce molecular multipole moments and charge asymmetry of full-length antibodies instead of only localized charges. This leads to great improvements in the agreement between CG and all-atom electrostatic fields. Second, the distinctive hydrophobic character of each antibody is incorporated through empirical adjustments to the short-range van der Waals terms dictated by cosolvent all-atom molecular dynamics simulations of antibody variable regions. CG simulations performed on a set of 15 different mabs reveal that diffusion coefficients in crowded environments are markedly impacted by intermolecular interactions. Diffusion coefficients computed from the simulations are in correlation with experimentally measured observables, including viscosities at a high concentration. Further, we show that the evaluation of electrostatic and hydrophobic characters of the mabs is useful in predicting the nonuniform effect of salt on the viscosity of mab solutions. This CG modeling approach is particularly applicable as a material-free screening tool for selecting antibody candidates with desirable viscosity properties.

Potential and limits of a colloid approach to protein solutions

Looking at globular proteins with the eyes of a colloid scientist has a long tradition, in fact a significant part of the early colloid literature was focused on protein solutions. However, it has also been recognized that proteins are much more complex than the typical hard sphere-like synthetic model colloids. Proteins are not perfect spheres, their interaction potentials are in general not isotropic, and using theories developed for such particles are thus clearly inadequate in many cases. In this perspective article, we now take a closer look at the field. In particular, we reflect on the fact that modern colloid science has been undergoing a tremendous development, where a multitude of novel systems have been developed in the lab and in silico. During the last decade we have seen a rapidly increasing number of reports on the synthesis of anisotropic, patchy and/or responsive synthetic colloids, that start to resemble their complex biological counterparts. This experimental development is also reflected in a corresponding theoretical and simulation effort. The experimental and theoretical toolbox of colloid science has thus rapidly expanded, and there is obviously an enormous potential for an application of these new concepts to protein solutions, which has already been realized and harvested in recent years. In this perspective article we make an attempt to critically discuss the exploitation of colloid science concepts to better understand protein solutions. We not only consider classical applications such as the attempt to understand and predict solution stability and phase behaviour, but also discuss new challenges related to the dynamics, flow behaviour and liquid-solid transitions found in concentrated or crowded protein solutions. It not only aims to provide an overview on the progress in experimental and theoretical (bio)colloid science, but also discusses current shortcomings in our ability to correctly reproduce and predict the structural and dynamic properties of protein solutions based on such a colloid approach.