Computational backbone mutagenesis of Abeta peptides: probing the role of backbone hydrogen bonds in aggregation

Human Neuronal Calcium Sensor-1 Protein Avoids Histidine Residues To Decrease pH Sensitivity

pH is highly regulated in mammalian central nervous systems. Neuronal calcium sensor-1 (NCS-1) can interact with numerous target proteins. Compared to that in the NCS-1 protein of Caenorhabditis elegans, evolution has avoided the placement of histidine residues at positions 102 and 83 in the NCS-1 protein of humans and Xenopus laevis, possibly to decrease the conformational sensitivity to pH gradients in synaptic processes. We used all-atom molecular dynamics simulations to investigate the effects of amino acid substitutions between species on human NCS-1 by substituting Arg102 and Ser83 for histidine at neutral (R102H and S83H) and acidic pHs (R102Hp and S83Hp). Our cumulative 5 μs simulations revealed that the R102H mutation slightly increases the structural flexibility of loop L2 and the R102Hp mutation decreases protein stability. Community network analysis illustrates that the R102H and S83H mutations weaken the interdomain and strengthen the intradomain communications. Secondary structure contents in the S83H and S83Hp mutants are similar to those in the wild type, whereas the global structural stabilities and salt-bridge probabilities decrease. This study highlights the conformational dynamics effects of the R102H and S83H mutations on the local structural flexibility and global stability of NCS-1, whereas protonated histidine decreases the stability of NCS-1. Thus, histidines at positions 102 and 83 may not be compatible with the function of NCS-1 whether in the neutral or protonated state.

Environmental polarity induces conformational transitions in a helical peptide sequence from bacteriophage T4 lysozyme and its tandem duplicate: a molecular dynamics simulation study

Our recent molecular dynamics (MD) simulation of an insertion/duplication mutant 'L20' of bacteriophage T4 lysozyme demonstrated a solvent induced α→β transition in a loosely held duplicate helical region, while α-helical conformation in the parent region was relatively stabilized by its tertiary interactions with the neighboring residues. The solution NMR of the parent helical sequence, sans its protein context, showed no inherent tendency to adopt a particular secondary structure in pure water but showed α-helical propensity in TFE/water and SDS micelles. In this study we investigate the conformational preference of the 'parent' and 'duplicate' sequences, sans the protein context, in pure water and an apolar TFE/water solution. Apolar TFE/water solution is a model for non-polar protein context. We performed MD simulations of the two peptides, in explicit water and 80% (v/v) TFE/water, using GROMOS 53a6 force field, at 300 K and 1 bar (under NPT conditions). We show that in TFE/water mixture, salt bridges are stabilized by apolar TFE molecules and main chain-main chain hydrogen bonds promote the α-helical conformation, particularly in the duplicate peptide. Solvent exposure, in pure water, resulted in an α→β transition to form a triple stranded β-sheet structure in the 'duplicate' sequence, with a rare psi-loop topology, while a mixture of turn/bend conformations were adopted by the 'parent' sequence. Thus the differences in conformational preference of the parent and duplicate sequence sans protein context, in pure water and TFE/water, implicate the importance of the environment polarity in dictating the peptide conformation. Mechanism of folding of the observed psi-loop in the duplicate sequence gives insights into folding of this rare β-sheet topology.

Molecular dynamics study of an insertion/duplication mutant of bacteriophage T4 lysozyme reveals the nature of α→β transition in full protein context

An α→β transition underlies the first step of disease causing amyloidogenesis in many proteins. In view of this, many studies have been carried out using peptide models to characterize these secondary structural transitions. In this paper we show that an insertion/duplication mutant 'L20' of bacteriophage T4 lysozyme (M. Sagermann, W. A. Baase and B. W. Matthews, Proc. Natl. Acad. Sci. U.S.A., 1999, 96, 6078) displays an α→β transition. We performed molecular dynamics (MD) simulation of L20, using the GROMACS package of programs and united atom GROMOS 53a6 force field for a time period of 600 ns at 300 K, in explicit water. Our MD simulation demonstrated that the transition occurs in a duplicated α-helical region inserted tandemly at the N-terminus of the 'parent' helix. We show that a C-terminal β-sheet anchors the parent helix while the loosely held N-terminal loop in the duplicate region is vulnerable to solvent attack and thus undergoes an α→β transition. Main chain-solvent interactions were seen to stabilize the observed β-structure. Thus L20 serves as a good protein model for characterization of α→β transition in a full length protein.

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.

Biochemistry of amyloid β-protein and amyloid deposits in Alzheimer disease

Progressive cerebral deposition of the amyloid β-protein (Aβ) in brain regions serving memory and cognition is an invariant and defining feature of Alzheimer disease. A highly similar but less robust process accompanies brain aging in many nondemented humans, lower primates, and some other mammals. The discovery of Aβ as the subunit of the amyloid fibrils in meningocerebral blood vessels and parenchymal plaques has led to innumerable studies of its biochemistry and potential cytotoxic properties. Here we will review the discovery of Aβ, numerous aspects of its complex biochemistry, and current attempts to understand how a range of Aβ assemblies, including soluble oligomers and insoluble fibrils, may precipitate and promote neuronal and glial alterations that underlie the development of dementia. Although the role of Aβ as a key molecular factor in the etiology of Alzheimer disease remains controversial, clinical trials of amyloid-lowering agents, reviewed elsewhere in this book, are poised to resolve the question of its pathogenic primacy.

Molecular dynamics simulations of low-ordered alzheimer β-amyloid oligomers from dimer to hexamer on self-assembled monolayers

Accumulation of small soluble oligomers of amyloid-β (Aβ) in the human brain is thought to play an important pathological role in Alzheimer's disease. The interaction of these Aβ oligomers with cell membrane and other artificial surfaces is important for the understanding of Aβ aggregation and toxicity mechanisms. Here, we present a series of exploratory molecular dynamics (MD) simulations to study the early adsorption and conformational change of Aβ oligomers from dimer to hexamer on three different self-assembled monolayers (SAMs) terminated with CH(3), OH, and COOH groups. Within the time scale of MD simulations, the conformation, orientation, and adsorption of Aβ oligomers on the SAMs is determined by complex interplay among the size of Aβ oligomers, the surface chemistry of the SAMs, and the structure and dynamics of interfacial waters. Energetic analysis of Aβ adsorption on the SAMs reveals that Aβ adsorption on the SAMs is a net outcome of different competitions between dominant hydrophobic Aβ-CH(3)-SAM interactions and weak CH(3)-SAM-water interactions, between dominant electrostatic Aβ-COOH-SAM interactions and strong COOH-SAM-water interactions, and between comparable hydrophobic and electrostatic Aβ-OH-SAM interactions and strong OH-SAM-water interactions. Atomic force microscopy images also confirm that all of three SAMs can induce the adsorption and polymerization of Aβ oligomers. Structural analysis of Aβ oligomers on the SAMs shows a dramatic increase in structural stability and β-sheet content from dimer to trimer, suggesting that Aβ trimer could act as seeds for Aβ polymerization on the SAMs. This work provides atomic-level understanding of Aβ peptides at interface.

Binding of inositol stereoisomers to model amyloidogenic peptides

The self-aggregation of proteins into amyloid fibrils is a pathological hallmark of numerous incurable diseases such as Alzheimer's disease. scyllo-Inositol is a stereochemistry-dependent in vitro inhibitor of amyloid formation. As the first step to elucidate its mechanism of action, we present molecular dynamics simulations of scyllo-inositol and its inactive stereoisomer, chiro-inositol, with simple peptide models, alanine dipeptide (ADP) and (Gly-Ala)(4). We characterize molecular interactions and compute equilibrium binding constants between inositol and ADP as well as, successively, monomers, amorphous aggregates, and fibril-like β-sheet aggregates of (Gly-Ala)(4). Inositol interacts weakly with all peptide systems considered, with millimolar to molar affinities, and displaces the conformational equilibria of ADP but not of the (Gly-Ala)(4) systems. However, scyllo- and chiro-inositol adopt different binding modes on the surface of β-sheet aggregates. These results suggest that inositol does not inhibit amyloid formation by breaking up preformed aggregates but rather by binding to the surface of prefibrillar aggregates.

Mapping conformational ensembles of aβ oligomers in molecular dynamics simulations

Although the oligomers formed by Aβ peptides appear to be the primary cytotoxic species in Alzheimer's disease, detailed information about their structures appears to be lacking. In this article, we use exhaustive replica exchange molecular dynamics and an implicit solvent united-atom model to study the structural properties of Aβ monomers, dimers, and tetramers. Our analysis suggests that the conformational ensembles of Aβ dimers and tetramers are very similar, but sharply distinct from those sampled by the monomers. The key conformational difference between monomers and oligomers is the formation of β-structure in the oligomers occurring together with the loss of intrapeptide interactions and helix structure. Our simulations indicate that, independent of oligomer order, the Aβ aggregation interface is largely confined to the sequence region 10-23, which forms the bulk of interpeptide interactions. We show that the fractions of β structure computed in our simulations and measured experimentally are in good agreement.