Publications

You can also access our individual websites (via the Members page) for further information about our research and lists of our publications.

Search or filter publications

Filter by type:

Filter by publication type

Filter by year:

to

Results

  • Showing results for:
  • Reset all filters

Search results

  • Journal article
    Aryaman J, Johnston IG, Jones NS, 2017,

    Mitochondrial DNA Density Homeostasis Accounts for a Threshold Effect in a Cybrid Model of a Human Mitochondrial Disease

    , Biochemical Journal, Vol: 474, Pages: 4019-4034, ISSN: 1470-8728

    Mitochondrial dysfunction is involved in a wide array of devastating diseases, but the heterogeneity and complexity of the symptoms of these diseases challenges theoretical understanding of their causation. With the explosion of omics data, we have the unprecedented opportunity to gain deep understanding of the biochemical mechanisms of mitochondrial dysfunction. This goal raises the outstanding need to make these complex datasets interpretable. Quantitative modelling allows us to translate such datasets into intuition and suggest rational biomedical treatments. Taking an interdisciplinary approach, we use a recently published large-scale dataset and develop a descriptive and predictive mathematical model of progressive increase in mutant load of the MELAS 3243A>G mtDNA mutation. The experimentally observed behaviour is surprisingly rich, but we find that our simple, biophysically motivated model intuitively accounts for this heterogeneity and yields a wealth of biological predictions. Our findings suggest that cells attempt to maintain wild-type mtDNA density through cell volume reduction, and thus power demand reduction, until a minimum cell volume is reached. Thereafter, cells toggle from demand reduction to supply increase, up-regulating energy production pathways. Our analysis provides further evidence for the physiological significance of mtDNA density and emphasizes the need for performing single-cell volume measurements jointly with mtDNA quantification. We propose novel experiments to verify the hypotheses made here to further develop our understanding of the threshold effect and connect with rational choices for mtDNA disease therapies.

  • Journal article
    Fulcher B, Jones NS, 2017,

    hctsa: A computational framework for automated timeseriesphenotyping using massive feature extraction

    , Cell Systems, Vol: 5, Pages: 527-531.e3, ISSN: 2405-4712

    Phenotype measurements frequently take the form of time series, but we currently lack a systematic method for relating these complex data streams to scientifically meaningful outcomes, such as relating the movement dynamics of organisms to their genotype or measurements of brain dynamics of a patient to their disease diagnosis. Previous work addressed this problem by comparing implementations of thousands of diverse scientific time-series analysis methods in an approach termed highly comparative time-series analysis. Here, we introduce hctsa, a software tool for applying this methodological approach to data. hctsa includes an architecture for computing over 7,700 time-series features and a suite of analysis and visualization algorithms to automatically select useful and interpretable time-series features for a given application. Using exemplar applications to high-throughput phenotyping experiments, we show how hctsa allows researchers to leverage decades of time-series research to quantify and understand informative structure in time-series data.

  • Journal article
    Liu Z, Barahona M, 2017,

    Geometric multiscale community detection: Markov stability and vector partitioning

    , Journal of Complex Networks, Vol: 6, Pages: 157-172, ISSN: 2051-1329

    Multiscale community detection can be viewed from a dynamical perspective within the Markov stability framework, which uses the diffusion of a Markov process on the graph to uncover intrinsic network substructures across all scales. Here we reformulate multiscale community detection as a max-sum length vector partitioning problem with respect to the set of time-dependent node vectors expressed in terms of eigenvectors of the transition matrix. This formulation provides a geometric interpretation of Markov stability in terms of a time-dependent spectral embedding, where the Markov time acts as an inhomogeneous geometric resolution factor that zooms the components of the node vectors at different rates. Our geometric formulation encompasses both modularity and the multi-resolution Potts model, which are shown to correspond to vector partitioning in a pseudo-Euclidean space, and is also linked to spectral partitioning methods, where the number of eigenvectors used corresponds to the dimensionality of the underlying embedding vector space. Inspired by the Louvain optimization for community detection, we then propose an algorithm based on a graph-theoretical heuristic for the vector partitioning problem. We apply the algorithm to the spectral optimization of modularity and Markov stability community detection. The spectral embedding based on the transition matrix eigenvectors leads to improved partitions with higher information content and higher modularity than the eigen-decomposition of the modularity matrix. We illustrate the results with random network benchmarks.

  • Journal article
    Shahrezaei V, Robertson B, Thomas P, Priestman Met al., 2017,

    Mycobacteria modify their cell size control under sub-optimal carbon sources

    , Frontiers in Cell and Developmental Biology, Vol: 5, ISSN: 2296-634X

    The decision to divide is the most important one that any cell must make. Recent single cell studies suggest that most bacteria follow an “adder” model of cell size control, incorporating a fixed amount of cell wall material before dividing. Mycobacteria, including the causative agent of tuberculosis Mycobacterium tuberculosis, are known to divide asymmetrically resulting in heterogeneity in growth rate, doubling time, and other growth characteristics in daughter cells. The interplay between asymmetric cell division and adder size control has not been extensively investigated. Moreover, the impact of changes in the environment on growth rate and cell size control have not been addressed for mycobacteria. Here, we utilize time-lapse microscopy coupled with microfluidics to track live Mycobacterium smegmatis cells as they grow and divide over multiple generations, under a variety of growth conditions. We demonstrate that, under optimal conditions, M. smegmatis cells robustly follow the adder principle, with constant added length per generation independent of birth size, growth rate, and inherited pole age. However, the nature of the carbon source induces deviations from the adder model in a manner that is dependent on pole age. Understanding how mycobacteria maintain cell size homoeostasis may provide crucial targets for the development of drugs for the treatment of tuberculosis, which remains a leading cause of global mortality.

  • Journal article
    Branch T, Barahona M, Dodson C, Ying Let al., 2017,

    Kinetic analysis reveals the identity of Aβ-metal complex responsible for the initial aggregation of Aβ in the synapse

    , ACS Chemical Neuroscience, Vol: 8, Pages: 1970-1979, ISSN: 1948-7193

    The mechanism of Aβ aggregation in the absence of metal ions is well established, yet the role that Zn2+ and Cu2+, the two most studied metal ions, released during neurotransmission, paly in promoting Aβ aggregation in the vicinity of neuronal synapses remains elusive. Here we report the kinetics of Zn2+ binding to Aβ and Zn2+/Cu2+ binding to Aβ-Cu to form ternary complexes under near physiological conditions (nM Aβ, μM metal ions). We find that these reactions are several orders of magnitude slower than Cu2+ binding to Aβ. Coupled reaction-diffusion simulations of the interactions of synaptically released metal ions with Aβ show that up to a third of Aβ is Cu2+-bound under repetitive metal ion release, while any other Aβ-metal complexes (including Aβ-Zn) are insignificant. We therefore conclude that Zn2+ is unlikely to play an important role in the very early stages (i.e., dimer formation) of Aβ aggregation, contrary to a widely held view in the subject. We propose that targeting the specific interactions between Cu2+ and Aβ may be a viable option in drug development efforts for early stages of AD.

  • Journal article
    Aryaman J, hoitzing H, burgstaller J, johnston I, Jones NSet al., 2017,

    Mitochondrial heterogeneity, metabolic scaling and cell death

    , Bioessays, Vol: 39, ISSN: 1521-1878

    Heterogeneity in mitochondrial content has been previously suggested as a major contributor to cellular noise, with multiple studies indicating its direct involvement in biomedically important cellular phenomena. A recently published dataset explored the connection between mitochondrial functionality and cell physiology, where a non-linearity between mitochondrial functionality and cell size was found. Using mathematical models, we suggest that a combination of metabolic scaling and a simple model of cell death may account for these observations. However, our findings also suggest the existence of alternative competing hypotheses, such as a non-linearity between cell death and cell size. While we find that the proposed non-linear coupling between mitochondrial functionality and cell size provides a compelling alternative to previous attempts to link mitochondrial heterogeneity and cell physiology, we emphasise the need to account for alternative causal variables, including cell cycle, size, mitochondrial density and death, in future studies of mitochondrial physiology.

  • Journal article
    Johnson S, Jones NS, 2017,

    Looplessness in networks is linked to trophic coherence

    , Proceedings of the National Academy of Sciences of USA, Vol: 114, Pages: 5618-5623, ISSN: 0027-8424

    Many natural, complex systems are remarkably stable thanks to anabsence of feedback acting on their elements. When described as net-works, these exhibit few or no cycles, and associated matrices have smallleading eigenvalues. It has been suggested that this architecture can con-fer advantages to the system as a whole, such as ‘qualitative stability’,but this observation does not in itself explain how a loopless structuremight arise. We show here that the number of feedback loops in a net-work, as well as the eigenvalues of associated matrices, are determined bya structural property called trophic coherence, a measure of how neatlynodes fall into distinct levels. Our theory correctly classifies a variety ofnetworks – including those derived from genes, metabolites, species, neu-rons, words, computers and trading nations – into two distinct regimesof high and low feedback, and provides a null model to gauge the signifi-cance of related magnitudes. Since trophic coherence suppresses feedback,whereas an absence of feedback alone does not lead to coherence, our worksuggests that the reasons for ‘looplessness’ in nature should be sought incoherence-inducing mechanisms.

  • Book chapter
    Andreychenko A, Bortolussi L, Grima R, Thomas P, Wolf Vet al., 2017,

    Distribution approximations for the chemical master equation: comparisonof the method of moments and the system size expansion

    , Modeling Cellular Systems, Editors: Graw, Matthaus, Pahle, Publisher: Springer, Pages: 39-39, ISBN: 978-3-319-45833-5

    The stochastic nature of chemical reactions involving randomly fluctuatingpopulation sizes has lead to a growing research interest in discrete-statestochastic models and their analysis. A widely-used approach is the descriptionof the temporal evolution of the system in terms of a chemical master equation(CME). In this paper we study two approaches for approximating the underlyingprobability distributions of the CME. The first approach is based on anintegration of the statistical moments and the reconstruction of thedistribution based on the maximum entropy principle. The second approach relieson an analytical approximation of the probability distribution of the CME usingthe system size expansion, considering higher-order terms than the linear noiseapproximation. We consider gene expression networks with unimodal andmultimodal protein distributions to compare the accuracy of the two approaches.We find that both methods provide accurate approximations to the distributionsof the CME while having different benefits and limitations in applications.

  • Journal article
    Kiselev V, Kirschner K, Schaub MT, Andrews T, Yiu A, Chandra T, Natarajan KN, Reik W, Barahona M, Green AR, Hemberg Met al., 2017,

    SC3: consensus clustering of single-cell RNA-seq data

    , Nature Methods, Vol: 14, Pages: 483-486, ISSN: 1548-7105

    Single-cell RNA-seq enables the quantitative characterization of cell types based on global transcriptome profiles. We present single-cell consensus clustering (SC3), a user-friendly tool for unsupervised clustering, which achieves high accuracy and robustness by combining multiple clustering solutions through a consensus approach (http://bioconductor.org/packages/SC3). We demonstrate that SC3 is capable of identifying subclones from the transcriptomes of neoplastic cells collected from patients.

  • Journal article
    Gosztolai A, Schumacher J, Behrends V, Bundy JG, Heydenreich F, Bennett MH, Buck M, Barahona Met al., 2017,

    GlnK Facilitates the Dynamic Regulation of Bacterial Nitrogen Assimilation.

    , Biophysical journal, Vol: 112, Pages: 2219-2230, ISSN: 0006-3495

    Ammonium assimilation in Escherichia coli is regulated by two paralogous proteins (GlnB and GlnK), which orchestrate interactions with regulators of gene expression, transport proteins, and metabolic pathways. Yet how they conjointly modulate the activity of glutamine synthetase, the key enzyme for nitrogen assimilation, is poorly understood. We combine experiments and theory to study the dynamic roles of GlnB and GlnK during nitrogen starvation and upshift. We measure time-resolved in vivo concentrations of metabolites, total and posttranslationally modified proteins, and develop a concise biochemical model of GlnB and GlnK that incorporates competition for active and allosteric sites, as well as functional sequestration of GlnK. The model predicts the responses of glutamine synthetase, GlnB, and GlnK under time-varying external ammonium level in the wild-type and two genetic knock-outs. Our results show that GlnK is tightly regulated under nitrogen-rich conditions, yet it is expressed during ammonium run-out and starvation. This suggests a role for GlnK as a buffer of nitrogen shock after starvation, and provides a further functional link between nitrogen and carbon metabolisms.

This data is extracted from the Web of Science and reproduced under a licence from Thomson Reuters. You may not copy or re-distribute this data in whole or in part without the written consent of the Science business of Thomson Reuters.

Request URL: http://www.imperial.ac.uk:80/respub/WEB-INF/jsp/search-t4-html.jsp Request URI: /respub/WEB-INF/jsp/search-t4-html.jsp Query String: id=219&limit=10&page=10&respub-action=search.html Current Millis: 1714046110624 Current Time: Thu Apr 25 12:55:10 BST 2024