Computational Analysis of SAM Analogs as Methyltransferase Inhibitors of nsp16/nsp10 Complex from SARS-CoV-2

Methyltransferases (MTases) enzymes, responsible for RNA capping into severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), are emerging important targets for the design of new anti-SARS-CoV-2 agents. Here, analogs of S-adenosylmethionine (SAM), obtained from the bioisosteric substitution of the sulfonium and amino acid groups, were evaluated by rigorous computational modeling techniques such as molecular dynamics (MD) simulations followed by relative binding free analysis against nsp16/nsp10 complex from SARS-CoV-2. The most potent inhibitor (2a) shows the lowest binding free energy (–58.75 Kcal/mol) and more potency than Sinefungin (SFG) (–39.8 Kcal/mol), a pan-MTase inhibitor, which agrees with experimental observations. Besides, our results suggest that the total binding free energy of each evaluated SAM analog is driven by van der Waals interactions which can explain their poor cell permeability, as observed in experimental essays. Overall, we provide a structural and energetic analysis for the inhibition of the nsp16/nsp10 complex involving the evaluated SAM analogs as potential inhibitors.     

Exploring the Catalytic Mechanism of the RNA Cap Modification by nsp16-nsp10 Complex of SARS-CoV-2 through a QM/MM Approach

The inhibition of key enzymes that may contain the viral replication of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have assumed central importance in drug discovery projects. Nonstructural proteins (nsps) are essential for RNA capping and coronavirus replication since it protects the virus from host innate immune restriction. In particular, nonstructural protein 16 (nsp16) in complex with nsp10 is a Cap-0 binding enzyme. The heterodimer formed by nsp16-nsp10 methylates the 5′-end of virally encoded mRNAs to mimic cellular mRNAs and thus it is one of the enzymes that is a potential target for antiviral therapy. In this study, we have evaluated the mechanism of the 2′-O methylation of the viral mRNA cap using hybrid quantum mechanics/molecular mechanics (QM/MM) approach. It was found that the calculated free energy barriers obtained at M062X/6-31+G(d,p) is in agreement with experimental observations. Overall, we provide a detailed molecular analysis of the catalytic mechanism involving the 2′-O methylation of the viral mRNA cap and, as expected, the results demonstrate that the TS stabilization is critical for the catalysis.     

Binding Properties of RNA Quadruplex of SARS-CoV-2 to Berberine Compared to Telomeric DNA Quadruplex

Previous studies suggest that berberine, an isoquinoline alkaloid, has antiviral potential and is a possible therapeutic candidate against SARS-CoV-2. The molecular underpinnings of its action are still unknown. Potential targets include quadruplexes (G4Q) in the viral genome as they play a key role in modulating the biological activity of viruses. While several DNA-G4Q structures and their binding properties have been elucidated, RNA-G4Qs such as RG-1 of the N-gene of SARS-CoV-2 are less explored. Using biophysical techniques, the berberine binding thermodynamics and the associated conformational and hydration changes of RG-1 could be characterized and compared with human telomeric DNA-G4Q 22AG. Berberine can interact with both quadruplexes. Substantial changes were observed in the interaction of berberine with 22AG and RG-1, which adopt different topologies that can also change upon ligand binding. The strength of interaction and the thermodynamic signatures were found to dependent not only on the initial conformation of the quadruplex, but also on the type of salt present in solution. Since berberine has shown promise as a G-quadruplex stabilizer that can modulate viral gene expression, this study may also contribute to the development of optimized ligands that can discriminate between binding to DNA and RNA G-quadruplexes.     

State-of-the-Art Molecular Dynamics Simulation Studies of RNA-Dependent RNA Polymerase of SARS-CoV-2

Molecular dynamics (MD) simulations are powerful theoretical methods that can reveal biomolecular properties, such as structure, fluctuations, and ligand binding, at the level of atomic detail. In this review article, recent MD simulation studies on these biomolecular properties of the RNA-dependent RNA polymerase (RdRp), which is a multidomain protein, of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are presented. Although the tertiary structures of RdRps in SARS-CoV-2 and SARS-CoV are almost identical, the RNA synthesis activity of RdRp of SARS-CoV is higher than SARS-CoV-2. Recent MD simulations observed a difference in the dynamic properties of the two RdRps, which may cause activity differences. RdRp is also a drug target for Coronavirus disease 2019 (COVID-19). Nucleotide analogs, such as remdesivir and favipiravir, are considered to be taken up by RdRp and inhibit RNA replication. Recent MD simulations revealed the recognition mechanism of RdRp for these drug molecules and adenosine triphosphate (ATP). The ligand-recognition ability of RdRp decreases in the order of remdesivir, favipiravir, and ATP. As a typical recognition process, it was found that several lysine residues of RdRp transfer these ligand molecules to the binding site such as a “bucket brigade.” This finding will contribute to understanding the mechanism of the efficient ligand recognition by RdRp. In addition, various simulation studies on the complexes of SARS-CoV-2 RdRp with several nucleotide analogs are reviewed, and the molecular mechanisms by which these compounds inhibit the function of RdRp are discussed. The simulation studies presented in this review will provide useful insights into how nucleotide analogs are recognized by RdRp and inhibit the RNA replication.     

Evolution of SARS-CoV-2 in Spain during the First Two Years of the Pandemic: Circulating Variants,Amino Acid Conservation,and Genetic Variability in Structural,Non-Structural,and Accessory Proteins

Monitoring SARS-CoV-2’s genetic diversity and emerging mutations in this ongoing pandemic is crucial to understanding its evolution and ensuring the performance of COVID-19 diagnostic tests, vaccines, and therapies. Spain has been one of the main epicenters of COVID-19, reaching the highest number of cases and deaths per 100,000 population in Europe at the beginning of the pandemic. This study aims to investigate the epidemiology of SARS-CoV-2 in Spain and its 18 Autonomous Communities across the six epidemic waves established from February 2020 to January 2022. We report on the circulating SARS-CoV-2 variants in each epidemic wave and Spanish region and analyze the mutation frequency, amino acid (aa) conservation, and most frequent aa changes across each structural/non-structural/accessory viral protein among the Spanish sequences deposited in the GISAID database during the study period. The overall SARS-CoV-2 mutation frequency was 1.24 × 10⁻⁵. The aa conservation was >99% in the three types of protein, being non-structural the most conserved. Accessory proteins had more variable positions, while structural proteins presented more aa changes per sequence. Six main lineages spread successfully in Spain from 2020 to 2022. The presented data provide an insight into the SARS-CoV-2 circulation and genetic variability in Spain during the first two years of the pandemic.     

Nucleoside Analogs and Perylene Derivatives Modulate Phase Separation of SARS-CoV-2 N Protein and Genomic RNA In Vitro

The life cycle of severe acute respiratory syndrome coronavirus 2 includes several steps that are supposedly mediated by liquid–liquid phase separation (LLPS) of the viral nucleocapsid protein (N) and genomic RNA. To facilitate the rational design of LLPS-targeting therapeutics, we modeled N-RNA biomolecular condensates in vitro and analyzed their sensitivity to several small-molecule antivirals. The model condensates were obtained and visualized under physiological conditions using an optimized RNA sequence enriched with N-binding motifs. The antivirals were selected based on their presumed ability to compete with RNA for specific N sites or interfere with non-specific pi–pi/cation–pi interactions. The set of antivirals included fleximers, 5′-norcarbocyclic nucleoside analogs, and perylene-harboring nucleoside analogs as well as non-nucleoside amphiphilic and hydrophobic perylene derivatives. Most of these antivirals enhanced the formation of N-RNA condensates. Hydrophobic perylene derivatives and 5′-norcarbocyclic derivatives caused up to 50-fold and 15-fold enhancement, respectively. Molecular modeling data argue that hydrophobic compounds do not hamper specific N-RNA interactions and may promote non-specific ones. These findings shed light on the determinants of potent small-molecule modulators of viral LLPS.     

Exploring Mitochondrial Localization of SARS-CoV-2 RNA by Padlock Assay: A Pilot Study in Human Placenta

The ongoing COVID-19 pandemic dictated new priorities in biomedicine research. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19, is a single-stranded positive-sense RNA virus. In this pilot study, we optimized our padlock assay to visualize genomic and subgenomic regions using formalin-fixed paraffin-embedded placental samples obtained from a confirmed case of COVID-19. SARS-CoV-2 RNA was localized in trophoblastic cells. We also checked the presence of the virion by immunolocalization of its glycoprotein spike. In addition, we imaged mitochondria of placental villi keeping in mind that the mitochondrion has been suggested as a potential residence of the SARS-CoV-2 genome. We observed a substantial overlapping of SARS-CoV-2 RNA and mitochondria in trophoblastic cells. This intriguing linkage correlated with an aberrant mitochondrial network. Overall, to the best of our knowledge, this is the first study that provides evidence of colocalization of the SARS-CoV-2 genome and mitochondria in SARS-CoV-2 infected tissue. These findings also support the notion that SARS-CoV-2 infection can reprogram mitochondrial activity in the highly specialized maternal–fetal interface.     

In Silico Characterization of Masitinib Interaction with SARS-CoV-2 Main Protease

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection continues to be a global health problem. Despite the current implementation of COVID-19 vaccination schedules, identifying effective antiviral drug treatments for this disease continues to be a priority. A recent study showed that masitinib (MST), a tyrosine kinase inhibitor, blocks the proteolytic activity of SARS-CoV-2 main protease (M^pro). Although MST is a potential candidate for COVID-19 treatment, a comprehensive analysis of its interaction with M^pro has not been done. In this work, we performed molecular dynamics simulations of the MST-M^pro complex crystal structure. The effect of the protonation states of M^pro H163 residue and MST titratable groups were studied. Furthermore, we identified the MST substituents and M^pro mutations that affect the stability of the complex. Our results provide valuable insights into the design of new MST analogs as potential treatments for COVID-19.     

SARS-CoV-2 Variants,RBD Mutations,Binding Affinity,and Antibody Escape

Since 2020, the receptor-binding domain (RBD) of the spike protein of the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been constantly mutating, producing most of the notable missense mutations in the context of “variants of concern”, probably in response to the vaccine-driven alteration of immune profiles of the human population. The Delta variant, in particular, has become the most prevalent variant of the epidemic, and it is spreading in countries with the highest vaccination rates, causing the world to face the risk of a new wave of the contagion. Understanding the physical mechanism responsible for the mutation-induced changes in the RBD’s binding affinity, its transmissibility, and its capacity to escape vaccine-induced immunity is the “urgent challenge” in the development of preventive measures, vaccines, and therapeutic antibodies against the coronavirus disease 2019 (COVID-19) pandemic. In this study, entropy–enthalpy compensation and the Gibbs free energy change were used to analyze the impact of the RBD mutations on the binding affinity of SARS-CoV-2 variants with the receptor angiotensin converting enzyme 2 (ACE2) and existing antibodies. Through the analysis, we found that the existing mutations have already covered almost all possible detrimental mutations that could result in an increase of transmissibility, and that a possible mutation in amino-acid position 498 of the RBD can potentially enhance its binding affinity. A new calculation method for the binding energies of protein–protein complexes is proposed based on the entropy–enthalpy compensation rule. All known structures of RBD–antibody complexes and the RBD–ACE2 complex comply with the entropy–enthalpy compensation rule in providing the driving force behind the spontaneous protein–protein docking. The variant-induced risk of breakthrough infections in vaccinated people is attributed to the L452R mutation’s reduction of the binding affinity of many antibodies. Mutations reversing the hydrophobic or hydrophilic performance of residues in the spike RBD potentially cause breakthrough infections of coronaviruses due to the changes in geometric complementarity in the entropy–enthalpy compensations between antibodies and the virus at the binding sites.     

SARS-CoV-2, Cardiovascular Diseases,and Noncoding RNAs: A Connected Triad

Coronavirus Disease 2019 (COVID-19), caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), is characterized by important respiratory impairments frequently associated with severe cardiovascular damages. Moreover, patients with pre-existing comorbidity for cardiovascular diseases (CVD) often present a dramatic increase in inflammatory cytokines release, which increases the severity and adverse outcomes of the infection and, finally, mortality risk. Despite this evident association at the clinical level, the mechanisms linking CVD and COVID-19 are still blurry and unresolved. Noncoding RNAs (ncRNAs) are functional RNA molecules transcribed from DNA but usually not translated into proteins. They play an important role in the regulation of gene expression, either in relatively stable conditions or as a response to different stimuli, including viral infection, and are therefore considered a possible important target in the design of specific drugs. In this review, we introduce known associations and interactions between COVID-19 and CVD, discussing the role of ncRNAs within SARS-CoV-2 infection from the perspective of the development of efficient pharmacological tools to treat COVID-19 patients and taking into account the equally dramatic associated consequences, such as those affecting the cardiovascular system.     

Electrostatic Map of the SARS-CoV-2 Virion Specifies Binding Sites of the Antiviral Cationic Photosensitizer

Electrostatics is an important part of virus life. Understanding the detailed distribution of charges over the surface of a virus is important to predict its interactions with host cells, antibodies, drugs, and different materials. Using a coarse-grained model of the entire viral envelope developed by D. Korkin and S.-J. Marrink’s scientific groups, we created an electrostatic map of the external surface of SARS-CoV-2 and found a highly heterogeneous distribution of the electrostatic potential field of the viral envelope. Numerous negative patches originate mainly from negatively charged lipid domains in the viral membrane and negatively charged areas on the “stalks” of the spike (S) proteins. Membrane (M) and envelope (E) proteins with the total positive charge tend to colocalize with the negatively charged lipids. In the E protein pentamer exposed to the outer surface, negatively charged glutamate residues and surrounding lipids form a negative electrostatic potential ring around the channel entrance. We simulated the interaction of the antiviral octacationic photosensitizer octakis(cholinyl)zinc phthalocyanine with the surface structures of the entire model virion using the Brownian dynamics computational method implemented in ProKSim software (version r661). All mentioned negatively charged envelope components attracted the photosensitizer molecules and are thus potential targets for reactive oxygen generated in photosensitized reactions.     

Studying SARS-CoV-2 with Fluorescence Microscopy

The COVID-19 pandemic caused by SARS-CoV-2 coronavirus deeply affected the world community. It gave a strong impetus to the development of not only approaches to diagnostics and therapy, but also fundamental research of the molecular biology of this virus. Fluorescence microscopy is a powerful technology enabling detailed investigation of virus–cell interactions in fixed and live samples with high specificity. While spatial resolution of conventional fluorescence microscopy is not sufficient to resolve all virus-related structures, super-resolution fluorescence microscopy can solve this problem. In this paper, we review the use of fluorescence microscopy to study SARS-CoV-2 and related viruses. The prospects for the application of the recently developed advanced methods of fluorescence labeling and microscopy—which in our opinion can provide important information about the molecular biology of SARS-CoV-2—are discussed.     

Computational Study of Helicase from SARS-CoV-2 in RNA-Free and Engaged Form

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the pandemic that broke out in 2020 and continues to be the cause of massive global upheaval. Coronaviruses are positive-strand RNA viruses with a genome of ~30 kb. The genome is replicated and transcribed by RNA-dependent RNA polymerase together with accessory factors. One of the latter is the protein helicase (NSP13), which is essential for viral replication. The recently solved helicase structure revealed a tertiary structure composed of five domains. Here, we investigated NSP13 from a structural point of view, comparing its RNA-free form with the RNA-engaged form by using atomistic molecular dynamics (MD) simulations at the microsecond timescale. Structural analyses revealed conformational changes that provide insights into the contribution of the different domains, identifying the residues responsible for domain–domain interactions in both observed forms. The RNA-free system appears to be more flexible than the RNA-engaged form. This result underlies the stabilizing role of the nucleic acid and the functional core role of these domains.     

Small Structural Proteins E and M Render the SARS-CoV-2 Pseudovirus More Infectious and Reveal the Phenotype of Natural Viral Variants

The SARS-CoV-2 pseudovirus is a commonly used strategy that mimics certain biological functions of the authentic virus by relying on biological legitimacy at the molecular level. Despite the fact that spike (S), envelope (E), and membrane (M) proteins together wrap up the SARS-CoV-2 virion, most of the reported pseudotype viruses consist of only the S protein. Here, we report that the presence of E and M increased the virion infectivity by promoting the S protein priming. The S, E, and M (SEM)-coated pseudovirion is spherical, containing crown-like spikes on the surface. Both S and SEM pseudoviruses packaged the same amounts of viral RNA, but the SEM virus bound more efficiently to cells stably expressing the viral receptor human angiotensin-converting enzyme II (hACE2) and became more infectious. Using this SEM pseudovirus, we examined the infectivity and antigenic properties of the natural SARS-CoV-2 variants. We showed that some variants have higher infectivity than the original virus and that some render the neutralizing plasma with lower potency. These studies thus revealed possible mechanisms of the dissemination advantage of these variants. Hence, the SEM pseudovirion provides a useful tool to evaluate the viral infectivity and capability of convalescent sera in neutralizing specific SARS-CoV-2 S dominant variants.     

Glycosylation of Receptor Binding Domain of SARS-CoV-2 S-Protein Influences on Binding to Immobilized DNA Aptamers

Nucleic acid aptamers specific to S-protein and its receptor binding domain (RBD) of SARS-CoV-2 (severe acute respiratory syndrome-related coronavirus 2) virions are of high interest as potential inhibitors of viral infection and recognizing elements in biosensors. Development of specific therapy and biosensors is complicated by an emergence of new viral strains bearing amino acid substitutions and probable differences in glycosylation sites. Here, we studied affinity of a set of aptamers to two Wuhan-type RBD of S-protein expressed in Chinese hamster ovary cell line and Pichia pastoris that differ in glycosylation patterns. The expression system for the RBD protein has significant effects, both on values of dissociation constants and relative efficacy of the aptamer binding. We propose glycosylation of the RBD as the main force for observed differences. Moreover, affinity of a several aptamers was affected by a site of biotinylation. Thus, the robustness of modified aptamers toward new virus variants should be carefully tested.     

Role of RNA Motifs in RNA Interaction with Membrane Lipid Rafts: Implications for Therapeutic Applications of Exosomal RNAs

RNA motifs may promote interactions with exosomes (EXO-motifs) and lipid rafts (RAFT-motifs) that are enriched in exosomal membranes. These interactions can promote selective RNA loading into exosomes. We quantified the affinity between RNA aptamers containing various EXO- and RAFT-motifs and membrane lipid rafts in a liposome model of exosomes by determining the dissociation constants. Analysis of the secondary structure of RNA molecules provided data about the possible location of EXO- and RAFT-motifs within the RNA structure. The affinity of RNAs containing RAFT-motifs (UUGU, UCCC, CUCC, CCCU) and some EXO-motifs (CCCU, UCCU) to rafted liposomes is higher in comparison to aptamers without these motifs, suggesting direct RNA-exosome interaction. We have confirmed these results through the determination of the dissociation constant values of exosome-RNA aptamer complexes. RNAs containing EXO-motifs GGAG or UGAG have substantially lower affinity to lipid rafts, suggesting indirect RNA-exosome interaction via RNA binding proteins. Bioinformatics analysis revealed RNA aptamers containing both raft- and miRNA-binding motifs and involvement of raft-binding motifs UCCCU and CUCCC. A strategy is proposed for using functional RNA aptamers (fRNAa) containing both RAFT-motif and a therapeutic motif (e.g., miRNA inhibitor) to selectively introduce RNAs into exosomes for fRNAa delivery to target cells for personalized therapy.     

A Systematic Study on the Optimal Nucleotide Analogue Concentration and Rate Limiting Nucleotide of the SARS-CoV-2 RNA-Dependent RNA Polymerase

The current COVID-19 pandemic has highlighted the necessity of more efficient antiviral compounds. The antiviral efficacy of adenosine-based analogs, the main repurposed drugs for SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) inhibition, is mainly assessed through in vitro or cell-free polymerization assays, under arbitrary conditions that do not reflect the physiological environment. We show that SARS-CoV-2 RdRp inhibition efficiency of remdesivir and cordycepin, two common adenosine analogs, is influenced by endogenous adenosine level, and that the current clinically approved concentrations for COVID-19 treatment are suboptimal for effective RdRp inhibition. Furthermore, we identified GTP as the rate-limiting nucleotide of SARS-CoV-2 replication. Our results demonstrate that nucleotide sensitivity of the RdRp complex and competition of nucleoside analog drugs against endogenous concentrations of nucleotides are crucial elements to be considered when designing new SARS-CoV-2 antiviral compounds.     

A Collection of Designed Peptides to Target SARS-CoV-2 Spike RBD—ACE2 Interaction

The angiotensin-converting enzyme 2 (ACE2) is the receptor used by SARS-CoV and SARS-CoV-2 coronaviruses to attach to cells via the receptor-binding domain (RBD) of their viral spike protein. Since the start of the COVID-19 pandemic, several structures of protein complexes involving ACE2 and RBD as well as monoclonal antibodies and nanobodies have become available. We have leveraged the structural data to design peptides to target the interaction between the RBD of SARS-CoV-2 and ACE2 and SARS-CoV and ACE2, as contrasting exemplar, as well as the dimerization surface of ACE2 monomers. The peptides were modelled using our original method: PiPreD that uses native elements of the interaction between the targeted protein and cognate partner(s) that are subsequently included in the designed peptides. These peptides recapitulate stretches of residues present in the native interface plus novel and highly diverse conformations surrogating key interactions at the interface. To facilitate the access to this information we have created a freely available and dedicated web-based repository, PepI-Covid19 database, providing convenient access to this wealth of information to the scientific community with the view of maximizing its potential impact in the development of novel therapeutic and diagnostic agents.