Kinetics of Nanomedicine in Tumor Spheroid as an In Vitro Model System for Efficient Tumor-Targeted Drug Delivery With Insights From Mathematical Models

Supramolecular approaches for the treatment of hypoxic regions in tumours

Supramolecular chemistry provides a range of 'weak' intermolecular interactions that allow drugs and prodrugs to self-assemble. In the complex biological setting of blood and tumours, these interactions must be stable enough for efficient and selective drug delivery to the tumour site, but weak enough to allow the release of the cytotoxic load. The non-covalent nature of supramolecular interactions enables the detachment of smaller (pro)drug monomers that can penetrate cancer cells differently to the original nanoparticles. Hypoxic tumours show low oxygen levels due to poor vascularization, which poses challenges for drug delivery and generates biological resistances. Supramolecular building blocks specifically designed for hypoxic tumours offer targeted activation of prodrug self-assemblies, enhancing effectiveness against hypoxic cancer cells and hypoxic regions in tumours. This Review explores how supramolecular chemistry can improve (pro)drug delivery and activation in hypoxic tumours.

Fluorescent Molecular Probe for Imaging Hypoxia in 2D Cell Culture Monolayers and 3D Tumor Spheroids: The Cell Membrane Partition Model for Predicting Probe Distribution in a Spheroid

Compared to cultured 2D cell monolayers, 3D multicellular spheroids are more realistic tumor models. Nonetheless, spheroids remain under-utilized in preclinical research, in part, because there is a lack of fluorescence sensors that can noninvasively interrogate all the individual cells within a spheroid. This present study describes a deep-red fluorogenic molecular probe for microscopic imaging of cells that contain a high level of nitroreductase enzyme activity as a biomarker of cell hypoxia. A first-generation version of the probe produced "turn-on" fluorescence in a 2D cell monolayer under hypoxic conditions; however, it was not useful in a 3D multicellular tumor spheroid because it only accumulated in the peripheral cells. To guide the probe structural optimization process, an intuitive theoretical membrane partition model was conceived to predict how a dosed probe will distribute within a 3D spheroid. The model identifies three limiting molecular diffusion pathways that are determined by a probe's membrane partition properties. A lipophilic probe with high membrane affinity rapidly becomes trapped in the membranes of the peripheral cells. In contrast, a very hydrophilic probe molecule with negligible membrane affinity diffuses rapidly through the spheroid intercellular space and rarely enters the cells. However, a probe molecule with intermediate membrane affinity undergoes sequential diffusion in and out of cells and distributes to all the cells within a spheroid. Using the model as a predictive tool, a second-generation fluorescent probe was prepared with a smaller and more hydrophilic molecular structure, and optical sectioning using structured illumination or light sheet microscopy revealed roughly even probe diffusion throughout a tumor spheroid. The membrane permeation model is likely to be broadly applicable for the structural optimization of various classes of molecules and nanoparticles to enable even distribution within a tumor spheroid.

Protocols for translocation processes of flexible polymers through a pore using LAMMPS

Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) is an open-source, powerful simulator with a customizable platform for extensive Langevin dynamics simulations. Here, we present a protocol for using LAMMPS to develop coarse-grained models of polymeric systems with macromolecular crowding, an integral part of any soft matter or biophysical system. We describe steps for installing software, using LAMMPS basic commands and code, and translocating polymers. This protocol has potential applications in developing mathematical models for DNA sequencing, controlled drug delivery, and cellular transport processes. For complete details on the use and execution of this protocol, please refer to Garg et al.1.

Lipid nanoparticles as nano-Trojan-horses for siRNA delivery and gene-knockdown

The therapeutic messenger RNA strategies, such as those using small interfering RNAs, take several advantages (versatility, efficiency and selectivity) over plasmid DNA-based strategies. However, the challenge remains to find nanovectors capable of properly loading the genetic material, transporting it through troublesome environments, like a tumoral site, and delivering it into the cytoplasm of target cells. Here, lipid nanoparticles, consisting of a gemini cationic/neutral helper lipid mixture, are proposed as siRNA nanovector. Cells from cervical and brain cancer overexpressing the green fluorescent protein (GFP) were chosen to analyse the biological response as well as the efficiency and safety of the siRNA-loaded nanovector according to the cell phenotype. Flow cytometry and epifluorescence or confocal microscopy were used to follow the gene knockdown in these overexpressed cells. The effect of the nanovector on cellular proliferation was evaluated with cytotoxicity assays while their potential oxidative stress generation was determined by quantifying the generation of reactive oxygen species. To explore the mechanism of cellular uptake, different inhibitors of endocytic pathways were used during incubation with cells. Finally, nanovectors were incubated in 3D-grown cells (spheroids) to see whether they can penetrate the complex tumoral microenvironments, their efficiency to knockdown GFP expression being monitored by confocal microscopy.

Targeted drug delivery in neurodegenerative diseases: the role of nanotechnology

Neurodegenerative diseases, characterized by progressive neuronal loss and cognitive impairments, pose a significant global health challenge. This study explores the potential of nanotherapeutics as a promising approach to enhance drug delivery across physiological barriers, particularly the blood-brain barrier (BBB) and blood-cerebrospinal fluid barrier (B-CSFB). By employing nanoparticles, this research aims to address critical challenges in the diagnosis and treatment of conditions such as Alzheimer's, Parkinson's, and Huntington's diseases. The multifactorial nature of these disorders necessitates innovative solutions that leverage nanomedicine to improve drug solubility, circulation time, and targeted delivery while minimizing off-target effects. The findings underscore the importance of advancing nanomedicine applications to develop effective therapeutic strategies that can alleviate the burden of neurodegenerative diseases on individuals and healthcare systems.

Modified Au:Fe-Ni magnetic micromotors improve drug delivery and diagnosis in MCF-7 cells and spheroids

Nano/micromotors hold immense potential for revolutionizing drug delivery and detection systems, especially in the realm of cancer diagnosis and treatment, owing to their distinctive features, including precise propulsion, maneuverability, and meticulously designed surface modifications. In this study, we explore the capabilities of modified and magnetically driven micromotors as active drug delivery systems within 2D and 3D cell culture environments and cancer diagnosis. We synthesized gold (Au) and iron-nickel (Fe-Ni) metallic-based magnetic micromotors (Au:Fe-Ni MMs) through electrochemical methods, equipping them with functionalities for controlled doxorubicin (DOX) release and cancer cell recognition. In 2D and spheroids of MCF-7 adenocarcinoma cells, the Au segment of these micromotors was utilized to help DOX loading through poly(sodium-4-styrenesulfonate) (PSS) functionalization, and the attachment of antiHER2 antibodies for specific recognition. This innovative approach enabled controlled drug release within the cancerous microenvironment, coupled with magnetic (Fe-Ni) propulsion for biocompatible drug delivery to MCF-7 cells. Furthermore, antiHER2 immobilized Au:Fe-Ni MMs effectively interacted with receptors, capitalizing on the overexpression of HER2 antigens on MCF-7 cells. Encouraging outcomes were observed, particularly in spheroid models, underscoring the remarkable potential of these multifunctional micromotors for advancing intelligent drug delivery methodologies and diagnostic purposes.

Recent Update on Nanocarrier(s) as the Targeted Therapy for Breast Cancer

Despite ongoing advances in cancer therapy, the results for the treatment of breast cancer are not satisfactory. The advent of nanotechnology promises to be an essential tool to improve drug delivery effectiveness in cancer therapy. Nanotechnology provides an opportunity to enhance the treatment modality by preventing degradation, improving tumour targeting, and controlling drug release. Recent advances have revealed several strategies to prevent cancer metastasis using nano-drug delivery systems (NDDS). These strategies include the design of appropriate nanocarriers loaded with anti-cancer drugs that target the optimization of physicochemical properties, modulate the tumour microenvironment, and target biomimetic techniques. Nanocarriers have emerged as a preferential approach in the chemotropic treatment for breast cancer due to their pivotal role in safeguarding the therapeutic agents against degradation. They facilitate efficient drug concentration in targeted cells, surmount the resistance of drugs, and possess a small size. Nevertheless, these nanocarrier(s) have some limitations, such as less permeability across the barrier and low bioavailability of loaded drugs. To overcome these challenges, integrating external stimuli has been employed, encompassing infrared light, thermal stimulation, microwaves, and X-rays. Among these stimuli, ultrasound-triggered nanocarriers have gained significant attention due to their cost-effectiveness, non-invasive nature, specificity, ability to penetrate tissues, and capacity to deliver elevated drug concentrations to intended targets. This article comprehensively reviews recent advancements in different nanocarriers for breast cancer chemotherapy. It also delves into the associated hurdles and offers valuable insights into the prospective directions for this innovative field.

Quantifying Imaging Agent Binding and Dissociation in 3-D Cancer Spheroid Tissue Culture Using Paired-Agent Principles

Binding kinetics play an important role in cancer diagnosis and therapeutics. However, current methods of quantifying binding kinetics fail to consider the three-dimensional environment that drugs and imaging agents experience in biological tissue. In response, a methodology to assay agent binding and dissociation in 3-D tissue culture was developed using paired-agent molecular imaging principles. To test the methodology, the uptakes of ABY-029 (an IRDye 800CW-labeled epidermal growth factor receptor (EGFR)-targeted antibody mimetic) and IRDye-700DX carboxylate in 3-D spheroids were measured in four different human cancer cell lines throughout staining and rinsing. A compartment model (optimized for the application) was then fit to the kinetic curves of both imaging agents to estimate binding and dissociation rate constants of the EGFR-targeted ABY-029 agent. A statistically significant correlation was observed between apparent association rate constant (k3) and the receptor concentration experimentally and in simulations (r = 0.99, p < 0.05). A statistically significant difference was found between effective k3 (apparent rate constant of ABY-029 binding to EGFR) values for cell lines with varying levels of EGFR expression (p < 0.05), with no significant difference found between cell lines and controls for other fit parameters. Additionally, a similar binding affinity profile compared to a gold standard method was determined by this model. This low-cost methodology to quantify imaging agent or drug binding affinity in clinically relevant 3-D tumor spheroid models can be used to guide timing of imaging in molecular guided surgery and could have implications in drug development.

Effect of Polyplex Size on Penetration into Tumor Spheroids

Ovarian cancer is one of the most lethal gynecological cancers in the world. In recent years, nucleic acid (NA)-based formulations have been shown to be promising treatments for ovarian cancer, including tumor nodules. However, gene therapy is not that far advanced in clinical reality due to unfavorable physicochemical properties of the NAs, such as high molecular weight, poor cellular uptake, rapid degradation by nucleases, etc. One of the strategies used to overcome these drawbacks is the complexation of anionic NAs via electrostatic interactions with cationic polymers, resulting in the formation of so-called polyplexes. In this work, the role of the size of pDNA and siRNA polyplexes on their penetration into ovarian-cancer-based tumor spheroids was investigated. For this, a methoxypoly(ethylene glycol) poly(2-(dimethylamino)ethyl methacrylate) (mPEG-pDMAEMA) diblock copolymer was synthesized as a polymeric carrier for NA binding and condensation with either plasmid DNA (pDNA) or short interfering RNA (siRNA). When prepared in HEPES buffer (10 mM, pH 7.4) at a nitrogen/phosphate (N/P) charge ratio of 5 and pDNA polyplexes were formed with a size of 162 ± 11 nm, while siRNA-based polyplexes displayed a size of 25 ± 2 nm. The polyplexes had a slightly positive zeta potential of +7-8 mV in the same buffer. SiRNA and pDNA polyplexes were tracked in vitro into tumor spheroids, resembling in vivo avascular ovarian tumor nodules. For this purpose, reproducible spheroids were obtained by coculturing ovarian carcinoma cells with primary mouse embryonic fibroblasts in different ratios (5:2, 1:1, and 2:5). Penetration studies revealed that after 24 h of incubation, siRNA polyplexes were able to penetrate deeper into the homospheroids (composed of only cancer cells) and heterospheroids (cancer cells cocultured with fibroblasts) compared to pDNA polyplexes which were mainly located in the rim. The penetration of the polyplexes was slowed when increasing the fraction of fibroblasts present in the spheroids. Furthermore, in the presence of serum siRNA polyplexes encoding for luciferase showed a high cellular uptake in 2D cells resulting in ∼50% silencing of luciferase expression. Taken together, these findings show that self-assembled small siRNA polyplexes have good potential as a platform to test ovarian tumor nodulus penetration..

Heterotypic tumor spheroids: a platform for nanomedicine evaluation

Nanomedicine has emerged as a promising therapeutic approach, but its translation to the clinic has been hindered by the lack of cellular models to anticipate how tumor cells will respond to therapy. Three-dimensional (3D) cell culture models are thought to more accurately recapitulate key features of primary tumors than two-dimensional (2D) cultures. Heterotypic 3D tumor spheroids, composed of multiple cell types, have become more popular than homotypic spheroids, which consist of a single cell type, as a superior model for mimicking in vivo tumor heterogeneity and physiology. The stromal interactions demonstrated in heterotypic 3D tumor spheroids can affect various aspects, including response to therapy, cancer progression, nanomedicine penetration, and drug resistance. Accordingly, to design more effective anticancer nanomedicinal therapeutics, not only tumor cells but also stromal cells (e.g., fibroblasts and immune cells) should be considered to create a more physiologically relevant in vivo microenvironment. This review aims to demonstrate current knowledge of heterotypic 3D tumor spheroids in cancer research, to illustrate current advances in utilizing these tumor models as a novel and versatile platform for in vitro evaluation of nanomedicine-based therapeutics in cancer research, and to discuss challenges, guidelines, and future directions in this field.

Gene Expression Profile of 3D Spheroids in Comparison with 2D Cell Cultures and Tissue Strains of Diffuse High-Grade Gliomas

The use of relevant, accessible, and easily reproducible preclinical models of diffuse gliomas is a prerequisite for the development of successful therapeutic approaches to their treatment. Here we studied the gene expression profile of 3D spheroids in a comparison with 2D cell cultures and tissue strains of diffuse high-grade gliomas. Using real time PCR, we evaluated the expression of Gfap, Cd44, Pten, S100b, Vegfa, Hif1a, Sox2, Melk, Gdnf, and Mgmt genes playing an important role in the progression of gliomas and regulating tumor cell proliferation, adhesion, invasion, plasticity, apoptosis, DNA repair, and recruitment of tumor-associated cells. Gene expression analysis showed that 3D spheroids are more similar to tumor tissue strains by the expression levels of Gfap, Cd44, and Pten, while the expression levels of Hif1a and Sox2 in 3D spheroids did not differ from those of 2D cell cultures, the expression levels S100b and Vegfa in 3D spheroids was higher than in other models, and the expression levels of Melk, Gdnf, and Mgmt genes changed diversely. Thus, 3D spheroid model more closely mimics the tumor tissue than 2D cell culture, but still is not the most relevant, probably due to too small size of spheroids, which does not allow reproducing hypoxia and apoptotic and necrotic processes in the tumor tissue.

Spatiotemporally controlled drug delivery via photothermally driven conformational change of self-integrated plasmonic hybrid nanogels

BACKGROUND

Spatiotemporal regulation is one of the major considerations for developing a controlled and targeted drug delivery system to treat diseases efficiently. Light-responsive plasmonic nanostructures take advantage due to their tunable optical and photothermal properties by changing size, shape, and spatial arrangement.

RESULTS

In this study, self-integrated plasmonic hybrid nanogels (PHNs) are developed for spatiotemporally controllable drug delivery through light-driven conformational change and photothermally-boosted endosomal escape. PHNs are easily synthesized through the simultaneous integration of gold nanoparticles (GNPs), thermo-responsive poly (N-isopropyl acrylamide), and linker molecules during polymerization. Wave-optic simulations reveal that the size of the PHNs and the density of the integrated GNPs are crucial factors in modulating photothermal conversion. Several linkers with varying molecular weights are inserted for the optimal PHNs, and the alginate-linked PHN (A-PHN) achieves more than twofold enhanced heat conversion compared with others. Since light-mediated conformational changes occur transiently, drug delivery is achieved in a spatiotemporally controlled manner. Furthermore, light-induced heat generation from cellular internalized A-PHNs enables pinpoint cytosolic delivery through the endosomal rupture. Finally, the deeper penetration for the enhanced delivery efficiency by A-PHNs is validated using multicellular spheroid.

CONCLUSION

This study offers a strategy for synthesizing light-responsive nanocarriers and an in-depth understanding of light-modulated site-specific drug delivery.

Controlling Cell Organization in 3D Coculture Spheroids Using DNA Interactions

The role of stromal and immune cells in transforming the tumor microenvironment is a key consideration in understanding tumor cell behavior and anticancer drug development. To better model these systems in vitro, 3D coculture tumor spheroids have been engineered using a variety of techniques including centrifugation to microwells, hanging drop, low adhesion cultures, and culture of cells in a microfluidic platform. Aside from using bioprinting, however, it has remained more challenging to direct the spatial organization of heterotypic cells in standalone 3D spheroids. To address this, we present an in vitro 3D coculture tumor model where we modulated the interactions between cancer cells and fibroblasts through DNA hybridization. When native heterotypic cells are simply mixed, the cell aggregates typically show cell sorting behavior to form phase separated structures composed of single cell types. In this work, we demonstrate that when MDA-MB-468 breast cancer and NIH/3T3 fibroblasts are directed to associate via complementary DNA, a uniform distribution of the two cell types within a single spheroid was observed. In contrast, in the absence of specific DNA interactions between the cancer cells and fibroblasts, individual clusters of the NIH/3T3 cells formed in each spheroid due to cell sorting. To better understand the effect of heterotypic cell organization on either cell-cell contacts or matrix protein production, the spheroids were further stained with anti-E-cadherin and antifibronectin antibodies. While the amounts of E-cadherin appeared to be similar between the spheroids, a significantly higher amount of fibronectin secretion was observed in the coculture spheroids with uniform mixing of two cell types. This result showed that different heterotypic cell distributions within 3D architecture can influence the ECM protein production that can again alter the properties of the tumor or tumor microenvironment. The present study thus describes the use of DNA templating to direct the organization of cells in coculture spheroids, which can provide mechanistic biological insight into how heterotypic distribution in tumor spheroids can influence tumor progression, metastasis, and drug resistance.

Efficient Synergistic Antibacterial Activity of α-MSH Using Chitosan-Based Versatile Nanoconjugates

The application of antimicrobial peptides has emerged as an alternative therapeutic tool to encounter against multidrug resistance of different pathogenic organisms. α-Melanocyte stimulating hormone (α-MSH), an endogenous neuropeptide, is found to be efficient in eradicating infection of various kinds of Staphylococcus aureus, including methicillin-resistant Staphylococcus aureus (MRSA). However, the chemical stability and efficient delivery of these biopharmaceuticals (i.e., α-MSH) to bacterial cells with a significant antibacterial effect remains a key challenge. To address this issue, we have developed a chitosan-cholesterol polymer using a single-step, one-pot, and simple chemical conjugation technique, where α-MSH is loaded with a significantly high amount (37.7%), and the final product is obtained as chitosan-cholesterol α-MSH polymer-drug nanoconjugates. A staphylococcal growth inhibition experiment was performed using chitosan-cholesterol α-MSH and individual controls. α-MSH and chitosan-cholesterol both show bacterial growth inhibition by a magnitude of 50 and 79%, respectively. The killing efficiency of polymer-drug nanoconjugates was very drastic, and almost no bacterial colony was observed (∼100% inhibition) after overnight incubation. Phenotypic alternation was observed in the presence of α-MSH causing changes in the cell structure and shape, indicating stress on Staphylococcus aureus. As a further consequence, vigorous cell lysis with concomitant release of the cellular material in the nearby medium was observed after treatment of chitosan-cholesterol α-MSH nanoconjugates. This vigorous lysis of the cell structure is associated with extensive aggregation of the bacterial cells evident in scanning electron microscopy (SEM). The dose-response experiment was performed with various concentrations of chitosan-cholesterol α-MSH nanoconjugates to decipher the degree of the bactericidal effect. The concentration of α-MSH as low as 1 pM also shows significant inhibition of bacterial growth (∼40% growth inhibition) of Staphylococcus aureus. Despite playing an important role in inhibiting bacterial growth, our investigation on hemolytic assay shows that chitosan-cholesterol α-MSH is significantly nontoxic at a wide range of concentrations. In a nutshell, our analysis demonstrated novel antimicrobial activity of nanoparticle-conjugated α-MSH, which could be used as future therapeutics against multidrug-resistant Staphylococcus aureus and other types of bacterial cells.

Glioblastoma immuno-endothelial multicellular microtissue as a 3D in vitro evaluation tool of anti-cancer nano-therapeutics

Despite being the most prevalent and lethal type of adult brain cancer, glioblastoma (GBM) remains intractable. Promising anti-GBM nanoparticle (NP) systems have been developed to improve the anti-cancer performance of difficult-to-deliver therapeutics, with particular emphasis on tumor targeting strategies. However, current disease modeling toolboxes lack close-to-native in vitro models that emulate GBM microenvironment and bioarchitecture, thus partially hindering translation due to poorly predicted clinical responses. Herein, human GBM heterotypic multicellular tumor microtissues (MCTMs) are generated through high-throughput 3D modeling of U-251 MG tumor cells, tissue differentiated macrophages isolated from peripheral monocytes, and brain microvascular primary endothelial cells. GBM MCTMs mimicked tumor spatial organization, extracellular matrix production and necrosis areas. The bioactivity of a model drug, docetaxel (DTX), and of tumor-targeted DTX-loaded polymeric NPs with a surface L-Histidine moiety (H-NPs), were assessed in the MCTMs. MCTMs cell uptake and anti-proliferative effect was 8- and 3-times higher for H-NPs, respectively, compared to the non-targeted NPs and to free DTX. H-NPs provided a decrease of MCTMs anti-inflammatory M2-macrophages, while increasing their pro-inflammatory M1 counterparts. Moreover, H-NPs showed a particular biomolecular signature through reduced secretion of an array of medium cytokines (IFN-γ, IL-1β, IL-1Ra, IL-6, IL-8, TGF-β). Overall, MCTMs provide an in vitro biomimetic model to recapitulate key cellular and structural features of GBM and improve in vivo drug response predictability, fostering future clinical translation of anti-GBM nano-therapeutic strategies.

Time-averaging and nonergodicity of reset geometric Brownian motion with drift

How do near-bankruptcy events in the past affect the dynamics of stock-market prices in the future? Specifically, what are the long-time properties of a time-local exponential growth of stock-market prices under the influence of stochastically occurring economic crashes? Here, we derive the ensemble- and time-averaged properties of the respective "economic" or geometric Brownian motion (GBM) with a nonzero drift exposed to a Poissonian constant-rate price-restarting process of "resetting." We examine-based both on thorough analytical calculations and on findings from systematic stochastic computer simulations-the general situation of reset GBM with a nonzero [positive] drift and for all special cases emerging for varying parameters of drift, volatility, and reset rate in the model. We derive and summarize all short- and long-time dependencies for the mean-squared displacement (MSD), the variance, and the mean time-averaged MSD (TAMSD) of the process of Poisson-reset GBM under the conditions of both rare and frequent resetting. We consider three main regions of model parameters and categorize the crossovers between different functional behaviors of the statistical quantifiers of this process. The analytical relations are fully supported by the results of computer simulations. In particular, we obtain that Poisson-reset GBM is a nonergodic stochastic process, with generally MSD(Δ)≠TAMSD(Δ) and Variance(Δ)≠TAMSD(Δ) at short lag times Δ and for long trajectory lengths T. We investigate the behavior of the ergodicity-breaking parameter in each of the three regions of parameters and examine its dependence on the rate of reset at Δ/T≪1. Applications of these theoretical results to the analysis of prices of reset-containing options are pertinent.

Intrinsic Differences in Spatiotemporal Organization and Stromal Cell Interactions Between Isogenic Lung Cancer Cells of Epithelial and Mesenchymal Phenotypes Revealed by High-Dimensional Single-Cell Analysis of Heterotypic 3D Spheroid Models

The lack of inadequate preclinical models remains a limitation for cancer drug development and is a primary contributor to anti-cancer drug failures in clinical trials. Heterotypic multicellular spheroids are three-dimensional (3D) spherical structures generated by self-assembly from aggregates of two or more cell types. Compared to traditional monolayer cell culture models, the organization of cells into a 3D tissue-like structure favors relevant physiological conditions with chemical and physical gradients as well as cell-cell and cell-extracellular matrix (ECM) interactions that recapitulate many of the hallmarks of cancer in situ. Epidermal growth factor receptor (EGFR) mutations are prevalent in non-small cell lung cancer (NSCLC), yet various mechanisms of acquired resistance, including epithelial-to-mesenchymal transition (EMT), limit the clinical benefit of EGFR tyrosine kinase inhibitors (EGFRi). Improved preclinical models that incorporate the complexity induced by epithelial-to-mesenchymal plasticity (EMP) are urgently needed to advance new therapeutics for clinical NSCLC management. This study was designed to provide a thorough characterization of multicellular spheroids of isogenic cancer cells of various phenotypes and demonstrate proof-of-principle for the applicability of the presented spheroid model to evaluate the impact of cancer cell phenotype in drug screening experiments through high-dimensional and spatially resolved imaging mass cytometry (IMC) analyses. First, we developed and characterized 3D homotypic and heterotypic spheroid models comprising EGFRi-sensitive or EGFRi-resistant NSCLC cells. We observed that the degree of EMT correlated with the spheroid generation efficiency in monocultures. In-depth characterization of the multicellular heterotypic spheroids using immunohistochemistry and high-dimensional single-cell analyses by IMC revealed intrinsic differences between epithelial and mesenchymal-like cancer cells with respect to self-sorting, spatiotemporal organization, and stromal cell interactions when co-cultured with fibroblasts. While the carcinoma cells harboring an epithelial phenotype self-organized into a barrier sheet surrounding the fibroblasts, mesenchymal-like carcinoma cells localized to the central hypoxic and collagen-rich areas of the compact heterotypic spheroids. Further, deep-learning-based single-cell segmentation of IMC images and application of dimensionality reduction algorithms allowed a detailed visualization and multiparametric analysis of marker expression across the different cell subsets. We observed a high level of heterogeneity in the expression of EMT markers in both the carcinoma cell populations and the fibroblasts. Our study supports further application of these models in pre-clinical drug testing combined with complementary high-dimensional single-cell analyses, which in turn can advance our understanding of the impact of cancer-stroma interactions and epithelial phenotypic plasticity on innate and acquired therapy resistance in NSCLC.