High throughput functional profiling of genes at intraocular pressure loci reveals distinct networks for glaucoma.

INTRODUCTION: Primary open angle glaucoma (POAG) is a leading cause of blindness globally. Characterized by progressive retinal ganglion cell degeneration, the precise pathogenesis remains unknown. Genome-wide association studies (GWAS) have uncovered many genetic variants associated with elevated intraocular pressure (IOP), one of the key risk factors for POAG. We aimed to identify genetic and morphological variation that can be attributed to trabecular meshwork cell (TMC) dysfunction and raised IOP in POAG. METHODS: 62 genes across 55 loci were knocked-out in a primary human TMC line. Each knockout group, including five non-targeting control groups, underwent single-cell RNA-sequencing (scRNA-seq) for differentially-expressed gene (DEG) analysis. Multiplexed fluorescence coupled with CellProfiler image analysis allowed for single-cell morphological profiling. RESULTS: Many gene knockouts invoked DEGs relating to matrix metalloproteinases and interferon-induced proteins. We have prioritized genes at four loci of interest to identify gene knockouts that may contribute to the pathogenesis of POAG, including ANGPTL2, LMX1B, CAV1, and KREMEN1. Three genetic networks of gene knockouts with similar transcriptomic profiles were identified, suggesting a synergistic function in trabecular meshwork cell physiology. TEK knockout caused significant upregulation of nuclear granularity on morphological analysis, while knockout of TRIOBP, TMCO1 and PLEKHA7 increased granularity and intensity of actin and the cell-membrane. CONCLUSION: High-throughput analysis of cellular structure and function through multiplex fluorescent single-cell analysis and scRNA-seq assays enabled the direct study of genetic perturbations at the single-cell resolution. This work provides a framework for investigating the role of genes in the pathogenesis of glaucoma and heterogenous diseases with a strong genetic basis.

Profiling Cellular Morphological Changes Induced by Dual-Targeting PROTACs of Aurora Kinase and RNA-Binding Protein YTHDF2.

Proteolysis targeting chimeras (PROTACs) are new chemical modalities that degrade proteins of interest, including established kinase targets and emerging RNA-binding proteins (RBPs). Whereas diverse sets of biochemical, biophysical and cellular assays are available for the evaluation and optimizations of PROTACs in understanding the involved ubiquitin-proteasome-mediated degradation mechanism and the structure-degradation relationship, a phenotypic method profiling the cellular morphological changes is rarely used. In this study, first, we reported the only examples of PROTACs degrading the mRNA-binding protein YTHDF2 via screening of multikinase PROTACs. Second, we reported the profiling of cellular morphological changes of the dual kinase- and RBP-targeting PROTACs using the unbiased cell painting assay (CPA). The CPA analysis revealed the high biosimilarity with the established aurora kinase cluster and annotated aurora kinase inhibitors, which reflected the association between YTHDF2 and the aurora kinase signaling network. Broadly, the results demonstrated that the cell painting assay can be a straightforward and powerful approach to evaluate PROTACs. Complementary to the existing biochemical, biophysical and cellular assays, CPA provided a new perspective in characterizing PROTACs at the cellular morphology.

Distinct profile of antiviral drugs effects in aortic and pulmonary endothelial cells revealed by high-content microscopy and cell painting assays.

Antiretroviral therapy have significantly improved the treatment of viral infections and reduced the associated mortality and morbidity rates. However, highly effective antiretroviral therapy (HAART) may lead to an increased risk of cardiovascular diseases, which could be related to endothelial toxicity. Here, seven antiviral drugs (remdesivir, PF-00835231, ritonavir, lopinavir, efavirenz, zidovudine and abacavir) were characterized against aortic (HAEC) and pulmonary (hLMVEC) endothelial cells, using high-content microscopy. The colourimetric study (MTS test) revealed similar toxicity profiles of all antiviral drugs tested in the concentration range of 1 nM-50 µM in aortic and pulmonary endothelial cells. Conversely, the drugs' effects on morphological parameters were more pronounced in HAECs as compared with hLMVECs. Based on the antiviral drugs' effects on the cytoplasmic and nuclei architecture (analyzed by multiple pre-defined parameters including SER texture and STAR morphology), the studied compounds were classified into five distinct morphological subgroups, each linked to a specific cellular response profile. In relation to morphological subgroup classification, antiviral drugs induced a loss of mitochondrial membrane potential, elevated ROS, changed lipid droplets/lysosomal content, decreased von Willebrand factor expression and micronuclei formation or dysregulated cellular autophagy. In conclusion, based on specific changes in endothelial cytoplasm, nuclei and subcellular morphology, the distinct endothelial response was identified for remdesivir, ritonavir, lopinavir, efavirenz, zidovudine and abacavir treatments. The effects detected in aortic endothelial cells were not detected in pulmonary endothelial cells. Taken together, high-content microscopy has proven to be a robust and informative method for endothelial drug profiling that may prove useful in predicting the organ-specific endothelial toxicity of various drugs.

Cell Painting unravels insecticidal modes of action on Spodoptera frugiperda insect cells.

The "Cell Painting" technology utilizes multiplexed fluorescent staining of various cell organelles, to produce high-content microscopy images of cells for multidimensional phenotype assessment. The phenotypic profiles extracted from those images can be analyzed upon perturbations with biologically active molecules to annotate the mode of action or biological activity by comparison with reference profiles of already known mechanisms of action, ultimately enabling the determination of on-target and off-target effects. This approach is already described in various human cell cultures, the most commonly used being the U2OS cell line, yet allows broad applications in additional areas of chemical-biological research. Here we describe for the first time the application and adaptation of Cell Painting to an insect cell line, the Sf9 cells from Spodoptera frugiperda. By adjusting image acquisition and analysis models, specific phenotypic profiles were obtained in a dose-dependent manner for 20 reference compounds, including representatives for the most relevant insecticidal modes of action categories (nerve & muscle, respiration and growth & development). Through a dimensionality-reduction method, both calculations of phenotypic half maximal inhibition concentration (IC50) values as well as similarity analysis of the obtained profiles by hierarchical clustering were performed. By Cell Painting effects on the phenotype could be obtained at higher sensitivity than in other assay formats, such as cytotoxicity assessments. More importantly, these analyses provide insight into mechanistic determinants of biological activity. Compounds with similar modes of action showed a high degree of proximity in a hierarchical clustering analysis while being distinct from actives with an unrelated mode of action. In essence, we provide strong evidence on the impact of Cell Painting mechanistic understanding of insecticides with regards to determinants of efficacy and safety utilizing an insect cell model system.

[Phenotypic drug discovery promotes research and development of innovative drugs based on traditional Chinese medicine].

Traditional Chinese medicine with rich resources in China and definite therapeutic effects on complex diseases demonstrates great development potential. However, the complex composition, the unclear pharmacodynamic substances and mechanisms of action, and the lack of reasonable methods for evaluating clinical safety and efficacy have limited the research and development of innovative drugs based on traditional Chinese medicine. The progress in cutting-edge disciplines such as artificial intelligence and biomimetics, especially the emergence of cell painting and organ-on-a-chip, helps to identify and characterize the active ingredients in traditional Chinese medicine based on the changes in model characteristics, thus providing more accurate guidance for the development and application of traditional Chinese medicine. The application of phenotypic drug discovery in the research and development of innovative drugs based on traditional Chinese medicine is gaining increasing attention. In recent years, the technology for phenotypic drug discovery keeps advancing, which improves the early discovery rate of new drugs and the success rate of drug research and development. Accordingly, phenotypic drug discovery gradually becomes a key tool for the research on new drugs. This paper discusses the enormous potential of traditional Chinese medicine in the discovery and development of innovative drugs and illustrates how the application of phenotypic drug discovery, supported by cutting-edge technologies such as cell painting, deep learning, and organ-on-a-chip, propels traditional Chinese medicine into a new stage of development.

Identification of Biologically Diverse Tetrahydronaphthalen-2-ols through the Synthesis and Phenotypic Profiling of Chemically Diverse, Estradiol-Inspired Compounds.

Combining natural product fragments to design new scaffolds with unprecedented bioactivity is a powerful strategy for the discovery of tool compounds and potential therapeutics. However, the choice of fragments to couple and the biological screens to employ remain open questions in the field. By choosing a primary fragment containing the A/B ring system of estradiol and fusing it to nine different secondary fragments, we were able to identify compounds that modulated four different phenotypes: inhibition of autophagy and osteoblast differentiation, as well as potassium channel and tubulin modulation. The latter two were uncovered by using unbiased morphological profiling with a cell-painting assay. The number of hits and variety in bioactivity discovered validates the use of recombining natural product fragments coupled to phenotypic screening for the rapid identification of biologically diverse compounds.

Assessing the performance of the Cell Painting assay across different imaging systems.

Quantitative microscopy is a powerful method for performing phenotypic screens from which image-based profiling can extract a wealth of information, termed profiles. These profiles can be used to elucidate the changes in cellular phenotypes across cell populations from different patient samples or following genetic or chemical perturbations. One such image-based profiling method is the Cell Painting assay, which provides morphological insight through the imaging of eight cellular compartments. Here, we examine the performance of the Cell Painting assay across multiple high-throughput microscope systems and find that all are compatible with this assay. Furthermore, we determine independently for each microscope system the best performing settings, providing those who wish to adopt this assay an ideal starting point for their own assays. We also explore the impact of microscopy setting changes in the Cell Painting assay and find that few dramatically reduce the quality of a Cell Painting profile, regardless of the microscope used.

Linking chemicals, genes and morphological perturbations to diseases.

The progress in image-based high-content screening technology has facilitated high-throughput phenotypic profiling notably the quantification of cell morphology perturbation by chemicals. However, understanding the mechanism of action of a chemical and linking it to cell morphology and phenotypes remains a challenge in drug discovery. In this study, we intended to integrate molecules that induced transcriptomic perturbations and cellular morphological changes into a biological network in order to assess chemical-phenotypic relationships in humans. Such a network was enriched with existing disease information to suggest molecular and cellular profiles leading to phenotypes. Two datasets were used for this study. Firstly, we used the "Cell Painting morphological profiling assay" dataset, composed of 30,000 compounds tested on human osteosarcoma cells (named U2OS). Secondly, we used the "L1000 mRNA profiling assay" dataset, a collection of transcriptional expression data from cultured human cells treated with approximately 20,000 bioactive small molecules from the Library of Integrated Network-based Cellular Signatures (LINCS). Furthermore, pathways, gene ontology terms and disease enrichments were performed on the transcriptomics data. Overall, our study makes it possible to develop a biological network combining chemical-gene-pathway-morphological perturbation and disease relationships. It contains an ensemble of 9989 chemicals, 732 significant morphological features and 12,328 genes. Through diverse examples, we demonstrated that some drugs shared similar genes, pathways and morphological profiles that, taken together, could help in deciphering chemical-phenotype observations.

Application of Cell Painting for chemical hazard evaluation in support of screening-level chemical assessments.

'Cell Painting' is an imaging-based high-throughput phenotypic profiling (HTPP) method in which cultured cells are fluorescently labeled to visualize subcellular structures (i.e., nucleus, nucleoli, endoplasmic reticulum, cytoskeleton, Golgi apparatus / plasma membrane and mitochondria) and to quantify morphological changes in response to chemicals or other perturbagens. HTPP is a high-throughput and cost-effective bioactivity screening method that detects effects associated with many different molecular mechanisms in an untargeted manner, enabling rapid in vitro hazard assessment for thousands of chemicals. Here, 1201 chemicals from the ToxCast library were screened in concentration-response up to ∼100 µM in human U-2 OS cells using HTPP. A phenotype altering concentration (PAC) was estimated for chemicals active in the tested range. PACs tended to be higher than lower bound potency values estimated from a broad collection of targeted high-throughput assays, but lower than the threshold for cytotoxicity. In vitro to in vivo extrapolation (IVIVE) was used to estimate administered equivalent doses (AEDs) based on PACs for comparison to human exposure predictions. AEDs for 18/412 chemicals overlapped with predicted human exposures. Phenotypic profile information was also leveraged to identify putative mechanisms of action and group chemicals. Of 58 known nuclear receptor modulators, only glucocorticoids and retinoids produced characteristic profiles; and both receptor types are expressed in U-2 OS cells. Thirteen chemicals with profile similarity to glucocorticoids were tested in a secondary screen and one chemical, pyrene, was confirmed by an orthogonal gene expression assay as a novel putative GR modulating chemical. Most active chemicals demonstrated profiles not associated with a known mechanism-of-action. However, many structurally related chemicals produced similar profiles, with exceptions such as diniconazole, whose profile differed from other active conazoles. Overall, the present study demonstrates how HTPP can be applied in screening-level chemical assessments through a series of examples and brief case studies.

Transcriptomics and cell painting analysis reveals molecular and morphological features associated with fed-batch production performance in CHO recombinant clones.

Stable, highly productive mammalian cells are critical for manufacturing affordable and effective biological medicines. Establishing a rational design of optimal biotherapeutic expression systems requires understanding how cells support the high demand for efficient biologics production. To that end, we performed transcriptomics and high-throughput imaging studies to identify putative genes and morphological features that underpin differences in antibody productivity among clones from a Chinese hamster ovary cell line. During log phase growth, we found that the expression of genes involved in biological processes related to cellular morphology varied significantly between clones with high specific productivity (qP > 35 pg/cell/day) and low specific productivity (qP < 20 pg/cell/day). At Day 10 of a fed-batch production run, near peak viable cell density, differences in gene expression related to metabolism, epigenetic regulation, and proliferation became prominent. Furthermore, we identified a subset of genes whose expression predicted overall productivity, including glutathione synthetase (Gss) and lactate dehydrogenase A (LDHA). Finally, we demonstrated the feasibility of cell painting coupled with high-throughput imaging to assess the morphological properties of intracellular organelles in relation to growth and productivity in fed-batch production. Our efforts lay the groundwork for systematic elucidation of clone performance using a multiomics approach that can guide future process design strategies.

Combining phenotypic profiling and targeted RNA-Seq reveals linkages between transcriptional perturbations and chemical effects on cell morphology: Retinoic acid as an example.

The United States Environmental Protection Agency has proposed a tiered testing strategy for chemical hazard evaluation based on new approach methods (NAMs). The first tier includes in vitro profiling assays applicable to many (human) cell types, such as high-throughput transcriptomics (HTTr) and high-throughput phenotypic profiling (HTPP). The goals of this study were to: (1) harmonize the seeding density of U-2 OS human osteosarcoma cells for use in both assays; (2) compare HTTr- versus HTPP-derived potency estimates for 11 mechanistically diverse chemicals; (3) identify candidate reference chemicals for monitoring assay performance in future screens; and (4) characterize the transcriptional and phenotypic changes in detail for all-trans retinoic acid (ATRA) as a model compound known for its adverse effects on osteoblast differentiation. The results of this evaluation showed that (1) HTPP conducted at low (400 cells/well) and high (3000 cells/well) seeding densities yielded comparable potency estimates and similar phenotypic profiles for the tested chemicals; (2) HTPP and HTTr resulted in comparable potency estimates for changes in cellular morphology and gene expression, respectively; (3) three test chemicals (etoposide, ATRA, dexamethasone) produced concentration-dependent effects on cellular morphology and gene expression that were consistent with known modes-of-action, demonstrating their suitability for use as reference chemicals for monitoring assay performance; and (4) ATRA produced phenotypic changes that were highly similar to other retinoic acid receptor activators (AM580, arotinoid acid) and some retinoid X receptor activators (bexarotene, methoprene acid). This phenotype was observed concurrently with autoregulation of the RARB gene. Both effects were prevented by pre-treating U-2 OS cells with pharmacological antagonists of their respective receptors. Thus, the observed phenotype could be considered characteristic of retinoic acid pathway activation in U-2 OS cells. These findings lay the groundwork for combinatorial screening of chemicals using HTTr and HTPP to generate complementary information for the first tier of a NAM-based chemical hazard evaluation strategy.

A phenomics approach for antiviral drug discovery.

BACKGROUND: The emergence and continued global spread of the current COVID-19 pandemic has highlighted the need for methods to identify novel or repurposed therapeutic drugs in a fast and effective way. Despite the availability of methods for the discovery of antiviral drugs, the majority tend to focus on the effects of such drugs on a given virus, its constituent proteins, or enzymatic activity, often neglecting the consequences on host cells. This may lead to partial assessment of the efficacy of the tested anti-viral compounds, as potential toxicity impacting the overall physiology of host cells may mask the effects of both viral infection and drug candidates. Here we present a method able to assess the general health of host cells based on morphological profiling, for untargeted phenotypic drug screening against viral infections. RESULTS: We combine Cell Painting with antibody-based detection of viral infection in a single assay. We designed an image analysis pipeline for segmentation and classification of virus-infected and non-infected cells, followed by extraction of morphological properties. We show that this methodology can successfully capture virus-induced phenotypic signatures of MRC-5 human lung fibroblasts infected with human coronavirus 229E (CoV-229E). Moreover, we demonstrate that our method can be used in phenotypic drug screening using a panel of nine host- and virus-targeting antivirals. Treatment with effective antiviral compounds reversed the morphological profile of the host cells towards a non-infected state. CONCLUSIONS: The phenomics approach presented here, which makes use of a modified Cell Painting protocol by incorporating an anti-virus antibody stain, can be used for the unbiased morphological profiling of virus infection on host cells. The method can identify antiviral reference compounds, as well as novel antivirals, demonstrating its suitability to be implemented as a strategy for antiviral drug repurposing and drug discovery.

Design, Synthesis, and Biological Evaluation of Chemically and Biologically Diverse Pyrroquinoline Pseudo Natural Products.

Natural product (NP) structures are a rich source of inspiration for the discovery of new biologically relevant chemical matter. In natural product inspired pseudo-NPs, NP-derived fragments are combined de novo in unprecedented arrangements. Described here is the design and synthesis of a 155-member pyrroquinoline pseudo-NP collection in which fragments characteristic of the tetrahydroquinoline and pyrrolidine NP classes are combined with eight different connectivities and regioisomeric arrangements. Cheminformatic analysis and biological evaluation of the compound collection by means of phenotyping in the morphological "cell painting" assay followed by principal component analysis revealed that the pseudo-NP classes are chemically diverse and that bioactivity patterns differ markedly, and are dependent on connectivity and regioisomeric arrangement of the fragments.

Multiple Chemical Features Impact Biological Performance Diversity of a Highly Active Natural Product-Inspired Library.

A systematic, diversity-oriented synthesis approach was employed to access a natural product-inspired flavonoid library with diverse chemical features, including chemical properties, scaffold, stereochemistry, and appendages. Using Cell Painting, the effects of these diversity elements were evaluated, and multiple chemical features that predict biological performance diversity were identified. Scaffold identity appears to be the dominant predictor of performance diversity, but stereochemistry and appendages also contribute to a lesser degree. In addition, the diversity of chemical properties contributed to performance diversity, and the driving chemical property was dependent on the scaffold. These results highlight the importance of key chemical features that may inform the creation of small-molecule, performance-diverse libraries to improve the efficiency and success of high-throughput screening campaigns.

Bioactivity screening of environmental chemicals using imaging-based high-throughput phenotypic profiling.

The present study adapted an existing high content imaging-based high-throughput phenotypic profiling (HTPP) assay known as "Cell Painting" for bioactivity screening of environmental chemicals. This assay uses a combination of fluorescent probes to label a variety of organelles and measures a large number of phenotypic features at the single cell level in order to detect chemical-induced changes in cell morphology. First, a small set of candidate phenotypic reference chemicals (n = 14) known to produce changes in the cellular morphology of U-2 OS cells were identified and screened at multiple time points in concentration-response format. Many of these chemicals produced distinct cellular phenotypes that were qualitatively similar to those previously described in the literature. A novel workflow for phenotypic feature extraction, concentration-response modeling and determination of in vitro thresholds for chemical bioactivity was developed. Subsequently, a set of 462 chemicals from the ToxCast library were screened in concentration-response mode. Bioactivity thresholds were calculated and converted to administered equivalent doses (AEDs) using reverse dosimetry. AEDs were then compared to effect values from mammalian toxicity studies. In many instances (68%), the HTPP-derived AEDs were either more conservative than or comparable to the in vivo effect values. Overall, we conclude that the HTPP assay can be used as an efficient, cost-effective and reproducible screening method for characterizing the biological activity and potency of environmental chemicals for potential use in in vitro-based safety assessments.

High content phenotypic screening identifies serotonin receptor modulators with selective activity upon breast cancer cell cycle and cytokine signaling pathways.

Heterogeneity in disease mechanisms between genetically distinct patients contributes to high attrition rates in late stage clinical drug development. New personalized medicine strategies aim to identify predictive biomarkers which stratify patients most likely to respond to a particular therapy. However, for complex multifactorial diseases not characterized by a single genetic driver, empirical approaches to identifying predictive biomarkers and the most promising therapies for personalized medicine are required. In vitro pharmacogenomics seeks to correlate in vitro drug sensitivity testing across panels of genetically distinct cell models with genomic, gene expression or proteomic data to identify predictive biomarkers of drug response. However, the vast majority of in vitro pharmacogenomic studies performed to date are limited to dose-response screening upon a single viability assay endpoint. In this article we describe the application of multiparametric high content phenotypic screening and the theta comparative cell scoring method to quantify and rank compound hits, screened at a single concentration, which induce a broad variety of divergent phenotypic responses between distinct breast cancer cell lines. High content screening followed by transcriptomic pathway analysis identified serotonin receptor modulators which display selective activity upon breast cancer cell cycle and cytokine signaling pathways correlating with inhibition of cell growth and survival. These methods describe a new evidence-led approach to rapidly identify compounds which display distinct response between different cell types. The results presented also warrant further investigation of the selective activity of serotonin receptor modulators upon breast cancer cell growth and survival as a potential drug repurposing opportunity.

Image-Based Morphological Profiling Identifies a Lysosomotropic, Iron-Sequestering Autophagy Inhibitor.

Chemical proteomics is widely applied in small-molecule target identification. However, in general it does not identify non-protein small-molecule targets, and thus, alternative methods for target identification are in high demand. We report the discovery of the autophagy inhibitor autoquin and the identification of its molecular mode of action using image-based morphological profiling in the cell painting assay. A compound-induced fingerprint representing changes in 579 cellular parameters revealed that autoquin accumulates in lysosomes and inhibits their fusion with autophagosomes. In addition, autoquin sequesters Fe2+ in lysosomes, resulting in an increase of lysosomal reactive oxygen species and ultimately cell death. Such a mechanism of action would have been challenging to unravel by current methods. This work demonstrates the potential of the cell painting assay to deconvolute modes of action of small molecules, warranting wider application in chemical biology.

A Scaffold-Diversity Synthesis of Biologically Intriguing Cyclic Sulfonamides.

A "branching-folding" synthetic strategy that affords a range of diverse cyclic benzo-sulfonamide scaffolds is presented. Whereas different annulation reactions on common ketimine substrates build the branching phase of the scaffold synthesis, a common hydrogenative ring-expansion method, facilitated by an increase of the ring-strain during the branching phase, led to sulfonamides bearing medium-sized rings in a folding pathway. Cell painting assay was successfully employed to identify tubulin targeting sulfonamides as novel mitotic inhibitors.

Design, Synthesis, and Phenotypic Profiling of Pyrano-Furo-Pyridone Pseudo Natural Products.

Natural products (NPs) inspire the design and synthesis of novel biologically relevant chemical matter, for instance through biology-oriented synthesis (BIOS). However, BIOS is limited by the partial coverage of NP-like chemical space by the guiding NPs. The design and synthesis of "pseudo NPs" overcomes these limitations by combining NP-inspired strategies with fragment-based compound design through de novo combination of NP-derived fragments to unprecedented compound classes not accessible through biosynthesis. We describe the development and biological evaluation of pyrano-furo-pyridone (PFP) pseudo NPs, which combine pyridone- and dihydropyran NP fragments in three isomeric arrangements. Cheminformatic analysis indicates that the PFPs reside in an area of NP-like chemical space not covered by existing NPs but rather by drugs and related compounds. Phenotypic profiling in a target-agnostic "cell painting" assay revealed that PFPs induce formation of reactive oxygen species and are structurally novel inhibitors of mitochondrial complex I.

Metabolic and phenotypic changes induced by PFAS exposure in two human hepatocyte cell models.

PFAS are ubiquitous industrial chemicals with known adverse health effects, particularly on the liver. The liver, being a vital metabolic organ, is susceptible to PFAS-induced metabolic dysregulation, leading to conditions such as hepatotoxicity and metabolic disturbances. In this study, we investigated the phenotypic and metabolic responses of PFAS exposure using two hepatocyte models, HepG2 (male cell line) and HepaRG (female cell line), aiming to define phenotypic alterations, and metabolic disturbances at the metabolite and pathway levels. The PFAS mixture composition was selected based on epidemiological data, covering a broad concentration spectrum observed in diverse human populations. Phenotypic profiling by Cell Painting assay disclosed predominant effects of PFAS exposure on mitochondrial structure and function in both cell models as well as effects on F-actin, Golgi apparatus, and plasma membrane-associated measures. We employed comprehensive metabolic characterization using liquid chromatography combined with high-resolution mass spectrometry (LC-HRMS). We observed dose-dependent changes in the metabolic profiles, particularly in lipid, steroid, amino acid and sugar and carbohydrate metabolism in both cells as well as in cell media, with HepaRG cell line showing a stronger metabolic response. In cells, most of the bile acids, acylcarnitines and free fatty acids showed downregulation, while medium-chain fatty acids and carnosine were upregulated, while the cell media showed different response especially in relation to the bile acids in HepaRG cell media. Importantly, we observed also nonmonotonic response for several phenotypic features and metabolites. On the pathway level, PFAS exposure was also associated with pathways indicating oxidative stress and inflammatory responses. Taken together, our findings on PFAS-induced phenotypic and metabolic disruptions in hepatocytes shed light on potential mechanisms contributing to the broader comprehension of PFAS-related health risks.