Introduction

The pandemic caused by Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has put immense pressure on the health care system and resulted in the indefinite future. Its genome consists of positive single-stranded RNA (genus betacoronavirus) comprising 14 open reading frames (Orfs) which encodes 16 non-structural proteins (nsP1-16), four structural (S, E, M and N) and nine accessory proteins (orf3a, orf3b, orf6, orf7a, orf7b, orf8, orf9b, orf9c and orf10). SARS-CoV-2 may confound seasons, and its transmission is not dependent on specific temperature conditions [21]. Dengue virus (DENV) and Chikungunya virus (CHIKV) infections become prevalent during monsoon season in tropical and subtropical countries. Mosquito species Aedes aegypti is the major vector for the transmission of DENV as well as CHIKV [45, 55]. There are reports that the severity of the disease is much more than the mono-infection during their co-infection [29, 45, 55]. SARS-CoV-2 and DENV infections have also been reported to share clinical and laboratory features [72]. Accumulated data from Singapore [72], India [3, 25], and Mayotte (Island in the Indian Ocean) [15] suggested that the patients have a probability of being co-infected with SARS-CoV-2 and DENV. A report from Thailand reported that the SARS-CoV-2 patients exhibited symptoms like skin rash with petechiae and low platelet count initially, which are common clinical manifestations during DENV infection [27]. Though the co-infection of CHIKV and SARS-CoV-2 is not reported yet, the probability in the future exists. Furthermore, the co-incidence of DENV, CHIKV, and SARS-CoV-2 at the same time could also be predicted with more severity in DENV/CHIKV predominant areas. SARS-CoV-2 is one of the greatest challenges of recent times, and its interaction with host proteins during co-infection needs to be deciphered.

Both DENV and CHIKV are single-stranded RNA viruses that belong to genus flavivirus and alphavirus, respectively [29]. The genome of CHIKV encodes four nsPs (nsP1-4) and five structural proteins (E1, E2, E3, 6 K and C) while DENV encodes seven nsPs (NS1, NS2A, NS2B, NS3, NS4A, NS4B, and NS5) and three structural proteins (C, prM, and E) [38, 56].

SARS-CoV-2, DENV and CHIKV infections have been strongly associated with marked increase in inflammatory cytokines. Higher concentrations of pro-inflammatory cytokines have been reported in the severely critically ill patients of coronavirus disease (Covid-19) [7, 74]. The cytokine storm was reported to be responsible for the tissue damage (heart, liver, and kidney) and multiple organ failure [14]. Studies have demonstrated the association of cardiovascular manifestations during CHIKV and DENV infections due to enhanced up-regulation of cytokines [1, 57]. Thus, the cytokine storm and the multi-organ dysfunction have been observed during the mono-infections of CHIKV, DENV and SARS-CoV-2.

The characterization of interactions between viral and host proteins help to understand the molecular mechanism of infection. Several studies have reported the interactions of human proteins with DENV and CHIKV during mono-infection [14, 33, 43]. However, the data about SARS-CoV-2 infection is limited [23, 55]. The study of viral-host protein interactome during co-infection will help to understand the residing mechanism adopted by the virus for its replication and pathogenesis.

In the present study, the viral-host protein–protein interactions were analysed during the co-infection of DENV, CHIKV, and SARS-CoV-2 using in silico approach. The cellular interactors during mono-infections were obtained from the published studies. The interactome network of each disease was constructed and merged to identify the cellular interactors during co-infection. The biological and molecular functions of identified proteins, as well as the associated pathways, were analysed. Moreover, the drug target proteins were predicted, which could be targeted to liquidate DENV, CHIKV, and SARS-CoV-2 co-infection.

Materials and methods

Collection of data

The data of human host interactors of DENV, CHIKV, and SARS-CoV-2 mono-infections were collected from the published studies [14, 21, 43]. Doolittle et al. and Rana et al. identified the host interactors using computationally approaches for DENV and CHIKV, respectively. They used structural similarity approach with the human proteins, followed by utilization of databases, Biological General Repository for Interaction Datasets (BIOGRID), and Human Protein Reference Database (HPRD), for identifying human interactors [14, 43]. Gordon et al. utilized the affinity purification mass spectrometry (AP-MS) approach to identify the cellular interactors of SARS-CoV-2 [21].

Identification of cellular interactors during co-infection

The network of accumulated host protein interactions data during mono-infections was constructed using the STRING online tool [60] and then imported in Cytoscape [52]. The in-built merge tool in Cytoscape was utilized to merge the individual three networks. The merging of data yielded a group of human proteins that were common in all three individual proteins networks.

Functional analysis

The biological and functional analyses of identified proteins were studied to understand the cellular response during co-infection. The Protein Analysis Through Evolutionary Relationships (PANTHER) tool, Kyoto Encyclopedia of Genes and Genomes (KEGG) database, and the Universal Protein Resource (UniProt) Knowledgebase (UniProtKB) were exploited to study the functions and the linked pathways [28, 62, 63]. The PANTHER tool facilitated the molecular and biological functions as well as processes of uncharacterized genes/proteins based on the evolutionary relationships [62]. UniProtKB is a comprehensive resource for protein sequence and annotation data [75] and KEGG database assisted in deciphering the functional behaviours of gene, biological pathways, and diseases [28].

Identification of drug targets protein

The identified human proteins were further analysed to be defined as drug targets. DrugBank [69], and Therapeutic Target Database (TTD) [67] were used in the study for the prediction of the druggable proteins. The DrugBank database is a unique bioinformatics and cheminformatics resource that combines detailed drug data with comprehensive drug target information. The DrugBank (version 5.1.6) comprises of 13,570 drug entries including 2631 approved small molecule drugs, 1377 approved biologics (proteins, peptides, vaccines, and allergenics), 131 nutraceuticals and over 6373 experimental (discovery-phase) drugs. Additionally, 5252 non-redundant protein (i.e. drug target/enzyme/transporter/carrier) sequences are linked to these drug entries. Each entry contains more than 200 data fields with half of the information being devoted to drug/chemical data, and the other half devoted to drug target or protein data [69]. Furthermore, TTD contains the known and explored therapeutic protein and nucleic acid targets, the targeted disease, pathway information, and the corresponding drugs directed at each of these targets [67].

Results and discussion

Decode host interactors during co-infection

The data of human protein interactions with DENV, CHIKV, and SARS-CoV-2 during mono-infections were collected for further analysis [14, 21, 43]. The networks of human proteins during viral mono-infections were constructed using Cytoscape and observed to be densely interconnected. The construction of proteome network and its analysis using Cytoscape have also been utilized to study protein–protein interactions for other viruses including the Hepatitis C virus (HCV) [17], Hepatitis B virus (HBV) [54], and gynaecological cancer [31]. The human interactors during the co-infection of DENV, CHIKV, and SARS-CoV-2 were elucidated by the merged constructed networks. It was observed that the total 17 human interactors shared in all the three networks could act as cellular interacting partners for viral proteins during co-infection (Table 1). The same strategy was also exploited to identify the common nodes and the novel candidate biomarkers during HBV and HCV co-infection [54]. Taz et al. generated the protein–protein interaction network during dengue, malaria, and chikungunya disease to elucidate the common genes present in all the three diseased conditions as drug targets. This computational approach is a powerful tool to study the hub of genes, proteins, and pathways related to the diseased conditions [61].

Table 1 The identified human proteins as interactors of DENV, CHIKV, and SARS-CoV-2 during co-infection

Moreover, it was observed that most elucidated host proteins mainly interacted with the nsPs and few structural (mainly envelope) proteins. It could be predicted that the interactions involved with nsPs might influence the RNA synthesis, replication, and viral life cycle. Host proteins interacting with the envelop proteins might affect the viral entry, virion formation, as well as virulence.

Functional analysis of identified proteins

The cellular functions and the associated pathways of the identified interactors of DENV, CHIKV, and SARS-CoV-2 during co-infection were analysed using the KEGG database, Panther tool, and UniProtKB (Table 2). It was found that the identified proteins were involved in signaling pathways, RNA transport and processing, ubiquitination, vesicle trafficking, innate and adaptive immune responses. These pathways help to understand viral biology in relation to cellular and molecular mechanisms.

Table 2 The functions and pathways linked with the identified host proteins

One of the identified proteins in the study was integrin beta 1 (ITGB1), as a cell adhesion receptor is exploited by viruses to make a successful entry into the host cell. The up-regulation of ITGB1 has been reported in CHIKV, DENV and SARS coronavirus (SARS-CoV) infections [10, 11, 58]. ITGBI was also reported to be crucial to maintain the infection of flaviviruses like HCV [75] and the West Nile virus (WNV) [50]. Therefore, it could be predicted that during DENV, CHIKV, and SARS-CoV-2 co-infection, their interaction with the host’s ITGB1 might assist the viral entry and disseminate infection.

The identified interactors, Ras homolog family member A (RHOA) and ITGB1, have been associated with the phosphatidylinositol-3-kinase (PI3K)/Akt and the mammalian target of rapamycin (mTOR) signaling pathways. CHIKV has been reported to utilize this pathway for its replication [13]. Liu et al. reported the modulation of the PI3K-Akt pathway to induce apoptosis during DENV infection [30].

The common interactor identified in the study, Poly(A)-binding protein cytoplasmic 1 (PABPC1), is a member of poly(A) binding proteins (PABP) family, which plays multiple roles in regulating the mRNA stability, mRNA metabolism and RNA transport. Suzuki et al. reported that the binding of PABPC1 with the A-rich region in the 3′UTR of DENV mRNA positively regulated its translation [59]. Its interaction with the nsP2 protein of CHIKV [6] and N protein of SARS-CoV-2 was also reported [21].

Another identified interactor, Ring box protein-1 (RBX1), is a component of ubiquitin–proteasome pathways (UPP) that has been shown to regulate the protein stability and protein trafficking [19]. RBX-1 is also associated with the Wnt signaling pathway, which is critical to control cell cycle, development, cell motility, and cell death. RBX-1 was reported to be differentially expressed during CHIKV infection. Though, the upregulation of its associated pathways was observed during DENV, WNV, Japanese encephalitis virus (JEV), and Human immunodeficiency virus (HIV) infections [19, 70, 76].

The nuclear pore complexes (NPC) present in the nuclear envelope are composed of nucleoporins (NUPs) and act as multiprotein channels. The NPC acts as a selective barrier for bidirectional transport of RNA, viral proteins, and molecules from the nucleus to the cytoplasm. The viruses could hijack this trafficking by remodelling the nuclear membrane, redirecting the host machinery for its own replication and combat anti-viral response [46]. In the current study, Nup214, NUP62, and NUP98 were identified as shared host interactors for the SARS-CoV-2, DENV, and CHIKV viruses. It is known that CHIKV, DENV, and SARS-CoV-2 primarily replicate in host cytoplasm, but intriguingly, some of the viral proteins get localised into the nucleus, which is essential for maintaining their infection [41, 44, 71]. NUP98 was reported to be accumulated in the nucleus with the nsP3 protein of CHIKV, but the relation is still unexplored [44]. The cleavage and degradation of NUP98 and NUP62 by the serine protease activity of NS2B-NS3 were observed during DENV infection. This degradation significantly inhibited the mRNA export through the NPC [8]. Nucleoporins (NUP98, NUP214, and NUP62) have been reported as interactors of SARS-CoV-2 and interfered in nuclear export, but the detailed mechanism is unknown [21]. Moreover, NUP98 and NUP214 were reported to be crucial in many other viral infections [24, 39, 48, 51].

Up-frame shift 1 (UPF1) is an ATP-dependent RNA and DNA helicase involved in nonsense-mediated mRNA decay (NMD), part of the eukaryotic RNA surveillance pathway [18]. UPF1 interacts with the initiation factors during translation and releases the ribosome to degrade mRNAs containing premature stop codons. Balistreri et al. observed that depletion of UPF1 with other co-factors resulted in the increased infection of the Semliki Forest virus (SFV, alphavirus) and hence, concluded that the UPF1 was involved in the early steps of infection to reduce the synthesis of viral RNA and proteins [4]. Moreover, the N protein of SARS-CoV and Middle East respiratory syndrome CoV (MERS-CoV) has been shown to suppress the NMD pathway to protect its mRNA from degradation and ensure efficient replication [65].

Another identified interactor, a histone deacetylase-2 (HDAC2), plays a crucial role in transcriptional regulation, cell cycle progression, and developmental events [12]. HDAC2, primarily located in the nucleus, is reported to have an essential role in regulating the inflammatory response of macrophages. HDAC2 is associated with the TNF signaling as well as longevity regulating pathway. It was reported that dengue hemorrhage with a high TNF level caused the endothelial cells apoptosis resulting in vascular damage [73]. Although the molecular insight of HDAC2 is not fully understood, it was reported to be crucial for viral infections like Influenza A virus (IAV) [37] and Human alphaherpesvirus 1 (HSV-1) [66].

Identified interactor, microtubule affinity-regulating kinase 3 (MAPK3) accounts for the specific phosphorylation of microtubule-associated proteins and linked with the intracellular signal transduction pathways. F2RL1 (coagulation factor II (thrombin) receptor-like 1), was found as a common interacting host protein during CHIKV, DENV, and SARS-CoV-2 co-infection in the current study. F2RL1 is also known as protease-activated receptor 2 (PAR2) and belongs to the family of G protein-coupled receptors. It accounts for the inflammatory process and modulates the immune response mediated by toll-like receptors (TLR) signaling and macrophages during viral infection [2]. Schanoski et al. reported the production of Granzyme A, which was correlated to PAR-1 and PAR-2, responsible for the inflammation and swelling induced during CHIKV, but the exact mechanism is not known [49]. Moreover, the activation of PAR2 was reported to modulate HIV and IAV infections [2].

Another identified protein, ADP Ribosylation Factor 6 (ARF6), is a membrane trafficking protein that coordinates the plasma membrane dynamics and regulates the internalization of cargo via endocytosis [64]. GTPase-activating proteins (GAP) are generally required for the regulation (activation and inactivation) of Arf6. In the case of DENV, the type I interferons (IFN-I) up-regulated the expression of ADAP2 (ArfGAP domain-containing protein 2), a GTPase-activating protein for Arf6, that suppressed its infection primarily at the stage of viral entry [53, 64]. Moreover, Radoshitzky et al. reported ARF6 in association with other host factors were involved in alphaviruses trafficking [42]. Furthermore, it was required for entry and early infection during HIV infection [20].

Identification of druggable proteins

The identified human interactors were further studied to predict drug target proteins using databases like DrugBank and TTD. The databases contain the necessary information regarding the drug targets, approved drugs, and experimental drugs that were exploited in the current study to elucidate the possible drug targets [67, 69]. Among the identified cellular interactors, seven human proteins were identified and predicted as druggable targets. The information regarding the identified target proteins, agonist, purposed molecules/ drugs, along with their indications, are summarised in Table 3. It was analysed that ITGB1, HDAC2, and MARK3 human proteins were targeted by FDA approved small molecules/ drugs. At the same time, the other predicted druggable proteins were targeted by the small molecules, which are present either in the experimental phase or in the clinical phase.

Table 3 The summary of predicted drug target proteins

Moreover, the human protein HDAC2 was observed to be targeted by nine FDA approved drugs emphasizing it to be a significant druggable protein. Additionally, the approved molecules/drugs targeting HDAC2 were observed to be effective for the inhibition of other viral infections. Vorinostat, an HDAC inhibitor used in the treatment of lymphoma, was reported to significantly abrogate the genome replication of Human papillomaviruses via induction of apoptosis, shown as a promising therapeutic agent [5]. Antiviral activity of lovastatin has been reported against Respiratory syncytial virus (RSV) [22], DENV Serotype 2 [34], Zika virus [16], and HIV-1 [32].

Montoya et al. suggested that the administration of lovastatin as anti-HIV-1 in chronically infected patients was safe [36]. Martinez-Gutierrez et al. suggested that lovastatin delayed the progression of DENV-2 in infected mice [34]. However, Whitehorn et al. did not observe any clinical benefits while treating DENV infected adults with lovastatin [68]. Gower et al. reported that lovastatin inhibited RSV by reducing the production of cholesterol and isoprenylation of RhoA [22]. Moreover, atorvastatin, simvastatin, and pravastatin were considered conventional therapeutics for IAV-H1N1, as these approved molecules effectively reduced the cytokines profile in the H1N1 infected cells [35]. Additionally, the use of statins (class of drugs) for the treatment of COVID-19 has been predicted as they help to reduce the level of cholesterol, anti‐inflammatory, and immunomodulatory activities [40]. Another approved molecule, valproic acid targeting HDAC2, has been reported to significantly inhibit the HSV infection in an oligodendrocyte cell line [9].

Recently, in silico study revealed that the fostamatinib (FDA approved drug) has the potential to treat COVID-19. Fostamatinib, approved for chronic immune thrombocytopenia, is an inhibitor of spleen tyrosine kinase (Syk), an important signalling component in immune cells. The docking of fostamatinib with the protease of SARS-CoV-2 had the highest docking score compared to other drugs like hydroxychloroquine, remdesivir, favipiravir, and darunavir, which have the potential to treat COVID-19 [47]. Additionally, fostamatinib is in phase 2 clinical trial to treat COVID-19 patients [26] as it can control dysregulated immune system responsible for the underlying symptoms of SARS-CoV-2. In the present study, MARK3 was found to be the possible host target of fostamatinib to be explored further.

The chances of co-infection of arboviruses (DENV and CHIKV) with SARS-CoV-2 predicted during the rainy season because of the favourable breeding conditions for the mosquitoes, and the current pandemic of SARS-CoV-2 could severely affect the human health. Deciphering the viral-host interactions assist in understanding the viral biology as well as the cellular mechanism. The present study elucidates the interacting cellular proteins and druggable targets that could play a crucial role during the co-infection of DENV, CHIKV, and SARS-CoV-2. The identified proteins were linked with several cellular and signaling pathways that might be crucial for sustained viral infection. Moreover, the approved therapeutic molecules against the potential drug target proteins were predicted in the study. Eventually, these interactions could be confirmed and utilized to develop effective therapeutics.