scholarly journals Computational Saturation Mutagenesis of SARS-CoV-1 Spike Glycoprotein: Stability, Binding Affinity, and Comparison With SARS-CoV-2

2021 ◽  
Vol 8 ◽  
Author(s):  
Adebiyi Sobitan ◽  
Vidhyanand Mahase ◽  
Raina Rhoades ◽  
Dejaun Williams ◽  
Dongxiao Liu ◽  
...  

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.

2021 ◽  
Author(s):  
Adebiyi Sobitan ◽  
Vidhyanand Mahase ◽  
Raina Rhoades ◽  
Dejaun Williams ◽  
Dongxiao Liu ◽  
...  

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 performed 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.


Author(s):  
Shaolei Teng ◽  
Adebiyi Sobitan ◽  
Raina Rhoades ◽  
Dongxiao Liu ◽  
Qiyi Tang

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.


Author(s):  
Shaolei Teng ◽  
Adebiyi Sobitan ◽  
Raina Rhoades ◽  
Dongxiao Liu ◽  
Qiyi Tang

ABSTRACTThe spike (S) glycoprotein of 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 3,705 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.


Author(s):  
Nabarun Chandra Das ◽  
Rajendra Kumar Labala ◽  
Ritwik Patra ◽  
Asamanja Chattoraj ◽  
Suprabhat Mukherjee

Background: The recent outbreak of novel coronavirus disease (COVID-19) pandemic caused by SARS-CoV-2 has posed a tremendous threat to mankind. The unavailability of a specific drug or vaccine has been the major concern to date. Spike (S) glycoprotein of SARS-CoV-2 plays the most crucial role in the viral infection and immunopathogenesis, and hence this protein appears to be an efficacious target for drug discovery. Objective: Identifying potent bioactive phytocompound that can target viral spike (S) glycoprotein and human TLR4 to reduce immunopathological manifestations of COVID-19. Method: A series of thirty (30) bioactive phytocompounds, previously documented for antiviral activity, were theoretically screened for their binding efficacy against key proteins related to pathogenesis of SARSCoV-2 namely viral spike (S) glycoprotein and human TLR4. MD simulation was employed to verify the postulations of molecular docking study and further ADME analysis was performed to predict the most effective one. Results: Studies hypothesized that two new phytochemicals viz. cajaninstilbene acid (-8.83 kcal/mol) and papaverine (-5.81 kcal/mol) might be the potent inhibitors of spike glycoprotein with stout binding affinity and favourable ADME attributes. MD simulation further ratified the stability of the docked complexes between the phytochemicals and S protein through strong hydrogen bonding. Our in silico data also indicated that cajaninstilbene acid and papaverine might block human TLR4 which could be useful to mitigate SARS-CoV-2-induced lethal proinflammatory responses. Conclusion: Experimental data collectively predict cajaninstilbene acid as the potential blocker of S protein which may be used as anti-viral against COVID-19 in the future. However, further experimental validations alongside toxicological detailing are needed for claiming the candidature of these molecules as future anti-corona therapeutics.


2021 ◽  
Vol 22 (12) ◽  
pp. 6490
Author(s):  
Olga A. Postnikova ◽  
Sheetal Uppal ◽  
Weiliang Huang ◽  
Maureen A. Kane ◽  
Rafael Villasmil ◽  
...  

The SARS-CoV-2 Spike glycoprotein (S protein) acquired a unique new 4 amino acid -PRRA- insertion sequence at amino acid residues (aa) 681–684 that forms a new furin cleavage site in S protein as well as several new glycosylation sites. We studied various statistical properties of the -PRRA- insertion at the RNA level (CCUCGGCGGGCA). The nucleotide composition and codon usage of this sequence are different from the rest of the SARS-CoV-2 genome. One of such features is two tandem CGG codons, although the CGG codon is the rarest codon in the SARS-CoV-2 genome. This suggests that the insertion sequence could cause ribosome pausing as the result of these rare codons. Due to population variants, the Nextstrain divergence measure of the CCU codon is extremely large. We cannot exclude that this divergence might affect host immune responses/effectiveness of SARS-CoV-2 vaccines, possibilities awaiting further investigation. Our experimental studies show that the expression level of original RNA sequence “wildtype” spike protein is much lower than for codon-optimized spike protein in all studied cell lines. Interestingly, the original spike sequence produces a higher titer of pseudoviral particles and a higher level of infection. Further mutagenesis experiments suggest that this dual-effect insert, comprised of a combination of overlapping translation pausing and furin sites, has allowed SARS-CoV-2 to infect its new host (human) more readily. This underlines the importance of ribosome pausing to allow efficient regulation of protein expression and also of cotranslational subdomain folding.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Shin Irumagawa ◽  
Kaito Kobayashi ◽  
Yutaka Saito ◽  
Takeshi Miyata ◽  
Mitsuo Umetsu ◽  
...  

AbstractThe stability of proteins is an important factor for industrial and medical applications. Improving protein stability is one of the main subjects in protein engineering. In a previous study, we improved the stability of a four-helix bundle dimeric de novo protein (WA20) by five mutations. The stabilised mutant (H26L/G28S/N34L/V71L/E78L, SUWA) showed an extremely high denaturation midpoint temperature (Tm). Although SUWA is a remarkably hyperstable protein, in protein design and engineering, it is an attractive challenge to rationally explore more stable mutants. In this study, we predicted stabilising mutations of WA20 by in silico saturation mutagenesis and molecular dynamics simulation, and experimentally confirmed three stabilising mutations of WA20 (N22A, N22E, and H86K). The stability of a double mutant (N22A/H86K, rationally optimised WA20, ROWA) was greatly improved compared with WA20 (ΔTm = 10.6 °C). The model structures suggested that N22A enhances the stability of the α-helices and N22E and H86K contribute to salt-bridge formation for protein stabilisation. These mutations were also added to SUWA and improved its Tm. Remarkably, the most stable mutant of SUWA (N22E/H86K, rationally optimised SUWA, ROSA) showed the highest Tm (129.0 °C). These new thermostable mutants will be useful as a component of protein nanobuilding blocks to construct supramolecular protein complexes.


2021 ◽  
Vol 9 (5) ◽  
pp. 1015
Author(s):  
Tianyu Zhang ◽  
Xin Gao ◽  
Dongqiang Wang ◽  
Jixue Zhao ◽  
Nan Zhang ◽  
...  

Cryptosporidium parvum is a globally recognized zoonotic parasite of medical and veterinary importance. This parasite mainly infects intestinal epithelial cells and causes mild to severe watery diarrhea that could be deadly in patients with weakened or defect immunity. However, its molecular interactions with hosts and pathogenesis, an important part in adaptation of parasitic lifestyle, remain poorly understood. Here we report the identification and characterization of a C. parvum T-cell immunomodulatory protein homolog (CpTIPH). CpTIPH is a 901-aa single-pass type I membrane protein encoded by cgd5_830 gene that also contains a short Vibrio, Colwellia, Bradyrhizobium and Shewanella (VCBS) repeat and relatively long integrin alpha (ITGA) N-terminus domain. Immunofluorescence assay confirmed the location of CpTIPH on the cell surface of C. parvum sporozoites. In congruence with the presence of VCBS repeat and ITGA domain, CpTIPH displayed high, nanomolar binding affinity to host cell surface (i.e., Kd(App) at 16.2 to 44.7 nM on fixed HCT-8 and CHO-K1 cells, respectively). The involvement of CpTIPH in the parasite invasion is partly supported by experiments showing that an anti-CpTIPH antibody could partially block the invasion of C. parvum sporozoites into host cells. These observations provide a strong basis for further investigation of the roles of CpTIPH in parasite-host cell interactions.


Genetics ◽  
2003 ◽  
Vol 163 (1) ◽  
pp. 91-101 ◽  
Author(s):  
Erin N Asleson ◽  
Dennis M Livingston

Abstract We investigated the stability of the Saccharomyces cerevisiae Rad52 protein to learn how a cell controls its quantity and longevity. We measured the cellular levels of wild-type and mutant forms of Rad52p when expressed from the RAD52 promoter and the half-lives of the various forms of Rad52p when expressed from the GAL1 promoter. The wild-type protein has a half-life of 15 min. rad52 mutations variably affect the cellular levels of the protein products, and these levels correlate with the measured half-lives. While missense mutations in the N terminus of the protein drastically reduce the cellular levels of the mutant proteins, two mutations—one a deletion of amino acids 210-327 and the other a missense mutation of residue 235—increase the cellular level and half-life more than twofold. These results suggest that Rad52p is subject to post-translational regulation. Proteasomal mutations have no effect on Rad52p half-life but increase the amount of RAD52 message. In contrast to Rad52p, the half-life of Rad51p is >2 hr, and RAD51 expression is unaffected by proteasomal mutations. These differences between Rad52p and Rad51p suggest differential regulation of two proteins that interact in recombinational repair.


2021 ◽  
Vol 22 (5) ◽  
pp. 2732
Author(s):  
Nadine Reichhart ◽  
Vladimir M. Milenkovic ◽  
Christian H. Wetzel ◽  
Olaf Strauß

The anoctamin (TMEM16) family of transmembrane protein consists of ten members in vertebrates, which act as Ca2+-dependent ion channels and/or Ca2+-dependent scramblases. ANO4 which is primarily expressed in the CNS and certain endocrine glands, has been associated with various neuronal disorders. Therefore, we focused our study on prioritizing missense mutations that are assumed to alter the structure and stability of ANO4 protein. We employed a wide array of evolution and structure based in silico prediction methods to identify potentially deleterious missense mutations in the ANO4 gene. Identified pathogenic mutations were then mapped to the modeled human ANO4 structure and the effects of missense mutations were studied on the atomic level using molecular dynamics simulations. Our data show that the G80A and A500T mutations significantly alter the stability of the mutant proteins, thus providing new perspective on the role of missense mutations in ANO4 gene. Results obtained in this study may help to identify disease associated mutations which affect ANO4 protein structure and function and might facilitate future functional characterization of ANO4.


Sign in / Sign up

Export Citation Format

Share Document