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.

Structure and sucrose binding mechanism of the plant SUC1 sucrose transporter

Sucrose import from photosynthetic tissues into the phloem is mediated by transporters from the low-affinity sucrose transporter family (SUC/SUT family). Furthermore, sucrose redistribution to other tissues is driven by phloem sap movement, the product of high turgor pressure created by this import activity. Additionally, sink organs such as fruits, cereals and seeds that accumulate high concentrations of sugar also depend on this active transport of sucrose. Here we present the structure of the sucrose-proton symporter, Arabidopsis thaliana SUC1, in an outward open conformation at 2.7 Å resolution, together with molecular dynamics simulations and biochemical characterization. We identify the key acidic residue required for proton-driven sucrose uptake and describe how protonation and sucrose binding are strongly coupled. Sucrose binding is a two-step process, with initial recognition mediated by the glucosyl moiety binding directly to the key acidic residue in a stringent pH-dependent manner. Our results explain how low-affinity sucrose transport is achieved in plants, and pinpoint a range of SUC binders that help define selectivity. Our data demonstrate a new mode for proton-driven symport with links to cation-driven symport and provide a broad model for general low-affinity transport in highly enriched substrate environments.

A pH-dependent cluster of charges in a conserved cryptic pocket on flaviviral envelopes

Flaviviruses are enveloped viruses which include human pathogens that are predominantly transmitted by mosquitoes and ticks. Some, such as dengue virus, exhibit the phenomenon of antibody-dependent enhancement (ADE) of disease, making vaccine-based routes of fighting infections problematic. The pH-dependent conformational change of the envelope (E) protein required for fusion between the viral and endosomal membranes is an attractive point of inhibition by antivirals as it has the potential to diminish the effects of ADE. We examined six flaviviruses by employing large-scale molecular dynamics (MD) simulations of raft systems that represent a substantial portion of the flaviviral envelope. We utilised a benzene-mapping approach that led to a discovery of shared hotspots and conserved cryptic sites. A cryptic pocket previously shown to bind a detergent molecule exhibited strain-specific characteristics. An alternative conserved cryptic site at the E protein domain interfaces showed a consistent dynamic behaviour across flaviviruses and contained a conserved cluster of ionisable residues. Constant-pH simulations revealed cluster and domain-interface disruption under low pH conditions. Based on this, we propose a cluster-dependent mechanism that addresses inconsistencies in the histidine-switch hypothesis and highlights the role of cluster protonation in orchestrating the domain dissociation pivotal for the formation of the fusogenic trimer.

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.

Bridging the N-terminal and middle domains in FliG of the flagellar rotor

Flagella are necessary for bacterial movement and contribute to various aspects of virulence. They are complex cylindrical structures built of multiple molecular rings with self-assembly properties. The flagellar rotor is composed of the MS-ring and the C-ring. The FliG protein of the C-ring is central to flagellar assembly and function due to its roles in linking the C-ring with the MS-ring and in torque transmission from stator to rotor. No high-resolution structure of an assembled C-ring has been resolved to date, and the conformation adopted by FliG within the ring is unclear due to variations in available crystallographic data. Here, we use molecular dynamics (MD) simulations to study the conformation and dynamics of FliG in different states of assembly, including both in physiologically relevant and crystallographic lattice environments. We conclude that the linker between the FliG N-terminal and middle domain likely adopts an extended helical conformation in vivo, in contrast with the contracted conformation observed in some previous X-ray studies. We further support our findings with integrative model building of full-length FliG and a FliG ring model that is compatible with cryo-electron tomography (cryo-ET) and electron microscopy (EM) densities of the C-ring. Collectively, our study contributes to a better mechanistic understanding of the flagellar rotor assembly and its function.

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.

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.

A Benzene-Mapping Approach for Uncovering Cryptic Pockets in Membrane-Bound Proteins

Molecular dynamics (MD) simulations in combination with small organic probes present in the solvent have previously been used as a method to reveal cryptic pockets that may not have been identified in experimental structures. We report such a method implemented within the CHARMM force field using the GROMACS simulation package to effectively explore cryptic pockets on the surfaces of membrane-embedded proteins using benzene as a probe molecule. This method, for which we have made implementation files freely available, relies on modified nonbonded parameters in addition to repulsive potentials between membrane lipids and benzene molecules. The method was tested on part of the outer shell of the dengue virus (DENV), for which research into a safe and effective neutralizing antibody or drug molecule is still ongoing. In particular, the envelope (E) protein, associated with the membrane (M) protein, is a lipid membrane-embedded complex which forms a dimer in the mature viral envelope. Solvent mapping was performed for the full, membrane-embedded EM protein complex and compared with similar calculations performed for the isolated, soluble E protein ectodomain dimer in the solvent. Ectodomain-only simulations with benzene exhibited unfolding effects not observed in the more physiologically relevant membrane-associated systems. A cryptic pocket which has been experimentally shown to bind n-octyl-β-d-glucoside detergent was consistently revealed in all benzene-containing simulations. The addition of benzene also enhanced the flexibility and hydrophobic exposure of cryptic pockets at a key, functional interface in the E protein and revealed a novel, potentially druggable pocket that may be targeted to prevent conformational changes associated with viral entry into the cell.