Modelling conformational state dynamics and its role on infection for SARS-CoV-2 Spike protein variants

Understanding and Predicting Ligand Efficacy in the μ-Opioid Receptor through Quantitative Dynamical Analysis of Complex Structures

The μ-opioid receptor (MOR) is a G-protein coupled receptor involved in nociception and the primary target of opioid drugs. Understanding the relationships among the ligand structure, receptor dynamics, and efficacy in activating MOR is crucial for drug discovery and development. Here, we use coarse-grained normal-mode analysis to predict ligand-induced changes in receptor dynamics with the Quantitative Dynamics Activity Relationship (QDAR) DynaSig-ML methodology, training a LASSO regression model on the entropic signatures (ESs) computed from ligand-receptor complexes. We train and validate the methodology using a data set of 179 MOR ligands with experimentally measured efficacies split into strictly chemically different cross-validation sets. By analyzing the coefficients of the ES LASSO model, we identified key residues involved in MOR activation, several of which have mutational data supporting their role in MOR activation. Additionally, we explored a contact-only LASSO model based on ligand-protein interactions. While the model showed predictive power, it failed at predicting efficacy for ligands with low structural similarity to the training set, emphasizing the importance of receptor dynamics for predicting ligand-induced receptor activation. Moreover, the low computational cost of our approach, at 3 CPU s per ligand-receptor complex, opens the door to its application in large-scale virtual screening contexts. Our work contributes to a better understanding of dynamics-function relationships in the μ-opioid receptor and provides a framework for predicting ligand efficacy based on ligand-induced changes in receptor dynamics.

A Mode Evolution Metric to Extract Reaction Coordinates for Biomolecular Conformational Transitions

The complex, multidimensional energy landscape of biomolecules makes the extraction of suitable, nonintuitive collective variables (CVs) that describe their conformational transitions challenging. At present, dimensionality reduction approaches and machine learning (ML) schemes are employed to obtain CVs from molecular dynamics (MD)/Monte Carlo (MC) trajectories or structural databanks for biomolecules. However, minimum sampling conditions to generate reliable CVs that accurately describe the underlying energy landscape remain unclear. Here, we address this issue by developing a Mode evolution Metric (MeM) to extract CVs that can pinpoint new states and describe local transitions in the vicinity of a reference minimum from nonequilibrated MD/MC trajectories. We present a general mathematical formulation of MeM for both statistical dimensionality reduction and machine learning approaches. Application of MeM to MC trajectories of model potential energy landscapes and MD trajectories of solvated alanine dipeptide reveals that the principal components which locate new states in the vicinity of a reference minimum emerge well before the trajectories locally equilibrate between the associated states. Finally, we demonstrate a possible application of MeM in designing efficient biased sampling schemes to construct accurate energy landscape slices that link transitions between states. MeM can help speed up the search for new minima around a biomolecular conformational state and enable the accurate estimation of thermodynamics for states lying on the energy landscape and the description of associated transitions.

COVID-19 mutations: An overview

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) belongs to the genus Beta coronavirus and the family of Coronaviridae. It is a positive-sense, non-segmented single-strand RNA virus. Four common types of human coronaviruses circulate globally, particularly in the fall and winter seasons. They are responsible for 10%-30% of all mild upper respiratory tract infections in adults. These are 229E, NL63 of the Alfacoronaviridae family, OC43, and HKU1 of the Betacoronaviridae family. However, there are three highly pathogenic human coronaviruses: SARS-CoV-2, Middle East respiratory syndrome coronavirus, and the latest pandemic caused by the SARS-CoV-2 infection. All viruses, including SARS-CoV-2, have the inherent tendency to evolve. SARS-CoV-2 is still evolving in humans. Additionally, due to the development of herd immunity, prior infection, use of medication, vaccination, and antibodies, the viruses are facing immune pressure. During the replication process and due to immune pressure, the virus may undergo mutations. Several SARS-CoV-2 variants, including the variants of concern (VOCs), such as B.1.1.7 (Alpha), B.1.351 (Beta), B.1.617/B.1.617.2 (Delta), P.1 (Gamma), and B.1.1.529 (Omicron) have been reported from various parts of the world. These VOCs contain several important mutations; some of them are on the spike proteins. These mutations may lead to enhanced infectivity, transmissibility, and decreased neutralization efficacy by monoclonal antibodies, convalescent sera, or vaccines. Mutations may also lead to a failure of detection by molecular diagnostic tests, leading to a delayed diagnosis, increased community spread, and delayed treatment. We searched PubMed, EMBASE, Covariant, the Stanford variant Database, and the CINAHL from December 2019 to February 2023 using the following search terms: VOC, SARS-CoV-2, Omicron, mutations in SARS-CoV-2, etc. This review discusses the various mutations and their impact on infectivity, transmissibility, and neutralization efficacy.

Nanomechanical collective vibration of SARS-CoV-2 spike proteins

The development of effective therapeutics against COVID-19 requires a thorough understanding of the receptor recognition mechanism of the SARS-CoV-2 spike (S) protein. Here the multidomain collective dynamics on the trimer of the spike protein has been analyzed using normal mode analysis (NMA). A common nanomechanical profile was identified in the spike proteins of SARS-CoV-2 and its variants. The profile involves collective vibrations of the receptor-binding domain (RBD) and the N-terminal domain (NTD), which may mediate the physical interaction process. Quantitative analysis of the collective modes suggests a nanomechanical property involving large-scale conformational changes, which explains the difference in receptor binding affinity among different variants. These results support the use of intrinsic global dynamics as a valuable perspective for studying the allosteric and functional mechanisms of the S protein. This approach also provides a low-cost theoretical toolkit for screening potential pathogenic mutations and drug targets.

Some mechanistic underpinnings of molecular adaptations of SARS-COV-2 spike protein by integrating candidate adaptive polymorphisms with protein dynamics

We integrate evolutionary predictions based on the neutral theory of molecular evolution with protein dynamics to generate mechanistic insight into the molecular adaptations of the SARS-COV-2 spike (S) protein. With this approach, we first identified candidate adaptive polymorphisms (CAPs) of the SARS-CoV-2 S protein and assessed the impact of these CAPs through dynamics analysis. Not only have we found that CAPs frequently overlap with well-known functional sites, but also, using several different dynamics-based metrics, we reveal the critical allosteric interplay between SARS-CoV-2 CAPs and the S protein binding sites with the human ACE2 (hACE2) protein. CAPs interact far differently with the hACE2 binding site residues in the open conformation of the S protein compared to the closed form. In particular, the CAP sites control the dynamics of binding residues in the open state, suggesting an allosteric control of hACE2 binding. We also explored the characteristic mutations of different SARS-CoV-2 strains to find dynamic hallmarks and potential effects of future mutations. Our analyses reveal that Delta strain-specific variants have non-additive (i.e., epistatic) interactions with CAP sites, whereas the less pathogenic Omicron strains have mostly additive mutations. Finally, our dynamics-based analysis suggests that the novel mutations observed in the Omicron strain epistatically interact with the CAP sites to help escape antibody binding.

Some mechanistic underpinnings of molecular adaptations of SARS-COV-2 spike protein by integrating candidate adaptive polymorphisms with protein dynamics

We integrate evolutionary predictions based on the neutral theory of molecular evolution with protein dynamics to generate mechanistic insight into the molecular adaptations of the SARS-COV-2 Spike (S) protein. With this approach, we first identified Candidate Adaptive Polymorphisms (CAPs) of the SARS-CoV-2 Spike protein and assessed the impact of these CAPs through dynamics analysis. Not only have we found that CAPs frequently overlap with well-known functional sites, but also, using several different dynamics-based metrics, we reveal the critical allosteric interplay between SARS-CoV-2 CAPs and the S protein binding sites with the human ACE2 (hACE2) protein. CAPs interact far differently with the hACE2 binding site residues in the open conformation of the S protein compared to the closed form. In particular, the CAP sites control the dynamics of binding residues in the open state, suggesting an allosteric control of hACE2 binding. We also explored the characteristic mutations of different SARS-CoV-2 strains to find dynamic hallmarks and potential effects of future mutations. Our analyses reveal that Delta strain-specific variants have non-additive (i.e., epistatic) interactions with CAP sites, whereas the less pathogenic Omicron strains have mostly additive mutations. Finally, our dynamics-based analysis suggests that the novel mutations observed in the Omicron strain epistatically interact with the CAP sites to help escape antibody binding.

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.

Raman Fingerprints of SARS-CoV-2 Omicron Subvariants: Molecular Roots of Virological Characteristics and Evolutionary Directions

The latest RNA genomic mutation of SARS-CoV-2 virus, termed the Omicron variant, has generated a stream of highly contagious and antibody-resistant strains, which in turn led to classifying Omicron as a variant of concern. We systematically collected Raman spectra from six Omicron subvariants available in Japan (i.e., BA.1.18, BA.2, BA.4, BA.5, XE, and BA.2.75) and applied machine-learning algorithms to decrypt their structural characteristics at the molecular scale. Unique Raman fingerprints of sulfur-containing amino acid rotamers, RNA purines and pyrimidines, tyrosine phenol ring configurations, and secondary protein structures clearly differentiated the six Omicron subvariants. These spectral characteristics, which were linked to infectiousness, transmissibility, and propensity for immune evasion, revealed evolutionary motifs to be compared with the outputs of genomic studies. The availability of a Raman "metabolomic snapshot", which was then translated into a barcode to enable a prompt subvariant identification, opened the way to rationalize in real-time SARS-CoV-2 activity and variability. As a proof of concept, we applied the Raman barcode procedure to a nasal swab sample retrieved from a SARS-CoV-2 patient and identified its Omicron subvariant by coupling a commercially available magnetic bead technology with our newly developed Raman analyses.

Protein conformational ensembles in function: roles and mechanisms

The sequence-structure-function paradigm has dominated twentieth century molecular biology. The paradigm tacitly stipulated that for each sequence there exists a single, well-organized protein structure. Yet, to sustain cell life, function requires (i) that there be more than a single structure, (ii) that there be switching between the structures, and (iii) that the structures be incompletely organized. These fundamental tenets called for an updated sequence-conformational ensemble-function paradigm. The powerful energy landscape idea, which is the foundation of modernized molecular biology, imported the conformational ensemble framework from physics and chemistry. This framework embraces the recognition that proteins are dynamic and are always interconverting between conformational states with varying energies. The more stable the conformation the more populated it is. The changes in the populations of the states are required for cell life. As an example, in vivo, under physiological conditions, wild type kinases commonly populate their more stable "closed", inactive, conformations. However, there are minor populations of the "open", ligand-free states. Upon their stabilization, e.g., by high affinity interactions or mutations, their ensembles shift to occupy the active states. Here we discuss the role of conformational propensities in function. We provide multiple examples of diverse systems, including protein kinases, lipid kinases, and Ras GTPases, discuss diverse conformational mechanisms, and provide a broad outlook on protein ensembles in the cell. We propose that the number of molecules in the active state (inactive for repressors), determine protein function, and that the dynamic, relative conformational propensities, rather than the rigid structures, are the hallmark of cell life.

Physicochemical effects of emerging exchanges on the spike protein's RBM of the SARS-CoV-2 Omicron subvariants BA.1-BA.5 and its influence on the biological properties and attributes developed by these subvariants

Emerging in South Africa, SARS-CoV-2 Omicron variant was marked by the expression of an exaggerated number of mutations throughout its genome and by the emergence of subvariants, whose attributes developed by them have been associated with amino acid exchanges that occur mainly in the RBM region of the spike protein. The RBM comprises a region within the RBD and is directly involved in the SARS-CoV-2 spike protein interaction with the host cell ACE2 receptor, during the infection mechanism and viral transmission. Defined as the region from aa 437 to aa 508, there are several residues in certain positions that interact directly with the human ACE-2 receptor during these processes. The occurrence of amino acid exchanges in these positions causes physicochemical alterations in the SARS-CoV-2 spike protein, which confer additional advantages and attributes to the agent. In addition, these exchanges serve as a basis for the characterization of new variants and subvariants of SARS-CoV-2. In this review, the amino acid exchanges that have occurred in the RBM of the subvariants BA.1 to BA.5 of SARS-CoV-2 that emerged from the Omicron are described. The physicochemical effects caused by them on spike protein are also described, as well as their influence on the biological properties and attributes developed by the subvariants BA.1, BA.2, BA.3, BA.4 and BA.5.

Interspecies transmission of SARS CoV-2 with special emphasis on viral mutations and ACE-2 receptor homology roles

COVID-19 outbreak was first reported in 2019, Wuhan, China. The spillover of the disease caused by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), to a wide range of pet, zoo, wild, and farm animals has emphasized potential zoonotic and reverse zoonotic viral transmission. Furthermore, it has evoked inquiries about susceptibility of different animal species to SARS-CoV-2 infection and role of these animals as viral reservoirs. Therefore, studying susceptible and non-susceptible hosts for SARS-CoV-2 infection could give a better understanding for the virus and will help in preventing further outbreaks. Here, we review structural aspects of SARS-CoV-2 spike protein, the effect of the different mutations observed in the spike protein, and the impact of ACE2 receptor variations in different animal hosts on inter-species transmission. Moreover, the SARS-CoV-2 spillover chain was reviewed. Combination of SARS-CoV-2 high mutation rate and homology of cellular ACE2 receptors enable the virus to transcend species barriers and facilitate its transmission between humans and animals.

The influence of single-point mutation D614G on the binding process between human angiotensin-converting enzyme 2 and the SARS-CoV-2 spike protein-an atomistic simulation study

SARS-CoV-2 has continuously evolved as changes in the genetic code occur during replication of the genome, with some of the mutations leading to higher transmission among human beings. The spike aspartic acid-614 to glycine (D614G) substitution in the spike represents a "more transmissible form of SARS-CoV-2" and occurs in all SARS-CoV-2 mutants. However, the underlying mechanism of the D614G substitution in virus infectivity has remained unclear. In this paper, we adopt molecular simulations to study the contact processes of the D614G mutant and wild-type (WT) spikes with hACE2. The interaction areas with hACE2 for the two spikes are completely different by visualizing the whole binding processes. The D614G mutant spike moves towards the hACE2 faster than the WT spike. We have also found that the receptor-binding domain (RBD) and N-terminal domain (NTD) of the D614G mutant extend more outwards than those of the WT spike. By analyzing the distances between the spikes and hACE2, the changes of number of hydrogen bonds and interaction energy, we suggest that the increased infectivity of the D614G mutant is not possibly related to the binding strength, but to the binding velocity and conformational change of the mutant spike. This work reveals the impact of D614G substitution on the infectivity of the SARS-CoV-2, and hopefully could provide a rational explanation of interaction mechanisms for all the SARS-CoV-2 mutants.

Internal Normal Mode Analysis Applied to RNA Flexibility and Conformational Changes

We investigated the capability of internal normal modes to reproduce RNA flexibility and predict observed RNA conformational changes and, notably, those induced by the formation of RNA-protein and RNA-ligand complexes. Here, we extended our iNMA approach developed for proteins to study RNA molecules using a simplified representation of the RNA structure and its potential energy. Three data sets were also created to investigate different aspects. Despite all the approximations, our study shows that iNMA is a suitable method to take into account RNA flexibility and describe its conformational changes opening the route to its applicability in any integrative approach where these properties are crucial.

The diverse role of heparan sulfate and other GAGs in SARS-CoV-2 infections and therapeutics

In December 2019, the global coronavirus disease 2019 (COVID-19) pandemic began in Wuhan, China. COVID-19 is caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which infects host cells primarily through the angiotensin-converting enzyme 2 (ACE2) receptor. In addition to ACE2, several studies have shown the importance of heparan sulfate (HS) on the host cell surface as a co-receptor for SARS-CoV-2-binding. This insight has driven research into antiviral therapies, aimed at inhibiting the HS co-receptor-binding, e.g., by glycosaminoglycans (GAGs), a family of sulfated polysaccharides that includes HS. Several GAGs, such as heparin (a highly sulfated analog of HS), are used to treat various health indications, including COVID-19. This review is focused on current research on the involvement of HS in SARS-CoV-2 infection, implications of viral mutations, as well as the use of GAGs and other sulfated polysaccharides as antiviral agents.

Evaluating the effect of SARS-CoV-2 spike mutations with a linear doubly robust learner

Driven by various mutations on the viral Spike protein, diverse variants of SARS-CoV-2 have emerged and prevailed repeatedly, significantly prolonging the pandemic. This phenomenon necessitates the identification of key Spike mutations for fitness enhancement. To address the need, this manuscript formulates a well-defined framework of causal inference methods for evaluating and identifying key Spike mutations to the viral fitness of SARS-CoV-2. In the context of large-scale genomes of SARS-CoV-2, it estimates the statistical contribution of mutations to viral fitness across lineages and therefore identifies important mutations. Further, identified key mutations are validated by computational methods to possess functional effects, including Spike stability, receptor-binding affinity, and potential for immune escape. Based on the effect score of each mutation, individual key fitness-enhancing mutations such as D614G and T478K are identified and studied. From individual mutations to protein domains, this paper recognizes key protein regions on the Spike protein, including the receptor-binding domain and the N-terminal domain. This research even makes further efforts to investigate viral fitness via mutational effect scores, allowing us to compute the fitness score of different SARS-CoV-2 strains and predict their transmission capacity based solely on their viral sequence. This prediction of viral fitness has been validated using BA.2.12.1, which is not used for regression training but well fits the prediction. To the best of our knowledge, this is the first research to apply causal inference models to mutational analysis on large-scale genomes of SARS-CoV-2. Our findings produce innovative and systematic insights into SARS-CoV-2 and promotes functional studies of its key mutations, serving as reliable guidance about mutations of interest.

Force-tuned avidity of spike variant-ACE2 interactions viewed on the single-molecule level

Recent waves of COVID-19 correlate with the emergence of the Delta and the Omicron variant. We report that the Spike trimer acts as a highly dynamic molecular caliper, thereby forming up to three tight bonds through its RBDs with ACE2 expressed on the cell surface. The Spike of both Delta and Omicron (B.1.1.529) Variant enhance and markedly prolong viral attachment to the host cell receptor ACE2, as opposed to the early Wuhan-1 isolate. Delta Spike shows rapid binding of all three Spike RBDs to three different ACE2 molecules with considerably increased bond lifetime when compared to the reference strain, thereby significantly amplifying avidity. Intriguingly, Omicron (B.1.1.529) Spike displays less multivalent bindings to ACE2 molecules, yet with a ten time longer bond lifetime than Delta. Delta and Omicron (B.1.1.529) Spike variants enhance and prolong viral attachment to the host, which likely not only increases the rate of viral uptake, but also enhances the resistance of the variants against host-cell detachment by shear forces such as airflow, mucus or blood flow. We uncover distinct binding mechanisms and strategies at single-molecule resolution, employed by circulating SARS-CoV-2 variants to enhance infectivity and viral transmission.

Sequence-sensitive elastic network captures dynamical features necessary for miR-125a maturation

The Elastic Network Contact Model (ENCoM) is a coarse-grained normal mode analysis (NMA) model unique in its all-atom sensitivity to the sequence of the studied macromolecule and thus to the effect of mutations. We adapted ENCoM to simulate the dynamics of ribonucleic acid (RNA) molecules, benchmarked its performance against other popular NMA models and used it to study the 3D structural dynamics of human microRNA miR-125a, leveraging high-throughput experimental maturation efficiency data of over 26 000 sequence variants. We also introduce a novel way of using dynamical information from NMA to train multivariate linear regression models, with the purpose of highlighting the most salient contributions of dynamics to function. ENCoM has a similar performance profile on RNA than on proteins when compared to the Anisotropic Network Model (ANM), the most widely used coarse-grained NMA model; it has the advantage on predicting large-scale motions while ANM performs better on B-factors prediction. A stringent benchmark from the miR-125a maturation dataset, in which the training set contains no sequence information in common with the testing set, reveals that ENCoM is the only tested model able to capture signal beyond the sequence. This ability translates to better predictive power on a second benchmark in which sequence features are shared between the train and test sets. When training the linear regression model using all available data, the dynamical features identified as necessary for miR-125a maturation point to known patterns but also offer new insights into the biogenesis of microRNAs. Our novel approach combining NMA with multivariate linear regression is generalizable to any macromolecule for which relatively high-throughput mutational data is available.

Novel coronavirus mutations: Vaccine development and challenges

The ongoing global pandemic of novel coronavirus pneumonia (COVID-19) caused by the SARS-CoV-2 has a significant impact on global health and economy system. In this context, there have been some landmark advances in vaccine development. Over 100 new coronavirus vaccine candidates have been approved for clinical trials, with ten WHO-approved vaccines including four inactivated virus vaccines, two mRNA vaccines, three recombinant viral vectored vaccines and one protein subunit vaccine on the "Emergency Use Listing". Although the SARS-CoV-2 has an internal proofreading mechanism, there have been a number of mutations emerged in the pandemic affecting its transmissibility, pathogenicity and immunogenicity. Of these, mutations in the spike (S) protein and the resultant mutant variants have posed new challenges for vaccine development and application. In this review article, we present an overview of vaccine development, the prevalence of new coronavirus variants and their impact on protective efficacy of existing vaccines and possible immunization strategies coping with the viral mutation and diversity.

Probing the mutational landscape of the SARS-CoV-2 spike protein via quantum mechanical modeling of crystallographic structures

We employ a recently developed complexity-reduction quantum mechanical (QM-CR) approach, based on complexity reduction of density functional theory calculations, to characterize the interactions of the SARS-CoV-2 spike receptor binding domain (RBD) with ACE2 host receptors and antibodies. QM-CR operates via ab initio identification of individual amino acid residue's contributions to chemical binding and leads to the identification of the impact of point mutations. Here, we especially focus on the E484K mutation of the viral spike protein. We find that spike residue 484 hinders the spike's binding to the human ACE2 receptor (hACE2). In contrast, the same residue is beneficial in binding to the bat receptor Rhinolophus macrotis ACE2 (macACE2). In agreement with empirical evidence, QM-CR shows that the E484K mutation allows the spike to evade categories of neutralizing antibodies like C121 and C144. The simulation also shows how the Delta variant spike binds more strongly to hACE2 compared to the original Wuhan strain, and predicts that a E484K mutation can further improve its binding. Broad agreement between the QM-CR predictions and experimental evidence supports the notion that ab initio modeling has now reached the maturity required to handle large intermolecular interactions central to biological processes.

The Comparison of Mutational Progression in SARS-CoV-2: A Short Updated Overview

The COVID-19 pandemic has impacted the world population adversely, posing a threat to human health. In the past few years, various strains of SARS-CoV-2, each with different mutations in its structure, have impacted human health in negative ways. The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) mutations influence the virulence, antibody evasion, and Angiotensin-converting enzyme 2 (ACE2) affinity of the virus. These mutations are essential to understanding how a new strain of SARS-CoV-2 has changed and its possible effects on the human body. This review provides an insight into the spike mutations of SARS-CoV-2 variants. As the current scientific data offer a scattered outlook on the various type of mutations, we aimed to categorize the mutations of Beta (B.1.351), Gamma (P.1), Delta (B.1.612.2), and Omicron (B.1.1.529) systematically according to their location in the subunit 1 (S1) and subunit 2 (S2) domains and summarized their consequences as a result. We also compared the miscellany of mutations that have emerged in all four variants to date. The comparison shows that mutations such as D614G and N501Y have emerged in all four variants of concern and that all four variants have multiple mutations within the N-terminal domain (NTD), as in the case of the Delta variant. Other mutations are scattered in the receptor binding domain (RBD) and subdomain 2 (SD2) of the S1 domain. Mutations in RBD or NTD are often associated with antibody evasion. Few mutations lie in the S2 domain in the Beta, Gamma, and Delta variants. However, in the Omicron variant many mutations occupy the S2 domain, hinting towards a much more evasive virus.

Human coronaviruses: The emergence of SARS-CoV-2 and management of COVID-19

To date, a total of seven human coronaviruses (HCoVs) have been identified, all of which are important respiratory pathogens. Recently, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of coronavirus disease 2019 (COVID-19), has led to a global pandemic causing millions of infections and deaths. Here, we summarize the discovery and fundamental virology of HCoVs, discuss their zoonotic transmission and highlight the weak species barrier of SARS-CoV-2. We also discuss the possible origins of SARS-CoV-2 variants of concern identified to date and discuss the experimental challenges in characterizing mutations of interest and propose methods to circumvent them. As the COVID-19 treatment and prevention landscape rapidly evolves, we summarize current therapeutics and vaccines, and their implications on SARS-CoV-2 variants. Finally, we explore how interspecies transmission of SARS-CoV-2 may drive the emergence of novel strains, how disease severity may evolve and how COVID-19 will likely continue to burden healthcare systems globally.

Applications of Molecular Dynamics Simulation in Protein Study

Molecular Dynamics (MD) Simulations is increasingly used as a powerful tool to study protein structure-related questions. Starting from the early simulation study on the photoisomerization in rhodopsin in 1976, MD Simulations has been used to study protein function, protein stability, protein-protein interaction, enzymatic reactions and drug-protein interactions, and membrane proteins. In this review, we provide a brief review for the history of MD Simulations application and the current status of MD Simulations applications in protein studies.

Raman Fingerprints of the SARS-CoV-2 Delta Variant and Mechanisms of Its Instantaneous Inactivation by Silicon Nitride Bioceramics

Raman spectroscopy uncovered molecular scale markers of the viral structure of the SARS-CoV-2 Delta variant and related viral inactivation mechanisms at the biological interface with silicon nitride (Si3N4) bioceramics. A comparison of Raman spectra collected on the TY11-927 variant (lineage B.1.617.2; simply referred to as the Delta variant henceforth) with those of the JPN/TY/WK-521 variant (lineage B.1.617.1; referred to as the Kappa variant or simply as the Japanese isolate henceforth) revealed the occurrence of key mutations of the spike receptor together with profound structural differences in the molecular structure/symmetry of sulfur-containing amino acid and altered hydrophobic interactions of the tyrosine residue. Additionally, different vibrational fractions of RNA purines and pyrimidines and dissimilar protein secondary structures were also recorded. Despite mutations, hydrolytic reactions at the surface of silicon nitride (Si3N4) bioceramics induced instantaneous inactivation of the Delta variant at the same rate as that of the Kappa variant. Contact between virions and micrometric Si3N4 particles yielded post-translational deimination of arginine spike residues, methionine sulfoxidation, tyrosine nitration, and oxidation of RNA purines to form formamidopyrimidines. Si3N4 bioceramics proved to be a safe and effective inorganic compound for instantaneous environmental sanitation.

A potent virucidal activity of functionalized TiO2 nanoparticles adsorbed with flavonoids against SARS-CoV-2

Abstract

The coronavirus SARS-CoV-2 has caused a pandemic with > 550 millions of cases and > 6 millions of deaths worldwide. Medical management of COVID-19 relies on supportive care as no specific targeted therapies are available yet. Given its devastating effects on the economy and mental health, it is imperative to develop novel antivirals. An ideal candidate will be an agent that blocks the early events of viral attachment and cell entry, thereby preventing viral infection and spread. This work reports functionalized titanium dioxide (TiO₂)-based nanoparticles adsorbed with flavonoids that block SARS-CoV-2 entry and fusion. Using molecular docking analysis, two flavonoids were chosen for their specific binding to critical regions of the SARS-CoV-2 spike glycoprotein that interacts with the host cell angiotensin-converting enzyme-2 (ACE-2) receptor. These flavonoids were adsorbed onto TiO₂ functionalized nanoparticles (FTNP). This new nanoparticulate compound was assayed in vitro against two different coronaviruses; HCoV 229E and SARS-CoV-2, in both cases a clear antiviral effect was observed. Furthermore, using a reporter-based cell culture model, a potent antiviral activity is demonstrated. The adsorption of flavonoids to functionalized TiO₂ nanoparticles induces a ~ threefold increase of that activity. These studies also indicate that FTNP interferes with the SARS-CoV-2 spike, impairing the cell fusion mechanism.

Key points/Highlights

• Unique TiO₂nanoparticles displaying flavonoid showed potent anti-SARS-CoV-2 activity.

• The nanoparticles precisely targeting SARS-CoV-2 were quantitatively verified by cell infectivity in vitro.

• Flavonoids on nanoparticles impair the interactions between the spike glycoprotein and ACE-2 receptor.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s00253-022-12112-9.

Independent acquisition of short insertions at the RIR1 site in the spike N-terminal domain of the SARS-CoV-2 BA.2 lineage

Although the major SARS‐CoV‐2 omicron lineages share over 30 non‐synonymous substitutions in the spike glycoprotein, they show several unique mutations that were acquired after their ancestral split. One of the most intriguing mutations associated with BA.1 is the presence of the inserted tripeptide Glu‐Pro‐Glu within the N‐terminal domain, at a site that had previously independently acquired short insertions in several other SARS‐CoV‐2 lineages. Although the functional implications of the small nucleotide sequences found at this insertion hotspot, named RIR1, are still unclear, we have previously hypothesized that they may play a compensatory role in counterbalancing minor fitness deficits associated with other co‐occurring spike non‐synonymous mutations.

Here we show that similar insertion events have independently occurred at RIR1 at least 20 times in early 2022 within the BA.2 lineage, being occasionally associated with significant community transmission. One of these omicron sublineages, characterized by a Ser‐Gly‐Arg insertion in position 212, has been responsible of over 4,000 documented covid‐19 cases worldwide between January and July 2022, for the most part concentrated in Denmark, where it reached a national prevalence close to 4% (10% in the Nordjylland region) in mid‐May. Although the concurrent spread of the BA.2.12.1, BA.4 and BA.5 lineages led to the rapid decline of this BA.2 sublineage, the independent acquisition of several other RIR1 insertions on a BA.2 genomic background suggests that these events may provide a slight fitness advantage. Therefore, we they should be carefully monitored in the upcoming months in other emerging omicron‐related lineages, including BA.5.

This article is protected by copyright. All rights reserved

Decreased Interfacial Dynamics Caused by the N501Y Mutation in the SARS-CoV-2 S1 Spike:ACE2 Complex

Coronavirus Disease of 2019 (COVID-19) caused by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has resulted in a massive health crisis across the globe, with some genetic variants gaining enhanced infectivity and competitive fitness, and thus significantly aggravating the global health concern. In this regard, the recent SARS-CoV-2 alpha, beta, and gamma variants (B.1.1.7, B.1.351, and P.1 lineages, respectively) are of great significance in that they contain several mutations that increase their transmission rates as evident from clinical reports. By the end of March 2021, these variants were accounting for about two-thirds of SARS-CoV-2 variants circulating worldwide. Specifically, the N501Y mutation in the S1 spike receptor binding domain (S1-RBD) of these variants have been reported to increase its affinity for ACE2, although the basis for this is not entirely clear yet. Here, we dissect the mechanism underlying the increased binding affinity of the N501Y mutant for ACE2 using molecular dynamics (MD) simulations of the available ACE2-S1-RBD complex structure (6M0J) and show a prolonged and stable interfacial interaction of the N501Y mutant S1-RBD with ACE2 compared to the wild type S1-RBD. Additionally, we find that the N501Y mutant S1-RBD displays altered dynamics that likely aids in its enhanced interaction with ACE2. By elucidating a mechanistic basis for the increased affinity of the N501Y mutant S1-RBD for ACE2, we believe that the results presented here will aid in developing therapeutic strategies against SARS-CoV-2 including designing of therapeutic agents targeting the ACE2-S1-RBD interaction.

An Electrostatically-steered Conformational Selection Mechanism Promotes SARS-CoV-2 Spike Protein Variation

After two years since the outbreak, the COVID-19 pandemic remains a global public health emergency. SARS-CoV-2 variants with substitutions on the spike (S) protein emerge increasing the risk of immune evasion and cross-species transmission. Here, we analyzed the evolution of the S protein as recorded in 276,712 samples collected before the start of vaccination efforts. Our analysis shows that most variants destabilize the S protein trimer, increase its conformational heterogeneity and improve the odds of the recognition by the host cell receptor. Most frequent substitutions promote overall hydrophobicity by replacing charged amino acids, reducing stabilizing local interactions in the unbound S protein trimer. Moreover, our results identify "forbidden" regions that rarely show any sequence variation, and which are related to conformational changes occurring upon fusion. These results are significant for understanding the structure and function of SARS-CoV-2 related proteins which is a critical step in vaccine development and epidemiological surveillance.

Methodology-Centered Review of Molecular Modeling, Simulation, and Prediction of SARS-CoV-2

Despite tremendous efforts in the past two years, our understanding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), virus-host interactions, immune response, virulence, transmission, and evolution is still very limited. This limitation calls for further in-depth investigation. Computational studies have become an indispensable component in combating coronavirus disease 2019 (COVID-19) due to their low cost, their efficiency, and the fact that they are free from safety and ethical constraints. Additionally, the mechanism that governs the global evolution and transmission of SARS-CoV-2 cannot be revealed from individual experiments and was discovered by integrating genotyping of massive viral sequences, biophysical modeling of protein-protein interactions, deep mutational data, deep learning, and advanced mathematics. There exists a tsunami of literature on the molecular modeling, simulations, and predictions of SARS-CoV-2 and related developments of drugs, vaccines, antibodies, and diagnostics. To provide readers with a quick update about this literature, we present a comprehensive and systematic methodology-centered review. Aspects such as molecular biophysics, bioinformatics, cheminformatics, machine learning, and mathematics are discussed. This review will be beneficial to researchers who are looking for ways to contribute to SARS-CoV-2 studies and those who are interested in the status of the field.

Physicochemical effect of the N501Y, E484K/Q, K417N/T, L452R and T478K mutations on the SARS-CoV-2 spike protein RBD and its influence on agent fitness and on attributes developed by emerging variants of concern

The spike protein comprises one of the main structural components of SARS-CoV-2 because it is directly involved in the infection process and viral transmission, and also because of its immunogenic properties, as an inducer of the protective antibodies production and as a vaccine component. The occurrence of mutations in this region or in other the virus genome regions, comprises a natural phenomenon in its evolution. However, they also occur due to the selective immune pressure, to which the agent is continuously subjected, especially in the spike protein immunodominant regions, such as the RBD. Mutations in the spike protein can change the virus' fitness, increasing its affinity for target cells, its transmissibility and its virulence. In addition, these mutations can giving it the potential ability to evade the protective antibodies action obtained from convalescent sera or vaccine origin, as well as those used in therapy, which may favor the virus expansion and compromise the infection control. Five mutations N501Y, E484K/Q, K417N/T, L452R and T478K, located in the spike protein RBD, have had a greater impact because they are associated with new attributes developed by the virus, which characterize the emerging variants of concern (VOCs) of SARS-Cov-2 identified so far. The occurrence of these mutations induces complex physicochemical effects that can alter the spike protein's structure and its function, which in turn, lead to changes in the agents' fitness. This manuscript discusses the attributes of VOCs associated with the physicochemical effects caused by the aforementioned mutations.

Molecular characteristics, immune evasion, and impact of SARS-CoV-2 variants

The persistent COVID-19 pandemic since 2020 has brought an enormous public health burden to the global society and is accompanied by various evolution of the virus genome. The consistently emerging SARS-CoV-2 variants harboring critical mutations impact the molecular characteristics of viral proteins and display heterogeneous behaviors in immune evasion, transmissibility, and the clinical manifestation during infection, which differ each strain and endow them with distinguished features during populational spread. Several SARS-CoV-2 variants, identified as Variants of Concern (VOC) by the World Health Organization, challenged global efforts on COVID-19 control due to the rapid worldwide spread and enhanced immune evasion from current antibodies and vaccines. Moreover, the recent Omicron variant even exacerbated the global anxiety in the continuous pandemic. Its significant evasion from current medical treatment and disease control even highlights the necessity of combinatory investigation of the mutational pattern and influence of the mutations on viral dynamics against populational immunity, which would greatly facilitate drug and vaccine development and benefit the global public health policymaking. Hence in this review, we summarized the molecular characteristics, immune evasion, and impacts of the SARS-CoV-2 variants and focused on the parallel comparison of different variants in mutational profile, transmissibility and tropism alteration, treatment effectiveness, and clinical manifestations, in order to provide a comprehensive landscape for SARS-CoV-2 variant research.

The past, current and future epidemiological dynamic of SARS-CoV-2

SARS-CoV-2, the agent of the COVID-19 pandemic, emerged in late 2019 in China, and rapidly spread throughout the world to reach all continents. As the virus expanded in its novel human host, viral lineages diversified through the accumulation of around two mutations a month on average. Different viral lineages have replaced each other since the start of the pandemic, with the most successful Alpha, Delta and Omicron variants of concern (VoCs) sequentially sweeping through the world to reach high global prevalence. Neither Alpha nor Delta was characterized by strong immune escape, with their success coming mainly from their higher transmissibility. Omicron is far more prone to immune evasion and spread primarily due to its increased ability to (re-)infect hosts with prior immunity. As host immunity reaches high levels globally through vaccination and prior infection, the epidemic is expected to transition from a pandemic regime to an endemic one where seasonality and waning host immunization are anticipated to become the primary forces shaping future SARS-CoV-2 lineage dynamics. In this review, we consider a body of evidence on the origins, host tropism, epidemiology, genomic and immunogenetic evolution of SARS-CoV-2 including an assessment of other coronaviruses infecting humans. Considering what is known so far, we conclude by delineating scenarios for the future dynamic of SARS-CoV-2, ranging from the good-circulation of a fifth endemic 'common cold' coronavirus of potentially low virulence, the bad-a situation roughly comparable with seasonal flu, and the ugly-extensive diversification into serotypes with long-term high-level endemicity.

SPEAR: Systematic ProtEin AnnotatoR

SUMMARY

We present SPEAR, a lightweight and rapid SARS-CoV-2 variant annotation and scoring tool, for identifying mutations contributing to potential immune escape and transmissibility (ACE2 binding) at point of sequencing. SPEAR can be used in the field to evaluate genomic surveillance results in real-time and features a powerful interactive data visualisation report.

AVAILABILITY AND IMPLEMENTATION

SPEAR and documentation are freely available on GitHub: https://github.com/m-crown/SPEAR and is implemented in Python and installable via Conda environment.

SUPPLEMENTAL

Supplementary data are available at Bioinformatics online.

Imaging Severity COVID-19 Assessment in Vaccinated and Unvaccinated Patients: Comparison of the Different Variants in a High Volume Italian Reference Center

Purpose: To analyze the vaccine effect by comparing five groups: unvaccinated patients with Alpha variant, unvaccinated patients with Delta variant, vaccinated patients with Delta variant, unvaccinated patients with Omicron variant, and vaccinated patients with Omicron variant, assessing the “gravity” of COVID-19 pulmonary involvement, based on CT findings in critically ill patients admitted to Intensive Care Unit (ICU). Methods: Patients were selected by ICU database considering the period from December 2021 to 23 March 2022, according to the following inclusion criteria: patients with proven Omicron variant COVID-19 infection with known COVID-19 vaccination with at least two doses and with chest Computed Tomography (CT) study during ICU hospitalization. Wee also evaluated the ICU database considering the period from March 2020 to December 2021, to select unvaccinated consecutive patients with Alpha variant, subjected to CT study, consecutive unvaccinated and vaccinated patients with Delta variant, subjected to CT study, and, consecutive unvaccinated patients with Omicron variant, subjected to CT study. CT images were evaluated qualitatively using a severity score scale of 5 levels (none involvement, mild: ≤25% of involvement, moderate: 26−50% of involvement, severe: 51−75% of involvement, and critical involvement: 76−100%) and quantitatively, using the Philips IntelliSpace Portal clinical application CT COPD computer tool. For each patient the lung volumetry was performed identifying the percentage value of aerated residual lung volume. Non-parametric tests for continuous and categorical variables were performed to assess statistically significant differences among groups. Results: The patient study group was composed of 13 vaccinated patients affected by the Omicron variant (Omicron V). As control groups we identified: 20 unvaccinated patients with Alpha variant (Alpha NV); 20 unvaccinated patients with Delta variant (Delta NV); 18 vaccinated patients with Delta variant (Delta V); and 20 unvaccinated patients affected by the Omicron variant (Omicron NV). No differences between the groups under examination were found (p value > 0.05 at Chi square test) in terms of risk factors (age, cardiovascular diseases, diabetes, immunosuppression, chronic kidney, cardiac, pulmonary, neurologic, and liver disease, etc.). A different median value of aerated residual lung volume was observed in the Delta variant groups: median value of aerated residual lung volume was 46.70% in unvaccinated patients compared to 67.10% in vaccinated patients. In addition, in patients with Delta variant every other extracted volume by automatic tool showed a statistically significant difference between vaccinated and unvaccinated group. Statistically significant differences were observed for each extracted volume by automatic tool between unvaccinated patients affected by Alpha variant and vaccinated patients affected by Delta variant of COVID-19. Good statistically significant correlations among volumes extracted by automatic tool for each lung lobe and overall radiological severity score were obtained (ICC range 0.71−0.86). GGO was the main sign of COVID-19 lesions on CT images found in 87 of the 91 (95.6%) patients. No statistically significant differences were observed in CT findings (ground glass opacities (GGO), consolidation or crazy paving sign) among patient groups. Conclusion: In our study, we showed that in critically ill patients no difference were observed in terms of severity of disease or exitus, between unvaccinated and vaccinated patients. The only statistically significant differences were observed, with regard to the severity of COVID-19 pulmonary parenchymal involvement, between unvaccinated patients affected by Alpha variant and vaccinated patients affected by Delta variant, and between unvaccinated patients with Delta variant and vaccinated patients with Delta variant.

Synthesis, crystal structure elucidation, DFT analysis, drug-likeness and ADMET evaluation and molecular docking studies of triazole derivatives: Binary inhibition of spike protein and ACE2 receptor protein of COVID-19

The recent incidence of terrible acute respiratory syndrome coronavirus 2 (SARS CoV-2) has presently experienced some noteworthy mutations since its discovery in 2019 in Wuhan, China. The present research work focuses on the synthesis of three triazole derivatives (BMTPP, BMTTP, and BMTIP) and their inhibition activities against SARS-Cov-2 spike and ACE2 receptor proteins. The crystal structure for BMTTP was determined by the SCXRD method and optimized geometrical parameters for the three triazole derivatives were obtained by DFT calculations. HOMO-LUMO, Global reactive descriptors [GRD], and Molecular electrostatic potential (MEP) investigations exposed that all three compounds have biological properties. The drug-likeness ability of the synthesized compounds was examined using Molinspiration and a pre-ADMET online Server. Further, to explore the binding nature of three synthesized compounds with SARS-Cov-2 spike proteins/ACE2 receptor molecular docking studies were executed. The outcomes we obtained from molecular docking simulation studies suggest that the synthesized triazole derivatives may be well utilized as curing medicines against COVID-19. Ultimately, animal tests and precise clinical tests are required to prove the potent nature of these compounds against COVID-19. Finally, the present outcomes must be proved to utilize in-vitro and in-vivo antiviral methods.

Probing structural basis for enhanced binding of SARS-CoV-2 P.1 variant spike protein with the human ACE2 receptor

The initial step of infection by severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) involves the binding of receptor binding domain (RBD) of the spike protein to the angiotensin converting enzyme 2 (ACE2) receptor. Each successive wave of SARS‐CoV‐2 reports emergence of many new variants, which is associated with mutations in the RBD as well as other parts of the spike protein. These mutations are reported to have enhanced affinity towards the ACE2 receptor as well as are also crucial for the virus transmission. Many computational and experimental studies have demonstrated the effect of individual mutation on the RBD‐ACE2 binding. However, the cumulative effect of mutations on the RBD and away from the RBD was not investigated in detail. We report here a comparative analysis on the structural communication and dynamics of the RBD and truncated S1 domain of spike protein in complex with the ACE2 receptor from SARS‐CoV‐2 wild type and its P.1 variant. Our integrative network and dynamics approaches highlighted a subtle conformational changes in the RBD as well as truncated S1 domain of spike protein at the protein contact level, responsible for the increased affinity with the ACE2 receptor. Moreover, our study also identified the commonalities and differences in the dynamics of the interactions between spike protein of SARS‐CoV‐2 wild type and its P.1 variant with the ACE2 receptor. Further, our investigation yielded an understanding towards identification of the unique RBD residues crucial for the interaction with the ACE2 host receptor. Overall, the study provides an insight for designing better therapeutics against the circulating P.1 variants as well as other future variants.

Resilience of Spike-Specific Immunity Induced by COVID-19 Vaccines against SARS-CoV-2 Variants

The outbreak of SARS-CoV-2 leading to the declaration of the COVID-19 global pandemic has led to the urgent development and deployment of several COVID-19 vaccines. Many of these new vaccines, including those based on mRNA and adenoviruses, are aimed to generate neutralizing antibodies against the spike glycoprotein, which is known to bind to the receptor angiotensin converting enzyme 2 (ACE2) in host cells via the receptor-binding domain (RBD). Antibodies binding to this domain can block the interaction with the receptor and prevent viral entry into the cells. Additionally, these vaccines can also induce spike-specific T cells which could contribute to providing protection against the virus. However, the emergence of new SARS-CoV-2 variants can impair the immunity generated by COVID-19 vaccines if mutations occur in cognate epitopes, precluding immune recognition. Here, we evaluated the chance of five SARS-CoV-2 variants of concern (VOCs), Alpha, Beta, Gamma, Delta and Omicron, to escape spike-specific immunity induced by vaccines. To that end, we examined the impact of the SARS-CoV-2 variant mutations on residues located on experimentally verified spike-specific epitopes, deposited at the Immune Epitope Database, that are targeted by neutralizing antibodies or recognized by T cells. We found about 300 of such B cell epitopes, which were largely overlapping, and could be grouped into 54 B cell epitope clusters sharing ≥ 7 residues. Most of the B cell epitope clusters map in the RBD domain (39 out of 54) and 20%, 50%, 37%, 44% and 57% of the total are mutated in SARS-CoV-2 Alpha, Beta, Gamma, Delta and Omicron variants, respectively. We also found 234 experimentally verified CD8 and CD4 T cell epitopes that were distributed evenly throughout the spike protein. Interestingly, in each SARS-CoV-2 VOC, over 87% and 79% of CD8 and CD4 T cell epitopes, respectively, are not mutated. These observations suggest that SARS-CoV-2 VOCs—particularly the Omicron variant—may be prone to escape spike-specific antibody immunity, but not cellular immunity, elicited by COVID-19 vaccines.

Severe Acute Respiratory Syndrome Coronavirus 2 Variants of Concern: A Perspective for Emerging More Transmissible and Vaccine-Resistant Strains

Novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOC) are constantly threatening global public health. With no end date, the pandemic persists with the emergence of novel variants that threaten the effectiveness of diagnostic tests and vaccines. Mutations in the Spike surface protein of the virus are regularly observed in the new variants, potentializing the emergence of novel viruses with different tropism from the current ones, which may change the severity and symptoms of the disease. Growing evidence has shown that mutations are being selected in favor of variants that are more capable of evading the action of neutralizing antibodies. In this context, the most important factor guiding the evolution of SARS-CoV-2 is its interaction with the host's immune system. Thus, as current vaccines cannot block the transmission of the virus, measures complementary to vaccination, such as the use of masks, hand hygiene, and keeping environments ventilated remain essential to delay the emergence of new variants. Importantly, in addition to the involvement of the immune system in the evolution of the virus, we highlight several chemical parameters that influence the molecular interactions between viruses and host cells during invasion and are also critical tools making novel variants more transmissible. In this review, we dissect the impacts of the Spike mutations on biological parameters such as (1) the increase in Spike binding affinity to hACE2; (2) bound time for the receptor to be cleaved by the proteases; (3) how mutations associate with the increase in RBD up-conformation state in the Spike ectodomain; (4) expansion of uncleaved Spike protein in the virion particles; (5) increment in Spike concentration per virion particles; and (6) evasion of the immune system. These factors play key roles in the fast spreading of SARS-CoV-2 variants of concern, including the Omicron.

Computer Simulations and Network-Based Profiling of Binding and Allosteric Interactions of SARS-CoV-2 Spike Variant Complexes and the Host Receptor: Dissecting the Mechanistic Effects of the Delta and Omicron Mutations

In this study, we combine all-atom MD simulations and comprehensive mutational scanning of S-RBD complexes with the angiotensin-converting enzyme 2 (ACE2) host receptor in the native form as well as the S-RBD Delta and Omicron variants to (a) examine the differences in the dynamic signatures of the S-RBD complexes and (b) identify the critical binding hotspots and sensitivity of the mutational positions. We also examined the differences in allosteric interactions and communications in the S-RBD complexes for the Delta and Omicron variants. Through the perturbation-based scanning of the allosteric propensities of the SARS-CoV-2 S-RBD residues and dynamics-based network centrality and community analyses, we characterize the global mediating centers in the complexes and the nature of local stabilizing communities. We show that a constellation of mutational sites (G496S, Q498R, N501Y and Y505H) correspond to key binding energy hotspots and also contribute decisively to the key interfacial communities that mediate allosteric communications between S-RBD and ACE2. These Omicron mutations are responsible for both favorable local binding interactions and long-range allosteric interactions, providing key functional centers that mediate the high transmissibility of the virus. At the same time, our results show that other mutational sites could provide a "flexible shield" surrounding the stable community network, thereby allowing the Omicron virus to modulate immune evasion at different epitopes, while protecting the integrity of binding and allosteric interactions in the RBD-ACE2 complexes. This study suggests that the SARS-CoV-2 S protein may exploit the plasticity of the RBD to generate escape mutants, while engaging a small group of functional hotspots to mediate efficient local binding interactions and long-range allosteric communications with ACE2.

Landscape-Based Protein Stability Analysis and Network Modeling of Multiple Conformational States of the SARS-CoV-2 Spike D614G Mutant: Conformational Plasticity and Frustration-Induced Allostery as Energetic Drivers of Highly Transmissible Spike Variants

The structural and functional studies of the SARS-CoV-2 spike protein variants revealed an important role of the D614G mutation that is shared across many variants of concern (VOCs), suggesting the effect of this mutation on the enhanced virus infectivity and transmissibility. The recent structural and biophysical studies provided important evidence about multiple conformational substates of the D614G spike protein. The development of a plausible mechanistic model that can explain the experimental observations from a more unified thermodynamic perspective is an important objective of the current work. In this study, we employed efficient and accurate coarse-grained simulations of multiple structural substates of the D614G spike trimers together with the ensemble-based mutational frustration analysis to characterize the dynamics signatures of the conformational landscapes. By combining the local frustration profiling of the conformational states with residue-based mutational scanning of protein stability and network analysis of allosteric interactions and communications, we determine the patterns of mutational sensitivity in the functional regions and sites of variants. We found that the D614G mutation may induce a considerable conformational adaptability of the open states in the SARS-CoV-2 spike protein without compromising the folding stability and integrity of the spike protein. The results suggest that the D614G mutant may employ a hinge-shift mechanism in which the dynamic couplings between the site of mutation and the interprotomer hinge modulate the interdomain interactions, global mobility change, and the increased stability of the open form. This study proposes that mutation-induced modulation of the conformational flexibility and energetic frustration at the interprotomer interfaces may serve as an efficient mechanism for allosteric regulation of the SARS-CoV-2 spike proteins.

Clinical Characteristics, Transmissibility, Pathogenicity, Susceptible Populations, and Re-infectivity of Prominent COVID-19 Variants

In addition to the rapid, global spread of SARS-CoV-2, new and comparatively more contagious variants are of considerable concern. These emerging mutations have become a threat to the global public health, creating COVID-19 surges in different countries. However, information on these emerging variants is limited and scattered. In this review, we discuss new variants that have emerged worldwide and identify several variants of concern, such as B.1.1.7, B.1.351, P.1, B.1.617.2 and B.1.1.529, and their basic characteristics. Other significant variants such as C.37, B.1.621, B.1.525, B.1.526, AZ.5, C.1.2, and B.1.617.1 are also discussed. This review highlights the clinical characteristics of these variants, including transmissibility, pathogenicity, susceptible population, and re-infectivity. It provides the latest information on the recent variants of SARS-CoV-2. The summary of this information will help researchers formulate reasonable strategies to curb the ongoing COVID-19 pandemic.

Structural and Computational Studies of the SARS-CoV-2 Spike Protein Binding Mechanisms with Nanobodies: From Structure and Dynamics to Avidity-Driven Nanobody Engineering

Nanobodies provide important advantages over traditional antibodies, including their smaller size and robust biochemical properties such as high thermal stability, high solubility, and the ability to be bioengineered into novel multivalent, multi-specific, and high-affinity molecules, making them a class of emerging powerful therapies against SARS-CoV-2. Recent research efforts on the design, protein engineering, and structure-functional characterization of nanobodies and their binding with SARS-CoV-2 S proteins reflected a growing realization that nanobody combinations can exploit distinct binding epitopes and leverage the intrinsic plasticity of the conformational landscape for the SARS-CoV-2 S protein to produce efficient neutralizing and mutation resistant characteristics. Structural and computational studies have also been instrumental in quantifying the structure, dynamics, and energetics of the SARS-CoV-2 spike protein binding with nanobodies. In this review, a comprehensive analysis of the current structural, biophysical, and computational biology investigations of SARS-CoV-2 S proteins and their complexes with distinct classes of nanobodies targeting different binding sites is presented. The analysis of computational studies is supplemented by an in-depth examination of mutational scanning simulations and identification of binding energy hotspots for distinct nanobody classes. The review is focused on the analysis of mechanisms underlying synergistic binding of multivalent nanobodies that can be superior to single nanobodies and conventional nanobody cocktails in combating escape mutations by effectively leveraging binding avidity and allosteric cooperativity. We discuss how structural insights and protein engineering approaches together with computational biology tools can aid in the rational design of synergistic combinations that exhibit superior binding and neutralization characteristics owing to avidity-mediated mechanisms.

Emergence of a recurrent insertion in the N-terminal domain of the SARS-CoV-2 spike glycoprotein

Tracking the evolution of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) through genomic surveillance programs is undoubtedly one of the key priorities in the current pandemic situation. Although the genome of SARS-CoV-2 acquires mutations at a slower rate compared with other RNA viruses, evolutionary pressures derived from the widespread circulation of SARS-CoV-2 in the human population have progressively favored the global emergence, though natural selection, of several variants of concern that carry multiple non-synonymous mutations in the spike glycoprotein. These are often placed in key sites within major antibody epitopes and may therefore confer resistance to neutralizing antibodies, leading to partial immune escape, or otherwise compensate infectivity deficits associated with other non-synonymous substitutions. As previously shown by other authors, several emerging variants carry recurrent deletion regions (RDRs) that display a partial overlap with antibody epitopes located in the spike N-terminal domain (NTD). Comparatively, very little attention had been directed towards spike insertion mutations prior to the emergence of the B.1.1.529 (omicron) lineage. This manuscript describes a single recurrent insertion region (RIR1) in the N-terminal domain of SARS-CoV-2 spike protein, characterized by at least 49 independent acquisitions of 1-8 additional codons between Val213 and Leu216 in different viral lineages. Even though RIR1 is unlikely to confer antibody escape, its association with two distinct formerly widespread lineages (A.2.5 and B.1.214.2), with the quickly spreading omicron and with other VOCs and VOIs warrants further investigation concerning its effects on spike structure and viral infectivity.

Allosteric Determinants of the SARS-CoV-2 Spike Protein Binding with Nanobodies: Examining Mechanisms of Mutational Escape and Sensitivity of the Omicron Variant

Structural and biochemical studies have recently revealed a range of rationally engineered nanobodies with efficient neutralizing capacity against the SARS-CoV-2 virus and resilience against mutational escape. In this study, we performed a comprehensive computational analysis of the SARS-CoV-2 spike trimer complexes with single nanobodies Nb6, VHH E, and complex with VHH E/VHH V nanobody combination. We combined coarse-grained and all-atom molecular simulations and collective dynamics analysis with binding free energy scanning, perturbation-response scanning, and network centrality analysis to examine mechanisms of nanobody-induced allosteric modulation and cooperativity in the SARS-CoV-2 spike trimer complexes with these nanobodies. By quantifying energetic and allosteric determinants of the SARS-CoV-2 spike protein binding with nanobodies, we also examined nanobody-induced modulation of escaping mutations and the effect of the Omicron variant on nanobody binding. The mutational scanning analysis supported the notion that E484A mutation can have a significant detrimental effect on nanobody binding and result in Omicron-induced escape from nanobody neutralization. Our findings showed that SARS-CoV-2 spike protein might exploit the plasticity of specific allosteric hotspots to generate escape mutants that alter response to binding without compromising activity. The network analysis supported these findings showing that VHH E/VHH V nanobody binding can induce long-range couplings between the cryptic binding epitope and ACE2-binding site through a broader ensemble of communication paths that is less dependent on specific mediating centers and therefore may be less sensitive to mutational perturbations of functional residues. The results suggest that binding affinity and long-range communications of the SARS-CoV-2 complexes with nanobodies can be determined by structurally stable regulatory centers and conformationally adaptable hotspots that are allosterically coupled and collectively control resilience to mutational escape.

Conformational Flexibility and Local Frustration in the Functional States of the SARS-CoV-2 Spike B.1.1.7 and B.1.351 Variants: Mutation-Induced Allosteric Modulation Mechanism of Functional Dynamics and Protein Stability

Structural and functional studies of the SARS-CoV-2 spike proteins have recently determined distinct functional states of the B.1.1.7 and B.1.351 spike variants, providing a molecular framework for understanding the mechanisms that link the effect of mutations with the enhanced virus infectivity and transmissibility. A detailed dynamic and energetic analysis of these variants was undertaken in the present work to quantify the effects of different mutations on functional conformational changes and stability of the SARS-CoV-2 spike protein. We employed the efficient and accurate coarse-grained (CG) simulations of multiple functional states of the D614G mutant, B.1.1.7 and B.1.351 spike variants to characterize conformational dynamics of the SARS-CoV-2 spike proteins and identify dynamic signatures of the functional regions that regulate transitions between the closed and open forms. By combining molecular simulations with full atomistic reconstruction of the trajectories and the ensemble-based mutational frustration analysis, we characterized how the intrinsic flexibility of specific spike regions can control functional conformational changes required for binding with the host-cell receptor. Using the residue-based mutational scanning of protein stability, we determined protein stability hotspots and identified potential energetic drivers favoring the receptor-accessible open spike states for the B.1.1.7 and B.1.351 spike variants. The results suggested that modulation of the energetic frustration at the inter-protomer interfaces can serve as a mechanism for allosteric couplings between mutational sites and the inter-protomer hinges of functional motions. The proposed mechanism of mutation-induced energetic frustration may result in greater adaptability and the emergence of multiple conformational states in the open form. This study suggested that SARS-CoV-2 B.1.1.7 and B.1.351 variants may leverage the intrinsic plasticity of functional regions in the spike protein for mutation-induced modulation of protein dynamics and allosteric regulation to control binding with the host cell receptor.

Emergence of a novel SARS-CoV-2 Pango lineage B.1.1.526 in West Bengal, India

BACKGROUND

Since its inception in late 2019, SARS-CoV-2 has been evolving continuously by procuring mutations, leading to emergence of numerous variants, causing second wave of pandemic in many countries including India in 2021. To control this pandemic continuous mutational surveillance and genomic epidemiology of circulating strains is very important to unveil the emergence of the novel variants and also monitor the evolution of existing variants.

METHODS

SARS-CoV-2 sequences were retrieved from GISAID database. Sequence alignment was performed with MAFT version 7. Phylogenetic tree was constructed by using MEGA (version X) and UShER.

RESULTS

In this study, we reported the emergence of a novel variant of SARS-CoV-2, named B.1.1.526, in India. This novel variant encompasses 129 SARS-CoV-2 strains which are characterized by the presence of 11 coexisting mutations including D614G, P681H, and V1230L in S glycoprotein. Out of these 129 sequences, 27 sequences also harbored E484K mutation in S glycoprotein. Phylogenetic analysis revealed strains of this novel variant emerged from the GR clade and formed a new cluster. Geographical distribution showed, out of 129 sequences, 126 were found in seven different states of India. Rest 3 sequences were observed in USA. Temporal analysis revealed this novel variant was first collected from Kolkata district of West Bengal, India.

CONCLUSIONS

The D614G, P618H and E484K mutations have previously been reported to favor increased transmissibility, enhanced infectivity, and immune invasion, respectively. The transmembrane domain (TM) of S2 subunit anchors S glycoprotein to the virus envelope. The V1230L mutation, present within the TM domain of S glycoprotein, might strengthen the interaction of S glycoprotein with the viral envelope and increase S glycoprotein deposition to the virion, resulting in more infectious virion. Therefore, the new variant having D614G, P618H, V1230L, and E484K may have higher infectivity, transmissibility, and immune invasion characteristics, and thus need to be monitored closely.

SARS-CoV-2 Variants and Their Relevant Mutational Profiles: Update Summer 2021

Since the beginning of the coronavirus disease 2019 (COVID-19) pandemic caused by it, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been undergoing a genetic diversification leading to the emergence of new variants. Nevertheless, a clear definition of the genetic signatures underlying the circulating variants is still missing. Here, we provide a comprehensive insight into mutational profiles characterizing each SARS-CoV-2 variant, focusing on spike mutations known to modulate viral infectivity and/or antigenicity. We focused on variants and on specific relevant mutations reported by GISAID, Nextstrain, Outbreak.info, Pango, and Stanford database websites that were associated with any clinical/diagnostic impact, according to published manuscripts. Furthermore, 1,223,338 full-length high-quality SARS-CoV-2 genome sequences were retrieved from GISAID and used to accurately define the specific mutational patterns in each variant. Finally, mutations were mapped on the three-dimensional structure of the SARS-CoV-2 spike protein to assess their localization in the different spike domains. Overall, this review sheds light and assists in defining the genetic signatures characterizing the currently circulating variants and their clinical relevance. IMPORTANCE Since the emergence of SARS-CoV-2, several recurrent mutations, particularly in the spike protein, arose during human-to-human transmission or spillover events between humans and animals, generating distinct worrisome variants of concern (VOCs) or of interest (VOIs), designated as such due to their clinical and diagnostic impacts. Characterizing these variants and their related mutations is important in tracking SAR-CoV-2 evolution and understanding the efficacy of vaccines and therapeutics based on monoclonal antibodies, convalescent-phase sera, and direct antivirals. Our study provides a comprehensive survey of the mutational profiles characterizing the important SARS-CoV-2 variants, focusing on spike mutations and highlighting other protein mutations.

RBD Double Mutations of SARS-CoV-2 Strains Increase Transmissibility through Enhanced Interaction between RBD and ACE2 Receptor

The COVID-19 pandemic, caused by SARS-CoV-2, has led to catastrophic damage for global human health. The initial step of SARS-CoV-2 infection is the binding of the receptor-binding domain (RBD) in its spike protein to the ACE2 receptor in the host cell membrane. Constant evolution of SARS-CoV-2 generates new mutations across its genome including the coding region for the RBD in the spike protein. In addition to the well-known single mutation in the RBD, the recent new mutation strains with an RBD "double mutation" are causing new outbreaks globally, as represented by the delta strain containing RBD L452R/T478K. Although it is considered that the increased transmissibility of double-mutated strains could be attributed to the altered interaction between the RBD and ACE2 receptor, the molecular details remain to be elucidated. Using the methods of molecular dynamics simulation, superimposed structural comparison, free binding energy estimation, and antibody escaping, we investigated the relationship between the ACE2 receptor and the RBD double mutants of L452R/T478K (delta), L452R/E484Q (kappa), and E484K/N501Y (beta, gamma). The results demonstrated that each of the three RBD double mutants altered the RBD structure and enhanced the binding of the mutated RBD to ACE2 receptor. Together with the mutations in other parts of the virus genome, the double mutations increase the transmissibility of SARS-CoV-2 to host cells.

Potential inhibitors of SARS-CoV-2 (COVID 19) spike protein of the delta and delta plus variant: In silico studies of medicinal plants of North-East India

Phytochemicals of 38 Medicinal plants of North-East India, with anti-viral, anti-oxidant or anti-bacterial properties were screened for properties of drug likeness. 231 phytochemicals were screened with LIPINSKI rule of five to obtain 131 candidates, which were further screened with SWISS-ADME, to obtain 50 phytochemicals. These phytochemicals were docked with the spike protein of the Delta variant (B.1.617.2) and Delta-Plus (AY.1) variant of SARS-CoV-2 using Autodock Vina and MOE 09. The target proteins were constructed by homology modeling using Swiss-Model. Hydroxychloroquine, taken as a standard in docking analysis, exhibited a binding energy of -6.5 ​kcal/mol and -6.1 ​kcal/mol with respect to the Delta variant and Delta-Plus variant respectively. Among the 50 docked results most flavones showed very good docking scores. 3,5,8-Trimethoxy-6,7,4,5-bis(methylenedioxy)flavone, a Poly-Methoxyflavone, produced a highest docking score of -8.7 ​kcal/mol with respect to both the spike protein targets. Poly-Methoxyflavones and Poly-Ethoxyflavones exhibited good binding affinity for the target spike protein of SARS-CoV-2, and can be potential anti-viral drug candidates against the existing Delta variant of the SARS-CoV-2.

A next generation sequencing (NGS) analysis to reveal genomic and proteomic mutation landscapes of SARS-CoV-2 in South Asia

•

48 SNPs were identified from the genome-wide analysis of 410 South Asian SARS-CoV-2 sequences.

•

About 85% SNPs are packaged in ORF1ab, spike protein, and nucleocapsid.

•

South Asian strains are highly related to the South American and European strains according to the phylogenetic analysis.

•

Unlike other countries, frequency of 1163A>T missense mutation is very high (78.80%) in Bangladeshi samples.

Counts for SARS-CoV-2 associated infections and fatalities are on the rise globally even in regions which contained the spread momentarily. The pattern of infections has been found to be controlled by the distinctive selection pressures exerted by fluctuating environmental nature and hosts. A total of 410 whole-genome sequences submitted by the South Asian countries were retrieved from the GISAID database and analyzed to assess the impact and pattern of mutations in this region. Most common and frequent mutations in the South Asian countries are 241C > T, 3037C > T, 14408C > T, and 23403A > G and about 85% SNPs are localized in ORF1ab, spike protein, and nucleocapsid. Among the identified mutations, the proportion of missense type (54.17%) was highest, followed by the synonymous (41.66%) and the non-coding types (4.17%). While analyzing transmission source in terms of geolocation, the largest clustered group from the South Asian countries was based on the G-clade (D614G) (81.7%; 335/410 samples), tracing the inception and transmission of SARS-CoV-2 infections in the South Asian countries from European regions. Phylogenetic analysis also revealed that the South Asian strains are highly related to the South American and European strains. We found that G-clade mutations are more prevalent (96.19%) in the samples of Bangladesh which were also prevalent in the European isolates. Surprisingly, one missense mutation (1163A > T) in ORF1ab gene became dominant only in Bangladesh (78.8%), which led to debates regarding effects on the pathogenicity and transmissibility of the virus. Overall, the findings of this study highlight the frequently mutated SARS-CoV-2 variants among the COVID-19 patients in the South Asian countries which might ease containment of the disease in this region through investigating the virulence reducing factors as the identified mutations are strongly correlated with low infection and mortality rate.