Content deleted Content added
m Open access bot: doi updated in citation with #oabot. |
|||
(219 intermediate revisions by more than 100 users not shown) | |||
Line 1:
{{Short description|Prediction by observation and computation}}'''Protein–protein interaction prediction''' is a field combining [[bioinformatics]] and [[structural biology]] in an attempt to identify and catalog physical interactions between pairs or groups of proteins. Understanding [[protein–protein interaction]]s is important for the investigation of intracellular signaling pathways, modelling of protein complex structures and for gaining insights into various biochemical processes.
''Experimentally'', physical interactions between pairs of proteins can be inferred from a variety of techniques, including yeast [[two-hybrid screening|two-hybrid]] systems, [[Protein-fragment complementation assay|protein-fragment complementation assays]] (PCA), affinity purification/[[mass spectrometry]], [[protein microarray]]s, fluorescence resonance energy transfer (FRET), and [[Microscale Thermophoresis]] (MST). Efforts to experimentally determine the [[interactome]] of numerous species are ongoing. Experimentally determined interactions usually provide the basis for ''computational methods'' to predict interactions, e.g. using [[Homology (biology)|homologous]] protein sequences across species. However, there are also methods that predict interactions ''de novo'', without prior knowledge of existing interactions.
== Methods ==
Proteins that interact are more likely to co-evolve,{{r|Dandekar}}{{r|Enright}}{{r|Marcotte}}{{r|Pazos}} therefore, it is possible to make inferences about interactions between pairs of proteins based on their phylogenetic distances. It has also been observed in some cases that pairs of interacting proteins have fused orthologues in other organisms. In addition, a number of bound protein complexes have been structurally solved and can be used to identify the residues that mediate the interaction so that similar motifs can be located in other organisms.
=== Phylogenetic profiling ===
[[File:Phylogenetic Profiling Method.png|thumb|''Figure A.''The phylogenetic profiles of four genes (A, B, C and D) are shown on the right. A '1' denotes presence of the gene in the genome and a '0' denotes absence. The two identical profiles of genes A and B are highlighted in yellow<ref name=":0" />.]]
''[[Phylogenetic profiling|The phylogenetic profile]] method'' is based on the hypothesis that if two or more proteins are concurrently present or absent across several genomes, then they are likely functionally related.<ref name=":0">{{Cite journal|last=Raman|first=Karthik|date=2010-02-15|title=Construction and analysis of protein–protein interaction networks|journal=Automated Experimentation|volume=2|issue=1|pages=2|doi=10.1186/1759-4499-2-2|issn=1759-4499|pmc=2834675|pmid=20334628 |doi-access=free }}</ref> ''Figure A'' illustrates a hypothetical situation in which proteins A and B are identified as functionally linked due to their identical phylogenetic profiles across 5 different genomes. The Joint Genome Institute provides an Integrated Microbial Genomes and Microbiomes database ([https://img.jgi.doe.gov JGI IMG]) that has a phylogenetic profiling tool for single genes and gene cassettes.
=== Prediction of co-evolved protein pairs based on similar phylogenetic trees===
It was observed that the phylogenetic trees of ligands and receptors were often more similar than due to random chance.{{r|Pazos}} This is likely because they faced similar selection pressures and co-evolved. This method{{r|Tan}} uses the phylogenetic trees of protein pairs to determine if interactions exist. To do this, homologs of the proteins of interest are found (using a sequence search tool such as [[BLAST (biotechnology)|BLAST]]) and multiple-sequence alignments are done (with alignment tools such as [[Clustal]]) to build distance matrices for each of the proteins of interest.{{r|Pazos}} The distance matrices should then be used to build phylogenetic trees. However, comparisons between phylogenetic trees are difficult, and current methods circumvent this by simply comparing distance matrices{{r|Pazos}}. The distance matrices of the proteins are used to calculate a correlation coefficient, in which a larger value corresponds to co-evolution. The benefit of comparing distance matrices instead of phylogenetic trees is that the results do not depend on the method of tree building that was used. The downside is that difference matrices are not perfect representations of phylogenetic trees, and inaccuracies may result from using such a shortcut.{{r|Pazos}} Another factor worthy of note is that there are background similarities between the phylogenetic trees of any protein, even ones that do not interact. If left unaccounted for, this could lead to a high false-positive rate. For this reason, certain methods construct a background tree using 16S rRNA sequences which they use as the canonical tree of life. The distance matrix constructed from this tree of life is then subtracted from the distance matrices of the proteins of interest.{{r|PazosRanea}} However, because RNA distance matrices and DNA distance matrices have different scale, presumably because RNA and DNA have different mutation rates, the RNA matrix needs to be rescaled before it can be subtracted from the DNA matrices.{{r|PazosRanea}} By using molecular clock proteins, the scaling coefficient for protein distance/RNA distance can be calculated.{{r|PazosRanea}} This coefficient is used to rescale the RNA matrix.
[[File:Alignment of the Rosette Stone protein Human Succinyl-CoA-Transferase with Acetate-CoA-Transferase subunits alpha and beta.png|thumb|''Figure B.'' The Human succinyl-CoA-Transferase enzyme is represented by the two joint blue and green bars at the top of the image. The alpha subunit of the Acetate-CoA-Transferase enzyme is homologous with the first half of the enzyme, represents by the blue bar. The beta subunit of the Acetate-CoA-Transferase enzyme is homologous with the second half of the enzyme, represents by the green bar. This mage was adapted from Uetz, P. & Pohl, E. (2018) ''Protein–Protein and Protein–DNA Interactions''. In: Wink, M. (ed.), Introduction to Molecular Biotechnology, 3rd ed. Wiley-VCH, ''in press''.]]
=== Rosetta stone (gene fusion) method ===
''The Rosetta Stone or Domain Fusion method'' is based on the hypothesis that interacting proteins are sometimes fused into a single protein.<ref name="Marcotte" /> For instance, two or more separate proteins in a genome may be identified as fused into one single protein in another genome. The separate proteins are likely to interact and thus are likely functionally related. An example of this is the ''Human Succinyl coA Transferase'' enzyme, which is found as one protein in humans but as two separate proteins, ''[[Acetate CoA-transferase|Acetate coA Transferase]] alpha'' and ''Acetate coA Transferase beta'', in ''Escherichia coli''.<ref name="Marcotte" /> In order to identify these sequences, a sequence similarity algorithm such as the one used by ''[[BLAST (biotechnology)|BLAST]]'' is necessary. For example, if we had the amino acid sequences of proteins A and B and the amino acid sequences of all proteins in a certain genome, we could check each protein in that genome for non-overlapping regions of sequence similarity to both proteins A and B. ''Figure B'' depicts the BLAST sequence alignment of Succinyl coA Transferase with its two separate homologs in E. coli. The two subunits have non-overlapping regions of sequence similarity with the human protein, indicated by the pink regions, with the alpha subunit similar to the first half of the protein and the beta similar to the second half. One limit of this method is that not all proteins that interact can be found fused in another genome, and therefore cannot be identified by this method. On the other hand, the fusion of two proteins does not necessitate that they physically interact. For instance, the [[SH2 ___domain|SH2]] and [[SH3 ___domain|SH3]] domains in the [[Src family kinase|src protein]] are known to interact. However, many proteins possess homologs of these domains and they do not all interact.<ref name="Marcotte" />
[[File:Trp operon organization across three different bacterial species.png|thumb|''FigureC.'' Organization of the trp operon in three different species of bacteria: ''Escherichia coli'', ''Haemophilus influenzae'', ''Helicobacter pylori''. Only the trpA and trpB genes are adjacent across all three organisms and are thus predicted to interact by the conserved gene neighborhood method. This image was adapted from Dandekar, T., Snel, B., Huynen, M., & Bork, P. (1998). Conservation of gene order: a fingerprint of proteins that physically interact. ''Trends in biochemical sciences'', ''23''(9), 324-328.<ref name="Dandekar" />]]
=== Conserved gene neighborhood ===
The conserved neighborhood method is based on the hypothesis that if genes encoding two proteins are neighbors on a chromosome in many genomes, then they are likely functionally related. The method is based on an observation by Bork et al. of gene pair conservation across nine bacterial and archaeal genomes. The method is most effective in prokaryotes with operons as the organization of genes in an operon is generally related to function.<ref name=":1">{{Cite journal|date=1998-09-01|title=Conservation of gene order: a fingerprint of proteins that physically interact|journal=Trends in Biochemical Sciences|language=en|volume=23|issue=9|pages=324–328|doi=10.1016/S0968-0004(98)01274-2|issn=0968-0004|last1=Dandekar|first1=T.|pmid=9787636}}</ref> For instance, the ''trpA'' and ''trpB'' genes in ''[[Escherichia coli]]'' encode the two subunits of the ''[[tryptophan synthase]]'' enzyme known to interact to catalyze a single reaction. The adjacency of these two genes was shown to be conserved across nine different bacterial and archaeal genomes.<ref name=":1" />
=== Classification methods ===
Classification methods use data to train a program (classifier) to distinguish positive examples of interacting protein/___domain pairs with negative examples of non-interacting pairs. Popular classifiers used are Random Forest Decision (RFD) and Support Vector Machines. RFD produces results based on the ___domain composition of interacting and non-interacting protein pairs. When given a protein pair to classify, RFD first creates a representation of the protein pair in a vector.{{r|Chen}} The vector contains all the ___domain types used to train RFD, and for each ___domain type the vector also contains a value of 0, 1, or 2. If the protein pair does not contain a certain ___domain, then the value for that ___domain is 0. If one of the proteins of the pair contains the ___domain, then the value is 1. If both proteins contain the ___domain, then the value is 2.{{r|Chen}} Using training data, RFD constructs a decision forest, consisting of many decision trees. Each decision tree evaluates several domains, and based on the presence or absence of interactions in these domains, makes a decision as to if the protein pair interacts. The vector representation of the protein pair is evaluated by each tree to determine if they are an interacting pair or a non-interacting pair. The forest tallies up all the input from the trees to come up with a final decision.{{r|Chen}} The strength of this method is that it does not assume that domains interact independent of each other. This makes it so that multiple domains in proteins can be used in the prediction.{{r|Chen}} This is a big step up from previous methods which could only predict based on a single ___domain pair. The limitation of this method is that it relies on the training dataset to produce results. Thus, usage of different training datasets could influence the results. A caveat of most methods is the lacks negative data, e.g non-interactions for proteins which can be overcome using topology-driven negative sampling.<ref>{{Citation |last=Chatterjee |first=Ayan |title=Topology-Driven Negative Sampling Enhances Generalizability in Protein-Protein Interaction Prediction |date=2024-04-29 |url=https://www.biorxiv.org/content/10.1101/2024.04.27.591478v1 |access-date=2024-05-04 |language=en |doi=10.1101/2024.04.27.591478 |last2=Ravandi |first2=Babak |last3=Philip |first3=Naomi H. |last4=Abdelmessih |first4=Mario |last5=Mowrey |first5=William R. |last6=Ricchiuto |first6=Piero |last7=Liang |first7=Yupu |last8=Ding |first8=Wei |last9=Mobarec |first9=Juan C.|doi-access=free }}</ref>
=== Inference of interactions from homologous structures===
This group of methods{{r|Aloy}}{{r|Chen}}{{r|Fukuhara}}{{r|Kittichotirat}}{{r|Ibis}} makes use of known protein complex structures to predict and structurally model interactions between query protein sequences. The prediction process generally starts by employing a sequence based method (e.g. [[Interolog]]) to search for protein complex structures that are homologous to the query sequences. These known complex structures are then used as templates to structurally model the interaction between query sequences. This method has the advantage of not only inferring protein interactions but also suggests models of how proteins interact structurally, which can provide some insights into the atomic level mechanism of that interaction. On the other hand, the ability for these methods to make a prediction is constrained by a limited number of known protein complex structures.
===
Association methods look for characteristic sequences or motifs that can help distinguish between interacting and non-interacting pairs. A classifier is trained by looking for sequence-signature pairs where one protein contains one sequence-signature, and its interacting partner contains another sequence-signature.{{r|Sprinzak}} They look specifically for sequence-signatures that are found together more often than by chance. This uses a log-odds score which is computed as log2(Pij/PiPj), where Pij is the observed frequency of domains i and j occurring in one protein pair; Pi and Pj are the background frequencies of domains i and j in the data. Predicted ___domain interactions are those with positive log-odds scores and also having several occurrences within the database.{{r|Sprinzak}} The downside with this method is that it looks at each pair of interacting domains separately, and it assumes that they interact independently of each other.
=== Identification of structural patterns ===
This method{{r|Aytuna}}{{r|Ogmen}} builds a library of known protein–protein interfaces from the [[Protein Data Bank|PDB]], where the interfaces are defined as pairs of polypeptide fragments that are below a threshold slightly larger than the [[Van der Waals radius]] of the atoms involved. The sequences in the library are then clustered based on structural alignment and redundant sequences are eliminated. The residues that have a high (generally >50%) level of frequency for a given position are considered hotspots.{{r|Keskin}} This library is then used to identify potential interactions between pairs of targets, providing that they have a known structure (i.e. present in the [[Protein Data Bank|PDB]]).
=== Bayesian network modelling ===
[[Bayesian method]]s{{r|Jansen}} integrate data from a wide variety of sources, including both experimental results and prior computational predictions, and use these features to assess the likelihood that a particular potential protein interaction is a true positive result. These methods are useful because experimental procedures, particularly the yeast two-hybrid experiments, are extremely noisy and produce many false positives, while the previously mentioned computational methods can only provide circumstantial evidence that a particular pair of proteins might interact.<ref>{{cite journal|last1=Zhang|first1=QC|last2=Petrey|first2=D|last3=Deng|first3=L|last4=Qiang|first4=L|last5=Shi|first5=Y|last6=Thu|first6=CA|last7=Bisikirska|first7=B|last8=Lefebvre|first8=C|last9=Accili|first9=D|last10=Hunter|first10=T|last11=Maniatis|first11=T|last12=Califano|first12=A|last13=Honig|first13=B|year=2012|title=Structure-based prediction of protein–protein interactions on a genome-wide scale|journal=Nature|volume=490|issue=7421|pages=556–60|doi=10.1038/nature11503|pmid=23023127|pmc=3482288 |bibcode=2012Natur.490..556Z}}</ref>
=== Domain-pair exclusion analysis ===
The ___domain-pair exclusion analysis{{r|Shoemaker}} detects specific ___domain interactions that are hard to detect using Bayesian methods. Bayesian methods are good at detecting nonspecific promiscuous interactions and not very good at detecting rare specific interactions. The ___domain-pair exclusion analysis method calculates an E-score which measures if two domains interact. It is calculated as log(probability that the two proteins interact given that the domains interact/probability that the two proteins interact given that the domains don’t interact). The probabilities required in the formula are calculated using an Expectation Maximization procedure, which is a method for estimating parameters in statistical models. High E-scores indicate that the two domains are likely to interact, while low scores indicate that other domains form the protein pair are more likely to be responsible for the interaction. The drawback with this method is that it does not take into account false positives and false negatives in the experimental data.
=== Supervised learning problem ===
The problem of PPI prediction can be framed as a supervised learning problem. In this paradigm the known protein interactions supervise the estimation of a function that can predict whether an interaction exists or not between two proteins given data about the proteins (e.g., expression levels of each gene in different experimental conditions, ___location information, phylogenetic profile, etc.).
== Relationship to docking methods ==
The field of protein–protein interaction prediction is closely related to the field of [[protein-protein docking|protein–protein docking]], which attempts to use geometric and steric considerations to fit two proteins of known structure into a bound complex. This is a useful mode of inquiry in cases where both proteins in the pair have known structures and are known (or at least strongly suspected) to interact, but since so many proteins do not have experimentally determined structures, sequence-based interaction prediction methods are especially useful in conjunction with experimental studies of an organism's [[interactome]].
== See also ==
*[[Interactome]]
*[[Protein–protein interaction]]
*[[Protein function prediction]]
*[[Protein structure prediction]]
*[[Protein structure prediction software]]
*[[Gene prediction]]
*[[Macromolecular docking]]
*[[Protein–DNA interaction site predictor]]
*[[Two-hybrid screening]]
*[[FastContact]]
== References ==
{{reflist|refs=
<ref name="Dandekar">Dandekar T., Snel B.,Huynen M. and Bork P. (1998) "Conservation of gene order: a fingerprint of proteins that physically interact." ''Trends Biochem. Sci.'' (23),324-328</ref>
<ref name="Enright">Enright A.J.,Iliopoulos I.,Kyripides N.C. and Ouzounis C.A. (1999) "Protein interaction maps for complete genomes based on gene fusion events." ''Nature'' (402), 86-90</ref>
<ref name="Marcotte">Marcotte E.M., Pellegrini M., Ng H.L., Rice D.W., Yeates T.O., Eisenberg D. (1999) "Detecting protein function and protein–protein interactions from genome sequences." ''Science'' (285), 751-753</ref>
<ref name="Pazos">{{cite journal|last1=Pazos|first1=F.|last2=Valencia|first2=A.|year=2001|title=Similarity of phylogenetic trees as indicator of protein–protein interaction|journal=Protein Engineering|volume=9|issue=14|pages=609–614 |doi=10.1093/protein/14.9.609|pmid=11707606|doi-access=free}}</ref>
<!-- not being used in article
<ref name="Pellegrini">{{cite journal|last1=Pellegrini|first1=M|last2=Marcotte|first2=EM|last3=Thompson|first3=MJ|last4=Eisenberg|first4=D|last5=Yeates|first5=TO|year=1999|title=Assigning protein functions by comparative genome analysis: protein phylogenetic profiles|journal=Proc Natl Acad Sci U S A|volume=96|pages=4285–8|doi=10.1073/pnas.96.8.4285|pmid=10200254|pmc=16324|bibcode=1999PNAS...96.4285P}}</ref>-->
<ref name="Tan">Tan S.H., Zhang Z., Ng S.K. (2004) "ADVICE: Automated Detection and Validation of Interaction by Co-Evolution." ''Nucleic Acids Res.'', ''32'' (Web Server issue):W69-72.</ref>
<ref name="Aytuna">{{cite journal|last1=Aytuna|first1=A. S.|last2=Keskin|first2=O.|last3=Gursoy|first3=A.|year=2005|title=Prediction of protein–protein interactions by combining structure and sequence conservation in protein interfaces|journal=Bioinformatics|volume=21|issue=12|pages=2850–2855|doi=10.1093/bioinformatics/bti443|pmid=15855251|doi-access=free}}</ref>
<ref name="Ogmen">{{cite journal|last1=Ogmen|first1=U.|last2=Keskin|first2=O.|last3=Aytuna|first3=A.S.|last4=Nussinov|first4=R.|last5=Gursoy|first5=A.|year=2005|title=PRISM: protein interactions by structural matching|journal=Nucleic Acids Res.|volume=33|issue=Web Server issue|pages=W331–336|doi=10.1093/nar/gki585|pmid=15991339|pmc=1160261|doi-access=free}}</ref>
<ref name="Keskin">{{cite journal|last1=Keskin|first1=O.|last2=Ma|first2=B.|last3=Nussinov|first3=R.|year=2004|title=Hot regions in protein–protein interactions: The organization and contribution of structurally conserved hot spot residues|journal=J. Mol. Biol.|volume=345|issue=5|pages=1281–1294|doi=10.1016/j.jmb.2004.10.077|pmid=15644221}}</ref>
<ref name="Jansen">{{cite journal|last1=Jansen|first1=R|last2=Yu|first2=H|last3=Greenbaum|first3=D|last4=Kluger|first4=Y|last5=Krogan|first5=NJ|last6=Chung|first6=S|last7=Emili|first7=A|last8=Snyder|first8=M|last9=Greenblatt|first9=JF|last10=Gerstein|first10=M|year=2003|title=A Bayesian networks approach for predicting protein–protein interactions from genomic data|journal=Science|volume=302|issue=5644|pages=449–53|doi=10.1126/science.1087361|pmid=14564010|bibcode=2003Sci...302..449J|citeseerx=10.1.1.217.8151|s2cid=5293611}}</ref>
<ref name="Aloy">{{cite journal|last1=Aloy|first1=P.|last2=Russell|first2=R. B.|year=2003|title=InterPreTS: protein Interaction Prediction through Tertiary Structure|journal=Bioinformatics|volume=19|issue=1|pages=161–162|doi=10.1093/bioinformatics/19.1.161|pmid=12499311|doi-access=free}}</ref>
<ref name="Chen">{{cite journal|last1=Chen|first1=XW|last2=Liu|first2=M|year=2005|title=Prediction of protein–protein interactions using random decision forest framework|journal=Bioinformatics|volume=21|issue=24|pages=4394–4400|doi=10.1093/bioinformatics/bti721|pmid=16234318|doi-access=free}}</ref>
<ref name="Fukuhara">Fukuhara, Naoshi, and Takeshi Kawabata. (2008) "HOMCOS: a server to predict interacting protein pairs and interacting sites by homology modeling of complex structures" ''Nucleic Acids Research'', ''36'' (S2): 185-.</ref>
<ref name="Kittichotirat">Kittichotirat W, M Guerquin, RE Bumgarner, and R Samudrala (2009) "Protinfo PPC: a web server for atomic level prediction of protein complexes" ''Nucleic Acids Research'', ''37'' (Web Server issue): 519-25.</ref>
<ref name="PazosRanea">{{cite journal|last1=Pazos|first1=F|last2=Ranea|first2=JA|last3=Juan|first3=D|last4=Sternberg|first4=MJ|year=2005|title=Assessing protein coevolution in the context of the tree of life assists in the prediction of the interactome|journal=J Mol Biol|volume=352|issue=4|pages=1002–1015|doi=10.1016/j.jmb.2005.07.005|pmid=16139301}}</ref>
<ref name="Sprinzak">{{cite journal|last1=Sprinzak|first1=E|last2=Margalit|first2=H|year=2001|title=Correlated sequence-signatures as markers of protein–protein interaction|journal=J Mol Biol|volume=311|issue=4|pages=681–692|doi=10.1006/jmbi.2001.4920|pmid=11518523}}</ref>
<ref name="Shoemaker">{{cite journal|last1=Shoemaker|first1=BA|last2=Panchenko|first2=AR|year=2007|title=Deciphering protein–protein interactions. Part II. Computational methods to predict protein and ___domain interaction partners|journal=PLOS Comput Biol|volume=3|issue=4|page=e43|doi=10.1371/journal.pcbi.0030043|pmid=17465672|pmc=1857810|bibcode=2007PLSCB...3...43S |doi-access=free }}</ref>
<ref name="Ibis">{{cite journal|last1=Shoemaker|first1=BA|last2=Zhang|first2=D|last3=Thangudu|first3=RR|last4=Tyagi|first4=M|last5=Fong|first5=JH|last6=Marchler-Bauer|first6=A|last7=Bryant|first7=SH|last8=Madej|first8=T|last9=Panchenko|first9=AR|date=Jan 2010|title=Inferred Biomolecular Interaction Server--a web server to analyze and predict protein interacting partners and binding sites|journal=Nucleic Acids Res|volume=38|issue=Database issue|pages=D518–24|pmid=19843613|doi=10.1093/nar/gkp842|pmc=2808861}}</ref>
<!-- not being used
<ref name="Marsh">{{cite journal|last1=Marsh|first1=J|last2=Hernandez|first2=H|last3=Hall|first3=Z|last4=Ahnert|first4=S|last5=Perica|first5=T|last6=Robinson|first6=C|last7=Teichmann|first7=S|year=2013|title=Protein complexes are under evolutionary selection to assemble via ordered pathways|journal=Cell|volume=153|issue=2|pages=461–70|doi=10.1016/j.cell.2013.02.044|pmid=23582331|pmc=4009401}}</ref>-->
}}
==External links==
*[http://openwetware.org/wiki/Protein–protein_interaction_databases Overview of protein interaction databases]
{{Prone to spam|date=December 2015}}
<!-- {{No more links}}
Please be cautious adding more external links.
Wikipedia is not a collection of links and should not be used for advertising.
Excessive or inappropriate links will be removed.
See [[Wikipedia:External links]] and [[Wikipedia:Spam]] for details.
If there are already suitable links, propose additions or replacements on
the article's talk page, or submit your link to the relevant category at
DMOZ (dmoz.org) and link there using {{Dmoz}}.
-->
{{Protein methods}}
{{DEFAULTSORT:Protein-protein interaction prediction}}
[[Category:Proteomics]]
|