Exposure to PFAS chemicals induces sex-dependent alterations in key rate-limiting steps of lipid metabolism in liver steatosis

Toxicants with the potential to bioaccumulate in humans and animals have long been a cause for concern, particularly due to their association with multiple diseases and organ injuries. Per- and polyfluoro alkyl substances (PFAS) and polycyclic aromatic hydrocarbons (PAH) are two such classes of chemicals that bioaccumulate and have been associated with steatosis in the liver. Although PFAS and PAH are classified as chemicals of concern, their molecular mechanisms of toxicity remain to be explored in detail. In this study, we aimed to identify potential mechanisms by which an acute exposure to PFAS and PAH chemicals can induce lipid accumulation and whether the responses depend on chemical class, dose, and sex. To this end, we analyzed mechanisms beginning with the binding of the chemical to a molecular initiating event (MIE) and the consequent transcriptomic alterations. We collated potential MIEs using predictions from our previously developed ToxProfiler tool and from published steatosis adverse outcome pathways. Most of the MIEs are transcription factors, and we collected their target genes by mining the TRRUST database. To analyze the effects of PFAS and PAH on the steatosis mechanisms, we performed a computational MIE-target gene analysis on high-throughput transcriptomic measurements of liver tissue from male and female rats exposed to either a PFAS or PAH. The results showed peroxisome proliferator-activated receptor (PPAR)-α targets to be the most dysregulated, with most of the genes being upregulated. Furthermore, PFAS exposure disrupted several lipid metabolism genes, including upregulation of fatty acid oxidation genes (Acadm, Acox1, Cpt2, Cyp4a1-3) and downregulation of lipid transport genes (Apoa1, Apoa5, Pltp). We also identified multiple genes with sex-specific behavior. Notably, the rate-limiting genes of gluconeogenesis (Pck1) and bile acid synthesis (Cyp7a1) were specifically downregulated in male rats compared to female rats, while the rate-limiting gene of lipid synthesis (Scd) showed a PFAS-specific upregulation. The results suggest that the PPAR signaling pathway plays a major role in PFAS-induced lipid accumulation in rats. Together, these results show that PFAS exposure induces a sex-specific multi-factorial mechanism involving rate-limiting genes of gluconeogenesis and bile acid synthesis that could lead to activation of an adverse outcome pathway for steatosis

A physicochemical descriptor-based scoring scheme for effective and rapid filtering of kinase-like chemical space

<p>Abstract</p> <p>Background</p> <p>The current chemical space of known small molecules is estimated to exceed 10⁶⁰structures. Though the largest physical compound repositories contain only a few tens of millions of unique compounds, virtual screening of databases of this size is still difficult. In recent years, the application of physicochemical descriptor-based profiling, such as Lipinski's rule-of-five for drug-likeness and Oprea's criteria of lead-likeness, as early stage filters in drug discovery has gained widespread acceptance. In the current study, we outline a kinase-likeness scoring function based on known kinase inhibitors.</p> <p>Results</p> <p>The method employs a collection of 22,615 known kinase inhibitors from the ChEMBL database. A kinase-likeness score is computed using statistical analysis of nine key physicochemical descriptors for these inhibitors. Based on this score, the kinase-likeness of four publicly and commercially available databases, i.e., National Cancer Institute database (NCI), the Natural Products database (NPD), the National Institute of Health's Molecular Libraries Small Molecule Repository (MLSMR), and the World Drug Index (WDI) database, is analyzed. Three of these databases, i.e., NCI, NPD, and MLSMR are frequently used in the virtual screening of kinase inhibitors, while the fourth WDI database is for comparison since it covers a wide range of known chemical space. Based on the kinase-likeness score, a kinase-focused library is also developed and tested against three different kinase targets selected from three different branches of the human kinome tree.</p> <p>Conclusions</p> <p>Our proposed methodology is one of the first that explores how the narrow chemical space of kinase inhibitors and its relevant physicochemical information can be utilized to build kinase-focused libraries and prioritize pre-existing compound databases for screening. We have shown that focused libraries generated by filtering compounds using the kinase-likeness score have, on average, better docking scores than an equivalent number of randomly selected compounds. Beyond library design, our findings also impact the broader efforts to identify kinase inhibitors by screening pre-existing compound libraries. Currently, the NCI library is the most commonly used database for screening kinase inhibitors. Our research suggests that other libraries, such as MLSMR, are more kinase-like and should be given priority in kinase screenings.</p

Molecular Structure-Based Large-Scale Prediction of Chemical-Induced Gene Expression Changes

The quantitative structure–activity relationship (QSAR) approach has been used to model a wide range of chemical-induced biological responses. However, it had not been utilized to model chemical-induced genomewide gene expression changes until very recently, owing to the complexity of training and evaluating a very large number of models. To address this issue, we examined the performance of a variable nearest neighbor (<i>v</i>-NN) method that uses information on near neighbors conforming to the principle that similar structures have similar activities. Using a data set of gene expression signatures of 13 150 compounds derived from cell-based measurements in the NIH Library of Integrated Network-based Cellular Signatures program, we were able to make predictions for 62% of the compounds in a 10-fold cross validation test, with a correlation coefficient of 0.61 between the predicted and experimentally derived signaturesa reproducibility rivaling that of high-throughput gene expression measurements. To evaluate the utility of the predicted gene expression signatures, we compared the predicted and experimentally derived signatures in their ability to identify drugs known to cause specific liver, kidney, and heart injuries. Overall, the predicted and experimentally derived signatures had similar receiver operating characteristics, whose areas under the curve ranged from 0.71 to 0.77 and 0.70 to 0.73, respectively, across the three organ injury models. However, detailed analyses of enrichment curves indicate that signatures predicted from multiple near neighbors outperformed those derived from experiments, suggesting that averaging information from near neighbors may help improve the signal from gene expression measurements. Our results demonstrate that the <i>v</i>-NN method can serve as a practical approach for modeling large-scale, genomewide, chemical-induced, gene expression changes

Predicting Rat and Human Pregnane X Receptor Activators Using Bayesian Classification Models

The pregnane X receptor (PXR) is a ligand-activated transcription factor that acts as a master regulator of metabolizing enzymes and transporters. To avoid adverse drug–drug interactions and diseases such as steatosis and cancers associated with PXR activation, identifying drugs and chemicals that activate PXR is of crucial importance. In this work, we developed ligand-based predictive computational models for both rat and human PXR activation, which allowed us to identify potentially harmful chemicals and evaluate species-specific effects of a given compound. We utilized a large publicly available data set of nearly 2000 compounds screened in cell-based reporter gene assays to develop Bayesian quantitative structure–activity relationship models using physicochemical properties and structural descriptors. Our analysis showed that PXR activators tend to be hydrophobic and significantly different from nonactivators in terms of their physicochemical properties such as molecular weight, logP, number of rings, and solubility. Our Bayesian models, evaluated by using 5-fold cross-validation, displayed a sensitivity of 75% (76%), specificity of 76% (75%), and accuracy of 89% (89%) for human (rat) PXR activation. We identified structural features shared by rat and human PXR activators as well as those unique to each species. We compared rat <i>in vitro</i> PXR activation data to <i>in vivo</i> data by using DrugMatrix, a large toxicogenomics database with gene expression data obtained from rats after exposure to diverse chemicals. Although <i>in vivo</i> gene expression data pointed to cross-talk between nuclear receptor activators that is captured only by <i>in vivo</i> assays, overall we found broad agreement between <i>in vitro</i> and <i>in vivo</i> PXR activation. Thus, the models developed here serve primarily as efficient initial high-throughput <i>in silico</i> screens of <i>in vitro</i> activity

Predicting Rat and Human Pregnane X Receptor Activators Using Bayesian Classification Models

The pregnane X receptor (PXR) is a ligand-activated transcription factor that acts as a master regulator of metabolizing enzymes and transporters. To avoid adverse drug–drug interactions and diseases such as steatosis and cancers associated with PXR activation, identifying drugs and chemicals that activate PXR is of crucial importance. In this work, we developed ligand-based predictive computational models for both rat and human PXR activation, which allowed us to identify potentially harmful chemicals and evaluate species-specific effects of a given compound. We utilized a large publicly available data set of nearly 2000 compounds screened in cell-based reporter gene assays to develop Bayesian quantitative structure–activity relationship models using physicochemical properties and structural descriptors. Our analysis showed that PXR activators tend to be hydrophobic and significantly different from nonactivators in terms of their physicochemical properties such as molecular weight, logP, number of rings, and solubility. Our Bayesian models, evaluated by using 5-fold cross-validation, displayed a sensitivity of 75% (76%), specificity of 76% (75%), and accuracy of 89% (89%) for human (rat) PXR activation. We identified structural features shared by rat and human PXR activators as well as those unique to each species. We compared rat <i>in vitro</i> PXR activation data to <i>in vivo</i> data by using DrugMatrix, a large toxicogenomics database with gene expression data obtained from rats after exposure to diverse chemicals. Although <i>in vivo</i> gene expression data pointed to cross-talk between nuclear receptor activators that is captured only by <i>in vivo</i> assays, overall we found broad agreement between <i>in vitro</i> and <i>in vivo</i> PXR activation. Thus, the models developed here serve primarily as efficient initial high-throughput <i>in silico</i> screens of <i>in vitro</i> activity

Probing Liver Injuries Induced by Thioacetamide in Human In Vitro Pooled Hepatocyte Experiments

Animal studies are typically utilized to understand the complex mechanisms associated with toxicant-induced hepatotoxicity. Among the alternative approaches to animal studies, in vitro pooled human hepatocytes have the potential to capture population variability. Here, we examined the effect of the hepatotoxicant thioacetamide on pooled human hepatocytes, divided into five lots, obtained from forty diverse donors. For 24 h, pooled human hepatocytes were exposed to vehicle, 1.33 mM (low dose), and 12 mM (high dose) thioacetamide, followed by RNA-seq analysis. We assessed gene expression variability using heat maps, correlation plots, and statistical variance. We used KEGG pathways and co-expression modules to identify underlying physiological processes/pathways. The co-expression module analysis showed that the majority of the lots exhibited activation for the bile duct proliferation module. Despite lot-to-lot variability, we identified a set of common differentially expressed genes across the lots with similarities in their response to amino acid, lipid, and carbohydrate metabolism. We also examined efflux transporters and found larger lot-to-lot variability in their expression patterns, indicating a potential for alteration in toxicant bioavailability within the cells, which could in turn affect the gene expression patterns between the lots. Overall, our analysis highlights the challenges in using pooled hepatocytes to understand mechanisms of toxicity

Identifying a predictive gene signature and signaling networks/pathways associated with acute kidney injury

Poster presented in SOT 2016<div><br></div><div>Understanding the molecular mechanisms and signaling networks of acute kidney injury (AKI) will aid in biomarker development. In this study, we carried out co-expression-based analyses of DrugMatrix, a toxicogenomics database with kidney gene expression data from rats after exposure to diverse chemicals. We used the iterative signature algorithm and exhaustively generated modules using 50 different parameter combinations. We clustered the modules using gene and condition overlap scores and obtained 16 module clusters. Two of the module clusters showed activation in chemical exposures causing kidney injury and mapped well-known AKI marker genes such as <i>Havcr1</i>, <i>Tff3,</i> and <i>Clu</i>. We used the genes in these AKI-relevant module clusters to develop a signature of 30 genes that could assess the potential of a chemical to cause kidney injury well before injury actually occurs. We integrated AKI-relevant module cluster genes with protein-protein interaction networks and identified the involvement of immunoproteasomes in AKI. To identify biological networks and processes linked to <em>Havcr1</em>, we determined genes within the modules that frequently co-express with <em>Havcr1</em>, including <em>Cd44</em>, <em>Plk2</em>, <em>Mdm2</em>, <em>Hnmt</em>, <em>Macrod1</em>, and <em>Gtpbp4</em>. In this gene set, CD44 is a potential non-invasive biomarker candidate as it is up-regulated during AKI, undergoes cleavage of its ectodomain, and is secreted in urine. Overall, our analysis shows data mining of toxicological big data and identification of new insights/biomarker candidates for acute kidney injury.</div

Assessing Kidney Injury Induced by Mercuric Chloride in Guinea Pigs with In Vivo and In Vitro Experiments

Acute kidney injury, which is associated with high levels of morbidity and mortality, affects a significant number of individuals, and can be triggered by multiple factors, such as medications, exposure to toxic chemicals or other substances, disease, and trauma. Because the kidney is a critical organ, understanding and identifying early cellular or gene-level changes can provide a foundation for designing medical interventions. In our earlier work, we identified gene modules anchored to histopathology phenotypes associated with toxicant-induced liver and kidney injuries. Here, using in vivo and in vitro experiments, we assessed and validated these kidney injury-associated modules by analyzing gene expression data from the kidneys of male Hartley guinea pigs exposed to mercuric chloride. Using plasma creatinine levels and cell-viability assays as measures of the extent of renal dysfunction under in vivo and in vitro conditions, we performed an initial range-finding study to identify the appropriate doses and exposure times associated with mild and severe kidney injuries. We then monitored changes in kidney gene expression at the selected doses and time points post-toxicant exposure to characterize the mechanisms of kidney injury. Our injury module-based analysis revealed a dose-dependent activation of several phenotypic cellular processes associated with dilatation, necrosis, and fibrogenesis that were common across the experimental platforms and indicative of processes that initiate kidney damage. Furthermore, a comparison of activated injury modules between guinea pigs and rats indicated a strong correlation between the modules, highlighting their potential for cross-species translational studies