Computational prediction of protein interactions related to the invasion of erythrocytes by malarial parasites
© Liu et al.; licensee BioMed Central Ltd. 2014
Received: 18 July 2014
Accepted: 19 November 2014
Published: 30 November 2014
The invasion of red blood cells (RBCs) by malarial parasites is an essential step in the life cycle of Plasmodium falciparum. Human-parasite surface protein interactions play a critical role in this process. Although several interactions between human and parasite proteins have been discovered, the mechanism related to invasion remains poorly understood because numerous human-parasite protein interactions have not yet been identified. High-throughput screening experiments are not feasible for malarial parasites due to difficulty in expressing the parasite proteins. Here, we performed computational prediction of the PPIs involved in malaria parasite invasion to elucidate the mechanism by which invasion occurs.
In this study, an expectation maximization algorithm was used to estimate the probabilities of domain-domain interactions (DDIs). Estimates of DDI probabilities were then used to infer PPI probabilities. We found that our prediction performance was better than that based on the information of D. melanogaster alone when information related to the six species was used. Prediction performance was assessed using protein interaction data from S. cerevisiae, indicating that the predicted results were reliable. We then used the estimates of DDI probabilities to infer interactions between 490 parasite and 3,787 human membrane proteins. A small-scale dataset was used to illustrate the usability of our method in predicting interactions between human and parasite proteins. The positive predictive value (PPV) was lower than that observed in S. cerevisiae. We integrated gene expression data to improve prediction accuracy and to reduce false positives. We identified 80 membrane proteins highly expressed in the schizont stage by fast Fourier transform method. Approximately 221 erythrocyte membrane proteins were identified using published mass spectral datasets. A network consisting of 205 interactions was predicted. Results of network analysis suggest that SNARE proteins of parasites and APP of humans may function in the invasion of RBCs by parasites.
We predicted a small-scale PPI network that may be involved in parasite invasion of RBCs by integrating DDI information and expression profiles. Experimental studies should be conducted to validate the predicted interactions. The predicted PPIs help elucidate the mechanism of parasite invasion and provide directions for future experimental investigations.
KeywordsPlasmodium falciparum Expectation maximization Fast Fourier transform Invasion Red blood cell
The ability of a malaria parasite to invade red blood cells (RBCs) is a key step for its survival and pathogenesis. The malarial parasites can successfully invade RBCs by several mechanisms, including rhoptry and microneme protein secretion, tight junction formation, and acto-myosin motor assembly ,. Previous investigations revealed that surface protein interactions between parasite and human surface proteins are necessary to allow parasites to invade RBCs. However, although several surface proteins were identified to be involved in this process -, the molecular basis of invasion remains poorly understood. Indeed, specific interactions between parasite proteins and human proteins remain largely undiscovered. Elucidation of the molecular basis of such interactions is critical to develop new intervention strategies against malaria.
Although high-throughput biological experiments such as yeast two-hybridization and avidity-based extracellular interaction screen can be conducted to reveal interactions between parasite proteins and human proteins ,, there are several technical difficulties in conducting such experiments with malarial parasites. A major difficulty is in expression of parasite proteins in yeasts or mammalian cells because of high AT content (approximately 80%) in the genome of P. falciparum . For this reason, large-scale experimental validation of protein interaction networks is often ruled out in the case of malaria parasite. Moreover, inherent shortcomings of high-throughput techniques, such as bias toward yeast homologs, large false positive rates, and low coverage, could contribute to the significant inaccuracies of projected interactions . Computational methods have been developed and used to predict protein interactions when experimental methods are not feasible. For example, Date et al. utilized a Bayesian framework to integrate different data sources, including polygenetic profiles, gene expression profiles, and Rosetta stone fusion genes . Positive and negative high confidence (HC) datasets derived from Gene Ontology and KEGG databases were used to determine the likelihood scores of each interaction. By setting a threshold, they generated a protein–protein interaction (PPI) network with sensitivity of about 21%. Based on the resulting network, several uncharacterized proteins were assigned to various biological processes. Another approach based on protein sequence similarity was developed and implemented to predict putative protein interactions between human and malaria parasite . Candidate interactions were then assessed by random forest classification and further filtered in terms of expression and molecular characteristics. The resulting network revealed that parasites possibly utilize their proteins in a combined manner by predominantly targeting hub proteins. Although several predicted protein networks have been constructed, predictions of membrane protein interactions related to parasite invasion have not been conducted before.
In this study, membrane protein interactions between human and P. falciparum were predicted to elucidate the protein interactions involved in parasite invasion of RBCs. Considering that a protein domain serves as a unit of protein-protein interactions and is evolutionally conserved, a model was developed to relate protein interaction probabilities with domain interaction probabilities. In the present study, an expectation maximization (EM) algorithm proposed by Liu et al. was used to estimate the probabilities of domain-domain interactions (DDIs) . The EM algorithm employs a likelihood-based approach and exhibits good performance in estimating DDI probabilities . In this approach, PPIs and DDIs were treated as random variables. The probabilities of DDIs were computed on the basis of information of PPIs after false positive rate (f p ) and false negative rate (f n ) were specified. DDI estimates were then used to infer plausible interactions between human and parasite membrane proteins.
Collection of physical protein-protein interaction data
The protein-domain relationships of each protein in six species, including Arabidopsis thaliana, Caenorhabditis elegans, Drosophila melanogaster, Homo sapiens, Mus musculus, and Schizosaccharomyces pombe, were extracted from the PFAM database . To reduce the false positive rate (FPR) of protein-domain relations, we defined a significant protein-domain relationship (PDR) as one with E-value less than 1E-5 for proteins in A. thaliana, C. elegans, D. melanogaster, H. sapiens, M. musculus, and less than 1E-4 for proteins in S. pombe. Proteins with significant PDRs were represented by their respective Entrez gene ID or gene symbol. The physical protein-protein interaction data of these six organisms were collected from the BioGRID database . We only considered interactions among proteins with significant PDRs. After removing redundant interactions, we obtained 9,478 protein interactions from A. thaliana, 3,052 interactions from C. elegans, 21,752 interactions from D. melanogaster, 94,396 interactions from H. sapiens, 7,409 interactions from M. musculus, and 3,828 interactions from S. pombe. Detailed information on the collected protein interactions for each species used in our study can be found in Additional file 1.
High-confidence dataset preparation
Three high confidence datasets were constructed and used to evaluate the reliability of prediction. The first HC dataset contained DDIs collected from iPFAM and 3DID databases ,. A total of 6,634 structurally verified physical domain interactions were found in these two databases. After removing domains not found in D. melanogaster, we constructed a DDI matrix containing information of interactions among 1,515 domains. The second high confidence dataset consisted of protein interactions from S. cerevisiae, in which physical protein interactions have been widely investigated. The PPIs in S. cerevisiae were also collected from BioGrid database. We defined PDR as significant when its E-value was less than 1E-4. After removing proteins containing domains not found in the six species studied here, we obtained 49,084 interactions among 3,960 proteins. We considered these 49,084 interactions as positive interactions and the remaining protein pairs were considered to be non-interacting (negative interactions). The last high confidence dataset was a small-scale dataset that contained 456 experimentally determined interactions between human and parasite proteins . After removing proteins that do not satisfy the PDR condition, we obtained 132 interactions between 66 parasite proteins and 107 human proteins. The first and second datasets were used to evaluate the reliabilities of our prediction in DDIs and PPIs, respectively. The last dataset was utilized to assess the prediction performance of protein interactions between humans and parasites.
where N is the number of protein pairs containing a domain pair (m, n). The initial value of λ mn was expressed as . We can estimate λ mn by iterating between E- and M-steps to obtain the maximum likelihood estimation of λ mn . If was 0, then the non-zero λ mn was updated in the EM algorithm; however, computational consumption was very high.
the value of λ mn was set to zero and not used in the E-step. In the present study, p = 0.00143 was the observed protein interaction probability and computed using the number of the experimentally confirmed PPIs to divide the total number of protein pairs across the six species. Therefore, we only considered λ mn with initial values that were significantly higher than the observed protein interaction probability (p).
The EM algorithm was implemented in Matlab 2012b for two months by using the parallel computing toolbox. EM steps were repeated until the difference between two consecutive steps was <0.01 or the number of repeats >40. Estimates of λ mn were then obtained and used for protein interaction prediction.
FFT analysis of microarray dataset
With Fourier analysis, periodic genes within a genome can be identified, and the peak expression time of each gene can be determined.
We performed FFT analysis on time course microarrays in the asexual intraerythrocytic developmental cycle (IDC) of P. falciparum . These expression profiles were obtained at intervals of 1 h for 48 h except on the 23rd and 29th h. The missing data were filled using the k nearest neighbor (KNN) method . Before FFT analysis was performed, the expression value of each gene was centered by subtracting the mean value.
Identification of membrane proteins
Proteins located in the membrane of human erythrocyte or malaria parasite should contain a transmembrane domain (TM) in their structures. They may or may not contain cleavable signal peptides in their sequences. To identify the human and malaria parasite membrane proteins that may interact, we used SignaIP and TMHMM to predict membrane proteins ,. The sequence regions where predicted signal peptide and TM domain coexist were excluded because signal peptides are likely to be incorrectly predicted as TM domains by TMHMM. Proteins with multiple TM domains in their sequences were considered to be membrane proteins. The SignaIP program was run under default parameters except the D-cutoff value was adjusted to 0.33. The parameters for TMHMM were retained.
Receiver operating characteristic (ROC) curve
If the predicted result showed high PPV, a high success rate could be obtained in experimental validation of predicted protein interactions.
Gene ontology (GO) enrichment analysis
To understand the biological meaning of a large list of genes identified using the FFT method, we performed GO term enrichment analysis on these genes via the DAVID web server (http://david.abcc.ncifcrf.gov/) . The GO terms of the gene molecular function were selected, and a gene set was considered to be enriched if the p-value was <0.05.
Estimation of DDI probabilities using the EM algorithm
Protein domains are the structural units of proteins, and PPIs are mostly achieved through DDIs. Therefore, we first needed to estimate the DDI probabilities before inferring PPI probabilities. In our analysis, an EM algorithm was introduced to estimate the probability of each DDI.
f n and f p used in the EM algorithm. t represents assumed average number of interaction partners
Considering that the computational consumption of estimating all DDI probabilities was very high, we assigned the initial low λ mn as zero (see Methods). The initial λ mn significantly higher than the observed PPI probability (p = 0.00143) was used in EM algorithm. This simplification indicated that only parts of non-zero λ mn were used for computation. To assess the effect of this alteration on the prediction accuracy of the EM algorithm, we fixed f n and f p , and the ROC curves obtained using totally non-zero λ mn were then compared with those obtained using partly non-zero λ mn . The ROC curves were nearly identical when the cut off value was set at ≥0.05 (Figure 1b). Slight differences were observed when the cut off value was <0.05. Similar trends were observed when f n and f p were assigned with different values (Additional file 2). Despite the slight differences in the ROC curves obtained using totally or partly non-zero λ mn , the time consumed by running the EM algorithm by using partly non-zero λ mn was reduced by almost half.
Previous investigation had established that the prediction performance based on multiple organisms outperformed the performance obtained from only a single organism. Therefore, we integrated the protein interaction data from six organisms to estimate DDI probabilities. Our computer memory was limited; as such, the program could not be run using totally non-zero λ mn . Therefore, we treated λ mn with a low value as zero and used the rest of the non-zero λ mn to estimate DDI probabilities. Indeed, the prediction accuracy based on six organisms was higher than that based on the protein interaction data of a single organism. The maximum TPR level increased from 0.209 to 0.616 (Figure 1c). With the the same FPR value, the TPR value observed from the information of six organisms was apparently higher than that observed from D. melanogaster alone. On the basis of this result, we used the estimates of DDI probabilities to compute PPI probabilities.
Before estimating the interaction probabilities between human and parasite membrane proteins by using estimates of DDI probabilities, we accessed the performance of those DDI probabilities in inferring the PPIs of S. cerevisiae. The probability of each PPI was computed using the formula of P inter (see Methods). At a threshold of 0.05, 117,805 protein pairs were predicted to be interacting; among these interactions, 7,117 protein interactions were identified (Figure 1d, upper panel). The TPR of this method was 0.145, and PPV reached 0.0604, which was nearly 42 times higher than the observed interaction probability (p = 0.00143). These results suggested that the PPIs inferred using the estimates of DDI probabilities were reliable. At a threshold of 0.01, 299,602 proteins were predicted to be interacting (Figure 1d, lower panel) and TPR and PPV were 0.235 and 0.0385, respectively. Although a decrease in the threshold value increased TPR, FPR was also found to be increased.
Prediction of interactions between parasite and human membrane proteins
Considering that significantly limited interactions were observed between human and parasite membrane proteins, we could not easily evaluate the prediction accuracy of our method. To demonstrate the viability of our method in predicting interactions between proteins from different species, 456 experimentally confirmed interactions between human and parasite proteins were used. After removing proteins not fulfilling the PDR conditions, we obtained a network consisting of 132 interactions between 66 parasite membrane proteins and 107 human membrane proteins. The observed interaction probability of this small-scale network was p = 0.0187. We then used this small-scale network to roughly assess the performance of our method in predicting interactions between human and parasite proteins. The interaction probabilities of these proteins were calculated using estimates of DDI probabilities. Figure 2b shows that 1,282 protein interaction probabilities were >0. The maximum TPR was 0.2270. At a threshold of 0.05, 11 of 317 protein interactions were predicted. TPR and PPV were 0.0833 and 0.0347, respectively, which were evidently lower than those obtained in S. cerevisiae. Even when the threshold was set at 0.01, only 15 interactions were predicted to be interacting. When TPR was increased to 0.114, PPV decreased to 0.0224, which was slightly higher than the observed interaction probability (p = 0.0187). The low TPR may be caused by several protein interactions in this network being false positives, since the observed interaction probability (p = 0.0187) was apparently larger than the interaction probability observed in the aforementioned six species (p = 0.00143). Low PPV suggests that integration of other information is necessary to improve our prediction accuracy.
We excluded falsely predicted PPIs by integrating the expression profiles of merozoite and human erythrocytes to improve the PPV of PPI prediction and to provide the direction for future experimental studies. The interactions between merozoite and human erythrocyte membrane proteins are mostly responsible for parasite invasion. Therefore, highly expressed proteins in merozoites and human erythrocytes should be identified to improve prediction accuracy.
FFT identification of proteins highly expressed on merozoite membrane
To identify highly expressed genes in merozoites, we only need to extract the genes highly expressed in the schizont stage because mature schizonts consist of tens of merozoites. We performed FFT analysis on a published microarray dataset . Probes whose spectral amplitudes were >0.5 at ω = 1 were considered to be significant periodic probes. The peak expression time of each periodic gene was calculated using a previously described method (T peak calculation, see Methods). The profiles of the periodic genes were organized by increasing the time of peak expression (Figure 3b). The results of Giemsa staining of the major morphological stages throughout the IDC had indicated that schizonts initially appeared at 25 h after the erythrocytes were invaded. In our study, genes with peak expression times at intervals of 30 h and 48 h were considered highly expressed in merozoites. A total of 3,442 genes contained significant frequency components, of which 800 were highly expressed in the schizont stage.
GO term enrichment analysis of genes highly expressed at schizont stage
GO:0004674 protein serine/threonine kinase activity
GO:0004672 protein kinase activity
GO:0016773 phosphotransferase activity, alcohol group as acceptor
GO:0016301 kinase activity
GO:0005509 calcium ion binding
GO:0003779 actin binding
GO:0008092 cytoskeletal protein binding
GO:0016772 transferase activity, transferring phosphorus-containing groups
GO:0016491 oxidoreductase activity
GO:0004197 cysteine-type endopeptidase activity
GO:0008234 cysteine-type peptidase activity
GO:0004252 serine-type endopeptidase activity
GO:0042578 phosphoric ester hydrolase activity
GO:0005200 structural constituent of cytoskeleton
To infer interactions between merozoite and human erythrocyte proteins, we focused on membrane proteins of these 800 proteins expressed highly in the schizont stage. We overlapped these 800 proteins and 490 parasite membrane proteins previously identified and found that 80 membrane proteins were highly expressed in the schizont stage (Figure 3c). To identify membrane proteins of RBCs, we used published mass spectral datasets ,. A total of 2,396 proteins were found to be expressed in RBCs; among them, 221 were predicted to be membrane proteins and had significant PDRs.
Prediction of protein interactions involved in malaria parasite invasion of human erythrocytes
On the basis of previously computed 62,197 probabilities between human and parasite membrane proteins, we obtained a network consisting of 80 parasite proteins and 221 human proteins after the proteins not expressed in the schizont stage and those not found in RBCs were filtered out. A total of 960 interactions with probabilities >0 were found. At a threshold of 0.01, 467 interactions were predicted. Considering that experimentally discovered surface protein interactions are limited, we could not evaluate the accuracy of our prediction. Six non-interacting protein pairs were manually collected from the third high confidence dataset. These non-interacting protein pairs consisted of two parasite membrane proteins and three human membrane proteins. Our prediction results indicated that all of the probabilities of these protein pairs were 0, suggesting that no interactions were found between those proteins. Thus, our predicted results are consistent with the experimental results. Further experimental validation is needed to fully assess the accuracy of our predictions.
Proteins with node degrees higher than six in the predicted network
a. Human proteins with node degrees higher than six.
amyloid beta (A4) precursor protein
VAMP (vesicle-associated membrane protein)-associated protein A
VAMP (vesicle-associated membrane protein)-associated protein B and C
motile sperm domain containing 1
b. Parasite proteins with node degrees higher than six.
syntaxin, Qa-SNARE family
LEM3/CDC50 family protein
Mitochondrial calcium-proton antiporter
signal peptidase complex subunit 2
oligosaccharyl transferase STT3 subunit
Eight parasite proteins in this network exhibited degrees of interaction more than six and were selected for further analysis (Table 3b). Some of these proteins, including PFF0170w, PFC0912w, and PF11_0173, were false positives because they are not localized on the plasma membrane of parasites and hence cannot interact with human proteins. The parasite SNARE proteins PfSyn17 and PfSyn6 were predicted to interact with human proteins. The peak expression times of these two proteins based on the microarray data were at 42 h (Figure 4b). Considering that SNARE proteins mainly regulate vesicle exocytosis and are involved in protein secretion ,, we speculate that the inferred interaction between SNARE proteins and human erythrocyte proteins may stimulate protein secretion by microneme and rhoptry. This may in turn facilitate the ability of malaria parasite to invade erythrocytes. Thus, our results suggest that SNARE proteins may possibly be involved in parasite invasion.
On the basis of these network analysis results, we chose parasite SNARE proteins and human APPs for further biological experimental studies.
In the present study, we sought to elucidate the molecular basis of invasion of RBCs by malarial parasites. We used the EM algorithm to estimate the DDI probabilities based on currently available large-scale protein interaction datasets, and then used the estimates of DDI probabilities to predict interactions between human RBCs and merozoite membrane proteins. We also integrated gene expression information to reduce false positives and thus improve our prediction accuracy. Using FFT, we identified the genes highly expressed in the schizont stage. By combining estimated PPI probabilities and gene expression profiles, we predicted a network consisting of 205 PPIs between parasite and human membrane proteins. The results of further analysis suggested that SNARE proteins in parasites and APP in humans may play an important role in the invasion of erythrocytes.
Although we inferred possible interactions between human and malaria parasite proteins, our prediction did not capture several discovered protein interactions related to parasite invasion of RBCs because we only considered possible interactions between proteins that satisfy the PDR condition. For instance, previous experimental investigations reported that several interactions between parasite proteins and human proteins, including MSP1-Band3 and EBA175-Glycophorin A ,, are involved in invasion. Our predicted interactions did not include these two interactions because the parasite proteins MSP1and EBA175 do not contain any domain found in the six species and were therefore not used for PPI probability estimation. Actually, a total of 815 parasite proteins that may be located in the parasite membrane were identified. However, only 490 of those contain domains that are conserved in all six species. This finding indicated that our method could not be used to predict whether the 325 remaining proteins interact with human membrane proteins. Furthermore, a total of 839 proteins were identified from the merozoite proteome , of which 62 were assumed to be surface proteins. Only 22 surface proteins were found to share recognizable domains across the six organisms, which implied that interactions mediated by the remaining 40 proteins could not be inferred from our study. Expansion of the PPI collection could overcome this limitation. Another limitation of our method is that, since we only considered the information of domain-domain interactions, PPIs mediated by posttranslational modifications, such as those in GYPC-EBA140 and GYPB-EBL1 ,, cannot be inferred from our results. To our knowledge, there is no computational method suitable for inferring PPIs mediated by posttranslational modifications.
In addition, interactions mediated by secretory or peri-membrane proteins are also not predicted because they were excluded from our analysis for the reason that most secretory parasite proteins contain domains not found in all six organisms and identification of human peri-membrane proteins is not feasible. In our study, we only predicted interactions between membrane proteins in parasites and human RBCs. However, several secretory parasite proteins and human perimembrane proteins are known to be utilized by parasites in the invasion of RBCs. For instance, PfRh5 proteins secreted by parasite rhoptry are involved in parasite invasion of erythrocytes by directly interacting with the surface protein Basigin in RBCs . The human perimembrane protein SEMA7A interacts with parasite membrane proteins PfMtrap and participates in parasite invasion . Therefore, extending the protein candidates to human perimembrane and parasite secretory proteins is the best solution to capture interactions in which these kinds of proteins are involved.
We found that the results of incorporating information from six species outperformed those obtained from D. melanogaster alone, which allowed us to make more informative inferences regarding protein interactions. A network composed of interactions between human and parasite proteins was then inferred. Assessment of our predictions using a small-scale protein interaction dataset indicated that our prediction accuracy was lower than that obtained in S. cerevisiae. This phenomenon may be caused by several protein interactions being false positives in the small-scale dataset, and suggested that other information, such as gene expression profile and gene ontology information, should be integrated to improve prediction accuracy . In our analysis, gene expression profile data was incorporated to increase the prediction accuracy. In addition, GO information, particularly information related to protein localization, could be integrated to further improve the accuracy of our predictions. For example, in our predicted PPI network, several proteins were found to be located on the membrane of mitochondria and the endoplasmic reticulum. This result suggested that proteins located on the membrane of these two cellular components may be falsely predicted to interact with human membrane proteins. The exclusion of interactions in which such proteins participate possibly increased our prediction accuracy. Another key factor affecting prediction performance is the selection of threshold value. Predicted interactions depend on a chosen threshold value. Although an increase in the threshold value likely improves prediction accuracy, the TPR is evidently reduced. Considering the incompleteness of domain information and the unknown number of definitely interacting protein pairs, we could not easily set a threshold value to determine whether or not a protein pair exhibits interaction.
Despite these limitations, our predicted interactions revealed that SNARE proteins that participate in protein secretion may be involved in malaria parasite invasion of erythrocytes. This result is consistent with previous reports that protein secretion via microneme and rhoptry are necessary to induce parasite invasion of RBCs . Predicted interactions also indicated that human APP may function in parasite invasion. The erythrocyte proteome data shows that APP is expressed at a medium level in RBCs . Analysis of dataset from the 1000 genome database  indicated that several SNPs of APP, such as rs216772, rs2829992 and rs216773, are more prevalent in among Africans. Understanding the significance of the relation between these SNPs and parasite invasion requires further investigation. Considering that the functions of APP and SNARE in parasite invasion have not been elucidated previously, we will conduct a series of biological experiments to study them. Although experimental validation of these interactions was hampered by difficulty in expression of parasite proteins in mammalian cells, some encouraging progress has been made by malariaologists ,.
In the present study, the EM algorithm was used to estimate probabilities of DDIs, which were then applied to infer possible protein interactions between human RBCs and parasite merozoite membrane proteins. Gene expression data was integrated to filter out false positives. FFT was introduced to identify genes expressed highly in the schizont stage, while published data was utilized to search for membrane proteins of RBCs. A network consisting of 205 PPIs, including several novel interactions, was obtained. Our results help elucidate the molecular mechanism of parasite invasion and provide promising research directions for further experimental studies.
Xuewu Liu, Yuxiao Huang and Jiao Liang are co-first authors.
This work was supported by the National Science Foundation of China (Grant No. 31172167).
- Sharma P, Chitnis CE: Key molecular events during host cell invasion by Apicomplexan pathogens. Curr Opin Microbiol. 2013, 16 (4): 432-437. 10.1016/j.mib.2013.07.004.View ArticlePubMedGoogle Scholar
- Singh S, Alam MM, Pal-Bhowmick I, Brzostowski JA, Chitnis CE: Distinct external signals trigger sequential release of apical organelles during erythrocyte invasion by malaria parasites. PLoS Pathog. 2010, 6 (2): e1000746-10.1371/journal.ppat.1000746.View ArticlePubMed CentralPubMedGoogle Scholar
- Crosnier C, Bustamante LY, Bartholdson SJ, Bei AK, Theron M, Uchikawa M, Mboup S, Ndir O, Kwiatkowski DP, Duraisingh MT, Rayner JC, Wright GJ: Basigin is a receptor essential for erythrocyte invasion by Plasmodium falciparum. Nature. 2011, 480 (7378): 534-537.PubMed CentralPubMedGoogle Scholar
- Rydzak J, Kmiecik AM, Jaskiewicz E: [Human erythrocyte glycophorin C as the receptor for EBA-140 Plasmodium falciparum merozoite ligand]. Postepy Hig Med Dosw (Online). 2013, 67: 1331-1339. 10.5604/17322693.1081865.View ArticleGoogle Scholar
- Tham WH, Wilson DW, Lopaticki S, Schmidt CQ, Tetteh-Quarcoo PB, Barlow PN, Richard D, Corbin JE, Beeson JG, Cowman AF: Complement receptor 1 is the host erythrocyte receptor for Plasmodium falciparum PfRh4 invasion ligand. Proc Natl Acad Sci U S A. 2010, 107 (40): 17327-17332. 10.1073/pnas.1008151107.View ArticlePubMed CentralPubMedGoogle Scholar
- LaCount DJ, Vignali M, Chettier R, Phansalkar A, Bell R, Hesselberth JR, Schoenfeld LW, Ota I, Sahasrabudhe S, Kurschner C, Fields S, Hughes RE: A protein interaction network of the malaria parasite Plasmodium falciparum. Nature. 2005, 438 (7064): 103-107. 10.1038/nature04104.View ArticlePubMedGoogle Scholar
- Kerr JS, Wright GJ: Avidity-based extracellular interaction screening (AVEXIS) for the scalable detection of low-affinity extracellular receptor-ligand interactions. J Vis Exp. 2012, 61: e3881-PubMedGoogle Scholar
- Gardner MJ, Hall N, Fung E, White O, Berriman M, Hyman RW, Carlton JM, Pain A, Nelson KE, Bowman S, Paulsen IT, James K, Eisen JA, Rutherford K, Salzberg SL, Craig A, Kyes S, Chan MS, Nene V, Shallom SJ, Suh B, Peterson J, Angiuoli S, Pertea M, Allen J, Selengut J, Haft D, Mather MW, Vaidya AB, Martin DM, et al: Genome sequence of the human malaria parasite Plasmodium falciparum. Nature. 2002, 419 (6906): 498-511. 10.1038/nature01097.View ArticlePubMedGoogle Scholar
- Von Mering C, Krause R, Snel B, Cornell M, Oliver SG, Fields S, Bork P: Comparative assessment of large-scale data sets of protein-protein interactions. Nature. 2002, 417 (6887): 399-403. 10.1038/nature750.View ArticlePubMedGoogle Scholar
- Date SV, Stoeckert CJ: Computational modeling of the Plasmodium falciparum interactome reveals protein function on a genome-wide scale. Genome Res. 2006, 16 (4): 542-549. 10.1101/gr.4573206.View ArticlePubMed CentralPubMedGoogle Scholar
- Wuchty S: Computational prediction of host-parasite protein interactions between P. falciparum and H. sapiens. PLoS One. 2011, 6 (11): e26960-10.1371/journal.pone.0026960.View ArticlePubMed CentralPubMedGoogle Scholar
- Liu Y, Liu N, Zhao H: Inferring protein-protein interactions through high-throughput interaction data from diverse organisms. Bioinformatics. 2005, 21 (15): 3279-3285. 10.1093/bioinformatics/bti492.View ArticlePubMedGoogle Scholar
- Lee H, Deng M, Sun F, Chen T: An integrated approach to the prediction of domain-domain interactions. BMC Bioinformatics. 2006, 7: 269-10.1186/1471-2105-7-269.View ArticlePubMed CentralPubMedGoogle Scholar
- Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, Heger A, Hetherington K, Holm L, Mistry J, Sonnhammer EL, Tate J, Punta M: Pfam: the protein families database. Nucleic Acids Res. 2014, 42 (Database issue): D222-D230. 10.1093/nar/gkt1223.View ArticlePubMed CentralPubMedGoogle Scholar
- Chatr-Aryamontri A, Breitkreutz BJ, Heinicke S, Boucher L, Winter A, Stark C, Nixon J, Ramage L, Kolas N, O'Donnell L, Reguly T, Breitkreutz A, Sellam A, Chen D, Chang C, Rust J, Livstone M, Oughtred R, Dolinski K, Tyers M: The BioGRID interaction database: 2013 update. Nucleic Acids Res. 2013, 41 (Database issue): D816-D823. 10.1093/nar/gks1158.View ArticlePubMed CentralPubMedGoogle Scholar
- Mosca R, Ceol A, Stein A, Olivella R, Aloy P: 3did: a catalog of domain-based interactions of known three-dimensional structure. Nucleic Acids Res. 2014, 42 (Database issue): D374-D379. 10.1093/nar/gkt887.View ArticlePubMed CentralPubMedGoogle Scholar
- Finn RD, Miller BL, Clements J, Bateman A: iPfam: a database of protein family and domain interactions found in the Protein Data Bank. Nucleic Acids Res. 2014, 42 (Database issue): D364-D373. 10.1093/nar/gkt1210.View ArticlePubMed CentralPubMedGoogle Scholar
- Vignali M, McKinlay A, LaCount DJ, Chettier R, Bell R, Sahasrabudhe S, Hughes RE, Fields S: Interaction of an atypical Plasmodium falciparum ETRAMP with human apolipoproteins. Malar J. 2008, 7: 211-10.1186/1475-2875-7-211.View ArticlePubMed CentralPubMedGoogle Scholar
- Kim I, Liu Y, Zhao H: Bayesian methods for predicting interacting protein pairs using domain information. Biometrics. 2007, 63 (3): 824-833. 10.1111/j.1541-0420.2007.00755.x.View ArticlePubMedGoogle Scholar
- Kim J, Ogden RT, Kim H: A method to identify differential expression profiles of time-course gene data with Fourier transformation. BMC Bioinformatics. 2013, 14: 310-10.1186/1471-2105-14-310.View ArticlePubMed CentralPubMedGoogle Scholar
- Kim BR, Zhang L, Berg A, Fan J, Wu R: A computational approach to the functional clustering of periodic gene-expression profiles. Genetics. 2008, 180 (2): 821-834. 10.1534/genetics.108.093690.View ArticlePubMed CentralPubMedGoogle Scholar
- Silva GJ, Ushizima MR, Lessa PS, Cardoso L, Drager LF, Atala MM, Consolim-Colombo FM, Lopes HF, Cestari IA, Krieger JE, Krieger EM: Critical analysis of autoregressive and fast Fourier transform markers of cardiovascular variability in rats and humans. Braz J Med Biol Res. 2009, 42 (4): 386-396. 10.1590/S0100-879X2009000400012.View ArticlePubMedGoogle Scholar
- Bozdech Z, Llinas M, Pulliam BL, Wong ED, Zhu J, DeRisi JL: The transcriptome of the intraerythrocytic developmental cycle of Plasmodium falciparum. PLoS Biol. 2003, 1 (1): E5-10.1371/journal.pbio.0000005.View ArticlePubMed CentralPubMedGoogle Scholar
- Sakka E, Prentza A, Koutsouris D: Classification algorithms for microcalcifications in mammograms (Review). Oncol Rep. 2006, 15 Spec no: 1049-1055.PubMedGoogle Scholar
- Petersen TN, Brunak S, Von Heijne G, Nielsen H: SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011, 8 (10): 785-786. 10.1038/nmeth.1701.View ArticlePubMedGoogle Scholar
- Chen Y, Yu P, Luo J, Jiang Y: Secreted protein prediction system combining CJ-SPHMM, TMHMM, and PSORT. Mamm Genome. 2003, 14 (12): 859-865. 10.1007/s00335-003-2296-6.View ArticlePubMedGoogle Scholar
- Linden A: Measuring diagnostic and predictive accuracy in disease management: an introduction to receiver operating characteristic (ROC) analysis. J Eval Clin Pract. 2006, 12 (2): 132-139. 10.1111/j.1365-2753.2005.00598.x.View ArticlePubMedGoogle Scholar
- Da Huang W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4 (1): 44-57. 10.1038/nprot.2008.211.View ArticlePubMedGoogle Scholar
- Srinivasan P, Yasgar A, Luci DK, Beatty WL, Hu X, Andersen J, Narum DL, Moch JK, Sun H, Haynes JD, Maloney DJ, Jadhav A, Simeonov A, Miller LH: Disrupting malaria parasite AMA1-RON2 interaction with a small molecule prevents erythrocyte invasion. Nat Commun. 2013, 4: 2261-10.1038/ncomms3261.View ArticlePubMed CentralPubMedGoogle Scholar
- Greenbaum DC, Baruch A, Grainger M, Bozdech Z, Medzihradszky KF, Engel J, DeRisi J, Holder AA, Bogyo M: A role for the protease falcipain 1 in host cell invasion by the human malaria parasite. Science. 2002, 298 (5600): 2002-2006. 10.1126/science.1077426.View ArticlePubMedGoogle Scholar
- Bachi A, Simo C, Restuccia U, Guerrier L, Fortis F, Boschetti E, Masseroli M, Righetti PG: Performance of combinatorial peptide libraries in capturing the low-abundance proteome of red blood cells. 2. Behavior of resins containing individual amino acids. Anal Chem. 2008, 80 (10): 3557-3565. 10.1021/ac8001353.View ArticlePubMedGoogle Scholar
- Pasini EM, Kirkegaard M, Mortensen P, Lutz HU, Thomas AW, Mann M: In-depth analysis of the membrane and cytosolic proteome of red blood cells. Blood. 2006, 108 (3): 791-801. 10.1182/blood-2005-11-007799.View ArticlePubMedGoogle Scholar
- Rivera-Escalera F, Matousek SB, Ghosh S, Olschowka JA, O'Banion MK: Interleukin-1beta mediated amyloid plaque clearance is independent of CCR2 signaling in the APP/PS1 mouse model of Alzheimer’s disease. Neurobiol Dis. 2014, 69: 124-133. 10.1016/j.nbd.2014.05.018.View ArticlePubMed CentralPubMedGoogle Scholar
- Muller UC, Zheng H: Physiological functions of APP family proteins. Cold Spring Harb Perspect Med. 2012, 2 (2): a006288-10.1101/cshperspect.a006288.View ArticlePubMed CentralPubMedGoogle Scholar
- Smoot ME, Ono K, Ruscheinski J, Wang PL, Ideker T: Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. 2011, 27 (3): 431-432. 10.1093/bioinformatics/btq675.View ArticlePubMed CentralPubMedGoogle Scholar
- Chaineau M, Danglot L, Galli T: Multiple roles of the vesicular-SNARE TI-VAMP in post-Golgi and endosomal trafficking. FEBS Lett. 2009, 583 (23): 3817-3826. 10.1016/j.febslet.2009.10.026.View ArticlePubMedGoogle Scholar
- Jena BP: Role of SNAREs in membrane fusion. Adv Exp Med Biol. 2011, 713: 13-32. 10.1007/978-94-007-0763-4_3.View ArticlePubMedGoogle Scholar
- Goel VK, Li X, Chen H, Liu SC, Chishti AH, Oh SS: Band 3 is a host receptor binding merozoite surface protein 1 during the Plasmodium falciparum invasion of erythrocytes. Proc Natl Acad Sci U S A. 2003, 100 (9): 5164-5169. 10.1073/pnas.0834959100.View ArticlePubMed CentralPubMedGoogle Scholar
- Bharara R, Singh S, Pattnaik P, Chitnis CE, Sharma A: Structural analogs of sialic acid interfere with the binding of erythrocyte binding antigen-175 to glycophorin A, an interaction crucial for erythrocyte invasion by Plasmodium falciparum. Mol Biochem Parasitol. 2004, 138 (1): 123-129. 10.1016/j.molbiopara.2004.07.012.View ArticlePubMedGoogle Scholar
- Florens L, Washburn MP, Raine JD, Anthony RM, Grainger M, Haynes JD, Moch JK, Muster N, Sacci JB, Tabb DL, Witney AA, Wolters D, Wu Y, Gardner MJ, Holder AA, Sinden RE, Yates JR, Carucci DJ: A proteomic view of the Plasmodium falciparum life cycle. Nature. 2002, 419 (6906): 520-526. 10.1038/nature01107.View ArticlePubMedGoogle Scholar
- Mayer DC, Cofie J, Jiang L, Hartl DL, Tracy E, Kabat J, Mendoza LH, Miller LH: Glycophorin B is the erythrocyte receptor of Plasmodium falciparum erythrocyte-binding ligand, EBL-1. Proc Natl Acad Sci U S A. 2009, 106 (13): 5348-5352. 10.1073/pnas.0900878106.View ArticlePubMed CentralPubMedGoogle Scholar
- Maier AG, Baum J, Smith B, Conway DJ, Cowman AF: Polymorphisms in erythrocyte binding antigens 140 and 181 affect function and binding but not receptor specificity in Plasmodium falciparum. Infect Immun. 2009, 77 (4): 1689-1699. 10.1128/IAI.01331-08.View ArticlePubMed CentralPubMedGoogle Scholar
- Wanaguru M, Liu W, Hahn BH, Rayner JC, Wright GJ: RH5-Basigin interaction plays a major role in the host tropism of Plasmodium falciparum. Proc Natl Acad Sci U S A. 2013, 110 (51): 20735-20740. 10.1073/pnas.1320771110.View ArticlePubMed CentralPubMedGoogle Scholar
- Bartholdson SJ, Bustamante LY, Crosnier C, Johnson S, Lea S, Rayner JC, Wright GJ: Semaphorin-7A is an erythrocyte receptor for P. falciparum merozoite-specific TRAP homolog, MTRAP. PLoS Pathog. 2012, 8 (11): e1003031-10.1371/journal.ppat.1003031.View ArticlePubMed CentralPubMedGoogle Scholar
- Xu Y, Hu W, Chang Z, Duanmu H, Zhang S, Li Z, Yu L, Li X: Prediction of human protein-protein interaction by a mixed Bayesian model and its application to exploring underlying cancer-related pathway crosstalk. J R Soc Interface. 2011, 8 (57): 555-567. 10.1098/rsif.2010.0384.View ArticlePubMed CentralPubMedGoogle Scholar
- Abecasis GR, Altshuler D, Auton A, Brooks LD, Durbin RM, Gibbs RA, Hurles ME, McVean GA: A map of human genome variation from population-scale sequencing. Nature. 2010, 467 (7319): 1061-1073. 10.1038/nature09534.View ArticlePubMedGoogle Scholar
- Crosnier C, Wanaguru M, McDade B, Osier FH, Marsh K, Rayner JC, Wright GJ: A library of functional recombinant cell-surface and secreted P. falciparum merozoite proteins. Mol Cell Proteomics. 2013, 12 (12): 3976-3986. 10.1074/mcp.O113.028357.View ArticlePubMed CentralPubMedGoogle Scholar
- Zenonos ZA, Rayner JC, Wright GJ: Towards a comprehensive Plasmodium falciparum merozoite cell surface and secreted recombinant protein library. Malar J. 2014, 13: 93-10.1186/1475-2875-13-93.View ArticlePubMed CentralPubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.