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

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.

Variant-Specific Viral Kinetics in Acute COVID-19

Understanding variant-specific differences in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral kinetics may explain differences in transmission efficiency and provide insights on pathogenesis and prevention. We evaluated SARS-CoV-2 kinetics from nasal swabs across multiple variants (Alpha, Delta, Epsilon, Gamma) in placebo recipients of the ACTIV-2/A5401 trial. Delta variant infection led to the highest maximum viral load and shortest time from symptom onset to viral load peak. There were no significant differences in time to viral clearance across the variants. Viral decline was biphasic with first- and second-phase decays having half-lives of 11 hours and 2.5 days, respectively, with differences among variants, especially in the second phase. These results suggest that while variant-specific differences in viral kinetics exist, post-peak viral load all variants appeared to be efficiently cleared by the host. Clinical Trials Registration. NCT04518410.

Enhancement of SARS-CoV-2 receptor-binding domain activity by two microbial defensins

Peptide binders are of great interest to both basic and biomedical research due to their unique properties in manipulating protein functions in a precise spatial and temporal manner. The receptor-binding domain (RBD) of the SARS-CoV-2 Spike protein is a ligand that captures human angiotensin-converting enzyme 2 (ACE2) to initiate infection. The development of binders of RBDs has value either as antiviral leads or as versatile tools to study the functional properties of RBDs dependent on their binding positions on the RBDs. In this study, we report two microbe-derived antibacterial defensins with RBD-binding activity. These two naturally occurring binders bind wild-type RBD (WT RBD) and RBDs from various variants with moderate-to-high affinity (7.6-1,450 nM) and act as activators that enhance the ACE2-binding activity of RBDs. Using a computational approach, we mapped an allosteric pathway in WT RBD that connects its ACE2-binding sites to other distal regions. The latter is targeted by the defensins, in which a cation-π interaction could trigger the peptide-elicited allostery in RBDs. The discovery of the two positive allosteric peptides of SARS-CoV-2 RBD will promote the development of new molecular tools for investigating the biochemical mechanisms of RBD allostery.

Novel Inhibitory Role of Fenofibric Acid by Targeting Cryptic Site on the RBD of SARS-CoV-2

The emergence of the recent pandemic causing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has created an alarming situation worldwide. It also prompted extensive research on drug repurposing to find a potential treatment for SARS-CoV-2 infection. An active metabolite of the hyperlipidemic drug fenofibrate (also called fenofibric acid or FA) was found to destabilize the receptor-binding domain (RBD) of the viral spike protein and therefore inhibit its binding to human angiotensin-converting enzyme 2 (hACE2) receptor. Despite being considered as a potential drug candidate for SARS-CoV-2, FA's inhibitory mechanism remains to be elucidated. We used molecular dynamics (MD) simulations to investigate the binding of FA to the RBD of the SARS-CoV-2 spike protein and revealed a potential cryptic FA binding site. Free energy calculations were performed for different FA-bound RBD complexes. The results suggest that the interaction of FA with the cryptic binding site of RBD alters the conformation of the binding loop of RBD and effectively reduces its binding affinity towards ACE2. Our study provides new insights for the design of SARS-CoV-2 inhibitors targeting cryptic sites on the RBD of SARS-CoV-2.