Attention:The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Thursday, October 15th until 12:00 AM ET on Friday, October 16th due to maintenance. We apologize for the inconvenience.


Title: In-silico investigation of systematic missense mutations of middle east respiratory coronavirus spike protein
Middle East Respiratory Syndrome Coronavirus (MERS-CoV) causes severe pneumonia-like symptoms and is still pose a significant threat to global public health. A key component in the virulence of MERS-CoV is the Spike (S) protein, which binds with the host membrane receptor dipeptidyl peptidase 4 (DPP4). The goal of the present investigation is to examine the effects of missense mutations in the MERS-CoV S protein on protein stability and binding affinity with DPP4 to provide insight that is useful in developing vaccines to prevent coronavirus infection. We utilized a saturation mutagenesis approach to simulate all possible mutations in the MERS-CoV full-length S, S Receptor Binding Domain (RBD) and DPP4. We found the mutations in MERS-CoV S protein residues, G552, C503, C526, N468, G570, S532, S451, S419, S465, and S435, affect protein stability. We identified key residues, G538, E513, V555, S557, L506, L507, R511, M452, D537, and S454 in the S protein RBD region are important in the binding of MERS-CoV S protein to the DPP4 receptor. We investigated the effects of MERS-CoV S protein viral mutations on protein stability and binding affinity. In addition, we studied all DPP4 mutations and found the functional substitution R336T weakens both DPP4 protein stability and S-DPP4 binding affinity. We compared the S protein structures of MERS-CoV, SARS-CoV, and SARS-CoV-2 viruses and identified the residues like C526, C383, and N468 located in equivalent positions of these viruses have effects on S protein structure. These findings provide further information on how mutations in coronavirus S proteins effect protein function.  more » « less
Award ID(s):
1924092 2000296
PAR ID:
10358831
Author(s) / Creator(s):
; ; ; ; ; ; ;
Date Published:
Journal Name:
Frontiers in Molecular Biosciences
Volume:
9
ISSN:
2296-889X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Liu, Shan-Lu (Ed.)
    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
  3. 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
  4. Infection by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) involves the attachment of the receptor-binding domain (RBD) of its spike proteins to the ACE2 receptors on the peripheral membrane of host cells. Binding is initiated by a down-to-up conformational change in the spike protein, the change that presents the RBD to the receptor. To date, computational and experimental studies that search for therapeutics have concentrated, for good reason, on the RBD. However, the RBD region is highly prone to mutations, and is therefore a hotspot for drug resistance. In contrast, we here focus on the correlations between the RBD and residues distant to it in the spike protein. This allows for a deeper understanding of the underlying molecular recognition events and prediction of the highest-effect key mutations in distant, allosteric sites, with implications for therapeutics. Also, these sites can appear in emerging mutants with possibly higher transmissibility and virulence, and preidentifying them can give clues for designing pan-coronavirus vaccines against future outbreaks. Our model, based on time-lagged independent component analysis (tICA) and protein graph connectivity network, is able to identify multiple residues that exhibit long-distance coupling with the RBD opening. Residues involved in the most ubiquitous D614G mutation and the A570D mutation of the highly contagious UK SARS-CoV-2 variant are predicted ab initio from our model. Conversely, broad-spectrum therapeutics like drugs and monoclonal antibodies can target these key distant-but-conserved regions of the spike protein. 
    more » « less
  5. Abstract A novel coronavirus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was identified as the cause of the COVID-19 pandemic that originated in China in December 2019. Although extensive research has been performed on SARS-CoV-2, the binding behavior of spike (S) protein and receptor binding domain (RBD) of SARS-CoV-2 at different environmental conditions have yet to be studied. The objective of this study is to investigate the effect of temperature, fatty acids, ions, and protein concentration on the binding behavior and rates of association and dissociation between the S protein and RBD of SARS-CoV-2 and the hydrophobic aminopropylsilane (APS) biosensors using biolayer interferometry (BLI) validated with molecular dynamics simulation. Our results suggest three conditions—high ionic concentration, presence of hydrophobic fatty acids, and low temperature—favor the attachment of S protein and RBD to hydrophobic surfaces. Increasing the temperature within an hour from 0 to 25 °C results in S protein detachment, suggesting that freezing can cause structural changes in the S protein, affecting its binding kinetics at higher temperature. At all the conditions, RBD exhibits lower dissociation capabilities than the full-length S trimer protein, indicating that the separated RBD formed stronger attachment to hydrophobic surfaces compared to when it was included in the S protein. 
    more » « less