Molecular insights into receptor binding energetics and neutralization of SARS-CoV-2 variants

GōMartini 3: From large conformational changes in proteins to environmental bias corrections

Coarse-grained modeling has become an important tool to supplement experimental measurements, allowing access to spatio-temporal scales beyond all-atom based approaches. The GōMartini model combines structure- and physics-based coarse-grained approaches, balancing computational efficiency and accurate representation of protein dynamics with the capabilities of studying proteins in different biological environments. This paper introduces an enhanced GōMartini model, which combines a virtual-site implementation of Gō models with Martini 3. The implementation has been extensively tested by the community since the release of the reparametrized version of Martini. This work demonstrates the capabilities of the model in diverse case studies, ranging from protein-membrane binding to protein-ligand interactions and AFM force profile calculations. The model is also versatile, as it can address recent inaccuracies reported in the Martini protein model. Lastly, the paper discusses the advantages, limitations, and future perspectives of the Martini 3 protein model and its combination with Gō models.

Conformational and Stability Analysis of SARS-CoV-2 Spike Protein Variants by Molecular Simulation

We performed a comprehensive structural analysis of the conformational space of several spike (S) protein variants using molecular dynamics (MD) simulations. Specifically, we examined four well-known variants (Delta, BA.1, XBB.1.5, and JN.1) alongside the wild-type (WT) form of SARS-CoV-2. The conformational states of each variant were characterized by analyzing their distributions within a selected space of collective variables (CVs), such as inter-domain distances between the receptor-binding domain (RBD) and the N-terminal domain (NTD). Our primary focus was to identify conformational states relevant to potential structural transitions and to determine the set of native contacts (NCs) that stabilize these conformations. The results reveal that genetically more distant variants, such as XBB.1.5, BA.1, and JN.1, tend to adopt more compact conformational states compared to the WT. Additionally, these variants exhibit novel NC profiles, characterized by an increased number of specific contacts distributed among ionic, polar, and nonpolar residues. We further analyzed the impact of specific mutations, including T478K, N500Y, and Y504H. These mutations not only enhance interactions with the human host receptor but also alter inter-chain stability by introducing additional NCs compared to the WT. Consequently, these mutations may influence the accessibility of certain protein regions to neutralizing antibodies. Overall, these findings contribute to a deeper understanding of the structural and functional variations among S protein variants.

Variant-Specific Interactions at the Plasma Membrane: Heparan Sulfate's Impact on SARS-CoV-2 Binding Kinetics

The spread of SARS-CoV-2 led to the emergence of several variants of concern (VOCs). The spike glycoprotein, responsible for engaging the viral receptor, exhibits the highest density of mutations, suggesting an ongoing evolution to optimize viral entry. This study characterizes the bond formed by virion mimics carrying the SARS-CoV-2 spike protein and the plasma membrane of host cells in the early stages of virus entry. Contrary to the traditional analysis of isolated ligand-receptor pairs, we utilized well-defined biomimetic models and biochemical and biophysical techniques to characterize the multivalent interaction of VOCs with the complex cell membrane. We observed an overall increase in the binding affinity for newer VOCs. By progressively reducing the system complexity, we identify heparan sulfate (HS) as a main driver of this variation, with a 10-fold increase in affinity for Omicron BA.1 over that of the original strain. These results demonstrate the essential role of coreceptors, particularly HS, in the modulation of SARS-CoV-2 infection and highlight the importance of multiscale biophysical and biochemical assays that account for membrane complexity to fully characterize and understand the role of molecular components and their synergy in viral attachment and entry.

From viral assembly to host interaction: AFM's contributions to virology

Viruses represent a diverse pool of obligate parasites that infect virtually every known organism, as such, their study is incredibly valuable for a range of fields including public health, medicine, agriculture, and ecology, and the development of biomedical technologies. Having evolved over millions of years, each virus has a unique and often complicated biology, that must be characterized on a case-by-case basis, even between strains of the same taxon. Owing to its nanoscale spatial resolution, atomic force microscopy (AFM) represents a powerful tool for exploring virus biology, including structural features, kinetics of binding to host cell ligands, virion self-assembly, and budding behaviors. Through the availability of numerous chemistries and advances in imaging modes, AFM is able to explore the complex web of host-virus interactions and life-cycle at a single virus level, exploring features at the level of individual bonds and molecules. Due to the wide array of techniques developed and data analysis approaches available, AFM can provide information that cannot be furnished by other modalities, especially at a single virus level. Here, we highlight the unique methods and information that can be obtained through the use of AFM, demonstrating both its utility and versatility in the study of viruses. As the technology continues to rapidly evolve, AFM is likely to remain an integral part of research, providing unique and important insight into many aspects of virology.

Atomic Force Microscopy: Mechanosensor and Mechanotransducer for Probing Biological System from Molecules to Tissues

Atomic Force Microscopy (AFM) is a powerful technique with widespread applications in various scientific fields, including biology. It operates by precisely detecting the interaction between a sharp tip and a sample surface, providing high-resolution topographical information and mechanical properties at a nanoscale. Through the years, a deeper understanding of this tip-sample interaction and the mechanisms by which it can be more precisely regulated have invariably led to improvements in AFM imaging. Additionally, AFM can serve not only as a sensor but also as a tool for actively manipulating the mechanical properties of biological systems. By applying controlled forces to the sample surface, AFM allows for a deeper understanding of mechanotransduction pathways, the intricate signaling cascades that convert physical cues into biochemical responses. This review, is an extensive overview of the current status of AFM working either as a mechanosensor or a mechanotransducer to probe biological systems across diverse scales, from individual molecules to entire tissues is presented. Challenges are discussed and potential future research directions are elaborated.

Probing living cell dynamics and molecular interactions using atomic force microscopy

Atomic force microscopy (AFM) has emerged as a powerful tool for studying biological interactions at the single-molecule level, offering unparalleled insights into receptor-ligand dynamics on living cells. This review discusses key developments in the application of AFM, highlighting its ability to capture nanomechanical properties of cellular surfaces and probe dynamic interactions, such as virus-host binding. AFM's versatility in measuring mechanical forces and mapping molecular interactions in near-physiological conditions is explored. The review also emphasizes how AFM provides critical insights into cell surface organization, receptor functionality, and viral entry mechanisms, advancing the understanding of cellular and molecular processes.

Nanomechanical footprint of SARS-CoV-2 variants in complex with a potent nanobody by molecular simulations

Rational design of novel antibody therapeutics against viral infections such as coronavirus relies on surface complementarity and high affinity for their effectiveness. Here, we explore an additional property of protein complexes, the intrinsic mechanical stability, in SARS-CoV-2 variants when complexed with a potent antibody. In this study, we utilized a recent implementation of the GōMartini 3 approach to investigate large conformational changes in protein complexes with a focus on the mechanostability of the receptor-binding domain (RBD) from WT, Alpha, Delta, and XBB.1.5 variants in complex with the H11-H4 nanobody. The analysis revealed moderate differences in mechanical stability among these variants. Also, we identified crucial residues in both the RBD and certain protein segments in the nanobody that contribute to this property. By performing pulling simulations and monitoring the presence of specific native and non-native contacts across the protein complex interface, we provided mechanistic insights into the dissociation process. Force-displacement profiles indicate a tensile force clamp mechanism associated with the type of protein complex. Our computational approach not only highlights the key mechanostable interactions that are necessary to maintain overall stability, but it also paves the way for the rational design of potent antibodies that are mechanostable and effective against emergent SARS-CoV-2 variants.

Widespread exposure to SARS-CoV-2 in wildlife communities

Pervasive SARS-CoV-2 infections in humans have led to multiple transmission events to animals. While SARS-CoV-2 has a potential broad wildlife host range, most documented infections have been in captive animals and a single wildlife species, the white-tailed deer. The full extent of SARS-CoV-2 exposure among wildlife communities and the factors that influence wildlife transmission risk remain unknown. We sampled 23 species of wildlife for SARS-CoV-2 and examined the effects of urbanization and human use on seropositivity. Here, we document positive detections of SARS-CoV-2 RNA in six species, including the deer mouse, Virginia opossum, raccoon, groundhog, Eastern cottontail, and Eastern red bat between May 2022-September 2023 across Virginia and Washington, D.C., USA. In addition, we found that sites with high human activity had three times higher seroprevalence than low human-use areas. We obtained SARS-CoV-2 genomic sequences from nine individuals of six species which were assigned to seven Pango lineages of the Omicron variant. The close match to variants circulating in humans at the time suggests at least seven recent human-to-animal transmission events. Our data support that exposure to SARS-CoV-2 has been widespread in wildlife communities and suggests that areas with high human activity may serve as points of contact for cross-species transmission.

Micro-second time-resolved X-ray single-molecule internal motions of SARS-CoV-2 spike variants

Single-molecule intramolecular dynamics were successfully measured for three variants of SARS-CoV-2 spike protein, alpha: B.1.1.7, delta: B.1.617, and omicron: B.1.1.529, with a time resolution of 100 μs using X-rays. The results were then compared with respect to the magnitude and directions of motions for the three variants. The largest 3-D intramolecular movement was observed for the omicron variant irrespective of ACE2 receptor binding. A more detailed analysis of the intramolecular motions revealed that the distribution state of intramolecular motion for the three variants was completely different with and without ACE2 receptor binding. The molecular dynamics for the trimeric spike protein of the omicron variant increased when ACE2 binding occurred. At that time, the diffusion constant increased from 71.0 [mrad2/ms] to 91.1 [mrad2/ms].

Single-Molecule Investigation of the Binding Interface Stability of SARS-CoV-2 Variants with ACE2

The SARS-CoV-2 pandemic spurred numerous research endeavors to comprehend the virus and mitigate its global severity. Understanding the binding interface between the virus and human receptors is pivotal to these efforts and paramount to curbing infection and transmission. Here we employ atomic force microscopy and steered molecular dynamics simulation to explore SARS-CoV-2 receptor binding domain (RBD) variants and angiotensin-converting enzyme 2 (ACE2), examining the impact of mutations at key residues upon binding affinity. Our results show that the Omicron and Delta variants possess strengthened binding affinity in comparison to the Mu variant. Further, using sera from individuals either vaccinated or with acquired immunity following Delta strain infection, we assess the impact of immunity upon variant RBD/ACE2 complex formation. Single-molecule force spectroscopy analysis suggests that vaccination before infection may provide stronger protection across variants. These results underscore the need to monitor antigenic changes in order to continue developing innovative and effective SARS-CoV-2 abrogation strategies.

Single-molecule force stability of the SARS-CoV-2-ACE2 interface in variants-of-concern

Mutations in SARS-CoV-2 have shown effective evasion of population immunity and increased affinity to the cellular receptor angiotensin-converting enzyme 2 (ACE2). However, in the dynamic environment of the respiratory tract, forces act on the binding partners, which raises the question of whether not only affinity but also force stability of the SARS-CoV-2-ACE2 interaction might be a selection factor for mutations. Using magnetic tweezers, we investigate the impact of amino acid substitutions in variants of concern (Alpha, Beta, Gamma and Delta) and on force-stability and bond kinetic of the receptor-binding domain-ACE2 interface at a single-molecule resolution. We find a higher affinity for all of the variants of concern (>fivefold) compared with the wild type. In contrast, Alpha is the only variant of concern that shows higher force stability (by 17%) compared with the wild type. Using molecular dynamics simulations, we rationalize the mechanistic molecular origins of this increase in force stability. Our study emphasizes the diversity of contributions to the transmissibility of variants and establishes force stability as one of the several factors for fitness. Understanding fitness advantages opens the possibility for the prediction of probable mutations, allowing a rapid adjustment of therapeutics, vaccines and intervention measures.

Molecular-Scale Investigations Reveal the Effect of Natural Polyphenols on BAX/Bcl-2 Interactions

Apoptosis signaling controls the cell cycle through the protein-protein interactions (PPIs) of its major B-cell lymphoma 2-associated x protein (BAX) and B-cell lymphoma 2 protein (Bcl-2). Due to the antagonistic function of both proteins, apoptosis depends on a properly tuned balance of the kinetics of BAX and Bcl-2 activities. The utilization of natural polyphenols to regulate the binding process of PPIs is feasible. However, the mechanism of this modulation has not been studied in detail. Here, we utilized atomic force microscopy (AFM) to evaluate the effects of polyphenols (kaempferol, quercetin, dihydromyricetin, baicalin, curcumin, rutin, epigallocatechin gallate, and gossypol) on the BAX/Bcl-2 binding mechanism. We demonstrated at the molecular scale that polyphenols quantitatively affect the interaction forces, kinetics, thermodynamics, and structural properties of BAX/Bcl-2 complex formation. We observed that rutin, epigallocatechin gallate, and baicalin reduced the binding affinity of BAX/Bcl-2 by an order of magnitude. Combined with surface free energy and molecular docking, the results revealed that polyphenols are driven by multiple forces that affect the orientation freedom of PPIs, with hydrogen bonding, hydrophobic interactions, and van der Waals forces being the major contributors. Overall, our work provides valuable insights into how molecules tune PPIs to modulate their function.

Deciphering the free energy landscapes of SARS-CoV-2 wild type and Omicron variant interacting with human ACE2

The binding of the receptor binding domain (RBD) of the SARS-CoV-2 spike protein to the host cell receptor angiotensin-converting enzyme 2 (ACE2) is the first step in human viral infection. Therefore, understanding the mechanism of interaction between RBD and ACE2 at the molecular level is critical for the prevention of COVID-19, as more variants of concern, such as Omicron, appear. Recently, atomic force microscopy has been applied to characterize the free energy landscape of the RBD-ACE2 complex, including estimation of the distance between the transition state and the bound state, xu. Here, using a coarse-grained model and replica-exchange umbrella sampling, we studied the free energy landscape of both the wild type and Omicron subvariants BA.1 and XBB.1.5 interacting with ACE2. In agreement with experiment, we find that the wild type and Omicron subvariants have similar xu values, but Omicron binds ACE2 more strongly than the wild type, having a lower dissociation constant KD.

Single-molecule scale quantification reveals interactions underlying protein-protein interface: from forces to non-covalent bonds

Protein-protein interactions (PPIs) between the B-cell lymphoma 2 (Bcl-2) family are considered a major driving force in cell cycle regulation and signaling. However, how this interfacial noncovalent interaction is achieved molecularly remains poorly understood. Herein, anti-apoptotic protein (Bcl-2) and pro-apoptotic protein (BAX) were used as models and their PPIs were explored for the first time using atomic force microscopy-based single-molecule force spectroscopy (SMFS) and in silico approaches. In addition, we used advanced analytical models, including multiple kinetic models, thermodynamic models, Poisson distributions, and contact angle molecular recognition to fully reveal the complexity of the BAX/Bcl-2 interaction interfaces. We propose that the binding kinetics between BAX/Bcl-2 are mainly mediated by specific (hydrogen bonding) and non-specific forces (hydrophobic interactions and electrostatic interactions) and show that the complicated multivalent binding interaction induces stable BAX/Bcl-2 complexes. This study enriches our understanding of the molecular mechanisms by which BAX interacts with Bcl-2. It provides valuable insights into the physical factors that need to be considered when designing PPI inhibitors.

Culinary spices and herbs in managing early and long-COVID-19 complications: A comprehensive review

Human race is preparing for the upsurge and aftermath of COVID-19 pandemic complicated by novel variants, new waves, variable mortality rate, and post-COVID complications. Despite use of repurposed drugs, symptomatic treatments and licensing of multiple vaccines, the daily number of cases and rate of transmission are significant. Culinary spices and herbs have been historically used in pandemic and non-pandemic times to reduce respiratory viral burden. Specific food items and culinary spices can boost the levels of protective immunity and also offer therapeutic benefits against impervious bugs via well-known as well as less-known but scientifically testable mechanisms. Here, we analyzed the phytochemicals profile of Ayurvedic herbs and inferred from the clinical trials/observational studies to provide a focused and succinct perspective on the relevance of "food-based" traditional decoction to moderate COVID-19 disease and long-COVID via modulation of immunity and reinstatement of homeostasis. We also underscore the druggable targets in pathogenesis of COVID-19 which are relevant to the ongoing clinical trials using spices and herbs. This information will provide a strong scientific rationale for standardization of the traditional herbs-based therapies and adopting the use of herbs, spices, and their formulations for reducing SARS-CoV-2 transmission, long-COVID symptoms, and COVID-19 disease progression.

Riding the Wave: Unveiling the Conformational Waves from RBD of SARS-CoV-2 Spike Protein to ACE2

The binding affinity between angiotensin-converting enzyme 2 (ACE2) and the receptor-binding domain (RBD) plays a crucial role in the transmission and reinfection of SARS-CoV2. Here, microsecond molecular dynamics simulations revealed that point mutations in the RBD domain induced conformational transitions that determined the binding affinity between ACE2 and RBD. These structural changes propagated through the RBD domain, altering the orientation of both ACE2 and RBD residues at the binding site. ACE2 receptor shows significant structural heterogeneity, whereas its binding to the RBD domain indicates a much greater degree of structural homogeneity. The receptor was more flexible in its unbound state with the binding of RBD domains inducing structural transitions. The structural heterogeneity observed in the ACE2 unbound form plays a role in the promiscuity of viral entry, as it may allow the receptor to interact with various related and unrelated ligands. Furthermore, rigidity may be important for stabilizing the complex and ensuring the proper orientation of the RBD-binding interface with ACE2. The greater structural homogeneity observed in the ACE2-RBD complex revealed the effectiveness of neutralizing antibodies and vaccines that are primarily directed toward the RBD-binding interface. The binding of the B38 monoclonal antibody revealed restricted conformational transitions in the RBD and ACE2 receptors, attributed to its potent binding interaction.

Interaction of SARS-CoV-2 with host cells and antibodies: experiment and simulation

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the devastating global COVID-19 pandemic announced by WHO in March 2020. Through unprecedented scientific effort, several vaccines, drugs and antibodies have been developed, saving millions of lives, but the fight against COVID-19 continues as immune escape variants of concern such as Delta and Omicron emerge. To develop more effective treatments and to elucidate the side effects caused by vaccines and therapeutic agents, a deeper understanding of the molecular interactions of SARS-CoV-2 with them and human cells is required. With special interest in computational approaches, we will focus on the structure of SARS-CoV-2 and the interaction of its spike protein with human angiotensin-converting enzyme-2 (ACE2) as a prime entry point of the virus into host cells. In addition, other possible viral receptors will be considered. The fusion of viral and human membranes and the interaction of the spike protein with antibodies and nanobodies will be discussed, as well as the effect of SARS-CoV-2 on protein synthesis in host cells.

Identification of conserved linear epitopes in the SARS-CoV-2 receptor-binding region using monoclonal antibodies

The ongoing coronavirus disease 2019 (COVID-19) pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has caused millions of cases of infections, leading to a global health emergency. The SARS-CoV-2 spike (S) protein plays the most important role in viral infection, and S1 subunit and its receptor-binding domain (RBD) are widely considered the most attractive vaccine targets. The RBD is highly immunogenic and its linear epitopes are important for vaccine development and therapy, but linear epitopes on the RBD have rarely been reported. In this study, 151 mouse monoclonal antibodies (mAbs) against the SARS-CoV-2 S1 protein were characterized and used to identify epitopes. Fifty-one mAbs reacted with eukaryotic SARS-CoV-2 RBD. Sixty-nine mAbs reacted with the S proteins of Omicron variants B.1.1.529 and BA.5, indicating their potential as rapid diagnostic materials. Three novel linear epitopes of RBD, R6 (³⁹¹CFTNVYADSFVIRGD⁴⁰⁵), R12 (⁴⁶³PFERDISTEIYQAGS⁴⁷⁷), and R16 (⁵¹⁰VVVLSFELLHAPAT⁵²³), were identified; these were highly conserved in SARS-CoV-2 variants of concern and could be detected in the convalescent serum of COVID-19 patients. From pseudovirus neutralization assays, some mAbs including one detecting R12 were found to possess neutralizing activity. Together, from the reaction of mAbs with eukaryotic RBD (N501Y), RBD (E484K), and S1 (D614G), we found that a single amino acid mutation in the SARS-CoV-2 S protein may cause a structural alteration, exerting substantial impact on mAb recognition. Our results could, therefore, help us better understand the function of the SARS-CoV-2 S protein and develop diagnostic tools for COVID-19.

Molecular docking and dynamics simulation approach of Camellia sinensis leaf extract derived compounds as potential cholinesterase inhibitors

The tea plant (Camellia sinensis) belongs to the family Theaceae and contains many phytochemicals that are effective against various diseases, including neurodegenerative disorders. In this study, we aimed to characterize the phytochemicals present in the methanolic and n-hexane leaf extracts of C. sinensis using GC-MS, FTIR, and UV-visible analysis. We detected a total of 19 compounds of different chemical classes. We also performed molecular docking studies using the GC-MS detected phytochemicals, targeting acetylcholinesterase (AChE, PBD ID: 4BDT) and butyrylcholinesterase (BChE, PDB ID: 6QAB), which are responsible for the breakdown of the neurotransmitter acetylcholine (ACh). This breakdown leads to dementia and cognitive decline in Alzheimer's patients. The compounds Ergosta-7,22-dien-3-ol, (3.beta.,5.alpha.,22E)- and Benzene, 1,3-bis(1,1-dimethylethyl) showed better binding affinity against AChE, while dl-.alpha.-Tocopherol and Ergosta-7,22-dien-3-ol, (3.beta.,5.alpha.,22E)- showed better binding affinity against BChE. We determined the stability and rigidity of these best docked complexes through molecular dynamics simulation for a period of 100 ns. All complexes showed stability in terms of SASA, Rg, and hydrogen bonds, but some variations were found in the RMSD values. Our ADMET analysis revealed that all lead compounds are non-toxic. Therefore, these compounds could be potential inhibitors of AChE and BChE.

Graphical abstract

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

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

Single-Molecule Analysis of SARS-CoV-2 Binding to C-Type Lectin Receptors

Despite intense scrutiny throughout the pandemic, development of efficacious drugs against SARS-CoV-2 spread remains hindered. Understanding the underlying mechanisms of viral infection is fundamental for developing novel treatments. While angiotensin converting enzyme 2 (ACE2) is accepted as the key entry receptor of the virus, other infection mechanisms exist. Dendritic cell-specific intercellular adhesion molecule-3 grabbing non-integrin (DC-SIGN) and its counterpart DC-SIGN-related (DC-SIGNR, also known as L-SIGN) have been recognized as possessing functional roles in COVID-19 disease and binding to SARS-CoV-2 has been demonstrated previously with ensemble and qualitative techniques. Here we examine the thermodynamic and kinetic parameters of the ligand-receptor interaction between these C-type lectins and the SARS-CoV-2 S1 protein using force-distance curve-based AFM and biolayer interferometry. We evidence that the S1 receptor binding domain is likely involved in this bond formation. Further, we employed deglycosidases and examined a nonglycosylated S1 variant to confirm the significance of glycosylation in this interaction. We demonstrate that the high affinity interactions observed occur through a mechanism distinct from that of ACE2.

Novel perspective for protein-drug interaction analysis: atomic force microscope

Proteins are major drug targets, and drug-target interaction identification and analysis are important factors for drug discovery. Atomic force microscopy (AFM) is a powerful tool making it possible to image proteins with nanometric resolution and probe intermolecular forces under physiological conditions. We review recent studies conducted in the field of target protein drug discovery using AFM-based analysis technology, including drug-driven changes in nanomechanical properties of protein morphology and interactions. Underlying mechanisms (including thermodynamic and kinetic parameters) of the drug-target interaction and drug-modulating protein-protein interaction (PPI) on the surfaces of models or living cells are discussed. Furthermore, challenges and the outlook for the field are likewise discussed. Overall, this insight into the mechanical properties of protein-drug interactions provides an unprecedented information framework for rational drug discovery in the pharmaceutical field.

Physical-Chemical Regulation of Membrane Receptors Dynamics in Viral Invasion and Immune Defense

Mechanical cues dynamically regulate membrane receptors functions to trigger various physiological and pathological processes from viral invasion to immune defense. These cues mainly include various types of dynamic mechanical forces and the spatial confinement of plasma membrane. However, the molecular mechanisms of how they couple with biochemical cues in regulating membrane receptors functions still remain mysterious. Here, we review recent advances in methodologies of single-molecule biomechanical techniques and in novel biomechanical regulatory mechanisms of critical ligand recognition of viral and immune receptors including SARS-CoV-2 spike protein, T cell receptor (TCR) and other co-stimulatory immune receptors. Furthermore, we provide our perspectives of the general principle of how force-dependent kinetics determine the dynamic functions of membrane receptors and of biomechanical-mechanism-driven SARS-CoV-2 neutralizing antibody design and TCR engineering for T-cell-based therapies.

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.

Pathogen-host adhesion between SARS-CoV-2 spike proteins from different variants and human ACE2 studied at single-molecule and single-cell levels

The binding of the receptor binding domain (RBD) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein onto human angiotensin-converting enzyme 2 (ACE2) is considered as the first step for the virus to adhere onto the host cells during the infection. Here, we investigated the adhesion of spike proteins from different variants and ACE2 using single-molecule and single-cell force spectroscopy. We found that the unbinding force and binding probability of the spike protein from Delta variant to the ACE2 were the highest among the variants tested in our study at both single-molecule and single-cell levels. As the most popular variants, the Omicron variants have slightly higher unbinding force to the ACE2 than wild type. Molecular dynamics simulation showed that ACE2-RBD (Omicron BA.1) complex is destabilized by the E484A and Y505H mutations and stabilized by S477N and N501Y mutations, when compared with Delta variant. In addition, a neutralizing antibody, produced by immunization with wild type spike protein, could effectively inhibit the binding of spike proteins from wild type, Delta and Omicron variants (BA.1 and BA.5) onto ACE2. Our results provide new insight for the molecular mechanism of the adhesive interactions between spike protein and ACE2 and suggest that effective monoclonal antibody can be prepared using wild type spike protein against different variants.

Spike protein of SARS-CoV-2 Omicron variant: An in-silico study evaluating spike interactions and immune evasion

Background

The fundamentals of the infectivity and immune evasion of the SARS-CoV-2 Omicron variant are not yet fully understood. Here, we carried out an in-silico study analyzing the spike protein, the protein electrostatic potential, and the potential immune evasion.

Methods

The analysis was based on the structure of the spike protein from two SARS-CoV-2 variants, the original Wuhan and the Botswana (Omicron). The full-length genome sequences and protein sequences were obtained from databanks. The interaction of the spike proteins with the human Angiotensin Converting Enzyme 2 (ACE2) receptor was evaluated through the open-source software. The Immune Epitope Database was used to analyze the potential immune evasion of the viruses.

Results

Our data show that the Omicron spike protein resulted in 37 amino acid changes. The physicochemical properties of the spike had changed, and the electrostatic potentials differed between both variants. This resulted in a decrease in protein interactions, which does not establish a greater interaction with the ACE2 receptor. These changes compromise key receptor-binding motif residues in the SARS-CoV-2 spike protein that interact with neutralizing antibodies and ACE2.

Conclusions

These mutations appear to confer enhanced properties of infectivity. The Omicron variant appears to be more effective at evading immune responses.

Electrostatic Interactions Are the Primary Determinant of the Binding Affinity of SARS-CoV-2 Spike RBD to ACE2: A Computational Case Study of Omicron Variants

To explore the mechanistic origin that determines the binding affinity of SARS-CoV-2 spike receptor binding domain (RBD) to human angiotensin converting enzyme 2 (ACE2), we constructed the homology models of RBD-ACE2 complexes of four Omicron subvariants (BA.1, BA.2, BA.3 and BA.4/5), and compared them with wild type complex (RBD_WT-ACE2) in terms of various structural dynamic properties by molecular dynamics (MD) simulations and binding free energy (BFE) calculations. The results of MD simulations suggest that the RBDs of all the Omicron subvariants (RBD_OMIs) feature increased global structural fluctuations when compared with RBD_WT. Detailed comparison of BFE components reveals that the enhanced electrostatic attractive interactions are the main determinant of the higher ACE2-binding affinity of RBD_OMIs than RBD_WT, while the weakened electrostatic attractive interactions determine RBD of BA.4/5 subvariant (RBD_BA.4/5) lowest ACE2-binding affinity among all Omicron subvariants. The per-residue BFE decompositions and the hydrogen bond (HB) networks analyses indicate that the enhanced electrostatic attractive interactions are mainly through gain/loss of the positively/negatively charged residues, and the formation or destruction of the interfacial HBs and salt bridges can also largely affect the ACE2-binding affinity of RBD. It is worth pointing out that since Q493R plays the most important positive contribution in enhancing binding affinity, the absence of this mutation in RBD_BA.4/5 results in a significantly weaker binding affinity to ACE2 than other Omicron subvariants. Our results provide insight into the role of electrostatic interactions in determining of the binding affinity of SARS-CoV-2 RBD to human ACE2.

Validation of Viral Inactivation Protocols for Therapeutic Blood Products against Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-COV-2)

Therapeutic blood products including convalescent plasma/serum and immunoglobulins concentrated from convalescent plasma, such as intravenous immunoglobulins or hyperimmune globulins, and monoclonal antibodies are passive immunotherapy options for novel coronavirus disease 2019 (COVID-19). They have been shown to improve the clinical status and biological and radiological parameters in some groups of COVID-19 patients. However, blood products are still potential sources of virus transmission in recipients. The use of pathogen reduction technology (PRT) should increase the safety of the products. The purpose of this study was to determine the impact of solvent/detergents (S/D) procedures on SARS-COV-2 infectivity elimination in the plasma of donors but also on COVID-19 convalescent serum (CCS) capacity to neutralize SARS-COV-2 infectivity. In this investigation, S/D treatment for all experiments was performed at a shortened process time (30 min). We first evaluated the impact of S/D treatments (1% TnBP/1% TritonX-45 and 1% TnBP/1% TritonX-100) on the inactivation of SARS-COV-2 pseudoparticles (SARS-COV-2pp)-spiked human plasma followed by S/D agent removal using a Sep-Pak Plus C18 cartridge. Both treatments were able to completely inactivate SARS-COV-2pp infectivity to an undetectable level. Moreover, the neutralizing activity of CCS against SARS-COV-2pp was preserved after S/D treatments. Our data suggested that viral inactivation methods using such S/D treatments could be useful in the implementation of viral inactivation/elimination processes of therapeutic blood products against SARS-COV-2.

When the Dust Has Settled: Calculation of Binding Affinities from First Principles for SARS-CoV-2 Variants with Quantitative Accuracy

Accurate determination of binding free energy is pivotal for the study of many biological processes and has been applied in a number of theoretical investigations to compare the affinity of severe acute respiratory syndrome coronavirus 2 variants toward the host cell. Diversity of these variants challenges the development of effective general therapies, their transmissibility relying either on an increased affinity toward their dedicated human receptor, the angiotensin-converting enzyme 2 (ACE2), or on escaping the immune response. Now that robust structural data are available, we have determined with utmost accuracy the standard binding free energy of the receptor-binding domain to the most widespread variants, namely, Alpha, Beta, Delta, and Omicron BA.2, as well as the wild type (WT) in complex either with ACE2 or with antibodies, namely, S2E12 and H11-D4, using a rigorous theoretical framework that combines molecular dynamics and potential-of-mean-force calculations. Our results show that an appropriate starting structure is crucial to ensure appropriate reproduction of the binding affinity, allowing the variants to be compared. They also emphasize the necessity to apply the relevant methodology, bereft of any shortcut, to account for all the contributions to the standard binding free energy. Our estimates of the binding affinities support the view that while the Alpha and Beta variants lean on an increased affinity toward the host cell, the Delta and Omicron BA.2 variants choose immune escape. Moreover, the S2E12 antibody, already known to be active against the WT (Starr et al., 2021; Mlcochova et al., 2021), proved to be equally effective against the Delta variant. In stark contrast, H11-D4 retains a low affinity toward the WT compared to that of ACE2 for the latter. Assuming robust structural information, the methodology employed herein successfully addresses the challenging protein-protein binding problem in the context of coronavirus disease 2019 while offering promising perspectives for predictive studies of ever-emerging variants.

Evolution of Stronger SARS-CoV-2 Variants as Revealed Through the Lens of Molecular Dynamics Simulations

Using molecular dynamics simulations, the protein-protein interactions of the receptor-binding domain of the wild-type and seven variants of the severe acute respiratory syndrome coronavirus 2 spike protein and the peptidase domain of human angiotensin-converting enzyme 2 were investigated. These variants are alpha, beta, gamma, delta, eta, kappa, and omicron. Using 100 ns simulation data, the residue interaction networks at the protein-protein interface were identified. Also, the impact of mutations on essential protein dynamics, backbone flexibility, and interaction energy of the simulated protein-protein complexes were studied. The protein-protein interface for the wild-type, delta, and omicron variants contained several stronger interactions, while the alpha, beta, gamma, eta, and kappa variants exhibited an opposite scenario as evident from the analysis of the inter-residue interaction distances and pair-wise interaction energies. The study reveals that two distinct residue networks at the central and right contact regions forge stronger binding affinity between the protein partners. The study provides a molecular-level insight into how enhanced transmissibility and infectivity by delta and omicron variants are most likely tied to a handful of interacting residues at the binding interface, which could potentially be utilized for future antibody constructs and structure-based antiviral drug design.

Mutation informatics: SARS-CoV-2 receptor-binding domain of the spike protein

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) undergoes mutations at a high rate and with frequent genetic reassortment (antigenic drift/shift), leading to variability in targets. The receptor-binding domain (RBD) of the spike (S) protein has a major role in the binding of SARS-CoV-2 with human angiotensin-converting enzyme 2 (ACE2). Mutations at the RBD influence the binding interaction at the SARS-CoV-2 S-ACE2 interface and impact viral pathogenicity. Here, we discuss different reported mutations of concern in RBD, physicochemical characteristic changes resulting from mutated amino acids and their effect on binding between the RBD and ACE2. Along with mutation informatics, we highlight recently developed small-molecule inhibitors of RBD and the ACE2 interface. This information provides a rational basis for the design of inhibitors against the multivariant strains of SARS-CoV-2.

Mechanosensing view of SARS-CoV-2 infection by a DNA nano-assembly

The mechanical force between a virus and its host cell plays a critical role in viral infection. However, characterization of the virus-cell mechanical force at the whole-virus level remains a challenge. Herein, we develop a platform in which the virus is anchored with multivalence-controlled aptamers to achieve transfer of the virus-cell mechanical force to a DNA tension gauge tether (Virus-TGT). When the TGT is ruptured, the complex of binding module-virus-cell is detached from the substrate, accompanied by decreased host cell-substrate adhesion, thus revealing the mechanical force between whole-virus and cell. Using Virus-TGT, direct evidence about the biomechanical force between SARS-CoV-2 and the host cell is obtained. The relative mechanical force gap (<10 pN) at the cellular level between the wild-type virus to cell and a variant virus to cell is measured, suggesting a possible positive correlation between virus-cell mechanical force and infectivity. Overall, this strategy provides a new perspective to probe the SARS-CoV-2 mechanical force.

Enhancing the Immunogenicity of RBD Protein Variants through Amino Acid E484 Mutation in SARS-CoV-2

In the context of the COVID-19 pandemic, conducting antibody testing and vaccination is critical. In particular, the continued evolution of SARS-CoV-2 raises concerns about the effectiveness of vaccines currently in use and the activity of neutralizing antibodies. Here, we used the Escherichia coli expression system to obtain nine different SARS-CoV-2 RBD protein variants, including six single-point mutants, one double-point mutant, and two three-point mutants. Western blotting results show that nine mutants of the RBD protein had strong antigenic activity in vitro. The immunogenicity of all RBD proteins was detected in mice to screen for protein mutants with high immunogenicity. The results show that the mutants E484K, E484Q, K417T-E484K-N501Y, and K417N-E484K-N501Y, especially the former two, had better immunogenicity than the wild type. This suggests that site E484 has a significant impact on the function of the RBD protein. Our results demonstrate that recombinant RBD protein expressed in E. coli can be an effective tool for the development of antibody detection methods and vaccines.

Exceptional diversity and selection pressure on coronavirus host receptors in bats compared to other mammals

Pandemics originating from non-human animals highlight the need to understand how natural hosts have evolved in response to emerging human pathogens and which groups may be susceptible to infection and/or potential reservoirs to mitigate public health and conservation concerns. Multiple zoonotic coronaviruses, such as severe acute respiratory syndrome-associated coronavirus (SARS-CoV), SARS-CoV-2 and Middle Eastern respiratory syndrome-associated coronavirus (MERS-CoV), are hypothesized to have evolved in bats. We investigate angiotensin-converting enzyme 2 (ACE2), the host protein bound by SARS-CoV and SARS-CoV-2, and dipeptidyl-peptidase 4 (DPP4 or CD26), the host protein bound by MERS-CoV, in the largest bat datasets to date. Both the ACE2 and DPP4 genes are under strong selection pressure in bats, more so than in other mammals, and in residues that contact viruses. Additionally, mammalian groups vary in their similarity to humans in residues that contact SARS-CoV, SARS-CoV-2 and MERS-CoV, and increased similarity to humans in binding residues is broadly predictive of susceptibility to SARS-CoV-2. This work augments our understanding of the relationship between coronaviruses and mammals, particularly bats, provides taxonomically diverse data for studies of how host proteins are bound by coronaviruses and can inform surveillance, conservation and public health efforts.

Integrating Conformational Dynamics and Perturbation-Based Network Modeling for Mutational Profiling of Binding and Allostery in the SARS-CoV-2 Spike Variant Complexes with Antibodies: Balancing Local and Global Determinants of Mutational Escape Mechanisms

In this study, we combined all-atom MD simulations, the ensemble-based mutational scanning of protein stability and binding, and perturbation-based network profiling of allosteric interactions in the SARS-CoV-2 spike complexes with a panel of cross-reactive and ultra-potent single antibodies (B1-182.1 and A23-58.1) as well as antibody combinations (A19-61.1/B1-182.1 and A19-46.1/B1-182.1). Using this approach, we quantify the local and global effects of mutations in the complexes, identify protein stability centers, characterize binding energy hotspots, and predict the allosteric control points of long-range interactions and communications. Conformational dynamics and distance fluctuation analysis revealed the antibody-specific signatures of protein stability and flexibility of the spike complexes that can affect the pattern of mutational escape. A network-based perturbation approach for mutational profiling of allosteric residue potentials revealed how antibody binding can modulate allosteric interactions and identified allosteric control points that can form vulnerable sites for mutational escape. The results show that the protein stability and binding energetics of the SARS-CoV-2 spike complexes with the panel of ultrapotent antibodies are tolerant to the effect of Omicron mutations, which may be related to their neutralization efficiency. By employing an integrated analysis of conformational dynamics, binding energetics, and allosteric interactions, we found that the antibodies that neutralize the Omicron spike variant mediate the dominant binding energy hotpots in the conserved stability centers and allosteric control points in which mutations may be restricted by the requirements of the protein folding stability and binding to the host receptor. This study suggested a mechanism in which the patterns of escape mutants for the ultrapotent antibodies may not be solely determined by the binding interaction changes but are associated with the balance and tradeoffs of multiple local and global factors, including protein stability, binding affinity, and long-range interactions.

Neutralization Activity against SARS-CoV-2 Variants after Booster Vaccination in Populations without COVID-19: A Meta-Analysis

A number of SARS-CoV-2 variants that have evolved to have significant immune escape have emerged worldwide since the COVID-19 outbreak. The efficacy of prime vaccination is waning with the evolution of SARS-CoV-2, and the necessity of booster doses is more and more prominent. Therefore, this study aimed to compare the neutralization activity against the wild type and variants (Beta, Delta, and Omicron) in different prime-boost vaccination regimens. Electronic databases including PubMed, the Cochrane Library, Embase, medRxiv, Wanfang and CNKI were used to retrieve original studies. A total of 16 studies, 9 prime-boost vaccination regimes, and 3134 subjects were included in the meta-analysis and random effect models were used to estimate pooled neutralization titers. The neutralization activity against SARS-CoV-2 showed a significant decline with the evolution of the virus, especially in the populations primed with inactivated vaccines. For homologous immunization, only the populations boosted with mRNA vaccines consistently had a significant rise in neutralization titers (Beta: MD = 0.97; Delta: MD = 1.33; Omicron: MD = 0.74). While the heterologous immunization was more effective, the increment of neutralization titers against wild type, Beta, Delta and Omicron was 1.65 (95% CI: 1.32-1.96), 1.03 (95% CI: 0.53-1.54), 1.46 (95% CI: 1.07-1.85) and 1.15 (95% CI: 0.68-1.61), respectively. With the evolution of SARS-CoV-2, the effectiveness of prime immunization is waning. Although the administration of the booster dose could ameliorate the neutralization titers, homologous immunization regimens were gradually losing their effectiveness. Therefore, a heterologous booster dose is required, especially in populations primed with inactivated vaccines.

Kinetics of the SARS-CoV-2 Antibody Avidity Response Following Infection and Vaccination

Serological assays capable of measuring antibody responses induced by previous infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have been critical tools in the response to the COVID-19 pandemic. In this study, we use bead-based multiplex assays to measure IgG and IgA antibodies and IgG avidity to five SARS-CoV-2 antigens (Spike (S), receptor-binding domain (RBD), Nucleocapsid (N), S subunit 2, and Membrane-Envelope fusion (ME)). These assays were performed in several cohorts of healthcare workers and nursing home residents, who were followed for up to eleven months after SARS-CoV-2 infection or up to six months after vaccination. Our results show distinct kinetic patterns of antibody quantity (IgG and IgA) and avidity. While IgG and IgA antibody levels waned over time, with IgA antibody levels waning more rapidly, avidity increased with time after infection or vaccination. These contrasting kinetic patterns allow for the estimation of time since previous SARS-CoV-2 infection. Including avidity measurements in addition to antibody levels in a classification algorithm for estimating time since infection led to a substantial improvement in accuracy, from 62% to 78%. The inclusion of antibody avidity in panels of serological assays can yield valuable information for improving serosurveillance during SARS-CoV-2 epidemics.

Rational design of AAVrh10-vectored ACE2 functional domain to broadly block the cell entry of SARS-CoV-2 variants

The frequently emerging SARS-CoV-2 variants have weakened the effectiveness of existing COVID-19 vaccines and neutralizing antibody therapy. Nevertheless, the infections of SARS-CoV-2 variants still depend on angiotensin-converting enzyme 2 (ACE2) receptor-mediated cell entry, and thus the soluble human ACE2 (shACE2) is a potential decoy for broadly blocking SARS-CoV-2 variants. In this study, we firstly generated the recombinant AAVrh10-vectored shACE2 constructs, a kind of adeno-associated virus (AAV) serotype with pulmonary tissue tropism, and then validated its inhibition capacity against SARS-CoV-2 infection. To further optimize the minimized ACE2 functional domain candidates, a comprehensive analysis was performed to clarify the interactions between the ACE2 orthologs from various species and the receptor binding domain (RBD) of SARS-CoV-2 spike (S) protein. Based on the key interface amino acids, we designed a series of truncated ACE2 orthologs, and then assessed their potential affinity to bind to SARS-CoV-2 variants RBD in silico. Of note, we found that the 24-83aa fragment of dog ACE2 (dACE2_24-83) had a higher affinity to the RBD of SARS-CoV-2 variants than that of human ACE2. Importantly, AAVrh10-vectored shACE2 or dACE2_24-83 constructs exhibited a broadly blockage breadth against SARS-CoV-2 prototype and variants in vitro and ex vivo. Collectively, these data highlighted a promising therapeutic strategy against SARS-CoV-2 variants.

Point mutations in SARS-CoV-2 variants induce long-range dynamical perturbations in neutralizing antibodies

Monoclonal antibodies are emerging as a viable treatment for the coronavirus disease 19 (COVID-19). However, newly evolved variants of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) can reduce the efficacy of currently available antibodies and can diminish vaccine-induced immunity. Here, we demonstrate that the microscopic dynamics of neutralizing monoclonal antibodies can be profoundly modified by the mutations present in the spike proteins of the SARS-COV-2 variants currently circulating in the world population. The dynamical perturbations within the antibody structure, which alter the thermodynamics of antigen recognition, are diverse and can depend both on the nature of the antibody and on the spatial location of the spike mutation. The correlation between the motion of the antibody and that of the spike receptor binding domain (RBD) can also be changed, modulating binding affinity. Using protein-graph-connectivity networks, we delineated the mutant-induced modifications in the information-flow along allosteric pathway throughout the antibody. Changes in the collective dynamics were spatially distributed both locally and across long-range distances within the antibody. On the receptor side, we identified an anchor-like structural element that prevents the detachment of the antibodies; individual mutations there can significantly affect the antibody binding propensity. Our study provides insight into how virus neutralization by monoclonal antibodies can be impacted by local mutations in the epitope via a change in dynamics. This realization adds a new layer of sophistication to the efforts for rational design of monoclonal antibodies against new variants of SARS-CoV2, taking the allostery in the antibody into consideration.

Mutations in the new variants of SARS-CoV-2 spike protein modulates the dynamics of the neutralizing antibodies. Capturing such modulations from MD simulations and graph network model identifies the role of mutations in facilitating immune evasion.

Developing Pseudovirus-Based Neutralization Assay against Omicron-Included SARS-CoV-2 Variants

The global spread of SARS-CoV-2 and its variants poses a serious threat to human health worldwide. Recently, the emergence of Omicron has presented a new challenge to the prevention and control of the COVID-19 pandemic. A convenient and reliable in vitro neutralization assay is an important method for validating the efficiency of antibodies, vaccines, and other potential drugs. Here, we established an effective assay based on a pseudovirus carrying a full-length spike (S) protein of SARS-CoV-2 variants in the HIV-1 backbone, with a luciferase reporter gene inserted into the non-replicate pseudovirus genome. The key parameters for packaging the pseudovirus were optimized, including the ratio of the S protein expression plasmids to the HIV backbone plasmids and the collection time for the Alpha, Beta, Gamma, Kappa, and Omicron pseudovirus particles. The pseudovirus neutralization assay was validated using several approved or developed monoclonal antibodies, underscoring that Omicron can escape some neutralizing antibodies, such as REGN10987 and REGN10933, while S309 and ADG-2 still function with reduced neutralization capability. The neutralizing capacity of convalescent plasma from COVID-19 convalescent patients in Wuhan was tested against these pseudoviruses, revealing the immune evasion of Omicron. Our work established a practical pseudovirus-based neutralization assay for SARS-CoV-2 variants, which can be conducted safely under biosafety level-2 (BSL-2) conditions, and this assay will be a promising tool for studying and characterizing vaccines and therapeutic candidates against Omicron-included SARS-CoV-2 variants.

Newly Emerged Antiviral Strategies for SARS-CoV-2: From Deciphering Viral Protein Structural Function to the Development of Vaccines, Antibodies, and Small Molecules

Coronavirus disease 2019 (COVID-19) caused by the infection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has become the most severe health crisis, causing extraordinary economic disruption worldwide. SARS-CoV-2 is a single-stranded RNA-enveloped virus. The process of viral replication and particle packaging is finished in host cells. Viral proteins, including both structural and nonstructural proteins, play important roles in the viral life cycle, which also provides the targets of treatment. Therefore, a better understanding of the structural function of virus proteins is crucial to speed up the development of vaccines and therapeutic strategies. Currently, the structure and function of proteins encoded by the SARS-CoV-2 genome are reviewed by several studies. However, most of them are based on the analysis of SARS-CoV-1 particles, lacking a systematic review update for SARS-CoV-2. Here, we specifically focus on the structure and function of proteins encoded by SARS-CoV-2. Viral proteins that contribute to COVID-19 infection and disease pathogenesis are reviewed according to the most recent research findings. The structure-function correlation of viral proteins provides a fundamental rationale for vaccine development and targeted therapy. Then, current antiviral vaccines are updated, such as inactive viral vaccines and protein-based vaccines and DNA, mRNA, and circular RNA vaccines. A summary of other therapeutic options is also reviewed, including monoclonal antibodies such as a cross-neutralizer antibody, a constructed cobinding antibody, a dual functional monoclonal antibody, an antibody cocktail, and an engineered bispecific antibody, as well as peptide-based inhibitors, chemical compounds, and clustered regularly interspaced short palindromic repeats (CRISPR) exploration. Overall, viral proteins and their functions provide the basis for targeted therapy and vaccine development.

BioExcel Building Blocks Workflows (BioBB-Wfs), an integrated web-based platform for biomolecular simulations

We present BioExcel Building Blocks Workflows, a web-based graphical user interface (GUI) offering access to a collection of transversal pre-configured biomolecular simulation workflows assembled with the BioExcel Building Blocks library. Available workflows include Molecular Dynamics setup, protein-ligand docking, trajectory analyses and small molecule parameterization. Workflows can be launched in the platform or downloaded to be run in the users’ own premises. Remote launching of long executions to user's available High-Performance computers is possible, only requiring configuration of the appropriate access credentials. The web-based graphical user interface offers a high level of interactivity, with integration with the NGL viewer to visualize and check 3D structures, MDsrv to visualize trajectories, and Plotly to explore 2D plots. The server requires no login but is recommended to store the users’ projects and manage sensitive information such as remote credentials. Private projects can be made public and shared with colleagues with a simple URL. The tool will help biomolecular simulation users with the most common and repetitive processes by means of a very intuitive and interactive graphical user interface. The server is accessible at https://mmb.irbbarcelona.org/biobb-wfs.

Autophagy and evasion of the immune system by SARS-CoV-2. Structural features of the non-structural protein 6 from wild type and Omicron viral strains interacting with a model lipid bilayer

The viral cycle of SARS-CoV-2 is based on a complex interplay with the cellular machinery, which is mediated by specific proteins eluding or hijacking the cellular defense mechanisms. Among the complex pathways induced by the viral infection, autophagy is particularly crucial and is strongly influenced by the action of the non-structural protein 6 (Nsp6) interacting with the endoplasmic reticulum membrane. Importantly, differently from other non-structural proteins, Nsp6 is mutated in the recently emerged Omicron variant, suggesting a possible different role of autophagy. In this contribution we explore, for the first time, the structural properties of Nsp6 thanks to long-timescale molecular dynamics simulations and machine learning analysis, identifying the interaction patterns with the lipid membrane. We also show how the mutation brought by the Omicron variant may indeed modify some of the specific interactions, and more particularly help anchor the viral protein to the lipid bilayer interface.

The viral cycle of SARS-CoV-2 is based on a complex interplay with the cellular machinery, which is mediated by specific proteins eluding or hijacking the cellular defense mechanisms.

Binding Interactions between Receptor-Binding Domain of Spike Protein and Human Angiotensin Converting Enzyme-2 in Omicron Variant

The emergence of new SARS-CoV-2 Omicron variant of concern (OV) has exacerbated the COVID-19 pandemic because of a large number of mutations in the spike protein, particularly in the receptor-binding domain (RBD), resulting in highly contagious and/or vaccine-resistant strains. Herein, we present a systematic analysis based on detailed molecular dynamics (MD) simulations in order to understand how the OV RBD mutations affect the ACE2 binding. We show that the OV RBD binds to ACE2 more efficiently and tightly predominantly because of strong electrostatic interactions, thereby promoting increased infectivity and transmissibility compared to other strains. Some of the OV RBD mutations are predicted to affect the antibody neutralization either through their role in the S-protein conformational changes, such as S371L, S373P, and S375F, or through changing its surface charge distribution, such as G339D, N440K, T478K, and E484A. Other mutations, such as K417N, G446S, and Y505H, decrease the ACE2 binding, whereas S447N, Q493R, G496S, Q498R, and N501Y tend to increase it.

Atomic force microscopy applied to interrogate nanoscale cellular chemistry and supramolecular bond dynamics for biomedical applications

Understanding biological interactions at a molecular level grants valuable information relevant to improving medical treatments and outcomes. Among the suite of technologies available, Atomic Force Microscopy (AFM) is unique in its ability to quantitatively probe forces and receptor-ligand interactions in real-time. The ability to assess the formation of supramolecular bonds and intermediates in real-time on surfaces and living cells generates important information relevant to understanding biological phenomena. Combining AFM with fluorescence-based techniques allows for an unprecedented level of insight not only concerning the formation and rupture of bonds, but understanding medically relevant interactions at a molecular level. As the ability of AFM to probe cells and more complex models improves, being able to assess binding kinetics, chemical topographies, and garner spectroscopic information will likely become key to developing further improvements in fields such as cancer, nanomaterials, and virology. The rapid response to the COVID-19 crisis, producing information regarding not just receptor affinities, but also strain-dependent efficacy of neutralizing nanobodies, demonstrates just how viable and integral to the pre-clinical development of information AFM techniques are in this era of medicine.

Mutational Effect of Some Major COVID-19 Variants on Binding of the S Protein to ACE2

COVID-19 is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which has many variants that accelerated the spread of the virus. In this study, we investigated the quantitative effect of some major mutants of the spike protein of SARS-CoV-2 binding to the human angiotensin-converting enzyme 2 (ACE2). These mutations are directly related to the Variant of Concern (VOC) including Alpha, Beta, Gamma, Delta and Omicron. Our calculations show that five major mutations (N501Y, E484K, L452R, T478K and K417N), first reported in Alpha, Beta, Gamma and Delta variants, all increase the binding of the S protein to ACE2 (except K417N), consistent with the experimental findings. We also studied an additional eight mutations of the Omicron variant that are located on the interface of the receptor binding domain (RDB) and have not been reported in other VOCs. Our study showed that most of these mutations (except Y505H and G446S) enhance the binding of the S protein to ACE2. The computational predictions helped explain why the Omicron variant quickly became dominant worldwide. Finally, comparison of several different computational methods for binding free energy calculation of these mutants was made. The alanine scanning method used in the current calculation helped to elucidate the residue-specific interactions responsible for the enhanced binding affinities of the mutants. The results show that the ASGB (alanine scanning with generalized Born) method is an efficient and reliable method for these binding free energy calculations due to mutations.

Functional Analysis of Spike from SARS-CoV-2 Variants Reveals the Role of Distinct Mutations in Neutralization Potential and Viral Infectivity

Enhanced viral transmission and escape from vaccine–elicited neutralizing antibodies drive worldwide spread of SARS-CoV-2 variants and promote disease progression. However, the impact of specific spike mutations that are carried by different viral variants on viral infectivity and neutralization sensitivity has not been completely defined. Here, we use pseudoviruses to assess the contribution of spike mutations within the Receptor Binding Domain (RBD) and the Furin Cleavage Site (FCS), and appear in circulating viral variants, on viral infectivity and neutralization potential against sera that was drawn from fully vaccinated individuals. Our functional analysis demonstrates that single, P681H, P681R or A701V–FCS mutations do not play a role in viral infectivity and neutralization potential. However, when in conjunction with the RBD–N501Y mutation, viral infectivity is enhanced. Similarly, combining the E484K–RBD mutation to the spike that carries FCS mutations reduces neutralization sensitivity with no effects on viral infectivity. Employing a similar approach onto the spike from Delta or Lota SARS-CoV-2 variants further reveals that specific RBD mutations affect neutralization sensitivity or viral infectivity differently. Our results validate the efficacy of the Pfizer third dose vaccine against Delta and Lota SARS-CoV-2 variants, and outline the significance of distinct RBD mutations in promoting viral infectivity and neutralization sensitivity to post–vaccination sera.

Molecular Dynamics and MM-PBSA Analysis of the SARS-CoV-2 Gamma Variant in Complex with the hACE-2 Receptor

The SARS-CoV-2 virus, since its appearance in 2019, has caused millions of cases and deaths. To date, there is no effective treatment or a vaccine that is fully protective. Despite the efforts made by governments and health institutions around the globe to control its propagation, the evolution of the virus has accelerated, diverging into hundreds of variants. However, not all of them are variants of concern (VoC’s). VoC’s have appeared in different regions and throughout the two years of the pandemic they have spread around the world. Specifically, in South America, the gamma variant (previously known as P.1) appeared in early 2021, bringing with it a second wave of infections. This variant contains the N501Y, E484K and K417T mutations in the receptor binding domain (RBD) of the spike protein. Although these mutations have been described experimentally, there is still no clarity regarding their role in the stabilization of the complex with the human angiotensin converting enzyme 2 (hACE-2) receptor. In this article we dissect the influence of mutations on the interaction with the hACE-2 receptor using molecular dynamics and estimations of binding affinity through a screened version of the molecular mechanics Poisson Boltzmann surface area (MM-PBSA) and interaction entropy. Our results indicate that mutations E484K and K417T compensate each other in terms of binding affinity, while the mutation N501Y promotes a more convoluted effect. This effect consists in the adoption of a cis configuration in the backbone of residue Y495 within the RBD, which in turn promotes polar interactions with the hACE-2 receptor. These results not only correlate with experimental observations and complement previous knowledge, but also expose new features associated with the specific contribution of concerned mutations. Additionally, we propose a recipe to assess the residue-specific contribution to the interaction entropy.

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

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

Martini 3 Model of Cellulose Microfibrils: On the Route to Capture Large Conformational Changes of Polysaccharides

High resolution data from all-atom molecular simulations is used to parameterize a Martini 3 coarse-grained (CG) model of cellulose I allomorphs and cellulose type-II fibrils. In this case, elementary molecules are represented by four effective beads centred in the positions of O2, O3, C6, and O6 atoms in the D-glucose cellulose subunit. Non-bonded interactions between CG beads are tuned according to a low statistical criterion of structural deviation using the Martini 3 type of interactions and are capable of being indistinguishable for all studied cases. To maintain the crystalline structure of each single cellulose chain in the microfibrils, elastic potentials are employed to retain the ribbon-like structure in each chain. We find that our model is capable of describing different fibril-twist angles associated with each type of cellulose fibril in close agreement with atomistic simulation. Furthermore, our CG model poses a very small deviation from the native-like structure, making it appropriate to capture large conformational changes such as those that occur during the self-assembly process. We expect to provide a computational model suitable for several new applications such as cellulose self-assembly in different aqueous solutions and the thermal treatment of fibrils of great importance in bioindustrial applications.