Abstract: Effective exploration and analysis tools are vital for the extraction of insights from single-cell data. However, current techniques for modeling single-cell studies performed across experimental conditions (e.g. samples) require restrictive assumptions or do not adequately deconvolute condition-to-condition variation from cell-to-cell variation. Here, we report that Reduction and Insight in Single-cell Exploration (RISE), an adaptation of the tensor decomposition method PARAFAC2, enables the dimensionality reduction and analysis of single-cell data across conditions. We demonstrate the benefits of RISE across distinct examples of single-cell RNA-sequencing experiments of peripheral immune cells: pharmacologic drug perturbations and systemic lupus erythematosus patient samples. RISE enables associations of gene variation patterns with patients or perturbations, while connecting each coordinated change to single cells without requiring cell type annotations. The theoretical grounding of RISE suggests a unified framework for many single-cell data modeling tasks, while providing an intuitive dimensionality reduction approach for multi-sample single-cell studies across biological contexts.

Abstract: Systems serology aims to broadly profile the antigen binding, Fc biophysical features, immune receptor engagement, and effector functions of antibodies. This experimental approach excels at identifying antibody functional features that are relevant to a particular disease. However, a crucial limitation of this approach is its incomplete description of what structural features of the antibodies are responsible for the observed immune receptor engagement and effector functions. Knowing these antibody features is important for both understanding how effector responses are naturally controlled through antibody Fc structure and designing antibody therapies with specific effector profiles. Here, we address this limitation by modeling the molecular interactions occurring in these assays and using this model to infer quantities of specific antibody Fc species among the antibodies being profiled. We used several validation strategies to show that the model accurately infers antibody properties and then applied the model to infer previously unavailable antibody fucosylation information from existing systems serology data. Using this capability, we find that COVID-19 vaccine efficacy is associated with the induction of afucosylated spike protein-targeting IgG. Our results also question an existing assumption that controllers of HIV exhibit gp120-targeting IgG that are less fucosylated than those of progressors. Additionally, we confirm that afucosylated IgG is associated with membrane-associated antigens for COVID-19 and HIV, and present new evidence indicating that this relationship is specific to the host cell membrane. Finally, we use the model to identify redundant assay measurements and subsets of information-rich measurements from which Fc properties can be inferred. In total, our modeling approach provides a quantitative framework for the reasoning typically applied in these studies, improving the ability to draw mechanistic conclusions from these data.

Abstract: Cytokines mediate cell-to-cell communication across the immune system and therefore are critical to immunosurveillance in cancer and other diseases. Several cytokines show dysregulated abundance or signaling responses in breast cancer, associated with the disease and differences in survival and progression. Cytokines operate in a coordinated manner to affect immune surveillance and regulate one another, necessitating a systems approach for a complete picture of this dysregulation. Here, we profiled cytokine signaling responses of peripheral immune cells from breast cancer patients as compared to healthy controls in a multidimensional manner across ligands, cell populations, and responsive pathways. We find alterations in cytokine responsiveness across pathways and cell types that are best defined by integrated signatures across dimensions. Alterations in the abundance of a cytokine’s cognate receptor do not explain differences in responsiveness. Rather, alterations in baseline signaling and receptor abundance suggesting immune cell reprogramming are associated with altered responses. These integrated features suggest a global reprogramming of immune cell communication in breast cancer.

Abstract: Recent biological studies have been revolutionized in scale and granularity by multiplex and high-throughput assays. Profiling cell responses across several experimental parameters, such as perturbations, time, and genetic contexts, leads to richer and more generalizable findings. However, these multidimensional datasets necessitate a reevaluation of the conventional methods for their representation and analysis. Traditionally, experimental parameters are merged to flatten the data into a two-dimensional matrix, sacrificing crucial experiment context reflected by the structure. As Marshall McLuhan famously stated, “the medium is the message.” In this work, we propose that the experiment structure is the medium in which subsequent analysis is performed, and the optimal choice of data representation must reflect the experiment structure. We review how tensor-structured analyses and decompositions can preserve this information. We contend that tensor methods are poised to become integral to the biomedical data sciences toolkit.

Abstract: The cytokine interleukin-2 (IL-2) has the potential to treat autoimmune disease but is limited by its modest specificity toward immunosuppressive regulatory T (Treg) cells. IL-2 receptors consist of combinations of α, β, and γ chains of variable affinity and cell specificity. Engineering IL-2 to treat autoimmunity has primarily focused on retaining binding to the relatively Treg-selective, high-affinity receptor while reducing binding to the less selective, low-affinity receptor. However, we found that refining the designs to focus on targeting the high-affinity receptor through avidity effects is key to optimizing Treg selectivity. We profiled the dynamics and dose dependency of signaling responses in primary human immune cells induced by engineered fusions composed of either wild-type IL-2 or mutant forms with altered affinity, valency, and fusion to the antibody Fc region for stability. Treg selectivity and signaling response variations were explained by a model of multivalent binding and dimer-enhanced avidity—a combined measure of the strength, number, and conformation of interaction sites—from which we designed tetravalent IL-2–Fc fusions that had greater Treg selectivity in culture than do current designs. Biasing avidity toward IL2Rα with an asymmetrical multivalent design consisting of one α/β chain–binding and one α chain–binding mutant further enhanced Treg selectivity. Comparative analysis revealed that IL2Rα was the optimal cell surface target for Treg selectivity, indicating that avidity for IL2Rα may be the optimal route to producing IL-2 variants that selectively target Tregs.

Abstract: Multivalent cell surface receptor binding is a ubiquitous biological phenomenon with functional and therapeutic significance. Predicting the amount of ligand binding for a cell remains an important question in computational biology as it can provide great insight into cell-to-cell communication and rational drug design toward specific targets. In this study, we extend a mechanistic, two-step multivalent binding model. This model predicts the behavior of a mixture of different multivalent ligand complexes binding to cells expressing various types of receptors. It accounts for the combinatorially large number of interactions between multiple ligands and receptors, optionally allowing a mixture of complexes with different valencies and complexes that contain heterogeneous ligand units. We derive the macroscopic predictions and demonstrate how this model enables large-scale predictions on mixture binding and the binding space of a ligand. This model thus provides an elegant and computationally efficient framework for analyzing multivalent binding.