Human adenovirus DNA polymerase is evolutionarily and functionally associated with human telomerase reverse transcriptase based on in silico molecular characterization that implicate abacavir and zidovudine

Abstract

Human adenoviruses (HAdVs) are non-enveloped, small double stranded DNA (dsDNA) viruses that cause asymptomatic infections, clinical syndromes and significant susceptibility to infections in immunocompromised people. The aim of the present study was to identify critical host proteins and HAdV hypothetical proteins that could be developed as potential host-viral targets for antiHAdV therapy. Here, the function of selected hypothetical proteins of HAdV based on phylogenetic relationship with the therapeutic targets of antiretroviral drugs of human immunodeficiency virus (HIV) was predicted computationally, and characterized the molecular dynamics and binding affinity of DNA polymerase of HAdV. Thirty-eight hypothetical proteins (HPs) of human adenovirus (HAdV) were used in this study. The results showed that HAdV DNA polymerase (P03261) is related to Human TERT (O14746) and HLA-B (P01889) genes. The protein-protein interaction of human five molecular targets (PNP, TERT, CCR5, HLA-B, and NR1I2) of ARVDs are well-coordinated/networked with CD4, AHR, FKBP4, NR3C1, HSP90AA1, and STUB1 proteins in the anti-HIV infection mechanism. The results showed that the free energy score of abacavir and zidovudine binding to HAdV DNA polymerase are −5.8 and −5.4 kcal mol^-1 respectively. Also, the control drug, cidofovir and ganciclovir have less binding affinity for DNA polymerase of HAdV when compare to that of abacavir and zidovudine. Similarity was observed in the binding of abacavir and zidovudine to HAdV DNA polymerase (ASP742, ALA743, LEU772, ARG773 and VAL776). In conclusion, combination of abacavir and zidovudine was predicted to be potential therapy for controlling HAdV infection targeting HAdV DNA polymerase.

Introduction

Human adenoviruses (HAdVs) are non-enveloped, small double stranded DNA (dsDNA) viruses that cause asymptomatic infection and clinical syndromes, including hepatitis, myocarditis, gastroenteritis, pneumonia, upper respiratory tract infection, conjunctivitis, and infections of the urinary tract (Curlin et al., 2010; Robinson et al., 2013; Naveed et al., 2017). Globally, about 5%–7% of respiratory tract infections has been attributed to HAdV in pediatrics population below age of 5 years, and that HAdVs cause significant susceptibility to infections in immunocompromised people (Hierholzer, 1992; James et al., 2007; Ghebremedhin, 2014; Scott et al., 2016). Adenoviruses can cause life-threatening infections in hematopoietic stem cell transplant recipients, including hepatitis, pneumonitis, colitis, nephritis, and encephalitis (La Rosa et al., 2001).

The adenovirus has a linear genome of about 35 kBp in size with a 55 kDa the terminal protein (TP) covalently attached to each 5’ end. The 140 kDa adenovirus-encoded DNA polymerase (Ad Pol) is required for viral DNA replication both in vitro and in vivo (Field et al., 1984). The in vitro replication of adenovirus DNA using HeLa cell, requires at three proteins from the virus (the 59 kDa DNA-binding protein (DBP), the 80 kDa precursor TP, and the 140 kDa Ad Pol); and two proteins from the host (HeLa cell) are nuclear proteins (factor I and factor II), which are both required for initiation and elongation of viral DNA synthesis (Field et al., 1984). A study has shown that nuclear factor II of HAdV possesses type I topoisomerase activity and could be substituted in the in vitro DNA replication system with purified type I topoisomerases from HeLa cells or calf thymus but not that Escherichia coli (Nagata et al., 1983).

Currently, there is no Food and Drug Administration (FDA) approved drug for the treatment of adenovirus infection, but the available regimen are the combinations of existing antiviral drugs such as cidofovir and its derivative (3-hexadecyloxy-1-propanol-cidofovir (also known as brincidofovir or CMX001)), ganciclovir, filociclovir, cytarabine, ivermectin, zalacitabine, arbidol (umifenovir), stavudine and ribavirin as well as anticancer drugs such as camptothecin, and other chemicals such as epiandrosterone (Wong and Hsu, 1990; Waye and Sing, 2010; Kumar et al., 2019; Saha and Parks, 2020; Dodge et al., 2021), RIDK34, digitoxigenin and digoxin (Grosso et al., 2017), Tazarotene (Wang et al., 2018) and Verdinexor (Widman et al., 2018).

Moreover, HAdVs have wide applications as vector systems for drug delivery and vaccines development (such as human immunodeficiency virus (HIV), malaria, influenza, tuberculosis, and Ebola virus) due to their ability to infect a range of mammalian host cells, flexibility to accommodate substantial genetic recombinant inserts, and excellent safety profile (Curlin et al., 2010). Molecular associations have been established between HIV, HAdV and hepatitis B virus (HBV), and the correlation of molecular targets of antiretroviral drugs (ARVDs) with the HPs of HAdV might open avenues for prioritizing novel drug targets (Ghebremedhin, 2014; Kolawole et al., 2014; Singh et al., 2017).

A deeper understanding of metabolic impact of the existing antiretroviral drugs (ARVDs) in the host cell gene expression at the protein and RNA levels could be used to infer treatment, for posttranscriptional regulation in adenovirus-infected cells (Zhao et al., 2017). The aim of the present study was to identify critical host proteins and HAdV hypothetical proteins that could be developed as potential host-viral targets for anti-HAdV therapy. Here, the function of selected hypothetical proteins of HAdV based on phylogenetic relationship with the therapeutic targets of antiretroviral drugs of human immunodeficiency virus (HIV) was predicted computationally, and characterized the molecular dynamics and binding affinity of DNA polymerase of HAdV.

Methods

HAdV hypothetical protein sequence analysis

The list of 38 hypothetical proteins (HPs) of human adenovirus (HAdV) along with their Uniprot ID, corresponding genome and protein length were obtained from the literature (Naveed et al., 2017), and their sequences were obtained from UniProt database. the sequence properties such as isoelectric point, net charge, extinction coefficients and improbability of expression in inclusion bodies, were analyzed on EMBOSS Pepstats (www.ebi.ac.uk/Tools/seqstats/emboss_pepstats).

Antiretroviral drugs and their molecular targets preparation

The structures of the antiretroviral drugs (ARVDs) were obtained from NCBI PubChem Compound database (https://pubchem.ncbi.nlm.nih.gov/) in SMILES, and they were reconstructed and subjected to 3D structure optimization using ACDLabs/ChemSketch software, and the structures were saved in MOL format. File conversion from MOL format to PDB format was done using PyMoL v2.0.7. The molecular targets of these antiretroviral drugs were obtained from DrugBank database (https://go.drugbank.com/drugs/).

Phylogenetic analysis and structural modeling

The phylogenetic analysis of the HAdV hypothetical proteins with and without the molecular targets of the ARVDs, were done using ClustalO server (www.ebi.ac.uk/tools/clustalo). The structure of selected one HAdV and HBV proteins were modelled using SwissModel server: https://swissmodel.expasy.org/(Waterhouse et al., 2018; Studer et al., 2021). The modeled structure was analyzed using SAVES webserver v6.0 (http://saves.mbi.ucla.edu), tools such as ERRAT and PROCHECK, were used to check for the error function and stereochemical quality of the modelled protein structure. Based on the phylogenetic result, the sequence comparison of HAdV DNA polymerase and Human TERT and HLA-B target proteins of ARVD, were done on EMBOSS DotMatcher (www.ebi.ac.uk/Tools/seqstats/emboss_dotmatcher).

Human host metabolic genes network

It is important to appreciate the intracellular metabolism of the ARVDS. The gene IDs of the predicted target proteins for the ARVDs were compiled and used for expression network analyses (transcription factor enrichment analysis and protein-protein interaction network expansion and kinase enrichment analysis), using eXpression2Kinases (X2K) Web server https://maayanlab.cloud/X2K/(Clarke et al., 2018), where human was selected as background organism.

ARVDs human targets protein-protein interaction analysis

The gene ID of ARVDS molecular target proteins in humans and HPs of HAdV, were respectively analyzed for protein-protein interaction (PPI) profile on the STRING webserver (https://string-db.org/, Szklarczyk et al., 2021).

Molecular docking studies

The molecular docking studies were carried out according to the method of Fatoki et al. (2022), using the ARVDs that have multiorganism (HIV, HBV or Human) targets based on documented reports on drugbank database as shown in Table 1, where cidofovir (DB DB00369) and ganciclovir (DB01004) were used as control. Briefly, the target proteins and ligands were prepared for docking using AutoDock Tools (ADT) v1.5.6 (Morris et al., 2009) at default settings, with only polar hydrogen atoms, Kollman charges and Gasteiger charges added, and the output file was saved in pdbqt format. Molecular docking program AutoDock Vina v1.2.3 (Trott and Olson, 2010; Eberhardt et al., 2021) was employed to perform the active site docking experiment. After docking, close interactions of binding of the target with the ligands were analyzed and visualized on PyMol and PoseView webserver available at https://proteins.plus/(Stierand et al., 2006).

Protein molecular dynamics simulation

The flexibility dynamics of HAdV DNA Pol proteins from the PPI were analyzed using CABSflex2 server (http://biocomp.chem.uw.edu.pl/CABSflex2/) to evaluate the root mean square fluctuation (RMSF) of the structures. Small- and wide-angle X-ray scattering (SWAXS) curves based on explicit-solvent all-atom molecular dynamics (MD) simulations of the selected core proteins were analyzed on WAXSiS server: http://waxsis.uni-goettingen.de (Knight and Hub, 2015) to assess average value of radius of gyration (Rg).

Protein-ligand molecular dynamics simulation

Molecular dynamics simulations were performed for 100 nanoseconds using Desmond, a Package of Schrödinger LLC (Bowers et al., 2006; Schrödinger, 2018; Fatoki et al., 2023). The initial stage of protein and ligand complexes for molecular dynamics simulation were obtained from docking studies. The protein–ligand complexes were preprocessed using maestro’s protein preparation wizard, which also included optimization and minimization of complexes. All systems were prepared by the System Builder tool. Solvent Model with an orthorhombic box was selected as TIP3P (Transferable Intermolecular Interaction Potential 3 Points). The Optimized Potential for Liquid Simulations (OPLS)-2005 force field was used in the simulation (Shivakumar et al., 2010). The models were made neutral by adding counter ions 0.15 M NaCl to mimic the physiological conditions (Fatoki, 2022). The NPT ensemble (Isothermal-Isobaric: moles (N), pressure (P), and temperature (T) are conserved) with 300 K temperature and 1 atm pressure was select for complete simulation. The models were relaxed before the simulation. The trajectories were saved after every 100 ps during simulation, and post-simulation analysis of the trajectories were done to determine the root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), radius of gyration (Rg), solvent accessibility surface area (SASA), protein-ligand interaction profile. Also, prime molecular mechanics/generalized Born surface area (MMGBSA) was evaluated and binding free energy was calculated as follows:where is the total Prime energy, Hbond denote hydrogen bonding energy, Lipo is lipophilic energy, Packing represents pi-pi packing correction. SolvGB is generalized Born electrostatic solvation energy, and vdW is Van der Waals energy (Zhang et al., 2017; Schrödinger, 2019).

Results

Thirty-eight hypothetical proteins (HPs) of human adenovirus (HAdV) were used in this study. Among these 38 HPs, only HAdV-2 protein P03261 has the largest amino acid residues (aa.1198) with isoelectric point near neutrality (pH 7.4) as shown in Supplementary Table S1. The molecular targets of 30 ARVDs obtained are classified into three based on organism status: human (telomerase reverse transcriptase (hTERT) [UniProt ID: O14746], HLA class I histocompatibility antigen, B-57 alpha chain (HLA-B) [UniProt ID: P01889]), purine nucleoside phosphorylase (PNP) [UniProt ID: P00491], nuclear receptor subfamily 1 group I member 2 (NR1I2) (UniProt ID: O75469) and C-C chemokine receptor type (CCR5) [UniProt ID: P51681]); hepatitis B visus (protein P [UniProt ID: Q05486] and DNA polymerase (UniProt ID: P24024]); and human immunodeficiency virus (reverse transcriptase/RNaseH [UniProt ID: Q72547], integrase [UniProt ID: Q7ZJM1], Gag-Pol polyprotein [UniProt ID: P04585], HIV-1 protease (Pol polyprotein) [UniProt ID: Q72874; O90777], and envelope glycoprotein [UniProt ID: Q53I07, P12488]), as shown in Table 1.

The phylogenetic results in Figure 1, showed that HAdV proteins Q2KS67, Q3ZKV7, and A0A0B4SI61 are related to HBV DNA polymerase (P24024 and Q05486); P03261 is related to Human TERT (O14746) and HLA-B (P01889) genes; Q3ZKV3 is related to human CCR5 protein (P51681); A6MLW9 is related to human NR1I2 (O75469) and PNP (P00491) proteins. A0A0B4SH32 and P03263 are related to HIV-1 protease (O90777). I1V173, A0A0B4SIA5, Q3ZKV2 and Q2KSC0 are related to HIV-1 envelope glycoprotein (Q53I07) in functions. The relationship between HAdV (P03261) and two human proteins (TERT, O14746 and HLA-B, P01889) is shown in Figure 2, where greater similarity was found between HAdV (P03261) and hTERT (O14746).

FIGURE 1

FIGURE 2

The results of protein-protein interaction of human five molecular targets (PNP, TERT, CCR5, HLA-B, and NR1I2) of ARVDs are well-coordinated/networked with CD4, AHR, FKBP4, NR3C1, HSP90AA1, and STUB1 proteins in the anti-HIV infection mechanism (Figure 3). The gene network analysis of the metabolism genes involved in ARVD reveal the role of transcription factors (MYC, HNF4A, MAX, TP63, GATA1, REST, PPARD, FOXA2, and NFE2L2) and kinases (CDK1, CDK2, CDK4, HIPK, DNAPK, MAPK3 and MAPK14) in antiviral mechanism (Figure 4).

FIGURE 3

FIGURE 4

The model structure of HAdV DNA polymerase (HAdV DNa pol) was done on swiss-modelling server using Phi29 DNA polymerase (PDB ID: 1XHX) with resolution of 2.35Å, which has sequence identity and similarity of 19.58 and 0.29 respectively, in a range of 410-825 amino acid residues, forming a monomeric structure. The model on Swissmodel showed MolProbility score of 1.90 and Ramachandra favoured 83.60% of the amino acids. The overall quality factor of the model by ERRAT plot was 81.70% (Supplementary Figure S1) while the Ramachandran plot generated by PROCHECK showed that 80.5% of the amino acid residues are in the most favoured region, 16.5% in the allowed region, 1.8% in the general region, and 1.3% in the disallowed region. This model structure of HAdV DNa pol was also used for molecular dynamic simulation and docking analyses.

The antiviral drugs that met the criteria of multi-organismal targets were used for molecular docking study and their chemical structures are shown in Figure 5. Molecular docking studies provide a prediction of ligand binding status in static conditions. Molecular docking of HAdV DNA polymerase was done based on the results of phylogenetic analysis, and the results are shown in Table 2. The results showed that the free energy score of abacavir and zidovudine binding to HAdV DNA polymerase (−5.8 and −5.4 kcal mol^-1 respectively) are less than that of HIV reverse transcriptase/RNaseH (−6.5 and −6.2 kcal mol^-1 respectively), human telomerase reverse transcriptase (hTERT) (−6.2 kcal mol^-1) and human HLA class I histocompatibility antigen, B-57 alpha chain (−6.3 kcal mol^-1). The 3D structures of some of the binding interactions are shown in Figure 6, while the 2D structures are indicated in Figure 7. Also, the control drug, cidofovir and ganciclovir have less binding affinity for DNA polymerase of HAdV when compare to that of abacavir and zidovudine. Similarity was observed in the binding of abacavir and zidovudine to HAdV DNA polymerase (ASP742, ALA743, LEU772, ARG773 and VAL776), with that of zidovudine binding to human telomerase reverse transcriptase (TYR137, TYR168, TYR772, MET773, and GLN775). Also, similarity was observed in the binding of cidofovir and ganciclovir to HAdV DNA polymerase.

FIGURE 5

FIGURE 6

FIGURE 7

Molecular dynamic simulations (MDS) were carried out to predict the structural stability of the protein as well as the ligand binding status in the physiological environment. The results of MDS study of the HAdV DNa pol is showed in Figure 8. The results showed radius of gyration (Rg) of 31.98Å (3.198 nm) from SWAXS analysis, also native and the residue wise RMSF value of the free protein and protein bound to the ligands. The residues showing higher peaks correspond to loop regions, as identified from MD trajectories. RMSD values for the protein and ligand in complex are 6 Å and 1 Å respectively, and the Rg of the protein in complex with the ligand is about 6 Å (Figure 9). The results of interactions of various binding forces, SASA and PCA results were indicated in Figure 10 and Figure 11, where the protein SASA values were about 6,000 Å. Val738, Phe739, Leu771, Arg773, Gly774, Leu794 and Val795 were the most interacted amino acid residues of HAdV DNA pol—abacavir complex while Leu429, Asn432, Ile573 and Glu769 were the most interacted amino acid residues of HAdV DNA pol—zidovudine complex. Also, the results of MMGBSA binding energy of protein-ligand interaction showed that HAdV DNA pol—zidovudine complex is higher than that of DNA pol—abacavir complex as shown in Table 3.

FIGURE 8

FIGURE 9

FIGURE 10

FIGURE 11

Discussion

HAdVs have been implicated as infectious agents which are responsible for numerous diseases, including respiratory tract infections, ocular and gastrointestinal tract disorders (Kinchington et al., 2005). The HAdV early regions (E1A, E1B, E2A, E2B, E3 and E4) are the first viral proteins to be expressed during infection (Saha and Parks, 2020). The early E2 region of the HAdV genome consists of two transcription units, E2A encode DNA binding protein (DBP), whereas E2B encode terminal protein precursor (pTP), and DNA polymerase) that are required for viral replication (Robinson et al., 2008). The E2B gene of HAdV encodes DNA polymerase, which plays a key role in the viral genome replication and serves as a potent target for developing antiviral agents (Kumar et al., 2019).

Predicting the structure and function of hypothetical protein provides an opportunity to discover novel potential drug targets (Singh et al., 2022). Computational analysis has shown that two conserved domains; the polymerase/reverse transcriptase (RT) domain, and the C-terminal RNase H (RH) domain, are presence in all P proteins (Xiong and Eickbush, 1990). Previous bioinformatics study has characterized 38 randomly retrieved hypothetical proteins and assessed their physiochemical properties, subcellular localization, and predicted the function of six hypothetical proteins with UniProt IDs: P03269, P03261, P03263, Q83127, Q1L4D7 and I6LEV1 to be DNA terminal protein, DNA polymerase, DNA binding protein, adenovirus E3 region protein CR1 and adenoviral protein L1, respectively (Naveed et al., 2017).

The result of this study showed that antiretroviral drugs mechanism of action is centralized on the simultaneous up-regulating CD4 gene and inhibiting the HAdV DNA polymerase activity by downregulation of hTERT. Telomerase which is highly expressed in the vast majority of human cancers but not in most host tissues, is a unique DNA polymerase that catalyzes the addition of telomeric repeats (T2AG3)_n to the ends of chromosomes to maintain chromosomal integrity (Mo et al., 2003; Irving et al., 2004). Over 85% of hTERT gene, which codes for the catalytic component of the telomerase ribonucleoprotein complex, are involve in cancers (Meyerson et al., 1997; Irving et al., 2004; Doloff et al., 2008). A research report showed that expression of the E1A gene under the control of the hTERT promoter and subsequently viral replication restricted to tumor cells derived from hepatic and prostate tissues, induced significant tumor reduction and sometimes resulted in complete tumor regression (Irving et al., 2004; Doloff et al., 2008).

A study has shown that adenovirus E3-19k inhibits the phosphorylation of major histocompatibility complex (MHC) class I proteins in the rough endoplasmic reticulum (RER) with ability to bind to HLA molecules, during the infection of a human cell line (Lippe et al., 1991). E3-19K proteins with MHC (class I and class II-restricted) activity binds with high affinity to HLA-A and HLA-B molecules, and counteract activation of CD8⁺ and CD4⁺ T cells (Scheppler et al., 1989; Peek et al., 1994; Li et al., 2016). Thus, vaccines will be unsuitable due to declining of CD4⁺ T cell counts, but antiviral therapies are a promising therapeutic strategy against HAdV infection (Widman et al., 2018).

Several CD4⁺ and CD8⁺ T-cell epitopes from HAdV hexon have been identified and found to be highly conserved among different serotypes (Kumar et al., 2019). The clearance of human adenovirus viremia in patients after allogeneic stem cell transplantation, concurred with advent of a coordinated CD4⁺ and CD8⁺ T-cell response against adenovirus hexon epitopes (Zandvliet et al., 2010). Adenoviruses produce factors that block the synthesis and expression of HLA class I molecules on the cell surface, thereby inhibit the presentation of viral antigens attacked by CD8⁺ T-lymphocytes; and the infected cells acquire increased resistance to interferons and TNF-α (Cook and Radke, 2017). Reduction in the level of surface-expressed HLA class I cause Immunosuppression and it has been reported to occur during infection by several viruses including HAdV through mechanism that differs from inhibition of host cell protein synthesis. This is possible by downregulation of the surface expression of the T cell specific CD3 complex and HLA class I, which resulted in impairment of the presentation of viral peptides to the immune system and thus prevent the initiation of cell-mediated immunity to the virus (Boshkov et al., 1992; Peek et al., 1994).

Insight on the DEGs during adenovirus infection could provide rich clues to understanding the host’s response to adenovirus infection at the different stages of their replicative cycle. This study showed the metabolism genes involved in ARVD, which include kinases (CDK1, CDK2, CDK4, HIPK, DNAPK, MAPK3 and MAPK14) and transcription factors (MYC, HNF4A, MAX, TP63, GATA1, REST, PPARD, FOXA2, and NFE2L2), which must be well regulated in order to inhibit HAdV patho-mechanism. A study on HBV replication in human hepatic cells has published that FOXA1 and FOXA2 suppressed the expression of nuclear hormone receptors including PPARα, HNF4α, and RXRα (Okumura et al., 2015). A study has identified CD24, CD44, TP63, AHR, MAP3K14, RARB, PIMI, NEDD9, and few others, to be DEGs suppressed in HAdV-6 and HAdV-12, while CBX7 is an activated DEGs found in HAdV-12 (Wang et al., 2018). AHR is one of the DEG found involved in cell cycle control in HAdVs while CDK2 is one of the key regulators for the progression from G1 to the S phase, and c-Myc:Max is one of the consensus transcription factor expressed in HAdVs (Zhao et al., 2019). In fowl adenovirus serotype 4 (FAdV-4 strain GX-1), FABP3 and PPARG were reported to be upregulated DEGs at 21dpi of KEGG pathway (Ren et al., 2019).

In respect to HAdV antivirals molecular targets, the CDKs have been found promising because they regulate key aspects of the replicative cycles of a number of DNA and RNA viruses. Specifically, CDK2 was shown to promote the replication of HAdV, while CDK9 has been suggested to enhance transcription of the HAdV E1A gene (Dodge et al., 2021). MYC is a transcription factor that regulates energy metabolism through direct activation of metabolic genes, thus the ARVDS mechanism would be to deactivates MYC in the virus-infected patients. During hepatitis C virus (HCV) infection, upregulation of MYC was shown to suppresses the expression of N-Myc downstream-regulated gene 1 (NDRG1), and lead to decreased expression of NDRG1-specific kinase serum/glucocorticoid-regulated kinase 1 (SGK1) (Schweitzer et al., 2018).

Moreover, several DEGs were found downregulated during the progression of an adenovirus infection including DEGs involved in cell cycle and proliferation (MYC, SGK), in immune and stress response (NFE2L2, MAP2K3), in cellular structure and cell communication (CD44, KIF23), but about 50% of DEGs involved in metabolism (e.g., FKBP4, HSD17B2) were upregulated (Granberg et al., 2005). Insightfully, HAdV could activates MYC gene of the host to upregulates the expression of host cell glucose metabolic enzymes (specifically, hexokinase 2 and phosphofructokinase 1) through a process that mimics Warburg effect in cancer, to promote biosynthesis of nucleotide from glucose intermediates and ensure sustainable environments for its replication in primary lung epithelial cells (Thai et al., 2014). Thus, downregulation of MYC by ARVDS combination therapy will trigger expression of NDRG1, which inhibits HCV and HAdV replication at the stage of viral assembly. The results of this study showed that rilpivirine and ritonavir targets the human nuclear receptor subfamily 1 group I member 2 (NR1I2). Also, nuclear receptor subfamily 4, group A, member 1 (NR4A1) have been reported to be among the DEGs that showed decrease expression during the progression of an adenovirus infection (Granberg et al., 2005).

The result of this study showed that cidofovir and ganciclovir have less binding affinity for DNA polymerase of HAdV when compare to zidovudine and abacavir. Curated information on DrugBank database showed that cidofovir (DB00369) is an inhibitor of DNA polymerase catalytic subunit (UniProt ID: P08547) of human herpesvirus 5 (HHV-5) in cytomegalovirus (CMV) family; Ganciclovir (DB01004) is an inhibitor of DNA polymerase catalytic subunit (UniProt ID: P04293) of HHV-1. Ganciclovir is of little value as a HAdV therapeutic due to the fact that it requires activation by a viral thymidine kinase to be most effective, which is lacking in HAdV (Dodge et al., 2021). However, study has shown that HAdV co-infection with human CMV improved anti-HAdV activity of ganciclovir, due to CMV viral kinase-mediated activation of ganciclovir (Aguilar-Guisado et al., 2020).

Abacavir and Zidovudine inhibit the HIV-1 reverse transcriptase enzyme competitively and act as a chain terminator of DNA synthesis. Abacavir is a carbocyclic synthetic nucleoside analogue and an antiviral agent. Intracellularly, abacavir is phosphorylated by cellular enzymes to the active metabolite carbovir triphosphate, an analogue of deoxyguanosine-5′-triphosphate (dGTP). Carbovir triphosphate inhibits the activity of HIV-1 reverse transcriptase (RT) both by competing with the natural substrate dGTP and by its incorporation into viral DNA. Viral DNA growth is terminated because the incorporated nucleotide lacks a 3′-OH group, which is needed to form the 5′to 3′phosphodiester linkage essential for DNA chain elongation (De Clercq, 2002). Zidovudine, a structural analog of thymidine, is a prodrug that is phosphorylated to produce active metabolite, zidovudine triphosphate (ZDV-TP), which inhibits the activity of HIV-1 reverse transcriptase (RT) via DNA chain termination by competing with the natural substrate dGTP and incorporates itself into viral DNA. zidovudine is also a weak inhibitor of cellular DNA polymerase α and γ (De Clercq, 2010; Leeansyah et al., 2013). A study has shown that zidovudine could inhibit human telomerase reverse transcriptase (hTERT) (Leeansyah et al., 2013), while abacavir could target human HLA class I histocompatibility antigen, B-57 alpha chain protein (Yang et al., 2009).

Computational approaches such as in silico molecular docking and molecular dynamic simulation, provide opportunities for target discovery and drug repositioning are useful in modern pharmaceutical processes. Drug repurposing or repositioning is now a significantly relevant strategy to circumvent the roadblocks in novel drug discovery (Pushpakom et al., 2019; Dodge et al., 2021).

Molecular dynamics simulation (MDS) was performed to understand the structural and conformational changes in the protein and protein-ligand complex. Maestro-Desmond was used for the MSD in this study, and it is an interoperable tool built for biomolecular (protein, RNA and DNA) simulation (Schrodinger, 2018). RMSD plots is used to indicate whether there are no huge changes in the Cα backbone of the protein. The smaller amount of difference in the value of Rg as well as RMSD indicate firmly folded protein (Saini et al., 2021). Radius of gyration (Rg) indicates the overall compactness of the protein during the molecular dynamics. It is the distance between the center of mass of all atoms of protein and its terminal in a particular time interval. Low RMSF values of binding site residues indicate the stability of ligand binding to the protein. The MMGBSA was used to determine the binding energy from the free energies of the reactants and product of the reaction (Kollman et al., 2000).

Conclusion

The results of this study pinpointed abacavir and zidovudine as potential repurposing FDA drugs for inhibiting HAdV DNA polymerase. Combination therapy for treating multisystemic diseases such as neurodegeneration (i.e., amyotrophic lateral sclerosis) and parasitic (i.e., HIV-AIDS, malaria) illnesses could be well implemented for HAdV treatment. Combined use of existing approved fosters synergistic efficacy, ameliorate resistance to single drug, and improved tolerance based on complementarity in the pharmacokinetics and pharmacodynamics of the drug. Future work will be necessary to validate the combine therapy of abacavir, zidovudine and other existing investigational drug for treating HAdV infection such as brincidofovir, tazarotene, and verdinexor.

References

• 1

  Aguilar-GuisadoM.Marrugal-LorenzoJ. A.Berastegui-CabreraJ.MerinoL.PachónJ.SánchezCéspedesJ. (2020). In vitro co-infection by cytomegalovirus improves the antiviral activity of ganciclovir against human adenovirus. Int. J. Antimicrob. Agents56, 106046. 10.1016/j.ijantimicag.2020.106046

  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  BoshkovL. K.MacenJ. L.McFaddenG. (1992). Virus-induced loss of class I MHC antigens from the surface of cells infected with myxoma virus and malignant rabbit fibroma virus. J. Immunol.148, 881–887. 10.4049/jimmunol.148.3.881

  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  BowersK. J.ChowD. E.XuH.DrorR. O.EastwoodM. P.GregersenB. A.et al (2006). “Molecular dynamics—scalable algorithms for molecular dynamics simulations on commodity clusters,” in Proceedings of the 2006 ACM/IEEE Conference on Supercomputing—SC’06, Tampa, FL, USA, November 11 - 17, 2006, 11–17.

  • Google Scholar
  • View reference in article

• 4

  ClarkeD. J. B.KuleshovM. V.SchilderB. M.TorreD.DuffyM. E.KeenanA. B.et al (2018). eXpression2Kinases (X2K) Web: linking expression signatures to upstream cell signaling networks. Nucleic Acids Res.46 (1), 171–179. 10.1093/nar/gky458

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  CookJ.RadkeJ. (2017). Mechanisms of pathogenesis of emerging adenoviruses. F1000 Res.6, 90. 10.12688/f1000research.10152.1

  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  CurlinM. E.HuangM-L.LuX.CelumC. L.SanchezJ.SelkeS.et al (2010). Frequent detection of human adenovirus from the lower gastrointestinal tract in men who have sex with men. PLoS ONE5 (6), e11321. 10.1371/journal.pone.0011321

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  De ClercqE. (2002). New anti-HIV agents and targets. Med. Res. Rev.22 (6), 531–565. 10.1002/med.10021

  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  De ClercqE. (2010). Antiretroviral drugs. Curr. Opin. Pharmacol.10 (5), 507–515. 10.1016/j.coph.2010.04.011

  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  DodgeM. J.MacNeilK. M.TessierT. M.WeinbergJ. B.MymrykJ. S. (2021). Emerging antiviral therapeutics for human adenovirus infection: Recent developments and novel strategies. Antivir. Res.188, 105034. 10.1016/j.antiviral.2021.105034

  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  DoloffJ. C.WaxmanD. J.JounaidiY. (2008). Human telomerase reverse transcriptase promoter-driven oncolytic adenovirus with E1B-19 kDa and E1B-55 kDa gene deletions. Hum. Gene Ther.19, 1383–1399. 10.1089/hum.2008.056

  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  EberhardtJ.Santos-MartinsD.TillackA. F.ForliS. (2021). AutoDock Vina 1.2.0: New docking methods, expanded force field, and Python bindings. J. Chem. Inf. Model.61, 3891–3898. 10.1021/acs.jcim.1c00203

  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  FatokiT.ChukwuejimS.IbraheemO.OkeC.EjimaduB.OlaoyeI.et al (2022). Harmine and 7,8-dihydroxyflavone synergistically suitable for amyotrophic lateral sclerosis management: An insilico study. Res. Results Pharmacol.8 (3), 49–61. 10.3897/rrpharmacology.8.83332

  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  FatokiT. H.AjiboyeB. O.AremuA. O. (2023). Silico evaluation of the antioxidant, anti-inflammatory, and dermatocosmetic activities of phytoconstituents in licorice (Glycyrrhiza glabra L). Cosmetics10, 69. 10.3390/cosmetics10030069

  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  FatokiT. H. (2022). Effect of pH on structural dynamics of HMGCoA reductase and binding affinity to β-sitosterol. J. Biomol. Struct. Dyn., 1–7. 10.1080/07391102.2022.2067240

  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  FieldJ.GronostajskiR. M.HurwitzJ. (1984). Properties of the adenovirus DNA polymerase. J. Biol. Chem.259 (15), 9487–9495. 10.1016/s0021-9258(17)42726-8

  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  GhebremedhinB. (2014). Human adenovirus: Viral pathogen with increasing importance. Eur. J. Microbiol. Immunol.4 (1), 26–33. 10.1556/eujmi.4.2014.1.2

  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  GranbergF.SvenssonC.PetterssonU.ZhaoH. (2005). Modulation of host cell gene expression during onset of the late phase of an adenovirus infection is focused on growth inhibition and cell architecture. Virology343, 236–245. 10.1016/j.virol.2005.08.023

  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  GrossoF.StoilovP.LingwoodC.BrownM.CochraneA. (2017). Suppression of adenovirus replication by cardiotonic steroids. J. Virol.91, e01623–16. 10.1128/jvi.01623-16

  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  HierholzerJ. C. (1992). Adenoviruses in the immunocompromised host. Clin. Microbiol. Rev.5, 262–274. 10.1128/cmr.5.3.262

  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  IrvingJ.WangZ.PowellS.O’SullivanC.MokM.MurphyB.et al (2004). Conditionally replicative adenovirus driven by the human telomerase promoter provides broad-spectrum antitumor activity without liver toxicity. Cancer Gene Ther.11, 174–185. 10.1038/sj.cgt.7700666

  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  JamesL.VernonM. O.JonesR. C.StewartA.LuX.ZollarL. M.et al (2007). Outbreak of human adenovirus type 3 infection in a pediatric long-term care facility - Illinois. Clin. Infect. Dis.45, 416–420. 10.1086/519938

  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  KinchingtonP. R.RomanowskiE. G.JeroldG. Y. (2005). Prospects for adenovirus antivirals. J. Antimicrob. Chemother.55, 424–429. 10.1093/jac/dki057

  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  KnightC. J.HubJ. S. (2015). WAXSiS: a web server for the calculation of SAXS/WAXS curves based on explicit-solvent molecular dynamics. Nucleic Acids Research171 (W1), W225–W230. 10.1093/nar/gkv309

  • CrossRef
  • Google Scholar

• 24

  KolawoleO. M.OladosuT. O.AbdulkarimA. A.OkohA. I. (2014). Prevalence of adenovirus respiratory tract and hiv co-infections in patients attending the University of Ilorin, teaching hospital, Ilorin, Nigeria. BMC Res. Notes7, 870. 10.1186/1756-0500-7-870

  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  KollmanP. A.MassovaI.ReyesC.KuhnB.HuoS.ChongL.et al (2000). Calculating structures and free energies of complex molecules: Combining molecular mechanics and continuum models. Acc. Chem. Res.33, 889–897. 10.1021/ar000033j

  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  KumarR.TiwariA.PandeyG.TrivediR. N.KhanA. S.SaxenaM. K. (2019). “Adenovirus DNA polymerase: Structure, function, and prospects in diagnostics and therapeutics,” in Viral polymerases (Amsterdam: Elsevier), 429–455. 10.1016/B978-0-12-815422-9.00015-2

  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  La RosaA. M.ChamplinR. E.MirzaN.GajewskiJ.GiraltS.RolstonK. V.et al (2001). Adenovirus infections in adult recipients of blood and marrow transplants. Clin. Infect. Dis.32, 871–876. 10.1086/319352

  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  LeeansyahE.CameronP. U.SolomonA.TennakoonS.VelayudhamP.GouillouM.et al (2013). Inhibition of telomerase activity by human immunodeficiency virus (HIV) nucleos(t)ide reverse transcriptase inhibitors: A potential factor contributing to HIV-associated accelerated aging. J. Infect. Dis.207 (7), 1157–1165. 10.1093/infdis/jit006

  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  LiL.SantarsieroB. D.BouvierM. (2016). Structure of the adenovirus type 4 (species E) E3-19K/HLA-A2 complex reveals species-specific features in MHC class I recognition. J. Immunol.197, 1399–1407. 10.4049/jimmunol.1600541

  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  LippeR.LukeE.KuahY. T.LomasC.JefferiesW. A. (1991). Adenovirus infection inhibits the phosphorylation of major histocompatibility complex class I proteins. J. Exp. Med.174, 1159–1166. 10.1084/jem.174.5.1159

  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  MeyersonM.CounterC. M.EatonE. N.EllisenL. W.SteinerP.CaddleS. D.et al (1997). hEST2, the putative human telomerase catalytic subunit gene is upregulated in tumor cells and during immortalization. Cell90, 785–795. 10.1016/s0092-8674(00)80538-3

  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  MoY.GanY.SongS.JohnstonJ.XiaoX.WientjesM. G.et al (2003). Simultaneous targeting of telomeres and telomerase as a cancer therapeutic approach. Cancer Res.63, 579–585.

  • Google Scholar
  • View reference in article

• 33

  MorrisG. M.HueyR.LindstromW.SannerM. F.BelewR. K.GoodsellD. S.et al (2009). AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J. Comput. Chem.30 (16), 2785–2791. 10.1002/jcc.21256

  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  NagataK.GuggenheimerR. A.HurwitzJ. (1983). Adenovirus DNA replication in vitro: Synthesis of full-length DNA with purified proteins. Proc. Natl. Acad. Sci. U. S. A.80, 4266–4270. 10.1073/pnas.80.14.4266

  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  NaveedM.TehreemS.UsmanM.ChaudhryZ.AbbasG. (2017). Structural and functional annotation of hypothetical proteins of human adenovirus: Prioritizing the novel drug targets. BMC Res. Notes10, 706. 10.1186/s13104-017-2992-z

  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  OkumuraN.IkedaM.SatohS.DansakoH.SugiyamaM.MizokamiM.et al (2015). Negative regulation of Hepatitis B virus replication by forkhead box protein A in human hepatoma cells. FEBS Lett.589, 1112–1118. 10.1016/j.febslet.2015.03.022

  • CrossRef
  • Google Scholar
  • View reference in article

• 37

  PeekR.WestphalJ. R.PruijnG. J. M.Van Der KempA. J. W.Van VenrooijW. J. (1994). Adenovirus infection induces loss of HLA class I and CD3 antigens, but does not induce cell surface presentation of the La (SS-B) autoantigen. Clin. Exp. Immunol.96, 395–402. 10.1111/j.1365-2249.1994.tb06041.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  PushpakomS.IorioF.EyersP. A.EscottK. J.HopperS.WellsA.et al (2019). Drug repurposing: Progress, challenges and recommendations. Nat. Rev. Drug Discov.18, 41–58. 10.1038/nrd.2018.168

  • CrossRef
  • Google Scholar
  • View reference in article

• 39

  RenG.WangH.HuangM.YanY.LiuF.ChenR., (2019). Transcriptome analysis of fowl adenovirus serotype 4 infection in chickens. Virus Genes55, 619–629. 10.1007/s11262-019-01676-w

  • CrossRef
  • Google Scholar

• 40

  RobinsonC. M.SinghG.LeeJ. Y.DehghanS.RajaiyaJ.LiuE. B.et al (2013). Molecular evolution of human adenoviruses. Sci. Rep.3, 1812. 10.1038/srep01812

  • CrossRef
  • Google Scholar

• 41

  RobinsonC. M.Fatemeh ShariatiF.GillaspyA. F.DyerD. W.ChodoshJ. (2008). Genomic and bioinformatics analysis of human adenovirus type 37: New insights into corneal tropism. BMC Genomics9, 213. 10.1186/1471-2164-9-213

  • CrossRef
  • Google Scholar

• 42

  SahaB.ParksR. J. (2020). Recent advances in novel antiviral therapies against human adenovirus. Microorganisms8, 1284. 10.3390/microorganisms8091284

  • CrossRef
  • Google Scholar

• 43

  SainiG.DalalV.GuptaD. N.SharmaN.KumarP.SharmaA. K. (2021). A molecular docking and dynamic approach to screen inhibitors against ZnuA1 of Candidatus Liberibacter asiaticus. Mol. Simul.47 (6), 510–525. 10.1080/08927022.2021.1888948

  • CrossRef
  • Google Scholar

• 44

  SchepplerJ. A.NicholsonJ. K. A.SwanD. C.Ahmed-AnsariA.McDougalJ. S. (1989). Down-modulation of MHC-I in a CD4+ T cell line, CEM-E5, after HIV-1 infection. J. Immunol.143, 2858–2866. 10.4049/jimmunol.143.9.2858

  • CrossRef
  • Google Scholar

• 45

  Schrödinger (2018). Schrödinger release 2018-3. Desmond molecular dynamics system. New York, NY: D.E. Shaw research. Maestro Desmond Interoperability Tools, Schrödinger, New York, NY, 2018.

  • Google Scholar

• 46

  Schrödinger (2019). What do all the Prime MM-GBSA energy properties mean?Available at: www.schrodinger.com/kb/1875.

  • Google Scholar

• 47

  SchweitzerC. J.ZhangF.BoyerA.ValdezK.CamM.LiangT. J. (2018). N-Myc downstream-regulated gene 1 restricts hepatitis C virus propagation by regulating lipid droplet biogenesis and viral assembly. J. Virol.92, e01166–17. 10.1128/JVI.01166-17

  • CrossRef
  • Google Scholar

• 48

  ScottM. K.ChommanardC.LuX.AppelgateD.GrenzL.SchneiderE.et al (2016). Human adenovirus associated with severe respiratory infection, Oregon, USA, 2013–2014. Emerg. Infect. Dis.22 (6), 1044–1051. 10.3201/eid2206.151898

  • CrossRef
  • Google Scholar

• 49

  ShivakumarD.WilliamsJ.WuY.DammW.ShelleyJ.ShermanW. (2010). Prediction of absolute solvation free energies using molecular dynamics free energy perturbation and the OPLS force field. J. Chem. Theory Comput.6, 1509–1519. 10.1021/ct900587b

  • CrossRef
  • Google Scholar

• 50

  SinghK. P.CraneM.AudsleyJ.LewinS. R.SasadeuszJ. (2017). HIV-hepatitis B virus co-infection: Epidemiology, pathogenesis and treatment. AIDS31 (15), 2035–2052. 10.1097/QAD.0000000000001574

  • CrossRef
  • Google Scholar

• 51

  SinghV.DhankharP.DalalV.TomarS.KumarP. (2022). Insilico functional and structural annotation of hypothetical protein from Klebsiella pneumonia: A potential drug target. J. Mol. Graph. Modell.116, 108262. 10.1016/j.jmgm.2022.108262

  • CrossRef
  • Google Scholar

• 52

  StierandK.MaaßP. C.RareyM. (2006). Molecular complexes at a glance: Automated generation of two-dimensional complex diagrams. Bioinformatics22, 1710–1716. 10.1093/bioinformatics/btl150

  • CrossRef
  • Google Scholar

• 53

  StuderG.TaurielloG.BienertS.BiasiniM.JohnerN.SchwedeT. (2021). ProMod3 - a versatile homology modelling toolbox. PLOS Comp. Biol.17 (1), e1008667. 10.1371/journal.pcbi.1008667

  • CrossRef
  • Google Scholar

• 54

  SzklarczykD.GableA. L.NastouK. C.LyonD.KirschR.PyysaloS.et al (2021). The STRING database in 2021: Customizable protein–protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res.49, D605–D612. 10.1093/nar/gkaa1074

  • CrossRef
  • Google Scholar

• 55

  ThaiM.GrahamN. A.BraasD.NehilM.KomisopoulouE.KurdistaniS. K.et al (2014). Adenovirus E4ORF1-induced MYC activation promotes host cell anabolic glucose metabolism and virus replication. Cell Metab.19 (4), 694–701. 10.1016/j.cmet.2014.03.009

  • CrossRef
  • Google Scholar

• 56

  TrottO.OlsonA. J. (2010). AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem.31 (2), 455–461. 10.1002/jcc.21334

  • CrossRef
  • Google Scholar

• 57

  WangX.ZhangQ.ZhouZ.LiuM.ChenY.LiJ.et al (2018). Retinoic acid receptor β, a potential therapeutic target in the inhibition of adenovirus replication. Antivir. Res.152, 84–93. 10.1016/j.antiviral.2018.01.014

  • CrossRef
  • Google Scholar

• 58

  WaterhouseA.BertoniM.BienertS.StuderG.TaurielloG.GumiennyR.et al (2018). SWISS-MODEL: Homology modelling of protein structures and complexes. Nucleic Acids Res.46 (1), W296–W303. 10.1093/nar/gky427

  • CrossRef
  • Google Scholar

• 59

  WayeM. M. Y.SingC. W. (2010). Anti-viral drugs for human adenoviruses. Pharmaceuticals3, 3343–3354. 10.3390/ph3103343

  • CrossRef
  • Google Scholar

• 60

  WidmanD. G.GornisiewiczS.ShachamS.TamirS. (2018). In vitro toxicity and efficacy of verdinexor, an exportin 1 inhibitor, on opportunistic viruses affecting immunocompromised individuals. PLoS ONE13 (10), e0200043. 10.1371/journal.pone.0200043

  • CrossRef
  • Google Scholar

• 61

  WongM-L.HsuM-T. (1990). Involvement of topoisomerases in replication, transcription, and packaging of the linear adenovirus genome. J. Virology64 (2), 691–699. 10.1128/jvi.64.2.691-699.1990

  • CrossRef
  • Google Scholar

• 62

  XiongY.EickbushT. H. (1990). Origin and evolution of retroelements based upon their reverse transcriptase sequences. EMBO J.9, 3353–3362. 10.1002/j.1460-2075.1990.tb07536.x

  • CrossRef
  • Google Scholar

• 63

  YangL.ChenJ.HeL. (2009). Harvesting candidate genes responsible for serious adverse drug reactions from a chemical-protein interactome. PLoS Comput. Biol.5 (7), e1000441. 10.1371/journal.pcbi.1000441

  • CrossRef
  • Google Scholar

• 64

  ZandvlietM. L.FalkenburgJ. H.van LiemptE.Veltrop-DuitsL. A.LankesterA. C.KalpoeJ. S.et al (2010). Combined CD8+ and CD4+ adenovirus hexon-specific T cells associated with viral clearance after stem cell transplantation as treatment for adenovirus infection. Haematologica95 (11), 1943–1951. 10.3324/haematol.2010.022947

  • CrossRef
  • Google Scholar

• 65

  ZhangX.Perez-SanchezH.LightstoneF. C. (2017). A comprehensive docking and MM/GBSA rescoring study of ligand recognition upon binding antithrombin. Curr. Top. Med. Chem.17, 1631–1639. 10.2174/1568026616666161117112604

  • CrossRef
  • Google Scholar

• 66

  ZhaoH.ChenM.ValdésA.LindS. B.PetterssonU. (2019). Transcriptomic and proteomic analyses reveal new insights into the regulation of immune pathways during adenovirus type 2 infection. BMC Microbiol.19, 15. 10.1186/s12866-018-1375-5

  • CrossRef
  • Google Scholar

• 67

  ZhaoH.KonzerA.MiJ.ChenM.PetterssonU.LindS. B. (2017). Posttranscriptional regulation in adenovirus infected cells. J. Proteome Res.16, 872–888. 10.1021/acs.jproteome.6b00834

  • CrossRef
  • Google Scholar