Atomic and dynamic insights into the beneficial effect of the 1,4-naphthoquinon-2-yl-L-tryptophan inhibitor on Alzheimer's Aβ1-42 dimer in terms of aggregation and toxicity

Virtual and In Vitro Screens Reveal a Potential Pharmacophore that Avoids the Fibrillization of Aβ1-42

Among the multiple factors that induce Alzheimer’s disease, aggregation of the amyloid β peptide (Aβ) is considered the most important due to the ability of the 42-amino acid Aβ peptides (Aβ_1–42) to form oligomers and fibrils, which constitute Aβ pathological aggregates. For this reason, the development of inhibitors of Aβ_1–42 pathological aggregation represents a field of research interest. Several Aβ_1–42 fibrillization inhibitors possess tertiary amine and aromatic moieties. In the present study, we selected 26 compounds containing tertiary amine and aromatic moieties with or without substituents and performed theoretical studies that allowed us to select four compounds according to their free energy values for Aβ_1–42 in α-helix (Aβ-α), random coil (Aβ-RC) and β-sheet (Aβ-β) conformations. Docking studies revealed that compound 5 had a higher affinity for Aβ-α and Aβ-RC than the other compounds. In vitro, this compound was able to abolish Thioflavin T fluorescence and favored an RC conformation of Aβ_1–42 in circular dichroism studies, resulting in the formation of amorphous aggregates as shown by atomic force microscopy. The results obtained from quantum studies allowed us to identify a possible pharmacophore that can be used to design Aβ_1–42 aggregation inhibitors. In conclusion, compounds with higher affinity for Aβ-α and Aβ-RC prevented the formation of oligomeric species.

Coupling of Zinc-Binding and Secondary Structure in Nonfibrillar Aβ40 Peptide Oligomerization

Nonfibrillar neurotoxic amyloid β (Aβ) oligomer structures are typically rich in β-sheets, which could be promoted by metal ions like Zn(2+). Here, using molecular dynamics (MD) simulations, we systematically examined combinations of Aβ40 peptide conformations and Zn(2+) binding modes to probe the effects of secondary structure on Aβ dimerization energies and kinetics. We found that random conformations do not contribute to dimerization either thermodynamically or kinetically. Zn(2+) couples with preformed secondary structures (α-helix and β-hairpin) to speed dimerization and stabilize the resulting dimer. Partial α-helices increase the dimerization speed, and dimers with α-helix rich conformations have the lowest energy. When Zn(2+) coordinates with residues D1, H6, H13, and H14, Aβ40 β-hairpin monomers have the fastest dimerization speed. Dimers with experimentally observed zinc coordination (E11, H6, H13, and H14) form with slower rate but have lower energy. Zn(2+) cannot stabilize fibril-like β-arch dimers. However, Zn(2+)-bound β-arch tetramers have the lowest energy. Collectively, zinc-stabilized β-hairpin oligomers could be important in the nucleation-polymerization of cross-β structures. Our results are consistent with experimental findings that α-helix to β-structural transition should accompany Aβ aggregation in the presence of zinc ions and that Zn(2+) stabilizes nonfibrillar Aβ oligomers and, thus, inhibits formation of less toxic Aβ fibrils.

Structural Conversion of Aβ17-42 Peptides from Disordered Oligomers to U-Shape Protofilaments via Multiple Kinetic Pathways

Discovering the mechanisms by which proteins aggregate into fibrils is an essential first step in understanding the molecular level processes underlying neurodegenerative diseases such as Alzheimer's and Parkinson's. The goal of this work is to provide insights into the structural changes that characterize the kinetic pathways by which amyloid-β peptides convert from monomers to oligomers to fibrils. By applying discontinuous molecular dynamics simulations to PRIME20, a force field designed to capture the chemical and physical aspects of protein aggregation, we have been able to trace out the entire aggregation process for a system containing 8 Aβ17-42 peptides. We uncovered two fibrillization mechanisms that govern the structural conversion of Aβ17-42 peptides from disordered oligomers into protofilaments. The first mechanism is monomeric conversion templated by a U-shape oligomeric nucleus into U-shape protofilament. The second mechanism involves a long-lived and on-pathway metastable oligomer with S-shape chains, having a C-terminal turn, en route to the final U-shape protofilament. Oligomers with this C-terminal turn have been regarded in recent experiments as a major contributing element to cell toxicity in Alzheimer's disease. The internal structures of the U-shape protofilaments from our PRIME20/DMD simulation agree well with those from solid state NMR experiments. The approach presented here offers a simple molecular-level framework to describe protein aggregation in general and to visualize the kinetic evolution of a putative toxic element in Alzheimer's disease in particular.

The mechanisms of flavonoids inhibiting conformational transition of amyloid-β42monomer: a comparative molecular dynamics simulation study

Flavonoids can bind Aβ₄₂to inhibit the aggregation of Aβ₄₂monomer.

Molecular mechanism of the inhibition and remodeling of human islet amyloid polypeptide (hIAPP(1-37)) oligomer by resveratrol from molecular dynamics simulation

Natural polyphenols are one of the most actively investigated categories of amyloid inhibitors, and resveratrol has recently been reported to inhibit and remodel the human islet amyloid polypeptide (hIAPP) oligomers and fibrils. However, the exact mechanism of its action is still unknown, especially for the full-length hIAPP1-37. To this end, we performed all-atom molecular dynamics simulations for hIAPP1-37 pentamer with and without resveratrol. The obtained results show that the binding of resveratrol is able to cause remarkable conformational changes of hIAPP1-37 pentamer, in terms of secondary structures, order degree, and morphology. By clustering analysis, two possible binding sites of resveratrol on the hIAPP1-37 pentamer were found, located at the grooves of the top and bottom surfaces of β-sheet layer, respectively. After the binding free energy calculation and residue energy decomposition, it can be concluded that the bottom site is the more possible one, and that the nonpolar interactions act as the driving force for the binding of hIAPP1-37 to resveratrol. In addition, Arg11 is the most important residue for the binding of resveratrol. The full understanding of inhibitory mechanism of resveratrol on the hIAPP1-37 oligomer, and the identification of its binding sites on this protein are helpful for the future design and discovery of new amyloid inhibitors.

Alzheimer's disease--a panorama glimpse

The single-mutation of genes associated with Alzheimer's disease (AD) increases the production of Aβ peptides. An elevated concentration of Aβ peptides is prone to aggregation into oligomers and further deposition as plaque. Aβ plaques and neurofibrillary tangles are two hallmarks of AD. In this review, we provide a broad overview of the diverses sources that could lead to AD, which include genetic origins, Aβ peptides and tau protein. We shall discuss on tau protein and tau accumulation, which result in neurofibrillary tangles. We detail the mechanisms of Aβ aggregation, fibril formation and its polymorphism. We then show the possible links between Aβ and tau pathology. Furthermore, we summarize the structural data of Aβ and its precursor protein obtained via Nuclear Magnetic Resonance (NMR) or X-ray crystallography. At the end, we go through the C-terminal and N-terminal truncated Aβ variants. We wish to draw reader's attention to two predominant and toxic Aβ species, namely Aβ4-42 and pyroglutamate amyloid-beta peptides, which have been neglected for more than a decade and may be crucial in Aβ pathogenesis due to their dominant presence in the AD brain.

The OPEP protein model: from single molecules, amyloid formation, crowding and hydrodynamics to DNA/RNA systems

The OPEP coarse-grained protein model has been applied to a wide range of applications since its first release 15 years ago. The model, which combines energetic and structural accuracy and chemical specificity, allows the study of single protein properties, DNA-RNA complexes, amyloid fibril formation and protein suspensions in a crowded environment. Here we first review the current state of the model and the most exciting applications using advanced conformational sampling methods. We then present the current limitations and a perspective on the ongoing developments.

Molecular structure of the NQTrp inhibitor with the Alzheimer Aβ1-28 monomer

The self-assembly of the amyloid-β (Aβ) peptide of various amino acid lengths into senile plaques is one hallmark of Alzheimer's disease pathology. In the past decade, many small molecules, including NQTrp, have been identified to reduce aggregation and toxicity. However, due to the heterogeneity of the conformational ensemble of Aβ with drugs, we lack detailed structures of the transient complexes. Following our previous simulation of the monomer of Aβ1-28, here we characterize the equilibrium ensemble of the Aβ1-28 monomer with NQTrp by means of extensive atomistic replica exchange molecular dynamics simulations using a force field known to fold diverse proteins correctly. While the secondary structure content and the intrinsic disorder of the whole peptides are very similar and the lifetimes of the salt-bridges remain constant, the population of β-hairpin is reduced by a factor of 1.5 and the population of α-helix in the region 17-24 is increased by a factor of two upon NQTrp binding. These two factors, which impact the free energy barrier for nucleation, provide a first explanation for the reported reduced Aβ1-40/1-42 aggregation kinetics in the presence of NQTrp. Backbone and side-chain interactions of Aβ with NQTrp may also inhibit Aβ-Aβ contacts. The fraction of free Aβ1-28 monomer is, however, on the order of 20-25% at 17.5 mM, and this shows that the affinity of NQTrp is low and hence its inhibitory activity is not very strong. This inhibitor can be improved to reduce the formation of dimer, a critical step in aggregation and toxicity.

Understanding amyloid fibril nucleation and aβ oligomer/drug interactions from computer simulations

Evolution has fine-tuned proteins to accomplish a variety of tasks. Yet, with aging, some proteins assemble into harmful amyloid aggregates associated with neurodegenerative diseases, such as Alzheimer's disease (AD), which presents a complex and costly challenge to our society. Thus, far, drug after drug has failed to slow the progression of AD, characterized by the self-assembly of the 39-43 amino acid β-amyloid (Aβ) protein into extracellular senile plaques that form a cross-β structure. While there is experimental evidence that the Aβ small oligomers are the primary toxic species, standard tools of biology have failed to provide structures of these transient, inhomogeneous assemblies. Despite extensive experimental studies, researchers have not successfully characterized the nucleus ensemble, the starting point for rapid fibril formation. Similarly scientists do not have atomic data to show how the compounds that reduce both fibril formation and toxicity in cells bind to Aβ42 oligomers. In this context, computer simulations are important tools for gaining insights into the self-assembly of amyloid peptides and the molecular mechanism of inhibitors. This Account reviews what analytical models and simulations at different time and length scales tell us about the dynamics, kinetics, and thermodynamics of amyloid fibril formation and, notably, the nucleation process. Though coarse-grained and mesoscopic protein models approximate atomistic details by averaging out unimportant degrees of freedom, they provide generic features of amyloid formation and insights into mechanistic details of the self-assembly process. The thermodynamics and kinetics vary from linear peptides adopting straight β-strands in fibrils to longer peptides adopting in parallel U shaped conformations in fibrils. In addition, these properties change with the balance between electrostatic and hydrophobic interactions and the intrinsic disorder of the system. However, simulations suggest that the critical nucleus size might be on the order of 20 chains under physiological conditions. The transition state might be characterized by a simultaneous change from mixed antiparallel/parallel β-strands with random side-chain packing to the final antiparallel or parallel states with the steric zipper packing of the side chains. Second, we review our current computer-based knowledge of the 3D structures of inhibitors with Aβ42 monomer and oligomers, a prerequisite for developing new drugs against AD. Recent extensive all-atom simulations of Aβ42 dimers with known inhibitors such as the green tea compound epigallocatechin-3-gallate and 1,4-naphthoquinon-2-yl-l-tryptophan provide a spectrum of initial Aβ42/inhibitor structures useful for screening and drug design. We conclude by discussing future directions that may offer opportunities to fully understand nucleation and further AD drug development.