Single-Cell Receptor Quantification of an In Vitro Coculture Angiogenesis Model Reveals VEGFR, NRP1, Tie2, and PDGFR Regulation and Endothelial Heterogeneity

Understanding the effects of oxytocin receptor variants on OXT-OXT receptor binding: A mathematical model

Approximately half of U.S. women giving birth annually receive Pitocin, the synthetic form of oxytocin (OXT), yet its effective dose can vary significantly. This variability presents safety concerns due to unpredictable responses, which may lead to adverse outcomes for both mother and baby. To address the need for improved dosing, we developed a data-driven mathematical model to predict OXT receptor (OXTR) binding. Our study focuses on five prevalent OXTR variants (V45L, P108A, L206V, V281M, and E339K) and their impact on OXT-OXTR binding dynamics in two distinct cell types: human embryonic kidney cells (HEK293T), commonly used in experimental systems, and human myometrial smooth muscle cells, containing endogenous OXTR. We parameterized the model with cell-specific OXTR surface localization measurements. To strengthen the robustness of our study, we conducted a comprehensive meta-analysis of OXT- OXTR binding, enabling parameterization of our model with cell-specific OXT-OXTR binding kinetics (myometrial OXT-OXTR K d = 1.6 nM, kon = 6.8 × 10 5 M -1 min -1 , and koff = 0.0011 min -1 ). Our meta-analysis revealed significant homogeneity in OXT-OXTR affinity across experiments and species with a K d = 0.52 - 9.32 nM and mean K d = 1.48 ± 0.36 nM. Our model achieves several valuable insights into designing dosage strategies. First, we predicted that the OXTR complex reaches maximum occupancy at 10 nM OXT in myometrial cells and at 1 µM in HEK293T cells. This information is pivotal for guiding experimental design and data interpretation when working with these distinct cell types, emphasizing the need to consider effects for specific cell types when choosing OXTR-transfected cell lines. Second, our model recapitulated the significant effects of genetic variants for both experimental and physiologically relevant systems, with V281M and E339K substantially compromising OXT-OXTR binding capacity. These findings suggest the need for personalized oxytocin dosing based on individual genetic profiles to enhance therapeutic efficacy and reduce risks, especially in the context of labor and delivery. Third, we demonstrated the potential for rescuing the attenuated cell response observed in V281M and E339K variants by increasing the OXT dosage at specific, early time points. Cellular responses to OXT, including Ca 2+ release, manifest within minutes. Our model indicates that providing V281M- and E339K-expressing cells with doubled OXT dose during the initial minute of binding can elevate OXT-OXTR complex formation to levels comparable to wild-type OXTR. In summary, our study provides a computational framework for precision oxytocin dosing strategies, paving the way for personalized medicine.

Cellular subpopulations identified using an ensemble average of multiple dielectrophoresis measurements

While the average value measurement approach can successfully analyze and predict the general behavior and biophysical properties of an isogenic cell population, it fails when significant differences among individual cells are generated in the population by intracellular changes such as the cell cycle, or different cellular responses to certain stimuli. Detecting such single-cell differences in a cell population has remained elusive. Here, we describe an easy-to-implement and generalizable platform that measures the dielectrophoretic cross-over frequency of individual cells by decreasing measurement noise with a stochastic method and computing ensemble average statistics. This platform enables multiple, real-time, label-free detection of individual cells with significant dielectric variations over time within an isogenic cell population. Using a stochastic method in combination with the platform, we distinguished cell subpopulations from a mixture of drug-untreated and -treated isogenic cells. Furthermore, we demonstrate that our platform can identify drug-treated isogenic cells with different recovery rates.

Quantification of surface-localized and total oxytocin receptor in myometrial smooth muscle cells

Oxytocin acts through the oxytocin receptor (OXTR) to modulate uterine contractility. We previously identified OXTR genetic variants and showed that, in HEK293T cells, two of the OXTR protein variants localized to the cell surface less than wild-type OXTR. Here, we sought to measure OXTR in the more native human myometrial smooth muscle cell (HMSMC) line on both the cell-surface and across the whole cell, and used CRISPR editing to add an HA tag to the endogenous OXTR gene for anti-HA measurement. Quantitative flow cytometry revealed that these cells possessed 55,000 ± 3200 total OXTRs and 4900 ± 390 cell-surface OXTRs per cell. To identify any differential wild-type versus variant localization, we transiently transfected HMSMCs to exogenously express wild-type or variant OXTR with HA and green fluorescent protein tags. Total protein expression of wild-type OXTR and all tested variants were similar. However, the two variants with lower surface localization in HEK293T cells also presented lower surface localization in HMSMCs. Overall, we confirm the differential surface localization of variant OXTR in a more native cell type, and further demonstrate that the quantitative flow cytometry technique is adaptable to whole-cell measurements.

Axl and Vascular Endothelial Growth Factor Receptors Exhibit Variations in Membrane Localization and Heterogeneity Across Monolayer and Spheroid High-Grade Serous Ovarian Cancer Models

Vascular endothelial growth factor receptors (VEGFRs) and Axl are receptor tyrosine kinases (RTK) that are targeted in ovarian cancer therapy. Two-dimensional monolayer culture and three-dimensional spheroids are common models for RTK-targeted drug screening: monolayers are simple and economical while spheroids include several genetic and histological tumor features. RTK membrane localization dictates RTK signaling and drug response, however, it is not characterized in these models. We quantify plasma membrane RTK concentrations and show differential RTK abundance and heterogeneity in monolayers versus spheroids. We show VEGFR1 concentrations on the plasma membrane to be 10 times higher in OVCAR8 spheroids than in monolayers; OVCAR8 spheroids are more heterogeneous than monolayers, exhibiting a bimodal distribution of a low-Axl (6200/cell) and a high-Axl subpopulation (25,000/cell). In addition, plasma membrane Axl concentrations differ by 100 times between chemosensitive (OVCAR3) and chemoresistant (OVCAR8) cells and by 10 times between chemoresistant cell lines (OVCAR5 vs. OVCAR8). These systematic findings can guide ovarian cancer model selection for drug screening.

Transcriptomic analysis of 3D vasculature-on-a-chip reveals paracrine factors affecting vasculature growth and maturation

In vitro models of vasculature are of great importance for modelling vascular physiology and pathology. However, there is usually a lack of proper spatial patterning of interacting heterotypic cells in conventional vasculature dish models, which might confound results between contact and non-contact interactions. We use a microfluidic platform with structurally defined separation between human microvasculature and fibroblasts to probe their dynamic, paracrine interactions. We also develop a novel, versatile technique to retrieve cells embedded in extracellular matrix from the microfluidic device for downstream transcriptomic analysis, and uncover growth factor and cytokine expression profiles associated with improved vasculature growth. Paired receptor-ligand analysis further reveals paracrine signaling molecules that could be supplemented into the medium for vasculatures models where fibroblast coculture is undesirable or infeasible. These findings also provide deeper insights into the molecular cues for more physiologically relevant vascular mimicry and vascularized organoid model for clinical applications such as drug screening and disease modeling.

Absolute Quantification of Plasma Membrane Receptors Via Quantitative Flow Cytometry

Plasma membrane receptors are transmembrane proteins that initiate cellular response following the binding of specific ligands (e.g., growth factors, hormones, and cytokines). The abundance of plasma membrane receptors can be a diagnostic or prognostic biomarker in many human diseases. One of the best techniques for measuring plasma membrane receptors is quantitative flow cytometry (qFlow). qFlow employs fluorophore-conjugated antibodies against the receptors of interest and corresponding fluorophore-loaded calibration beads offers standardized and reproducible measurements of plasma membrane receptors. More importantly, qFlow can achieve absolute quantification of plasma membrane receptors when phycoerythrin (PE) is the fluorophore of choice. Here we describe a detailed qFlow protocol to obtain absolute receptor quantities on the basis of PE calibration. This protocol is foundational for many previous and ongoing studies in quantifying tyrosine kinase receptors and G-protein-coupled receptors with in vitro cell models and ex vivo cell samples.

A quantitative view on multivalent nanomedicine targeting

Although the concept of selective delivery has been postulated over 100 years ago, no targeted nanomedicine has been clinically approved so far. Nanoparticles modified with targeting ligands to promote the selective delivery of therapeutics towards a specific cell population have been extensively reported. However, the rational design of selective particles is still challenging. One of the main reasons for this is the lack of quantitative theoretical and experimental understanding of the interactions involved in cell targeting. In this review, we discuss new theoretical models and experimental methods that provide a quantitative view of targeting. We show the new advancements in multivalency theory enabling the rational design of super-selective nanoparticles. Furthermore, we present the innovative approaches to obtain key targeting parameters at the single-cell and single molecule level and their role in the design of targeting nanoparticles. We believe that the combination of new theoretical multivalent design and experimental methods to quantify receptors and ligands aids in the rational design and clinical translation of targeted nanomedicines.

Systems Biology Will Direct Vascular-Targeted Therapy for Obesity

Healthy adipose tissue expansion and metabolism during weight gain require coordinated angiogenesis and lymphangiogenesis. These vascular growth processes rely on the vascular endothelial growth factor (VEGF) family of ligands and receptors (VEGFRs). Several studies have shown that controlling vascular growth by regulating VEGF:VEGFR signaling can be beneficial for treating obesity; however, dysregulated angiogenesis and lymphangiogenesis are associated with several chronic tissue inflammation symptoms, including hypoxia, immune cell accumulation, and fibrosis, leading to obesity-related metabolic disorders. An ideal obesity treatment should minimize adipose tissue expansion and the advent of adverse metabolic consequences, which could be achieved by normalizing VEGF:VEGFR signaling. Toward this goal, a systematic investigation of the interdependency of vascular and metabolic systems in obesity and tools to predict personalized treatment ranges are necessary to improve patient outcomes through vascular-targeted therapies. Systems biology can identify the critical VEGF:VEGFR signaling mechanisms that can be targeted to regress adipose tissue expansion and can predict the metabolic consequences of different vascular-targeted approaches. Establishing a predictive, biologically faithful platform requires appropriate computational models and quantitative tissue-specific data. Here, we discuss the involvement of VEGF:VEGFR signaling in angiogenesis, lymphangiogenesis, adipogenesis, and macrophage specification – key mechanisms that regulate adipose tissue expansion and metabolism. We then provide useful computational approaches for simulating these mechanisms, and detail quantitative techniques for acquiring tissue-specific parameters. Systems biology, through computational models and quantitative data, will enable an accurate representation of obese adipose tissue that can be used to direct the development of vascular-targeted therapies for obesity and associated metabolic disorders.