Researcher profile

Abhyudai Singh

Abhyudai Singh contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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.

Published work

15 published item(s)

preprint2025arXiv

Stochastic multi-step cell size homeostasis model for cycling human cells

Measurements of cell size dynamics have established the adder principle as a robust mechanism of cell size homeostasis. In this framework, cells add a nearly constant amount of size during each cell cycle, independent of their size at birth. Theoretical studies have shown that the adder principle can be achieved when cell-cycle progression is coupled to cell size. Here, we extend this framework by considering a general growth law modeled as a Hill-type function of cell size. This assumption introduces growth saturation to the model, such that very large cells grow approximately linearly rather than exponentially. Additionally, to capture the sequential nature of division, we implement a stochastic multi-step adder model in which cells progress through internal regulatory stages before dividing. From this model, we derive exact analytical expressions for the moments of cell size distributions. Our results show that stronger growth saturation increases the mean cell size in steady state, while slightly reducing fluctuations compared to exponential growth. Importantly, despite these changes, the adder property is preserved. This emphasizes that the reduction in size variability is a consequence of~the growth law rather than simple scaling with mean size. Finally, we analyze stochastic clonal proliferation and find that growth saturation influences both single-cell size statistics and variability across populations. Our results provide a generalized framework for connecting multi-step adder mechanisms with proliferation dynamics, extending size control theory beyond exponential growth.

preprint2019arXiv

Global redistribution and local migration in semi-discrete host-parasitoid population dynamic models

Host-parasitoid population dynamics is often probed using a semi-discrete/hybrid modeling framework. Here, the update functions in the discrete-time model connecting year-to-year changes in the population densities are obtained by solving ordinary differential equations that mechanistically describe interactions when hosts become vulnerable to parasitoid attacks. We use this semi-discrete formalism to study two key spatial effects: local movement (migration) of parasitoids between patches during the vulnerable period; and yearly redistribution of populations across patches outside the vulnerable period. Our results show that in the absence of any redistribution, constant density-independent migration and parasitoid attack rates are unable to stabilize an otherwise unstable host-parasitoid population dynamics. Interestingly, inclusion of host redistribution (but not parasitoid redistribution) before the start of the vulnerable period can lead to stable coexistence of both species. Next, we consider a Type-III functional response (parasitoid attack rate increases with host density), where the absence of any spatial effects leads to a neutrally stable host-parasitoid equilibrium. As before, density-independent parasitoid migration by itself is again insufficient to stabilize the population dynamics and host redistribution provides a stabilizing influence. Finally, we show that a Type-III functional response combined with density-dependent parasitoid migration leads to stable coexistence, even in the absence of population redistributions. In summary, we have systematically characterized parameter regimes leading to stable/unstable population dynamics with different forms of spatial heterogeneity coupled to the parasitoid's functional response using mechanistically formulated semi-discrete models.

preprint2016arXiv

Cell-cycle coupled expression minimizes random fluctuations in gene product levels

Expression of many genes varies as a cell transitions through different cell-cycle stages. How coupling between stochastic expression and cell cycle impacts cell-to-cell variability (noise) in the level of protein is not well understood. We analyze a model, where a stable protein is synthesized in random bursts, and the frequency with which bursts occur varies within the cell cycle. Formulas quantifying the extent of fluctuations in the protein copy number are derived and decomposed into components arising from the cell cycle and stochastic processes. The latter stochastic component represents contributions from bursty expression and errors incurred during partitioning of molecules between daughter cells. These formulas reveal an interesting trade-off: cell-cycle dependencies that amplify the noise contribution from bursty expression also attenuate the contribution from partitioning errors. We investigate existence of optimum strategies for coupling expression to the cell cycle that minimize the stochastic component. Intriguingly, results show that a zero production rate throughout the cell cycle, with expression only occurring just before cell division minimizes noise from bursty expression for a fixed mean protein level. In contrast, the optimal strategy in the case of partitioning errors is to make the protein just after cell division. We provide examples of regulatory proteins that are expressed only towards the end of cell cycle, and argue that such strategies enhance robustness of cell-cycle decisions to the intrinsic stochasticity of gene expression.

preprint2016arXiv

Conditions for cell size homeostasis: A stochastic hybrid systems approach

A ubiquitous feature of living cells is their growth over time followed by division into daughter cells. How isogenic cell populations maintain size homeostasis, i.e., a narrow distribution of cell size, is an intriguing fundamental problem. We model cell size using a stochastic hybrid system, where a cell grows exponentially in size (volume) over time and probabilistic division events are triggered at discrete time intervals. Moreover, whenever division events occur, size is randomly partitioned among daughter cells. We first consider a scenario, where a timer (i.e., cell-cycle clock) that measures the time since the last division event regulates both the cellular growth and division rates. Analysis reveals that such a timer-controlled system cannot achieve size homeostasis, in the sense that, the cell-to-cell size variation grows unboundedly with time. To explore biologically meaningful mechanisms for controlling size we consider two classes of regulation: a size-dependent growth rate and a size-dependent division rate. Our results show that these strategies can provide bounded intercellular variation in cell size, and exact mathematical conditions on the form of regulation needed for size homeostasis are derived. Different known forms of size control strategies, such as, the adder and the sizer are shown to be consistent with these results. Interestingly, for timer-based division mechanisms, the mean cell size depends on the noise in the cell-cycle duration but independent of errors incurred in partitioning of volume among daughter cells. In contrast, the mean cell size decreases with increasing partitioning errors for size-based division mechanisms. Finally, we discuss how organisms ranging from bacteria to mammalian cells have adopted different control approaches for maintaining size homeostasis.

preprint2016arXiv

Effect of gene-expression bursts on stochastic timing of cellular events

Gene expression is inherently a noisy process which manifests as cell-to-cell variability in time evolution of proteins. Consequently, events that trigger at critical threshold levels of regulatory proteins exhibit stochasticity in their timing. An important contributor to the noise in gene expression is translation bursts which correspond to randomness in number of proteins produced in a single mRNA lifetime. Modeling timing of an event as a first-passage time (FPT) problem, we explore the effect of burst size distribution on event timing. Towards this end, the probability density function of FPT is computed for a gene expression model with burst size drawn from a generic non-negative distribution. Analytical formulas for FPT moments are provided in terms of known vectors and inverse of a matrix. The effect of burst size distribution is investigated by looking at how the feedback regulation strategy that minimizes noise in timing around a given time deviates from the case when burst is deterministic. Interestingly, results show that the feedback strategy for deterministic burst case is quite robust to change in burst size distribution, and deviations from it are confined to about 20% of the optimal value. These findings facilitate an improved understanding of noise regulation in event timing.

preprint2016arXiv

Gene expression noise is affected differentially by feedback in burst frequency and burst size

Inside individual cells, expression of genes is stochastic across organisms ranging from bacterial to human cells. A ubiquitous feature of stochastic expression is burst-like synthesis of gene products, which drives considerable intercellular variability in protein levels across an isogenic cell population. One common mechanism by which cells control such stochasticity is negative feedback regulation, where a protein inhibits its own synthesis. For a single gene that is expressed in bursts, negative feedback can affect the burst frequency or the burst size. In order to compare these feedback types, we study a piecewise deterministic model for gene expression of a self-regulating gene. Mathematically tractable steady-state protein distributions are derived and used to compare the noise suppression abilities of the two feedbacks. Results show that in the low noise regime, both feedbacks are similar in term of their noise buffering abilities. Intriguingly, feedback in burst size outperforms the feedback in burst frequency in the high noise regime. Finally, we discuss various regulatory strategies by which cells implement feedback to control burst sizes of expressed proteins at the level of single cells.

preprint2016arXiv

Moment dynamics for a class of time-triggered stochastic hybrid systems

Stochastic Hybrid Systems (SHS) constitute an important class of mathematical models that integrate discrete stochastic events with continuous dynamics. The time evolution of statistical moments is generally not closed for SHS, in the sense that the time derivative of the lower-order moments depends on higher-order moments. Here, we identify an important class of SHS where moment dynamics is automatically closed, and hence moments can be computed exactly by solving a system of coupled differential equations. This class is referred to as linear time-triggered SHS (TTSHS), where the state evolves according to a linear dynamical system. Stochastic events occur at discrete times and the intervals between them are independent random variables that follow a general class of probability distributions. Moreover, whenever the event occurs, the state of the SHS changes randomly based on a probability distribution. Our approach relies on embedding a Markov chain based on phase-type processes to model timing of events, and showing that the resulting system has closed moment dynamics. Interestingly, we identify a subclass of linear TTSHS, where the first and second-order moments depend only on the mean time interval between events and invariant of their higher-order statistics. Finally, we discuss applicability of our results to different application areas such as network control systems and systems biology.

preprint2015arXiv

A mechanistic first--passage time framework for bacterial cell-division timing

How exponentially growing cells maintain size homeostasis is an important fundamental problem. Recent single-cell studies in prokaryotes have uncovered the adder principle, where cells on average, add a fixed size (volume) from birth to division. Interestingly, this added volume differs considerably among genetically-identical newborn cells with similar sizes suggesting a stochastic component in the timing of cell-division. To mechanistically explain the adder principle, we consider a time-keeper protein that begins to get stochastically expressed after cell birth at a rate proportional to the volume. Cell-division time is formulated as the first-passage time for protein copy numbers to hit a fixed threshold. Consistent with data, the model predicts that while the mean cell-division time decreases with increasing size of newborns, the noise in timing increases with size at birth. Intriguingly, our results show that the distribution of the volume added between successive cell-division events is independent of the newborn cell size. This was dramatically seen in experimental studies, where histograms of the added volume corresponding to different newborn sizes collapsed on top of each other. The model provides further insights consistent with experimental observations: the distributions of the added volume and the cell-division time when scaled by their respective means become invariant of the growth rate. Finally, we discuss various modifications to the proposed model that lead to deviations from the adder principle. In summary, our simple yet elegant model explains key experimental findings and suggests a mechanism for regulating both the mean and fluctuations in cell-division timing for size control.

preprint2015arXiv

Decomposing variability in protein levels from noisy expression, genome duplication and partitioning errors during cell-divisions

Inside individual cells, expression of genes is inherently stochastic and manifests as cell-to-cell variability or noise in protein copy numbers. Since proteins half-lives can be comparable to the cell-cycle length, randomness in cell-division times generates additional intercellular variability in protein levels. Moreover, as many mRNA/protein species are expressed at low-copy numbers, errors incurred in partitioning of molecules between the mother and daughter cells are significant. We derive analytical formulas for the total noise in protein levels for a general class of cell-division time and partitioning error distributions. Using a novel hybrid approach the total noise is decomposed into components arising from i) stochastic expression; ii) partitioning errors at the time of cell-division and iii) random cell-division events. These formulas reveal that random cell-division times not only generate additional extrinsic noise but also critically affect the mean protein copy numbers and intrinsic noise components. Counter intuitively, in some parameter regimes noise in protein levels can decrease as cell-division times become more stochastic. Computations are extended to consider genome duplication, where the gene dosage is increased by two-fold at a random point in the cell-cycle. We systematically investigate how the timing of genome duplication influences different protein noise components. Intriguingly, results show that noise contribution from stochastic expression is minimized at an optimal genome duplication time. Our theoretical results motivate new experimental methods for decomposing protein noise levels from single-cell expression data. Characterizing the contributions of individual noise mechanisms will lead to precise estimates of gene expression parameters and techniques for altering stochasticity to change phenotype of individual cells.

preprint2015arXiv

Host-feeding enhances stability of discrete-time host-parasitoid population dynamic models

Discrete-time models are the traditional approach for capturing population dynamics of a host-parasitoid system. Recent work has introduced a semi-discrete framework for obtaining model update functions that connect host-parasitoid population levels from year-to-year. In particular, this framework uses differential equations to describe the hosts-parasitoid interaction during the time of year where they come in contact, allowing specific behaviors to be mechanistically incorporated into the model. We use the semi-discrete approach to study the effects of host-feeding, which occurs when a parasitoid consumes a potential host larva without ovipositing. Our results show that host-feeding by itself cannot stabilize the system, and both the host and parasitoid populations exhibit diverging oscillations similar to the Nicholson-Bailey model. However, when combined with other stabilizing mechanisms such as density-dependent host mortality or density-dependent parasitoid attack rate, host-feeding expands the region of parameter space that allows for a stable host-parasitoid equilibrium. Finally, our results show that host-feeding causes inefficiency in the parasitoid population, which yields a higher population of hosts per generation. This suggests that host-feeding may have limited long-term impact in terms of suppressing host levels for biological control applications.

preprint2015arXiv

Nonspecific transcription factor binding reduces variability in transcription factor and target protein expression

Transcription factors (TFs) interact with a multitude of binding sites on DNA and partner proteins inside cells. We investigate how nonspecific binding/unbinding to such decoy binding sites affects the magnitude and time-scale of random fluctuations in TF copy numbers arising from stochastic gene expression. A stochastic model of TF gene expression, together with decoy site interactions is formulated. Distributions for the total (bound and unbound) and free (unbound) TF levels are derived by analytically solving the chemical master equation under physiologically relevant assumptions. Our results show that increasing the number of decoy binding sides considerably reduces stochasticity in free TF copy numbers. The TF autocorrelation function reveals that decoy sites can either enhance or shorten the time-scale of TF fluctuations depending on model parameters. To understand how noise in TF abundances propagates downstream, a TF target gene is included in the model. Intriguingly, we find that noise in the expression of the target gene decreases with increasing decoy sites for linear TF-target protein dose-responses, even in regimes where decoy sites enhance TF autocorrelation times. Moreover, counterintuitive noise transmissions arise for nonlinear dose-responses. In summary, our study highlights the critical role of molecular sequestration by decoy binding sites in regulating the stochastic dynamics of TFs and target proteins at the single-cell level.

preprint2015arXiv

Stochastic Analysis Of An Incoherent Feedforward Genetic Motif

Gene products (RNAs, proteins) often occur at low molecular counts inside individual cells, and hence are subject to considerable random fluctuations (noise) in copy number over time. Not surprisingly, cells encode diverse regulatory mechanisms to buffer noise. One such mechanism is the incoherent feedforward circuit. We analyze a simplistic version of this circuit, where an upstream regulator X affects both the production and degradation of a protein Y. Thus, any random increase in X's copy numbers would increase both production and degradation, keeping Y levels unchanged. To study its stochastic dynamics, we formulate this network into a mathematical model using the Chemical Master Equation formulation. We prove that if the functional dependence of Y's production and degradation on X is similar, then the steady-distribution of Y's copy numbers is independent of X. To investigate how fluctuations in Y propagate downstream, a protein Z whose production rate only depend on Y is introduced. Intriguingly, results show that the extent of noise in Z increases with noise in X, in spite of the fact that the magnitude of noise in Y is invariant of X. Such counter intuitive results arise because X enhances the time-scale of fluctuations in Y, which amplifies fluctuations in downstream processes. In summary, while feedforward systems can buffer a protein from noise in its upstream regulators, noise can propagate downstream due to changes in the time-scale of fluctuations.

preprint2014arXiv

Moment Closure Approximations in a Genetic Negative Feedback Circuit

Auto-regulation, a process wherein a protein negatively regulates its own production, is a common motif in gene expression networks. Negative feedback in gene expression plays a critical role in buffering intracellular fluctuations in protein concentrations around optimal value. Due to the nonlinearities present in these feedbacks, moment dynamics are typically not closed, in the sense that the time derivative of the lower-order statistical moments of the protein copy number depends on high-order moments. Moment equations are closed by expressing higher-order moments as nonlinear functions of lower-order moments, a technique commonly referred to as moment closure. Here, we compare the performance of different moment closure techniques. Our results show that the commonly used closure method, which assumes a priori that the protein population counts are normally distributed, performs poorly. In contrast, conditional derivative matching, a novel closure scheme proposed here provides a good approximation to the exact moments across different parameter regimes. In summary our study provides a new moment closure method for studying stochastic dynamics of genetic negative feedback circuits, and can be extended to probe noise in more complex gene networks.

preprint2014arXiv

Optimal first-passage time in gene regulatory networks

The inherent probabilistic nature of the biochemical reactions, and low copy number of species can lead to stochasticity in gene expression across identical cells. As a result, after induction of gene expression, the time at which a specific protein count is reached is stochastic as well. Therefore events taking place at a critical protein level will see stochasticity in their timing. First-passage time (FPT), the time at which a stochastic process hits a critical threshold, provides a framework to model such events. Here, we investigate stochasticity in FPT. Particularly, we consider events for which controlling stochasticity is advantageous. As a possible regulatory mechanism, we also investigate effect of auto-regulation, where the transcription rate of gene depends on protein count, on stochasticity of FPT. Specifically, we investigate for an optimal auto-regulation which minimizes stochasticity in FPT, given fixed mean FPT and threshold. For this purpose, we model the gene expression at a single cell level. We find analytic formulas for statistical moments of the FPT in terms of model parameters. Moreover, we examine the gene expression model with auto-regulation. Interestingly, our results show that the stochasticity in FPT, for a fixed mean, is minimized when the transcription rate is independent of protein count. Further, we discuss the results in context of lysis time of an \textit{E. coli} cell infected by a $λ$ phage virus. An optimal lysis time provides evolutionary advantage to the $λ$ phage, suggesting a possible regulation to minimize its stochasticity. Our results indicate that there is no auto-regulation of the protein responsible for lysis. Moreover, congruent to experimental evidences, our analysis predicts that the expression of the lysis protein should have a small burst size.

preprint2014arXiv

Transcriptional bursting in gene expression: analytical results for general stochastic models

Gene expression in individual cells is highly variable and sporadic, often resulting in the synthesis of mRNAs and proteins in bursts. Bursting in gene expression is known to impact cell-fate in diverse systems ranging from latency in HIV-1 viral infections to cellular differentiation. It is generally assumed that bursts are geometrically distributed and that they arrive according to a Poisson process. On the other hand, recent single-cell experiments provide evidence for complex burst arrival processes, highlighting the need for more general stochastic models. To address this issue, we invoke a mapping between general models of gene expression and systems studied in queueing theory to derive exact analytical expressions for the moments associated with mRNA/protein steady-state distributions. These moments are then used to derive explicit conditions, based entirely on experimentally measurable quantities, that determine if the burst distributions deviate from the geometric distribution or if burst arrival deviates from a Poisson process. For non-Poisson arrivals, we develop approaches for accurate estimation of burst parameters.