1. Introduction
Alzheimer’s (AD) and Parkinson’s (PD) diseases are neurodegenerative disorders that are associated with the formation of insoluble protein deposits. These deposits disrupt the normal functioning of neurons, blocking neurotransmitters and leading to apoptosis. The main protein responsible for the damage in Alzheimer’s disease is beta-amyloid (A
), which is formed from an amyloidogenic process, where two enzymes,
-secretase (BACE1) and
-secretase, cleave the amyloid precursor protein (APP) into A
monomers [
1,
2]. These parts can aggregate to form fibrils, and then plaques. The main protein responsible for the damage in Parkinson’s disease is alpha-synuclein (
-syn), which is involved in the process of storing dopamine neurotransmitter reserves. The cause of its aggregation is still unknown. Scientific evidence points to the centrality of
-syn, an intrinsically disordered protein in the cytosol that is considered to be essential for the development of the disease. The etiology of PD is varied, and depending on the triggering factor of the pathology, several changes can occur in the physiological pathways and the monomer
-syn. Changes in the protein can lead to damage to cellular pathways, and damage to these pathways can also cause changes in the protein. Ultimately, neuronal death occurs, and in resilient neurons, Lewy bodies composed mainly of
-syn are found. These bodies are formed because, in an initial phase, changes in cellular homeostasis induce disturbances in the equilibrium between the different forms of
-syn, originating in the formation of toxic oligomers and fibrils. In addition to the damage they cause to cellular pathways, these species can spread in a prion-like manner, interfering with the homeostasis of neighboring cells. It is believed that their accumulation is due to combinations of genetic (mutations in Park1) and environmental factors, which can reduce the ability of this protein to cleave
-syn [
3,
4,
5].
The most common A filaments found in AD are peptides of 38-43 amino acids, with the A(1-42) and A(1-40) forms being the most common. Aggregation depends primarily on the hydrophobic amino acid portion, which, due to this property, makes it prone to clumping together, forming a hydrophobic structure. This portion represents the amino acids from 16-20, consisting of Lysine 16 (K), Leucine 17 (L), Valine 18 (V), Phenylalanine 19 (F), and Phenylalanine 20 (F), and from 39-42, consisting of Valine 39 (V), Valine 40 (V), Isoleucine 41 (Ile), and Alanine (42). The latter is known as the hydrophobic C-terminus. The most common -syn filaments found in PD are peptides of 121 amino acids. The main amino acids, also hydrophobic, responsible for the stability of the -syn protofilaments are those from 50-57, found in the interface region between two -syn fibrils. These are Histidine (50), Glycine (51), Valine (52), Alanine (53), Threonine (54), Valine (55), Alanine (56), and Glutamate (57).
Finding drugs that promote the disaggregation of A
and
-syn is as important as obtaining drugs that can act as antiaggregators, since oligomers, which are smaller aggregates of these proteins, are even more toxic than fibrils [
6]. Oligomers can penetrate cells and cause more damage than fibrils [
7,
8]. Therefore, it is important to develop therapies that can disaggregate fibrils and prevent them from re-aggregating into oligomers. Based on this knowledge, some research groups have developed some experimental and computational strategies to find inhibitors that can block the aggregation process, intervene in metal complexation, or even induce the disaggregation of A
and
-syn fibrils [
9,
10,
11,
12].
Following these same steps, our research group has been investigating the properties of some peptides derived from wasp venom intending to find a compound that can inhibit the action of the peptides A
and
-syn. Among them, we chose to investigate the derivatives of the venom of the wasps
Polybia occidentalis and
Parachartergus fraternus. From the first, it was obtained the compound called Fraternine-10 [
13,
14] and, from the second, Octovespin [
15]. The molecular structure of Octovespin and Fraternine-10 is presented in
Figure 1, Fraternine-10 has the peptide sequence: Leu-Ser-Phe-Gln-Lys-Val-Lys-Glu-Lys-Val-NH2, and Octovespin: Gln-Met-Trp-Ala-Val-Phe-Trp-Met-NH2. Fraternine-10 showed a relative antiparkinsonian effect in a model of neurodegeneration in mice [
16].
In vitro assays showed that Octovespin decreased the aggregation of A
fibrils, and
in vivo assays in mice demonstrated the ability of this compound to prevent the aggregation of A
and reduce induced negative cognitive effects [
17]. No information on Fraternine-10 efficacy as an anti-aggregator to A
and
-syn is available at the moment.
Based on these reports, the present work aims to investigate, using computational methods, the properties and capabilities of the peptides Octovespin and Fraternine-10 to promote the disaggregation of the fibrils of A and -syn. To this end, molecular docking and molecular dynamics techniques were used using A(18-42) and -syn(38-95) fibrils as targets. Additionally, experimental results of Fraternine-10 are presented regarding its ability to prevent or delay the formation of oligomers that may be formed after the disaggregation of fibrils, to complement the results obtained computationally.
1.1. Target Proteins
To carry out this study, the proteins A
and
-syn (represented in the form of fibrils) were considered as targets, as shown in
Figure 2. For beta-amyloid, the structure with PDB ID: 2BEG was chosen, which is a 3D structure of A
(18-42) fibrils determined by Luhrs et al. [
18] in 2005 using solid-state nuclear magnetic resonance spectroscopy. The fibrils are composed of five twisted beta sheets, named A, B, C, D, E, and F, with the A
(18-42) peptides stacked on top of each other in a ladder-like fashion. The beta sheets are stabilized by a network of hydrogen bonds and hydrophobic interactions. The fibrils are highly polymorphic, meaning that they can exist in different shapes and sizes. However, all fibrils share the same basic structure of twisted beta sheets.
For
-syn, the structure with PDB ID: 6H6B was chosen, which is a 3D structure of
-syn(38-95) fibrils determined by Guerrero-Ferreira et al. [
19] in 2005 using cryo-electron microscopy. This structure shows that
-syn(38-95) fibrils are formed by a parallel arrangement of two sets of beta sheets. The beta sheets are composed of long, thin protein filaments arranged in a zigzag pattern. The chains are held together by hydrogen bonds, which are weak bonds that form between oxygen atoms in the protein structure. The beta sheets are stacked on top of each other and the strands of one sheet run perpendicular to the strands of the next sheet. One set is named by the sequence of sheets C, B, A, H, and G, and the other by the sequence F, D, E, J, and I.
3. Results and Discussion
The results of the molecular docking showed that Octovespin and Fraternine-10 bind to the same region of both A
(17-42) and
-syn(38-95) fibrils, as shown in
Figure 3. In the case of A
(17-42), the ligands docked along the axis of the fibril between one of the lateral fragments of the
- sheet. Fraternine-10 docked along the axis of the fibril with the lateral fragment of the A sheet of the fibril in the most central part, interacting with the amino acids Leu17(A), Leu34(A), Phe19(A), Phe20(A), Gly25(A), Gly38(A), Glu22(A), Asp23(A), and Val24(A). Octovespin docked in a similar region, but shifted towards the edge of the A filament, with interactions with the A, B, and C sheets, and interacting with the amino acids Asp23(A), Ala21(A), Val18(A), Val36(A), Phe20(A), Ala21(B), Leu34(B), Val36(B), and Val36(C). The score function, which measures the binding affinity, of the docking of A
(17-42) with Fraternine-10 had a value of -5.1 kcal/mol and that of Octovespin of -6.2 kcal/mol. In the case of
-syn(38-95), both ligands docked near the region of interaction between the two sets of fibrils, in one of the pockets, near the intersection region of the set of fibrils, formed by the sheets C, B, A, H, and G. Fraternine-10 showed interaction with the amino acids Lys60(C), Thr59(C), Val63(C), Glu61(C), Gly73(C), Val63(C), Ala56(B), Glu61(B), Thr72(B), Gly73(B), Thr75(B), Gly73(A), Thr75(A), Glu61(A), Thr59(A), Thr59(H), Glu61(H), Val55(H), and Gly73(H). Octovespin showed interaction with the amino acids Gly73(C), Thr59(C), Glu61(C), Thr75(B), Gly73(B), Ala56(B), Thr59(B), Thr75(A), Thr75(H), Thr59(H), Ala56(H), and Glu61(H). The score function, which measures the binding affinity, of the docking of
-syn(38-95) with Fraternine-10 had a value of -6.4 kcal/mol and that of Octovespin of -7.9 kcal/mol.
Subsequently, molecular dynamics simulations were performed to understand how Fraternine-10 and Octovespin interact with the target fibrils. The changes in secondary structure were analyzed during the dynamic simulation for complexes of each compound, as well as for isolated target fibrils. Differences were observed in the conformation of the fibrils in the presence of the compounds, mainly about the A
(17-42) fibril. The conformational changes are presented visually in
Figure 4 and
Figure 5, with poses of the isolated fibrils and in the presence of the compounds at times of 0ns, 50ns, and 100ns of simulation.
Based on these results, the root-mean-square deviation (RMSD) of the fibrils was calculated to understand how the structural changes occurred concerning the initial conformation. Numerically, the average RMSD values obtained were 4.536Å for A
(17-42) alone, 5.661Å for A
(17-42) with Fraternine-10, 4.906Å for A
(17-42) with Octovespin, 2.823Å for
-syn(38-95) alone, 2.834Å for
-syn(38-95) with Fraternine-10, and 3.663Å for
-syn(38-95) with Octovespin. For the A
(17-42) fibril, the interaction with Fraternine-10 showed considerable changes in the first 20 ns of simulation, decreasing the RMSD value after that time and up to 100 ns. Between 60-70 ns there was an increase in the peak, as shown in
Figure 6. With the interaction with Octovespin, small conformational changes occurred until approximately 85 ns, where from that moment, the RMSD showed a sharp trend of increasing in its value, indicating a greater change in the conformation of the fibril after that time. About the
-syn(38-95) fibril, the RMSD of the fibril alone and in the presence of the compounds showed small variations, but with emphasis on its value when in the presence of Octovespin, which showed a more considerable change in the RMSD. This is because
-syn(38-95) fibrils are relatively larger than A
(17-42) fibrils and, therefore, have less mobility and greater resistance to conformational change. In addition to RMSD, the number of hydrogen bonds in the intermolecular structure of the fibrils was also calculated, presented employing the graph in
Figure 6. According to what was observed, for the A
(17-42) fibril, the decrease in bonds is noticeable in the presence of Fraternine-10, mainly after approximately 60 ns of simulation, indicating a disaggregating effect of this compound. To the
-syn(38-95) fibril, there can also be noted a decrease in the number of hydrogen bonds in the intermolecular structure, mainly around 80 ns and after, but in smaller quantity.
The root-mean-square fluctuation (RMSF) of the fibrils was also calculated in order to understand how the change in the flexibility of the fibrils occurs in the presence of peptides. The obtained results are shown in
Figure 6. For A
(17-42) fibrils, residues 1-26 represent the A sheet, 27-52 the B sheet, 53-78 the C sheet, 79-104 the D sheet, and 105-130 the E sheet. The results indicate a increase in fibril flexibility in the presence of Fraternine-10, with RMSF peaks pronounced close to the central region of the fibril cavity. For Octovespin, the results indicated little change in the flexibility of the fibril, following the same behavior of the isolated fibril. For
-syn(38-95) fibrils, residues 1-58 represent the C sheet, 59-116 the B sheet, 117-174 the A sheet, 175-232 the H sheet, 233-290 the sheet G, 291-348 the F sheet, 349-406 the D sheet, 407-464 the E sheet, 465-522 the J sheet, and 523-580 the I sheet. Our findings suggest that both peptides slightly influence the flexibility of
-syn(38-95). Notably, Octovespin exhibited a more pronounced effect, transiently increasing the fibrils RMSF close to the interface region during the simulation.
Subsequently, an analysis of the interactions of the compounds with A
(17-42) and
-syn(38-95) fibrils at 0, 50, and 100 ns of the simulation was performed, the interactions of which are presented in
Figure 7 and
Figure 8. For Fraternine-10 and A
(17-42), during the 100 ns of simulation, the interactions occurred primarily with the amino acids Val36 (
-alkyl), Ala21 (hydrogen bond) and Leu17 (hydrogen bond). It is worth noting that a large part of the interactions also occurred with the amino acids Val39 (
-alkyl) and Val40 (
-alkyl), which are amino acids responsible for the aggregation of A
(17-42) fibrils. Concerning Fraternine-10 and
-syn(38-95), during the 100 ns of simulation, the interactions occur primarily with the amino acids Thr59 (hydrogen bond), Glu61 (
-anion) and Gly73 (carbon hydrogen bond). To Octovespin and A
(17-42) during the 100 ns of simulation, the interactions occurred primarily with the amino acid Phe19 (hydrogen bond). For Octovespin and
-syn(38-95), during the 100 ns of simulation, the interactions occurred primarily with the amino acids Ala56 (
-alkyl) and Thr75 (hydrogen bond). It is worth noting that a large part of the interactions occurred with the amino acid Ala56 (
-alkyl) which is one of the amino acids responsible for the aggregation of
-syn fibrils.
Considering the ThT assay, Fraternine-10 intermediate (37.5
M) and low dose (18.8
M) had ThT signal decrease compared to A
(1-42) alone (two-way ANOVA; F (864, 9216) = 2.048; p<0.0001), as shown in
Figure 9. However, the signal started to increase after 12h which indicated the ability of Fraternine-10 to delay amyloid aggregation but not to prevent it. Finally, Fraternine-10 high dose (75
M) did not affect amyloid aggregation. According to Camargo and collaborators [
17], Octovespin also prevented amyloid aggregation in the ThT assay.
Certain natural peptides also showed promising results both
in silico and
in vitro in the A
aggregation. Carnosine, a peptide naturally produced in muscles and brain, has demonstrated a reduction in the formation of A
fibrils, as evidenced by atomic force microscopy (AFM) and the ThT assay. Moreover,
in silico data showed that carnosine and homocarnosine interact in D23 and K28 of the A
sequence, which are amino acid residues involved in the self-association process [
32,
33].
While several peptides have been specifically designed to prevent amyloid aggregation, a category called “
-sheet breakers” had limitations throughout the drug development. Synthetic pseudo-peptides with around 8 amino acid residues (Sg inhibitor) were designed to bind in the aggregation-prone region A
(16–23). The evaluation of MD simulations showed that SGC1 had more favorable binding to A
although demonstrated inefficacy against A
toxicity in cells. Interestingly, other peptides with less favorable binding had a better protective effect [
34]. In the present study, both Fraternine-10 and octovespin showed good binding to A
and were effective in the ThT assay and in vivo [
17]. Similar to fraternine-10, a dipeptide D-Trp-Aib showed also promising results in the MD simulation and probable interaction with Lis17 and Phe19 of A
hydrophobic core [
35].
Regarding
in silico assessment of peptide interactions with
-syn, only a limited number of studies are available. Other molecules have been studied, such as Copper (Cu (II)). Cu(II) binds to
-syn seems to increase conformation stability [
36]. For instance, a study showed that 5,7-dihydroxy-2-[(2R,3R)-3-(hydroxymethyl)-2-(4-hydroxyphenyl)-2,3-dihydro-1,4-benzodioxin-6-yl]chromen-4-one can alter
-syn conformation by molecular dynamics evaluation and interaction with Glu57 [
37]. Similarly, derivatives of 7,8-Dihydroxyflavone had high binding affinity and interacted with Gln99, Ala90, Ser87, Val95, Gly93, Leu100 and Lys97 amoung other amino acids residues [
38]. In the present study, the impact of Octovespin and Fraternine-10 on
-syn conformation was assessed through RMSD and molecular dynamics. Octovespin induced conformational changes and exhibited interactions with Ala56 and Thr75 amino acid residues. Fraternine-10, on the other hand, induced subtle changes in
-synuclein conformation and interacted with Thr59, Glu61, and Gly73. In summary, both Octovespin and Fraternine-10 demonstrated good interaction with
-syn and A
, which shows their potential as compounds of interest in further exploration and development.
Author Contributions
“Conceptualization, Y.A.O.S., L.C.C., M.R.M., and R.G.; methodology, Y.A.O.S., C.V.S.C., L.C.C., L.G.V, L.A.R.J, M.R.M, and R.G.; software, Y.A.O.S., C.V.S.C., L.C.C., and L.A.R.J; validation, Y.A.O.S., C.V.S.C., L.C.C., M.R.M, and R.G.; formal analysis, Y.A.O.S., C.V.S.C., L.C.C., L.G.V, M.R.M, and R.G.; investigation, Y.A.O.S., L.C.C., L.G.V, M.R.M, and R.G.; resources, L.A.R.J, M.R.M, and R.G.; data curation, Y.A.O.S., C.V.S.C., L.C.C., and L.G.V.; writing—original draft preparation, Y.A.O.S., L.C.C., L.A.R.J, M.R.M, and R.G.; writing—review and editing, Y.A.O.S., L.C.C., L.A.R.J, M.R.M, and R.G.; visualization, Y.A.O.S. and L.C.C; supervision, L.A.R.J, M.R.M, and R.G.; project administration, L.A.R.J, M.R.M, and R.G.; funding acquisition, L.A.R.J, M.R.M, and R.G. All authors have read and agreed to the published version of the manuscript.”