Polycystic kidney disease (PKD) is a prevalent disorder characterized by renal cysts that lead to kidney failure. Various signaling pathways have been targeted to stop disease progression, but most... Show morePolycystic kidney disease (PKD) is a prevalent disorder characterized by renal cysts that lead to kidney failure. Various signaling pathways have been targeted to stop disease progression, but most interventions still focus on alleviating PKD-associated symptoms. The mechanistic complexity of the disease, as well as the lack of functional in vitro assays for compound testing, has made drug discovery for PKD challenging. To identify modulators of PKD, Pkd1–/– kidney tubule epithelial cells were applied to a scalable and automated 3D cyst culture model for compound screening, followed by phenotypic profiling to determine compound efficacy. We used this screening platform to screen a library of 273 kinase inhibitors to probe various signaling pathways involved in cyst growth. We show that inhibition of several targets, including aurora kinase, CDK, Chk, IGF-1R, Syk, and mTOR, but, surprisingly, not PI3K, prevented forskolin-induced cyst swelling. Additionally, we show that multiparametric phenotypic classification discriminated potentially undesirable (i.e., cytotoxic) compounds from molecules inducing the desired phenotypic change, greatly facilitating hit selection and validation. Our findings show that a pathophysiologically relevant 3D cyst culture model of PKD coupled to phenotypic profiling can be used to identify potentially therapeutic compounds and predict and validate molecular targets for PKD. Show less
Human and animal studies have converged to suggest that caffeine consumption prevents memory deficits in aging and Alzheimer’s disease through the antagonism of adenosine A2A receptors (A2ARs). To... Show moreHuman and animal studies have converged to suggest that caffeine consumption prevents memory deficits in aging and Alzheimer’s disease through the antagonism of adenosine A2A receptors (A2ARs). To test if A2AR activation in the hippocampus is actually sufficient to impair memory function and to begin elucidating the intracellular pathways operated by A2AR, we have developed a chimeric rhodopsin-A2AR protein (optoA2AR), which retains the extracellular and transmembrane domains of rhodopsin (conferring light responsiveness and eliminating adenosine-binding pockets) fused to the intracellular loop of A2AR to confer specific A2AR signaling. The specificity of the optoA2AR signaling was confirmed by light-induced selective enhancement of cAMP and phospho-mitogen-activated protein kinase (p-MAPK) (but not cGMP) levels in human embryonic kidney 293 (HEK293) cells, which was abolished by a point mutation at the C terminal of A2AR. Supporting its physiological relevance, optoA2AR activation and the A2AR agonist CGS21680 produced similar activation of cAMP and p-MAPK signaling in HEK293 cells, of p-MAPK in the nucleus accumbens and of c-Fos/phosphorylated-CREB (p-CREB) in the hippocampus, and similarly enhanced long-term potentiation in the hippocampus. Remarkably, optoA2AR activation triggered a preferential p-CREB signaling in the hippocampus and impaired spatial memory performance, while optoA2AR activation in the nucleus accumbens triggered MAPK signaling and modulated locomotor activity. This shows that the recruitment of intracellular A2AR signaling in the hippocampus is sufficient to trigger memory dysfunction. Furthermore, the demonstration that the biased A2AR signaling and functions depend on intracellular A2AR loops prompts the possibility of targeting the intracellular A2AR-interacting partners to selectively control different neuropsychiatric behaviors. Show less
Cyclooxygenases (COX) are present in the body in two isoforms, namely: COX-1, constitutively expressed, and COX-2, induced in physiopathological conditions such as cancer or chronic inflammation.... Show moreCyclooxygenases (COX) are present in the body in two isoforms, namely: COX-1, constitutively expressed, and COX-2, induced in physiopathological conditions such as cancer or chronic inflammation. The inhibition of COX with non-steroideal anti-inflammatory drugs (NSAIDs) is the most widely used treatment for chronic inflammation despite the adverse effects associated to prolonged NSAIDs intake. Although selective COX-2 inhibition has been shown not to palliate all adverse effects (e.g.cardiotoxicity), there are still niche populations which can benefit from selective COX-2 inhibition. Thus, capitalizing on bioactivity data from both isoforms simultaneously would contribute to develop COX inhibitors with better safety profiles. We applied ensemble proteochemometric modeling (PCM) for the prediction of the potency of 3,228 distinct COX inhibitors on 11 mammalian cyclooxygenases. Ensemble PCM models (R20test=0.65R20test=0.65R20test=0.65, and RMSEtest= 0.71) outperformed models exclusively trained on compound (R20test=0.17R20test=0.17R20test=0.17, and RMSEtest= 1.09) or protein descriptors (R20test=0.16R20test=0.16R20test=0.16and RMSEtest= 1.10) on the test set. Moreover, PCM predicted COX potency for 1,086 selective and non-selective COX inhibitors withR20test=0.59R20test=0.59R20test=0.59and RMSEtest= 0.76. These values are in agreement with the maximum and minimum achievableR20testR20testR20testand RMSEtestvalues of approximately 0.68 for both metrics. Confidence intervals for individual predictions were calculated from the standard deviation of the predictions from the individual models composing the ensembles. Finally, two substructure analysis pipelines singled out chemical substructures implicated in both potency and selectivity in agreement with the literature. Show less
Proteochemometric (PCM) modelling is a computational method to model the bioactivity of multiple ligands against multiple related protein targets simultaneously. Hence it has been found to be... Show moreProteochemometric (PCM) modelling is a computational method to model the bioactivity of multiple ligands against multiple related protein targets simultaneously. Hence it has been found to be particularly useful when exploring the selectivity and promiscuity of ligands on different proteins. In this review, we will firstly provide a brief introduction to the main concepts of PCM for readers new to the field. The next part focuses on recent technical advances, including the application of support vector machines (SVMs) using different kernel functions, random forests, Gaussian processes and collaborative filtering. The subsequent section will then describe some novel practical applications of PCM in the medicinal chemistry field, including studies on GPCRs, kinases, viral proteins (e.g.from HIV) and epigenetic targets such as histone deacetylases. Finally, we will conclude by summarizing novel developments in PCM, which we expect to gain further importance in the future. These developments include adding three-dimensional protein target information, application of PCM to the prediction of binding energies, and application of the concept in the fields of pharmacogenomics and toxicogenomics. This review is an update to a related publication in 2011 and it mainly focuses on developments in the field since then. Show less
Proteochemometrics (PCM) is an approach for bioactivity predictive modeling which models the relationship between protein and chemical information. Gaussian Processes (GP), based on Bayesian... Show moreProteochemometrics (PCM) is an approach for bioactivity predictive modeling which models the relationship between protein and chemical information. Gaussian Processes (GP), based on Bayesian inference, provide the most objective estimation of the uncertainty of the predictions, thus permitting the evaluation of the applicability domain (AD) of the model. Furthermore, the experimental error on bioactivity measurements can be used as input for this probabilistic model. In this study, we apply GP implemented with a panel of kernels on three various (and multispecies) PCM datasets. The first dataset consisted of information from 8 human and rat adenosine receptors with 10,999 small molecule ligands and their binding affinity. The second consisted of the catalytic activity of four dengue virus NS3 proteases on 56 small peptides. Finally, we have gathered bioactivity information of small molecule ligands on 91 aminergic GPCRs from 9 different species, leading to a dataset of 24,593 datapoints with a matrix completeness of only 2.43%. GP models trained on these datasets are statistically sound, at the same level of statistical significance as Support Vector Machines (SVM), with R20R20R20 values on the external dataset ranging from 0.68 to 0.92, and RMSEP values close to the experimental error. Furthermore, the best GP models obtained with the normalized polynomial and radial kernels provide intervals of confidence for the predictions in agreement with the cumulative Gaussian distribution. GP models were also interpreted on the basis of individual targets and of ligand descriptors. In the dengue dataset, the model interpretation in terms of the amino-acid positions in the tetra-peptide ligands gave biologically meaningful results. Show less
Resistance to pesticides is an increasing problem in agriculture. Despite practices such as phased use and cycling of ‘orthogonally resistant’ agents, resistance remains a major risk to national... Show moreResistance to pesticides is an increasing problem in agriculture. Despite practices such as phased use and cycling of ‘orthogonally resistant’ agents, resistance remains a major risk to national and global food security. To combat this problem, there is a need for both new approaches for pesticide design, as well as for novel chemical entities themselves. As summarized in this opinion article, a technique termed ‘proteochemometric modelling’ (PCM), from the field of chemoinformatics, could aid in the quantification and prediction of resistance that acts via point mutations in the target proteins of an agent. The technique combines information from both the chemical and biological domain to generate bioactivity models across large numbers of ligands as well as protein targets. PCM has previously been validated in prospective, experimental work in the medicinal chemistry area, and it draws on the growing amount of bioactivity information available in the public domain. Here, two potential applications of proteochemometric modelling to agrochemical data are described, based on previously published examples from the medicinal chemistry literature. Show less
Allosteric modulators are ligands for proteins that exert their effects via a different binding site than the natural (orthosteric) ligand site and hence form a conceptually distinct class of... Show moreAllosteric modulators are ligands for proteins that exert their effects via a different binding site than the natural (orthosteric) ligand site and hence form a conceptually distinct class of ligands for a target of interest. Here, the physicochemical and structural features of a large set of allosteric and non-allosteric ligands from the ChEMBL database of bioactive molecules are analyzed. In general allosteric modulators are relatively smaller, more lipophilic and more rigid compounds, though large differences exist between different targets and target classes. Furthermore, there are differences in the distribution of targets that bind these allosteric modulators. Allosteric modulators are over-represented in membrane receptors, ligand-gated ion channels and nuclear receptor targets, but are underrepresented in enzymes (primarily proteases and kinases). Moreover, allosteric modulators tend to bind to their targets with a slightly lower potency (5.96 log units versus 6.66 log units, p<0.01). However, this lower absolute affinity is compensated by their lower molecular weight and more lipophilic nature, leading to similar binding efficiency and surface efficiency indices. Subsequently a series of classifier models are trained, initially target class independent models followed by finer-grained target (architecture/functional class) based models using the target hierarchy of the ChEMBL database. Applications of these insights include the selection of likely allosteric modulators from existing compound collections, the design of novel chemical libraries biased towards allosteric regulators and the selection of targets potentially likely to yield allosteric modulators on screening. All data sets used in the paper are available for download. Show less
Kufareva, I.; Katritch, V.; Westen, G.J.P. van; Lenselink, E.B.; Overington, J.P.; Participants, of GPCR Dock 2013; Stevens, RC Abagyan R 2014
Despite tremendous successes of GPCR crystallography, the receptors with available structures represent only a small fraction of human GPCRs. An important role of the modeling community is to... Show moreDespite tremendous successes of GPCR crystallography, the receptors with available structures represent only a small fraction of human GPCRs. An important role of the modeling community is to maximize structural insights for the remaining receptors and complexes. The community-wide GPCR Dock assessment was established to stimulate and monitor the progress in molecular modeling and ligand docking for GPCRs. The four targets in the present third assessment round presented new and diverse challenges for modelers, including prediction of allosteric ligand interaction and activation states in 5-hydroxytryptamine receptors 1B and 2B, and modeling by extremely distant homology for smoothened receptor. Forty-four modeling groups participated in the assessment. State-of-the-art modeling approaches achieved close-to-experimental accuracy for small rigid orthosteric ligands and models built by close homology, and they correctly predicted protein fold for distant homology targets. Predictions of long loops and GPCR activation states remain unsolved problems. Show less
Proteochemometric (PCM) modelling is a computational method to model the bioactivity of multiple ligands against multiple related protein targets simultaneously. Hence it has been found to be... Show moreProteochemometric (PCM) modelling is a computational method to model the bioactivity of multiple ligands against multiple related protein targets simultaneously. Hence it has been found to be particularly useful when exploring the selectivity and promiscuity of ligands on different proteins. In this review, we will firstly provide a brief introduction to the main concepts of PCM for readers new to the field. The next part focuses on recent technical advances, including the application of support vector machines (SVMs) using different kernel functions, random forests, Gaussian processes and collaborative filtering. The subsequent section will then describe some novel practical applications of PCM in the medicinal chemistry field, including studies on GPCRs, kinases, viral proteins (e.g.from HIV) and epigenetic targets such as histone deacetylases. Finally, we will conclude by summarizing novel developments in PCM, which we expect to gain further importance in the future. These developments include adding three-dimensional protein target information, application of PCM to the prediction of binding energies, and application of the concept in the fields of pharmacogenomics and toxicogenomics. This review is an update to a related publication in 2011 and it mainly focuses on developments in the field since then. Show less
Westen, G.J.P. van; Swier, R.F.; Wegner, J.K.; IJzerman, A.P.; Vlijmen, H. van; Bender, A. 2013
Background While a large body of work exists on comparing and benchmarking of descriptors of molecular structures, a similar comparison of protein descriptor sets is lacking. Hence, in the... Show moreBackground While a large body of work exists on comparing and benchmarking of descriptors of molecular structures, a similar comparison of protein descriptor sets is lacking. Hence, in the current work a total of 13 different protein descriptor sets have been compared with respect to their behavior in perceiving similarities between amino acids. The descriptor sets included in the study are Z-scales (3 variants), VHSE, T-scales, ST-scales, MS-WHIM, FASGAI and BLOSUM, and a novel protein descriptor set termed ProtFP (4 variants). We investigate to which extent descriptor sets show collinear as well as orthogonal behavior via principal component analysis (PCA). Results In describing amino acid similarities, MSWHIM, T-scales and ST-scales show related behavior, as do the VHSE, FASGAI, and ProtFP (PCA3) descriptor sets. Conversely, the ProtFP (PCA5), ProtFP (PCA8), Z-Scales (Binned), and BLOSUM descriptor sets show behavior that is distinct from one another as well as both of the clusters above. Generally, the use of more principal components (>3 per amino acid, per descriptor) leads to a significant differences in the way amino acids are described, despite that the later principal components capture less variation per component of the original input data. Conclusion In this work a comparison is provided of how similar (and differently) currently available amino acids descriptor sets behave when converting structure to property space. The results obtained enable molecular modelers to select suitable amino acid descriptor sets for structure-activity analyses, e.g. those showing complementary behavior. Show less
Infection with HIV cannot currently be cured; however it can be controlled by combination treatment with multiple anti-retroviral drugs. Given different viral genotypes for virtually each... Show moreInfection with HIV cannot currently be cured; however it can be controlled by combination treatment with multiple anti-retroviral drugs. Given different viral genotypes for virtually each individual patient, the question now arises which drug combination to use to achieve effective treatment. With the availability of viral genotypic data and clinical phenotypic data, it has become possible to create computational models able to predict an optimal treatment regimen for an individual patient. Current models are based only on sequence data derived from viral genotyping; chemical similarity of drugs is not considered. To explore the added value of chemical similarity inclusion we applied proteochemometric models, combining chemical and protein target properties in a single bioactivity model. Our dataset was a large scale clinical database of genotypic and phenotypic information (in total ca. 300,000 drug-mutant bioactivity data points, 4 (NNRTI), 8 (NRTI) or 9 (PI) drugs, and 10,700 (NNRTI) 10,500 (NRTI) or 27,000 (PI) mutants). Our models achieved a prediction error below 0.5 Log Fold Change. Moreover, when directly compared with previously published sequence data, derived models PCM performed better in resistance classification and prediction of Log Fold Change (0.76 log units versus 0.91). Furthermore, we were able to successfully confirm both known and identify previously unpublished, resistance-conferring mutations of HIV Reverse Transcriptase (e.g. K102Y, T216M) and HIV Protease (e.g. Q18N, N88G) from our dataset. Finally, we applied our models prospectively to the public HIV resistance database from Stanford University obtaining a correct resistance prediction rate of 84% on the full set (compared to 80% in previous work on a high quality subset). We conclude that proteochemometric models are able to accurately predict the phenotypic resistance based on genotypic data even for novel mutants and mixtures. Furthermore, we add an applicability domain to the prediction, informing the user about the reliability of predictions. Show less
Over the last decades several disciplines relevant to medicinal chemistry and preclinical drug discovery have made gigantic leaps; this includes chemistry, biology and measurement of bioactivity.... Show moreOver the last decades several disciplines relevant to medicinal chemistry and preclinical drug discovery have made gigantic leaps; this includes chemistry, biology and measurement of bioactivity. Better techniques have led to massive amounts of data. Moreover, sources of chemical and bioactivity data have become available in the public domain. Hence there is a need for new techniques combining and mining these data sources. This thesis focuses on computational methods combining data from these disciplines and demonstrates that the sum of these methods leads to better quality predictions than models using the individual data sources. One of the techniques central in this thesis is proteochemometric modeling, a machine learning approach linking chemical descriptors and protein descriptors to a biologically relevant output variable. This output variable describes the activity of molecules on biological macromolecules and hence proteochemometric models can make relevant predictions for both unseen molecules and unseen macromolecules (e.g. novel viral mutants). Secondly we present a novel technique that is able to combine information from multiple crystal structures in such a way that shared and unique pharmacophoric features can be isolated and visualized. Approaches presented here have been validated prospectively and have been shown to be widely applicable. Show less
