Spatial patterning of different cell types is crucial for tissue engineering and is characterized by the formation of sharp boundary between segregated groups of cells of different lineages. The cell−cell boundary layers, depending on the relative adhesion forces, can result in kinks in the border, similar to fingering patterns between two viscous partially miscible fluids which can be characterized by its fractal dimension. This suggests that mathematical models used to analyze the fingering patterns can be applied to cell migration data as a metric for intercellular adhesion forces. In this study, we develop a novel computational analysis method to characterize the interactions between blood endothelial cells (BECs) and lymphatic endothelial cells (LECs), which form segregated vasculature by recognizing each other through podoplanin. We observed indiscriminate mixing with LEC−LEC and BEC−BEC pairs and a sharp boundary between LEC−BEC pair, and fingering-like patterns with pseudo-LEC−BEC pairs. We found that the box counting method yields fractal dimension between 1 for sharp boundaries and 1.3 for indiscriminate mixing, and intermediate values for fingering-like boundaries. We further verify that these results are due to differential affinity by performing random walk simulations with differential attraction to nearby cells and generate similar migration pattern, confirming that higher differential attraction between different cell types result in lower fractal dimensions. We estimate the characteristic velocity and interfacial tension for our simulated and experimental data to show that the fractal dimension negatively correlates with capillary number (
Zebrafish (
- Award ID(s):
- 1714429
- NSF-PAR ID:
- 10153582
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- Nature Communications
- Volume:
- 9
- Issue:
- 1
- ISSN:
- 2041-1723
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract Ca ), further indicating that the mathematical models used to study viscous fingering pattern can be used to characterize cell−cell mixing. Taken together, these results indicate that the fractal analysis of segregation boundaries can be used as a simple metric to estimate relative cell−cell adhesion forces between different cell types. -
more » « less
Abstract Proper cell–material interactions are critical to remain cell function and thus successful tissue regeneration. Many fabrication processes have been developed to create microenvironments to control cell attachment and organization on a three‐dimensional (3D) scaffold. However, these approaches often involve heavy engineering and only the surface layer can be patterned. We found that 3D extrusion based printing at high temperature and pressure will result an aligned effect on the polymer molecules, and this molecular arrangement will further induce the cell alignment and different differentiation capacities. In particular, articular cartilage tissue is known to have zonal collagen fiber and cell orientation to support different functions, where collagen fibers and chondrocytes align parallel, randomly, and perpendicular, respectively, to the surface of the joint. Therefore, cell alignment was evaluated in a cartilage model in this study. We used small angle X‐ray scattering analysis to substantiate the polymer molecule alignment phenomenon. The cellular response was evaluated both
in vitro andin vivo . Seeded mesenchymal stem cells (MSCs) showed different morphology and orientation on scaffolds, as a combined result of polymer molecule alignment and printed scaffold patterns. Gene expression results showed improved superficial zonal chondrogenic marker expression in parallel‐aligned group. The cell alignment was successfully maintained in the animal model after 7 days with distinct MSC morphology between the casted and parallel printed scaffolds. This 3D printing induced polymer and cell alignment will have a significant impact on developing scaffold with controlled cell–material interactions for complex tissue engineering while avoiding complicated surface treatment, and therefore provides new concept for effective tissue repairing in future clinical applications. © 2018 Wiley Periodicals, Inc. J Biomed Mater Res Part A: 106A: 2190‐2199, 2018. -
more » « less
Abstract HIV-1 viral transcription persists in patients despite antiretroviral treatment, potentially due to intermittent HIV-1 LTR activation. While several mathematical models have been explored in the context of LTR-protein interactions, in this work for the first time HIV-1 LTR model featuring repressed, intermediate, and activated LTR states is integrated with generation of long (
env ) and short (TAR) RNAs and proteins (Tat, Pr55, and p24) in T-cells and macrophages using both cell lines and infected primary cells. This type of extended modeling framework allows us to compare and contrast behavior of these two cell types. We demonstrate that they exhibit unique LTR dynamics, which ultimately results in differences in the magnitude of viral products generated. One of the distinctive features of this work is that it relies on experimental data in reaction rate computations. Two RNA transcription rates from the activated promoter states are fit by comparison of experimental data to model predictions. Fitting to the data also provides estimates for the degradation/exit rates for long and short viral RNA. Our experimentally generated data is in reasonable agreement for the T-cell as well macrophage population and gives strong evidence in support of using the proposed integrated modeling paradigm. Sensitivity analysis performed using Latin hypercube sampling method confirms robustness of the model with respect to small parameter perturbations. Finally, incorporation of a transcription inhibitor (F07#13) into the governing equations demonstrates how the model can be used to assess drug efficacy. Collectively, our model indicates transcriptional differences between latently HIV-1 infected T-cells and macrophages and provides a novel platform to study various transcriptional dynamics leading to latency or activation in numerous cell types and physiological conditions. -
more » « less
Abstract Mathematical and computational modeling enables biologists to integrate data from observations and experiments into a theoretical framework. In this review, we describe how developmental processes associated with stem‐cell‐driven growth of tissue in both the embryonic and adult nervous system can be modeled using cellular automata (CA). A cellular automaton is defined by its discrete nature in time, space, and state. The discrete space is represented by a uniform grid or lattice containing agents that interact with other agents within their local neighborhood. This possibility of local interactions of agents makes the cellular automata approach particularly well suited for studying through modeling how complex patterns at the tissue level emerge from fundamental developmental processes (such as proliferation, migration, differentiation, and death) at the single‐cell level. As part of this review, we provide a primer for how to define biologically inspired rules governing these processes so that they can be implemented into a CA model. We then demonstrate the power of the CA approach by presenting simulations (in the form of figures and movies) based on building models of three developmental systems: the formation of the enteric nervous system through invasion by neural crest cells; the growth of normal and tumorous neurospheres induced by proliferation of adult neural stem/progenitor cells; and the neural fate specification through lateral inhibition of embryonic stem cells in the neurogenic region of
Drosophila . -
more » « less
Abstract “Supermodeling” climate by allowing different models to assimilate data from one another in run time has been shown to give results superior to those of any one model and superior to any weighted average of model outputs. The only free parameters, connection strengths between corresponding variables in each pair of models, are determined using some form of machine learning. It is demonstrated that supermodeling succeeds because near critical states, interscale interactions are important but unresolved processes cannot be effectively represented diagnostically in any single parameterization scheme. In two examples, a pair of toy quasigeostrophic (QG) channel models of the midlatitudes and a pair of ECHAM5 models of the tropical Pacific atmosphere with a common ocean, supermodels dynamically combine parameterization schemes so as to capture criticality, associated critical structures, and the supporting scale interactions. The QG supermodeling scheme extends a previous configuration in which two such models synchronize with intermodel connections only between medium-scale components of the flow; here the connections are trained against a third “real” model. Intermittent blocking patterns characterize the critical behavior thus obtained, even where such patterns are missing in the constituent models. In the ECHAM-based climate supermodel, the corresponding critical structure is the single ITCZ pattern, a pattern that occurs in neither of the constituent models. For supermodels of both types, power spectra indicate enhanced interscale interactions in frequency or energy ranges of physical interest, in agreement with observed data, and supporting a generalized form of the self-organized criticality hypothesis.
Significance Statement In a “supermodel” of Earth’s climate, alternative models (climate simulations), which differ in the way they represent processes on the smallest scales, are trained to exchange information as they run, adjusting to one another much as weather prediction models adjust to new observations. They form a consensus, capturing atmospheric behaviors that have eluded all the separate models. We demonstrate that simplified supermodels succeed, where no single approach can, by correctly representing
critical phenomena involving sudden qualitative transitions, such as occur in El Niño events, that depend on interactions among atmospheric processes on many different scales in space and time. The correct reproduction of critical phenomena is vital both for predicting weather and for projecting the effects of climate change.
