How persistent infection overcomes peripheral tolerance mechanisms to cause T cell-mediated autoimmune disease.

T cells help orchestrate immune responses to pathogens, and their aberrant regulation can trigger autoimmunity. Recent studies highlight that a threshold number of T cells (a quorum) must be activated in a tissue to mount a functional immune response. These collective effects allow the T cell repertoire to respond to pathogens while suppressing autoimmunity due to circulating autoreactive T cells. Our computational studies show that increasing numbers of pathogenic peptides targeted by T cells during persistent or severe viral infections increase the probability of activating T cells that are weakly reactive to self-antigens (molecular mimicry). These T cells are easily re-activated by the self-antigens and contribute to exceeding the quorum threshold required to mount autoimmune responses. Rare peptides that activate many T cells are sampled more readily during severe/persistent infections than in acute infections, which amplifies these effects. Experiments in mice to test predictions from these mechanistic insights are suggested.

Polymer folding through active processes recreates features of genome organization.

From proteins to chromosomes, polymers fold into specific conformations that control their biological function. Polymer folding has long been studied with equilibrium thermodynamics, yet intracellular organization and regulation involve energy-consuming, active processes. Signatures of activity have been measured in the context of chromatin motion, which shows spatial correlations and enhanced subdiffusion only in the presence of adenosine triphosphate. Moreover, chromatin motion varies with genomic coordinate, pointing toward a heterogeneous pattern of active processes along the sequence. How do such patterns of activity affect the conformation of a polymer such as chromatin? We address this question by combining analytical theory and simulations to study a polymer subjected to sequence-dependent correlated active forces. Our analysis shows that a local increase in activity (larger active forces) can cause the polymer backbone to bend and expand, while less active segments straighten out and condense. Our simulations further predict that modest activity differences can drive compartmentalization of the polymer consistent with the patterns observed in chromosome conformation capture experiments. Moreover, segments of the polymer that show correlated active (sub)diffusion attract each other through effective long-ranged harmonic interactions, whereas anticorrelations lead to effective repulsions. Thus, our theory offers nonequilibrium mechanisms for forming genomic compartments, which cannot be distinguished from affinity-based folding using structural data alone. As a first step toward exploring whether active mechanisms contribute to shaping genome conformations, we discuss a data-driven approach.

Adenovirus-vectored vaccine containing multidimensionally conserved parts of the HIV proteome is immunogenic in rhesus macaques.

An effective vaccine that can protect against HIV infection does not exist. A major reason why a vaccine is not available is the high mutability of the virus, which enables it to evolve mutations that can evade human immune responses. This challenge is exacerbated by the ability of the virus to evolve compensatory mutations that can partially restore the fitness cost of immune-evading mutations. Based on the fitness landscapes of HIV proteins that account for the effects of coupled mutations, we designed a single long peptide immunogen comprising parts of the HIV proteome wherein mutations are likely to be deleterious regardless of the sequence of the rest of the viral protein. This immunogen was then stably expressed in adenovirus vectors that are currently in clinical development. Macaques immunized with these vaccine constructs exhibited T-cell responses that were comparable in magnitude to animals immunized with adenovirus vectors with whole HIV protein inserts. Moreover, the T-cell responses in immunized macaques strongly targeted regions contained in our immunogen. These results suggest that further studies aimed toward using our vaccine construct for HIV prophylaxis and cure are warranted.

Deconvolving mutational patterns of poliovirus outbreaks reveals its intrinsic fitness landscape.

Vaccination has essentially eradicated poliovirus. Yet, its mutation rate is higher than that of viruses like HIV, for which no effective vaccine exists. To investigate this, we infer a fitness model for the poliovirus viral protein 1 (vp1), which successfully predicts in vitro fitness measurements. This is achieved by first developing a probabilistic model for the prevalence of vp1 sequences that enables us to isolate and remove data that are subject to strong vaccine-derived biases. The intrinsic fitness constraints derived for vp1, a capsid protein subject to antibody responses, are compared with those of analogous HIV proteins. We find that vp1 evolution is subject to tighter constraints, limiting its ability to evade vaccine-induced immune responses. Our analysis also indicates that circulating poliovirus strains in unimmunized populations serve as a reservoir that can seed outbreaks in spatio-temporally localized sub-optimally immunized populations.

CD45 functions as a signaling gatekeeper in T cells.

T cells require the protein tyrosine phosphatase CD45 to detect and respond to antigen because it activates the Src family kinase Lck, which phosphorylates the T cell antigen receptor (TCR) complex. CD45 activates Lck by opposing the negative regulatory kinase Csk. Paradoxically, CD45 has also been implicated in suppressing TCR signaling by dephosphorylating the same signaling motifs within the TCR complex upon which Lck acts. We sought to reconcile these observations using chemical and genetic perturbations of the Csk/CD45 regulatory axis incorporated with computational analyses. Specifically, we titrated the activities of Csk and CD45 and assessed their influence on Lck activation, TCR-associated ζ-chain phosphorylation, and more downstream signaling events. Acute inhibition of Csk revealed that CD45 suppressed ζ-chain phosphorylation and was necessary for a regulatable pool of active Lck, thereby interconnecting the activating and suppressive roles of CD45 that tune antigen discrimination. CD45 suppressed signaling events that were antigen independent or induced by low-affinity antigen but not those initiated by high-affinity antigen. Together, our findings reveal that CD45 acts as a signaling "gatekeeper," enabling graded signaling outputs while filtering weak or spurious signaling events.

The low spike density of HIV may have evolved because of the effects of T helper cell depletion on affinity maturation.

The spikes on virus surfaces bind receptors on host cells to propagate infection. High spike densities (SDs) can promote infection, but spikes are also targets of antibody-mediated immune responses. Thus, diverse evolutionary pressures can influence virus SDs. HIV's SD is about two orders of magnitude lower than that of other viruses, a surprising feature of unknown origin. By modeling antibody evolution through affinity maturation, we find that an intermediate SD maximizes the affinity of generated antibodies. We argue that this leads most viruses to evolve high SDs. T helper cells, which are depleted during early HIV infection, play a key role in antibody evolution. We find that T helper cell depletion results in high affinity antibodies when SD is high, but not if SD is low. This special feature of HIV infection may have led to the evolution of a low SD to avoid potent immune responses early in infection.

Fitness landscape of the human immunodeficiency virus envelope protein that is targeted by antibodies.

HIV is a highly mutable virus, and over 30 years after its discovery, a vaccine or cure is still not available. The isolation of broadly neutralizing antibodies (bnAbs) from HIV-infected patients has led to renewed hope for a prophylactic vaccine capable of combating the scourge of HIV. A major challenge is the design of immunogens and vaccination protocols that can elicit bnAbs that target regions of the virus's spike proteins where the likelihood of mutational escape is low due to the high fitness cost of mutations. Related challenges include the choice of combinations of bnAbs for therapy. An accurate representation of viral fitness as a function of its protein sequences (a fitness landscape), with explicit accounting of the effects of coupling between mutations, could help address these challenges. We describe a computational approach that has allowed us to infer a fitness landscape for gp160, the HIV polyprotein that comprises the viral spike that is targeted by antibodies. We validate the inferred landscape through comparisons with experimental fitness measurements, and various other metrics. We show that an effective antibody that prevents immune escape must selectively bind to high escape cost residues that are surrounded by those where mutations incur a low fitness cost, motivating future applications of our landscape for immunogen design.

Rational design of vaccine targets and strategies for HIV: a crossroad of statistical physics, biology, and medicine.

Vaccination has saved more lives than any other medical procedure. Pathogens have now evolved that have not succumbed to vaccination using the empirical paradigms pioneered by Pasteur and Jenner. Vaccine design strategies that are based on a mechanistic understanding of the pertinent immunology and virology are required to confront and eliminate these scourges. In this perspective, we describe just a few examples of work aimed to achieve this goal by bringing together approaches from statistical physics with biology and clinical research.

Phosphorylation of a Tyrosine Residue on Zap70 by Lck and Its Subsequent Binding via an SH2 Domain May Be a Key Gatekeeper of T Cell Receptor Signaling In Vivo.

The initiation of signaling in T lymphocytes in response to the binding of the T cell receptor (TCR) to cognate ligands is a key step in the emergence of adaptive immune responses. Conventional models posit that TCR signaling is initiated by the phosphorylation of receptor-associated immune receptor activation motifs (ITAMs). The cytoplasmic tyrosine kinase Zap70 binds to phosphorylated ITAMs, is subsequently activated, and then propagates downstream signaling. While evidence for such models is provided by experiments with cell lines, in vivo, Zap70 is bound to phosphorylated ITAMs in resting T cells. However, Zap70 is activated only upon TCR binding to cognate ligand. We report the results of computational studies of a new model for the initiation of TCR signaling that incorporates these in vivo observations. Importantly, the new model is shown to allow better and faster TCR discrimination between self-ligands and foreign ligands. The new model is consistent with many past experimental observations, and experiments that could further test the model are proposed.

Relative rate and location of intra-host HIV evolution to evade cellular immunity are predictable.

Human immunodeficiency virus (HIV) evolves within infected persons to escape being destroyed by the host immune system, thereby preventing effective immune control of infection. Here, we combine methods from evolutionary dynamics and statistical physics to simulate in vivo HIV sequence evolution, predicting the relative rate of escape and the location of escape mutations in response to T-cell-mediated immune pressure in a cohort of 17 persons with acute HIV infection. Predicted and clinically observed times to escape immune responses agree well, and we show that the mutational pathways to escape depend on the viral sequence background due to epistatic interactions. The ability to predict escape pathways and the duration over which control is maintained by specific immune responses open the door to rational design of immunotherapeutic strategies that might enable long-term control of HIV infection. Our approach enables intra-host evolution of a human pathogen to be predicted in a probabilistic framework.

Identification of drug resistance mutations in HIV from constraints on natural evolution.

Human immunodeficiency virus (HIV) evolves with extraordinary rapidity. However, its evolution is constrained by interactions between mutations in its fitness landscape. Here we show that an Ising model describing these interactions, inferred from sequence data obtained prior to the use of antiretroviral drugs, can be used to identify clinically significant sites of resistance mutations. Successful predictions of the resistance sites indicate progress in the development of successful models of real viral evolution at the single residue level and suggest that our approach may be applied to help design new therapies that are less prone to failure even where resistance data are not yet available.

The catalytic activity of the kinase ZAP-70 mediates basal signaling and negative feedback of the T cell receptor pathway.

T cell activation by antigens binding to the T cell receptor (TCR) must be properly regulated to ensure normal T cell development and effective immune responses to pathogens and transformed cells while avoiding autoimmunity. The Src family kinase Lck and the Syk family kinase ZAP-70 (ζ chain-associated protein kinase of 70 kD) are sequentially activated in response to TCR engagement and serve as critical components of the TCR signaling machinery that leads to T cell activation. We performed a mass spectrometry-based phosphoproteomic study comparing the quantitative differences in the temporal dynamics of phosphorylation in stimulated and unstimulated T cells with or without inhibition of ZAP-70 catalytic activity. The data indicated that the kinase activity of ZAP-70 stimulates negative feedback pathways that target Lck and thereby modulate the phosphorylation patterns of the immunoreceptor tyrosine-based activation motifs (ITAMs) of the CD3 and ζ chain components of the TCR and of signaling molecules downstream of Lck, including ZAP-70. We developed a computational model that provides a mechanistic explanation for the experimental findings on ITAM phosphorylation in wild-type cells, ZAP-70-deficient cells, and cells with inhibited ZAP-70 catalytic activity. This model incorporated negative feedback regulation of Lck activity by the kinase activity of ZAP-70 and predicted the order in which tyrosines in the ITAMs of TCR ζ chains must be phosphorylated to be consistent with the experimental data.

Scaling laws describe memories of host-pathogen riposte in the HIV population.

The enormous genetic diversity and mutability of HIV has prevented effective control of this virus by natural immune responses or vaccination. Evolution of the circulating HIV population has thus occurred in response to diverse, ultimately ineffective, immune selection pressures that randomly change from host to host. We show that the interplay between the diversity of human immune responses and the ways that HIV mutates to evade them results in distinct sets of sequences defined by similar collectively coupled mutations. Scaling laws that relate these sets of sequences resemble those observed in linguistics and other branches of inquiry, and dynamics reminiscent of neural networks are observed. Like neural networks that store memories of past stimulation, the circulating HIV population stores memories of host-pathogen combat won by the virus. We describe an exactly solvable model that captures the main qualitative features of the sets of sequences and a simple mechanistic model for the origin of the observed scaling laws. Our results define collective mutational pathways used by HIV to evade human immune responses, which could guide vaccine design.

Understanding immunology: fun at an intersection of the physical, life, and clinical sciences.

Understanding how the immune system works is a grand challenge in science with myriad direct implications for improving human health. The immune system protects us from infectious pathogens and cancer, and maintains a harmonious steady state with essential microbiota in our gut. Vaccination, the medical procedure that has saved more lives than any other, involves manipulating the immune system. Unfortunately, the immune system can also go awry to cause autoimmune diseases. Immune responses are the product of stochastic collective dynamic processes involving many interacting components. These processes span multiple scales of length and time. Thus, statistical mechanics has much to contribute to immunology, and the oeuvre of biological physics will be further enriched if the number of physical scientists interested in immunology continues to increase. I describe how I got interested in immunology and provide a glimpse of my experiences working on immunology using approaches from statistical mechanics and collaborating closely with immunologists.

Coreceptor scanning by the T cell receptor provides a mechanism for T cell tolerance.

In the thymus, high-affinity, self-reactive thymocytes are eliminated from the pool of developing T cells, generating central tolerance. Here, we investigate how developing T cells measure self-antigen affinity. We show that very few CD4 or CD8 coreceptor molecules are coupled with the signal-initiating kinase, Lck. To initiate signaling, an antigen-engaged T cell receptor (TCR) scans multiple coreceptor molecules to find one that is coupled to Lck; this is the first and rate-limiting step in a kinetic proofreading chain of events that eventually leads to TCR triggering and negative selection. MHCII-restricted TCRs require a shorter antigen dwell time (0.2 s) to initiate negative selection compared to MHCI-restricted TCRs (0.9 s) because more CD4 coreceptors are Lck-loaded compared to CD8. We generated a model (Lck come&stay/signal duration) that accurately predicts the observed differences in antigen dwell-time thresholds used by MHCI- and MHCII-restricted thymocytes to initiate negative selection and generate self-tolerance.

Molecular kinetics. Ras activation by SOS: allosteric regulation by altered fluctuation dynamics.

Activation of the small guanosine triphosphatase H-Ras by the exchange factor Son of Sevenless (SOS) is an important hub for signal transduction. Multiple layers of regulation, through protein and membrane interactions, govern activity of SOS. We characterized the specific activity of individual SOS molecules catalyzing nucleotide exchange in H-Ras. Single-molecule kinetic traces revealed that SOS samples a broad distribution of turnover rates through stochastic fluctuations between distinct, long-lived (more than 100 seconds), functional states. The expected allosteric activation of SOS by Ras-guanosine triphosphate (GTP) was conspicuously absent in the mean rate. However, fluctuations into highly active states were modulated by Ras-GTP. This reveals a mechanism in which functional output may be determined by the dynamical spectrum of rates sampled by a small number of enzymes, rather than the ensemble average.

Statistical linkage analysis of substitutions in patient-derived sequences of genotype 1a hepatitis C virus nonstructural protein 3 exposes targets for immunogen design.

UNLABELLED: Chronic hepatitis C virus (HCV) infection is one of the leading causes of liver failure and liver cancer, affecting around 3% of the world's population. The extreme sequence variability of the virus resulting from error-prone replication has thwarted the discovery of a universal prophylactic vaccine. It is known that vigorous and multispecific cellular immune responses, involving both helper CD4(+) and cytotoxic CD8(+) T cells, are associated with the spontaneous clearance of acute HCV infection. Escape mutations in viral epitopes can, however, abrogate protective T-cell responses, leading to viral persistence and associated pathologies. Despite the propensity of the virus to mutate, there might still exist substitutions that incur a fitness cost. In this paper, we identify groups of coevolving residues within HCV nonstructural protein 3 (NS3) by analyzing diverse sequences of this protein using ideas from random matrix theory and associated methods. Our analyses indicate that one of these groups comprises a large percentage of residues for which HCV appears to resist multiple simultaneous substitutions. Targeting multiple residues in this group through vaccine-induced immune responses should either lead to viral recognition or elicit escape substitutions that compromise viral fitness. Our predictions are supported by published clinical data, which suggested that immune genotypes associated with spontaneous clearance of HCV preferentially recognized and targeted this vulnerable group of residues. Moreover, mapping the sites of this group onto the available protein structure provided insight into its functional significance. An epitope-based immunogen is proposed as an alternative to the NS3 epitopes in the peptide-based vaccine IC41. IMPORTANCE: Despite much experimental work on HCV, a thorough statistical study of the HCV sequences for the purpose of immunogen design was missing in the literature. Such a study is vital to identify epistatic couplings among residues that can provide useful insights for designing a potent vaccine. In this work, ideas from random matrix theory were applied to characterize the statistics of substitutions within the diverse publicly available sequences of the genotype 1a HCV NS3 protein, leading to a group of sites for which HCV appears to resist simultaneous substitutions possibly due to deleterious effect on viral fitness. Our analysis leads to completely novel immunogen designs for HCV. In addition, the NS3 epitopes used in the recently proposed peptide-based vaccine IC41 were analyzed in the context of our framework. Our analysis predicts that alternative NS3 epitopes may be worth exploring as they might be more efficacious.

Coreceptor affinity for MHC defines peptide specificity requirements for TCR interaction with coagonist peptide-MHC.

Recent work has demonstrated that nonstimulatory endogenous peptides can enhance T cell recognition of antigen, but MHCI- and MHCII-restricted systems have generated very different results. MHCII-restricted TCRs need to interact with the nonstimulatory peptide-MHC (pMHC), showing peptide specificity for activation enhancers or coagonists. In contrast, the MHCI-restricted cells studied to date show no such peptide specificity for coagonists, suggesting that CD8 binding to noncognate MHCI is more important. Here we show how this dichotomy can be resolved by varying CD8 and TCR binding to agonist and coagonists coupled with computer simulations, and we identify two distinct mechanisms by which CD8 influences the peptide specificity of coagonism. Mechanism 1 identifies the requirement of CD8 binding to noncognate ligand and suggests a direct relationship between the magnitude of coagonism and CD8 affinity for coagonist pMHCI. Mechanism 2 describes how the affinity of CD8 for agonist pMHCI changes the requirement for specific coagonist peptides. MHCs that bind CD8 strongly were tolerant of all or most peptides as coagonists, but weaker CD8-binding MHCs required stronger TCR binding to coagonist, limiting the potential coagonist peptides. These findings in MHCI systems also explain peptide-specific coagonism in MHCII-restricted cells, as CD4-MHCII interaction is generally weaker than CD8-MHCI.

Activation of extracellular signal-regulated kinase but not of p38 mitogen-activated protein kinase pathways in lymphocytes requires allosteric activation of SOS.

Thymocytes convert graded T cell receptor (TCR) signals into positive selection or deletion, and activation of extracellular signal-related kinase (ERK), p38, and Jun N-terminal protein kinase (JNK) mitogen-activated protein kinases (MAPKs) has been postulated to play a discriminatory role. Two families of Ras guanine nucleotide exchange factors (RasGEFs), SOS and RasGRP, activate Ras and the downstream RAF-MEK-ERK pathway. The pathways leading to lymphocyte p38 and JNK activation are less well defined. We previously described how RasGRP alone induces analog Ras-ERK activation while SOS and RasGRP cooperate to establish bimodal ERK activation. Here we employed computational modeling and biochemical experiments with model cell lines and thymocytes to show that TCR-induced ERK activation grows exponentially in thymocytes and that a W729E allosteric pocket mutant, SOS1, can only reconstitute analog ERK signaling. In agreement with RasGRP allosterically priming SOS, exponential ERK activation is severely decreased by pharmacological or genetic perturbation of the phospholipase Cγ (PLCγ)-diacylglycerol-RasGRP1 pathway. In contrast, p38 activation is not sharply thresholded and requires high-level TCR signal input. Rac and p38 activation depends on SOS1 expression but not allosteric activation. Based on computational predictions and experiments exploring whether SOS functions as a RacGEF or adaptor in Rac-p38 activation, we established that the presence of SOS1, but not its enzymatic activity, is critical for p38 activation.