Ke Wang
Wensheng Cai- Research Center for Analytical Sciences, College of Chemistry, Tianjin Key Laboratory of Biosensing and Molecular Recognition, Nankai University, Tianjin, China
The destructive effect of Aβ peptides on membranes is an important source of its cytotoxicity in the pathogenesis of Alzheimer’s disease. We have investigated the binding mechanism between the Aβ42 peptide and bilayer in our former work. However, as another abundant form of Aβ peptides in the physiological environment, the binding mechanism between Aβ40 peptide and the lipid bilayer still remains ambiguous. Hence, we performed all-atom simulations on the Aβ40 peptides with the lipid bilayer herein using replica exchange with the solute tempering 2 method. We obtained four major binding models with the hydrophobic C-terminus as the most preferable binding region. Hydrophobic residues and positively charged residues are the principal residues involved in the peptide-bilayer interactions. Aβ40 peptides in our simulation mainly adopt a β-rich conformation in both bound and unbound states. Besides, we determined peptide-water interactions and found that bound peptides prefer forming hydrogen bonds with water molecules than unbound peptides. Our findings herein may provide new insights for the in-depth understanding of the membrane-destructive mechanism of Aβ peptides.
1 Introduction
The formation of senile plaques composed of Aβ peptides is one of the main hallmarks of Alzheimer’s disease (AD) (Guillozet et al., 2003; Chabrier et al., 2012; Bloom, 2014). According to the mainstream view of the “amyloid cascade hypothesis,” the excessive accumulation and abnormal aggregation of Aβ peptides are recognized as the key factors in the onset and progression of AD (Hardy and Selkoe, 2002; Shankar et al., 2007; Sondag et al., 2009; Ferreira and Klein, 2011; Barage and Sonawane, 2015; Awasthi et al., 2016; Minter et al., 2016). Aβ is an intrinsically disordered peptide of 36–43 residues cleaved from amyloid precursor protein, with Aβ40 composed of 40 residues and Aβ42 composed of 42 residues as two predominant forms in physiological conditions (Murphy and LeVine, 2010). Aβ peptides and their products of low aggregation level, such as dimers and trimers, are toxic to the neuron and can cause neuroinflammation and further induce synaptic plasticity impairment and synapse loss (Shankar et al., 2007; Selkoe, 2008; Ferreira and Klein, 2011; Bartolotti et al., 2016; Minter et al., 2016; Mueller-Schiffmann et al., 2016). Abundant evidence shows that many pathogenic features, such as tau hyperphosphorylation, acetylcholine deficiency, oxidative stress, and inflammation, are in the downstream of Aβ pathway (Sondag et al., 2009; Chabrier et al., 2012; Lesne et al., 2013; Bloom, 2014; Amar et al., 2017; Cheignon et al., 2018). Therefore, the investigation of the toxic mechanism of Aβ peptides is vital for the understanding of the pathological mechanism of AD.
Recent studies have demonstrated that the interaction between Aβ peptides and lipid membrane is an important source of its cytotoxicity to AD (Kremer and Murphy, 2003; Arispe et al., 2007; Serra-Batiste et al., 2016; Lindberg et al., 2017; Sparr and Linse, 2019; Ciudad et al., 2020). On the one hand, the membrane can accelerate the aggregation of Aβ peptides by acting as a catalytic site for Aβ nucleation (Kremer and Murphy, 2003; Bokvist et al., 2004; Bokvist and Grobner, 2007; Lindberg et al., 2017; Sparr and Linse, 2019). The aggregation rates of Aβ in a membrane environment are faster than those in bulk solution (Kremer and Murphy, 2003; Bokvist et al., 2004; Bokvist and Grobner, 2007; Banerjee et al., 2020). On the other hand, Aβ peptides can insert into the membrane, forming destructive channels allowing water or ions to flow (Williams and Serpell, 2011; Fantini et al., 2014; Serra-Batiste et al., 2016; Osterlund et al., 2019; Ciudad et al., 2020). Full-length Aβ(1–42) or truncated β-amyloid peptide Aβ(9–42) and Aβ(17–42) peptides can form channels or channel-like structures inside the lipid bilayer (Lin et al., 2001; Arispe et al., 2007; Jang et al., 2010; Serra-Batiste et al., 2016; Bode et al., 2017; Osterlund et al., 2019; Ciudad et al., 2020). The ion-channel-like structures are found to be toxic by inducing neurite degeneration or neuritic abnormality irrespective of their size and morph (Lin et al., 2001; Jang et al., 2010). The channels formed can elicit ion-channel currents, allow calcium uptake, and disrupt the homeostasis of calcium ions (Lin et al., 2001; Quist et al., 2005; Jang et al., 2010; Williams and Serpell, 2011; Fantini et al., 2014). Moreover, the direct interactions between Aβ peptide and bilayer are destructive to the membrane, causing membrane thinning and curvature (Wong et al., 2009; Williams and Serpell, 2011; Garcia-Vinuales et al., 2021). Aβ monomer forms α-helix structure in the membrane-like environment and binding induces a coil-to-helix structure change (Utsumi et al., 2009; Wong et al., 2009). Replica exchange molecular dynamics (REMD) simulations have been utilized to study the Aβ(10–40) monomer in the dimyristoylphosphatidylcholine (DMPC) bilayer environment and found that peptides bound with the bilayer favor the structure with central hydrophobic cluster inserted inside the bilayer (Lockhart et al., 2020). Replica exchange with solute tempering (REST) has also been performed to study the Aβ(25–35) peptide in the DMPC bilayer environment (Smith and Klimov, 2018; Smith et al., 2019; Khayat et al., 2020; Khayat et al., 2021). Interactions between Aβ peptide and bilayer can happen in residues at diverse regions, and both helix, coil, or β-strand structures have been found to exist in peptides at different membrane environments (Brown and Bevan, 2017; Fatafta et al., 2020; Fatafta et al., 2022). Besides, all-atom (AA) simulations and coarse-grained (CG) simulations have also been applied to study the interactions between trimeric or hexameric Aβ fibrils with different bilayers, finding that their binding affinity with bilayer increases with increasing cholesterol content (Agrawal et al., 2023). CG models can be performed in a large timescale due to the simplification of the system, whereas compared with AA models, they lose some key interactions such as hydrogen bonds, salt bridges, etc.
Two Aβ species, Aβ40 and Aβ42, are majorly found in physiological conditions (Murphy and LeVine, 2010; Gu and Guo, 2013; Qiu et al., 2015). The contents of Aβ40 are significantly higher than Aβ42 peptides, whereas Aβ42 peptides are more toxic and more ready to aggregate (Gu and Guo, 2013; Qiu et al., 2015). Previously, we determined the binding mechanism of the Aβ42 peptide with a mixed bilayer using all-atom conventional molecular dynamics (cMD) simulation (Wang et al., 2022). Herein, as a comparison, we carried out all-atom simulations for the full-length Aβ40 peptide with an identical mixed bilayer using the enhanced sampling method of replica exchange with solute tempering 2 (REST2). The sequence for the Aβ40 peptide was divided into four function regions according to the hydrophobicity and charges of residues analogous to our former work of Aβ42: the hydrophilic N-terminus of residues D1-K16 (NT), the central hydrophobic core of residues L17-A21 (CHC), the hydrophilic central loop region of residues E22-G29 (CL), and the hydrophobic C-terminus of residues A30-V40 (CT) (Man Hoang et al., 2014; Cao et al., 2017; Owen et al., 2018; Wang et al., 2022). We found four binding models for the Aβ40 peptide binding to the bilayer with the hydrophobic CT as the most preferable interacting region. We determined the structure features of the Aβ40 peptide of both bound and unbound states. Moreover, we investigated the interactions such as hydrogen bonds formed between peptides and lipids and explored the role of water molecules in peptide-bilayer binding.
2 Methods
2.1 System setup
Before the REST2 simulation, we performed short conventional molecular dynamic simulations for the Aβ40 peptide in solution to acquire ththe onset and progression of ADe pre-equilibrated initial conformations of the peptide for REST2. The initial structure of the full-length Aβ40 peptide is obtained from the Protein Data Bank (PDB ID: 2lfm; Figure 1A) with a helical peptide determined by nuclear magnetic resonance (Vivekanandan et al., 2011). The sequence of the Aβ40 peptide compared with the Aβ42 peptide is shown in Figure 1B. Compared to Aβ42, Aβ40 only lacks the last two hydrophobic residues I41 and A42 at the C-terminal (Tomaselli et al., 2006; Vivekanandan et al., 2011). This peptide was then under the calculation of the protonation state through the H++ web server (Gordon et al., 2005; Anandakrishnan et al., 2012) and put into a rectangle box filled with 0.15 M NaCl and TIP3P (Jorgensen et al., 1983) water molecules. The system underwent a 100 ns cMD at 343 K, and the final conformation was used as the initial peptide structure for the REST2 simulations. The peptide was placed 5.01 nm center of mass (COM) distance above the lipid bilayer and filled with 150 mM NaCl and TIP3P (Jorgensen et al., 1983) water molecules in a rectangular box. The lipid bilayer was constructed in CHARMM-GUI (Jo et al., 2009; Lee et al., 2016) with each leaflet containing 18 cholesterols (CHOL), 18 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-L-serine (POPS) lipids, and 54 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipids, with a total of 36 CHOL, 36 POPS, and 108 POPC in the system. This POPC/POPS/CHOL bilayer in a ratio of POPC: POPS: CHOL = 3: 1: 1 is identical to our previous work (Wang et al., 2022).
FIGURE 1. (A) Initial Aβ40 peptide structure obtained from Protein Data Bank for the pre-equilibrium of the peptide (left) and initial Aβ40 peptide and bilayer structure for the REST2 simulation at 343 K (right). (B) Sequence for the Aβ40 and Aβ42 peptides, where the two additional residues compared to Aβ40 at CT are highlighted in red in Aβ42.
2.2 Simulation details
All the simulations were carried out employing the GROMACS software, version 2020.6 (Van der Spoel et al., 2005), with the CHARMM36m force field (Huang et al., 2017) and a time step of 4 fs with hydrogen-mass repartitioning (Hopkins et al., 2015). The temperature coupling was dealt with using the Nosé–Hoover method while the Parrinello-Rahman barostat (Parrinello and Rahman, 1981; Nos’e and Klein, 1983) was employed for pressure coupling. The periodic boundary condition was used in all directions with the electrostatic interactions using the particle-mesh Ewald method (Essmann et al., 1995). The length of all bonds was constrained by the LINCS algorithm (Hess, 2008). The distance cutoff for the Lennard-Jones potentials and electrostatic interactions were all 1.2 nm. We also performed a 5000-step energy minimization with the steepest descent algorithm followed by 1 ns NVT and 1 ns NPT equilibration before final production runs.
The REST2 method coupled with PLUMED version 2.7.5 (Tribello et al., 2014) was used herein to explore the binding mechanism between full-length Aβ40 peptide and lipid bilayer, which provided an enhanced sampling method analogous to REMD, whereas fewer replicas were needed to achieve equal performance (Wang et al., 2011). The following expression was utilized to determine the temperature distribution in our REST2 (Jo and Jiang, 2015).
Peptides and ions were dealt with as the “hot” region while water molecules and membranes remained cold. Exchanges between the adjacent replicas were attempted every 2 ps. Eighteen replicas were used with effective temperatures ranging between 343 K and 500 K. Each replica runs for 500 ns, resulting in a total of 9 μs simulation times for this entire work.
2.3 Analysis
Trajectories were analyzed using in-house codes and the GROMACS built-in programs. A free energy landscape (FEL) describing the binding process was constructed along the number of contacts and peptide-bilayer distance. Contacts between the Aβ40 peptide and lipid bilayer were determined between any heavy atoms of the peptide and bilayer within 0.5 nm. Peptide-bilayer distance was defined as the perpendicular COM distance between the peptide and the bilayer. For a closer look, the FEL describing each basin in the binding process was also constructed along the root mean square deviation (RMSD) of the peptide with respect to its initial conformation and β-sheet content. Distances between residues and membrane surfaces were generated by computing the perpendicular COM distance between each residue and bilayer and further subtracting half of the membrane thickness (2 nm). Hydrogen bond (H-bond) was defined using the criteria of 0.3 nm donor-acceptor distance cutoff and 20° angle cutoff. Free energies were defined using the following expression (1).
where
3 Results and discussion
3.1 Binding models explored by REST2
The REST2 method was used herein to explore the binding mechanism between Aβ40 and the bilayer. The lipid bilayer is the POPC/POPS/CHOL membrane used in our previous work (Wang et al., 2022). Similarly, this membrane was used here based on its characteristic of being extensively studied and close to the authentic cell membrane, as well as its simplicity (van Meer et al., 2008; Lemkul and Bevan, 2011; Lai et al., 2017; Banerjee et al., 2021). As our previous study revealed, Aβ42 peptides showed little interaction with the bilayer, and the obtained binding models were rare at the physiological temperature of 310 K, while the high temperature of 343 K could make the membranes ready to be penetrated without modifying their basic architecture (Wang et al., 2022). The interactions raised a lot at 343 K and the analyses of Aβ42 peptide-bilayer interactions were all based on the results at 343 K. Hence, with the aim of promoting Aβ40 binding to the lipid bilayer and for the convenience of direct comparison, we performed the REST2 simulation at 343 K for the Aβ40-bilayer system herein. The temperature distribution is calculated using Eq. 1a. Exchange rates between adjacent replicas were all larger than 17% (Supplementary Figure S1). The average exchange rate was 20.80% averaged over 18 replicas. The temperature trajectory for each replica sufficiently visits all the temperatures from 343 K to 500 K (Supplementary Figure S2). All the above outcomes demonstrated the sufficient sampling of our REST2 simulations. Consistent with our expectations, the peptides showed a high tendency to form interactions with bilayers at 343 K (Supplementary Figure S3). For most residues, the contact fractions were larger than 10% and the average contact fraction was 25.29%.
First, the free energy landscape describing the position and interaction of Aβ40 with respect to the bilayer is mapped (Figure 2A). The two-dimensional (2D) FEL is constructed based on the peptide-bilayer distance and the number of peptide-bilayer contacts. Not surprisingly, as the peptide approaches the bilayer, the interactions between the peptide and bilayer increase, leading to the FEL in an L-form. Seven free energy basins, a to g, are identified. Basin c corresponds to the lowest free energy in the 2D FEL, where ΔG = 0. By projecting the 2D FEL onto the peptide-bilayer distance, we obtain the one-dimensional (1D) free energy profile (FEP, Figure 2B). From the FEP, it can be found that Aβ40 is most populated in the region of 2.65 nm (e, f, and g) and 3.33–3.56 nm (a, b, c, and d) distance. These results indicate that Aβ40 energetically favors binding with the bilayer.
FIGURE 2. (A) FEL as a function of the distance and contacts between the Aβ40 peptide and bilayer with the minima are marked. (B) FEL in (A) is projected into the one dimension of distance.
To further identify the bound structures distributed at the seven minima a to g characterized in Figure 2, we extracted the structural ensemble of the peptides at each minimum to construct the 2D FEL of the ensemble for each minimum to describe their structural features (Figure 3A). The FEL is mapped along the positional RMSD and β-sheet content of the peptide. As can be seen from Figure 3A, as the peptides are getting closer to the bilayer (from a to g), the β-sheet content is progressively reduced. Conformations at each minimum can be explicitly classified into two groups, and their representative structures are shown in Figure 3B. One group, labeled as a1, b1, c1, d1, e1, f1, and g1 of Figure 3B, contains bound peptides with predominantly β-sheet structure. Another group, labeled as a2, b2, c2, d2, e2, f2, and g2 of Figure 3B, contains bound peptides with predominantly helix or coil structure.
FIGURE 3. (A) FES for each minimum in Figure 2A as a function of RMSD between heavy atoms of the peptide with respect to its initial conformation and the β-sheet content of the peptide. Microstates in FES of each minimum are labeled at the right. (B) Representative structure of the microstates in (A). Membrane atoms are shown in gray lines. Β-sheet structures are shown in yellow. Helix structures, including α-helix, 3–10 helix, and Pi-helix are shown in red. The remaining turn, bridge, and coil structures are shown in blue.
Snapshots in Figure 3B can accurately represent the microstates in Figure 3A except for a2. a2 microstate contains peptides with few residues interacting with the membrane surface and their conformations are disordered and distinct from each other. The remaining 13 microstates can be classified into four binding models according to their structural features and membrane-interacting regions. Distances between each residue and bilayer surface of the four interacting models are shown in Figure 4A. Snapshots of peptides and bilayers are depicted in Figure 4B. Microstates a1, b1, c1, and d1 are classified into model 1, where barely several residues of CT are inserted inside the bilayer with other residues remaining in the solution. Conformations in this model adopt a β-rich structure of 30%–50% β-sheet content with three parallel β-strands as NT, CHC, and CT region each possessing one β-strand. Microstates b2, c2, and d2 are classified into model 2, where several residues of NT are lying on the membrane surface. Peptides in this model adopt two helical fragments at NT and residues at these two fragments are slightly touching the bilayer surface. Model 3 contains microstates e1, f1, and g1, where the peptides are mainly β-sheet structures with most residues of CT buried inside the bilayer, and several residues in NT and CHC are lying on the membrane surface. Compared to the β-rich structure in model 1, peptides of this model also adopt the structure with three parallel β-strands, whereas the β-sheet structures are of 10%–30% content lower than model 1 with NT, CHC, and CL region, each possessing one β-strand. Besides, residues at CT in model 3 were buried significantly deeper than in model 1. Microstates e2, f2, and g2 are classified into the most deeply buried model 4. Peptides in this model are unstructured and dominated by a random coil with the CHC and CT buried below the membrane surface while NT and CL stretch into the solution. Two to 15 residues of peptides in model 1 and model 2 contact with the bilayer in the 3.0 nm–4.0 nm peptide-bilayer distance region, whereas 16 to 25 residues in model 3 and model 4 contact and insert deeper in the 2.5 nm to the 3.0 nm distance region.
FIGURE 4. Four binding models were acquired from the REST2 simulation. (A) The distance between each residue and membrane surface. Negative values represent residues buried below the bilayer surface. (B) Corresponding structures of peptide and bilayer of each model, where the four function regions NT, CHC, CL, and CT of the peptides are colored in blue, red, yellow, and magenta, respectively.
These binding models of Aβ40 were further compared with those of Aβ42 reported in our previous work (Wang et al., 2022). The results show two significant differences in their binding models. From their representative structures in Figure 4, it can be seen that Aβ40 can adopt more diverse structures of bound peptide including β-strand, helix, or unstructured structure, whereas bound peptides in Aβ42 are mostly unstructured. This difference can be attributed to the highly enhanced sampling efficiency of REST2 compared to cMD as the latter was used in the work of Aβ42. Another difference is that residues at CT frequently interact with the bilayer in Aβ40 (Models one to four in Figure 4) while showing little tendency to bind with the bilayer in Aβ42. This can be explained by the two additional residues Ile41 and Ala42 in CT of Aβ42. Due to the absence of these two residues in Aβ40, two Val residues with stronger hydrophobicity are exposed at the CT region, which can interact with the bilayer more frequently driven by the hydrophobic interactions.
To give a deeper view of the conformation for the binding peptides, the H-bonds formed between the peptide and membrane/water for minima a to g were obtained (Figure 5A). From minima a to g, with the contacts increased, the number of H-bonds formed by Aβ40 with bilayer also shows a roughly rising trend. However, there is no obvious rising or declining trend for the H-bond formed with water. We then computed the H-bonds formed between binding models and the membrane/water (Figure 5B). The peptide-bilayer H-bonds exhibit no correlation with the peptide-water H-bonds. Conformations in model 2 form the highest number of H-bonds with the membrane but also show a high tendency to form H-bonds with water. Conformations in model 1 form the least H-bonds with the membrane. The average number of H-bonds formed by each residue with membrane was then acquired to explore the role of residues in peptide-membrane interactions in four binding models (Figure 5C). Two residues forming the highest number of H-bond with bilayer were pointed in each model. Most of these residues were charged especially for the deepest binding model of model 4, and positively charged K16 and K28 form strong H-bonds with the bilayer. Positively charged R5 residue at NT also plays a key role in forming H-bonds in models 2 and 3. We have labeled these residues intuitively in the structure of binding models in Figure 5D. They are all distributed in the peptide-bilayer interacting regions. The free energy data were computed using Eq. 1b to compare the relative stability of each microstate (Supplementary Table S1). Microstates d1, g2, and c1 showed the lowest free energy, indicating the high stability of these microstates. Free energies of other microstates were relatively high, especially for b2, c2, and d2 of model 2, suggesting the instability of this binding model.
FIGURE 5. (A) The average number of H-bonds formed between peptide and membrane/water for minima a to g. (B) The average number of H-bonds formed between peptide and membrane/water for four binding models. (C) The average number of H-bonds formed between each residue and membrane for four binding models. The two residues forming the highest number of H-bonds are pointed in each model. (D) The representative conformations for four binding models with the two residues forming the highest H-bond are highlighted in green balls. The color scheme for peptide and bilayer is identical to Figure 3.
3.2 Structural features of the Aβ40 peptide
Aβ peptides, as numerous works previously characterized, are intrinsically disordered in bulk water (Viola and Klein, 2015; Ono, 2018). The initial Aβ40 conformation used here is a partially folded structure with a helix structure spanning residues H13-D23. To explore the conformational transition of Aβ40 induced by binding to the lipid bilayer, an analysis of the secondary structure was carried out, as delineated in Figure 6. Aβ40 peptides were considered bound when their peptide-bilayer contacts in Figure 2A are larger than zero; otherwise, they are considered unbound. As can be seen from Figures 6A, B, for both the bound and unbound peptides, β-sheet (E), turn (T), and coil (C) structures were dominant. These three structures took up at least 79.50% and 90.80% of the bound and unbound peptides, respectively, and alternatively exist in each residue. In Figure 6C, the seven structures obtained were classified into three groups: the sum of the fraction for coil, turn, and bridge structure (C + T + B) was taken as the unstructured structure fraction; the fraction for β-sheet structure (E) was taken as the β-sheet structure fraction; the sum of 3–10 helix, α-helix, and Pi-helix structure (G + H + I) was taken as the helix structure fraction. Fractions of helix structure were low in both bound and unbound peptides, whereas residues E3-R5 and H13-Q15 at NT could form helix structures with fractions larger than 10% in the bound peptides. Peptides possessed more helix structure in the bound state than in the unbound state. β-sheet content (E) of both bound and unbound peptides was much higher than that of helix, especially for unbound peptides. Residues D7, E11, Q15, K16, L17, and F20 in unbound peptides adopted the β-sheet structure with fractions larger than 80%, which were all distributed at the NT and CHC regions. In the bound state, the highest fraction of β-sheet was 49.84%, which was evidently lower than the unbound state, and the residues with β-sheet structure were widely distributed in NT, CHC, CL, and CT regions. In summary, the β-sheet structure is dominant in both bound and unbound Aβ40 peptides.
FIGURE 6. (A) Secondary structure probability for each residue of bound peptide. (B) Secondary structure probability for each residue of unbound peptide. (C) Secondary structure probability for each residue for the sum of coil, turn, bridge structure (C + T + B); the sum of the β-sheet (E); and the sum of 3–10 helix, α-helix, and Pi-helix structure (G + H + I) for bound and unbound Aβ40 peptides. E: extended β-sheet structure, B: bridge structure, H: α-helix structure, G: 3–10 helix structure, I: Pi-helix structure, T: turn structure, and C: coil structure.
3.3 Peptide-bilayer interactions
To explore the driving force for the peptide binding, the peptide-bilayer contacts and H-bonds were analyzed. It is clear from Figure 7A; Supplementary Figure S4 that binding events occurred mainly in the NT, CHC, and CT regions. The contact fractions of residues in hydrophilic CL were significantly lower, especially for residues E22-N27. Residues at CT showed the highest tendency to interact with the bilayer. The top 10 residues with the highest contact fractions were G38, A40, G37, V39, Q15, V36, K16, H14, M35, and L34, which were all hydrophobic or positively charged residues. We compared the average fractions of residue-bilayer contact between Aβ40 and Aβ42 in the four sequence regions in Supplementary Figure S4 (Wang et al., 2022). It is remarkable that residues at CT in Aβ40 frequently interacted with the bilayer, whereas they showed little tendency to bind with the bilayer in Aβ42. From Figure 7B, R5 formed the highest number of H-bonds much larger than other residues, followed by K16 and K28 residues. These three residues also exhibited the same behavior in Aβ42 (Wang et al., 2022). The highest 10 types of H-bonds formed between R5, K16, and K28 residues and bilayer are provided in Supplementary Figure S5. The most populated H-bonds are formed with the side chain of POPC or POPS lipids. Like in the Aβ42-bilayer system, cholesterols also do not favor forming H-bonds with these residues (Wang et al., 2022). Overall, residues in the CT region of the Aβ40 peptide showed the highest tendency to interact with the bilayer, and H-bonds formed between positively charged residues and bilayer may drive the binding. Both hydrophobic interactions and electrostatic interactions contributed to the binding of Aβ40 to the bilayer.
FIGURE 7. (A) Contact fractions of residues with the membrane in the bound peptides. When there is a contact between any of the heavy atoms of one residue and the bilayer, then this residue is considered to be in contact. (B) The average number of H-bonds formed between each residue and the bilayer.
3.4 Intra-peptide and peptide-water interactions
Besides peptide-bilayer interactions, intra-peptide and peptide-water interactions may also drive the peptide bind to the bilayer. Residue-residue contact maps for the bound and unbound Aβ40 peptides were plotted to investigate the intra-peptide interactions (Figure 8A). Unlike Aβ42, the conformations of Aβ40 in the unbound state were rather different from the bound state (Wang et al., 2022). The evident cross-diagonal formed by residues F4-V12 and residues H14-D23 in the contact map of the unbound state corresponded to the anti-parallel β-sheet structure formed at these regions. The contact probabilities in these regions were extremely high, corresponding to the high β-sheet structure of these regions in Figure 6C. A cross-diagonal represented that the β-sheet structure was also formed by residues H14-A21 and V24-V36 of the bound peptide, whereas the probabilities were not as high as the β-sheet structure in the unbound peptide. Besides, there was a band indicating the contacts formed between residues A2-H6 and residues G29-G33, which corresponded to the β-sheet interactions at these two fragments. The highest 10 contact fractions for the residue pairs (at least three residues apart) were used to characterize long-range interactions (Supplementary Figure S6).
FIGURE 8. (A) Residue-residue contact maps of peptides in bound and unbound states. The color bar corresponding to the contact probability from 0 to 1 is shown on the right. (B) The average number of H-bonds formed between each residue and water molecules for bound and unbound peptides.
To explore the peptide-water interactions, the average number of H-bonds formed between peptides and water molecules was counted (Figure 8B). Surprisingly, for most residues in the bound state, the average number of H-bonds formed with water was higher than residues in the unbound state except for a few residues at CT. This phenomenon occurred due to the fact that bound peptides were more exposed to water than unbound peptides and water molecules could drive the binding of peptide to the bilayer. Another reasonable explanation is that Aβ40 in extended solvent-accessible structure preferred binding with the bilayer than in the globular structure. These results suggested that water molecules played an important role in the binding of Aβ40 to membranes.
4 Conclusion
In the present work, we investigated the interactions between full-length Aβ40 peptide and POPC/POPS/CHOL bilayer using the REST2 method. We first explored the binding mechanism by extracting the conformations in the free energy landscape. The conformations could be classified into four binding models. Peptides in model 1 adopted a β-rich structure with only several residues of CT inserted inside the bilayer. Model 2 contained peptides with NT lying on the membrane surface and adopted two helical fragments at NT. Peptides in model 3 were also β-rich structures with most residues of CT buried inside the bilayer and several residues in NT and CHC lying on the bilayer surface. The most deeply buried model 4 contains peptides dominated by the random coil with the CHC and CT buried deeply inside the membrane. Hydrophobic CT was the region showing the highest tendency to interact with the bilayer. Residues most preferably forming H-bonds with bilayer were positively charged R5, K16, and K28 residues. Aβ40 peptides in both bound and unbound states mainly adopted the β-rich structure, whereas bound peptides showed slightly higher fractions of the helix structure than unbound peptides. We also computed the H-bonds formed between peptide and water molecules to unveil the role of water in peptide-bilayer binding. Peptides in the bound state form more H-bonds with water than in the unbound state, which showed strong proof of the vital role of water molecules in driving the peptide-membrane binding.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
KW: Writing–original draft, Writing–review and editing, Data curation, Formal Analysis, Methodology, Software, Visualization. WC: Writing–original draft, Writing–review and editing, Funding acquisition, Project administration, Supervision.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. The work was supported by the National Natural Science Foundation of China (22073050).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2024.1367793/full#supplementary-material
References
Agrawal, N., Skelton, A. A., and Parisini, E. (2023). A coarse-grained molecular dynamics investigation on spontaneous binding of Aβ1–40 fibrils with cholesterol-mixed DPPC bilayers. Comp. Struct. Biotechnol. J. 21, 2688–2695. doi:10.1016/j.csbj.2023.04.013
Amar, F., Sherman, M. A., Rush, T., Larson, M., Boyle, G., Chang, L., et al. (2017). The amyloid-β oligomer Aβ*56 induces specific alterations in neuronal signaling that lead to tau phosphorylation and aggregation. Sci. Signal. 10, eaal2021. doi:10.1126/scisignal.aal2021
Anandakrishnan, R., Aguilar, B., and Onufriev, A. V. (2012). H++3.0: automating pK prediction and the preparation of biomolecular structures for atomistic molecular modeling and simulations. Nucleic Acids Res. 40, W537–W541. doi:10.1093/nar/gks375
Arispe, N., Diaz, J. C., and Simakova, O. (2007). Aβ ion channels. Prospects for treating Alzheimer's disease with Aβ channel blockers. Biochim. Biophys. Acta-Biomembr. 1768, 1952–1965. doi:10.1016/j.bbamem.2007.03.014
Awasthi, M., Singh, S., Pandey, V. P., and Dwivedi, U. N. (2016). Alzheimer's disease: an overview of amyloid beta dependent pathogenesis and its therapeutic implications along with in silico approaches emphasizing the role of natural products. J. Neurol. Sci. 361, 256–271. doi:10.1016/j.jns.2016.01.008
Banerjee, S., Hashemi, M., Zagorski, K., and Lyubchenko, Y. L. (2020). Interaction of Aβ42 with membranes triggers the self-assembly into oligomers. Int. J. Mol. Sci. 21, 1129. doi:10.3390/ijms21031129
Banerjee, S., Hashemi, M., Zagorski, K., and Lyubchenko, Y. L. (2021). Cholesterol in membranes facilitates aggregation of amyloid β protein at physiologically relevant concentrations. ACS Chem. Neurosci. 12, 506–516. doi:10.1021/acschemneuro.0c00688
Barage, S. H., and Sonawane, K. D. (2015). Amyloid cascade hypothesis: pathogenesis and therapeutic strategies in Alzheimer's disease. Neuropeptides 52, 1–18. doi:10.1016/j.npep.2015.06.008
Bartolotti, N., Bennett, D. A., and Lazarov, O. (2016). Reduced pCREB in Alzheimer's disease prefrontal cortex is reflected in peripheral blood mononuclear cells. Mol. Psychiatr. 21, 1158–1166. doi:10.1038/mp.2016.111
Bloom, G. S. (2014). Amyloid-β and tau: the trigger and bullet in alzheimer disease pathogenesis. JAMA Neurol. 71, 505–508. doi:10.1001/jamaneurol.2013.5847
Bode, D. C., Baker, M. D., and Viles, J. H. (2017). Ion Channel formation by amyloid-β42 oligomers but not amyloid-β40 in cellular membranes. J. Biol. Chem. 292, 1404–1413. doi:10.1074/jbc.M116.762526
Bokvist, M., and Grobner, G. (2007). Misfolding of amyloidogenic proteins at membrane surfaces: the impact of macromolecular crowding. J. Am. Chem. Soc. 129, 14848–14849. doi:10.1021/ja076059o
Bokvist, M., Lindstrom, F., Watts, A., and Grobner, G. (2004). Two types of alzheimer's β-amyloid (1–40) peptide membrane interactions: aggregation preventing transmembrane anchoring versus accelerated surface fibril formation. J. Mol. Biol. 335, 1039–1049. doi:10.1016/j.jmb.2003.11.046
Brown, A. M., and Bevan, D. R. (2017). Influence of sequence and lipid type on membrane perturbation by human and rat amyloid β-peptide (1–42). Arch. Biochem. Biophys. 614, 1–13. doi:10.1016/j.abb.2016.11.006
Cao, Y., Jiang, X., and Han, W. (2017). Self-assembly pathways of β-sheet-rich amyloid-β(1–40) dimers: Markov state model analysis on millisecond hybrid-resolution simulations. J. Chem. Theory Comput. 13, 5731–5744. doi:10.1021/acs.jctc.7b00803
Chabrier, M. A., Blurton-Jones, M., Agazaryan, A. A., Nerhus, J. L., Martinez-Coria, H., and LaFerla, F. M. (2012). Soluble Aβ promotes wild-type tau pathologyin vivo. J. Neurosci. 32, 17345–17350. doi:10.1523/jneurosci.0172-12.2012
Cheignon, C., Tomas, M., Bonnefont-Rousselot, D., Faller, P., Hureau, C., and Collin, F. (2018). Oxidative stress and the amyloid beta peptide in Alzheimer’s disease. Redox Biol. 14, 450–464. doi:10.1016/j.redox.2017.10.014
Ciudad, S., Puig, E., Botzanowski, T., Meigooni, M., Arango, A. S., Do, J., et al. (2020). Aβ(1-42) tetramer and octamer structures reveal edge conductivity pores as a mechanism for membrane damage. Nat. Commun. 11, 3014. doi:10.1038/s41467-020-16566-1
Essmann, U., Perera, L., Berkowitz, M. L., Darden, T., Lee, H., and Pedersen, L. G. (1995). A smooth particle mesh ewald method. J. Chem. Phys. 103, 8577–8593. doi:10.1063/1.470117
Fantini, J., Di Scala, C., Yahi, N., Troadec, J.-D., Sadelli, K., Chahinian, H., et al. (2014). Bexarotene blocks calcium-permeable ion channels formed by neurotoxic Alzheimer’s β-amyloid peptides. ACS Chem. Neurosci. 5, 216–224. doi:10.1021/cn400183w
Fatafta, H., Kav, B., Bundschuh, B. F., Loschwitz, J., and Strodel, B. (2022). Disorder-to-order transition of the amyloid-β peptide upon lipid binding. Biophys. Chem. 280, 106700. doi:10.1016/j.bpc.2021.106700
Fatafta, H., Poojari, C., Sayyed-Ahmad, A., Strodel, B., and Owen, M. C. (2020). Role of oxidized Gly25, Gly29, and Gly33 residues on the interactions of aβ1–42 with lipid membranes. ACS Chem. Neurosci. 11, 535–548. doi:10.1021/acschemneuro.9b00558
Ferreira, S. T., and Klein, W. L. (2011). The Aβ oligomer hypothesis for synapse failure and memory loss in Alzheimer’s disease. Neurobiol. Learn. Mem. 96, 529–543. doi:10.1016/j.nlm.2011.08.003
Garcia-Vinuales, S., Sciacca, M. F. M., Lanza, V., Santoro, A. M., Grasso, G., Tundo, G. R., et al. (2021). The interplay between lipid and Aβ amyloid homeostasis in Alzheimer’s Disease: risk factors and therapeutic opportunities. Chem. Phys. Lipids 236, 105072. doi:10.1016/j.chemphyslip.2021.105072
Gordon, J. C., Myers, J. B., Folta, T., Shoja, V., Heath, L. S., and Onufriev, A. (2005). H++: a server for estimating pK(a)s and adding missing hydrogens to macromolecules. Nucleic Acids Res. 33, W368–W371. doi:10.1093/nar/gki464
Gu, L., and Guo, Z. (2013). Alzheimer's Aβ42 and Aβ40 peptides form interlaced amyloid fibrils. J. Neurochem. 126, 305–311. doi:10.1111/jnc.12202
Guillozet, A. L., Weintraub, S., Mash, D. C., and Mesulam, M. M. (2003). Neurofibrillary tangles, amyloid, and memory in aging and mild cognitive impairment. Arch. Neurol. 60, 729–736. doi:10.1001/archneur.60.5.729
Hardy, J., and Selkoe, D. J. (2002). The amyloid hypothesis of alzheimer's disease: progress and problems on the road to therapeutics. Science 297, 353–356. doi:10.1126/science.1072994
Hess, B. (2008). P-LINCS: a parallel linear constraint solver for molecular simulation. J. Chem. Theory Comput. 4, 116–122. doi:10.1021/ct700200b
Hopkins, C. W., Le Grand, S., Walker, R. C., and Roitberg, A. E. (2015). Long-time-step molecular dynamics through hydrogen mass repartitioning. J. Chem. Theory Comput. 11, 1864–1874. doi:10.1021/ct5010406
Huang, J., Rauscher, S., Nawrocki, G., Ran, T., Feig, M., de Groot, B. L., et al. (2017). CHARMM36m: an improved force field for folded and intrinsically disordered proteins. Nat. Methods. 14, 71–73. doi:10.1038/nmeth.4067
Humphrey, W., Dalke, A., and Schulten, K. (1996). VMD: visual molecular dynamics. J. Mol. Graph. 14 (33-38), 33–38. doi:10.1016/0263-7855(96)00018-5
Jang, H., Arce, F. T., Ramachandran, S., Capone, R., Azimova, R., Kagan, B. L., et al. (2010). Truncated β-amyloid peptide channels provide an alternative mechanism for Alzheimer’s Disease and Down syndrome. Proc. Natl. Acad. Sci. U. S. A. 107, 6538–6543. doi:10.1073/pnas.0914251107
Jo, S., and Jiang, W. (2015). A generic implementation of replica exchange with solute tempering (REST2) algorithm in NAMD for complex biophysical simulations. Comput. Phys. Commun. 197, 304–311. doi:10.1016/j.cpc.2015.08.030
Jo, S., Lim, J. B., Klauda, J. B., and Im, W. (2009). CHARMM-GUI membrane builder for mixed bilayers and its application to yeast membranes. Biophys. J. 97, 50–58. doi:10.1016/j.bpj.2009.04.013
Jorgensen, W. L., Chandrasekhar, J., Madura, J. D., Impey, R. W., and Klein, M. L. (1983). Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 79, 926–935. doi:10.1063/1.445869
Khayat, E., Klimov, D. K., and Smith, A. K. (2020). Phosphorylation promotes aβ25–35 peptide aggregation within the DMPC bilayer. ACS Chem. Neurosci. 11, 3430–3441. doi:10.1021/acschemneuro.0c00541
Khayat, E., Lockhart, C., Delfing, B. M., Smith, A. K., and Klimov, D. K. (2021). Met35 oxidation hinders aβ25-35 peptide aggregation within the dimyristoylphosphatidylcholine bilayer. ACS Chem. Neurosci. 12, 3225–3236. doi:10.1021/acschemneuro.1c00407
Kremer, J. J., and Murphy, R. M. (2003). Kinetics of adsorption of β-amyloid peptide Aβ(1–40) to lipid bilayers. J. Biochem. Biophys. Methods. 57, 159–169. doi:10.1016/s0165-022x(03)00103-9
Lai, A. L., Millet, J. K., Daniel, S., Freed, J. H., and Whittaker, G. R. (2017). The SARS-CoV fusion peptide forms an extended bipartite fusion platform that perturbs membrane order in a calcium-dependent manner. J. Mol. Biol. 429, 3875–3892. doi:10.1016/j.jmb.2017.10.017
Lee, J., Cheng, X., Swails, J. M., Yeom, M. S., Eastman, P. K., Lemkul, J. A., et al. (2016). CHARMM-GUI input generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM simulations using the CHARMM36 additive force field. J. Chem. Theory Comput. 12, 405–413. doi:10.1021/acs.jctc.5b00935
Lemkul, J. A., and Bevan, D. R. (2011). Lipid composition influences the release of Alzheimer's amyloid β-peptide from membranes. Protein Sci. 20, 1530–1545. doi:10.1002/pro.678
Lesne, S. E., Sherman, M. A., Grant, M., Kuskowski, M., Schneider, J. A., Bennett, D. A., et al. (2013). Brain amyloid-β oligomers in ageing and Alzheimer’s disease. Brain 136, 1383–1398. doi:10.1093/brain/awt062
Lin, H., Bhatia, R., and Lal, R. (2001). Amyloid β protein forms ion channels: implications for Alzheimer's disease pathophysiology. Faseb J. 15, 2433–2444. doi:10.1096/fj.01-0377com
Lindberg, D. J., Wesen, E., Bjorkeroth, J., Rocha, S., and Esbjorner, E. K. (2017). Lipid membranes catalyse the fibril formation of the amyloid-β (1–42) peptide through lipid-fibril interactions that reinforce secondary pathways. Biochim. Biophys. Acta-Biomembr. 1859, 1921–1929. doi:10.1016/j.bbamem.2017.05.012
Lockhart, C., Smith, A. K., and Klimov, D. K. (2020). Three popular force fields predict consensus mechanism of amyloid β peptide binding to the dimyristoylgylcerophosphocholine bilayer. J. Chem. Inf. Model. 60, 2282–2293. doi:10.1021/acs.jcim.0c00096
Man Hoang, V., Nguyen, P. H., Derreumaux, P., and Li, M. S. (2014). Effect of the English familial disease mutation (H6R) on the monomers and dimers of Aβ40 and Aβ42. ACS Chem. Neurosci. 5, 646–657. doi:10.1021/cn500007j
Minter, M. R., Taylor, J. M., and Crack, P. J. (2016). The contribution of neuroinflammation to amyloid toxicity in Alzheimer's disease. J. Neurochem. 136, 457–474. doi:10.1111/jnc.13411
Mueller-Schiffmann, A., Herring, A., Abdel-Hafiz, L., Chepkova, A. N., Schaeble, S., Wedel, D., et al. (2016). Amyloid-β dimers in the absence of plaque pathology impair learning and synaptic plasticity. Brain 139, 509–525. doi:10.1093/brain/awv355
Murphy, M. P., and LeVine, H. (2010). Alzheimer's disease and the amyloid-β peptide. J. Alzheimers Dis. 19, 311–323. doi:10.3233/jad-2010-1221
Nos’e, S., and Klein, M. L. (1983). Constant pressure molecular dynamics for molecular systems. Mol. Phys. 50, 1055–1076. doi:10.1080/00268978300102851
Ono, K. (2018). Alzheimer's disease as oligomeropathy. Neurochem. Int. 119, 57–70. doi:10.1016/j.neuint.2017.08.010
Osterlund, N., Moons, R., Ilag, L. L., Sobott, F., and Graslund, A. (2019). Native ion mobility-mass spectrometry reveals the formation of β-barrel shaped amyloid-β hexamers in a membrane-mimicking environment. J. Am. Chem. Soc. 141, 10440–10450. doi:10.1021/jacs.9b04596
Owen, M. C., Kulig, W., Poojari, C., Rog, T., and Strodel, B. (2018). Physiologically-relevant levels of sphingomyelin, but not GM1, induces a β-sheet-rich structure in the amyloid-β(1-42) monomer. Biochim. Biophys. Acta-Biomembr. 1860, 1709–1720. doi:10.1016/j.bbamem.2018.03.026
Parrinello, M., and Rahman, A. (1981). Polymorphic transitions in single crystals: a new molecular dynamics method. J. Appl. Phys. 52, 7182–7190. doi:10.1063/1.328693
Qiu, T., Liu, Q., Chen, Y.-X., Zhao, Y.-F., and Li, Y.-M. (2015). Aβ42 and Aβ40: similarities and differences. J. Pept. Sci. 21, 522–529. doi:10.1002/psc.2789
Quist, A., Doudevski, L., Lin, H., Azimova, R., Ng, D., Frangione, B., et al. (2005). Amyloid ion channels: a common structural link for protein-misfolding disease. Proc. Natl. Acad. Sci. U. S. A. 102, 10427–10432. doi:10.1073/pnas.0502066102
Selkoe, D. J. (2008). Soluble oligomers of the amyloid β-protein impair synaptic plasticity and behavior. Behav. Brain Res. 192, 106–113. doi:10.1016/j.bbr.2008.02.016
Serra-Batiste, M., Ninot-Pedrosa, M., Bayoumi, M., Gairi, M., Maglia, G., and Carulla, N. (2016). Aβ42 assembles into specific β-barrel pore-forming oligomers in membrane-mimicking environments. Proc. Natl. Acad. Sci. U. S. A. 113, 10866–10871. doi:10.1073/pnas.1605104113
Shankar, G. M., Bloodgood, B. L., Townsend, M., Walsh, D. M., Selkoe, D. J., and Sabatini, B. L. (2007). Natural oligomers of the alzheimer amyloid-β protein induce reversible synapse loss by modulating an NMDA-type glutamate receptor-dependent signaling pathway. J. Neurosci. 27, 2866–2875. doi:10.1523/jneurosci.4970-06.2007
Smith, A. K., Khayat, E., Lockhart, C., and Klimov, D. K. (2019). Do cholesterol and sphingomyelin change the mechanism of Aβ(25-35) peptide binding to zwitterionic bilayer? J. Chem. Inf. Model. 59, 5207–5217. doi:10.1021/acs.jcim.9b00763
Smith, A. K., and Klimov, D. K. (2018). Binding of cytotoxic aβ25–35 peptide to the dimyristoylphosphatidylcholine lipid bilayer. J. Chem. Inf. Model. 58, 1053–1065. doi:10.1021/acs.jcim.8b00045
Sondag, C. M., Dhawan, G., and Combs, C. K. (2009). Beta amyloid oligomers and fibrils stimulate differential activation of primary microglia. J. Neuroinflamm. 6, 1. doi:10.1186/1742-2094-6-1
Sparr, E., and Linse, S. (2019). Lipid-protein interactions in amyloid formation. BBA-Proteins Proteomics 1867, 455–457. doi:10.1016/j.bbapap.2019.03.006
Tomaselli, S., Esposito, V., Vangone, P., van Nuland, N. A. J., Bonvin, A., Guerrini, R., et al. (2006). The α-to-β conformational transition of Alzheimer's Aβ-(1-42) peptide in aqueous media is reversible: a step by step conformational analysis suggests the location of β conformation seeding. Chembiochem 7, 257–267. doi:10.1002/cbic.200500223
Tribello, G. A., Bonomi, M., Branduardi, D., Camilloni, C., and Bussi, G. (2014). PLUMED 2: new feathers for an old bird. Comput. Phys. Commun. 185, 604–613. doi:10.1016/j.cpc.2013.09.018
Utsumi, M., Yamaguchi, Y., Sasakawa, H., Yamamoto, N., Yanagisawa, K., and Kato, K. (2009). Up-and-down topological mode of amyloid β-peptide lying on hydrophilic/hydrophobic interface of ganglioside clusters. Glycoconj. J. 26, 999–1006. doi:10.1007/s10719-008-9216-7
Van der Spoel, D., Lindahl, E., Hess, B., Groenhof, G., Mark, A. E., and Berendsen, H. J. C. (2005). GROMACS: fast, flexible, and free. J. Comput. Chem. 26, 1701–1718. doi:10.1002/jcc.20291
van Meer, G., Voelker, D. R., and Feigenson, G. W. (2008). Membrane lipids: where they are and how they behave. Nat. Rev. Mol. Cell Biol. 9, 112–124. doi:10.1038/nrm2330
Viola, K. L., and Klein, W. L. (2015). Amyloid β oligomers in Alzheimer’s disease pathogenesis, treatment, and diagnosis. Acta Neuropathol. 129, 183–206. doi:10.1007/s00401-015-1386-3
Vivekanandan, S., Brender, J. R., Lee, S. Y., and Ramamoorthy, A. (2011). A partially folded structure of amyloid-beta(1–40) in an aqueous environment. Biochem. Biophys. Res. Commun. 411, 312–316. doi:10.1016/j.bbrc.2011.06.133
Wang, K., Shao, X., and Cai, W. (2022). Binding models of Aβ42 peptide with membranes explored by molecular simulations. J. Chem. Inf. Model. 62, 6482–6493. doi:10.1021/acs.jcim.2c00444
Wang, L., Friesner, R. A., and Berne, B. J. (2011). Replica exchange with solute scaling: a more efficient version of replica exchange with solute tempering (REST2). J. Phys. Chem. B 115, 9431–9438. doi:10.1021/jp204407d
Williams, T. L., and Serpell, L. C. (2011). Membrane and surface interactions of Alzheimer’s Aβ peptide – insights into the mechanism of cytotoxicity. Febs J. 278, 3905–3917. doi:10.1111/j.1742-4658.2011.08228.x
Keywords: Aβ40 peptide, lipid bilayer, binding mechanism, binding models, peptide-bilayer interactions
Citation: Wang K and Cai W (2024) Binding mechanism of full-length Aβ40 peptide to a mixed lipid bilayer. Front. Chem. 12:1367793. doi: 10.3389/fchem.2024.1367793
Received: 09 January 2024; Accepted: 05 February 2024;
Published: 21 February 2024.
Edited by:
Angelo Spinello, University of Palermo, ItalyReviewed by:
Daniel Mejía Rodríguez, Pacific Northwest National Laboratory (DOE), United StatesLixue Cheng, Microsoft Research, Germany
Copyright © 2024 Wang and Cai. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wensheng Cai, wscai@nankai.edu.cn