In silico investigation and targeting of amyloid β oligomers of different size
Click here to load reader
Transcript of In silico investigation and targeting of amyloid β oligomers of different size
Registered Charity Number 207890
Accepted Manuscript
This is an Accepted Manuscript, which has been through the RSC Publishing peer review process and has been accepted for publication.
Accepted Manuscripts are published online shortly after acceptance, which is prior to technical editing, formatting and proof reading. This free service from RSC Publishing allows authors to make their results available to the community, in citable form, before publication of the edited article. This Accepted Manuscript will be replaced by the edited and formatted Advance Article as soon as this is available.
To cite this manuscript please use its permanent Digital Object Identifier (DOI®), which is identical for all formats of publication.
More information about Accepted Manuscripts can be found in the Information for Authors.
Please note that technical editing may introduce minor changes to the text and/or graphics contained in the manuscript submitted by the author(s) which may alter content, and that the standard Terms & Conditions and the ethical guidelines that apply to the journal are still applicable. In no event shall the RSC be held responsible for any errors or omissions in these Accepted Manuscript manuscripts or any consequences arising from the use of any information contained in them.
www.rsc.org/molecularbiosystems
ISSN 1742-206X
1742-206X(2010)6:1;1-G
Top Quality Bioscience Journals
Take a look today!
www.rsc.org/journalsRegistered Charity Number 207890
New for 2010
MedChemComm - focusing on medicinal chemistry research, including new studies related to
biologically-active chemical or biochemical entities that can act as pharmacological agents with therapeutic
potential or relevance. www.rsc.org/medchemcomm
New for 2009
Integrative Biology - a journal focusing on quantitative multi-scale biology using enabling technologies and tools to exploit the convergence of biology with physics, chemistry, engineering, imaging and informatics. www.rsc.org/ibiology
Metallomics - a journal covering the research fields related to metals in biological, environmental and clinical systems. www.rsc.org/metallomics
Molecular BioSystems - a journal with a focus on the interface between chemistry and the -omic sciences and systems biology. www.molecularbiosystems.org
Organic & Biomolecular Chemistry - an international journal covering the breadth of synthetic, physical and biomolecular chemistry. www.rsc.org/obc
Natural Product Reports (NPR) - a critical review journal which stimulates progress in all areas of natural products research. www.rsc.org/npr
Photochemical & Photobiological Sciences - publishing high quality research on all aspects of photochemistry and photobiology, encouraging synergism between the two areas. www.rsc.org/pps
Volume 6 | N
umber 1 | 2010
Molecular B
ioSystems
Pages 1–276
www.molecularbiosystems.org Volume 6 | Number 1 | January 2010 | Pages 1–276
PAPERDieter Willbold et al.Competitively selected protein ligands pay their increase in specificity by a decrease in affinity
METHODThomas Kodadek et al.Rapid identification of orexin receptor binding ligands using cell-based screening accelerated with magnetic beads
Indexed in
MEDLINE!
mbs006001_Cover.indd 5-1 26/11/2009 15:05:22
MolecularBiosystems
View Article OnlineView Journal
Structure and dynamics of Aβ-fibril models used as targets to investigate the effect of putative inhibitors at
different stages of Aβ−related disease
Page 1 of 19 Molecular BioSystems
Mo
lecu
lar
Bio
Sys
tem
s A
ccep
ted
Man
usc
rip
t
Dow
nloa
ded
by M
ount
Alli
son
Uni
vers
ity o
n 03
/05/
2013
11:
21:4
1.
Publ
ishe
d on
23
Apr
il 20
13 o
n ht
tp://
pubs
.rsc
.org
| do
i:10.
1039
/C3M
B70
086K
View Article Online
In silico investigation and targeting of amyloid β
oligomers of different size
Ida Autiero,a Michele Savianob* and Emma Langellaa*
*Correspondence to: Michele Saviano, National Research Council, Institute of
Crystallography, 70126 Bari, Italy. E-mail: [email protected]. Emma Langella,
National Research Council, Institute of Biostructures and Bioimaging, 80138 Naples,
Italy. E-mail: [email protected]
a National Research Council, Institute of Biostructures and Bioimaging, 80138 Naples,
Italy
b National Research Council, Institute of Crystallography, 70126 Bari, Italy
1
2
3
4
5
6
7
8
9
10
Page 2 of 19Molecular BioSystems
Mo
lecu
lar
Bio
Sys
tem
s A
ccep
ted
Man
usc
rip
t
Dow
nloa
ded
by M
ount
Alli
son
Uni
vers
ity o
n 03
/05/
2013
11:
21:4
1.
Publ
ishe
d on
23
Apr
il 20
13 o
n ht
tp://
pubs
.rsc
.org
| do
i:10.
1039
/C3M
B70
086K
View Article Online
ABSTRACT
Aggregation of Amyloid β (Aβ) peptide into fibrils has been implicated in the
pathogenesis of Alzheimer's disease (AD). As a result, in recent years, substantial
efforts have been expended in the study of the mechanism of aggregation of the Aβ
peptide as well as of its inhibition by potential drug molecules. In this context, we have
built a model of the Aβ (17-42) deca-oligomer using the solid-state NMR (ssNMR)
structure of the Aβ (17-42) penta-oligomer as a reference. Both the penta- and deca-
oligomer systems have been studied by all-atom molecular dynamics (MD) simulations
and used as target systems for the investigation of the mechanism of action of a
trehalose-derived Aβ aggregation inhibitor. In the deca-oligomer all the main structural
features of the putative fibrillar state are retained. Moreover, the simulations reveal a
remarkable gain in stability as the oligomer grows. MD studies of the inhibitor in
complex with the penta- and deca-oligomers indicate a significant destabilization of
the structure beyond the hampering of the addition of successive Aβ peptides at the
ends of the fibril due to the presence of the inhibitor molecule. Our work provides an
easy and effective approach which could be useful for the in silico development of
potential drug molecules acting at different stages of the progression of Aβ–related
disease.
INTRODUCTION
Alzheimer's disease (AD) is currently one of the most common and devastating forms
of dementia correlated with β-amyloid peptides (Aβ) accumulation in human brain
tissue (1-3). These peptides are proteolytic byproducts of the Aβ protein precursor and
are most commonly composed of 40 (Aβ1-40) and 42 (Aβ1-42) amino acids (4-11). Aβ
peptides easily aggregate to form fibrils (12-14) which finally accumulate into plaques.
Important insights into the structural characteristics of (Aβ1-40) as well as Aβ(1–42)
fibrils have been determined, establishing that amyloid fibrils show a cross-β structure
(9) that contains parallel, in-register β-sheets (15-16). The Aβ(1–42) fragment is the
dominant Aβ species in the amyloid plaques of AD patients, and, compared with Aβ(1–
40), its neurological toxicity is stronger (17-20).
Preventing the deposition of Aβ plaques in brain tissue remains a viable approach to
treating this disorder. Thus, it is of fundamental importance to elucidate the
mechanism of nucleation of Aβ fibrils and increase the ability of small ligands to
interfere with plaque formation at the molecular level. However, recently, increasing
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
Page 3 of 19 Molecular BioSystems
Mo
lecu
lar
Bio
Sys
tem
s A
ccep
ted
Man
usc
rip
t
Dow
nloa
ded
by M
ount
Alli
son
Uni
vers
ity o
n 03
/05/
2013
11:
21:4
1.
Publ
ishe
d on
23
Apr
il 20
13 o
n ht
tp://
pubs
.rsc
.org
| do
i:10.
1039
/C3M
B70
086K
View Article Online
evidence suggests that small soluble oligomers other than fibrils may be the primary
cause of synaptic dysfunction in AD (21,22). Thus, much effort has been invested in
studying the composition as well as the stability of the possible oligomeric states and
in the understanding of the amyloid fibril growth process and its modulation by
putative drug molecules through either experimental or computational studies (23-43).
Important explanations have derived from the solid state NMR structure of the
pentameric Aβ(17-42) oligomer as determined by Lhurs et al (16), which evince the
key roles behind the quaternary structure of the amyloid oligomers, although little is
still known about the structural features of the growing fibrils.
In this context, we modeled a decameric structure of the Aβ fibril composed of 10
peptides layers representing a successive step in the process of fibril growth compared
to the experimentally available structure of the Aβ(17-42) penta-oligomer (16). Here
we report the results of the pentameric as well as decameric system, subjected to a
full atomistic molecular dynamics simulation of 100 ns, alone and in complex with an
aggregate inhibitor proposed by De Bona et al. (44) to investigate if our models could
represent target systems for the design of new inhibitor compounds. Our work could
be useful for the development of new drugs able to act not only on the first phases of
the Aβ-correlated diseases, but also able to target the successive state of grown fibril.
METHODS
The solid state NMR structure obtained by Lührs et al (pdb entry 2BEG) has been used
for the pentameric system. The model of the decameric oligomer was built using two
consecutive modules of the pentameric structure. The two oligomeric systems were
subjected to 100ns molecular dynamics (MD) simulations with explicit waters.
MD simulations were performed and analysed using the GROMACS package program
and GROMOS43a1 force field (45-47). The Aβ oligomer and each Aβ-inhibitor complex
were solvated in a cubic box of simple point charge (SPC) water with at least 11 Å
distance to the border adding counter-ions to neutralize the system. Periodic boundary
conditions were employed and the LINCS algorithm (48) was used to constrain all bond
lengths. Simulations were run under NPT conditions, using Berendsen's coupling
algorithm (49) to retain the temperature and pressure constant (300 K with a P = 1
bar). Electrostatic interactions were handled by means of the particle mesh Ewald
method (PME) (50) whereas for the Lennard-Jones potential a non-bonded cutoff of 1.4
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
Page 4 of 19Molecular BioSystems
Mo
lecu
lar
Bio
Sys
tem
s A
ccep
ted
Man
usc
rip
t
Dow
nloa
ded
by M
ount
Alli
son
Uni
vers
ity o
n 03
/05/
2013
11:
21:4
1.
Publ
ishe
d on
23
Apr
il 20
13 o
n ht
tp://
pubs
.rsc
.org
| do
i:10.
1039
/C3M
B70
086K
View Article Online
nm was used. The systems were then gradually heated up to 300 K in a five-step
process, from 50K to 300 K, and then simulated under standard NPT conditions for 100
ns without restraints.
The CT ligand was built using the Insight BUILDER module (51) with the covalent
valence force field (CVFF), then minimized using a conjugate gradient algorithm (52).
The conformational space was sampled through simulated annealing with distance-
dependent dielectrics of 80, in order to mimic the aqueous environment. The ligand
energy was first minimized using the steepest-descent algorithm until the maximum
derivative was less than 0.001 kcal/A. Structures were thus heated in 1000 fs up to
1000 K for 1000 fs and 200 structures were chosen on the basis of energy and
geometric criteria. Then the CT best conformation was docked on both the systems
under investigation following the complementarity between the peptidic LPFFD
scaffold and the homologous sequence of Aβ(17-21) and leaving the sugar without
structural restraints. Subsequently, the CT ligand was parametrized using the PRODRG
server. Partial charges for the sugar portion were adjusted taking into account those
reported in reference (53). Finally, the obtained complexes were subjected to MD
simulations following the above described protocol, and subsequently both MD
ensembles have been clustered, following the algorithm described by Dura et al. (54)
with a cut-off of 0.2 nm. For the simulated complexes, the representative models of
the most populated clusters have been subjected to the MMPBSA (molecular
mechanics (MM) Poisson Boltzmann/Generalized Born surface area MM-GBSA) analysis
to obtain the relative binding free energy, using MMPBSA.py program of the AmberTool
suite. (55-56).
RESULTS AND DISCUSSION
Starting from the solid-state NMR structure of Aβ(17-42) as determined by Luhrs at al.
(16), which is formed by 5 peptide layers (or chains) in parallel in register β-sheets
packing and which is considered the repeat unit of the mature fibril, we created a
model of decameric oligomer of Aβ(17-42) amyloid (Fig 1). Following Luhr's
nomenclature, B1 refers to the residues from 18 to 26, B2 to those from 31 to 42, and
U-turn those from 27 to 30, where the turn region connects the two β-strands B1 and
B2 (Fig 1).
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
Page 5 of 19 Molecular BioSystems
Mo
lecu
lar
Bio
Sys
tem
s A
ccep
ted
Man
usc
rip
t
Dow
nloa
ded
by M
ount
Alli
son
Uni
vers
ity o
n 03
/05/
2013
11:
21:4
1.
Publ
ishe
d on
23
Apr
il 20
13 o
n ht
tp://
pubs
.rsc
.org
| do
i:10.
1039
/C3M
B70
086K
View Article Online
The analysis of the root mean square (RMS) deviation (RMSD) computed for the
pentameric as well as for the decameric system during the whole simulations reveals
that the pentameric system is less stable with respect to the decameric one (Fig 2). In
greater detail, the extreme chains are less stable in terms of secondary structure,
RMSD and RMS fluctuation (RMSF) for both systems studied as shown in Figure 3. The
closeness of each layer was monitored by measuring the distance between the center
of mass (COM) of the F19 phenyl ring of each layer and the successive one. During the
simulation, the penta-oligomer displays a comparable trend to those shown by the
deca-oligomer. Indeed, in both systems only the last chain, chain 5 and chain 10
respectively, shows an increase of the distance from the previous chain, suggesting a
detachment from the rest of the system (Fig 4).
The pentameric NMR structure (16) is principally taken into its characteristic
arrangement by important hydrophobic interactions between B1 and B2, (mainly
between the residues F19 and G38, as well as A21 and V36) and by a fundamental salt
bridge involving D23 and K28 which connects the two strands B1 and B2. These key
interactions are kept into the same layer, as intra-layer contacts, and between
adjacent layers as inter-layer contacts. In order to evaluate the structural and dynamic
features of the studied systems, these key interactions have been analyzed over the
entire course of the simulations. The stabilizing hydrophobic intra-layer contacts F19-
G38 and A21-V36 are substantially retained in both simulated systems, apart from a
sporadic loss of the A21-V36 interaction in the decamer (Supplementary Figure
S1). In contrast, the inter-layer contacts are continually lost between the last two
chains in the pentamer in a more significant way then in the decamer (Fig 5). This
fact is mirrored by the decrease of the hydrophobic contact surface between the last
two chains of the pentameric system during the simulation, as shown Figure 6. The
inter-layer D23-K28 salt bridge is often lost between the last two chains in a
comparable way for the pentameric and the decameric systems (Supplementary
Figure S2). Thus, our simulations reveal that the smaller oligomer maintains the
intra-layer interactions but shows higher vulnerability for the inter-layer contacts.
These data suggest a modest detachment of each layer from the previous one, most
pronounced in the last chain, which is attenuated with oligomer growth as suggested
by the simulation of the decameric system. It is worth noting that the higher structural
stability of the decameric system with respect to the pentamer as suggested by the
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
Page 6 of 19Molecular BioSystems
Mo
lecu
lar
Bio
Sys
tem
s A
ccep
ted
Man
usc
rip
t
Dow
nloa
ded
by M
ount
Alli
son
Uni
vers
ity o
n 03
/05/
2013
11:
21:4
1.
Publ
ishe
d on
23
Apr
il 20
13 o
n ht
tp://
pubs
.rsc
.org
| do
i:10.
1039
/C3M
B70
086K
View Article Online
analysis of our results is in agreement with previous structural studies focusing on the
stability differences between Aβ(1-42) oligomers of different size (23, 34-37, 57).
It should be noted, that our model of deca-oligomer gives represents a step forward in
the process of fibril growth, along the direction of the fibril axis (longitudinal direction
of fibril growth). In recent years, by means of MD simulations, other models of
amyloid protofibrils have been studied, which mainly focused on the fibril growth by
lateral association (direction perpendicular to the fibril axis), evaluating different
possible association interfaces (58, 35-37). Even though a direct comparison with our
results is difficult due to the significantly different systems topologies, there is a
substantial agreement in predicting the higher instability of fibril ends and a key role
of the D23-K28 salt bridge as well as of hydrophobic interactions in stabilizing the
overall structure. Though our MD analysis provides insights into the structural stability
and conformational dynamics of Aβ (17-42) oligomeric models, we cannot make any
hypothesis on the mechanism of amyloid fibril elongation and further investigations
would be needed as reported in previous studies (40, 60-61).
Finally we evaluated the effect of an aggregate inhibitor on the two oligomeric
systems. To this end, an inhibitor, CT, proposed by De Bona et al. (44) has been
chosen, this is a trehalose derivative of the well known Soto’s peptide (LPFFD) (62).
This peptide retains high affinity toward the self-complementary LVFFA region of Aβ(17-
21) and was thought to block elongation by binding to the ends of the growing fibril.
So, the initial complex of inhibitor-amyloid oligomer was built following the original
idea of Soto and coworkers stated above by manually docking the peptide portion of
CT ligand to the 17-21 region of the first chain of the Aβ oligomers in a self
complementary manner, leaving the sugar portion free (Fig 7). The resulting
complexes of CT inhibitor with both the penta- and deca-oligomers were subjected to
100ns MD simulations without any positional restraint in order to avoid bias, allowing
the ligand to freely move and eventually go away from the system. Indeed, the ligand
significantly moves during the simulations, but it never detaches from the oligomers
neither from the starting binding site (Supplementary Figure S3).
The secondary structure of both systems is not affected in a very significant way by
the presence of the inhibitor, however, there is a not negligible increase of the
residues in coil structure simultaneously with a decrease of the residues in β-sheet and
bend structure, most pronounced for the penta-oligomeric system (Table 1). Moreover,
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
Page 7 of 19 Molecular BioSystems
Mo
lecu
lar
Bio
Sys
tem
s A
ccep
ted
Man
usc
rip
t
Dow
nloa
ded
by M
ount
Alli
son
Uni
vers
ity o
n 03
/05/
2013
11:
21:4
1.
Publ
ishe
d on
23
Apr
il 20
13 o
n ht
tp://
pubs
.rsc
.org
| do
i:10.
1039
/C3M
B70
086K
View Article Online
the stability of the global structure is differently influenced by CT ligand binding: i) the
RMSD value over the evolution of the simulations increases due to the presence of the
ligand in both systems. Nonetheless, for the decameric system the RMSD gap is
constant and of only about 0.1 nm, while for the pentamer this growth becomes
significant at the end of the simulation, reaching a value of 0.5 nm, thus suggesting a
most important destabilization in the latter case (Fig 8). ii) The RMSF value in the
decameric system grows only for the chain interacting with the ligand, whereas all the
chains composing the pentameric oligomer show an increase in fluctuations when the
system interacts with the ligand, in a relevant and more remarkable manner for the
first chain (Fig 9). Therefore, our results underline a most pronounced destabilizing
effect due to the interaction with the inhibitor in the case of the pentameric oligomer.
In order to gain insights into the energetics of the two complexes, we have performed
binding free energy calculations using the MMPBSA method (55). The computed
energies for ligand binding to the penta- and the deca-oligomers are -16.6 kcal/mol
and -13.6 kcal/mol, respectively. These data clearly indicate that the ligand binding to
the fibril edge is energetically favorite. Although it cannot be excluded that other CT
binding sites along Aβ fibril exist, it should be noted that the binding to the amyloid
fibril edge is emerged as an important player into inhibition, also in recent MD studies
of two nonsteroidal anti-inflammatory drugs as well as of Congo red (41-43). In these
regards, a more comprehensive analysis of the inhibitor-oligomer interactions at the
atomic level is needed.
Table 1: Percentage of the residue showing the Coil, B-Sheet and Bend Secondary Structures during the last 50 ns of Molecular Dynamics simulations.
Secondary Structure
Pentamer
Pentamer-CT
Decamer
Decamer-CT
Coil (%) 26.7 28.1 20.0 21.0
B-Sheet (%) 63.1 62.9 71.7 70.6
Bend (%) 8.4 7.8 7.3 7.2
CONCLUSIONS
We have built a model of the deca-oligomer Aβ(17-42) which represents a step forward
in the process of fibril growth compared to the experimentally available structure of
the penta-oligomer (16). The structure and dynamics of both the penta- and the deca-
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
Page 8 of 19Molecular BioSystems
Mo
lecu
lar
Bio
Sys
tem
s A
ccep
ted
Man
usc
rip
t
Dow
nloa
ded
by M
ount
Alli
son
Uni
vers
ity o
n 03
/05/
2013
11:
21:4
1.
Publ
ishe
d on
23
Apr
il 20
13 o
n ht
tp://
pubs
.rsc
.org
| do
i:10.
1039
/C3M
B70
086K
View Article Online
oligomers of Aβ alone and in complex with an agglomerate inhibitor (44) have been
studied with MD simulations. Our simulations reveal that in the deca-oligomer all the
structural elements are well retained; moreover it is more stable compared to the
pentameric system in agreement with previous studies (23, 34, 35). Interestingly, MD
studies of the interaction between the inhibitor and the two systems indicate that the
inhibitor acts not only by hampering the addition of successive layers at the ends of
the oligomers but also by affecting the structure and stability of the oligomers. It is
worthy of note that the strategy we employed in our study could be fruitfully applied
to build models of the elongated Aβ fibril in order to study the effect of putative
inhibitors at different stages of Aβ−related disease.
ACKNOWLEDGEMENTS
The authors would like to give their special thanks to Dr. Giuseppe Pappalardo for his
help, support, interest and valuable hints. The technical assistance of Mr. Luca de Luca
and Mr. Giovanni Filograsso are gratefully acknowledged. This work was supported by
the M.I.U.R, Ministero dell'Istruzione, dell'Università e della Ricerca (FIRB-MERIT
RBNE08HWLZ_002).
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
Page 9 of 19 Molecular BioSystems
Mo
lecu
lar
Bio
Sys
tem
s A
ccep
ted
Man
usc
rip
t
Dow
nloa
ded
by M
ount
Alli
son
Uni
vers
ity o
n 03
/05/
2013
11:
21:4
1.
Publ
ishe
d on
23
Apr
il 20
13 o
n ht
tp://
pubs
.rsc
.org
| do
i:10.
1039
/C3M
B70
086K
View Article Online
REFERENCES
1 C. Adessi and C. Soto, Drug Dev Res., 2002, 56, 184-193
2 J. Hardy and J.S. Dennis, Science, 2002, 297, 353-356
3 J. Clippingdale, Pept.Sci., 2001, 7, 227-249
4 W. Annaert and B. De Strooper, Annu. Rev. Cell Dev. Biol., 2002, 18, 25-51
5 M. Hutton, J. Perez-Tur and J. Hardy, Essays Biochem. 1998, 33, 117-131
6 D. J. Selkoe, Nat. Cell Biol., 2004, 6, 1054-1061
7 J. P. Taylor, J. Hardy and K. H. Fischbeck, Science, 2002, 296, 1991-1995
8 W.L. Klein, W.B. Stine and D.B. Teplow. Neurobiol Aging., 2004, 25, 569-580
9 J.T. Jarrett, E.P. Berger and P.T. Lansbury, Biochemistry, 1993, 32, 4693-4697
10 J.T. Jarrett, Berger EP and Lansbury PT, Ann NY Acad Sci., 1993, 695, 144-148
11 A. T. Petkova, W. M. Yau, R. Tycko, Biochemistry, 2006, 45, 498-512
12 E. D. Eanes and G. G. J Glenner, Histochem. Cytochem., 1968, 16, 673-677
13 D. A. Kirschner, C. Abraham and D. J. Selkoe, Proc. Natl. Acad. Sci., U.S.A. 1986, 83,
503-507
14 A. T. Petkova, Y. Ishii, J. Balbach, O. Antzutkin, R. Leapman, F. Delaglio and R. Tycko,
Proc.Natl. Acad. Sci.U.S.A., 2002, 99, 16742-16747
15 J.J. Balbach, A. T .Petkova, N. A. Oyler, O. N. Antzutkin, D. J, Gordon, S.C. Meredith
and R. Tycko, Biophys. J., 2002, 83, 1205-1216
16 T. Luhrs, C. Ritter, M. Adrian, D. Riek-Loher, B. Bohrmann, H. Doeli, D.Schubert and
R. Riek, Proc. Natl. Acad. Sci. U.S.A, 2005, 102, 17342-17347
17 J. T. Jarrett and P.T.J. Lansbury, Cell, 1993, 73, 1055-1058.
18 L. Mucke, E. Masliah, G. Q. Yu, M. Mallory, E. M Rockenstein, G. Tatsuno,K. Hu, D.
Kholodenko and K. Johnson-Wood, L. McConlogue, J Neurosci., 2000, 20, 4050-4058.
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
Page 10 of 19Molecular BioSystems
Mo
lecu
lar
Bio
Sys
tem
s A
ccep
ted
Man
usc
rip
t
Dow
nloa
ded
by M
ount
Alli
son
Uni
vers
ity o
n 03
/05/
2013
11:
21:4
1.
Publ
ishe
d on
23
Apr
il 20
13 o
n ht
tp://
pubs
.rsc
.org
| do
i:10.
1039
/C3M
B70
086K
View Article Online
19 A. Schmechel, H. Zentgraf, S. Scheuermann, G. Fritz, R.D. Pipkorn, J. Reed, K.
Beyreuther, T.A. Bayer, and G. Multhaup, J Biol Chem, 2003, 278, 35317-35324.
20 Y. Zhang, R. McLaughlin, C. Goodyer and A. LeBlanc, J Cell Biol, 2002, 156- 519
21 C.A. McLean, R.A. Cherny, F.W. Fraser, S.J. Fuller, M.J. Smith, K. Beyreuther, A.I. Bush
and C.L. Masters, Ann Neurol, 1999, 46, 860-866
22 C. Haassand and D.J. Selkoe DJ, Nat Rev Mol Cell Biol, 2007, 8, 101-112
23 A. H.C. Horn and H. Sticht, J.Phys.Chem., B 2010, 1141, 2219-222
24 N.J. Bruce, D. Chen, S.G. Dastidar, G.E. Marks, C.H. Schein and R. A. Bryce,
Peptides, 2010, 31, 2100-2108
25 C. Soto, G.P. Saborio and B. Permanne, Acta Neurol Scandinavica, 2000, 176, 90-95
26 A.J. Lemkul and D.R. Bevan, J. Phys. Chem., 2010, 114, 1652-1660
27 B. Ma and R. Nussinov, PNAS, 2002, 99, 14126-14131
28 L.O. Tjernberg, D.J.E. Callaway, A. Tjernberg, S. Hahne, C. Lilliehook, L. Terenius, J.
Thyberg and C Nordstedt, J. Biol. Chem, 1999, 274, 12619-12625,
29 K. Hochdörffer, J. März-Berberich, L. Nagel-Steger, M. Epple, W. Meyer-Zaika, A.H.
Horn, H. Sticht, S. Sinha, G. Bitan and T. Schrader, J Am Chem Soc, 2011, 133, 4348-
4358
30 P. Rzepecki and T. Schrader, J Am Chem Soc, 2005, 127, 3016-3025
31 B. Urbanc, L. Cruz, F. Ding, D. Sammond, S. Khare, S.V. Buldyrev, H.E. Stanley and
N.V. Dokholyan, Biophys J., 2004, 87, 2310-2321
32 B. Urbanc, L. Cruz, S. Yun, S.V. Buldyrev, G. Bitan, D.B. Teplow and H.E. Stanley,
Proc Natl Acad Sci USA, 2004, 101, 17345-17350.
33 S. Yun, B. Urbanc, L. Cruz, G. Bitan, D.B. Teplow and H.E. Stanley, Biophys J., 2007,
92, 4064-4077
34 F.M. Masman, U.l.L.M. Eisel, I. G. Csizmadia, B. Penke, R. D. Enriz, S. J. Marrink and P.
G.M. Luiten, J. Phys. Chem. 2009, 113, 11710-11719.
35 J. Zheng, H. Jang, B. Ma, C.J. Tsai, R. and Nussinov, Biophys J., 2007, 93, 3046-57.
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
Page 11 of 19 Molecular BioSystems
Mo
lecu
lar
Bio
Sys
tem
s A
ccep
ted
Man
usc
rip
t
Dow
nloa
ded
by M
ount
Alli
son
Uni
vers
ity o
n 03
/05/
2013
11:
21:4
1.
Publ
ishe
d on
23
Apr
il 20
13 o
n ht
tp://
pubs
.rsc
.org
| do
i:10.
1039
/C3M
B70
086K
View Article Online
36 Y. Miller, B. Ma and R. Nussinov, Biophys J., 2009, 97, 1168-77.
37 J. Zheng, H. Jang, B. Ma and R. Nussinov, J Phys Chem B., 2008, 112, 6856-65.
38 R. Friedman, R. Pellarin and A. Caflisch, J Mol Biol., 2009, 407-415.
39 F. Tofoleanu and N.V. Buchete, J Mol Biol., 2012, 572-586.
40 T. Takeda and D. K. Klimov, Biophys. J. 2008, 95,1758-1772.
41 E. Raman, P.T. Takeda and D. K. Klimov, Biophys. J., 2009, 97, 2070-2079.
42 T. Takeda, W.L.E. Chang and D.K. Klimov, Proteins, 2010, 78, 2849-2860.
43 C. Wu, J. Scott and J-E Shea, Biophys J., 2012,103, 550-557.
44 P. De Bona, M.L. Giuffrida, F. Caraci, A. Copani, B. Pignataro, F. Attanasio, S. Cataldo,
G. Pappalardo and E. Rizzarelli, J. Pept Scie., 2009 , 15, 220-228
45 H.J.C. Berendsen, H. J. C., van der Spoel and D. van Drunen, R. Comp. Phys. Comm.,
1995, 91, 43–56,.
46 D. van der Spoel, E. Lindahl, B. Hess, G. Groenhof, A.E. Mark and H.J.C. Berendsen,
J. Comp. Chem., 2005 26, 1701–1718.
47 B. Hess, C. Kutzner, D. van der Spoel and E. Lindahl, J. Chem. Theory Comp, 2008 4,
435–447
48 B. Hess, H. Bekker, H.J.C. Berendsen and J.G.E.M. Fraaije, J Comp Chem., 1997, 18,
1463–1472
49 H. J. C. Berendsen, J. P. M. Postma, W. F. Van Gunsteren, A. DiNola, and J. R. Haak,
J.Chem. Phy, 1984, 81,3648-3690
50 T. Darden, D. York and L. Pedersen, J. ChemPhys, 1993, 98, 10089-10092
51 Insight2000, Accelrys, San Diego
52 R. Fletcher, John Wiley & Sons, 1980, Vol. 1
53 E. Fadda and R.J. Woods, Drug Discovery today, 2010, 15, 596-609
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
Page 12 of 19Molecular BioSystems
Mo
lecu
lar
Bio
Sys
tem
s A
ccep
ted
Man
usc
rip
t
Dow
nloa
ded
by M
ount
Alli
son
Uni
vers
ity o
n 03
/05/
2013
11:
21:4
1.
Publ
ishe
d on
23
Apr
il 20
13 o
n ht
tp://
pubs
.rsc
.org
| do
i:10.
1039
/C3M
B70
086K
View Article Online
54 X. Dura, K. Gademann, J. Bernhard, D. Seebach, F. Wilfred van Gunsteren and E.
Mark, Angev.Chem.Int.Ed., 1999, 38, 236-240.
55 B.R. Miller, T.D. McGee, Jr. J.M. Swails, Homeyer,H. Gohlke and A.E. Roitberg, J.
Chem. Theory Comput., 2012, 8, 3314-3321.
56 D.A. Case, T.A. Darden, T.E. Cheatham, III, C.L. Simmerling, J. Wang, R.E. Duke, R.
Luo, R.C. Walker, W. Zhang, K.M. Merz, B. Roberts, S. Hayik, A. Roitberg, G. Seabra, J.
Swails, A.W. Goetz, I. Kolossváry, K.F. Wong, F. Paesani, J. Vanicek, R.M. Wolf, J. Liu, X.
Wu, S.R. Brozell, T. Steinbrecher, H. Gohlke, Q. Cai, X. Ye, J. Wang, M.-J. Hsieh, G. Cui,
D.R. Roe, D.H. Mathews, M.G. Seetin, R. Salomon-Ferrer, C. Sagui, V. Babin, T. Luchko,
S. Gusarov, A. Kovalenko, and P.A. Kollman (2012), AMBER 12, University of California,
San Francisco
57 U.F. Rohrig, A. Laio, N. Tantalo, M. Parrinello and R. Petronzio, Biophys J. 2006, 91,
3217-3229
58 N.V. Buchete and G. Hummer, Biophys. J., 2007, 92, 3032-3039
59 N.V. Buchete and G. Hummer, J. Mol. Biol., 2005, 353, 804-821
60 M. Schor, J. Vreede and P. G. Bolhuis, Biophys. J., 2012, 103, 1296-1304
61 E.P. O’Brien, Y. Okamoto, J.E. Straub, B.R. Brooks and D. Thirumala, J. Phys. Chem. B
2009, 113, 14421–14430.
62 C. Soto, E. M. Sigurdsson, L. Morelli, R.A. Kumar, E.M. Castaño and B. Frangione,
Nat. Med. 1998, 4, 822-826.
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
Page 13 of 19 Molecular BioSystems
Mo
lecu
lar
Bio
Sys
tem
s A
ccep
ted
Man
usc
rip
t
Dow
nloa
ded
by M
ount
Alli
son
Uni
vers
ity o
n 03
/05/
2013
11:
21:4
1.
Publ
ishe
d on
23
Apr
il 20
13 o
n ht
tp://
pubs
.rsc
.org
| do
i:10.
1039
/C3M
B70
086K
View Article Online
FIGURE CAPTIONS
Figure 1: Structural details of the pentameric and decameric oligomer systems. A)
Top-view of the Aβ(17-42) pentamer (16). Upper panel Labels used to indicate
secondary structure elements (B1, B2) and connecting regions (U-turn) are shown.
Residues belonging to the different secondary structure regions are displayed in sticks
of different color. Lower panel Pairs of residues involved in fibril stabilizing interactions
are displayed in sticks of different color. B) Side-view of the pentameric (top) and
decameric (bottom) Aβ oligomers colored by chain.
Figure 2: Root mean square deviation (RMSD) of the Cα carbon atoms positions in the
pentameric (black line) and decameric (red line) systems as a function of time.
Figure 3: Single chain RMS deviation and RMS fluctuation for the pentameric (A,C)
and decameric (B, D) systems.
Figure 4: Interchain distance between the centers of mass (COM) of the F19 phenyl
rings as a function of time for the penta- and deca-oligomer (top and bottom,
respectively).
Figure 5: F19-G38 and A21-V36 inter-layer Cα distances as a function of time for the
pentamer (A and C) and decamer (B and D), respectively.
Figure 6: Hydrophobic interacting area among the chains composing the penta- and
deca-oligomeric systems (top and bottom, respectively).
Figure 7: Top view of the binding mode of CT with the Aβ first layer. Oligomer residues
are depicted in blue, CT peptidic and sugar portions are in red and yellow, respectively.
Figure 8: Root mean-square deviation of the Cα carbon atom positions in the
pentameric (black line) and decameric (red line) complexes as a function of time.
Figure 9: Single chain RMSF for the pentameric (top) and decameric (bottom)
complexes.
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
Page 14 of 19Molecular BioSystems
Mo
lecu
lar
Bio
Sys
tem
s A
ccep
ted
Man
usc
rip
t
Dow
nloa
ded
by M
ount
Alli
son
Uni
vers
ity o
n 03
/05/
2013
11:
21:4
1.
Publ
ishe
d on
23
Apr
il 20
13 o
n ht
tp://
pubs
.rsc
.org
| do
i:10.
1039
/C3M
B70
086K
View Article Online
Figure 1.
Figure 2.
343
344
345
346
347
348
Page 15 of 19 Molecular BioSystems
Mo
lecu
lar
Bio
Sys
tem
s A
ccep
ted
Man
usc
rip
t
Dow
nloa
ded
by M
ount
Alli
son
Uni
vers
ity o
n 03
/05/
2013
11:
21:4
1.
Publ
ishe
d on
23
Apr
il 20
13 o
n ht
tp://
pubs
.rsc
.org
| do
i:10.
1039
/C3M
B70
086K
View Article Online
Figure 3.
Figure 4.
350
351
352
353
354
355
356
357
358
Page 16 of 19Molecular BioSystems
Mo
lecu
lar
Bio
Sys
tem
s A
ccep
ted
Man
usc
rip
t
Dow
nloa
ded
by M
ount
Alli
son
Uni
vers
ity o
n 03
/05/
2013
11:
21:4
1.
Publ
ishe
d on
23
Apr
il 20
13 o
n ht
tp://
pubs
.rsc
.org
| do
i:10.
1039
/C3M
B70
086K
View Article Online
Figure 5.
Figure 6.
360
361
362
363
364
365
366
367
368
Page 17 of 19 Molecular BioSystems
Mo
lecu
lar
Bio
Sys
tem
s A
ccep
ted
Man
usc
rip
t
Dow
nloa
ded
by M
ount
Alli
son
Uni
vers
ity o
n 03
/05/
2013
11:
21:4
1.
Publ
ishe
d on
23
Apr
il 20
13 o
n ht
tp://
pubs
.rsc
.org
| do
i:10.
1039
/C3M
B70
086K
View Article Online
Figure 7.
Figure 8.
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
Page 18 of 19Molecular BioSystems
Mo
lecu
lar
Bio
Sys
tem
s A
ccep
ted
Man
usc
rip
t
Dow
nloa
ded
by M
ount
Alli
son
Uni
vers
ity o
n 03
/05/
2013
11:
21:4
1.
Publ
ishe
d on
23
Apr
il 20
13 o
n ht
tp://
pubs
.rsc
.org
| do
i:10.
1039
/C3M
B70
086K
View Article Online
Figure 9.
384
385
386
Page 19 of 19 Molecular BioSystems
Mo
lecu
lar
Bio
Sys
tem
s A
ccep
ted
Man
usc
rip
t
Dow
nloa
ded
by M
ount
Alli
son
Uni
vers
ity o
n 03
/05/
2013
11:
21:4
1.
Publ
ishe
d on
23
Apr
il 20
13 o
n ht
tp://
pubs
.rsc
.org
| do
i:10.
1039
/C3M
B70
086K
View Article Online