Spatially resolved in silico modeling of NKG2D signaling kinetics suggests a key role of NKG2D and Vav1 Co-clustering in generating natural killer cell activation

Abstract

Natural Killer (NK) cells provide key resistance against viral infections and tumors. A diverse set of activating and inhibitory NK cell receptors (NKRs) interact with cognate ligands presented by target host cells, where integration of dueling signals initiated by the ligand-NKR interactions determines NK cell activation or tolerance. Imaging experiments over decades have shown micron and sub-micron scale spatial clustering of activating and inhibitory NKRs. The mechanistic roles of these clusters in affecting downstream signaling and activation are often unclear. To this end, we developed a preDISCLOSUREdictive in silico framework by combining spatially resolved mechanistic agent based modeling, published TIRF imaging data, and parameter estimation to determine mechanisms by which formation and spatial movements of activating NKG2D microclusters affect early time NKG2D signaling kinetics in a human cell line NKL. We show co-clustering of NKG2D and the guanosine nucleotide exchange factor Vav1 in NKG2D microclusters plays a dominant role over ligand (ULBP3) rebinding in increasing production of phospho-Vav1(pVav1), an activation marker of early NKG2D signaling. The in silico model successfully predicts several scenarios of inhibition of NKG2D signaling and time course of NKG2D spatial clustering over a short (~3 min) interval. Modeling shows the presence of a spatial positive feedback relating formation and centripetal movements of NKG2D microclusters, and pVav1 production offers flexibility towards suppression of activating signals by inhibitory KIR ligands organized in inhomogeneous spatial patterns (e.g., a ring). Our in silico framework marks a major improvement in developing spatiotemporal signaling models with quantitatively estimated model parameters using imaging data.

PASCAR: a multiscale framework to explore the design space of constitutive and inducible CAR T cells

“…This could be a reason why our model underpredicted percentage lysis for the constitutive CAR T cells ( Fig 2G ) at low and high CAR abundances where PASCAR used the same signaling rates estimated for intermediate CAR abundances. A potential extension of PASCAR could be to include such details of signaling and spatial clustering that can be developed by building on several existing CAR T cell ( Rohrs et al, 2018 ) and signaling models in lymphocytes ( Čemerski et al, 2008 ; Grewal & Das, 2022 ).…”

BioNetGMMFit: estimating parameters of a BioNetGen model from time-stamped snapshots of single cells

“…BioNetGMMFit is currently limited to parameter estimation for well-mixed mechanistic models describing deterministic and stochastic kinetics. Unlike PyBioNetFit 18 , other forms of simulation such as NFsim 31 and spatial modeling 32 , 33 are not supported. Furthermore, BioNetGMMFit does not support any optimization algorithms other than PSO.…”

Killer instincts: natural killer cells as multifactorial cancer immunotherapy

“…A subset of KIRs, the Fc receptor CD16a, and the natural cytotoxicity receptors (NCRs) convey activating signals via immunoreceptor tyrosine activating motifs (ITAMs), which activate SFKs, and other signaling intermediaries including ZAP-70 and Vav-1, then phospholipase C, PI3K, Rho-family GTPases ( 58 , 59 ). Additional activating receptors include the NKG2 family members NKG2C and NKG2D, which signal via DAP12 and DAP10, respectively, and therefore bypass the need for SFKs, instead shunting directly to activating downstream mediators including Vav-1, PI3K, Rho-family GTPases and phospholipase C ( 60 – 63 ). Immunoglobulin tail tyrosine motifs (ITTs), used by DNAM-1, are phosphorylated by SFKs and similarly drive downstream activation ( 64 , 65 ).…”