2019新型冠状病毒S蛋白的结构和功能分析

运用生物信息学,预测急性呼吸系统综合症冠状病毒2(SARS-CoV-2/2019-nCoV)的基本理化性质、结构、功能和抗原表位等,为新型冠状病毒肺炎(COVID-19)的防治提供思路。应用ExPASy分析S蛋白的消光系数、不稳定系数和半衰期等理化性质;利用SignaIP v5.0分析S蛋白的信号肽;应用TMHMM分析S蛋白的跨膜区;利用NetPhos3.1在线工具预测S蛋白的磷酸化位点;应用Pfam预测S蛋白的结构域;应用PSIPRED分析S蛋白的二级结构特征;利用SWISS-MODEL构建S蛋白的三级结构;利用BLAST分析SARS-CoV-2的S蛋白与其他物种的相似性;利用MEGA软件分析2019-nCoV的S蛋白与其他物种的进化关系。S蛋白由1 273个氨基酸组成,其相对分子质量为141 178.47,等电点为6.24,含有一个跨膜区,是低亲水性分泌蛋白;S蛋白的基本组成单位为纤突蛋白,其二级结构中以无规则卷曲和螺旋结构为主,三级结构中纤突糖蛋白和ACE2复合体具有重要的意义;2019-nCoV与蝙蝠冠状病毒和SARS-CoV同源;S蛋白存在多个潜在的线性T细胞和B细胞表位,1 202～1 210位氨基酸区域的抗原性和应答频率最高。生物信息学技术有利于了解S蛋白的理化性质、结构、功能和潜在的线性T细胞表位等,可为新型冠状肺炎的研究和防治提供参考依据。     

Enhanced sampling protocol to elucidate fusion peptide opening of SARS-CoV-2 spike protein

Large-scale conformational transitions in the spike protein S2 domain are required during host-cell infection of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus. Although conventional molecular dynamics simulations have been extensively used to study therapeutic targets of SARS-CoV-2, it is still challenging to gain molecular insight into the key conformational changes because of the size of the spike protein and the long timescale required to capture these transitions. In this work, we have developed an efficient simulation protocol that leverages many short simulations, a dynamic selection algorithm, and Markov state models to interrogate the structural changes of the S2 domain. We discovered that the conformational flexibility of the dynamic region upstream of the fusion peptide in S2 is coupled to the proteolytic cleavage state of the spike protein. These results suggest that opening of the fusion peptide likely occurs on a submicrosecond timescale after cleavage at the S2′ site. Building on the structural and dynamical information gained to date about S2 domain dynamics, we provide proof of principle that a small molecule bound to a seam neighboring the fusion peptide can slow the opening of the fusion peptide, leading to a new inhibition strategy for experiments to confirm. In aggregate, these results will aid the development of drug cocktails to inhibit infections caused by SARS-CoV-2 and other coronaviruses.     

Natural selection in the evolution of SARS-CoV-2 in bats created a generalist virus and highly capable human pathogen

Virus host shifts are generally associated with novel adaptations to exploit the cells of the new host species optimally. Surprisingly, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has apparently required little to no significant adaptation to humans since the start of the Coronavirus Disease 2019 (COVID-19) pandemic and to October 2020. Here we assess the types of natural selection taking place in Sarbecoviruses in horseshoe bats versus the early SARS-CoV-2 evolution in humans. While there is moderate evidence of diversifying positive selection in SARS-CoV-2 in humans, it is limited to the early phase of the pandemic, and purifying selection is much weaker in SARS-CoV-2 than in related bat Sarbecoviruses. In contrast, our analysis detects evidence for significant positive episodic diversifying selection acting at the base of the bat virus lineage SARS-CoV-2 emerged from, accompanied by an adaptive depletion in CpG composition presumed to be linked to the action of antiviral mechanisms in these ancestral bat hosts. The closest bat virus to SARS-CoV-2, RmYN02 (sharing an ancestor about 1976), is a recombinant with a structure that includes differential CpG content in Spike; clear evidence of coinfection and evolution in bats without involvement of other species. While an undiscovered “facilitating” intermediate species cannot be discounted, collectively, our results support the progenitor of SARS-CoV-2 being capable of efficient human–human transmission as a consequence of its adaptive evolutionary history in bats, not humans, which created a relatively generalist virus.

A study of the natural origins of SARS-CoV-2 reveals very little adaptive evolution occurring since it emerged in humans, but strong evolutionary signals in the bat virus lineage from which SARS-CoV-2 arose. Evolution in bats involved lineage-specific depletion of CpG nucleotides (linked to host anti-viral molecules), and clear evidence of recombination across these lineages, supporting bat host species’ influence on the ancestral viruses.     

The polybasic insert,the RBD of the SARS-CoV-2 spike protein,and the feline coronavirus – evolved or yet to evolve

Recent research on the SARS-CoV-2 pandemic has exploded around the furin-cleavable polybasic insert PRRAR↓S, found within the spike protein. The insert and the receptor-binding domain, (RBD), are vital clues in the Sherlock Holmes-like investigation into the origin of the virus and in its zoonotic crossover. Based on comparative analysis of the whole genome and the sequence features of the insert and the RBD domain, the bat and the pangolin have been proposed as very likely intermediary hosts. In this study, using the various databases, in-house developed tools, sequence comparisons, structure-guided docking, and molecular dynamics simulations, we cautiously present a fresh, theoretical perspective on the SARS-CoV-2 virus activation and its intermediary host. They are a) the SARS-CoV-2 has not yet acquired a fully optimal furin binding site or this seemingly less optimal sequence, PRRARS, has been selected for survival; b) in structural models of furin complexed with peptides, PRRAR↓S binds less well and with distinct differences as compared to the all basic RRKRR↓S; c) these differences may be exploited for the design of virus-specific inhibitors; d) the novel polybasic insert of SARS-CoV-2 may be promiscuous enough to be cleaved by multiple enzymes of the human airway epithelium and tissues which may explain its unexpected broad tropism; e) the RBD domain of the feline coronavirus spike protein carries residues that are responsible for high-affinity binding of the SARS-CoV-2 to the ACE 2 receptor; f) en route zoonotic transfer, the virus may have passed through the domestic cat whose very human-like ACE2 receptor and furin may have played some role in optimizing the traits required for zoonotic transfer.     

Contrasting Patterns in the Early Stage of SARS-CoV-2 Evolution between Humans and Minks

One of the unique features of SARS-CoV-2 is its apparent neutral evolution during the early pandemic (before February 2020). This contrasts with the preceding SARS-CoV epidemics, where viruses evolved adaptively. SARS-CoV-2 may exhibit a unique or adaptive feature which deviates from other coronaviruses. Alternatively, the virus may have been cryptically circulating in humans for a sufficient time to have acquired adaptive changes before the onset of the current pandemic. To test the scenarios above, we analyzed the SARS-CoV-2 sequences from minks (Neovision vision) and parental humans. In the early phase of the mink epidemic (April to May 2020), nonsynonymous to synonymous mutation ratio per site in the spike protein is 2.93, indicating a selection process favoring adaptive amino acid changes. Mutations in the spike protein were concentrated within its receptor-binding domain and receptor-binding motif. An excess of high-frequency derived variants produced by genetic hitchhiking was found during the middle (June to July 2020) and late phase I (August to September 2020) of the mink epidemic. In contrast, the site frequency spectra of early SARS-CoV-2 in humans only show an excess of low-frequency mutations, consistent with the recent outbreak of the virus. Strong positive selection in the mink SARS-CoV-2 implies that the virus may not be preadapted to a wide range of hosts and illustrates how a virus evolves to establish a continuous infection in a new host. Therefore, the lack of positive selection signal during the early pandemic in humans deserves further investigation.     

The binding of heparin to spike glycoprotein inhibits SARS-CoV-2 infection by three mechanisms

Heparin, a naturally occurring glycosaminoglycan, has been found to have antiviral activity against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative virus of COVID-19. To elucidate the mechanistic basis for the antiviral activity of heparin, we investigated the binding of heparin to the SARS-CoV-2 spike glycoprotein by means of sliding window docking, molecular dynamics simulations, and biochemical assays. Our simulations show that heparin binds at long, positively charged patches on the spike glycoprotein, thereby masking basic residues of both the receptor-binding domain (RBD) and the multifunctional S1/S2 site. Biochemical experiments corroborated the simulation results, showing that heparin inhibits the furin-mediated cleavage of spike by binding to the S1/S2 site. Our simulations showed that heparin can act on the hinge region responsible for motion of the RBD between the inactive closed and active open conformations of the spike glycoprotein. In simulations of the closed spike homotrimer, heparin binds the RBD and the N-terminal domain of two adjacent spike subunits and hinders opening. In simulations of open spike conformations, heparin induces stabilization of the hinge region and a change in RBD motion. Our results indicate that heparin can inhibit SARS-CoV-2 infection by three mechanisms: by allosterically hindering binding to the host cell receptor, by directly competing with binding to host heparan sulfate proteoglycan coreceptors, and by preventing spike cleavage by furin. Furthermore, these simulations provide insights into how host heparan sulfate proteoglycans can facilitate viral infection. Our results will aid the rational optimization of heparin derivatives for SARS-CoV-2 antiviral therapy.     

新型冠状病毒进化与变异

本文分析了新型冠状病毒(SARS-CoV-2,新冠病毒)的进化来源及刺突蛋白(spike protein,S)基因的突变情况.从GenBank数据库中下载相关病毒全基因组序列及S基因序列,运用DNAMAN9.0、MEGAX等生物信息学软件,进行多序列比对,构建系统进化树,并统计S基因位点突变情况.分析结果提示,新冠病毒...     

High-throughput,single-copy sequencing reveals SARS-CoV-2 spike variants coincident with mounting humoral immunity during acute COVID-19

Tracking evolution of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) within infected individuals will help elucidate coronavirus disease 2019 (COVID-19) pathogenesis and inform use of antiviral interventions. In this study, we developed an approach for sequencing the region encoding the SARS-CoV-2 virion surface proteins from large numbers of individual virus RNA genomes per sample. We applied this approach to the WA-1 reference clinical isolate of SARS-CoV-2 passaged in vitro and to upper respiratory samples from 7 study participants with COVID-19. SARS-CoV-2 genomes from cell culture were diverse, including 18 haplotypes with non-synonymous mutations clustered in the spike NH₂-terminal domain (NTD) and furin cleavage site regions. By contrast, cross-sectional analysis of samples from participants with COVID-19 showed fewer virus variants, without structural clustering of mutations. However, longitudinal analysis in one individual revealed 4 virus haplotypes bearing 3 independent mutations in a spike NTD epitope targeted by autologous antibodies. These mutations arose coincident with a 6.2-fold rise in serum binding to spike and a transient increase in virus burden. We conclude that SARS-CoV-2 exhibits a capacity for rapid genetic adaptation that becomes detectable in vivo with the onset of humoral immunity, with the potential to contribute to delayed virologic clearance in the acute setting.     

Population Dynamics and Structural Effects at Short and Long Range Support the Hypothesis of the Selective Advantage of the G614 SARS-CoV-2 Spike Variant

SARS-CoV-2 epidemics quickly propagated worldwide, sorting virus genomic variants in newly established propagules of infections. Stochasticity in transmission within and between countries or an actual selective advantage could explain the global high frequency reached by some genomic variants. Using statistical analyses, demographic reconstructions, and molecular dynamics simulations, we show that the globally invasive G614 spike variant 1) underwent a significant demographic expansion in most countries explained neither by stochastic effects nor by overrepresentation in clinical samples, 2) increases the spike S1/S2 furin-like site conformational plasticity (short-range effect), and 3) modifies the internal motion of the receptor-binding domain affecting its cross-connection with other functional domains (long-range effect). Our results support the hypothesis of a selective advantage at the basis of the spread of the G614 variant, which we suggest may be due to structural modification of the spike protein at the S1/S2 proteolytic site, and provide structural information to guide the design of variant-specific drugs.     

Molecular Cloning and Evolutionary Analysis of Hemoglobin α-Chain Genes in Bats

Bats are the only mammals with the capacity for powered flight. When flying, they need abundant energy and oxygen. According to previous works, the hemoglobin (Hb) oxygen loading function of bats is insensitive to variations in body temperature, although different bat species have different heat sensitivity. We cloned Hb α-chain sequences from eight bat species to investigate whether they have different characteristics. We found that Hb in the bat lineages is under purifying selection, which accords with the importance of its function in bats. Three turn regions in bat Hb, however, have distinct evolutionary rates compared with those of other mammals, and the codons in these regions have an accelerated rate of evolution. These codons are under divergent selection in bats. These changes in Hb may have occurred in response to the physiological requirements of the species concerned, as adaptations to different lifestyles.     

Evidence for a mouse origin of the SARS-CoV-2 Omicron variant

The rapid accumulation of mutations in the SARS-CoV-2 Omicron variant that enabled its outbreak raises questions as to whether its proximal origin occurred in humans or another mammalian host. Here, we identified 45 point mutations that Omicron acquired since divergence from the B.1.1 lineage. We found that the Omicron spike protein sequence was subjected to stronger positive selection than that of any reported SARS-CoV-2 variants known to evolve persistently in human hosts, suggesting a possibility of host-jumping. The molecular spectrum of mutations (i.e., the relative frequency of the 12 types of base substitutions) acquired by the progenitor of Omicron was significantly different from the spectrum for viruses that evolved in human patients but resembled the spectra associated with virus evolution in a mouse cellular environment. Furthermore, mutations in the Omicron spike protein significantly overlapped with SARS-CoV-2 mutations known to promote adaptation to mouse hosts, particularly through enhanced spike protein binding affinity for the mouse cell entry receptor. Collectively, our results suggest that the progenitor of Omicron jumped from humans to mice, rapidly accumulated mutations conducive to infecting that host, then jumped back into humans, indicating an inter-species evolutionary trajectory for the Omicron outbreak.     

Engineered disulfide reveals structural dynamics of locked SARS-CoV-2 spike

The spike (S) protein of SARS-CoV-2 has been observed in three distinct pre-fusion conformations: locked, closed and open. Of these, the function of the locked conformation remains poorly understood. Here we engineered a SARS-CoV-2 S protein construct “S-R/x3” to arrest SARS-CoV-2 spikes in the locked conformation by a disulfide bond. Using this construct we determined high-resolution structures confirming that the x3 disulfide bond has the ability to stabilize the otherwise transient locked conformations. Structural analyses reveal that wild-type SARS-CoV-2 spike can adopt two distinct locked-1 and locked-2 conformations. For the D614G spike, based on which all variants of concern were evolved, only the locked-2 conformation was observed. Analysis of the structures suggests that rigidified domain D in the locked conformations interacts with the hinge to domain C and thereby restrains RBD movement. Structural change in domain D correlates with spike conformational change. We propose that the locked-1 and locked-2 conformations of S are present in the acidic high-lipid cellular compartments during virus assembly and egress. In this model, release of the virion into the neutral pH extracellular space would favour transition to the closed or open conformations. The dynamics of this transition can be altered by mutations that modulate domain D structure, as is the case for the D614G mutation, leading to changes in viral fitness. The S-R/x3 construct provides a tool for the further structural and functional characterization of the locked conformations of S, as well as how sequence changes might alter S assembly and regulation of receptor binding domain dynamics.     

Prefusion spike protein conformational changes are slower in SARS-CoV-2 than in SARS-CoV-1

Within the last 2 decades, severe acute respiratory syndrome coronaviruses 1 and 2 (SARS-CoV-1 and SARS-CoV-2) have caused two major outbreaks; yet, for reasons not fully understood, the coronavirus disease 2019 pandemic caused by SARS-CoV-2 has been significantly more widespread than the 2003 SARS epidemic caused by SARS-CoV-1, despite striking similarities between these two viruses. The SARS-CoV-1 and SARS-CoV-2 spike proteins, both of which bind to host cell angiotensin-converting enzyme 2, have been implied to be a potential source of their differential transmissibility. However, the mechanistic details of prefusion spike protein binding to angiotensin-converting enzyme 2 remain elusive at the molecular level. Here, we performed an extensive set of equilibrium and nonequilibrium microsecond-level all-atom molecular dynamics simulations of SARS-CoV-1 and SARS-CoV-2 prefusion spike proteins to determine their differential dynamic behavior. Our results indicate that the active form of the SARS-CoV-2 spike protein is more stable than that of SARS-CoV-1 and the energy barrier associated with the activation is higher in SARS-CoV-2. These results suggest that not only the receptor-binding domain but also other domains such as the N-terminal domain could play a crucial role in the differential binding behavior of SARS-CoV-1 and SARS-CoV-2 spike proteins.     

The mutational dynamics of the SARS-CoV-2 virus in serial passages in vitro

Since its outbreak in 2019, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) keeps surprising the medical community by evolving diverse immune escape mutations in a rapid and effective manner. To gain deeper insight into mutation frequency and dynamics, we isolated ten ancestral strains of SARS-CoV-2 and performed consecutive serial incubation in ten replications in a suitable and common cell line and subsequently analysed them using RT-qPCR and whole genome sequencing. Along those lines we hoped to gain fundamental insights into the evolutionary capacity of SARS-CoV-2 in vitro. Our results identified a series of adaptive genetic changes, ranging from unique convergent substitutional mutations and hitherto undescribed insertions. The region coding for spike proved to be a mutational hotspot, evolving a number of mutational changes including the already known substitutions at positions S:484 and S:501. We discussed the evolution of all specific adaptations as well as possible reasons for the seemingly inhomogeneous potential of SARS-CoV-2 in the adaptation to cell culture. The combination of serial passage in vitro with whole genome sequencing uncovers the immense mutational potential of some SARS-CoV-2 strains. The observed genetic changes of SARS-CoV-2 in vitro could not be explained solely by selectively neutral mutations but possibly resulted from the action of directional selection accumulating favourable genetic changes in the evolving variants, along the path of increasing potency of the strain. Competition among a high number of quasi-species in the SARS-CoV-2 in vitro population gene pool may reinforce directional selection and boost the speed of evolutionary change.     

Positive selection on transposase genes of insertion sequences in the Crocosphaera watsonii genome

Insertion sequences (ISs) are mobile elements that are commonly found in bacterial genomes. Here, the structural and functional diversity of these mobile elements in the genome of the cyanobacterium Crocosphaera watsonii WH8501 is analyzed. The number, distribution, and diversity of nucleotide and amino acid stretches with similarity to the transposase gene of this IS family suggested that this genome harbors many functional as well as truncated IS fragments. The selection pressure acting on full-length transposase open reading frames of these ISs suggested (i) the occurrence of positive selection and (ii) the presence of one or more positively selected codons. These results were obtained using three data sets of transposase genes from the same IS family that were collected based on the level of amino acid similarity, the presence of an inverted repeat, and the number of sequences in the data sets. Neither recombination nor ribosomal frameshifting, which may interfere with the selection analyses, appeared to be important forces in the transposase gene family. Some positively selected codons were located in a conserved domain, suggesting that these residues are functionally important. The finding that this type of selection acts on IS-carried genes is intriguing, because although ISs have been associated with the adaptation of the bacterial host to new environments, this has typically been attributed to transposition or transformation, thus involving different genomic locations. Intragenic adaptation of IS-carried genes identified here may constitute a novel mechanism associated with bacterial diversification and adaptation.     

Viral surface geometry shapes influenza and coronavirus spike evolution through antibody pressure

The evolution of circulating viruses is shaped by their need to evade antibody response, which mainly targets the viral spike. Because of the high density of spikes on the viral surface, not all antigenic sites are targeted equally by antibodies. We offer here a geometry-based approach to predict and rank the probability of surface residues of SARS spike (S protein) and influenza H1N1 spike (hemagglutinin) to acquire antibody-escaping mutations utilizing in-silico models of viral structure. We used coarse-grained MD simulations to estimate the on-rate (targeting) of an antibody model to surface residues of the spike protein. Analyzing publicly available sequences, we found that spike surface sequence diversity of the pre-pandemic seasonal influenza H1N1 and the sarbecovirus subgenus highly correlates with our model prediction of antibody targeting. In particular, we identified an antibody-targeting gradient, which matches a mutability gradient along the main axis of the spike. This identifies the role of viral surface geometry in shaping the evolution of circulating viruses. For the 2009 H1N1 and SARS-CoV-2 pandemics, a mutability gradient along the main axis of the spike was not observed. Our model further allowed us to identify key residues of the SARS-CoV-2 spike at which antibody escape mutations have now occurred. Therefore, it can inform of the likely functional role of observed mutations and predict at which residues antibody-escaping mutation might arise.     

新型冠状病毒PLpro蛋白酶结构与功能的生物信息学分析

2019年12月,由新型冠状病毒(SARS-CoV-2)引起的新型冠状病毒肺炎(COVID-19)在中国武汉暴发。SARS-CoV-2的基因组编码2种病毒蛋白酶,即木瓜样蛋白酶(Papain-like protease,PLpro)和3C样蛋白酶(3C-like protease)。其中PLpro是SARS-CoV-2复制酶复合体(RC)形成的重要调节蛋白分子,对于病毒基因组转录和复制至关重要。因此,将SARS-CoV-2 PLpro作为药物的靶点对COVID-19的治疗具有积极意义。本研究应用生物信息学工具分析新型冠状病毒的木瓜样蛋白酶的结构和功能,首先利用BLAST和BioEdit获取SARS-CoV-2 PLpro蛋白酶(SC2-PLpro)及其同源蛋白的氨基酸序列,并利用BLAST和MEGA 6.0进行同源性分析。之后,利用ProParam和Proscale分别对SC2-PLpro蛋白酶的理化性质、亲水性和疏水性进行分析。然后,通过SOMPA、ScanProsite和InterPro分别预测SC2-PLpro蛋白酶的二级结构和功能区域,进一步利用SignalP 4.0和TMHMM对SC2-PLpro蛋白酶的信号肽和跨膜区进行分析。最后,通过SWISS-MODEL对SARS-CoV-2 PLpro蛋白酶进行三级结构同源建模。结果显示,对SARS-CoV-2 PLpro蛋白酶与已报道的PLpro蛋白酶进行多序列比对后,发现SARS-CoV-2 nsp3的746~1063段氨基酸与多种冠状病毒PLpro蛋白酶氨基酸序列高度相似。同时,同源性分析发现SARS-CoV-2与蝙蝠冠状病毒的PLpro蛋白酶具有同源性,其中与QHR63299、AVP78030相似性最高。对SC2-PLpro进行理化性质预测结果显示,其由318个氨基酸所组成,为稳定亲水性蛋白。二级结构预测结果显示SC2-PLpro主要含有α-螺旋、延伸链、β-转角、无规卷曲,四种结构贯穿整条氨基酸链。进一步进行功能分析,发现其具有完整的催化三联体、锌结合域、泛素样N末端结构域,故推测该蛋白具有去泛素化的功能。然后,信号肽假说和跨膜结构域分析结果表明,SC2-PLpro既不是分泌蛋白,也不属于跨膜蛋白。本研究提示,生物信息学分析SC2-PLpro为稳定性亲水蛋白,属于非跨膜蛋白,比较保守,具有去泛素化的功能,利用此功能可以进一步规避宿主的固有免疫反应。通过制备PLpro蛋白酶小分子抑制剂,可能有助于治疗新型冠状病毒肺炎。     

Screening of inhibitors against SARS-CoV-2 spike protein and their capability to block the viral entry mechanism: A viroinformatics study

SARS-CoV-2, previously named 2019 novel coronavirus (2019-nCoV), has been associated with the global pandemic of acute respiratory distress syndrome. First reported in December 2019 in the Wuhan province of China, this new RNA virus has several folds higher transmission among humans than its other family member (SARS-CoV and MERS-CoV). The SARS-CoV-2 spike receptor-binding domain (RBD) is the region mediating the binding of the virus to host cells via Angiotensin-converting enzyme 2 (ACE2), a critical step of viral. Here in this study, we have utilized in silico approach for the virtual screening of antiviral library extracted from the Asinex database against the Receptor binding domain (RBD) of the S1 subunit of the SARS-CoV-2 spike glycoprotein. Further, the molecules were ranked based on their binding affinity against RBD, and the top 15 molecules were selected. The affinity of these selected molecules to interrupt the ACE2-Spike interaction was also studied. It was found that the chosen molecules were demonstrating excellent binding affinity against spike protein, and these molecules were also very effectively interrupting the ACE2-RBD interaction.Furthermore, molecular dynamics (MD) simulation studies were utilized to investigate the top 3 selected molecules' stability in the ACE2-RBD complexes. To the best of our knowledge, this is the first study where molecules' inhibitory potential against the Receptor binding domain (RBD) of the S1 subunit of the SARS-CoV-2 spike glycoprotein and their inhibitory potential against the ACE2-Spike has been studied. We believe that these compounds can be further tested as a potential therapeutic option against COVID-19.     

Diverse high-affinity DNA aptamers for wild-type and B.1.1.7 SARS-CoV-2 spike proteins from a pre-structured DNA library

We performed in vitro selection experiments to identify DNA aptamers for the S1 subunit of the SARS-CoV-2 spike protein (S1 protein). Using a pool of pre-structured random DNA sequences, we obtained over 100 candidate aptamers after 13 cycles of enrichment under progressively more stringent selection pressure. The top 10 sequences all exhibited strong binding to the S1 protein. Two aptamers, named MSA1 (K_d = 1.8 nM) and MSA5 (K_d = 2.7 nM), were assessed for binding to the heat-treated S1 protein, untreated S1 protein spiked into 50% human saliva and the trimeric spike protein of both the wildtype and the B.1.1.7 variant, demonstrating comparable affinities in all cases. MSA1 and MSA5 also recognized the pseudotyped lentivirus of SARS-CoV-2 with respective K_d values of 22.7 pM and 11.8 pM. Secondary structure prediction and sequence truncation experiments revealed that both MSA1 and MSA5 adopted a hairpin structure, which was the motif pre-designed into the original library. A colorimetric sandwich assay was developed using MSA1 as both the recognition element and detection element, which was capable of detecting the pseudotyped lentivirus in 50% saliva with a limit of detection of 400 fM, confirming the potential of these aptamers as diagnostic tools for COVID-19 detection.     

The landscape of antibody binding in SARS-CoV-2 infection

The search for potential antibody-based diagnostics, vaccines, and therapeutics for pandemic severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has focused almost exclusively on the spike (S) and nucleocapsid (N) proteins. Coronavirus membrane (M), ORF3a, and ORF8 proteins are humoral immunogens in other coronaviruses (CoVs) but remain largely uninvestigated for SARS-CoV-2. Here, we use ultradense peptide microarray mapping to show that SARS-CoV-2 infection induces robust antibody responses to epitopes throughout the SARS-CoV-2 proteome, particularly in M, in which 1 epitope achieved excellent diagnostic accuracy. We map 79 B cell epitopes throughout the SARS-CoV-2 proteome and demonstrate that antibodies that develop in response to SARS-CoV-2 infection bind homologous peptide sequences in the 6 other known human CoVs. We also confirm reactivity against 4 of our top-ranking epitopes by enzyme-linked immunosorbent assay (ELISA). Illness severity correlated with increased reactivity to 9 SARS-CoV-2 epitopes in S, M, N, and ORF3a in our population. Our results demonstrate previously unknown, highly reactive B cell epitopes throughout the full proteome of SARS-CoV-2 and other CoV proteins.

Profiling of antibody binding from naïve and COVID-19 convalescent human sera to the entire proteome of SARS-CoV-2 and other human, bat and pangolin coronaviruses identifies 79 B cell epitopes throughout the SARS-CoV-2 proteome, finding that the most sensitive and specific binding occurred in the membrane (M) protein, and revealing cross-reactivity patterns.