Transcription Factor Binding Activity Research Articles

Weakened tanning ability is an important mechanism for evolutionary skin lightening in East Asians

The evolution of light-skin pigmentation among Eurasians is considered as an adaptation to the high-latitude environments. East Asians are ideal populations for studying skin color evolution because of the complex environment of East Asia. Here, we report a strong selection signal for the pigmentation gene phenylalanine hydroxylase (PAH) in light-skinned Han Chinese individuals. The intron mutation rs10778203 in PAH is enriched in East Asians and is significantly associated with skin color of the back of the hand in Han Chinese males (P < 0.05). In vitro luciferase and transcription factor binding assays show that the ancestral allele of rs10778203 could bind to SMAD2 and has a significant enhancer activity for PAH. However, the derived T allele (the major allele in East Asians) of rs10778203 decreases the binding activity of transcription factors and enhancer activity. Meanwhile, the derived T allele of rs10778203 shows a weaker ultraviolet radiation response in A375 cells and zebrafish embryos. Furthermore, rs10778203 decreases melanin production in transgenic zebrafish embryos after ultraviolet B (UVB) treatment. Collectively, PAH is a potential pigmentation gene that regulates skin tanning ability. Natural selection has enriched the adaptive allele, resulting in weakened tanning ability in East Asians, suggesting a unique genetic mechanism for evolutionary skin lightening in East Asians.

Bisphenol A regulates bladder cells responses via control of G2/M-phase cell cycle, apoptotic signaling, MAPK pathway, and transcription factor-associated MMP modulation.

Bisphenol A (BPA), an exogenous endocrine-disrupting chemical,is widely used to produce polycarbonate plastics. The widely used BPA has been detected in human urine samples, raising public anxiety about the detrimental effects of BPA on the bladder. In this study, we explored regulatory mechanisms for the adverse effects of BPA in human bladder BdFCand T24 cells. BPA induced extrinsic and intrinsic apoptosis and G2/M cell cycle arrest caused by the ATM-CHK1/CHK2-CDC25c-CDC2 signaling, which ultimately inhibited the growth of human bladder cells. We also found that BPA decreased the binding activity of AP-1 and NF-κB transcription factors in human bladder cells, which inhibited migration and invasion through matrix metallopeptidase-2 and -9 inactivation. Phosphorylation of MAPKs was implicated with BPA-mediated detrimental effects in human bladder cells. Collectively, our results provide a novel explanation for the underlying molecular mechanisms that BPA induces cytotoxicity in human bladder cells.

Bioinformatics identification of potential biomarkers and therapeutic targets for ischemic stroke and vascular dementia

Ischemic stroke and vascular dementia, as common cerebrovascular diseases, with the former causing irreversible neurological damage and the latter causing cognitive and memory impairment, are closely related and have long received widespread attention. Currently, the potential causative genes of these two diseases have yet to be investigated, and effective early diagnostic tools for the diseases have not yet emerged. In this study, we screened new potential biomarkers and analyzed new therapeutic targets for both diseases from the perspective of immune infiltration. Two gene expression profiles on ischemic stroke and vascular dementia were obtained from the NCBI GEO database, and key genes were identified by LASSO regression and SVM-RFE algorithms, and key genes were analyzed by GO and KEGG enrichment. The CIBERSORT algorithm was applied to the gene expression profile species of the two diseases to quantify the 24 subpopulations of immune cells. Moreover, logistic regression modeling analysis was applied to illustrate the stability of the key genes in the diagnosis. Finally, the key genes were validated using RT-PCR assay. A total of 105 intersecting DEGs genes were obtained in the 2 sets of GEO datasets, and bioinformatics functional analysis of the intersecting DEGs genes showed that GO was mainly involved in the purine ribonucleoside triphosphate metabolic process，respiratory chain complex，DNA−binding transcription factor binding and active transmembrane transporter activity. KEGG is mainly involved in the Oxidative phosphorylation, cAMP signaling pathway. The LASSO regression algorithm and SVM-RFE algorithm finally obtained three genes, GAS2L1, ARHGEF40 and PFKFB3, and the logistic regression prediction model determined that the three genes, GAS2L1 (AUC: 0.882), ARHGEF40 (AUC: 0.867) and PFKFB3 (AUC: 0.869), had good diagnostic performance. Meanwhile, the two disease core genes and immune infiltration were closely related, GAS2L1 and PFKFB3 had the highest positive correlation with macrophage M1 (p < 0.001) and the highest negative correlation with mast cell activation (p = 0.0017); ARHGEF40 had the highest positive correlation with macrophage M1 and B cells naive (p < 0.001), the highest negative correlation with B cell memory highest correlation (p = 0.0047). RT-PCR results showed that the relative mRNA expression levels of GAS2L1, ARHGEF40, and PFKFB3 were significantly elevated in the populations of both disease groups (p < 0.05). Immune infiltration-based models can be used to predict the diagnosis of patients with ischemic stroke and vascular dementia and provide a new perspective on the early diagnosis and treatment of both diseases.

Genomic transcription factor binding site selection is edited by the chromatin remodeling factor CHD4.

Biologically precise enhancer licensing by lineage-determining transcription factors enables activation of transcripts appropriate to biological demand and prevents deleterious gene activation. This essential process is challenged by the millions of matches to most transcription factor binding motifs present in many eukaryotic genomes, leading to questions about how transcription factors achieve the exquisite specificity required. The importance of chromatin remodeling factors to enhancer activation is highlighted by their frequent mutation in developmental disorders and in cancer. Here, we determine the roles of CHD4 inenhancer licensing and maintenance in breast cancer cells and during cellular reprogramming. In unchallenged basal breast cancer cells, CHD4 modulates chromatin accessibility. Its depletion leads to redistribution of transcription factors to previously unoccupied sites. During cellular reprogramming induced by the pioneer factor GATA3, CHD4 activity is necessary to prevent inappropriate chromatin opening. Mechanistically, CHD4 promotes nucleosome positioning over GATA3 binding motifs to compete with transcription factor-DNA interaction. We propose that CHD4 acts as a chromatin proof-reading enzyme that prevents unnecessary gene expression by editing chromatin binding activities of transcription factors.

Genome-wide association study of susceptibility to hospitalised respiratory infections

Background: Globally, respiratory infections contribute to significant morbidity and mortality. However, genetic determinants of respiratory infections are understudied and remain poorly understood. Methods: We conducted a genome-wide association study in 19,459 hospitalised respiratory infection cases and 101,438 controls from UK Biobank (Stage 1). We followed-up well-imputed top signals from our Stage 1 analysis in 50,912 respiratory infection cases and 150,442 controls from 11 cohorts (Stage 2). We aggregated effect estimates across studies using inverse variance-weighted meta-analyses. Additionally, we investigated the function of the top signals in order to gain understanding of the underlying biological mechanisms. Results: From our Stage 1 analysis, we report 56 signals at P&lt;5 ×10 -6, one of which was genome-wide significant ( P&lt;5 ×10 -8). The genome-wide significant signal was in an intron of PBX3, a gene that encodes pre-B-cell leukaemia transcription factor 3, a homeodomain-containing transcription factor. Further, the genome-wide significant signal was found to colocalise with gene-specific expression quantitative trait loci (eQTLs) affecting expression of PBX3 in lung tissue, where the respiratory infection risk alleles were associated with decreased PBX3 expression in lung tissue, highlighting a possible biological mechanism. Of the 56 signals, 40 were well-imputed in UK Biobank and were investigated in Stage 2. None of the 40 signals replicated, with effect estimates attenuated. Conclusions: Our Stage 1 analysis implicated PBX3 as a candidate causal gene and suggests a possible role of transcription factor binding activity in respiratory infection susceptibility. However, the PBX3 signal, and the other well-imputed signals, did not replicate in the meta-analysis of Stages 1 and 2. Significant phenotypic heterogeneity and differences in study ascertainment may have contributed to this lack of statistical replication. Overall, our study highlighted putative associations and possible biological mechanisms that may provide insight into respiratory infection susceptibility.

Benzimidazole based hybrids as privileged candidates for topoII inhibition: Design, synthesis and molecular modeling studies

Two series of new benzimidazole derivatives were designed and synthesized. The antiproliferative activity of the new derivatives was screened against five human cancer cell lines namely; HepG-2, MCF-7, PC-3, HCT-116 and Hela. Relative to Doxorubicin (IC50 = 6.72 µM), compounds 18 and 19 showed broad spectrum anticancer effect against all the tested cancer cell lines, with IC50 range of 5.16–20.35 µM, however they were safe to WI-38 normal cells with IC50 of 32.61 µM and 50.62 µM, respectively. Inspection of human TopoII inhibitory activity disclosed that compound 18 was the most potent (IC50 = 6.9 µM) compared to Staurosporine (IC50 = 4.64 µM). Further mechanistic studies implied that compound 18 triggered the apoptotic death of MCF-7 cells via decreasing Bcl-2, activating caspases 8 and 9, and increasing the accumulation of ROS. In addition, it caused cell cycle arrest at the G1-phase and lowered the binding activity of transcription factors’ protein to DNA. The biological results were assisted by molecular modeling studies. The physicochemical parameters of compounds 18 and 19 conferred convenient drug-like properties.

Characterizing Different Modes of Interplay Between Rap1 and H3 Using Inducible H3-depletion Yeast

Histones and transcription factors (TFs) are two important DNA-binding proteins that interact, compete, and together regulate transcriptional processes in response to diverse internal and external stimuli. Condition-specific depletion of histones in Saccharomyces cerevisiae using a galactose-inducible H3 promoter provides a suitable framework for examining transcriptional alteration resulting from reduced nucleosome content. However, the effect on DNA binding activities of TFs is yet to be fully explored. In this work, we combine ChIP-seq of H3 with RNA-seq to elucidate the genome-scale relationships between H3 occupancy patterns and transcriptional dynamics before and after global H3 depletion. ChIP-seq of Rap1 is also conducted in the H3-depletion and control treatments, to investigate the interplay between this master regulator TF and nucleosomal H3, and to explore the impact on diverse transcriptional responses of different groups of target genes and functions. Ultimately, we propose a working model and testable hypotheses regarding the impact of global and local H3 depletion on transcriptional changes. We also demonstrate different potential modes of interaction between Rap1 and H3, which sheds light on the potential multifunctional regulatory capabilities of Rap1 and potentially other pioneer factors.

Genome-wide census of ATF4 binding sites and functional profiling of trait-associated genetic variants overlapping ATF4 binding motifs.

Activating Transcription Factor 4 (ATF4) is an important regulator of gene expression in stress responses and developmental processes in many cell types. Here, we catalogued ATF4 binding sites in the human genome and identified overlaps with trait-associated genetic variants. We probed these genetic variants for allelic regulatory activity using a massively parallel reporter assay (MPRA) in HepG2 hepatoma cells exposed to tunicamycin to induce endoplasmic reticulum stress and ATF4 upregulation. The results revealed that in the majority of cases, the MPRA allelic activity of these SNPs was in agreement with the nucleotide preference seen in the ATF4 binding motif from ChIP-Seq. Luciferase and electrophoretic mobility shift assays in additional cellular models further confirmed ATF4-dependent regulatory effects for the SNPs rs532446 (GADD45A intronic; linked to hematological parameters), rs7011846 (LPL upstream; myocardial infarction), rs2718215 (diastolic blood pressure), rs281758 (psychiatric disorders) and rs6491544 (educational attainment). CRISPR-Cas9 disruption and/or deletion of the regulatory elements harboring rs532446 and rs7011846 led to the downregulation of GADD45A and LPL, respectively. Thus, these SNPs could represent examples of GWAS genetic variants that affect gene expression by altering ATF4-mediated transcriptional activation.

Transcription factors direct epigenetic reprogramming at specific loci in human cancers.

The characterization of epigenetic changes during cancer development and progression led to notable insights regarding the roles of cancer-specific epigenetic reprogramming. Recent studies showed that transcription factors (TFs) are capable to regulate epigenetic reprogramming at specific loci in different cancer types through their DNA-binding activities. However, the causal association of dynamic histone modification change mediated by TFs is still not well elucidated. Here we evaluated the impacts of 636 transcription factor binding activities on histone modification in 24 cancer types. We performed Instrumental Variables analysis by using genetic lesions of TFs as our instrumental proxies, which previously discovered to be associated with histone mark activities. As a result, we showed a total of 6 EpiTFs as strong directors of epigenetic reprogramming of histone modification in cancers, which alters the molecular and clinical phenotypes of cancer. Together our findings highlight a causal mechanism driven by the TFs and genome-wide histone modification, which is relevant to multiple status of oncogenesis.

A sticky-end probe biosensor for homogeneous detection of transcription factor binding activity.

Transcription factors are essential regulators of various physiological and pathological processes. However, detecting transcription factor-DNA binding activities is often time-consuming and labor-intensive. Homogeneous biosensors that are compatible with mix-and-measure protocols have the potential to simplify the workflow for therapeutic screening and disease diagnostics. In this study, we apply a combined computational-experimental approach to investigate the design of a sticky-end probe biosensor, where the transcription factor-DNA complex stabilizes the fluorescence resonance energy transfer signal of the donor-acceptor pair. We design a sticky-end biosensor for the SOX9 transcription factor based on the consensus sequence and characterize its sensing performance. A systems biology model is also developed to investigate the reaction kinetics and optimize the operating conditions. Taken together, our study provides a conceptual framework for the design and optimization of sticky-end probe biosensors for homogeneous detection of transcription factor-DNA binding activity.

Tomatidine targets ATF4-dependent signaling and induces ferroptosis to limit pancreatic cancer progression

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive cancer with high metastasis and therapeutic resistance. Activating transcription factor 4 (ATF4), a master regulator of cellular stress, is exploited by cancer cells to survive. Prior research and data reported provide evidence that high ATF4 expression correlates with worse overall survival in PDAC. Tomatidine, a natural steroidal alkaloid, is associated with inhibition of ATF4 signaling in multiple diseases. Here, we discovered that invitro and invivo tomatidine treatment of PDAC cells inhibits tumor growth. Tomatidine inhibited nuclear translocation of ATF4 and reduced the transcriptional binding of ATF4 with downstream promoters. Tomatidine enhanced gemcitabine chemosensitivity in 3D ECM-hydrogels and invivo. Tomatidine treatment was associated with induction of ferroptosis signaling validated by increased lipid peroxidation, mitochondrial biogenesis, and decreased GPX4 expression in PDAC cells. This study highlights a possible therapeutic approach utilizing a plant-derived metabolite, tomatidine, to target ATF4 activity in PDAC.

Identification of genes involved in regulating the development of feathered feet in chicken embryo

The genetic and developmental factors driving the diverse distribution and morphogenesis of feathers and scales on bird feet are yet unclear. Within a single species, Guangxi domestic chickens exhibit dramatic variety in feathered feet, making them an accessible model for research into the molecular basis of variations in skin appendages. In this study, we used H&E staining to observe the morphogenesis of feathered feet, scaled feet and wings skin at different embryonic stages in Longsheng-Feng chickens and Guangxi Partridge chickens. We selected 4 periods (E6, E7, E8, and E12) that play an important role in feather development and performed transcriptome sequencing to screen for candidate genes associated with feathered feet. Through comparison and analysis of transcriptome data, we identified a set of differently expressed genes (DGEs), which were enriched in appendage organ development, hindlimb morphogenesis, activation of transcription factor binding, and binding of sequence-specific DNA in the cis-regulatory region. In addition, we identified some feathered feet-related genes by analyzing the classical signaling pathways that regulate feather development. Finally, we identified candidate genes that regulate feathered feet formation, which include TBX5, PITX1, ZIC1, FGF20, WNT11, WNT7A, WNT16, and SHH. Interestingly, we found that TBX5 was significantly overexpressed in the skin of the feathered feet and had the highest expression at E7 (P < 0.01), whereas PITX1 expression was significantly reduced at E7(P < 0.01). It is hypothesized that TBX5 and PITX1 regulate the development of hair follicles through the Wnt/β-catenin signaling pathway at E7. Our results provide a theoretical basis for investigating the molecular regulatory mechanisms underlying the formation of chicken feathered feet.

Using Synthetic DNA Libraries to Investigate Chromatin and Gene Regulation.

Despite the recent explosion in genome-wide studies in chromatin and gene regulation, we are still far from extracting a set of genetic rules that can predict the function of theregulatory genome. One major reason for this deficiency is that gene regulation is a multi-layered process that involves an enormous variable space, which cannot be fully explored using native genomes. This problem can be partially solved by introducing synthetic DNA libraries into cells, a method that can test the regulatory roles of thousands to millions of sequences with limited variables. Here, we review recent applications of this method to study transcription factor (TF) binding, nucleosome positioning, and transcriptional activity. We discuss the design principles, experimental procedures, and major findings from these studies and compare the pros and cons of different approaches.

Functional analysis of selenok, selenot and selenop promoters and their regulation by selenium in yellow catfish Pelteobagrus fulvidraco

The selenok, selenot and selenop are three key selenoproteins involved in stress response. Our study, using the yellow catfish Pelteobagrus fulvidraco as the experimental animal, obtained the 1993-bp, 2000-bp and 1959-bp sequences of selenok, selenot and selenop promoters, respectively, and predicted the binding sites of several transcriptional factors on their promoters, such as Forkhead box O 4 (FoxO4), activating transcription factor 4 (ATF4), Kruppel-like factor 4 (KLF4) and nuclear factor erythroid 2-related factor 2 (NRF2). Selenium (Se) increased the activities of the selenok, selenot and selenop promoters. FoxO4 and Nrf2 can directly bind with selenok promoter and controlled selenok promoter activities positively; KLF4 and Nrf2 can directly bind with selenot promoter and controlled selenot promoter activities positively; FoxO4 and ATF4 can directly bind to selenop promoter and regulated selenop promoter activities positively. Se promoted FoxO4 and Nrf2 binding to selenok promoter, KLF4 and Nrf2 binding to selenot promoter, and FoxO4 and ATF4 binding to selenop promoter. Thus, we provide the first evidence for FoxO4 and Nrf2 bindnig elements in selenok promoter, KLF4 and Nrf2 binding elements in selenot promoter, and FoxO4 and ATF4 binding elements in selenop promoter, and offer novel insight into regulatory mechanism of these selenoproteins induced by Se.

A potential histone-chaperoneactivity for the MIER1 histone deacetylase complex.

Histone deacetylases 1 and 2 (HDAC1/2) serve as the catalytic subunit of six distinct families of nuclear complexes. These complexes repress gene transcription through removing acetyl groups from lysine residues in histone tails. In addition to the deacetylase subunit, these complexes typically contain transcription factor and/or chromatin binding activities. The MIER:HDAC complex has hitherto been poorly characterized. Here, we show that MIER1 unexpectedly co-purifies with an H2A:H2B histone dimer. We show that MIER1 is also able to bind a complete histone octamer. Intriguingly, we found that a larger MIER1:HDAC1:BAHD1:C1QBP complex additionally co-purifies with an intact nucleosome on which H3K27 is either di- or tri-methylated. Together this suggests that the MIER1 complex acts downstream of PRC2 to expand regions of repressed chromatin and could potentially deposit histone octamer onto nucleosome-depleted regions of DNA.

Correction of transposase sequence bias in ATAC-seq data with rule ensemble modeling.

Chromatin accessibility assays have revolutionized the field of transcription regulation by providing single-nucleotide resolution measurements of regulatory features such as promoters and transcription factor binding sites. ATAC-seq directly measures how well the Tn5 transposase accesses chromatinized DNA. Tn5 has a complex sequence bias that is not effectively scaled with traditional bias-correction methods. We model this complex bias using a rule ensemble machine learning approach that integrates information from many input k-mers proximal to the ATAC sequence reads. We effectively characterize and correct single-nucleotide sequence biases and regional sequence biases of the Tn5 enzyme. Correction of enzymatic sequence bias is an important step in interpreting chromatin accessibility assays that aim to infer transcription factor binding and regulatory activity of elements in the genome.

The role of pyrethroid derivatives in autophagy and apoptosis crosstalk signaling and potential risk for malignancies.

Pyrethroids and its derivatives widespread and uncontrolled continuous use has influenced multiple deleterious effects resulting in as a potential risk factor causing damage to the organ systems. Allethrin and prallethrin are extensively used yet their influences on human primary cells are very limited or under reported. The potential mechanisms by which allethrin and prallethrin modulates human primary cells, especially the molecular mechanisms or interconnectivity of autophagy-apoptosis, their clinical relevance in human subjects or patients are not well defined. In this current study, we've furnished the evidence that both allethrin and prallethrin user samples significantly induced Ccl2 mRNA expression, increased amount of reactive oxygen intermediate, inhibited membrane bound enzymes and altered membrane fluidity. Pyrethroid derivative users had induced levels of lipid peroxidation and induced binding activities of transcription factors(tfs) like CEBP-β and NF-AT. Pyrethroid derivatives induced autophagy, elicited intracellular Ca2+ concentration, calcineurin and regulated proapoptotic genes, DAPK1, Bim. Our current study presumably comprises the initial investigation of a very new mechanism of pyrethroid derivatives-moderated programed cell death in various cell sets or types, like human primary cells where-in this is a late event, is documented. Hence, current research-study might be significant in the various pyrethroid derivatives-allied hematological-related cancers and immunosuppressant or auto-immune disorders. In the foremost instance, we present data stating that pyrethroid derivatives induces multiple cell signaling cascades, like CEBP-β, NF-AT, ERK and MAPK having a role in autophagy thereby; synchronously effectively impact on the apoptosis, therefore causing hematological tumors and toxic or immune related disorders.

In Situ Hybridization as a Method to Examine Gene Regulatory Activity In Vivo.

Transcription factor-enhancer binding events are among the most well-studied protein-DNA interactions, allowing researchers to determine mechanisms of transcriptional activation or repression during development. While large-scale ChIP-sequence datasets, together with computational predictions and chromatin accessibility data, yield information on potential transcription factor binding activities, reporter gene assays provide measurable information on whether these binding activities are functional in particular cell types during development. Here, we present a detailed protocol to examine enhancer activity in Drosophila embryos using cloning, transgenesis, and in situ hybridization.

A leukemia-protective germline variant mediates chromatin module formation via transcription factor nucleation

Non-coding variants coordinate transcription factor (TF) binding and chromatin mark enrichment changes over regions spanning >100 kb. These molecularly coordinated regions are named “variable chromatin modules” (VCMs), providing a conceptual framework of how regulatory variation might shape complex traits. To better understand the molecular mechanisms underlying VCM formation, here, we mechanistically dissect a VCM-modulating noncoding variant that is associated with reduced chronic lymphocytic leukemia (CLL) predisposition and disease progression. This common, germline variant constitutes a 5-bp indel that controls the activity of an AXIN2 gene-linked VCM by creating a MEF2 binding site, which, upon binding, activates a super-enhancer-like regulatory element. This triggers a large change in TF binding activity and chromatin state at an enhancer cluster spanning >150 kb, coinciding with subtle, long-range chromatin compaction and robust AXIN2 up-regulation. Our results support a model in which the indel acts as an AXIN2 VCM-activating TF nucleation event, which modulates CLL pathology.

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.