Chi Cheng Chiu

Two-dimensional infrared spectroscopy reveals the complex behaviour of an amyloid fibril inhibitor

While amyloid formation has been implicated in the pathology of over twenty human diseases, the rational design of amyloid inhibitors is hampered by a lack of structural information about amyloid-inhibitor complexes. We use isotope labeling and two-dimensional infrared spectroscopy to obtain a residue-specific structure for the complex of human amylin, the peptide responsible for islet amyloid formation in type 2 diabetes, with a known inhibitor, rat amylin. Based on its sequence, rat amylin should block formation of the C-terminal β-sheet, but at 8 hours after mixing rat amylin blocks the N-terminal β-sheet instead. At 24 hours after mixing, rat amylin blocks neither β-sheet and forms its own β-sheet most likely on the outside of the human fibrils. This is striking because rat amylin is natively disordered and not previously known to form amyloid β-sheets. The results show that even seemingly intuitive inhibitors may function by unforeseen and complex structural processes.

Mechanism of IAPP amyloid fibril formation involves an intermediate with a transient β-sheet

Significance There is an enormous interest in the mechanism by which proteins misfold and aggregate into amyloid fibrils. Amyloid has been implicated in many human diseases, but the mechanism of aggregation is not understood. Intermediates have been postulated to play an important role in the process, but there have been very few direct measurements that provide specific structural details. The use of isotope labeling and 2D IR methods has allowed the characterization of a critical intermediate generated during amyloid formation by islet amyloid polypeptide, the peptide responsible for amyloid formation in type 2 diabetes. Identification of this intermediate provides a structural explanation for the lag phase and may explain why some species develop amyloid deposits of hIAPP while others do not. Amyloid formation is implicated in more than 20 human diseases, yet the mechanism by which fibrils form is not well understood. We use 2D infrared spectroscopy and isotope labeling to monitor the kinetics of fibril formation by human islet amyloid polypeptide (hIAPP or amylin) that is associated with type 2 diabetes. We find that an oligomeric intermediate forms during the lag phase with parallel β-sheet structure in a region that is ultimately a partially disordered loop in the fibril. We confirm the presence of this intermediate, using a set of homologous macrocyclic peptides designed to recognize β-sheets. Mutations and molecular dynamics simulations indicate that the intermediate is on pathway. Disrupting the oligomeric β-sheet to form the partially disordered loop of the fibrils creates a free energy barrier that is the origin of the lag phase during aggregation. These results help rationalize a wide range of previous fragment and mutation studies including mutations in other species that prevent the formation of amyloid plaques.

Parallel β-sheet vibrational couplings revealed by 2D IR spectroscopy of an isotopically labeled macrocycle: quantitative benchmark for the interpretation of amyloid and protein infrared spectra.

Infrared spectroscopy is playing an important role in the elucidation of amyloid fiber formation, but the coupling models that link spectra to structure are not well tested for parallel β-sheets. Using a synthetic macrocycle that enforces a two stranded parallel β-sheet conformation, we measured the lifetimes and frequency for six combinations of doubly (13)C═(18)O labeled amide I modes using 2D IR spectroscopy. The average vibrational lifetime of the isotope labeled residues was 550 fs. The frequencies of the labels ranged from 1585 to 1595 cm(-1), with the largest frequency shift occurring for in-register amino acids. The 2D IR spectra of the coupled isotope labels were calculated from molecular dynamics simulations of a series of macrocycle structures generated from replica exchange dynamics to fully sample the conformational distribution. The models used to simulate the spectra include through-space coupling, through-bond coupling, and local frequency shifts caused by environment electrostatics and hydrogen bonding. The calculated spectra predict the line widths and frequencies nearly quantitatively. Historically, the characteristic features of β-sheet infrared spectra have been attributed to through-space couplings such as transition dipole coupling. We find that frequency shifts of the local carbonyl groups due to nearest neighbor couplings and environmental factors are more important, while the through-space couplings dictate the spectral intensities. As a result, the characteristic absorption spectra empirically used for decades to assign parallel β-sheet secondary structure arises because of a redistribution of oscillator strength, but the through-space couplings do not themselves dramatically alter the frequency distribution of eigenstates much more than already exists in random coil structures. Moreover, solvent exposed residues have amide I bands with >20 cm(-1) line width. Narrower line widths indicate that the amide I backbone is solvent protected inside the macrocycle. This work provides calculated and experimentally verified couplings for parallel β-sheets that can be used in structure-based models to simulate and interpret the infrared spectra of β-sheet containing proteins and protein assemblies, such as amyloid fibers.

Molecular dynamics study of a nanotube-binding amphiphilic helical peptide at different water/hydrophobic interfaces.

Many potential applications of single-walled carbon nanotubes (SWNTs) require that they be isolated from one another. This may be accomplished through covalent or noncovalent SWNT functionalization. The noncovalent approach preserves the intrinsic electrical, optical, and mechanical properties of SWNTs and can be achieved by dispersing SWNTs in aqueous solution using surfactants, polymers, or biomacromolecules like DNA or polypeptides. The designed amphiphilic helical peptide nano-1, which contains hydrophobic valine and aromatic phenylalanine residues for interaction with SWNTs and glutamic acid and lysine residues for water solubility, has been shown to debundle and disperse SWNTs, although the details of the peptide-SWNT interactions await elucidation. Here we use fully atomistic molecular dynamics simulations to investigate the nano-1 peptide at three different water/hydrophobic interfaces: water/oil, water/graphite, and water/SWNT. The amphiphilic nature of the peptide is characterized by its secondary structure, peptide-water hydrogen bonding, and peptide-hydrophobic surface van der Waals energy. We show that nano-1 has reduced amphiphilic character at the water/oil interface because the peptide helix penetrates into the hydrophobic phase. The peptide alpha-helix cannot match its hydrophobic face to the rigid planar graphite surface without partially unfolding. In contrast, nano-1 can curve on the SWNT surface in an alpha-helical conformation to simultaneously maximize its hydrophobic contacts with the SWNT and its hydrogen bonds with water. The molecular insight into the peptide conformation at the various hydrophobic surfaces provides guidelines for future peptide design.

Coarse-grained potential models for phenyl-based molecules: I. Parametrization using experimental data.

A coarse-grained intermolecular potential has been parametrized for phenyl-based molecules. The parametrization was accomplished by fitting to experimental thermodynamic data. Specifically, the intermolecular potentials, which were based on Lennard-Jones functional forms, were parametrized and validated using experimental surface tension, density, and partitioning data. This approach has been used herein to develop parameters for coarse-grained interaction sites that are applicable to a variety of phenyl-based molecules, including analogues of the amino acid side chains of phenylalanine and tyrosine. Comparison of the resulting coarse-grain model to atomistic simulations shows a high level of structural and thermodynamic agreement between the two models, despite the fact that no atomistic simulation data was used in the parametrization of the coarse-grain intermolecular potentials.

Effect of proline mutations on the monomer conformations of amylin

The formation of human islet amyloid polypeptide (hIAPP) is implicated in the loss of pancreatic β-cells in type II diabetes. Rat amylin, which differs from human amylin at six residues, does not lead to formation of amyloid fibrils. Pramlintide is a synthetic analog of human amylin that shares three proline substitutions with rat amylin. Pramlintide has a much smaller propensity to form amyloid aggregates and has been widely prescribed in amylin replacement treatment. It is known that the three prolines attenuate β-sheet formation. However, the detailed effects of these proline substitutions on full-length hIAPP remain poorly understood. In this work, we use molecular simulations and bias-exchange metadynamics to investigate the effect of proline substitutions on the conformation of the hIAPP monomer. Our results demonstrate that hIAPP can adopt various β-sheet conformations, some of which have been reported in experiments. The proline substitutions perturb the formation of long β-sheets and reduce their stability. More importantly, we find that all three proline substitutions of pramlintide are required to inhibit β conformations and stabilize the α-helical conformation. Fewer substitutions do not have a significant inhibiting effect.

Parametrization and Application of a Coarse Grained Force Field for Benzene/Fullerene Interactions with Lipids

Recently, we reported new coarse grain (CG) force fields for lipids and phenyl/fullerene based molecules. Here, we developed the cross parameters necessary to unite those force fields and then applied the model to investigate the nature of benzene and C(60) interactions with lipid bilayers. The interaction parameters between the phenyl and lipid CG sites are based on experimental and all atom (AA) molecular dynamics (MD) data. The resulting force field was tested on benzene rich lipid bilayers and shown to reproduce general behavior expected from experiments. The parameters were then applied to C(60) interactions with lipid bilayers. Overall, the results showed excellent agreement with AA MD and experimental observations. In the C(60) lipid systems, the fullerenes were shown to aggregate even at the lowest concentrations investigated.

Coarse-grained potential models for phenyl-based molecules: II. Application to fullerenes.

The interaction of fullerenes with biological systems and the environment is a topic of current interest. Coarse-grained molecular dynamics (CGMD) simulations are well-suited to investigate some of the factors involved because they provide access to time and length scales that are not accessible using fully atomistic simulation methods. In the case of buckyballs (C(60)) and single-walled carbon nanotubes (SWNTs), it is necessary to parametrize a CG force field that accurately captures the balance between fullerene-fullerene and fullerene-solvent interactions. Herein, we derive CG force field parameters for C(60) and SWNTs by using the optimized benzene parameters from part I [DeVane, R.; Chiu, C.-c.; Nielsen, S. O.; Shinoda, W.; Moore, P. B.; Klein, M. L. J. Phys. Chem. B 2010, doi: 10.1021/jp9117369]. Solubility, transfer free energy, and dimerization free-energy data for C(60) and SWNTs obtained using the proposed models show excellent agreement with experimental and fully atomistic MD data. In particular, cluster analysis of C(60) aggregation in a hydrocarbon melt corroborates the force field parameters. The aggregation behavior of the present CG force field differs considerably from that of models currently in widespread use. The combined results provide a strong basis for applying this model for further large-scale MD studies involving C(60) and SWNTs.

Role of peptide-peptide interactions in stabilizing peptide-wrapped single-walled carbon nanotubes: A molecular dynamics study

Single‐walled carbon nanotubes (SWNTs) have unique properties and are projected to have a major impact in nanoscale electronics, materials science, and nanomedicine. Yet, these potential applications are hindered by the need for sample purification to separate SWNTs from each other and from metallic catalyst and amorphous carbon present in as‐synthesized samples. Common purification strategies involve dispersing SWNTs as individual tubes in aqueous solution. Towards this end, a designed helical peptide was shown to be excellent at dispersing SWNTs. However, the molecular details of the peptide‐SWNT and peptide‐peptide interactions await elucidation. Here we explore these molecular interactions using fully atomistic molecular dynamics simulations of peptide‐wrapped SWNTs. We characterize the interactions by measuring the aromatic residue‐to‐SWNT surface distance, the peptide amphiphilicity, the peptide‐SWNT crossing angle, the peptide‐SWNT contact area, the peptide helix axis‐to‐axis distance, and the inter‐peptide hydrogen bonding. We find that the peptides collectively tilt with respect to the SWNT long axis, are α‐helical, and form interpeptide hydrogen bonds through their lysine and glutamate residues, which helps to stabilize the multipeptide/SWNT complex. All hydrophobic residues interact with the SWNT and are sequestered from water. The picture that emerges from this study gives insight into subsequent peptide design.

α-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.