Impact of the mutation A21G (Flemish variant) on Alzheimer's beta-amyloid dimers by molecular dynamics simulations

Soluble oligomers of the amyloid beta-protein (Abeta) are linked to Alzheimer's disease. Irrespective of the nature of the nucleus before fibril growth, dimers are essential species in Abeta assembly, but their transient character has precluded, thus far, high-resolution structure determination. We have investigated the effects of the point mutation A21G on Abeta dimers by performing high temperature all-atom molecular dynamics simulations of Abeta(40), Abeta(42), and their Flemish variants (A21G) starting from their fibrillar conformations. Abeta dimers are found in equilibrium between various topologies, and the absence of common structural features shared by the four species makes problematic the design of a unique inhibitor for blocking dimers. We also show that the impact of the point mutation A21G on Abeta structure and dynamics varies from Abeta(40) to Abeta(42). Finally, we provide a possible structural explanation for the reduced aggregation rate of Abeta fibrils containing the Flemish disease-causing mutation.     

Investigating the mechanism of peptide aggregation: insights from mixed monte carlo-molecular dynamics simulations

The early stages of peptide aggregation are currently not accessible by experimental techniques at atomic resolution. In this article, we address this problem through the application of a mixed simulation scheme in which a preliminary coarse-grained Monte Carlo analysis of the free-energy landscape is used to identify representative conformations of the aggregates and subsequent all-atom molecular dynamics simulations are used to analyze in detail possible pathways for the stabilization of oligomers. This protocol was applied to systems consisting of multiple copies of the model peptide GNNQQNY, whose detailed structures in the aggregated state have been recently solved in another study. The analysis of the various trajectories provides dynamical and structural insight into the details of aggregation. In particular, the simulations suggest a hierarchical mechanism characterized by the initial formation of stable parallel β-sheet dimers and identify the formation of the polar zipper motif as a fundamental feature for the stabilization of initial oligomers. Simulation results are consistent with experimentally derived observations and provide an atomically detailed view of the putative initial stages of fibril formation.     

Folding Trp-cage to NMR resolution native structure using a coarse-grained protein model

We develop a coarse-grained protein model with a simplified amino acid interaction potential. Using this model, we perform discrete molecular dynamics folding simulations of a small 20-residue protein--Trp-cage--from a fully extended conformation. We demonstrate the ability of the Trp-cage model to consistently reach conformations within 2-angstroms backbone root-mean-square distance from the corresponding NMR structures. The minimum root-mean-square distance of Trp-cage conformations in simulations can be <1 angstroms. Our findings suggest that, at least in the case of Trp-cage, a detailed all-atom protein model with a molecular mechanics force field is not necessary to reach the native state of a protein. Our results also suggest that the success of folding Trp-cage in our simulations and in the reported all-atom molecular mechanics simulation studies may be mainly due to the special stabilizing features specific to this miniprotein.     

Spontaneous beta-barrel formation: an all-atom Monte Carlo study of Abeta16-22 oligomerization

Using all-atom Monte Carlo simulations with implicit water, combined with a cluster size analysis, we study the aggregation of Abeta(16) (-22), a peptide capable of forming amyloid fibrils. We consider a system of six initially randomly oriented Abeta(16) (-22) peptides, and investigate the thermodynamics and structural properties of aggregates formed by this system. The system is unaggregated without ordered secondary structure at high temperature, and forms beta-sheet rich aggregates at low temperature. At the crossover between these two regimes, we find that clusters of all sizes occur, whereas the beta-strand content is low. In one of several runs, we observe the spontaneous formation of a beta-barrel with six antiparallel strands. The beta-barrel stands out as the by far most long-lived aggregate seen in our simulations.     

Dimer formation enhances structural differences between amyloid β-protein (1-40) and (1-42): an explicit-solvent molecular dynamics study

Amyloid β-protein (Aβ) is central to the pathology of Alzheimer's disease. A 5% difference in the primary structure of the two predominant alloforms, Aβ(1-40) and Aβ(1-42), results in distinct assembly pathways and toxicity properties. Discrete molecular dynamics (DMD) studies of Aβ(1-40) and Aβ(1-42) assembly resulted in alloform-specific oligomer size distributions consistent with experimental findings. Here, a large ensemble of DMD-derived Aβ(1-40) and Aβ(1-42) monomers and dimers was subjected to fully atomistic molecular dynamics (MD) simulations using the OPLS-AA force field combined with two water models, SPCE and TIP3P. The resulting all-atom conformations were slightly larger, less compact, had similar turn and lower β-strand propensities than those predicted by DMD. Fully atomistic Aβ(1-40) and Aβ(1-42) monomers populated qualitatively similar free energy landscapes. In contrast, the free energy landscape of Aβ(1-42) dimers indicated a larger conformational variability in comparison to that of Aβ(1-40) dimers. Aβ(1-42) dimers were characterized by an increased flexibility in the N-terminal region D1-R5 and a larger solvent exposure of charged amino acids relative to Aβ(1-40) dimers. Of the three positively charged amino acids, R5 was the most and K16 the least involved in salt bridge formation. This result was independent of the water model, alloform, and assembly state. Overall, salt bridge propensities increased upon dimer formation. An exception was the salt bridge propensity of K28, which decreased upon formation of Aβ(1-42) dimers and was significantly lower than in Aβ(1-40) dimers. The potential relevance of the three positively charged amino acids in mediating the Aβ oligomer toxicity is discussed in the light of available experimental data.     

The role of DNA twist in the packaging of viral genomes

We performed molecular dynamics simulations of the genome packaging of bacteriophage P4 using two coarse-grained models of DNA. The first model, 1DNA6 (one pseudo-atom per six DNA basepairs), represents DNA as a string of beads, for which DNA torsions are undefined. The second model, 3DNA6 (three pseudo-atoms per six DNA basepairs), represents DNA as a series of base planes with torsions defined by the angles between successive planes. Bacteriophage P4 was packaged with 1DNA6, 3DNA6 in a torsionally relaxed state, and 3DNA6 in a torsionally strained state. We observed good agreement between the packed conformation of 1DNA6 and the packed conformations of 3DNA6. The free energies of packaging were in agreement, as well. Our results suggest that DNA torsions can be omitted from coarse-grained bacteriophage packaging simulations without significantly altering the DNA conformations or free energies of packaging that the simulations predict.     

Folding landscapes of the Alzheimer amyloid-beta(12-28) peptide

The energy landscape for folding of the 12-28 fragment of the Alzheimer amyloid beta (Abeta) peptide is characterized using replica-exchange molecular dynamics simulations with an all-atom peptide model and explicit solvent. At physiological temperatures, the peptide exists mostly as a collapsed random coil, populating a small fraction (less than 10%) of hairpins with a beta-turn at position V18F19, with another 10% of hairpin-like conformations possessing a bend rather than a turn in the central VFFA positions. A small fraction of the populated states, approximately 14%, adopt polyproline II (PPII) conformations. Folding of the structured hairpin states proceeds through the assembly of two locally stable segments, VFFAE and EDVGS. The interactions stabilizing these locally folded structural motifs are in conflict with those stabilizing the global fold of A12-28, a signature of underlying residual frustration in this peptide. At increased temperature, the population of both beta-strand and PPII conformations diminishes in favor of beta-turn and random-coil states. On the basis of the conformational preferences of Abeta 12-28 monomers, two models for the molecular structure of amyloid fibrils formed by this peptide are proposed.     

Structures and dynamics of β-barrel oligomer intermediates of amyloid-beta16-22 aggregation

Accumulating evidence suggests that soluble oligomers are more toxic than final fibrils of amyloid aggregations. Among the mixture of inter-converting intermediates with continuous distribution of sizes and secondary structures, oligomers in the β-barrel conformation – a common class of protein folds with a closed β-sheet – have been postulated as the toxic species with well-defined three-dimensional structures to perform pathological functions. A common mechanism for amyloid toxicity, therefore, implies that all amyloid peptides should be able to form β-barrel oligomers as the aggregation intermediates. Here, we applied all-atom discrete molecular dynamics (DMD) simulations to evaluate the formation of β-barrel oligomers and characterize their structures and dynamics in the aggregation of a seven-residue amyloid peptide, corresponding to the amyloid core of amyloid-β with a sequence of ¹⁶KLVFFAE²² (Aβ16-22). We carried out aggregation simulations with various numbers of peptides to study the size dependence of aggregation dynamics and assembly structures. Consistent with previous computational studies, we observed the formation of β-barrel oligomers in all-atom DMD simulations. Using a network-based approach to automatically identify β-barrel conformations, we systematically characterized β-barrels of various sizes. Our simulations revealed the conformational inter-conversion between β-barrels and double-layer β-sheets due to increased structural strains upon forming a closed β-barrel while maximizing backbone hydrogen bonds. The potential of mean force analysis further characterized the free energy barriers between these two states. The obtained structural and dynamic insights of β-barrel oligomers may help better understand the molecular mechanism of oligomer toxicities and design novel therapeutics targeting the toxic β-barrel oligomers. This article is part of a Special Issue entitled: Protein Aggregation and Misfolding at the Cell Membrane Interface edited by Ayyalusamy Ramamoorthy     

Millisecond molecular dynamics simulations of KRas-dimer formation and interfaces

Ras dimers have been proposed as building blocks for initiating the extracellular signal-regulated kinase (ERK)/mitogen-activated protein kinase (MAPK) cellular signaling pathway. To better examine the structure of possible dimer interfaces, the dynamics of Ras dimerization, and its potential signaling consequences, we performed molecular dynamics simulations totaling 1 ms of sampling, using an all-atom model of two full-length, farnesylated, guanosine triphosphate (GTP)-bound, wild-type KRas4b proteins diffusing on 29%POPS (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-L-serine)-mixed POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) membranes. Our simulations unveil an ensemble of thermodynamically weak KRas dimers spanning multiple conformations. The most stable conformations, having the largest interface areas, involve helix α₂ and a hypervariable region (HVR). Among the dimer conformations, we found that the HVR of each KRas has frequent interactions with various parts of the dimer, thus potentially mediating the dimerization. Some dimer configurations have one KRas G-domain elevated above the lipid bilayer surface by residing on top of the other G-domain, thus likely contributing to the recruitment of cytosolic Raf kinases in the context of a stably formed multi-protein complex. We identified a variant of the α₄-α₅ KRas-dimer interface that is similar to the interfaces obtained with fluorescence resonance energy transfer (FRET) data of HRas on lipid bilayers. Interestingly, we found two arginine fingers, R68 and R149, that directly interact with the beta-phosphate of the GTP bound in KRas, in a manner similar to what is observed in a crystal structure of GAP-HRas complex, which can facilitate the GTP hydrolysis via the arginine finger of GTPase-activating protein (GAP).     

Mapping the conformational dynamics and pathways of spontaneous steric zipper Peptide oligomerization

The process of protein misfolding and self-assembly into various, polymorphic aggregates is associated with a number of important neurodegenerative diseases. Only recently, crystal structures of several short peptides have provided detailed structural insights into -sheet rich aggregates, known as amyloid fibrils. Knowledge about early events of the formation and interconversion of small oligomeric states, an inevitable step in the cascade of peptide self-assembly, however, remains still limited. We employ molecular dynamics simulations in explicit solvent to study the spontaneous aggregation process of steric zipper peptide segments from the tau protein and insulin in atomistic detail. Starting from separated chains with random conformations, we find a rapid formation of structurally heterogeneous, -sheet rich oligomers, emerging from multiple bimolecular association steps and diverse assembly pathways. Furthermore, our study provides evidence that aggregate intermediates as small as dimers can be kinetically trapped and thus affect the structural evolution of larger oligomers. Alternative aggregate structures are found for both peptide sequences in the different independent simulations, some of which feature characteristics of the known steric zipper conformation (e.g., -sheet bilayers with a dry interface). The final aggregates interconvert with topologically distinct oligomeric states exclusively via internal rearrangements. The peptide oligomerization was analyzed through the perspective of a minimal oligomer, i.e., the dimer. Thereby all observed multimeric aggregates can be consistently mapped onto a space of reduced dimensionality. This novel method of conformational mapping reveals heterogeneous association and reorganization dynamics that are governed by the characteristics of peptide sequence and oligomer size.     

Effects of ATP and Actin-Filament Binding on the Dynamics of the Myosin II S1 Domain

Actin and myosin interact with one another to perform a variety of cellular functions. Central to understanding the processive motion of myosin on actin is the characterization of the individual states along the mechanochemical cycle. We present an all-atom molecular dynamics simulation of the myosin II S1 domain in the rigor state interacting with an actin filament. We also study actin-free myosin in both rigor and post-rigor conformations. Using all-atom level and coarse-grained analysis methods, we investigate the effects of myosin binding on actin, and of actin binding on myosin. In particular, we determine the domains of actin and myosin that interact strongly with one another at the actomyosin interface using a highly coarse-grained level of resolution, and we identify a number of salt bridges and hydrogen bonds at the interface of myosin and actin. Applying coarse-grained analysis, we identify differences in myosin states dependent on actin-binding, or ATP binding. Our simulations also indicate that the actin propeller twist-angle and nucleotide cleft-angles are influenced by myosin at the actomyosin interface. The torsional rigidity of the myosin-bound filament is also calculated, and is found to be increased compared to previous simulations of the free filament.     

Pathway complexity of Alzheimer's beta-amyloid Abeta16-22 peptide assembly

Recent studies suggest that both soluble oligomers and insoluble fibrils have toxic effects in cell cultures, raising the interest in determining the first steps of the assembly process. We have determined the aggregation mechanisms of Abeta(16-22) dimer using the activation-relaxation technique and an approximate free energy model. Consistent with the NMR solid-state analysis, the dimer is predicted to prefer an antiparallel beta sheet structure with the expected registry of intermolecular hydrogen bonds. The simulations, however, locate three other antiparallel minima with nonnative beta sheet registries and one parallel beta sheet structure, slightly destabilized with respect to the ground state. This result is significant because it can explain the observed dependency of beta sheet registry on pH conditions. We also find that assembly of Abeta(16-22) into dimers follows multiple routes, but alpha-helical intermediates are not obligatory. This indicates that destabilization of alpha-helical intermediates is unlikely to abolish oligomerization of Abeta peptides.     

Effects of ATP and Actin-Filament Binding on the Dynamics of the Myosin II S1 Domain

Actin and myosin interact with one another to perform a variety of cellular functions. Central to understanding the processive motion of myosin on actin is the characterization of the individual states along the mechanochemical cycle. We present an all-atom molecular dynamics simulation of the myosin II S1 domain in the rigor state interacting with an actin filament. We also study actin-free myosin in both rigor and post-rigor conformations. Using all-atom level and coarse-grained analysis methods, we investigate the effects of myosin binding on actin, and of actin binding on myosin. In particular, we determine the domains of actin and myosin that interact strongly with one another at the actomyosin interface using a highly coarse-grained level of resolution, and we identify a number of salt bridges and hydrogen bonds at the interface of myosin and actin. Applying coarse-grained analysis, we identify differences in myosin states dependent on actin-binding, or ATP binding. Our simulations also indicate that the actin propeller twist-angle and nucleotide cleft-angles are influenced by myosin at the actomyosin interface. The torsional rigidity of the myosin-bound filament is also calculated, and is found to be increased compared to previous simulations of the free filament.     

Probing the initial stage of aggregation of the Abeta(10-35)-protein: assessing the propensity for peptide dimerization

Characterization of the early stages of peptide aggregation is of fundamental importance in elucidating the mechanism of the formation of deposits associated with amyloid disease. The initial step in the pathway of aggregation of the Abeta-protein, whose monomeric NMR structure is known, was studied through the simulation of the structure and stability of the peptide dimer in aqueous solution. A protocol based on shape complementarity was used to generate an assortment of possible dimer structures. The structures generated based on shape complementarity were evaluated using rapidly computed estimates of the desolvation and electrostatic interaction energies to identify a putative stable dimer structure. The potential of mean force associated with the dimerization of the peptides in aqueous solution was computed for both the hydrophobic and the electrostatic driven forces using umbrella sampling and classical molecular dynamics simulation at constant temperature and pressure with explicit solvent and periodic boundary conditions. The comparison of the two free energy profiles suggests that the structure of the peptide dimer is determined by the favorable desolvation of the hydrophobic residues at the interface. Molecular dynamics trajectories originating from two putative dimer structures indicate that the peptide dimer is stabilized primarily through hydrophobic interactions, while the conformations of the peptide monomers undergo substantial structural reorganization in the dimerization process. The finding that the phi-dimer may constitute the ensemble of stable Abeta(10-35) dimer has important implications for fibril formation. In particular, the expulsion of water molecules at the interface might be a key event, just as in the oligomerization of Abeta(16-22) fragments. We conjecture that events prior to the nucleation process themselves might involve crossing free energy barriers which depend on the peptide-peptide and peptide-water interactions. Consistent with existing experimental studies, the peptides within the ensemble of aggregated states show no signs of formation of secondary structure.     

Insights into subunit interactions in the heterotetrameric structure of potato ADP-glucose pyrophosphorylase

ADP-glucose pyrophosphorylase, a key allosteric enzyme involved in higher plant starch biosynthesis, is composed of pairs of large (LS) and small subunits (SS). Current evidence indicates that the two subunit types play distinct roles in enzyme function. The LS is involved in mainly allosteric regulation through its interaction with the catalytic SS. Recently the crystal structure of the SS homotetramer has been solved, but no crystal structure of the native heterotetrameric enzyme is currently available. In this study, we first modeled the three-dimensional structure of the LS to construct the heterotetrameric enzyme. Because the enzyme has a 2-fold symmetry, six different dimeric (either up-down or side-by-side) interactions were possible. Molecular dynamics simulations were carried out for each of these possible dimers. Trajectories obtained from molecular dynamics simulations of each dimer were then analyzed by the molecular mechanics/Poisson-Boltzmann surface area method to identify the most favorable dimers, one for up-down and the other for side-by-side. Computational results combined with site directed mutagenesis and yeast two hybrid experiments suggested that the most favorable heterotetramer is formed by LS-SS (side-by-side), and LS-SS (up-down). We further determined the order of assembly during the heterotetrameric structure formation. First, side-by-side LS-SS dimers form followed by the up-down tetramerization based on the relative binding free energies.     

Free energy landscapes for amyloidogenic tetrapeptides dimerization

The oligomerization of four peptide sequences, KFFE, KVVE, KLLE, and KAAE is studied using replica-exchange molecular dynamics simulations with an atomically detailed peptide model. Previous experimental studies reported that of these four peptides, only those containing phenylalanine and valine residues form fibrils. We show that the fibrillogenic propensities of these peptides can be rationalized in terms of the equilibrium thermodynamics of their early oligomers. Thermodynamic stability of dimers, as measured by the temperature of monomer association, is seen to be higher for those peptides that are able to form fibrils. Although the relative high and low stabilities of the KFFE and KAAE dimers arise from their respective high and low interpeptide interaction energies, the higher stability of the KVVE dimer over the KLLE system results from the smaller loss of configurational entropy accompanying the dimerization of KVVE. Free energy landscapes for dimerization are found to be strongly sequence-dependent, with a high free energy barrier separating the monomeric and dimeric states for KVVE, KLLE, and KAAE sequences. In contrast, the most fibrillogenic peptide, KFFE, displayed downhill assembly, indicating enhanced kinetic accessibility of its dimeric states. The dimeric phase for all peptide sequences is found to be heterogeneous, containing both antiparallel beta-sheet structures that can grow into full fibrils as well as disordered dimers acting as on- or off-pathway intermediates for fibrillation.     

Folding, misfolding, and amyloid protofibril formation of WW domain FBP28

We study the folding mechanism of a triple beta-strand WW domain from the Formin binding protein 28 (FBP28) at atomic resolution with explicit water model using replica exchange molecular dynamics computer simulations. Extended sampling over a wide range of temperatures to obtain the free energy, enthalpy, and entropy surfaces as a function of structural coordinates has been performed. Simulations were started from different configurations covering the folded and unfolded states. In the free energy landscape a transition state is identified and its structures and -values are compared with experimental data from a homologous protein, the prolyl-isomerase Pin1 WW domain. A stable intermediate state is found to accumulate during the simulation characterized by the carboxyl-terminal beta-strand 3 having misregistered hydrogen bonds and where the structural heterogeneity is due to nonnative turn II formation. Furthermore, the aggregation behavior of the FBP28 WW domain may be related to one such misfolded structure, which has a much lower free energy of dimer formation than that of the native dimer. Based on the misfolded dimer, aggregation to form protofibril structure is discussed.     

Mechanisms of the self-assembly of EAK16-family peptides into fibrillar and globular structures: molecular dynamics simulations from nano- to micro-seconds

The self-assembly of EAK16-family peptides in a bulk solution was studied using a combination of all-atom and coarse-grained molecular dynamics simulations. In addition, specified concentrations of EAK16 peptides were induced to form fibrillary or globular assemblies in vitro. The results show that the combination of all-atom molecular dynamics simulations on the single- and double-chain levels and coarse-grained simulations on the many-chain level predicts the experimental observations reasonably well. At neutral pH conditions, EAK16-I and EAK16-II assemble into fibrillary structures, whereas EAK16-IV aggregates into globular assemblies. Mechanisms of the formation of fibrillar and globular assemblies are described using the simulation results.     

Targeting the early steps of Abeta16-22 protofibril disassembly by N-methylated inhibitors: a numerical study

Aggregation of the Abeta1-40/Abeta1-42 peptides is a key factor in Alzheimer's disease. Though the inhibitory effect of N-methylated Abeta16-22 (mAbeta16-22) peptides is well characterized in vitro, there is little information on how they disassemble full-length Abeta fibrils or block fibril formation. Here, we report coarse-grained implicit solvent molecular dynamics (MD) and replica exchange molecular dynamics (REMD) simulations on Abeta16-22 and mAbeta16-22 monomers, and then a preformed six-chain Abeta16-22 bilayer with either four copies of Abeta16-22 or four copies of mAbeta16-22. Our simulations show that the effect of N-methylation on mAbeta16-22 monomer is to reduce the density of compact forms. While 100 ns MD trajectories do not reveal any significant differences between the two ten-chain systems, the REMD simulations totaling 1 micros help understand the first steps of Abeta16-22 protofibril disassembly by N-methylated inhibitors. Notably, we find that mAbeta16-22 preferentially interacts with Abeta16-22 by blocking both beta-sheet extension and lateral association of layers, but also by intercalation of the inhibitors allowing sequestration of Abeta16-22 peptides. This third binding mode is particularly appealing for blocking Abeta fibrillogenesis.