Researcher profile

Jaroslav Albert

Jaroslav Albert contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Molecular Networks Biomolecules Cell Behavior cond-mat.soft Quantitative Methods Subcellular Processes

Trust snapshot

Quick read

Trust 19 - Baseline

5works

0followers

6topics

2close collaborators

Actions

Decide how to stay connected

Follow researcher0

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

preprint2023arXiv

A tau-leaping method for computing joint probability distributions of the first-passage time and position of a Brownian particle

First passage time (FPT) is the time a particle, subject to some stochastic process, hits or crosses a closed surface for the very first time. $τ$-leaping methods are a class of stochastic algorithms in which, instead of simulating every single reaction, many reactions are ``leaped" over in order to shorten the computing time. In this paper we developed a $τ$-leaping method for computing the FPT and position in arbitrary volumes for a Brownian particle governed by the Langevin equation. The $τ$-leaping method proposed here works as follows. A sphere is inscribed within the volume of interest (VOI) centered at the initial particle's location. On this sphere, the FPT is sampled, as well as the position, which becomes the new initial position. Then, another sphere, centered at this new location, is inscribed. This process continues until the sphere becomes smaller than some minimal radius $R_{\text{min}}$. When this occurs, the $τ$-leaping switches to the conventional Monte Carlo, which runs until the particle either crosses the surface of the VOI or finds its way to a position where a sphere of radius $>R_{\text{min}}$ can be inscribed. The switching between $τ$-leaping and MC continues until the particle crosses the surface of the VOI. The size of this radius depends on the system parameters and on one's notion of accuracy: the larger this radius the more accurate the $τ$-leaping method, but also less efficient. This trade off between accuracy and efficiency is discussed. For two VOI, the $τ$-leaping method is shown to be accurate and more efficient than MC by at least a factor of 10 and up to a factor of about 110. However, while MC becomes exponentially slower with increasing VOI, the efficiency of the $τ$-leaping method remains relatively unchanged. Thus, the $τ$-leaping method can potentially be many orders of magnitude more efficient than MC.

preprint2022arXiv

A detailed model of gene promoter dynamics reveals the entry into productive elongation to be a highly punctual process

Gene transcription is a stochastic process that involves thousands of reactions. The first set of these reactions, which happen near a gene promoter, are considered to be the most important in the context of stochastic noise. The most common models of transcription are primarily concerned with the effect of activators/repressors on the overall transcription rate and approximate the basal transcription processes as a one step event. According to such effective models, the Fano factor of mRNA copy distributions is always greater than (super-Poissonian) or equal to 1 (Poissonian), and the only way to go below this limit (sub-Poissonian) is via a negative feedback. It is partly due to this limit that the first stage of transcription is held responsible for most of the stochastic noise in mRNA copy numbers. However, by considering all major reactions that build and drive the basal transcription machinery, from the first protein that binds a promoter to the entrance of the transcription complex (TC) into productive elongation, it is shown that the first two stages of transcription, namely the pre-initiation complex (PIC) formation and the promoter proximal pausing (PPP), is a highly punctual process. In other words, the time between the first and the last step of this process is narrowly distributed, which gives rise to sub-Poissonian distributions for the number of TCs that have entered productive elongation. In fact, having simulated the PIC formation and the PPP via the Gillespie algorithm using 2000 distinct parameter sets and 4 different reaction network topologies, it is shown that only 4.4% give rise to a Fano factor that is > 1 with the upper bound of 1.7, while for 31% of cases the Fano factor is below 0.5, with 0.19 as the lower bound. These results cast doubt on the notion that most of the stochastic noise observed in mRNA distributions always originates at the promoter.

preprint2020arXiv

Dimensionality reduction via path integration for computing mRNA distributions

Inherent stochasticity in gene expression leads to distributions of mRNA copy numbers in a population of identical cells. These distributions are determined primarily by the multitude of states of a gene promoter, each driving transcription at a different rate. In an era where single-cell mRNA copy number data are more and more available, there is an increasing need for fast computations of mRNA distributions. In this paper, we present a method for computing separate distributions for each species of mRNA molecules, i. e. mRNAs that have been either partially or fully processed post-transcription. The method involves the integration over all possible realizations of promoter states, which we cast into a set of linear ordinary differential equations of dimension $M\times n_j$, where $M$ is the number of available promoter states and $n_j$ is the mRNA copy number of species $j$ up to which one wishes to compute the probability distribution. This approach is superior to solving the Master equation (ME) directly in two ways: a) the number of coupled differential equations in the ME approach is $M\timesΛ_1\timesΛ_2\times ...\timesΛ_L$, where $Λ_j$ is the cutoff for the probability of the $j^{\text{th}}$ species of mRNA; and b) the ME must be solved up to the cutoffs $Λ_j$, which are {\it ad hoc} and must be selected {\it a priori}. In our approach, the equation for the probability to observe $n$ mRNAs of any species depends only on the the probability of observing $n-1$ mRNAs of that species, thus yielding a correct probability distribution up to an arbitrary $n$. To demonstrate the validity of our derivations, we compare our results with Gillespie simulations for ten randomly selected system parameters.

preprint2014arXiv

Stochastic noise reduction upon complexification: positively correlated birth-death type systems

Cell systems consist of a huge number of various molecules that display specific patterns of interactions, which have a determining influence on the cell's functioning. In general, such complexity is seen to increase with the complexity of the organism, with a concomitant increase of the accuracy and specificity of the cellular processes. The question thus arises how the complexification of systems - modeled here by simple interacting birth-death type processes - can lead to a reduction of the noise - described by the variance of the number of molecules. To gain understanding of this issue, we investigated the difference between a single system containing molecules that are produced and degraded, and the same system - with the same average number of molecules - connected to a buffer. We modeled these systems using Ito stochastic differential equations in discrete time, as they allow straightforward analytical developments. In general, when the molecules in the system and the buffer are positively correlated, the variance on the number of molecules in the system is found to decrease compared to the equivalent system without a buffer. Only buffers that are too noisy by themselves tend to increase the noise in the main system. We tested this result on two model cases, in which the system and the buffer contain proteins in their active and inactive state, or protein monomers and homodimers. We found that in the second test case, where the interconversion terms are non-linear in the number of molecules, the noise reduction is much more pronounced; it reaches up to 20% reduction of the Fano factor with the parameter values tested in numerical simulations on an unperturbed birth-death model. We extended our analysis to two arbitrary interconnected systems.

preprint2011arXiv

Dynamic modeling of gene expression in prokaryotes: application to glucose-lactose diauxie in Escherichia coli

Coexpression of genes or, more generally, similarity in the expression profiles poses an unsurmountable obstacle to inferring the gene regulatory network (GRN) based solely on data from DNA microarray time series. Clustering of genes with similar expression profiles allows for a course-grained view of the GRN and a probabilistic determination of the connectivity among the clusters. We present a model for the temporal evolution of a gene cluster network which takes into account interactions of gene products with genes and, through a non-constant degradation rate, with other gene products. The number of model parameters is reduced by using polynomial functions to interpolate temporal data points. In this manner, the task of parameter estimation is reduced to a system of linear algebraic equations, thus making the computation time shorter by orders of magnitude. To eliminate irrelevant networks, we test each GRN for stability with respect to parameter variations, and impose restrictions on its behavior near the steady state. We apply our model and methods to DNA microarray time series' data collected on Escherichia coli during glucose-lactose diauxie and infer the most probable cluster network for different phases of the experiment.