Systems biology in immunology: a computational modeling perspective.

Systems biology is an emerging discipline that combines high-content, multiplexed measurements with informatic and computational modeling methods to better understand biological function at various scales. Here we present a detailed review of the methods used to create computational models and to conduct simulations of immune function. We provide descriptions of the key data-gathering techniques employed to generate the quantitative and qualitative data required for such modeling and simulation and summarize the progress to date in applying these tools and techniques to questions of immunological interest, including infectious disease. We include comments on what insights modeling can provide that complement information obtained from the more familiar experimental discovery methods used by most investigators and the reasons why quantitative methods are needed to eventually produce a better understanding of immune system operation in health and disease.

Author Correction: Commensal-driven immune zonation of the liver promotes host defence.

A Correction to this paper has been published: https://doi.org/10.1038/s41586-021-03346-0.

Commensal-driven immune zonation of the liver promotes host defence.

The liver connects the intestinal portal vasculature with the general circulation, using a diverse array of immune cells to protect from pathogens that translocate from the gut1. In liver lobules, blood flows from portal triads that are situated in periportal lobular regions to the central vein via a polarized sinusoidal network. Despite this asymmetry, resident immune cells in the liver are considered to be broadly dispersed across the lobule. This differs from lymphoid organs, in which immune cells adopt spatially biased positions to promote effective host defence2,3. Here we used quantitative multiplex imaging, genetic perturbations, transcriptomics, infection-based assays and mathematical modelling to reassess the relationship between the localization of immune cells in the liver and host protection. We found that myeloid and lymphoid resident immune cells concentrate around periportal regions. This asymmetric localization was not developmentally controlled, but resulted from sustained MYD88-dependent signalling induced by commensal bacteria in liver sinusoidal endothelial cells, which in turn regulated the composition of the pericellular matrix involved in the formation of chemokine gradients. In vivo experiments and modelling showed that this immune spatial polarization was more efficient than a uniform distribution in protecting against systemic bacterial dissemination. Together, these data reveal that liver sinusoidal endothelial cells sense the microbiome, actively orchestrating the localization of immune cells, to optimize host defence.

Anti-Inflammatory Chromatinscape Suggests Alternative Mechanisms of Glucocorticoid Receptor Action.

Despite the widespread use of glucocorticoids (GCs), their anti-inflammatory effects are not understood mechanistically. Numerous investigations have examined the effects of glucocorticoid receptor (GR) activation prior to inflammatory challenges. However, clinical situations are emulated by a GC intervention initiated in the midst of rampant inflammatory responses. To characterize the effects of a late GC treatment, we profiled macrophage transcriptional and chromatinscapes with Dexamethasone (Dex) treatment before or after stimulation by lipopolysaccharide (LPS). The late activation of GR had a similar gene-expression profile as from GR pre-activation, while ameliorating the disruption of metabolic genes. Chromatin occupancy of GR was not predictive of Dex-regulated gene expression, contradicting the "trans-repression by tethering" model. Rather, GR activation resulted in genome-wide blockade of NF-κB interaction with chromatin and directly induced inhibitors of NF-κB and AP-1. Our investigation using GC treatments with clinically relevant timing highlights mechanisms underlying GR actions for modulating the "inflamed epigenome."

High-throughput imaging of mRNA at the single-cell level in human primary immune cells.

Measurement of gene expression at the single-cell level has advanced the study of transcriptional regulation programs in healthy and disease states. In particular, single-cell approaches have shed light on the high level of transcriptional heterogeneity of individual cells, both at baseline and in response to experimental or environmental perturbations. We have developed a method for high-content imaging (HCI)-based quantification of relative changes in transcript abundance at the single-cell level in human primary immune cells and have validated its performance under multiple experimental conditions to demonstrate its general applicability. This method, named hcHCR, combines the sensitivity of the hybridization chain reaction (HCR) for the visualization of RNA in single cells, with the speed, scalability, and reproducibility of HCI. We first tested eight cell attachment substrates for short-term culture of primary human B cells, T cells, monocytes, or neutrophils. We then miniaturized HCR in 384-well format and documented the ability of the method to detect changes in transcript abundance at the single-cell level in thousands of cells for each experimental condition by HCI. Furthermore, we demonstrated the feasibility of multiplexing gene expression measurements by simultaneously assaying the abundance of three transcripts per cell at baseline and in response to an experimental stimulus. Finally, we tested the robustness of the assay to technical and biological variation. We anticipate that hcHCR will be suitable for low- to medium-throughput chemical or functional genomics screens in primary human cells, with the possibility of performing screens on cells obtained from patients with a specific disease.

Lipid regulation of NLRP3 inflammasome activity through organelle stress.

Inflammation driven by the NLRP3 inflammasome in macrophages is an important contributor to chronic metabolic diseases that affect growing numbers of individuals. Many of these diseases involve the pathologic accumulation of endogenous lipids or their oxidation products, which can activate NLRP3. Other endogenous lipids, however, can inhibit the activation of NLRP3. The intracellular mechanisms by which these lipids modulate NLRP3 activity are now being identified. This review discusses emerging evidence suggesting that organelle stress, particularly involving mitochondria, lysosomes, and the endoplasmic reticulum, may be key in lipid-induced modification of NLRP3 inflammasome activity.

Type I IFNs facilitate innate immune control of the opportunistic bacteria Burkholderia cenocepacia in the macrophage cytosol.

The mammalian immune system is constantly challenged by signals from both pathogenic and non-pathogenic microbes. Many of these non-pathogenic microbes have pathogenic potential if the immune system is compromised. The importance of type I interferons (IFNs) in orchestrating innate immune responses to pathogenic microbes has become clear in recent years. However, the control of opportunistic pathogens-and especially intracellular bacteria-by type I IFNs remains less appreciated. In this study, we use the opportunistic, Gram-negative bacterial pathogen Burkholderia cenocepacia (Bc) to show that type I IFNs are capable of limiting bacterial replication in macrophages, preventing illness in immunocompetent mice. Sustained type I IFN signaling through cytosolic receptors allows for increased expression of autophagy and linear ubiquitination mediators, which slows bacterial replication. Transcriptomic analyses and in vivo studies also show that LPS stimulation does not replicate the conditions of intracellular Gram-negative bacterial infection as it pertains to type I IFN stimulation or signaling. This study highlights the importance of type I IFNs in protection against opportunistic pathogens through innate immunity, without the need for damaging inflammatory responses.

Correction: Type I IFNs facilitate innate immune control of the opportunistic bacteria Burkholderia cenocepacia in the macrophage cytosol.

[This corrects the article DOI: 10.1371/journal.ppat.1009395.].

Age influences susceptibility of brain capillary endothelial cells to La Crosse virus infection and cell death.

BACKGROUND: A key factor in the development of viral encephalitis is a virus crossing the blood-brain barrier (BBB). We have previously shown that age-related susceptibility of mice to the La Crosse virus (LACV), the leading cause of pediatric arbovirus encephalitis in the USA, was associated with the ability of the virus to cross the BBB. LACV infection in weanling mice (aged around 3 weeks) results in vascular leakage in the olfactory bulb/tract (OB/OT) region of the brain, which is not observed in adult mice aged > 6-8 weeks. Thus, we studied age-specific differences in the response of brain capillary endothelial cells (BCECs) to LACV infection. METHODS: To examine mechanisms of LACV-induced BBB breakdown and infection of the CNS, we analyzed BCECs directly isolated from weanling and adult mice as well as established a model where these cells were infected in vitro and cultured for a short period to determine susceptibility to virus infection and cell death. Additionally, we utilized correlative light electron microscopy (CLEM) to examine whether changes in cell morphology and function were also observed in BCECs in vivo. RESULTS: BCECs from weanling, but not adult mice, had detectable infection after several days in culture when taken ex vivo from infected mice suggesting that these cells could be infected in vitro. Further analysis of BCECs from uninfected mice, infected in vitro, showed that weanling BCECs were more susceptible to virus infection than adult BCECs, with higher levels of infected cells, released virus as well as cytopathic effects (CPE) and cell death. Although direct LACV infection is not detected in the weanling BCECs, CLEM analysis of brain tissue from weanling mice indicated that LACV infection induced significant cerebrovascular damage which allowed virus-sized particles to enter the brain parenchyma. CONCLUSIONS: These findings indicate that BCECs isolated from adult and weanling mice have differential viral load, infectivity, and susceptibility to LACV. These age-related differences in susceptibility may strongly influence LACV-induced BBB leakage and neurovascular damage allowing virus invasion of the CNS and the development of neurological disease.

Navigating the network: signaling cross-talk in hematopoietic cells.

Recent studies in hematopoietic cells have led to a growing appreciation of the diverse modes of molecular and functional cross-talk between canonical signaling pathways. However, these intersections represent only the tip of the iceberg. Emerging global analytical methods are providing an even richer and more complete picture of the many components that measurably interact in a network manner to produce cellular responses. Here we highlight the pieces in this Focus, emphasize the limitations of the present canonical pathway paradigm, and discuss the value of a systems biology approach using more global, quantitative experimental design and data analysis strategies. Lastly, we urge caution about overly facile interpretation of genome- and proteome-level studies.

Dynamic ADP-Ribosylome, Phosphoproteome, and Interactome in LPS-Activated Macrophages.

We have used mass spectrometry (MS) to characterize protein signaling in lipopolysaccharide (LPS)-stimulated macrophages from human blood, human THP1 cells, mouse bone marrow, and mouse Raw264.7 cells. Protein ADP-ribosylation was truncated down to phosphoribose, allowing for enrichment and identification of the resulting phosphoribosylated peptides alongside phosphopeptides. Size exclusion chromatography-MS (SEC-MS) was used to separate proteoforms by size; protein complexes were then identified by weighted correlation network analysis (WGCNA) based on their correlated movement into or out of SEC fractions following stimulation, presenting an analysis method for SEC-MS that does not rely on established databases. We highlight two modules of interest: one linked to the apoptosis signal-regulating kinase (ASK) signalosome and the other containing poly(ADP-ribose) polymerase 9 (PARP9). Finally, PARP inhibition was used to perturb the characterized systems, demonstrating the importance of ADP-ribosylation for the global interactome. All post-translational modification (PTM) and interactome data have been aggregated into a meta-database of 6729 proteins, with ADP-ribosylation characterized on 2905 proteins and phosphorylation characterized on 2669 proteins. This database-titled MAPCD, for Macrophage ADP-ribosylation, Phosphorylation, and Complex Dynamics-serves as an invaluable resource for studying crosstalk between the ADP-ribosylome, phosphoproteome, and interactome.

A Deep Learning Pipeline for Nucleus Segmentation.

Deep learning is rapidly becoming the technique of choice for automated segmentation of nuclei in biological image analysis workflows. In order to evaluate the feasibility of training nuclear segmentation models on small, custom annotated image datasets that have been augmented, we have designed a computational pipeline to systematically compare different nuclear segmentation model architectures and model training strategies. Using this approach, we demonstrate that transfer learning and tuning of training parameters, such as the composition, size, and preprocessing of the training image dataset, can lead to robust nuclear segmentation models, which match, and often exceed, the performance of existing, off-the-shelf deep learning models pretrained on large image datasets. We envision a practical scenario where deep learning nuclear segmentation models trained in this way can be shared across a laboratory, facility, or institution, and continuously improved by training them on progressively larger and varied image datasets. Our work provides computational tools and a practical framework for deep learning-based biological image segmentation using small annotated image datasets. Published [2020]. This article is a U.S. Government work and is in the public domain in the USA.

Dual Roles for Ikaros in Regulation of Macrophage Chromatin State and Inflammatory Gene Expression.

Macrophage activation by bacterial LPS leads to induction of a complex inflammatory gene program dependent on numerous transcription factor families. The transcription factor Ikaros has been shown to play a critical role in lymphoid cell development and differentiation; however, its function in myeloid cells and innate immune responses is less appreciated. Using comprehensive genomic analysis of Ikaros-dependent transcription, DNA binding, and chromatin accessibility, we describe unexpected dual repressor and activator functions for Ikaros in the LPS response of murine macrophages. Consistent with the described function of Ikaros as transcriptional repressor, Ikzf1-/- macrophages showed enhanced induction for select responses. In contrast, we observed a dramatic defect in expression of many delayed LPS response genes, and chromatin immunoprecipitation sequencing analyses support a key role for Ikaros in sustained NF-κB chromatin binding. Decreased Ikaros expression in Ikzf1+/- mice and human cells dampens these Ikaros-enhanced inflammatory responses, highlighting the importance of quantitative control of Ikaros protein level for its activator function. In the absence of Ikaros, a constitutively open chromatin state was coincident with dysregulation of LPS-induced chromatin remodeling, gene expression, and cytokine responses. Together, our data suggest a central role for Ikaros in coordinating the complex macrophage transcriptional program in response to pathogen challenge.

Pharmacokinetic Interactions between the Hepatitis C Virus Inhibitors Elbasvir and Grazoprevir and HIV Protease Inhibitors Ritonavir, Atazanavir, Lopinavir, and Darunavir in Healthy Volunteers.

The combination of the hepatitis C virus (HCV) nonstructural protein 5A (NS5A) inhibitor elbasvir and the NS3/4A protease inhibitor grazoprevir is a potent, once-daily therapy indicated for the treatment of chronic HCV infection in individuals coinfected with human immunodeficiency virus (HIV). We explored the pharmacokinetic interactions of elbasvir and grazoprevir with ritonavir and ritonavir-boosted HIV protease inhibitors in three phase 1 trials. Drug-drug interaction trials with healthy participants were conducted to evaluate the effect of ritonavir on the pharmacokinetics of grazoprevir (n = 10) and the potential two-way pharmacokinetic interactions of elbasvir (n = 30) or grazoprevir (n = 39) when coadministered with ritonavir-boosted atazanavir, lopinavir, or darunavir. Coadministration of ritonavir with grazoprevir increased grazoprevir exposure; the geometric mean ratio (GMR) for grazoprevir plus ritonavir versus grazoprevir alone area under the concentration-time curve from 0 to 24 h (AUC0-24) was 1.91 (90% confidence interval [CI]; 1.31 to 2.79). Grazoprevir exposure was markedly increased with coadministration of atazanavir-ritonavir, lopinavir-ritonavir, and darunavir-ritonavir, with GMRs for grazoprevir AUC0-24 of 10.58 (90% CI, 7.78 to 14.39), 12.86 (90% CI, 10.25 to 16.13), and 7.50 (90% CI, 5.92 to 9.51), respectively. Elbasvir exposure was increased with coadministration of atazanavir-ritonavir, lopinavir-ritonavir, and darunavir-ritonavir, with GMRs for elbasvir AUC0-24 of 4.76 (90% CI, 4.07 to 5.56), 3.71 (90% CI, 3.05 to 4.53), and 1.66 (90% CI, 1.35 to 2.05), respectively. Grazoprevir and elbasvir had little effect on atazanavir, lopinavir, and darunavir pharmacokinetics. Coadministration of elbasvir-grazoprevir with atazanavir-ritonavir, lopinavir-ritonavir, or darunavir-ritonavir is contraindicated, owing to an increase in grazoprevir exposure. Therefore, HIV treatment regimens without HIV protease inhibitors should be considered for HCV/HIV-coinfected individuals who are being treated with elbasvir-grazoprevir.

Assessment of drug interaction potential between the HCV direct-acting antiviral agents elbasvir/grazoprevir and the HIV integrase inhibitors raltegravir and dolutegravir.

BACKGROUND: Elbasvir/grazoprevir is a once-daily fixed-dose combination therapy for the treatment of chronic HCV infection, including HCV/HIV coinfection. OBJECTIVES: To evaluate the pharmacokinetic interaction of elbasvir and grazoprevir with raltegravir or dolutegravir. METHODS: Three open-label trials in healthy adult participants were conducted. In the raltegravir trials, participants received a single dose of raltegravir 400 mg, a single dose of elbasvir 50 mg or grazoprevir 200 mg, and raltegravir with either elbasvir or grazoprevir. In the dolutegravir trial, participants received a single dose of dolutegravir 50 mg alone or co-administered with once-daily elbasvir 50 mg and grazoprevir 200 mg. RESULTS: The raltegravir AUC0-∞ geometric mean ratio (GMR) (90% CI) was 1.02 (0.81-1.27) with elbasvir and 1.43 (0.89-2.30) with grazoprevir. Dolutegravir AUC0-∞ GMR (90% CI) was 1.16 (1.00-1.34) with elbasvir and grazoprevir. The elbasvir AUC0-∞ GMR (90% CI) was 0.81 (0.57-1.17) with raltegravir and 0.98 (0.93-1.04) with dolutegravir. The grazoprevir AUC0-24 GMR (90% CI) was 0.89 (0.72-1.09) with raltegravir and 0.81 (0.67-0.97) with dolutegravir. CONCLUSIONS: Elbasvir or grazoprevir co-administered with raltegravir or dolutegravir resulted in no clinically meaningful drug-drug interactions and was generally well tolerated. These results support the assertion that no dose adjustments for elbasvir, grazoprevir, raltegravir or dolutegravir are needed for co-administration in HCV/HIV-coinfected people.

A randomized clinical trial to evaluate the single-dose pharmacokinetics, pharmacodynamics, and safety of sitagliptin in pediatric patients with type 2 diabetes.

OBJECTIVE: To evaluate the single-dose pharmacokinetics (PK), pharmacodynamics (PD), and safety of sitagliptin in pediatric patients with type 2 diabetes mellitus (T2DM). STUDY DESIGN: This was a randomized, placebo-controlled, double-blind evaluation of sitagliptin in 35 patients 10 to 17 years old with T2DM at 7 clinical research sites. The safety, tolerability, PK, and PD (dipeptidyl peptidase-4 [DPP-4] inhibition and aspects of glucose metabolism) of single doses of 50, 100, and 200 mg were assessed. Appropriate transformations on the PK parameters were used and back-transformed summary statistics are reported. RESULTS: Adverse experiences were reported by eight study participants; all were of mild intensity except one (intravenous site pain of moderate intensity). PK characteristics in the young patients were comparable to reference adult data, with geometric mean ratios (youths/adults) for AUC0-∞ , Cmax , and C24hr of 0.82, 1.04, and 0.74, respectively. Single doses of 50, 100, and 200 mg sitagliptin inhibited 67.2%, 73.8%, and 81.2% of plasma DPP-4 activity over 24 hours, respectively. Least squares (LS) mean glucose concentrations 2 hours after an oral glucose tolerance test or a meal tolerance test decreased in study participants treated with sitagliptin, compared to placebo, while active LS mean glucagon-like peptide 1 concentrations increased significantly at all sitagliptin doses in both tests. CONCLUSIONS: Single doses of sitagliptin as high as 200 mg were generally well tolerated in 10- to 17-year-old male and female study participants with T2DM, and a daily sitagliptin dose of 100 mg is appropriate for evaluation in Phase III safety and efficacy studies in pediatric patients with T2DM. (ClinicalTrials.gov: NCT00730275).

Proteome and Secretome Analysis Reveals Differential Post-transcriptional Regulation of Toll-like Receptor Responses.

The innate immune system is the organism's first line of defense against pathogens. Pattern recognition receptors (PRRs) are responsible for sensing the presence of pathogen-associated molecules. The prototypic PRRs, the membrane-bound receptors of the Toll-like receptor (TLR) family, recognize pathogen-associated molecular patterns (PAMPs) and initiate an innate immune response through signaling pathways that depend on the adaptor molecules MyD88 and TRIF. Deciphering the differences in the complex signaling events that lead to pathogen recognition and initiation of the correct response remains challenging. Here we report the discovery of temporal changes in the protein signaling components involved in innate immunity. Using an integrated strategy combining unbiased proteomics, transcriptomics and macrophage stimulations with three different PAMPs, we identified differences in signaling between individual TLRs and revealed specifics of pathway regulation at the protein level.

The cAMP Pathway Amplifies Early MyD88-Dependent and Type I Interferon-Independent LPS-Induced Interleukin-10 Expression in Mouse Macrophages.

Interleukin-10 (IL-10) is a key anti-inflammatory cytokine, secreted by macrophages and other immune cells to attenuate inflammation. Autocrine type I interferons (IFNs) largely mediate the delayed expression of IL-10 by LPS-stimulated macrophages. We have previously shown that IL-10 is synergistically expressed in macrophages following a costimulus of a TLR agonist and cAMP. We now show that the cAMP pathway directly upregulates IL-10 transcription and plays an important permissive and synergistic role in early, but not late, LPS-stimulated IL-10 mRNA and protein expression in mouse macrophages and in a mouse septic shock model. Our results suggest that the loss of synergism is not due to desensitization of the cAMP inducing signal, and it is not mediated by a positive crosstalk between the cAMP and type I IFN pathways. First, cAMP elevation in LPS-treated cells decreased the secretion of type I IFN. Second, autocrine/paracrine type I IFNs induce IL-10 promoter reporter activity only additively, but not synergistically, with the cAMP pathway. IL-10 promoter reporter activity was synergistically induced by cAMP elevation in macrophages stimulated by an agonist of either TLR4, TLR2/6, or TLR7, receptors which signal via MyD88, but not by an agonist of TLR3 which signals independently of MyD88. Moreover, MyD88 knockout largely reduced the synergistic IL-10 expression, indicating that MyD88 is required for the synergism displayed by LPS with cAMP. This report delineates the temporal regulation of early cAMP-accelerated vs. late type I IFN-dependent IL-10 transcription in LPS-stimulated murine macrophages that can limit inflammation at its onset.

Monoacylglycerol Lipase Inhibition in Human and Rodent Systems Supports Clinical Evaluation of Endocannabinoid Modulators.

Monoacylglycerol lipase (MGLL) is the primary degradative enzyme for the endocannabinoid 2-arachidonoylglycerol (2-AG). The first MGLL inhibitors have recently entered clinical development for the treatment of neurologic disorders. To support this clinical path, we report the pharmacological characterization of the highly potent and selective MGLL inhibitor ABD-1970 [1,1,1,3,3,3-hexafluoropropan-2-yl 4-(2-(8-oxa-3-azabicyclo[3.2.1]octan-3-yl)-4-chlorobenzyl)piperazine-1-carboxylate]. We used ABD-1970 to confirm the role of MGLL in human systems and to define the relationship between MGLL target engagement, brain 2-AG concentrations, and efficacy. Because MGLL contributes to arachidonic acid metabolism in a subset of rodent tissues, we further used ABD-1970 to evaluate whether selective MGLL inhibition would affect prostanoid production in several human assays known to be sensitive to cyclooxygenase inhibitors. ABD-1970 robustly elevated brain 2-AG content and displayed antinociceptive and antipruritic activity in a battery of rodent models (ED50 values of 1-2 mg/kg). The antinociceptive effects of ABD-1970 were potentiated when combined with analgesic standards of care and occurred without overt cannabimimetic effects. ABD-1970 also blocked 2-AG hydrolysis in human brain tissue and elevated 2-AG content in human blood without affecting stimulated prostanoid production. These findings support the clinical development of MGLL inhibitors as a differentiated mechanism to treat pain and other neurologic disorders.