Signal execution modes emerge in biochemical reaction networks calibrated to experimental data.

Mathematical models of biomolecular networks are commonly used to study cellular processes; however, their usefulness to explain and predict dynamic behaviors is often questioned due to the unclear relationship between parameter uncertainty and network dynamics. In this work, we introduce PyDyNo (Python dynamic analysis of biochemical networks), a non-equilibrium reaction-flux based analysis to identify dominant reaction paths within a biochemical reaction network calibrated to experimental data. We first show, in a simplified apoptosis execution model, that despite the thousands of parameter vectors with equally good fits to experimental data, our framework identifies the dynamic differences between these parameter sets and outputs three dominant execution modes, which exhibit varying sensitivity to perturbations. We then apply our methodology to JAK2/STAT5 network in colony-forming unit-erythroid (CFU-E) cells and provide previously unrecognized mechanistic explanation for the survival responses of CFU-E cell population that would have been impossible to deduce with traditional protein-concentration based analyses.

Mapping the proteogenomic landscape enables prediction of drug response in acute myeloid leukemia.

Acute myeloid leukemia is a poor-prognosis cancer commonly stratified by genetic aberrations, but these mutations are often heterogeneous and fail to consistently predict therapeutic response. Here, we combine transcriptomic, proteomic, and phosphoproteomic datasets with ex vivo drug sensitivity data to help understand the underlying pathophysiology of AML beyond mutations. We measure the proteome and phosphoproteome of 210 patients and combine them with genomic and transcriptomic measurements to identify four proteogenomic subtypes that complement existing genetic subtypes. We build a predictor to classify samples into subtypes and map them to a "landscape" that identifies specific drug response patterns. We then build a drug response prediction model to identify drugs that target distinct subtypes and validate our findings on cell lines representing various stages of quizartinib resistance. Our results show how multiomics data together with drug sensitivity data can inform therapy stratification and drug combinations in AML.

CDK4/6 inhibition enhances SHP2 inhibitor efficacy and is dependent upon RB function in malignant peripheral nerve sheath tumors.

Malignant peripheral nerve sheath tumors (MPNSTs) are highly aggressive soft tissue sarcomas with limited treatment options, and new effective therapeutic strategies are desperately needed. We observe antiproliferative potency of genetic depletion of PTPN11 or pharmacological inhibition using the SHP2 inhibitor (SHP2i) TNO155. Our studies into the signaling response to SHP2i reveal that resistance to TNO155 is partially mediated by reduced RB function, and we therefore test the addition of a CDK4/6 inhibitor (CDK4/6i) to enhance RB activity and improve TNO155 efficacy. In combination, TNO155 attenuates the adaptive response to CDK4/6i, potentiates its antiproliferative effects, and converges on enhancement of RB activity, with greater suppression of cell cycle and inhibitor-of-apoptosis proteins, leading to deeper and more durable antitumor activity in in vitro and in vivo patient-derived models of MPNST, relative to either single agent. Overall, our study provides timely evidence to support the clinical advancement of this combination strategy in patients with MPNST and other tumors driven by loss of NF1.

Proteome Mapping of the Human Pancreatic Islet Microenvironment Reveals Endocrine-Exocrine Signaling Sphere of Influence.

The need for a clinically accessible method with the ability to match protein activity within heterogeneous tissues is currently unmet by existing technologies. Our proteomics sample preparation platform, named microPOTS (Microdroplet Processing in One pot for Trace Samples), can be used to measure relative protein abundance in micron-scale samples alongside the spatial location of each measurement, thereby tying biologically interesting proteins and pathways to distinct regions. However, given the smaller pixel/voxel number and amount of tissue measured, standard mass spectrometric analysis pipelines have proven inadequate. Here we describe how existing computational approaches can be adapted to focus on the specific biological questions asked in spatial proteomics experiments. We apply this approach to present an unbiased characterization of the human islet microenvironment comprising the entire complex array of cell types involved while maintaining spatial information and the degree of the islet's sphere of influence. We identify specific functional activity unique to the pancreatic islet cells and demonstrate how far their signature can be detected in the adjacent tissue. Our results show that we can distinguish pancreatic islet cells from the neighboring exocrine tissue environment, recapitulate known biological functions of islet cells, and identify a spatial gradient in the expression of RNA processing proteins within the islet microenvironment.

Ex vivo to in vivo model of malignant peripheral nerve sheath tumors for precision oncology.

BACKGROUND: Malignant peripheral nerve sheath tumors (MPNST) are aggressive soft tissue sarcomas that often develop in patients with neurofibromatosis type 1 (NF1). To address the critical need for novel therapeutics in MPNST, we aimed to establish an ex vivo 3D platform that accurately captured the genomic diversity of MPNST and could be utilized in a medium-throughput manner for drug screening studies to be validated in vivo using patient-derived xenografts (PDX). METHODS: Genomic analysis was performed on all PDX-tumor pairs. Selected PDX were harvested for assembly into 3D microtissues. Based on prior work in our labs, we evaluated drugs (trabectedin, olaparib, and mirdametinib) ex vivo and in vivo. For 3D microtissue studies, cell viability was the endpoint as assessed by Zeiss Axio Observer. For PDX drug studies, tumor volume was measured twice weekly. Bulk RNA sequencing was performed to identify pathways enriched in cells. RESULTS: We developed 13 NF1-associated MPNST-PDX and identified mutations or structural abnormalities in NF1 (100%), SUZ12 (85%), EED (15%), TP53 (15%), CDKN2A (85%), and chromosome 8 gain (77%). We successfully assembled PDX into 3D microtissues, categorized as robust (>90% viability at 48 h), good (>50%), or unusable (<50%). We evaluated drug response to "robust" or "good" microtissues, namely MN-2, JH-2-002, JH-2-079-c, and WU-225. Drug response ex vivo predicted drug response in vivo, and enhanced drug effects were observed in select models. CONCLUSIONS: These data support the successful establishment of a novel 3D platform for drug discovery and MPNST biology exploration in a system representative of the human condition.

Coupling microdroplet-based sample preparation, multiplexed isobaric labeling, and nanoflow peptide fractionation for deep proteome profiling of tissue microenvironment.

There is increasing interest in developing in-depth proteomic approaches for mapping tissue heterogeneity at a cell-type-specific level to better understand and predict the function of complex biological systems, such as human organs. Existing spatially resolved proteomics technologies cannot provide deep proteome coverages due to limited sensitivity and poor sample recovery. Herein, we seamlessly combined laser capture microdissection with a low-volume sample processing technology that includes a microfluidic device named microPOTS (Microdroplet Processing in One pot for Trace Samples), the multiplexed isobaric labelling, and a nanoflow peptide fractionation approach. The integrated workflow allowed to maximize proteome coverage of laser-isolated tissue samples containing nanogram proteins. We demonstrated the deep spatial proteomics can quantify more than 5,000 unique proteins from a small-sized human pancreatic tissue pixel (∼60,000 µm2) and reveal unique islet microenvironments.

CDK4/6 inhibition enhances SHP2 inhibitor efficacy and is dependent upon restoration of RB function in malignant peripheral nerve sheath tumors.

Malignant peripheral nerve sheath tumors (MPNST) are highly aggressive soft tissue sarcomas with limited treatment options, and novel effective therapeutic strategies are desperately needed. We observe anti-proliferative efficacy of genetic depletion or pharmacological inhibition using the clinically available SHP2 inhibitor (SHP2i) TNO155. Our studies into the signaling response to SHP2i reveal that resistance to TNO155 is partially mediated by reduced RB function, and we therefore test the addition of a CDK4/6 inhibitor (CDK4/6i) to enhance RB activity and improve TNO155 efficacy. In combination, TNO155 attenuates the adaptive response to CDK4/6i, potentiates its anti-proliferative effects, and converges on enhancement of RB activity, with greater suppression of cell cycle and inhibitor-of-apoptosis proteins, leading to deeper and more durable anti-tumor activity in in vitro and in vivo patient-derived models of MPNST, relative to either single agent. Overall, our study provides timely evidence to support the clinical advancement of this combination strategy in patients with MPNST and other tumors driven by loss of NF1.

Processes in DNA damage response from a whole-cell multi-omics perspective.

Technological advances have made it feasible to collect multi-condition multi-omic time courses of cellular response to perturbation, but the complexity of these datasets impedes discovery due to challenges in data management, analysis, visualization, and interpretation. Here, we report a whole-cell mechanistic analysis of HL-60 cellular response to bendamustine. We integrate both enrichment and network analysis to show the progression of DNA damage and programmed cell death over time in molecular, pathway, and process-level detail using an interactive analysis framework for multi-omics data. Our framework, Mechanism of Action Generator Involving Network analysis (MAGINE), automates network construction and enrichment analysis across multiple samples and platforms, which can be integrated into our annotated gene-set network to combine the strengths of networks and ontology-driven analysis. Taken together, our work demonstrates how multi-omics integration can be used to explore signaling processes at various resolutions and demonstrates multi-pathway involvement beyond the canonical bendamustine mechanism.

Proteomic and phosphoproteomic measurements enhance ability to predict ex vivo drug response in AML.

Acute Myeloid Leukemia (AML) affects 20,000 patients in the US annually with a five-year survival rate of approximately 25%. One reason for the low survival rate is the high prevalence of clonal evolution that gives rise to heterogeneous sub-populations of leukemic cells with diverse mutation spectra, which eventually leads to disease relapse. This genetic heterogeneity drives the activation of complex signaling pathways that is reflected at the protein level. This diversity makes it difficult to treat AML with targeted therapy, requiring custom patient treatment protocols tailored to each individual's leukemia. Toward this end, the Beat AML research program prospectively collected genomic and transcriptomic data from over 1000 AML patients and carried out ex vivo drug sensitivity assays to identify genomic signatures that could predict patient-specific drug responses. However, there are inherent weaknesses in using only genetic and transcriptomic measurements as surrogates of drug response, particularly the absence of direct information about phosphorylation-mediated signal transduction. As a member of the Clinical Proteomic Tumor Analysis Consortium, we have extended the molecular characterization of this cohort by collecting proteomic and phosphoproteomic measurements from a subset of these patient samples (38 in total) to evaluate the hypothesis that proteomic signatures can improve the ability to predict response to 26 drugs in AML ex vivo samples. In this work we describe our systematic, multi-omic approach to evaluate proteomic signatures of drug response and compare protein levels to other markers of drug response such as mutational patterns. We explore the nuances of this approach using two drugs that target key pathways activated in AML: quizartinib (FLT3) and trametinib (Ras/MEK), and show how patient-derived signatures can be interpreted biologically and validated in cell lines. In conclusion, this pilot study demonstrates strong promise for proteomics-based patient stratification to assess drug sensitivity in AML.

Selective Inhibition of JAK1 Primes STAT5-Driven Human Leukemia Cells for ATRA-Induced Differentiation.

BACKGROUND: All-trans retinoic acid (ATRA), a derivate of vitamin A, has been successfully used as a therapy to induce differentiation in M3 acute promyelocytic leukemia (APML), and has led to marked improvement in outcomes. Previously, attempts to use ATRA in non-APML in the clinic, however, have been underwhelming, likely due to persistent signaling through other oncogenic drivers. Dysregulated JAK/STAT signaling is known to drive several hematologic malignancies, and targeting JAK1 and JAK2 with the JAK1/JAK2 inhibitor ruxolitinib has led to improvement in survival in primary myelofibrosis and alleviation of vasomotor symptoms and splenomegaly in polycythemia vera and myelofibrosis. OBJECTIVE: While dose-dependent anemia and thrombocytopenia limit the use of JAK2 inhibition, selectively targeting JAK1 has been explored as a means to suppress inflammation and STAT-associated pathologies related to neoplastogenesis. The objective of this study is to employ JAK1 inhibition (JAK1i) in the presence of ATRA as a potential therapy in non-M3 acute myeloid leukemia (AML). METHODS: Efficacy of JAK1i using INCB52793 was assessed by changes in cell cycle and apoptosis in treated AML cell lines. Transcriptomic and proteomic analysis evaluated effects of JAK1i. Synergy between JAK1i+ ATRA was assessed in cell lines in vitro while efficacy in vivo was assessed by tumor reduction in MV-4-11 cell line-derived xenografts. RESULTS: Here we describe novel synergistic activity between JAK1i inhibition and ATRA in non-M3 leukemia. Transcriptomic and proteomic analysis confirmed structural and functional changes related to maturation while in vivo combinatory studies revealed significant decreases in leukemic expansion. CONCLUSIONS: JAK1i+ ATRA lead to decreases in cell cycle followed by myeloid differentiation and cell death in human leukemias. These findings highlight potential uses of ATRA-based differentiation therapy of non-M3 human leukemia.

GPU-powered model analysis with PySB/cupSODA.

SUMMARY: A major barrier to the practical utilization of large, complex models of biochemical systems is the lack of open-source computational tools to evaluate model behaviors over high-dimensional parameter spaces. This is due to the high computational expense of performing thousands to millions of model simulations required for statistical analysis. To address this need, we have implemented a user-friendly interface between cupSODA, a GPU-powered kinetic simulator, and PySB, a Python-based modeling and simulation framework. For three example models of varying size, we show that for large numbers of simulations PySB/cupSODA achieves order-of-magnitude speedups relative to a CPU-based ordinary differential equation integrator. AVAILABILITY AND IMPLEMENTATION: The PySB/cupSODA interface has been integrated into the PySB modeling framework (version 1.4.0), which can be installed from the Python Package Index (PyPI) using a Python package manager such as pip. cupSODA source code and precompiled binaries (Linux, Mac OS/X, Windows) are available at github.com/aresio/cupSODA (requires an Nvidia GPU; developer.nvidia.com/cuda-gpus). Additional information about PySB is available at pysb.org. CONTACT: paolo.cazzaniga@unibg.it or c.lopez@vanderbilt.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Activated Oncogenic Pathway Modifies Iron Network in Breast Epithelial Cells: A Dynamic Modeling Perspective.

Dysregulation of iron metabolism in cancer is well documented and it has been suggested that there is interdependence between excess iron and increased cancer incidence and progression. In an effort to better understand the linkages between iron metabolism and breast cancer, a predictive mathematical model of an expanded iron homeostasis pathway was constructed that includes species involved in iron utilization, oxidative stress response and oncogenic pathways. The model leads to three predictions. The first is that overexpression of iron regulatory protein 2 (IRP2) recapitulates many aspects of the alterations in free iron and iron-related proteins in cancer cells without affecting the oxidative stress response or the oncogenic pathways included in the model. This prediction was validated by experimentation. The second prediction is that iron-related proteins are dramatically affected by mitochondrial ferritin overexpression. This prediction was validated by results in the pertinent literature not used for model construction. The third prediction is that oncogenic Ras pathways contribute to altered iron homeostasis in cancer cells. This prediction was validated by a combination of simulation experiments of Ras overexpression and catalase knockout in conjunction with the literature. The model successfully captures key aspects of iron metabolism in breast cancer cells and provides a framework upon which more detailed models can be built.

Integrated, High-Throughput, Multiomics Platform Enables Data-Driven Construction of Cellular Responses and Reveals Global Drug Mechanisms of Action.

An understanding of how cells respond to perturbation is essential for biological applications; however, most approaches for profiling cellular response are limited in scope to pre-established targets. Global analysis of molecular mechanism will advance our understanding of the complex networks constituting cellular perturbation and lead to advancements in areas, such as infectious disease pathogenesis, developmental biology, pathophysiology, pharmacology, and toxicology. We have developed a high-throughput multiomics platform for comprehensive, de novo characterization of cellular mechanisms of action. Platform validation using cisplatin as a test compound demonstrates quantification of over 10 000 unique, significant molecular changes in less than 30 days. These data provide excellent coverage of known cisplatin-induced molecular changes and previously unrecognized insights into cisplatin resistance. This proof-of-principle study demonstrates the value of this platform as a resource to understand complex cellular responses in a high-throughput manner.

Folding propensity of intrinsically disordered proteins by osmotic stress.

Proteins imparted with intrinsic disorder conduct a range of essential cellular functions. To better understand the folding and hydration properties of intrinsically disordered proteins (IDPs), we used osmotic stress to induce conformational changes in nuclear co-activator binding domain (NCBD) and activator for thyroid hormone and retinoid receptor (ACTR) separate from their mutual binding. Osmotic stress was applied by the addition of small and polymeric osmolytes, where we discovered that water contributions to NCBD folding always exceeded those for ACTR. Both NCBD and ACTR were found to gain α-helical structure with increasing osmotic stress, consistent with their folding upon NCBD/ACTR complex formation. Using small-angle neutron scattering (SANS), we further characterized NCBD structural changes with the osmolyte ethylene glycol. Here a large reduction in overall size initially occurred before substantial secondary structural change. By focusing on folding propensity, and linked hydration changes, we uncover new insights that may be important for how IDP folding contributes to binding.

CASP11--An Evaluation of a Modular BCL::Fold-Based Protein Structure Prediction Pipeline.

In silico prediction of a protein's tertiary structure remains an unsolved problem. The community-wide Critical Assessment of Protein Structure Prediction (CASP) experiment provides a double-blind study to evaluate improvements in protein structure prediction algorithms. We developed a protein structure prediction pipeline employing a three-stage approach, consisting of low-resolution topology search, high-resolution refinement, and molecular dynamics simulation to predict the tertiary structure of proteins from the primary structure alone or including distance restraints either from predicted residue-residue contacts, nuclear magnetic resonance (NMR) nuclear overhauser effect (NOE) experiments, or mass spectroscopy (MS) cross-linking (XL) data. The protein structure prediction pipeline was evaluated in the CASP11 experiment on twenty regular protein targets as well as thirty-three 'assisted' protein targets, which also had distance restraints available. Although the low-resolution topology search module was able to sample models with a global distance test total score (GDT_TS) value greater than 30% for twelve out of twenty proteins, frequently it was not possible to select the most accurate models for refinement, resulting in a general decay of model quality over the course of the prediction pipeline. In this study, we provide a detailed overall analysis, study one target protein in more detail as it travels through the protein structure prediction pipeline, and evaluate the impact of limited experimental data.