# Effect of Cholesterol Molecules on Aβ1-42 wild-type and mutants trimers

Trung Hai Nguyen, Phuong H. Nguyen, Son Tung Ngo, Philippe Derreumaux

► To cite this version:

Trung Hai Nguyen, Phuong H. Nguyen, Son Tung Ngo, Philippe Derreumaux. Effect of Cholesterol Molecules on Aβ1-42 wild-type and mutants trimers. Molecules, 2022, Special Issue Frontiers in Protein Folding and Related Areas – in Memory of Professor Sir Christopher M. Dobson (1949–2019), 27 (4), pp.1395. 10.3390/molecules27041395. [https://doi.org/10.3390/molecules27041395](https://doi.org/10.3390/molecules27041395)

HAL Id: hal-03853301

https://hal.science/hal-03853301

Submitted on 15 Nov 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Effect of Cholesterol Molecules on Aβ1-42 wild-type and mutants trimers

Trung Hai Nguyen,<sup>1,2</sup> Phuong H Nguyen,<sup>3,4</sup> Son Tung Ngo,<sup>1,2</sup> and Philippe Derreumaux <sup>3,4,5</sup>

1 Laboratory of Theoretical and Computational Biophysics, Ton Duc Thang University, Ho Chi Minh City, Vietnam.

2 Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam.

3 CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, 13 rue Pierre et Marie Curie, 75005 Paris, France.

4 Institut de Biologie Physico-Chimique, Fondation Edmond de Rothschild, PSL Research University, Paris, France.

5 Institut Universitaire de France (IUF).

ABSTRACT

Alzheimer's disease displays aggregates of the amyloid-beta (Aβ) peptide in the brain and there is increasing evidence that cholesterol may contribute to the pathogenesis of the disease. Though many experimental and theoretical studies have focused on the interactions of Aβ oligomers with membrane models containing cholesterol, an understanding of the effect of free cholesterol on small Aβ42 oligomers is not fully established. To address this question, we report on replica exchange with solute tempering simulation of Aβ42 trimer with cholesterol and compare with a previous replica exchange molecular dynamics simulation. We show that the binding hot spots of cholesterol are rather complex involving essentially the hydrophobic residues L17-F20 and L30-M35 with a non-negligible contribution of the loop residues D22-K28 and the N-terminus residues. We also examine the effects of cholesterol on the trimers of the disease-causing A21G and disease-protective A2T mutations by molecular dynamics simulations. We show that these two mutations impact moderately the cholesterol binding modes. In our REST2 simulations, we find that cholesterol is rarely inserted into the aggregates, but rather attached as dimers and trimers at the surface of Aβ42 oligomers. We propose that cholesterol acts as a glue to speed up the formation of larger aggregates, this provides a mechanistic link between cholesterol and Alzheimer's disease.
Introduction

The two hallmarks of Alzheimer's disease (AD) are extracellular amyloid-beta (Aβ) plaques of Aβ42 and Aβ40 peptides, the 42-residue species being the most toxic, and intracellular neurofibrillary tangles built from hyperphosphorylated tau protein [1]. Despite extensive research, all drugs targeting Aβ and tau oligomers have failed in AD [2,3]. Among many cellular factors contributing to AD development, cholesterol plays a critical role by different actions.

First, cholesterol is present in micro-dissected AD senile plaques with a molar ratio of 1:1 [4] and purified AD paired helical fragments of tangles [5]. The level of plasma cholesterol is 10% higher in AD patients than in normal individuals [6]. Cholesterol levels in the brain positively correlate with the severity of dementia in AD patients [7] in contrast to Aβ plaque burden which correlates weakly with disease severity [8].

Second, cholesterol impacts Aβ production through the cleavage of the amyloid protein precursor (APP) [9]. Based on nuclear magnetic resonance and electron paramagnetic resonance, the residues E22-N27 in the bulk solution and the membrane-buried residues G29-G33 (using Aβ42 amino acid numbering) play a key role in cholesterol binding to APP [10]. Recently, it has been proposed from atomistic molecular dynamics (MD) simulations that cholesterol modulates the conformation and activity of the C-terminal domain of APP directly through hydrogen bonding and indirectly through induction of the liquid ordered phase [11]. Furthermore, Aβ accumulation in neurons is tightly regulated by cholesterol production in astrocytes [12] and elevated membrane cholesterol alters lysosomal degradation to induce Aβ degradation [13]. Finally, the apolipoprotein E gene, a major transporter of cholesterol in the brain, and in particular its e4 allele, the most genetic risk factor for AD, is associated with higher total cholesterol [14].

Third, numerous *in vitro* experiments have explored Aβ aggregation in a water-phospholipid-cholesterol membrane environment. They showed that cholesterol favors formation of Aβ pores in the membrane of brain cells [15]. Cholesterol inserted into membranes facilitates aggregation of Aβ at the surface of the membranes at physiological concentrations, as reported by atomic force microscopy images [16]. Lipid membranes containing cholesterol also enhance the primary nucleation rate of Aβ42 aggregation by up to 20-fold, as reported by aggregation kinetics experiments [17].
At the theoretical level, many all-atom MD simulations have been performed to study the interactions between Aβ42 peptides and membranes with various contents of cholesterol, revealing the mechanisms by which cholesterol changes the mechanisms of Aβ monomers and dimers binding and affinity to the lipid bilayer [16,18-20]. Binding of Aβ42 amyloid fibrils to lipid bilayer has also been investigated by using coarse-grained MD simulations and showed that the addition of cholesterol increases the binding frequency and alters the binding interface and contacts [21].

Little is known on the interactions of Aβ42 peptides with free cholesterol in the bulk solution. The effect of unmodified cholesterol and charged cholesterol derivatives on Aβ40 fibril formation was explored by kinetic aggregation assay Thioflavin T, atomic force microscopy and dynamic light scattering. At a concentration lower than the critical micelle concentration, and therefore assuming that cholesterol exists mainly as monomeric molecules, unmodified cholesterol and positively and negatively charged cholesterol derivatives accelerate the aggregation rate of Aβ40, the aggregation half time being reduced by 15 to 25% compared to Aβ40 peptide in phosphate buffer [22].

Using transmission electron microscopy, Harris *et al*. demonstrated the binding of soluble 10 nm diameter cholesterol-PEG 600 micelles to Aβ42 fibrils and proposed that binding involves the central hydrophobic core CHC spanning residues L17-A21 [23]. Atomistic 20 ns MD simulation of a system of 4 Aβ42 - 12 cholesterol molecules in a box of water showed a stable contact between the benzyl group of F19 and the cholesterol micelle, which forms a flat surface and, notably the steroid group of cholesterol [24]. It was suggested that the carboxyl terminus spanning residues have a higher tendency to interact with cholesterol and form β-sheet conformation [24].

It is well-established that lipids can be removed from the membrane because of oligomers - membrane interactions [25-27]. Recent experimental studies indicate the critical role of free phospholipids at a nanomolar to micromolar concentration in equilibrium with the membrane (large unilamellar vesicles) in forming Aβ-lipids complexes that favor amyloid-membrane poration [28,29]. We recently performed atomistic replica exchange molecular dynamics (REMD) simulations of Aβ42 dimer and trimer with cholesterol using the protein AMBER99sb-ildn force field (also called AMBER ff99sb-ildn) with TIP3P water force field. We reported on the drastic
effect of cholesterol on the conformational ensemble of both Aβ42 species and showed multiple transient binding modes involving essentially the residues L17-A21 and L30-M35 [30].

In this study, we have further studied the impact of cholesterol with a ratio 1:1 on Aβ42 trimer in the bulk solution. First, replica exchange with solute tempering (REST2) was chosen for sampling, as this procedure converges faster to equilibrium than REMD, while requiring a smaller number of replicas [31-33]. Second, we used CHARMM36m, as this force field has been shown to be more relevant for intrinsically disordered proteins [34] and small oligomers of Aβ [35-39]. This allows comparison with our previous REMD study using AMBER99sb-ildn TIP3P [30] Third, we have examined the impact of the A2T mutation, known to be AD protective, and the disease-causing A21G mutation [3] on the interactions of Aβ trimer - cholesterol by MD simulations. It is to be noted that the A21G and A2T Aβ42 oligomers have been extensively studied by experimental and theoretical means in the bulk solution [40-46], but none study reports on their interactions with free cholesterol.

## Materials and Methods

As in our previous REMD study [30], the cholesterol molecule was geometrically optimized by quantum mechanics with the B3LYP functional and the 6-31G(d,p) basis set, and the cholesterol force field was parameterized using the general Amber force field (GAFF) [42], the Restrained Electrostatic Potential (RESP) method [47] and the B3LYP functional and the 6-31G(d,p) basis set [48].

In all simulations, Aβ42 peptides at pH 7 have NH3⁺ and CO2⁻ termini, deprotonated Glu and Asp, protonated Arg and Lys, and neutral His with a protonated N<sub>e</sub> atom. The Aβ42/cholesterol system was neutralized by sodium ions. Prior to MD and REST2 simulations, all systems were minimized with harmonic restraints on the positions of the peptide Cα atoms, and then equilibrated by MD at 310 K for 1 ns in the NVT ensemble, followed by 2 ns in the NPT ensemble.

In REST2, we used the same initial structure of the Aβ42 trimer – cholesterol system as in our REMD simulation inserted into a cubic water box of 7.26 nm size and 382.66 nm³ volume (10869 water molecules). This structure with high β-hairpin propensity of each chain has already been discussed in many simulations of Aβ monomer [3,22] and oligomers [28,37,38,49,50] and in exploring amyloid oligomers with peptide model system [51]. The equilibrated structure shown on Figure 1 started from three separated β-hairpins with one hairpin perpendicular to the other two
peptides and three cholesterol molecules randomly positioned and orientated with respect to Aβ42
peptides with a minimal distance of 1 nm from all Aβ42 atoms.

REST2 simulation with the CHARMM36m-TIP3P force field was performed using NAMD [52]. We used a time step of 2 fs, a cutoff of 1.2 nm for Van der Waals interactions and a cutoff of 1.1 nm for electrostatic interactions using the particle mesh Ewald (PME) method [53]. REST2 scales the solute interactions by λ with the solute consisting of Aβ42 and cholesterol molecules, the solute-water interactions by λ¹/² and leaves the water-water interactions unaltered. Using T_min = 310 K and T_max = 500 K with the number of replicas set to 16, we performed REST2 simulation in the NPT ensemble at the temperature of 310 K for 250 ns on 16 replicas exchanging the following solute-solute corresponding temperatures of 310, 320, 330.4, 341.1, 352.1, 363.5, 375.3, 387.5, 400, 413, 426, 440.1, 454.4, 469.1, 484.3 and 500 K, i.e., for λ varying between 1 and 0.75. Exchanges of configurations between neighboring replicas were attempted every 1 ps. In REST2 simulation, we used Nose-Hoover Langevin pressure control and Langevin temperature control as described in NAMD.

Note that MD simulations of the selected WT, A2T and A21G systems were performed at 310 K using GROMACS [54] with the velocity-rescaling thermostat [55], a cutoff of 1.2 nm for Van der Waals interactions and a cutoff of 1.1 nm for electrostatic interactions using the PME method.

To determine the interface between Aβ42 chains and cholesterol, we calculated the distances between heavy atoms of cholesterol and non-hydrogen atoms of the side chain residues of each Aβ42 chain. We considered a contact formed between a side chain residue and cholesterol if there is at least one distance below 0.45 nm. For clarity, we report on the probability of side chain contacts per Aβ42 chain and cholesterol molecule. We also calculated the percentage of cholesterol monomers, dimers+monomer and trimers. Monomers of cholesterol exist if they are no contacts between any two cholesterol, and trimers of cholesterol are formed if there is at least one contact between all cholesterol, a contact being defined if the intermolecular distance between any two heavy atoms is less than 0.45 nm. The DSSP protocol was used to determine the secondary structure of Aβ42 peptides [56].
The full energy landscape was approximated by 2D energy landscape. First, the trajectory was projected on the first two principal components obtained after diagonalization of the positional fluctuation covariance matrix of the backbone peptides and cholesterol atoms. Then the free energy landscape was constructed from the previous data using the formula -RT × log(H(x,y)), where H(x,y) is the histogram of the two selected order parameters x and y [57]. The population of each minimum was determined by counting all conformations around each minimum. The representative structure or center of each cluster was obtained by the Daura clustering method [58].

In what follows, we used for comparison the time intervals 50-250 ns of the present REST2 simulation at 310 K, and of our previous REMD simulation at 315 K. Note that if death is imminent for humans at 315 K, atomistic protein simulations at 310 and 315 K lead to very similar thermodynamic and structure properties. Also, a pure cholesterol bilayer has already been explored at 310 K by all-atom MD simulations [59], and REST2 has been successfully shown to simulate very well the weak binding of Aβ40 peptide on a lipid bilayer [32], and the lateral equilibration in mixed cholesterol-DPPC bilayers [60].

## Results and Discussion

The sampling of the REST2 simulation using 16 replicas is first illustrated by the good overlap of the distributions of the potential energy between neighboring replicas (Figure 2A). The exchanges of coordinates as a function of simulation time for replicas 1 (Figure 2B) and 16 (Figure 2C) indicate exploration of the full replica space in the nanosecond time scale. Overall, the average acceptance probability at different replicas is 0.3. The convergence of the REST2 simulation at 310 K is also assessed by the very high similarity of the two FELs using the single REST2 trajectory at the time intervals 50-220 ns and 50-250 ns (Figure 3A and 3B).

The secondary structure compositions of Aβ42 slightly vary between the two simulations, reaching β-sheet and coil contents of 37% and 36% in REMD vs. 36% and 43% in REST2. There are, however, differences along the amino acid sequence. Both methods give the same β-strand content for the residues 15-21 and residues 28-37 (Figure 4A) and the same turn probability for residues 22-27 (Figure 4B). The β-strand probability of residues 3-5 and the turn character of residues 6-10 are, however, reduced by a factor of two from REMD to REST2 (18% β-strand and
30% turn by REST2), leading to an enhance of coil character of the N-terminus (residues 1-10) in REST2 compared to REMD (Figure 4C). Also, REST2 explores more β-strands at positions 36-41 (Figure 4A).

Differences in the conformational space explored by the REST2 and REMD simulations is further analyzed by comparing the two FEL's using the combined REMD and REST2 trajectories to compute the first two principal components. Though direct comparison is not possible because we use two distinct force fields, there is some overlap, but overall, the free energy landscapes are different. In contrast to the REMD FEL (Figure 3C), the REST2 FEL is divided into well separated regions (Figure 3D). Also, the amplitudes of fluctuation along the PC1 and PC2 are also much larger in REST2, and the REST2 FEL is dominated by two states S1 and S2, representing 40% and 22% of the full ensemble. In contrast, the REMD FEL is dominated by two states representing 12% and 28% of the conformational ensemble.

Deviation between the REST2 and REMD conformational ensembles is observed for the probability of side chain contacts per cholesterol molecule and Aβ42 chain (Figure 5). Using REST2, the highest binding spot with cholesterol molecules involves the side-chains of the CHC region (residues 17-21) with a probability of 33% for the aromatic interaction with F19 and 17% for the interaction with L17, and then the side-chains of residues L30, I31 and I32 and residues L34 and M35 (average probability of 12%). The same residues are identified with REMD, but the probability decreases notably for F19 to 25% and increases moderately for the residues H6, Y10, V12, H13 and N27, most probabilities remaining < 10% however. There is almost no difference between the REMD and REST2 populations of (free monomers of cholesterol, dimers + monomer of cholesterol and trimers of cholesterol), as they reach (0.5%, 43% and 56%) in REST2 vs. (0.6%, 47% and 52%) in REMD. It is found that the network of interactions between the side-chains of Aβ42 varies with the aggregated forms of cholesterol. Using REST2, there are 15 contacts per Aβ42 chain and per cholesterol with a probability > 7.5% when cholesterol is in dimers + monomer form, while there are 8 contacts formed when cholesterol is in trimer form, and this is accompanied by a substantial probability reduction with the CHC region and residues 30-41 (Figure 6).

Overall, there are non-negligible differences between the REST2 and REMD results. REST2 results emphasize the role of the CHC region L17-A21 and the hydrophobic residues L30-M35 and V39-I41 in the binding with the dimer form of cholesterol, and the residues L17, F19,
A21, V24, N27, I31 and I32 in the binding with the trimer forms of cholesterol. This binding mechanistic view between the hydrophobic region of Aβ42 and cholesterol explains why the nonvesicle forms of the negatively charged cholesterol sulfate and the cationic cholesterol derivative, 3β(N-dimethylaminoethane)carbaloyl)-cholesterol, moderately change the aggregation kinetics of Aβ40 [22]. Our cholesterol binding mechanism with an average probability of contacts of 16% averaged over the residues 17-21 is clearly more complex than that previously described by short MD simulations which only emphasized on the critical role of the side-chain of F19 in the potentiation effect of cholesterol on Aβ40 fibril formation [24].

The representative structures of the eight free energy minima designated as S1-S8 on the FEL with decreasing populations are shown in Figure 7. The first minimum, S1, with a population of 40%, is characterized by three chains forming well defined β-hairpins and chains A and B forming a short-twisted four-stranded antiparallel β-sheet. The β-strands cover residues L17-A21 and N27-I31 in chain A, residues L17-A21 and I32-V36 in chain B and residues H14-A21 and K28-M35 in chain C, indicating that the β-hairpins are distinct and extend beyond the CHC region by including some residues in the loop region (residues E22-K28). The second minimum (S2, population of 22%) displays two highly bent and twisted β-hairpins with strands covering residues L17-F20 and A30-I41 in chain A, residues L17-F20, K28-L34 and V39-I41 in chain C and a chain B adopting a five stranded β-sheet at R5-S8, Y10-Q15, L17-E22, I31-V36, and V39-I41. The orientation of the three chains is complex as the C-terminus residues (30-41) of chain A is antiparallel with the residues K28-L34 of chain C, and the CHC of chain C is parallel to residues V39-I41 of chain C.

The S3, S4, S6 and S7 minima, with a total population of 29%, have a β-strand content of 36% and share the same topological features: two peptides with highly twisted β-hairpins (chains A and B) covering distinct amino acids perpendicular to each other, and a third disordered and compact peptide (chain C) with no preferred interface with the other two chains. For instance, there is a small intermolecular β-sheet between residues V40-I41 and K28-G38 in S3, and between residues L34-V36 and V39-I41 in S4. The S5 state with a population of 7%, dominated by coil (49%) with low β-strand content (29%), forms a twisted four-stranded antiparallel β-sheet (chains A and B) with a chain C having very deformed β-strands. Finally, the S8 state with a population of 2% is amorphous and has high turn and coil contents of 36% and 38%. S8 has a low β-strand
content of 24% with strands at residues Q15-V18, L30-G33 and V40-I41 in chain A and residues V18-E22, I31-M35 and V39-V40 in chain B, and forms a short intermolecular parallel β-sheet between the two termini of chains A and B. In this state chain C is almost devoid of any secondary structure with the exception of a small helix covering residues A30-G33.

The various globular shapes of β-strand contents varying between 24% (S8 state) and 42% (S2) are all characterized by different orientations and packings of the chains, but also different binding interfaces and hot spots with cholesterol (Figure 8). In the S1 state, the residue hot spots with a trimer of cholesterol involve residues H14, K16, N27, K28 and I31 of chain A, and residues A2, F19, A21, D23, L30, I32 and L34 of chain B. In this state, cholesterol is located at the surface of the Aβ trimer (Figure 7, S1 state). The S3 state has almost the same residues for binding and a trimer of cholesterol located at the surface. The S5 state is also characterized by a trimer of cholesterol at the surface, but binding involves many residues in the region F4-L34 of chain B, the residues V24, K28, I31, I41 and A42 of chain A, and residues E22 and A42 of chain C (Figure 8). In total 57% of Aβ states have a trimer of cholesterol at the surface.

The S4, S6 and S7 states representing 22% of the ensemble are characterized by a dimer of cholesterol at the surface and one cholesterol inserted into the complex (Figure 7). In these states, we find that the residues 16-22 are essential for binding but many residues from the region K28 to A42 and a few residues in the N-terminus (residues F4, H6, S8 and V12) participate as well (Figure 8). Finally, the S2 state reveals a dimer of cholesterol inserted in Aβ trimer and one cholesterol at the surface with binding residues D1-F4, Q15-V24 and K28-I41, while the S8 state reveals a trimer of cholesterol in the interior of oligomers with binding residues D1-Q15 and K16-V36 (Figure 7 panel S8 and Figure 8).

Overall, these results are very different from a previous MD simulation which only reported on interactions of Aβ42 peptides bound to the surface cholesterol micelle [22]. Our cholesterol-binding sites of Aβ42 also differ from combined docking modelling and surface pressure studies of Langmuir monolayers that suggested the region E22-M35 as the minimal cholesterol binding site. It must be stressed that the docking procedure used the NMR structure of Aβ monomer mixed with detergent micelles, characterized by an alpha-helix at residues Q15-V36 with a hinge at residues G25-N27 [61]. It is interesting however that our simulations report on the contribution of the loop residues E22-K28 in the binding process.
To further understand the binding mechanism of cholesterol, Figure 9 reports the time-
averaged probability of contacts between the cholesterol molecules and the side chains of wild-
type (WT) Aβ42, A21G Aβ42 and A2T Aβ42 obtained from a total of 3 microseconds per system
(namely 1 microsecond MD simulation starting from the S1, S2 and S4 structures at 310 K shown
on Figure 7). Though the S1, S2 and S4 structures of WT represent 70% of the full conformation
space of the wild-type sequence, they were selected because they display hot spots involving the
N-terminus, the CHC, the loop region and the C-terminus. It is important to stress that these two
mutations were also selected because A2T introduces a hydrophilic residue and thus potentially
reduces the hydrophobic surface with cholesterol, and A21G because it reduces the total
hydrophobic character of the CHC region. Our results show that both mutations do not change the
profile of interactions per Aβ42 residue and per cholesterol in all regions of Aβ42, including the
N-terminus and the CHC region (Figure 9). These results suggest that cholesterol should
moderately alter the aggregation kinetics of Aβ42 A2T and Aβ42 A21G in the bulk solution, with
an enhancement that should be comparable to that observed for WT Aβ42 [22].

**Conclusions**

In summary, we have determined the Aβ42/cholesterol trimeric states by means of REST2
simulations and a force field designed for intrinsically disordered proteins. Consistent with our
previous REMD simulation using a force field for well-structured proteins [30], we found that the
conformational space does not contain the aggregation prone state of the parallel U- and S-shape
Aβ42 fibrils [2], but displays some antiparallel dimers with short intermolecular β-sheet built on
strands located at different regions that go beyond the CHC region. Based on our previous REMD
study of Aβ42/cholesterol dimers showing an increase of the population of β-hairpin and β-sheet
contents upon cholesterol addition compared to pure bulk solution [30], it is likely that all Aβ42
oligomers mixed with free cholesterol will show more β-structural and antiparallel β-sheet features
of the peptides than in pure bulk solution as the oligomer size augments [50].

Our simulations also show that the formation of Aβ42 trimers in the presence of cholesterol
involves two binding interfaces. A first highly populated one, where cholesterols are mainly
located at the surface of Aβ oligomers (either as a trimer or a dimer of cholesterol), and a second
much less populated where cholesterol is fully inserted in the interior of the Aβ oligomers. The
predominance of the first binding interface implies that cholesterol should act as a glue to speed
up the formation of larger aggregates and therefore catalyzes primary nucleation. This explains the
modest but non-negligible experimentally-observed acceleration of the aggregation rate of Aβ40
and Aβ42 in the presence of cholesterol [22].

Finally, we find that the binding hot spots with Aβ are very complex and go much beyond
the CHC region (residues L17-A21) and the hydrophobic residues L30-M35. Our mechanism
cannot be generalized to all amyloid proteins. Indeed, it has been shown that free cholesterol has
an inhibitory effect on aggregation of the 37-residue amylin protein which has two hydrophobic
regions LANFLV and FGAIL separated by HSSNN, in both solution and on model membranes
[62]. Clearly, a better understanding of the interactions of free cholesterols or free phospholipids
on amyloid aggregates in the brain either by computational [63] or experimental [25,64] means
must be further explored. To this end, we are coupling a coarse-grained protein force field in
aqueous solution [65,66] with coarse-grained cholesterols and phospholipids [67] to explore larger
aggregates and the impact of other disease-causing and disease-protecting mutations [68-71].

AUTHOR INFORMATION

Corresponding Author *E-mail : philippe.derreumaux@ibpc.fr

ORCID and E-mails

Phuong H. Nguyen: 0000-0003-1284-967X, phuong.nguyen@ibpc.fr

Trung Hai Nguyen: 0000-0003-1848-3963, nguyentrunghai@tdtu.edu.vn

Son Tung Ngo: 0000-0003-1034-1768, ngosontung@tdtu.edu.vn

Philippe Derreumaux: 0000-0001-9110-5585

Notes The authors declare no competing financial interest.

Author Contributions :

Conceptualization, Philippe Derreumaux; Data curation, Trung Hai Nguyen and Phuong H. Nguyen; Formal analysis, Trung Hai Nguyen, Phuong H. Nguyen, Son Tung Ngo and Philippe Derreumaux; Investigation, Phuong H. Nguyen; Methodology, Trung Hai Nguyen, Phuong H. Nguyen and Son Tung Ngo; Writing – review & editing, Philippe Derreumaux.

ACKNOWLEDGMENTS
This work was supported by grants from the French IDRIS and CINES computer centers, and the "Initiative d'Excellence" program from the French State, Grants "DYNAMO", ANR-11-LABX-0011, and "CACSICE", ANR-11-EQPX-0008.

References

1. Benson, M.D.; Buxbaum, J.N.; Eisenberg, D.S., Merlini, G.; Saraiva, M.J.M.; Sekijima, Y.; Sipe, J.D., Westermark, P. Amyloid Nomenclature 2018: Recommendations by the International Society of Amyloidosis (ISA) Nomenclature Committee. *Amyloid* **2018**, *25*, 215-219.

2. Doig, A.J.; Del Castillo-Frias, M.P.; Berthoumieu, O.; Tarus, B.; Nasica-Labouze, J.; Sterpone, F.; Nguyen, P.H.; Hooper, N.M.; Faller, P.; Derreumaux, P. Why Is Research on Amyloid-β Failing to Give New Drugs for Alzheimer's Disease? *ACS Chem. Neurosci.* **2017**, *8*, 1435-1437.

3. Nasica-Labouze, J.; Nguyen, P.H.; Sterpone, F.; Berthoumieu, O.; Buchete, N.V.; Coté, S.; De Simone, A.; Doig, A.J.; Faller, P.; Garcia, A. et al. Amyloid β Protein and Alzheimer's Disease: When Computer Simulations Complement Experimental Studies. *Chem. Rev.* **2015**, *115*, 3518-3563.

4. Panchal, M.; Loeper, J.; Cossec, J.C.; Perruchini, C.; Lazar, A.; Pompon, D.; Duyckaerts, C. Enrichment of Cholesterol in Microdissected Alzheimer's Disease Senile Plaques as Assessed by Mass Spectrometry. *J. Lipid Res.* **2010**, *51*, 598-605.

5. Gellermann, G.P.; Appel, T.R.; Davies, P.; Diekmann, S. Paired helical filaments contain small amounts of cholesterol, phosphatidylcholine and sphingolipids. *Biol Chem.* **2006**, *387*, 1267-1274.

6. Popp, J.; Meichsner, S.; Kölsch, H.; Lewczuk, P.; Maier, W.; Kornhuber, J.; Jessen, F.; Lütjohann, D. Cerebral and Extracerebral Cholesterol Metabolism and CSF Markers of Alzheimer's Disease. *Biochem. Pharmacol.* **2013**, *86*, 37-42.

7. Cutler, R.G.; Kelly, J.; Storie, K.; Pedersen, W.A.; Tammara, A.; Hatanpaa, K.; Troncoso, J.C.; Mattson, M.P. Involvement of oxidative stress-induced abnormalities in ceramide and cholesterol metabolism in brain aging and Alzheimer's disease. *Proc Natl Acad Sci U S A.* **2004**, *101*, 2070-2075.

8. Makin, S. The amyloid hypothesis on trial. *Nature*. **2018**, *559*, S4-S7.

9. Xiong, H.; Callaghan, D.; Jones, A.; Walker, D.G.; Lue, L.F.; Beach, T.G.; Sue, L.I.; Woulfe, J.; Xu, H.; Stanimirovic, D.B. et al. Cholesterol retention in Alzheimer's brain is responsible for high beta- and gamma-secretase activities and Abeta production. *Neurobiol Dis.* **2008**, *29*, 422-437.

10. Barrett, P.J.; Song, Y.; Van Horn, W.D.; Hustedt, E.J.; Schafer, J.M.; Hadziselimovic, A.; Beel, A.J.; Sanders, C.R. The amyloid precursor protein has a flexible transmembrane domain and binds cholesterol. *Science.* **2012**, *336*, 1168-1171.

11. Pantelopulos, G.A.; Panahi, A.; Straub, J.E. Impact of Cholesterol Concentration and Lipid Phase on Structure and Fluctuation of Amyloid Precursor Protein. *J Phys Chem B.* **2020**, *124*, 10173-10185.
12. Wang, H.; Kulas, J.A.; Wang, C.; Holtzman, D.M.; Ferris, H.A.; Hansen, S.B. Regulation of beta-amyloid production in neurons by astrocyte-derived cholesterol. *Proc Natl Acad Sci USA.* **2021**, *118*, e2102191118.

13. Takeuchi, S.; Ueda, N.; Suzuki, K.; Shimozawa, N.; Yasutomi, Y.; Kimura, N. Elevated Membrane Cholesterol Disrupts Lysosomal Degradation to Induce β-Amyloid Accumulation: The Potential Mechanism Underlying Augmentation of β-Amyloid Pathology by Type 2 Diabetes Mellitus. *Am J Pathol.* **2019**, *189*, 391-404.

14. Corder, E.H.; Saunders, A.M.; Strittmatter, W.J.; Schmechel, D.E.; Gaskell, P.C.; Small, G.W.; Roses, A.D.; Haines, J.L.; Pericak-Vance, M.A. Gene dose of apolipoprotein E type 4 allele and the risk of Alzheimer's disease in late onset families. *Science.* **1993**, *261*, 921-923.

15. Di Scala, C.; Yahi, N.; Boutemeur, S.; Flores, A.; Rodriguez, L.; Chahinian, H.; Fantini, J. Common Molecular Mechanism of Amyloid Pore Formation by Alzheimer's β-amyloid Peptide and α-synuclein. *Sci. Rep.* **2016**, *6*, 28781.

16. Banerjee, S.; Hashemi, M.; Zagorski, K.; Lyubchenko, Y.L. Cholesterol in Membranes Facilitates Aggregation of Amyloid β Protein at Physiologically Relevant Concentrations. *ACS Chem Neurosci.* **2021**, *12*, 506-516.

17. Habchi, J.; Chia, S.; Galvagnion, C.; Michaels, T.C.T.; Bellaiche, M.M.J.; Ruggeri, F.S.; Sanguanini, M.; Idini, I.; Kumita, J.R.; Sparr, E et al. Cholesterol Catalyses Aβ42 Aggregation through a Heterogeneous Nucleation Pathway in the Presence of Lipid Membranes. *Nat. Chem.* **2018**, *10*, 673-683.

18. Lockhart, C.; Klimov, D.K. Cholesterol Changes the Mechanisms of Aβ Peptide Binding to the DMPC Bilayer. *J. Chem. Inf. Model.* **2017**, *57*, 2554-2565.

19. Grouleff, J.; Irudayam, S.J.; Skeby, K.K.; Schiøtt, B. The Influence of Cholesterol on Membrane Protein Structure, Function, and Dynamics Studied by Molecular Dynamics Simulations. *Biochim. Biophys. Acta.* **2015**, *1848*, 1783-1795.

20. Yu, X.; Zheng, J. Cholesterol Promotes the Interaction of Alzheimer β-amyloid Monomer with Lipid Bilayer. *J. Mol. Biol.* **2012**, *421*, 561-571.

21. Dias, C.L.; Jalali, S.; Yang, Y.; Cruz, L. Role of Cholesterol on Binding of Amyloid Fibrils to Lipid Bilayers. *J. Phys. Chem. B.* **2020**, *124*, 3036-3042.

22. Elbassal, E.A.; Liu, H.; Morris, C.; Wojcikiewicz, E.P.; Du, D. Effects of Charged Cholesterol Derivatives on Aβ40 Amyloid Formation. *J Phys Chem B.* **2016**, *120*, 59-68.

23. Harris, J.R.; Milton, N.G. Cholesterol in Alzheimer's Disease and Other Amyloidogenic Disorders. *Subcell. Biochem.* **2010**, *51*, 47-75.

24. Zhou, X.; Xu, J. Free Cholesterol Induces Higher β-sheet Content in Aβ Peptide Oligomers by Aromatic Interaction with Phe19. *PLoS One* **2012**, *7*, e46245.

25. Ait-Bouziad, N.; Lv, G.; Mahul-Mellier, A.L.; Xiao, S.; Zorludemir, G.; Eliezer, D.; Walz, T.; Lashuel, H.A. Discovery and Characterization of Stable and Toxic Tau/phospholipid Oligomeric Complexes. *Nat. Commun.* **2017**, *8*, 1678.

26. Sciacca, M.F.; Kotler, S.A.; Brender, J.R.; Chen, J.; Lee, D.K.; Ramamoorthy, A. Two-step mechanism of membrane disruption by Aβ through membrane fragmentation and pore formation. *Biophys J.* **2012**, *103*, 702-710.
27. Nguyen, P.H.; Ramamoorthy, A.; Sahoo, B.R.; Zheng, J.; Faller, P.; Straub, J.E., Dominguez, L.; Shea, J.E.; Dokholyan, N.V.; De Simone, A. et al. Amyloid Oligomers: A Joint Experimental/Computational Perspective on Alzheimer's Disease, Parkinson's Disease, Type II Diabetes, and Amyotrophic Lateral Sclerosis. *Chem Rev* **2021**, *121*, 2545-2647.

28. Scollo, F.; Tempra, C.; Lolicato, F.; Sciacca, M.F.M.; Raudino, A.; Milardi, D.; La Rosa, C. Phospholipids Critical Micellar Concentrations Trigger Different Mechanisms of Intrinsically Disordered Proteins Interaction with Model Membranes. *J Phys Chem Lett.* **2018**, *9*, 5125-5129.

29. Sciacca, M.F.; Lolicato, F.; Tempra, C.; Scollo, F.; Sahoo, B.R.; Watson, M.D.; García-Viñuales, S.; Milardi, D.; Raudino, A., Lee, J.C. et al. Lipid-Chaperone Hypothesis: A Common Molecular Mechanism of Membrane Disruption by Intrinsically Disordered Proteins. *ACS Chem Neurosci.* **2020**, *11*, 4336-4350.

30. Ngo, S.T.; Nguyen, P.H.; Derreumaux, P. Cholesterol Molecules Alter the Energy Landscape of Small Aβ1-42 Oligomers. *J Phys Chem B.* **2021**, *125*, 2299-2307.

31. Wang, L.; Friesner, R.A.; Berne, B.J. Replica exchange with solute scaling: a more efficient version of replica exchange with solute tempering (REST2). *J Phys Chem B.* **2011**, *115*, 9431-9438.

32. Smith, A.K.; Lockhart, C.; Klimov, D.K. Does Replica Exchange with Solute Tempering Efficiently Sample Aβ Peptide Conformational Ensembles? *J Chem Theory Comput.* **2016**, *12*, 5201-5214.

33. Noda, K.; Tachi, Y.; Okamoto, Y. Structural Characteristics of Monomeric Aβ42 on Fibril in the Early Stage of Secondary Nucleation Process. *ACS Chem Neurosci.* **2020**, *11*, 2989-2998.

34. Huang, J.; Rauscher, S.; Nawrocki, G.; Ran, T.; Feig, M.; de Groot, B.L.; Grubmüller, H.; MacKerell, A.D. Jr. CHARMM36m: An Improved Force Field for Folded and Intrinsically Disordered Proteins. *Nat. Methods* **2017**, *14*, 71-73.

35. Man, V.H.; He, X.; Derreumaux, P.; Ji, B.; Xie, X.Q.; Nguyen, P.H.; Wang J. Effects of All-Atom Molecular Mechanics Force Fields on Amyloid Peptide Assembly: The Case of Aβ<sub>16-22</sub> Dimer. *J. Chem. Theory Comput.* **2019**, *15*, 1440-1452.

36. Samantray, S.; Yin, F.; Kav, B; Strodel, B. Different Force Fields Give Rise to Different Amyloid Formation Pathways in Molecular Dynamics Simulations. *J. Chem. Inf. Model.* **2020**, *60*, 6462-6475.

37. Nguyen, P.H.; Derreumaux, P. Structures of the intrinsically disordered Aβ, tau and α-synuclein proteins in aqueous solution from computer simulations. *Biophys Chem.* **2020**, *264*, 106421.

38. Nguyen, P.H.; Campanera, J.M.; Ngo, S.T.; Loquet, A.; Derreumaux, P. Tetrameric Aβ40 and Aβ42 β-Barrel Structures by Extensive Atomistic Simulations. II. In Aqueous Solution. *J Phys Chem B.* **2019**, *123*, 6750-6756.

39. Strodel, B. Amyloid aggregation simulations: challenges, advances and perspectives. *Curr Opin Struct Biol.* **2021**, *67*, 145-152.

40. Huet, A.; Derreumaux, P. Impact of the mutation A21G (Flemish variant) on Alzheimer's beta-amyloid dimers by molecular dynamics simulations. *Biophys J.* **2006**, *91*, 3829-3840.
41. Benilova, I.; Gallardo, R.; Ungureanu, A.A.; Castillo Cano, V.; Snellinx, A.; Ramakers, M., Bartic, C.; Rousseau, F.; Schymkowitz, J.; De Strooper, B. The Alzheimer disease protective mutation A2T modulates kinetic and thermodynamic properties of amyloid-β (Aβ) aggregation. *J Biol Chem.* **2014**, *289*, 30977-30989.

42. Zheng, X.; Liu, D.; Roychaudhuri, R.; Teplow, D.B.; Bowers, M.T. Amyloid β-Protein Assembly: Differential Effects of the Protective A2T Mutation and Recessive A2V Familial Alzheimer's Disease Mutation. *ACS Chem Neurosci.* **2015**, *6*, 1732-1740.

43. Nguyen, P.H.; Sterpone, F.; Pouplana, R.; Derreumaux, P.; Campanera, J.M. Dimerization Mechanism of Alzheimer Aβ₄₀ Peptides: The High Content of Intrapeptide-Stabilized Conformations in A2V and A2T Heterozygous Dimers Retards Amyloid Fibril Formation. *J Phys Chem B.* **2016**, *120*, 12111-12126.

44. Sharma, B.; Ranganathan, S.V.; Belfort, G. Weaker N-Terminal Interactions for the Protective over the Causative Aβ Peptide Dimer Mutants. *ACS Chem Neurosci.* **2018**, *9*, 1247-1253.

45. Aggarwal, L.; Biswas, P. Hydration Thermodynamics of the N-Terminal FAD Mutants of Amyloid-β. *J Chem Inf Model.* **2021**, *61*, 298-310.

46. Yang, X.; Meisl, G.; Frohm, B.; Thulin, E.; Knowles, T.P.J; Linse, S. On the role of sidechain size and charge in the aggregation of Aβ42 with familial mutations. *Proc Natl Acad Sci U S A.* **2018**, *115*, E5849-E5858.

47. Wang, J.; Wang, W.; Kollman, P. A.; Case, D. A., Automatic Atom Type and Bond Type Perception in Molecular Mechanical Calculations. *J. Mol. Graph. Model.* **2006**, *25*, 247-260.

48. Wang, J.; Wolf, R. M.; Caldwell, J. W.; Kollman, P. A.; Case, D. A., Development and Testing of a General Amber Force Field. *J. Comput. Chem.* **2004**, *25*, 1157-1174.

49. Cao, Y.; Jiang, X.; Han, W. Self-Assembly Pathways of β-Sheet-Rich Amyloid-β(1-40) Dimers: Markov State Model Analysis on Millisecond Hybrid-Resolution Simulations. *J. Chem. Theory Comput.* **2017**, *13*, 5731-5744.

50. Zheng, W.; Tsai, M.Y.; Wolynes, P.G. Comparing the Aggregation Free Energy Landscapes of Amyloid Beta(1-42) and Amyloid Beta(1-40). *J. Am. Chem. Soc.* **2017**, *139*, 16666-16676.

51. Kreutzer, A.G.; Nowick, J.S. Elucidating the Structures of Amyloid Oligomers with Macrocyclic β-Hairpin Peptides: Insights into Alzheimer's Disease and Other Amyloid Diseases. *Acc. Chem. Res.* **2018**, *51*, 706-718.

52. Phillips, J.C., Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R.D.; Kalé, L.; Schulten, K. Scalable molecular dynamics with NAMD. *J Comput Chem.* **2005**, *26*, 1781-1802.

53. Man, V.H.; Nguyen, P.H.; Derreumaux, P. High-Resolution Structures of the Amyloid-β 1-42 Dimers from the Comparison of Four Atomistic Force Fields. *J. Phys. Chem. B.* **2017**, *121*, 5977-5987.

54. van der Spoel, D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A. E.; Berendsen, H. J. GROMACS: Fast, Flexible, and Free. *J. Comput. Chem.* **2005**, *26*, 1701-1718.

55. Bussi, G.; Donadio, D.; Parrinello, M. Canonical sampling through velocity rescaling. *J Chem Phys.* **2007**, *126*, 014101.
56. Kabsch, W.; Sander, C. Dictionary of Protein Secondary Structure: Pattern Recognition of Hydrogen-bonded and Geometrical Features. *Biopolymers* **1983**, *22*, 2577-2637.

57. Lu, Y.; Derreumaux, P.; Guo, Z.; Mousseau, N.; Wei, G. Thermodynamics and dynamics of amyloid peptide oligomerization are sequence dependent. *Proteins*. **2009**, *75*, 954-63.

58. Smith, L.J.; Daura, X.; van Gunsteren, W.F. Assessing Equilibration and Convergence in Biomolecular Simulations. *Proteins* **2002**, *48*, 487-96.

59. Plesnar, E., Subczynski, W.K., Pasenkiewicz-Gierula, M. Comparative computer simulation study of cholesterol in hydrated unary and binary lipid bilayers and in an anhydrous crystal. *J Phys Chem B.* **2013**, *117*, 8758-8769.

60. Huang K, García AE. Acceleration of Lateral Equilibration in Mixed Lipid Bilayers Using Replica Exchange with Solute Tempering. *J Chem Theory Comput.* **2014**, *10*, 4264-4272.

61. Di Scala, C.; Yahi, N.; Lelièvre, C.; Garmy, N.; Chahinian, H.; Fantini, J. Biochemical identification of a linear cholesterol-binding domain within Alzheimer's β amyloid peptide. *ACS Chem Neurosci.* **2013**, *4*, 509-517.

62. Cho, W.J.; Trikha, S.; Jeremic, A.M. Cholesterol regulates assembly of human islet amyloid polypeptide on model membranes. *J Mol Biol.* **2009**, *393*, 765-775.

63. Fatafta, H.; Kav, B.; Bundschuh, B.F.; Loschwitz, J.; Strodel, B. Disorder-to-order transition of the amyloid-β peptide upon lipid binding. *Biophys Chem.* **2022**, *280*, 106700.

64. Confer, M.P.; Holcombe, B.M.; Foes, A.G.; Holmquist, J.M.; Walker, S.C.; Deb, S.; Ghosh, A. Label-Free Infrared Spectroscopic Imaging Reveals Heterogeneity of β-Sheet Aggregates in Alzheimer's Disease. *J Phys Chem Lett.* **2021**, *12*, 9662-9671.

65. Maupetit, J., Tuffery, P., Derreumaux, P. A coarse-grained protein force field for folding and structure prediction. *Proteins.* **2007**, *69*, 394-408.

66. Sterpone, F.; Melchionna, S.; Tuffery, P.; Pasquali, S.; Mousseau, N.; Cragnolini, T.; Chebaro, Y.; St-Pierre, J.F.; Kalimeri, M.; Barducci A. et al. The OPEP protein model: from single molecules, amyloid formation, crowding and hydrodynamics to DNA/RNA systems. *Chem Soc Rev.* **2014**, *43*, 4871-4893.

67. Brandner, A.F.; Timr, S.; Melchionna, S.; Derreumaux, P.; Baaden, M.; Sterpone, F. Modelling lipid systems in fluid with Lattice Boltzmann Molecular Dynamics simulations and hydrodynamics. *Sci Rep.* **2019**, *9*, 16450.

68. Viet, M.H.; Nguyen, P.H.; Derreumaux, P., Li, M.S. Effect of the English familial disease mutation (H6R) on the monomers and dimers of Aβ40 and Aβ42. *ACS Chem Neurosci.* **2014**, *5*, 646-657.

69. Nguyen, P.H., Tarus, B.; Derreumaux, P. Familial Alzheimer A2 V mutation reduces the intrinsic disorder and completely changes the free energy landscape of the Aβ1-28 monomer. *J Phys Chem B.* **2014**, *118*, 501-510.

70. Viet, M.H.; Nguyen, P.H., Ngo, S.T.; Li, M.S.; Derreumaux, P. Effect of the Tottori familial disease mutation (D7N) on the monomers and dimers of Aβ40 and Aβ42. *ACS Chem Neurosci.* **2013**, *4*, 1446-1457.

71. Côté, S.; Laghaei, R.; Derreumaux, P.; Mousseau, N. Distinct dimerization for various alloforms of the amyloid-beta protein: Aβ(1-40), Aβ(1-42), and Aβ(1-40)(D23N). *J Phys Chem B.* **2012**, *116*, 4043-4055.
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Initial REST2 conformation of Aβ42 trimer - 3 cholesterols with Aβ42 forming β-sheets at residues 13-21 and 29-35 of chain A (red), residues 16-22 and 31-37 of chain B (green), and residues 13-21 and 27-35 of chain C (yellow). Each N-terminus is shown by a blue ball. The three all-atom cholesterol molecules are shown in green.</figcaption></figure>
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** REST2 Sampling. (A) Overlap of the total potential energies for replicas 1 to 16. (B and C) Exchanges of coordinates as a function of simulation time of replica 1 and replica 16, respectively.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Free energy landscapes of Aβ42 trimer + 3 cholesterol from PCA. REST2 FEL by using the single REST2 trajectory at 310 K and time intervals 50-220 ns (A) and 50-250 ns (B). REMD FEL (C) and REST2 FEL (D) by using the combined REST2 (310 K) and REMD (315 K) trajectories and the time interval 50-250 ns. Note that the PC1 and PC2 differ from REMD to REST2.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Secondary structure composition along the amino acid sequence using the time interval 50-250 ns of the REMD simulation at 315 K and the time interval 50-250 ns of the REST2 simulation at 310 K. (A) beta-strand, (B) turn and (C) coil content.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Probability of side chain contacts between each cholesterol molecule and each amino acid of each Aβ42 peptide using the time interval 50-250 ns of REMD at 315 K, and REST2 at 310 K.</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Probability of contacts between the cholesterol molecules in dimer+monomer form and trimer form and the side-chains of Aβ42 using the time intervals 50-250 ns of REST2 at 310 K.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Representative structures and populations of eight free energy minima from REST2 simulation at 310 K. Chain A is in red, chain C is in white, and chain B is in pastel pink. The ball shows the N-terminus of each chain. The cholesterols are visualised in all-atom and van der Waals representations.</figcaption></figure>
<figure><img src="image_11.png" /><figcaption>**Figure 8.** Contact between each Aβ42 chain (A, B, C) and cholesterol. Shown are results for each free energy minimum Si (i=1...8). One Aβ42 side-chain and one cholesterol molecule are considered in contact if there is at least one intermolecular distance below 0.45 nm.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Figure 9.** Time-averaged probability of contacts between the cholesterol molecules and the side-chain of Aβ42 obtained from 1 microsecond MD trajectory at 310 K of the wild-type (WT, black), A2T mutation (red) and A21G mutation (green) sequence. Shown are the averaged values over the three simulations starting from the S1, S2 and S4 structures displayed in Figure 7.</figcaption></figure>