Engineering Papers⌕ Search

Engineering topics

Lythe, Grant

Publications and source records attributed to Lythe, Grant.

Stochastic journeys of cell progenies through compartments and the role of self-renewal, symmetric and asymmetric division

Abstract Division and differentiation events by which cell populations with specific functions are generated often take place as part of a developmental programme, which can be represented by a sequence of compartments. A compartment is the set of cells with common characteristics; sharing, for instance, a spatial location or a phenotype. Differentiation events are transitions from one compartment to the next. Cells may also die or divide. We consider three different types of division events: (i) where both daughter cells inherit the mother’s phenotype (self-renewal), (ii) where only one of the daughters changes phenotype (asymmetric division), and (iii) where both daughters change phenotype (symmetric division). The self-renewal probability in each compartment determines whether the progeny of a single cell, moving through the sequence of compartments, is finite or grows without bound. We analyse the progeny stochastic dynamics with probability generating functions. In the case of self-renewal, by following one of the daughters after any division event, we may construct lifelines containing only one cell at any time. We analyse the number of divisions along such lines, and the compartment where lines terminate with a death event. Analysis and numerical simulations are applied to a five-compartment model of the gradual differentiation of hematopoietic stem cells and to a model of thymocyte development: from pre-double positive to single positive (SP) cells with a bifurcation to either SP4 or SP8 in the last compartment of the sequence.

Science & Technology - Other Topics↗

The reproduction number and its probability distribution for stochastic viral dynamics

We consider stochastic models of individual infected cells. The reproduction number, R, is understood as a random variable representing the number of new cells infected by one initial infected cell in an otherwise susceptible (target cell) population. Variability in R results partly from heterogeneity in the viral burst size (the number of viral progeny generated from an infected cell during its lifetime), which depends on the distribution of cellular lifetimes and on the mechanism of virion release. We analyse viral dynamics models with an eclipse phase: the period of time after a cell is infected but before it is capable of releasing virions. The duration of the eclipse, or the subsequent infectious, phase is non-exponential, but composed of stages. We derive the probability distribution of the reproduction number for these viral dynamics models, and show it is a negative binomial distribution in the case of constant viral release from infectious cells, and under the assumption of an excess of target cells. In a deterministic model, the ultimate in-host establishment or extinction of the viral infection depends entirely on whether the mean reproduction number is greater than, or less than, one, respectively. Here, the probability of extinction is determined by the probability distribution of R, not simply its mean value. In particular, we show that in some cases the probability of infection is not an increasing function of the mean reproduction number.

59 BASIC BIOLOGICAL SCIENCES↗

Multi-variate model of T cell clonotype competition and homeostasis

Abstract Diversity of the naive T cell repertoire is maintained by competition for stimuli provided by self-peptides bound to major histocompatibility complexes (self-pMHCs). We extend an existing bi-variate competition model to a multi-variate model of the dynamics of multiple T cell clonotypes which share stimuli. In order to understand the late-time behaviour of the system, we analyse: (i) the dynamics until the extinction of the first clonotype, (ii) the time to the first extinction event, (iii) the probability of extinction of each clonotype, and (iv) the size of the surviving clonotypes when the first extinction event takes place. We also find the probability distribution of the number of cell divisions per clonotype before its extinction. The mean size of a new clonotype at quasi-steady state is an increasing function of the stimulus available to it, and a decreasing function of the fraction of stimuli it shares with other clonotypes. Thus, the probability of, and time to, extinction of a new clonotype entering the pool of T cell clonotypes is determined by the extent of competition for stimuli it experiences and by its initial number of cells.

59 BASIC BIOLOGICAL SCIENCES↗

Why are cell populations maintained via multiple compartments?

We consider the maintenance of ‘product’ cell populations from ‘progenitor’ cells via a sequence of one or more cell types, or compartments, where each cell’s fate is chosen stochastically. If there is only one compartment then large amplification, that is, a large ratio of product cells to progenitors comes with disadvantages. The product cell population is dominated by large families (cells descended from the same progenitor) and many generations separate, on average, product cells from progenitors. These disadvantages are avoided using suitably constructed sequences of compartments: the amplification factor of a sequence is the product of the amplification factors of each compartment, while the average number of generations is a sum over contributions from each compartment. Passing through multiple compartments is, in fact, an efficient way to maintain a product cell population from a small flux of progenitors, avoiding excessive clonality and minimizing the number of rounds of division en route. We use division, exit and death rates, estimated from measurements of single-positive thymocytes, to choose illustrative parameter values in the single-compartment case. We also consider a five-compartment model of thymocyte differentiation, from double-negative precursors to single-positive product cells.

59 BASIC BIOLOGICAL SCIENCES↗