BMC Genomics Protein disorder in the human diseasome: unfoldomics of human genetic diseases

Uros Midic
Center for Information Science and Technology
Temple University
19122PhiladelphiaPAUSA

Christopher J Oldfield email@example.com
Center for Computational Biology and Bioinformatics
Indiana University School of Informatics
46202IndianapolisINUSA

A Keith Dunker
Center for Computational Biology and Bioinformatics
Department of Biochemistry and Molecular Biology
Indiana University School of Medicine
46202IndianapolisINUSA

Zoran Obradovic
Center for Information Science and Technology
Temple University
19122PhiladelphiaPAUSA

Vladimir N Uversky firstname.lastname@example.com
Center for Computational Biology and Bioinformatics
Department of Biochemistry and Molecular Biology
Indiana University School of Medicine
46202IndianapolisINUSA

Institute for Intrinsically Disordered Protein Research
Indiana University School of Medicine
46202IndianapolisINUSA

Institute for Biological Instrumentation
Russian Academy of Sciences
142290Pushchino, Moscow RegionRussia

BMC Genomics Protein disorder in the human diseasome: unfoldomics of human genetic diseases
10.1186/1471-2164-10-S1-S12BioMed Central Page 1 of 24 (page number not for citation purposes) Research * Corresponding authors
Background: Intrinsically disordered proteins lack stable structure under physiological conditions, yet carry out many crucial biological functions, especially functions associated with regulation, recognition, signaling and control. Recently, human genetic diseases and related genes were organized into a bipartite graph(Goh KI, Cusick ME, Valle D, Childs B, Vidal M, et al. (2007)The human disease network. Proc Natl Acad Sci U S A 104: 8685-8690). This diseasome network revealed several significant features such as the common genetic origin of many diseases.Methods and findings:We analyzed the abundance of intrinsic disorder in these diseasome network proteins by means of several prediction algorithms, and we analyzed the functional repertoires of these proteins based on prior studies relating disorder to function. Our analyses revealed that (i) Intrinsic disorder is common in proteins associated with many human genetic diseases; (ii) Different disease classes vary in the IDP contents of their associated proteins; (iii) Molecular recognition features, which are relatively short loosely structured protein regions within mostly disordered sequences and which gain structure upon binding to partners, are common in the diseasome, and their abundance correlates with the intrinsic disorder level; (iv) Some disease classes have a significant fraction of genes affected by alternative splicing, and the alternatively spliced regions in the corresponding proteins are predicted to be highly disordered; and (v) Correlations were found among the various diseasome graph-related properties and intrinsic disorder.Conclusion: These observations provide the basis for the construction of the human-geneticdisease-associated unfoldome.Author summaryMany proteins with important biological functions lack stable structure under physiological conditions. These proteins, being known as intrinsically disordered, are very common in regulation, recognition, signaling and control, and play crucial roles in protein-protein interaction networks. Many of such intrinsically disordered proteins are associated with various human diseases such as cancer, cardiovascular disease, amyloidoses, neurodegenerative diseases, diabetes and others. Recently, human genetic diseases and related genes were organized into a specific network, diseasome. Previous analysis of this diseasome revealed several significant features including the common genetic origin of many diseases. However, the abundance of intrinsically disordered proteins involved in human genetic diseases and the functional repertoire of these proteins have never been before. We filled this gap by performing the thorough bioinformatics analysis of all the proteins form the diseasome utilizing several disorder predictors and by performing the intensive text mining.Here we show that intrinsic disorder is common in diseasome, and that proteins from different diseases possess different levels of intrinsic disorder. Many disordered regions are subjected to alternative splicing and contain specific molecular recognition features responsible for the protein-protein interactions. We also show that many hub proteins are generally more disordered than non-hub proteins. Our study provides the basis for the construction of the human-genetic-disease-associated unfoldome; i.e., a part of the diseasome dealing with the intrinsically disordered proteins.

Introduction

Significant experimental and computational data show that many biologically active proteins lack rigid 3-D structure, remaining unstructured, or incompletely structured, under physiological conditions, and, thus, these proteins exist as dynamic ensembles of interconverting structures. These proteins are known by different names, including intrinsically disordered [1], natively denatured [2], natively unfolded [3], intrinsically unstructured [4], and natively disordered [5] among others. The terms intrinsic disorder (ID), intrinsically disordered protein (IDP), and intrinsically disordered region (IDR) will be used here.

The manifestation of ID is manifold, and functional disordered segments can be as short as only a few amino acid residues or can occupy rather long loop regions and/or protein ends. Proteins, even large ones, can be partially or even wholly disordered. Some IDPs and IDRs exhibit collapsed disordered conformations with pronounced residual structure (thus, resembling a molten globule), others can stay in extended highly disordered states (such as the random coil), while still others form collapsed random coils or semi-collapsed premolten globules [1,[5][6][7][8]. The relationships among the different ID forms needs further study.

There are several crucial differences between amino acid sequences of IDPs/IDRs and structured globular proteins and domains. These differences include divergence in amino acid composition, sequence complexity, hydrophobicity, aromaticity, charge, flexibility index value, and type and rate of amino acid substitutions over evolutionary time. For example, IDPs are significantly depleted in bulky hydrophobic (Ile, Leu, and Val) and aromatic amino acid residues (Trp, Tyr, and Phe), which form and stabilize the hydrophobic cores of folded globular proteins. IDPs also possess a low content of Asn and of the cross-linking Cys residues. The residues that are less abundant in IDPs, and that are more abundant in structured proteins, have been called order-promoting amino acids. On the other hand, IDPs/IDRs are substantially enriched in polar and charged amino acids: Arg, Gln, Ser, Glu, and Lys and in structure-breaking Gly and Pro residues, collectively called disorder-promoting amino acid residues [1,9,10]. Thus, in addition to the well-known "protein folding code" stating that all the information necessary for a given protein to fold is encoded in its amino acid sequence [11], we have proposed that there exists a "protein non-folding code", according to which the propensity of a protein to stay intrinsically disordered is likewise encoded in its amino acid sequence [12,13].

Amino acid differences between IDPs and ordered proteins have been utilized to develop numerous disorder predictors, including PONDR ® (Predictor of Naturally Disordered Regions) [9], charge-hydropathy plots (CHplots) [14] and IUPred [15] to name a few. Intrinsic disorder predictors fall into two general groups. Per-residue predictors (such as the PONDR ® group of predictors) output a score for each residue in a protein and are especially useful when applied to proteins having both structured and disordered regions. The other type of algorithm gives a single prediction value for the entire protein. This type is useful when the objective is to identify mostly or wholly disordered or structured proteins. The charge-hydrophobicity (CH)-plot and the cumulative distribution function (CDF) are the two main predictors of this type [16].

The current state of the art in the field of IDP predictions, including advantages and drawbacks, has been summarized recently [17]. Links to many of the servers for these predictors, when available, can be found in the Disordered Protein Database, DisProt http://www.dis prot.org [18].

Although experimentally characterized IDPs have been discussed in the literature over at least four decades, these proteins have not been viewed as a group but rather as a collection of unusual protein outliers. Bioinformatics is playing a major role in transforming this collection of examples into a sub-field of protein science. For example, soon after the first disorder predictor was developed [19], it was shown that 25% of proteins in Swiss-Prot contained predicted ID regions longer than 40 consecutive residues and that about 11% of residues in Swiss-Prot were likely to be disordered [20]. Subsequent analyses confirmed these trends and revealed that eukaryotic proteomes are significantly more enriched in IDPs in comparison to bacterial and archaeal proteomes [16,21]. This increased utilization of IDPs in higher organisms was attributed to the greater need for signaling and coordination among the various organelles in the more complex eukaryotic domain [1,22].

IDPs carry out numerous biological functions, many of which obviously rely on high flexibility and lack of stable structure. These functions are diverse and complement those of ordered proteins and protein regions. While structured proteins are mainly involved in molecular recognition leading to catalysis or transport, disordered proteins and regions are typically involved in signaling, recognition, regulation, and control by a diversity of mechanisms [23][24][25].

IDPs play crucial roles in protein-protein interaction networks, which generally involve a few proteins binding to many partners (called hub proteins or hubs) and many proteins interacting with just a few partners. Consideration of structure data revealed that several hub proteins are entirely disordered, from one end to the other, and to be capable of binding large numbers of partners, other hubs contain both ordered and disordered regions, and some hubs are structured throughout [26]. Fully disordered hubs can serve as scaffolds for organizing the components of multi-step pathways [27]. For the mixedstructure hubs, many, but not all, of the interactions map to the regions of disorder. For the highly structured hubs (such as 14-3-3 [28] and calmodulin [29]), the binding regions of their partner proteins are intrinsically disordered [30]. Overall, these observations support two previously proposed mechanisms by which ID is utilized in protein-protein interactions: namely, one disordered region binding to many partners and many disordered regions binding to one partner [30,31].

The binding diversity of IDPs plays important roles in the establishment, regulation and control of various signaling networks. Such disorder-based signaling is further modulated in multicellular eukaryotes by posttranslational modification and by alternative splicing, both of which very likely occur much more often in IDRs compared to structured regions of proteins [25,32]. Locating alternative splicing in disordered regions avoids the folding problems that arise upon removal of segments of from structured domains. The flexibility of IDRs facilitates the binding of the enzymes that bring about the disorderassociated posttranslational modifications. We have suggested that the intersection of binding sites, posttranslational modifications, and alternative splicing variants within IDRs provide a powerful combination to bring about signaling diversity in different cell types [25,28,32] Many IDPs and IDRs fold upon binding with their specific partners. Said partners include other proteins, nucleic acids, membranes or small molecules [33]. The concept of the "molecular recognition feature," abbreviated as MoRF, was introduced to describe short, intrinsically disordered regions that "morph" from disorder-to-order upon partner recognition [34][35][36]. Based on several specific features in the disorder prediction scores, a predictor of helixforming MoRFs was elaborated [34,37]. The application of this predictor to several proteomes revealed that such foldable recognition features are especially abundant among eukaryotic proteins [34,37]. MoRFs that form sheet or irregular structure also exist [35,36]. Predictors of these non-helical MoRFs have not yet been developed, so the predictions of helix-forming MoRFs should be regarded as providing lower-bound estimates of binding sites in disordered regions.

Proteins are involved in virtually all cellular and in many extracellular processes. Protein dysfunction can therefore cause development of various pathological conditions and a broad range of human diseases known as proteinconformation or protein-misfolding diseases. Such diseases arise from the failure of a specific peptide or protein to adopt its functional conformational state; i.e., from protein misfolding and malfunctioning.

Misfolding diseases can affect a single organ or be spread through multiple tissues. Consequences of misfolding include protein aggregation, loss of normal function, and gain of toxic function. Misfolding and misfunction can originate from point mutation(s) or result from an exposure to internal or external toxins, from impaired posttranslational modification (phosphorylation, advanced glycation, deamidation, racemization, etc.), from an increased probability of degradation, from impaired trafficking, from lost binding partners or from oxidative damage among other causes. These factors can act independently or in complex associations with one another [38]. Furthermore, numerous IDPs are associated with human diseases such as cancer [22], cardiovascular disease [39], amyloidoses [40], neurodegenerative diseases [41], diabetes and others [38]. Based on these intriguing links among intrinsic disorder, cell signaling and human diseases, suggesting that protein conformational diseases may result not only from protein misfold-ing, but also from misidentification and missignaling [30], the "disorder in disorders" or D 2 concept was recently introduced [38].

As a result of decades-long efforts, impressive lists of disease-gene association pairs were generated [42,43]. In parallel, analysis of protein-protein interactions in humans produced detailed maps of the relationships between different genes including those related to disease [44,45]. To gain a better understanding of the relationship between the genes implicated in a selected disease, network-based tools were successfully utilized for a single disease, e.g., human inherited ataxias and disorders of Purkinje cell degeneration [46].

Recently, to estimate whether human genetic diseases and the corresponding disease genes are related to each other at a higher level of cellular and organism organization, a bipartite graph was utilized in a dual way: to represent a network of genetic diseases, the "human disease network", HDN, where two diseases are directly linked if there is a gene that is directly related to both of them, and a network of disease genes, the "disease gene network", DGN, where two genes are directly linked if there is a disease to which they are both directly related [47]. This framework, called the human diseasome, systematically linked the human disease phenome (which includes all the human genetic diseases) with the human disease genome (which contains all the disease-related genes). This diseaseome opens a new avenue for the analysis and understanding of human genetic diseases, moving from single gene-single disease viewpoint to a frameworkbased approach [47].

The analysis of the HDN and DGN properties revealed that these networks are significantly different in many aspects from randomly generated networks of the same size. By these analyses the various diseases became classified into 20 types, some diseases were unclassified, and several diseases were annotated as belonging to multiple classes. Similarly, genes were clustered into classes via their associations with specific diseases [47]. Analysis of this network of genetic diseases and disease genes linked by known disease-gene associations revealed the common genetic origin of many diseases. The vast majority of these disease genes was non-essential and showed no tendency to encode hub proteins. Overall, the expression pattern of these disease-related genes indicated that they are localized in the functional periphery of the network [47].

In the present study, we started from the disease-related classification of genes from [47] and then performed a large-scale analysis of the abundance of intrinsic disorder in transcripts of the various disease-related genes. Since structural information was available only for a limited number of these proteins, we used intrinsic disorder predictions. We also analyzed the correlation between various HDN/DGN graph-related properties of genes and intrinsic disorder. We compared the occurrence of alternative splicing in various disease classes and analyzed the relationship between alternative splicing and intrinsic disorder. In essence, the aim of our study was to build an unfoldome, which we define as the IDP-containing subset of a given genome, associated with human genetic diseases.

Overall, our findings indicate that there are significant differences in occurrence of intrinsic disorder in the proteins arising from genes related to diseases as compared to proteins arising from genes unrelated to specific diseases. Furthermore, there are significant differences with respect to intrinsic disorder among the various disease classes. Our analysis shows noticeable positive trends that link intrinsic disorder to graph-related features of genes, such as the number of other genes that are directly linked to a given gene via the diseasome network. Certain disease classes have a significantly greater fraction of genes involved in alternative splicing, and these alternative splicing regions are predicted to be highly disordered. In summary, disorder analysis provides interesting new insights regarding the human diseasome.

Methods

The basis for our experimental dataset is the dual Human Disease Network/Disease Gene Network (HDN/DGN) [47]. It consists of two types of nodes that represent human genes (1,777) and diseases (1,284), and links that connect diseases with related genes. A disease and a gene were connected by a link if mutation(s) in the corresponding gene were implicated in the given disease [47]. The network is dual, because it can be observed as both a Human Disease Network (two diseases are linked if they are both related to the same gene), or as a Disease Gene Network (two disease genes are linked if they are both related to the same disease).

We augmented the set of disease genes from DGN with human genes with known protein sequences. Protein sequences for all human genes were collected from NCBI Gene database; we excluded all model proteins obtained solely with automated genome annotation processing. After this exclusion, our dataset consists of 1,751 human disease related genes and 16,358 other human genes with known protein sequences. If several protein sequences were collected for a single gene; i.e., for genes with multiple alternatively spliced isoforms, then any duplicate sequences were discarded.

The diseases in DGN were grouped into twenty classes. In addition to these twenty classes we introduced sets of unclassified diseases and diseases belonging to multiple classes as two separate disease classes. We used this approach to classify genes as well. In our model, a gene belongs to all classes to which its related diseases also belong. Furthermore, since a gene can be related to multiple diseases that belong to various classes, we defined an additional multiple class gene group. Thus, overall, this approach defined 22 gene classes: the twenty original classes, as well as classes of unclassified genes (related to unclassified diseases) and multi-class disease genes (genes related to diseases that belong to multiple classes). Note that the 22 gene classes were not necessarily disjoint, and that all genes from multiple class gene class also belonged to at least two more classes. Two more sets were used for comparison: disease genes (this set included all genes from DGN with known protein sequences; i.e., genes from all 22 previously defined classes), and human genes (this was the whole dataset that included the disease genes set). Table 1 contains preliminary statistics for 22 disease/gene classes and 3 additional classes of genes, namely multiple class genes, disease genes, and human genes.

Intrinsic disorder prediction

Three predictors of intrinsic disorder were used on the protein sequences: PONDR ® VSL2B, CH and CDF. VSL2B is a variant of VSL2 predictor described in [48]. For an amino acid sequence, VSL2B outputs ID prediction in [0, 1] range per residue. These outputs were then compared to a threshold (we used the default threshold 0.5) and residues with prediction value greater than the threshold were predicted to be ID. In the case of multiple sequences for one gene, sequences were aligned using our own multiple alignment algorithm, which was aimed at rediscovering identical exons in multiple sequences by only matching identical amino acids and optimizing the alignment for long contiguous matched subsequences. A sequence obtained from such multiple alignments included all exons from individual sequences, and was considered to represent the whole gene sequence. For each position in the alignment sequence, we obtained a single prediction by averaging predictions for all residues from protein sequences that are aligned at that position.

CH and CDF give outputs that predict disorder on the level of whole proteins. The CH (Charge-Hydrophobicity) The same order of classes is used in graphs in the Results section; the first 22 classes are sorted in descending order with respect to the median of disorder content (defined in Experimental procedures). The difference between "multi-class diseases" and "multiple class genes" is that "multi-class diseases" set includes genes that are only related to diseases that are classified as "multiple" in [47], whereas "multiple class genes" includes genes that are related to several diseases that belong to different classes.

predictor is based on the finding [14] that two sets of proteins -a set of natively unfolded proteins and a set of small globular folded proteins -occupy two distinct regions in the charge-hydrophobicity phase space, and can be almost perfectly separated with a straight line. The CH predictor calculates the mean hydrophobicity and the mean net charge for a protein sequence, identifies the part of the charge-hydrophobicity plane that the corresponding point belongs to, and calculates its distance from the separating line. The CDF predictor [16,49] compiles the predictions of a per-residue predictor to a single binary predictor per protein, by observing the cumulative distribution function (CDF) of per-residue predictions, and comparing it to a set of 7 boundary CDF points obtained from a training set [16]. In the case of multiple sequences for one gene, we used weighted voting to determine a single prediction for the gene. For the CH predictor, we calculate the mean of signed distances (distance is multiplied by -1 if prediction is negative, i.e. protein is predicted to be ordered). The prediction for the gene depends on the sign of the weighted mean (disorder if the weighted mean is positive, order otherwise). Similarly to the CH predictor, CDF predictor has a parameter (CDF count), the mean of which over all proteins sequences for a gene is compared to the threshold to determine a single prediction for the gene.

Since VSL2B provides per-residue predictions, we measure the disorder content, which is defined as the fraction of residues in a protein sequence (or sequence alignment in case of alternative splicing) that is predicted to be disordered. This provides a single prediction value for a given gene. Note that, unlike the CH and CDF predictors, this prediction can take any value in the range [0, 1].

-MoRF predictions

The predictor of an -helix forming Molecular Recognition Feature (-MoRF) is based on observations that predictions of order in otherwise highly disordered proteins corresponds to protein regions that mediate interaction with other proteins or DNA. This predictor focuses on short binding regions within long regions of disorder that are likely to form helical structure upon binding [34]. It uses a stacked architecture, where PONDR ® VLXT is used to identify short predictions of order within long predictions of disorder and then a second level predictor determines whether the order prediction is likely to be a binding site based on attributes of both the predicted ordered region and the predicted surrounding disordered region. An -MoRF prediction indicates the presence of a relatively short (~20 residues), loosely structured helical region within a largely disordered sequence [34]. Such regions gain functionality upon a disorder-to-helix transition induced by binding to partner sequences [35,36].

Recently it has been indicated that the -MoRF predictor has a poor sensitivity, i.e., misses many -MoRF regions [37], due to the small set of -MoRF regions used in its development. In this study, the modified -MoRF predictor, -MoRF-PredII, was used [37]. This algorithm was improved by including additional -MoRF examples and their cross species homologues in the positive training set, carefully extracting monomer structure chains from PDB as the negative training set and including attributes from recently developed disorder predictors, secondary structure predictions, and amino acid indices as attributes [37].

Alternative splicing analysis

For genes with multiple isoforms, the multiple alignments provide the information on the alternative splicing regions. We define the alternative splicing regions (AS regions) as exons that are expressed in some, but not all protein sequences for a given gene. Similarly as for a whole gene, we define disorder content for an AS region as the fraction of its residues that are predicted to be disordered.

Statistical analysis of the data

When disorder content measurements -as predicted by VSL2B predictor -for all genes in a disease class were observed as a sample, we used statistical tests to compare the samples arising from different disease classes. Since we cannot make any assumptions on the distributions for disorder content in disease classes, we used the nonparametric Mann-Whitney U test (Wilcoxon rank-sum test) [50,51] to test whether two samples of observations (i.e. disorder content for two classes) came from the same distribution. The Mann-Whitney U test was not appropriate for similar comparison in the case of CH and CDF predictors, as their predictions were binary. For these two predictors, we counted the number of positive (disordered) and negative (ordered) observation in two samples (classes) and then used the  2 test to estimate the likelihood of whether the two samples come from the same distribution.

We dealt with the possible problems of multiple hypotheses testing by controlling false discovery rate (FDR) with the Benjamini-Hochberg (for independent tests) [52] or with the Benjamini-Yekutieli method [53].

Several of our hypotheses dealt with the dependency between graph-related numeric properties of nodes representing genes and disorder content. The numeric properties were defined as:

• number of related diseases: number of diseases the gene is directly related to (as provided in [47])

• number of related disease classes: number of distinct disease classes that diseases related to the gene belong to

• degree: number of other genes that are related to the diseases the gene is related to; or defined in the terms of DGN graph: the number of other genes that are directly linked to the gene (through some disease node).

For such hypotheses we used (first-order) linear regression to model the relationship, and then we used the corresponding F-statistic to assess the validity of the linear model.

The HDN/DGN graph was not completely connected.

Using the usual definition of connectivity in graphs, we identified the connected components. One of the components was large and included 516 disease nodes and 903 gene nodes. All of the remaining components contained 15 genes or less; for example, 399 components contain only one gene each. We split the set of disease genes (DIS) into the set of 896 disease genes that belong to the large component (LARGECOMP) and the set of 855 disease genes that belong to one of the smaller components (SMALLCOMPS). Note that although the 16 disease genes with no available protein sequences were not included in the DIS set, and therefore neither in LARGECOMP nor the SMALLCOMPS set, these 16 genes were still included in the HDN/DGN graph for the purpose of identification of connected components.

Results

Analysis of ID in human diseasome

Prediction of intrinsic disorder using PONDR ® VSL2B predictor on all 30053 initially collected protein sequences showed significant differences in predicted ID content for the 7525 (25.04%) model protein sequences obtained with automated genome annotation processing, and the 22528 (74.96%) protein sequences with additional experimental support. The medians of disorder content for model protein sequences was much higher (68.6% vs. 37.5%), as well as the first quartile (37.9% vs. 21.4%) and the third quartile (96.5% vs. 61.7%). Furthermore, 40.6% of model protein sequences were predicted to have disorder content above 80%, compared to only 11.3% for remaining sequences.

The boxplot in Figure 1 depicts the distributions of disorder content for genes in 25 classes. The 22 disease classes are sorted according to their medians of disorder content. The distributions for the majority of classes appear to be positively skewed. The ranges of disorder content between the first and the third quartile differ greatly between classes. For example, connective tissue disorder (CTD) class is ranked eleventh in disorder content median among the 22 disease classes, but has the highest third quartile.

The distributions of disorder content in disease classes are further compared in histograms in Figure 2. The various classes have irregular disorder content distributions that can hardly be fit by any of the standard distributions. Furthermore, the distributions associated with the different disease classes are dissimilar both in shape and size. For these reasons we use a nonparametric test, Wilcoxon ranksum test [50,51], to compare the distributions by comparing their medians. Figure 3 shows an overview of pair-wise comparisons of disorder content medians. We used Benjamini-Yekutieli (BY) method of false discovery rate (FDR) control [53], as the family-wise error rate multiple comparisons methods, such as the Tukey-Kramer method [54,55], are much more conservative. With an FDR rate of 0.05, it is expected that 2.8 of 56 class pairs reported to have significantly different disorder content medians were false discoveries. The BY method is still quite conservative as it does not make any assumption on the independence of the pair-wise com- parisons. Therefore we included Table 2 [52], which makes an assumption that the tests are independent.

Comparison of disorder content distributions in disease classes and human gene class using boxplots

We continued with the investigation of the relationship between disorder content and several HDN/DGN graphrelated properties. We used linear regression to model disorder content as a linear function of number of related diseases for a gene (Figure 4), number of related disease classes for a gene ( Figure 5), and gene degree in DGN ( Figure 6).

For all three cases, the F-test gave p-values that were smaller than 0.05; for the number of related diseases and gene degree the p-values were smaller than 0.01. Although it is not likely that the observed linear trends were obtained by pure chance, they explained only a very small amount of variation in the disorder content; the respective R 2 values were 6.12·10 -3 , 3.51·10 -3 , and 6.10·10 -3 .

The disease genes set DIS is split almost evenly between LARGECOMP, the 896 (51.17%) disease genes in the large DGN component, and SMALLCOMPS, the 855 (48.83%) disease genes in the remaining small DGN components. This split can be further observed in individual disease classes. The histogram in Figure 7 shows the split between LARGECOMP and SMALLCOMPS for all disease gene classes. Using the  2 test to compare the split in each class to the overall split in the disease gene set, we identified classes of disease genes that were significantly overrepresented or underrepresented in LARGECOMP. For example, 85.99% of genes related to cancer diseases belonged to the large component, while only 19.03% of genes related to metabolic diseases belonged to the large Pairwise comparison of disorder content medians for disease classes and human gene class Comparison of disorder content distributions in disease classes and human gene class using stacked histograms Figure 2 Comparison of disorder content distributions in disease classes and human gene class using stacked histograms. The histograms are stacked horizontally to save space. They show what fraction of genes in each class has disorder content within various ranges. Each of the five major ranges, that cover 20% each, is further split into two smaller 10% ranges (they use the same color, but are divided with a line). Distributions can be visually compared by observing the balance between darker and lighter shades of gray; the class with a darker histogram has on average more disorder content.

component. We then compared the medians of disorder content for genes from LARGECOMP and SMALLCOMPS for each class individually, as well as for the whole disease genes set. The median of disorder content for LARGE-COMP genes was significantly greater than for SMALL-COMPS genes, with an adjusted p-value of 7.56·10 -7 on the rank sum test. Similarly, the median of disorder content for LARGECOMP genes related to metabolic diseases was significantly greater than for the SMALLCOMPS genes related to metabolic diseases, with an adjusted p-value of 0.0112. These comparisons are illustrated in Figure 8. Substantial differences between disorder content medians for genes in LARGECOMP and genes in SMALLCOMPS can also be observed for several other classes; in the majority of cases, the median for the LARGECOMP genes is greater than the median for the SMALLCOMPS genes. However, none of these differences were statistically significant; which was partially due to the small numbers of genes in subsets compared.

Alternative splicing and ID in human diseasome

We applied similar methodology to analyze alternative splicing. We divided the set of all genes (HUM) into the set of genes with multiple isoforms and the set of genes with a single isoform. The same division can also be applied to all disease classes, and the disease gene set. The comparison of fractions of genes with multiple isoforms is shown in Figure 9.

The disease gene set DIS had significantly higher fraction of genes with multiple isoforms than the human gene set HUM. Out of 1751 disease related genes, 410 genes (23.4%) had multiple isoforms (average of 2.77 for disease related genes with multiple isoforms), and they included 991 alternatively spliced regions (2.41 AS regions per disease-related gene with multiple isoforms).

Linear regression of disorder content with respect to number of related disease classes (for genes) Figure 5 Linear regression of disorder content with respect to number of related disease classes (for genes). The genes with number of related disease classes between 1 and 3 are represented as a boxplot, while the remaining genes are represented as points. Note that the disorder content means (inverted triangles) for subsets are greater than the respective medians, because the disorder content distributions in these subsets are positively skewed. Linear regression of disorder content with respect to number of related diseases (for genes) Figure 4 Linear regression of disorder content with respect to number of related diseases (for genes). The genes with number of related diseases up to 4 are represented as a boxplot, while the remaining genes are represented as points.

Note that the disorder content means (inverted triangles) for subsets are greater than the respective medians, because the disorder content distributions in these subsets are positively skewed.

Out of 16358 non-disease genes, 2445 (14.95%) had multiple isoforms (average of 2.51 for non-disease genes with multiple isoforms), and they included 4954 AS regions (2.02 AS regions per non-disease gene with multiple isoforms).

Furthermore, all the disease classes but one (unclassified diseases) had higher fraction of genes with multiple isoforms than the HUM set, and for several classes this difference in fractions was statistically significant. The highest fraction of genes with multiple isoforms was 40.10% for the cancer disease gene class.

The comparisons of distributions of disorder content for genes with multiple isoforms with genes with single isoform showed that for three sets the medians of disorder content for genes with multiple isoforms were significantly greater than for genes with single isoform: human genes set HUM (BY adjusted p = 1.50·10 -7 ), disease genes set DIS (BY adjusted p = 5.08·10 -7 ) and multiple class genes set MULT (adjusted p = 0.0176). Individual tests for three disease classes also returned low p-values (hematological, p = 0.0196; renal, p = 0.0283; bone, p = 0.0291), but the corresponding BY adjusted p-values were above  = 0.05. Figure 10 shows the distributions of disorder content for genes with multiple isoforms (disease, non-disease, and all genes) and for all human genes. Although there are significant differences in medians, the distributions have similar shapes; the peaks are in the 20-40% range, and the fractions decrease with the increase in the disorder content. Figure 10 also shows the disorder content distributions for AS regions in disease related and non-disease genes. These two distributions have different shape than shape of the disorder content distributions for whole pro-teins; the fractions decrease with the increase in the disorder content, but then suddenly increase in the 80-100% range. The rank sum test for medians shows that the distribution of disorder content in AS regions is significantly different from distributions of disorder content in whole proteins for all genes (p ~10 -142 ), as well as for subset of genes with multiple isoforms (p ~10 -48 ). However, as is clearly seen in Figure 10, the distributions of disorder content for AS regions in disease genes and non-disease genes were not significantly different (p = 0.5278). We compared the disorder content distributions for AS regions for genes from individual classes to the overall distribution for AS regions from all human genes. The distributions for classes with significant statistical results are shown in Figure 11. For developmental and neurological disease classes, the fraction of AS regions in 80-100% range is significantly increased. Similarly, there is an increase in 0-20% range for hematological disease class. Metabolic disease class is an extreme case, as there is both a big increase in 0-20% range and decrease in 80-100% range; the AS regions in metabolic disease genes have significantly less disorder when compared to whole sequences in human genes. Linear regression of disorder content with respect to gene degree in DGN Figure 6 Linear regression of disorder content with respect to gene degree in DGN.

-MoRFs in the human diseasome

genes with predicted -MoRFs are highly correlated to fractions of disordered residues (corr. coefficient ~0.89).  Table 3 lists the quotients of the MoRF density in AS regions over the overall MoRF density, as well as corresponding p-value for comparison of these densities.

-MoRFs and alternative splicing in the human diseasome

Evaluation of ID by binary classifiers

We compared the fractions of genes predicted to be disordered by per-protein predictors CDF and CH in Figure 14.

Overall, the CDF predictor identified more genes to be disordered than the CH predictor. The relationship between alternative splicing and intrinsic disorder, as predicted by CDF and CH predictors, can only be observed at the level of whole proteins. For the CDF predictor, classes with significantly different fractions of predicted disordered genes in genes with single isoform and in genes with multiple isoforms were: DIS, HUM, RENA, MULT, and (borderline significance) HEMA; in all cases, fraction of predicted disordered genes for genes with multiple isoforms was greater than for genes with a single isoform. For the CH predictor, significant difference of fractions of predicted disordered genes in genes with a Comparison of fractions of disease genes with multiple iso-forms (i.e. with alternative splicing) and with single isoform Figure 9 Comparison of fractions of disease genes with multiple isoforms (i.e. with alternative splicing) and with single isoform. The classes with the + signs after their acronyms have significantly higher fraction of genes with multiple isoforms than the human gene set. The error bars represent one standard deviation or 68.2% confidence interval.

Comparison of distributions of disorder content in LARGE-COMP and SMALLCOMP for genes related to metabolic dis-eases, and for the whole disease gene set Figure 8 Comparison of distributions of disorder content in LARGECOMP and SMALLCOMP for genes related to metabolic diseases, and for the whole disease gene set. Distribution of disorder content for human gene set is included for comparison. The error bars represent one standard deviation or 68.2% confidence interval.

single isoform and in genes with multiple isoforms was only observed in the HUM set.

In general, CDF predicts a much higher fraction of genes to be disordered than CH. Vectors of fractions of predicted disordered genes in various classes for CDF and CH predictors are fairly correlated, though there are several classes with substantial differences. For example, IMMU, RESP, RENA, and UNCL have very low (or even zero) fractions of disordered genes for CH predictor. The relative difference between HUM and DIS sets is much larger for the CH predictor (approximately two-fold) than for the CDF predictor. For the CDF predictor, the fractions of predicted disordered genes for several classes are higher than (although not strictly significantly higher) or similar to the same fraction in the HUM set, while this is not the case for CH predictor.

Overall, the fractions of predicted disordered genes for both binary predictors are correlated to medians of disorder content for VSL2B predictors, but there are some striking differences. For example, low fractions of disordered genes in IMMU class for both predictors, or relatively high fraction of PSYC and RESP classes for CDF predictor. Looking only at the medians without at least comparing whole distributions is not a good way to compare prevalence of intrinsic disorder in two classes/sets of genes.

The difference between these two methods in the magnitude of predicted disorder is generally similar to previously published data [16,39,56]. This difference was explaine by the fact that the CH-plot is a linear classifier that takes into account only two parameters of the particular sequence -charge and hydropathy, whereas CDF analysis is dependent upon the output of the PONDR ® VLXT predictor, a nonlinear neural network classifier, which was trained to distinguish order and disorder based on a significantly larger feature space that explicitly includes net charge and hydropathy. According to these methodological differences, CH-plot analysis is predisposed to discriminate proteins with substantial amounts of extended disorder (random coils and pre-molten globules) from proteins with globular conformations (molten globule-like and rigid well-structured proteins). On the other hand, PONDR-based CDF analysis may discriminate all disordered conformations including molten globules from rigid well-folded proteins. Therefore, this discrepancy in the disorder prediction by CDF and CHplot might provide a computational tool to discriminate proteins with extended disorder from native molten globules, which might be predicted to be disordered by CDF, but compact by CH-plot. This model is consistent with the behavior of several intrinsically disordered proteins (e.g., [57]). Comparison of disorder content distributions for whole pro-teins and for alternative splicing (AS) regions Figure 10 Comparison of disorder content distributions for whole proteins and for alternative splicing (AS) regions. Series #1-4 represent disorder content distributions in 1) disease genes with multiple isoforms, 2) non-disease genes with multiple isoforms, 3) human genes with multiple isoforms, 4) all human genes. Series #5 and #6 represent disorder content distributions for AS regions in disease genes, and AS regions in non-disease genes. The error bars represent one standard deviation or 68.2% confidence interval.

disease within the CH-CDF phase space. In these plots, each spot corresponds to a single protein and its coordinates are calculated as a distance of this protein from the boundary in the corresponding CH-plot (Y-coordinate) and an averaged distance of the corresponding CDF curve from the boundary (X-coordinate). Positive and negative Y values correspond to proteins which, according to CHplot analysis, are predicted to be natively unfolded or compact, respectively. Whereas positive and negative X values are attributed to proteins that, by the CDF analysis, are predicted to be ordered or intrinsically disordered, respectively. Therefore, each plot contains four quadrants: (-, -) contains proteins predicted to be disordered by CDF, but compact by CH-plot (i.e., potential native molten globules); (-, +) includes proteins predicted to be disordered by both methods (i.e., proteins with extended disorder); (+, -) contains ordered proteins; (+, +) includes proteins predicted to be disordered by CH-plot, but ordered by the CDF analysis. A sharp cut-off at the right side of each plot is due to the upper limit of a difference between the CDF curve (which might have a maximum value of 1.0) and a boundary separating IDPs and ordered proteins in CDF plots. Figure 15 suggests that the majority of the wholly disordered proteins could possibly be native molten globules.

Discussion

Intrinsic disorder in predicted and experimentally identified proteins

The difference in predicted disorder content between sets of predicted model protein sequences and (partially) confirmed protein sequences is greater in both magnitude and statistical significance than difference between any two classes in our final data set. The simplest explanation for this is that automated annotation procedure has a high error rate that introduces a large number of incorrect amino acid sequences. Alternatively, this dramatic difference in the level of predicted ID between the experimentally and automatically identified proteins could be due the bias of the existing identification techniques toward the ordered proteins. To some extent this resembles a problem the Structural Genomics Initiative Centers are facing, where the use of the traditional target search criteria (mostly based on the sequence identity) and protein purification and isolation methods generated mostly ordered targets, whereas alternatively identified and purified proteins awaiting structure determination were richer in disorder than an average protein in PDB [58,59]. It has been pointed out that this bottleneck was determined by the strategy chosen were in efforts to identify proteins with novel folds researchers started with proteins having amino acid sequences unlike those of proteins with known 3D structures [58,59]. In a similar manner, tradi- tional experimental approaches developed for protein identification could be biased toward order (as ordered well-folded proteins where at the research focus for many years), whereas predictive tools are mostly dealing with the remaining part of the proteomes and therefore are inevitably identifying more disordered proteins.

Comparison of overall density of predicted MoRFs vs density of predicted MoRFs in AS regions for 25 classes/sets

The predicted sequences were unevenly distributed between disease-related and disease-unrelated proteins. In fact, the majority of the predicted "model" sequences were products of the non-disease genes. Therefore, including such sequences into the data set would introduce significant bias for disorder in the non-disease gene part of the data set. Based on these observations, we decided to exclude such sequences from the final datasets.

An important assumption that we made was that ID predictors have no bias towards any class of genes. Although errors are unavoidable in prediction of disorder, we assumed that both false positive and false negative errors occur equally likely in all gene classes. Under this assumption we can expect that any observed variations in predicted disorder content between disease classes are due to real variations in disorder content and not due to bias introduced by prediction. Although we have not found any obvious reason for questioning these assumptions, more structural data are needed to test such biases.

Fractions of genes predicted to be disordered by CDF and CH predictors Figure 14 Fractions of genes predicted to be disordered by CDF and CH predictors. The error bars represent one standard deviation or 68.2% confidence interval.

Intrinsic disorder in human genetic diseases

Contrary to our initial expectations based on known abundance of ID in such diseases as cancer [22], cardiovascular disease [39], amyloidoses [40], neurodegenerative diseases [41], diabetes and others [38], the disease genes have in general slightly lower disorder content than the non-disease genes. This can be explained by the fact that the human disease network (HDN) and disease gene network (DGN) are based on the genetic diseases and genes, mutations in which were associated with disease development, respectively. Based on the expression pattern analyses of the DGN genes it has been concluded that they are mostly localized in the functional periphery of the protein-protein interaction network [47]. This peripheral localization of most disease genes was explained assuming that mutations in topologically central, highly connected, and widely expressed genes were more likely to result in severe impairment of normal development, leading to early lethality and therefore to deletion from the population, whereas mutations compatible with survival into the reproductive years were more likely to be maintained in a population [47]. Overall, the vast majority of disease genes in DGN was non-essential and showed no tendency to encode hub proteins [47]. On the other hand, the above mentioned studies on various individual diseases [22,[38][39][40][41] dealt with all proteins known to be associated with a given disease and not just those proteins bearing the disease-promoting mutations. Therefore, the various datasets of proteins associated with individual diseases contained wider variety of proteins, including hubs. It is important to remember that hub proteins were shown to be highly enriched in intrinsic disorder [26][27][28][60][61][62][63][64]. In fact, hubs were shown to have multiple interactions, either being intrinsically disordered and serving as an anchor, or acting as a stable globular scaffold that interacts with intrinsically disordered regions of its targets [26][27][28][60][61][62][63][64]. Therefore a systematic depletion of hub proteins in HDN and DGN can in part explain their slightly lower disorder contents.

Functional analysis of the disorder predictions for several specific diseases

Our data revealed that there was a large variation in disorder content in various disease classes. Several disease classes had median disorder content higher or comparable with the human gene set HUM (SKEL, BONE, DERM, Figure 15 Comparison of CDF and CH predictions in various disease gene classes and gene sets. Each spot represents a gene whose coordinates were calculated as the distance of the corresponding point in the CH-plot from the boundary (x-coordinate) and the averaged distance of the corresponding CDF-curve from the CDF boundary (y-coordinate). Four quadrants in each plot correspond to the following predictions: (-,-) proteins predicted to be disordered by CDF, but compact by CH, (-,+) proteins predicted to be disordered by both methods, (+,-) ordered proteins, (+,+) proteins predicted to be disordered by CH, but ordered by CDF. This is further illustrated by an explanatory plot at the bottom right corner. Percentages represent the fractions of genes in the corresponding quadrants.

Comparison of CDF and CH predictions in various disease gene classes and gene sets

CANC, DEVE, MCD), and/or several classes had a higher or comparable fraction of highly disordered genes (SKEL, BONE, CTD, DERM, DEVE, CTD, PSYC, MCD).

Unfortunately, structural information on proteins associated with various genetic diseases is sparse. Therefore, in the analyses given below we used an established earlier correlation between various protein functions and ID [23][24][25]. The strongest correlation with ID was observed for regulatory functions involved in biological processes such as those described by the following functional keywords: differentiation, transcription, transcription regulation, spermatogenesis, DNA condensation, cell cycle, mRNA processing, mRNA splicing, mitosis, apoptosis, protein transport, meiosis, cell division, Ubl conjugation pathway, Wnt signaling pathway, chromosome partition, neurogenesis, ribosome biogenesis, chondrogenesis, growth regulation [23][24][25]. The major ID-associated functional keywords covered a wide spectrum of protein activities including ribonucleoprotein, ribosomal protein, develop-mental protein, chromatin regulator, hormone, growth factor, GTPase activation, cytokine, GAP protein, repressor, cyclin, activator protein phosphatase inhibitor.

Many of these processes and functions are performed by proteins associated with various human genetic diseases. Based on this correlation between ID and protein function, a given protein was assigned either to the ordered protein family or to the IDR class, assuming that if this protein possesses ID-associated function then it likely contains at least one long IDR. A brief overview of IDenriched classes of human genetic diseases is presented below. In addition, several examples of disease-related genes encoding for important well-characterized IDPs are discussed. This includes disordered proteins, mutations in which were associated with particular diseases, -synuclein (one of the major players in the Parkinson's diseases pathogenesis), p53 (a key tumor-supressor protein), huntingtin (a protein involved in the Huntington's disease pathogenesis), BRCA1 (a breast and/or ovarian can-  The class identifier is followed by the number of genes with predicted -MoRFs, the total number of genes, and the fraction of genes with predicted -MoRFs in that class. This is followed by the list of genes with multiple (>= 2) predicted -MoRFs. The numbers following each gene symbol are the number of predicted -MoRFs, the number of amino acids and the disorder content (note that in the case of alternative splicing the number of amino acids includes all exons in a gene and may be larger than the lengths of individual isoforms).

cer-associated protein) and EWS-FLI1 fusion protein (a protein associated with the Ewing's sarcoma family of tumors).

Illustrative examples of well-studied IDPs encoded by disease-related genes

-Synuclein is a typical IDP that links various synucleinopathies, a group of neurodegenerative disorders characterized by deposition of aggregated -synuclein in the cytoplasm of selective populations of neurons and glia [65][66][67][68]. Clinically, synucleinopathies are characterized by a chronic and progressive decline in motor, cognitive, behavioral, and autonomic functions, depending on the distribution of the lesions. Some of the most common synucleinopathies include Parkinson's disease, dementia with Lewy bodies, Alzheimer's disease, Down's syndrome, multiple system atrophy, and neurodegeneration with brain iron accumulation type 1. Different diseases are characterized by the morphologically different -synuclein-containing inclusions. In Parkinson's diseases and various Lewy body diseases these inclusions are Lewy bodies and Lewy; multiple system atrophy is characterized by the accumulation of glial cytoplasmic inclusions and neuronal cytoplasmic inclusions; whereas axonal spheroids are frequently found in neurodegeneration with brain iron accumulation type 1 [69,70].

Several observations implicate -synuclein in the pathogenesis the pathogenesis of Parkinson's disease. For example, a direct role for -synuclein in the neurodegenerative processes in PD and Lewy body dementia was demonstrated by genetic evidence. Autosomal dominant earlyonset Parkinson's disease and Lewy body dementia was shown to be induced in a small number of kindreds as a result of three different missense mutations in the -synuclein gene, corresponding to A30P, E46K, and A53T substitutions in -synuclein [71][72][73] or as a result of the hyper-expression of the wild type -synuclein protein due to gene triplication [74][75][76][77]. Besides this genetic evidence many other observations correlate -synuclein and PD pathogenesis (reviewed in [67][68][69][70][78][79][80][81][82]. Some of these observations are briefly outlined below. The recombinant -synuclein easily assembles into amyloid-like fibrils in vitro and this process is modulated by familial point mutations. Characteristic depositions in various synucleinopathies invariably contain aggregated -synuclein. -Synuclein is abnormally phosphorylated, ubiquitinated, and nitrated in pathology-related inclusions. Co-expression of chaperones or -synuclein with -synuclein in transgenic animals was shown to suppress the neurodegeneration. -Synuclein-positive proteinaceous deposits were shown to accumulate in several animal models where Parkinsonism was induced by exposure to different neurotoxicants. All this indicates that -synuclein is a key player in the pathogenesis of several neurodegenerative disorders Conformational behavior of -synuclein under a variety of environments has been extensively analyzed (for recent reviews see [69,70,82]). This analysis has revealed that the structure of -synuclein is extremely sensitive to the environment and can be easily modified. As a result, -synuclein was shown to possess a remarkable conformational plasticity, being able to adopt structurally unrelated conformations including the substantially unfolded state, an amyloidogenic partially folded conformation, different helical or -structural species folded to a different degree, both monomeric and oligomeric, several morphologically different types of aggregates, including various oligomers, amorphous aggregates, and amyloid-like fibrils [69,70,82]. Based on this astonishing conformational behavior the concept of a protein-chameleon was proposed, according to which the structure of -synuclein to a dramatic degree depends on the environment: the choice between all the mentioned above conformations is determined by the peculiarities of protein surroundings [82]. Functionally, -synuclein is an example of disordered hub, as in a case-by-case studies, this protein was shown to interact with at least 50 ligands and other proteins [78], whereas a recent proteomic analysis identified 587 proteins involved in the formation of complexes with -synuclein in the dopaminergic cells, with 141 proteins displaying significant changes in their relative abundance (increase or decrease) after these cell were treated with rotenone [83].

Huntingtin is a large protein with an estimated molecular mass of 350 kDa, which contains a polyglutamine tract near its N terminus expansion of which causes Huntington's disease [84]. Huntington's disease is a member of the family of neurodegenerative diseases associated with the expansion of a CAG repeat in the gene which is translated into the extension of the polyglutamine (polyQ) tract in the corresponding protein. The polyQ repeat varies between 16 and 37 residues in healthy individuals, and individuals who are afflicted by disease have repeats of >38 residues. The mechanistic hypothesis linking CAG repeat expansion to toxicity involves the tendency of longer polyQ sequences, regardless of protein context, to form insoluble aggregates [85][86][87][88][89][90][91][92][93]. The far-UV CD spectra of polyQ peptides with various repeat lengths were shown to be nearly identical and were consistent with a high degree of random coil structure, suggesting that the length-dependence of disease is not related to a conformational change in the monomeric states of expanded polyQ sequences [92]. In contrast, there was a dramatic acceleration in the spontaneous formation of ordered, amyloidlike aggregates for poly(Gln) peptides with repeat lengths of greater than 37 residues.

The N terminus of wild-type huntingtin interacts with proteins involved in nuclear functions, including HYPA/ FBP-11, which functions in pre-mRNA processing (splicesome function) [94], nuclear receptor co-repressor protein (NCoR) [95], which plays a role in the repression of gene activity, and p53 [96], a tumor suppressor involved in regulation of the cell cycle. Full-length huntingtin contains candidate binding sites for other proteins with nuclear functions. Huntingtin contains a PXDLS motif, a candidate-binding site for the transcriptional corepressor C-terminal binding protein (CtBP) [97], suggesting that huntingtin may play a role in transcriptional repression.

The p53 protein is a transcription factor that targets genes involved in cell cycle regulation and apoptosis, among other functions [98]. p53 is at the center of a large signaling network, regulating expression of genes involved in such cellular processes as cell cycle progression, apoptosis induction, DNA repair, response to cellular stress, etc. [99]. When p53 function is lost, either directly through mutations or indirectly through several other mechanisms, the cell often undergoes cancerous transformation [100]. In fact, it is believed that all human cancers exhibit defects in the p53-signaling pathway [101]. p53 is considered as the most commonly mutated tumour-suppressor gene in human cancers [102]. In roughly half of all cancer cases the p53 gene is mutated [100]. Typically these are missense mutations within the DNA-binding core domain resulting in the expression of a protein with aberrant function. Among missense mutations, there are hotspot mutations at four codons (175, 248, 249 and 273), which together account for over 25% of all missense mutations identified in human cancers [103]. Cancers showing mutations in p53 are found in colon, lung, esophagus, breast, liver, brain, reticuloendothelial tissues and hemopoietic tissues [100]. For these reasons, a loss of p53 function is believed to be a major factor in cancer development [100] and this protein has attracted significant attention of cancer researchers. A database of p53 point mutations was created http://www.iarc.fr/ p53Ulomepage.htm, which currently is the largest singlelocus mutation database, containing more than 10,000 somatic mutations identified by sequencing [103].

There are three structural domains in p53: N-terminal translational activation domain, central DNA binding domain, and C-terminal tetramerization and regulatory domain. The analysis of the intrinsic order-disorder state in these revealed that the DNA binding domain is intrinsically structured, whereas the terminal domains are intrinsically disordered [104,105]. It has been shown that p53 induces or inhibits over 150 genes, including p21, GADD45, MDM2, IGFBP3, and BAX [106]. At the transactivation region, p53 interacts with TFIID, TFIIH, Mdm2, RPA, CBP/p300 and CSN5/Jab1 [99]. At the C-terminal domain, it interacts with GSK3, PARP-1, TAF1, TRRAP, hGcn5, TAF, 14-3-3, S100B() and many other proteins [99]. Overall, ~70% of the interactions between p53 and its binding partners are mediated by IDRs in p53 [28]. A bias toward intrinsic disorder is even more pronounced in the sites of posttranslational modifications, with 86%, 90%, and 100% of observed acetylation, phosphorylation, and protein conjugation sites, respectively, found in IDRs [28]. This concentration of functional elements within IRDs comparing to just 29% of the residues being disordered [26,28] clearly shows that p53 extensively utilizes IDRs to mediate and modulate interactions with other proteins.

BRCA1

About 5%-10% of breast cancer and ovarian cancer are hereditary and 30%-50% of these are due to the autosomal dominant mutations in the susceptibility genes, BRCA1 and BRCA2 [107]. In both cases the variants are distributed uniformly along the entire coding region and intronic sequences flanking each exon [108]. Women with the BRCA1 mutations are susceptible to the development of a breast cancer before age 35-40 and of an ovarian cancer with a probability rate of, respectively, 45%-60% and 20%-40%. Women carrying BRCA2 mutations present a 25%-40% risk of breast cancer development and a 10%-20% risk of an ovarian cancer development [108].

BRCA1 participates in many different cellular pathways, including transcription, apoptosis and DNA repair, through direct or indirect interaction with a variety of partners [109]. It has multiple alternatively spliced isoforms. One of the most studied BRCA1 isoform has 1863 amino acids and comprises a long highly disordered central region flanked by ordered domains at the two termini. At the N-terminus is a RING finger domain of 103 residues. This domain is reported to form a heterodimer with BARD1 (BRCA1 associated RING domain 1) and to bind to the ubiquitin carboxy-terminal hydrolase BAP1. At the C-terminus are two tandem copies of the BRCA1 C-terminal domain (BRCT) with 218 total residues for the two domains. These two domains are reported to bind with transcriptional activators and repressors like CtlP.

The structural characterization by various spectroscopic techniques revealed that the 1500 amino acid long central region of BRCA1 is completely disordered [110]. However this disordered central region contains molecular recognition domains for both DNA and several protein binding partners, including tumor suppressors such as p53, retinoblastoma protein (RB) and BRCA2; oncogenes like c-Myc and JunB; DNA damage repair proteins such as Rad50 and Rad51; and the Fanconi anemia protein (FANCA) [110]. Importantly, BRCA1 was shown to have at least 24 alternatively spliced isoforms [111]. Alternative splicing was shown to affect mostly central IDR of BRCA1 modulating its functionality by removing different functional domains [32].

EWS-FLI1 fusion protein

Ewing's sarcoma family of tumors is a set of highly malignant tumors of bone and soft tissue that occur in children, adolescents, and young adults. These tumors share a recurrent and specific t (

Hubness and intrinsic disorder in human diseasome

Linear regression of disorder content with respect to number of related diseases, number of related disease classes, and gene degree, shows that the correlation between disorder content and these graph-related gene features are positive and significant. The very low R 2 coefficient tells us that disorder content cannot be predicted from these features (which was never our intention), but that the positive correlations should be observed as trends. Two genes with the highest number of related diseases are PAX6, encoding a developmentally regulated transcription factor paired box protein 6 (Pax-6), which is related to 9 ophthalmological and one developmental disorder, and TP53, a well-studied gene encoding another transcription factor p53, that is involved in 11 different forms of cancer. Some peculiarities of the p53 structure and functions as well as a role of ID in function of this protein were already discussed (see above), whereas a brief overview of the Pax-6 protein is presented below. Pax-6 is a member of a family of developmentally regulated transcription factors that includes at least 8 members expressed in temporally and spatially restricted patterns during development and have been implicated in a number of human congenital disorders, as well as in tum-origenesis [129]. These proteins are characterized by the presence of a specific DNA-binding domain, termed the paired domain. They are highly conserved across millions of years of evolution and human PAX-6 gene is identical to that of axolotol [129]. CD and NMR structural analyses of the purified Pax-6 reveal that it is largely unstuctured in solution. However, upon binding to the recognition DNA sequence, the Pax-6 folds and displays CD spectroscopic evidence of significant -helical structure [129].

A number of related disease classes and gene degree are features related to whether a gene/protein is a hub. The observed trends in predicted disorder content provide additional support for the hypothesis that hub proteins are more likely to be disordered, to accommodate the various interactions and functions they are involved with [26]. All three graph-related gene features are related to the partition of the HDN/DGN graph into one large connected component and a series of small connected components. Genes for which any of the three graph-related features is a high number belong to the large component.

Since such genes are more likely to be disordered, they contribute to the difference in disorder content between large component and small components. This difference is particularly significant for genes related to metabolic diseases. More than 60% of metabolic disease genes that belong to the small components have disorder content in the 0-20% range, and further 30+% have disorder content in the 20-30% range. On the other hand, 25% metabolic disease proteins that belong to the large component have disorder content higher than 40%, which is lower when compared to other disease proteins, but substantially higher than the level of ID in metabolic disease proteins in the small component. Of note, most of metabolic disease genes in the large component are also related to disease from other classes.

The difference in disorder content between one large connected component of HDN/DGN and remaining small connected components has to be observed with caution. The connectivity of HDN/DGN is influenced heavily by small components. Only one link between a gene/disease in the large component and a disease/gene in some small components that has not yet been established, but is discovered in the future can change the partition completely, by leading to inclusion of that whole small component into the large component.

Alternative splicing, intrinsic disorder and human genetic diseases

Prediction of intrinsic disorder in proteins encoded by genes with alternative splicing shows that AS regions have a much higher predicted disorder content than the whole protein sequences. This is in agreement with previous observations [32]. No difference was observed in disorder content for AS regions in disease and non-disease genes/ proteins; the distributions were almost identical. However, alternative splicing can be observed as an important link between diseases and intrinsic disorder, as several disease classes have significantly higher fraction of genes with multiple isoforms; i.e., with AS regions. The presence of AS regions in such genes is associated with increased disorder content. Distributions of disorder content in AS regions were fairly similar across various genes, except for three classes. DEVE and NEUR have a very high fraction of highly disordered AS regions (disorder content 80-100%). This fact might be related to the functionality of proteins involved in these diseases (see above). AS regions in META genes are much less disordered than AS regions in other disease classes, just like whole META gene sequences are much less disordered than other disease genes.

Abundance of -MoRFs in proteins associated with human genetic diseases

IDRs frequently participate in protein-protein interactions and molecular recognitions [1,5,10,22,30,34,130]. Many IDPs and IDRs undergo disorder-to-order transitions upon binding, which is crucial for recognition, regulation, and signaling [1,4,14,28,[34][35][36][37][131][132][133]. A recent confounding observation is that not all specific interactions between intrinsically disordered proteins and their partners are necessarily accompanied by the disorder-to-order transitions, but may somehow remain unstructured even after binding [134][135][136][137][138]. Nevertheless, a correlation has been established between the specific pattern in the PONDR ® VL-XT curve and the ability of a given short disordered regions to undergo a disorder-to-helix transition upon binding [139]. Based on these specific features, a predictor helix-forming MoRFs was recently developed [34,37]. Not all helix forming MoRF regions share these same features, and some MoRFs form or irregular structure rather than the -helix [35,36]. A further complication is that MoRFs can exhibit partner-dependent structures, with at least one example morphing into helix, sheet, or irregular structure, depending on the partner [28]. Overall, therefore, these predicted MoRFs represent only fractions of the total numbers of MoRFs for each organism.

The application of the -MoRF predictor to various datasets reveals that helix forming molecular recognition features are highly abundant in proteins associated with all human genetic diseases as well as in proteins encoded by disease genes and by all human genes, suggesting the existence of extensive interaction networks. In the HUM set, 57.9% of human genes contain -MoRFs. In the DIS set, 54.4% of all disease-associated genes contain -MoRFs, with significant variation between various disease classes, ranging from 26.0% in metabolic diseases and 27.3% in nutritional disorders to 73.4% in cancer and 78.6% in skeletal diseases. In most disease classes some long, highly disordered proteins have multiple predicted -MoRF regions ( Interestingly, fractions of proteins with predicted -MoRF regions were highly correlated with the content of predicted disorder in a given dataset (correlation coefficient is ~0.89). This suggests that the major function of IDRs in the proteins from analyzed datasets is protein-protein interaction. -MoRFs, being disordered in the unbound state and gaining -helical structure upon interaction with binding partners, suit ideally this function. In fact, it has been proposed that that the involvement of IDRs in protein-protein interactions have several advantages [27], including: (i) Decoupled specificity and strength of binding (high-specificity-low-affinity interactions); (ii) Increased speed of interaction due to greater capture radius and the ability to spatially search interaction space; (iii) Efficient regulation via rapid degradation; (iv) Increased interaction (surface) area per residue; (v) Strengthened encounter complex (less stringent spatial orientation requirements); (vi) A single disordered region may bind to several structurally diverse partners; (vii) Many (structured) proteins may bind a single disordered region; (viii) Less sterically restricted to allow elongation of binding area; (ix) Efficient regulation via posttranslational modification; (x) Ease of regulation/redirection by alternative splicing; (xi) Overlapping binding sites due to extended linear conformation; (xii) High evolutionary rate; (xiii) Flexibility that allows masking (or not) of interaction sites or allow interaction between bound partners. Many of these features are specific properties of -MoRFs.

Abundance of -MoRFs in alternative spliced regions of proteins from human diseasome

Interestingly, our analysis revealed that -MoRFs are abundantly present in alternatively spliced regions of proteins from some human genetic diseases. This observation is very important as it sheds some light on the potential functional repertoire of alternatively spliced regions. In several diseases, these regions play a crucial role in pro-tein-protein interaction, as they are enriched in molecular recognition features.

Concluding remarks

Intrinsically disordered proteins are highly abundant in nature. Although they lack stable tertiary and/or secondary structure under physiological conditions in vitro, IDPs carry out a number of crucial biological functions, being involved in regulation, recognition, signaling and control. The functional repertoire of IDPs complements the functions of ordered proteins. Earlier studies revealed that many IDPs are associated with various human diseases, including cancer, cardiovascular disease, amyloidoses, neurodegenerative diseases, diabetes and others, emphasizing the existence of intriguing interconnections between IDPs, cell signaling and human diseases. Based on these observations, the "disorder in disorders" or D 2 concept was introduced.

Here, a large-scale analysis of the abundance of intrinsic disorder in human genetic diseases joined into the human diseasome [47] was performed. This analysis uncovered an unfoldome (an IDP-containing subset of a proteome) associated with human genetic diseases and revealed several interesting peculiarities. Particularly, we are showing here that proteins associated with various human genetic diseases are enriched in intrinsic disorder with the IDP content being markedly different for different genetic diseases. The diseasome possesses a high level of MoRFs, whose abundance correlates with the intrinsic disorder level. Alternative splicing is commonly present in several genetic diseases. Alternatively spliced regions in corresponding proteins are predicted to be highly disordered and in some diseases contain a significant number of MoRFs. The various diseasome graph-related properties are correlated with the levels of intrinsic disorder (hub proteins are generally more disordered). These data were used to build the unfoldome for the diseaseome.

Figure 1
1Comparison of disorder content distributions in disease classes and human gene class using boxplots. The 22 disease classes are sorted according to their disorder content medians. The boxes in the boxplot represent the first quartile (left edge), median (line in the middle), and third quartile (right edge); the whiskers extend to the lowest/highest values within the 1.5 IQR interval from the box (IQR is the range between the first and the third quartile), while the + signs represent the outliers. Medians for two classes can be compared by looking at the notches at their median lines; if the notches do not overlap, the medians are different at the significance level  = 0.05.

Figure 3
3Pairwise comparison of disorder content medians for disease classes and human gene class. Filled squares represent pairs for which adjusted Wilcoxon rank sum test p-values are smaller than  = 0.05 (p-values are adjusted for false discovery rate control with Benjamini-Yekutieli method). Squares are filled black if the median for the row class is greater than the median for the column class, or gray if the median for the row class is smaller than the median for the column class.

Figure 12
12compiles the -MoRF prediction data and shows the fractions of genes with predicted -MoRFs and the densities -MoRFs (number of -MoRFs per residue) for all disease classes, as well as for sets of all disease genes and all human genes. The overall fractions of disordered residues are included for comparison. The fractions of Comparison of fractions of disease genes in the large compo-nent and the small components of the DGN Figure 7 Comparison of fractions of disease genes in the large component and the small components of the DGN. The classes with the + signs after their acronyms are significantly overrepresented in the big component; the classes with the -signs after their acronyms are significantly underrepresented in the big component. The error bars represent one standard deviation or 68.2% confidence interval.

Figure 13 compares
13the overall density of predicted -MoRFs versus density of predicted -MoRFs in AS regions for the 25 classes. The differences between densities of MoRFs (overall vs. AS regions) are significant for the majority of classes (listed by increasing p-values: HUM, NEUR, META, CANC, DIS, GI, DEVE, IMMU, ENDO, RESP, BONE, DERM, MUSC, CARD, HEMA, ENT), borderline significant for NUTR and OPHT, and not significant for the remaining classes (RENA, MCD, UNCL, SKEL, MULT). Two classes (PSYC, UNCL) have no -MoRFs predicted in AS regions (while genes in both classes have very small number of residues in AS regions, for PSYC class this difference in densities is still statistically significant). Finally,

Figure 15
15compares the results of the CH-plot and CDF analyses by showing the distributions of proteins in each Comparison of disorder content distributions for AS regions in various classes of human genes Figure 11 Comparison of disorder content distributions for AS regions in various classes of human genes. Series #1 and #2 represent disorder content distribution for whole human gene sequences and AS regions in human genes. Series #3-6 represent disorder content distributions for AS regions in: 3) developmental, 4) neurological, 5) hematological, and 6) metabolic disease classes. The error bars represent one standard deviation or 68.2% confidence interval.

Figure 13
13Comparison of overall density of predicted MoRFs vs density of predicted MoRFs in AS regions for 25 classes/sets. The error bars represent one standard deviation or 68.2% confidence interval. Comparison of fractions of genes with predicted -MoRFs and densities of -MoRFs with fractions of disordered resi-dues Figure 12 Comparison of fractions of genes with predicted -MoRFs and densities of -MoRFs with fractions of disordered residues. The plot on the left compares fractions of genes with predicted -MoRFs (top/first series) with fractions of disordered residues (bottom/second series). The plot on the right compares densities of -MoRFs (top/first series) with fractions of disordered residues (bottom/second series). In both plots the series are shown with different scales, such that the values for HUM set are aligned. The error bars represent one standard deviation or 68.2% confidence interval.

6%), GNAS (3, 1323, 74.9%), FAS (3, 376, 59.3%), MXI1 (3, 320, 88.8%), ATM (2, 3056, 21.1%), DLC1 (2, 1591, 62.7%), PML (2, 1225, 63.9%), ABL1 (2, 1175, 59.2%), AR (2, 927, 60.0%), CHEK2 (2, 586, 37.2%), CASP8 (2, 548, 38.5%), PARK2 (2, 465, 37.0%), SMARCB1 (2, 385, 37.1%), SSX2 (2, 255, 77.3%) DEVE (36/53 ~67.9%) NSD1 (2, 2706, 77.3%), UBE3A (2, 882, 27.8%), PVRL1 (2, 658, 37.2%), TGIF1 (2, 424, 77.4%) MCD (142/209 ~67.9%) MITF (5, 598, 75.6%), GNAS (3, 1323, 74.9%), PITX2 (3, 385, 81.0%), NSD1 (2, 2706, 77.3%), ATRX (2, 2492, 72.2%), COL11A1 (2, 1857, 81.3%), COL18A1 (2, 1551, 74.4%), L1CAM (2, 1257, 27.0%), USH1C (2, 926, 58.4%), FGFR1 (2, 901, 38.5%), HPS4 (2, 783, 42.9%), KCNQ1 (2, 718, 42.3%), PVRL1 (2, 658, 37.2%), DTNBP1 (2, 383, 81.7%) CARD (53/96 ~55.2%) DMD (7, 3771, 54.4%), DTNA (4, 767, 61.0%), KCNH2 (2, 1283, 46.0%), KCNQ1 (2, 718, 42.3%), EYA4 (2, 665, 61.7%), TPM1 (2, 443, 100.0%) MUSC (45/68 ~66.2%) DMD (7, 3771, 54.4%), PLEC1 (3, 4904, 56.6%), COL6A3 (2, 3177, 28.6%), AR (2, 927, 60.0%), CHAT (2, 748, 34.1%), TPM3 (2, 378, 97.6%) IMMU (59/115 ~51.3%) FAS (3, 376, 59.3%), ATM (2, 3056, 21.1%), PTPRC (2, 1307, 42.2%), CASP8 (2, 548, 38.5%), PARK2 (2, 465, 37.0%), UNG (2, 348, 44.5%) OPHT (68/120 ~56.7%) OPA1 (3, 1015, 36.1%), PITX2 (3, 385, 81.0%), PIP5K3 (2, 2108, 52.0%), EYA1 (2, 600, 60.5%) CTD (32/51 ~62.7%) --ENDO (54/96 ~56.3%) GNAS (3, 1323, 74.9%), AR (2, 927, 60.0%), HNF4A (2, 531, 51.6%), GCK (2, 495, 35.6%) NEUR (154/254 ~60.6%) COLQ (3, 622, 76.0%), PTPRC (2, 1307, 42.2%), L1CAM (2, 1257, 27.0%), KCNQ2 (2, 892, 55.5%), FOXP2 (2, 740, 78.9%), MTMR2 (2, 643, 31.6%), SPAST (2, 616, 51.1%), EYA1 (2, 600, 60.5%), NR4A2 (2, 599, 55.8%), EIF2B4 (2, 555, 47.2%), CACNB4 (2, 538, 60.0%), OPRM1 (2, 492, 34.6%), CCM2 (2, 475, 54.1%), PARK2 (2, 465, 37.0%), DCX (2, 446, 57.8%), DRD2 (2, 443, 39.1%), PNKD (2, 440, 39.5%), ATXN3 (2, 370, 61.6%), FGF14 (2, 316, 46.2%) PSYC (18/30 ~60.0%) --ENT (26/44 ~59.1%) OTOF (2, 2100, 33.9%), USH1C (2, 926, 58.4%), KCNQ4 (2, 695, 41.3%), EYA4 (2, 665, 61.7%) RESP (12/34 ~35.3%) GDNF (2, 230, 51.3%) RENA (35/58 ~60.3%) --HEMA (58/146 ~39.7%) CD44 (3, 807, 76.0%), ATRX (2, 2492, 72.2%), ANK1 (2, 2001, 40.9%), ADAMTS13 (2, 1497, 48.6%), EPB41 (2, 850, 66.2%), AMPD3 (2, 781, 36.9%), IGLL1 (2, 228, 86.4%) NUTR (6/22 ~27.3%) --GI (16/34 ~47.1%) GDNF (2, 230, 51.3%) UNCL (12/29 ~41.4%) --META (75/289 ~26.0%) GCK (2, 495, 35.6%), HFE2 (2, 426, 42.0%)

Table 1 :
1Disease class names and acronyms, number of diseases and number of genes related to disease classes.Class name
Acronym
Number of diseases
%
(of 1284)

Number of genes
%
(of 1751)

Skeletal
SKEL
64
4.98%
56
3.20%
Bone
BONE
30
2.34%
44
2.51%
Dermatological
DERM
48
3.74%
80
4.57%
Cancer
CANC
113
8.80%
207
11.82%
Developmental
DEVE
37
2.88%
53
3.03%
Multi-class disease
MCD
155
12.07%
209
11.94%
Cardiovascular
CARD
41
3.19%
96
5.48%
Muscular
MUSC
31
2.41%
68
3.88%
Immunological
IMMU
69
5.37%
115
6.57%
Ophthamological
OPHT
62
4.83%
120
6.85%
Connective tissue disorder
CTD
28
2.18%
51
2.91%
Endocrine
ENDO
56
4.36%
96
5.48%
Neurological
NEUR
117
9.11%
254
14.51%
Psychiatric
PSYC
17
1.32%
30
1.71%
Ear, Nose, Throat
ENT
6
0.47%
44
2.51%
Respiratory
RESP
13
1.01%
34
1.94%
Renal
RENA
36
2.80%
58
3.31%
Hematological
HEMA
88
6.85%
146
8.34%
Nutritional
NUTR
4
0.31%
22
1.26%
Gastrointestinal
GI
23
1.79%
34
1.94%
Unclassified
UNCL
31
2.41%
29
1.66%
Metabolic
META
215
16.74%
289
16.50%
Multiple class genes
MULT
295
16.85%
Disease genes
DIS
1751
100.00%
Human genes
HUM
18109

content median than the human gene set with a BY false discovery rate of 0.05. Several other classes have low p-values in comparison with human gene set, but the adjustment for the BY method pushes them above the 0.05 limit. Note that adjusted p-values would be ~3.7 times smaller if we used the Benjamini-Hochberg false discovery methodwhich shows the
top 15 p-values and BY adjusted p-values for comparison
of disease classes with disease gene (DIS) set, as well as for
comparison of disease classes with human gene (HUM)
set. Several other classes, besides the one indicated in Fig-
ure 2, can be considered to have disorder content medians
significantly different from the DIS and HUM classes,
depending on how strict the comparisons are to be. For
example, cancer gene class has (borderline) significantly
different disorder

Table 2 :
2Comparison of disorder content medians of disease classes with disease gene set (DIS) and with human gene set (HUM).Both the p-values and the adjusted p-values (for Benjamini-Yekutieli FDR control method) are listed in the table.Comparison with DIS
Comparison with HUM

p-value
BY p-value
p-value
BY p-value

META 9.10·10 -31 7.81·10 -29
META 1.38·10 -50 1.25·10 -48
CANC 9.76·10 -09 4.19·10 -07
DIS
6.16·10 -15 2.79·10 -13
SKEL
3.92·10 -05 0.001123
HEMA 7.13·10 -08 2.15·10 -06
MCD
0.000548 0.011771
UNCL 0.000192 0.004349
DERM 0.001852 0.031810
CANC 0.002397 0.043445
HEMA 0.003684 0.052740
NUTR 0.007141 0.107855
UNCL 0.004152 0.050941
SKEL
0.011080 0.143441
DEVE
0.008386 0.090036
GI
0.015816 0.179167
NEUR 0.033455 0.319267
IMMU 0.016768 0.168843
BONE 0.042742 0.367102
RENA 0.026136 0.236856
NUTR 0.063282 0.494113
RESP
0.093824 0.772967
MULT 0.090375 0.646849
MULT 0.105644 0.797813
MUSC 0.122463 0.809091
DERM 0.178919 1.247247
GI
0.130811 0.802516
ENT
0.195823 1.267578
ENDO 0.164391 0.941288
DEVE
0.208293 1.258409

Table 3 :
3Comparison of densities of predicted -MoRFs in AS regions and complete genes.The quotients of density of predicted -MoRFs in AS regions over overall density of predicted -MoRFs and the p-values for comparison of corresponding densities.Class acronym
Density of MoRFs in AS/Overall density of MoRFs
p-value for comparison of densities of MoRFs

GI
4.68
3.84·10 -021
META
3.71
1.16·10 -061
RESP
2.65
1.19·10 -010
NEUR
2.45
2.54·10 -076
DEVE
2.37
4.76·10 -017
BONE
2.10
1.28·10 -010
NUTR
2.06
0.017694
DERM
1.75
5.11·10 -010
ENDO
1.71
2.93·10 -011
IMMU
1.69
1.92·10 -013
CANC
1.68
1.07·10 -031
ENT
1.50
0.00079422
HEMA
1.48
4.65·10 -007
HUM
1.30
6.12·10 -233
DIS
1.27
1.10·10 -030
RENA
1.21
0.55174
MCD
1.08
0.55638
SKEL
1.06
2.8735
MULT
0.99
3.0621
OPHT
0.78
0.049566
CARD
0.72
3.30·10 -007
MUSC
0.66
8.89·10 -009
CTD
0.32
3.78·10 -006
PSYC
0
2.07·10 -005
UNCL
0
0.56375

Table 4 :
4Disease-related genes with multiple predicted -MoRFs.SKEL (44/56 ~78.6%)
FGFR1 (2, 901, 38.5%)
BONE (29/44 ~65.9%)
GNAS (3, 1323, 74.9%), COL9A1 (2, 945, 76.5%), AMELX (2, 205, 66.8%)
DERM (55/80 ~68.8%)
EDA
(4, 460, 63.9%), PLEC1 (3, 4904, 56.6%), ADAR (2, 1226, 53.6%), SLC39A4 (2, 686, 36.3%), PVRL1 (2, 658, 37.2%)
CANC (152/207 ~73.4%) BRCA1 (3, 1864, 80.

to form EWS-FLI1 fusion protein, a chimeric transcription factor. EWS-FLI1 is expressed only in tumor cells and its function is required for the malignant phenotype of Ewing's sarcoma family of tumors [113]. EWS-FLI1 retains the Ets DNA binding domain from FLI1 and modulates a diverse group of target genes by binding to specific promoters including transforming growth factor- receptor type-II [114], p21 (WAF1/CIP1) [115], PTPL1 [116], Id2 [117], andtenascin-C [118], EAT-2 [119], mE2C [120], manic fringe [121], c-myc [122], platelet-derived growth factor C [123], p57KIP [124], and PIM-3 [125]. EWS-FLI1 also regulates gene expression by modulating RNA splicing as shown by alteration of an E1A splice site and interaction with U1C [126,127]. Despite these numerous activities, the EWS-FLI1 fusion protein was shown to approach a largely unfolded conformation under native conditions [128].11;22) (q24;q12) chromosome
translocation [112], which combines the N-terminus of
EWS (residues 1-264) from chromosome 22 with the C-
terminus of FLI1 (232 carboxy-terminal residues) from
chromosome 11

Table 4 )
4that may potentially serve as binding sites for multiple proteins. For example, DMD from CARD/MUSC (7 predicted -MoRFs, 3771 amino acids, 54.4% disorder content); MITF from MCD (5, 598,

(page number not for citation purposes)
(page number not for citation purposes) for multiple protein-protein and protein-DNA interactions? J Mol Biol 2005, 345(2):275-287. 111. Orban TI, Olah E: Emerging roles of BRCA1 alternative splic-
AcknowledgementsThis work was supported in part by the grants R01 LM007688-01A1 (to A.K.D and V.N.U.) and GM071714-01A2 (to A.K.D and V.N.U.) from the National Institutes of Health and the Program of the Russian Academy of Sciences for the "Molecular and cellular biology" (to V. N. U.). We gratefully acknowledge the support of the IUPUI Signature Centers Initiative.This article has been published as part of BMC Genomics Volume 10 Supplement 1, 2009: The 2008 International Conference on Bioinformatics & Computational Biology (BIOCOMP'08). The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2164/ 10?issue=S1.Competing interestsThe authors declare that they have no competing interests.Authors' contributionsVNU was involved in design and planning of all the experiments, contributed to the manuscript writing, revised the final version and headed the project. UM performed the computational analysis, designed figures and drafted the manuscript. CJO performed the computational analysis and contributed to the manuscript writing. ZO and AKD were involved in design and planning of all the experiments and contributed to the manuscript writing. All authors have read and approved the final manuscript.
. A K Dunker, J D Lawson, C J Brown, R M Williams, P Romero, J S Oh, C J Oldfield, A M Campen, C M Ratliff, K W Hipps, Intrinsically disordered protein. J Mol Graph Model. 191Dunker AK, Lawson JD, Brown CJ, Williams RM, Romero P, Oh JS, Oldfield CJ, Campen AM, Ratliff CM, Hipps KW, et al.: Intrinsically disordered protein. J Mol Graph Model 2001, 19(1):26-59.

Structural studies of tau protein and Alzheimer paired helical filaments show no evidence for beta-structure. O Schweers, E Schonbrunn-Hanebeck, A Marx, E Mandelkow, J Biol Chem. 26939Schweers O, Schonbrunn-Hanebeck E, Marx A, Mandelkow E: Struc- tural studies of tau protein and Alzheimer paired helical fila- ments show no evidence for beta-structure. J Biol Chem 1994, 269(39):24290-24297.

NACP, a protein implicated in Alzheimer's disease and learning, is natively unfolded. P H Weinreb, W Zhen, A W Poon, K A Conway, P T LansburyJr, Biochemistry. 3543Weinreb PH, Zhen W, Poon AW, Conway KA, Lansbury PT Jr: NACP, a protein implicated in Alzheimer's disease and learning, is natively unfolded. Biochemistry 1996, 35(43):13709-13715.

Intrinsically unstructured proteins: reassessing the protein structure-function paradigm. P E Wright, H J Dyson, J Mol Biol. 2932Wright PE, Dyson HJ: Intrinsically unstructured proteins: re- assessing the protein structure-function paradigm. J Mol Biol 1999, 293(2):321-331.

Natively disordered proteins. G W Daughdrill, G J Pielak, V N Uversky, M S Cortese, A K Dunker, Handbook of Protein Folding Edited by: Buchner J, Kiefhaber T. Weinheim, GermanyWiley-VCH, Verlag GmbH & Co. KGaADaughdrill GW, Pielak GJ, Uversky VN, Cortese MS, Dunker AK: Natively disordered proteins. In Handbook of Protein Folding Edited by: Buchner J, Kiefhaber T. Weinheim, Germany: Wiley-VCH, Verlag GmbH & Co. KGaA; 2005:271-353.

Protein folding revisited. A polypeptide chain at the folding-misfolding-nonfolding cross-roads: which way to go?. V N Uversky, Cell Mol Life Sci. 609Uversky VN: Protein folding revisited. A polypeptide chain at the folding-misfolding-nonfolding cross-roads: which way to go? Cell Mol Life Sci 2003, 60(9):1852-1871.

Fluorescence correlation spectroscopy shows that monomeric polyglutamine molecules form collapsed structures in aqueous solutions. S L Crick, M Jayaraman, C Frieden, R Wetzel, R V Pappu, Proc Natl Acad Sci. 10345Crick SL, Jayaraman M, Frieden C, Wetzel R, Pappu RV: Fluores- cence correlation spectroscopy shows that monomeric pol- yglutamine molecules form collapsed structures in aqueous solutions. Proc Natl Acad Sci USA 2006, 103(45):16764-16769.

Role of backbone-solvent interactions in determining conformational equilibria of intrinsically disordered proteins. H T Tran, A Mao, R V Pappu, J Am Chem Soc. 13023Tran HT, Mao A, Pappu RV: Role of backbone-solvent interac- tions in determining conformational equilibria of intrinsi- cally disordered proteins. J Am Chem Soc 2008, 130(23):7380-7392.

Sequence complexity of disordered protein. P Romero, Z Obradovic, X Li, E C Garner, C J Brown, A K Dunker, Proteins. 421Romero P, Obradovic Z, Li X, Garner EC, Brown CJ, Dunker AK: Sequence complexity of disordered protein. Proteins 2001, 42(1):38-48.

Intrinsic disorder and functional proteomics. Biophys J. P Radivojac, L M Iakoucheva, C J Oldfield, Z Obradovic, V N Uversky, A K Dunker, 92Radivojac P, Iakoucheva LM, Oldfield CJ, Obradovic Z, Uversky VN, Dunker AK: Intrinsic disorder and functional proteomics. Bio- phys J 2007, 92(5):1439-1456.

The protein folding problem. T E Creighton, Science. 2404850Creighton TE: The protein folding problem. Science 1988, 240(4850):267-344.

The protein non-folding problem: amino acid determinants of intrinsic order and disorder. R M Williams, Z Obradovic, V Mathura, W Braun, E C Garner, J Young, S Takayama, C J Brown, A K Dunker, Pac Symp Biocomput. Williams RM, Obradovic Z, Mathura V, Braun W, Garner EC, Young J, Takayama S, Brown CJ, Dunker AK: The protein non-folding problem: amino acid determinants of intrinsic order and dis- order. Pac Symp Biocomput 2001:89-100.

Natively unfolded proteins: a point where biology waits for physics. V N Uversky, Protein Sci. 114Uversky VN: Natively unfolded proteins: a point where biol- ogy waits for physics. Protein Sci 2002, 11(4):739-756.

Why are "natively unfolded" proteins unstructured under physiologic conditions? Proteins. V N Uversky, J R Gillespie, A L Fink, 41Uversky VN, Gillespie JR, Fink AL: Why are "natively unfolded" proteins unstructured under physiologic conditions? Proteins 2000, 41(3):415-427.

IUPred: web server for the prediction of intrinsically unstructured regions of proteins based on estimated energy content. Z Dosztanyi, V Csizmok, P Tompa, I Simon, Bioinformatics. 2116Dosztanyi Z, Csizmok V, Tompa P, Simon I: IUPred: web server for the prediction of intrinsically unstructured regions of proteins based on estimated energy content. Bioinformatics 2005, 21(16):3433-3434.

Comparing and combining predictors of mostly disordered proteins. C J Oldfield, Y Cheng, M S Cortese, C J Brown, V N Uversky, A K Dunker, Biochemistry. 446Oldfield CJ, Cheng Y, Cortese MS, Brown CJ, Uversky VN, Dunker AK: Comparing and combining predictors of mostly disor- dered proteins. Biochemistry 2005, 44(6):1989-2000.

A practical overview of protein disorder prediction methods. F Ferron, S Longhi, B Canard, D Karlin, Proteins. 651Ferron F, Longhi S, Canard B, Karlin D: A practical overview of protein disorder prediction methods. Proteins 2006, 65(1):1-14.

DisProt: the Database of Disordered Proteins. M Sickmeier, J A Hamilton, T Legall, V Vacic, M S Cortese, A Tantos, B Szabo, P Tompa, J Chen, V N Uversky, Nucleic Acids Res. Sickmeier M, Hamilton JA, LeGall T, Vacic V, Cortese MS, Tantos A, Szabo B, Tompa P, Chen J, Uversky VN, et al.: DisProt: the Data- base of Disordered Proteins. Nucleic Acids Res 2007:D786-793.

Identifying disordered regions in proteins from amino acid sequence. P Romero, Z Obradovic, C Kissinger, J E Villafranca, A K Dunker, Proceedings of International Conference on Neural Networks 1997. International Conference on Neural Networks 19971Romero P, Obradovic Z, Kissinger C, Villafranca JE, Dunker AK: Identifying disordered regions in proteins from amino acid sequence. 1997 Proceedings of International Conference on Neural Net- works 1997, 1:90-95.

Thousands of proteins likely to have long disordered regions. P Romero, Z Obradovic, C R Kissinger, J E Villafranca, E Garner, S Guilliot, A K Dunker, Pac Symp Biocomput. Romero P, Obradovic Z, Kissinger CR, Villafranca JE, Garner E, Guil- liot S, Dunker AK: Thousands of proteins likely to have long disordered regions. Pac Symp Biocomput 1998:437-448.

Prediction and functional analysis of native disorder in proteins from the three kingdoms of life. J J Ward, J S Sodhi, L J Mcguffin, B F Buxton, D T Jones, J Mol Biol. 3373Ward JJ, Sodhi JS, McGuffin LJ, Buxton BF, Jones DT: Prediction and functional analysis of native disorder in proteins from the three kingdoms of life. J Mol Biol 2004, 337(3):635-645.

Intrinsic disorder in cell-signaling and cancer-associated proteins. L M Iakoucheva, C J Brown, J D Lawson, Z Obradovic, A K Dunker, J Mol Biol. 3233Iakoucheva LM, Brown CJ, Lawson JD, Obradovic Z, Dunker AK: Intrinsic disorder in cell-signaling and cancer-associated pro- teins. J Mol Biol 2002, 323(3):573-584.

Functional anthology of intrinsic disorder. 1. Biological processes and functions of proteins with long disordered regions. H Xie, S Vucetic, L M Iakoucheva, C J Oldfield, A K Dunker, V N Uversky, Z Obradovic, J Proteome Res. 65Xie H, Vucetic S, Iakoucheva LM, Oldfield CJ, Dunker AK, Uversky VN, Obradovic Z: Functional anthology of intrinsic disorder. 1. Biological processes and functions of proteins with long dis- ordered regions. J Proteome Res 2007, 6(5):1882-1898.

Functional anthology of intrinsic disorder. 2. Cellular components, domains, technical terms, developmental processes, and coding sequence diversities correlated with long disordered regions. S Vucetic, H Xie, L M Iakoucheva, C J Oldfield, A K Dunker, Z Obradovic, V N Uversky, J Proteome Res. 65Vucetic S, Xie H, Iakoucheva LM, Oldfield CJ, Dunker AK, Obradovic Z, Uversky VN: Functional anthology of intrinsic disorder. 2. Cellular components, domains, technical terms, develop- mental processes, and coding sequence diversities corre- lated with long disordered regions. J Proteome Res 2007, 6(5):1899-1916.

Functional anthology of intrinsic disorder. 3. Ligands, post-translational modifications, and diseases associated with intrinsically disordered proteins. H Xie, S Vucetic, L M Iakoucheva, C J Oldfield, A K Dunker, Z Obradovic, V N Uversky, J Proteome Res. 65Xie H, Vucetic S, Iakoucheva LM, Oldfield CJ, Dunker AK, Obradovic Z, Uversky VN: Functional anthology of intrinsic disorder. 3. Ligands, post-translational modifications, and diseases asso- ciated with intrinsically disordered proteins. J Proteome Res 2007, 6(5):1917-1932.

Flexible nets. The roles of intrinsic disorder in protein interaction networks. A K Dunker, M S Cortese, P Romero, L M Iakoucheva, V N Uversky, Febs J. 27220Dunker AK, Cortese MS, Romero P, Iakoucheva LM, Uversky VN: Flexible nets. The roles of intrinsic disorder in protein inter- action networks. Febs J 2005, 272(20):5129-5148.

Keith Dunker A: Intrinsic disorder in scaffold proteins: Getting more from less. M S Cortese, V N Uversky, Prog Biophys Mol Biol. 981Cortese MS, Uversky VN, Keith Dunker A: Intrinsic disorder in scaffold proteins: Getting more from less. Prog Biophys Mol Biol 2008, 98(1):85-106.

Flexible nets: disorder and induced fit in the associations of p53 and 14-3-3 with their partners. C J Oldfield, J Meng, J Y Yang, M Q Yang, V N Uversky, A K Dunker, BMC Genomics. 91Suppl 1Oldfield CJ, Meng J, Yang JY, Yang MQ, Uversky VN, Dunker AK: Flexible nets: disorder and induced fit in the associations of p53 and 14-3-3 with their partners. BMC Genomics 2008, 9(Suppl 1):S1.

Calmodulin signaling: analysis and prediction of a disorder-dependent molecular recognition. P Radivojac, S Vucetic, T R O'connor, V N Uversky, Z Obradovic, A K Dunker, Proteins. 632Radivojac P, Vucetic S, O'Connor TR, Uversky VN, Obradovic Z, Dunker AK: Calmodulin signaling: analysis and prediction of a disorder-dependent molecular recognition. Proteins 2006, 63(2):398-410.

Showing your ID: intrinsic disorder as an ID for recognition, regulation and cell signaling. V N Uversky, C J Oldfield, A K Dunker, J Mol Recognit. 185Uversky VN, Oldfield CJ, Dunker AK: Showing your ID: intrinsic disorder as an ID for recognition, regulation and cell signal- ing. J Mol Recognit 2005, 18(5):343-384.

Protein disorder and the evolution of molecular recognition: theory, predictions and observations. A K Dunker, E Garner, S Guilliot, P Romero, K Albrecht, J Hart, Z Obradovic, C Kissinger, J E Villafranca, Pac Symp Biocomput. Dunker AK, Garner E, Guilliot S, Romero P, Albrecht K, Hart J, Obra- dovic Z, Kissinger C, Villafranca JE: Protein disorder and the evo- lution of molecular recognition: theory, predictions and observations. Pac Symp Biocomput 1998:473-484.

Alternative splicing in concert with protein intrinsic disorder enables increased functional diversity in multicellular organisms. P R Romero, S Zaidi, Y Y Fang, V N Uversky, P Radivojac, C J Oldfield, M S Cortese, M Sickmeier, T Legall, Z Obradovic, Proc Natl Acad Sci. 10322Romero PR, Zaidi S, Fang YY, Uversky VN, Radivojac P, Oldfield CJ, Cortese MS, Sickmeier M, LeGall T, Obradovic Z, et al.: Alternative splicing in concert with protein intrinsic disorder enables increased functional diversity in multicellular organisms. Proc Natl Acad Sci USA 2006, 103(22):8390-8395.

Intrinsic disorder and protein function. A K Dunker, C J Brown, J D Lawson, L M Iakoucheva, Z Obradovic, Biochemistry. 4121Dunker AK, Brown CJ, Lawson JD, Iakoucheva LM, Obradovic Z: Intrinsic disorder and protein function. Biochemistry 2002, 41(21):6573-6582.

Coupled folding and binding with alpha-helix-forming molecular recognition elements. C J Oldfield, Y Cheng, M S Cortese, P Romero, V N Uversky, A K Dunker, Biochemistry. 4437Oldfield CJ, Cheng Y, Cortese MS, Romero P, Uversky VN, Dunker AK: Coupled folding and binding with alpha-helix-forming molecular recognition elements. Biochemistry 2005, 44(37):12454-12470.

Analysis of molecular recognition features (MoRFs). A Mohan, C J Oldfield, P Radivojac, V Vacic, M S Cortese, A K Dunker, V N Uversky, J Mol Biol. 3625Mohan A, Oldfield CJ, Radivojac P, Vacic V, Cortese MS, Dunker AK, Uversky VN: Analysis of molecular recognition features (MoRFs). J Mol Biol 2006, 362(5):1043-1059.

Characterization of molecular recognition features, MoRFs, and their binding partners. V Vacic, C J Oldfield, A Mohan, P Radivojac, M S Cortese, V N Uversky, A K Dunker, J Proteome Res. 66Vacic V, Oldfield CJ, Mohan A, Radivojac P, Cortese MS, Uversky VN, Dunker AK: Characterization of molecular recognition fea- tures, MoRFs, and their binding partners. J Proteome Res 2007, 6(6):2351-2366.

Mining -helix-forming molecular recognition features -MoRFs with cross species sequence alignments. Y Cheng, C J Oldfield, P Romero, V N Uversky, A K Dunker, Biochemistry. 4647Cheng Y, Oldfield CJ, Romero P, Uversky VN, Dunker AK: Mining -helix-forming molecular recognition features -MoRFs with cross species sequence alignments. Biochemistry 2007, 46(47):13468-13477.

Intrinsically disordered proteins in human diseases: Introducing the D2 concept. V N Uversky, C J Oldfield, A K Dunker, Ann Rev Biophys Biomol Structure. 37Uversky VN, Oldfield CJ, Dunker AK: Intrinsically disordered proteins in human diseases: Introducing the D2 concept. Ann Rev Biophys Biomol Structure 2008, 37:215-246.

Abundance of intrinsic disorder in protein associated with cardiovascular disease. Y Cheng, T Legall, C J Oldfield, A K Dunker, V N Uversky, Biochemistry. 4535Cheng Y, LeGall T, Oldfield CJ, Dunker AK, Uversky VN: Abun- dance of intrinsic disorder in protein associated with cardio- vascular disease. Biochemistry 2006, 45(35):10448-10460.

Amyloidogenesis of natively unfolded proteins. V N Uversky, Current Alzheimer Research. 53Uversky VN: Amyloidogenesis of natively unfolded proteins. Current Alzheimer Research 2008, 5(3):260-287.

Intrinsic disorder in proteins associated with neurodegenerative diseases. V N Uversky, Protein Folding and Misfolding: Neurodegenerative Diseases Edited by: Ovádi J, Orosz F. New York, USASpringerUversky VN: Intrinsic disorder in proteins associated with neurodegenerative diseases. In Protein Folding and Misfolding: Neu- rodegenerative Diseases Edited by: Ovádi J, Orosz F. New York, USA: Springer; 2008:21-75.

Human disease genes. G Jimenez-Sanchez, B Childs, D Valle, Nature. 4096822Jimenez-Sanchez G, Childs B, Valle D: Human disease genes. Nature 2001, 409(6822):853-855.

Genomics and medicine. Dissecting human disease in the postgenomic era. L Peltonen, V A Mckusick, Science. 2915507Peltonen L, McKusick VA: Genomics and medicine. Dissecting human disease in the postgenomic era. Science 2001, 291(5507):1224-1229.

Towards a proteome-scale map of the human protein-protein interaction network. J F Rual, K Venkatesan, T Hao, T Hirozane-Kishikawa, A Dricot, N Li, G F Berriz, F D Gibbons, M Dreze, N Ayivi-Guedehoussou, Nature. 4377062Rual JF, Venkatesan K, Hao T, Hirozane-Kishikawa T, Dricot A, Li N, Berriz GF, Gibbons FD, Dreze M, Ayivi-Guedehoussou N, et al.: Towards a proteome-scale map of the human protein-pro- tein interaction network. Nature 2005, 437(7062):1173-1178.

A human protein-protein interaction network: a resource for annotating the proteome. U Stelzl, U Worm, M Lalowski, C Haenig, F H Brembeck, H Goehler, M Stroedicke, M Zenkner, A Schoenherr, S Koeppen, Cell. 1226Stelzl U, Worm U, Lalowski M, Haenig C, Brembeck FH, Goehler H, Stroedicke M, Zenkner M, Schoenherr A, Koeppen S, et al.: A human protein-protein interaction network: a resource for annotat- ing the proteome. Cell 2005, 122(6):957-968.

A protein-protein interaction network for human inherited ataxias and disorders of Purkinje cell degeneration. J Lim, T Hao, C Shaw, A J Patel, G Szabo, J F Rual, C J Fisk, N Li, A Smolyar, D E Hill, Cell. 1254Lim J, Hao T, Shaw C, Patel AJ, Szabo G, Rual JF, Fisk CJ, Li N, Smolyar A, Hill DE, et al.: A protein-protein interaction network for human inherited ataxias and disorders of Purkinje cell degeneration. Cell 2006, 125(4):801-814.

The human disease network. K I Goh, M E Cusick, D Valle, B Childs, M Vidal, A L Barabasi, Proc Natl Acad Sci. 10421Goh KI, Cusick ME, Valle D, Childs B, Vidal M, Barabasi AL: The human disease network. Proc Natl Acad Sci USA 2007, 104(21):8685-8690.

Lengthdependent prediction of protein intrinsic disorder. K Peng, P Radivojac, S Vucetic, A K Dunker, Z Obradovic, BMC Bioinformatics. 7208Peng K, Radivojac P, Vucetic S, Dunker AK, Obradovic Z: Length- dependent prediction of protein intrinsic disorder. BMC Bioin- formatics 2006, 7:208.

Intrinsic protein disorder in complete genomes. A K Dunker, Z Obradovic, P Romero, E C Garner, C J Brown, Genome Inform Ser Workshop Genome Inform. 11Dunker AK, Obradovic Z, Romero P, Garner EC, Brown CJ: Intrin- sic protein disorder in complete genomes. Genome Inform Ser Workshop Genome Inform 2000, 11:161-171.

On a test of whether one of two random variables is stochastically larger than the other. H B Mann, D R Whitney, Annals of Mathematical Statistics. 18Mann HB, Whitney DR: On a test of whether one of two ran- dom variables is stochastically larger than the other. Annals of Mathematical Statistics 1947, 18:50-60.

F Wilcoxon, Individual comparisons by ranking methods. Biometrics Bulletin. 1Wilcoxon F: Individual comparisons by ranking methods. Bio- metrics Bulletin 1945, 1:80-83.

Controlling the false discovery rate: a practical and powerful approach to multiple testing. Y Benjamini, Y , Journal of the Royal Statistical Society Series B (Methodological). 571Benjamini Y, Y H: Controlling the false discovery rate: a practi- cal and powerful approach to multiple testing. Journal of the Royal Statistical Society Series B (Methodological) 1995, 57(1):289-300.

The control of the false discovery rate in multiple testing under dependency. Y Benjamini, D Yekutieli, Annals of Statistics. 29Benjamini Y, Yekutieli D: The control of the false discovery rate in multiple testing under dependency. Annals of Statistics 2001, 29:1165-1188.

The problem of multiple comparisons. Princeton University. J W Tukey, unpublished reportTukey JW: The problem of multiple comparisons. Princeton Uni- versity; 1953. unpublished report.

Extension of multiple range tests to group means with unequal numbers of replications. C Y Kramer, Biometrics Bulletin. 12Kramer CY: Extension of multiple range tests to group means with unequal numbers of replications. Biometrics Bulletin 1956, 12:307-310.

Intrinsic disorder in pathogenic and non-pathogenic microbes: discovering and analyzing the unfoldomes of early-branching eukaryotes. A Mohan, W J SullivanJr, P Radivojac, A K Dunker, V N Uversky, Mol Biosyst. 44Mohan A, Sullivan WJ Jr, Radivojac P, Dunker AK, Uversky VN: Intrinsic disorder in pathogenic and non-pathogenic microbes: discovering and analyzing the unfoldomes of early-branching eukaryotes. Mol Biosyst 2008, 4(4):328-340.

Structural characterization of the native NH2-terminal transactivation domain of the human androgen receptor: a collapsed disordered conformation underlies structural plasticity and protein-induced folding. D N Lavery, I J Mcewan, Biochemistry. 4711Lavery DN, McEwan IJ: Structural characterization of the native NH2-terminal transactivation domain of the human androgen receptor: a collapsed disordered conformation underlies structural plasticity and protein-induced folding. Biochemistry 2008, 47(11):3360-3369.

Addressing the intrinsic disorder bottleneck in structural proteomics. C J Oldfield, E L Ulrich, Y Cheng, A K Dunker, J L Markley, Proteins. 593Oldfield CJ, Ulrich EL, Cheng Y, Dunker AK, Markley JL: Addressing the intrinsic disorder bottleneck in structural proteomics. Proteins 2005, 59(3):444-453.

Proteomics of Mycoplasma genitalium: identification and characterization of unannotated and atypical proteins in a small model genome. S Balasubramanian, T Schneider, M Gerstein, L Regan, Nucleic Acids Res. 2816Balasubramanian S, Schneider T, Gerstein M, Regan L: Proteomics of Mycoplasma genitalium: identification and characteriza- tion of unannotated and atypical proteins in a small model genome. Nucleic Acids Res 2000, 28(16):3075-3082.

Intrinsic disorder is a common feature of hub proteins from four eukaryotic interactomes. C Haynes, C J Oldfield, Ji F Klitgord, N Cusick, M E Radivojac, P Uversky, V N Vidal, M Iakoucheva, L M , PLoS Comput Biol. 28100Haynes C, Oldfield CJ, Ji F, Klitgord N, Cusick ME, Radivojac P, Uver- sky VN, Vidal M, Iakoucheva LM: Intrinsic disorder is a common feature of hub proteins from four eukaryotic interactomes. PLoS Comput Biol 2006, 2(8):e100.

Disordered domains and high surface charge confer hubs with the ability to interact with multiple proteins in interaction networks. A Patil, H Nakamura, FEBS Lett. 5808Patil A, Nakamura H: Disordered domains and high surface charge confer hubs with the ability to interact with multiple proteins in interaction networks. FEBS Lett 2006, 580(8):2041-2045.

Disorder and sequence repeats in hub proteins and their implications for network evolution. Z Dosztanyi, J Chen, A K Dunker, I Simon, P Tompa, J Proteome Res. 511Dosztanyi Z, Chen J, Dunker AK, Simon I, Tompa P: Disorder and sequence repeats in hub proteins and their implications for network evolution. J Proteome Res 2006, 5(11):2985-2995.

What properties characterize the hub proteins of the protein-protein interaction network of Saccharomyces cerevisiae?. D Ekman, S Light, A K Bjorklund, A Elofsson, Genome Biol. 7645Ekman D, Light S, Bjorklund AK, Elofsson A: What properties characterize the hub proteins of the protein-protein interac- tion network of Saccharomyces cerevisiae? Genome Biol 2006, 7(6):R45.

Intrinsic unstructuredness and abundance of PEST motifs in eukaryotic proteomes. G P Singh, M Ganapathi, K S Sandhu, D Dash, Proteins. 622Singh GP, Ganapathi M, Sandhu KS, Dash D: Intrinsic unstruc- turedness and abundance of PEST motifs in eukaryotic pro- teomes. Proteins 2006, 62(2):309-315.

M Goedert, Filamentous nerve cell inclusions in neurodegenerative diseases: tauopathies and alpha-synucleinopathies. Goedert M: Filamentous nerve cell inclusions in neurodegen- erative diseases: tauopathies and alpha-synucleinopathies.

. Philos Trans R Soc Lond B Biol Sci. 354Philos Trans R Soc Lond B Biol Sci 1999, 354(1386):1101-1118.

The alpha-synucleinopathies: Parkinson's disease, dementia with Lewy bodies, and multiple system atrophy. M G Spillantini, M Goedert, Ann N Y Acad Sci. 920Spillantini MG, Goedert M: The alpha-synucleinopathies: Parkin- son's disease, dementia with Lewy bodies, and multiple sys- tem atrophy. Ann N Y Acad Sci 2000, 920:16-27.

Parkinson's disease and related alphasynucleinopathies are brain amyloidoses. J Q Trojanowski, V M Lee, Ann N Y Acad Sci. 991Trojanowski JQ, Lee VM: Parkinson's disease and related alpha- synucleinopathies are brain amyloidoses. Ann N Y Acad Sci 2003, 991:107-110.

Synucleinopathies: clinical and pathological implications. J E Galvin, V M Lee, J Q Trojanowski, Arch Neurol. 582Galvin JE, Lee VM, Trojanowski JQ: Synucleinopathies: clinical and pathological implications. Arch Neurol 2001, 58(2):186-190.

Uversky VN: alpha-Synuclein misfolding and neurodegenerative diseases. Current Protein and Peptide Science. 95Uversky VN: alpha-Synuclein misfolding and neurodegenera- tive diseases. Current Protein and Peptide Science 2008, 9(5):507-540.

Uversky VN: Neuropathology, biochemistry, and biophysics of alpha-synuclein aggregation. J Neurochem. 1031Uversky VN: Neuropathology, biochemistry, and biophysics of alpha-synuclein aggregation. J Neurochem 2007, 103(1):17-37.

Ala30Pro mutation in the gene encoding alpha-synuclein in Parkinson's disease. R Kruger, W Kuhn, T Muller, D Woitalla, M Graeber, S Kosel, H Przuntek, J T Epplen, L Schols, O Riess, Nat Genet. 182Kruger R, Kuhn W, Muller T, Woitalla D, Graeber M, Kosel S, Przuntek H, Epplen JT, Schols L, Riess O: Ala30Pro mutation in the gene encoding alpha-synuclein in Parkinson's disease. Nat Genet 1998, 18(2):106-108.

Mutation in the alphasynuclein gene identified in families with Parkinson's disease. M H Polymeropoulos, C Lavedan, E Leroy, S E Ide, A Dehejia, A Dutra, B Pike, H Root, J Rubenstein, R Boyer, Science. 2765321Polymeropoulos MH, Lavedan C, Leroy E, Ide SE, Dehejia A, Dutra A, Pike B, Root H, Rubenstein J, Boyer R, et al.: Mutation in the alpha- synuclein gene identified in families with Parkinson's disease. Science 1997, 276(5321):2045-2047.

The new mutation, E46K, of alpha-synuclein causes Parkinson and Lewy body dementia. J J Zarranz, J Alegre, J C Gomez-Esteban, E Lezcano, Ros R Ampuero, I Vidal, L Hoenicka, J Rodriguez, O Atares, B , Ann Neurol. 552Zarranz JJ, Alegre J, Gomez-Esteban JC, Lezcano E, Ros R, Ampuero I, Vidal L, Hoenicka J, Rodriguez O, Atares B, et al.: The new muta- tion, E46K, of alpha-synuclein causes Parkinson and Lewy body dementia. Ann Neurol 2004, 55(2):164-173.

Comparison of kindreds with parkinsonism and alpha-synuclein genomic multiplications. M Farrer, J Kachergus, L Forno, S Lincoln, D S Wang, M Hulihan, D Maraganore, K Gwinn-Hardy, Z Wszolek, D Dickson, Ann Neurol. 552Farrer M, Kachergus J, Forno L, Lincoln S, Wang DS, Hulihan M, Mara- ganore D, Gwinn-Hardy K, Wszolek Z, Dickson D, et al.: Compari- son of kindreds with parkinsonism and alpha-synuclein genomic multiplications. Ann Neurol 2004, 55(2):174-179.

A B Singleton, M Farrer, J Johnson, A Singleton, S Hague, J Kachergus, M Hulihan, T Peuralinna, A Dutra, R Nussbaum, alpha-Synuclein locus triplication causes Parkinson's disease. 302841Singleton AB, Farrer M, Johnson J, Singleton A, Hague S, Kachergus J, Hulihan M, Peuralinna T, Dutra A, Nussbaum R, et al.: alpha-Synu- clein locus triplication causes Parkinson's disease. Science 2003, 302(5646):841.

Association between cardiac denervation and parkinsonism caused by alpha-synuclein gene triplication. A Singleton, K Gwinn-Hardy, Y Sharabi, S T Li, C Holmes, R Dendi, J Hardy, A Crawley, D S Goldstein, Brain. 127Pt 4Singleton A, Gwinn-Hardy K, Sharabi Y, Li ST, Holmes C, Dendi R, Hardy J, Crawley A, Goldstein DS: Association between cardiac denervation and parkinsonism caused by alpha-synuclein gene triplication. Brain 2004, 127(Pt 4):768-772.

Alpha-synuclein in blood and brain from familial Parkinson disease with SNCA locus triplication. D W Miller, S M Hague, J Clarimon, M Baptista, K Gwinn-Hardy, M R Cookson, A B Singleton, Neurology. 6210Miller DW, Hague SM, Clarimon J, Baptista M, Gwinn-Hardy K, Cook- son MR, Singleton AB: Alpha-synuclein in blood and brain from familial Parkinson disease with SNCA locus triplication. Neu- rology 2004, 62(10):1835-1838.

Part II: alpha-synuclein and its molecular pathophysiological role in neurodegenerative disease. K K Dev, K Hofele, S Barbieri, V L Buchman, H Putten, Van Der, Neuropharmacology. 451Dev KK, Hofele K, Barbieri S, Buchman VL, Putten H van der: Part II: alpha-synuclein and its molecular pathophysiological role in neurodegenerative disease. Neuropharmacology 2003, 45(1):14-44.

Alpha-synuclein and the Lewy body disorders. D W Dickson, Curr Opin Neurol. 144Dickson DW: Alpha-synuclein and the Lewy body disorders. Curr Opin Neurol 2001, 14(4):423-432.

Alpha-synuclein and neurodegenerative diseases. M Goedert, Nat Rev Neurosci. 27Goedert M: Alpha-synuclein and neurodegenerative diseases. Nat Rev Neurosci 2001, 2(7):492-501.

Parkinson's disease and other alpha-synucleinopathies. M Goedert, Clin Chem Lab Med. 394Goedert M: Parkinson's disease and other alpha-synucleinopa- thies. Clin Chem Lab Med 2001, 39(4):308-312.

A protein-chameleon: conformational plasticity of alpha-synuclein, a disordered protein involved in neurodegenerative disorders. V N Uversky, J Biomol Struct Dyn. 212Uversky VN: A protein-chameleon: conformational plasticity of alpha-synuclein, a disordered protein involved in neurode- generative disorders. J Biomol Struct Dyn 2003, 21(2):211-234.

J Jin, G J Li, J Davis, D Zhu, Y Wang, C Pan, J Zhang, Identification of novel proteins interacting with both a-synuclein and DJ-1. Jin J, Li GJ, Davis J, Zhu D, Wang Y, Pan C, Zhang J: Identification of novel proteins interacting with both a-synuclein and DJ-1.

. Mol Cell Proteomics. Mol Cell Proteomics 2006.

Glutamine repeats and inherited neurodegenerative diseases: molecular aspects. M F Perutz, Curr Opin Struct Biol. 66Perutz MF: Glutamine repeats and inherited neurodegenera- tive diseases: molecular aspects. Curr Opin Struct Biol 1996, 6(6):848-858.

Polyglutamine diseases: protein cleavage and aggregation. H Y Zoghbi, H T Orr, Curr Opin Neurobiol. 95Zoghbi HY, Orr HT: Polyglutamine diseases: protein cleavage and aggregation. Curr Opin Neurobiol 1999, 9(5):566-570.

Polyglutamine pathogenesis. C A Ross, J D Wood, G Schilling, M F Peters, F C NuciforaJr, J K Cooper, A H Sharp, R L Margolis, D R Borchelt, Philos Trans R Soc Lond B Biol Sci. 354Ross CA, Wood JD, Schilling G, Peters MF, Nucifora FC Jr, Cooper JK, Sharp AH, Margolis RL, Borchelt DR: Polyglutamine patho- genesis. Philos Trans R Soc Lond B Biol Sci 1999, 354(1386):1005-1011.

Evidence for a recruitment and sequestration mechanism in Huntington's disease. E Preisinger, B M Jordan, A Kazantsev, D Housman, Philos Trans R Soc Lond B Biol Sci. 354Preisinger E, Jordan BM, Kazantsev A, Housman D: Evidence for a recruitment and sequestration mechanism in Huntington's disease. Philos Trans R Soc Lond B Biol Sci 1999, 354(1386):1029-1034.

Protein aggregation and pathogenesis of Huntington's disease: mechanisms and correlations. E E Wanker, Biol Chem. 3819Wanker EE: Protein aggregation and pathogenesis of Hunting- ton's disease: mechanisms and correlations. Biol Chem 2000, 381(9-10):937-942.

Aggregation of proteins with expanded glutamine and alanine repeats of the glutamine-rich and asparagine-rich domains of Sup35 and of the amyloid beta-peptide of amyloid plaques. M F Perutz, B J Pope, D Owen, E E Wanker, E Scherzinger, Proc Natl Acad Sci. 998Perutz MF, Pope BJ, Owen D, Wanker EE, Scherzinger E: Aggrega- tion of proteins with expanded glutamine and alanine repeats of the glutamine-rich and asparagine-rich domains of Sup35 and of the amyloid beta-peptide of amyloid plaques. Proc Natl Acad Sci USA 2002, 99(8):5596-5600.

CREB-binding protein sequestration by expanded polyglutamine. A Mccampbell, J P Taylor, A A Taye, J Robitschek, M Li, J Walcott, D Merry, Y Chai, H Paulson, G Sobue, Hum Mol Genet. 914McCampbell A, Taylor JP, Taye AA, Robitschek J, Li M, Walcott J, Merry D, Chai Y, Paulson H, Sobue G, et al.: CREB-binding protein sequestration by expanded polyglutamine. Hum Mol Genet 2000, 9(14):2197-2202.

Polyglutamine and CBP: fatal attraction. A Mccampbell, K H Fischbeck, Nat Med. 75McCampbell A, Fischbeck KH: Polyglutamine and CBP: fatal attraction? Nat Med 2001, 7(5):528-530.

Polyglutamine aggregation behavior in vitro supports a recruitment mechanism of cytotoxicity. S Chen, V Berthelier, W Yang, R Wetzel, J Mol Biol. 3111Chen S, Berthelier V, Yang W, Wetzel R: Polyglutamine aggrega- tion behavior in vitro supports a recruitment mechanism of cytotoxicity. J Mol Biol 2001, 311(1):173-182.

Amyloidlike features of polyglutamine aggregates and their assembly kinetics. S Chen, V Berthelier, J B Hamilton, B O'nuallain, R Wetzel, Biochemistry. 4123Chen S, Berthelier V, Hamilton JB, O'Nuallain B, Wetzel R: Amyloid- like features of polyglutamine aggregates and their assembly kinetics. Biochemistry 2002, 41(23):7391-7399.

Huntingtin interacts with a family of WW domain proteins. P W Faber, G T Barnes, J Srinidhi, J Chen, J F Gusella, M E Macdonald, Hum Mol Genet. 79Faber PW, Barnes GT, Srinidhi J, Chen J, Gusella JF, MacDonald ME: Huntingtin interacts with a family of WW domain proteins. Hum Mol Genet 1998, 7(9):1463-1474.

Aberrant interactions of transcriptional repressor proteins with the Huntington's disease gene product, huntingtin. J M Boutell, P Thomas, J W Neal, V J Weston, J Duce, P S Harper, A L Jones, Hum Mol Genet. 89Boutell JM, Thomas P, Neal JW, Weston VJ, Duce J, Harper PS, Jones AL: Aberrant interactions of transcriptional repressor pro- teins with the Huntington's disease gene product, hunting- tin. Hum Mol Genet 1999, 8(9):1647-1655.

The Huntington's disease protein interacts with p53 and CREBbinding protein and represses transcription. J S Steffan, A Kazantsev, O Spasic-Boskovic, M Greenwald, Y Z Zhu, H Gohler, E E Wanker, G P Bates, D E Housman, L M Thompson, Proc Natl Acad Sci. 9712Steffan JS, Kazantsev A, Spasic-Boskovic O, Greenwald M, Zhu YZ, Gohler H, Wanker EE, Bates GP, Housman DE, Thompson LM: The Huntington's disease protein interacts with p53 and CREB- binding protein and represses transcription. Proc Natl Acad Sci USA 2000, 97(12):6763-6768.

Molecular cloning and characterization of a cellular phosphoprotein that interacts with a conserved C-terminal domain of adenovirus E1A involved in negative modulation of oncogenic transformation. U Schaeper, J M Boyd, S Verma, E Uhlmann, T Subramanian, G Chinnadurai, Proc Natl Acad Sci. 9223Schaeper U, Boyd JM, Verma S, Uhlmann E, Subramanian T, Chinna- durai G: Molecular cloning and characterization of a cellular phosphoprotein that interacts with a conserved C-terminal domain of adenovirus E1A involved in negative modulation of oncogenic transformation. Proc Natl Acad Sci USA 1995, 92(23):10467-10471.

Live or let die: the cell's response to p53. K H Vousden, X Lu, Nat Rev Cancer. 28Vousden KH, Lu X: Live or let die: the cell's response to p53. Nat Rev Cancer 2002, 2(8):594-604.

Signaling to the p53 tumor suppressor through pathways activated by genotoxic and nongenotoxic stress. C W Anderson, E Appella, R A Bradshaw, E A Dennis, Handbook of Cell Signaling. New YorkAcademic PressAnderson CW, Appella E: Signaling to the p53 tumor suppres- sor through pathways activated by genotoxic and nongeno- toxic stress. In Handbook of Cell Signaling Edited by: Bradshaw RA, Dennis EA. New York: Academic Press; 2003:237-247.

p53 mutations in human cancers. M Hollstein, D Sidransky, B Vogelstein, C C Harris, Science. 2535015Hollstein M, Sidransky D, Vogelstein B, Harris CC: p53 mutations in human cancers. Science 1991, 253(5015):49-53.

Surfing the p53 network. B Vogelstein, D Lane, A J Levine, Nature. 4086810Vogelstein B, Lane D, Levine AJ: Surfing the p53 network. Nature 2000, 408(6810):307-310.

The p53 mutation "gradient effect" and its clinical implications. G P Zambetti, J Cell Physiol. 2132Zambetti GP: The p53 mutation "gradient effect" and its clin- ical implications. J Cell Physiol 2007, 213(2):370-373.

p53 and human cancer: the first ten thousand mutations. P Hainaut, M Hollstein, Adv Cancer Res. 77Hainaut P, Hollstein M: p53 and human cancer: the first ten thousand mutations. Adv Cancer Res 2000, 77:81-137.

The N-terminal domain of p53 is natively unfolded. R Dawson, L Muller, A Dehner, C Klein, H Kessler, J Buchner, J Mol Biol. 3325Dawson R, Muller L, Dehner A, Klein C, Kessler H, Buchner J: The N-terminal domain of p53 is natively unfolded. J Mol Biol 2003, 332(5):1131-1141.

Local structural elements in the mostly unstructured transcriptional activation domain of human p53. H Lee, K H Mok, R Muhandiram, K H Park, J E Suk, D H Kim, J Chang, Y C Sung, K Y Choi, K H Han, J Biol Chem. 27538Lee H, Mok KH, Muhandiram R, Park KH, Suk JE, Kim DH, Chang J, Sung YC, Choi KY, Han KH: Local structural elements in the mostly unstructured transcriptional activation domain of human p53. J Biol Chem 2000, 275(38):29426-29432.

Analysis of p53-regulated gene expression patterns using oligonucleotide arrays. R Zhao, K Gish, M Murphy, Y Yin, D Notterman, W H Hoffman, E Tom, D H Mack, A J Levine, Genes Dev. 148Zhao R, Gish K, Murphy M, Yin Y, Notterman D, Hoffman WH, Tom E, Mack DH, Levine AJ: Analysis of p53-regulated gene expres- sion patterns using oligonucleotide arrays. Genes Dev 2000, 14(8):981-993.

Inherited breast and ovarian cancer. C I Szabo, M C King, Hum Mol Genet. 4Szabo CI, King MC: Inherited breast and ovarian cancer. Hum Mol Genet 1995, 4(Spec No):1811-1817.

Founder mutations in BRCA1 and BRCA2 genes. R Ferla, V Calo, S Cascio, G Rinaldi, G Badalamenti, I Carreca, E Surmacz, G Colucci, V Bazan, A Russo, Ann Oncol. 186SupplFerla R, Calo V, Cascio S, Rinaldi G, Badalamenti G, Carreca I, Sur- macz E, Colucci G, Bazan V, Russo A: Founder mutations in BRCA1 and BRCA2 genes. Ann Oncol 2007, 18(Suppl 6):vi93-98.

Roles of BRCA1 and its interacting proteins. C X Deng, S G Brodie, Bioessays. 228Deng CX, Brodie SG: Roles of BRCA1 and its interacting pro- teins. Bioessays 2000, 22(8):728-737.

Characterization of segments from the central region of BRCA1: an intrinsically disordered scaffold Submit your manuscript here. W Y Mark, J C Liao, Y Lu, A Ayed, R Laister, B Szymczyna, A Chakrabartty, C H Arrowsmith, BioMedcentral BMC Genomics. 101SupplMark WY, Liao JC, Lu Y, Ayed A, Laister R, Szymczyna B, Chakrabar- tty A, Arrowsmith CH: Characterization of segments from the central region of BRCA1: an intrinsically disordered scaffold Submit your manuscript here: http://www.biomedcentral.com/info/publishing_adv.asp BioMedcentral BMC Genomics 2009, 10(Suppl 1):S12 http://www.biomedcentral.com/1471-2164/10/S1/S12

The protein trinity -linking function and disorder. A K Dunker, Z Obradovic, Nat Biotechnol. 199Dunker AK, Obradovic Z: The protein trinity -linking function and disorder. Nat Biotechnol 2001, 19(9):805-806.

Coupling of folding and binding for unstructured proteins. H J Dyson, P E Wright, Curr Opin Struct Biol. 121Dyson HJ, Wright PE: Coupling of folding and binding for unstructured proteins. Curr Opin Struct Biol 2002, 12(1):54-60.

Intrinsically unstructured proteins and their functions. H J Dyson, P E Wright, Nat Rev Mol Cell Biol. 63Dyson HJ, Wright PE: Intrinsically unstructured proteins and their functions. Nat Rev Mol Cell Biol 2005, 6(3):197-208.

Homooligomerization of the cytoplasmic domain of the T cell receptor zeta chain and of other proteins containing the immunoreceptor tyrosine-based activation motif. A Sigalov, D Aivazian, L Stern, Biochemistry. 437Sigalov A, Aivazian D, Stern L: Homooligomerization of the cyto- plasmic domain of the T cell receptor zeta chain and of other proteins containing the immunoreceptor tyrosine-based activation motif. Biochemistry 2004, 43(7):2049-2061.

Multichain immune recognition receptor signaling: different players, same game?. A B Sigalov, Trends Immunol. 2511Sigalov AB: Multichain immune recognition receptor signal- ing: different players, same game? Trends Immunol 2004, 25(11):583-589.

Immune cell signaling: a novel mechanistic model reveals new therapeutic targets. A B Sigalov, Trends Pharmacol Sci. 2710Sigalov AB: Immune cell signaling: a novel mechanistic model reveals new therapeutic targets. Trends Pharmacol Sci 2006, 27(10):518-524.

Lipid-binding activity of intrinsically unstructured cytoplasmic domains of multichain immune recognition receptor signaling subunits. A B Sigalov, D A Aivazian, V N Uversky, L J Stern, Biochemistry. 4551Sigalov AB, Aivazian DA, Uversky VN, Stern LJ: Lipid-binding activ- ity of intrinsically unstructured cytoplasmic domains of mul- tichain immune recognition receptor signaling subunits. Biochemistry 2006, 45(51):15731-15739.

Binding of intrinsically disordered proteins is not necessarily accompanied by a structural transition to a folded form. A B Sigalov, A V Zhuravleva, V Y Orekhov, Biochimie. 893Sigalov AB, Zhuravleva AV, Orekhov VY: Binding of intrinsically disordered proteins is not necessarily accompanied by a structural transition to a folded form. Biochimie 2007, 89(3):419-421.

Predicting Binding Regions within Disordered Proteins. E Garner, P Romero, A K Dunker, C Brown, Z Obradovic, Genome Inform Ser Workshop Genome Inform. 10Garner E, Romero P, Dunker AK, Brown C, Obradovic Z: Predict- ing Binding Regions within Disordered Proteins. Genome Inform Ser Workshop Genome Inform 1999, 10:41-50.