G Protein-coupled receptors (GPCRs) are the most frequently exploited drug target family, moreover they are often found mutated in cancer. Here we used a dataset of mutations found in patient... Show moreG Protein-coupled receptors (GPCRs) are the most frequently exploited drug target family, moreover they are often found mutated in cancer. Here we used a dataset of mutations found in patient samples derived from the Genomic Data Commons and compared it to the natural human variance as exemplified by data from the 1000 genomes project. We explored cancer-related mutation patterns in all GPCR classes combined and individually. While the location of the mutations across the protein domains did not differ significantly in the two datasets, a mutation enrichment in cancer patients was observed among class-specific conserved motifs in GPCRs such as the Class A "DRY" motif. A Two-Entropy Analysis confirmed the correlation between residue conservation and cancer-related mutation frequency. We subsequently created a ranking of high scoring GPCRs, using a multi-objective approach (Pareto Front Ranking). Our approach was confirmed by re-discovery of established cancer targets such as the LPA and mGlu receptor families, but also discovered novel GPCRs which had not been linked to cancer before such as the P2Y Receptor 10 (P2RY10). Overall, this study presents a list of GPCRs that are amenable to experimental follow up to elucidate their role in cancer. Show less
Hypoxia and HIF signaling drive cancer progression and therapy resistance and have been demonstrated in breast cancer. To what extent breast cancer subtypes differ in their response to hypoxia has... Show moreHypoxia and HIF signaling drive cancer progression and therapy resistance and have been demonstrated in breast cancer. To what extent breast cancer subtypes differ in their response to hypoxia has not been resolved. Here, we show that hypoxia similarly triggers HIF1 stabilization in luminal and basal A triple negative breast cancer cells and we use high throughput targeted RNA sequencing to analyze its effects on gene expression in these subtypes. We focus on regulation of YAP/TAZ/TEAD targets and find overlapping as well as distinct target genes being modulated in luminal and basal A cells under hypoxia. We reveal a HIF1 mediated, basal A specific response to hypoxia by which TAZ, but not YAP, is phosphorylated at Ser89. While total YAP/TAZ localization is not affected by hypoxia, hypoxia drives a shift of [p-TAZ(Ser89)/p-YAP(Ser127)] from the nucleus to the cytoplasm in basal A but not luminal breast cancer cells. Cell fractionation and YAP knock-out experiments confirm cytoplasmic sequestration of TAZ(Ser89) in hypoxic basal A cells. Pharmacological and genetic interference experiments identify c-Src and CDK3 as kinases involved in such phosphorylation of TAZ at Ser89 in hypoxic basal A cells. Hypoxia attenuates growth of basal A cells and the effect of verteporfin, a disruptor of YAP/TAZ-TEAD-mediated transcription, is diminished under those conditions, while expression of a TAZ-S89A mutant does not confer basal A cells with a growth advantage under hypoxic conditions, indicating that other hypoxia regulated pathways suppressing cell growth are dominant. Show less
Zwol, W. van; Rimbert, A.; Wolters, J.C.; Smit, M.; Bloks, V.W.; Kloosterhuis, N.J.; ... ; Kuivenhoven, J.A. 2022
The assembly and secretion of VLDL from the liver, a pathway that affects hepatic and plasma lipids, remains incompletely understood. We set out to identify players in the VLDL biogenesis pathway... Show moreThe assembly and secretion of VLDL from the liver, a pathway that affects hepatic and plasma lipids, remains incompletely understood. We set out to identify players in the VLDL biogenesis pathway by identifying genes that are co-expressed with the MTTP gene that encodes for microsomal triglyceride transfer protein, key to the lipidation of apolipoprotein B, the core protein of VLDL. Using human and murine transcriptomic data sets, we identified small leucine-rich protein 1 (SMLR1), encoding for small leucine-rich protein 1, a protein of unknown function that is exclusively expressed in liver and small intestine. To assess the role of SMLR1 in the liver, we used somatic CRISPR/CRISPR-associated protein 9 gene editing to silence murine Smlr1 in hepatocytes (Smlr1-LKO). When fed a chow diet, male and female mice show hepatic steatosis, reduced plasma apolipoprotein B and triglycerides, and reduced VLDL secretion without affecting microsomal triglyceride transfer protein activity. Immunofluorescence studies show that SMLR1 is in the endoplasmic reticulum and Cis-Golgi complex. The loss of hepatic SMLR1 in female mice protects against diet-induced hyperlipidemia and atherosclerosis but causes NASH. On a high-fat, high-cholesterol diet, insulin and glucose tolerance tests did not reveal differences in male Smlr1-LKO mice versus controls. We propose a role for SMLR1 in the trafficking of VLDL from the endoplasmic reticulum to the Cis-Golgi complex. While this study uncovers SMLR1 as a player in the VLDL assembly, trafficking, and secretion pathway, it also shows that NASH can occur with undisturbed glucose homeostasis and atheroprotection. Show less
Aims During phase I study conduct, blinded data are reviewed to predict the safety of increasing the dose level. The aim of the present study was to describe the probability that effects are... Show moreAims During phase I study conduct, blinded data are reviewed to predict the safety of increasing the dose level. The aim of the present study was to describe the probability that effects are observed in blinded evaluations of data in a simulated phase I study design. Methods An application was created to simulate blinded pharmacological response curves over time for 6 common safety/efficacy measurements in phase I studies for 1 or 2 cohorts (6 active, 2 placebo per cohort). Effect sizes between 0 and 3 between-measurement standard deviations (SDs) were simulated. Each set of simulated graphs contained the individual response and mean +/- SD over time. Reviewers (n = 34) reviewed a median of 100 simulated datasets and indicated whether an effect was present. Results Increasing effect sizes resulted in a higher chance of the effect being identified by the blinded reviewer. On average, 6% of effect sizes of 0.5 between-measurement SD were correctly identified, increasing to 72% in 3.0 between-measurement SD effect sizes. In contrast, on average 92-95% of simulations with no effect were correctly identified, with little effect of between-measurement variability in single cohort simulations. Adding a dataset of a second cohort at half the simulated dose did not appear to improve the interpretation. Conclusion Our analysis showed that effect sizes <2x the between-measurement SD of the investigated outcome frequently go unnoticed by blinded reviewers, indicating that the weight given to these blinded analyses in current phase I practice is inappropriate and should be re-evaluated. Show less
The adenosine A(2A) receptor (A(2A)AR) is a class A G-protein-coupled receptor (GPCR). It is an immune checkpoint in the tumor micro-environment and has become an emerging target for cancer... Show moreThe adenosine A(2A) receptor (A(2A)AR) is a class A G-protein-coupled receptor (GPCR). It is an immune checkpoint in the tumor micro-environment and has become an emerging target for cancer treatment. In this study, we aimed to explore the effects of cancer-patient-derived A(2A)AR mutations on ligand binding and receptor functions. The wild-type A(2A)AR and 15 mutants identified by Genomic Data Commons (GDC) in human cancers were expressed in HEK293T cells. Firstly, we found that the binding affinity for agonist NECA was decreased in six mutants but increased for the V275A mutant. Mutations A165V and A265V decreased the binding affinity for antagonist ZM241385. Secondly, we found that the potency of NECA (EC50) in an impedance-based cell-morphology assay was mostly correlated with the binding affinity for the different mutants. Moreover, S132L and H278N were found to shift the A(2A)AR towards the inactive state. Importantly, we found that ZM241385 could not inhibit the activation of V275A and P285L stimulated by NECA. Taken together, the cancer-associated mutations of A(2A)AR modulated ligand binding and receptor functions. This study provides fundamental insights into the structure-activity relationship of the A(2A)AR and provides insights for A(2A)AR-related personalized treatment in cancer. Show less
Gorostiola Gonzalez, M.; Janssen, A.P.A.; IJzerman, A.P.; Heitman, L.H.; Westen, G.J.P. van 2022
The integration of machine learning and structure-based methods has proven valuable in the past as a way to prioritize targets and compounds in early drug discovery. In oncological research, these... Show moreThe integration of machine learning and structure-based methods has proven valuable in the past as a way to prioritize targets and compounds in early drug discovery. In oncological research, these methods can be highly beneficial in addressing the diversity of neoplastic diseases portrayed by the different hallmarks of cancer. Here, we review six use case scenarios for integrated computational methods, namely driver prediction, computational mutagenesis, (off)-target prediction, binding site prediction, virtual screening, and allosteric modulation analysis. We address the heterogeneity of integration approaches and individual methods, while acknowledging their current limitations and highlighting their potential to bring drugs for personalized oncological therapies to the market faster. Show less
Overexpression of the adenosine A1 receptor (A1AR) has been detected in various cancer cell lines. However, the role of A1AR in tumor development is still unclear. Thirteen A1AR mutations were... Show moreOverexpression of the adenosine A1 receptor (A1AR) has been detected in various cancer cell lines. However, the role of A1AR in tumor development is still unclear. Thirteen A1AR mutations were identified in the Cancer Genome Atlas from cancer patient samples. We have investigated the pharmacology of the mutations located at the 7-transmembrane domain using a yeast system. Concentration-growth curves were obtained with the full agonist CPA and compared to the wild type hA1AR. H78L3.23 and S246T6.47 showed increased constitutive activity, while only the constitutive activity of S246T6.47 could be reduced to wild type levels by the inverse agonist DPCPX. Decreased constitutive activity was observed on five mutant receptors, among which A52V2.47 and W188C5.46 showed a diminished potency for CPA. Lastly, a complete loss of activation was observed in five mutant receptors. A selection of mutations was also investigated in a mammalian system, showing comparable effects on receptor activation as in the yeast system, except for residues pointing toward the membrane. Taken together, this study will enrich the view of the receptor structure and function of A1AR, enlightening the consequences of these mutations in cancer. Ultimately, this may provide an opportunity for precision medicine for cancer patients with pathological phenotypes involving these mutations. Show less
ObjectiveThe aim of the present study was to develop a neural network to characterize the effect of aging on the ECG in healthy volunteers. Moreover, the impact of the various ECG features on aging... Show moreObjectiveThe aim of the present study was to develop a neural network to characterize the effect of aging on the ECG in healthy volunteers. Moreover, the impact of the various ECG features on aging was evaluated.Methods & resultsA total of 6228 healthy subjects without structural heart disease were included in this study. A neural network regression model was created to predict age of the subjects based on their ECG; 577 parameters derived from a 12‑lead ECG of each subject were used to develop and validate the neural network; A tenfold cross-validation was performed, using 118 subjects for validation each fold. Using SHapley Additive exPlanations values the impact of the individual features on the prediction of age was determined. Of 6228 subjects tested, 1808 (29%) were females and mean age was 34 years, range 18-75 years. Physiologic age was estimated as a continuous variable with an average error of 6.9 ± 5.6 years (R2 = 0.72 ± 0.04). The correlation was slightly stronger for men (R2 = 0.74) than for women (R2 = 0.66). The most important features on the prediction of physiologic age were T wave morphology indices in leads V4 and V5, and P wave amplitude in leads AVR and II.ConclusionThe application of machine learning to the ECG using a neural network regression model, allows accurate estimation of physiologic cardiac age. This technique could be used to pick up subtle age-related cardiac changes, but also estimate the reversing of these age-associated effects by administered treatments. Show less
G protein-coupled receptors (GPCRs) are known to be involved in tumor progression and metastasis. The adenosine A1 receptor (A1 AR) has been detected to be over-expressed in various cancer cell... Show moreG protein-coupled receptors (GPCRs) are known to be involved in tumor progression and metastasis. The adenosine A1 receptor (A1 AR) has been detected to be over-expressed in various cancer cell lines. However, the role of A1 AR in tumor development is not yet well characterized. A series of A1 AR mutations were identified in the Cancer Genome Atlas from cancer patient samples. In this study, we have investigated the pharmacology of mutations located outside of the 7-transmembrane domain by using a "single-GPCR-one-G protein" yeast system. Concentration-growth curves were obtained with the full agonist CPA for 12 mutant receptors and compared to the wild-type hA1 AR. Most mutations located at the extracellular loops (EL) reduced the levels of constitutive activity of the receptor and agonist potency. For mutants at the intracellular loops (ILs) of the receptor, an increased constitutive activity was found for mutant receptor L211R5.69 , while a decreased constitutive activity and agonist response were found for mutant receptor L113F34.51 . Lastly, mutations identified on the C-terminus did not significantly influence the pharmacological function of the receptor. A selection of mutations was also investigated in a mammalian system. Overall, similar effects on receptor activation compared to the yeast system were found with mutations located at the EL, but some contradictory effects were observed for mutations located at the IL. Taken together, this study will enrich the insight of A1 AR structure and function, enlightening the consequences of these mutations in cancer. Ultimately, this may provide potential precision medicine in cancer treatment. Show less
Wall, H.E.C. van der; Hassing, G.J.; Doll, R.J.; Westen, G.J.P. van; Cohen, A.F.; Selder, J.L.; ... ; Gal, P. 2022
ObjectiveThe aim of the present study was to develop a neural network to characterize the effect of aging on the ECG in healthy volunteers. Moreover, the impact of the various ECG features on aging... Show moreObjectiveThe aim of the present study was to develop a neural network to characterize the effect of aging on the ECG in healthy volunteers. Moreover, the impact of the various ECG features on aging was evaluated.Methods & resultsA total of 6228 healthy subjects without structural heart disease were included in this study. A neural network regression model was created to predict age of the subjects based on their ECG; 577 parameters derived from a 12‑lead ECG of each subject were used to develop and validate the neural network; A tenfold cross-validation was performed, using 118 subjects for validation each fold. Using SHapley Additive exPlanations values the impact of the individual features on the prediction of age was determined. Of 6228 subjects tested, 1808 (29%) were females and mean age was 34 years, range 18–75 years. Physiologic age was estimated as a continuous variable with an average error of 6.9 ± 5.6 years (R2 = 0.72 ± 0.04) . The correlation was slightly stronger for men (R2 = 0.74) than for women (R2 = 0.66). The most important features on the prediction of physiologic age were T wave morphology indices in leads V4 and V5, and P wave amplitude in leads AVR and II.ConclusionThe application of machine learning to the ECG using a neural network regression model, allows accurate estimation of physiologic cardiac age. This technique could be used to pick up subtle age-related cardiac changes, but also estimate the reversing of these age-associated effects by administered treatments. Show less
Integration of statistical learning methods with structure-based modeling approaches is a contemporary strategy to identify novel lead compounds in drug discovery. Hepatic organic anion... Show moreIntegration of statistical learning methods with structure-based modeling approaches is a contemporary strategy to identify novel lead compounds in drug discovery. Hepatic organic anion transporting polypeptides (OATP1B1, OATP1B3, and OATP2B1) are classical off-targets, and it is well recognized that their ability to interfere with a wide range of chemically unrelated drugs, environmental chemicals, or food additives can lead to unwanted adverse effects like liver toxicity and drug-drug or drug-food interactions. Therefore, the identification of novel (tool) compounds for hepatic OATPs by virtual screening approaches and subsequent experimental validation is a major asset for elucidating structure-function relationships of (related) transporters: they enhance our understanding about molecular determinants and structural aspects of hepatic OATPs driving ligand binding and selectivity. In the present study, we performed a consensus virtual screening approach by using different types of machine learning models (proteochemometric models, conformal prediction models, and XGBoost models for hepatic OATPs), followed by molecular docking of preselected hits using previously established structural models for hepatic OATPs. Screening the diverse REAL drug-like set (Enamine) shows a comparable hit rate for OATP1B1 (36% actives) and OATP1B3 (32% actives), while the hit rate for OATP2B1 was even higher (66% actives). Percentage inhibition values for 44 selected compounds were determined using dedicated in vitro assays and guided the prioritization of several highly potent novel hepatic OATP inhibitors: six (strong) OATP2B1 inhibitors (IC50 values ranging from 0.04 to 6 μM), three OATP1B1 inhibitors (2.69 to 10 μM), and five OATP1B3 inhibitors (1.53 to 10 μM) were identified. Strikingly, two novel OATP2B1 inhibitors were uncovered (C7 and H5) which show high affinity (IC50 values: 40 nM and 390 nM) comparable to the recently described estrone-based inhibitor (IC50 = 41 nM). A molecularly detailed explanation for the observed differences in ligand binding to the three transporters is given by means of structural comparison of the detected binding sites and docking poses. Show less
Toxicology has been an active research field for many decades, with academic, industrial and government involvement. Modern omics and computational approaches are changing the field, from merely... Show moreToxicology has been an active research field for many decades, with academic, industrial and government involvement. Modern omics and computational approaches are changing the field, from merely disease-specific observational models into target-specific predictive models. Traditionally, toxicology has strong links with other fields such as biology, chemistry, pharmacology and medicine. With the rise of synthetic and new engineered materials, alongside ongoing prioritisation needs in chemical risk assessment for existing chemicals, early predictive evaluations are becoming of utmost importance to both scientific and regulatory purposes. ELIXIR is an intergovernmental organisation that brings together life science resources from across Europe. To coordinate the linkage of various life science efforts around modern predictive toxicology, the establishment of a new ELIXIR Community is seen as instrumental. In the past few years, joint efforts, building on incidental overlap, have been piloted in the context of ELIXIR. For example, the EU-ToxRisk, diXa, HeCaToS, transQST, and the nanotoxicology community have worked with the ELIXIR TeSS, Bioschemas, and Compute Platforms and activities. In 2018, a core group of interested parties wrote a proposal, outlining a sketch of what this new ELIXIR Toxicology Community would look like. A recent workshop (held September 30th to October 1st, 2020) extended this into an ELIXIR Toxicology roadmap and a shortlist of limited investment-high gain collaborations to give body to this new community. This Whitepaper outlines the results of these efforts and defines our vision of the ELIXIR Toxicology Community and how it complements other ELIXIR activities. Show less
Transmembranal G Protein-Coupled Receptors (GPCRs) transduce extracellular chemical signals to the cell, via conformational change from a resting (inactive) to an active (canonically bound to a G... Show moreTransmembranal G Protein-Coupled Receptors (GPCRs) transduce extracellular chemical signals to the cell, via conformational change from a resting (inactive) to an active (canonically bound to a G-protein) conformation. Receptor activation is normally modulated by extracellular ligand binding, but mutations in the receptor can also shift this equilibrium by stabilizing different conformational states. In this work, we built structure-energetic relationships of receptor activation based on original thermodynamic cycles that represent the conformational equilibrium of the prototypical A2A adenosine receptor (AR). These cycles were solved with efficient free energy perturbation (FEP) protocols, allowing to distinguish the pharmacological profile of different series of A2AAR agonists with different efficacies. The modulatory effects of point mutations on the basal activity of the receptor or on ligand efficacies could also be detected. This methodology can guide GPCR ligand design with tailored pharmacological properties, or allow the identification of mutations that modulate receptor activation with potential clinical implications. Show less
Liu, X.; Ye, K.; Vlijmen, H.W.T. van; Emmerich, M.T.M.; IJzerman, A.P.; Westen, G.J.P. van 2021
In polypharmacology drugs are required to bind to multiple specific targets, for example to enhance efficacy or to reduce resistance formation. Although deep learning has achieved a breakthrough in... Show moreIn polypharmacology drugs are required to bind to multiple specific targets, for example to enhance efficacy or to reduce resistance formation. Although deep learning has achieved a breakthrough in de novo design in drug discovery, most of its applications only focus on a single drug target to generate drug-like active molecules. However, in reality drug molecules often interact with more than one target which can have desired (polypharmacology) or undesired (toxicity) effects. In a previous study we proposed a new method named DrugEx that integrates an exploration strategy into RNN-based reinforcement learning to improve the diversity of the generated molecules. Here, we extended our DrugEx algorithm with multi-objective optimization to generate drug-like molecules towards multiple targets or one specific target while avoiding off-targets (the two adenosine receptors, A1AR and A2AAR, and the potassium ion channel hERG in this study). In our model, we applied an RNN as the agent and machine learning predictors as the environment. Both the agent and the environment were pre-trained in advance and then interplayed under a reinforcement learning framework. The concept of evolutionary algorithms was merged into our method such that crossover and mutation operations were implemented by the same deep learning model as the agent. During the training loop, the agent generates a batch of SMILES-based molecules. Subsequently scores for all objectives provided by the environment are used to construct Pareto ranks of the generated molecules. For this ranking a non-dominated sorting algorithm and a Tanimoto-based crowding distance algorithm using chemical fingerprints are applied. Here, we adopted GPU acceleration to speed up the process of Pareto optimization. The final reward of each molecule is calculated based on the Pareto ranking with the ranking selection algorithm. The agent is trained under the guidance of the reward to make sure it can generate desired molecules after convergence of the training process. All in all we demonstrate generation of compounds with a diverse predicted selectivity profile towards multiple targets, offering the potential of high efficacy and low toxicity. Show less
Many contemporary cheminformatics methods, including computer-aided de novo drug design, hold promise to significantly accelerate and reduce the cost of drug discovery. Thanks to this attractive... Show moreMany contemporary cheminformatics methods, including computer-aided de novo drug design, hold promise to significantly accelerate and reduce the cost of drug discovery. Thanks to this attractive outlook, the field has thrived and in the past few years has seen an especially significant growth, mainly due to the emergence of novel methods based on deep neural networks. This growth is also apparent in the development of novel de novo drug design methods with many new generative algorithms now available. However, widespread adoption of new generative techniques in the fields like medicinal chemistry or chemical biology is still lagging behind the most recent developments. Upon taking a closer look, this fact is not surprising since in order to successfully integrate the most recent de novo drug design methods in existing processes and pipelines, a close collaboration between diverse groups of experimental and theoretical scientists needs to be established. Therefore, to accelerate the adoption of both modern and traditional de novo molecular generators, we developed Generator User Interface (GenUI), a software platform that makes it possible to integrate molecular generators within a feature-rich graphical user interface that is easy to use by experts of diverse backgrounds. GenUI is implemented as a web service and its interfaces offer access to cheminformatics tools for data preprocessing, model building, molecule generation, and interactive chemical space visualization. Moreover, the platform is easy to extend with customizable frontend React.js components and backend Python extensions. GenUI is open source and a recently developed de novo molecular generator, DrugEx, was integrated as a proof of principle. In this work, we present the architecture and implementation details of GenUI and discuss how it can facilitate collaboration in the disparate communities interested in de novo molecular generation and computer-aided drug discovery. Show less
