Computational modelling of potentially emerging SARS-CoV-2 spike protein RBDs mutations with higher binding affinity towards ACE2: A structural modelling study.
Ontology highlight
ABSTRACT:
SUBMITTER: Khan A
PROVIDER: S-EPMC8716154 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has caused substantially more infections, deaths, and economic disruptions than the 2002-2003 SARS-CoV. The key to understanding SARS-CoV-2's higher infectivity lies partly in its host receptor recognition mechanism. Experiments show that the human angiotensin converting enzyme 2 (ACE2) protein, which serves as the primary receptor for both CoVs, binds to the receptor binding domain (RBD) of CoV-2's spike protein stronger than SARS-CoV's spike RBD. The molecular basis for this difference in binding affinity, however, remains unexplained from X-ray structures. To go beyond insights gained from X-ray structures and investigate the role of thermal fluctuations in structure, we employ all-atom molecular dynamics simulations. Microseconds-long simulations reveal that while CoV and CoV-2 spike-ACE2 interfaces have similar conformational binding modes, CoV-2 spike interacts with ACE2 via a larger combinatorics of polar contacts, and on average, makes 45% more polar contacts. Correlation analysis and thermodynamic calculations indicate that these differences in the density and dynamics of polar contacts arise from differences in spatial arrangements of interfacial residues, and dynamical coupling between interfacial and non-interfacial residues. These results recommend that ongoing efforts to design spike-ACE2 peptide blockers will benefit from incorporating dynamical information as well as allosteric coupling effects.
| S-EPMC8250905 | biostudies-literature
Project description:The receptor-binding domain (RBD) of the SARS-CoV-2 spike protein (RBDCoV2) has a higher binding affinity to the human receptor angiotensin-converting enzyme 2 (ACE2) than the SARS-CoV RBD (RBDCoV). Here, we performed molecular dynamics (MD) simulations, binding free energy (BFE) calculations, and interface residue contact network (IRCN) analysis to explore the mechanistic origin of different ACE2-binding affinities of the two RBDs. The results demonstrate that, when compared to the RBDCoV2-ACE2 complex, RBDCoV-ACE2 features enhanced dynamicsand inter-protein positional movements and increased conformational entropy and conformational diversity. Although the inter-protein electrostatic attractive interactions are the primary determinant for the high ACE2-binding affinities of both RBDs, the significantly enhanced electrostatic attractive interactions between ACE2 and RBDCoV2 determine the higher ACE2-binding affinity of RBDCoV2 than of RBDCoV. Comprehensive comparative analyses of the residue BFE components and IRCNs between the two complexes reveal that it is the residue changes at the RBD interface that lead to the overall stronger inter-protein electrostatic attractive force in RBDCoV2-ACE2, which not only tightens the interface packing and suppresses the dynamics of RBDCoV2-ACE2, but also enhances the ACE2-binding affinity of RBDCoV2. Since the RBD residue changes involving gain/loss of the positively/negatively charged residues can greatly enhance the binding affinity, special attention should be paid to the SARS-CoV-2 variants carrying such mutations, particularly those near or at the binding interfaces with the potential to form hydrogen bonds and/or salt bridges with ACE2.
| S-EPMC9032924 | biostudies-literature
Project description:Coronaviruses (CoVs) pose a major risk to global public health due to their ability to infect diverse animal species and potential for emergence in humans. The CoV spike protein mediates viral entry into the cell and plays a crucial role in determining the binding affinity to host cell receptors. With particular emphasis on α- and β-coronaviruses that infect humans and domestic animals, current research on CoV receptor use suggests that the exploitation of the angiotensin-converting enzyme 2 (ACE2) receptor poses a significant threat for viral emergence with pandemic potential. This review summarizes the approaches used to study binding interactions between CoV spike proteins and the human ACE2 (hACE2) receptor. Solid-phase enzyme immunoassays and cell binding assays allow qualitative assessment of binding but lack quantitative evaluation of affinity. Surface plasmon resonance, Bio-layer interferometry, and Microscale Thermophoresis on the other hand, provide accurate affinity measurement through equilibrium dissociation constants (KD). In silico modeling predicts affinity through binding structure modeling, protein-protein docking simulations, and binding energy calculations but reveals inconsistent results due to the lack of a standardized approach. Machine learning and deep learning models utilize simulated and experimental protein-protein interaction data to elucidate the critical residues associated with CoV binding affinity to hACE2. Further optimization and standardization of existing approaches for studying binding affinity could aid pandemic preparedness. Specifically, prioritizing surveillance of CoVs that can bind to human receptors stands to mitigate the risk of zoonotic spillover.
| S-EPMC10831124 | biostudies-literature
Project description:Initial phase of COVID-19 infection is associated with the binding of viral spike protein S1 receptor binding domain (RBD) with the host cell surface receptor, ACE2. Peptide inhibitors typically interact with spike proteins in order to block its interaction with ACE2, and this knowledge would promote the use of such peptides as therapeutic scaffolds. The present study examined the competitive inhibitor activity of a broad spectrum antimicrobial peptide, Dermaseptin-S4 (S4) and its analogues. Three structural S4 analogues viz., S4 (K4), S4 (K20) and S4 (K4K20) were modelled by substituting charged lysine for non-polar residues in S4 and subsequently, docked with S1. Further, the comparative analysis of inter-residue contacts and non-covalent intermolecular interactions among S1-S4 (K4), S1-S4 (K4K20) and S1-ACE2 complexes were carried out to explore their mode of binding with S1. Interestingly, S1-S4 (K4) established more inter-molecular interactions compared to S4 (K4K20) and S1-ACE2. In order to substantiate this study, the normal mode analysis (NMA) was conducted to show how the structural stability of the flexible loop region in S1 is affected by atomic displacements in unbound S1 and docked complexes. Markedly, the strong interactions consistently maintained by S1-S4 (K4) complex revealed their conformational transition over the harmonic motion period. Moreover, S1-S4 (K4) peptide complex showed a higher energy deformation profile compared to S1-S4 (K4K20), where the higher energy deformation suggests the rigidity of the docked complex and thus it's harder deformability, which is also substantiated by molecular dynamics simulation. In conclusion, S1-S4 (K4) complex has definitely exhibited a functionally significant dynamics compared to S1-ACE2 complex; this peptide inhibitor, S4 (K4) will need to be considered as the best therapeutic scaffold to block SARS-CoV-2 infection.
| S-EPMC7811342 | biostudies-literature
Project description:The association of the receptor binding domain (RBD) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein with human angiotensin-converting enzyme 2 (hACE2) represents the first required step for cellular entry. SARS-CoV-2 has continued to evolve with the emergence of several novel variants, and amino acid changes in the RBD have been implicated with increased fitness and potential for immune evasion. Reliably predicting the effect of amino acid changes on the ability of the RBD to interact more strongly with the hACE2 can help assess the implications for public health and the potential for spillover and adaptation into other animals. Here, we introduce a two-step framework that first relies on 48 independent 4-ns molecular dynamics (MD) trajectories of RBD-hACE2 variants to collect binding energy terms decomposed into Coulombic, covalent, van der Waals, lipophilic, generalized Born solvation, hydrogen bonding, π-π packing, and self-contact correction terms. The second step implements a neural network to classify and quantitatively predict binding affinity changes using the decomposed energy terms as descriptors. The computational base achieves a validation accuracy of 82.8% for classifying single-amino acid substitution variants of the RBD as worsening or improving binding affinity for hACE2 and a correlation coefficient of 0.73 between predicted and experimentally calculated changes in binding affinities. Both metrics are calculated using a fivefold cross-validation test. Our method thus sets up a framework for screening binding affinity changes caused by unknown single- and multiple-amino acid changes offering a valuable tool to predict host adaptation of SARS-CoV-2 variants toward tighter hACE2 binding.
| S-EPMC8594574 | biostudies-literature
Project description:A novel coronavirus (SARS-CoV-2) has emerged to a global pandemic and caused significant damages to public health. Human angiotensin-converting enzyme 2(ACE2) was identified as the entry receptor for SARS-CoV-2. As a carboxypeptidase, ACE2 cleaves many biological substrates besides Ang II to control vasodilatation and permeability. Given the nanomolar high affinity between ACE2 and SARS-CoV-2 spike protein, we wonder how this interaction would affect the enzymatic activity of ACE2. Surprisingly, SARS-CoV-2 trimeric spike protein increased ACE2 proteolytic activity ~3-10 fold when fluorogenic caspase-1 substrate and Bradykinin-analog peptides were used to characterize ACE2 activity. In addition, the enhancement was mediated by ACE2 binding of RBD domain of SARS-CoV-2 spike. These results highlighted the altered activity of ACE2 during SARS-CoV-2 infection and would shed new lights on the pathogenesis of COVID-19 and its complications for better treatments.
| S-EPMC7337377 | biostudies-literature
Project description:The novel severe acute respiratory syndrome coronavirus (SARS-CoV-2) has emerged to a pandemic and caused global public health crisis. Human angiotensin-converting enzyme 2(ACE2) was identified as the entry receptor for SARS-CoV-2. As a carboxypeptidase, ACE2 cleaves many biological substrates besides angiotensin II to control vasodilatation and vascular permeability. Given the nanomolar high affinity between ACE2 and SARS-CoV-2 spike protein, we investigated how this interaction would affect the enzymatic activity of ACE2. Surprisingly, SARS-CoV-2 trimeric spike protein increased ACE2 proteolytic activity ∼3-10 fold against model peptide substrates, such as caspase-1 substrate and Bradykinin-analog. The enhancement in ACE2 enzymatic function was mediated by the binding of SARS-CoV-2 spike RBD domain. These results highlighted the potential for SARS-CoV-2 infection to enhance ACE2 activity, which may be relevant to the cardiovascular symptoms associated with COVID-19.
| S-EPMC7833600 | biostudies-literature
Project description:The current pandemic of COVID-19 is caused by a novel coronavirus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The SARS-CoV-2 spike protein receptor-binding domain (RBD) is the critical determinant of viral tropism and infectivity. To investigate whether naturally occurring RBD mutations during the early transmission phase have altered the receptor binding affinity and infectivity, we first analyzed in silico the binding dynamics between SARS-CoV-2 RBD mutants and the human angiotensin-converting enzyme 2 (ACE2) receptor. Among 32,123 genomes of SARS-CoV-2 isolates (December 2019 through March 2020), 302 nonsynonymous RBD mutants were identified and clustered into 96 mutant types. The six dominant mutations were analyzed applying molecular dynamics simulations (MDS). The mutant type V367F continuously circulating worldwide displayed higher binding affinity to human ACE2 due to the enhanced structural stabilization of the RBD beta-sheet scaffold. The MDS also indicated that it would be difficult for bat SARS-like CoV to infect humans. However, the pangolin CoV is potentially infectious to humans. The increased infectivity of V367 mutants was further validated by performing receptor-ligand binding enzyme-linked immunosorbent assay (ELISA), surface plasmon resonance, and pseudotyped virus assays. Phylogenetic analysis of the genomes of V367F mutants showed that during the early transmission phase, most V367F mutants clustered more closely with the SARS-CoV-2 prototype strain than the dual-mutation variants (V367F+D614G), which may derivate from recombination. The analysis of critical RBD mutations provides further insights into the evolutionary trajectory of early SARS-CoV-2 variants of zoonotic origin under negative selection pressure and supports the continuing surveillance of spike mutations to aid in the development of new COVID-19 drugs and vaccines. IMPORTANCE A novel coronavirus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has caused the pandemic of COVID-19. The origin of SARS-CoV-2 was associated with zoonotic infections. The spike protein receptor-binding domain (RBD) is identified as the critical determinant of viral tropism and infectivity. Thus, whether mutations in the RBD of the circulating SARS-CoV-2 isolates have altered the receptor binding affinity and made them more infectious has been the research hot spot. Given that SARS-CoV-2 is a novel coronavirus, the significance of our research is in identifying and validating the RBD mutant types emerging during the early transmission phase and increasing human angiotensin-converting enzyme 2 (ACE2) receptor binding affinity and infectivity. Our study provides insights into the evolutionary trajectory of early SARS-CoV-2 variants of zoonotic origin. The continuing surveillance of RBD mutations with increased human ACE2 affinity in human or other animals is critical to the development of new COVID-19 drugs and vaccines against these variants during the sustained COVID-19 pandemic.
| S-EPMC8373230 | biostudies-literature
Project description:AimsSevere acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a novel member of the betacoronaviruses family affecting the lower respiratory tract mainly through binding to angiotensin converting enzyme 2 (ACE2) via its S-protein. Genetic analysis of (ACE2) gene revealed several variants that have been suggested to regulate the interaction with S protein. This study investigates the N720D variant, positioned in the collectrin-like domain (CLD) at proximity to type II transmembrane serine protease (TMPRSS2) cleavage site.Main methodsThe effect of N720D variant on ACE2 structure and thermodynamic stability was studied by DynaMut. HDOCK was utilised to model TMPRSS2 protease binding to ACE2 WT and D720 variant cleavage site. PRODIGY was used to calculate binding affinities and MD simulation tools calculated the at 100 ns for ACE2 apo structure and the ACE2-TMPRSS2 complex.Key findingsThe N720D variant is a more dynamic structure with a free energy change (ΔΔG): -0.470 kcal/mol. As such, introducing a tighter binding affinity of Kd = 3.2 × 10-10 M between TMPRSS2 and N720D variant. RMSD, RMSF calculations showed the N720D variant is less stable, however, RMSF values of the D720-TMPRSS2 complex reflected a slower dynamic motion.SignificanceThe hotspot N720D variant in the CLD of ACE2 affected the stability and flexibility of ACE2 by increasing the level of motion in the loop region, resulting in a more favourable site for TMPRSS2 binding and cleavage. Consequently, this would facilitate S-protein binding and can potentially increase viral entry highlighting the importance of variants affecting the ACE2-TMPRSS2 complex.
| S-EPMC7405906 | biostudies-literature
Project description:The causative agent of the COVID-19 pandemic, SARS-CoV-2, is steadily mutating during continuous transmission among humans. Such mutations can occur in the spike (S) protein that binds to the ACE2 receptor and is cleaved by TMPRSS2. However, whether S mutations affect SARS-CoV-2 cell entry remains unknown. Here, we show that naturally occurring S mutations can reduce or enhance cell entry via ACE2 and TMPRSS2. A SARS-CoV-2 S-pseudotyped lentivirus exhibits substantially lower entry than that of SARS-CoV S. Among S variants, the D614G mutant shows the highest cell entry, as supported by structural and binding analyses. Nevertheless, the D614G mutation does not affect neutralization by antisera against prototypic viruses. Taken together, we conclude that the D614G mutation increases cell entry by acquiring higher affinity to ACE2 while maintaining neutralization susceptibility. Based on these findings, further worldwide surveillance is required to understand SARS-CoV-2 transmissibility among humans.
| S-EPMC7870668 | biostudies-literature