Atomistic simulation reveals structural mechanisms underlying D614G spike glycoprotein-enhanced fitness in SARS-COV-2

D614G spike glycoprotein (sgp) mutation in rapidly spreading severe acute respiratory syndrome coronavirus-2 (SARS-COV-2) is associated with enhanced fitness and higher transmissibility in new cases of COVID-19 but the underlying mechanism is unknown. Here, using atomistic simulation, a plausible mechanism has been delineated. In G614 sgp but not wild type, increased D(G)614-T859 Cα-distance within 65 ns is interpreted as S1/S2 protomer dissociation. Overall, ACE2-binding, post-fusion core, open-state and sub-optimal antibody-binding conformations were preferentially sampled by the G614 mutant, but not wild type. Furthermore, in the wild type, only one of the three sgp chains has optimal communication route between residue 614 and the receptor-binding domain (RBD); whereas, two of the three chains communicated directly in G614 mutant. These data provide evidence that D614G sgp mutant is more available for receptor binding, cellular invasion and reduced antibody interaction; thus, providing framework for enhanced fitness and higher transmissibility in D614G SARS-COV-2 mutant.     

In silico study of novel niclosamide derivatives,SARS-CoV-2 nonstructural proteins catalytic residue-targeting small molecules drug candidates

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2)-mediated coronavirus disease 2019 (COVID-19) infection remains a global pandemic and health emergency with overwhelming social and economic impacts throughout the world. Therapeutics for COVID-19 are limited to only remdesivir; therefore, there is a need for combined, multidisciplinary efforts to develop new therapeutic molecules and explore the effectiveness of existing drugs against SARS-CoV-2. In the present study, we reported eight (SCOV-L-02, SCOV-L-09, SCOV-L-10, SCOV-L-11, SCOV-L-15, SCOV-L-18, SCOV-L-22, and SCOV-L-23) novel structurally related small-molecule derivatives of niclosamide (SCOV-L series) for their targeting potential against angiotensin-converting enzyme-2 (ACE2), type II transmembrane serine protease (TMPRSS2), and SARS-COV-2 nonstructural proteins (NSPs) including NSP5 (3CLpro), NSP3 (PLpro), and RdRp. Our correlation analysis suggested that ACE2 and TMPRSS2 modulate host immune response via regulation of immune-infiltrating cells at the site of tissue/organs entries. In addition, we identified some TMPRSS2 and ACE2 microRNAs target regulatory networks in SARS-CoV-2 infection and thus open up a new window for microRNAs-based therapy for the treatment of SARS-CoV-2 infection. Our in vitro study revealed that with the exception of SCOV-L-11 and SCOV-L-23 which were non-active, the SCOV-L series exhibited strict antiproliferative activities and non-cytotoxic effects against ACE2- and TMPRSS2-expressing cells. Our molecular docking for the analysis of receptor-ligand interactions revealed that SCOV-L series demonstrated high ligand binding efficacies (at higher levels than clinical drugs) against the ACE2, TMPRSS2, and SARS-COV-2 NSPs. SCOV-L-18, SCOV-L-15, and SCOV-L-09 were particularly found to exhibit strong binding affinities with three key SARS-CoV-2’s proteins: 3CLpro, PLpro, and RdRp. These compounds bind to the several catalytic residues of the proteins, and satisfied the criteria of drug-like candidates, having good adsorption, distribution, metabolism, excretion, and toxicity (ADMET) pharmacokinetic profile. Altogether, the present study suggests the therapeutic potential of SCOV-L series for preventing and managing SARs-COV-2 infection and are currently under detailed investigation in our lab.     

血管紧张素转换酶与抑制肽结合模式的分子动力学研究

应用分子模拟方法研究了血管紧张素转换酶(Angiotensin-converting enzyme,ACE)C端结构域(C-domain)与两种抑制肽(RIGLF/AHEPVK)的结合机制,预测了两个体系的结合模式,提出在C-domain-RIGLF中His353,Asp377,Asp453,Phe457,His513,Tyr523和Phe527为RIGLF主要结合残基,而在C-domainAHEPVK中Gln281,His353,Ser355,Glu384,Lys511,His513和Tyr523等残基起关键作用.应用结合自由能计算比较了两个体系的结合能力,结果表明,RIGLF和AHEPVK均与C-domain活性位点残基存在较强作用,且AHEPVK对C-domain的结合能力较强,与实验结果一致.     

In Silico Discovery of Antimicrobial Peptides as an Alternative to Control SARS-CoV-2

A serious pandemic has been caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The interaction between spike surface viral protein (Sgp) and the angiotensin-converting enzyme 2 (ACE2) cellular receptor is essential to understand the SARS-CoV-2 infectivity and pathogenicity. Currently, no drugs are available to treat the infection caused by this coronavirus and the use of antimicrobial peptides (AMPs) may be a promising alternative therapeutic strategy to control SARS-CoV-2. In this study, we investigated the in silico interaction of AMPs with viral structural proteins and host cell receptors. We screened the antimicrobial peptide database (APD3) and selected 15 peptides based on their physicochemical and antiviral properties. The interactions of AMPs with Sgp and ACE2 were performed by docking analysis. The results revealed that two amphibian AMPs, caerin 1.6 and caerin 1.10, had the highest affinity for Sgp proteins while interaction with the ACE2 receptor was reduced. The effective AMPs interacted particularly with Arg995 located in the S2 subunits of Sgp, which is key subunit that plays an essential role in viral fusion and entry into the host cell through ACE2. Given these computational findings, new potentially effective AMPs with antiviral properties for SARS-CoV-2 were identified, but they need experimental validation for their therapeutic effectiveness.     

Coronaviruses Use ACE2 Monomers as Entry-Receptors

The angiotensin-converting enzyme 2 (ACE2) has been identified as entry receptor on cells enabling binding and infection with the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) via trimeric spike (S) proteins protruding from the viral surface. It has been suggested that trimeric S proteins preferably bind to plasma membrane areas with high concentrations of possibly multimeric ACE2 receptors to achieve a higher binding and infection efficiency. Here we used direct stochastic optical reconstruction microscopy (dSTORM) in combination with different labeling approaches to visualize the distribution and quantify the expression of ACE2 on different cells. Our results reveal that endogenous ACE2 receptors are present as monomers in the plasma membrane with densities of only 1–2 receptors μm⁻². In addition, binding of trimeric S proteins does not induce the formation of ACE2 oligomers in the plasma membrane. Supported by infection studies using vesicular stomatitis virus (VSV) particles bearing S proteins our data demonstrate that a single S protein interaction per virus particle with a monomeric ACE2 receptor is sufficient for infection, which provides SARS-CoV-2 a high infectivity.     

A molecular docking model of SARS-CoV S1 protein in complex with its receptor, human ACE2

The exact residues within severe acute respiratory syndrome coronavirus (SARS-CoV) S1 protein and its receptor, human ACE2, involved in their interaction still remain largely undetermined. Identification of exact amino acid residues that are crucial for the interaction of S1 with ACE2 could provide working hypotheses for experimental studies and might be helpful for the development of antiviral inhibitor. In this paper, a molecular docking model of SARS-CoV S1 protein in complex with human ACE2 was constructed. The interacting residue pairs within this complex model and their contact types were also identified. Our model, supported by significant biochemical evidence, suggested receptor-binding residues were concentrated in two segments of S1 protein. In contrast, the interfacial residues in ACE2, though close to each other in tertiary structure, were found to be widely scattered in the primary sequence. In particular, the S1 residue ARG453 and ACE2 residue LYS341 might be the key residues in the complex formation.     

Blocking Effect of Demethylzeylasteral on the Interaction between Human ACE2 Protein and SARS-CoV-2 RBD Protein Discovered Using SPR Technology

The novel coronavirus disease (2019-nCoV) has been affecting global health since the end of 2019, and there is no sign that the epidemic is abating. Targeting the interaction between the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein and the human angiotensin-converting enzyme 2 (ACE2) receptor is a promising therapeutic strategy. In this study, surface plasmon resonance (SPR) was used as the primary method to screen a library of 960 compounds. A compound 02B05 (demethylzeylasteral, CAS number: 107316-88-1) that had high affinities for S-RBD and ACE2 was discovered, and binding affinities (K_D, μM) of 02B05-ACE2 and 02B05-S-RBD were 1.736 and 1.039 μM, respectively. The results of a competition experiment showed that 02B05 could effectively block the binding of S-RBD to ACE2 protein. Furthermore, pseudovirus infection assay revealed that 02B05 could inhibit entry of SARS-CoV-2 pseudovirus into 293T cells to a certain extent at nontoxic concentration. The compoundobtained in this study serve as references for the design of drugs which have potential in the treatment of COVID-19 and can thus accelerate the process of developing effective drugs to treat SARS-CoV-2 infections.     

Glycyrrhizic Acid Inhibits SARS-CoV-2 Infection by Blocking Spike Protein-Mediated Cell Attachment

Glycyrrhizic acid (GA), also known as glycyrrhizin, is a triterpene glycoside isolated from plants of Glycyrrhiza species (licorice). GA possesses a wide range of pharmacological and antiviral activities against enveloped viruses including severe acute respiratory syndrome (SARS) virus. Since the S protein (S) mediates SARS coronavirus 2 (SARS-CoV-2) cell attachment and cell entry, we assayed the GA effect on SARS-CoV-2 infection using an S protein-pseudotyped lentivirus (Lenti-S). GA treatment dose-dependently blocked Lenti-S infection. We showed that incubation of Lenti-S virus, but not the host cells with GA prior to the infection, reduced Lenti-S infection, indicating that GA targeted the virus for infection. Surface plasmon resonance measurement showed that GA interacted with a recombinant S protein and blocked S protein binding to host cells. Autodocking analysis revealed that the S protein has several GA-binding pockets including one at the interaction interface to the receptor angiotensin-converting enzyme 2 (ACE2) and another at the inner side of the receptor-binding domain (RBD) which might impact the close-to-open conformation change of the S protein required for ACE2 interaction. In addition to identifying GA antiviral activity against SARS-CoV-2, the study linked GA antiviral activity to its effect on virus cell binding.     

Effect of peptide-based captopril analogues on angiotensin converting enzyme activity and peroxynitrite-mediated tyrosine nitration

Angiotensin converting enzyme (ACE) regulates the blood pressure by converting angiotensin I to angiotensin II and bradykinin to bradykinin 1-7. These two reactions elevate the blood pressure as angiotensin II and bradykinin are vasoconstrictory and vasodilatory hormones, respectively. Therefore, inhibition of ACE is an important strategy for the treatment of hypertension. The natural substrates of ACE, i.e., angiotensin II and bradykinin, contain a Pro-Phe motif near the site of hydrolysis. Therefore, there may be a Pro-Phe binding pocket at the active site of ACE, which may facilitate the substrate binding. In view of this, we have synthesized a series of thiol- and selenol-containing dipeptides and captopril analogues and studied their ACE inhibition activities. This study reveals that both the selenol or thiol moiety and proline residues are essential for ACE inhibition. Although the introduction of a Phe residue to captopril and its selenium analogue considerably reduces the inhibitory effect, there appears to be a Phe binding pocket at the active site of ACE.     

Virtual Screening of Natural Chemical Databases to Search for Potential ACE2 Inhibitors

The angiotensin-converting enzyme II (ACE2) is a multifunctional protein in both health and disease conditions, which serves as a counterregulatory component of RAS function in a cardioprotective role. ACE2 modulation may also have relevance to ovarian cancer, diabetes, acute lung injury, fibrotic diseases, etc. Furthermore, since the outbreak of the coronavirus disease in 2019 (COVID-19), ACE2 has been recognized as the host receptor of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The receptor binding domain of the SARS-CoV-2 S-protein has a strong interaction with ACE2, so ACE2 may be a potent drug target to prevent the virus from invading host cells for anti-COVID-19 drug discovery. In this study, structure- and property-based virtual screening methods were combined to filter natural product databases from ChemDiv, TargetMol, and InterBioScreen to find potential ACE2 inhibitors. The binding affinity between protein and ligands was predicted using both Glide SP and XP scoring functions and the MM-GBSA method. ADME properties were also calculated to evaluate chemical drug-likeness. Then, molecular dynamics (MD) simulations were performed to further explore the binding modes between the highest-potential compounds and ACE2. Results showed that the compounds 154-23-4 and STOCK1N-07141 possess potential ACE2 inhibition activities and deserve further study.     

In silico studies on phytochemicals to combat the emerging COVID-19 infection

The current COVID-19 pandemic, caused by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) and its variants, remains a serious health hazard globally. The SARS-CoV-2 Mpro and spike proteins, as well as the human ACE2 receptor, have previously been reported as good targets for the development of new drug leads to combat COVID-19. Various ligands, including synthetic and plant-derived small molecules, can interact with the aforementioned proteins. In this study, we investigated the interaction of eight phytochemicals, from selected medicinal plants (Aegle marmelos, Azadirachta indica, and Ocimum sanctum) commonly used in Indian traditional medicine, with SARS-CoV-2 Mpro (PDBID: 6LU7), SARS-CoV-2S spike protein (PDB ID: 6M0J) and the human ACE2 receptor (PDB ID: 6M18). All compounds were subjected to density functional theory (DFT) and frontier molecular orbitals (FMO) analysis to determine their geometry, and key electronic and energetic properties. Upon examining the interactions of the phytochemicals with the human ACE2 receptor and the SARS-CoV-2 Mpro, spike protein targets, two compounds (C-5 and C-8) were identified as the best binding ligands. These were further examined in MD simulation studies to determine the stability of the ligand–protein interactions. QSAR, pharmacokinetic and drug-likeness properties studies revealed that C-5 may be the best candidate to serve as a template for the design and development of new drugs to combat COVID-19.     

Evidence for Multiple Binding Modes in the Initial Contact Between SARS-CoV-2 Spike S1 Protein and Cell Surface Glycans**

Infection of host cells by SARS-CoV-2 begins with recognition by the virus S (spike) protein of cell surface heparan sulfate (HS), tethering the virus to the extracellular matrix environment, and causing the subunit S1-RBD to undergo a conformational change into the ‘open’ conformation. These two events promote the binding of S1-RBD to the angiotensin converting enzyme 2 (ACE2) receptor, a preliminary step toward viral-cell membrane fusion. Combining ligand-based NMR spectroscopy with molecular dynamics, oligosaccharide analogues were used to explore the interactions between S1-RBD of SARS CoV-2 and HS, revealing several low-specificity binding modes and previously unidentified potential sites for the binding of extended HS polysaccharide chains. The evidence for multiple binding modes also suggest that highly specific inhibitors will not be optimal against protein S but, rather, diverse HS-based structures, characterized by high affinity and including multi-valent compounds, may be required.     

Possible Transmission Flow of SARS-CoV-2 Based on ACE2 Features

Angiotensin-converting enzyme 2 (ACE2) is the cellular receptor for the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) that is engendering the severe coronavirus disease 2019 (COVID-19) pandemic. The spike (S) protein receptor-binding domain (RBD) of SARS-CoV-2 binds to the three sub-domains viz. amino acids (aa) 22–42, aa 79–84, and aa 330–393 of ACE2 on human cells to initiate entry. It was reported earlier that the receptor utilization capacity of ACE2 proteins from different species, such as cats, chimpanzees, dogs, and cattle, are different. A comprehensive analysis of ACE2 receptors of nineteen species was carried out in this study, and the findings propose a possible SARS-CoV-2 transmission flow across these nineteen species.     

Search of ligands for the amyloidogenic protein beta2-microglobulin by capillary electrophoresis and other techniques

Beta2-microglobulin (beta2-m) is a small amyloidogenic protein normally present on the surface of most nucleated cells and responsible for dialysis-related amyloidosis, which represents a severe complication of long-term hemodialysis. A therapeutic approach for this amyloidosis could be based on the stabilization of beta2-m through the binding to a small molecule, and consequent inhibition of protein misfolding and amyloid fibril formation. A few compounds have been described to weakly bind beta2-m, including the drug suramin. The lack of a binding site for nonpolypeptidic ligands on the beta2-m structure makes it difficult for both the identification of functional groups responsible for the binding and the search of hits to be optimized. The characterization of the binding properties of suramin for beta2-m by using three different techniques (surface plasmon resonance, affinity CE (ACE), ultrafiltration) is here described and the results obtained are compared. The common features of the chemical structures of the compounds known to bind the protein led us to select 200 sulfonated/suramin-like molecules from a wider chemical library on the basis of similarity rules, so as to possibly single out some interesting hits and to gain more information on the functional groups involved in the binding. The development of screening methods to test the compounds by using ultrafiltration and ACE is described.     

Design of an adenosine analogue that selectively improves the affinity of a mutant U1A protein for RNA

The RNA recognition motif (RRM), one of the most common RNA binding domains, contains three highly conserved aromatic amino acids that participate in stacking interactions with RNA bases. We have investigated the contribution of these highly conserved aromatic amino acids to the affinity of the complex formed between the N-terminal RRM of the U1A protein and stem loop 2 of U1 snRNA. Previously, we found that substitution of one of these conserved aromatic amino acids, Phe56, with Ala resulted in a large destabilization of the complex. Here, we have modified A6, the base in stem loop 2 RNA that stacks with Phe56, to compensate for a portion of the destabilization caused by the Phe56Ala mutation. We have designed two modified adenosines, A-3CPh and A-4CPh, in which a phenyl group is linked to the adenosine such that it may replace the phenyl group that is eliminated by the Phe56Ala mutation in the complex. We have found that incorporation of A-3CPh into stem loop 2 RNA stabilizes the complex formed with Phe56Ala by 0.6 kcal/mol, while incorporation of A-4CPh into stem loop 2 RNA stabilizes this complex by 1.8 kcal/mol. Either base modification destabilizes the wild-type complex by 0.8-0.9 kcal/mol. Experiments with other U1A mutant proteins suggest that the stabilization of the complex between the Phe56Ala U1A protein and stem loop 2 RNA is due to a specific interaction between the Phe56Ala U1A protein and A6-4CPh stem loop 2 RNA.     

Carnosine to Combat Novel Coronavirus (nCoV): Molecular Docking and Modeling to Cocrystallized Host Angiotensin-Converting Enzyme 2 (ACE2) and Viral Spike Protein

Aims: Angiotensin-converting enzyme 2 (ACE2) plays an important role in the entry of coronaviruses into host cells. The current paper described how carnosine, a naturally occurring supplement, can be an effective drug candidate for coronavirus disease (COVID-19) on the basis of molecular docking and modeling to host ACE2 cocrystallized with nCoV spike protein. Methods: First, the starting point was ACE2 inhibitors and their structure–activity relationship (SAR). Next, chemical similarity (or diversity) and PubMed searches made it possible to repurpose and assess approved or experimental drugs for COVID-19. Parallel, at all stages, the authors performed bioactivity scoring to assess potential repurposed inhibitors at ACE2. Finally, investigators performed molecular docking and modeling of the identified drug candidate to host ACE2 with nCoV spike protein. Results: Carnosine emerged as the best-known drug candidate to match ACE2 inhibitor structure. Preliminary docking was more optimal to ACE2 than the known typical angiotensin-converting enzyme 1 (ACE1) inhibitor (enalapril) and quite comparable to known or presumed ACE2 inhibitors. Viral spike protein elements binding to ACE2 were retained in the best carnosine pose in SwissDock at 1.75 Angstroms. Out of the three main areas of attachment expected to the protein–protein structure, carnosine bound with higher affinity to two compared to the known ACE2 active site. LibDock score was 92.40 for site 3, 90.88 for site 1, and inside the active site 85.49. Conclusion: Carnosine has promising inhibitory interactions with host ACE2 and nCoV spike protein and hence could offer a potential mitigating effect against the current COVID-19 pandemic.     

Deciphering the molecular details for the binding of the prion protein to main ganglioside GM1 of neuronal membranes

The prion protein (PrP) resides in lipid rafts in?vivo, and lipids modulate misfolding of the protein to infectious isoforms. Here we demonstrate that binding of recombinant PrP to model raft membranes requires the presence of ganglioside GM1. A combination of liquid- and solid-state NMR revealed the binding sites of PrP to the saccharide head group of GM1. The binding epitope for GM1 was mapped to the folded C-terminal domain of PrP, and docking simulations identified key residues in the C-terminal region of helix C and the loop between strand S2 and helix B. Crucially, this region of PrP is linked to prion resistance in?vivo, and structural changes caused by lipid binding in this region may explain the requirement for lipids in the generation of infectious prions in?vitro.     

The Two-Way Switch Role of ACE2 in the Treatment of Novel Coronavirus Pneumonia and Underlying Comorbidities

December 2019 saw the emergence of the coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), which has spread across the globe. The high infectivity and ongoing mortality of SARS-CoV-2 emphasize the demand of drug discovery. Angiotensin-converting enzyme II (ACE2) is the functional receptor for SARS-CoV-2 entry into host cells. ACE2 exists as a membrane-bound protein on major viral target pulmonary epithelial cells, and its peptidase domain (PD) interacts SARS-CoV-2 spike protein with higher affinity. Therefore, targeting ACE2 is an important pharmacological intervention for a SARS-CoV-2 infection. In this review, we described the two-way switch role of ACE2 in the treatment of novel coronavirus pneumonia and underlying comorbidities, and discussed the potential effect of the ACE inhibitor and angiotensin receptor blocker on a hypertension patient with the SARS-CoV-2 infection. In addition, we analyzed the S-protein-binding site on ACE2 and suggested that blocking hot spot-31 and hot spot-353 on ACE2 could be a therapeutic strategy for preventing the spread of SARS-CoV-2. Besides, the recombinant ACE2 protein could be another potential treatment option for SARS-CoV-2 induced acute severe lung failure. This review could provide beneficial information for the development of anti-SARS-CoV-2 agents via targeting ACE2 and the clinical usage of renin-angiotensin system (RAS) drugs for novel coronavirus pneumonia treatment.     

Native Structure-Based Peptides as Potential Protein–Protein Interaction Inhibitors of SARS-CoV-2 Spike Protein and Human ACE2 Receptor

Severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2) is a positive-strand RNA virus that causes severe respiratory syndrome in humans, which is now referred to as coronavirus disease 2019 (COVID-19). Since December 2019, the new pathogen has rapidly spread globally, with over 65 million cases reported to the beginning of December 2020, including over 1.5 million deaths. Unfortunately, currently, there is no specific and effective treatment for COVID-19. As SARS-CoV-2 relies on its spike proteins (S) to bind to a host cell-surface receptor angiotensin-converting enzyme-2(ACE2), and this interaction is proved to be responsible for entering a virus into host cells, it makes an ideal target for antiviral drug development. In this work, we design three very short peptides based on the ACE2 sequence/structure fragments, which may effectively bind to the receptor-binding domain (RBD) of S protein and may, in turn, disrupt the important virus-host protein–protein interactions, blocking early steps of SARS-CoV-2 infection. Two of our peptides bind to virus protein with affinity in nanomolar range, and as very short peptides have great potential for drug development.     

Prediction and mitigation of mutation threats to COVID-19 vaccines and antibody therapies

Antibody therapeutics and vaccines are among our last resort to end the raging COVID-19 pandemic. They, however, are prone to over 5000 mutations on the spike (S) protein uncovered by a Mutation Tracker based on over 200 000 genome isolates. It is imperative to understand how mutations will impact vaccines and antibodies in development. In this work, we first study the mechanism, frequency, and ratio of mutations on the S protein which is the common target of most COVID-19 vaccines and antibody therapies. Additionally, we build a library of 56 antibody structures and analyze their 2D and 3D characteristics. Moreover, we predict the mutation-induced binding free energy (BFE) changes for the complexes of S protein and antibodies or ACE2. By integrating genetics, biophysics, deep learning, and algebraic topology, we reveal that most of the 462 mutations on the receptor-binding domain (RBD) will weaken the binding of S protein and antibodies and disrupt the efficacy and reliability of antibody therapies and vaccines. A list of 31 antibody disrupting mutants is identified, while many other disruptive mutations are detailed as well. We also unveil that about 65% of the existing RBD mutations, including those variants recently found in the United Kingdom (UK) and South Africa, will strengthen the binding between the S protein and human angiotensin-converting enzyme 2 (ACE2), resulting in more infectious COVID-19 variants. We discover the disparity between the extreme values of RBD mutation-induced BFE strengthening and weakening of the bindings with antibodies and angiotensin-converting enzyme 2 (ACE2), suggesting that SARS-CoV-2 is at an advanced stage of evolution for human infection, while the human immune system is able to produce optimized antibodies. This discovery, unfortunately, implies the vulnerability of current vaccines and antibody drugs to new mutations. Our predictions were validated by comparison with more than 1400 deep mutations on the S protein RBD. Our results show the urgent need to develop new mutation-resistant vaccines and antibodies and to prepare for seasonal vaccinations.

Antibody therapeutics and vaccines are among our last resort to end the raging COVID-19 pandemic.