Source author record

Michael C. Mackey

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

ResearcherUnclaimed source record

Molecular Networks math.PR nlin.CD math-ph math.CA math.DS math.MP Subcellular Processes

Catalog footprint

What is connected

7works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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

Building this map preview

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

preprint2021arXiv

Operon Dynamics with State Dependent Transcription and/or Translation Delays

Transcription and translation retrieve and operationalize gene encoded information in cells. These processes are not instantaneous and incur significant delays. In this paper we study Goodwin models of both inducible and repressible operons with state-dependent delays. The paper provides justification and derivation of the model, detailed analysis of the appropriate setting of the corresponding dynamical system, and extensive numerical analysis of its dynamics. Comparison with constant delay models shows significant differences in dynamics that include existence of stable periodic orbits in inducible systems and multistability in repressible systems. A combination of parameter space exploration, numerics, analysis of steady state linearization and bifurcation theory indicates the likely presence of Shilnikov-type homoclinic bifurcations in the repressible operon model.

preprint2016arXiv

Mathematical model of galactose regulation and metabolic consumption in yeast

The galactose network is a complex system responsible for galactose metabolism. It has been extensively studied experimentally and mathematically at the unicellular level to broaden our understanding of its regulatory mechanisms at higher order species. Although the key molecular players involved in the metabolic and regulatory processes underlying this system have been known for decades, their interactions and chemical kinetics remain incompletely understood. Mathematical models can provide an alternative method to study the dynamics of this network from a quantitative and a qualitative perspective. Here, we employ such approaches to unravel the main properties of the galactose network, including equilibrium binary and temporal responses, as a way to decipher its adaptation to actively-changing inputs. We combine the two main components of the network; namely, the genetic branch which allows for bistable responses, and a metabolic branch, encompassing the relevant metabolic processes and glucose repressive reactions. We use both computational tools to estimate model parameters based on published experimental data, as well as bifurcation analysis to decipher the properties of the system in various parameter regimes. Our model analysis reveals that the interplay between the inducer (galactose) and the repressor (glucose) creates the bistability regime which dictates the temporal responses of the system. Based on the same bifurcation techniques, we can also explain why the system is robust to genetic mutations and molecular instabilities. These findings may provide experimentalists with a theoretical framework upon which they can determine how the galactose network functions under various conditions.

preprint2015arXiv

Random Numbers from a Delay Equation

Delay differential equations (DDE) can have "chaotic" solutions that can be used to mimic Brownian motion. Since a Brownian motion is random in its velocity, it is reasonable to think that a random number generator (RNG) might be constructed from such a model. In this preliminary study, we consider one specific example of this and show that it satisfies criteria commonly employed in the testing of random number generators (from TestU01's very stringent "Big Crush" battery of tests). A technique termed digit discarding, commonly used in both this generator and physical RNG's using laser feedback systems, is discussed with regard to the maximal Lyapunov exponent. Also, we benchmark the generator to a contemporary common method: the multiple recursive generator, MRG32k3a. Although our method is about 7 times slower than MRG32k3a, there is in principle no apparent limit on the number of possible values that can be generated from the scheme we present here.

preprint2015arXiv

The Limiting Dynamics of a Bistable Molecular Switch With and Without Noise

We consider the dynamics of a population of organisms containing two mutually inhibitory gene regulatory networks, that can result in a bistable switch-like behaviour. We completely characterize their local and global dynamics in the absence of any noise, and then go on to consider the effects of either noise coming from bursting (transcription or translation), or Gaussian noise in molecular degradation rates when there is a dominant slow variable in the system. We show analytically how the steady state distribution in the population can range from a single unimodal distribution through a bimodal distribution and give the explicit analytic form for the invariant stationary density which is globally asymptotically stable. Rather remarkably, the behaviour of the stationary density with respect to the parameters characterizing the molecular behaviour of the bistable switch is qualitatively identical in the presence of noise coming from bursting as well as in the presence of Gaussian noise in the degradation rate. This implies that one cannot distinguish between either the dominant source or nature of noise based on the stationary molecular distribution in a population of cells. We finally show that the switch model with bursting but two dominant slow genes has an asymptotically stable stationary density.

preprint2013arXiv

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

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

preprint2011arXiv

Deterministic Brownian motion generated from differential delay equations

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

preprint2010arXiv

Molecular Distributions in Gene Regulatory Dynamics

We show how one may analytically compute the stationary density of the distribution of molecular constituents in populations of cells in the presence of noise arising from either bursting transcription or translation, or noise in degradation rates arising from low numbers of molecules. We have compared our results with an analysis of the same model systems (either inducible or repressible operons) in the absence of any stochastic effects, and shown the correspondence between behaviour in the deterministic system and the stochastic analogs. We have identified key dimensionless parameters that control the appearance of one or two steady states in the deterministic case, or unimodal and bimodal densities in the stochastic systems, and detailed the analytic requirements for the occurrence of different behaviours. This approach provides, in some situations, an alternative to computationally intensive stochastic simulations. Our results indicate that, within the context of the simple models we have examined, bursting and degradation noise cannot be distinguished analytically when present alone.