Modeling Stem Cell Induction Processes

Technology for converting human cells to pluripotent stem cell using induction processes has the potential to revolutionize regenerative medicine. However, the production of these so called iPS cells is still quite inefficient and may be dominated by stochastic effects. In this work we build mass-action models of the core regulatory elements controlling stem cell induction and maintenance. The models include not only the network of transcription factors NANOG, OCT4, SOX2, but also important epigenetic regulatory features of DNA methylation and histone modification. We show that the network topology reported in the literature is consistent with the observed experimental behavior of bistability and inducibility. Based on simulations of stem cell generation protocols, and in particular focusing on changes in epigenetic cellular states, we show that cooperative and independent reaction mechanisms have experimentally identifiable differences in the dynamics of reprogramming, and we analyze such differences and their biological basis. It had been argued that stochastic and elite models of stem cell generation represent distinct fundamental mechanisms. Work presented here suggests an alternative possibility that they represent differences in the amount of information we have about the distribution of cellular states before and during reprogramming protocols. We show further that unpredictability and variation in reprogramming decreases as the cell progresses along the induction process, and that identifiable groups of cells with elite-seeming behavior can come about by a stochastic process. Finally we show how different mechanisms and kinetic properties impact the prospects of improving the efficiency of iPS cell generation protocols.Fundação para a Ciência e a Tecnologia (BD 42942)MIT-Portugal ProgramNational Institutes of Health (U.S.) (CA112967)Singapore–MIT Alliance for Research and TechnologyIntel Corporatio

A single molecule assay to probe monovalent and multivalent bonds between hyaluronan and its key leukocyte receptor CD44 under force

Glycosaminoglycans (GAGs), a category of linear, anionic polysaccharides, are ubiquitous in the extracellular space, and important extrinsic regulators of cell function. Despite the recognized significance of mechanical stimuli in cellular communication, however, only few single molecule methods are currently available to study how monovalent and multivalent GAG•protein bonds respond to directed mechanical forces. Here, we have devised such a method, by combining purpose-designed surfaces that afford immobilization of GAGs and receptors at controlled nanoscale organizations with single molecule force spectroscopy (SMFS). We apply the method to study the interaction of the GAG polymer hyaluronan (HA) with CD44, its receptor in vascular endothelium. Individual bonds between HA and CD44 are remarkably resistant to rupture under force in comparison to their low binding affinity. Multiple bonds along a single HA chain rupture sequentially and independently under load. We also demonstrate how strong non-covalent bonds, which are versatile for controlled protein and GAG immobilization, can be effectively used as molecular anchors in SMFS. We thus establish a versatile method for analyzing the nanomechanics of GAG•protein interactions at the level of single GAG chains, which provides new molecular-level insight into the role of mechanical forces in the assembly and function of GAG-rich extracellular matrices

The route to transcription initiation determines the mode of transcriptional bursting in E. coli

Transcription is fundamentally noisy, leading to significant heterogeneity across bacterial populations. Noise is often attributed to burstiness, but the underlying mechanisms and their dependence on the mode of promotor regulation remain unclear. Here, we measure E. coli single cell mRNA levels for two stress responses that depend on bacterial sigma factors with different mode of transcription initiation (σ70 and σ54). By fitting a stochastic model to the observed mRNA distributions, we show that the transition from low to high expression of the σ70-controlled stress response is regulated via the burst size, while that of the σ54-controlled stress response is regulated via the burst frequency. Therefore, transcription initiation involving σ54 differs from other bacterial systems, and yields bursting kinetics characteristic of eukaryotic systems

Phenotypic Signatures Arising from Unbalanced Bacterial Growth

Fluctuations in the growth rate of a bacterial culture during unbalanced growth are generally considered undesirable in quantitative studies of bacterial physiology. Under well-controlled experimental conditions, however, these fluctuations are not random but instead reflect the interplay between intra-cellular networks underlying bacterial growth and the growth environment. Therefore, these fluctuations could be considered quantitative phenotypes of the bacteria under a specific growth condition. Here, we present a method to identify “phenotypic signatures” by time-frequency analysis of unbalanced growth curves measured with high temporal resolution. The signatures are then applied to differentiate amongst different bacterial strains or the same strain under different growth conditions, and to identify the essential architecture of the gene network underlying the observed growth dynamics. Our method has implications for both basic understanding of bacterial physiology and for the classification of bacterial strains

Insights on the Neuromagnetic Representation of Temporal Asymmetry in Human Auditory Cortex.

Communication sounds are typically asymmetric in time and human listeners are highly sensitive to this short-term temporal asymmetry. Nevertheless, causal neurophysiological correlates of auditory perceptual asymmetry remain largely elusive to our current analyses and models. Auditory modelling and animal electrophysiological recordings suggest that perceptual asymmetry results from the presence of multiple time scales of temporal integration, central to the auditory periphery. To test this hypothesis we recorded auditory evoked fields (AEF) elicited by asymmetric sounds in humans. We found a strong correlation between perceived tonal salience of ramped and damped sinusoids and the AEFs, as quantified by the amplitude of the N100m dynamics. The N100m amplitude increased with stimulus half-life time, showing a maximum difference between the ramped and damped stimulus for a modulation half-life time of 4 ms which is greatly reduced at 0.5 ms and 32 ms. This behaviour of the N100m closely parallels psychophysical data in a manner that: i) longer half-life times are associated with a stronger tonal percept, and ii) perceptual differences between damped and ramped are maximal at 4 ms half-life time. Interestingly, differences in evoked fields were significantly stronger in the right hemisphere, indicating some degree of hemispheric specialisation. Furthermore, the N100m magnitude was successfully explained by a pitch perception model using multiple scales of temporal integration of auditory nerve activity patterns. This striking correlation between AEFs, perception, and model predictions suggests that the physiological mechanisms involved in the processing of pitch evoked by temporal asymmetric sounds are reflected in the N100m

Effects of Noise Bandwidth and Amplitude Modulation on Masking in Frog Auditory Midbrain Neurons

Natural auditory scenes such as frog choruses consist of multiple sound sources (i.e., individual vocalizing males) producing sounds that overlap extensively in time and spectrum, often in the presence of other biotic and abiotic background noise. Detection of a signal in such environments is challenging, but it is facilitated when the noise shares common amplitude modulations across a wide frequency range, due to a phenomenon called comodulation masking release (CMR). Here, we examined how properties of the background noise, such as its bandwidth and amplitude modulation, influence the detection threshold of a target sound (pulsed amplitude modulated tones) by single neurons in the frog auditory midbrain. We found that for both modulated and unmodulated masking noise, masking was generally stronger with increasing bandwidth, but it was weakened for the widest bandwidths. Masking was less for modulated noise than for unmodulated noise for all bandwidths. However, responses were heterogeneous, and only for a subpopulation of neurons the detection of the probe was facilitated when the bandwidth of the modulated masker was increased beyond a certain bandwidth – such neurons might contribute to CMR. We observed evidence that suggests that the dips in the noise amplitude are exploited by TS neurons, and observed strong responses to target signals occurring during such dips. However, the interactions between the probe and masker responses were nonlinear, and other mechanisms, e.g., selective suppression of the response to the noise, may also be involved in the masking release

In vivo single cell analysis reveals Gata2 dynamics in cells transitioning to hematopoietic fate

Cell fate is established through coordinated gene expression programs in individual cells. Regulatory networks that include the Gata2 transcription factor play central roles in hematopoietic fate establishment. Although Gata2 is essential to the embryonic development and function of hematopoietic stem cells that form the adult hierarchy, little is known about the in vivo expression dynamics of Gata2 in single cells. Here, we examine Gata2 expression in single aortic cells as they establish hematopoietic fate in Gata2Venus mouse embryos. Time-lapse imaging reveals rapid pulsatile level changes in Gata2 reporter expression in cells undergoing endothelial-to-hematopoietic transition. Moreover, Gata2 reporter pulsatile expression is dramatically altered in Gata2+/- aortic cells, which undergo fewer transitions and are reduced in hematopoietic potential. Our novel finding of dynamic pulsatile expression of Gata2 suggests a highly unstable genetic state in single cells concomitant with their transition to hematopoietic fate. This reinforces the notion that threshold levels of Gata2 influence fate establishment and has implications for transcription factor-related hematologic dysfunctions

Inference for stochastic chemical kinetics using moment equations and system size expansion

Quantitative mechanistic models are valuable tools for disentangling biochemical pathways and for achieving a comprehensive understanding of biological systems. However, to be quantitative the parameters of these models have to be estimated from experimental data. In the presence of significant stochastic fluctuations this is a challenging task as stochastic simulations are usually too time-consuming and a macroscopic description using reaction rate equations (RREs) is no longer accurate. In this manuscript, we therefore consider moment-closure approximation (MA) and the system size expansion (SSE), which approximate the statistical moments of stochastic processes and tend to be more precise than macroscopic descriptions. We introduce gradient-based parameter optimization methods and uncertainty analysis methods for MA and SSE. Efficiency and reliability of the methods are assessed using simulation examples as well as by an application to data for Epo-induced JAK/STAT signaling. The application revealed that even if merely population-average data are available, MA and SSE improve parameter identifiability in comparison to RRE. Furthermore, the simulation examples revealed that the resulting estimates are more reliable for an intermediate volume regime. In this regime the estimation error is reduced and we propose methods to determine the regime boundaries. These results illustrate that inference using MA and SSE is feasible and possesses a high sensitivity