Binding of Aβ peptide creates lipid density depression in DMPC bilayer

A Small Molecule Impedes the Aβ1-42 Tetramer Neurotoxicity by Preserving Membrane Integrity: Microsecond Multiscale Simulations

Amyloid-β (Aβ1-42) peptides aggregated into plaques deposited in the brain are the main hallmark of Alzheimer's disease (AD), a social and economic burden worldwide. In this context, insoluble Aβ1-42 fibrils are the main components of plaques. The recent trials that used approved AD drugs show that they can remove the fibrils from AD patients' brains, but they did not halt the course of the disease. Mounting evidence envisages that the soluble Aβ1-42 oligomers' interactions with the neuronal membrane trigger higher cell death than Aβ1-42 fibril interactions. Developing a compound that can alleviate the oligomer's toxicity is one of the most demanding tasks for curing the disease. We performed two molecular dynamics (MD) simulations in an explicit solvent model. In the first case, 55-μs of multiscale all-atom (AA)/coarse-grained (CG) MD simulations were carried out to decipher the impact of a previously described small anti-Aβ molecule, termed M30 (2-octahydroisoquinolin-2(1H)-ylethanamine), on an Aβ1-42 tetramer structure in close contact with a DMPC bilayer. In the second case, 15-μs AA/CG MD simulations were performed to rationalize the dynamics between Aβ1-42 and Aβ1-42-M30 tetramer complexes embedded in DMPC. On the membrane bilayer, we found that the Aβ1-42 tetramer penetrates the bilayer surface due to unrestricted conformational flexibility and many contacts with the membrane phosphate groups. In contrast, no Aβ1-42-M30 tetramer penetration was observed during the entire course of the simulation. In the case of the membrane-embedded Aβ1-42 tetramer, the integrity of the bottom bilayer leaflet was severely affected by the interactions between the negatively charged phosphate groups and the positively charged residues of the Aβ1-42 tetramer, resulting in a deep tetramer penetration into the bilayer hydrophobic region. These contacts were not observed in the case of the membrane-embedded Aβ1-42-M30 tetramer. It was noted that M30 molecules bind to Aβ1-42 tetramer through hydrogen bonds, resulting in a conformational stable Aβ1-42-M30 complex. The associated complex has reduced conformational changes and an enhanced rigidity that prevents the tetramer dissociation by interfering with the tetramer-membrane contacts. Our findings suggest that the M30 molecules could bind to Aβ1-42 tetramer resulting in a rigid structure, and that such complexes do not significantly perturb the membrane bilayer organization. These observations support the in vitro and in vivo experimental evidence that the M30 molecules prevent synaptotocity, improving AD-affected mice memory.

Binding and dimerization of PGLa peptides in anionic lipid bilayer studied by replica exchange molecular dynamics

The 21-residue PGLa peptide is well known for antimicrobial activity attributed to its ability to compromize bacterial membranes. Using all-atom explicit solvent replica exchange molecular dynamics with solute tempering, we studied PGLa binding to a model anionic DMPC/DMPG bilayer at the high peptide:lipid ratio that promotes PGLa dimerization (a two peptides per leaflet system). As a reference we used our previous simulations at the low peptide:lipid ratio (a one peptide per leaflet system). We found that the increase in the peptide:lipid ratio suppresses PGLa helical propensity, tilts the bound peptide toward the bilayer hydrophobic core, and forces it deeper into the bilayer. Surprisingly, at the high peptide:lipid ratio PGLa binding induces weaker bilayer thinning, but deeper water permeation. We explain these effects by the cross-correlations between lipid shells surrounding PGLa that leads to a much diminished efflux of DMPC lipids from the peptide proximity at the high peptide:lipid ratio. Consistent with the experimental data the propensity for PGLa dimerization was found to be weak resulting in coexistence of monomers and dimers with distinctive properties. PGLa dimers assemble via apolar criss-cross interface and become partially expelled from the bilayer residing at the bilayer-water boundary. We rationalize their properties by the dimer tendency to preserve favorable electrostatic interactions between lysine and phosphate lipid groups as well as to avoid electrostatic repulsion between lysines in the low dielectric environment of the bilayer core. PGLa homedimer interface is predicted to be distinct from that involved in PGLa-magainin heterodimers.

Replica Exchange with Hybrid Tempering Efficiently Samples PGLa Peptide Binding to Anionic Bilayer

We evaluated the utility of a variant of the replica exchange method, a replica exchange with hybrid tempering (REHT), for all-atom explicit water biomolecular simulations and compared it with a more traditional replica exchange with the solute tempering (REST) algorithm. As a test system, we selected a 21-mer antimicrobial peptide PGLa binding to an anionic DMPC/DMPG lipid bilayer. Application of REHT revealed the following binding mechanism. Due to the strong hydrophobic moment, the bound PGLa adopts an extensive helical structure. The binding free energy landscape identifies two major bound states, a metastable surface bound state and a dominant inserted state. In both states, positively charged PGLa amino acids maintain electrostatic interactions with anionic phosphate groups by rotating the PGLa helix around its axis. PGLa binding causes an influx of anionic DMPG and an efflux of zwitterionic DMPC lipids from the peptide proximity. PGLa thins the bilayer and disorders the adjacent fatty acid tails. Deep invasion of water wires into the bilayer hydrophobic core is detected in the inserted peptide state. The analysis of charge density distributions indicated that peptide positive charges are nearly compensated for by lipid negative charges and water dipole ordering, whereas ions play no role in peptide binding. Thus, electrostatic interactions are the key energetic factor in binding cationic PGLa to an anionic DMPC/DMPG bilayer. Comparison of REHT and REST shows that due to exclusion of lipids from tempered partition, REST lags behind REHT in peptide equilibration, particularly, with respect to peptide insertion and helix acquisition. As a result, REST struggles to provide accurate details of PGLa binding, although it still qualitatively maps the bimodal binding mechanism. Importantly, REHT not only equilibrates PGLa in the bilayer faster than REST, but also with less computational effort. We conclude that REHT is a preferable choice for studying interfacial biomolecular systems.

Binding Models of Aβ42 Peptide with Membranes Explored by Molecular Simulations

One of the factors contributing to the toxicity of amyloid-β (Aβ) peptides is the destruction of membrane integrity through Aβ peptide-membrane interactions. The binding of Aβ peptides to membranes has been studied by experiments and theoretical simulations extensively. The exact binding mechanism, however, still remains elusive. In the present study, the molecular basis of the peptide-bilayer binding mechanism of the full-length Aβ42 monomer with POPC/POPS/CHOL bilayers is investigated by all-atom (AA) simulations. Three main binding models in coil, bend, and turn structures are obtained. Model 1 of the three models with the central hydrophobic core (CHC) buried inside the membrane is the dominant binding model. The structural features of the peptide, the peptide-bilayer interacting regions, the intrapeptide interactions, and peptide-water interactions are studied. The binding of the Aβ42 monomer to the POPC/POPS/CHOL bilayer is also explored by coarse-grained (CG) simulations as a complement. Both the AA and CG simulations show that residues in CHC prefer forming interactions with the bilayer, indicating the crucial role of CHC in peptide-bilayer binding. Our results can provide new insights for the investigation of the peptide-bilayer binding mechanism of the Aβ peptide.

De Novo Transmembrane Aggregation of Aβ10-40 Peptides in an Anionic Lipid Bilayer

Using the all-atom model and 10 μs serial replica-exchange molecular dynamics (SREMD), we investigated the binding of Alzheimer's Aβ10-40 peptides to the anionic dimyristoylphosphatidylcholine/dimyristoylphosphatidylglycerol (DMPC/DMPG) lipid bilayer. Our objective was to probe de novo transmembrane Aβ10-40 aggregation and to test the utility of SREMD. Our results are threefold. First, upon binding, Aβ10-40 adopts a helical structure in the C-terminus and deeply inserts into the bilayer. Binding is primarily controlled by electrostatic interactions of the peptides with water, ions, and lipids, particularly, anionic DMPG. Second, Aβ-bilayer interactions reorganize lipids in the proximity of the bound peptides, causing an influx of DMPG lipids into the Aβ binding footprint. Third and most important, computed free energy landscapes reveal that Aβ10-40 peptides partition into monomeric and dimeric species. The dimers result from transmembrane aggregation of the peptides and induce a striking lipid density void throughout both leaflets in the bilayer. There are multiple factors stabilizing transmembrane dimers, including van der Waals and steric interactions, electrostatic interactions, and hydrogen bonding, hydration, and entropic gains originating from dimer conformations and lipid disorder. We argue that helix dipole-dipole interactions underestimated in the all-atom force field must be a contributing factor to stabilizing antiparallel transmembrane dimers. We propose that transmembrane aggregates serve as mechanistic links between the populations of extra- and intracellular Aβ peptides. From the computational perspective, SREMD is found to be a viable alternative to traditional replica-exchange simulations.

Calcium inhibits penetration of Alzheimer's Aβ1 -42 monomers into the membrane

Calcium ion regulation plays a crucial role in maintaining neuronal functions such as neurotransmitter release and synaptic plasticity. Copper (Cu²⁺ ) coordination to amyloid-β (Aβ) has accelerated Aβ_1-42 aggregation that can trigger calcium dysregulation by enhancing the influx of calcium ions by extensive perturbing integrity of the membranes. Aβ_1-42 aggregation, calcium dysregulation, and membrane damage are Alzheimer disease (AD) implications. To gain a detail of calcium ions' role in the full-length Aβ_1-42 and Aβ_1-42 -Cu²⁺ monomers contact, the cellular membrane before their aggregation to elucidate the neurotoxicity mechanism, we carried out 2.5 μs extensive molecular dynamics simulation (MD) to rigorous explorations of the intriguing feature of the Aβ_1-42 and Aβ_1-42 -Cu²⁺ interaction with the dimyristoylphosphatidylcholine (DMPC) bilayer in the presence of calcium ions. The outcome of the results compared to the same simulations without calcium ions. We surprisingly noted robust binding energies between the Aβ_1-42 and membrane observed in simulations containing without calcium ions and is two and a half fold lesser in the simulation with calcium ions. Therefore, in the case of the absence of calcium ions, N-terminal residues of Aβ_1-42 deeply penetrate from the surface to the center of the bilayer; in contrast to calcium ions presence, the N- and C-terminal residues are involved only in surface contacts through binding phosphate moieties. On the other hand, Aβ_1-42 -Cu²⁺ actively participated in surface bilayer contacts in the absence of calcium ions. These contacts are prevented by forming a calcium bridge between Aβ_1-42 -Cu²⁺ and the DMPC bilayer in the case of calcium ions presence. In a nutshell, Calcium ions do not allow Aβ_1-42 penetration into the membranes nor contact of Aβ_1-42 -Cu²⁺ with the membranes. These pieces of information imply that the calcium ions mediate the membrane perturbation via the monomer interactions but do not damage the membrane; they agree with the western blot experimental results of a higher concentration of calcium ions inhibit the membrane pore formation by Aβ peptides.

Mechanisms of Binding of Antimicrobial Peptide PGLa to DMPC/DMPG Membrane

PGLa belongs to a class of antimicrobial peptides showing strong affinity to anionic bacterial membranes. Using all-atom explicit solvent replica exchange molecular dynamics with solute tempering, we studied binding of PGLa to a model anionic dimyristoylphosphatidylcholine/dimyristoylphosphatidylglycerol (DMPC/DMPG) bilayer. Due to a strong hydrophobic moment, PGLa upon binding adopts a helical structure and two distinct bound states separated by a significant free energy barrier. In these states, the C-terminus helix is either surface bound or inserted into the bilayer, whereas the N-terminus remains anchored in the bilayer. Analysis of the free energy landscape indicates that the transition between the two states involves a C-terminus helix rotation permitting the peptide to preserve the interactions between cationic Lys amino acids and anionic lipid phosphorus groups. We calculated the free energy of PGLa binding and showed that it is mostly governed by the balance between desolvation of PGLa positive charges and formation of electrostatic PGLa-lipid interactions. PGLa binding induces minor bilayer thinning but causes pronounced lipid redistribution resulting from an influx of DMPG lipids into the binding footprint and efflux of DMPC lipids. Our in silico results rationalize the S-state detected in NMR experiments.

Partitioning of Benzoic Acid into 1,2-Dimyristoyl-sn-glycero-3-phosphocholine and Blood-Brain Barrier Mimetic Bilayers

Using an all-atom explicit water model and replica exchange umbrella sampling simulations, we investigated the molecular mechanisms of benzoic acid partitioning into two model lipid bilayers. The first was formed of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) lipids, whereas the second was composed of an equimolar mixture of DMPC, 1,2-dimyristoyl-sn-glycero-3-phosphoethanolamine, palmitoylsphingomyelin, and cholesterol to constitute a blood-brain barrier (BBB) mimetic bilayer. Comparative analysis of benzoic acid partitioning into the two bilayers has revealed qualitative similarities. Partitioning into the DMPC and BBB bilayers is thermodynamically favorable although insertion into the former lowers the free energy of benzoic acid by approximately an additional 1 kcal mol^-1. The partitioning energetics for the two bilayers is also largely similar based on the balance of benzoic acid interactions with apolar fatty acid tails, polar lipid headgroups, and water. In both bilayers, benzoic acid retains a considerable number of residual water molecules until reaching the bilayer midplane where it experiences nearly complete dehydration. Upon insertion into the bilayers, benzoic acid undergoes several rotations primarily determined by the interactions with the lipid headgroups. Nonetheless, in addition to the depth of the free energy minimum, the BBB bilayer differs from the DMPC counterpart by a much deeper location of the free energy minimum and the appearance of a high free energy barrier and positioning of benzoic acid near the midplane. Furthermore, DMPC and BBB bilayers exhibit different structural responses to benzoic acid insertion. Taken together, the BBB mimetic bilayer is preferable for an accurate description of benzoic acid partitioning.

Label-free, chronological and selective detection of aggregation and fibrillization of amyloid β protein in serum by microcantilever sensor immobilizing cholesterol-incorporated liposome

To facilitate the early diagnosis of Alzheimer's disease and mild cognitive impairment patients, we developed a cantilever-based microsensor that immobilized liposomes of various phospholipids to detect a trace amount of amyloid β (Aβ) protein, and investigated its aggregation and fibrillization on model cell membranes in human serum. Three species of liposomes composed of different phospholipids of 1,2-dipalmtoyl-sn-glycero-3-phosphocholine (DPPC), DPPC/phosphatidyl ethanolamine and 1,2-dipalmitoyl-sn-glycero-3-phosphorylglycerol having varied hydrophilic groups were applied, which showed different chronological interactions with Aβ(1-40) protein and varied sensitivities of the cantilever sensor, depending on their specific electrostatic charged conditions, hydrophilicity, and membrane fluidity. 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) having short hydrophobic carbon chains confirmed to show a large interaction with Aβ(1-40) and a high sensitivity. Furthermore, the incorporation of cholesterol into DMPC was effective to selectively detect Aβ(1-40) in human serum, which effect was also checked by quartz crystal microbalance. Finally, Aβ detection of 100-pM order was expected selectively in the serum by using the developed biosensor.

The monomers, oligomers, and fibrils of amyloid-β inhibit the activity of mitoBKCa channels by a membrane-mediated mechanism

A causative agent of Alzheimer's disease (AD) is a short amphipathic peptide called amyloid beta (Aβ). Aβ monomers undergo structural changes leading to their oligomerization or fibrillization. The monomers as well as all aggregated forms of Aβ, i.e., oligomers, and fibrils, can bind to biological membranes, thereby modulating membrane mechanical properties. It is also known that some isoforms of the large-conductance calcium-activated potassium (BK_Ca) channel, including the mitochondrial BK_Ca (mitoBK_Ca) channel, respond to mechanical changes in the membrane. Here, using the patch-clamp technique, we investigated the impact of full-length Aβ (Aβ_1-42) and its fragment, Aβ_25-35, on the activity of mitoBK_Ca channels. We found that all forms of Aβ inhibited the activity of the mitoBK_Ca channel in a concentration-dependent manner. Since monomers, oligomers, and fibrils of Aβ exhibit different molecular characteristics and structures, we hypothesized that the inhibition was not due to direct peptide-protein interactions but rather to membrane-binding of the Aβ peptides. Our findings supported this hypothesis by showing that Aβ peptides block mitoBK_Ca channels irrespective of the side of the membrane to which they are applied. In addition, we found that the enantiomeric peptide, D-Aβ_1-42, demonstrated similar inhibitory activity towards mitoBK_Ca channels. As a result, we proposed a general model in which all Aβ forms i.e., monomers, oligomers, and amyloid fibrils, contribute to the progression of AD by exerting a modulatory effect on mechanosensitive membrane components.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Methionine Oxidation Changes the Mechanism of Aβ Peptide Binding to the DMPC Bilayer

Using all-atom explicit solvent replica exchange molecular dynamics simulations with solute tempering, we study the effect of methionine oxidation on Aβ10–40 peptide binding to the zwitterionic DMPC bilayer. By comparing oxidized and reduced peptides, we identified changes in the binding mechanism caused by this modification. First, Met35 oxidation unravels C-terminal helix in the bound peptides. Second, oxidation destabilizes intrapeptide interactions and expands bound peptides. We explain these outcomes by the loss of amphiphilic character of the C-terminal helix due to oxidation. Third, oxidation “polarizes” Aβ binding to the DMPC bilayer by strengthening the interactions of the C-terminus with lipids while largely releasing the rest of the peptide from bilayer. Fourth, in contrast to the wild-type peptide, oxidized Aβ induces significantly smaller bilayer thinning and drop in lipid density within the binding footprint. These observations are the consequence of mixing oxidized peptide amino acids with lipids promoted by enhanced Aβ conformational fluctuations. Fifth, methionine oxidation reduces the affinity of Aβ binding to the DMPC bilayer by disrupting favorable intrapeptide interactions upon binding, which offset the gains from better hydration. Reduced binding affinity of the oxidized Aβ may represent the molecular basis for its reduced cytotoxicity.

Cholesterol Changes the Mechanisms of Aβ Peptide Binding to the DMPC Bilayer

Using isobaric-isothermal all-atom replica-exchange molecular dynamics (REMD) simulations, we investigated the equilibrium binding of Aβ10-40 monomers to the zwitterionic dimyristoylphosphatidylcholine (DMPC) bilayer containing cholesterol. Our previous REMD simulations, which studied binding of the same peptide to the cholesterol-free DMPC bilayer, served as a control, against which we measured the impact of cholesterol. Our findings are as follows. First, addition of cholesterol to the DMPC bilayer partially expels the Aβ peptide from the hydrophobic core and promotes its binding to bilayer polar headgroups. Using thermodynamic and energetics analyses, we argued that Aβ partial expulsion is not related to cholesterol-induced changes in lateral pressure within the bilayer but is caused by binding energetics, which favors Aβ binding to the surface of the densely packed cholesterol-rich bilayer. Second, cholesterol has a protective effect on the DMPC bilayer structure against perturbations caused by Aβ binding. More specifically, cholesterol reduces bilayer thinning and overall depletion of bilayer density beneath the Aβ binding footprint. Third, we found that the Aβ peptide contains a single cholesterol binding site, which involves hydrophobic C-terminal amino acids (Ile31-Val36), Phe19, and Phe20 from the central hydrophobic cluster, and cationic Lys28 from the turn region. This binding site accounts for about 76% of all Aβ-cholesterol interactions. Because cholesterol binding site in the Aβ10-40 peptide does not contain the GXXXG motif featured in cholesterol interactions with the transmembrane domain C99 of the β-amyloid precursor protein, we argued that the binding mechanisms for Aβ and C99 are distinct reflecting their different conformations and positions in the lipid bilayer. Fourth, cholesterol sharply reduces the helical propensity in the bound Aβ peptide. As a result, cholesterol largely eliminates the emergence of helical structure observed upon Aβ transition from a water environment to the cholesterol-free DMPC bilayer. We explain this effect by the formation of hydrogen bonds between cholesterol and the Aβ backbone, which prevent helix formation. Taken together, we expect that our simulations will advance understanding of a molecular-level mechanism behind the role of cholesterol in Alzheimer's disease.

Reduced Lipid Bilayer Thickness Regulates the Aggregation and Cytotoxicity of Amyloid-β

The aggregation of amyloid-β (Aβ) on lipid bilayers has been implicated as a mechanism by which Aβ exerts its toxicity in Alzheimer's disease (AD). Lipid bilayer thinning has been observed during both oxidative stress and protein aggregation in AD, but whether these pathological modifications of the bilayer correlate with Aβ misfolding is unclear. Here, we studied peptide-lipid interactions in synthetic bilayers of the short-chain lipid dilauroyl phosphatidylcholine (DLPC) as a simplified model for diseased bilayers to determine their impact on Aβ aggregate, protofibril, and fibril formation. Aβ aggregation and fibril formation in membranes composed of dioleoyl phosphatidylcholine (DOPC) or 1- palmitoyl-2-oleoyl phosphatidylcholine mimicking normal bilayers served as controls. Differences in aggregate formation and stability were monitored by a combination of thioflavin-T fluorescence, circular dichroism, atomic force microscopy, transmission electron microscopy, and NMR. Despite the ability of all three lipid bilayers to catalyze aggregation, DLPC accelerates aggregation at much lower concentrations and prevents the fibrillation of Aβ at low micromolar concentrations. DLPC stabilized globular, membrane-associated oligomers, which could disrupt the bilayer integrity. DLPC bilayers also remodeled preformed amyloid fibrils into a pseudo-unfolded, molten globule state, which resembled on-pathway, protofibrillar aggregates. Whereas the stabilized, membrane-associated oligomers were found to be nontoxic, the remodeled species displayed toxicity similar to that of conventionally prepared aggregates. These results provide mechanistic insights into the roles that pathologically thin bilayers may play in Aβ aggregation on neuronal bilayers, and pathological lipid oxidation may contribute to Aβ misfolding.

Is the Conformational Ensemble of Alzheimer's Aβ10-40 Peptide Force Field Dependent?

By applying REMD simulations we have performed comparative analysis of the conformational ensembles of amino-truncated Aβ10-40 peptide produced with five force fields, which combine four protein parameterizations (CHARMM36, CHARMM22*, CHARMM22/cmap, and OPLS-AA) and two water models (standard and modified TIP3P). Aβ10-40 conformations were analyzed by computing secondary structure, backbone fluctuations, tertiary interactions, and radius of gyration. We have also calculated Aβ10-40 ³J_HNHα-coupling and RDC constants and compared them with their experimental counterparts obtained for the full-length Aβ1-40 peptide. Our study led us to several conclusions. First, all force fields predict that Aβ adopts unfolded structure dominated by turn and random coil conformations. Second, specific TIP3P water model does not dramatically affect secondary or tertiary Aβ10-40 structure, albeit standard TIP3P model favors slightly more compact states. Third, although the secondary structures observed in CHARMM36 and CHARMM22/cmap simulations are qualitatively similar, their tertiary interactions show little consistency. Fourth, two force fields, OPLS-AA and CHARMM22* have unique features setting them apart from CHARMM36 or CHARMM22/cmap. OPLS-AA reveals moderate β-structure propensity coupled with extensive, but weak long-range tertiary interactions leading to Aβ collapsed conformations. CHARMM22* exhibits moderate helix propensity and generates multiple exceptionally stable long- and short-range interactions. Our investigation suggests that among all force fields CHARMM22* differs the most from CHARMM36. Fifth, the analysis of ³J_HNHα-coupling and RDC constants based on CHARMM36 force field with standard TIP3P model led us to an unexpected finding that in silico Aβ10-40 and experimental Aβ1-40 constants are generally in better agreement than these quantities computed and measured for identical peptides, such as Aβ1-40 or Aβ1-42. This observation suggests that the differences in the conformational ensembles of Aβ10-40 and Aβ1-40 are small and the former can be used as proxy of the full-length peptide. Based on this argument, we concluded that CHARMM36 force field with standard TIP3P model produces the most accurate representation of Aβ10-40 conformational ensemble.

Dependence of protein conformational ensembles on force field parameterizations limits the predictive power of molecular dynamics simulations. To address this problem, we evaluated five all-atom force fields for their consistency in reproducing the conformational ensemble of Alzheimer’s Aβ10-40 peptide. To generate conformational ensembles, we have used replica exchange molecular dynamics and computed Aβ10-40 secondary and tertiary structures. We found that, although all force fields predict Aβ10-40 unfolded structure, they strongly disagree on helix and β propensities and tertiary structure distributions. We have also calculated Aβ10-40 J-coupling and residual dipolar coupling constants and compared them with the experimental data for the full-length Aβ1-40 peptide. Unexpectedly, we determined that in silico Aβ10-40 and experimental Aβ1-40 constants are in better agreement than these quantities computed and measured previously for identical peptides, such as Aβ1-40 or Aβ1-42. We then concluded that the conformational ensembles of Aβ10-40 and Aβ1-40 are similar and on this basis argue that CHARMM36 force field with standard TIP3P water model provides the most accurate description of Aβ10-40. Although our objective was not to evaluate the biomolecular force fields in general, our study is expected to facilitate their proper selection for the simulations of Alzheimer’s peptides.

Inclusion of lipopeptides into the DMPC lipid bilayers prevents Aβ peptide insertion

Lipopeptides prevent penetration of Alzheimer's Aβ peptides into lipid bilayers.

Greedy replica exchange algorithm for heterogeneous computing grids

Replica exchange molecular dynamics (REMD) has become a valuable tool in studying complex biomolecular systems. However, its application on distributed computing grids is limited by the heterogeneity of this environment. In this study, we propose a REMD implementation referred to as greedy REMD (gREMD) suitable for computations on heterogeneous grids. To decentralize replica management, gREMD utilizes a precomputed schedule of exchange attempts between temperatures. Our comparison of gREMD against standard REMD suggests four main conclusions. First, gREMD accelerates grid REMD simulations by as much as 40 %. Second, gREMD increases CPU utilization rates in grid REMD by up to 60 %. Third, we argue that gREMD is expected to maintain approximately constant CPU utilization rates and simulation wall-clock times with the increase in the number of replicas. Finally, we show that gREMD correctly implements the REMD algorithm and reproduces the conformational ensemble of a short peptide sampled in our previous standard REMD simulations. We believe that gREMD can find its place in large-scale REMD simulations on heterogeneous computing grids.

Calcium enhances binding of Aβ monomer to DMPC lipid bilayer

Using isobaric-isothermal replica-exchange molecular dynamics and the all-atom explicit-solvent model, we studied the equilibrium binding of Aβ monomers to a zwitterionic dimyristoylphosphatidylcholine (DMPC) bilayer coincubated with calcium ions. Using our previous replica-exchange molecular dynamics calcium-free simulations as a control, we reached three conclusions. First, calcium ions change the tertiary structure of the bound Aβ monomer by destabilizing several long-range intrapeptide interactions, particularly the salt bridge Asp(23)-Lys(28). Second, calcium strengthens Aβ peptide binding to the DMPC bilayer by enhancing electrostatic interactions between charged amino acids and lipid polar headgroups. As a result, Aβ monomer penetrates deeper into the bilayer, making disorder in proximal lipids and bilayer thinning more pronounced. Third, because calcium ions demonstrate strong affinity to negatively charged amino acids, a considerable influx of calcium into the area proximal to the bound Aβ monomer is observed. Consequently, the localizations of negatively charged amino acids and calcium ions in the Aβ binding footprint overlap. Based on our data, we propose a mechanism by which calcium ions strengthen Aβ-bilayer interactions. This mechanism involves two factors: 1) calcium ions make the DMPC bilayer partially cationic and thus attractive to the anionic Aβ peptide; and 2) destabilization of the Asp(23)-Lys(28) salt bridge makes Lys(28) available for interactions with the bilayer. Finally, we conclude that a single Aβ monomer does not promote permeation of calcium ions through the zwitterionic bilayer.