The role of dimerisation and nuclear transport in the Hes1 gene regulatory network

Hes1 is a member of the family of basic helix-loop-helix transcription factors and the Hes1 gene regulatory network (GRN) may be described as the canonical example of transcriptional control in eukaryotic cells, since it involves only the Hes1 protein and its own mRNA. Recently, the Hes1 protein has been established as an excellent target for an anti-cancer drug treatment, with the design of a small molecule Hes1 dimerisation inhibitor representing a promising if challenging approach to therapy. In this paper, we extend a previous spatial stochastic model of the Hes1 GRN to include nuclear transport and dimerisation of Hes1 monomers. Initially, we assume that dimerisation occurs only in the cytoplasm, with only dimers being imported into the nucleus. Stochastic simulations of this novel model using the URDME software show that oscillatory dynamics in agreement with experimental studies are retained. Furthermore, we find that our model is robust to changes in the nuclear transport and dimerisation parameters. However, since the precise dynamics of the nuclear import of Hes1 and the localisation of the dimerisation reaction are not known, we consider a second modelling scenario in which we allow for both Hes1 monomers and dimers to be imported into the nucleus, and we allow dimerisation of Hes1 to occur everywhere in the cell. Once again, computational solutions of this second model produce oscillatory dynamics in agreement with experimental studies. We also explore sensitivity of the numerical solutions to nuclear transport and dimerisation parameters. Finally, we compare and contrast the two different modelling scenarios using numerical experiments that simulate dimer disruption, and suggest a biological experiment that could distinguish which model more faithfully captures the Hes1 GRN.PostprintPeer reviewe

Data-driven spatio-temporal modelling of glioblastoma

Mathematical oncology provides unique and invaluable insights into tumour growth on both the microscopic and macroscopic levels. This review presents state-of-the-art modelling techniques and focuses on their role in understanding glioblastoma, a malignant form of brain cancer. For each approach, we summarise the scope, drawbacks, and assets. We highlight the potential clinical applications of each modelling technique and discuss the connections between the mathematical models and the molecular and imaging data used to inform them. By doing so, we aim to prime cancer researchers with current and emerging computational tools for understanding tumour progression. Finally, by providing an in-depth picture of the different modelling techniques, we also aim to assist researchers who seek to build and develop their own models and the associated inference frameworks.Comment: 30 pages, 3 figures, 3 table

Hopf bifurcation in a gene regulatory network model:Molecular movement causes oscillations

Gene regulatory networks, i.e. DNA segments in a cell which interact with each other indirectly through their RNA and protein products, lie at the heart of many important intracellular signal transduction processes. In this paper we analyse a mathematical model of a canonical gene regulatory network consisting of a single negative feedback loop between a protein and its mRNA (e.g. the Hes1 transcription factor system). The model consists of two partial differential equations describing the spatio-temporal interactions between the protein and its mRNA in a 1-dimensional domain. Such intracellular negative feedback systems are known to exhibit oscillatory behaviour and this is the case for our model, shown initially via computational simulations. In order to investigate this behaviour more deeply, we next solve our system using Green's functions and then undertake a linearized stability analysis of the steady states of the model. Our results show that the diffusion coefficient of the protein/mRNA acts as a bifurcation parameter and gives rise to a Hopf bifurcation. This shows that the spatial movement of the mRNA and protein molecules alone is sufficient to cause the oscillations. This has implications for transcription factors such as p53, NF-$\kappa$B and heat shock proteins which are involved in regulating important cellular processes such as inflammation, meiosis, apoptosis and the heat shock response, and are linked to diseases such as arthritis and cancer

A New Method to Reconstruct Recombination Events at a Genomic Scale

Recombination is one of the main forces shaping genome diversity, but the information it generates is often overlooked. A recombination event creates a junction between two parental sequences that may be transmitted to the subsequent generations. Just like mutations, these junctions carry evidence of the shared past of the sequences. We present the IRiS algorithm, which detects past recombination events from extant sequences and specifies the place of each recombination and which are the recombinants sequences. We have validated and calibrated IRiS for the human genome using coalescent simulations replicating standard human demographic history and a variable recombination rate model, and we have fine-tuned IRiS parameters to simultaneously optimize for false discovery rate, sensitivity, and accuracy in placing the recombination events in the sequence. Newer recombinations overwrite traces of past ones and our results indicate more recent recombinations are detected by IRiS with greater sensitivity. IRiS analysis of the MS32 region, previously studied using sperm typing, showed good concordance with estimated recombination rates. We also applied IRiS to haplotypes for 18 X-chromosome regions in HapMap Phase 3 populations. Recombination events detected for each individual were recoded as binary allelic states and combined into recotypes. Principal component analysis and multidimensional scaling based on recotypes reproduced the relationships between the eleven HapMap Phase III populations that can be expected from known human population history, thus further validating IRiS. We believe that our new method will contribute to the study of the distribution of recombination events across the genomes and, for the first time, it will allow the use of recombination as genetic marker to study human genetic variation

Defining the scope for altering rice leaf anatomy to improve photosynthesis: a modelling approach.

Leaf structure plays an important role in photosynthesis. However, the causal relationship and the quantitative importance of any single structural parameter to the overall photosynthetic performance of a leaf remains open to debate. In this paper, we report on a mechanistic model, eLeaf, which successfully captures rice leaf photosynthetic performance under varying environmental conditions of light and CO2. We developed a 3D reaction-diffusion model for leaf photosynthesis parameterised using a range of imaging data and biochemical measurements from plants grown under ambient and elevated CO2 and then interrogated the model to quantify the importance of these elements. The model successfully captured leaf-level photosynthetic performance in rice. Photosynthetic metabolism underpinned the majority of the increased carbon assimilation rate observed under elevated CO2 levels, with a range of structural elements making positive and negative contributions. Mesophyll porosity could be varied without any major outcome on photosynthetic performance, providing a theoretical underpinning for experimental data. eLeaf allows quantitative analysis of the influence of morphological and biochemical properties on leaf photosynthesis. The analysis highlights a degree of leaf structural plasticity with respect to photosynthesis of significance in the context of attempts to improve crop photosynthesis

Identification of genetic variants associated with Huntington's disease progression: a genome-wide association study

Background Huntington's disease is caused by a CAG repeat expansion in the huntingtin gene, HTT. Age at onset has been used as a quantitative phenotype in genetic analysis looking for Huntington's disease modifiers, but is hard to define and not always available. Therefore, we aimed to generate a novel measure of disease progression and to identify genetic markers associated with this progression measure. Methods We generated a progression score on the basis of principal component analysis of prospectively acquired longitudinal changes in motor, cognitive, and imaging measures in the 218 indivduals in the TRACK-HD cohort of Huntington's disease gene mutation carriers (data collected 2008–11). We generated a parallel progression score using data from 1773 previously genotyped participants from the European Huntington's Disease Network REGISTRY study of Huntington's disease mutation carriers (data collected 2003–13). We did a genome-wide association analyses in terms of progression for 216 TRACK-HD participants and 1773 REGISTRY participants, then a meta-analysis of these results was undertaken. Findings Longitudinal motor, cognitive, and imaging scores were correlated with each other in TRACK-HD participants, justifying use of a single, cross-domain measure of disease progression in both studies. The TRACK-HD and REGISTRY progression measures were correlated with each other (r=0·674), and with age at onset (TRACK-HD, r=0·315; REGISTRY, r=0·234). The meta-analysis of progression in TRACK-HD and REGISTRY gave a genome-wide significant signal (p=1·12 × 10−10) on chromosome 5 spanning three genes: MSH3, DHFR, and MTRNR2L2. The genes in this locus were associated with progression in TRACK-HD (MSH3 p=2·94 × 10−8 DHFR p=8·37 × 10−7 MTRNR2L2 p=2·15 × 10−9) and to a lesser extent in REGISTRY (MSH3 p=9·36 × 10−4 DHFR p=8·45 × 10−4 MTRNR2L2 p=1·20 × 10−3). The lead single nucleotide polymorphism (SNP) in TRACK-HD (rs557874766) was genome-wide significant in the meta-analysis (p=1·58 × 10−8), and encodes an aminoacid change (Pro67Ala) in MSH3. In TRACK-HD, each copy of the minor allele at this SNP was associated with a 0·4 units per year (95% CI 0·16–0·66) reduction in the rate of change of the Unified Huntington's Disease Rating Scale (UHDRS) Total Motor Score, and a reduction of 0·12 units per year (95% CI 0·06–0·18) in the rate of change of UHDRS Total Functional Capacity score. These associations remained significant after adjusting for age of onset. Interpretation The multidomain progression measure in TRACK-HD was associated with a functional variant that was genome-wide significant in our meta-analysis. The association in only 216 participants implies that the progression measure is a sensitive reflection of disease burden, that the effect size at this locus is large, or both. Knockout of Msh3 reduces somatic expansion in Huntington's disease mouse models, suggesting this mechanism as an area for future therapeutic investigation