Decoding information in cell shape

Insights gained from computational modeling of YAP/TAZ signaling for cellular mechanotransduction

YAP/TAZ signaling pathway is regulated by a multiplicity of feedback loops, crosstalk with other pathways, and both mechanical and biochemical stimuli. Computational modeling serves as a powerful tool to unravel how these different factors can regulate YAP/TAZ, emphasizing biophysical modeling as an indispensable tool for deciphering mechanotransduction and its regulation of cell fate. We provide a critical review of the current state-of-the-art of computational models focused on YAP/TAZ signaling.

A bypass flow model to study endothelial cell mechanotransduction across diverse flow environments

Disturbed flow is one of the pathological initiators of endothelial dysfunction in intimal hyperplasia (IH) which is commonly seen in vascular bypass grafts, and arteriovenous fistulas. Various in vitro disease models have been designed to simulate the hemodynamic conditions found in the vasculature. Nonetheless, prior investigations have encountered challenges in establishing a robust disturbed flow model, primarily attributed to the complex bifurcated geometries and distinctive flow dynamics. In the present study, we aim to address this gap by introducing an in vitro bypass flow model capable of inducing disturbed flow and other hemodynamics patterns through a pulsatile flow in the same model. To assess the model's validity, we employed computational fluid dynamics (CFD) to simulate hemodynamics and compared the morphology and functions of human umbilical venous endothelial cells (HUVECs) under disturbed flow conditions to those in physiological flow or stagnant conditions. CFD analysis revealed the generation of disturbed flow within the model, pinpointing the specific location in the channel where the effects of disturbed flow were observed. High-content screening, a single-cell morphological profile assessment, demonstrated that HUVECs in the disturbed flow area exhibited random orientation, and morphological features were significantly distinct compared to cells in the physiological flow or stagnant condition after a two days of flow exposure. Furthermore, HUVECs exposed to disturbed flow underwent extensive remodeling of the adherens junctions and expressed higher levels of endothelial cell activation markers compared to other hemodynamic conditions. In conclusion, our in vitro bypass flow model provides a robust platform for investigating the associations between disturbed flow pattern and vascular diseases.

Biophysical Modeling of Synaptic Plasticity

Dendritic spines are small, bulbous compartments that function as postsynaptic sites and undergo intense biochemical and biophysical activity. The role of the myriad signaling pathways that are implicated in synaptic plasticity is well studied. A recent abundance of quantitative experimental data has made the events associated with synaptic plasticity amenable to quantitative biophysical modeling. Spines are also fascinating biophysical computational units because spine geometry, signal transduction, and mechanics work in a complex feedback loop to tune synaptic plasticity. In this sense, ideas from modeling cell motility can inspire us to develop multiscale approaches for predictive modeling of synaptic plasticity. In this article, we review the key steps in postsynaptic plasticity with a specific focus on the impact of spine geometry on signaling, cytoskeleton rearrangement, and membrane mechanics. We summarize the main experimental observations and highlight how theory and computation can aid our understanding of these complex processes.

Direct GPCR-EGFR interaction enables synergistic membrane-to-nucleus information transfer

We addressed the heteromerization of the epidermal growth factor receptor (EGFR) with G-protein coupled receptors (GPCR) on the basis of angiotensin-II-receptor-subtype-1(AT1R)-EGFR interaction as proof-of-concept and show its functional relevance during synergistic nuclear information transfer, beyond ligand-dependent EGFR transactivation. Following in silico modelling, we generated EGFR-interaction deficient AT1R-mutants and compared them to AT1R-wildtype. Receptor interaction was assessed by co-immunoprecipitation (CoIP), Förster resonance energy transfer (FRET) and fluorescence-lifetime imaging microscopy (FLIM). Changes in cell morphology, ERK1/2-phosphorylation (ppERK1/2), serum response factor (SRF)-activation and cFOS protein expression were determined by digital high content microscopy at the single cell level. FRET, FLIM and CoIP confirmed the physical interaction of AT1R-wildtype with EGFR that was strongly reduced for the AT1R-mutants. Responsiveness of cells transfected with AT1R-WT or -mutants to angiotensin II or EGF was similar regarding changes in cell circularity, ppERK1/2 (direct and by ligand-dependent EGFR-transactivation), cFOS-expression and SRF-activity. By contrast, the EGFR-AT1R-synergism regarding these parameters was completely absent for in the interaction-deficient AT1R mutants. The results show that AT1R-EGFR heteromerisation enables AT1R-EGFR-synergism on downstream gene expression regulation, modulating the intensity and the temporal pattern of nuclear AT1R/EGFR-information transfer. Furthermore, remote EGFR transactivation, via ligand release or cytosolic tyrosine kinases, is not sufficient for the complete synergistic control of gene expression.

Novel Kraft Softwood Lignin-Derived Carbon Quantum Dots: Synthesis, Characterization, and In Vitro Cytocompatibility

Carbon quantum dots (CQDs) have been investigated for biomedical applications in medical imaging due to their fluorescent properties, overall long-term stability, and excellent cytocompatibility and biocompatibility. Lignin is an organic polymer in the tissues of woody plants. It is also considered a byproduct of the wood and pulp industries. Hence, it presents as a renewable source of carbon nanoparticles. In this study, we report the synthesis and material and biological characterization of two colloidal suspensions of CQDs in water derived from lignin-based carbon. One was the native form of CQDs derived from lignin carbon, and the second was doped with nitrogen to evaluate material differences. Material characterization was carried out using various commonly used techniques, including Fourier transform infrared spectroscopy (FTIR), emission and absorbance spectra, zeta potential, transmission electron microscopy (TEM), and scanning electron microscopy (SEM). Thin films of CQDs were formed on glass and silicon substrates to assess the in vitro cytocompatibility with human mesenchymal stem cells (hMSCs). Observations suggest that the two forms of CQDs promote cell attachment within 24 h and sustain it for at least 7 days. The overall structure and shape of cells suggest a lack of any adverse or toxic effects of CQDs. The data lay down the novel foundation to support the use of lignin-derived CQDs in tissue engineering applications.

A Computational Model of Mechanical Stretching of Cultured Cells on a Flexible Membrane

Mechanical forces applied to cells are known to regulate a wide variety of biological processes. Recent studies have supported that mechanical forces can cause nuclear deformation, leading to significant alterations in the gene expression and chromatin landscape of the cell. While the stresses and strains applied to cells is it is often known or controlled experimentally on a macroscopic length scale, it is often unclear what the actual forces and displacements are at the microscopic level of the cell. In this work, we created a model of cell deformation during application of mechanical stretch to cultured cells growth on a flexible membrane. This configuration is commonly used is in experimental studies as a means to apply controlled mechanical strains to adherent cultured cells. The parameters used in the study were used for application of strain to a mesenchymal stem cell stretched on a membrane. computational model was created to simulate the stresses and strains within the cell under a variety of stain amplitudes, waveforms and frequencies of mechanical loading with the range of commonly used experimental systems. The results demonstrate the connection between mechanical loading parameters applied through the flexible membrane and the resulting stresses and strains within the cell and nucleus. Using a viscoelastic model of chromatin, we connected the results provide to a rough model of resulting deformation within chromatin from the forces applied to the nucleus. Overall, the model is useful in providing insight between experimentally applied mechanical forces and the actual forces within the cell to better interpret the results of experimental studies.

Statement of Significance

In this work, we created a computational model of the mechanical stretching of cell on a flexible membrane under cyclic mechanical loading. This model provides insight into the forces and displacements inside of cell that result from that application of stretch. As many experiments use this set up, our work is relevant to interpreting many studies that use mechanical stretch to stimulate mechanotransduction.

Hierarchical Surface Pattern on Ni-Free Ti-Based Bulk Metallic Glass to Control Cell Interactions

Ni-free Ti-based bulk metallic glasses (BMGs) are exciting materials for biomedical applications because of their outstanding biocompatibility and advantageous mechanical properties. The glassy nature of BMGs allows them to be shaped and patterned via thermoplastic forming (TPF). This work demonstrates the versatility of the TPF technique to create micro- and nano-patterns and hierarchical structures on Ti40Zr10Cu34Pd14Sn2 BMG. Particularly, a hierarchical structure fabricated by a two-step TPF process integrates 400 nm hexagonal close-packed protrusions on 2.5 µm square protuberances while preserving the advantageous mechanical properties from the as-cast material state. The correlations between thermal history, structure, and mechanical properties are explored. Regarding biocompatibility, Ti40Zr10Cu34Pd14Sn2 BMGs with four surface topographies (flat, micro-patterned, nano-patterned, and hierarchical-structured surfaces) are investigated using Saos-2 cell lines. Alamar Blue assay and live/dead analysis show that all tested surfaces have good cell proliferation and viability. Patterned surfaces are observed to promote the formation of longer filopodia on the edge of the cytoskeleton, leading to star-shaped and dendritic cell morphologies compared with the flat surface. In addition to potential implant applications, TPF-patterned Ti-BMGs enable a high level of order and design flexibility on the surface topography, expanding the available toolbox for studying cell behavior on rigid and ordered surfaces.

Cell-TIMP: Cellular Trajectory Inference based on Morphological Parameter

Cellular morphology, shaped by various genetic and environmental influences, is pivotal to studying experimental cell biology, necessitating precise measurement and analysis techniques. Traditional approaches, which rely on geometric metrics derived from stained images, encounter obstacles stemming from both the imaging and analytical domains. Staining processes can disrupt the cell's natural state and diminish accuracy due to photobleaching, while conventional analysis techniques, which categorize cells based on shape to discern pathophysiological conditions, often fail to capture the continuous and asynchronous nature of biological processes such as cell differentiation, immune responses, and cancer progression. In this work, we propose the use of quantitative phase imaging for morphological assessment due to its label-free nature. For analysis, we repurposed the genomic analysis toolbox to perform trajectory inference analysis purely based on morphology information. We applied the developed framework to study the progression of leukemia and breast cancer metastasis. Our approach revealed a clear pattern of morphological evolution tied to the diseases' advancement, highlighting the efficacy of our method in identifying functionally significant shape changes where conventional techniques falter. This advancement offers a fresh perspective on analyzing cellular morphology and holds significant potential for the broader research community, enabling a deeper understanding of complex biological dynamics.

Understanding Macrophage-Tumor Interactions: Insights from Single-Cell Behavior Monitoring in a Sessile Microdroplet System

Interaction between tumor-associated macrophages and tumor cells is crucial for tumor development, metastasis, and the related immune process. However, the macrophages are highly heterogeneous spanning from anti-tumorigenic to pro-tumorigenic, which needs to be understood at the single-cell level. Herein, a sessile microdroplet system designed for monitoring cellular behavior and analyzing intercellular interaction, demonstrated with macrophage-tumor cell pairs is presented. An automatic procedure based on the inkjet printing method is utilized for the precise pairing and co-encapsulation of heterotypic cells within picoliter droplets. The sessile nature of microdroplets ensures controlled fusion and provides stable environments conducive to adherent cell culture. The nitric oxide generation and morphological changes over incubation are explored to reveal the complicated interactions from a single-cell perspective. The immune response of macrophages under distinct cellular microenvironments is recorded. The results demonstrate that the tumor microenvironment displays a modulating role in polarizing macrophages from anti-tumorigenic into pro-tumorigenic phenotype. The approach provides a versatile and compatible platform to investigate intercellular interaction at the single-cell level, showing promising potential for advancing single-cell behavior studies.

"Patchiness" in mechanical stiffness across a tumor as an early-stage marker for malignancy

Mechanical phenotyping of tumors, either at an individual cell level or tumor cell population level is gaining traction as a diagnostic tool. However, the extent of diagnostic and prognostic information that can be gained through these measurements is still unclear. In this work, we focus on the heterogeneity in mechanical properties of cells obtained from a single source such as a tissue or tumor as a potential novel biomarker. We believe that this heterogeneity is a conventionally overlooked source of information in mechanical phenotyping data. We use mechanics-based in-silico models of cell-cell interactions and cell population dynamics within 3D environments to probe how heterogeneity in cell mechanics drives tissue and tumor dynamics. Our simulations show that the initial heterogeneity in the mechanical properties of individual cells and the arrangement of these heterogenous sub-populations within the environment can dictate overall cell population dynamics and cause a shift towards the growth of malignant cell phenotypes within healthy tissue environments. The overall heterogeneity in the cellular mechanotype and their spatial distributions is quantified by a "patchiness" index, which is the ratio of the global to local heterogeneity in cell populations. We observe that there exists a threshold value of the patchiness index beyond which an overall healthy population of cells will show a steady shift towards a more malignant phenotype. Based on these results, we propose that the "patchiness" of a tumor or tissue sample, can be an early indicator for malignant transformation and cancer occurrence in benign tumors or healthy tissues. Additionally, we suggest that tissue patchiness, measured either by biochemical or biophysical markers, can become an important metric in predicting tissue health and disease likelihood just as landscape patchiness is an important metric in ecology.

Order from chaos: cellular asymmetries explained with modelling

Molecules inside cells are subject to physical forces and undergo biochemical interactions, continuously changing their physical properties and dynamics. Despite this, cells achieve highly ordered molecular patterns that are crucial to regulate various cellular functions and to specify cell fate. In the Caenorhabditis elegans one-cell embryo, protein asymmetries are established in the narrow time window of a cell division. What are the mechanisms that allow molecules to establish asymmetries, defying the randomness imposed by Brownian motion? Mathematical and computational models have paved the way to the understanding of protein dynamics up to the 'single-molecule level' when resolution represents an issue for precise experimental measurements. Here we review the models that interpret cortical and cytoplasmic asymmetries in the one-cell C. elegans embryo.

Cristae formation is a mechanical buckling event controlled by the inner mitochondrial membrane lipidome

Cristae are high-curvature structures in the inner mitochondrial membrane (IMM) that are crucial for ATP production. While cristae-shaping proteins have been defined, analogous lipid-based mechanisms have yet to be elucidated. Here, we combine experimental lipidome dissection with multi-scale modeling to investigate how lipid interactions dictate IMM morphology and ATP generation. When modulating phospholipid (PL) saturation in engineered yeast strains, we observed a surprisingly abrupt breakpoint in IMM topology driven by a continuous loss of ATP synthase organization at cristae ridges. We found that cardiolipin (CL) specifically buffers the inner mitochondrial membrane against curvature loss, an effect that is independent of ATP synthase dimerization. To explain this interaction, we developed a continuum model for cristae tubule formation that integrates both lipid and protein-mediated curvatures. This model highlighted a snapthrough instability, which drives IMM collapse upon small changes in membrane properties. We also showed that cardiolipin is essential in low-oxygen conditions that promote PL saturation. These results demonstrate that the mechanical function of cardiolipin is dependent on the surrounding lipid and protein components of the IMM.

Small Plasma Membrane-Targeted Fluorescent Dye for Long-Time Imaging and Protein Degradation Analyses

Plasma membrane (PM)-targeted fluorescent dyes have become an important tool to visualize morphological and dynamic changes in the cell membrane. However, most of these PM dyes are either too large and thus might potentially perturb the membrane and affect its functions or exhibit a short retention time on the cell membrane. The rapid internalization problem is particularly severe for PM dyes based on cationic and neutral hydrophobic fluorescent dyes, which can be easily transported into the cells by transmembrane potential and passive diffusion mechanisms. In this paper, we report a small but highly specific PM fluorescent dye, PM-1, which exhibits a very long retention time on the plasma membrane with a half-life of approximately 15 h. For biological applications, we demonstrated that PM-1 can be used in combination with protein labeling probes to study ectodomain shedding and endocytosis processes of cell surface proteins and successfully demonstrated that native transmembrane human carbonic anhydrase IX (hCAIX) is degraded via the ectodomain shedding mechanism. In contrast, hCAIX undergoes endocytic degradation in the presence of sheddase inhibitors. We believe that PM-1 can be a versatile tool to provide detailed insights into the dynamic processes of the cell surface proteins.

Heterogeneous off-target impact of ion-channel deletion on intrinsic properties of hippocampal model neurons that self-regulate calcium

How do neurons that implement cell-autonomous self-regulation of calcium react to knockout of individual ion-channel conductances? To address this question, we used a heterogeneous population of 78 conductance-based models of hippocampal pyramidal neurons that maintained cell-autonomous calcium homeostasis while receiving theta-frequency inputs. At calcium steady-state, we individually deleted each of the 11 active ion-channel conductances from each model. We measured the acute impact of deleting each conductance (one at a time) by comparing intrinsic electrophysiological properties before and immediately after channel deletion. The acute impact of deleting individual conductances on physiological properties (including calcium homeostasis) was heterogeneous, depending on the property, the specific model, and the deleted channel. The underlying many-to-many mapping between ion channels and properties pointed to ion-channel degeneracy. Next, we allowed the other conductances (barring the deleted conductance) to evolve towards achieving calcium homeostasis during theta-frequency activity. When calcium homeostasis was perturbed by ion-channel deletion, post-knockout plasticity in other conductances ensured resilience of calcium homeostasis to ion-channel deletion. These results demonstrate degeneracy in calcium homeostasis, as calcium homeostasis in knockout models was implemented in the absence of a channel that was earlier involved in the homeostatic process. Importantly, in reacquiring homeostasis, ion-channel conductances and physiological properties underwent heterogenous plasticity (dependent on the model, the property, and the deleted channel), even introducing changes in properties that were not directly connected to the deleted channel. Together, post-knockout plasticity geared towards maintaining homeostasis introduced heterogenous off-target effects on several channels and properties, suggesting that extreme caution be exercised in interpreting experimental outcomes involving channel knockouts.

Cristae formation is a mechanical buckling event controlled by the inner membrane lipidome

Cristae are high curvature structures in the inner mitochondrial membrane (IMM) that are crucial for ATP production. While cristae-shaping proteins have been defined, analogous mechanisms for lipids have yet to be elucidated. Here we combine experimental lipidome dissection with multi-scale modeling to investigate how lipid interactions dictate IMM morphology and ATP generation. When modulating phospholipid (PL) saturation in engineered yeast strains, we observed a surprisingly abrupt breakpoint in IMM topology driven by a continuous loss of ATP synthase organization at cristae ridges. We found that cardiolipin (CL) specifically buffers the IMM against curvature loss, an effect that is independent of ATP synthase dimerization. To explain this interaction, we developed a continuum model for cristae tubule formation that integrates both lipid and protein-mediated curvatures. The model highlighted a snapthrough instability, which drives IMM collapse upon small changes in membrane properties. We also showed that CL is essential in low oxygen conditions that promote PL saturation. These results demonstrate that the mechanical function of CL is dependent on the surrounding lipid and protein components of the IMM.

CellRemorph: A Toolkit for Transforming, Selecting, and Slicing 3D Cell Structures on the Road to Morphologically Detailed Astrocyte Simulations

Understanding functions of astrocytes can be greatly enhanced by building and simulating computational models that capture their morphological details. Novel computational tools enable utilization of existing morphological data of astrocytes and building models that have appropriate level of details for specific simulation purposes. In addition to analyzing existing computational tools for constructing, transforming, and assessing astrocyte morphologies, we present here the CellRemorph toolkit implemented as an add-on for Blender, a 3D modeling platform increasingly recognized for its utility for manipulating 3D biological data. To our knowledge, CellRemorph is the first toolkit for transforming astrocyte morphologies from polygonal surface meshes into adjustable surface point clouds and vice versa, precisely selecting nanoprocesses, and slicing morphologies into segments with equal surface areas or volumes. CellRemorph is an open-source toolkit under the GNU General Public License and easily accessible via an intuitive graphical user interface. CellRemorph will be a valuable addition to other Blender add-ons, providing novel functionality that facilitates the creation of realistic astrocyte morphologies for different types of morphologically detailed simulations elucidating the role of astrocytes both in health and disease.

Surface-guided computing to analyze subcellular morphology and membrane-associated signals in 3D

Signal transduction and cell function are governed by the spatiotemporal organization of membrane-associated molecules. Despite significant advances in visualizing molecular distributions by 3D light microscopy, cell biologists still have limited quantitative understanding of the processes implicated in the regulation of molecular signals at the whole cell scale. In particular, complex and transient cell surface morphologies challenge the complete sampling of cell geometry, membrane-associated molecular concentration and activity and the computing of meaningful parameters such as the cofluctuation between morphology and signals. Here, we introduce u-Unwrap3D, a framework to remap arbitrarily complex 3D cell surfaces and membrane-associated signals into equivalent lower dimensional representations. The mappings are bidirectional, allowing the application of image processing operations in the data representation best suited for the task and to subsequently present the results in any of the other representations, including the original 3D cell surface. Leveraging this surface-guided computing paradigm, we track segmented surface motifs in 2D to quantify the recruitment of Septin polymers by blebbing events; we quantify actin enrichment in peripheral ruffles; and we measure the speed of ruffle movement along topographically complex cell surfaces. Thus, u-Unwrap3D provides access to spatiotemporal analyses of cell biological parameters on unconstrained 3D surface geometries and signals.

Surface-guided computing to analyze subcellular morphology and membrane-associated signals in 3D

Signal transduction and cell function are governed by the spatiotemporal organization of membrane-associated molecules. Despite significant advances in visualizing molecular distributions by 3D light microscopy, cell biologists still have limited quantitative understanding of the processes implicated in the regulation of molecular signals at the whole cell scale. In particular, complex and transient cell surface morphologies challenge the complete sampling of cell geometry, membrane-associated molecular concentration and activity and the computing of meaningful parameters such as the cofluctuation between morphology and signals. Here, we introduce u-Unwrap3D, a framework to remap arbitrarily complex 3D cell surfaces and membrane-associated signals into equivalent lower dimensional representations. The mappings are bidirectional, allowing the application of image processing operations in the data representation best suited for the task and to subsequently present the results in any of the other representations, including the original 3D cell surface. Leveraging this surface-guided computing paradigm, we track segmented surface motifs in 2D to quantify the recruitment of Septin polymers by blebbing events; we quantify actin enrichment in peripheral ruffles; and we measure the speed of ruffle movement along topographically complex cell surfaces. Thus, u-Unwrap3D provides access to spatiotemporal analyses of cell biological parameters on unconstrained 3D surface geometries and signals.

Understanding the interplay of membrane trafficking, cell surface mechanics, and stem cell differentiation

Stem cells can generate a diversity of cell types during development, regeneration and adult tissue homeostasis. Differentiation changes not only the cell fate in terms of gene expression but also the physical properties and functions of cells, e.g. the secretory activity, cell shape, or mechanics. Conversely, these activities and properties can also regulate differentiation itself. Membrane trafficking is known to modulate signal transduction and thus has the potential to control stem cell differentiation. On the other hand, membrane trafficking, particularly from and to the plasma membrane, depends on the mechanical properties of the cell surface such as tension within the plasma membrane or the cortex. Indeed, recent findings demonstrate that cell surface mechanics can also control cell fate. Here, we review the bidirectional relationships between these three fundamental cellular functions, i.e. membrane trafficking, cell surface mechanics, and stem cell differentiation. Furthermore, we discuss commonly used methods in each field and how combining them with new tools will enhance our understanding of their interplay. Understanding how membrane trafficking and cell surface mechanics can guide stem cell fate holds great potential as these concepts could be exploited for directed differentiation of stem cells for the fields of tissue engineering and regenerative medicine.

Measuring dimensionality of cell-scaffold contacts of primary human bone marrow stromal cells cultured on electrospun fiber scaffolds

The properties and structure of the cellular microenvironment can influence cell behavior. Sites of cell adhesion to the extracellular matrix (ECM) initiate intracellular signaling that directs cell functions such as proliferation, differentiation, and apoptosis. Electrospun fibers mimic the fibrous nature of native ECM proteins and cell culture in fibers affects cell shape and dimensionality, which can drive specific functions, such as the osteogenic differentiation of primary human bone marrow stromal cells (hBMSCs), by. In order to probe how scaffolds affect cell shape and behavior, cell-fiber contacts were imaged to assess their shape and dimensionality through a novel approach. Fluorescent polymeric fiber scaffolds were made so that they could be imaged by confocal fluorescence microscopy. Fluorescent polymer films were made as a planar control. hBSMCs were cultured on the fluorescent substrates and the cells and substrates were imaged. Two different image analysis approaches, one having geometrical assumptions and the other having statistical assumptions, were used to analyze the 3D structure of cell-scaffold contacts. The cells cultured in scaffolds contacted the fibers in multiple planes over the surface of the cell, while the cells cultured on films had contacts confined to the bottom surface of the cell. Shape metric analysis indicated that cell-fiber contacts had greater dimensionality and greater 3D character than the cell-film contacts. These results suggest that cell adhesion site-initiated signaling could emanate from multiple planes over the cell surface during culture in fibers, as opposed to emanating only from the cell's basal surface during culture on planar surfaces.

An Approach for Systems-Level Understanding of Prostate Cancer from High-Throughput Data Integration to Pathway Modeling and Simulation

To understand complex diseases, high-throughput data are generated at large and multiple levels. However, extracting meaningful information from large datasets for comprehensive understanding of cell phenotypes and disease pathophysiology remains a major challenge. Despite tremendous advances in understanding molecular mechanisms of cancer and its progression, current knowledge appears discrete and fragmented. In order to render this wealth of data more integrated and thus informative, we have developed a GECIP toolbox to investigate the crosstalk and the responsible genes'/proteins' connectivity of enriched pathways from gene expression data. To implement this toolbox, we used mainly gene expression datasets of prostate cancer, and the three datasets were GSE17951, GSE8218, and GSE1431. The raw samples were processed for normalization, prediction of differentially expressed genes, and the prediction of enriched pathways for the differentially expressed genes. The enriched pathways have been processed for crosstalk degree calculations for which number connections per gene, the frequency of genes in the pathways, sharing frequency, and the connectivity have been used. For network prediction, protein-protein interaction network database FunCoup2.0 was used, and cytoscape software was used for the network visualization. In our results, we found that there were enriched pathways 27, 45, and 22 for GSE17951, GSE8218, and GSE1431, respectively, and 11 pathways in common between all of them. From the crosstalk results, we observe that focal adhesion and PI3K pathways, both experimentally proven central for cellular output upon perturbation of numerous individual/distinct signaling pathways, displayed highest crosstalk degree. Moreover, we also observe that there were more critical pathways which appear to be highly significant, and these pathways are HIF1a, hippo, AMPK, and Ras. In terms of the pathways' components, GSK3B, YWHAE, HIF1A, ATP1A3, and PRKCA are shared between the aforementioned pathways and have higher connectivity with the pathways and the other pathway components. Finally, we conclude that the focal adhesion and PI3K pathways are the most critical pathways, and since for many other pathways, high-rank enrichment did not translate to high crosstalk degree, the global impact of one pathway on others appears distinct from enrichment.

Frequency-dependent signaling in cardiac myocytes

Background: Recent experimental data support the view that signaling activity at the membrane depends on its geometric parameters such as surface area and curvature. However, a mathematical, biophysical concept linking shape to receptor signaling is missing. The membranes of cardiomyocytes are constantly reshaped due to cycles of contraction and relaxation. According to constant-volume behavior of cardiomyocyte contraction, the length shortening is compensated by Z-disc myofilament lattice expansion and dynamic deformation of membrane between two adjacent Z-discs. Both morphological changes are strongly dependent on the frequency of contraction. Here, we developed the hypothesis that dynamic geometry of cardiomyocytes could be important for their plasticity and signaling. This effect may depend on the frequency of the beating heart and may represent a novel concept to explain how changes in frequency affect cardiac signaling.

Methods: This hypothesis is almost impossible to answer with experiments, as the in-vitro cardiomyocytes are almost two-dimensional and flattened rather than being in their real in-vivo shape. Therefore, we designed a COMSOL multiphysics program to mathematically model the dynamic geometry of a human cardiomyocyte and explore whether the beating frequency can modulate membrane signal transduction. Src kinase is an important component of cardiac mechanotransduction. We first presented that Src mainly localizes at costameres. Then, the frequency-dependent signaling effect was studied mathematically by numerical simulation of Src-mediated PDGFR signaling pathway. The reaction-convection-diffusion partial differential equation was formulated to simulate PDGFR pathway in a contracting sarcomeric disc for a range of frequencies from 1 to 4 Hz. Results: Simulations exhibits higher concentration of phospho-Src when a cardiomyocyte beats with higher rates. The calculated phospho-Src concentration at 4, 2, and 1 Hz beat rates, comparing to 0 Hz, was 21.5%, 9.4%, and 4.7% higher, respectively.

Conclusion: Here we provide mathematical evidence for a novel concept in biology. Cell shape directly translates into signaling, an effect of importance particularly for the myocardium, where cells continuously reshape their membranes. The concept of locality of surface-to-volume ratios is demonstrated to lead to changes in membrane-mediated signaling and may help to explain the remarkable plasticity of the myocardium in response to biomechanical stress.

Cell Architecture-Dependent Constraints: Critical Safeguards to Carcinogenesis

Animal cells display great diversity in their shape. These morphological characteristics result from crosstalk between the plasma membrane and the force-generating capacities of the cytoskeleton macromolecules. Changes in cell shape are not merely byproducts of cell fate determinants, they also actively drive cell fate decisions, including proliferation and differentiation. Global and local changes in cell shape alter the transcriptional program by a multitude of mechanisms, including the regulation of physical links between the plasma membrane and the nuclear envelope and the mechanical modulation of cation channels and signalling molecules. It is therefore not surprising that anomalies in cell shape contribute to several diseases, including cancer. In this review, we discuss the possibility that the constraints imposed by cell shape determine the behaviour of normal and pro-tumour cells by organizing the whole interconnected regulatory network. In turn, cell behaviour might stabilize cells into discrete shapes. However, to progress towards a fully transformed phenotype and to acquire plasticity properties, pro-tumour cells might need to escape these cell shape restrictions. Thus, robust controls of the cell shape machinery may represent a critical safeguard against carcinogenesis.

Sensing the shape of a cell with reaction diffusion and energy minimization

Some dividing cells sense their shape by becoming polarized along their long axis. Cell polarity is controlled in part by polarity proteins, like Rho GTPases, cycling between active membrane-bound forms and inactive cytosolic forms, modeled as a "wave-pinning" reaction-diffusion process. Does shape sensing emerge from wave pinning? We show that wave pinning senses the cell's long axis. Simulating wave pinning on a curved surface, we find that high-activity domains migrate to peaks and troughs of the surface. For smooth surfaces, a simple rule of minimizing the domain perimeter while keeping its area fixed predicts the final position of the domain and its shape. However, when we introduce roughness to our surfaces, shape sensing can be disrupted, and high-activity domains can become localized to locations other than the global peaks and valleys of the surface. On rough surfaces, the domains of the wave-pinning model are more robust in finding the peaks and troughs than the minimization rule, although both can become trapped in steady states away from the peaks and valleys. We can control the robustness of shape sensing by altering the Rho GTPase diffusivity and the domain size. We also find that the shape-sensing properties of cell polarity models can explain how domains localize to curved regions of deformed cells. Our results help to understand the factors that allow cells to sense their shape-and the limits that membrane roughness can place on this process.

Dendritic spine morphology regulates calcium-dependent synaptic weight change

Dendritic spines act as biochemical computational units and must adapt their responses according to their activation history. Calcium influx acts as the first signaling step during postsynaptic activation and is a determinant of synaptic weight change. Dendritic spines also come in a variety of sizes and shapes. To probe the relationship between calcium dynamics and spine morphology, we used a stochastic reaction-diffusion model of calcium dynamics in idealized and realistic geometries. We show that despite the stochastic nature of the various calcium channels, receptors, and pumps, spine size and shape can modulate calcium dynamics and subsequently synaptic weight updates in a deterministic manner. Through a series of exhaustive simulations and analyses, we found that the calcium dynamics and synaptic weight change depend on the volume-to-surface area of the spine. The relationships between calcium dynamics and spine morphology identified in idealized geometries also hold in realistic geometries, suggesting that there are geometrically determined deterministic relationships that may modulate synaptic weight change.

Localization of signaling receptors maximizes cellular information acquisition in spatially structured natural environments

Cells in natural environments, such as tissue or soil, sense and respond to extracellular ligands with intricately structured and non-monotonic spatial distributions, sculpted by processes such as fluid flow and substrate adhesion. In this work, we show that spatial sensing and navigation can be optimized by adapting the spatial organization of signaling pathways to the spatial structure of the environment. We develop an information-theoretic framework for computing the optimal spatial organization of a sensing system for a given signaling environment. We find that receptor localization previously observed in cells maximizes information acquisition in simulated natural contexts, including tissue and soil. Specifically, information acquisition is maximized when receptors form localized patches at regions of maximal ligand concentration. Receptor localization extends naturally to produce a dynamic protocol for continuously redistributing signaling receptors, which when implemented using simple feedback, boosts cell navigation efficiency by 30-fold.

Differences in cell death and division rules can alter tissue rigidity and fluidization

Tissue mechanical properties such as rigidity and fluidity, and changes in these properties driven by jamming-unjamming transitions (UJT), have come under recent highlight as mechanical markers of health and disease in various biological processes including cancer. However, most analyses of these mechanical properties and UJT have sidestepped the effect of cellular death and division in these systems. Cellular apoptosis (programmed cell death) and mitosis (cell division) can drive significant changes in tissue properties. The balance between the two is crucial in maintaining tissue function, and an imbalance between the two is seen in situations such as cancer progression, wound healing and necrosis. In this work we investigate the impact of cell death and division on tissue mechanical properties, by incorporating specific mechanosensitive triggers of cell death and division based on the size and geometry of the cell within in silico models of tissue dynamics. Specifically, we look at cell migration, tissue response to external stress, tissue extrusion propensity and self-organization of different cell types within the tissue, as a function of cell death and division and the rules that trigger these events. We find that not only do cell death and division events significantly alter tissue mechanics when compared to systems without these events, but that the choice of triggers driving these cell death and division events also alters the predicted tissue mechanics and overall system behavior.

The Cell Physiome: What Do We Need in a Computational Physiology Framework for Predicting Single-Cell Biology?

Modern biology and biomedicine are undergoing a big data explosion, needing advanced computational algorithms to extract mechanistic insights on the physiological state of living cells. We present the motivation for the Cell Physiome project: a framework and approach for creating, sharing, and using biophysics-based computational models of single-cell physiology. Using examples in calcium signaling, bioenergetics, and endosomal trafficking, we highlight the need for spatially detailed, biophysics-based computational models to uncover new mechanisms underlying cell biology. We review progress and challenges to date toward creating cell physiome models. We then introduce bond graphs as an efficient way to create cell physiome models that integrate chemical, mechanical, electromagnetic, and thermal processes while maintaining mass and energy balance. Bond graphs enhance modularization and reusability of computational models of cells at scale. We conclude with a look forward at steps that will help fully realize this exciting new field of mechanistic biomedical data science. Expected final online publication date for the Annual Review of Biomedical Data Science, Volume 5 is August 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Deformation of the nucleus by TGFβ1 via the remodeling of nuclear envelope and histone isoforms

The cause of nuclear shape abnormalities which are often seen in pre-neoplastic and malignant tissues is not clear. In this study we report that deformation of the nucleus can be induced by TGFβ1 stimulation in several cell lines including Huh7. In our results, the upregulated histone H3.3 expression downstream of SMAD signaling contributed to TGFβ1-induced nuclear deformation, a process of which requires incorporation of the nuclear envelope (NE) proteins lamin B1 and SUN1. During this process, the NE constitutively ruptured and reformed. Contrast to lamin B1 which was relatively stationary around the nucleus, the upregulated lamin A was highly mobile, clustering at the nuclear periphery and reintegrating into the nucleoplasm. The chromatin regions that lost NE coverage formed a supra-nucleosomal structure characterized by elevated histone H3K27me3 and histone H1, the formation of which depended on the presence of lamin A. These results provide evidence that shape of the nucleus can be modulated through TGFβ1-induced compositional changes in the chromatin and nuclear lamina.

Three-Dimensional Visualization of the Podocyte Actin Network Using Integrated Membrane Extraction, Electron Microscopy, and Machine Learning

BACKGROUND

Actin stress fibers are abundant in cultured cells, but little is known about them in vivo. In podocytes, much evidence suggests that mechanobiologic mechanisms underlie podocyte shape and adhesion in health and in injury, with structural changes to actin stress fibers potentially responsible for pathologic changes to cell morphology. However, this hypothesis is difficult to rigorously test in vivo due to challenges with visualization. A technology to image the actin cytoskeleton at high resolution is needed to better understand the role of structures such as actin stress fibers in podocytes.

METHODS

We developed the first visualization technique capable of resolving the three-dimensional cytoskeletal network in mouse podocytes in detail, while definitively identifying the proteins that comprise this network. This technique integrates membrane extraction, focused ion-beam scanning electron microscopy, and machine learning image segmentation.

RESULTS

Using isolated mouse glomeruli from healthy animals, we observed actin cables and intermediate filaments linking the interdigitated podocyte foot processes to newly described contractile actin structures, located at the periphery of the podocyte cell body. Actin cables within foot processes formed a continuous, mesh-like, electron-dense sheet that incorporated the slit diaphragms.

CONCLUSIONS

Our new technique revealed, for the first time, the detailed three-dimensional organization of actin networks in healthy podocytes. In addition to being consistent with the gel compression hypothesis, which posits that foot processes connected by slit diaphragms act together to counterbalance the hydrodynamic forces across the glomerular filtration barrier, our data provide insight into how podocytes respond to mechanical cues from their surrounding environment.

Stochasticity and positive feedback enable enzyme kinetics at the membrane to sense reaction size

Cellular membranes span a wide range of spatial dimensions, from the plasma membrane with a scale of microns to vesicles on the nanometer scale. The work presented here identifies a molecular mechanism, based on common features of cellular signaling enzymes, that causes the average enzymatic catalytic rate to exhibit reaction size dependency. This effect stems from stochastic variation, but the final results can be essentially deterministic. In competitive enzymatic reaction cycles, the final product can depend on the size of the reaction system. The simplicity of the mechanism suggests that size-dependent reaction rates may be widespread among signaling enzymes and thus enable reaction size to be an important factor in signal regulation at the membrane.

Here, we present detailed kinetic analyses of a panel of soluble lipid kinases and phosphatases, as well as Ras activating proteins, acting on their respective membrane surface substrates. The results reveal that the mean catalytic rate of such interfacial enzymes can exhibit a strong dependence on the size of the reaction system—in this case membrane area. Experimental measurements and kinetic modeling reveal how stochastic effects stemming from low molecular copy numbers of the enzymes alter reaction kinetics based on mechanistic characteristics of the enzyme, such as positive feedback. For the competitive enzymatic cycles studied here, the final product—consisting of a specific lipid composition or Ras activity state—depends on the size of the reaction system. Furthermore, we demonstrate how these reaction size dependencies can be controlled by engineering feedback mechanisms into the enzymes.

On the Three-Dimensional Correlation Between Myofibroblast Shape and Contraction

Myofibroblasts are responsible for wound healing and tissue repair across all organ systems. In periods of growth and disease, myofibroblasts can undergo a phenotypic transition characterized by an increase in extracellular matrix (ECM) deposition rate, changes in various protein expression (e.g., alpha-smooth muscle actin (αSMA)), and elevated contractility. Cell shape is known to correlate closely with stress-fiber geometry and function and is thus a critical feature of cell biophysical state. However, the relationship between myofibroblast shape and contraction is complex, even as well in regards to steady-state contractile level (basal tonus). At present, the relationship between myofibroblast shape and basal tonus in three-dimensional (3D) environments is poorly understood. Herein, we utilize the aortic valve interstitial cell (AVIC) as a representative myofibroblast to investigate the relationship between basal tonus and overall cell shape. AVICs were embedded within 3D poly(ethylene glycol) (PEG) hydrogels containing degradable peptide crosslinkers, adhesive peptide sequences, and submicron fluorescent microspheres to track the local displacement field. We then developed a methodology to evaluate the correlation between overall AVIC shape and basal tonus induced contraction. We computed a volume averaged stretch tensor ⟨U⟩ for the volume occupied by the AVIC, which had three distinct eigenvalues (λ1,2,3=1.08,0.99, and 0.89), suggesting that AVIC shape is a result of anisotropic contraction. Furthermore, the direction of maximum contraction correlated closely with the longest axis of a bounding ellipsoid enclosing the AVIC. As gel-imbedded AVICs are known to be in a stable state by 3 days of incubation used herein, this finding suggests that the overall quiescent AVIC shape is driven by the underlying stress-fiber directional structure and potentially contraction level.

Describing the movement of molecules in reduced-dimension models

When addressing spatial biological questions using mathematical models, symmetries within the system are often exploited to simplify the problem by reducing its physical dimension. In a reduced-dimension model molecular movement is restricted to the reduced dimension, changing the nature of molecular movement. This change in molecular movement can lead to quantitatively and even qualitatively different results in the full and reduced systems. Within this manuscript we discuss the condition under which restricted molecular movement in reduced-dimension models accurately approximates molecular movement in the full system. For those systems which do not satisfy the condition, we present a general method for approximating unrestricted molecular movement in reduced-dimension models. We will derive a mathematically robust, finite difference method for solving the 2D diffusion equation within a 1D reduced-dimension model. The methods described here can be used to improve the accuracy of many reduced-dimension models while retaining benefits of system simplification.

Systems modeling predicts that mitochondria ER contact sites regulate the postsynaptic energy landscape

Spatiotemporal compartmentation of calcium dynamics is critical for neuronal function, particularly in postsynaptic spines. This exquisite level of Ca²⁺ compartmentalization is achieved through the storage and release of Ca²⁺ from various intracellular organelles particularly the endoplasmic reticulum (ER) and the mitochondria. Mitochondria and ER are established storage organelles controlling Ca²⁺ dynamics in neurons. Mitochondria also generate a majority of energy used within postsynaptic spines to support the downstream events associated with neuronal stimulus. Recently, high resolution microscopy has unveiled direct contact sites between the ER and the mitochondria (MERCs), which directly channel Ca²⁺ release from the ER into the mitochondrial membrane. In this study, we develop a computational 3D reaction-diffusion model to investigate the role of MERCs in regulating Ca²⁺ and ATP dynamics. This spatiotemporal model accounts for Ca²⁺ oscillations initiated by glutamate stimulus of metabotropic and ionotropic glutamate receptors and Ca²⁺ changes in four different compartments: cytosol, ER, mitochondria, and the MERC microdomain. Our simulations predict that the organization of these organelles and inter-organellar contact sites play a key role in modulating Ca²⁺ and ATP dynamics.

We further show that the crosstalk between geometry (mitochondria and MERC) and metabolic parameters (cytosolic ATP hydrolysis, ATP generation) influences the neuronal energy state. Our findings shed light on the importance of organelle interactions in predicting Ca²⁺ dynamics in synaptic signaling. Overall, our model predicts that a combination of MERC linkage and mitochondria size is necessary for optimal ATP production in the cytosol.

Wnt- and glutamate-receptors orchestrate stem cell dynamics and asymmetric cell division

The Wnt-pathway is part of a signalling network that regulates many aspects of cell biology. Recently, we discovered crosstalk between AMPA/Kainate-type ionotropic glutamate receptors (iGluRs) and the Wnt-pathway during the initial Wnt3a-interaction at the cytonemes of mouse embryonic stem cells (ESCs). Here, we demonstrate that this crosstalk persists throughout the Wnt3a-response in ESCs. Both AMPA and Kainate receptors regulate early Wnt3a-recruitment, dynamics on the cell membrane, and orientation of the spindle towards a Wnt3a-source at mitosis. AMPA receptors specifically are required for segregating cell fate components during Wnt3a-mediated asymmetric cell division (ACD). Using Wnt-pathway component knockout lines, we determine that Wnt co-receptor Lrp6 has particular functionality over Lrp5 in cytoneme formation, and in facilitating ACD. Both Lrp5 and 6, alongside pathway effector β-catenin act in concert to mediate the positioning of the dynamic interaction with, and spindle orientation to, a localised Wnt3a-source. Wnt-iGluR crosstalk may prove pervasive throughout embryonic and adult stem cell signalling.

A spatial model of YAP/TAZ signaling reveals how stiffness, dimensionality, and shape contribute to emergent outcomes

YAP/TAZ is a master regulator of mechanotransduction whose functions rely on translocation from the cytoplasm to the nucleus in response to diverse physical cues. Substrate stiffness, substrate dimensionality, and cell shape are all input signals for YAP/TAZ, and through this pathway, regulate critical cellular functions and tissue homeostasis. Yet, the relative contributions of each biophysical signal and the mechanisms by which they synergistically regulate YAP/TAZ in realistic tissue microenvironments that provide multiplexed input signals remain unclear. For example, in simple two-dimensional culture, YAP/TAZ nuclear localization correlates strongly with substrate stiffness, while in three-dimensional (3D) environments, YAP/TAZ translocation can increase with stiffness, decrease with stiffness, or remain unchanged. Here, we develop a spatial model of YAP/TAZ translocation to enable quantitative analysis of the relationships between substrate stiffness, substrate dimensionality, and cell shape. Our model couples cytosolic stiffness to nuclear mechanics to replicate existing experimental trends, and extends beyond current data to predict that increasing substrate activation area through changes in culture dimensionality, while conserving cell volume, forces distinct shape changes that result in nonlinear effect on YAP/TAZ nuclear localization. Moreover, differences in substrate activation area versus total membrane area can account for counterintuitive trends in YAP/TAZ nuclear localization in 3D culture. Based on this multiscale investigation of the different system features of YAP/TAZ nuclear translocation, we predict that how a cell reads its environment is a complex information transfer function of multiple mechanical and biochemical factors. These predictions reveal a few design principles of cellular and tissue engineering for YAP/TAZ mechanotransduction.

Ionic Liquid-Controlled Shape Transformation of Spherical to Nonspherical Polymersomes via Hierarchical Self-Assembly of a Diblock Copolymer

Here, we report the self-assembly of poly(ethylene glycol) methyl ether-block-poly(ε-caprolactone) (PEG-b-PCL) copolymer in three ionic liquids (ILs) possessing different cations with common bis(trifluoromethylsulfonyl)imide anion. The observed polymeric nanostructures in ILs were directly visualized by room temperature conventional transmission and field emission scanning electron microscopy and were further examined for their size and shape by dynamic light scattering technique. The results show that through changes in the concentration of PEG-b-PCL and/or changing the solvent by using a different IL, we can effectively induce shape transformation of self-assembled PEG-b-PCL nanostructures in order to generate nonspherical polymersomes, such as worm-like aggregates, stomatocytes, nanotubes, large hexagonal and tubular-shaped polymersomes. These findings provide a promising platform for the design of biodegradable soft dynamic systems in the micro-/nano-motor field for cancer-targeted delivery, diagnosis and imaging-guided therapy, and controlled release of therapeutic drugs for treatment of many diseases. Non-spherical polymersome-based vaccines may be taken up more efficiently, especially against viruses for pulmonary drug delivery than the spherical polymersomes-based.

Active perception during angiogenesis: filopodia speed up Notch selection of tip cells in silico and in vivo

How do cells make efficient collective decisions during tissue morphogenesis? Humans and other organisms use feedback between movement and sensing known as 'sensorimotor coordination' or 'active perception' to inform behaviour, but active perception has not before been investigated at a cellular level within organs. Here we provide the first proof of concept in silico/in vivo study demonstrating that filopodia (actin-rich, dynamic, finger-like cell membrane protrusions) play an unexpected role in speeding up collective endothelial decisions during the time-constrained process of 'tip cell' selection during blood vessel formation (angiogenesis). We first validate simulation predictions in vivo with live imaging of zebrafish intersegmental vessel growth. Further simulation studies then indicate the effect is due to the coupled positive feedback between movement and sensing on filopodia conferring a bistable switch-like property to Notch lateral inhibition, ensuring tip selection is a rapid and robust process. We then employ measures from computational neuroscience to assess whether filopodia function as a primitive (basal) form of active perception and find evidence in support. By viewing cell behaviour through the 'basal cognitive lens' we acquire a fresh perspective on the tip cell selection process, revealing a hidden, yet vital time-keeping role for filopodia. Finally, we discuss a myriad of new and exciting research directions stemming from our conceptual approach to interpreting cell behaviour. This article is part of the theme issue 'Basal cognition: multicellularity, neurons and the cognitive lens'.

A self-organized synthetic morphogenic liposome responds with shape changes to local light cues

Reconstituting artificial proto-cells capable of transducing extracellular signals into cytoskeletal changes can reveal fundamental principles of how non-equilibrium phenomena in cellular signal transduction affect morphogenesis. Here, we generated a Synthetic Morphogenic Membrane System (SynMMS) by encapsulating a dynamic microtubule (MT) aster and a light-inducible signaling system driven by GTP/ATP chemical potential into cell-sized liposomes. Responding to light cues in analogy to morphogens, this biomimetic design embodies basic principles of localized Rho-GTPase signal transduction that generate an intracellular MT-regulator signaling gradient. Light-induced signaling promotes membrane-deforming growth of MT-filaments by dynamically elevating the membrane-proximal tubulin concentration. The resulting membrane deformations enable recursive coupling of the MT-aster with the signaling system, which generates global self-organized morphologies that reorganize towards local external cues in dependence on prior shape. SynMMS thereby signifies a step towards bio-inspired engineering of self-organized cellular morphogenesis.

The authors generated a Synthetic Morphogenic Membrane System by encapsulating a dynamic microtubule aster and a light-inducible signaling system driven by GTP/ATP chemical potential into cell-sized liposomes. This reconstitution of artificial proto-cells reveals how non-equilibrium phenomena affect cellular information processing in morphogenesis.

Local and Systemic In Vivo Responses to Osseointegrative Titanium Nanotube Surfaces

Orthopedic implants requiring osseointegration are often surface modified; however, implants may shed these coatings and generate wear debris leading to complications. Titanium nanotubes (TiNT), a new surface treatment, may promote osseointegration. In this study, in vitro (rat marrow-derived bone marrow cell attachment and morphology) and in vivo (rat model of intramedullary fixation) experiments characterized local and systemic responses of two TiNT surface morphologies, aligned and trabecular, via animal and remote organ weight, metal ion, hematologic, and nondecalcified histologic analyses. In vitro experiments showed total adherent cells on trabecular and aligned TiNT surfaces were greater than control at 30 min and 4 h, and cells were smaller in diameter and more eccentric. Control animals gained more weight, on average; however, no animals met the institutional trigger for weight loss. No hematologic parameters (complete blood count with differential) were significantly different for TiNT groups vs. control. Inductively coupled plasma mass spectrometry (ICP-MS) showed greater aluminum levels in the lungs of the trabecular TiNT group than in those of the controls. Histologic analysis demonstrated no inflammatory infiltrate, cytotoxic, or necrotic conditions in proximity of K-wires. There were significantly fewer eosinophils/basophils and neutrophils in the distal region of trabecular TiNT-implanted femora; and, in the midshaft of aligned TiNT-implanted femora, there were significantly fewer foreign body giant/multinucleated cells and neutrophils, indicating a decreased immune response in aligned TiNT-implanted femora compared to controls.

The nuclear import of the transcription factor MyoD is reduced in mesenchymal stem cells grown in a 3D micro-engineered niche

Smart biomaterials are increasingly being used to control stem cell fate in vitro by the recapitulation of the native niche microenvironment. By integrating experimental measurements with numerical models, we show that in mesenchymal stem cells grown inside a 3D synthetic niche both nuclear transport of a myogenic factor and the passive nuclear diffusion of a smaller inert protein are reduced. Our results also suggest that cell morphology modulates nuclear proteins import through a partition of the nuclear envelope surface, which is a thin but extremely permeable annular portion in cells cultured on 2D substrates. Therefore, our results support the hypothesis that in stem cell differentiation, the nuclear import of gene-regulating transcription factors is controlled by a strain-dependent nuclear envelope permeability, probably related to the reorganization of stretch-activated nuclear pore complexes.

Machine-Learning-Based Approach to Decode the Influence of Nanomaterial Properties on Their Interaction with Cells

In an in vitro nanotoxicity system, cell-nanoparticle (NP) interaction leads to the surface adsorption, uptake, and changes into nuclei/cell phenotype and chemistry, as an indicator of oxidative stress, genotoxicity, and carcinogenicity. Different types of nanomaterials and their chemical composition or "corona" have been widely studied in context with nanotoxicology. However, rare reports are available, which delineate the details of the cell shape index (CSI) and nuclear area factors (NAFs) as a descriptor of the type of nanomaterials. In this paper, we propose a machine-learning-based graph modeling and correlation-establishing approach using tight junction protein ZO-1-mediated alteration in the cell/nuclei phenotype to quantify and propose it as indices of cell-NP interactions. We believe that the phenotypic variation (CSI and NAF) in the epithelial cell is governed by the physicochemical descriptors (e.g., shape, size, zeta potential, concentration, diffusion coefficients, polydispersity, and so on) of the different classes of nanomaterials, which critically determines the intracellular uptake or cell membrane interactions when exposed to the epithelial cells at sub-lethal concentrations. The intrinsic and extrinsic physicochemical properties of the representative nanomaterials (NMs) were measured using optical (dynamic light scattering, NP tracking analysis) methods to create a set of nanodescriptors contributing to cell-NM interactions via phenotype adjustments. We used correlation function as a machine-learning algorithm to successfully predict cell and nuclei shapes and polarity functions as phenotypic markers for five different classes of nanomaterials studied herein this report. The CSI and NAF as nanodescriptors can be used as intuitive cell phenotypic parameters to define the safety of nanomaterials extensively used in consumer products and nanomedicine.

Biophysical phenotyping of mesenchymal stem cells along the osteogenic differentiation pathway

Mesenchymal stem cells represent an important resource, for bone regenerative medicine and therapeutic applications. This review focuses on new advancements and biophysical tools which exploit different physical and chemical markers of mesenchymal stem cell populations, to finely characterize phenotype changes along their osteogenic differentiation process. Special attention is paid to recently developed label-free methods, which allow monitoring cell populations with minimal invasiveness. Among them, quantitative phase imaging, suitable for single-cell morphometric analysis, and nanoindentation, functional to cellular biomechanics investigation. Moreover, the pool of ion channels expressed in cells during differentiation is discussed, with particular interest for calcium homoeostasis.Altogether, a biophysical perspective of osteogenesis is proposed, offering a valuable tool for the assessment of the cell stage, but also suggesting potential physiological links between apparently independent phenomena.

Tools for computational analysis of moving boundary problems in cellular mechanobiology

A cell's ability to change shape is one of the most fundamental biological processes and is essential for maintaining healthy organisms. When the ability to control shape goes awry, it often results in a diseased system. As such, it is important to understand the mechanisms that allow a cell to sense and respond to its environment so as to maintain cellular shape homeostasis. Because of the inherent complexity of the system, computational models that are based on sound theoretical understanding of the biochemistry and biomechanics and that use experimentally measured parameters are an essential tool. These models involve an inherent feedback, whereby shape is determined by the action of regulatory signals whose spatial distribution depends on the shape. To carry out computational simulations of these moving boundary problems requires special computational techniques. A variety of alternative approaches, depending on the type and scale of question being asked, have been used to simulate various biological processes, including cell motility, division, mechanosensation, and cell engulfment. In general, these models consider the forces that act on the system (both internally generated, or externally imposed) and the mechanical properties of the cell that resist these forces. Moving forward, making these techniques more accessible to the non-expert will help improve interdisciplinary research thereby providing new insight into important biological processes that affect human health. This article is categorized under: Cancer > Cancer>Computational Models Cancer > Cancer>Molecular and Cellular Physiology.

Development of a kinetic model expressing anomalous phenomena in human induced pluripotent stem cell culture

During culture with feeder cells, deviation from the undifferentiated state of human induced pluripotent stem cells (hiPSCs) occurs at a very low frequency. Anomalous cell migration in central and peripheral regions of hiPSC colonies has been suggested to be the trigger for this phenomenon. To confirm this hypothesis, sequential cell migration prior to deviation must be demonstrated. This has been difficult using in vitro methods. We therefore developed a kinetic model with a proposed definition of anomalous cell migration as continuous relatively fast or slow cell migration. The developed model was validated via in silico reproduction of deviation phenomenon observed in vitro, such as the positions of deviated cells in a colony and the frequency of deviation in culture. This model suggests that anomalous cell migration-driven hiPSC deviation can be explained by two factors: a mechanical stimulus, represented by cell migration, and duration of the mechanical stimulus. The factor "duration of mechanical stimulus" sets our model apart from others, and helps to realize the ultra-rare trigger (approximately 10^-5) of deviation from the undifferentiated state in hiPSC culture.

Cadmium-Induced Cell Homeostasis Impairment is Suppressed by the Tor1 Deficiency in Fission Yeast

Cadmium has no known physiological function in the body; however, its adverse effects are associated with cancer and many types of organ system damage. Although much has been shown about Cd toxicity, the underlying mechanisms of its responses to the organism remain unclear. In this study, the role of Tor1, a catalytic subunit of the target of rapamycin complex 2 (TORC2), in Cd-mediated effects on cell proliferation, the antioxidant system, morphology, and ionome balance was investigated in the eukaryotic model organism Schizosaccharomyces pombe. Surprisingly, spectrophotometric and biochemical analyses revealed that the growth rate conditions and antioxidant defense mechanisms are considerably better in cells lacking the Tor1 signaling. The malondialdehyde (MDA) content of Tor1-deficient cells upon Cd treatment represents approximately half of the wild-type content. The microscopic determination of the cell morphological parameters indicates the role for Tor1 in cell shape maintenance. The ion content, determined by inductively coupled plasma optical emission spectroscopy (ICP-OES), showed that the Cd uptake potency was markedly lower in Tor1-depleted compared to wild-type cells. Conclusively, we show that the cadmium-mediated cell impairments in the fission yeast significantly depend on the Tor1 signaling. Additionally, the data presented here suggest the yet-undefined role of Tor1 in the transport of ions.

Cell shape regulates subcellular organelle location to control early Ca2+ signal dynamics in vascular smooth muscle cells

The shape of the cell is connected to its function; however, we do not fully understand underlying mechanisms by which global shape regulates a cell's functional capabilities. Using theory, experiments and simulation, we investigated how physiologically relevant cell shape changes affect subcellular organization, and consequently intracellular signaling, to control information flow needed for phenotypic function. Vascular smooth muscle cells going from a proliferative and motile circular shape to a contractile fusiform shape show changes in the location of the sarcoplasmic reticulum, inter-organelle distances, and differential distribution of receptors in the plasma membrane. These factors together lead to the modulation of signals transduced by the M3 muscarinic receptor/Gq/PLCβ pathway at the plasma membrane, amplifying Ca2+ dynamics in the cytoplasm, and the nucleus resulting in phenotypic changes, as determined by increased activity of myosin light chain kinase in the cytoplasm and enhanced nuclear localization of the transcription factor NFAT. Taken together, our observations show a systems level phenomenon whereby global cell shape affects subcellular organization to modulate signaling that enables phenotypic changes.

Effect of cell imprinting on viability and drug susceptibility of breast cancer cells to doxorubicin

This study demonstrates the effect of substrate's geometrical cues on viability and the efficacy of an anti-cancer drug, doxorubicin (DOX), on breast cancer cells. It is hypothesized that the surface topographical properties can mediate the cellular drug intake. Pseudo-three dimensional (3D) platforms were fabricated using imprinting technique from polydimethylsiloxane (PDMS) and gelatin methacryloyl (GelMA) hydrogel to recapitulate topography of cells' membranes. The cells exhibited higher viability on the cell-imprinted platforms for both PDMS and GelMA materials compared to the plain/flat counterparts. For instance, MCF7 cells showed a higher metabolic activity (11.9%) on MCF7-imprinted PDMS substrate than plain PDMS. The increased metabolic activity for the imprinted GelMA was about 44.2% compared to plain hydrogel. The DOX response of cells was monitored for 24 h. Although imprinted substrates demonstrated enhanced biocompatibility, the cultured cells were more susceptible to the drug compared to the plain substrates. In particular, MCF7 cells on imprinted PDMS and GelMA substrates showed 37% and 50% higher in cell death compared to the corresponding plain PDMS and GelMA, respectively. Interestingly, the drug susceptibility of the cells on the imprinted hydrogel was about 70% higher than the cells cultured on imprinted PDMS substrates. Having MCF7 cell-imprinted substrates, DOX responses of two other breast cancer cell lines, SKBR3 and ZR-75-1, were also evaluated. The results support that cell membrane curvature developed by multiscale topography is able to mediate intracellular signaling and drug intake. STATEMENT OF SIGNIFICANCE: Research in biological sciences and drug discovery mostly rely on two dimensional (2D) cell culture techniques which cannot provide a reliable physiologically relevant environment. Lack of extracellular matrix and a large shift in physicochemical properties of conventional 2D substrates can induce aberrant cellular behaviors. While chemical composition, topographical, and mechanical properties of substrates have remarkable impacts on drug susceptibility, gene expression, and protein synthesis, the most cell culture plates are from rigid and plain substrates. A number of (bio)polymeric 3D-platforms have been introduced to resemble innate cell microenvironment. However, their intricate culture protocols restrain their applications in demanding high-throughput drug screening. To address the above concerns, in the present study, a hydrogel-based pseudo-3D substrate with imprinted cell features has been introduced.

In vitro characterization of osteoblast cells on polyelectrolyte multilayers containing detonation nanodiamonds

Nanoparticle-enhanced coatings of bone implants are a promising method to facilitate sustainable wound healing, leading to an increase in patient well-being. This article describes the in vitro characterization of osteoblast cells interacting with polyelectrolyte multilayers, which contain detonation nanodiamonds (NDs), as a novel class of carbon-based coating material, which presents a unique combination of photoluminescence and drug-binding properties. The cationic polyelectrolyte, namely polydiallyldimethylammonium chloride (PDDA), has been used to immobilize NDs on silica glass. The height of ND-PDDA multilayers varies from a minimum of 10 nm for one bilayer to a maximum of 90 nm for five bilayers of NDs and PDDA. Human fetal osteoblasts (hFOBs) cultured on ND-PDDA multilayers show a large number of focal adhesions, which were studied via quantitative fluorescence imaging analysis. The influence of the surface roughness on the filopodia formation was assessed via scanning electron microscopy and atomic force microscopy. The nano-rough surface of five bilayers constrained the filopodia formation. The hFOBs grown on NDs tend to show not only a similar cell morphology compared to cells cultured on extracellular matrix protein-coated silica glass substrates, but also increased cell viability by about 40%. The high biocompatibility of the ND-PDDA multilayers, indicated via high cell proliferation and sound cell adhesion, shows their potential for biomedical applications such as drug-eluting coatings and biomaterials in general.