Tracking Changes in SARS-CoV-2 Spike: Evidence that D614G Increases Infectivity of the COVID-19 Virus.

A SARS-CoV-2 variant carrying the Spike protein amino acid change D614G has become the most prevalent form in the global pandemic. Dynamic tracking of variant frequencies revealed a recurrent pattern of G614 increase at multiple geographic levels: national, regional, and municipal. The shift occurred even in local epidemics where the original D614 form was well established prior to introduction of the G614 variant. The consistency of this pattern was highly statistically significant, suggesting that the G614 variant may have a fitness advantage. We found that the G614 variant grows to a higher titer as pseudotyped virions. In infected individuals, G614 is associated with lower RT-PCR cycle thresholds, suggestive of higher upper respiratory tract viral loads, but not with increased disease severity. These findings illuminate changes important for a mechanistic understanding of the virus and support continuing surveillance of Spike mutations to aid with development of immunological interventions.

High throughput analysis of B cell dynamics and neutralizing antibody development during immunization with a novel clade C HIV-1 envelope.

A protective HIV-1 vaccine has been hampered by a limited understanding of how B cells acquire neutralizing activity. Our previous vaccines expressing two different HIV-1 envelopes elicited robust antigen specific serum IgG titers in 20 rhesus macaques; yet serum from only two animals neutralized the autologous virus. Here, we used high throughput immunoglobulin receptor and single cell RNA sequencing to characterize the overall expansion, recall, and maturation of antigen specific B cells longitudinally over 90 weeks. Diversification and expansion of many B cell clonotypes occurred broadly in the absence of serum neutralization. However, in one animal that developed neutralization, two neutralizing B cell clonotypes arose from the same immunoglobulin germline and were tracked longitudinally. Early antibody variants with high identity to germline neutralized the autologous virus while later variants acquired somatic hypermutation and increased neutralization potency. The early engagement of precursors capable of neutralization with little to no SHM followed by prolonged affinity maturation allowed the two neutralizing lineages to successfully persist despite many other antigen specific B cells. The findings provide new insight into B cells responding to HIV-1 envelope during heterologous prime and boost immunization in rhesus macaques and the development of selected autologous neutralizing antibody lineages.

Acoustic force spectroscopy reveals subtle differences in cellulose unbinding behavior of carbohydrate-binding modules.

Protein adsorption to solid carbohydrate interfaces is critical to many biological processes, particularly in biomass deconstruction. To engineer more-efficient enzymes for biomass deconstruction into sugars, it is necessary to characterize the complex protein-carbohydrate interfacial interactions. A carbohydrate-binding module (CBM) is often associated with microbial surface-tethered cellulosomes or secreted cellulase enzymes to enhance substrate accessibility. However, it is not well known how CBMs recognize, bind, and dissociate from polysaccharides to facilitate efficient cellulolytic activity, due to the lack of mechanistic understanding and a suitable toolkit to study CBM-substrate interactions. Our work outlines a general approach to study the unbinding behavior of CBMs from polysaccharide surfaces using a highly multiplexed single-molecule force spectroscopy assay. Here, we apply acoustic force spectroscopy (AFS) to probe a Clostridium thermocellum cellulosomal scaffoldin protein (CBM3a) and measure its dissociation from nanocellulose surfaces at physiologically relevant, low force loading rates. An automated microfluidic setup and method for uniform deposition of insoluble polysaccharides on the AFS chip surfaces are demonstrated. The rupture forces of wild-type CBM3a, and its Y67A mutant, unbinding from nanocellulose surfaces suggests distinct multimodal CBM binding conformations, with structural mechanisms further explored using molecular dynamics simulations. Applying classical dynamic force spectroscopy theory, the single-molecule unbinding rate at zero force is extrapolated and found to agree with bulk equilibrium unbinding rates estimated independently using quartz crystal microbalance with dissipation monitoring. However, our results also highlight critical limitations of applying classical theory to explain the highly multivalent binding interactions for cellulose-CBM bond rupture forces exceeding 15 pN.

A neutralizing antibody target in early HIV-1 infection was recapitulated in rhesus macaques immunized with the transmitted/founder envelope sequence.

Transmitted/founder (T/F) HIV-1 envelope proteins (Envs) from infected individuals that developed neutralization breadth are likely to possess inherent features desirable for vaccine immunogen design. To explore this premise, we conducted an immunization study in rhesus macaques (RM) using T/F Env sequences from two human subjects, one of whom developed potent and broad neutralizing antibodies (Z1800M) while the other developed little to no neutralizing antibody responses (R66M) during HIV-1 infection. Using a DNA/MVA/protein immunization protocol, 10 RM were immunized with each T/F Env. Within each T/F Env group, the protein boosts were administered as either monomeric gp120 or stabilized trimeric gp140 protein. All vaccination regimens elicited high titers of antigen-specific IgG, and two animals that received monomeric Z1800M Env gp120 developed autologous neutralizing activity. Using early Env escape variants isolated from subject Z1800M as guides, the serum neutralizing activity of the two immunized RM was found to be dependent on the gp120 V5 region. Interestingly, the exact same residues of V5 were also targeted by a neutralizing monoclonal antibody (nmAb) isolated from the subject Z1800M early in infection. Glycan profiling and computational modeling of the Z1800M Env gp120 immunogen provided further evidence that the V5 loop is exposed in this T/F Env and was a dominant feature that drove neutralizing antibody targeting during infection and immunization. An expanded B cell clonotype was isolated from one of the neutralization-positive RM and nmAbs corresponding to this group demonstrated V5-dependent neutralization similar to both the RM serum and the human Z1800M nmAb. The results demonstrate that neutralizing antibody responses elicited by the Z1800M T/F Env in RM converged with those in the HIV-1 infected human subject, illustrating the potential of using immunogens based on this or other T/F Envs with well-defined immunogenicity as a starting point to drive breadth.

Visualization of the HIV-1 Env glycan shield across scales.

The dense array of N-linked glycans on the HIV-1 envelope glycoprotein (Env), known as the "glycan shield," is a key determinant of immunogenicity, yet intrinsic heterogeneity confounds typical structure-function analysis. Here, we present an integrated approach of single-particle electron cryomicroscopy (cryo-EM), computational modeling, and site-specific mass spectrometry (MS) to probe glycan shield structure and behavior at multiple levels. We found that dynamics lead to an extensive network of interglycan interactions that drive the formation of higher-order structure within the glycan shield. This structure defines diffuse boundaries between buried and exposed protein surface and creates a mapping of potentially immunogenic sites on Env. Analysis of Env expressed in different cell lines revealed how cryo-EM can detect subtle changes in glycan occupancy, composition, and dynamics that impact glycan shield structure and epitope accessibility. Importantly, this identified unforeseen changes in the glycan shield of Env obtained from expression in the same cell line used for vaccine production. Finally, by capturing the enzymatic deglycosylation of Env in a time-resolved manner, we found that highly connected glycan clusters are resistant to digestion and help stabilize the prefusion trimer, suggesting the glycan shield may function beyond immune evasion.

Exploring CRD mobility during RAS/RAF engagement at the membrane.

During the activation of mitogen-activated protein kinase (MAPK) signaling, the RAS-binding domain (RBD) and cysteine-rich domain (CRD) of RAF bind to active RAS at the plasma membrane. The orientation of RAS at the membrane may be critical for formation of the RAS-RBDCRD complex and subsequent signaling. To explore how RAS membrane orientation relates to the protein dynamics within the RAS-RBDCRD complex, we perform multiscale coarse-grained and all-atom molecular dynamics (MD) simulations of KRAS4b bound to the RBD and CRD domains of RAF-1, both in solution and anchored to a model plasma membrane. Solution MD simulations describe dynamic KRAS4b-CRD conformations, suggesting that the CRD has sufficient flexibility in this environment to substantially change its binding interface with KRAS4b. In contrast, when the ternary complex is anchored to the membrane, the mobility of the CRD relative to KRAS4b is restricted, resulting in fewer distinct KRAS4b-CRD conformations. These simulations implicate membrane orientations of the ternary complex that are consistent with NMR measurements. While a crystal structure-like conformation is observed in both solution and membrane simulations, a particular intermolecular rearrangement of the ternary complex is observed only when it is anchored to the membrane. This configuration emerges when the CRD hydrophobic loops are inserted into the membrane and helices α3-5 of KRAS4b are solvent exposed. This membrane-specific configuration is stabilized by KRAS4b-CRD contacts that are not observed in the crystal structure. These results suggest modulatory interplay between the CRD and plasma membrane that correlate with RAS/RAF complex structure and dynamics, and potentially influence subsequent steps in the activation of MAPK signaling.

Molecular origins of reduced activity and binding commitment of processive cellulases and associated carbohydrate-binding proteins to cellulose III.

Efficient enzymatic saccharification of cellulosic biomass into fermentable sugars can enable production of bioproducts like ethanol. Native crystalline cellulose, or cellulose I, is inefficiently processed via enzymatic hydrolysis but can be converted into the structurally distinct cellulose III allomorph that is processed via cellulase cocktails derived from Trichoderma reesei up to 20-fold faster. However, characterization of individual cellulases from T. reesei, like the processive exocellulase Cel7A, shows reduced binding and activity at low enzyme loadings toward cellulose III. To clarify this discrepancy, we monitored the single-molecule initial binding commitment and subsequent processive motility of Cel7A enzymes and associated carbohydrate-binding modules (CBMs) on cellulose using optical tweezers force spectroscopy. We confirmed a 48% lower initial binding commitment and 32% slower processive motility of Cel7A on cellulose III, which we hypothesized derives from reduced binding affinity of the Cel7A binding domain CBM1. Classical CBM-cellulose pull-down assays, depending on the adsorption model fitted, predicted between 1.2- and 7-fold reduction in CBM1 binding affinity for cellulose III. Force spectroscopy measurements of CBM1-cellulose interactions, along with molecular dynamics simulations, indicated that previous interpretations of classical binding assay results using multisite adsorption models may have complicated analysis, and instead suggest simpler single-site models should be used. These findings were corroborated by binding analysis of other type-A CBMs (CBM2a, CBM3a, CBM5, CBM10, and CBM64) on both cellulose allomorphs. Finally, we discuss how complementary analytical tools are critical to gain insight into the complex mechanisms of insoluble polysaccharides hydrolysis by cellulolytic enzymes and associated carbohydrate-binding proteins.

Binding and movement of individual Cel7A cellobiohydrolases on crystalline cellulose surfaces revealed by single-molecule fluorescence imaging.

The efficient catalytic conversion of biomass to bioenergy would meet a large portion of energy requirements in the near future. A crucial step in this process is the enzyme-catalyzed hydrolysis of cellulose to glucose that is then converted into fuel such as ethanol by fermentation. Here we use single-molecule fluorescence imaging to directly monitor the movement of individual Cel7A cellobiohydrolases from Trichoderma reesei (TrCel7A) on the surface of insoluble cellulose fibrils to elucidate molecular level details of cellulase activity. The motion of multiple, individual TrCel7A cellobiohydrolases was simultaneously recorded with ∼15-nm spatial resolution. Time-resolved localization microscopy provides insights on the activity of TrCel7A on cellulose and informs on nonproductive binding and diffusion. We measured single-molecule residency time distributions of TrCel7A bound to cellulose both in the presence of and absence of cellobiose the major product and a potent inhibitor of Cel7A activity. Combining these results with a kinetic model of TrCel7A binding provides microscopic insight into interactions between TrCel7A and the cellulose substrate.

Influence of membrane on the antigen presentation of the HIV-1 envelope membrane proximal external region (MPER).

The membrane proximal external region (MPER) of the HIV envelope glycoproteins has generated renewed interest after a recent phase I vaccine trial that presented MPER lipid-peptide epitopes demonstrated promise to elicit a broad neutralization response. The antigenicity of MPER is intimately associated with the membrane, and its presentation relies significantly on the lipid composition. This review brings together recent findings on the influence of membranes on the conformation of MPER and its recognition by broadly neutralizing antibodies. Specifically, the review highlights the importance of properly accounting for the balance between protein-protein and membrane-protein interactions in vaccine design.

Membrane lipids drive formation of KRAS4b-RAF1 RBDCRD nanoclusters on the membrane.

The oncogene RAS, extensively studied for decades, presents persistent gaps in understanding, hindering the development of effective therapeutic strategies due to a lack of precise details on how RAS initiates MAPK signaling with RAF effector proteins at the plasma membrane. Recent advances in X-ray crystallography, cryo-EM, and super-resolution fluorescence microscopy offer structural and spatial insights, yet the molecular mechanisms involving protein-protein and protein-lipid interactions in RAS-mediated signaling require further characterization. This study utilizes single-molecule experimental techniques, nuclear magnetic resonance spectroscopy, and the computational Machine-Learned Modeling Infrastructure (MuMMI) to examine KRAS4b and RAF1 on a biologically relevant lipid bilayer. MuMMI captures long-timescale events while preserving detailed atomic descriptions, providing testable models for experimental validation. Both in vitro and computational studies reveal that RBDCRD binding alters KRAS lateral diffusion on the lipid bilayer, increasing cluster size and decreasing diffusion. RAS and membrane binding cause hydrophobic residues in the CRD region to penetrate the bilayer, stabilizing complexes through ß-strand elongation. These cooperative interactions among lipids, KRAS4b, and RAF1 are proposed as essential for forming nanoclusters, potentially a critical step in MAP kinase signal activation.

Dynamic electrophoretic fingerprinting of the HIV-1 envelope glycoprotein.

BACKGROUND: Interactions between the HIV-1 envelope glycoprotein (Env) and its primary receptor CD4 are influenced by the physiological setting in which these events take place. In this study, we explored the surface chemistry of HIV-1 Env constructs at a range of pH and salinities relevant to mucosal and systemic compartments through electrophoretic mobility (EM) measurements. Sexual transmission events provide a more acidic environment for HIV-1 compared to dissemination and spread of infection occurring in blood or lymph node. We hypothesize functional, trimeric Env behaves differently than monomeric forms. RESULTS: The dynamic electrophoretic fingerprint of trimeric gp140 revealed a change in EM from strongly negative to strongly positive as pH increased from that of the lower female genital tract (pHx) to that of the blood (pHy). Similar findings were observed using a trimeric influenza Haemagglutinin (HA) glycoprotein, indicating that this may be a general attribute of trimeric viral envelope glycoproteins. These findings were supported by computationally modeling the surface charge of various gp120 and HA crystal structures. To identify the behavior of the infectious agent and its target cells, EM measurements were made on purified whole HIV-1 virions and primary T-lymphocytes. Viral particles had a largely negative surface charge, and lacked the regions of positivity near neutral pH that were observed with trimeric Env. T cells changed their surface chemistry as a function of activation state, becoming more negative over a wider range of pH after activation. Soluble recombinant CD4 (sCD4) was found to be positively charged under a wide range of conditions. Binding studies between sCD4 and gp140 show that the affinity of CD4-gp140 interactions depends on pH. CONCLUSIONS: Taken together, these findings allow a more complete model of the electrochemical forces involved in HIV-1 Env functionality. These results indicate that the influence of the localized environment on the interactions of HIV with target cells are more pronounced than previously appreciated. There is differential chemistry of trimeric, but not monomeric, Env under conditions which mimic the mucosa compared to those found systemically. This should be taken into consideration during design of immunogens which targets virus at mucosal portals of entry.

Molecular insights into antibody-mediated protection against the prototypic simian immunodeficiency virus.

SIVmac239 infection of macaques is a favored model of human HIV infection. However, the SIVmac239 envelope (Env) trimer structure, glycan occupancy, and the targets and ability of neutralizing antibodies (nAbs) to protect against SIVmac239 remain unknown. Here, we report the isolation of SIVmac239 nAbs that recognize a glycan hole and the V1/V4 loop. A high-resolution structure of a SIVmac239 Env trimer-nAb complex shows many similarities to HIV and SIVcpz Envs, but with distinct V4 features and an extended V1 loop. Moreover, SIVmac239 Env has a higher glycan shield density than HIV Env that may contribute to poor or delayed nAb responses in SIVmac239-infected macaques. Passive transfer of a nAb protects macaques from repeated intravenous SIVmac239 challenge at serum titers comparable to those described for protection of humans against HIV infection. Our results provide structural insights for vaccine design and shed light on antibody-mediated protection in the SIV model.

Asynchronous Reciprocal Coupling of Martini 2.2 Coarse-Grained and CHARMM36 All-Atom Simulations in an Automated Multiscale Framework.

The appeal of multiscale modeling approaches is predicated on the promise of combinatorial synergy. However, this promise can only be realized when distinct scales are combined with reciprocal consistency. Here, we consider multiscale molecular dynamics (MD) simulations that combine the accuracy and macromolecular flexibility accessible to fixed-charge all-atom (AA) representations with the sampling speed accessible to reductive, coarse-grained (CG) representations. AA-to-CG conversions are relatively straightforward because deterministic routines with unique outcomes are achievable. Conversely, CG-to-AA conversions have many solutions due to a surge in the number of degrees of freedom. While automated tools for biomolecular CG-to-AA transformation exist, we find that one popular option, called Backward, is prone to stochastic failure and the AA models that it does generate frequently have compromised protein structure and incorrect stereochemistry. Although these shortcomings can likely be circumvented by human intervention in isolated instances, automated multiscale coupling requires reliable and robust scale conversion. Here, we detail an extension to Multiscale Machine-learned Modeling Infrastructure (MuMMI), including an improved CG-to-AA conversion tool called sinceCG. This tool is reliable (∼98% weakly correlated repeat success rate), automatable (no unrecoverable hangs), and yields AA models that generally preserve protein secondary structure and maintain correct stereochemistry. We describe how the MuMMI framework identifies CG system configurations of interest, converts them to AA representations, and simulates them at the AA scale while on-the-fly analyses provide feedback to update CG parameters. Application to systems containing the peripheral membrane protein RAS and proximal components of RAF kinase on complex eight-component lipid bilayers with ∼1.5 million atoms is discussed in the context of MuMMI.

Exploring the Role of Glycans in the Interaction of SARS-CoV-2 RBD and Human Receptor ACE2.

COVID-19 is a highly infectious respiratory disease caused by the novel coronavirus SARS-CoV-2. It has become a global pandemic and its frequent mutations may pose new challenges for vaccine design. During viral infection, the Spike RBD of SARS-CoV-2 binds the human host cell receptor ACE2, enabling the virus to enter the host cell. Both the Spike and ACE2 are densely glycosylated, and it is unclear how distinctive glycan types may modulate the interaction of RBD and ACE2. Detailed understanding of these determinants is key for the development of novel therapeutic strategies. To this end, we perform extensive all-atom simulations of the (i) RBD-ACE2 complex without glycans, (ii) RBD-ACE2 with oligomannose MAN9 glycans in ACE2, and (iii) RBD-ACE2 with complex FA2 glycans in ACE2. These simulations identify the key residues at the RBD-ACE2 interface that form contacts with higher probabilities, thus providing a quantitative evaluation that complements recent structural studies. Notably, we find that this RBD-ACE2 contact signature is not altered by the presence of different glycoforms, suggesting that RBD-ACE2 interaction is robust. Applying our simulated results, we illustrate how the recently prevalent N501Y mutation may alter specific interactions with host ACE2 that facilitate the virus-host binding. Furthermore, our simulations reveal how the glycan on Asn90 of ACE2 can play a distinct role in the binding and unbinding of RBD. Finally, an energetics analysis shows that MAN9 glycans on ACE2 decrease RBD-ACE2 affinity, while FA2 glycans lead to enhanced binding of the complex. Together, our results provide a more comprehensive picture of the detailed interplay between virus and human receptor, which is much needed for the discovery of effective treatments that aim at modulating the physical-chemical properties of this virus.

SARS-CoV-2 variant B.1.1.7 is susceptible to neutralizing antibodies elicited by ancestral spike vaccines.

All current vaccines for COVID-19 utilize ancestral SARS-CoV-2 spike with the goal of generating protective neutralizing antibodies. The recent emergence and rapid spread of several SARS-CoV-2 variants carrying multiple spike mutations raise concerns about possible immune escape. One variant, first identified in the United Kingdom (B.1.1.7, also called 20I/501Y.V1), contains eight spike mutations with potential to impact antibody therapy, vaccine efficacy, and risk of reinfection. Here, we show that B.1.1.7 remains sensitive to neutralization, albeit at moderately reduced levels (∼sim;2-fold), by serum samples from convalescent individuals and recipients of an mRNA vaccine (mRNA-1273, Moderna) and a protein nanoparticle vaccine (NVX-CoV2373, Novavax). A subset of monoclonal antibodies to the receptor binding domain (RBD) of spike are less effective against the variant, while others are largely unaffected. These findings indicate that variant B.1.1.7 is unlikely to be a major concern for current vaccines or for an increased risk of reinfection.

HIV-1 and SARS-CoV-2: Patterns in the evolution of two pandemic pathogens.

Humanity is currently facing the challenge of two devastating pandemics caused by two very different RNA viruses: HIV-1, which has been with us for decades, and SARS-CoV-2, which has swept the world in the course of a single year. The same evolutionary strategies that drive HIV-1 evolution are at play in SARS-CoV-2. Single nucleotide mutations, multi-base insertions and deletions, recombination, and variation in surface glycans all generate the variability that, guided by natural selection, enables both HIV-1's extraordinary diversity and SARS-CoV-2's slower pace of mutation accumulation. Even though SARS-CoV-2 diversity is more limited, recently emergent SARS-CoV-2 variants carry Spike mutations that have important phenotypic consequences in terms of both antibody resistance and enhanced infectivity. We review and compare how these mutational patterns manifest in these two distinct viruses to provide the variability that fuels their evolution by natural selection.

SARS-CoV-2 variant B.1.1.7 is susceptible to neutralizing antibodies elicited by ancestral Spike vaccines.

The SARS-CoV-2 Spike glycoprotein mediates virus entry and is a major target for neutralizing antibodies. All current vaccines are based on the ancestral Spike with the goal of generating a protective neutralizing antibody response. Several novel SARS-CoV-2 variants with multiple Spike mutations have emerged, and their rapid spread and potential for immune escape have raised concerns. One of these variants, first identified in the United Kingdom, B.1.1.7 (also called VUI202012/01), contains eight Spike mutations with potential to impact antibody therapy, vaccine efficacy and risk of reinfection. Here we employed a lentivirus-based pseudovirus assay to show that variant B.1.1.7 remains sensitive to neutralization, albeit at moderately reduced levels (~2-fold), by serum samples from convalescent individuals and recipients of two different vaccines based on ancestral Spike: mRNA-1273 (Moderna), and protein nanoparticle NVX-CoV2373 (Novavax). Some monoclonal antibodies to the receptor binding domain (RBD) of Spike were less effective against the variant while others were largely unaffected. These findings indicate that B.1.1.7 is not a neutralization escape variant that would be a major concern for current vaccines, or for an increased risk of reinfection.

Improving pandemic influenza risk assessment.

Assessing the pandemic risk posed by specific non-human influenza A viruses is an important goal in public health research. As influenza virus genome sequencing becomes cheaper, faster, and more readily available, the ability to predict pandemic potential from sequence data could transform pandemic influenza risk assessment capabilities. However, the complexities of the relationships between virus genotype and phenotype make such predictions extremely difficult. The integration of experimental work, computational tool development, and analysis of evolutionary pathways, together with refinements to influenza surveillance, has the potential to transform our ability to assess the risks posed to humans by non-human influenza viruses and lead to improved pandemic preparedness and response.