Algorithmic Mapping and Characterization of the Drug-Induced Phenotypic-Response Space of Parasites Causing Schistosomiasis.

Neglected tropical diseases, especially those caused by helminths, constitute some of the most common infections of the world's poorest people. Amongst these, schistosomiasis (bilharzia or 'snail fever'), caused by blood flukes of the genus Schistosoma, ranks second only to malaria in terms of human impact: two hundred million people are infected and close to 800 million are at risk of infection. Drug screening against helminths poses unique challenges: the parasite cannot be cloned and is difficult to target using gene knockouts or RNAi. Consequently, both lead identification and validation involve phenotypic screening, where parasites are exposed to compounds whose effects are determined through the analysis of the ensuing phenotypic responses. The efficacy of leads thus identified derives from one or more or even unknown molecular mechanisms of action. The two most immediate and significant challenges that confront the state-of-the-art in this area are: the development of automated and quantitative phenotypic screening techniques and the mapping and quantitative characterization of the totality of phenotypic responses of the parasite. In this paper, we investigate and propose solutions for the latter problem in terms of the following: (1) mathematical formulation and algorithms that allow rigorous representation of the phenotypic response space of the parasite, (2) application of graph-theoretic and network analysis techniques for quantitative modeling and characterization of the phenotypic space, and (3) application of the aforementioned methodology to analyze the phenotypic space of S. mansoni - one of the etiological agents of schistosomiasis, induced by compounds that target its polo-like kinase 1 (PLK 1) gene - a recently validated drug target. In our approach, first, bio-image analysis algorithms are used to quantify the phenotypic responses of different drugs. Next, these responses are linearly mapped into a low- dimensional space using Principle Component Analysis (PCA). The phenotype space is modeled using neighborhood graphs which are used to represent the similarity amongst the phenotypes. These graphs are characterized and explored using network analysis algorithms. We present a number of results related to both the nature of the phenotypic space of the S. mansoni parasite as well as algorithmic issues encountered in constructing and analyzing the phenotypic-response space. In particular, the phenotype distribution of the parasite was found to have a distinct shape and topology. We have also quantitatively characterized the phenotypic space by varying critical model parameters. Finally, these maps of the phenotype space allows visualization and reasoning about complex relationships between putative drugs and their system-wide effects and can serve as a highly efficient paradigm for assimilating and unifying information from phenotypic screens both during lead identification and lead optimization.

Structure-Bioactivity Relationship for Benzimidazole Thiophene Inhibitors of Polo-Like Kinase 1 (PLK1), a Potential Drug Target in Schistosoma mansoni.

BACKGROUND: Schistosoma flatworm parasites cause schistosomiasis, a chronic and debilitating disease of poverty in developing countries. Praziquantel is employed for treatment and disease control. However, its efficacy spectrum is incomplete (less active or inactive against immature stages of the parasite) and there is a concern of drug resistance. Thus, there is a need to identify new drugs and drug targets. METHODOLOGY/PRINCIPAL FINDINGS: We show that RNA interference (RNAi) of the Schistosoma mansoni ortholog of human polo-like kinase (huPLK)1 elicits a deleterious phenotypic alteration in post-infective larvae (schistosomula or somules). Phenotypic screening and analysis of schistosomula and adult S. mansoni with small molecule inhibitors of huPLK1 identified a number of potent anti-schistosomals. Among these was a GlaxoSmithKline (GSK) benzimidazole thiophene inhibitor that has completed Phase I clinical trials for treatment of solid tumor malignancies. We then obtained GSKs Published Kinase Inhibitor Sets (PKIS) 1 and 2, and phenotypically screened an expanded series of 38 benzimidazole thiophene PLK1 inhibitors. Computational analysis of controls and PLK1 inhibitor-treated populations of somules demonstrated a distinctive phenotype distribution. Using principal component analysis (PCA), the phenotypes exhibited by these populations were mapped, visualized and analyzed through projection to a low-dimensional space. The phenotype distribution was found to have a distinct shape and topology, which could be elicited using cluster analysis. A structure-activity relationship (SAR) was identified for the benzimidazole thiophenes that held for both somules and adult parasites. The most potent inhibitors produced marked phenotypic alterations at 1-2 µM within 1 h. Among these were compounds previously characterized as potent inhibitors of huPLK1 in cell assays. CONCLUSIONS/SIGNIFICANCE: The reverse genetic and chemical SAR data support a continued investigation of SmPLK1 as a possible drug target and/or the prosecution of the benzimidazole thiophene chemotype as a source of novel anti-schistosomals.