Structural modeling of protein ensembles between E3 RING ligases and SARS-CoV-2: The role of zinc binding domains

BackgroundThe ubiquitin system is a modification process with many different cellular functions including immune signaling and antiviral functions. E3 ubiquitin ligases are enzymes that recruit an E2 ubiquitin-conjugating enzyme bound to ubiquitin in order to catalyze the transfer of ubiquitin from the E2 to a protein substrate. The RING E3s, the most abundant type of ubiquitin ligases, are characterized by a zinc (II)-binding domain called RING (Really Interesting New Gene). Viral replication requires modifying and hijacking key cellular pathways within host cells such as cellular ubiquitination. There are well-established examples where a viral proteins bind to RING E3s, redirecting them to degrade otherwise long-lived host proteins or inhibiting E3’s ubiquitination activity. Recently, three binary interactions between SARS-CoV-2 proteins and innate human immune signaling Ε3 RING ligases: NSP15-RNF41, ORF3a-TRIM59 and NSP9-MIB1 have been experimentally established.MethodsIn this work, we have investigated the mode of the previous experimentally supported NSP15-RNF41, ORF3a,-TRIM59 and NSP9-MIB1 binary interactions by in silico methodologies intending to provide structural insights of E3-virus interplay that can help identify potential inhibitors that could block SARS-CoV-2 infection of immune cells.ConclusionIn silico methodologies have shown that the above human E3 ligases interact with viral partners through their Zn(II) binding domains. This RING mediated formation of stable SARS-CoV-2-E3 complexes indicates a critical structural role of RING domains in immune system disruption by SARS-CoV-2-infection.Data AvailabilityThe data used to support the findings of this research are included within the article and are labeled with references.     

The proximal proteome of 17 SARS-CoV-2 proteins links to disrupted antiviral signaling and host translation

Viral proteins localize within subcellular compartments to subvert host machinery and promote pathogenesis. To study SARS-CoV-2 biology, we generated an atlas of 2422 human proteins vicinal to 17 SARS-CoV-2 viral proteins using proximity proteomics. This identified viral proteins at specific intracellular locations, such as association of accessary proteins with intracellular membranes, and projected SARS-CoV-2 impacts on innate immune signaling, ER-Golgi transport, and protein translation. It identified viral protein adjacency to specific host proteins whose regulatory variants are linked to COVID-19 severity, including the TRIM4 interferon signaling regulator which was found proximal to the SARS-CoV-2 M protein. Viral NSP1 protein adjacency to the EIF3 complex was associated with inhibited host protein translation whereas ORF6 localization with MAVS was associated with inhibited RIG-I 2CARD-mediated IFNB1 promoter activation. Quantitative proteomics identified candidate host targets for the NSP5 protease, with specific functional cleavage sequences in host proteins CWC22 and FANCD2. This data resource identifies host factors proximal to viral proteins in living human cells and nominates pathogenic mechanisms employed by SARS-CoV-2.     

Staphylococcus aureus protects its immune‐evasion proteins against degradation by neutrophil serine proteases

Neutrophils store large quantities of neutrophil serine proteases (NSPs) that contribute, via multiple mechanisms, to antibacterial immune defences. Even though neutrophils are indispensable in fighting Staphylococcus aureus infections, the importance of NSPs in anti‐staphylococcal defence is yet unknown. However, the fact that S. aureus produces three highly specific inhibitors for NSPs [the extracellular adherence proteins (EAPs) Eap, EapH1 and EapH2], suggests that these proteases are important for host defences against this bacterium. In this study we demonstrate that NSPs can inactivate secreted virulence factors of S. aureus and that EAP proteins function to prevent this degradation. Specifically, we find that a large group of S. aureus immune‐evasion proteins is vulnerable to proteolytic inactivation by NSPs. In most cases, NSP cleavage leads to functional inactivation of virulence proteins. Interestingly, proteins with similar immune‐escape functions appeared to have differential cleavage sensitivity towards NSPs. Using targeted mutagenesis and complementation analyses in S. aureus, we demonstrate that all EAP proteins can protect other virulence factors from NSP degradation in complex bacterial supernatants. These findings show that NSPs inactivate S. aureus virulence factors. Moreover, the protection by EAP proteins can explain why this antibacterial function of NSPs was masked in previous studies. Furthermore, our results indicate that therapeutic inactivation of EAP proteins can help to restore the natural host immune defences against S. aureus.     

Identification of SARS-CoV-2 RNA-dependent RNA polymerase inhibitors from the major phytochemicals of Nigella sativa: An in silico approach

The coronavirus disease 2019 (COVID-19), which emerged in December 2019, continues to be a serious health concern worldwide. There is an urgent need to develop effective drugs and vaccines to control the spread of this disease. In the current study, the main phytochemical compounds of Nigella sativa were screened for their binding affinity for the active site of the RNA-dependent RNA polymerase (RdRp) enzyme of the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2). The binding affinity was investigated using molecular docking methods, and the interaction of phytochemicals with the RdRp active site was analyzed and visualized using suitable software. Out of the nine phytochemicals of N. sativa screened in this study, a significant docking score was observed for four compounds, namely α-hederin, dithymoquinone, nigellicine, and nigellidine. Based on the findings of our study, we report that α-hederin, which was found to possess the lowest binding energy (–8.6 kcal/mol) and hence the best binding affinity, is the best inhibitor of RdRp of SARS-CoV-2, among all the compounds screened here. Our results prove that the top four potential phytochemical molecules of N. sativa, especially α-hederin, could be considered for ongoing drug development strategies against SARS-CoV-2. However, further in vitro and in vivo testing are required to confirm the findings of this study.     

SARS-CoV-2 suppresses IFNβ production mediated by NSP1, 5, 6, 15, ORF6 and ORF7b but does not suppress the effects of added interferon

Type I Interferons (IFN-Is) are a family of cytokines which play a major role in inhibiting viral infection. Resultantly, many viruses have evolved mechanisms in which to evade the IFN-I response. Here we tested the impact of expression of 27 different SARS-CoV-2 genes in relation to their effect on IFN production and activity using three independent experimental methods. We identified six gene products; NSP6, ORF6, ORF7b, NSP1, NSP5 and NSP15, which strongly (>10-fold) blocked MAVS-induced (but not TRIF-induced) IFNβ production. Expression of the first three of these SARS-CoV-2 genes specifically blocked MAVS-induced IFNβ-promoter activity, whereas all six genes induced a collapse in IFNβ mRNA levels, corresponding with suppressed IFNβ protein secretion. Five of these six genes furthermore suppressed MAVS-induced activation of IFNλs, however with no effect on IFNα or IFNγ production. In sharp contrast, SARS-CoV-2 infected cells remained extremely sensitive to anti-viral activity exerted by added IFN-Is. None of the SARS-CoV-2 genes were able to block IFN-I signaling, as demonstrated by robust activation of Interferon Stimulated Genes (ISGs) by added interferon. This, despite the reduced levels of STAT1 and phospho-STAT1, was likely caused by broad translation inhibition mediated by NSP1. Finally, we found that a truncated ORF7b variant that has arisen from a mutant SARS-CoV-2 strain harboring a 382-nucleotide deletion associating with mild disease (Δ382 strain identified in Singapore & Taiwan in 2020) lost its ability to suppress type I and type III IFN production. In summary, our findings support a multi-gene process in which SARS-CoV-2 blocks IFN-production, with ORF7b as a major player, presumably facilitating evasion of host detection during early infection. However, SARS-CoV-2 fails to suppress IFN-I signaling thus providing an opportunity to exploit IFN-Is as potential therapeutic antiviral drugs.     

Evidence for multiple modes of neutrophil serine protease recognition by the EAP family of Staphylococcal innate immune evasion proteins

Neutrophils contain high levels of chymotrypsin‐like serine proteases (NSPs) within their azurophilic granules that have a multitude of functions within the immune system. In response, the pathogen Staphylococcus aureus has evolved three potent inhibitors (Eap, EapH1, and EapH2) that protect the bacterium as well as several of its secreted virulence factors from the degradative action of NSPs. We previously showed that these so‐called EAP domain proteins represent a novel class of NSP inhibitors characterized by a non‐covalent inhibitory mechanism and a distinct target specificity profile. Based upon high levels of structural homology amongst the EAP proteins and the NSPs, as well as supporting biochemical data, we predicted that the inhibited complex would be similar for all EAP/NSP pairs. However, we present here evidence that EapH1 and EapH2 bind the canonical NSP, Neutrophil Elastase (NE), in distinct orientations. We discovered that alteration of EapH1 residues at the EapH1/NE interface caused a dramatic loss of affinity and inhibition of NE, while mutation of equivalent positions in EapH2 had no effect on NE binding or inhibition. Surprisingly, mutation of residues in an altogether different region of EapH2 severely impacted both the NE binding and inhibitory properties of EapH2. Even though EapH1 and EapH2 bind and inhibit NE and a second NSP, Cathepsin G, equally well, neither of these proteins interacts with the structurally related, but non‐proteolytic granule protein, azurocidin. These studies expand our understanding of EAP/NSP interactions and suggest that members of this immune evasion protein family are capable of diverse target recognition modes.     

Genome composition and genetic characterization of SARS-CoV-2

SARS-CoV-2 is a type of Betacoronaviruses responsible for COVID-19 pandemic disease, with more than 1.745 million fatalities globally as of December-2020. Genetically, it is considered the second largest genome of all RNA viruses with a 5′ cap and 3′ poly-A tail. Phylogenetic analyses of coronaviruses reveal that SARS-CoV-2 is genetically closely related to the Bat-SARS Like-Corona virus (Bat-SL-Cov) with 96% whole-genome identity. SARS-CoV-2 genome consists of 15 ORFs coded into 29 proteins. At the 5′ terminal of the genome, we have ORF1ab and ORF1a, which encode the 1ab and 1a polypeptides that are proteolytically cleaved into 16 different nonstructural proteins (NSPs). The 3′ terminal of the genome represents four structural (spike, envelope, matrix, and nucleocapsid) and nine accessory (3a, 3b, 6, 7a, 7b, 8b, 9a, 9b, and orf10) proteins. As the number of COVID-19 patients increases dramatically worldwide, there is an urgent need to find a quick and sensitive diagnostic tool for controlling the outbreak of SARS-CoV-2 in the community. Today, molecular testing methods utilizing viral genetic material (e.g., PCR) represent the crucial diagnostic tool for the SARS-CoV-2 virus despite its low sensitivity in the early stage of viral infection. This review summarizes the genome composition and genetic characterization of the SARS-CoV-2.     

Potential role of nicotinamide analogues against SARS-COV-2 target proteins

Background and objectiveCoronavirus 2019 (COVID-19) is caused by ‘severe acute respiratory syndrome coronavirus 2′ (SARS-CoV-2), first reported in Wuhan, China in December 2019, which eventually became a global disaster. Various key mediators have been reported in the pathogenesis of COVID-19. However, no effective pharmacological intervention has been available to combat COVID-19 complications. The present study screens nicotinamide riboside (NR) and nicotinamide mononucleotide (NMN) as potential inhibitors of this present generation coronavirus infection using an in-silico approach.Materials and methodsThe SARS-CoV-2 proteins (nucleocapsid, proteases, post-fusion core, phosphatase, endoriboruclease) and ACE-2 protein were selected. The 2D structure of nicotinamide ribonucleoside and nicotinamide ribonucleotide was drawn using ChemDraw 14.0 and saved in .cdx format. The results were analyzed using two parameters: full fitness energy and binding free energy (ΔG).ResultsThe full fitness energy and estimated ΔG values from docking of NM, and NMN with selected SARS-CoV-2 target proteins, ADMET prediction and Target prediction indicate the interaction of NR and NMN in the treatment of COVID-19.ConclusionsBased on full fitness energy and estimated ΔG values from docking studies of NM and NAM with selected SARS-CoV-2 target proteins, ADME prediction, target prediction and toxicity prediction, we expect a possible therapeutic efficacy of NR in the treatment of COVID-19.Keyword: COVID-19, SARS-CoV-2, Molecular docking, Poly (ADP-ribose) polymerase enzyme, Nicotinamide     

Activation of the SARS-CoV-2 NSP14 3′–5′ exoribonuclease by NSP10 and response to antiviral inhibitors

Understanding the core replication complex of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is essential to the development of novel coronavirus-specific antiviral therapeutics. Among the proteins required for faithful replication of the SARS-CoV-2 genome are nonstructural protein 14 (NSP14), a bifunctional enzyme with an N-terminal 3′-to-5′ exoribonuclease (ExoN) and a C-terminal N7-methyltransferase, and its accessory protein, NSP10. The difficulty in producing pure and high quantities of the NSP10/14 complex has hampered the biochemical and structural study of these important proteins. We developed a straightforward protocol for the expression and purification of both NSP10 and NSP14 from Escherichia coli and for the in vitro assembly and purification of a stoichiometric NSP10/14 complex with high yields. Using these methods, we observe that NSP10 provides a 260-fold increase in k_cat/K_m in the exoribonucleolytic activity of NSP14 and enhances protein stability. We also probed the effect of two small molecules on NSP10/14 activity, remdesivir monophosphate and the methyltransferase inhibitor S-adenosylhomocysteine. Our analysis highlights two important factors for drug development: first, unlike other exonucleases, the monophosphate nucleoside analog intermediate of remdesivir does not inhibit NSP14 activity; and second, S-adenosylhomocysteine modestly activates NSP14 exonuclease activity. In total, our analysis provides insights for future structure–function studies of SARS-CoV-2 replication fidelity for the treatment of coronavirus disease 2019.     

新型冠状病毒相关靶点小分子抑制剂虚拟筛选

为确定治疗新型冠状病毒(SARS-CoV-2)感染的候选药物,开展了针对SARS-CoV-2的药物虚拟筛选研究。以SARS-CoV-2的刺突蛋白(S蛋白)和3CL蛋白酶(主蛋白酶)作为药物靶点,以美国食品药品监督管理局(FDA)批准上市的2 726个小分子药物作为候选,通过分子对接方法,筛选出了3种(阿巴瑞克(Abarelix)、西曲瑞克(Cetrorelix)、鞣酸(Tannic acid))与S蛋白具有较强结合能力的小分子药物,1种(戈舍瑞林(Goserelin))与3CL蛋白酶具有较好结合能力的小分子药物,它们理论上都具有抑制新型冠状病毒复制的效果。将靶向3CL蛋白酶的候选药物与辉瑞公司开发的药物Paxlovid进行比较,发现其作用位点均集中于3CL蛋白酶的第130-200位的残基周围,具有相似的结合位点与相互作用。此外也对候选药物的物理与化学性质及与基因相互作用进行了分析。本研究可为开发新型冠状病毒感染的治疗药物提供参考。     

Clinical and genomic data of sars-cov-2 detected in maternal–fetal interface during the first wave of infection in Brazil

Brazil has the highest SARS-CoV-2 case-fatality rate in pregnant women in the Americas. In this study, clinical and virological findings of five mildly symptomatic pregnant women and their infected fetuses/newborns treated at a referral hospital for COVID19-pregnant women in Midwestern Brazil are reported. Mother and fetal samples were tested by RT-qPCR, ECLIA and Illumina MiSeq sequencing. From the five cases, one resulted in spontaneous abortion, one was stillborn, two were preterm births and one full-term birth. Maternal and fetal placenta, newborn and stillborn secretions were SARS-CoV-2+; one neonate developed ground-glass opacities in his lungs. One neonate's umbilical cord was IgG+ and all were IgM negative upon hospital discharge. Genomes recovered from two placentas belong to the B.1.1.28 and B.1.1.33 lineages and present nonsynonymous mutations associated with virus fitness and infectivity; other not frequently reported mutations (B.1.1.33: NSP3 V2090G, M A2S and ORF3ab S253P and Y264N; B.1.1.28: NSP3 E995D, NSP12 R240K, NSP14H1897Y and in ORF7b V21F) were found in proteins involved in viral replication, viral induction of apoptosis, viral interference on interferon and on NF-Κβ pathways. Phylogeny indicates the south of Brazil as the possible origin of these lineages circulating in Mato Grosso State. These findings contribute to describe SARS-CoV-2 infection and outcomes in pregnant women and their fetuses, at any stage of gestation and even in mild symptomatic cases.     

Caffeic acid derivatives (CAFDs) as inhibitors of SARS-CoV-2: CAFDs-based functional foods as a potential alternative approach to combat COVID-19

BackgroundSARS-CoV-2, an emerging strain of coronavirus, has affected millions of people from all the continents of world and received worldwide attention. This emerging health crisis calls for the urgent development of specific therapeutics against COVID-19 to potentially reduce the burden of this emerging pandemic.PurposeThis study aims to evaluate the anti-viral efficacy of natural bioactive entities against COVID-19 via molecular docking and molecular dynamics simulation.MethodsA library of 27 caffeic-acid derivatives was screened against 5 proteins of SARS-CoV-2 by using Molegro Virtual Docker 7 to obtain the binding energies and interactions between compounds and SARS-CoV-2 proteins. ADME properties and toxicity profiles were investigated via www.swissadme.ch web tools and Toxtree respectively. Molecular dynamics simulation was performed to determine the stability of the lead-protein interactions.ResultsOur obtained results has uncovered khainaoside C, 6-O-Caffeoylarbutin, khainaoside B, khainaoside C and vitexfolin A as potent modulators of COVID-19 possessing more binding energies than nelfinavir against COVID-19 M^pro, Nsp15, SARS-CoV-2 spike S2 subunit, spike open state and closed state structure respectively. While Calceolarioside B was identified as pan inhibitor, showing strong molecular interactions with all proteins except SARS-CoV-2 spike glycoprotein closed state. The results are supported by 20 ns molecular dynamics simulations of the best complexes.ConclusionThis study will hopefully pave a way for development of phytonutrients-based antiviral therapeutic for treatment or prevention of COVID-19 and further studies are recommended to evaluate the antiviral effects of these phytochemicals against SARS-CoV-2 in in vitro and in vivo models.     

基于处方挖掘与分子动力学模拟筛选严重急性呼吸综合征冠状病毒2潜在抑制剂分子的研究

目的 从中药筛选具有潜在抑制严重急性呼吸综合征冠状病毒2 （SARS-CoV-2） 活性的成分，进一步从原子水平揭 示其抑制SARS-CoV-2 表面刺突蛋白（S 蛋白） 受体结合域（RBD） 与血管紧张素转化酶2 （ACE2） 结合的内在机制。 方法 检索新型冠状病毒（简称“新冠肺炎”） 治疗中药处方，构建“新冠肺炎中药候选活性成分数据库”。用具有ACE2 抑制活性的小分子化合物构建HipHop药效团模型，并对“新冠肺炎中药候选活性成分数据库”中活性成分筛选。采用分子 对接和分子动力学模拟方法研究候选活性成分与ACE2 的结合方式及其对SARS-CoV-2 S 蛋白与ACE2 识别的影响。 结果 本文通过中药处方挖掘和分子动力学模拟，从143 个新冠肺炎治疗中药处方中筛选出10 种可与SARS-CoV-2 S 蛋白/ 人源ACE2 识别位点结合的中药成分。其中，枇杷叶主要活性成分23-trans-p-coumaryhormentic acid 与ACE2 具有最高的亲和 力，且23-trans-p-coumaryhormentic acid 的结合可有效阻断SARS-CoV-2 S蛋白与宿主细胞ACE2 的结合。结论 本文通过虚 拟筛选发现了SARS-CoV-2 潜在抑制剂分子23-trans-p-coumaryhormentic acid，同时从原子水平预测了其抑制SARS-CoV-2 S 蛋白与ACE2 结合的内在机制，这将为SARS-CoV-2 特异性抗病毒药物的研发提供理论依据。     

Structural analysis of protein‐ligand interactions: the binding of endogenous compounds and of synthetic drugs

The large number of macromolecular structures deposited with the Protein Data Bank (PDB) describing complexes between proteins and either physiological compounds or synthetic drugs made it possible a systematic analysis of the interactions occurring between proteins and their ligands. In this work, the binding pockets of about 4000 PDB protein‐ligand complexes were investigated and amino acid and interaction types were analyzed. The residues observed with lowest frequency in protein sequences, Trp, His, Met, Tyr, and Phe, turned out to be the most abundant in binding pockets. Significant differences between drug‐like and physiological compounds were found. On average, physiological compounds establish with respect to drugs about twice as many hydrogen bonds with protein atoms, whereas drugs rely more on hydrophobic interactions to establish target selectivity. The large number of PDB structures describing homologous proteins in complex with the same ligand made it possible to analyze the conservation of binding pocket residues among homologous protein structures bound to the same ligand, showing that Gly, Glu, Arg, Asp, His, and Thr are more conserved than other amino acids. Also in the cases in which the same ligand is bound to unrelated proteins, the binding pockets showed significant conservation in the residue types. In this case, the probability of co‐occurrence of the same amino acid type in the binding pockets could be up to thirteen times higher than that expected on a random basis. The trends identified in this study may provide an useful guideline in the process of drug design and lead optimization. Copyright © 2014 John Wiley & Sons, Ltd.     

Investigation of Binding Phenomenon of NSP3 and p130Cas Mutants and Their Effect on Cell Signalling

Members of the novel SH2-containing protein (NSP3) and Crk-associated substrate (p130Cas) protein families form a multi-domain signalling platforms that mediate cell signalling process. We analysed the damaging consequences of three mutations, each from NSP3 (NSP3^L469R, NSP3^L623E, NSP3^R627E) and p130Cas (p130Cas^F794R, p130Cas^L787E, p130Cas^D797R) protein with respect to their native biological partners. Mutations depicted notable loss in interaction affinity towards their corresponding biological partners. NSP3^L469R and p130Cas^D797R mutations were predicted as most prominent in docking analysis. Molecular dynamics (MD) studies were conducted to evaluate structural consequences of most prominent mutation in NSP3 and p130Cas obtained from the docking analysis. MD analysis confirmed that mutation in NSP3^L469R and p130Cas^D797R showed significant structural deviation, changes in conformations and increased flexibility, which in turn affected the binding affinity with their biological partners. Moreover, the root mean square fluctuation has indicated a rise in fluctuation of residues involved in moderate interaction acquired between the NSP3 and p130Cas. It has significantly affected the binding interaction in mutant complexes. The results obtained in this work present a detailed overview of molecular mechanisms involved in the loss of cell signalling associated with NSP3 and p130Cas protein.