Identifying potential drug targets and candidate drugs for COVID-19: biological networks and structural modeling approaches

Using genetics and proteomics data to identify proteins causally related to COVID-19, healthspan and lifespan: a Mendelian randomization study

BACKGROUND

COVID-19 pandemic poses a heavy burden on public health and accounts for substantial mortality and morbidity. Proteins are building blocks of life, but specific proteins causally related to COVID-19, healthspan and lifespan have not been systematically examined.

METHODS

We conducted a Mendelian randomization study to assess the effects of 1,361 plasma proteins on COVID-19, healthspan and lifespan, using large GWAS of severe COVID-19 (up to 13,769 cases and 1,072,442 controls), COVID-19 hospitalization (32,519 cases and 2,062,805 controls) and SARS-COV2 infection (122,616 cases and 2,475,240 controls), healthspan (n = 300,477) and parental lifespan (~0.8 million of European ancestry).

RESULTS

We identified 35, 43, and 63 proteins for severe COVID, COVID-19 hospitalization, and SARS-COV2 infection, and 4, 32, and 19 proteins for healthspan, father's attained age, and mother's attained age. In addition to some proteins reported previously, such as SFTPD related to severe COVID-19, we identified novel proteins involved in inflammation and immunity (such as ICAM-2 and ICAM-5 which affect COVID-19 risk, CXCL9, HLA-DRA and LILRB4 for healthspan and lifespan), apoptosis (such as FGFR2 and ERBB4 which affect COVID-19 risk and FOXO3 which affect lifespan) and metabolism (such as PCSK9 which lowers lifespan). We found 2, 2 and 3 proteins shared between COVID-19 and healthspan/lifespan, such as CXADR and LEFTY2, shared between severe COVID-19 and healthspan/lifespan. Three proteins affecting COVID-19 and seven proteins affecting healthspan/lifespan are targeted by existing drugs.

CONCLUSIONS

Our study provided novel insights into protein targets affecting COVID-19, healthspan and lifespan, with implications for developing new treatment and drug repurposing.

Single Nucleus RNA Sequencing of Remnant Kidney Biopsies and Urine Cell RNA Sequencing Reveal Cell Specific Markers of Covid-19 Acute Kidney Injury

Acute kidney injury (AKI) in COVID-19 patients is associated with high mortality and morbidity. Critically ill COVID-19 patients are at twice the risk of in-hospital mortality compared to non-COVID AKI patients. We know little about the cell-specific mechanism in the kidney that contributes to worse clinical outcomes in these patients. New generation single cell technologies have the potential to provide insights into physiological states and molecular mechanisms in COVID-AKI. One of the key limitations is that these patients are severely ill posing significant risks in procuring additional biopsy tissue. We recently generated single nucleus RNA-sequencing data using COVID-AKI patient biopsy tissue as part of the human kidney atlas. Here we describe this approach in detail and report deeper comparative analysis of snRNAseq of 4 COVID-AKI, 4 reference, and 6 non-COVID-AKI biopsies. We also generated and analyzed urine transcriptomics data to find overlapping COVID-AKI-enriched genes and their corresponding cell types in the kidney from snRNA-seq data. We identified all major and minor cell types and states by using by using less than a few cubic millimeters of leftover tissue after pathological workup in our approach. Differential expression analysis of COVID-AKI biopsies showed pathways enriched in viral response, WNT signaling, kidney development, and cytokines in several nephron epithelial cells. COVID-AKI profiles showed a much higher proportion of altered TAL cells than non-COVID AKI and the reference samples. In addition to kidney injury and fibrosis markers indicating robust remodeling we found that, 17 genes overlap between urine cell COVID-AKI transcriptome and the snRNA-seq data from COVID-AKI biopsies. A key feature was that several of the distal nephron and collecting system cell types express these markers. Some of these markers have been previously observed in COVID-19 studies suggesting a common mechanism of injury and potentially the kidney as one of the sources of soluble factors with a potential role in disease progression.

Translational Statement

The manuscript describes innovation, application and discovery that impact clinical care in kidney disease. First, the approach to maximize use of remnant frozen clinical biopsies to inform on clinically relevant molecular features can augment existing pathological workflow for any frozen tissue without much change in the protocol. Second, this approach is transformational in medical crises such as pandemics where mechanistic insights are needed to evaluate organ injury, targets for drug therapy and diagnostic and prognostic markers. Third, the cell type specific and soluble markers identified and validated can be used for diagnoses or prognoses in AKI due to different etiologies and in multiorgan injury.

Integration of Omics Data and Network Models to Unveil Negative Aspects of SARS-CoV-2, from Pathogenic Mechanisms to Drug Repurposing

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) caused the COVID-19 health emergency, affecting and killing millions of people worldwide. Following SARS-CoV-2 infection, COVID-19 patients show a spectrum of symptoms ranging from asymptomatic to very severe manifestations. In particular, bronchial and pulmonary cells, involved at the initial stage, trigger a hyper-inflammation phase, damaging a wide range of organs, including the heart, brain, liver, intestine and kidney. Due to the urgent need for solutions to limit the virus' spread, most efforts were initially devoted to mapping outbreak trajectories and variant emergence, as well as to the rapid search for effective therapeutic strategies. Samples collected from hospitalized or dead COVID-19 patients from the early stages of pandemic have been analyzed over time, and to date they still represent an invaluable source of information to shed light on the molecular mechanisms underlying the organ/tissue damage, the knowledge of which could offer new opportunities for diagnostics and therapeutic designs. For these purposes, in combination with clinical data, omics profiles and network models play a key role providing a holistic view of the pathways, processes and functions most affected by viral infection. In fact, in addition to epidemiological purposes, networks are being increasingly adopted for the integration of multiomics data, and recently their use has expanded to the identification of drug targets or the repositioning of existing drugs. These topics will be covered here by exploring the landscape of SARS-CoV-2 survey-based studies using systems biology approaches derived from omics data, paying particular attention to those that have considered samples of human origin.

Role of E3 ubiquitin ligases and deubiquitinating enzymes in SARS-CoV-2 infection

Ever since its emergence in 2019, COVID-19 has rapidly disseminated worldwide, engendering a pervasive pandemic that has profoundly impacted healthcare systems and the socio-economic milieu. A plethora of studies has been conducted targeting its pathogenic virus, SARS-CoV-2, to find ways to combat COVID-19. The ubiquitin-proteasome system (UPS) is widely recognized as a crucial mechanism that regulates human biological activities by maintaining protein homeostasis. Within the UPS, the ubiquitination and deubiquitination, two reversible modifications, of substrate proteins have been extensively studied and implicated in the pathogenesis of SARS-CoV-2. The regulation of E3 ubiquitin ligases and DUBs(Deubiquitinating enzymes), which are key enzymes involved in the two modification processes, determines the fate of substrate proteins. Proteins associated with the pathogenesis of SARS-CoV-2 may be retained, degraded, or even activated, thus affecting the ultimate outcome of the confrontation between SARS-CoV-2 and the host. In other words, the clash between SARS-CoV-2 and the host can be viewed as a battle for dominance over E3 ubiquitin ligases and DUBs, from the standpoint of ubiquitin modification regulation. This review primarily aims to clarify the mechanisms by which the virus utilizes host E3 ubiquitin ligases and DUBs, along with its own viral proteins that have similar enzyme activities, to facilitate invasion, replication, escape, and inflammation. We believe that gaining a better understanding of the role of E3 ubiquitin ligases and DUBs in COVID-19 can offer novel and valuable insights for developing antiviral therapies.

Prospects of Novel and Repurposed Immunomodulatory Drugs against Acute Respiratory Distress Syndrome (ARDS) Associated with COVID-19 Disease

Acute respiratory distress syndrome (ARDS) is intricately linked with SARS-CoV-2-associated disease severity and mortality, especially in patients with co-morbidities. Lung tissue injury caused as a consequence of ARDS leads to fluid build-up in the alveolar sacs, which in turn affects oxygen supply from the capillaries. ARDS is a result of a hyperinflammatory, non-specific local immune response (cytokine storm), which is aggravated as the virus evades and meddles with protective anti-viral innate immune responses. Treatment and management of ARDS remain a major challenge, first, because the condition develops as the virus keeps replicating and, therefore, immunomodulatory drugs are required to be used with caution. Second, the hyperinflammatory responses observed during ARDS are quite heterogeneous and dependent on the stage of the disease and the clinical history of the patients. In this review, we present different anti-rheumatic drugs, natural compounds, monoclonal antibodies, and RNA therapeutics and discuss their application in the management of ARDS. We also discuss on the suitability of each of these drug classes at different stages of the disease. In the last section, we discuss the potential applications of advanced computational approaches in identifying reliable drug targets and in screening out credible lead compounds against ARDS.

Identification of host genomic biomarkers from multiple transcriptomics datasets for diagnosis and therapies of SARS-CoV-2 infections

The pandemic of COVID-19 is a severe threat to human life and the global economy. Despite the success of vaccination efforts in reducing the spread of the virus, the situation remains largely uncontrolled due to the random mutation in the RNA sequence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which demands different variants of effective drugs. Disease-causing gene-mediated proteins are usually used as receptors to explore effective drug molecules. In this study, we analyzed two different RNA-Seq and one microarray gene expression profile datasets by integrating EdgeR, LIMMA, weighted gene co-expression network and robust rank aggregation approaches, which revealed SARS-CoV-2 infection causing eight hub-genes (HubGs) including HubGs; REL, AURKA, AURKB, FBXL3, OAS1, STAT4, MMP2 and IL6 as the host genomic biomarkers. Gene Ontology and pathway enrichment analyses of HubGs significantly enriched some crucial biological processes, molecular functions, cellular components and signaling pathways that are associated with the mechanisms of SARS-CoV-2 infections. Regulatory network analysis identified top-ranked 5 TFs (SRF, PBX1, MEIS1, ESR1 and MYC) and 5 miRNAs (hsa-miR-106b-5p, hsa-miR-20b-5p, hsa-miR-93-5p, hsa-miR-106a-5p and hsa-miR-20a-5p) as the key transcriptional and post-transcriptional regulators of HubGs. Then, we conducted a molecular docking analysis to determine potential drug candidates that could interact with HubGs-mediated receptors. This analysis resulted in the identification of top-ranked ten drug agents, including Nilotinib, Tegobuvir, Digoxin, Proscillaridin, Olysio, Simeprevir, Hesperidin, Oleanolic Acid, Naltrindole and Danoprevir. Finally, we investigated the binding stability of the top-ranked three drug molecules Nilotinib, Tegobuvir and Proscillaridin with the three top-ranked proposed receptors (AURKA, AURKB, OAS1) by using 100 ns MD-based MM-PBSA simulations and observed their stable performance. Therefore, the findings of this study might be useful resources for diagnosis and therapies of SARS-CoV-2 infections.

The pathogenetic influence of smoking on SARS-CoV-2 infection: Integrative transcriptome and regulomics analysis of lung epithelial cells

Corona virus disease (COVID-19) has been emerged as pandemic infectious disease. The recent epidemiological data suggest that the smokers are more vulnerable to infection with COVID-19; however, the influence of smoking (SMK) on the COVID-19 infected patients and the mortality is not known yet. In this study, we aimed to discern the influence of SMK on COVID-19 infected patients utilizing the transcriptomics data of COVID-19 infected lung epithelial cells and transcriptomics data smoking matched with controls from lung epithelial cells. The bioinformatics based analysis revealed the molecular insights into the level of transcriptional changes and pathways which are important to identify the impact of smoking on COVID-19 infection and prevalence. We compared differentially expressed genes (DEGs) between COVID-19 and SMK and 59 DEGs were identified as consistently dysregulated at transcriptomics levels. The correlation network analyses were constructed for these common genes using WGCNA R package to see the relationship among these genes. Integration of DEGs with network analysis (protein-protein interaction) showed the presence of 9 hub proteins as key so called "candidate hub proteins" overlapped between COVID-19 patients and SMK. The Gene Ontology and pathways analysis demonstrated the enrichment of inflammatory pathway such as IL-17 signaling pathway, Interleukin-6 signaling, TNF signaling pathway and MAPK1/MAPK3 signaling pathways that might be the therapeutic targets in COVID-19 for smoking persons. The identified genes, pathways, hubs genes, and their regulators might be considered for establishment of key genes and drug targets for SMK and COVID-19.

Identification of host transcriptome-guided repurposable drugs for SARS-CoV-1 infections and their validation with SARS-CoV-2 infections by using the integrated bioinformatics approaches

Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) is one of the most severe global pandemic due to its high pathogenicity and death rate starting from the end of 2019. Though there are some vaccines available against SAER-CoV-2 infections, we are worried about their effectiveness, due to its unstable sequence patterns. Therefore, beside vaccines, globally effective supporting drugs are also required for the treatment against SARS-CoV-2 infection. To explore commonly effective repurposable drugs for the treatment against different variants of coronavirus infections, in this article, an attempt was made to explore host genomic biomarkers guided repurposable drugs for SARS-CoV-1 infections and their validation with SARS-CoV-2 infections by using the integrated bioinformatics approaches. At first, we identified 138 differentially expressed genes (DEGs) between SARS-CoV-1 infected and control samples by analyzing high throughput gene-expression profiles to select drug target key receptors. Then we identified top-ranked 11 key DEGs (SMAD4, GSK3B, SIRT1, ATM, RIPK1, PRKACB, MED17, CCT2, BIRC3, ETS1 and TXN) as hub genes (HubGs) by protein-protein interaction (PPI) network analysis of DEGs highlighting their functions, pathways, regulators and linkage with other disease risks that may influence SARS-CoV-1 infections. The DEGs-set enrichment analysis significantly detected some crucial biological processes (immune response, regulation of angiogenesis, apoptotic process, cytokine production and programmed cell death, response to hypoxia and oxidative stress), molecular functions (transcription factor binding and oxidoreductase activity) and pathways (transcriptional mis-regulation in cancer, pathways in cancer, chemokine signaling pathway) that are associated with SARS-CoV-1 infections as well as SARS-CoV-2 infections by involving HubGs. The gene regulatory network (GRN) analysis detected some transcription factors (FOXC1, GATA2, YY1, FOXL1, TP53 and SRF) and micro-RNAs (hsa-mir-92a-3p, hsa-mir-155-5p, hsa-mir-106b-5p, hsa-mir-34a-5p and hsa-mir-19b-3p) as the key transcriptional and post- transcriptional regulators of HubGs, respectively. We also detected some chemicals (Valproic Acid, Cyclosporine, Copper Sulfate and arsenic trioxide) that may regulates HubGs. The disease-HubGs interaction analysis showed that our predicted HubGs are also associated with several other diseases including different types of lung diseases. Then we considered 11 HubGs mediated proteins and their regulatory 6 key TFs proteins as the drug target proteins (receptors) and performed their docking analysis with the SARS-CoV-2 3CL protease-guided top listed 90 anti-viral drugs out of 3410. We found Rapamycin, Tacrolimus, Torin-2, Radotinib, Danoprevir, Ivermectin and Daclatasvir as the top-ranked 7 candidate-drugs with respect to our proposed target proteins for the treatment against SARS-CoV-1 infections. Then, we validated these 7 candidate-drugs against the already published top-ranked 11 target proteins associated with SARS-CoV-2 infections by molecular docking simulation and found their significant binding affinity scores with our proposed candidate-drugs. Finally, we validated all of our findings by the literature review. Therefore, the proposed candidate-drugs might play a vital role for the treatment against different variants of SARS-CoV-2 infections with comorbidities, since the proposed HubGs are also associated with several comorbidities.

Computational identification of host genomic biomarkers highlighting their functions, pathways and regulators that influence SARS-CoV-2 infections and drug repurposing

The pandemic threat of COVID-19 has severely destroyed human life as well as the economy around the world. Although, the vaccination has reduced the outspread, but people are still suffering due to the unstable RNA sequence patterns of SARS-CoV-2 which demands supplementary drugs. To explore novel drug target proteins, in this study, a transcriptomics RNA-Seq data generated from SARS-CoV-2 infection and control samples were analyzed. We identified 109 differentially expressed genes (DEGs) that were utilized to identify 10 hub-genes/proteins (TLR2, USP53, GUCY1A2, SNRPD2, NEDD9, IGF2, CXCL2, KLF6, PAG1 and ZFP36) by the protein–protein interaction (PPI) network analysis. The GO functional and KEGG pathway enrichment analyses of hub-DEGs revealed some important functions and signaling pathways that are significantly associated with SARS-CoV-2 infections. The interaction network analysis identified 5 TFs proteins and 6 miRNAs as the key regulators of hub-DEGs. Considering 10 hub-proteins and 5 key TFs-proteins as drug target receptors, we performed their docking analysis with the SARS-CoV-2 3CL protease-guided top listed 90 FDA approved drugs. We found Torin-2, Rapamycin, Radotinib, Ivermectin, Thiostrepton, Tacrolimus and Daclatasvir as the top ranked seven candidate drugs. We investigated their resistance performance against the already published COVID-19 causing top-ranked 11 independent and 8 protonated receptor proteins by molecular docking analysis and found their strong binding affinities, which indicates that the proposed drugs are effective against the state-of-the-arts alternatives independent receptor proteins also. Finally, we investigated the stability of top three drugs (Torin-2, Rapamycin and Radotinib) by using 100 ns MD-based MM-PBSA simulations with the two top-ranked proposed receptors (TLR2, USP53) and independent receptors (IRF7, STAT1), and observed their stable performance. Therefore, the proposed drugs might play a vital role for the treatment against different variants of SARS-CoV-2 infections.

Structure-Based Virtual Screening Reveals Ibrutinib and Zanubrutinib as Potential Repurposed Drugs against COVID-19

Coronavirus disease (COVID)-19 is the leading global health threat to date caused by a severe acute respiratory syndrome coronavirus (SARS-CoV-2). Recent clinical trials reported that the use of Bruton's tyrosine kinase (BTK) inhibitors to treat COVID-19 patients could reduce dyspnea and hypoxia, thromboinflammation, hypercoagulability and improve oxygenation. However, the mechanism of action remains unclear. Thus, this study employs structure-based virtual screening (SBVS) to repurpose BTK inhibitors acalabrutinib, dasatinib, evobrutinib, fostamatinib, ibrutinib, inositol 1,3,4,5-tetrakisphosphate, spebrutinib, XL418 and zanubrutinib against SARS-CoV-2. Molecular docking is conducted with BTK inhibitors against structural and nonstructural proteins of SARS-CoV-2 and host targets (ACE2, TMPRSS2 and BTK). Molecular mechanics-generalized Born surface area (MM/GBSA) calculations and molecular dynamics (MD) simulations are then carried out on the selected complexes with high binding energy. Ibrutinib and zanubrutinib are found to be the most potent of the drugs screened based on the results of computational studies. Results further show that ibrutinib and zanubrutinib could exploit different mechanisms at the viral entry and replication stage and could be repurposed as potential inhibitors of SARS-CoV-2 pathogenesis.