The fatty acid site is coupled to functional motifs in the SARS-CoV-2 spike protein and modulates spike allosteric behaviour

Allosteric modulation by the fatty acid site in the glycosylated SARS-CoV-2 spike

The spike protein is essential to the SARS-CoV-2 virus life cycle, facilitating virus entry and mediating viral-host membrane fusion. The spike contains a fatty acid (FA) binding site between every two neighbouring receptor-binding domains. This site is coupled to key regions in the protein, but the impact of glycans on these allosteric effects has not been investigated. Using dynamical nonequilibrium molecular dynamics (D-NEMD) simulations, we explore the allosteric effects of the FA site in the fully glycosylated spike of the SARS-CoV-2 ancestral variant. Our results identify the allosteric networks connecting the FA site to functionally important regions in the protein, including the receptor-binding motif, an antigenic supersite in the N-terminal domain, the fusion peptide region, and another allosteric site known to bind heme and biliverdin. The networks identified here highlight the complexity of the allosteric modulation in this protein and reveal a striking and unexpected link between different allosteric sites. Comparison of the FA site connections from D-NEMD in the glycosylated and non-glycosylated spike revealed that glycans do not qualitatively change the internal allosteric pathways but can facilitate the transmission of the structural changes within and between subunits.

A Reflection on the Use of Molecular Simulation to Respond to SARS-CoV-2 Pandemic Threats

Molecular simulations play important roles in understanding the lifecycle of the SARS-CoV-2 virus and contribute to the design and development of antiviral agents and diagnostic tests for COVID. Here, we discuss the insights that such simulations have provided and the challenges involved, focusing on the SARS-CoV-2 main protease (Mpro) and the spike glycoprotein. Mpro is the leading target for antivirals, while the spike glycoprotein is the target for vaccine design. Finally, we reflect on lessons from this pandemic for the simulation community. Data sharing initiatives and collaborations across the international research community contributed to advancing knowledge and should be built on to help in future pandemics and other global challenges such as antimicrobial resistance.

Design and Test of Molecules that Interfere with the Recognition Mechanisms between the SARS-CoV-2 Spike Protein and Its Host Cell Receptors

The disruptive impact of the COVID-19 pandemic has led the scientific community to undertake an unprecedented effort to characterize viral infection mechanisms. Among these, interactions between the viral glycosylated Spike and the human receptors ACE2 and TMPRSS2 are key to allowing virus invasion. Here, we report and test a fully rational methodology to design molecules that are capable of perturbing the interactions between these critical players in SARS-CoV-2 pathogenicity. To this end, we computationally identify substructures on the fully glycosylated Spike protein that are not intramolecularly optimized and are thus prone to being stabilized by forming complexes with ACE2 and TMPRSS2. With the aim of competing with the Spike-mediated cell entry mechanisms, we have engineered the predicted putative interaction regions in the form of peptide mimics that could compete with Spike for interaction with ACE2 and/or TMPRSS2. Experimental models of viral entry demonstrate that the designed molecules are able to interfere with viral entry into ACE2/TMPRSS2 expressing cells, while they have no effects on the entry of control viral particles that do not harbor the Spike protein or on the entry of Spike-presenting viral particles into cells that do not display its receptors on their surface.

Valproic acid use is associated with diminished risk of contracting COVID-19, and diminished disease severity: Epidemiologic and in vitro analysis reveal mechanistic insights

The SARS-CoV-2 pandemic has caused unprecedented worldwide infections from persistent mutant variants with various degrees of infectivity and virulence. The elusiveness of a highly penetrant, worldwide vaccination strategy suggests that the complete eradication of SARS-CoV-2 is unlikely. Even with the advent of new antiviral agents, the disease burden worldwide continues to exceed current preventative and therapeutic strategies. Greater interest has been placed towards the development of affordable,broadly effective antiviral therapeutics. Here, we report that the small branched-chain fatty acid Valproic acid (VPA), approved for maintenance of seizure and bipolar disorder, has a novel anti- coronavirus activity that can be augmented with the addition of a long-chain, polyunsaturated omega-3 fatty acid, Docosahexaenoic acid (DHA). An EMR-based epidemiological study of patients tested for COVID-19 demonstrated a correlation exists between a reduced infection rate in patients treated withVPA of up to 25%, as well as a decreased risk of emergency room visits, hospitalization, ICU admission,and use of mechanical ventilation. In vitro studies have demonstrated that VPA modifies gene expression in MRC5 cells. Interestingly, VPA correlates with the inhibition of several SARS-CoV2 interacting genes and the greater inhibition of alpha-coronavirus HCoV-229E (a "common cold" virus) and SARS-CoV2. The VPA-DHA combination activates pre-existing intracellular antiviral mechanisms normally repressed by coronaviruses. Gene expression profiles demonstrate subtle differences in overall gene expression between VPA-treated and VPA-DHA-treated cells. HCoV-229E infection caused an intensely different response with a marked induction of multiple intracellular inflammatory genes. Changes in gene expression took at least 24 hours to manifest and most likely why prior drug screens failed to identify any antiviral VPA activity despite in silico predictions. This report demonstrates an interaction between HDAC inhibition and the potent activation of cellular antiviral responses. A foundation now exists for a low-cost, highly effective antiviral strategy when supplemented with DHA.

Signal Propagation in the ATPase Domain of Mycobacterium tuberculosis DNA Gyrase from Dynamical-Nonequilibrium Molecular Dynamics Simulations

DNA gyrases catalyze negative supercoiling of DNA, are essential for bacterial DNA replication, transcription, and recombination, and are important antibacterial targets in multiple pathogens, including Mycobacterium tuberculosis, which in 2021 caused >1.5 million deaths worldwide. DNA gyrase is a tetrameric (A2B2) protein formed from two subunit types: gyrase A (GyrA) carries the breakage-reunion active site, whereas gyrase B (GyrB) catalyzes ATP hydrolysis required for energy transduction and DNA translocation. The GyrB ATPase domains dimerize in the presence of ATP to trap the translocated DNA (T-DNA) segment as a first step in strand passage, for which hydrolysis of one of the two ATPs and release of the resulting inorganic phosphate is rate-limiting. Here, dynamical-nonequilibrium molecular dynamics (D-NEMD) simulations of the dimeric 43 kDa N-terminal fragment of M. tuberculosis GyrB show how events at the ATPase site (dissociation/hydrolysis of bound nucleotides) are propagated through communication pathways to other functionally important regions of the GyrB ATPase domain. Specifically, our simulations identify two distinct pathways that respectively connect the GyrB ATPase site to the corynebacteria-specific C-loop, thought to interact with GyrA prior to DNA capture, and to the C-terminus of the GyrB transduction domain, which in turn contacts the C-terminal GyrB topoisomerase-primase (TOPRIM) domain responsible for interactions with GyrA and the centrally bound G-segment DNA. The connection between the ATPase site and the C-loop of dimeric GyrB is consistent with the unusual properties of M. tuberculosis DNA gyrase relative to those from other bacterial species.

Mechanistic insights into ligand dissociation from the SARS-CoV-2 spike glycoprotein

The COVID-19 pandemic, driven by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has spurred an urgent need for effective therapeutic interventions. The spike glycoprotein of the SARS-CoV-2 is crucial for infiltrating host cells, rendering it a key candidate for drug development. By interacting with the human angiotensin-converting enzyme 2 (ACE2) receptor, the spike initiates the infection of SARS-CoV-2. Linoleate is known to bind the spike glycoprotein, subsequently reducing its interaction with ACE2. However, the detailed mechanisms underlying the protein-ligand interaction remain unclear. In this study, we characterized the pathways of ligand dissociation and the conformational changes associated with the spike glycoprotein by using ligand Gaussian accelerated molecular dynamics (LiGaMD). Our simulations resulted in eight complete ligand dissociation trajectories, unveiling two distinct ligand unbinding pathways. The preference between these two pathways depends on the gate distance between two α-helices in the receptor binding domain (RBD) and the position of the N-linked glycan at N343. Our study also highlights the essential contributions of K417, N121 glycan, and N165 glycan in ligand unbinding, which are equally crucial in enhancing spike-ACE2 binding. We suggest that the presence of the ligand influences the motions of these residues and glycans, consequently reducing accessibility for spike-ACE2 binding. These findings enhance our understanding of ligand dissociation from the spike glycoprotein and offer significant implications for drug design strategies in the battle against COVID-19.

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.

Exploring the Fatty Acid Binding Pocket in the SARS-CoV-2 Spike Protein - Confirmed and Potential Ligands

Severe Acute Respiratory syndrome 2 (SARS-CoV-2) is a respiratory virus responsible for coronavirus disease 19 (COVID-19) and the still ongoing and unprecedented global pandemic. The key viral protein for cell infection is the spike glycoprotein, a surface-exposed fusion protein that both recognizes and mediates entry into host cells. Within the spike glycoprotein, a fatty acid binding pocket (FABP) was confirmed, with the crystallization of linoleic acid (LA) occupying a well-defined site. Importantly, when the pocket is occupied by a fatty acid, an inactive conformation is stabilized, and cell recognition is hindered. In this review, we discuss ligands reported so far for this site, correlating their activity predicted through in silico studies with antispike experimental activity, assessed by either binding assays or cell-infection assays. LA was the first confirmed ligand, cocrystallized in a cryo-EM structure of the spike protein, resulting in increased stability of the inactive conformation of the spike protein. The next identified ligand, lifitegrast, was also experimentally confirmed as a ligand with antiviral activity, suggesting the potential for diverse chemical scaffolds to bind this site. Finally, SPC-14 was also confirmed as a ligand, although no inhibition assays were performed. In this review, we identified 20 studies describing small-molecule compounds predicted to bind the pocket in in silico studies and with confirmed binding or in vitro activity, either inhibitory activity against the spike-ACE2 interaction or antiviral activity in cell-based assays. When considering all ligands confirmed with in vitro assays, a good overall occupation of the pocket should be complemented with the ability to make direct interactions, both hydrophilic and hydrophobic, with key amino acid residues defining the pocket surface. Among the active compounds, long flexible carbon chains are recurrent, with retinoids capable of binding the FABP, although bulkier systems are also capable of affecting viral fitness. Compounds able to bind this site with high affinity have the potential to stabilize the inactive conformation of the SARS-CoV-2 spike protein and therefore reduce the virus's ability to infect new cells. Since this pocket is conserved in highly pathogenic human coronaviruses, including MERS-CoV and SARS-CoV, this effect could be exploited for the development of new antiviral agents, with broad-spectrum anticoronavirus activity.

Cooperative and structural relationships of the trimeric Spike with infectivity and antibody escape of the strains Delta (B.1.617.2) and Omicron (BA.2, BA.5, and BQ.1)

Herein, we conducted simulations of trimeric Spike from several SARS-CoV-2 variants of concern (Delta and Omicron sub-variants BA.2, BA.5, and BQ.1) and investigated the mechanisms by which specific mutations confer resistance to neutralizing antibodies. We observed that the mutations primarily affect the cooperation between protein domains within and between protomers. The substitutions K417N and L452R expand hydrogen bonding interactions, reducing their interaction with neutralizing antibodies. By interacting with nearby residues, the K444T and N460K mutations in the SpikeBQ.1 variant potentially reduces solvent exposure, thereby promoting resistance to antibodies. We also examined the impact of D614G, P681R, and P681H substitutions on Spike protein structure that may be related to infectivity. The D614G substitution influences communication between a glycine residue and neighboring domains, affecting the transition between up- and -down RBD states. The P681R mutation, found in the Delta variant, enhances correlations between protein subunits, while the P681H mutation in Omicron sub-variants weakens long-range interactions that may be associated with reduced fusogenicity. Using a multiple linear regression model, we established a connection between inter-protomer communication and loss of sensitivity to neutralizing antibodies. Our findings underscore the importance of structural communication between protein domains and provide insights into potential mechanisms of immune evasion by SARS-CoV-2. Overall, this study deepens our understanding of how specific mutations impact SARS-CoV-2 infectivity and shed light on how the virus evades the immune system.

SARS-CoV-2 spike variants differ in their allosteric responses to linoleic acid

The SARS-CoV-2 spike protein contains a functionally important fatty acid (FA) binding site, which is also found in some other coronaviruses, e.g. SARS-CoV and MERS-CoV. The occupancy of the FA site by linoleic acid (LA) reduces infectivity by 'locking' the spike in a less infectious conformation. Here, we use dynamical-nonequilibrium molecular dynamics (D-NEMD) simulations to compare the allosteric responses of spike variants to LA removal. D-NEMD simulations show that the FA site is coupled to other functional regions of the protein, e.g. the receptor-binding motif (RBM), N-terminal domain (NTD), furin cleavage site, and regions surrounding the fusion peptide. D-NEMD simulations also identify the allosteric networks connecting the FA site to these functional regions. The comparison between the wild-type spike and four variants (Alpha, Delta, Delta plus, and Omicron BA.1) shows that the variants differ significantly in their responses to LA removal. The allosteric connections to the FA site on Alpha are generally similar to those on the wild-type protein, with the exception of the RBM and the S71-R78 region, which show a weaker link to the FA site. In contrast, Omicron is the most different variant, exhibiting significant differences in the RBM, NTD, V622-L629, and furin cleavage site. These differences in the allosteric modulation may be of functional relevance, potentially affecting transmissibility and virulence. Experimental comparison of the effects of LA on SARS-CoV-2 variants, including emerging variants, is warranted.

Dynamical Nonequilibrium Molecular Dynamics Simulations Identify Allosteric Sites and Positions Associated with Drug Resistance in the SARS-CoV-2 Main Protease

The SARS-CoV-2 main protease (Mpro) plays an essential role in the coronavirus lifecycle by catalyzing hydrolysis of the viral polyproteins at specific sites. Mpro is the target of drugs, such as nirmatrelvir, though resistant mutants have emerged that threaten drug efficacy. Despite its importance, questions remain on the mechanism of how Mpro binds its substrates. Here, we apply dynamical nonequilibrium molecular dynamics (D-NEMD) simulations to evaluate structural and dynamical responses of Mpro to the presence and absence of a substrate. The results highlight communication between the Mpro dimer subunits and identify networks, including some far from the active site, that link the active site with a known allosteric inhibition site, or which are associated with nirmatrelvir resistance. They imply that some mutations enable resistance by altering the allosteric behavior of Mpro. More generally, the results show the utility of the D-NEMD technique for identifying functionally relevant allosteric sites and networks including those relevant to resistance.

Metabolic alterations upon SARS-CoV-2 infection and potential therapeutic targets against coronavirus infection

The coronavirus disease 2019 (COVID-19) caused by coronavirus SARS-CoV-2 infection has become a global pandemic due to the high viral transmissibility and pathogenesis, bringing enormous burden to our society. Most patients infected by SARS-CoV-2 are asymptomatic or have mild symptoms. Although only a small proportion of patients progressed to severe COVID-19 with symptoms including acute respiratory distress syndrome (ARDS), disseminated coagulopathy, and cardiovascular disorders, severe COVID-19 is accompanied by high mortality rates with near 7 million deaths. Nowadays, effective therapeutic patterns for severe COVID-19 are still lacking. It has been extensively reported that host metabolism plays essential roles in various physiological processes during virus infection. Many viruses manipulate host metabolism to avoid immunity, facilitate their own replication, or to initiate pathological response. Targeting the interaction between SARS-CoV-2 and host metabolism holds promise for developing therapeutic strategies. In this review, we summarize and discuss recent studies dedicated to uncovering the role of host metabolism during the life cycle of SARS-CoV-2 in aspects of entry, replication, assembly, and pathogenesis with an emphasis on glucose metabolism and lipid metabolism. Microbiota and long COVID-19 are also discussed. Ultimately, we recapitulate metabolism-modulating drugs repurposed for COVID-19 including statins, ASM inhibitors, NSAIDs, Montelukast, omega-3 fatty acids, 2-DG, and metformin.

Conformational Behavior of SARS-Cov-2 Spike Protein Variants: Evolutionary Jumps in Sequence Reverberate in Structural Dynamic Differences

SARS-CoV-2 has evolved rapidly in the first 3 years of pandemic diffusion. The initial evolution of the virus appeared to proceed through big jumps in sequence changes rather than through the stepwise accumulation of point mutations on already established variants. Here, we examine whether this nonlinear mutational process reverberates in variations of the conformational dynamics of the SARS-CoV-2 Spike protein (S-protein), the first point of contact between the virus and the human host. We run extensive microsecond-scale molecular dynamics simulations of seven distinct variants of the protein in their fully glycosylated state and set out to elucidate possible links between the mutational spectrum of the S-protein and the structural dynamics of the respective variant, at global and local levels. The results reveal that mutation-dependent structural and dynamic modulations mostly consist of increased coordinated motions in variants that acquire stability and in an increased internal flexibility in variants that are less stable. Importantly, a limited number of functionally important substructures (the receptor binding domain, in particular) share the same time of movements in all variants, indicating efficient preorganization for functional regions dedicated to host interactions. Our results support a model in which the internal dynamics of the S-proteins from different strains varies in a way that reflects the observed random and non-stepwise jumps in sequence evolution, while conserving the functionally oriented traits of conformational dynamics necessary to support productive interactions with host receptors.

Cryo-EM reveals binding of linoleic acid to SARS-CoV-2 spike glycoprotein, suggesting an antiviral treatment strategy

The COVID-19 pandemic and concomitant lockdowns presented a global health challenge and triggered unprecedented research efforts to elucidate the molecular mechanisms and pathogenicity of SARS-CoV-2. The spike glycoprotein decorating the surface of SARS-CoV-2 virions is a prime target for vaccine development, antibody therapy and serology as it binds the host cell receptor and is central for viral cell entry. The electron cryo-microscopy structure of the spike protein revealed a hydrophobic pocket in the receptor-binding domain that is occupied by an essential fatty acid, linoleic acid (LA). The LA-bound spike protein adopts a non-infectious locked conformation which is more stable than the infectious form and shields important immunogenic epitopes. Here, the impact of LA binding on viral infectivity and replication, and the evolutionary conservation of the pocket in other highly pathogenic coronaviruses, including SARS-CoV-2 variants of concern (VOCs), are reviewed. The importance of LA metabolic products, the eicosanoids, in regulating the human immune response and inflammation is highlighted. Lipid and fatty-acid binding to a hydrophobic pocket in proteins on the virion surface appears to be a broader strategy employed by viruses, including picornaviruses and Zika virus. Ligand binding stabilizes their protein structure and assembly, and downregulates infectivity. In the case of rhinoviruses, this has been exploited to develop small-molecule antiviral drugs that bind to the hydrophobic pocket. The results suggest a COVID-19 antiviral treatment based on the LA-binding pocket.

Progressive membrane-binding mechanism of SARS-CoV-2 variant spike proteins

Membrane recognition by viral spike proteins is critical for infection. Here we show the host cell membrane-binding surfaces of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike variants Alpha, Beta, Gamma, Delta, Epsilon, Kappa, and Omicron as well as SARS-CoV-1 and pangolin and bat relatives. They show increases in membrane binding propensities over time, with all spike head mutations in variants, and particularly BA.1, impacting the protein's affinity to cell membranes. Comparison of hundreds of structures yields a progressive model of membrane docking in which spike protein trimers shift from initial perpendicular stances to increasingly tilted positions that draw viral particles alongside host cell membranes before optionally engaging angiotensin-converting enzyme 2 (ACE2) receptors. This culminates in the assembly of the symmetric fusion apparatus, with enhanced membrane interactions of variants explaining their unique cell fusion capacities and COVID-19 disease transmission rates.

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.

Multifaceted membrane binding head of the SARS-CoV-2 spike protein

The SARS-CoV-2 spike protein presents a surface with enormous membrane binding potential to host tissues and organelles of infected cells. Its exposed trimeric head binds not only the angiotensin-converting enzyme 2 (ACE2), but also host phospholipids which are missing from all existing structures. Hence, the membrane interaction surfaces that mediate viral fusion, entry, assembly and egress remain unclear. Here the spike:membrane docking sites are identified based on membrane optimal docking area (MODA) analysis of 3D structures of spike proteins in closed and open conformations at endocytic and neutral pH levels as well as ligand complexes. This reveals multiple membrane binding sites in the closed spike head that together prefer convex membranes and are modulated by pH, fatty acids and post-translational modifications including glycosylation. The exposure of the various membrane interaction sites adjusts upon domain repositioning within the trimer, allowing formation of intermediate bilayer complexes that lead to the prefusion state while also enabling ACE2 receptor recognition. In contrast, all antibodies that target the spike head would block the membrane docking process that precedes ACE2 recognition. Together this illuminates the engagements of the spike protein with plasma, endocytic, ER or exocytic vesicle membranes that help to drive the cycle of viral infection, and offers novel sites for intervention.

Binding behavior of spike protein and receptor binding domain of the SARS-CoV-2 virus at different environmental conditions

A novel coronavirus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was identified as the cause of the COVID-19 pandemic that originated in China in December 2019. Although extensive research has been performed on SARS-CoV-2, the binding behavior of spike (S) protein and receptor binding domain (RBD) of SARS-CoV-2 at different environmental conditions have yet to be studied. The objective of this study is to investigate the effect of temperature, fatty acids, ions, and protein concentration on the binding behavior and rates of association and dissociation between the S protein and RBD of SARS-CoV-2 and the hydrophobic aminopropylsilane (APS) biosensors using biolayer interferometry (BLI) validated with molecular dynamics simulation. Our results suggest three conditions-high ionic concentration, presence of hydrophobic fatty acids, and low temperature-favor the attachment of S protein and RBD to hydrophobic surfaces. Increasing the temperature within an hour from 0 to 25 °C results in S protein detachment, suggesting that freezing can cause structural changes in the S protein, affecting its binding kinetics at higher temperature. At all the conditions, RBD exhibits lower dissociation capabilities than the full-length S trimer protein, indicating that the separated RBD formed stronger attachment to hydrophobic surfaces compared to when it was included in the S protein.