Exploring the Association Between Sialic Acid and SARS-CoV-2 Spike Protein Through a Molecular Dynamics-Based Approach

Multi-combined QSAR, molecular docking, molecular dynamics simulation, and ADMET of Flavonoid derivatives as potent cholinesterase inhibitors

In searching for a new and efficient therapeutic agent against Alzheimer's disease, a Quantitative structure-activity relationship (QSAR) was derived for 45 Flavonoid derivatives recently synthesized and evaluated as cholinesterase inhibitors. The multiple linear regression method (MLR) was adopted to develop an adequate mathematical model that describes the relationship between a variety of molecular descriptors of the studied compounds and their biological activities (cholinesterase inhibitors). Golbraikh and Tropsha criteria were applied to verify the validity of the built model. The built MLR model was statistically reliable, robust, and predictive (R2 = 0.801, Q2cv = 0.876, R2test = 0.824). Dreiding energy and Molar Refractivity were the major factors that govern the Anti-cholinesterase activity. These results were further exploited to design a new series of Flavonoid derivatives with higher Anti-cholinesterase activities than the existing ones. Thereafter, molecular docking and molecular dynamic studies were performed to predict the binding types of the designed compounds and to investigate their stability at the active site of the Butyrylcholinestérase BuChE protein. The negative and low binding affinity calculated for all designed compounds shows that designed compound 1 has a favorable affinity for the 4TPK. Moreover, molecular dynamics simulation studies confirmed the stability of designed compound 1 in the active pocket of 4TPK over 100 ns. Finally, the ADMET analysis was incorporated to analyze the pharmacokinetics and toxicity parameters. The designed compounds were found to meet the ADMET descriptor criteria at an acceptable level having respectable intestinal permeability and water solubility and can reach the intended destinations.Communicated by Ramaswamy H. Sarma.

Examining sialic acid derivatives as potential inhibitors of SARS-CoV-2 spike protein receptor binding domain

Severe acute respiratory syndrome coronavirus 2 (SARS CoV-2) has been the primary reason behind the COVID-19 global pandemic which has affected millions of lives worldwide. The fundamental cause of the infection is the molecular binding of the viral spike protein receptor binding domain (SP-RBD) with the human cell angiotensin-converting enzyme 2 (ACE2) receptor. The infection can be prevented if the binding of RBD-ACE2 is resisted by utilizing certain inhibitors or drugs that demonstrate strong binding affinity towards the SP RBD. Sialic acid based glycans found widely in human cells and tissues have notable propensity of binding to viral proteins of the coronaviridae family. Recent experimental literature have used N-acetyl neuraminic acid (Sialic acid) to create diagnostic sensors for SARS-CoV-2, but a detailed interrogation of the underlying molecular mechanisms is warranted. Here, we perform all atom molecular dynamics (MD) simulations for the complexes of certain Sialic acid-based molecules with that of SP RBD of SARS CoV-2. Our results indicate that Sialic acid not only reproduces a binding affinity comparable to the RBD-ACE2 interactions, it also assumes the longest time to dissociate completely from the protein binding pocket of SP RBD. Our predictions corroborate that a combination of electrostatic and van der Waals energies as well the polar hydrogen bond interactions between the RBD residues and the inhibitors influence free energy of binding.Communicated by Ramaswamy H. Sarma.

Single-molecule imaging reveals allosteric stimulation of SARS-CoV-2 spike receptor binding domain by host sialic acid

Conformational dynamics of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike glycoprotein (S) mediate exposure of the binding site for the cellular receptor, angiotensin-converting enzyme 2 (ACE2). The N-terminal domain (NTD) of S binds terminal sialic acid (SA) moieties on the cell surface, but the functional role of this interaction in virus entry is unknown. Here, we report that NTD-SA interaction enhances both S-mediated virus attachment and ACE2 binding. Through single-molecule Förster resonance energy transfer imaging of individual S trimers, we demonstrate that SA binding to the NTD allosterically shifts the S conformational equilibrium, favoring enhanced exposure of the ACE2-binding site. Antibodies that target the NTD block SA binding, which contributes to their mechanism of neutralization. These findings inform on mechanisms of S activation at the cell surface.

Two Receptor Binding Strategy of SARS-CoV-2 Is Mediated by Both the N-Terminal and Receptor-Binding Spike Domain

It is not well understood why severe acute respiratory syndrome (SARS)-CoV-2 spreads much faster than other β-coronaviruses such as SARS-CoV and Middle East respiratory syndrome (MERS)-CoV. In a previous publication, we predicted the binding of the N-terminal domain (NTD) of SARS-CoV-2 spike to sialic acids (SAs). Here, we experimentally validate this interaction and present simulations that reveal a second possible interaction between SAs and the spike protein via a binding site located in the receptor-binding domain (RBD). The predictions from molecular-dynamics simulations and the previously-published 2D-Zernike binding-site recognition approach were validated through flow-induced dispersion analysis (FIDA)─which reveals the capability of the SARS-CoV-2 spike to bind to SA-containing (glyco)lipid vesicles, and flow-cytometry measurements─which show that spike binding is strongly decreased upon inhibition of SA expression on the membranes of angiotensin converting enzyme-2 (ACE2)-expressing HEK cells. Our analyses reveal that the SA binding of the NTD and RBD strongly enhances the infection-inducing ACE2 binding. Altogether, our work provides in silico, in vitro, and cellular evidence that the SARS-CoV-2 virus utilizes a two-receptor (SA and ACE2) strategy. This allows the SARS-CoV-2 spike to use SA moieties on the cell membrane as a binding anchor, which increases the residence time of the virus on the cell surface and aids in the binding of the main receptor, ACE2, via 2D diffusion.

Triterpene Derivatives as Potential Inhibitors of the RBD Spike Protein from SARS-CoV-2: An In Silico Approach

The appearance of a new coronavirus, SARS-CoV-2, in 2019 kicked off an international public health emergency. Although rapid progress in vaccination has reduced the number of deaths, the development of alternative treatments to overcome the disease is still necessary. It is known that the infection begins with the interaction of the spike glycoprotein (at the virus surface) and the angiotensin-converting enzyme 2 cell receptor (ACE2). Therefore, a straightforward solution for promoting virus inhibition seems to be the search for molecules capable of abolishing such attachment. In this work, we tested 18 triterpene derivatives as potential inhibitors of SARS-CoV-2 against the receptor-binding domain (RBD) of the spike protein by means of molecular docking and molecular dynamics simulations, modeling the RBD S1 subunit from the X-ray structure of the RBD-ACE2 complex (PDB ID: 6M0J). Molecular docking revealed that at least three triterpene derivatives of each type (i.e., oleanolic, moronic and ursolic) present similar interaction energies as the reference molecule, i.e., glycyrrhizic acid. Molecular dynamics suggest that two compounds from oleanolic and ursolic acid, OA5 and UA2, can induce conformational changes capable of disrupting the RBD-ACE2 interaction. Finally, physicochemical and pharmacokinetic properties simulations revealed favorable biological activity as antivirals.

Involvement of sialoglycans in SARS-COV-2 infection: Opportunities and challenges for glyco-based inhibitors

Viral infections have been the causes of global pandemics, including the ongoing coronavirus disease 2019, which prompted the investigation into the infection mechanisms to find treatment and aid the vaccine design. Betacoronaviruses use spike glycoprotein on their surface to bind to host receptors, aiding their host attachment and cell fusion. Protein-glycan interaction has been implicated in the viral entry mechanism of many viruses and has recently been shown in SARS-CoV-2. Here, we reviewed the current knowledge on protein-glycan interactions that facilitate SARS-CoV-2 host entry, with special interest in sialoglycans present on both the virions and host cell surfaces. We also analyze how such information provides opportunities and challenges in glyco-based inhibitors.

Cell Entry and Unusual Replication of SARS-CoV-2

BACKGROUND

SARS-CoV-2 is the causative virus for the CoVID-19 pandemic that has frequently mutated to continue to infect and resist available vaccines. Emerging new variants of the virus have complicated notions of immunity conferred by vaccines versus immunity that results from infection. While we continue to progress from epidemic to endemic as a result of this collective immunity, the pandemic remains a morbid and mortal problem.

OBJECTIVE

The SARS-CoV-2 virus has a very complex manner of replication. The spike protein, one of the four structural proteins of the encapsulated virus, is central to the ability of the virus to penetrate cells to replicate. The objective of this review is to summarize these complex features of viral replication.

METHODS

A review of the recent literature was performed on the biology of SARS-CoV-2 infection from published work from PubMed and works reported to preprint servers, e.g., bioRxiv and medRxiv.

RESULTS AND CONCLUSION

The complex molecular and cellular biology involved in SARS-CoV-2 replication and the origination of >30 proteins from a single open reading frame (ORF) have been summarized, as well as the structural biology of spike protein, a critical factor in the cellular entry of the virus, which is a necessary feature for it to replicate and cause disease.

Molecular Docking and Dynamics Simulation Revealed the Potential Inhibitory Activity of New Drugs against Human Topoisomerase I Receptor

Human Topoisomerase I (hTop1p) is a ubiquitous enzyme that relaxes supercoiled DNA through a conserved mechanism involving transient breakage, rotation, and binding. Htop1p is the molecular target of the chemotherapeutic drug camptothecin (CPT). It causes the hTop1p-DNA complex to slow down the binding process and clash with the replicative machinery during the S phase of the cell cycle, forcing cells to activate the apoptotic response. This gives hTop1p a central role in cancer therapy. Recently, two artesunic acid derivatives (compounds c6 and c7) have been proposed as promising inhibitors of hTop1p with possible antitumor activity. We used several computational approaches to obtain in silico confirmations of the experimental data and to form a comprehensive dynamic description of the ligand-receptor system. We performed molecular docking analyses to verify the ability of the two new derivatives to access the enzyme-DNA interface, and a classical molecular dynamics simulation was performed to assess the capacity of the two compounds to maintain a stable binding pose over time. Finally, we calculated the noncovalent interactions between the two new derivatives and the hTop1p receptor in order to propose a possible inhibitory mechanism like that adopted by CPT.

Structural and energetic analyses of SARS-CoV-2 N-terminal domain characterise sugar binding pockets and suggest putative impacts of variants on COVID-19 transmission

Coronavirus disease 2019 (COVID-19) caused by SARS-CoV-2 is an ongoing pandemic that causes significant health/socioeconomic burden. Variants of concern (VOCs) have emerged affecting transmissibility, disease severity and re-infection risk. Studies suggest that the - N-terminal domain (NTD) of the spike protein may have a role in facilitating virus entry via sialic-acid receptor binding. Furthermore, most VOCs include novel NTD variants. Despite global sequence and structure similarity, most sialic-acid binding pockets in NTD vary across coronaviruses. Our work suggests ongoing evolutionary tuning of the sugar-binding pockets and recent analyses have shown that NTD insertions in VOCs tend to lie close to loops. We extended the structural characterisation of these sugar-binding pockets and explored whether variants could enhance sialic acid-binding. We found that recent NTD insertions in VOCs (i.e., Gamma, Delta and Omicron variants) and emerging variants of interest (VOIs) (i.e., Iota, Lambda and Theta variants) frequently lie close to sugar-binding pockets. For some variants, including the recent Omicron VOC, we find increases in predicted sialic acid-binding energy, compared to the original SARS-CoV-2, which may contribute to increased transmission. These binding observations are supported by molecular dynamics simulations (MD). We examined the similarity of NTD across Betacoronaviruses to determine whether the sugar-binding pockets are sufficiently similar to be exploited in drug design. Whilst most pockets are too structurally variable, we detected a previously unknown highly structurally conserved pocket which can be investigated in pursuit of a generic pan-Betacoronavirus drug. Our structure-based analyses help rationalise the effects of VOCs and provide hypotheses for experiments. Our findings suggest a strong need for experimental monitoring of changes in NTD of VOCs.

Electrostatic Features for the Receptor Binding Domain of SARS-COV-2 Wildtype and Its Variants. Compass to the Severity of the Future Variants with the Charge-Rule

Electrostatic intermolecular interactions are important in many aspects of biology. We have studied the main electrostatic features involved in the interaction of the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein with the human receptor Angiotensin-converting enzyme 2 (ACE2). As the principal computational tool, we have used the FORTE approach, capable to model proton fluctuations and computing free energies for a very large number of protein-protein systems under different physical-chemical conditions, here focusing on the RBD-ACE2 interactions. Both the wild-type and all critical variants are included in this study. From our large ensemble of extensive simulations, we obtain, as a function of pH, the binding affinities, charges of the proteins, their charge regulation capacities, and their dipole moments. In addition, we have calculated the pK_as for all ionizable residues and mapped the electrostatic coupling between them. We are able to present a simple predictor for the RBD-ACE2 binding based on the data obtained for Alpha, Beta, Gamma, Delta, and Omicron variants, as a linear correlation between the total charge of the RBD and the corresponding binding affinity. This "RBD charge rule" should work as a quick test of the degree of severity of the coming SARS-CoV-2 variants in the future.

Lactoferrin Inhibition of the Complex Formation between ACE2 Receptor and SARS CoV-2 Recognition Binding Domain

The present investigation focuses on the analysis of the interactions among human lactoferrin (LF), SARS-CoV-2 receptor-binding domain (RBD) and human angiotensin-converting enzyme 2 (ACE2) receptor in order to assess possible mutual interactions that could provide a molecular basis of the reported preventative effect of lactoferrin against CoV-2 infection. In particular, kinetic and thermodynamic parameters for the pairwise interactions among the three proteins were measured via two independent techniques, biolayer interferometry and latex nanoparticle-enhanced turbidimetry. The results obtained clearly indicate that LF is able to bind the ACE2 receptor ectodomain with significantly high affinity, whereas no binding to the RBD was observed up to the maximum “physiological” lactoferrin concentration range. Lactoferrin, above 1 µM concentration, thus appears to directly interfere with RBD–ACE2 binding, bringing about a measurable, up to 300-fold increase of the K_D value relative to RBD–ACE2 complex formation.

The SARS-CoV-2 Spike Glycoprotein Directly Binds Exogeneous Sialic Acids: A NMR View

The interaction of the SARS CoV2 spike glycoprotein with two sialic acid-containing trisaccharides (α2,3 and α2,6 sialyl N-acetyllactosamine) has been demonstrated by NMR. The NMR-based distinction between the signals of those sialic acids in the glycans covalently attached to the spike protein and those belonging to the exogenous α2,3 and α2,6 sialyl N-acetyllactosamine ligands has been achieved by synthesizing uniformly 13 C-labelled trisaccharides at the sialic acid and galactose moieties. STD-1 H,13 C-HSQC NMR experiments elegantly demonstrate the direct interaction of the sialic acid residues of both trisaccharides with additional participation of the galactose moieties, especially for the α2,3-linked analogue. Additional experiments with the spike protein in the presence of a specific antibody for the N-terminal domain and with the isolated receptor binding and N-terminal domains of the spike protein unambiguously show that the sialic acid binding site is located at the N-terminal domain.

Mucins Inhibit Coronavirus Infection in a Glycan-Dependent Manner

Mucins are a diverse and heterogeneous family of glycoproteins that comprise the bulk of mucus and the epithelial glycocalyx. Mucins are intimately involved in viral transmission. Mucin and virus laden particles can be expelled from the mouth and nose to later infect others. Viruses must also penetrate the mucus layer before cell entry and replication. The role of mucins and their molecular structure have not been well-characterized in coronavirus transmission studies. Laboratory studies predicting high rates of fomite transmission have not translated to real-world infections, and mucins may be one culprit. Here, we probed both surface and direct contact transmission scenarios for their dependence on mucins and their structure. We utilized disease-causing, bovine-derived, human coronavirus OC43. We found that bovine mucins could inhibit the infection of live cells in a concentration- and glycan-dependent manner. The effects were observed in both mock fomite and direct contact transmission experiments and were not dependent upon surface material or time-on-surface. However, the effects were abrogated by removal of the glycans or in a cross-species infection scenario where bovine mucin could not inhibit the infection of a murine coronavirus. Together, our data indicate that the mucin molecular structure plays a complex and important role in host defense.

A Computational Approach to Investigate TDP-43 RNA-Recognition Motif 2 C-Terminal Fragments Aggregation in Amyotrophic Lateral Sclerosis

Many of the molecular mechanisms underlying the pathological aggregation of proteins observed in neurodegenerative diseases are still not fully understood. Among the aggregate-associated diseases, Amyotrophic Lateral Sclerosis (ALS) is of relevant importance. In fact, although understanding the processes that cause the disease is still an open challenge, its relationship with protein aggregation is widely known. In particular, human TDP-43, an RNA/DNA binding protein, is a major component of the pathological cytoplasmic inclusions observed in ALS patients. Indeed, the deposition of the phosphorylated full-length TDP-43 in spinal cord cells has been widely studied. Moreover, it has also been shown that the brain cortex presents an accumulation of phosphorylated C-terminal fragments (CTFs). Even if it is debated whether the aggregation of CTFs represents a primary cause of ALS, it is a hallmark of TDP-43 related neurodegeneration in the brain. Here, we investigate the CTFs aggregation process, providing a computational model of interaction based on the evaluation of shape complementarity at the molecular interfaces. To this end, extensive Molecular Dynamics (MD) simulations were conducted for different types of protein fragments, with the aim of exploring the equilibrium conformations. Adopting a newly developed approach based on Zernike polynomials, able to find complementary regions in the molecular surface, we sampled a large set of solvent-exposed portions of CTFs structures as obtained from MD simulations. Our analysis proposes and assesses a set of possible association mechanisms between the CTFs, which could drive the aggregation process of the CTFs. To further evaluate the structural details of such associations, we perform molecular docking and additional MD simulations to propose possible complexes and assess their stability, focusing on complexes whose interacting regions are both characterized by a high shape complementarity and involve β3 and β5 strands at their interfaces.

Inhibition of Severe Acute Respiratory Syndrome Coronavirus 2 Replication by Hypertonic Saline Solution in Lung and Kidney Epithelial Cells

An unprecedented global health crisis has been caused by a new virus called severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). We performed experiments to test if a hypertonic saline solution was capable of inhibiting virus replication. Our data show that 1.2% NaCl inhibited virus replication by 90%, achieving 100% of inhibition at 1.5% in the nonhuman primate kidney cell line Vero, and 1.1% of NaCl was sufficient to inhibit the virus replication by 88% in human epithelial lung cell line Calu-3. Furthermore, our results indicate that the inhibition is due to an intracellular mechanism and not to the dissociation of the spike SARS-CoV-2 protein and its human receptor. NaCl depolarizes the plasma membrane causing a low energy state (high ADP/ATP concentration ratio) without impairing mitochondrial function, supposedly associated with the inhibition of the SARS-CoV-2 life cycle. Membrane depolarization and intracellular energy deprivation are possible mechanisms by which the hypertonic saline solution efficiently prevents virus replication in vitro assays.

Up State of the SARS-COV-2 Spike Homotrimer Favors an Increased Virulence for New Variants

The COVID-19 pandemic has spread worldwide. However, as soon as the first vaccines-the only scientifically verified and efficient therapeutic option thus far-were released, mutations combined into variants of SARS-CoV-2 that are more transmissible and virulent emerged, raising doubts about their efficiency. This study aims to explain possible molecular mechanisms responsible for the increased transmissibility and the increased rate of hospitalizations related to the new variants. A combination of theoretical methods was employed. Constant-pH Monte Carlo simulations were carried out to quantify the stability of several spike trimeric structures at different conformational states and the free energy of interactions between the receptor-binding domain (RBD) and angiotensin-converting enzyme II (ACE2) for the most worrying variants. Electrostatic epitopes were mapped using the PROCEEDpKa method. These analyses showed that the increased virulence is more likely to be due to the improved stability to the S trimer in the opened state, in which the virus can interact with the cellular receptor, ACE2, rather than due to alterations in the complexation RBD-ACE2, since the difference observed in the free energy values was small (although more attractive in general). Conversely, the South African/Beta variant (B.1.351), compared with the SARS-CoV-2 wild type (wt), is much more stable in the opened state with one or two RBDs in the up position than in the closed state with three RBDs in the down position favoring the infection. Such results contribute to understanding the natural history of disease and indicate possible strategies for developing new therapeutic molecules and adjusting the vaccine doses for higher B-cell antibody production.