Site-Specific Steric Control of SARS-CoV-2 Spike Glycosylation.

A central tenet in the design of vaccines is the display of native-like antigens in the elicitation of protective immunity. The abundance of N-linked glycans across the SARS-CoV-2 spike protein is a potential source of heterogeneity among the many different vaccine candidates under investigation. Here, we investigate the glycosylation of recombinant SARS-CoV-2 spike proteins from five different laboratories and compare them against S protein from infectious virus, cultured in Vero cells. We find patterns that are conserved across all samples, and this can be associated with site-specific stalling of glycan maturation that acts as a highly sensitive reporter of protein structure. Molecular dynamics simulations of a fully glycosylated spike support a model of steric restrictions that shape enzymatic processing of the glycans. These results suggest that recombinant spike-based SARS-CoV-2 immunogen glycosylation reproducibly recapitulates signatures of viral glycosylation.

Atomistic and Coarse Grain Simulations of the Cell Envelope of Gram-Negative Bacteria: What Have We Learned?

Bacterial membranes, and those of Gram-negative bacteria in particular, are some of the most biochemically diverse membranes known. They incorporate a wide range of lipid types and proteins of varying sizes, architectures, and functions. While simpler biological membranes have been the focus of myriad simulation studies over the years that have yielded invaluable details to complement, and often to direct, ongoing experimental studies, simulations of complex bacterial membranes have been slower to emerge. However, the past few years have seen tremendous activity in this area, leading to advances such as the development of atomistic and coarse-grain models of the lipopolysaccharide (LPS) component of the outer membrane that are compatible with widely used simulation codes. In this Account, we review our contributions to the field of molecular simulations of the bacterial cell envelope, including the development of models of both membranes and the cell wall of Gram-negative bacteria, with a predominant focus on E. coli. At the atomistic level, simulations of chemically accurate models of both membranes have revealed the tightly cross-linked nature of the LPS headgroups and have shown that penetration of solutes through these regions is not as straightforward as the route through phospholipids. The energetic differences between the two routes have been calculated. Simulations of native outer membrane proteins in LPS-containing membranes have shown that the conformational dynamics of the proteins is not only slower in LPS but also different compared to in simpler models of phospholipid bilayers. These chemically more complex and consequently biologically more relevant models are leading to details of conformational dynamics that were previously inaccessible from simulations. Coarse-grain models have enabled simulations of multiprotein systems on time scales of microseconds, leading to insights not only into the rates of protein and lipid diffusion but also into the trends in their respective directions of flow. We find that the motions of LPS molecules are highly correlated with each other but also with outer membrane proteins embedded within the membrane. We have shown that the two leaflets of the outer membrane exhibit communication, whereby regions of low disorder in one leaflet correspond to regions of high disorder in the other. The cell wall remains a comparatively neglected component, although models of the E. coli peptidoglycan are now emerging, particularly at the atomistic level. Our simulations of Braun's lipoprotein have shown that bending and tilting of this protein afford a degree of variability in the gap between the cell wall and the OM. The noncovalent interactions with the cell wall of proteins such as OmpA can further influence the width of this gap by extension or contraction of their linker domains. Overall we have shown that the dynamics of proteins, lipids, and other molecular species within the outer membrane cannot be approximated using simpler phospholipid bilayers, if one is addressing questions regarding the in vivo behavior of Gram-negative bacteria. These membranes have their own unique chemical characteristics that cannot be decoupled from their biological functions.

Directional Porin Binding of Intrinsically Disordered Protein Sequences Promotes Colicin Epitope Display in the Bacterial Periplasm.

Protein bacteriocins are potent narrow spectrum antibiotics that exploit outer membrane porins to kill bacteria by poorly understood mechanisms. Here, we determine how colicins, bacteriocins specific for Escherichia coli, engage the trimeric porin OmpF to initiate toxin entry. The N-terminal ∼80 residues of the nuclease colicin ColE9 are intrinsically unstructured and house two OmpF binding sites (OBS1 and OBS2) that reside within the pores of OmpF and which flank an epitope that binds periplasmic TolB. Using a combination of molecular dynamics simulations, chemical trimerization, isothermal titration calorimetry, fluorescence microscopy, and single channel recording planar lipid bilayer measurements, we show that this arrangement is achieved by OBS2 binding from the extracellular face of OmpF, while the interaction of OBS1 occurs from the periplasmic face of OmpF. Our study shows how the narrow pores of oligomeric porins are exploited by colicin disordered regions for direction-specific binding, which ensures the constrained presentation of an activating signal within the bacterial periplasm.

Braun's Lipoprotein Facilitates OmpA Interaction with the Escherichia coli Cell Wall.

Gram-negative bacteria such as Escherichia coli are protected by a complex cell envelope. The development of novel therapeutics against these bacteria necessitates a molecular level understanding of the structure-dynamics-function relationships of the various components of the cell envelope. We use atomistic MD simulations to reveal the details of covalent and noncovalent protein interactions that link the outer membrane to the aqueous periplasmic region. We show that the Braun's lipoprotein tilts and bends, and thereby lifts the cell wall closer to the outer membrane. Both monomers and dimers of the outer membrane porin OmpA can interact with peptidoglycan in the presence of Braun's lipoprotein, but in the absence of the latter, only dimers of OmpA show a propensity to form contacts with peptidoglycan. Our study provides a glimpse of how the molecular components of the bacterial cell envelope interact with each other to mediate cell wall attachment in E. coli.

Full-Length OmpA: Structure, Function, and Membrane Interactions Predicted by Molecular Dynamics Simulations.

OmpA is a multidomain protein found in the outer membranes of most Gram-negative bacteria. Despite a wealth of reported structural and biophysical studies, the structure-function relationships of this protein remain unclear. For example, it is still debated whether it functions as a pore, and the precise molecular role it plays in attachment to the peptidoglycan of the periplasm is unknown. The absence of a consensus view is partly due to the lack of a complete structure of the full-length protein. To address this issue, we performed molecular-dynamics simulations of the full-length model of the OmpA dimer proposed by Robinson and co-workers. The N-terminal domains were embedded in an asymmetric model of the outer membrane, with lipopolysaccharide molecules in the outer leaflet and phospholipids in the inner leaflet. Our results reveal a large dimerization interface within the membrane environment, ensuring that the dimer is stable over the course of the simulations. The linker is flexible, expanding and contracting to pull the globular C-terminal domain up toward the membrane or push it down toward the periplasm, suggesting a possible mechanism for providing mechanical stability to the cell. The external loops were more stabilized than was observed in previous studies due to the extensive dimerization interface and presence of lipopolysaccharide molecules in our outer-membrane model, which may have functional consequences in terms of OmpA adhesion to host cells. In addition, the pore-gating behavior of the protein was modulated compared with previous observations, suggesting a possible role for dimerization in channel regulation.

Mechanisms of allostery at the viral surface through the eyes of molecular simulation.

The outermost surface layer of any virus is formed by either a capsid shell or envelope. Such layers have traditionally been thought of as immovable structures, but it is becoming apparent that they cannot be viewed exclusively as static architectures protecting the viral genome. A limited number of proteins on the virion surface must perform a multitude of functions in order to orchestrate the viral life cycle, and allostery can regulate their structures at multiple levels of organization, spanning individual molecules, protomers, large oligomeric assemblies, or entire viral surfaces. Here, we review recent contributions from the molecular simulation field to viral surface allostery, with a particular focus on the trimeric spike glycoprotein emerging from the coronavirus surface, and the icosahedral flaviviral envelope complex. As emerging viral pathogens continue to pose a global threat, an improved understanding of viral dynamics and allosteric regulation will prove crucial in developing novel therapeutic strategies.

Membrane platform protein PulF of the Klebsiella type II secretion system forms a trimeric ion channel essential for endopilus assembly and protein secretion.

IMPORTANCE: Type IV pili and type II secretion systems are members of the widespread type IV filament (T4F) superfamily of nanomachines that assemble dynamic and versatile surface fibers in archaea and bacteria. The assembly and retraction of T4 filaments with diverse surface properties and functions require the plasma membrane platform proteins of the GspF/PilC superfamily. Generally considered dimeric, platform proteins are thought to function as passive transmitters of the mechanical energy generated by the ATPase motor, to somehow promote insertion of pilin subunits into the nascent pilus fibers. Here, we generate and experimentally validate structural predictions that support the trimeric state of a platform protein PulF from a type II secretion system. The PulF trimers form selective proton or sodium channels which might energize pilus assembly using the membrane potential. The conservation of the channel sequence and structural features implies a common mechanism for all T4F assembly systems. We propose a model of the oligomeric PulF-PulE ATPase complex that provides an essential framework to investigate and understand the pilus assembly mechanism.

SARS-CoV-2 spike protein as a bacterial lipopolysaccharide delivery system in an overzealous inflammatory cascade.

Accumulating evidence indicates a potential role for bacterial lipopolysaccharide (LPS) in the overactivation of the immune response during SARS-CoV-2 infection. LPS is recognized by Toll-like receptor 4, mediating proinflammatory effects. We previously reported that LPS directly interacts with SARS-CoV-2 spike (S) protein and enhances proinflammatory activities. Using native gel electrophoresis and hydrogen-deuterium exchange mass spectrometry, we showed that LPS binds to multiple hydrophobic pockets spanning both the S1 and S2 subunits of the S protein. Molecular simulations validated by a microscale thermophoresis binding assay revealed that LPS binds to the S2 pocket with a lower affinity compared to S1, suggesting a role as an intermediate in LPS transfer. Congruently, nuclear factor-kappa B (NF-κB) activation in monocytic THP-1 cells is strongly boosted by S2. Using NF-κB reporter mice followed by bioimaging, a boosting effect was observed for both S1 and S2, with the former potentially facilitated by proteolysis. The Omicron S variant binds to LPS, but with reduced affinity and LPS boosting in vitro and in vivo. Taken together, the data provide a molecular mechanism by which S protein augments LPS-mediated hyperinflammation.

Bioactive Suture with Added Innate Defense Functionality for the Reduction of Bacterial Infection and Inflammation.

Surgical site infections (SSI) are a clinical and economic burden. Suture-associated SSI may develop when bacteria colonize the suture surface and form biofilms that are resistant to antibiotics. Thrombin-derived C-terminal peptide (TCP)-25 is a host defense peptide with a unique dual mode of action that can target both bacteria and the excessive inflammation induced by bacterial products. The peptide demonstrates therapeutic potential in preclinical in vivo wound infection models. In this study, the authors set out to explore whether TCP-25 can provide a new bioactive innate immune feature to hydrophilic polyglactin sutures (Vicryl). Using a combination of biochemical, biophysical, antibacterial, biofilm, and anti-inflammatory assays in vitro, in silico molecular modeling studies, along with experimental infection and inflammation models in mice, a proof-of-concept that TCP-25 can provide Vicryl sutures with a previously undisclosed host defense capacity, that enables targeting of bacteria, biofilms, and the accompanying inflammatory response, is shown.

Spike-Independent Infection of Human Coronavirus 229E in Bat Cells.

Bats are the reservoir for numerous human pathogens, including coronaviruses. Despite many coronaviruses having descended from bat ancestors, little is known about virus-host interactions and broader evolutionary history involving bats. Studies have largely focused on the zoonotic potential of coronaviruses with few infection experiments conducted in bat cells. To determine genetic changes derived from replication in bat cells and possibly identify potential novel evolutionary pathways for zoonotic virus emergence, we serially passaged six human 229E isolates in a newly established Rhinolophus lepidus (horseshoe bat) kidney cell line. Here, we observed extensive deletions within the spike and open reading frame 4 (ORF4) genes of five 229E viruses after passaging in bat cells. As a result, spike protein expression and infectivity of human cells was lost in 5 of 6 viruses, but the capability to infect bat cells was maintained. Only viruses that expressed the spike protein could be neutralized by 229E spike-specific antibodies in human cells, whereas there was no neutralizing effect on viruses that did not express the spike protein inoculated on bat cells. However, one isolate acquired an early stop codon, abrogating spike expression but maintaining infection in bat cells. After passaging this isolate in human cells, spike expression was restored due to acquisition of nucleotide insertions among virus subpopulations. Spike-independent infection of human coronavirus 229E may provide an alternative mechanism for viral maintenance in bats that does not rely on the compatibility of viral surface proteins and known cellular entry receptors. IMPORTANCE Many viruses, including coronaviruses, originated from bats. Yet, we know little about how these viruses switch between hosts and enter human populations. Coronaviruses have succeeded in establishing in humans at least five times, including endemic coronaviruses and the recent severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). In an approach to identify requirements for host switches, we established a bat cell line and adapted human coronavirus 229E viruses by serial passage. The resulting viruses lost their spike protein but maintained the ability to infect bat cells, but not human cells. Maintenance of 229E viruses in bat cells appears to be independent of a canonical spike receptor match, which in turn might facilitate cross-species transmission in bats.

Targeting Toll-like receptor-driven systemic inflammation by engineering an innate structural fold into drugs.

There is a clinical need for conceptually new treatments that target the excessive activation of inflammatory pathways during systemic infection. Thrombin-derived C-terminal peptides (TCPs) are endogenous anti-infective immunomodulators interfering with CD14-mediated TLR-dependent immune responses. Here we describe the development of a peptide-based compound for systemic use, sHVF18, expressing the evolutionarily conserved innate structural fold of natural TCPs. Using a combination of structure- and in silico-based design, nuclear magnetic resonance spectroscopy, biophysics, mass spectrometry, cellular, and in vivo studies, we here elucidate the structure, CD14 interactions, protease stability, transcriptome profiling, and therapeutic efficacy of sHVF18. The designed peptide displays a conformationally stabilized, protease resistant active innate fold and targets the LPS-binding groove of CD14. In vivo, it shows therapeutic efficacy in experimental models of endotoxin shock in mice and pigs and increases survival in mouse models of systemic polymicrobial infection. The results provide a drug class based on Nature´s own anti-infective principles.

Defining neutralization and allostery by antibodies against COVID-19 variants.

The changing landscape of SARS-CoV-2 Spike protein is linked to the emergence of variants, immune-escape and reduced efficacy of the existing repertoire of anti-viral antibodies. The functional activity of neutralizing antibodies is linked to their quaternary changes occurring as a result of antibody-Spike trimer interactions. Here, we reveal the conformational dynamics and allosteric perturbations linked to binding of novel human antibodies and the viral Spike protein. We identified epitope hotspots, and associated changes in Spike dynamics that distinguish weak, moderate and strong neutralizing antibodies. We show the impact of mutations in Wuhan-Hu-1, Delta, and Omicron variants on differences in the antibody-induced conformational changes in Spike and illustrate how these render certain antibodies ineffective. Antibodies with similar binding affinities may induce destabilizing or stabilizing allosteric effects on Spike, with implications for neutralization efficacy. Our results provide mechanistic insights into the functional modes and synergistic behavior of human antibodies against COVID-19 and may assist in designing effective antiviral strategies.

SARS-CoV-2 spike protein aggregation is triggered by bacterial lipopolysaccharide.

SARS-CoV-2 spike (S) protein is crucial for virus invasion in COVID-19. Here, we showed that lipopolysaccharide (LPS) can trigger S protein aggregation at high doses of LPS and S protein. We demonstrated the formation of S protein aggregates by microscopy analyses, aggregation and gel shift assays. LPS at high levels boosts the formation of S protein aggregates as detected by amytracker and thioflavin T dyes that specifically bind to aggregating proteins. We validated the role of LPS by blocking the formation of aggregates by the endotoxin-scavenging thrombin-derived peptide TCP-25. Aggregation-prone sequences in S protein are predicted to be nearby LPS binding sites, while molecular simulations showed stable formation of S protein-LPS higher-order oligomers. Collectively, our results provide evidence of LPS-induced S protein aggregation.

Allosteric perspective on the mutability and druggability of the SARS-CoV-2 Spike protein.

Recent developments in the SARS-CoV-2 pandemic point to its inevitable transformation into an endemic disease, urging both refinement of diagnostics for emerging variants of concern (VOCs) and design of variant-specific drugs in addition to vaccine adjustments. Exploring the structure and dynamics of the SARS-CoV-2 Spike protein, we argue that the high-mutability characteristic of RNA viruses coupled with the remarkable flexibility and dynamics of viral proteins result in a substantial involvement of allosteric mechanisms. While allosteric effects of mutations should be considered in predictions and diagnostics of new VOCs, allosteric drugs advantageously avoid escape mutations via non-competitive inhibition originating from alternative distal locations. The exhaustive allosteric signaling and probing maps presented herein provide a comprehensive picture of allostery in the spike protein, making it possible to locate potential mutations that could work as new VOC "drivers" and to determine binding patches that may be targeted by newly developed allosteric drugs.

Uncovering cryptic pockets in the SARS-CoV-2 spike glycoprotein.

The COVID-19 pandemic has prompted a rapid response in vaccine and drug development. Herein, we modeled a complete membrane-embedded SARS-CoV-2 spike glycoprotein and used molecular dynamics simulations with benzene probes designed to enhance discovery of cryptic pockets. This approach recapitulated lipid and host metabolite binding sites previously characterized by cryo-electron microscopy, revealing likely ligand entry routes, and uncovered a novel cryptic pocket with promising druggable properties located underneath the 617-628 loop. A full representation of glycan moieties was essential to accurately describe pocket dynamics. A multi-conformational behavior of the 617-628 loop in simulations was validated using hydrogen-deuterium exchange mass spectrometry experiments, supportive of opening and closing dynamics. The pocket is the site of multiple mutations associated with increased transmissibility found in SARS-CoV-2 variants of concern including Omicron. Collectively, this work highlights the utility of the benzene mapping approach in uncovering potential druggable sites on the surface of SARS-CoV-2 targets.

Glycosylation and Serological Reactivity of an Expression-enhanced SARS-CoV-2 Viral Spike Mimetic.

Extensive glycosylation of viral glycoproteins is a key feature of the antigenic surface of viruses and yet glycan processing can also be influenced by the manner of their recombinant production. The low yields of the soluble form of the trimeric spike (S) glycoprotein from SARS-CoV-2 has prompted advances in protein engineering that have greatly enhanced the stability and yields of the glycoprotein. The latest expression-enhanced version of the spike incorporates six proline substitutions to stabilize the prefusion conformation (termed SARS-CoV-2 S HexaPro). Although the substitutions greatly enhanced expression whilst not compromising protein structure, the influence of these substitutions on glycan processing has not been explored. Here, we show that the site-specific N-linked glycosylation of the expression-enhanced HexaPro resembles that of an earlier version containing two proline substitutions (2P), and that both capture features of native viral glycosylation. However, there are site-specific differences in glycosylation of HexaPro when compared to 2P. Despite these discrepancies, analysis of the serological reactivity of clinical samples from infected individuals confirmed that both HexaPro and 2P protein are equally able to detect IgG, IgA, and IgM responses in all sera analysed. Moreover, we extend this observation to include an analysis of glycan engineered S protein, whereby all N-linked glycans were converted to oligomannose-type and conclude that serological activity is not impacted by large scale changes in glycosylation. These observations suggest that variations in glycan processing will not impact the serological assessments currently being performed across the globe.

The impact of Gag non-cleavage site mutations on HIV-1 viral fitness from integrative modelling and simulations.

The high mutation rate in retroviruses is one of the leading causes of drug resistance. In human immunodeficiency virus type-1 (HIV-1), synergistic mutations in its protease and the protease substrate - the Group-specific antigen (Gag) polyprotein - work together to confer drug resistance against protease inhibitors and compensate the mutations affecting viral fitness. Some Gag mutations can restore Gag-protease binding, yet most Gag-protease correlated mutations occur outside of the Gag cleavage site. To investigate the molecular basis for this, we now report multiscale modelling approaches to investigate various sequentially cleaved Gag products in the context of clinically relevant mutations that occur outside of the cleavage sites, including simulations of the largest Gag proteolytic product in its viral membrane-bound state. We found that some mutations, such as G123E and H219Q, involve direct interaction with cleavage site residues to influence their local environment, while certain mutations in the matrix domain lead to the enrichment of lipids important for Gag targeting and assembly. Collectively, our results reveal why non-cleavage site mutations have far-reaching implications outside of Gag proteolysis, with important consequences for drugging Gag maturation intermediates and tackling protease inhibitor resistance.

SARS-CoV-2 S protein:ACE2 interaction reveals novel allosteric targets.

The spike (S) protein is the main handle for SARS-CoV-2 to enter host cells via surface angiotensin-converting enzyme 2 (ACE2) receptors. How ACE2 binding activates proteolysis of S protein is unknown. Here, using amide hydrogen-deuterium exchange mass spectrometry and molecular dynamics simulations, we have mapped the S:ACE2 interaction interface and uncovered long-range allosteric propagation of ACE2 binding to sites necessary for host-mediated proteolysis of S protein, critical for viral host entry. Unexpectedly, ACE2 binding enhances dynamics at a distal S1/S2 cleavage site and flanking protease docking site ~27 Å away while dampening dynamics of the stalk hinge (central helix and heptad repeat [HR]) regions ~130 Å away. This highlights that the stalk and proteolysis sites of the S protein are dynamic hotspots in the prefusion state. Our findings provide a dynamics map of the S:ACE2 interface in solution and also offer mechanistic insights into how ACE2 binding is allosterically coupled to distal proteolytic processing sites and viral-host membrane fusion. Thus, protease docking sites flanking the S1/S2 cleavage site represent alternate allosteric hotspot targets for potential therapeutic development.

Site-specific steric control of SARS-CoV-2 spike glycosylation.

A central tenet in the design of vaccines is the display of native-like antigens in the elicitation of protective immunity. The abundance of N-linked glycans across the SARS-CoV-2 spike protein is a potential source of heterogeneity between the many different vaccine candidates under investigation. Here, we investigate the glycosylation of recombinant SARS-CoV-2 spike proteins from five different laboratories and compare them against infectious virus S protein. We find patterns which are conserved across all samples and this can be associated with site-specific stalling of glycan maturation which act as a highly sensitive reporter of protein structure. Molecular dynamics (MD) simulations of a fully glycosylated spike support s a model of steric restrictions that shape enzymatic processing of the glycans. These results suggest that recombinant spike-based SARS-CoV-2 immunogen glycosylation reproducibly recapitulates signatures of viral glycosylation.