Structural and biochemical rationale for enhanced spike protein fitness in delta and kappa SARS-CoV-2 variants

Liquid-Electron Microscopy and the Real-Time Revolution

Advances in imaging technology enable striking views of life's most minute details. A missing piece of the puzzle, however, is the direct atomic observation of biomolecules in action. Liquid-phase transmission electron microscopy (liquid-EM) is the room-temperature correlate to cryo-electron microscopy, which is leading the resolution revolution in biophysics. This article reviews current challenges and opportunities in the liquid-EM field while discussing technical considerations for specimen enclosures, devices and systems, and scientific data management. Since liquid-EM is gaining traction in the life sciences community, cross talk among the disciplines of materials and life sciences is needed to disseminate knowledge of best practices along with high-level user engagement. How liquid-EM technology is inspiring the real-time revolution in molecular microscopy is also discussed. Looking ahead, the new movement can be better supported through open resource sharing and partnerships among academic, industry, and federal organizations, which may benefit from the scientific equity foundational to the technique.

A bivalent spike-targeting nanobody with anti-sarbecovirus activity

The continued emergence and zoonotic threat posed by coronaviruses highlight the urgent need for effective antiviral strategies with broad reactivity to counter new emerging strains. Nanobodies (or single-domain antibodies) are promising alternatives to traditional monoclonal antibodies, due to their small size, cost-effectiveness and ease of bioengineering. Here, we describe 7F, a llama-derived nanobody, targeting the spike receptor binding domain of sarbecoviruses and SARS-like coronaviruses. 7F demonstrates potent neutralization against SARS-CoV-2 and cross-neutralizing activity against SARS-CoV and SARS-like CoV WIV16 pseudoviruses. Structural analysis reveals 7F's ability to induce the formation of spike trimer dimers by engaging with two SARS-CoV-2 spike RBDs, targeting the highly conserved class IV region, though concentration dependent. Bivalent 7F constructs substantially enhance neutralization potency and breadth, up to more recent SARS-CoV-2 variants of concern. Furthermore, we demonstrate the therapeutic potential of bivalent 7F against SARS-CoV-2 in the fully differentiated 3D tissue cultures mirroring the epithelium of the human airway ex vivo. The broad sarbecovirus activity and distinctive structural features of bivalent 7F underscore its potential as promising antiviral against emerging and evolving sarbecoviruses.

Differential Impact of Spike Protein Mutations on SARS-CoV-2 Infectivity and Immune Evasion: Insights from Delta and Kappa Variants

SARS-CoV-2 continues to pose a global health challenge due to its high transmissibility and mutability, with new variants emerging that potentially undermine vaccination and therapeutic efforts. Mutations in the spike protein, particularly in the receptor-binding domain (RBD), significantly influence viral transmissibility and immune escape. However, the complex interplay of these mutations and their combined effects on viral fitness remain to be analyzed. In this study, we investigated the functional impact of key mutations found in the Delta and Kappa variants of SARS-CoV-2. Using pseudovirus assays, we demonstrated that the T478K and L452R mutations characteristic of the Delta variant primarily enhance viral infectivity, with minimal effect on antibody-mediated neutralization. Conversely, the E484Q mutation of the Kappa variant, alone or in combination with L452R, significantly improved evasion of antibody-mediated neutralization but appeared to compromise viral fitness and infectivity. Notably, contrary to previous reports, we found that the P681R mutation contributed neither to increased infectivity nor immune evasion at least in the assay system employed in this study. Our findings suggest that the Delta variant's global dominance over the Kappa variant may be attributed to its superior infectivity and transmissibility rather than enhanced immune evasion capabilities. These results provide valuable insights into the functional consequences of spike protein mutations and may aid in predicting the emergence and spread of future SARS-CoV-2 variants. Such understanding is crucial for enhancing public health preparedness and informing the development of next-generation vaccines and therapeutics.

Structure and dynamics of the interaction of Delta and Omicron BA.1 SARS-CoV-2 variants with REGN10987 Fab reveal mechanism of antibody action

Study of mechanisms by which antibodies recognize different viral strains is necessary for the development of new drugs and vaccines to treat COVID-19 and other infections. Here, we report 2.5 Å cryo-EM structure of the SARS-CoV-2 Delta trimeric S-protein in complex with Fab of the recombinant analog of REGN10987 neutralizing antibody. S-protein adopts "two RBD-down and one RBD-up" conformation. Fab interacts with RBDs in both conformations, blocking the recognition of angiotensin converting enzyme-2. Three-dimensional variability analysis reveals high mobility of the RBD/Fab regions. Interaction of REGN10987 with Wuhan, Delta, Omicron BA.1, and mutated variants of RBDs is analyzed by microscale thermophoresis, molecular dynamics simulations, and ΔG calculations with umbrella sampling and one-dimensional potential of mean force. Variability in molecular dynamics trajectories results in a large scatter of calculated ΔG values, but Boltzmann weighting provides an acceptable correlation with experiment. REGN10987 evasion of the Omicron variant is found to be due to the additive effect of the N440K and G446S mutations located at the RBD/Fab binding interface with a small effect of Q498R mutation. Our study explains the influence of known-to-date SARS-CoV-2 RBD mutations on REGN10987 recognition and highlights the importance of dynamics data beyond the static structure of the RBD/Fab complex.

Exploring distinct modes of inter-spike cross-linking for enhanced neutralization by SARS-CoV-2 antibodies

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and its Omicron subvariants drastically amplifies transmissibility, infectivity, and immune escape, mainly due to their resistance to most neutralizing antibodies. Thus, exploring the mechanisms underlying antibody evasion is crucial. Although the full-length native form of antibody, immunoglobulin G (IgG), offers valuable insights into the neutralization, structural investigations primarily focus on the fragment of antigen-binding (Fab). Here, we employ single-particle cryo-electron microscopy (cryo-EM) to characterize a W328-6H2 antibody, in its native IgG form complexed with severe acute respiratory syndrome (SARS), severe acute respiratory syndrome coronavirus 2 wild-type (WT) and Omicron variant BA.1 spike protein (S). Three high-resolution structures reveal that the full-length IgG forms a centered head-to-head dimer of trimer when binds fully stoichiometrically with both SARS and WT S, while adopting a distinct offset configuration with Omicron BA.1 S. Combined with functional assays, our results suggest that, beyond the binding affinity between the RBD epitope and Fab, the higher-order architectures of S trimer and full-length IgG play an additional role in neutralization, enriching our understanding of enhanced neutralization by SARS-CoV-2 antibodies.

Virion morphology and on-virus spike protein structures of diverse SARS-CoV-2 variants

The evolution of SARS-CoV-2 variants with increased fitness has been accompanied by structural changes in the spike (S) proteins, which are the major target for the adaptive immune response. Single-particle cryo-EM analysis of soluble S protein from SARS-CoV-2 variants has revealed this structural adaptation at high resolution. The analysis of S trimers in situ on intact virions has the potential to provide more functionally relevant insights into S structure and virion morphology. Here, we characterized B.1, Alpha, Beta, Gamma, Delta, Kappa, and Mu variants by cryo-electron microscopy and tomography, assessing S cleavage, virion morphology, S incorporation, "in-situ" high-resolution S structures, and the range of S conformational states. We found no evidence for adaptive changes in virion morphology, but describe multiple different positions in the S protein where amino acid changes alter local protein structure. Taken together, our data are consistent with a model where amino acid changes at multiple positions from the top to the base of the spike cause structural changes that can modulate the conformational dynamics of the S protein.

Dynamic expedition of leading mutations in SARS-CoV-2 spike glycoproteins

The continuous evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which caused the recent pandemic, has generated countless new variants with varying fitness. Mutations of the spike glycoprotein play a particularly vital role in shaping its evolutionary trajectory, as they have the capability to alter its infectivity and antigenicity. We present a time-resolved statistical method, Dynamic Expedition of Leading Mutations (deLemus), to analyze the evolutionary dynamics of the SARS-CoV-2 spike glycoprotein. The proposed L -index of the deLemus method is effective in quantifying the mutation strength of each amino acid site and outlining evolutionarily significant sites, allowing the comprehensive characterization of the evolutionary mutation pattern of the spike glycoprotein.

AI Promoted Virtual Screening, Structure-Based Hit Optimization, and Synthesis of Novel COVID-19 S-RBD Domain Inhibitors

Coronavirus disease 2019 (COVID-19) is caused by a new, highly pathogenic severe-acute-respiratory syndrome coronavirus 2 (SARS-CoV-2) that infects human cells through its transmembrane spike (S) glycoprotein. The receptor-binding domain (RBD) of the S protein interacts with the angiotensin-converting enzyme II (ACE2) receptor of the host cells. Therefore, pharmacological targeting of this interaction might prevent infection or spread of the virus. Here, we performed a virtual screening to identify small molecules that block S-ACE2 interaction. Large compound libraries were filtered for drug-like properties, promiscuity and protein-protein interaction-targeting ability based on their ADME-Tox descriptors and also to exclude pan-assay interfering compounds. A properly designed AI-based virtual screening pipeline was applied to the remaining compounds, comprising approximately 10% of the starting data sets, searching for molecules that could bind to the RBD of the S protein. All molecules were sorted according to their screening score, grouped based on their structure and postfiltered for possible interaction patterns with the ACE2 receptor, yielding 31 hits. These hit molecules were further tested for their inhibitory effect on Spike RBD/ACE2 (19-615) interaction. Six compounds inhibited the S-ACE2 interaction in a dose-dependent manner while two of them also prevented infection of human cells from a pseudotyped virus whose entry is mediated by the S protein of SARS-CoV-2. Of the two compounds, the benzimidazole derivative CKP-22 protected Vero E6 cells from infection with SARS-CoV-2, as well. Subsequent, hit-to-lead optimization of CKP-22 was effected through the synthesis of 29 new derivatives of which compound CKP-25 suppressed the Spike RBD/ACE2 (19-615) interaction, reduced the cytopathic effect of SARS-CoV-2 in Vero E6 cells (IC50 = 3.5 μM) and reduced the viral load in cell culture supernatants. Early in vitro ADME-Tox studies showed that CKP-25 does not possess biodegradation or liver metabolism issues, while isozyme-specific CYP450 experiments revealed that CKP-25 was a weak inhibitor of the CYP450 system. Moreover, CKP-25 does not elicit mutagenic effect on Escherichia coli WP2 uvrA strain. Thus, CKP-25 is considered a lead compound against COVID-19 infection.

The evaluation of SARS-CoV-2 mutations at the early stage of the pandemic in Istanbul population

BACKGROUND

Determination of SARS-CoV-2 variant is significant to prevent the spreads of COVID-19 disease.

METHODS

We aimed to evaluate the variants of SARS-CoV-2 rate in positive patients in Kanuni Sultan Suleyman Training and Research Hospital (KSS-TRH), Istanbul, Türkiye between 1st January and 30th November 2021 by using RT-PCR method.

RESULTS

Herein, 825,169 patients were evaluated (male:58.53% and female:41.47%) whether COVID-19 positive or not [( +):21.3% and (-):78.7%] and 175,367 patient was described as positive (53.2%-female and 46.8%-male) by RT-PCR. COVID-19 positive rate is observed highest in the 6-15- and 66-75-year age range. The frequencies were obtained as SARS-CoV-2 positive (without mutation of B.1.1.7 [B.1.1.7 (U.K), E484K, L452R, B.1.351 (S. Africa/Brazil) spike mutations] as 66.1% (n: 115,899), B.1.1.7 Variant as 23.2% (n:40,686), Delta mutation (L452R) variant as 9.8% (n:17,182), B.1.351 variant as 0.8% (n:1370) and E484K as 0.1% (n: 230). In April 2021, general SARS-CoV-2 and B.1.1.7 variant were dominantly observed. Up to July 2021, B.1.617.2 (Delta variant/ Indian variant) and E484K has been not observed. B.1.351 variant of SARS-CoV-2 has been started in February 2021 at the rarest ratio and March 2021 is the top point. September 2021 is the pick point of E484K. African/Brazil variant of SARS-CoV-2 has been started in February 2021 at the rarest ratio and March 2021 is the top point. September 2021 is the pick point of E484K. When the gender type is compared within the variants, women were found to be more prevalent in all varieties.

CONCLUSIONS

The meaning of these mutations is very important to understand the transmission capacity of the COVID-19 disease, pandemic episode, and diagnosis of the virus with mutation types. Understanding the variant type is important for monitoring herd immunity and the spread of the disease.

Comparison of the Neutralization Power of Sotrovimab Against SARS-CoV-2 Variants: Development of a Rapid Computational Method

BACKGROUND

The rapid evolution of SARS-CoV-2 imposed a huge challenge on disease control. Immune evasion caused by genetic variations of the SARS-CoV-2 spike protein's immunogenic epitopes affects the efficiency of monoclonal antibody-based therapy of COVID-19. Therefore, a rapid method is needed to evaluate the efficacy of the available monoclonal antibodies against the new emerging variants or potential novel variants.

OBJECTIVE

The aim of this study is to develop a rapid computational method to evaluate the neutralization power of anti-SARS-CoV-2 monoclonal antibodies against new SARS-CoV-2 variants and other potential new mutations.

METHODS

The amino acid sequence of the extracellular domain of the spike proteins of the severe acute respiratory syndrome coronavirus (GenBank accession number YP_009825051.1) and SARS-CoV-2 (GenBank accession number YP_009724390.1) were used to create computational 3D models for the native spike proteins. Specific mutations were introduced to the curated sequence to generate the different variant spike models. The neutralization potential of sotrovimab (S309) against these variants was evaluated based on its molecular interactions and Gibbs free energy in comparison to a reference model after molecular replacement of the reference receptor-binding domain with the variant's receptor-binding domain.

RESULTS

Our results show a loss in the binding affinity of the neutralizing antibody S309 with both SARS-CoV and SARS-CoV-2. The binding affinity of S309 was greater to the Alpha, Beta, Gamma, and Kappa variants than to the original Wuhan strain of SARS-CoV-2. However, S309 showed a substantially decreased binding affinity to the Delta and Omicron variants. Based on the mutational profile of Omicron subvariants, our data describe the effect of the G339H and G339D mutations and their role in escaping antibody neutralization, which is in line with published clinical reports.

CONCLUSIONS

This method is rapid, applicable, and of interest to adapt the use of therapeutic antibodies to the treatment of emerging variants. It could be applied to antibody-based treatment of other viral infections.

Structural basis for receptor-binding domain mobility of the spike in SARS-CoV-2 BA.2.86 and JN.1

Since 2019, SARS-CoV-2 has undergone mutations, resulting in pandemic and epidemic waves. The SARS-CoV-2 spike protein, crucial for cellular entry, binds to the ACE2 receptor exclusively when its receptor-binding domain (RBD) adopts the up-conformation. However, whether ACE2 also interacts with the RBD in the down-conformation to facilitate the conformational shift to RBD-up remains unclear. Herein, we present the structures of the BA.2.86 and the JN.1 spike proteins bound to ACE2. Notably, we successfully observed the ACE2-bound down-RBD, indicating an intermediate structure before the RBD-up conformation. The wider and mobile angle of RBDs in the up-state provides space for ACE2 to interact with the down-RBD, facilitating the transition to the RBD-up state. The K356T, but not N354-linked glycan, contributes to both of infectivity and neutralizing-antibody evasion in BA.2.86. These structural insights the spike-protein dynamics would help understand the mechanisms underlying SARS-CoV-2 infection and its neutralization.

Exploring conformational landscapes and binding mechanisms of convergent evolution for the SARS-CoV-2 spike Omicron variant complexes with the ACE2 receptor using AlphaFold2-based structural ensembles and molecular dynamics simulations

In this study, we combined AlphaFold-based approaches for atomistic modeling of multiple protein states and microsecond molecular simulations to accurately characterize conformational ensembles evolution and binding mechanisms of convergent evolution for the SARS-CoV-2 spike Omicron variants BA.1, BA.2, BA.2.75, BA.3, BA.4/BA.5 and BQ.1.1. We employed and validated several different adaptations of the AlphaFold methodology for modeling of conformational ensembles including the introduced randomized full sequence scanning for manipulation of sequence variations to systematically explore conformational dynamics of Omicron spike protein complexes with the ACE2 receptor. Microsecond atomistic molecular dynamics (MD) simulations provide a detailed characterization of the conformational landscapes and thermodynamic stability of the Omicron variant complexes. By integrating the predictions of conformational ensembles from different AlphaFold adaptations and applying statistical confidence metrics we can expand characterization of the conformational ensembles and identify functional protein conformations that determine the equilibrium dynamics for the Omicron spike complexes with the ACE2. Conformational ensembles of the Omicron RBD-ACE2 complexes obtained using AlphaFold-based approaches for modeling protein states and MD simulations are employed for accurate comparative prediction of the binding energetics revealing an excellent agreement with the experimental data. In particular, the results demonstrated that AlphaFold-generated extended conformational ensembles can produce accurate binding energies for the Omicron RBD-ACE2 complexes. The results of this study suggested complementarities and potential synergies between AlphaFold predictions of protein conformational ensembles and MD simulations showing that integrating information from both methods can potentially yield a more adequate characterization of the conformational landscapes for the Omicron RBD-ACE2 complexes. This study provides insights in the interplay between conformational dynamics and binding, showing that evolution of Omicron variants through acquisition of convergent mutational sites may leverage conformational adaptability and dynamic couplings between key binding energy hotspots to optimize ACE2 binding affinity and enable immune evasion.

Phylogeny and evolution of SARS-CoV-2 during Delta and Omicron variant waves in India

SARS-CoV-2 evolution has continued to generate variants, responsible for new pandemic waves locally and globally. Varying disease presentation and severity has been ascribed to inherent variant characteristics and vaccine immunity. This study analyzed genomic data from 305 whole genome sequences from SARS-CoV-2 patients before and through the third wave in India. Delta variant was reported in patients without comorbidity (97%), while Omicron BA.2 was reported in patients with comorbidity (77%). Tissue adaptation studies brought forth higher propensity of Omicron variants to bronchial tissue than lung, contrary to observation in Delta variants from Delhi. Study of codon usage pattern distinguished the prevalent variants, clustering them separately, Omicron BA.2 isolated in February grouped away from December strains, and all BA.2 after December acquired a new mutation S959P in ORF1b (44.3% of BA.2 in the study) indicating ongoing evolution. Loss of critical spike mutations in Omicron BA.2 and gain of immune evasion mutations including G142D, reported in Delta but absent in BA.1, and S371F instead of S371L in BA.1 could explain very brief period of BA.1 in December 2021, followed by complete replacement by BA.2. Higher propensity of Omicron variants to bronchial tissue, probably ensured increased transmission while Omicron BA.2 became the prevalent variant possibly due to evolutionary trade-off. Virus evolution continues to shape the epidemic and its culmination.Communicated by Ramaswamy H. Sarma.

AlphaFold2 Predictions of Conformational Ensembles and Atomistic Simulations of the SARS-CoV-2 Spike XBB Lineages Reveal Epistatic Couplings between Convergent Mutational Hotspots that Control ACE2 Affinity

In this study, we combined AlphaFold-based atomistic structural modeling, microsecond molecular simulations, mutational profiling, and network analysis to characterize binding mechanisms of the SARS-CoV-2 spike protein with the host receptor ACE2 for a series of Omicron XBB variants including XBB.1.5, XBB.1.5+L455F, XBB.1.5+F456L, and XBB.1.5+L455F+F456L. AlphaFold-based structural and dynamic modeling of SARS-CoV-2 Spike XBB lineages can accurately predict the experimental structures and characterize conformational ensembles of the spike protein complexes with the ACE2. Microsecond molecular dynamics simulations identified important differences in the conformational landscapes and equilibrium ensembles of the XBB variants, suggesting that combining AlphaFold predictions of multiple conformations with molecular dynamics simulations can provide a complementary approach for the characterization of functional protein states and binding mechanisms. Using the ensemble-based mutational profiling of protein residues and physics-based rigorous calculations of binding affinities, we identified binding energy hotspots and characterized the molecular basis underlying epistatic couplings between convergent mutational hotspots. Consistent with the experiments, the results revealed the mediating role of the Q493 hotspot in the synchronization of epistatic couplings between L455F and F456L mutations, providing a quantitative insight into the energetic determinants underlying binding differences between XBB lineages. We also proposed a network-based perturbation approach for mutational profiling of allosteric communications and uncovered the important relationships between allosteric centers mediating long-range communication and binding hotspots of epistatic couplings. The results of this study support a mechanism in which the binding mechanisms of the XBB variants may be determined by epistatic effects between convergent evolutionary hotspots that control ACE2 binding.

Cryo2StructData: A Large Labeled Cryo-EM Density Map Dataset for AI-based Modeling of Protein Structures

The advent of single-particle cryo-electron microscopy (cryo-EM) has brought forth a new era of structural biology, enabling the routine determination of large biological molecules and their complexes at atomic resolution. The high-resolution structures of biological macromolecules and their complexes significantly expedite biomedical research and drug discovery. However, automatically and accurately building atomic models from high-resolution cryo-EM density maps is still time-consuming and challenging when template-based models are unavailable. Artificial intelligence (AI) methods such as deep learning trained on limited amount of labeled cryo-EM density maps generate inaccurate atomic models. To address this issue, we created a dataset called Cryo2StructData consisting of 7,600 preprocessed cryo-EM density maps whose voxels are labelled according to their corresponding known atomic structures for training and testing AI methods to build atomic models from cryo-EM density maps. Cryo2StructData is larger than existing, publicly available datasets for training AI methods to build atomic protein structures from cryo-EM density maps. We trained and tested deep learning models on Cryo2StructData to validate its quality showing that it is ready for being used to train and test AI methods for building atomic models.

Ensemble-Based Mutational Profiling and Network Analysis of the SARS-CoV-2 Spike Omicron XBB Lineages for Interactions with the ACE2 Receptor and Antibodies: Cooperation of Binding Hotspots in Mediating Epistatic Couplings Underlies Binding Mechanism and Immune Escape

In this study, we performed a computational study of binding mechanisms for the SARS-CoV-2 spike Omicron XBB lineages with the host cell receptor ACE2 and a panel of diverse class one antibodies. The central objective of this investigation was to examine the molecular factors underlying epistatic couplings among convergent evolution hotspots that enable optimal balancing of ACE2 binding and antibody evasion for Omicron variants BA.1, BA2, BA.3, BA.4/BA.5, BQ.1.1, XBB.1, XBB.1.5, and XBB.1.5 + L455F/F456L. By combining evolutionary analysis, molecular dynamics simulations, and ensemble-based mutational scanning of spike protein residues in complexes with ACE2, we identified structural stability and binding affinity hotspots that are consistent with the results of biochemical studies. In agreement with the results of deep mutational scanning experiments, our quantitative analysis correctly reproduced strong and variant-specific epistatic effects in the XBB.1.5 and BA.2 variants. It was shown that Y453W and F456L mutations can enhance ACE2 binding when coupled with Q493 in XBB.1.5, while these mutations become destabilized when coupled with the R493 position in the BA.2 variant. The results provided a molecular rationale of the epistatic mechanism in Omicron variants, showing a central role of the Q493/R493 hotspot in modulating epistatic couplings between convergent mutational sites L455F and F456L in XBB lineages. The results of mutational scanning and binding analysis of the Omicron XBB spike variants with ACE2 receptors and a panel of class one antibodies provide a quantitative rationale for the experimental evidence that epistatic interactions of the physically proximal binding hotspots Y501, R498, Q493, L455F, and F456L can determine strong ACE2 binding, while convergent mutational sites F456L and F486P are instrumental in mediating broad antibody resistance. The study supports a mechanism in which the impact on ACE2 binding affinity is mediated through a small group of universal binding hotspots, while the effect of immune evasion could be more variant-dependent and modulated by convergent mutational sites in the conformationally adaptable spike regions.

Predicting Functional Conformational Ensembles and Binding Mechanisms of Convergent Evolution for SARS-CoV-2 Spike Omicron Variants Using AlphaFold2 Sequence Scanning Adaptations and Molecular Dynamics Simulations

In this study, we combined AlphaFold-based approaches for atomistic modeling of multiple protein states and microsecond molecular simulations to accurately characterize conformational ensembles and binding mechanisms of convergent evolution for the SARS-CoV-2 Spike Omicron variants BA.1, BA.2, BA.2.75, BA.3, BA.4/BA.5 and BQ.1.1. We employed and validated several different adaptations of the AlphaFold methodology for modeling of conformational ensembles including the introduced randomized full sequence scanning for manipulation of sequence variations to systematically explore conformational dynamics of Omicron Spike protein complexes with the ACE2 receptor. Microsecond atomistic molecular dynamic simulations provide a detailed characterization of the conformational landscapes and thermodynamic stability of the Omicron variant complexes. By integrating the predictions of conformational ensembles from different AlphaFold adaptations and applying statistical confidence metrics we can expand characterization of the conformational ensembles and identify functional protein conformations that determine the equilibrium dynamics for the Omicron Spike complexes with the ACE2. Conformational ensembles of the Omicron RBD-ACE2 complexes obtained using AlphaFold-based approaches for modeling protein states and molecular dynamics simulations are employed for accurate comparative prediction of the binding energetics revealing an excellent agreement with the experimental data. In particular, the results demonstrated that AlphaFold-generated extended conformational ensembles can produce accurate binding energies for the Omicron RBD-ACE2 complexes. The results of this study suggested complementarities and potential synergies between AlphaFold predictions of protein conformational ensembles and molecular dynamics simulations showing that integrating information from both methods can potentially yield a more adequate characterization of the conformational landscapes for the Omicron RBD-ACE2 complexes. This study provides insights in the interplay between conformational dynamics and binding, showing that evolution of Omicron variants through acquisition of convergent mutational sites may leverage conformational adaptability and dynamic couplings between key binding energy hotspots to optimize ACE2 binding affinity and enable immune evasion.

AlphaFold2-Enabled Atomistic Modeling of Structure, Conformational Ensembles, and Binding Energetics of the SARS-CoV-2 Omicron BA.2.86 Spike Protein with ACE2 Host Receptor and Antibodies: Compensatory Functional Effects of Binding Hotspots in Modulating Mechanisms of Receptor Binding and Immune Escape

The latest wave of SARS-CoV-2 Omicron variants displayed a growth advantage and increased viral fitness through convergent evolution of functional hotspots that work synchronously to balance fitness requirements for productive receptor binding and efficient immune evasion. In this study, we combined AlphaFold2-based structural modeling approaches with atomistic simulations and mutational profiling of binding energetics and stability for prediction and comprehensive analysis of the structure, dynamics, and binding of the SARS-CoV-2 Omicron BA.2.86 spike variant with ACE2 host receptor and distinct classes of antibodies. We adapted several AlphaFold2 approaches to predict both the structure and conformational ensembles of the Omicron BA.2.86 spike protein in the complex with the host receptor. The results showed that the AlphaFold2-predicted structural ensemble of the BA.2.86 spike protein complex with ACE2 can accurately capture the main conformational states of the Omicron variant. Complementary to AlphaFold2 structural predictions, microsecond molecular dynamics simulations reveal the details of the conformational landscape and produced equilibrium ensembles of the BA.2.86 structures that are used to perform mutational scanning of spike residues and characterize structural stability and binding energy hotspots. The ensemble-based mutational profiling of the receptor binding domain residues in the BA.2 and BA.2.86 spike complexes with ACE2 revealed a group of conserved hydrophobic hotspots and critical variant-specific contributions of the BA.2.86 convergent mutational hotspots R403K, F486P, and R493Q. To examine the immune evasion properties of BA.2.86 in atomistic detail, we performed structure-based mutational profiling of the spike protein binding interfaces with distinct classes of antibodies that displayed significantly reduced neutralization against the BA.2.86 variant. The results revealed the molecular basis of compensatory functional effects of the binding hotspots, showing that BA.2.86 lineage may have evolved to outcompete other Omicron subvariants by improving immune evasion while preserving binding affinity with ACE2 via through a compensatory effect of R493Q and F486P convergent mutational hotspots. This study demonstrated that an integrative approach combining AlphaFold2 predictions with complementary atomistic molecular dynamics simulations and robust ensemble-based mutational profiling of spike residues can enable accurate and comprehensive characterization of structure, dynamics, and binding mechanisms of newly emerging Omicron variants.

Altered receptor binding, antibody evasion and retention of T cell recognition by the SARS-CoV-2 XBB.1.5 spike protein

The XBB.1.5 variant of SARS-CoV-2 has rapidly achieved global dominance and exhibits a high growth advantage over previous variants. Preliminary reports suggest that the success of XBB.1.5 stems from mutations within its spike glycoprotein, causing immune evasion and enhanced receptor binding. We present receptor binding studies that demonstrate retention of binding contacts with the human ACE2 receptor and a striking decrease in binding to mouse ACE2 due to the revertant R493Q mutation. Despite extensive evasion of antibody binding, we highlight a region on the XBB.1.5 spike protein receptor binding domain (RBD) that is recognized by serum antibodies from a donor with hybrid immunity, collected prior to the emergence of the XBB.1.5 variant. T cell assays reveal high frequencies of XBB.1.5 spike-specific CD4+ and CD8+ T cells amongst donors with hybrid immunity, with the CD4+ T cells skewed towards a Th1 cell phenotype and having attenuated effector cytokine secretion as compared to ancestral spike protein-specific cells. Thus, while the XBB.1.5 variant has retained efficient human receptor binding and gained antigenic alterations, it remains susceptible to recognition by T cells induced via vaccination and previous infection.

The Diverse Nature of the Molecular Interactions That Govern the COV-2 Variants' Cell Receptor Affinity Ranking and Its Experimental Variability

A critical determinant of infectivity and virulence of the most infectious and or lethal variants of concern (VOCs): Wild Type, Delta and Omicron is related to the binding interactions between the receptor-binding domain of the spike and its host receptor, the initial step in cell infection. It is of the utmost importance to understand how mutations of a viral strain, especially those that are in the viral spike, affect the resulting infectivity of the emerging VOC, knowledge that could help us understand the variant virulence and inform the therapies applied or the vaccines developed. For this sake, we have applied a battery of computational protocols of increasing complexity to the calculation of the spike binding affinity for three variants of concern to the ACE2 cell receptor. The results clearly illustrate that the attachment of the spikes of the Delta and Omicron variants to the receptor originates through different molecular interaction mechanisms. All our protocols unanimously predict that the Delta variant has the highest receptor-binding affinity, while the Omicron variant displays a substantial variability in the binding affinity of the spike that relates to the structural plasticity of the Omicron spike-receptor complex. We suggest that the latter result could explain (at least in part) the variability of the in vitro binding results for this VOC and has led us to suggest a reason for the lower virulence of the Omicron variant as compared to earlier strains. Several hypotheses have been developed around this subject.

Cryo2StructData: A Large Labeled Cryo-EM Density Map Dataset for AI-based Modeling of Protein Structures

The advent of single-particle cryo-electron microscopy (cryo-EM) has brought forth a new era of structural biology, enabling the routine determination of large biological molecules and their complexes at atomic resolution. The high-resolution structures of biological macromolecules and their complexes significantly expedite biomedical research and drug discovery. However, automatically and accurately building atomic models from high-resolution cryo-EM density maps is still time-consuming and challenging when template-based models are unavailable. Artificial intelligence (AI) methods such as deep learning trained on limited amount of labeled cryo-EM density maps generate inaccurate atomic models. To address this issue, we created a dataset called Cryo2StructData consisting of 7,600 preprocessed cryo-EM density maps whose voxels are labelled according to their corresponding known atomic structures for training and testing AI methods to build atomic models from cryo-EM density maps. It is larger and of higher quality than any existing, publicly available dataset. We trained and tested deep learning models on Cryo2StructData to make sure it is ready for the large-scale development of AI methods for building atomic models from cryo-EM density maps.

ACE2-Fc and DPP4-Fc decoy receptors against SARS-CoV-2 and MERS-CoV variants: a quick therapeutic option for current and future coronaviruses outbreaks

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) and the Middle East respiratory syndrome coronavirus (MERS-CoV) are highly pathogenic human coronaviruses (CoVs). Anti-CoVs mAbs and vaccines may be effective, but the emergence of neutralization escape variants is inevitable. Angiotensin-converting enzyme 2 and dipeptidyl peptidase 4 enzyme are the getaway receptors for SARS-CoV-2 and MERS-CoV, respectively. Thus, we reformatted these receptors as Fc-fusion decoy receptors. Then, we tested them in parallel with anti-SARS-CoV (ab1-IgG) and anti-MERS-CoV (M336-IgG) mAbs against several variants using pseudovirus neutralization assay. The generated Fc-based decoy receptors exhibited a strong inhibitory effect against all pseudotyped CoVs. Results showed that although mAbs can be effective antiviral drugs, they might rapidly lose their efficacy against highly mutated viruses. We suggest that receptor traps can be engineered as Fc-fusion proteins for highly mutating viruses with known entry receptors, for a faster and effective therapeutic response even against virus harboring antibodies escape mutations.

AlphaFold2-Enabled Atomistic Modeling of Epistatic Binding Mechanisms for the SARS-CoV-2 Spike Omicron XBB.1.5, EG.5 and FLip Variants: Convergent Evolution Hotspots Cooperate to Control Stability and Conformational Adaptability in Balancing ACE2 Binding and Antibody Resistance

In this study, we combined AI-based atomistic structural modeling and microsecond molecular simulations of the SARS-CoV-2 Spike complexes with the host receptor ACE2 for XBB.1.5+L455F, XBB.1.5+F456L(EG.5) and XBB.1.5+L455F/F456L (FLip) lineages to examine the mechanisms underlying the role of convergent evolution hotspots in balancing ACE2 binding and antibody evasion. Using the ensemble-based mutational scanning of the spike protein residues and physics-based rigorous computations of binding affinities, we identified binding energy hotspots and characterized molecular basis underlying epistatic couplings between convergent mutational hotspots. Consistent with the experiments, the results revealed the mediating role of Q493 hotspot in synchronization of epistatic couplings between L455F and F456L mutations providing a quantitative insight into the mechanism underlying differences between XBB lineages. Mutational profiling is combined with network-based model of epistatic couplings showing that the Q493, L455 and F456 sites mediate stable communities at the binding interface with ACE2 and can serve as stable mediators of non-additive couplings. Structure-based mutational analysis of Spike protein binding with the class 1 antibodies quantified the critical role of F456L and F486P mutations in eliciting strong immune evasion response. The results of this analysis support a mechanism in which the emergence of EG.5 and FLip variants may have been dictated by leveraging strong epistatic effects between several convergent revolutionary hotspots that provide synergy between the improved ACE2 binding and broad neutralization resistance. This interpretation is consistent with the notion that functionally balanced substitutions which simultaneously optimize immune evasion and high ACE2 affinity may continue to emerge through lineages with beneficial pair or triplet combinations of RBD mutations involving mediators of epistatic couplings and sites in highly adaptable RBD regions.

Accurate Characterization of Conformational Ensembles and Binding Mechanisms of the SARS-CoV-2 Omicron BA.2 and BA.2.86 Spike Protein with the Host Receptor and Distinct Classes of Antibodies Using AlphaFold2-Augmented Integrative Computational Modeling

The latest wave SARS-CoV-2 Omicron variants displayed a growth advantage and the increased viral fitness through convergent evolution of functional hotspots that work synchronously to balance fitness requirements for productive receptor binding and efficient immune evasion. In this study, we combined AlphaFold2-based structural modeling approaches with all-atom MD simulations and mutational profiling of binding energetics and stability for prediction and comprehensive analysis of the structure, dynamics, and binding of the SARS-CoV-2 Omicron BA.2.86 spike variant with ACE2 host receptor and distinct classes of antibodies. We adapted several AlphaFold2 approaches to predict both structure and conformational ensembles of the Omicron BA.2.86 spike protein in the complex with the host receptor. The results showed that AlphaFold2-predicted conformational ensemble of the BA.2.86 spike protein complex can accurately capture the main dynamics signatures obtained from microscond molecular dynamics simulations. The ensemble-based dynamic mutational scanning of the receptor binding domain residues in the BA.2 and BA.2.86 spike complexes with ACE2 dissected the role of the BA.2 and BA.2.86 backgrounds in modulating binding free energy changes revealing a group of conserved hydrophobic hotspots and critical variant-specific contributions of the BA.2.86 mutational sites R403K, F486P and R493Q. To examine immune evasion properties of BA.2.86 in atomistic detail, we performed large scale structure-based mutational profiling of the S protein binding interfaces with distinct classes of antibodies that displayed significantly reduced neutralization against BA.2.86 variant. The results quantified specific function of the BA.2.86 mutations to ensure broad resistance against different classes of RBD antibodies. This study revealed the molecular basis of compensatory functional effects of the binding hotspots, showing that BA.2.86 lineage may have primarily evolved to improve immune escape while modulating binding affinity with ACE2 through cooperative effect of R403K, F486P and R493Q mutations. The study supports a hypothesis that the impact of the increased ACE2 binding affinity on viral fitness is more universal and is mediated through cross-talk between convergent mutational hotspots, while the effect of immune evasion could be more variant-dependent.

Structural understanding of SARS-CoV-2 virus entry to host cells

Coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is a major global health concern associated with millions of fatalities worldwide. Mutant variants of the virus have further exacerbated COVID-19 mortality and infection rates, emphasizing the urgent need for effective preventive strategies. Understanding the viral infection mechanism is crucial for developing therapeutics and vaccines. The entry of SARS-CoV-2 into host cells is a key step in the infection pathway and has been targeted for drug development. Despite numerous reviews of COVID-19 and the virus, there is a lack of comprehensive reviews focusing on the structural aspects of viral entry. In this review, we analyze structural changes in Spike proteins during the entry process, dividing the entry process into prebinding, receptor binding, proteolytic cleavage, and membrane fusion steps. By understanding the atomic-scale details of viral entry, we can better target the entry step for intervention strategies. We also examine the impacts of mutations in Spike proteins, including the Omicron variant, on viral entry. Structural information provides insights into the effects of mutations and can guide the development of therapeutics and vaccines. Finally, we discuss available structure-based approaches for the development of therapeutics and vaccines. Overall, this review provides a detailed analysis of the structural aspects of SARS-CoV-2 viral entry, highlighting its significance in the development of therapeutics and vaccines against COVID-19. Therefore, our review emphasizes the importance of structural information in combating SARS-CoV-2 infection.

Infection and Vaccine Induced Spike Antibody Responses Against SARS-CoV-2 Variants of Concern in COVID-19-Naïve Children and Adults

Although a more efficient adaptive humoral immune response has been proposed to underlie the usually favorable outcome of pediatric COVID-19, the breadth of viral and vaccine cross-reactivity toward the ever-mutating Spike protein among variants of concern (VOCs) has not yet been compared between children and adults. We assessed antibodies to conformational Spike in COVID-19-naïve children and adults vaccinated by BNT162b2 and ChAdOx1, and naturally infected with SARS-CoV-2 Early Clade, Delta, and Omicron. Sera were analyzed against Spike including naturally occurring VOCs Alpha, Beta, Gamma, Delta, and Omicron BA.1, BA.2, BA.5, BQ.1.1, BA2.75.2, and XBB.1, and variants of interest Epsilon, Kappa, Eta, D.2, and artificial mutant Spikes. There was no notable difference between breadth and longevity of antibody against VOCs in children and adults. Vaccinated individuals displayed similar immunoreactivity profiles across variants compared with naturally infected individuals. Delta-infected patients had an enhanced cross-reactivity toward Delta and earlier VOCs compared to patients infected by Early Clade SARS-CoV-2. Although Omicron BA.1, BA.2, BA.5, BQ.1.1, BA2.75.2, and XBB.1 antibody titers were generated after Omicron infection, cross-reactive binding against Omicron subvariants was reduced across all infection, immunization, and age groups. Some mutations, such as 498R and 501Y, epistatically combined to enhance cross-reactive binding, but could not fully compensate for antibody-evasive mutations within the Omicron subvariants tested. Our results reveal important molecular features central to the generation of high antibody titers and broad immunoreactivity that should be considered in future vaccine design and global serosurveillance in the context of limited vaccine boosters available to the pediatric population.

Comparative Analysis of Conformational Dynamics and Systematic Characterization of Cryptic Pockets in the SARS-CoV-2 Omicron BA.2, BA.2.75 and XBB.1 Spike Complexes with the ACE2 Host Receptor: Confluence of Binding and Structural Plasticity in Mediating Networks of Conserved Allosteric Sites

In the current study, we explore coarse-grained simulations and atomistic molecular dynamics together with binding energetics scanning and cryptic pocket detection in a comparative examination of conformational landscapes and systematic characterization of allosteric binding sites in the SARS-CoV-2 Omicron BA.2, BA.2.75 and XBB.1 spike full-length trimer complexes with the host receptor ACE2. Microsecond simulations, Markov state models and mutational scanning of binding energies of the SARS-CoV-2 BA.2 and BA.2.75 receptor binding domain complexes revealed the increased thermodynamic stabilization of the BA.2.75 variant and significant dynamic differences between these Omicron variants. Molecular simulations of the SARS-CoV-2 Omicron spike full-length trimer complexes with the ACE2 receptor complemented atomistic studies and enabled an in-depth analysis of mutational and binding effects on conformational dynamic and functional adaptability of the Omicron variants. Despite considerable structural similarities, Omicron variants BA.2, BA.2.75 and XBB.1 can induce unique conformational dynamic signatures and specific distributions of the conformational states. Using conformational ensembles of the SARS-CoV-2 Omicron spike trimer complexes with ACE2, we conducted a comprehensive cryptic pocket screening to examine the role of Omicron mutations and ACE2 binding on the distribution and functional mechanisms of the emerging allosteric binding sites. This analysis captured all experimentally known allosteric sites and discovered networks of inter-connected and functionally relevant allosteric sites that are governed by variant-sensitive conformational adaptability of the SARS-CoV-2 spike structures. The results detailed how ACE2 binding and Omicron mutations in the BA.2, BA.2.75 and XBB.1 spike complexes modulate the distribution of conserved and druggable allosteric pockets harboring functionally important regions. The results are significant for understanding the functional roles of druggable cryptic pockets that can be used for allostery-mediated therapeutic intervention targeting conformational states of the Omicron variants.

Cryo-EM and cryo-ET of the spike, virion, and antibody neutralization of SARS-CoV-2 and VOCs

Since the outbreak of the COVID-19 pandemic, cryo-electron microscopy (cryo-EM) and cryo-electron tomography (cryo-ET) have been demonstrated to be powerful and efficient tools in structural studies of distinct conformational states of the trimeric spike protein of SARS-CoV-2 and the VOCs as well as the intact virion. Cryo-EM has also contributed greatly to revealing the molecular basis of receptor recognition and antibody neutralization of the S trimer. Additionally, it has provided structural insights into the enhanced transformation and immune evasion of the VOCs, thus facilitating antiviral antibody and drug discovery. In this review, we summarize the contributions of cryo-EM and cryo-ET in revealing the structures of SARS-CoV-2 S trimer and intact virion and the mechanisms of receptor binding and antibody neutralization. We also highlight their prospective utilities in the development of vaccines and future therapeutics against emerging SARS-CoV-2 variants and other epidemic viruses.

Exploring Conformational Landscapes and Cryptic Binding Pockets in Distinct Functional States of the SARS-CoV-2 Omicron BA.1 and BA.2 Trimers: Mutation-Induced Modulation of Protein Dynamics and Network-Guided Prediction of Variant-Specific Allosteric Binding Sites

A significant body of experimental structures of SARS-CoV-2 spike trimers for the BA.1 and BA.2 variants revealed a considerable plasticity of the spike protein and the emergence of druggable binding pockets. Understanding the interplay of conformational dynamics changes induced by the Omicron variants and the identification of cryptic dynamic binding pockets in the S protein is of paramount importance as exploring broad-spectrum antiviral agents to combat the emerging variants is imperative. In the current study, we explore conformational landscapes and characterize the universe of binding pockets in multiple open and closed functional spike states of the BA.1 and BA.2 Omicron variants. By using a combination of atomistic simulations, a dynamics network analysis, and an allostery-guided network screening of binding pockets in the conformational ensembles of the BA.1 and BA.2 spike conformations, we identified all experimentally known allosteric sites and discovered significant variant-specific differences in the distribution of binding sites in the BA.1 and BA.2 trimers. This study provided a structural characterization of the predicted cryptic pockets and captured the experimentally known allosteric sites, revealing the critical role of conformational plasticity in modulating the distribution and cross-talk between functional binding sites. We found that mutational and dynamic changes in the BA.1 variant can induce the remodeling and stabilization of a known druggable pocket in the N-terminal domain, while this pocket is drastically altered and may no longer be available for ligand binding in the BA.2 variant. Our results predicted the experimentally known allosteric site in the receptor-binding domain that remains stable and ranks as the most favorable site in the conformational ensembles of the BA.2 variant but could become fragmented and less probable in BA.1 conformations. We also uncovered several cryptic pockets formed at the inter-domain and inter-protomer interface, including functional regions of the S2 subunit and stem helix region, which are consistent with the known role of pocket residues in modulating conformational transitions and antibody recognition. The results of this study are particularly significant for understanding the dynamic and network features of the universe of available binding pockets in spike proteins, as well as the effects of the Omicron-variant-specific modulation of preferential druggable pockets. The exploration of predicted druggable sites can present a new and previously underappreciated opportunity for therapeutic interventions for Omicron variants through the conformation-selective and variant-specific targeting of functional sites involved in allosteric changes.

From Alpha to Omicron: How Different Variants of Concern of the SARS-Coronavirus-2 Impacted the World

SARS-CoV-2, the virus that causes COVID-19, is prone to mutations and the generation of genetic variants. Since its first outbreak in 2019, SARS-CoV-2 has continually evolved, resulting in the emergence of several lineages and variants of concern (VOC) that have gained more efficient transmission, severity, and immune evasion properties. The World Health Organization has given these variants names according to the letters of the Greek Alphabet, starting with the Alpha (B.1.1.7) variant, which emerged in 2020, followed by the Beta (B.1.351), Gamma (P.1), Delta (B.1.617.2), and Omicron (B.1.1.529) variants. This review explores the genetic variation among different VOCs of SARS-CoV-2 and how the emergence of variants made a global impact on the pandemic.

Examining Functional Linkages Between Conformational Dynamics, Protein Stability and Evolution of Cryptic Binding Pockets in the SARS-CoV-2 Omicron Spike Complexes with the ACE2 Host Receptor: Recombinant Omicron Variants Mediate Variability of Conserved Allosteric Sites and Binding Epitopes

In the current study, we explore coarse-grained simulations and atomistic molecular dynamics together with binding energetics scanning and cryptic pocket detection in a comparative examination of conformational landscapes and systematic characterization of allosteric binding sites in the SARS-CoV-2 Omicron BA.2, BA.2.75 and XBB.1 spike full-length trimer complexes with the host receptor ACE2. Microsecond simulations, Markov state models and mutational scanning of binding energies of the SARS-CoV-2 BA.2 and BA.2.75 receptor binding domain complexes revealed the increased thermodynamic stabilization of the BA.2.75 variant and significant dynamic differences between these Omicron variants. Molecular simulations of the SARS-CoV-2 Omicron spike full length trimer complexes with the ACE2 receptor complemented atomistic studies and enabled an in-depth analysis of mutational and binding effects on conformational dynamic and functional adaptability of the Omicron variants. Despite considerable structural similarities, Omicron variants BA.2, BA.2.75 and XBB.1 can induce unique conformational dynamic signatures and specific distributions of the conformational states. Using conformational ensembles of the SARS-CoV-2 Omicron spike trimer complexes with ACE2, we conducted a comprehensive cryptic pocket screening to examine the role of Omicron mutations and ACE2 binding on the distribution and functional mechanisms of the emerging allosteric binding sites. This analysis captured all experimentally known allosteric sites and discovered networks of inter-connected and functionally relevant allosteric sites that are governed by variant-sensitive conformational adaptability of the SARS-CoV-2 spike structures. The results detailed how ACE2 binding and Omicron mutations in the BA.2, BA.2.75 and XBB.1 spike complexes modulate the distribution of conserved and druggable allosteric pockets harboring functionally important regions. The results of are significant for understanding functional roles of druggable cryptic pockets that can be used for allostery-mediated therapeutic intervention targeting conformational states of the Omicron variants.

Knowledge and attitude regarding the COVID-19 pandemic among undergraduate health science students of Nepal: An online survey

Objectives

This study aims to assess health science students' knowledge and attitude about COVID-19 epidemiology, management, and prevention; and the association of knowledge and attitude with various sociodemographic characteristics.

Methods

An online survey was done among 524 undergraduate health science students using a pre-tested questionnaire across 19 health science institutions in Nepal from 30 June to 11 August 2021. All subjects were enrolled in the study after informed consent. Outcomes were Knowledge level, attitude level, and predictors of knowledge level and attitude level. Bivariate analysis was done to determine the association between variables.

Results

Of 524 students, 42.9% were male, and 57.1% were female. More than half (54.6%) and the majority (85.1%) participants had good knowledge and attitude, respectively. Approximately three-fifths (59.4%) of the participants were from the B. P. Koirala Institute of Health Sciences (BPKIHS). Students in BPKIHS (odds ratio = 1.774; 95% confidence interval = 1.243-2.533), junior years (odds ratio = 8.892; 95% confidence interval = 5.814-13.599), age less than 23 years (odds ratio = 2.985; 95% confidence interval = 2.089-4.266) were more likely to have good knowledge. Students under 23 years (odds ratio = 24.160; 95% confidence interval: 9.570-60.992) and those in junior years (odds ratio = 4.460; 95% confidence interval = 3.753-5.300) were likely to have a good attitude level. Students in BPKIHS (odds ratio = 0.443; 95% confidence interval = 0.272-0.722) were less likely to have a good attitude.

Conclusions

Overall, health science students had adequate knowledge and a good attitude regarding COVID-19. However, students lacked knowledge regarding infectiousness, transmission, post-vaccination observation period, remdesivir use, convalescent plasma therapy, and awake-prone positioning. Knowledge and attitude scores were associated with age, stream, and study institution.

Markov State Models and Perturbation-Based Approaches Reveal Distinct Dynamic Signatures and Hidden Allosteric Pockets in the Emerging SARS-Cov-2 Spike Omicron Variant Complexes with the Host Receptor: The Interplay of Dynamics and Convergent Evolution Modulates Allostery and Functional Mechanisms

The new generation of SARS-CoV-2 Omicron variants displayed a significant growth advantage and increased viral fitness by acquiring convergent mutations, suggesting that the immune pressure can promote convergent evolution leading to the sudden acceleration of SARS-CoV-2 evolution. In the current study, we combined structural modeling, microsecond molecular dynamics simulations, and Markov state models to characterize conformational landscapes and identify specific dynamic signatures of the SARS-CoV-2 spike complexes with the host receptor ACE2 for the recently emerged highly transmissible XBB.1, XBB.1.5, BQ.1, and BQ.1.1 Omicron variants. Microsecond simulations and Markovian modeling provided a detailed characterization of the functional conformational states and revealed the increased thermodynamic stabilization of the XBB.1.5 subvariant, which can be contrasted to more dynamic BQ.1 and BQ.1.1 subvariants. Despite considerable structural similarities, Omicron mutations can induce unique dynamic signatures and specific distributions of the conformational states. The results suggested that variant-specific changes of the conformational mobility in the functional interfacial loops of the receptor-binding domain in the SARS-CoV-2 spike protein can be fine-tuned through crosstalk between convergent mutations which could provide an evolutionary path for modulation of immune escape. By combining atomistic simulations and Markovian modeling analysis with perturbation-based approaches, we determined important complementary roles of convergent mutation sites as effectors and receivers of allosteric signaling involved in modulation of conformational plasticity and regulation of allosteric communications. This study also revealed hidden allosteric pockets and suggested that convergent mutation sites could control evolution and distribution of allosteric pockets through modulation of conformational plasticity in the flexible adaptable regions.

Probing conformational landscapes of binding and allostery in the SARS-CoV-2 omicron variant complexes using microsecond atomistic simulations and perturbation-based profiling approaches: hidden role of omicron mutations as modulators of allosteric signaling and epistatic relationships

In this study, we systematically examine the conformational dynamics, binding and allosteric communications in the Omicron BA.1, BA.2, BA.3 and BA.4/BA.5 spike protein complexes with the ACE2 host receptor using molecular dynamics simulations and perturbation-based network profiling approaches. Microsecond atomistic simulations provided a detailed characterization of the conformational landscapes and revealed the increased thermodynamic stabilization of the BA.2 variant which can be contrasted with the BA.4/BA.5 variants inducing a significant mobility of the complexes. Using the dynamics-based mutational scanning of spike residues, we identified structural stability and binding affinity hotspots in the Omicron complexes. Perturbation response scanning and network-based mutational profiling approaches probed the effect of the Omicron mutations on allosteric interactions and communications in the complexes. The results of this analysis revealed specific roles of Omicron mutations as conformationally plastic and evolutionary adaptable modulators of binding and allostery which are coupled to the major regulatory positions through interaction networks. Through perturbation network scanning of allosteric residue potentials in the Omicron variant complexes performed in the background of the original strain, we characterized regions of epistatic couplings that are centered around the binding affinity hotspots N501Y and Q498R. Our results dissected the vital role of these epistatic centers in regulating protein stability, efficient ACE2 binding and allostery which allows for accumulation of multiple Omicron immune escape mutations at other sites. Through integrative computational approaches, this study provides a systematic analysis of the effects of Omicron mutations on thermodynamics, binding and allosteric signaling in the complexes with ACE2 receptor.

Mechanism and evolution of human ACE2 binding by SARS-CoV-2 spike

Spike glycoprotein of SARS-CoV-2 mediates viral entry into host cells by facilitating virus attachment and membrane fusion. ACE2 is the main receptor of SARS-CoV-2 and its interaction with spike has shaped the virus' emergence from an animal reservoir and subsequent evolution in the human host. Many structural studies on the spike:ACE2 interaction have provided insights into mechanisms driving viral evolution during the on-going pandemic. This review describes the molecular basis of spike binding to ACE2, outlines mechanisms that have optimised this interaction during viral evolution, and suggests directions for future research.

Free Energy Perturbation Calculations of Mutation Effects on SARS-CoV-2 RBD::ACE2 Binding Affinity

The strength of binding between human angiotensin converting enzyme 2 (ACE2) and the receptor binding domain (RBD) of viral spike protein plays a role in the transmissibility of the SARS-CoV-2 virus. In this study we focus on a subset of RBD mutations that have been frequently observed in infected individuals and probe binding affinity changes to ACE2 using surface plasmon resonance (SPR) measurements and free energy perturbation (FEP) calculations. Our SPR results are largely in accord with previous studies but discrepancies do arise due to differences in experimental methods and to protocol differences even when a single method is used. Overall, we find that FEP performance is superior to that of other computational approaches examined as determined by agreement with experiment and, in particular, by its ability to identify stabilizing mutations. Moreover, the calculations successfully predict the observed cooperative stabilization of binding by the Q498R N501Y double mutant present in Omicron variants and offer a physical explanation for the underlying mechanism. Overall, our results suggest that despite the significant computational cost, FEP calculations may offer an effective strategy to understand the effects of interfacial mutations on protein-protein binding affinities and, hence, in a variety of practical applications such as the optimization of neutralizing antibodies.

Emergence of COVID-19 Variants: An Update

Severe acute respiratory disease virus-2 (SARS CoV-2) is one of the deadliest global threats faced by mankind to date. Despite the colossal efforts, the viral pandemic swept across all boundaries. Besides the virulence and susceptible population, the low proofreading capacity and error-prone mechanism of RNA-dependent RNA polymerase (RdRp) have contributed to new variants and reinfections. The World Health Organization has officially categorized these variants as variants of concern or variants of interest. This nomenclature is not merely to suffice the surveillance but also to have effective treatment and vaccine options in place. Coronavirus disease 2019 (COVID-19) variants have the propensity to render the available treatment strategies futile owing to the mutations they acquire. The futility of treatment strategies can be attributed either to the ineffectiveness or the shortage of supply given the skyrocketing increase in the number of cases. Presently, the Omicron variant is the most widespread one and is known to escape the protection, be it immune-derived, vaccination-derived, or hybrid. WHO has recommended modification in vaccine development policies and few companies have introduced Omicron-adapted vaccine jabs. Keeping in view the unending tale of COVID-19 variants and the huge data available on the same, this review focuses on providing insight into the emergence and ongoing dynamics of these new COVID-19 variants.

Site specific N- and O-glycosylation mapping of the spike proteins of SARS-CoV-2 variants of concern

The glycosylation on the spike (S) protein of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes COVID-19, modulates the viral infection by altering conformational dynamics, receptor interaction and host immune responses. Several variants of concern (VOCs) of SARS-CoV-2 have evolved during the pandemic, and crucial mutations on the S protein of the virus have led to increased transmissibility and immune escape. In this study, we compare the site-specific glycosylation and overall glycomic profiles of the wild type Wuhan-Hu-1 strain (WT) S protein and five VOCs of SARS-CoV-2: Alpha, Beta, Gamma, Delta and Omicron. Interestingly, both N- and O-glycosylation sites on the S protein are highly conserved among the spike mutant variants, particularly at the sites on the receptor-binding domain (RBD). The conservation of glycosylation sites is noteworthy, as over 2 million SARS-CoV-2 S protein sequences have been reported with various amino acid mutations. Our detailed profiling of the glycosylation at each of the individual sites of the S protein across the variants revealed intriguing possible association of glycosylation pattern on the variants and their previously reported infectivity. While the sites are conserved, we observed changes in the N- and O-glycosylation profile across the variants. The newly emerged variants, which showed higher resistance to neutralizing antibodies and vaccines, displayed a decrease in the overall abundance of complex-type glycans with both fucosylation and sialylation and an increase in the oligomannose-type glycans across the sites. Among the variants, the glycosylation sites with significant changes in glycan profile were observed at both the N-terminal domain and RBD of S protein, with Omicron showing the highest deviation. The increase in oligomannose-type happens sequentially from Alpha through Delta. Interestingly, Omicron does not contain more oligomannose-type glycans compared to Delta but does contain more compared to the WT and other VOCs. O-glycosylation at the RBD showed lower occupancy in the VOCs in comparison to the WT. Our study on the sites and pattern of glycosylation on the SARS-CoV-2 S proteins across the VOCs may help to understand how the virus evolved to trick the host immune system. Our study also highlights how the SARS-CoV-2 virus has conserved both N- and O- glycosylation sites on the S protein of the most successful variants even after undergoing extensive mutations, suggesting a correlation between infectivity/ transmissibility and glycosylation.

Markov State Models and Perturbation-Based Approaches Reveal Distinct Dynamic Signatures and Hidden Allosteric Pockets in the Emerging SARS-Cov-2 Spike Omicron Variants Complexes with the Host Receptor: The Interplay of Dynamics and Convergent Evolution Modulates Allostery and Functional Mechanisms

The new generation of SARS-CoV-2 Omicron variants displayed a significant growth advantage and the increased viral fitness by acquiring convergent mutations, suggesting that the immune pressure can promote convergent evolution leading to the sudden acceleration of SARS-CoV-2 evolution. In the current study, we combined structural modeling, extensive microsecond MD simulations and Markov state models to characterize conformational landscapes and identify specific dynamic signatures of the SARS-CoV-2 spike complexes with the host receptor ACE2 for the recently emerged highly transmissible XBB.1, XBB.1.5, BQ.1, and BQ.1.1 Omicron variants. Microsecond simulations and Markovian modeling provided a detailed characterization of the conformational landscapes and revealed the increased thermodynamic stabilization of the XBB.1.5 subvariant which is contrasted to more dynamic BQ.1 and BQ.1.1 subvariants. Despite considerable structural similarities, Omicron mutations can induce unique dynamic signatures and specific distributions of conformational states. The results suggested that variant-specific changes of conformational mobility in the functional interfacial loops of the spike receptor binding domain can be fine-tuned through cross-talk between convergent mutations thereby providing an evolutionary path for modulation of immune escape. By combining atomistic simulations and Markovian modeling analysis with perturbation-based approaches, we determined important complementary roles of convergent mutation sites as effectors and receivers of allosteric signaling involved in modulating conformational plasticity at the binding interface and regulating allosteric responses. This study also characterized the dynamics-induced evolution of allosteric pockets in the Omicron complexes that revealed hidden allosteric pockets and suggested that convergent mutation sites could control evolution and distribution of allosteric pockets through modulation of conformational plasticity in the flexible adaptable regions. Through integrative computational approaches, this investigation provides a systematic analysis and comparison of the effects of Omicron subvariants on conformational dynamics and allosteric signaling in the complexes with the ACE2 receptor.

For Table of Contents Use Only

Balancing Functional Tradeoffs between Protein Stability and ACE2 Binding in the SARS-CoV-2 Omicron BA.2, BA.2.75 and XBB Lineages: Dynamics-Based Network Models Reveal Epistatic Effects Modulating Compensatory Dynamic and Energetic Changes

Evolutionary and functional studies suggested that the emergence of the Omicron variants can be determined by multiple fitness trade-offs including the immune escape, binding affinity for ACE2, conformational plasticity, protein stability and allosteric modulation. In this study, we systematically characterize conformational dynamics, structural stability and binding affinities of the SARS-CoV-2 Spike Omicron complexes with the host receptor ACE2 for BA.2, BA.2.75, XBB.1 and XBB.1.5 variants. We combined multiscale molecular simulations and dynamic analysis of allosteric interactions together with the ensemble-based mutational scanning of the protein residues and network modeling of epistatic interactions. This multifaceted computational study characterized molecular mechanisms and identified energetic hotspots that can mediate the predicted increased stability and the enhanced binding affinity of the BA.2.75 and XBB.1.5 complexes. The results suggested a mechanism driven by the stability hotspots and a spatially localized group of the Omicron binding affinity centers, while allowing for functionally beneficial neutral Omicron mutations in other binding interface positions. A network-based community model for the analysis of epistatic contributions in the Omicron complexes is proposed revealing the key role of the binding hotspots R498 and Y501 in mediating community-based epistatic couplings with other Omicron sites and allowing for compensatory dynamics and binding energetic changes. The results also showed that mutations in the convergent evolutionary hotspot F486 can modulate not only local interactions but also rewire the global network of local communities in this region allowing the F486P mutation to restore both the stability and binding affinity of the XBB.1.5 variant which may explain the growth advantages over the XBB.1 variant. The results of this study are consistent with a broad range of functional studies rationalizing functional roles of the Omicron mutation sites that form a coordinated network of hotspots enabling a balance of multiple fitness tradeoffs and shaping up a complex functional landscape of virus transmissibility.

Probing Conformational Landscapes of Binding and Allostery in the SARS-CoV-2 Omicron Variant Complexes Using Microsecond Atomistic Simulations and Perturbation-Based Profiling Approaches: Hidden Role of Omicron Mutations as Modulators of Allosteric Signaling and Epistatic Relationships

In this study, we systematically examine the conformational dynamics, binding and allosteric communications in the Omicron BA.1, BA.2, BA.3 and BA.4/BA.5 complexes with the ACE2 host receptor using molecular dynamics simulations and perturbation-based network profiling approaches. Microsecond atomistic simulations provided a detailed characterization of the conformational landscapes and revealed the increased thermodynamic stabilization of the BA.2 variant which is contrasted with the BA.4/BA.5 variants inducing a significant mobility of the complexes. Using ensemble-based mutational scanning of binding interactions, we identified binding affinity and structural stability hotspots in the Omicron complexes. Perturbation response scanning and network-based mutational profiling approaches probed the effect of the Omicron variants on allosteric communications. The results of this analysis revealed specific roles of Omicron mutations as "plastic and evolutionary adaptable" modulators of binding and allostery which are coupled to the major regulatory positions through interaction networks. Through perturbation network scanning of allosteric residue potentials in the Omicron variant complexes, which is performed in the background of the original strain, we identified that the key Omicron binding affinity hotspots N501Y and Q498R could mediate allosteric interactions and epistatic couplings. Our results suggested that the synergistic role of these hotspots in controlling stability, binding and allostery can enable for compensatory balance of fitness tradeoffs with conformationally and evolutionary adaptable immune-escape Omicron mutations. Through integrative computational approaches, this study provides a systematic analysis of the effects of Omicron mutations on thermodynamics, binding and allosteric signaling in the complexes with ACE2 receptor. The findings support a mechanism in which Omicron mutations can evolve to balance thermodynamic stability and conformational adaptability in order to ensure proper tradeoff between stability, binding and immune escape.

Coarse-Grained Molecular Simulations and Ensemble-Based Mutational Profiling of Protein Stability in the Different Functional Forms of the SARS-CoV-2 Spike Trimers: Balancing Stability and Adaptability in BA.1, BA.2 and BA.2.75 Variants

Evolutionary and functional studies have suggested that the emergence of Omicron variants can be determined by multiple fitness tradeoffs including immune escape, binding affinity, conformational plasticity, protein stability, and allosteric modulation. In this study, we embarked on a systematic comparative analysis of the conformational dynamics, electrostatics, protein stability, and allostery in the different functional states of spike trimers for BA.1, BA.2, and BA.2.75 variants. Using efficient and accurate coarse-grained simulations and atomistic reconstruction of the ensembles, we examined the conformational dynamics of the spike trimers that agree with the recent functional studies, suggesting that BA.2.75 trimers are the most stable among these variants. A systematic mutational scanning of the inter-protomer interfaces in the spike trimers revealed a group of conserved structural stability hotspots that play a key role in the modulation of functional dynamics and are also involved in the inter-protomer couplings through local contacts and interaction networks with the Omicron mutational sites. The results of mutational scanning provided evidence that BA.2.75 trimers are more stable than BA.2 and comparable in stability to the BA.1 variant. Using dynamic network modeling of the S Omicron BA.1, BA.2, and BA.2.75 trimers, we showed that the key network mediators of allosteric interactions are associated with the major stability hotspots that are interconnected along potential communication pathways. The network analysis of the BA.1, BA.2, and BA.2.75 trimers suggested that the increased thermodynamic stability of the BA.2.75 variant may be linked with the organization and modularity of the residue interaction network that allows for allosteric communications between structural stability hotspots and Omicron mutational sites. This study provided a plausible rationale for a mechanism in which Omicron mutations may evolve by targeting vulnerable sites of conformational adaptability to elicit immune escape while maintaining their control on balancing protein stability and functional fitness through robust allosteric communications with the stability hotspots.

Epitope mapping of SARS-CoV-2 spike protein differentiates the antibody binding activity in vaccinated and infected individuals

Previous studies have attempted to characterize the antibody response of individuals to the SARS-CoV-2 virus on a linear peptide level by utilizing peptide microarrays. These studies have helped to identify epitopes that have potential to be used for diagnostic tests to identify infected individuals. The immunological responses of individuals who have received the two most popular vaccines available in the US, the Moderna mRNA-1273 or the Pfizer BNT162b2 mRNA vaccines, have not been characterized. We aimed to identify linear peptides of the SARS-CoV-2 spike protein that elicited high IgG or IgA binding activity and to compare the immunoreactivity of infected individuals to those who received both doses of either vaccine by utilizing peptide microarrays. Our results revealed peptide epitopes of significant IgG binding among recently infected individuals. Some of these peptides are located near variable regions of the receptor binding domains as well as the conserved region in the c-terminal of the spike protein implicated in the high infectivity of SARS-CoV-2. Vaccinated individuals lacked a response to these distinct markers despite the overall antibody binding activity being similar.

Animal Models, Zoonotic Reservoirs, and Cross-Species Transmission of Emerging Human-Infecting Coronaviruses

Over the past three decades, coronavirus (CoV) diseases have impacted humans more than any other emerging infectious disease. The recent emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19 (coronavirus disease 2019), has resulted in huge economic disruptions and loss of human lives. The SARS-CoV-2 genome was found to mutate more rapidly due to sustained transmission in humans and potentially animals, resulting in variants of concern (VOCs) that threaten global human health. However, the primary difficulties are filling in the current knowledge gaps in terms of the origin and modalities of emergence for these viruses. Because many CoVs threatening human health are suspected to have a zoonotic origin, identifying the animal hosts implicated in the spillover or spillback events would be beneficial for current pandemic management and to prevent future outbreaks. In this review, wesummarize the animal models, zoonotic reservoirs, and cross-species transmission of the emerging human CoVs. Finally, we comment on potential sources of SARS-CoV-2 Omicron VOCs and the new SARS-CoV-2 recombinants currently under investigation.

SARS-CoV-2 variant evasion of monoclonal antibodies based on in vitro studies

Monoclonal antibodies (mAbs) offer a treatment option for individuals with severe COVID-19 and are especially important in high-risk individuals where vaccination is not an option. Given the importance of understanding the evolution of resistance to mAbs by SARS-CoV-2, we reviewed the available in vitro neutralization data for mAbs against live variants and viral constructs containing spike mutations of interest. Unfortunately, evasion of mAb-induced protection is being reported with new SARS-CoV-2 variants. The magnitude of neutralization reduction varied greatly among mAb-variant pairs. For example, sotrovimab retained its neutralization capacity against Omicron BA.1 but showed reduced efficacy against BA.2, BA.4 and BA.5, and BA.2.12.1. At present, only bebtelovimab has been reported to retain its efficacy against all SARS-CoV-2 variants considered here. Resistance to mAb neutralization was dominated by the action of epitope single amino acid substitutions in the spike protein. Although not all observed epitope mutations result in increased mAb evasion, amino acid substitutions at non-epitope positions and combinations of mutations also contribute to evasion of neutralization. This Review highlights the implications for the rational design of viral genomic surveillance and factors to consider for the development of novel mAb therapies.

Structural analysis of receptor engagement and antigenic drift within the BA.2 spike protein

The BA.2 sub-lineage of the Omicron (B.1.1.529) severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant rapidly supplanted the original BA.1 sub-lineage in early 2022. Both lineages threatened the efficacy of vaccine-elicited antibodies and acquired increased binding to several mammalian ACE2 receptors. Cryoelectron microscopy (cryo-EM) analysis of the BA.2 spike (S) glycoprotein in complex with mouse ACE2 (mACE2) identifies BA.1- and BA.2-mutated residues Q493R, N501Y, and Y505H as complementing non-conserved residues between human and mouse ACE2, rationalizing the enhanced S protein-mACE2 interaction for Omicron variants. Cryo-EM structures of the BA.2 S-human ACE2 complex and of the extensively mutated BA.2 amino-terminal domain (NTD) reveal a dramatic reorganization of the highly antigenic N1 loop into a β-strand, providing an explanation for decreased binding of the BA.2 S protein to antibodies isolated from BA.1-convalescent patients. Our analysis reveals structural mechanisms underlying the antigenic drift in the rapidly evolving Omicron variant landscape.

Valproate-coenzyme A conjugate blocks opening of receptor binding domains in the spike trimer of SARS-CoV-2 through an allosteric mechanism

The receptor-binding domains (RBDs) of the SARS-CoV-2 spike trimer exhibit "up" and "down" conformations often targeted by neutralizing antibodies. Only in the "up" configuration can RBDs bind to the ACE2 receptor of the host cell and initiate the process of viral multiplication. Here, we identify a lead compound (3-oxo-valproate-coenzyme A conjugate or Val-CoA) that stabilizes the spike trimer with RBDs in the down conformation. Val-CoA interacts with three R408 residues, one from each RBD, which significantly reduces the inter-subunit R408-R408 distance by ∼ 13 Å and closes the central pore formed by the three RBDs. Experimental evidence is presented that R408 is part of a triggering mechanism that controls the prefusion to postfusion state transition of the spike trimer. By stabilizing the RBDs in the down configuration, this and other related compounds can likely attenuate viral transmission. The reported findings for binding of Val-CoA to the spike trimer suggest a new approach for the design of allosteric antiviral drugs that do not have to compete for specific virus-receptor interactions but instead hinder the conformational motion of viral membrane proteins essential for interaction with the host cell. Here, we introduce an approach to target the spike protein by identifying lead compounds that stabilize the RBDs in the trimeric "down" configuration. When these compounds trimerize monomeric RBD immunogens as co-immunogens, they could also induce new types of non-ACE2 blocking antibodies that prevent local cell-to-cell transmission of the virus, providing a novel approach for inhibition of SARS-CoV-2.

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.

Effect of Double Mutation (L452R and E484Q) on the Binding Affinity of Monoclonal Antibodies (mAbs) against the RBD-A Target for Vaccine Development

The COVID-19 pandemic, caused by SARS-CoV-2, emerges as a global health problem, as the viral genome is evolving rapidly to form several variants. Advancement and progress in the development of effective vaccines and neutralizing monoclonal antibodies are promising to combat viral infections. In the current scenario, several lineages containing "co-mutations" in the receptor-binding domain (RBD) region of the spike (S) protein are imposing new challenges. Co-occurrence of some co-mutations includes delta (L452R/T478K), kappa (L452R/E484Q), and a common mutation in both beta and gamma variants (E484K/N501Y). The effect of co-mutants (L452R/E484Q) on human angiotensin-converting enzyme 2 (hACE2) binding has already been elucidated. Here, for the first time, we investigated the role of these RBD co-mutations (L452R/E484Q) on the binding affinity of mAbs by adopting molecular dynamics (MD) simulation and free-energy binding estimation. The results obtained from our study suggest that the structural and dynamic changes introduced by these co-mutations reduce the binding affinity of the viral S protein to monoclonal antibodies (mAbs). The structural changes imposed by L452R create a charged patch near the interfacial surface that alters the affinity towards mAbs. In E484Q mutation, polar negatively charged E484 helps in the formation of electrostatic interaction, while the neutrally charged Q residue affects the interaction by forming repulsive forces. MD simulations along with molecular mechanics-generalized Born surface area (MMGBSA) studies revealed that the REGN 10933, BD-368-2, and S2M11 complexes have reduced binding affinity towards the double-mutant RBD. This indicates that their mutant (MT) structures have a stronger ability to escape from most antibodies than the wild type (WT). However, EY6A Ab showed higher affinity towards the double MT-RBD complex as compared to the WT. However, no significant effect of the per-residue contribution of double-mutated residues was observed, as this mAb does not interact with the region harboring L452 and E484 residues.

Electron microscopy holdings of the Protein Data Bank: the impact of the resolution revolution, new validation tools, and implications for the future

As a discipline, structural biology has been transformed by the three-dimensional electron microscopy (3DEM) "Resolution Revolution" made possible by convergence of robust cryo-preservation of vitrified biological materials, sample handling systems, and measurement stages operating a liquid nitrogen temperature, improvements in electron optics that preserve phase information at the atomic level, direct electron detectors (DEDs), high-speed computing with graphics processing units, and rapid advances in data acquisition and processing software. 3DEM structure information (atomic coordinates and related metadata) are archived in the open-access Protein Data Bank (PDB), which currently holds more than 11,000 3DEM structures of proteins and nucleic acids, and their complexes with one another and small-molecule ligands (~ 6% of the archive). Underlying experimental data (3DEM density maps and related metadata) are stored in the Electron Microscopy Data Bank (EMDB), which currently holds more than 21,000 3DEM density maps. After describing the history of the PDB and the Worldwide Protein Data Bank (wwPDB) partnership, which jointly manages both the PDB and EMDB archives, this review examines the origins of the resolution revolution and analyzes its impact on structural biology viewed through the lens of PDB holdings. Six areas of focus exemplifying the impact of 3DEM across the biosciences are discussed in detail (icosahedral viruses, ribosomes, integral membrane proteins, SARS-CoV-2 spike proteins, cryogenic electron tomography, and integrative structure determination combining 3DEM with complementary biophysical measurement techniques), followed by a review of 3DEM structure validation by the wwPDB that underscores the importance of community engagement.