Volume 12 Supplement 10

## Proceedings of the Eighth Annual MCBIOS Conference. Computational Biology and Bioinformatics for a New Decade

# A modified Stokes-Einstein equation for A*β* aggregation

- Srisairam Achuthan
^{1}, - Bong Jae Chung
^{2}, - Preetam Ghosh
^{3}, - Vijayaraghavan Rangachari
^{4}and - Ashwin Vaidya
^{5}Email author

**12(Suppl 10)**:S13

https://doi.org/10.1186/1471-2105-12-S10-S13

© Achuthan et al; licensee BioMed Central Ltd. 2011

**Published: **18 October 2011

## Abstract

### Background

In all amyloid diseases, protein aggregates have been implicated fully or partly, in the etiology of the disease. Due to their significance in human pathologies, there have been unprecedented efforts towards physiochemical understanding of aggregation and amyloid formation over the last two decades. An important relation from which hydrodynamic radii of the aggregate is routinely measured is the classic Stokes-Einstein equation. Here, we report a modification in the classical Stokes-Einstein equation using a mixture theory approach, in order to accommodate the changes in viscosity of the solvent due to the changes in solute size and shape, to implement a more realistic model for A*β* aggregation involved in Alzheimer’s disease. Specifically, we have focused on validating this model in protofibrill lateral association reactions along the aggregation pathway, which has been experimentally well characterized.

### Results

The modified Stokes-Einstein equation incorporates an effective viscosity for the mixture consisting of the macromolecules and solvent where the lateral association reaction occurs. This effective viscosity is modeled as a function of the volume fractions of the different species of molecules. The novelty of our model is that in addition to the volume fractions, it incorporates previously published reports on the dimensions of the protofibrils and their aggregates to formulate a more appropriate shape rather than mere spheres. The net result is that the diffusion coefficient which is inversely proportional to the viscosity of the system is now dependent on the concentration of the different molecules as well as their proper shapes. Comparison with experiments for variations in diffusion coefficients over time reveals very similar trends.

### Conclusions

We argue that the standard Stokes-Einstein’s equation is insufficient to understand the temporal variations in diffusion when trying to understand the aggregation behavior of A*β* 42 proteins. Our modifications also involve inclusion of improved shape factors of molecules and more appropriate viscosities. The modification we are reporting is not only useful in A*β* aggregation but also will be important for accurate measurements in all protein aggregation systems.

## Introduction

Aberrant misfolding and aggregation of proteins have been implicated in over 40 different human pathologies including Alzheimer’s disease (AD), Parkinson’s disease, type 2 diabetes, transmissible spongiform encephalopathies (TSE) or Prion diseases, Huntington’s disease and Cruedfelt-Jacob disease (CJD). A biochemical commonality in these diseases is that the protein involved forms pathogenic aggregates, irrespective of whether the monomeric protein is in the misfolded or intrinsically disordered forms, that have a consensus structural moiety, commonly known as “amyloids”. In all amyloid diseases, such protein aggregates have been implicated fully or partly, in the etiology of the disease. Due to their significance in human pathologies, there have been unprecedented efforts towards physiochemical understanding of aggregation and amyloid formation over the last two decades. It is well understood that the process of aggregation towards amyloid fibrils from a monomeric state is a nucleation-dependant mechanism, analogous to crystal growth. In such a process, it is believed that there is a conformational change during the nucleation event followed by rapid aggregation to a state that forms large insoluble or soluble aggregates, which in turn is accompanied by a structural convergence to cross *β*-sheet conformation [1].

In Alzheimer’s disease (AD), a protein called amyloid-*β* (A*β*) peptide forms aggregates that deposit as senile plaques in brains. The nucleation-dependent process of A*β* aggregation is inferred by the occurrence of a ’lag-phase’ prior to fibril growth that shows a sigmoidal pattern [2]. In the aggregation pathway, one important intermediate called “protofibrils” has been isolated and characterized by several groups [3][4][5][6][7][8]. Protofibrils mainly differ from fibrils in their size and solubility; while fibrils can be sedimented with relatively smaller forces (19000g, 10min)-where *g* is the acceleration due to gravity- protofibrils require substantially high sedimentation forces and have smaller diameters than the fibrils [9]. Protofibrils have propensities to both elongate (by monomer addition) as well as to laterally associate (protofibril-protofibril association) to grow into mature fibrils and have been well characterized [9]. These two mechanisms of fibril growth depend on the structure and stability of protofibrils, which in turn depend on the factors affecting nucleation.

*in vitro*A

*β*aggregation process is well known to be affected by several environmental factors such as concentration, pH, ionic strength and temperature. In addition, differential effects of solvents on the aggregation process and size has also been documented [10]. Several biophysical methods are used to monitor the aggregation process among which, light scattering and centrifugation techniques are widely used to measure the hydrodynamic radii (a) as well as the diffusion coefficient (D) of the protein sample [5][11][12][13][14]. For spherical particles, the relation between

*a*and

*D*is given by the Stokes-Einstein (SE) equation;

*k*

_{ B }is Boltzmann’s constant,

*T*is the absolute temperature and

*η*

_{0}is solution viscosity. However, the aggregates of A

*β*peptide cannot be assumed to be spherical and hence determination of

*a*can be erroneous. This difficulty is overcome by the use of multi-angle light scattering (MALS) technique which measures the intensity of the scattered light from multiple angles that in turn helps to calculate the molecular weight without assuming the shape of the solute [9]. A more sophisticated method of measuring the molecular size is by analytical ultracentrifugation (AUC), which utilizes the rate of the movement of sedimentation boundary [15]. For a spherical solute with a molecular mass,

*MA*and radius

*a*, the value of sedimentation coefficient measured at temperature

*T*in buffer

*b*, the rate of migration per unit field,

*s*

_{ A }is given by;

where *N* is the Avogadro’s number and
is the partial specific volume of the solute, and *ρ*, *η* are the density and viscosity of the buffer medium respectively.

An important aspect of the Stokes-Einstein’s equation and its utility in protein aggregation systems in particular is that the viscosity of the bulk solvent in a given buffer solution is considered to be a constant.

*β*aggregation. Specifically, we have focused on validating this model in protofibril lateral association reaction (Figure 1) along the aggregation pathway, which has been experimentally well characterized [9, 16]. The modification we are reporting is not only useful in A

*β*aggregation but also will be important for accurate measurement in all protein aggregation systems.

*η*

_{0}is written as equation (1) where the denominator is the effective drag coefficient for the molecule which is treated as a rigid sphere of effective radius

*a*. We point to several shortcomings of this model with respect to the A

*β*aggregation and suggest modifications to improve it which include:

- (i)
firstly, molecules of concern to us are not necessarily spheres [17] although this is the simplest assumption one can make. In general for our problem, where we consider the aggregation of molecules, several classes of molecules (or n-mers) are assumed to populate the solvent. While it is not clear what the specific shape of each class (n-mer) is, it is most certainly not a sphere [18] for the most part.

- (ii)
The second point we make is that the SE relation takes the viscosity to be a constant, based upon the solvent viscosity. However we argue that it should depend upon the concentration of the various species of molecules present in the system. Any single molecule (individual or aggregated), during diffusion will experience the viscous forces of the liquid as well as due to the remaining molecules around it. As a result, the diffusion coefficient,

*k*, is also a function of the concentration of the molecules. In fact, more generally, we will consider an effective diffusion coefficient*k*_{ e }where*k*_{ e }=*k*_{ e }(*ϕ*_{ i },*η*_{0}) where*i*= 1, 2, …,*n*.

In this paper, we have considered these two important conditions and incorporated modifications to the Stokes-Einstein equation that may faithfully represent the changes occurring in amyloid aggregation systems, especially A*β*. The significance of non-sphericity in computing the diffusion coefficient has been realized in various contexts and the biophysics, chemistry, geological and chemical engineering literature is replete with discussions of the role of non-spherical particles, specifically spheroidal particles [19–23]. While the biophysics community has been well aware of the importance of non-sphericity, the importance of concentration upon the viscosity and hence diffusivity, has been largely neglected. The chemical engineering literature is however a very rich source of discussion on this matter; see [24] and references cited therein for more information on this subject.

## Results

### Theoretical models

*n*-species of molecules, we define

*V*

_{0}to be the volume of the solvent liquid and

*V*

_{ i },

*i*= 1, 2,…,

*n*to be the volume of the

*i-th*species. Then the volume fraction of each element of this mixture is given by:

*ϕ*

_{ i }and must be of the form:

*j*= 1, 2,…,

*n*,

*f*

_{ i }refers to the drag coefficient of the

*i-th*species and

*β*

_{ i }, its respective shape coefficient. As is seen from the last term on the right hand side of the equation, the volume fraction enters the diffusion coefficient through the viscosity. In fact, in this paper, we suggest that the viscosity term itself be treated as an average of the viscosity due to its constituents as a result of which

*η*

_{ i }(

*ϕ*

_{ j }) is replaced by a scalar, the effective viscosity,

*η*

_{ e }(

*ϕ*

_{ j }) where

*j*= 1, 2,…,

*n*. Hence:

While we are following the vast majority of physics and engineering literature in following this approach of modeling a mixture through its overall material parameters, we are also conscious that it not a rigorous mathematical approach. In his paper on the “viscosities of mixtures” Massoudi [24] provides an overview of the historical approach to the modeling of effective viscosities but more importantly, he points to the mathematically rigorous approach of modeling a mixture by examining the stress tensors. However, this is recognized to be a rather formidable challenge to verify experimentally and very tedious to implement in computations such as these. Even within this approach of modeling viscosities, it has been argued [25] that the effective viscosity of say two liquids would be not just a linear superposition of the two individual viscosities, but also include an interaction term, referred to “mutual viscosity” which again is very difficult to verify experimentally. All these point to the complexity of rigorously modeling even a seemingly simple parameter as viscosity. In the rest of the paper, we consider two specific and successful models applied in the engineering literature to our problem of A*β* aggregation. With regards to the shape of the molecules, there has been much discussed about the validity of spherical-shape assumption of biomolecules( [17, 18] and references therein) and while it is recognized that molecular shapes can be fairly complex and random, it still remains worthwhile to represent these arbitrarily composed molecules using some standard shapes that can be imagined to envelope the molecules in a reasonable manner. Our attempt in particular is to go beyond the assumption of the sphere. To do so we need at our disposal, the drag coefficients corresponding to arbitrarily shaped bodies. This is a very daunting challenge in fluid mechanics and at best relatively simpler shapes such as spheres, oblate and prolate spheroids (of any eccentricity) have been analyzed mathematically [26, 27]. However there is some literature in terms of empirical or computational relations which provide the drag coefficients for other shapes.

### Case 1: a sphere-rod empirical model

We begin with an empirical model. Based upon the work of [28–31] we will consider the modification of the SE law by taking a suspension composing of spheres and rods with the former particles being considerably smaller than the latter such that the ratio *a*_{
r
}/*a*_{
s
} > 20, where *a*_{
r
} is the major axis of the rod and *a*_{
s
} is the radius of the sphere. This particular assumption can also be very meaningful in our case where we consider our sample primarily to be composed of smaller monomers and much larger protofibrils. In particular we follow the specific approach of [31] in obtaining an effective viscosity for such a suspension. In the words of [31] the primary assumption in the computation of viscosity is that:

... *for each fraction of a given particle size*, *the smaller particles in suspension have the same effect as a homogeneous fluid with Newtonian viscosity similar to the effective viscosity of a suspension made up of the fraction of smaller spheres. In other words*, *the smaller suspended particles do not interact with the larger particles and are ’sensed’ by the large particles as part of the continuous suspending fluid*.

*η*

_{rel,r}(

*ϕ*

_{ r }) refers to the viscosity of the rods relative to the continuum of a homogeneously distributed spheres in a solvent, , similarly is the viscosity of the spheres relative to the continuum of the pure solvent in the absence of the rods. Also, Φ

_{0}=

*ϕ*

_{ r }+

*ϕ*

_{ s },

*ϕ*

_{ r }refers to the concentration of the rods,

*ϕ*

_{ s }the concentration of spheres and:

*L*,

*b*refer to the length and diameter respectively, depending on the orientation of the rod, we have the relation for the diffusion coefficient, namely:

*ϕ*

_{ r }< 0.125. A similar relation can be written for

*ϕ*

_{ r }> 0.125 which is given by:

If the net volume fraction of all the solute molecules is held fixed (i.e. Φ), then the diffusion coefficient can be written purely in terms of *ϕ*_{
r
} or *ϕ*_{
s
}.

### Case 2: a mixture theory approach

*η*

_{0}is the viscosity of the background solvent,

*ϕ*is the volume fraction of the suspensions and

*ϕ*

_{ m }is the maximum volume fraction at which the viscosity becomes singular, i.e. the critical limit of concentration beyond which flow ceases. The literature estimates 0.58 <

*ϕ*

_{ m }< 0.69 [36, 37] (As in past literature, let us take the average value, namely

*ϕ*

_{ m }≈ 0.6). In classical mixture theory the effective viscosity of two liquids of different viscosities has been taken to be the weighted average viscosity of the components. In attempting to do the same for suspensions, one must exercise a little caution. In the model proposed here, we adopt a mixture approach by trying to combine the solvent viscosity appropriately with equation (15). We schematically display our recipe for obtaining a suspension with solvent volume

*V*

_{0}and such that the volume of each constituent particle (molecule) is

*V*

_{ i }(

*i*= 1, 2,

*…*,

*n*) in Figure 2, while requiring that equations (3) and (4) be satisfied.

where as in the sphere-rod model treated above, *η*_{
i
}(*i* ≠ 0) where *i* ranges over all cyclic permutations of (1, *…*, *j*, *j* + 1, *…*,*n*), represents the relative viscosity of the suspension containing the *j* th solute in a “solvent” composed of *the liquid of viscosity η*_{0}*and all other solute particles*. Note that the equation (16) is not unique. Other formulations of the effective viscosity are possible. The one provided here seems mathematically and physically reasonable and consistent with some models provided in the literature. The term
represents the relative volume fraction of the *j*’th solute relative to the “solvent” that contains it and is independent for different *j*’s. Therefore, for an n-species suspension, the viscosity relation for the *j* solute in a background containing all other species and the pure solvent can be inferred from the algorithm: *η*_{
j
} *→ η*_{
j
}_{+1}*→* …*η*_{
n
}*→ η*_{1}*→* …*η*_{
j
}_{–1} where *η*_{
a
} *→ η*_{
b
} indicates that *η*_{
a
} is measured relative to *η*_{
b
} in the sense of equation (8). In particular, the cyclic permutations of *i* are taken to ensure evenness in the definition of the volume fractions,
for all *j*’s. This is clearly illustrated in the cases of 2 and 3-species which are examined in detail. We illustrate this by means of a few examples.

#### (i) The case of two species

#### (ii) The case of three species

*k*

_{ e }gives:

where *η*_{
e
} is defined by the equation (16) and
where *i*, *j* ≠ 0 that is, the fluid is discounted. Since the diffusion that we are interested in only pertains to the solute species embedded in the surrounding solvent, we need only consider the volume fraction of each solute with respect to the others in the evaluation of the effective diffusion. In particular let us assume, for instance, that two classes (or species) or molecules populate the solvent; the shape are so selected since the shape factors are only known for a small class of shapes. The first are the molecules at the protofibril stage or beyond which are treated as being spheroidal in shape(diffusion of spheroidal bodies had been considered as far back as 1936 by Perrin [38]) and the second class are those smaller than the protofibril stage which will be treated as spheres, even though Zwanzig and Harrison [18] rightly argue that the shapes of molecules are not correctly described by either of these shapes, the latter being far more realistic than the former. We however feel that this is still an improvement over the current approach of treating all the molecules as hard spheres.

_{0}= 1, the volume fraction of the spheroidal particles to be

*ψ*, then the volume of spherical molecules can be given by 1 –

*ψ*. For the two species case we can write

*ψ*

_{1}=

*ψ*and ψ

_{2}= 1 –

*ψ*. Therefore for this example, we can write the effective viscosity as:

*c*and

*b*refer to the major and minor axis of the spheroid. Hence combining the above results allows us to formulate a specific model for the two-species case where the diffusion coefficient can be given by:

*ψ*= 1 and the effective diffusion coefficient can be written in the form:

*ϕ*

_{0}+

*ϕ*= 1 based on the constraint equation (4). For the one-species suspension case when

*ϕ*/ϕ

_{0}= ∈ << 1, then equation (33) can be approximately written at

*O*(

*∈*) as:

where *ϕ*^{(1)} = *ϕ*. As ∈ *→* 0, the diffusion coefficient approaches the classical Stokes Einstein equation since then *ϕ*_{0}*→* 1.

### Protofibril-protofibril lateral association experiment

*µ*M) that was monitored by dynamic light scattering (DLS) instrument. Protofibrils without salt were used as a negative control. The data is plotted in Figure 4 against diffusion coefficient (D) which is inversely proportional to the hydrodynamic radius of the sample.

### Molecular level simulation of protofibril-protofibril lateral association

Next, we used some simple molecular level reaction models to predict the concentration of protofibrils of different size in the lateral association stage that can correspond to the diffusion coefficient estimates from DLS reported above. First, we assumed that the protofibrils are composed of 1600 monomers based on the previous reports [39, 40]. Here, *F*_{1600} denote the protofibrils (comprising 1600-mers and average length 64 nm) [39, 40]. In general, we will let *F*_{
i
} denote a protofibril with *i* number of *Aβ* molecules binding to it during the lateral association phase.

We present separate reaction models considering different number of species formed during the lateral association reaction. The reactions involving these species were modeled under the assumption of a homogeneous mixture of protofibrils (with different initial concentrations) as used in the in vitro experiments. In particular, we considered a maximum of five different molecular species involved in this phase as follows: *F*_{1600}, *F*_{3200}, *F*_{4800} , *F*_{6400} and *F*_{8000}. We have assumed that lateral association will result in at most a 5-fold increase in size of the protofibrils within the first 1.5 hrs of observation such that the system of reactions will not involve the formation of protofibrils beyond *F*_{8000}. This assumption is purely theoretical although it is experimentally well known that three protofibrils associate together to form mature fibrils. We wanted to consider upto 5 species for the interest of modeling and calculations. The reaction fluxes in this are denoted by *R*(*i*, *j*) corresponding to the reaction between *F*_{
i
}_{*1600} and *F*_{
j
}_{*1600} forming *F*_{(}_{
i
}_{+}_{
j
}_{)*1600} for *i* = 1, …, 5 and *j* = *i*, *i* + 1,…, 5 – *i*. We have also assumed that the initial concentration of protofibrils comprise of only *F*_{1600}’s and that of the other species are zero in the reaction models as only minute levels of associated protofibrils are expected to be present initially in the in vitro system.

In the following, we have considered three case studies to estimate the concentration change of the different molecular species involved in lateral association: 2-species (*F*_{1600}, *F*_{3200}), 3-species (*F*_{1600}, *F*_{3200}, *F*_{4800}) and 5 species (*F*_{1600}, *F*_{3200}, *F*_{4800}, *F*_{6400}, *F*_{8000}). Each reaction was considered reversible with *k*_{+} and *k*_{–} being the forward and backward rate constants respectively. Following our previously validated lateral association stage model [16], we have considered *k*_{+} = 9.0 × 10^{–1}*h*^{–1}*mM*^{–1} and *k*_{–} = 6.0 × 10^{–3}*h*^{–1} respectively. While the 2-species and 3-species models are less realistic, we have seen that protofibrils do not grow beyond *F*_{8000} in the first 1.5 hrs of lateral association, and hence the 5-species model should be the best approximation of the in vitro system under study.

### Case 1: 2-species reaction model

### Case 2: 3-species reaction model

### Case 3: 5-speciesreaction model

After the kinetic schemes are established in the reaction models above, we computed the corresponding reaction fluxes, and hence the differential equations that govern the temporal change of these species can be derived from material balances and reaction kinetics. The reaction flux and differential equations for each of the models reported here have also been shown alongside. The initial concentration of protofibrils ([*F*_{1600}]) is set to 1*µ* M, 3*µ* M and 5*µ* M respectively following the in vitro experiments while the concentrations of the other species are assumed to be zero at the start (i.e., [*F*_{
i
}_{*1600}] = 0, *i* = 2,…, 5). Thus, the system of differential equations for each model is properly defined and were solved using Matlab’s ODE toolbox to estimate the concentration of the different species at various time points between 0 – 90 mins as shown in the previous section.

## Case studies

In this section, we consider some special cases to test our proposed models. We consider the cases of two, three and five species embedded in a solvent, introduced earlier. In these examples, we consider each species to be an ellipsoid since the length and width of the molecules considered here are significantly different so they cannot be correctly represented as spheres. As is indicated in Figure 3(c), associations of these spheroids occur laterally; therefore according to our model, the result of n-associations of spheroids of major axis c and minor axis *b* results in a new spheroid with major axis c and minor axis *n* × *b*. Our objective is to employ our model to estimate the effective diffusion coefficient of various systems and compare them to the experimentally measured values of diffusion.

^{–3}Pa-s. The results of our computations are plotted in Figure 5. The right columns of this plot (Figures 5D, E and F) show the results of the experiments which yield the diffusion coefficient and the concentrations alike. The initial concentration of the protofibril is taken to be 5

*µM*in each of the cases. The black points in the curves refer to the 1600 mers while the red points refer to the final stage in each case, i.e. 3200mers in the 2-species case, 4800 mers in the 3 species and the 8000 mers in the 5 species case. Note that the diffusion coefficient always decreases with the diminishing 1600 mers on the other hand growing with the diminishing concentrations of the bigger, associated molecule. The panels on the left hand side (Figures 5A, B and C) indicate the results of the numerical computations for the same cases as the experiments and indicate very similar trends. For ease of comprehension of the plots the values of alone are shown in the graphs here while those of show very similar features. It is to be noted that the black and red points do not actually intersect; the intersection is an artefact of the two different scales represented on the y-axis. While the increasing and decreasing trends remain consistent with observations the actual value of the diffusion coefficient seems to differ. Our model puts the lower and upper limit values of the diffusion coefficient consistently between 1.7 × 10

^{–12}

*m*

^{2}/s and about 2.4 × 10

^{–12}

*m*

^{2}/

*s*while experimental values are higher, ranging upto four times more that our upper limit at some concentrations. Firstly the large variability in the experimental data when compared to our model can be attributed to the fact that the experimental study is based upon DLS observations which provides a time varying effective radius for the system, corresponding to the production of larger molecules which form as a result of the association. These experimentally evaluated radii which range between 30-80nm are then inserted into the classical Stokes-Einstein formula (equation (1)) to get the diffusion coefficient. The experimental estimates for the effective length scale are substantially lower than the more accurate values used in our model, resulting in the deviation from our values of the diffusion coefficient.

Certain features of the computations merit immediate attention. In particular, the profile of the 1600, 3200 and 8000mers is striking. These include the declining values of the apparent intersection point with increasing number of species and also the slope of the red curves whose declining trend in the graph as plotted, gets more rapid with the larger molecule.

## Discussion and future work

*η*

_{0}) in the Stokes-Einstein formula is the one pertaining to the solvent alone. However if one wants to think of the diffusion of any molecule, the effective medium becomes the solvent as well as the remaining molecules embedded in the system. Therefore the solvent viscosity of the classical SE formula needs to be replaced appropriately with one of an effective viscosity which is no longer just the solvent viscosity alone. In fact the concentration terms are as important if not more than the non-sphericity requirement. It is also important to note that while the apparent form of the modified SE formula is similar to that of the classical SE formula, our model has grave implications in understanding the dynamics of the protofibrils and aggregate molecules since the diffusion coefficient is no longer a constant but a function of the shape factors and concentrations of the different species of molecules that populate the suspension. While the change from a sphere to a spheroid might appear to be a minor point, its full impact can only be realized when one studies the problem at the smaller scale where this change can display its full impact as we plan to do in the future; this paper is focused on the larger scale continuum properties of the system as experimental measurements were hard to achieve at a smaller scale for the protein aggregation problem. To see the complexity of this model, one need only look at the diffusion equation that our model generates. For instance, in the case of n-species of n-mers embedded in a solvent,

where *i* = 1,2,…,*n*. The new *k*_{
e
} renders the equation highly nonlinear. The modified SE equation provided in this paper now allows for a more appropriate coupling of the solute-solvent dynamics. Under the realistic conditions of flow we have a two way coupling; the dynamics of the solute is described by the concentration equation (35), while the dynamics of the solvent would be given by the Navier-Stokes equation. Our next goal is to take on this multiscale approach coupled with a reaction term to understand the temporal and spatial variations in the solute (protofibrils and aggregates) distribution. Due to the contributions of this paper, this would be a very new and rigorous addition to the biophysics literature. The model proposed here, therefore merits serious attention since it provides a systemic treatment of the biophysics involved. The current paper is therefore a first step in this long term approach of understanding the problem of protein aggregation thoroughly and rigorously.

## Conclusion

In this report, we have shown that the use of classical Stokes-Einstein equation may not be accurate for the measurement of hydrodynamic radii and consequently molecular weights of amyloid aggregates. Although it is known that appropriate assumption of the shape of solute molecules is essential for accurate size determination while using Stokes-Einstein equation, this is the first time we have shown that it is incorrect to assume solvent viscosity be a constant during the reaction. This is particularly true for reactions involving A*β* protofibrils, as large size of the solute molecules indicate significant effect on bulk solvent viscosity. We have attempted to rectify the problem by incorporating modifications to the Stokes-Einstein equation (Equation 16)that seem to be valid for the protofibril lateral association reaction. It is noteworthy that the modifications reflecting the viscosity changes can lead to realistic models not only for A*β* aggregation but also for a majority of amyloid aggregation systems.

## Methods

### A*β* protofibril preparation

Synthetic A*β* 42 peptide was obtained from synthesis facility at Mayo Clinic, Rochester, MN as a lyophilized powder. A*β* 42 protofibrils were generated and isolated as previously reported (9). Freshly purified A*β* 42 monomers (100 *µ* M) from size exclusion chromatography buffered in 10 mM Tris, 50 mM NaCl, pH 8.0 was agitated at room temperature for 48 h. The aggregation was monitored using ThT fluorescence. The sample was then centrifuged at 19000g for 12 min to spin out any fibril that may have formed. The supernatant was then fractionated by Superdex-75 size exclusion column to isolate protofibrils from unreacted monomers and smaller oligomers. The concentration of protofibrils was measured by UV-Vis with a molar extinction coefficient of 1450 *cm*^{–1}*M*^{–1} corresponding to A*β* 42.

### Dynamic light scattering (DLS)

The DLS experiments were performed using Zetasizer nanoseries instrument ( Malvern Instruments, Inc.). The samples were equilibrated for 120s and a total of 90 measurements of 10s each were made with a time interval of 1min after each measurement. The data was exported as diffusion coefficient using Zetasizer software version 6.20 and processed using origin 6.0 software.

### Reaction flux and differential equations for the reaction pathway

*F*

_{1600}whose rate of change is expressed in terms of its disappearance:

*F*

_{1600}enlisted in 5-species reaction model. Similarly, the differential equations of each of the other species involved can be written based on the fluxes of the reactions where the species is produced (positive terms) or consumed (negative terms):

The initial concentration of *F*_{1600}s is equal to the amount of protofibrils added initially and is the main driving force for the downstream reactions. The concentrations of the other species are assumed to be zero at the start (i.e., [*F*_{
i
}_{×1600}] = 0, *i* = 1, …, 5). This set of differential equations is properly defined and can be solved using Matlab’s ode15s.

## Authors contributions

The study was conceptualized by PG,VR and AV. PG designed the computational model for the reaction pathway and VG designed the ThT fluorescence and DLS experiments. AV,SA and BJC performed the mathematical modeling and the related calculations. All authors read and approved the final manuscript.

## Declarations

### Acknowledgments

The authors thank the following agencies for their financial support: NSF EAGER(1049962) and Mississippi INBRE (P20RR016476) (for PG and VR); American Heart Association (10GRNT4190124) (for VR). The authors thank Amit Kumar for his help with biophysical experiments and Dr. Cannon for letting us use the DLS instrument.

This article has been published as part of *BMC Bioinformatics* Volume 12 Supplement 10, 2011: Proceedings of the Eighth Annual MCBIOS Conference. Computational Biology and Bioinformatics for a New Decade. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/12?issue=S10.

## Authors’ Affiliations

## References

- Tycko R:
**Molecular structure of amyloid fibrils: insights from solid-state NMR.***Q Rev Biophys*2006, 1–55.Google Scholar - Harper JD, Lansbury PT Jr.:
**Models of amyloid seeding in Alzheimer’s disease and scrapie: mechanistic truths and physiological consequences of the time-dependent solubility of amyloid proteins.***Annu Rev Biochem*1997,**66:**385–407. 10.1146/annurev.biochem.66.1.385View ArticlePubMedGoogle Scholar - Walsh DM, Lomakin A, Benedek GB, Condron MM, Teplow DB:
**Amyloid beta-protein fibrillogenesis. Detection of a protofibrillar intermediate.***J. Biol. Chem*1997,**272:**22364–22372. 10.1074/jbc.272.35.22364View ArticlePubMedGoogle Scholar - Ward RV, Jennings KH, Jepras R, Neville W, Owen DE, Hawkins J, Christie G, Davis JB, George A, Karran EH, Howlett DR:
**Fractionation and characterization of oligomeric, protofibrillar and fibrillar forms of beta-amyloid peptide.***Biochem. J*2000,**348**(1):137–144. 10.1042/0264-6021:3480137PubMed CentralView ArticlePubMedGoogle Scholar - Murphy RM, Pallitto MM:
**Probing the kinetics of beta-amyloid self-association.***J Struct Biol*2000,**130:**109–122. 10.1006/jsbi.2000.4253View ArticlePubMedGoogle Scholar - Harper JD, Wong SS, Lieber CM, Lansbury PT Jr.:
**Assembly of A beta amyloid protofibrils: an in vitro model for a possible early event in Alzheimer’s disease.***Biochemistry*1999,**38:**8972–8980. 10.1021/bi9904149View ArticlePubMedGoogle Scholar - Blackley HK, Sanders GH, Davies MC, Roberts CJ, Tendler SJ, Wilkinson MJ:
**In-situ atomic force microscopy study of beta-amyloid fibrillization.***J. Mol. Biol*2000,**298:**833–840. 10.1006/jmbi.2000.3711View ArticlePubMedGoogle Scholar - Harper JD, Wong SS, Lieber CM, Lansbury PT:
**Observation of metastable Abeta amyloid protofibrils by atomic force microscopy.***Chem. Biol*1997,**4:**119–125. 10.1016/S1074-5521(97)90255-6View ArticlePubMedGoogle Scholar - Nichols MR, Moss MA, Reed DK, Lin WL, Mukhopadhyay R, Hoh JH, Rosenberry TL:
**Growth of beta-amyloid(1–40) protofibrils by monomer elongation and lateral association. Characterization of distinct products by light scattering and atomic force microscopy.***Biochemistry*2002,**41:**6115–6127. 10.1021/bi015985rView ArticlePubMedGoogle Scholar - Shen CL, Murphy RM:
**Solvent effects on self-assembly of beta-amyloid peptide.***Biophys. J*1995,**69:**640–651. 10.1016/S0006-3495(95)79940-4PubMed CentralView ArticlePubMedGoogle Scholar - Bitan G, Kirkitadze MD, Lomakin A, Vollers SS, Benedek GB, Teplow DB:
**Amyloid beta -protein (Abeta) assembly: Abeta 40 and Abeta 42 oligomerize through distinct pathways.**In*Proc. Natl. Acad. Sci*.*Volume 100*. U.S.A; 2003:330–335. 10.1073/pnas.222681699Google Scholar - Pallitto MM, Murphy RM:
**A mathematical model of the kinetics of beta-amyloid fibril growth from the denatured state.***Biophys. J*2001,**81:**1805–1822. 10.1016/S0006-3495(01)75831-6PubMed CentralView ArticlePubMedGoogle Scholar - Hepler RW, Grimm KM, Nahas DD, Breese R, Dodson EC, Acton P, Keller PM, Yeager M, Wang H, Shughrue P, Kinney G, Joyce JG:
**Solution state characterization of amyloid beta-derived diffusible ligands.***Biochemistry*2006,**45:**15157–15167. 10.1021/bi061850fView ArticlePubMedGoogle Scholar - Smith AM, Jahn TR, Ashcroft AE, Radford SE:
**Direct observation of oligomeric species formed in the early stages of amyloid fibril formation using electrospray ionisation mass spectrometry.***J. Mol. Biol*2006,**364:**9–19. 10.1016/j.jmb.2006.08.081View ArticlePubMedGoogle Scholar - Schuck P, Perugini MA, Gonzales NR, Howlett GJ, Schubert D:
**Size-distribution analysis of proteins by analytical ultracentrifugation: strategies and application to model systems.***Biophys. J*2002,**82:**1096–1111. 10.1016/S0006-3495(02)75469-6PubMed CentralView ArticlePubMedGoogle Scholar - Ghosh P, Kumar A, Datta B, Rangachari V:
**Dynamics of protofibril elongation and association involved in Abeta42 peptide aggregation in Alzheimer’s disease.***BMC Bioinformatics*2010,**11**(Suppl 6):S24. 10.1186/1471-2105-11-S6-S24PubMed CentralView ArticlePubMedGoogle Scholar - Cantor CR, Shimmel PR:
*Biophysical chemistry: Part II: Techniques for the study of biological structure and function.*W.H. Freeman; 1980.Google Scholar - Zwanzig R, Harrison AK:
**Modifications of the Stokes-Einstein formula.***Journal of Chemical Physics*1985,**83**(11):5861–5863. 10.1063/1.449616View ArticleGoogle Scholar - Gregorova E, Pabst W, Bouchet JB:
**Influence of particle shape on the viscosity of Kaolin suspensions.***Acta Geodyn. Geomater*2009,**6**(1):101–109.Google Scholar - Shenhar R, Wang H, Hoffman RE, Frish L, Avram L, Willner I, Rajca A, Rabinovitz M:
**Self-Assembled, helically stacked anionic aggregates of 2,5,8,11 tetra tert butylcycloocta [1,2,3,4-def;5,6,7,8-def]bisbiphenylene, stabilized by electrostatic interactions.***J. Am. Chem. Soc*2002,**124:**4685–4692. 10.1021/ja012140uView ArticlePubMedGoogle Scholar - Vasanthi R, Ravichandran S, Bagchi B:
**Needlelike motion of prolate ellipsoids in the sea of spheres.***J. Chem. Phys*2001,**114:**7989. 10.1063/1.1363674View ArticleGoogle Scholar - Stein A, Hoffmann DA, Marcus AH, Leezenberg PB, Frank CW, Fayer MD:
**Dynamics in poly(dimethylsi1oxane) melts: fluorescence depolarization measurements of probe chromophore orientational relaxation.***J. Phys. Chem*1992,**96:**5255–5263. 10.1021/j100192a017View ArticleGoogle Scholar - Kim MS:
**Study of molecular reorientation in liquid with Raman spectroscopy (III). Temperature dependence of molecular rotation of C6Fc in neat liquid.***J. of Korean Chem. Soc*1984,**28**(1):34–40.Google Scholar - Massoudi M:
**A note on the meaning of mixture viscosity using the classical continuum theories of mixtures.***International Journal of Engineering Science*2008,**46:**677–689. 10.1016/j.ijengsci.2008.01.008View ArticleGoogle Scholar - Sampaio R, Williams WO:
**Viscosities of liquid-mixtures, Zeitschrift Angew.**In*Math. Phys*.*Volume 28*. ZAMP; 1977:607–614.Google Scholar - Chwang AT, Wu TY:
**Hydromechanics of low-reynolds-number flow. Part 2. Singularity method for Stokes flows.***Journal Fluid Mechanics*1972,**45:**105–145.Google Scholar - Happel V, Brenner H:
*Low Reynolds number hydrodynamics.*Prentice Hall; 1965.Google Scholar - Chang C, Powell RL:
**Dynamic simulation of bimodal suspensions of hydrodynamically interacting spherical particles.***Journal of Fluid Mechanics*1993,**253:**1–25. 10.1017/S0022112093001697View ArticleGoogle Scholar - Farris RI:
**Prediction of the viscosity of multimodal suspensions from unimodal viscosity data.***Journal of Rheology*1968,**12:**281. 10.1122/1.549109View ArticleGoogle Scholar - Fidleris V, Whitmore RL:
**The physical interaction of spherical particles in suspensions.***Rheologica Acta*1961,**I:**573.View ArticleGoogle Scholar - Morr R, Gottlieb M, Graham A, Mondy L:
**Viscosity of concentrated suspensions of sphere/rod mixtures.***Chemical Engineering Communications*1996,**148**(1):421–430.Google Scholar - Milliken WJ, Mondy LA, Gottlieb M, Graham AL, Powell RL:
**The effect of the diameter of falling balls on the apparent viscosity of suspensions of spheres and rods.***Phys. Chem. Hydro*1989,**11:**341.Google Scholar - Thomas DG:
**Transport characteristics of suspension: VIII. A note on the viscosity of Newtonian suspensions of uniform spherical particles.***Journal of Colloid Science*1965,**20:**267. 10.1016/0095-8522(65)90016-4View ArticleGoogle Scholar - Quemada D:
**Rheology of concentrated disperse systems. I. Minimum energy dissipation principle and viscosity-concentration relationship.***Rheologica Acta*1977,**16**(1):82–94. 10.1007/BF01516932View ArticleGoogle Scholar - Quemada D:
**Rheological modelling of complex fuids. I. The concept of effective volume fraction revisited.***European Physical Journal of Applied Physics*1998,**1**(1):119–127. 10.1051/epjap:1998125View ArticleGoogle Scholar - Lee JD, So JH, Yang SM:
**Rheological behavior and stability of concentrated silica suspensions.***Journal of Rheology*1999,**43**(5):1117–1140. 10.1122/1.551018View ArticleGoogle Scholar - Wolthers W, Van den Ende D, Duits MHG, Mellema J:
**The viscosity and sedimentation of aggregating colloidal dispersions in a Couette flow.***Journal of Rheology*1996,**40**(1):55–67. 10.1122/1.550784View ArticleGoogle Scholar - Perrin F:
**Brownian movement of an ellipsoid.***J. Phys. Radium*1936,**7:**1–11. 10.1051/jphysrad:01936007010100View ArticleGoogle Scholar - Blackley HG, Sanders GH, Davies MC, Roberts CJ, Tendler Sj, Wilkinson MJ:
**In-situ atomic force microscopy study of beta-amyloid fibrillization.***J Mol Biol*2000,**298:**833–840. 10.1006/jmbi.2000.3711View ArticlePubMedGoogle Scholar - Nichols MR, Moss MA, Reed Dk, Lin WL, Mukhopadhyay R, Hoh JH, Rosenberry TL:
**Growth of beta-amyloid(1–40) protofibrils by monomer elongation and lateral association. Characterization of distinct products by light scattering and atomic force microscopy.***Biochemistry*2002,**41:**6115–6127. 10.1021/bi015985rView ArticlePubMedGoogle Scholar - Goldbury C, Frey P, Olivieri V, Aebi U, Muller S:
**Multiple assembly pathways underlie amyloid-beta fibrill polymorphisms.***J. Mol. Biol*2005,**352:**292–298.Google Scholar - Barthelmesa SEG., Pratsinisa H, Buggisch B:
**Particle size distributions and viscosity of suspensions undergoing shear-induced coagulation and fragmentation.***Chemical Engineering Science*2003,**58:**2893–2902. 10.1016/S0009-2509(03)00133-7View ArticleGoogle Scholar - Naoya H, Ikuro S:
**Experiments on granular rheology: Effects of particle size and fluid viscosity.***Journal of Geophysical Research*2009,**114:**B04413.Google Scholar - Phan-Thien N, Pham DC:
**Differential multiphase models for polydispersed pheroidal inclusions: thermal conductivity and effective viscosity.***International Journal of Engineering Science*2008,**38:**73–88.View ArticleGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.