Next Article in Journal
Calcium Phosphate Particles Coated with Humic Substances: A Potential Plant Biostimulant from Circular Economy
Next Article in Special Issue
Cheminformatic Profiling and Hit Prioritization of Natural Products with Activities against Methicillin-Resistant Staphylococcus aureus (MRSA)
Previous Article in Journal
Structural Relevance of Intramolecular H-Bonding in Ortho-Hydroxyaryl Schiff Bases: The Case of 3-(5-bromo-2-hydroxybenzylideneamino) Phenol
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Curcumin and Ferulic Acid on the Folding of Amyloid-β Peptide

Faculty of Pharmacy, Medical University of Sofia, 1000 Sofia, Bulgaria
*
Author to whom correspondence should be addressed.
Submission received: 29 April 2021 / Revised: 6 May 2021 / Accepted: 8 May 2021 / Published: 10 May 2021
(This article belongs to the Special Issue In Silico Activity Profiling of Natural Products)

Abstract

:
The polyphenols curcumin (CU) and ferulic acid (FA) are able to inhibit the aggregation of amyloid-β (Aβ) peptide with different strengths. CU is a strong inhibitor while FA is a weaker one. In the present study, we examine the effects of CU and FA on the folding process of an Aβ monomer by 1 µs molecular dynamics (MD) simulations. We found that both inhibitors increase the helical propensity and decrease the non-helical propensity of Aβ peptide. They prevent the formation of a dense bulk core and shorten the average lifetime of intramolecular hydrogen bonds in Aβ. CU makes more and longer-lived hydrogen bonds, hydrophobic, π–π, and cation–π interactions with Aβ peptide than FA does, which is in a good agreement with the observed stronger inhibitory activity of CU on Aβ aggregation.

1. Introduction

Polyphenols are a group of natural compounds known for their numerous beneficial effects on different functions of the human body. In addition to their well-known antioxidant properties [1], new effects have recently been identified. Several studies have found that flavonoids and phenolic acids are able to reduce the risk of type 2 diabetes by inhibition of glucose absorption [2,3], stimulation of insulin secretion [4,5] and reduction of hepatic glucose output [6,7]. Similar mechanisms underline the antiobesity effects of polyphenols [8]. The antioxidant activity of polyphenols is associated with anti-inflammatory effects [9] and reduced risk of cardiovascular disease [10]. It has been found that polyphenol-rich plant extracts have anticoagulant properties [11] and anti-cancer effects [12,13]. Green tea polyphenols have digestion-improving activity by promoting the growth of probiotic bacteria and inhibiting the growth of pathogenic ones [14]. The polyphenols in grape juice [15], cocoa [16] and Ginkgo biloba [17] are able to improve brain activity and boost memory and concentration.
Curcumin (CU, Figure 1) is a natural polyphenol originating from the plant turmeric (Curcuma longa) [18]. Most of the beneficial effects of CU are due to its antioxidant [19] and anti-inflammatory properties [20]. CU is a free radical scavenger [19] and chain-breaking antioxidant [21]. It is an inhibitor of several enzymes involved in the generation of reactive oxygen species like lipoxygenase, cyclooxygenase and xanthine hydrogenase [22]. Based on its strong antioxidant properties, CU acts as an anti-inflammatory agent and prevents the development of many chronic diseases like neurodegeneration, cancer, osteoarthritis [18]. It was found that CU prevents the aggregation of amyloid-beta (Aβ) peptide [23] and tau proteins [24] by direct binding and prevention of their assemblage into neurotoxic species in the brain. The IC50 value of CU for in vitro inhibition of Aβ aggregation is 0.81 µM [23]. Moreover, CU is able to disintegrate preformed Aβ fibrils with EC50 of 1 µM [23].
Ferulic acid (FA, Figure 1) is a phenolic acid with antioxidant [25], hypotensive [26] and anti-inflammatory properties [27] found in various fruits and vegetables. It has been shown that FA inhibits aldose reductase, a key enzyme involved in the development of insulin resistance and metabolic syndrome [28]. Additionally, FA has neuroprotective properties like antidepressant [29,30,31], antinociceptive [32] and antiepileptic effects [33], protection from ischemia [34,35], Parkinson’s disease [36] and inflammation [37,38]. FA is a weak inhibitor of the Aβ aggregation as well with IC50 of 5.5 µM [39].
The effects of CU on the structure of Aβ oligomers and fibrils have been studied by molecular dynamics (MD) simulations. MD is a computational method for in silico mimicking the movement of a molecule in a given medium [40]. The effects of CU on the stability of dimers and protofibrils have been simulated and the mechanisms of interactions have been investigated at atomistic level [41,42,43,44,45,46,47]. Recently, we simulated the primary nucleation process of 12 Aβ peptides and studied the effects of CU and FA on the process [48]. We found that CU intercalates among the peptide chains, binds tightly to Aβ by hydrogen bonds, hydrophobic, π–π, and cation–π interactions and prevents the formation of a compact primary nucleus. The interactions of FA with Aβ are weak and short-living in accordance to the weaker inhibitory effect of FA on the Aβ aggregation.
Here, we examine by MD simulation the effect of CU on the secondary structure and folding of a single Aβ molecule and identify the interactions between the two molecules. The process of Aβ packing in saline is conducted for a duration of 1 µs (1000 ns) in the presence and absence of CU. For comparison, the same process for the same time and at the same conditions is conducted in the presence of the weaker inhibitor FA.

2. Models and Methods

2.1. Modeled Ligands and Systems

The structures of the two ligands CU (CID 969516) and FA (CID 445858) were downloaded from PubChem [49]. The CU was taken in tautomeric keto-enol form as it has been found that this form is able to prevent Aβ fibril formation [50,51]. The FA acid was simulated in anionic form as it is a weak acid with pKa = 4.58. Structurally, FA mirrors the half-cutoff molecule of CU.
The point charges of CU and FA were derived as AM1-BCC charges using the AMBER module antechamber [52]. The topology of the Aβ peptide was built using the ff14SB force field. CU and FA were parameterized with GAFF atom type [53].
The Aβ1-42 monomer structure was obtained from PDB (pdb code 1IYT) [54]. Both ends were capped and all hydrogens were added. Three systems were arranged: Aβ peptide, Aβ peptide and CU, Aβ peptide and FA (Figure 2). The ligands were placed randomly near the peptide and the systems were solvated in a truncated octahedron water box using periodic boundary conditions (PBC) with TIP3P water explicit solvent molecules. Na+ and Cl- were added to neutralize the structures and achieve a salt concentration of 0.1 mol/L to mimic the physiological pH.

2.2. Molecular Dynamics Protocol

The MD simulations in the present study were performed using AMBER 18 [56]. As a first step in the MD protocol, the three systems were minimized for 10,000 steps using 10 kcal/molÅ2 harmonic restraints on solute heavy atoms and 10 Å cutoff for the non-bonded van der Waals and electrostatic interactions. The structures were heated to 300 K in a timestep of 2 fs with releasing the restraints to 3 kcal/molÅ2 for 100 ps. The next step was a system density equilibration for 100 ps under the same force constant on the heavy atoms. The restraints were released in 1 ns equilibration at constant temperature maintained with Langevin thermostat (300 K) and constant pressure using the Berendsen barostat (1 bar). The production run was held for 1 µs for any of the three systems. Frames were saved every 1 ns for a total of 1000 per trajectory.
The trajectories were analyzed by cpptraj V4.24.0 [57]. The following parameters were calculated: root mean square deviation (RMSD) accounts for changes in atomic coordinates of input frames to a reference frame (frame #1 by default); RMSF (root mean square fluctuations) accounts for atomic positional fluctuations; secondary structural propensities for residues (α-, 310 and 314 helices, parallel and anti-parallel β-sheets, turns and bends); solvent accessible surface area (SASA) accounts for the surface area in Å2 of atoms; hydrogen bond formation. The RMSD and RMSF were calculated for the backbone atoms of Aβ peptide. The calculation of secondary structure propensities was based on the DPSS method of Kabsch and Sander [58]. The SASA was calculated according to thr Weiser, Shenkin, and Still LCPO approximation method [59]. The electrostatic potentials on the SASA were calculated by Poisson-Boltzmann equation and visualized by PBEQ Solver [60]. The hydrogen bonds were calculated by applying simple geometric criteria for Acceptor-H-Donor distance cutoff of 3.0 Å and dihedral cutoff of 135°.

3. Results

The three systems–Aβ peptide, Aβ peptide and CU, Aβ peptide and FA–were solvated in water, energy-minimized, heated to 300 K, equilibrated and simulated for 1 µs (1000 ns) as described in the Models and Methods section. The topologies of the modelled systems before and after the production phase are given in Figure 3. The derived trajectories were analysed to assess the effects of CU and FA on the folding of Aβ peptide.

3.1. Curcumin Stabilizes the Aβ Structure

The backbone RMSDs of the Aβ peptide in the three studied systems as a function of time are given in Figure 4 (left panels). Initially, the single Aβ molecule reaches equilibrium for 100 ns, stays stable for the next 200 ns, between 300 and 500 ns undergoes a distortion and after 500 ns it gradually stabilizes again. The difference between maximum and minimum Aβ backbone RMSDs over the last 100 frames of the simulation is 4.762 Å. In the presence of CU, the peptide is stabilized for the first 50 ns, stays stable for the next 300 ns, between 400 and 600 ns undergoes a distortion and after 700 ns is stable again till the end of simulation. The difference between maximum and minimum Aβ backbone RMSDs over the last 100 frames of the simulation is only 1.621 Å. In the presence of FA, the peptide distortions continue until 700 ns and then it stabilizes with difference between maximum and minimum backbone RMSDs over the last 100 frames of 4.383 Å. Obviously, the addition of CU to the solvated Aβ stabilizes the peptide structure.
The backbone RMSF values per residue showed significant differences between the three simulated systems (Figure 4, right panels). The RMSF of the single Aβ peptide ranges from 4.90 to 8.38 Å with average value of 6.32 (±0.80) Å. The fluctuations are spread almost evenly along the sequence. In the presence of CU, the peptide structure is stabilized in the middle part between residues 9 and 32, while the both ends remain flanking. The RMSF of the peptide ranges from 2.70 to 9.39 Å with average value of 5.11 (±1.43) Å. FA has a weaker stabilizing effect focused mainly on the fragments 9–21 and 28–34 with RMSF values in the range 3.57–9.84 Å and average of 5.96 (±1.43) Å.

3.2. Curcumin and Ferulic Acid Increase the Helical Structure Propensity of Aβ Peptide

The propensities of secondary structures of the modelled systems are given as average values over the total 1000 frames (1000 ns) in Table 1. A-helix and parallel β-turn are the dominating structures in the single Aβ peptide, followed by anti-parallel β-turn. 310 helices and bends are rare. In the presence of CU, the propensity of α-helix increases, the parallel β-turn decreases, the anti-parallel β-turn almost disappears, bends increase slightly and 310 helix does not change. Similar changes were observed in the presence of FA.
The propensities of α- and 310-helices were summed and presented as helical structure propensity, while the sum of propensities of parallel and anti-parallel β-turns and bends gives the non-helical propensity (Table 1). Both CU and FA increase the helical propensity and decrease the non-helical propensity of Aβ peptide.
The helical and non-helical propensities averaged over the initial 100 ns and over the final 100 ns of the production phase of the MD simulation are presented along the peptide sequence in Figure 5. Initially, three helices exist in the structure of Aβ peptide (Figure 3A). They include residues 2–5, 11–23 and 28–36. At the end of simulation, the first helix has disappeared, the second has shrinnked between positions 12–21 and the third has split into two shorter helices between position 25–28 and 31–34. Between the helices non-helical domains are spread. Initially, they involve positions 6–10, 24–27 and 37–39. At the final, the non-helical structures are localized at the two peptide ends including positions 4–10 and 32–40. Two anti-parallel β-turns are formed here as is shown in Figure 3D.
During the initial 100 frames (100 ns) of the production phase in the presence of CU, two well-defined α-helices cover almost the whole structure of Aβ peptide and two short non-helical domains exist including residues 26, 27, 38 and 39 (Figure 3B). At the end of simulation, an α-helix exists only between positions 13–24, while the remaining part of the peptide contains non-helical structures (Figure 3E).
In the presence of FA, the initial structure of Aβ consists of three helices and three non-helical structures (Figure 3C). At the final 100 frames, the peptide still contains three shorter helical structures but at the end of simulation the last one is disordered (Figure 3F).

3.3. Curcumin and Ferulic Acid Increase the Solvent-Accessible Surface Area (SASA)

SASA accounts for the density of packing of the peptide molecule in saline. The average SASA over the production trajectory of a single Aβ molecule is 3237 Å2. In the presence of CU, the average SASA of Aβ increases to 3586 Å2. In the presence of FA, it reaches 3679 Å2. Both ligands prevent the formation of a dense bulk core and disfavor the growth of packing.
The level of packing of the Aβ peptide is visualized by SASAs averaged over every 100 ns for the three analysed systems (Figure 6). The single Aβ molecule formes in a two-step manner a compact nucleus which slightly relaxes at the end of simulation. In the precense of CU, the peptide initially shrinks, then relaxes, followed by small decrease and increase in SASA. Similar behaivior has Aβ in the presence of FA.
The electrostatic potentials on the SASA calculated by Poisson-Boltzmann equation and visualized by PBEQ Solver showed that the compact nucleus of Aβ peptide formed in the three systems consists of a hydrophobic core and a hydrophilic polar surface (Figure 7). As an entropy-driven process, the packing of Aβ in saline forms nuclei of different shapes and surface arrangements.

3.4. Curcumin Forms More and Longer-Lived Hydrogen Bonds with Aβ Peptide

During the 1 µs simulation, 444 intramolecular hydrogen bonds in total were recorded between the Aβ residues with average lifetime of 30.53 ns (Table 2). Among the longest-lived hydrogen bonds are those between Gln15 and Phe19, Glu11 and Gln15, and Val24 and Lys28 with lifetimes >400 ns.
In the presence of CU, the intramolecular hydrogen bonds in Aβ are preserved as a number but decreased as an average lifetime to 26.59 ns. Ninety four short-living intermolecular hydrogen bonds are formed between Aβ and CU with average lifetime of 5.46 ns. In 65% of these bonds, CU acts as a donor through the phenolic OH groups and the enol group. The longest-lived hydrogen bonds between Aβ and CU are given in Figure 8. As corresponding acceptors in these bonds act the backbone O-atoms of Phe4, His13, Ala2, Gly33 and Leu34, the side-chain O-atoms of Glu11, Asp7, Glu22 and Asp23, the imidazole N-atom (Nδ) of His13.
The presence of FA in the vicinity of Aβ increases the number of intramolecular hydrogen bonds but also descreases the average lifetime to 25.91 ns (Table 2). Only 55 hydrogen bonds between Aβ and FA were detected during the simulation with ultra-short average lifetime of 2.67 ns. In 75% of them FA acts as an acceptor mainly through the carboxy O-atoms (Figure 8). The longest-lived hydrogen bonds involve the guanidino group of Arg5.

4. Discussion

The single monomer Aβ peptide is a water-soluble non-toxic molecule. In an apolar microevironment, it exists as two α-helices connected by a flexible kink [54] as is shown in Figure 2A. In water, the peptide adopts a collapsed compact coil structure [61]. This structure is meta-stable and is able to be easily arranged into oligomers, protofibrils and fibrils containing predominantly intermolecular β-sheet structures which are water insoluble and neurotoxic [62]. Thus, the prevention of nascent Aβ monomer from misfolding and conversion to toxic conformers can serve as a rational therapeutic goal.
It was proven by in vitro studies that the polyphenol CU can block the formation of Aβ oligomers [23]. Additionally, CU is able to cross the blood-brain barrier and to bind to and disintegrate preformed amyloid plaques in mice brains [23]. However, the submolecular mechanisms of the interactions between CU and Aβ and the effects of CU on the process of folding of a single Aβ molecule are still unclear. Even more unclear are the effects of the phenolic acid FA on the folding of Aβ. FA could be considered as a half-cutoff molecule of CU with almost 7-fold weaker inhibitory activity on Aβ aggregation [39].
The aim of the present study is to analyze the effects of CU and FA on the folding process of a single Aβ peptide molecule in a 1:1 ratio in saline. CU is a strong inhibitor of Aβ aggregation, while FA is a weak inhibitor. Additionally, CU is a neutral molecule, while FA is an anion at pH 7.4. The two molecules were positioned randomly in a close proximity to Aβ peptide and the interactions between them were conducted in an octahedron box of water molecules and NaCl for a time of 1 µs. The movements and interactions were described quantitatively by parameters calculated over the trajectories. For comparison, the movement of a single Aβ peptide in the same box for the same time was also conducted.
The quantitative parameters derived from the MD simulation reveal several similarities and dissimilarities in the effects of the two inhibitors on the folding of Aβ. CU has a stronger stabilizing effect on Aβ structure than FA. Stabilization is achieved for 50 ns. Except for both ends, the whole main structure of Aβ is stable in the presence of CU. In the presence of FA, two peptide fragments remain stable: residues 9–21 and 28–34. The middle part of the peptide (residues 22–27) and both ends are rather flexible.
Both ligands have similar effects on the secondary structure of Aβ. They increase the propensity of α-helices and decrease the propensities of β-turns. In the presence of CU, the helical structure involves the middle peptide part (residues 13–24). In the presence of FA, this helical structure splits into two short helices including residues 10–16 and 18–23 with a bend at position 17. In the absence of ligands, the helical structure of Aβ gradually decreases and converts either into random coil or into a hairpin β-strand. At the end of simulation, two β-strands are formed (residues 4–10 and 32–40). The formation of such regions of β-turns is considered as the first step, from which the subsequent assembly of misfolded peptides proceeds [63].
Both ligands increase the SASA of Aβ peptide with 11% (CU) and 14% (FA). The Aβ nuclei formed in the absence or presence of ligands are similar–they consist of hydrophobic core and hydrophilic polar surface of different shape and arrangement. The number of intramolecular hydrogen bonds in Aβ peptide is not affected by the ligands but both shorten the average lifetimes of the bonds. CU binds to Aβ by 94 hydrogen bonds with average lifetime of 5.46 ns. In most of them, CU acts as a hydrogen-bond donor mainly through the phenolic OH and enol groups. The longest-lived bonds exist between 30 and 40 ns. On the contrary, FA forms half of the CU bonds with shorter average lifetime of only 2.67 ns. In these bonds, FA acts dominantly as a hydrogen-bond acceptor through the carboxy O-atoms. The longest-living hydrogen bonds of FA last up to 15 ns.
The top three longest-lived interactions of CU and FA with Aβ peptide are given in Figure 9. CU makes hydrogen bonds with Phe4 (Figure 9A, yellow dashes), His13 (Figure 9B) and Glu11 (Figure 9C). FA binds bidentately to Arg5 (Figure 9D) making two parallel hydrogen bonds between the carboxy O-atoms of FA and the N-atoms of guanidino group of Arg5. Apart from the hydrogen bonds, the two ligands are involved in a variety of intermolecular interactions with the peptide. CU makes hydrophobic contacts with Glu3, Phe4, Arg5, His6, Asp7, Met35 (Figure 9A, green lines), with His13, Lys16, Phe20, Val36, Val40 (Figure 9B), Tyr10, Leu17, Ile41 (Figure 9C). Additionally, CU is involved in π-π interactions with Phe4 (Figure 9A, red lines), His13 and Phe20 (Figure 9B), Tyr10 and His14 (Figure 9C) and in cation-π interactions with Arg5 (Figure 9A,C, blue lines) and Lys16 (Figure 9B).
In comparison with CU, FA makes fewer interactions with the Aβ peptide (Figure 9D). Apart from the bidentate hydrogen bond, FA interacts hydrophobically with Lys16 and Val36, participates in π-π stacking with Phe19 and cation-π bond with Arg5.

5. Conclusions

In conclusion, the polyphenols CU and FA affect the folding of Aβ peptide by increasing the helical and decreasing the non-helical propensities. CU forms numerous long-living hydrogen bonds, hydrophobic, π-π and cation-π interactions with the peptide, while the interactions of FA are fewer and short-living. These discrepancies in the interactions with Aβ are in a good agreement with the observed stronger inhibitory activity of CU on Aβ aggregation and the weaker activity of FA.

Author Contributions

Conceptualization, I.D. (Irini Doytchinova) and E.S.; methodology, E.S.; investigation, E.S., M.A. and I.D. (Ivan Dimitrov); writing—original draft preparation, E.S. and I.D. (Irini Doytchinova); and writing—review and editing, E.S., M.A., I.D. (Ivan Dimitrov) and I.D. (Irini Doytchinova); visualization, I.D. (Irini Doytchinova); supervision, I.D. (Irini Doytchinova); project administration, I.D. (Irini Doytchinova); funding acquisition, I.D. (Irini Doytchinova). All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Bulgarian National Science Fund (Grant DN03/9/2016) and by the Bulgarian National Roadmap for Research Infrastructure (Grant No. D01-271/2019).

Institutional Review Board Statement

Not available.

Informed Consent Statement

Not available.

Data Availability Statement

N/A.

Acknowledgments

This work was performed in the Centre of Excellence for Informatics and ICT supported by the Science and Education for Smart Growth Operational Program and co-financed by the European Union through the European Structural and Investment funds (Grant No BG05M2OP001-1.001-0003).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Sample Availability

Samples of the compounds are not available from the authors.

References

  1. Swallah, M.S.; Sun, H.; Affoh, R.; Fu, H.; Yu, H. Antioxidant potential overviews of secondary metabolites (polyphenols) in fruits. Int. J. Food. Sci. 2020, 2020, 9081686. [Google Scholar] [CrossRef] [PubMed]
  2. Roder, P.V.; Geillinger, K.E.; Zietek, T.S.; Thorens, B.; Koepsell, H.; Daniel, H. The role of sglt1 and GLUT2 in intestinal glucose transport and sensing. PLoS ONE 2014, 9, e89977. [Google Scholar] [CrossRef] [PubMed]
  3. McDougall, G.J.; Shpiro, F.; Dobson, P.; Smith, P.; Blake, A.; Stewart, D. Different polyphenolic components of soft fruits inhibit alpha-amylase and alpha-glucosidase. J. Agric. Food Chem. 2005, 53, 2760–2766. [Google Scholar] [CrossRef] [PubMed]
  4. Henquin, J.C. Triggering and amplifying pathways of regulation of insulin secretion by glucose. Diabetes 2000, 49, 1751–1760. [Google Scholar] [CrossRef] [Green Version]
  5. Cai, E.P.; Lin, J.K. Epigallocatechin gallate (EGCG) and rutin suppress the glucotoxicity through activating IRS2 and AMPK signaling in rat pancreatic beta cells. J. Agric. Food Chem. 2009, 57, 9817–9827. [Google Scholar] [CrossRef]
  6. Waltner-Law, M.E.; Wang, X.L.; Law, B.K.; Hall, R.K.; Nawano, M.; Granner, D.K. Epigallocatechin gallate, a constituent of green tea, represses hepatic glucose production. J. Biol. Chem. 2002, 277, 34933–34940. [Google Scholar] [CrossRef] [Green Version]
  7. Collins, Q.F.; Liu, H.Y.; Pi, J.; Liu, Z.; Quon, M.J.; Cao, W. Epigallocatechin-3-gallate (EGCG), a green tea polyphenol, suppresses hepatic gluconeogenesis through 51-AMP-activated protein kinase. J. Biol. Chem. 2007, 282, 30143–30149. [Google Scholar] [CrossRef] [Green Version]
  8. Azzini, E.; Giacometti, J.; Russo, G.L. Antiobesity Effects of Anthocyanins in Preclinical and Clinical Studies. Oxid. Med. Cell Longev. 2017, 2740364. [Google Scholar] [CrossRef]
  9. Hussain, T.; Tan, B.; Yin, Y.; Blachier, F.; Tossou, M.C.; Rahu, N. Oxidative Stress and Inflammation: What Polyphenols Can Do for Us? Oxid. Med. Cell Longev. 2016, 7432797. [Google Scholar] [CrossRef] [Green Version]
  10. Tangney, C.C.; Rasmussen, H.E. Polyphenols, inflammation, and cardiovascular disease. Curr. Atheroscler. Rep. 2013, 15, 324. [Google Scholar] [CrossRef]
  11. Bijak, M.; Saluk, J.; Szelenberger, R.; Nowak, P. Popular naturally occurring antioxidants as potential anticoagulant drugs. Chem. Biol. Interact. 2016, 257, 35–45. [Google Scholar] [CrossRef]
  12. Wang, D.; Chen, Q.; Liu, B.; Li, Y.; Tan, Y.; Yang, B. Ellagic acid inhibits proliferation and induces apoptosis in human glioblastoma cells. Acta Cir. Bras. 2016, 31, 143–149. [Google Scholar] [CrossRef] [Green Version]
  13. Shailasree, S.; Venkataramana, M.; Niranjana, S.R. Cytotoxic effect of p-coumaric acid on neuroblastoma, N2a cell via generation of reactive oxygen species leading to dysfunction of mitochondria inducing apoptosis and autophagy. Mol. Neurobiol. 2015, 51, 119–130. [Google Scholar] [CrossRef]
  14. Pacheco-Ordaz, R.; Wall-Medrano, A.; Goñi, M.G.; Ramos-Clamont-Montfort, G.; Ayala-Zavala, J.F.; González-Aguilar, G.A. Effect of phenolic compounds on the growth of selected probiotic and pathogenic bacteria. Lett. Appl. Microbiol. 2018, 66, 25–31. [Google Scholar] [CrossRef]
  15. Krikorian, R.; Nash, T.A.; Shidler, M.D.; Shukitt-Hale, B.; Joseph, J.A. Concord grape juice supplementation improves memory function in older adults with mild cognitive impairment. Br. J. Nutr. 2010, 103, 730–734. [Google Scholar] [CrossRef] [Green Version]
  16. Field, D.T.; Williams, C.M.; Butler, L.T. Consumption of cocoa flavanols results in an acute improvement in visual and cognitive functions. Physiol. Behav. 2011, 103, 255–260. [Google Scholar] [CrossRef] [Green Version]
  17. Vauzour, D. Dietary polyphenols as modulators of brain functions: Biological actions and molecular mechanisms underpinning their beneficial effects. Oxid. Med. Cell Longev. 2012, 2012, 914273. [Google Scholar] [CrossRef] [Green Version]
  18. Hewlings, S.J.; Kalman, D.S. Curcumin: A Review of Its Effects on Human Health. Foods 2017, 6, 92. [Google Scholar] [CrossRef]
  19. Sahebkar, A.; Serbanc, M.C.; Ursoniuc, S.; Banach, M. Effect of curcuminoids on oxidative stress: A systematic review and meta-analysis of randomized controlled trials. J. Funct. Foods 2015, 18, 898–909. [Google Scholar] [CrossRef]
  20. Aggarwal, B.B.; Harikumar, K.B. Potential therapeutic effects of curcumin, the anti-inflammatory agent, against neurodegenerative, cardiovascular, pulmonary, metabolic, autoimmune and neoplastic diseases. Int. J. Biochem. Cell Biol. 2009, 41, 40–59. [Google Scholar] [CrossRef] [Green Version]
  21. Priyadarsini, K.I.; Maity, D.K.; Naik, G.H.; Kumar, M.S.; Unnikrishnan, M.K.; Satav, J.G.; Mohan, H. Role of phenolic O-H and methylene hydrogen on the free radical reactions and antioxidant activity of curcumin. Free Radic. Biol. Med. 2003, 35, 475–484. [Google Scholar] [CrossRef]
  22. Marchiani, A.; Rozzo, C.; Fadda, A.; Delogu, G.; Ruzza, P. Curcumin and curcumin-like molecules: From spice to drugs. Curr. Med. Chem. 2014, 21, 204–222. [Google Scholar] [CrossRef] [PubMed]
  23. Yang, F.; Lim, G.P.; Begum, A.N.; Ubeda, O.J.; Simmons, M.R.; Ambegaokar, S.S.; Chen, P.P.; Kayed, R.; Glabe, C.G.; Frautschy, S.A.; et al. Curcumin inhibits formation of amyloid beta oligomers and fibrils, binds plaques, and reduces amyloid in vivo. J. Biol. Chem. 2005, 280, 5892–5901. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Ma, Q.L.; Zuo, X.; Yang, F.; Ubeda, O.J.; Gant, D.J.; Alaverdyan, M.; Teng, E.; Hu, S.; Chen, P.P.; Maiti, P.; et al. Curcumin suppresses soluble tau dimers and corrects molecular chaperone, synaptic, and behavioral deficits in aged human tau transgenic mice. J. Biol. Chem. 2013, 288, 4056–4065. [Google Scholar] [CrossRef] [Green Version]
  25. Kanski, J.; Aksenova, M.; Stoyanova, A.; Butterfield, D.A. Ferulic acid antioxidant protection against hydroxyl and peroxyl radical oxidation in synaptosomal and neuronal cell culture systems in vitro: Structure-activity studies. J. Nutr. Biochem. 2002, 13, 273–281. [Google Scholar] [CrossRef]
  26. Suzuki, A.; Kagawa, D.; Fujii, A.; Ochiai, R.; Tokimitsu, I.; Saito, I. Short- and long-term effects of ferulic acid on blood pressure in spontaneously hypertensive rats. Am. J. Hypertens. 2002, 15, 351–357. [Google Scholar] [CrossRef] [Green Version]
  27. Hirabayashi, T.; Ochiai, H.; Sakai, S.; Nakajima, K.; Terasawa, K. Inhibitory effect of ferulic acid and isoferulic acid on murine interleukin-8 production in response to influenza virus infections in vitro and in vivo. Planta Med. 1995, 61, 221–226. [Google Scholar] [CrossRef]
  28. El-Bassossy, H.; Badawy, D.; Neamatallah, T.; Fahmy, A. Ferulic acid, a natural polyphenol, alleviates insulin resistance and hypertension in fructose fed rats: Effect on endothelial-dependent relaxation. Chemico-Biol. Interact. 2016, 254, 191–197. [Google Scholar] [CrossRef]
  29. Chen, J.; Lin, D.; Zhang, C.; Li, G.; Zhang, N.; Ruan, L.; Yan, Q.; Li, J.; Yu, X.; Xie, X.; et al. Antidepressant-like effects of ferulic acid: Involvement of serotonergic and norepinergic systems. Metab. Brain Dis. 2015, 30, 129–136. [Google Scholar] [CrossRef]
  30. Lenzi, J.; Rodriguez, A.F.; Rós Ade, S.; de Castro, A.B.; de Lima, D.D.; Magro, D.D.; Zeni, A.L. Ferulic acid chronic treatment exerts antidepressant-like effect: Role of antioxidant defense system. Metab. Brain Dis. 2015, 30, 1453–1463. [Google Scholar] [CrossRef]
  31. Liu, Y.M.; Hu, C.Y.; Shen, J.D.; Wu, S.H.; Li, Y.C.; Yi, L.T. Elevation of synaptic protein is associated with the antidepressant-like effects of ferulic acid in a chronic model of depression. Physiol. Behav. 2017, 169, 184–188. [Google Scholar] [CrossRef]
  32. Xu, Y.; Lin, D.; Yu, X.; Xie, X.; Wang, L.; Lian, L.; Fei, N.; Chen, J.; Zhu, N.; Wang, G.; et al. The antinociceptive effects of ferulic acid on neuropathic pain: Involvement of descending monoaminergic system and opioid receptors. Oncotarget 2016, 7, 20455–20468. [Google Scholar] [CrossRef] [Green Version]
  33. Hassanzadeh, P.; Arbabi, E.; Atyabi, F.; Dinarvand, R. Ferulic acid exhibits antiepileptogenic effect and prevents oxidative stress and cognitive impairment in the kindling model of epilepsy. Life Sci. 2017, 179, 9–14. [Google Scholar] [CrossRef]
  34. Zhang, L.; Wang, H.; Wang, T.; Jiang, N.; Yu, P.; Chong, Y.; Fu, F. Ferulic acid ameliorates nerve injury induced by cerebral ischemia in rats. Exp. Ther. Med. 2015, 9, 972–976. [Google Scholar] [CrossRef]
  35. Koh, P.O. Ferulic acid attenuates the down-regulation of MEK/ERK/p90RSK signaling pathway in focal cerebral ischemic injury. Neurosci. Lett. 2015, 588, 18–23. [Google Scholar] [CrossRef]
  36. Ojha, S.; Javed, H.; Azimullah, S.; Khair, S.B.A.; Haque, M.E. Neuroprotective potential of ferulic acid in the rotenone model of Parkinson’s disease. Drug Des. Devel. Ther. 2015, 9, 5499–5510. [Google Scholar]
  37. Wu, J.L.; Shen, M.M.; Yang, S.X.; Wang, X.; Ma, Z.C. Inhibitory effect of ferulic acid on neuroinflammation in LPS-activated microglia. Chin. Pharm. Bull. 2015, 31, 97–102. [Google Scholar]
  38. Szwajgier, D.; Borowiec, K.; Pustelniak, K. The Neuroprotective Effects of Phenolic Acids: Molecular Mechanism of Action. Nutrients 2017, 9, 477. [Google Scholar] [CrossRef] [Green Version]
  39. Ono, K.; Hirohata, M.; Yamada, M. Ferulic acid destabilizes preformed β-amyloid fibrils in vitro. Biochem. Biophys. Res. Commun. 2005, 336, 444–449. [Google Scholar] [CrossRef]
  40. Nair, P.C.; Miners, J.O. Molecular dynamics simulations: From structure function relationships to drug discovery. Silico Pharmacol. 2014, 2, 4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Zhao, L.N.; Chiu, S.W.; Benoit, J.; Chew, L.Y.; Mu, Y. The effect of curcumin on the stability of Aβ dimers. J. Phys. Chem. B 2012, 116, 7428–7435. [Google Scholar] [CrossRef] [PubMed]
  42. Awasthi, M.; Singh, S.; Pandey, V.P.; Dwivedi, U.N. Modulation in the conformational and stability attributes of the Alzheimer’s disease associated amyloid-beta mutants and their favorable stabilization by curcumin: Molecular dynamics simulation analysis. J. Biomol. Struct. Dyn. 2018, 36, 407–422. [Google Scholar] [CrossRef] [PubMed]
  43. Ngo, S.T.; Li, M.S. Curcumin binds to Aβ1-40 peptides and fibrils stronger than ibuprofen and naproxen. J. Phys. Chem. B 2012, 116, 10165–10175. [Google Scholar] [CrossRef]
  44. Kundaikar, H.S.; Degani, M.S. Insights into the Interaction Mechanism of Ligands with Aβ42 Based on Molecular Dynamics Simulations and Mechanics: Implications of Role of Common Binding Site in Drug Design for Alzheimer’s Disease. Chem. Biol. Drug Des. 2015, 86, 805–812. [Google Scholar] [CrossRef]
  45. Tavanti, F.; Pedone, A.; Menziani, M.C. Computational Insight into the Effect of Natural Compounds on the Destabilization of Preformed Amyloid-β(1–40) Fibrils. Molecules 2018, 23, 1320. [Google Scholar] [CrossRef] [Green Version]
  46. Jakubowski, J.M.; Orr, A.A.; Le, D.A.; Tamamis, P. Interactions between Curcumin Derivatives and Amyloid-β Fibrils: Insights from Molecular Dynamics Simulations. J. Chem. Inf. Model. 2020, 60, 289–305. [Google Scholar] [CrossRef]
  47. Bajda, M.; Filipek, S. Computational approach for the assessment of inhibitory potency against beta-amyloid aggregation. Bioorg. Med. Chem. Lett. 2017, 27, 212–216. [Google Scholar] [CrossRef]
  48. Doytchinova, I.; Atanasova, M.; Salamanova, E.; Ivanov, S.; Dimitrov, I. Curcumin Inhibits the Primary Nucleation of Amyloid-Beta Peptide: A Molecular Dynamics Study. Biomolecules 2020, 10, 1323. [Google Scholar] [CrossRef]
  49. Kim, S.; Chen, J.; Cheng, T.; Gindulyte, A.; He, J.; He, S.; Li, Q.; Shoemaker, B.A.; Thiessen, P.A.; Yu, B.; et al. PubChem 2019 update: Improved access to chemical data. Nucleic Acids Res. 2019, 47, D1102–D1109. [Google Scholar] [CrossRef] [Green Version]
  50. Yanagisawa, D.; Shirai, N.; Amatsubo, T.; Taguchi, H.; Hirao, K.; Urushitani, M.; Morikawa, S.; Inubushi, T.; Kato, M.; Kato, F.; et al. Relationship between the tautomeric structures of curcumin derivatives and their Abeta-binding activities in the context of therapies for Alzheimer’s disease. Biomaterials 2010, 31, 4179–4185. [Google Scholar] [CrossRef]
  51. Rao, P.P.; Mohamed, T.; Teckwani, K.; Tin, G. Curcumin Binding to Beta Amyloid: A Computational Study. Chem. Biol. Drug Des. 2015, 86, 813–820. [Google Scholar] [CrossRef]
  52. Jakalian, A.; Bush, B.L.; Jack, D.B.; Bayly, C.I.F. Efcient generation of high-quality atomic charges. AM1-BCC model: I. Method. J. Comput. Chem. 2000, 21, 132–146. [Google Scholar] [CrossRef]
  53. Wang, J.M.; Wolf, R.M.; Caldwell, J.W.; Kollman, P.A.; Case, D.A. Development and Testing of a General Amber Force Field. J. Comput. Chem. 2004, 25, 1157–1174. [Google Scholar] [CrossRef]
  54. Crescenzi, O.; Tomaselli, S.; Guerrini, R.; Salvadori, S.; D’Ursi, A.M.; Temussi, P.A.; Picone, D. Solution structure of the Alzheimer amyloid beta-peptide (1–42) in an apolar microenvironment. Similarity with a virus fusion domain. Eur. J. Biochem. 2002, 269, 5642–5648. [Google Scholar] [CrossRef]
  55. Krieger, E.; Vriend, G. YASARA View—molecular graphics for all devices—from smartphones to workstations. Bioinformatics 2014, 30, 2981–2982. [Google Scholar] [CrossRef] [Green Version]
  56. Case, D.A.; Cheatham, T.E., 3rd; Darden, T.; Gohlke, H.; Luo, R.; Merz, K.M., Jr.; Onufriev, A.; Simmerling, C.; Wang, B.; Woods, R.J. The Amber biomolecular simulation programs. J. Comput. Chem. 2005, 26, 1668–1688. [Google Scholar] [CrossRef] [Green Version]
  57. Roe, D.R.; Cheatham, T.E. PTRAJ and CPPTRAJ: Sofware for Processing and Analysis of Molecular Dynamics Trajectory Data. J. Chem. Theory Comput. 2013, 9, 3084–3095. [Google Scholar] [CrossRef]
  58. Kabsch, W.; Sander, C. Dictionary of protein secondary structure: Pattern recognition of hydrogen-bonded and geometrical features. Biopolymers 1983, 22, 2577–2637. [Google Scholar] [CrossRef]
  59. Weiser, J.; Shenkin, P.S.; Still, W.C. Approximate solvent-accessible surface areas from tetrahedrally directed neighbor densities. Biopolymers 1999, 50, 373–380. [Google Scholar] [CrossRef]
  60. Jo, S.; Vargyas, M.; Vasko-Szedlar, J.; Roux, B.; Im, W. PBEQ-Solver for Online Visualization of Electrostatic Potential of Biomolecules. Nucleic Acids Res. 2008, 36, W270–W275. [Google Scholar] [CrossRef] [Green Version]
  61. Zhang, S.; Iwata, K.; Lachenmann, M.J.; Peng, J.W.; Li, S.; Stimson, E.R.; Lu, Y.; Felix, A.M.; Maggio, J.E.; Lee, J.P. The Alzheimer’s peptide a beta adopts a collapsed coil structure in water. J. Struct. Biol. 2000, 130, 130–141. [Google Scholar] [CrossRef] [PubMed]
  62. Chen, G.F.; Xu, T.H.; Yan, Y.; Zhou, Y.R.; Jiang, Y.; Melcher, K.; Xu, H.E. Amyloid beta: Structure, biology and structure-based therapeutic development. Acta Pharmacol. Sin. 2017, 38, 1205–1235. [Google Scholar] [CrossRef] [PubMed]
  63. Lazo, N.D.; Grant, M.A.; Condron, M.C.; Rigby, A.C.; Teplow, D.B. On the nucleation of amyloid beta-protein monomer folding. Protein Sci. 2005, 14, 1581–1596. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Chemical structures of the amyloid-beta (Aβ) inhibitors (A) curcumin (CU) and (B) ferulic acid (FA). The O-atoms in the molecules are numbered.
Figure 1. Chemical structures of the amyloid-beta (Aβ) inhibitors (A) curcumin (CU) and (B) ferulic acid (FA). The O-atoms in the molecules are numbered.
Molecules 26 02815 g001
Figure 2. Initial topology of the systems modelled in the present study: (A) Aβ peptide, (B) Aβ peptide and CU, (C) Aβ peptide and FA. The peptide is given in cartoon, the ligand–in ball-and-stick. Structures are visualized by YASARA [55].
Figure 2. Initial topology of the systems modelled in the present study: (A) Aβ peptide, (B) Aβ peptide and CU, (C) Aβ peptide and FA. The peptide is given in cartoon, the ligand–in ball-and-stick. Structures are visualized by YASARA [55].
Molecules 26 02815 g002
Figure 3. Topologies of the modelled systems at the beginning (frame #1) and at the end (frame #100,000) of the production phase: (A) Aβ peptide at frame #1, (B) Aβ peptide and CU at frame #1, (C) Aβ peptide and FA at frame #1, (D) Aβ peptide at frame #100,000, (E) Aβ peptide and CU at frame #100,000, (F) Aβ peptide and FA at frame #100,000. The peptide is given in cartoon, the ligand–in ball-and-stick.
Figure 3. Topologies of the modelled systems at the beginning (frame #1) and at the end (frame #100,000) of the production phase: (A) Aβ peptide at frame #1, (B) Aβ peptide and CU at frame #1, (C) Aβ peptide and FA at frame #1, (D) Aβ peptide at frame #100,000, (E) Aβ peptide and CU at frame #100,000, (F) Aβ peptide and FA at frame #100,000. The peptide is given in cartoon, the ligand–in ball-and-stick.
Molecules 26 02815 g003aMolecules 26 02815 g003b
Figure 4. Aβ backbone RMSDs (left panels) and RMSFs per residue (right panels) of the modelled systems: Aβ peptide (blue), Aβ peptide and CU (red), Aβ peptide and FA (green).
Figure 4. Aβ backbone RMSDs (left panels) and RMSFs per residue (right panels) of the modelled systems: Aβ peptide (blue), Aβ peptide and CU (red), Aβ peptide and FA (green).
Molecules 26 02815 g004aMolecules 26 02815 g004b
Figure 5. Propensities of helical and non-helical structures of the modelled systems averaged over the initial 100 ns (dash lines) and the final 100 ns (solid lines) of the production phase.
Figure 5. Propensities of helical and non-helical structures of the modelled systems averaged over the initial 100 ns (dash lines) and the final 100 ns (solid lines) of the production phase.
Molecules 26 02815 g005
Figure 6. SASA averaged over every 100 ns of the production trajectory of Aβ peptide in the three modelled systems.
Figure 6. SASA averaged over every 100 ns of the production trajectory of Aβ peptide in the three modelled systems.
Molecules 26 02815 g006
Figure 7. 3D electrostatic potential map on SASA visualized by PBEQ Solver. (A) Aβ peptide; (B) Aβ peptide and CU; (C) Aβ peptide and FA. Areas with negative potential (−2 kcal/(mol.e)) are given in red, areas with positive potential (+2 kcal/(mol.e))–in blue, areas with neutral potential–in white.
Figure 7. 3D electrostatic potential map on SASA visualized by PBEQ Solver. (A) Aβ peptide; (B) Aβ peptide and CU; (C) Aβ peptide and FA. Areas with negative potential (−2 kcal/(mol.e)) are given in red, areas with positive potential (+2 kcal/(mol.e))–in blue, areas with neutral potential–in white.
Molecules 26 02815 g007
Figure 8. Intermolecular hydrogen bonds between Aβ peptide and ligands CU and FA with lifetime ≥10 ns.
Figure 8. Intermolecular hydrogen bonds between Aβ peptide and ligands CU and FA with lifetime ≥10 ns.
Molecules 26 02815 g008
Figure 9. Interactions between CU and Aβ peptide (AC) and between FA and Aβ peptide (D). The peptide is given in cartoon, the ligand–in ball-and-stick. Hydrogen bonds are given by yellow dashes, hydrophobic interactions–by green lines, π-π interactions–by red lines and cation-π interaction–by blue lines. Interactions are visualized by YASARA [55].
Figure 9. Interactions between CU and Aβ peptide (AC) and between FA and Aβ peptide (D). The peptide is given in cartoon, the ligand–in ball-and-stick. Hydrogen bonds are given by yellow dashes, hydrophobic interactions–by green lines, π-π interactions–by red lines and cation-π interaction–by blue lines. Interactions are visualized by YASARA [55].
Molecules 26 02815 g009
Table 1. Propensities of secondary structures of the modelled systems given as average values over the total 1000 ns (1 µs). Helical propensity is the sum of α- and 310- helices. Non-helical propensity summarizes parallel and anti-parallel β-turns and bending structures.
Table 1. Propensities of secondary structures of the modelled systems given as average values over the total 1000 ns (1 µs). Helical propensity is the sum of α- and 310- helices. Non-helical propensity summarizes parallel and anti-parallel β-turns and bending structures.
Secondary StructureAβ + CUAβ + FA
Alpha helix0.2750.3550.397
310 helix0.0670.0590.069
Anti-parallel β-turn0.1110.0190.009
Parallel β-turn0.2110.1920.205
Bend0.0810.1090.090
Helical0.3420.4140.466
Non-helical0.4030.3200.304
Table 2. Number, percentage and lifetime of intramolecular and intermolecular hydrogen bonds, calculated over 1000 ns.
Table 2. Number, percentage and lifetime of intramolecular and intermolecular hydrogen bonds, calculated over 1000 ns.
Hydrogen BondsAβ + CUAβ + FA
Total444535533
Intramolecular444 (100%)441 (82%)478 (90%)
Average lifetime (ns)30.5326.5925.91
Intermolecular-94 (18%)55 (10%)
Average lifetime (ns)-5.462.67
Ligand is a donor-61 (65%)14 (25%)
Ligand in an acceptor-33 (35%)41 (75%)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Salamanova, E.; Atanasova, M.; Dimitrov, I.; Doytchinova, I. Effects of Curcumin and Ferulic Acid on the Folding of Amyloid-β Peptide. Molecules 2021, 26, 2815. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules26092815

AMA Style

Salamanova E, Atanasova M, Dimitrov I, Doytchinova I. Effects of Curcumin and Ferulic Acid on the Folding of Amyloid-β Peptide. Molecules. 2021; 26(9):2815. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules26092815

Chicago/Turabian Style

Salamanova, Evdokiya, Mariyana Atanasova, Ivan Dimitrov, and Irini Doytchinova. 2021. "Effects of Curcumin and Ferulic Acid on the Folding of Amyloid-β Peptide" Molecules 26, no. 9: 2815. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules26092815

Article Metrics

Back to TopTop