Computational Hot-Spot Analysis of the SARS-CoV-2 Receptor Binding Domain/ACE2 Complex*

Deciphering ACE2-RBD binding affinity through peptide scanning: A molecular dynamics simulation approach

Rapid discovery of target information for protein-protein interactions (PPIs) is significant in drug design, diagnostics, vaccine development, antibody therapy, etc. Peptide microarray is an ideal tool for revealing epitope information of PPIs. In this work, the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) spike receptor-binding domain (RBD) and the host cell receptor angiotensin-converting enzyme 2 (ACE2) were introduced as a model to study the epitope information of RBD-specific binding to ACE2 via a combination of theoretical calculations and experimental validation. Through dock and molecular dynamics simulations, it was found that among the 22 peptide fragments that consist of RBD, #14 (YNYLYRLFRKSNLKP) has the highest binding strength. Subsequently, the experiments of peptide microarray constructed based on plasmonic materials chip also confirmed the theoretical calculation data. Compared to other methods, such as phage display technology and surface plasmon resonance (SPR), this method is rapid and cost-effective, providing insights into the investigation of pathogen invasion processes and the timely development of peptide drugs and other fields.

A tool for the cheap and rapid screening of SARS-CoV-2 variants of concern (VoCs) by Sanger sequencing

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an emerging virus that, since March 2020, has been responsible for a global and ongoing pandemic. Its rapid spread over the past nearly 3 years has caused novel variants to arise. To monitor the circulation and emergence of SARS-CoV-2 variants, surveillance systems based on nucleotide mutations are required. In this regard, we searched in the spike, ORF8, and nucleocapsid genes to detect variable sites among SARS-CoV-2 variants. We describe polymorphic genetic regions that enable us to differentiate between the Alpha, Beta, Gamma, Delta, and Omicron variants of concern (VoCs). We found 21 relevant mutations, 13 of which are unique for Omicron lineages BA.1/BA.1.1, BA.2, BA.3, BA.4, and BA.5. This genetic profile enables the discrimination between VoCs using only four reverse transcription PCR fragments and Sanger sequencing, offering a cheaper and faster alternative to whole-genome sequencing for SARS-CoV-2 surveillance. IMPORTANCE Our work describes a new (Sanger sequencing-based) screening methodology for SARS-CoV-2, performing PCR amplifications of a few target regions to detect diagnostic mutations between virus variants. Using the methodology developed in this work, we were able to discriminate between the following VoCs: Alpha, Beta, Gamma, Delta, and Omicron (BA.1/BA.1.1, BA.2, BA.3, BA.4, and BA.5). This becomes important, especially in low-income countries where current methodologies like next-generation sequencing have prohibitive costs. Furthermore, rapid detection would allow sanitary authorities to take rapid measures to limit the spread of the virus and therefore reduce the probability of new virus dispersion. With this methodological approach, 13 previously unreported diagnostic mutations among several Omicron lineages were found.

Mechanistic Origin of Different Binding Affinities of SARS-CoV and SARS-CoV-2 Spike RBDs to Human ACE2

The receptor-binding domain (RBD) of the SARS-CoV-2 spike protein (RBD_CoV2) has a higher binding affinity to the human receptor angiotensin-converting enzyme 2 (ACE2) than the SARS-CoV RBD (RBD_CoV). 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 RBD_CoV2-ACE2 complex, RBD_CoV-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 RBD_CoV2 determine the higher ACE2-binding affinity of RBD_CoV2 than of RBD_CoV. 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 RBD_CoV2-ACE2, which not only tightens the interface packing and suppresses the dynamics of RBD_CoV2-ACE2, but also enhances the ACE2-binding affinity of RBD_CoV2. 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.

Epidemiological associations with genomic variation in SARS-CoV-2

SARS-CoV-2 (CoV) is the etiological agent of the COVID-19 pandemic and evolves to evade both host immune systems and intervention strategies. We divided the CoV genome into 29 constituent regions and applied novel analytical approaches to identify associations between CoV genomic features and epidemiological metadata. Our results show that nonstructural protein 3 (nsp3) and Spike protein (S) have the highest variation and greatest correlation with the viral whole-genome variation. S protein variation is correlated with nsp3, nsp6, and 3'-to-5' exonuclease variation. Country of origin and time since the start of the pandemic were the most influential metadata associated with genomic variation, while host sex and age were the least influential. We define a novel statistic-coherence-and show its utility in identifying geographic regions (populations) with unusually high (many new variants) or low (isolated) viral phylogenetic diversity. Interestingly, at both global and regional scales, we identify geographic locations with high coherence neighboring regions of low coherence; this emphasizes the utility of this metric to inform public health measures for disease spread. Our results provide a direction to prioritize genes associated with outcome predictors (e.g., health, therapeutic, and vaccine outcomes) and to improve DNA tests for predicting disease status.