Development of constrained fuzzy logic for modeling biological regulatory networks and predicting contextual therapeutic effects

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Biological Engineering, 2012.Cataloged from PDF version of thesis.Includes bibliographical references (p. 199-213).Upon exposure to environmental cues, protein modifications form a complex signaling network that dictates cellular response. In this thesis, we develop methods for using continuous logic-based models to aide our understanding of these signaling networks and facilitate data interpretation. We present a novel modeling framework called constrained fuzzy logic (cFL) that maintains a simple logic-based description of interactions with AND, OR, and NOT gates, but allows for intermediate species activities with mathematical functions relating input and output values (transfer functions). We first train a prior knowledge network (PKN) to data with cFL, which reveals what aspects of the dataset agree or disagree with prior knowledge. The cFL models are trained to a dataset describing signaling proteins in a hepatocellular carcinoma cell line after exposure to ligand cues in the presence or absence of small molecule inhibitors. We find that multiple models with differing topology and parameters explain the data equally well, and it is crucial to consider this non-identifiability during model training and subsequence analysis. Our trained models generate new biological understanding of network crosstalk as well as quantitative predictions of signaling protein activation. In our next applications of cFL, we explore the ability of models either constructed based solely on prior knowledge or trained to dedicated biochemical data to make predictions that answer the following questions: 1) What perturbations to species in the system are effective at accomplishing a clinical goal? and 2) In what environmental conditions are these perturbations effective? We find that we are able to make accurate predictions in both cases. Thus, we offer cFL as a flexible modeling methodology to assist data interpretation and hypothesis generation for choice of therapeutic targets.by Melody K. Morris.Ph.D

The Implications and Concerns of Algorithms in the Modern World

In an ever emerging world of technology and advancements, there should be an increase in the importance and awareness of what each user is allowing and agreeing to on the Internet. Algorithms are becoming a part of everyday life, and many individuals don’t know just how much of their data, personal information, and privacy is being violated and collected from different companies and applications that they use. There needs to be an increase in the awareness and acknowledgement of algorithms and the “F.E.A.T” principles that companies should be following. New technologies from companies, like Google and Nimo, that are considered exciting and interesting in the technology world bring along their own concerns over privacy and the data that is being collected on users. In the next twenty years, technology will change as everyone knows it. How do users keep up with the new technologies that are coming up, and what are the implications of these algorithms and technologies in the modern world

Systematic Analysis of Quantitative Logic Model Ensembles Predicts Drug Combination Effects on Cell Signaling Networks

A major challenge in developing anticancer therapies is determining the efficacies of drugs and their combinations in physiologically relevant microenvironments. We describe here our application of “constrained fuzzy logic” (CFL) ensemble modeling of the intracellular signaling network for predicting inhibitor treatments that reduce the phospho-levels of key transcription factors downstream of growth factors and inflammatory cytokines representative of hepatocellular carcinoma (HCC) microenvironments. We observed that the CFL models successfully predicted the effects of several kinase inhibitor combinations. Furthermore, the ensemble predictions revealed ambiguous predictions that could be traced to a specific structural feature of these models, which we resolved with dedicated experiments, finding that IL-1α activates downstream signals through TAK1 and not MEKK1 in HepG2 cells. We conclude that CFL-Q2LM (Querying Quantitative Logic Models) is a promising approach for predicting effective anticancer drug combinations in cancer-relevant microenvironments.United States. Army Research Office (W911NF-09-0001

Logic-Based Models for the Analysis of Cell Signaling Networks

Computational models are increasingly used to analyze the operation of complex biochemical networks, including those involved in cell signaling networks. Here we review recent advances in applying logic-based modeling to mammalian cell biology. Logic-based models represent biomolecular networks in a simple and intuitive manner without describing the detailed biochemistry of each interaction. A brief description of several logic-based modeling methods is followed by six case studies that demonstrate biological questions recently addressed using logic-based models and point to potential advances in model formalisms and training procedures that promise to enhance the utility of logic-based methods for studying the relationship between environmental inputs and phenotypic or signaling state outputs of complex signaling networks.National Institutes of Health (U.S.) (Grant P50- GM68762)National Institutes of Health (U.S.) (Grant U54-CA112967)United States. Dept. of Defense (Institute for Collaborative Biotechnologies

Comparing Signaling Networks between Normal and Transformed Hepatocytes Using Discrete Logical Models

Substantial effort in recent years has been devoted to constructing and analyzing large-scale gene and protein networks on the basis of “omic” data and literature mining. These interaction graphs provide valuable insight into the topologies of complex biological networks but are rarely context specific and cannot be used to predict the responses of cell signaling proteins to specific ligands or drugs. Conversely, traditional approaches to analyzing cell signaling are narrow in scope and cannot easily make use of network-level data. Here, we combine network analysis and functional experimentation by using a hybrid approach in which graphs are converted into simple mathematical models that can be trained against biochemical data. Specifically, we created Boolean logic models of immediate-early signaling in liver cells by training a literature-based prior knowledge network against biochemical data obtained from primary human hepatocytes and 4 hepatocellular carcinoma cell lines exposed to combinations of cytokines and small-molecule kinase inhibitors. Distinct families of models were recovered for each cell type, and these families clustered topologically into normal and diseased sets.National Institutes of Health (U.S.) (Grant GM68762)National Institutes of Health (U.S.) (Grant CA112967

Oncogenic D816V-KIT signaling in mast cells causes persistent IL-6 production

Persistent dysregulation of IL-6 production and signaling have been implicated in the pathology of various cancers. In systemic mastocytosis, increased serum levels of IL-6 associate with disease severity and progression, although the mechanisms involved are not well understood. Since systemic mastocytosis often associates with the presence in hematopoietic cells of a somatic gain-of-function variant in KIT, D816V-KIT, we examined its potential role in IL-6 upregulation. Bone marrow mononuclear cultures from patients with greater D816V allelic burden released increased amounts of IL-6 which correlated with the percentage of mast cells in the cultures. Intracellular IL-6 staining by flow cytometry and immunofluorescence was primarily associated with mast cells and suggested a higher percentage of IL-6 positive mast cells in patients with higher D816V allelic burden. Furthermore, mast cell lines expressing D816V-KIT, but not those expressing normal KIT or other KIT variants, produced constitutively high IL-6 amounts at the message and protein levels. We further demonstrate that aberrant KIT activity and signaling are critical for the induction of IL-6 and involve STAT5 and PI3K pathways but not STAT3 or STAT4. Activation of STAT5A and STATB downstream of D816V-KIT was mediated by JAK2 but also by MEK/ERK1/2, which not only promoted STAT5 phosphorylation but also its long-term transcription. Our study thus supports a role for mast cells and D816V-KIT activity in IL-6 dysregulation in mastocytosis and provides insights into the intracellular mechanisms. The findings contribute to a better understanding of the physiopathology of mastocytosis and suggest the importance of therapeutic targeting of these pathwaysThis work was supported by the Division of Intramural Research within the National Institute of Allergy and Infectious Diseases (NIAID), at the National Institutes of Health.S

Normalization and Statistical Analysis of Multiplexed Bead-Based Immunoassay Data Using Mixed-Effects Modeling

Multiplexed bead-based flow cytometric immunoassays are a powerful experimental tool for investigating cellular communication networks, yet their widespread adoption is limited in part by challenges in robust quantitative analysis of the measurements. Here we report our application of mixed-effects modeling for the normalization and statistical analysis of bead-based immunoassay data. Our data set consisted of bead-based immunoassay measurements of 16 phospho-proteins in lysates of HepG2 cells treated with ligands that regulate acute-phase protein secretion. Mixed-effects modeling provided estimates for the effects of both the technical and biological sources of variance, and normalization was achieved by subtracting the technical effects from the measured values. This approach allowed us to detect ligand effects on signaling with greater precision and sensitivity and to more accurately characterize the HepG2 cell signaling network using constrained fuzzy logic. Mixed-effects modeling analysis of our data was vital for ascertaining that IL-1α and TGF-α treatment increased the activities of more pathways than IL-6 and TNF-α and that TGF-α and TNF-α increased p38 MAPK and c-Jun N-terminal kinase (JNK) phospho-protein levels in a synergistic manner. Moreover, we used mixed-effects modeling-based technical effect estimates to reveal the substantial variance contributed by batch effects along with the absence of loading order and assay plate position effects. We conclude that mixed-effects modeling enabled additional insights to be gained from our data than would otherwise be possible and we discuss how this methodology can play an important role in enhancing the value of experiments employing multiplexed bead-based immunoassays.United States. Army Research Office (Contract W911NF-09-D-0001)National Institutes of Health (U.S.) (NIH P50-GM68762

Protein signatures of centenarians and their offspring suggest centenarians age slower than other humans

Using samples from the New England Centenarian Study (NECS), we sought to characterize the serum proteome of 77 centenarians, 82 centenarians\u27 offspring, and 65 age-matched controls of the offspring (mean ages: 105, 80, and 79 years). We identified 1312 proteins that significantly differ between centenarians and their offspring and controls (FDR \u3c 1%), and two different protein signatures that predict longer survival in centenarians and in younger people. By comparing the centenarian signature with 2 independent proteomic studies of aging, we replicated the association of 484 proteins of aging and we identified two serum protein signatures that are specific of extreme old age. The data suggest that centenarians acquire similar aging signatures as seen in younger cohorts that have short survival periods, suggesting that they do not escape normal aging markers, but rather acquire them much later than usual. For example, centenarian signatures are significantly enriched for senescence-associated secretory phenotypes, consistent with those seen with younger aged individuals, and from this finding, we provide a new list of serum proteins that can be used to measure cellular senescence. Protein co-expression network analysis suggests that a small number of biological drivers may regulate aging and extreme longevity, and that changes in gene regulation may be important to reach extreme old age. This centenarian study thus provides additional signatures that can be used to measure aging and provides specific circulating biomarkers of healthy aging and longevity, suggesting potential mechanisms that could help prolong health and support longevity

Convective and Wave Signatures in Ozone Profiles Over the Equatorial Americas: Views from TC4 (2007) and SHADOZ

During the months of July-August 2007 NASA conducted a research campaign called the Tropical Composition, Clouds and Climate Coupling (TC4) experiment. Vertical profiles of ozone were measured daily using an instrument known as an ozonesonde, which is attached to a weather balloon and launch to altitudes in excess of 30 km. These ozone profiles were measured over coastal Las Tablas, Panama (7.8N, 80W) and several times per week at Alajuela, Costa Rica (ION, 84W). Meteorological systems in the form of waves, detected most prominently in 100- 300 in thick ozone layer in the tropical tropopause layer, occurred in 50% (Las Tablas) and 40% (Alajuela) of the soundings. These layers, associated with vertical displacements and classified as gravity waves ("GW," possibly Kelvin waves), occur with similar stricture and frequency over the Paramaribo (5.8N, 55W) and San Cristobal (0.925, 90W) sites of the Southern Hemisphere Additional Ozonesondes (SHADOZ) network. The gravity wave labeled layers in individual soundings correspond to cloud outflow as indicated by the tracers measured from the NASA DC-8 and other aircraft data, confirming convective initiation of equatorial waves. Layers representing quasi-horizontal displacements, referred to as Rossby waves, are robust features in soundings from 23 July to 5 August. The features associated with Rossby waves correspond to extra-tropical influence, possibly stratospheric, and sometimes to pollution transport. Comparison of Las Tablas and Alajuela ozone budgets with 1999-2007 Paramaribo and San Cristobal soundings shows that TC4 is typical of climatology for the equatorial Americas. Overall during TC4, convection and associated meteorological waves appear to dominate ozone transport in the tropical tropopause layer