intrinsically disordered proteins ( idps ) and regions ( idrs ) have biological activities that , at least for part of the time , require the absence of stable 3-dimensional or secondary structure under physiological conditions . numbers estimating the amount of intrinsic disorder in proteins are stunning ; about 43% of known mammalian protein sequence is in predicted intrinsic disorder ( pid ) . between 35 and 51% of eukaryotic proteins have been predicted to contain idrs that span 40 or more residues . between 25 and 30% of eukaryotic proteins more than 70% of signaling proteins , and most of the cancer - associated proteins have been predicted to contain long disordered regions . sequence / structure information in disprot has been used by members of the idp community to build more than 50 different methods for predicting regions of intrinsic disorder in proteins , and to estimate statistics for the accuracy of these methods . nine of these predictors , along with predictions for 10,429,761 sequences in 1,765 proteomes from 1,256 distinct species , are available on the dp site . direct experimental evidence for the existence of idps and idrs comes primarily from the protein data bank ( pdb ) where nmr solution structures show conformational ensembles that clearly indicate dynamic disorder . many pdb entries contain segments of protein sequence that are completely missing from x - ray and neutron diffraction crystal structures , but as is well - known by crystallographers , these segments can correspond to structured regions that are unobserved for a variety of technical reasons . for this reason , intrinsic disorder can not be assigned to residues in crystallized proteins solely because they are not located . nevertheless , 7% of the crystal structures in the pdb , which is highly selective for ordered proteins , have been assigned to idrs longer than 10 residues . idps and proteins with significant idrs continue to resist crystallization , a situation that is unlikely to change . with few exceptions , crystal structures of largely disordered proteins have been obtained only from relatively small isolated sections of biologically active idps , many of them co - crystallized with , or covalently bound to , much larger structured molecular complexes . even under these conditions many residues in the idrs can not be located in the electron density , and the crystallized complex may not represent the interaction in solution . solution nmr data from a variety of chemical shift , relaxation , and heteronuclear noe measurements yield unique information about the spectrum of conformational disorder and dynamics in proteins that is less specific but more accurate than that obtained from diffraction measurements . taken together , and supported by circular dichroism , vibrational spectroscopy , chromatography , and small angle scattering methods , these measurements provide certainty about the extent of static and dynamic disorder in idps and idrs , and this information is recorded in the disprot database . there are many biophysical techniques that can be used to characterize dynamic structure of idps , and many of these methods have been the subjects of focused reviews and books . the p53 and 143 - 3 proteins provide 2 particularly striking examples of the biological activity of intrinsically disordered proteins . the p53 protein has several different idrs that bind to different partners , and some idrs that each bind in different conformations to several different partners , an association termed  one - to - many .  14 - 3 - 3 on the other hand is a structured protein that binds many different intrinsically disordered partners in associations termed  many - to - one  . the pfam - a database is a curated collection of biologically conserved , and for many  functional , regions in proteins . pfam sequences are grouped in large part by function and used to train hidden markov models ( hmms ) that are used to find similar regions in proteins where there is no protein based evidence of biological activity . the training set of sequences , pfam-a.seed , contains regions from proteins that have been experimentally validated . it is this set that has been used in the present study to relate biological function with intrinsic disorder . recent work , looked at 71,974 version 22.0 pfam - a seed members of 6,857 unique domains , limited to those that included go annotations or had at least one literature citation . 12.14% of the domains had greater than 50% predicted disorder , and 4.15% were fully ( 95100% ) disordered . the high percentage of fully disordered domains was attributed to the uneven length distribution of domains , with somewhat shorter domains dominating at high percentage of disorder . earlier , 40% of pfam domains were shown to contain conserved protein fragments that were predicted to be disordered ( conserved disorder predictions , cdps ) . these cdps were found in proteins from all domains / kingdoms of life , including viruses , with eukaryota having one order of magnitude more proteins containing long disordered regions than did archaea and bacteria . functional analyses revealed that cdp regions frequently participate in signaling , regulation , and interaction with dna / rna and other proteins , common in ribosomal proteins . in the present work here we analyze the distribution of intrinsic disorder in pfam domain sequences using the 10 dimensional space provided by the kidera factors for the 20 naturally occurring amino acids combined with pondr predictions of intrinsic disorder . the kidera factors have been developed expressly to describe properties of the amino acid residues that may be related to protein folding with a minimum number of parameters . they are derived from a multivariate statistical analysis beginning with 188 quantitative measurements of the amino acids available in 1985 . because this set of factors contains most of the measurable information relating the amino acids , it is possible to estimate the numerical values for amino acids where measurements may be missing . four of the kidera factors , helix / bend preference , side - chain size , extended structure preference , and hydrophobicity are essentially pure factors ( table  1 ) . for example , the cluster for hydrophobicity contains only the relatively large set of measurements related to amino acid solubility . these factors are principle components of a large set of experimental measurements , scaled from 1 to 1.3 the first 4 : hel , siz , ext , and hph , are close to being pure physical properties . kidera factors , being orthonormal , contain 10 times the information contained in the sequence alone . an analysis of short range interactions in sequences benefits when these factors are included , potentially increasing the information available to an hmm analysis many fold . the remaining 6 factors consist of weighted linear combinations of different measurements , labeled for convenience by the name of the most heavily weighted component . where factors appear to have names related to similar properties , the similarity is in name only . for example , the vectors composed of the 10 factors for extended structure preference ( ext ) and flat extended preference ( fle ) for each of the 20 amino acids are themselves orthonormal . likewise , the similarly named pairs hel / alp and siz / psb are also orthonormal ( table  1 ) ; there is no correlation between these factors . as indicated by the authors , the kidera factors are orthonormal by design to avoid problems arising from incompleteness and correlation ( they are normalized and their inner product is zero . ) they do not contain information about interactions as may arise from an analysis of a length of sequence for periodicity or interactions . information about interactions can be derived from a sliding window analysis of any of the 10 factors . as we show here , an average over sets of orthonormal factors can yield correlated results with reduced information content . some of the factors included in the original kidera data set have been included in the set used to train the pondr predictors used to make predictions here . these include the hydrophobicity scales from kyte and doolittle and from rose ( for predicting turns in globular proteins ) , and side chain volume . also included in pondr training were charge ( k+rd - e ) , aromatic count ( w + f + y ) , and coordination number , which correlate strongly with hydrophobicity and side chain volume . side chain volume , related to the convenience named partial specific volume kidera factor , does not come close to being a pure physical property ( table  1 and ref . other factors included in pondr training depend on the kidera information in lengths of sequence and include a flexibility index calculated from a sliding window , the hydrophobic periodicity moments from eisenberg , weiss , and terwilliger , codon number , and alphabet size . these scales are derived from a small part of the information contained in the kidera data . the pondr predictors continue to be among the most accurate available , suggesting that hydrophobicity , may be among the principle physical properties of the amino acids determining the tendency of proteins to evolve with functions in intrinsically disordered states . however , we show here that all 10 of the kidera factors contribute to a clustering of different types of predicted intrinsic disorder in pfam domains , and that plots of pondr vsl2b predictor scores against hydrophobicity have a spread and appearance that is similar to that seen in plots of the vsl2b scores against any of the other kidera factors . the inclusion of the kidera factors in the training of disorder predictors could hypothetically increase the information content available to an analysis of intrinsic protein disorder by several fold . that while figure 1 includes the 1520% of pfam members that contained no pid , subsequent figures exclude members with 0% pid by simply excluding data below 2% from the plot . end effects do not contribute significantly to disorder distributions , as is shown here and in table  2 in a comparison of the distribution of percent predicted disorder , as a function of the percentage of predicted disorder in 2% wide bins , in all pfam seed proteins where domain members start or end within 19 , 29 , and 39 residues of the whole seed protein ends . all predictions in this work were performed on whole proteins , not on isolated domains , so sequence end prediction artifacts are restricted to pfam domains at the n or c - terminus of proteins . when all proteins where end effects may affect prediction are removed , the prevalence of 100% predicted disordered domain members , shown here at 100% and in previous work , does not significantly decrease ( table  2 ) . the inset shows the distribution of lengths for 100% predicted intrinsically disordered ( pid ) pfam sequences . the mean length for pid regions in whole mammalian proteins in the pfam seed set is 16 residues ( table  7 ) . testing the hypothesis that intrinsic protein disorder preference , like helix preference , is a fundamental physical property of the amino acids that was not included when the kidera factors were calculated , a singular value decomposition of the 20 by 11 matrix composed of the kidera factors and the frequencies of occurrence of the amino acids in idps yielded only 10 non - zero eigenvalues ; the disorder preferences of the amino acids are linear combinations of the kidera factors . of the 618 100% pid members of pfam domains , the vsl2b predictor is biased to assign disordered structure to the first 20 residues of a protein or isolated segment . we tested the hypothesis that about 176 members may be incorrectly predicted to be 100% disordered and that this error biased our calculations . the results , shown in figure 1 and table  2 , indicate that end effects do not contribute significantly to the 100% pid group . we test the null hypothesis that the distribution for ( 1 ) the object under percent disorder statistics are given as the quartiles of the distribution of pid calculated with r from a table where each row lists a domain or protein followed by the % disorder in that domain or protein . compared with : this test , and the 99.9% confidence interval , are calculated using the r wilcoxon rank sum ( mann - whitney ) test with continuity correction . the mann - whitney test is appropriate and accurate for comparing the medians of 2 large sample non - paired non - normal distributions , when those distributions are the same . mammals include only human , mouse , rat , bovine , rabbit , pig , and horse . the sequence length mean of 100% pid pfam sequences is 82 , the median is just below 50 , and includes many above 150 . this broad and skewed distribution distorts the cluster analysis performed here to some extent where averages of kidera factors taken over long sequences tend toward their central value of zero , while in shorter sequences there is a greater chance that a particular type of intrinsic disorder may be isolated , in a way that is analogous to searching for segments of helix , -strand , or turn in folded structures . earlier work finding 3 distinct albeit overlapping flavors of disorder evaluated amino acid composition in windows 41 residues long . the effects of our use of a homogeneous distribution of lengths here are discussed below where we evaluate the high dimensional analysis of our data . figure  2 and table  3 show the distributions of pid in several proteomes compared with that in human pfam seed proteins . there are more long pid regions in the proteome than in the pfam seeds for humans . about 4% of the chimp , mouse , and human proteome whole proteins are 100% pid . lower phylogenetic domains have smaller proportions of 100% disordered pfam domains , and of disorder overall , in the following order : mammals > other eukaryota > viruses > bacteria > archaea . distributions of disorder in each of the 5 phylogenetic domains have positions that differ at the 0.999 level with p values less than 1  10 figure  2 . in whole proteins there appears to be a predominance of 100% pid over what is found in the pfam seed database for humans , with more ( here ) in mouse and chimp than in human proteins . this can not be attributed to an uneven distribution of pfam domain or pid lengths , and suggests that the pfam seed database excludes some idps . this is shown here in comparisons of the distribution of percent predicted disorder in the human proteome with the pantroglodytes chimpanzee and musmusculus mouse proteomes , and with human proteins chosen as sources for pfam seed members ( red ) . we note that the distribution for mouse is shifted 1% to the left of that for chimp and human proteomes while 100% pid is highest for the mouse . these figures derive from whole proteins from entire proteomes , and not pfam domains . overall , 28% of pfam members with 100% pid are derived from whole proteins that are 100% pid , and mammalian domains in the pfam seed database contain 27% pid . this is 13.8% less than in their source proteins and 16% less than in mammalian proteomes . this can be seen , with the additional 3.1% , in tables  3 and 4 . we hypothesize that this is an artifact of low complexity sequence filtering in the selection of pfam seeds , by accident or design . the null hypothesis , that each distribution does not differ from any other , can be rejected at the 0.999 level . percent disorder statistics are given here as the quartiles of the distribution of intrinsic disorder calculated with r from a table where each row lists a domain or protein followed by the % disorder in that domain or protein . mammals included only human , mouse , rat , bovine , rabbit , pig , and horse . eukaryota here have only these removed.99.9% confidence interval calculated using the r wilcoxon rank sum test with continuity correction . the null hypothesis is that the distributions for mammals and the other domains listed here differ by a location shift of zero . the alternative is that they differ by some other one or two sided location shift . figure  3 , and tables  4 and 5 , show the distributions of pid for pfam seed members ( regions of proteins ) and the whole proteins from which seeds are derived . figure  3 . left : a comparison of the distribution of percent pid in pfam seed domain members and pfam whole proteins for mammals and other eukaryota in 2% wide bins show 2 things . about 1.7% of whole proteins and 6% of pfam members are predicted to be 100% pid , but whole proteins contain significantly more pid overall . again , whole proteins contain 28% of the pfam sequences predicted to be 100% disordered , and do not fall into the category where uneven length distribution of domains accounts for 100% disordered pfam domains . also , predicted disorder is estimated to be 13% lower in pfam domains than in pfam whole proteins , shown in table  4 ,  est dif . right : here , the median length of predicted 100% disordered whole proteins is about 70 residues longer than that of predicted 100% disordered pfam domains ( table  5 ) , and in figure 4 the median length of pfam domains is much larger than the median length of predicted intrinsically disordered regions . clearly , some 100% pid pfam sequences derive from whole proteins where pid extends beyond the ends of the pfam segments as proposed earlier , but there is no obvious reason why this classifies the significant category of predicted entirely disordered pfam sequences as an artifact .   as in table  4 . figure  4 , and tables  6 and 7 , show the distributions for percent disorder in pfam members , pfam member length , predicted intrinsic disorder ( i d ) length in whole proteins , and numbers of i d regions in whole proteins . distributions of pfam seed sequences for mammals , other eukaryota , viruses , bacteria , and archaea , showing , upper left : the percent of pfam domains as a function of the percentage of predicted disorder in 2% wide bins , and the spike in the number of domains that are 100% disordered ( see also table  4 ) . upper right : pfam domain sequence lengths for mammals are 12 to 40 residues shorter than they are in other eukaryota , viruses , bacteria , and archaea , but the order of domain sequence lengths does not follow the order of 100% pid , as also shown in table  7 under  shift .  the mean domain length is 145 residues ( table  7 ) . lower left : intrinsically disordered region ( idr ) lengths ( distinct from pfam domain lengths ) for mammals are 15 residues longer than in other domains / kingdoms of life , also shown in table  6 , and the order of i d length follows the order of 100% pid . most predicted intrinsically disordered regions are shorter than 10 residues , much shorter than the median pfam domain sequence length . lower right : there are significantly fewer idrs in mammals than in viruses , and marginally fewer than in other eukaryota . the median is near 2 , and some proteins are predicted to have more than 20 idrs .   as in table  4 . shown here for comparison , this line includes allmmammalian predicted intrinsically disordered region lengths , and not pfam domains ; see figure 4 lower left . most striking is the spike at the right side of figure 4 top left indicating that there are more mammalian pfam members 100% pid than there are in any other 2% wide bin of the data . there are 618 100% pid mammalian pfam members in this set . our initial hypothesis here  that this set of pid sequences was characterized by factors that differed from those in other pid sequences  could not be supported . also shown in figure 4 top right and bottom left , while most pid sequences are shorter than 10 residues , very few pfam members are this short , and the mean pfam member length is 145 residues ( table  7 ) . table 7 shows that each 100% pid pfam member clearly derives from a pid region that is longer than the pfam member . however , there is no compelling evidence here that the presence of entirely disordered pfam members is an artifact , and does not have a special evolutionary significance , perhaps conferring an advantage to  higher  or  lower  phylogenetic domains . the evidence against artifact is also not especially compelling : the mean length of pfam domain members , 145 residues , is much longer here than the mean length of pid regions , 16 residues . pfam member lengths for mammals are 1240 residues shorter than they are in other eukaryota , viruses , bacteria , and archaea . there appears to be a predominance of 100% pid in the eukaryotic proteome that is also seen in pfam domain members . predicted disorder is estimated here to be 13% lower in pfam domains than in pfam whole proteins . twenty - eight percent of pfam members that are 100% pid are derived from whole proteins that are 100% pid . figure  5 shows ch / cdf scatter plots for the distributions of predicted disorder in pfam - a version 23.0 members in mammals ( 10,660 pfam members ) , other eukaryota ( 71,765 ) , viruses ( 6,360 ) , bacteria ( 101,959 ) , and archaea ( 12,721 ) . the ch / cdf plot shows , for all practical purposes , the vsl2b disorder prediction score on the x axis and hydropathy on the y axis . this can be seen in the scatter plot matrix in the lower right where ch , hydropathy , cdf , and the vsl2b scores for mammals are all plotted against each other . ch / cdf plots for mammals , other eukaryota , viruses , bacteria , and archaea . the scatterplot matrix ( bottom right ) shows density ( number of sequences ) as a function of ch , hydrophobicity , cdf , and vsl2b , while the off - diagonal plots show the correlations between each of these parameters . the off - diagonal plots : ( 1 ) corresponding to ch / cdf on the diagonal duplicate the plot for all pfam seed sequences for mammals ( top left ) , and ( 2 ) corresponding to ch or cdf on the one hand and hydrophobicity or vsl2b on the other , show the high correlation between ch and hydrophobicity , and between cdf and vsl2b . only 2 or 3 clearly separate clusters are evident here , but the cluster analyses of the mammalian data , shown in figures 915 reveal much more information . density in the lower and middle right hand sections of the ch / cdf plots corresponds to mostly ordered pfam members , while density in the middle and upper left hand sections represents mostly disordered pfam members . as can be seen , density shifts from the lower right to the upper left as plots progress from archaea , bacteria , viruses , and other eukaryota , to mammals . the grouping in the lower right hand side of the plots for archaea and bacteria in figure 5 , corresponding to proteins with no predicted intrinsic disorder , is not included in the cluster analyses below . figure  6 shows more detailed structure than can be seen in the figure 5 . each plot here contrasts data from all of the phylogenetic domains , and shows the differences between domains more clearly than can be seen in figure 4 . here the percent of pfam members is on the y axis , and the percent of predicted disorder in each region is on the x axis in 5% wide bins . we note that in each sub - plot the lines for all of the domains tend to cross at 1 point , giving the appearance of an isosbestic point , suggesting that there are at least 2 independent states in each of the 5 sections . a cluster analysis of these sections ( not shown ) indicates 2 components in each section , with the exception of the  spot  representing ordered proteins . second , in the upper left quadrant , there are more eukaryota pfam members in the 95100% disordered group than in any other 5% wide group pfam sequences plotted in each section of figure 5 are plotted here as a function of percent disorder in each sequence , correlated with but not the same as the vsl2b parameter . we note that the ordering of phylogenetic domains , with respect to increasing pid , is preserved with the exception of the upper right quadrant where the sample size is too small to be significant . the shifts in pid from 1 domain to the next are more quantitative here , and each  quadrant  appears to represent 2 distinct states . as discussed above , some methods of predicting intrinsically disordered structure are subject to end effects , such that , for vsl2b here , n - and c - terminals are predicted to be disordered when they are possibly not . this problem probably accounts for part of the approximately 20 percent false positive rate observed for some disorder predictors . as can be seen in the top left frame of figure 6 , corresponding to the lower - right spot in figure 5 bottom center , in that group of proteins that are the least disordered , only 24% figures  7 and 8 show structure - function feature maps for 3 examples of the 618 mammalian proteins found here with 100% predicted disordered pfam members , chosen on the basis of their known involvement in human disease and presented in order of length . the program that generates these figures was written for this work and is now part of the disprot pondr predictors available on the disprot site . structure - function maps of 2 100% pid whole proteins , each dominated by 1 pfam domain , are shown above with pfam sequences marked in green . many similar examples exist where parts of these proteins have been crystallized , but in the presence of sds or in co - crystals with partner molecules , illustrating the induced conformation nature of idps . pf05160 dss1 human has been co - crystallized in complexes : 1iyj , 1mje , and 1miu ) . figure  8 . a structure - function map of apc human shows an example of a protein with a mix of pid and structure with multiple pfam domains ( green ) , most predicted to be mostly disordered ( blue ) . x - ray crystal structures ( black ) also include many residues not observed in these structured segments ( red ) . the longest fragment of apc human to be obtained in large quantities is seen to be 100% disordered in the nmr solution structure ( also red ) . many of the x - ray structures for apc human have been co - crystalized with other molecules , suggesting that these conformations are induced by the binding of an idp with a structured partner . included in each figure are the idr prediction profiles from the pondr vsl2b , vl3 , and vxlt methods , markers for : pfam family and domain members , predicted idrs based on vsl2b , predicted molecular recognition features based on the vxlt profiles , regions of these proteins represented in the protein data bank by x - ray crystal diffraction and nmr evidence both for order and disorder , sites of phosphorylation and methylation , and for sequence variations related to human disease . these are described below with a comparison to other members of the same pfam family . pf05160 : dss1 human plays a role in ubiquitin - dependent proteolysis , interacts with the c - terminal of brca2 , and is involved in split hand - split foot malformation . other members of this family predicted here to be 100% disordered are dss1mouse-3 - 63 , sem1-drome-15 - 75 , and sem11arath-8 - 70 . the following members of the pfama.seed set predicted to be less disordered are : q9liy2orysj-9 - 97 , 86% , sem1-yeast-19 - 84 , 67% , and dss1-schpo-2 - 66 , 91% . in a limited sample some non - seed members show even lower amounts of pid , but we can not say that it is generally true that some family members have widely differing amounts of pid . pf04440 : dbnd1 human-14 - 153 binds to - and -dystrobrevin in muscle and brain , and genetic variation is thought to be associated with schizophrenia . other members of this family in the seed set : dbnd1-mouse-14 - 155 , 94% pid , and dbnd2human-100 - 254 , 99% pid . pf05923 : apc human , found repeated in the mid region of the adenomatous polyposis proteins ( apcs ) , near many cancer - linked snps . most other v23.0 seed members of this family in human are 100% pid : 13691394 , 1840 - 1866 , 2006 - 2031 , and 1948 - 1973 , but 1485 - 1510 is 85% pid and 1636 - 1661 is 27% pid . the pfam v23.0 seed set appears to be accurate with respect to most of the idps we have sampled . however , there are , as with all prediction methods , some inconsistencies , and we note that we are applying a prediction to a prediction here . figure  9 shows , on the right , a parallel analysis plot indicating that there are probably 3 independent types of mostly or entirely disordered pfam members where differences are based only on the averaged physical properties of their amino acids , and vsl2b scores are greater than 0.5 . likewise , on the left , we can not reject the null hypothesis that there are 5 types of pfam members in mammals based on the averaged physical characteristics . retention and rejection in the glorfeld principle component analysis of pfam sequence information represented by the 10 kidera factors . right : 3 components are retained when the sample is restricted to 1913 members here where the vsl2b parameter is greater than 0.5 ( more than half of each member is predicted to be disordered ) . as can be seen in the figures below , these components form well separated groups with very little overlap . in an earlier study aimed specifically at finding different flavors of intrinsic protein disorder in the swissprot database , vucetic et al . found 3 groups characterized by their amino acid composition . these groups , composed of predicted intrinsically disordered regions divided into windows 41 residues long , were related to function . it is not possible to compare those groups to the ones seen in this work with respect to physical properties . the question here is : how do the clusters found here relate to pfam families and function . there is substantial evidence that functional intrinsic disorder , particularly in long sequences , is composed of multiple shorter structural features , such as those found in morfs . to be meaningful , the mapping of types of disorder to function or pfam family should include a windowing analysis , similar to that performed earlier , of shorter segments that includes the kind of periodic features used to train pondr in this respect . figure  10 shows a pairs graph ( a scatterplot of all factors against all others ) of the mean kidera factors for pfam members where the average vsl2b score is greater than 0.5 , indicating that the sequence is likely to be more than half disordered . as in figure 9 , red and green delineate 100% pid members from others that are between 50 and 100% pid . this is where the vsl2b score , the predictor of intrinsic disorder , is highest . the 100% pid members here appear above the correlation line in the ext vs fle plot , and in the lower left of the dbe vs alp plot . however , these are not particularly special cases ; there is some degree of separation in the predictions of partly and completely disordered structure by all pairs of factors , indicating again that all factors contribute to this distinction . also , in the scatterplot matrix figure below , where the vsl2b score is included and the color distinction is made for clusters and not for degrees of disorder , the vsl2b score does not appear to be well correlated with either color mapping or with the overall trend of the scatter . another way of viewing the clustering , and non - clustering , of 100% pid pfam members is shown in the dendrogram plots below . a pairs graph shows scatterplot , regression lines , and typographically scaled absolute values of the correlation coefficients between the mean kidera factors for the 1913 members where the mean vsl2b parameters are greater than 0.5 ( half disordered or more ) . points for domains that are 100% pid are red , while all others are green , showing several visually distinct clusters in each of the 2-dimensional plots in a way that is almost impossible with 2-dimensional projections of high - dimensional objects , and the corresponding dendrogram shown in figure 13 , right , does not reveal the level of structure shown here . helix / bend ( hel ) and extended structure preferences ( ext and fle ) are negatively correlated , and the 2 extended structure preferences are positively correlated , as can be seen also in figure 11 . note that although the kidera factors themselves have zero correlation , some factors averaged over the entire lengths of pfam sequences do show substantial correlation , albeit with large variance . there is a positive correlation here between average ext ( extended ) and fle ( flat extended ) factors , and a negative correlation between these and the average hel ( helix ) factors for entire pfam sequences . a biplot of the data calculated from 10572 pfam sequences and 10 mean kidera factors shows a 2-dimensional projection of 3 of the 5 principle components , vectors of equal length representing the kidera factors ( table  1 ) , and 100% pid sequences colored in red . note that the apparent correlation between hph and sur is caused by the projection into 2 dimensions , and that the actual correlation is 0.25 ( fig . these distinctions become even more evident when only shorter pfam sequences , or windows of uniformly shorter segments , are chosen ( not shown ) . as pfam sequences become longer , the means of the kidera factors tend toward the central value , zero , masking the diversity of idp types in pfam domains . we anticipate that , for example , a cluster will appear along the hel ( helix ) or alp ( occurrence in  region ) axis in an analysis that includes windows of 20 residues each , and that shorter windows may yield a cluster along the dbe ( double - bend preference ) axis . figure  11 shows a pca biplot of all pfam members here scored by the principle components of the 10 mean kidera factors for the members . each point represents an identifiable member with eigenvalues or coefficients corresponding to how much each component contributes to the variance in the original data . the 10 labeled vectors are projections of the mean kidera axes onto the 2-dimensional plot shown here , indicating their contributions , however difficult to see in only 2 dimensions , to the 3 principle components represented here . the positively correlated average extended- and flat - extended factors point away from the average helix factor to which they are negatively correlated , and these are somewhat aligned along the pc1 and ( for fle ) pc3 axes . note that 100% pid members , in red , are clustered primarily along the fle axis , but that some occupy positions at the far ends of other factors . all of the kidera factors appear to contribute to the preference for types of order or disorder . a better visualization of how all 10 factors contribute to the principle components , particularly when there are more than 3 pcs , can be obtained by rotating the pca biplot in 2- and 3-dimensional projections in real time using the rggobi and ( for windows only ) biplotgui r packages . both there , as well as in static plots shown here , individual points can be labeled with their accession numbers to better explore these relationships . figure  12 shows a partial least squares and principal component regression of the mean vsl2b score against the 10 mean kidera factors for all pfam members included here . each point on the left represents a pfam member scaled by 6 eigenvectors of the pca decomposition . plots of the regression coefficients on the right show that only 5 components show noticeable differences , so the regression here is essentially in 5 dimensions . the distribution of lengths in pfam members here is broad , causing the mean values of the kidera factors to have a limited value , tending toward zero a central value where lengths are long . a partial least squares and principal component regression of the mean vsl2b factor against the 10 mean kidera factors for 10572 pfam sequences , left , shows a prediction of the mean vsl2b factor for pfam sequences using  leave one out  validation , with 100% pid sequences in red . here the x axis represents mean vsl2b scores calculated directly from the vsl2b predictions for each sequence ( these are just mean vsl2b scores ) and the y axis represents predictions from the multivariate linear regression on 6 principle components of the kidera factors ( see methods ) . we hypothesize that the nonlinearity and spread of the data here is partly due to errors in the vsl2b predictions themselves , that the kidera factors more accurately represent the tendency to disorder , and that the spread will narrow in an analysis of smaller uniform windows of sequence . five components yield the same results , consistent with the results shown in figure 9 , and here to the right . the plot on the right shows the regression coefficients for these 6 sets of principle components . here the sets containing 5 and 6 components have nearly the same coefficients , indicating , as was also shown using a different analysis in figure 9 , that the sixth component contributes little to the information contained in the kidera factor averages . we hypothesize here that convergence will shift to more coefficients when smaller uniform windows are analyzed . figure  13 shows dendrograms of the mean kidera and vsl2b factors for the complete 10532 member set , and for the 1913 member set that includes only members with mean vsl2b scores above 0.5 . it appears to be clear that 100% pid members are in distinct groups . depending on where the y axis line is drawn , it is possible to find between 5 and 7 clusters here . however , in the 1,913 member subset , although we know that there are only 3 distinct clusters ( fig . the dendrogram on the left shows the mean kidera and vsl2b factors for 10532 pfam members , with clustered 100% pid members in red . five clusters are obtained here by drawing a horizontal line crossing six of the vertical markers near y = 2.1 . the marker on the far left is an outlier . the dendrogram on the right shows the same mean factors for 1913 pfam members where the vsl2b mean is above 0.5 and sequences are half pid or more . the 100% pid members plotted in red are not well clustered , indicating that here we can not say they are different from those that are half pid . when sequences are chosen to be on the order of 10 residues , distinct clusters do appear ( not shown ) . we now turn to a principle component analysis ( pca ) of the kidera factors averaged over entire pfam sequence lengths , in figure 14 below where the mean vsl2b scores are excluded , and in figure 15 where they are included to show an almost complete lack of correlation with any single mean kidera factor . figure  14 shows a bivariate cluster plot , showing only 2 components , of the 1,913 member subset where 4 components are retained to show the tendency of long members to have the same kidera factor means , reducing the dimensionality available for discriminating between disorder types . the 4 clusters have the same colors in both the cluster and silhouette plots . while values around 0.2 in silhouette plots are low , indicating weak clustering , these values become larger as longer members are excluded . figure  14 . left : a bivariate cluster plot of a partitioning around medoids ( pam ) for k = 4 clusters shows 3 clusters with almost no overlap , and 1 cluster with considerable overlap ( black ) . the observations near zero tend to be from long pfam sequences where means in the distance matrix tend to the central value . this calculation included the mean vsl2b factors and the 10 mean kidera factors for 1913 pfam members where the mean vsl2b factor is greater than 0.5 . when k is set to 3 , the black and blue points are combined into one group . right : a silhouette plot of the same pam object using the same colors also shows 3 mostly non - overlapping groups and one with significant overlap . silhouette widths : near 1 indicate well clustered groups , near 0 indicates that observations lay between 2 clusters , and negative means that observations overlap or are in the wrong cluster . scatterplot matrix of the mean vsl2b factors and 10 mean kidera factors for 1913 pfam members where the mean vsl2b factor is greater than 0.5 . points plotted here have the same colors used in the pam cluster analysis plotted in figure 14 . the scatter plot shown here for vsl2b ( vsl ) and hydrophobicity ( hph ) is essentially the same plot shown in figure 5 , with the exception that here only the left side of figure 5 is represented . note that the cluster identified in red is most prominent where size is plotted against the other parameters , and that it runs parallel to the vsl2b axis , indicating that sequences in this cluster containing large amino acids are present in both ordered and pid pfam members . the cluster identified in green is most prominent when helix preference is plotted against other parameters , marking sequences with low mean helix preference , and for blue sequences with high mean helix preference . the distribution of residue lengths has a distorting effect on cluster plots . while residue lengths of most pid sequences are less than 10 ( 4 ) , the mean length of 100% pid pfam sequences is 82 , the median is less than 50 , and many are longer than 150 ( fig . 1 ) . these much longer sequences tend to have mean kidera factors near zero , even though they may contain regions that vary from this and from each other significantly . presently not well supported by the statistics , both the blue and red clusters on the left ( fig . 14 ) can be divided again to form 5 well separated clusters , and we think these will become more significant as a windowing scheme is imposed on the size of cluster members and the members colored in black become redistributed . figure  15 shows a scatterplot matrix for the 1913 member set where mean vsl2b factors greater than 0.5 . members plotted here have the same colors used in the pam cluster analysis plotted in figure 14 , but no pca is performed here . again , these plots take on different characteristics as sequence lengths become more uniform ( not shown ) . perhaps the most notable feature of this figure : the cross - plot for vsl and hph , showing the relationship between vsl2b score and hydrophobicity , is essentially the same as is shown in the top left of figure 5 where cdf is less than zero and ch is between -0.4 and 0.6 . note that the cross plots for vsl against all of the other kidera factors show a similar spread of data , indicating that hydrophobicity is not the only factor important in predicting intrinsic disorder . as in figure 10 , a positive correlation can be seen here between extended and flat extended structure , but here there is also a clear separation of clusters ( colors ) along the axis of the correlation . also , the negative correlation between the 2 extended factors and helix shows a separation of clusters . note in this figure that amino acid size ( siz ) and partial specific volume ( psv ) have very different distributions plotted against predicted disorder ( vsl ) . at least from this perspective , size divides both folded and disordered pfam sequences along the vsl axis in a very clear way , while volume , one of the factors used in pondr training , does not appear to discriminate . likewise , hydrophobicity ( hph ) , while showing more of a separation of the three primary clusters , does not show a direct relationship with vsl . the single kidera factor that appears to be in the most direct relationship with vsl here is pkc where , on the top line of plots , blue predominates in the upper right and green somewhat in the lower left . the kidera factors contain most of what is known about the statistically countable and experimentally measurable physical characteristics of the amino acids , effectively increasing the information content of the sequence by 10-fold . principle component and cluster analyses of the averages of these factors for each pfam member , included those that are folded , reveal 6 separate groups . in the subset of pfam members where the vsl2b disorder predictor is greater than 0.5 at least 3 groups are found , despite the homogeneous distribution of pfam member lengths . in the earlier study aimed specifically at finding different flavors of intrinsic protein disorder , vucetic et al . found three groups characterized by their amino acid composition in windows 41 residues long in predicted intrinsically disordered regions in 80,000 sequences of the swissprot database . it is not possible to compare those groups to the ones seen in this work with respect to physical properties . however , we can observe that the overlap of clusters in the previous work is considerable , while here in figure 14 the overlap of 1913 sequences is very small but the clustering is weak in the usual interpretation of the silhouette plot . where sequences are long the average of the kidera factors tend toward zero , limiting the capacity of the analyses to discriminate . it is hypothesized that grouping sequences into smaller subsets with more uniform lengths will reveal more information , and lead to more finely structured and perhaps more accurate predictions of intrinsic disorder . adding additional information derived from sliding window calculations based on the physical properties , related to interactions such as is provided by hydrophobic moment , and flexibility , strengthens predictor training , particularly for features important for molecular recognition . tables 8 and 9 show the percent of non - mammal single and multiple cell eukaryotes , respectively , with unique pfam members that are 100% pid . it has been shown previously that the average fractions of disordered residues in unicellular and multicellular eukaryotes are about the same , with unicellular eukaryotes having more scatter ( ref . 1 ) . we observe here , in a relatively small sample size , that single cell eukaryotes are predicted to have fewer species with 100% intrinsically disordered pfam domains than multiple cell eukaryota . this may be an artifact of the selection of pfam seed sequences in these 2 groups , or it may have some evolutionary significance . only species represented by 100 or more pfam domains in the seed database are included . ndis : number of proteins ; 100% pid ; tot : total number of proteins analyzed ;    ...  indicates where not all levels of taxonomy are listed . a standing hypothesis among idp investigators proposes that intrinsically disordered protein ( idp ) has evolved into different classes that can be identified by physical characteristics or functions , analogous to those identified by hmms in the pfam database . we analyze distributions and clusters of pid in 193024 members of the version 23.0 pfam seed database , representing 12456 unique domains , families , and repeats . of these sequences , 616 mammalian members associated with 315 biological functions , and 120 of the parent whole protein set , are 100% pid . exploring ways to find the maximum information content in intrinsically disordered proteins independent of sequence information used to train hmms , we applied the 10 linearly independent kidera factors containing most of the variance in the physical properties of the amino acids to a transformation and analysis of sequence in pfam training set members . the hypothesis that intrinsic protein disorder exists in multiple forms with preferences analogous to those for protein secondary structure , as is evidenced in the success of morf predictors is supported by this analysis . there are 3 principle conclusions : ( 1 ) each of the 10 orthonormal kidera factors contributes to intrinsic protein disorder . the pondr predictors used here , while still among the most accurate , included less than 2 of the kidera factors in training . the inclusion of all 10 kidera factors could hypothetically increase the information available to a predictor by more than 5-fold . ( 2 ) we identify 3 clearly separate and non - overlapping clusters of intrinsically disordered mammalian pfam sequences where the vsl2b score for disorder is not clearly correlated with cluster parameters , and where hmm training or sliding window calculations are not involved . the distributions of length in these pfam sequences are broad and skewed , so these 3 clusters are groupings of intrinsically disordered pfam sequences ( or families ) , and not of types of intrinsic disorder . a sliding window analysis limiting sequence size to uniform lengths characteristic of morfs and other features important to idp function could hypothetically increase the number of non - overlapping clusters several fold and reveal a meaningful correlation between disorder cluster type and function . the next step is to train an intrinsic disorder predictor against the kidera transformed disprot database , with the addition of these other methods for recognizing patterns in sequence . ( 3 ) while mammalian proteins in the pfam seed database contain 40% pid sequence here , only 27% of the pfam members in this set contain pid regions , suggesting that the pfam database is missing functionally important pid protein domains . the kidera factors consist of 10 linearly independent vectors derived from an eigenvalue  eigenvector decomposition of statistically and experimentally determined physical properties of the amino acids . the frequencies of occurrence of each amino acid in experimentally known intrinsic disorder , from the disprot database , can be added to the kidera matrix to make 11 columns , but a singular value decomposition of this matrix shows only 10 linearly independent vectors ; the preference for intrinsic disorder is contained in the 10 kidera factors . this analysis included 193024 members of the version 23.0 pfam seed database , representing 12456 unique domains , families , and repeats . here refers to individual pfam domain member sequences and not to entire pfam domains . the term  domains  here is used for both pfam domains and phylogenetic domains , but the distinction is spelled out . the term  regions  here refers to parts of proteins , or parts of pfam members . the term  sequence  here refers to protein in whole or part , and if in part , whether it is isolated or not is spelled out . whole protein sequences represented in the pfam seed set , and the human , mouse , and chimp proteomes were obtained from uniprot . predictions of intrinsic disorder were performed on 193024 whole proteins using predictors vsl2b , vl3 , and a method that is predictive of molecular recognition features , vlxt . the vsl2b predictor used for the statistical results here compares well with other methods , available on the dp site . statistics , cluster , factor , and principle component analyses , were calculated using r core and contributed packages . the r implementation of the wilcoxon rank sum test with continuity correction is used to compare distributions , testing the null hypothesis that the distributions of x and y differ by a location shift of  mu .  the alternative is that they differ by some other 1- or 2-sided location shift . rejection of a null hypothesis in some cases sampled here may have no practical significance if other larger effects are present , such as the possibility that populations of pfam domains may contain artifacts of hmm training methods , as it appears to be a factor here with respect to low complexity sequences . confidence intervals on the results of a singular value or eigenvector  eigenvalue decomposition can be calculated when the standard error in the data noise is known because the scale of the original data are maintained through the transformation and eigenvectors can be compared with noise . again , even if the level of the noise were known here , other larger effects may render it meaningless . to obtain some statistical confidence in our results , a parallel analysis was performed using an eigenvalue - eigenvector decomposition and glorfelds principle component analysis of the 10532-pfam - member by 10-mean - kidera - factor matrix . each pfam family or domain member was represented by a vector of the means of the 10 kidera factors listed in table 1 over each entire pfam region . paran options included 1000 and 5000 iterations , giving the same results at the 99th percentile . a common factor analysis ( not shown ) retained seven factors from the 10532 sample set . no real difference was seen in the number of eigenvalues retained when vsl2b scores are added to the kidera set . a partial least squares and principal component regression , plsr , of the mean vsl2b factor against the 10 mean kidera factors for 10572 pfam sequences were calculated using the following model : vsl hel + ext + dbe + e + pkc + siz + hph + psv + alp + sur . structure and function information extracted from uniprot database or fasta format files is combined with predictions of intrinsic disorder here to provide graphical representation of the relationships between these features . supported by department of defense high performance computing modernization program project number odefn26263101 and usu project c02928 .
the pfam database groups regions of proteins by how well hidden markov models ( hmms ) can be trained to recognize similarities among them . 
 conservation pressure is probably in play here . 
 the pfam seed training set includes sequence and structure information , being drawn largely from the pdb . 
 a long standing hypothesis among intrinsically disordered protein ( idp ) investigators has held that conservation pressures are also at play in the evolution of different kinds of intrinsic disorder , but we find that predicted intrinsic disorder ( pid ) is not always conserved across pfam domains . here 
 we analyze distributions and clusters of pid regions in 193024 members of the version 23.0 pfam seed database . to include the maximum information available for proteins that remain unfolded in solution 
 , we employ the 10 linearly independent kidera factors13 for the amino acids , combined with pondr4 predictions of disorder tendency , to transform the sequences of these pfam members into an 11 column matrix where the number of rows is the length of each pfam region . 
 cluster analyses of the set of all regions , including those that are folded , show 6 groupings of domains . 
 cluster analyses of domains with mean vsl2b scores greater than 0.5 ( half predicted disorder or more ) show at least 3 separated groups . 
 it is hypothesized that grouping sets into shorter sequences with more uniform length will reveal more information about intrinsic disorder and lead to more finely structured and perhaps more accurate predictions . 
 hmms could be trained to include this information .