Gain of function of a metalloproteinase associated with multiple myeloma, bicuspid aortic valve, and Von Hippel-Lindau syndrome.

A patient diagnosed with multiple myeloma, bicuspid aortic valve, and Von Hippel-Lindau syndrome underwent whole-exome sequencing seeking a unified genetic cause for these three pathologies. The patient possessed a single-point mutation of arginine to cysteine (R24C) in the N-terminal region(pro-domain) of matrix metalloproteinase 9 (MMP-9). The pro-domain interacts with the catalytic site of this enzyme rendering it inactive. MMP-9 has previously been associated with all three pathologies suffered by the patient. We hypothesized that the observed mutation in the pro-domain would influence the activity of this enzyme. We expressed recombinant versions of MMP-9 and an investigation of their biochemical properties revealed that MMP-9 R24C is a constitutively active zymogen. To our knowledge, this is the first example of a mutation that discloses catalytic activity in the pro-form in any of the 24 human MMPs.

Polled Digital Cell Sorter (p-DCS): Automatic identification of hematological cell types from single cell RNA-sequencing clusters.

BACKGROUND: Single cell RNA sequencing (scRNA-seq) brings unprecedented opportunities for mapping the heterogeneity of complex cellular environments such as bone marrow, and provides insight into many cellular processes. Single cell RNA-seq has a far larger fraction of missing data reported as zeros (dropouts) than traditional bulk RNA-seq, and unsupervised clustering combined with Principal Component Analysis (PCA) can be used to overcome this limitation. After clustering, however, one has to interpret the average expression of markers on each cluster to identify the corresponding cell types, and this is normally done by hand by an expert curator. RESULTS: We present a computational tool for processing single cell RNA-seq data that uses a voting algorithm to automatically identify cells based on approval votes received by known molecular markers. Using a stochastic procedure that accounts for imbalances in the number of known molecular signatures for different cell types, the method computes the statistical significance of the final approval score and automatically assigns a cell type to clusters without an expert curator. We demonstrate the utility of the tool in the analysis of eight samples of bone marrow from the Human Cell Atlas. The tool provides a systematic identification of cell types in bone marrow based on a list of markers of immune cell types, and incorporates a suite of visualization tools that can be overlaid on a t-SNE representation. The software is freely available as a Python package at https://github.com/sdomanskyi/DigitalCellSorter . CONCLUSIONS: This methodology assures that extensive marker to cell type matching information is taken into account in a systematic way when assigning cell clusters to cell types. Moreover, the method allows for a high throughput processing of multiple scRNA-seq datasets, since it does not involve an expert curator, and it can be applied recursively to obtain cell sub-types. The software is designed to allow the user to substitute the marker to cell type matching information and apply the methodology to different cellular environments.

Cell cycle time series gene expression data encoded as cyclic attractors in Hopfield systems.

Modern time series gene expression and other omics data sets have enabled unprecedented resolution of the dynamics of cellular processes such as cell cycle and response to pharmaceutical compounds. In anticipation of the proliferation of time series data sets in the near future, we use the Hopfield model, a recurrent neural network based on spin glasses, to model the dynamics of cell cycle in HeLa (human cervical cancer) and S. cerevisiae cells. We study some of the rich dynamical properties of these cyclic Hopfield systems, including the ability of populations of simulated cells to recreate experimental expression data and the effects of noise on the dynamics. Next, we use a genetic algorithm to identify sets of genes which, when selectively inhibited by local external fields representing gene silencing compounds such as kinase inhibitors, disrupt the encoded cell cycle. We find, for example, that inhibiting the set of four kinases AURKB, NEK1, TTK, and WEE1 causes simulated HeLa cells to accumulate in the M phase. Finally, we suggest possible improvements and extensions to our model.

Evolutionary and Topological Properties of Genes and Community Structures in Human Gene Regulatory Networks.

The diverse, specialized genes present in today's lifeforms evolved from a common core of ancient, elementary genes. However, these genes did not evolve individually: gene expression is controlled by a complex network of interactions, and alterations in one gene may drive reciprocal changes in its proteins' binding partners. Like many complex networks, these gene regulatory networks (GRNs) are composed of communities, or clusters of genes with relatively high connectivity. A deep understanding of the relationship between the evolutionary history of single genes and the topological properties of the underlying GRN is integral to evolutionary genetics. Here, we show that the topological properties of an acute myeloid leukemia GRN and a general human GRN are strongly coupled with its genes' evolutionary properties. Slowly evolving ("cold"), old genes tend to interact with each other, as do rapidly evolving ("hot"), young genes. This naturally causes genes to segregate into community structures with relatively homogeneous evolutionary histories. We argue that gene duplication placed old, cold genes and communities at the center of the networks, and young, hot genes and communities at the periphery. We demonstrate this with single-node centrality measures and two new measures of efficiency, the set efficiency and the interset efficiency. We conclude that these methods for studying the relationships between a GRN's community structures and its genes' evolutionary properties provide new perspectives for understanding evolutionary genetics.

Naturally occurring combinations of receptors from single cell transcriptomics in endothelial cells.

VEGF inhibitor drugs are part of standard care in oncology and ophthalmology, but not all patients respond to them. Combinations of drugs are likely to be needed for more effective therapies of angiogenesis-related diseases. In this paper we describe naturally occurring combinations of receptors in endothelial cells that might help to understand how cells communicate and to identify targets for drug combinations. We also develop and share a new software tool called DECNEO to identify them. Single-cell gene expression data are used to identify a set of co-expressed endothelial cell receptors, conserved among species (mice and humans) and enriched, within a network, of connections to up-regulated genes. This set includes several receptors previously shown to play a role in angiogenesis. Multiple statistical tests from large datasets, including an independent validation set, support the reproducibility, evolutionary conservation and role in angiogenesis of these naturally occurring combinations of receptors. We also show tissue-specific combinations and, in the case of choroid endothelial cells, consistency with both well-established and recent experimental findings, presented in a separate paper. The results and methods presented here advance the understanding of signaling to endothelial cells. The methods are generally applicable to the decoding of intercellular combinations of signals.

LIF, a mitogen for choroidal endothelial cells, protects the choriocapillaris: implications for prevention of geographic atrophy.

In the course of our studies aiming to discover vascular bed-specific endothelial cell (EC) mitogens, we identified leukemia inhibitory factor (LIF) as a mitogen for bovine choroidal EC (BCE), although LIF has been mainly characterized as an EC growth inhibitor and an anti-angiogenic molecule. LIF stimulated growth of BCE while it inhibited, as previously reported, bovine aortic EC (BAE) growth. The JAK-STAT3 pathway mediated LIF actions in both BCE and BAE cells, but a caspase-independent proapoptotic signal mediated by cathepsins was triggered in BAE but not in BCE. LIF administration directly promoted activation of STAT3 and increased blood vessel density in mouse eyes. LIF also had protective effects on the choriocapillaris in a model of oxidative retinal injury. Analysis of available single-cell transcriptomic datasets shows strong expression of the specific LIF receptor in mouse and human choroidal EC. Our data suggest that LIF administration may be an innovative approach to prevent atrophy associated with AMD, through protection of the choriocapillaris.

Digital Cell Sorter (DCS): a cell type identification, anomaly detection, and Hopfield landscapes toolkit for single-cell transcriptomics.

MOTIVATION: Analysis of singe cell RNA sequencing (scRNA-seq) typically consists of different steps including quality control, batch correction, clustering, cell identification and characterization, and visualization. The amount of scRNA-seq data is growing extremely fast, and novel algorithmic approaches improving these steps are key to extract more biological information. Here, we introduce: (i) two methods for automatic cell type identification (i.e., without expert curator) based on a voting algorithm and a Hopfield classifier, (ii) a method for cell anomaly quantification based on isolation forest, and (iii) a tool for the visualization of cell phenotypic landscapes based on Hopfield energy-like functions. These new approaches are integrated in a software platform that includes many other state-of-the-art methodologies and provides a self-contained toolkit for scRNA-seq analysis. RESULTS: We present a suite of software elements for the analysis of scRNA-seq data. This Python-based open source software, Digital Cell Sorter (DCS), consists in an extensive toolkit of methods for scRNA-seq analysis. We illustrate the capability of the software using data from large datasets of peripheral blood mononuclear cells (PBMC), as well as plasma cells of bone marrow samples from healthy donors and multiple myeloma patients. We test the novel algorithms by evaluating their ability to deconvolve cell mixtures and detect small numbers of anomalous cells in PBMC data. AVAILABILITY: The DCS toolkit is available for download and installation through the Python Package Index (PyPI). The software can be deployed using the Python import function following installation. Source code is also available for download on Zenodo: DOI 10.5281/zenodo.2533377. SUPPLEMENTARY INFORMATION: Supplemental Materials are available at PeerJ online.

Integrating metabolomics and phenomics with systems models of cardiac hypoxia.

Hypoxia is the major cause of necrotic cell death in myocardial infarction. Cellular energy supply and demand under hypoxic conditions is regulated by many interacting signaling and transcriptional networks, which complicates studies on individual proteins and pathways. We apply an integrated systems approach to understand the metabolic and functional response to hypoxia in muscle cells of the fruit fly Drosophila melanogaster. In addition to its utility as a hypoxia-tolerant model organism, Drosophila also offers advantages due to its small size, fecundity, and short life cycle. These traits, along with a large library of single-gene mutations, motivated us to develop new, computer-automated technology for gathering in vivo measurements of heart function under hypoxia for a large number of mutant strains. Phenotype data can be integrated with in silico cellular networks, metabolomic data, and microarrays to form qualitative and quantitative network models for prediction and hypothesis generation. Here we present a framework for a systems approach to hypoxia in the cardiac myocyte, starting from nuclear magnetic resonance (NMR) metabolomics, a constraint-based metabolic model, and phenotypic profiles.

Metabolomic and flux-balance analysis of age-related decline of hypoxia tolerance in Drosophila muscle tissue.

The fruitfly Drosophila melanogaster is increasingly used as a model organism for studying acute hypoxia tolerance and for studying aging, but the interactions between these two factors are not well known. Here we show that hypoxia tolerance degrades with age in post-hypoxic recovery of whole-body movement, heart rate and ATP content. We previously used (1)H NMR metabolomics and a constraint-based model of ATP-generating metabolism to discover the end products of hypoxic metabolism in flies and generate hypotheses for the biological mechanisms. We expand the reactions in the model using tissue- and age-specific microarray data from the literature, and then examine metabolomic profiles of thoraxes after 4 h at 0.5% O(2) and after 5 min of recovery in 40- versus 3-day-old flies. Model simulations were constrained to fluxes calculated from these data. Simulations suggest that the decreased ATP production during reoxygenation seen in aging flies can be attributed to reduced recovery of mitochondrial respiration pathways and concomitant overdependence on the acetate production pathway as an energy source.

Search algorithms as a framework for the optimization of drug combinations.

Combination therapies are often needed for effective clinical outcomes in the management of complex diseases, but presently they are generally based on empirical clinical experience. Here we suggest a novel application of search algorithms -- originally developed for digital communication -- modified to optimize combinations of therapeutic interventions. In biological experiments measuring the restoration of the decline with age in heart function and exercise capacity in Drosophila melanogaster, we found that search algorithms correctly identified optimal combinations of four drugs using only one-third of the tests performed in a fully factorial search. In experiments identifying combinations of three doses of up to six drugs for selective killing of human cancer cells, search algorithms resulted in a highly significant enrichment of selective combinations compared with random searches. In simulations using a network model of cell death, we found that the search algorithms identified the optimal combinations of 6-9 interventions in 80-90% of tests, compared with 15-30% for an equivalent random search. These findings suggest that modified search algorithms from information theory have the potential to enhance the discovery of novel therapeutic drug combinations. This report also helps to frame a biomedical problem that will benefit from an interdisciplinary effort and suggests a general strategy for its solution.

Modeling disease progression in Multiple Myeloma with Hopfield networks and single-cell RNA-seq.

Associative memories in Hopfield's neural networks are mapped to gene expression pattern to model different paths of disease progression towards Multiple Myeloma (MM). The model is built using single cell RNA-seq data from bone marrow aspirates of MM patients as well as patients diagnosed with Monoclonal Gammopathy of Undetermined Significance (MGUS) and Smoldering Multiple Myeloma (SMM), two medical conditions that often progress to full MM. Results: We identify different clusters of MGUS, SMM, and MM cells, map them to Hopfield associative memory patterns, and model the dynamics of transition between the different patterns. The model is then used to identify genes that are differentialy expressed across different MM stages and whose simultaneous inhibition is associated to a delayed disease progression.

Flexibility in energy metabolism supports hypoxia tolerance in Drosophila flight muscle: metabolomic and computational systems analysis.

The fruitfly Drosophila melanogaster offers promise as a genetically tractable model for studying adaptation to hypoxia at the cellular level, but the metabolic basis for extreme hypoxia tolerance in flies is not well known. Using (1)H NMR spectroscopy, metabolomic profiles were collected under hypoxia. Accumulation of lactate, alanine, and acetate suggested that these are the major end products of anaerobic metabolism in the fly. A constraint-based model of ATP-producing pathways was built using the annotated genome, existing models, and the literature. Multiple redundant pathways for producing acetate and alanine were added and simulations were run in order to find a single optimal strategy for producing each end product. System-wide adaptation to hypoxia was then investigated in silico using the refined model. Simulations supported the hypothesis that the ability to flexibly convert pyruvate to these three by-products might convey hypoxia tolerance by improving the ATP/H(+) ratio and efficiency of glucose utilization.

Discovering regulators of the Drosophila cardiac hypoxia response using automated phenotyping technology.

Necrosis and apoptosis during acute myocardial infarction result in part from the inability of hypoxic cardiac myocytes to match ATP supply and demand. In contrast, hypoxia-tolerant organisms, such as Drosophila, can rapidly regulate cellular metabolism to survive large oxygen fluctuations. A genetic screen of fly heart function during acute hypoxia can be an unbiased way to discover essential enzymes and novel signaling proteins involved in this response. We have developed a prototype to show proof of concept for a genome-scale screen, using computer automation to rapidly gather in vivo hypoxic heart data in adult Drosophila. Our system automatically anesthetizes flies, deposits them on a microscope slide, and locates the heart organ of each fly. The system then applies a hypoxia stimulus, acquires time-space (M-mode) images of the heart walls, and analyzes heart rate and rhythm. The prototype can produce highly controlled measurements of up to 55 flies per hour, which we demonstrated by characterizing the effect of temperature, oxygen content, and genetic background on the hypoxia response. We discuss the possible applications of a genome-wide cardiac phenotype data set in systems biology analyses of hypoxic metabolism, using genome-scale interaction networks and constraint-based metabolic models.

The nitric oxide scavenger cobinamide profoundly improves survival in a Drosophila melanogaster model of bacterial sepsis.

Septic shock has an extremely high mortality rate, with approximately 200,000 people dying from sepsis annually in the U.S. The high mortality results in part from severe hypotension secondary to high serum NO concentrations. Reducing NO levels should be beneficial in sepsis, but NOS inhibitors have had a checkered history in animal models, and one such agent increased mortality in a clinical trial. An alternative approach to reduce NO levels in sepsis is to use an NO scavenger, which should leave sufficient free NO for normal physiological functions. Using a well-established model of bacterial sepsis in Drosophila melanogaster, we found that cobinamide, a B(12) analog and an effective NO scavenger in vitro, dramatically improved fly survival. Cobinamide augmented the effect of an antibiotic and was beneficial even in immune-deficient flies. Cobinamide's mechanism of action appeared to be from reducing NO levels and improving cardiac function.

High-throughput screening and bioinformatic analysis to ascertain compounds that prevent saturated fatty acid-induced ß-cell apoptosis.

Pancreatic ß-cell lipotoxicity is a central feature of the pathogenesis of type 2 diabetes. To study the mechanism by which fatty acids cause ß-cell death and develop novel approaches to prevent it, a high-throughput screen on the ß-cell line INS1 was carried out. The cells were exposed to palmitate to induce cell death and compounds that reversed palmitate-induced cytotoxicity were ascertained. Hits from the screen were analyzed by an increasingly more stringent testing funnel, ending with studies on primary human islets treated with palmitate. MAP4K4 inhibitors, which were not part of the screening libraries but were ascertained by a bioinformatics analysis, and the endocannabinoid anandamide were effective at inhibiting palmitate-induced apoptosis in INS1 cells as well as primary rat and human islets. These targets could serve as the starting point for the development of therapeutics for type 2 diabetes.

Cardiac systems biology.

As more detailed molecular information accumulates on the biology of the heart and other complex systems in health and disease, the need for new integrative analyses and tools is growing. Systems biology and bioengineering seek to use high-throughput technologies and integrative computational analysis to construct networks of the interactions between molecular components in the system, to develop systems models of their functionally integrated biological properties, and to incorporate these systems models into structurally integrated multi-scale models for predicting clinical phenotypes. This review gives examples of recent applications using these approaches to elucidate the electromechanical function of the heart in aging and disease.

Combinations of Kinase Inhibitors Protecting Myoblasts against Hypoxia.

Cell-based therapies to treat skeletal muscle disease are limited by the poor survival of donor myoblasts, due in part to acute hypoxic stress. After confirming that the microenvironment of transplanted myoblasts is hypoxic, we screened a kinase inhibitor library in vitro and identified five kinase inhibitors that protected myoblasts from cell death or growth arrest in hypoxic conditions. A systematic, combinatorial study of these compounds further improved myoblast viability, showing both synergistic and additive effects. Pathway and target analysis revealed CDK5, CDK2, CDC2, WEE1, and GSK3ß as the main target kinases. In particular, CDK5 was the center of the target kinase network. Using our recently developed statistical method based on elastic net regression we computationally validated the key role of CDK5 in cell protection against hypoxia. This method provided a list of potential kinase targets with a quantitative measure of their optimal amount of relative inhibition. A modified version of the method was also able to predict the effect of combinations using single-drug response data. This work is the first step towards a broadly applicable system-level strategy for the pharmacology of hypoxic damage.

A scalable method for molecular network reconstruction identifies properties of targets and mutations in acute myeloid leukemia.

A key aim of systems biology is the reconstruction of molecular networks. We do not yet, however, have networks that integrate information from all datasets available for a particular clinical condition. This is in part due to the limited scalability, in terms of required computational time and power, of existing algorithms. Network reconstruction methods should also be scalable in the sense of allowing scientists from different backgrounds to efficiently integrate additional data. We present a network model of acute myeloid leukemia (AML). In the current version (AML 2.1), we have used gene expression data (both microarray and RNA-seq) from 5 different studies comprising a total of 771 AML samples and a protein-protein interactions dataset. Our scalable network reconstruction method is in part based on the well-known property of gene expression correlation among interacting molecules. The difficulty of distinguishing between direct and indirect interactions is addressed by optimizing the coefficient of variation of gene expression, using a validated gold-standard dataset of direct interactions. Computational time is much reduced compared to other network reconstruction methods. A key feature is the study of the reproducibility of interactions found in independent clinical datasets. An analysis of the most significant clusters, and of the network properties (intraset efficiency, degree, betweenness centrality, and PageRank) of common AML mutations demonstrated the biological significance of the network. A statistical analysis of the response of blast cells from 11 AML patients to a library of kinase inhibitors provided an experimental validation of the network. A combination of network and experimental data identified CDK1, CDK2, CDK4, and CDK6 and other kinases as potential therapeutic targets in AML.

Control of asymmetric Hopfield networks and application to cancer attractors.

The asymmetric Hopfield model is used to simulate signaling dynamics in gene regulatory networks. The model allows for a direct mapping of a gene expression pattern into attractor states. We analyze different control strategies aimed at disrupting attractor patterns using selective local fields representing therapeutic interventions. The control strategies are based on the identification of signaling bottlenecks, which are single nodes or strongly connected clusters of nodes that have a large impact on the signaling. We provide a theorem with bounds on the minimum number of nodes that guarantee control of bottlenecks consisting of strongly connected components. The control strategies are applied to the identification of sets of proteins that, when inhibited, selectively disrupt the signaling of cancer cells while preserving the signaling of normal cells. We use an experimentally validated non-specific and an algorithmically-assembled specific B cell gene regulatory network reconstructed from gene expression data to model cancer signaling in lung and B cells, respectively. Among the potential targets identified here are TP53, FOXM1, BCL6 and SRC. This model could help in the rational design of novel robust therapeutic interventions based on our increasing knowledge of complex gene signaling networks.

Cellular protection using Flt3 and PI3Kα inhibitors demonstrates multiple mechanisms of oxidative glutamate toxicity.

Glutamate-induced oxidative stress is a major contributor to neurodegenerative diseases. Here, we identify small-molecule inhibitors of this process. We screen a kinase inhibitor library on neuronal cells and identify Flt3 and PI3Kα inhibitors as potent protectors against glutamate toxicity. Both inhibitors prevented reactive oxygen species (ROS) generation, mitochondrial hyperpolarization and lipid peroxidation in neuronal cells, but they do so by distinct molecular mechanisms. The PI3Kα inhibitor protects cells by inducing partial restoration of depleted glutathione levels and accumulation of intracellular amino acids, whereas the Flt3 inhibitor prevents lipid peroxidation, a key mechanism of glutamate-mediated toxicity. We also demonstrate that glutamate toxicity involves a combination of ferroptosis, necrosis and AIF-dependent apoptosis. We confirm the protective effect by using multiple inhibitors of these kinases and multiple cell types. Our results not only identify compounds that protect against glutamate-stimulated oxidative stress, but also provide new insights into the mechanisms of glutamate toxicity in neurons.