Synthesizing and tuning chemical reaction networks with specified behaviours

We consider how to generate chemical reaction networks (CRNs) from functional specifications. We propose a two-stage approach that combines synthesis by satisfiability modulo theories and Markov chain Monte Carlo based optimisation. First, we identify candidate CRNs that have the possibility to produce correct computations for a given finite set of inputs. We then optimise the reaction rates of each CRN using a combination of stochastic search techniques applied to the chemical master equation, simultaneously improving the of correct behaviour and ruling out spurious solutions. In addition, we use techniques from continuous time Markov chain theory to study the expected termination time for each CRN. We illustrate our approach by identifying CRNs for majority decision-making and division computation, which includes the identification of both known and unknown networks.Comment: 17 pages, 6 figures, appeared the proceedings of the 21st conference on DNA Computing and Molecular Programming, 201

Mechanics and dynamics of X-chromosome pairing at X inactivation

At the onset of X-chromosome inactivation, the vital process whereby female mammalian cells equalize X products with respect to males, the X chromosomes are colocalized along their Xic (X-inactivation center) regions. The mechanism inducing recognition and pairing of the X’s remains, though, elusive. Starting from recent discoveries on the molecular factors and on the DNA sequences (the so-called "pairing sites") involved, we dissect the mechanical basis of Xic colocalization by using a statistical physics model. We show that soluble DNA-specific binding molecules, such as those experimentally identified, can be indeed sufficient to induce the spontaneous colocalization of the homologous chromosomes but only when their concentration, or chemical affinity, rises above a threshold value as a consequence of a thermodynamic phase transition. We derive the likelihood of pairing and its probability distribution. Chromosome dynamics has two stages: an initial independent Brownian diffusion followed, after a characteristic time scale, by recognition and pairing. Finally, we investigate the effects of DNA deletion/insertions in the region of pairing sites and compare model predictions to available experimental data

Arabidopsis plants perform arithmetic division to prevent starvation at night

Photosynthetic starch reserves that accumulate in Arabidopsis leaves during the day decrease approximately linearly with time at night to support metabolism and growth. We find that the rate of decrease is adjusted to accommodate variation in the time of onset of darkness and starch content, such that reserves last almost precisely until dawn. Generation of these dynamics therefore requires an arithmetic division computation between the starch content and expected time to dawn. We introduce two novel chemical kinetic models capable of implementing analog arithmetic division. Predictions from the models are successfully tested in plants perturbed by a night-time light period or by mutations in starch degradation pathways. Our experiments indicate which components of the starch degradation apparatus may be important for appropriate arithmetic division. Our results are potentially relevant for any biological system dependent on a food reserve for survival over a predictable time period

Normalizing single-cell RNA sequencing data: challenges and opportunities

Single-cell transcriptomics is becoming an important component of the molecular biologist's toolkit. A critical step when analyzing data generated using this technology is normalization. However, normalization is typically performed using methods developed for bulk RNA sequencing or even microarray data, and the suitability of these methods for single-cell transcriptomics has not been assessed. We here discuss commonly used normalization approaches and illustrate how these can produce misleading results. Finally, we present alternative approaches and provide recommendations for single-cell RNA sequencing users

A pivotal role for starch in the reconfiguration of 14C-partitioning and allocation in Arabidopsis thaliana under short-term abiotic stress.

Plant carbon status is optimized for normal growth but is affected by abiotic stress. Here, we used 14C-labeling to provide the first holistic picture of carbon use changes during short-term osmotic, salinity, and cold stress in Arabidopsis thaliana. This could inform on the early mechanisms plants use to survive adverse environment, which is important for efficient agricultural production. We found that carbon allocation from source to sinks, and partitioning into major metabolite pools in the source leaf, sink leaves and roots showed both conserved and divergent responses to the stresses examined. Carbohydrates changed under all abiotic stresses applied; plants re-partitioned 14C to maintain sugar levels under stress, primarily by reducing 14C into the storage compounds in the source leaf, and decreasing 14C into the pools used for growth processes in the roots. Salinity and cold increased 14C-flux into protein, but as the stress progressed, protein degradation increased to produce amino acids, presumably for osmoprotection. Our work also emphasized that stress regulated the carbon channeled into starch, and its metabolic turnover. These stress-induced changes in starch metabolism and sugar export in the source were partly accompanied by transcriptional alteration in the T6P/SnRK1 regulatory pathway that are normally activated by carbon starvation

Active and poised promoter states drive folding of the extended HoxB locus in mouse embryonic stem cells

Gene expression states influence the three-dimensional conformation of the genome through poorly understood mechanisms. Here, we investigate the conformation of the murine HoxB locus, a gene-dense genomic region containing closely spaced genes with distinct activation states in mouse embryonic stem (ES) cells. To predict possible folding scenarios, we performed computer simulations of polymer models informed with different chromatin occupancy features, which define promoter activation states or CTCF binding sites. Single cell imaging of the locus folding was performed to test model predictions. While CTCF occupancy alone fails to predict the in vivo folding at genomic length scale of 10 kb, we found that homotypic interactions between active and Polycomb-repressed promoters co-occurring in the same DNA fibre fully explain the HoxB folding patterns imaged in single cells. We identify state-dependent promoter interactions as major drivers of chromatin folding in gene-dense regions

Conformation Regulation of the X Chromosome Inactivation Center: A Model

X-Chromosome Inactivation (XCI) is the process whereby one, randomly chosen X becomes transcriptionally silenced in female cells. XCI is governed by the Xic, a locus on the X encompassing an array of genes which interact with each other and with key molecular factors. The mechanism, though, establishing the fate of the X's, and the corresponding alternative modifications of the Xic architecture, is still mysterious. In this study, by use of computer simulations, we explore the scenario where chromatin conformations emerge from its interaction with diffusing molecular factors. Our aim is to understand the physical mechanisms whereby stable, non-random conformations are established on the Xic's, how complex architectural changes are reliably regulated, and how they lead to opposite structures on the two alleles. In particular, comparison against current experimental data indicates that a few key cis-regulatory regions orchestrate the organization of the Xic, and that two major molecular regulators are involved