Ion-induced free energy landscapes of Aβ33–42 peptide dimer in wet ionic liquids

Investigating the Restricted Dynamical Environment in and Around Aβ Peptide Oligomers in Aqueous Ionic Liquid Solutions

It is widely believed that the aggregation of amyloid β (Aβ) peptides into soluble oligomers is the root cause behind Alzheimer's disease. In this study, we have performed room-temperature molecular dynamics (MD) simulations of aggregated Aβ oligomers of different sizes (pentamer (O(5)), decamer (O(10)), and hexadecamer (O(16))) in binary aqueous solutions containing 1-butyl-3-methylimidazolium tetrafluoroborate ([BMIM][BF4]) ionic liquid (IL). Investigations have been carried out to obtain a microscopic understanding of the effects of the IL on the dynamic environment around the exterior surfaces and within the confined nanocores of the oligomers. The calculations revealed that in contrast to nearly uniform dynamics near the exterior surface, heterogeneous structural distortions of oligomers of varying sizes and nonuniform distributions of water and IL components within their core volumes modify the core dynamics in a differential manner. It is demonstrated that increasingly restricted mobility of water and IL components is the origin behind the longer time scale of dynamic heterogeneity in and around the oligomers. Importantly, due to the equivalent nondirectional nature of the B-F bonds, BF4- anions are found to reorient on a time scale faster than that of water molecules. Further, the structural relaxation of protein-anion (PA) hydrogen bonds around the oligomers has been found to be correlated with sluggish translational motions of the anions but anticorrelated with their reorientational time scale. In addition, it is quantified that compared to the pure aqueous medium, strengthening of protein-water (PW) hydrogen bonds in the presence of the IL leads to their longer lifetimes.

Conformational Properties of Aβ Peptide Oligomers in Aqueous Ionic Liquid Solution: Insights from Molecular Simulation Studies

Alzheimer's disease is a progressive irreversible neurological disorder with abnormal extracellular deposition of amyloid β (Aβ) peptides in the brain. We have carried out atomistic molecular dynamics simulations to investigate the size-dependent conformational properties of aggregated Aβ oligomers of different orders, namely, pentamer [O(5)], decamer [O(10)], and hexadecamer [O(16)] in aqueous solutions containing the ionic liquid (IL) 1-butyl-3-methylimidazolium tetrafluoroborate ([BMIM][BF4]). The calculations revealed reduced peptide conformational fluctuations in O(5) and O(10) in the presence of the IL. In contrast, the higher order oligomer [O(16)] has been found to exhibit greater structural distortion due to enhanced flexibilities of its peptide units in the presence of the IL. Based on the distributions of the solvent (water) and the cosolvent (IL) components, it is demonstrated that exchange of water by the IL ion pairs at the exterior surface of the oligomers primarily occurs beyond the first layer of surface-bound water molecules. Importantly, a reduced number of relatively weaker peptide salt bridges have been found in O(16) in binary water-IL solution as compared to the other two smaller-sized oligomers [O(5) and O(10)]. Such differential influence of the IL on peptide salt bridges results in less favorable binding free energies of peptide monomers to O(16), which leads to its greater structural distortion and reduced stability compared to those of O(5) and O(10).

Molecular Modelling of Ionic Liquids: Situations When Charge Scaling Seems Insufficient

Charge scaling as an effective solution to the experiment-computation disagreement in molecular modelling of ionic liquids (ILs) could bring the computational results close to the experimental reference for various thermodynamic properties. According to the large-scale benchmark calculations of mass density, solvation, and water-ILs transfer-free energies in our series of papers, the charge-scaling factor of 0.8 serves as a near-optimal option generally applicable to most ILs, although a system-dependent parameter adjustment could be attempted for further improved performance. However, there are situations in which such a charge-scaling treatment would fail. Namely, charge scaling cannot really affect the simulation outcome, or minimally perturbs the results that are still far from the experimental value. In such situations, the vdW radius as an additional adjustable parameter is commonly tuned to minimize the experiment-calculation deviation. In the current work, considering two ILs from the quinuclidinium family, we investigate the impacts of this vdW-scaling treatment on the mass density and the solvation/partition thermodynamics in a fashion similar to our previous charge-scaling works, i.e., scanning the vdW-scaling factor and computing physical properties under these parameter sets. It is observed that the mass density exhibits a linear response to the vdW-scaling factor with slopes close to -1.8 g/mL. By further investigating a set of physiochemically relevant temperatures between 288 K and 348 K, we confirm the robustness of the vdW-scaling treatment in the estimation of bulk properties. The best vdW-scaling parameter for mass density would worsen the computation of solvation/partition thermodynamics, and a marginal decrease in the vdW-scaling factor is considered as an intermediate option balancing the reproductions of bulk properties and solvation thermodynamics. These observations could be understood in a way similar to the charge-scaling situation. i.e., overfitting some properties (e.g., mass density) would degrade the accuracy of the other properties (e.g., solvation free energies). Following this principle, the general guideline for applying this vdW-tuning protocol is by using values between the density-derived choice and the solvation/partition-derived solution. The charge and current vdW scaling treatments cover commonly encountered ILs, completing the protocol for accurate modelling of ILs with fixed-charge force fields.

Helix-coil transition and conformational deformity in Aβ42-monomer: a case study using the Zn2+ cation

The metal ions (like Fe2+, Zn²⁺, Cu2+) are known to influence the amyloid beta (Aβ) aggregation. In this study, we have examined the conformational and dynamical changes during the coordination of Aβ-monomer with the Zn2+ ion using all-atom molecular dynamics (MD) simulations using explicit solvent models. We have probed the unfolding of the full-length Aβ₄₂ monomer both inclusive and exclusive of the Zn2+ cation, with 1:1 ratio of the peptide and the Zn2+ cation. The inclusion of the Zn2+ cation shows differential intra-peptide interactions which has been probed using various analyses. The Helix - Coil transition of the wild type Aβ42 monomer is studied using the steered molecular dynamics simulations by taking the end-to-end C-α distance across the peptide. This gives an idea of the unequal intra - peptide and peptide - water interactions being found across the length of the Aβ monomer. The transition of an α-helix dominated wild-type (WT) Aβ structure to the unfolded coil structure gives significant evidence of the intra-peptide hydrogen bonding shifts in the presence of the Zn2+ cation. This accounts for the structural and the dynamical variations that take place in the Aβ monomer in the presence of the Zn2+ cation to mimic the conditions/environment at the onset of fibrillation.Communicated by Ramaswamy H. Sarma.

Conformational dynamics of amyloid-β (16–22) peptide in aqueous ionic liquids

Molecular dynamics simulations of amyloid-β (16–22) peptide dimer in water as well as at two different experimentally studied concentrations of hydrated ionic liquids (ILs), ethylammonium mesylate (EAM), ethylammonium nitrate (EAN), and triethylammonium mesylate (TEAM), were carried out employing an umbrella sampling method. We used the average Ψ angle of the peptide backbone as the reaction coordinate to observe the conformational changes of a peptide dimer. Secondary structural element values were calculated for the peptide dimer along the reaction coordinate to see the transition of the peptide dimer between β-sheet and α-helix conformations. We observe the β-sheet conformation as the global minimum on the free energy surfaces in both EAM and EAN ILs at both the concentrations and at a low concentration of TEAM. However, we observe α-helix conformation as the global minimum at a high concentration of TEAM. Our results are in good correlation with the experimental findings. We calculated the average number of intramolecular and intermolecular hydrogen bonds of α-helix and β-sheet conformations in all solutions, and they are in correlation with the secondary structure element values. To understand the peptide–IL interactions, atom–atom radial distribution functions of cation, anion, and water around amide oxygen and hydrogen atoms were calculated. The solvent-accessible surface area of the peptide dimer was calculated to understand the exposure of the peptide towards the solvent during conformational changes. Finally, van der Waals (vdW) and Coulomb interaction energies were calculated between peptide–cation, peptide–anion, and peptide–water to understand the stability of conformations in different concentrations. We find that the TEA cation has more vdW interaction energy compared to Coulomb interaction energy with peptide in 70% (w/w) TEAM, which mimics a membrane-like environment to induce α-helix conformation rather than β-sheet conformation.

Molecular dynamics simulations of amyloid-β (16–22) peptide dimer at two different experimentally studied concentrations of hydrated ethylammonium mesylate, ethylammonium nitrate, and triethylammonium mesylate were carried out employing an umbrella sampling method.