BACKGROUND
virus-host protein-protein interactions  are essential to understand viral attack mechanisms.
vhppis are used by viruses to disrupt or modulate host pathways in order to achieve goals like the evasion of the complement system  <cit> , modulation of the cytokine system  <cit>  and abrogation of apoptosis  <cit> .
some of these ppis are based on mimicry: a viral protein mimicking a host protein might interact with the host protein binding partners.
the mimicry is achieved through protein sequence or structural similarity  <cit> .
we focus our study on predicting a subset of ppis, the ones mediated by mimicked short linear motifs .
slim-mediated ppi predictions, conveniently ranked, might help researchers to postulate hypothesis to elucidate viral attack mechanisms, design antivirals and vaccines .
a short linear motif   is a short region of a protein,  <dig> to  <dig> residues long, with functions like controlling the assembly of protein complexes, marking proteolytic cleavage, tagging protein localization and enzyme recruiting  <cit> .
slims are structurally compact and participate in transitory low-affinity interactions  <cit> .
slims in eukaryotic proteins are curated in the elm database  <cit> .
slims might evolve rapidly in viral disordered regions through insertions, deletions and mutations  <cit> .
the new slims can change the ppi networks creating new advantageous ppis that can alter the cell cycle  <cit> , form protein complexes and mediate conformational changes  <cit> .
a recent analysis of the experimentally inferred human-virus ppis concludes that human proteins interacting with viruses are enriched in slims and binding interfaces  <cit> .
viruses use vhppis mediated by host-mimicked slims to hijack cell regulation  <cit>  and execute their viral cycle  <cit> .
an example of this strategy is the set of slim-mediated interactions of human papilloma virus  protein e <dig> with members of the 14-3- <dig> protein family and proteins containing pdz domains  <cit> .
experimental determination of vhppis is expensive since the number of proteins for some host organisms is large, more than  <dig>  in humans.
there are many viral protein sequences available but few corresponding three-dimensional structures resolved to use structure-based interaction prediction methods.
these are reasons for developing a general method for predicting mimetic host-virus ppis based solely in sequence data.
a bioinformatic approach to predict slim-mediated vhppis might be an inexpensive alternative to experimentation or can guide experimental design.
slims are represented computationally as regular expressions.
a slim instance is a protein subsequence that matches the regular expression.
for instance, a slim represented by the regular expression r.r.
have several instances like rvrre in ebola virus  <cit>  and rkrrf in human respiratory syncytial virus a <dig>  <cit> .
an algorithm for predicting virus-host ppis consist in finding viral instances of slims located in host proteins.
the viral instances found need to be filtered by some criteria that increase the probability of inferring real interactions.
if a slim is conserved in a small viral genome it probably could be used to interact with a host protein.
evans et al.
find that common slims between hiv- <dig> and humans are significantly conserved in hiv- <dig> proteins  <cit> .
they propose a criterion to filter slims if they are conserved above a 70% in the available viral sequences.
viral genomes have high mutation rates and are not too thermodynamically stable.
this seems to favor protein structures with a small number of inter-residue interactions and a high number of polar residues that account for the abundance of disordered protein regions  <cit> .
slims occur more frequently in viral protein disordered regions  <cit> , in different amounts between viral families  <cit> .
viral hubs, proteins that have many interactions with host proteins, tend to have more disordered regions  <cit> .
with these antecedents hagai et al.
propose a criterion to filter slims based on location in protein disordered regions  <cit> .
hagai et al.
also propose another criterion to filter slims based on rarity in a big set of randomized proteins  <cit> .
a slim is judged as rare, or hard to form by pure chance, if it is counted in less than a percentage of the sequences in the set of randomized proteins, e.g.
1% of the sequences.
hagai et al.
find that rare slims located in disordered regions have a significant enrichment in functional slims i.e.
with experimental evidence for interaction with host proteins  <cit> .
to our knowledge, there is no comparison of slim filtering methods in the literature.
for that reason, we implement and compare the three criteria introduced above for slim filtering: conservation above a threshold of the available viral sequences, localization in a protein disordered region and rarity, or difficulty to form by pure chance.
each filtering method produces a set of slims – conserved , disordered  and rare .
with sets c,d,r we form derived union and intersection sets: c∪d, c∪r, d∪r, c∪d∪r and c∩d, c∩r, d∩r, c∩r∩d.
each of these sets allow us to predict interactions between the viral protein containing the slim and the host proteins that interact with the slim.
all the sets generated are compared by filtering strength.
they also are compared by the number of vhppis derived from the set that have supporting evidence in a database –i.e.
correctly predicted.
the comparison by number of vhppis correctly predicted by set allow us to rank the vhppis partially.
to conduct the comparison of the sets we choose the human immunodeficiency virus .
it is the virus with more bioinformatic data available, with the niaid databases for sequences and alignments  <cit>  and for interactions with human proteins  <cit> .
we also use the hiv-1-human ppis mediated by slims as reported in the lmpid database  <cit> .
methods
disorder prediction
protein preprocessing
we download alignments for hiv- <dig> proteins env, gag, nef, pol, rev, vif, vpr, tat, vpu for the year  <dig> and an alignment of gag-pol dna sequences with years previous to  <dig> from the niaid hiv- <dig> sequence database  <cit> .
gag-pol sequences were translated following reference  <cit> : of  <dig> sequences,  <dig> containing the slippery subsequence tttttta were used to perform a computational translation considering the frameshift at the given subsequence.
we filter all protein sequences by hiv- <dig> subtypes b and c for their worldwide dominance and computationally cleave some of the alignments in the following manner: env into gp <dig>  gp <dig>  pol into pr, rt, rtp <dig>  in and gag into ma, ca, p <dig>  nc, p <dig>  p <dig>  <cit> .
after the cleavage we eliminate the gaps and asterisks in the resulting alignments in order to reinterpret the files as sets of sequences, fig.
<dig>  disorder panel.
the number of sequences per hiv- <dig> protein is in additional file 1: table s <dig>
fig.
<dig> methodology.
the methods are divided in three parts: 1) disorder: sequence preprocessing and prediction of disordered regions, 2) counting: counting of slim patterns and instances, and 3) comparison: analysis of the overlap between predicted interactions against interactions in niaid and lmpid databases
protein disorder prediction with iupred
among several disorder prediction algorithms for proteins  <cit>  we use iupred  <cit> .
this predictor implements a physical model based on force fields between residues statistically calibrated with a set of globular proteins in pdb  <cit> .
its performance is comparable to other predictors  <cit>  and can be installed locally.
iupred is enhanced with a sliding window addition proposed by hagai et al.
that allows to define disordered regions  <cit> .
residues with iupred computed values higher that  <dig>  are considered disordered.
for each residue an average disorder value is computed considering the iupred values for surrounding residues in a window of size  <dig>  this averaging is justified because the disorder tendency of the neighbors of a residue influence its disorder tendency.
residue windows with average disorder value higher than  <dig>  are considered as disordered.
as iupred receives as input a fasta file with only one sequence, we split fasta files with multiple sequences, call iupred on every split sequence-file, compute the sliding window based average values and give as output a list of disordered regions per protein sequence id.
we set the parameter long when calling iupred, see fig.
<dig>  disorder panel.
protein randomization
we randomize the hiv- <dig> proteins to create a big data set.
for each sequence in a protein file we create  <dig> shuffled versions randomizing the residues located in disordered regions of the sequence, as computed with iupred.
all disordered residues in a protein are joined together in a temporary list, shuffled with the modern fisher–yates algorithm, and put back in the disordered regions, leaving the ordered residues intact.
slim counting
we download all the slims, instances and interactions from the elm database  <cit>  and create an in-memory elm data structure with each slim identifier, its regular expression, its instances and its interactions with protein domains.
we wrote scripts to compute: the number of sequences with a given slim, the number of slim instances per protein, the number of slims conserved above a percentage of sequences  and the number of slims in disordered regions .
after randomizing as described above, we count the rare  slims in these shuffled data set, i.e.
the slims that are found in 1% of the randomized sequences or less .
based on c, d, r we create the union sets c∪d, c∪r, d∪r, c∪d∪r and intersection sets c∩d, c∩r, d∩r, c∩r∩d.
see fig.
<dig>  panels counting and comparison.
prediction of protein-protein interactions
we download the niaid human–hiv- <dig> ppi database  <cit> .
as the proteins in the database are identified by refseq records and the slim-domain interactions given by the elm database are given by uniprot records, we map refseq to uniprot identifiers for human proteins using uniprot id mappings.
we also download the lmpid database that curates virus-host elm-mediated interactions  <cit> .
for each slim set  obtained per hiv- <dig> protein we create vhppis based on the elm database interactions and interacting domains.
for each interaction reported in elm we add the human protein interacting with the slim located in the viral protein.
we also add the proteins that contain the domains listed as interacting with the slim.
to map domains to human proteins we used the domain-protein mapping for the human proteome in the pfam ftp server  <cit> .
figure  <dig>  comparison panel.
comparison of filtering methods
to validate a prediction we use two sets: the niaid hiv-1-human interactions and the set of elm mediated hiv-1-human interactions, as identified in lmpid  <cit> .
we count the number of correctly predicted interactions, when an interaction deduced with one of the slim sets is in the niaid database.
for all the slim sets obtained, and all the hiv- <dig> proteins, we analyze the overlap between the set of predicted human proteins interacting with hiv- <dig> and the set human proteins in niaid interactions.
we compute p-values for this overlap using the hyper-geometric distribution from the scipy python library, table  <dig>  the total number of human proteins was estimated as  <dig>  from reference  <cit> .p-values for the overlap between predicted interactions and niaid ppis
the p-value indicates the probability that the overlap between our sets of predicted ppis and the ppis with literature support in the niaid database takes place under the null hypothesis, that our sets were formed by random sampling.
red values are not significant at a level of  <dig>
RESULTS
a general method to identify slim-mediated ppis in eukaryotes
as slims are computationally represented by regular expressions there is always a possibility of finding instances in viral sequences by pure chance.
for this reason, it is important to develop slim filtering methods.
three filtering methods are implemented and systematically compared: conservation, location in disordered regions and rarity.
the combination of filters produces a method to predict virus-host ppis and rank them.
the comparison of filtering methods performance is conducted with the virus with more abundant data, hiv- <dig>  in fig.
<dig> there is an overview of the methods used.
the developed method only use protein sequences as input and do not depend on protein 3d structures, for this reason it can be used with any sequenced eukaryotic virus to infer candidate vhppis.
the restriction to eukaryotic viruses is based on the higher number of slims in eukaryotes and the use of the elm database, because the elm slim classes mod  and trg  are less used in prokaryotes  <cit> .
candidate interactions
the lists of predicted human-hiv- <dig> interactions that are not in the niaid database are in the .
disordered regions and slims in hiv- <dig> proteins
the disordered regions for hiv- <dig> proteins are in the .
they are depicted in fig.
<dig>  subfigures a to u show the predicted disordered content in hiv- <dig> proteins and polyproteins.
each protein sequence is represented as a yellow line and disordered regions are depicted as red segments.
fig.
<dig> disordered regions for hiv- <dig> proteins.
each subfigure from a to u contains an hiv- <dig> protein or protein precursor.
for all subfigures, each yellow line represents a protein sequence.
the red segments represent disordered regions as deduced with iupred with the sliding window addition explained in the “methods” section
we find that predicted disordered regions for hiv- <dig> proteins are relatively conserved.
perhaps the virus must keep flexibility in their proteins in order to interact with several partners.
in fig.
<dig>  we show the percentage of slims conserved above a 70% of the input sequences that are also located inside a disordered region.
most of the conserved slims in hiv- <dig> are located in protein disordered regions.
fig.
<dig> percentage of conserved slims that are located in disordered regions in hiv- <dig>  we plot the percentage of conserved slims, present in 70% or more of the input sequences, that are localized in a predicted disordered region
the proteins that deviate the most from this tendency are vpr, vpu, gp <dig>  in, and pr, with a percentage of conserved motifs that are located in disordered regions of  <dig> ,  <dig> ,  <dig> ,  <dig>  and 0% respectively.
the reason for this discrepancy lies in the few disordered regions predicted in the five proteins.
indeed, pr, in and gp <dig> are considered mostly ordered, while vpr and vpu are considered moderately disordered  <cit> .
a similar correlation between evolutionary conservation and location in disordered regions was found for the slims that bind to sh <dig>  sh <dig> and ser/thr kinase domains  <cit> .
we use iupred as disorder predictor only because its performance finding the disordered regions of the vif protein is outstanding compared to other  <dig> disorder predictors  <cit> .
one procedure that could be used to avoid structured regions entirely is a blast query against hiv- <dig> proteins in the protein data bank excluding hit regions.
however, it seems that disorder is a viral strategy to buffer mutations and increase interactions with host proteins  <cit> .
in this perspective, small disordered regions could be located inside structured protein regions to allow some interactions with the host, and not excluding the structured regions opens the possibility of finding these regions.
analysis of slim sets obtained
a ranking of slim sets by filtering strength
the slim set sizes are in additional file 3: table s <dig> and the slim sets for hiv- <dig> proteins are in the .
in fig.
<dig> we plot the the number of slim regular expressions that were found in the hiv- <dig> proteins identified by set.
the intersection slim sets c∩r  and c∩d∩r  were discarded for being almost empty for all proteins.
fig.
<dig> number of slims by set.
we plot the number of slims  that were found in hiv- <dig> proteins.
the intersection slim sets c∩r  and c∩d∩r  were discarded for being almost null in all entries
considering the sizes of slim sets we can rank them by the filtering strength; from low to high filtering.
the obtained ranking is r,d,c,c∩d,d∩r.
the criterion that filters the most is location in a disordered region and rarity.
it is followed by location in a disordered region and conservation.
the sets d∩r  studied by hagai et al.
<cit> , tend to have a smaller size than sets c∩d, of slims conserved and located in protein disordered regions, fig.
<dig>  the intersection slim sets c∩r  and c∩d∩r  are almost empty so they can be discarded as useful filtering criteria –data in additional file 3: table s <dig>
protein-protein interactions predicted with the slim sets are enriched in experimentally validated hiv-1-human protein-protein interactions
we validate against two virus-host ppis databases: niaid  <cit>  and lmpid  <cit> .
the niaid contains  <dig> ppis at the moment of writing while lmpid contains  <dig> ppis between several viruses and hosts, with  <dig> ppis between hiv- <dig> and human proteins.
the validation of the predicted ppis with the niaid database is not the best way to gauge the proportion of slim-based interactions.
this database contains ppis of all kinds, not only slim-mediated ones.
however, it is the most complete virus-host ppi dataset.
a better validation set, conceptually, is constructed with pairs deemed to interact through a slim with the lmpid database.
nevertheless, this dataset is too small.
we do the comparison with both databases, selecting the niaid database to compare the sets prediction performance and check the statistical significance of the results.
although we are suggesting a partial ranking of slim-based predicted ppis, another addition would be to rank totally the interactions with a score representing the probability that the interaction takes place based on experimental data  <cit>  or other techniques  <cit> .
for the moment, a total ranking is difficult to achieve given the scarcity of data about slim-mediated ppis  <cit> .
in the niaid database
in fig.
<dig> we plot the percentage of correctly predicted interactions, i.e.
stored in the niaid database and predicted with base on our slim sets.
in fig.
<dig> we plot the number of interactions predicted against the total number of interactions in the niaid database per hiv- <dig> protein.
the number of correctly predicted interactions is in additional file 4: table s <dig> and the number of novel interactions found with our method is in additional file 4: table s <dig>
fig.
<dig> percentage of validated interactions per slim set.
each subfigure plots the percentage of validated interactions with a slim set.
interactions were validated with the niaid hiv- <dig> human interaction database.
the percentage of predicted interactions is represented with a magenta bar.
the hiv- <dig> proteins are sorted by percentage in each subfigure
fig.
<dig> number of predicted interactions per slim set.
the number of predicted interactions by slim set, contrasted with the number of interactions in the niaid hiv- <dig> human interaction database.
the number of predicted interactions is represented with a red bar.
the green bar represents the total number of interactions in the niaid database
we use the hyper-geometric distribution to measure the statistical significance of the sets of interactions we found.
the p-values for the overlap between the ppis predicted with base on each slim set and the ppis in the niaid database are in table  <dig>  the sensitivity and specificity for the slim sets as ppi predictors is in additional file 4: table s <dig> and additional file 4: table s <dig>
in the lmpid database
using the literature curated lmpid database  <cit> , we find that the motif sets c,c∩d,c∩r,c∩d∩r capture half of the interactions in lmpid, while the sets c∪d,c∪r,d∪r,c∪d∪r,d∪r allow to infer all of them.
all the interactions between hiv- <dig> and human extracted from lmpid are in additional file 4: table s <dig>
the small number of human-hiv- <dig> interactions in this database , leaves open two possibilities: the number is really small, or the number is larger but few experiments have been performed to detect them.
to estimate the number of human-hiv- <dig> slim-mediated ppis more work is needed, perhaps an approach based on combining expert opinions  <cit> .
ppis correctly predicted serve as a ranking of filtering methods
in fig.
<dig> we plot the number of predicted interactions correctly validated against the niaid database identified by the slim set used to infer them.
we find that the slim sets have an almost general tendency with respect to the number of ppis correctly predicted across all hiv- <dig> proteins.
for this reason we propose to rank the ppis predicted according to the set used to deduce them.
fig.
<dig> number of validated interactions per set.
the number of interactions as validated with the niaid hiv- <dig> human interaction database.
each set is represented by a dashed line with a different color
the ranking of the sets we found by its capacity to infer real interactions is: d∪r,c∪r,c∪d,c∩d,d∩r.
this ranking allow to present the ppis predicted to researchers in a partial order: first the set of interactions deduced with d∪r –slims located in disordered regions or hard to form by pure chance, then the set deduced with c –conserved slims.
most used slims in hiv- <dig> proteins suggest hiv- <dig> extensive use of human protein signaling and other processes
we consider the set c∩d of slims conserved and located in disordered regions for their biological relevance to analyze their human counter domains.
in table  <dig> we include the most used slims from this set, i.e the slims that are present in  <dig> or more hiv- <dig> proteins, are conserved, and localize in a disordered region.
in general, most of these slims would interfere with host signaling.
the most used counter domains are the protein kinase domain  that interact with  <dig> of the most used slims and the peptidase_s <dig> subtilase family  that interacts with  <dig> of the most used cleavage slims.
however, the list of counter domains in table  <dig> suggest that hiv- <dig> slims interfere with transcription regulation, autophagy, cell cycle control, apoptosis and cellular transport.
from the slims that are conserved in more than 70% of the hiv- <dig> protein sequences, and are located in disordered regions too we counted how many hiv- <dig> proteins include them.
in this table we report the slims more commonly used in hiv- <dig> proteins, the ones that are included in  <dig> or more of the hiv- <dig> proteins.
the table includes the counter domain for every slim
CONCLUSIONS
we develop a method to predict virus-host slim-mediated ppis and rank them.
it is applicable to any eukaryotic virus and host with available protein sequences.
using data for the most studied virus, hiv- <dig>  we find a partial ordering of the ppis obtained based on the set used to infer the interactions.
this order is descending in the expected probability of inferring real interactions.
we expect that the method gives interesting candidate interactions with other eukaryotic viruses and hosts.
the call for using high-throughput methods to detect slim-mediated ppis illustrates the benefits of a bioinformatic method that predicts slim-mediated ppis and might guide experimental design  <cit> .although the number of slim-mediated ppis might be small, there is evidence that these ppis are used by several viruses, in contrast to virus-host domain-domain ppis, that are virus-specific  <cit> .
this kind of interactions can help to analyze common viral strategies for infection.
indeed, in a previous work we used the method with the viruses in the ncbi virus variation resource to predict interactions with the proteins from the human protein synthesis machinery  <cit> .
we found evidence that viruses interact with eukaryotic initiation factors  <dig> and  <dig>  and the poly-binding proteins using slims.
even though the method developed is not a strong predictor, by using several viruses, interesting interactions with host subsystems can be uncovered.
in a future work we want to scale the approach considering all the human proteome and more human viruses.
in future work we could also incorporate structural information in the prediction and analysis of slim-mediated vhppis in order to create otter slim filtering methods and compare them with the filters obtained in this work.
one possibility is the study of fuzziness and slim flanking regions  <cit> , another one is the use of disordered binding region prediction methods, like anchor  <cit> .
additional files
additional file  <dig>
table s <dig>  candidate interactions between human and hiv- <dig>  available at https://figshare.com/articles/interactions_zip/ <dig>
additional file  <dig>
table s <dig>  disordered regions for hiv- <dig> proteins , using filenames layout explained in table s <dig> available at https://figshare.com/articles/disordered_regions_predicted_in_hiv-1/ <dig>
additional file  <dig>
table s <dig>  slim sets c,d,r and derived for hiv- <dig> proteins.
available at https://figshare.com/articles/short_linear_motif_sets_common_to_human_and_hiv-1/ <dig>
additional file  <dig>
table s <dig>  supplement information describing the previous files.
available at https://figshare.com/articles/supplement_-_prediction_of_virus-host_protein-protein_interactions_based_on_short_linear_motifs/ <dig>
abbreviations
elmeukaryotic linear motifs
hiv-1human immunodeficiency virus - 1
lmpidlinear motif mediated protein interaction database
niaidnational institute of allergy and infectious diseases
pdbprotein data bank
pfamprotein families database
ppiprotein-protein interaction
slimshort linear motif
vhppivirus-host protein-protein interaction
the authors would like to thank dr. aydin tozeren for receiving ab in the biomedical engineering department at drexel university, and drs.
irene tischer and angel garcía for their valuable input and comments at initial stages of the research.
funding
this study was supported by a colciencias scholarship granted to ab and the escuela de ingeniería de sistemas y computación de la facultad de ingeniería de la universidad del valle.
availability of data and materials
the datasets generated during the current study are available in the figshare repository, with links in the additional files section.
authors’ contributions
this study was conceived by ab with significant input to intellectual content at all stages from pm and vb.
ab executed the model operations.
all authors worked on and approved the final manuscript.
competing interests
the authors declare that they have no competing interests.
consent for publication
not applicable.
ethics approval and consent to participate
not applicable.

 virus-host protein-protein interactions are essential to understand viral attack mechanisms. vhppis are used by viruses to disrupt or modulate host pathways in order to achieve goals like the evasion of the complement system, modulation of the cytokine system and abrogation of apoptosis. some of these ppis are based on mimicry: a viral protein mimicking a host protein might interact with the host protein binding partners. the mimicry is achieved through protein sequence or structural similarity. we focus our study on predicting a subset of ppis, the ones mediated by mimicked short linear motifs. slim-mediated ppi predictions, conveniently ranked, might help researchers to postulate hypothesis to elucidate viral attack mechanisms, design antivirals and vaccines. a short linear motif is a short region of a protein, to residues long, with functions like controlling the assembly of protein complexes, marking proteolytic cleavage, tagging protein localization and enzyme recruiting. slims are structurally compact and participate in transitory low-affinity interactions. slims in eukaryotic proteins are curated in the elm database. slims might evolve rapidly in viral disordered regions through insertions, deletions and mutations. the new slims can change the ppi networks creating new advantageous ppis that can alter the cell cycle, form protein complexes and mediate conformational changes. a recent analysis of the experimentally inferred human-virus ppis concludes that human proteins interacting with viruses are enriched in slims and binding interfaces. viruses use vhppis mediated by host-mimicked slims to hijack cell regulation and execute their viral cycle. an example of this strategy is the set of slim-mediated interactions of human papilloma virus protein e with members of the 14-3- protein family and proteins containing pdz domains. experimental determination of vhppis is expensive since the number of proteins for some host organisms is large, more than in humans. there are many viral protein sequences available but few corresponding three-dimensional structures resolved to use structure-based interaction prediction methods. these are reasons for developing a general method for predicting mimetic host-virus ppis based solely in sequence data. a bioinformatic approach to predict slim-mediated vhppis might be an inexpensive alternative to experimentation or can guide experimental design. slims are represented computationally as regular expressions. a slim instance is a protein subsequence that matches the regular expression. for instance, a slim represented by the regular expression r.r. have several instances like rvrre in ebola virus and rkrrf in human respiratory syncytial virus a. an algorithm for predicting virus-host ppis consist in finding viral instances of slims located in host proteins. the viral instances found need to be filtered by some criteria that increase the probability of inferring real interactions. if a slim is conserved in a small viral genome it probably could be used to interact with a host protein. evans et al. find that common slims between hiv- and humans are significantly conserved in hiv- proteins. they propose a criterion to filter slims if they are conserved above a 70% in the available viral sequences. viral genomes have high mutation rates and are not too thermodynamically stable. this seems to favor protein structures with a small number of inter-residue interactions and a high number of polar residues that account for the abundance of disordered protein regions. slims occur more frequently in viral protein disordered regions, in different amounts between viral families. viral hubs, proteins that have many interactions with host proteins, tend to have more disordered regions. with these antecedents hagai et al. propose a criterion to filter slims based on location in protein disordered regions. hagai et al. also propose another criterion to filter slims based on rarity in a big set of randomized proteins. a slim is judged as rare, or hard to form by pure chance, if it is counted in less than a percentage of the sequences in the set of randomized proteins, e.g. 1% of the sequences. hagai et al. find that rare slims located in disordered regions have a significant enrichment in functional slims i.e. with experimental evidence for interaction with host proteins. to our knowledge, there is no comparison of slim filtering methods in the literature. for that reason, we implement and compare the three criteria introduced above for slim filtering: conservation above a threshold of the available viral sequences, localization in a protein disordered region and rarity, or difficulty to form by pure chance. each filtering method produces a set of slims – conserved, disordered and rare. with sets c,d,r we form derived union and intersection sets: c∪d, c∪r, d∪r, c∪d∪r and c∩d, c∩r, d∩r, c∩r∩d. each of these sets allow us to predict interactions between the viral protein containing the slim and the host proteins that interact with the slim. all the sets generated are compared by filtering strength. they also are compared by the number of vhppis derived from the set that have supporting evidence in a database –i.e. correctly predicted. the comparison by number of vhppis correctly predicted by set allow us to rank the vhppis partially. to conduct the comparison of the sets we choose the human immunodeficiency virus. it is the virus with more bioinformatic data available, with the niaid databases for sequences and alignments and for interactions with human proteins. we also use the hiv-1-human ppis mediated by slims as reported in the lmpid database. methods disorder prediction protein preprocessing we download alignments for hiv- proteins env, gag, nef, pol, rev, vif, vpr, tat, vpu for the year and an alignment of gag-pol dna sequences with years previous to from the niaid hiv- sequence database. gag-pol sequences were translated following reference: of sequences, containing the slippery subsequence tttttta were used to perform a computational translation considering the frameshift at the given subsequence. we filter all protein sequences by hiv- subtypes b and c for their worldwide dominance and computationally cleave some of the alignments in the following manner: env into gp gp pol into pr, rt, rtp in and gag into ma, ca, p nc, p p. after the cleavage we eliminate the gaps and asterisks in the resulting alignments in order to reinterpret the files as sets of sequences, fig. disorder panel. the number of sequences per hiv- protein is in additional file 1: table s fig. methodology. the methods are divided in three parts: 1) disorder: sequence preprocessing and prediction of disordered regions, 2) counting: counting of slim patterns and instances, and 3) comparison: analysis of the overlap between predicted interactions against interactions in niaid and lmpid databases protein disorder prediction with iupred among several disorder prediction algorithms for proteins we use iupred. this predictor implements a physical model based on force fields between residues statistically calibrated with a set of globular proteins in pdb. its performance is comparable to other predictors and can be installed locally. iupred is enhanced with a sliding window addition proposed by hagai et al. that allows to define disordered regions. residues with iupred computed values higher that are considered disordered. for each residue an average disorder value is computed considering the iupred values for surrounding residues in a window of size this averaging is justified because the disorder tendency of the neighbors of a residue influence its disorder tendency. residue windows with average disorder value higher than are considered as disordered. as iupred receives as input a fasta file with only one sequence, we split fasta files with multiple sequences, call iupred on every split sequence-file, compute the sliding window based average values and give as output a list of disordered regions per protein sequence id. we set the parameter long when calling iupred, see fig. disorder panel. protein randomization we randomize the hiv- proteins to create a big data set. for each sequence in a protein file we create shuffled versions randomizing the residues located in disordered regions of the sequence, as computed with iupred. all disordered residues in a protein are joined together in a temporary list, shuffled with the modern fisher–yates algorithm, and put back in the disordered regions, leaving the ordered residues intact. slim counting we download all the slims, instances and interactions from the elm database and create an in-memory elm data structure with each slim identifier, its regular expression, its instances and its interactions with protein domains. we wrote scripts to compute: the number of sequences with a given slim, the number of slim instances per protein, the number of slims conserved above a percentage of sequences and the number of slims in disordered regions. after randomizing as described above, we count the rare slims in these shuffled data set, i.e. the slims that are found in 1% of the randomized sequences or less. based on c, d, r we create the union sets c∪d, c∪r, d∪r, c∪d∪r and intersection sets c∩d, c∩r, d∩r, c∩r∩d. see fig. panels counting and comparison. prediction of protein-protein interactions we download the niaid human–hiv- ppi database. as the proteins in the database are identified by refseq records and the slim-domain interactions given by the elm database are given by uniprot records, we map refseq to uniprot identifiers for human proteins using uniprot id mappings. we also download the lmpid database that curates virus-host elm-mediated interactions. for each slim set obtained per hiv- protein we create vhppis based on the elm database interactions and interacting domains. for each interaction reported in elm we add the human protein interacting with the slim located in the viral protein. we also add the proteins that contain the domains listed as interacting with the slim. to map domains to human proteins we used the domain-protein mapping for the human proteome in the pfam ftp server. figure comparison panel. comparison of filtering methods to validate a prediction we use two sets: the niaid hiv-1-human interactions and the set of elm mediated hiv-1-human interactions, as identified in lmpid. we count the number of correctly predicted interactions, when an interaction deduced with one of the slim sets is in the niaid database. for all the slim sets obtained, and all the hiv- proteins, we analyze the overlap between the set of predicted human proteins interacting with hiv- and the set human proteins in niaid interactions. we compute p-values for this overlap using the hyper-geometric distribution from the scipy python library, table the total number of human proteins was estimated as from reference.p-values for the overlap between predicted interactions and niaid ppis the p-value indicates the probability that the overlap between our sets of predicted ppis and the ppis with literature support in the niaid database takes place under the null hypothesis, that our sets were formed by random sampling. red values are not significant at a level of  a general method to identify slim-mediated ppis in eukaryotes as slims are computationally represented by regular expressions there is always a possibility of finding instances in viral sequences by pure chance. for this reason, it is important to develop slim filtering methods. three filtering methods are implemented and systematically compared: conservation, location in disordered regions and rarity. the combination of filters produces a method to predict virus-host ppis and rank them. the comparison of filtering methods performance is conducted with the virus with more abundant data, hiv- in fig. there is an overview of the methods used. the developed method only use protein sequences as input and do not depend on protein 3d structures, for this reason it can be used with any sequenced eukaryotic virus to infer candidate vhppis. the restriction to eukaryotic viruses is based on the higher number of slims in eukaryotes and the use of the elm database, because the elm slim classes mod and trg are less used in prokaryotes. candidate interactions the lists of predicted human-hiv- interactions that are not in the niaid database are in the. disordered regions and slims in hiv- proteins the disordered regions for hiv- proteins are in the. they are depicted in fig. subfigures a to u show the predicted disordered content in hiv- proteins and polyproteins. each protein sequence is represented as a yellow line and disordered regions are depicted as red segments. fig. disordered regions for hiv- proteins. each subfigure from a to u contains an hiv- protein or protein precursor. for all subfigures, each yellow line represents a protein sequence. the red segments represent disordered regions as deduced with iupred with the sliding window addition explained in the “methods” section we find that predicted disordered regions for hiv- proteins are relatively conserved. perhaps the virus must keep flexibility in their proteins in order to interact with several partners. in fig. we show the percentage of slims conserved above a 70% of the input sequences that are also located inside a disordered region. most of the conserved slims in hiv- are located in protein disordered regions. fig. percentage of conserved slims that are located in disordered regions in hiv- we plot the percentage of conserved slims, present in 70% or more of the input sequences, that are localized in a predicted disordered region the proteins that deviate the most from this tendency are vpr, vpu, gp in, and pr, with a percentage of conserved motifs that are located in disordered regions of,,, and 0% respectively. the reason for this discrepancy lies in the few disordered regions predicted in the five proteins. indeed, pr, in and gp are considered mostly ordered, while vpr and vpu are considered moderately disordered. a similar correlation between evolutionary conservation and location in disordered regions was found for the slims that bind to sh sh and ser/thr kinase domains. we use iupred as disorder predictor only because its performance finding the disordered regions of the vif protein is outstanding compared to other disorder predictors. one procedure that could be used to avoid structured regions entirely is a blast query against hiv- proteins in the protein data bank excluding hit regions. however, it seems that disorder is a viral strategy to buffer mutations and increase interactions with host proteins. in this perspective, small disordered regions could be located inside structured protein regions to allow some interactions with the host, and not excluding the structured regions opens the possibility of finding these regions. analysis of slim sets obtained a ranking of slim sets by filtering strength the slim set sizes are in additional file 3: table s and the slim sets for hiv- proteins are in the. in fig. we plot the the number of slim regular expressions that were found in the hiv- proteins identified by set. the intersection slim sets c∩r and c∩d∩r were discarded for being almost empty for all proteins. fig. number of slims by set. we plot the number of slims that were found in hiv- proteins. the intersection slim sets c∩r and c∩d∩r were discarded for being almost null in all entries considering the sizes of slim sets we can rank them by the filtering strength; from low to high filtering. the obtained ranking is r,d,c,c∩d,d∩r. the criterion that filters the most is location in a disordered region and rarity. it is followed by location in a disordered region and conservation. the sets d∩r studied by hagai et al., tend to have a smaller size than sets c∩d, of slims conserved and located in protein disordered regions, fig. the intersection slim sets c∩r and c∩d∩r are almost empty so they can be discarded as useful filtering criteria –data in additional file 3: table s protein-protein interactions predicted with the slim sets are enriched in experimentally validated hiv-1-human protein-protein interactions we validate against two virus-host ppis databases: niaid and lmpid. the niaid contains ppis at the moment of writing while lmpid contains ppis between several viruses and hosts, with ppis between hiv- and human proteins. the validation of the predicted ppis with the niaid database is not the best way to gauge the proportion of slim-based interactions. this database contains ppis of all kinds, not only slim-mediated ones. however, it is the most complete virus-host ppi dataset. a better validation set, conceptually, is constructed with pairs deemed to interact through a slim with the lmpid database. nevertheless, this dataset is too small. we do the comparison with both databases, selecting the niaid database to compare the sets prediction performance and check the statistical significance of the results. although we are suggesting a partial ranking of slim-based predicted ppis, another addition would be to rank totally the interactions with a score representing the probability that the interaction takes place based on experimental data or other techniques. for the moment, a total ranking is difficult to achieve given the scarcity of data about slim-mediated ppis. in the niaid database in fig. we plot the percentage of correctly predicted interactions, i.e. stored in the niaid database and predicted with base on our slim sets. in fig. we plot the number of interactions predicted against the total number of interactions in the niaid database per hiv- protein. the number of correctly predicted interactions is in additional file 4: table s and the number of novel interactions found with our method is in additional file 4: table s fig. percentage of validated interactions per slim set. each subfigure plots the percentage of validated interactions with a slim set. interactions were validated with the niaid hiv- human interaction database. the percentage of predicted interactions is represented with a magenta bar. the hiv- proteins are sorted by percentage in each subfigure fig. number of predicted interactions per slim set. the number of predicted interactions by slim set, contrasted with the number of interactions in the niaid hiv- human interaction database. the number of predicted interactions is represented with a red bar. the green bar represents the total number of interactions in the niaid database we use the hyper-geometric distribution to measure the statistical significance of the sets of interactions we found. the p-values for the overlap between the ppis predicted with base on each slim set and the ppis in the niaid database are in table the sensitivity and specificity for the slim sets as ppi predictors is in additional file 4: table s and additional file 4: table s in the lmpid database using the literature curated lmpid database, we find that the motif sets c,c∩d,c∩r,c∩d∩r capture half of the interactions in lmpid, while the sets c∪d,c∪r,d∪r,c∪d∪r,d∪r allow to infer all of them. all the interactions between hiv- and human extracted from lmpid are in additional file 4: table s the small number of human-hiv- interactions in this database, leaves open two possibilities: the number is really small, or the number is larger but few experiments have been performed to detect them. to estimate the number of human-hiv- slim-mediated ppis more work is needed, perhaps an approach based on combining expert opinions. ppis correctly predicted serve as a ranking of filtering methods in fig. we plot the number of predicted interactions correctly validated against the niaid database identified by the slim set used to infer them. we find that the slim sets have an almost general tendency with respect to the number of ppis correctly predicted across all hiv- proteins. for this reason we propose to rank the ppis predicted according to the set used to deduce them. fig. number of validated interactions per set. the number of interactions as validated with the niaid hiv- human interaction database. each set is represented by a dashed line with a different color the ranking of the sets we found by its capacity to infer real interactions is: d∪r,c∪r,c∪d,c∩d,d∩r. this ranking allow to present the ppis predicted to researchers in a partial order: first the set of interactions deduced with d∪r –slims located in disordered regions or hard to form by pure chance, then the set deduced with c –conserved slims. most used slims in hiv- proteins suggest hiv- extensive use of human protein signaling and other processes we consider the set c∩d of slims conserved and located in disordered regions for their biological relevance to analyze their human counter domains. in table we include the most used slims from this set, i.e the slims that are present in or more hiv- proteins, are conserved, and localize in a disordered region. in general, most of these slims would interfere with host signaling. the most used counter domains are the protein kinase domain that interact with of the most used slims and the peptidase_s subtilase family that interacts with of the most used cleavage slims. however, the list of counter domains in table suggest that hiv- slims interfere with transcription regulation, autophagy, cell cycle control, apoptosis and cellular transport. from the slims that are conserved in more than 70% of the hiv- protein sequences, and are located in disordered regions too we counted how many hiv- proteins include them. in this table we report the slims more commonly used in hiv- proteins, the ones that are included in or more of the hiv- proteins. the table includes the counter domain for every slim  we develop a method to predict virus-host slim-mediated ppis and rank them. it is applicable to any eukaryotic virus and host with available protein sequences. using data for the most studied virus, hiv- we find a partial ordering of the ppis obtained based on the set used to infer the interactions. this order is descending in the expected probability of inferring real interactions. we expect that the method gives interesting candidate interactions with other eukaryotic viruses and hosts. the call for using high-throughput methods to detect slim-mediated ppis illustrates the benefits of a bioinformatic method that predicts slim-mediated ppis and might guide experimental design.although the number of slim-mediated ppis might be small, there is evidence that these ppis are used by several viruses, in contrast to virus-host domain-domain ppis, that are virus-specific. this kind of interactions can help to analyze common viral strategies for infection. indeed, in a previous work we used the method with the viruses in the ncbi virus variation resource to predict interactions with the proteins from the human protein synthesis machinery. we found evidence that viruses interact with eukaryotic initiation factors and and the poly-binding proteins using slims. even though the method developed is not a strong predictor, by using several viruses, interesting interactions with host subsystems can be uncovered. in a future work we want to scale the approach considering all the human proteome and more human viruses. in future work we could also incorporate structural information in the prediction and analysis of slim-mediated vhppis in order to create otter slim filtering methods and compare them with the filters obtained in this work. one possibility is the study of fuzziness and slim flanking regions, another one is the use of disordered binding region prediction methods, like anchor. additional files additional file table s candidate interactions between human and hiv- available at https://figshare.com/articles/interactions_zip/ additional file table s disordered regions for hiv- proteins, using filenames layout explained in table s available at https://figshare.com/articles/disordered_regions_predicted_in_hiv-1/ additional file table s slim sets c,d,r and derived for hiv- proteins. available at https://figshare.com/articles/short_linear_motif_sets_common_to_human_and_hiv-1/ additional file table s supplement information describing the previous files. available at https://figshare.com/articles/supplement_-_prediction_of_virus-host_protein-protein_interactions_based_on_short_linear_motifs/ abbreviations elmeukaryotic linear motifs hiv-1human immunodeficiency virus - 1 lmpidlinear motif mediated protein interaction database niaidnational institute of allergy and infectious diseases pdbprotein data bank pfamprotein families database ppiprotein-protein interaction slimshort linear motif vhppivirus-host protein-protein interaction the authors would like to thank dr. aydin tozeren for receiving ab in the biomedical engineering department at drexel university, and drs. irene tischer and angel garcía for their valuable input and comments at initial stages of the research. funding this study was supported by a colciencias scholarship granted to ab and the escuela de ingeniería de sistemas y computación de la facultad de ingeniería de la universidad del valle. availability of data and materials the datasets generated during the current study are available in the figshare repository, with links in the additional files section. authors’ contributions this study was conceived by ab with significant input to intellectual content at all stages from pm and vb. ab executed the model operations. all authors worked on and approved the final manuscript. competing interests the authors declare that they have no competing interests. consent for publication not applicable. ethics approval and consent to participate not applicable.
with respect to hiv- the sets c∪r, d∪r, c∪d∪r infer all known interactions between hiv and human proteins mediated by linear motifs. we have developed a method for predicting protein-protein interactions mediated by linear motifs between hiv- and human proteins. given that viral linear motif instances of host motif regular expressions can be found by chance, it is necessary to develop filtering methods of functional linear motifs. short linear motifs in host organisms proteins can be mimicked by viruses to create protein-protein interactions that disable or control metabolic pathways. we conduct a systematic comparison of linear motifs filtering methods to develop a computational approach for predicting motif-mediated protein-protein interactions between human and the human immunodeficiency virus. the resulting sets are compared by the number of protein-protein interactions correctly inferred with them – with experimental validation. we can extend the software developed to any other eukaryotic virus and host in order to find and rank candidate interactions. we implemented three filtering methods to obtain linear motif sets: 1) conserved in viral proteins, 2) located in disordered regions and 3) rare or scarce in a set of randomized viral sequences. the method only use protein sequences as inputs. we found that the majority of conserved linear motifs in the virus are located in disordered regions.
 short linear motifs in host organisms proteins can be mimicked by viruses to create protein-protein interactions that disable or control metabolic pathways. given that viral linear motif instances of host motif regular expressions can be found by chance, it is necessary to develop filtering methods of functional linear motifs. we conduct a systematic comparison of linear motifs filtering methods to develop a computational approach for predicting motif-mediated protein-protein interactions between human and the human immunodeficiency virus.  we implemented three filtering methods to obtain linear motif sets: 1) conserved in viral proteins, 2) located in disordered regions and 3) rare or scarce in a set of randomized viral sequences. the sets c,d,r are united and intersected. the resulting sets are compared by the number of protein-protein interactions correctly inferred with them – with experimental validation. the comparison is done with hiv- sequences and interactions from the national institute of allergy and infectious diseases. the number of correctly inferred interactions allows to rank the interactions by the sets used to deduce them: d∪r and c. the ordering of the sets is descending on the probability of capturing functional interactions. with respect to hiv- the sets c∪r, d∪r, c∪d∪r infer all known interactions between hiv and human proteins mediated by linear motifs. we found that the majority of conserved linear motifs in the virus are located in disordered regions.  we have developed a method for predicting protein-protein interactions mediated by linear motifs between hiv- and human proteins. the method only use protein sequences as inputs. we can extend the software developed to any other eukaryotic virus and host in order to find and rank candidate interactions. in future works we will use it to explore possible viral attack mechanisms based on linear motif mimicry. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords virushosteukaryotesproteininteractionpredictionshortlinearmotifdisorderescuela de ingeniería de sistemas y computaciónissue-copyright-statement© the author 2017 