Dramatic Differences between the Structural Susceptibility of the S1 Pre- and S2 Postfusion States of the SARS-CoV-2 Spike Protein to External Electric Fields Revealed by Molecular Dynamics Simulations

In its prefusion state, the SARS-CoV-2 spike protein (similarly to other class I viral fusion proteins) is metastable, which is considered to be an important feature for optimizing or regulating its functions. After the binding process of its S1 subunit (S1) with ACE2, the spike protein (S) undergoes a dramatic conformational change where S1 splits from the S2 subunit, which then penetrates the membrane of the host cell, promoting the fusion of the viral and cell membranes. This results in the infection of the host cell. In a previous work, we showed-using large-scale molecular dynamics simulations-that the application of external electric fields (EFs) induces drastic changes and damage in the receptor-binding domain (RBD) of the wild-type spike protein, as well of the Alpha, Beta, and Gamma variants, leaving a structure which cannot be recognized anymore by ACE2. In this work, we first extend the study to the Delta and Omicron variants and confirm the high sensitivity and extreme vulnerability of the RBD of the prefusion state of S to moderate EF (as weak as 104 V/m), but, more importantly, we also show that, in contrast, the S2 subunit of the postfusion state of the spike protein does not suffer structural damage even if electric field intensities four orders of magnitude higher are applied. These results provide a solid scientific basis to confirm the connection between the prefusion-state metastability of the SARS-CoV-2 spike protein and its susceptibility to be damaged by EF. After the virus docks to the ACE2 receptor, the stable and robust postfusion conformation develops, which exhibits a similar resistance to EF (damage threshold higher than 108 V/m) like most globular proteins.

Mutation and codon bias analysis of the spike protein of Omicron, the recent variant of SARS-CoV-2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant is a heavily mutated virus and designated as a variant of concern. To investigate the codon usage pattern of this new variant, we performed mutation and codon bias analysis for Omicron as well as for its sub-lineages BA.1 and BA.2 and compared them with the original SARS-CoV-2 and the Delta variant sequences obtained in this study. Our results indicate that the sub-lineage BA.1 and BA.2 have up to 23 sites of difference on the spike protein, which have minimal impact on function. The Omicron variant and its sub-lineages have similar codon usage patterns and A/U ending codons appear to be preferred over G/C ending codons. The Omicron has a lower degree of codon usage bias in spite of evidence that natural selection, mutation pressure and dinucleotide abundance shape the codon usage bias of Omicron, with natural selection being more significant on BA.2 than the other sub-lineages of Omicron. The codon usage pattern of Omicron variant that we explored provides valid information for a clearer understanding of Omicron and its sub-lineages, which could find application in vaccine development and optimization.

Comprehensive analysis unveils altered binding kinetics of 5-/6-methylCytosine/adenine modifications in R2R3-DNA system

Recent studies have shown that DNA methylation is an important epigenetic marker. Two prominent forms are methylation of the C5 position of cytosine and methylation of the C6 position of adenine. Given the vital significance of DNA methylation, investigating the mechanisms that influence protein binding remains a compelling pursuit. This study used molecular dynamics simulations to investigate the binding patterns of R2R3 protein and four differentially methylated DNAs. The alanine scanning combined with interaction entropy method was used to identify key residues that respond to different methylation patterns. The order of protein binding ability to DNA is as follows: unmethylated DNA > A11 methylation (5'-A6mAC-3') (6m2A system) > A10 methylation (5'-6mAAC-3') (6m1A system) > both A10 and A11 methylation (5'-6mA6mAC-3') (6mAA system) > C12 methylation (5'-AA5mC-3') (5mC system). All methylation systems lead to the sixth α helix (H6) (residues D105 to L116) moving away from the binding interface, and in the 5mC and 6m1A systems, the third α helix (H3) (residues G54 to L65) exhibits a similar trend. When the positively charged amino acids in H3 and H6 move away from the binding interface, their electrostatic and van der Waals interactions with the negatively charged DNA are weakened. Structural changes induced by methylation contributed to the destabilization of the hydrogen bond network near the original binding site, except for the 6m2A system. Moreover, there is a positive correlation between the number of methylated sites and the probability of distorting the DNA structure. Our study explores how different methylation patterns affect binding and structural adaptability, and have implications for drug discovery and understanding diseases related to abnormal methylation.

RBD spatial orientation of the spike protein and its binding to ACE2: insight into the high infectivity of the SARS-CoV-2 Delta variant from MD simulations

The spike glycoprotein on the surface of the SARS-CoV-2 envelope plays an important role in its invasion into host cells. The binding of the spike glycoprotein RBD to the angiotensin-converting enzyme 2 (ACE2) receptor as a critical step in the spread of the virus has been explored intensively since the outbreak of COVID-19, but the high transmissibility of the virus such as the Delta variant is still not fully understood. Here, molecular simulations on the binding interactions of the wild-type spike protein and its four variants (Beta, Kappa, Delta, and Mu) with ACE2 and the antibody were performed, and the present results reveal that the residue mutations will not strengthen the binding affinity of the variant for ACE2, but remarkably influences the spatial orientation of the spike protein. Only the up-right conformational receptor binding domain (RBD) can bind ACE2, which is stabilized by the nearby RBDs in the down state, revealing that the RBD bears dual functional characteristics. The present results provide new insights into plausible mechanisms for high infectivity of the virus variants and their immune escape.

Molecular adaptations during viral epidemics

In 1977, the world witnessed both the eradication of smallpox and the beginning of the modern age of genomics. Over the following half-century, 7 epidemic viruses of international concern galvanized virologists across the globe and led to increasingly extensive virus genome sequencing. These sequencing efforts exerted over periods of rapid adaptation of viruses to new hosts, in particular, humans provide insight into the molecular mechanisms underpinning virus evolution. Investment in virus genome sequencing was dramatically increased by the unprecedented support for phylogenomic analyses during the COVID-19 pandemic. In this review, we attempt to piece together comprehensive molecular histories of the adaptation of variola virus, HIV-1 M, SARS, H1N1-SIV, MERS, Ebola, Zika, and SARS-CoV-2 to the human host. Disruption of genes involved in virus-host interaction in animal hosts, recombination including genome segment reassortment, and adaptive mutations leading to amino acid replacements in virus proteins involved in host receptor binding and membrane fusion are identified as the key factors in the evolution of epidemic viruses.

Impact of the temperature on the interactions between common variants of the SARS-CoV-2 receptor binding domain and the human ACE2

Several key mutations in the Spike protein receptor binding domain (RBD) have been identified to influence its affinity for the human Angiotensin-Converting Enzyme 2 (ACE2). Here, we perform a comparative study of the ACE2 binding to the wild type (Wuhan) RBD and some of its variants: Alpha B.1.1.7, Beta B.1.351, Delta B.1.617.2, Kappa B.1.617.1, B.1.1.7 + L452R and Omicron B.1.1.529. Using a coiled-coil mediated tethering approach of ACE2 in a novel surface plasmon resonance (SPR)-based assay, we measured interactions at different temperatures. Binding experiments at 10 °C enhanced the kinetic dissimilarities between the RBD variants and allowed a proper fit to a Langmuir 1:1 model with high accuracy and reproducibility, thus unraveling subtle differences within RBD mutants and ACE2 glycovariants. Our study emphasizes the importance of SPR-based assay parameters in the acquisition of biologically relevant data and offers a powerful tool to deepen our understanding of the role of the various RBD mutations in ACE2 interaction binding parameters.