Deciphering the miRNA-mRNA Interaction Landscape between Breast Cancer and Triple-Negative Breast Cancer: An Integrated Bioinformatics Approach.

Breast cancer (BC) is globally recognized as the second most prevalent form of cancer. It predominantly affects women and can be categorized into distinct types based on the overexpression of specific cancer receptors.The key receptors implicated in this context are the human epidermal growth factor receptor-2 (HER2), estrogen receptor (ER), and progesterone receptor (PR), alongside a particularly intricate subclass known as triple-negative breast cancer (TNBC). This subclassification is critical for the stratification of breast cancer and informs therapeutic decision-making processes. Due to a lack of therapeutic targets, such as growth factor receptors, TNBC is the most aggressive type. Hence, identifying targetable regulators such as miRNAs could pave the way for potential therapeutic interventions. To identify common differentially expressed mRNAs (DE-mRNAs) in BC, including TNBC, we leveraged two data sets from the GEO collection and The Cancer Genome Atlas (TCGA). Significant DE-mRNAs were identified through PPI, MCODE, CytoNCA, and CytoHubba analyses. Following this, miRNAs were predicted using mirDIP. We utilized GSE42568, GSE185645, and TCGA and identified 159 common DE-mRNAs. Using Cytoscape plug-ins, we identified the 10 most significant DE-mRNAs in BC. Using mirDIP, target miRNAs for 10 DE-mRNAs were identified. We conducted an advanced analysis on the TNBC GEO data set (GSE45498) to corroborate the significance of shared DE-mRNAs and DE-miRNAs in TNBC. We identified four downregulated DE-miRNAs, including hsa-miR-802, hsa-miR-1258, hsa-miR-548a-3p, and hsa-miR-2053, significantly associated with TNBC. Our study revealed significant miRNA-mRNA interactions, specifically hsa-miR-802/MELK, hsa-miR-1258/NCAPG, miR-548a-3p/CCNA2, and hsa-miR-2053/NUSAP1, in both BC and TNBC. The observed downregulation of hsa-miR-548a-3p is associated with diminished survival rates in BC patients, emphasizing their potential utility as prognostic indicators. Furthermore, the differential expression of mRNAs, including CCNB2, UBE2C, MELK, and KIF2C, correlates with reduced survival outcomes, signifying their critical role as potential targets for therapeutic intervention in both BC and TNBC. These findings highlight specific regulatory mechanisms that are potentially crucial for understanding and treating these cancer types.

Deciphering the Impact of Rare Missense Variants in EGFR-TKI-Resistant Non-Small-Cell Lung Cancer through Whole Exome Sequencing: A Computational Approach.

Targeted therapy revolutionizes the treatment of non-small-cell lung cancer (NSCLC), harboring molecular change. Epidermal growth factor receptor(EGFR) mutations play a crucial role in the development of NSCLC, serving as a pivotal factor in its pathogenesis. We elucidated the mechanisms of resistance and potential therapeutic strategies in NSCLC resistant to the EGFR-tyrosine kinase inhibitor (EGFR-TKI). This is achieved by identifying rare missense variants through whole exome sequencing (WES). The goal is to enhance our understanding, identify biomarkers, and lay the groundwork for targeted interventions, thereby offering hope for an improved NSCLC treatment landscape. We conducted WES analysis on 16 NSCLC samples with EGFR-TKI-resistant NSCLC obtained from SRA-NCBI (PRJEB50602) to reveal genomic profiles within the EGFR-TKI. Our findings showed that 48% of the variants were missense, and after filtering with the Ensembl variant effect predictor, 53 rare missense variants in 23 genes were identified as highly deleterious. Further examination using pathogenic tools like PredictSNP revealed 12 deleterious rare missense variants in 7 genes: ZNF717, PSPH, ESRRA, SEMA3G, PTPN7, CAVIN4, and MYBBP1A. Molecular dynamics simulation (MDS) suggested that the L385P variant alters the structural flexibility of ESRRA, potentially leading to unfolding of ERRα proteins. This could impact their function and alter ERRα expression. These insights from MDS enhance our understanding of the structural and dynamic consequences of the L385P ESRRA variant and provide valuable implications for subsequent therapeutic considerations and targeted interventions.

Comparative Atomistic Insights on Apo and ATP-I1171N/S/T in Nonsmall-Cell Lung Cancer.

Anaplastic lymphoma kinase (ALK) rearrangements occur in about 5% of nonsmall cell lung cancer (NSCLC) patients. Despite being first recognized as EML4-ALK, fusions with several additional genes have been identified, all of which cause constitutive activation of the ALK kinase and subsequently lead to tumor development. ALK inhibitors first-line crizotinib, second-line ceritinib, and alectinib are effective against NSCLC patients with these rearrangements. Patients progressing on crizotinib had various mutations in the ALK kinase domain. ALK fusion proteins are activated by oligomerization through the fusion partner, which leads to the autophosphorylation of the kinase's domain and consequent downstream activation. The proposed computational study focuses on understanding the activation mechanism of ALK and ATP binding of wild-type (WT) and I1171N/S/T mutations. We analyzed the conformational change of ALK I1171N/S/T mutations and ATP binding using molecular docking and molecular dynamics simulation approaches. According to principal component analysis and free energy landscape, it is clear that I1171N/S/T mutations in Apo and ATP showed different energy minima/unstable structures compared to WT-Apo. The results revealed that I1171N/S/T mutations and ATP binding significantly supported a change toward an active-state conformation, whereas WT-Apo remained inactive. We demonstrated that I1171N/S/T mutations are persistent in an active state and independent of ATP. The I1171S/T mutations showed greater intermolecular H-bonds with ATP than WT-ATP. The molecular mechanics Poisson-Boltzmann surface area analysis revealed that the I1171N/S/T mutation binding energy was similar to that of WT-ATP. This study shows that I1171N/S/T can form stable bonds with ATP and may contribute to a constitutively active kinase. Based on the Y1278-C1097 H-bond and E1167-K1150 salt bridge interaction, I1171N strongly promotes the constitutively active kinase independent of ATP. This structural mechanism study will aid in understanding the oncogenic activity of ALK and the basis for improving the ALK inhibitors.

The targeted next-generation sequence revealed SMAD4, AKT1, and TP53 mutations from circulating cell-free DNA of breast cancer and its effect on protein structure - A computational approach.

Breast cancer biomarkers that detect marginally advanced stages are still challenging. The detection of specific abnormalities, targeted therapy selection, prognosis, and monitoring of treatment effectiveness over time are all made possible by circulating free DNA (cfDNA) analysis. The proposed study will detect specific genetic abnormalities from the plasma cfDNA of a female breast cancer patient by sequencing a cancer-related gene panel (MGM455 - Oncotrack Ultima), including 56 theranostic genes (SNVs and small INDELs). Initially, we determined the pathogenicity of the observed mutations using PredictSNP, iStable, Align-GVGD, and ConSurf servers. As a next step, molecular dynamics (MD) was implemented to determine the functional significance of SMAD4 mutation (V465M). Lastly, the mutant gene relationships were examined using the Cytoscape plug-in GeneMANIA. Using ClueGO, we determined the gene's functional enrichment and integrative analysis. The structural characteristics of SMAD4 V465M protein by MD simulation analysis further demonstrated that the mutation was deleterious. The simulation showed that the native structure was more significantly altered by the SMAD4 (V465M) mutation. Our findings suggest that SMAD4 V465M mutation might be significantly associated with breast cancer, and other patient-found mutations (AKT1-E17K and TP53-R175H) are synergistically involved in the process of SMAD4 translocate to nuclease, which affects the target gene translation. Therefore, this combination of gene mutations could alter the TGF-ß signaling pathway in BC. We further proposed that the SMAD4 protein loss may contribute to an aggressive phenotype by inhibiting the TGF-ß signaling pathway. Thus, breast cancer's SMAD4 (V465M) mutation might increase their invasive and metastatic capabilities.Communicated by Ramaswamy H. Sarma.

Controlling cell proliferation by targeting cyclin-dependent kinase 6 using drug repurposing approach.

Cyclin-dependent kinase 6 (CDK6) is an essential kinase in cell cycle progression, which is a viable target for inhibitors in various malignancies, including breast cancer. This study aimed to virtually screen efficient compounds as new leads in treating breast cancer using a drug repurposing approach. Apoptosis regulatory compounds were taken from the seleckchem database. Molecular docking experiments were carried out in the presence of abemaciclib, a routinely used FDA drug. Compared to conventional drugs, the two compounds demonstrated a higher binding affinity for CDK6. Compounds (N-benzyl-6-[(4-hydroxyphenyl)methyl]-8-(naphthalen-1-ylmethyl)-4,7-dioxo-3,6,9,9a-tetrahydro-2H-pyrazino[1,2-a]pyrimidine-1-carboxamide) and (1'-[4-[1-(4-fluorophenyl)indol-3-yl]butyl]spiro[1H-2-benzofuran-3,4'-piperidine]) were discovered to have an inhibitory effect against CDK6 at -8.49 and -6.78kcal/mol, respectively, compared to -8.09kcal/mol of the control molecule, the interacting residues of these two new compounds were found to fall within the binding site of the CDK6 molecule. Both compounds exhibited equal ADME features compared with abemaciclib and would be well distributed and metabolized by the body with an appropriate druglikeness range. Lastly, molecular dynamics was initiated for 200ns for the selected potent inhibitors and abemaciclib as complexed with CDK6. The RMSD, RMSF, Rg, H-Bond interactions, SASA, PCA, FEL, and MM/PBSA analysis were performed for the complexes to assess the stability, fluctuations, radius of gyration, hydrogen bond interaction, solvent accessibility, essential dynamics, free energy landscape, and MM/PBSA. The selected two compounds are small molecules in the appropriate druglikeness range. The results observed in molecular docking and molecular dynamics simulations were most promising for two compounds, suggesting their potent inhibitory effect against CDK6. We propose that these candidate compounds can undergo in vitro validation and in vivo testing for their further use against cancer.

HPV and molecular mimicry in systemic lupus erythematosus and an impact of compiling B-cell epitopes and MHC-class II binding profiles with in silico evidence.

Epidemiological link between HPV and SLE is evolving. The possibility of HPV infection-induced molecular mimicry and systemic lupus erythematosus (SLE) was elucidated through detailed in silico analyses. Conserved regions in the structural protein sequences of high-risk HPV types were inferred, and sequence homologies between viral and human peptides were identified to delineate proteins implicated in SLE. B-cell epitopes and MHC-class II binding were compiled using Immune Epitope Database and ProPred II analysis tool. Molecular modeling and molecular dynamics/simulation (MDS) were performed using AutoDock Vina and GROMACS, respectively. Sequence alignment revealed 32 conserved regions, and 27/32 viral peptides showed varying similarities to human peptides, rich in B-cell epitopes with superior accessibility, high hydrophilicity, antigenicity and disposition to bind many class-II HLA alleles. Molecular docking of 13 viral peptides homologous (100%) to human peptides implicated in SLE showed that VIR-PEP1 (QLFNKPYWL) and VIR-PEP2 (DTYRFVTS) exhibited higher binding affinities than corresponding human peptides to SLE predisposing HLA-DRB1 allele. MDS of these peptides showed that the viral peptides had superior folding, compactness, and a higher number of hydrogen bonds than human peptides throughout the simulation period. SASA analysis revealed that the VIR-PEP1&2 fluctuated less frequently than corresponding human peptides. MM-PBSA revealed that the VIR-PEP2 complex exhibited higher binding energy than the human peptide complex. This suggests that highly conserved structural peptides of high-risk HPV types homologous to human peptides could compete and bind avidly to the HLA allele associated with SLE and predispose HPV-infected individuals to SLE through molecular mimicry.Communicated by Ramaswamy H. Sarma.

In silico analysis revealed the potential circRNA-miRNA-mRNA regulative network of non-small cell lung cancer (NSCLC).

BACKGROUND: The primary source of death in the world is non-small cell lung cancer (NSCLC). However, NSCLCs pathophysiology is still not completely understood. The current work sought to study the differential expression of mRNAs involved in NSCLC and their interactions with miRNAs and circRNAs. METHODS: We utilized three microarray datasets (GSE21933, GSE27262, and GSE33532) from the GEO NCBI database to identify the differentially expressed genes (DEGs) in NSCLC. We employed DAVID Functional annotation tool to investigate the underlying GO biological process, molecular functions, and KEGG pathways involved in NSCLC. We performed the Protein-protein interaction (PPI) network, MCODE, and CytoHubba analysis from Cytoscape software to identify the significant DEGs in NSCLC. We utilized miRnet to anticipate and build interaction between miRNAs and mRNAs in NSCLC and ENCORI to predict the miRNA-circRNA relationships and build the ceRNA regulatory network. Finally, we executed the gene expression and Kaplan-Meier survival analysis to validate the significant DEGs in the ceRNA network utilizing TCGA NSCLC and GEPIA data. RESULTS: We revealed a total of 156 overlapped DEGs (47 upregulated and 109 downregulated genes) in NSCLC. The PPI network, MCODE, and CytoHubba analysis revealed 12 hub genes (cdkn3, rrm2, ccnb1, aurka, nuf2, tyms, kif11, hmmr, ccnb2, nek2, anln, and birc5) that are associated with NSCLC. We identified that these 12 genes encode 12 mRNAs that are strongly linked with 8 miRNAs, and further, we revealed that 1 circRNA was associated with this 5 miRNA. We constructed the ceRNAs network that contained 1circRNA-5miRNAs-7mRNAs. The expression of these seven significant genes in LUAD & LUSC (NSCLC) was considerably higher in the TCGA database than in normal tissues. Kaplan-Meier survival plot reveals that increased expression of these hub genes was related to a poor survival rate in LUAD. CONCLUSION: Overall, we developed a circRNA-miRNA-mRNA regulation network to study the probable mechanism of NSCLC.

A computational examination of the therapeutic advantages of fourth-generation ALK inhibitors TPX-0131 and repotrectinib over third-generation lorlatinib for NSCLC with ALK F1174C/L/V mutations.

Background: In non-small-cell lung cancer (NSCLC), a pivotal factor in promoting cancer development is the rearrangement in the anaplastic lymphoma kinase ALK gene, resulting in elevated ALK protein expression. F1174C/L/V is the acquired secondary resistant mutation in ALK. Significant survival improvements have been seen while tyrosine kinase inhibitors specifically target ALK. Nevertheless, the emergence of drug resistance hinders the clinical effectiveness of these drugs. Objective: This research sought to find the binding affinity/inhibitory effects of the existing drug lorlatinib (LOR) and upcoming TPX-0131 (zotizalkib/TPX) and repotrectinib (TPX-0005/REP) inhibitors against ALK F1174C/L/V mutations using computational approaches to identify potential strategies over resistance. Methods: We conducted molecular docking, molecular dynamics simulation, and MMPBSA calculations to investigate how compact macrocyclic inhibitors, such as TPX-0131 and repotrectinib, fit within the ATP-binding boundary and differ from LOR. Results: Our results demonstrated that TPX-0131 and repotrectinib contributed to higher binding energy in F1174C and F1174L mutations than LOR. Repotrectinib showed greater binding energy in the F1174V mutation, whereas LOR and TPX-0131 exhibited similar binding energy. However, all three inhibitors showed significant binding energy toward F1174C/L/V mutations found in NSCLC. Conclusion: This comparative study of the potential binding effects of fourth-generation inhibitors TPX-0131 and repotrectinib and third-generation inhibitor LOR for ALK F1174C/L/V mutations revealed the atomistic insights of the binding mechanism. These computational findings enable us to carry out further research for the clinical implementation of fourth-generation ALK inhibitors on ALK-positive NSCLC.

Unraveling the Structural Changes in the DNA-Binding Region of Tumor Protein p53 (TP53) upon Hotspot Mutation p53 Arg248 by Comparative Computational Approach.

The vital tissue homeostasis regulator p53 forms a tetramer when it binds to DNA and regulates the genes that mediate essential biological processes such as cell-cycle arrest, senescence, DNA repair, and apoptosis. Missense mutations in the core DNA-binding domain (109-292) simultaneously cause the loss of p53 tumor suppressor function and accumulation of the mutant p53 proteins that are carcinogenic. The most common p53 hotspot mutation at codon 248 in the DNA-binding region, where arginine (R) is substituted by tryptophan (W), glycine (G), leucine (L), proline (P), and glutamine (Q), is reported in various cancers. However, it is unclear how the p53 Arg248 mutation with distinct amino acid substitution affects the structure, function, and DNA binding affinity. Here, we characterized the pathogenicity and protein stability of p53 hotspot mutations at codon 248 using computational tools PredictSNP, Align GVGD, HOPE, ConSurf, and iStable. We found R248W, R248G, and R248P mutations highly deleterious and destabilizing. Further, we subjected all five R248 mutant-p53-DNA and wt-p53-DNA complexes to molecular dynamics simulation to investigate the structural stability and DNA binding affinity. From the MD simulation analysis, we observed increased RMSD, RMSF, and Rg values and decreased protein-DNA intermolecular hydrogen bonds in the R248-p53-DNA than the wt-p53-DNA complexes. Likewise, due to high SASA values, we observed the shrinkage of proteins in R248W, R248G, and R248P mutant-p53-DNA complexes. Compared to other mutant p53-DNA complexes, the R248W, R248G, and R248P mutant-p53-DNA complexes showed more structural alteration. MM-PBSA analysis showed decreased binding energies with DNA in all five R248-p53-DNA mutants than the wt-p53-DNA complexes. Henceforth, we conclude that the amino acid substitution of Arginine with the other five amino acids at codon 248 reduces the p53 protein's affinity for DNA and may disrupt cell division, resulting in a gain of p53 function. The proposed study influences the development of rationally designed molecular-targeted treatments that improve p53-based therapeutic outcomes in cancer.

Integrative ontology and pathway-based approach identifies distinct molecular signatures in transcriptomes of esophageal squamous cell carcinoma.

Esophageal squamous cell carcinoma (ESCC) remains a serious concern globally due to many factors that including late diagnosis, lack of an ideal biomarker for diagnosis and prognosis, and high rate of mortality. In this study, we aimed to identify the essential dysregulated genes and molecular signatures associated with the progression and development of ESCC. The dataset with 15 ESCCs and the 15 adjacent normal tissue samples from the surrounding histopathologically tumor-free mucosa was selected. We applied bioinformatics pipelines including various topological parameters from MCODE, CytoNCA, and cytoHubba to prioritize the most significantly associated DEGs with ESCC. We performed functional enrichment annotation for the identified DEGs using DAVID and MetaCore™ GeneGo platforms. Furthermore, we validated the essential core genes in TCGA and GTEx datasets between the normal mucosa and ESCC for their expression levels. These DEGs were primarily enriched in positive regulation of transferase activity, negative regulation of organelle organization, cell cycle mitosis/S-phase transition, spindle organization/assembly, development, and regulation of angiogenesis. Subsequently, the DEGs were associated with the pathways such as oocyte meiosis, cell cycle, and DNA replication. Our study identified the eight-core genes (AURKA, AURKB, MCM2, CDC20, TPX2, PLK1, FOXM1, and MCM7) that are highly expressed among the ESCC, and TCGA dataset. The multigene comparison and principal component analysis resulted in elevated signals for the AURKA, MCM2, CDC20, TPX2, PLK1, and FOXM1. Overall, our study reported GO profiles and molecular signatures that might help researchers to grasp the pathological mechanisms underlying ESCC development and eventually provide novel therapeutic and diagnostic strategies.

Whole-exome sequencing analysis of NSCLC reveals the pathogenic missense variants from cancer-associated genes.

BACKGROUND: Non-small-cell lung cancer (NSCLC) is the most common type of lung cancer. NSCLC accounts for 84% of all lung cancer cases. In recent years, advances in pathway understanding, methods for discovering novel genetic biomarkers, and new drugs designed to inhibit the signaling cascades have enabled clinicians to personalize therapy for NSCLC. OBJECTIVES: The primary aim of this study is to identify the genes associated with NSCLC that harbor pathogenic variants that could be causative for NSCLC. The second aim is to investigate their roles in different pathways that lead to NSCLC. METHODS: We examined exome-sequencing datasets from 54 NSCLC patients to characterize the variants associated with NSCLC. RESULTS: Our findings revealed that 17 variants in 14 genes were considered highly pathogenic, including CDKN2A, ERBB2, FOXP1, IDH1, JAK3, KMT2D, K-Ras, MSH3, MSH6, POLE, RNF43, TCF7L2, TP53, and TSC1. Gene set enrichment analysis revealed the involvement of transmembrane receptor protein tyrosine kinase activity, protein binding, ATP binding, phosphatidylinositol-4,5-bisphosphate 3-kinase, and Ras guanyl-nucleotide exchange factor activity. Pathway analysis of these genes yielded different cancer-related pathways, including colorectal, prostate, endometrial, pancreatic, PI3K-Akt signaling pathways, and signaling pathways regulating pluripotency of stem cells. Module 1 from protein-protein interactions (PPIs) identified genes that harbor pathogenic SNPs. Three of the most deleterious SNPs are ERBB2 (rs1196929947), K-Ras (rs121913529), and POLE (rs751425952). Interestingly, one patient has a pathogenic K-Ras variant (rs121913529) co-occurred with the missense variant (rs752054698) inTSC1 gene. CONCLUSION: This study maps highly pathogenic variants associated with NSCLC and investigates their contributions to the pathogenesis of NSCLC. This study sheds light on the potential applications of precision medicine in patients with NSCLC.

ZFTA (Zinc Finger Translocation Associated) Fusion in Supratentorial Ependymomas: Low Prevalence in South Asians and No Correlation with Survival.

BACKGROUND: Supratentorial ependymomas (STEs) are an aggressive group of ependymomas, topographically distinct from their posterior fossa and spinal counterparts. Zinc finger translocation associated (ZFTA) fusion-positive cases have been reported to account for the majority of STEs, although data on its association with poorer outcomes are inconsistent. MATERIALS AND METHODS: We assessed the prevalence of the ZFTA fusion by reverse-transcription polymerase chain reaction and fluorescence in situ hybridization in a cohort of 61 patients (68 samples) with STE. Our primary outcome was to determine the role of the ZFTA fusion on progression-free and overall survival of patients with STE. Our secondary objectives were to assess the impact of ZFTA fusion on nuclear factor (NF)-kB pathway signaling via surrogate markers of this pathway, namely COX-2, CCND1, and L1 cell adhesion molecule. RESULTS: ZFTA fusion was noted in 21.3% of STEs in our cohort. The presence of this rearrangement did not significantly impact the progression-free or overall survival of patients with STEs and was not associated with upregulation of markers of the NF-kB pathway. Only gross total resection was significantly associated with better progression-free survival. CONCLUSIONS: In contradiction to previous reports from across the world, the ZFTA fusion is far less prevalent among our population. It does not appear to drive NF-kB signaling or significantly affect outcomes. Gross total resection must be attempted in all cases of STE and adjuvant radiation and/or chemotherapy employed when gross total resection is not achieved.