Evolution of Stronger SARS-CoV-2 Variants as Revealed Through the Lens of Molecular Dynamics Simulations.

Using molecular dynamics simulations, the protein-protein interactions of the receptor-binding domain of the wild-type and seven variants of the severe acute respiratory syndrome coronavirus 2 spike protein and the peptidase domain of human angiotensin-converting enzyme 2 were investigated. These variants are alpha, beta, gamma, delta, eta, kappa, and omicron. Using 100 ns simulation data, the residue interaction networks at the protein-protein interface were identified. Also, the impact of mutations on essential protein dynamics, backbone flexibility, and interaction energy of the simulated protein-protein complexes were studied. The protein-protein interface for the wild-type, delta, and omicron variants contained several stronger interactions, while the alpha, beta, gamma, eta, and kappa variants exhibited an opposite scenario as evident from the analysis of the inter-residue interaction distances and pair-wise interaction energies. The study reveals that two distinct residue networks at the central and right contact regions forge stronger binding affinity between the protein partners. The study provides a molecular-level insight into how enhanced transmissibility and infectivity by delta and omicron variants are most likely tied to a handful of interacting residues at the binding interface, which could potentially be utilized for future antibody constructs and structure-based antiviral drug design.

PCN-Miner: an open-source extensible tool for the analysis of Protein Contact Networks.

MOTIVATION: Protein Contact Network (PCN) is a powerful method for analysing the structure and function of proteins, with a specific focus on disclosing the molecular features of allosteric regulation through the discovery of modular substructures. The importance of PCN analysis has been shown in many contexts, such as the analysis of SARS-CoV-2 Spike protein and its complexes with the Angiotensin Converting Enzyme 2 (ACE2) human receptors. Even if there exist many software tools implementing such methods, there is a growing need for the introduction of tools integrating existing approaches. RESULTS: We present PCN-Miner, a software tool implemented in the Python programming language, able to (i) import protein structures from the Protein Data Bank; (ii) generate the corresponding PCN; (iii) model, analyse and visualize PCNs and related protein structures by using a set of known algorithms and metrics. The PCN-Miner can cover a large set of applications: from clustering to embedding and subsequent analysis. AVAILABILITY AND IMPLEMENTATION: The PCN-Miner tool is freely available at the following GitHub repository: https://github.com/hguzzi/ProteinContactNetworks. It is also available in the Python Package Index (PyPI) repository.

PEPPI: Whole-proteome Protein-protein Interaction Prediction through Structure and Sequence Similarity, Functional Association, and Machine Learning.

Proteome-wide identification of protein-protein interactions is a formidable task which has yet to be sufficiently addressed by experimental methodologies. Many computational methods have been developed to predict proteome-wide interaction networks, but few leverage both the sensitivity of structural information and the wide availability of sequence data. We present PEPPI, a pipeline which integrates structural similarity, sequence similarity, functional association data, and machine learning-based classification through a naïve Bayesian classifier model to accurately predict protein-protein interactions at a proteomic scale. Through benchmarking against a set of 798 ground truth interactions and an equal number of non-interactions, we have found that PEPPI attains 4.5% higher AUROC than the best of other state-of-the-art methods. As a proteomic-scale application, PEPPI was applied to model the interactions which occur between SARS-CoV-2 and human host cells during coronavirus infection, where 403 high-confidence interactions were identified with predictions covering 73% of a gold standard dataset from PSICQUIC and demonstrating significant complementarity with the most recent high-throughput experiments. PEPPI is available both as a webserver and in a standalone version and should be a powerful and generally applicable tool for computational screening of protein-protein interactions.

Baicalin Targets HSP70/90 to Regulate PKR/PI3K/AKT/eNOS Signaling Pathways.

Baicalin is a major active ingredient of traditional Chinese medicine Scutellaria baicalensis, and has been shown to have antiviral, anti-inflammatory, and antitumor activities. However, the protein targets of baicalin have remained unclear. Herein, a chemical proteomics strategy was developed by combining baicalin-functionalized magnetic nanoparticles (BCL-N3@MNPs) and quantitative mass spectrometry to identify the target proteins of baicalin. Bioinformatics analysis with the use of Gene Ontology, STRING and Ingenuity Pathway Analysis, was performed to annotate the biological functions and the associated signaling pathways of the baicalin targeting proteins. Fourteen proteins in human embryonic kidney cells were identified to interact with baicalin with various binding affinities. Bioinformatics analysis revealed these proteins are mainly ATP-binding and/or ATPase activity proteins, such as CKB, HSP86, HSP70-1, HSP90, ATPSF1ß and ACTG1, and highly associated with the regulation of the role of PKR in interferon induction and the antiviral response signaling pathway (P = 10-6), PI3K/AKT signaling pathway (P = 10-5) and eNOS signaling pathway (P = 10-4). The results show that baicalin exerts multiply pharmacological functions, such as antiviral, anti-inflammatory, antitumor, and antioxidant functions, through regulating the PKR and PI3K/AKT/eNOS signaling pathways by targeting ATP-binding and ATPase activity proteins. These findings provide a fundamental insight into further studies on the mechanism of action of baicalin.

Competing endogenous RNA network mediated by circ_3205 in SARS-CoV-2 infected cells.

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is a new member of the Betacoronaviridae family, responsible for the recent pandemic outbreak of COVID-19. To start exploring the molecular events that follow host cell infection, we queried VirusCircBase and identified a circular RNA (circRNA) predicted to be synthesized by SARS-CoV-2, circ_3205, which we used to probe: (i) a training cohort comprised of two pools of cells from three nasopharyngeal swabs of SARS-CoV-2 infected (positive) or uninfected (negative, UCs) individuals; (ii) a validation cohort made up of 12 positive and 3 negative samples. The expression of circRNAs, miRNAs and miRNA targets was assayed through real-time PCR. CircRNA-miRNA interactions were predicted by TarpMiR, Analysis of Common Targets for circular RNAs (ACT), and STarMir tools. Enrichment of the biological processes and the list of predicted miRNA targets were retrieved from DIANA miRPath v3.0. Our results showed that the predicted SARS-CoV-2 circ_3205 was expressed only in positive samples and its amount positively correlated with that of SARS-CoV-2 Spike (S) mRNA and the viral load (r values = 0.80952 and 0.84867, Spearman's correlation test, respectively). Human (hsa) miR-298 was predicted to interact with circ_3205 by all three predictive tools. KCNMB4 and PRKCE were predicted as hsa-miR-298 targets. Interestingly, the function of both is correlated with blood coagulation and immune response. KCNMB4 and PRKCE mRNAs were upregulated in positive samples as compared to UCs (6 and 8.1-fold, p values = 0.049 and 0.02, Student's t test, respectively) and their expression positively correlated with that of circ_3205 (r values = 0.6 and 0.25, Spearman's correlation test, respectively). We propose that our results convincingly suggest that circ_3205 is a circRNA synthesized by SARS-CoV-2 upon host cell infection and that it may behave as a competitive endogenous RNA (ceRNA), sponging hsa-miR-298 and contributing to the upregulation of KCNMB4 and PRKCE mRNAs.

Inflammatory activation and immune cell infiltration are main biological characteristics of SARS-CoV-2 infected myocardium.

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) can target cardiomyocytes (CMs) to directly invade the heart resulting in high mortality. This study aims to explore the biological characteristics of SARS-CoV-2 infected myocardium based on omics by collecting transcriptome data and analyzing them with a series of bioinformatics tools. Totally, 86 differentially expressed genes (DEGs) were discovered in SARS-CoV-2 infected CMs, and 15 miRNAs were discovered to target 60 genes. Functional enrichment analysis indicated that these DEGs were mainly enriched in the inflammatory signaling pathway. After the protein-protein interaction (PPI) network was constructed, several genes including CCL2 and CXCL8 were regarded as the hub genes. SRC inhibitor saracatinib was predicted to potentially act against the cardiac dysfunction induced by SARS-CoV-2. Among the 86 DEGs, 28 were validated to be dysregulated in SARS-CoV-2 infected hearts. Gene Set Enrichment Analysis (GSEA) analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) showed that malaria, IL-17 signaling pathway, and complement and coagulation cascades were significantly enriched. Immune infiltration analysis indicated that 'naive B cells' was significantly increased in the SARS-CoV-2 infected heart. The above results may help to improve the prognosis of patients with COVID-19.

SARS-CoV-2 NSP13 Inhibits Type I IFN Production by Degradation of TBK1 via p62-Dependent Selective Autophagy.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which causes coronavirus disease 2019 (COVID-19), has seriously threatened global public health. Severe COVID-19 has been reported to be associated with an impaired IFN response. However, the mechanisms of how SARS-CoV-2 antagonizes the host IFN response are poorly understood. In this study, we report that SARS-CoV-2 helicase NSP13 inhibits type I IFN production by directly targeting TANK-binding kinase 1 (TBK1) for degradation. Interestingly, inhibition of autophagy by genetic knockout of Beclin1 or pharmacological inhibition can rescue NSP13-mediated TBK1 degradation in HEK-293T cells. Subsequent studies revealed that NSP13 recruits TBK1 to p62, and the absence of p62 can also inhibit TBK1 degradation in HEK-293T and HeLa cells. Finally, TBK1 and p62 degradation and p62 aggregation were observed during SARS-CoV-2 infection in HeLa-ACE2 and Calu3 cells. Overall, our study shows that NSP13 inhibits type I IFN production by recruiting TBK1 to p62 for autophagic degradation, enabling it to evade the host innate immune response, which provides new insights into the transmission and pathogenesis of SARS-CoV-2 infection.

Super.Complex: A supervised machine learning pipeline for molecular complex detection in protein-interaction networks.

Characterization of protein complexes, i.e. sets of proteins assembling into a single larger physical entity, is important, as such assemblies play many essential roles in cells such as gene regulation. From networks of protein-protein interactions, potential protein complexes can be identified computationally through the application of community detection methods, which flag groups of entities interacting with each other in certain patterns. Most community detection algorithms tend to be unsupervised and assume that communities are dense network subgraphs, which is not always true, as protein complexes can exhibit diverse network topologies. The few existing supervised machine learning methods are serial and can potentially be improved in terms of accuracy and scalability by using better-suited machine learning models and parallel algorithms. Here, we present Super.Complex, a distributed, supervised AutoML-based pipeline for overlapping community detection in weighted networks. We also propose three new evaluation measures for the outstanding issue of comparing sets of learned and known communities satisfactorily. Super.Complex learns a community fitness function from known communities using an AutoML method and applies this fitness function to detect new communities. A heuristic local search algorithm finds maximally scoring communities, and a parallel implementation can be run on a computer cluster for scaling to large networks. On a yeast protein-interaction network, Super.Complex outperforms 6 other supervised and 4 unsupervised methods. Application of Super.Complex to a human protein-interaction network with ~8k nodes and ~60k edges yields 1,028 protein complexes, with 234 complexes linked to SARS-CoV-2, the COVID-19 virus, with 111 uncharacterized proteins present in 103 learned complexes. Super.Complex is generalizable with the ability to improve results by incorporating domain-specific features. Learned community characteristics can also be transferred from existing applications to detect communities in a new application with no known communities. Code and interactive visualizations of learned human protein complexes are freely available at: https://sites.google.com/view/supercomplex/super-complex-v3-0.

SCoV2-MD: a database for the dynamics of the SARS-CoV-2 proteome and variant impact predictions.

SCoV2-MD (www.scov2-md.org) is a new online resource that systematically organizes atomistic simulations of the SARS-CoV-2 proteome. The database includes simulations produced by leading groups using molecular dynamics (MD) methods to investigate the structure-dynamics-function relationships of viral proteins. SCoV2-MD cross-references the molecular data with the pandemic evolution by tracking all available variants sequenced during the pandemic and deposited in the GISAID resource. SCoV2-MD enables the interactive analysis of the deposited trajectories through a web interface, which enables users to search by viral protein, isolate, phylogenetic attributes, or specific point mutation. Each mutation can then be analyzed interactively combining static (e.g. a variety of amino acid substitution penalties) and dynamic (time-dependent data derived from the dynamics of the local geometry) scores. Dynamic scores can be computed on the basis of nine non-covalent interaction types, including steric properties, solvent accessibility, hydrogen bonding, and other types of chemical interactions. Where available, experimental data such as antibody escape and change in binding affinities from deep mutational scanning experiments are also made available. All metrics can be combined to build predefined or custom scores to interrogate the impact of evolving variants on protein structure and function.

HumanNet v3: an improved database of human gene networks for disease research.

Network medicine has proven useful for dissecting genetic organization of complex human diseases. We have previously published HumanNet, an integrated network of human genes for disease studies. Since the release of the last version of HumanNet, many large-scale protein-protein interaction datasets have accumulated in public depositories. Additionally, the numbers of research papers and functional annotations for gene-phenotype associations have increased significantly. Therefore, updating HumanNet is a timely task for further improvement of network-based research into diseases. Here, we present HumanNet v3 (https://www.inetbio.org/humannet/, covering 99.8% of human protein coding genes) constructed by means of the expanded data with improved network inference algorithms. HumanNet v3 supports a three-tier model: HumanNet-PI (a protein-protein physical interaction network), HumanNet-FN (a functional gene network), and HumanNet-XC (a functional network extended by co-citation). Users can select a suitable tier of HumanNet for their study purpose. We showed that on disease gene predictions, HumanNet v3 outperforms both the previous HumanNet version and other integrated human gene networks. Furthermore, we demonstrated that HumanNet provides a feasible approach for selecting host genes likely to be associated with COVID-19.

Contemporary Approaches to the Discovery and Development of Broad-Spectrum Natural Product Prototypes for the Control of Coronaviruses.

The pressing need for SARS-CoV-2 controls has led to a reassessment of strategies to identify and develop natural product inhibitors of zoonotic, highly virulent, and rapidly emerging viruses. This review article addresses how contemporary approaches involving computational chemistry, natural product (NP) and protein databases, and mass spectrometry (MS) derived target-ligand interaction analysis can be utilized to expedite the interrogation of NP structures while minimizing the time and expense of extraction, purification, and screening in BioSafety Laboratories (BSL)3 laboratories. The unparalleled structural diversity and complexity of NPs is an extraordinary resource for the discovery and development of broad-spectrum inhibitors of viral genera, including Betacoronavirus, which contains MERS, SARS, SARS-CoV-2, and the common cold. There are two key technological advances that have created unique opportunities for the identification of NP prototypes with greater efficiency: (1) the application of structural databases for NPs and target proteins and (2) the application of modern MS techniques to assess protein-ligand interactions directly from NP extracts. These approaches, developed over years, now allow for the identification and isolation of unique antiviral ligands without the immediate need for BSL3 facilities. Overall, the goal is to improve the success rate of NP-based screening by focusing resources on source materials with a higher likelihood of success, while simultaneously providing opportunities for the discovery of novel ligands to selectively target proteins involved in viral infection.

COVID19 Disease Map, a computational knowledge repository of virus-host interaction mechanisms.

We need to effectively combine the knowledge from surging literature with complex datasets to propose mechanistic models of SARS-CoV-2 infection, improving data interpretation and predicting key targets of intervention. Here, we describe a large-scale community effort to build an open access, interoperable and computable repository of COVID-19 molecular mechanisms. The COVID-19 Disease Map (C19DMap) is a graphical, interactive representation of disease-relevant molecular mechanisms linking many knowledge sources. Notably, it is a computational resource for graph-based analyses and disease modelling. To this end, we established a framework of tools, platforms and guidelines necessary for a multifaceted community of biocurators, domain experts, bioinformaticians and computational biologists. The diagrams of the C19DMap, curated from the literature, are integrated with relevant interaction and text mining databases. We demonstrate the application of network analysis and modelling approaches by concrete examples to highlight new testable hypotheses. This framework helps to find signatures of SARS-CoV-2 predisposition, treatment response or prioritisation of drug candidates. Such an approach may help deal with new waves of COVID-19 or similar pandemics in the long-term perspective.

Drug repurposing improves disease targeting 11-fold and can be augmented by network module targeting, applied to COVID-19.

This analysis presents a systematic evaluation of the extent of therapeutic opportunities that can be obtained from drug repurposing by connecting drug targets with disease genes. When using FDA-approved indications as a reference level we found that drug repurposing can offer an average of an 11-fold increase in disease coverage, with the maximum number of diseases covered per drug being increased from 134 to 167 after extending the drug targets with their high confidence first neighbors. Additionally, by network analysis to connect drugs to disease modules we found that drugs on average target 4 disease modules, yet the similarity between disease modules targeted by the same drug is generally low and the maximum number of disease modules targeted per drug increases from 158 to 229 when drug targets are neighbor-extended. Moreover, our results highlight that drug repurposing is more dependent on target proteins being shared between diseases than on polypharmacological properties of drugs. We apply our drug repurposing and network module analysis to COVID-19 and show that Fostamatinib is the drug with the highest module coverage.

Identification of Novel COVID-19 Biomarkers by Multiple Feature Selection Strategies.

Coronavirus disease 2019 (COVID-19) arising from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in a global pandemic since its first report in December 2019. So far, SARS-CoV-2 nucleic acid detection has been deemed as the golden standard of COVID-19 diagnosis. However, this detection method often leads to false negatives, thus triggering missed COVID-19 diagnosis. Therefore, it is urgent to find new biomarkers to increase the accuracy of COVID-19 diagnosis. To explore new biomarkers of COVID-19 in this study, expression profiles were firstly accessed from the GEO database. On this basis, 500 feature genes were screened by the minimum-redundancy maximum-relevancy (mRMR) feature selection method. Afterwards, the incremental feature selection (IFS) method was used to choose a classifier with the best performance from different feature gene-based support vector machine (SVM) classifiers. The corresponding 66 feature genes were set as the optimal feature genes. Lastly, the optimal feature genes were subjected to GO functional enrichment analysis, principal component analysis (PCA), and protein-protein interaction (PPI) network analysis. All in all, it was posited that the 66 feature genes could effectively classify positive and negative COVID-19 and work as new biomarkers of the disease.

SARS-CoV-2 structural coverage map reveals viral protein assembly, mimicry, and hijacking mechanisms.

We modeled 3D structures of all SARS-CoV-2 proteins, generating 2,060 models that span 69% of the viral proteome and provide details not available elsewhere. We found that ˜6% of the proteome mimicked human proteins, while ˜7% was implicated in hijacking mechanisms that reverse post-translational modifications, block host translation, and disable host defenses; a further ˜29% self-assembled into heteromeric states that provided insight into how the viral replication and translation complex forms. To make these 3D models more accessible, we devised a structural coverage map, a novel visualization method to show what is-and is not-known about the 3D structure of the viral proteome. We integrated the coverage map into an accompanying online resource (https://aquaria.ws/covid) that can be used to find and explore models corresponding to the 79 structural states identified in this work. The resulting Aquaria-COVID resource helps scientists use emerging structural data to understand the mechanisms underlying coronavirus infection and draws attention to the 31% of the viral proteome that remains structurally unknown or dark.

Mass spectrometry-based proteomic platforms for better understanding of SARS-CoV-2 induced pathogenesis and potential diagnostic approaches.

While protein-protein interaction is the first step of the SARS-CoV-2 infection, recent comparative proteomic profiling enabled the identification of over 11,000 protein dynamics, thus providing a comprehensive reflection of the molecular mechanisms underlying the cellular system in response to viral infection. Here we summarize and rationalize the results obtained by various mass spectrometry (MS)-based proteomic approaches applied to the functional characterization of proteins and pathways associated with SARS-CoV-2-mediated infections in humans. Comparative analysis of cell-lines versus tissue samples indicates that our knowledge in proteome profile alternation in response to SARS-CoV-2 infection is still incomplete and the tissue-specific response to SARS-CoV-2 infection can probably not be recapitulated efficiently by in vitro experiments. However, regardless of the viral infection period, sample types, and experimental strategies, a thorough cross-comparison of the recently published proteome, phosphoproteome, and interactome datasets led to the identification of a common set of proteins and kinases associated with PI3K-Akt, EGFR, MAPK, Rap1, and AMPK signaling pathways. Ephrin receptor A2 (EPHA2) was identified by 11 studies including all proteomic platforms, suggesting it as a potential future target for SARS-CoV-2 infection mechanisms and the development of new therapeutic strategies. We further discuss the potentials of future proteomics strategies for identifying prognostic SARS-CoV-2 responsive age-, gender-dependent, tissue-specific protein targets.

Structural Insight into the Binding of Cyanovirin-N with the Spike Glycoprotein, Mpro and PLpro of SARS-CoV-2: Protein-Protein Interactions, Dynamics Simulations and Free Energy Calculations.

The emergence of COVID-19 continues to pose severe threats to global public health. The pandemic has infected over 171 million people and claimed more than 3.5 million lives to date. We investigated the binding potential of antiviral cyanobacterial proteins including cyanovirin-N, scytovirin and phycocyanin with fundamental proteins involved in attachment and replication of SARS-CoV-2. Cyanovirin-N displayed the highest binding energy scores (-16.8 ± 0.02 kcal/mol, -12.3 ± 0.03 kcal/mol and -13.4 ± 0.02 kcal/mol, respectively) with the spike protein, the main protease (Mpro) and the papainlike protease (PLpro) of SARS-CoV-2. Cyanovirin-N was observed to interact with the crucial residues involved in the attachment of the human ACE2 receptor. Analysis of the binding affinities calculated employing the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) approach revealed that all forms of energy, except the polar solvation energy, favourably contributed to the interactions of cyanovirin-N with the viral proteins. With particular emphasis on cyanovirin-N, the current work presents evidence for the potential inhibition of SARS-CoV-2 by cyanobacterial proteins, and offers the opportunity for in vitro and in vivo experiments to deploy the cyanobacterial proteins as valuable therapeutics against COVID-19.

Network pharmacology study on the mechanism of the Chinese medicine Radix Isatidis (Banlangen) for COVID-19.

ABSTRACT: Radix Isatidis (Banlangen) is a well-known traditional Chinese medicine for the treatment of different diseases and prevention of many body disorders. Besides, it also plays a pivotal role in novel coronavirus pneumonia, coronavirus disease 2019 (COVID-19). However, few researchers know its active ingredients and mechanism of action for COVID-19. To find whether Banlangen has a pharmacological effect on COVID-19. In this research, we systematically analyze Banlangen and COVID-19 through network pharmacology technology. A total of 33 active ingredients in Banlangen, 92 targets of the active ingredients, and 259 appropriate targets of COVID-19 were obtained, with 11 common targets. The analysis of the biological process of gene ontology and the enrichment of Kyoto Encyclopedia of Genes and Genomes signaling pathway suggests that Banlangen participated in the biological processes of protein phosphatase binding, tetrapyrrole binding, the apoptotic process involving cysteine-type endopeptidase activity, etc. The COVID-19 may be treated by regulating advanced glycation end products/a receptor for advanced glycation end products signaling pathway, interleukin-17 signaling pathway, tumor necrosis factor signaling pathway, sphingolipid signaling pathway, and p53 signaling pathway. Banlangen has a potential pharmacological effect on COVID-19, which has the value of further exploration in the following experiment and clinical application.

Clinical significance and molecular mechanism of angiotensin-converting enzyme 2 in hepatocellular carcinoma tissues.

During the pandemic of the coronavirus disease 2019, there exist quite a few studies on angiotensin-converting enzyme 2 (ACE2) and SARS-CoV-2 infection, while little is known about ACE2 in hepatocellular carcinoma (HCC). The detailed mechanism among ACE2 and HCC still remains unclear, which needs to be further investigated. In the current study with a total of 6,926 samples, ACE2 expression was downregulated in HCC compared with non-HCC samples (standardized mean difference = -0.41). With the area under the curve of summary receiver operating characteristic = 0.82, ACE2 expression showed a better ability to differentiate HCC from non-HCC. The mRNA expression of ACE2 was related to the age, alpha-fetoprotein levels and cirrhosis of HCC patients, and it was identified as a protected factor for HCC patients via Kaplan-Meier survival, Cox regression analyses. The potential molecular mechanism of ACE2 may be relevant to catabolic and cell division. In all, decreasing ACE2 expression can be seen in HCC, and its protective role for HCC patients and underlying mechanisms were explored in the study.