α-helix to ß-hairpin transition of human amylin monomer.

The human islet amylin polypeptide is produced along with insulin by pancreatic islets. Under some circumstances, amylin can aggregate to form amyloid fibrils, whose presence in pancreatic cells is a common pathological feature of Type II diabetes. A growing body of evidence indicates that small, early stage aggregates of amylin are cytotoxic. A better understanding of the early stages of the amylin aggregation process and, in particular, of the nucleation events leading to fibril growth could help identify therapeutic strategies. Recent studies have shown that, in dilute solution, human amylin can adopt an α-helical conformation, a ß-hairpin conformation, or an unstructured coil conformation. While such states have comparable free energies, the ß-hairpin state exhibits a large propensity towards aggregation. In this work, we present a detailed computational analysis of the folding pathways that arise between the various conformational states of human amylin in water. A free energy surface for amylin in explicit water is first constructed by resorting to advanced sampling techniques. Extensive transition path sampling simulations are then employed to identify the preferred folding mechanisms between distinct minima on that surface. Our results reveal that the α-helical conformer of amylin undergoes a transformation into the ß-hairpin monomer through one of two mechanisms. In the first, misfolding begins through formation of specific contacts near the turn region, and proceeds via a zipping mechanism. In the second, misfolding occurs through an unstructured coil intermediate. The transition states for these processes are identified. Taken together, the findings presented in this work suggest that the inter-conversion of amylin between an α-helix and a ß-hairpin is an activated process and could constitute the nucleation event for fibril growth.

2DIR spectroscopy of human amylin fibrils reflects stable ß-sheet structure.

The aggregation of human amylin to form amyloid contributes to islet ß-cell dysfunction in type 2 diabetes. Studies of amyloid formation have been hindered by the low structural resolution or relatively modest time resolution of standard methods. Two-dimensional infrared (2DIR) spectroscopy, with its sensitivity to protein secondary structures and its intrinsic fast time resolution, is capable of capturing structural changes during the aggregation process. Moreover, isotope labeling enables the measurement of residue-specific information. The diagonal line widths of 2DIR spectra contain information about dynamics and structural heterogeneity of the system. We illustrate the power of a combined atomistic molecular dynamics simulation and theoretical and experimental 2DIR approach by analyzing the variation in diagonal line widths of individual amide I modes in a series of labeled samples of amylin amyloid fibrils. The theoretical and experimental 2DIR line widths suggest a "W" pattern, as a function of residue number. We show that large line widths result from substantial structural disorder and that this pattern is indicative of the stable secondary structure of the two ß-sheet regions. This work provides a protocol for bridging MD simulation and 2DIR experiments for future aggregation studies.

Structural effects and translocation of doxorubicin in a DPPC/Chol bilayer: the role of cholesterol.

We use molecular dynamics simulations to characterize the influence of cholesterol (Chol) on the interaction between the anticancer drug doxorubicin (DOX) and a dipalmitoyl phosphatidylcholine/Chol lipid bilayer. We calculate the potential of mean force, which gives us an estimate of the free energy barrier for DOX translocation across the membrane. We find free energy barriers of 23.1 ± 3.1 k(B)T, 36.8 ± 5.1 k(B)T, and 54.5 ± 4.7 k(B)T for systems composed of 0%, 15%, and 30% Chol, respectively. Our predictions agree with Arrhenius activation energies from experiments using phospholipid membranes, including 20 k(B)T for 0% Chol and 37.2 k(B)T for 20% Chol. The location of the free energy barrier for translocation across the bilayer is dependent on composition. As Chol concentration increases, this barrier changes from the release of DOX into the water to flip-flop over the membrane center. The drug greatly affects local membrane structure by attracting dipalmitoyl phosphatidylcholine headgroups, curving the membrane, and allowing water penetration. Despite its hydrophobicity, DOX facilitates water transport via its polar groups.

Effect of trehalose on the interaction of Alzheimer's Aß-peptide and anionic lipid monolayers.

The interaction of amyloid ß-peptide (Aß) with cell membranes is believed to play a central role in the pathogenesis of Alzheimer's disease. In particular, recent experimental evidence indicates that bilayer and monolayer membranes accelerate the aggregation and amyloid fibril formation rate of Aß. Understanding that interaction could help develop therapeutic strategies for treatment of the disease. Trehalose, a disaccharide of glucose, has been shown to be effective in preventing the aggregation of numerous proteins. It has also been shown to delay the onset of certain amyloid-related diseases in a mouse model. Using Langmuir monolayers and molecular simulations of the corresponding system, we study several thermodynamic and kinetic aspects of the insertion of Aß peptide into DPPG monolayers in water and trehalose subphases. In the water subphase, the insertion of the Aß peptide into the monolayer exhibits a lag time which decreases with increasing temperature of the subphase. In the presence of trehalose, the lag time is completely eliminated and peptide insertion is completed within a shorter time period compared to that observed in pure water. Molecular simulations show that more peptide is inserted into the monolayer in the water subphase, and that such insertion is deeper. The peptide at the monolayer interface orients itself parallel to the monolayer, while it inserts with an angle of 50° in the trehalose subphase. Simulations also show that trehalose reduces the conformational change that the peptide undergoes when it inserts into the monolayer. This observation helps explain the experimentally observed elimination of the lag time by trehalose and the temperature dependence of the lag time in the water subphase.

Stable and metastable states of human amylin in solution.

Patients with type II diabetes exhibit fibrillar deposits of human amylin protein in the pancreas. It has been proposed that amylin oligomers arising along the aggregation or fibril-formation pathways are important in the genesis of the disease. In a step toward understanding these aggregation pathways, in this work we report the conformational preferences of human amylin monomer in solution using molecular simulations and infrared experiments. In particular, we identify a stable conformer that could play a key role in aggregation. We find that amylin adopts three stable conformations: one with an α-helical segment comprising residues 9-17 and a short antiparallel ß-sheet comprising residues 24-28 and 31-35; one with an extended antiparallel ß-hairpin with the turn region comprising residues 20-23; and one with no particular structure. Using detailed calculations, we determine the relative stability of these various conformations, finding that the ß-hairpin conformation is the most stable, followed by the α-helical conformation, and then the unstructured coil. To test our predicted structure, we calculate its infrared spectrum in the amide I stretch regime, which is sensitive to secondary structure through vibrational couplings and linewidths, and compare it to experiment. We find that theoretically predicted spectra are in good agreement with the experimental line shapes presented herein. The implications of the monomer secondary structures on its aggregation pathway and on its interaction with cell membranes are discussed.

GNNQQNY--investigation of early steps during amyloid formation.

Protein aggregation has been implicated in the pathology of several neurodegenerative diseases, and a better understanding of how it proceeds is essential for the development of therapeutic strategies. Recently, the amyloidogenic heptapeptide GNNQQNY has emerged as a molecule of choice for fundamental studies of protein aggregation. A number of experimental and computational studies have examined the structure of the GNNQQNY aggregate. Less work, however, has been aimed at understanding its aggregation pathway. In this study, we present a detailed computational analysis of such a pathway. To that end, transition path sampling Monte Carlo simulations are used to examine the dimerization process. A statistical analysis of the reaction pathways shows that the dimerization reaction proceeds via a zipping mechanism, initiated with the formation of distinct contacts at the third residue (N). Asparagine residues are found to play a key role in the early stages of aggregation. And, contrary to previous belief, it is also shown that the tyrosine terminal group is not required to stabilize the dimer. In fact, an asparagine residue leads to faster aggregation of the peptide.

2D IR line shapes probe ovispirin peptide conformation and depth in lipid bilayers.

We report a structural study on the membrane binding of ovispirin using 2D IR line shape analysis, isotope labeling, and molecular dynamics simulations. Ovispirin is an antibiotic polypeptide that binds to the surfaces of membranes as an alpha-helix. By resolving individual backbone vibrational modes (amide I) using 1-(13)C=(18)O labeling, we measured the 2D IR line shapes for 15 of the 18 residues in this peptide. A comparison of the line shapes reveals an oscillation in the inhomogeneous line width that has a period equal to that of an alpha-helix (3.6 amino acids). The periodic trend is caused by the asymmetric environment of the membrane bilayer that exposes one face of the alpha-helix to much stronger environmental electrostatic forces than the other. We compare our experimental results to 2D IR line shapes calculated using the lowest free energy structure identified from molecular dynamics simulations. These simulations predict a periodic trend similar to the experiment and lead us to conclude that ovispirin lies in the membrane just below the headgroups, is tilted, and may be kinked. Besides providing insight into the antibiotic mechanism of ovispirin, our procedure provides an infrared method for studying peptide and protein structures that relies on the natural vibrational modes of the backbone. It is a complementary method to other techniques that utilize line shapes, such as fluorescence, NMR, and ESR spectroscopies, because it does not require mutations, the spectra can be quantitatively simulated using molecular dynamics, and the technique can be applied to difficult-to-study systems like ion channels, aggregated proteins, and kinetically evolving systems.

Solution structures of rat amylin peptide: simulation, theory, and experiment.

Amyloid deposits of amylin in the pancreas are an important characteristic feature found in patients with Type-2 diabetes. The aggregate has been considered important in the disease pathology and has been studied extensively. However, the secondary structures of the individual peptide have not been clearly identified. In this work, we present detailed solution structures of rat amylin using a combination of Monte Carlo and molecular dynamics simulations. A new Monte Carlo method is presented to determine the free energy of distinct biomolecular conformations. Both folded and random-coil conformations of rat amylin are observed in water and their relative stability is examined in detail. The former contains an alpha-helical segment comprised of residues 7-17. We find that at room temperature the folded structure is more stable, whereas at higher temperatures the random-coil structure predominates. From the configurations and weights we calculate the alpha-carbon NMR chemical shifts, with results that are in reasonable agreement with experiments of others. We also calculate the infrared spectrum in the amide I stretch regime, and the results are in fair agreement with the experimental line shape presented herein.

Effect of trehalose on amyloid beta (29-40)-membrane interaction.

A growing body of experimental evidence indicates that the interaction between amyloid beta peptide and lipid bilayer membranes plays an important role in the development of Alzheimer disease. Recent experimental evidence also suggests that trehalose, a simple disaccharide, reduces the toxicity of amyloid beta peptide. Molecular simulations are used to examine the effect of trehalose on the conformational stability of amyloid beta peptide in aqueous solution and its effect on the interaction between amyloid beta peptide and a model phospholipid bilayer membrane. It is found that, in aqueous solution, the peptide exhibits a random coil conformation but, in the presence of trehalose, it adopts an alpha helical conformation. It is then shown that the insertion of amyloid beta peptide into a membrane is more favorable when the peptide is folded into an alpha-helix than in a random coil conformation, thereby suggesting that trehalose promotes the insertion of alpha-helical amyloid beta into biological membranes.

Folding of polyglutamine chains.

Long polyglutamine chains have been associated with a number of neurodegenerative diseases. These include Huntington's disease, where expanded polyglutamine (PolyQ) sequences longer than 36 residues are correlated with the onset of symptoms. In this paper we study the folding pathway of a 54-residue PolyQ chain into a beta-helical structure. Transition path sampling Monte Carlo simulations are used to generate unbiased reactive pathways between unfolded configurations and the folded beta-helical structure of the polyglutamine chain. The folding process is examined in both explicit water and an implicit solvent. Both models reveal that the formation of a few critical contacts is necessary and sufficient for the molecule to fold. Once the primary contacts are formed, the fate of the protein is sealed and it is largely committed to fold. We find that, consistent with emerging hypotheses about PolyQ aggregation, a stable beta-helical structure could serve as the nucleus for subsequent polymerization of amyloid fibrils. Our results indicate that PolyQ sequences shorter than 36 residues cannot form that nucleus, and it is also shown that specific mutations inferred from an analysis of the simulated folding pathway exacerbate its stability.

Improved transition path sampling methods for simulation of rare events.

The free energy surfaces of a wide variety of systems encountered in physics, chemistry, and biology are characterized by the existence of deep minima separated by numerous barriers. One of the central aims of recent research in computational chemistry and physics has been to determine how transitions occur between deep local minima on rugged free energy landscapes, and transition path sampling (TPS) Monte-Carlo methods have emerged as an effective means for numerical investigation of such transitions. Many of the shortcomings of TPS-like approaches generally stem from their high computational demands. Two new algorithms are presented in this work that improve the efficiency of TPS simulations. The first algorithm uses biased shooting moves to render the sampling of reactive trajectories more efficient. The second algorithm is shown to substantially improve the accuracy of the transition state ensemble by introducing a subset of local transition path simulations in the transition state. The system considered in this work consists of a two-dimensional rough energy surface that is representative of numerous systems encountered in applications. When taken together, these algorithms provide gains in efficiency of over two orders of magnitude when compared to traditional TPS simulations.