Does Replica Exchange with Solute Tempering Efficiently Sample Aβ Peptide Conformational Ensembles?

Surface Hydrophobicity Strongly Influences Adsorption and Conformation of Amyloid Beta Derived Peptides

The formation of amyloid fibrils is a common feature of many protein systems. It has implications in both health, as amyloid fibrils are implicated in over 30 degenerative diseases, and in the biological functions of proteins. Surfaces have long been known to affect the formation of fibrils but the specific effect depends on the details of both the surface and protein. Fully understanding the role of surfaces in fibrillization requires microscopic information on protein conformation on surfaces. In this paper replica exchange molecular dynamics simulation is used to investigate the model fibril forming protein, Aβ(10-40) (a 31-residue segment of the amyloid-beta protein) on surfaces of different hydrophobicity. Similar to other proteins Aβ(10-40) is found to adsorb strongly onto hydrophobic surfaces. It also adopts significantly different sets of conformations on hydrophobic and polar surfaces, as well as in bulk solution. On hydrophobic surfaces, it adopts partially helical structures, with the helices overlapping with beta-strand regions in the mature fibril. These may be helical intermediates on the fibril formation pathway, suggesting a mechanism for the enhanced fibril formation seen on hydrophobic surfaces.

Competitive Binding of Viral Nuclear Localization Signal Peptide and Inhibitor Ligands to Importin-α Nuclear Transport Protein

Venezuelan equine encephalitis virus (VEEV) is a highly virulent pathogen whose nuclear localization signal (NLS) sequence from capsid protein binds to the host importin-α transport protein and blocks nuclear import. We studied the molecular mechanisms by which two small ligands, termed I1 and I2, interfere with the binding of VEEV's NLS peptide to importin-α protein. To this end, we performed all-atom replica exchange molecular dynamics simulations probing the competitive binding of the VEEV coreNLS peptide and I1 or I2 ligand to the importin-α major NLS binding site. As a reference, we used our previous simulations, which examined noncompetitive binding of the coreNLS peptide or the inhibitors to importin-α. We found that both inhibitors completely abrogate the native binding of the coreNLS peptide, forcing it to adopt a manifold of nonnative loosely bound poses within the importin-α major NLS binding site. Both inhibitors primarily destabilize the native coreNLS binding by masking its amino acids rather than competing with it for binding to importin-α. Because I2, in contrast to I1, binds off-site localizing on the edge of the major NLS binding site, it inhibits fewer coreNLS native binding interactions than I1. Structural analysis is supported by computations of the free energies of the coreNLS peptide binding to importin-α with or without competition from the inhibitors. Specifically, both inhibitors reduce the free energy gain from coreNLS binding, with I1 causing significantly larger loss than I2. To test our simulations, we performed AlphaScreen experiments measuring IC50 values for both inhibitors. Consistent with in silico results, the IC50 value for I1 was found to be lower than that for I2. We hypothesize that the inhibitory action of I1 and I2 ligands might be specific to the NLS from VEEV's capsid protein.

Dissecting the geometric and hydrophobic constraints of stapled peptides

Stapled peptides are a promising class of molecules with potential as highly specific probes of protein-protein interactions and as therapeutics. Hydrocarbon stapling affects the peptide properties through the interplay of two factors: enhancing the overall hydrophobicity and constraining the conformational flexibility. By constructing a series of virtual peptides, we study the role of each factor in modulating the structural properties of a hydrocarbon-stapled peptide PM2, which has been shown to enter cells, engage its target Mouse Double Minute 2 (MDM2), and activate p53. Hamiltonian replica exchange molecular dynamics (HREMD) simulations suggest that hydrocarbon stapling favors helical populations of PM2 through a combination of the geometric constraints and the enhanced hydrophobicity of the peptide. To further understand the conformational landscape of the stapled peptides along the binding pathway, we performed HREMD simulations by restraining the peptide at different distances from MDM2. When the peptide approaches MDM2, the binding pocket undergoes dehydration which appears to be greater in the presence of the stapled peptide compared with the linear peptide. In the binding pocket, the helicity of the stapled peptide is increased due to the favorable interactions between the peptide residues as well as the staple and the microenvironment of the binding pocket, contributing to enhanced affinity. The dissection of the multifaceted mechanism of hydrocarbon stapling into individual factors not only deepens fundamental understanding of peptide stapling, but also provides guidelines for the design of new stapled peptides.

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 of viral nuclear localization signal peptides to importin-α nuclear transport protein

Using all-atom replica-exchange molecular dynamics simulations, we mapped the mechanisms of binding of the nuclear localization signal (NLS) sequence from Venezuelan equine encephalitis virus (VEEV) capsid protein to importin-α (impα) transport protein. Our objective was to identify the VEEV NLS sequence fragment that confers native, experimentally resolved binding to impα as well as to study associated binding energetics and conformational ensembles. The two selected VEEV NLS peptide fragments, KKPK and KKPKKE, show strikingly different binding mechanisms. The minNLS peptide KKPK binds non-natively and nonspecifically by adopting five diverse conformational clusters with low similarity to the x-ray structure 3VE6 of NLS-impα complex. Despite the prevalence of non-native interactions, the minNLS peptide still largely binds to the impα major NLS binding site. In contrast, the coreNLS peptide KKPKKE binds specifically and natively, adopting a largely homogeneous binding ensemble with a dominant, highly native-like conformational cluster. The coreNLS peptide retains most of native binding interactions, including π-cation contacts and a tryptophan cage. While KKPK binding is governed by a complex multistate free energy landscape featuring transitions between multiple binding poses, the coreNLS peptide free energy map is simple, exhibiting a single dominant native-like bound basin. We argue that the origin of the coreNLS peptide binding specificity is several electrostatic interactions formed by the two C-terminal amino acids, Lys10 and Glu11, with impα. The coreNLS sequence is then sufficient for native binding, but none of the amino acids flanking minNLS, including Lys10 and Glu11, are strictly necessary for the native pose. Our analyses indicate that the VEEV coreNLS sequence is virtually unique among human and viral proteins interacting with impα making it a potential target for VEEV-specific inhibitors.

Can Free Energy Perturbation Simulations Coupled with Replica-Exchange Molecular Dynamics Study Ligands with Distributed Binding Sites?

Free energy perturbation coupled with replica exchange with solute tempering (FEP/REST) offers a rigorous approach to compute relative free energy changes for ligands. To determine the applicability of FEP/REST for the ligands with distributed binding poses, we considered two alchemical transformations involving three putative inhibitors I0, I1, and I2 of the Venezuelan equine encephalitis virus nuclear localization signal sequence binding to the importin-α (impα) transporter protein. I0 → I1 and I0 → I2 transformations, respectively, increase or decrease the polarity of the parent molecule. Our objective was three-fold─(i) to verify FEP/REST technical performance and convergence, (ii) to estimate changes in binding free energy ΔΔG, and (iii) to determine the utility of FEP/REST simulations for conformational binding analysis. Our results are as follows. First, our FEP/REST implementation properly follows FEP/REST formalism and produces converged ΔΔG estimates. Due to ligand inherent unbinding, the better FEP/REST strategy lies in performing multiple independent trajectories rather than extending their length. Second, I0 → I1 and I0 → I2 transformations result in overall minor changes in inhibitor binding free energy, slightly strengthening the affinity of I1 and weakening that of I2. Electrostatic interactions dominate binding interactions, determining the enthalpic changes. The two transformations cause opposite entropic changes, which ultimately govern binding affinities. Importantly, we confirm the validity of FEP/REST free energy estimates by comparing them with our previous REST simulations, directly probing binding of three ligands to impα. Third, we established that FEP/REST simulations can sample binding ensembles of ligands. Thus, FEP/REST can be applied (i) to study the energetics of the ligand binding without defined poses and showing minor differences in affinities |ΔΔG| ≲ 0.5 kcal/mol and (ii) to collect ligand binding conformational ensembles.

Binding of Venezuelan Equine Encephalitis Virus Inhibitors to Importin-α Receptors Explored with All-Atom Replica Exchange Molecular Dynamics

Although Venezuelan equine encephalitis virus (VEEV) is a life-threatening pathogen with a capacity for epidemic outbreaks, there are no FDA-approved VEEV antivirals for humans. VEEV cytotoxicity is partially attributed to the formation of a tetrameric complex between the VEEV capsid protein, the nuclear import proteins importin-α and importin-β, and the nuclear export protein CRM1, which together block trafficking through the nuclear pore complex. Experimental studies have identified small molecules from the CL6662 scaffold as potential inhibitors of the viral nuclear localization signal (NLS) sequence binding to importin-α. However, little is known about the molecular mechanism of CL6662 inhibition. To address this issue, we employed all-atom replica exchange molecular dynamics simulations to probe, in atomistic detail, the binding mechanism of CL6662 ligands to importin-α. Three ligands, including G281-1485 and two congeners with varying hydrophobicities, were considered. We investigated the distribution of ligand binding poses, their locations, and ligand specificities measured by the strength of binding interactions. We found that G281-1485 binds nonspecifically without forming well-defined binding poses throughout the NLS binding site. Binding of the less hydrophobic congener becomes strongly on-target with respect to the NLS binding site but remains nonspecific. However, a more hydrophobic congener is a strongly specific binder and the only ligand out of three to form a well-defined binding pose, while partially overlapping with the NLS binding site. On the basis of free energy estimates, we argue that all three ligands weakly compete with the viral NLS sequence for binding to importin-α in an apparent compromise to preserve host NLS binding. We further show that all-atom replica exchange binding simulations are a viable tool for studying ligands binding nonspecifically without forming well-defined binding poses.

Fuzzy Drug Targets: Disordered Proteins in the Drug-Discovery Realm

Intrinsically disordered proteins (IDPs) and regions (IDRs) form a large part of the eukaryotic proteome. Contrary to the structure-function paradigm, the disordered proteins perform a myriad of functions in vivo. Consequently, they are involved in various disease pathways and are plausible drug targets. Unlike folded proteins, that have a defined structure and well carved out drug-binding pockets that can guide lead molecule selection, the disordered proteins require alternative drug-development methodologies that are based on an acceptable picture of their conformational ensemble. In this review, we discuss various experimental and computational techniques that contribute toward understanding IDP "structure" and describe representative pursuances toward IDP-targeting drug development. We also discuss ideas on developing rational drug design protocols targeting IDPs.

Ingenuity in performing replica permutation: How to order the state labels for improving sampling efficiency

In the replica-permutation method, an advanced version of the replica-exchange method, all combinations of replicas and parameters are considered for parameter permutation, and a list of all the combinations is prepared. Here, we report that the temperature transition probability depends on how the list is created, especially in replica permutation with solute tempering (RPST). We found that the transition probabilities decrease at large replica indices when the combinations are sequentially assigned to the state labels as in the originally proposed list. To solve this problem, we propose to modify the list by randomly assigning the combinations to the state labels. We performed molecular dynamics simulations of amyloid-β(16-22) peptides using RPST with the "randomly assigned" list (RPST-RA) and RPST with the "sequentially assigned" list (RPST-SA). The results show the decreases in the transition probabilities in RPST-SA are eliminated, and the sampling efficiency is improved in RPST-RA.

Lysine Acetylation Changes the Mechanism of Aβ25-35 Peptide Binding and Dimerization in the DMPC Bilayer

The impact of Lys28 acetylation on Alzheimer's Aβ peptide binding to the lipid bilayer has not been previously studied, either experimentally or computationally. To probe this common post-translational modification, we performed all-atom replica exchange molecular dynamics simulations targeting binding and aggregation of acetylated acAβ25-35 peptide within the DMPC bilayer. Using the unmodified Aβ25-35 studied previously as a reference, our results can be summarized as follows. First, Lys28 acetylation strengthens the Aβ25-35 hydrophobic moment and consequently promotes the helical structure across the peptide extending it into the N-terminus. Second, because Lys28 acetylation disrupts electrostatic contact between Lys28 and lipid phosphate groups, it reduces the binding affinity of acAβ25-35 peptides to the DMPC bilayer. Accordingly, although acetylation preserves the bimodal binding featuring a preferred inserted state and a less probable surface bound state, it decreases the stability of the former. Third, acetylation promotes acAβ25-35 aggregation and eliminates monomers as thermodynamically viable species. More importantly, acAβ25-35 retains as the most thermodynamically stable the inserted dimer with unique head-to-tail helical aggregation interface. However, due to enhanced helix structure, this dimer state becomes less stable and is less likely to propagate into higher order aggregates. Thus, acetylation is predicted to facilitate the formation of low-molecular-weight oligomers. Other post-translational modifications, including phosphorylation and oxidation, reduce helical propensity and have divergent impact on aggregation. Consequently, acetylation, when considered in its totality, has distinct consequences on Aβ25-35 binding and aggregation in the lipid bilayer.

ACES: Optimized Alchemically Enhanced Sampling

We present an alchemical enhanced sampling (ACES) method implemented in the GPU-accelerated AMBER free energy MD engine. The methods hinges on the creation of an "enhanced sampling state" by reducing or eliminating selected potential energy terms and interactions that lead to kinetic traps and conformational barriers while maintaining those terms that curtail the need to otherwise sample large volumes of phase space. For example, the enhanced sampling state might involve transforming regions of a ligand and/or protein side chain into a noninteracting "dummy state" with internal electrostatics and torsion angle terms turned off. The enhanced sampling state is connected to a real-state end point through a Hamiltonian replica exchange (HREMD) framework that is facilitated by newly developed alchemical transformation pathways and smoothstep softcore potentials. This creates a counterdiffusion of real and enhanced-sampling states along the HREMD network. The effect of a differential response of the environment to the real and enhanced-sampling states is minimized by leveraging the dual topology framework in AMBER to construct a counterbalancing HREMD network in the opposite alchemical direction with the same (or similar) real and enhanced sampling states at inverted end points. The method has been demonstrated in a series of test cases of increasing complexity where traditional MD, and in several cases alternative REST2-like enhanced sampling methods, are shown to fail. The hydration free energy for acetic acid was shown to be independent of the starting conformation, and the values for four additional edge case molecules from the FreeSolv database were shown to have a significantly closer agreement with experiment using ACES. The method was further able to handle different rotamer states in a Cdk2 ligand identified as fractionally occupied in crystal structures. Finally, ACES was applied to T4-lysozyme and demonstrated that the side chain distribution of V111χ1 could be reliably reproduced for the apo state, bound to p-xylene, and in p-xylene→ benzene transformations. In these cases, the ACES method is shown to be highly robust and superior to a REST2-like enhanced sampling implementation alone.

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.

Insights into the Atomistic Mechanisms of Phosphorylation in Disrupting Liquid-Liquid Phase Separation and Aggregation of the FUS Low-Complexity Domain

Fused in sarcoma (FUS), a nuclear RNA binding protein, can not only undergo liquid-liquid phase separation (LLPS) to form dynamic biomolecular condensates but also aggregate into solid amyloid fibrils which are associated with the pathology of amyotrophic lateral sclerosis and frontotemporal lobar degeneration diseases. Phosphorylation in the FUS low-complexity domain (FUS-LC) inhibits FUS LLPS and aggregation. However, it remains largely elusive what are the underlying atomistic mechanisms of this inhibitory effect and whether phosphorylation can disrupt preformed FUS fibrils, reversing the FUS gel/solid phase toward the liquid phase. Herein, we systematically investigate the impacts of phosphorylation on the conformational ensemble of the FUS_37-97 monomer and dimer and the structure of the FUS_37-97 fibril by performing extensive all-atom molecular dynamics simulations. Our simulations reveal three key findings: (1) phosphorylation shifts the conformations of FUS_37-97 from the β-rich, fibril-competent state toward a helix-rich, fibril-incompetent state; (2) phosphorylation significantly weakens protein-protein interactions and enhances protein-water interactions, which disfavor FUS-LC LLPS as well as aggregation and facilitate the dissolution of the preformed FUS-LC fibril; and (3) the FUS_37-97 peptide displays a high β-strand probability in the region spanning residues 52-67, and phosphorylation at S54 and S61 residues located in this region is crucial for the disruption of LLPS and aggregation of FUS-LC. This study may pave the way for ameliorating phase-separation-related pathologies via site-specific phosphorylation.

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.

Intrinsically disordered proteins: Ensembles at the limits of Anfinsen's dogma

Intrinsically disordered proteins (IDPs) are proteins that lack rigid 3D structure. Hence, they are often misconceived to present a challenge to Anfinsen's dogma. However, IDPs exist as ensembles that sample a quasi-continuum of rapidly interconverting conformations and, as such, may represent proteins at the extreme limit of the Anfinsen postulate. IDPs play important biological roles and are key components of the cellular protein interaction network (PIN). Many IDPs can interconvert between disordered and ordered states as they bind to appropriate partners. Conformational dynamics of IDPs contribute to conformational noise in the cell. Thus, the dysregulation of IDPs contributes to increased noise and "promiscuous" interactions. This leads to PIN rewiring to output an appropriate response underscoring the critical role of IDPs in cellular decision making. Nonetheless, IDPs are not easily tractable experimentally. Furthermore, in the absence of a reference conformation, discerning the energy landscape representation of the weakly funneled IDPs in terms of reaction coordinates is challenging. To understand conformational dynamics in real time and decipher how IDPs recognize multiple binding partners with high specificity, several sophisticated knowledge-based and physics-based in silico sampling techniques have been developed. Here, using specific examples, we highlight recent advances in energy landscape visualization and molecular dynamics simulations to discern conformational dynamics and discuss how the conformational preferences of IDPs modulate their function, especially in phenotypic switching. Finally, we discuss recent progress in identifying small molecules targeting IDPs underscoring the potential therapeutic value of IDPs. Understanding structure and function of IDPs can not only provide new insight on cellular decision making but may also help to refine and extend Anfinsen's structure/function paradigm.

Replica permutation with solute tempering for molecular dynamics simulation and its application to the dimerization of amyloid-β fragments

We propose the replica permutation with solute tempering (RPST) by combining the replica-permutation method (RPM) and the replica exchange with solute tempering (REST). Temperature permutations are performed among more than two replicas in RPM, whereas temperature exchanges are performed between two replicas in the replica-exchange method (REM). The temperature transition in RPM occurs more efficiently than in REM. In REST, only the temperatures of the solute region, the solute temperatures, are exchanged to reduce the number of replicas compared to REM. Therefore, RPST is expected to be an improved method taking advantage of these methods. For comparison, we applied RPST, REST, RPM, and REM to two amyloid-β(16–22) peptides in explicit water. We calculated the transition ratio and the number of tunneling events in the temperature space and the number of dimerization events of amyloid-β(16–22) peptides. The results indicate that, in RPST, the number of replicas necessary for frequent random walks in the temperature and conformational spaces is reduced compared to the other three methods. In addition, we focused on the dimerization process of amyloid-β(16–22) peptides. The RPST simulation with a relatively small number of replicas shows that the two amyloid-β(16–22) peptides form the intermolecular antiparallel β-bridges due to the hydrophilic side-chain contact between Lys and Glu and hydrophobic side-chain contact between Leu, Val, and Phe, which stabilizes the dimer of the peptides.

Effect of Cholesterol Molecules on Aβ1-42 Wild-Type and Mutants Trimers

Alzheimer’s disease displays aggregates of the amyloid-beta (Aβ) peptide in the brain, and there is increasing evidence that cholesterol may contribute to the pathogenesis of the disease. Though many experimental and theoretical studies have focused on the interactions of Aβ oligomers with membrane models containing cholesterol, an understanding of the effect of free cholesterol on small Aβ42 oligomers is not fully established. To address this question, we report on replica exchange with a solute tempering simulation of an Aβ42 trimer with cholesterol and compare it with a previous replica exchange molecular dynamics simulation. We show that the binding hot spots of cholesterol are rather complex, involving hydrophobic residues L17–F20 and L30–M35 with a non-negligible contribution of loop residues D22–K28 and N-terminus residues. We also examine the effects of cholesterol on the trimers of the disease-causing A21G and disease-protective A2T mutations by molecular dynamics simulations. We show that these two mutations moderately impact cholesterol-binding modes. In our REST2 simulations, we find that cholesterol is rarely inserted into aggregates but rather attached as dimers and trimers at the surface of Aβ42 oligomers. We propose that cholesterol acts as a glue to speed up the formation of larger aggregates; this provides a mechanistic link between cholesterol and Alzheimer’s disease.

Computer Simulations Aimed at Exploring Protein Aggregation and Dissociation

Protein aggregation can lead to well-defined structures that are functional, but is also the cause of the death of neuron cells in many neurodegenerative diseases. The complexity of the molecular events involved in the aggregation kinetics of amyloid proteins and the transient and heterogeneous characters of all oligomers prevent high-resolution structural experiments. As a result, computer simulations have been used to determine the atomic structures of amyloid proteins at different association stages as well as to understand fibril dissociation. In this chapter, we first review the current computer simulation methods used for aggregation with some atomistic and coarse-grained results aimed at better characterizing the early formed oligomers and amyloid fibril formation. Then we present the applications of non-equilibrium molecular dynamics simulations to comprehend the dissociation of protein assemblies.

Met35 Oxidation Hinders Aβ25-35 Peptide Aggregation within the Dimyristoylphosphatidylcholine Bilayer

Using all-atom explicit solvent replica exchange molecular dynamics simulations, we studied the aggregation of oxidized (ox) Aβ25-35 peptides into dimers mediated by the zwitterionic dimyristoylphosphatidylcholine (DMPC) lipid bilayer. By comparing oxAβ25-35 aggregation with that observed for reduced and phosphorylated Aβ25-35 peptides, we elucidated plausible impact of post-translational modifications on cytotoxicity of Aβ peptides involved in Alzheimer's disease. We found that Met35 oxidation reduces helical propensity in oxAβ25-35 peptides bound to the lipid bilayer and enhances backbone fluctuations. These factors destabilize the wild-type head-to-tail dimer interface and lower the aggregation propensity. Met35 oxidation diversifies aggregation pathways by adding monomeric species to the bound conformational ensemble. The oxAβ25-35 dimer becomes partially expelled from the DMPC bilayer and as a result inflicts limited disruption to the bilayer structure compared to wild-type Aβ25-35. Interestingly, the effect of Ser26 phosphorylation is largely opposite, as it preserves the wild-type head-to-tail aggregation interface and strengthens, not weakens, aggregation propensity. The differing effects can be attributed to the sequence locations of these post-translational modifications, since in contrast to Ser26 phosphorylation, Met35 oxidation directly affects the wild-type C-terminal aggregation interface. A comparison with experimental data is provided.

Albumin Alters the Conformational Ensemble of Amyloid-β by Promiscuous Interactions: Implications for Amyloid Inhibition

Human serum albumin (HSA) is a key endogenous inhibitor of amyloid-β (Αβ) aggregation. In vitro HSA inhibits Aβ fibrillization and targets multiple species along the aggregation pathway including monomers, oligomers, and protofibrils. Amyloid inhibition by HSA has both pathological implications and therapeutic potential, but the underlying molecular mechanism remains elusive. As a first step towards addressing this complex question, we studied the interactions of an Aβ42 monomer with HSA by molecular dynamics simulations. To adequately sample the conformational space, we adapted the replica exchange with solute tempering (REST2) method to selectively heat the Aβ42 peptide in the absence and presence of HSA. Aβ42 binds to multiple sites on HSA with a preference to domain III and adopts various conformations that all differ from the free state. The β-sheet abundances of H14-E22 and A30-M33 regions are significantly reduced by HSA, so are the β-sheet lengths. HSA shifts the conformational ensemble towards more disordered states and alters the β-sheet association patterns. In particular, the frequent association of Q15-V24 and N27-V36 regions into β-hairpin which is critical for aggregation is impeded. HSA primarily interacts with the latter β-region and the N-terminal charged residues. They form promiscuous interactions characterized by salt bridges at the edge of the peptide-protein interface and hydrophobic cores at the center. Consequently, intrapeptide interactions crucial for β-sheet formation are disrupted. Our work builds the bridge between the modification of Aβ conformational ensemble and amyloid inhibition by HSA. It also illustrates the potential of the REST2 method in studying interactions between intrinsically disordered peptides and globular proteins.

High resolution ensemble description of metamorphic and intrinsically disordered proteins using an efficient hybrid parallel tempering scheme

Mapping free energy landscapes of complex multi-funneled metamorphic proteins and weakly-funneled intrinsically disordered proteins (IDPs) remains challenging. While rare-event sampling molecular dynamics simulations can be useful, they often need to either impose restraints or reweigh the generated data to match experiments. Here, we present a parallel-tempering method that takes advantage of accelerated water dynamics and allows efficient and accurate conformational sampling across a wide variety of proteins. We demonstrate the improved sampling efficiency by benchmarking against standard model systems such as alanine di-peptide, TRP-cage and β-hairpin. The method successfully scales to large metamorphic proteins such as RFA-H and to highly disordered IDPs such as Histatin-5. Across the diverse proteins, the calculated ensemble averages match well with the NMR, SAXS and other biophysical experiments without the need to reweigh. By allowing accurate sampling across different landscapes, the method opens doors for sampling free energy landscape of complex uncharted proteins.

An Implementation of Replica Exchange with Dynamical Scaling for Efficient Large-Scale Simulations

An implementation of the replica exchange with dynamical scaling (REDS) method in the commonly used molecular dynamics program GROMACS is presented. REDS is a replica exchange method that requires fewer replicas than conventional replica exchange while still providing data over a range of temperatures and can be used in either constant volume or constant pressure ensembles. Details for running REDS simulations are given, and an application to the human islet amyloid polypeptide (hIAPP) 11-25 fragment shows that the model efficiently samples conformational space.

In Silico Prediction of the Binding, Folding, Insertion, and Overall Stability of Membrane-Active Peptides

Membrane-active peptides (MAPs) are short-length peptides used for potential biomedical applications in diagnostic imaging of tissues, targeted drug delivery, gene delivery, and antimicrobials and antibiotics. The broad appeal of MAPs is that they are infinitely variable, relatively low cost, and biocompatible. However, experimentally characterizing the specific properties of a MAP or its many variants is a low-resolution and potentially time-consuming endeavor; molecular dynamics (MD) simulations have emerged as an invaluable tool in identifying the biophysical interactions that are fundamental to the function of MAPs. In this chapter, a step-by-step approach to discreetly model the binding, folding, and insertion of a membrane-active peptide to a model lipid bilayer using MD simulations is described. Detailed discussion is devoted to the critical aspects of running these types of simulations: prior knowledge of the system, understanding the strengths and weaknesses of molecular mechanics force fields, proper construction and equilibration of the system, realistically estimating both experimental and computational timescales, and leveraging analysis to make direct comparisons to experimental results as often as possible.

Disclosing the Interaction of Gold Nanoparticles with Aβ(1-40) Monomers through Replica Exchange Molecular Dynamics Simulations

Amyloid-β aggregation is one of the principal causes of amyloidogenic diseases that lead to the loss of neuronal cells and to cognitive impairments. The use of gold nanoparticles treating amyloidogenic diseases is a promising approach, because the chemistry of the gold surface can be tuned in order to have a specific binding, obtaining effective tools to control the aggregation. In this paper, we show, by means of Replica Exchange Solute Tempering Molecular Simulations, how electrostatic interactions drive the absorption of Amyloid-β monomers onto citrates-capped gold nanoparticles. Importantly, upon binding, amyloid monomers show a reduced propensity in forming β-sheets secondary structures that are characteristics of mature amyloid fibrils.

Folding Free-Energy Landscape of α-Synuclein (35-97) Via Replica Exchange Molecular Dynamics

The misfolding and aggregation of α-synuclein (α-syn) in Lewy bodies are implicated in the pathogenesis of various neurodegenerative disorders, such as Parkinson's disease and dementia. The formation of α-syn fibrils is a complex process, involving various intermediates and oligomeric forms. These intermediates establish at an early stage of aggregation and subsequently lead to fibrillation. Determining which conformations are accessible to monomeric α-syn and especially, as shown in a recent work, to the central amino acids from residue 35 to residue 97 (63 residues) is thus crucial to understand the formation of these oligomers. Here, we carry out extensive replica exchange molecular dynamics (total time-18 μs) with an all-atom model and explicit solvent to characterize the free-energy landscape of human α-syn (residue 35 to residue 97). The simulation results lead us to identify two free-energy basins. Clustering analysis for the deepest free-energy minimum reveals a compact structure, with a secondary structure predominantly α-helix, while the shallower minimum corresponds to an elongated conformation, also predominantly α-helix. Furthermore, at physiological temperature, we find that conformational rearrangements happen via helix breaks due to the presence of glycine. We also show that the most likely conformations are characterized by the α-helix structure rather than the β-hairpin structure (for residue 38 to residue 53), in contrast with prior simulation studies using coarse-grained models or an implicit solvent. For higher temperatures, we observe a shift in secondary structure with a decrease in the population of α-helix in favor of random coils, β-bend, and β-turns.

Effect of the air-water interface on the conformation of amyloid beta

It has long been recognized that liquid interfaces, such as the air-water interface (AWI), can enhance the formation of protein fibrils. This makes liquid interfaces attractive templates for fibril formation but fully realizing this requires knowledge of protein behavior at interfaces, which is currently lacking. To address this, molecular dynamics simulation is used to investigate fragments of amyloid beta, a model fibril forming protein, at the air-water interface. At the air-water interface, the enrichment of aggregation-prone helical conformations provides a mechanism for the enhancement of fibrillation at interfaces. The conformational ensemble at the air-water interface was also considerably reduced compared to bulk solution due to the tendency of hydrophobic side chains partitioning into the air restricting the range of conformations. Little overlap between the conformational ensembles at the AWI and in the bulk solution was found, suggesting that AWI induces the formation of a different set of structures compared to bulk solution. The smaller Aβ(16-22) and Aβ(25-35) fragments show an increase in the propensity for an ordered secondary structure at the air-water interface but with a increased propensity for turn over other motifs, illustrating the importance of intra-protein interactions for stabilizing helical and extended conformations.

Structural Analysis of (p)ppGpp Reveals Its Versatile Binding Pattern for Diverse Types of Target Proteins

(p)ppGpp, highly phosphorylated guanosine, are global regulatory nucleotides that modulate several biochemical events in bacterial physiology ranging from core central dogma to various metabolic pathways. Conventionally, (p)ppGpp collectively refers to two nucleotides, ppGpp, and pppGpp in the literature. Initially, (p)ppGpp has been discovered as a transcription regulatory molecule as it binds to RNA polymerase and regulates transcriptional gene regulation. During the past decade, several other target proteins of (p)ppGpp have been discovered and as of now, more than 30 proteins have been reported to be regulated by the binding of these two signaling nucleotides. The regulation of diverse biochemical activities by (p)ppGpp requires fine-tuned molecular interactions with various classes of proteins so that it can moderate varied functions. Here we report a structural dynamics of (p)ppGpp in the unbound state using well-defined computational tools and its interactions with target proteins to understand the differential regulation by (p)ppGpp at the molecular level. We carried out replica exchange molecular dynamics simulation studies to enhance sampling of conformations during (p)ppGpp simulation. The detailed comparative analysis of torsion angle conformation of ribose sugar of unbound (p)ppGpp and bound states of (p)ppGpp was carried out. The structural dynamics shows that two linear phosphate chains provide plasticity to (p)ppGpp nucleotides for the binding to diverse proteins. Moreover, the intermolecular interactions between (p)ppGpp and target proteins were characterized through various physicochemical parameters including, hydrogen bonds, van der Waal’s interactions, aromatic stacking, and side chains of interacting residues of proteins. Surprisingly, we observed that interactions of (p)ppGpp to target protein have a consensus binding pattern for a particular functional class of enzymes. For example, the binding of (p)ppGpp to RNA polymerase is significantly different from the binding of (p)ppGpp to the proteins involved in the ribosome biogenesis pathway. Whereas, (p)ppGpp binding to enzymes involved in nucleotide metabolism facilitates the functional regulation through oligomerization. Analysis of these datasets revealed that guanine base-specific contacts are key determinants to discriminate functional class of protein. Altogether, our studies provide significant information to understand the differential interaction pattern of (p)ppGpp to its target and this information may be useful to design antibacterial compounds based on (p)ppGpp analogs.

Phosphorylation Promotes Aβ25-35 Peptide Aggregation within the DMPC Bilayer

The consequences of phosphorylation of the Aβ25-35 peptide at the position Ser26 on its aggregation have not been examined. To investigate them, we performed all-atom replica exchange simulations probing the binding of phosphorylated Aβ25-35 (pAβ25-35) peptides to the dimyristoyl phosphatidylcholine (DMPC) bilayer and their subsequent aggregation. As a control, we used our previous study of unmodified peptides. We found that phosphorylation moderately reduces the helical propensity in pAβ25-35 and its binding affinity to the DMPC bilayer. Phosphorylation preserves the bimodal binding observed for unmodified Aβ25-35, which features a preferred inserted state and a less probable surface bound state. Phosphorylation also retains the inserted dimer with a head-to-tail helical aggregation interface as the most thermodynamically stable state. Importantly, this post-translation modification strengthens interpeptide interactions by adding a new aggregation "hot spot" created by cross-bridging phosphorylated Ser26 with water, cationic ions, or Lys28. The cross-bridging constitutes the molecular mechanism behind most structural phosphorylation effects. In addition, phosphorylation eliminates pAβ25-35 monomers and diversifies the pool of aggregated species. As a result, it changes the binding and aggregation mechanism by multiplying pathways leading to stable inserted dimers. These findings offer a plausible molecular rationale for experimental observations, including enhanced production of low molecular weight oligomers and cytotoxicity of phosphorylated Aβ peptides.

Aggregation of disease-related peptides

Protein misfolding and aggregation of amyloid proteins is the fundamental cause of more than 20 diseases. Molecular mechanisms of the self-assembly and the formation of the toxic aggregates are still elusive. Computer simulations have been intensively used to study the aggregation of amyloid peptides of various amino acid lengths related to neurodegenerative diseases. We review atomistic and coarse-grained simulations of short amyloid peptides aimed at determining their transient oligomeric structures and the early and late aggregation steps.

Do Cholesterol and Sphingomyelin Change the Mechanism of Aβ25-35 Peptide Binding to Zwitterionic Bilayer?

Using replica exchange with solute tempering all-atom molecular dynamics, we studied the equilibrium binding of Aβ_25-35 peptide to the ternary bilayer composed of an equimolar mixture of dimyristoylphosphatidylcholine (DMPC), N-palmitoylsphingomyelin (PSM), and cholesterol. Binding of the same peptide to the pure DMPC bilayer served as a control. Due to significant C-terminal hydrophobic moment, binding to the ternary and DMPC bilayers promotes helical structure in the peptide. For both bilayers a polarized binding profile is observed, in which the N-terminus anchors to the bilayer surface, whereas the C-terminus alternates between unbound and inserted states. Both ternary and DMPC bilayers feature two Aβ_25-35 bound states, surface bound, S, and inserted, I, separated by modest free energy barriers. Experimental data are in agreement with our results but indicate that cholesterol impact is Aβ fragment dependent. For Aβ_25-35, we predict that its binding mechanism is independent of the inclusion of PSM and cholesterol into the bilayer.

Role of the N-Terminal Transmembrane Helix Contacts in the Activation of FGFR3

Fibroblast growth factor receptor 3 (FGFR3) is a member of receptor tyrosine kinases, which is involved in skeletal cell growth, differentiation, and migration. FGFR3 transduces biochemical signals from the extracellular ligand-binding domain to the intracellular kinase domain through the conformational changes of the transmembrane (TM) helix dimer. Here, we apply generalized replica exchange with solute tempering method to wild type (WT) and G380R mutant (G380R) of FGFR3. The dimer interface in G380R is different from WT and the simulation results are in good agreement with the solid-state nuclear magnetic resonance (NMR) spectroscopy. TM helices in G380R are extended more than WT, and thereby, G375 in G380R contacts near the N-termini of the TM helix dimer. Considering that both G380R and G375C show the constitutive activation, the formation of the N-terminal contacts of the TM helices can be generally important for the activation mechanism. © 2019 Wiley Periodicals, Inc.

De novo aggregation of Alzheimer's Aβ25-35 peptides in a lipid bilayer

A potential mechanism of cytotoxicity attributed to Alzheimer’s Aβ peptides postulates that their aggregation disrupts membrane structure causing uncontrollable permeation of Ca²⁺ ions. To gain molecular insights into these processes, we have performed all-atom explicit solvent replica exchange with solute tempering molecular dynamics simulations probing aggregation of the naturally occurring Aβ fragment Aβ25-35 within the DMPC lipid bilayer. To compare the impact produced on the bilayer by Aβ25-35 oligomers and monomers, we used as a control our previous simulations, which explored binding of Aβ25-35 monomers to the same bilayer. We found that compared to monomeric species aggregation results in much deeper insertion of Aβ25-35 peptides into the bilayer hydrophobic core causing more pronounced disruption in its structure. Aβ25-35 peptides aggregate by incorporating monomer-like structures with stable C-terminal helix. As a result the Aβ25-35 dimer features unusual helix head-to-tail topology supported by a parallel off-registry interface. Such topology affords further growth of an aggregate by recruiting additional peptides. Free energy landscape reveals that inserted dimers represent the dominant equilibrium state augmented by two metastable states associated with surface bound dimers and inserted monomers. Using the free energy landscape we propose the pathway of Aβ25-35 binding, aggregation, and insertion into the lipid bilayer.

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.

Self-Assembly of Nanocubic Molecular Capsules via Solvent-Guided Formation of Rectangular Blocks

We investigate the mechanism underlying the self-assembly of gear-shaped amphiphilic molecules into a highly ordered nanocubic capsule ("nanocube") in aqueous methanol. Simulation results show that the solvent molecules play a significant role in the assembly process by directing the primitive intermediates to orthogonal/rectangular shapes, thus creating appropriate building blocks for cubic assembly while avoiding off-pathway stacked aggregates. Free-energy analyses reveal that the interplay of the direct intermonomer interaction and the solvent-mediated repulsion between large aromatic cores (via preferential solvation of methanol on hydrophobic surfaces) leads to the strong trend for perpendicular binding of monomers and hence the solvent-guided formation of rectangular blocks. Furthermore, we report the self-assembly simulation of the nanocube using replica exchange with solute tempering and demonstrate that the simulation can predict a highly ordered nanocapsule structure, assembly intermediates, and encapsulated molecules, which helps promote computer-aided design of functional molecular self-assemblies in explicit solvent.

Temperature-induced collapse of a disordered peptide observed by three sampling methods in molecular dynamics simulations

The conformational ensembles of a disordered peptide, polyglutamine Q15, over a wide temperature range were sampled using multiple replicates of conventional molecular dynamics (cMD) simulations as well as two enhanced sampling methods, temperature replica exchange (TREMD) and replica exchange with solute tempering (REST). The radius of gyration, asphericity, secondary structure, and hydrogen bonding patterns were used for the comparison of the sampling methods. Overall, the three sampling methods generated similar conformational ensembles, with progressive collapse at higher temperatures. Although accumulating the longest simulation time (90 μs), cMD at room temperature missed a small subspace that was sampled by both TREMD and REST. This subspace was high in α-helical content and separated from the main conformational space by an energy barrier. REST used less simulation time than TREMD (36 μs versus 42 μs), and this gap is expected to widen significantly for larger disordered proteins. We conclude that REST is the method of choice for conformational sampling of intrinsically disordered proteins.

Flexible selection of the solute region in replica exchange with solute tempering: Application to protein-folding simulations

Replica-exchange molecular dynamics (REMD) and their variants have been widely used in simulations of the biomolecular structure and dynamics. Replica exchange with solute tempering (REST) is one of the methods where temperature of a pre-defined solute molecule is exchanged between replicas, while solvent temperatures in all the replicas are kept constant. REST greatly reduces the number of replicas compared to the temperature REMD, while replicas at low temperatures are often trapped under their conditions, interfering with the conformational sampling. Here, we introduce a new scheme of REST, referred to as generalized REST (gREST), where the solute region is defined as a part of a molecule or a part of the potential energy terms, such as the dihedral-angle energy term or Lennard-Jones energy term. We applied this new method to folding simulations of a β-hairpin (16 residues) and a Trp-cage (20 residues) in explicit water. The protein dihedral-angle energy term is chosen as the solute region in the simulations. gREST reduces the number of replicas necessary for good random walks in the solute-temperature space and covers a wider conformational space compared to the conventional REST2. Considering the general applicability, gREST should become a promising tool for the simulations of protein folding, conformational dynamics, and an in silico drug design.

The length but not the sequence of peptide linker modules exerts the primary influence on the conformations of protein domains in cellulosome multi-enzyme complexes

Cellulosomes are large multi-protein catalysts produced by various anaerobic microorganisms to efficiently degrade plant cell-wall polysaccharides down into simple sugars. X-ray and physicochemical structural characterisations show that cellulosomes are composed of numerous protein domains that are connected by unstructured polypeptide segments, yet the properties and possible roles of these 'linker' peptides are largely unknown. We have performed coarse-grained and all-atom molecular dynamics computer simulations of a number of cellulosomal linkers of different lengths and compositions. Our data demonstrates that the effective stiffness of the linker peptides, as quantified by the equilibrium fluctuations in the end-to-end distances, depends primarily on the length of the linker and less so on the specific amino acid sequence. The presence of excluded volume - provided by the domains that are connected - dampens the motion of the linker residues and reduces the effective stiffness of the linkers. Simultaneously, the presence of the linkers alters the conformations of the protein domains that are connected. We demonstrate that short, stiff linkers induce significant rearrangements in the folded domains of the mini-cellulosome composed of endoglucanase Cel8A in complex with scaffoldin ScafT (Cel8A-ScafT) of Clostridium thermocellum as well as in a two-cohesin system derived from the scaffoldin ScaB of Acetivibrio cellulolyticus. We give experimentally testable predictions on structural changes in protein domains that depend on the length of linkers.

Communication: Self-assembly of a model supramolecular polymer studied by replica exchange with solute tempering

Conventional molecular-dynamics (cMD) simulation has a well-known limitation in accessible time and length scales, and thus various enhanced sampling techniques have been proposed to alleviate the problem. In this paper, we explore the utility of replica exchange with solute tempering (REST) (i.e., a variant of Hamiltonian replica exchange methods) to simulate the self-assembly of a supramolecular polymer in explicit solvent and compare the performance with temperature-based replica exchange MD (T-REMD) as well as cMD. As a test system, we consider a relatively simple all-atom model of supramolecular polymerization (namely, benzene-1,3,5-tricarboxamides in methylcyclohexane solvent). Our results show that both REST and T-REMD are able to predict highly ordered polymer structures with helical H-bonding patterns, in contrast to cMD which completely fails to obtain such a structure for the present model. At the same time, we have also experienced some technical challenge (i.e., aggregation-dispersion transition and the resulting bottleneck for replica traversal), which is illustrated numerically. Since the computational cost of REST scales more moderately than T-REMD, we expect that REST will be useful for studying the self-assembly of larger systems in solution with enhanced rearrangement of monomers.

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.

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

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