The ability to engineer complex multicellular systems has enormous potential to inform our understanding of biological processes and disease and alter the drug development process. Engineering living systems to emulate natural processes or to incorporate new functions relies on a detailed understanding of the biochemical, mechanical, and other cues between cells and between cells and their environment that result in the coordinated action of multicellular systems. On April 3–6, 2022, experts in the field met at the Keystone symposium “Engineering Multicellular Living Systems” to discuss recent advances in understanding how cells cooperate within a multicellular system, as well as recent efforts to engineer systems like organ-on-a-chip models, biological robots, and organoids. Given the similarities and common themes, this meeting was held in conjunction with the symposium “Organoids as Tools for Fundamental Discovery and Translation”.

While microscopy-based cellular assays, including microfluidics, have significantly advanced over the last several decades, there has not been concurrent development of widely-accessible techniques to analyze time-dependent microscopy data incorporating phenomena such as fluid flow and dynamic cell adhesion. As such, experimentalists typically rely on error-prone and time-consuming manual analysis, resulting in lost resolution and missed opportunities for innovative metrics. We present a user-adaptable toolkit packaged into the open-source, standalone Interactive Cellular assay Labeled Observation and Tracking Software (iCLOTS). We benchmark cell adhesion, single-cell tracking, velocity profile, and multiscale microfluidic-centric applications with blood samples, the prototypical biofluid specimen. Moreover, machine learning algorithms characterize previously imperceptible data groupings from numerical outputs. Free to download/use, iCLOTS addresses a need for a field stymied by a lack of analytical tools for innovative, physiologically-relevant assays of any design, democratizing use of well-validated algorithms for all end-user biomedical researchers who would benefit from advanced computational methods.

Understanding intracellular redox chemistry requires new tools for the site-specific visualization of intracellular oxidation. We have developed a spatially-resolved intracellular sensor of hydrogen peroxide, HyPer-Tau, for time-resolved imaging in live cells. This sensor consists of a hydrogen peroxide-sensing protein tethered to microtubules. We demonstrate the use of the HyPer-Tau sensor for three applications; dose-dependent response of human cells to exogenous hydrogen peroxide, a model immune response of mouse macrophages to stimulation by bacterial toxin, and a spatially-resolved response to localized delivery of hydrogen peroxide. These results demonstrate that HyPer-Tau can be used as an effective tool for tracking changes in spatially localized intracellular hydrogen peroxide and for future applications in redox signaling.

Quantifying the magnitude and dynamics of protein oxidation during cell signaling is technically challenging. Computational modeling provides tractable, quantitative methods to test hypotheses of redox mechanisms that may be simultaneously operative during signal transduction. The interleukin-4 (IL-4) pathway, which has previously been reported to induce reactive oxygen species and oxidation of PTP1B, may be controlled by several other putative mechanisms of redox regulation; widespread proteomic thiol oxidation observed via 2D redox differential gel electrophoresis upon IL-4 treatment suggests more than one redox-sensitive protein implicated in this pathway. Through computational modeling and a model selection strategy that relied on characteristic STAT6 phosphorylation dynamics of IL-4 signaling, we identified reversible protein tyrosine phosphatase (PTP) oxidation as the primary redox regulatory mechanism in the pathway. A systems-level model of IL-4 signaling was developed that integrates synchronous pan-PTP oxidation with ROS-independent mechanisms. The model quantitatively predicts the dynamics of IL-4 signaling over a broad range of new redox conditions, offers novel hypotheses about regulation of JAK/STAT signaling, and provides a framework for interrogating putative mechanisms involving receptor-initiated oxidation.

The side population (SP) assay, a technique used in cancer and stem cell research, assesses the activity of ABC transporters on Hoechst staining in the presence and absence of transporter inhibition, identifying SP and non-SP cell (NSP) subpopulations by differential staining intensity. The interpretation of the assay is complicated because the transporter-mediated mechanisms fail to account for cell-to-cell variability within a population or adequately control the direct role of transporter activity on staining intensity. We hypothesized that differences in dye kinetics at the single-cell level, such as ABCG2 transporter-mediated efflux and DNA binding, are responsible for the differential cell staining that demarcates SP/NSP identity. We report changes in A549 phenotype during time in culture and with TGFβ treatment that correlate with SP size. Clonal expansion of individually sorted cells re-established both SP and NSPs, indicating that SP membership is dynamic. To assess the validity of a purely kinetics-based interpretation of SP/NSP identity, we developed a computational approach that simulated cell staining within a heterogeneous cell population; this exercise allowed for the direct inference of the role of transporter activity and inhibition on cell staining. Our simulated SP assay yielded appropriate SP responses for kinetic scenarios in which high transporter activity existed in a portion of the cells and little differential staining occurred in the majority of the population. With our approach for single-cell analysis, we observed SP and NSP cells at both ends of a transporter activity continuum, demonstrating that features of transporter activity as well as DNA content are determinants of SP/NSP identity.

The initiation of heterogeneity within a population of phenotypically identical progenitors is a critical event for the onset of morphogenesis and differentiation patterning. Gap junction communication within multicellular systems produces complex networks of intercellular connectivity that result in heterogeneous distributions of intracellular signaling molecules. In this study, we investigate emergent systems-level behavior of the intercellular network within embryonic stem cell (ESC) populations and corresponding spatial organization during early neural differentiation. An agent-based model incorporates experimentally-determined parameters to yield complex transport networks for delivery of pro-differentiation cues between neighboring cells, reproducing the morphogenic trajectories during retinoic acid-accelerated mouse ESC differentiation. Furthermore, the model correctly predicts the delayed differentiation and preserved spatial features of the morphogenic trajectory that occurs in response to intercellular perturbation. These findings suggest an integral role of gap junction communication in the temporal coordination of emergent patterning during early differentiation and neural commitment of pluripotent stem cells.

Cell response to extracellular ligand is affected not only by ligand availability, but also by pre-existing cell-to-cell variability that enables a range of responses within a cell population. We developed a computational model that incorporates cell heterogeneity in order to investigate Jurkat T cell response to time dependent extracellular IL-2 stimulation. Our model predicted preferred timing of IL-2 oscillatory input for maximizing downstream intracellular STAT5 nuclear translocation. The modeled cytokine exposure was replicated experimentally through the use of a microfluidic platform that enabled the parallelized capture of dynamic single cell response to precisely delivered pulses of IL-2 stimulus. The in vitro results demonstrate that single cell response profiles vary with pulsatile IL-2 input at pre-equilibrium levels. These observations confirmed our model predictions that Jurkat cells have a preferred range of extracellular IL-2 fluctuations, in which downstream response is rapidly initiated. Further investigation into this filtering behavior could increase our understanding of how pre-existing cellular states within immune cell populations enable a systems response within a preferred range of ligand fluctuations, and whether the observed cytokine range corresponds to in vivo conditions.

Single-cell analysis has the potential to provide us with a host of new knowledge about biological systems, but it comes with the challenge of correctly interpreting the biological information. While emerging techniques have made it possible to measure inter-cellular variability at the transcriptome level, no consensus yet exists on the most appropriate method of data analysis of such single cell data. Methods for analysis of transcriptional data at the population level are well established but are not well suited to single cell analysis due to their dependence on population averages. In order to address this question, we have systematically tested combinations of methods for primary data analysis on single cell transcription data generated from two types of primary immune cells, neutrophils and T lymphocytes. Cells were obtained from healthy individuals, and single cell transcript expression data was obtained by a combination of single cell sorting and nanoscale quantitative real time PCR (qRT-PCR) for markers of cell type, intracellular signaling, and immune functionality. Gene expression analysis was focused on hierarchical clustering to determine the existence of cellular subgroups within the populations. Nine combinations of criteria for data exclusion and normalization were tested and evaluated. Bimodality in gene expression indicated the presence of cellular subgroups which were also revealed by data clustering. We observed evidence for two clearly defined cellular subtypes in the neutrophil populations and at least two in the T lymphocyte populations. When normalizing the data by different methods, we observed varying outcomes with corresponding interpretations of the biological characteristics of the cell populations. Normalization of the data by linear standardization taking into account technical effects such as plate effects, resulted in interpretations that most closely matched biological expectations. Single cell transcription profiling provides evidence of cellular subclasses in neutrophils and leukocytes that may be independent of traditional classifications based on cell surface markers. The choice of primary data analysis method had a substantial effect on the interpretation of the data. Adjustment for technical effects is critical to prevent misinterpretation of single cell transcript data.

The complexity of morphogenesis poses a fundamental challenge to understanding the mechanisms governing the formation of biological patterns and structures. Over the past century, numerous processes have been identified as critically contributing to morphogenetic events, but the interplay between the various components and aspects of pattern formation have been much harder to grasp. The combination of traditional biology with mathematical and computational methods has had a profound effect on our current understanding of morphogenesis and led to significant insights and advancements in the field. In particular, the theoretical concepts of reaction-diffusion systems and positional information, proposed by Alan Turing and Lewis Wolpert, respectively, dramatically influenced our general view of morphogenesis, although typically in isolation from one another. In recent years, agent-based modeling has been emerging as a consolidation and implementation of the two theories within a single framework. Agent-based models (ABMs) are unique in their ability to integrate combinations of heterogeneous processes and investigate their respective dynamics, especially in the context of spatial phenomena. In this review, we highlight the benefits and technical challenges associated with ABMs as tools for examining morphogenetic events. These models display unparalleled flexibility for studying various morphogenetic phenomena at multiple levels and have the important advantage of informing future experimental work, including the targeted engineering of tissues and organs.

Titanium dioxide nanoparticles (TiO 2 NPs), used as pigments and photocatalysts, are ubiquitous in our daily lives. Previous work has observed cellular oxidative stress in response to the UV-excitation of photocatalytic TiO 2 NPs. In comparison, most human exposure to TiO 2 NPs takes place in the dark, in the lung following inhalation or in the gut following consumption of TiO 2 NP food pigment. Our spectroscopic characterization shows that both photocatalytic and food grade TiO 2 NPs, in the dark, generate low levels of reactive oxygen species (ROS), specifically hydroxyl radicals and superoxides. These ROS oxidize serum proteins that form a corona of proteins on the NP surface. This protein layer is the interface between the NP and the cell. An oxidized protein corona triggers an oxidative stress response, detected with PCR and western blotting. Surface modification of TiO 2 NPs to increase or decrease surface defects correlates with ROS generation and oxidative stress, suggesting that NP surface defects, likely oxygen vacancies, are the underlying cause of TiO 2 NP-induced oxidative stress.