Mutations of SARS-CoV-2 RBD May Alter Its Molecular Structure to Improve Its Infection Efficiency

Molecular dynamics of SARS-CoV-2 omicron variants from Philippine isolates against hesperidin as spike protein inhibitor

SARS-CoV-2 remains a global threat with new sublineages posing challenges, particularly in the Philippines. Hesperidin (HD) is being studied as a potential prophylactic for COVID-19. However, the virus's rapid evolution could alter how HD binds to it, affecting its effectiveness. Here, we study the mutation-induced variabilities of HD dynamics and their effects on molecular energetics in SARS-CoV-2 spike receptor complex systems. We considered eight different point mutations present in the Omicron variant. Root-mean-square deviation and binding energy analysis showed that S477N and Omicron did not eject HD throughout the simulation. Hydrogen bond distribution analysis highlighted the involvement of hydrogen bonding in mutant-HD stabilization, especially for S477N and Omicron. Root-mean-square fluctuation analysis revealed evidence of Y505H destabilization on complex systems, while distal-end loop mutations increased loop flexibility for all models bearing the three mutations. Per-residue energy decomposition demonstrated that Q493R substitution increased HD interaction. Free energy landscape and essential dynamics through principal component analysis provided insights into the conformational subspace distribution of mutant model molecular dynamics trajectories. In conclusion, significant mutations contributed to the HD interaction in different ways. S477N has shown significant binding contributions through favorable ligand interaction, while other mutations contribute via conformational modifications, increased affinity due to sidechain mutations, and increased loop flexibility.

Mechanistic insights into the orthogonal functionality of an AHL-mediated quorum-sensing circuit in Yersinia pseudotuberculosis

YpsR, a pivotal regulatory protein in the quorum-sensing (QS) of Yersinia pseudotuberculosis(Y. pstb), is essential for molecular signaling, yet its molecular mechanisms remain poorly understood. Herein, this study systematically investigates the interactions between YpsR and acyl-homoserine lactones (AHLs), shedding light on the selective mechanism of YpsR to various AHL molecules. Using molecular docking and surface plasmon resonance (SPR) analysis, we confirmed YpsR's binding affinities, with the strongest observed for 3OC6-HSL, which notably inhibited Y. pstb growth. Additionally, we engineered a whole-cell biosensor based on YpsR-AHL interaction, which exhibited sensitivity to the signal molecule 3OC6-HSL produced by Y. pstb. Furthermore, key YpsR residues (S32, Y50, W54, D67) involved in AHL binding were identified and validated. Overall, this research elucidates the mechanisms of QS signal recognition in Y. pstb, providing valuable insights that support the development of diagnostic tools for detecting Y. pstb infections.

Deep learning for discriminating non-trivial conformational changes in molecular dynamics simulations of SARS-CoV-2 spike-ACE2

Molecular dynamics (MD) simulations produce a substantial volume of high-dimensional data, and traditional methods for analyzing these data pose significant computational demands. Advances in MD simulation analysis combined with deep learning-based approaches have led to the understanding of specific structural changes observed in MD trajectories, including those induced by mutations. In this study, we model the trajectories resulting from MD simulations of the SARS-CoV-2 spike protein-ACE2, specifically the receptor-binding domain (RBD), as interresidue distance maps, and use deep convolutional neural networks to predict the functional impact of point mutations, related to the virus's infectivity and immunogenicity. Our model was successful in predicting mutant types that increase the affinity of the S protein for human receptors and reduce its immunogenicity, both based on MD trajectories (precision = 0.718; recall = 0.800; [Formula: see text] = 0.757; MCC = 0.488; AUC = 0.800) and their centroids. In an additional analysis, we also obtained a strong positive Pearson's correlation coefficient equal to 0.776, indicating a significant relationship between the average sigmoid probability for the MD trajectories and binding free energy (BFE) changes. Furthermore, we obtained a coefficient of determination of 0.602. Our 2D-RMSD analysis also corroborated predictions for more infectious and immune-evading mutants and revealed fluctuating regions within the receptor-binding motif (RBM), especially in the [Formula: see text] loop. This region presented a significant standard deviation for mutations that enable SARS-CoV-2 to evade the immune response, with RMSD values of 5Å in the simulation. This methodology offers an efficient alternative to identify potential strains of SARS-CoV-2, which may be potentially linked to more infectious and immune-evading mutations. Using clustering and deep learning techniques, our approach leverages information from the ensemble of MD trajectories to recognize a broad spectrum of multiple conformational patterns characteristic of mutant types. This represents a strategic advantage in identifying emerging variants, bypassing the need for long MD simulations. Furthermore, the present work tends to contribute substantially to the field of computational biology and virology, particularly to accelerate the design and optimization of new therapeutic agents and vaccines, offering a proactive stance against the constantly evolving threat of COVID-19 and potential future pandemics.

Antibody evasion associated with the RBD significant mutations in several emerging SARS-CoV-2 variants and its subvariants

Since the origin of the wild strain of SARS-CoV-2, several variants have emerged, which were designated as VOC, VOI, and VUM from time to time. The Omicron variant is noted as the recent VOC. After the origin of the Omicron variant on November 2021, several subvariants of Omicron have originated subsequently, like BA.1/2, BA.2.75/2.75.2, BA.4/5, BF.7, BQ.1/1.1, XBB.1/1.5, etc. which are circulated throughout the globe. Scientists reported that antibody escape is a common phenomenon observed in all the previous VOCs, VOIs, including Omicron and its subvariants. The mutations in the NTD (N-terminal domain) and RBD (Receptor-binding domain) of the spike of these variants and subvariants are responsible for antibody escape. At the same time, it has been noted that spike RBD mutations have been increasing in the last few months. This review illustrates significant RBD mutations namely R346T, K417N/T, L452R, N460K E484A/K/Q, and N501Y found in the previous emerging SARS-CoV-2 variants, including Omicron and its subvariants in high frequency and their role in antibody evasion and immune evasion. The review also describes the different classes of nAb responsible for antibody escape in SARS-CoV-2 variants and the molecular perspective of the mutation in nAb escape. It will help the future researchers to develop efficient vaccines which can finally prevent the pandemic.

Dynamics of SARS-CoV-2 Seroprevalence in a Large US population Over a Period of 12 Months

Due to a combination of asymptomatic or undiagnosed infections, the proportion of the United States population infected with SARS-CoV-2 was unclear from the beginning of the pandemic. We previously established a platform to screen for SARS-CoV-2 positivity across a representative proportion of the US population, from which we reported that almost 17 million Americans were estimated to have had undocumented infections in the Spring of 2020. Since then, vaccine rollout and prevalence of different SARS-CoV-2 variants have further altered seropositivity trends within the United States population. To explore the longitudinal impacts of the pandemic and vaccine responses on seropositivity, we re-enrolled participants from our baseline study in a 6- and 12- month follow-up study to develop a longitudinal antibody profile capable of representing seropositivity within the United States during a critical period just prior to and during the initiation of vaccine rollout. Initial measurements showed that, since July 2020, seropositivity elevated within this population from 4.8% at baseline to 36.2% and 89.3% at 6 and 12 months, respectively. We also evaluated nucleocapsid seropositivity and compared to spike seropositivity to identify trends in infection versus vaccination relative to baseline. These data serve as a window into a critical timeframe within the COVID-19 pandemic response and serve as a resource that could be used in subsequent respiratory illness outbreaks.

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.

Understanding the role of conserved proline and serine residues in the SARS-CoV-2 spike cleavage sites in the virus entry, fusion, and infectivity

The spike (S) glycoprotein of the SARS-CoV-2 virus binds to the host cell receptor and promotes the virus's entry into the target host cell. This interaction is primed by host cell proteases like furin and TMPRSS2, which act at the S1/S2 and S2´ cleavage sites, respectively. Both cleavage sites have serine or proline residues flanking either the single or polybasic region and were found to be conserved in coronaviruses. Unravelling the effects of these conserved residues on the virus entry and infectivity might facilitate the development of novel therapeutics. Here, we have investigated the role of the conserved serine and proline residues in the SARS-CoV-2 spike mediated entry, fusogenicity, and viral infectivity by using the HIV-1/spike-based pseudovirus system. A conserved serine residue mutation to alanine (S2´S-A) at the S2´ cleavage site resulted in the complete loss of spike cleavage. Exogenous treatment with trypsin or overexpression of TMPRSS2 protease could not rescue the loss of spike cleavage and biological activity. The S2´S-A mutant showed no significant responses against E-64d, TMPRSS2 or other relevant inhibitors. Taken together, serine at the S2´ site in the spike protein was indispensable for spike protein cleavage and virus infectivity. Thus, novel interventions targeting the conserved serine at the S2´ cleavage site should be explored to reduce severe disease caused by SARS-CoV-2-and novel emerging variants.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03749-y.

Physicochemical effects of emerging exchanges on the spike protein's RBM of the SARS-CoV-2 Omicron subvariants BA.1-BA.5 and its influence on the biological properties and attributes developed by these subvariants

Emerging in South Africa, SARS-CoV-2 Omicron variant was marked by the expression of an exaggerated number of mutations throughout its genome and by the emergence of subvariants, whose attributes developed by them have been associated with amino acid exchanges that occur mainly in the RBM region of the spike protein. The RBM comprises a region within the RBD and is directly involved in the SARS-CoV-2 spike protein interaction with the host cell ACE2 receptor, during the infection mechanism and viral transmission. Defined as the region from aa 437 to aa 508, there are several residues in certain positions that interact directly with the human ACE-2 receptor during these processes. The occurrence of amino acid exchanges in these positions causes physicochemical alterations in the SARS-CoV-2 spike protein, which confer additional advantages and attributes to the agent. In addition, these exchanges serve as a basis for the characterization of new variants and subvariants of SARS-CoV-2. In this review, the amino acid exchanges that have occurred in the RBM of the subvariants BA.1 to BA.5 of SARS-CoV-2 that emerged from the Omicron are described. The physicochemical effects caused by them on spike protein are also described, as well as their influence on the biological properties and attributes developed by the subvariants BA.1, BA.2, BA.3, BA.4 and BA.5.

Adsorption of SARS-CoV-2 Spike (N501Y) RBD to Human Angiotensin-Converting Enzyme 2 at a Lipid/Water Interface

The receptor binding domain (RBD) of spike proteins plays a crucial role in the process of severe acute respiratory syndrome corona virus 2 (SARS-CoV-2) attachment to the human angiotensin-converting enzyme 2 (ACE2). The N501Y mutation and later mutations introduced extra positive charges on the spike RBD and resulted in higher transmissibility, likely due to stronger binding with the highly negatively charged ACE2. Consequently, many studies have been devoted to understanding the molecular mechanism of spike protein binding with the ACE2 receptor. Most of the theoretical studies, however, have been done on isolated proteins. ACE2 is a transmembrane protein; thus, it is important to understand the interaction of spike proteins with ACE2 in a lipid matrix. In this study, the adsorption of ACE2 and spike (N501Y) RBD at a lipid/water interface was studied using the heterodyne-detected vibrational sum frequency generation (HD-VSFG) technique. The technique is a non-linear optical spectroscopy which measures vibrational spectra of molecules at an interface and provides information on their structure and orientation. It is found that ACE2 is effectively adsorbed at the positively charged 1,2-dipalmitoyl-3-trimethylammonium-propane (DPTAP) lipid monolayer via electrostatic interactions. The adsorption of ACE2 at the DPTAP monolayer causes a reorganization of interfacial water (D2O) from the D-down to the D-up orientation, indicating that the originally positively charged DPTAP interface becomes negatively charged due to ACE2 adsorption. The negatively charged interface (DPTAP/ACE2) allows further adsorption of positively charged spike RBD. HD-VSFG spectra in the amide I region show differences for spike (N501Y) RBD adsorbed at D2O, DPTAP, and DPTAP/ACE2 interfaces. A red shift observed for the spectra of spike RBD/DPTAP suggests that spike RBD oligomers are formed upon contact with DPTAP lipids.

Digestive system infection by SARS‑CoV‑2: Entry mechanism, clinical symptoms and expression of major receptors (Review)

Besides causing severe acute respiratory syndrome (SARS), SARS‑coronavirus 2 (SARS‑CoV‑2) also harms the digestive system. Given the appearance of numerous cases of SARS‑CoV‑2, it has been demonstrated that SARS‑CoV‑2 is able to harm target organs such as the gastrointestinal tract, liver and pancreas, and either worsen the condition of patients with basic digestive illnesses or make their prognosis poor. According to several previously published studies, angiotensin‑converting enzyme II (ACE2) and transmembrane serine protease II (TMPRSS2) are expressed either singly or in combination in the digestive system and in other regions of the human body. In order to change the viral conformation, create a fusion hole and release viral RNA into the host cell for replication and transcription, SARS‑CoV‑2 is capable of binding to these two proteins through the spike protein on its surface. As a result, the body experiences an immune reaction and an inflammatory reaction, which may lead to nausea, diarrhea, abdominal pain and even gastrointestinal bleeding, elevated levels of liver enzymes, acute liver injury, pancreatitis and other serious lesions. In order to provide possible strategies for the clinical diagnosis and treatment of digestive system diseases during the COVID‑19 pandemic, the molecular structure of SARS‑CoV‑2 and the mechanism via which SARS‑CoV‑2 enters the human body through ACE2 and TMPRSS2 were discussed in the present review, and the clinical manifestations of SARS‑CoV‑2 infection in the digestive system were also summarized. Finally, the expression characteristics of ACE2 and TMPRSS2 in the main target organs of the digestive system were described.

The role of S477N mutation in the molecular behavior of SARS-CoV-2 spike protein: An in-silico perspective

The attachment of SARA-CoV-2 happens between ACE2 and the receptor binding domain (RBD) on the spike protein. Mutations in this domain can affect the binding affinity of the spike protein for ACE2. S477N, one of the most common mutations reported in the recent variants, is located in the RBD. Today's computational approaches in biology, especially during the SARS-CoV-2 pandemic, assist researchers in predicting a protein's behavior in contact with other proteins in more detail. In this study, we investigated the interactions of the S477N-hACE2 in silico to find the impact of this mutation on its binding affinity for ACE2 and immunity responses using dynamics simulation, protein-protein docking, and immunoinformatics methods. Our computational analysis revealed an increased binding affinity of N477 for ACE2. Four new hydrogen and hydrophobic bonds in the mutant RBD-ACE2 were formed (with S19 and Q24 of ACE2), which do not exist in the wild type. Also, the protein spike structure in this mutation was associated with an increase in stabilization and a decrease in its fluctuations at the atomic level. N477 mutation can be considered as the cause of increased escape from the immune system through MHC-II.

Insight into free energy and dynamic cross-correlations of residue for binding affinity of antibody and receptor binding domain SARS-CoV-2

SARS-CoV-2 virus continues to evolve and mutate causing most of the mutated variants resist to many of the therapeutic monoclonal antibodies (mAbs). Despite several mAbs retained neutralizing capability for Omicron BA.1 and BA.2, reduction in neutralization potency was reported. Hence, effort of searching for mAb that is broader in neutralization breadth without losing the neutralizing ability is continued. MW06 was reported with capability in neutralizing most of the variants of concern (VOC) and it binds to the conserved region (left flank) near epitope mAb sotrovimab (S309). In this study, binding affinity of mAb MW06 and its cocktail formulation with MW05 for receptor binding domain (RBD) SARS-CoV-2 virus was investigated under molecular dynamics simulations (MDs). Binding free energies computed by Molecular Mechanics Generalised Born Surface Area (MM-GBSA) algorithm predicted the binding affinity of MW06 for RBD BA.1 (-53 kcal/mol) as strong as RBD wildtype (-58 kcal/mol) while deterioration was observed for RBD BA.2 (-43 kcal/mol). Alike S309 and MW06, simulated cocktail mAb (MW05 and MW06)-RBD interactions suggested the neutralizing capability of the cocktail formulation for RBD BA.1 and BA.2 reduced. Meanwhile, residue pairs that favour the communication between the mAb and RBD have been identified by decomposing the free energy per pairwise residue basis. Apart from understanding the effects of mutation occurred in the RBD region on human angiotensin-converting enzyme 2 (hACE2) binding, impact of heavily mutated RBD on mAb-RBD interactions was investigated in this study as well. In addition to energetic profile obtained from MDs, plotting the dynamics cross-correlation map of the mAb-RBD complex under elastic network model (ENM) was aimed to understand the cross-correlations between residue fluctuations. It allows simple and rapid analysis on the motions or dynamics of the protein residues of mAbs and RBD in complex. Protein residues having correlated motions are normally part of the structural domains of the protein and their respective motions and protein function are related. Motion of mutated RBD residues and mAb residues was less correlated while their respective interactions energy computed to be higher. The combined techniques of MDs and ENM offered simplicity in understanding dynamics and energy contribution that explain binding affinity of mAb-RBD complexes.

Rapid and Sensitive Genotyping of SARS-CoV-2 Key Mutation L452R with an RPA-PfAgo Method

In the two years of COVID-19 pandemic, the SARS-CoV-2 variants have caused waves of infections one after another, and the pandemic is not ending. The key mutations on the S protein enable the variants with enhanced viral infectivity, immune evasion, and/or antibody neutralization resistance, bringing difficulties to epidemic prevention and control. In support of precise epidemic control and precision medicine of the virus, a fast and simple genotyping method for the key mutations of SARS-CoV-2 variants needs to be developed. By utilizing the specific recognition and cleavage property of the nuclease Argonaute from Pyrococcus furiosus (PfAgo), we developed a recombinase polymerase amplification (RPA) and PfAgo combined method for a rapid and sensitive genotyping of SARS-CoV-2 key mutation L452R. With a delicate design of the strategy, careful screening of the RPA primers and PfAgo gDNA, and optimization of the reaction, the method achieves a high sensitivity of a single copy per reaction, which is validated with the pseudovirus. This is the highest sensitivity that can be achieved theoretically and the highest sensitivity as compared to the available SARS-CoV-2 genotyping assays. Using RPA, the procedure of the method is finished within 1.5 h and only needs a minimum laboratorial support, suggesting that the method can be easily applied locally or on-site. The RPA-PfAgo method established in this study provides a strong support to the precise epidemic control and precision medicine of SARS-CoV-2 variants and can be readily developed for the simultaneous genotyping of multiple SARS-CoV-2 mutations.

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).

Identification of Doxorubicin as Repurposing Inhibitory Drug for MERS-CoV PLpro

Middle East respiratory syndrome coronavirus (MERS-CoV), belonging to the betacoronavirus genus can cause severe respiratory illnesses, accompanied by pneumonia, multiorgan failure, and ultimately death. CoVs have the ability to transgress species barriers and spread swiftly into new host species, with human-to-human transmission causing epidemic diseases. Despite the severe public health threat of MERS-CoV, there are currently no vaccines or drugs available for its treatment. MERS-CoV papain-like protease (PLpro) is a key enzyme that plays an important role in its replication. In the present study, we evaluated the inhibitory activities of doxorubicin (DOX) against the recombinant MERS-CoV PLpro by employing protease inhibition assays. Hydrolysis of fluorogenic peptide from the Z-RLRGG-AMC-peptide bond in the presence of DOX showed an IC50 value of 1.67 μM at 30 min. Subsequently, we confirmed the interaction between DOX and MERS-CoV PLpro by thermal shift assay (TSA), and DOX increased ΔTm by ~20 °C, clearly indicating a coherent interaction between the MERS-CoV PL protease and DOX. The binding site of DOX on MERS-CoV PLpro was assessed using docking techniques and molecular dynamic (MD) simulations. DOX bound to the thumb region of the catalytic domain of the MERS-CoV PLpro. MD simulation results showed flexible BL2 loops, as well as other potential residues, such as R231, R233, and G276 of MERS-CoV PLpro. Development of drug repurposing is a remarkable opportunity to quickly examine the efficacy of different aspects of treating various diseases. Protease inhibitors have been found to be effective against MERS-CoV to date, and numerous candidates are currently undergoing clinical trials to prove this. Our effort follows a in similar direction.

IL-2 Flexible Loops Might Play a Role in IL-2 Interaction with the High-Affinity IL-2 Receptor: A Molecular Dynamics (MD) Study

The clinical use of high-dose IL-2 in cancer immunotherapy faces several drawbacks such as toxicity and unfavorable pharmacokinetic profile. These drawbacks can be avoided by inhibiting IL-2 interaction with the CD25 subunit, which is a component of the high-affinity IL-2 receptor (IL-2Rαβγ). Several studies showed mutations of potential IL-2 residues such as R38, F42, Y45, and Y72 would produce IL-2 that is CD25-independent. In essence, structural comparison between wild-type (WT) IL-2 and CD25-independent IL-2 can be very insightful to assess the role of IL-2 flexibility and conformation in the IL-2 receptor interactions. Here, we investigated the flexibility loops and conformation of IL-2m (F24A, Y45A, and L72G), which is known to be CD25-independent, and IL-2m2 (F42Y and L72R) mutants along with WT IL-2 using MD simulations. Despite residue mutations, both IL-2m and IL-2m2 showed comparable conformational compactness and better stability than WT IL-2. Interestingly, IL-2m and IL-2m2 mutants showed rigid BC and CD loops in comparison to WT IL-2 . Also, the AB loop conformation of IL-2m was a bent structure compared to the WT IL-2 and IL-2m2. Principal component analysis (PCA) and free-energy landscape results suggested IL-2m and IL-2m2 have stable conformations compared to the WT IL-2. Therefore, these mutation sites of IL-2 produced stable and rigid loops that might prevent IL-2 from binding to the CD25 subunit. Our results can help to assess IL-2 flexibility loops to design new CD25-independent IL-2 mutants without compromising the IL-2 structure.

Evolutionary trajectory of receptor binding specificity and promiscuity of the spike protein of SARS-CoV-2

SARS-CoV-2 infects cells by attachment to its receptor-the angiotensin converting enzyme 2 (ACE2). Regardless of the wealth of structural data, little is known about the physicochemical mechanism of interactions of the viral spike (S) protein with ACE2 and how this mechanism has evolved during the pandemic. Here, we applied experimental and computational approaches to characterize the molecular interaction of S proteins from SARS-CoV-2 variants of concern (VOC). Data on kinetics, activation-, and equilibrium thermodynamics of binding of the receptor binding domain (RBD) from VOC with ACE2 as well as data from computational protein electrostatics revealed a profound remodeling of the physicochemical characteristics of the interaction during the evolution. Thus, as compared to RBDs from Wuhan strain and other VOC, Omicron RBD presented as a unique protein in terms of conformational dynamics and types of non-covalent forces driving the complex formation with ACE2. Viral evolution resulted in a restriction of the RBD structural dynamics, and a shift to a major role of polar forces for ACE2 binding. Further, we investigated how the reshaping of the physicochemical characteristics of interaction affects the binding specificity of S proteins. Data from various binding assays revealed that SARS-CoV-2 Wuhan and Omicron RBDs manifest capacity for promiscuous recognition of unrelated human proteins, but they harbor distinct reactivity patterns. These findings might contribute for mechanistic understanding of the viral tropism and capacity to evade immune responses during evolution.

When the Dust Has Settled: Calculation of Binding Affinities from First Principles for SARS-CoV-2 Variants with Quantitative Accuracy

Accurate determination of binding free energy is pivotal for the study of many biological processes and has been applied in a number of theoretical investigations to compare the affinity of severe acute respiratory syndrome coronavirus 2 variants toward the host cell. Diversity of these variants challenges the development of effective general therapies, their transmissibility relying either on an increased affinity toward their dedicated human receptor, the angiotensin-converting enzyme 2 (ACE2), or on escaping the immune response. Now that robust structural data are available, we have determined with utmost accuracy the standard binding free energy of the receptor-binding domain to the most widespread variants, namely, Alpha, Beta, Delta, and Omicron BA.2, as well as the wild type (WT) in complex either with ACE2 or with antibodies, namely, S2E12 and H11-D4, using a rigorous theoretical framework that combines molecular dynamics and potential-of-mean-force calculations. Our results show that an appropriate starting structure is crucial to ensure appropriate reproduction of the binding affinity, allowing the variants to be compared. They also emphasize the necessity to apply the relevant methodology, bereft of any shortcut, to account for all the contributions to the standard binding free energy. Our estimates of the binding affinities support the view that while the Alpha and Beta variants lean on an increased affinity toward the host cell, the Delta and Omicron BA.2 variants choose immune escape. Moreover, the S2E12 antibody, already known to be active against the WT (Starr et al., 2021; Mlcochova et al., 2021), proved to be equally effective against the Delta variant. In stark contrast, H11-D4 retains a low affinity toward the WT compared to that of ACE2 for the latter. Assuming robust structural information, the methodology employed herein successfully addresses the challenging protein-protein binding problem in the context of coronavirus disease 2019 while offering promising perspectives for predictive studies of ever-emerging variants.

Mutation informatics: SARS-CoV-2 receptor-binding domain of the spike protein

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) undergoes mutations at a high rate and with frequent genetic reassortment (antigenic drift/shift), leading to variability in targets. The receptor-binding domain (RBD) of the spike (S) protein has a major role in the binding of SARS-CoV-2 with human angiotensin-converting enzyme 2 (ACE2). Mutations at the RBD influence the binding interaction at the SARS-CoV-2 S-ACE2 interface and impact viral pathogenicity. Here, we discuss different reported mutations of concern in RBD, physicochemical characteristic changes resulting from mutated amino acids and their effect on binding between the RBD and ACE2. Along with mutation informatics, we highlight recently developed small-molecule inhibitors of RBD and the ACE2 interface. This information provides a rational basis for the design of inhibitors against the multivariant strains of SARS-CoV-2.

Spike protein of SARS-CoV-2 variants: a brief review and practical implications

The scientific community has been alarmed by the possible immunological evasion, higher infectivity, and severity of disease caused by the newest variants of SARS-CoV-2. The spike protein has an important role in the cellular invasion of viruses and is the target of several vaccines and therapeutic resources, such as monoclonal antibodies. In addition, some of the most relevant mutations in the different variants are on the spike (S) protein gene sequence that leads to structural alterations in the predicted protein, thus causing concern about the protection mediated by vaccines against these new strains. The present review highlights the most recent knowledge about COVID-19 and vaccines, emphasizing the different spike protein structures of SARS-CoV-2 and updating the reader about the emerging viral variants and their classifications, the more common viral mutations described and their distribution in Brazil. It also compiles a table with the most recent knowledge about all of the Omicron spike mutations.

Impact of SARS-CoV-2 RBD Mutations on the Production of a Recombinant RBD Fusion Protein in Mammalian Cells

SARS-CoV-2 receptor-binding domain (RBD) is a major target for the development of diagnostics, vaccines and therapeutics directed against COVID-19. Important efforts have been dedicated to the rapid and efficient production of recombinant RBD proteins for clinical and diagnostic applications. One of the main challenges is the ongoing emergence of SARS-CoV-2 variants that carry mutations within the RBD, resulting in the constant need to design and optimise the production of new recombinant protein variants. We describe here the impact of naturally occurring RBD mutations on the secretion of a recombinant Fc-tagged RBD protein expressed in HEK 293 cells. We show that mutation E484K of the B.1.351 variant interferes with the proper disulphide bond formation and folding of the recombinant protein, resulting in its retention into the endoplasmic reticulum (ER) and reduced protein secretion. Accumulation of the recombinant B.1.351 RBD-Fc fusion protein in the ER correlated with the upregulation of endogenous ER chaperones, suggestive of the unfolded protein response (UPR). Overexpression of the chaperone and protein disulphide isomerase PDIA2 further impaired protein secretion by altering disulphide bond formation and increasing ER retention. This work contributes to a better understanding of the challenges faced in producing mutant RBD proteins and can assist in the design of optimisation protocols.

Understanding gilteritinib resistance to FLT3-F691L mutation through an integrated computational strategy

FMS-like tyrosine kinase 3 (FLT3) serves as an important drug target for acute myeloid leukemia (AML), and gene mutations of FLT3 have been closely associated with AML patients with an incidence rate of ~ 30%. However, the mechanism of the clinically relevant F691L gatekeeper mutation conferred resistance to the drug gilteritinib remained poorly understood. In this study, multiple microsecond molecular dynamics (MD) simulations, end-point free energy calculations, and dynamic correlated and network analyses were performed to investigate the molecular basis of gilteritinib resistance to the FLT3-F691L mutation. The simulations revealed that the resistant mutation largely induced the conformational changes of the activation loop (A-loop), the phosphate-binding loop, and the helix αC of the FLT3 protein. The binding abilities of the gilteritinib to the wild-type and the F691L mutant were different through the binding free energy prediction. The simulation results further indicated that the driving force to determine the binding affinity of gilteritinib was derived from the differences in the energy terms of electrostatic and van der Waals interactions. Moreover, the per-residue free energy decomposition suggested that the four residues (Phe803, Gly831, Leu832, and Ala833) located at the A-loop of FLT3 had a significant impact on the binding affinity of gilteritinib to the F691L mutant. This study may provide useful information for the design of novel FLT3 inhibitors specially targeting the F691L gatekeeper mutant.

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.

Physicochemical effect of the N501Y, E484K/Q, K417N/T, L452R and T478K mutations on the SARS-CoV-2 spike protein RBD and its influence on agent fitness and on attributes developed by emerging variants of concern

The spike protein comprises one of the main structural components of SARS-CoV-2 because it is directly involved in the infection process and viral transmission, and also because of its immunogenic properties, as an inducer of the protective antibodies production and as a vaccine component. The occurrence of mutations in this region or in other the virus genome regions, comprises a natural phenomenon in its evolution. However, they also occur due to the selective immune pressure, to which the agent is continuously subjected, especially in the spike protein immunodominant regions, such as the RBD. Mutations in the spike protein can change the virus' fitness, increasing its affinity for target cells, its transmissibility and its virulence. In addition, these mutations can giving it the potential ability to evade the protective antibodies action obtained from convalescent sera or vaccine origin, as well as those used in therapy, which may favor the virus expansion and compromise the infection control. Five mutations N501Y, E484K/Q, K417N/T, L452R and T478K, located in the spike protein RBD, have had a greater impact because they are associated with new attributes developed by the virus, which characterize the emerging variants of concern (VOCs) of SARS-Cov-2 identified so far. The occurrence of these mutations induces complex physicochemical effects that can alter the spike protein's structure and its function, which in turn, lead to changes in the agents' fitness. This manuscript discusses the attributes of VOCs associated with the physicochemical effects caused by the aforementioned mutations.

The Glu143 Residue Might Play a Significant Role in T20 Peptide Binding to HIV-1 Receptor gp41: An In Silico Study

Despite the enormous efforts made to develop other fusion inhibitors for HIV, the enfuvirtide (known as T20) peptide is the only approved HIV-1 inhibitory drug so far. Investigating the role of potential residues of the T20 peptide’s conformational dynamics could help us to understand the role of potential residues of the T20 peptide. We investigated T20 peptide conformation and binding interactions with the HIV-1 receptor (i.e., gp41) using MD simulations and docking techniques, respectively. Although the mutation of E143 into alanine decreased the flexibility of the E143A mutant, the conformational compactness of the mutant was increased. This suggests a potential role of E143 in the T20 peptide’s conformation. Interestingly, the free energy landscape showed a significant change in the wild-type T20 minimum, as the E143A mutant produced two observed minima. Finally, the docking results of T20 to the gp41 receptor showed a different binding interaction in comparison to the E143A mutant. This suggests that E143 residue can influence the binding interaction with the gp41 receptor. Overall, the E143 residue showed a significant role in conformation and binding to the HIV-1 receptor. These findings can be helpful in optimizing and developing HIV-1 inhibitor peptides.

Insights into the structure and dynamics of SARS-CoV-2 spike glycoprotein double mutant L452R-E484Q

The Receptor Binding Domain (RBD) of SARS-CoV-2, located on the S1 subunit, plays a vital role in the virus binding and its entry into the host cell through angiotensin-converting enzyme 2 (ACE2) receptor. Therefore, understanding the dynamic effects of mutants on the SARS-CoV-2 RBD is essential for discovering drugs to inhibit the virus binding and disrupt its entry into the host cells. A recent study reported a double mutant of SARS-CoV-2, L452R-E484Q, located in the RBD region. Thus, this study employed various computational algorithms and methods to understand the structural impact of both individual variants L452R, E484Q, and the double mutant L452R-E484Q on the native RBD of spike glycoprotein. The effects of the mutations on native RBD structure were predicted by in silico algorithms, which predicted changes in the protein structure and function upon the mutations. Subsequently, molecular dynamics (MD) simulations were employed to understand the conformational stability and functional changes on the RBD upon the mutations. The comparative results of MD simulation parameters displayed that the double mutant induces significant conformational changes in the spike glycoprotein RBD, which may alter its biological functions.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03151-0.

Host adaptation of codon usage in SARS-CoV-2 from mammals indicates potential natural selection and viral fitness

SARS-CoV-2 infection, which is the cause of the COVID-19 pandemic, has expanded across various animal hosts, and the virus can be transmitted particularly efficiently in minks. It is still not clear how SARS-CoV-2 is selected and evolves in its hosts, or how mutations affect viral fitness. In this report, sequences of SARS-CoV-2 isolated from human and animal hosts were analyzed, and the binding energy and capacity of the spike protein to bind human ACE2 and the mink receptor were compared. Codon adaptation index (CAI) analysis indicated the optimization of viral codons in some animals such as bats and minks, and a neutrality plot demonstrated that natural selection had a greater influence on some SARS-CoV-2 sequences than mutational pressure. Molecular dynamics simulation results showed that the mutations Y453F and N501T in mink SARS-CoV-2 could enhance the binding of the viral spike to the mink receptor, indicating the involvement of these mutations in natural selection and viral fitness. Receptor binding analysis revealed that the mink SARS-CoV-2 spike interacted more strongly with the mink receptor than the human receptor. Tracking the variations and codon bias of SARS-CoV-2 is helpful for understanding the fitness of the virus in virus transmission, pathogenesis, and immune evasion.