Network analysis uncovers the communication structure of SARS-CoV-2 spike protein identifying sites for immunogen design.

The COVID-19 pandemic, caused by the SARS-CoV-2 virus, has triggered myriad efforts to understand the structure and dynamics of this complex pathogen. The spike glycoprotein of SARS-CoV-2 is a significant target for immunogens as it is the means by which the virus enters human cells, while simultaneously sporting mutations responsible for immune escape. These functional and escape processes are regulated by complex molecular-level interactions. Our study presents quantitative insights on domain and residue contributions to allosteric communication, immune evasion, and local- and global-level control of functions through the derivation of a weighted graph representation from all-atom MD simulations. Focusing on the ancestral form and the D614G-variant, we provide evidence of the utility of our approach by guiding the selection of a mutation that alters the spike's stability. Taken together, the network approach serves as a valuable tool to evaluate communication "hot-spots" in proteins to guide design of stable immunogens.

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.

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.

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.

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.

The SARS-CoV-2 Spike variant D614G favors an open conformational state.

The COVID-19 (coronavirus disease 2019) pandemic underwent a rapid transition with the emergence of a dominant viral variant (from the "D-form" to the "G-form") that carried an amino acid substitution D614G in its "Spike" protein. The G-form is more infectious in vitro and is associated with increased viral loads in the upper airways. To gain insight into the molecular-level underpinnings of these characteristics, we used microsecond all-atom simulations. We show that changes in the protein energetics favor a higher population of infection-capable states in the G-form through release of asymmetry present in the D-form inter-protomer interactions. Thus, the increased infectivity of the G-form is likely due to a higher rate of profitable binding encounters with the host receptor. It is also predicted to be more neutralization sensitive owing to enhanced exposure of the receptor binding domain, a key target region for neutralizing antibodies. These results are critical for vaccine design.

Development of Martini 2.2 parameters for N-glycans: a case study of the HIV-1 Env glycoprotein dynamics.

N-linked glycans are ubiquitous in nature and play key roles in biology. For example, glycosylation of pathogenic proteins is a common immune evasive mechanism, hampering the development of successful vaccines. Due to their chemical variability and complex dynamics, an accurate molecular understanding of glycans is still limited by the lack of effective resolution of current experimental approaches. Here, we have developed and implemented a reductive model based on the popular Martini 2.2 coarse-grained force field for the computational study of N-glycosylation. We used the HIV-1 Env as a direct applied example of a highly glycosylated protein. Our results indicate that the model not only reproduces many observables in very good agreement with a fully atomistic force field but also can be extended to study large amount of glycosylation variants, a fundamental property that can aid in the development of drugs and vaccines.

Quantification of the Resilience and Vulnerability of HIV-1 Native Glycan Shield at Atomistic Detail.

Dense surface glycosylation on the HIV-1 envelope (Env) protein acts as a shield from the adaptive immune system. However, the molecular complexity and flexibility of glycans make experimental studies a challenge. Here we have integrated high-throughput atomistic modeling of fully glycosylated HIV-1 Env with graph theory to capture immunologically important features of the shield topology. This is the first complete all-atom model of HIV-1 Env SOSIP glycan shield that includes both oligomannose and complex glycans, providing physiologically relevant insights of the glycan shield. This integrated approach including quantitative comparison with cryo-electron microscopy data provides hitherto unexplored details of the native shield architecture and its difference from the high-mannose glycoform. We have also derived a measure to quantify the shielding effect over the antigenic protein surface that defines regions of relative vulnerability and resilience of the shield and can be harnessed for rational immunogen design.

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.

The SARS-CoV-2 Spike Variant D614G Favors an Open Conformational State.

The COVID-19 pandemic underwent a rapid transition with the emergence of a SARS-CoV-2 variant that carried the amino acid substitution D614G in the Spike protein that became globally prevalent. The G-form is both more infectious in vitro and associated with increased viral loads in infected people. To gain insight into the mechanism underlying these distinctive characteristics, we employed multiple replicas of microsecond all-atom simulations to probe the molecular-level impact of this substitution on Spike closed and open states. The open state enables Spike interactions with its human cellular receptor, ACE2. Here we show that changes in the inter-protomer energetics due to the D614G substitution favor a higher population of infection-capable (open) states. The inter-protomer interactions between S1 and S2 subunits in the open state of the D-form are asymmetric. This asymmetry is resolved in the G-form due to the release of tensile hydrogen bonds resulting in an increased population of open conformations. Thus, the increased infectivity of the G-form is likely due to a higher rate of profitable binding encounters with the host receptor. It is also predicted to be more neutralization sensitive due to enhanced exposure of the receptor binding domain, a key target region for neutralizing antibodies.

Graph-Directed Approach for Downselecting Toxins for Experimental Structure Determination.

Conotoxins are short, cysteine-rich peptides of great interest as novel therapeutic leads and of great concern as lethal biological agents due to their high affinity and specificity for various receptors involved in neuromuscular transmission. Currently, of the approximately 6000 known conotoxin sequences, only about 3% have associated structural characterization, which leads to a bottleneck in rapid high-throughput screening (HTS) for identification of potential leads or threats. In this work, we combine a graph-based approach with homology modeling to expand the library of conotoxin structures and to identify those conotoxin sequences that are of the greatest value for experimental structural characterization. The latter would allow for the rapid expansion of the known structural space for generating high quality template-based models. Our approach generalizes to other evolutionarily-related, short, cysteine-rich venoms of interest. Overall, we present and validate an approach for venom structure modeling and experimental guidance and employ it to produce a 290%-larger library of approximate conotoxin structures for HTS. We also provide a set of ranked conotoxin sequences for experimental structure determination to further expand this library.

Emergence of Alternative Structures in Amyloid Beta 1-42 Monomeric Landscape by N-terminal Hexapeptide Amyloid Inhibitors.

Alzheimer's disease (AD) is characterized by deposition of amyloid beta (Aß) peptides into senile plaques in the brain. While most familial mutations are associated with early-onset AD, recent studies report the AD-protective nature of two genetic human Aß variants, i.e. A2T and A2V, in the heterozygous state. The mixture of A2V Aß1-6 (Aß6) hexapeptide and WT Aß1-42 (Αß42) is also found neuroprotective. Motivated by these findings, in this study we investigate the effects of WT, A2V, and A2T Aß6 hexapeptide binding on the monomeric WT Aß42 landscape. For this purpose, we have performed extensive atomistic Replica Exchange Molecular Dynamics simulations, elucidating preferential binding of Aß42 with the A2V and A2T hexapeptides compared to WT Aß6. A notable reorganization of the Aß42 landscape is revealed due to hexapeptide association, as manifested by lowering of transient interactions between the central and C-terminal hydrophobic patches. Concurrently, Aß6-bound Aß42 monomer exhibits alternative structural features that are strongly dependent on the hexapeptide sequence. For example, a central helix is more frequently populated within the A2T-bound monomer, while A2V-bound Aß42 is often enhanced in overall disorder. Taken together, the present simulations offer novel molecular insights onto the effect of the N-terminal hexapeptide binding on the Aß42 monomer structure, which might help in explaining their reported amyloid inhibition properties.

A mechanism for acetylcholine receptor gating based on structure, coupling, phi, and flip.

Nicotinic acetylcholine receptors are allosteric proteins that generate membrane currents by isomerizing ("gating") between resting and active conformations under the influence of neurotransmitters. Here, to explore the mechanisms that link the transmitter-binding sites (TBSs) with the distant gate, we use mutant cycle analyses to measure coupling between residue pairs, phi value analyses to sequence domain rearrangements, and current simulations to reproduce a microsecond shut component ("flip") apparent in single-channel recordings. Significant interactions between amino acids separated by >15 Å are rare; an exception is between the αM2-M3 linkers and the TBSs that are ∼30 Å apart. Linker residues also make significant, local interactions within and between subunits. Phi value analyses indicate that without agonists, the linker is the first region in the protein to reach the gating transition state. Together, the phi pattern and flip component suggest that a complete, resting↔active allosteric transition involves passage through four brief intermediate states, with brief shut events arising from sojourns in all or a subset. We derive energy landscapes for gating with and without agonists, and propose a structure-based model in which resting→active starts with spontaneous rearrangements of the M2-M3 linkers and TBSs. These conformational changes stabilize a twisted extracellular domain to promote transmembrane helix tilting, gate dilation, and the formation of a "bubble" that collapses to initiate ion conduction. The energy landscapes suggest that twisting is the most energetically unfavorable step in the resting→active conformational change and that the rate-limiting step in the reverse process is bubble formation.

Structural correlates of affinity in fetal versus adult endplate nicotinic receptors.

Adult-type nicotinic acetylcholine receptors (AChRs) mediate signalling at mature neuromuscular junctions and fetal-type AChRs are necessary for proper synapse development. Each AChR has two neurotransmitter binding sites located at the interface of a principal and a complementary subunit. Although all agonist binding sites have the same core of five aromatic amino acids, the fetal site has ∼30-fold higher affinity for the neurotransmitter ACh. Here we use molecular dynamics simulations of adult versus fetal homology models to identify complementary-subunit residues near the core that influence affinity, and use single-channel electrophysiology to corroborate the results. Four residues in combination determine adult versus fetal affinity. Simulations suggest that at lower-affinity sites, one of these unsettles the core directly and the others (in loop E) increase backbone flexibility to unlock a key, complementary tryptophan from the core. Swapping only four amino acids is necessary and sufficient to exchange function between adult and fetal AChRs.

Function of the M1 π-helix in endplate receptor activation and desensitization.

KEY POINTS: A conserved proline in M1 causes a kink between α and π helical segments. The kink is under greater tension in the resting versus active conformation. The kink and the agonist do not interact directly. The π-helix separates the gating functions of the extracellular and transmembrane domains. Mutations of the conserved proline and propofol increase desensitization. ABSTRACT: Nicotinic acetylcholine receptors (AChRs) switch on/off to generate transient membrane currents (C↔O; closed-open 'gating') and enter/recover from long-lived, refractory states (O↔D; 'desensitization'). The M1 transmembrane helix of the muscle endplate AChR is linked to a ß-strand of the extracellular domain that extends to a neurotransmitter binding site. We used electrophysiology to measure the effects of mutations of amino acids that are located at a proline kink in M1 that separates π and α helices, in both α (N217, V218 and P221) and non-α subunits. In related receptors, the kink is straighter and more stable in O vs. C structures (gating is 'spring-loaded'). None of the AChR kink mutations had a measureable effect on agonist affinity but many influenced the allosteric gating constant substantially. Side chains in the M1 α-helix experience extraordinarily large energy differences between C and O structures, probably because of a ∼2 Å displacement and tilt of M2 relative to M1. There is a discrete break in the character of the gating transition state between αN217 and αV218, indicating that the π-helix is a border between extracellular- and transmembrane-domain function. Mutations of the conserved M1 proline, and the anaesthetic propofol, increase a rate constant for desensitization. The results suggest that straightening of the M1 proline kink triggers AChR desensitization.

Decrypting the structural, dynamic, and energetic basis of a monomeric kinesin interacting with a tubulin dimer in three ATPase states by all-atom molecular dynamics simulation.

We have employed molecular dynamics (MD) simulation to investigate, with atomic details, the structural dynamics and energetics of three major ATPase states (ADP, APO, and ATP state) of a human kinesin-1 monomer in complex with a tubulin dimer. Starting from a recently solved crystal structure of ATP-like kinesin-tubulin complex by the Knossow lab, we have used flexible fitting of cryo-electron-microscopy maps to construct new structural models of the kinesin-tubulin complex in APO and ATP state, and then conducted extensive MD simulations (total 400 ns for each state), followed by flexibility analysis, principal component analysis, hydrogen bond analysis, and binding free energy analysis. Our modeling and simulation have revealed key nucleotide-dependent changes in the structure and flexibility of the nucleotide-binding pocket (featuring a highly flexible and open switch I in APO state) and the tubulin-binding site, and allosterically coupled motions driving the APO to ATP transition. In addition, our binding free energy analysis has identified a set of key residues involved in kinesin-tubulin binding. On the basis of our simulation, we have attempted to address several outstanding issues in kinesin study, including the possible roles of ß-sheet twist and neck linker docking in regulating nucleotide release and binding, the structural mechanism of ADP release, and possible extension and shortening of α4 helix during the ATPase cycle. This study has provided a comprehensive structural and dynamic picture of kinesin's major ATPase states, and offered promising targets for future mutational and functional studies to investigate the molecular mechanism of kinesin motors.

Functional differences between neurotransmitter binding sites of muscle acetylcholine receptors.

A muscle acetylcholine receptor (AChR) has two neurotransmitter binding sites located in the extracellular domain, at αδ and either αε (adult) or αγ (fetal) subunit interfaces. We used single-channel electrophysiology to measure the effects of mutations of five conserved aromatic residues at each site with regard to their contribution to the difference in free energy of agonist binding to active versus resting receptors (ΔGB1). The two binding sites behave independently in both adult and fetal AChRs. For four different agonists, including ACh and choline, ΔGB1 is ∼-2 kcal/mol more favorable at αγ compared with at αε and αδ. Only three of the aromatics contribute significantly to ΔGB1 at the adult sites (αY190, αY198, and αW149), but all five do so at αγ (as well as αY93 and γW55). γW55 makes a particularly large contribution only at αγ that is coupled energetically to those contributions of some of the α-subunit aromatics. The hydroxyl and benzene groups of loop C residues αY190 and αY198 behave similarly with regard to ΔGB1 at all three kinds of site. ACh binding energies estimated from molecular dynamics simulations are consistent with experimental values from electrophysiology and suggest that the αγ site is more compact, better organized, and less dynamic than αε and αδ. We speculate that the different sensitivities of the fetal αγ site versus the adult αε and αδ sites to choline and ACh are important for the proper maturation and function of the neuromuscular synapse.