COMO: a pipeline for multi-omics data integration in metabolic modeling and drug discovery.

Identifying potential drug targets using metabolic modeling requires integrating multiple modeling methods and heterogeneous biological datasets, which can be challenging without efficient tools. We developed Constraint-based Optimization of Metabolic Objectives (COMO), a user-friendly pipeline that integrates multi-omics data processing, context-specific metabolic model development, simulations, drug databases and disease data to aid drug discovery. COMO can be installed as a Docker Image or with Conda and includes intuitive instructions within a Jupyter Lab environment. It provides a comprehensive solution for the integration of bulk and single-cell RNA-seq, microarrays and proteomics outputs to develop context-specific metabolic models. Using public databases, open-source solutions for model construction and a streamlined approach for predicting repurposable drugs, COMO enables researchers to investigate low-cost alternatives and novel disease treatments. As a case study, we used the pipeline to construct metabolic models of B cells, which simulate and analyze them to predict metabolic drug targets for rheumatoid arthritis and systemic lupus erythematosus, respectively. COMO can be used to construct models for any cell or tissue type and identify drugs for any human disease where metabolic inhibition is relevant. The pipeline has the potential to improve the health of the global community cost-effectively by providing high-confidence targets to pursue in preclinical and clinical studies. The source code of the COMO pipeline is available at https://github.com/HelikarLab/COMO. The Docker image can be pulled at https://github.com/HelikarLab/COMO/pkgs/container/como.

A multiscale mechanistic model of human dendritic cells for in-silico investigation of immune responses and novel therapeutics discovery.

Dendritic cells (DCs) are professional antigen-presenting cells (APCs) with the unique ability to mediate inflammatory responses of the immune system. Given the critical role of DCs in shaping immunity, they present an attractive avenue as a therapeutic target to program the immune system and reverse immune disease disorders. To ensure appropriate immune response, DCs utilize intricate and complex molecular and cellular interactions that converge into a seamless phenotype. Computational models open novel frontiers in research by integrating large-scale interaction to interrogate the influence of complex biological behavior across scales. The ability to model large biological networks will likely pave the way to understanding any complex system in more approachable ways. We developed a logical and predictive model of DC function that integrates the heterogeneity of DCs population, APC function, and cell-cell interaction, spanning molecular to population levels. Our logical model consists of 281 components that connect environmental stimuli with various layers of the cell compartments, including the plasma membrane, cytoplasm, and nucleus to represent the dynamic processes within and outside the DC, such as signaling pathways and cell-cell interactions. We also provided three sample use cases to apply the model in the context of studying cell dynamics and disease environments. First, we characterized the DC response to Sars-CoV-2 and influenza co-infection by in-silico experiments and analyzed the activity level of 107 molecules that play a role in this co-infection. The second example presents simulations to predict the crosstalk between DCs and T cells in a cancer microenvironment. Finally, for the third example, we used the Kyoto Encyclopedia of Genes and Genomes enrichment analysis against the model's components to identify 45 diseases and 24 molecular pathways that the DC model can address. This study presents a resource to decode the complex dynamics underlying DC-derived APC communication and provides a platform for researchers to perform in-silico experiments on human DC for vaccine design, drug discovery, and immunotherapies.

Tim-3 mediates T cell trogocytosis to limit antitumor immunity.

T cell immunoglobulin mucin domain-containing protein 3 (Tim-3) negatively regulates innate and adaptive immunity in cancer. To identify the mechanisms of Tim-3 in cancer immunity, we evaluated the effects of Tim-3 blockade in human and mouse melanoma. Here, we show that human programmed cell death 1-positive (PD-1+) Tim-3+CD8+ tumor-infiltrating lymphocytes (TILs) upregulate phosphatidylserine (PS), a receptor for Tim-3, and acquire cell surface myeloid markers from antigen-presenting cells (APCs) through transfer of membrane fragments called trogocytosis. Tim-3 blockade acted on Tim-3+ APCs in a PS-dependent fashion to disrupt the trogocytosis of activated tumor antigen-specific CD8+ T cells and PD-1+Tim-3+ CD8+ TILs isolated from patients with melanoma. Tim-3 and PD-1 blockades cooperated to disrupt trogocytosis of CD8+ TILs in 2 melanoma mouse models, decreasing tumor burden and prolonging survival. Deleting Tim-3 in dendritic cells but not in CD8+ T cells impeded the trogocytosis of CD8+ TILs in vivo. Trogocytosed CD8+ T cells presented tumor peptide-major histocompatibility complexes and became the target of fratricide T cell killing, which was reversed by Tim-3 blockade. Our findings have uncovered a mechanism Tim-3 uses to limit antitumor immunity.

The follicular lymphoma epigenome regulates its microenvironment.

Follicular lymphoma (FL) is a B-cell non-Hodgkin lymphoma of germinal center (GC) origin with a distinctive tumor microenvironment (TME) and a unique spectrum of mutations. Despite the important therapeutic advances, FL is still incurable. During B-cell development, the GC reaction is a complex multistep process in which epigenetic regulators dynamically induce or suppress transcriptional programs. In FL, epigenetic gene mutations perturb the regulation of these programs, changing GC B-cell function and skewing differentiation towards tumor cells and altering the microenvironment interactions. FL pathogenesis and malignant transformation are promoted by epigenetic reprogramming of GC B cells that alters the immunological synapse and niche. Despite the extensive characterization of FL epigenetic signature and TME, the functional consequences of epigenetic dysregulation on TME and niche plasticity need to be better characterized. In this review, first we describe the most frequent epigenomic alterations in FL (KMT2D, CREBBP and EZH2) that affect the immunological niche, and their potential consequences on the informational transfer between tumor B cells and their microenvironment. Then, we discuss the latest progress to harness epigenetic targets for inhibiting the FL microenvironment. Finally, we highlight unexplored research areas and outstanding questions that should be considered for a successful long-term treatment of FL.

Recent applications of quantitative systems pharmacology and machine learning models across diseases.

Quantitative systems pharmacology (QSP) is a quantitative and mechanistic platform describing the phenotypic interaction between drugs, biological networks, and disease conditions to predict optimal therapeutic response. In this meta-analysis study, we review the utility of the QSP platform in drug development and therapeutic strategies based on recent publications (2019-2021). We gathered recent original QSP models and described the diversity of their applications based on therapeutic areas, methodologies, software platforms, and functionalities. The collection and investigation of these publications can assist in providing a repository of recent QSP studies to facilitate the discovery and further reusability of QSP models. Our review shows that the largest number of QSP efforts in recent years is in Immuno-Oncology. We also addressed the benefits of integrative approaches in this field by presenting the applications of Machine Learning methods for drug discovery and QSP models. Based on this meta-analysis, we discuss the advantages and limitations of QSP models and propose fields where the QSP approach constitutes a valuable interface for more investigations to tackle complex diseases and improve drug development.

Follicular lymphoma triggers phenotypic and functional remodeling of the human lymphoid stromal cell landscape.

Lymphoid stromal cells (LSCs) are essential organizers of immune responses. We analyzed tonsillar tissue by combining flow cytometry, in situ imaging, RNA sequencing, and functional assays, defining three distinct human LSC subsets. The integrin CD49a designated perivascular stromal cells exhibiting features of local committed LSC precursors and segregated cytokine and chemokine-producing fibroblastic reticular cells (FRCs) supporting B and T cell survival. The follicular dendritic cell transcriptional profile reflected active responses to B cell and non-B cell stimuli. We therefore examined the effect of B cell stimuli on LSCs in follicular lymphoma (FL). FL B cells interacted primarily with CD49a+ FRCs. Transcriptional analyses revealed LSC reprogramming in situ downstream of the cytokines tumor necrosis factor (TNF) and transforming growth factor ß (TGF-ß), including increased expression of the chemokines CCL19 and CCL21. Our findings define human LSC populations in healthy tissue and reveal bidirectional crosstalk between LSCs and malignant B cells that may present a targetable axis in lymphoma.

Integrative computational approach identifies drug targets in CD4+ T-cell-mediated immune disorders.

CD4+ T cells provide adaptive immunity against pathogens and abnormal cells, and they are also associated with various immune-related diseases. CD4+ T cells' metabolism is dysregulated in these pathologies and represents an opportunity for drug discovery and development. Genome-scale metabolic modeling offers an opportunity to accelerate drug discovery by providing high-quality information about possible target space in the context of a modeled disease. Here, we develop genome-scale models of naïve, Th1, Th2, and Th17 CD4+ T-cell subtypes to map metabolic perturbations in rheumatoid arthritis, multiple sclerosis, and primary biliary cholangitis. We subjected these models to in silico simulations for drug response analysis of existing FDA-approved drugs and compounds. Integration of disease-specific differentially expressed genes with altered reactions in response to metabolic perturbations identified 68 drug targets for the three autoimmune diseases. In vitro experimental validation, together with literature-based evidence, showed that modulation of fifty percent of identified drug targets suppressed CD4+ T cells, further increasing their potential impact as therapeutic interventions. Our approach can be generalized in the context of other diseases, and the metabolic models can be further used to dissect CD4+ T-cell metabolism.

Nuclear Heparanase Regulates Chromatin Remodeling, Gene Expression and PTEN Tumor Suppressor Function.

Heparanase (HPSE) is an endoglycosidase that cleaves heparan sulfate and has been shown in various cancers to promote metastasis, angiogenesis, osteolysis, and chemoresistance. Although heparanase is thought to act predominantly extracellularly or within the cytoplasm, it is also present in the nucleus, where it may function in regulating gene transcription. Using myeloma cell lines, we report here that heparanase enhances chromatin accessibility and confirm a previous report that it also upregulates the acetylation of histones. Employing the Multiple Myeloma Research Foundation CoMMpass database, we demonstrate that patients expressing high levels of heparanase display elevated expression of proteins involved in chromatin remodeling and several oncogenic factors compared to patients expressing low levels of heparanase. These signatures were consistent with the known function of heparanase in driving tumor progression. Chromatin opening and downstream target genes were abrogated by inhibition of heparanase. Enhanced levels of heparanase in myeloma cells led to a dramatic increase in phosphorylation of PTEN, an event known to stabilize PTEN, leading to its inactivity and loss of tumor suppressor function. Collectively, this study demonstrates that heparanase promotes chromatin opening and transcriptional activity, some of which likely is through its impact on diminishing PTEN tumor suppressor activity.

Heparanase promotes myeloma stemness and in vivo tumorigenesis.

Heparanase is known to enhance the progression of many cancer types and is associated with poor patient prognosis. We recently reported that after patients with multiple myeloma were treated with high dose chemotherapy, the tumor cells that emerged upon relapse expressed a much higher level of heparanase than was present prior to therapy. Because tumor cells having stemness properties are thought to seed tumor relapse, we investigated whether heparanase had a role in promoting myeloma stemness. When plated at low density and grown in serum-free conditions that support survival and expansion of stem-like cells, myeloma cells expressing a low level of heparanase formed tumor spheroids poorly. In contrast, cells expressing a high level of heparanase formed significantly more and larger spheroids than did the heparanase low cells. Importantly, heparanase-low expressing cells exhibited plasticity and were induced to exhibit stemness properties when exposed to recombinant heparanase or to exosomes that contained a high level of heparanase cargo. The spheroid-forming heparanase-high cells had elevated expression of GLI1, SOX2 and ALDH1A1, three genes known to be associated with myeloma stemness. Inhibitors that block the heparan sulfate degrading activity of heparanase significantly diminished spheroid formation and expression of stemness genes implying a direct role of the enzyme in regulating stemness. Blocking the NF-κB pathway inhibited spheroid formation and expression of stemness genes demonstrating a role for NF-κB in heparanase-mediated stemness. Myeloma cells made deficient in heparanase exhibited decreased stemness properties in vitro and when injected into mice they formed tumors poorly compared to the robust tumorigenic capacity of cells expressing higher levels of heparanase. These studies reveal for the first time a role for heparanase in promoting cancer stemness and provide new insight into its function in driving tumor progression and its association with poor prognosis in cancer patients.

CD226 opposes TIGIT to disrupt Tregs in melanoma.

CD4+ Tregs impede T cell responses to tumors. They express multiple inhibitory receptors that support their suppressive functions, including T cell Ig and ITIM domain (TIGIT). In melanoma patients, we show that Tregs exhibit increased TIGIT expression and decreased expression of its competing costimulatory receptor CD226 as compared with CD4+ effector T cells, resulting in an increased TIGIT/CD226 ratio. Tregs failed to upregulate CD226 upon T cell activation. TIGIT+ Tregs are highly suppressive, stable, and enriched in tumors. TIGIT and CD226 oppose each other to augment or disrupt, respectively, Treg suppression and stability. A high TIGIT/CD226 ratio in Tregs correlates with increased Treg frequencies in tumors and poor clinical outcome upon immune checkpoint blockade. Altogether, our findings show that a high TIGIT/CD226 ratio in Tregs regulates their suppressive function and stability in melanoma. They provide the rationale for novel immunotherapies to activate CD226 in Tregs together with TIGIT blockade to counteract Treg suppression in cancer patients.

Loss of the HVEM Tumor Suppressor in Lymphoma and Restoration by Modified CAR-T Cells.

The HVEM (TNFRSF14) receptor gene is among the most frequently mutated genes in germinal center lymphomas. We report that loss of HVEM leads to cell-autonomous activation of B cell proliferation and drives the development of GC lymphomas in vivo. HVEM-deficient lymphoma B cells also induce a tumor-supportive microenvironment marked by exacerbated lymphoid stroma activation and increased recruitment of T follicular helper (TFH) cells. These changes result from the disruption of inhibitory cell-cell interactions between the HVEM and BTLA (B and T lymphocyte attenuator) receptors. Accordingly, administration of the HVEM ectodomain protein (solHVEM(P37-V202)) binds BTLA and restores tumor suppression. To deliver solHVEM to lymphomas in vivo, we engineered CD19-targeted chimeric antigen receptor (CAR) T cells that produce solHVEM locally and continuously. These modified CAR-T cells show enhanced therapeutic activity against xenografted lymphomas. Hence, the HVEM-BTLA axis opposes lymphoma development, and our study illustrates the use of CAR-T cells as "micro-pharmacies" able to deliver an anti-cancer protein.

Loss of IL-22 inhibits autoantibody formation in collagen-induced arthritis in mice.

Interleukin 22 (IL-22) expression is associated with increased joint destruction and disease progression in rheumatoid arthritis (RA). Although IL-22 is considered a pro-inflammatory cytokine, its mechanism of action in RA remains incompletely understood. Here, we used the collagen-induced arthritis model in IL-22 deficient (IL-22(-/-) ) mice to study the role of IL-22 in RA. In spite of normal disease incidence, disease severity is significantly diminished in IL-22(-/-) mice. Moreover, pathogenicity of Th17 cells and development and function of B cells are unaffected. In contrast, splenic plasma cells, as well as serum autoantibody titers, are reduced in the absence of IL-22. At the peak of disease, germinal centers (GCs) are severely reduced in the spleens of IL-22(-/-) mice, correlating with a decline in GC B-cell numbers. Within the GC, we identified IL-22R1 expressing follicular dendritic cell-like stromal cells. Human lymphoid stromal cells respond to IL-22 ex vivo by inducing transcription of CXCL12 and CXCL13. We therefore postulate IL-22 as an important enhancer of the GC reaction, maintaining chemokine levels for the persistence of GC reactions, essential for the production of autoantibody-secreting plasma cells. Blocking IL-22 might therefore prevent immune-complex deposition and destruction of joints in RA patients.

DC-SIGN-expressing macrophages trigger activation of mannosylated IgM B-cell receptor in follicular lymphoma.

Follicular lymphoma (FL) results from the accumulation of malignant germinal center (GC) B cells leading to the development of an indolent and largely incurable disease. FL cells remain highly dependent on B-cell receptor (BCR) signaling and on a specific cell microenvironment, including T cells, macrophages, and stromal cells. Importantly, FL BCR is characterized by a selective pressure to retain surface immunoglobulin M (IgM) BCR despite an active class-switch recombination process, and by the introduction, in BCR variable regions, of N-glycosylation acceptor sites harboring unusual high-mannose oligosaccharides. However, the relevance of these 2 FL BCR features for lymphomagenesis remains unclear. In this study, we demonstrated that IgM(+) FL B cells activated a stronger BCR signaling network than IgG(+) FL B cells and normal GC B cells. BCR expression level and phosphatase activity could both contribute to such heterogeneity. Moreover, we underlined that a subset of IgM(+) FL samples, displaying highly mannosylated BCR, efficiently bound dendritic cell-specific intercellular adhesion molecule-3-grabbing nonintegrin (DC-SIGN), which could in turn trigger delayed but long-lasting BCR aggregation and activation. Interestingly, DC-SIGN was found within the FL cell niche in situ. Finally, M2 macrophages induced a DC-SIGN-dependent adhesion of highly mannosylated IgM(+) FL B cells and triggered BCR-associated kinase activation. Interestingly, pharmacologic BCR inhibitors abolished such crosstalk between macrophages and FL B cells. Altogether, our data support an important role for DC-SIGN-expressing infiltrating cells in the biology of FL and suggest that they could represent interesting therapeutic targets.

The class-specific BCR tonic signal modulates lymphomagenesis in a c-myc deregulation transgenic model.

Deregulation of c-myc by translocation onto immunoglobulin (Ig) loci can promote B cell malignant proliferations with phenotypes as diverse as acute lymphoid leukemia, Burkitt lymphoma, diffuse large B cell lymphoma, myeloma... The B cell receptor (BCR) normally providing tonic signals for cell survival and mitogenic responses to antigens, can also contribute to lymphomagenesis upon sustained ligand binding or activating mutations. BCR signaling varies among cell compartments and BCR classes. For unknown reasons, some malignancies associate with expression of either IgM or class-switched Ig. We explored whether an IgA BCR, with strong tonic signaling, would affect lymphomagenesis in c-myc IgH 3'RR transgenic mice prone to lymphoproliferations. Breeding c-myc transgenics in a background where IgM expression was replaced with IgA delayed lymphomagenesis. By comparison to single c-myc transgenics, lymphomas from double mutant animals were more differentiated and less aggressive, with an altered transcriptional program. Larger tumor cells more often expressed CD43 and CD138, which culminated in a plasma cell phenotype in 10% of cases. BCR class-specific signals thus appear to modulate lymphomagenesis and may partly explain the observed association of specific Ig classes with human B cell malignancies of differential phenotype, progression and prognosis.

Deletion of the α immunoglobulin chain membrane-anchoring region reduces but does not abolish IgA secretion.

Class switching and plasma cell differentiation occur at a high level within all mucosa-associated lymphoid tissues. The different classes of membrane immunoglobulin heavy chains are associated with the Igα/Igß heterodimer within the B-cell receptor (BCR). Whether BCR isotypes convey specific signals adapted to the corresponding differentiation stages remains debated but IgG and IgA membranes have been suggested to promote plasma cell differentiation. We investigated the impact of blocking expression of the IgA-class BCR through a 'αΔtail' targeted mutation, deleting the Cα immunoglobulin gene membrane exon. This allowed us to evaluate to what extent class switching and plasma cell differentiation can be concurrent processes, allowing some αΔtail(+/+) B cells with an IgM BCR to directly differentiate into IgA plasma cells and yield serum secreted IgA in spite of the absence of membrane IgA(+) B lymphocytes. By contrast, in secretions the secretory IgA was very low, indicating that J-chain-positive plasma cells producing secretory IgA overwhelmingly differentiate from previously class-switched membrane IgA(+) memory B cells. In addition, although mucosa-associated lymphoid tissues are a major site for plasma cell accumulation, αΔtail(+/+) mice showed that the gut B-cell lineage homeostasis is not polarized toward plasma cell differentiation through a specific influence of the membrane IgA BCR.

Premature replacement of mu with alpha immunoglobulin chains impairs lymphopoiesis and mucosal homing but promotes plasma cell maturation.

Sequentially along B cell differentiation, the different classes of membrane Ig heavy chains associate with the Ig alpha/Ig beta heterodimer within the B cell receptor (BCR). Whether each Ig class conveys specific signals adapted to the corresponding differentiation stage remains debated. We investigated the impact of the forced expression of an IgA-class receptor throughout murine B cell differentiation by knocking in the human C alpha Ig gene in place of the S mu region. Despite expression of a functional BCR, homozygous mutant mice showed a partial developmental blockade at the pro-B/pre-BI and large pre-BII cell stages, with decreased numbers of small pre-BII cells. Beyond this stage, peripheral B cell compartments of reduced size developed and allowed specific antibody responses, whereas mature cells showed constitutive activation and a strong commitment to plasma cell differentiation. Secreted IgA correctly assembled into polymers, associated with the murine J chain, and was transported into secretions. In heterozygous mutants, cells expressing the IgA allele competed poorly with those expressing IgM from the wild-type allele and were almost undetectable among peripheral B lymphocytes, notably in gut-associated lymphoid tissues. Our data indicate that the IgM BCR is more efficient in driving early B cell education and in mucosal site targeting, whereas the IgA BCR appears particularly suited to promoting activation and differentiation of effector plasma cells.

A myeloma translocation-like model associating CCND1 with the immunoglobulin heavy-chain locus 3' enhancers does not promote by itself B-cell malignancies.

Cyclin D1 overexpression is associated with mantle cell lymphoma and multiple myeloma. In myeloma, it often results from chromosomal translocations linking the CCND1 gene to the 3' part of the IgH locus constant region. This region includes a single and potent transcriptional regulatory region (RR) 3' of the Calpha gene mostly active in mature B-cells. To check whether this RR alone was sufficient to deregulate CCND1, we generated mice carrying a 3'IgH RR-driven human CCND1 transgene and specifically up-regulating cyclin D1 expression in B-cells. In transgenic B-cells, cyclin D1 enforced cell cycle entry in response to various stimuli (LPS, anti-IgM, anti-CD40) but also increased cell death, so that exaggerated proliferation did not result in peripheral lymphocytosis. Despite exaggerated B-cell entry into G(1) phase, malignant lymphoproliferation did not occur either. Crossing of CCND1-3'IgH RR mice with c-myc-3'IgH RR mice did not reveal accelerated tumorigenesis as compared with c-myc-3'IgH RR mice alone. The data presented here demonstrate that the 3'IgH RR-mediated deregulation of CCND1 in mature B-cells cannot by itself trigger the development of lymphomas and strengthen the concept that cyclin D1 per se is not an armful proto-oncogene. Rather its overexpression in several malignancies might be only a stigma of lymphomagenesis or represent a single hit within a multiple hit process.

ABCA2 is a marker of neural progenitors and neuronal subsets in the adult rodent brain.

The notion that the ATP-binding cassette transporter-A2 (ABCA2) may be involved in brain sterol homeostasis and is associated with early onset Alzheimer's disease led us to explore its neural expression. Our data support and extend the previous reports on ABCA2 expression by oligodendrocytes. They evidence that ABCA2 (i) is located in intracellular vesicles, identified in transfected cells as lysosome-related organelles only partially overlapping with classical endolysosomes; (ii) is a marker of neural progenitors as it is expressed in the subventricular zone of the lateral ventricle and the dentate gyrus of the hippocampal formation, sites of continual neurogenesis in the adult brain, and in nestin(+) cells differentiated in vitro from embryonic stem cells; (iii) persists, in the adult rodent brain, in a subset of GABAergic and glutamatergic neurons. Considering that the latter are targets of Alzheimer's lesions, these data provide a new rationale to explore the neuropathological consequences of ABCA2 functional dysregulations.