Nascent RNA kinetics with complex promoter architecture: Analytic results and parameter inference.

Transcription is a stochastic process that involves several downstream operations which make it difficult to model and infer transcription kinetics from mature RNA numbers in individual cell. However, recent advances in single-cell technologies have enabled a more precise measurement of the fluctuations of nascent RNA that closely reflect transcription kinetics. In this paper we introduce a general stochastic model to mimic nascent RNA kinetics with complex promoter architecture.

Exact burst-size distributions for gene-expression models with complex promoter structure.

In prokaryotic and eukaryotic cells, most genes are transcribed in a bursty fashion on one hand and complex gene regulations may lead to complex promoter structure on the other hand. This raises an unsolved issue: how does promoter structure shape transcriptional bursting kinetics characterized by burst size and frequency? Here we analyze stochastic models of gene transcription, which consider complex regulatory mechanisms. Notably, we develop an efficient method to derive exact burst-size distributions.

Reconstructing the evolution history of networked complex systems.

The evolution processes of complex systems carry key information in the systems' functional properties. Applying machine learning algorithms, we demonstrate that the historical formation process of various networked complex systems can be extracted, including protein-protein interaction, ecology, and social network systems. The recovered evolution process has demonstrations of immense scientific values, such as interpreting the evolution of protein-protein interaction network, facilitating structure prediction, and particularly revealing the key co-evolution features of network structures such as preferential attachment, community structure, local clustering, degree-degree correlation that could not be explained collectively by previous theories.

Exact results for gene-expression models with general waiting-time distributions.

Complex molecular details of transcriptional regulation can be coarse-grained by assuming that reaction waiting times for promoter-state transitions, the mRNA synthesis, and the mRNA degradation follow general distributions. However, how such a generalized two-state model is analytically solved is a long-standing issue. Here we first present analytical formulas of burst-size distributions for this model.

Balanced implicit Patankar-Euler methods for positive solutions of stochastic differential equations of biological regulatory systems.

Stochastic differential equations (SDEs) are a powerful tool to model fluctuations and uncertainty in complex systems. Although numerical methods have been designed to simulate SDEs effectively, it is still problematic when numerical solutions may be negative, but application problems require positive simulations. To address this issue, we propose balanced implicit Patankar-Euler methods to ensure positive simulations of SDEs.

Power-law behavior of transcriptional bursting regulated by enhancer-promoter communication.

Revealing how transcriptional bursting kinetics are genomically encoded is challenging because genome structures are stochastic at the organization level and are suggestively linked to gene transcription. To address this challenge, we develop a generic theoretical framework that integrates chromatin dynamics, enhancer-promoter (E-P) communication, and gene-state switching to study transcriptional bursting. The theory predicts that power law can be a general rule to quantitatively describe bursting modulations by E-P spatial communication.

4D nucleome equation predicts gene expression controlled by long-range enhancer-promoter interaction.

Recent experimental evidence strongly supports that three-dimensional (3D) long-range enhancer-promoter (E-P) interactions have important influences on gene-expression dynamics, but it is unclear how the interaction information is translated into gene expression over time (4D). To address this question, we developed a general theoretical framework (named as a 4D nucleome equation), which integrates E-P interactions on chromatin and biochemical reactions of gene transcription. With this equation, we first present the distribution of mRNA counts as a function of the E-P genomic distance and then reveal a power-law scaling of the expression level in this distance.

Phase separation reduces cell-to-cell variability of transcriptional bursting.

Gene expression is a stochastic and noisy process often occurring in "bursts". Experiments have shown that the compartmentalization of proteins by liquid-liquid phase separation is conducive to reducing the noise of gene expression. Therefore, an important goal is to explore the role of bursts in phase separation noise reduction processes.

Machine learning study of the extended drug-target interaction network informed by pain related voltage-gated sodium channels.

Pain is a significant global health issue, and the current treatment options for pain management have limitations in terms of effectiveness, side effects, and potential for addiction. There is a pressing need for improved pain treatments and the development of new drugs. Voltage-gated sodium channels, particularly Nav1.

Stochastic modeling of the mRNA life process: A generalized master equation.

The mRNA life cycle is a complex biochemical process, involving transcription initiation, elongation, termination, splicing, and degradation. Each of these molecular events is multistep and can create a memory. The effect of this molecular memory on gene expression is not clear, although there are many related yet scattered experimental reports.

Machine Learning Study of the Extended Drug-target Interaction Network informed by Pain Related Voltage-Gated Sodium Channels.

Pain is a significant global health issue, and the current treatment options for pain management have limitations in terms of effectiveness, side effects, and potential for addiction. There is a pressing need for improved pain treatments and the development of new drugs. Voltage-gated sodium channels, particularly Nav1.

Composite Patankar-Euler methods for positive simulations of stochastic differential equation models for biological regulatory systems.

Stochastic differential equations (SDE) are a powerful tool to model biological regulatory processes with intrinsic and extrinsic noise. However, numerical simulations of SDE models may be problematic if the values of noise terms are negative and large, which is not realistic for biological systems since the molecular copy numbers or protein concentrations should be non-negative. To address this issue, we propose the composite Patankar-Euler methods to obtain positive simulations of SDE models.

Inferring transcriptional bursting kinetics from single-cell snapshot data using a generalized telegraph model.

Gene expression has inherent stochasticity resulting from transcription's burst manners. Single-cell snapshot data can be exploited to rigorously infer transcriptional burst kinetics, using mathematical models as blueprints. The classical telegraph model (CTM) has been widely used to explain transcriptional bursting with Markovian assumptions.

TIDAL: Topology-Inferred Drug Addiction Learning.

Drug addiction is a global public health crisis, and the design of antiaddiction drugs remains a major challenge due to intricate mechanisms. Since experimental drug screening and optimization are too time-consuming and expensive, there is urgent need to develop innovative artificial intelligence (AI) methods for addressing the challenge. We tackle this challenge by topology-inferred drug addiction learning (TIDAL) built from integrating multiscale topological Laplacians, deep bidirectional transformer, and ensemble-assisted neural networks (EANNs).

Genome-wide inference reveals that feedback regulations constrain promoter-dependent transcriptional burst kinetics.

Gene expression in mammalian cells is highly variable and episodic, resulting in a series of discontinuous bursts of mRNAs. A challenge is to understand how static promoter architecture and dynamic feedback regulations dictate bursting on a genome-wide scale. Although single-cell RNA sequencing (scRNA-seq) provides an opportunity to address this challenge, effective analytical methods are scarce.

Silent transcription intervals and translational bursting lead to diverse phenotypic switching.

Gene-expression bimodality, as a potential mechanism generating phenotypic cell diversity, can enhance the survival of cells in a fluctuating environment. Previous studies have shown that intrinsic or extrinsic regulations could induce bimodal gene expressions, but it is unclear whether this bimodality can occur without regulation. Here we develop an interpretable and tractable model, namely a generalized telegraph model (GTM), which considers silent transcription intervals and translational bursting, each being characterized by a general distribution.

Generalized fluctuation-dissipation theorem for non-Markovian reaction networks.

Intracellular biochemical networks often display large fluctuations in the molecule numbers or the concentrations of reactive species, making molecular approaches necessary for system descriptions. For Markovian reaction networks, the fluctuation-dissipation theorem (FDT) has been well established and extensively used in fast evaluation of fluctuations in reactive species. For non-Markovian reaction networks, however, the similar FDT has not been established so far.

Analysis of a non-Markov transcription model with nuclear RNA export and RNA nuclear retention.

Transcription involves gene activation, nuclear RNA export (NRE) and RNA nuclear retention (RNR). All these processes are multistep and biochemical. A multistep reaction process can create memories between reaction events, leading to non-Markovian kinetics.

Kinetic characteristics of transcriptional bursting in a complex gene model with cyclic promoter structure.

While transcription often occurs in a bursty manner, various possible regulations can lead to complex promoter patterns such as promoter cycles, giving rise to an important question: How do promoter kinetics shape transcriptional bursting kinetics? Here we introduce and analyze a general model of the promoter cycle consisting of multi-OFF states and multi-ON states, focusing on the effects of multi-ON mechanisms on transcriptional bursting kinetics. The derived analytical results indicate that burst size follows a mixed geometric distribution rather than a single geometric distribution assumed in previous studies, and ON and OFF times obey their own mixed exponential distributions. In addition, we find that the multi-ON mechanism can lead to bimodal burst-size distribution, antagonistic timing of ON and OFF, and diverse burst frequencies, each further contributing to cell-to-cell variability in the mRNA expression level.

A robust method for designing multistable systems by embedding bistable subsystems.

Although multistability is an important dynamic property of a wide range of complex systems, it is still a challenge to develop mathematical models for realising high order multistability using realistic regulatory mechanisms. To address this issue, we propose a robust method to develop multistable mathematical models by embedding bistable models together. Using the GATA1-GATA2-PU.

Exact distributions for stochastic gene expression models with arbitrary promoter architecture and translational bursting.

Gene expression in individual cells is inherently variable and sporadic, leading to cell-to-cell variability in mRNA and protein levels. Recent single-cell and single-molecule experiments indicate that promoter architecture and translational bursting play significant roles in controlling gene expression noise and generating the phenotypic diversity that life exhibits. To quantitatively understand the impact of these factors, it is essential to construct an accurate mathematical description of stochastic gene expression and find the exact analytical results, which is a formidable task.

Exploring the role of eRNA in regulating gene expression.

eRNAs as the products of enhancers can regulate gene expression via various possible ways, but which regulation way is more reasonable is debatable in biology, and in particular, how eRNAs impact gene expression remains unclear. Here we introduce a mechanistic model of gene expression to address these issues. This model considers three possible regulation ways of eRNA: Type-I by which eRNA regulates transcriptional activity by facilitating the formation of enhancer-promoter (E-P) loop, Type-II by which eRNA directly promotes the mRNA production rate, and mixed regulation (i.

Exact results for queuing models of stochastic transcription with memory and crosstalk.

Gene transcription is a complex multistep biochemical process, which can create memory between individual reaction events. On the other hand, many inducible genes, when activated by external cues, are often coregulated by several competitive pathways with crosstalk. This raises an unexplored question: how do molecular memory and crosstalk together affect gene expressions? To address this question, we introduce a queuing model of stochastic transcription, where two crossing signaling pathways are used to direct gene activation in response to external signals and memory functions to model multistep reaction processes involved in transcription.

Correlation between external regulators governs the mean-noise relationship in stochastic gene expression.

Gene transcription in single cells is inherently a probabilistic process. The relationship between variance ($ \sigma^{2} $) and mean expression ($ \mu $) is of paramount importance for investigations into the evolutionary origins and consequences of noise in gene expression. It is often formulated as $ \log \left({{{\sigma}^{2}}}/{{{\mu}^{2}}}\; \right) = \beta\log\mu+\log\alpha $, where $ \beta $ is a key parameter since its sign determines the qualitative dependence of noise on mean.

Exploring dissipative sources of non-Markovian biochemical reaction systems.

Many biological processes including important intracellular processes are governed by biochemical reaction networks. Usually, these reaction systems operate far from thermodynamic equilibrium, implying free-energy dissipation. On the other hand, single reaction events happen often in a memory manner, leading to non-Markovian kinetics.