Computational assessment of MCM2 transcriptional expression and identification of the prognostic biomarker for human breast cancer

Pixantrone as a novel MCM2 inhibitor for ovarian cancer treatment

BACKGROUND

Mini-chromosome maintenance protein 2 (MCM2) is a potential target for the development of cancer therapeutics. However, small molecule inhibitors targeting MCM2 need further investigation.

METHODS

Molecular dynamics simulation was performed to identify active pockets in the MCM2 protein structure (6EYC). The active pocket was used as a docking model to discover MCM2 inhibitors by using structure-based virtual screening and surface plasmon resonance (SPR) assay. Furthermore, the efficacy of pixantrone targeting MCM2 in ovarian cancer was evaluated in vitro and in vivo.

RESULTS

Pixantrone was identified as a novel inhibitor of MCM2 by virtual screening. SPR binding affinity analysis confirmed the direct binding of pixantrone to MCM2 protein. Pixantrone significantly reduced the viability of ovarian cancer cells A2780 and SKOV3 in a dose- and time-dependent manner. In addition, pixantrone inhibited DNA replication, and induced cell cycle arrest and apoptosis in ovarian cancer cells via targeting MCM2. Knockdown of MCM2 could attenuate the inhibitory activity of pixantrone in ovarian cancer cells. Furthermore, pixantrone significantly suppressed ovarian cancer growth in the A2780 cell xenograft mouse model and showed favorable safety.

CONCLUSION

These findings suggest that pixantrone may be a promising drug for ovarian cancer patients by targeting MCM2 in the clinic.

Pharmacokinetics and drug-likeness of anti-cancer traditional Chinese medicine: molecular docking and molecular dynamics simulation study

MCM7 (Minichromosome Maintenance Complex Component 7) is a component of the DNA replication licensing factor, which controls DNA replication. The MCM7 protein is linked to tumor cell proliferation and has a function in the development of several human cancers. Several types of cancer may be treated by inhibiting the protein, as it is strongly produced throughout this process. Significantly, Traditional Chinese Medicine (TCM), which has a long history of clinical adjuvant use against cancer, is rapidly gaining traction as a valuable medical resource for the development of novel cancer therapies, including immunotherapy. Therefore, the goal of the research was to find small molecular therapeutic candidates against the MCM7 protein that may be used to treat human cancers. A computational-based virtual screening of 36,000 natural TCM libraries is carried out for this goal using a molecular docking and dynamic simulation technique. Thereby, ∼8 novel potent compounds i.e., ZINC85542762, ZINC95911541, ZINC85542617, ZINC85542646, ZINC85592446, ZINC85568676, ZINC85531303, and ZINC95914464 were successfully shortlisted, each having the capacity to penetrate the cell as potent inhibitors for MCM7 to curb this disorder. These selected compounds were found to have high binding affinities compared to the reference (AGS compound) i.e. < -11.0 kcal/mol. ADMET and pharmacological properties showed that none of these 8 compounds poses any toxic property (carcinogenicity) and have anti-metastatic, and anticancer activity. Additionally, MD simulations were run to assess the compounds' stability and dynamic behavior with the MCM7 complex for about 100 ns. Finally, ZINC95914464, ZINC95911541, ZINC85568676, ZINC85592446, ZINC85531303, and ZINC85542646 are identified as highly stable within the complex throughout the 100 ns simulations. Moreover, the results of binding free energy suggested that the selected virtual hits significantly bind to the MCM7 which implied these compounds may act as a potential MCM7 inhibitor. However, in vitro testing protocols are required to further support these results. Further, assessment through various lab-based trial methods can assist with deciding the action of the compound that will give options in contrast to human cancer immunotherapy.Communicated by Ramaswamy H. Sarma.

Unveiling promising breast cancer biomarkers: an integrative approach combining bioinformatics analysis and experimental verification

BACKGROUND

Breast cancer remains a significant health challenge worldwide, necessitating the identification of reliable biomarkers for early detection, accurate prognosis, and targeted therapy.

MATERIALS AND METHODS

Breast cancer RNA expression data from the TCGA database were analyzed to identify differentially expressed genes (DEGs). The top 500 up-regulated DEGs were selected for further investigation using random forest analysis to identify important genes. These genes were evaluated based on their potential as diagnostic biomarkers, their overexpression in breast cancer tissues, and their low median expression in normal female tissues. Various validation methods, including online tools and quantitative Real-Time PCR (qRT-PCR), were used to confirm the potential of the identified genes as breast cancer biomarkers.

RESULTS

The study identified four overexpressed genes (CACNG4, PKMYT1, EPYC, and CHRNA6) among 100 genes with higher importance scores. qRT-PCR analysis confirmed the significant upregulation of these genes in breast cancer patients compared to normal samples.

CONCLUSIONS

These findings suggest that CACNG4, PKMYT1, EPYC, and CHRNA6 may serve as valuable biomarkers for breast cancer diagnosis, and PKMYT1 may also have prognostic significance. Furthermore, CACNG4, CHRNA6, and PKMYT1 show promise as potential therapeutic targets. These findings have the potential to advance diagnostic methods and therapeutic approaches for breast cancer.

Multi-omics prognostic signatures of IPO11 mRNA expression and clinical outcomes in colorectal cancer using bioinformatics approaches

The most prevalent malignant illness of the gastrointestinal system, colorectal cancer, is the third most prevalent cancer in males and the second most prevalent cancer in women. Importin-11 is a protein that acts as a regulator of cancer cell proliferation in colorectal tumours by conveying β -catenin to the cell nucleus. However, the IPO11 gene was found to encode a protein called Importin-11, which functions as a nucleus importer for the cell. As a result, preventing β -catenin from entering the nucleus requires blocking Importin-11. As a result, we conducted a multi-omics investigation to assess IPO11 gene potential as a therapeutic biomarker for human colorectal cancer (CC). Oncomine, GEPIA2, immunohisto-chemistry, and UALCAN databases were used to analyses the mRNA expression profiles of IPO11 in CC. The investigation has yielded clear evidence of the increase of IPO11 expression in CC subtypes, as indicated by the data acquired. Analysing CC research from the cBioPortal database, the study discovered three new missense mutations in the importin-11 protein sequence at a frequency of 0.00-1.50% copy number changes. Additionally, the Kaplan-Meier plots demonstrated a strong connection concerning IPO11 downregulation and a poorer CC patient survival rate. The co-expressed gene profile of IPO11 was likewise associated with the onset of CC. IPO11 co-expressed gene profile was also linked to CC development. Moreover, the correlation analysis using bc-GenExMiner and the UCSC Xena server identified KIF2A as the most positively co-expressed gene. The study found that KIF2A and its co-expressed genes were involved in a wide variety of cancer progression pathways using the Enrichr database. Cumulatively, this result will not only provide new information about the expression of IPO11 associated with CC progression and patient survival, but could also serve as a therapeutic biomarker for treating CC in a significant and worthwhile manner.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13755-023-00259-2.

Identification of natural antiviral drug candidates against Tilapia Lake Virus: Computational drug design approaches

Tilapia Lake Virus (TiLV) is a disease that affects tilapia fish, causing a high rate of sudden death at any stage in their life cycle. Unfortunately, there are currently no effective antiviral drugs or vaccines to prevent or control the progression of this disease. Researchers have discovered that the CRM1 protein plays a critical function in the development and spreading of animal viruses. By inhibiting CRM1, the virus's spread in commercial fish farms can be suppressed. With this in mind, this study intended to identify potential antiviral drugs from two different tropical mangrove plants from tropical regions: Heritiera fomes and Ceriops candolleana. To identify promising compounds that target the CRM1 protein, a computer-aided drug discovery approach is employed containing molecular docking, ADME (absorption, distribution, metabolism and excretion) analysis, toxicity assessment as well as molecular dynamics (MD) simulation. To estimate binding affinities of all phytochemicals, molecular docking is used and the top three candidate compounds with the highest docking scores were selected, which are CID107876 (-8.3 Kcal/mol), CID12795736 (-8.2 Kcal/mol), and CID12303662 (-7.9 Kcal/mol). We also evaluated the ADME and toxicity properties of these compounds. Finally, MD simulation was conducted to analyze the stability of the protein-ligand complex structures and confirm the suitability of these compounds. The computational study demonstrated that the phytochemicals found in H. fomes and C. candolleana could potentially serve as important inhibitors of TiLV, offering practical utility. However, further in vivo investigations are necessary to investigate and potentially confirm the effectiveness of these compounds as antiviral drugs against the virus TiLV.

MCM-2 Levels as a Potential Biomarker for Predicting High-Risk Breast Cancer Patients According to TAILORx Classification

Background

The minichromosome maintenance protein-2 (MCM-2) is a more sensitive proliferation marker than Ki-67. This study aimed to evaluate the relationship between MCM-2 and Oncotype DX recurrence score (ODX-RS) and determine an MCM-2 cutoff value in high-risk patients according to TAILORx risk categorization.

Methods

Hormone receptor (HR) positive HER-2 negative early-stage breast cancer patients (pT1-2, pN0-N1, M0) who had ODX-RS were included in the study. According to the TAILORx trial, patients were divided into two groups with high (ODX-RS ≥26) and low risk (ODX-RS <26) in terms of ODX-RS. Formalin-fixed-paraffin-embedded tissues of patients were re-evaluated, and 3 µm sections were prepared for MCM-2 immuno-histochemical staining. The relationship between ODX-RS and the percentage of MCM-2 staining was evaluated in two groups. The ROC curve analysis was performed to determine the MCM-2 cut-off value for the TAILORx high-risk group (ODX-RS ≥26).

Results

The mean MCM-2 value was significantly higher in the high-risk group [(60.2 ± 11.2 vs 34.4 ± 13.8, p < 0.001)]. In the multivariate analysis, MCM-2 (OR: 1.27, 95% CI: 1.08-1.49, p = 0.003) and progesterone receptor (PR) levels ≤10% (OR: 60.9, 95% CI: 4.1-89.7, p = 0.003) were found to be independent factors indicating a high-risk group. A one-unit increase in MCM-2 level increased the likelihood of being in the high-risk group by 1.27 times. In the ROC curve analysis, the optimal MCM-2 cut-off level was 50 (AUC: 0.921, sensitivity: 86.7%, specificity: 96.0%, p < 0.001).

Conclusion

Our study is the first study in the literature to investigate the relationship between ODX-RS and MCM-2 levels in HR-positive HER-2 negative early breast-cancer patients. In this study, MCM-2 was an independent risk factor in identifying high-risk patients according to TAILORx risk classification. MCM 2 cut-off value (50) may help the decision on adjuvant chemotherapy in patients where the Oncotype DX test cannot be performed.

Marine-derived sea urchin compounds as potential anti-cancer drug candidate against colorectal cancer: In silico and in vitro studies

Sea urchin-derived compounds are potential candidates for the development of effective drugs for the treatment of cancer diseases. In this study, 19 compounds derived from sea urchin (Diadema savignyi) were used to treat colorectal cancer using the HCT116 cell line. However, molecular docking, ADME (absorption, distribution, metabolism, and excretion), toxicity, molecular dynamic (MD) simulation, and molecular mechanics generalized Born surface area (MM-GBSA) were used to confirm the ligand-protein interaction. Interactions of Importin-11 receptor with sea urchin compounds reveal that four compounds have higher binding affinities (ranging from -8.6 to -7.1 kcal/mol). In vitro testing revealed that the CID 6432458 compound was effective (docking score of -8.6 kcal/mol) against the HCT116 cell line. The cytotoxicity of HCT116 has been documented, with an IC50 value of 6.885 ± 4. MTT assay, apoptosis analysis, and cell cycle assay were utilized to examine cell death in colorectal cancer. In the MTT experiment, 15 µM and 20 µM dosages were associated with 77% cell death; however, flow cytometry analysis using the IC50 value revealed that the selected chemical induced greater apoptosis in the HCT116 cell line (58.5%). The gene expression data revealed that the apoptotic gene BAX is expressed at a higher level than the BCL-2 gene. The IPO11 gene was downregulated during treatment. In the experiment involving the cell cycle, the S phase for the 30  µM dose showed 75.1% apoptosis, which was greater than the other concentrations used alone. These in silico and in vitro analysis will not only provide new information about Importin-11 receptor and insight into colorectal cancer but will also facilitate the development of natural compounds in a significant and worthwhile manner.

Impaired histone inheritance promotes tumor progression

Faithful inheritance of parental histones is essential to maintain epigenetic information and cellular identity during cell division. Parental histones are evenly deposited onto the replicating DNA of sister chromatids in a process dependent on the MCM2 subunit of DNA helicase. However, the impact of aberrant parental histone partition on human disease such as cancer is largely unknown. In this study, we construct a model of impaired histone inheritance by introducing MCM2-2A mutation (defective in parental histone binding) in MCF-7 breast cancer cells. The resulting impaired histone inheritance reprograms the histone modification landscapes of progeny cells, especially the repressive histone mark H3K27me3. Lower H3K27me3 levels derepress the expression of genes associated with development, cell proliferation, and epithelial to mesenchymal transition. These epigenetic changes confer fitness advantages to some newly emerged subclones and consequently promote tumor growth and metastasis after orthotopic implantation. In summary, our results indicate that impaired inheritance of parental histones can drive tumor progression.

Molecular Docking and Simulation-Binding Analysis of Plant Phytochemicals with the Hepatocellular Carcinoma Targets Epidermal Growth Factor Receptor and Caspase-9

Among primary liver cancers, hepatocellular carcinoma (HCC) is one of the most common forms and it has been categorized as the joint-fourth largest reason of cancer-related deaths globally. Different factors such as alcohol abuse, hepatitis B and C, viral infections, and fatty liver diseases are mainly related to the pathogenesis of HCC. In the current study, 1000 total various plant phytochemicals were docked to proteins involved in HCC. The compounds were docked to the active site amino acids of epidermal growth factor receptor and caspase-9 as receptor proteins in order to explore their inhibiting potential. The top five compounds against each receptor protein were explored as potential drug candidates on the basis of their binding affinity and root-mean square deviation values. The top two compounds against each protein were found to be liquoric acid (S-score -9.8 kcal/mol) and madecassic acid (S-score -9.3 kcal/mol) against EGFR, and limonin (S-score -10.5 kcal/mol) and obamegine (S-score -9.3 kcal/mol) against the caspase-9 protein. The selected phytochemicals were further assessed through drug scanning using Lipinski's rule of five to explore their molecular properties and druggability. According to the ADMET analysis, the selected phytochemicals were found to be non-toxic and non-carcinogenic. Finally, the molecular dynamics simulation study revealed that liquoric acid and limonin were stabilized within the binding pockets of EGFR and capase-9, respectively, and stayed firmly bound throughout the simulation. In light of the current findings, the phytochemicals reported in this study, especially liquoric acid and limonin, could be used as potential drugs for the treatment of HCC in the future.

Integrative Ligand-Based Pharmacophore Modeling, Virtual Screening, and Molecular Docking Simulation Approaches Identified Potential Lead Compounds against Pancreatic Cancer by Targeting FAK1

Pancreatic cancer is a very deadly disease with a 5-year survival rate, making it one of the leading causes of cancer-related deaths globally. Focal adhesion kinase 1 (FAK1) is a ubiquitously expressed protein in pancreatic cancer. FAK, a tyrosine kinase that is overexpressed in cancer cells, is crucial for the development of tumors into malignant phenotypes. FAK functions in response to extracellular signals by triggering transmembrane receptor signaling, which enhances focal adhesion turnover, cell adhesion, cell migration, and gene expression. The ligand-based drug design approach was used to identify potential compounds against the target protein, which included molecular docking: ADME (absorption, distribution, metabolism, and excretion), toxicity, molecular dynamics (MD) simulation, and molecular mechanics generalized born surface area (MM-GBSA). Following the retrieval of twenty hits, four compounds were selected for further evaluation based on a molecular docking approach. Three newly discovered compounds, including PubChem CID24601203, CID1893370, and CID16355541, with binding scores of -10.4, -10.1, and -9.7 kcal/mol, respectively, may serve as lead compounds for the treatment of pancreatic cancer associated with FAK1. The ADME (absorption, distribution, metabolism, and excretion) and toxicity analyses demonstrated that the compounds were effective and nontoxic. However, further wet laboratory investigations are required to evaluate the activity of the drugs against the cancer.

Distinctive expression of DNA replication factors in squamous cell carcinomas of the lip, face and oral cavity

OBJECTIVE

Uncontrolled proliferation and aberrations in cell-cycle progression are fundamental issues in cancer. In this study we aimed to determine and compare deoxyribonucleic acid (DNA) replication licensing factors at the mRNA and protein levels among squamous cell carcinomas (SCCs) of the lip, facial-skin and oral cavity.

MATERIALS AND METHODS

A total of 103 lip, oral and face SCCs were immunohistochemically stained with MCM2 (mini-chromosome maintenance 2), geminin, and ki67, and their labeling-indices were calculated. Also, 57 SCCs from the same regions along with their adjacent normal tissues underwent quantitative reverse transcription-polymerase chain reaction analysis.

RESULTS

All three proteins were overexpressed in the studied SCCs, but only geminin (P = 0.004) showed significant difference among the three regions, with higher levels in oral SCCs compared to lip (P = 0.005) and skin (P = 0.024) tumors. Geminin expression did not differ between skin- and lip-SCCs (P = 0.822). MCM2/ki67 ratio was higher in oral- compared to skin-neoplasms (P = 0.039), but no difference was found in geminin/ki67 among the SCC-subsites. There were significant differences in MCM2 and geminin mRNA between carcinomatous- and normal-tissues in all tumors, but not among the three locations.

CONCLUSION

MCM2 and geminin are involved in the tumorigenesis of lip, face and oral SCC at both mRNA- and protein-levels. Geminin may have a role in the site-specific biologic behavior of SCC. Skin SCCs had the highest proportion of licensed non-proliferating cells, while actively proliferating cells were more prominent in oral tumors. Regarding DNA replication, lip SCCs seem to be closer to skin tumors compared to their oral counterparts.

In silico investigation of the role of vitamins in cancer therapy through inhibition of MCM7 oncoprotein

An overabundance of MCM7 protein, a component of the minichromosome maintenance complex that normally initiates DNA replication, has been reported to cause different types of cancers with aggressive malignancy. Inhibition of MCM7 may lead to a significant reduction in cancer-associated cell proliferation. Despite such significance of MCM7 in cancer, the protein structure is yet to be resolved experimentally. This significantly halts the structure-guided ligand designing for cancer therapy targeting the MCM7. The present study aims to resolve the tertiary structure of MCM7 and repurpose the FDA-approved clinically used drugs for cancer therapy by targeting MCM7 protein. The secondary and 3D structures of MCM7 were generated using multiple bioinformatics tools, including the Self-Optimized Prediction Method with Alignment (SOPMA), SWISS-MODEL, and I-TASSER. The reliability of the modeled structure was assessed using PROCHECK. Initially, a structure-guided virtual screening was performed on the approved drug library to identify potential hits against MCM7. The detailed molecular mechanism of receptor interactions of the identified hits was evaluated using extensive molecular dynamics simulation. The results from this study reveal an intriguing discovery of the potential of ergocalciferol (vitamin D2), cholecalciferol (vitamin D3), ergosterol (precursor of vitamin D2) and menaquinone (vitamin K2) as oncoprotein inhibitors for cancer therapy via inhibition of MCM7.

High-Throughput Screening and Molecular Dynamics Simulation of Natural Products for the Identification of Anticancer Agents against MCM7 Protein

Minichromosome maintenance complex component 7 (MCM7) belongs to the minichromosome maintenance family that is necessary for the initiation of eukaryotic DNA replication. Overexpression of the MCM7 protein is linked to cellular proliferation and is accountable for critical malignancy in many cancers. Mechanistically, the suppression of MCM7 greatly lowers the cellular proliferation associated with cancer. Advances in immunotherapy have revolutionized treatments for many types of cancer. To date, no effective small molecular candidate has been found that can stop the advancement of cancer produced by the MCM7 protein. Here, we present the findings of methods that used a combination of structure-assisted drug design, high-throughput virtual screening, and simulations studies to swiftly generate lead compounds against MCM7 protein. In the current study, we designed efficient compounds that may combat all emerging cancer targeting the common MCM7 protein. For this objective, a molecular docking and molecular dynamics (MD) simulation-based virtual screening of 29,000 NPASS library was carried out. As a consequence of using specific pharmacological, physiological, and ADMET criteria, four new prevailing compounds, NPA000018, NPA000111, NPA00305, and NPA014826, were successfully selected. The MD simulations were also used for a time period of 50 ns to evaluate for stability and dynamics behavior of the compounds. Eventually, compounds NPA000111 and NPA014826 were found to be highly potent against MCM7 protein. According to our results, the selected compounds may be effective in treating certain cancer subtypes, for which additional follow-up experimental validation is recommended.

Integrated analysis of the clinical consequence and associated gene expression of ALK in ALK-positive human cancers

Anaplastic lymphoma kinase (ALK) is a tyrosine kinase receptor that is genetically altered in several cancers, including NSCLC, melanoma, lymphoma, and other tumors. Although ALK is associated with various cancers, the relationship between ALK expression and patient prognosis in different cancers is poorly understood. Here, using multidimensional approaches, we revealed the correlation between ALK expression and the clinical outcomes of patients with LUAD, melanoma, OV, DLBC, AML, and BC. We analyzed ALK transcriptional expression, patient survival rate, genetic alteration, protein network, and gene and microRNA (miRNA) co-expression. Compared to that in normal tissues, higher ALK expression was found in LUAD, melanoma, and OV, which are associated with poor patient survival rates. In contrast, lower transcriptional expression was found to decrease the survival rate of patients with DLBC, AML, and BC. A total of 202 missense mutations, 17 truncating mutations, 7 fusions, and 3 in-frame mutations were identified. Further, 17 genes and 19 miRNAs were found to be exclusively co-expressed and echinoderm microtubule-associated protein-like 4 (EML4) was identified as the most positively correlated gene (log odds ratio >3). The gene ontology and signaling pathways of the genes co-expressed with ALK in these six cancers were also identified. Our findings offer a basis for ALK as a prognostic biomarker and therapeutic target in cancers, which will potentially contribute to precision oncology and assist clinicians in identifying suitable treatment options.

25 (S)-Hydroxycholesterol acts as a possible dual enzymatic inhibitor of SARS-CoV-2 Mpro and RdRp-: an insight from molecular docking and dynamics simulation approaches

The coronavirus disease (COVID-19) pandemic has rapidly extended globally and killed approximately 5.83 million people all over the world. But, to date, no effective therapeutic against the disease has been developed. The disease is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and enters the host cell through the spike glycoprotein (S protein) of the virus. Subsequently, RNA-dependent RNA polymerase (RdRp) and main protease (M^pro) of the virus mediate viral transcription and replication. Mechanistically inhibition of these proteins can hinder the transcription as well as replication of the virus. Recently oxysterols and its derivative, such as 25 (S)-hydroxycholesterol (25-HC) has shown antiviral activity against SARS-CoV-2. But the exact mechanisms and their impact on RdRp and M^pro have not been explored yet. Therefore, the study aimed to identify the inhibitory activity of 25-HC against the viral enzymes RdRp and M^pro simultaneously. Initially, a molecular docking simulation was carried out to evaluate the binding activity of the compound against the two proteins. The pharmacokinetics (PK) and toxicity parameters were analyzed to observe the 'drug-likeness' properties of the compound. Additionally, molecular dynamics (MD) simulation was performed to confirm the binding stability of the compound to the targeted protein. Furthermore, molecular mechanics generalized Born surface area (MM-GBSA) was used to predict the binding free energies of the compound to the targeted protein. Molecular docking simulation identified low glide energy -51.0 kcal/mol and -35.0 kcal/mol score against the RdRp and M^pro, respectively, where MD simulation found good binding stability of the compound to the targeted proteins. In addition, the MM/GBSA approach identified a good value of binding free energies (ΔG bind) of the compound to the targeted proteins. Therefore, the study concludes that the compound 25-HC could be developed as a treatment and/or prevention option for SARS-CoV-2 disease-related complications. Although, experimental validation is suggested for further evaluation of the work.Communicated by Ramaswamy H. Sarma.

Toward the Identification of Natural Antiviral Drug Candidates against Merkel Cell Polyomavirus: Computational Drug Design Approaches

Merkel cell carcinoma (MCC) is a rare form of aggressive skin cancer mainly caused by Merkel cell polyomavirus (MCPyV). Most MCC tumors express MCPyV large T (LT) antigens and play an important role in the growth-promoting activities of oncoproteins. Truncated LT promotes tumorigenicity as well as host cell proliferation by activating the viral replication machinery, and inhibition of this protein in humans drastically lowers cellular growth linked to the corresponding cancer. Our study was designed with the aim of identifying small molecular-like natural antiviral candidates that are able to inhibit the proliferation of malignant tumors, especially those that are aggressive, by blocking the activity of viral LT protein. To identify potential compounds against the target protein, a computational drug design including molecular docking, ADME (absorption, distribution, metabolism, and excretion), toxicity, molecular dynamics (MD) simulation, and molecular mechanics generalized Born surface area (MM-GBSA) approaches were applied in this study. Initially, a total of 2190 phytochemicals isolated from 104 medicinal plants were screened using the molecular docking simulation method, resulting in the identification of the top five compounds having the highest binding energy, ranging between −6.5 and −7.6 kcal/mol. The effectiveness and safety of the selected compounds were evaluated based on ADME and toxicity features. A 250 ns MD simulation confirmed the stability of the selected compounds bind to the active site (AS) of the target protein. Additionally, MM-GBSA analysis was used to determine the high values of binding free energy (ΔG bind) of the compounds binding to the target protein. The five compounds identified by computational approaches, Paulownin (CID: 3084131), Actaealactone (CID: 11537736), Epigallocatechin 3-O-cinnamate (CID: 21629801), Cirsilineol (CID: 162464), and Lycoricidine (CID: 73065), can be used in therapy as lead compounds to combat MCPyV-related cancer. However, further wet laboratory investigations are required to evaluate the activity of the drugs against the virus.

Bioinformatics approaches identified dasatinib and bortezomib inhibit the activity of MCM7 protein as a potential treatment against human cancer

Minichromosome Maintenance Complex Component 7 (MCM7) is a key component of the DNA replication licensing factor and hexamer MCM (MCM2-7) complex that regulates the DNA replication process. The MCM7 protein is associated with tumor cell proliferation that plays an important role in different human cancer progression. As the protein is highly expressed during the cancer development process, therefore, inhibition of the protein can be utilized as a treatment option for different human cancer. However, the study aimed to identify potential small molecular drug candidates against the MCM7 protein that can utilize treatment options for human cancer. Initially, the compounds identified from protein-drugs network analysis have been retrieved from NetworkAnalyst v3.0 server and screened through molecular docking, MM-GBSA, DFT, pharmacokinetics, toxicity, and molecular dynamics (MD) simulation approach. Two compounds namely Dasatinib (CID_3062316) and Bortezomib (CID_387447) have been identified throughout the screening process, which have the highest negative binding affinity (Kcal/mol) and binding free energy (Kcal/mol). The pharmacokinetics and toxicity analysis identified drug-like properties and no toxicity properties of the compounds, where 500 ns MD simulation confirmed structural stability of the two compounds to the targeted proteins. Therefore, we can conclude that the compounds dasatinib and bortezomib can inhibit the activity of the MCM7 and can be developed as a treatment option against human cancer.

Spike protein recognizer receptor ACE2 targeted identification of potential natural antiviral drug candidates against SARS-CoV-2

Angiotensin-converting enzyme 2 (ACE2), also known as peptidyl-dipeptidase A, belongs to the dipeptidyl carboxydipeptidases family has emerged as a potential antiviral drug target against SARS-CoV-2. Most of the ACE2 inhibitors discovered till now are chemical synthesis; suffer from many limitations related to stability and adverse side effects. However, natural, and selective ACE2 inhibitors that possess strong stability and low side effects can be replaced instead of those chemicals' inhibitors. To envisage structurally diverse natural entities as an ACE2 inhibitor with better efficacy, a 3D structure-based-pharmacophore model (SBPM) has been developed and validated by 20 known selective inhibitors with their correspondence 1166 decoy compounds. The validated SBPM has excellent goodness of hit score and good predictive ability, which has been appointed as a query model for further screening of 11,295 natural compounds. The resultant 23 hits compounds with pharmacophore fit score 75.31 to 78.81 were optimized using in-silico ADMET and molecular docking analysis. Four potential natural inhibitory molecules namely D-DOPA (Amb17613565), L-Saccharopine (Amb6600091), D-Phenylalanine (Amb3940754), and L-Mimosine (Amb21855906) have been selected based on their binding affinity (−7.5, −7.1, −7.1, and −7.0 kcal/mol), respectively. Moreover, 250 ns molecular dynamics (MD) simulations confirmed the structural stability of the ligands within the protein. Additionally, MM/GBSA approach also used to support the stability of molecules to the binding site of the protein that also confirm the stability of the selected four natural compounds. The virtual screening strategy used in this study demonstrated four natural compounds that can be utilized for designing a future class of potential natural ACE2 inhibitor that will block the spike (S) protein dependent entry of SARS-CoV-2 into the host cell.

Knockout of Calcyclin Binding Protein Impedes the Growth of Breast Cancer Cells by Regulating Cell Apoptosis and β-Catenin Signaling

Breast invasive carcinoma (BRCA) is becoming the most common malignant disease worldwide, and there is intense interest in identifying diagnostic biomarkers that can be targeted for treatment of BRCA. Recent evidence has shown that calcyclin binding protein (CacyBP) can function as either a tumor promoter or suppressor during carcinogenesis. Data in The Cancer Genome Atlas (TCGA) database show that CacyBP is overexpressed in human BRCA tissues, and high levels of CacyBP are associated with shorter overall survival. Immunohistochemical staining has shown that CacyBP levels are high in cancer tissue samples and associated with a higher likelihood of disease progression. We, therefore, conducted a knockout assay to determine the role of CacyBP in the development of BRCA. Knockout of CacyBP significantly inhibited MCF7 cell proliferation and colony formation. Apoptosis was higher in CacyBP knockout cells compared with control cells. Microarray analysis showed that the CacyBP knockout caused dysregulation of numerous genes closely related to β-catenin signaling, whereas quantitative reverse-transcription PCR and immunoblotting showed that it to be inactivated. In summary, we conclude that when overexpressed, CacyBP acts as a potential oncogene for BRCA by regulating β-catenin signaling.

Transporter associated with antigen processing 1 (TAP1) expression and prognostic analysis in breast, lung, liver, and ovarian cancer

Abstract

Transporter associated with antigen processing 1 (TAP1) is a transporter protein that represent tumor antigen in the MHC I or HLA complex. Any defect in the TAP1 gene resulting in inadequate tumor tracking. TAP1 influences multidrug resistance (MDR) in human cancer cell lines and hinders the treatment during chemotherapeutic. The association of TAP1 in cancer progression remains mostly unknown and further study of the gene in relation with cancer need to conduct. Thus, the study has designed to analyze the association between the TAP1 with cancer by computationally. The expression pattern of the gene has determined by using ONCOMINE, GENT2, and GEPIA2 online platforms. The protein level of TAP1 was examined by the help of Human Protein Atlas. Samples with different clinical outcomes were investigated to evaluate the expression and promoter methylation in cancer vs. normal tissues by using UALCAN server. The copy number alteration, mutation frequency, and expression level of the gene in different cancer were analyzed by using cBioPortal server. The PrognoScan and KM plotter platforms were used to perform the survival analysis and represented graphically. Additionally, pathway and gene ontology (GO) features correlated to the TAP1 gene were analyzed and presented by bar charts. After arranging the data in a single panel like correlating expression to prognosis, mutational and alterations characteristic, and pathways analysis, we observed some interesting insights that emphasized the importance of the gene in cancer progression. The study found the relationship between the TAP1 expression pattern and prognosis in different cancer tissues and shows how TAP1 affects the clinical characteristics. The analytical data presented in the study is vital to learn about the effect of TAP1 in tumor tissue, where previously studies showing contradicting expression of TAP1 in cancer tissue. The analyzed data can also be utilized further to evade the threats against chemotherapy. Overall, the study provided a new aspect to consider the role of TAP1 gene in cancer progression and survival status.

Key messages

• This study demonstrated, for the first time, a correlation between the TAP1 gene and tumor progression.

• An upregulation of TAP1 mRNA was demonstrated in various cancer types.

• This study reported a significant negative correlation for TAP1 gene expression and the survival rate in different cancer types.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00109-021-02088-w.

Pharmacoinformatics and molecular dynamics simulation-based phytochemical screening of neem plant (Azadiractha indica) against human cancer by targeting MCM7 protein

Minichromosome maintenance complex component 7 (MCM7) belongs to the minichromosome maintenance family that is important for the initiation of eukaryotic DNA replication. Overexpression of the MCM7 protein is relative to cellular proliferation and responsible for aggressive malignancy in various cancers. Mechanistically, inhibition of MCM7 significantly reduces the cellular proliferation associated with cancer. To date, no effective small molecular candidate has been identified that can block the progression of cancer induced by the MCM7 protein. Therefore, the study has been designed to identify small molecular-like natural drug candidates against aggressive malignancy associated with various cancers by targeting MCM7 protein. To identify potential compounds against the targeted protein a comprehensive in silico drug design including molecular docking, ADME (Absorption, Distribution, Metabolism and Excretion), toxicity, and molecular dynamics (MD) simulation approaches has been applied. Seventy phytochemicals isolated from the neem tree (Azadiractha indica) were retrieved and screened against MCM7 protein by using the molecular docking simulation method, where the top four compounds have been chosen for further evaluation based on their binding affinities. Analysis of ADME and toxicity properties reveals the efficacy and safety of the selected four compounds. To validate the stability of the protein-ligand complex structure MD simulations approach has also been performed to the protein-ligand complex structure, which confirmed the stability of the selected three compounds including CAS ID:105377-74-0, CID:12308716 and CID:10505484 to the binding site of the protein. In the study, a comprehensive data screening process has performed based on the docking, ADMET properties, and MD simulation approaches, which found a good value of the selected four compounds against the targeted MCM7 protein and indicates as a promising and effective human anticancer agent.

CRISPR-cas9 Screening Identified Lethal Genes Enriched in Cell Cycle Pathway and of Prognosis Significance in Breast Cancer

Background

Lethal genes have not been systematically analyzed in breast cancer which may have significant prognostic value. The current study aims to investigate vital genes related to cell viability by analyzing the CRISPR-cas9 screening data, which may provide novel therapeutic target for patients.

Methods

Genes differentially expressed between tumor and normal tissue from the Cancer Genome Atlas (TCGA) and genes related to cell viability by CRISPR-cas9 screening from Depmap (Cancer Dependency Map) were overlapped. Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analysis was conducted to identify which pathways of overlapped genes were enriched. GSE21653 set was randomized into training and internal validation dataset at a ratio of 3:1, and external validation was performed in GSE20685 set. The least absolute shrinkage and selection operator (LASSO) regression was used to construct a signature to predict recurrence-free survival (RFS) of breast cancer patients. Univariate and multivariate Cox regression were used to evaluate the prognostic value of this signature. Differentially expressed genes (DEGs) between high-risk and low-risk patients were then analyzed to identify the main pathways regulated by this signature. Weighted correlation network analysis (WGCNA) was conducted to recognize modules correlated with high risk. Enrichment analysis was then used to identify pathways regulated by genes shared in the overlapped genes, DEGs, and WGCNA.

Results

A total of 86 oncogenes were upregulated in TCGA database and overlapped with lethal genes in Depmap database, which were enriched in cell cycle pathway. A total of 51 genes were included in the gene signature based on LASSO regression, and the median risk score of 2.36 was used as cut-off to separate low-risk patients from high-risk patients. High-risk patients showed worse RFS compared with low-risk patients in internal training, internal validation, and external validation dataset. Time-dependent receiver operating characteristic curves of 3 and 5 years indicated that risk score was superior to tumor stage, age, and PAM50 in both entire and external validation datasets. Cell cycle was the main different pathway between the high-risk and low-risk groups. Meanwhile, cell cycle was also the main pathway enriched in the 25 genes which were shared among 86 genes, DEGs, and WGCNA.

Conclusion

Cell cycle pathway, identified by CRISPR-cas9 screening, was a key pathway regulating cell viability, which has significant prognostic values and can serve as a new target for breast cancer patient treatment.

Structure based pharmacophore modeling, virtual screening, molecular docking and ADMET approaches for identification of natural anti-cancer agents targeting XIAP protein

X-linked inhibitor of apoptosis protein (XIAP) is a member of inhibitor of apoptosis protein (IAP) family responsible for neutralizing the caspases-3, caspases-7, and caspases-9. Overexpression of the protein decreased the apoptosis process in the cell and resulting development of cancer. Different types of XIAP antagonists are generally used to repair the defective apoptosis process that can eliminate carcinoma from living bodies. The chemically synthesis compounds discovered till now as XIAP inhibitors exhibiting side effects, which is making difficulties during the treatment of chemotherapy. So, the study has design to identifying new natural compounds that are able to induce apoptosis by freeing up caspases and will be low toxic. To identify natural compound, a structure-based pharmacophore model to the protein active site cavity was generated following by virtual screening, molecular docking and molecular dynamics (MD) simulation. Initially, seven hit compounds were retrieved and based on molecular docking approach four compounds has chosen for further evaluation. To confirm stability of the selected drug candidate to the target protein the MD simulation approach were employed, which confirmed stability of the three compounds. Based on the finding, three newly obtained compounds namely Caucasicoside A (ZINC77257307), Polygalaxanthone III (ZINC247950187), and MCULE-9896837409 (ZINC107434573) may serve as lead compounds to fight against the treatment of XIAP related cancer, although further evaluation through wet lab is necessary to measure the efficacy of the compounds.

Validation of CSN1S1 transcriptional expression, promoter methylation, and prognostic power in breast cancer using independent datasets

Breast cancer ranked second among most frequent cancer in the world playing a significant role in mortality rate. Having prior knowledge on differentially expressed genes in breast cell carcinoma elucidated important indications to understand the molecular mechanism underneath breast carcinogenesis. In this study we have investigated the distinguished CSN1S1 expression in human breast cancer. We have analyzed CSN1S1 mRNA expression between cancer and normal tissues using TCGA datasets. Moreover, analysis including promoter methylation, mutations, prognosis, co-expression, gene ontology, and pathways of CSN1S1 were performed by the TCGA Wanderer, UCSC Xena, cBioPortal, PrognoScan, UALCAN, and Enricher server. We have observed low mRNA expression and high promoter methylation of CSN1S1 in cancer tissues compared to normal tissues. Furthermore, we have also identified low mRNA expression in clinicopathological patients, as well as 9 deleterious mutations with highly co-expressed protein MRC1, and significantly related signaling pathways. We have found a positive correlation between the lower expression of CSN1S1 and patients surviving with breast cancer. Here we have concluded that CSN1S1 acts as a biomarker for the surveillance and prognosis of breast cancer, and also works as a novel therapeutic target at the molecular and pathway levels.

•

Low transcriptional expression and low survival rate of CSN1S1 in breast cancer.

•

The investigation of clinical profiles and mutational positions of CSN1S1 in breast cancer.

•

The investigation of gene ontology and signaling pathway of CSN1S1 and their co-expressed genes.

•

We identified CSN1S1 and also their co-expressed proteins are the potential biomarkers in breast cancer.

High expression of bone morphogenetic protein 1 (BMP1) is associated with a poor survival rate in human gastric cancer, a dataset approaches

Bone morphogenetic protein 1 (BMP1) is a secreted metalloprotease of the astacin M12A family of bone morphogenetic proteins (BMPs). BMP1 activates transforming growth factor-β (TGF-β) and BMP signaling pathways by proteolytic cleavage, which has dual roles in gastrointestinal tumor development and progression.TGF-β promotes invasion and metastasis of gastric cancer (GC) by the help of BMP1, so upregulation of the BMP1 may increase cancer invasiveness in GC. In this study,the transcriptional expression, mutations, survival rate, TFs, miRNAs, gene ontology, and signaling pathways of BMP1 were analyzed by using different web servers. We found higher transcriptional and clinicopathological characteristics expression compared to normal tissues, worsening survival rate in GC. We detected 25 missenses, 15 truncating mutations, 23 TFs, and 8 miRNAs. Finally, we identified and analyzed the co-expressed genes and found that the leukemia inhibitory factor is the most positively correlated gene. The gene ontological features and signaling pathways involved in GC development were evaluated as well. We believe that this study will provide a basis for BMP1 to be a significant biomarker for human GC prognosis.