Abstract:To explain the differentiation of stem cells in terms of dynamical systems theory, models of interacting cells with intracellular protein expression dynamics are analyzed and simulated. Simulations were carried out for all possible protein expression networks consisting of two genes under cell--cell interactions mediated by the diffusion of a protein. Networks that show cell differentiation are extracted and two forms of symmetric differentiation based on Turing's mechanism and asymmetric differentiation are identified. In the latter network, the intracellular protein levels show oscillatory dynamics at a single-cell level, while cell-to-cell synchronicity of the oscillation is lost with an increase in the number of cells. Differentiation to a fixed-point type behavior follows with a further increase in the number of cells. The cell type with oscillatory dynamics corresponds to a stem cell that can both proliferate and differentiate, while the latter fixed-point type only proliferates. This differentiation is analyzed as a saddle-node bifurcation on an invariant circle, while the number ratio of each cell type is shown to be robust against perturbations due to self-consistent determination of the effective bifurcation parameter as a result of the cell--cell interaction. Complex cell differentiation is designed by combing these simple two-gene networks. The generality of the present differentiation mechanism, as well as its biological relevance, is discussed.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to explain the process of stem cell differentiation through the theory of dynamical systems. Specifically, the author constructed a model of cell - cell interactions and the dynamics of intracellular protein expression, and carried out simulation analysis. The main objectives of the paper can be summarized as follows:
1. **Understand how stem cells simultaneously support the two seemingly contradictory functions of proliferation and differentiation**:
- Stem cells can self - replicate (i.e., keep their components unchanged), and at the same time can also differentiate into other types of cells (i.e., transform into states with different components). How to explain this phenomenon within the framework of dynamical systems is a key issue.
2. **Determine how initial conditions select different attractor states**:
- During development, how cells transition from an initial state to different stable states (i.e., different cell types) is another question that needs to be answered.
3. **Explore the stability of the cell - differentiation process**:
- Including the stability of the time points of cell differentiation and the distribution of the number of each cell type, and whether these stability mechanisms are related to cell - cell interactions.
To achieve the above objectives, the author developed a coupled dynamical - system model, which includes changes in intracellular protein concentrations, cell - cell interactions, and an increase in the number of cells due to cell division. By simulating all possible protein - expression networks composed of two genes, the author extracted the networks that can exhibit cell differentiation and further analyzed the dynamical characteristics of these networks.
### Main Findings
1. **Turing - type Differentiation (Fixed Point → Fixed Point)**:
- The dynamics of a single cell has a unique stable fixed point. As the number of cells increases, the cell population splits into two groups, each with different fixed - point values. This mechanism can be explained by the classical Turing pattern.
2. **Oscillatory Death (Oscillation → Fixed Point)**:
- The dynamics of a single cell has a limit - cycle attractor. As the number of cells increases, the cell population again splits into two groups, and each group shows different fixed points through cell - cell interactions. This phenomenon of losing oscillation is called "oscillatory death".
3. **Asymmetric Differentiation (Oscillation → Oscillation Transition)**:
- The dynamics of a single cell is still a limit cycle, but as the number of cells increases, the synchronous oscillation becomes unstable, and some cells enter an oscillation state with a small amplitude, while other cells remain on the original limit cycle. These two oscillation states stably co - exist with each other and are not affected by the initial conditions.
### Dynamical Mechanisms
The mechanism of asymmetric differentiation is based on two stages: oscillation desynchronization and saddle - node - on - invariant - circle (SNIC) bifurcation.
1. **Desynchronization**:
- When there are cell - cell interactions, the oscillations of cells lose synchronization. This can be verified by calculating the synchronization - stability index, and the results show that the magnitude of the tangent vector representing the deviation between cells grows exponentially with time.
2. **Differentiation**:
- After desynchronization, cell - cell interactions lead to differentiation. Through the diffusion term \( \alpha(t) \equiv D\left(\frac{1}{N} \sum_{k = 1}^N y_k - y_i\right) \), the dynamical equations become:
\[
\frac{dx(t)}{dt} = \frac{1}{1 + \exp(-\beta(x(t) - y(t) - g_x))} - x(t)
\]
\[
\frac{dy(t)}{dt} = \frac{1}{1 + \exp(-\beta(x(t) - g_y))} - y(t) + \alpha
\]
- If the oscillations are synchronous, then \( \alpha \equiv 0 \), and the equations degenerate into single - cell dynamics. Due to desynchronization, \( \alpha \) is not zero and, as a time - dependent bifurcation parameter, causes the single - cell dynamics to bifurcate. When \( \alpha \) decreases, the zero - line of \( \frac{dy}{dt} = 0 \) moves downward and intersects with the zero - line of \( \frac{dx}{dt} = 0 \), forming an SNIC bifurcation