BACKGROUND
intrinsically disordered proteins  and intrinsically disordered regions  within structured proteins are defined by the lack of a stable tertiary structure and a corresponding high degree of flexibility under physiological conditions  <cit> .
the importance of conformational flexibility is reflected in the observation that idps and proteins with idrs are often involved in essential cellular processes, such as cell-cycle regulation, transcription, and translation  <cit> .
additionally, they often play major roles in pathologies associated with aggregation and misfolding  <cit> , making them attractive potential drug targets  <cit> .
genes encoding such amino acid sequences are under reduced selective pressure, which is manifest in a higher sequence diversity compared to genes of structured proteins/domains  <cit> .
whereas the functionality of a protein segment is often approached by investigating the evolutionary history of its primary sequence  <cit> , this is often difficult to achieve with idp/idr sequences, due to their generally high sequence diversity  <cit> .
on the other hand, combining the information derived from analyzing the conservation of both sequence and disorder can be much more useful, and this idea has recently been suggested to partition disordered residue positions into three separate groups of potentially different functional attributes: i) ‘constrained’, if both features  are conserved; ii) ‘flexible’, if only disorder is conserved; and finally, iii) ‘non-conserved’ positions where disorder is not conserved.
these specific evolutionary behaviors have been shown to correlate with distinct disorder-related functional categories  <cit> .
in general, segments of constrained disorder are often associated with protein binding and molecular recognition, whereas flexible disorder is prevalent in linker segments acting as entropic chains.
non-conserved disorder has not been associated with specific protein function so far.
here, we present discons, a novel web application and downloadable, stand-alone source code that offers a description of the conservation of both the amino acid sequence and of the feature of structural disorder, and performs the classification of disordered positions into these three categories.
thus, discons provides an additional  layer of information that together with other sequence-based tools, such as the paml software package  <cit> , morfpred  <cit>  and anchor  <cit> , should facilitate the effective identification of functionally important disordered regions in proteins.
implementation
both the web application and the downloadable version of discons hosted at the website are freely available without registration.
the discons website is divided into four functional sections that are accessible both through the menu and via the options shown on the welcome page.
these sections correspond to the three running mode interfaces: ‘quick’, ‘advanced’ and ‘from alignment’.
the fourth functional section is ‘help’, which offers the complete documentation of the server and source code, in addition to a user guide.
the ‘quick’ running mode requires a single protein sequence  as input, or alternately, a uniprot  <cit>  accession id.
in this mode, the default parameters are used through all the calculation steps of the discons workflow.
although this calculation is the easiest to set up, experienced users might prefer to use the ‘advanced’ tool enabling a better understanding of the results leading to more fine-tuned functional interpretations.
the ‘advanced’ mode also accepts a single protein sequence or uniprot id in a manner similar to the ‘quick’ calculation, but in this mode users can manually set all the parameters of the underlying calculation, allowing for a detailed optimization of the protocol pipeline, and a better overall command of the final results.
finally, the ‘from alignment’ mode is best suited if the user already has a custom made, reliable multiple sequence alignment that can be used for the calculations.
the main advantage of this mode is speed, since the need for running a blast search and constructing the multiple sequence alignment with mafft  is circumvented.
therefore, this mode is significantly faster than the others which start from a single sequence .
by default, the stand-alone source code is also running ‘from alignment’; however if the necessary dependencies, namely blast+  <cit>  and mafft  <cit>  are available locally in the user’s computer, the full pipeline can be utilized in a straightforward manner.
depending on the running mode, the workflow of the calculations has a different starting point .
in ‘quick’ and ‘advanced’ modes, the procedure starts with a blastp or psi-blast search  <cit>  to collect sequences similar to the query sequence.
in ‘advanced’ mode, the search dataset  or pdb  <cit> ) and the blast threshold values can be specified.
next, a multiple sequence alignment  is created from the set of identified homologous sequences using mafft  <cit> .
this is the most crucial part of the procedure, since aligning disordered regions is non-trivial due to the potential diversity of related sequences.
because an incorrect alignment will compromise the subsequent calculations, it is advised either to use the ‘advanced’ mode to fine-tune the alignment procedure or to use a reliable, user-defined multiple alignment in the ‘from alignment’ mode.figure  <dig> discons work flow.
the schematic work flow of the discons sequence analysis pipeline is displayed along with a sample output table, showing the discons profile of human calpastatin domain i .
the procedure starts either from a query sequence, followed by a blastp or psi-blast search and the aligning of the retrieved sequences, or from a user-provided multiple sequence alignment.
next, position-specific sequence- and disorder conservation scores are calculated and finally these pairwise scores are combined to create the  discons profile.
the result of one such classification is displayed at the bottom of the figure.
in the next step , the msa is used to construct an aligned disorder profile by running iupred   <cit> , vsl <dig>  <cit> , espritz  <cit>  or foldindex  <cit>  on each of the aligned sequences.
the disorder scores are first transformed to a binary scale  for each sequence.
details about the transformation procedure for the three different disorder predictors are given in the help section of the website.
next the fraction of disordered residues at every position in the msa across the different sequences is calculated, thereby effectively quantifying the position-wise conservation of disorder.
next, sequence conservation scores for each position in the alignment are calculated using the algorithm developed by capra et al.
<cit> .
in ‘advanced’ and ‘from alignment’ modes, a number of parameters such as the algorithm of choice or the background distribution can be specified to make the calculation more robust.
finally, positions in the msa are scored by both their sequence- and disorder conservation, ranging from  <dig>  to  <dig>  .
based on these pairs of scores, each position will fall into one of four distinct categories.
as suggested by bellay et al.
<cit> , positions with a higher degree of disorder can be ‘constrained’ , if both the sequence and disorder conservation scores are  <dig> or greater; ‘flexible’ , if the sequence conservation is lower than  <dig> but the disorder conservation is  <dig> or greater; or ‘non-conserved’ , if the disorder conservation is lower than  <dig>  but higher than  <dig>  positions with a disorder conservation score of  <dig> are completely lacking disorder and therefore are considered as ‘structured’ .
thus, regions of constrained disorder show strong conservation both at the amino acid level, and also of the disorder feature, while flexible disordered regions are variable in terms of amino acid sequence, but retain a significant level of disorder in evolution.
lastly, non-conserved disordered regions lack disorder as a conserved feature, and are generally thought not to be associated with functions  <cit> .
on the results page, the position-specific conservation profile is provided on the output screen, and the fractions of residues falling into each of these distinct categories are also displayed in a tabular format at the bottom of the page, effectively quantifying the conservation of disorder in the query sequence .
the sequences of consecutive ‘constrained’ disordered regions are also recorded and are available for download in fasta format along with the profiles and fractions in text format using the links that are provided.
such segments of consecutive stretches of constrained disorder are most likely to correspond to functionally important idrs such as linear motifs or morfs, as we describe below.
RESULTS
bench-marking on molecular recognition features
we evaluated the performance of discons on a set of molecular recognition features, or morfs, that are short peptide segments mostly found within longer disordered regions  and involved in the binding to protein partners via disorder-to-order transition  <cit> .
morfs have been implicated in functions involving regulation and signaling, among other cellular processes.
these recognition features are enriched in disordered residues, however, they may also have some residual  structure, and their sequences are relatively more conserved than their flanking disordered regions  <cit> .
to estimate the efficiency of the discons protocol in distinguishing between such functionally important disordered segments on a large scale, we retrieved three morf datasets from morfpred that are available at their website  <cit>  and combined them into a single benchmarking dataset.
the three datasets were the ‘test dataset’ containing morfs deposited in the protein data bank  before 2008; the ‘experimental dataset’ with morfs identified between  <dig> and 2012; and the ‘test 2012’ dataset with morfs from  <dig>  the combined dataset contained  <dig> morf instances.
after applying a sequence redundancy filter on the full length sequences using cd-hit  <cit> ,  <dig> unique sequences remained.
morf sequences were extracted from the full length protein sequences along with up to  <dig> residue long flanking segments on both sides.
disorder propensity scores were then calculated for the extracted morfs, the flanking residues, the full length proteins, and the complete uniprot/swissprot database, using iupred.
figure 2a displays the distribution of disorder scores, comparing these four datasets, demonstrating that morf residues  and morf-flanking regions  are more disordered than morf-containing proteins .
the difference is even more pronounced when compared to proteins from the complete uniprot/swissprot database .
since the distributions did not follow gaussian or even symmetric distribution, we chose the non-parametric kolmogorov-smirnov  test, which only assumes that the compared variables are continuous.
the distributions of morfs and flanking residues are significantly different from both their proteins and the uniprot/swissprot dataset, according to ks tests with p-values less than the precision limit of r , indicating very strong significance.figure  <dig> analyzing morfs with discons.
morfs are known to be enriched in disordered residues.
panel a shows normalized density distributions of morf regions, morf-flanking segments, full-length sequences of the morf-containing proteins and all proteins from swiss-prot.
in these distributions, the area under each curve adds up to one.
morf residues  are predicted to be more disordered on average than the proteins they are found in, and especially more than the proteins of uniprot/swissprot.
the flanking regions  of the morfs also have significantly higher disorder content, compared to the full-length proteins they are found in.
panel b shows the combined disorder- and sequence conservation scores, which range from  <dig>  to  <dig> .
by comparing these residue-specific score pairs, discons further supports the idea that the sequences of morfs  are more conserved than that of their flanking regions ; however, even in the flanking regions intrinsic disorder as a feature is highly conserved, indicating that these segment are required to be flexible in order for the protein segment to function.
after performing psi-blast  <cit>  searches against each morf and their flanking regions, the conservation of the aligned positions were quantified in terms of sequence- and disorder conservation according to the protocol explained in the ‘implementation’ section.
the binned conservation score pairs are displayed on figure 2b, comparing morfs  and flanking residues .
in order to statistically compare the conservation scores, we used the welch t-test, since this test only assumes gaussian distribution of the variables, and does not require equal variances.
the conservation of disorder in morf  and flanking  residues is similar , while the underlying amino acid sequence is significantly more conserved in the case of morf residues .
therefore, the comparison of morfs and morf-flanking sequences shows a trend that is in agreement with the literature  <cit> , namely that while both sequence and disorder are rather conserved in the morfs, their neighboring protein segments are less conserved sequence-wise.
bench-marking on short linear motifs
following the analysis of the morf dataset, we applied the discons procedure on all the  <dig> known instances of short linear motifs  of the elm database  <cit> .
generally, these motifs are enriched in disorder, and their sequences show higher than average conservation on the amino acid level  <cit> .
we compared these motifs to all the available idp sequences retrieved from disprot  <cit>  .
the average disorder content of slims is significantly higher than that of the full length idps , and in fact even more so than in the morf dataset  .
as expected and demonstrated in figure 3b, both the sequence and disorder conservation scores of the slim sites are significantly higher than that of the full length idps , with 52% of all the slim residues being of constrained disorder compared to only 12% in the full length idps, indicating the importance of structural disorder in slims.
in comparison to morfs, where 36% of the residues are of constrained disorder, structural flexibility seems to play a larger role in slims, and indeed, morfs are known to often have some residual pre-formed structural elements, while slims are more disordered overall.
additionally, the majority of slims are localized on consecutive segments of constrained disorder.
concretely, 88% of the slims were found within consecutive segments of  <dig> or more constrained residues,  <dig> % in segments that are at least  <dig> residues long, and  <dig> % in segments of at least  <dig> constrained residues in length.
since slim segments often form ligand/protein binding sites, their high sequence conservation is necessitated by the formation of interface contacts with partner proteins.figure  <dig> analyzing slims with discons.
short linear motifs  are peptide motifs with characteristic sequence patterns and are generally enriched in disorder.
as displayed on panel a with the help of normalized density curves, slims are significantly more disordered than the full length proteins they are found in or compared to the proteins found in disprot.
not only are slims more disordered, but this flexibility is highly conserved across homologs as well .
when compared to the idps of disprot  the difference in conservation is striking.
as in figure 2b, disorder- and sequence conservation scores range from  <dig>  to  <dig> .
case studies of discons uncovering constrained and flexible disorder
finally, we provide two case studies using two different protein segments; one exemplifying “constrained” disorder and the other being an example of the “flexible” disorder class.
figure 4a shows a morf region of constrained disorderfound in the c-terminal negative regulatory domain of the p <dig> protein  bound to the s <dig> calcium-binding protein , along with the discons profile of the c-terminal part of the p <dig> sequence .
all the residues forming the morf are constrained based on the conservation profile of sequence- and of structural disorder, and, not surprisingly, these are the only residues of the disordered segment that appear in the crystal structure.
the interaction with s <dig> restricts access to phosphorylation and acetylation sites on p <dig> that are important for transcription activation  <cit> .
thus, this region that is important for mediating a critical interaction is clearly identified by our protocol as a conserved disordered segment.figure  <dig> example of a constrained disordered morf.
the intrinsically disordered c-terminus of p <dig> adopts a helical fold upon partner binding .
this is the only segment of the disordered region that appears in this crystal structure.
upon quantifying the conservation of sequence and of disorder in the c-terminus, the full length morf segment is classified as ‘constrained’ disordered, where both features are highly conserved .
since ‘flexible’ disordered segments generally function as entropic chains linking structured  segments, finding structural data for them is less straightforward.
these residues are often missing from the structures found in the protein data bank   <cit> , however, ensemble descriptions of such regions are available from the protein ensemble database   <cit> .
in order to demonstrate “flexible” disorder, we retrieved the “fuzzy” complex formed between sic <dig> and the cdc <dig> subunit of an scf ubiquitin ligase .
as seen generally in “fuzzy” complexes, sic <dig> is a fully disordered protein that remains disordered even when bound to its partner.
sic <dig> has multiple binding segments along its disordered chain that compete for binding to the same pocket on the receptor protein.
figures 5a and b display two conformations  present in the ensemble description of the complex, while figure 5c provides the discons profile of the corresponding segment of the sic <dig> protein.
the profile clearly shows that a “flexible” disordered linker connects the two “constrained” disordered binding regions.
in several of the  <dig> different conformations that constitute the ensemble of the “fuzzy” complex, both binding regions are found to contact the binding pocket, whereas the linker segment does not bind to cdc <dig> in any of the conformations.
this indicates that the ‘unstructured-ness’ of the linker is more important for the function of the protein than the corresponding amino acid sequence of this segment.figure  <dig> example of a flexible disordered linker.
the intrinsically disordered protein sic <dig> binds to the cdc <dig> subunit of an scf ubiquitin ligase forming a fuzzy complex .
two conformations of the same sic <dig> fragment from the ensemble of this complex are displayed on panels a and b. below, the discons profile on panel c corresponds to the same fragment, and shows the segment that is flexible disordered between two constrained disordered regions  that bind to cdc <dig>  this former segment corresponds to a flexible linker , which is found between two morfs  that bind to the same binding pocket in the fuzzy complex.
the flexible conservation of the linker indicates that the sequence is of less importance than its conformational flexibility.
CONCLUSIONS
discons is a novel and freely available online and downloadable tool that combines the quantitative description of the position-specific evolutionary conservation of the amino acid sequence with predictions of the conservation of its disordered/flexible state, providing meaningful information on the evolutionary context of a disordered protein segment.
furthermore, discons uses this combined information to classify each disordered position into one of three categories, namely: constrained, flexible and non-conserved.
these classes have been suggested to correlate with distinct functions that arise from the disordered state; therefore discons may provide information orthogonal to those obtained by other methods, which potentially enhances the reliability of the identification of functionally relevant disordered segments within proteins.
we demonstrated that discons can be used to investigate both sequence- and disorder conservation in a functionally meaningful manner by bench-marking our procedure on morf and slim datasets, which are known to be conserved functional units enriched in structural disorder.
it is important to emphasize that the success of calculation with discons strongly depends on the quality of the underlying multiple sequence alignment; therefore it is advised to review and optimize each msa to maximize the information of the output.
taken this into consideration, discons can be used as an online or stand-alone tool for quantifying the conservation of both sequence and structural disorder by analyzing large-scale protein datasets and individual proteins.
as such, discons might provide an additional layer of information for the investigation of protein disorder, and could serve to enhance the performance of prediction software such as morfpred  <cit> , or provide descriptive information for disorder related databases such as d2p <dig>  <cit> , mobidb  <cit>  or ped  <cit> .
availability and requirements
discons is available as a web application, and as source code, both hosted at http://pedb.vib.be/discons.
the source code is written in python, and has two versions: the multiple sequence alignment -based script, and the complete pipeline.
the multiple alignment-based version has no requirements , while deployment of the full pipeline locally requires the following software: blast+  <cit>  and mafft  <cit> .
the software is distributed under the gnu gpl license.
abbrevations
idpintrinsically disordered protein
ldrlong disordered region
morfmolecular recognition feature
msamultiple sequence alignment
slimshort linear motif
competing interests
the authors declare that they have no competing interests.
authors' contributions
vm conceived the project and developed the web application of discons.
vm developed the underlying source code.
vm, fzs and mg tested the code and the web application.
vm and mg carried out the bench-marking analyses.
vm, mg, fzs and pt wrote the paper.
all authors read and approved the final manuscript.

 intrinsically disordered proteins and intrinsically disordered regions within structured proteins are defined by the lack of a stable tertiary structure and a corresponding high degree of flexibility under physiological conditions. the importance of conformational flexibility is reflected in the observation that idps and proteins with idrs are often involved in essential cellular processes, such as cell-cycle regulation, transcription, and translation. additionally, they often play major roles in pathologies associated with aggregation and misfolding, making them attractive potential drug targets. genes encoding such amino acid sequences are under reduced selective pressure, which is manifest in a higher sequence diversity compared to genes of structured proteins/domains. whereas the functionality of a protein segment is often approached by investigating the evolutionary history of its primary sequence, this is often difficult to achieve with idp/idr sequences, due to their generally high sequence diversity. on the other hand, combining the information derived from analyzing the conservation of both sequence and disorder can be much more useful, and this idea has recently been suggested to partition disordered residue positions into three separate groups of potentially different functional attributes: i) ‘constrained’, if both features are conserved; ii) ‘flexible’, if only disorder is conserved; and finally, iii) ‘non-conserved’ positions where disorder is not conserved. these specific evolutionary behaviors have been shown to correlate with distinct disorder-related functional categories. in general, segments of constrained disorder are often associated with protein binding and molecular recognition, whereas flexible disorder is prevalent in linker segments acting as entropic chains. non-conserved disorder has not been associated with specific protein function so far. here, we present discons, a novel web application and downloadable, stand-alone source code that offers a description of the conservation of both the amino acid sequence and of the feature of structural disorder, and performs the classification of disordered positions into these three categories. thus, discons provides an additional layer of information that together with other sequence-based tools, such as the paml software package, morfpred and anchor, should facilitate the effective identification of functionally important disordered regions in proteins. implementation both the web application and the downloadable version of discons hosted at the website are freely available without registration. the discons website is divided into four functional sections that are accessible both through the menu and via the options shown on the welcome page. these sections correspond to the three running mode interfaces: ‘quick’, ‘advanced’ and ‘from alignment’. the fourth functional section is ‘help’, which offers the complete documentation of the server and source code, in addition to a user guide. the ‘quick’ running mode requires a single protein sequence as input, or alternately, a uniprot accession id. in this mode, the default parameters are used through all the calculation steps of the discons workflow. although this calculation is the easiest to set up, experienced users might prefer to use the ‘advanced’ tool enabling a better understanding of the results leading to more fine-tuned functional interpretations. the ‘advanced’ mode also accepts a single protein sequence or uniprot id in a manner similar to the ‘quick’ calculation, but in this mode users can manually set all the parameters of the underlying calculation, allowing for a detailed optimization of the protocol pipeline, and a better overall command of the final results. finally, the ‘from alignment’ mode is best suited if the user already has a custom made, reliable multiple sequence alignment that can be used for the calculations. the main advantage of this mode is speed, since the need for running a blast search and constructing the multiple sequence alignment with mafft is circumvented. therefore, this mode is significantly faster than the others which start from a single sequence. by default, the stand-alone source code is also running ‘from alignment’; however if the necessary dependencies, namely blast+ and mafft are available locally in the user’s computer, the full pipeline can be utilized in a straightforward manner. depending on the running mode, the workflow of the calculations has a different starting point. in ‘quick’ and ‘advanced’ modes, the procedure starts with a blastp or psi-blast search to collect sequences similar to the query sequence. in ‘advanced’ mode, the search dataset or pdb ) and the blast threshold values can be specified. next, a multiple sequence alignment is created from the set of identified homologous sequences using mafft. this is the most crucial part of the procedure, since aligning disordered regions is non-trivial due to the potential diversity of related sequences. because an incorrect alignment will compromise the subsequent calculations, it is advised either to use the ‘advanced’ mode to fine-tune the alignment procedure or to use a reliable, user-defined multiple alignment in the ‘from alignment’ mode.figure discons work flow. the schematic work flow of the discons sequence analysis pipeline is displayed along with a sample output table, showing the discons profile of human calpastatin domain i. the procedure starts either from a query sequence, followed by a blastp or psi-blast search and the aligning of the retrieved sequences, or from a user-provided multiple sequence alignment. next, position-specific sequence- and disorder conservation scores are calculated and finally these pairwise scores are combined to create the discons profile. the result of one such classification is displayed at the bottom of the figure. in the next step, the msa is used to construct an aligned disorder profile by running iupred, vsl, espritz or foldindex on each of the aligned sequences. the disorder scores are first transformed to a binary scale for each sequence. details about the transformation procedure for the three different disorder predictors are given in the help section of the website. next the fraction of disordered residues at every position in the msa across the different sequences is calculated, thereby effectively quantifying the position-wise conservation of disorder. next, sequence conservation scores for each position in the alignment are calculated using the algorithm developed by capra et al.. in ‘advanced’ and ‘from alignment’ modes, a number of parameters such as the algorithm of choice or the background distribution can be specified to make the calculation more robust. finally, positions in the msa are scored by both their sequence- and disorder conservation, ranging from to. based on these pairs of scores, each position will fall into one of four distinct categories. as suggested by bellay et al., positions with a higher degree of disorder can be ‘constrained’, if both the sequence and disorder conservation scores are or greater; ‘flexible’, if the sequence conservation is lower than but the disorder conservation is or greater; or ‘non-conserved’, if the disorder conservation is lower than but higher than positions with a disorder conservation score of are completely lacking disorder and therefore are considered as ‘structured’. thus, regions of constrained disorder show strong conservation both at the amino acid level, and also of the disorder feature, while flexible disordered regions are variable in terms of amino acid sequence, but retain a significant level of disorder in evolution. lastly, non-conserved disordered regions lack disorder as a conserved feature, and are generally thought not to be associated with functions. on the results page, the position-specific conservation profile is provided on the output screen, and the fractions of residues falling into each of these distinct categories are also displayed in a tabular format at the bottom of the page, effectively quantifying the conservation of disorder in the query sequence. the sequences of consecutive ‘constrained’ disordered regions are also recorded and are available for download in fasta format along with the profiles and fractions in text format using the links that are provided. such segments of consecutive stretches of constrained disorder are most likely to correspond to functionally important idrs such as linear motifs or morfs, as we describe below.  bench-marking on molecular recognition features we evaluated the performance of discons on a set of molecular recognition features, or morfs, that are short peptide segments mostly found within longer disordered regions and involved in the binding to protein partners via disorder-to-order transition. morfs have been implicated in functions involving regulation and signaling, among other cellular processes. these recognition features are enriched in disordered residues, however, they may also have some residual structure, and their sequences are relatively more conserved than their flanking disordered regions. to estimate the efficiency of the discons protocol in distinguishing between such functionally important disordered segments on a large scale, we retrieved three morf datasets from morfpred that are available at their website and combined them into a single benchmarking dataset. the three datasets were the ‘test dataset’ containing morfs deposited in the protein data bank before 2008; the ‘experimental dataset’ with morfs identified between and 2012; and the ‘test 2012’ dataset with morfs from the combined dataset contained morf instances. after applying a sequence redundancy filter on the full length sequences using cd-hit, unique sequences remained. morf sequences were extracted from the full length protein sequences along with up to residue long flanking segments on both sides. disorder propensity scores were then calculated for the extracted morfs, the flanking residues, the full length proteins, and the complete uniprot/swissprot database, using iupred. figure 2a displays the distribution of disorder scores, comparing these four datasets, demonstrating that morf residues and morf-flanking regions are more disordered than morf-containing proteins. the difference is even more pronounced when compared to proteins from the complete uniprot/swissprot database. since the distributions did not follow gaussian or even symmetric distribution, we chose the non-parametric kolmogorov-smirnov test, which only assumes that the compared variables are continuous. the distributions of morfs and flanking residues are significantly different from both their proteins and the uniprot/swissprot dataset, according to ks tests with p-values less than the precision limit of r, indicating very strong significance.figure analyzing morfs with discons. morfs are known to be enriched in disordered residues. panel a shows normalized density distributions of morf regions, morf-flanking segments, full-length sequences of the morf-containing proteins and all proteins from swiss-prot. in these distributions, the area under each curve adds up to one. morf residues are predicted to be more disordered on average than the proteins they are found in, and especially more than the proteins of uniprot/swissprot. the flanking regions of the morfs also have significantly higher disorder content, compared to the full-length proteins they are found in. panel b shows the combined disorder- and sequence conservation scores, which range from to. by comparing these residue-specific score pairs, discons further supports the idea that the sequences of morfs are more conserved than that of their flanking regions; however, even in the flanking regions intrinsic disorder as a feature is highly conserved, indicating that these segment are required to be flexible in order for the protein segment to function. after performing psi-blast searches against each morf and their flanking regions, the conservation of the aligned positions were quantified in terms of sequence- and disorder conservation according to the protocol explained in the ‘implementation’ section. the binned conservation score pairs are displayed on figure 2b, comparing morfs and flanking residues. in order to statistically compare the conservation scores, we used the welch t-test, since this test only assumes gaussian distribution of the variables, and does not require equal variances. the conservation of disorder in morf and flanking residues is similar, while the underlying amino acid sequence is significantly more conserved in the case of morf residues. therefore, the comparison of morfs and morf-flanking sequences shows a trend that is in agreement with the literature, namely that while both sequence and disorder are rather conserved in the morfs, their neighboring protein segments are less conserved sequence-wise. bench-marking on short linear motifs following the analysis of the morf dataset, we applied the discons procedure on all the known instances of short linear motifs of the elm database. generally, these motifs are enriched in disorder, and their sequences show higher than average conservation on the amino acid level. we compared these motifs to all the available idp sequences retrieved from disprot. the average disorder content of slims is significantly higher than that of the full length idps, and in fact even more so than in the morf dataset. as expected and demonstrated in figure 3b, both the sequence and disorder conservation scores of the slim sites are significantly higher than that of the full length idps, with 52% of all the slim residues being of constrained disorder compared to only 12% in the full length idps, indicating the importance of structural disorder in slims. in comparison to morfs, where 36% of the residues are of constrained disorder, structural flexibility seems to play a larger role in slims, and indeed, morfs are known to often have some residual pre-formed structural elements, while slims are more disordered overall. additionally, the majority of slims are localized on consecutive segments of constrained disorder. concretely, 88% of the slims were found within consecutive segments of or more constrained residues, % in segments that are at least residues long, and % in segments of at least constrained residues in length. since slim segments often form ligand/protein binding sites, their high sequence conservation is necessitated by the formation of interface contacts with partner proteins.figure analyzing slims with discons. short linear motifs are peptide motifs with characteristic sequence patterns and are generally enriched in disorder. as displayed on panel a with the help of normalized density curves, slims are significantly more disordered than the full length proteins they are found in or compared to the proteins found in disprot. not only are slims more disordered, but this flexibility is highly conserved across homologs as well. when compared to the idps of disprot the difference in conservation is striking. as in figure 2b, disorder- and sequence conservation scores range from to. case studies of discons uncovering constrained and flexible disorder finally, we provide two case studies using two different protein segments; one exemplifying “constrained” disorder and the other being an example of the “flexible” disorder class. figure 4a shows a morf region of constrained disorderfound in the c-terminal negative regulatory domain of the p protein bound to the s calcium-binding protein, along with the discons profile of the c-terminal part of the p sequence. all the residues forming the morf are constrained based on the conservation profile of sequence- and of structural disorder, and, not surprisingly, these are the only residues of the disordered segment that appear in the crystal structure. the interaction with s restricts access to phosphorylation and acetylation sites on p that are important for transcription activation. thus, this region that is important for mediating a critical interaction is clearly identified by our protocol as a conserved disordered segment.figure example of a constrained disordered morf. the intrinsically disordered c-terminus of p adopts a helical fold upon partner binding. this is the only segment of the disordered region that appears in this crystal structure. upon quantifying the conservation of sequence and of disorder in the c-terminus, the full length morf segment is classified as ‘constrained’ disordered, where both features are highly conserved. since ‘flexible’ disordered segments generally function as entropic chains linking structured segments, finding structural data for them is less straightforward. these residues are often missing from the structures found in the protein data bank, however, ensemble descriptions of such regions are available from the protein ensemble database. in order to demonstrate “flexible” disorder, we retrieved the “fuzzy” complex formed between sic and the cdc subunit of an scf ubiquitin ligase. as seen generally in “fuzzy” complexes, sic is a fully disordered protein that remains disordered even when bound to its partner. sic has multiple binding segments along its disordered chain that compete for binding to the same pocket on the receptor protein. figures 5a and b display two conformations present in the ensemble description of the complex, while figure 5c provides the discons profile of the corresponding segment of the sic protein. the profile clearly shows that a “flexible” disordered linker connects the two “constrained” disordered binding regions. in several of the different conformations that constitute the ensemble of the “fuzzy” complex, both binding regions are found to contact the binding pocket, whereas the linker segment does not bind to cdc in any of the conformations. this indicates that the ‘unstructured-ness’ of the linker is more important for the function of the protein than the corresponding amino acid sequence of this segment.figure example of a flexible disordered linker. the intrinsically disordered protein sic binds to the cdc subunit of an scf ubiquitin ligase forming a fuzzy complex. two conformations of the same sic fragment from the ensemble of this complex are displayed on panels a and b. below, the discons profile on panel c corresponds to the same fragment, and shows the segment that is flexible disordered between two constrained disordered regions that bind to cdc this former segment corresponds to a flexible linker, which is found between two morfs that bind to the same binding pocket in the fuzzy complex. the flexible conservation of the linker indicates that the sequence is of less importance than its conformational flexibility.  discons is a novel and freely available online and downloadable tool that combines the quantitative description of the position-specific evolutionary conservation of the amino acid sequence with predictions of the conservation of its disordered/flexible state, providing meaningful information on the evolutionary context of a disordered protein segment. furthermore, discons uses this combined information to classify each disordered position into one of three categories, namely: constrained, flexible and non-conserved. these classes have been suggested to correlate with distinct functions that arise from the disordered state; therefore discons may provide information orthogonal to those obtained by other methods, which potentially enhances the reliability of the identification of functionally relevant disordered segments within proteins. we demonstrated that discons can be used to investigate both sequence- and disorder conservation in a functionally meaningful manner by bench-marking our procedure on morf and slim datasets, which are known to be conserved functional units enriched in structural disorder. it is important to emphasize that the success of calculation with discons strongly depends on the quality of the underlying multiple sequence alignment; therefore it is advised to review and optimize each msa to maximize the information of the output. taken this into consideration, discons can be used as an online or stand-alone tool for quantifying the conservation of both sequence and structural disorder by analyzing large-scale protein datasets and individual proteins. as such, discons might provide an additional layer of information for the investigation of protein disorder, and could serve to enhance the performance of prediction software such as morfpred, or provide descriptive information for disorder related databases such as d2p, mobidb or ped. availability and requirements discons is available as a web application, and as source code, both hosted at http://pedb.vib.be/discons. the source code is written in python, and has two versions: the multiple sequence alignment -based script, and the complete pipeline. the multiple alignment-based version has no requirements, while deployment of the full pipeline locally requires the following software: blast+ and mafft. the software is distributed under the gnu gpl license. abbrevations idpintrinsically disordered protein ldrlong disordered region morfmolecular recognition feature msamultiple sequence alignment slimshort linear motif competing interests the authors declare that they have no competing interests. authors' contributions vm conceived the project and developed the web application of discons. vm developed the underlying source code. vm, fzs and mg tested the code and the web application. vm and mg carried out the bench-marking analyses. vm, mg, fzs and pt wrote the paper. all authors read and approved the final manuscript. 
discons is an openly accessible sequence analysis tool that identifies and highlights structurally disordered segments of proteins where the conformational flexibility is conserved across homologs, and therefore potentially functional. however, in the case of many idps, the lack of evolutionary conservation of the primary sequence can hamper the study of functionality, because the conservation of their disorder profile and ensuing function may not appear in a traditional analysis of the evolutionary history of the protein. discons therefore helps with elucidating the function arising from the disordered state by analyzing individual proteins as well as large-scale proteomics datasets. here we present discons, a novel pipelined tool that combines the quantification of sequence- and disorder conservation to classify disordered residue positions. the tool is freely available both as a web application and as stand-alone source code hosted at http://pedb.vib.be/discons. analyzing the amino acid sequence of an intrinsically disordered protein in an evolutionary context can yield novel insights on the functional role of disordered regions and sequence element.
 analyzing the amino acid sequence of an intrinsically disordered protein in an evolutionary context can yield novel insights on the functional role of disordered regions and sequence element. however, in the case of many idps, the lack of evolutionary conservation of the primary sequence can hamper the study of functionality, because the conservation of their disorder profile and ensuing function may not appear in a traditional analysis of the evolutionary history of the protein.  here we present discons, a novel pipelined tool that combines the quantification of sequence- and disorder conservation to classify disordered residue positions. according to this scheme, the most interesting categories are constrained disordered residues and flexible disordered residues. the former residues show conservation of both the sequence and the property of disorder and are associated mainly with specific binding functionalities, whereas the latter class correspond to segments where disorder as a feature is important for function as opposed to the identity of the underlying sequence. discons therefore helps with elucidating the function arising from the disordered state by analyzing individual proteins as well as large-scale proteomics datasets.  discons is an openly accessible sequence analysis tool that identifies and highlights structurally disordered segments of proteins where the conformational flexibility is conserved across homologs, and therefore potentially functional. the tool is freely available both as a web application and as stand-alone source code hosted at http://pedb.vib.be/discons. keywords intrinsic protein disorderlarge-scale sequence analysismolecular recognition features short linear motifs issue-copyright-statement© the author 2015 