Predicting potential SARS-CoV-2 mutations of concern via full quantum mechanical modelling.

Ab initio quantum mechanical models can characterize and predict intermolecular binding, but only recently have models including more than a few hundred atoms gained traction. Here, we simulate the electronic structure for approximately 13 000 atoms to predict and characterize binding of SARS-CoV-2 spike variants to the human ACE2 (hACE2) receptor using the quantum mechanics complexity reduction (QM-CR) approach. We compare four spike variants in our analysis: Wuhan, Omicron, and two Omicron-based variants. To assess binding, we mechanistically characterize the energetic contribution of each amino acid involved, and predict the effect of select single amino acid mutations. We validate our computational predictions experimentally by comparing the efficacy of spike variants binding to cells expressing hACE2. At the time we performed our simulations (December 2021), the mutation A484K which our model predicted to be highly beneficial to ACE2 binding had not been identified in epidemiological surveys; only recently (August 2023) has it appeared in variant BA.2.86. We argue that our computational model, QM-CR, can identify mutations critical for intermolecular interactions and inform the engineering of high-specificity interactors.

Reverse transcriptase inhibitors prevent liver abscess formation during Escherichia coli bloodstream infection.

The presence of bacteria in the bloodstream is associated with severe clinical outcomes. In mice, intravenous inoculation of Escherichia coli can lead to the formation of macroscopic abscesses in the liver. Abscesses are regions of severe necrosis and consist of millions of bacteria surrounded by inflammatory immune cells. Liver abscess susceptibility varies widely across strains of mice, but the host factors governing this variation are unknown. Here, we profiled hepatic transcriptomes in mice with varying susceptibility to liver abscess formation. We found that transcripts from endogenous retroviruses (ERVs) are robustly induced in the liver by E. coli infection and ERV expression positively correlates with the frequency of abscess formation. Hypothesizing that ERV-encoded reverse transcriptase may generate cytoplasmic DNA and heighten inflammatory responses, we tested whether nucleoside/nucleotide reverse transcriptase inhibitors (NRTIs) influence abscess formation. Strikingly, a single NRTI dose administered immediately following E. coli inoculation prevented abscess formation, leading to a concomitant 100,000-fold reduction in bacterial burden. We provide evidence that NRTIs inhibit abscess formation by preventing the tissue necrosis that facilitates bacterial replication. Together, our findings suggest that endogenous reverse transcriptases drive inflammatory responses during bacterial bloodstream infection to drive abscess formation. The high efficacy of NRTIs in preventing abscess formation suggests that the consequences of reverse transcription on inflammation should be further examined, particularly in infectious diseases where inflammation drives negative clinical outcomes, such as sepsis.

Unique Structure and Distinctive Properties of the Ancient and Ubiquitous Gamma-Type Envelope Glycoprotein.

After the onset of the AIDS pandemic, HIV-1 (genus Lentivirus) became the predominant model for studying retrovirus Env glycoproteins and their role in entry. However, HIV Env is an inadequate model for understanding entry of viruses in the Alpharetrovirus, Gammaretrovirus and Deltaretrovirus genera. For example, oncogenic model system viruses such as Rous sarcoma virus (RSV, Alpharetrovirus), murine leukemia virus (MLV, Gammaretrovirus) and human T-cell leukemia viruses (HTLV-I and HTLV-II, Deltaretrovirus) encode Envs that are structurally and functionally distinct from HIV Env. We refer to these as Gamma-type Envs. Gamma-type Envs are probably the most widespread retroviral Envs in nature. They are found in exogenous and endogenous retroviruses representing a broad spectrum of vertebrate hosts including amphibians, birds, reptiles, mammals and fish. In endogenous form, gamma-type Envs have been evolutionarily coopted numerous times, most notably as placental syncytins (e.g., human SYNC1 and SYNC2). Remarkably, gamma-type Envs are also found outside of the Retroviridae. Gp2 proteins of filoviruses (e.g., Ebolavirus) and snake arenaviruses in the genus Reptarenavirus are gamma-type Env homologs, products of ancient recombination events involving viruses of different Baltimore classes. Distinctive hallmarks of gamma-type Envs include a labile disulfide bond linking the surface and transmembrane subunits, a multi-stage attachment and fusion mechanism, a highly conserved (but poorly understood) "immunosuppressive domain", and activation by the viral protease during virion maturation. Here, we synthesize work from diverse retrovirus model systems to illustrate these distinctive properties and to highlight avenues for further exploration of gamma-type Env structure and function.

Probing the mutational landscape of the SARS-CoV-2 spike protein via quantum mechanical modeling of crystallographic structures.

We employ a recently developed complexity-reduction quantum mechanical (QM-CR) approach, based on complexity reduction of density functional theory calculations, to characterize the interactions of the SARS-CoV-2 spike receptor binding domain (RBD) with ACE2 host receptors and antibodies. QM-CR operates via ab initio identification of individual amino acid residue's contributions to chemical binding and leads to the identification of the impact of point mutations. Here, we especially focus on the E484K mutation of the viral spike protein. We find that spike residue 484 hinders the spike's binding to the human ACE2 receptor (hACE2). In contrast, the same residue is beneficial in binding to the bat receptor Rhinolophus macrotis ACE2 (macACE2). In agreement with empirical evidence, QM-CR shows that the E484K mutation allows the spike to evade categories of neutralizing antibodies like C121 and C144. The simulation also shows how the Delta variant spike binds more strongly to hACE2 compared to the original Wuhan strain, and predicts that a E484K mutation can further improve its binding. Broad agreement between the QM-CR predictions and experimental evidence supports the notion that ab initio modeling has now reached the maturity required to handle large intermolecular interactions central to biological processes.

ICTV Virus Taxonomy Profile: Retroviridae 2021.

Viruses in the family Retroviridae are found in a wide variety of vertebrate hosts. Enveloped virions are 80-100 nm in diameter with an inner core containing the viral genome and replicative enzymes. Core morphology is often characteristic for viruses within the same genus. Replication involves reverse transcription and integration into host cell DNA, resulting in a provirus. Integration into germline cells can result in a heritable provirus known as an endogenous retrovirus. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the family Retroviridae, which is available at ictv.global/report/retroviridae.

APOBEC3F Constitutes a Barrier to Successful Cross-Species Transmission of Simian Immunodeficiency Virus SIVsmm to Humans.

Simian immunodeficiency virus infecting sooty mangabeys (SIVsmm) has been transmitted to humans on at least nine occasions, giving rise to human immunodeficiency virus type 2 (HIV-2) groups A to I. SIVsmm isolates replicate in human T cells and seem capable of overcoming major human restriction factors without adaptation. However, only groups A and B are responsible for the HIV-2 epidemic in sub-Saharan Africa, and it is largely unclear whether adaptive changes were associated with spread in humans. To address this, we examined the sensitivity of infectious molecular clones (IMCs) of five HIV-2 strains and representatives of five different SIVsmm lineages to various APOBEC3 proteins. We confirmed that SIVsmm strains replicate in human T cells, albeit with more variable replication fitness and frequently lower efficiency than HIV-2 IMCs. Efficient viral propagation was generally dependent on intact vif genes, highlighting the need for counteraction of APOBEC3 proteins. On average, SIVsmm was more susceptible to inhibition by human APOBEC3D, -F, -G, and -H than HIV-2. For example, human APOBEC3F reduced infectious virus yield of SIVsmm by ∼80% but achieved only ∼40% reduction in the case of HIV-2. Functional and mutational analyses of human- and monkey-derived alleles revealed that an R128T polymorphism in APOBEC3F contributes to species-specific counteraction by HIV-2 and SIVsmm Vifs. In addition, a T84S substitution in SIVsmm Vif increased its ability to counteract human APOBEC3F. Altogether, our results confirm that SIVsmm Vif proteins show intrinsic activity against human APOBEC3 proteins but also demonstrate that epidemic HIV-2 strains evolved an increased ability to counteract this class of restriction factors during human adaptation. IMPORTANCE Viral zoonoses pose a significant threat to human health, and it is important to understand determining factors. SIVs infecting great apes gave rise to HIV-1. In contrast, SIVs infecting African monkey species have not been detected in humans, with one notable exception. SIVsmm from sooty mangabeys has crossed the species barrier to humans on at least nine independent occasions and seems capable of overcoming many innate defense mechanisms without adaptation. Here, we confirmed that SIVsmm Vif proteins show significant activity against human APOBEC3 proteins. Our analyses also revealed, however, that different lineages of SIVsmm are significantly more susceptible to inhibition by various human APOBEC3 proteins than HIV-2 strains. Mutational analyses suggest that an R128T substitution in APOBEC3F and a T84S change in Vif contribute to species-specific counteraction by HIV-2 and SIVsmm. Altogether, our results support that epidemic HIV-2 strains acquired increased activity against human APOBEC3 proteins to clear this restrictive barrier.

Structure of the Receptor Binding Domain of EnvP(b)1, an Endogenous Retroviral Envelope Protein Expressed in Human Tissues.

EnvP(b)1 is an endogenous retroviral envelope gene found in human and other primate genomes. We report EnvP(b)1 sequences in primate genomes consistent with an integration event between 40 and 71 million years ago. Using a highly specific polyclonal antiserum raised against the putative receptor binding domain (RBD) of human EnvP(b)1, we detected expression in human placenta, ovaries, and thymus. We found that EnvP(b)1 is proteolytically processed, and using cell-cell fusion assays in multiple primate cell lines, we demonstrated that extant EnvP(b)1 proteins from a variety of primate genomes are fusogenic. This work supports the idea that EnvP(b)1 is under purifying selection and its fusogenic activity has been maintained for over 40 million years. We determined the structure of the RBD of human EnvP(b)1, which defines structural similarities with extant leukemia viruses, despite little sequence conservation. This structure highlights a common scaffold from which novel receptor binding specificities likely evolved. The evolutionary plasticity of this domain may underlie the diversity of related Envs in circulating viruses.IMPORTANCE Organisms can access genetic and functional novelty by capturing viral elements within their genomes, where they can evolve to drive new cellular or organismal processes. We demonstrate that a retroviral envelope gene, EnvP(b)1, has been maintained and its fusion activity preserved for 40 to 71 million years. It is expressed as a protein in multiple healthy human tissues. We determined the structure of its inferred receptor binding domain and compared it with the same domain in modern viruses. We found a common conserved architecture that underlies the varied receptor binding activity of divergent Env genes. The modularity and versatility of this domain may underpin the evolutionary success of this clade of fusogens.

Not Your Typical Anti-CRISPR.

In this issue of Cell Host & Microbe, two papers by Osuna et al. describe the characterization of AcrIIA1, an anti-CRISPR protein distributed widely among Listeria phages. AcrIIA1 functions as an anti-CRISPR and as a dynamic repressor of acr loci, suggesting it may play an important role in lysogeny.

Complex Problems Require Ternary Solutions: Another Lesson from SIV Nef.

In this issue of Cell Host & Microbe, Buffalo et al. describe a cryo-EM structure of SIV Nef complexed with AP-2 and tetherin. The structure helps explain why human tetherin is Nef-resistant and why lentiviruses that successfully emerged in humans (HIV-1 and HIV-2) had to evolve novel anti-tetherin strategies.

Origins and evolutionary consequences of ancient endogenous retroviruses.

Retroviruses infect a broad range of vertebrate hosts that includes amphibians, reptiles, fish, birds and mammals. In addition, a typical vertebrate genome contains thousands of loci composed of ancient retroviral sequences known as endogenous retroviruses (ERVs). ERVs are molecular remnants of ancient retroviruses and proof that the ongoing relationship between retroviruses and their vertebrate hosts began hundreds of millions of years ago. The long-term impact of retroviruses on vertebrate evolution is twofold: first, as with other viruses, retroviruses act as agents of selection, driving the evolution of host genes that block viral infection or that mitigate pathogenesis, and second, through the phenomenon of endogenization, retroviruses contribute an abundance of genetic novelty to host genomes, including unique protein-coding genes and cis-acting regulatory elements. This Review describes ERV origins, their diversity and their relationships to retroviruses and discusses the potential for ERVs to reveal virus-host interactions on evolutionary timescales. It also describes some of the many examples of cellular functions, including protein-coding genes and regulatory elements, that have evolved from ERVs.

Nomenclature for endogenous retrovirus (ERV) loci.

Retroviral integration into germline DNA can result in the formation of a vertically inherited proviral sequence called an endogenous retrovirus (ERV). Over the course of their evolution, vertebrate genomes have accumulated many thousands of ERV loci. These sequences provide useful retrospective information about ancient retroviruses, and have also played an important role in shaping the evolution of vertebrate genomes. There is an immediate need for a unified system of nomenclature for ERV loci, not only to assist genome annotation, but also to facilitate research on ERVs and their impact on genome biology and evolution. In this review, we examine how ERV nomenclatures have developed, and consider the possibilities for the implementation of a systematic approach for naming ERV loci. We propose that such a nomenclature should not only provide unique identifiers for individual loci, but also denote orthologous relationships between ERVs in different species. In addition, we propose that-where possible-mnemonic links to previous, well-established names for ERV loci and groups should be retained. We show how this approach can be applied and integrated into existing taxonomic and nomenclature schemes for retroviruses, ERVs and transposable elements.

Spumaretroviruses: Updated taxonomy and nomenclature.

Spumaretroviruses, commonly referred to as foamy viruses, are complex retroviruses belonging to the subfamily Spumaretrovirinae, family Retroviridae, which naturally infect a variety of animals including nonhuman primates (NHPs). Additionally, cross-species transmissions of simian foamy viruses (SFVs) to humans have occurred following exposure to tissues of infected NHPs. Recent research has led to the identification of previously unknown exogenous foamy viruses, and to the discovery of endogenous spumaretrovirus sequences in a variety of host genomes. Here, we describe an updated spumaretrovirus taxonomy that has been recently accepted by the International Committee on Taxonomy of Viruses (ICTV) Executive Committee, and describe a virus nomenclature that is generally consistent with that used for other retroviruses, such as lentiviruses and deltaretroviruses. This taxonomy can be applied to distinguish different, but closely related, primate (e.g., human, ape, simian) foamy viruses as well as those from other hosts. This proposal accounts for host-virus co-speciation and cross-species transmission.

High-Resolution Sequencing of Viral Populations during Early Simian Immunodeficiency Virus Infection Reveals Evolutionary Strategies for Rapid Escape from Emerging Env-Specific Antibody Responses.

Primate lentiviruses, including the human and simian immunodeficiency viruses (HIV and SIV), produce infections marked by persistent, ongoing viral replication. This occurs despite the presence of virus-specific adaptive immune responses, including antibodies targeting the viral envelope glycoprotein (Env), and evolution of antibody-escape variants is a well-documented feature of lentiviral infection. Here, we examined the evolutionary dynamics of the SIV env gene during early infection (≤29 weeks postinfection) in a cohort of four SIVmac251-infected rhesus macaques. We tracked env evolution during acute and early infection using frequent sampling and ultradeep sequencing of viral populations, capturing a transmission bottleneck and the subsequent reestablishment of Env diversity. A majority of changes in the gp120 subunit mapped to two short clusters, one in the first variable region (V1) and one in V4, while most changes in the gp41 subunit appeared in the cytoplasmic domain. Variation in V1 was dominated by short duplications and deletions of repetitive sequence, while variation in V4 was marked by short in-frame deletions and closely overlapping substitutions. The most common substitutions in both patches did not alter viral replicative fitness when tested using a highly sensitive, deep-sequencing-based competition assay. Our results, together with the observation that very similar or identical patterns of sequence evolution also occur in different macaque species infected with related but divergent strains of SIV, suggest that resistance to early, strain-specific anti-Env antibodies is the result of temporally and mutationally predictable pathways of escape that occur during the early stages of infection.IMPORTANCE The envelope glycoprotein (Env) of primate lentiviruses mediates entry by binding to host cell receptors followed by fusion of the viral membrane with the cell membrane. The exposure of Env complexes on the surface of the virion results in targeting by antibodies, leading to selection for virus escape mutations. We used the SIV/rhesus macaque model to track in vivo evolution of variation in Env during acute/early infection in animals with and without antibody responses to Env, uncovering remarkable variation in animals with antibody responses within weeks of infection. Using a deep-sequencing-based fitness assay, we found substitutions associated with antibody escape had little to no effect on inherent replicative capacity. The ability to readily propagate advantageous changes that incur little to no replicative fitness costs may be a mechanism to maintain continuous replication under constant immune selection, allowing the virus to persist for months to years in the infected host.

Retroviruses of the RDR superinfection interference group: ancient origins and broad host distribution of a promiscuous Env gene.

Due to recombination, different regions of a retrovirus genome can have distinct phylogenetic histories. The RD114-and-D-type-retrovirus (RDR) interference group provides an extreme example: the RDR group comprises a variety of taxonomically distinct retroviruses, isolated from diverse mammalian and avian hosts, that share a homologous env gene and use the same cell-surface entry receptor. RDR env homologs are also found among ancient endogenous retrovirus (ERV) sequences, including the syncytin genes of humans and rabbits, indicating that RDR Env glycoproteins have likely mediated endogenization on multiple occasions in diverse vertebrate lineages. The distribution of RDR env among exogenous and endogenous retroviruses indicates that it has been swapped between viruses many times, and that it likely facilitated multiple cross-species transmission events spanning millions of years of vertebrate evolution.

A penalized regression approach to haplotype reconstruction of viral populations arising in early HIV/SIV infection.

MOTIVATION: Next generation sequencing (NGS) has been increasingly applied to characterize viral evolution during HIV and SIV infections. In particular, NGS datasets sampled during the initial months of infection are characterized by relatively low levels of diversity as well as convergent evolution at multiple loci dispersed across the viral genome. Consequently, fully characterizing viral evolution from NGS datasets requires haplotype reconstruction across large regions of the viral genome. Existing haplotype reconstruction algorithms have not been developed with the particular characteristics of early HIV/SIV infection in mind, raising the possibility that better performance could be achieved through a specifically designed algorithm. RESULTS: Here, we introduce a haplotype reconstruction algorithm, RegressHaplo, specifically designed for low diversity and convergent evolution regimes. The algorithm uses a penalized regression that balances a data fitting term with a penalty term that encourages solutions with few haplotypes. The regression covariates are a large set of potential haplotypes and fitting the regression is made computationally feasible by the low diversity setting. Using simulated and in vivo datasets, we compare RegressHaplo to PredictHaplo and QuRe, two existing haplotype reconstruction algorithms. RegressHaplo performs better than these algorithms on simulated datasets with relatively low diversity levels. We suggest RegressHaplo as a novel tool for the investigation of early infection HIV/SIV datasets and, more generally, low diversity viral NGS datasets. CONTACT: sr286@georgetown.edu. AVAILABILITY AND IMPLEMENTATION: https://github.com/SLeviyang/RegressHaplo.

Analysis of SIVmac Envelope-Specific Antibodies Selected Through Phage Display.

We have constructed a single chain fragment variable (scFv) phage display library from a simian immunodeficiency virus (SIV)-infected rhesus macaque that developed unusually high-titer neutralizing antibody responses against tier-3, neutralization-resistant SIVmac239. The library was screened using trimeric (gp140) and monomeric (gp120) forms of the SIVmac239 envelope (Env) glycoprotein. We also cloned variable-heavy and variable-light (VH-VL) antibody fragments from seven previously described rhesus macaque B-cell lines (BLCLs) that produce SIV gp120-specific monoclonal antibodies (mAbs). Thirty-two gp140-specific mAbs were selected along with 20 gp120-specific ones. gp120-specific mAbs were only from the VH4 family, while gp41-specific mAbs were primarily from VH1, followed by VH4 and VH3. Rhesus macaque BLCL-derived mAbs belonged primarily to the VH4 family of antibodies followed by VH3 and a smaller number of VH1s. A preferential VH combination with Vλ light chain was observed with phage display-selected SIV Env-specific mAbs (gp120 and gp140), but not with BLCL-derived antibodies or the unpanned library. None of the tested antibodies had detectable neutralizing activity against tier-3 SIVmac239. The majority of gp120-specifc mAbs potently neutralized tier-1 SIVmac316 with 50% inhibitory concentration (IC50) values below 1 µg/ml. For gp140-specific antibodies, which were all specific for the gp41-subunit, 2 out of 11 tested neutralized SIVmac316 (IC50 of 7 and 5 µg/ml, respectively). These data suggest an order of preferential VH segment usage for SIV-specific antibodies in rhesus macaques. These antibodies will be useful in assessing the contribution of non-neutralizing antibodies to inhibition of SIV infection in vitro and in vivo.

An Intact Retroviral Gene Conserved in Spiny-Rayed Fishes for over 100 My.

We have identified a retroviral envelope gene with a complete, intact open reading frame (ORF) in 20 species of spiny-rayed fishes (Acanthomorpha). The taxonomic distribution of the gene, "percomORF", indicates insertion into the ancestral lineage >110 Ma, making it the oldest known conserved gene of viral origin in a vertebrate genome. Underscoring its ancient provenence, percomORF exists as an isolated ORF within the intron of a widely conserved host gene, with no discernible proviral sequence nearby. Despite its remarkable age, percomORF retains canonical features of a retroviral glycoprotein, and tests for selection strongly suggest cooption for a host function. Retroviral envelope genes have been coopted for a role in placentogenesis by numerous lineages of mammals, including eutherians and marsupials, representing a variety of placental structures. Therefore percomORF's presence within the group Percomorpha-unique among spiny-finned fishes in having evolved placentation and live birth-is especially intriguing.

TRIM5α Resistance Escape Mutations in the Capsid Are Transferable between Simian Immunodeficiency Virus Strains.

TRIM5α polymorphism limits and complicates the use of simian immunodeficiency virus (SIV) for evaluation of human immunodeficiency virus (HIV) vaccine strategies in rhesus macaques. We previously reported that the TRIM5α-sensitive SIV from sooty mangabeys (SIVsm) clone SIVsmE543-3 acquired amino acid substitutions in the capsid that overcame TRIM5α restriction when it was passaged in rhesus macaques expressing restrictive TRIM5α alleles. Here we generated TRIM5α-resistant clones of the related SIVsmE660 strain without animal passage by introducing the same amino acid capsid substitutions. We evaluated one of the variants in rhesus macaques expressing permissive and restrictive TRIM5α alleles. The SIVsmE660 variant infected and replicated in macaques with restrictive TRIM5α genotypes as efficiently as in macaques with permissive TRIM5α genotypes. These results demonstrated that mutations in the SIV capsid can confer SIV resistance to TRIM5α restriction without animal passage, suggesting an applicable method to generate more diverse SIV strains for HIV vaccine studies. IMPORTANCE: Many strains of SIV from sooty mangabey monkeys are susceptible to resistance by common rhesus macaque TRIM5α alleles and result in reduced virus acquisition and replication in macaques that express these restrictive alleles. We previously observed that spontaneous variations in the capsid gene were associated with improved replication in macaques, and the introduction of two amino acid changes in the capsid transfers this improved replication to the parent clone. In the present study, we introduced these mutations into a related but distinct strain of SIV that is commonly used for challenge studies for vaccine trials. These mutations also improved the replication of this strain in macaques with the restrictive TRIM5α genotype and thus will eliminate the confounding effects of TRIM5α in vaccine studies.