Chiral Polymerization in Open Systems From Chiral-Selective Reaction Rates

We investigate the possibility that prebiotic homochirality can be achieved exclusively through chiral-selective reaction rate parameters without any other explicit mechanism for chiral bias. Specifically, we examine an open network of polymerization reactions, where the reaction rates can have chiral-selective values. The reactions are neither autocatalytic nor do they contain explicit enantiomeric cross-inhibition terms. We are thus investigating how rare a set of chiral-selective reaction rates needs to be in order to generate a reasonable amount of chiral bias. We quantify our results adopting a statistical approach: varying both the mean value and the rms dispersion of the relevant reaction rates, we show that moderate to high levels of chiral excess can be achieved with fairly small chiral bias, below 10%. Considering the various unknowns related to prebiotic chemical networks in early Earth and the dependence of reaction rates to environmental properties such as temperature and pressure variations, we argue that homochirality could have been achieved from moderate amounts of chiral selectivity in the reaction rates.Comment: 15 pages, 6 figures, accepted for publication in Origins of Life and Evolution of Biosphere

A mathematical model of the human metabolic system and metabolic flexibility

In healthy subjects some tissues in the human body display metabolic flexibility, by this we mean the ability for the tissue to switch its fuel source between predominantly carbohydrates in the post prandial state and predominantly fats in the fasted state. Many of the pathways involved with human metabolism are controlled by insulin, and insulin- resistant states such as obesity and type-2 diabetes are characterised by a loss or impairment of metabolic flexibility. In this paper we derive a system of 12 first-order coupled differential equations that describe the transport between and storage in different tissues of the human body. We find steady state solutions to these equations and use these results to nondimensionalise the model. We then solve the model numerically to simulate a healthy balanced meal and a high fat meal and we discuss and compare these results. Our numerical results show good agreement with experimental data where we have data available to us and the results show behaviour that agrees with intuition where we currently have no data with which to compare

Modelling the regulation of telomere length: the effects of telomerase and G-quadruplex stabilising drugs

Telomeres are guanine-rich sequences at the end of chromosomes which shorten during each replication event and trigger cell cycle arrest and/or controlled death (apoptosis) when reaching a threshold length. The enzyme telomerase replenishes the ends of telomeres and thus prolongs the life span of cells, but also causes cellular immortalisation in human cancer. G-quadruplex (G4) stabilising drugs are a potential anticancer treatment which work by changing the molecular structure of telomeres to inhibit the activity of telomerase. We investigate the dynamics of telomere length in different conformational states, namely t-loops, G-quadruplex structures and those being elongated by telomerase. By formulating deterministic differential equation models we study the effects of various levels of both telomerase and concentrations of a G4-stabilising drug on the distribution of telomere lengths, and analyse how these effects evolve over large numbers of cell generations. As well as calculating numerical solutions, we use quasicontinuum methods to approximate the behaviour of the system over time, and predict the shape of the telomere length distribution. We find those telomerase and G4-concentrations where telomere length maintenance is successfully regulated. Excessively high levels of telomerase lead to continuous telomere lengthening, whereas large concentrations of the drug lead to progressive telomere erosion. Furthermore, our models predict a positively skewed distribution of telomere lengths, that is, telomeres accumulate over lengths shorter than the mean telomere length at equilibrium. Our model results for telomere length distributions of telomerase-positive cells in drug-free assays are in good agreement with the limited amount of experimental data available

Mathematical modelling of telomere length dynamics

Telomeres are repetitive DNA sequences located at the ends of chromosomes. During cell division, an incomplete copy of each chromosome’s DNA is made, causing telomeres to shorten on successive generations. When a threshold length is reached replication ceases and the cell becomes ‘senescent’. In this paper, we consider populations of telomeres and, from discrete models, we derive partial differential equations which describe how the distribution of telomere lengths evolves over many generations. We initially consider a population of cells each containing just a single telomere. We use continuum models to compare the effects of various mechanisms of telomere shortening and rates of cell division during normal ageing. For example, the rate (or probability) of cell replication may be fixed or it may decrease as the telomeres shorten. Furthermore, the length of telomere lost on each replication may be constant, or may decrease as the telomeres shorten. Where possible, explicit solutions for the evolution of the distribution of telomere lengths are presented. In other cases, expressions for the mean of the distribution are derived. We extend the models to describe cell populations in which each cell contains a distinct subpopulation of chromosomes. As for the simpler models, constant telomere shortening leads to a linear reduction in telomere length over time, whereas length-dependent shortening results in initially rapid telomere length reduction, slowing at later times. Our analysis also reveals that constant telomere loss leads to a Gaussian (normal) distribution of telomere lengths, whereas length-dependent loss leads to a log-normal distribution. We show that stochastic models, which include a replication probability, also lead to telomere length distributions which are skewed

DNA charge neutralization by linear polymers: Irreversible binding

We develop a deterministic mathematical model to describe the way in which polymers bind to DNA by considering the dynamics of the gap distribution that forms when polymers bind to a DNA plasmid. In so doing, we generalize existing theory to account for overlaps and binding cooperativity whereby the polymer binding rate depends on the size of the overlap. The proposed mean-field models are then solved using a combination of numerical and asymptotic methods. We find that overlaps lead to higher coverage and hence higher charge neutralizations, results which are more in line with recent experimental observations. Our work has applications to gene therapy where polymers are used to neutralize the negative charges of the DNA phosphate backbone, allowing condensation prior to delivery into the nucleus of an abnormal cell. © 2006 The American Physical Society