Jag1-Notch cis-interaction determines cell fate segregation in pancreatic development.

The Notch ligands Jag1 and Dll1 guide differentiation of multipotent pancreatic progenitor cells (MPCs) into unipotent pro-acinar cells (PACs) and bipotent duct/endocrine progenitors (BPs). Ligand-mediated trans-activation of Notch receptors induces oscillating expression of the transcription factor Hes1, while ligand-receptor cis-interaction indirectly represses Hes1 activation. Despite Dll1 and Jag1 both displaying cis- and trans-interactions, the two mutants have different phenotypes for reasons not fully understood. Here, we present a mathematical model that recapitulates the spatiotemporal differentiation of MPCs into PACs and BPs. The model correctly captures cell fate changes in Notch pathway knockout mice and small molecule inhibitor studies, and a requirement for oscillatory Hes1 expression to maintain the multipotent state. Crucially, the model entails cell-autonomous attenuation of Notch signaling by Jag1-mediated cis-inhibition in MPC differentiation. The model sheds light on the underlying mechanisms, suggesting that cis-interaction is crucial for exiting the multipotent state, while trans-interaction is required for adopting the bipotent fate.

A possible universal role for mRNA secondary structure in bacterial translation revealed using a synthetic operon.

In bacteria, translation re-initiation is crucial for synthesizing proteins encoded by genes that are organized into operons. The mechanisms regulating translation re-initiation remain, however, poorly understood. We now describe the ribosome termination structure (RTS), a conserved and stable mRNA secondary structure localized immediately downstream of stop codons, and provide experimental evidence for its role in governing re-initiation efficiency in a synthetic Escherichia coli operon. We further report that RTSs are abundant, being associated with 18%-65% of genes in 128 analyzed bacterial genomes representing all phyla, and are selectively depleted when translation re-initiation is advantageous yet selectively enriched so as to insulate translation when re-initiation is deleterious. Our results support a potentially universal role for the RTS in controlling translation termination-insulation and re-initiation across bacteria.

Expanding the Genetic Code of a Photoautotrophic Organism.

The photoautotrophic freshwater cyanobacterium Synechococcus elongatus is widely used as a chassis for biotechnological applications as well as a photosynthetic bacterial model. In this study, a method for expanding the genetic code of this cyanobacterium has been established, thereby allowing the incorporation of unnatural amino acids into proteins. This was achieved through UAG stop codon suppression, using an archaeal pyrrolysyl orthogonal translation system. We demonstrate incorporation of unnatural amino acids into green fluorescent protein with 20 ± 3.5% suppression efficiency. The introduced components were shown to be orthogonal to the host translational machinery. In addition, we observed that no significant growth impairment resulted from the integration of the system. To interpret the observations, we modeled and investigated the competition over the UAG codon between release factor 1 and pyl-tRNACUA. On the basis of the model results, and the fact that 39.6% of the stop codons in the S. elongatus genome are UAG stop codons, the suppression efficiency in S. elongatus is unexpectedly high. The reason for this unexpected suppression efficiency has yet to be determined.

Tuning of Recombinant Protein Expression in Escherichia coli by Manipulating Transcription, Translation Initiation Rates, and Incorporation of Noncanonical Amino Acids.

Protein synthesis in cells has been thoroughly investigated and characterized over the past 60 years. However, some fundamental issues remain unresolved, including the reasons for genetic code redundancy and codon bias. In this study, we changed the kinetics of the Eschrichia coli transcription and translation processes by mutating the promoter and ribosome binding domains and by using genetic code expansion. The results expose a counterintuitive phenomenon, whereby an increase in the initiation rates of transcription and translation lead to a decrease in protein expression. This effect can be rescued by introducing slow translating codons into the beginning of the gene, by shortening gene length or by reducing initiation rates. On the basis of the results, we developed a biophysical model, which suggests that the density of co-transcriptional-translation plays a role in bacterial protein synthesis. These findings indicate how cells use codon bias to tune translation speed and protein synthesis.

A Monte Carlo Study of the Early Steps of Functional Amyloid Formation.

In addition to their well-known roles in neurodegenerative diseases and amyloidoses, amyloid structures also assume important functional roles in the cell. Although functional amyloid shares many physiochemical properties with its pathogenic counterpart, it is evolutionarily optimized to avoid cytotoxicity. This makes it an interesting study case for aggregation phenomenon in general. One of the most well-known examples of a functional amyloid, E. coli curli, is an essential component in the formation of bacterial biofilm, and is primarily formed by aggregates of the protein CsgA. Previous studies have shown that the minor sequence variations observed in the five different subrepeats (R1-R5), which comprise the CsgA primary sequence, have a substantial influence on their individual aggregation propensities. Using a recently described diffusion-optimized enhanced sampling approach for Monte Carlo simulations, we here investigate the equilibrium properties of the monomeric and dimeric states of these subrepeats, to probe whether structural properties observed in these early stage oligomers are decisive for the characteristics of the resulting aggregate. We show that the dimerization propensities of these peptides have strong correlations with their propensity for amyloid formation, and provide structural insights into the inter- and intramolecular contacts that appear to be essential in this process.

Coupled positive and negative feedbacks produce diverse gene expression patterns in colonies.

UNLABELLED: Formation of patterns is a common feature in the development of multicellular organism as well as of microbial communities. To investigate the formation of gene expression patterns in colonies, we build a mathematical model of two-dimensional colony growth, where cells carry a coupled positive-and-negative-feedback circuit. We demonstrate that the model can produce sectored, target (concentric), uniform, and scattered expression patterns of regulators, depending on gene expression dynamics and nutrient diffusion. We reconstructed the same regulatory structure in Escherichia coli cells and found gene expression patterns on the surface of colonies similar to the ones produced by the computer simulations. By comparing computer simulations and experimental results, we observed that very simple rules of gene expression can yield a spectrum of well-defined patterns in a growing colony. Our results suggest that variations of the protein content among cells lead to a high level of heterogeneity in colonies. IMPORTANCE: Formation of patterns is a common feature in the development of microbial communities. In this work, we show that a simple genetic circuit composed of a positive-feedback loop and a negative-feedback loop can produce diverse expression patterns in colonies. We obtained similar sets of gene expression patterns in the simulations and in the experiments. Because the combination of positive feedback and negative feedback is common in intracellular molecular networks, our results suggest that the protein content of cells is highly diversified in colonies.

Multiple roles of heparin in the aggregation of p25α.

The 219-residue protein p25α stimulates the fibrillation of α-synuclein (αSN) in vitro and colocalizes with it in several α-synucleinopathies. Although p25α does not fibrillate by itself under native conditions in vitro, αSN-free p25α aggregates have also been observed in vivo in, for example, multiple system atrophy. To investigate which environmental conditions might trigger this aggregation, we investigated the effect of polyanionic biomolecules on p25α aggregation. Heparin, polyglutamate, arachidonic acid micelles, and RNA all induce p25α aggregation. More detailed studies using heparin as template for aggregation reveal that a minimum of 10-14 heparin monosaccharide units per heparin polymer are required. Bona fide fibrils are only formed at intermediate heparin concentrations, possibly because an excess of heparin binding sites blocks the inter-p25α contacts required for amyloid formation. Other polyanions also show an optimum for amyloid formation. Aggregation involves only modest structural changes according to both spectroscopic and proteolytic experiments. The aggregates do not seed aggregation of heparin-free p25α, suggesting that heparin is required in stoichiometric amounts to form organized structures. We are able to reproduce these observations in a model involving two levels of binding of p25α to heparin. We conclude that the modest structural changes that p25α undergoes can promote weak intermolecular contacts and that polyanions such as heparin play a central role in stabilizing these aggregates but in multiple ways, leading to different types of aggregates. This highlights the role of non-protein components in promoting protein aggregation in vivo.

Analyzing inflammatory response as excitable media.

The regulatory system of the transcription factor NF-κB plays a great role in many cell functions, including inflammatory response. Interestingly, the NF-κB system is known to up-regulate production of its own triggering signal-namely, inflammatory cytokines such as TNF, IL-1, and IL-6. In this paper we investigate a previously presented model of the NF-κB, which includes both spatial effects and the positive feedback from cytokines. The model exhibits the properties of an excitable medium and has the ability to propagate waves of high cytokine concentration. These waves represent an optimal way of sending an inflammatory signal through the tissue as they create a chemotactic signal able to recruit neutrophils to the site of infection. The simple model displays three qualitatively different states; low stimuli leads to no or very little response. Intermediate stimuli leads to reoccurring waves of high cytokine concentration. Finally, high stimuli leads to a sustained high cytokine concentration, a scenario which is toxic for the tissue cells and corresponds to chronic inflammation. Due to the few variables of the simple model, we are able to perform a phase-space analysis leading to a detailed understanding of the functional form of the model and its limitations. The spatial effects of the model contribute to the robustness of the cytokine wave formation and propagation.

Dickkopf1--a new player in modelling the Wnt pathway.

The Wnt signaling pathway transducing the stabilization of ß-catenin is essential for metazoan embryo development and is misregulated in many diseases such as cancers. In recent years models have been proposed for the Wnt signaling pathway during the segmentation process in developing embryos. Many of these include negative feedback loops where Axin2 plays a key role. However, Axin2 null mice show no segmentation phenotype. We therefore propose a new model where the negative feedback involves Dkk1 rather than Axin2. We show that this model can exhibit the same type of oscillations as the previous models with Axin2 and as observed in experiments. We show that a spatial Wnt gradient can consistently convert this temporal periodicity into the spatial periodicity of somites, provided the oscillations in new cells arising in the presomitic mesoderm are synchronized with the oscillations of older cells. We further investigate the hypothesis that a change in the Wnt level in the tail bud during the later stages of somitogenesis can lengthen the time period of the oscillations and hence the size and separation of the later somites.

Genetic regulation of fluxes: iron homeostasis of Escherichia coli.

Iron is an essential trace-element for most organisms. However, because high concentration of free intracellular iron is cytotoxic, cells have developed complex regulatory networks that keep free intracellular iron concentration at optimal range, allowing the incorporation of the metal into iron-using enzymes and minimizing damage to the cell. We built a mathematical model of the network that controls iron uptake and usage in the bacterium Escherichia coli to explore the dynamics of iron flow. We simulate the effect of sudden decrease or increase in the extracellular iron level on intracellular iron distribution. Based on the results of simulations we discuss the possible roles of the small RNA RyhB and the Fe-S cluster assembly systems in the optimal redistribution of iron flows. We suggest that Fe-S cluster assembly is crucial to prevent the accumulation of toxic levels of free intracellular iron when the environment suddenly becomes iron rich.