1
|
Klink GV, Kalinina OV, Bazykin GA. Changing selection on amino acid substitutions in Gag protein between major HIV-1 subtypes. Virus Evol 2024; 10:veae036. [PMID: 38808036 PMCID: PMC11131029 DOI: 10.1093/ve/veae036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2023] [Revised: 12/27/2023] [Accepted: 04/28/2024] [Indexed: 05/30/2024] Open
Abstract
Amino acid preferences at a protein site depend on the role of this site in protein function and structure as well as on external constraints. All these factors can change in the course of evolution, making amino acid propensities of a site time-dependent. When viral subtypes divergently evolve in different host subpopulations, such changes may depend on genetic, medical, and sociocultural differences between these subpopulations. Here, using our previously developed phylogenetic approach, we describe sixty-nine amino acid sites of the Gag protein of human immunodeficiency virus type 1 (HIV-1) where amino acids have different impact on viral fitness in six major subtypes of the type M. These changes in preferences trigger adaptive evolution; indeed, 32 (46 per cent) of these sites experienced strong positive selection at least in one of the subtypes. At some of the sites, changes in amino acid preferences may be associated with differences in immune escape between subtypes. The prevalence of an amino acid in a protein site within a subtype is only a poor predictor for whether this amino acid is preferred in this subtype according to the phylogenetic analysis. Therefore, attempts to identify the factors of viral evolution from comparative genomics data should integrate across multiple sources of information.
Collapse
Affiliation(s)
- Galya V Klink
- Laboratory of Molecular Evolution, Institute for Information Transmission Problems (Kharkevich Institute) of the Russian Academy of Sciences, Bolshoy Karetny per. 19, build.1, Moscow 127051, Russia
- Center for Molecular and Cellular Biology, Skolkovo Institute of Science and Technology, Bolshoy Boulevard, 30, p.1, Skolkovo 121205, Russia
| | - Olga V Kalinina
- Drug Bioinformatics, Helmholtz Institute for Pharmaceutical Research Saarland (HIPS)/Helmholtz Centre for Infection Research (HZI), Campus E8.1, Saarbrücken 66123, Germany
- Center for Bioinformatics, Saarland University, Campus E2.1, Saarbrücken 66123, Germany
- Medical Faculty, Saarland University, Kirrberger Str. 100, Homburg 66421, Germany
| | - Georgii A Bazykin
- Laboratory of Molecular Evolution, Institute for Information Transmission Problems (Kharkevich Institute) of the Russian Academy of Sciences, Bolshoy Karetny per. 19, build.1, Moscow 127051, Russia
| |
Collapse
|
2
|
Meger AT, Spence MA, Sandhu M, Matthews D, Chen J, Jackson CJ, Raman S. Rugged fitness landscapes minimize promiscuity in the evolution of transcriptional repressors. Cell Syst 2024; 15:374-387.e6. [PMID: 38537640 DOI: 10.1016/j.cels.2024.03.002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2023] [Revised: 09/08/2023] [Accepted: 03/05/2024] [Indexed: 04/20/2024]
Abstract
How a protein's function influences the shape of its fitness landscape, smooth or rugged, is a fundamental question in evolutionary biochemistry. Smooth landscapes arise when incremental mutational steps lead to a progressive change in function, as commonly seen in enzymes and binding proteins. On the other hand, rugged landscapes are poorly understood because of the inherent unpredictability of how sequence changes affect function. Here, we experimentally characterize the entire sequence phylogeny, comprising 1,158 extant and ancestral sequences, of the DNA-binding domain (DBD) of the LacI/GalR transcriptional repressor family. Our analysis revealed an extremely rugged landscape with rapid switching of specificity, even between adjacent nodes. Further, the ruggedness arises due to the necessity of the repressor to simultaneously evolve specificity for asymmetric operators and disfavors potentially adverse regulatory crosstalk. Our study provides fundamental insight into evolutionary, molecular, and biophysical rules of genetic regulation through the lens of fitness landscapes.
Collapse
Affiliation(s)
- Anthony T Meger
- Department of Biochemistry, University of Wisconsin-Madison, Madison, WI 53706, USA
| | - Matthew A Spence
- Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia
| | - Mahakaran Sandhu
- Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia
| | - Dana Matthews
- Research School of Biology, Australian National University, Canberra, ACT 2601, Australia
| | - Jackie Chen
- Department of Biochemistry, University of Wisconsin-Madison, Madison, WI 53706, USA
| | - Colin J Jackson
- Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia; ARC Centre of Excellence for Innovations in Peptide & Protein Science, Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia; ARC Centre of Excellence for Innovations in Synthetic Biology, Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia.
| | - Srivatsan Raman
- Department of Biochemistry, University of Wisconsin-Madison, Madison, WI 53706, USA; Department of Bacteriology, University of Wisconsin-Madison, Madison, WI 53706, USA; Department of Chemical and Biological Engineering, University of Wisconsin-Madison, Madison, WI 53706, USA.
| |
Collapse
|
3
|
Gao Y, Barton JP. A binary trait model reveals the fitness effects of HIV-1 escape from T cell responses. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.03.03.583183. [PMID: 38464239 PMCID: PMC10925374 DOI: 10.1101/2024.03.03.583183] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/12/2024]
Abstract
Natural selection often acts on multiple traits simultaneously. For example, the virus HIV-1 faces pressure to evade host immunity while also preserving replicative fitness. While past work has studied selection during HIV-1 evolution, it is challenging to quantitatively separate different contributions to fitness. This task is made more difficult because a single mutation can affect both immune escape and replication. Here, we develop an evolutionary model that disentangles the effects of escaping CD8+ T cell-mediated immunity, which we model as a binary trait, from other contributions to fitness. After validation in simulations, we applied this model to study within-host HIV-1 evolution in a clinical data set. We observed strong selection for immune escape, sometimes greatly exceeding past estimates, especially early in infection. Conservative estimates suggest that roughly half of HIV-1 fitness gains during the first months to years of infection can be attributed to T cell escape. Our approach is not limited to HIV-1 or viruses, and could be adapted to study the evolution of quantitative traits in other contexts.
Collapse
Affiliation(s)
- Yirui Gao
- Department of Physics and Astronomy, University of California, Riverside, USA
| | - John P. Barton
- Department of Physics and Astronomy, University of California, Riverside, USA
- Department of Computational and Systems Biology, University of Pittsburgh School of Medicine, USA
- Department of Physics and Astronomy, University of Pittsburgh, USA
| |
Collapse
|
4
|
Wang Y, Stebe KJ, de la Fuente-Nunez C, Radhakrishnan R. Computational Design of Peptides for Biomaterials Applications. ACS APPLIED BIO MATERIALS 2024; 7:617-625. [PMID: 36971822 DOI: 10.1021/acsabm.2c01023] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/29/2023]
Abstract
Computer-aided molecular design and protein engineering emerge as promising and active subjects in bioengineering and biotechnological applications. On one hand, due to the advancing computing power in the past decade, modeling toolkits and force fields have been put to use for accurate multiscale modeling of biomolecules including lipid, protein, carbohydrate, and nucleic acids. On the other hand, machine learning emerges as a revolutionary data analysis tool that promises to leverage physicochemical properties and structural information obtained from modeling in order to build quantitative protein structure-function relationships. We review recent computational works that utilize state-of-the-art computational methods to engineer peptides and proteins for various emerging biomedical, antimicrobial, and antifreeze applications. We also discuss challenges and possible future directions toward developing a roadmap for efficient biomolecular design and engineering.
Collapse
Affiliation(s)
- Yiming Wang
- Department of Bioengineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
- Penn Institute for Computational Science, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Kathleen J Stebe
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Cesar de la Fuente-Nunez
- Department of Bioengineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
- Penn Institute for Computational Science, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
- Machine Biology Group, Department of Psychiatry and Microbiology, Institute for Biomedical Informatics, Institute for Translational Medicine and Therapeutics, Perelman School of Medicine University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Ravi Radhakrishnan
- Department of Bioengineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
- Penn Institute for Computational Science, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
5
|
Alvarez S, Nartey CM, Mercado N, de la Paz JA, Huseinbegovic T, Morcos F. In vivo functional phenotypes from a computational epistatic model of evolution. Proc Natl Acad Sci U S A 2024; 121:e2308895121. [PMID: 38285950 PMCID: PMC10861889 DOI: 10.1073/pnas.2308895121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2023] [Accepted: 12/19/2023] [Indexed: 01/31/2024] Open
Abstract
Computational models of evolution are valuable for understanding the dynamics of sequence variation, to infer phylogenetic relationships or potential evolutionary pathways and for biomedical and industrial applications. Despite these benefits, few have validated their propensities to generate outputs with in vivo functionality, which would enhance their value as accurate and interpretable evolutionary algorithms. We demonstrate the power of epistasis inferred from natural protein families to evolve sequence variants in an algorithm we developed called sequence evolution with epistatic contributions (SEEC). Utilizing the Hamiltonian of the joint probability of sequences in the family as fitness metric, we sampled and experimentally tested for in vivo [Formula: see text]-lactamase activity in Escherichia coli TEM-1 variants. These evolved proteins can have dozens of mutations dispersed across the structure while preserving sites essential for both catalysis and interactions. Remarkably, these variants retain family-like functionality while being more active than their wild-type predecessor. We found that depending on the inference method used to generate the epistatic constraints, different parameters simulate diverse selection strengths. Under weaker selection, local Hamiltonian fluctuations reliably predict relative changes to variant fitness, recapitulating neutral evolution. SEEC has the potential to explore the dynamics of neofunctionalization, characterize viral fitness landscapes, and facilitate vaccine development.
Collapse
Affiliation(s)
- Sophia Alvarez
- Department of Biological Sciences, University of Texas at Dallas, Richardson, TX75080
| | - Charisse M. Nartey
- Department of Biological Sciences, University of Texas at Dallas, Richardson, TX75080
| | - Nicholas Mercado
- Department of Biological Sciences, University of Texas at Dallas, Richardson, TX75080
| | | | - Tea Huseinbegovic
- Department of Biological Sciences, University of Texas at Dallas, Richardson, TX75080
| | - Faruck Morcos
- Department of Biological Sciences, University of Texas at Dallas, Richardson, TX75080
- Department of Bioengineering, University of Texas at Dallas, Richardson, TX75080
- Center for Systems Biology, University of Texas at Dallas, Richardson, TX75080
| |
Collapse
|
6
|
Rajendra D, Maroli N, Dixit NM, Maiti PK. Molecular dynamics simulations show how antibodies may rescue HIV-1 mutants incapable of infecting host cells. J Biomol Struct Dyn 2023:1-11. [PMID: 38111161 DOI: 10.1080/07391102.2023.2294835] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2023] [Accepted: 12/01/2023] [Indexed: 12/20/2023]
Abstract
High mutation and replication rates of HIV-1 result in the continuous generation of variants, allowing it to adapt to changing host environments. Mutations often have deleterious effects, but variants carrying them are rapidly purged. Surprisingly, a particular variant incapable of entering host cells was found to be rescued by host antibodies targeting HIV-1. Understanding the molecular mechanism of this rescue is important to develop and improve antibody-based therapies. To unravel the underlying mechanisms, we performed fully atomistic molecular dynamics simulations of the HIV-1 gp41 trimer responsible for viral entry into host cells, its entry-deficient variant, and its complex with the rescuing antibody. We find that the Q563R mutation, which the entry-deficient variant carries, prevents the native conformation of the gp41 6-helix bundle required for entry and stabilizes an alternative conformation instead. This is the consequence of substantial changes in the secondary structure and interactions between the domains of gp41. Binding of the antibody F240 to gp41 reverses these changes and re-establishes the native conformation, resulting in rescue. To test the generality of this mechanism, we performed simulations with the entry-deficient L565A variant and antibody 3D6. We find that 3D6 binding was able to reverse structural and interaction changes introduced by the mutation and restore the native gp41 conformation. Viral variants may not only escape antibodies but be aided by them in their survival, potentially compromising antibody-based therapies, including vaccination and passive immunization. Our simulation framework could serve as a tool to assess the likelihood of such resistance against specific antibodies.Communicated by Ramaswamy H. SarmaCommunicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Dharanish Rajendra
- Centre for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bengaluru, India
| | - Nikhil Maroli
- Centre for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bengaluru, India
| | - Narendra M Dixit
- Department of Chemical Engineering, Indian Institute of Science, Bengaluru, India
- Department of Bioengineering, Indian Institute of Science, Bengaluru, India
| | - Prabal K Maiti
- Centre for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bengaluru, India
| |
Collapse
|
7
|
Zhang H, Bull RA, Quadeer AA, McKay MR. HCV E1 influences the fitness landscape of E2 and may enhance escape from E2-specific antibodies. Virus Evol 2023; 9:vead068. [PMID: 38107333 PMCID: PMC10722114 DOI: 10.1093/ve/vead068] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2023] [Revised: 09/27/2023] [Accepted: 11/16/2023] [Indexed: 12/19/2023] Open
Abstract
The Hepatitis C virus (HCV) envelope glycoprotein E1 forms a non-covalent heterodimer with E2, the main target of neutralizing antibodies. How E1-E2 interactions influence viral fitness and contribute to resistance to E2-specific antibodies remain largely unknown. We investigate this problem using a combination of fitness landscape and evolutionary modeling. Our analysis indicates that E1 and E2 proteins collectively mediate viral fitness and suggests that fitness-compensating E1 mutations may accelerate escape from E2-targeting antibodies. Our analysis also identifies a set of E2-specific human monoclonal antibodies that are predicted to be especially resilient to escape via genetic variation in both E1 and E2, providing directions for robust HCV vaccine development.
Collapse
Affiliation(s)
- Hang Zhang
- Department of Electronic and Computer Engineering, The Hong Kong University of Science and Technology, Clear Water Bay, Hong Kong, SAR, China
| | - Rowena A Bull
- School of Biomedical Sciences, Faculty of Medicine and Health, University of New South Wales, Sydney, NSW 2052, Australia
- The Kirby Institute for Infection and Immunity, Sydney, NSW 2052, Australia
| | - Ahmed Abdul Quadeer
- Department of Electronic and Computer Engineering, The Hong Kong University of Science and Technology, Clear Water Bay, Hong Kong, SAR, China
- Department of Electrical and Electronic Engineering, University of Melbourne, Parkville, VIC 3010, Australia
| | - Matthew R McKay
- Department of Electrical and Electronic Engineering, University of Melbourne, Parkville, VIC 3010, Australia
- Department of Microbiology and Immunology, The Peter Doherty Institute for Infection and Immunity, University of Melbourne, Melbourne, VIC 3000, Australia
| |
Collapse
|
8
|
Zhang H, Quadeer AA, McKay MR. Direct-acting antiviral resistance of Hepatitis C virus is promoted by epistasis. Nat Commun 2023; 14:7457. [PMID: 37978179 PMCID: PMC10656532 DOI: 10.1038/s41467-023-42550-6] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2023] [Accepted: 10/13/2023] [Indexed: 11/19/2023] Open
Abstract
Direct-acting antiviral agents (DAAs) provide efficacious therapeutic treatments for chronic Hepatitis C virus (HCV) infection. However, emergence of drug resistance mutations (DRMs) can greatly affect treatment outcomes and impede virological cure. While multiple DRMs have been observed for all currently used DAAs, the evolutionary determinants of such mutations are not currently well understood. Here, by considering DAAs targeting the nonstructural 3 (NS3) protein of HCV, we present results suggesting that epistasis plays an important role in the evolution of DRMs. Employing a sequence-based fitness landscape model whose predictions correlate highly with experimental data, we identify specific DRMs that are associated with strong epistatic interactions, and these are found to be enriched in multiple NS3-specific DAAs. Evolutionary modelling further supports that the identified DRMs involve compensatory mutational interactions that facilitate relatively easy escape from drug-induced selection pressures. Our results indicate that accounting for epistasis is important for designing future HCV NS3-targeting DAAs.
Collapse
Affiliation(s)
- Hang Zhang
- Department of Electronic and Computer Engineering, The Hong Kong University of Science and Technology, Clear Water Bay, Hong Kong SAR, China
| | - Ahmed Abdul Quadeer
- Department of Electronic and Computer Engineering, The Hong Kong University of Science and Technology, Clear Water Bay, Hong Kong SAR, China.
| | - Matthew R McKay
- Department of Electrical and Electronic Engineering, University of Melbourne, Melbourne, VIC, Australia.
- Department of Microbiology and Immunology, University of Melbourne, at The Peter Doherty Institute for Infection and Immunity, Melbourne, VIC, Australia.
| |
Collapse
|
9
|
Li M, Oliveira Passos D, Shan Z, Smith SJ, Sun Q, Biswas A, Choudhuri I, Strutzenberg TS, Haldane A, Deng N, Li Z, Zhao XZ, Briganti L, Kvaratskhelia M, Burke TR, Levy RM, Hughes SH, Craigie R, Lyumkis D. Mechanisms of HIV-1 integrase resistance to dolutegravir and potent inhibition of drug-resistant variants. SCIENCE ADVANCES 2023; 9:eadg5953. [PMID: 37478179 DOI: 10.1126/sciadv.adg5953] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/09/2023] [Accepted: 06/16/2023] [Indexed: 07/23/2023]
Abstract
HIV-1 infection depends on the integration of viral DNA into host chromatin. Integration is mediated by the viral enzyme integrase and is blocked by integrase strand transfer inhibitors (INSTIs), first-line antiretroviral therapeutics widely used in the clinic. Resistance to even the best INSTIs is a problem, and the mechanisms of resistance are poorly understood. Here, we analyze combinations of the mutations E138K, G140A/S, and Q148H/K/R, which confer resistance to INSTIs. The investigational drug 4d more effectively inhibited the mutants compared with the approved drug Dolutegravir (DTG). We present 11 new cryo-EM structures of drug-resistant HIV-1 intasomes bound to DTG or 4d, with better than 3-Å resolution. These structures, complemented with free energy simulations, virology, and enzymology, explain the mechanisms of DTG resistance involving E138K + G140A/S + Q148H/K/R and show why 4d maintains potency better than DTG. These data establish a foundation for further development of INSTIs that potently inhibit resistant forms in integrase.
Collapse
Affiliation(s)
- Min Li
- National Institute of Diabetes and Digestive Diseases, National Institutes of Health, Bethesda, MD, 20892, USA
| | | | - Zelin Shan
- The Salk Institute for Biological Studies, La Jolla, CA, 92037, USA
| | - Steven J Smith
- Center for Cancer Research, National Cancer Institute, National Institutes of Health, Frederick, MD, 21702, USA
| | - Qinfang Sun
- Center for Biophysics and Computational Biology, and Department of Chemistry, Temple University, Philadelphia, PA 19122, USA
| | - Avik Biswas
- The Salk Institute for Biological Studies, La Jolla, CA, 92037, USA
- Center for Biophysics and Computational Biology and Department of Physics, Temple University, Philadelphia, PA 19122, USA
| | - Indrani Choudhuri
- Center for Biophysics and Computational Biology, and Department of Chemistry, Temple University, Philadelphia, PA 19122, USA
| | | | - Allan Haldane
- Center for Biophysics and Computational Biology and Department of Physics, Temple University, Philadelphia, PA 19122, USA
| | - Nanjie Deng
- Department of Chemistry and Physical Sciences, Pace University, New York, NY, 10038, USA
| | - Zhaoyang Li
- National Institute of Diabetes and Digestive Diseases, National Institutes of Health, Bethesda, MD, 20892, USA
| | - Xue Zhi Zhao
- Center for Cancer Research, National Cancer Institute, National Institutes of Health, Frederick, MD, 21702, USA
| | - Lorenzo Briganti
- Division of Infectious Diseases, University of Colorado Anschutz Medical Campus, Aurora, CO 80045, USA
| | - Mamuka Kvaratskhelia
- Division of Infectious Diseases, University of Colorado Anschutz Medical Campus, Aurora, CO 80045, USA
| | - Terrence R Burke
- Center for Cancer Research, National Cancer Institute, National Institutes of Health, Frederick, MD, 21702, USA
| | - Ronald M Levy
- Center for Biophysics and Computational Biology and Department of Physics, Temple University, Philadelphia, PA 19122, USA
| | - Stephen H Hughes
- Center for Cancer Research, National Cancer Institute, National Institutes of Health, Frederick, MD, 21702, USA
| | - Robert Craigie
- National Institute of Diabetes and Digestive Diseases, National Institutes of Health, Bethesda, MD, 20892, USA
| | - Dmitry Lyumkis
- The Salk Institute for Biological Studies, La Jolla, CA, 92037, USA
- Department of Integrative Structural and Computational Biology, The Scripps Research Institute, La Jolla, CA, 92037, USA
- Graduate School of Biological Sciences, Section of Molecular Biology, University of California San Diego, La Jolla, CA 92093, USA
| |
Collapse
|
10
|
Fili M, Hu G, Han C, Kort A, Trettin J, Haim H. A classification algorithm based on dynamic ensemble selection to predict mutational patterns of the envelope protein in HIV-infected patients. Algorithms Mol Biol 2023; 18:4. [PMID: 37337202 DOI: 10.1186/s13015-023-00228-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2022] [Accepted: 06/04/2023] [Indexed: 06/21/2023] Open
Abstract
BACKGROUND Therapeutics against the envelope (Env) proteins of human immunodeficiency virus type 1 (HIV-1) effectively reduce viral loads in patients. However, due to mutations, new therapy-resistant Env variants frequently emerge. The sites of mutations on Env that appear in each patient are considered random and unpredictable. Here we developed an algorithm to estimate for each patient the mutational state of each position based on the mutational state of adjacent positions on the three-dimensional structure of the protein. METHODS We developed a dynamic ensemble selection algorithm designated k-best classifiers. It identifies the best classifiers within the neighborhood of a new observation and applies them to predict the variability state of each observation. To evaluate the algorithm, we applied amino acid sequences of Envs from 300 HIV-1-infected individuals (at least six sequences per patient). For each patient, amino acid variability values at all Env positions were mapped onto the three-dimensional structure of the protein. Then, the variability state of each position was estimated by the variability at adjacent positions of the protein. RESULTS The proposed algorithm showed higher performance than the base learner and a panel of classification algorithms. The mutational state of positions in the high-mannose patch and CD4-binding site of Env, which are targeted by multiple therapeutics, was predicted well. Importantly, the algorithm outperformed other classification techniques for predicting the variability state at multi-position footprints of therapeutics on Env. CONCLUSIONS The proposed algorithm applies a dynamic classifier-scoring approach that increases its performance relative to other classification methods. Better understanding of the spatiotemporal patterns of variability across Env may lead to new treatment strategies that are tailored to the unique mutational patterns of each patient. More generally, we propose the algorithm as a new high-performance dynamic ensemble selection technique.
Collapse
Affiliation(s)
- Mohammad Fili
- Department of Industrial and Manufacturing Systems Engineering, Iowa State University, 3014 Black Engineering, 2529 Union Drive, Ames, IA, 50011, USA
| | - Guiping Hu
- Department of Industrial and Manufacturing Systems Engineering, Iowa State University, 3014 Black Engineering, 2529 Union Drive, Ames, IA, 50011, USA.
| | - Changze Han
- Department of Microbiology and Immunology, Carver College of Medicine, University of Iowa, 51 Newton Rd, 3-770 BSB, Iowa City, IA, 52242, USA
| | - Alexa Kort
- Department of Microbiology and Immunology, Carver College of Medicine, University of Iowa, 51 Newton Rd, 3-770 BSB, Iowa City, IA, 52242, USA
| | - John Trettin
- Department of Industrial and Manufacturing Systems Engineering, Iowa State University, 3014 Black Engineering, 2529 Union Drive, Ames, IA, 50011, USA
| | - Hillel Haim
- Department of Microbiology and Immunology, Carver College of Medicine, University of Iowa, 51 Newton Rd, 3-770 BSB, Iowa City, IA, 52242, USA.
| |
Collapse
|
11
|
Alvarez S, Nartey CM, Mercado N, de la Paz A, Huseinbegovic T, Morcos F. In vivo functional phenotypes from a computational epistatic model of evolution. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.05.24.542176. [PMID: 37292895 PMCID: PMC10245989 DOI: 10.1101/2023.05.24.542176] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Computational models of evolution are valuable for understanding the dynamics of sequence variation, to infer phylogenetic relationships or potential evolutionary pathways and for biomedical and industrial applications. Despite these benefits, few have validated their propensities to generate outputs with in vivo functionality, which would enhance their value as accurate and interpretable evolutionary algorithms. We demonstrate the power of epistasis inferred from natural protein families to evolve sequence variants in an algorithm we developed called Sequence Evolution with Epistatic Contributions. Utilizing the Hamiltonian of the joint probability of sequences in the family as fitness metric, we sampled and experimentally tested for in vivo β -lactamase activity in E. coli TEM-1 variants. These evolved proteins can have dozens of mutations dispersed across the structure while preserving sites essential for both catalysis and interactions. Remarkably, these variants retain family-like functionality while being more active than their WT predecessor. We found that depending on the inference method used to generate the epistatic constraints, different parameters simulate diverse selection strengths. Under weaker selection, local Hamiltonian fluctuations reliably predict relative changes to variant fitness, recapitulating neutral evolution. SEEC has the potential to explore the dynamics of neofunctionalization, characterize viral fitness landscapes and facilitate vaccine development.
Collapse
Affiliation(s)
- Sophia Alvarez
- Department of Biological Sciences, University of Texas at Dallas, Richardson, TX 75080, USA
| | - Charisse M. Nartey
- Department of Biological Sciences, University of Texas at Dallas, Richardson, TX 75080, USA
| | - Nicholas Mercado
- Department of Biological Sciences, University of Texas at Dallas, Richardson, TX 75080, USA
| | - Alberto de la Paz
- Department of Biological Sciences, University of Texas at Dallas, Richardson, TX 75080, USA
| | - Tea Huseinbegovic
- School of Natural Sciences and Mathematics, University of Texas at Dallas, Richardson, TX 75080, USA
| | - Faruck Morcos
- Department of Biological Sciences, University of Texas at Dallas, Richardson, TX 75080, USA
- Department of Bioengineering, University of Texas at Dallas, Richardson, TX 75080, USA
- Center for Systems Biology, University of Texas at Dallas, Richardson, TX 75080, USA
| |
Collapse
|
12
|
Dichio V, Zeng HL, Aurell E. Statistical genetics in and out of quasi-linkage equilibrium. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2023; 86:052601. [PMID: 36944245 DOI: 10.1088/1361-6633/acc5fa] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/08/2022] [Accepted: 03/21/2023] [Indexed: 06/18/2023]
Abstract
This review is about statistical genetics, an interdisciplinary topic between statistical physics and population biology. The focus is on the phase ofquasi-linkage equilibrium(QLE). Our goals here are to clarify under which conditions the QLE phase can be expected to hold in population biology and how the stability of the QLE phase is lost. The QLE state, which has many similarities to a thermal equilibrium state in statistical mechanics, was discovered by M Kimura for a two-locus two-allele model, and was extended and generalized to the global genome scale byNeher&Shraiman (2011). What we will refer to as the Kimura-Neher-Shraiman theory describes a population evolving due to the mutations, recombination, natural selection and possibly genetic drift. A QLE phase exists at sufficiently high recombination rate (r) and/or mutation ratesµwith respect to selection strength. We show how in QLE it is possible to infer the epistatic parameters of the fitness function from the knowledge of the (dynamical) distribution of genotypes in a population. We further consider the breakdown of the QLE regime for high enough selection strength. We review recent results for the selection-mutation and selection-recombination dynamics. Finally, we identify and characterize a new phase which we call the non-random coexistence where variability persists in the population without either fixating or disappearing.
Collapse
Affiliation(s)
- Vito Dichio
- Sorbonne Université, Paris Brain Institute-ICM, CNRS, Inria, Inserm, AP-HP, Hôpital de la Pitié Salpêtrière, F-75013 Paris, France
| | - Hong-Li Zeng
- School of Science, Nanjing University of Posts and Telecommunications, New Energy Technology Engineering Laboratory of Jiangsu Province, Nanjing 210023, People's Republic of China
| | - Erik Aurell
- Department of Computational Science and Technology, KTH-Royal Institute of Technology, AlbaNova University Center, SE-106 91 Stockholm, Sweden
| |
Collapse
|
13
|
Poulsen VM, DeDeo S. Inferring Cultural Landscapes with the Inverse Ising Model. ENTROPY (BASEL, SWITZERLAND) 2023; 25:264. [PMID: 36832631 PMCID: PMC9955041 DOI: 10.3390/e25020264] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/13/2022] [Revised: 01/22/2023] [Accepted: 01/25/2023] [Indexed: 06/18/2023]
Abstract
The space of possible human cultures is vast, but some cultural configurations are more consistent with cognitive and social constraints than others. This leads to a "landscape" of possibilities that our species has explored over millennia of cultural evolution. However, what does this fitness landscape, which constrains and guides cultural evolution, look like? The machine-learning algorithms that can answer these questions are typically developed for large-scale datasets. Applications to the sparse, inconsistent, and incomplete data found in the historical record have received less attention, and standard recommendations can lead to bias against marginalized, under-studied, or minority cultures. We show how to adapt the minimum probability flow algorithm and the Inverse Ising model, a physics-inspired workhorse of machine learning, to the challenge. A series of natural extensions-including dynamical estimation of missing data, and cross-validation with regularization-enables reliable reconstruction of the underlying constraints. We demonstrate our methods on a curated subset of the Database of Religious History: records from 407 religious groups throughout human history, ranging from the Bronze Age to the present day. This reveals a complex, rugged, landscape, with both sharp, well-defined peaks where state-endorsed religions tend to concentrate, and diffuse cultural floodplains where evangelical religions, non-state spiritual practices, and mystery religions can be found.
Collapse
Affiliation(s)
- Victor Møller Poulsen
- Department of Social and Decision Sciences, Carnegie Mellon University, 5000 Forbes Avenue, Pittsburgh, PA 15213, USA
| | - Simon DeDeo
- Department of Social and Decision Sciences, Carnegie Mellon University, 5000 Forbes Avenue, Pittsburgh, PA 15213, USA
- Santa Fe Institute, 1399 Hyde Park Road, Santa Fe, NM 87501, USA
| |
Collapse
|
14
|
Sappington A, Mohanty V. Probabilistic Genotype-Phenotype Maps Reveal Mutational Robustness of RNA Folding, Spin Glasses, and Quantum Circuits. ARXIV 2023:arXiv:2301.01847v1. [PMID: 36713233 PMCID: PMC9882568] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
Recent studies of genotype-phenotype (GP) maps have reported universally enhanced phenotypic robustness to genotype mutations, a feature essential to evolution. Virtually all of these studies make a simplifying assumption that each genotype maps deterministically to a single phenotype. Here, we introduce probabilistic genotype-phenotype (PrGP) maps, where each genotype maps to a vector of phenotype probabilities, as a more realistic framework for investigating robustness. We study three model systems to show that our generalized framework can handle uncertainty emerging from various physical sources: (1) thermal fluctuation in RNA folding, (2) external field disorder in spin glass ground state finding, and (3) superposition and entanglement in quantum circuits, which are realized experimentally on a 7-qubit IBM quantum computer. In all three cases, we observe a novel biphasic robustness scaling which is enhanced relative to random expectation for more frequent phenotypes and approaches random expectation for less frequent phenotypes.
Collapse
Affiliation(s)
- Anna Sappington
- Harvard-MIT Health Sciences and Technology, Harvard Medical School, Boston, MA 02115 and Massachusetts Institute of Technology, Cambridge, MA 02139
| | - Vaibhav Mohanty
- Harvard-MIT Health Sciences and Technology, Harvard Medical School, Boston, MA 02115 and Massachusetts Institute of Technology, Cambridge, MA 02139
| |
Collapse
|
15
|
Delgado MS, López-Galíndez C, Moran F. Viral Fitness Landscapes Based on Self-organizing Maps. Curr Top Microbiol Immunol 2023; 439:95-119. [PMID: 36592243 DOI: 10.1007/978-3-031-15640-3_2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
The creation of fitness maps from viral populations especially in the case of RNA viruses, with high mutation rates producing quasispecies, is complex since the mutant spectrum is in a very high-dimensional space. In this work, a new approach is presented using a class of neural networks, Self-Organized Maps (SOM), to represent realistic fitness landscapes in two RNA viruses: Human Immunodeficiency Virus type 1 (HIV-1) and Hepatitis C Virus (HCV). This methodology has proven to be very effective in the classification of viral quasispecies, using as criterium the mutant sequences in the population. With HIV-1, the fitness landscapes are constructed by representing the experimentally determined fitness on the sequence map. This approach permitted the depiction of the evolutionary paths of the variants subjected to processes of fitness loss and gain in cell culture. In the case of HCV, the efficiency was measured as a function of the frequency of each haplotype in the population by ultra-deep sequencing. The fitness landscapes obtained provided information on the efficiency of each variant in the quasispecies environment, that is, in relation to the entire spectrum of mutants. With the SOM maps, it is possible to determine the evolutionary dynamics of the different haplotypes.
Collapse
Affiliation(s)
- M Soledad Delgado
- Departamento de Sistemas Informáticos, Escuela Técnica Superior de Ingeniería de Sistemas Informáticos (ETSISI), Universidad Politécnica de Madrid, 28031, Madrid, Spain.
| | - Cecilio López-Galíndez
- Unidad de Virología Molecular, Laboratorio de Referencia e Investigación en Retrovirus, Centro Nacional de Microbiología, Instituto de Salud Carlos III, Majadahonda, Madrid, Spain
| | - Federico Moran
- Departamento de Bioquímica y Biología Molecular, Universidad Complutense de Madrid, 28040, Madrid, Spain
| |
Collapse
|
16
|
Mohanty V, Louis AA. Robustness and stability of spin-glass ground states to perturbed interactions. Phys Rev E 2023; 107:014126. [PMID: 36797942 DOI: 10.1103/physreve.107.014126] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2020] [Accepted: 12/16/2022] [Indexed: 06/18/2023]
Abstract
Across many problems in science and engineering, it is important to consider how much the output of a given system changes due to perturbations of the input. Here, we investigate the glassy phase of ±J spin glasses at zero temperature by calculating the robustness of the ground states to flips in the sign of single interactions. For random graphs and the Sherrington-Kirkpatrick model, we find relatively large sets of bond configurations that generate the same ground state. These sets can themselves be analyzed as subgraphs of the interaction domain, and we compute many of their topological properties. In particular, we find that the robustness, equivalent to the average degree, of these subgraphs is much higher than one would expect from a random model. Most notably, it scales in the same logarithmic way with the size of the subgraph as has been found in genotype-phenotype maps for RNA secondary structure folding, protein quaternary structure, gene regulatory networks, as well as for models for genetic programming. The similarity between these disparate systems suggests that this scaling may have a more universal origin.
Collapse
Affiliation(s)
- Vaibhav Mohanty
- Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford, OX1 3NP, United Kingdom
- MD-PhD Program and Program in Health Sciences and Technology, Harvard Medical School, Boston, Massachusetts 02125, USA and Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Ard A Louis
- Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford, OX1 3NP, United Kingdom
| |
Collapse
|
17
|
Abstract
The landscape paradigm is revisited in the light of evolution in simple systems. A brief overview of different classes of fitness landscapes is followed by a more detailed discussion of the RNA model, which is currently the only evolutionary model that allows for a comprehensive molecular analysis of a fitness landscape. Neutral networks of genotypes are indispensable for the success of evolution. Important insights into the evolutionary mechanism are gained by considering the topology of sequence and shape spaces. The dynamic concept of molecular quasispecies is viewed in the light of the landscape paradigm. The distribution of fitness values in state space is mirrored by the population structures of mutant distributions. Two classes of thresholds for replication error or mutations are important: (i) the-conventional-genotypic error threshold, which separates ordered replication from random drift on neutral networks, and (ii) a phenotypic error threshold above which the molecular phenotype is lost. Empirical landscapes are reviewed and finally, the implications of the landscape concept for virus evolution are discussed.
Collapse
Affiliation(s)
- Peter Schuster
- Institut für Theoretische Chemie der Universität Wien, Währingerstraße 17, 1090, Wien, Austria.
| | - Peter F Stadler
- Institut für Informatik der Universität Leipzig, Härtelstraße 16-18, 04107, Leipzig, Germany.,The Santa Fe Institute, 1399 Hyde Park Road, Santa Fe, NM 87501, USA
| |
Collapse
|
18
|
Choudhuri I, Biswas A, Haldane A, Levy RM. Contingency and Entrenchment of Drug-Resistance Mutations in HIV Viral Proteins. J Phys Chem B 2022; 126:10622-10636. [PMID: 36493468 PMCID: PMC9841799 DOI: 10.1021/acs.jpcb.2c06123] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
The ability of HIV-1 to rapidly mutate leads to antiretroviral therapy (ART) failure among infected patients. Drug-resistance mutations (DRMs), which cause a fitness penalty to intrinsic viral fitness, are compensated by accessory mutations with favorable epistatic interactions which cause an evolutionary trapping effect, but the kinetics of this overall process has not been well characterized. Here, using a Potts Hamiltonian model describing epistasis combined with kinetic Monte Carlo simulations of evolutionary trajectories, we explore how epistasis modulates the evolutionary dynamics of HIV DRMs. We show how the occurrence of a drug-resistance mutation is contingent on favorable epistatic interactions with many other residues of the sequence background and that subsequent mutations entrench DRMs. We measure the time-autocorrelation of fluctuations in the likelihood of DRMs due to epistatic coupling with the sequence background, which reveals the presence of two evolutionary processes controlling DRM kinetics with two distinct time scales. Further analysis of waiting times for the evolutionary trapping effect to reverse reveals that the sequences which entrench (trap) a DRM are responsible for the slower time scale. We also quantify the overall strength of epistatic effects on the evolutionary kinetics for different mutations and show these are much larger for DRM positions than polymorphic positions, and we also show that trapping of a DRM is often caused by the collective effect of many accessory mutations, rather than a few strongly coupled ones, suggesting the importance of multiresidue sequence variations in HIV evolution. The analysis presented here provides a framework to explore the kinetic pathways through which viral proteins like HIV evolve under drug-selection pressure.
Collapse
Affiliation(s)
| | | | - Allan Haldane
- Center for Biophysics and Computational Biology, Temple University, Philadelphia, Pennsylvania 19122, United States; Department of Physics, Temple University, Philadelphia, Pennsylvania 19122-6008, United States
| | - Ronald M. Levy
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States; Center for Biophysics and Computational Biology, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
19
|
Zeng HL, Liu Y, Dichio V, Aurell E. Temporal epistasis inference from more than 3 500 000 SARS-CoV-2 genomic sequences. Phys Rev E 2022; 106:044409. [PMID: 36397507 DOI: 10.1103/physreve.106.044409] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2022] [Accepted: 09/19/2022] [Indexed: 06/16/2023]
Abstract
We use direct coupling analysis (DCA) to determine epistatic interactions between loci of variability of the SARS-CoV-2 virus, segmenting genomes by month of sampling. We use full-length, high-quality genomes from the GISAID repository up to October 2021 for a total of over 3 500 000 genomes. We find that DCA terms are more stable over time than correlations but nevertheless change over time as mutations disappear from the global population or reach fixation. Correlations are enriched for phylogenetic effects, and in particularly statistical dependencies at short genomic distances, while DCA brings out links at longer genomic distance. We discuss the validity of a DCA analysis under these conditions in terms of a transient auasilinkage equilibrium state. We identify putative epistatic interaction mutations involving loci in spike.
Collapse
Affiliation(s)
- Hong-Li Zeng
- School of Science, Nanjing University of Posts and Telecommunications, New Energy Technology Engineering Laboratory of Jiangsu Province, Nanjing 210023, China
| | - Yue Liu
- School of Science, Nanjing University of Posts and Telecommunications, New Energy Technology Engineering Laboratory of Jiangsu Province, Nanjing 210023, China
| | - Vito Dichio
- Inria Paris, Aramis Project Team, Paris 75013, France
- Institut du Cerveau, ICM, Inserm U 1127, CNRS UMR 7225, Sorbonne Université, Paris, France
| | - Erik Aurell
- Department of Computational Science and Technology, AlbaNova University Center, SE-106 91 Stockholm, Sweden
| |
Collapse
|
20
|
LaMont C, Otwinowski J, Vanshylla K, Gruell H, Klein F, Nourmohammad A. Design of an optimal combination therapy with broadly neutralizing antibodies to suppress HIV-1. eLife 2022; 11:76004. [PMID: 35852143 PMCID: PMC9467514 DOI: 10.7554/elife.76004] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2021] [Accepted: 07/04/2022] [Indexed: 11/13/2022] Open
Abstract
Infusion of broadly neutralizing antibodies (bNAbs) has shown promise as an alternative to anti-retroviral therapy against HIV. A key challenge is to suppress viral escape, which is more effectively achieved with a combination of bNAbs. Here, we propose a computational approach to predict the efficacy of a bNAb therapy based on the population genetics of HIV escape, which we parametrize using high-throughput HIV sequence data from bNAb-naive patients. By quantifying the mutational target size and the fitness cost of HIV-1 escape from bNAbs, we predict the distribution of rebound times in three clinical trials. We show that a cocktail of three bNAbs is necessary to effectively suppress viral escape, and predict the optimal composition of such bNAb cocktail. Our results offer a rational therapy design for HIV, and show how genetic data can be used to predict treatment outcomes and design new approaches to pathogenic control.
Collapse
Affiliation(s)
- Colin LaMont
- Max Planck Institute for Dynamics and Self-Organization
| | | | | | | | | | | |
Collapse
|
21
|
Multiscale affinity maturation simulations to elicit broadly neutralizing antibodies against HIV. PLoS Comput Biol 2022; 18:e1009391. [PMID: 35442968 PMCID: PMC9020693 DOI: 10.1371/journal.pcbi.1009391] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2021] [Accepted: 03/24/2022] [Indexed: 11/29/2022] Open
Abstract
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. Vaccination has saved more lives than any other medical procedure. 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 the in silico investigation 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.
Collapse
|
22
|
Biswas A, Haldane A, Levy RM. Limits to detecting epistasis in the fitness landscape of HIV. PLoS One 2022; 17:e0262314. [PMID: 35041711 PMCID: PMC8765623 DOI: 10.1371/journal.pone.0262314] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2021] [Accepted: 12/20/2021] [Indexed: 02/05/2023] Open
Abstract
The rapid evolution of HIV is constrained by interactions between mutations which affect viral fitness. In this work, we explore the role of epistasis in determining the mutational fitness landscape of HIV for multiple drug target proteins, including Protease, Reverse Transcriptase, and Integrase. Epistatic interactions between residues modulate the mutation patterns involved in drug resistance, with unambiguous signatures of epistasis best seen in the comparison of the Potts model predicted and experimental HIV sequence “prevalences” expressed as higher-order marginals (beyond triplets) of the sequence probability distribution. In contrast, experimental measures of fitness such as viral replicative capacities generally probe fitness effects of point mutations in a single background, providing weak evidence for epistasis in viral systems. The detectable effects of epistasis are obscured by higher evolutionary conservation at sites. While double mutant cycles in principle, provide one of the best ways to probe epistatic interactions experimentally without reference to a particular background, we show that the analysis is complicated by the small dynamic range of measurements. Overall, we show that global pairwise interaction Potts models are necessary for predicting the mutational landscape of viral proteins.
Collapse
Affiliation(s)
- Avik Biswas
- Department of Physics, Temple University, Philadelphia, PA, United States of America
- Center for Biophysics and Computational Biology, Temple University, Philadelphia, PA, United States of America
| | - Allan Haldane
- Department of Physics, Temple University, Philadelphia, PA, United States of America
- Center for Biophysics and Computational Biology, Temple University, Philadelphia, PA, United States of America
| | - Ronald M. Levy
- Department of Physics, Temple University, Philadelphia, PA, United States of America
- Center for Biophysics and Computational Biology, Temple University, Philadelphia, PA, United States of America
- Department of Chemistry, Temple University, Philadelphia, PA, United States of America
- * E-mail:
| |
Collapse
|
23
|
Doelger J, Kardar M, Chakraborty AK. Inferring the intrinsic mutational fitness landscape of influenzalike evolving antigens from temporally ordered sequence data. Phys Rev E 2022; 105:024401. [PMID: 35291059 DOI: 10.1103/physreve.105.024401] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2021] [Accepted: 01/19/2022] [Indexed: 06/14/2023]
Abstract
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. To find 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 influenzalike 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 the design of immunization protocols by identifying intrinsically vulnerable immune target combinations on antigens that evolve under immune-driven selection. In the future, this approach may 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.
Collapse
Affiliation(s)
- Julia Doelger
- Institute for Medical Engineering and Science, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Mehran Kardar
- Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Arup K Chakraborty
- Institute for Medical Engineering and Science, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA; Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA; Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA; Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA; and Ragon Institute of MGH, MIT, and Harvard, Cambridge, Massachusetts 02139, USA
| |
Collapse
|
24
|
Epistatic models predict mutable sites in SARS-CoV-2 proteins and epitopes. Proc Natl Acad Sci U S A 2022; 119:2113118119. [PMID: 35022216 PMCID: PMC8795541 DOI: 10.1073/pnas.2113118119] [Citation(s) in RCA: 44] [Impact Index Per Article: 22.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/13/2021] [Indexed: 12/21/2022] Open
Abstract
During the COVID pandemic, new severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants emerge and spread, some being of major concern due to their increased infectivity or capacity to reduce vaccine efficiency. Anticipating mutations, which might give rise to new variants, would be of great interest. We construct sequence models predicting how mutable SARS-CoV-2 positions are, using a single SARS-CoV-2 sequence and databases of other coronaviruses. Predictions are tested against available mutagenesis data and the observed variability of SARS-CoV-2 proteins. Interestingly, predictions agree increasingly with observations, as more SARS-CoV-2 sequences become available. Combining predictions with immunological data, we find an overrepresentation of mutations in current variants of concern. The approach may become relevant for potential outbreaks of future viral diseases. The emergence of new variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a major concern given their potential impact on the transmissibility and pathogenicity of the virus as well as the efficacy of therapeutic interventions. Here, we predict the mutability of all positions in SARS-CoV-2 protein domains to forecast the appearance of unseen variants. Using sequence data from other coronaviruses, preexisting to SARS-CoV-2, we build statistical models that not only capture amino acid conservation but also more complex patterns resulting from epistasis. We show that these models are notably superior to conservation profiles in estimating the already observable SARS-CoV-2 variability. In the receptor binding domain of the spike protein, we observe that the predicted mutability correlates well with experimental measures of protein stability and that both are reliable mutability predictors (receiver operating characteristic areas under the curve ∼0.8). Most interestingly, we observe an increasing agreement between our model and the observed variability as more data become available over time, proving the anticipatory capacity of our model. When combined with data concerning the immune response, our approach identifies positions where current variants of concern are highly overrepresented. These results could assist studies on viral evolution and future viral outbreaks and, in particular, guide the exploration and anticipation of potentially harmful future SARS-CoV-2 variants.
Collapse
|
25
|
Evolutionary modeling reveals enhanced mutational flexibility of HCV subtype 1b compared with 1a. iScience 2022; 25:103569. [PMID: 34988406 PMCID: PMC8704487 DOI: 10.1016/j.isci.2021.103569] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2021] [Revised: 11/19/2021] [Accepted: 12/02/2021] [Indexed: 11/24/2022] Open
Abstract
Hepatitis C virus (HCV) is a leading cause of liver-associated disease and liver cancer. Of the major HCV subtypes, patients infected with subtype 1b have been associated with having a higher risk of developing chronic infection and hepatocellular carcinoma. However, underlying reasons for this increased disease severity remain unknown. Here, we provide an evolutionary rationale, based on a comparative study of fitness landscape and in-host evolutionary models of the E2 glycoprotein of HCV subtypes 1a and 1b. Our analysis demonstrates that a higher chronicity rate of 1b may be attributed to lower fitness constraints, enabling 1b viruses to more easily escape antibody responses. More generally, our results suggest that differences in evolutionary constraints between HCV subtypes may be an important factor in mediating distinct disease outcomes. Our analysis also identifies antibodies that appear escape-resistant against both subtypes 1a and 1b, providing directions for designing HCV vaccines having cross-subtype protection. Comparative analysis of the fitness landscapes of HCV subtypes 1a and 1b Subtype 1b evolution is subject to less constraints than 1a Subtype 1b appears to evade antibodies more easily compared with 1a Antibodies are identified that are difficult to escape for both subtypes 1a and 1b
Collapse
|
26
|
Abstract
Antibody immunodominance refers to the preferential and asymmetric elicitation of antibodies against specific epitopes on a complex protein antigen. Traditional vaccination approaches for rapidly evolving pathogens have had limited success in part because of this phenomenon, as elicited antibodies preferentially target highly variable regions of antigens, and thus do not confer long lasting protection. While antibodies targeting functionally conserved epitopes have the potential to be broadly protective, they often make up a minority of the overall repertoire. Here, we discuss recent protein engineering strategies used to favorably alter patterns of immunodominance, and selectively focus antibody responses toward broadly protective epitopes in the pursuit of next-generation vaccines for rapidly evolving pathogens.
Collapse
|
27
|
Viral surface geometry shapes influenza and coronavirus spike evolution through antibody pressure. PLoS Comput Biol 2021; 17:e1009664. [PMID: 34898597 PMCID: PMC8699686 DOI: 10.1371/journal.pcbi.1009664] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2020] [Revised: 12/23/2021] [Accepted: 11/19/2021] [Indexed: 01/02/2023] Open
Abstract
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. The immune system responds to viruses by making neutralizing antibodies to regions of the viral spike protein, which mutates to escape. To inform vaccine design and understand how the fitness landscape of the viral spike changes over time, it is necessary to identify and quantify the factors directing its evolution. Based on the 3D structure of the viral surface and spike as captured with Cryo-EM and crystallography, we aimed to create a coarse-grained model for the effect of antibodies in forcing surface residues of the spike to mutate. We found that for pre-pandemic influenza (hemagglutinin) and the corona sarbecovirus subgenus (S protein), the location of a residue on the spike protein, which modulates its accessibility to antibodies, highly correlates with its propensity to mutate. Hence, a mechanistic approach can be used to identify aspects of viral spike sequence diversity related to antibody escape.
Collapse
|
28
|
Laine E, Eismann S, Elofsson A, Grudinin S. Protein sequence-to-structure learning: Is this the end(-to-end revolution)? Proteins 2021; 89:1770-1786. [PMID: 34519095 DOI: 10.1002/prot.26235] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2021] [Revised: 08/16/2021] [Accepted: 09/03/2021] [Indexed: 01/08/2023]
Abstract
The potential of deep learning has been recognized in the protein structure prediction community for some time, and became indisputable after CASP13. In CASP14, deep learning has boosted the field to unanticipated levels reaching near-experimental accuracy. This success comes from advances transferred from other machine learning areas, as well as methods specifically designed to deal with protein sequences and structures, and their abstractions. Novel emerging approaches include (i) geometric learning, that is, learning on representations such as graphs, three-dimensional (3D) Voronoi tessellations, and point clouds; (ii) pretrained protein language models leveraging attention; (iii) equivariant architectures preserving the symmetry of 3D space; (iv) use of large meta-genome databases; (v) combinations of protein representations; and (vi) finally truly end-to-end architectures, that is, differentiable models starting from a sequence and returning a 3D structure. Here, we provide an overview and our opinion of the novel deep learning approaches developed in the last 2 years and widely used in CASP14.
Collapse
Affiliation(s)
- Elodie Laine
- Sorbonne Université, CNRS, IBPS, Laboratoire de Biologie Computationnelle et Quantitative (LCQB), Paris, France
| | - Stephan Eismann
- Department of Computer Science and Applied Physics, Stanford University, Stanford, California, USA
| | - Arne Elofsson
- Department of Biochemistry and Biophysics and Science for Life Laboratory, Stockholm University, Solna, Sweden
| | - Sergei Grudinin
- Univ. Grenoble Alpes, CNRS, Grenoble INP, LJK, Grenoble, France
| |
Collapse
|
29
|
Abstract
Broadly neutralizing antibodies (bnAbs) show promise for antibody-mediated prevention or treatment of HIV-1 infection. Recent clinical trials, however, indicate that the virus can evolve resistance mutations that escape neutralization by bnAbs. Here, we establish a fitness model for the escape dynamics of HIV in humans. We show that this model can be applied universally across different human hosts, providing a proof of principle that the in vivo response of HIV to bnAb therapies is predictable. Our analysis identifies a fitness trade-off in viral evolution that can inform antibody design and therapy protocols. Broadly neutralizing antibodies are promising candidates for treatment and prevention of HIV-1 infections. Such antibodies can temporarily suppress viral load in infected individuals; however, the virus often rebounds by escape mutants that have evolved resistance. In this paper, we map a fitness model of HIV-1 interacting with broadly neutralizing antibodies using in vivo data from a recent clinical trial. We identify two fitness factors, antibody dosage and viral load, that determine viral reproduction rates reproducibly across different hosts. The model successfully predicts the escape dynamics of HIV-1 in the course of an antibody treatment, including a characteristic frequency turnover between sensitive and resistant strains. This turnover is governed by a dosage-dependent fitness ranking, resulting from an evolutionary trade-off between antibody resistance and its collateral cost in drug-free growth. Our analysis suggests resistance–cost trade-off curves as a measure of antibody performance in the presence of resistance evolution.
Collapse
|
30
|
Zhang Y, Chowdhury S, Rodrigues JV, Shakhnovich E. Development of antibacterial compounds that constrain evolutionary pathways to resistance. eLife 2021; 10:64518. [PMID: 34279221 PMCID: PMC8331180 DOI: 10.7554/elife.64518] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2020] [Accepted: 07/13/2021] [Indexed: 01/27/2023] Open
Abstract
Antibiotic resistance is a worldwide challenge. A potential approach to block resistance is to simultaneously inhibit WT and known escape variants of the target bacterial protein. Here, we applied an integrated computational and experimental approach to discover compounds that inhibit both WT and trimethoprim (TMP) resistant mutants of E. coli dihydrofolate reductase (DHFR). We identified a novel compound (CD15-3) that inhibits WT DHFR and its TMP resistant variants L28R, P21L and A26T with IC50 50–75 µM against WT and TMP-resistant strains. Resistance to CD15-3 was dramatically delayed compared to TMP in in vitro evolution. Whole genome sequencing of CD15-3-resistant strains showed no mutations in the target folA locus. Rather, gene duplication of several efflux pumps gave rise to weak (about twofold increase in IC50) resistance against CD15-3. Altogether, our results demonstrate the promise of strategy to develop evolution drugs - compounds which constrain evolutionary escape routes in pathogens.
Collapse
Affiliation(s)
- Yanmin Zhang
- School of Science, China Pharmaceutical University, Nanjing, China.,Department of Chemistry and Chemical Biology, Harvard University, Cambridge, United States
| | - Sourav Chowdhury
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, United States
| | - João V Rodrigues
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, United States
| | - Eugene Shakhnovich
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, United States
| |
Collapse
|
31
|
Sheng J, Wang S. Coevolutionary transitions emerging from flexible molecular recognition and eco-evolutionary feedback. iScience 2021; 24:102861. [PMID: 34401660 PMCID: PMC8353512 DOI: 10.1016/j.isci.2021.102861] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2020] [Revised: 05/16/2021] [Accepted: 07/13/2021] [Indexed: 01/16/2023] Open
Abstract
Highly mutable viruses evolve to evade host immunity that exerts selective pressure and adapts to viral dynamics. Here, we provide a framework for identifying key determinants of the mode and fate of viral-immune coevolution by linking molecular recognition and eco-evolutionary dynamics. We find that conservation level and initial diversity of antigen jointly determine the timing and efficacy of narrow and broad antibody responses, which in turn control the transition between viral persistence, clearance, and rebound. In particular, clearance of structurally complex antigens relies on antibody evolution in a larger antigenic space than where selection directly acts; viral rebound manifests binding-mediated feedback between ecology and rapid evolution. Finally, immune compartmentalization can slow viral escape but also delay clearance. This work suggests that flexible molecular binding allows a plastic phenotype that exploits potentiating neutral variations outside direct contact, opening new and shorter paths toward highly adaptable states. A scale-crossing framework identifies key determinants of viral-immune coevolution Fast specific response influences slow broad response by shaping antigen dynamics Antibody footprint shift enables breadth acquisition and viral clearance Model explains divergent kinetics and outcomes of HCV infection in humans
Collapse
Affiliation(s)
- Jiming Sheng
- Department of Physics and Astronomy, University of California, Los Angeles, Los Angeles, CA 90095, USA
| | - Shenshen Wang
- Department of Physics and Astronomy, University of California, Los Angeles, Los Angeles, CA 90095, USA
| |
Collapse
|
32
|
Pedruzzi G, Rouzine IM. An evolution-based high-fidelity method of epistasis measurement: Theory and application to influenza. PLoS Pathog 2021; 17:e1009669. [PMID: 34153082 PMCID: PMC8248644 DOI: 10.1371/journal.ppat.1009669] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2020] [Revised: 07/01/2021] [Accepted: 05/25/2021] [Indexed: 12/18/2022] Open
Abstract
Linkage effects in a multi-locus population strongly influence its evolution. The models based on the traveling wave approach enable us to predict the average speed of evolution and the statistics of phylogeny. However, predicting statistically the evolution of specific sites and pairs of sites in the multi-locus context remains a mathematical challenge. In particular, the effects of epistasis, the interaction of gene regions contributing to phenotype, is difficult to predict theoretically and detect experimentally in sequence data. A large number of false-positive interactions arises from stochastic linkage effects and indirect interactions, which mask true epistatic interactions. Here we develop a proof-of-principle method to filter out false-positive interactions. We start by demonstrating that the averaging of haplotype frequencies over multiple independent populations is necessary but not sufficient for epistatic detection, because it still leaves high numbers of false-positive interactions. To compensate for the residual stochastic noise, we develop a three-way haplotype method isolating true interactions. The fidelity of the method is confirmed analytically and on simulated genetic sequences evolved with a known epistatic network. The method is then applied to a large sequence database of neurominidase protein of influenza A H1N1 obtained from various geographic locations to infer the epistatic network responsible for the difference between the pre-pandemic virus and the pandemic strain of 2009. These results present a simple and reliable technique to measure epistatic interactions of any sign from sequence data. Interactions between genomic sites create a fitness landscape. The knowledge of topology and strength of interactions is vital for predicting the escape of viruses from drugs and immune response and their passing through fitness valleys. Many efforts have been invested into measuring these interactions from DNA sequence sets. Unfortunately, reproducibility of the results remains low due partly to a very small fraction of interaction pairs and partly to stochastic linkage noise masking true interactions. Here we propose a method to separate stochastic linkage and indirect interactions from epistatic interactions and apply it to influenza virus sequence data.
Collapse
Affiliation(s)
- Gabriele Pedruzzi
- Sorbonne Université, Institute de Biologie Paris-Seine, Laboratoire de Biologie Computationelle et Quantitative LCQB, Paris, France
| | - Igor M. Rouzine
- Sorbonne Université, Institute de Biologie Paris-Seine, Laboratoire de Biologie Computationelle et Quantitative LCQB, Paris, France
- * E-mail:
| |
Collapse
|
33
|
Fernandez-de-Cossio-Diaz J, Uguzzoni G, Pagnani A. Unsupervised Inference of Protein Fitness Landscape from Deep Mutational Scan. Mol Biol Evol 2021; 38:318-328. [PMID: 32770229 PMCID: PMC7783173 DOI: 10.1093/molbev/msaa204] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
The recent technological advances underlying the screening of large combinatorial libraries in high-throughput mutational scans deepen our understanding of adaptive protein evolution and boost its applications in protein design. Nevertheless, the large number of possible genotypes requires suitable computational methods for data analysis, the prediction of mutational effects, and the generation of optimized sequences. We describe a computational method that, trained on sequencing samples from multiple rounds of a screening experiment, provides a model of the genotype-fitness relationship. We tested the method on five large-scale mutational scans, yielding accurate predictions of the mutational effects on fitness. The inferred fitness landscape is robust to experimental and sampling noise and exhibits high generalization power in terms of broader sequence space exploration and higher fitness variant predictions. We investigate the role of epistasis and show that the inferred model provides structural information about the 3D contacts in the molecular fold.
Collapse
Affiliation(s)
- Jorge Fernandez-de-Cossio-Diaz
- Systems Biology Department, Center of Molecular Immunology, Havana, Cuba.,Laboratory of Physics of the Ecole Normale Superieure, CNRS UMR 8023 & PSL Research, Paris, France
| | | | - Andrea Pagnani
- Politecnico di Torino, Torino, Italy.,Italian Institute for Genomic Medicine, IRCCS Candiolo, Candiolo, TO, Italy.,INFN, Sezione di Torino, Torino, Italy
| |
Collapse
|
34
|
Adenovirus-vectored vaccine containing multidimensionally conserved parts of the HIV proteome is immunogenic in rhesus macaques. Proc Natl Acad Sci U S A 2021; 118:2022496118. [PMID: 33514660 DOI: 10.1073/pnas.2022496118] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
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.
Collapse
|
35
|
McCann CD, van Dorp CH, Danesh A, Ward AR, Dilling TR, Mota TM, Zale E, Stevenson EM, Patel S, Brumme CJ, Dong W, Jones DS, Andresen TL, Walker BD, Brumme ZL, Bollard CM, Perelson AS, Irvine DJ, Jones RB. A participant-derived xenograft model of HIV enables long-term evaluation of autologous immunotherapies. J Exp Med 2021; 218:212105. [PMID: 33988715 PMCID: PMC8129803 DOI: 10.1084/jem.20201908] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2020] [Revised: 02/15/2021] [Accepted: 04/12/2021] [Indexed: 12/13/2022] Open
Abstract
HIV-specific CD8+ T cells partially control viral replication and delay disease progression, but they rarely provide lasting protection, largely due to immune escape. Here, we show that engrafting mice with memory CD4+ T cells from HIV+ donors uniquely allows for the in vivo evaluation of autologous T cell responses while avoiding graft-versus-host disease and the need for human fetal tissues that limit other models. Treating HIV-infected mice with clinically relevant HIV-specific T cell products resulted in substantial reductions in viremia. In vivo activity was significantly enhanced when T cells were engineered with surface-conjugated nanogels carrying an IL-15 superagonist, but it was ultimately limited by the pervasive selection of a diverse array of escape mutations, recapitulating patterns seen in humans. By applying mathematical modeling, we show that the kinetics of the CD8+ T cell response have a profound impact on the emergence and persistence of escape mutations. This “participant-derived xenograft” model of HIV provides a powerful tool for studying HIV-specific immunological responses and facilitating the development of effective cell-based therapies.
Collapse
Affiliation(s)
- Chase D McCann
- Infectious Diseases Division, Department of Medicine, Weill Cornell Medical College, New York, NY.,Immunology & Microbial Pathogenesis Graduate Program, Weill Cornell Graduate School of Medical Sciences, New York, NY
| | | | - Ali Danesh
- Infectious Diseases Division, Department of Medicine, Weill Cornell Medical College, New York, NY
| | - Adam R Ward
- Department of Microbiology, Immunology & Tropical Medicine, The George Washington University School of Medicine and Health Sciences, Washington, DC.,PhD Program in Epidemiology, Milken Institute School of Public Health, The George Washington University, Washington, DC
| | - Thomas R Dilling
- Infectious Diseases Division, Department of Medicine, Weill Cornell Medical College, New York, NY
| | - Talia M Mota
- Infectious Diseases Division, Department of Medicine, Weill Cornell Medical College, New York, NY
| | - Elizabeth Zale
- Immunology & Microbial Pathogenesis Graduate Program, Weill Cornell Graduate School of Medical Sciences, New York, NY
| | - Eva M Stevenson
- Infectious Diseases Division, Department of Medicine, Weill Cornell Medical College, New York, NY
| | - Shabnum Patel
- Center for Cancer and Immunology Research, Children's National Health System, Washington, DC.,George Washington University Cancer Center, George Washington University, Washington, DC
| | - Chanson J Brumme
- British Columbia Centre for Excellence in HIV/AIDS, Vancouver, British Columbia, Canada
| | - Winnie Dong
- British Columbia Centre for Excellence in HIV/AIDS, Vancouver, British Columbia, Canada
| | | | | | - Bruce D Walker
- Ragon Institute of Massachusetts General Hospital, Massachusetts Institute of Technology, and Harvard University, Boston, MA.,Institute for Medical and Engineering Sciences, Massachusetts Institute of Technology, Cambridge, MA.,Howard Hughes Medical Institute, Chevy Chase, MD
| | - Zabrina L Brumme
- British Columbia Centre for Excellence in HIV/AIDS, Vancouver, British Columbia, Canada.,Faculty of Health Sciences, Simon Fraser University, Burnaby, British Columbia, Canada
| | - Catherine M Bollard
- Center for Cancer and Immunology Research, Children's National Health System, Washington, DC.,George Washington University Cancer Center, George Washington University, Washington, DC
| | - Alan S Perelson
- Theoretical Biology and Biophysics, Los Alamos National Laboratory, Los Alamos, NM
| | - Darrell J Irvine
- Howard Hughes Medical Institute, Chevy Chase, MD.,Department of Material Science and Engineering, Massachusetts Institute of Technology, Cambridge, MA.,Koch Institute for Integrative Cancer Research, Massachusetts Institute of Technology, Cambridge, MA.,Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, MA
| | - R Brad Jones
- Infectious Diseases Division, Department of Medicine, Weill Cornell Medical College, New York, NY.,Immunology & Microbial Pathogenesis Graduate Program, Weill Cornell Graduate School of Medical Sciences, New York, NY
| |
Collapse
|
36
|
Garrett ME, Galloway J, Chu HY, Itell HL, Stoddard CI, Wolf CR, Logue JK, McDonald D, Weight H, Matsen FA, Overbaugh J. High-resolution profiling of pathways of escape for SARS-CoV-2 spike-binding antibodies. Cell 2021; 184:2927-2938.e11. [PMID: 34010620 PMCID: PMC8096189 DOI: 10.1016/j.cell.2021.04.045] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2020] [Revised: 02/26/2021] [Accepted: 04/23/2021] [Indexed: 01/02/2023]
Abstract
Defining long-term protective immunity to SARS-CoV-2 is one of the most pressing questions of our time and will require a detailed understanding of potential ways this virus can evolve to escape immune protection. Immune protection will most likely be mediated by antibodies that bind to the viral entry protein, spike (S). Here, we used Phage-DMS, an approach that comprehensively interrogates the effect of all possible mutations on binding to a protein of interest, to define the profile of antibody escape to the SARS-CoV-2 S protein using coronavirus disease 2019 (COVID-19) convalescent plasma. Antibody binding was common in two regions, the fusion peptide and the linker region upstream of the heptad repeat region 2. However, escape mutations were variable within these immunodominant regions. There was also individual variation in less commonly targeted epitopes. This study provides a granular view of potential antibody escape pathways and suggests there will be individual variation in antibody-mediated virus evolution.
Collapse
Affiliation(s)
- Meghan E Garrett
- Division of Human Biology, Fred Hutchinson Cancer Research Center, Seattle, WA 98102, USA; Molecular and Cellular Biology Graduate Program, University of Washington and Fred Hutchinson Cancer Research Center, Seattle, WA 98195, USA
| | - Jared Galloway
- Division of Public Health Sciences, Fred Hutchinson Cancer Research Center, Seattle, WA 98102, USA
| | - Helen Y Chu
- Department of Medicine, University of Washington, Seattle, WA 98195, USA
| | - Hannah L Itell
- Division of Human Biology, Fred Hutchinson Cancer Research Center, Seattle, WA 98102, USA; Molecular and Cellular Biology Graduate Program, University of Washington and Fred Hutchinson Cancer Research Center, Seattle, WA 98195, USA
| | - Caitlin I Stoddard
- Division of Human Biology, Fred Hutchinson Cancer Research Center, Seattle, WA 98102, USA
| | - Caitlin R Wolf
- Department of Medicine, University of Washington, Seattle, WA 98195, USA
| | - Jennifer K Logue
- Department of Medicine, University of Washington, Seattle, WA 98195, USA
| | - Dylan McDonald
- Department of Medicine, University of Washington, Seattle, WA 98195, USA
| | - Haidyn Weight
- Division of Human Biology, Fred Hutchinson Cancer Research Center, Seattle, WA 98102, USA
| | - Frederick A Matsen
- Division of Public Health Sciences, Fred Hutchinson Cancer Research Center, Seattle, WA 98102, USA; Computational Biology Program, Fred Hutchinson Cancer Research Center, Seattle, WA 98102, USA
| | - Julie Overbaugh
- Division of Human Biology, Fred Hutchinson Cancer Research Center, Seattle, WA 98102, USA; Division of Public Health Sciences, Fred Hutchinson Cancer Research Center, Seattle, WA 98102, USA.
| |
Collapse
|
37
|
Gao A, Chen Z, Amitai A, Doelger J, Mallajosyula V, Sundquist E, Pereyra Segal F, Carrington M, Davis MM, Streeck H, Chakraborty AK, Julg B. Learning from HIV-1 to predict the immunogenicity of T cell epitopes in SARS-CoV-2. iScience 2021; 24:102311. [PMID: 33748696 PMCID: PMC7956900 DOI: 10.1016/j.isci.2021.102311] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2021] [Revised: 02/22/2021] [Accepted: 03/10/2021] [Indexed: 12/18/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Ang Gao
- Institute for Medical Engineering & Science, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Zhilin Chen
- Ragon Institute of Massachusetts General Hospital, Massachusetts Institute of Technology, and Harvard, 400 Technology Sq., Cambridge, MA 02139, USA
| | - Assaf Amitai
- Institute for Medical Engineering & Science, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA
| | - Julia Doelger
- Institute for Medical Engineering & Science, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA
| | - Vamsee Mallajosyula
- Institute for Immunity, Transplantation, and Infection, Stanford University School of Medicine, Stanford, CA 94305, USA
| | - Emily Sundquist
- Ragon Institute of Massachusetts General Hospital, Massachusetts Institute of Technology, and Harvard, 400 Technology Sq., Cambridge, MA 02139, USA
| | | | - Mary Carrington
- Ragon Institute of Massachusetts General Hospital, Massachusetts Institute of Technology, and Harvard, 400 Technology Sq., Cambridge, MA 02139, USA
- Basic Science Program, Frederick National Laboratory for Cancer Research, Frederick, MD 21702, USA
| | - Mark M. Davis
- Institute for Immunity, Transplantation, and Infection, Stanford University School of Medicine, Stanford, CA 94305, USA
- Department of Microbiology and Immunology, Stanford University School of Medicine, Stanford, CA 94305, USA
- Howard Hughes Medical Institute, Stanford University School of Medicine, Stanford, CA 94305, USA
| | - Hendrik Streeck
- Institut für Virologie, Universitätsklinikum Bonn, 53127 Bonn, Germany
| | - Arup K. Chakraborty
- Institute for Medical Engineering & Science, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
- Ragon Institute of Massachusetts General Hospital, Massachusetts Institute of Technology, and Harvard, 400 Technology Sq., Cambridge, MA 02139, USA
- Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Boris Julg
- Ragon Institute of Massachusetts General Hospital, Massachusetts Institute of Technology, and Harvard, 400 Technology Sq., Cambridge, MA 02139, USA
| |
Collapse
|
38
|
Narayanan KK, Procko E. Deep Mutational Scanning of Viral Glycoproteins and Their Host Receptors. Front Mol Biosci 2021; 8:636660. [PMID: 33898517 PMCID: PMC8062978 DOI: 10.3389/fmolb.2021.636660] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2020] [Accepted: 03/18/2021] [Indexed: 11/17/2022] Open
Abstract
Deep mutational scanning or deep mutagenesis is a powerful tool for understanding the sequence diversity available to viruses for adaptation in a laboratory setting. It generally involves tracking an in vitro selection of protein sequence variants with deep sequencing to map mutational effects based on changes in sequence abundance. Coupled with any of a number of selection strategies, deep mutagenesis can explore the mutational diversity available to viral glycoproteins, which mediate critical roles in cell entry and are exposed to the humoral arm of the host immune response. Mutational landscapes of viral glycoproteins for host cell attachment and membrane fusion reveal extensive epistasis and potential escape mutations to neutralizing antibodies or other therapeutics, as well as aiding in the design of optimized immunogens for eliciting broadly protective immunity. While less explored, deep mutational scans of host receptors further assist in understanding virus-host protein interactions. Critical residues on the host receptors for engaging with viral spikes are readily identified and may help with structural modeling. Furthermore, mutations may be found for engineering soluble decoy receptors as neutralizing agents that specifically bind viral targets with tight affinity and limited potential for viral escape. By untangling the complexities of how sequence contributes to viral glycoprotein and host receptor interactions, deep mutational scanning is impacting ideas and strategies at multiple levels for combatting circulating and emergent virus strains.
Collapse
Affiliation(s)
| | - Erik Procko
- Department of Biochemistry and Cancer Center at Illinois, University of Illinois, Urbana, IL, United States
| |
Collapse
|
39
|
Sohail MS, Ahmed SF, Quadeer AA, McKay MR. In silico T cell epitope identification for SARS-CoV-2: Progress and perspectives. Adv Drug Deliv Rev 2021; 171:29-47. [PMID: 33465451 PMCID: PMC7832442 DOI: 10.1016/j.addr.2021.01.007] [Citation(s) in RCA: 43] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2020] [Revised: 12/31/2020] [Accepted: 01/07/2021] [Indexed: 02/06/2023]
Abstract
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.
Collapse
Affiliation(s)
- Muhammad Saqib Sohail
- Department of Electronic and Computer Engineering, The Hong Kong University of Science and Technology, Hong Kong, China
| | - Syed Faraz Ahmed
- Department of Electronic and Computer Engineering, The Hong Kong University of Science and Technology, Hong Kong, China
| | - Ahmed Abdul Quadeer
- Department of Electronic and Computer Engineering, The Hong Kong University of Science and Technology, Hong Kong, China.
| | - Matthew R McKay
- Department of Electronic and Computer Engineering, The Hong Kong University of Science and Technology, Hong Kong, China; Department of Chemical and Biological Engineering, The Hong Kong University of Science and Technology, Hong Kong, China.
| |
Collapse
|
40
|
Zhen A, Carrillo MA, Mu W, Rezek V, Martin H, Hamid P, Chen ISY, Yang OO, Zack JA, Kitchen SG. Robust CAR-T memory formation and function via hematopoietic stem cell delivery. PLoS Pathog 2021; 17:e1009404. [PMID: 33793675 PMCID: PMC8016106 DOI: 10.1371/journal.ppat.1009404] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2020] [Accepted: 02/17/2021] [Indexed: 12/16/2022] Open
Abstract
Due to the durability and persistence of reservoirs of HIV-1-infected cells, combination antiretroviral therapy (ART) is insufficient in eradicating infection. Achieving HIV-1 cure or sustained remission without ART treatment will require the enhanced and persistent effective antiviral immune responses. Chimeric Antigen Receptor (CAR) T-cells have emerged as a powerful immunotherapy and show promise in treating HIV-1 infection. Persistence, trafficking, and maintenance of function remain to be a challenge in many of these approaches, which are based on peripheral T cell modification. To overcome many of these issues, we have previously demonstrated successful long-term engraftment and production of anti-HIV CAR T cells in modified hematopoietic stem cells (HSCs) in vivo. Here we report the development and in vivo testing of second generation CD4-based CARs (CD4CAR) against HIV-1 infection using a HSCs-based approach. We found that a modified, truncated CD4-based CAR (D1D2CAR) allows better CAR-T cell differentiation from gene modified HSCs, and maintains similar CTL activity as compared to the full length CD4-based CAR. In addition, D1D2CAR does not mediate HIV infection or stimulation mediated by IL-16, suggesting lower risk of off-target effects. Interestingly, stimulatory domains of 4-1BB but not CD28 allowed successful hematopoietic differentiation and improved anti-viral function of CAR T cells from CAR modified HSCs. Addition of 4-1BB to CD4 based CARs led to faster suppression of viremia during early untreated HIV-1 infection. D1D2CAR 4-1BB mice had faster viral suppression in combination with ART and better persistence of CAR T cells during ART. In summary, our data indicate that the D1D2CAR-41BB is a superior CAR, showing better HSC differentiation, viral suppression and persistence, and less deleterious functions compared to the original CD4CAR, and should continue to be pursued as a candidate for clinical study. Engineering T cells with anti-HIV chimeric antigen receptors (CAR) has emerged as a promising strategy to control HIV infection through a genetic vaccination strategy. Here we report a novel CAR-based approach targeting HIV infection using the genetic modification of blood forming hematopoietic stem cells (HSCs). This novel CAR approach uses a modified HIV receptor molecule (the primary HIV receptor CD4) as well as anti-HIV agents to modify HSCs to allow them to develop into cells that are protected from HIV infection and target HIV infected cells for the life of the individual. We found this latest generation of CARs successfully modified and allowed in vivo engraftment that resulted in the development of effective anti-HIV CAR T cells with robust memory formation and viral control. Our study highlights the identification of a next-generation CAR molecule that protected cells from infection, targeted and reduced HIV burdens, and serves as an ideal developmental candidate for further clinical studies.
Collapse
Affiliation(s)
- Anjie Zhen
- Division of Hematology/Oncology, Department of Medicine, David Geffen School of Medicine at UCLA, Los Angeles, California, United States of America
- UCLA AIDS Institute and the Eli and Edythe Broad Center of Regenerative Medicine and Stem Cell Research, David Geffen School of Medicine at UCLA, Los Angeles, California, United States of America
| | - Mayra A. Carrillo
- Division of Hematology/Oncology, Department of Medicine, David Geffen School of Medicine at UCLA, Los Angeles, California, United States of America
- UCLA AIDS Institute and the Eli and Edythe Broad Center of Regenerative Medicine and Stem Cell Research, David Geffen School of Medicine at UCLA, Los Angeles, California, United States of America
| | - Wenli Mu
- Division of Hematology/Oncology, Department of Medicine, David Geffen School of Medicine at UCLA, Los Angeles, California, United States of America
- UCLA AIDS Institute and the Eli and Edythe Broad Center of Regenerative Medicine and Stem Cell Research, David Geffen School of Medicine at UCLA, Los Angeles, California, United States of America
| | - Valerie Rezek
- Division of Hematology/Oncology, Department of Medicine, David Geffen School of Medicine at UCLA, Los Angeles, California, United States of America
- UCLA AIDS Institute and the Eli and Edythe Broad Center of Regenerative Medicine and Stem Cell Research, David Geffen School of Medicine at UCLA, Los Angeles, California, United States of America
| | - Heather Martin
- Division of Hematology/Oncology, Department of Medicine, David Geffen School of Medicine at UCLA, Los Angeles, California, United States of America
- UCLA AIDS Institute and the Eli and Edythe Broad Center of Regenerative Medicine and Stem Cell Research, David Geffen School of Medicine at UCLA, Los Angeles, California, United States of America
| | - Philip Hamid
- Division of Hematology/Oncology, Department of Medicine, David Geffen School of Medicine at UCLA, Los Angeles, California, United States of America
- UCLA AIDS Institute and the Eli and Edythe Broad Center of Regenerative Medicine and Stem Cell Research, David Geffen School of Medicine at UCLA, Los Angeles, California, United States of America
| | - Irvin S. Y. Chen
- UCLA AIDS Institute and the Eli and Edythe Broad Center of Regenerative Medicine and Stem Cell Research, David Geffen School of Medicine at UCLA, Los Angeles, California, United States of America
- Department of Microbiology, Immunology, and Molecular Genetics, University of California, Los Angeles, Los Angeles, California, United States of America
| | - Otto O. Yang
- UCLA AIDS Institute and the Eli and Edythe Broad Center of Regenerative Medicine and Stem Cell Research, David Geffen School of Medicine at UCLA, Los Angeles, California, United States of America
- Department of Infectious Disease, David Geffen School of Medicine at UCLA, Los Angeles, California, United States of America
| | - Jerome A. Zack
- Division of Hematology/Oncology, Department of Medicine, David Geffen School of Medicine at UCLA, Los Angeles, California, United States of America
- UCLA AIDS Institute and the Eli and Edythe Broad Center of Regenerative Medicine and Stem Cell Research, David Geffen School of Medicine at UCLA, Los Angeles, California, United States of America
- Department of Microbiology, Immunology, and Molecular Genetics, University of California, Los Angeles, Los Angeles, California, United States of America
| | - Scott G. Kitchen
- Division of Hematology/Oncology, Department of Medicine, David Geffen School of Medicine at UCLA, Los Angeles, California, United States of America
- UCLA AIDS Institute and the Eli and Edythe Broad Center of Regenerative Medicine and Stem Cell Research, David Geffen School of Medicine at UCLA, Los Angeles, California, United States of America
- * E-mail:
| |
Collapse
|
41
|
Sohail MS, Louie RHY, McKay MR, Barton JP. MPL resolves genetic linkage in fitness inference from complex evolutionary histories. Nat Biotechnol 2021; 39:472-479. [PMID: 33257862 PMCID: PMC8044047 DOI: 10.1038/s41587-020-0737-3] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2019] [Accepted: 10/14/2020] [Indexed: 12/13/2022]
Abstract
Genetic linkage causes the fate of new mutations in a population to be contingent on the genetic background on which they appear. This makes it challenging to identify how individual mutations affect fitness. To overcome this challenge, we developed marginal path likelihood (MPL), a method to infer selection from evolutionary histories that resolves genetic linkage. Validation on real and simulated data sets shows that MPL is fast and accurate, outperforming existing inference approaches. We found that resolving linkage is crucial for accurately quantifying selection in complex evolving populations, which we demonstrate through a quantitative analysis of intrahost HIV-1 evolution using multiple patient data sets. Linkage effects generated by variants that sweep rapidly through the population are particularly strong, extending far across the genome. Taken together, our results argue for the importance of resolving linkage in studies of natural selection.
Collapse
Affiliation(s)
- Muhammad Saqib Sohail
- Department of Electronic and Computer Engineering, Hong Kong University of Science and Technology, Hong Kong, China
| | - Raymond H Y Louie
- Department of Electronic and Computer Engineering, Hong Kong University of Science and Technology, Hong Kong, China
- Institute for Advanced Study, Hong Kong University of Science and Technology, Hong Kong, China
- The Kirby Institute, University of New South Wales, Sydney, New South Wales, Australia
- School of Medical Sciences, University of New South Wales, Sydney, New South Wales, Australia
| | - Matthew R McKay
- Department of Electronic and Computer Engineering, Hong Kong University of Science and Technology, Hong Kong, China.
- Department of Chemical and Biological Engineering, Hong Kong University of Science and Technology, Hong Kong, China.
| | - John P Barton
- Department of Physics and Astronomy, University of California, Riverside, Riverside, CA, USA.
| |
Collapse
|
42
|
Ferguson AL, Ranganathan R. 100th Anniversary of Macromolecular Science Viewpoint: Data-Driven Protein Design. ACS Macro Lett 2021; 10:327-340. [PMID: 35549066 DOI: 10.1021/acsmacrolett.0c00885] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
The design of synthetic proteins with the desired function is a long-standing goal in biomolecular science, with broad applications in biochemical engineering, agriculture, medicine, and public health. Rational de novo design and experimental directed evolution have achieved remarkable successes but are challenged by the requirement to find functional "needles" in the vast "haystack" of protein sequence space. Data-driven models for fitness landscapes provide a predictive map between protein sequence and function and can prospectively identify functional candidates for experimental testing to greatly improve the efficiency of this search. This Viewpoint reviews the applications of machine learning and, in particular, deep learning as part of data-driven protein engineering platforms. We highlight recent successes, review promising computational methodologies, and provide an outlook on future challenges and opportunities. The article is written for a broad audience comprising both polymer and protein scientists and computer and data scientists interested in an up-to-date review of recent innovations and opportunities in this rapidly evolving field.
Collapse
Affiliation(s)
- Andrew L. Ferguson
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
| | - Rama Ranganathan
- Pritzker School of Molecular Engineering, University of Chicago, Chicago, Illinois 60637, United States
- Center for Physics of Evolving Systems, University of Chicago, Chicago, Illinois 60637, United States
- Biochemistry and Molecular Biology, University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
43
|
Conti S, Kaczorowski KJ, Song G, Porter K, Andrabi R, Burton DR, Chakraborty AK, Karplus M. Design of immunogens to elicit broadly neutralizing antibodies against HIV targeting the CD4 binding site. Proc Natl Acad Sci U S A 2021; 118:e2018338118. [PMID: 33637649 PMCID: PMC7936365 DOI: 10.1073/pnas.2018338118] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Simone Conti
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA 02138
| | - Kevin J Kaczorowski
- Institute for Medical Engineering & Science, Massachusetts Institute of Technology, Cambridge, MA 02139
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139
| | - Ge Song
- Scripps Consortium for HIV/AIDS Vaccine Development, The Scripps Research Institute, La Jolla, CA 92037
- IAVI Neutralizing Antibody Center, The Scripps Research Institute, La Jolla, CA 92037
- Department of Immunology and Microbiology, The Scripps Research Institute, La Jolla, CA 92037
| | - Katelyn Porter
- Scripps Consortium for HIV/AIDS Vaccine Development, The Scripps Research Institute, La Jolla, CA 92037
- IAVI Neutralizing Antibody Center, The Scripps Research Institute, La Jolla, CA 92037
- Department of Immunology and Microbiology, The Scripps Research Institute, La Jolla, CA 92037
| | - Raiees Andrabi
- Scripps Consortium for HIV/AIDS Vaccine Development, The Scripps Research Institute, La Jolla, CA 92037
- IAVI Neutralizing Antibody Center, The Scripps Research Institute, La Jolla, CA 92037
- Department of Immunology and Microbiology, The Scripps Research Institute, La Jolla, CA 92037
| | - Dennis R Burton
- Scripps Consortium for HIV/AIDS Vaccine Development, The Scripps Research Institute, La Jolla, CA 92037
- IAVI Neutralizing Antibody Center, The Scripps Research Institute, La Jolla, CA 92037
- Department of Immunology and Microbiology, The Scripps Research Institute, La Jolla, CA 92037
- Ragon Institute of Massachusetts General Hospital, Massachusetts Institute of Technology and Harvard University, Cambridge, MA 02139
| | - Arup K Chakraborty
- Institute for Medical Engineering & Science, Massachusetts Institute of Technology, Cambridge, MA 02139;
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139
- Ragon Institute of Massachusetts General Hospital, Massachusetts Institute of Technology and Harvard University, Cambridge, MA 02139
- Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139
| | - Martin Karplus
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA 02138;
- Laboratoire de Chimie Biophysique, Institut de Science et d'Ingénierie Supramoléculaires, Université de Strasbourg, 67000 Strasbourg, France
| |
Collapse
|
44
|
Akand EH, Murray JM. NGlyAlign: an automated library building tool to align highly divergent HIV envelope sequences. BMC Bioinformatics 2021; 22:54. [PMID: 33557755 PMCID: PMC7869453 DOI: 10.1186/s12859-020-03901-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2020] [Accepted: 11/23/2020] [Indexed: 08/29/2023] Open
Abstract
BACKGROUND The high variability in envelope regions of some viruses such as HIV allow the virus to establish infection and to escape subsequent immune surveillance. This variability, as well as increasing incorporation of N-linked glycosylation sites, is fundamental to this evasion. It also creates difficulties for multiple sequence alignment methods (MSA) that provide the first step in their analysis. Existing MSA tools often fail to properly align highly variable HIV envelope sequences requiring extensive manual editing that is impractical with even a moderate number of these variable sequences. RESULTS We developed an automated library building tool NGlyAlign, that organizes similar N-linked glycosylation sites as block constraints and statistically conserved global sites as single site constraints to automatically enforce partial columns in consistency-based MSA methods such as Dialign. This combined method accurately aligns variable HIV-1 envelope sequences. We tested the method on two datasets: a set of 156 founder and chronic gp160 HIV-1 subtype B sequences as well as a set of reference sequences of gp120 in the highly variable region 1. On measures such as entropy scores, sum of pair scores, column score, and similarity heat maps, NGlyAlign+Dialign proved superior against methods such as T-Coffee, ClustalOmega, ClustalW, Praline, HIValign and Muscle. The method is scalable to large sequence sets producing accurate alignments without requiring manual editing. As well as this application to HIV, our method can be used for other highly variable glycoproteins such as hepatitis C virus envelope. CONCLUSIONS NGlyAlign is an automated tool for mapping and building glycosylation motif libraries to accurately align highly variable regions in HIV sequences. It can provide the basis for many studies reliant on single robust alignments. NGlyAlign has been developed as an open-source tool and is freely available at https://github.com/UNSW-Mathematical-Biology/NGlyAlign_v1.0 .
Collapse
Affiliation(s)
- Elma H Akand
- School of Mathematics and Statistics, UNSW, Sydney, NSW, Australia.
| | - John M Murray
- School of Mathematics and Statistics, UNSW, Sydney, NSW, Australia
| |
Collapse
|
45
|
Amitai A. Viral surface geometry shapes influenza and coronavirus spike evolution through antibody pressure. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2020:2020.10.20.347641. [PMID: 33106808 PMCID: PMC7587782 DOI: 10.1101/2020.10.20.347641] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Abstract
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.
Collapse
Affiliation(s)
- Assaf Amitai
- Institute for Medical Engineering and Science, Massachusetts Institute of Technology, Cambridge, MA 02139
- Ragon Institute of MGH, MIT, and Harvard, Cambridge, MA 02139
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139
| |
Collapse
|
46
|
Akand EH, Maher SJ, Murray JM. Mutational networks of escape from transmitted HIV-1 infection. PLoS One 2020; 15:e0243391. [PMID: 33284837 PMCID: PMC7721145 DOI: 10.1371/journal.pone.0243391] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2020] [Accepted: 11/19/2020] [Indexed: 02/08/2023] Open
Abstract
Human immunodeficiency virus (HIV) is subject to immune selective pressure soon after it establishes infection at the founder stage. As an individual progresses from the founder to chronic stage of infection, immune pressure forces a history of mutations that are embedded in envelope sequences. Determining this pathway of coevolving mutations can assist in understanding what is different with the founder virus and the essential pathways it takes to maintain infection. We have combined operations research and bioinformatics methods to extract key networks of mutations that differentiate founder and chronic stages for 156 subtype B and 107 subtype C envelope (gp160) sequences. The chronic networks for both subtypes revealed strikingly different hub-and-spoke topologies compared to the less structured transmission networks. This suggests that the hub nodes are impacted by the immune response and the resulting loss of fitness is compensated by mutations at the spoke positions. The major hubs in the chronic C network occur at positions 12, 137 (within the N136 glycan), and 822, and at position 306 for subtype B. While both founder networks had a more heterogeneous connected network structure, interestingly founder B subnetworks around positions 640 and 837 preferentially contained CD4 and coreceptor binding domains. Finally, we observed a differential effect of glycosylation between founder and chronic subtype B where the latter had mutational pathways significantly driven by N-glycosylation. Our study provides insights into the mutational pathways HIV takes to evade the immune response, and presents features more likely to establish founder infection, valuable for effective vaccine design.
Collapse
Affiliation(s)
- Elma H. Akand
- School of Mathematics and Statistics, UNSW Sydney, Kensington, NSW, Australia
| | - Stephen J. Maher
- College of Engineering, Mathematical and Physical Sciences, University of Exeter, Exeter, United Kingdom
| | - John M. Murray
- School of Mathematics and Statistics, UNSW Sydney, Kensington, NSW, Australia
| |
Collapse
|
47
|
Zhang TH, Dai L, Barton JP, Du Y, Tan Y, Pang W, Chakraborty AK, Lloyd-Smith JO, Sun R. Predominance of positive epistasis among drug resistance-associated mutations in HIV-1 protease. PLoS Genet 2020; 16:e1009009. [PMID: 33085662 PMCID: PMC7605711 DOI: 10.1371/journal.pgen.1009009] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2019] [Revised: 11/02/2020] [Accepted: 07/24/2020] [Indexed: 12/12/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Tian-hao Zhang
- Molecular Biology Institute, University of California, Los Angeles, CA 90095, USA
| | - Lei Dai
- CAS Key Laboratory of Quantitative Engineering Biology, Shenzhen Institute of Synthetic Biology, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen 518055, China
| | - John P. Barton
- Department of Physics and Astronomy, University of California, Riverside, CA 92521, USA
| | - Yushen Du
- School of Medicine, ZheJiang University, Hangzhou, 210000, China
- Molecular and Medical Pharmacology, University of California, Los Angeles, CA 90095, USA
| | - Yuxiang Tan
- CAS Key Laboratory of Quantitative Engineering Biology, Shenzhen Institute of Synthetic Biology, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen 518055, China
| | - Wenwen Pang
- Department of Public Health Laboratory Science, West China School of Public Health, Sichuan University, Chengdu 610041, China
| | - Arup K. Chakraborty
- Institute for Medical Engineering and Science, Departments of Chemical Engineering, Physics, & Chemistry, Massachusetts Institute of Technology, MA 21309, USA
- Ragon Institute of MGH, MIT, & Harvard, Cambridge, MA 21309, USA
| | - James O. Lloyd-Smith
- Department of Ecology and Evolutionary Biology, University of California, Los Angeles, CA 90095, USA
| | - Ren Sun
- Molecular and Medical Pharmacology, University of California, Los Angeles, CA 90095, USA
| |
Collapse
|
48
|
A de novo approach to inferring within-host fitness effects during untreated HIV-1 infection. PLoS Pathog 2020; 16:e1008171. [PMID: 32492061 PMCID: PMC7295245 DOI: 10.1371/journal.ppat.1008171] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2019] [Revised: 06/15/2020] [Accepted: 05/11/2020] [Indexed: 12/15/2022] Open
Abstract
In the absence of effective antiviral therapy, HIV-1 evolves in response to the within-host environment, of which the immune system is an important aspect. During the earliest stages of infection, this process of evolution is very rapid, driven by a small number of CTL escape mutations. As the infection progresses, immune escape variants evolve under reduced magnitudes of selection, while competition between an increasing number of polymorphic alleles (i.e., clonal interference) makes it difficult to quantify the magnitude of selection acting upon specific variant alleles. To tackle this complex problem, we developed a novel multi-locus inference method to evaluate the role of selection during the chronic stage of within-host infection. We applied this method to targeted sequence data from the p24 and gp41 regions of HIV-1 collected from 34 patients with long-term untreated HIV-1 infection. We identify a broad distribution of beneficial fitness effects during infection, with a small number of variants evolving under strong selection and very many variants evolving under weaker selection. The uniquely large number of infections analysed granted a previously unparalleled statistical power to identify loci at which selection could be inferred to act with statistical confidence. Our model makes no prior assumptions about the nature of alleles under selection, such that any synonymous or non-synonymous variant may be inferred to evolve under selection. However, the majority of variants inferred with confidence to be under selection were non-synonymous in nature, and in most cases were have previously been associated with either CTL escape in p24 or neutralising antibody escape in gp41. We also identified a putative new CTL escape site (residue 286 in gag), and a region of gp41 (including residues 644, 648, 655 in env) likely to be associated with immune escape. Sites inferred to be under selection in multiple hosts have high within-host and between-host diversity although not all sites with high between-host diversity were inferred to be under selection at the within-host level. Our identification of selection at sites associated with resistance to broadly neutralising antibodies (bNAbs) highlights the need to fully understand the role of selection in untreated individuals when designing bNAb based therapies.
Collapse
|
49
|
Gao A, Chen Z, Segal FP, Carrington M, Streeck H, Chakraborty AK, Julg B. Predicting the Immunogenicity of T cell epitopes: From HIV to SARS-CoV-2. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2020:2020.05.14.095885. [PMID: 32511339 PMCID: PMC7241102 DOI: 10.1101/2020.05.14.095885] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
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.
Collapse
Affiliation(s)
- Ang Gao
- Institute for Medical Engineering & Science, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
- Department of Chemical Engineering, Massachusetts Institute of Technology Cambridge, MA 02139, USA
| | - Zhilin Chen
- Ragon Insitute of Massachusetts General Hospital, Massachusetts Institute of Technology, and Harvard, Cambridge, MA 02139, USA
| | | | - Mary Carrington
- Ragon Insitute of Massachusetts General Hospital, Massachusetts Institute of Technology, and Harvard, Cambridge, MA 02139, USA
- Basic Science Program, Frederick National Laboratory for Cancer Research, Frederick, MD 21702, USA
| | - Hendrik Streeck
- Institut für Virologie, Universitätsklinikum Bonn, 53127 Bonn, Germany
| | - Arup K. Chakraborty
- Institute for Medical Engineering & Science, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
- Ragon Insitute of Massachusetts General Hospital, Massachusetts Institute of Technology, and Harvard, Cambridge, MA 02139, USA
- Department of Chemical Engineering, Massachusetts Institute of Technology Cambridge, MA 02139, USA
- Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Boris Julg
- Ragon Insitute of Massachusetts General Hospital, Massachusetts Institute of Technology, and Harvard, Cambridge, MA 02139, USA
| |
Collapse
|
50
|
A model for the interplay between plastic tradeoffs and evolution in changing environments. Proc Natl Acad Sci U S A 2020; 117:8934-8940. [PMID: 32245811 DOI: 10.1073/pnas.1915537117] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
Performance tradeoffs are ubiquitous in both ecological and evolutionary modeling, yet they are usually postulated and built into fitness and ecological landscapes. However, tradeoffs depend on genetic background and evolutionary history and can themselves evolve. We present a simple model capable of capturing the key feedback loop: evolutionary history shapes tradeoff strength, which, in turn, shapes evolutionary future. One consequence of this feedback is that genomes with identical fitness can have different evolutionary properties shaped by prior environmental exposure. Another is that, generically, the best adaptations to one environment may evolve in another. Our simple framework bridges the gap between the phenotypic Fisher's Geometric Model and the genotypic properties, such as modularity and evolvability, and can serve as a rich playground for investigating evolution in multiple or changing environments.
Collapse
|