Optimizing information flow in small genetic networks. I

In order to survive, reproduce and (in multicellular organisms) differentiate, cells must control the concentrations of the myriad different proteins that are encoded in the genome. The precision of this control is limited by the inevitable randomness of individual molecular events. Here we explore how cells can maximize their control power in the presence of these physical limits; formally, we solve the theoretical problem of maximizing the information transferred from inputs to outputs when the number of available molecules is held fixed. We start with the simplest version of the problem, in which a single transcription factor protein controls the readout of one or more genes by binding to DNA. We further simplify by assuming that this regulatory network operates in steady state, that the noise is small relative to the available dynamic range, and that the target genes do not interact. Even in this simple limit, we find a surprisingly rich set of optimal solutions. Importantly, for each locally optimal regulatory network, all parameters are determined once the physical constraints on the number of available molecules are specified. Although we are solving an over--simplified version of the problem facing real cells, we see parallels between the structure of these optimal solutions and the behavior of actual genetic regulatory networks. Subsequent papers will discuss more complete versions of the problem

Challenges in quantifying changes in the global water cycle

Human influences have likely already impacted the large-scale water cycle but natural variability and observational uncertainty are substantial. It is essential to maintain and improve observational capabilities to better characterize changes. Understanding observed changes to the global water cycle is key to predicting future climate changes and their impacts. While many datasets document crucial variables such as precipitation, ocean salinity, runoff, and humidity, most are uncertain for determining long-term changes. In situ networks provide long time-series over land but are sparse in many regions, particularly the tropics. Satellite and reanalysis datasets provide global coverage, but their long-term stability is lacking. However, comparisons of changes among related variables can give insights into the robustness of observed changes. For example, ocean salinity, interpreted with an understanding of ocean processes, can help cross-validate precipitation. Observational evidence for human influences on the water cycle is emerging, but uncertainties resulting from internal variability and observational errors are too large to determine whether the observed and simulated changes are consistent. Improvements to the in situ and satellite observing networks that monitor the changing water cycle are required, yet continued data coverage is threatened by funding reductions. Uncertainty both in the role of anthropogenic aerosols, and due to large climate variability presently limits confidence in attribution of observed changes

Exploring the Contextual Sensitivity of Factors that Determine Cell-to-Cell Variability in Receptor-Mediated Apoptosis

Stochastic fluctuations in gene expression give rise to cell-to-cell variability in protein levels which can potentially cause variability in cellular phenotype. For TRAIL (TNF-related apoptosis-inducing ligand) variability manifests itself as dramatic differences in the time between ligand exposure and the sudden activation of the effector caspases that kill cells. However, the contribution of individual proteins to phenotypic variability has not been explored in detail. In this paper we use feature-based sensitivity analysis as a means to estimate the impact of variation in key apoptosis regulators on variability in the dynamics of cell death. We use Monte Carlo sampling from measured protein concentration distributions in combination with a previously validated ordinary differential equation model of apoptosis to simulate the dynamics of receptor-mediated apoptosis. We find that variation in the concentrations of some proteins matters much more than variation in others and that precisely which proteins matter depends both on the concentrations of other proteins and on whether correlations in protein levels are taken into account. A prediction from simulation that we confirm experimentally is that variability in fate is sensitive to even small increases in the levels of Bcl-2. We also show that sensitivity to Bcl-2 levels is itself sensitive to the levels of interacting proteins. The contextual dependency is implicit in the mathematical formulation of sensitivity, but our data show that it is also important for biologically relevant parameter values. Our work provides a conceptual and practical means to study and understand the impact of cell-to-cell variability in protein expression levels on cell fate using deterministic models and sampling from parameter distributions

Allosteric Indole Amide Inhibitors of p97: Identification of a Novel Probe of the Ubiquitin Pathway

A high-throughput screen to discover inhibitors of p97 ATPase activity identified an indole amide that bound to an allosteric site of the protein. Medicinal chemistry optimization led to improvements in potency and solubility. Indole amide 3 represents a novel uncompetitive inhibitor with excellent physical and pharmaceutical properties that can be used as a starting point for drug discovery efforts

Massively parallel fitness profiling reveals multiple novel enzymes in Pseudomonas putida lysine metabolism

P. putida lysine metabolism can produce multiple commodity chemicals, conferring great biotechnological value. Despite much research, the connection of lysine catabolism to central metabolism in P. putida remained undefined. Here, we used random barcode transposon sequencing to fill the gaps of lysine metabolism in P. putida. We describe a route of 2-oxoadipate (2OA) catabolism, which utilizes DUF1338-containing protein P. putida 5260 (PP_5260) in bacteria. Despite its prevalence in many domains of life, DUF1338-containing proteins have had no known biochemical function. We demonstrate that PP_5260 is a metalloenzyme which catalyzes an unusual route of decarboxylation of 2OA to d-2-hydroxyglutarate (d-2HG). Our screen also identified a recently described novel glutarate metabolic pathway. We validate previous results and expand the understanding of glutarate hydroxylase CsiD by showing that can it use either 2OA or 2KG as a cosubstrate. Our work demonstrated that biological novelty can be rapidly identified using unbiased experimental genetics and that RB-TnSeq can be used to rapidly validate previous results.Despite intensive study for 50 years, the biochemical and genetic links between lysine metabolism and central metabolism in Pseudomonas putida remain unresolved. To establish these biochemical links, we leveraged random barcode transposon sequencing (RB-TnSeq), a genome-wide assay measuring the fitness of thousands of genes in parallel, to identify multiple novel enzymes in both l- and d-lysine metabolism. We first describe three pathway enzymes that catabolize l-2-aminoadipate (l-2AA) to 2-ketoglutarate (2KG), connecting d-lysine to the TCA cycle. One of these enzymes, P. putida 5260 (PP_5260), contains a DUF1338 domain, representing a family with no previously described biological function. Our work also identified the recently described coenzyme A (CoA)-independent route of l-lysine degradation that results in metabolization to succinate. We expanded on previous findings by demonstrating that glutarate hydroxylase CsiD is promiscuous in its 2-oxoacid selectivity. Proteomics of selected pathway enzymes revealed that expression of catabolic genes is highly sensitive to the presence of particular pathway metabolites, implying intensive local and global regulation. This work demonstrated the utility of RB-TnSeq for discovering novel metabolic pathways in even well-studied bacteria, as well as its utility a powerful tool for validating previous research