Effects of Small Random Perturbations in the Extended Glass–Kauffman Model of Gene Regulatory Networks

Arcady Ponosov,Irina Shlykova,Ramazan I. Kadiev
DOI: https://doi.org/10.3390/math12081223
IF: 2.4
2024-04-19
Mathematics
Abstract:A mathematical justification of some basic structural properties of stochastically perturbed gene regulatory networks, including those with autoregulation and delay, is offered in this paper. By using the theory of stochastic differential equations, it is, in particular, shown how to control the asymptotic behavior of the diffusion terms in order to not destroy certain qualitative features of the networks, for instance, their sliding modes. The results also confirm that the level of randomness is gradually reduced if the gene activation times become much smaller than the time of interaction of genes. Finally, the suggested analysis explains why the deterministic numerical schemes based on replacing smooth, steep response functions by the simpler yet discontinuous Heaviside function, the well-known simplification algorithm, are robust with respect to uncertainties in data. The main technical difficulties of the analysis are handled by applying the uniform version of the stochastic Tikhonov theorem in singular perturbation analysis suggested by Yu. Kabanov and S. Pergamentshchikov.
mathematics
What problem does this paper attempt to address?
### Problems the paper attempts to solve This paper aims to explore the impact after introducing small random perturbations in the extended Glass - Kauffman gene regulatory network model. Specifically, the paper focuses on the following aspects: 1. **The influence of random perturbations on the structural characteristics of the gene regulatory network**: - Using the theory of stochastic differential equations, study how to control the asymptotic behavior of the diffusion term to maintain certain qualitative features of the network, such as sliding modes. 2. **The change in the level of randomness**: - Study the phenomenon that the level of randomness gradually decreases when the gene activation time is much smaller than the gene interaction time. 3. **The robustness of deterministic numerical methods**: - Explain why the simplified algorithm based on replacing the smooth and steep response function with the discontinuous Heaviside function is robust in the face of data uncertainty. 4. **Sliding - mode analysis**: - Provide a method to calculate the sliding - mode trajectories and prove the stability of these sliding - modes under small perturbations, which helps to understand the regulation efficiency of gene expression and thus contributes to the overall function of the network. ### Main technical difficulties - **Application of the stochastic Tikhonov theorem**: - Deal with the main technical problems by applying the uniform version of the Tikhonov theorem in stochastic singular perturbation analysis proposed by Yu. Kabanov and S. Pergamentshchikov. - **Handling of the delay effect**: - Convert the system with delay into a system of ordinary differential equations by the modified linear chain trick for analysis. ### Key assumptions - **A single variable is close to the threshold**: - Assume that at any time, at most only one variable \(x_i\) is close to its threshold \(\theta_i\), while other variables are far from their thresholds. ### Conclusions Through the above analysis, the paper verifies that under small random perturbations, the deterministic gene regulatory network model dominates when the response function is steep enough. In addition, the paper also provides a method for dealing with the stochastic gene regulatory network with delay effects, ensuring the robustness of the deterministic numerical method under small random perturbations.