Ligand-specific c-Fos expression emerges from the spatiotemporal control of ErbB network dynamics.

Activation of ErbB receptors by epidermal growth factor (EGF) or heregulin (HRG) determines distinct cell-fate decisions, although signals propagate through shared pathways. Using mathematical modeling and experimental approaches, we unravel how HRG and EGF generate distinct, all-or-none responses of the phosphorylated transcription factor c-Fos. In the cytosol, EGF induces transient and HRG induces sustained ERK activation. In the nucleus, however, ERK activity and c-fos mRNA expression are transient for both ligands. Knockdown of dual-specificity phosphatases extends HRG-stimulated nuclear ERK activation, but not c-fos mRNA expression, implying the existence of a HRG-induced repressor of c-fos transcription. Further experiments confirmed that this repressor is mainly induced by HRG, but not EGF, and requires new protein synthesis. We show how a spatially distributed, signaling-transcription cascade robustly discriminates between transient and sustained ERK activities at the c-Fos system level. The proposed control mechanisms are general and operate in different cell types, stimulated by various ligands.

Increasing Signal Intensity of Fluorescent Oligo-Labeled Antibodies to Enable Combination Multiplexing.

Full-spectrum flow cytometry has increased antibody-based multiplexing, yet further increases remain potentially impactful. We recently proposed how fluorescence multiplexing using spectral imaging and combinatorics (MuSIC) could do so using tandem dyes and an oligo-based antibody labeling method. In this work, we found that such labeled antibodies had significantly lower signal intensities than conventionally labeled antibodies in human cell experiments. To improve signal intensity, we tested moving the fluorophores from the original external (ext.) 5' or 3' end-labeled orientation to internal (int.) fluorophore modifications. Cell-free spectrophotometer measurements showed a ∼6-fold signal intensity increase of the new int. configuration compared to the previous ext. configuration. Time-resolved fluorescence and fluorescence correlation spectroscopy showed that the ∼3-fold brightness difference is due to static quenching most likely by the oligo or solution in the ext. configuration. Spectral flow cytometry experiments using peripheral blood mononuclear cells show int. MuSIC probe-labeled antibodies (i) retained increased signal intensity while having no significant difference in the estimated % of CD8+ lymphocytes and (ii) labeled with Atto488, Atto647, and Atto488/647 combinations can be demultiplexed in triple-stained samples. The antibody labeling approach is general and can be broadly applied to many biological and diagnostic applications where spectral detection is available.

Predicting anti-cancer drug combination responses with a temporal cell state network model.

Cancer chemotherapy combines multiple drugs, but predicting the effects of drug combinations on cancer cell proliferation remains challenging, even for simple in vitro systems. We hypothesized that by combining knowledge of single drug dose responses and cell state transition network dynamics, we could predict how a population of cancer cells will respond to drug combinations. We tested this hypothesis here using three targeted inhibitors of different cell cycle states in two different cell lines in vitro. We formulated a Markov model to capture temporal cell state transitions between different cell cycle phases, with single drug data constraining how drug doses affect transition rates. This model was able to predict the landscape of all three different pairwise drug combinations across all dose ranges for both cell lines with no additional data. While further application to different cell lines, more drugs, additional cell state networks, and more complex co-culture or in vivo systems remain, this work demonstrates how currently available or attainable information could be sufficient for prediction of drug combination response for single cell lines in vitro.

Protein structure-based gene expression signatures.

Gene expression signatures (GES) connect phenotypes to differential messenger RNA (mRNA) expression of genes, providing a powerful approach to define cellular identity, function, and the effects of perturbations. The use of GES has suffered from vague assessment criteria and limited reproducibility. Because the structure of proteins defines the functional capability of genes, we hypothesized that enrichment of structural features could be a generalizable representation of gene sets. We derive structural gene expression signatures (sGES) using features from multiple levels of protein structure (e.g., domain and fold) encoded by the mRNAs in GES. Comprehensive analyses of data from the Genotype-Tissue Expression Project (GTEx), the all RNA-seq and ChIP-seq sample and signature search (ARCHS4) database, and mRNA expression of drug effects on cardiomyocytes show that sGES are useful for characterizing biological phenomena. sGES enable phenotypic characterization across experimental platforms, facilitates interoperability of expression datasets, and describe drug action on cells.

SysMod: the ISCB community for data-driven computational modelling and multi-scale analysis of biological systems.

Computational models of biological systems can exploit a broad range of rapidly developing approaches, including novel experimental approaches, bioinformatics data analysis, emerging modelling paradigms, data standards and algorithms. A discussion about the most recent advances among experts from various domains is crucial to foster data-driven computational modelling and its growing use in assessing and predicting the behaviour of biological systems. Intending to encourage the development of tools, approaches and predictive models, and to deepen our understanding of biological systems, the Community of Special Interest (COSI) was launched in Computational Modelling of Biological Systems (SysMod) in 2016. SysMod's main activity is an annual meeting at the Intelligent Systems for Molecular Biology (ISMB) conference, which brings together computer scientists, biologists, mathematicians, engineers, computational and systems biologists. In the five years since its inception, SysMod has evolved into a dynamic and expanding community, as the increasing number of contributions and participants illustrate. SysMod maintains several online resources to facilitate interaction among the community members, including an online forum, a calendar of relevant meetings and a YouTube channel with talks and lectures of interest for the modelling community. For more than half a decade, the growing interest in computational systems modelling and multi-scale data integration has inspired and supported the SysMod community. Its members get progressively more involved and actively contribute to the annual COSI meeting and several related community workshops and meetings, focusing on specific topics, including particular techniques for computational modelling or standardisation efforts.

Method for Improved Fluorescence Corrections for Molar Mass Characterization by Multiangle Light Scattering.

Multiangle light scattering (MALS) was used to determine the absolute molar mass of fluorescent macromolecules. It is standard protocol to install bandwidth filters before MALS detectors to suppress detection of fluorescent emissions. Fluorescence can introduce tremendous error in light scattering measurements and is a formidable challenge in accurately characterizing fluorescent macromolecules and particles. However, we show that for some systems, bandwidth filters alone are insufficient for blocking fluorescence in molar mass determinations. For these systems, we have devised a correction procedure to calculate the amount of fluorescence interference in the filtered signal. By determining the intensity of fluorescent emission not blocked by the bandwidth filters, we can correct the filtered signal accordingly and accurately determine the true molar mass. The transmission rates are calculated before MALS experimentation using emission data from standard fluorimetry techniques, allowing for the characterization of unknown samples. To validate the correction procedure, we synthesized fluorescent dye-conjugated proteins using an IR800CW (LI-COR) fluorophore and Bovine Serum Albumin protein. We successfully eliminated fluorescence interference in MALS measurements using this approach. This correction procedure has potential application toward more accurate molar mass characterizations of macromolecules with intrinsic fluorescence, such as lignins, fluorescent proteins, fluorescence-tagged proteins, and optically active nanoparticles.

Protocol for Creating Antibodies with Complex Fluorescence Spectra.

Fluorescent antibodies are a workhorse of biomedical science, but fluorescence multiplexing has been notoriously difficult due to spectral overlap between fluorophores. We recently established proof-of-principal for fluorescence Multiplexing using Spectral Imaging and Combinatorics (MuSIC), which uses combinations of existing fluorophores to create unique spectral signatures for increased multiplexing. However, a method for labeling antibodies with MuSIC probes has not yet been developed. Here, we present a method for labeling antibodies with MuSIC probes. We conjugate a DBCO-Peg5-NHS ester linker to antibodies and a single-stranded DNA "docking strand" to the linker and, finally, hybridize two MuSIC-compatible, fluorescently labeled oligos to the docking strand. We validate the labeling protocol with spin-column purification and absorbance measurements. We demonstrate the approach using (i) Cy3, (ii) Tex615, and (iii) a Cy3-Tex615 combination as three different MuSIC probes attached to three separate batches of antibodies. We created single-, double-, and triple-positive beads that are analogous to single cells by incubating MuSIC probe-labeled antibodies with protein A beads. Spectral flow cytometry experiments demonstrate that each MuSIC probe can be uniquely distinguished, and the fraction of beads in a mixture with different staining patterns are accurately inferred. The approach is general and might be more broadly applied to cell-type profiling or tissue heterogeneity studies in clinical, biomedical, and drug discovery research.

A mechanistic pan-cancer pathway model informed by multi-omics data interprets stochastic cell fate responses to drugs and mitogens.

Most cancer cells harbor multiple drivers whose epistasis and interactions with expression context clouds drug and drug combination sensitivity prediction. We constructed a mechanistic computational model that is context-tailored by omics data to capture regulation of stochastic proliferation and death by pan-cancer driver pathways. Simulations and experiments explore how the coordinated dynamics of RAF/MEK/ERK and PI-3K/AKT kinase activities in response to synergistic mitogen or drug combinations control cell fate in a specific cellular context. In this MCF10A cell context, simulations suggest that synergistic ERK and AKT inhibitor-induced death is likely mediated by BIM rather than BAD, which is supported by prior experimental studies. AKT dynamics explain S-phase entry synergy between EGF and insulin, but simulations suggest that stochastic ERK, and not AKT, dynamics seem to drive cell-to-cell proliferation variability, which in simulations is predictable from pre-stimulus fluctuations in C-Raf/B-Raf levels. Simulations suggest MEK alteration negligibly influences transformation, consistent with clinical data. Tailoring the model to an alternate cell expression and mutation context, a glioma cell line, allows prediction of increased sensitivity of cell death to AKT inhibition. Our model mechanistically interprets context-specific landscapes between driver pathways and cell fates, providing a framework for designing more rational cancer combination therapy.

Cell size assays for mass cytometry.

Mass cytometry offers the advantage of allowing the simultaneous measurement of a greater number parameters than conventional flow cytometry. However, to date, mass cytometry has lacked a reliable alternative to the light scatter properties that are commonly used as a cell size metric in flow cytometry (forward scatter intensity-FSC). Here, we report the development of two plasma membrane staining assays to evaluate mammalian cell size in mass cytometry experiments. One is based on wheat germ agglutinin (WGA) staining and the other on Osmium tetroxide (OsO4 ) staining, both of which have preferential affinity for cell membranes. We first perform imaging and flow cytometry experiments to establish a relationship between WGA staining intensity and traditional measures of cell size. We then incorporate WGA staining in mass cytometry analysis of human whole blood and show that WGA staining intensity has reproducible patterns within and across immune cell subsets that have distinct cell sizes. Lastly, we stain PBMCs or dissociated lung tissue with both WGA and OsO4 ; mass cytometry analysis demonstrates that the two staining intensities correlate well with one another. We conclude that both WGA and OsO4 may be used to acquire cell size-related parameters in mass cytometry experiments, and expect these stains to be broadly useful in expanding the range of parameters that can be measured in mass cytometry experiments. © 2016 International Society for Advancement of Cytometry.

ERK2 drives tumour cell migration in three-dimensional microenvironments by suppressing expression of Rab17 and liprin-ß2.

Upregulation of the extracellular signal-regulated kinase (ERK) pathway has been shown to contribute to tumour invasion and progression. Because the two predominant ERK isoforms (ERK1 and ERK2, also known as MAPK3 and MAPK1, respectively) are highly homologous and have indistinguishable kinase activities in vitro, both enzymes were believed to be redundant and interchangeable. To challenge this view, we show that ERK2 silencing inhibits invasive migration of MDA-MB-231 cells, and re-expression of ERK2 but not ERK1 restores the normal invasive phenotype. A detailed quantitative analysis of cell movement on 3D matrices indicates that ERK2 knockdown impairs cellular motility by decreasing the migration velocity as well as increasing the time that cells spend not moving. Using gene expression arrays we found that the expression of the genes for Rab17 and liprin-ß2 was increased by knockdown of ERK2 and restored to normal levels following re-expression of ERK2, but not ERK1. Both play inhibitory roles in the invasive behaviour of three independent cancer cell lines. Importantly, knockdown of either Rab17 or liprin-ß2 restores invasiveness of ERK2-depleted cells, indicating that ERK2 drives invasion of MDA-MB-231 cells by suppressing expression of these genes.

Network Analyses of Brain Tumor Patients' Multiomic Data Reveals Pharmacological Opportunities to Alter Cell State Transitions.

Glioblastoma Multiforme (GBM) remains a particularly difficult cancer to treat, and survival outcomes remain poor. In addition to the lack of dedicated drug discovery programs for GBM, extensive intratumor heterogeneity and epigenetic plasticity related to cell-state transitions are major roadblocks to successful drug therapy in GBM. To study these phenomenon, publicly available snRNAseq and bulk RNAseq data from patient samples were used to categorize cells from patients into four cell states (i.e. phenotypes), namely: (i) neural progenitor-like (NPC-like), (ii) oligodendrocyte progenitor-like (OPC-like), (iii) astrocyte-like (AC-like), and (iv) mesenchymal-like (MES-like). Patients were subsequently grouped into subpopulations based on which cell-state was the most dominant in their respective tumor. By incorporating phosphoproteomic measurements from the same patients, a protein-protein interaction network (PPIN) was constructed for each cell state. These four-cell state PPINs were pooled to form a single Boolean network that was used for in silico protein knockout simulations to investigate mechanisms that either promote or prevent cell state transitions. Simulation results were input into a boosted tree machine learning model which predicted the cell states or phenotypes of GBM patients from an independent public data source, the Glioma Longitudinal Analysis (GLASS) Consortium. Combining the simulation results and the machine learning predictions, we generated hypotheses for clinically relevant causal mechanisms of cell state transitions. For example, the transcription factor TFAP2A can be seen to promote a transition from the NPC-like to the MES-like state. Such protein nodes and the associated signaling pathways provide potential drug targets that can be further tested in vitro and support cell state-directed (CSD) therapy.

Network Analyses of Brain Tumor Patients' Multiomic Data Reveals Pharmacological Opportunities to Alter Cell State Transitions.

Glioblastoma Multiforme (GBM) remains a particularly difficult cancer to treat, and survival outcomes remain poor. In addition to the lack of dedicated drug discovery programs for GBM, extensive intratumor heterogeneity and epigenetic plasticity related to cell-state transitions are major roadblocks to successful drug therapy in GBM. To study these phenomenon, publicly available snRNAseq and bulk RNAseq data from patient samples were used to categorize cells from patients into four cell states (i.e. phenotypes), namely: (i) neural progenitor-like (NPC-like), (ii) oligodendrocyte progenitor-like (OPC-like), (iii) astrocyte- like (AC-like), and (iv) mesenchymal-like (MES-like). Patients were subsequently grouped into subpopulations based on which cell-state was the most dominant in their respective tumor. By incorporating phosphoproteomic measurements from the same patients, a protein-protein interaction network (PPIN) was constructed for each cell state. These four-cell state PPINs were pooled to form a single Boolean network that was used for in silico protein knockout simulations to investigate mechanisms that either promote or prevent cell state transitions. Simulation results were input into a boosted tree machine learning model which predicted the cell states or phenotypes of GBM patients from an independent public data source, the Glioma Longitudinal Analysis (GLASS) Consortium. Combining the simulation results and the machine learning predictions, we generated hypotheses for clinically relevant causal mechanisms of cell state transitions. For example, the transcription factor TFAP2A can be seen to promote a transition from the NPC-like to the MES-like state. Such protein nodes and the associated signaling pathways provide potential drug targets that can be further tested in vitro and support cell state-directed (CSD) therapy.

Low-frequency ERK and Akt activity dynamics are predictive of stochastic cell division events.

Understanding the dynamics of intracellular signaling pathways, such as ERK1/2 (ERK) and Akt1/2 (Akt), in the context of cell fate decisions is important for advancing our knowledge of cellular processes and diseases, particularly cancer. While previous studies have established associations between ERK and Akt activities and proliferative cell fate, the heterogeneity of single-cell responses adds complexity to this understanding. This study employed a data-driven approach to address this challenge, developing machine learning models trained on a dataset of growth factor-induced ERK and Akt activity time courses in single cells, to predict cell division events. The most predictive models were developed by applying discrete wavelet transforms (DWTs) to extract low-frequency features from the time courses, followed by using Ensemble Integration, a data integration and predictive modeling framework. The results demonstrated that these models effectively predicted cell division events in MCF10A cells (F-measure=0.524, AUC=0.726). ERK dynamics were found to be more predictive than Akt, but the combination of both measurements further enhanced predictive performance. The ERK model`s performance also generalized to predicting division events in RPE cells, indicating the potential applicability of these models and our data-driven methodology for predicting cell division across different biological contexts. Interpretation of these models suggested that ERK dynamics throughout the cell cycle, rather than immediately after growth factor stimulation, were associated with the likelihood of cell division. Overall, this work contributes insights into the predictive power of intra-cellular signaling dynamics for cell fate decisions, and highlights the potential of machine learning approaches in unraveling complex cellular behaviors.

Mechanistic modeling of cell viability assays with in silico lineage tracing.

Data from cell viability assays, which measure cumulative division and death events in a population and reflect substantial cellular heterogeneity, are widely available. However, interpreting such data with mechanistic computational models is hindered because direct model/data comparison is often muddled. We developed an algorithm that tracks simulated division and death events in mechanistically detailed single-cell lineages to enable such a model/data comparison and suggest causes of cell-cell drug response variability. Using our previously developed model of mammalian single-cell proliferation and death signaling, we simulated drug dose response experiments for four targeted anti-cancer drugs (alpelisib, neratinib, trametinib and palbociclib) and compared them to experimental data. Simulations are consistent with data for strong growth inhibition by trametinib (MEK inhibitor) and overall lack of efficacy for alpelisib (PI-3K inhibitor), but are inconsistent with data for palbociclib (CDK4/6 inhibitor) and neratinib (EGFR inhibitor). Model/data inconsistencies suggest (i) the importance of CDK4/6 for driving the cell cycle may be overestimated, and (ii) that the cellular balance between basal (tonic) and ligand-induced signaling is a critical determinant of receptor inhibitor response. Simulations show subpopulations of rapidly and slowly dividing cells in both control and drug-treated conditions. Variations in mother cells prior to drug treatment all impinging on ERK pathway activity are associated with the rapidly dividing phenotype and trametinib resistance. This work lays a foundation for the application of mechanistic modeling to large-scale cell viability assay datasets and better understanding determinants of cellular heterogeneity in drug response.

A 96-Well Polyacrylamide Gel for Electrophoresis and Western Blotting.

Western blotting is a stalwart technique for analyzing specific proteins and/or their post-translational modifications. However, it remains challenging to accommodate more than ~10 samples per experiment without substantial departure from trusted, established protocols involving accessible instrumentation. Here, we describe a 96-sample western blot that conforms to standard 96-well plate dimensional constraints and has little operational deviation from standard western blotting. The main differences are that (i) submerged polyacrylamide gel electrophoresis is operated horizontally (similar to agarose gels) as opposed to vertically, and (ii) a 6 mm thick gel is used, with 2 mm most relevant for membrane transfer (vs ~1 mm typical). Results demonstrate both wet and semi-dry transfer are compatible with this gel thickness. The major tradeoff is reduced molecular weight resolution, due primarily to less available migration distance per sample. We demonstrate proof-of-principle using gels loaded with molecular weight ladder, recombinant protein, and cell lysates. We expect the 96-well western blot will increase reproducibility, efficiency (cost and time ~8-fold), and capacity for biological characterization relative to established western blots.

Multiscale mapping of transcriptomic signatures for cardiotoxic drugs.

Drug-induced gene expression profiles can identify potential mechanisms of toxicity. We focus on obtaining signatures for cardiotoxicity of FDA-approved tyrosine kinase inhibitors (TKIs) in human induced-pluripotent-stem-cell-derived cardiomyocytes, using bulk transcriptomic profiles. We use singular value decomposition to identify drug-selective patterns across cell lines obtained from multiple healthy human subjects. Cellular pathways affected by cardiotoxic TKIs include energy metabolism, contractile, and extracellular matrix dynamics. Projecting these pathways to published single cell expression profiles indicates that TKI responses can be evoked in both cardiomyocytes and fibroblasts. Integration of transcriptomic outlier analysis with whole genomic sequencing of our six cell lines enables us to correctly reidentify a genomic variant causally linked to anthracycline-induced cardiotoxicity and predict genomic variants potentially associated with TKI-induced cardiotoxicity. We conclude that mRNA expression profiles when integrated with publicly available genomic, pathway, and single cell transcriptomic datasets, provide multiscale signatures for cardiotoxicity that could be used for drug development and patient stratification.

Bayesian multivariate Poisson abundance models for T-cell receptor data.

A major feature of an adaptive immune system is its ability to generate B- and T-cell clones capable of recognizing and neutralizing specific antigens. These clones recognize antigens with the help of the surface molecules, called antigen receptors, acquired individually during the clonal development process. In order to ensure a response to a broad range of antigens, the number of different receptor molecules is extremely large, resulting in a huge clonal diversity of both B- and T-cell receptor populations and making their experimental comparisons statistically challenging. To facilitate such comparisons, we propose a flexible parametric model of multivariate count data and illustrate its use in a simultaneous analysis of multiple antigen receptor populations derived from mammalian T-cells. The model relies on a representation of the observed receptor counts as a multivariate Poisson abundance mixture (m PAM). A Bayesian parameter fitting procedure is proposed, based on the complete posterior likelihood, rather than the conditional one used typically in similar settings. The new procedure is shown to be considerably more efficient than its conditional counterpart (as measured by the Fisher information) in the regions of m PAM parameter space relevant to model T-cell data.

NFκB and HIF display synergistic behaviour during hypoxic inflammation.

The oxygen-sensitive transcription factor hypoxia inducible factor (HIF) is a key regulator of gene expression during adaptation to hypoxia. Crucially, inflamed tissue often displays regions of prominent hypoxia. Recent studies have shown HIF signalling is intricately linked to that of the pro-inflammatory transcription factor nuclear factor kappa B (NFκB) during hypoxic inflammation. We describe the relative temporal contributions of each to hypoxia-induced inflammatory gene expression and investigate the level of crosstalk between the two pathways using a novel Gaussia princeps luciferase (Gluc) reporter system. Under the control of an active promoter, Gluc is expressed and secreted into the cell culture media, where it can be sampled and measured over time. Thus, Gluc constructs under the control of either HIF or NFκB were used to resolve their temporal transcriptional dynamics in response to hypoxia and to cytokine stimuli, respectively. We also investigated the interactions between HIF and NFκB activities using a construct containing the sequence from the promoter of the inflammatory gene cyclooxygenase 2 (COX-2), which includes functionally active binding sites for both HIF and NFκB. Finally, based on our experimental data, we constructed a mathematical model of the binding affinities of HIF and NFκB to their respective response elements to analyse transcriptional crosstalk. Taken together, these data reveal distinct temporal HIF and NFκB transcriptional activities in response to hypoxic inflammation. Furthermore, we demonstrate synergistic activity between these two transcription factors on the regulation of the COX-2 promoter, implicating a co-ordinated role for both HIF and NFκB in the expression of COX-2 in hypoxic inflammation.

Modeling the Dynamics of Eukaryotic DNA Synthesis in Remembrance of Tunde Ogunnaike.

Two things Tunde loved were dynamics and probability. The work described herein combined them both, which explains why Tunde invariably asked me each time we talked how this work was proceeding. However, as I've come to appreciate as reminiscent of a surprisingly large amount of work in almost any researcher's career, I did not complete a peer-reviewed article on the matter while he could see it. We were broadly motivated by analysis of data for total DNA content in single cells, across thousands of cells. From such data one can estimate the proportions of cells in different phases of the cell cycle by fitting a mixture model for subpopulations of G0/G1 phase cells (1 relative copy of the genome), S phase cells (between 1 and 2 relative copies of the genome), and G2/M phase cells (2 relative copies of the genome). Given an asynchronously cycling population, Gaussian models are reasonable for the G0/G1 and G2/M subpopulations, but an appropriate functional form for the S-phase subpopulation was unclear. Since the probability of observing an S-phase cell is intimately related to the dynamics of DNA replication, we worked to derive a model for DNA replication dynamics from first principles, resulting in a closed-form, analytic expression for the dynamics of DNA synthesis. While quite arguably a somewhat superfluous effort, there is a certain satisfaction and academic beauty to modeling systems from a first-principles approach, and it can sometimes lead to unexpected scientific insights. Yet, while mathematically elegant, there was a fundamental issue with a key assumption that the so-called inter-origin distance distribution (distances between DNA replication initiation sites) was time-invariant. First, I present the model as developed previously. Then, to address the time-invariant inter-origin distance distribution issue, I provide a treatment of time-varying inter-origin distance distributions that, while mathematically simple, provides (i) mechanistic predictions for how all the DNA in a fertilized frog egg can be replicated "on time" despite some inter-origin distances initially exceeding the corresponding amount of allowable time and (ii) evidence that, based only on data from DNA content versus time and average inter-origin distances, somatic cell DNA is parsed into distinct regions whose replication is temporally separate.