AUTHOR SUMMARY Phosphatidylserine (PS) receptors are PS binding proteins that mediate uptake of apoptotic bodies. Many enveloped viruses utilize this PS/PS receptor mechanism to adhere to and internalize into the endosomal compartment of cells and this is termed apoptotic mimicry. For viruses that have a mechanism(s) of endosomal escape, apoptotic mimicry is a productive route of virus entry. We evaluated if PS receptors serve as cell surface receptors for SARS-CoV-2 and found that the PS receptors, AXL, TIM-1 and TIM-4, facilitated virus infection when low concentrations of the SARS-CoV-2 cognate receptor, ACE2, was present. Consistent with the established mechanism of PS receptor utilization by other viruses, PS liposomes competed with SARS-CoV-2 for binding and entry. We demonstrated that this PS receptor enhances SARS-CoV-2 binding to and infection of an array of human lung cell lines and is an under-appreciated but potentially important host factor facilitating SARS-CoV-2 entry.

ABSTRACT Therapeutics that target the envelope glycoproteins (Envs) of human immunodeficiency virus type 1 (HIV-1) effectively reduce virus levels in patients. However, due to mutations, new Env variants are frequently generated, which may be resistant to the treatments. The appearance of such sequence variance at any Env position is seemingly random. A better understanding of the spatiotemporal patterns of variance across Env may lead to the development of new therapeutic strategies. We hypothesized that, at any time point in a patient, positions with sequence variance are clustered on the three-dimensional structure of Env. To test this hypothesis, we examined whether variance at any Env position can be predicted by the variance measured at adjacent positions. Sequences from 300 HIV-infected patients were applied to a new algorithm we developed. The k-best classifiers (KBC) method is a dynamic ensemble selection technique that identifies the best classifier(s) within the neighborhood of a new observation. It applies bootstrap resampling to generate out-of-bag samples that are used with the resampled set to evaluate each classifier. For many positions of Env, primarily in the CD4-binding site, KBC accurately predicted variance based on the variance at their adjacent positions. KBC improved performance compared to the initial learners, static ensemble, and other baseline models. KBC also outperformed other algorithms for predicting variance at multi-position footprints of therapeutics on Env. These understandings can be applied to refine models that predict future changes in HIV-1 Env. More generally, we propose KBC as a new high-performance dynamic ensemble selection technique.

New sublineages of SARS-CoV-2 variants-of-concern (VOCs) continuously emerge with mutations in the spike glycoprotein. In most cases, the sublineage-defining mutations vary between the VOCs. It is unclear whether these differences reflect lineage-specific likelihoods for mutations at each spike position or the stochastic nature of their appearance. Here we show that SARS-CoV-2 lineages have distinct evolutionary spaces (a probabilistic definition of the sequence states that can be occupied by expanding virus subpopulations). This space can be accurately inferred from the patterns of amino acid variability at the whole-protein level. Robust networks of co-variable sites identify the highest-likelihood mutations in new VOC sublineages and predict remarkably well the emergence of subvariants with resistance mutations to COVID-19 therapeutics. Our studies reveal the contribution of low frequency variant patterns at heterologous sites across the protein to accurate prediction of the changes at each position of interest.

ABSTRACT The error-prone replication machinery of HIV-1 continuously generates new variants of the envelope glycoproteins (Envs). Antibody selection pressures applied in the host can limit their persistence. The target specificity of antibodies elicited in different hosts varies considerably. Whether some specificities are shared and have affected the population-level evolution of Env structure is still unclear. We examined the historical changes in amino acid sequence of the gp41 fusion peptide proximal region (FPPR), which is not exposed on the Env trimer. For three FPPR positions, the residue found in the clade B ancestor was mainly replaced by alanine. However, the changes in alanine frequency at these positions between 1979 and 2016 followed different patterns; two positions maintained a historically-constant frequency whereas the third showed a gradual increase. To understand these patterns, we introduced alanine substitutions in the FPPR of primary HIV-1 strains and examined their fitness and antigenicity relative to the clade-ancestral form. The evolutionary patterns could not be explained by effects on Env fitness. Instead, the FPPR variants with a historically-constant alanine frequency exhibited a unique open-at-the-base conformation of the trimer that exposes partially-cryptic epitopes. These Envs were modestly but significantly more sensitive to poorly-neutralizing sera from HIV-infected individuals than the clade-ancestral form. Our findings suggest that weakly-neutralizing antibodies targeting the base of the trimer are commonly elicited. Such low-level antibody pressures do not exert catastrophic effects on the emerging variants but rather determine their set-point frequency in the population and historical patterns of change. IMPORTANCE HIV-1 infection elicits antibodies that target the Env proteins of the virus. The specific targets of these antibodies vary between infected individuals. It is unclear whether some target specificities are shared between the antibody responses of different individuals. Our data suggest that antibodies against the base of the Env protein are commonly elicited during infection and are contained in sera with low neutralization efficacy. Such antibody pressures are weak. As a result, they do not completely eliminate the sensitive Env forms from the population, but rather maintain their frequency at a low level that has not increased during the past 40 years.

ABSTRACT The spike protein of SARS-CoV-2 is arranged as a trimer on the virus surface, composed of three S1 and three S2 subunits. Infected and vaccinated individuals generate antibodies against spike, which can neutralize the virus. Most antibodies target the receptor-binding domain (RBD) and N-terminal domain (NTD) of S1; however, antibodies against other regions of spike have also been isolated. The variation between infected individuals in domain specificity of the antibodies and in their relative neutralization efficacy is still poorly characterized. To this end, we tested serum and plasma samples from 85 COVID-19 convalescent subjects using 7 immunoassays that employ different domains, subunits and oligomeric forms of spike to capture the antibodies. Samples were also tested for their neutralization of pseudovirus containing SARS-CoV-2 spike and of replication-competent SARS-CoV-2. We observed strong correlations between the levels of NTD- and RBD-specific antibodies, with a fixed ratio of each type to all anti-spike antibodies. The relative potency of the response (defined as the measured neutralization efficacy relative to the total level of spike-targeting antibodies) also exhibited limited variation between subjects, and was not associated with the overall amount of anti-spike antibodies produced. Accordingly, the ability of immunoassays that use RBD, NTD and different forms of S1 or S1/S2 as capture antigens to estimate the neutralizing efficacy of convalescent samples was largely similar. These studies suggest that host-to-host variation in the polyclonal response elicited against SARS-CoV-2 spike is primarily limited to the quantity of antibodies generated rather than their domain specificity or relative neutralization potency. IMPORTANCE Infection by SARS-CoV-2 elicits antibodies against various domains of the spike protein, including the RBD, NTD and S2. Different infected individuals generate vastly different amounts of anti-spike antibodies. By contrast, as we show here, there is a remarkable similarity in the properties of the antibodies produced. Different individuals generate the same proportions of antibodies against each domain of the spike protein. Furthermore, the relationship between the amount of anti-spike antibodies produced and their neutralization efficacy of SARS-CoV-2 is highly conserved. Therefore, the observed variation in the neutralizing activity of the antibody response in COVID-19 convalescent subjects is caused by differences in the amounts of antibodies rather than their recognition properties or relative antiviral activity. These findings suggest that COVID-19 vaccine strategies that focus on enhancing the overall level of the antibodies will likely elicit a more uniformly efficacious protective response.

HIV-1 group M was transmitted to humans nearly one century ago. The virus has since diversified to form distinct clades, which spread to multiple regions worldwide. Of the different proteins encoded by HIV-1, the envelope glycoproteins (Envs) have diversified most rapidly in all infected populations. We compared the range of variants that emerged during the AIDS pandemic in diverse HIV-1 clades and distinct geographic regions. Our analyses focused on two components of Env that contain multiple epitopes of broadly-neutralizing antibodies: the glycan shield and apex domain. Interestingly, at each Env position, the amino acid in the inferred clade ancestor was replaced by a unique combination of emerging variants. Key antigenic sites and genetic signatures of vaccine protection have gradually evolved toward conserved frequency distributions (FDs) of all amino acids. FDs are specific for position and clade and are highly conserved in populations from different regions. Remarkably, founder effects of Env mutations in distinct clades and recently-infected regions were significantly reduced during the epidemic by evolution of each site toward the position-specific FD. These findings suggest that the selective pressures that guide evolution of Env are conserved in different populations. They are sufficiently strong to reduce founder effects at the clade and regional levels and have significantly altered the distribution of Env forms that circulate worldwide. Consequently, the intra-population diversity of the Env protein continues to increase whereas the inter-population diversity is gradually decreasing.

ABSTRACT Mutations in the spike glycoprotein of SARS-CoV-2 allow the virus to probe the sequence space in search of higher-fitness states. New sublineages of SARS-CoV-2 variants-of-concern (VOCs) continuously emerge with such mutations. Interestingly, the sites of mutation in these sublineages vary between the VOCs. Whether such differences reflect the random nature of mutation appearance or distinct evolutionary spaces of spike in the VOCs is unclear. Here we show that each position of spike has a lineage-specific likelihood for mutations to appear and dominate descendent sublineages. This likelihood can be accurately estimated from the lineage-specific mutational profile of spike at a protein-wide level. The mutability environment of each position, including adjacent sites on the protein structure and neighboring sites on the network of comutability, accurately forecast changes in descendent sublineages. Mapping of imminent changes within the VOCs can contribute to the design of immunogens and therapeutics that address future forms of SARS-CoV-2.