On the origins of the weak folding cooperativity of a designed ββα ultrafast protein FSD-1

Hetero-oligomeric Amyloid Assembly and Mechanism: Prion Fragment PrP(106-126) Catalyzes the Islet Amyloid Polypeptide β-Hairpin

Protein aggregation is typically attributed to the association of homologous amino acid sequences between monomers of the same protein. Coaggregation of heterogeneous peptide species can occur, however, and is implicated in the proliferation of seemingly unrelated protein diseases in the body. The prion protein fragment (PrP_106-126) and human islet amyloid polypeptide (hIAPP) serve as an interesting model of nonhomologous protein assembly as they coaggregate, despite a lack of sequence homology. We have applied ion-mobility mass spectrometry, atomic force microscopy, circular dichroism, and high-level molecular modeling to elucidate this important assembly process. We found that the prion fragment not only forms pervasive hetero-oligomeric aggregates with hIAPP but also promotes the transition of hIAPP into its amyloidogenic β-hairpin conformation. Further, when PrP_106-126 was combined with non-amyloidogenic rIAPP, the two formed nearly identical hetero-oligomers to those seen with hIAPP, despite rIAPP containing β-sheet breaking proline substitutions. Additionally, while rIAPP does not natively form the amyloidogenic β-hairpin structure, it did so in the presence of PrP_106-126 and underwent a conformational transition to β-sheet in solution. We also find that PrP_106-126 forms hetero-oligomers with the IAPP_8-20 fragment but not with the "aggregation hot spot" IAPP_20-29 fragment. PrP_106-126 apparently induces IAPP into a β-hairpin structure within the PrP:IAPP heterodimer complex and then, through ligand exchange, catalytically creates the amyloidogenic β-hairpin dimer of IAPP in significantly greater abundance than IAPP does on its own. This is a new mechanistic model that provides a critical foundation for the detailed study of hetero-oligomerization and prion-like proliferation in amyloid systems.

Validation of Molecular Dynamics Simulations for Prediction of Three-Dimensional Structures of Small Proteins

Although various higher-order protein structure prediction methods have been developed, almost all of them were developed based on the three-dimensional (3D) structure information of known proteins. Here we predicted the short protein structures by molecular dynamics (MD) simulations in which only Newton’s equations of motion were used and 3D structural information of known proteins was not required. To evaluate the ability of MD simulationto predict protein structures, we calculated seven short test protein (10–46 residues) in the denatured state and compared their predicted and experimental structures. The predicted structure for Trp-cage (20 residues) was close to the experimental structure by 200-ns MD simulation. For proteins shorter or longer than Trp-cage, root-mean square deviation values were larger than those for Trp-cage. However, secondary structures could be reproduced by MD simulations for proteins with 10–34 residues. Simulations by replica exchange MD were performed, but the results were similar to those from normal MD simulations. These results suggest that normal MD simulations can roughly predict short protein structures and 200-ns simulations are frequently sufficient for estimating the secondary structures of protein (approximately 20 residues). Structural prediction method using only fundamental physical laws are useful for investigating non-natural proteins, such as primitive proteins and artificial proteins for peptide-based drug delivery systems.

Finding multiple reaction pathways via global optimization of action

Global searching for reaction pathways is a long-standing challenge in computational chemistry and biology. Most existing approaches perform only local searches due to computational complexity. Here we present a computational approach, Action-CSA, to find multiple diverse reaction pathways connecting fixed initial and final states through global optimization of the Onsager-Machlup action using the conformational space annealing (CSA) method. Action-CSA successfully overcomes large energy barriers via crossovers and mutations of pathways and finds all possible pathways of small systems without initial guesses on pathways. The rank order and the transition time distribution of multiple pathways are in good agreement with those of long Langevin dynamics simulations. The lowest action folding pathway of FSD-1 is consistent with recent experiments. The results show that Action-CSA is an efficient and robust computational approach to study the multiple pathways of complex reactions and large-scale conformational changes.

Folding of a Zinc-Finger ββα-Motif Investigated Using Two-Dimensional and Time-Resolved Vibrational Spectroscopy

Small proteins provide good model systems for studying the fundamental forces that control protein folding. Here, we investigate the folding dynamics of the 28-residue zinc-finger mutant FSD-1, which is designed to form a metal-independent folded ββα-motif, and which provides a testing ground for proteins containing a mixed α/β fold. Although the folding of FSD-1 has been actively studied, the folding mechanism remains largely unclear. In particular, it is unclear in what stage of folding the α-helix is formed. To address this issue we investigate the folding mechanism of FSD-1 using a combination of temperature-dependent UV circular dichroism (UV-CD), Fourier transform infrared (FTIR) spectroscopy, two-dimensional infrared (2D-IR) spectroscopy, and temperature-jump (T-jump) transient-IR spectroscopy. Our UV-CD and FTIR data show different thermal melting transitions, indicating multistate folding behavior. Temperature-dependent 2D-IR spectra indicate that the α-helix is the most stable structural element of FSD-1. To investigate the folding/unfolding re-equilibration dynamics of FSD-1, the conformational changes induced by a nanosecond T-jump are probed with transient-IR and transient dispersed-pump-probe (DPP) IR spectroscopy. We observe biexponential T-jump relaxation kinetics (with time constants of 80 ± 13 ns and 1300 ± 100 ns at 322 K), confirming that the folding involves an intermediate state. The IR and dispersed-pump-probe IR spectra associated with the two kinetic components suggest that the folding of FSD-1 involves early formation of the α-helix, followed by the formation of the β-hairpin and hydrophobic contacts.

The Folding of de Novo Designed Protein DS119 via Molecular Dynamics Simulations

As they are not subjected to natural selection process, de novo designed proteins usually fold in a manner different from natural proteins. Recently, a de novo designed mini-protein DS119, with a βαβ motif and 36 amino acids, has folded unusually slowly in experiments, and transient dimers have been detected in the folding process. Here, by means of all-atom replica exchange molecular dynamics (REMD) simulations, several comparably stable intermediate states were observed on the folding free-energy landscape of DS119. Conventional molecular dynamics (CMD) simulations showed that when two unfolded DS119 proteins bound together, most binding sites of dimeric aggregates were located at the N-terminal segment, especially residues 5-10, which were supposed to form β-sheet with its own C-terminal segment. Furthermore, a large percentage of individual proteins in the dimeric aggregates adopted conformations similar to those in the intermediate states observed in REMD simulations. These results indicate that, during the folding process, DS119 can easily become trapped in intermediate states. Then, with diffusion, a transient dimer would be formed and stabilized with the binding interface located at N-terminals. This means that it could not quickly fold to the native structure. The complicated folding manner of DS119 implies the important influence of natural selection on protein-folding kinetics, and more improvement should be achieved in rational protein design.

Self-aggregation and coaggregation of the p53 core fragment with its aggregation gatekeeper variant

Recent studies suggested that p53 aggregation can lead to loss-of-function (LoF), dominant-negative (DN) and gain-of-function (GoF) effects, with adverse cancer consequences. The p53 aggregation-nucleating (251)ILTIITL(257) fragment is a key segment in wild-type p53 aggregation; however, an I254R mutation can prevent it. It was suggested that self-assembly of wild-type p53 and its cross-interaction with mutants differ from the classical amyloid nucleation-growth mechanism. Here, using replica exchange molecular dynamics (REMD) simulations, we studied the cross-interactions of this p53 core fragment and its aggregation rescue I254R mutant. We found that the core fragment displays strong aggregation propensity, whereas the gatekeeper I254R mutant tends to be disordered, consistent with experiments. Our cross-interaction results reveal that the wild-type p53 fragment promotes β-sheet formation of the I254R mutant by shifting the disordered mutant peptides into aggregating states. As a result, the system has similar oligomeric structures, inter-peptide interactions and free energy landscape as the wild type fragment does, revealing a prion-like process. We also found that in the cross-interaction system, the wild-type species has higher tendency to interact with the mutant than with itself. This phenomenon illustrates synergistic effects between the p53 (251)ILTIITL(257) fragment and the mutant resembling prion cross-species propagation, cautioning against exploiting it in drug discovery.

Polymorphic cross-seeding amyloid assemblies of amyloid-β and human islet amyloid polypeptide

Epidemiological studies have shown that the development of Alzheimer's disease (AD) is associated with type 2 diabetes (T2D), but it still remains unclear how AD and T2D are connected. Heterologous cross-seeding between the causative peptides of Aβ and hIAPP may represent a molecular link between AD and T2D. Here, we computationally modeled and simulated a series of cross-seeding double-layer assemblies formed by Aβ and hIAPP peptides using all-atom and coarse-gained molecular dynamics (MD) simulations. The cross-seeding Aβ-hIAPP assemblies showed a wide range of polymorphic structures via a combination of four β-sheet-to-β-sheet interfaces and two packing orientations, focusing on a comparison of different matches of β-sheet layers. Two cross-seeding Aβ-hIAPP assemblies with different interfacial β-sheet packings exhibited high structural stability and favorable interfacial interactions in both oligomeric and fibrillar states. Both Aβ-hIAPP assemblies displayed interfacial dehydration to different extents, which in turn promoted Aβ-hIAPP association depending on interfacial polarity and geometry. Furthermore, computational mutagenesis studies revealed that disruption of interfacial salt bridges largely disfavor the β-sheet-to-β-sheet association, highlighting the importance of salt bridges in the formation of cross-seeding assemblies. This work provides atomic-level information on the cross-seeding interactions between Aβ and hIAPP, which may be involved in the interplay between these two disorders.

Defining the molecular basis of amyloid inhibitors: human islet amyloid polypeptide-insulin interactions

Human islet amyloid polypeptide (hIAPP or Amylin) is a 37 residue hormone that is cosecreted with insulin from the pancreatic islets. The aggregation of hIAPP plays a role in the progression of type 2 diabetes and contributes to the failure of islet cell grafts. Despite considerable effort, little is known about the mode of action of IAPP amyloid inhibitors, and this has limited rational drug design. Insulin is one of the most potent inhibitors of hIAPP fibril formation, but its inhibition mechanism is not understood. In this study, the aggregation of mixtures of hIAPP with insulin, as well as with the separate A and B chains of insulin, were characterized using ion mobility spectrometry-based mass spectrometry and atomic force microscopy. Insulin and the insulin B chain target the hIAPP monomer in its compact isoform and shift the equilibrium away from its extended isoform, an aggregation-prone conformation, and thus inhibit hIAPP from forming β-sheets and subsequently amyloid fibrils. All-atom molecular modeling supports these conclusions.

Structural similarities and differences between amyloidogenic and non-amyloidogenic islet amyloid polypeptide (IAPP) sequences and implications for the dual physiological and pathological activities of these peptides

IAPP, a 37 amino-acid peptide hormone belonging to the calcitonin family, is an intrinsically disordered protein that is coexpressed and cosecreted along with insulin by pancreatic islet β-cells in response to meals. IAPP plays a physiological role in glucose regulation; however, in certain species, IAPP can aggregate and this process is linked to β-cell death and Type II Diabetes. Using replica exchange molecular dynamics with extensive sampling (16 replicas per sequence and 600 ns per replica), we investigate the structure of the monomeric state of two species of aggregating peptides (human and cat IAPP) and two species of non-aggregating peptides (pig and rat IAPP). Our simulations reveal that the pig and rat conformations are very similar, and consist of helix-coil and helix-hairpin conformations. The aggregating sequences, on the other hand, populate the same helix-coil and helix-hairpin conformations as the non-aggregating sequence, but, in addition, populate a hairpin structure. Our exhaustive simulations, coupled with available peptide-activity data, leads us to a structure-activity relationship (SAR) in which we propose that the functional role of IAPP is carried out by the helix-coil conformation, a structure common to both aggregating and non-aggregating species. The pathological role of this peptide may have multiple origins, including the interaction of the helical elements with membranes. Nonetheless, our simulations suggest that the hairpin structure, only observed in the aggregating species, might be linked to the pathological role of this peptide, either as a direct precursor to amyloid fibrils, or as part of a cylindrin type of toxic oligomer. We further propose that the helix-hairpin fold is also a possible aggregation prone conformation that would lead normally non-aggregating variants of IAPP to form fibrils under conditions where an external perturbation is applied. The SAR relationship is used to suggest the rational design of therapeutics for treating diabetes.

IAPP, a 37 amino-acid peptide hormone belonging to the calcitonin family, is an intrinsically disordered peptide produced along with insulin by pancreatic islet β-cells in response to meals. In its functional form, IAPP acts as a synergic partner of insulin to reduce blood glucose. IAPP can, however, also play a pathological role, contributing to Type II diabetes (T2D). Knowledge of the structural nature of the physiological and pathological forms of IAPP will facilitate the rational design of novel drugs for therapeutic treatment of T2D. However, because IAPP does not fold to a single structure, but rather co-exists between multiple functional (and toxic) structures, it is extremely challenging for experimental methods to gain detailed structural information. Using a computational approach, we were able to obtain detailed structures of four IAPP variants and propose a novel structural hypothesis for the two opposing roles of this peptide.

Folding models of mini-protein FSD-1

We have carried out all-atom action-derived molecular dynamics (ADMD) folding simulations of the full-size FSD-1. FSD-1 is a designed mini-protein of 28 residues containing both α and β secondary structure elements. Multiple folding pathways are found for FSD-1, which is consistent with existing computational studies. Hydrophobic collapse is observed first, and then subsequent folding events proceeds by forming either α-helix or β-hairpin. Concurrent formation of the full tertiary structure and the secondary structure elements of α-helix and β-hairpin is observed. The folding pathway of FSD-1 elucidated by ADMD simulations does not follow the scenario of the framework model. ADMD simulations provide significant insights for the general mechanisms of protein folding and conformational changes.

The amyloid formation mechanism in human IAPP: dimers have β-strand monomer-monomer interfaces

Early oligomerization of human IAPP (hIAPP) is responsible for β-cell death in the pancreas and is increasingly considered a primary pathological process linked to Type II Diabetes (T2D). Yet, the assembly mechanism remains poorly understood, largely due to the inability of conventional techniques to probe distributions or detailed structures of early oligomeric species. Here, we describe the first experimental data on the isolated and unmodified dimers of human (hIAPP) and nonamyloidogenic rat IAPP (rIAPP). The experiments reveal that the human IAPP dimers are more extended than those formed by rat IAPP and likely descend from extended monomers. Independent all-atom molecular dynamics simulations show that rIAPP forms compact helix and coil rich dimers, whereas hIAPP forms β-strand rich dimers that are generally more extended. Also, the simulations reveal that the monomer-monomer interfaces of the hIAPP dimers are dominated by β-strands and that β-strands can recruit coil or helix structured regions during the dimerization process. Our β-rich interface contrasts with an N-terminal helix-to-helix interface proposed in the literature but is consistent with existing experimental data on the self-interaction pattern of hIAPP, mutation effects, and inhibition effects of the N-methylation in the mutation region.