Mutational Effect of Some Major COVID-19 Variants on Binding of the S Protein to ACE2

Molecular Dynamics Simulations on Spike Protein Mutants Binding with Human β Defensin Type 2

Human β defensin type 2 (hBD-2), a cationic cysteine-rich peptide secreted from the human innate immune system, can bind Spike-RBD at the same site as receptor ACE2, thus blocking viral entry into ACE2-expressing cells. In order to find out the impact of CoV-2 mutations on hBD-2's antiviral activity, it is important to investigate the binding and interaction of hBD-2 with RBD mutants. All-atom molecular dynamics simulations were conducted on typical RBD mutants, including N501Y, E484K, P479S, T478I, S477N, N439K, K417N, and N501Y-E484K-K417N, binding with hBD-2. Starting from the stable binding structure of hBD-2 and wt-RBD and ClusPro and HADDOCK docking-predicted initial structures, the RBD variants bound with hBD-2 simulations were set up, and NAMD simulations were conducted. Based on the structure and dynamics analysis, it was found that most RBD variants can still form a similar number of hydrogen bonds with hBD-2, in addition to having a similar-sized buried surface area (BSA) and a similar binding interface to the RBD wildtype. However, the RBD triple mutant formed a less stable binding structure with hBD-2 than other variants. Additionally, the free energy perturbation (FEP) method was applied to calculate the contribution of key mutant residues to the binding and the free energy change caused by the mutations. The result shows that N439K, K417N, and the trimutation increase the binding free energy of RBD with hBD-2; thus, RBD should bind less stably with hBD-2. E484K decreases the binding free energy, thus it should bind more stably with hBD-2, while N501Y, S477N, T478I, and P479S almost do not change the binding free energy with hBD-2. The MM-GBSA method predicted the binding interaction energy which shows that the trimutant should be able to escape the binding with hBD-2 but N501Y should not. The result can provide insight into understanding the functional mechanism of hBD-2 combating SARS-CoV-2 mutants.

Structure, dynamics and free energy studies on the effect of point mutations on SARS-CoV-2 spike protein binding with ACE2 receptor

The ongoing COVID-19 pandemic continues to infect people worldwide, and the virus continues to evolve in significant ways which can pose challenges to the efficiency of available vaccines and therapeutic drugs and cause future pandemic. Therefore, it is important to investigate the binding and interaction of ACE2 with different RBD variants. A comparative study using all-atom MD simulations was conducted on ACE2 binding with 8 different RBD variants, including N501Y, E484K, P479S, T478I, S477N, N439K, K417N and N501Y-E484K-K417N on RBD. Based on the RMSD, RMSF, and DSSP results, overall the binding of RBD variants with ACE2 is stable, and the secondary structure of RBD and ACE2 are consistent after the point mutation. Besides that, a similar buried surface area, a consistent binding interface and a similar amount of hydrogen bonds formed between RBD and ACE2 although the exact residue pairs on the binding interface were modified. The change of binding free energy from point mutation was predicted using the free energy perturbation (FEP) method. It is found that N501Y, N439K, and K417N can strengthen the binding of RBD with ACE2, while E484K and P479S weaken the binding, and S477N and T478I have negligible effect on the binding. Point mutations modified the dynamic correlation of residues in RBD based on the dihedral angle covariance matrix calculation. Doing dynamic network analysis, a common intrinsic network community extending from the tail of RBD to central, then to the binding interface region was found, which could communicate the dynamics in the binding interface region to the tail thus to the other sections of S protein. The result can supply unique methodology and molecular insight on studying the molecular structure and dynamics of possible future pandemics and design novel drugs.

SARS-CoV-2 Delta Variant: Interplay between Individual Mutations and Their Allosteric Synergy

Since its first appearance in April 2021, B.1.617.2, also termed variant Delta, catalyzed one major worldwide wave dominating the second year of coronavirus disease 2019 (COVID-19) pandemic. Despite its quick disappearance worldwide, the strong virulence caused by a few point mutations remains an unsolved problem largely. Along with the other two sublineages, the Delta variant harbors an accumulation of Spike protein mutations, including the previously identified L452R, E484Q, and the newly emerged T478K on its receptor binding domain (RBD). We used molecular dynamics (MD) simulations, in combination with free energy perturbation (FEP) calculations, to examine the effects of two combinative mutation sets, L452R + E484Q and L452R + T478K. Our dynamic trajectories reveal an enhancement in binding affinity between mutated RBD and the common receptor protein angiotensin converting enzyme 2 (ACE2) through a net increase in the buried molecular surface area of the binary complex. This enhanced binding, mediated through Gln493, sets the same stage for all three sublineages due to the presence of L452R mutation. The other mutation component, E484Q or T478K, was found to impact the RBD-ACE2 binding and help the variant to evade several monoclonal antibodies (mAbs) in a distinct manner. Especially for L452R + T478K, synergies between mutations are mediated through a complex residual and water interaction network and further enhance its binding to ACE2. Taking together, this study demonstrates that new variants of SARS-CoV-2 accomplish both "attack" (infection) and "defense" (antibody neutralization escape) with the same "polished sword" (mutated Spike RBD).

A comprehensive update on CIDO: the community-based coronavirus infectious disease ontology

Background

The current COVID-19 pandemic and the previous SARS/MERS outbreaks of 2003 and 2012 have resulted in a series of major global public health crises. We argue that in the interest of developing effective and safe vaccines and drugs and to better understand coronaviruses and associated disease mechenisms it is necessary to integrate the large and exponentially growing body of heterogeneous coronavirus data. Ontologies play an important role in standard-based knowledge and data representation, integration, sharing, and analysis. Accordingly, we initiated the development of the community-based Coronavirus Infectious Disease Ontology (CIDO) in early 2020.

Results

As an Open Biomedical Ontology (OBO) library ontology, CIDO is open source and interoperable with other existing OBO ontologies. CIDO is aligned with the Basic Formal Ontology and Viral Infectious Disease Ontology. CIDO has imported terms from over 30 OBO ontologies. For example, CIDO imports all SARS-CoV-2 protein terms from the Protein Ontology, COVID-19-related phenotype terms from the Human Phenotype Ontology, and over 100 COVID-19 terms for vaccines (both authorized and in clinical trial) from the Vaccine Ontology. CIDO systematically represents variants of SARS-CoV-2 viruses and over 300 amino acid substitutions therein, along with over 300 diagnostic kits and methods. CIDO also describes hundreds of host-coronavirus protein-protein interactions (PPIs) and the drugs that target proteins in these PPIs. CIDO has been used to model COVID-19 related phenomena in areas such as epidemiology. The scope of CIDO was evaluated by visual analysis supported by a summarization network method. CIDO has been used in various applications such as term standardization, inference, natural language processing (NLP) and clinical data integration. We have applied the amino acid variant knowledge present in CIDO to analyze differences between SARS-CoV-2 Delta and Omicron variants. CIDO's integrative host-coronavirus PPIs and drug-target knowledge has also been used to support drug repurposing for COVID-19 treatment.

Conclusion

CIDO represents entities and relations in the domain of coronavirus diseases with a special focus on COVID-19. It supports shared knowledge representation, data and metadata standardization and integration, and has been used in a range of applications.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13326-022-00279-z.

CoVM2: Molecular Biological Data Integration of SARS-CoV-2 Proteins in a Macro-to-Micro Method

The COVID-19 pandemic has been a major public health event since 2020. Multiple variant strains of SARS-CoV-2, the causative agent of COVID-19, were detected based on the mutation sites in their sequences. These sequence mutations may lead to changes in the protein structures and affect the binding states of SARS-CoV-2 and human proteins. Experimental research on SARS-CoV-2 has accumulated a large amount of structural data and protein-protein interactions (PPIs), but the studies on the SARS-CoV-2–human PPI networks lack integration of physical associations with possible protein docking information. In addition, the docking structures of variant viral proteins with human receptor proteins are still insufficient. This study constructed SARS-CoV-2–human protein–protein interaction network with data integration methods. Crystal structures were collected to map the interaction pairs. The pairs of direct interactions and physical associations were selected and analyzed for variant docking calculations. The study examined the structures of spike (S) glycoprotein of variants Delta B.1.617.2, Omicron BA.1, and Omicron BA.2. The calculated docking structures of S proteins and potential human receptors were obtained. The study integrated binary protein interactions with 3D docking structures to fulfill an extended view of SARS-CoV-2 proteins from a macro- to micro-scale.