A new multicompartmental reaction-diffusion modeling method links transient membrane attachment of E. coli MinE to E-ring formation

Many important cellular processes are regulated by reaction-diffusion (RD) of molecules that takes place both in the cytoplasm and on the membrane. To model and analyze such multicompartmental processes, we developed a lattice-based Monte Carlo method, Spatiocyte that supports RD in volume and surface compartments at single molecule resolution. Stochasticity in RD and the excluded volume effect brought by intracellular molecular crowding, both of which can significantly affect RD and thus, cellular processes, are also supported. We verified the method by comparing simulation results of diffusion, irreversible and reversible reactions with the predicted analytical and best available numerical solutions. Moreover, to directly compare the localization patterns of molecules in fluorescence microscopy images with simulation, we devised a visualization method that mimics the microphotography process by showing the trajectory of simulated molecules averaged according to the camera exposure time. In the rod-shaped bacterium _Escherichia coli_, the division site is suppressed at the cell poles by periodic pole-to-pole oscillations of the Min proteins (MinC, MinD and MinE) arising from carefully orchestrated RD in both cytoplasm and membrane compartments. Using Spatiocyte we could model and reproduce the _in vivo_ MinDE localization dynamics by accounting for the established properties of MinE. Our results suggest that the MinE ring, which is essential in preventing polar septation, is largely composed of MinE that is transiently attached to the membrane independently after recruited by MinD. Overall, Spatiocyte allows simulation and visualization of complex spatial and reaction-diffusion mediated cellular processes in volumes and surfaces. As we showed, it can potentially provide mechanistic insights otherwise difficult to obtain experimentally

Bayesian inference of biochemical kinetic parameters using the linear noise approximation

Background Fluorescent and luminescent gene reporters allow us to dynamically quantify changes in molecular species concentration over time on the single cell level. The mathematical modeling of their interaction through multivariate dynamical models requires the deveopment of effective statistical methods to calibrate such models against available data. Given the prevalence of stochasticity and noise in biochemical systems inference for stochastic models is of special interest. In this paper we present a simple and computationally efficient algorithm for the estimation of biochemical kinetic parameters from gene reporter data. Results We use the linear noise approximation to model biochemical reactions through a stochastic dynamic model which essentially approximates a diffusion model by an ordinary differential equation model with an appropriately defined noise process. An explicit formula for the likelihood function can be derived allowing for computationally efficient parameter estimation. The proposed algorithm is embedded in a Bayesian framework and inference is performed using Markov chain Monte Carlo. Conclusion The major advantage of the method is that in contrast to the more established diffusion approximation based methods the computationally costly methods of data augmentation are not necessary. Our approach also allows for unobserved variables and measurement error. The application of the method to both simulated and experimental data shows that the proposed methodology provides a useful alternative to diffusion approximation based methods

Geometry-controlled kinetics

It has long been appreciated that transport properties can control reaction kinetics. This effect can be characterized by the time it takes a diffusing molecule to reach a target -- the first-passage time (FPT). Although essential to quantify the kinetics of reactions on all time scales, determining the FPT distribution was deemed so far intractable. Here, we calculate analytically this FPT distribution and show that transport processes as various as regular diffusion, anomalous diffusion, diffusion in disordered media and in fractals fall into the same universality classes. Beyond this theoretical aspect, this result changes the views on standard reaction kinetics. More precisely, we argue that geometry can become a key parameter so far ignored in this context, and introduce the concept of "geometry-controlled kinetics". These findings could help understand the crucial role of spatial organization of genes in transcription kinetics, and more generally the impact of geometry on diffusion-limited reactions.Comment: Submitted versio

Community-based adaptation research in the Canadian Arctic

Community-based adaptation (CBA) has emerged over the last decade as an approach to empowering communities to plan for and cope with the impacts of climate change. While such approaches have been widely advocated, few have critically examined the tensions and challenges that CBA brings. Responding to this gap, this article critically examines the use of CBA approaches with Inuit communities in Canada. We suggest that CBA holds significant promise to make adaptation research more democratic and responsive to local needs, providing a basis for developing locally appropriate adaptations based on local/indigenous and Western knowledge. Yet, we argue that CBA is not a panacea, and its common portrayal as such obscures its limitations, nuances, and challenges. Indeed, if uncritically adopted, CBA can potentially lead to maladaptation, may be inappropriate in some instances, can legitimize outside intervention and control, and may further marginalize communities. We identify responsibilities for researchers engaging in CBA work to manage these challenges, emphasizing the centrality of how knowledge is generated, the need for project flexibility and openness to change, and the importance of ensuring partnerships between researchers and communities are transparent. Researchers also need to be realistic about what CBA can achieve, and should not assume that research has a positive role to play in community adaptation just because it utilizes participatory approaches

Analytic philosophy for biomedical research: the imperative of applying yesterday's timeless messages to today's impasses

The mantra that "the best way to predict the future is to invent it" (attributed to the computer scientist Alan Kay) exemplifies some of the expectations from the technical and innovative sides of biomedical research at present. However, for technical advancements to make real impacts both on patient health and genuine scientific understanding, quite a number of lingering challenges facing the entire spectrum from protein biology all the way to randomized controlled trials should start to be overcome. The proposal in this chapter is that philosophy is essential in this process. By reviewing select examples from the history of science and philosophy, disciplines which were indistinguishable until the mid-nineteenth century, I argue that progress toward the many impasses in biomedicine can be achieved by emphasizing theoretical work (in the true sense of the word 'theory') as a vital foundation for experimental biology. Furthermore, a philosophical biology program that could provide a framework for theoretical investigations is outlined

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

The potential for sand dams to increase the adaptive capacity of East African drylands to climate change

Drylands are home to more than two billion people and are characterised by frequent, severe droughts. Such extreme events are expected to be exacerbated in the near future by climate change. A potentially simple and cost-effective mitigation measure against drought periods is sand dams. This little-known technology aims to promote subsoil rainwater storage to support dryland agro-ecosystems. To date, there is little long-term empirical analysis that tests the effectiveness of this approach during droughts. This study addresses this shortcoming by utilising multi-year satellite imagery to monitor the effect of droughts at sand dam locations. A time series of satellite images was analysed to compare vegetation at sand dam sites and control sites over selected periods of drought, using the normalised difference vegetation index. The results show that vegetation biomass was consistently and significantly higher at sand dam sites during periods of extended droughts. It is also shown that vegetation at sand dam sites recovers more quickly from drought. The observed findings corroborate modelling-based research which identified related impacts on ground water, land cover, and socio-economic indicators. Using past periods of drought as an analogue to future climate change conditions, this study indicates that sand dams have potential to increase adaptive capacity and resilience to climate change in drylands. It therefore can be concluded that sand dams enhance the resilience of marginal environments and increase the adaptive capacity of drylands. Sand dams can therefore be a promising adaptation response to the impacts of future climate change on drylands

MicroRNAs can generate thresholds in target gene expression

MicroRNAs (miRNAs) are short, highly conserved noncoding RNA molecules that repress gene expression in a sequence-dependent manner. We performed single-cell measurements using quantitative fluorescence microscopy and flow cytometry to monitor a target gene's protein expression in the presence and absence of regulation by miRNA. We find that although the average level of repression is modest, in agreement with previous population-based measurements, the repression among individual cells varies dramatically. In particular, we show that regulation by miRNAs establishes a threshold level of target mRNA below which protein production is highly repressed. Near this threshold, protein expression responds sensitively to target mRNA input, consistent with a mathematical model of molecular titration. These results show that miRNAs can act both as a switch and as a fine-tuner of gene expression.National Institutes of Health (U.S.). Director's Pioneer Award (1DP1OD003936)National Cancer Institute (U.S.). Physical Sciences-Oncology Center (U54CA143874)United States. Public Health Service (Grant R01-CA133404)United States. Public Health Service (Grant R01-GM34277)National Cancer Institute (U.S.) (PO1-CA42063)National Cancer Institute (U.S.) Cancer Center Support (Grant P30-CA14051)Howard Hughes Medical Institute. Predoctoral FellowshipCleo and Paul Schimmel Foundation. FellowshipNatural Sciences and Engineering Research Council of Canada PGS Scholarshi

The Dynamics of Supply and Demand in mRNA Translation

We study the elongation stage of mRNA translation in eukaryotes and find that, in contrast to the assumptions of previous models, both the supply and the demand for tRNA resources are important for determining elongation rates. We find that increasing the initiation rate of translation can lead to the depletion of some species of aa-tRNA, which in turn can lead to slow codons and queueing. Particularly striking “competition” effects are observed in simulations of multiple species of mRNA which are reliant on the same pool of tRNA resources. These simulations are based on a recent model of elongation which we use to study the translation of mRNA sequences from the Saccharomyces cerevisiae genome. This model includes the dynamics of the use and recharging of amino acid tRNA complexes, and we show via Monte Carlo simulation that this has a dramatic effect on the protein production behaviour of the system