Proteomic Analysis of Exosomes from Mutant KRAS Colon Cancer Cells Identifies Intercellular Transfer of Mutant KRAS

Activating mutations in KRAS occur in 30% to 40% of colorectal cancers. How mutant KRAS alters cancer cell behavior has been studied intensively, but non-cell autonomous effects of mutant KRAS are less understood. We recently reported that exosomes isolated from mutant KRAS-expressing colon cancer cells enhanced the invasiveness of recipient cells relative to exosomes purified from wild-type KRAS-expressing cells, leading us to hypothesize mutant KRAS might affect neighboring and distant cells by regulating exosome composition and behavior. Herein, we show the results of a comprehensive proteomic analysis of exosomes from parental DLD-1 cells that contain both wild-type and G13D mutant KRAS alleles and isogenically matched derivative cell lines, DKO-1 (mutant KRAS allele only) and DKs-8 (wild-type KRAS allele only). Mutant KRAS status dramatically affects the composition of the exosome proteome. Exosomes from mutant KRAS cells contain many tumor-promoting proteins, including KRAS, EGFR, SRC family kinases, and integrins. DKs-8 cells internalize DKO-1 exosomes, and, notably, DKO-1 exosomes transfer mutant KRAS to DKs-8 cells, leading to enhanced three-dimensional growth of these wild-type KRAS-expressing non-transformed cells. These results have important implications for non-cell autonomous effects of mutant KRAS, such as field effect and tumor progression.K-RAS (KRAS) is a small, monomeric GTPase whose biological activity is specified by its nucleotide binding state. Multiple lines of evidence highlight the importance of KRAS in colorectal cancer (CRC).¹ For example, activating missense mutations in KRAS, which lock the protein into the GTP-bound state, occur in 30% to 40% of CRCs and are strongly associated with poor prognosis (1, 2). Also, mutant KRAS negatively predicts responsiveness to anti-EGF receptor (EGFR) therapy (3).Early attempts to decipher the neoplastic consequences of mutant KRAS relied on overexpression studies. A drawback of these studies is their failure to simulate the genetic conditions present in human tumors, where there is often one wild-type (WT) and one mutant KRAS allele (1). More recently, KRAS mutant CRC cell lines have been engineered to selectively contain either the wild-type or the mutant KRAS allele (4), and a single mutant Kras allele has been activated in the intestine using genetically engineered mice (5). Detailed studies using these complementary approaches demonstrate a wide range of tumor-promoting effects of mutant KRAS (reviewed in Ref. 6). Much of what is known about mutant KRAS pertains to its ability to alter the behavior of a transformed cell in a cell autonomous manner. With the exception of increased tumor vascularity via increased tumor-derived VEGF expression (7, 8), non-cell autonomous effects of mutant KRAS have been much less studied.Exosomes are 30- to 100-nm secreted vesicles that have emerged as a novel mode of intercellular communication (9). We recently reported that exosomes purified from conditioned medium of mutant KRAS CRC cells contained higher levels of the EGFR ligand amphiregulin (AREG) and enhanced invasiveness of recipient cancer cells relative to exosomes from isogenically matched wild-type KRAS cells (10). These results prompted us to perform a comprehensive analysis of exosomes purified from these cells. Herein, we show that mutant KRAS induces many changes in exosomal protein composition. Notably, we show that (i) KRAS is contained within exosomes, (ii) exosomes can transfer mutant KRAS to cells expressing only wild-type KRAS, and (iii) mutant KRAS-containing exosomes enhance wild-type KRAS cell growth in collagen matrix and soft agar. These results have important implications for the progression of CRC tumors by providing a mechanism by which the tumor microenvironment may be influenced by non-cell autonomous signals released by mutant KRAS-expressing tumor cells.     

Proteomic and Genetic Analyses Demonstrate that Plasmodium berghei Blood Stages Export a Large and Diverse Repertoire of Proteins

Malaria parasites actively remodel the infected red blood cell (irbc) by exporting proteins into the host cell cytoplasm. The human parasite Plasmodium falciparum exports particularly large numbers of proteins, including proteins that establish a vesicular network allowing the trafficking of proteins onto the surface of irbcs that are responsible for tissue sequestration. Like P. falciparum, the rodent parasite P. berghei ANKA sequesters via irbc interactions with the host receptor CD36. We have applied proteomic, genomic, and reverse-genetic approaches to identify P. berghei proteins potentially involved in the transport of proteins to the irbc surface. A comparative proteomics analysis of P. berghei non-sequestering and sequestering parasites was used to determine changes in the irbc membrane associated with sequestration. Subsequent tagging experiments identified 13 proteins (Plasmodium export element (PEXEL)-positive as well as PEXEL-negative) that are exported into the irbc cytoplasm and have distinct localization patterns: a dispersed and/or patchy distribution, a punctate vesicle-like pattern in the cytoplasm, or a distinct location at the irbc membrane. Members of the PEXEL-negative BIR and PEXEL-positive Pb-fam-3 show a dispersed localization in the irbc cytoplasm, but not at the irbc surface. Two of the identified exported proteins are transported to the irbc membrane and were named erythrocyte membrane associated proteins. EMAP1 is a member of the PEXEL-negative Pb-fam-1 family, and EMAP2 is a PEXEL-positive protein encoded by a single copy gene; neither protein plays a direct role in sequestration. Our observations clearly indicate that P. berghei traffics a diverse range of proteins to different cellular locations via mechanisms that are analogous to those employed by P. falciparum. This information can be exploited to generate transgenic humanized rodent P. berghei parasites expressing chimeric P. berghei/P. falciparum proteins on the surface of rodent irbc, thereby opening new avenues for in vivo screening adjunct therapies that block sequestration.Malaria parasites invade and develop inside red blood cells, and extensive remodeling of the host cell is required in order for the parasite to take up nutrients and grow (1). In addition, infected red blood cells (irbcs)¹ of several Plasmodium species adhere to endothelium lining blood capillaries, and this is achieved through modification of the irbc, specifically, alteration of the irbc membrane (2, 3). This active remodeling of the erythrocyte requires the export of parasite proteins into the host cell cytoplasm and their incorporation in the irbc membrane of the host cell (1, 2). Schizont-infected red blood cells of the rodent parasite P. berghei ANKA adhere to endothelial cells of the microvasculature, leading to the sequestration of irbcs in organs such as the lungs and adipose tissue (4–6). P. berghei irbcs adhere to the class II scavenger receptor CD36 (7), which is highly conserved in mammals and is the receptor most commonly used by irbcs infected with the human parasite P. falciparum (8). These observations suggest that P. berghei may export proteins onto the surface of irbcs in a fashion analogous to the processes employed by P. falciparum that expresses PfEMP1, the protein known to be responsible for P. falciparum irbc sequestration. However, P. berghei does not contain Pfemp1 orthologs or proteins with domains with clear homology to the domains of PfEMP1 (9), and the P. berghei proteins responsible for irbc cytoadherence and proteins involved in the transport of these proteins to the irbc membrane remain largely unknown. Recently we used a proteomic analysis of P. berghei ANKA irbc membranes to identify parasite proteins associated with the erythrocyte membrane, and we have demonstrated that the deletion of a single-copy gene of P. berghei that encodes a small exported protein known as SMAC results in strongly reduced irbc sequestration (6). No evidence was found for the presence of SMAC on the irbc surface, and therefore this protein is most likely involved in the transport or anchoring of other P. berghei proteins that directly interact with host receptors on endothelial cells.For P. falciparum, a large number of exported proteins have been predicted based on the presence of an N-terminal motif known as the Plasmodium export element (PEXEL) motif (10, 11). Many of these PEXEL-positive proteins belong to species-specific gene families. Comparison of PEXEL-positive proteins in different Plasmodium species suggested that P. falciparum expresses a significantly higher number of exported proteins than other Plasmodium species, which in part can be attributed to the expansion of P. falciparum–specific protein families, including those containing DnaJ or PHIST domains (12–17). One explanation for the elevated number of exported proteins in P. falciparum is that they are necessary for export of the P. falciparum–specific protein PfEMP1 to the irbc surface (10). Comparisons of different Plasmodium exportomes have mainly focused on identifying orthologs of the PEXEL-positive proteins of P. falciparum in the other species (14, 15, 18). For example, of the >500 PEXEL-positive P. falciparum proteins, only between 11 and 33 had orthologs in P. berghei (14, 15, 19). However, such an approach might underestimate the total number of exported proteins. A recent hidden Markov model (HMM) analysis of the PEXEL motif for P. berghei proteins identified at least 75 PEXEL-positive P. berghei proteins (6). Moreover, in different Plasmodium species, a number of exported proteins have been described that are PEXEL-negative, indicating that alternative export pathways might exist that are independent of the presence of a PEXEL motif (20, 21). It has been suggested that in species with a small number of PEXEL-positive proteins, PEXEL-negative exported proteins play a more prominent role in host cell remodeling (21). An example of a PEXEL-negative exported protein family is the large PIR family of proteins, which are expressed by rodent Plasmodium species (9, 22), the monkey parasite P. knowlesi (23), and the human parasite P. vivax (24, 25).To date, export to the irbc cytosol has been shown for only a few P. berghei proteins (i.e. several members of the BIR family; TIGR01590) (6), two members of the ETRAMP family (26), and two proteins encoded by a single copy gene, SMAC and IBIS1 (6, 27). In this study, comparative proteomic, genomic, and reverse-genetic approaches have been used to identify novel exported proteins of P. berghei. We report proteome analyses of samples enriched for proteins associated with membranes of irbcs from both sequestering P. berghei ANKA and non-sequestering P. berghei K173 parasites, and we also present analyses of the full genome sequence of a non-sequestering P. berghei K173 line. Fluorescent tagging of parasite proteins selected from the proteome and genome analyses identified a number of novel P. berghei ANKA proteins that are exported into the irbc cytoplasm. We report for the first time the export of members of the PEXEL-negative Pb-fam-1 gene family (pyst-a; TIGR01599) and show that two proteins are transported to the P. berghei ANKA irbc membrane. This is the first comprehensive study of exported proteins of P. berghei that has been validated via the generation of a large number of transgenic P. berghei ANKA parasites expressing tagged proteins and has shown the export of both PEXEL-positive and PEXEL-negative proteins to the irbc cytoplasm. The identification of P. berghei ANKA proteins exported to the irbc membrane and proteins involved in sequestration suggests the possibility of developing “humanized” small animal models for the in vivo analysis of the sequestration properties of P. falciparum proteins that would express (domains of) P. falciparum proteins on the surface of rodent irbcs (4, 6).     

Proteomic Analysis of the Relationship between Metabolism and Nonhost Resistance in Soybean Exposed to Bipolaris maydis

Nonhost resistance (NHR) pertains to the most common form of plant resistance against pathogenic microorganisms of other species. Bipolaris maydis is a non-adapted pathogen affecting soybeans, particularly of maize/soybean intercropping systems. However, no experimental evidence has described the immune response of soybeans against B. maydis. To elucidate the molecular mechanism underlying NHR in soybeans, proteomics analysis based on two-dimensional polyacrylamide gel electrophoresis (2-DE) was performed to identify proteins involved in the soybean response to B. maydis. The spread of B. maydis spores across soybean leaves induced NHR throughout the plant, which mobilized almost all organelles and various metabolic processes in response to B. maydis. Some enzymes, including ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO), mitochondrial processing peptidase (MPP), oxygen evolving enhancer (OEE), and nucleoside diphosphate kinase (NDKs), were found to be related to NHR in soybeans. These enzymes have been identified in previous studies, and STRING analysis showed that most of the protein functions related to major metabolic processes were induced as a response to B. maydis, which suggested an array of complex interactions between soybeans and B. maydis. These findings suggest a systematic NHR against non-adapted pathogens in soybeans. This response was characterized by an overlap between metabolic processes and response to stimulus. Several metabolic processes provide the soybean with innate immunity to the non-adapted pathogen, B. maydis. This research investigation on NHR in soybeans may foster a better understanding of plant innate immunity, as well as the interactions between plant and non-adapted pathogens in intercropping systems.     

大豆种子萌发过程中的差异蛋白质组研究

运用蛋白质组学技术对大豆(Glycinemax)N2899种子萌发0h、8h、36h、60h4个时期蛋白质的差异表达情况进行了研究.结果发现,在考马斯亮蓝染色的双向电泳pH3～10胶上,PDQuest图像分析软件可识别的点约350个,其中表达量变化2.5倍以上的蛋白质点有24个,而绝大部分大豆种子贮藏蛋白在萌发期尚未降解.在萌发的第一阶段,24个差异表达蛋白中有10个蛋白质的丰度发生变化.第二阶段,差异表达蛋白的种类和量增加,其中15个蛋白质是动态变化的,14个蛋白质在胚根突破种皮时表达量达到峰值,表明吸胀后种子内的生命活动越来越强.对这24个蛋白质点进行胶内酶解,用基质辅助激光解析电离飞行时间质谱测定均获得肽质量指纹图谱.搜索大豆的UniGene库初步鉴定出6个蛋白质,分别是核苷二磷酸激酶、热激蛋白、硫氧还蛋白、35ku种子成熟蛋白及种子成熟蛋白PM36.对这些蛋白质在种子萌发过程中可能的作用进行了讨论.     

Proteomic Analysis of Kunitz-Type Trypsin Inhibitor Deleted Soybean

Russian Journal of Plant Physiology - To explore proteomic characters of Kunitz-type trypsin inhibitors (KTIs) deleted soybean (Glycine max (L.) Merr.), seeds without KTIs and its female parent...     

The Temperature Dependent Proteomic Analysis of Thermotoga maritima

Thermotoga maritima (T. maritima) is a typical thermophile, and its proteome response to environmental temperature changes has yet to be explored. This study aims to uncover the temperature-dependent proteins of T. maritima using comparative proteomic approach. T. maritima was cultured under four temperatures, 60°C, 70°C, 80°C and 90°C, and the bacterial proteins were extracted and electrophoresed in two-dimensional mode. After analysis of gel images, a total of 224 spots, either cytoplasm or membrane, were defined as temperature-dependent. Of these spots, 75 unique bacterial proteins were identified using MALDI TOF/TOF MS. As is well known, the chaperone proteins such as heat shock protein 60 and elongation factor Tu, were up-regulated in abundance due to increased temperature. However, several temperature-dependent proteins of T. maritima responded very differently when compared to responses of the thermophile T. tengcongensis. Intriguingly, a number of proteins involved in central carbohydrate metabolism were significantly up-regulated at higher temperature. Their corresponding mRNA levels were elevated accordingly. The increase in abundance of several key enzymes indicates that a number of central carbohydrate metabolism pathways of T. maritima are activated at higher temperatures.     

Proteomic Analysis Identifies Ribosome Reduction as an Effective Proteotoxic Stress Response

Stress responses are adaptive cellular programs that identify and mitigate potentially dangerous threats. Misfolded proteins are a ubiquitous and clinically relevant stress. Trivalent metalloids, such as arsenic, have been proposed to cause protein misfolding. Using tandem mass tag-based mass spectrometry, we show that trivalent arsenic results in widespread reorganization of the cell from an anabolic to a catabolic state. Both major pathways of protein degradation, the proteasome and autophagy, show increased abundance of pathway components and increased functional output, and are required for survival. Remarkably, cells also showed a down-regulation of ribosomes at the protein level. That this represented an adaptive response and not an adverse toxic effect was indicated by enhanced survival of ribosome mutants after arsenic exposure. These results suggest that a major source of toxicity of trivalent arsenic derives from misfolding of newly synthesized proteins and identifies ribosome reduction as a rapid, effective, and reversible proteotoxic stress response.     

Dose-Dependent Regulation of Alternative Splicing by MBNL Proteins Reveals Biomarkers for Myotonic Dystrophy

Alternative splicing is a regulated process that results in expression ofspecific mRNA and protein isoforms. Alternative splicing factors determine therelative abundance of each isoform. Here we focus on MBNL1, a splicing factormisregulated in the disease myotonic dystrophy. By altering the concentration ofMBNL1 in cells across a broad dynamic range, we show that different splicingevents require different amounts of MBNL1 for half-maximal response, and respondmore or less steeply to MBNL1. Motifs around MBNL1 exon 5 were studied to assesshow cis-elements mediate the MBNL1 dose-dependent splicingresponse. A framework was developed to estimate MBNL concentration usingsplicing responses alone, validated in the cell-based model, and applied tomyotonic dystrophy patient muscle. Using this framework, we evaluated theability of individual and combinations of splicing events to predict functionalMBNL concentration in human biopsies, as well as their performance as biomarkersto assay mild, moderate, and severe cases of DM.     

Mechanism-based Proteomic Screening Identifies Targets of Thioredoxin-like Proteins

Thioredoxin (Trx)-fold proteins are protagonists of numerous cellular pathways that are subject to thiol-based redox control. The best characterized regulator of thiols in proteins is Trx1 itself, which together with thioredoxin reductase 1 (TR1) and peroxiredoxins (Prxs) comprises a key redox regulatory system in mammalian cells. However, there are numerous other Trx-like proteins, whose functions and redox interactors are unknown. It is also unclear if the principles of Trx1-based redox control apply to these proteins. Here, we employed a proteomic strategy to four Trx-like proteins containing CXXC motifs, namely Trx1, Rdx12, Trx-like protein 1 (Txnl1) and nucleoredoxin 1 (Nrx1), whose cellular targets were trapped in vivo using mutant Trx-like proteins, under conditions of low endogenous expression of these proteins. Prxs were detected as key redox targets of Trx1, but this approach also supported the detection of TR1, which is the Trx1 reductant, as well as mitochondrial intermembrane proteins AIF and Mia40. In addition, glutathione peroxidase 4 was found to be a Rdx12 redox target. In contrast, no redox targets of Txnl1 and Nrx1 could be detected, suggesting that their CXXC motifs do not engage in mixed disulfides with cellular proteins. For some Trx-like proteins, the method allowed distinguishing redox and non-redox interactions. Parallel, comparative analyses of multiple thiol oxidoreductases revealed differences in the functions of their CXXC motifs, providing important insights into thiol-based redox control of cellular processes.     

Comparative Proteomic Analysis Identifies Key Metabolic Regulators of Gemcitabine Resistance in Pancreatic Cancer

Pancreatic adenocarcinoma (PDAC) is highly refractory to treatment. Standard-of-care gemcitabine (Gem) provides only modest survival benefits, and development of Gem resistance (GemR) compromises its efficacy. Highly GemR clones of Gem-sensitive MIAPaCa-2 cells were developed to investigate the molecular mechanisms of GemR and implemented global quantitative differential proteomics analysis with a comprehensive, reproducible ion-current–based MS1 workflow to quantify ～6000 proteins in all samples. In GemR clone MIA-GR8, cellular metabolism, proliferation, migration, and ‘drug response’ mechanisms were the predominant biological processes altered, consistent with cell phenotypic alterations in cell cycle and motility. S100 calcium binding protein A4 was the most downregulated protein, as were proteins associated with glycolytic and oxidative energy production. Both responses would reduce tumor proliferation. Upregulation of mesenchymal markers was prominent, and cellular invasiveness increased. Key enzymes in Gem metabolism pathways were altered such that intracellular utilization of Gem would decrease. Ribonucleoside-diphosphate reductase large subunit was the most elevated Gem metabolizing protein, supporting its critical role in GemR. Lower Ribonucleoside-diphosphate reductase large subunit expression is associated with better clinical outcomes in PDAC, and its downregulation paralleled reduced MIAPaCa-2 proliferation and migration and increased Gem sensitivity. Temporal protein-level Gem responses of MIAPaCa-2 versus GemR cell lines (intrinsically GemR PANC-1 and acquired GemR MIA-GR8) implicate adaptive changes in cellular response systems for cell proliferation and drug transport and metabolism, which reduce cytotoxic Gem metabolites, in DNA repair, and additional responses, as key contributors to the complexity of GemR in PDAC. These findings additionally suggest targetable therapeutic vulnerabilities for GemR PDAC patients.     

Diversification of Two Lineages of Symbiotic Photobacterium

Understanding of processes driving bacterial speciation requires examination of closely related, recently diversified lineages. To gain an insight into diversification of bacteria, we conducted comparative genomic analysis of two lineages of bioluminescent symbionts, Photobacterium leiognathi and ‘P. mandapamensis’. The two lineages are evolutionary and ecologically closely related. Based on the methods used in bacterial taxonomy for classification of new species (DNA-DNA hybridization and ANI), genetic relatedness of the two lineages is at a cut-off point for species delineation. In this study, we obtained the whole genome sequence of a representative P. leiognathi strain lrivu.4.1, and compared it to the whole genome sequence of ‘P. mandapamensis’ svers.1.1. Results of the comparative genomic analysis suggest that P. leiognathi has a more plastic genome and acquired genes horizontally more frequently than ‘P. mandapamensis’. We predict that different rates of recombination and gene acquisition contributed to diversification of the two lineages. Analysis of lineage-specific sequences in 25 strains of P. leiognathi and ‘P. mandapamensis’ found no evidence that bioluminescent symbioses with specific host animals have played a role in diversification of the two lineages.     

Different Epidermal Growth Factor Receptor (EGFR) Agonists Produce Unique Signatures for the Recruitment of Downstream Signaling Proteins

The EGF receptor can bind seven different agonist ligands. Although each agonist appears to stimulate the same suite of downstream signaling proteins, different agonists are capable of inducing distinct responses in the same cell. To determine the basis for these differences, we used luciferase fragment complementation imaging to monitor the recruitment of Cbl, CrkL, Gab1, Grb2, PI3K, p52 Shc, p66 Shc, and Shp2 to the EGF receptor when stimulated by the seven EGF receptor ligands. Recruitment of all eight proteins was rapid, dose-dependent, and inhibited by erlotinib and lapatinib, although to differing extents. Comparison of the time course of recruitment of the eight proteins in response to a fixed concentration of each growth factor revealed differences among the growth factors that could contribute to their differing biological effects. Principal component analysis of the resulting data set confirmed that the recruitment of these proteins differed between agonists and also between different doses of the same agonist. Ensemble clustering of the overall response to the different growth factors suggests that these EGF receptor ligands fall into two major groups as follows: (i) EGF, amphiregulin, and EPR; and (ii) betacellulin, TGFα, and epigen. Heparin-binding EGF is distantly related to both clusters. Our data identify differences in network utilization by different EGF receptor agonists and highlight the need to characterize network interactions under conditions other than high dose EGF.     

The Ubiquitous Distribution of Late Embryogenesis Abundant Proteins across Cell Compartments in Arabidopsis Offers Tailored Protection against Abiotic Stress

Late embryogenesis abundant (LEA) proteins are hydrophilic, mostly intrinsically disordered proteins, which play major roles in desiccation tolerance. In Arabidopsis thaliana, 51 genes encoding LEA proteins clustered into nine families have been inventoried. To increase our understanding of the yet enigmatic functions of these gene families, we report the subcellular location of each protein. Experimental data highlight the limits of in silico predictions for analysis of subcellular localization. Thirty-six LEA proteins localized to the cytosol, with most being able to diffuse into the nucleus. Three proteins were exclusively localized in plastids or mitochondria, while two others were found dually targeted to these organelles. Targeting cleavage sites could be determined for five of these proteins. Three proteins were found to be endoplasmic reticulum (ER) residents, two were vacuolar, and two were secreted. A single protein was identified in pexophagosomes. While most LEA protein families have a unique subcellular localization, members of the LEA_4 family are widely distributed (cytosol, mitochondria, plastid, ER, and pexophagosome) but share the presence of the class A α-helix motif. They are thus expected to establish interactions with various cellular membranes under stress conditions. The broad subcellular distribution of LEA proteins highlights the requirement for each cellular compartment to be provided with protective mechanisms to cope with desiccation or cold stress.     

In Silico Repositioning-Chemogenomics Strategy Identifies New Drugs with Potential Activity against Multiple Life Stages of Schistosoma mansoni

Morbidity and mortality caused by schistosomiasis are serious public health problems in developing countries. Because praziquantel is the only drug in therapeutic use, the risk of drug resistance is a concern. In the search for new schistosomicidal drugs, we performed a target-based chemogenomics screen of a dataset of 2,114 proteins to identify drugs that are approved for clinical use in humans that may be active against multiple life stages of Schistosoma mansoni. Each of these proteins was treated as a potential drug target, and its amino acid sequence was used to interrogate three databases: Therapeutic Target Database (TTD), DrugBank and STITCH. Predicted drug-target interactions were refined using a combination of approaches, including pairwise alignment, conservation state of functional regions and chemical space analysis. To validate our strategy, several drugs previously shown to be active against Schistosoma species were correctly predicted, such as clonazepam, auranofin, nifedipine, and artesunate. We were also able to identify 115 drugs that have not yet been experimentally tested against schistosomes and that require further assessment. Some examples are aprindine, gentamicin, clotrimazole, tetrabenazine, griseofulvin, and cinnarizine. In conclusion, we have developed a systematic and focused computer-aided approach to propose approved drugs that may warrant testing and/or serve as lead compounds for the design of new drugs against schistosomes.