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.This is the final published version of the paper. It's also available from PLOS at http://www.plosone.org/article/info%3Adoi%2F10.1371%2Fjournal.pone.0106035

Identification of novel modifiers of Aβ toxicity by transcriptomic analysis in the fruitfly.

The strongest risk factor for developing Alzheimer's Disease (AD) is age. Here, we study the relationship between ageing and AD using a systems biology approach that employs a Drosophila (fruitfly) model of AD in which the flies overexpress the human Aβ42 peptide. We identified 712 genes that are differentially expressed between control and Aβ-expressing flies. We further divided these genes according to how they change over the animal's lifetime and discovered that the AD-related gene expression signature is age-independent. We have identified a number of differentially expressed pathways that are likely to play an important role in the disease, including oxidative stress and innate immunity. In particular, we uncovered two new modifiers of the Aβ phenotype, namely Sod3 and PGRP-SC1b

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

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

Intrinsic determinants of neurotoxic aggregate formation by the amyloid β 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 β peptide (Aβ) 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 Aβ40 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 Aβ42 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 Aβ peptide. © 2010 by the Biophysical Society

Role of bulk and of interface contacts in the behaviour of model dimeric proteins

Some dimeric proteins first fold and then dimerize (three--state dimers) while others first dimerize and then fold (two--state dimers). Within the framework of a minimal lattice model, we can distinguish between sequences obeying to one or to the other mechanism on the basis of the partition of the ground state energy between bulk than for interface contacts. The topology of contacts is very different for the bulk than for the interface: while the bulk displays a rich network of interactions, the dimer interface is built up a set of essentially independent contacts. Consequently, the two sets of interactions play very different roles both in the the folding and in the evolutionary history of the protein. Three--state dimers, where a large fraction of the energy is concentrated in few contacts buried in the bulk, and where the relative contact energy of interface contacts is considerably smaller than that associated with bulk contacts, fold according to a hierarchycal pathway controlled by local elementary structures, as also happens in the folding of single--domain monomeric proteins. On the other hand, two--state dimers display a relative contact energy of interface contacts which is larger than the corresponding quantity associated with the bulk. In this case, the assembly of the interface stabilizes the system and lead the two chains to fold. The specific properties of three--state dimers acquired through evolution are expected to be more robust than those of two--state dimers, a fact which has consequences on proteins connected with viral diseases

An effective all-atom potential for proteins

We describe and test an implicit solvent all-atom potential for simulations of protein folding and aggregation. The potential is developed through studies of structural and thermodynamic properties of 17 peptides with diverse secondary structure. Results obtained using the final form of the potential are presented for all these peptides. The same model, with unchanged parameters, is furthermore applied to a heterodimeric coiled-coil system, a mixed alpha/beta protein and a three-helix-bundle protein, with very good results. The computational efficiency of the potential makes it possible to investigate the free-energy landscape of these 49--67-residue systems with high statistical accuracy, using only modest computational resources by today's standards

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

Multi-Scale Simulations Provide Supporting Evidence for the Hypothesis of Intramolecular Protein Translocation in GroEL/GroES Complexes

The biological function of chaperone complexes is to assist the folding of non-native proteins. The widely studied GroEL chaperonin is a double-barreled complex that can trap non-native proteins in one of its two barrels. The ATP-driven binding of a GroES cap then results in a major structural change of the chamber where the substrate is trapped and initiates a refolding attempt. The two barrels operate anti-synchronously. The central region between the two barrels contains a high concentration of disordered protein chains, the role of which was thus far unclear. In this work we report a combination of atomistic and coarse-grained simulations that probe the structure and dynamics of the equatorial region of the GroEL/GroES chaperonin complex. Surprisingly, our simulations show that the equatorial region provides a translocation channel that will block the passage of folded proteins but allows the passage of secondary units with the diameter of an alpha-helix. We compute the free-energy barrier that has to be overcome during translocation and find that it can easily be crossed under the influence of thermal fluctuations. Hence, strongly non-native proteins can be squeezed like toothpaste from one barrel to the next where they will refold. Proteins that are already fairly close to the native state will not translocate but can refold in the chamber where they were trapped. Several experimental results are compatible with this scenario, and in the case of the experiments of Martin and Hartl, intra chaperonin translocation could explain why under physiological crowding conditions the chaperonin does not release the substrate protein