An official website of the United States government
We integrate evolutionary predictions based on the neutral theory of molecular evolution with protein dynamics to generate mechanistic insight into the molecular adaptations of the SARS-COV-2 spike (S) protein. With this approach, we first identified candidate adaptive polymorphisms (CAPs) of the SARS-CoV-2 S protein and assessed the impact of these CAPs through dynamics analysis. Not only have we found that CAPs frequently overlap with well-known functional sites, but also, using several different dynamics-based metrics, we reveal the critical allosteric interplay between SARS-CoV-2 CAPs and the S protein binding sites with the human ACE2 (hACE2) protein. CAPs interact far differently with the hACE2 binding site residues in the open conformation of the S protein compared to the closed form. In particular, the CAP sites control the dynamics of binding residues in the open state, suggesting an allosteric control of hACE2 binding. We also explored the characteristic mutations of different SARS-CoV-2 strains to find dynamic hallmarks and potential effects of future mutations. Our analyses reveal that Delta strain-specific variants have non-additive (i.e., epistatic) interactions with CAP sites, whereas the less pathogenic Omicron strains have mostly additive mutations. Finally, our dynamics-based analysis suggests that the novel mutations observed in the Omicron strain epistatically interact with the CAP sites to help escape antibody binding.
more »
« less
This content will become publicly available on March 24, 2027
Computational and experimental identification of potential neutralizing peptides derived from human ACE2 against SARS-CoV-2 infection
ABSTRACT The human angiotensin-converting enzyme 2 (hACE2) is the primary receptor for the entry of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Some human alleles of ACE2 exhibit an improved affinity for the SARS-CoV-2 Spike protein. However, the impact of ACE2 polymorphisms on SARS-CoV-2 infection remains unclear. Our previous study predicted that G431 and S514 in the receptor-binding domain (RBD) of SARS-CoV-2 S1 domain are important for S protein stability, and that S protein residues G496 and F497 and ACE2 residues D355 and Y41 are critical for the RBD-ACE2 interaction. In this study, we explored the potential of hACE2-derived neutralizing peptides as a therapeutic strategy against SARS-CoV-2 and investigated how ACE2 polymorphisms affect RBD-ACE2 binding affinity. We applied computational saturation mutagenesis to systematically screen the binding affinity changes among all possible ACE2 missense mutations within the ACE2-Wuhan-S1 complex. Mutations at ACE2 residues D355 and Y41 were predicted to weaken binding affinity, whereas those at N330 and D30 enhanced it. We identified six ACE2 regions (19–49, 65–102, 320–333, 348–359, 378–395, and 552–563) to be vital for ACE2-RBD interaction. We synthesized peptides corresponding to these six regions and tested them using a pseudotyped viral particle system and dot blot assay. Three peptides were confirmed to bind with the S protein, and four exhibited inhibitory effects. We aligned ACE2-Wuhan-S1 and ACE2-Omicron-S1 complexes, conducted correlation analysis, and observed similar binding patterns, suggesting that these peptides also have the potential to neutralize Omicron strains.IMPORTANCESARS-CoV-2 continues its global spread. In this research, we identified six regions within ACE2 that are vital for interaction with the viral S receptor-binding domain and have the potential to neutralize SARS-CoV-2 infection. Among the six peptides derived from ACE2, three were confirmed to bind with the S protein of the Wuhan strain, and four exhibited inhibitory effects on the Wuhan strain SARS-CoV-2. We also found ACE2 residues D355 and Y41 as weakening affinity, and N330 and D30 as enhancing it. We also aligned this complex with the ACE2-Omicron-S1 complex, performed correlation analyses, and compared their patterns of stability changes upon mutations and obtained similar results, indicating that these peptides may also be effective against Omicron variants. These results provide insight into the role of ACE2 polymorphism in viral entry and suggest that hACE2-derived peptides may offer a promising therapeutic strategy against SARS-CoV-2, demonstrating strong consistency between our computational predictions and experimental outcomes.
more »
« less
- PAR ID:
- 10695055
- Editor(s):
- Liu, Shan-Lu
- Publisher / Repository:
- American Society for Microbiology
- Date Published:
- Journal Name:
- Journal of Virology
- Volume:
- 100
- Issue:
- 3
- ISSN:
- 0022-538X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Severe Acute respiratory syndrome coronavirus (SARS-CoV-1) attaches to the host cell surface to initiate the interaction between the receptor-binding domain (RBD) of its spike glycoprotein (S) and the human Angiotensin-converting enzyme (hACE2) receptor. SARS-CoV-1 mutates frequently because of its RNA genome, which challenges the antiviral development. Here, we per-formed computational saturation mutagenesis of the S protein of SARS-CoV-1 to identify the residues crucial for its functions. We used the structure-based energy calculations to analyze the effects of the missense mutations on the SARS-CoV-1 S stability and the binding affinity with hACE2. The sequence and structure alignment showed similarities between the S proteins of SARS-CoV-1 and SARS-CoV-2. Interestingly, we found that target mutations of S protein amino acids generate similar effects on their stabilities between SARS-CoV-1 and SARS-CoV-2. For example, G839W of SARS-CoV-1 corresponds to G857W of SARS-CoV-2, which decrease the stability of their S glycoproteins. The viral mutation analysis of the two different SARS-CoV-1 isolates showed that mutations, T487S and L472P, weakened the S-hACE2 binding of the 2003–2004 SARS-CoV-1 isolate. In addition, the mutations of L472P and F360S destabilized the 2003–2004 viral isolate. We further predicted that many mutations on N-linked glycosylation sites would increase the stability of the S glycoprotein. Our results can be of therapeutic importance in the design of antivirals or vaccines against SARS-CoV-1 and SARS-CoV-2.more » « less
-
Abstract The Omicron BA.1 variant emerged in late 2021 and quickly spread across the world. Compared to the earlier SARS-CoV-2 variants, BA.1 has many mutations, some of which are known to enable antibody escape. Many of these antibody-escape mutations individually decrease the spike receptor-binding domain (RBD) affinity for ACE2, but BA.1 still binds ACE2 with high affinity. The fitness and evolution of the BA.1 lineage is therefore driven by the combined effects of numerous mutations. Here, we systematically map the epistatic interactions between the 15 mutations in the RBD of BA.1 relative to the Wuhan Hu-1 strain. Specifically, we measure the ACE2 affinity of all possible combinations of these 15 mutations (2 15 = 32,768 genotypes), spanning all possible evolutionary intermediates from the ancestral Wuhan Hu-1 strain to BA.1. We find that immune escape mutations in BA.1 individually reduce ACE2 affinity but are compensated by epistatic interactions with other affinity-enhancing mutations, including Q498R and N501Y. Thus, the ability of BA.1 to evade immunity while maintaining ACE2 affinity is contingent on acquiring multiple interacting mutations. Our results implicate compensatory epistasis as a key factor driving substantial evolutionary change for SARS-CoV-2 and are consistent with Omicron BA.1 arising from a chronic infection.more » « less
-
null (Ed.)Abstract The spike (S) glycoprotein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the binding to the permissive cells. The receptor-binding domain (RBD) of SARS-CoV-2 S protein directly interacts with the human angiotensin-converting enzyme 2 (ACE2) on the host cell membrane. In this study, we used computational saturation mutagenesis approaches, including structure-based energy calculations and sequence-based pathogenicity predictions, to quantify the systemic effects of missense mutations on SARS-CoV-2 S protein structure and function. A total of 18 354 mutations in S protein were analyzed, and we discovered that most of these mutations could destabilize the entire S protein and its RBD. Specifically, residues G431 and S514 in SARS-CoV-2 RBD are important for S protein stability. We analyzed 384 experimentally verified S missense variations and revealed that the dominant pandemic form, D614G, can stabilize the entire S protein. Moreover, many mutations in N-linked glycosylation sites can increase the stability of the S protein. In addition, we investigated 3705 mutations in SARS-CoV-2 RBD and 11 324 mutations in human ACE2 and found that SARS-CoV-2 neighbor residues G496 and F497 and ACE2 residues D355 and Y41 are critical for the RBD–ACE2 interaction. The findings comprehensively provide potential target sites in the development of drugs and vaccines against COVID-19.more » « less
-
The COVID-19 pandemic was declared due to the spread of the novel coronavirus, SARS-CoV-2. Viral infection is caused by the interaction between the SARS-CoV-2 receptor binding domain (RBD) and the human ACE2 receptor (hACE2). Previous computational studies have identified repurposed small molecules that target the RBD, but very few have screened drugs in the RBD–hACE2 interface. When studies focus solely on the binding affinity between the drug and the RBD, they ignore the effect of hACE2, resulting in an incomplete analysis. We screened ACE inhibitors and previously identified SARS-CoV-2 inhibitors for binding to the RBD—hACE2 interface, and then conducted 500 ns of unrestrained molecular dynamics (MD) simulations of fosinopril, fosinoprilat, lisinopril, emodin, diquafosol, and physcion bound to the interface to assess the binding characteristics of these ligands. Based on MM-GBSA analysis, all six ligands bind favorably in the interface and inhibit the RBD–hACE2 interaction. However, when we repeat our simulation by first binding the drug to the RBD before interacting with hACE2, we find that fosinopril, fosinoprilat, and lisinopril result in a strongly interacting trimeric complex (RBD-drug-hACE2). Hydrogen bonding and pairwise decomposition analyses further suggest that fosinopril is the best RBD inhibitor. However, when lisinopril is bound, it stabilizes the trimeric complex and, therefore, is not an ideal potential drug candidate. Overall, these results reveal important atomistic interactions critical to the binding of the RBD to hACE2 and highlight the significance of including all protein partners in the evaluation of a potential drug candidate.more » « less
