Genome-wide analysis, identification, evolution and genomic organization of dehydration responsive element-binding (DREB) gene family in Solanum tuberosum

Background The dehydration responsive element-binding (DREB) gene family plays a crucial role as transcription regulators and enhances plant tolerance to abiotic stresses. Although the DREB gene family has been identified and characterized in many plants, knowledge about it in Solanum tuberosum (Potato) is limited. Results In the present study, StDREB gene family was comprehensively analyzed using bioinformatics approaches. We identified 66 StDREB genes through genome wide screening of the Potato genome based on the AP2 domain architecture and amino acid conservation analysis (Valine at position 14th). Phylogenetic analysis divided them into six distinct subgroups (A1–A6). The categorization of StDREB genes into six subgroups was further supported by gene structure and conserved motif analysis. Potato DREB genes were found to be distributed unevenly across 12 chromosomes. Gene duplication proved that StDREB genes experienced tandem and segmental duplication events which led to the expansion of the gene family. The Ka/Ks ratios of the orthologous pairs also demonstrated the StDREB genes were under strong purification selection in the course of evolution. Interspecies synteny analysis revealed 45 and 36 StDREB genes were orthologous to Arabidopsis and Solanum lycopersicum, respectively. Moreover, subcellular localization indicated that StDREB genes were predominantly located within the nucleus and the StDREB family’s major function was DNA binding according to gene ontology (GO) annotation. Conclusions This study provides a comprehensive and systematic understanding of precise molecular mechanism and functional characterization of StDREB genes in abiotic stress responses and will lead to improvement in Solanum tuberosum

A Computational Systems Analyses to Identify Biomarkers and Mechanistic Link in Psoriasis and Cutaneous Squamous Cell Carcinoma.

Psoriasis is the most common and chronic skin disease that affects individuals from every age group. The rate of psoriasis is increasing over the time in both developed and developing countries. Studies have revealed the possibility of association of psoriasis with skin cancers, particularly non-melanoma skin cancers (NMSC), which, include basal cell carcinoma and cutaneous squamous cell carcinoma (cSCC). There is a need to analyze the disease at molecular level to propose potential biomarkers and therapeutic targets in comparison to cSCC. Therefore, the second analyzed disease of this study is cSCC. It is the second most common prevalent skin cancer all over the world with the potential to metastasize and recur. There is an urge to validate the proposed biomarkers and discover new potential biomarkers as well. In order to achieve the goals and objectives of the study, microarray and RNA-sequencing data analyses were performed followed by network analysis. Afterwards, quantitative systems biology was implemented to analyze the results at a holistic level. The aim was to predict the molecular patterns that can lead psoriasis to cancer. The current study proposed potential biomarkers and therapeutic targets for psoriasis and cSCC. IL-17 signaling pathway is also identified as significant pathway in both diseases. Moreover, the current study proposed that autoimmune pathology, neutrophil recruitment, and immunity to extracellular pathogens are sensitive towards MAPKs (MAPK13 and MAPK14) and genes for AP-1 (FOSL1 and FOS). Therefore, these genes should be further studied in gene knock down based studies as they may play significant role in leading psoriasis towards cancer

Modelling and analysis of the complement system signalling pathways: roles of C3, C5a and pro-inflammatory cytokines in SARS-CoV-2 infection

The complement system is an essential part of innate immunity. It is activated by invading pathogens causing inflammation, opsonization, and lysis via complement anaphylatoxins, complement opsonin’s and membrane attack complex (MAC), respectively. However, in SARS-CoV-2 infection overactivation of complement system is causing cytokine storm leading to multiple organs damage. In this study, the René Thomas kinetic logic approach was used for the development of biological regulatory network (BRN) to model SARS-CoV-2 mediated complement system signalling pathways. Betweenness centrality analysis in cytoscape was adopted for the selection of the most biologically plausible states in state graph. Among the model results, in strongly connected components (SCCs) pro-inflammatory cytokines (PICyts) oscillatory behaviour between recurrent generation and downregulation was found as the main feature of SARS-CoV-2 infection. Diversion of trajectories from the SCCs leading toward hyper-inflammatory response was found in agreement with in vivo studies that overactive innate immunity response caused PICyts storm during SARS-CoV-2 infection. The complex of negative regulators FI, CR1 and DAF in the inhibition of complement peptide (C5a) and PICyts was found desirable to increase immune responses. In modelling role of MAC and PICyts in lowering of SARS-CoV-2 titre was found coherent with experimental studies. Intervention in upregulation of C5a and PICyts by C3 was found helpful in back-and-forth variation of signalling pattern linked with the levels of PICyts. Moreover, intervention in upregulation of PICyts by C5a was found productive in downregulation of all activating factors in the normal SCCs. However, the computational model predictions require experimental studies to be validated by exploring the activation role of C3 and C5a which could change levels of PICyts at various phases of SARS-CoV-2 infection

Structural insights and characterization of human Npas4 protein

Npas4 is an activity dependent transcription factor which is responsible for gearing the expression of target genes involved in neuro-transmission. Despite the importance of Npas4 in many neuronal diseases, the tertiary structure of Npas4 protein along with its physico-chemical properties is limited. In the current study, first we perfomed the phylogenetic analysis of Npas4 and determined the content of hydrophobic, flexible and order-disorder promoting amino acids. The protein binding regions, post-translational modifications and crystallization propensity of Npas4 were predicted through different in-silico methods. The three dimensional model of Npas4 was predicted through LOMET, SPARSKS-X, I-Tasser, RaptorX, MUSTER and Pyhre and the best model was selected on the basis of Ramachandran plot, PROSA, and Qmean scores. The best model was then subjected to further refinement though MODREFINER. Finally the interacting partners of Npas4 were identified through STRING database. The phylogenetic analysis showed the human Npas4 gene to be closely related to other primates such as chimpanzees, monkey, gibbon. The physiochemical properties of Npas4 showed that it is an intrinsically disordered protein with N-terminal ordered region. The post-translational modification analyses indicated absence of acetylation and mannosylation sites. Three potential phosphorylation sites (S108, T130 and T136) were found in PAS A domain whilst a single phosphorylation site (S273) was present in PAS B domain. The predicted tertiary structure of Npas4 showed that bHLH domain and PAS domain possess tertiary structures while the rest of the protein exhibited disorder property. Protein-protein interaction analysis revealed NPas4 interaction with various proteins which are mainly involved in nuclear trafficking of proteins to cytoplasm, activity regulated gene transcription and neurodevelopmental disorders. Moreover the analysis also highlighted the direct relation to proteins involved in promoting neuronal survival, plasticity and cAMP responsive element binding protein proteins. The current study helps in understanding the physicochemical properties and reveals the neuro-modulatory role of Npas4 in crucial pathways involved in neuronal survival and neural signalling hemostasis

Long non-coding RNAs and their targets as potential biomarkers in breast cancer.

Breast cancer is among the lethal types of cancer with a high mortality rate, globally. Its high prevalence can be controlled through improved analysis and identification of disease-specific biomarkers. Recently, long non-coding RNAs (lncRNAs) have been reported as key contributors of carcinogenesis and regulate various cellular pathways through post-transcriptional regulatory mechanisms. The specific aim of this study was to identify the novel interactions of aberrantly expressed genetic components in breast cancer by applying integrative analysis of publicly available expression profiles of both lncRNAs and mRNAs. Differential expression patterns were identified by comparing the breast cancer expression profiles of samples with controls. Significant co-expression networks were identified through WGCNA analysis. WGCNA is a systems biology approach used to elucidate the pattern of correlation between genes across microarray samples. It is also used to identify the highly correlated modules. The results obtained from this study revealed significantly differentially expressed and co-expressed lncRNAs and their cis- and trans-regulating mRNA targets which include RP11-108F13.2 targeting TAF5L, RPL23AP2 targeting CYP4F3, CYP4F8 and AL022324.2 targeting LRP5L, AL022324.3, and Z99916.3, respectively. Moreover, pathway analysis revealed the involvement of identified mRNAs and lncRNAs in major cell signalling pathways, and target mRNAs expression is also validated through cohort data. Thus, the identified lncRNAs and their target mRNAs represent novel biomarkers that could serve as potential therapeutics for breast cancer and their roles could also be further validated through wet labs to employ them as potential therapeutic targets in future

A Computational Systems Analyses to Identify Biomarkers and Mechanistic Link in Psoriasis and Cutaneous Squamous Cell Carcinoma

Psoriasis is the most common and chronic skin disease that affects individuals from every age group. The rate of psoriasis is increasing over the time in both developed and developing countries. Studies have revealed the possibility of association of psoriasis with skin cancers, particularly non-melanoma skin cancers (NMSC), which, include basal cell carcinoma and cutaneous squamous cell carcinoma (cSCC). There is a need to analyze the disease at molecular level to propose potential biomarkers and therapeutic targets in comparison to cSCC. Therefore, the second analyzed disease of this study is cSCC. It is the second most common prevalent skin cancer all over the world with the potential to metastasize and recur. There is an urge to validate the proposed biomarkers and discover new potential biomarkers as well. In order to achieve the goals and objectives of the study, microarray and RNA-sequencing data analyses were performed followed by network analysis. Afterwards, quantitative systems biology was implemented to analyze the results at a holistic level. The aim was to predict the molecular patterns that can lead psoriasis to cancer. The current study proposed potential biomarkers and therapeutic targets for psoriasis and cSCC. IL-17 signaling pathway is also identified as significant pathway in both diseases. Moreover, the current study proposed that autoimmune pathology, neutrophil recruitment, and immunity to extracellular pathogens are sensitive towards MAPKs (MAPK13 and MAPK14) and genes for AP-1 (FOSL1 and FOS). Therefore, these genes should be further studied in gene knock down based studies as they may play significant role in leading psoriasis towards cancer

Integrated Analysis of Microarray and RNA-Seq Data for the Identification of Hub Genes and Networks Involved in the Pancreatic Cancer

Pancreatic cancer (PaCa) is the seventh most fatal malignancy, with more than 90% mortality rate within the first year of diagnosis. Its treatment can be improved the identification of specific therapeutic targets and their relevant pathways. Therefore, the objective of this study is to identify cancer specific biomarkers, therapeutic targets, and their associated pathways involved in the PaCa progression. RNA-seq and microarray datasets were obtained from public repositories such as the European Bioinformatics Institute (EBI) and Gene Expression Omnibus (GEO) databases. Differential gene expression (DE) analysis of data was performed to identify significant differentially expressed genes (DEGs) in PaCa cells in comparison to the normal cells. Gene co-expression network analysis was performed to identify the modules co-expressed genes, which are strongly associated with PaCa and as well as the identification of hub genes in the modules. The key underlaying pathways were obtained from the enrichment analysis of hub genes and studied in the context of PaCa progression. The significant pathways, hub genes, and their expression profile were validated against The Cancer Genome Atlas (TCGA) data, and key biomarkers and therapeutic targets with hub genes were determined. Important hub genes identified included ITGA1, ITGA2, ITGB1, ITGB3, MET, LAMB1, VEGFA, PTK2, and TGFβ1. Enrichment analysis characterizes the involvement of hub genes in multiple pathways. Important ones that are determined are ECM–receptor interaction and focal adhesion pathways. The interaction of overexpressed surface proteins of these pathways with extracellular molecules initiates multiple signaling cascades including stress fiber and lamellipodia formation, PI3K-Akt, MAPK, JAK/STAT, and Wnt signaling pathways. Identified biomarkers may have a strong influence on the PaCa early stage development and progression. Further, analysis of these pathways and hub genes can help in the identification of putative therapeutic targets and development of effective therapies for PaCa

In silico evaluation of molecular virus–virus interactions taking place between Cotton leaf curl Kokhran virus- Burewala strain and Tomato leaf curl New Delhi virus

Background Cotton leaf curl disease (CLCuD) is a disease of cotton caused by begomoviruses, leading to a drastic loss in the annual yield of the crop. Pakistan has suffered two epidemics of this disease leading to the loss of billions in annual exports. The speculation that a third epidemic of CLCuD may result as consequence of the frequent occurrence of Tomato leaf curl New Delhi virus (ToLCNDV) and Cotton leaf curl Kokhran Virus-Burewala Strain (CLCuKoV-Bu) in CLCuD infected samples, demand that the interactions taking between the two viruses be properly evaluated. This study is designed to assess virus-virus interactions at the molecular level and determine the type of co-infection taking place. Methods Based on the amino acid sequences of the gene products of both CLCuKoV-Bu and ToLCNDV, protein structures were generated using different software, i.e., MODELLER, I-TASSER, QUARKS, LOMETS and RAPTORX. A consensus model for each protein was selected after model quality assessment using ERRAT, QMEANDisCo, PROCHECK Z-Score and Ramachandran plot analysis. The active and passive residues in the protein structures were identified using the CPORT server. Protein–Protein Docking was done using the HADDOCK webserver, and 169 Protein–Protein Interaction (PPIs) were performed between the proteins of the two viruses. The docked complexes were submitted to the PRODIGY server to identify the interacting residues between the complexes. The strongest interactions were determined based on the HADDOCK Score, Desolvation energy, Van der Waals Energy, Restraint Violation Energy, Electrostatic Energy, Buried Surface Area and Restraint Violation Energy, Binding Affinity and Dissociation constant (Kd). A total of 50 ns Molecular Dynamic simulations were performed on complexes that exhibited the strongest affinity in order to validate the stability of the complexes, and to remove any steric hindrances that may exist within the structures. Results Our results indicate significant interactions taking place between the proteins of the two viruses. Out of all the interactions, the strongest were observed between the Replication Initiation protein (Rep) of CLCuKoV-Bu with the Movement protein (MP), Nuclear Shuttle Protein (NSP) of ToLCNDV (DNA-B), while the weakest were seen between the Replication Enhancer protein (REn) of CLCuKoV-Bu with the REn protein of ToLCNDV. The residues identified to be taking a part in interaction belonged to domains having a pivotal role in the viral life cycle and pathogenicity. It maybe deduced that the two viruses exhibit antagonistic behavior towards each other, and the type of infection may be categorised as a type of Super Infection Exclusion (SIE) or homologous interference. However, further experimentation, in the form of transient expression analysis, is needed to confirm the nature of these interactions and increase our understanding of the direct interactions taking place between two viruses

Antibody designing against IIIabc junction (JIIIabc) of HCV IRES through affinity maturation; RNA-Antibody docking and interaction analysis.

Hepatitis C virus is a single-stranded RNA based virus which can cause chronic HCV and hepatocellular carcinoma. HCV genotype 3a has relatively higher rate of fibrosis progression, prevalence of steatosis and incidence of HCC. Despite HCVs variation in genomic sequence, the 5' untranslated region containing internal ribosome entry site (IRES) is highly conserved among all genotypes. It is responsible for translation and initiation of the viral protein. In present study, IRES was targeted by designing variants of reported antigen binding fragment (Fab) through affinity maturation approach. Affinity maturation strategy allowed the rational antibody designing with better biophysical properties and antibody-antigen binding interactions. Complementarity determining regions of reported Fab (wild type) were assessed and docked with IRES. Best generated model of Fab was selected and subjected to alanine scanning Three sets of insilico mutations for variants (V) designing were selected; single (1-71), double (a-j) and triple (I-X). Redocking of IRES-Fab variants consequently enabled the discovery of three variants exhibiting better docking score as compared to the wild type Fab. V1, V39 and V4 exhibited docking scores of -446.51, -446.52 and-446.29 kcal/mol respectively which is better as compared to the wild type Fab that exhibited the docking score of -351.23 kcal/mol. Variants exhibiting better docking score were screened for aggregation propensity by assessing the aggregation prone regions in Fab structure. Total A3D scores of wild type Fab, V1, V4 and V39 were predicted as -315.325, -312.727, -316.967 and -317.545 respectively. It is manifested that solubility of V4 and V39 is comparable to wild type Fab. In future, development and invitro assessment of these promising Fab HCV3 variants is aimed