Network motifs: structure does not determine function

BACKGROUND: A number of publications have recently examined the occurrence and properties of the feed-forward motif in a variety of networks, including those that are of interest in genome biology, such as gene networks. The present work looks in some detail at the dynamics of the bi-fan motif, using systems of ordinary differential equations to model the populations of transcription factors, mRNA and protein, with the aim of extending our understanding of what appear to be important building blocks of gene network structure. RESULTS: We develop an ordinary differential equation model of the bi-fan motif and analyse variants of the motif corresponding to its behaviour under various conditions. In particular, we examine the effects of different steady and pulsed inputs to five variants of the bifan motif, based on evidence in the literature of bifan motifs found in Saccharomyces cerevisiae (commonly known as baker's yeast). Using this model, we characterize the dynamical behaviour of the bi-fan motif for a wide range of biologically plausible parameters and configurations. We find that there is no characteristic behaviour for the motif, and with the correct choice of parameters and of internal structure, very different, indeed even opposite behaviours may be obtained. CONCLUSION: Even with this relatively simple model, the bi-fan motif can exhibit a wide range of dynamical responses. This suggests that it is difficult to gain significant insights into biological function simply by considering the connection architecture of a gene network, or its decomposition into simple structural motifs. It is necessary to supplement such structural information by kinetic parameters, or dynamic time series experimental data, both of which are currently difficult to obtain

Understanding the Slow Depletion of Memory CD4+ T Cells in HIV Infection

Using a simple mathematical model, Andrew Yates and colleagues show that a runaway cycle of T cell activation and infection cannot explain the slow rate of CD4 decline during chronic HIV infection

Dissection of a complex transcriptional response using genome-wide transcriptional modelling

Modern genomics technologies generate huge data sets creating a demand for systems level, experimentally verified, analysis techniques. We examined the transcriptional response to DNA damage in a human T cell line (MOLT4) using microarrays. By measuring both mRNA accumulation and degradation over a short time course, we were able to construct a mechanistic model of the transcriptional response. The model predicted three dominant transcriptional activity profiles—an early response controlled by NFκB and c-Jun, a delayed response controlled by p53, and a late response related to cell cycle re-entry. The method also identified, with defined confidence limits, the transcriptional targets associated with each activity. Experimental inhibition of NFκB, c-Jun and p53 confirmed that target predictions were accurate. Model predictions directly explained 70% of the 200 most significantly upregulated genes in the DNA-damage response. Genome-wide transcriptional modelling (GWTM) requires no prior knowledge of either transcription factors or their targets. GWTM is an economical and effective method for identifying the main transcriptional activators in a complex response and confidently predicting their targets

Abnormal Cortical Development after Premature Birth Shown by Altered Allometric Scaling of Brain Growth

BACKGROUND: We postulated that during ontogenesis cortical surface area and cerebral volume are related by a scaling law whose exponent gives a quantitative measure of cortical development. We used this approach to investigate the hypothesis that premature termination of the intrauterine environment by preterm birth reduces cortical development in a dose-dependent manner, providing a neural substrate for functional impairment. METHODS AND FINDINGS: We analyzed 274 magnetic resonance images that recorded brain growth from 23 to 48 wk of gestation in 113 extremely preterm infants born at 22 to 29 wk of gestation, 63 of whom underwent neurodevelopmental assessment at a median age of 2 y. Cortical surface area was related to cerebral volume by a scaling law with an exponent of 1.29 (95% confidence interval, 1.25–1.33), which was proportional to later neurodevelopmental impairment. Increasing prematurity and male gender were associated with a lower scaling exponent (p < 0.0001) independent of intrauterine or postnatal somatic growth. CONCLUSIONS: Human brain growth obeys an allometric scaling relation that is disrupted by preterm birth in a dose-dependent, sexually dimorphic fashion that directly parallels the incidence of neurodevelopmental impairments in preterm infants. This result focuses attention on brain growth and cortical development during the weeks following preterm delivery as a neural substrate for neurodevelopmental impairment after premature delivery

Nonidentifiability of the Source of Intrinsic Noise in Gene Expression from Single-Burst Data

Over the last few years, experimental data on the fluctuations in gene activity between individual cells and within the same cell over time have confirmed that gene expression is a “noisy” process. This variation is in part due to the small number of molecules taking part in some of the key reactions that are involved in gene expression. One of the consequences of this is that protein production often occurs in bursts, each due to a single promoter or transcription factor binding event. Recently, the distribution of the number of proteins produced in such bursts has been experimentally measured, offering a unique opportunity to study the relative importance of different sources of noise in gene expression. Here, we provide a derivation of the theoretical probability distribution of these bursts for a wide variety of different models of gene expression. We show that there is a good fit between our theoretical distribution and that obtained from two different published experimental datasets. We then prove that, irrespective of the details of the model, the burst size distribution is always geometric and hence determined by a single parameter. Many different combinations of the biochemical rates for the constituent reactions of both transcription and translation will therefore lead to the same experimentally observed burst size distribution. It is thus impossible to identify different sources of fluctuations purely from protein burst size data or to use such data to estimate all of the model parameters. We explore methods of inferring these values when additional types of experimental data are available

Reconstruction of cell population dynamics using CFSE

Background: Quantifying cell division and death is central to many studies in the biological sciences. The fluorescent dye CFSE allows the tracking of cell division in vitro and in vivo and provides a rich source of information with which to test models of cell kinetics. Cell division and death have a stochastic component at the single-cell level, and the probabilities of these occurring in any given time interval may also undergo systematic variation at a population level. This gives rise to heterogeneity in proliferating cell populations. Branching processes provide a natural means of describing this behaviour. Results: We present a likelihood-based method for estimating the parameters of branching process models of cell kinetics using CFSE-labeling experiments, and demonstrate its validity using synthetic and experimental datasets. Performing inference and model comparison with real CFSE data presents some statistical problems and we suggest methods of dealing with them. Conclusion: The approach we describe here can be used to recover the (potentially variable) division and death rates of any cell population for which division tracking information is available

Measurement of the B0-anti-B0-Oscillation Frequency with Inclusive Dilepton Events

The $B^0$-$\bar B^0$ oscillation frequency has been measured with a sample of 23 million \B\bar B pairs collected with the BABAR detector at the PEP-II asymmetric B Factory at SLAC. In this sample, we select events in which both B mesons decay semileptonically and use the charge of the leptons to identify the flavor of each B meson. A simultaneous fit to the decay time difference distributions for opposite- and same-sign dilepton events gives $\Delta m_d = 0.493 \pm 0.012{(stat)}\pm 0.009{(syst)}$ ps$^{-1}$.Comment: 7 pages, 1 figure, submitted to Physical Review Letter