Developing a Fully Glycosylated Full-Length SARS-CoV-2 Spike Protein Model in a Viral Membrane

Significance of the RBD mutations in the SARS-CoV-2 Omicron: from spike opening to antibody escape and cell attachment

We computationally investigated the role of the Omicron RBD mutations on its structure and interactions with surrounding domains in the spike trimer as well as with ACE2. Our results suggest...

Using Molecular Visualisation Techniques to Explain the Molecular Biology of SARS-CoV-2 Spike Protein Mutations to a General Audience

Since the COVID-19 pandemic started in 2019, the virus responsible for the outbreak-SARS-CoV-2-has continued to evolve. Mutations of the virus' spike protein, the main protein driving infectivity and transmissibility, are especially concerning as they may allow the virus to improve its infectivity, transmissibility, and ability to evade the immune system. Understanding how specific molecular changes can alter the behaviour of a virus is challenging for non-experts, but this information helps us to understand the pandemic we are living through and the public health measures and interventions needed to bring it under control. In response to communication challenges arising from the COVID-19 pandemic, we recently developed an online educational application to explain the molecular biology of SARS-CoV-2 spike protein mutations to the general public. We used visualisation techniques such as 3D modelling and animation, which have been shown to be highly effective teaching tools in molecular biology, allowing the viewer to better understand protein structure, function, and dynamics. We also included interactive elements for users to learn actively by engaging with the digital content, and consequently improve information retention.This chapter presents the methodological and technological framework which we used to create this resource, the 'SARS-CoV-2 Spike Protein Mutation Explorer' (SSPME). It explains how molecular visualisation and 3D modelling software were used to develop accurate models of relevant proteins; how 3D animation software was used to accurately visualise the dynamic molecular processes of SARS-CoV-2 infection, transmission, and antibody evasion; and how game development software was used to compile the 3D models and animations into a comprehensive, informative interactive application on SARS-CoV-2 spike protein mutations. This chapter indicates how cutting-edge visualisation techniques and technologies can be used to improve science communication about complex topics in molecular biology and infection biology to the general public, something that is critical to gaining control of the continuing COVID-19 pandemic.

Hotspot Mutations in SARS-CoV-2

Since its emergence in Wuhan, China, severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has spread very rapidly around the world, resulting in a global pandemic. Though the vaccination process has started, the number of COVID-affected patients is still quite large. Hence, an analysis of hotspot mutations of the different evolving virus strains needs to be carried out. In this regard, multiple sequence alignment of 71,038 SARS-CoV-2 genomes of 98 countries over the period from January 2020 to June 2021 is performed using MAFFT followed by phylogenetic analysis in order to visualize the virus evolution. These steps resulted in the identification of hotspot mutations as deletions and substitutions in the coding regions based on entropy greater than or equal to 0.3, leading to a total of 45 unique hotspot mutations. Moreover, 10,286 Indian sequences are considered from 71,038 global SARS-CoV-2 sequences as a demonstrative example that gives 52 unique hotspot mutations. Furthermore, the evolution of the hotspot mutations along with the mutations in variants of concern is visualized, and their characteristics are discussed as well. Also, for all the non-synonymous substitutions (missense mutations), the functional consequences of amino acid changes in the respective protein structures are calculated using PolyPhen-2 and I-Mutant 2.0. In addition to this, SSIPe is used to report the binding affinity between the receptor-binding domain of Spike protein and human ACE2 protein by considering L452R, T478K, E484Q, and N501Y hotspot mutations in that region.

Differential Interactions between Human ACE2 and Spike RBD of SARS-CoV-2 Variants of Concern

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the current coronavirus disease 2019 (COVID-19) pandemic. It is known that the receptor-binding domain (RBD) of the spike protein of SARS-CoV-2 interacts with the human angiotensin-converting enzyme 2 (ACE2) receptor, initiating the entry of SARS-CoV-2. Since its emergence, a number of SARS-CoV-2 variants have been reported, and the variants that show high infectivity are classified as variants of concern according to the United States Centers for Disease Control and Prevention. In this study, we performed both all-atom steered molecular dynamics (SMD) simulations and microscale thermophoresis (MST) experiments to characterize the binding interactions between ACE2 and RBD of all current variants of concern (Alpha, Beta, Gamma, and Delta) and two variants of interest (Epsilon and Kappa). We report that RBD of the Alpha (N501Y) variant requires the highest amount of force initially to be detached from ACE2 due to the N501Y mutation in addition to the role of N90-glycan, followed by Beta/Gamma (K417N/T, E484 K, and N501Y) or Delta (L452R and T478 K) variants. Among all variants investigated in this work, RBD of the Epsilon (L452R) variant is relatively easily detached from ACE2. Our results from both SMD simulations and MST experiments indicate what makes each variant more contagious in terms of RBD and ACE2 interactions. This study could shed light on developing new drugs to inhibit SARS-CoV-2 entry effectively.

Mutational landscape of SARS-CoV-2 genome in Turkey and impact of mutations on spike protein structure

The Coronavirus Disease 2019 (COVID-19) was declared a pandemic in March 2020 by the World Health Organization (WHO). As of May 25th, 2021 there were 2.059.941 SARS-COV2 genome sequences that have been submitted to the GISAID database, with numerous variations. Here, we aim to analyze the SARS-CoV-2 genome data submitted to the GISAID database from Turkey and to determine the variant and clade distributions by the end of May 2021, in accordance with their appearance timeline. We compared these findings to USA, Europe, and Asia data as well. We have also evaluated the effects of spike protein variations, detected in a group of genome sequences of 13 patients who applied to our clinic, by using 3D modeling algorithms. For this purpose, we analyzed 4607 SARS-CoV-2 genome sequences submitted by different lab centers from Turkey to the GISAID database between March 2020 and May 2021. Described mutations were also introduced in silico to the spike protein structure to analyze their isolated impacts on the protein structure. The most abundant clade was GR followed by G, GH, and GRY and we did not detect any V clade. The most common variant was B.1, followed by B.1.1, and the UK variant, B.1.1.7. Our results clearly show a concordance between the variant distributions, the number of cases, and the timelines of different variant accumulations in Turkey. The 3D simulations indicate an increase in the surface hydrophilicity of the reference spike protein and the detected mutations. There was less surface hydrophilicity increase in the Asp614Gly mutation, which exhibits a more compact conformation around the ACE-2 receptor binding domain region, rendering the structure in a “down” conformation. Our genomic findings can help to model vaccination programs and protein modeling may lead to different approaches for COVID-19 treatment strategies.

SARS-CoV-2 spike binding to ACE2 is stronger and longer ranged due to glycan interaction

Highly detailed steered molecular dynamics simulations are performed on differently glycosylated receptor binding domains of the severe acute respiratory syndrome coronavirus-2 spike protein. The binding strength and the binding range increase with glycosylation. The interaction energy rises very quickly when pulling the proteins apart and only slowly drops at larger distances. We see a catch-slip-type behavior whereby interactions during pulling break and are taken over by new interactions forming. The dominant interaction mode is hydrogen bonds, but Lennard-Jones and electrostatic interactions are relevant as well.

N-glycosylation profiles of the SARS-CoV-2 spike D614G mutant and its ancestral protein characterized by advanced mass spectrometry

N-glycosylation plays an important role in the structure and function of membrane and secreted proteins. The spike protein on the surface of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes COVID-19, is heavily glycosylated and the major target for developing vaccines, therapeutic drugs and diagnostic tests. The first major SARS-CoV-2 variant carries a D614G substitution in the spike (S-D614G) that has been associated with altered conformation, enhanced ACE2 binding, and increased infectivity and transmission. In this report, we used mass spectrometry techniques to characterize and compare the N-glycosylation of the wild type (S-614D) or variant (S-614G) SARS-CoV-2 spike glycoproteins prepared under identical conditions. The data showed that half of the N-glycosylation sequons changed their distribution of glycans in the S-614G variant. The S-614G variant showed a decrease in the relative abundance of complex-type glycans (up to 45%) and an increase in oligomannose glycans (up to 33%) on all altered sequons. These changes led to a reduction in the overall complexity of the total N-glycosylation profile. All the glycosylation sites with altered patterns were in the spike head while the glycosylation of three sites in the stalk remained unchanged between S-614G and S-614D proteins.

Evolutionary and Phenotypic Characterization of Two Spike Mutations in European Lineage 20E of SARS-CoV-2

We have detected two mutations in the spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) at amino acid positions 1163 and 1167 that appeared independently in multiple transmission clusters and different genetic backgrounds. Furthermore, both mutations appeared together in a cluster of 1,627 sequences belonging to clade 20E. This cluster is characterized by 12 additional single nucleotide polymorphisms but no deletions. The available structural information on the S protein in the pre- and postfusion conformations predicts that both mutations confer rigidity, which could potentially decrease viral fitness. Accordingly, we observed reduced infectivity of this spike genotype relative to the ancestral 20E sequence in vitro, and the levels of viral RNA in nasopharyngeal swabs were not significantly higher. Furthermore, the mutations did not impact thermal stability or antibody neutralization by sera from vaccinated individuals but moderately reduce neutralization by convalescent-phase sera from the early stages of the pandemic. Despite multiple successful appearances of the two spike mutations during the first year of SARS-CoV-2 evolution, the genotype with both mutations was displaced upon the expansion of the 20I (Alpha) variant. The midterm fate of the genotype investigated was consistent with the lack of advantage observed in the clinical and experimental data.

Advanced descriptors for long-range noncovalent interactions between SARS-CoV-2 spikes and polymer surfaces

The recent pandemic triggered numerous societal efforts aimed to control and limit the spread of SARS-CoV-2. One of these aspects is related on how the virion interacts with inanimate surfaces, which might be the source of secondary infection. Although recent works address the adsorption of the spike protein on surfaces, there is no information concerning the long-range interactions between spike and surfaces, experimented by the virion when is dispersed in the droplet before its possible adsorption. Some descriptors, namely the interaction potentials per single protein and global potentials, were calculated in this work. These descriptors, evaluated for the closed and open states of the spike protein, are correlated to the long-range noncovalent interactions between the SARS-CoV-2 spikes and polymeric surfaces. They are associated with the surface's affinity towards SARS-CoV-2 dispersed in respiratory droplets or water solutions. Molecular-Dynamics simulations were performed to model the surface of three synthetic polymeric materials: Polypropylene (PP), Polyethylene Terephthalate (PET), and Polylactic Acid (PLA), used in Molecular Mechanics simulations to define the above potentials. The descriptors show a similar trend for the three surfaces, highlighting a greater affinity towards the spikes of PP and PLA over PET. For closed and open structures, the long-range interactions with the surfaces decreased in the following order PP ∼ PLA > PET and PLA > PP > PET, respectively. Thus, PLA and PP interact with the virion quite distant from these surfaces to a greater extent concerning the PET surface, however, the differences among the considered surfaces were small. The global potentials show that the long-range interactions are weak compared to classic binding energy of covalent or ionic bonds. The proposed descriptors are useful most of all for a comparative study aimed at quickly preliminary screening of polymeric surfaces. The obtained results should be validated by more accurate method which will be subject of a subsequent work.

Distant residues modulate conformational opening in SARS-CoV-2 spike protein

The novel coronavirus (SARS-CoV-2) pandemic resulted in the largest public health crisis in recent times. Significant drug design effort against SARS-CoV-2 is focused on the receptor-binding domain (RBD) of the spike protein, although this region is highly prone to mutations causing therapeutic resistance. We applied deep data analysis methods on all-atom molecular dynamics simulations to identify key non-RBD residues that play a crucial role in spike−receptor binding and infection. Because the non-RBD residues are typically conserved across multiple coronaviruses, they can be targeted by broad-spectrum antibodies and drugs to treat infections from new strains that might appear during future epidemics.

Infection by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) involves the attachment of the receptor-binding domain (RBD) of its spike proteins to the ACE2 receptors on the peripheral membrane of host cells. Binding is initiated by a down-to-up conformational change in the spike protein, the change that presents the RBD to the receptor. To date, computational and experimental studies that search for therapeutics have concentrated, for good reason, on the RBD. However, the RBD region is highly prone to mutations, and is therefore a hotspot for drug resistance. In contrast, we here focus on the correlations between the RBD and residues distant to it in the spike protein. This allows for a deeper understanding of the underlying molecular recognition events and prediction of the highest-effect key mutations in distant, allosteric sites, with implications for therapeutics. Also, these sites can appear in emerging mutants with possibly higher transmissibility and virulence, and preidentifying them can give clues for designing pan-coronavirus vaccines against future outbreaks. Our model, based on time-lagged independent component analysis (tICA) and protein graph connectivity network, is able to identify multiple residues that exhibit long-distance coupling with the RBD opening. Residues involved in the most ubiquitous D614G mutation and the A570D mutation of the highly contagious UK SARS-CoV-2 variant are predicted ab initio from our model. Conversely, broad-spectrum therapeutics like drugs and monoclonal antibodies can target these key distant-but-conserved regions of the spike protein.

Structural and Functional Analysis of Female Sex Hormones against SARS-CoV-2 Cell Entry

Emerging evidence suggests that males are more susceptible to severe infection by the SARS-CoV-2 virus than females. A variety of mechanisms may underlie the observed gender-related disparities including differences in sex hormones. However, the precise mechanisms by which female sex hormones may provide protection against SARS-CoV-2 infectivity remains unknown. Here we report new insights into the molecular basis of the interactions between the SARS-CoV-2 spike (S) protein and the human ACE2 receptor. We further report that glycosylation of the ACE2 receptor enhances SARS-CoV-2 infectivity. Importantly, estrogens can disrupt glycan-glycan interactions and glycan-protein interactions between the human ACE2 and the SARS-CoV-2 thereby blocking its entry into cells. In a mouse model of COVID-19, estrogens reduced ACE2 glycosylation and thereby alveolar uptake of the SARS-CoV-2 spike protein. These results shed light on a putative mechanism whereby female sex hormones may provide protection from developing severe infection and could inform the development of future therapies against COVID-19.

S-acylation controls SARS-CoV-2 membrane lipid organization and enhances infectivity

SARS-CoV-2 virions are surrounded by a lipid bilayer that contains membrane proteins such as spike, responsible for target-cell binding and virus fusion. We found that during SARS-CoV-2 infection, spike becomes lipid modified, through the sequential action of the S-acyltransferases ZDHHC20 and 9. Particularly striking is the rapid acylation of spike on 10 cytosolic cysteines within the ER and Golgi. Using a combination of computational, lipidomics, and biochemical approaches, we show that this massive lipidation controls spike biogenesis and degradation, and drives the formation of localized ordered cholesterol and sphingolipid-rich lipid nanodomains in the early Golgi, where viral budding occurs. Finally, S-acylation of spike allows the formation of viruses with enhanced fusion capacity. Our study points toward S-acylating enzymes and lipid biosynthesis enzymes as novel therapeutic anti-viral targets.

Allosteric Control of Structural Mimicry and Mutational Escape in the SARS-CoV-2 Spike Protein Complexes with the ACE2 Decoys and Miniprotein Inhibitors: A Network-Based Approach for Mutational Profiling of Binding and Signaling

We developed a computational framework for comprehensive and rapid mutational scanning of binding energetics and residue interaction networks in the SARS-CoV-2 spike protein complexes. Using this approach, we integrated atomistic simulations and conformational landscaping of the SARS-CoV-2 spike protein complexes with ensemble-based mutational screening and network modeling to characterize mechanisms of structure-functional mimicry and resilience toward mutational escape by the ACE2 protein decoy and de novo designed miniprotein inhibitors. A detailed analysis of structural plasticity of the SARS-CoV-2 spike proteins obtained from atomistic simulations of conformational landscapes and sequence-based profiling of the disorder propensities revealed the intrinsically flexible regions that harbor key functional sites targeted by circulating variants. The conservation of collective dynamics in the SARS-CoV-2 spike protein complexes showed that mutational escape positions are important for modulation of functional motions and that mutational changes in these sites can alter allosteric interaction networks. Through mutational profiling of binding and allosteric propensities in the SARS-CoV-2 spike protein complexes, we identified the key binding and regulatory hotspots that collectively determine functional response and resilience of miniproteins to mutational variants. The results suggest that binding affinities and allosteric signatures of the SARS-CoV-2 complexes can be determined by dynamic crosstalk between structurally stable regulatory centers and conformationally adaptable allosteric hotspots that collectively control the resilience toward mutational escape. This may underlie a mechanism in which moderate perturbations in the mutational escape positions can induce global allosteric changes and alter functional protein response by modulating signaling in the residue interaction networks.

Assembly of Biomolecular Gigastructures and Visualization with the Vulkan Graphics API

Building and displaying all-atom models of biomolecular structures with millions or billions of atoms, like virus particles or cells, remain a challenge due to the sheer size of the data, the required levels of automated building, and the visualization limits of today’s graphics hardware. Based on concepts introduced with the CellPack program, we report new algorithms to create such large-scale models using an intermediate coarse-grained “pet representation” of biomolecules with 1/10th the normal size. Pet atoms are placed such that they optimally trace the surface of the original molecule with just ∼1/50th the original atom number and are joined with covalent bonds. Molecular dynamics simulations of pet molecules allow for efficient packing optimization, as well as the generation of realistic DNA/RNA conformations. This pet world can be expanded back to the all-atom representation to be explored and visualized with full details. Essential for the efficient interactive visualization of gigastructures is the use of multiple levels of detail (LODs), where distant molecules are drawn with a heavily reduced polygon count. We present a grid-based algorithm to create such LODs for all common molecular graphics styles (including ball-and-sticks, ribbons, and cartoons) that do not require monochrome molecules to hide LOD transitions. As a practical application, we built all-atom models of SARS-CoV-2, HIV, and an entire presynaptic bouton with 1 μm diameter and 3.6 billion atoms, using modular building blocks to significantly reduce GPU memory requirements through instancing. We employ the Vulkan graphics API to maximize performance on consumer grade hardware and describe how to use the mmCIF format to efficiently store such giant models. An implementation is available as part of the YASARA molecular modeling and simulation program from www.YASARA.org. The free YASARA View program can be used to explore the presented models, which can be downloaded from www.YASARA.org/petworld, a Creative Commons platform for sharing giant biomolecular structures.

Molecular Modeling of the HR2 and Transmembrane Domains of the SARS-CoV-2 S Protein in the Prefusion State

SARS-CoV-2, the causative agent of COVID-19, remains the focus of research worldwide. SARS-CoV-2 entry into the cell starts with its S protein binding to the angiotensin-converting enzyme-2 (ACE2) expressed on the cell surface. The knowledge of the S protein’s spatial structure is indispensable for understanding the molecular principles of its work. The S protein structure has been almost fully described using experimental approaches with the only exception for the protein’s endodomain, the transmembrane domain, and the ectodomain parts adjacent to the latter. The paper reports molecular modelling of the S protein fragment corresponding to its coiled coil HR2 domain and fully palmitoylated transmembrane domain. Model stability in lipid bilayer was confirmed by all-atom and coarse-grained molecular dynamics simulations. It has been demonstrated that palmitoylation leads to a significant decrease in transmembrane domain mobility and local bilayer thickening, which may be relevant for protein trimerization.

Mutation-induced changes in the receptor-binding interface of the SARS-CoV-2 Delta variant B.1.617.2 and implications for immune evasion

Following the initial surges of the Alpha (B.1.1.7) and the Beta (B.1.351) variants, a more infectious Delta variant (B.1.617.2) is now surging, further deepening the health crises caused by the pandemic. The sharp rise in cases attributed to the Delta variant has made it especially disturbing and is a variant of concern. Fortunately, current vaccines offer protection against known variants of concern, including the Delta variant. However, the Delta variant has exhibited some ability to dodge the immune system as it is found that neutralizing antibodies from prior infections or vaccines are less receptive to binding with the Delta spike protein. Here, we investigated the structural changes caused by the mutations in the Delta variant's receptor-binding interface and explored the effects on binding with the ACE2 receptor as well as with neutralizing antibodies. We find that the receptor-binding β-loop-β motif adopts an altered but stable conformation causing separation in some of the antibody binding epitopes. Our study shows reduced binding of neutralizing antibodies and provides a possible mechanism for the immune evasion exhibited by the Delta variant.

Dynamic Interactions of Fully Glycosylated SARS-CoV-2 Spike Protein with Various Antibodies

The spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) presents a public health crisis, and the vaccines that can induce highly potent neutralizing antibodies are essential for ending the pandemic. The spike (S) protein on the viral envelope mediates human angiotensin-converting enzyme 2 binding and thus is the target of a variety of neutralizing antibodies. In this work, we built various S trimer-antibody complex structures on the basis of the fully glycosylated S protein models described in our previous work and performed all-atom molecular dynamics simulations to gain insight into the structural dynamics and interactions between S protein and antibodies. Investigation of the residues critical for S-antibody binding allows us to predict the potential influence of mutations in SARS-CoV-2 variants. Comparison of the glycan conformations between S-only and S-antibody systems reveals the roles of glycans in S-antibody binding. In addition, we explored the antibody binding modes and the influences of antibody on the motion of S protein receptor binding domains. Overall, our analyses provide a better understanding of S-antibody interactions, and the simulation-based S-antibody interaction maps could be used to predict the influences of S mutation on S-antibody interactions, which will be useful for the development of vaccine and antibody-based therapy.

Identifying Distinct Structural Features of the SARS-CoV-2 Spike Protein Fusion Domain Essential for Membrane Interaction

The SARS-CoV-2 spike protein is the primary antigenic determinant of the virus and has been studied extensively, yet the process of membrane fusion remains poorly understood. The fusion domain (FD) of viral glycoproteins is well established as facilitating the initiation of membrane fusion. An improved understanding of the structural plasticity associated with these highly conserved regions aids in our knowledge of the molecular mechanisms that drive viral fusion. Within the spike protein, the FD of SARS-CoV-2 exists immediately following S2' cleavage at the N-terminus of the S2 domain. Here we have shown that following the introduction of a membrane at pH 7.4, the FD undergoes a transition from a random coil to a more structurally well-defined postfusion state. Furthermore, we have classified the domain into two distinct regions, a fusion peptide (FP, S816-G838) and a fusion loop (FL, D839-F855). The FP forms a helix-turn-helix motif upon association with a membrane, and the favorable entropy gained during this transition from a random coil is likely the driving force behind membrane insertion. Membrane depth experiments then revealed the FP is found inserted within the membrane below the lipid headgroups, while the interaction of the FL with the membrane is shallower in nature. Thus, we propose a structural model relevant to fusion at the plasma membrane in which the FP inserts itself just below the phospholipid headgroups and the FL lays upon the lipid membrane surface.

Atomistic Simulations and In Silico Mutational Profiling of Protein Stability and Binding in the SARS-CoV-2 Spike Protein Complexes with Nanobodies: Molecular Determinants of Mutational Escape Mechanisms

Structure-functional studies have recently revealed a spectrum of diverse high-affinity nanobodies with efficient neutralizing capacity against SARS-CoV-2 virus and resilience against mutational escape. In this study, we combine atomistic simulations with the ensemble-based mutational profiling of binding for the SARS-CoV-2 S-RBD complexes with a wide range of nanobodies to identify dynamic and binding affinity fingerprints and characterize the energetic determinants of nanobody-escaping mutations. Using an in silico mutational profiling approach for probing the protein stability and binding, we examine dynamics and energetics of the SARS-CoV-2 complexes with single nanobodies Nb6 and Nb20, VHH E, a pair combination VHH E + U, a biparatopic nanobody VHH VE, and a combination of the CC12.3 antibody and VHH V/W nanobodies. This study characterizes the binding energy hotspots in the SARS-CoV-2 protein and complexes with nanobodies providing a quantitative analysis of the effects of circulating variants and escaping mutations on binding that is consistent with a broad range of biochemical experiments. The results suggest that mutational escape may be controlled through structurally adaptable binding hotspots in the receptor-accessible binding epitope that are dynamically coupled to the stability centers in the distant binding epitope targeted by VHH U/V/W nanobodies. This study offers a plausible mechanism in which through cooperative dynamic changes, nanobody combinations and biparatopic nanobodies can elicit the increased binding affinity response and yield resilience to common escape mutants.

Local topology and bifurcation hot-spots in proteins with SARS-CoV-2 spike protein as an example

Novel topological methods are introduced to protein research. The aim is to identify hot-spot sites where a bifurcation can alter the local topology of the protein backbone. Since the shape of a protein is intimately related to its biological function, a substitution that causes a bifurcation should have an enhanced capacity to change the protein's function. The methodology applies to any protein but it is developed with the SARS-CoV-2 spike protein as a timely example. First, topological criteria are introduced to identify and classify potential bifurcation hot-spot sites along the protein backbone. Then, the expected outcome of asubstitution, if it occurs, is estimated for a general class of hot-spots, using a comparative analysis of the surrounding backbone segments. The analysis combines the statistics of structurally commensurate amino acid fragments in the Protein Data Bank with general stereochemical considerations. It is observed that the notorious D614G substitution of the spike protein is a good example of a bifurcation hot-spot. A number of topologically similar examples are then analyzed in detail, some of them are even better candidates for a bifurcation hot-spot than D614G. The local topology of the more recently observed N501Y substitution is also inspected, and it is found that this site is proximal to a different kind of local topology changing bifurcation.

Computational Study of Potential Galectin-3 Inhibitors in the Treatment of COVID-19

Galectin-3 is a carbohydrate-binding protein and the most studied member of the galectin family. It regulates several functions throughout the body, among which are inflammation and post-injury remodelling. Recent studies have highlighted the similarity between Galectin-3's carbohydrate recognition domain and the so-called "galectin fold" present on the N-terminal domain of the S1 sub-unit of the SARS-CoV-2 spike protein. Sialic acids binding to the N-terminal domain of the Spike protein are known to be crucial for viral entry into humans, and the role of Galectin-3 as a mediator of lung fibrosis has long been the object of study since its levels have been found to be abnormally high in alveolar macrophages following lung injury. In this context, the discovery of a double inhibitor may both prevent viral entry and reduce post-infection pulmonary fibrosis. In this study, we use a database of 56 compounds, among which 37 have known experimental affinity with Galectin-3. We carry out virtual screening of this database with respect to Galectin-3 and Spike protein. Several ligands are found to exhibit promising binding affinity and interaction with the Spike protein's N-terminal domain as well as with Galectin-3. This finding strongly suggests that existing Galectin-3 inhibitors possess dual-binding capabilities to disrupt Spike-ACE2 interactions. Herein we identify the most promising inhibitors of Galectin-3 and Spike proteins, of which five emerge as potential dual effective inhibitors. Our preliminary results warrant further in vitro and in vivo testing of these putative inhibitors against SARS-CoV-2 with the hope of being able to halt the spread of the virus in the future.

A Novel Multicolour-thresholding Auto-detection Method to Detect the Location and Severity of Inflammation in Confirmed SARS-COV-2 Cases Using Chest X-ray Images

OBJECTIVES

Since late 2019, Coronavirus Disease 2019 (COVID-19) has spread around the world. It has been determined that the disease is very contagious and can cause acute respiratory distress (ARD). Medical imaging has the potential to help identify, detect, and quantify the severity of this infection. This work seeks to develop a novel auto-detection technique for verified COVID-19 cases that can detect aberrant alterations in traditional X-ray pictures.

METHODS

Nineteen separate-colored layers were created from X ray scans of patients diagnosed with COVID-19. Each layer represents objects that have a similar contrast and can be represented by a single color. On a single layer, objects with similar contrasts are formed. A single color image was created by extracting all the objects from all the layers. The prototype model could recognize a wide range of abnormal changes in the image texture based on color differentiation. This was true even when the contrast values of the detected uncleared abnormalities varied a little.

RESULTS

The results indicate that the proposed novel method is 91% accurate in detecting and grading COVID-19 lung infection when compared to the opinions of three experienced radiologists evaluating chest X-ray images. Additionally, the method can be used to determine the infection site and severity of the disease by categorizing the X-rays into five severity levels.

CONCLUSION

By comparing affected tissue to healthy tissue, the proposed COVID-19 auto-detection method can identify locations and indicate the severity of the disease, as well as predict where the disease may spread.

Computational insights into the membrane fusion mechanism of SARS-CoV-2 at the cellular level

The membrane fusion mechanism of SARS-CoV-2 offers an attractive target for the development of small molecule antiviral inhibitors. Fusion involves an initial binding of the crown-like trimeric spike glycoproteins of SARS-CoV-2 to the receptor angiotensin II-converting enzyme 2 (ACE2) on the permissive host cellular membrane and a prefusion to post-fusion conversion of the spike trimer. During this conversion, the fusion peptides of the spike trimer are inserted into the host membrane to bring together the host and viral membranes for membrane fusion in highly choreographic events. However, geometric constraints due to interactions with the membranes remain poorly understood. In this study, we build structural models of super-complexes of spike trimer/ACE2 dimers based on the molecular structures of the ACE2/neutral amino acid transporter B(0)AT heterodimer. We determine the conformational constraints due to the membrane geometry on the enzymatic activity of ACE2 and on the viral fusion process. Furthermore, we find that binding three ACE2 dimers per spike trimer is essential to open the central pore as necessary for triggering productive membrane fusion through an elongation of the central stalk. The reported findings thus provide valuable insights for targeting the membrane fusion mechanism for drug design at the molecular level.

A monoclonal antibody against staphylococcal enterotoxin B superantigen inhibits SARS-CoV-2 entry in vitro

We recently discovered a superantigen-like motif sequentially and structurally similar to a staphylococcal enterotoxin B (SEB) segment, near the S1/S2 cleavage site of the SARS-CoV-2 spike protein, which might explain the multisystem inflammatory syndrome (MIS-C) observed in children and the cytokine storm in severe COVID-19 patients. We show here that an anti-SEB monoclonal antibody (mAb), 6D3, can bind this viral motif at its polybasic (PRRA) insert to inhibit infection in live virus assays. The overlap between the superantigenic site of the spike and its proteolytic cleavage site suggests that the mAb prevents viral entry by interfering with the proteolytic activity of cell proteases (furin and TMPRSS2). The high affinity of 6D3 for this site originates from a polyacidic segment at its heavy chain CDR2. The study points to the potential utility of 6D3 for possibly treating COVID-19, MIS-C, or common colds caused by human coronaviruses that also possess a furin-like cleavage site.

Determination of the interaction between the receptor binding domain of 2019-nCoV spike protein, TMPRSS2, cathepsin B and cathepsin L, and glycosidic and aglycon forms of some flavonols

The novel coronavirus (COVID-19, SARS-CoV-2) is a rapidly spreading disease with a high mortality. In this research, the interactions between specific flavonols and the 2019-nCoV receptor binding domain (RBD), transmembrane protease, serine 2 (TMPRSS2), and cathepsins (CatB and CatL) were analyzed. According to the relative binding capacity index (RBCI) calculated based on the free energy of binding and calculated inhibition constants, it was determined that robinin (ROB) and gossypetin (GOS) were the most effective flavonols on all targets. While the binding free energy of ROB with the spike glycoprotein RBD, TMPRSS2, CatB, and CatL were -5.02, -7.57, -10.10, and -6.11 kcal/mol, the values for GOS were -4.67, -5.24, -8.31, and -6.76, respectively. Furthermore, both compounds maintained their stability for at least 170 ns on respective targets in molecular dynamics simulations. The molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) calculations also corroborated these data. Considering Lipinski's rule of five, ROB and GOS exhibited 3 (MW>500, N or O>10, NH or OH>5), and 1 (NH or OH>5) violations, respectively. Neither ROB nor GOS showed AMES toxicity or hepatotoxicity. The LD50 of these compounds in rats were 2.482 and 2.527 mol/kg, respectively. Therefore, we conclude that these compounds could be considered as alternative therapeutic agents in the treatment of COVID-19. However, the possible inhibitory effects of these compounds on cytochromes (CYPs) should be verified by in vitro or in vivo tests and their adverse effects on cellular energy metabolism should be minimized by performing molecular modifications if necessary.

Peptide Platform as a Powerful Tool in the Fight against COVID-19

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in a global pandemic causing over 195 million infections and more than 4 million fatalities as of July 2021.To date, it has been demonstrated that a number of mutations in the spike glycoprotein (S protein) of SARS-CoV-2 variants of concern abrogate or reduce the neutralization potency of several therapeutic antibodies and vaccine-elicited antibodies. Therefore, the development of additional vaccine platforms with improved supply and logistic profile remains a pressing need. In this work, we have validated the applicability of a peptide-based strategy focused on a preventive as well as a therapeutic purpose. On the basis of the involvement of the dipeptidyl peptidase 4 (DPP4), in addition to the angiotensin converting enzyme 2 (ACE2) receptor in the mechanism of virus entry, we analyzed peptides bearing DPP4 sequences by protein-protein docking and assessed their ability to block pseudovirus infection in vitro. In parallel, we have selected and synthetized peptide sequences located within the highly conserved receptor-binding domain (RBD) of the S protein, and we found that RBD-based vaccines could better promote elicitation of high titers of neutralizing antibodies specific against the regions of interest, as confirmed by immunoinformatic methodologies and in vivo studies. These findings unveil a key antigenic site targeted by broadly neutralizing antibodies and pave the way to the design of pan-coronavirus vaccines.

A critical overview of computational approaches employed for COVID-19 drug discovery

COVID-19 has resulted in huge numbers of infections and deaths worldwide and brought the most severe disruptions to societies and economies since the Great Depression. Massive experimental and computational research effort to understand and characterize the disease and rapidly develop diagnostics, vaccines, and drugs has emerged in response to this devastating pandemic and more than 130 000 COVID-19-related research papers have been published in peer-reviewed journals or deposited in preprint servers. Much of the research effort has focused on the discovery of novel drug candidates or repurposing of existing drugs against COVID-19, and many such projects have been either exclusively computational or computer-aided experimental studies. Herein, we provide an expert overview of the key computational methods and their applications for the discovery of COVID-19 small-molecule therapeutics that have been reported in the research literature. We further outline that, after the first year the COVID-19 pandemic, it appears that drug repurposing has not produced rapid and global solutions. However, several known drugs have been used in the clinic to cure COVID-19 patients, and a few repurposed drugs continue to be considered in clinical trials, along with several novel clinical candidates. We posit that truly impactful computational tools must deliver actionable, experimentally testable hypotheses enabling the discovery of novel drugs and drug combinations, and that open science and rapid sharing of research results are critical to accelerate the development of novel, much needed therapeutics for COVID-19.

What Binds Cationic Photosensitizers Better: Brownian Dynamics Reveals Key Interaction Sites on Spike Proteins of SARS-CoV, MERS-CoV, and SARS-CoV-2

We compared the electrostatic properties of the spike proteins (S-proteins) of three coronaviruses, SARS-CoV, MERS-CoV, and SARS-CoV-2, and their interactions with photosensitizers (PSs), octacationic octakis(cholinyl)zinc phthalocyanine (Zn-PcChol₈₊) and monocationic methylene blue (MB). We found a major common PS binding site at the connection of the S-protein stalk and head. The molecules of Zn-PcChol₈₊ and MB also form electrostatic encounter complexes with large area of negative electrostatic potential at the head of the S-protein of SARS-CoV-2, between fusion protein and heptad repeat 1 domain. The top of the SARS-CoV spike head demonstrates a notable area of electrostatic contacts with Zn-PcChol₈₊ and MB that corresponds to the N-terminal domain. The S-protein protomers of SARS-CoV-2 in “open” and “closed” conformations demonstrate different ability to attract PS molecules. In contrast with Zn-PcChol₈₊, MB possesses the ability to penetrate inside the pocket formed as a result of SARS-CoV-2 receptor binding domain transition into the “open” state. The existence of binding site for cationic PSs common to the S-proteins of SARS-CoV, SARS-CoV-2, and MERS-CoV creates prospects for the wide use of this type of PSs to combat the spread of coronaviruses.

DNA Aptamers Block the Receptor Binding Domain at the Spike Protein of SARS-CoV-2

DNA aptamers are versatile molecular species obtained by the folding of short single-stranded nucleotide sequences, with highly specific recognition capabilities against proteins. Here we test the ability of DNA aptamers to interact with the spike (S-)protein of the SARS-CoV-2 viral capsid. The S-protein, a trimer made up of several subdomains, develops the crucial function of recognizing the ACE2 receptors on the surface of human cells, and subsequent fusioning of the virus membrane with the host cell membrane. In order to achieve this, the S1 domain of one protomer switches between a closed conformation, in which the binding site is inaccessible to the cell receptors, and an open conformation, in which ACE2 can bind, thereby initiating the entry process of the viral genetic material in the host cell. Here we show, by means of state-of-the-art molecular simulations, that small DNA aptamers experimentally identified can recognize the S-protein of SARS-CoV-2, and characterize the details of the binding process. We find that their interaction with different subdomains of the S-protein can effectively block, or at least considerably slow down the opening process of the S1 domain, thereby significantly reducing the probability of virus-cell binding. We provide evidence that, as a consequence, binding of the human ACE2 receptor may be crucially affected under such conditions. Given the facility and low cost of fabrication of specific aptamers, the present findings could open the way to both an innovative viral screening technique with sub-nanomolar sensitivity, and to an effective and low impact curative strategy.

Effect of natural mutations of SARS-CoV-2 on spike structure, conformation, and antigenicity

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants with multiple spike mutations enable increased transmission and antibody resistance. We combined cryo-electron microscopy (cryo-EM), binding, and computational analyses to study variant spikes, including one that was involved in transmission between minks and humans, and others that originated and spread in human populations. All variants showed increased angiotensin-converting enzyme 2 (ACE2) receptor binding and increased propensity for receptor binding domain (RBD)-up states. While adaptation to mink resulted in spike destabilization, the B.1.1.7 (UK) spike balanced stabilizing and destabilizing mutations. A local destabilizing effect of the RBD E484K mutation was implicated in resistance of the B.1.1.28/P.1 (Brazil) and B.1.351 (South Africa) variants to neutralizing antibodies. Our studies revealed allosteric effects of mutations and mechanistic differences that drive either interspecies transmission or escape from antibody neutralization.

Differential Interactions Between Human ACE2 and Spike RBD of SARS-CoV-2 Variants of Concern

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the current coronavirus disease 2019 (COVID-19) pandemic. It is known that the receptor-binding domain (RBD) of the spike protein of SARS-CoV-2 interacts with the human angiotensin-converting enzyme 2 (ACE2) receptor, initiating the entry of SARS-CoV-2. Since its emergence, a number of SARS-CoV-2 variants have been reported, and the variants that show high infectivity are classified as the variants of concern according to the US CDC. In this study, we performed both all-atom steered molecular dynamics (SMD) simulations and microscale thermophoresis (MST) experiments to characterize the binding interactions between ACE2 and RBD of all current variants of concern (Alpha, Beta, Gamma, and Delta) and two variants of interest (Epsilon and Kappa). We report that the RBD of the Alpha (N501Y) variant requires the highest amount of force initially to be detached from ACE2 due to the N501Y mutation in addition to the role of N90-glycan, followed by Beta/Gamma (K417N/T, E484K, and N501Y) or Delta (L452R and T478K) variant. Among all variants investigated in this work, the RBD of the Epsilon (L452R) variant is relatively easily detached from ACE2. Our results combined SMD simulations and MST experiments indicate what makes each variant more contagious in terms of RBD and ACE2 interactions. This study could help develop new drugs to inhibit SARS-CoV-2 entry effectively.

ABSTRACT FIGURE

Identification of Potential Binding Sites of Sialic Acids on the RBD Domain of SARS-CoV-2 Spike Protein

COVID-19, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is still an emergent pandemic for humans. The virus infection is achieved by penetrating its spike protein to host cells via binding with ACE2. Moreover, recent studies show that SARS-CoV-2 may have multiple receptors that need to be further revealed. SARS-CoV-2 shares similar sequences of the spike protein with the Middle East Respiratory Syndrome Coronavirus (MERS-CoV), which can invade host cells by binding to either DPP4 or sialic acids. Sialic acids can be linked to the terminal of glycoproteins and gangliosides are used as one of the receptors of many types of viruses. Therefore, it is very interesting to determine whether sialic acid is a potential receptor of SARS-CoV-2. To address this question, we took N-Acetylneuraminic acid (Neu5Ac), a type of predominant sialic acid found in human cells, as the molecular probe to computationally search the surface of the spike protein to locate the potential binding sites of Neu5Ac. SPR analysis and mass spectrum analysis confirmed the interaction between Neu5Ac and spike protein. This study shows that sialic acids can moderately interact with the spike protein of SARS-CoV-2 by binding between the two RBDs of the spike protein, indicating it could be a potential secondary or auxiliary receptor of SARS-CoV-2.

Enhanced sampling protocol to elucidate fusion peptide opening of SARS-CoV-2 spike protein

Large-scale conformational transitions in the spike protein S2 domain are required during host-cell infection of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus. Although conventional molecular dynamics simulations have been extensively used to study therapeutic targets of SARS-CoV-2, it is still challenging to gain molecular insight into the key conformational changes because of the size of the spike protein and the long timescale required to capture these transitions. In this work, we have developed an efficient simulation protocol that leverages many short simulations, a dynamic selection algorithm, and Markov state models to interrogate the structural changes of the S2 domain. We discovered that the conformational flexibility of the dynamic region upstream of the fusion peptide in S2 is coupled to the proteolytic cleavage state of the spike protein. These results suggest that opening of the fusion peptide likely occurs on a submicrosecond timescale after cleavage at the S2' site. Building on the structural and dynamical information gained to date about S2 domain dynamics, we provide proof of principle that a small molecule bound to a seam neighboring the fusion peptide can slow the opening of the fusion peptide, leading to a new inhibition strategy for experiments to confirm. In aggregate, these results will aid the development of drug cocktails to inhibit infections caused by SARS-CoV-2 and other coronaviruses.

Neuropilin-1 assists SARS-CoV-2 infection by stimulating the separation of Spike protein S1 and S2

The cell surface receptor Neuropilin-1 (Nrp1) was recently identified as a host factor for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) entry. The Spike protein of SARS-CoV-2 is cleaved into two segments, the S1 (residues (res.) 1-685) and the S2 (res. 686-1273) domains by furin protease. Nrp1 predominantly binds to the C-terminal RRAR amino acid motif (res. 682-685) of the S1 domain. In this study, we firstly modeled the association of an Nrp1 protein (consisting of domains a2-b1-b2) with the Spike protein. Next, we studied the separation of S2 from the S1 domain, with and without Nrp1 bound, by utilizing molecular dynamics pulling simulations. During the separation, Nrp1 stabilizes the S1 C-terminal region (res. 640-685) and thereby assists the detachment of S2 N-terminal region (res. 686-700). Without Nrp1 bound, S1 tends to become stretched, whereas the bound Nrp1 stimulates an earlier separation of S2 from the S1 domain. The liberated S2 domain is known to mediate the fusion of virus and host membranes; thus, Nrp1 likely increases virus infectivity by facilitating the S1 and S2 separation. We further analyzed the possible topological structure of the SARS-CoV-2 Spike protein when bound with Nrp1 and angiotensin-converting enzyme 2 (ACE2). Understanding of such an Nrp1-assisted viral infection opens the gate for the generation of protein-protein inhibitors, such as antibodies, which could attenuate the infection mechanism and protect certain cells in a future Nrp1-ACE2 targeted combination therapy.

Structural and Dynamical Differences in the Spike Protein RBD in the SARS-CoV-2 Variants B.1.1.7 and B.1.351

The novel coronavirus (SARS-CoV-2) pandemic that started in late 2019 is responsible for hundreds of millions of cases worldwide and millions of fatalities. Though vaccines are available, the virus is mutating to form new strains among which are the variants B.1.1.7 and B.1.351 that demonstrate increased transmissivity and infectivity. In this study, we performed molecular dynamics simulations to explore the role of the mutations in the interaction of the virus spike protein receptor binding domain (RBD) with the host receptor ACE2. We find that the hydrogen bond networks are rearranged in the variants and also that new hydrogen bonds are established between the RBD and ACE2 as a result of mutations. We investigated three variants: the wild-type (WT), B.1.1.7, and B.1.351. We find that the B.1.351 variant (also known as 501Y.V2) shows larger flexibility in the RBD loop segment involving residue K484, yet the RBD-ACE2 complex shows higher stability. Mutations that allow a more flexible interface that can result in a more stable complex may be a factor contributing to the increased infectivity of the mutated variants.

On-Site Detection of SARS-CoV-2 Antigen by Deep Learning-Based Surface-Enhanced Raman Spectroscopy and Its Biochemical Foundations

A rapid, on-site, and accurate SARS-CoV-2 detection method is crucial for the prevention and control of the COVID-19 epidemic. However, such an ideal screening technology has not yet been developed for the diagnosis of SARS-CoV-2. Here, we have developed a deep learning-based surface-enhanced Raman spectroscopy technique for the sensitive, rapid, and on-site detection of the SARS-CoV-2 antigen in the throat swabs or sputum from 30 confirmed COVID-19 patients. A Raman database based on the spike protein of SARS-CoV-2 was established from experiments and theoretical calculations. The corresponding biochemical foundation for this method is also discussed. The deep learning model could predict the SARS-CoV-2 antigen with an identification accuracy of 87.7%. These results suggested that this method has great potential for the diagnosis, monitoring, and control of SARS-CoV-2 worldwide.

Landscape-Based Mutational Sensitivity Cartography and Network Community Analysis of the SARS-CoV-2 Spike Protein Structures: Quantifying Functional Effects of the Circulating D614G Variant

We developed and applied a computational approach to simulate functional effects of the global circulating mutation D614G of the SARS-CoV-2 spike protein. All-atom molecular dynamics simulations are combined with deep mutational scanning and analysis of the residue interaction networks to investigate conformational landscapes and energetics of the SARS-CoV-2 spike proteins in different functional states of the D614G mutant. The results of conformational dynamics and analysis of collective motions demonstrated that the D614 site plays a key regulatory role in governing functional transitions between open and closed states. Using mutational scanning and sensitivity analysis of protein residues, we identified the stability hotspots in the SARS-CoV-2 spike structures of the mutant trimers. The results suggest that the D614G mutation can induce the increased stability of the open form acting as a driver of conformational changes, which may result in the increased exposure to the host receptor and promote infectivity of the virus. The network community analysis of the SARS-CoV-2 spike proteins showed that the D614G mutation can enhance long-range couplings between domains and strengthen the interdomain interactions in the open form, supporting the reduced shedding mechanism. This study provides the landscape-based perspective and atomistic view of the allosteric interactions and stability hotspots in the SARS-CoV-2 spike proteins, offering a useful insight into the molecular mechanisms underpinning functional effects of the global circulating mutations.

Dynamic Profiling of Binding and Allosteric Propensities of the SARS-CoV-2 Spike Protein with Different Classes of Antibodies: Mutational and Perturbation-Based Scanning Reveals the Allosteric Duality of Functionally Adaptable Hotspots

The functional adaptability and conformational plasticity of SARS-CoV-2 spike proteins allow for the efficient modulation of complex phenotypic responses to the host receptor and antibodies. In this study, we combined atomistic simulations with mutational and perturbation-based scanning approaches to examine binding mechanisms of the SARS-CoV-2 spike proteins with three different classes of antibodies. The ensemble-based profiling of binding and allosteric propensities of the SARS-CoV-2 spike protein residues showed that these proteins can work as functionally adaptable and allosterically regulated machines. Conformational dynamics analysis revealed that binding-induced modulation of soft modes can elicit the unique protein response to different classes of antibodies. Mutational scanning heatmaps and sensitivity analysis revealed the binding energy hotspots for different classes of antibodies that are consistent with the experimental deep mutagenesis, showing that differences in the binding affinity caused by global circulating variants in spike positions K417, E484, and N501 are relatively moderate and may not fully account for the observed antibody resistance effects. Through functional dynamics analysis and perturbation-response scanning of the SARS-CoV-2 spike protein residues in the unbound form and antibody-bound forms, we examine how antibody binding can modulate allosteric propensities of spike protein residues and determine allosteric hotspots that control signal transmission and global conformational changes. These results show that residues K417, E484, and N501 targeted by circulating mutations correspond to a group of versatile allosteric centers in which small perturbations can modulate collective motions, alter the global allosteric response, and elicit binding resistance. We suggest that the SARS-CoV-2 S protein may exploit the plasticity of specific allosteric hotspots to generate escape mutants that alter the response to antibody binding without compromising the activity of the spike protein.

Proteomics-Based Insights Into the SARS-CoV-2-Mediated COVID-19 Pandemic: A Review of the First Year of Research

In late 2019, a virus subsequently named severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged in China and led to a worldwide pandemic of the disease termed coronavirus disease 2019. The global health threat posed by this pandemic led to an extremely rapid and robust mobilization of the scientific and medical communities as evidenced by the publication of more than 10,000 peer-reviewed articles and thousands of preprints in the first year of the pandemic alone. With the publication of the initial genome sequence of SARS-CoV-2, the proteomics community immediately joined this effort publishing, to date, more than 100 peer-reviewed proteomics studies and submitting many more preprints to preprint servers. In this review, we focus on peer-reviewed articles published on the proteome, glycoproteome, and glycome of SARS-CoV-2. At a basic level, proteomic studies provide valuable information on quantitative aspects of viral infection course; information on the identities, sites, and microheterogeneity of post-translational modifications; and, information on protein-protein interactions. At a biological systems level, these studies elucidate host cell and tissue responses, characterize antibodies and other immune system factors in infection, suggest biomarkers that may be useful for diagnosis and disease-course monitoring, and help in the development or repurposing of potential therapeutics. Here, we summarize results from selected early studies to provide a perspective on the current rapidly evolving literature.

Determination of intermediate state structures in the opening pathway of SARS-CoV-2 spike using cryo-electron microscopy

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the cause of COVID-19, a highly infectious disease that is severely affecting our society and welfare systems. In order to develop therapeutic interventions against this condition, one promising strategy is to target spike, the trimeric transmembrane glycoprotein that the virus uses to recognise and bind its host cells. Here we use a metainference cryo-electron microscopy approach to determine the opening pathway that brings spike from its inactive (closed) conformation to its active (open) one. The knowledge of the structures of the intermediate states of spike along this opening pathway enables us to identify a cryptic pocket that is not exposed in the open and closed states. These results underline the opportunities offered by the determination of the structures of the transient intermediate states populated during the dynamics of proteins to allow the therapeutic targeting of otherwise invisible cryptic binding sites.

A structural ensemble derived from cryo-electron microscopy reveals a cryptic pocket site in intermediate states along the opening pathway of the SARS-CoV-2 spike protein.

Comparative Perturbation-Based Modeling of the SARS-CoV-2 Spike Protein Binding with Host Receptor and Neutralizing Antibodies: Structurally Adaptable Allosteric Communication Hotspots Define Spike Sites Targeted by Global Circulating Mutations

In this study, we used an integrative computational approach to examine molecular mechanisms and determine functional signatures underlying the role of functional residues in the SARS-CoV-2 spike protein that are targeted by novel mutational variants and antibody-escaping mutations. Atomistic simulations and functional dynamics analysis are combined with alanine scanning and mutational sensitivity profiling of the SARS-CoV-2 spike protein complexes with the ACE2 host receptor and the REGN-COV2 antibody cocktail(REG10987+REG10933). Using alanine scanning and mutational sensitivity analysis, we have shown that K417, E484, and N501 residues correspond to key interacting centers with a significant degree of structural and energetic plasticity that allow mutants in these positions to afford the improved binding affinity with ACE2. Through perturbation-based network modeling and community analysis of the SARS-CoV-2 spike protein complexes with ACE2, we demonstrate that E406, N439, K417, and N501 residues serve as effector centers of allosteric interactions and anchor major intermolecular communities that mediate long-range communication in the complexes. The results provide support to a model according to which mutational variants and antibody-escaping mutations constrained by the requirements for host receptor binding and preservation of stability may preferentially select structurally plastic and energetically adaptable allosteric centers to differentially modulate collective motions and allosteric interactions in the complexes with the ACE2 enzyme and REGN-COV2 antibody combination. This study suggests that the SARS-CoV-2 spike protein may function as a versatile and functionally adaptable allosteric machine that exploits the plasticity of allosteric regulatory centers to fine-tune response to antibody binding without compromising the activity of the spike protein.

Integrated Biophysical Modeling of the SARS-CoV-2 Spike Protein Binding and Allosteric Interactions with Antibodies

Structural and biochemical studies of the severe acute respiratory syndrome (SARS)-CoV-2 spike glycoproteins and complexes with highly potent antibodies have revealed multiple conformation-dependent epitopes highlighting conformational plasticity of spike proteins and capacity for eliciting specific binding and broad neutralization responses. In this study, we used coevolutionary analysis, molecular simulations, and perturbation-based hierarchical network modeling of the SARS-CoV-2 spike protein complexes with a panel of antibodies targeting distinct epitopes to explore molecular mechanisms underlying binding-induced modulation of dynamics and allosteric signaling in the spike proteins. Through coevolutionary analysis of the SARS-CoV-2 spike proteins, we identified highly coevolving hotspots and functional clusters that enable a functional cross-talk between distant allosteric regions in the SARS-CoV-2 spike complexes with antibodies. Coarse-grained and all-atom molecular dynamics simulations combined with mutational sensitivity mapping and perturbation-based profiling of the SARS-CoV-2 receptor-binding domain (RBD) complexes with CR3022 and CB6 antibodies enabled a detailed validation of the proposed approach and an extensive quantitative comparison with the experimental structural and deep mutagenesis scanning data. By combining in silico mutational scanning, perturbation-based modeling, and network analysis of the SARS-CoV-2 spike trimer complexes with H014, S309, S2M11, and S2E12 antibodies, we demonstrated that antibodies can incur specific and functionally relevant changes by modulating allosteric propensities and collective dynamics of the SARS-CoV-2 spike proteins. The results provide a novel insight into regulatory mechanisms of SARS-CoV-2 S proteins showing that antibody-escaping mutations can preferentially target structurally adaptable energy hotspots and allosteric effector centers that control functional movements and allosteric communication in the complexes.

Energy Landscape of the SARS-CoV-2 Reveals Extensive Conformational Heterogeneity

Cryo-electron microscopy (cryo-EM) has produced a number of structural models of the SARS-CoV-2 spike, already prompting biomedical outcomes. However, these reported models and their associated electrostatic potential maps represent an unknown admixture of conformations stemming from the underlying energy landscape of the spike protein. As for any protein, some of the spike's conformational motions are expected to be biophysically relevant, but cannot be interpreted only by static models. Using experimental cryo-EM images, we present the energy landscape of the spike protein conformations, and identify molecular rearrangements along the most-likely conformational path in the vicinity of the open (so called 1RBD-up) state. The resulting global and local atomic refinements reveal larger movements than those expected by comparing the reported 1RBD-up and 1RBD-down cryo-EM models. Here we report greater degrees of "openness" in global conformations of the 1RBD-up state, not revealed in the single-model interpretations of the density maps, together with conformations that overlap with the reported models. We discover how the glycan shield contributes to the stability of these conformations along the minimum free-energy pathway. A local analysis of seven key binding pockets reveals that six out them, including those for engaging ACE2, therapeutic mini-proteins, linoleic acid, two different kinds of antibodies, and protein-glycan interaction sites, switch conformations between their known apo- and holo-conformations, even when the global spike conformation is 1RBD-up. This is reminiscent of a conformational pre-equilibrium. We found only one binding pocket, namely antibody AB-C135 to remain closed along the entire minimum free energy path, suggesting an induced fit mechanism for this enzyme.

Dynamic Interactions of Fully Glycosylated SARS-CoV-2 Spike Protein with Various Antibodies

The spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) presents a public health crisis, and the vaccines that can induce highly potent neutralizing antibodies are essential for ending the pandemic. The spike (S) protein on the viral envelope mediates human angiotensin-converting enzyme 2 (ACE2) binding and thus is the target of a variety of neutralizing antibodies. In this work, we built various S trimer-antibody complex structures on the basis of the fully glycosylated S protein models described in our previous work, and performed all-atom molecular dynamics simulations to get insight into the structural dynamics and interactions between S protein and antibodies. Investigation of the residues critical for S-antibody binding allows us to predict the potential influence of mutations in SARS-CoV-2 variants. Comparison of the glycan conformations between S-only and S-antibody systems reveals the roles of glycans in S-antibody binding. In addition, we explored the antibody binding modes, and the influences of antibody on the motion of S protein receptor binding domains. Overall, our analyses provide a better understanding of S-antibody interactions, and the simulation-based S-antibody interaction maps could be used to predict the influences of S mutation on S-antibody interactions, which will be useful for the development of vaccine and antibody-based therapy.

Fighting viruses with computers, right now

The synergistic conjunction of various technological revolutions with the accumulated knowledge and workflows is rapidly transforming several scientific fields. Particularly, Virology can now feed from accurate physical models, polished computational tools, and massive computational power to readily integrate high-resolution structures into biological representations of unprecedented detail. That preparedness allows for the first time to get crucial information for vaccine and drug design from in-silico experiments against emerging pathogens of worldwide concern at relevant action windows. The present work reviews some of the main milestones leading to these breakthroughs in Computational Virology, providing an outlook for future developments in capacity building and accessibility to computational resources.

Mass spectrometry-based proteomic platforms for better understanding of SARS-CoV-2 induced pathogenesis and potential diagnostic approaches

While protein–protein interaction is the first step of the SARS‐CoV‐2 infection, recent comparative proteomic profiling enabled the identification of over 11,000 protein dynamics, thus providing a comprehensive reflection of the molecular mechanisms underlying the cellular system in response to viral infection. Here we summarize and rationalize the results obtained by various mass spectrometry (MS)‐based proteomic approaches applied to the functional characterization of proteins and pathways associated with SARS‐CoV‐2‐mediated infections in humans. Comparative analysis of cell‐lines versus tissue samples indicates that our knowledge in proteome profile alternation in response to SARS‐CoV‐2 infection is still incomplete and the tissue‐specific response to SARS‐CoV‐2 infection can probably not be recapitulated efficiently by in vitro experiments. However, regardless of the viral infection period, sample types, and experimental strategies, a thorough cross‐comparison of the recently published proteome, phosphoproteome, and interactome datasets led to the identification of a common set of proteins and kinases associated with PI3K‐Akt, EGFR, MAPK, Rap1, and AMPK signaling pathways. Ephrin receptor A2 (EPHA2) was identified by 11 studies including all proteomic platforms, suggesting it as a potential future target for SARS‐CoV‐2 infection mechanisms and the development of new therapeutic strategies. We further discuss the potentials of future proteomics strategies for identifying prognostic SARS‐CoV‐2 responsive age‐, gender‐dependent, tissue‐specific protein targets.

Evidence of a putative glycosaminoglycan binding site on the glycosylated SARS-CoV-2 spike protein N-terminal domain

SARS-CoV-2 has rapidly spread throughout the world's population since its initial discovery in 2019. The virus infects cells via a glycosylated spike protein located on its surface. The protein primarily binds to the angiotensin-converting enzyme-2 (ACE2) receptor, using glycosaminoglycans (GAGs) as co-receptors. Here, we performed bioinformatics and molecular dynamics simulations of the spike protein to investigate the existence of additional GAG binding sites on the receptor-binding domain (RBD), separate from previously reported heparin-binding sites. A putative GAG binding site in the N-terminal domain (NTD) of the protein was identified, encompassing residues 245-246. We hypothesized that GAGs of a sufficient length might bridge the gap between this site and the PRRARS furin cleavage site, including the mutation S247R. Docking studies using GlycoTorch Vina and subsequent MD simulations of the spike trimer in the presence of dodecasaccharides of the GAGs heparin and heparan sulfate supported this possibility. The heparan sulfate chain bridged the gap, binding the furin cleavage site and S247R. In contrast, the heparin chain bound the furin cleavage site and surrounding glycosylation structures, but not S247R. These findings identify a site in the spike protein that favors heparan sulfate binding that may be particularly pertinent for a better understanding of the recent UK and South African strains. This will also assist in future targeted therapy programs that could include repurposing clinical heparan sulfate mimetics.

Structural biology in the time of COVID-19: perspectives on methods and milestones

The global COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has wreaked unprecedented havoc on global society, in terms of a huge loss of life and burden of morbidity, economic upheaval and social disruption. Yet the sheer magnitude and uniqueness of this event has also spawned a massive mobilization of effort in the scientific community to investigate the virus, to develop therapeutics and vaccines, and to understand the public health impacts. Structural biology has been at the center of these efforts, and so it is advantageous to take an opportunity to reflect on the status of structural science vis-à-vis its role in the fight against COVID-19, to register the unprecedented response and to contemplate the role of structural biology in addressing future outbreak threats. As the one-year anniversary of the World Health Organization declaration that COVID-19 is a pandemic has just passed, over 1000 structures of SARS-CoV-2 biomolecules have been deposited in the Worldwide Protein Data Bank (PDB). It is rare to obtain a snapshot of such intense effort in the structural biology arena and is of special interest as the 50th anniversary of the PDB is celebrated in 2021. It is additionally timely as it overlaps with a period that has been termed the 'resolution revolution' in cryoelectron microscopy (CryoEM). CryoEM has recently become capable of producing biomolecular structures at similar resolutions to those traditionally associated with macromolecular X-ray crystallo-graphy. Examining SARS-CoV-2 protein structures that have been deposited in the PDB since the virus was first identified allows a unique window into the power of structural biology and a snapshot of the advantages of the different techniques available, as well as insight into the complementarity of the structural methods.