Signatures in SARS-CoV-2 spike protein conferring escape to neutralizing antibodies

Understanding SARS-CoV-2 evolution and host immunity is critical to control COVID-19 pandemics. At the core is an arms-race between SARS-CoV-2 antibody and angiotensin-converting enzyme 2 (ACE2) recognition, a function of the viral protein spike. Mutations in spike impacting antibody and/or ACE2 binding are appearing worldwide, imposing the need to monitor SARS-CoV2 evolution and dynamics in the population. Determining signatures in SARS-CoV-2 that render the virus resistant to neutralizing antibodies is critical. We engineered 25 spike-pseudotyped lentiviruses containing individual and combined mutations in the spike protein, including all defining mutations in the variants of concern, to identify the effect of single and synergic amino acid substitutions in promoting immune escape. We confirmed that E484K evades antibody neutralization elicited by infection or vaccination, a capacity augmented when complemented by K417N and N501Y mutations. In silico analysis provided an explanation for E484K immune evasion. E484 frequently engages in interactions with antibodies but not with ACE2. Importantly, we identified a novel amino acid of concern, S494, which shares a similar pattern. Using the already circulating mutation S494P, we found that it reduces antibody neutralization of convalescent and post-immunization sera, particularly when combined with E484K and with mutations able to increase binding to ACE2, such as N501Y. Our analysis of synergic mutations provides a signature for hotspots for immune evasion and for targets of therapies, vaccines and diagnostics.     

In-silico genomic landscape characterization and evolution of SARS-CoV-2 variants isolated in India shows significant drift with high frequency of mutations

In-silico studies on SARS-CoV-2 genome are considered important to identify the significant pattern of variations and its possible effects on the structural and functional characteristics of the virus. The current study determined such genetic variations and their possible impact among SARS-CoV-2 variants isolated in India. A total of 546 SARS-CoV-2 genomic sequences (India) were retrieved from the gene bank (NCBI) and subjected to alignment against the Wuhan variant (NC_045512.2), the corresponding amino acid changes were analyzed using NCBI Protein-BLAST. These 546 variants revealed 841 mutations; most of these were non-synonymous 464/841 (55.1%), there was no identical variant compared to the original strain. All genes; coding and non-coding showed nucleotide changes, most of the structural genes showed frequent nonsynonymous mutations. The most affected genes were ORF1a/b followed by the S gene which showed 515/841 (61.2%) and 120/841 (14.3%) mutations, respectively. The most frequent non-synonymous mutation 486/546 (89.01%) occurred in the S gene (structural gene) at position 23,403 where A changed to G leading to the replacement of aspartic acid by glycine in position (D614G). Interestingly, four variants also showed deletion. The variants MT800923 and MT800925 showed 12 consecutive nucleotide deletion in position 21982–21993 resulting in 4 consecutive amino acid deletions that were leucine, glycine, valine, and tyrosine in positions 141, 142, 143, and 144 respectively. The present study exhibited a higher mutations rate per variant compared to other studies carried out in India.     

Analyzing the effect of mutations in SARS-CoV2 papain-like protease from Saudi isolates on protein structure and drug-protein binding: Molecular modelling and dynamics studies

The continuous and rapid development of the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) virus remains a health concern especially with the emergence of numerous variants and mutations worldwide. As with other RNA viruses, SARS-CoV-2 has a genetically high mutation rate. These mutations have an impact on the virus characteristics, including transmissibility, antigenicity and development of drug and vaccine resistance. This work was pursued to identify the differences that exist in the papain-like protease (PL^Pro) from 58 Saudi isolates in comparison to the first reported sequence from Wuhan, China and determine their implications on protein structure and the inhibitor binding. PL^pro is a key protease enzyme for the host cells invasion and viral proteolytic cleavage, hence, it emerges as a valuable antiviral therapeutic target. Two mutations were identified including D108G and A249V and shown to increase the molecular flexibility of PL^Pro protein and alter the protein stability, particularly with D108G mutation. The effect of these mutations on the stability and dynamic behavior of PL^Pro structures as well as their effect on the binding of a known inhibitor; GRL0617 were further investigated by molecular docking and dynamic simulation.     

Experimental Evidence for Enhanced Receptor Binding by Rapidly Spreading SARS-CoV-2 Variants

Rapidly spreading new variants of SARS-CoV-2 carry multiple mutations in the viral spike protein which attaches to the angiotensin converting enzyme 2 (ACE2) receptor on host cells. Among these mutations are amino acid changes N501Y (lineage B.1.1.7, first identified in the UK), and the combination N501Y, E484K, K417N (B.1.351, first identified in South Africa), all located at the interface on the receptor binding domain (RBD). We experimentally establish that RBD containing the N501Y mutation results in 7-fold stronger binding to the hACE2 receptor than wild type RBD. The E484K mutation only slightly enhances the affinity for the receptor, while K417N attenuates affinity. As a result, RBD from B.1.351 containing all three mutations binds 3-fold stronger to hACE2 than wild type RBD but 2-fold weaker than N501Y. However, the recently emerging double mutant E484K/N501Y binds even stronger than N501Y. The independent evolution of lineages containing mutations with different effects on receptor binding affinity, viral transmission and immune evasion underscores the importance of global viral genome surveillance and functional characterization.     

A Comparative Study between Spanish and British SARS-CoV-2 Variants

The study of the interaction between the SARS-CoV-2 spike protein and the angiotensin-converting enzyme 2 (ACE2) receptor is key to understanding binding affinity and stability. In the present report, we sought to investigate the differences between two already sequenced genome variants (Spanish and British) of SARS-CoV-2. Methods: In silico model evaluating the homology, identity and similarity in the genome sequence and the structure and alignment of the predictive spike by computational docking methods. Results: The identity results between the Spanish and British variants of the Spike protein were 28.67%. This close correspondence in the results between the Spanish and British SARS-CoV-2 variants shows that they are very similar (99.99%). The alignment obtained results in four deletions. There were 23 nucleotide substitutions also predicted which could affect the functionality of the proteins produced from this sequence. The interaction between the binding receptor domain from the spike protein and the ACE2 receptor produces some of the mutations found and, therefore, the energy of this ligand varies. However, the estimated antigenicity of the British variant is higher than its Spanish counterpart. Conclusions: Our results indicate that minimal mutations could interfere in the infectivity of the virus due to changes in the fitness between host cell recognition and interaction proteins. In particular, the N501Y substitution, situated in the RBD of the spike of the British variant, might be the reason for its extraordinary infective potential.     

Vaccine-escape and fast-growing mutations in the United Kingdom,the United States,Singapore, Spain,India, and other COVID-19-devastated countries

Recently, the SARS-CoV-2 variants from the United Kingdom (UK), South Africa, and Brazil have received much attention for their increased infectivity, potentially high virulence, and possible threats to existing vaccines and antibody therapies. The question remains if there are other more infectious variants transmitted around the world. We carry out a large-scale study of 506,768 SARS-CoV-2 genome isolates from patients to identify many other rapidly growing mutations on the spike (S) protein receptor-binding domain (RBD). We reveal that essentially all 100 most observed mutations strengthen the binding between the RBD and the host angiotensin-converting enzyme 2 (ACE2), indicating the virus evolves toward more infectious variants. In particular, we discover new fast-growing RBD mutations N439K, S477N, S477R, and N501T that also enhance the RBD and ACE2 binding. We further unveil that mutation N501Y involved in United Kingdom (UK), South Africa, and Brazil variants may moderately weaken the binding between the RBD and many known antibodies, while mutations E484K and K417N found in South Africa and Brazilian variants, L452R and E484Q found in India variants, can potentially disrupt the binding between the RBD and many known antibodies. Among these RBD mutations, L452R is also now known as part of the California variant B.1.427. Finally, we hypothesize that RBD mutations that can simultaneously make SARS-CoV-2 more infectious and disrupt the existing antibodies, called vaccine escape mutations, will pose an imminent threat to the current crop of vaccines. A list of most likely vaccine escape mutations is given, including S494P, Q493L, K417N, F490S, F486L, R403K, E484K, L452R, K417T, F490L, E484Q, and A475S. Mutation T478K appears to make the Mexico variant B.1.1.222 the most infectious one. Our comprehensive genetic analysis and protein-protein binding study show that the genetic evolution of SARS-CoV-2 on the RBD, which may be regulated by host gene editing, viral proofreading, random genetic drift, and natural selection, gives rise to more infectious variants that will potentially compromise existing vaccines and antibody therapies.     

Genotyping coronavirus SARS-CoV-2: methods and implications

The emerging global infectious COVID-19 disease by novel Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) presents critical threats to global public health and the economy since it was identified in late December 2019 in China. The virus has gone through various pathways of evolution. To understand the evolution and transmission of SARS-CoV-2, genotyping of virus isolates is of great importance. This study presents an accurate method for effectively genotyping SARS-CoV-2 viruses using complete genomes. The method employs the multiple sequence alignments of the genome isolates with the SARS-CoV-2 reference genome. The single-nucleotide polymorphism (SNP) genotypes are then measured by Jaccard distances to track the relationship of virus isolates. The genotyping analysis of SARS-CoV-2 isolates from the globe reveals that specific multiple mutations are the predominated mutation type during the current epidemic. The proposed method serves an effective tool for monitoring and tracking the epidemic of pathogenic viruses in their global and local genetic variations. The genotyping analysis shows that the genes encoding the S proteins and RNA polymerase, RNA primase, and nucleoprotein, undergo frequent mutations. These mutations are critical for vaccine development in disease control.     

Spatio-temporal dynamics of intra-host variability in SARS-CoV-2 genomes

During the course of the COVID-19 pandemic, large-scale genome sequencing of SARS-CoV-2 has been useful in tracking its spread and in identifying variants of concern (VOC). Viral and host factors could contribute to variability within a host that can be captured in next-generation sequencing reads as intra-host single nucleotide variations (iSNVs). Analysing 1347 samples collected till June 2020, we recorded 16 410 iSNV sites throughout the SARS-CoV-2 genome. We found ∼42% of the iSNV sites to be reported as SNVs by 30 September 2020 in consensus sequences submitted to GISAID, which increased to ∼80% by 30th June 2021. Following this, analysis of another set of 1774 samples sequenced in India between November 2020 and May 2021 revealed that majority of the Delta (B.1.617.2) and Kappa (B.1.617.1) lineage-defining variations appeared as iSNVs before getting fixed in the population. Besides, mutations in RdRp as well as RNA-editing by APOBEC and ADAR deaminases seem to contribute to the differential prevalence of iSNVs in hosts. We also observe hyper-variability at functionally critical residues in Spike protein that could alter the antigenicity and may contribute to immune escape. Thus, tracking and functional annotation of iSNVs in ongoing genome surveillance programs could be important for early identification of potential variants of concern and actionable interventions.     

Exhaustive mutational analysis of severe acute respiratory syndrome coronavirus 2 ORF3a: An essential component in the pathogen's infectivity cycle

Detailed knowledge of a protein's key residues may assist in understanding its function and designing inhibitors against it. Consequently, such knowledge of one of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)'s proteins is advantageous since the virus is the etiological agent behind one of the biggest health crises of recent times. To that end, we constructed an exhaustive library of bacteria differing from each other by the mutated version of the virus's ORF3a viroporin they harbor. Since the protein is harmful to bacterial growth due to its channel activity, genetic selection followed by deep sequencing could readily identify mutations that abolish the protein's function. Our results have yielded numerous mutations dispersed throughout the sequence that counteract ORF3a's ability to slow bacterial growth. Comparing these data with the conservation pattern of ORF3a within the coronavirinae provided interesting insights: Deleterious mutations obtained in our study corresponded to conserved residues in the protein. However, despite the comprehensive nature of our mutagenesis coverage (108 average mutations per site), we could not reveal all of the protein's conserved residues. Therefore, it is tempting to speculate that our study unearthed positions in the protein pertinent to channel activity, while other conserved residues may correspond to different functionalities of ORF3a. In conclusion, our study provides important information on a key component of SARS-CoV-2 and establishes a procedure to analyze other viroporins comprehensively.     

Detection of R.1 lineage severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with spike protein W152L/E484K/G769V mutations in Japan

We aimed to investigate novel emerging severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineages in Japan that harbor variants in the spike protein receptor-binding domain (RBD). The total nucleic acid contents of samples from 159 patients with coronavirus disease 2019 (COVID-19) were subjected to whole genome sequencing. The SARS-CoV-2 genome sequences from these patients were examined for variants in spike protein RBD. In January 2021, three family members (one aged in their 40s and two aged under 10 years old) were found to be infected with SARS-CoV-2 harboring W152L/E484K/G769V mutations. These three patients were living in Japan and had no history of traveling abroad. After identifying these cases, we developed a TaqMan assay to screen for the above hallmark mutations and identified an additional 14 patients with the same mutations. The associated virus strain was classified into the GR clade (Global Initiative on Sharing Avian Influenza Data [GISAID]), 20B clade (Nextstrain), and R.1 lineage (Phylogenetic Assignment of Named Global Outbreak [PANGO] Lineages). As of April 22, 2021, R.1 lineage SARS-CoV-2 has been identified in 2,388 SARS-CoV-2 entries in the GISAID database, many of which were from Japan (38.2%; 913/2,388) and the United States (47.1%; 1,125/2,388). Compared with that in the United States, the percentage of SARS-CoV-2 isolates belonging to the R.1 lineage in Japan increased more rapidly over the period from October 24, 2020 to April 18, 2021. R.1 lineage SARS-CoV-2 has potential escape mutations in the spike protein RBD (E484K) and N-terminal domain (W152L); therefore, it will be necessary to continue to monitor the R.1 lineage as it spreads around the world.     

Insights on cross-species transmission of SARS-CoV-2 from structural modeling

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the ongoing global pandemic that has infected more than 31 million people in more than 180 countries worldwide. Like other coronaviruses, SARS-CoV-2 is thought to have been transmitted to humans from wild animals. Given the scale and widespread geographical distribution of the current pandemic and confirmed cases of cross-species transmission, the question of the extent to which this transmission is possible emerges, as well as what molecular features distinguish susceptible from non-susceptible animal species. Here, we investigated the structural properties of several ACE2 orthologs bound to the SARS-CoV-2 spike protein. We found that species known not to be susceptible to SARS-CoV-2 infection have non-conservative mutations in several ACE2 amino acid residues that disrupt key polar and charged contacts with the viral spike protein. Our models also allow us to predict affinity-enhancing mutations that could be used to design ACE2 variants for therapeutic purposes. Finally, our study provides a blueprint for modeling viral-host protein interactions and highlights several important considerations when designing these computational studies and analyzing their results.     

MSH2 missense mutations alter cisplatin cytotoxicity and promote cisplatin-induced genome instability

Defects in the mismatch repair protein MSH2 cause tolerance to DNA damage. We report how cancer-derived and polymorphic MSH2 missense mutations affect cisplatin cytotoxicity. The chemotolerance phenotype was compared with the mutator phenotype in a yeast model system. MSH2 missense mutations display a strikingly different effect on cell death and genome instability. A mutator phenotype does not predict chemotolerance or vice versa. MSH2 mutations that were identified in tumors (Y109C) or as genetic variations (L402F) promote tolerance to cisplatin, but leave the initial mutation rate of cells unaltered. A secondary increase in the mutation rate is observed after cisplatin exposure in these strains. The mutation spectrum of cisplatin-resistant mutators identifies persistent cisplatin adduction as the cause for this acquired genome instability. Our results demonstrate that MSH2 missense mutations that were identified in tumors or as polymorphic variations can cause increased cisplatin tolerance independent of an initial mutator phenotype. Cisplatin exposure promotes drug-induced genome instability. From a mechanistical standpoint, these data demonstrate functional separation between MSH2-dependent cisplatin cytotoxicity and repair. From a clinical standpoint, these data provide valuable information on the consequences of point mutations for the success of chemotherapy and the risk for secondary carcinogenesis.     

Hydrolysis of oxidized nucleotides by the Escherichia coli Orf135 protein.

To examine the possibility that the Orf135 protein of Escherichia coli functions as a hydrolyzing enzyme for a damaged DNA precursor (deoxyribonucleoside 5'-triphosphate), we purified the recombinant Orf135 protein and incubated it with oxidized deoxynucleotides. Of the nucleotides tested, 2-hydroxydeoxyadenosine 5'-triphosphate, and somewhat less efficiently, 8-hydroxydeoxyguanosine 5'-triphosphate, were hydrolyzed by this protein. These damaged deoxynucleotides elicit transversion mutations in E. coli (Inoue, M., Kamiya, H., Fujikawa, K., Ootsuyama, Y., Murata-Kamiya, N., Osaki, T., Yasumoto, K., Kasai, H. (1998) J. Biol. Chem. 273, 11069-11074). These results suggest that this protein may be involved in the prevention of mutations induced by these oxidized deoxynucleotides.     

High-resolution melting analysis to discriminate between the SARS-CoV-2 Omicron variants BA.1 and BA.2

High-resolution melting (HRM) analysis was conducted to discriminate between SARS-CoV-2 Omicron variant BA.1 (B.1.1.529.1) and subvariant BA.2 (B.1.1.529.2). We performed two-step PCR consisting of the first PCR and the second nested PCR to prepare the amplicon for HRM analysis, which detected G339D, N440K, G446S and D796Y variations in the SARS-CoV-2 spike protein. The melting temperatures (Tms) of the amplicons from the cDNA of the Omicron variant BA.1 and subvariant BA.2 receptor binding domain (RBD) in spike protein were the same: 75.2 °C (G339D variation) and 73.4 °C (D796Y variation). These Tms were distinct from those of SARS-CoV-2 isolate Wuhan-Hu-1, and were specific to the Omicron variant. In HRM analyses that detected the N440K and G446S variations, the Tms of amplicons from the cDNA of the Omicron variant BA.1 and subvariant BA.2 RBDs were 73.0 °C (N440K and G446S variations) and 73.5 °C (G446S variation). This difference indicates that the SARS-CoV-2 Omicron variants BA.1 and BA.2 can be clearly discriminated. Our study demonstrates the usefulness of HRM analysis after two-step PCR for the discrimination of SARS-CoV-2 variants.     

The C-terminus of the phage λ Orf recombinase is involved in DNA binding

Phage λ Orf substitutes for the activities of the Escherichia coli RecFOR proteins in vivo and is therefore implicated as a recombination mediator, encouraging the assembly of bacterial RecA onto single-stranded DNA (ssDNA) coated with SSB. Orf exists as a dimer in solution, associates with E. coli SSB and binds preferentially to ssDNA. To help identify interacting domains we analysed Orf and SSB proteins carrying mutations or truncations in the C-terminal region. A cluster of acidic residues at the carboxy-terminus of SSB is known to attract multiple protein partners to assist in DNA replication and repair. In this case an alternative domain must be utilized since Orf association with SSB was unaffected by an SSB113 point mutant (P176S) or removal of the last ten residues (ΔC10). Structurally the Orf C-terminus consists of a helix with a flexible tail that protrudes from each side of the dimer and could serve as a binding site for either SSB or DNA. Eliminating the six residue flexible tail (ΔC6) or the entire helix (ΔC19) had no significant impact on the Orf-SSB interaction. However, the OrfΔC6 protein exhibited reduced DNA binding, a feature shared by single amino acid substitutions within (W141F) or adjacent (R140A) to this region. The OrfΔC19 mutant bound poorly to DNA and secondary structure analysis in solution revealed that this truncation induces protein misfolding and aggregation. The results show that the carboxy-terminus of Orf is involved in nucleic acid recognition and also plays an unexpected role in maintaining structural integrity.     

SARS冠状病毒未知蛋白Orf 9b的克隆、表达和结构分析

从SARS病毒基因组中克隆Orf 9b基因并构建pET32-Orf 9b原核表达质粒,在原核细胞中表达目的重组蛋白并将其纯化,经肠激酶酶切后,得到分子量为11 kD的Orf 9b蛋白。ELISA检测表明重组Orf 9b蛋白可以与SARS康复病人血清反应而不与健康人血清反应。用圆二色光谱和红外光谱初步分析了重组Orf 9b蛋白的二级结构组成。圆二色光谱显示重组Orf 9b蛋白中α-螺旋占12.5%,β-折叠占40%,无规则卷曲占47.5%,红外光谱显示重组Orf 9b蛋白中α-螺旋占13.7%,β-折叠占47.5%,无规则卷曲占37.9%,两种检测方法对重组Orf 9b蛋白结构分析的结果基本一致,提示Orf 9b蛋白富含β折叠,而α-螺旋含量较少。获得该重组蛋白及对其结构的初步分析为进一步研究Orf 9b蛋白的功能奠定基础。     

Discovery of re-purposed drugs that slow SARS-CoV-2 replication in human cells

COVID-19 vaccines based on the Spike protein of SARS-CoV-2 have been developed that appear to be largely successful in stopping infection. However, therapeutics that can help manage the disease are still required until immunity has been achieved globally. The identification of repurposed drugs that stop SARS-CoV-2 replication could have enormous utility in stemming the disease. Here, using a nano-luciferase tagged version of the virus (SARS-CoV-2-ΔOrf7a-NLuc) to quantitate viral load, we evaluated a range of human cell types for their ability to be infected and support replication of the virus, and performed a screen of 1971 FDA-approved drugs. Hepatocytes, kidney glomerulus, and proximal tubule cells were particularly effective in supporting SARS-CoV-2 replication, which is in-line with reported proteinuria and liver damage in patients with COVID-19. Using the nano-luciferase as a measure of virus replication we identified 35 drugs that reduced replication in Vero cells and human hepatocytes when treated prior to SARS-CoV-2 infection and found amodiaquine, atovaquone, bedaquiline, ebastine, LY2835219, manidipine, panobinostat, and vitamin D3 to be effective in slowing SARS-CoV-2 replication in human cells when used to treat infected cells. In conclusion, our study has identified strong candidates for drug repurposing, which could prove powerful additions to the treatment of COVID.     

SARS冠状病毒基因结构中未知蛋白Orf7的研究

为探讨SARS-CoV未知蛋白在基因表达调控中的作用,本研究以鉴定的SARS-CoVBJ01毒株cDNA为材料,采用PCR方法扩增SARS-CoV未知蛋白基因Orf7,Orf8,Orf9,引入GFP和pDsRed作为报告基因,构建出5种表达质粒:pEGFP-Orf7,pEGFP-Orf8,pEGFP-Orf79,pcDNA3-Orf7,pcDNA3-Orf8。进行体外转染实验后,采用荧光酶标仪测量其细胞中荧光蛋白的表达量,再采用流式细胞仪进一步进行瞬时表达的分析比较。结果发现,pcDNA3-Orf7与pEGFP共转染时,pcDNA3-Orf7能够增强绿色荧光蛋白(GFP)的表达,pEGFP-Orf7与pDsRed共转染时pEGFP-Orf7可以增强红色荧光蛋白(RFP)的表达,且同时表达红色和绿色荧光的细胞数增加。本文首次证明了ORF7编码的63个氨基酸的蛋白具有反式激活作用,即为SARS-CoV自身的反式激活因子。