Viral surface geometry shapes influenza and coronavirus spike evolution through antibody pressure.

The evolution of circulating viruses is shaped by their need to evade antibody response, which mainly targets the viral spike. Because of the high density of spikes on the viral surface, not all antigenic sites are targeted equally by antibodies. We offer here a geometry-based approach to predict and rank the probability of surface residues of SARS spike (S protein) and influenza H1N1 spike (hemagglutinin) to acquire antibody-escaping mutations utilizing in-silico models of viral structure. We used coarse-grained MD simulations to estimate the on-rate (targeting) of an antibody model to surface residues of the spike protein. Analyzing publicly available sequences, we found that spike surface sequence diversity of the pre-pandemic seasonal influenza H1N1 and the sarbecovirus subgenus highly correlates with our model prediction of antibody targeting. In particular, we identified an antibody-targeting gradient, which matches a mutability gradient along the main axis of the spike. This identifies the role of viral surface geometry in shaping the evolution of circulating viruses. For the 2009 H1N1 and SARS-CoV-2 pandemics, a mutability gradient along the main axis of the spike was not observed. Our model further allowed us to identify key residues of the SARS-CoV-2 spike at which antibody escape mutations have now occurred. Therefore, it can inform of the likely functional role of observed mutations and predict at which residues antibody-escaping mutation might arise.

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.

Damage-induced chromatome dynamics link Ubiquitin ligase and proteasome recruitment to histone loss and efficient DNA repair.

Eukaryotic cells package their genomes around histone octamers. In response to DNA damage, checkpoint activation in yeast induces core histone degradation resulting in 20%-40% reduction in nucleosome occupancy. To gain insight into this process, we developed a new approach to analyze the chromatin-associated proteome comprehensively before and after damage. This revealed extensive changes in protein composition after Zeocin-induced damage. First, core histones and the H1 homolog Hho1 were partially lost from chromatin along with replication, transcription, and chromatin remodeling machineries, while ubiquitin ligases and the proteasome were recruited. We found that the checkpoint- and INO80C-dependent recruitment of five ubiquitin-conjugating factors (Rad6, Bre1, Pep5, Ufd4, and Rsp5) contributes to core and linker histone depletion, reducing chromatin compaction and enhancing DNA locus mobility. Importantly, loss of Rad6/Bre1, Ufd4/TRIP12, and Pep5/VPS11 compromise DNA strand invasion kinetics during homology-driven repair. Thus we provide a comprehensive overview of a functionally relevant genome-wide chromatin response to DNA damage.

Defining and Manipulating B Cell Immunodominance Hierarchies to Elicit Broadly Neutralizing Antibody Responses against Influenza Virus.

The antibody repertoire possesses near-limitless diversity, enabling the adaptive immune system to accommodate essentially any antigen. However, this diversity explores the antigenic space unequally, allowing some pathogens like influenza virus to impose complex immunodominance hierarchies that distract antibody responses away from key sites of virus vulnerability. We developed a computational model of affinity maturation to map the patterns of immunodominance that evolve upon immunization with natural and engineered displays of hemagglutinin (HA), the influenza vaccine antigen. Based on this knowledge, we designed immunization protocols that subvert immune distraction and focus serum antibody responses upon a functionally conserved, but immunologically recessive, target of human broadly neutralizing antibodies. We tested in silico predictions by vaccinating transgenic mice in which antibody diversity was humanized to mirror clinically relevant humoral output. Collectively, our results demonstrate that complex patterns in antibody immunogenicity can be rationally defined and then manipulated to elicit engineered immunity.

Viral surface geometry shapes influenza and coronavirus spike evolution through antibody pressure.

The evolution of circulating viruses is shaped by their need to evade antibody response, which mainly targets the glycoprotein (spike). However, not all antigenic sites are targeted equally by antibodies, leading to complex immunodominance patterns. We used 3D computational models to estimate antibody pressure on the seasonal influenza H1N1 and SARS spikes. Analyzing publically available sequences, we show that antibody pressure, through the geometrical organization of spikes on the viral surface, shaped their mutability. Studying the mutability patterns of SARS-CoV-2 and the 2009 H1N1 pandemic spikes, we find that they are not predominantly shaped by antibody pressure. However, for SARS-CoV-2, we find that over time, it acquired mutations at antibody-accessible positions, which could indicate possible escape as define by our model. We offer a geometry-based approach to predict and rank the probability of surface resides of SARS-CoV-2 spike to acquire antibody escaping mutations.

DNA Damage-Induced Nucleosome Depletion Enhances Homology Search Independently of Local Break Movement.

To determine whether double-strand break (DSB) mobility enhances the physical search for an ectopic template during homology-directed repair (HDR), we tested the effects of factors that control chromatin dynamics, including cohesin loading and kinetochore anchoring. The former but not the latter is altered in response to DSBs. Loss of the nonhistone high-mobility group protein Nhp6 reduces histone occupancy and increases chromatin movement, decompaction, and ectopic HDR. The loss of nucleosome remodeler INO80-C did the opposite. To see whether enhanced HDR depends on DSB mobility or the global chromatin response, we tested the ubiquitin ligase mutant uls1Δ, which selectively impairs local but not global movement in response to a DSB. Strand invasion occurs in uls1Δ cells with wild-type kinetics, arguing that global histone depletion rather than DSB movement is rate limiting for HDR. Impaired break movement in uls1Δ correlates with elevated MRX and cohesin loading, despite normal resection and checkpoint activation.

Advances in Chromatin and Chromosome Research: Perspectives from Multiple Fields.

Nucleosomes package genomic DNA into chromatin. By regulating DNA access for transcription, replication, DNA repair, and epigenetic modification, chromatin forms the nexus of most nuclear processes. In addition, dynamic organization of chromatin underlies both regulation of gene expression and evolution of chromosomes into individualized sister objects, which can segregate cleanly to different daughter cells at anaphase. This collaborative review shines a spotlight on technologies that will be crucial to interrogate key questions in chromatin and chromosome biology including state-of-the-art microscopy techniques, tools to physically manipulate chromatin, single-cell methods to measure chromatin accessibility, computational imaging with neural networks and analytical tools to interpret chromatin structure and dynamics. In addition, this review provides perspectives on how these tools can be applied to specific research fields such as genome stability and developmental biology and to test concepts such as phase separation of chromatin.

Guided nuclear exploration increases CTCF target search efficiency.

The enormous size of mammalian genomes means that for a DNA-binding protein the number of nonspecific, off-target sites vastly exceeds the number of specific, cognate sites. How mammalian DNA-binding proteins overcome this challenge to efficiently locate their target sites is not known. Here, through live-cell single-molecule tracking, we show that CCCTC-binding factor, CTCF, is repeatedly trapped in small zones that likely correspond to CTCF clusters, in a manner that is largely dependent on an internal RNA-binding region (RBRi). We develop a new theoretical model called anisotropic diffusion through transient trapping in zones to explain CTCF dynamics. Functionally, transient RBRi-mediated trapping increases the efficiency of CTCF target search by ~2.5-fold. Overall, our results suggest a 'guided' mechanism where CTCF clusters concentrate diffusing CTCF proteins near cognate binding sites, thus increasing the local ON-rate. We suggest that local guiding may allow DNA-binding proteins to more efficiently locate their target sites.

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.

A Population Dynamics Model for Clonal Diversity in a Germinal Center.

Germinal centers (GCs) are micro-domains where B cells mature to develop high affinity antibodies. Inside a GC, B cells compete for antigen and T cell help, and the successful ones continue to evolve. New experimental results suggest that, under identical conditions, a wide spectrum of clonal diversity is observed in different GCs, and high affinity B cells are not always the ones selected. We use a birth, death and mutation model to study clonal competition in a GC over time. We find that, like all evolutionary processes, diversity loss is inherently stochastic. We study two selection mechanisms, birth-limited and death limited selection. While death limited selection maintains diversity and allows for slow clonal homogenization as affinity increases, birth limited selection results in more rapid takeover of successful clones. Finally, we qualitatively compare our model to experimental observations of clonal selection in mice.

Visualization of Chromatin Decompaction and Break Site Extrusion as Predicted by Statistical Polymer Modeling of Single-Locus Trajectories.

Chromatin moves with subdiffusive and spatially constrained dynamics within the cell nucleus. Here, we use single-locus tracking by time-lapse fluorescence microscopy to uncover information regarding the forces that influence chromatin movement following the induction of a persistent DNA double-strand break (DSB). Using improved time-lapse imaging regimens, we monitor trajectories of tagged DNA loci at a high temporal resolution, which allows us to extract biophysical parameters through robust statistical analysis. Polymer modeling based on these parameters predicts chromatin domain expansion near a DSB and damage extrusion from the domain. Both phenomena are confirmed by live imaging in budding yeast. Calculation of the anomalous exponent of locus movement allows us to differentiate forces imposed on the nucleus through the actin cytoskeleton from those that arise from INO80 remodeler-dependent changes in nucleosome organization. Our analytical approach can be applied to high-density single-locus trajectories obtained in any cell type.

Histone degradation in response to DNA damage enhances chromatin dynamics and recombination rates.

Nucleosomes are essential for proper chromatin organization and the maintenance of genome integrity. Histones are post-translationally modified and often evicted at sites of DNA breaks, facilitating the recruitment of repair factors. Whether such chromatin changes are localized or genome-wide is debated. Here we show that cellular levels of histones drop 20-40% in response to DNA damage. This histone loss occurs from chromatin, is proteasome-mediated and requires both the DNA damage checkpoint and the INO80 nucleosome remodeler. We confirmed reductions in histone levels by stable isotope labeling of amino acids in cell culture (SILAC)-based mass spectrometry, genome-wide nucleosome mapping and fluorescence microscopy. Chromatin decompaction and increased fiber flexibility accompanied histone degradation, both in response to DNA damage and after artificial reduction of histone levels. As a result, recombination rates and DNA-repair focus turnover were enhanced. Thus, we propose that a generalized reduction in nucleosome occupancy is an integral part of the DNA damage response in yeast that provides mechanisms for enhanced chromatin mobility and homology search.

Structural Fluctuations of the Chromatin Fiber within Topologically Associating Domains.

Experiments based on chromosome conformation capture have shown that mammalian genomes are partitioned into topologically associating domains (TADs), within which the chromatin fiber preferentially interacts. TADs may provide three-dimensional scaffolds allowing genes to contact their appropriate distal regulatory DNA sequences (e.g., enhancers) and thus to be properly regulated. Understanding the cell-to-cell and temporal variability of the chromatin fiber within TADs, and what determines them, is thus of great importance to better understand transcriptional regulation. We recently described an equilibrium polymer model that can accurately predict cell-to-cell variation of chromosome conformation within single TADs, from chromosome conformation capture-based data. Here we further analyze the conformational and energetic properties of our model. We show that the chromatin fiber within TADs can easily fluctuate between several conformational states, which are hierarchically organized and are not separated by important free energy barriers, and that this is facilitated by the fact that the chromatin fiber within TADs is close to the onset of the coil-globule transition. We further show that in this dynamic state the properties of the chromatin fiber, and its contact probabilities in particular, are determined in a nontrivial manner not only by site-specific interactions between strongly interacting loci along the fiber, but also by nonlocal correlations between pairs of contacts. Finally, we use live-cell experiments to measure the dynamics of the chromatin fiber in mouse embryonic stem cells, in combination with dynamical simulations, and predict that conformational changes within one TAD are likely to occur on timescales that are much shorter than the duration of one cell cycle. This suggests that genes and their regulatory elements may come together and disassociate several times during a cell cycle. These results have important implications for transcriptional regulation as they support the concept of highly dynamic interactions driven by a complex interplay between site-specific interactions and the intrinsic biophysical properties of the chromatin fiber.

Analysis of Single Locus Trajectories for Extracting In Vivo Chromatin Tethering Interactions.

Is it possible to extract tethering forces applied on chromatin from the statistics of a single locus trajectories imaged in vivo? Chromatin fragments interact with many partners such as the nuclear membrane, other chromosomes or nuclear bodies, but the resulting forces cannot be directly measured in vivo. However, they impact chromatin dynamics and should be reflected in particular in the motion of a single locus. We present here a method based on polymer models and statistics of single trajectories to extract the force characteristics and in particular when they are generated by the gradient of a quadratic potential well. Using numerical simulations of a Rouse polymer and live cell imaging of the MAT-locus located on the yeast Saccharomyces cerevisiae chromosome III, we recover the amplitude and the distance between the observed and the interacting monomer. To conclude, the confined trajectories we observed in vivo reflect local interaction on chromatin.

First-passage distributions in a collective model of anomalous diffusion with tunable exponent.

We consider a model system in which anomalous diffusion is generated by superposition of underlying linear modes with a broad range of relaxation times. In the language of Gaussian polymers, our model corresponds to Rouse (Fourier) modes whose friction coefficients scale as wave number to the power 2-z. A single (tagged) monomer then executes subdiffusion over a broad range of time scales, and its mean square displacement increases as t(alpha) with alpha = 1/z. To demonstrate nontrivial aspects of the model, we numerically study the absorption of the tagged particle in one dimension near an absorbing boundary or in the interval between two such boundaries. We obtain absorption probability densities as a function of time, as well as the position-dependent distribution for unabsorbed particles, at several values of alpha. Each of these properties has features characterized by exponents that depend on alpha. Characteristic distributions found for different values of alpha have similar qualitative features, but are not simply related quantitatively. Comparison of the motion of translocation coordinate of a polymer moving through a pore in a membrane with the diffusing tagged monomer with identical alpha also reveals quantitative differences.