Search

The structure of the 21-30 fragment of the amyloid beta-protein (Abeta) was investigated by ion mobility mass spectrometry and replica exchange dynamics simulations. Mutations associated with familial Alzheimer's disease (E22G, E22Q, E22K, and D23N) of Abeta(21-30) were also studied, in order to understand any structural changes that might occur with these substitutions. The structure of the WT peptide shows a bend and a perpendicular turn in the backbone which is maintained by a network of D23 hydrogen bonding. Results for the mutants show that substitutions at E22 do little to alter the overall structure of the fragment. A substitution at D23 resulted in a change of structure for Abeta(21-30). A comparison of these gas-phase studies to previous solution-phase studies reveals that the peptide can fold in the absence of solvent to a structure also seen in solution, highlighting the important role of the D23 hydrogen bonding network in stabilizing the fragment's folded structure.

The structure of the 21-30 fragment of the amyloid beta-protein (Abeta) was investigated by ion mobility mass spectrometry and replica exchange dynamics simulations. Mutations associated with familial Alzheimer's disease (E22G, E22Q, E22K, and D23N) of Abeta(21-30) were also studied, in order to understand any structural changes that might occur with these substitutions. The structure of the WT peptide shows a bend and a perpendicular turn in the backbone which is maintained by a network of D23 hydrogen bonding. Results for the mutants show that substitutions at E22 do little to alter the overall structure of the fragment. A substitution at D23 resulted in a change of structure for Abeta(21-30). A comparison of these gas-phase studies to previous solution-phase studies reveals that the peptide can fold in the absence of solvent to a structure also seen in solution, highlighting the important role of the D23 hydrogen bonding network in stabilizing the fragment's folded structure.

The effect of single amino acid substitutions associated with the Italian (E22K), Arctic (E22G), Dutch (E22Q) and Iowa (D23N) familial forms of Alzheimer's disease and cerebral amyloid angiopathy on the structure of the 21-30 fragment of the Alzheimer amyloid beta-protein (Abeta) is investigated by replica-exchange molecular dynamics simulations. The 21-30 segment has been shown in our earlier work to adopt a bend structure in solution that may serve as the folding nucleation site for Abeta. Our simulations reveal that the 24-28 bend motif is retained in all E22 mutants, suggesting that mutations involving residue E22 may not affect the structure of the folding nucleation site of Abeta. Enhanced aggregation in Abeta with familial Alzheimer's disease substitutions may result from the depletion of the E22-K28 salt bridge, which destabilizes the bend structure. Alternately, the E22 mutations may affect longer-range interactions outside the 21-30 segment that can impact the aggregation of Abeta. Substituting at residue D23, on the other hand, leads to the formation of a turn rather than a bend motif, implying that in contrast to E22 mutants, the D23N mutant may affect monomer Abeta folding and subsequent aggregation. Our simulations suggest that the mechanisms by which E22 and D23 mutations affect the folding and aggregation of Abeta are fundamentally different.

The effect of single amino acid substitutions associated with the Italian (E22K), Arctic (E22G), Dutch (E22Q) and Iowa (D23N) familial forms of Alzheimer's disease and cerebral amyloid angiopathy on the structure of the 21-30 fragment of the Alzheimer amyloid beta-protein (Abeta) is investigated by replica-exchange molecular dynamics simulations. The 21-30 segment has been shown in our earlier work to adopt a bend structure in solution that may serve as the folding nucleation site for Abeta. Our simulations reveal that the 24-28 bend motif is retained in all E22 mutants, suggesting that mutations involving residue E22 may not affect the structure of the folding nucleation site of Abeta. Enhanced aggregation in Abeta with familial Alzheimer's disease substitutions may result from the depletion of the E22-K28 salt bridge, which destabilizes the bend structure. Alternately, the E22 mutations may affect longer-range interactions outside the 21-30 segment that can impact the aggregation of Abeta. Substituting at residue D23, on the other hand, leads to the formation of a turn rather than a bend motif, implying that in contrast to E22 mutants, the D23N mutant may affect monomer Abeta folding and subsequent aggregation. Our simulations suggest that the mechanisms by which E22 and D23 mutations affect the folding and aggregation of Abeta are fundamentally different.

Oligomeric, neurotoxic amyloid protein assemblies are thought to be causative agents in Alzheimer's and other neurodegenerative diseases. Development of oligomer-specific therapeutic agents requires a mechanistic understanding of the oligomerization process. This is a daunting task because amyloidogenic protein oligomers often are metastable and comprise structurally heterogeneous populations in equilibrium with monomers and fibrils. A single methodological approach cannot elucidate the entire protein assembly process. An integrated multidisciplinary program is required. We discuss here the synergistic application of in hydro, in vacuo, and in silico methods to the study of the amyloid beta-protein, the key pathogenetic agent in Alzheimer's disease.

Oligomeric, neurotoxic amyloid protein assemblies are thought to be causative agents in Alzheimer's and other neurodegenerative diseases. Development of oligomer-specific therapeutic agents requires a mechanistic understanding of the oligomerization process. This is a daunting task because amyloidogenic protein oligomers often are metastable and comprise structurally heterogeneous populations in equilibrium with monomers and fibrils. A single methodological approach cannot elucidate the entire protein assembly process. An integrated multidisciplinary program is required. We discuss here the synergistic application of in hydro, in vacuo, and in silico methods to the study of the amyloid beta-protein, the key pathogenetic agent in Alzheimer's disease.

Folding and self-assembly of the 42-residue amyloid beta-protein (Abeta) are linked to Alzheimer's disease (AD). The 21-30 region of Abeta, Abeta(21-30), is resistant to proteolysis and is believed to nucleate the folding of full-length Abeta. The conformational space accessible to the Abeta(21-30) peptide is investigated by using replica exchange molecular dynamics simulations in explicit solvent. Conformations belonging to the global free energy minimum (the "native" state) from simulation are in good agreement with reported NMR structures. These conformations possess a bend motif spanning the central residues V24-K28. This bend is stabilized by a network of hydrogen bonds involving the side chain of residue D23 and the amide hydrogens of adjacent residues G25, S26, N27, and K28, as well as by a salt bridge formed between side chains of K28 and E22. The non-native states of this peptide are compact and retain a native-like bend topology. The persistence of structure in the denatured state may account for the resistance of this peptide to protease degradation and aggregation, even at elevated temperatures.

Folding and self-assembly of the 42-residue amyloid beta-protein (Abeta) are linked to Alzheimer's disease (AD). The 21-30 region of Abeta, Abeta(21-30), is resistant to proteolysis and is believed to nucleate the folding of full-length Abeta. The conformational space accessible to the Abeta(21-30) peptide is investigated by using replica exchange molecular dynamics simulations in explicit solvent. Conformations belonging to the global free energy minimum (the "native" state) from simulation are in good agreement with reported NMR structures. These conformations possess a bend motif spanning the central residues V24-K28. This bend is stabilized by a network of hydrogen bonds involving the side chain of residue D23 and the amide hydrogens of adjacent residues G25, S26, N27, and K28, as well as by a salt bridge formed between side chains of K28 and E22. The non-native states of this peptide are compact and retain a native-like bend topology. The persistence of structure in the denatured state may account for the resistance of this peptide to protease degradation and aggregation, even at elevated temperatures.

The structural properties of the Abeta42 peptide, a main constituent of the amyloid plaques formed in Alzheimer's disease, were investigated through a combination of ion-mobility mass spectrometry and theoretical modeling. Replica exchange molecular dynamics simulations using a fully atomic description of the peptide and implicit water solvent were performed on the -3 charge state of the peptide, its preferred state under experimental conditions. Equilibrated structures at 300 K were clustered into three distinct families with similar structural features within a family and with significant root mean square deviations between families. An analysis of secondary structure indicates the Abeta42 peptide conformations are dominated by loops and turns but show some helical structure in the C-terminal hydrophobic tail. A second calculation on Abeta42 in a solvent-free environment yields compact structures turned "inside out" from the solution structures (hydrophobic parts on the outside, polar parts on the inside). Ion mobility experiments on the Abeta42 -3 charge state electrosprayed from solution yield a bimodal arrival time distribution. This distribution can be quantitatively fit using cross-sections from dehydrated forms of the three families of calculated solution structures and the calculated solvent-free family of structures. Implications of the calculations on the early stages of aggregation of Abeta42 are discussed.

The structural properties of the Abeta42 peptide, a main constituent of the amyloid plaques formed in Alzheimer's disease, were investigated through a combination of ion-mobility mass spectrometry and theoretical modeling. Replica exchange molecular dynamics simulations using a fully atomic description of the peptide and implicit water solvent were performed on the -3 charge state of the peptide, its preferred state under experimental conditions. Equilibrated structures at 300 K were clustered into three distinct families with similar structural features within a family and with significant root mean square deviations between families. An analysis of secondary structure indicates the Abeta42 peptide conformations are dominated by loops and turns but show some helical structure in the C-terminal hydrophobic tail. A second calculation on Abeta42 in a solvent-free environment yields compact structures turned "inside out" from the solution structures (hydrophobic parts on the outside, polar parts on the inside). Ion mobility experiments on the Abeta42 -3 charge state electrosprayed from solution yield a bimodal arrival time distribution. This distribution can be quantitatively fit using cross-sections from dehydrated forms of the three families of calculated solution structures and the calculated solvent-free family of structures. Implications of the calculations on the early stages of aggregation of Abeta42 are discussed.

The amyloid beta-protein (Abeta) is a seminal neuropathic agent in Alzheimer's disease (AD). Recent evidence points to soluble Abeta oligomers as the probable neurotoxic species. Among the naturally occurring Abeta peptides, the 42-residue form Abeta42 is linked particularly strongly with AD, even though it is produced at approximately 10% of the levels of the more abundant 40-residue form Abeta40. Here, we apply mass spectrometry and ion mobility to the study of Abeta42 and its Pro19 alloform. The Phe19 --> Pro19 substitution blocks fibril formation by [Pro19]Abeta42. Evidence indicates that solution-like structures of Abeta monomers are electrosprayed and characterized. Unfiltered solutions of Abeta42 produce only monomers and large oligomers, whereas [Pro19]Abeta42 solutions produce abundant monomers, dimers, trimers, and tetramers but no large oligomers. When passed through a 10,000 amu filter and immediately sampled, Abeta42 solutions produce monomers, dimers, tetramers, hexamers, and an aggregate of two hexamers that may be the first step in protofibril formation. These results are consistent with recently published photochemical cross-linking data and lend support to recent aggregation mechanisms proposed by Bitan, Teplow, and co-workers [J. Biol. Chem. 2003, 278, 34882-34889].

The amyloid beta-protein (Abeta) is a seminal neuropathic agent in Alzheimer's disease (AD). Recent evidence points to soluble Abeta oligomers as the probable neurotoxic species. Among the naturally occurring Abeta peptides, the 42-residue form Abeta42 is linked particularly strongly with AD, even though it is produced at approximately 10% of the levels of the more abundant 40-residue form Abeta40. Here, we apply mass spectrometry and ion mobility to the study of Abeta42 and its Pro19 alloform. The Phe19 --> Pro19 substitution blocks fibril formation by [Pro19]Abeta42. Evidence indicates that solution-like structures of Abeta monomers are electrosprayed and characterized. Unfiltered solutions of Abeta42 produce only monomers and large oligomers, whereas [Pro19]Abeta42 solutions produce abundant monomers, dimers, trimers, and tetramers but no large oligomers. When passed through a 10,000 amu filter and immediately sampled, Abeta42 solutions produce monomers, dimers, tetramers, hexamers, and an aggregate of two hexamers that may be the first step in protofibril formation. These results are consistent with recently published photochemical cross-linking data and lend support to recent aggregation mechanisms proposed by Bitan, Teplow, and co-workers [J. Biol. Chem. 2003, 278, 34882-34889].

Amyloid fibrils, large ordered aggregates of amyloid [beta] proteins (A[beta]), are clinical hallmarks of Alzheimer's disease (AD). The aggregation properties of amyloid [beta] proteins can be strongly affected by single-point mutations at positions 22 and 23. The Dutch mutation involves a substitution at position 22 (E22Q) and leads to increased deposition rates of the A[beta]E22Q peptide onto preseeded fibrils. We investigate the effect of the E22Q mutation on two key regions involved in the folding and aggregation of the A[beta] peptide through replica exchange molecular dynamics simulations of the 15-28 fragment of the A[beta] peptide. The A[beta]15-28 peptide encompasses the 22-28 region that constitutes the most structured part of the A[beta] peptide (the E22-K28 bend), as well as the central hydrophobic cluster (CHC) (segment 17-21), the primary docking site for A[beta] monomers depositing onto fibrils. Our simulations show that the 22-28 bend is preserved in the A[beta](15-28) peptide and that the CHC, which is mostly unstructured, interacts with this bend region. The E22Q mutation does not affect the structure of the bend but weakens the interactions between the CHC and the bend. This leads to an increased population of [beta]-structure in the CHC. Our analysis of the fibril elongation reaction reveals that the CHC adopts a [beta]-strand conformation in the transition state ensemble, and that the E22Q mutation increases aggregation rates by lowering the barrier for A[beta] monomer deposition onto a fibril. Thermodynamic signatures of this enhanced fibrillization process from our simulations are in good agreement with experimental observations. Alzheimer's disease | explicit water | familial type | molecular dynamics simulations | replica exchange

Increasingly, theoretical studies of proteins focus on large systems. This trend demands the development of computational models that are fast, to overcome the growing complexity, and accurate, to capture the physically relevant features. To address this demand, we introduce a protein model that uses all-atom architecture to ensure the highest level of chemical detail while employing effective pair potentials to represent the effect of solvent to achieve the maximum speed. The effective potentials are derived for amino acid residues based on the condition that the solvent-free model matches the relevant pair-distribution functions observed in explicit solvent simulations. As a test, the model is applied to alanine polypeptides. For the chain with 10 amino acid residues, the model is found to reproduce properly the native state and its population. Small discrepancies are observed for other folding properties and can be attributed to the approximations inherent in the model. The transferability of the generated effective potentials is investigated in simulations of a longer peptide with 25 residues. A minimal set of potentials is identified that leads to qualitatively correct results in comparison with the explicit solvent simulations. Further tests, conducted for multiple peptide chains, show that the transferable model correctly reproduces the experimentally observed tendency of polyalanines to aggregate into β-sheets more strongly with the growing length of the peptide chain. Taken together, the reported results suggest that the proposed model could be used to succesfully simulate folding and aggregation of small peptides in atomic detail. Further tests are needed to assess the strengths and limitations of the model more thoroughly. [ABSTRACT FROM AUTHOR]

In a preceding paper [J. Chem. Phys. 131, 154103 (2009)], we introduced a new, hybrid explicit/implicit method to treat electrostatic interactions in computer simulations, and tested its performance for liquid water. In this paper, we report further tests of this method, termed the image-charge solvation model (ICSM), in simulations of ions solvated in water. We find that our model can faithfully reproduce known solvation properties of sodium and chloride ions. The charging free energy of a single sodium ion is in excellent agreement with the estimates by other electrostatics methods, while offering much lower finite-size errors. Similarly, the potentials of mean force computed for Na-Cl, Na-Na, and Cl-Cl pairs closely reproduce those reported previously. Collectively, our results demonstrate the superior accuracy of the proposed ICSM method for simulations of mixed media. [ABSTRACT FROM AUTHOR]

In this paper, a new solvation model is proposed for simulations of biomolecules in aqueous solutions that combines the strengths of explicit and implicit solvent representations. Solute molecules are placed in a spherical cavity filled with explicit water, thus providing microscopic detail where it is most needed. Solvent outside of the cavity is modeled as a dielectric continuum whose effect on the solute is treated through the reaction field corrections. With this explicit/implicit model, the electrostatic potential represents a solute molecule in an infinite bath of solvent, thus avoiding unphysical interactions between periodic images of the solute commonly used in the lattice-sum explicit solvent simulations. For improved computational efficiency, our model employs an accurate and efficient multiple-image charge method to compute reaction fields together with the fast multipole method for the direct Coulomb interactions. To minimize the surface effects, periodic boundary conditions are employed for nonelectrostatic interactions. The proposed model is applied to study liquid water. The effect of model parameters, which include the size of the cavity, the number of image charges used to compute reaction field, and the thickness of the buffer layer, is investigated in comparison with the particle-mesh Ewald simulations as a reference. An optimal set of parameters is obtained that allows for a faithful representation of many structural, dielectric, and dynamic properties of the simulated water, while maintaining manageable computational cost. With controlled and adjustable accuracy of the multiple-image charge representation of the reaction field, it is concluded that the employed model achieves convergence with only one image charge in the case of pure water. Future applications to pKa calculations, conformational sampling of solvated biomolecules and electrolyte solutions are briefly discussed. [ABSTRACT FROM AUTHOR]

The performance of reaction-field methods to treat electrostatic interactions is tested in simulations of ions solvated in water. The potential of mean force between sodium chloride pair of ions and between side chains of lysine and aspartate are computed using umbrella sampling and molecular dynamics simulations. It is found that in comparison with lattice sum calculations, the charge-group-based approaches to reaction-field treatments produce a large error in the association energy of the ions that exhibits strong systematic dependence on the size of the simulation box. The atom-based implementation of the reaction field is seen to (i) improve the overall quality of the potential of mean force and (ii) remove the dependence on the size of the simulation box. It is suggested that the atom-based truncation be used in reaction-field simulations of mixed media. [ABSTRACT FROM AUTHOR]

Surface-tethered proteins are increasingly being used in a variety of experimental situations, and they are the basis for many new technologies. Nevertheless, a thorough understanding of how a surface can impact the native state stability of an attached protein is lacking. In this work, the authors use molecular dynamics simulations of a model β-barrel protein to investigate how surface tethering influences native state stability. They find that stability, as measured by the folding temperature Tf, can be either increased, decreased, or remain unchanged as a result of tethering. Observed shifts are highly dependent on the location of residue used as the tether point, and stability is influenced by a number of factors, both energetic and entropic. These factors include native state vibrations, loss of bulk unfolded conformations, changes to the unfolded state ensemble, and the emergence of an entropic term not present for the bulk protein. They discuss each of these contributions in detail and comment on their relative importance and connection to experiment. [ABSTRACT FROM AUTHOR]

The dimerization of polyalanine peptides in a hydrophobic environment was explored using replica exchange molecular dynamics simulations. A nonpolar solvent (cyclohexane) was used to mimic, among other hydrophobic environments, the hydrophobic interior of a membrane in which the peptides are fully embedded. Our simulations reveal that while the polyalanine monomer preferentially adopts a β-hairpin conformation, dimeric phases exist in an equilibrium between random coil, α-helical, β-sheet, and β-hairpin states. A thermodynamic characterization of the dimeric phases reveals that electric dipole-dipole interactions and optimal side-chain packing stabilize α-helical conformations, while hydrogen bond interactions favor β-sheet conformations. Possible pathways leading to the formation of α-helical and β-sheet dimers are discussed. [ABSTRACT FROM AUTHOR]

Replica exchange molecular dynamics simulations were performed to investigate the effects of different electrostatic treatments on the structure and thermodynamics of a small beta-hairpin forming peptide. Simulations on the beta-hairpin forming peptide and on additional control peptides indicate that generalized reaction field treatment of electrostatics offers a satisfactory compromise between accuracy and computational speed for the identification of low-energy conformations.

Mass spectrometry and ion mobility is applied to the study of A(beta)42 and its Pro(super 19) alloform. Unfiltered solutions of A(beta)42 produce only monomers and large oligomers, whereas [Pro(super 19)]A(beta)42 solutions produce abundant monomers, dimers, trimers, and tetramers but no large oligomers.

Foldingand aggregation lie on competing reaction pathways in proteins.Altering the occupancy of one pathway is automatically relayed tothe other pathway, leading to a shift in the balance between the twoprocesses. In particular, it is known that the stabilization of thenative state through mutations or solvent alterations is able to haltaggregation. In this work, we explore the feasibility of using externalelectric field as an agent preventing aggregation through the promotionof folding. We use an atomically accurate protein model and computersimulations to investigate folding and aggregation of alanine polypeptidesin electric field of varying strength. The studied peptides are mostlyunstructured in the absence of the field but experience a transitioninto α-helical states when the field is applied. The transitionis accompanied by the disassembly of preseeded stacked β-sheets,which are used as a model of amyloid fibrils, suggesting that electricfield can be employed to control aggregation propensity of intrinsicallydisordered peptides. According to our calculations, the strength ofthe field required for the disaggregation could be suitable for bothcontrolled in vitroexperiments as well as for experimentson live cells. Additionally, our estimates suggest that endogenouselectric fields may have a significant effect on in vivoamyloid formation. [ABSTRACT FROM AUTHOR]