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

Co-Mutations and Possible Variation Tendency of the Spike RBD and Membrane Protein in SARS-CoV-2 by Machine Learning

Since the onset of the coronavirus disease 2019 (COVID-19) pandemic, SARS-CoV-2 variants capable of breakthrough infections have attracted global attention. These variants have significant mutations in the receptor-binding domain (RBD) of the spike protein and the membrane (M) protein, which may imply an enhanced ability to evade immune responses. In this study, an examination of co-mutations within the spike RBD and their potential correlation with mutations in the M protein was conducted. The EVmutation method was utilized to analyze the distribution of the mutations to elucidate the relationship between the mutations in the spike RBD and the alterations in the M protein. Additionally, the Sequence-to-Sequence Transformer Model (S2STM) was employed to establish mapping between the amino acid sequences of the spike RBD and M proteins, offering a novel and efficient approach for streamlined sequence analysis and the exploration of their interrelationship. Certain mutations in the spike RBD, G339D-S373P-S375F and Q493R-Q498R-Y505, are associated with a heightened propensity for inducing mutations at specific sites within the M protein, especially sites 3 and 19/63. These results shed light on the concept of mutational synergy between the spike RBD and M proteins, illuminating a potential mechanism that could be driving the evolution of SARS-CoV-2.

Molecular Imaging of ACE2 Expression in Infectious Disease and Cancer

Angiotensin-converting enzyme 2 (ACE2) is a cell-surface receptor that plays a critical role in the pathogenesis of SARS-CoV-2 infection. Through the use of ligands engineered for the receptor, ACE2 imaging has emerged as a valuable tool for preclinical and clinical research. These can be used to visualize the expression and distribution of ACE2 in tissues and cells. A variety of techniques including optical, magnetic resonance, and nuclear medicine contrast agents have been developed and employed in the preclinical setting. Positron-emitting radiotracers for highly sensitive and quantitative tomography have also been translated in the context of SARS-CoV-2-infected and control patients. Together this information can be used to better understand the mechanisms of SARS-CoV-2 infection, the potential roles of ACE2 in homeostasis and disease, and to identify potential therapeutic modulators in infectious disease and cancer. This review summarizes the tools and techniques to detect and delineate ACE2 in this rapidly expanding field.

Promising Schiff bases in antiviral drug design and discovery

Emerging and re-emerging illnesses will probably present a new hazard of infectious diseases and have fostered the urge to research new antiviral agents. Most of the antiviral agents are analogs of nucleosides and only a few are non-nucleoside antiviral agents. There is quite a less percentage of marketed/clinically approved non-nucleoside antiviral medications. Schiff bases are organic compounds that possess a well-demonstrated profile against cancer, viruses, fungus, and bacteria, as well as in the management of diabetes, chemotherapy-resistant cases, and malarial infections. Schiff bases resemble aldehydes or ketones with an imine/azomethine group instead of a carbonyl ring. Schiff bases have a broad application profile not only in therapeutics/medicine but also in industrial applications. Researchers have synthesized and screened various Schiff base analogs for their antiviral potential. Some of the important heterocyclic compounds like istatin, thiosemicarbazide, quinazoline, quinoyl acetohydrazide, etc. have been used to derive novel Schiff base analogs. Keeping in view the outbreak of viral pandemics and epidemics, this manuscript compiles a review of Schiff base analogs concerning their antiviral properties and structural-activity relationship analysis.

Extracellular Vesicle-Based SARS-CoV-2 Vaccine

Messenger ribonucleic acid (RNA) vaccines are mainly used as SARS-CoV-2 vaccines. Despite several issues concerning storage, stability, effective period, and side effects, viral vector vaccines are widely used for the prevention and treatment of various diseases. Recently, viral vector-encapsulated extracellular vesicles (EVs) have been suggested as useful tools, owing to their safety and ability to escape from neutral antibodies. Herein, we summarize the possible cellular mechanisms underlying EV-based SARS-CoV-2 vaccines.

Molecular Evolution of SARS-CoV-2 during the COVID-19 Pandemic

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) produced diverse molecular variants during its recent expansion in humans that caused different transmissibility and severity of the associated disease as well as resistance to monoclonal antibodies and polyclonal sera, among other treatments. In order to understand the causes and consequences of the observed SARS-CoV-2 molecular diversity, a variety of recent studies investigated the molecular evolution of this virus during its expansion in humans. In general, this virus evolves with a moderate rate of evolution, in the order of 10-3-10-4 substitutions per site and per year, which presents continuous fluctuations over time. Despite its origin being frequently associated with recombination events between related coronaviruses, little evidence of recombination was detected, and it was mostly located in the spike coding region. Molecular adaptation is heterogeneous among SARS-CoV-2 genes. Although most of the genes evolved under purifying selection, several genes showed genetic signatures of diversifying selection, including a number of positively selected sites that affect proteins relevant for the virus replication. Here, we review current knowledge about the molecular evolution of SARS-CoV-2 in humans, including the emergence and establishment of variants of concern. We also clarify relationships between the nomenclatures of SARS-CoV-2 lineages. We conclude that the molecular evolution of this virus should be monitored over time for predicting relevant phenotypic consequences and designing future efficient treatments.

High-Throughput Molecular Dynamics-Based Alchemical Free Energy Calculations for Predicting the Binding Free Energy Change Associated with the Selected Omicron Mutations in the Spike Receptor-Binding Domain of SARS-CoV-2

The ongoing pandemic caused by SARS-CoV-2 has gone through various phases. Since the initial outbreak, the virus has mutated several times, with some lineages showing even stronger infectivity and faster spread than the original virus. Among all the variants, omicron is currently classified as a variant of concern (VOC) by the World Health Organization, as the previously circulating variants have been replaced by it. In this work, we have focused on the mutations observed in omicron sub lineages BA.1, BA.2, BA.4 and BA.5, particularly at the receptor-binding domain (RBD) of the spike protein that is responsible for the interactions with the host ACE2 receptor and binding of antibodies. Studying such mutations is particularly important for understanding the viral infectivity, spread of the disease and for tracking the escape routes of this virus from antibodies. Molecular dynamics (MD) based alchemical free energy calculations have been shown to be very accurate in predicting the free energy change, due to a mutation that could have a deleterious or a stabilizing effect on either the protein itself or its binding affinity to another protein. Here, we investigated the significance of five spike RBD mutations on the stability of the spike protein binding to ACE2 by free energy calculations using high throughput MD simulations. For comparison, we also used conventional MD simulations combined with a Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) based approach, and compared our results with the available experimental data. Overall, the alchemical free energy calculations performed far better than the MM-GBSA approach in predicting the individual impact of the mutations. When considering the experimental variation, the alchemical free energy method was able to produce a relatively accurate prediction for N501Y, the mutant that has previously been reported to increase the binding affinity to hACE2. On the other hand, the other individual mutations seem not to have a significant effect on the spike RBD binding affinity towards hACE2.

Kinetics of Drug Molecule Interactions with a Newly Developed Nano-Gold-Modified Spike Protein Electrochemical Receptor Sensor

In March 2020, the World Health Organization (WHO) declared COVID-19 a pandemic, and the spike protein has been reported to be an important drug target for anti-COVID-19 treatment. As such, in this study, we successfully developed a novel electrochemical receptor biosensor by immobilizing the SARS-CoV-2 spike protein and using AuNPs-HRP as an electrochemical signal amplification system. Moreover, the time-current method was used to quantify seven antiviral drug compounds, such as arbidol and chloroquine diphosphate. The results show that the spike protein and the drugs are linearly correlated within a certain concentration range and that the detection sensitivity of the sensor is extremely high. In the low concentration range of linear response, the kinetics of receptor-ligand interactions are similar to that of an enzymatic reaction. Among the investigated drug molecules, bromhexine exhibits the smallest Ka value, and thus, is most sensitively detected by the sensor. Hydroxychloroquine exhibits the largest Ka value. Molecular docking simulations of the spike protein with six small-molecule drugs show that residues of this protein, such as Asp, Trp, Asn, and Gln, form hydrogen bonds with the -OH or -NH₂ groups on the branched chains of small-molecule drugs. The electrochemical receptor biosensor can directly quantify the interaction between the spike protein and drugs such as abidor and hydroxychloroquine and perform kinetic studies with a limit of detection 3.3 × 10^-20 mol/L, which provides a new research method and idea for receptor-ligand interactions and pharmacodynamic evaluation.

Cheminformatics-Based Discovery of Potential Chemical Probe Inhibitors of Omicron Spike Protein

During the past two decades, the world has witnessed the emergence of various SARS-CoV-2 variants with distinct mutational profiles influencing the global health, economy, and clinical aspects of the COVID-19 pandemic. These variants or mutants have raised major concerns regarding the protection provided by neutralizing monoclonal antibodies and vaccination, rates of virus transmission, and/or the risk of reinfection. The newly emerged Omicron, a genetically distinct lineage of SARS-CoV-2, continues its spread in the face of rising vaccine-induced immunity while maintaining its replication fitness. Efforts have been made to improve the therapeutic interventions and the FDA has issued Emergency Use Authorization for a few monoclonal antibodies and drug treatments for COVID-19. However, the current situation of rapidly spreading Omicron and its lineages demands the need for effective therapeutic interventions to reduce the COVID-19 pandemic. Several experimental studies have indicated that the FDA-approved monoclonal antibodies are less effective than antiviral drugs against the Omicron variant. Thus, in this study, we aim to identify antiviral compounds against the Spike protein of Omicron, which binds to the human angiotensin-converting enzyme 2 (ACE2) receptor and facilitates virus invasion. Initially, docking-based virtual screening of the in-house database was performed to extract the potential hit compounds against the Spike protein. The obtained hits were optimized by DFT calculations to determine the electronic properties and molecular reactivity of the compounds. Further, MD simulation studies were carried out to evaluate the dynamics of protein–ligand interactions at an atomistic level in a time-dependent manner. Collectively, five compounds (AKS-01, AKS-02, AKS-03, AKS-04, and AKS-05) with diverse scaffolds were identified as potential hits against the Spike protein of Omicron. Our study paves the way for further in vitro and in vivo studies.

Evaluation of humoral immune responses induced by different SARS-CoV-2 spike trimers from wild-type and emerging variants with individual, sequential, and combinational delivered strategies

The spike trimer of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an effective target for inducing neutralizing antibodies by coronavirus disease 2019 (COVID-19) vaccines. However, the diversity of spike protein from emerging SASR-CoV-2 variants has become the major challenge for development of a universal vaccine. To investigate the immunogenicity of spike proteins from various circulating strains including wild type, Delta, and Omicron variants, we produced various natural spike trimers and designed three vaccination strategies, that is, individual, sequential, and bivalent regimens to assess autologous and heterogenous antibody responses in a mouse model. The results indicated that monovalent vaccine strategy with individual spike trimer could only induce binding and neutralizing antibodies against homologous viruses. However, sequential and bivalent immunization with Delta and Omicron spike trimers could induce significantly broader neutralizing antibody responses against heterogenous SARS-CoV-2. Interestingly, the spike trimer from Omicron variant showed superior immunogenicity in inducing antibody response against recently emerging XE variant. Taken together, our data supported the development of novel vaccination strategies or multivalent vaccine against emerging variants.