Single cell classification of macrophage subtypes by label-free cell signatures and machine learning

Pro-inflammatory (M1) and anti-inflammatory (M2) macrophage phenotypes play a fundamental role in the immune response. The interplay and consequently the classification between these two functional subtypes is significant for many therapeutic applications. Albeit, a fast classification of macrophage phenotypes is challenging. For instance, image-based classification systems need cell staining and coloration, which is usually time- and cost-consuming, such as multiple cell surface markers, transcription factors and cytokine profiles are needed. A simple alternative would be to identify such cell types by using single-cell, label-free and high throughput light scattering pattern analyses combined with a straightforward machine learning-based classification. Here, we compared different machine learning algorithms to classify distinct macrophage phenotypes based on their optical signature obtained from an ad hoc developed wide-angle static light scattering apparatus. As the main result, we were able to identify unpolarized macrophages from M1- and M2-polarized phenotypes and distinguished them from naive monocytes with an average accuracy above 85%. Therefore, we suggest that optical single-cell signatures within a lab-on-a-chip approach along with machine learning could be used as a fast, affordable, non-invasive macrophage phenotyping tool to supersede resource-intensive cell labelling.

Temporal responsiveness of adipose-derived stem/stromal cell immune plasticity

We determined the role of time in adipose-derived stem/stromal cell (ASC) response to a model inflammatory environment. ASCs and other mesenchymal stem/stromal cells exhibit immune plasticity. We evaluated the persistence of pro- and anti-inflammatory phenotypes for ASCs exposed to a sustained or pulse inflammatory stimulus. Using qPCR, flow cytometry, and immunocytochemistry, we monitored the temporal expression and up-regulation patterns of a pro-inflammatory gene (caspase 1), a pleiotropic gene/protein (interleukin 6, IL-6), and an anti-inflammatory gene/protein (indoleamine 2, 3-dioxygenase, IDO1) after exposing ASCs to the cytokines tumor necrosis factor-α and interferon-γ. In response to sustained cytokine stimulation, we discovered that time played a role in the balance of pro- and anti-inflammatory ASC phenotypes. IL-6 was present at all time points for both cytokine-stimulated and non-stimulated conditions, whereas IDO1 was heterogeneously up-regulated in stimulated conditions at later time points. After a pulse stimulus, ASC immunoresponse remained consistent for 96-168 h. As a final measure of immune plasticity, we cultured cytokine-stimulated ASCs with blood-derived macrophages to observe macrophage polarization. While the presence of ASCs altered macrophage phenotype, there was no dependency on the length of ASC cytokine exposure time.

Molecular mechanisms of detection and discrimination of dynamic signals

Many molecules decode not only the concentration of cellular signals, but also their temporal dynamics. However, little is known about the mechanisms that underlie the detection and discrimination of dynamic signals. We used computational modelling of the interaction of a ligand with multiple targets to investigate how kinetic and thermodynamic parameters regulate their capabilities to respond to dynamic signals. Our results demonstrated that the detection and discrimination of temporal features of signal inputs occur for reactions proceeding outside mass-action equilibrium. For these reactions, thermodynamic parameters such as affinity do not predict their outcomes. Additionally, we showed that, at non-equilibrium, the association rate constants determine the amount of product formed in reversible reactions. In contrast, the dissociation rate constants regulate the time interval required for reversible reactions to achieve equilibrium and, consequently, control their ability to detect and discriminate dynamic features of cellular signals.

A data-driven modeling approach to identify disease-specific multi-organ networks driving physiological dysregulation

Multiple physiological systems interact throughout the development of a complex disease. Knowledge of the dynamics and connectivity of interactions across physiological systems could facilitate the prevention or mitigation of organ damage underlying complex diseases, many of which are currently refractory to available therapeutics (e.g., hypertension). We studied the regulatory interactions operating within and across organs throughout disease development by integrating in vivo analysis of gene expression dynamics with a reverse engineering approach to infer data-driven dynamic network models of multi-organ gene regulatory influences. We obtained experimental data on the expression of 22 genes across five organs, over a time span that encompassed the development of autonomic nervous system dysfunction and hypertension. We pursued a unique approach for identification of continuous-time models that jointly described the dynamics and structure of multi-organ networks by estimating a sparse subset of ∼12,000 possible gene regulatory interactions. Our analyses revealed that an autonomic dysfunction-specific multi-organ sequence of gene expression activation patterns was associated with a distinct gene regulatory network. We analyzed the model structures for adaptation motifs, and identified disease-specific network motifs involving genes that exhibited aberrant temporal dynamics. Bioinformatic analyses identified disease-specific single nucleotide variants within or near transcription factor binding sites upstream of key genes implicated in maintaining physiological homeostasis. Our approach illustrates a novel framework for investigating the pathogenesis through model-based analysis of multi-organ system dynamics and network properties. Our results yielded novel candidate molecular targets driving the development of cardiovascular disease, metabolic syndrome, and immune dysfunction.

Complex diseases such as hypertension often involve maladaptive autonomic nervous system control over the cardiovascular, renal, hepatic, immune, and endocrine systems. We studied the pathogenesis of physiological homeostasis by examining the temporal dynamics of gene expression levels from multiple organs in an animal model of autonomic dysfunction characterized by cardiovascular disease, metabolic dysregulation, and immune system aberrations. We employed a data-driven modeling approach to jointly predict continuous gene expression dynamics and gene regulatory interactions across organs in the disease and control phenotypes. We combined our analyses of multi-organ gene regulatory network dynamics and connectivity with bioinformatic analyses of genetic mutations that could regulate gene expression. Our multi-organ modeling approach to investigate the mechanisms of complex disease pathogenesis revealed novel candidates for therapeutic interventions against the development and progression of complex diseases involving autonomic nervous system dysfunction.

Acetylation of MKL1 by PCAF regulates pro-inflammatory transcription

Inflammation is considered a fundamental host defense mechanism and, when aberrantly activated, contributes to a host of human diseases. Previously we have reported that the transcriptional regulator megakaryocytic leukemia 1 (MKL1) plays a role programming cellular inflammatory response by modulating NF-κB activity. Here we report that MKL1 was acetylated in vivo and pro-inflammatory stimuli (TNF-α and LPS) augmented MKL1 acetylation accompanying increased MKL1 binding to NF-κB target promoters. Further analysis revealed that the lysine acetyltransferase PCAF mediated MKL1 acetylation: TNF-α and LPS promoted the interaction between MKL1 and PCAF whereas depletion of PCAF abrogated the induction of MKL1 acetylation by TNF-α and LPS. Acetylation of MKL1 was necessary for MKL1 to activate the transcription of pro-inflammatory genes because mutation of four conserved lysine residues in MKL1 attenuated its capacity as a trans-activator of NF-κB target genes. Mechanistically, MKL1 acetylation served to promote MKL1 nuclear enrichment, to enhance the MKL1-NF-κB interaction, and to stabilize the binding of MKL1 on target promoters. In conclusion, our data unveil an important pathway that contributes to the transcriptional regulation of inflammatory response.

Understanding Physiology in the Continuum: Integration of Information from Multiple -Omics Levels

In this paper, we discuss approaches for integrating biological information reflecting diverse physiologic levels. In particular, we explore statistical and model-based methods for integrating transcriptomic, proteomic and metabolomics data. Our case studies reflect responses to a systemic inflammatory stimulus and in response to an anti-inflammatory treatment. Our paper serves partly as a review of existing methods and partly as a means to demonstrate, using case studies related to human endotoxemia and response to methylprednisolone (MPL) treatment, how specific questions may require specific methods, thus emphasizing the non-uniqueness of the approaches. Finally, we explore novel ways for integrating -omics information with PKPD models, toward the development of more integrated pharmacology models.

Quantitative Systems Pharmacology: A Framework for Context

Quantitative Systems Pharmacology (QSP) is receiving increased attention. As the momentum builds and the expectations grow it is important to (re)assess and formalize the basic concepts and approaches. In this short review, I argue that QSP, in addition to enabling the rational integration of data and development of complex models, maybe more importantly, provides the foundations for developing an integrated framework for the assessment of drugs and their impact on disease within a broader context expanding the envelope to account in great detail for physiology, environment and prior history. I articulate some of the critical enablers, major obstacles and exciting opportunities manifesting themselves along the way. Charting such overarching themes will enable practitioners to identify major and defining factors as the field progressively moves towards personalized and precision health care delivery.

Self-Catalyzed Assembly of Peptide Scaffolded Nanozyme as a Dynamic Biosensing System

In this work, a new strategy of biosensor design is developed based on the assembly of amyloid beta and its multiple interactions with other bioactive species. These interactions can enable amyloid beta peptide as a multifunctional sensing element, so the immobilization of sensing probe and the step-by-step modification of the sensing interface have all been dispensed with. Instead, the kinetics of the assembly of a peptide-based catalytic network serves to convert the quantity of analyte into amplified signal readout. The designed dynamic assembling and biosensing system has also been successfully applied in detecting the activity of polyglutamylation, an essential post translation modification controlling cell skeleton and cell cycle, in biological complex samples. Further studies reveal that the serum abundance of a polyglutamylase, tubulin tyrosine ligase-like protein 12, may show parallel with the degree of development of prostate cancer and the discrimination between early cancerous development and benign conditions. And the obtained result is more distinct than that based on PSA detection, the current gold standard. This study may also point to the prospective of extending this design strategy to broader range of biosensing applications in the future.

Quantitative analysis of receptor tyrosine kinase-effector coupling at functionally relevant stimulus levels

A major goal of current signaling research is to develop a quantitative understanding of how receptor activation is coupled to downstream signaling events and to functional cellular responses. Here, we measure how activation of the RET receptor tyrosine kinase on mouse neuroblastoma cells by the neurotrophin artemin (ART) is quantitatively coupled to key downstream effectors. We show that the efficiency of RET coupling to ERK and Akt depends strongly on ART concentration, and it is highest at the low (∼100 pM) ART levels required for neurite outgrowth. Quantitative discrimination between ERK and Akt pathway signaling similarly is highest at this low ART concentration. Stimulation of the cells with 100 pM ART activated RET at the rate of ∼10 molecules/cell/min, leading at 5-10 min to a transient peak of ∼150 phospho-ERK (pERK) molecules and ∼50 pAkt molecules per pRET, after which time the levels of these two signaling effectors fell by 25-50% while the pRET levels continued to slowly rise. Kinetic experiments showed that signaling effectors in different pathways respond to RET activation with different lag times, such that the balance of signal flux among the different pathways evolves over time. Our results illustrate that measurements using high, super-physiological growth factor levels can be misleading about quantitative features of receptor signaling. We propose a quantitative model describing how receptor-effector coupling efficiency links signal amplification to signal sensitization between receptor and effector, thereby providing insight into design principles underlying how receptors and their associated signaling machinery decode an extracellular signal to trigger a functional cellular outcome.

Dynamic multiprotein assemblies shape the spatial structure of cell signaling

Cell signaling underlies critical cellular decisions. Coordination, efficiency as well as fail-safe mechanisms are key elements. How the cell ensures that these hallmarks are at play are important questions. Cell signaling is often viewed as taking place through discrete and cross-talking pathways; oftentimes these are modularized to emphasize distinct functions. While simple, convenient and clear, such models largely neglect the spatial structure of cell signaling; they also convey inter-modular (or inter-protein) spatial separation that may not exist. Here our thesis is that cell signaling is shaped by a network of multiprotein assemblies. While pre-organized, the assemblies and network are loose and dynamic. They contain transiently-associated multiprotein complexes which are often mediated by scaffolding proteins. They are also typically anchored in the membrane, and their continuum may span the cell. IQGAP1 scaffolding protein which binds proteins including Raf, calmodulin, Mek, Erk, actin, and tens more, with actin shaping B-cell (and likely other) membrane-anchored nanoclusters and allosterically polymerizing in dynamic cytoskeleton formation, and Raf anchoring in the membrane along with Ras, provides a striking example. The multivalent network of dynamic proteins and lipids, with specific interactions forming and breaking, can be viewed as endowing gel-like properties. Collectively, this reasons that efficient, productive and reliable cell signaling takes place primarily through transient, preorganized and cooperative protein-protein interactions spanning the cell rather than stochastic, diffusion-controlled processes.