Insights into the Inhibitory Mechanism of Dicyanovinyl-Substituted J147 Derivative against Aβ42 Aggregation and Protofibril Destabilization: A Molecular Dynamics Simulation Study

Molecular insights into the structure destabilization effects of ECG and EC on the Aβ protofilament: An all-atom molecular dynamics simulation study

The formation of Aβ into amyloid fibrils was closely connected to AD, therefore, the Aβ aggregates were the primary therapeutic targets against AD. Previous studies demonstrated that epicatechin-3-gallate (ECG), which possessed a gallate moiety, exhibited a greater ability to disrupt the preformed Aβ amyloid fibrils than epicatechin (EC), indicating that the gallate moiety was crucial. In the present study, the molecular mechanisms were investigated. Our results demonstrated that ECG had more potent disruptive impacts on the β-sheet structure and K28-A42 salt bridges than EC. We found that ECG significantly interfered the interactions between Peptide-4 and Peptide-5. However, EC could not. The disruption of K28-A42 salt bridges by ECG was mainly due to the interactions between ECG and the hydrophobic residues located at C-terminus. Interestingly, EC disrupted the K28-A42 salt bridges by the interactions with C-terminal hydrophobic residues and the cation-π interactions with K28. Moreover, our results indicated that hydrophobic interactions, H-bonds, π-π interactions and cation-π interactions between ECG and the bend of L-shaped region caused the disaggregation of interactions between Peptide-4 and Peptide-5. Significantly, gallate moiety in ECG had contributed tremendously to the disaggregation. We believed that our findings could be useful for designing prospective drug candidates targeting AD.

Identification of new pentapeptides as potential inhibitors of amyloid-β42 aggregation using virtual screening and molecular dynamics simulations

Alzheimer's disease (AD) is a multifactorial neurodegenerative disease mainly characterized by extracellular accumulation of amyloid-β (Aβ) peptide. Previous studies reported pentapeptide RIIGL as an effective inhibitor of Aβ aggregation and neurotoxicity induced by Aβ aggregates. In this work, a library of 912 pentapeptides based on RIIGL has been designed and assessed for their efficacy to inhibit Aβ42 aggregation using computational techniques. The top hit pentapeptides revealed by molecular docking were further assessed for their binding affinity with Aβ42 monomer using MM-PBSA (molecular mechanics Poisson-Boltzmann surface area) method. The MM-PBSA analysis identified RLAPV, RVVPI, and RIAPA, which bind to Aβ42 monomer with a higher binding affinity -55.80, -46.32, and -44.26 kcal/mol, respectively, as compared to RIIGL (ΔGbinding = -41.29 kcal/mol). The residue-wise binding free energy predicted hydrophobic contacts between Aβ42 monomer and pentapeptides. The secondary structure analysis of the conformational ensembles generated by molecular dynamics (MD) depicted remarkably enhanced sampling of helical and no β-sheet conformations in Aβ42 monomer on the incorporation of RVVPI and RIAPA. Notably, RVVPI and RIAPA destabilized the D23-K28 salt bridge in Aβ42 monomer, which plays a crucial role in Aβ42 oligomer stability and fibril formation. The MD simulations highlighted that the incorporation of proline and arginine in pentapeptides contributed to their strong binding with Aβ42 monomer. Furthermore, RVVPI and RIAPA prevented conformational conversion of Aβ42 monomer to aggregation-prone structures, which, in turn, resulted in a lower aggregation tendency of Aβ42 monomer.

In silico design and identification of new peptides for mitigating hIAPP aggregation in type 2 diabetes

The aberrant misfolding and self-aggregation of human islet amyloid polypeptide (hIAPP or amylin) into cytotoxic aggregates are implicated in the pathogenesis of type 2 diabetes (T2D). Among various inhibitors, short peptides derived from the amyloidogenic regions of hIAPP have been employed as hIAPP aggregation inhibitors due to their low immunogenicity, biocompatibility, and high chemical diversity. Recently, hIAPP fragment HSSNN18-22 was identified as an amyloidogenic sequence and displayed higher antiproliferative activity to RIN-5F cells. Various hIAPP aggregation inhibitors have been designed by chemical modifications of the highly amyloidogenic sequence (NFGAIL) of hIAPP. In this work, a library of pentapeptides based on fragment HSSNN18-22 was designed and assessed for their efficacy in blocking hIAPP aggregation using an integrated computational screening approach. The binding free energy calculations by molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method identified HSSQN and HSSNQ that bind to hIAPP monomer with a binding affinity of -21.25 ± 4.90 and -19.73 ± 3.10 kcal/mol, respectively, which is notably higher as compared to HSSNN (-11.90 ± 4.12 kcal/mol). The sampling of the non aggregation-prone helical conformation was notably increased from 23.5 ± 3.0 in the hIAPP monomer to 38.1 ± 3.6, and 33.8 ± 3.0% on the incorporation of HSSQN, and HSSNQ, respectively, which indicate reduced aggregation propensity of hIAPP monomer. The pentapeptides, HSSQN and HSSNQ, identified as hIAPP aggregation inhibitors in this work can be further conjugated with various metal chelating peptides to yield more efficacious and clinically relevant multifunctional modulators for targeting various pathological hallmarks of T2D.

Unveiling How Hydroxytyrosol Destabilizes α-Syn Oligomers Using Molecular Simulations

The etiology of Parkinson's disease (PD) is mainly linked to the α-synuclein (α-Syn) fibrillogenesis. Hydroxytyrosol (HT), also known as 3,4-dihydroxyphenylethanol, is a naturally occurring polyphenol, found in extra virgin olive oil, and has been shown to have cardioprotective, anticancer, antiobesity, and antidiabetic properties. HT has neuroprotective benefits in neurodegenerative diseases and lessens the severity of PD by reducing the aggregation of α-Syn and destabilizing the preformed toxic α-Syn oligomers. However, the molecular mechanism by which HT destabilizes α-Syn oligomers and alleviates the accompanying cytotoxicity remains unexplored. The impact of HT on the α-Syn oligomer structure and its potential binding mechanism was examined in this work by employing molecular dynamics (MD) simulations. The secondary structure analysis depicted that HT significantly reduces the β-sheet and concomitantly increases the coil content of α-Syn trimer. Visualization of representative conformations from the clustering analysis depicted the hydrogen bond interactions of the hydroxyl groups in HT with the N-terminal and nonamyloid-β component (NAC) region residues of α-Syn trimer, which, in turn, leads to the weakening of interchain interactions in α-Syn trimer and resulted in the disruption of the α-Syn oligomer. The binding free energy calculations depict that HT binds favorably to α-Syn trimer (ΔG_binding = -23.25 ± 7.86 kcal/mol) and a notable reduction in the interchain binding affinity of α-Syn trimer on the incorporation of HT, which, in turn, highlights its potential to disrupt α-Syn oligomers. The current research provided mechanistic insights into the destabilization of α-Syn trimer by HT, which, in turn, will provide new clues for developing therapeutics against PD.

Enhanced stability of a disaggregated Aβ fibril on removal of ligand inhibits refibrillation: An all atom Molecular Dynamics simulation study

The extraneuronally deposited senile plaques, composed of neurotoxic aggregates of Aβ fibril, define Alzheimer's disease (AD). Natural compounds have been tested for their destabilization potential on Aβ fibril, thereby curing AD. However, the resultant destabilized Aβ fibril, needs to be checked for its irreversibility to the native organized state after removal of the ligand. Herein, we assessed the stability of a destabilized fibril after the ligand (ellagic acid represented as REF) is removed from the complex. The study has been conducted via Molecular Dynamics (MD) simulation of 1 μs for both Aβ-Water (control) and Aβ-REF″ (test or REF removed) system. The increased value of RMSD, R_g, SASA, lower β-sheet content and reduced number of H-bonds explains enhanced destabilization observed in Aβ-REF″ system. The increased inter-chain distance demonstrates breaking of the residual contacts, testifying the drift of terminal chains from the pentamer. The increased SASA along with the ∆G_ps(polar solvation energy) accounts for the reduced interaction amongst residues, and more with solvent molecules, governing irreversibility to native state. The higher Gibb's free energy of the misaligned structure of Aβ-REF″ ensures irreversibility to the organized structure due to its inability to cross such high energy barrier. The observed stability of the disaggregated structure, despite ligand elimination, establishes the effectiveness of the destabilization technique as a promising therapeutic approach towards treating AD.

Structural and molecular insights into tacrine-benzofuran hybrid induced inhibition of amyloid-β peptide aggregation and BACE1 activity

Amyloid-β (Aβ) aggregation and β-amyloid precursor protein cleaving enzyme 1 (BACE1) are the potential therapeutic drug targets for Alzheimer's disease (AD). A recent study highlighted that tacrine-benzofuran hybrid C1 displayed anti-aggregation activity against Aβ42 peptide and inhibit BACE1 activity. However, the inhibition mechanism of C1 against Aβ42 aggregation and BACE1 activity remains unclear. Thus, molecular dynamics (MD) simulations of Aβ42 monomer and BACE1 with and without C1 were performed to inspect the inhibitory mechanism of C1 against Aβ42 aggregation and BACE1 activity. In addition, a ligand-based virtual screening followed by MD simulations was employed to explore potent new small-molecule dual inhibitors of Aβ42 aggregation and BACE1 activity. MD simulations highlighted that C1 promotes the non aggregating helical conformation in Aβ42 and destabilizes D23-K28 salt bridge that plays a vital role in the self-aggregation of Aβ42. C1 displays a favourable binding free energy (-50.7 ± 7.3 kcal/mol) with Aβ42 monomer and preferentially binds to the central hydrophobic core (CHC) residues. MD simulations highlighted that C1 strongly interacted with the BACE1 active site (Asp32 and Asp228) and active pockets. The scrutiny of interatomic distances among key residues of BACE1 highlighted the close flap (non-active) position in BACE1 on the incorporation of C1. The MD simulations explain the observed high inhibitory activity of C1 against Aβ aggregation and BACE1 in the in vitro studies. The ligand-based virtual screening followed by MD simulations identified CHEMBL2019027 (C2) as a promising dual inhibitor of Aβ42 aggregation and BACE1 activity.Communicated by Ramaswamy H. Sarma.

Unravelling the destabilization potential of ellagic acid on α-synuclein fibrils using molecular dynamics simulations

The aberrant deposition of α-synuclein (α-Syn) protein into the intracellular neuronal aggregates termed Lewy bodies and Lewy neurites characterizes the devastating neurodegenerative condition known as Parkinson's disease (PD). The disruption of pre-existing disease-relevant α-Syn fibrils is recognized as a viable therapeutic approach for PD. Ellagic acid (EA), a natural polyphenolic compound, is experimentally proven as a potential candidate that prevents or reverses the α-Syn fibrillization process. However, the detailed inhibitory mechanism of EA against the destabilization of α-Syn fibril remains largely unclear. In this work, the influence of EA on α-Syn fibril and its putative binding mechanism were explored using molecular dynamics (MD) simulations. EA interacted primarily with the non-amyloid-β component (NAC) of α-Syn fibril, disrupting its β-sheet content and thereby increasing the coil content. The E46-K80 salt bridge, critical for the stability of Greek-key-like α-Syn fibril, was disrupted in the presence of EA. The binding free energy analysis using the MM-PBSA method demonstrates the favourable binding of EA to α-Syn fibril (ΔG_binding = -34.62 ± 11.33 kcal mol^-1). Interestingly, the binding affinity between chains H and J of the α-Syn fibril was significantly reduced on the incorporation of EA, which highlights the disruptive ability of EA towards α-Syn fibril. The MD simulations provide mechanistic insights into the α-Syn fibril disruption by EA, which gives a valuable direction for the development of potential inhibitors of α-Syn fibrillization and its associated cytotoxicity.

Mechanistic insights into the mitigation of Aβ aggregation and protofibril destabilization by a D–enantiomeric decapeptide rk10

According to clinical studies, the development of Alzheimer’s disease (AD) is linked to the abnormal aggregation of amyloid-β (Aβ) peptides into toxic soluble oligomers, protofibrils as well as mature fibrils....

Insights into β-amyloid transition prevention by cucurbit[7]uril from molecular modeling

In this study, comparable molecular dynamic (MD) simulations of 1.2 microseconds were performed to clarify the prevention of the β-amyloid peptide (Aβ_1-42) aggregation by cucurbit[7]uril (CB[7]). The accumulation of this peptide in the brain is one of the most harmful in Alzheimer's disease. The inhibition mechanism of Aβ_1-42 aggregation by different molecules is attributed to preventing of Aβ_1-42 conformational transition from α-helix to the β-sheet structure. However, our structural analysis shows that the pure water and aqueous solution of the CB[7] denature the native Aβ_1-42 α-helix structure forming different compactness and unfolded conformations, not in β-sheet form. On the other hand, in the three CB[7]@Aβ_1-42 complexes, it was observed the encapsulation of N-terminal (Asp1), Lys16, and Val36 by CB[7] along the MD trajectory, and not with aromatic residues as suggested by the literature. Only in one CB[7]@Aβ_1-42 complex was observed stable Asp23-Lys28 salt bridge with an average distance of 0.36 nm. All CB[7]@Aβ_1-42 complexes are very stable with binding free energy lowest than ∼-50 kcal/mol between the CB[7] and Aβ_1-42 monomer from MM/PBSA calculation. Therefore, herein we show that the mechanism of the prevention of elongation protofibril by CB[7] is due to the disruption of the Asp23-Lys28 salt bridge and steric effects of CB[7]@Aβ_1-42 complex with the fibril lattice, and not due to the transition from α-helix to β-sheet following the dock-lock mechanism.Communicated by Ramaswamy H. Sarma.

How Does the Mono-Triazole Derivative Modulate Aβ42 Aggregation and Disrupt a Protofibril Structure: Insights from Molecular Dynamics Simulations

Clinical studies have identified that abnormal self-assembly of amyloid-β (Aβ) peptide into toxic fibrillar aggregates is associated with the pathology of Alzheimer's disease (AD). The most acceptable therapeutic approach to stop the progression of AD is to inhibit the formation of β-sheet-rich structures. Recently, we designed and evaluated a series of novel mono-triazole derivatives 4(a-x), where compound 4v was identified as the most potent inhibitor of Aβ₄₂ aggregation and disaggregates preformed Aβ₄₂ fibrils significantly. Moreover, 4v strongly averts the Cu²⁺-induced Aβ₄₂ aggregation and disaggregates the preformed Cu²⁺-induced Aβ₄₂ fibrils, halts the generation of reactive oxygen species, and shows neuroprotective effects in SH-SY5Y cells. However, the underlying molecular mechanism of inhibition of Aβ₄₂ aggregation by 4v and disaggregation of preformed Aβ₄₂ fibrils remains obscure. In this work, molecular dynamics (MD) simulations have been performed to explore the conformational ensemble of the Aβ₄₂ monomer and a pentameric protofibril structure of Aβ₄₂ in the presence of 4v. The MD simulations highlighted that 4v binds preferentially at the central hydrophobic core region of the Aβ₄₂ monomer and chains D and E of the Aβ₄₂ protofibril. The dictionary of secondary structure of proteins analysis indicated that 4v retards the conformational conversion of the helix-rich structure of the Aβ₄₂ monomer into the aggregation-prone β-sheet conformation. The binding free energy calculated by the molecular mechanics Poisson-Boltzmann surface area method revealed an energetically favorable process with ΔG _binding = -44.9 ± 3.3 kcal/mol for the Aβ₄₂ monomer-4v complex. The free energy landscape analysis highlighted that the Aβ₄₂ monomer-4v complex sampled conformations with significantly higher helical contents (35 and 49%) as compared to the Aβ₄₂ monomer alone (17%). Compound 4v displayed hydrogen bonding with Gly37 (chain E) and π-π interactions with Phe19 (chain D) of the Aβ₄₂ protofibril. Further, the per-residue binding free energy analysis also highlighted that Phe19 (chain D) and Gly37 (chain E) of the Aβ₄₂ protofibril showed the maximum contribution in the binding free energy. The decreased binding affinity and residue-residue contacts between chains D and E of the Aβ₄₂ protofibril in the presence of 4v indicate destabilization of the Aβ₄₂ protofibril structure. Overall, the structural information obtained through MD simulations indicated that 4v stabilizes the native helical conformation of the Aβ₄₂ monomer and persuades a destabilization in the protofibril structure of Aβ₄₂. The results of the study will be useful in the rational design of potent inhibitors against amyloid aggregation.

Molecular insights into the inhibitory mechanisms of gallate moiety on the Aβ1-40 amyloid aggregation: A molecular dynamics simulation study

Alzheimer's disease is the most common form of neurodegenerative disease and the formation of Aβ amyloid aggregates has been widely demonstrated to be the principal cause of Alzheimer's disease. Our previous study and other studies suggested that the gallate moiety played an obligatory role in the inhibition process of naturally occurring polyphenols on Aβ amyloid fibrils formation. However, the detailed mechanisms were still unknown. Thus, in the present study, the gallic acid (GA) was specially selected and the molecular recognition mechanisms between GA molecules and Aβ_1-40 monomer were examined and analyzed by molecular dynamics simulation. The in silico experiments revealed that GA significantly prevented the conformational changes of Aβ_1-40 monomer with no β-sheet structure during the whole 100 ns. By analyzing the binding sites of GA molecules to Aβ_1-40 monomer, we found that both hydrophilic and hydrophobic amino acid residues were participated in the binding of GA molecules to Aβ_1-40 monomer. Moreover, results from the binding free energy analysis further demonstrated that the strength of polar interactions was significantly stronger than that of nonpolar interactions. We believed that our results could help to elucidate the underlying mechanisms of gallate moiety on the anti-amyloidogenic effects of polyphenols at the atomic level.

Caffeine destabilizes preformed Aβ protofilaments: insights from all atom molecular dynamics simulations

The aggregation and deposition of neurotoxic Aβ fibrils are key in the etiology of Alzheimer's disease (AD). It has been clinically recognized as a major form of dementia across the globe. Finding and testing various natural compounds to target Aβ fibrils to disrupt their stable structures seems to be a promising and attractive therapeutic strategy. The destabilization effects of caffeine on Aβ fibrils are investigated via in silico studies, where a series of molecular dynamics (MD) simulations, each of 100 ns, was conducted. The simulation outcomes obtained henceforth clearly indicated the drift of the terminal chains from the protofibrils, leading to disorganization of the characteristically organized cross-β structures of Aβ fibrils. The structural instability of Aβ17-42 protofibrils is explained through enhanced fluctuations in the RMSD, radius of gyration and RMSF values in the presence of caffeine. The key interactions providing stability, comprising D23-K28 salt bridges, intra- and inter-chain hydrogen bonding and hydrophobic interactions involving interchain A21-V36 and F19-G38 and intrachain L34-V36, were found to be disrupted due to increases in the distances between the participating components. The loss of β-sheet structure with the introduction of turns and α-helices in terminal chains may further inhibit the formation of higher order aggregates, which is necessary to stop the progression of the disease. The atomistic details obtained via MD studies relating to the mechanism behind the underlying destabilization of Aβ17-42 protofibrils by caffeine encourage further investigations exploring the potency of natural compounds to treat AD via disrupting preformed neurotoxic Aβ protofibrils.

Edaravone inhibits the conformational transition of amyloid-β42: insights from molecular dynamics simulations

Previous work has shown that edaravone inhibits fibrillogenesis of amyloid-β protein (Aβ). However, the detailed mechanism by which edaravone inhibits the conformational transition of the Aβ42 monomer is not known at the molecular level. Here, explicit-solvent molecular dynamics (MD) simulations were coupled with molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method to address the issue. MD simulations confirmed that edaravone inhibits the conformational transition of the Aβ42 monomer in a dose-dependent manner. It was found that the direct interactions between edaravone and Aβ42 are responsible for its inhibiting effects. The analysis of binding free energy using the MM-PBSA method demonstrated that the nonpolar interactions provide favourable contributions (about -71.7 kcal/mol). Conversely, the polar interactions are unfavourable for the binding process. A total of 14 residues were identified as greatly contributing to the binding free energy between edaravone and the Aβ42 monomer. In addition, the intra-peptide hydrophobic interactions were weakened and the salt bridge D23-K28 was interrupted by edaravone. Therefore, the conformational transition was inhibited. Our studies provide molecular-level insights into how edaravone molecules inhibit the conformational transition of the Aβ42 monomer, which may be useful for designing amyloid inhibitors.Communicated by Ramaswamy H. Sarma.

Molecular Simulations Reveal Terminal Group Mediated Stabilization of Helical Conformers in Both Amyloid-β42 and α-Synuclein

The presence of partially structured helices in natively unfolded amyloid-β42 (Aβ42) and α-synuclein (αS) has been shown to accelerate fibrillation in the onset of Alzheimer's and Parkinson's disease, respectively. At the other extreme, folded stable helical conformers have also been reported to resist amyloid formation. Recent studies indicate that amyloidogenic aggregation can be impeded using small molecules that stabilize the α-helical monomers and switch off the neurotoxic pathway. We predict a common intrapeptide route to stabilization based on the plasticity of helical conformations of Aβ42 and αS as assessed through extensive atomistic molecular dynamics (MD) computer simulations (∼36 μs) across ten distinct protein force field and water model combinations. Computed free energies and interaction maps (not obtainable from experiments alone) show that flexible terminal groups (N-terminus of Aβ42 and C-terminus of αS) show a tendency to stabilize folded helical conformations in both peptides via primary hydrophobic interactions with central hydrophobic domains, and secondary salt bridges with other domains. These interactions confer aggregation resistance by decreasing the population of partially structured helices and are absent in control simulations of complete unfolding. Computed helical stability is also significantly reduced in terminal-deleted variants. The models suggest new strategies to tackle neurodegeneration by rationally re-engineering terminal groups to optimize their predicted ability to deactivate helical monomers.

Conformational Dynamics and Stability of U-Shaped and S-Shaped Amyloid β Assemblies

Alzheimer’s disease is the most fatal neurodegenerative disorder characterized by the aggregation and deposition of Amyloid β (Aβ) oligomers in the brain of patients. Two principal variants of Aβ exist in humans: Aβ_1–40 and Aβ_1–42. The former is the most abundant in the plaques, while the latter is the most toxic species and forms fibrils more rapidly. Interestingly, fibrils of Aβ_1–40 peptides can only assume U-shaped conformations while Aβ_1–42 can also arrange as S-shaped three-stranded chains, as recently discovered. As alterations in protein conformational arrangement correlate with cell toxicity and speed of disease progression, it is important to characterize, at molecular level, the conformational dynamics of amyloid fibrils. In this work, Replica Exchange Molecular Dynamics simulations were carried out to compare the conformational dynamics of U-shaped and S-shaped Aβ_17–42 small fibrils. Our computational results provide support for the stability of the recently proposed S-shaped model due to the maximized interactions involving the C-terminal residues. On the other hand, the U-shaped motif is characterized by significant distortions resulting in a more disordered assembly. Outcomes of our work suggest that the molecular architecture of the protein aggregates might play a pivotal role in formation and conformational stability of the resulting fibrils.

Molecular modeling studies of quinazolinone derivatives as novel PI3Kδ selective inhibitors

The main molecular modeling method, the docking results of newly designed compoundD04and the best pharmacophore model are reported herein.