Spinal motor neuron protein supersaturation patterns are associated with inclusion body formation in ALS.

Amyotrophic lateral sclerosis (ALS) is a heterogeneous degenerative motor neuron disease linked to numerous genetic mutations in apparently unrelated proteins. These proteins, including SOD1, TDP-43, and FUS, are highly aggregation-prone and form a variety of intracellular inclusion bodies that are characteristic of different neuropathological subtypes of the disease. Contained within these inclusions are a variety of proteins that do not share obvious characteristics other than coaggregation. However, recent evidence from other neurodegenerative disorders suggests that disease-affected biochemical pathways can be characterized by the presence of proteins that are supersaturated, with cellular concentrations significantly greater than their solubilities. Here, we show that the proteins that form inclusions of mutant SOD1, TDP-43, and FUS are not merely a subset of the native interaction partners of these three proteins, which are themselves supersaturated. To explain the presence of coaggregating proteins in inclusions in the brain and spinal cord, we observe that they have an average supersaturation even greater than the average supersaturation of the native interaction partners in motor neurons, but not when scores are generated from an average of other human tissues. These results suggest that inclusion bodies in various forms of ALS result from a set of proteins that are metastable in motor neurons, and thus prone to aggregation upon a disease-related progressive collapse of protein homeostasis in this specific setting.

Expression in drosophila of tandem amyloid ß peptides provides insights into links between aggregation and neurotoxicity.

The generation and subsequent aggregation of amyloid ß (Aß) peptides play a crucial initiating role in the pathogenesis of Alzheimer disease (AD). The two main isoforms of these peptides have 40 (Aß(40)) or 42 residues (Aß(42)), the latter having a higher propensity to aggregate in vitro and being the main component of the plaques observed in vivo in AD patients. We have designed a series of tandem dimeric constructs of these Aß peptides to probe the manner in which changes in the aggregation kinetics of Aß affect its deposition and toxicity in a Drosophila melanogaster model system. The levels of insoluble aggregates were found to be substantially elevated in flies expressing the tandem constructs of both Aß(40) and Aß(42) compared with the equivalent monomeric peptides, consistent with the higher effective concentration, and hence increased aggregation rate, of the peptides in the tandem repeat. A unique feature of the Aß(42) constructs, however, is the appearance of high levels of soluble oligomeric aggregates and a corresponding dramatic increase in their in vivo toxicity. The toxic nature of the Aß(42) peptide in vivo can therefore be attributed to the higher kinetic stability of the oligomeric intermediate states that it populates relative to those of Aß(40) rather than simply to its higher rate of aggregation.

Intrinsic determinants of neurotoxic aggregate formation by the amyloid beta peptide.

The extent to which proteins aggregate into distinct structures ranging from prefibrillar oligomers to amyloid fibrils is key to the pathogenesis of many age-related degenerative diseases. We describe here for the Alzheimer's disease-related amyloid beta peptide (Abeta) an investigation of the sequence-based determinants of the balance between the formation of prefibrillar aggregates and amyloid fibrils. We show that by introducing single-point mutations, it is possible to convert the normally harmless Abeta40 peptide into a pathogenic species by increasing its relative propensity to form prefibrillar but not fibrillar aggregates, and, conversely, to abolish the pathogenicity of the highly neurotoxic E22G Abeta42 peptide by reducing its relative propensity to form prefibrillar species rather than mature fibrillar ones. This observation can be rationalized by the demonstration that whereas regions of the sequence of high aggregation propensity dominate the overall tendency to aggregate, regions with low intrinsic aggregation propensities exert significant control over the balance of the prefibrillar and fibrillar species formed, and therefore play a major role in determining the neurotoxicity of the Abeta peptide.

Correction: esyN: Network Building, Sharing and Publishing.

[This corrects the article DOI: 10.1371/journal.pone.0106035.].

Calculation of the free energy barriers in the oligomerisation of Abeta peptide fragments.

Protein misfolding and aggregation are associated with a range of severe human neurodegenerative conditions. We use all-atom simulations to describe the process of assembly of the Abeta(16-22) and Abeta(25-35) fragments of Abeta, a peptide associated with Alzheimer's disease. Our results indicate that the pathways of aggregation of these two peptides depend predominantly on the relative strength of hydrophobic interactions and hydrogen bonding. In the Abeta(25-35) peptide, which is weakly hydrophobic, the tendency to form hydrogen bonds drives the crossing of a single major free energy barrier for the formation of a cross-beta structure. By contrast, in the more hydrophobic Abeta(16-22) peptide, the process of ordered assembly is preceded by an initial collapse into disordered oligomers. These results provide support for a recently proposed two-step mechanism of amyloid formation. We have also found that the barriers for reordering are lower for large oligomers than for small oligomers, a result that provides an explanation of the recent experimental observation that the efficiency of the seeding reaction depends on the size of the seeds themselves.

Structural reorganisation and potential toxicity of oligomeric species formed during the assembly of amyloid fibrils.

Increasing evidence indicates that oligomeric protein assemblies may represent the molecular species responsible for cytotoxicity in a range of neurological disorders including Alzheimer and Parkinson diseases. We use all-atom computer simulations to reveal that the process of oligomerization can be divided into two steps. The first is characterised by a hydrophobic coalescence resulting in the formation of molten oligomers in which hydrophobic residues are sequestered away from the solvent. In the second step, the oligomers undergo a process of reorganisation driven by interchain hydrogen bonding interactions that induce the formation of beta sheet rich assemblies in which hydrophobic groups can become exposed. Our results show that the process of aggregation into either ordered or amorphous species is largely determined by a competition between the hydrophobicity of the amino acid sequence and the tendency of polypeptide chains to form arrays of hydrogen bonds. We discuss how the increase in solvent-exposed hydrophobic surface resulting from such a competition offers an explanation for recent observations concerning the cytotoxicity of oligomeric species formed prior to mature amyloid fibrils.

Cholinergic neuron gene expression differences captured by translational profiling in a mouse model of Alzheimer's disease.

Cholinergic neurotransmission is impaired in Alzheimer's disease (AD), and loss of basal forebrain cholinergic neurons is a key component of disease pathogenicity and symptomatology. To explore the molecular basis of this cholinergic dysfunction, we paired translating ribosome affinity purification (TRAP) with RNA sequencing (TRAP-Seq) to identify the actively translating mRNAs in anterior forebrain cholinergic neurons in the TgCRND8 mouse model of AD. Bioinformatic analyses revealed the downregulation of 67 of 71 known cholinergic-related transcripts, consistent with cholinergic neuron dysfunction in TgCRND8 mice, as well as transcripts related to oxidative phosphorylation, neurotrophins, and ribosomal processing. Upregulated transcripts included those related to axon guidance, glutamatergic synapses and kinase activity and included AD-risk genes Sorl1 and Ptk2b. In contrast, the total transcriptome of the anterior forebrain showed upregulation in cytokine signaling, microglia, and immune system pathways, including Trem2, Tyrobp, and Inpp5d. Hence, TRAP-Seq clearly distinguished the differential gene expression alterations occurring in cholinergic neurons of TgCRND8 mice compared with wild-type littermates, providing novel candidate pathways to explore for therapeutic development in AD.

Network Approaches to the Understanding of Alzheimer's Disease: From Model Organisms to Humans.

It is becoming increasingly evident that Alzheimer's disease cannot be considered as the outcome of a single pathway, but rather we should view it as a system, that is, a network of interactions between large numbers of different protein molecules. In the last few years, probably because of the inherent limitations of traditional methods and because of the great increase in availability of sequencing data, this type of approach is being used more and more. In the following, we will discuss what constitutes a "network approach," what are its pros and cons, a number of recent case studies and finally what are the future perspectives of this type of analysis.

ALS/FTD Mutation-Induced Phase Transition of FUS Liquid Droplets and Reversible Hydrogels into Irreversible Hydrogels Impairs RNP Granule Function.

The mechanisms by which mutations in FUS and other RNA binding proteins cause ALS and FTD remain controversial. We propose a model in which low-complexity (LC) domains of FUS drive its physiologically reversible assembly into membrane-free, liquid droplet and hydrogel-like structures. ALS/FTD mutations in LC or non-LC domains induce further phase transition into poorly soluble fibrillar hydrogels distinct from conventional amyloids. These assemblies are necessary and sufficient for neurotoxicity in a C. elegans model of FUS-dependent neurodegeneration. They trap other ribonucleoprotein (RNP) granule components and disrupt RNP granule function. One consequence is impairment of new protein synthesis by cytoplasmic RNP granules in axon terminals, where RNP granules regulate local RNA metabolism and translation. Nuclear FUS granules may be similarly affected. Inhibiting formation of these fibrillar hydrogel assemblies mitigates neurotoxicity and suggests a potential therapeutic strategy that may also be applicable to ALS/FTD associated with mutations in other RNA binding proteins.

Folding of a small helical protein using hydrogen bonds and hydrophobicity forces.

A reduced protein model with five to six atoms per amino acid and five amino acid types is developed and tested on a three-helix-bundle protein, a 46-amino acid fragment from staphylococcal protein A. The model does not rely on the widely used Go approximation, which ignores non-native interactions. We find that the collapse transition is considerably more abrupt for the protein A sequence than for random sequences with the same composition. The chain collapse is found to be at least as fast as helix formation. Energy minimization restricted to the thermodynamically favored topology gives a structure that has a root-mean-square deviation of 1.8 A from the native structure. The sequence-dependent part of our potential is pairwise additive. Our calculations suggest that fine-tuning this potential by parameter optimization is of limited use.

Sequence-based study of two related proteins with different folding behaviors.

Z(SPA-1) is an engineered protein that binds to its parent, the three-helix-bundle Z domain of staphylococcal protein A. Uncomplexed Z(SPA-1) shows a reduced helix content and a melting behavior that is less cooperative, compared with the wild-type Z domain. Here we show that the difference in folding behavior between these two sequences can be partly understood in terms of an off-lattice model with 5-6 atoms per amino acid and a minimalistic potential, in which folding is driven by backbone hydrogen bonding and effective hydrophobic attraction.

esyN: network building, sharing and publishing.

The construction and analysis of networks is increasingly widespread in biological research. We have developed esyN ("easy networks") as a free and open source tool to facilitate the exchange of biological network models between researchers. esyN acts as a searchable database of user-created networks from any field. We have developed a simple companion web tool that enables users to view and edit networks using data from publicly available databases. Both normal interaction networks (graphs) and Petri nets can be created. In addition to its basic tools, esyN contains a number of logical templates that can be used to create models more easily. The ability to use previously published models as building blocks makes esyN a powerful tool for the construction of models and network graphs. Users are able to save their own projects online and share them either publicly or with a list of collaborators. The latter can be given the ability to edit the network themselves, allowing online collaboration on network construction. esyN is designed to facilitate unrestricted exchange of this increasingly important type of biological information. Ultimately, the aim of esyN is to bring the advantages of Open Source software development to the construction of biological networks.

The TRiC/CCT chaperone is implicated in Alzheimer's disease based on patient GWAS and an RNAi screen in Aß-expressing Caenorhabditis elegans.

The human Aß peptide causes progressive paralysis when expressed in the muscles of the nematode worm, C. elegans. We have exploited this model of Aß toxicity by carrying out an RNAi screen to identify genes whose reduced expression modifies the severity of this locomotor phenotype. Our initial finding was that none of the human orthologues of these worm genes is identical with the genome-wide significant GWAS genes reported to date (the "white zone"); moreover there was no identity between worm screen hits and the longer list of GWAS genes which included those with borderline levels of significance (the "grey zone"). This indicates that Aß toxicity should not be considered as equivalent to sporadic AD. To increase the sensitivity of our analysis, we then considered the physical interactors (+1 interactome) of the products of the genes in both the worm and the white+grey zone lists. When we consider these worm and GWAS gene lists we find that 4 of the 60 worm genes have a +1 interactome overlap that is larger than expected by chance. Two of these genes form a chaperonin complex, the third is closely associated with this complex and the fourth gene codes for actin, the major substrate of the same chaperonin.

ANS binding reveals common features of cytotoxic amyloid species.

Oligomeric assemblies formed from a variety of disease-associated peptides and proteins have been strongly associated with toxicity in many neurodegenerative conditions, such as Alzheimer's disease. The precise nature of the toxic agents, however, remains still to be established. We show that prefibrillar aggregates of E22G (arctic) variant of the Abeta(1-42) peptide bind strongly to 1-anilinonaphthalene 8-sulfonate and that changes in this property correlate significantly with changes in its cytotoxicity. Moreover, we show that this phenomenon is common to other amyloid systems, such as wild-type Abeta(1-42), the I59T variant of human lysozyme and an SH3 domain. These findings are consistent with a model in which the exposure of hydrophobic surfaces as a result of the aggregation of misfolded species is a crucial and common feature of these pathogenic species.

Finite size effects in simulations of protein aggregation.

It is becoming increasingly clear that the soluble protofibrillar species that proceed amyloid fibril formation are associated with a range of neurodegenerative disorders such as Alzheimer's and Parkinson diseases. Computer simulations of the processes that lead to the formation of these oligomeric species are starting to make significant contributions to our understanding of the determinants of protein aggregation. We simulate different systems at constant concentration but with a different number of peptides and we study the how the finite number of proteins affects the underlying free energy of the system and therefore the relative stability of the species involved in the process. If not taken into account, this finite size effect can undermine the validity of theoretical predictions regarding the relative stability of the species involved and the rates of conversion from one to the other. We discuss the reasons that give rise to this finite size effect form both a probabilistic and energy fluctuations point of view and also how this problem can be dealt by a finite size scaling analysis.

Oligomerization of amyloid Abeta16-22 peptides using hydrogen bonds and hydrophobicity forces.

The 16-22 amino-acid fragment of the beta-amyloid peptide associated with the Alzheimer's disease, Abeta, is capable of forming amyloid fibrils. Here we study the aggregation mechanism of Abeta16-22 peptides by unbiased thermodynamic simulations at the atomic level for systems of one, three, and six Abeta16-22 peptides. We find that the isolated Abeta16-22 peptide is mainly a random coil in the sense that both the alpha-helix and beta-strand contents are low, whereas the three- and six-chain systems form aggregated structures with a high beta-sheet content. Furthermore, in agreement with experiments on Abeta16-22 fibrils, we find that large parallel beta-sheets are unlikely to form. For the six-chain system, the aggregated structures can have many different shapes, but certain particularly stable shapes can be identified.

Two-state folding over a weak free-energy barrier.

We present a Monte Carlo study of a model protein with 54 amino acids that folds directly to its native three-helix-bundle state without forming any well-defined intermediate state. The free-energy barrier separating the native and unfolded states of this protein is found to be weak, even at the folding temperature. Nevertheless, we find that melting curves to a good approximation can be described in terms of a simple two-state system, and that the relaxation behavior is close to single exponential. The motion along individual reaction coordinates is roughly diffusive on timescales beyond the reconfiguration time for a single helix. A simple estimate based on diffusion in a square-well potential predicts the relaxation time within a factor of two.