Computational identification of antibody-binding epitopes from mimotope datasets.

Introduction: A fundamental challenge in computational vaccinology is that most B-cell epitopes are conformational and therefore hard to predict from sequence alone. Another significant challenge is that a great deal of the amino acid sequence of a viral surface protein might not in fact be antigenic. Thus, identifying the regions of a protein that are most promising for vaccine design based on the degree of surface exposure may not lead to a clinically relevant immune response. Methods: Linear peptides selected by phage display experiments that have high affinity to the monoclonal antibody of interest ("mimotopes") usually have similar physicochemical properties to the antigen epitope corresponding to that antibody. The sequences of these linear peptides can be used to find possible epitopes on the surface of the antigen structure or a homology model of the antigen in the absence of an antigen-antibody complex structure. Results and Discussion: Herein we describe two novel methods for mapping mimotopes to epitopes. The first is a novel algorithm named MimoTree that allows for gaps in the mimotopes and epitopes on the antigen. More specifically, a mimotope may have a gap that does not match to the epitope to allow it to adopt a conformation relevant for binding to an antibody, and residues may similarly be discontinuous in conformational epitopes. MimoTree is a fully automated epitope detection algorithm suitable for the identification of conformational as well as linear epitopes. The second is an ensemble approach, which combines the prediction results from MimoTree and two existing methods.

Looking under the lamp-post: quantifying the performance of contact tracing in the United States during the SARS-CoV-2 pandemic.

Contact tracing forms a crucial part of the public-health toolbox in mitigating and understanding emergent pathogens and nascent disease outbreaks. Contact tracing in the United States was conducted during the pre-Omicron phase of the ongoing COVID-19 pandemic. This tracing relied on voluntary reporting and responses, often using rapid antigen tests due to lack of accessibility to PCR tests. These limitations, combined with SARS-CoV-2's propensity for asymptomatic transmission, raise the question "how reliable was contact tracing for COVID-19 in the United States"? We answered this question using a Markov model to examine the efficiency with which transmission could be detected based on the design and response rates of contact tracing studies in the United States. Our results suggest that contact tracing protocols in the U.S. are unlikely to have identified more than 1.65% (95% uncertainty interval: 1.62-1.68%) of transmission events with PCR testing and 1.00% (95% uncertainty interval 0.98-1.02%) with rapid antigen testing. When considering a more robust contact tracing scenario, based on compliance rates in East Asia with PCR testing, this increases to 62.7% (95% uncertainty interval: 62.6-62.8%). We did not assume presence of asymptomatic transmission or superspreading, making our estimates upper bounds on the actual percentages traced. These findings highlight the limitations in interpretability for studies of SARS-CoV-2 disease spread based on U.S. contact tracing and underscore the vulnerability of the population to future disease outbreaks, for SARS-CoV-2 and other pathogens.

Antibody escape, the risk of serotype formation, and rapid immune waning: Modeling the implications of SARS-CoV-2 immune evasion.

As the COVID-19 pandemic progresses, widespread community transmission of SARS-CoV-2 has ushered in a volatile era of viral immune evasion rather than the much-heralded stability of "endemicity" or "herd immunity." At this point, an array of viral strains has rendered essentially all monoclonal antibody therapeutics obsolete and strongly undermined the impact of vaccinal immunity on SARS-CoV-2 transmission. In this work, we demonstrate that antibody escape resulting in evasion of pre-existing immunity is highly evolutionarily favored and likely to cause waves of short-term transmission. In the long-term, invading strains that induce weak cross-immunity against pre-existing strains may co-circulate with those pre-existing strains. This would result in the formation of serotypes that increase disease burden, complicate SARS-CoV-2 control, and raise the potential for increases in viral virulence. Less durable immunity does not drive positive selection as a trait, but such strains may transmit at high levels if they establish. Overall, our results draw attention to the importance of inter-strain cross-immunity as a driver of transmission trends and the importance of early immune evasion data to predict the trajectory of the pandemic.

Vaccines Alone Cannot Slow the Evolution of SARS-CoV-2.

The rapid emergence of immune-evading viral variants of SARS-CoV-2 calls into question the practicality of a vaccine-only public-health strategy for managing the ongoing COVID-19 pandemic. It has been suggested that widespread vaccination is necessary to prevent the emergence of future immune-evading mutants. Here, we examined that proposition using stochastic computational models of viral transmission and mutation. Specifically, we looked at the likelihood of emergence of immune escape variants requiring multiple mutations and the impact of vaccination on this process. Our results suggest that the transmission rate of intermediate SARS-CoV-2 mutants will impact the rate at which novel immune-evading variants appear. While vaccination can lower the rate at which new variants appear, other interventions that reduce transmission can also have the same effect. Crucially, relying solely on widespread and repeated vaccination (vaccinating the entire population multiple times a year) is not sufficient to prevent the emergence of novel immune-evading strains, if transmission rates remain high within the population. Thus, vaccines alone are incapable of slowing the pace of evolution of immune evasion, and vaccinal protection against severe and fatal outcomes for COVID-19 patients is therefore not assured.

No magic bullet: Limiting in-school transmission in the face of variable SARS-CoV-2 viral loads.

In the face of a long-running pandemic, understanding the drivers of ongoing SARS-CoV-2 transmission is crucial for the rational management of COVID-19 disease burden. Keeping schools open has emerged as a vital societal imperative during the pandemic, but in-school transmission of SARS-CoV-2 can contribute to further prolonging the pandemic. In this context, the role of schools in driving SARS-CoV-2 transmission acquires critical importance. Here we model in-school transmission from first principles to investigate the effectiveness of layered mitigation strategies on limiting in-school spread. We examined the effect of masks and air quality (ventilation, filtration and ionizers) on steady-state viral load in classrooms, as well as on the number of particles inhaled by an uninfected person. The effectiveness of these measures in limiting viral transmission was assessed for variants with different levels of mean viral load (ancestral, Delta, Omicron). Our results suggest that a layered mitigation strategy can be used effectively to limit in-school transmission, with certain limitations. First, poorly designed strategies (insufficient ventilation, no masks, staying open under high levels of community transmission) will permit in-school spread even if some level of mitigation is present. Second, for viral variants that are sufficiently contagious, it may be difficult to construct any set of interventions capable of blocking transmission once an infected individual is present, underscoring the importance of other measures. Our findings provide practical recommendations; in particular, the use of a layered mitigation strategy that is designed to limit transmission, with other measures such as frequent surveillance testing and smaller class sizes (such as by offering remote schooling options to those who prefer it) as needed.

Genomic analysis of early-stage lung cancer reveals a role for TP53 mutations in distant metastasis.

Patients with non-small cell lung cancer (NSCLC) who have distant metastases have a poor prognosis. To determine which genomic factors of the primary tumor are associated with metastasis, we analyzed data from 759 patients originally diagnosed with stage I-III NSCLC as part of the AACR Project GENIE Biopharma Collaborative consortium. We found that TP53 mutations were significantly associated with the development of new distant metastases. TP53 mutations were also more prevalent in patients with a history of smoking, suggesting that these patients may be at increased risk for distant metastasis. Our results suggest that additional investigation of the optimal management of patients with early-stage NSCLC harboring TP53 mutations at diagnosis is warranted in light of their higher likelihood of developing new distant metastases.

Transcriptional differences between JAK2-V617F and wild-type bone marrow cells in patients with myeloproliferative neoplasms.

The JAK2-V617F mutation is the most common cause of myeloproliferative neoplasms. Although experiments have revealed that this gain-of-function mutation is associated with myeloid blood cell expansion and increased production of white cells, red cells, and platelets, the transcriptional consequences of the JAK2-V617F mutation in different cellular compartments of the bone marrow have not yet been fully elucidated. To study the direct effects of JAK2-V617F on bone marrow cells in patients with myeloproliferative neoplasms, we performed joint single-cell RNA sequencing and JAK2 genotyping on CD34+-enriched cells from eight patients with newly diagnosed essential thrombocythemia or polycythemia vera. We found that the JAK2-V617F mutation increases the expression of interferon-response genes (e.g., HLAs) and the leptin receptor in hematopoietic progenitor cells. Furthermore, we sequenced a population of CD34- bone marrow monocytes and found that the JAK2 mutation increased expression of intermediate monocyte genes and the fibrocyte-associated surface protein SLAMF7 in these cells.

Rapid relaxation of pandemic restrictions after vaccine rollout favors growth of SARS-CoV-2 variants: A model-based analysis.

The development and deployment of several SARS-CoV-2 vaccines in a little over a year is an unprecedented achievement of modern medicine. The high levels of efficacy against transmission for some of these vaccines makes it feasible to use them to suppress SARS-CoV-2 altogether in regions with high vaccine acceptance. However, viral variants with reduced susceptibility to vaccinal and natural immunity threaten the utility of vaccines, particularly in scenarios where a return to pre-pandemic conditions occurs before the suppression of SARS-CoV-2 transmission. In this work we model the situation in the United States in May-June 2021, to demonstrate how pre-existing variants of SARS-CoV-2 may cause a rebound wave of COVID-19 in a matter of months under a certain set of conditions. A high burden of morbidity (and likely mortality) remains possible, even if the vaccines are partially effective against new variants and widely accepted. Our modeling suggests that variants that are already present within the population may be capable of quickly defeating the vaccines as a public health intervention, a serious potential limitation for strategies that emphasize rapid reopening before achieving control of SARS-CoV-2.

Controlling long-term SARS-CoV-2 infections can slow viral evolution and reduce the risk of treatment failure.

The rapid emergence and expansion of novel SARS-CoV-2 variants threatens our ability to achieve herd immunity for COVID-19. These novel SARS-CoV-2 variants often harbor multiple point mutations, conferring one or more evolutionarily advantageous traits, such as increased transmissibility, immune evasion and longer infection duration. In a number of cases, variant emergence has been linked to long-term infections in individuals who were either immunocompromised or treated with convalescent plasma. In this paper, we used a stochastic evolutionary modeling framework to explore the emergence of fitter variants of SARS-CoV-2 during long-term infections. We found that increased viral load and infection duration favor emergence of such variants. While the overall probability of emergence and subsequent transmission from any given infection is low, on a population level these events occur fairly frequently. Targeting these low-probability stochastic events that lead to the establishment of novel advantageous viral variants might allow us to slow the rate at which they emerge in the patient population, and prevent them from spreading deterministically due to natural selection. Our work thus suggests practical ways to achieve control of long-term SARS-CoV-2 infections, which will be critical for slowing the rate of viral evolution.

Individually optimal choices can be collectively disastrous in COVID-19 disease control.

BACKGROUND: The word 'pandemic' conjures dystopian images of bodies stacked in the streets and societies on the brink of collapse. Despite this frightening picture, denialism and noncompliance with public health measures are common in the historical record, for example during the 1918 Influenza pandemic or the 2015 Ebola epidemic. The unique characteristics of SARS-CoV-2-its high basic reproduction number (R0), time-limited natural immunity and considerable potential for asymptomatic spread-exacerbate the public health repercussions of noncompliance with interventions (such as vaccines and masks) to limit disease transmission. Our work explores the rationality and impact of noncompliance with measures aimed at limiting the spread of SARS-CoV-2. METHODS: In this work, we used game theory to explore when noncompliance confers a perceived benefit to individuals. We then used epidemiological modeling to predict the impact of noncompliance on control of SARS-CoV-2, demonstrating that the presence of a noncompliant subpopulation prevents suppression of disease spread. RESULTS: Our modeling demonstrates that noncompliance is a Nash equilibrium under a broad set of conditions and that the existence of a noncompliant population can result in extensive endemic disease in the long-term after a return to pre-pandemic social and economic activity. Endemic disease poses a threat for both compliant and noncompliant individuals; all community members are protected if complete suppression is achieved, which is only possible with a high degree of compliance. For interventions that are highly effective at preventing disease spread, however, the consequences of noncompliance are borne disproportionately by noncompliant individuals. CONCLUSIONS: In sum, our work demonstrates the limits of free-market approaches to compliance with disease control measures during a pandemic. The act of noncompliance with disease intervention measures creates a negative externality, rendering suppression of SARS-CoV-2 spread ineffective. Our work underscores the importance of developing effective strategies for prophylaxis through public health measures aimed at complete suppression and the need to focus on compliance at a population level.

Risk of rapid evolutionary escape from biomedical interventions targeting SARS-CoV-2 spike protein.

The spike protein receptor-binding domain (RBD) of SARS-CoV-2 is the molecular target for many vaccines and antibody-based prophylactics aimed at bringing COVID-19 under control. Such a narrow molecular focus raises the specter of viral immune evasion as a potential failure mode for these biomedical interventions. With the emergence of new strains of SARS-CoV-2 with altered transmissibility and immune evasion potential, a critical question is this: how easily can the virus escape neutralizing antibodies (nAbs) targeting the spike RBD? To answer this question, we combined an analysis of the RBD structure-function with an evolutionary modeling framework. Our structure-function analysis revealed that epitopes for RBD-targeting nAbs overlap one another substantially and can be evaded by escape mutants with ACE2 affinities comparable to the wild type, that are observed in sequence surveillance data and infect cells in vitro. This suggests that the fitness cost of nAb-evading mutations is low. We then used evolutionary modeling to predict the frequency of immune escape before and after the widespread presence of nAbs due to vaccines, passive immunization or natural immunity. Our modeling suggests that SARS-CoV-2 mutants with one or two mildly deleterious mutations are expected to exist in high numbers due to neutral genetic variation, and consequently resistance to vaccines or other prophylactics that rely on one or two antibodies for protection can develop quickly -and repeatedly- under positive selection. Predicted resistance timelines are comparable to those of the decay kinetics of nAbs raised against vaccinal or natural antigens, raising a second potential mechanism for loss of immunity in the population. Strategies for viral elimination should therefore be diversified across molecular targets and therapeutic modalities.

Reconstructing the Lineage Histories and Differentiation Trajectories of Individual Cancer Cells in Myeloproliferative Neoplasms.

Some cancers originate from a single mutation event in a single cell. Blood cancers known as myeloproliferative neoplasms (MPNs) are thought to originate when a driver mutation is acquired by a hematopoietic stem cell (HSC). However, when the mutation first occurs in individuals and how it affects the behavior of HSCs in their native context is not known. Here we quantified the effect of the JAK2-V617F mutation on the self-renewal and differentiation dynamics of HSCs in treatment-naive individuals with MPNs and reconstructed lineage histories of individual HSCs using somatic mutation patterns. We found that JAK2-V617F mutations occurred in a single HSC several decades before MPN diagnosis-at age 9 ± 2 years in a 34-year-old individual and at age 19 ± 3 years in a 63-year-old individual-and found that mutant HSCs have a selective advantage in both individuals. These results highlight the potential of harnessing somatic mutations to reconstruct cancer lineages.

Fitness variation in isogenic populations leads to a novel evolutionary mechanism for crossing fitness valleys.

Individuals in a population often have different fitnesses even when they have identical genotypes, but the effect of this variation on the evolution of a population through complicated fitness landscapes is unknown. Here, we investigate how populations with non-genetic fitness variation cross fitness valleys, common barriers to adaptation in rugged fitness landscapes in which a population must pass through a deleterious intermediate to arrive at a final advantageous stage. We develop a stochastic computational model describing the dynamics of an asexually reproducing population crossing a fitness valley, in which individuals of the same evolutionary stage can have variable fitnesses. We find that fitness variation that persists over multiple generations increases the rate of valley crossing through a novel evolutionary mechanism different from previously characterized mechanisms such as stochastic tunneling. By reducing the strength of selection against deleterious intermediates, persistent fitness variation allows for faster adaptation through rugged fitness landscapes.

Microfluidic squeezing for intracellular antigen loading in polyclonal B-cells as cellular vaccines.

B-cells are promising candidate autologous antigen-presenting cells (APCs) to prime antigen-specific T-cells both in vitro and in vivo. However to date, a significant barrier to utilizing B-cells as APCs is their low capacity for non-specific antigen uptake compared to "professional" APCs such as dendritic cells. Here we utilize a microfluidic device that employs many parallel channels to pass single cells through narrow constrictions in high throughput. This microscale "cell squeezing" process creates transient pores in the plasma membrane, enabling intracellular delivery of whole proteins from the surrounding medium into B-cells via mechano-poration. We demonstrate that both resting and activated B-cells process and present antigens delivered via mechano-poration exclusively to antigen-specific CD8(+)T-cells, and not CD4(+)T-cells. Squeezed B-cells primed and expanded large numbers of effector CD8(+)T-cells in vitro that produced effector cytokines critical to cytolytic function, including granzyme B and interferon-γ. Finally, antigen-loaded B-cells were also able to prime antigen-specific CD8(+)T-cells in vivo when adoptively transferred into mice. Altogether, these data demonstrate crucial proof-of-concept for mechano-poration as an enabling technology for B-cell antigen loading, priming of antigen-specific CD8(+)T-cells, and decoupling of antigen uptake from B-cell activation.

Structure-based programming of lymph-node targeting in molecular vaccines.

In cancer patients, visual identification of sentinel lymph nodes (LNs) is achieved by the injection of dyes that bind avidly to endogenous albumin, targeting these compounds to LNs, where they are efficiently filtered by resident phagocytes. Here we translate this 'albumin hitchhiking' approach to molecular vaccines, through the synthesis of amphiphiles (amph-vaccines) comprising an antigen or adjuvant cargo linked to a lipophilic albumin-binding tail by a solubility-promoting polar polymer chain. Administration of structurally optimized CpG-DNA/peptide amph-vaccines in mice resulted in marked increases in LN accumulation and decreased systemic dissemination relative to their parent compounds, leading to 30-fold increases in T-cell priming and enhanced anti-tumour efficacy while greatly reducing systemic toxicity. Amph-vaccines provide a simple, broadly applicable strategy to simultaneously increase the potency and safety of subunit vaccines.