Rules for biological regulation based on error minimization

The control of gene expression involves complex mechanisms that show large variation in design. For example, genes can be turned on either by the binding of an activator (positive control) or the unbinding of a repressor (negative control). What determines the choice of mode of control for each gene? This study proposes rules for gene regulation based on the assumption that free regulatory sites are exposed to nonspecific binding errors, whereas sites bound to their cognate regulators are protected from errors. Hence, the selected mechanisms keep the sites bound to their designated regulators for most of the time, thus minimizing fitness-reducing errors. This offers an explanation of the empirically demonstrated Savageau demand rule: Genes that are needed often in the natural environment tend to be regulated by activators, and rarely needed genes tend to be regulated by repressors; in both cases, sites are bound for most of the time, and errors are minimized. The fitness advantage of error minimization appears to be readily selectable. The present approach can also generate rules for multi-regulator systems. The error-minimization framework raises several experimentally testable hypotheses. It may also apply to other biological regulation systems, such as those involving protein-protein interactions.Comment: biological physics, complex networks, systems biology, transcriptional regulation http://www.weizmann.ac.il/complex/tlusty/papers/PNAS2006.pdf http://www.pnas.org/content/103/11/3999.ful

The Role of Regulated mRNA Stability in Establishing Bicoid Morphogen Gradient in Drosophila Embryonic Development

The Bicoid morphogen is amongst the earliest triggers of differential spatial pattern of gene expression and subsequent cell fate determination in the embryonic development of Drosophila. This maternally deposited morphogen is thought to diffuse in the embryo, establishing a concentration gradient which is sensed by downstream genes. In most model based analyses of this process, the translation of the bicoid mRNA is thought to take place at a fixed rate from the anterior pole of the embryo and a supply of the resulting protein at a constant rate is assumed. Is this process of morphogen generation a passive one as assumed in the modelling literature so far, or would available data support an alternate hypothesis that the stability of the mRNA is regulated by active processes? We introduce a model in which the stability of the maternal mRNA is regulated by being held constant for a length of time, followed by rapid degradation. With this more realistic model of the source, we have analysed three computational models of spatial morphogen propagation along the anterior-posterior axis: (a) passive diffusion modelled as a deterministic differential equation, (b) diffusion enhanced by a cytoplasmic flow term; and (c) diffusion modelled by stochastic simulation of the corresponding chemical reactions. Parameter estimation on these models by matching to publicly available data on spatio-temporal Bicoid profiles suggests strong support for regulated stability over either a constant supply rate or one where the maternal mRNA is permitted to degrade in a passive manner

Regulatory control and the costs and benefits of biochemical noise

Experiments in recent years have vividly demonstrated that gene expression can be highly stochastic. How protein concentration fluctuations affect the growth rate of a population of cells, is, however, a wide open question. We present a mathematical model that makes it possible to quantify the effect of protein concentration fluctuations on the growth rate of a population of genetically identical cells. The model predicts that the population's growth rate depends on how the growth rate of a single cell varies with protein concentration, the variance of the protein concentration fluctuations, and the correlation time of these fluctuations. The model also predicts that when the average concentration of a protein is close to the value that maximizes the growth rate, fluctuations in its concentration always reduce the growth rate. However, when the average protein concentration deviates sufficiently from the optimal level, fluctuations can enhance the growth rate of the population, even when the growth rate of a cell depends linearly on the protein concentration. The model also shows that the ensemble or population average of a quantity, such as the average protein expression level or its variance, is in general not equal to its time average as obtained from tracing a single cell and its descendants. We apply our model to perform a cost-benefit analysis of gene regulatory control. Our analysis predicts that the optimal expression level of a gene regulatory protein is determined by the trade-off between the cost of synthesizing the regulatory protein and the benefit of minimizing the fluctuations in the expression of its target gene. We discuss possible experiments that could test our predictions.Comment: Revised manuscript;35 pages, 4 figures, REVTeX4; to appear in PLoS Computational Biolog

Modeling recursive RNA interference.

An important application of the RNA interference (RNAi) pathway is its use as a small RNA-based regulatory system commonly exploited to suppress expression of target genes to test their function in vivo. In several published experiments, RNAi has been used to inactivate components of the RNAi pathway itself, a procedure termed recursive RNAi in this report. The theoretical basis of recursive RNAi is unclear since the procedure could potentially be self-defeating, and in practice the effectiveness of recursive RNAi in published experiments is highly variable. A mathematical model for recursive RNAi was developed and used to investigate the range of conditions under which the procedure should be effective. The model predicts that the effectiveness of recursive RNAi is strongly dependent on the efficacy of RNAi at knocking down target gene expression. This efficacy is known to vary highly between different cell types, and comparison of the model predictions to published experimental data suggests that variation in RNAi efficacy may be the main cause of discrepancies between published recursive RNAi experiments in different organisms. The model suggests potential ways to optimize the effectiveness of recursive RNAi both for screening of RNAi components as well as for improved temporal control of gene expression in switch off-switch on experiments

Outcome of ATP-based tumor chemosensitivity assay directed chemotherapy in heavily pre-treated recurrent ovarian carcinoma

BACKGROUND: We wished to evaluate the clinical response following ATP-Tumor Chemosensitivity Assay (ATP-TCA) directed salvage chemotherapy in a series of UK patients with advanced ovarian cancer. The results are compared with that of a similar assay used in a different country in terms of evaluability and clinical endpoints. METHODS: From November 1998 to November 2001, 46 patients with pre-treated, advanced ovarian cancer were given a total of 56 courses of chemotherapy based on in-vitro ATP-TCA responses obtained from fresh tumor samples or ascites. Forty-four patients were evaluable for results. Of these, 18 patients had clinically platinum resistant disease (relapse < 6 months after first course of chemotherapy). There was evidence of cisplatin resistance in 31 patients from their first ATP-TCA. Response to treatment was assessed by radiology, clinical assessment and tumor marker level (CA 125). RESULTS: The overall response rate was 59% (33/56) per course of chemotherapy, including 12 complete responses, 21 partial responses, 6 with stable disease, and 15 with progressive disease. Two patients were not evaluable for response having received just one cycle of chemotherapy: if these were excluded the response rate is 61%. Fifteen patients are still alive. Median progression free survival (PFS) was 6.6 months per course of chemotherapy; median overall survival (OAS) for each patient following the start of TCA-directed therapy was 10.4 months (95% confidence interval 7.9-12.8 months). CONCLUSION: The results show similar response rates to previous studies using ATP-TCA directed therapy in recurrent ovarian cancer. The assay shows high evaluability and this study adds weight to the reproducibility of results from different centre

A Genome-Wide Analysis of Promoter-Mediated Phenotypic Noise in Escherichia coli

Gene expression is subject to random perturbations that lead to fluctuations in the rate of protein production. As a consequence, for any given protein, genetically identical organisms living in a constant environment will contain different amounts of that particular protein, resulting in different phenotypes. This phenomenon is known as “phenotypic noise.” In bacterial systems, previous studies have shown that, for specific genes, both transcriptional and translational processes affect phenotypic noise. Here, we focus on how the promoter regions of genes affect noise and ask whether levels of promoter-mediated noise are correlated with genes' functional attributes, using data for over 60% of all promoters in Escherichia coli. We find that essential genes and genes with a high degree of evolutionary conservation have promoters that confer low levels of noise. We also find that the level of noise cannot be attributed to the evolutionary time that different genes have spent in the genome of E. coli. In contrast to previous results in eukaryotes, we find no association between promoter-mediated noise and gene expression plasticity. These results are consistent with the hypothesis that, in bacteria, natural selection can act to reduce gene expression noise and that some of this noise is controlled through the sequence of the promoter region alon

The interplay of intrinsic and extrinsic bounded noises in genetic networks

After being considered as a nuisance to be filtered out, it became recently clear that biochemical noise plays a complex role, often fully functional, for a genetic network. The influence of intrinsic and extrinsic noises on genetic networks has intensively been investigated in last ten years, though contributions on the co-presence of both are sparse. Extrinsic noise is usually modeled as an unbounded white or colored gaussian stochastic process, even though realistic stochastic perturbations are clearly bounded. In this paper we consider Gillespie-like stochastic models of nonlinear networks, i.e. the intrinsic noise, where the model jump rates are affected by colored bounded extrinsic noises synthesized by a suitable biochemical state-dependent Langevin system. These systems are described by a master equation, and a simulation algorithm to analyze them is derived. This new modeling paradigm should enlarge the class of systems amenable at modeling. We investigated the influence of both amplitude and autocorrelation time of a extrinsic Sine-Wiener noise on: $(i)$ the Michaelis-Menten approximation of noisy enzymatic reactions, which we show to be applicable also in co-presence of both intrinsic and extrinsic noise, $(ii)$ a model of enzymatic futile cycle and $(iii)$ a genetic toggle switch. In $(ii)$ and $(iii)$ we show that the presence of a bounded extrinsic noise induces qualitative modifications in the probability densities of the involved chemicals, where new modes emerge, thus suggesting the possibile functional role of bounded noises

Effect of promoter architecture on the cell-to-cell variability in gene expression

According to recent experimental evidence, the architecture of a promoter, defined as the number, strength and regulatory role of the operators that control the promoter, plays a major role in determining the level of cell-to-cell variability in gene expression. These quantitative experiments call for a corresponding modeling effort that addresses the question of how changes in promoter architecture affect noise in gene expression in a systematic rather than case-by-case fashion. In this article, we make such a systematic investigation, based on a simple microscopic model of gene regulation that incorporates stochastic effects. In particular, we show how operator strength and operator multiplicity affect this variability. We examine different modes of transcription factor binding to complex promoters (cooperative, independent, simultaneous) and how each of these affects the level of variability in transcription product from cell-to-cell. We propose that direct comparison between in vivo single-cell experiments and theoretical predictions for the moments of the probability distribution of mRNA number per cell can discriminate between different kinetic models of gene regulation.Comment: 35 pages, 6 figures, Submitte

Clinical outcome and risk factors for recurrence in borderline ovarian tumours

We investigated the long-term prognosis of borderline ovarian tumours and determined risk factors for recurrence. One hundred and twenty-one borderline ovarian tumours treated between 1994 and 2003 at the participating institutions in the Tohoku Gynecologic Cancer Unit were retrospectively investigated for clinical stage, histopathological subtype, surgical technique, postoperative chemotherapy, the presence or absence of recurrence, and prognosis. The median follow-up period was 57 months (1–126 months). One hundred and nine cases (90.6%) were at clinical stage I. The histopathological subtypes consisted of 91 cases of mucinous tumour (75.2%), 27 cases of serous tumour (22.3%), and three cases of endometrioid tumour. Conservative surgery was used in 53 cases (43.8%), radical surgery in 68 cases (56.2%), a staging laparotomy in 43 cases (35.5%), and postoperative adjuvant therapy in 30 cases (24.8%). Recurrence was found in eight cases, but no tumour-related deaths were reported. Although no significant difference in disease-free survival rate was seen between different clinical stages, the difference in disease-free survival rate between serous and nonserous (mucinous and endometrioid) types was significant (P<0.05). The 10-year disease-free survival rate was 89.1% for the radical surgery group and 57.4% for the conservative surgery group – this difference was significant (P<0.05). In the conservative surgery group, cystectomy and serous tumour were independent risk factors for recurrence. Although recurrence was observed, the long-term prognosis of borderline ovarian tumour was favourable, without tumour-related deaths. Considering the favourable prognosis, conservative surgery can be chosen as far as the patient has a nonserous tumour and receive adnexectomy. However, in cases of serous type and/or receiving cystectomy special care should be given as relative risk rates of recurrence elevate by 2–4-folds