Constructing conformational library for amyloid-β42 dimers as the smallest toxic oligomers using two CHARMM force fields

Disrupting Dimeric β-Amyloid by Electric Fields

The early oligomers of the amyloid Aβ peptide are implicated in Alzheimer's disease, but their transient nature complicates the characterization of their structure and toxicity. Here, we investigate the stability of the minimal toxic species, i.e., β-amyloid dimers, in the presence of an oscillating electric field. We first use deep learning (AlphaFold-multimer) for generating initial models of Aβ42 dimers. The flexibility and secondary structure content of the models are then analyzed by multiple runs of molecular dynamics (MD). Structurally stable models are similar to ensemble representatives from microsecond-long MD sampling. Finally, we employ the validated model as the starting structure of MD simulations in the presence of an external oscillating electric field and observe a fast decay of β-sheet content at high field strengths. Control simulations using the helical dimer of the 42-residue leucine zipper peptide show higher structural stability than the Aβ42 dimer. The simulation results provide evidence that an external electric field (oscillating at 1 GHz) can disrupt amyloid oligomers which should be further investigated by experiments with brain organoids in vitro and eventually in vivo.

Identification of Catechins' Binding Sites in Monomeric Aβ42 through Ensemble Docking and MD Simulations

The assembly of the amyloid-β peptide (Aβ) into toxic oligomers and fibrils is associated with Alzheimer's disease and dementia. Therefore, disrupting amyloid assembly by direct targeting of the Aβ monomeric form with small molecules or antibodies is a promising therapeutic strategy. However, given the dynamic nature of Aβ, standard computational tools cannot be easily applied for high-throughput structure-based virtual screening in drug discovery projects. In the current study, we propose a computational pipeline-in the framework of the ensemble docking strategy-to identify catechins' binding sites in monomeric Aβ₄₂. It is shown that both hydrophobic aromatic interactions and hydrogen bonding are crucial for the binding of catechins to Aβ₄₂. Additionally, it has been found that all the studied ligands, especially EGCG, can act as potent inhibitors against amyloid aggregation by blocking the central hydrophobic region of Aβ. Our findings are evaluated and confirmed with multi-microsecond MD simulations. Finally, it is suggested that our proposed pipeline, with low computational cost in comparison with MD simulations, is a suitable approach for the virtual screening of ligand libraries against Aβ.

Interactions of Curcumin's Degradation Products with the Aβ42 Dimer: A Computational Study

Amyloid-β (Aβ) dimers are the smallest toxic species along the amyloid-aggregation pathway and among the most populated oligomeric accumulations present in the brain affected by Alzheimer's disease (AD). A proposed therapeutic strategy to avoid the aggregation of Aβ into higher-order structures is to develop molecules that inhibit the early stages of aggregation, i.e., dimerization. Under physiological conditions, the Aβ dimer is highly dynamic and does not attain a single well-defined structure but is rather characterized by an ensemble of conformations. In a recent study, a highly heterogeneous library of conformers of the Aβ dimer was generated by an efficient sampling method with constraints based on ion mobility mass spectrometry data. Here, we make use of the Aβ dimer library to study the interaction with two curcumin degradation products, ferulic aldehyde and vanillin, by molecular dynamics (MD) simulations. Ensemble docking and MD simulations are used to provide atomistic detail of the interactions between the curcumin degradation products and the Aβ dimer. The simulations show that the aromatic residues of Aβ, and in particular ¹⁹FF²⁰, interact with ferulic aldehyde and vanillin through π-π stacking. The binding of these small molecules induces significant changes on the ¹⁶KLVFF²⁰ region.