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

Tom Chou

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

ResearcherUnclaimed source record
Quantitative Methodscond-mat.stat-mechPopulations and EvolutionBiomoleculescond-mat.softSubcellular ProcessesCell Behaviorcond-mat.mes-hallmath.OCphysics.data-anTissues and OrgansApplicationsBiological PhysicsComputer Visioncond-mat.mtrl-scihep-thmath-phmath.APmath.MPmath.NAmath.PRmath.STNeurons and CognitionNumerical Analysisphysics.chem-phphysics.comp-phStatistics Theory

Catalog footprint

What is connected

19works
27topics
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

19 published item(s)

preprint2022arXiv

Adaptive Hermite Spectral Methods in Unbounded Domains

Recently, new adaptive techniques were developed that greatly improved the efficiency of solving PDEs using spectral methods. These adaptive spectral techniques are especially suited for accurately solving problems in unbounded domains and require the monitoring and dynamic adjustment of three key tunable parameters: the scaling factor, the displacement of the basis functions, and the spectral expansion order. There have been few analyses of numerical methods for unbounded domain problems. Specifically, there is no analysis of adaptive spectral methods to provide insight into how to increase efficiency and accuracy through dynamical adjustment of parameters. In this paper, we perform the first numerical analysis of the adaptive spectral method using generalized Hermite functions in both one- and multi-dimensional problems. Our analysis reveals why adaptive spectral methods work well when a "frequency indicator" of the numerical solution is controlled. We then investigate how the implementation of the adaptive spectral methods affects numerical results, thereby providing guidelines for the proper tuning of parameters. Finally, we further improve performance by extending the adaptive methods to allow bidirectional basis function translation, and the prospect of carrying out similar numerical analysis to solving PDEs arising from realistic difficult-to-solve unbounded models with adaptive spectral methods is also briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Density- and elongation speed-dependent error correction in RNA polymerization

Backtracking of RNA polymerase (RNAP) is an important pausing mechanism during DNA transcription that is part of the error correction process that enhances transcription fidelity. We model the backtracking mechanism of RNA polymerase, which usually happens when the polymerase tries to incorporate a mismatched nucleotide triphosphate. Previous models have made simplifying assumptions such as neglecting the trailing polymerase behind the backtracking polymerase or assuming that the trailing polymerase is stationary. We derive exact analytic solutions of a stochastic model that includes locally interacting RNAPs by explicitly showing how a trailing RNAP influences the probability that an error is corrected or incorporated by the leading backtracking RNAP. We also provide two related methods for computing the mean times to error correction or incorporation given an initial local RNAP configuration.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Diversity in Biology: definitions, quantification, and models

Diversity indices are useful single-number metrics for characterizing a complex distribution of a set of attributes across a population of interest. The utility of these different metrics or sets of metrics depend on the context and application, and whether a predictive mechanistic model exists. In this topical review, we first summarize the relevant mathematical principles underlying heterogeneity in a large population before outlining the various definitions of `diversity' and providing examples of scientific topics in which its quantification plays an important role. We then review how diversity has been a ubiquitous concept across multiple fields including ecology, immunology, cellular barcoding experiments, and socioeconomic studies. Since many of these applications involve sampling of populations, we also review how diversity in small samples is related to the diversity in the entire population. Features that arise in each of these applications are highlighted.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

How heterogeneous thymic output and homeostatic proliferation shape naive T cell receptor clone abundance distributions

The set of T cells that express the same T cell receptor (TCR) sequence represents a T cell clone. The number of different naive T cell clones in an organism reflects the number of different T cell receptors (TCRs) arising from recombination of the V(D)J gene segments during T cell development in the thymus. TCR diversity and more specifically, the clone abundance distribution is an important factor in immune function. Specific recombination patterns occur more frequently than others while subsequent interactions between TCRs and self-antigens are known to trigger proliferation and sustain naive T cell survival. These processes are TCR-dependent, leading to clone-dependent thymic export and naive T cell proliferation rates. Using a mean-field approximation to the solution of a regulated birth-death-immigration model and a modification arising from sampling, we systematically quantify how TCR-dependent heterogeneities in immigration and proliferation rates affect the shape of clone abundance distributions (the number of different clones that are represented by a specific number of cells, or "clone counts"). By comparing predicted clone counts derived from our heterogeneous birth-death-immigration model with experimentally sampled clone abundances, we show that although heterogeneity in immigration rates causes very little change to predicted clone-counts, significant heterogeneity in proliferation rates is necessary to generate the observed abundances with reasonable physiological parameter values. Our analysis provides constraints among physiological parameters that are necessary to yield predictions that qualitatively match the data. Assumptions of the model and potentially other important mechanistic factors are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

PDE models of adder mechanisms in cellular proliferation

Cell division is a process that involves many biochemical steps and complex biophysical mechanisms. To simplify the understanding of what triggers cell division, three basic models that subsume more microscopic cellular processes associated with cell division have been proposed. Cells can divide based on the time elapsed since their birth, their size, and/or the volume added since their birth -- the timer, sizer, and adder models, respectively. Here, we propose unified adder-sizer models and investigate some of the properties of different adder processes arising in cellular proliferation. Although the adder-sizer model provides a direct way to model cell population structure, we illustrate how it is mathematically related to the well-known model in which cell division depends on age and size. Existence and uniqueness of weak solutions to our 2+1-dimensional PDE model are proved, leading to the convergence of the discretized numerical solutions and allowing us to numerically compute the dynamics of cell population densities. We then generalize our PDE model to incorporate recent experimental findings of a system exhibiting mother-daughter correlations in cellular growth rates. Numerical experiments illustrating possible average cell volume blowup and the dynamical behavior of cell populations with mother-daughter correlated growth rates are carried out. Finally, motivated by new experimental findings, we extend our adder model cases where the controlling variable is the added size between DNA replication initiation points in the cell cycle.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Uniformly accurate effective equations for disease transmission mediated by pair formation dynamics

We derive and asymptotically analyze mass-action models for disease spread that include transient pair formation and dissociation. Populations of unpaired susceptibles and infecteds are distinguished from the population of three types of pairs of individuals; both susceptible, one susceptible and one infected, and both infected. Disease transmission can occur only within a pair consisting of one susceptible individual and one infected individual. By considering the fast pair formation and fast pair dissociation limits, we use a perturbation expansion to formally derive a uniformly valid approximation for the dynamics of the total infected and susceptible populations. Under different parameter regimes, we derive uniformly valid effective equations for the total infected population and compare their results to those of the full mass-action model. Our results are derived from the fundamental mass-action system without implicitly imposing transmission mechanisms such as that used in frequency-dependent models. They provide a new formulation for effective pairing models and are compared with previous models.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A hierarchical kinetic theory of birth, death, and fission in age-structured interacting populations

We study mathematical models describing the evolution of stochastic age-structured populations. After reviewing existing approaches, we present a full kinetic framework for age-structured interacting populations undergoing birth, death and fission processes, in spatially dependent environments. We define the complete probability density for the population-size-age-chart and find results under specific conditions. Connections with more classical models are also explicitly derived. In particular, we show that factorial moments for non-interacting processes are described by a natural generalization of the McKendrick-von Foerster equation, which describes mean-field deterministic behaviour. Our approach utilizes mixed type, multi-dimensional probability distributions similar to those employed in the study of gas kinetics, with terms that satisfy BBGKY-like equation hierarchies.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Onset, timing, and exposure therapy of stress disorders: mechanistic insight from a mathematical model of oscillating neuroendocrine dynamics

The hypothalamic-pituitary-adrenal (HPA) axis is a neuroendocrine system that regulates numerous physiological processes. Disruptions in the activity of the HPA axis are correlated with many stress-related diseases such as post-traumatic stress disorder (PTSD) and major depressive disorder. In this paper, we characterize "normal" and "diseased" states of the HPA axis as basins of attraction of a dynamical system describing the inhibition of peptide hormones such as corticotropin-releasing hormone (CRH) and adrenocorticotropic hormone (ACTH) by circulating glucocorticoids such as cortisol (CORT). In addition to including key physiological features such as ultradian oscillations in cortisol levels and self-upregulation of CRH neuron activity, our model distinguishes the relatively slow process of cortisol-mediated CRH biosynthesis from rapid trans-synaptic effects that regulate the CRH secretion process. Crucially, we find that the slow regulation mechanism mediates external stress-driven transitions between the stable states in novel, intensity, duration, and timing-dependent ways. These results indicate that the timing of traumatic events may be an important factor in determining if and how patients will exhibit hallmarks of stress disorders. Our model also suggests a mechanism whereby exposure therapy of stress disorders such as PTSD may act to normalize downstream dysregulation of the HPA axis.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A path-integral approach to Bayesian inference for inverse problems using the semiclassical approximation

We demonstrate how path integrals often used in problems of theoretical physics can be adapted to provide a machinery for performing Bayesian inference in function spaces. Such inference comes about naturally in the study of inverse problems of recovering continuous (infinite dimensional) coefficient functions from ordinary or partial differential equations (ODE, PDE), a problem which is typically ill-posed. Regularization of these problems using $L^2$ function spaces (Tikhonov regularization) is equivalent to Bayesian probabilistic inference, using a Gaussian prior. The Bayesian interpretation of inverse problem regularization is useful since it allows one to quantify and characterize error and degree of precision in the solution of inverse problems, as well as examine assumptions made in solving the problem -- namely whether the subjective choice of regularization is compatible with prior knowledge. Using path-integral formalism, Bayesian inference can be explored through various perturbative techniques, such as the semiclassical approximation, which we use in this manuscript. Perturbative path-integral approaches, while offering alternatives to computational approaches like Markov-Chain-Monte-Carlo (MCMC), also provide natural starting points for MCMC methods that can be used to refine approximations. In this manuscript, we illustrate a path-integral formulation for inverse problems and demonstrate it on an inverse problem in membrane biophysics as well as inverse problems in potential theories involving the Poisson equation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fixation times in differentiation and evolution in the presence of bottlenecks, deserts, and oases

Cellular differentiation and evolution are stochastic processes that can involve multiple types (or states) of particles moving on a complex, high-dimensional state-space or "fitness" landscape. Cells of each specific type can thus be quantified by their population at a corresponding node within a network of states. Their dynamics across the state-space network involve genotypic or phenotypic transitions that can occur upon cell division, such as during symmetric or asymmetric cell differentiation, or upon spontaneous mutation. Waiting times between transitions can be nonexponentially distributed and reflect e.g., the cell cycle. Here, we use a multi-type branching processes to study first passage time statistics for a single cell to appear in a specific state. We present results for a sequential evolutionary process in which $L$ successive transitions propel a population from a "wild-type" state to a given "terminally differentiated," "resistant," or "cancerous" state. Analytic and numeric results are also found for first passage times across an evolutionary chain containing a node with increased death or proliferation rate, representing a desert/bottleneck or an oasis. Processes involving cell proliferation are shown to be "nonlinear" (even though mean-field equations for the expected particle numbers are linear) resulting in first passage time statistics that depend on the position of the bottleneck or oasis. Our results highlight the sensitivity of stochastic measures to cell division fate and quantify the limitations of using certain approximations and assumptions (such as fixed-population and mean-field assumptions) in evaluating fixation times.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Iterative graph cuts for image segmentation with a nonlinear statistical shape prior

Shape-based regularization has proven to be a useful method for delineating objects within noisy images where one has prior knowledge of the shape of the targeted object. When a collection of possible shapes is available, the specification of a shape prior using kernel density estimation is a natural technique. Unfortunately, energy functionals arising from kernel density estimation are of a form that makes them impossible to directly minimize using efficient optimization algorithms such as graph cuts. Our main contribution is to show how one may recast the energy functional into a form that is minimizable iteratively and efficiently using graph cuts.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

First passage times in homogeneous nucleation and self-assembly

Motivated by nucleation and molecular aggregation in physical, chemical and biological settings, we present a thorough analysis of the general problem of stochastic self-assembly of a fixed number of identical particles in a finite volume. We derive the Backward Kolmogorov equation (BKE) for the cluster probability distribution. From the BKE we study the distribution of times it takes for a single maximal cluster to be completed, starting from any initial particle configuration. In the limits of slow and fast self-assembly, we develop analytical approaches to calculate the mean cluster formation time and to estimate the first assembly time distribution. We find, both analytically and numerically, that faster detachment can lead to a shorter mean time to first completion of a maximum-sized cluster. This unexpected effect arises from a redistribution of trajectory weights such that upon increasing the detachment rate, paths that take a shorter time to complete a cluster become more likely.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Reconstruction of potential energy profiles from multiple rupture time distributions

We explore the mathematical and numerical aspects of reconstructing a potential energy profile of a molecular bond from its rupture time distribution. While reliable reconstruction of gross attributes, such as the height and the width of an energy barrier, can be easily extracted from a single first passage time (FPT) distribution, the reconstruction of finer structure is ill-conditioned. More careful analysis shows the existence of optimal bond potential amplitudes (represented by an effective Peclet number) and initial bond configurations that yield the most efficient numerical reconstruction of simple potentials. Furthermore, we show that reconstruction of more complex potentials containing multiple minima can be achieved by simultaneously using two or more measured FPT distributions, obtained under different physical conditions. For example, by changing the effective potential energy surface by known amounts, additional measured FPT distributions improve the reconstruction. We demonstrate the possibility of reconstructing potentials with multiple minima, motivate heuristic rules-of-thumb for optimizing the reconstruction, and discuss further applications and extensions.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Accelerated search kinetics mediated by redox reactions of DNA repair enzymes

A Charge Transport (CT) mechanism has been proposed in several papers (e.g., Yavin, et al. PNAS, v102 3546 (2005)) to explain the localization of Base Excision Repair (BER) enzymes to lesions on DNA. The CT mechanism relies on redox reactions of iron-sulfur cofactors that modify the enzyme's binding affinity. These redox reactions are mediated by the DNA strand and involve the exchange of electrons between BER enzymes along DNA. We propose a mathematical model that incorporates enzyme binding/unbinding, electron transport, and enzyme diffusion along DNA. Analysis of our model within a range of parameter values suggests that the redox reactions can increase desorption of BER enzymes not already bound to their targets, allowing the enzymes to be recycled, thus accelerating the overall search process. This acceleration mechanism is most effective when enzyme copy numbers and enzyme diffusivity along the DNA are small. Under such conditions, we find that CT BER enzymes find their targets more quickly than simple "passive" enzymes that simply attach to the DNA without desorbing.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Enhancement of charged macromolecule capture by nanopores in a salt gradient

Nanopores spanning synthetic membranes have been used as key components in proof-of-principle nanofluidic applications, particularly those involving manipulation of biomolecules or sequencing of DNA. The only practical way of manipulating charged macromolecules near nanopores is through a voltage difference applied across the nanopore-spanning membrane. However, recent experiments have shown that salt concentration gradients applied across nanopores can also dramatically enhance charged particle capture from a low concentration reservoir of charged molecules at one end of the nanopore. This puzzling effect has hitherto eluded a physically consistent theoretical explanation. Here, we propose an electrokinetic mechanism of this enhanced capture that relies on the electrostatic potential near the pore mouth. For long pores with diameter much greater than the local screening length, we obtain accurate analytic expressions showing how salt gradients control the local conductivity which can lead to increased local electrostatic potentials and charged analyte capture rates. We also find that the attractive electrostatic potential may be balanced by an outward, repulsive electroosmotic flow (EOF) that can in certain cases conspire with the salt gradient to further enhance the analyte capture rate.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Fluctuation Theorems for Entropy Production and Heat Dissipation in Periodically Driven Markov Chains

Asymptotic fluctuation theorems are statements of a Gallavotti-Cohen symmetry in the rate function of either the time-averaged entropy production or heat dissipation of a process. Such theorems have been proved for various general classes of continuous-time deterministic and stochastic processes, but always under the assumption that the forces driving the system are time independent, and often relying on the existence of a limiting ergodic distribution. In this paper we extend the asymptotic fluctuation theorem for the first time to inhomogeneous continuous-time processes without a stationary distribution, considering specifically a finite state Markov chain driven by periodic transition rates. We find that for both entropy production and heat dissipation, the usual Gallavotti-Cohen symmetry of the rate function is generalized to an analogous relation between the rate functions of the original process and its corresponding backward process, in which the trajectory and the driving protocol have been time-reversed. The effect is that spontaneous positive fluctuations in the long time average of each quantity in the forward process are exponentially more likely than spontaneous negative fluctuations in the backward process, and vice-versa, revealing that the distributions of fluctuations in universes in which time moves forward and backward are related. As an additional result, the asymptotic time-averaged entropy production is obtained as the integral of a periodic entropy production rate that generalizes the constant rate pertaining to homogeneous dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

The stochastic entry of enveloped viruses: Fusion vs. endocytosis

Viral infection requires the binding of receptors on the target cell membrane to glycoproteins, or ``spikes,'' on the viral membrane. The initial entry is usually classified as fusogenic or endocytotic. However, binding of viral spikes to cell surface receptors not only initiates the viral adhesion and the wrapping process necessary for internalization, but can simultaneously initiate direct fusion with the cell membrane. Both fusion and internalization have been observed to be viable pathways for many viruses. We develop a stochastic model for viral entry that incorporates a competition between receptor mediated fusion and endocytosis. The relative probabilities of fusion and endocytosis of a virus particle initially nonspecifically adsorbed on the host cell membrane are computed as functions of receptor concentration, binding strength, and number of spikes. We find different parameter regimes where the entry pathway probabilities can be analytically expressed. Experimental tests of our mechanistic hypotheses are proposed and discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Interplay of chemotaxis and chemokinesis mechanisms in bacterial dynamics

Motivated by observations of the dynamics of {\it Myxococcus xanthus}, we present a self-interacting random walk model that describes the competition between chemokinesis and chemotaxis. Cells are constrained to move in one dimension, but release a chemical chemoattractant at a steady state. The bacteria senses the chemical that it produces. The probability of direction reversals is modeled as a function of both the absolute level of chemoattractant sensed directly under each cell as well as the gradient sensed across the length of the cell. If the chemical does not degrade or diffuse rapidly, the one dimensional trajectory depends on the entire past history of the trajectory. We derive the corresponding Fokker-Planck equations, use an iterative mean field approach that we solve numerically for short times, and perform extensive Monte-Carlo simulations of the model. Cell positional distributions and the associated moments are computed in this feedback system. Average drift and mean squared displacements are found. Crossover behavior among different diffusion regimes are found.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

A length-dynamic Tonks gas theory of histone isotherms

We find exact solutions to a new one-dimensional (1D) interacting particle theory and apply the results to the adsorption and wrapping of polymers (such as DNA) around protein particles (such as histones). Each adsorbed protein is represented by a Tonks gas particle. The length of each particle is a degree of freedom that represents the degree of DNA wrapping around each histone. Thermodynamic quantities are computed as functions of wrapping energy, adsorbed histone density, and bulk histone concentration (or chemical potential); their experimental signatures are also discussed. Histone density is found to undergo a two-stage adsorption process as a function of chemical potential, while the mean coverage by high affinity proteins exhibits a maximum as a function of the chemical potential. However, {\it fluctuations} in the coverage are concurrently maximal. Histone-histone correlation functions are also computed and exhibit rich two length scale behavior.

0 citations0 reviewsNo code linkedOpen access