Aggregate Size Dependence of Amyloid Adsorption onto Charged

Dec 28, 2017 - ⊥Central European Institute of Technology and #Faculty of Science, Masaryk University, 625 00 Brno, Czech Republic ... Abstract Image...
0 downloads 11 Views 14MB Size
Article Langmuir Aggregate is published Sizeby the American Dependence Chemical Society. of 1155 Subscriber access Sixteenth

provided Street by Grand

Amyloid Adsorption Langmuir is published onto by the American Charged Chemical Society. Interfaces 1155

Subscriber access Sixteenth provided Street by Grand

Giulio Tesei, Erik Langmuir Hellstrand, is published Kalyani by the Sanagavarapu, American SaraChemical S. Society. Linse, 1155 Emma Subscriber access Sixteenth provided Street by Grand

Sparr, Robert Vacha, and Langmuir MikaelisLund published

Langmuir, Just by the Accepted American Chemical Manuscript • Society. DOI: 10.1021/ 1155 acs.langmuir.7b03155 Subscriber access Sixteenth provided Street by Grand

• Publication Date (Web): 28 DecLangmuir 2017 is

Downloaded published from http:// by the American pubs.acs.org Chemical on December Society. 30, 2017 1155

Subscriber access Sixteenth provided Street by Grand

Just Accepted

Langmuir “Just Accepted” manus is online published prior to technica the Societyby provides “Jus American dissemination of scient Chemical appearSociety. in full in PDF fo 1155 fully peer reviewed, but Subscriber access Sixteenth

provided Street by Grand

readers and citable by to authors. Therefore, Langmuir After a m in the journal. is Accepted” Web site an published by the changes to the manus American and ethical guidelines Chemical Society. or consequences arisin 1155

Subscriber access Sixteenth provided Street by Grand

64

16

Page 1 of 36 Langmuir

8

32

2

128

1 2 segr ai ned l i nesegment 3coar ACS Paragon mo del Plus Environment model 4 5

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

cross section

A

side view

1 42

N = 64

A

1 40

N = 96

A

Page 2 of 36

17 42

N = 64 ACS Paragon Plus Environment

3 0 3 6

C

E

3 0

3 6

0 5 10 15 0 20 40 60 0 5 10 15

B

D

F

0 5 10 15 0 5 10 15 0 5 10 15 t (min) t (min) t (min) ACS Paragon Plus Environment

6

A

D / 10

3 0 3 6

POPC:DOTAP 3:1 cs = 0.05 M cs = 0.2 M

3 0

6

Monomers F/n (Hz)

Aggregates F/n (Hz)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36

POPC cs = 0.05 M

Langmuir

D / 10

Page 3 of 36

3 6

Free Energy, w(r)

1

Henry's Law Constant, KHA (nm)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

A

Page 4 of 36

0 1 2 8 16

2 4 107 105

0

4 8 12 16 Aggregate-Surface Separation, r (nm)

Attraction

101 0

Repulsion

100

101 102 Aggregate Size, N

ACS Paragon Plus Environment

20

B c (M) s

103

101

32 64 128

0.15 0.2 0.25 0.3 0.35 0.4 0.45 n. s. 103

107 105

Langmuir

A

Attraction

cs (M)

103

0.2 0.25 0.3 0.35 0.4 0.45 n. s.

101 0 101

Repulsion

100

Henry's Law Constant, KHA (nm)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Henry's Law Constant, KHA (nm)

Page 5 of 36

107 105

101

102

Attraction

103

cs (M)

101 0 101

B

103

Repulsion

100

101 102 Aggregate Size, N

ACS Paragon Plus Environment

0.1 0.15 0.2 0.25 0.3 n. s. 103

Langmuir

1

A

B

C

0

2 1 0

Total Charge

2 Total Charge

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 36

1

1

2

2

3

3

0 2 4 6 0 2 4 6 0 2 4 6 Rc (nm) Rc (nm) Rc (nm)

ACS Paragon Plus Environment

4 3 2 1

0 0.35

0

2)

0.25

Surface Excess, × 103 (nm 2)

Repulsion

0 0.15 0.25 0.35 Salt Concentration, cs (M) ACS Paragon Plus Environment

80

Adsorbed Amount, m (ng cm

Aggregate Size, N

Attraction

80

B

Attraction

600 500 400 300 200 100 0 0.05

03 - 0.0

Aggregate Size, N

A

600 500 400 300 0.5 200 100 1.51.0 0 0.05 0.15

03 -0.0 10

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56

Langmuir

Repulsion

Page 7 of 36

60 40 20 0

A

4

B

cs (M) 0.05 0.10 0.15 0.20

3 2

cs (M) 0.05 0.10 0.15 0.20

Page 8 of 36

120 100 80 60 40

1

20

0

0

0

40

80 120 0 L (nm)

40

ACS Paragon Plus Environment

80 120 L (nm)

m (ng cm 2)

× 103 (nm 2)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32

5

Langmuir

Page 9 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Aggregate Size-Dependence of Amyloid Adsorption onto Charged Interfaces Giulio Tesei,∗,† Erik Hellstrand,‡ Kalyani Sanagavarapu,¶ Sara Linse,¶ Emma Sparr,§ Robert Vacha,∗,k,⊥ and Mikael Lund∗,† Theoretical Chemistry, Lund University, 221 00 Lund, Sweden, Biophysical Chemistry, Lund University, 221 00 Lund, Sweden, Biochemistry & Structural Biology, Lund University, 221 00 Lund, Sweden, Physical Chemistry, Lund University, 221 00 Lund, Sweden, Central European Institute of Technology, Masaryk University, 625 00 Brno, Czech Republic, and Faculty of Science, Masaryk University, 625 00 Brno, Czech Republic E-mail: [email protected]; [email protected]; [email protected]

Abstract Amyloid aggregates are associated with a range of human neurodegenerative disorders, and it has previously been shown that neurotoxicity is dependent on aggregate size. Combining molecular simulation with analytical theory, a predictive model is proposed for the adsorption of amyloid aggregates onto oppositely charged surfaces, where the interaction is governed by an interplay between electrostatic attraction and ∗

To whom correspondence should be addressed Theoretical Chemistry, Lund University, 221 00 Lund, Sweden ‡ Biophysical Chemistry, Lund University, 221 00 Lund, Sweden ¶ Biochemistry & Structural Biology, Lund University, 221 00 Lund, Sweden § Physical Chemistry, Lund University, 221 00 Lund, Sweden k Central European Institute of Technology, Masaryk University, 625 00 Brno, Czech Republic ⊥ Faculty of Science, Masaryk University, 625 00 Brno, Czech Republic †

1

ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

entropic repulsion. The theoretical predictions are experimentally validated against quartz crystal microbalance-dissipation experiments of amyloid beta peptides and fragmented fibrils in the presence of a supported lipid bilayer. Assuming amyloids as rigid, elongated particles, we observe non-monotonic trends for the extent of adsorption with respect to aggregate size, and preferential adsorption of smaller aggregates over larger ones. Our findings describe a general phenomenon with implications for stiff-polyions and rod-like particles that are electrostatically attracted to a surface.

Introduction A large number of peptides and proteins self-assemble into elongated and highly ordered structures, rich in β-sheets, which are generally termed amyloid fibrils. 1–3 The formation process is controlled by amino acid sequence; solution conditions such as pH, and concentration of salts or co-solutes; as well as by the presence of surfaces in contact with the solution. In biology, protein aggregation occurs in environments containing proteins and protein complexes, glyco-proteins, nucleic acids, ribosome particles, and lipid membranes. The protein self-assembly process and the final composition of the amyloid aggregates are affected by this environment. 4–7 Likewise, protein aggregation influences properties of existing self-assembled entities, for example the structure and integrity of biological membranes. 4,8,9 Amyloid formation is associated with many human diseases. 2 To understand cellular toxicity, it is important to resolve the role played by the interactions between proteins, in different aggregation states, and other components of the complex environment. Several amyloid proteins, e.g. Aβ, α-synuclein, and IAPP are surface active and adsorb to solid surfaces, air-liquid interfaces, and to lipid bilayers. 10–14 Adsorption of protein aggregates to a membrane interface is the first step towards cell permeabilization and in this work we investigate the interaction between elongated aggregates of charged peptides and an oppositely charged surface. Aggregate size is one possible determinant of the observed cytotoxicity. Several studies have 2

ACS Paragon Plus Environment

Page 10 of 36

Page 11 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

shown that oligomers and short, fragmented fibrillar aggregates are more potent in permeabilizing cell membranes, and in reducing cell viability, compared to longer fibrils. 15–18 Fragmentation of β2 -microglobulin, α-synuclein, and lysozyme fibrils into shorter aggregates has been shown to cause increased cellular damage, 19 and the enhanced cytotoxicity has been related to the enrichment in fibril extremities upon fragmentation. 19,20 Indeed, interactions between the extremities of short β2 -microglobulin fibrils and negatively charged liposomes has been reported to induce membrane distortions via lipid extraction from the bilayer. 20 Analogously, short α-synuclein fibrils have been shown to associate with negatively charged lipid bilayers and form protein-lipid co-aggregates. 4 Moreover, it has been reported that long (100–400 nm) and short (10–100 nm) α-synuclein aggregates have comparable affinity to the plasma membrane of mammalian cells and they bind with predominantly parallel orientation to liposomes composed of brain lipids. 21 Lateral binding has also been observed for huntingtin exon 1 fibrils of length 40–120 nm, while shorter fragments displayed low affinity to both liposomes and cells. 21 We use the quartz crystal microbalance-dissipation (QCM-D) 22 technique to measure the interaction of Aβ1−40 in either monomeric or fibrillar form with fluid lipid bilayers. Our model system consists of a combination of Aβ1−40 and a positively charged POPC:DOTAP 3:1 lipid bilayer. Experimental results are interpreted using molecular simulations and an analytical theory, based on the representation of aggregates as charged line segments. Molecular simulations show how the adsorption is influenced by aggregate size, and decay length of the electrostatic attraction to the surface, which is controlled by the bulk ionic strength. The analytical line segment theory captures the dependence of surface–aggregate interaction on bulk ionic strength and aggregate size. Using generalized van the Waals theory, we calculate the surface excess of line segments, a quantity which can be directly related to the QCM-D data.

3

ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 36

Materials and Methods 1-palmitoyl-2-oleoyl−sn−glycero-3-phosphocholine (POPC) and 1,2-dioleoyl-3-trimethylammoniumpropane (DOTAP) were purchased lyophilized from Avanti Polar Lipids Inc. (Alabaster, USA). All chemicals were of analytical grade.

Amyloyd-β peptides and fibrils. The amiloyd-β peptide Aβ1−40 was expressed in Escherichia Coli and purified according to Walsh. 23 Aβ1−40 monomers were isolated by size exclusion chromatography in 20 mM sodium phosphate buffer, 0.02% NaN3 , 0.2 mM EDTA, pH 7.4. For samples with salt concentration 150 mM, NaCl was added to the monomer solution on ice from a 30× concentrated stock. Fibrils were prepared by keeping the monomer solution at 37◦ C with moderate shaking (IKA-VIBAX-VXR motor operating at 200 rpm) in low binding tubes (Axygen). Fibrils were sonicated for 30 minutes with a tip sonicator using cycles of 1-s pulses and 1-s waiting time. The length of sonicated fibrils was characterized using atomic force microscopy (AFM) (Supporting Information, Figure S4). The peptide concentration of the samples used for QCM-D experiments was 4 µM and for AFM 10 µM. The contribution of the sodium phosphate buffer to the ionic strength of the solutions is 50 mM.

Bilayer formation. For the positively charged lipid bilayer, a POPC:DOTAP 3:1 mixture was prepared in chloroform:methanol 2:1 v:v and deposited as a thin film on glass under flow of N2 gas and dried under vacuum overnight. The lipid film was dispersed in 200 mM NaCl aqueous solution to obtain a concentration of 0.5 mg/ml. For the neutral lipid bilayer, POPC was dispersed in water to obtain a concetration of 0.5 mg/ml. Lipid dispersions were probe-sonicated in an ice bath for 10 minutes alternating 5-second pulses with 5 seconds of cooling. The clear vesicle dispersions were centrifuged for 20 minutes at 2000 rcf to remove debris from the probe tip.

4

ACS Paragon Plus Environment

Page 13 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

The supernatant which contains small unilamellar vesicles was used for the deposition of the lipid bilayer.

QCM-D Experiments. The Q-sense E4 from Qsense (Göteborg, Sweden) is the QCM-D instrument used for the experiments. Measurements were performed in parallel in four cells thermostated at 25◦ C with quartz crystals covered by a thin gold surface coated with 50 nm SiO2 (QSX 303, Q-sense). The crystals were stored in 2% w/v SDS solution for at least one hour, rinsed in MilliQ water followed by ethanol and dried with N2 gas. Finally the crystals were treated in the low pressure (0.02 mbar) chamber of a plasma cleaner (Harrik Scientific Corp., Pleasentville, NY, USA, model PDC-3XG) for 10 minutes. The crystals were mounted in the cells and equilibrated with water until steady frequency and dissipation responses were observed. In the experiments with the positively charged lipid bilayer, an injection of a 200 mM NaCl aqueous solution, with flow-rate of 300 µl/min, preceded the injection of the vesicle dispersion, with flow rate of 100 µl/min. The high ionic strength of the dispersion facilitates vesicle formation and the deposition of the lipid bilayer by screening the electrostatic repulsions between the positively charged DOTAP head-groups. The supported lipid bilayer was equilibrated in the buffer solution and finally the freshly isolated Aβ1−40 monomers or the sonicated Aβ1−40 fibrils were injected. Both buffer and protein solution flowed through the cells at a rate of 50 µl/min. 24–26

MC Simulations. We performed Metropolis Monte Carlo (MC) simulations using the Faunus framework. 27 A single aggregate (Figure 1) is placed in a rectangular box of volume 200×200×100 nm3 , and a wall of uniform charge density is located at z=0 nm. The potential energy contribution of the ith bead of radius Ri at a distance ri from the surface is given by the linearized

5

ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 36

Gouy-Chapman potential 28

ui (r) =

   − zi eφ exp (−ri /λD ), ri > Ri

(1)

  ∞, ri ≤ Ri where zi is the charge number of the amino acid, e is the elementary charge, and λD = √ 3.04/ cs Å is the Debye screening length, at T = 298 K, for the surrounding monovalent (1:1) salt of molar concentration, cs . The surface potential is βφe = 2 asinh (2πσλB λD /e), where λB = 7.1 Å is the Bjerrum length of water, β = 1/kB T is the inverse thermal energy, while σ = 1/266 eÅ−2 is the surface charge density of a 3:1 POPC:DOTAP lipid bilayer. 29 The internal degrees of freedom of the aggregates are frozen, and the configurational space is sampled by attempting roto-translational moves of the rigid body. Sampling is enhanced using the Wang-Landau method, 30 where the reaction coordinate is the separation of the center-of-mass of the aggregate from the surface, r. To study the effect of counterion condensation, a single coarse-grained aggregate is placed in a cylindrical box with the fibril principal axis aligned to the cylinder axis. We performed MC simulations of the static aggregate surrounded by explicit mobile monovalent (cs =0.25 mM, 0.5 mM, 5 mM) or divalent (cs =2 mM) counterions of radius 1.9 Å. 31 We calculated the average total charge within coaxial cylinders of increasing radii enclosing the aggregate and the ions. The effective charge of the aggregate is determined as the total charge at the radial distance corresponding to the peak in the counterion density profile.

Results and Discussion QCM-D Experiments. QCM-D provides measurements of the temporal variation in frequency (∆F ) and energy dissipation (∆D) of a quartz crystal resonator when the lipid or protein and associated solvent molecules adsorb on the sensor surface. Figure 2 shows the measured |∆F |/n and ∆D 6

ACS Paragon Plus Environment

Page 15 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

cross section

A

side view

1 42

N = 64

A

1 40

N = 96

A

17 42

N = 64

Figure 1: Cross sections and side views of coarse-grained aggregates of N Aβ1−42 , Aβ17−42 and Aβ1−40 peptides. Each bead corresponds to an amino acid that can be neutral (white), cationic (blue), or anionic (red). Fibrils of the Aβ17−42 fragment present a smooth rod-like structure. The additional 16 N-terminal residues of the Aβ1−42 peptide are mostly charged or titratable, and they decorate the core of Aβ1−42 fibrils with flexible chains which extend in the surrounding solution, 32,33 resembling a polyelectrolyte brush. Aβ17−42 and Aβ1−42 structures have two-fold symmetry, whereas the fibril made of Aβ1−40 peptides has three peptides per cross section and three-fold symmetry. for solutions of Aβ1−40 monomers and fragmented fibrils in the presence of a supported POPC or POPC:DOTAP 3:1 lipid bilayer. Signals are recorded for various overtone numbers, n, of the fundamental frequency of the quartz crystal. The zero value in Figure 2 corresponds to the stationary response after deposition of the lipid bilayer and equilibration with the buffer solution. Figure 2 A and B show that negatively charged Aβ1−40 monomers and fragmented fibrils in solutions of pH 7.4 and 50 mM ionic strength, cs , do not adsorb to neutral POPC lipid bilayers. Therefore, despite being surface active, the association of Aβ to fluid lipid bilayers is not controlled by hydrophobic interactions. 34 Panels d, e, and f of Figure 2 show that in the presence of the positively charged POPC:DOTAP 3:1 lipid bilayer, no significant changes in |∆F |/n and ∆D are detected upon addition of Aβ40 monomers, neither at high nor at low ionic strength. Similarly, when sonicated fibrils are added at high salt condition, no adsorption on the deposited bilayer is detected. However, adsorption clearly occurs after the injection of sonicated Aβ1−40 fibrils in the buffer with no added NaCl, as directly inferred from the decrease in ∆F/n by around 5-10 Hz and from the simultaneous increase in ∆D 7

ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

in Figure 2C. The dissipation shift and the n-dependent response observed upon injection of Aβ1−40 fibrils at low ionic strength conditions, indicate that the adsorbed fibril layer on the positively charged lipid bilayer is viscoelastic and acoustically coupled to the aqueous solution. Hence, the experimental data were analyzed with a one-layer extended viscoelastic Voigt-based model 35,36 (Supporting Information, QCM-D Data Analysis). In the continuum model, the adsorbed fibrils are represented by a viscoelastic layer, while the supported lipid bilayer is considered to be unaffected by fibril adsorption. The fibril layer is in contact with a semi-infinite solution modelled as a Newtonian fluid. The extended model accounts for a linear frequency dependence of the viscosity and shear modulus of the fibril layer, 37 expressed in terms of frequency factors. The parameters evaluated in the fitting procedure are the frequency factors and, more importantly, the time evolutions of the shear modulus, viscosity, and mass per area of the fibril layer (Supporting Information, Table S1 and Figure S3). The fitted quantities significantly increase upon injection of Aβ1−40 fibrils, and the equilibrium wet mass of adsorbed hydrated fibrils is estimated to be of 660 ± 30 ng cm−2 .

Molecular Simulations. Computer simulations are used to calculate the interaction free energy between Aβ elongated aggregates and a planar surface as a function of aggregate size and salt concentration. The strategy is to first use all-atom molecular dynamics (MD) simulations to relax the structure of fibrillar assemblies (Supporting Information), which are then coarse-grained to the amino acid level and used in Metropolis MC simulations of surface–aggregate interactions. Previous studies provided evidence of a difference in morphology between amyloid oligomers, protofibrils and fibrils. 38–40 To take polymorphism of amyloid aggregates into account, our in silico investigation extends over three fibril architectures (Figure 1): two structures with two-fold symmetry for Aβ1−42 peptides 33 and Aβ17−42 fragments, 41,42 as well as a structure with three-fold symmetry for Aβ1−40 peptides, which is predominant in Alzheimer’s disease 8

ACS Paragon Plus Environment

Page 16 of 36

E

3 0

3 6

0 5 10 15 0 20 40 60 0 5 10 15

B

D

F

0 5 10 15 0 5 10 15 0 5 10 15 t (min) t (min) t (min)

6

C

D / 10

3 0 3 6

A

POPC:DOTAP 3:1 cs = 0.05 M cs = 0.2 M

3 0

6

3 0 3 6

POPC cs = 0.05 M

D / 10

Aggregates F/n (Hz)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Monomers F/n (Hz)

Page 17 of 36

3 6

Figure 2: Frequency (blue lines) and dissipation changes (red lines) for overtone-numbers 5, 7, 9 after injection of (A, C, E) fibrils solutions and (B, D, F) monomer solutions of peptide concentration 4 µM. The zero values of ∆F/n and ∆D correspond to the supported lipid bilayer equilibrated with the buffer solution. No adsorption is observed for monomers in both low (cs =0.05 M) and high salt conditions (cs =0.2 M) on POPC as well as on POPC:DOTAP 3:1 bilayers. For fibrils, significant adsorption is detected only at cs =0.05 M on the POPC:DOTAP 3:1 bilayer (C). The corresponding ∆F/n and ∆D are fitted to the one-layer extended Voigt-based model (C, black lines). The yellow symbols facilitate the comparison with the predictions of the line segment model shown in Figure 6B. human brain tissues. 43,44 Assemblies of around 200 peptides are generated by stacking fibril segments obtained from the protein data bank (PDB entries: 5KK3, 2M4J and 2BEG, after Zheng et al. 42 ), and simulated annealing 45 is used to minimize the energy of the large structures. At pH 7.4, Aβ1−42 and Aβ1−40 have the same net charge number, zm , of -3, while the Aβ17−42 fragment has zm = −1. All-atom elongated aggregates are coarse-grained to rigid models where each amino acid is represented by a bead that can be neutral, cationic or anionic. Approximating aggregates as rigid bodies is justified by the stiffness of Aβ fibrils, which have persistence lengths of micrometers. 43,46,47 While the distinction between oligomer and fibril may be operational, structural, or size-based, 48 for our simulations we generate aggregates of various numbers of monomers, N , by cross-sectional slicing of the long fibril. This approach is analogous to the physical fibril fragmentation process that has previously been used to show the enhanced cytotoxic potential of small-sized amyloid fragments. 19 In MC simulations of the rigid aggregates in the presence of a charged, planar interface, we

9

ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 36

evaluate the surface–aggregate interaction free energy using *

+

βw(r) = − ln exp(−β

X

ui (ri ))

i

(2) Ω

where the angular brackets denote a Canonical ensemble average over the orientational degrees of freedom of the aggregate, Ω; ui is the interaction of the ith amino acid with the surface at a distance ri (Equation 1); and r is the aggregate mass center to surface separation. The relative protein concentration in the interfacial volume, with respect to the bulk, is quantified by the surface excess, Γ. 28 An aggregate surplus (Γ > 0) or depletion (Γ < 0) is determined by surface–protein and protein–protein interactions. For a weakly interacting system, Γ can be described by Henry’s law: Γ ≈ KHA ρ/N , where ρ is the bulk peptide number density, and KHA is Henry’s law constant, related to w(r) through the Mayer integral Z



dr [exp(−βw(r)) − 1] .

KHA =

(3)

0

Figure 3A shows free energy profiles obtained from Equation 2 for Aβ1−42 assemblies of various N at cs =0.4 M. For smaller aggregates, the interaction is short-ranged and repulsive. Conversely, free energy profiles for larger aggregates are characterized by a minimum for the fibril in proximity of the surface, and by a free energy barrier between surface and bulk solution. These profiles result from the interplay between surface–aggregate electrostatic attraction and the entropic repulsion due to the decrease in rotational degrees of freedom as the aggregate approaches the surface. Both attractive and repulsive forces are heightened with increasing fibril length. The free energy barrier increases with N , suggesting that, in real systems, larger aggregates might be kinetically trapped in solution. The subtle kinks occurring at r around 5, 9 and 18 nm for aggregates of 32, 64 and 128 peptides, respectively, reflect the possibility for the aggregate to be oriented perpendicularly to the surface. This orientation is favored by end-point electrostatic attraction (Supporting Information, Figure 10

ACS Paragon Plus Environment

Page 19 of 36

S7), as well as by the absence of steric hindrance from the surface. At cs =0.4 M, aggregates of N = 64 preferentially bind to the surface through their extremities, while they bind laterally at cs =0.35 M. Aggregates of N = 128 bind laterally at all explored cs values. This illustrates that cs and N modulate the preferential orientation of the adsorbed aggregate. Figure 3B shows Henry’s law constants, KHA , calculated from MC simulations for various

A

Free Energy, w(r)

1 0 1 2 8 16

2 4 Henry's Law Constant, KHA (nm)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

107 105

0

4 8 12 16 Aggregate-Surface Separation, r (nm)

cs (M) 0.15 0.2 0.25 0.3 0.35 0.4 0.45 n. s.

101 0

Repulsion

100

20

B

Attraction

103

101

32 64 128

101 102 Aggregate Size, N

103

Figure 3: (A) Angularly averaged interaction free energy, w(r), as a function of surface– aggregate separation for Aβ1−42 assemblies of various size, N , and 0.4 M salt concentration, cs . (B) Henry’s law constants, KHA , calculated from MC simulations (circles) and from the line segment model (lines) for aggregates of Aβ1−42 of increasing N and at different cs . Data points labeled n.s. are from calculations with a neutral surface. N and cs . At conditions where the electrostatic interactions are negligible, i.e. at high cs and for neutral surfaces (n.s.), the surface–aggregate interaction is repulsive at all separations, and more so for longer fibrils. This is expected as the interaction is controlled solely by the orientational entropy loss, when fibrils approach the interface. At low cs , the trend reverses, and long fibrils adsorb to the interface with interaction free energies of several kB T . Here, 11

ACS Paragon Plus Environment

Langmuir

the electrostatic attraction dominates the interaction, greatly exceeding the entropic cost of aligning the rod-like fibrils parallel to the surface. For intermediate cs , the interplay between entropic loss and electrostatic attraction results in a more complex behavior: KHA varies non-monotonically with N (brown line in Figure 3B), indicating that fibrillar assemblies of certain lengths may be repelled while others attracted. The same conclusion is drawn by inspecting the free energy profiles in Figure 3A. We now compare KHA values for aggregates of the Aβ1−42 fragment with the corresponding

Henry's Law Constant, KHA (nm)

107 105

Attraction

A

Repulsion

cs (M) 0.2 0.25 0.3 0.35 0.4 0.45 n. s.

103 101 0 101

100 107 Henry's Law Constant, KHA (nm)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 36

105

101

102

Attraction

B

Repulsion

cs (M) 0.1 0.15 0.2 0.25 0.3 n. s.

103 101 0 101

100

101 102 Aggregate Size, N

103

103

Figure 4: Henry’s law constants, KHA , calculated from MC simulations for (A) Aβ1−40 and (B) Aβ17−42 aggregates of increasing N and at different cs . Data points labeled n.s. are from calculations with a neutral surface. Lines are calculated using the line segment model. data for rigid aggregates of Aβ17−42 and Aβ1−40 peptides. The three amyloids have architectures differing in cross-sectional area and symmetry, line charge density, and extent of exposure of charged residues to the surface. Nonetheless, Figures 3 and 4 show that aggregates of Aβ1−42 , Aβ1−40 , and Aβ17−42 display similarities in the dependence of surface 12

ACS Paragon Plus Environment

Page 21 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

interaction on cs and N . Notably, we observe fair agreement between the dependence of KHA on N and cs for Aβ1−42 and Aβ1−40 aggregates which have different symmetry but similar line charge densities (Figure 3B and 4A). Therefore, it seems reasonable to approximate the elongated rigid aggregates by line segments with a characteristic line charge density. We begin by studying the adsorption of a negatively charged aggregate on a positively charged surface, to conclude that the underlying interactions are independent of the sign of the net charge of the interacting entities, as long as they are oppositely charged. As a consequence, the line segment model developed below is equally applicable to the case of a positively charged fibril interacting with a cell membrane.

Line Segment Model. In the following, we describe how the conceptual line segment model for surface–aggregate interactions is constructed. Consider a freely rotating line segment of length L, with the center point located at a distance r away from a planar surface. The entropy change corresponding to the reduced number of available rotational states is related to the relative area of the spherical belt spanned by the ends of the line segment with respect to a sphere of diameter L,  2r Abelt  = kB ln , r < L/2  kB ln Asphere L ∆S(r) =   0, r ≥ L/2.

(4)

Using the expression for the electrostatic interaction between a line charge and the unperturbed double layer of a charged surface, 49 the total surface–line interaction free energy is approximated by " βw(r)rod ≈ 2z ln

1 + tanh(βφe/4)e

− λr

1 − tanh(βφe/4)e

D

− λr

D

# (5)

−∆S(r)/kB + (s/r)9 where s is the distance of closest approach between aggregate and surface, while z is the net charge number of the aggregate. Each monomer contributes to the aggregate length

13

ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 36

and charge number by lm = L/N and zm = z/N , respectively. The r−9 -term introduces a soft repulsion between rod and surface. 50 To allow for the difference in dimensionality between oligomers and longer fibrils, s is modeled by a smooth function varying between s0 and s1 , s = s0 + (s1 − s0 ) tanh (N/ν). s0 , s1 , and ν are determined from global leastsquare fits to the simulated KHA values for high cs (Supporting Information, Table S2) while lm values of 2.8 Å, 2.2 Å, and 2.6 Å are obtained from long Aβ1−42 , Aβ1−40 , and Aβ17−42 coarse-grained aggregates, respectively. Figure 3B and 4 show KHA as a function of N as obtained from Equations 3 and 5 for the three different amyloids. The line segment model closely reproduces the MC simulation results: adsorption is enhanced by increasing N , and decreases with increasing cs . Further, the model captures the oscillating trends at cs values where entropic repulsion and electrostatic attraction are of comparable magnitude. The agreement between molecular simulation results and the analytical model indicates that the main features of the adsorption behavior are independent of molecular-level structural details, such as the discrete charge distribution on the fibril surface. Hence, the line segment model is also valid for positively charged elongated aggregates adsorbing onto a negatively charged surface. The calculation of Γ is based on the description of the adsorbed particles as a two-dimensional fluid, using a generalized van der Waals approach. The excluded area per monomer is set to 3σ −1 , implying that the maximum coverage occurs when the surface charge is neutralized by the adsorbed particles. Adsorbed fibrils laterally repel each other as cylinders of line charge density z/L and diameter s 51

βw(r)surface

λB z 2 K0 =2 L

14



r λD



+ (s/r)12

ACS Paragon Plus Environment

(6)

Page 23 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

where K0 is the zero order modified Bessel function of the second kind. Γ is given by the implicit equation 52 ρ Γ = KHA N

    3N Γ 3N Γ 1− exp 2Γˆ a− σ σ − 3N Γ

(7)

where KHA was defined in Equation 3, while a ˆ is the mean-field energy constant for the R∞ adsorbed aggregates a ˆ = − 12 s dr2πrβw(r)surface . The charge density of rod-like polyelectrolytes in solution is compensated by counterion condensation, 53,54 and highly charged rods are likely to release only a small fraction of condensed ions upon adsorption onto a weakly charged surface of opposite sign. 55 Therefore, for the calculation of Γ, we consider an effec∗ , estimated from Monte Carlo simulations of a single tive charge number per monomer, zm

fibril surrounded by its counterions. Figure 5 shows the cumulative sum of aggregate and counterion charges as a function of the radial distance from the aggregate axis, Rc . Even with monovalent counterions, the effective charge of Aβ1−42 and Aβ1−40 aggregates is significantly reduced by counterion condensation. In contrast, for the less densely charged Aβ17−42 aggregate, counterion condensation occurs only with divalent ions. For a wide range of N and cs , Figure 6A displays Γ values predicted by the line segment model using Equation 7 and the parameters derived for Aβ1−40 aggregates (Supporting Information, Table S2). Solid contour lines connect conditions of N and cs yielding constant Γ values. The line of Γ = −0.003 nm−2 indicates that, for high cs , the surface excess of aggregates varies non-monotonically with N . Points on the left-hand side of the line of zeroΓ correspond to conditions where the aggregates are attracted to the surface. The contour lines of Γ =0.5, 1.0, and 1.5 nm−2 highlight that our model predicts larger Γ values for oligomers than for fibrils, at low-to-intermediate cs . This stems in part from the decrease, with increasing N , of the number of aggregates required to neutralize the surface. Figure 5B displays line segment model estimates of the adsorbed amount, ∆m = N ΓMW /NA , where MW is the molecular weight of the Aβ1−40 peptide. This conversion allows for direct

15

ACS Paragon Plus Environment

Langmuir

1

A

B

C

0

2 1 0

Total Charge

2 Total Charge

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 36

1

1

2

2

3

0 2 4 6 0 2 4 6 0 2 4 6 Rc (nm) Rc (nm) Rc (nm)

3

Figure 5: Total charge of aggregate plus counterions (red lines) within coaxial cylinders as a function of the cylinder radius, Rc . Grey and blue lines are the total charge of the aggregate and counterions, respectively. Counterions are monovalent at cs = 0.25 mM (solid lines) or divalent at cs = 2 mM (dashed lines). Aggregates consist of (A) N = 120 Aβ1−42 peptides, (B) N = 150 Aβ1−40 peptides, and (C) N = 130 Aβ17−42 peptides. The dashed vertical lines indicate the positions of the peaks in the counterion radial density profiles which are used to determine the effective charge of the aggregates. comparison with QCM-D experimental results. Contrary to what we observe for Γ, at low cs , ∆m is higher for long fibrils than for oligomers. However, Figure 7 shows that at physiological ionic strength both Γ and ∆m are the highest for a range of relatively short aggregates with maximum values for lengths of 20 nm and 70 nm, respectively.

For crystals with fundamental frequency of 5 MHz, the Sauerbrey equation (Supporting Information, Equation S1) gives an estimate of around 10 ng cm−2 for the smallest adsorbed wet mass yielding a ∆F/n signal significantly lower than zero. Sonicated fibril samples at low and high cs have length distributions between 40 and 130 nm (Supporting Information, Figure S4), corresponding to N between 180 and 590. Figure 6B shows that ∆m is large in a range of cs values which extends to cs =0.17 for N =250. Beyond the upper bound of this cs -interval, the screened surface–aggregate electrostatic force and the entropic repulsion counterbalance. As indicated by the yellow symbols in Figure 2 and Figure 6B, the line segment model reproduces the QCM-D results, predicting significant adsorption for aggregates at low ionic strengths, while the binding affinity is considerably lower for monomers, and for sonicated fibrils at cs ≥ 0.2 M. 16

ACS Paragon Plus Environment

Surface Excess, × 103 (nm 2)

Repulsion

3 2

0.35

0 0.15 0.25 0.35 Salt Concentration, cs (M)

0

Adsorbed Amount, m (ng cm 2)

0 0.25

4

1

Repulsion

Attraction Attraction

80

B

03

600 500 400 300 200 100 0 0.05

-0.0

Aggregate Size, N

A

600 500 400 300 0.5 200 100 1.51.0 0 0.05 0.15

03 -0.0 10

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

Aggregate Size, N

Page 25 of 36

80 60 40 20 0

Figure 6: (A) Surface excess of aggregates, Γ, and (B) adsorbed amount, ∆m, as a function of salt concentration, cs , and aggregate size, N . Values are calculated using the line segment ∗ model for Aβ1−40 aggregates (zm = −1.35, lm = 2.2 Å, s0 = 19.3 Å, s1 = 72.7 Å, and ν = 1556). The bulk peptide number density is ρ=2.4 10−6 nm−3 . The yellow symbols in (B) highlight the conditions of N and cs explored in the QCM-D experiments of Figure 2. Contour lines connect conditions of fibril length and ionic strength yielding same values of Γ or ∆m. At cs = 0.05 M, the predicted adsorbed dry peptide mass is ∆m ≈90 ng cm−2 , i.e. 13.6% of the wet mass determined from QCM-D experiments. Assuming this composition, the protein layer has density of 1.05 g cm−3 and a thickness of 6.3 ± 0.3 nm, implying that, on average, the fibril principal axis forms angles smaller than 10◦ with the surface. The QCM-D signals displayed in Figure 2 are characteristic of a highly hydrated protein layer, and the quasiparallel fibril orientation is expected for a strong surface–aggregate attraction.

17

ACS Paragon Plus Environment

Langmuir

4 3 2

A

B

cs (M) 0.05 0.10 0.15 0.20

cs (M) 0.05 0.10 0.15 0.20

1 0

120 100 80 60 40

m (ng cm 2)

5

× 103 (nm 2)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 36

20 0

40

80 120 0 L (nm)

40

80 120 L (nm)

0

Figure 7: (A) Surface excess of aggregates, Γ, and (B) adsorbed amount, ∆m, calculated as a function of aggregate length, L, from the line segment model for Aβ1−40 aggregates ∗ = −1.35, lm = 2.2 Å, s0 = 19.3 Å, s1 = 72.7 Å, and ν = 1556) and four values of cs . (zm The surface charge density is σ=1/266 eÅ−2 , while the bulk peptide number density is ρ=2.4 10−6 nm−3 .

Summary and Conclusions In summary, the proposed model offers a description of the adsorption of rigid rod-like amyloid aggregates to oppositely charged surfaces, and predicts that small aggregates can lead to higher surface excess and adsorbed amount than long fibrils at physiological conditions. This theoretical result offers a feasible explanation for why longer fibrils might be less cytotoxic than shorter ones, 16,17 albeit oligomer flexibility and membrane structure may also play significant roles. 56 The model predicts that at low ionic strength, cs , the number of adsorbed molecules at the surface is larger for short than for long aggregates, while the dependence of the extent of adsorption on N is less pronounced at higher cs . Moreover, large N and low cs favor lateral surface binding while end-point binding is favored by shorter length and higher cs . We show that small changes in solution ionic strength as well as fibril line charge density and length have a large impact on amyloid adsorption. This marked sensitivity on system conditions contributes to explain the seemingly contradictory experimental evidence regarding the length-dependence of the affinity and preferential orientation of positively charged

18

ACS Paragon Plus Environment

Page 27 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

amyloid aggregates to plasma membranes. 19–21 In particular, amyloid adsorption is likely to be affected by the low concentrations of divalent cations present in the extracellular space in vivo as well as in cell culture media. Compared to monovalent ions, divalent counterions reduce to a larger extent the effective line charge density of the aggregate 54 (Figure 5). However, they also mediate the interaction between adsorbed like-charged aggregates and may lead to a net lateral attraction at the interface. 57 Additionally, the presented adsorption mechanisms described in this work is of interest for a variety of colloids and molecules that can be likened to rigid charged elongated particles– DNA strands and polyelectrolytes, 58–61 rod-like particles, 62–64 as well as amyloid aggregates in adhesive bio-films, spider silk, aggregation of milk proteins, and new functional materials. 65–67 For systems where such particles interact with an oppositely charged surface, the presented results provide insight into the complex dependence of the adsorbed amount on particle size and solution ionic strength.

Supporting Information Available The following files are available free of charge. • SI.pdf: detailed information about QCM-D data analysis and MD simulations; as well as AFM images and additional MC results. • SI.zip: experimental data; simulation data; and Jupyter Notebooks which can be used to reproduce MD and MC simulations; the line segment model; and all presented plots. This material is available free of charge via the Internet at http://pubs.acs.org/.

Acknowledgement We are grateful to Maria Valldeperas Badell for experimental assistance. This work was supported by the Czech Science Foundation (grant 14-12598S), the Ministry of Education, Youth and Sports of the Czech Republic under the National Sustainability Program II. CEITEC 19

ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

2020 (LQ1601), the Swedish Research Council, and the Swedish Strategic Research Foundation. Computational resources were provided by Lunarc in Lund and by the MetaCentrum under the program LM2010005 and the CERIT-SC under the program Centre CERIT Scientific Cloud, part of the Operational Program Research and Development for Innovations, Reg. no. CZ.1.05/3.2.00/08.0144.

References (1) Fändrich, M. On the structural definition of amyloid fibrils and other polypeptide aggregates. Cellular and Molecular Life Sciences 2007, 64, 2066–2078. (2) Pulawski, W.; Ghoshdastider, U.; Andrisano, V.; Filipek, S. Ubiquitous Amyloids. Applied Biochemistry and Biotechnology 2012, 166, 1626–1643. (3) Salvatella, X. Progress in Molecular Biology and Translational Science; Elsevier BV, 2013; pp 73–101. (4) Grey, M.; Linse, S.; Nilsson, H.; Brundin, P.; Sparr, E. Membrane Interaction of αSynuclein in Different Aggregation States. Journal of Parkinson’s Disease 2011, 1, 359–371. (5) Hellstrand, E.; Nowacka, A.; Topgaard, D.; Linse, S.; Sparr, E. Membrane Lipid CoAggregation with α-Synuclein Fibrils. PLoS ONE 2013, 8, e77235. (6) Tycko, R. Physical and structural basis for polymorphism in amyloid fibrils. Protein Science 2014, 23, 1528–1539. (7) Tycko, R. Amyloid Polymorphism: Structural Basis and Neurobiological Relevance. Neuron 2015, 86, 632–645. (8) Williams, T. L.; Serpell, L. C. Membrane and surface interactions of Alzheimer’s Aβ peptide - insights into the mechanism of cytotoxicity. FEBS Journal 2011, 278, 3905– 3917. 20

ACS Paragon Plus Environment

Page 28 of 36

Page 29 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

(9) Jiang, Z.; de Messieres, M.; Lee, J. C. Membrane Remodeling by α-Synuclein and Effects on Amyloid Formation. Journal of the American Chemical Society 2013, 135, 15970–15973. (10) Schladitz, C.; Vieira, E. P.; Hermel, H.; Möhwald, H. Amyloid–β-Sheet Formation at the Air-Water Interface. Biophysical Journal 1999, 77, 3305–3310. (11) Engel, M. F. M.; Khemtemourian, L.; Kleijer, C. C.; Meeldijk, H. J. D.; Jacobs, J.; Verkleij, A. J.; de Kruijff, B.; Killian, J. A.; Hoppener, J. W. M. Membrane damage by human islet amyloid polypeptide through fibril growth at the membrane. Proceedings of the National Academy of Sciences 2008, 105, 6033–6038. (12) Chi, E. Y.; Frey, S. L.; Winans, A.; Lam, K. L. H.; Kjaer, K.; Majewski, J.; Lee, K. Y. C. Amyloid-β Fibrillogenesis Seeded by Interface-Induced Peptide Misfolding and SelfAssembly. Biophysical Journal 2010, 98, 2299–2308. (13) Pfefferkorn, C. M.; Heinrich, F.; Sodt, A. J.; Maltsev, A. S.; Pastor, R. W.; Lee, J. C. Depth of α-Synuclein in a Bilayer Determined by Fluorescence, Neutron Reflectometry, and Computation. Biophysical Journal 2012, 102, 613–621. (14) Campioni, S.; Carret, G.; Jordens, S.; Nicoud, L.; Mezzenga, R.; Riek, R. The Presence of an Air–Water Interface Affects Formation and Elongation of α-Synuclein Fibrils. J. Am. Chem. Soc. 2014, 136, 2866–2875. (15) Walsh, D. M.; Selkoe, D. J. Aβ Oligomers - a decade of discovery. Journal of Neurochemistry 2007, 101, 1172–1184. (16) Xue, W.-F.; Hellewell, A. L.; Hewitt, E. W.; Radford, S. E. Fibril fragmentation in amyloid assembly and cytotoxicity. Prion 2010, 4, 20–25. (17) Winner, B. et al. In vivo demonstration that α-synuclein oligomers are toxic. Proceedings of the National Academy of Sciences 2011, 108, 4194–4199. 21

ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(18) Cohen, S. I. A.; Linse, S.; Luheshi, L. M.; Hellstrand, E.; White, D. A.; Rajah, L.; Otzen, D. E.; Vendruscolo, M.; Dobson, C. M.; Knowles, T. P. J. Proliferation of amyloid- 42 aggregates occurs through a secondary nucleation mechanism. Proceedings of the National Academy of Sciences 2013, 110, 9758–9763. (19) Xue, W.-F.; Hellewell, A. L.; Gosal, W. S.; Homans, S. W.; Hewitt, E. W.; Radford, S. E. Fibril Fragmentation Enhances Amyloid Cytotoxicity. Journal of Biological Chemistry 2009, 284, 34272–34282. (20) Milanesi, L.; Sheynis, T.; Xue, W.-F.; Orlova, E. V.; Hellewell, A. L.; Jelinek, R.; Hewitt, E. W.; Radford, S. E.; Saibil, H. R. Direct three-dimensional visualization of membrane disruption by amyloid fibrils. Proceedings of the National Academy of Sciences 2012, 109, 20455–20460. (21) Monsellier, E.; Bousset, L.; Melki, R. α-Synuclein and huntingtin exon 1 amyloid fibrils bind laterally to the cellular membrane. Scientific Reports 2016, 6 . (22) Höök, F.; Rodahl, M.; Brzezinski, P.; ; Kasemo, B. Energy Dissipation Kinetics for Protein and Antibody-Antigen Adsorption under Shear Oscillation on a Quartz Crystal Microbalance. Langmuir 1998, 14, 729–734. (23) Walsh, D. M.; Thulin, E.; Minogue, A. M.; Gustavsson, N.; Pang, E.; Teplow, D. B.; Linse, S. A facile method for expression and purification of the Alzheimer’s diseaseassociated amyloid β-peptide. FEBS Journal 2009, 276, 1266–1281. (24) Hellstrand, E.; Grey, M.; Ainalem, M.-L.; Ankner, J.; Forsyth, V. T.; Fragneto, G.; Haertlein, M.; Dauvergne, M.-T.; Nilsson, H.; Brundin, P.; Linse, S.; Nylander, T.; Sparr, E. Adsorption of α-Synuclein to Supported Lipid Bilayers: Positioning and Role of Electrostatics. ACS Chem. Neurosci. 2013, 4, 1339–1351. (25) Wang, J.; Cao, Y.; Li, Q.; Liu, L.; Dong, M. Size Effect of Graphene Oxide on Modulating Amyloid Peptide Assembly. Chemistry – A European Journal 2015, 21, 9632–9637. 22

ACS Paragon Plus Environment

Page 30 of 36

Page 31 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

(26) Wang, J.; Zhu, Z.; Bortolini, C.; Hoffmann, S. V.; Amari, A.; Zhang, H. X.; Liu, L.; Dong, M. D. Dimensionality of carbon nanomaterial impacting on the modulation of amyloid peptide assembly. Nanotechnology 2016, 27, 304001. (27) Stenqvist, B.; Thuresson, A.; Kurut, A.; Vácha, R.; Lund, M. Faunus– a flexible framework for Monte Carlo simulation. Molecular Simulation 2013, 39, 1233–1239. (28) Evans, D. F.; Wennerström, H. The Colloidal Domain: Where Physics, Chemistry, Biology, and Technology Meet; Advances in Interfacial Engineering; Wiley-VCH, New York, 1999. (29) Heberle, F. A.; Pan, J.; Standaert, R. F.; Drazba, P.; Ku˘cerka, N.; Katsaras, J. Modelbased approaches for the determination of lipid bilayer structure from small-angle neutron and X-ray scattering data. European Biophysics Journal 2012, 41, 875–890. (30) Wang, F.; Landau, D. P. Determining the density of states for classical statistical models: A random walk algorithm to produce a flat histogram. Physical Review E 2001, 64 . (31) Pasquier, C.; Vazdar, M.; Forsman, J.; Jungwirth, P.; Lund, M. Anomalous ProteinâĂŞProtein Interactions in Multivalent Salt Solution. The Journal of Physical Chemistry B 2017, 121, 3000–3006. (32) Wälti, M. A.; Ravotti, F.; Arai, H.; Glabe, C. G.; Wall, J. S.; Böckmann, A.; Güntert, P.; Meier, B. H.; Riek, R. Atomic-resolution structure of a disease-relevant Aβ(1–42) amyloid fibril. Proceedings of the National Academy of Sciences 2016, 113, E4976–E4984. (33) Colvin, M. T.; Silvers, R.; Ni, Q. Z.; Can, T. V.; Sergeyev, I.; Rosay, M.; Donovan, K. J.; Michael, B.; Wall, J.; Linse, S.; Griffin, R. G. Atomic Resolution Structure of Monomorphic Aβ42 Amyloid Fibrils. Journal of the American Chemical Society 2016, 138, 9663–9674. 23

ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(34) Hellstrand, E. Protein-lipid association and aggregation: from neurodegenerative disease to nanosafety; Doctoral Thesis, Division of Biophysical Chemistry, Lund University: Lund, 2012. (35) Höök, F.; ; Kasemo, B.; Nylander, T.; Fant, C.; Sott, K.; ; Elwing, H. Variations in Coupled Water, Viscoelastic Properties, and Film Thickness of a Mefp-1 Protein Film during Adsorption and Cross-Linking: A Quartz Crystal Microbalance with Dissipation Monitoring, Ellipsometry, and Surface Plasmon Resonance Study. Analytical Chemistry 2001, 73, 5796–5804. (36) Patel, A. R.; Kanazawa, K. K.; Frank, C. W. Antibody Binding to a Tethered Vesicle Assembly Using QCM-D. Analytical Chemistry 2009, 81, 6021–6029. (37) Larsson, C.; Rodahl, M.; Höök, F. Characterization of DNA Immobilization and Subsequent Hybridization on a 2D Arrangement of Streptavidin on a Biotin-Modified Lipid Bilayer Supported on SiO2. Analytical Chemistry 2003, 75, 5080–5087. (38) Ahmed, M.; Davis, J.; Aucoin, D.; Sato, T.; Ahuja, S.; Aimoto, S.; Elliott, J. I.; Nostrand, W. E. V.; Smith, S. O. Structural conversion of neurotoxic amyloid-β1–42 oligomers to fibrils. Nature Structural & Molecular Biology 2010, 17, 561–567. (39) Dubnovitsky, A.; Sandberg, A.; Rahman, M. M.; Benilova, I.; Lendel, C.; Härd, T. Amyloid-β Protofibrils: Size, Morphology and Synaptotoxicity of an Engineered Mimic. PLoS ONE 2013, 8, e66101. (40) Lendel, C.; Bjerring, M.; Dubnovitsky, A.; Kelly, R. T.; Filippov, A.; Antzutkin, O. N.; Nielsen, N. C.; Härd, T. A Hexameric Peptide Barrel as Building Block of Amyloid-β Protofibrils. Angewandte Chemie International Edition 2014, 53, 12756–12760. (41) Luhrs, T.; Ritter, C.; Adrian, M.; Riek-Loher, D.; Bohrmann, B.; Dobeli, H.; Schubert, D.; Riek, R. 3D structure of Alzheimer’s amyloid- (1-42) fibrils. Proceedings of the National Academy of Sciences 2005, 102, 17342–17347. 24

ACS Paragon Plus Environment

Page 32 of 36

Page 33 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

(42) Zheng, J.; Jang, H.; Ma, B.; Tsai, C.-J.; Nussinov, R. Modeling the Alzheimer Aβ17-42 Fibril Architecture: Tight Intermolecular Sheet-Sheet Association and Intramolecular Hydrated Cavities. Biophysical Journal 2007, 93, 3046–3057. (43) Xu, Z.; Paparcone, R.; Buehler, M. J. Alzheimer’s Aβ(1-40) Amyloid Fibrils Feature Size-Dependent Mechanical Properties. Biophysical Journal 2010, 98, 2053–2062. (44) Lu, J.-X.; Qiang, W.; Yau, W.-M.; Schwieters, C. D.; Meredith, S. C.; Tycko, R. Molecular Structure of β-Amyloid Fibrils in Alzheimer’s Disease Brain Tissue. Cell 2013, 154, 1257–1268. (45) Kirkpatrick, S.; Gelatt, C. D.; Vecchi, M. P. Optimization by Simulated Annealing. Science 1983, 220, 671–680. (46) Manning, G. S. Correlation of polymer persistence length with Euler buckling fluctuations. Phys. Rev. A 1986, 34, 4467–4468. (47) Knowles, T. P.; Fitzpatrick, A. W.; Meehan, S.; Mott, H. R.; Vendruscolo, M.; Dobson, C. M.; Welland, M. E. Role of Intermolecular Forces in Defining Material Properties of Protein Nanofibrils. Science 2007, 318, 1900–1903. (48) Arosio, P.; Knowles, T. P. J.; Linse, S. On the lag phase in amyloid fibril formation. Phys. Chem. Chem. Phys. 2015, 17, 7606–7618. (49) Fleck, C.; Netz, R. R.; von Grünberg, H. H. Poisson-Boltzmann theory for membranes with mobile charged lipids and the pH-dependent interaction of a DNA molecule with a membrane. Biophysical Journal 2002, 82, 76–92. (50) Abraham, F. F.; Singh, Y. The structure of a hard-sphere fluid in contact with a soft repulsive wall. The Journal of Chemical Physics 1977, 67, 2384.

25

ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(51) Brenner, S. L.; Parsegian, V. A Physical Method for Deriving the Electrostatic Interaction between Rod-Like Polyions at All Mutual Angles. Biophysical Journal 1974, 14, 327–334. (52) Nordholm, S. Properties of Molecular Fluids in Equilibrium. Lecture notes; University of Gothenburg, 2004. (53) Deserno, M.; Holm, C.; May, S. Fraction of Condensed Counterions around a Charged Rod: Comparison of Poisson-Boltzmann Theory and Computer Simulations. Macromolecules 2000, 33, 199–206. (54) Manning, G. S. Counterion Condensation on Charged Spheres, Cylinders, and Planes. The Journal of Physical Chemistry B 2007, 111, 8554–8559. (55) Cheng, H.; de la Cruz, M. O. Adsorption of rod-like polyelectrolytes onto weakly charged surfaces. The Journal of Chemical Physics 2003, 119, 12635–12644. (56) Campioni, S.; Mannini, B.; Zampagni, M.; Pensalfini, A.; Parrini, C.; Evangelisti, E.; Relini, A.; Stefani, M.; Dobson, C. M.; Cecchi, C.; Chiti, F. A causative link between the structure of aberrant protein oligomers and their toxicity. Nature Chemical Biology 2010, 6, 140–147. (57) Koltover, I.; Wagner, K.; Safinya, C. R. DNA condensation in two dimensions. Proceedings of the National Academy of Sciences 2000, 97, 14046–14051. (58) Dias, R. S.; Pais, A. A. C. C.; Linse, P.; Miguel, M. G.; Lindman, B. Polyion Adsorption onto Catanionic Surfaces. A Monte Carlo Study. The Journal of Physical Chemistry B 2005, 109, 11781–11788. (59) Tzlil, S.; Ben-Shaul, A. Flexible Charged Macromolecules on Mixed Fluid Lipid Membranes: Theory and Monte Carlo Simulations. Biophysical Journal 2005, 89, 2972– 2987. 26

ACS Paragon Plus Environment

Page 34 of 36

Page 35 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Langmuir

(60) Duan, X.; Ding, M.; Zhang, R.; Li, L.; Shi, T.; An, L.; Huang, Q.; Xu, W.-S. Effects of Chain Rigidity on the Adsorption of a Polyelectrolyte Chain on Mixed Lipid Monolayer: A Monte Carlo Study. The Journal of Physical Chemistry B 2015, 119, 6041–6049. (61) Caetano, D. L. Z.; de Carvalho, S. J.; Metzler, R.; Cherstvy, A. G. Critical adsorption of periodic and random polyampholytes onto charged surfaces. Phys. Chem. Chem. Phys. 2017, 19, 23397–23413. (62) Britt, D. W.; Buijs, J.; Hlady, V. Tobacco mosaic virus adsorption on self-assembled and Langmuir–Blodgett monolayers studied by TIRF and SFM. Thin Solid Films 1998, 327-329, 824–828. (63) Ghosh, S. K.; Cherstvy, A. G.; Petrov, E. P.; Metzler, R. Interactions of rod-like particles on responsive elastic sheets. Soft Matter 2016, 12, 7908–7919. (64) Petrova, A. B.; Herold, C.; Petrov, E. P. Conformations and membrane-driven selforganization of rodlike fd virus particles on freestanding lipid membranes. Soft Matter 2017, 13, 7172–7187. (65) Fowler, D. M.; Koulov, A. V.; Balch, W. E.; Kelly, J. W. Functional amyloid - from bacteria to humans. Trends in Biochemical Sciences 2007, 32, 217–224. (66) Slotta, U.; Hess, S.; Spieß, K.; Stromer, T.; Serpell, L.; Scheibel, T. Spider Silk and Amyloid Fibrils: A Structural Comparison. Macromolecular Bioscience 2007, 7, 183– 188. (67) Larsen, P.; Nielsen, J. L.; Dueholm, M. S.; Wetzel, R.; Otzen, D.; Nielsen, P. H. Amyloid adhesins are abundant in natural biofilms. Environmental Microbiology 2007, 9, 3077– 3090.

27

ACS Paragon Plus Environment

Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Graphical TOC Entry 64

16 8

32

2

128 coar segr ai ned model

l i nesegment model

28

ACS Paragon Plus Environment

Page 36 of 36