Curcumin Inhibits the Primary Nucleation of Amyloid-Beta Peptide: A Molecular Dynamics Study.

The amyloid plaques are a key hallmark of neurodegenerative diseases such as Alzheimer's disease and Parkinson's disease. Amyloidogenesis is a complex long-lasting multiphase process starting with the formation of nuclei of amyloid peptides: a process assigned as a primary nucleation. Curcumin (CU) is a well-known inhibitor of the aggregation of amyloid-beta (Aβ) peptides. Even more, CU is able to disintegrate preformed Aβ firbils and amyloid plaques. Here, we simulate by molecular dynamics the primary nucleation process of 12 Aβ peptides and investigate the effects of CU on the process. We found that CU molecules intercalate among the Aβ chains and bind tightly to them by hydrogen bonds, hydrophobic, π-π, and cation-π interactions. In the presence of CU, the Aβ peptides form a primary nucleus of a bigger size. The peptide chains in the nucleus become less flexible and more disordered, and the number of non-native contacts and hydrogen bonds between them decreases. For comparison, the effects of the weaker Aβ inhibitor ferulic acid (FA) on the primary nucleation are also examined. Our study is in good agreement with the observation that taken regularly, CU is able to prevent or at least delay the onset of neurodegenerative disorders.

1. Introduction

Amyloidogenic peptides and proteins are soluble structures that are able to aggregate into insoluble amyloid fibrils with a characteristic core structure rich in β-sheets [1]. Evolutionary, these proteins have appeared in the microorganisms as necessary components of physiological functions such as biofilm formation [2], cell adhesion [3], interactions with host membranes [4], and permeabilization [5]. In humans, some of them also show physiological functions. For example, amylin is involved in the regulation of insulin and glucagon secretion [6], tau protein stabilizes the microtubules in the cells [7], the amyloid-beta peptide controls the synaptic activity [8] and acts as a protective agent against brain infections [9,10], α-synuclein regulates the neurotransmitter release [11], and premalanosome protein (PMEL) forms fibrils in the melanocytes required for the optimal cell function [12]. At the same time, the ability of amyloidogenic proteins to form fibrils is a key feature of more than 50 diseases such as Alzheimer’s disease [13], Parkinson’s disease [14], prion disease [15], type 2 diabetes [16], and sickle-cell anemia [17]. The triggers that convert the functional amyloids into pathological ones still remain unclear.

The amyloid fibril formation is a complex process consisting of three phases: nucleation, elongation, and saturation [18]. The nucleation starts with the aggregation of two monomers into a dimer, which further attracts other monomers and evolves into an oligomer. The oligomer might arrange into toxic fibrillar aggregate (on-pathway) or into non-toxic disordered aggregate (off-pathway). The toxic fibrils elongate by attaching novel monomers in an appropriate cross-β-conformation. In the saturation phase, the fibril formation reaches equilibrium as monomers detach and attach continuously at the ends of the fibril. The nucleation could be primary or secondary. A primary nucleation is assigned to the creation of a folding nucleus, from which oligomeric and fibrillar assemblies emanate [19], while the secondary nucleation describes the formation of nuclei around pre-existing crystals [20].

As the nucleation phase is the initial step in the process of fibril formation, inhibition of the aggregation of amyloidogenic peptides and proteins is one of the therapeutic strategies in the development of anti-amyloid drug candidates [21]. It was reported that curcumin (CU), rosmarinic acid, nordihydroquaiaretic acid, rifampicin, and wine-related polyphenols, such as myricetin, dose-dependently inhibit amyloid-beta (Aβ) fibril formation from Aβ1-40 and Aβ1-42 and destabilize preformed Aβ fibrils in vitro [22].

CU, a polyphenol from Curcuma longa (turmeric plant), is known to treat many diseases, including neurodegenerative ones such as Alzheimer’s disease, Parkinson’s disease, prion disease, multiple sclerosis, stroke, anxiety, and depression [23]. It has been observed that the death rate of Alzheimer’s/dementia per 100,000 population standardized by age in India is significantly lower than in the Middle East [24]. One of the reasons could be the regular intake of turmeric in the routine diet of Indians [25]. CU is a widely studied compound, and many good reviews on its properties and applications are available in the literature [23,26,27,28,29,30]. It was found in vitro that CU blocks the toxicity of Aβ oligomers (IC50 = 0.8 µM) [31] and disintegrates preformed Aβ fibrils (EC50 = 1 µM) [32]. In animals, it was shown that CU reduces insoluble Aβ deposits [33] and amyloid plaques [31]. CU formulations with improved bioavailability have been demonstrated to improve memory, mood, and alertness in clinical trials with non-demented elderly patients [34,35].

The interactions of CU with Aβ peptides and fibrils have been simulated by molecular dynamics (MD). Zhao et al. have studied the effects of CU on the stability of Aβ dimers and have found that CU acts as β-sheet breaker, reducing the β-sheet content within the Aβ oligomers [36]. Awasthi et al. have analyzed the effects of CU on the stability of two familial Aβ42 mutations—A2V (harmful) and A2T (protective)—and have found that the mutant A2V is more stable than A2T in the absence of CU, while A2T is more stable in the presence of CU [37]. Apart from oligomers, CU binds tightly to preform Aβ fibrils occupying the binding pocket inside fibrils [38]. It forms hydrogen bonds and hydrophobic interactions with protofibrils [38,39] and causes structural distortion and perturbation [39,40], leading to the dissociation of the outermost peptides [41]. The number of intermolecular backbone hydrogen bonds between adjacent β-threats in fibrils correlates with the anti-aggregation activity of Aβ inhibitors [42].

The aim of the present study was to simulate by MD the primary nucleation of an ensemble of 12 Aβ1-42 peptides and to analyze the effects of the inhibitor curcumin. For comparison, the effects of the weaker inhibitor ferulic acid (FA) (IC50 = 5.5 µM) [43] were also conducted.

2. Models and Method

2.1. Modeled Ligands and Systems

The chemical structures of CU (CID 969516) and FA (CID 445858) were retrieved from PubChem [44] (Figure 1). The keto-enol form of CU was used in the present study instead of β-diketone, as it has been proven that the keto-enol form exhibits anti-aggregatation properties and is responsible for the ability of CU to prevent Aβ oligomer and fibril assembly [45,46]. FA is a weak acid with pKa 4.58 and at physiological pH 7.4, the anionic form dominates. This form was used in the simulations. The molecules were parameterized by GAFF 2.11 [47] with AM1–BCC charges [48] before being used in Amber 18 [49]. The Aβ1-42 monomer structure was obtained from PDB (Protein Database); (pdb code 1IYT) [50]. Both ends were capped, and all hydrogens were added.

Figure 1

Chemical structures of the amyloid-beta (Aβ) inhibitor (A) curcumin (CU) and (B). non-inhibitor ferulic acid (FA).

Five systems were modeled (Figure 2). The first one contained 12 Aβ1-42 monomers randomly positioned in a truncated octahedron cell with a volume of 1,927,025 Å3 and solvated with TIP3P water yielding 10.4 mM solution. NaCl was added to maintain a physiological salt concentration and adjusted to neutralize the peptides’ charge. This system was used for the construction of the other four systems. Two of them contained 12 and 36 CU molecules spread randomly around the Aβ monomers corresponding to 9.3 mM and 30.4 mM solutions, respectively. The last two systems included 12 and 36 FA anions placed around the Aβ monomers, giving 10.3 mM and 31.2 mM solutions, respectively. The salt concentration was adjusted additionally to neutralize the FA charge.

Figure 2

Modeled systems in the present study: (A) 12 Aβ monomers; (B) 12 Aβ monomers and 12 CU molecules; (C) 12 Aβ monomers and 36 CU molecules; (D) 12 Aβ monomers and 12 FA anions; (E) 12 Aβ monomers and 36 FA anions.

2.2. MD Protocol

Initially, the solvated systems were energy minimized for 5000 steps with 3 kcal/molÅ2 harmonic restraints on solute heavy atoms, followed by heating them from 0 to 300 K over a period of 1 ns at constant volume with identical restraints. Next, the systems were equilibrated for 1 ns of constant pressure density with restraints and for 100 ns without any restraints. Finally, they were simulated for 1000 ns of production dynamics under constant temperature (310 K) and pressure (1 bar), maintained with the Langevin thermostat [51] and Berendsen barostat [52], respectively. The systems were simulated with the ff14SB [53] force field under periodic boundary conditions. A 12.0 Å cutoff was used for both van der Waals and electrostatic interactions; long-range electrostatics beyond the real-space cutoff were evaluated with the particle-mesh Ewald (PME) scheme [54]. During heating, density equilibration, preproduction, and production dynamics, covalent bonds involving hydrogen were constrained using the SHAKE algorithm [55], allowing a 2 fs time step; only during energy minimization were the bonds to hydrogen were not constrained. During production dynamics, frames were saved every 1 ns for a total of 1000 per trajectory to be used in subsequent analysis. The MD simulations were performed by Amber 18 [56] and were analyzed by cpptraj V4.24.0 [57].

3. Results

The five modeled Aβ ensembles—12 Aβ1-42 monomers (10.4 mM), 12 Aβ1-42 monomers + 12 CU molecules (9.3 mM), 12 Aβ1-42 monomers + 36 CU molecules (30.4 mM), 12 Aβ1-42 monomers + 12 FA anions (10.3 mM), and 12 Aβ1-42 monomers + 36 FA anions (31.2 mM) (Figure 2)—were simulated in saliva for 1 µs according to the MD protocol described in Models and Methods. The modeled Aβ ensembles formed stable nuclei for a different period of time. The coordinate trajectories were processed and used to analyze the effects of CU and FA on the process of primary nucleation of Aβ peptides.

3.1. Curcumin Stabilizes Immediately the Aβ1-42 Ensemble

The averaged backbone RMSDs (root mean square deviations) of the ensembles of 12 Aβ1-42 peptides with and without ligands were calculated as a function of time (Figure S1). The ensemble without a ligand underwent conformational distortions during the initial 300 frames (300 ns) and then was stabilized (Figure S1A). The presence of a ligand stabilizes the Aβ ensemble earlier. In the presence of 12 CU molecules, the ensemble was stabilized for 8 ns (Figure S1C); in the presence of 36 CU, it was stabilized for 2 ns (Figure S1E); in the presence of 12 FA, it was stabilized for 20 ns (Figure S1G); and in the presence of 36 FA, it was stabilized for 15 ns (Figure S1I). It is clearly evident that the presence of CU molecules stabilizes immediately the ensemble in a dose-dependent manner. Two additional 50 ns productions on the five systems were run. They confirmed the immediate stabilizing effect of CU (Figure S2).

The backbone RMSFs (root mean square fluctuations) indicate the most fluctuating Aβ residues during the MD simulation. Typically, the terminals of the Aβ peptides fluctuate more than the middle parts, with the N-terminals being the most flexible. The averaged RMSFs of the ensembles of 12 Aβ peptides with and without ligands are summarized in Figure 3A. The presence of any ligand reduces the fluctuations of the Aβ ensemble. The effects of CU are more prominent than the effects of FA. CU decreases mainly the fluctuations of the C-terminals (Figure S1D,F).

Figure 3

Root mean square fluctuations (RMSF) (A), solvent-accessible surface area (SASA) (B), number of non-native contacts (C), and number of H-bonds (D) averaged over 1000 frames (1000 ns) for 12 Aβ monomers (gray), 12 Aβ monomers and 12 CU molecules (pink), 12 Aβ monomers and 36 CU molecules (red), 12 Aβ monomers and 12 FA anions (light green), 12 Aβ monomers and 36 FA anions (dark green).

3.2. Curcumin Binds inside the Aβ Nucleus

The coordinates of the modeled systems after 1 µs MD simulations are given in Figure 4. In 300 ns, the 12 Aβ monomers collapsed into a compact nucleus losing part of the ordered structures (α-helix, 3-10 helix, β-turn) (Figure S2). The Aβ1-42 peptide is a water-soluble molecule containing 20 hydrophobic residues (Ala, Ile, Leu, Met, Val, Phe, Tyr, and Trp) clustered mainly in the middle part of the molecule and in the C-terminal. The Aβ1-42 peptide has an overall logP −12.2 [44] and isoelectric point pI 5.5 [58]. In polar medium, the flexible peptide chains adopt spontaneously a globular conformation with a polar surface and non-polar core (Figure 4A). As a typical hydrophobic interaction, the nucleation is an entropy-driven process governed by an increase of solvent entropy.

Figure 4

The modeled systems after 1 µs molecular dynamics (MD) simulation: (A) 12 Aβ monomers; (B) 12 Aβ monomers and 12 CU molecules; (C) 12 Aβ monomers and 36 CU molecules; (D) 12 Aβ monomers and 12 FA anions; (E) 12 Aβ monomers and 36 FA anions.

Regardless of the concentrations, all molecules of the strong inhibitor CU were bound to Aβ peptides (Figure 4B,C). Some of them were intercalated and bound inside the core, while others covered the surface. The main interactions between CU and Aβ residues were hydrophobic, but π–π interactions with the side chains of Phe4, Tyr10, Phe19, and Phe20 and cation–π interactions with the side chains of Arg5, Lys16, and Lys28 also occurred. Plenty of H-bonds exist between CU molecules and Aβ residues. The hydrogen bonds in the ensembles are considered in detail in Section 3.6.

The Aβ ensemble has quite a different look in the presence of the weak inhibitor FA (Figure 4D,E). In the presence of 12 FA, after 1 µs simulation, nine of them (75%) were bound to Aβ peptides, while in the presence of 36 FA, this number was only 18 (50%). Here again, the main interactions between the FA molecules and Aβ peptides were hydrophobic. Most of the bound molecules made π–π stacking with the aromatic side chains of Phe4, Tyr10, Phe19, and Phe20. Cation–π interactions appeared between the phenyl ring of FA and the side chains of Arg5, Lys16, and Lys28. The hydrogen bonds between Aβ and FA anions are considered in detail in Section 3.6.

Both ligands disrupted the salt bridges formed between the Aβ chains. Their number decreases from 142 (without ligands) to 121 (with 12 CU), 119 (with 36 CU), 123 (with 12 FA), and 118 (with 36 FA). The salt bridges were formed between Asp1, Asp7, Asp23, Glu3, Glu22, Arg5, Lys16, and Lys28.

3.3. Curcumin Does Not Affect the Secondary Structures of the Peptides in the Aβ1-42 Ensemble

The analysis of the secondary structures of the Aβ peptides indicates that the ensembles with and without ligands underwent similar changes during the simulations. The total structural propensities over the initial and final 100 ns are given in Figure S3. The propensities of helices (including 3–10 helices) decreased by 8% (12 Aβ) to 18% (12 Aβ + 12 FA). The β-turn propensities decreased for the ensembles of 12 Aβ peptides, 12 Aβ + 12 CU, and 12 Aβ + 36 FA, for the ensemble of 12 Aβ + 12 FA, they increased, and the ensemble of 12 Aβ + 36 CU remained almost unchanged. The bending propensities in all ensembles increased by 13% (12 Aβ + 36 FA) to 23% (12 Aβ + 36 CU). Indeed, the nucleation proceeds by bending and losing the helical and β-turn structures (Figure 4). The final ensembles are nuclei of randomly folded chains. A distinguishable effect of CU on the secondary structures of Aβ peptides during the 1 µs simulation of the process of nucleation was not observed.

3.4. Curcumin Increases the Solvent-Accessible Surface Area (SASA) of the Aβ1-42 Ensemble

The solvent-accessible surface area (SASA) was calculated using the algorithm of Weiser et al. [59] and expressed in Å2. The SASAs of the five ensembles averaged over 1000 frames are given in Figure 3B. CU at higher concentration (36 CU) increases the SASA of the 12 Aβ ensemble by 20%. At lower concentration (12 CU), the increase is only 4% comparable to the increase caused by 12 FA and 36 FA. The increased surface area is associated by the many CU molecules bound among the folded Aβ chains (Figure 4).

3.5. Curcumin Decreases the Number of Non-Native Contacts in the Aβ1-42 Ensemble

Amber defines as a non-native contact any contact between a pair of atoms satisfying the distance cutoff (by default 7 Å), which is not already a native contact [56]. The number of the non-native contacts between the Aβ monomers in the five ensembles averaged over 1000 frames are given at Figure 3C. The CU molecules at a higher concentration decrease the number of the non-native contacts in the Aβ ensemble by 16%, while the decrease caused by other ligands varied between 2% for 12 FA and 4% for 12 CU and 36 FA. The reduced contacts between the Aβ monomers are due to the intercalation of many CU molecules bound inside the Aβ nucleus (Figure 4).

3.6. Curcumin Decreases the Number of Hydrogen Bonds in the Aβ1-42 Ensemble

The number of hydrogen bonds averaged over 1000 frames of the five Aβ ensembles are plotted at Figure 3D. Here again, CU at higher concentration decreases the number of H-bonds the most. The CU molecules bound internally disrupt the pre-existing H-bonds between the Aβ monomers and/or prevent the formation of new ones (Figure 4B,C). As expected, FA as an anion and H-bond acceptor also decreases the number of H-bonds by replacing the peptide residues in some of them (Figure 4D,E).

In total, between 2191 and 2840 H-bonds were formed in the ensembles during the 1 µs MD simulation (Table 1). Of them, 73–90% were between the Aβ chains, 10–26% were between Aβ and ligands, and 0.15–4% were between ligands. In the presence of 12 and 36 CU molecules, the average lifetimes of the H-bonds between the Aβ peptides were 41.7 ns and 36.4 ns, respectively, 39.05 ns on average. The corresponding lifetimes for FA anions were 33.9 ns and 38 ns, on average 35.95 ns. The lifetimes of the H-bonds between Aβ peptides and ligands were on average 14.55 ns for CU and 5.7 ns for FA. The H-bonds formed between ligands lived on average 9.8 ns for CU and 2.05 ns for FA. It is clearly evident that CU decreases the number of H-bonds between Aβs without affecting their lifetimes, and it stays bound longer to Aβ chains and to other CU molecules than FA anions do. The lifetimes of H-bonds in the Aβ ensembles are given in detail in Figure S4.

Table 1

Number, percentage, and lifetime of hydrogen bonds formed between Aβ peptides, between Aβ peptides and ligands, and between ligands calculated over 1000 frames (1000 ns).

Hydrogen Bonds	12 Aβ	12 Aβ + 12 CU	12 Aβ + 36 CU	12 Aβ + 12 FA	12 Aβ + 36 FA
Total	2191	2304	2692	2632	2840
Between Aβ peptides	2191 (100%)	2010 (87%)	2041 (76%)	2356 (90%)	2074 (73%)
Average lifetime (ns)	38	41.7	36.4	33.9	38
Between Aβ and ligands	-	277 (12%)	545 (20%)	272 (10%)	738 (26%)
Average lifetime (ns)	-	14.3	14.8	7.2	4.2
Ligand is a donor	-	197 (71%)	386 (71%)	46 (17%)	139 (19%)
Ligand is an acceptor	-	80 (29%)	159 (29%)	226 (83%)	599 (81%)
Between ligands	-	17 (1%)	106 (4%)	4 (0.15%)	28 (1%)
Average lifetime (ns)	-	8.3	11.3	2.5	1.6

In 71% of the H-bonds with Aβ peptides, CU acts as an H-bond donor (Table 1) through the phenolic OH groups and the enol group, which interacts with the side chains of Asp1, Arg5, Asp7, Asp23 and the backbone carbonyl oxygens of Gly9, Lys28, and Leu34. In 29% of the H-bonds, CU is a H-bond acceptor through the keto O-atom such as in the bond with the guanidino group of Arg5 or in the bond with the backbone N-atom of Lys28. In contrast, the FA acts mainly as an H-bond acceptor (82% in average) through the phenolic and carboxy O-atoms.

4. Discussion

In the present study, the primary nucleation of 12 Aβ1-42 peptides was modeled by molecular dynamics in saline at 37 °C for 1 μs. CU is a well-known inhibitor of Aβ self-assembly [58]. At pH 7.4, it is a neutral hydrophobic molecule that is slightly soluble in water (Table 2). The enol form of CU has 11 rotatable bonds, 3 HB donors, and 6 HB acceptors. CU has poor bioavailability (up to 1%) [59] requiring high doses (from 2 to 12 g) to achieve therapeutic concentrations [60]. It crosses the blood–brain barrier (BBB), binds to Aβ peptides, and blocks the aggregation and fibril formation in vitro and in vivo [31]. Even more, CU increases the phagocytosis of Aβ peptides, leading to the effective clearance of plaques from the brain of patients with Alzheimer’s disease [31]. As a diferuloylmethane, CU could be considered as consisting of two FA molecules. FA is the most abundant phenolic acid in whole grain wheat [61,62]. It has strong antioxidant capacity [63,64] and positively affects inflammation, diabetes, cancer, aging, and neurodegeneration [65,66,67]. FA is a weak acid that is freely soluble in water with four rotatable bonds, two HB donors, and four HB acceptors (Table 2). At pH 7.4, the anionic form dominates. As in CU, FA has poor bioavailability [68] but crosses the BBB [69]. Reinke et al. [70] have found that two aromatic rinds are required for the inhibition of Aβ aggregation, and the optimal linker between them should be rigid (up to two rotatable bonds) with a length between 8 and 16 Å. Accordingly, the linker in the enol form of CU consists of seven carbon atoms; five of them are in the sp2 state and two are in sp3, corresponding to two rotatable bonds. The linker length is 12 Å. FA contains only one aromatic ring and is a weaker inhibitor of Aβ aggregation than CU.

Table 2

Physicochemical and ADME (adsorption, distribution, metabolism and excretion) properties of curcumin (CU) and ferulic acid (FA).

Property	Curcumin (CU)	Ferulic Acid (FA)
Molecular weight (g/mol) a	370.38	194.18
Solublity (mg/L) b	3.12	780
LogP a	2.52	1.64
pKa a	9.73; 10.02; 13.3	4.04; 10.22
Polar surface area (Å) a	96.22	66.76
Free rotatable bonds a	11	4
H-bond donors a	3	2
H-bond acceptors a	6	4
IC50 µM (inhibition of Aβ1-42 aggregation)	0.8 c	5.5 d
Bioavailability % e	0.47–1	1
Half-life (h)	12.85 f	1 g
BBB permeability	Yes f	Yes g

a Calculated by ACD/Labs (ACD/ChemSketch, version 2020.1.0); b PubChem [44]; c [31]; d [43]; e DrugBank [71]; f [68]; g [72].

It was found that the IC50 values for the inhibition of Aβ aggregation are 0.8 µM for CU and 5.5 µM for FA (Table 2). In the present study, CU was used as a referent inhibitor of Aβ nucleation, while FA was used as a referent weak inhibitor, and the effects of both compounds on the primary nucleation of Aβ1-42 peptides were studied. As the process of primary nucleation takes several days in vitro and decades in vivo [73], in order to accelerate the in silico simulation, the systems in the present study were modeled to contain molecules in high concentrations: from 9.3 to 31.2 mM. In 10.4 mM solution of Aβ peptides, the nucleation proceeded in 300 ns. A compact nucleus was formed with a micelle-like architecture. A similar, micelle-like architecture with a hydrophobic core and polar surface was observed in the Monte Carlo simulations of Aβ1-40 and Aβ1-42 monomers [74]. Based on NMR analysis, Roche et al. [75] found that the self-association of Aβ1-42 peptides into toxic oligomers is driven by intermolecular interactions between the hydrophobic regions of the peptides instead of elevated propensities of the monomeric species to adopt β-strand-like conformations.

The addition of ligands in ratios 1:1 and 1:3 accelerated the process dose-dependently. However, the cores that were formed were bigger with up to a 20% increase in SASA. In them, the Aβ chains were less flexible and more disordered. In addition, the ligands decreased the number of non-native contacts and hydrogen bonds between the peptides. Interestingly, the presence of ligands in the ensembles did not affect the longevity of the H-bonds between the peptide chains. More than half of them were short-living (less than 10 ns), one third lived between 10 and 100 ns, and only 12% were long-living (more than 100 ns).

CU binds tightly to Aβ peptides by a wide variety of intermolecular interactions: hydrogen bonds, hydrophobic interactions, π–π stacking, and cation–π attraction. Up to 20% of the hydrogen bonds in the ensembles were formed between the CU molecules and the Aβ peptides. Although CU contains more HB acceptors than HB donors, in 71% of the HB, it acted as a donor. The average lifetimes of the CU–Aβ H-bonds were between 14.3 and 14.8 ns. The longest-living HBs (454–506 ns) were with Asp1. CU made π–π interactions with the aromatic residues in Aβ (Phe4, Tyr10, Phe19, and Phe20) and cation–π interactions with the cationic residues (Arg5, Lys16, and Lys28).

FA binds loosely to Aβ peptides. At the end of the 1 µs simulation, between 25% and 50% of the FA molecules were flanking freely around the peptides. The rest were bound to Aβ by short-living H-bonds (lifetimes between 4.2 and 7.2 ns), hydrophobic interactions, π–π stacking with the aromatic residues, and cation–π interactions with the cationic residues.

5. Conclusions

In conclusion, both ligands CU and FA affect the process of primary nucleation of Aβ1-42 peptides, but the effects of CU are more prominent. The CU molecules intercalate among the Aβ chains and bind tightly to them by H-bonds, hydrophobic, π–π, and cation–π interactions. In the presence of CU, the primary nucleus is bigger, while the Aβ peptides are less flexible and more disordered, making fewer non-native contacts and H-bonds. The observed interactions act synergistically and explain at an atomistic level the better inhibitory activity of CU on Aβ aggregation than that of FA. Taken regularly, CU is able to inhibit the process of primary nucleation of Aβ peptides and thus to prevent or at least to delay the onset of neurodegenerative disorders such as Alzheimer’s disease.