BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Andrew Mugler

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

ResearcherUnclaimed source record
Molecular NetworksBiological PhysicsCell BehaviorQuantitative MethodsSubcellular Processescond-mat.stat-mechnlin.AO

Catalog footprint

What is connected

27works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

27 published item(s)

preprint2026arXiv

Optimal information transmission in a sequential model for cell division

In proliferating cell populations, adaptive changes to biochemical reactions can change a cell's division time, which in turn can change the population size. However, biochemical reactions are subject to noise, and therefore the conditions for optimal information transmission from the molecular to the population scale are poorly understood. Here, we model cell proliferation as a Bellman-Harris branching process with age-dependent division times. We identify a class of division time distributions, built from a series of Markovian steps, for which the population size distribution at all times is hierarchically calculable. We use this feature to characterize the amount of influence that a given reaction step has on the population size via the mutual information. We find that information transmission is optimal for a characteristic number of steps until division: too few and the population size is unpredictable; too many and any given step has vanishing influence on the population size. Our work reveals the potential tradeoffs involved in adaptive decision making at the sub-cellular, cellular and population scales.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Restoring information in aged gene regulatory networks by single knock-ins

A hallmark of aging is loss of information in gene regulatory networks. These networks are tightly connected, raising the question of whether information could be restored by perturbing single genes. We develop a simple theoretical framework for information transmission in gene regulatory networks that describes the information gained or lost when a gene is "knocked in" (exogenously expressed). Applying the framework to gene expression data from muscle cells in young and old mice, we find that single knock-ins can restore network information by up to 10%. Our work advances the study of information flow in networks and identifies potential gene targets for rejuvenation.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Autologous chemotaxis at high cell density

Autologous chemotaxis, in which cells secrete and detect molecules to determine the direction of fluid flow, is thwarted at high cell density because molecules from other cells interfere with a given cell's signal. Using a minimal model of autologous chemotaxis, we determine the cell density at which sensing fails and find that it agrees with experimental observations of metastatic cancer cells. To understand this agreement, we derive a physical limit to autologous chemotaxis in terms of the cell density, the Péclet number, and the length scales of the cell and its environment. Surprisingly, in an environment that is uniformly oversaturated in the signaling molecule, we find that sensing not only can fail, but can be reversed, causing backwards cell motion. Our results get to the heart of the competition between chemical and mechanical cellular sensing and shed light on a sensory strategy employed by cancer cells in dense tumor environments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Cell-to-cell information at a feedback-induced bifurcation point

A ubiquitous way that cells share information is by exchanging molecules. Yet, the fundamental ways that this information exchange is influenced by intracellular dynamics remain unclear. Here we use information theory to investigate a simple model of two interacting cells with internal feedback. We show that cell-to-cell molecule exchange induces a collective two-cell critical point and that the mutual information between the cells peaks at this critical point. Information can remain large far from the critical point on a manifold of cellular states, but scales logarithmically with the correlation time of the system, resulting in an information-correlation time tradeoff. This tradeoff is strictly imposed, suggesting the correlation time as a proxy for the mutual information.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Precision of flow sensing by self-communicating cells

Metastatic cancer cells detect the direction of lymphatic flow by self-communication: they secrete and detect a chemical which, due to the flow, returns to the cell surface anisotropically. The secretion rate is low, meaning detection noise may play an important role, but the sensory precision of this mechanism has not been explored. Here we derive the precision of flow sensing for two ubiquitous detection methods: absorption vs.\ reversible binding to surface receptors. We find that binding is more precise due to the fact that absorption distorts the signal that the cell aims to detect. Comparing to experiments, our results suggest that the cancer cells operate remarkably close to the physical detection limit. Our prediction that cells should bind the chemical reversibly, not absorb it, is supported by endocytosis data for this ligand-receptor pair.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Spiral wave propagation in communities with spatially correlated heterogeneity

Many multicellular communities propagate signals in a directed manner via excitable waves. Cell-to-cell heterogeneity is a ubiquitous feature of multicellular communities, but the effects of heterogeneity on wave propagation are still unclear. Here we use a minimal FitHugh-Nagumo-type model to investigate excitable wave propagation in a two-dimensional heterogeneous community. The model shows three dynamic regimes in which waves either propagate directionally, die out, or spiral indefinitely, and we characterize how these regimes depend on the heterogeneity parameters. We find that in some parameter regimes, spatial correlations in the heterogeneity enhance directional propagation and suppress spiraling. However, in other regimes, spatial correlations promote spiraling, a surprising feature that we explain by demonstrating that these spirals form by a second, distinct mechanism. Finally, we characterize the dependence of the spiral period on the degree of heterogeneity in the system by using techniques from percolation theory. Our results reveal that the spatial structure of cell-to-cell heterogeneity can have important consequences for signal propagation in cellular communities.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Statistics of correlated percolation in a bacterial community

Signal propagation over long distances is a ubiquitous feature of multicellular communities. In biofilms of the bacterium Bacillus subtilis, we recently discovered that some, but not all, cells participate in the propagation of an electrical signal, and the ones that do are organized in a way that is statistically consistent with percolation theory. However, two key assumptions of percolation theory are violated in this system. First, we find here that the probability for a cell to signal is not independent from other cells but instead is correlated with its nearby neighbors. We develop a mechanistic model, in which correlated signaling emerges from cell division, phenotypic inheritance, and cell displacement, that reproduces the experimental results. Second, we observed previously that the fraction of signaling cells is not constant but instead varies from biofilm to biofilm. We use our model to understand why percolation theory remains a valid description of the system despite these two violations of its assumptions. We find that the first violation does not significantly affect the spatial statistics, which we rationalize using a renormalization argument. We find that the second violation widens the range of signaling fractions around the percolation threshold at which one observes the characteristic power-law statistics of cluster sizes, consistent with our previous experimental results. We validate our model using a mutant biofilm whose signaling probability decays along the propagation direction. Our results identify key statistical features of a correlated percolating system and demonstrate their functional utility for a multicellular community.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Temporal precision of molecular events with regulation and feedback

Cellular behaviors such as migration, division, and differentiation rely on precise timing, and yet the molecular events that govern these behaviors are highly stochastic. We investigate regulatory strategies that decrease the timing noise of molecular events. Autoregulatory feedback increases noise. Yet, we find that in the presence of regulation by a second species, autoregulatory feedback decreases noise. To explain this finding, we develop a method to calculate the optimal regulation function that minimizes the timing noise. The method reveals that the combination of feedback and regulation minimizes noise by maximizing the number of molecular events that must happen in sequence before a threshold is crossed. We compute the optimal timing precision for all two-node networks with regulation and feedback, derive a generic lower bound on timing noise, and discuss our results in the context of neuroblast migration during Caenorhabditis elegans development.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Collective chemotaxis through noisy multicellular gradient sensing

Collective cell migration in response to a chemical cue occurs in many biological processes such as morphogenesis and cancer metastasis. Clusters of migratory cells in these systems are capable of responding to gradients of less than 1% difference in chemical concentration across a cell length. Multicellular systems are extremely sensitive to their environment and while the limits to multicellular sensing are becoming known, how this information leads to coherent migration remains poorly understood. We develop a computational model of multicellular sensing and migration in which groups of cells collectively measure noisy chemical gradients. The output of the sensing process is coupled to individual cells polarization to model migratory behavior. Through the use of numerical simulations, we find that larger clusters of cells detect the gradient direction with higher precision and thus achieve stronger polarization bias, but larger clusters also induce more drag on collective motion. The trade-off between these two effects leads to an optimal cluster size for most efficient migration. We discuss how our model could be validated using simple, phenomenological experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Sense and sensitivity: physical limits to multicellular sensing, migration and drug response

Metastasis is a process of cell migration that can be collective and guided by chemical cues. Viewing metastasis in this way, as a physical phenomenon, allows one to draw upon insights from other studies of collective sensing and migration in cell biology. Here we review recent progress in the study of cell sensing and migration as collective phenomena, including in the context of metastatic cells. We describe simple physical models that yield the limits to the precision of cell sensing, and we review experimental evidence that cells operate near these limits. Models of collective migration are surveyed in order understand how collective metastatic invasion can occur. We conclude by contrasting cells' sensory abilities with their sensitivity to drugs, and suggesting potential alternatives to cell-death-based cancer therapies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cell-cell communication enhances the capacity of cell ensembles to sense shallow gradients during morphogenesis

Collective cell responses to exogenous cues depend on cell-cell interactions. In principle, these can result in enhanced sensitivity to weak and noisy stimuli. However, this has not yet been shown experimentally, and, little is known about how multicellular signal processing modulates single cell sensitivity to extracellular signaling inputs, including those guiding complex changes in the tissue form and function. Here we explored if cell-cell communication can enhance the ability of cell ensembles to sense and respond to weak gradients of chemotactic cues. Using a combination of experiments with mammary epithelial cells and mathematical modeling, we find that multicellular sensing enables detection of and response to shallow Epidermal Growth Factor (EGF) gradients that are undetectable by single cells. However, the advantage of this type of gradient sensing is limited by the noisiness of the signaling relay, necessary to integrate spatially distributed ligand concentration information. We calculate the fundamental sensory limits imposed by this communication noise and combine them with the experimental data to estimate the effective size of multicellular sensory groups involved in gradient sensing. Functional experiments strongly implicated intercellular communication through gap junctions and calcium release from intracellular stores as mediators of collective gradient sensing. The resulting integrative analysis provides a framework for understanding the advantages and limitations of sensory information processing by relays of chemically coupled cells.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Limits to the precision of gradient sensing with spatial communication and temporal integration

Gradient sensing requires at least two measurements at different points in space. These measurements must then be communicated to a common location to be compared, which is unavoidably noisy. While much is known about the limits of measurement precision by cells, the limits placed by the communication are not understood. Motivated by recent experiments, we derive the fundamental limits to the precision of gradient sensing in a multicellular system, accounting for communication and temporal integration. The gradient is estimated by comparing a "local" and a "global" molecular reporter of the external concentration, where the global reporter is exchanged between neighboring cells. Using the fluctuation-dissipation framework, we find, in contrast to the case when communication is ignored, that precision saturates with the number of cells independently of the measurement time duration, since communication establishes a maximum lengthscale over which sensory information can be reliably conveyed. Surprisingly, we also find that precision is improved if the local reporter is exchanged between cells as well, albeit more slowly than the global reporter. The reason is that while exchange of the local reporter weakens the comparison, it decreases the measurement noise. We term such a model "regional excitation--global inhibition" (REGI). Our results demonstrate that fundamental sensing limits are necessarily sharpened when the need to communicate information is taken into account.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Noise expands the response range of the Bacillus subtilis competence circuit

Gene regulatory circuits must contend with intrinsic noise that arises due to finite numbers of proteins. While some circuits act to reduce this noise, others appear to exploit it. A striking example is the competence circuit in Bacillus subtilis, which exhibits much larger noise in the duration of its competence events than a synthetically constructed analog that performs the same function. Here, using stochastic modeling and fluorescence microscopy, we show that this larger noise allows cells to exit terminal phenotypic states, which expands the range of stress levels to which cells are responsive and leads to phenotypic heterogeneity at the population level. This is an important example of how noise confers a functional benefit in a genetic decision-making circuit.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Optimal Prediction by Cellular Signaling Networks

Living cells can enhance their fitness by anticipating environmental change. We study how accurately linear signaling networks in cells can predict future signals. We find that maximal predictive power results from a combination of input-noise suppression, linear extrapolation, and selective readout of correlated past signal values. Single-layer networks generate exponential response kernels, which suffice to predict Markovian signals optimally. Multilayer networks allow oscillatory kernels that can optimally predict non-Markovian signals. At low noise, these kernels exploit the signal derivative for extrapolation, while at high noise, they capitalize on signal values in the past that are strongly correlated with the future signal. We show how the common motifs of negative feedback and incoherent feed-forward can implement these optimal response functions. Simulations reveal that E. coli can reliably predict concentration changes for chemotaxis, and that the integration time of its response kernel arises from a trade-off between rapid response and noise suppression.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Role of spatial averaging in multicellular gradient sensing

Gradient sensing underlies important biological processes including morphogenesis, polarization, and cell migration. The precision of gradient sensing increases with the length of a detector (a cell or group of cells) in the gradient direction, since a longer detector spans a larger range of concentration values. Intuition from analyses of concentration sensing suggests that precision should also increase with detector length in the direction transverse to the gradient, since then spatial averaging should reduce the noise. However, here we show that, unlike for concentration sensing, the precision of gradient sensing decreases with transverse length for the simplest gradient sensing model, local excitation--global inhibition (LEGI). The reason is that gradient sensing ultimately relies on a subtraction of measured concentration values. While spatial averaging indeed reduces the noise in these measurements, which increases precision, it also reduces the covariance between the measurements, which results in the net decrease in precision. We demonstrate how a recently introduced gradient sensing mechanism, regional excitation--global inhibition (REGI), overcomes this effect and recovers the benefit of transverse averaging. Using a REGI-based model, we compute the optimal two- and three-dimensional detector shapes, and argue that they are consistent with the shapes of naturally occurring gradient-sensing cell populations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stochastic modeling of gene expression, protein modification, and polymerization

Many fundamental cellular processes involve small numbers of molecules. When numbers are small, fluctuations dominate, and stochastic models, which account for these fluctuations, are required. In this chapter, we describe minimal stochastic models of three fundamental cellular processes: gene expression, protein modification, and polymerization. We introduce key analytic tools for solving each model, including the generating function, eigenfunction expansion, and operator methods, and we discuss how these tools are extended to more complicated models. These analytic tools provide an elegant, efficient, and often insightful alternative to stochastic simulation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Prediction and Dissipation in Biochemical Sensing

Cells sense and predict their environment via energy-dissipating pathways. However, it is unclear whether dissipation helps or harms prediction. Here we study dissipation and prediction for a minimal sensory module of receptors that reversibly bind ligand. We find that the module performs short-term prediction optimally when operating in an adiabatic regime where dissipation vanishes. In contrast, beyond a critical forecast interval, prediction becomes most precise in a regime of maximal dissipation, suggesting that dissipative sensing in biological systems can serve to enhance prediction performance.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spatial partitioning improves the reliability of biochemical signaling

Spatial heterogeneity is a hallmark of living systems, even at the molecular scale in individual cells. A key example is the partitioning of membrane-bound proteins via lipid domain formation or cytoskeleton-induced corralling. Yet the impact of this spatial heterogeneity on biochemical signaling processes is poorly understood. Here we demonstrate that partitioning improves the reliability of biochemical signaling. We exactly solve a stochastic model describing a ubiquitous motif in membrane signaling. The solution reveals that partitioning improves signaling reliability via two effects: it moderates the non-linearity of the switching response, and it reduces noise in the response by suppressing correlations between molecules. An optimal partition size arises from a trade-off between minimizing the number of proteins per partition to improve signaling reliability and ensuring sufficient proteins per partition to maintain signal propagation. The predicted optimal partition size agrees quantitatively with experimentally observed systems. These results persist in spatial simulations with explicit diffusion barriers. Our findings suggest that molecular partitioning is not merely a consequence of the complexity of cellular substructures, but also plays an important functional role in cell signaling.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The importance of crowding in signaling, genetic, and metabolic networks

It is now well established that the cell is a highly crowded environment. Yet, the effects of crowding on the dynamics of signaling pathways, gene regulation networks and metabolic networks are still largely unknown. Crowding can alter both molecular diffusion and the equilibria of biomolecular reactions. In this review, we first discuss how diffusion can affect biochemical networks. Diffusion of transcription factors can increase noise in gene expression, while diffusion of proteins between intracellular compartments or between cells can reduce concentration fluctuations. In push-pull networks diffusion can impede information transmission, while in multi-site protein modification networks diffusion can qualitatively change the macroscopic response of the system, such as the loss or emergence of bistability. Moreover, diffusion can directly change the metabolic flux. We describe how crowding affects diffusion, and thus how all these phenomena are influenced by crowding. Yet, a potentially more important effect of crowding on biochemical networks is mediated via the shift in the equilibria of bimolecular reactions, and we provide computational evidence that supports this idea. Finally, we discuss how the effects of crowding can be incorporated in models of biochemical networks.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Protein logic: a statistical mechanical study of signal integration at the single-molecule level

Information processing and decision making is based upon logic operations, which in cellular networks has been well characterized at the level of transcription. In recent years however, both experimentalists and theorists have begun to appreciate that cellular decision making can also be performed at the level of a single protein, giving rise to the notion of protein logic. Here we systematically explore protein logic using a well known statistical mechanical model. As an example system, we focus on receptors which bind either one or two ligands, and their associated dimers. Notably, we find that a single heterodimer can realize any of the 16 possible logic gates, including the XOR gate, by variation of biochemical parameters. We then introduce the novel idea that a set of receptors with fixed parameters can encode functionally unique logic gates simply by forming different dimeric combinations. An exhaustive search reveals that the simplest set of receptors (two single-ligand receptors and one double-ligand receptor) can realize several different groups of three unique gates, a result for which the parametric analysis of single receptors and dimers provides a clear interpretation. Both results underscore the surprising functional freedom readily available to cells at the single-protein level.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The macroscopic effects of microscopic heterogeneity

Over the past decade, advances in super-resolution microscopy and particle-based modeling have driven an intense interest in investigating spatial heterogeneity at the level of single molecules in cells. Remarkably, it is becoming clear that spatiotemporal correlations between just a few molecules can have profound effects on the signaling behavior of the entire cell. While such correlations are often explicitly imposed by molecular structures such as rafts, clusters, or scaffolds, they also arise intrinsically, due strictly to the small numbers of molecules involved, the finite speed of diffusion, and the effects of macromolecular crowding. In this chapter we review examples of both explicitly imposed and intrinsic correlations, focusing on the mechanisms by which microscopic heterogeneity is amplified to macroscopic effect.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Membrane clustering and the role of rebinding in biochemical signaling

In many cellular signaling pathways, key components form clusters at the cell membrane. Although much work has focused on the mechanisms behind such cluster formation, the implications for downstream signaling remain poorly understood. Here, motivated by recent experiments, we study via particle-based simulation a covalent modification network in which the activating component is either clustered or randomly distributed on the membrane. We find that while clustering reduces the response of a single-modification network, clustering can enhance the response of a double-modification network. The reduction is a bulk effect: a cluster presents a smaller effective target to a substrate molecule in the bulk. The enhancement, on the other hand, is a local effect: a cluster promotes the rapid rebinding and second activation of singly active substrate molecules. As such, the enhancement relies upon frequent collisions on a short timescale, which leads to a diffusion coefficient at which the enhancement is optimal. We complement simulation with analytic results at both the mean-field and first-passage distribution levels. Our results emphasize the importance of spatially resolved models, showing that significant effects of spatial correlations persist even in spatially averaged quantities such as response curves.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A statistical method for revealing form-function relations in biological networks

Over the past decade, a number of researchers in systems biology have sought to relate the function of biological systems to their network-level descriptions -- lists of the most important players and the pairwise interactions between them. Both for large networks (in which statistical analysis is often framed in terms of the abundance of repeated small subgraphs) and for small networks which can be analyzed in greater detail (or even synthesized in vivo and subjected to experiment), revealing the relationship between the topology of small subgraphs and their biological function has been a central goal. We here seek to pose this revelation as a statistical task, illustrated using a particular setup which has been constructed experimentally and for which parameterized models of transcriptional regulation have been studied extensively. The question "how does function follow form" is here mathematized by identifying which topological attributes correlate with the diverse possible information-processing tasks which a transcriptional regulatory network can realize. The resulting method reveals one form-function relationship which had earlier been predicted based on analytic results, and reveals a second for which we can provide an analytic interpretation. Resulting source code is distributed via http://formfunction.sourceforge.net.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Analytic methods for modeling stochastic regulatory networks

The past decade has seen a revived interest in the unavoidable or intrinsic noise in biochemical and genetic networks arising from the finite copy number of the participating species. That is, rather than modeling regulatory networks in terms of the deterministic dynamics of concentrations, we model the dynamics of the probability of a given copy number of the reactants in single cells. Most of the modeling activity of the last decade has centered on stochastic simulation of individual realizations, i.e., Monte-Carlo methods for generating stochastic time series. Here we review the mathematical description in terms of probability distributions, introducing the relevant derivations and illustrating several cases for which analytic progress can be made either instead of or before turning to numerical computation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Telling time with an intrinsically noisy clock

Intracellular transmission of information via chemical and transcriptional networks is thwarted by a physical limitation: the finite copy number of the constituent chemical species introduces unavoidable intrinsic noise. Here we provide a method for solving for the complete probabilistic description of intrinsically noisy oscillatory driving. We derive and numerically verify a number of simple scaling laws. Unlike in the case of measuring a static quantity, response to an oscillatory driving can exhibit a resonant frequency which maximizes information transmission. Further, we show that the optimal regulatory design is dependent on the biophysical constraints (i.e., the allowed copy number and response time). The resulting phase diagram illustrates under what conditions threshold regulation outperforms linear regulation.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Network motifs come in sets: correlations in the randomization process

The identification of motifs--subgraphs that appear significantly more often in a particular network than in an ensemble of randomized networks--has become a ubiquitous method for uncovering potentially important subunits within networks drawn from a wide variety of fields. We find that the most common algorithms used to generate the ensemble from the real network change subgraph counts in a highly correlated manner, so that one subgraph's status as a motif may not be independent from the statuses of the other subgraphs. We demonstrate this effect for the problem of 3- and 4-node motif identification in the transcriptional regulatory networks of E. coli and S. cerevisiae in which randomized networks are generated via an edge-swapping algorithm (Milo et al., Science 298:824, 2002). We show that correlations among 3-node subgraphs are easily interpreted, and we present an information-theoretic tool that may be used to identify correlations among subgraphs of any size.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Statistical properties of multistep enzyme-mediated reactions

Enzyme-mediated reactions may proceed through multiple intermediate conformational states before creating a final product molecule, and one often wishes to identify such intermediate structures from observations of the product creation. In this paper, we address this problem by solving the chemical master equations for various enzymatic reactions. We devise a perturbation theory analogous to that used in quantum mechanics that allows us to determine the first (<n>) and the second (variance) cumulants of the distribution of created product molecules as a function of the substrate concentration and the kinetic rates of the intermediate processes. The mean product flux V=d<n>/dt (or "dose-response" curve) and the Fano factor F=variance/<n> are both realistically measurable quantities, and while the mean flux can often appear the same for different reaction types, the Fano factor can be quite different. This suggests both qualitative and quantitative ways to discriminate between different reaction schemes, and we explore this possibility in the context of four sample multistep enzymatic reactions. We argue that measuring both the mean flux and the Fano factor can not only discriminate between reaction types, but can also provide some detailed information about the internal, unobserved kinetic rates, and this can be done without measuring single-molecule transition events.

0 citations0 reviewsNo code linkedOpen access