Universality in the morphology and mechanics of coarsening amyloid fibril networks

Microrheology of semiflexible filament solutions based on relaxation simulations

We present an efficient computational methodology to obtain the viscoelastic response of dilute solutions of semiflexible filaments. By considering an approach based on the fluctuation-dissipation theorem, we were able to evaluate the dynamical properties of probe particles immersed in solutions of semiflexible filaments from relaxation simulations with a relatively low computational cost and higher precision in comparison to those based on stochastic dynamics. We used a microrheological approach to obtain the complex shear modulus and the complex viscosity of the solution through its compliance which was obtained directly from the dynamical properties of a probe particle attached to an effective medium described by a mesoscopic model, i.e., an effective filament model (EFM). The relaxation simulations were applied to assess the effects of the bending energy on the viscoelasticity of the semiflexible filament solutions, and our methodology was validated by comparing the numerical results to the experimental data on DNA and collagen solutions.

On the relationship between the plateau modulus and the threshold frequency in peptide gels

Relations between static and dynamic viscoelastic responses in gels can be very elucidating and may provide useful tools to study the behavior of bio-materials such as protein hydrogels. An important example comes from the viscoelasticity of semisolid gel-like materials, which is characterized by two regimes: a low-frequency regime, where the storage modulus G'(ω) displays a constant value G_eq, and a high-frequency power-law stiffening regime, where G'(ω) ∼ ωⁿ. Recently, by considering Monte Carlo simulations to study the formation of peptides networks, we found an intriguing and somewhat related power-law relationship between the plateau modulus and the threshold frequency, i.e., G_eq∼(ω^*)^Δ with Δ = 2/3. Here we present a simple theoretical approach to describe that relationship and test its validity by using experimental data from a β-lactoglobulin gel. We show that our approach can be used even in the coarsening regime where the fractal model fails. Remarkably, the very same exponent Δ is found to describe the experimental data.

Dilute Semiflexible Polymers with Attraction: Collapse, Folding and Aggregation

We review the current state on the thermodynamic behavior and structural phases of self- and mutually-attractive dilute semiflexible polymers that undergo temperature-driven transitions. In extreme dilution, polymers may be considered isolated, and this single polymer undergoes a collapse or folding transition depending on the internal structure. This may go as far as to stable knot phases. Adding polymers results in aggregation, where structural motifs again depend on the internal structure. We discuss in detail the effect of semiflexibility on the collapse and aggregation transition and provide perspectives for interesting future investigations.

Dynamics of protein aggregation and oligomer formation governed by secondary nucleation

The formation of aggregates in many protein systems can be significantly accelerated by secondary nucleation, a process where existing assemblies catalyse the nucleation of new species. In particular, secondary nucleation has emerged as a central process controlling the proliferation of many filamentous protein structures, including molecular species related to diseases such as sickle cell anemia and a range of neurodegenerative conditions. Increasing evidence suggests that the physical size of protein filaments plays a key role in determining their potential for deleterious interactions with living cells, with smaller aggregates of misfolded proteins, oligomers, being particularly toxic. It is thus crucial to progress towards an understanding of the factors that control the sizes of protein aggregates. However, the influence of secondary nucleation on the time evolution of aggregate size distributions has been challenging to quantify. This difficulty originates in large part from the fact that secondary nucleation couples the dynamics of species distant in size space. Here, we approach this problem by presenting an analytical treatment of the master equation describing the growth kinetics of linear protein structures proliferating through secondary nucleation and provide closed-form expressions for the temporal evolution of the resulting aggregate size distribution. We show how the availability of analytical solutions for the full filament distribution allows us to identify the key physical parameters that control the sizes of growing protein filaments. Furthermore, we use these results to probe the dynamics of the populations of small oligomeric species as they are formed through secondary nucleation and discuss the implications of our work for understanding the factors that promote or curtail the production of these species with a potentially high deleterious biological activity.

Importance of non-affine viscoelastic response in disordered fibre networks

Disordered fibre networks are ubiquitous in nature and have a wide range of industrial applications as novel biomaterials. Predicting their viscoelastic response is straightforward for affine deformations that are uniform over all length scales, but when affinity fails, as has been observed experimentally, modelling becomes challenging. Here we present a numerical methodology, related to an existing framework for amorphous packings, to predict the steady-state viscoelastic spectra and degree of affinity for disordered fibre networks driven at arbitrary frequencies. Applying this method to a peptide gel model reveals a monotonic increase of the shear modulus as the soft, non-affine normal modes are successively suppressed as the driving frequency increases. In addition to being dominated by fibril bending, these low frequency network modes are also shown to be delocalised. The presented methodology provides insights into the importance of non-affinity in the viscoelastic response of peptide gels, and is easily extendible to all types of fibre networks.

Amyloid Fibril Solubility

It is well established that amyloid fibril solubility is protein specific, but how solubility depends on the interactions between the fibril building blocks is not clear. Here we use a simple protein model and perform Monte Carlo simulations to directly measure the solubility of amyloid fibrils as a function of the interaction between the fibril building blocks. Our simulations confirms that the fibril solubility depends on the fibril thickness and that the relationship between the interactions and the solubility can be described by a simple analytical formula. The results presented in this study reveal general rules how side-chain-side-chain interactions, backbone hydrogen bonding, and temperature affect amyloid fibril solubility, which might prove to be a powerful tool to design protein fibrils with desired solubility and aggregation properties in general.