Shibin Mathew

Analysis of alternative signaling pathways of endoderm induction of human embryonic stem cells identifies context specific differences

BackgroundLineage specific differentiation of human embryonic stem cells (hESCs) is largely mediated by specific growth factors and extracellular matrix molecules. Growth factors initiate a cascade of signals which control gene transcription and cell fate specification. There is a lot of interest in inducing hESCs to an endoderm fate which serves as a pathway towards more functional cell types like the pancreatic cells. Research over the past decade has established several robust pathways for deriving endoderm from hESCs, with the capability of further maturation. However, in our experience, the functional maturity of these endoderm derivatives, specifically to pancreatic lineage, largely depends on specific pathway of endoderm induction. Hence it will be of interest to understand the underlying mechanism mediating such induction and how it is translated to further maturation. In this work we analyze the regulatory interactions mediating different pathways of endoderm induction by identifying co-regulated transcription factors.ResultshESCs were induced towards endoderm using activin A and 4 different growth factors (FGF2 (F), BMP4 (B), PI3KI (P), and WNT3A (W)) and their combinations thereof, resulting in 15 total experimental conditions. At the end of differentiation each condition was analyzed by qRT-PCR for 12 relevant endoderm related transcription factors (TFs). As a first approach, we used hierarchical clustering to identify which growth factor combinations favor up-regulation of different genes. In the next step we identified sets of co-regulated transcription factors using a biclustering algorithm. The high variability of experimental data was addressed by integrating the biclustering formulation with bootstrap re-sampling to identify robust networks of co-regulated transcription factors. Our results show that the transition from early to late endoderm is favored by FGF2 as well as WNT3A treatments under high activin. However, induction of late endoderm markers is relatively favored by WNT3A under high activin.ConclusionsUse of FGF2, WNT3A or PI3K inhibition with high activin A may serve well in definitive endoderm induction followed by WNT3A specific signaling to direct the definitive endoderm into late endodermal lineages. Other combinations, though still feasible for endoderm induction, appear less promising for pancreatic endoderm specification in our experiments.

Potential for Pancreatic Maturation of Differentiating Human Embryonic Stem Cells Is Sensitive to the Specific Pathway of Definitive Endoderm Commitment

This study provides a detailed experimental and mathematical analysis of the impact of the initial pathway of definitive endoderm (DE) induction on later stages of pancreatic maturation. Human embryonic stem cells (hESCs) were induced to insulin-producing cells following a directed-differentiation approach. DE was induced following four alternative pathway modulations. DE derivatives obtained from these alternate pathways were subjected to pancreatic progenitor (PP) induction and maturation and analyzed at each stage. Results indicate that late stage maturation is influenced by the initial pathway of DE commitment. Detailed quantitative analysis revealed WNT3A and FGF2 induced DE cells showed highest expression of insulin, are closely aligned in gene expression patterning and have a closer resemblance to pancreatic organogenesis. Conversely, BMP4 at DE induction gave most divergent differentiation dynamics with lowest insulin upregulation, but highest glucagon upregulation. Additionally, we have concluded that early analysis of PP markers is indicative of its potential for pancreatic maturation.

Regulatory interactions maintaining self-renewal of human embryonic stem cells as revealed through a systems analysis of PI3K/AKT pathway.

MOTIVATION Maintenance of the self-renewal state in human embryonic stem cells (hESCs) is the foremost critical step for regenerative therapy applications. The insulin-mediated PI3K/AKT pathway is well appreciated as being the central pathway supporting hESC self-renewal; however, the regulatory interactions in the pathway that maintain cell state are not yet known. Identification of these regulatory pathway components will be critical for designing targeted interventions to facilitate a completely defined platform for hESC propagation and differentiation. Here, we have developed a systems analysis approach to identify regulatory components that control PI3K/AKT pathway in self-renewing hESCs. RESULTS A detailed mathematical model was adopted to explain the complex regulatory interactions in the PI3K/AKT pathway. We evaluated globally sensitive processes of the pathway in a computationally efficient manner by replacing the detailed model by a surrogate meta-model. Our mathematical analysis, supported by experimental validation, reveals that negative regulators of the molecules IRS1 and PIP3 primarily govern the steady state of the pathway in hESCs. Among the regulators, negative feedback via IRS1 reduces the sensitivity of various reactions associated with direct trunk of the pathway and also constraints the propagation of parameter uncertainty to the levels of post receptor signaling molecules. Furthermore, our results suggest that inhibition of negative feedback can significantly increase p-AKT levels and thereby, better support hESC self-renewal. Our integrated mathematical modeling and experimental workflow demonstrates the significant advantage of computationally efficient meta-model approaches to detect sensitive targets from signaling pathways. AVAILABILITY AND IMPLEMENTATION FORTRAN codes for the PI3K/AKT pathway and the RS-HDMR implementation are available from the authors upon request.

A Neutrophil Phenotype Model for Extracorporeal Treatment of Sepsis

Neutrophils play a central role in eliminating bacterial pathogens, but may also contribute to end-organ damage in sepsis. Interleukin-8 (IL-8), a key modulator of neutrophil function, signals through neutrophil specific surface receptors CXCR-1 and CXCR-2. In this study a mechanistic computational model was used to evaluate and deploy an extracorporeal sepsis treatment which modulates CXCR-1/2 levels. First, a simplified mechanistic computational model of IL-8 mediated activation of CXCR-1/2 receptors was developed, containing 16 ODEs and 43 parameters. Receptor level dynamics and systemic parameters were coupled with multiple neutrophil phenotypes to generate dynamic populations of activated neutrophils which reduce pathogen load, and/or primed neutrophils which cause adverse tissue damage when misdirected. The mathematical model was calibrated using experimental data from baboons administered a two-hour infusion of E coli and followed for a maximum of 28 days. Ensembles of parameters were generated using a Bayesian parallel tempering approach to produce model fits that could recreate experimental outcomes. Stepwise logistic regression identified seven model parameters as key determinants of mortality. Sensitivity analysis showed that parameters controlling the level of killer cell neutrophils affected the overall systemic damage of individuals. To evaluate rescue strategies and provide probabilistic predictions of their impact on mortality, time of onset, duration, and capture efficacy of an extracorporeal device that modulated neutrophil phenotype were explored. Our findings suggest that interventions aiming to modulate phenotypic composition are time sensitive. When introduced between 3–6 hours of infection for a 72 hour duration, the survivor population increased from 31% to 40–80%. Treatment efficacy quickly diminishes if not introduced within 15 hours of infection. Significant harm is possible with treatment durations ranging from 5–24 hours, which may reduce survival to 13%. In severe sepsis, an extracorporeal treatment which modulates CXCR-1/2 levels has therapeutic potential, but also potential for harm. Further development of the computational model will help guide optimal device development and determine which patient populations should be targeted by treatment.

Development of an Alginate Array Platform to Decouple the Effect of Multiparametric Perturbations on Human Pluripotent Stem Cells During Pancreatic Differentiation

Human embryonic stem cells (hESC)-derived functional cells hold great promise for regenerative cell therapy. Currently approved strategies for clinical translation requires the isolation of the hESCs-derived cells in materials allowing transfer of reagents but preventing integration with the host. However, hESC fate is known to be sensitive to its local microenvironment, both chemical and physical. Given the complexity of hESC response to environmental parameters, it will be important to evaluate the cell response to multiple combinatorial perturbations. Such complex perturbations are best enabled by exploiting high-throughput screening platforms. In this study, the authors report the effect of multivariate perturbations on hESC differentiation, enabled by the development of high throughput 3D alginate array platform. Specifically, the sensitivity of hESC propagation and pancreatic differentiation to substrate properties and cell culture configuration is analyzed. Cellular response to array perturbations is analyzed by quantitative imaging, and cell sensitivity was determined through statistical modeling. The results indicate that configuration is the stronger determinant of hESC proliferation and differentiation, while substrate properties fine-tune the expression around the average levels. This platform allowed for multiparametric perturbations, and in combination with statistical modeling, allows to identify the sensitivity of hESC proliferation and fate to multiparametric modulation.

Finding the Optimal Tradeoffs

Computational analyses of a half-million circuit topologies provide a rationale for why certain fold-change detection topologies are more prevalent in nature.