
This study proposed an interpretable multi-scale infrared small object detection network (IMD-Net) design method to improve the precision of infrared small object detection and contour segmentation in complex backgrounds. To this end, a multi-scale object enhancement module was constructed, which converted artificially designed features into network structures. The network structure was used to enhance actual objects and extract shallow detail and deep semantic features of images. Next, a global object response, channel attention, and multilayer feature fusion modules were introduced, combining context and channel information and aggregated information, selected data, and decoded objects. Finally, the multiple loss constraint module was constructed, which effectively constrained the network output using multiple losses and solved the problems of high false alarms and high missed detections. Experimental results showed that the proposed network model outperformed local energy factor (LEF), self-regularized weighted sparse model (SRWS), asymmetric contextual modulation (ACM), and other state of the art methods in the intersection-over-union (IoU) and Fmeasure values by 10.8% and 11.3%, respectively. The proposed method performed best on the currently available datasets, achieving accurate detection and effective segmentation of dim and small objects in various infrared complex background images.
Citation: Dawei Li, Suzhen Lin, Xiaofei Lu, Xingwang Zhang, Chenhui Cui, Boran Yang. IMD-Net: Interpretable multi-scale detection network for infrared dim and small objects[J]. Mathematical Biosciences and Engineering, 2024, 21(1): 1712-1737. doi: 10.3934/mbe.2024074
[1] | Jihong Yang, Hao Xu, Congshu Li, Zhenhao Li, Zhe Hu . An explorative study for leveraging transcriptomic data of embryonic stem cells in mining cancer stemness genes, regulators, and networks. Mathematical Biosciences and Engineering, 2022, 19(12): 13949-13966. doi: 10.3934/mbe.2022650 |
[2] | Katrine O. Bangsgaard, Morten Andersen, Vibe Skov, Lasse Kjær, Hans C. Hasselbalch, Johnny T. Ottesen . Dynamics of competing heterogeneous clones in blood cancers explains multiple observations - a mathematical modeling approach. Mathematical Biosciences and Engineering, 2020, 17(6): 7645-7670. doi: 10.3934/mbe.2020389 |
[3] | E.V. Presnov, Z. Agur . The Role Of Time Delays, Slow Processes And Chaos In Modulating The Cell-Cycle Clock. Mathematical Biosciences and Engineering, 2005, 2(3): 625-642. doi: 10.3934/mbe.2005.2.625 |
[4] | Samantha L Elliott, Emek Kose, Allison L Lewis, Anna E Steinfeld, Elizabeth A Zollinger . Modeling the stem cell hypothesis: Investigating the effects of cancer stem cells and TGF−β on tumor growth. Mathematical Biosciences and Engineering, 2019, 16(6): 7177-7194. doi: 10.3934/mbe.2019360 |
[5] | Awatif Jahman Alqarni, Azmin Sham Rambely, Sana Abdulkream Alharbi, Ishak Hashim . Dynamic behavior and stabilization of brain cell reconstitution after stroke under the proliferation and differentiation processes for stem cells. Mathematical Biosciences and Engineering, 2021, 18(5): 6288-6304. doi: 10.3934/mbe.2021314 |
[6] | Heyrim Cho, Ya-Huei Kuo, Russell C. Rockne . Comparison of cell state models derived from single-cell RNA sequencing data: graph versus multi-dimensional space. Mathematical Biosciences and Engineering, 2022, 19(8): 8505-8536. doi: 10.3934/mbe.2022395 |
[7] | R. A. Everett, Y. Zhao, K. B. Flores, Yang Kuang . Data and implication based comparison of two chronic myeloid leukemia models. Mathematical Biosciences and Engineering, 2013, 10(5&6): 1501-1518. doi: 10.3934/mbe.2013.10.1501 |
[8] | Qiaojun Situ, Jinzhi Lei . A mathematical model of stem cell regeneration with epigenetic state transitions. Mathematical Biosciences and Engineering, 2017, 14(5&6): 1379-1397. doi: 10.3934/mbe.2017071 |
[9] | J. Ignacio Tello . On a mathematical model of tumor growth based on cancer stem cells. Mathematical Biosciences and Engineering, 2013, 10(1): 263-278. doi: 10.3934/mbe.2013.10.263 |
[10] | Fabien Crauste . Global Asymptotic Stability and Hopf Bifurcation for a Blood Cell Production Model. Mathematical Biosciences and Engineering, 2006, 3(2): 325-346. doi: 10.3934/mbe.2006.3.325 |
This study proposed an interpretable multi-scale infrared small object detection network (IMD-Net) design method to improve the precision of infrared small object detection and contour segmentation in complex backgrounds. To this end, a multi-scale object enhancement module was constructed, which converted artificially designed features into network structures. The network structure was used to enhance actual objects and extract shallow detail and deep semantic features of images. Next, a global object response, channel attention, and multilayer feature fusion modules were introduced, combining context and channel information and aggregated information, selected data, and decoded objects. Finally, the multiple loss constraint module was constructed, which effectively constrained the network output using multiple losses and solved the problems of high false alarms and high missed detections. Experimental results showed that the proposed network model outperformed local energy factor (LEF), self-regularized weighted sparse model (SRWS), asymmetric contextual modulation (ACM), and other state of the art methods in the intersection-over-union (IoU) and Fmeasure values by 10.8% and 11.3%, respectively. The proposed method performed best on the currently available datasets, achieving accurate detection and effective segmentation of dim and small objects in various infrared complex background images.
Embryonic stem cells (ESCs) are pluripotent stem cells that are derived from the inner cell mass of a blastocyst which is an early-stage preimplantation embryo that can be propagated by culturing in an undifferentiated state while maintaining the capacity to generate any cell type in the body[1]. Due to their plasticity and potentially unlimited capacity for self-renewal, human embryonic stem cells (hESCs) play an important role in animal cloning, organ transplants, etc. For instance, embryonic stem cell therapies have been proposed for regenerative medicine and tissue replacement after injury or disease[2,3].
The transition from embryonic stem cells to three different types of cells in the ectoderm, mesoderm and endoderm has been the first step in studying how pluripotent cells exit the pluripotent state and give rise to lineage-specific progenitors. From the perspective of molecular biology, the transition is determined by different levels of gene expression. Thus, understanding the regulatory mechanism of the genes underlying differentiation and the cell fate decision of hESCs has become necessary.
In the last few years, genome-wide profiling approaches have started to uncover the molecular programs that drive the developmental processes of hESCs. A number of research works have revealed the key genes or mechanisms underlying the differentiation of hESCs[4,5,6,7,8,9,10]. For example, Adamo et al. revealed that LSD1 could regulate the balance between self-renewal and differentiation in hESCs [4]. Cheryle et al. demonstrated that stable endoderm progenitors can be established from hESCs by the constitutive expression of SOX7 or SOX17, producing extraembryonic endoderm and definitive endoderm progenitors, respectively[9]. Ivanova et al. found that OCT4 regulates and interacts with the BNP4 pathway to specify four developmental fates and identified general and cell-line-specific requirements for NANOG, OCT4, and SOX2 in the hESC[10].
Recently, due to the development of single-cell sequencing technology, a new type of data named "single-cell RNA-seq data (scRNA-seq data)" is available and we can now quantify the gene expression of individual cells and analyze the heterogeneity among cells[11,12]. Through analyzing the scRNA-seq data of hESCs, Chu et al. revealed that the Kr¨uppel-like factor 8 (KLF8) plays a key role in modulating the mesendoderm, which is an intermediate state before the definitive endoderm and mesoderm are formed, to DE state differentiation[5].
However, the dynamic molecular mechanisms underlying the differentiation of individual human embryonic stem cells are still poorly understood. Moreover, how KLF8 dynamically interacts with other important genes in the core regulatory network and governs the transition from ESCs to the three different cell types in the ectoderm, mesoderm and endoderm remains unclear[13,14,15].
Mathematical modeling has been demonstrated as a powerful tool for investigating the dynamic mechanisms underlying signal regulatory network. To systematically analyze the dynamical regulatory mechanism underlying hESC differentiation, in this paper, we developed a mathematical model that is based on a possible core regulatory network underlying the differentiation process of hESCs. scRNA-seq data and the particle swarm optimization (PSO) algorithm are used to identify the parameters of the mathematical model. Then, we performed the simulation and dynamic bifurcation analysis that is based on the proposed model and investigated the molecular mechanisms underlying the core regulatory module. More importantly, we proved the reasonability of the proposed regulatory mechanisms through the dynamic analysis results and introduced several strategies that may be feasible for controlling the process of hESC differentiation.
It has been documented that epithelial-mesenchymal transition (EMT) underlies cell fate conversions in both reprogramming and differentiation along an endoderm cell fate[16]. HESCs begin differentiation with a near-synchronous EMT, and the differentiation of hESCs is considered an EMT process[17]. In fact, hESCs are epithelial cells with a high expression of E-cadherin (CDH1), which plays a central role in maintaining epithelial cell-cell adhesion and polarity. Previous work[18,19] has shown that down regulation of CDH1's transcription is thought to be a primary mechanism that contributes to the onset of EMT. The CDH1 has also been demonstrated to be the the paradigm of epithelial genes. Loss of CDH1 is normally determined to be the hallmark of EMT[20].
In the differentiation of hESCs, KLF8 is recognized to be an essential regulator in modulating mesendoderm to endoderm differentiation. It was observed that overexpression of KLF8 increases the mobility, which was evident by the up regulation of TWIST1, an EMT marker, which indicated that KLF8 plays a specific role in promoting the transition from mesendoderm to endoderm state[5]. Indeed, KLF8 induces EMT by directly binding to the CDH1 promoter through GT boxes and repressing the expression of CDH1[21].
Zinc finger E-box-binding homeobox1 (ZEB1) has been well studied in cancerous tissues, such as cervical cancer, gastric cancer, lung cancer[22,23,24]. In addition, it was demonstrated that ZEB1 is an important transcription factor of EMT by inhibiting an E-box gene, including CDH1[23,25,26].
To obtain the appropriate model, we first used single-cell RNA-seq data to compute Spearman's correlation coefficients among three genes. The calculated Spearman's correlation coefficients between ZEB1 and KLF8, ZEB1 and CDH1, and CDH1 and KLF8 are 0.4879, -0.4664 and -0.4851, respectively.
According to the above descriptions, it is obvious that ZEB1 inhibits CDH1 and KLF8 inhibits CDH1. No evidence shows how KLF8 influences ZEB1 and CDH1. It has also been suggested that KLF8 functions as a promoter of DE cells by upregulating many other DE cell marker genes, such as CXCR4[5]. We assumed several models to describe their different relations and determined that good hypotheses are that ZEB1 is promoted by KLF8 and CDH1 inhibits the expression of KLF8. The connections among the three genes are contained in the following Figure 1.
The single-cell RNA-seq data we used in this research is obtained from NCBIs Gene Expression Omnibus and are accessible through GEO series accession number GSE75748, which characterizes differentiation from the mesendoderm(a state before the formation of DE cells and mesoderm cells) to definitive endoderm cells[5]. This data set contains a large number of genes expression data collected at 6 time points with 758 cells.
Based on the single cell RNA-seq data, we first reconstructed a single-cell order using the Wave-Crest toolbox under R software. The cell order assigns a specific time point to each cell within the time interval [0 h, 96 h]. There are 758 time points evenly distributed from 0 h to 96 h (each pair of adjacent time points is 96 h/757 apart according to the Wave-Crest algorithm). These 758 time points correspond to the 758 reordered cells. The cell order produces a time-series expression of genes, which was used to calculate the Spearman's correlation coefficients (Table 1). Then, we obtained a pseudo trajectory of the gene-specific expression through polynomial fitting. The normalized scRNA-seq data is obtained by evenly selecting 758 points of the polynomial curve with the same time points. It has been demonstrated that the time series expression of the three considered genes reflect the dynamical differentiation well[5]. After further investigation, we discovered that the parameters estimated by the data not only capture the differentiation from the mesendoderm state to the DE state, but, in fact, changing some of the parameters can simulate the differentiation from ESCs to all three germ layers. We speculate that the function of genes and the interactions between genes at different stages of differentiation remain unchanged. The differentiation trend is mainly determined by specific biological processes.
Genes | CDH1&ZEB1 | CDH1&KLF8 | ZEB1&KLF8 |
Spearman's correlations | –0.4664 | –0.4851 | 0.4879 |
We built a mathematical model based on mass action law and Michaelis-Menten equation[27,28], which is described in the following ordinary differential equations (ODEs). The relationships among the three selected genes has been thoroughly discussed in the network construction, and each term of the equations has been marked accordingly.
d[CDH1]dτ=k1k2+[ZEB1]n1⏟ZEB1 inhibits the transcription of CDH1 +k3k4+[KLF8]n2⏟KLF8 inhibits the transcription of CDH1−d1[CDH1]⏟degradation of CDH1 | (1) |
d[ZEB1]dτ=a×k5[KLF8]n3k6+[KLF8]n3⏟assume KLF8 promotes transcription of ZEB1−d2[ZEB1]⏟degradation of ZEB1 | (2) |
d[KLF8]dτ=r×k7k8+[CDH1]n4⏟assume CDH1 inhibits transcription of KLF8−d3[KLF8]⏟degradation of KLF8 | (3) |
Here k1,k3,k5 and k7 represent the maximum production rate of CDH1 regulated by ZEB1, the maximum production rate of CDH1 regulated by KLF8, the maximum production rate of ZEB1 regulated by KLF8 and the maximum production rate of KLF8 regulated by CDH1, respectively. k2,k4,k6 and k8 represent the Michaelis constants that correspond to the production rates k1,k3,k5 and k7, respectively. d1,d2 and d3 represent the degradation rates of CDH1, ZEB1 and KLF8, respectively. a and r represent weight coefficients to reflect the degree of activation and inhibition, respectively.
To reduce the parameters and make the following analysis more convenient, we nondimensionalized the model by making the following substitutions, whereby we assume that the Hill coefficients are equal to 2.
y1=k2d1k1[CDH1],y2=k6d1k4k5[ZEB1],y3=[KLF8]√k4 | (4) |
Therefore, the simplified ODE model is as following, and all the analysis is based on the nondimensionalized ODEs.
dy1dt=11+p1y22+p21+y23−y1 | (5) |
dy2dt=a×y231+p3y23−p4y2 | (6) |
dy3dt=r×p51+p6y21−p7y3 | (7) |
where
p1=k24k25k26d21,p2=k2k3k1k4,p3=k4k6,p4=d2d1,p5=k22k7d1k21√k4,p6=k22k8d21k21,p7=d3d1 | (8) |
The model includes three variables and night parameters. To identify the parameters in the model, we converted the parameter identification into the problem of optimization by minimizing the following objective function. Mathematically, the objective function is defined as the error between the simulation results and the time series experimental data. The formulation can be expressed as
minPJ(P)=αK∑k=1J∑j=1(yDk(tj)−yk(tj,P)maxj(yDk(tj)))2 | (9) |
yDk(tj) represents the measured data of component k at time-point tj, which, in our case, is the normalized data obtained from the polynomial curves. yk(tj,P) represents the kth component of the solutions of ODE at time tj with parameter set P. Here J=758, since there are 758 pseudo time points. The numerical solution of the ODEs is solved by MATLAB, with the initial value set to be the approximation of the first data point of the normalized data and a,r assumed to be 1. The particle swarm optimization(PSO)[29] algorithm is used to optimize this object function within the optimization intervals that are carefully chosen. The obtained optimization parameters are listed in Table S1.
Processes | Parameters | Values | Remarks |
Associate constant | p1 | 49.9611 | Fitted |
Associate constant | p2 | 0.5661 | Fitted |
Relative dissociation constant | p3 | 9.4701 | Fitted |
Relative degradation constant | p4 | 0.0951 | Fitted |
Associate constant | p5 | 2.3591 | Fitted |
Associate constant | p6 | 11.9229 | Fitted |
Relative degradation constant | p7 | 0.9982 | Fitted |
Hill coefficient | n1 | 2 | assumed |
Hill coefficient | n2 | 2 | assumed |
Hill coefficient | n3 | 2 | assumed |
Hill coefficient | n4 | 2 | assumed |
Intensity of impact | a | 1 | assumed |
Intensity of impact | r | 1 | assumed |
We used the proposed model to simulate the dynamic process with time. Figure 2a, c, e demonstrates that the simulation results are consistent with the obtained trajectory of single-cell RNA-seq data, which shows that our model could reflect the main biological process underlying the differentiation of ESC. The fitting of 758 scRNA-seq data is shown in Appendix D.
Since the RNA-seq data we obtained was collected from large amount of cells which may lead to the errors between the real initial gene expression level and the initial expression level we estimated, we gave the initial value random perturbations within 0.3 (Figure 2b, d, f). The results show that the model well reflects the fluctuation of the expression data with different initial values. The simulation is independent of the selection of initial values to a certain extent.
We presented dynamic simulations of the system with different values of p4. The curves in Figure 3a indicate that the dimensionless concentration of CDH1 gradually decreases after a short increase until a steady state is attained during the time evolution. The steady state remains the same with an increase of p4. Intriguingly, the steady state significantly increases into an extremely high state when p4 is larger than 0.133. Conversely, a relatively lower state is maintained when p4 is less than 0.133. It is obvious that the system will sustain more steady states with the change of p4.
Through the nondimensionalization process, we know that p4 is the ratio of the degradation rate of ZEB1 to the degradation rate of CDH1, which we named the "relative degradation rate". We hypothesize the existence of thresholds for the relative degradation rate, which eventually distinguishes different levels of CDH1 expression.
The threshold of the relative degradation rate is clearly indicated by Figure 3b. Through qualitative and quantitative analyses, we suggest that a high expression of CDH1, in which relative degradation rate is greater than 0.225, corresponds to the ectoderm state of ESCs. With the relative degradation rate decreasing, the inhibition of ZEB1 targets on CDH1 is enhanced, at least partly, because of the relatively lower degradation rate of ZEB1, which results in weaker adhesion among cells. Thus, EMT is activated to a certain extent, and the system undergoes SN bifurcation and another stable steady state occurs, which we believe corresponds to the mesoderm state. As the relative degradation rate decreases further (<0.096), EMT is further activated and the system undergoes another SN bifurcation. Then, the lower expression state of CDH1 occurs, which corresponds to the formation of DE cells. In short, the threshold we discovered is the necessary condition for cell formation in different germ layers. The results also indicate that mesendoderm cells can repossess pluripotency and differentiate into ectodermal cells, and the differentiation process is reversible under certain conditions. There are reasons to believe that the proposed model is able to characterize the differentiation from ESCs to the three embryonic layers.
We investigated the effect of inhibition intensity on the dynamic behavior of the system. Figure 4a shows that a bistable phenomenon exists when p3 varies in the region [9.65, 23]. According to the nondimensionalization process (Eq 8), p3 is the ratio of k4 and k6, which are the dissociation constants that KLF8 targets on CDH1 and ZEB1, respectively. Since p3 is equal to the ratio of the dissociation rate of the ligand-receptor complex, we named it as the relative dissociation constant, and we found that there is a threshold of the relative dissociation constant. When the dissociation constant of the ligand-receptor complex of CDH1 and KLF8 is approximately 23 times greater than that of ZEB1 and KLF8, the system undergoes a bistable switch from a low-expression state to a high-expression state of CDH1. The bifurcation diagram indicates that the change of p3 might not be a good way to obtain all three different types of cells during ESC differentiation.
However, we suggest that influencing the relative dissociation constant might be a feasible way to obtain DE cells from cells in other germ layers.
To further determine how p3 affects the whole system, we analyzed the two-parameter bifurcation. According to Figure 4b, when the relative dissociation constant p3 decreases to a small enough amount (for example, p3=0.3), the tristable phenomenon disappears, and the system undergoes an irreversible bistable switch, which means that it is impossible to obtain mesoblastema and hard to change the fate of cells from differentiated cell into DE cells. In this case, CDH1 maintains a low level of expression, and KLF8 maintains a high level of expression. Interestingly, when the relative dissociation constant becomes greater, the system exhibits only one stable steady state of high concentration of CDH1.
Qualitatively we find that a large dissociation constant of the ligand-receptor complex of CDH1 and KLF8 and a relatively smaller dissociation constant of the ligand-receptor complex of ZEB1 and KLF8 resulted in the activation of CDH1 (the case where p3=35), which is partly due to the situation that the inhibitory effect of CDH1 by KLF8 is weakened by a large dissociation constant. Instead, with a smaller dissociation constant of ligand-receptor complex of CDH1 and KLF8 compared to that of ZEB1 and KLF8, the direct inhibition of CDH1 by KLF8 is enhanced (the case where p3=0.2), which leads to a state of low expression of CDH1, although ZEB1 is promoted by KLF8.
Furthermore we have the following equation according to Eq 8.
p1=p23k25d21 |
Under the circumstance that p3 and d1 are constants, the decrease of p1 results from the the decrease of k5, which characterizes the promotion of ZEB1 by KLF8. It reveals that the system may exhibit only one steady state if the relative dissociation constant p3 is too large or too small, no matter how the production rate of ZEB1 changes.
Because p3 and p4 play important roles in the ESC fate decision, we start to wonder how they synergistically affect the system. Figure 4c, d shows the varying steady states of the model with the two constants in different domains. From the parameter values of p3,p4 that correspond to different steady states, we suggest that it is difficult to obtain mesodermal cells with the relatively weakened inhibition effect that KLF8 targets on CDH1 or the weakened promotion that KLF8 targets on ZEB1. Meanwhile, a high relative degradation rate of ZEB1 to CDH1 determines the irreversible steady state of epithelial cell formation, which corresponds to the differentiation of the ectoderm according to the yellow district of Figure 4d.
In conclusion, the relative dissociation constant p3 should be neither too large nor too small in case of an irreversible situation of the steady states dominated by KLF8 if three types of embryonic cells are needed.
Previous experimental and dynamical analyses show that KLF8 plays an important role in modulating cell differentiation from mesendoderm to DE. To systematically analyze the role of KLF8 in the differentiation, we plotted the bifurcation diagram (Figure 5) with respect to the inhibition coefficient that CDH1 targets on KLF8 (p6). With an increase of the inhibition coefficient p6, CDH1 switches from the monostability of low level to tristable states in a narrow region, then to the high level of bistable states. It reveals that if we want to switch the DE state to other states of hESC, we may enhance the inhibition strength of CDH1 on KLF8, or artificially, we inhibit the expression of KLF8 through medication. Combined with the fact that KLF8 inhibit the expression of CDH1, we assume that KLF8 promotes the formation of DE cells possibly by mutual inhibition of KLF8 and CDH1.
The above analysis is based on the assumption that a=1,r=1. To further investigate the function of KLF8 in the system, we decreased parameter a to characterize different intensities of CDH1 on KLF8 and plotted Figure 6a to analyze the behavior of KLF8.
With the decrease of a, the tristability phenomenon gradually disappears, and the system undergoes an irreversible trend towards a high expression level of CDH1 (when a reaches the threshold 0.15), which marks the formation of epithelial cells. Similarly, the decrease of r leads to the disappearance of the tristability phenomenon of the system (Figure 6b), which finally resulted in the formation of the ectoderm. The simulation results show similar conclusions (see Figure S2). We reasonably suggest that the promotion of ZEB1 by KLF8 and the inhibition of KLF8 by CDH1 is indispensable for DE cell formation which is consistent with the complex biological process. More specifically, we assume that the promotion of ZEB1 by KLF8 does not exist, which results in changes of the promotion term from a nonlinear function of KLF8 to a linear constant c. The simulation results (Figure 6c) indicate that when the KLF8-related production rate of ZEB1 degenerates into constants c, the expression of CDH either maintains a high level of expression or reduces it directly, which is inconsistent with experimental observations. Figure 6d shows the same conclusion that the production rate of KLF8 is a nonlinear function of CDH1, as was indicated by our model.
In fact, the intrinsically complex processes require the coordinated dynamic expression of hundreds of genes and proteins in precise response to external signaling cues. Our work only focuses on a fraction of the complex regulatory network, which raises another question about whether the model really reflects the biological differentiation process. To verify the generality and reliability of the proposed model, we analyzed it from three aspects. First, we tried numbers of model settings, e.g., we set the Hill coefficients to 1 and found that the numerical simulations of the model did not fit the data well. We have also tried the linear functions in the model as a regulatory method between two genes. The fitting results remain poor. Secondly, to testify whether the parameters are sensitive to the single-cell data we used in the parameter estimation, we performed 10-fold cross-validation (Appendix E). The parameters that were estimated by 10 data subsets vary slightly and are similar to the parameters listed in Table S1. Third, since cell differentiation is typically heterogeneous and is spatially disordered, the intrinsic fluctuations and extrinsic signal fluctuations may play important roles in modulating the state switch in the differentiation of hESCs. To explore the influence of intrinsic and extrinsic noise on the transition between multiple steady states of the proposed model, we developed corresponding Master equation and Langevin equation models (Appendix F) and performed stochastic simulations using the Gillespie algorithm and Euler-Maruyama algorithm[30], respectively. According to the simulation results (Figures S4 and S5), the Master equation and Langevin equation models displayed similar dynamic behavior as that by the deterministic model.
The above evidences proved that the proposed model can characterize the differentiation from hESCs to three germ layers and we suggest that KLF8 promotes the formation of DE cells possibly with the promotion of ZEB1 and the mutual inhibition of KLF8 and CDH1, which partially answers the unsolved question[5].
In this study, based on the single-cell RNA-seq data, we reconstructed a core regulatory network underlying the stem-cell differentiation process, and demonstrated that the core regulatory module shows various behaviors, including the numerical fitting and the three states switch, which is in correspondence with the three types of cell in three germ layers. Thus, we proved that the novel important gene KLF8[5] affects the differentiation process by up-regulating ZEB1 and down-regulating the expression of CDH1, which forms a coupled feedback loop. In addition, we defined two indexes including the relative degradation rate and the relative disassociation rate which are tightly associated with the complex dynamic behaviors.
We also proposed some possible methods to realize the switch of the cell fate.
● The small relative dissociation constant ensures the formation of DE cells, and the large relative dissociation constant ensures the formation of epithelial cells;
● The inhibition of KLF8 contributes to the formation of epithelial cells by blocking EMT;
● The high relative degradation rate of ZEB1 to CDH1 determines the steady state of epithelial cell formation.
Altogether, we believe that the combination of scRNA-seq analysis and mathematical modeling can well reveal the molecular mechanism of cell fate decisions.
This work was supported by the key project of the National Natural Science Foundation of China (No. 11831015) and the Chinese National Natural Science Foundation (No. 61672388 and No. 61802125) and the Natural Science Foundation of Jiangxi Province (No. 20181BAB202006).
All authors declare no conflicts of interest in this paper.
The global sensitivity of parameters reflects how the system responds to the perturbation of parameters in the model. To obtain the sensitivity of the input parameters to all variables in the model, the sensitivity function sj(t) of the parameter pj at time t was defined as follows[31,32]:
sj=∂O(t)O(t)/∂Pj(t)Pj(t)≈‖O(Pj+ΔPj,t)−O(Pj−ΔPj,t)‖O(Pj,t)/2ΔPjPj |
where O(t) is the model output at time t (T is the total time length), ΔPj is a small perturbation which is 10% in our situation. Sj=∫T0sj(t)dt is the sensitivity value of parameter Pj. From Figure S1, we can see that the perturbation of parameter p2,p5 and p7 have relative large effects on the expression of CDH1 and KLF8.
We simulated the dynamics of CDH1 with different weights that characterize different interactive intensity of CDH1, ZEB1 and KLF8. The variation tendency of the curve indicates that the promotion of ZEB1 by KLF8 and the inhibition of KLF8 by CDH1 is necessary for the formation of DE cells in the differentiation process of ESCs, which is marked by a low expression of CDH1.
To clearly observe the original 758 single cell gene expression data and to thoroughly investigate the fitting of the model to the data, we plotted the following three graphs. The results indicate that the model reflects the main variation tendency of the gene expression levels. However, the simulation results do not fit the expression levels of ZEB1 and KLF8 well at approximately 90h. We assumed that the cells completed the differentiation process at that time, and the expression of ZEB1 and KLF8 were down-regulated for some reason that are not captured by the model. New regulatory genes need to be considered to study the regulatory mechanisms at the late stage of differentiation. Nevertheless, the present model can serve as a significant foundation for further investigation.
Through 10-fold validation, the scRNA-seq data were randomly and evenly divided into 10 parts. 9 of the data subsets were selected to estimate the parameters each time. After 10 optimizations, 10 sets of estimated parameters were listed in Table S2. The parameters estimated by 10 data subsets vary slightly and are similar to the parameters listed in Table S1.
Validation | p1 | p2 | p3 | p4 | p5 | p6 | p7 | Error |
1 | 49.99 | 0.44 | 17.24 | 0.06 | 48.88 | 15.96 | 22.47 | 12.76 |
2 | 29.98 | 0.44 | 14.06 | 0.06 | 5.08 | 12.51 | 2.56 | 12.63 |
3 | 50.00 | 0.46 | 16.98 | 0.06 | 38.76 | 15.35 | 17.48 | 12.63 |
4 | 49.98 | 0.48 | 15.06 | 0.07 | 44.83 | 14.48 | 20.13 | 12.51 |
5 | 50.00 | 0.47 | 11.68 | 0.11 | 31.98 | 12.80 | 14.41 | 12.88 |
6 | 49.95 | 0.44 | 17.82 | 0.06 | 49.99 | 16.21 | 22.83 | 12.71 |
7 | 50.00 | 0.45 | 18.01 | 0.06 | 47.61 | 15.95 | 21.04 | 12.63 |
8 | 50.00 | 0.44 | 17.00 | 0.06 | 23.69 | 16.05 | 10.81 | 12.75 |
9 | 50.00 | 0.44 | 18.94 | 0.05 | 47.20 | 16.23 | 21.81 | 12.74 |
10 | 49.99 | 0.46 | 16.57 | 0.06 | 2.27 | 14.49 | 1.04 | 12.65 |
Average | 47.99 | 0.45 | 16.34 | 0.07 | 34.03 | 15.00 | 15.46 | 12.69 |
We supposed that the state of the system is X(t0)=x0 and defined the conditional probability density function P(x,t|x0,t0), which is the probability of the system state satisfying X(t)=x at time t. The corresponding Master equation is developed according to the following probability equation:
∂∂tP(x,t|x0,t0)=6∑j=1(aj(x−vj)P(x−vj,t|x0,t0)−aj(x)P(x,t|x0,t0)) |
where aj is the propensity function corresponding to the chemical reaction channel j, and vj is the state-change vector when reaction j occurred. For simplicity, we used the Hill function that was defined in the deterministic model to depict the propensities instead of decomposing this into a set of elementary reaction steps. Following the modeling approaches for stochastic models[33], the propensities were defined in Table S3, and the parameters were set to be the same as those that are defined in Table S1.
NO. | Reaction | Propensity function |
1 | → CDH1 | Ω3/(Ω2+p1×[ZEB1]2)+p2×Ω3/(Ω2+[KLF8]2) |
2 | CDH1 → | [CDH1] |
3 | → ZEB1 | Ω×[KLF8]2/(Ω2+p3×[KLF8]2) |
4 | ZEB1 → | p4×[ZEB1] |
5 | → KLF8 | p5×Ω3/(Ω2+p6×[CDH1]2) |
6 | KLF8 → | p7×[KLF8] |
The stochastic simulation was performed using the Gillespie algorithm[30], and a steady state was obtained by dividing the system size by the molecular numbers. According to the simulation results (Figure S4), the system displayed similar dynamic behavior as that of the deterministic model.
The Langevin equations were built as follows, where ζi is the Gaussian white noises with ζi(t)>0,(ζi(t),ζj(t′))≥δijδ(t−t′), and similar results were obtained under a simulation of the chemical Langevin equations (Figure S5).
[CDH1](t+dt)=[CDH1](t)+11+p1[ZEB1]2+p21+[KLF8]2−[CDH1]+√11+p1[ZEB1]2+p21+[KLF8]2ζ1(t)(dt)12−[CDH1]ζ2(t)(dt)12 |
[ZEB1](t+dt)=[ZEB1](t)+a[KLF8]21+p3[KLF8]2−p4[ZEB1]+√a[KLF8]21+p3[KLF8]2ζ3(t)(dt)12−√p4[ZEB1]ζ4(t)(dt)12 |
[KLF8](t+dt)=[KLF8](t)+rp51+p6[CDH1]2−p7[KLF8]+√rp51+p6[CDH1]2ζ5(t)(dt)12−√p7[KLF8]ζ6(dt)12 |
Previous studies have demonstrated that the bistable systems can generate bimodal expression patterns of corresponding genes[34,35]. To explore how the negative feedback strength shapes the KLF8 expression distribution, we simulated the distribution of KLF8 in a collection of 10,000 cells under different values of p6 for the Master equation. Figure S6 gives the distribution diagram of KLF8 switching from a unimodal distribution of the high state (p6=13) via the bimodal distribution in the tristable region (p6=15), to the bimodal distribution in the bistable region (p6=17) and, then, to the low state with unimodal distribution (p6=22). Since the overexpression of KLF8 could accelerate the transition from the mesendoderm cell to DE cell under differentiation[5], the distribution of KLF8 become biomodal (p6=15) distribution, which may correspond to the undifferentiated state or differentiate from mesendoderm to DE state. Or the cell may directional differentiate from mesendoderm state to endoderm state (p6=13).
[1] |
S Wu, K Zhang, S Li, J Yan, Joint feature embedding learning and correlation filters for aircraft tracking with infrared imagery, Neurocomputing, 450 (2021), 104–118. https://doi.org/10.1016/j.neucom.2021.04.018 doi: 10.1016/j.neucom.2021.04.018
![]() |
[2] | N. Zou, J. W. Tian, Research on multi feature fusion infrared ship wake detection method, Comput. Sci., 45 (2018), 172–175. |
[3] |
C. Deng, S. He, Y. Han, B. Zhao, Learning dynamic spatial-temporal regularization for UAV object tracking, IEEE Signal Process. Lett., 6 (2021), 1230–1234. https://doi.org/10.1109/LSP.2021.3086675 doi: 10.1109/LSP.2021.3086675
![]() |
[4] |
Y. Han, H. Liu, Y. Wang, C. Liu, A comprehensive review for typical applications based upon unmanned aerial vehicle platform, IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens., 15 (2022), 9654–9666. https://doi.org/10.1109/JSTARS.2022.3216564 doi: 10.1109/JSTARS.2022.3216564
![]() |
[5] |
Y. Han, H. Wang, Z. Zhang, W. Wang, Boundary-aware vehicle tracking upon UAV, Electron. Lett., 8 (2020), 873–876. https://doi.org/10.1049/el.2020.1170 doi: 10.1049/el.2020.1170
![]() |
[6] |
X. Liang, L. Liu, M. Luo, Z. Yan, Y. Xin, Robust infrared small object detection using hough line suppression and rank-hierarchy in complex backgrounds, Infrared Phys. Technol., 120 (2022), 103893. https://doi.org/10.1016/j.infrared.2021.103893 doi: 10.1016/j.infrared.2021.103893
![]() |
[7] | L. Zhong, Y. He, J. W. Zhang, Small object detection algorithm based on the fusion of context and semantic features, Comput. Appl., 42 (2022), 6. |
[8] |
L. Deng, J. Zhang, G. Xu, H. Zhu, Infrared small object detection via adaptive M-estimator ring top-hat transformation, Patt. Recognit., 112 (2021), 107729. https://doi.org/10.1016/j.patcog.2020.107729 doi: 10.1016/j.patcog.2020.107729
![]() |
[9] |
Y. Li, Z. Li, C. Zhang, Z. Luo, Y. Zhu, Z. Ding, Infrared maritime dim small object detection based on spatiotemporal cues and directional morphological filtering, Infrared Phys. Technol., 115 (2021), 103657. https://doi.org/10.1016/j.infrared.2021.103657 doi: 10.1016/j.infrared.2021.103657
![]() |
[10] |
Y. Lu, S. Huang, W. Zhao, Sparse representation based infrared small object detection via an online-learned double sparse background dictionary, Infrared Phys. Technol., 99 (2019), 14–27. https://doi.org/10.1016/j.infrared.2019.04.001 doi: 10.1016/j.infrared.2019.04.001
![]() |
[11] |
C. Chen, H. Li, Y. Wei, T. Xia, Y. Tang, A local contrast method for small infrared object detection, IEEE Trans. Geosci. Remote Sens., 52 (2013), 574–581. https://doi.org/10.1109/TGRS.2013.2242477 doi: 10.1109/TGRS.2013.2242477
![]() |
[12] |
L. Zhang, Z. Peng, Infrared small object detection based on partial sum of the tensor nuclear norm, Remote Sens., 11 (2019), 382. https://doi.org/10.3390/rs11040382 doi: 10.3390/rs11040382
![]() |
[13] |
T. Zhang, Z. Peng, H. Wu, Y. He, C. Li, C. Yang, Infrared small object detection via self-regularized weighted sparse model, Neurocomputing, 420 (2021), 124–148. https://doi.org/10.1016/j.neucom.2020.08.065 doi: 10.1016/j.neucom.2020.08.065
![]() |
[14] |
Y. Han, C. Deng, B. Zhao, D. Tao, State-aware anti-drift object tracking, IEEE Trans. Image Process., 5 (2019), 4075–4086. https://doi.org/10.1109/TIP.2019.2905984 doi: 10.1109/TIP.2019.2905984
![]() |
[15] |
I. V. Pustokhina, D. A. Pustokhin, T. Vaiyapuri, D. Gupta, S. Kumar, K. Shankar, An automated deep learning based anomaly detection in pedestrian walkways for vulnerable road users safety, Saf. Sci., 142 (2021), 105356. https://doi.org/10.1016/j.ssci.2021.105356 doi: 10.1016/j.ssci.2021.105356
![]() |
[16] | S. H. Xie, W. Z. Zhang, P. Cheng, YOLOv4 fire and smoke detection model with embedded channel attention, Chin. J. Liquid Crystal Displ., 36 (2021), 1445–1453. |
[17] | Y. Dai, Y. Wu, F. Zhou, K. Barnard, Asymmetric contextual modulation for infrared small object detection, in Proceedings of the IEEE/CVF Winter Conference on Applications of Computer Vision, (2021). https://doi.org/10.1109/WACV48630.2021.00099 |
[18] |
Y. Dai, Y. Wu, F. Zhou, K. Barnard, Attentional local contrast networks for infrared small object detection, IEEE Trans. Geosci. Remote Sens., 59 (2021), 9813–9823. https://doi.org/10.1109/TGRS.2020.3044958 doi: 10.1109/TGRS.2020.3044958
![]() |
[19] | H. Wang, L. Zhou, L. Wang, Miss detection vs. false alarm: Adversarial learning for small object segmentation in infrared images, in Proceedings of the IEEE/CVF International Conference on Computer Vision, (2019), 8509–8518. https://doi.org/10.1109/ICCV.2019.00860 |
[20] |
B. Li, C. Xiao, L. Wang, Y. Wang, Z. Lin, M. Li, Dense nested attention network for infrared small target detection, IEEE Trans. Image Process., 32 (2022), 1745–1758. https://doi.org/10.1109/TIP.2022.3199107 doi: 10.1109/TIP.2022.3199107
![]() |
[21] |
L. Huang, S. Dai, T. Huang, X. Huang, H. Wang, Infrared small object segmentation with multiscale feature representation, Infrared Phys. Technol., 116 (2021), 103755. https://doi.org/10.1016/j.infrared.2021.103755 doi: 10.1016/j.infrared.2021.103755
![]() |
[22] |
T. Zhang, S. Cao, T. Pu, Z. Peng, AGPCNet: Attention-guided pyramid context networks for infrared small target detection, IEEE Trans. Aerosp. Electron. Syst., 59 (2023), 4250–4261. https://doi.org/10.1109/TAES.2023.3238703 doi: 10.1109/TAES.2023.3238703
![]() |
[23] | M. Zhang, K. Yue, J. Zhang, Y. Li, X. Gao, Exploring feature compensation and cross-level correlation for infrared small target detection, in Proceedings of the 30th ACM International Conference on Multimedia, (2022), 1857–1865. https://doi.org/10.1145/3503161.3548264 |
[24] |
Y. Zhang, P. Tiňo, A. Leonardis, K. Tang, A survey on neural network interpretability, IEEE Trans. Emerg. Top Comput. Intell., (2021), 1–17. https://doi.org/10.1109/TETCI.2021.3100641 doi: 10.1109/TETCI.2021.3100641
![]() |
[25] | Y. Cao, J. Xu, S. Lin, F. Wei, H. Hu, Gcnet: Non-local networks meet squeeze-excitation networks and beyond, in Proceedings of the IEEE/CVF International Conference on Computer Vision Workshops, 2019. https://doi.org/10.48550/arXiv.1904.11492 |
[26] |
S. K. Ghosh, A. Ghosh, ENResNet: A novel residual neural network for chest X-ray enhancement based COVID-19 detection, Biomed. Signal Process. Control, 72 (2022), 103286. https://doi.org/10.1016/j.bspc.2021.103286 doi: 10.1016/j.bspc.2021.103286
![]() |
[27] |
W. Li, J. Li, J. Li, Z. Huang, D. Zhou, A lightweight multi-scale channel attention network for image super-resolution, Neurocomputing, 456 (2021), 327–337. https://doi.org/10.1016/j.neucom.2021.05.090 doi: 10.1016/j.neucom.2021.05.090
![]() |
[28] | X. Wang, K. Yu, S. Wu, J. Gu, Y. Liu, C. Dong, Esrgan: Enhanced super-resolution generative adversarial networks, in Proceedings of the European conference on computer vision (ECCV) workshops, 2018. https://doi.org/10.48550/arXiv.1809.00219 |
[29] | T. Y. Lin, P. Goyal, R. Girshick, K. He, P. Dollár, Focal loss for dense object detection, in Proceedings of the IEEE international conference on computer vision, 2017. https://doi.org/10.48550/arXiv.1708.02002 |
[30] |
G. Chen, W. Wang, X. Li, Designing and learning a lightweight network for infrared small target detection via dilated pyramid and semantic distillation, Infrared Phys. Technol., 131 (2023), 104671. https://doi.org/10.1016/j.infrared.2023.104671 doi: 10.1016/j.infrared.2023.104671
![]() |
[31] |
C. Xia, X. Li, L. Zhao, R. Shu, Infrared small object detection based on multiscale local contrast measure using local energy factor, IEEE Trans. Geosci. Remote Sens., 17 (2019), 157–161. https://doi.org/10.1109/LGRS.2019.2914432 doi: 10.1109/LGRS.2019.2914432
![]() |
[32] |
J. Han, S. Moradi, I. Faramarzi, C. Liu, H. Zhang, Q. Zhao, A local contrast method for infrared small-target detection utilizing a tri-layer window, IEEE Trans. Geosci. Remote Sens., 17 (2019), 1822–1826. https://doi.org/10.1109/LGRS.2019.2954578 doi: 10.1109/LGRS.2019.2954578
![]() |
[33] |
S. Huang, Y. Liu, Y. He, T. Zhang, Z. Peng, Structure-adaptive clutter suppression for infrared small object detection: Chain-growth filtering, Remote Sens., 12 (2020), 47. https://doi.org/10.3390/rs12010047 doi: 10.3390/rs12010047
![]() |
1. | Ben Lambert, David J. Gavaghan, Simon J. Tavener, A Monte Carlo method to estimate cell population heterogeneity from cell snapshot data, 2021, 511, 00225193, 110541, 10.1016/j.jtbi.2020.110541 | |
2. | Juan Shen, Xiao Tu, Yuanyuan Li, Mathematical Modeling Reveals Mechanisms of Cancer-Immune Interactions Underlying Hepatocellular Carcinoma Development, 2023, 11, 2227-7390, 4261, 10.3390/math11204261 | |
3. | Hao Jiang, Jingxin Liu, You Song, Jinzhi Lei, Quantitative Modeling of Stemness in Single-Cell RNA Sequencing Data: A Nonlinear One-Class Support Vector Machine Method, 2023, 1557-8666, 10.1089/cmb.2022.0484 | |
4. | Shenbageshwaran Rajendiran, Francisco Galdos, Carissa Anne Lee, Sidra Xu, Justin Harvell, Shireen Singh, Sean M. Wu, Elizabeth A. Lipke, Selen Cremaschi, 2024, 3, 9781777940324, 344, 10.69997/sct.152564 |
Genes | CDH1&ZEB1 | CDH1&KLF8 | ZEB1&KLF8 |
Spearman's correlations | –0.4664 | –0.4851 | 0.4879 |
Processes | Parameters | Values | Remarks |
Associate constant | p1 | 49.9611 | Fitted |
Associate constant | p2 | 0.5661 | Fitted |
Relative dissociation constant | p3 | 9.4701 | Fitted |
Relative degradation constant | p4 | 0.0951 | Fitted |
Associate constant | p5 | 2.3591 | Fitted |
Associate constant | p6 | 11.9229 | Fitted |
Relative degradation constant | p7 | 0.9982 | Fitted |
Hill coefficient | n1 | 2 | assumed |
Hill coefficient | n2 | 2 | assumed |
Hill coefficient | n3 | 2 | assumed |
Hill coefficient | n4 | 2 | assumed |
Intensity of impact | a | 1 | assumed |
Intensity of impact | r | 1 | assumed |
Validation | p1 | p2 | p3 | p4 | p5 | p6 | p7 | Error |
1 | 49.99 | 0.44 | 17.24 | 0.06 | 48.88 | 15.96 | 22.47 | 12.76 |
2 | 29.98 | 0.44 | 14.06 | 0.06 | 5.08 | 12.51 | 2.56 | 12.63 |
3 | 50.00 | 0.46 | 16.98 | 0.06 | 38.76 | 15.35 | 17.48 | 12.63 |
4 | 49.98 | 0.48 | 15.06 | 0.07 | 44.83 | 14.48 | 20.13 | 12.51 |
5 | 50.00 | 0.47 | 11.68 | 0.11 | 31.98 | 12.80 | 14.41 | 12.88 |
6 | 49.95 | 0.44 | 17.82 | 0.06 | 49.99 | 16.21 | 22.83 | 12.71 |
7 | 50.00 | 0.45 | 18.01 | 0.06 | 47.61 | 15.95 | 21.04 | 12.63 |
8 | 50.00 | 0.44 | 17.00 | 0.06 | 23.69 | 16.05 | 10.81 | 12.75 |
9 | 50.00 | 0.44 | 18.94 | 0.05 | 47.20 | 16.23 | 21.81 | 12.74 |
10 | 49.99 | 0.46 | 16.57 | 0.06 | 2.27 | 14.49 | 1.04 | 12.65 |
Average | 47.99 | 0.45 | 16.34 | 0.07 | 34.03 | 15.00 | 15.46 | 12.69 |
NO. | Reaction | Propensity function |
1 | → CDH1 | Ω3/(Ω2+p1×[ZEB1]2)+p2×Ω3/(Ω2+[KLF8]2) |
2 | CDH1 → | [CDH1] |
3 | → ZEB1 | Ω×[KLF8]2/(Ω2+p3×[KLF8]2) |
4 | ZEB1 → | p4×[ZEB1] |
5 | → KLF8 | p5×Ω3/(Ω2+p6×[CDH1]2) |
6 | KLF8 → | p7×[KLF8] |
Genes | CDH1&ZEB1 | CDH1&KLF8 | ZEB1&KLF8 |
Spearman's correlations | –0.4664 | –0.4851 | 0.4879 |
Processes | Parameters | Values | Remarks |
Associate constant | p1 | 49.9611 | Fitted |
Associate constant | p2 | 0.5661 | Fitted |
Relative dissociation constant | p3 | 9.4701 | Fitted |
Relative degradation constant | p4 | 0.0951 | Fitted |
Associate constant | p5 | 2.3591 | Fitted |
Associate constant | p6 | 11.9229 | Fitted |
Relative degradation constant | p7 | 0.9982 | Fitted |
Hill coefficient | n1 | 2 | assumed |
Hill coefficient | n2 | 2 | assumed |
Hill coefficient | n3 | 2 | assumed |
Hill coefficient | n4 | 2 | assumed |
Intensity of impact | a | 1 | assumed |
Intensity of impact | r | 1 | assumed |
Validation | p1 | p2 | p3 | p4 | p5 | p6 | p7 | Error |
1 | 49.99 | 0.44 | 17.24 | 0.06 | 48.88 | 15.96 | 22.47 | 12.76 |
2 | 29.98 | 0.44 | 14.06 | 0.06 | 5.08 | 12.51 | 2.56 | 12.63 |
3 | 50.00 | 0.46 | 16.98 | 0.06 | 38.76 | 15.35 | 17.48 | 12.63 |
4 | 49.98 | 0.48 | 15.06 | 0.07 | 44.83 | 14.48 | 20.13 | 12.51 |
5 | 50.00 | 0.47 | 11.68 | 0.11 | 31.98 | 12.80 | 14.41 | 12.88 |
6 | 49.95 | 0.44 | 17.82 | 0.06 | 49.99 | 16.21 | 22.83 | 12.71 |
7 | 50.00 | 0.45 | 18.01 | 0.06 | 47.61 | 15.95 | 21.04 | 12.63 |
8 | 50.00 | 0.44 | 17.00 | 0.06 | 23.69 | 16.05 | 10.81 | 12.75 |
9 | 50.00 | 0.44 | 18.94 | 0.05 | 47.20 | 16.23 | 21.81 | 12.74 |
10 | 49.99 | 0.46 | 16.57 | 0.06 | 2.27 | 14.49 | 1.04 | 12.65 |
Average | 47.99 | 0.45 | 16.34 | 0.07 | 34.03 | 15.00 | 15.46 | 12.69 |
NO. | Reaction | Propensity function |
1 | → CDH1 | Ω3/(Ω2+p1×[ZEB1]2)+p2×Ω3/(Ω2+[KLF8]2) |
2 | CDH1 → | [CDH1] |
3 | → ZEB1 | Ω×[KLF8]2/(Ω2+p3×[KLF8]2) |
4 | ZEB1 → | p4×[ZEB1] |
5 | → KLF8 | p5×Ω3/(Ω2+p6×[CDH1]2) |
6 | KLF8 → | p7×[KLF8] |