Source author record

Jinzhi Lei

Jinzhi Lei appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
17topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2022arXiv

Power law dependence in a random differential equation

This paper studies a random differential equation with random switch perturbations. We explore how the maximum displacement from the equilibrium state depends on the statistical properties of time series of {the} random switches. We show a power law dependence between the upper bound of displacement and the frequency of random perturbation switches, and the slope of power law dependence is dependent on the specific distribution of the intervals between switching times. This result {suggests} a quantitative connection between frequency modulation and amplitude modulation under random perturbations.

preprint2020arXiv

DNA methylation heterogeneity induced by collaborations between enhancers

During mammalian embryo development, reprogramming of DNA methylation plays important roles in the erasure of parental epigenetic memory and the establishment of naïve pluripogent cells. Multiple enzymes that regulate the processes of methylation and demethylation work together to shape the pattern of genome-scale DNA methylation and guid the process of cell differentiation. Recent availability of methylome information from single-cell whole genome bisulfite sequencing (scBS-seq) provides an opportunity to study DNA methylation dynamics in the whole genome in individual cells, which reveal the heterogeneous methylation distributions of enhancers in embryo stem cells (ESCs). In this study, we developed a computational model of enhancer methylation inheritance to study the dynamics of genome-scale DNA methylation reprogramming during exit from pluripotency. The model enables us to track genome-scale DNA methylation reprogramming at single-cell level during the embryo development process, and reproduce the DNA methylation heterogeneity reported by scBS-seq. Model simulations show that DNA methylation heterogeneity is an intrinsic property driven by cell division along the development process, and the collaboration between neighboring enhancers is required for heterogeneous methylation. Our study suggest that the mechanism of genome-scale oscillation proposed by Rulands et al. (2018) might not necessary to the DNA methylation during exit from pluripotency.

preprint2020arXiv

Single-cell entropy to quantify the cellular transcription from single-cell RNA-seq data

We present the use of single-cell entropy (scEntropy) to measure the order of the cellular transcriptome profile from single-cell RNA-seq data, which leads to a method of unsupervised cell type classification through scEntropy followed by the Gaussian mixture model (scEGMM). scEntropy is straightforward in defining an intrinsic transcriptional state of a cell. scEGMM is a coherent method of cell type classification that includes no parameters and no clustering; however, it is comparable to existing machine learning-based methods in benchmarking studies and facilitates biological interpretation.

preprint2019arXiv

Evolutionary dynamics of cancer: from epigenetic regulation to cell population dynamics -- mathematical model framework, applications, and open problems

Predictive modeling of the evolutionary dynamics of cancer is a challenge issue in computational cancer biology. In this paper, we propose a general mathematical model framework for the evolutionary dynamics of cancer with plasticity and heterogeneity in cancer cells. Cancer is a group of diseases involving abnormal cell growth, during which abnormal regulations in stem cell regeneration are essential for the dynamics of cancer development. In general, the dynamics of stem cell regeneration can be simplified as a $\mathrm{G_0}$ phase cell cycle model, which lead to a delay differentiation equation. When cell heterogeneity and plasticity are considered, we establish a differential-integral equation based on the random transition of epigenetic states of stem cells during cell division. The proposed model highlights cell heterogeneity and plasticity, and connects the heterogeneity with cell-to-cell variance in cellular behaviors, e.g. proliferation, apoptosis, and differentiation/senescence, and can be extended to include gene mutation-induced tumor development. Hybrid computations models are developed based on the mathematical model framework, and are applied to the process of inflammation-induced tumorigenesis and tumor relapse after CAR-T therapy. Finally, we give rise to several mathematical problems related to the proposed differential-integral equation. Answers to these problems are crucial for the understanding of the evolutionary dynamics of cancer.

preprint2015arXiv

aa-tRNA competition is crucial for the effective translation efficiency

Translation is a central biological process by which proteins are synthesized from genetic information contained within mRNAs. Here we study the kinetics of translation at molecular level through a stochastic simulation model. The model explicitly include RNA sequences, ribosome dynamics, tRNA pool and biochemical reactions in the translation elongation. The results show that the translation efficiency is mainly limited by the available ribosome number, translation initiation and the translation elongation time. The elongation time is log-normal distribution with mean and variance determined by both the codon saturation and the process of aa-tRNA selection at each codon binding. Moreover, our simulations show that the translation accuracy exponentially decreases with the sequence length. These results suggest that aa-tRNA competition is crucial for both translation elongation, translation efficiency and the accuracy, which in turn determined the effective protein production rate of correct proteins. Our results improve the dynamical equation of protein production with a delay differential equation which is dependent on sequence informations through both the effective production rate and the distribution of elongation time.

preprint2015arXiv

Bifurcation analysis and potential landscape of the p53-Mdm2 oscillator regulated by the co-activator PDCD5

Dynamics of p53 is known to play important roles in the regulation of cell fate decisions in response to various stresses, and PDCD5 functions as a co-activator of p53 to modulate the p53 dynamics. In the present paper, we investigate how p53 dynamics are modulated by PDCD5 during the DNA damage response using methods of bifurcation analysis and potential landscape. Our results reveal that p53 activities can display rich dynamics under different PDCD5 levels, including monostability, bistability with two stable steady states, oscillations, and co-existence of a stable steady state and an oscillatory state. Physical properties of the p53 oscillations are further shown by the potential landscape, in which the potential force attracts the system state to the limit cycle attractor, and the curl flux force drives the coherent oscillation along the cyclic. We also investigate the effect of PDCD5 efficiency on inducing the p53 oscillations. We show that Hopf bifurcation is induced by increasing the PDCD5 efficiency, and the system dynamics show clear transition features in both barrier height and energy dissipation when the efficiency is close to the bifurcation point. This study provides a global picture of how PDCD5 regulates p53 dynamics via the interaction with the p53-Mdm2 oscillator and can be helpful in understanding the complicate p53 dynamics in a more complete p53 pathway.

preprint2015arXiv

PDCD5 interacts with p53 and functions as a regulator of p53 dynamics in the DNA damage response

The tumor suppressor p53 plays a central role in cell fate decisions after DNA damage. Programmed Cell Death 5 (PDCD5) is known to interact with the p53 pathway to promote cell apoptosis. Recombinant human PDCD5 can significantly sensitize different cancers to chemotherapies. In the present paper, we construct a computational model that includes PDCD5 interactions in the p53 signaling network and study the effects of PDCD5 on p53-mediated cell fate decisions during the DNA damage response. Our results revealed that PDCD5 functions as a co-activator of p53 that regulates p53-dependent cell fate decisions via the mediation of p53 dynamics. The effects of PDCD5 are dose-dependent such that p53 can display either sustained or pulsed dynamics at different PDCD5 levels. Moreover, PDCD5 regulates caspase-3 activation via two mechanisms during the two phases of sustained and pulsed p53 dynamics. This study provides insights regarding how PDCD5 functions as a regulator of the p53 pathway and might be helpful for increasing our understanding of the molecular mechanisms by which PDCD5 can be used to treat cancers.

preprint2013arXiv

Adiabatic reduction of a model of stochastic gene expression with jump Markov process

This paper considers adiabatic reduction in a model of stochastic gene expression with bursting transcription considered as a jump Markov process. In this model, the process of gene expression with auto-regulation is described by fast/slow dynamics. The production of mRNA is assumed to follow a compound Poisson process occurring at a rate depending on protein levels (the phenomena called bursting in molecular biology) and the production of protein is a linear function of mRNA numbers. When the dynamics of mRNA is assumed to be a fast process (due to faster mRNA degradation than that of protein) we prove that, with appropriate scalings in the burst rate, jump size or translational rate, the bursting phenomena can be transmitted to the slow variable. We show that, depending on the scaling, the reduced equation is either a stochastic differential equation with a jump Poisson process or a deterministic ordinary differential equation. These results are significant because adiabatic reduction techniques seem to have not been rigorously justified for a stochastic differential system containing a jump Markov process. We expect that the results can be generalized to adiabatic methods in more general stochastic hybrid systems.

preprint2012arXiv

Moment Boundedness of Linear Stochastic Delay Differential Equation with Distributed Delay

This paper studies the moment boundedness of solutions of linear stochastic delay differential equations with distributed delay. For a linear stochastic delay differential equation, the first moment stability is known to be identical to that of the corresponding deterministic delay differential equation. However, boundedness of the second moment is complicated and depends on the stochastic terms. In this paper, the characteristic function of the equation is obtained through techniques of Laplace transform. From the characteristic equation, sufficient conditions for the second moment to be bounded or unbounded are proposed.

preprint2012arXiv

On positive solutions and the Omega limit set for a class of delay differential equations

This paper studies the positive solutions of a class of delay differential equations with two delays. These equations originate from the modeling of hematopoietic cell populations. We give a sufficient condition on the initial function for $t\leq 0$ such that the solution is positive for all time $t>0$. The condition is "optimal". We also discuss the long time behavior of these positive solutions through a dynamical system on the space of continuous functions. We give a characteristic description of the $ω$ limit set of this dynamical system, which can provide informations about the long time behavior of positive solutions of the delay differential equation.

preprint2012arXiv

Second Moment Boundedness of Linear Stochastic Delay Differential Equations

This paper studies the second moment boundedness of solutions of linear stochastic delay differential equations. First, we give a framework, for general $\mathrm{N}$-dimensional linear stochastic differential equations with a single discrete delay, of calculating the characteristic function for the second moment boundedness. Next, we apply the proposed framework to a special case of a type of 2-dimensional equation that the stochastic terms are decoupled. For the 2-dimensional equation, we obtain the characteristic function explicitly given by equation coefficients, the characteristic function gives sufficient conditions for the second moment to be bounded or unbounded.

preprint2011arXiv

Deterministic Brownian motion generated from differential delay equations

This paper addresses the question of how Brownian-like motion can arise from the solution of a deterministic differential delay equation. To study this we analytically study the bifurcation properties of an apparently simple differential delay equation and then numerically investigate the probabilistic properties of chaotic solutions of the same equation. Our results show that solutions of the deterministic equation with randomly selected initial conditions display a Gaussian-like density for long time, but the densities are supported on an interval of finite measure. Using these chaotic solutions as velocities, we are able to produce Brownian-like motions, which show statistical properties akin to those of a classical Brownian motion over both short and long time scales. Several conjectures are formulated for the probabilistic properties of the solution of the differential delay equation. Numerical studies suggest that these conjectures could be "universal" for similar types of "chaotic" dynamics, but we have been unable to prove this.

preprint2011arXiv

Nonlinear differential Galois theory

Differential Galois theory has played important roles in the theory of integrability of linear differential equation. In this paper we will extend the theory to nonlinear case and study the integrability of the first order nonlinear differential equation. We will define for the differential equation the differential Galois group, will study the structure of the group, and will prove the equivalent between the existence of the Liouvillian first integral and the solvability of the corresponding differential Galois group.

preprint2011arXiv

On the classification of polynomial differential operators

This paper gives a classification of first order polynomial differential operators of form $\mathscr{X} = X_1(x_1,x_2)δ_1 + X_2(x_1,x_2)δ_2$, $(δ_i = \partial/\partial x_i)$. The classification is given through the order of an operator that is defined in this paper. Let $X=\mathscr{X}y$ to be the differential polynomial associated with $\mathscr{X}$, the order of $\mathscr{X}$, $\mathrm{ord}(\mathscr{X})$, is defined as the order of a differential ideal $Λ$ of differential polynomials that is a nontrivial expansion of the ideal $\{X\}$ and with the lowest order. In this paper, we prove that there are only four possible values for the order of a differential operator, 0, 1, 2, 3, or $\infty$. Furthermore, when the order is finite, the expansion $Λ$ is generated by $X$ and a differential polynomial $A$, which can be obtained through a rational solution of a partial differential equation that is given explicitly in this paper. When the order is infinite, the expansion $Λ$ is just the unit ideal. In additional, if, and only if, the order of $\mathscr{X}$ is 0, 1, or 2, the polynomial differential equation associating with $\mathscr{X}$ has Liouvillian first integrals. Examples for each class of differential operators are given at the end of this paper.

preprint2011arXiv

Stochastic Modeling in Systems Biology

Many cellular behaviors are regulated by gene regulation networks, kinetics of which is one of the main subjects in the study of systems biology. Because of the low number molecules in these reacting systems, stochastic effects are significant. In recent years, stochasticity in modeling the kinetics of gene regulation networks have been drawing the attention of many researchers. This paper is a self contained review trying to provide an overview of stochastic modeling. I will introduce the derivation of the main equations in modeling the biochemical systems with intrinsic noise (chemical master equation, Fokker-Plan equation, reaction rate equation, chemical Langevin equation), and will discuss the relations between these formulations. The mathematical formulations for systems with fluctuations in kinetic parameters are also discussed. Finally, I will introduce the exact stochastic simulation algorithm and the approximate explicit tau-leaping method for making numerical simulations.

preprint2010arXiv

Elastic energy of proteins and the stages of protein folding

We propose a universal elastic energy for proteins, which depends only on the radius of gyration $R_{g}$ and the residue number $N$. It is constructed using physical arguments based on the hydrophobic effect and hydrogen bonding. Adjustable parameters are fitted to data from the computer simulation of the folding of a set of proteins using the CSAW (conditioned self-avoiding walk) model. The elastic energy gives rise to scaling relations of the form $R_{g}\sim N^ν$ in different regions. It shows three folding stages characterized by the progression with exponents $ν= 3/5, 3/7, 2/5$, which we identify as the unfolded stage, pre-globule, and molten globule, respectively. The pre-globule goes over to the molten globule via a break in behavior akin to a first-order phase transition, which is initiated by a sudden acceleration of hydrogen bonding.

preprint2010arXiv

Protein Folding: A Perspective From Statistical Physics

In this paper, we introduce an approach to the protein folding problem from the point of view of statistical physics. Protein folding is a stochastic process by which a polypeptide folds into its characteristic and functional 3D structure from random coil. The process involves an intricate interplay between global geometry and local structure, and each protein seems to present special problems. We introduce CSAW (conditioned self-avoiding walk), a model of protein folding that combines the features of self-avoiding walk (SAW) and the Monte Carlo method. In this model, the unfolded protein chain is treated as a random coil described by SAW. Folding is induced by hydrophobic forces and other interactions, such as hydrogen bonding, which can be taken into account by imposing conditions on SAW. Conceptually, the mathematical basis is a generalized Langevin equation. To illustrate the flexibility and capabilities of the model, we consider several examples, including helix formation, elastic properties, and the transition in the folding of myoglobin. From the CSAW simulation and physical arguments, we find a universal elastic energy for proteins, which depends only on the radius of gyration $R_{g}$ and the residue number $N$. The elastic energy gives rise to scaling laws $R_{g}\sim N^ν$ in different regions with exponents $ν=3/5,3/7,2/5$, consistent with the observed unfolded stage, pre-globule, and molten globule, respectively. These results indicate that CSAW can serve as a theoretical laboratory to study universal principles in protein folding.

preprint2007arXiv

Dynamics of Alpha-Helix Formation in the CSAW Model

We study the folding dynamics of polyalanine (Ala$_{20}$), a protein fragment with 20 residues whose native state is a single alpha helix. We use the CSAW model (conditioned self-avoiding walk), which treats the protein molecule as a chain in Brownian motion, with interactions that include hydrophobic forces and internal hydrogen bonding. We find that large scale structures form before small scale structures, and obtain the relevant relaxation times. We find that helix nucleation occurs at two separate points on the protein chain. The evolution of small and large scale structures involve different mechanisms. While the former can be describe by rate equations governing the growth of helical content, the latter is akin to the relaxation of an elastic solid.