Effect of Ionic Aqueous Environments on the Structure and Dynamics of the Aβ21–30 Fragment: A Molecular-Dynamics Study

Effects of ion type and concentration on the structure and aggregation of the amyloid peptide A β 16 - 22 $$ {\boldsymbol{\beta}}_{16-22} $$

Among the various factors controlling the amyloid aggregation process, the influences of ions on the aggregation rate and the resulting structures are important aspects to consider, which can be studied by molecular simulations. There is a wide variety of protein force fields and ion models, raising the question of which model to use in such studies. To address this question, we perform molecular dynamics simulations of Aβ16-22 , a fragment of the Alzheimer's amyloid β peptide, using different protein force fields, AMBER99SB-disp (A99-d) and CHARMM36m (C36m), and different ion parameters. The influences of NaCl and CaCl2 at various concentrations are studied and compared with the systems without the addition of ions. Our results indicate a sensitivity of the peptide-ion interactions to the different ion models. In particular, we observe a strong binding of Ca2+ to residue E22 with C36m and also with the Åqvist ion model used together with A99-d, which slightly affects the monomeric Aβ16-22 structures and the aggregation rate, but significantly affects the oligomer structures formed in the aggregation simulations. For example, at high Ca2+ concentrations, there was a switch from an antiparallel to a parallel β-sheet. Such ionic influences are of biological relevance because local ion concentrations can change in vivo and could help explain the polymorphism of amyloid fibrils.

Ins and Outs of Rocker Switch Mechanism in Major Facilitator Superfamily of Transporters

The major facilitator superfamily (MFS) of transporters consists of three classes of membrane transporters: symporters, uniporters, and antiporters. Despite such diverse functions, MFS transporters are believed to undergo similar conformational changes within their distinct transport cycles, known as the rocker-switch mechanism. While the similarities between conformational changes are noteworthy, the differences are also important since they could potentially explain the distinct functions of symporters, uniporters, and antiporters of the MFS superfamily. We reviewed a variety of experimental and computational structural data on a select number of antiporters, symporters, and uniporters from the MFS family to compare the similarities and differences of the conformational dynamics of three different classes of transporters.

How do salt and lipids affect conformational dynamics of Aβ42 monomers in water?

It is well established that amyloid β-protein (Aβ) self-assembly is involved in triggering of Alzheimer's disease. On the other hand, evidence of physiological function of Aβ interacting with lipids has only begun to emerge. Details of Aβ-lipid interactions, which may underlie physiological and pathological activities of Aβ, are not well understood. Here, the effects of salt and 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) lipids on conformational dynamics of Aβ42 monomer in water are examined by all-atom molecular dynamics (MD). We acquired six sets of 250 ns long MD trajectories for each of the three lipid concentrations (0, 27, and 109 mM) in the absence and presence of 150 mM salt. Ten replica trajectories per set are used to enhance sampling of Aβ42 conformational space. We show that salt facilitates long-range tertiary contacts in Aβ42, resulting in more compact Aβ42 conformations. By contrast, addition of lipids results in lipid-concentration dependent Aβ42 unfolding concomitant with enhanced stability of the turn in the A21-A30 region. At the high lipid concentration, salt enables the N-terminal region of Aβ42 to form long-range tertiary contacts and interact with lipids, which results in formation of a parallel β-strand. Aβ42 forms stable lipid-protein complexes whereby the protein is adhered to the lipid cluster rather than embedded into it. We propose that the inability of Aβ42 monomer to get embedded into the lipid cluster may be important for facilitating repair of leaks in the blood-brain barrier without penetrating and damaging cellular membranes.

Influence of ion specificity and concentration on the conformational transition of intrinsically disordered sheep prion peptide

The structural sensitivity of the IDPs with the ions has been observed experimentally; however, it is still unclear how the presence of different metal ions affects structural stability. We performed atomistic molecular dynamics simulation of sheep prion peptide (142-167) in presence of different monovalent, divalent ions at various concentrations to find out the effect of the size, charge, and ionic concentration on the structure of the peptide. It is found that Li + ions have a higher survival probability compared to Na + , K + and Mg2 + affecting the solvation structure of the protein leading to the alpha-helix structure. At high concentration, due to the increase in the ion-solvent and ion-counter interactions, the effect of the ions is screened on the surface of the protein and hence no ion specificity is observed. This study demonstrates how ions can be used to regulate the protein structure and function that can help in designing drugs.

Does the inclusion of electronic polarisability lead to a better modelling of peptide aggregation?

Simulating the process of amyloid aggregation with atomic detail is a challenging task for various reasons. One of them is that it is difficult to parametrise a force field such that all protein states ranging from the folded through the unfolded to the aggregated state are represented with the same level of accuracy. Here, we test whether the consideration of electronic polarisability improves the description of the different states of Aβ_16–22. Surprisingly, the CHARMM Drude polarisable force field is found to perform worse than its unpolarisable counterpart CHARMM36m. Sources for this failure of the Drude model are discussed.

Simulating the process of amyloid aggregation is a hard task. We test whether the inclusion of electronic polarisability as done in CHARMM-Drude improves the modelling of Aβ_16–22 aggregation and find it does not. Reasons for the failure are given.

Destabilisation of the structure of transthyretin is driven by Ca2

Tetrameric transthyretin (TTR) transports thyroid hormones and retinol in plasma and cerebrospinal fluid and performs protective functions under stress conditions. Ageing and mutations result in TTR destabilisation and the formation of the amyloid deposits that dysregulate Ca²⁺ homeostasis. Our aim was to determine whether Ca²⁺ affects the structural stability of TTR. We show, using multiple techniques, that Ca²⁺ does not induce prevalent TTR dissociation and/or oligomerisation. However, in the presence of Ca²⁺, TTR exhibits altered conformational flexibility and different interactions with the solvent molecules. These structural changes lead to the formation of the sub-populations of non-native TTR conformers and to the destabilisation of the structure of TTR. Moreover, the sub-population of TTR molecules undergoes fragmentation that is augmented by Ca²⁺. We postulate that Ca²⁺ constitutes the structural and functional switch between the native and non-native forms of TTR, and therefore tip the balance towards age-dependent pathological calcification.

Symmetry-breaking transitions in the early steps of protein self-assembly

Protein misfolding and subsequent self-association are complex, intertwined processes, resulting in development of a heterogeneous population of aggregates closely related to many chronic pathological conditions including Type 2 Diabetes Mellitus and Alzheimer's disease. To address this issue, here, we develop a theoretical model in the general framework of linear stability analysis. According to this model, self-assemblies of peptides with pronounced conformational flexibility may become, under particular conditions, unstable and spontaneously evolve toward an alternating array of partially ordered and disordered monomers. The predictions of the theory were verified by atomistic molecular dynamics (MD) simulations of islet amyloid polypeptide (IAPP) used as a paradigm of aggregation-prone polypeptides (proteins). Simulations of dimeric, tetrameric, and hexameric human-IAPP self-assemblies at physiological electrolyte concentration reveal an alternating distribution of the smallest domains (of the order of the peptide mean length) formed by partially ordered (mainly β-strands) and disordered (turns and coil) arrays. Periodicity disappears upon weakening of the inter-peptide binding, a result in line with the predictions of the theory. To further probe the general validity of our hypothesis, we extended the simulations to other peptides, the Aβ(1-40) amyloid peptide, and the ovine prion peptide as well as to other proteins (SOD1 dimer) that do not belong to the broad class of intrinsically disordered proteins. In all cases, the oligomeric aggregates show an alternate distribution of partially ordered and disordered monomers. We also carried out Surface Enhanced Raman Scattering (SERS) measurements of hIAPP as an experimental validation of both the theory and in silico simulations.

Current Status of AMOEBA-IL: A Multipolar/Polarizable Force Field for Ionic Liquids

Computational simulations of ionic liquid solutions have become a useful tool to investigate various physical, chemical and catalytic properties of systems involving these solvents. Classical molecular dynamics and hybrid quantum mechanical/molecular mechanical (QM/MM) calculations of IL systems have provided significant insights at the atomic level. Here, we present a review of the development and application of the multipolar and polarizable force field AMOEBA for ionic liquid systems, termed AMOEBA–IL. The parametrization approach for AMOEBA–IL relies on the reproduction of total quantum mechanical (QM) intermolecular interaction energies and QM energy decomposition analysis. This approach has been used to develop parameters for imidazolium– and pyrrolidinium–based ILs coupled with various inorganic anions. AMOEBA–IL has been used to investigate and predict the properties of a variety of systems including neat ILs and IL mixtures, water exchange reactions on lanthanide ions in IL mixtures, IL–based liquid–liquid extraction, and effects of ILs on an aniline protection reaction.

Interactions of a multifunctional di-triazole derivative with Alzheimer's Aβ42monomer and Aβ42protofibril: a systematic molecular dynamics study

The molecular dynamics simulations results highlighted that the multi-target-directed ligand6nstabilizes the native α-helix conformation of the Aβ₄₂monomer and induces a sizable destabilization in the Aβ₄₂protofibril structure.

Alkali ion influence on structure and stability of fibrillar amyloid-β oligomers

Alzheimer's disease is characterized by the aggregation of Amyloid-β (Aβ) peptide into oligomers, fibrils and plaques. Many factors influencing this process as well as the stability of the various Aβ aggregates are known to date, and include the concentration and type of metal ions. Most experimental and theoretical studies have concentrated on heavy metal ions, like Fe²⁺, Zn²⁺, or Cu²⁺, while the smaller alkali ions Li⁺, Na⁺, and K⁺ have not gained much attention notwithstanding their role and ubiquity in physiological environments. In this work, we applied atomistic molecular dynamics simulations to investigate the potential role of these alkali ions in stabilizing fibrillar Aβ oligomers of different size and topology, i.e., single and double filament systems comprising 3-24 peptide chains per filament. We find a pronounced difference on the molecular level in the interaction behavior with free carboxylate groups of the Aβ oligomer: Li⁺ forms stable bridged interactions, whereas K⁺ interacts more transiently and lacks bridging. The behavior of Na⁺ is in between, so that this ion-protein interaction obeys the renowned Hofmeister series. These differences are also reflected in the ability of the alkali ions to stabilize the oligomer secondary structure. The stabilizing effect is most pronounced for the smaller fibrillar oligomers, suggesting that the type of alkali ion critically affects the initial stages of fibril formation. Our findings thus offer a molecular explanation for the observation that the polymorphisms of Aβ fibril structures are caused by differences in the surrounding ionic environment. Graphical abstract Influence of alkali ions on the structure and stability of fibrillar amyloid-β oligomers.

Salt Bridge in Aqueous Solution: Strong Structural Motifs but Weak Enthalpic Effect

Salt bridges are elementary motifs of protein secondary and tertiary structure and are commonly associated with structural driving force that increases stability. Often found on the interface to the solvent, they are highly susceptible to solvent–solute interactions, primarily with water but also with other cosolvents (especially ions). We have investigated the interplay of an Arginine–Aspartic acid salt bridge with simple salt ions in aqueous solution by means of molecular dynamics simulations. Besides structural and dynamical features at equilibrium, we have computed the mean force along the dissociation pathway of the salt bridge. We demonstrate that solvated ions influence the behavior of the salt bridge in a very specific and local way, namely the formation of tight ionic pairs Li⁺/Na⁺–Asp⁻. Moreover, our findings show that the enthalpic relevance of the salt bridge is minor, regardless of the presence of solvated ions.

Force-Field Induced Bias in the Structure of Aβ21-30: A Comparison of OPLS, AMBER, CHARMM, and GROMOS Force Fields

In this work we examine the dynamics of an intrinsically disordered protein fragment of the amyloid β, the Aβ21-30, under seven commonly used molecular dynamics force fields (OPLS-AA, CHARMM27-CMAP, AMBER99, AMBER99SB, AMBER99SB-ILDN, AMBER03, and GROMOS53A6), and three water models (TIP3P, TIP4P, and SPC/E). We find that the tested force fields and water models have little effect on the measures of radii of gyration and solvent accessible surface area (SASA); however, secondary structure measures and intrapeptide hydrogen-bonding are significantly modified, with AMBER (99, 99SB, 99SB-ILDN, and 03) and CHARMM22/27 force-fields readily increasing helical content and the variety of intrapeptide hydrogen bonds. On the basis of a comparison between the population of helical and β structures found in experiments, our data suggest that force fields that suppress the formation of helical structure might be a better choice to model the Aβ21-30 peptide.

Spontaneous dimer states of the Aβ21–30decapeptide

Computational examination of the spontaneous dimerization of Aβ_21–30and stability measures of the resulting parallel and anti-parallel aligned dimers.