Flux-dependent graphs for metabolic networks

Cells adapt their metabolic fluxes in response to changes in the environment. We present a framework for the systematic construction of flux-based graphs derived from organism-wide metabolic networks. Our graphs encode the directionality of metabolic flows via edges that represent the flow of metabolites from source to target reactions. The methodology can be applied in the absence of a specific biological context by modelling fluxes probabilistically, or can be tailored to different environmental conditions by incorporating flux distributions computed through constraint-based approaches such as Flux Balance Analysis. We illustrate our approach on the central carbon metabolism of Escherichia coli and on a metabolic model of human hepatocytes. The flux-dependent graphs under various environmental conditions and genetic perturbations exhibit systemic changes in their topological and community structure, which capture the re-routing of metabolic flows and the varying importance of specific reactions and pathways. By integrating constraint-based models and tools from network science, our framework allows the study of context-specific metabolic responses at a system level beyond standard pathway descriptions.

Cellular metabolism is the result of a highly enmeshed set of biochemical reactions that is naturally amenable to graph-based analyses. Yet there are multiple ways to construct a graph representation from any given metabolic model. Here an international research team of UK and Spain scientists presents a principled approach to study metabolic models through the lens of network science. They propose a framework to construct graphs for genome-scale metabolic models that resolve various challenges, such as the incorporation of pool metabolites, the preservation of the directionality of metabolic flows, and the capability to incorporate specific flux information. The method can be integrated into pipelines based on flux balance analysis and provides a systematic framework to explore changes in network connectivity as a result of environmental shifts or genetic perturbations. The framework thus allows to interrogate context-specific metabolic responses beyond standard pathway descriptions. The authors illustrate the approach through the analysis of Escherichia coli's core metabolism in different growth conditions, as well as a rare metabolic disease affecting human hepatocytes.

Allostery and cooperativity in multimeric proteins: bond-to-bond propensities in ATCase

Aspartate carbamoyltransferase (ATCase) is a large dodecameric enzyme with six active sites that exhibits allostery: its catalytic rate is modulated by the binding of various substrates at distal points from the active sites. A recently developed method, bond-to-bond propensity analysis, has proven capable of predicting allosteric sites in a wide range of proteins using an energy-weighted atomistic graph obtained from the protein structure and given knowledge only of the location of the active site. Bond-to-bond propensity establishes if energy fluctuations at given bonds have significant effects on any other bond in the protein, by considering their propagation through the protein graph. In this work, we use bond-to-bond propensity analysis to study different aspects of ATCase activity using three different protein structures and sources of fluctuations. First, we predict key residues and bonds involved in the transition between inactive (T) and active (R) states of ATCase by analysing allosteric substrate binding as a source of energy perturbations in the protein graph. Our computational results also indicate that the effect of multiple allosteric binding is non linear: a switching effect is observed after a particular number and arrangement of substrates is bound suggesting a form of long range communication between the distantly arranged allosteric sites. Second, cooperativity is explored by considering a bisubstrate analogue as the source of energy fluctuations at the active site, also leading to the identification of highly significant residues to the T ↔ R transition that enhance cooperativity across active sites. Finally, the inactive (T) structure is shown to exhibit a strong, non linear communication between the allosteric sites and the interface between catalytic subunits, rather than the active site. Bond-to-bond propensity thus offers an alternative route to explain allosteric and cooperative effects in terms of detailed atomistic changes to individual bonds within the protein, rather than through phenomenological, global thermodynamic arguments.

Computational Re-design of Synthetic Genetic Oscillators for Independent Amplitude and Frequency Modulation

To perform well in biotechnology applications, synthetic genetic oscillators must be engineered to allow independent modulation of amplitude and period. This need is currently unmet. Here, we demonstrate computationally how two classic genetic oscillators, the dual-feedback oscillator and the repressilator, can be re-designed to provide independent control of amplitude and period and improve tunability-that is, a broad dynamic range of periods and amplitudes accessible through the input "dials." Our approach decouples frequency and amplitude modulation by incorporating an orthogonal "sink module" where the key molecular species are channeled for enzymatic degradation. This sink module maintains fast oscillation cycles while alleviating the translational coupling between the oscillator's transcription factors and output. We characterize the behavior of our re-designed oscillators over a broad range of physiologically reasonable parameters, explain why this facilitates broader function and control, and provide general design principles for building synthetic genetic oscillators that are more precisely controllable.

Stochastic models of gene transcription with upstream drives: exact solution and sample path characterization

Gene transcription is a highly stochastic and dynamic process. As a result, the mRNA copy number of a given gene is heterogeneous both between cells and across time. We present a framework to model gene transcription in populations of cells with time-varying (stochastic or deterministic) transcription and degradation rates. Such rates can be understood as upstream cellular drives representing the effect of different aspects of the cellular environment. We show that the full solution of the master equation contains two components: a model-specific, upstream effective drive, which encapsulates the effect of cellular drives (e.g. entrainment, periodicity or promoter randomness) and a downstream transcriptional Poissonian part, which is common to all models. Our analytical framework treats cell-to-cell and dynamic variability consistently, unifying several approaches in the literature. We apply the obtained solution to characterize different models of experimental relevance, and to explain the influence on gene transcription of synchrony, stationarity, ergodicity, as well as the effect of time scales and other dynamic characteristics of drives. We also show how the solution can be applied to the analysis of noise sources in single-cell data, and to reduce the computational cost of stochastic simulations.

Kinetic Analysis Reveals the Identity of Aβ-Metal Complex Responsible for the Initial Aggregation of Aβ in the Synapse

The mechanism of Aβ aggregation in the absence of metal ions is well established, yet the role that Zn²⁺ and Cu²⁺, 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 Zn²⁺ binding to Aβ and Zn²⁺/Cu²⁺ 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 Cu²⁺ 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 Cu²⁺-bound under repetitive metal ion release, while any other Aβ-metal complexes (including Aβ-Zn) are insignificant. We therefore conclude that Zn²⁺ 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 Cu²⁺ and Aβ may be a viable option in drug development efforts for early stages of AD.

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

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.

[Effect of sports training in angular compartment of the lower limbs in children footballers aged 11 to 12 years old]

Anterior Cruciate Ligament (ACL) injury is an important cause of days lost in athletes. Most ACL injuries are non-contact and are associated with biomechanical risk factors that increase tension in the ACL: increased knee valgus (KV) and hip flexion (HF) and decreased flexion of knee (KF). Muscle around the knee contributes to knee stability, so fatigue produced by exercise could alter knee balance, increasing LCA tension. The aim of the study is to determine the angular behavior before and after a physical load for CF, RR and RV in children born in 2002-2003. A non-randomized clinical trial was conducted. The sample consisted of 50 students from soccer schools born between 2002 and 2003. The angular behavior of CF, RR and VR was compared, before and after performing standardized training. The angular behavior was measured by performing the DJ test with data obtained by inertial sensors. After exercise, the 3 variables increased, but only HF reached significant difference. Other important finding was the difference found in KV between the dominant leg and the support limb, at both times: rest and post exercise. It was concluded that the angular behavior of CF increases significantly in both limbs post-exercise and that preventive measures should be applied for the management of valgus in the supporting limb.

Toward Precision Healthcare: Context and Mathematical Challenges

Precision medicine refers to the idea of delivering the right treatment to the right patient at the right time, usually with a focus on a data-centered approach to this task. In this perspective piece, we use the term "precision healthcare" to describe the development of precision approaches that bridge from the individual to the population, taking advantage of individual-level data, but also taking the social context into account. These problems give rise to a broad spectrum of technical, scientific, policy, ethical and social challenges, and new mathematical techniques will be required to meet them. To ensure that the science underpinning "precision" is robust, interpretable and well-suited to meet the policy, ethical and social questions that such approaches raise, the mathematical methods for data analysis should be transparent, robust, and able to adapt to errors and uncertainties. In particular, precision methodologies should capture the complexity of data, yet produce tractable descriptions at the relevant resolution while preserving intelligibility and traceability, so that they can be used by practitioners to aid decision-making. Through several case studies in this domain of precision healthcare, we argue that this vision requires the development of new mathematical frameworks, both in modeling and in data analysis and interpretation.

The 'who' and 'what' of #diabetes on Twitter

Social media are being increasingly used for health promotion, yet the landscape of users, messages and interactions in such fora is poorly understood. Studies of social media and diabetes have focused mostly on patients, or public agencies addressing it, but have not looked broadly at all of the participants or the diversity of content they contribute. We study Twitter conversations about diabetes through the systematic analysis of 2.5 million tweets collected over 8 months and the interactions between their authors. We address three questions. (1) What themes arise in these tweets? (2) Who are the most influential users? (3) Which type of users contribute to which themes? We answer these questions using a mixed-methods approach, integrating techniques from anthropology, network science and information retrieval such as thematic coding, temporal network analysis and community and topic detection. Diabetes-related tweets fall within broad thematic groups: health information, news, social interaction and commercial. At the same time, humorous messages and references to popular culture appear consistently, more than any other type of tweet. We classify authors according to their temporal 'hub' and 'authority' scores. Whereas the hub landscape is diffuse and fluid over time, top authorities are highly persistent across time and comprise bloggers, advocacy groups and NGOs related to diabetes, as well as for-profit entities without specific diabetes expertise. Top authorities fall into seven interest communities as derived from their Twitter follower network. Our findings have implications for public health professionals and policy makers who seek to use social media as an engagement tool and to inform policy design.

Graph partitions and cluster synchronization in networks of oscillators

Synchronization over networks depends strongly on the structure of the coupling between the oscillators. When the coupling presents certain regularities, the dynamics can be coarse-grained into clusters by means of External Equitable Partitions of the network graph and their associated quotient graphs. We exploit this graph-theoretical concept to study the phenomenon of cluster synchronization, in which different groups of nodes converge to distinct behaviors. We derive conditions and properties of networks in which such clustered behavior emerges and show that the ensuing dynamics is the result of the localization of the eigenvectors of the associated graph Laplacians linked to the existence of invariant subspaces. The framework is applied to both linear and non-linear models, first for the standard case of networks with positive edges, before being generalized to the case of signed networks with both positive and negative interactions. We illustrate our results with examples of both signed and unsigned graphs for consensus dynamics and for partial synchronization of oscillator networks under the master stability function as well as Kuramoto oscillators.

Linear models of activation cascades: analytical solutions and coarse-graining of delayed signal transduction

Cellular signal transduction usually involves activation cascades, the sequential activation of a series of proteins following the reception of an input signal. Here, we study the classic model of weakly activated cascades and obtain analytical solutions for a variety of inputs. We show that in the special but important case of optimal gain cascades (i.e. when the deactivation rates are identical) the downstream output of the cascade can be represented exactly as a lumped nonlinear module containing an incomplete gamma function with real parameters that depend on the rates and length of the cascade, as well as parameters of the input signal. The expressions obtained can be applied to the non-identical case when the deactivation rates are random to capture the variability in the cascade outputs. We also show that cascades can be rearranged so that blocks with similar rates can be lumped and represented through our nonlinear modules. Our results can be used both to represent cascades in computational models of differential equations and to fit data efficiently, by reducing the number of equations and parameters involved. In particular, the length of the cascade appears as a real-valued parameter and can thus be fitted in the same manner as Hill coefficients. Finally, we show how the obtained nonlinear modules can be used instead of delay differential equations to model delays in signal transduction.

Prediction of allosteric sites and mediating interactions through bond-to-bond propensities

Allostery is a fundamental mechanism of biological regulation, in which binding of a molecule at a distant location affects the active site of a protein. Allosteric sites provide targets to fine-tune protein activity, yet we lack computational methodologies to predict them. Here we present an efficient graph-theoretical framework to reveal allosteric interactions (atoms and communication pathways strongly coupled to the active site) without a priori information of their location. Using an atomistic graph with energy-weighted covalent and weak bonds, we define a bond-to-bond propensity quantifying the non-local effect of instantaneous bond fluctuations propagating through the protein. Significant interactions are then identified using quantile regression. We exemplify our method with three biologically important proteins: caspase-1, CheY, and h-Ras, correctly predicting key allosteric interactions, whose significance is additionally confirmed against a reference set of 100 proteins. The almost-linear scaling of our method renders it suitable for high-throughput searches for candidate allosteric sites.

Flow-Based Network Analysis of the Caenorhabditis elegans Connectome

We exploit flow propagation on the directed neuronal network of the nematode C. elegans to reveal dynamically relevant features of its connectome. We find flow-based groupings of neurons at different levels of granularity, which we relate to functional and anatomical constituents of its nervous system. A systematic in silico evaluation of the full set of single and double neuron ablations is used to identify deletions that induce the most severe disruptions of the multi-resolution flow structure. Such ablations are linked to functionally relevant neurons, and suggest potential candidates for further in vivo investigation. In addition, we use the directional patterns of incoming and outgoing network flows at all scales to identify flow profiles for the neurons in the connectome, without pre-imposing a priori categories. The four flow roles identified are linked to signal propagation motivated by biological input-response scenarios.

One of the goals of systems neuroscience is to elucidate the relationship between the structure of neuronal networks and the functional dynamics that they implement. An ideal model organism to study such interactions is the roundworm C. elegans, which not only has a fully mapped connectome, but has also been the object of extensive behavioural, genetic and neurophysiological experiments. Here we present an analysis of the neuronal network of C. elegans from a dynamical flow perspective. Our analysis reveals a multi-scale organisation of the signal flow in the network linked to anatomical and functional features of neurons, as well as identifying different neuronal roles in relation to signal propagation. We use our computational framework to explore biological input-response scenarios as well as exhaustive in silico ablations, which we relate to experimental findings reported in the literature.

Great cities look small

Great cities connect people; failed cities isolate people. Despite the fundamental importance of physical, face-to-face social ties in the functioning of cities, these connectivity networks are not explicitly observed in their entirety. Attempts at estimating them often rely on unrealistic over-simplifications such as the assumption of spatial homogeneity. Here we propose a mathematical model of human interactions in terms of a local strategy of maximizing the number of beneficial connections attainable under the constraint of limited individual travelling-time budgets. By incorporating census and openly available online multi-modal transport data, we are able to characterize the connectivity of geometrically and topologically complex cities. Beyond providing a candidate measure of greatness, this model allows one to quantify and assess the impact of transport developments, population growth, and other infrastructure and demographic changes on a city. Supported by validations of gross domestic product and human immunodeficiency virus infection rates across US metropolitan areas, we illustrate the effect of changes in local and city-wide connectivities by considering the economic impact of two contemporary inter- and intra-city transport developments in the UK: High Speed 2 and London Crossrail. This derivation of the model suggests that the scaling of different urban indicators with population size has an explicitly mechanistic origin.

Interest communities and flow roles in directed networks: the Twitter network of the UK riots

Directionality is a crucial ingredient in many complex networks in which information, energy or influence are transmitted. In such directed networks, analysing flows (and not only the strength of connections) is crucial to reveal important features of the network that might go undetected if the orientation of connections is ignored. We showcase here a flow-based approach for community detection through the study of the network of the most influential Twitter users during the 2011 riots in England. Firstly, we use directed Markov Stability to extract descriptions of the network at different levels of coarseness in terms of interest communities, i.e. groups of nodes within which flows of information are contained and reinforced. Such interest communities reveal user groupings according to location, profession, employer and topic. The study of flows also allows us to generate an interest distance, which affords a personalized view of the attention in the network as viewed from the vantage point of any given user. Secondly, we analyse the profiles of incoming and outgoing long-range flows with a combined approach of role-based similarity and the novel relaxed minimum spanning tree algorithm to reveal that the users in the network can be classified into five roles. These flow roles go beyond the standard leader/follower dichotomy and differ from classifications based on regular/structural equivalence. We then show that the interest communities fall into distinct informational organigrams characterized by a different mix of user roles reflecting the quality of dialogue within them. Our generic framework can be used to provide insight into how flows are generated, distributed, preserved and consumed in directed networks.

Emergence of Slow-Switching Assemblies in Structured Neuronal Networks

Unraveling the interplay between connectivity and spatio-temporal dynamics in neuronal networks is a key step to advance our understanding of neuronal information processing. Here we investigate how particular features of network connectivity underpin the propensity of neural networks to generate slow-switching assembly (SSA) dynamics, i.e., sustained epochs of increased firing within assemblies of neurons which transition slowly between different assemblies throughout the network. We show that the emergence of SSA activity is linked to spectral properties of the asymmetric synaptic weight matrix. In particular, the leading eigenvalues that dictate the slow dynamics exhibit a gap with respect to the bulk of the spectrum, and the associated Schur vectors exhibit a measure of block-localization on groups of neurons, thus resulting in coherent dynamical activity on those groups. Through simple rate models, we gain analytical understanding of the origin and importance of the spectral gap, and use these insights to develop new network topologies with alternative connectivity paradigms which also display SSA activity. Specifically, SSA dynamics involving excitatory and inhibitory neurons can be achieved by modifying the connectivity patterns between both types of neurons. We also show that SSA activity can occur at multiple timescales reflecting a hierarchy in the connectivity, and demonstrate the emergence of SSA in small-world like networks. Our work provides a step towards understanding how network structure (uncovered through advancements in neuroanatomy and connectomics) can impact on spatio-temporal neural activity and constrain the resulting dynamics.

Amplification of small molecule-inducible gene expression via tuning of intracellular receptor densities

Ligand-responsive transcription factors in prokaryotes found simple small molecule-inducible gene expression systems. These have been extensively used for regulated protein production and associated biosynthesis of fine chemicals. However, the promoter and protein engineering approaches traditionally used often pose significant restrictions to predictably and rapidly tune the expression profiles of inducible expression systems. Here, we present a new unified and rational tuning method to amplify the sensitivity and dynamic ranges of versatile small molecule-inducible expression systems. We employ a systematic variation of the concentration of intracellular receptors for transcriptional control. We show that a low density of the repressor receptor (e.g. TetR and ArsR) in the cell can significantly increase the sensitivity and dynamic range, whereas a high activator receptor (e.g. LuxR) density achieves the same outcome. The intracellular concentration of receptors can be tuned in both discrete and continuous modes by adjusting the strength of their cognate driving promoters. We exemplified this approach in several synthetic receptor-mediated sensing circuits, including a tunable cell-based arsenic sensor. The approach offers a new paradigm to predictably tune and amplify ligand-responsive gene expression with potential applications in synthetic biology and industrial biotechnology.

PDGFRα demarcates the cardiogenic clonogenic Sca1+ stem/progenitor cell in adult murine myocardium

Cardiac progenitor/stem cells in adult hearts represent an attractive therapeutic target for heart regeneration, though (inter)-relationships among reported cells remain obscure. Using single-cell qRT-PCR and clonal analyses, here we define four subpopulations of cardiac progenitor/stem cells in adult mouse myocardium all sharing stem cell antigen-1 (Sca1), based on side population (SP) phenotype, PECAM-1 (CD31) and platelet-derived growth factor receptor-α (PDGFRα) expression. SP status predicts clonogenicity and cardiogenic gene expression (Gata4/6, Hand2 and Tbx5/20), properties segregating more specifically to PDGFRα(+) cells. Clonal progeny of single Sca1(+) SP cells show cardiomyocyte, endothelial and smooth muscle lineage potential after cardiac grafting, augmenting cardiac function although durable engraftment is rare. PDGFRα(-) cells are characterized by Kdr/Flk1, Cdh5, CD31 and lack of clonogenicity. PDGFRα(+)/CD31(-) cells derive from cells formerly expressing Mesp1, Nkx2-5, Isl1, Gata5 and Wt1, distinct from PDGFRα(-)/CD31(+) cells (Gata5 low; Flk1 and Tie2 high). Thus, PDGFRα demarcates the clonogenic cardiogenic Sca1(+) stem/progenitor cell.

Hhex and Cer1 mediate the Sox17 pathway for cardiac mesoderm formation in embryonic stem cells

Cardiac muscle differentiation in vivo is guided by sequential growth factor signals, including endoderm-derived diffusible factors, impinging on cardiogenic genes in the developing mesoderm. Previously, by RNA interference in AB2.2 mouse embryonic stem cells (mESCs), we identified the endodermal transcription factor Sox17 as essential for Mesp1 induction in primitive mesoderm and subsequent cardiac muscle differentiation. However, downstream effectors of Sox17 remained to be proven functionally. In this study, we used genome-wide profiling of Sox17-dependent genes in AB2.2 cells, RNA interference, chromatin immunoprecipitation, and luciferase reporter genes to dissect this pathway. Sox17 was required not only for Hhex (a second endodermal transcription factor) but also for Cer1, a growth factor inhibitor from endoderm that, like Hhex, controls mesoderm patterning in Xenopus toward a cardiac fate. Suppressing Hhex or Cer1 blocked cardiac myogenesis, although at a later stage than induction of Mesp1/2. Hhex was required but not sufficient for Cer1 expression. Over-expression of Sox17 induced endogenous Cer1 and sequence-specific transcription of a Cer1 reporter gene. Forced expression of Cer1 was sufficient to rescue cardiac differentiation in Hhex-deficient cells. Thus, Hhex and Cer1 are indispensable components of the Sox17 pathway for cardiopoiesis in mESCs, acting at a stage downstream from Mesp1/2.

Engineering modular and tunable genetic amplifiers for scaling transcriptional signals in cascaded gene networks

Synthetic biology aims to control and reprogram signal processing pathways within living cells so as to realize repurposed, beneficial applications. Here we report the design and construction of a set of modular and gain-tunable genetic amplifiers in Escherichia coli capable of amplifying a transcriptional signal with wide tunable-gain control in cascaded gene networks. The devices are engineered using orthogonal genetic components (hrpRS, hrpV and P_hrpL) from the hrp (hypersensitive response and pathogenicity) gene regulatory network in Pseudomonas syringae. The amplifiers can linearly scale up to 21-fold the transcriptional input with a large output dynamic range, yet not introducing significant time delay or significant noise during signal amplification. The set of genetic amplifiers achieves different gains and input dynamic ranges by varying the expression levels of the underlying ligand-free activator proteins in the device. As their electronic counterparts, these engineered transcriptional amplifiers can act as fundamental building blocks in the design of biological systems by predictably and dynamically modulating transcriptional signal flows to implement advanced intra- and extra-cellular control functions.

Uncovering allosteric pathways in caspase-1 using Markov transient analysis and multiscale community detection

Atomistic graph–theoretical analysis of caspase-1 reveals details of intra-protein communication pathways.

Observability and coarse graining of consensus dynamics through the external equitable partition

Using the intrinsic relationship between the external equitable partition (EEP) and the spectral properties of the graph Laplacian, we characterize convergence and observability properties of consensus dynamics on networks. In particular, we establish the relationship between the original consensus dynamics and the associated consensus of the quotient graph under varied initial conditions, and characterize the asymptotic convergence to the synchronization manifold under nonuniform input signals. We also show that the EEP with respect to a node can reveal nodes in the graph with an increased rate of asymptotic convergence to the consensus value, as characterized by the second smallest eigenvalue of the quotient Laplacian. Finally, we show that the quotient graph preserves the observability properties of the full graph and how the inheritance by the quotient graph of particular aspects of the eigenstructure of the full Laplacian underpins the observability and convergence properties of the system.

Risk factor-dependent dynamics of atopic dermatitis: modelling multi-scale regulation of epithelium homeostasis

Epithelial tissue provides the body with its first layer of protection against harmful environmental stimuli by enacting the regulatory interplay between a physical barrier preventing the influx of external stimuli and an inflammatory response to the infiltrating stimuli. Importantly, this interdependent regulation occurs on different time scales: the tissue-level barrier permeability is regulated over the course of hours, whereas the cellular-level enzymatic reactions leading to inflammation take place within minutes. This multi-scale regulation is key to the epithelium's function and its dysfunction leads to various diseases. This paper presents a mathematical model of regulatory mechanisms in the epidermal epithelium that includes processes on two different time scales at the cellular and tissue levels. We use this model to investigate the essential regulatory interactions between epidermal barrier integrity and skin inflammation and how their dysfunction leads to atopic dermatitis (AD). Our model exhibits a structure of dual (positive and negative) control at both cellular and tissue levels. We also determined how the variation induced by well-known risk factors for AD can break the balance of the dual control. Our model analysis based on time-scale separation suggests that each risk factor leads to qualitatively different dynamic behaviours of different severity for AD, and that the coincidence of multiple risk factors dramatically increases the fragility of the epithelium's function. The proposed mathematical framework should also be applicable to other inflammatory diseases that have similar time-scale separation and control architectures.

Tuning the dials of Synthetic Biology

Synthetic Biology is the ‘Engineering of Biology’ – it aims to use a forward-engineering design cycle based on specifications, modelling, analysis, experimental implementation, testing and validation to modify natural or design new, synthetic biology systems so that they behave in a predictable fashion. Motivated by the need for truly plug-and-play synthetic biological components, we present a comprehensive review of ways in which the various parts of a biological system can be modified systematically. In particular, we review the list of ‘dials’ that are available to the designer and discuss how they can be modelled, tuned and implemented. The dials are categorized according to whether they operate at the global, transcriptional, translational or post-translational level and the resolution that they operate at. We end this review with a discussion on the relative advantages and disadvantages of some dials over others.

Effect of whole linseed and rumen-protected conjugated linoleic acid enriched diets on feedlot performance, carcass characteristics, and adipose tissue development in young Holstein bulls

Forty-eight young Holstein bulls (slaughtered at 458.6±9.79 kg body weight) were used to evaluate the effect of whole linseed and conjugated linoleic acid (CLA) supplementation on animal performance, adipose tissue development, and carcass characteristics. The animals were fed with one of four isoenergetic and isoproteic diets: control (0% linseed, 0% CLA), linseed (10% linseed, 0% CLA), CLA (0% linseed, 2% CLA), and linseed plus CLA (10% linseed, 2% CLA). Animal performance and carcass characteristics were unaffected by diet composition. Adding linseed or CLA to the concentrate diet did not result in significant differences in adipocyte size and number or lipogenic enzyme activity. However, while the frequency distribution of subcutaneous adipocyte diameters followed a normal distribution, the frequency distribution of intramuscular adipocyte diameters was not normal in any dietary group (skewness coefficients: 0.8, 1.2, 0.9, 0.8 for control, linseed, CLA, and linseed plus CLA, respectively; P<0.05), indicative of adipocyte proliferation in the intramuscular adipose tissue.

Robustness of random graphs based on graph spectra

It has been recently proposed that the robustness of complex networks can be efficiently characterized through the natural connectivity, a spectral property of the graph which corresponds to the average Estrada index. The natural connectivity corresponds to an average eigenvalue calculated from the graph spectrum and can also be interpreted as the Helmholtz free energy of the network. In this article, we explore the use of this index to characterize the robustness of Erdős-Rényi (ER) random graphs, random regular graphs, and regular ring lattices. We show both analytically and numerically that the natural connectivity of ER random graphs increases linearly with the average degree. It is also shown that ER random graphs are more robust than the corresponding random regular graphs with the same number of vertices and edges. However, the relative robustness of ER random graphs and regular ring lattices depends on the average degree and graph size: there is a critical graph size above which regular ring lattices are more robust than random graphs. We use our analytical results to derive this critical graph size as a function of the average degree.

Compound stress response in stomatal closure: a mathematical model of ABA and ethylene interaction in guard cells

BACKGROUND

Stomata are tiny pores in plant leaves that regulate gas and water exchange between the plant and its environment. Abscisic acid and ethylene are two well-known elicitors of stomatal closure when acting independently. However, when stomata are presented with a combination of both signals, they fail to close.

RESULTS

Toshed light on this unexplained behaviour, we have collected time course measurements of stomatal aperture and hydrogen peroxide production in Arabidopsis thaliana guard cells treated with abscisic acid, ethylene, and a combination of both. Our experiments show that stomatal closure is linked to sustained high levels of hydrogen peroxide in guard cells. When treated with a combined dose of abscisic acid and ethylene, guard cells exhibit increased antioxidant activity that reduces hydrogen peroxide levels and precludes closure. We construct a simplified model of stomatal closure derived from known biochemical pathways that captures the experimentally observed behaviour.

CONCLUSIONS

Our experiments and modelling results suggest a distinct role for two antioxidant mechanisms during stomatal closure: a slower, delayed response activated by a single stimulus (abscisic acid 'or' ethylene) and another more rapid 'and' mechanism that is only activated when both stimuli are present. Our model indicates that the presence of this rapid 'and' mechanism in the antioxidant response is key to explain the lack of closure under a combined stimulus.

A modular cell-based biosensor using engineered genetic logic circuits to detect and integrate multiple environmental signals

Cells perceive a wide variety of cellular and environmental signals, which are often processed combinatorially to generate particular phenotypic responses. Here, we employ both single and mixed cell type populations, pre-programmed with engineered modular cell signalling and sensing circuits, as processing units to detect and integrate multiple environmental signals. Based on an engineered modular genetic AND logic gate, we report the construction of a set of scalable synthetic microbe-based biosensors comprising exchangeable sensory, signal processing and actuation modules. These cellular biosensors were engineered using distinct signalling sensory modules to precisely identify various chemical signals, and combinations thereof, with a quantitative fluorescent output. The genetic logic gate used can function as a biological filter and an amplifier to enhance the sensing selectivity and sensitivity of cell-based biosensors. In particular, an Escherichia coli consortium-based biosensor has been constructed that can detect and integrate three environmental signals (arsenic, mercury and copper ion levels) via either its native two-component signal transduction pathways or synthetic signalling sensors derived from other bacteria in combination with a cell-cell communication module. We demonstrate how a modular cell-based biosensor can be engineered predictably using exchangeable synthetic gene circuit modules to sense and integrate multiple-input signals. This study illustrates some of the key practical design principles required for the future application of these biosensors in broad environmental and healthcare areas.

Encoding dynamics for multiscale community detection: Markov time sweeping for the map equation

The detection of community structure in networks is intimately related to finding a concise description of the network in terms of its modules. This notion has been recently exploited by the map equation formalism [Rosvall and Bergstrom, Proc. Natl. Acad. Sci. USA 105, 1118 (2008)] through an information-theoretic description of the process of coding inter- and intracommunity transitions of a random walker in the network at stationarity. However, a thorough study of the relationship between the full Markov dynamics and the coding mechanism is still lacking. We show here that the original map coding scheme, which is both block-averaged and one-step, neglects the internal structure of the communities and introduces an upper scale, the "field-of-view" limit, in the communities it can detect. As a consequence, map is well tuned to detect clique-like communities but can lead to undesirable overpartitioning when communities are far from clique-like. We show that a signature of this behavior is a large compression gap: The map description length is far from its ideal limit. To address this issue, we propose a simple dynamic approach that introduces time explicitly into the map coding through the analysis of the weighted adjacency matrix of the time-dependent multistep transition matrix of the Markov process. The resulting Markov time sweeping induces a dynamical zooming across scales that can reveal (potentially multiscale) community structure above the field-of-view limit, with the relevant partitions indicated by a small compression gap.

Engineering and ethical perspectives in synthetic biology. Rigorous, robust and predictable designs, public engagement and a modern ethical framework are vital to the continued success of synthetic biology

The applications of synthetic biology will involve the release of artificial life forms into the environment. These organisms will present unique safety challenges that need to be addressed by researchers and regulators to win public engagement and support.

Combinatorial stresses kill pathogenic Candida species

Pathogenic microbes exist in dynamic niches and have evolved robust adaptive responses to promote survival in their hosts. The major fungal pathogens of humans, Candida albicans and Candida glabrata, are exposed to a range of environmental stresses in their hosts including osmotic, oxidative and nitrosative stresses. Significant efforts have been devoted to the characterization of the adaptive responses to each of these stresses. In the wild, cells are frequently exposed simultaneously to combinations of these stresses and yet the effects of such combinatorial stresses have not been explored. We have developed a common experimental platform to facilitate the comparison of combinatorial stress responses in C. glabrata and C. albicans. This platform is based on the growth of cells in buffered rich medium at 30°C, and was used to define relatively low, medium and high doses of osmotic (NaCl), oxidative (H ₂O₂) and nitrosative stresses (e.g., dipropylenetriamine (DPTA)-NONOate). The effects of combinatorial stresses were compared with the corresponding individual stresses under these growth conditions. We show for the first time that certain combinations of combinatorial stress are especially potent in terms of their ability to kill C. albicans and C. glabrata and/or inhibit their growth. This was the case for combinations of osmotic plus oxidative stress and for oxidative plus nitrosative stress. We predict that combinatorial stresses may be highly signif cant in host defences against these pathogenic yeasts.

Sensory experience modifies spontaneous state dynamics in a large-scale barrel cortical model

Experimental evidence suggests that spontaneous neuronal activity may shape and be shaped by sensory experience. However, we lack information on how sensory experience modulates the underlying synaptic dynamics and how such modulation influences the response of the network to future events. Here we study whether spike-timing-dependent plasticity (STDP) can mediate sensory-induced modifications in the spontaneous dynamics of a new large-scale model of layers II, III and IV of the rodent barrel cortex. Our model incorporates significant physiological detail, including the types of neurons present, the probabilities and delays of connections, and the STDP profiles at each excitatory synapse. We stimulated the neuronal network with a protocol of repeated sensory inputs resembling those generated by the protraction-retraction motion of whiskers when rodents explore their environment, and studied the changes in network dynamics. By applying dimensionality reduction techniques to the synaptic weight space, we show that the initial spontaneous state is modified by each repetition of the stimulus and that this reverberation of the sensory experience induces long-term, structured modifications in the synaptic weight space. The post-stimulus spontaneous state encodes a memory of the stimulus presented, since a different dynamical response is observed when the network is presented with shuffled stimuli. These results suggest that repeated exposure to the same sensory experience could induce long-term circuitry modifications via 'Hebbian' STDP plasticity.

Quantitative measure of hysteresis for memristors through explicit dynamics

We introduce a mathematical framework for the analysis of the input–output dynamics of externally driven memristors. We show that, under general assumptions, their dynamics comply with a Bernoulli differential equation and hence can be nonlinearly transformed into a formally solvable linear equation. The Bernoulli formalism, which applies to both charge- and flux-controlled memristors when either current or voltage driven, can, in some cases, lead to expressions of the output of the device as an explicit function of the input. We apply our framework to obtain analytical solutions of the i – v characteristics of the recently proposed model of the Hewlett–Packard memristor under three different drives without the need for numerical simulations. Our explicit solutions allow us to identify a dimensionless lumped parameter that combines device-specific parameters with properties of the input drive. This parameter governs the memristive behaviour of the device and, consequently, the amount of hysteresis in the i – v . We proceed further by defining formally a quantitative measure for the hysteresis of the device, for which we obtain explicit formulas in terms of the aforementioned parameter, and we discuss the applicability of the analysis for the design and analysis of memristor devices.

Markov dynamics as a zooming lens for multiscale community detection: non clique-like communities and the field-of-view limit

In recent years, there has been a surge of interest in community detection algorithms for complex networks. A variety of computational heuristics, some with a long history, have been proposed for the identification of communities or, alternatively, of good graph partitions. In most cases, the algorithms maximize a particular objective function, thereby finding the 'right' split into communities. Although a thorough comparison of algorithms is still lacking, there has been an effort to design benchmarks, i.e., random graph models with known community structure against which algorithms can be evaluated. However, popular community detection methods and benchmarks normally assume an implicit notion of community based on clique-like subgraphs, a form of community structure that is not always characteristic of real networks. Specifically, networks that emerge from geometric constraints can have natural non clique-like substructures with large effective diameters, which can be interpreted as long-range communities. In this work, we show that long-range communities escape detection by popular methods, which are blinded by a restricted 'field-of-view' limit, an intrinsic upper scale on the communities they can detect. The field-of-view limit means that long-range communities tend to be overpartitioned. We show how by adopting a dynamical perspective towards community detection [1], [2], in which the evolution of a Markov process on the graph is used as a zooming lens over the structure of the network at all scales, one can detect both clique- or non clique-like communities without imposing an upper scale to the detection. Consequently, the performance of algorithms on inherently low-diameter, clique-like benchmarks may not always be indicative of equally good results in real networks with local, sparser connectivity. We illustrate our ideas with constructive examples and through the analysis of real-world networks from imaging, protein structures and the power grid, where a multiscale structure of non clique-like communities is revealed.

Squeeze-and-breathe evolutionary Monte Carlo optimization with local search acceleration and its application to parameter fitting

Estimating parameters from data is a key stage of the modelling process, particularly in biological systems where many parameters need to be estimated from sparse and noisy datasets. Over the years, a variety of heuristics have been proposed to solve this complex optimization problem, with good results in some cases yet with limitations in the biological setting. In this work, we develop an algorithm for model parameter fitting that combines ideas from evolutionary algorithms, sequential Monte Carlo and direct search optimization. Our method performs well even when the order of magnitude and/or the range of the parameters is unknown. The method refines iteratively a sequence of parameter distributions through local optimization combined with partial resampling from a historical prior defined over the support of all previous iterations. We exemplify our method with biological models using both simulated and real experimental data and estimate the parameters efficiently even in the absence of a priori knowledge about the parameters.

Flow graphs: interweaving dynamics and structure

The behavior of complex systems is determined not only by the topological organization of their interconnections but also by the dynamical processes taking place among their constituents. A faithful modeling of the dynamics is essential because different dynamical processes may be affected very differently by network topology. A full characterization of such systems thus requires a formalization that encompasses both aspects simultaneously, rather than relying only on the topological adjacency matrix. To achieve this, we introduce the concept of flow graphs, namely weighted networks where dynamical flows are embedded into the link weights. Flow graphs provide an integrated representation of the structure and dynamics of the system, which can then be analyzed with standard tools from network theory. Conversely, a structural network feature of our choice can also be used as the basis for the construction of a flow graph that will then encompass a dynamics biased by such a feature. We illustrate the ideas by focusing on the mathematical properties of generic linear processes on complex networks that can be represented as biased random walks and their dual consensus dynamics, and show how our framework improves our understanding of these processes.

Transient dynamics around unstable periodic orbits in the generalized repressilator model

We study the temporal dynamics of the generalized repressilator, a network of coupled repressing genes arranged in a directed ring topology, and give analytical conditions for the emergence of a finite sequence of unstable periodic orbits that lead to reachable long-lived oscillating transients. Such transients dominate the finite time horizon dynamics that is relevant in confined, noisy environments such as bacterial cells (see our previous work [Strelkowa and Barahona, J. R. Soc. Interface 7, 1071 (2010)]), and are therefore of interest for bioengineering and synthetic biology. We show that the family of unstable orbits possesses spatial symmetries and can also be understood in terms of traveling wave solutions of kink-like topological defects. The long-lived oscillatory transients correspond to the propagation of quasistable two-kink configurations that unravel over a long time. We also assess the similarities between the generalized repressilator model and other unidirectionally coupled electronic systems, such as magnetic flux gates, which have been implemented experimentally.

Crowding-induced anisotropic transport modulates reaction kinetics in nanoscale porous media

We quantify the emergence of persistent anisotropy in the diffusion of spherical tracer particles through a nanoscale porous medium composed of a uniform distribution of purely symmetric crowding particles. We focus on the interior of a biological cell as an example of such a medium and find that diffusion is highly directional for distances comparable to the size of some organelles. We use a geometrical procedure that avoids the standard orientational averaging to quantify the anisotropy of diffusive paths and show that the point source distributions are predominantly of prolate ellipsoidal shape as a result of local volume exclusion. This geometrical symmetry breaking strongly skews the distribution of kinetic rates of diffusion-limited reactions toward small values, leading to the result that, for short to intermediate times, almost 80% of the rates measured in an ensemble of heterogeneous media are smaller than the expected rate in an ideal homogeneous medium of similar excluded volume fraction. This crowding-induced modulation may have implications for our understanding and measurement of diffusion-controlled intracellular reaction kinetics and for experimental nanotechnology applications, such as nanoparticle-based bioimaging and drug delivery, where diffusion plays an important role.

Switchable genetic oscillator operating in quasi-stable mode

Ring topologies of repressing genes have qualitatively different long-term dynamics if the number of genes is odd (they oscillate) or even (they exhibit bistability). However, these attractors may not fully explain the observed behaviour in transient and stochastic environments such as the cell. We show here that even repressilators possess quasi-stable, travelling wave periodic solutions that are reachable, long-lived and robust to parameter changes. These solutions underlie the sustained oscillations observed in even rings in the stochastic regime, even if these circuits are expected to behave as switches. The existence of such solutions can also be exploited for control purposes: operation of the system around the quasi-stable orbit allows us to turn on and off the oscillations reliably and on demand. We illustrate these ideas with a simple protocol based on optical interference that can induce oscillations robustly both in the stochastic and deterministic regimes.