Evaluation of Interactions between SARS-CoV-2 RBD and Full-Length ACE2 with Coarse-Grained Molecular Dynamics Simulations

Unraveling the binding mechanisms of SARS-CoV-2 variants through molecular simulations

The emergence of SARS-CoV-2 variants like Delta (AY.29) and Omicron (EG.5) poses continued challenges for vaccines and therapeutics. Mutations in the viral spike protein are key in altering infectivity and immune evasion. This study uses computational modeling to investigate the molecular binding mechanisms between spike protein variants and the ACE2 host receptor. Using the MARTNI force field, coarse-grained molecular dynamics (CGMD) simulations and nudged elastic band (NEB) calculations explore spike-ACE2 interactions for the wild type, Delta variant, and Omicron variant. The simulations reveal Omicron has the strongest binding affinity at -128.35 ± 10.91 kcal/mol, followed by Delta and wild type. Key mutations in Delta and Omicron, like Q493R and Q498R, optimize electrostatic contacts, enhancing ACE2 interactions. The wild-type spike has the highest transition state energy barrier at 17.87 kcal/mol, while Delta has the lowest barrier at 9.21 kcal/mol. Despite slightly higher dual barriers, Omicron's increased binding energy lowers its overall barrier to rapidly bind ACE2. These findings provide residue-level insights into mutation effects on SARS-CoV-2 infectivity. The computational modeling elucidates mechanisms underlying spike-ACE2 binding kinetics, aiding the development of vaccines and therapies targeting emerging viral strains.

Two-Stage Recognition Mechanism of the SARS-CoV-2 Receptor-Binding Domain to Angiotensin-Converting Enzyme-2 (ACE2)

The SARS-CoV-2 virus, commonly known as COVID-19, occurred in 2019. It is a highly contagious illness with effects ranging from mild symptoms to severe illness. It is also one of the best-known pathogens since more than 200,000 scientific papers occurred in the last few years. With the publication of the SARS-CoV-2 (SARS-CoV-2-CTD) spike (S) protein in a complex with human ACE2 (hACE2) (PDB (6LZG)), the molecular analysis of one of the most crucial steps on the infection pathway was possible. The aim of this manuscript is to simulate the most widely spread mutants of SARS-CoV-2, namely Alpha, Beta, Gamma, Delta, Omicron, and the first recognized variant (natural wild type). With the wide search of the hypersurface of the potential energy performed using the UNRES force field, the intermediate state of the ACE2-RBD complex was found. R403, K/N/T417, L455, F486, Y489, F495, Y501, and Y505 played a crucial role in the protein recognition mechanism. The intermediate state cannot be very stable since it will prevent the infection cascade.

Little Brown Bats (Myotis lucifugus) Support the Binding of SARS-CoV-2 Spike and Are Likely Susceptible to SARS-CoV-2 Infection

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), believed to have originated from a bat species, can infect a wide range of non-human hosts. Bats are known to harbor hundreds of coronaviruses capable of spillover into human populations. Recent studies have shown a significant variation in the susceptibility among bat species to SARS-CoV-2 infection. We show that little brown bats (LBB) express angiotensin-converting enzyme 2 receptor and the transmembrane serine protease 2, which are accessible to and support SARS-CoV-2 binding. All-atom molecular dynamics (MD) simulations revealed that LBB ACE2 formed strong electrostatic interactions with the RBD similar to human and cat ACE2 proteins. In summary, LBBs, a widely distributed North American bat species, could be at risk of SARS-CoV-2 infection and potentially serve as a natural reservoir. Finally, our framework, combining in vitro and in silico methods, is a useful tool to assess the SARS-CoV-2 susceptibility of bats and other animal species.

Short Peptide Self-Assembly in the Martini Coarse-Grain Force Field Family

ConspectusPivotal to the success of any computational experiment is the ability to make reliable predictions about the system under study and the time required to yield these results. Biomolecular interactions is one area of research that sits in every camp of resolution vs the time required, from the quantum mechanical level to in vivo studies. At an approximate midpoint, there is coarse-grained molecular dynamics, for which the Martini force fields have become the most widely used, fast enough to simulate the entire membrane of a mitochondrion though lacking atom-specific precision. While many force fields have been parametrized to account for a specific system under study, the Martini force field has aimed at casting a wider net with more generalized bead types that have demonstrated suitability for broad use and reuse in applications from protein-graphene oxide coassembly to polysaccharides interactions.In this Account, the progressive (Martini versions 1 through 3) and peripheral (Sour Martini, constant pH, Martini Straight, Dry Martini, etc.) developmental trajectory of the Martini force field will be analyzed in terms of self-assembling systems with a focus on short (two to three amino acids) peptide self-assembly in aqueous environments. In particular, this will focus on the effects of the Martini solvent model and compare how changes in bead definitions and mapping have effects on different systems. Considerable effort in the development of Martini has been expended to reduce the "stickiness" of amino acids to better simulate proteins in bilayers. We have included in this Account a short study of dipeptide self-assembly in water, using all mainstream Martini force fields, to examine their ability to reproduce this behavior. The three most recently released versions of Martini and variations in their solvents are used to simulate in triplicate all 400 dipeptides of the 20 gene-encoded amino acids. The ability of the force fields to model the self-assembly of the dipeptides in aqueoues environments is determined by the measurement of the aggregation propensity, and additional descriptors are used to gain further insight into the dipeptide aggregates.

Conformational Dynamics of the Receptor-Binding Domain of the SARS-CoV-2 Spike Protein

Variants of SARS-CoV-2 keep emerging and causing new waves of COVID-19 around the world. Effective new approaches in drug development are based on the binding of agents, such as neutralizing monoclonal antibodies to a receptor-binding domain (RBD) of SARS-CoV-2 spike protein. However, mutations in RBD may lower the affinity of previously developed antibodies. Therefore, rapid analysis of new variants and selection of a binding partner with high affinity is of great therapeutic importance. Here, we explore a computational approach based on molecular dynamics simulations and conformational clusterization techniques for the wild-type and omicron variants of RBD. Biochemical experiments support the hypothesis of the presence of several conformational states within the RBD assembly. The development of such an approach will facilitate the selection of neutralization drugs with higher affinity based on the primary structure of the target antigen.

Effect of Disulfide Bridge on the Binding of SARS-CoV-2 Fusion Peptide to Cell Membrane: A Coarse-Grained Study

In this paper, we present the parameterization of the CAVS coarse-grained (CG) force field for 20 amino acids, and our CG simulations show that the CAVS force field could accurately predict the amino acid tendency of the secondary structure. Then, we used the CAVS force field to investigate the binding of a severe acute respiratory syndrome-associated coronavirus fusion peptide (SARS-CoV-2 FP) to a phospholipid bilayer: a long FP (FP-L) containing 40 amino acids and a short FP (FP-S) containing 26 amino acids. Our CAVS CG simulations displayed that the binding affinity of the FP-L to the bilayer is higher than that of the FP-S. We found that the FP-L interacted more strongly with membrane cholesterol than the FP-S, which should be attributed to the stable helical structure of the FP-L at the C-terminus. In addition, we discovered that the FP-S had one major and two minor membrane-bound states, in agreement with previous all-atom molecular dynamics (MD) studies. However, we found that both the C-terminal and N-terminal amino acid residues of the FP-L can strongly interact with the bilayer membrane. Furthermore, we found that the disulfide bond formed between Cys840 and Cys851 stabilized the helices of the FP-L at the C-terminus, enhancing the interaction between the FP-L and the bilayer membrane. Our work indicates that the stable helical structure is crucial for binding the SARS-CoV-2 FP to cell membranes. In particular, the helical stability of FP should have a significant influence on the FP-membrane binding.

Human Fc-Conjugated Receptor Binding Domain-Based Recombinant Subunit Vaccines with Short Linker Induce Potent Neutralizing Antibodies against Multiple SARS-CoV-2 Variants

The coronavirus disease-19 (COVID-19) pandemic has been ongoing since December 2019, with more than 6.3 million deaths reported globally as of August 2022. Despite the success of several SARS-CoV-2 vaccines, the rise in variants, some of which are resistant to the effects of vaccination, highlights the need for a so-called pan-coronavirus (universal) vaccine. Here, we performed an immunogenicity comparison of prototype vaccines containing spike protein receptor-binding domain (RBD) residues 319–541, or spike protein regions S1, S2 and S fused to a histidine-tagged or human IgG1 Fc (hFC) fragment with either a longer (six residues) or shorter (three residues) linker. While all recombinant protein vaccines developed were effective in eliciting humoral immunity, the RBD-hFc vaccine was able to generate a potent neutralizing antibody response as well as a cellular immune response. We then compared the effects of recombinant protein length and linker size on immunogenicity in vivo. We found that a longer recombinant RBD protein (residues 319–583; RBD-Plus-hFc) containing a small alanine linker (AAA) was able to trigger long-lasting, high-titer neutralizing antibodies in mice. Finally, we evaluated cross-neutralization of wild-type and mutant RBD-Plus-hFc vaccines against wild-type, Alpha, Beta, Delta and Omicron SARS-CoV-2 variants. Significantly, at the same antigen dose, wild-type RBD-Plus-hFc immune sera induced broadly neutralizing antibodies against wild-type, Alpha, Beta, Delta and Omicron variants. Taken together, our findings provide valuable information for the continued development of recombinant protein-based SARS-CoV-2 vaccines and a basic foundation for booster vaccinations to avoid reinfection with SARS-CoV-2 variants.

Genomic surveillance of SARS-CoV-2 by sequencing the RBD region using Sanger sequencing from North Kerala

Next Generation Sequencing (NGS) is the gold standard for the detection of new variants of SARS-CoV-2 including those which have immune escape properties, high infectivity, and variable severity. This test is helpful in genomic surveillance, for planning appropriate and timely public health interventions. But labs with NGS facilities are not available in small or medium research settings due to the high cost of setting up such a facility. Transportation of samples from many places to few centers for NGS testing also produces delays due to transportation and sample overload leading in turn to delays in patient management and community interventions. This becomes more important for patients traveling from hotspot regions or those suspected of harboring a new variant. Another major issue is the high cost of NGS-based tests. Thus, it may not be a good option for an economically viable surveillance program requiring immediate result generation and patient follow-up. The current study used a cost-effective facility which can be set up in a common research lab and which is replicable in similar centers with expertise in Sanger nucleotide sequencing. More samples can be processed at a time and can generate the results in a maximum of 2 days (1 day for a 24 h working lab). We analyzed the nucleotide sequence of the Receptor Binding Domain (RBD) region of SARS-CoV-2 by the Sanger sequencing using in-house developed methods. The SARS-CoV-2 variant surveillance was done during the period of March 2021 to May 2022 in the Northern region of Kerala, a state in India with a population of 36.4 million, for implementing appropriate timely interventions. Our findings broadly agree with those from elsewhere in India and other countries during the period.

All-Atom Simulations of Human ACE2-Spike Protein RBD Complexes for SARS-CoV-2 and Some of its Variants: Nature of Interactions and Free Energy Diagrams for Dissociation of the Protein Complexes

The spike protein of SARS-CoV-2 is known to interact with the human ACE2 protein via its receptor binding domain (RBD). We have investigated the molecular nature of this interprotein interaction and the associated free energy diagrams for the unbinding of the two proteins for SARS-CoV-2 and some of its known variants through all-atom simulations. The present work involves generation and analysis of 2.5 μs of unbiased and 4.2 μs of biased molecular dynamics trajectories in total for five explicitly solvated RBD-ACE2 systems at full atomic level. First, we have made a comparative analysis of the details of residue-wise specific interactions of the spike protein with ACE2 for SARS-CoV-1 and SARS-CoV-2. It is found that the average numbers of both direct interprotein and water-bridged hydrogen bonds between the RBD and ACE2 are higher for SARS-CoV-2 than SARS-CoV-1. These higher hydrogen bonded interactions are further aided by enhanced nonspecific electrostatic attractions between the two protein surfaces for SARS-CoV-2. The free energy calculations reveal that there is an increase in the free energy barrier by ∼1.5 kcal/mol for the unbinding of RBD from ACE2 for SARS-CoV-2 compared to that for SARS-CoV-1. Subsequently, we considered the RBDs of three variants of SARS-CoV-2, namely N501Y, E484Q/L452R, and N440K. The free energy barrier of protein unbinding for the N501Y variant is found to be ∼4 kcal/mol higher than the wild type SARS-CoV-2 which can be attributed to additional specific interactions involving Tyr501 of RBD and Lys353 and Tyr42 of ACE2 and also enhanced nonspecific electrostatic interaction between the protein surfaces. For the other two mutant variants of E484Q/L452R and N440K, the free energy barrier for protein unbinding increases by ∼2 and ∼1 kcal/mol, respectively, compared with the wild type SARS-CoV-2, which can be attributed to an increase in the number of interprotein hydrogen bonds for the former and also to enhanced positive electrostatic potential on the RBD surfaces for both of the variants. The successive breaking of interprotein hydrogen bonds along the free energy pathway of the unbinding process is also found out for all five systems studied here.