Characterizing amyloid‐beta protein misfolding from molecular dynamics simulations with explicit water

Extracellular deposition of amyloid‐beta (Aβ) protein, a fragment of membrane glycoprotein called β‐amyloid precursor transmembrane protein (βAPP), is the major characteristic for the Alzheimer's disease (AD). However, the structural and mechanistic information of forming Aβ protein aggregates in a lag phase in cell exterior has been still limited. Here, we have performed multiple all‐atom molecular dynamics simulations for physiological 42‐residue amyloid‐beta protein (Aβ42) in explicit water to characterize most plausible aggregation‐prone structure (APS) for the monomer and the very early conformational transitions for Aβ42 protein misfolding process in a lag phase. Monitoring the early sequential conformational transitions of Aβ42 misfolding in water, the APS for Aβ42 monomer is characterized by the observed correlation between the nonlocal backbone H‐bond formation and the hydrophobic side‐chain exposure. Characteristics on the nature of the APS of Aβ42 allow us to provide new insight into the higher aggregation propensity of Aβ42 over Aβ40, which is in agreement with the experiments. On the basis of the structural features of APS, we propose a plausible aggregation mechanism from APS of Aβ42 to form fibril. The structural and mechanistic observations based on these simulations agree with the recent NMR experiments and provide the driving force and structural origin for the Aβ42 aggregation process to cause AD. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011     

A Hexameric Peptide Barrel as Building Block of Amyloid‐β Protofibrils

Oligomeric and protofibrillar aggregates formed by the amyloid‐β peptide (Aβ) are believed to be involved in the pathology of Alzheimer’s disease. Central to Alzheimer pathology is also the fact that the longer Aβ₄₂ peptide is more prone to aggregation than the more prevalent Aβ₄₀. Detailed structural studies of Aβ oligomers and protofibrils have been impeded by aggregate heterogeneity and instability. We previously engineered a variant of Aβ that forms stable protofibrils and here we use solid‐state NMR spectroscopy and molecular modeling to derive a structural model of these. NMR data are consistent with packing of residues 16 to 42 of Aβ protomers into hexameric barrel‐like oligomers within the protofibril. The core of the oligomers consists of all residues of the central and C‐terminal hydrophobic regions of Aβ, and hairpin loops extend from the core. The model accounts for why Aβ₄₂ forms oligomers and protofibrils more easily than Aβ₄₀.     

An Original Monomer Sampling from a Ready-Made Aβ42 NMR Fibril Suggests a Turn-β-Strand Synergetic Seeding Mechanism

Aggregation of amyloid beta (Aβ) is a central step of Alzheimer's disease. Aβ42 monomers are building blocks in the formation of both “on pathway” intermediate structures and “off pathway” oligomers. How to sample an Aβ monomer becomes a problem however due to the instinct of Aβ high flexibility and diversity as well as aggregation propensity. Currently, (1) most samplings focus on either the ready-made helix-rich 1Z0Q/1IYT NMR structure, or the completely extended conformation, but (2) few on a ready-made Aβ NMR fibril (i. e., 2BEG). Here we compare the simulation results from sampling in scheme (1) with that in scheme (2), and find that the coil and β-sheet contents in the 1Z0Q-sampled system are comparable to the counterparts in the 2BEG-sampled system, but with a large difference in simulation time and dynamics character. 1Z0Q-sampled system not only takes several times longer than the 2BEG-sampled one, and only β1-seeding dynamics characteristic is observed probably due to far insufficient conformation transition in the limited simulation time. Two dynamics characteristics of Aβ42 folding observed experimentally, that either β1 region or β2 region aggregates first, reproduce in the present simulations for 2BEG-sampled system however, suggesting a preferential sampling in the future simulation. In addition, a turn-β-strand synergetic seeding mechanism of aggregation is first proposed based on the trajectory analyses on the four regions of Aβ42 chain.     

Exploring the inter-molecular interactions in amyloid-β protofibril with molecular dynamics simulations and molecular mechanics Poisson-Boltzmann surface area free energy calculations

Aggregation of amyloid-β (Aβ) peptides correlates with the pathology of Alzheimer's disease. However, the inter-molecular interactions between Aβ protofibril remain elusive. Herein, molecular mechanics Poisson-Boltzmann surface area analysis based on all-atom molecular dynamics simulations was performed to study the inter-molecular interactions in Aβ(17-42) protofibril. It is found that the nonpolar interactions are the important forces to stabilize the Aβ(17-42) protofibril, while electrostatic interactions play a minor role. Through free energy decomposition, 18 residues of the Aβ(17-42) are identified to provide interaction energy lower than -2.5 kcal/mol. The nonpolar interactions are mainly provided by the main chain of the peptide and the side chains of nine hydrophobic residues (Leu17, Phe19, Phe20, Leu32, Leu34, Met35, Val36, Val40, and Ile41). However, the electrostatic interactions are mainly supplied by the main chains of six hydrophobic residues (Phe19, Phe20, Val24, Met35, Val36, and Val40) and the side chains of the charged residues (Glu22, Asp23, and Lys28). In the electrostatic interactions, the overwhelming majority of hydrogen bonds involve the main chains of Aβ as well as the guanidinium group of the charged side chain of Lys28. The work has thus elucidated the molecular mechanism of the inter-molecular interactions between Aβ monomers in Aβ(17-42) protofibril, and the findings are considered critical for exploring effective agents for the inhibition of Aβ aggregation.     

Why Is the C-terminus of Abeta(1-42) more unfolded than that of Abeta(1-40)? Clues from hydrophobic interaction

Abeta(1-40) and Abeta(1-42) are the main forms of amyloid beta (Abeta) peptides in the brain of Alzheimer's patients; however, the latter possesses much stronger aggregation and deposition propensity than the former, which is partially attributed to the more unfolded C-terminus of Abeta(1-42) than that of Abeta(1-40). To explore the physical basis underlying the different dynamic behaviors of both Abeta peptides, parallel molecular dynamics (MD) simulations on Abeta(1-40) and Abeta(1-42) were performed to investigate their thermal unfolding processes. It is revealed that the addition of residues 41 and 42 in Abeta(1-42) disrupts the C-terminal hydrophobic core, which triggers the unraveling of the C-terminal helix of Abeta(1-42). This conclusion is supported by the MD simulation on the I41A mutant of Abeta(1-42), in which the C-terminal helix possesses relatively higher conformational stability than that of wild type Abeta(1-42) owing to the change in hydrophobic interaction patterns.     

Monomer Dynamics of Alzheimer Peptides and Kinetic Control of Early Aggregation in Alzheimer's Disease

The rate of reconfiguration—or intramolecular diffusion—of monomeric Alzheimer (Aβ) peptides is measured and, under conditions that aggregation is more likely, peptide diffusion slows down significantly, which allows bimolecular associations to be initiated. By using the method of Trp–Cys contact quenching, the rate of reconfiguration is observed to be about five times faster for Aβ₄₀, which aggregates slowly, than that for Aβ₄₂, which aggregates quickly. Furthermore, the rate of reconfiguration for Aβ₄₂ speeds up at higher pH, which slows aggregation, and in the presence of the aggregation inhibitor curcumin. The measured reconfiguration rates are able to predict the early aggregation behavior of the Aβ peptide and provide a kinetic basis for why Aβ₄₂ is more prone to aggregation than Aβ₄₀, despite a difference of only two amino acids.     

Structural and Material Properties of Amyloid Aβ40/42 Fibrils

In this study, structural and mechanical properties of a series of models of Aβ₄₂ (one‐ and two‐fold) and Aβ₄₀ (two‐ and three‐fold) fibrils have been computed by using all‐atom molecular dynamics simulations. Based on calculations of the twist angle (θ) and periodicity (v=360d/θ), oligomers formed by 20, 11, and 13 monomers were found to be the smallest realistic models of three‐fold Aβ₄₀, one‐fold Aβ₄₂, and two‐fold Aβ₄₂ fibrils, respectively. Our results predict that the Aβ₄₀ fibrils initially exist in two staggered conformations [STAG(+2) and STAG(+1)] and then undergo a [STAG(+2)→STAG(+1)] transformation in a size‐dependent manner. The length of the loop region consisting of the residues 23–29 shrinks with the elongation of both Aβ₄₀ and Aβ₄₂ fibrils. A comparison of the computed potential energy suggests that a two‐fold Aβ₄₀ aggregate is more stable than its three‐fold counterpart, and that Aβ₄₂ oligomers can exist only in one‐fold conformation for aggregates of more than 11 monomers in length. The computed Young′s modulus and yield strengths of 50 GPa and 0.95 GPa, respectively, show that these aggregates possess excellent material properties.     

Formation and stabilization model of the 42-mer Abeta radical: implications for the long-lasting oxidative stress in Alzheimer's disease

Amyloid fibrils mainly consist of 40-mer and 42-mer peptides (Abeta40, Abeta42). Abeta42 is believed to play a crucial role in the pathogenesis of Alzheimer's disease because its aggregative ability and neurotoxicity are considerably greater than those of Abeta40. The neurotoxicity of Abeta peptides involving the generation of free radicals is closely related to the S-oxidized radical cation of Met-35. However, the cation's origin and mechanism of stabilization remain unclear. Recently, structural models of fibrillar Abeta42 and Abeta40 based on systematic proline replacement have been proposed by our group [Morimoto, A.; et al. J. Biol. Chem. 2004, 279, 52781] and Wetzel's group [Williams, A. D.; et al. J. Mol. Biol. 2004, 335, 833], respectively. A major difference between these models is that our model of Abeta42 has a C-terminal beta-sheet region. Our biophysical study on Abeta42 using electron spin resonance (ESR) suggests that the S-oxidized radical cation of Met-35 could be generated by the reduction of the tyrosyl radical at Tyr-10 through a turn structure at positions 22 and 23, and stabilized by a C-terminal carboxylate anion through an intramolecular beta-sheet at positions 35-37 and 40-42 to form a C-terminal core that would lead to aggregation. A time-course analysis of the generation of radicals using ESR suggests that stabilization of the radicals by aggregation might be a main reason for the long-lasting oxidative stress of Abeta42. In contrast, the S-oxidized radical cation of Abeta40 is too short-lived to induce potent neurotoxicity because no such stabilization of radicals occurs in Abeta40.     

Taylor dispersion analysis discloses the impairment of Aβ peptide aggregation by the presence of a fluorescent tag

The use of fluorescently tagged amyloid peptides, implicated in Alzheimer's disease, to study their aggregation at low concentrations is a common method; however, the fluorescent tag should not introduce a bias in the aggregation process. In this work, native amyloid peptides Aβ(1–40) and Aβ(1–42) and fluorescein-5-isothiocyanate (FITC), tagged ones, were studied using Taylor dispersion analysis coupled with a simultaneous UV and light-emitting diode-induced fluorescence detection, to unravel the effect of FITC on the aggregation process. For that, a total concentration of 100 µM of peptides consisting of a mixture of native and tagged ones (up to 10% in moles) was applied. Results demonstrated that FITC had a strong inhibition effect upon the aggregation behaviour of Aβ(1–42), whereas for Aβ(1–40), only a retardation in kinetics was observed. It was also shown that when mixed solutions of Aβ(1–40) and Aβ(1–42) are used, the Aβ(1–42) alloform was the leading peptide in the aggregation process, and when the latter was tagged, the aggregation kinetics decreased but the lifetime of potentially toxic oligomers was drastically increased. These results confirmed that the hydrophilicity of the N-terminus part of the peptide plays a major role in the aggregation process.     

Z_(Aβ3)和Aβ_(16–40)亲和作用的分子机理解析

淀粉样多肽(amyloid-βpeptide,Aβ)聚集是引起阿尔兹海默症(Alzheimer's disease,AD)的主要原因。开发Aβ聚集抑制剂是治疗AD的最有效手段之一。利用噬菌体展示技术筛选出来的Z_(Aβ3)蛋白质能够有效抑制Aβ聚集,但Z_(Aβ3)和Aβ之间的作用区域和关键氨基酸残基尚不清楚。针对此问题,本研究利用分子动力学模拟、MM-PBSA自由能计算和分解方法研究了Z_(Aβ3)-Aβ_(16–40)复合物之间的相互作用机制。结果表明,Z_(Aβ3)的β-股和Aβ_(16–40)之间的亲和作用占主导,而Z_(Aβ3)的α-螺旋贡献很小。利用分子力学-帕松波尔茨曼溶剂可及化表面积方法(MM-PBSA)自由能分解发现Z_(Aβ3)的热点残基为E15、I16、V17、Y18、L19、P20、N21和L22,而Aβ_(16–40)的热点残基为F19、F20、A21、E22、D23、K28、I31、I32、G33、L34、M35、V36、G38和V40。Z_(Aβ3)通过将发夹型Aβ单体包埋在α-螺旋围成的疏水性腔体内来阻碍Aβ聚集。这种结合模式为设计高效的Aβ蛋白质类抑制剂提供了三个基本要素:高亲和性的结合片段(β-股)、附属结构(α-螺旋)和通过二硫键形成的稳定构象。高亲和性结合片段能竞争性地与Aβ单体结合,附属结构α-螺旋可以阻碍其它Aβ单体靠近,而稳定的构象是上述两种要素发挥作用的基础,三者协同作用可以有效地抑制Aβ聚集。     

Solid‐Phase Synthesis and Characterization of N‐Terminally Elongated Aβ−3–x‐Peptides

In addition to the prototypic amyloid‐β (Aβ) peptides Aβ_1–40 and Aβ_1–42, several Aβ variants differing in their amino and carboxy termini have been described. Synthetic availability of an Aβ variant is often the key to study its role under physiological or pathological conditions. Herein, we report a protocol for the efficient solid‐phase peptide synthesis of the N‐terminally elongated Aβ‐peptides Aβ_?3–38, Aβ_?3–40, and Aβ_?3–42. Biophysical characterization by NMR spectroscopy, CD spectroscopy, an aggregation assay, and electron microscopy revealed that all three peptides were prone to aggregation into amyloid fibrils. Immunoprecipitation, followed by mass spectrometry, indicated that Aβ_?3–38 and Aβ_?3–40 are generated by transfected cells even in the presence of a tripartite β‐site amyloid precursor protein cleaving enzyme 1 (BACE1) inhibitor. The elongated Aβ peptides starting at Val(?3) can be separated from N‐terminally‐truncated Aβ forms by high‐resolution isoelectric‐focusing techniques, despite virtually identical isoelectric points. The synthetic Aβ variants and the methods presented here are providing tools to advance our understanding of the potential roles of N‐terminally elongated Aβ variants in Alzheimer's disease.     

Residue‐Specific Dynamics and Local Environmental Changes in Aβ40 Oligomer and Fibril Formation

Elucidating local dynamics of protein aggregation is crucial for understanding the mechanistic details of protein amyloidogenesis. Herein, we studied the residue‐specific dynamics and local environmental changes of Aβ40 along the course of aggregation by using para‐cyanophenylalanine (Phe_CN) as a fluorescent and vibrational probe. Our results show that the Phe_CN residues introduced at various positions all exhibited an immediate decay of fluorescence intensity, indicating a relatively synergistic process in early oligomer formation. The fast decreases in the fluorescence intensities of residues 19 and 20 in the central hydrophobic core region and residue 10 in the N‐terminal region suggest that they play crucial roles in the formation of the oligomeric core. The Phe_CN4 residue exhibits a remarkably slower decrease in fluorescence intensity, implicating its dynamic conformational characteristics in oligomer and fibril formation. Our results also suggest that the N‐terminal residues in fibrils are surrounded by a relatively hydrophobic local environment, as opposed to being solvated.     

Structure of gramicidin a in a lipid bilayer environment determined using molecular dynamics simulations and solid-state NMR data

Two different high-resolution structures recently have been proposed for the membrane-spanning gramicidin A channel: one based on solid-state NMR experiments in oriented phospholipid bilayers (Ketchem, R. R.; Roux, B.; Cross, T. A. Structure 1997, 5, 1655-1669; Protein Data Bank, PDB:1MAG); and one based on two-dimensional NMR in detergent micelles (Townsley, L. E.; Tucker, W. A.; Sham, S.; Hinton, J. F. Biochemistry 2001, 40, 11676-11686; PDB:1JNO). Despite overall agreement, the two structures differ in peptide backbone pitch and the orientation of several side chains; in particular that of the Trp at position 9. Given the importance of the peptide backbone and Trp side chains for ion permeation, we undertook an investigation of the two structures using molecular dynamics simulation with an explicit lipid bilayer membrane, similar to the system used for the solid-state NMR experiments. Based on 0.1 micros of simulation, both backbone structures converge to a structure with 6.25 residues per turn, in agreement with X-ray scattering, and broad agreement with SS backbone NMR observables. The side chain of Trp 9 is mobile, more so than Trp 11, 13, and 15, and undergoes spontaneous transitions between the orientations in 1JNO and 1MAG. Based on empirical fitting to the NMR results, and umbrella sampling calculations, we conclude that Trp 9 spends 80% of the time in the 1JNO orientation and 20% in the 1MAG orientation. These results underscore the utility of molecular dynamics simulations in the analysis and interpretation of structural information from solid-state NMR.     

All-atom molecular dynamics studies of the full-length β-amyloid peptides

β-Amyloid peptides are believed to play an essential role in Alzheimer’s disease (AD), due to their sedimentation in the form of β-amyloid aggregates in the brain of AD-patients, and the in vitro neurotoxicity of oligomeric aggregates. The monomeric peptides come in different lengths of 39–43 residues, of which the 42 alloform seems to be most strongly associated with AD-symptoms. Structural information on these peptides to date comes from NMR studies in acidic solutions, organic solvents, or on shorter fragments of the peptide. In addition X-ray and solid-state NMR investigations of amyloid fibrils yield insight into the structure of the final aggregate and therefore define the endpoint of any conformational change of an Aβ-monomer along the aggregation process. The conformational changes necessary to connect the experimentally known conformations are not yet understood and this process is an active field of research.     

Polymeric complements to the Alzheimer's disease biomarker β-amyloid isoforms Aβ1-40 and Aβ1-42 for blood serum analysis under denaturing conditions

Treatment of Alzheimer's disease (AD) is plagued by a lack of practical and reliable methods allowing early diagnosis of the disease. We here demonstrate that robust receptors prepared by molecular imprinting successfully address current limitations of biologically derived receptors in displaying affinity for hydrophobic peptide biomarkers for AD under denaturing conditions. C-terminal epitope-imprinted polymers showing enhanced binding affinity for Aβ1-42 were first identified from a 96-polymer combinatorial library. This information was then used to synthesize molecularly imprinted polymers for both of the β-amyloid (Aβ) isoforms and a corresponding nonimprinted polymer. A solid-phase extraction method was developed to be compatible with sample loading under conditions of complete protein denaturation. This resulted in a method capable of quantitatively and selectively enriching a shorter C-terminal peptide corresponding to the sequences Aβ33-40 and Aβ33-42 as well as the full-length sequence Aβ1-40 and Aβ1-42 from a 4 M guanidinum chloride solution. Application of the method to serum allowed selective, high-recovery extraction of both biomarkers at spiking levels marginally higher than clinically relevant concentrations found in cerebrospinal fluid.     

Natural Compounds against Alzheimer’s Disease: Molecular Recognition of Aβ1–42 Peptide by Salvia sclareoides Extract and its Major Component,Rosmarinic Acid,as Investigated by NMR

Amyloid peptides, Aβ1–40 and Aβ1–42, represent major molecular targets to develop potential drugs and diagnostic tools for Alzheimer’s Disease (AD). In fact, oligomeric and fibrillar aggregates generated by these peptides are amongst the principal components of amyloid plaques found post mortem in patients suffering from AD. Rosmarinic acid has been demonstrated to be effective in preventing the aggregation of amyloid peptides in vitro and to delay the progression of the disease in animal models. Nevertheless, no information is available about its molecular mechanism of action. Herein, we report the NMR characterization of the interaction of Salvia sclareoides extract and that of its major component, rosmarinic acid, with Aβ1–42 peptide, whose oligomers have been described as the most toxic Aβ species in vivo. Our data shed light on the structural determinants of rosmarinic acid–Aβ1–42 oligomers interaction, thus allowing the elucidation of its mechanism of action. They also provide important information for the rational design of new compounds with higher affinity for Aβ peptides to generate new anti‐amyloidogenic molecules and/or molecular tools for the specific targeting of amyloid aggregates in vivo. In addition, we identified methyl caffeate, another natural compound present in different plants and human diet, as a good ligand of Aβ1–42 oligomers, which also shows anti‐amyloidogenic activity. Finally, we demonstrated the possibility to exploit STD‐NMR and trNOESY experiments to screen extracts from natural sources for the presence of Aβ peptide ligands.     

突变对嗜热古菌蛋白[P62A]Ssh10b主链构象影响的核磁共振研究

嗜热古菌蛋白Ssh10b突变体[P62A]Ssh10b具有良好的热稳定性. [P62A]Ssh10b的结构测定结果显示, 其α2螺旋上残基K48和D51之间形成了一对盐键, 而且, 突变D51将影响蛋白的热稳定性. 为了探索D51的突变对蛋白热稳定性的影响, 构建了突变体[D51N/P62A]Ssh10b的质粒, 并获得了高纯度的15N和13C双标记[D51N/P62A]Ssh10b. 通过对异核三共振NMR实验数据的解析, 完成了对[D51N/P62A]Ssh10b的主链共振近乎完全的指认. 比较突变以及未突变蛋白质的主链1HN和15N化学位移, 在[P62A]Ssh10b的结构基础上进行分析发现, D51N突变显著地影响了α2螺旋骨架构象, 并进一步影响到古菌Lβ2α2loop区域、β4 的N端区域、Lβ3β4 loop的C端区域以及β3与Lβ3β4 loop的交界区域. 结果表明, 由于D51N突变破坏了α2螺旋上K48和D51之间的盐键, 影响了[D51N/P62A]Ssh10b的α2螺旋构象, 并影响到蛋白的其它相关部位的局部构象, 说明 [P62A]Ssh10b 的高热稳定性可能与其溶液构象密切相关. 为进一步运用NMR研究[D51N/P62A]Ssh10b分子结构特性与耐热机制间的关系奠定了基础.     

Aβ40 and Aβ42 amyloid fibrils exhibit distinct molecular recycling properties

A critical aspect to understanding the molecular basis of Alzheimer's disease (AD) is the characterization of the kinetics of interconversion between the different species present during amyloid-β protein (Aβ) aggregation. By monitoring hydrogen/deuterium exchange in Aβ fibrils using electrospray ionization mass spectrometry, we demonstrate that the Aβ molecules comprising the fibril continuously dissociate and reassociate, resulting in molecular recycling within the fibril population. Investigations on Aβ40 and Aβ42 amyloid fibrils reveal that molecules making up Aβ40 fibrils recycle to a much greater extent than those of Aβ42. By examining factors that could influence molecular recycling and by running simulations, we show that the rate constant for dissociation of molecules from the fibril (k(off)) is much greater for Aβ40 than that for Aβ42. Importantly, the k(off) values obtained for Aβ40 and Aβ42 reveal that recycling occurs on biologically relevant time scales. These results have implications for understanding the role of Aβ fibrils in neurotoxicity and for designing therapeutic strategies against AD.     

Comparative molecular dynamics studies of wild-type and oxidized forms of full-length Alzheimer amyloid beta-peptides Abeta(1-40) and Abeta(1-42)

In this study, all-atom 50 ns molecular dynamics simulations are performed on the full-length amyloid beta (Abeta) monomers (WT-Abeta(1-40) and WT-Abeta(1-42)) and their oxidized forms (Met35(O)-Abeta(1-40) and Met35(O)-Abeta(1-42)) in aqueous solution. The effects of the oxidation state of Met35 and the presence of dipeptide (Ile41-Ala42) on the secondary structures of the three distinct regions (the central hydrophobic core region 17-21 (LVFFA), the loop 23-28 (DVGSNK), and the second hydrophobic domain 29-35 (GAIIGLM)) of all monomers have been analyzed in detail, and results are compared with the available experimental information. Our simulations indicate that the WT-Abeta(1-40) monomer adopts an overall beta-hairpin-like structure, which is promoted by the turn region (24-27). This turn region is stabilized through salt-bridge formation between the Asp23 and Lys28 residues. In contrast, the overall structure of the oxidized (Met35(O)-Abeta(1-40)) monomer can be divided into three well-defined bend regions separated by coil segments. These structural differences may be critical for the measured decrease in the rate of aggregation of Met35(O)-Abeta(1-40) peptide. In the WT-Abeta(1-42) monomer, in comparison to the WT-Abeta(1-40), the Asp23-Lys28 salt bridge is absent, and consequently, the turn in the middle (24-27) region has a smaller curvature. The observed difference in the aggregation rates of these two peptides may be related to the opening of the turn (24-27) stabilized by the Asp23-Lys28 salt bridge. For WT-Abeta(1-42), in the absence of this salt bridge, the unfolding and aggregation events may be more favorable than for WT-Abeta(1-40).     

Combining conformational sampling and selection to identify the binding mode of zinc-bound amyloid peptides with bifunctional molecules

The pathogenesis of Alzheimer's disease (AD) has been suggested to be related with the aggregation of amyloid β (Aβ) peptides. Metal ions (e.g. Cu, Fe, and Zn) are supposed to induce the aggregation of Aβ. Recent development of bifunctional molecules that are capable of interacting with Aβ and chelating biometal ions provides promising therapeutics to AD. However, the molecular mechanism for how Aβ, metal ions, and bifunctional molecules interact with each other is still elusive. In this study, the binding mode of Zn(2+)-bound Aβ with bifunctional molecules was investigated by the combination of conformational sampling of full-length Aβ peptides using replica exchange molecular dynamics simulations (REMD) and conformational selection using molecular docking and classical MD simulations. We demonstrate that Zn(2+)-bound Aβ((1-40)) and Aβ((1-42)) exhibit different conformational ensemble. Both Aβ peptides can adopt various conformations to recognize typical bifunctional molecules with different binding affinities. The bifunctional molecules exhibit their dual functions by first preferentially interfering with hydrophobic residues 17-21 and/or 30-35 of Zn(2+)-bound Aβ. Additional interactions with residues surrounding Zn(2+) could possibly disrupt interactions between Zn(2+) and Aβ, which then facilitate these small molecules to chelate Zn(2+). The binding free energy calculations further demonstrate that the association of Aβ with bifunctional molecules is driven by enthalpy. Our results provide a feasible approach to understand the recognition mechanism of disordered proteins with small molecules, which could be helpful to the design of novel AD drugs.