Microsecond Dynamics and Network Analysis of the HIV-1 SOSIP Env Trimer Reveal Collective Behavior and Conserved Microdomains of the Glycan Shield

Vaccine-elicited and naturally elicited antibodies differ in their recognition of the HIV-1 fusion peptide

Broadly neutralizing antibodies have been proposed as templates for HIV-1 vaccine design, but it has been unclear how similar vaccine-elicited antibodies are to their naturally elicited templates. To provide insight, here we compare the recognition of naturally elicited and vaccine-elicited antibodies targeting the HIV-1 fusion peptide, which comprises envelope (Env) residues 512-526, with the most common sequence being AVGIGAVFLGFLGAA. Naturally elicited antibodies bound peptides with substitutions to negatively charged amino acids at residue positions 517-520 substantially better than the most common sequence, despite these substitutions rarely appearing in HIV-1; by contrast, vaccine-elicited antibodies were less tolerant of sequence variation, with no substitution of residues 512-516 showing increased binding. Molecular dynamics analysis and cryo-EM structural analysis of the naturally elicited ACS202 antibody in complex with the HIV-1 Env trimer with an alanine 517 to glutamine substitution suggested enhanced binding to result from electrostatic interactions with positively charged antibody residues. Overall, vaccine-elicited antibodies appeared to be more fully optimized to bind the most common fusion peptide sequence, perhaps reflecting the immunization with fusion peptide of the vaccine-elicited antibodies.

Glycosylation and the global virome

The sugars that coat the outsides of viruses and host cells are key to successful disease transmission, but they remain understudied compared to other molecular features. Understanding the comparative zoology of glycosylation - and harnessing it for predictive science - could help close the molecular gap in zoonotic risk assessment.

Principles of SARS-CoV-2 Glycosylation

The structure and post-translational processing of the SARS-CoV-2 spike glycoprotein (S) is intimately associated with the function of the virus and of sterilising vaccines. The surface of the S protein is extensively modified by glycans, and their biosynthesis is driven by both the wider cellular context, and importantly, the underlining protein structure and local glycan density. Comparison of virally derived S protein with both recombinantly derived and adenovirally induced proteins, reveal hotspots of protein-directed glycosylation that drive conserved glycosylation motifs. Molecular dynamics simulations revealed that, while the S surface is extensively shielded by N-glycans, it presents regions vulnerable to neutralising antibodies. Furthermore, glycans have been shown to influence the accessibility of the receptor binding domain and the binding to the cellular receptor. The emerging picture is one of unifying, principles of S protein glycosylation and an intimate role of glycosylation in immunogen structure and efficacy.

Conformational plasticity of the HIV-1 gp41 immunodominant region is recognized by multiple non-neutralizing antibodies

The early humoral immune response to acute HIV-1 infection is largely non-neutralizing. The principal target of these antibodies is the primary immunodominant region (PID) on the gp41 fusion protein. The PID is a highly conserved 15-residue region displayed on the surface of HIV-1 virions. In this study, we analyzed the humoral determinants of HIV-1 gp41 PID binding using biophysical, structural, and computational methods. In complex with a patient-derived near-germline antibody fragment, the PID motif adopts an elongated random coil, whereas the PID bound to affinity-matured Fab adopts a strand-turn-helix conformation. Molecular dynamics simulations showed that the PID is structurally plastic suggesting that the PID can form an ensemble of structural states recognized by various non-neutralizing antibodies, facilitating HIV-1 immunodominance observed in acute and chronic HIV-1 infections. An improved understanding of how the HIV-1 gp41 PID misdirects the early humoral response should guide the development of an effective HIV-1 vaccine.

The 15-amino-acid primary immunodominant (PID) region on HIV-1 gp41 adopts an ensemble of conformational states. This conformational plasticity is suggested to misdirect the early humoral immune response.

GLYCO: a tool to quantify glycan shielding of glycosylated proteins

MOTIVATION

Glycans play important roles in protein folding and cell-cell interactions-and, furthermore, glycosylation of protein antigens can dramatically impact immune responses. While there have been attempts to quantify the glycan shielding or coverage of a protein surface, none of the publicly available tools analyzes glycan shielding computationally at an atomistic level.

RESULTS

Here, we developed an in silico approach, GLYCO (GLYcan COverage), to quantify the glycan shielding of a protein surface. The software provides insights into glycan-dense/sparse regions of the entire protein surface or a subset of the protein surface. GLYCO calculates glycan shielding from a single coordinate file or from multiple coordinate files, for instance, as obtained from molecular dynamics simulations or by nuclear magnetic resonance spectroscopy structure determination, enabling analysis of glycan dynamics. Overall, GLYCO provides fundamental insights into the glycan shielding of glycosylated proteins.

AVAILABILITY AND IMPLEMENTATION

GLYCO is freely available at GitHub (https://github.com/myungjinlee/GLYCO).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Unsupervised Learning Methods for Molecular Simulation Data

Unsupervised learning is becoming an essential tool to analyze the increasingly large amounts of data produced by atomistic and molecular simulations, in material science, solid state physics, biophysics, and biochemistry. In this Review, we provide a comprehensive overview of the methods of unsupervised learning that have been most commonly used to investigate simulation data and indicate likely directions for further developments in the field. In particular, we discuss feature representation of molecular systems and present state-of-the-art algorithms of dimensionality reduction, density estimation, and clustering, and kinetic models. We divide our discussion into self-contained sections, each discussing a specific method. In each section, we briefly touch upon the mathematical and algorithmic foundations of the method, highlight its strengths and limitations, and describe the specific ways in which it has been used-or can be used-to analyze molecular simulation data.

Hitting the sweet spot: exploiting HIV-1 glycan shield for induction of broadly neutralizing antibodies

PURPOSE OF REVIEW

The surface of the HIV-1 Env glycoprotein, the target of neutralizing antibodies, is extensively covered by N-linked glycans that create a glycan shield. Broadly neutralizing antibodies (bNAbs), the primary targets of HIV-1 vaccine design, have to negotiate this glycan shield. Here, we review the barriers and opportunities that the HIV-1 glycan shield presents for vaccine induction of bNAbs.

RECENT FINDINGS

Glycan shields can impact the nature of the antibody response and influence the development of neutralization breadth in HIV-1 infections. The architecture of the glycan shield arising from glycan interactions and dynamics have been modeled, and its fine structure, that is, the site-wise glycan heterogeneity, has been determined for some isolates. Although the extent of glycan shielding is conserved, the precise number, location and processing of glycans, however, is strain-dependent. New insights continue to reveal how such differences can impact bNAb activity and development. Novel approaches have exploited the glycan shield for designing immunogens that bind the germline precursors of bNAbs, a critical roadblock for vaccine-induction of bNAbs.

SUMMARY

The HIV-1 glycan shield can significantly impact the induction and maturation of bNAbs, and a better understanding of how to manipulate it will improve immunogen design.

Protein and glycan molecular weight determination of highly glycosylated HIV-1 envelope trimers by HPSEC-MALS

High Performance Size-Exclusion Chromatography coupled with Multi-Angle Light Scattering detection (HPSEC-MALS) is an important tool to provide a reliable molecular weight measurement for a large complex biomolecule. A recent HIV-1 soluble envelope trimer vaccine candidate, BG505 DS-SOSIP.664, is among the most glycosylated proteins to enter a clinical trial to date, and determination of its protein and glycan molecular weight is one of the key attributes in pre-clinical characterization. However, protein and glycans possess disparate dndcvalues making molecular weight measurement inaccurate in conventional SEC-MALS. To overcome these challenges, a simple mathematically guided experiment was explored, and a composite dndcvalue was established by utilizing protein and glycan mass contributions for the HIV-1 envelope trimer. This establishment was further verified by an orthogonal mass spectrometry analysis. This innovative, simple, and quick analytical approach can be applied broadly to measuring the molecular weight of various composite molecular structures, such as complex glycoconjugates.

Glycans in Virus-Host Interactions: A Structural Perspective

Many interactions between microbes and their hosts are driven or influenced by glycans, whose heterogeneous and difficult to characterize structures have led to an underappreciation of their role in these interactions compared to protein-based interactions. Glycans decorate microbe glycoproteins to enhance attachment and fusion to host cells, provide stability, and evade the host immune system. Yet, the host immune system may also target these glycans as glycoepitopes. In this review, we provide a structural perspective on the role of glycans in host-microbe interactions, focusing primarily on viral glycoproteins and their interactions with host adaptive immunity. In particular, we discuss a class of topological glycoepitopes and their interactions with topological mAbs, using the anti-HIV mAb 2G12 as the archetypical example. We further offer our view that structure-based glycan targeting strategies are ready for application to viruses beyond HIV, and present our perspective on future development in this area.

HIV-1 Envelope Conformation, Allostery, and Dynamics

The HIV-1 envelope glycoprotein (Env) mediates host cell fusion and is the primary target for HIV-1 vaccine design. The Env undergoes a series of functionally important conformational rearrangements upon engagement of its host cell receptor, CD4. As the sole target for broadly neutralizing antibodies, our understanding of these transitions plays a critical role in vaccine immunogen design. Here, we review available experimental data interrogating the HIV-1 Env conformation and detail computational efforts aimed at delineating the series of conformational changes connecting these rearrangements. These studies have provided a structural mapping of prefusion closed, open, and transition intermediate structures, the allosteric elements controlling rearrangements, and state-to-state transition dynamics. The combination of these investigations and innovations in molecular modeling set the stage for advanced studies examining rearrangements at greater spatial and temporal resolution.

Exploring dynamics and network analysis of spike glycoprotein of SARS-COV-2

The ongoing pandemic caused by severe acute respiratory syndrome coronavirus 2 continues to rage with devastating consequences on human health and global economy. The spike glycoprotein on the surface of coronavirus mediates its entry into host cells and is the target of all current antibody design efforts to neutralize the virus. The glycan shield of the spike helps the virus to evade the human immune response by providing a thick sugar-coated barrier against any antibody. To study the dynamic motion of glycans in the spike protein, we performed microsecond-long molecular dynamics simulation in two different states that correspond to the receptor binding domain in open or closed conformations. Analysis of this microsecond-long simulation revealed a scissoring motion on the N-terminal domain of neighboring monomers in the spike trimer. The roles of multiple glycans in shielding of spike protein in different regions were uncovered by a network analysis, in which the high betweenness centrality of glycans at the apex revealed their importance and function in the glycan shield. Microdomains of glycans were identified featuring a high degree of intracommunication in these microdomains. An antibody overlap analysis revealed the glycan microdomains as well as individual glycans that inhibit access to the antibody epitopes on the spike protein. Overall, the results of this study provide detailed understanding of the spike glycan shield, which may be utilized for therapeutic efforts against this crisis.

HIV-1 Envelope Glycosylation and the Signal Peptide

The RV144 trial represents the only vaccine trial to demonstrate any protective effect against HIV-1 infection. While the reason(s) for this protection are still being evaluated, it serves as justification for widespread efforts aimed at developing new, more effective HIV-1 vaccines. Advances in our knowledge of HIV-1 immunogens and host antibody responses to these immunogens are crucial to informing vaccine design. While the envelope (Env) protein is the only viral protein present on the surface of virions, it exists in a complex trimeric conformation and is decorated with an array of variable N-linked glycans, making it an important but difficult target for vaccine design. Thus far, efforts to elicit a protective humoral immune response using structural mimics of native Env trimers have been unsuccessful. Notably, the aforementioned N-linked glycans serve as a component of many of the epitopes crucial for the induction of potentially protective broadly neutralizing antibodies (bnAbs). Thus, a greater understanding of Env structural determinants, most critically Env glycosylation, will no doubt be of importance in generating effective immunogens. Recent studies have identified the HIV-1 Env signal peptide (SP) as an important contributor to Env glycosylation. Further investigation into the mechanisms by which the SP directs glycosylation will be important, both in the context of understanding HIV-1 biology and in order to inform HIV-1 vaccine design.

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.

Glycan Positioning Impacts HIV-1 Env Glycan-Shield Density, Function, and Recognition by Antibodies

HIV-1 envelope (Env) N-glycosylation impact virus-cell entry and immune evasion. How each glycan interacts to shape the Env-protein-sugar complex and affects Env function is not well understood. Here, analysis of two Env variants from the same donor, with differing functional characteristics and N-glycosylation-site composition, revealed that changes to key N-glycosylation sites affected the Env structure at distant locations and had a ripple effect on Env-wide glycan processing, virus infectivity, antibody recognition, and virus neutralization. Specifically, the N262 glycan, although not in the CD4-binding site, modulated Env binding to the CD4 receptor, affected Env recognition by several glycan-dependent neutralizing antibodies, and altered site-specific glycosylation heterogeneity, with, for example, N448 displaying limited glycan processing. Molecular-dynamic simulations visualized differences in glycan density and how specific oligosaccharide positions can move to compensate for a glycan loss. This study demonstrates how changes in individual glycans can alter molecular dynamics, processing, and function of the Env-glycan shield.

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 dynamics and network analysis of spike glycoprotein of SARS-COV-2

The ongoing pandemic caused by coronavirus SARS-COV-2 continues to rage with devastating consequences on human health and global economy. The spike glycoprotein on the surface of coronavirus mediates its entry into host cells and is the target of all current antibody design efforts to neutralize the virus. The glycan shield of the spike helps the virus to evade the human immune response by providing a thick sugar-coated barrier against any antibody. To study the dynamic motion of glycans in the spike protein, we performed microsecond-long MD simulation in two different states that correspond to the receptor binding domain in open or closed conformations. Analysis of this microsecond-long simulation revealed a scissoring motion on the N-terminal domain of neighboring monomers in the spike trimer. Role of multiple glycans in shielding of spike protein in different regions were uncovered by a network analysis, where the high betweenness centrality of glycans at the apex revealed their importance and function in the glycan shield. Microdomains of glycans were identified featuring a high degree of intra-communication in these microdomains. An antibody overlap analysis revealed the glycan microdomains as well as individual glycans that inhibit access to the antibody epitopes on the spike protein. Overall, the results of this study provide detailed understanding of the spike glycan shield, which may be utilized for therapeutic efforts against this crisis.

Networks of HIV-1 Envelope Glycans Maintain Antibody Epitopes in the Face of Glycan Additions and Deletions

Numerous broadly neutralizing antibodies (bnAbs) have been identified that target the glycans of the HIV-1 envelope spike. Neutralization breadth is notable given that glycan processing can be substantially influenced by the presence or absence of neighboring glycans. Here, using a stabilized recombinant envelope trimer, we investigate the degree to which mutations in the glycan network surrounding an epitope impact the fine glycan processing of antibody targets. Using cryo-electron microscopy and site-specific glycan analysis, we reveal the importance of glycans in the formation of the 2G12 bnAb epitope and show that the epitope is only subtly impacted by variations in the glycan network. In contrast, we show that the PG9 and PG16 glycan-based epitopes at the trimer apex are dependent on the presence of the highly conserved surrounding glycans. Glycan networks underpin the conservation of bnAb epitopes and are an important parameter in immunogen design.

Mass spectrometry for the identification and analysis of highly complex glycosylation of therapeutic or pathogenic proteins

INTRODUCTION

Protein glycosylation influences characteristics such as folding, stability, protein interactions, and solubility. Therefore, glycan moieties of therapeutic proteins and proteins that are likely associated with disease pathogenesis should be analyzed in-depth, including glycan heterogeneity and modification sites. Recent advances in analytical methods and instrumentation have enabled comprehensive characterization of highly complex glycosylated proteins.

AREA COVERED

The following aspects should be considered when analyzing glycosylated proteins: sample preparation, chromatographic separation, mass spectrometry (MS) and fragmentation methods, and bioinformatics, such as software solutions for data analyses. Notably, analysis of glycoproteins with heavily sialylated glycans or multiple glycosylation sites requires special considerations. Here, we discuss recent methodological advances in MS that provide detailed characterization of heterogeneous glycoproteins.

EXPERT OPINION

As characterization of complex glycosylated proteins is still analytically challenging, the function or pathophysiological significance of these proteins is not fully understood. To reproducibly produce desired forms of therapeutic glycoproteins or to fully elucidate disease-specific patterns of protein glycosylation, a highly reproducible and robust analytical platform(s) should be established. In addition to advances in MS instrumentation, optimization of analytical and bioinformatics methods and utilization of glycoprotein/glycopeptide standards is desirable. Ultimately, we envision that an automated high-throughput MS analysis will provide additional power to clinical studies and precision medicine.

Vulnerabilities in coronavirus glycan shields despite extensive glycosylation

Severe acute respiratory syndrome (SARS) and Middle East respiratory syndrome (MERS) coronaviruses (CoVs) are zoonotic pathogens with high fatality rates and pandemic potential. Vaccine development focuses on the principal target of the neutralizing humoral immune response, the spike (S) glycoprotein. Coronavirus S proteins are extensively glycosylated, encoding around 66-87 N-linked glycosylation sites per trimeric spike. Here, we reveal a specific area of high glycan density on MERS S that results in the formation of oligomannose-type glycan clusters, which were absent on SARS and HKU1 CoVs. We provide a comparison of the global glycan density of coronavirus spikes with other viral proteins including HIV-1 envelope, Lassa virus glycoprotein complex, and influenza hemagglutinin, where glycosylation plays a known role in shielding immunogenic epitopes. Overall, our data reveal how organisation of glycosylation across class I viral fusion proteins influence not only individual glycan compositions but also the immunological pressure across the protein surface.

Disruption of the HIV-1 Envelope allosteric network blocks CD4-induced rearrangements

The trimeric HIV-1 Envelope protein (Env) mediates viral-host cell fusion via a network of conformational transitions, with allosteric elements in each protomer orchestrating host receptor-induced exposure of the co-receptor binding site and fusion elements. To understand the molecular details of this allostery, here, we introduce Env mutations aimed to prevent CD4-induced rearrangements in the HIV-1 BG505 Env trimer. Binding analysis and single-molecule Förster Resonance Energy Transfer confirm that these mutations prevent CD4-induced transitions of the HIV-1 Env. Structural analysis by single-particle cryo-electron microscopy performed on the BG505 SOSIP mutant Env proteins shows rearrangements in the gp120 topological layer contacts with gp41. Displacement of a conserved tryptophan (W571) from its typical pocket in these Env mutants renders the Env insensitive to CD4 binding. These results reveal the critical function of W571 as a conformational switch in Env allostery and receptor-mediated viral entry and provide insights on Env conformation that are relevant for vaccine design.

Systematic Interaction Analysis of Anti-Human Immunodeficiency Virus Type-1 Neutralizing Antibodies with High Mannose Glycans Using Fragment Molecular Orbital and Molecular Dynamics Methods

A series of broadly neutralizing antibodies called PGT have been shown to be bound directly to human immunodeficiency virus type-1 via high mannose glycans on glycoprotein gp120. Despite the sequence similarities of amino acids of the antibodies, their affinities to the glycan differ. Glycan-antibody interactions among these antibodies are systematically compared with quantum chemical fragment molecular orbital calculations and molecular dynamics simulations. The differences among structural stability of the glycan in the active site of the complexes and total interaction energies as well as binding free energies between the glycan and antibodies agree well with the experimentally shown affinities of the glycan to the antibodies. The terminal saccharide, Man D3, is structurally stable and responsible for the glycan-antibody binding through electrostatic and dispersion interactions. The structural stability of nonterminal saccharides such as Man 4 or Man C plays substantial roles in the interaction via direct hydrogen bonds. © 2019 Wiley Periodicals, Inc.

Griffithsin Retains Anti-HIV-1 Potency with Changes in gp120 Glycosylation and Complements Broadly Neutralizing Antibodies PGT121 and PGT126

Griffithsin (Grft) is an antiviral lectin that has been shown to potently inhibit HIV-1 by binding high-mannose N-linked glycosylation sites on HIV-1 gp120. A key factor for Grft potency is glycosylation at N295 of gp120, which is directly adjacent to N332, a target glycan for an entire class of broadly neutralizing antibodies (bNAbs). Here, we unify previous work on the importance of other glycans to Grft potency against HIV-1 and Grft's role in mediating the conformational change of gp120 by mutating nearly every glycosylation site in gp120. In addition to a significant loss of Grft activity by the removal of glycosylation at N295, glycan absence at N332 or N448 was found to have moderate effects on Grft potency. Interestingly, in the absence of N295, Grft effectiveness could be improved by a mutation that results in the glycan at N448 shifting to N446, indicating that the importance of individual glycans may be related to their effect on glycosylation density. Grft's ability to alter the structure of gp120, exposing the CD4 binding site, correlated with the presence of glycosylation at N295 only in clade B strains, not clade C strains. We further demonstrate that Grft can rescue the activity of the bNAbs PGT121 and PGT126 in the event of a loss or a shift of glycosylation at N332, where the bNAbs suffer a drastic loss of potency. Despite targeting the same region, Grft in combination with PGT121 and PGT126 produced additive effects. This indicates that Grft could be an important combinational therapeutic.

Completeness of HIV-1 Envelope Glycan Shield at Transmission Determines Neutralization Breadth

Densely arranged N-linked glycans shield the HIV-1 envelope (Env) trimer from antibody recognition. Strain-specific breaches in this shield (glycan holes) can be targets of vaccine-induced neutralizing antibodies that lack breadth. To understand the interplay between glycan holes and neutralization breadth in HIV-1 infection, we developed a sequence-and structure-based approach to identify glycan holes for individual Env sequences that are shielded in most M-group viruses. Applying this approach to 12 longitudinally followed individuals, we found that transmitted viruses with more intact glycan shields correlated with development of greater neutralization breadth. Within 2 years, glycan acquisition filled most glycan holes present at transmission, indicating escape from hole-targeting neutralizing antibodies. Glycan hole filling generally preceded the time to first detectable breadth, although time intervals varied across hosts. Thus, completely glycan-shielded viruses were associated with accelerated neutralization breadth development, suggesting that Env immunogens with intact glycan shields may be preferred components of AIDS vaccines.

Wagh et al. show that transmitted viruses with more intact glycan shields are correlated with development of neutralization breadth in HIV-1-infected individuals. This is consistent with previous findings that glycan holes in Env immunogens are targeted by strain-specific neutralizing responses, and suggests that immunogens with intact glycan Shields may be advantageous.

Glycosylator: a Python framework for the rapid modeling of glycans

BACKGROUND

Carbohydrates are a class of large and diverse biomolecules, ranging from a simple monosaccharide to large multi-branching glycan structures. The covalent linkage of a carbohydrate to the nitrogen atom of an asparagine, a process referred to as N-linked glycosylation, plays an important role in the physiology of many living organisms. Most software for glycan modeling on a personal desktop computer requires knowledge of molecular dynamics to interface with specialized programs such as CHARMM or AMBER. There are a number of popular web-based tools that are available for modeling glycans (e.g., GLYCAM-WEB (http:// https://dev.glycam.org/gp/ ) or Glycosciences.db ( http://www.glycosciences.de/ )). However, these web-based tools are generally limited to a few canonical glycan conformations and do not allow the user to incorporate glycan modeling into their protein structure modeling workflow.

RESULTS

Here, we present Glycosylator, a Python framework for the identification, modeling and modification of glycans in protein structure that can be used directly in a Python script through its application programming interface (API) or through its graphical user interface (GUI). The GUI provides a straightforward two-dimensional (2D) rendering of a glycoprotein that allows for a quick visual inspection of the glycosylation state of all the sequons on a protein structure. Modeled glycans can be further refined by a genetic algorithm for removing clashes and sampling alternative conformations. Glycosylator can also identify specific three-dimensional (3D) glycans on a protein structure using a library of predefined templates.

CONCLUSIONS

Glycosylator was used to generate models of glycosylated protein without steric clashes. Since the molecular topology is based on the CHARMM force field, new complex sugar moieties can be generated without modifying the internals of the code. Glycosylator provides more functionality for analyzing and modeling glycans than any other available software or webserver at present. Glycosylator will be a valuable tool for the glycoinformatics and biomolecular modeling communities.

Illuminating the virus life cycle with single-molecule FRET imaging

Single-molecule Förster resonance energy transfer (smFRET) imaging has emerged as a powerful tool to probe conformational dynamics of viral proteins, identify novel structural intermediates that are hiding in averaging population-based measurements, permit access to the energetics of transitions and as such to the precise molecular mechanisms of viral replication. One strength of smFRET is the capability of characterizing biological molecules in their fully hydrated/native state, which are not necessarily available to other structural methods. Elegant experimental design for physiologically relevant conditions, such as intact virions, has permitted the detection of previously unknown conformational states of viral glycoproteins, revealed asymmetric intermediates, and allowed access to the real-time imaging of conformational changes during viral fusion. As more laboratories are applying smFRET, our understanding of the molecular mechanisms and the dynamic nature of viral proteins throughout the virus life cycle are predicted to improve and assist the development of novel antiviral therapies and vaccine design.

Glycosylation of viral surface proteins probed by mass spectrometry

Glycosylation is a common and biologically significant post-translational modification that is found on numerous virus surface proteins (VSPs). Many of these glycans affect virulence through modulating virus receptor binding, masking antigenic sites, or by stimulating the host immune response. Mass spectrometry (MS) has arisen as a pivotal technique for the characterization of VSP glycosylation. This review will cover how MS-based analyses, such as released glycan profiles, glycan site localization, site-occupancy, and site-specific heterogeneity, are being utilized to map VSP glycosylation. Furthermore, this review will provide information on how MS glycoprofiling data are being used in conjunction with molecular and structural experiments to provide a better understanding of the role of specific glycans in VSP function.

Defining HIV-1 Envelope N-Glycan Microdomains through Site-Specific Heterogeneity Profiles

The HIV-1 envelope (Env) glycans shield the surface of Env from the immune system and form integral interactions important for a functional Env. To understand how individual N-glycosylation sites (NGS) coordinate to form a dynamic shield and evade the immune system through mutations, we tracked 20 NGS in Env from HIV-transmitted/founder (T/F) and immune escape variants and their mutants involving the N262 glycan. NGS were profiled in a site-specific manner using a high-resolution mass spectrometry (MS)-based workflow. Using this site-specific quantitative heterogeneity profiling, we empirically characterized the interdependent NGS of a microdomain in the high-mannose patch (HMP). The changes (shifts) in NGS heterogeneity between the T/F and immune escape variants defined a range of NGS that we further probed for exclusive combinations of sequons in the HMP microdomain using the Los Alamos National Laboratory HIV sequence database. The resultant sequon combinations, including the highly conserved NGS N262, N448, and N301, created an immune escape map of the conserved and variable sequons in the HMP microdomain. This report provides details on how some clustered NGS form microdomains that can be identified and tracked across Env variants. These microdomains have a limited number of N-glycan-sequon combinations that may allow the anticipation of immune escape variants.IMPORTANCE The Env protein of HIV is highly glycosylated, and the sites of glycosylation can change as the virus mutates during immune evasion. Due to these changes, the glycan location and heterogeneity of surrounding N-glycosylation sites can be altered, resulting in exposure of different glycan or proteoglycan surfaces while still producing a viable HIV variant. These changes present a need for vaccine developers to identify Env variants with epitopes most likely to induce durable protective responses. Here we describe a means of anticipating HIV-1 immune evasion by dividing Env into N-glycan microdomains that have a limited number of N-glycan sequon combinations.

Structural Rearrangements Maintain the Glycan Shield of an HIV-1 Envelope Trimer After the Loss of a Glycan

The HIV-1 envelope (Env) glycoprotein is the primary target of the humoral immune response and a critical vaccine candidate. However, Env is densely glycosylated and thereby substantially protected from neutralisation. Importantly, glycan N301 shields V3 loop and CD4 binding site epitopes from neutralising antibodies. Here, we use molecular dynamics techniques to evaluate the structural rearrangements that maintain the protective qualities of the glycan shield after the loss of glycan N301. We examined a naturally occurring subtype C isolate and its N301A mutant; the mutant not only remained protected against neutralising antibodies targeting underlying epitopes, but also exhibited an increased resistance to the VRC01 class of broadly neutralising antibodies. Analysis of this mutant revealed several glycans that were responsible, independently or through synergy, for the neutralisation resistance of the mutant. These data provide detailed insight into the glycan shield’s ability to compensate for the loss of a glycan, as well as the cascade of glycan movements on a protomer, starting at the point mutation, that affects the integrity of an antibody epitope located at the edge of the diminishing effect. These results present key, previously overlooked, considerations for HIV-1 Env glycan research and related vaccine studies.

Predicting the Structures of Glycans, Glycoproteins, and Their Complexes

Complex carbohydrates are ubiquitous in nature, and together with proteins and nucleic acids they comprise the building blocks of life. But unlike proteins and nucleic acids, carbohydrates form nonlinear polymers, and they are not characterized by robust secondary or tertiary structures but rather by distributions of well-defined conformational states. Their molecular flexibility means that oligosaccharides are often refractory to crystallization, and nuclear magnetic resonance (NMR) spectroscopy augmented by molecular dynamics (MD) simulation is the leading method for their characterization in solution. The biological importance of carbohydrate-protein interactions, in organismal development as well as in disease, places urgency on the creation of innovative experimental and theoretical methods that can predict the specificity of such interactions and quantify their strengths. Additionally, the emerging realization that protein glycosylation impacts protein function and immunogenicity places the ability to define the mechanisms by which glycosylation impacts these features at the forefront of carbohydrate modeling. This review will discuss the relevant theoretical approaches to studying the three-dimensional structures of this fascinating class of molecules and interactions, with reference to the relevant experimental data and techniques that are key for validation of the theoretical predictions.

All-atom virus simulations

The constant threat of viral disease can be combated by the development of novel vaccines and therapeutics designed to disrupt key features of virus structure or infection cycle processes. Such development relies on high-resolution characterization of viruses and their dynamical behaviors, which are often challenging to obtain solely by experiment. In response, all-atom molecular dynamics simulations are widely leveraged to study the structural components of viruses, leading to some of the largest simulation endeavors undertaken to date. The present work reviews exemplary all-atom simulation work on viruses, as well as progress toward simulating entire virions.

A Neutralizing Antibody Recognizing Primarily N-Linked Glycan Targets the Silent Face of the HIV Envelope

Virtually the entire surface of the HIV-1-envelope trimer is recognized by neutralizing antibodies, except for a highly glycosylated region at the center of the "silent face" on the gp120 subunit. From an HIV-1-infected donor, #74, we identified antibody VRC-PG05, which neutralized 27% of HIV-1 strains. The crystal structure of the antigen-binding fragment of VRC-PG05 in complex with gp120 revealed an epitope comprised primarily of N-linked glycans from N262, N295, and N448 at the silent face center. Somatic hypermutation occurred preferentially at antibody residues that interacted with these glycans, suggesting somatic development of glycan recognition. Resistance to VRC-PG05 in donor #74 involved shifting of glycan-N448 to N446 or mutation of glycan-proximal residue E293. HIV-1 neutralization can thus be achieved at the silent face center by glycan-recognizing antibody; along with other known epitopes, the VRC-PG05 epitope completes coverage by neutralizing antibody of all major exposed regions of the prefusion closed trimer.

Integrity of Glycosylation Processing of a Glycan-Depleted Trimeric HIV-1 Immunogen Targeting Key B-Cell Lineages

Broadly neutralizing antibodies (bNAbs) that target the trimeric HIV-1 envelope glycoprotein spike (Env) are tools that can guide the design of recombinant Env proteins intended to engage the predicted human germline precursors of bNAbs (gl-bNAbs). The protein components of gl-bNAb epitopes are often masked by glycans, while mature bNAbs can evolve to accommodate or bypass these shielding glycans. The design of germline-targeting Env immunogens therefore includes the targeted deletion of specific glycan sites. However, the processing of glycans on Env trimers can be influenced by the density with which they are packed together, a highly relevant point given the essential contributions under-processed glycans make to multiple bNAb epitopes. We sought to determine the impact of the removal of 15 potential N-glycan sites (5 per protomer) from the germline-targeting soluble trimer, BG505 SOSIP.v4.1-GT1, using quantitative, site-specific N-glycan mass spectrometry analysis. We find that, compared with SOSIP.664, there was little overall change in the glycan profile but only subtle increases in the extent of processing at sites immediately adjacent to where glycans had been deleted. We conclude that multiple glycans can be deleted from BG505 SOSIP trimers without perturbing the overall integrity of the glycan shield.

Proper balance of solvent-solute and solute-solute interactions in the treatment of the diffusion of glucose using the Drude polarizable force field

Motivated by underestimation of the diffusion constant of glucose in aqueous solution at high glucose concentrations we performed additional optimization of the Drude polarizable hexopyranose monosaccharide force field. This indicated aggregation of the glucose at higher concentrations, which is a concern for studies of complex glycan systems such as the HIV Envelope where high effective concentrations of sugars are present. High-level quantum mechanical calculations were undertaken on water monohydrate-glucose interactions, on water cluster-glucose interactions and on glucose-glucose dimers in stacked (parallel) and perpendicular orientations. Optimization of the nonbond and dihedral parameters targeting these data yielded a revised model that showed improved agreement with experimental aqueous diffusion data. However, limitations in the diffusion constants were still present. These were due to the SWM4-NDP inherently overestimating the diffusion constant of water, a problem that was validated by calculation of the aqueous diffusion constants using the SWM6-NDP water model. In addition, results show the water diffusion constant to be significantly overestimated at high glucose concentrations though the glucose diffusion is in satisfactory agreement with experiment. These results indicate the subtle balance of water-sugar, water-water and sugar-sugar interactions that needs to be properly modeled to account for the full range of aqueous behavior of sugars in aqueous solution.