Thymic selection of T-cell receptors as an extreme value problem

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.

Quantitative approaches for decoding the specificity of the human T cell repertoire

T cell receptor (TCR)-peptide-major histocompatibility complex (pMHC) interactions play a vital role in initiating immune responses against pathogens, and the specificity of TCRpMHC interactions is crucial for developing optimized therapeutic strategies. The advent of high-throughput immunological and structural evaluation of TCR and pMHC has provided an abundance of data for computational approaches that aim to predict favorable TCR-pMHC interactions. Current models are constructed using information on protein sequence, structures, or a combination of both, and utilize a variety of statistical learning-based approaches for identifying the rules governing specificity. This review examines the current theoretical, computational, and deep learning approaches for identifying TCR-pMHC recognition pairs, placing emphasis on each method's mathematical approach, predictive performance, and limitations.

Quantifying changes in the T cell receptor repertoire during thymic development

One of the feats of adaptive immunity is its ability to recognize foreign pathogens while sparing the self. During maturation in the thymus, T cells are selected through the binding properties of their antigen-specific T-cell receptor (TCR), through the elimination of both weakly (positive selection) and strongly (negative selection) self-reactive receptors. However, the impact of thymic selection on the TCR repertoire is poorly understood. Here, we use transgenic Nur77-mice expressing a T-cell activation reporter to study the repertoires of thymic T cells at various stages of their development, including cells that do not pass selection. We combine high-throughput repertoire sequencing with statistical inference techniques to characterize the selection of the TCR in these distinct subsets. We find small but significant differences in the TCR repertoire parameters between the maturation stages, which recapitulate known differentiation pathways leading to the CD4⁺ and CD8⁺ subtypes. These differences can be simulated by simple models of selection acting linearly on the sequence features. We find no evidence of specific sequences or sequence motifs or features that are suppressed by negative selection. These results favour a collective or statistical model for T-cell self non-self discrimination, where negative selection biases the repertoire away from self recognition, rather than ensuring lack of self-reactivity at the single-cell level.

Contact map dependence of a T-cell receptor binding repertoire

The T-cell arm of the adaptive immune system provides the host protection against unknown pathogens by discriminating between host and foreign material. This discriminatory capability is achieved by the creation of a repertoire of cells each carrying a T-cell receptor (TCR) specific to non-self-antigens displayed as peptides bound to the major histocompatibility complex (pMHC). The understanding of the dynamics of the adaptive immune system at a repertoire level is complex, due to both the nuanced interaction of a TCR-pMHC pair and to the number of different possible TCR-pMHC pairings, making computationally exact solutions currently unfeasible. To gain some insight into this problem, we study an affinity-based model for TCR-pMHC binding in which a crystal structure is used to generate a distance-based contact map that weights the pairwise amino acid interactions. We find that the TCR-pMHC binding energy distribution strongly depends both on the number of contacts and the repeat structure allowed by the topology of the contact map of choice; this in turn influences T-cell recognition probability during negative selection, with higher variances leading to higher survival probabilities. In addition, we quantify the degree to which neoantigens with mutations in sites with higher contacts are recognized at a higher rate.

Tregs self-organize into a computing ecosystem and implement a sophisticated optimization algorithm for mediating immune response

Regulatory T cells (Tregs) play a crucial role in mediating immune response. Yet an algorithmic understanding of the role of Tregs in adaptive immunity remains lacking. Here, we present a biophysically realistic model of Treg-mediated self-tolerance in which Tregs bind to self-antigens and locally inhibit the proliferation of nearby activated T cells. By exploiting a duality between ecological dynamics and constrained optimization, we show that Tregs tile the potential antigen space while simultaneously minimizing the overlap between Treg activation profiles. We find that for sufficiently high Treg diversity, Treg-mediated self-tolerance is robust to fluctuations in self-antigen concentrations but lowering the Treg diversity results in a sharp transition-related to the Gardner transition in perceptrons-to a regime where changes in self-antigen concentrations can result in an autoimmune response. We propose an experimental test of this transition in immune-deficient mice and discuss potential implications for autoimmune diseases.

Learning from HIV-1 to predict the immunogenicity of T cell epitopes in SARS-CoV-2

We describe a physics-based learning model for predicting the immunogenicity of cytotoxic T lymphocyte (CTL) epitopes derived from diverse pathogens including SARS-CoV-2. The model was trained and optimized on the relative immunodominance of CTL epitopes in human immunodeficiency virus infection. Its accuracy was tested against experimental data from patients with COVID-19. Our model predicts that only some SARS-CoV-2 epitopes predicted to bind to HLA molecules are immunogenic. The immunogenic CTL epitopes across all SARS-CoV-2 proteins are predicted to provide broad population coverage, but those from the SARS-CoV-2 spike protein alone are unlikely to do so. Our model also predicts that several immunogenic SARS-CoV-2 CTL epitopes are identical to seasonal coronaviruses circulating in the population and such cross-reactive CD8+ T cells can indeed be detected in prepandemic blood donors, suggesting that some level of CTL immunity against COVID-19 may be present in some individuals before SARS-CoV-2 infection.

In silico T cell epitope identification for SARS-CoV-2: Progress and perspectives

Growing evidence suggests that T cells may play a critical role in combating severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Hence, COVID-19 vaccines that can elicit a robust T cell response may be particularly important. The design, development and experimental evaluation of such vaccines is aided by an understanding of the landscape of T cell epitopes of SARS-CoV-2, which is largely unknown. Due to the challenges of identifying epitopes experimentally, many studies have proposed the use of in silico methods. Here, we present a review of the in silico methods that have been used for the prediction of SARS-CoV-2 T cell epitopes. These methods employ a diverse set of technical approaches, often rooted in machine learning. A performance comparison is provided based on the ability to identify a specific set of immunogenic epitopes that have been determined experimentally to be targeted by T cells in convalescent COVID-19 patients, shedding light on the relative performance merits of the different approaches adopted by the in silico studies. The review also puts forward perspectives for future research directions.

Predicting the Immunogenicity of T cell epitopes: From HIV to SARS-CoV-2

We describe a physics-based learning model for predicting the immunogenicity of Cytotoxic T Lymphocyte (CTL) epitopes derived from diverse pathogens, given a Human Leukocyte Antigen (HLA) genotype. The model was trained and tested on experimental data on the relative immunodominance of CTL epitopes in Human Immunodeficiency Virus infection. The method is more accurate than publicly available models. Our model predicts that only a fraction of SARS-CoV-2 epitopes that have been predicted to bind to HLA molecules is immunogenic. The immunogenic CTL epitopes across all SARS-CoV-2 proteins are predicted to provide broad population coverage, but the immunogenic epitopes in the SARS-CoV-2 spike protein alone are unlikely to do so. Our model predicts that several immunogenic SARS-CoV-2 CTL epitopes are identical to those contained in low-pathogenicity coronaviruses circulating in the population. Thus, we suggest that some level of CTL immunity against COVID-19 may be present in some individuals prior to SARS-CoV-2 infection.

Broad cross-reactivity of the T-cell repertoire achieves specific and sufficiently rapid target searching

The molecular recognition of T-cell receptors is the hallmark of the adaptive immunity. Given the finiteness of the T-cell repertoire, individual T-cell receptors are necessary to be cross-reactive to multiple antigenic peptides. In this study, we quantify the variability of the cross-reactivity by using a string model that estimates the binding affinity between two sequences of amino acids. We examine sequences of 10,000 human T-cell receptors and 10,000 antigenic peptides, and obtain a full spectrum of cross-reactivity of the receptor-peptide binding. Then, we find that the cross-reactivity spectrum is broad. Some T-cells are reactive to 1000 peptides, but some T-cells are reactive to only one or two peptides. Since the degree of cross-reactivity has a correlation with the (un)binding affinity of receptors, we further investigate how the broad cross-reactivity affects the target searching of T-cells. High cross-reactive T-cells may not require many trials for searching correct targets, but they may spend long time to unbind from incorrect targets. In contrast, low cross-reactive T-cells may not spend long time to ignore incorrect targets, but they require many trials for screening correct targets. We evaluate this hypothesis, and show that the broad cross-reactivity of the natural T-cell repertoire can balance the trade-off between the rapid screening and unbinding penalty.

Influence of correlated antigen presentation on T-cell negative selection in the thymus

The thymus is the primary organ for the generation of naive T cells, a key component of the immune system. Tolerance of T cells to self is achieved primarily in the thymic medulla, where immature T cells (thymocytes) sample self-peptides presented by medullary thymic epithelial cells (mTECs). A sufficiently strong interaction activates the thymocytes leading to negative selection. A key question of current interest is whether there is any structure in the manner in which mTECs present peptides: can any mTEC present any peptide at any time, or are there particular patterns of correlated peptide presentation? We investigate this question using a mathematical model of negative selection. We find that correlated patterns of peptide presentation may be advantageous in negatively selecting low-degeneracy thymocytes (that is, those thymocytes which respond to relatively few peptides). We also quantify the probability that an auto-reactive thymocyte exits the thymus before it encounters a cognate antigen. The results suggest that heterogeneity of gene co-expression in mTECs has an effect on the probability of escape of autoreactive thymocytes.

How nonuniform contact profiles of T cell receptors modulate thymic selection outcomes

T cell receptors (TCRs) bind foreign or self-peptides attached to major histocompatibility complex (MHC) molecules, and the strength of this interaction determines T cell activation. Optimizing the ability of T cells to recognize a diversity of foreign peptides yet be tolerant of self-peptides is crucial for the adaptive immune system to properly function. This is achieved by selection of T cells in the thymus, where immature T cells expressing unique, stochastically generated TCRs interact with a large number of self-peptide-MHC; if a TCR does not bind strongly enough to any self-peptide-MHC, or too strongly with at least one self-peptide-MHC, the T cell dies. Past theoretical work cast thymic selection as an extreme value problem and characterized the statistical enrichment or depletion of amino acids in the postselection TCR repertoire, showing how T cells are selected to be able to specifically recognize peptides derived from diverse pathogens yet have limited self-reactivity. Here, we investigate how the diversity of the postselection TCR repertoire is modified when TCRs make nonuniform contacts with peptide-MHC. Specifically, we were motivated by recent experiments showing that amino acids at certain positions of a TCR sequence have large effects on thymic selection outcomes, and crystal structure data that reveal a nonuniform contact profile between a TCR and its peptide-MHC ligand. Using a representative TCR contact profile as an illustration, we show via simulations that the statistical enrichment or depletion of amino acids now varies by position according to the contact profile, and, importantly, it depends on the implementation of nonuniform contacts during thymic selection. We explain these nontrivial results analytically. Our study has implications for understanding the selection forces that shape the functionality of the postselection TCR repertoire.

Organs on chip approach: a tool to evaluate cancer -immune cells interactions

In this paper we discuss the applicability of numerical descriptors and statistical physics concepts to characterize complex biological systems observed at microscopic level through organ on chip approach. To this end, we employ data collected on a microfluidic platform in which leukocytes can move through suitably built channels toward their target. Leukocyte behavior is recorded by standard time lapse imaging. In particular, we analyze three groups of human peripheral blood mononuclear cells (PBMC): heterozygous mutants (in which only one copy of the FPR1 gene is normal), homozygous mutants (in which both alleles encoding FPR1 are loss-of-function variants) and cells from 'wild type' donors (with normal expression of FPR1). We characterize the migration of these cells providing a quantitative confirmation of the essential role of FPR1 in cancer chemotherapy response. Indeed wild type PBMC perform biased random walks toward chemotherapy-treated cancer cells establishing persistent interactions with them. Conversely, heterozygous mutants present a weaker bias in their motion and homozygous mutants perform rather uncorrelated random walks, both failing to engage with their targets. We next focus on wild type cells and study the interactions of leukocytes with cancerous cells developing a novel heuristic procedure, inspired by Lyapunov stability in dynamical systems.

Effects of thymic selection on T cell recognition of foreign and tumor antigenic peptides

The advent of cancer immunotherapy has generated renewed hope for the treatment of many malignancies by introducing a number of novel strategies that exploit various properties of the immune system. These therapies are based on the idea that cytotoxic T lymphocytes (CTLs) directly recognize and respond to tumor-associated neoantigens (TANs) in much the same way as they would to foreign peptides presented on cell surfaces. To date, however, nearly all attempts to optimize immunotherapeutic strategies have been empirical. Here, we develop a model of T cell selection based on the assumption of random interaction strengths between a self-peptide and the various T cell receptors. The model enables the analytical study of the effects of selection on the CTL recognition of TANs and completely foreign peptides and can estimate the number of CTLs that can detect donor-matched transplants. We show that negative selection thresholds chosen to reflect experimentally observed thymic survival rates result in near-optimal production of T cells that are capable of surviving selection and recognizing foreign antigen. These analytical results are confirmed by simulation.

A Perspective on the Role of Computational Models in Immunology

This is an exciting time for immunology because the future promises to be replete with exciting new discoveries that can be translated to improve health and treat disease in novel ways. Immunologists are attempting to answer increasingly complex questions concerning phenomena that range from the genetic, molecular, and cellular scales to that of organs, whole animals or humans, and populations of humans and pathogens. An important goal is to understand how the many different components involved interact with each other within and across these scales for immune responses to emerge, and how aberrant regulation of these processes causes disease. To aid this quest, large amounts of data can be collected using high-throughput instrumentation. The nonlinear, cooperative, and stochastic character of the interactions between components of the immune system as well as the overwhelming amounts of data can make it difficult to intuit patterns in the data or a mechanistic understanding of the phenomena being studied. Computational models are increasingly important in confronting and overcoming these challenges. I first describe an iterative paradigm of research that integrates laboratory experiments, clinical data, computational inference, and mechanistic computational models. I then illustrate this paradigm with a few examples from the recent literature that make vivid the power of bringing together diverse types of computational models with experimental and clinical studies to fruitfully interrogate the immune system.

Can the Immune System Perform a t-Test?

The self-nonself discrimination hypothesis remains a landmark concept in immunology. It proposes that tolerance breaks down in the presence of nonself antigens. In strike contrast, in statistics, occurrence of nonself elements in a sample (i.e., outliers) is not obligatory to violate the null hypothesis. Very often, what is crucial is the combination of (self) elements in a sample. The two views on how to detect a change seem challengingly different and it could seem difficult to conceive how immunological cellular interactions could trigger responses with a precision comparable to some statistical tests. Here it is shown that frustrated cellular interactions reconcile the two views within a plausible immunological setting. It is proposed that the adaptive immune system can be promptly activated either when nonself ligands are detected or self-ligands occur in abnormal combinations. In particular we show that cellular populations behaving in this way could perform location statistical tests, with performances comparable to t or KS tests, or even more general data mining tests such as support vector machines or random forests. In more general terms, this work claims that plausible immunological models should provide accurate detection mechanisms for host protection and, furthermore, that investigation on mechanisms leading to improved detection in “in silico” models can help unveil how the real immune system works.

Coevolutionary immune system dynamics driving pathogen speciation

We introduce and analyze a within-host dynamical model of the coevolution between rapidly mutating pathogens and the adaptive immune response. Pathogen mutation and a homeostatic constraint on lymphocytes both play a role in allowing the development of chronic infection, rather than quick pathogen clearance. The dynamics of these chronic infections display emergent structure, including branching patterns corresponding to asexual pathogen speciation, which is fundamentally driven by the coevolutionary interaction. Over time, continued branching creates an increasingly fragile immune system, and leads to the eventual catastrophic loss of immune control.

Anergy in self-directed B lymphocytes: A statistical mechanics perspective

Self-directed lymphocytes may evade clonal deletion at ontogenesis but still remain harmless due to a mechanism called clonal anergy. For B-lymphocytes, two major explanations for anergy developed over the last decades: according to Varela theory, anergy stems from a proper orchestration of the whole B-repertoire, such that self-reactive clones, due to intensive feed-back from other clones, display strong inertia when mounting a response. Conversely, according to the model of cognate response, self-reacting cells are not stimulated by helper lymphocytes and the absence of such signaling yields anergy. Through statistical mechanics we show that helpers do not prompt activation of a sub-group of B-cells: remarkably, the latter are just those broadly interacting in the idiotypic network. Hence Varela theory can finally be reabsorbed into the prevailing framework of the cognate response model. Further, we show how the B-repertoire architecture may emerge, where highly connected clones are self-directed as a natural consequence of ontogenetic learning.

Theories and quantification of thymic selection

The peripheral T cell repertoire is sculpted from prototypic T cells in the thymus bearing randomly generated T cell receptors (TCR) and by a series of developmental and selection steps that remove cells that are unresponsive or overly reactive to self-peptide–MHC complexes. The challenge of understanding how the kinetics of T cell development and the statistics of the selection processes combine to provide a diverse but self-tolerant T cell repertoire has invited quantitative modeling approaches, which are reviewed here.

The past, present, and future of immune repertoire biology - the rise of next-generation repertoire analysis

T and B cell repertoires are collections of lymphocytes, each characterized by its antigen-specific receptor. We review here classical technologies and analysis strategies developed to assess immunoglobulin (IG) and T cell receptor (TR) repertoire diversity, and describe recent advances in the field. First, we describe the broad range of available methodological tools developed in the past decades, each of which answering different questions and showing complementarity for progressive identification of the level of repertoire alterations: global overview of the diversity by flow cytometry, IG repertoire descriptions at the protein level for the identification of IG reactivities, IG/TR CDR3 spectratyping strategies, and related molecular quantification or dynamics of T/B cell differentiation. Additionally, we introduce the recent technological advances in molecular biology tools allowing deeper analysis of IG/TR diversity by next-generation sequencing (NGS), offering systematic and comprehensive sequencing of IG/TR transcripts in a short amount of time. NGS provides several angles of analysis such as clonotype frequency, CDR3 diversity, CDR3 sequence analysis, V allele identification with a quantitative dimension, therefore requiring high-throughput analysis tools development. In this line, we discuss the recent efforts made for nomenclature standardization and ontology development. We then present the variety of available statistical analysis and modeling approaches developed with regards to the various levels of diversity analysis, and reveal the increasing sophistication of those modeling approaches. To conclude, we provide some examples of recent mathematical modeling strategies and perspectives that illustrate the active rise of a "next-generation" of repertoire analysis.

Quorum sensing allows T cells to discriminate between self and nonself

T cells orchestrate pathogen-specific adaptive immune responses by identifying peptides derived from pathogenic proteins that are displayed on the surface of infected cells. Host cells also display peptide fragments from the host's own proteins. Incorrectly identifying peptides derived from the body's own proteome as pathogenic can result in autoimmune disease. To minimize autoreactivity, immature T cells that respond to self-peptides are deleted in the thymus by a process called negative selection. However, negative selection is imperfect, and autoreactive T cells exist in healthy individuals. To understand how autoimmunity is yet avoided, without loss of responsiveness to pathogens, we have developed a model of T-cell training and response. Our model shows that T cells reliably respond to infection and avoid autoimmunity because collective decisions made by the T-cell population, rather than the responses of individual T cells, determine biological outcomes. The theory is qualitatively consistent with experimental data and yields a criterion for thymic selection to be adequate for suppressing autoimmunity.

Parallel processing in immune networks

In this work, we adopt a statistical-mechanics approach to investigate basic, systemic features exhibited by adaptive immune systems. The lymphocyte network made by B cells and T cells is modeled by a bipartite spin glass, where, following biological prescriptions, links connecting B cells and T cells are sparse. Interestingly, the dilution performed on links is shown to make the system able to orchestrate parallel strategies to fight several pathogens at the same time; this multitasking capability constitutes a remarkable, key property of immune systems as multiple antigens are always present within the host. We also define the stochastic process ruling the temporal evolution of lymphocyte activity and show its relaxation toward an equilibrium measure allowing statistical-mechanics investigations. Analytical results are compared with Monte Carlo simulations and signal-to-noise outcomes showing overall excellent agreement. Finally, within our model, a rationale for the experimentally well-evidenced correlation between lymphocytosis and autoimmunity is achieved; this sheds further light on the systemic features exhibited by immune networks.

Positive feedback produces broad distributions in maximum activation attained within a narrow time window in stochastic biochemical reactions

How do single cell fate decisions induced by activation of key signaling proteins above threshold concentrations within a time interval are affected by stochastic fluctuations in biochemical reactions? We address this question using minimal models of stochastic chemical reactions commonly found in cell signaling and gene regulatory systems. Employing exact solutions and semi-analytical methods we calculate distributions of the maximum value (N) of activated species concentrations (P(max)(N)) and the time (t) taken to reach the maximum value (P(max)(t)) within a time interval in the minimal models. We find, the presence of positive feedback interactions make P(max)(N) more spread out with a higher "peakedness" in P(max)(t). Thus positive feedback interactions may help single cells to respond sensitively to a stimulus when cell decision processes require upregulation of activated forms of key proteins to a threshold number within a time window.

Organization and evolution of synthetic idiotypic networks

We introduce a class of weighted graphs whose properties are meant to mimic the topological features of idiotypic networks, namely, the interaction networks involving the B core of the immune system. Each node is endowed with a bit string representing the idiotypic specificity of the corresponding B cell, and the proper distance between any couple of bit strings provides the coupling strength between the two nodes. We show that a biased distribution of the entries in bit strings can yield fringes in the (weighted) degree distribution, small-world features, and scaling laws, in agreement with experimental findings. We also investigate the role of aging, thought of as a progressive increase in the degree of bias in bit strings, and we show that it can possibly induce mild percolation phenomena, which are investigated too.

The effect of mutations on the alloreactive T cell receptor/peptide-MHC interface structure: a molecular dynamics study

T cells orchestrate adaptive, pathogen-specific immune responses. T cells have a surface receptor (called TCR) whose ligands are complexes (pMHCs) of peptides (derived from pathogens or host proteins) and major histocompatibility complex proteins (MHCs). MHC proteins vary between hosts. During organ transplants, host TCRs interact with peptides present in complex with genetically different MHCs. This usually causes a vigorous immune response: alloreactivity. Studies of alloreactive protein interactions have yielded results that present a puzzle. Some crystallographic studies concluded that the alloreactive TCR/MHC interface is essentially unaffected by changing the TCR peptide-binding region, suggesting that the peptide does not influence the interface. Another biochemical study concluded from mutation data that different peptides can alter the binding interface with the same TCR. To explore the origin of this puzzle, we used molecular dynamics simulations to study the dependence of the TCR/pMHC interface on changes in both the peptide and the TCR. Our simulations show that the footprint of the TCR on the pMHC is insensitive to mutations of the TCR peptide-binding loops, but peptide mutations can make multiple local changes to TCR/pMHC contacts. Therefore, our results demonstrate that the structural and mutation data do not conflict and reveal how subtle, but important, characteristics of the alloreactive TCR/pMHC interface are influenced by the TCR and the peptide.

Predicting interactions between T cell receptors and MHC-peptide complexes

Conserved interactions between T cell receptors (TCRs) and major histocompatibility complex (MHC) proteins with bound peptide antigens are not well understood. In order to gain a better understanding of the interaction modes of human TCR variable (V) regions, we have performed a structural analysis of the TCRs bound to their MHC-peptide ligands in human, using the available structural models determined by X-ray crystallography. We identified important differences to previous studies in which such interactions were evaluated. Based on the interactions found in the actual experimental structures we developed the first rule-based approach for predicting the ability of TCR residues in the complementarity-determining region (CDR) 1, CDR2, and CDR3 loops to interact with the MHC-peptide antigen complex. Two relatively simple algorithms show good performance under cross validation.

Statistical mechanical concepts in immunology

Higher organisms, such as humans, have an adaptive immune system that usually enables them to successfully combat diverse (and evolving) microbial pathogens. The adaptive immune system is not preprogrammed to respond to prescribed pathogens. Yet it mounts pathogen-specific responses against diverse microbes and establishes memory of past infections (the basis of vaccination). Although major advances have been made in understanding pertinent molecular and cellular phenomena, the mechanistic principles that govern many aspects of an immune response are not known. We illustrate how complementary approaches from the physical and life sciences can help confront this challenge. Specifically, we describe work that brings together statistical mechanics and cell biology to shed light on how key molecular/cellular components of the adaptive immune system are selected to enable pathogen-specific responses. We hope these examples encourage physical chemists to work at this crossroad of disciplines where fundamental discoveries with implications for human health might be made.

Effects of thymic selection of the T-cell repertoire on HLA class I-associated control of HIV infection

Without therapy, most people infected with human immunodeficiency virus (HIV) ultimately progress to AIDS. Rare individuals ('elite controllers') maintain very low levels of HIV RNA without therapy, thereby making disease progression and transmission unlikely. Certain HLA class I alleles are markedly enriched in elite controllers, with the highest association observed for HLA-B57 (ref. 1). Because HLA molecules present viral peptides that activate CD8(+) T cells, an immune-mediated mechanism is probably responsible for superior control of HIV. Here we describe how the peptide-binding characteristics of HLA-B57 molecules affect thymic development such that, compared to other HLA-restricted T cells, a larger fraction of the naive repertoire of B57-restricted clones recognizes a viral epitope, and these T cells are more cross-reactive to mutants of targeted epitopes. Our calculations predict that such a T-cell repertoire imposes strong immune pressure on immunodominant HIV epitopes and emergent mutants, thereby promoting efficient control of the virus. Supporting these predictions, in a large cohort of HLA-typed individuals, our experiments show that the relative ability of HLA-B alleles to control HIV correlates with their peptide-binding characteristics that affect thymic development. Our results provide a conceptual framework that unifies diverse empirical observations, and have implications for vaccination strategies.