The Aquila Digital Community The Aquila Digital Community

As the primary toxic species in the etiology of Alzheimer disease (AD) are low molecular weight oligomers of A β , it is crucial to understand the structure of A β oligomers for gaining molecular insights into AD pathology. We have earlier demonstrated that in the presence of fatty acids, A β 42 peptides assemble as 12-24mer oligomers. These Large Fatty Acid-derived Oligomers (LFAOs) exist predominantly as 12mers at low and as 24mers at high concentrations. The 12mers are more neurotoxic than the 24mers and undergo self-replication, while the latter propagate to morphologically distinct ﬁbrils with succinct pathological consequences. In order to glean into their functional differences and similarities, we have determined their structures in greater detail by combining molecular dynamic simulations with biophysical measurements. We conjecture that the LFAO are made of A β units in an S-shaped conformation, with the 12mers forming a double-layered hexamer ring (6 × 2) while the structure of 24mers is a double-layered dodecamer ring (12 × 2). A closer inspection of the (6 × 2) and (12 × 2) structures reveals a concentration and pH dependent molecular reorganization in the assembly of 12 to 24mers, which seems to be the underlying mechanism for the observed biophysical and cellular properties of LFAOs.


I. INTRODUCTION
One of the hallmarks of Alzheimer disease (AD) pathology is the deposition of amyloid-β (Aβ) peptide fibrillar aggregates as plaques in brains of patients.The neuronal loss, however, seems to be triggered primarily by low-molecular weight (LMW) oligomers that are formed earlier than the high-molecular weight fibrils during the aggregation process. 1 Therefore, there is a growing interest in isolating LMW oligomers and deriving their atomistic structures and dynamics.However, the transient nature and heterogeneity of the oligomers makes their isolation and characterization, either from endogenous or exogenous sources, difficult.As a consequence, deriving structural models, as needed for understanding their toxicity mechanism and mode of propagation, poses a challenge.
We have developed a method for generating distinct Aβ42 12-24mer assemblies, called large fatty acid-derived oligomers (LFAOs), [2][3][4][5][6][7] as their generation is catalyzed by saturated fatty acids that are toxic to neuroblastoma cells.At higher concentrations, LFAOs convert from the 12mer species to more disperse distribution of 12-24mer oligomers. 2This concentration-dependent transition is significant because the 12mers self-replicate in the presence of monomers and are more apoptotic to neuroblastoma cells than the 24mers. 2n the other hand, the 24mers faithfully propagate towards morphologically-unique fibrils and induce acute cerebral amyloid angiopathy (CAA) in transgenic mice. 3Therefore, it is imperative to obtain atomistic insights into the differences between the 12mer and 24mer LFAOs, as well as the transition from one form to the other, in order to better understand their unique properties.Unfortunately, due to some of the The Journal of Chemical Physics ARTICLE scitation.org/journal/jcpaforementioned reasons, the structures of these LFAOs have so far not been resolved.Guided by computational investigations of the Hansmann laboratory, we present in this paper structural models for 12mer and 24mer LFAOs and a mechanism for the transition from one to the other assembly.Derivation of these models relies on the computational techniques developed in the Hansmann lab.These methods have been already used in previous work to study the various polymorphs seen in amyloids [8][9][10][11] and to probe the factors that modulate conformational switching in amyloids. 11Combining computational investigations with novel biophysical experiments, we conclude that the 12mers prefer to form two-layered rings, each ring a hexamer (6 × 2), while 24mers transition to another species of two-layered assemblies, here each ring a dodecamer (12 × 2).We also eliminate the possibility of a single dodecamer ring (12 × 1) structure for 12mers and of four stacked hexamer rings (6 × 4) as the structure for 24mers.Our models allow us to explain the experimentally observed conformational-dynamics of LFAOs to identify the key residues involved in conformational switching, and provide hints at the structural basis for the different pathogenies of LFAO 12mers and 24mers.

B. Molecular dynamics simulation
The oligomer models presented in this article are based on previous work where we constructed a series of N-fold ring-like Aβ42 oligomer models, 8 including six-fold and twelve-fold models of Aβ (11-42), that are characterized by S-shaped chain configurations forming three β-strands liked by two turn regions.The chains in a ring are kept together by hydrophobic contacts in the region of residues 20-28 and an inter-chain salt bridge K16-D23.These models are the starting point in our construction of the (6 × 2), (6 × 4), and (12 × 2) models described later.
Our simulations rely on the software package GROMACS (GROningen MAchine for Chemical Simulations) 5.1.5 12and employ the AMBER ff99SB-ildn 13 force field for proteins and TIP3P water 14 as solvent, a choice also employed by us in our previous work. 15The temperature of 300 K and a pressure of 1 bar are controlled by v-rescale thermostat 16 and Parrinello-Rahman barostat. 17The bond-lengths are restrained by the LINCE algorithm 18 and the SETTLE algorithm 19 allowing us to use a time step of 2 fs for integration.The protein and solvent are put into a box with the side length and periodic boundary conditions of 13.28 nm (for 6 × 2), 18.86 nm (for 12 × 2), and 13.31 nm (for 6 × 4), and electrostatic interactions are calculated by the particle mesh Ewald (PME) method. 20Stability of our structures is probed by following our molecular dynamics trajectories over 20 ns, with only the last 10 ns used for analysis.
Most of our analysis is carried out within the tool set provided by GROMACS, with snapshots of configurations visualized by VMD (Visual Molecular Dynamics). 21The distance between residues are defined as the average distance between heavy atoms in the side chains of each residue; for example, the NZ atom on K16 and the OD1 atom on D23 are used to calculate the inter-chain salt-bridge between residues K16 and D23.We measure the solvent accessible surface area (SASA) by both the g_sasa and POPS 22 software tools, the later allowing one to separate the hydrophobic and hydrophilic areas.Since the two methods use different definitions of surface area, values may differ slightly.The binding energy are approximated with MMPBSA.py in AmberTools, 23 with the setting igb = 8 for the GBSA (Generalized Born and Surface Area continuum solvation) part.

C. Aβ monomer and oligomer purification
Aβ monomers and oligomers were purified as described previously. 4Briefly, 0.5-1 mg of the synthetic peptide was weighed into a sterile microcentrifuge tube and resuspended in 500 µL of 10 mM NaOH.After incubation at 25 • C for 30 min, the sample was loaded onto a Superdex-75 HR 10/30 size exclusion column pre-equilibrated in 20 mM Tris, pH 8.0 using either an AKTA FPLC (GE Healthcare) or a BioLogic DuoFlow (BioRad) purification system.Fractions of 500 µL were collected at a constant flow rate of 0.5 mL/min.Aβ concentrations were determined using intrinsic tyrosine absorbance (ε = 1450 cm −1 M −1 at 276 nm) on a Cary 50 Ultra violet-Visible (UV-Vis) spectrometer (Agilent Technologies).To generate LFAOs, Aβ monomers (50-60 µM) were incubated with 5 mM C12:0 NEFA and 50 mM NaCl at 37 • C for 48 h.LFAOs were then centrifuged at 20 000 g for 20 min before being purified via size exclusion chromatography (SEC) as described above.

D. ANS binding assay
ANS binding experiments were done as described previously. 2 In a set of separate experiments, varying pH or ionic concentration, LFAOs were exchanged into 20 mM Tris at the appropriate pH using a 3.5 kDa molecular weight cut-off Slide-A-Lyzer MINI Dialysis Device (Thermo Fisher Scientific) following the manufactures' protocol.Upon the addition of 100 µM ANS followed by 1 min of equilibration, fluorescence measurements of LFAOs (8, 6, 4, 2, 1, or 0.5 µM) were collected on a Cary Eclipse instrument (Agilent Technologies) by scanning the emission spectrum between 400 and 650 nm upon excitation at 388 nm.The area under the curve for each respective pH or NaCl titration experiment were then normalized and plotted as shown.The data presented are representative of three independent experiments for each specific pH or NaCl concentration.

E. Circular dichroism spectroscopy
Data were collected on a Jasco J-815 spectropolarimeter attached with a Peltier temperature controller.To a solution of LFAOs (1 or 8 µM), sodium dodecyl sulfate (SDS) was added to a final concentration of 1% (wt/v), followed by heating from 10 to 90 • C (ramp rate = 0.5 • C/min) while monitoring the signal at 206 nm every 1 min (hold time = 30 s, DIT = 32 s, bandwidth = 5 nm).The data was processed by normalizing from 0 to 1 [as shown in Fig. 6(a)], followed by determining the apparent standard free energy using the Gibbs-Helmholtz equation [as shown in Fig. 6(b)].The data presented are representative of three independent experiments.

F. Atomic force microscopy
Atomic force microscopy (AFM) samples were prepared following a previously published procedure. 3Freshly cleaved mica substrates were first treated with 150 µL of APTES (3-aminopropyltriethoxysilane) solution (500 µl in 50 ml of 1 mM acetic acid) for 20 min.The APTES solution was then decanted and rinsed three times with 150 µL DI H 2 O.The substrates were dried under a stream of N 2 and stored in the desiccator for 1 hour.A 150 µL aliquot of the Aβ solution (either 1 or 5 µM in 20 mM Tris-HCl, pH 8.0) was deposited onto the amine-treated mica substrates for 30 min to adsorb the proteins.The Aβ solution was then decanted, and the samples were rinsed three times with 150 µL DI H 2 O.The samples were dried under a stream of N 2 and stored in the desiccator until imaging.
AFM analysis of LFAOs was conducted using a Dimension Icon atomic force microscope (Bruker) in Peak-Force Tapping mode.AFM scanning was performed using NanoScope 8.15r3sr8 software and the images were analyzed in NanoScope Analysis 1.50 software.Imaging was performed using a sharp silicon nitride cantilever (SNL-C, nominal tip radius of 2 nm; nominal resonance frequency of 56 kHz; nominal spring constant of 0.24 N/m) in a standard probe holder under ambient conditions with 512 × 512 data point resolution.The AFM height image was deconvoluted and corrected for diameter manually by calculating the dead space between the base and the sample created by the curvature the cantilever tip.

A. LFAOs are two-layered rings
[4][5]7 LFAOs display the presence of two aggregate distributions corresponding to sedimentation coefficients 5S and 7S [Fig. 1 In a recent paper, 8 we have shown that unlike the more common, but less toxic, Aβ40, Aβ42 polypeptide innately can form pore-like trimers and larger oligomers.This is because Aβ42 chains are able to assume a S-shaped three-stranded motif, 24 while Aβ40 peptides are not stable in this form and instead take U-shaped conformations in fibrils. 15Building on our previously presented models 8 and guided by the experimental size measurements, we conjecture that the 12mer is organized as two ring-like hexamers stacked on top of one another (6 × 2) [Fig.2(a)]. 8In such an arrangement, the hydrated diameter of the oligomer measures 14-15 nm with a height of 3-4 nm, giving a flattened disc-like appearance.6][27] In order to model the heavier LFAO 24mers, two possible models were considered: a tetramer of hexamers (6 × 4) or a dimer of dodecamers (12 × 2).Both models are consistent with previous experimental results in their dimensions and size as seen by AFM, which range between 14 and 20 nm of the height observed 2,3,6  Our earlier AFM analysis on LFAOs showed a bimodal distribution of oligomers [Fig.1(c)], but we did not interrogate the morphology of individual aggregates.In this study, we used 1 and 5-8 µM LFAO samples to generate 12 and 12-24mers, respectively, throughout this manuscript.For detailed AFM analysis, we have employed a sharp silicon nitride cantilever with a lower resonant frequency and smaller spring constant to obtain higher resolution images.Height and phase images were obtained over small scan areas containing LFAOs of different diameters (Fig. 3).Height images indicate spherical particles with average diameters of 20 nm for the 12mers and average diameters of 40 nm for the 24mers.Average height for the 12mers is 2.5-3.5 nm and that of the 24mers is 5-8 nm.Phase images for both oligomers show a characteristic "donut" shape, indicative of differences in modulus and/or adhesion from the outer to the inner edge of the material, which we attribute to the existence of a cavity.Similar height and phase images have been reported for hollow nanoparticles. 28o see the correspondence between the morphology obtained from AFM and the simulated models, the height analysis and various structural models were overlaid (Fig. 4 27 Note that our structural models are for residues 11-42 in the Aβ42 chains as residues 1-10 are flexible in all resolved fibril structures.To evaluate the effect of first 10 residues, we also constructed the full-length Aβ42 oligomers for (6 × 2) and (12 × 2) models with the N-terminal segment in the beginning in a random coil configuration that is allowed to relax in a molecular dynamics simulation over 10 ns.The orientation of the first ten residues stayed random, and their sole effect was that the diameter of (6 × 2) oligomers changed from 10.5 to 14.6 nm.
In order to understand why the heavier LFAO 24mers appear to form a dimer of dodecamers (12 × 2) instead of a tetramer of hexamers (6 × 4), we have simulated all three Aβ11-42 oligomer models, the (6 × 2) 12mer and the two 24mer models (6 × 4) and (12 × 2), by atomistic molecular dynamics (MD).As we are neither modeling the association into the 12mers nor the transition between 12mers and 24mers, only relative short simulations are needed to explore energetics and stability of these models.Note that the experimental measurements were obtained at a pH = 8, i.e., under neutral/alkaline conditions, which were modeled in our simulations by changing the H13 and H14 residues into a deprotonated state (named by us the HIE state).Each system is followed in two independent runs.In Table I, we list the solvent accessible surface area (SASA) and the binding energy as approximated by the Molecular Mechanics/Generalized Born and Surface Area (MMGBSA) continuum solvation approach.This approximation is justified because we are not interested in absolute values for the binding energies of the three models but only in the qualitative differences between them.
Comparing the SASA of the three models, we note that the values for the (12 × 2) model of the 24mers is about double that of the (6 × 2) 12mer, while the corresponding SASA of the (6 × 4) 24mer model is less than 1.5 times that of the 12mer.Hence, based on the ANS binding data, which showed that the solvent accessible hydrophobic surface area doubles from LFAO 12mers to 24mers [Fig.1(e)].We conclude that the (12 × 2) model is the more likely structure for the 24mers.This model has a diameter of 15.7 nm, which extends to 19.5 nm for the full-sized model (including the first ten N-terminal residues), also in agreement with the experimentally measured dimensions.
The (12 × 2) model for the 24mer is also favored by the binding energies shown in Table I.The (12 × 2) structure has binding energies more favorable than the binding energies of where A is interaction between chains within a ring, B the interaction between chains of neighboring rings, and C the interaction of chains with surrounding water (i.e., proportional to exposed surface).With the same definitions, one finds that the binding energy of a ( 12  The solvent accessible surface area (SASA) and the binding energy of the three considered oligomer models under neutral pH conditions.Shown are for all quantities the averages (Avg) as obtained from two runs of 20 ns and their standard deviations (StD).The binding energy was calculated in a MMGBSA approximation and is composed by four terms: the van der Waals energy E vdW , the electrostatic energy E elect , a generalized Born approximation of the solvation energy E GB , and solvent surface tension interaction term E surf .The SASA values are calculated by two different tools: g_sasa and POPS.In POPS, the SASA values can be further separated into hydrophobic and hydrophilic contributions.The binding energy are defined as follows: The histidine side chains at positions 13 and 14 orient themselves towards one side of the hexamer ring in the upper layer, while similar histidines on the congruent Aβ42 unit on the lower layer hexamer ring do so on the opposite side.Hence, the positioning of the histidine chains hinders further association of hexamer rings (leading to steric clashes), which would make a (6 × 4) structure thermodynamically expensive.
The above results suggest that the transition between 12mers and 24mers is not a simple stacking of two 12mers.As the repulsive interaction between (6 × 2) 12mer prohibits stacking and leads to a large unfavorable binding energy, a reorganization of chains is needed, leading to a (12 × 2) structure where the histidines (H13 and H14) are moved ∼95 • perpendicular to the axis of the oligomer, thus preventing potential charge repulsion due to protonation.Furthermore, this transition also results in the exposure of hydrophobic residues along either side of oligomer face, see Table I, evident from the increase in ANS binding [Fig.1(e)].
This scenario is also consistent with the observation that the transition between 12mer and 24mer depends on concentration, see [Figs.1(b) and 1(e)], 2 and can be understood from the energetics of our (6 × 2) 12mer and (12 × 2) 24mer models.At low concentrations, the 12mers are separated and experience little repulsive forces.At the same time, the gain in energy from forming a 24mer is negligible, and therefore the equilibrium between the two forms shifts towards the 12mer.On the other hand, above a critical threshold, the distance between (6 × 2) 12mer rings becomes so close that the repulsive interaction between them becomes noticeable.At that point, it becomes energetically more favorable to rearrange and form (12 × 2) assemblies, i.e., the equilibrium is shifted toward the 24mers.

B. Dynamics based on stability
While the above reasoning is plausible, more evidence is needed to support our models for 12mer and 24mer and for the implied mechanism of the concentration-dependent transition between the two forms.Since the transition from 12 to 24mer involves in our model exposure of charged histidines, we reason that the transition is caused by electrostatic interactions, and that it should be modulated by pH and ion concentration in a way that can be understood from our proposed models.In an effort to obtain molecular details on such a possibility, we investigated the effect of buffer pH and salinity on 12 to 24mer transition using the ANS binding assay.This transition was previously established at pH 8.0 in low ionic strength conditions (absence of salt) [ ; Figs. 5(a .Specifically, this consolidates the idea derived from the structures regarding the involvement of protonation/deprotonation events in such a transition.To further investigate differences in 12mers and 24mers, thermodynamic stability analysis was performed in the presence of a denaturant (sodium dodecyl sulfate, SDS).In these experiments, the conversion of β-sheet structure adopted by LFAOs to an α-helix in the presence of SDS was monitored by far-UV CD at 206 nm as a function of temperature [Fig.6(a)].

TABLE II.
The solvent accessible surface area (SASA) and the binding energy of the three considered oligomer models under neutral pH conditions.Shown are for all quantities the averages (Avg) as obtained from two runs of 20 ns and their standard deviations (StD).The binding energy was calculated in a MMGBSA approximation and is composed by four terms: the van der Waals energy E vdW , the electrostatic energy E elect , a generalized Born approximation of the solvation energy E GB , and solvent surface tension interaction term E surf .The SASA values are calculated by two different tools: g_sasa and POPS.In POPS, the SASA values can be further separated into hydrophobic and hydrophilic contributions.The binding energy are as follows: The Journal of Chemical Physics

ARTICLE scitation.org/journal/jcp
The equilibrium data were processed (as described in the experimental) to obtain the apparent Gibbs free energy at 37 • C, which were found to be −0.879 and −0.489 kcal/mol for LFAO 12mer (1 µM) and 24mer (8 µM), respectively [Fig.6(b)].This difference of 0.4 kcal/mol is minimal and insignificant, revealing that both isoforms have similar thermodynamic stability.
In our simulations, we model the switch from neutral to acidic conditions by repeating our simulations of the three systems (the (6 × 2) for the 12mers, and the (6 × 4) and ( 12× 2) model for the 24mers), with protonated histidine 13 and 14, i.e., setting the charge state to positive (called HIP in this paper).All other parameters in the two sets of simulations are identical.Data for the pH < 7 simulations are summarized in Table II.Comparing binding energies between charged (pH < 7) and neutral forms, it appears that the energy differences between 24mer (12 × 2) and 12mer (6 × 2) favor the 24mers in the charged state more than in the neutral state [(∆∆G) = ∼13 kcal/mol].This is unexpected as the experimental data shows a faster transition toward 24 mer at pH ≥ 7. Hence, the sharper transition between 12mer and 24mer is not because at pH ≥ 7 the 24mer is energetically more favored over the 12mer than in acidic conditions.Instead, the sharper transition at pH ≥ 7 is because the repulsion between 6 × 2 dodecamers is much higher for the charged forms than the neutral forms (Tables I and II).In other words, the repulsion between the 12 mers is larger at low pH values than in the neutral range where a faster 12 to 24mer transition was observed.Hence, we conjecture that the less pronounced transition at low pH is because larger concentrations are needed to overcome the stronger repulsion between the 12mer than at neutral or higher pH.Note also that while the solvent accessible surface area does not differ between the charged and neutral forms in the (6 × 2) model, a difference in SASA is observed for the (12 × 2) structure.For neutral conditions (HIE), the (12 × 2) structure of a 24mer exposes roughly two times (332 Å 2 ; Table I) the hydrophobic surfaces of the (6 × 2) 12mer (164 Å 2 ; Table I), an observation that the experimental results concur with Fig. 1(e).On the other hand, under acidic conditions (HIP), the 24mer exposes with 347 Å 2 (Table II) more than double the hydrophobic surface of two 12mers (164 Å 2 ; Table II), making formation of 24mer less favorable), which further agrees with the experimentalZ observations.

IV. DISCUSSION
What is the cause for the above described changes in the way the 12mer and 24mer differ at low pH and neutral pH, changes that in turn modulate the concentration-dependent transition between 12mer and 24mer?These changes have to be connected with differences in the two geometries and with the way these differences are modulated by pH.One possibility is the way residues are exposed to the solvent.Figure 7 shows the per-residue differences in the solvent exposed surface area between HIE and HIP states for (6 × 2) and (12 × 2) structures.While on average the SASA does not differ between charged and neutral forms in (6 × 2) 12mer, the turn region t1 between β1 and β2 (see Fig. 2) is more exposed under neutral conditions (HIE) than under acidic conditions (HIP), while most other residues are more exposed in the HIP state than in the HIE state.A similar picture is seen for the (12 × 2) 24mer, only that here the difference in exposure of residues in the turn region is much smaller than for the 12mer; while the overall exposure of surface to the solvent is larger under acidic conditions (HIP) than under neutral conditions (HIE).Hence, our solvent accessible surface differences indicate that the pH-modulation of the transition between 12mer and 24mer involves this turn region.
This observation is confirmed by Fig. 8, where the perresidue contributions to the binding energy are shown.The contribution of binding energy for each residue was calculated in three different ways: first, the difference between HIE and HIP states for (6 × 2) and (12 × 2) models [Fig.8 and other residues located on the packing surface that directly interact with residues on the neighboring fold.For this segment, the binding energy contribution is for neutral pH similar between two isolated (6 × 2) rings and when fully associated as a (6 × 4) assembly, while under acidic conditions the binding energy contribution from these segments favor isolated (6 × 2) assemblies.This is consistent with Fig. 5(a) which shows that under acidic conditions, the residues in this segment are less exposed to solvent than under neutral conditions.Visual inspection of this segment in the (6 × 2) model for both HIE and HIP systems shows that this region is slightly more distorted in HIP than in HIE (Fig. 9).This distortion is related with (and can be quantified by) a weakening of the saltbridge between K16 and E22 or D23 (mainly the K16-D23, see Ref. 29) that is formed between neighboring chains (Fig. 9).In order to demonstrate this point, we have calculated the average distance of all corresponding salt-bridge pairs between the O δ1 /O δ2 atoms on D23 and N ζ atom on K16, and the distributions of such distances for the HIP and HIE models are shown in Fig. 10.
Compared with the neutral state (HIE), the acid state decreases the stability of the inter-chain salt-bridge in the (6 × 2) model.How do the different charge states of the histidine lead to this effect?In HIP, residues 13 and 16 are both positively charged and the repulsive interaction between them distorts the geometry of the salt-bridge between residue K16 and either D23 or E22.In order to quantify the repulsive interaction between the two positively charged H13 and K16, we have calculated the average distance between the mass center of imidazole on H13 and N ζ atoms on K16 on the same chain and drawn the distribution of distances between H13 and K16 in Fig. 10.The geometry in the 24mer (a two-layer dodecamer ring) is such that the different charge states of H13 do not change the average distance between H13 and K16 and therefore also does not weaken the salt bridge K16-D23 (E22), see Fig. 9. On the other hand, for the 12mer (a two-layer hexamer ring), the distribution of distances between residues H13 and K16 is shifted towards larger values for acidic conditions reflecting the repulsive interaction between the two residues under these conditions.This is not seen for the 24mer because for a (12 × 2) structure, the same histidines (H13 and H14) are moved ∼95 • perpendicular to the axis of oligomer, and thus preventing potential charge repulsion due to protonation.The net-effect of this repulsive interaction between the charged histidine and the lysine K16 is a weakening of the salt bridge K16-D23 (E22) (Fig. 10), which in turn reduces the stability of the turn region between the β1 and β2 strands and the hydrophobic core region for the peptides in the 12mer under acidic conditions.The binding energies [Fig.8(c)] indicate that this distortion leads to a larger repulsive interaction between the (6 × 2) structure, shifting the equilibrium toward the 12mer.

V. CONCLUSIONS
Combining a variety of biophysical measurements and molecular dynamics simulations, we put forward models for the fatty acid catalyzed Aβ42 assemblies called LFAOs.The experimentally observed existence of 12 and 24mer LFAOs in a concentration-dependent manner agrees with the models generated by simulations.First, ring-like two-layered assemblies were derived for both 12 and 24mers (two stacked hexamer rings in the case of the 12mer, and two stacked dodecamer rings in the case of the 24mer).The diameter and height of the two oligomers agree well with their morphology in AFM images that appear as punctate spherical particles.More importantly, the most noticeable feature of the oligomer structure, the ring structure with a cavity in the middle, is also in agreement with the phase changes observed in AFM images.Furthermore, the determination that protonation of histidine side chains (H13 and H14) are the key in determining the structural conversion of (6 × 2) to (12 × 2) is also supported by the experiments on both wild-type and double alanine mutant of the histidines.Finally, the energetic contributions calculated from the models are also in agreement with the free energy changes observed experimentally.Together, both simulations and experiments point out to the proposed models for the LFAO structure and dynamics that remained elusive thus far.6][27] Based on AFM imaging, ring-like, spherical low molecular oligomers are observed to be transiently formed before the formation of high molecular weight oligomers, which then laterally associate to form protofibrils. 25 Economou and colleagues observed that even at low concentrations, Aβ42 but not Aβ40, form ring-like hexamers that convert to dodecamers, which consequently seed protofibril formation. 26Our observations on the LFAO structure and propagation support our previous reports and add much needed detail.The studies undertaken in our labs have confirmed that LFAOs are ring-like dodecamers at low concentrations. 6,8Investigations on the mechanism of LFAO The Journal of Chemical Physics ARTICLE scitation.org/journal/jcppropagation suggested that LFAO 24mers are formed at higher concentrations, which grow larger to form a key intermediate on route to fibril formation in a three-step mechanism. 4The (12 × 2) structure observed for 24mers explains the fact that 24mers, and not (6 × 2) 12mers, are able to faithfully propagate by associating with one another mediated by increased hydrophobic surface interactions with 24mer units of LFAOs as observed previously. 25,27Overall, these observations are in agreement with those described previously.Perhaps the intriguing and enigmatic properties of LFAOs are (i) the ability of 12mers to self-replicate in the presence of monomers, (ii) to convert to 24mers in a concentration dependent manner, and (iii) the striking differences in the pathogenicity of 12 and 24mers.While both oligomeric forms are pathogenic, LFAO 12mers are more apoptotic to neuroblastoma cells than the 24mers. 2LFAOs also induce acute CAA in transgenic mice brains selectively, although it remains unclear which form of LFAO is responsible for this phenotype. 3The results presented here bring out the molecular signatures that are responsible for the structural and functional differences between LFAO 12 and 24mers and provide insights into the structure and mechanism by which LFAOs behave and become neurotoxic.At elevated concentrations, (6 × 2) LFAOs form 24mers by adopting a (12 × 2) structure accompanied by reorganization in the assembly, which exposes the hydrophobic residues along either side of oligomer face, an observation also supported by an increase in ANS binding.Such a reorganization increases the susceptibility of 24mers for further oligomer associations mediated largely by hydrophobic interactions.This is indeed supported by the fact that LFAO 24mers are able to propagate morphologically distinct fibrils made by repeats of LFAO units. 4nsights into the molecular underpinnings of oligomer behavior is much needed to understand AD pathology and for future therapeutic interventions.This work is a step toward advancing our knowledge into this critical area, providing insights into the structure and mechanism by which LFAOs behave and become neurotoxic.

FIG. 4 .
FIG. 2. Models of LFAO 12 and 24mers.(a)The backbone ribbon representation of our (6 × 2) model for LFAO 12mers built from two hexameric rings stacked on top of one another.β-sheets and turns are indicated as β and t, respectively.(b) and (c) Two possible models for LFAO 24mer assembly with (6 × 4) and (12 × 2), respectively.The structure of the individual chains is derived from the ones the chains take in the fibril model24 as adapted from Ref. 8.

FIG. 9 .FIG. 10 .
FIG. 9. Detailed structure of LFAO 12mers.(a) Final structures of the (6 × 2) model for neutral pH (HIE, blue) and acidic conditions (HIP, red) states after a 20 ns molecular dynamics trajectory.The side chains of residues 20 to 28, where the contributions to binding energy differed mostly with pH, are shown in bond representation.The histidine residues H13/H14 are colored in green (HIE) and mauve (HIP).(b) a close-up view of the stabilizing salt bridge interactions between D23 and K16 the interface of two Aβ units are shown along with the histidines 13 and 14.
Table I], while the binding energy of the 6 × 4 model is substantially higher than that of two isolated 6 × 2 models [see the difference ∆(B) = E_(6 × 4) − 2 × E_(6 × 2) in Table I], disadvantaging formation of this assembly.In order to show why the 12 × 2 model is more favorable than the (6 × 4) model, we neglect entropic contributions and approximate the binding energy of a 12mer by