Showing 1–26 of 26 results for author: Lei, J

AutoRG-Brain: Grounded Report Generation for Brain MRI

Authors: Jiayu Lei, Xiaoman Zhang, Chaoyi Wu, Lisong Dai, Ya Zhang, Yanyong Zhang, Yanfeng Wang, Weidi Xie, Yuehua Li

Abstract: Radiologists are tasked with interpreting a large number of images in a daily base, with the responsibility of generating corresponding reports. This demanding workload elevates the risk of human error, potentially leading to treatment delays, increased healthcare costs, revenue loss, and operational inefficiencies. To address these challenges, we initiate a series of work on grounded Automatic Re… ▽ More Radiologists are tasked with interpreting a large number of images in a daily base, with the responsibility of generating corresponding reports. This demanding workload elevates the risk of human error, potentially leading to treatment delays, increased healthcare costs, revenue loss, and operational inefficiencies. To address these challenges, we initiate a series of work on grounded Automatic Report Generation (AutoRG), starting from the brain MRI interpretation system, which supports the delineation of brain structures, the localization of anomalies, and the generation of well-organized findings. We make contributions from the following aspects, first, on dataset construction, we release a comprehensive dataset encompassing segmentation masks of anomaly regions and manually authored reports, termed as RadGenome-Brain MRI. This data resource is intended to catalyze ongoing research and development in the field of AI-assisted report generation systems. Second, on system design, we propose AutoRG-Brain, the first brain MRI report generation system with pixel-level grounded visual clues. Third, for evaluation, we conduct quantitative assessments and human evaluations of brain structure segmentation, anomaly localization, and report generation tasks to provide evidence of its reliability and accuracy. This system has been integrated into real clinical scenarios, where radiologists were instructed to write reports based on our generated findings and anomaly segmentation masks. The results demonstrate that our system enhances the report-writing skills of junior doctors, aligning their performance more closely with senior doctors, thereby boosting overall productivity. △ Less

Submitted 29 July, 2024; v1 submitted 23 July, 2024; originally announced July 2024.

Exploring Biomarker Relationships in Both Type 1 and Type 2 Diabetes Mellitus Through a Bayesian Network Analysis Approach

Authors: Yuyang Sun, Jingyu Lei, Panagiotis Kosmas

Abstract: Understanding the complex relationships of biomarkers in diabetes is pivotal for advancing treatment strategies, a pressing need in diabetes research. This study applies Bayesian network structure learning to analyze the Shanghai Type 1 and Type 2 diabetes mellitus datasets, revealing complex relationships among key diabetes-related biomarkers. The constructed Bayesian network presented notable pr… ▽ More Understanding the complex relationships of biomarkers in diabetes is pivotal for advancing treatment strategies, a pressing need in diabetes research. This study applies Bayesian network structure learning to analyze the Shanghai Type 1 and Type 2 diabetes mellitus datasets, revealing complex relationships among key diabetes-related biomarkers. The constructed Bayesian network presented notable predictive accuracy, particularly for Type 2 diabetes mellitus, with root mean squared error (RMSE) of 18.23 mg/dL, as validated through leave-one-domain experiments and Clarke error grid analysis. This study not only elucidates the intricate dynamics of diabetes through a deeper understanding of biomarker interplay but also underscores the significant potential of integrating data-driven and knowledge-driven methodologies in the realm of personalized diabetes management. Such an approach paves the way for more custom and effective treatment strategies, marking a notable advancement in the field. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Comments: Paper is accepted by EMBC 2024

A computational scheme connecting gene regulatory network dynamics with heterogeneous stem cell regeneration

Authors: Yakun Li, Xiyin Liang, Jinzhi Lei

Abstract: Stem cell regeneration is a vital biological process in self-renewing tissues, governing development and tissue homeostasis. Gene regulatory network dynamics are pivotal in controlling stem cell regeneration and cell type transitions. However, integrating the quantitative dynamics of gene regulatory networks at the single-cell level with stem cell regeneration at the population level poses signifi… ▽ More Stem cell regeneration is a vital biological process in self-renewing tissues, governing development and tissue homeostasis. Gene regulatory network dynamics are pivotal in controlling stem cell regeneration and cell type transitions. However, integrating the quantitative dynamics of gene regulatory networks at the single-cell level with stem cell regeneration at the population level poses significant challenges. This study presents a computational framework connecting gene regulatory network dynamics with stem cell regeneration through a data-driven formulation of the inheritance function. The inheritance function captures epigenetic state transitions during cell division in heterogeneous stem cell populations. Our scheme allows the derivation of the inheritance function based on a hybrid model of cross-cell-cycle gene regulation network dynamics. The proposed scheme enables us to derive the inheritance function based on the hybrid model of cross-cell-cycle gene regulation network dynamics. By explicitly incorporating gene regulatory network structure, it replicates cross-cell-cycling gene regulation dynamics through individual-cell-based modeling. The numerical scheme holds the potential for extension to diverse gene regulatory networks, facilitating a deeper understanding of the connection between gene regulation dynamics and stem cell regeneration. △ Less

Submitted 17 April, 2024; originally announced April 2024.

Comments: 27 pages, 9 figures

Modelling infectious disease transmission dynamics in conference environments: An individual-based approach

Authors: Xue Liu, Yue Deng, Jingying Huang, Yuhong Zhang, Jinzhi Lei

Abstract: The global public health landscape is perpetually challenged by the looming threat of infectious diseases. Central to addressing this concern is the imperative to prevent and manage disease transmission during pandemics, particularly in unique settings. This study addresses the transmission dynamics of infectious diseases within conference venues, presenting a computational model designed to simul… ▽ More The global public health landscape is perpetually challenged by the looming threat of infectious diseases. Central to addressing this concern is the imperative to prevent and manage disease transmission during pandemics, particularly in unique settings. This study addresses the transmission dynamics of infectious diseases within conference venues, presenting a computational model designed to simulate transmission processes within a condensed timeframe (one day), beginning with sporadic cases. Our model intricately captures the activities of individual attendees within the conference venue, encompassing meetings, rest intervals, and meal breaks. While meetings entail proximity seating, rest and lunch periods allow attendees to interact with diverse individuals. Moreover, the restroom environment poses an additional avenue for potential infection transmission. Employing an individual-based model, we meticulously replicated the transmission dynamics of infectious diseases, with a specific emphasis on close-contact interactions between infected and susceptible individuals. Through comprehensive analysis of model simulations, we elucidated the intricacies of disease transmission dynamics within conference settings and assessed the efficacy of control strategies to curb disease dissemination. Ultimately, our study proffers a numerical framework for assessing the risk of infectious disease transmission during short-duration conferences, furnishing conference organizers with valuable insights to inform the implementation of targeted prevention and control measures. △ Less

Submitted 17 April, 2024; originally announced April 2024.

Comments: 25 pages; 8 figures

scRDiT: Generating single-cell RNA-seq data by diffusion transformers and accelerating sampling

Authors: Shengze Dong, Zhuorui Cui, Ding Liu, Jinzhi Lei

Abstract: Motivation: Single-cell RNA sequencing (scRNA-seq) is a groundbreaking technology extensively utilized in biological research, facilitating the examination of gene expression at the individual cell level within a given tissue sample. While numerous tools have been developed for scRNA-seq data analysis, the challenge persists in capturing the distinct features of such data and replicating virtual d… ▽ More Motivation: Single-cell RNA sequencing (scRNA-seq) is a groundbreaking technology extensively utilized in biological research, facilitating the examination of gene expression at the individual cell level within a given tissue sample. While numerous tools have been developed for scRNA-seq data analysis, the challenge persists in capturing the distinct features of such data and replicating virtual datasets that share analogous statistical properties. Results: Our study introduces a generative approach termed scRNA-seq Diffusion Transformer (scRDiT). This method generates virtual scRNA-seq data by leveraging a real dataset. The method is a neural network constructed based on Denoising Diffusion Probabilistic Models (DDPMs) and Diffusion Transformers (DiTs). This involves subjecting Gaussian noises to the real dataset through iterative noise-adding steps and ultimately restoring the noises to form scRNA-seq samples. This scheme allows us to learn data features from actual scRNA-seq samples during model training. Our experiments, conducted on two distinct scRNA-seq datasets, demonstrate superior performance. Additionally, the model sampling process is expedited by incorporating Denoising Diffusion Implicit Models (DDIM). scRDiT presents a unified methodology empowering users to train neural network models with their unique scRNA-seq datasets, enabling the generation of numerous high-quality scRNA-seq samples. Availability and implementation: https://github.com/DongShengze/scRDiT △ Less

Submitted 9 April, 2024; originally announced April 2024.

Comments: 11 pages, 4 figures,

Mathematical modeling of heterogeneous stem cell regeneration: from cell division to Waddington's epigenetic landscape

Authors: Jinzhi Lei

Abstract: Stem cell regeneration is a crucial biological process for most self-renewing tissues during the development and maintenance of tissue homeostasis. In developing the mathematical models of stem cell regeneration and tissue development, cell division is the core process connecting different scale biological processes and leading to changes in cell population number and the epigenetic state of cells… ▽ More Stem cell regeneration is a crucial biological process for most self-renewing tissues during the development and maintenance of tissue homeostasis. In developing the mathematical models of stem cell regeneration and tissue development, cell division is the core process connecting different scale biological processes and leading to changes in cell population number and the epigenetic state of cells. This chapter focuses on the primary strategies for modeling cell division in biological systems. The Lagrange coordinate modeling approach considers gene network dynamics within each cell and random changes in cell states and model parameters during cell division. In contrast, the Euler coordinate modeling approach formulates the evolution of cell population numbers with the same epigenetic state via a differential-integral equation. These strategies focus on different scale dynamics, respectively, and result in two methods of modeling Waddington's epigenetic landscape: the Fokker-Planck equation and the differential-integral equation approaches. The differential-integral equation approach formulates the evolution of cell population density based on simple assumptions in cell proliferation, apoptosis, differentiation, and epigenetic state transitions during cell division. Moreover, machine learning methods can establish low-dimensional macroscopic measurements of a cell based on single-cell RNA sequencing data. The low dimensional measurements can quantify the epigenetic state of cells and become connections between static single-cell RNA sequencing data with dynamic equations for tissue development processes. The differential-integral equation presented in this chapter provides a reasonable approach to understanding the complex biological processes of tissue development and tumor progression. △ Less

Submitted 16 January, 2024; v1 submitted 14 September, 2023; originally announced September 2023.

Comments: 46pages, 1 figures

Dynamics of cell-type transition mediated by epigenetic modifications

Authors: Rongsheng Huang, Qiaojun Situ, Jinzhi Lei

Abstract: Maintaining tissue homeostasis requires appropriate regulation of stem cell differentiation. The Waddington landscape posits that gene circuits in a cell form a potential landscape of different cell types, wherein cells follow attractors of the probability landscape to develop into distinct cell types. However, how adult stem cells achieve a delicate balance between self-renewal and differentiatio… ▽ More Maintaining tissue homeostasis requires appropriate regulation of stem cell differentiation. The Waddington landscape posits that gene circuits in a cell form a potential landscape of different cell types, wherein cells follow attractors of the probability landscape to develop into distinct cell types. However, how adult stem cells achieve a delicate balance between self-renewal and differentiation remains unclear. We propose that random inheritance of epigenetic states plays a pivotal role in stem cell differentiation and present a hybrid model of stem cell differentiation induced by epigenetic modifications. Our comprehensive model integrates gene regulation networks, epigenetic state inheritance, and cell regeneration, encompassing multi-scale dynamics ranging from transcription regulation to cell population. Through model simulations, we demonstrate that random inheritance of epigenetic states during cell divisions can spontaneously induce cell differentiation, dedifferentiation, and transdifferentiation. Furthermore, we investigate the influences of interfering with epigenetic modifications and introducing additional transcription factors on the probabilities of dedifferentiation and transdifferentiation, revealing the underlying mechanism of cell reprogramming. This \textit{in silico} model provides valuable insights into the intricate mechanism governing stem cell differentiation and cell reprogramming and offers a promising path to enhance the field of regenerative medicine. △ Less

Submitted 13 September, 2023; originally announced September 2023.

Comments: 34 pages, 12 figures

Activate index: an integrated index to reveal disrupted brain network organizations of major depressive disorder patients

Authors: Yu Fu, Yanyan Huang, Meng Niu, Le Xue, Shunjie Dong, Shunlin Guo, Junqiang Lei, Cheng Zhuo

Abstract: Altered functional brain networks have been a typical manifestation that distinguishes major depressive disorder (MDD) patients from healthy control (HC) subjects in functional magnetic resonance imaging (fMRI) studies. Recently, rich club and diverse club metrics have been proposed for network or network neuroscience analyses. The rich club defines a set of nodes that tend to be the hubs of speci… ▽ More Altered functional brain networks have been a typical manifestation that distinguishes major depressive disorder (MDD) patients from healthy control (HC) subjects in functional magnetic resonance imaging (fMRI) studies. Recently, rich club and diverse club metrics have been proposed for network or network neuroscience analyses. The rich club defines a set of nodes that tend to be the hubs of specific communities, and the diverse club defines the nodes that span more communities and have edges diversely distributed across different communities. Considering the heterogeneity of rich clubs and diverse clubs, combining them and on the basis to derive a novel indicator may reveal new evidence of brain functional integration and separation, which might provide new insights into MDD. This study for the first time discussed the differences between MDD and HC using both rich club and diverse club metrics and found the complementarity of them in analyzing brain networks. Besides, a novel index, termed "active index", has been proposed in this study. The active index defines a group of nodes that tend to be diversely distributed across communities while avoiding being a hub of a community. Experimental results demonstrate the superiority of active index in analyzing MDD brain mechanisms. △ Less

Submitted 14 February, 2022; originally announced February 2022.

ChemRL-GEM: Geometry Enhanced Molecular Representation Learning for Property Prediction

Authors: Xiaomin Fang, Lihang Liu, Jieqiong Lei, Donglong He, Shanzhuo Zhang, Jingbo Zhou, Fan Wang, Hua Wu, Haifeng Wang

Abstract: Effective molecular representation learning is of great importance to facilitate molecular property prediction, which is a fundamental task for the drug and material industry. Recent advances in graph neural networks (GNNs) have shown great promise in applying GNNs for molecular representation learning. Moreover, a few recent studies have also demonstrated successful applications of self-supervise… ▽ More Effective molecular representation learning is of great importance to facilitate molecular property prediction, which is a fundamental task for the drug and material industry. Recent advances in graph neural networks (GNNs) have shown great promise in applying GNNs for molecular representation learning. Moreover, a few recent studies have also demonstrated successful applications of self-supervised learning methods to pre-train the GNNs to overcome the problem of insufficient labeled molecules. However, existing GNNs and pre-training strategies usually treat molecules as topological graph data without fully utilizing the molecular geometry information. Whereas, the three-dimensional (3D) spatial structure of a molecule, a.k.a molecular geometry, is one of the most critical factors for determining molecular physical, chemical, and biological properties. To this end, we propose a novel Geometry Enhanced Molecular representation learning method (GEM) for Chemical Representation Learning (ChemRL). At first, we design a geometry-based GNN architecture that simultaneously models atoms, bonds, and bond angles in a molecule. To be specific, we devised double graphs for a molecule: The first one encodes the atom-bond relations; The second one encodes bond-angle relations. Moreover, on top of the devised GNN architecture, we propose several novel geometry-level self-supervised learning strategies to learn spatial knowledge by utilizing the local and global molecular 3D structures. We compare ChemRL-GEM with various state-of-the-art (SOTA) baselines on different molecular benchmarks and exhibit that ChemRL-GEM can significantly outperform all baselines in both regression and classification tasks. For example, the experimental results show an overall improvement of 8.8% on average compared to SOTA baselines on the regression tasks, demonstrating the superiority of the proposed method. △ Less

Submitted 22 February, 2022; v1 submitted 10 June, 2021; originally announced June 2021.

Comments: Nature Machine Intelligence, 2022

Journal ref: Nature Machine Intelligence, 2022

Object Based Attention Through Internal Gating

Authors: Jordan Lei, Ari S. Benjamin, Konrad P. Kording

Abstract: Object-based attention is a key component of the visual system, relevant for perception, learning, and memory. Neurons tuned to features of attended objects tend to be more active than those associated with non-attended objects. There is a rich set of models of this phenomenon in computational neuroscience. However, there is currently a divide between models that successfully match physiological d… ▽ More Object-based attention is a key component of the visual system, relevant for perception, learning, and memory. Neurons tuned to features of attended objects tend to be more active than those associated with non-attended objects. There is a rich set of models of this phenomenon in computational neuroscience. However, there is currently a divide between models that successfully match physiological data but can only deal with extremely simple problems and models of attention used in computer vision. For example, attention in the brain is known to depend on top-down processing, whereas self-attention in deep learning does not. Here, we propose an artificial neural network model of object-based attention that captures the way in which attention is both top-down and recurrent. Our attention model works well both on simple test stimuli, such as those using images of handwritten digits, and on more complex stimuli, such as natural images drawn from the COCO dataset. We find that our model replicates a range of findings from neuroscience, including attention-invariant tuning, inhibition of return, and attention-mediated scaling of activity. Understanding object based attention is both computationally interesting and a key problem for computational neuroscience. △ Less

Submitted 8 June, 2021; originally announced June 2021.

Using single-cell entropy to describe the dynamics of reprogramming and differentiation of induced pluripotent stem cells

Authors: Yusong Ye, Zhuoqin Yang, Jinzhi Lei

Abstract: Induced pluripotent stem cells (iPSCs) provide a great model to study the process of reprogramming and differentiation of stem cells. Single-cell RNA sequencing (scRNA-seq) enables us to investigate the reprogramming process at single-cell level. Here, we introduce single-cell entropy (scEntropy) as a macroscopic variable to quantify the cellular transcriptome from scRNA-seq data during reprogramm… ▽ More Induced pluripotent stem cells (iPSCs) provide a great model to study the process of reprogramming and differentiation of stem cells. Single-cell RNA sequencing (scRNA-seq) enables us to investigate the reprogramming process at single-cell level. Here, we introduce single-cell entropy (scEntropy) as a macroscopic variable to quantify the cellular transcriptome from scRNA-seq data during reprogramming and differentiation of iPSCs. scEntropy measures the relative order parameter of genomic transcriptions at single cell level during the cell fate change process, which shows increasing during differentiation, and decreasing upon reprogramming. Moreover, based on the scEntropy dynamics, we construct a phenomenological stochastic differential equation model and the corresponding Fokker-Plank equation for cell state transitions during iPSC differentiation, which provide insights to infer cell fates changes and stem cell differentiation. This study is the first to introduce the novel concept of scEntropy to the biological process of iPSC, and suggests that the scEntropy can provide a suitable quantify to describe cell fate transition in differentiation and reprogramming of stem cells. △ Less

Submitted 16 April, 2020; originally announced April 2020.

Comments: 12 pages, 5 figures

DNA methylation heterogeneity induced by collaborations between enhancers

Authors: Yusong Ye, Zhuoqin Yang, Jinzhi Lei

Abstract: 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… ▽ More 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. △ Less

Submitted 15 February, 2020; originally announced February 2020.

Comments: 25 pages, 4 figures

Journal ref: Journal of Computational Biology, 2020

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

Authors: Jingxin Liu, You Song, Jinzhi Lei

Abstract: 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 t… ▽ More 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. △ Less

Submitted 15 February, 2020; originally announced February 2020.

Comments: 7 pages, 5 figures

Journal ref: Biophysical Reviews and Letters, 2020

Accurate Protein Structure Prediction by Embeddings and Deep Learning Representations

Authors: Iddo Drori, Darshan Thaker, Arjun Srivatsa, Daniel Jeong, Yueqi Wang, Linyong Nan, Fan Wu, Dimitri Leggas, Jinhao Lei, Weiyi Lu, Weilong Fu, Yuan Gao, Sashank Karri, Anand Kannan, Antonio Moretti, Mohammed AlQuraishi, Chen Keasar, Itsik Pe'er

Abstract: Proteins are the major building blocks of life, and actuators of almost all chemical and biophysical events in living organisms. Their native structures in turn enable their biological functions which have a fundamental role in drug design. This motivates predicting the structure of a protein from its sequence of amino acids, a fundamental problem in computational biology. In this work, we demonst… ▽ More Proteins are the major building blocks of life, and actuators of almost all chemical and biophysical events in living organisms. Their native structures in turn enable their biological functions which have a fundamental role in drug design. This motivates predicting the structure of a protein from its sequence of amino acids, a fundamental problem in computational biology. In this work, we demonstrate state-of-the-art protein structure prediction (PSP) results using embeddings and deep learning models for prediction of backbone atom distance matrices and torsion angles. We recover 3D coordinates of backbone atoms and reconstruct full atom protein by optimization. We create a new gold standard dataset of proteins which is comprehensive and easy to use. Our dataset consists of amino acid sequences, Q8 secondary structures, position specific scoring matrices, multiple sequence alignment co-evolutionary features, backbone atom distance matrices, torsion angles, and 3D coordinates. We evaluate the quality of our structure prediction by RMSD on the latest Critical Assessment of Techniques for Protein Structure Prediction (CASP) test data and demonstrate competitive results with the winning teams and AlphaFold in CASP13 and supersede the results of the winning teams in CASP12. We make our data, models, and code publicly available. △ Less

Submitted 8 November, 2019; originally announced November 2019.

Journal ref: Machine Learning in Computational Biology, 2019

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

Authors: Jinzhi Lei

Abstract: 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 essen… ▽ More 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. △ Less

Submitted 19 August, 2019; originally announced August 2019.

Comments: 19 pages, 3 figures

A general mathematical framework for understanding the behavior of heterogeneous stem cell regeneration

Authors: Jinzhi Lei

Abstract: Stem cell heterogeneity is essential for the homeostasis in tissue development. This paper established a general formulation for understanding the dynamics of stem cell regeneration with cell heterogeneity and random transitions of epigenetic states. The model generalizes the classical G0 cell cycle model, and incorporates the epigenetic states of stem cells that are represented by a continuous mu… ▽ More Stem cell heterogeneity is essential for the homeostasis in tissue development. This paper established a general formulation for understanding the dynamics of stem cell regeneration with cell heterogeneity and random transitions of epigenetic states. The model generalizes the classical G0 cell cycle model, and incorporates the epigenetic states of stem cells that are represented by a continuous multidimensional variable and the kinetic rates of cell behaviors, including proliferation, differentiation, and apoptosis, that are dependent on their epigenetic states. Moreover, the random transition of epigenetic states is represented by an inheritance probability that can be described as a conditional beta distribution. This model can be extended to investigate gene mutation-induced tumor development. The proposed formula is a generalized formula that helps us to understand various dynamic processes of stem cell regeneration, including tissue development, degeneration, and abnormal growth. △ Less

Submitted 11 July, 2024; v1 submitted 27 March, 2019; originally announced March 2019.

Comments: 36 pages, 7 figures

Journal ref: Journal of Theoretical Biology, 2020

aa-tRNA competition is crucial for the effective translation efficiency

Authors: Wenjun Xia, Jinzhi Lei

Abstract: 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 effic… ▽ More 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. △ Less

Submitted 30 March, 2015; originally announced April 2015.

Comments: 8 pages, 10 figures

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

Authors: Changjing Zhuge, Xiaojuan Sun, Yingyu Chen, Jinzhi Lei

Abstract: 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… ▽ More 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. △ Less

Submitted 28 March, 2015; originally announced March 2015.

Comments: 12 pages, 7 figures

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

Authors: Yuanhong Bi, Zhuoqin Yang, Changjing Zhuge, Jinzhi Lei

Abstract: 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 act… ▽ More 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. △ Less

Submitted 27 March, 2015; originally announced March 2015.

Comments: 11 pages, 8 figures

Stochastic Modeling in Systems Biology

Authors: Jinzhi Lei

Abstract: 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 se… ▽ More 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. △ Less

Submitted 22 April, 2011; originally announced April 2011.

Comments: 25 pages, 4 figures

Elastic energy of proteins and the stages of protein folding

Authors: Jinzhi Lei, Kerson Huang

Abstract: 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 ela… ▽ More 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. △ Less

Submitted 26 February, 2010; originally announced February 2010.

Journal ref: Europhys Lett, 88 (2009), 68004

Protein Folding: A Perspective From Statistical Physics

Authors: Jinzhi Lei, Kerson Huang

Abstract: 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.… ▽ More 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. △ Less

Submitted 26 February, 2010; originally announced February 2010.

Journal ref: Protein Folding, Eric C. Walters ed. NOVA Publishers, 2010, pp. 579-604

A Unified Model of $α$-Helix/$β$-Sheet/Random-Coil Transition in Proteins

Authors: Liu Hong, JinZhi Lei

Abstract: The theory of transition between $α$-helix, $β$-sheet and random coil conformation of a protein is discussed through a simple model, that includes both short and long-range interactions. Besides the bonding parameter and helical initiation factor in Zimm-Bragg model, three new parameters are introduced to describe beta structure: the local constraint factor for a single residue to be contained i… ▽ More The theory of transition between $α$-helix, $β$-sheet and random coil conformation of a protein is discussed through a simple model, that includes both short and long-range interactions. Besides the bonding parameter and helical initiation factor in Zimm-Bragg model, three new parameters are introduced to describe beta structure: the local constraint factor for a single residue to be contained in a $β$-strand, the long-range bonding parameter that accounts for the interaction between a pair of bonded $β$-strands, and a correction factor for the initiation of a $β$-sheet. Either increasing local constraint factor or long-range bonding parameter can cause a transition from $α$-helix or random coil conformation to $β$-sheet conformation. The sharpness of transition depends on the competition between short and long-range interactions. Other effective factors, such as the chain length and temperature, are also discussed. In this model, the entropy due to different ways to group $β$-strands into different $β$-sheets gives rise to significant contribution to partition function, and makes major differences between beta structure and helical structure. △ Less

Submitted 6 November, 2007; originally announced November 2007.

Comments: 8 pages, 7 figures

Sequential recruitment and combinatorial assembling of multiprotein complexes in transcriptional activation

Authors: Vincent Lemaire, Chiu Fan Lee, Jinzhi Lei, Raphael Metivier, Leon Glass

Abstract: In human cells, estrogenic signals induce cyclical association and dissociation of specific proteins with the DNA in order to activate transcription of estrogen-responsive genes. These oscillations can be modeled by assuming a large number of sequential reactions represented by linear kinetics with random kinetic rates. Application of the model to experimental data predicts robust binding sequen… ▽ More In human cells, estrogenic signals induce cyclical association and dissociation of specific proteins with the DNA in order to activate transcription of estrogen-responsive genes. These oscillations can be modeled by assuming a large number of sequential reactions represented by linear kinetics with random kinetic rates. Application of the model to experimental data predicts robust binding sequences in which proteins associate with the DNA at several different phases of the oscillation. Our methods circumvent the need to derive detailed kinetic graphs, and are applicable to other oscillatory biological processes involving a large number of sequential steps. △ Less

Submitted 30 May, 2006; originally announced May 2006.

Journal ref: Phys. Rev. Lett. 96, 198102 (2006)

Spectral Representation Theory for Dielectric Behavior of Nonspherical Cell Suspensions

Authors: J. P. Huang, K. W. Yu, Jun Lei, Hong Sun

Abstract: Recent experiments revealed that the dielectric dispersion spectrum of fission yeast cells in a suspension was mainly composed of two sub-dispersions. The low-frequency sub-dispersion depended on the cell length, while the high-frequency one was independent of it. The cell shape effect was simulated by an ellipsoidal cell model but the comparison between theory and experiment was far from being… ▽ More Recent experiments revealed that the dielectric dispersion spectrum of fission yeast cells in a suspension was mainly composed of two sub-dispersions. The low-frequency sub-dispersion depended on the cell length, while the high-frequency one was independent of it. The cell shape effect was simulated by an ellipsoidal cell model but the comparison between theory and experiment was far from being satisfactory. Prompted by the discrepancy, we proposed the use of spectral representation to analyze more realistic cell models. We adopted a shell-spheroidal model to analyze the effects of the cell membrane. It is found that the dielectric property of the cell membrane has only a minor effect on the dispersion magnitude ratio and the characteristic frequency ratio. We further included the effect of rotation of dipole induced by an external electric field, and solved the dipole-rotation spheroidal model in the spectral representation. Good agreement between theory and experiment has been obtained. △ Less

Submitted 23 April, 2001; originally announced April 2001.

Comments: 19 pages, 5 eps figures

Journal ref: Comm. Theor. Phys. 38, 113 (2002)

Dielectric Behavior of Nonspherical Cell Suspensions

Authors: Jun Lei, Jones T. K. Wan, K. W. Yu, Hong Sun

Abstract: Recent experiments revealed that the dielectric dispersion spectrum of fission yeast cells in a suspension was mainly composed of two sub-dispersions. The low-frequency sub-dispersion depended on the cell length, whereas the high-frequency one was independent of it. The cell shape effect was qualitatively simulated by an ellipsoidal cell model. However, the comparison between theory and experime… ▽ More Recent experiments revealed that the dielectric dispersion spectrum of fission yeast cells in a suspension was mainly composed of two sub-dispersions. The low-frequency sub-dispersion depended on the cell length, whereas the high-frequency one was independent of it. The cell shape effect was qualitatively simulated by an ellipsoidal cell model. However, the comparison between theory and experiment was far from being satisfactory. In an attempt to close up the gap between theory and experiment, we considered the more realistic cells of spherocylinders, i.e., circular cylinders with two hemispherical caps at both ends. We have formulated a Green function formalism for calculating the spectral representation of cells of finite length. The Green function can be reduced because of the azimuthal symmetry of the cell. This simplification enables us to calculate the dispersion spectrum and hence access the effect of cell structure on the dielectric behavior of cell suspensions. △ Less

Submitted 23 March, 2001; originally announced March 2001.

Comments: Preliminary results have been reported in the 2001 March Meeting of the American Physical Society. Accepted for publications in J. Phys.: Condens. Matter

Journal ref: J. Phys.: Condens. Matter 13, 3583 (2001).