Arrested dynamics in a model peptide hydrogel system

We report here on a peptide hydrogel system, which in contrast to most other such systems, is made up of relatively short ﬁbrillar aggregates, discussing resemblance with colloidal rods. The synthetic model peptides A 8 K and A 10 K, where A denotes alanine and K lysine, self-assemble in aqueous solutions into ribbon-like aggregates having an average length (cid:104) L (cid:105) on the order of 100 nm and with a diameter d ≈ 6 nm. The aggregates can be seen as weakly charged rigid rods and they undergo an isotropic to nematic phase transition at higher concentrations. Translational motion perpendicular to the rod axis gets strongly hindered when the concentration is increased above the overlap concentration. Similarly, the rotational motion is hindered, leading to very long stress relaxation times. The peptide self-assembly is driven by hydrophobic interactions and due to a net peptide charge the system is colloidally stable. However, at the same time short range, presumably hydrophobic, attractive interactions appear to affect the rheology of the system. Upon screening the long range electrostatic repulsion, with the addition of salt, the hydrophobic attraction becomes more dominant and we observe a transition from a repulsive glassy state to an attractive gel-state of the rod-like peptide aggregates.


Introduction
Fiber dispersions are interesting due to the strong rheological effects that can be observed already at very low concentrations 1 . In fact, gelation or glass formation can be obtained at a solid content of only a fraction of percent [2][3][4][5][6] . Materials of such low molecular weight content, where the continuous solvent is water, are often referred to as hydrogels. A specific class of anisotropic colloidal particles forming hydrogels are self-assembled peptide fibers, where peptides are the individual building blocks which spontaneously assemble into long one-dimensional fiber-like structures. These fibers may further interact with one another to form higher level assemblies or networks with varying mechanical properties. The biocompatibility and often high tunability make these fibrous networks of great interest for biomedical applications, such as drug delivery, implantable materials or as scaffolds for tissue regeneration [7][8][9][10] .
Fiber forming peptides can be found in nature. They are typical hallmarks of neurodegenerative diseases such as Alzheimer's and Parkison's 11,12 but are also found in many other protein systems 5,[13][14][15][16] . All these fibers appear to contain some sort of laminated β -sheet structure. Inspiration has been taken from these naturally occuring systems and the same or similar structural mo- tifs have been used when designing peptides for controlled selfassembly 3,6,17,18 . It has been found that the mechanical properties of the fibers as well as the macroscopical mechanical properties of their assemblies is strongly correlated to the structural motifs that are present within the fibers 19 .
Most peptide hydrogel systems consist of a network of many micrometer long fibers having contour lengths, L, that are much longer than their persistence length l p , which in turn can vary from tens of nanometers to tens of micrometers 5,20 . The elasticity of these hydrogel systems can often successfully be described using theory developed for flexible and semi-flexible polymers 21 , with fiber rigidity and network mesh size being the important parameters. However, fiber-fiber interactions also affect the stiffness of the network 1,22,23 . Peptides typically contain both charged and hydrophobic amino acid, giving rise to simultaneous long range electrostatic repulsion and short range attractive interactions. Recent work indicate that increasing the relative strength of attractive interactions leads to a stiffer network due to an increase in fiber-fiber contacts, that decreases the entanglement length in the network 22 .
gregates are crystalline 26 , making them highly rigid, and they therefore offer an interesting model system for studying the behavior of weakly charged colloidal rods that also exhibit a short range attractive interaction due to their hydrophobicity. Tuning the electrostatic repulsion by the addition of electrolyte allows for exploring this peptide hydrogel system from a repulsive glassy state at low salt content, to a heterogeneous gel state at a higher salt content where interactions are predominantly attractive. We here adopt the terminology that a gel is a material with arrested dynamics due to attractive "sticky" interactions, while in a glass the dynamics has been arrested because of repulsive excluded volume interactions 1,27,28 .

Materials
The model peptides A 8 K and A 10 K were purchased from CPC Scientific as a lyophilized powder. The peptides were delivered as trifluoroacetate (TFA) salts at a purity of roughly 98 % and were used without further purification. The peptides have one TFA counterion and can therefore be described for example as (NH 2 −A n K−COOH) + TFA − . Samples were prepared by mixing the amorphous powder in either water or D 2 O, the latter chosen for the possibility to perform nuclear magnetic resonance (NMR) measurements. Sample volume fractions were calculated from mass fractions and previously reported peptide densities 24 .

Dynamic Lights Scattering
Dynamic light scattering was measured on an CGF-8F compact goniometer system from ALV Gmbh, Langen, Germany. The system includes a CW He-Ne gas laser with a wavelength λ = 632.8 nm and an output power of 22 mW. The laser intensity is automatically attenuated and detected by a fiber optical near-monomodal detection system splitting up the signal on two avalanche photodiodes in pseudo-cross correlation arrangement. An ALV-7004 multiple tau digital correlator is used for the time correlation function of the scattered intensity. 8 samples in water in the concentration range 0.00034 < φ < 0.028 were measured at four different angles. The 4 lowest concentration were measured for 300 s while the higher concentrations, showing slower dynamics, were measured for 600 s.
The experiment measures the intensity auto correlation function g (2) (t) from which we construct C(t) = (g (2) Here, β is an instrumental constant, close to unity, g (1) (t) is the correlation function of the electric field, and the second equality corresponds to the so called Siegert relation 29 . C(t) is then defined in the inteval 0 ≤ C(t) ≤ 1. All experiments were performed at 25 • C.

Rheology
Rheology was performed on a stress-controlled Anton Paar Physica MCR 301 rheometer with a Peltier system for temperature control. A cone-plate geometry Anton Paar CP25-1 tool with a diameter of 25 mm and an angle of 1 • was used together with a water vapor solvent trap. All measurements were performed at 25 • C. Gel-like samples were loaded using a spatula whereas more viscous samples were loaded using a syringe. Before the start of each measurement a 20 min time sweep was performed at ω = 1 rads −1 and γ = 1 % to ensure that apparent steady state was reached after loading.
Oscillatory measurements were performed at 1 % strain in the linear viscoelastic regime, determined from amplitude sweeps at a frequency of 1 rads −1 . Steady shear flow curves were acquired by step-wise increasing the shear rate and awaiting steady state before the measurement point was collected. To determine the recovery time of the gel-like material after high shear a shear step measurement was performed where a half hour time sweep at γ = 1 % and ω = 1 rads −1 was followed by a 1 minute shear step atγ = 1000 s −1 and a 1 h time sweep with the same parameters as previously mentioned.
All samples used for rheological measurements including various salt concentrations were diluted from a concentration φ ≥ 0.035 and further aged approximately 20 h.

Small angle X-ray scattering
Small and wide angle X-ray scattering (SAXS/WAXS) experiments were performed on a Saxslab Ganesha instrument, equipped with JJ X-Ray system Aps pinhole, an X-Ray microsource with a wavelength λ = 1.54 nm, and a movable two-dimensional 300k Pilatus detector (Dectris Ltd., Switzerland). Three detector sampleto-detector distances were used for data collection and the azimuthally averaged scattering intensities, I(q), as a function of the scattering vector q = (4π/λ ) sin θ /2 where θ is the angle between the scattered beam and the incoming beam. Measured scattering intensities were put to absolute scale using a built-in calibration procedure of the instrument. All experiments were performed at 25 • C.

Ribbon phase behavior
A 8 K and A 10 K form twisted ribbon aggregates of laminated βsheets when dissolved in water. These ribbons have a crosssection of circa 4 × 8 nm 2 , and are typically about 100 nm long for A 10 K, and roughly 180 nm, for A 8 K 24,25 . The cross-section is strikingly monodisperse, a fact that can be explained as a compromise between hydrophobic interactions and a torsional deformation that stretches the β -sheet hydrogen bonds 25,30 . The ribbons form at a critical aggregation concentration, φ cac , which is 0.001 for A 8 K and 5 × 10 −6 for A 10 K, measured in volume fraction 24 . Due to the fact that the ribbons are crystalline [24][25][26] , they are also highly rigid, and can be viewed as charged rigid rods. At lower concentrations, the solution phase is isotropic. However, when increasing the concentration a transition to an anisotropic phase, presumably nematic, occurs. Fig. 1 shows a number of samples of A 10 K in water, viewed in transmitted light between crossed polarizers. A transition from an isotropic solution to an anisotropic birefringent phase is observed as the volume fraction, φ , is increased above approximately 0.02.
Surprisingly the same isotropic-to-nematic phase transition is estimated to φ IN ≈ 0.1 in the case of A 8 K, in spite of the fact that the A 10 K ribbons are shorter. These values can be compared with Onsager's predicition for the hard rods, where the isotropic-tonematic transition occurs at φ IN ≈ 4d/L 31 , d being the rod diameter and L its length. For A 8 K the average L/d ≈ 180/6 ≈ 30 from which we expect the transition at circa φ ≈ 0.13 only slightly higher than the observed value. For A 10 K, which has a smaller aspect ratio L/d ≈ 100/6 ≈ 15, the Onsager model predicts the transition to occur at a higher concentration φ IN ≈ 0.25, while we experimentally observe the transition already at φ IN ≈ 0.02.
This can be understood from the fact that the ribbons carry a net charge. Due to the electrostatic interparticle repulsion, the rods will act as if they have an effectively larger diameter, implying a decrease in the effective aspect ratio [31][32][33] . In addition, the electrostatic interaction destabilizes the nematic phase 32 . For two rods arranged side by side, the electrostatic potential is highest when they are parallel and lowest when the angle between them is 90 • . The increased effective diameter however also corresponds to a larger effective volume fraction and the overall effect of the electrostatic interaction is that the isotropic-to-nematic transition moves to a lower concentration 32 .
The situation regarding electrostatic interactions is here however somewhat complicated. The peptides have one TFA counterion and can therefore be described for example as (NH 2 -A n K-COOH) + TFA − . When the peptides self-assemble, a certain fraction (≈ 50 %) of the peptide ions dissociate a proton and become electroneutral 34 . When doing so, the solution becomes acidic, as trifluoroacetic acid is a relatively strong acid, pK a ≈ 1 35 . This implies that the ionic strength of the solution varies with the peptide concentration, as [H + ]≈ 0.5 (c − c s ), where c is the peptide molar concentration and c s is the monomer solubility. The pH reaches ≈ 2 at φ ≈ 0.028. In addition, there is the screening from the peptide monomers. In the case of A 10 K, c s = 7 µM 24 , corresponding to φ cac ≈ 5 × 10 −6 , which is very low. c s has however been found found to be an order of magnitude higher because of unaggregated monomers 36 . For A 8 K on the other hand, the solubility is then almost one order of magnitude higher, c s = 2 mM 24 . Thus, there is a higher ionic strength in the A 8 K system compared to A 10 K, resulting in a shift of the isotropic-to-nematic transition to higher concentrations.
It is further important to keep in mind that the background screening electrolyte concentration, due to [H + ] increases with c. Thus, electrostatic interactions become increasingly screened the higher the peptide concentration. Because of this, and because the aggregate lengths are polydisperse, the isotropic-nematic two phase region is expected to be wide. We refrain from a more detailed analysis of the electrostatics, but conclude that with A 8 K the interactions are more screened compared to A 10 K. We also note that the dynamics at higher concentrations are very slow (vide infra), and it is very difficult to establish an equilibrium state. The system can be trapped in a more or less isotropic glass.

Arrested diffusion
Rod-like particles in solution have translational and rotational degrees of freedom. In dilute conditions, diffusion is unhindered. However, as the concentration is increased, in particular above the overlap concentration, φ * ≈ 6 (d/L) 2 , certain degrees of freedom become significantly hindered. For translational diffusion we may consider diffusion parallel and perpendicular to the rod axis, as statistically independent diffusion modes, with diffusion coefficients D and D ⊥ , respectively. At infinite dilution the friction coefficients of these two diffusion modes only differ by a factor of 2 (D /D ⊥ = 2). Calculations of the diffusion coefficients were first made by Broersma [37][38][39][40] . Slightly different values of some coefficients were later proposed by Tirado et al. 41 . Following Broersma, we consider and Here, η S is the solvent viscosity, k B is Boltzmann's constant, T is the absolute temperature and the functions δ , γ and γ ⊥ are defined as In Figure 2A we present a normalized correlation function C(t), obtained from DLS experiments, for different concentrations of A 10 K. For the lowest concentration, the separation of modes is not clearly seen and the measured correlation function is consistent with only a single diffusion mode being only an average of the parallel and perpendicular modes. In this case, the diffusion coefficient is given by 42 From a second order cumulant analysis of the φ = 3.6 × 10 −4 data we estimate D t ≈ 7 × 10 −12 m 2 s −1 . Using the viscosity of water η S = 0.89 mPas, T = 298 K and setting the effective diameter d = 6 nm we obtain L = 300 nm, and a polydispersity index of 0.3. We note that L here corresponds to a Z-average, where larger particles in the ensemble are weighted more.
With increasing concentration, we see the appearance of a slow mode, clearly visible for φ > 0.01. At the highest concentration, φ = 0.028, the correlation function still decays to zero, but very slowly requiring several seconds. At these highest concentrations, the sampling time was 600 s. This indicates that the system approaches an arrested or jammed state. At the highest concentrations, φ > 0.02, the intercept, g (2) (0) − 1, shown in the insert in Fig. 2A, decreases. This further indicates that the system is becoming non-ergodic 43 . All modes are diffusive as relaxation rates are proportional to q 2 , Fig. S1, q being the magnitude of the scattering vector. Some deviation is seen at high concentration and low q, but this is likely due to too short measurement times when the decay time becomes very long. We see no sign of rotational diffusion in the present data. But this is not expected as qL < 1. Attempts to detect the rotational diffusion using depolarized (vertical-horizontal, VH) light scattering were made, but the signal was too weak. Jamming and dynamic arrest have been studied extensively with spherical colloids 44,45 . Similar bimodal decays are observed in the vicinity of a glass transition, that are interpreted in terms of α and β structural relaxation, respectively. The fast α relaxation corresponds to the rattling motion of particles trapped within the cage of nearest neighboring particles. The slow β -relaxation corresponds to the occasional opening of the cage due to the collective motion of the neighboring particles.
There are, however, some particular differences. Spheres jam at much higher volume fractions 44,45 . The hard sphere glass transition occurs at φ = 0.58. While this can be considered merely as a quantitative difference, there is also a significant qualitative difference. For the spherical geometry, both the local short time diffusion as well as the long time diffusion is isotropic, and when the system jams, all three translational degrees of freedom become arrested simultaneously. For anisotropic rod-like particles, on the other hand, the perpendicular diffusion is strongly hindered in crowded conditions, while the parallel diffusion may be essentially unhindered. For rod-like particles, rotation, with diffusion coefficient D r , is also important, and like the perpendicular diffusion, the rotation also become strongly hindered above the overlap concentration 46 .
In a previous study 24 , the fast and slow decays were interpreted as being associated with D and D ⊥ , respectively. A similar interpretation was also made recently in a study of cellulose nanocrystals 47 . Another possible source of slow mode could be the self-diffusion, as the particles are highly polydisperse 48 . Selfdiffusion modes have been observed clearly for example for polydisperse microemulsion droplets 49 . Kang and Dhont studied the dynamics of charged fd virus particles (L/d ≈ 100) approaching the glass transition 50,51 . A significant slowing down was observed with increasing concentration in the liquid state, followed by an abrupt vitrification when crossing the glass transition. In their experiments, however, no clear time-scale separation of two modes were observed.
To summarize, we conclude that the A 10 K dispersions approach a state of non-ergodicity and dynamic arrest at φ ≈ 0.03. SAXS experiments have shown that interactions are significantly repulsive 25 , and the arrested state should be a repulsive glass. There is still significant motion in the glassy state as seen by the apparent β -relaxation. Interestingly, the systems shows a non-ergodic behavior even though the particles are essentially free to diffuse parallel to their axis. Perpendicular diffusion and rotation, however, are expected to be strongly hindered.
In Fig. 3 we show the linear viscoelastic spectrum obtained for φ = 0.028 of A 10 K ribbons. We observe that G G , where G is the storage modulus and G the loss modulus. This is indicative of a system that behaves as an elastic solid rather than a viscous liquid. Furthermore, but both the moduli are independent of the deformation frequency, ω, within the frequency window of the experiment, covering three orders of magnitude. The frequency independence implies that τ ω −1 min = 100 s. This should be compared to the time it takes to diffuse its contour length, L 2 /D ≈ 10 −5 s, where we have used L ≈ 100 nm and D ≈ 10 −11 m 2 s −1 where D is taken from the data of Fig. 2A. A striking difference.
When an isotropic system of homogeneously distributed rigid rods is subjected to an applied shear stress, an anisotropic distribution of rod orientations is induced. Relaxation of this anisotropy, back to the isotropic state is done by particle rota- tion alone, as the translational degrees of freedom, parallel and perpendicular to the rod axis, do not change orientation 52 . The characteristic stress relaxation time is therefore τ ∼ 1/D r . This clearly differentiates rigid rod systems from systems of flexible chains with internal degrees of freedom. In melts or solutions of entangled flexible chains the stress is typically relaxed by reptation 53 , where chains diffuse along their contour (parallel diffusion) within a tube defined by entanglements with neighboring chains. Typically, flexible chains need only to escape the original tube in order to completely relax the stress and the induced anisotropy and the τ is essentially the time it takes for a tube renewal, referred to as the reptation time, τ rep ≈ L 2 /D rep , where D rep is the reptation diffusion coefficient, similar to D discussed above. Rigid rods, on the other hand, must diffuse many times their own contour length in order to relax the stress, as for each tube renewal there is only a minor rotation, depending on the concentration (crowding). For this reason, very long τ, due to very low D r , are obtained even for relatively short rods. A striking difference compared to flexible chains where corresponding glassy states are not obtained in semi-dilute solutions.

Linear Viscoelastic Regime
The onset of the dynamical arrest can also be studied by probing the viscoelastic properties of the system. Measurements of A 8 K and A 10 K were performed with small amplitude oscillatory shear (SAOS) rheology in the linear viscoelastic regime, where no plastic deformation occurs. By studying the dynamic moduli G and G as a function of concentration we obtain information about when the system transforms from a viscous liquid to an elastic solid 54,55 . In viscous liquids the input energy is dissipated as heat and G < G . Elastic solids on the other hand, where the energy is stored shows G > G . At the transition between the two regimes, G ≈ G , the characteristic viscoelastic relaxation time is typically found. In Fig. 4 values of G are obtained at 1 % strain and 1 rads −1 are shown versus φ . At low enough concentration, roughly φ < 0.005, the A 8 K and A 10 K aggregate dispersions show G ≈ G , however the absolute values of G are noisy and close to the instrumental limit. Therefore, we classify the system in this concentration regime as a viscous liquid. At φ ≈ 0.005 the response starts to become predominantly elastic (G > G ), and we identify this as the overlap concentration φ * ≈ (d/L) 2 . Here, the system is till ergodic and we do not yet see any significant slow mode in DLS, Fig. 2A. Increasing the concentration further, the ratio G /G increases. For φ > 0.02, we have G G and the system is also becoming non-ergodic and show the presence of very slow dynamics in DLS, Fig 2A. φ ≈ 0.02 also corresponds with the visually observed liquid to gel transition. At the transition concentration (φ = 0.02) we further observe the proportionality G ∝ G ∝ ω (∼0.5) , reminiscent of the criteria of the Winter-Chambon principle for when a system transitions into a fully percolated network 56 . The viscoelastic relaxation time in this high concentrated regime, φ > 0.02, becomes practically infinite and G is above 10 3 Pa, not uncommon for peptide hydrogels 15,57 . Representative amplitude and frequency sweeps from this regime can be found in Figs. 5C and D.
An interesting observation is the particularly strong concentration dependence G ∝ φ n . A power of n = 4.8 was fitted to the high concentration data and is shown as a dashed line in Fig. 4. Models describing the scaling behavior of the plateau modulus as a function of concentration for semi-flexible and rigid polymers and rod dispersions show a power-law behavior where n ≈ 2 − 2.5 depending on fiber properties 21,58 . These models have also been shown to hold for other fibril forming peptide model systems, such as the MAX1 system 3 and the RADA16-II 4 , for amyloid fibrils formed from proteins such as β -lactoglobulin 16 and lysozyme 14 and for other biologically relevant fibrillar systems such as actin 59 .
One of few reports of a peptide system deviating from the previously mentioned theoretical scaling behavior is the octapeptide FEFEFKFK 60 , where a scaling behavior of G ∝ C 3.7 was found for systems at pH= 4 whereas higher and lower pH showed a lower power. When the peptides fibrils became less charged at a specific pH, they could laterally aggregate and become more rigid. This would lead to a higher concentration dependence of G .
A strong concentration dependence of the storage modulus has also been found in carboxylated nano-crystalline cellulose (NCC) systems 61,62 where scaling behaviors with n ≥ 3 are found. The NCC systems studied are very similar to the peptide aggregates found in the A 8 K and the A 10 K system. They are not only charged stabilized and show a similar aspect ratio, but are also highly crystalline and thus very rigid. Inspired by the cellulose work, 61 , Hill proposed a model for the storage modulus of rigid rod networks, allowing for the effects of repulsive or attractive interparticle interacions 23 . In this model G is given by where E is the Young's modulus of the rods, β is an order-one constant, and interactions are accounted for by the parameter α.
Attractive interactions lead to an increased density of rod contacts in the density of rod contacts, and α > 0. For repulsive interactions, Hill considered α < 0. A similar idea with attractive ("sticky") interactions results in an increased contact density and an increase in G was recently proposed by Semerdzhiev et al. 22 The mechanical properties of β -sheet containing fibers have a strong correlation to the density of hydrogen bonds in the crosssection of the fiber. The more regular and tightly packed the hydrogen bonds are, the higher E the fiber will have 63 . A rough method of determining the E has been given by Knowles et al. 19 , where the E is estimated by a product of the hydrogen bond lengths of within the β -sheet, δ β , the hydrogen bond density, ζ , and a hydrogen bond spring constant, k H . The parameters regarding the geometry of the hydrogen bond have been investigated previously and have been determined to δ β = 0.45 nm and ζ = 1/(0.51 × 0.35) nm −2 . Using k H = 13 Nm −1 , as suggested by Knowles et al. 19 , results in E ≈ 32 GPa for both the A n K peptides in this work. This value is similar to the Young's modulus estimated for spider silk 64 , but larger when compared to amyloid fibrils where different values of Young's mudulus have been reported (0.3 − 4 GPa) 65,66 . We also note that an estimate of E for the A 8 K and A 10 K fibres is higher than that for NCC used as the example system by Hill (5.5 GPa) 23,61 .
Finally we note that from a simple fit of G = Aφ 4.8 the high concentration data in Fig. 4 we obtain A ≈ 30 GPa, which is similar to the estimated value of the Young's modulus presented above. This is consistent with the ideas of Hill 23 and Equation 4 that the dominating local fibril deformation is elongation.
As seen in Figs. 4 and 5, the rheological response is essentially the same for A 8 K and A 10 K for a given φ . This in spite that A 10 K samples are birefringent for φ > 0.02. A 8 K samples showed no birefringence within the concentration range studied using rheology. This shows that the effect of ordering in A 10 K is small. Possibly, because we are dealing with a two phase region where the amount of ordered domains is small.

Repulsive and Attractive Fibril Interactions
As suggested above, the A 8 K and A 10 K aggregates are hydrophobic, which influences the rheological behavior. In pure water they are sufficiently charged to provide colloidal stability 24,25 , but when screening the long range electrostatic interaction by the addition of salt, the short range attractive interactions become predominant and the rods aggregate into clusters. The samples then strongly scatter light and become turbid. With pure water, the samples are highly transparent. This change is visualized in Fig. 6A comparing two samples of A 10 K φ ≈ 0.025 in pure water and 0.1 M NaCl. It is also possible to observe this change in interactions upon the addition of salt from small angle X-ray scattering (SAXS) measurements. In Fig. 6B we present SAXS patterns of A 10 K φ = 0.025 recorded in different aqueous concentrations of NaCl, varying from 0 to 0.2 M.
The scattered intensity in a small angle scattering experiment, I(q), is proportional to a product of the average single particle scattering function, often referred to as the form factor P(q) , and the effective structure factor, S eff(q) , that contains information on solution structure and hence interparticle interactions. We use Fig. 6 (A) Comparison of the turbidity of a φ ≈ 0.025 A 10 K sample in pure water, and in an aqueous salt solution of 0.1 M NaCl. In pure water the system is colloidally stable due to repulsive interactions. Upon the addition of salt peptide fibers cluster resulting in a turbid sample. In (B) the same A 10 K concentration is analyzed with SAXS at different aqueous concentrations of NaCl. No change in the observed scattering curve at q < 0.07 Å −1 shows that the individual peptide fibrils do not change. An increase is observed in the low q region which is shown closer in the insert indicates a transition from repulsive to attractive interactions. The scattering pattern at high salt concentration shows a power law scattering of I(q) ∼ q −1.5 which could possibly be interpreted as a fractal dimension of the fiber clusters.
the word effective here, because a system with size polydispersity should strictly be considered as a multicomponent system where all different cross correlations should be taken into account as a sum of partial structure factors. Using this approximation we get 67 Here ∆ρ is the contrast, i.e. the difference in scattering length density between particles and solvent, V p is the particle volume and the brackets ... denote an ensemble average. In the so called thermodynamic limit of q = 0, the effective structure factor can be interpreted in terms of an effective osmotic compressibility, S eff (0) = k B T V p (∂ π/∂ φ ) −1 , where π denotes the osmotic pressure. For repulsive interparticle interactions S eff (0) < 1, and the effect of the repulsive interaction is to decrease the scattering intensity at smaller q-values. If interactions, on the other hand are predominantly attractive, S eff (0) > 1, resulting in an increased scattering at smaller q-values. S eff (0) = 1 corresponds to an ideal solution.
The A 8 K and A 10 K systems have been quantitatively character-ized in water by SAXS previously, both with respect to P(q) and S eff (q) 24,25 . From those experiments it was concluded that interactions are repulsive (S eff (0) < 1). From P(q) it was concluded that the fibril cross section was biaxial and that it can be described as elliptical with semi-axes of 1.8 nm and 4 nm 24,25 . With this information at hand we compare the scattering curves obtained at different NaCl concentrations, that are shown in Fig. 6B. The data cover a q-range of roughly 3 orders of magnitude including both the small angle regime with information of the aggregate morphology and interactions as well as the wide angle regime (high q) where the diffraction peaks show that the aggregates are crystalline. The diffraction peaks can be indexed to a two dimensional oblique unit cell 26 .
First of all we note that the different scattering curves completely superimpose for q > 0.07 Å −1 . From this we can conclude that the individual fibrils are not affected by the addition of salt. Changes in I(q) are only observed at lower q, which we interpret as resulting from changes in S eff (q). The low q scattering is highlighted in the insert of Fig. 6B. The addition of salt leads to an increase in the small q scattering, indicating a crossover to attractive interactions. We do not have access to sufficiently low q-values to determine S eff (0). For the highest salt concentrations, 0.1 and 0.2 M, the small angle scattering shows a power law with I(q) ∼ q −1.5 . These samples are also very turbid in visible light suggesting that the attractive interactions have become predominant and the fibrils have aggregated further into clusters, resulting in a heterogeneous fibrillar network. The power −1.5 can possibly be interpreted as a fractal dimension of the fiber clusters 1,68 .
Also the non-linear rheology of a φ = 0.025 A 10 K sample depends on the electrolyte concentration. In Fig. 7 we present the variation of the stress amplitude, σ 0 , as a function of the strain amplitude, γ 0 , in oscillatory shear experiments at the different NaCl concentrations. The oscillation frequency used was ω = 1 rads −1 . In the linear regime, at lower γ, we can extract the dynamic shear modulus |G * | = σ 0 /γ 0 . The variation of the modulus with the NaCl concentration is plotted as an insert in Fig. 7. |G * | first increases with increasing salt concentration and then essentially levels off. A similar value is obtained at 0.1 and 0.2 M salt. This value is approximately an order of magnitude larger compared to the value at zero salt. The increase in |G * | with increasing salt concentration is understood as an increase in the entanglement density due to the increase of hydrophobic contacts when the electrostatic repulsion gets screened 22,69 . A similar increase in network stiffness is also observed in α-Synuclein hydrogels when decreasing the fibril charge density by pH 70 . Moreover, similar salt effects are observed for nanocrystalline cellulose 69 .
Frequency sweeps were also performed with the samples containing salt, and the results are presented in Fig. S2. The samples with salt showed the same features as the sample without salt, Fig. 3, with essentially no frequency dependence within the frequency window. However, G increases with increasing salt concentration similar to the behavior of G * , shown as an insert in Fig. 7.
We also see in Fig. 7 that the addition of salt leads to the build up of a stress overshoot at γ 0 ≈ 10 %. This is consistent with Fig. 7 The stress amplitude, σ 0 , as a function of the strain amplitude, γ 0 , and salt concentration for a φ = 0.025 A 10 K sample and at an oscillation frequency ω = 1 rads −1 . Upon the addition of salt a yield stress is built up due to the transition from repulsive to attractive interactions. The insert shows |G * | versus the salt concentration as determined from the σ 0 /γ 0 in the linear viscoelastic regime. |G * | increases roughly an order of magnitude with the salt concentration and eventually levels off. The line in the insert is added as a guide to the eye. the transition from predominantly repulsive to predominantly attractive interactions. In the attractive gel, an additional stress is needed to overcome the attractive hydrophobic force and pull particles apart inducing a flow 55 .
As has been shown above, we observe a transition from a homogeneous network in water to a more heterogeneous state with voids in the overlapping network giving rise to an increased turbidity. We can understand this as follows. When screening the electrostatic interactions, the short range attraction becomes the dominating interaction and the system starts to undergo a liquidliquid (i.e. "liquid-gas") phase separation where the rods get concentrated in one phase coexisting with a dilute or excess solvent phase. However, a full macroscopic phase separation is not obtained. The process gets at some point arrested resulting in the kinetically trapped gel state 28 .

Flow properties
Knowing how peptide hydrogels respond to shear is important to correctly assess future applicability. Materials that are injected by syringe are typically subjected to extremely high shear rates, which can have a high impact on the material. By rheological flow curves and shear rate step measurements it is possible to a certain degree to simulate the process of injection through a needle 71,72 . In Fig. 8A we present the viscosity η as a function of shear rateγ for two samples, A 8 K φ = 0.033 and A 10 K φ = 0.039 respectively. As can be seen, both samples show a strong shearthinning behavior, with essentially η ∝γ −1 over the whole range of shear rates investigated. This corresponds to the maximum possible shear thinning and implies a constant shear stress, σ , within the range ofγ. These types of systems are shear thinning because the rod-like aggregates align in the shear flow 26,73 . A question that remains is whether this alignment is continuous, with an order parameter describing the degree of alignment, that is homogeneous throughout the gap, or whether the alignment is discontinuous. The latter case implies shear banding in the gradient direction 74 , where a fully aligned nematic state, N, of relatively low viscosity, coexisting with a highly viscous isotropic network state, I. The present experiments where performed on a stress controlled rheometer but in a simulated rate controlled mode. In this case the setγ is a weighted average (lever rule),γ = (1 − f N )γ I + f NγN , where f N is the fraction of aligned nematic state in the gap, andγ I andγ N , withγ γ N , are the shear rates of the isotropic and aligned nematic states, respectively. Withγ I andγ N being fixed, f N increases from 0 to 1 in the rangeγ I <γ <γ N . Such a scenario with a shear induced nematic state has been found and studied extensively in wormlike micelle systems [75][76][77] . In the present system, it is known that the aggregates align in shear flow 26 , but the detailed nature of the alignment remains to be investigated.
Shearing the samples at a high shear rate result in a low viscosity, presumably an aligned (nematic) state. After cessation of flow, the sample is expected to relax back to its initial state of a stiff network. The kinetics of this recovery is of interest for some hydrogel applications, for example for injected subcutaneous drug deposition. In Fig. 8B we have plotted the elastic and storage modulus obtained from an oscillatory experiment as a function of time after pre-shearing the samples atγ = 1000 s −1 for 60 s. The elastic modulus measured directly before the pre-shear is shown as the horizontal lines. Roughly 75 − 80 % of the elastic modulus is recovered within 1 h. This reversibility is also shown by the fact that the flow curve in Fig. 8A for A 10 K φ = 0.039 is near to identical both asγ is increased and whenγ is decreased. Interestingly, the moduli are G ≈ G ≈ 0 just after the pre-shear (t = 0) indicating that the network is completely disentangled and a fully aligned nematic state has been achieved.
Shear-thinning as a response to increased shear rate is, as previously stated, typical behavior for rod dispersions as the rods align along the direction of the flow. This behavior has recently been explained thoroughly for various viruses 73 . Many studied peptide fibril systems however have long fibres where alignment is not possible due to the large number of entanglement points. For the MAX1 peptide system the shear thinning instead occurs by the fracturing of the network and not of alignment of fibres 71 . Due to the shorter aggregate lengths in these systems, alignment has been shown to take place 26 , this does however not exclude network fracturing.

Conclusions
Self-assembled peptide hydrogels are interesting biomaterials with a range of potential applications [7][8][9][10] . Typically, such systems consist of micrometer long fibrils forming a transient network 3,6,15 . Here, we have shown that peptide hydrogels can also be obtained with relatively short, 100−200 nm fibrils, corresponding to an aspect ratio L/d of roughly 20. These peptide aggregates can be viewed as rigid rods with a sticky hydrophobic surface that are still colloidally stable due to a net charge. Thus in water, these hydrophobic but charged aggregates interact through a combination of a long range repulsion and a short range attraction like in other similar systems 22,69,78 , and where both interactions are significant. There are strong similarities to the classical DLVO potential although here the attraction is mainly due to sticky hydrophobic interactions rather than van der Waals interactions.
With increasing concentration, above the overlap concentration φ * , the rotational diffusion becomes increasingly hindered and D r gradually approaches zero. In water, the repulsive electrostatic interaction still provides colloidal stability, and we identify the system as a repulsive glass, where the rotational degree of freedom has been arrested. In the present system, the short range attraction still seems to be significant, increasing the density of effective aggregate contacts leading to an increased network stiff-ness and a strong concentration dependence of the elastic modulus (G ∼ φ 4 .8) 23 . When screening the electrostatic repulsion by the addition of salt, the attractive hydrophobic interactions becomes predominant and the system crosses over from a repulsive glass to a gel state. A similar behavior has previously been reported for nanocrystalline cellulose having a similar aspect ratio as the peptide aggregates here 69 .

Conflicts of interest
There are no conflicts to declare.