Search

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 (ATP). Moreover, chromatin motion varies with genomic coordinate, pointing towards 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 non-equilibrium mechanisms for forming genomic compartments, which cannot be distinguished from affinity-based folding using structural data alone. As a first step toward disentangling active and passive mechanisms of folding, we discuss a data-driven approach to discern if and how active processes affect genome organization.

The molecules of the cell are compartmentalized into membrane- and non-membrane-bound organelles. Many non-membrane-bound organelles are phase-separated biomolecular condensates with distinct physicochemical properties that can absorb and concentrate specific proteins and nucleic acids involved in discrete biochemical processes. We reasoned that selective condensate partitioning might also occur with small molecule drugs whose targets occur within condensates, and that the therapeutic index and efficacy of such compounds might therefore relate to their ability to partition into condensates. To test this idea, we focused our study on nuclear condensates reported in cell lines, demonstrated they occur in normal human and malignant breast cancer, and developed assays to test clinically active antineoplastic small molecule drugs relative to these condensates.To study the behavior of drugs within condensates, these were modeled in vitro with purified proteins and visualized by fluorescent confocal microscopy. We found that cisplatin, tamoxifen, JQ1, THZ1, and mitoxantrone are concentrated in specific protein condensates in vitro, and that this occurs through physicochemical properties independent of the drug target. For each drug, the small molecule partitioned into the same condensate in vitro in which its established target resides in vivo. A screen of a chemically diverse fluorescent probes and mutant-protein condensates demonstrated that pi-system interactions between aromatic moieties in the protein and small molecule govern concentration in condensates. These results show that clinically important drugs partition into specific protein condensates in vitro by virtue of defined chemical properties, thereby altering their local concentration.Alkylating agents are a class of commonly used antineoplastic compounds, of which cisplatin is a prominent example. In vitro droplet assays revealed that cisplatin is selectively concentrated in transcriptional condensates, and that this ability is required for efficient platination of target DNA. In cell studies revealed that cisplatin preferentially targets DNA contained within MED1 condensates, and disrupts the genetic regulatory elements that compose phase-separated transcriptional condensates. Live cell imaging demonstrated that transcriptional condensates are dissolved by cisplatin, whereas other condensates remain intact. Thus, we conclude that cisplatin preferentially modifies transcriptional condensate-associated DNA in cells, and that this causes selective condensate disruption. The mechanisms that produce drug resistance can provide clues to drug activity in the clinical setting. Investigating the behavior of tamoxifen within ER transcriptional condensates demonstrated that it disrupts these condensates in vitro and on oncogenes in cells; hormonal therapy resistant ESR1 mutations render these condensates resistant. MED1 overexpression, a poorly understood mechanism of tamoxifen resistance, increased the size of ER-MED1 condensates, thereby rending tamoxifen more dilute and ineffective when concentrated therein. This suggest that altering the size and nature of transcription condensates in breast cancer can mediate drug resistance in the clinical setting.Our results show that antineoplastic drugs partition selectively into condensates, that this can occur through physicochemical properties independent of their molecular targets, and that resistance to drugs may occur through condensate altering mechanisms. These results have implications for development of efficacious cancer therapeutics; effective target engagement will depend on factors such as drug partitioning in condensates. Assays of the type described here may thus help optimize condensate partitioning, target engagement, and the therapeutic index of drugs for cancer treatment. Citation Format: Isaac Klein, Ann Boija, Lena Afeyan, Susana Wilson Hawken, Mengyang Fan, Alessandra Dall'Agnese, Ozgur Oksuz, Jonathan Henninger, Krishna Shrinivas, Benjamin Sabari, Ido Sagi, Victoria Clark, Jesse Platt, Mrityunjoy Kar, Patrick McCall, Alicia Zamudio, John Mantiega, Eliot Coffey, Charles Li, Nancy Hannett, Yang Guo, Tim-Michael Decker, Tong Lee, Tinghu Zhang, Jing-Ke Weng, Dylan Taatjes, Arup Chakraborty, Phillip Sharp, Young Tae Chang, Anthony Hyman, Nathanael Gray, Richard Young. Partitioning of cancer therapeutics in nuclear condensates [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr GS3-10.

T cells exhibit remarkable sensitivity and selectivity in detecting and responding to agonist peptides (p) bound to MHC molecules in a sea of self pMHC molecules. Despite much work, understanding of the underlying mechanisms of distinguishing such ligands remains incomplete. Here, we quantify T cell discriminatory capacity using channel capacity, a direct measure of the signaling network's ability to discriminate between antigen-presenting cells (APCs) displaying either self ligands or a mixture of self and agonist ligands. This metric shows how differences in information content between these two types of peptidomes are decoded by the topology and rates of kinetic proofreading signaling steps inside T cells. Using channel capacity, we constructed numerically substantiated hypotheses to explain the discriminatory role of a recently identified slow LAT Y132 phosphorylation step. Our results revealed that in addition to the number and kinetics of sequential signaling steps, a key determinant of discriminatory capability is spatial localization of a minimum number of these steps to the engaged TCR. Biochemical and imaging experiments support these findings. Our results also reveal the discriminatory role of early negative feedback and necessary amplification conferred by late positive feedback.

David Chandler, a major figure in statistical mechanics, spent his career at the University of Illinois Urbana-Champaign, the University of Pennsylvania and the University of California Berkeley. Starting in his graduate work, he made significant advances in many areas of statistical mechanics theory, such as the structure and thermodynamics of simple liquids and nonpolar molecular liquids, the nature of hydrophobic hydration and hydrophobic interactions in aqueous systems, chemical reaction rates, quantum processes in liquids such as electron transfer and the solvation of an excess electron in water, ‘transition path sampling’ (a method for using computer simulations to study chemical reaction rates and other dynamic processes in liquids), and the ‘dynamic facilitation’ theory of the properties of supercooled liquids and the glass transition. He received the Hildebrand Award and the Theoretical Chemistry Award from the American Chemical Society and the Irving Langmuir Chemical Physics Award from the American Physical Society. He was elected to membership in the US National Academy of Sciences and the American Academy of Arts and Sciences.

The rise of SARS-CoV-2 variants and the history of outbreaks caused by zoonotic coronaviruses point to the need for next-generation vaccines that confer protection against variant strains. Here, we combined analyses of diverse sequences and structures of coronavirus spikes with data from deep mutational scanning to design SARS-CoV-2 variant antigens containing the most significant mutations that may emerge. We trained a neural network to predict RBD expression and ACE2 binding from sequence, which allowed us to determine that these antigens are stable and bind to ACE2. Thus, they represent viable variants. We then used a computational model of affinity maturation (AM) to study the antibody response to immunization with different combinations of the designed antigens. The results suggest that immunization with a cocktail of the antigens is likely to promote evolution of higher titers of antibodies that target SARS-CoV-2 variants than immunization or infection with the wildtype virus alone. Finally, our analysis of 12 coronaviruses from different genera identified the S2’ cleavage site and fusion peptide as potential pan-coronavirus vaccine targets.Author SummarySARS-CoV-2 variants have already emerged and future variants may pose greater threats to the efficacy of current vaccines. Rather than using a reactive approach to vaccine development that would lag behind the evolution of the virus, such as updating the sequence in the vaccine with a current variant, we sought to use a proactive approach that predicts some of the mutations that could arise that could evade current immune responses. Then, by including these mutations in a new vaccine antigen, we might be able to protect against those potential variants before they appear. Toward this end, we used various computational methods including sequence analysis and machine learning to design such antigens. We then used simulations of antibody development, and the results suggest that immunization with our designed antigens is likely to result in an antibody response that is better able to target SARS-CoV-2 variants than current vaccines. We also leveraged our sequence analysis to suggest that a particular site on the spike protein could serve as a useful target for a pan-coronavirus vaccine.

Macroscopic membraneless organelles containing RNA such as the nucleoli, germ granules, and the Cajal body have been known for decades. These biomolecular condensates are liquid-like bodies that can be formed by a phase transition. Recent evidence has revealed the presence of similar microscopic condensates associated with the transcription of genes. This brief article summarizes thoughts about the importance of condensates in the regulation of transcription and how RNA molecules, as components of such condensates, control the synthesis of RNA. Models and experimental data suggest that RNAs from enhancers facilitate the formation of a condensate that stabilizes the binding of transcription factors and accounts for a burst of transcription at the promoter. Termination of this burst is pictured as a nonequilibrium feedback loop where additional RNA destabilizes the condensate.

The design of vaccines against highly mutable pathogens, such as HIV and influenza, requires a detailed understanding of how the adaptive immune system responds to encountering multiple variant antigens (Ags). Here, we describe a multiscale model of B cell receptor (BCR) affinity maturation that employs actual BCR nucleotide sequences and treats BCR/Ag interactions in atomistic detail. We apply the model to simulate the maturation of a broadly neutralizing Ab (bnAb) against HIV. Starting from a germline precursor sequence of the VRC01 anti-HIV Ab, we simulate BCR evolution in response to different vaccination protocols and different Ags, which were previously designed by us. The simulation results provide qualitative guidelines for future vaccine design and reveal unique insights into bnAb evolution against the CD4 binding site of HIV. Our model makes possible direct comparisons of simulated BCR populations with results of deep sequencing data, which will be explored in future applications.Author SummaryVaccination has saved more lives than any other medical procedure, and the impending end of the COVID-19 pandemic is also due to the rapid development of highly efficacious vaccines. But, we do not have robust ways to develop vaccines against highly mutable pathogens. For example, there is no effective vaccine against HIV, and a universal vaccine against diverse strains of influenza is also not available. The development of immunization strategies to elicit antibodies that can neutralize diverse strains of highly mutable pathogens (so-called ‘broadly neutralizing antibodies’, or bnAbs) would enable the design of universal vaccines against such pathogens, as well as other viruses that may emerge in the future. In this paper, we present an agent-based model of affinity maturation – the Darwinian process by which antibodies evolve against a pathogen – that, for the first time, enables thein silicoinvestigation of real germline nucleotide sequences of antibodies known to evolve into potent bnAbs, evolving against real amino acid sequences of HIV-based vaccine-candidate proteins. Our results provide new insights into bnAb evolution against HIV, and can be used to qualitatively guide the future design of vaccines against highly mutable pathogens.

The design of vaccines against highly mutable pathogens, such as HIV and influenza, requires a detailed understanding of how the adaptive immune system responds to encountering multiple variant antigens (Ags). Here, we describe a multiscale model of B cell receptor (BCR) affinity maturation that employs actual BCR nucleotide sequences and treats BCR/Ag interactions in atomistic detail. We apply the model to simulate the maturation of a broadly neutralizing Ab (bnAb) against HIV. Starting from a germline precursor sequence of the VRC01 anti-HIV Ab, we simulate BCR evolution in response to different vaccination protocols and different Ags, which were previously designed by us. The simulation results provide qualitative guidelines for future vaccine design and reveal unique insights into bnAb evolution against the CD4 binding site of HIV. Our model makes possible direct comparisons of simulated BCR populations with results of deep sequencing data, which will be explored in future applications.

There still are no effective long-term protective vaccines against viruses that continuously evolve under immune pressure such as seasonal influenza, which has caused, and can cause, devastating epidemics in the human population. For finding such a broadly protective immunization strategy it is useful to know how easily the virus can escape via mutation from specific antibody responses. This information is encoded in the fitness landscape of the viral proteins (i.e., knowledge of the viral fitness as a function of sequence). Here we present a computational method to infer the intrinsic mutational fitness landscape of influenza-like evolving antigens from yearly sequence data. We test inference performance with computer-generated sequence data that are based on stochastic simulations mimicking basic features of immune-driven viral evolution. Although the numerically simulated model does create a phylogeny based on the allowed mutations, the inference scheme does not use this information. This provides a contrast to other methods that rely on reconstruction of phylogenetic trees. Our method just needs a sufficient number of samples over multiple years. With our method we are able to infer single-as well as pairwise mutational fitness effects from the simulated sequence time series for short antigenic proteins. Our fitness inference approach may have potential future use for design of immunization protocols by identifying intrinsically vulnerable immune target combinations on antigens that evolve under immune-driven selection. This approach may in the future be applied to influenza and other novel viruses such as SARS-CoV-2, which evolves and, like influenza, might continue to escape the natural and vaccine-mediated immune pressures.

A vaccine which is effective against the HIV virus is considered to be the best solution to the ongoing global HIV/AIDS epidemic. In the past thirty years, numerous attempts to develop an effective vaccine have been made with little or no success, due, in large part, to the high mutability of the virus. More recent studies showed that a vaccine able to elicit broadly neutralizing antibodies (bnAbs), that is, antibodies that can neutralize a high fraction of global virus variants, has promise to protect against HIV. Such a vaccine has been proposed to involve at least three separate stages: First, activate the appropriate precursor B cells; second, shepherd affinity maturation along pathways toward bnAbs; and, third, polish the Ab response to bind with high affinity to diverse HIV envelopes (Env). This final stage may require immunization with a mixture of Envs. In this paper, we set up a framework based on theory and modeling to design optimal panels of antigens to use in such a mixture. The designed antigens are characterized experimentally and are shown to be stable and to be recognized by known HIV antibodies.

Regulation of biological processes typically incorporate mechanisms that both initiate and terminate the process and, where understood, these mechanisms often involve feedback control. Regulation of transcription is a fundamental cellular process where the mechanisms involved in initiation have been studied extensively but those involved in arresting the process are poorly understood. Modeling of the potential roles of RNA in transcriptional control suggested a non-equilibrium feedback control mechanism wherein low levels of RNA promote condensates formed by electrostatic interactions whereas relatively high levels promote dissolution of these condensates. Evidence from both in vitro and in vivo experiments support the model that RNAs produced during early steps in transcription initiation stimulate condensate formation whereas the burst of RNAs produced during elongation stimulate condensate dissolution. We propose that transcriptional regulation incorporates a feedback mechanism whereby transcribed RNAs initially stimulate but then ultimately arrest the process.

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.

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 COVID-19 patients. 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 prior to SARS-CoV-2 infection., National Science Foundation (U.S.) (Grant PHY-2026995), Frederick National Laboratory for Cancer Research (Contract HHSN261200800001E), National Institutes of Health (U.S.) (Grant AI138790)

Effective prophylactic vaccines usually induce the immune system to generate potent antibodies that can bind to an antigen and thus prevent it from infecting host cells. B cells produce antibodies by a Darwinian evolutionary process called affinity maturation (AM). During AM, the B cell population evolves in response to the antigen to produce antibodies that bind specifically and strongly to the antigen. Highly mutable pathogens pose a major challenge to the development of effective vaccines because antibodies that are effective against one strain of the virus may not protect against a mutant strain. Antibodies that can protect against diverse strains of a mutable pathogen have high “breadth” and are called broadly neutralizing antibodies (bnAbs). In spite of extensive studies, an effective vaccination strategy that can generate bnAbs in humans does not exist for any highly mutable pathogen. Here we study a minimal model to explore the mechanisms underlying how the selection forces imposed by antigens can be optimally chosen to guide AM to maximize the evolution of bnAbs. For logistical reasons, only a finite number of antigens can be administered in a finite number of vaccinations; that is, guiding the non-equilibrium dynamics of AM to produce bnAbs must be accomplished non-adiabatically. The time-varying Kullback-Leibler divergence (KLD) between the existing B cell population distribution and the fitness landscape imposed by antigens is a quantitative metric of the thermodynamic force acting on B cells. If this force is too small, adaptation is minimal. If the force is too large, contrary to expectations, adaptation is not faster; rather, the B cell population is extinguished for reasons that we describe. We define the conditions necessary for the force to be set optimally such that the flux of B cells from low to high breadth states is maximized. Even in this case we show why the dynamics of AM prevent perfect adaptation. If two shots of vaccination are allowed, the optimal protocol is characterized by a relatively low optimal KLD during the first shot that appropriately increases the diversity of the B cell population so that the surviving B cells have a high chance of evolving into bnAbs upon subsequently increasing the KLD during the second shot. Phylogenetic tree analysis further reveals the evolutionary pathways that lead to bnAbs. The connections between the mechanisms revealed by our analyses and recent simulation studies of bnAb evolution, the problem of generalist versus specialist evolution, and learning theory are discussed.

Natural infections and vaccination with a pathogen typically stimulate the production of potent antibodies specific for the pathogen through a Darwinian evolutionary process known as affinity maturation. Such antibodies provide protection against reinfection by the same strain of a pathogen. A highly mutable virus, like HIV or influenza, evades recognition by these strain-specific antibodies via the emergence of new mutant strains. A vaccine that elicits antibodies that can bind to many diverse strains of the virus-known as broadly neutralizing antibodies (bnAbs)-could protect against highly mutable pathogens. Despite much work, the mechanisms by which bnAbs emerge remain uncertain. Using a computational model of affinity maturation, we studied a wide variety of vaccination strategies. Our results suggest that an effective strategy to maximize bnAb evolution is through a sequential immunization protocol, wherein each new immunization optimally increases the pressure on the immune system to target conserved antigenic sites, thus conferring breadth. We describe the mechanisms underlying why sequentially driving the immune system increasingly further from steady state, in an optimal fashion, is effective. The optimal protocol allows many evolving B cells to become bnAbs via diverse evolutionary paths.

Drug partitioning in nuclear condensates There is increasing interest in the function of phase-separated biomolecular condensates in cells because of their distinct properties and expanding roles in important biological processes. Klein et al. considered the fate of small-molecule therapeutics in the context of nuclear condensates (see the Perspective by Viny and Levine). They show that certain antineoplastic drugs have physicochemical properties that cause them to concentrate preferentially in condensates, both in vitro and in cancer cells. This property influences drug activity, and protein mutations that alter condensate formation can lead to drug resistance. Optimizing condensate partitioning may be valuable in developing improved therapeutics. Science , this issue p. 1386 ; see also p. 1314

Biological organisms experience constantly changing environments, from sudden changes in physiology brought about by feeding, to the regular rising and setting of the Sun, to ecological changes over evolutionary timescales. Living organisms have evolved to thrive in this changing world but the general principles by which organisms shape and are shaped by time varying environments remain elusive. Our understanding is particularly poor in the intermediate regime with no separation of timescales, where the environment changes on the same timescale as the physiological or evolutionary response. Experiments to systematically characterize the response to dynamic environments are challenging since such environments are inherently high dimensional. This roadmap deals with the unique role played by time varying environments in biological phenomena across scales, from physiology to evolution, seeking to emphasize the commonalities and the challenges faced in this emerging area of research.

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.

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., Poliovirus has a higher mutation rate than HIV, yet has been almost eradicated by vaccination while an effective vaccine against HIV does not exist. Here, the authors develop a fitness model for poliovirus viral protein 1 to show that it is subject to stringent evolutionary constraints that limit its ability to avoid vaccine-induced immune responses.

Significance An effective vaccine for HIV is still not available, although recent hope has emerged through the discovery of antibodies capable of neutralizing diverse HIV strains. Nonetheless, there exist mutational pathways through which HIV can evade known broadly neutralizing antibody responses. An ideal vaccine would elicit broadly neutralizing antibodies that target parts of the virus’s spike proteins where mutations severely compromise the virus’s fitness. Here, we employ a computational approach that allows estimation of the fitness landscape (fitness as a function of sequence) of the polyprotein that comprises HIV’s spike. We validate the inferred landscape through comparisons with diverse experimental measurements. The availability of this fitness landscape will aid the rational design of immunogens for effective vaccines., 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.

Drug-resistant mutations often have deleterious impacts on replication fitness, posing a fitness cost that can only be overcome by compensatory mutations. However, the role of fitness cost in the evolution of drug resistance has often been overlooked in clinical studies or in vitro selection experiments, as these observations only capture the outcome of drug selection. In this study, we systematically profile the fitness landscape of resistance-associated sites in HIV-1 protease using deep mutational scanning. We construct a mutant library covering combinations of mutations at 11 sites in HIV-1 protease, all of which are associated with resistance to protease inhibitors in clinic. Using deep sequencing, we quantify the fitness of thousands of HIV-1 protease mutants after multiple cycles of replication in human T cells. Although the majority of resistance-associated mutations have deleterious effects on viral replication, we find that epistasis among resistance-associated mutations is predominantly positive. Furthermore, our fitness data are consistent with genetic interactions inferred directly from HIV sequence data of patients. Fitness valleys formed by strong positive epistasis reduce the likelihood of reversal of drug resistance mutations. Overall, our results support the view that strong compensatory effects are involved in the emergence of clinically observed resistance mutations and provide insights to understanding fitness barriers in the evolution and reversion of drug resistance., Author summary Antiretroviral drugs have achieved great success in controlling the HIV pandemic. However, the therapy fails sometimes owing to the low drug adherence and/or the emergence of resistance associated mutations on viral genome. The persistence of drug resistance poses challenges in using antiretroviral drugs for long term control or pre-exposure prophylaxis. To understand the mechanisms of resistance evolution and persistence, we profiled the replication fitness of over 1000 HIV-1 mutants with combinations of resistance associated mutations on its protease gene. We found that although resistance associated mutations greatly reduce replication fitness, they interact positively to alleviate the mutational load. These genetic interactions, termed epistasis, increase the ruggedness along the evolution paths, restricting resistance associated mutations from reversal. Our data support the clinical observations that drug resistance mutations tend to persist even when antiretroviral drug is discontinued.

Copyright © 2019 The Authors 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.

An effective vaccine is urgently required to curb the HIV-1 epidemic. We have previously described an approach to model the fitness landscape of several HIV-1 proteins, and have validated the results against experimental and clinical data. The fitness landscape may be used to identify mutation patterns harmful to virus viability, and consequently inform the design of immunogens that can target such regions for immunological control. Here we apply such an analysis and complementary experiments to HIV-1 Nef, a multifunctional protein which plays a key role in HIV-1 pathogenesis. We measured Nef-driven replication capacities as well as Nef-mediated CD4 and HLA-I down-modulation capacities of thirty-two different Nef mutants, and tested model predictions against these results. Furthermore, we evaluated the models using 448 patient-derived Nef sequences for which several Nef activities were previously measured. Model predictions correlated significantly with Nef-driven replication and CD4 down-modulation capacities, but not HLA-I down-modulation capacities, of the various Nef mutants. Similarly, in our analysis of patient-derived Nef sequences, CD4 down-modulation capacity correlated the most significantly with model predictions, suggesting that of the tested Nef functions, this is the most important in vivo. Overall, our results highlight how the fitness landscape inferred from patient-derived sequences captures, at least in part, the in vivo functional effects of mutations to Nef. However, the correlation between predictions of the fitness landscape and measured parameters of Nef function is not as accurate as the correlation observed in past studies for other proteins. This may be because of the additional complexity associated with inferring the cost of mutations on the diverse functions of Nef.

An effective vaccine is urgently required to curb the HIV-1 epidemic. We have previously described an approach to model the fitness landscape of several HIV-1 proteins, and have validated the results against experimental and clinical data. The fitness landscape may be used to identify mutation patterns harmful to virus viability, and consequently inform the design of immunogens that can target such regions for immunological control. Here we apply such an analysis and complementary experiments to HIV-1 Nef, a multifunctional protein which plays a key role in HIV-1 pathogenesis. We measured Nef-driven replication capacities as well as Nef-mediated CD4 and HLA-I down-modulation capacities of thirty-two different Nef mutants, and tested model predictions against these results. Furthermore, we evaluated the models using 448 patient-derived Nef sequences for which several Nef activities were previously measured. Model predictions correlated significantly with Nef-driven replication and CD4 down-modulation capacities, but not HLA-I down-modulation capacities, of the various Nef mutants. Similarly, in our analysis of patient-derived Nef sequences, CD4 down-modulation capacity correlated the most significantly with model predictions, suggesting that of the tested Nef functions, this is the most important in vivo. Overall, our results highlight how the fitness landscape inferred from patient-derived sequences captures, at least in part, the in vivo functional effects of mutations to Nef. However, the correlation between predictions of the fitness landscape and measured parameters of Nef function is not as accurate as the correlation observed in past studies for other proteins. This may be because of the additional complexity associated with inferring the cost of mutations on the diverse functions of Nef.

3131 Background: The molecules of the cell are compartmentalized into membrane- and non-membrane-bound organelles. Many non-membrane-bound organelles are phase-separated biomolecular condensates with distinct physicochemical properties that can absorb and concentrate specific proteins and nucleic acids involved in discrete biochemical processes. We reasoned that selective condensate partitioning might also occur with small molecule drugs whose targets occur within condensates, and that the therapeutic index and efficacy of such compounds might therefore relate to their ability to partition into condensates. Methods: To study the behavior of drugs within condensates, these were modeled in cells and in vitro with purified proteins and visualized by fluorescent confocal microscopy. The functional outcomes of condensate partitioning were queried in cells. Results: We found that cisplatin, tamoxifen, JQ1, THZ1, and mitoxantrone are concentrated in specific protein condensates in vitro, and that this occurs through physicochemical properties independent of the drug target. A screen of a chemically diverse fluorescent probes and mutant-protein condensates demonstrated that pi-system interactions between aromatic moieties in the protein and small molecule govern concentration in condensates. These results show that clinically important drugs partition into specific protein condensates in vitro by virtue of defined chemical properties, thereby altering their local concentration. In vitro droplet assays revealed that cisplatin is selectively concentrated in transcriptional condensates, and that this ability is required for efficient platination of target DNA. In cell studies revealed that cisplatin preferentially targets DNA contained within MED1 condensates, and disrupts the genetic regulatory elements that compose phase-separated transcriptional condensates. Live cell imaging demonstrated that transcriptional condensates are dissolved by cisplatin, whereas other condensates remain intact. Conclusions: Our results show that antineoplastic drugs partition selectively into biomolecular condensates, that this can occur through physicochemical properties independent of their molecular targets, and that drug activity may occur through condensate-related mechanisms. These results have implications for development of efficacious cancer therapeutics; effective target engagement will depend on factors such as drug partitioning in condensates. Assays of the type described here may thus help optimize condensate partitioning, target engagement, and the therapeutic index of drugs for cancer treatment.

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.

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., National Institute of Allergy and Infectious Diseases (U.S.) (Grant UM1-AI100663)