Novel DNA Bis-Intercalator XR5944 as a Potent Anticancer Drug-Design and Mechanism of Action

Interactions between DNA and the acridine intercalator: A computational study

Cancer is a global public health problem characterized by deviations in the mechanisms that control cell proliferation, resulting in mutations and variations in the structure of DNA. The mechanisms of action of chemotherapeutic drugs are related to their interactions and binding with DNA; consequently, the development of antineoplastic agents that target DNA has extensively focused on use of acridine, a heterocyclic molecule that binds to deoxyribonucleic acid via intercalation, a process that modifies DNA and makes replication impossible. In this context, this study aimed to computationally investigate how acridine intercalators interact with DNA by evaluating the mechanism of interactions, binding, and interaction energies using quantum mechanics calculations. Molecular electrostatic potential (MEP) analysis revealed that acridine has well- distributed negative charges in the center of the molecule, indicative of a dominant electron-rich region. Acridine exhibits well-defined π orbitals (HOMO and LUMO) on the aromatic rings, suggesting that charge transfer occurs within the molecule and may be responsible for the pharmacological activity of the compound. Structural analysis revealed that acridine interacts with DNA mainly through hydrogen bonds between HAcridine… ODNA with bond lengths ranging from 2.370 Å to 3.472 Å. The Binding energy (ΔEBind) showed that acridine interacts with DNA effectively for all complexes and the electronic energy results (E+ZPE) for complexes revealed that the complexes are more stable when the DNA-centered acridine molecule. The Laplacian-analysis topological QTAIM parameter (∇2ρ(r)) and total energy (H(r)) categorized the interactions as being non-covalent in nature. The RGD peak distribution in the NCI analysis reveals the presence of van der Waals interactions, predominantly between the intercalator and DNA. Accordingly, we confirm that acridine/DNA interactions are relevant for understanding how the intercalator acts within nucleic acids.

A DNA Replication Stress-Based Prognostic Model for Lung Adenocarcinoma

Tumor cells endure continuous DNA replication stress, which opens the way to cancer development. Despite previous research, the prognostic implications of DNA replication stress on lung adenocarcinoma (LUAD) have yet to be investigated. Here, we aimed to investigate the potential of DNA replication stress-related genes (DNARSs) in predicting the prognosis of individuals with LUAD. Differentially expressed genes (DEGs) originated from the TCGA-LUAD dataset, and we constructed a 10-gene LUAD prognostic model based on DNARSs-related DEGs (DRSDs) using Cox regression analysis. The receiver operating characteristic (ROC) curve demonstrated excellent predictive capability for the LUAD prognostic model, while the Kaplan-Meier survival curve indicated a poorer prognosis in a high-risk (HR) group. Combined with clinical data, the Riskscore was found to be an independent predictor of LUAD prognosis. By incorporating Riskscore and clinical data, we developed a nomogram that demonstrated a capacity to predict overall survival and exhibited clinical utility, which was validated through the calibration curve, ROC curve, and decision curve analysis curve tests, confirming its effectiveness in prognostic evaluation. Immune analysis revealed that individuals belonging to the low-risk (LR) group exhibited a greater abundance of immune cell infiltration and higher levels of immune function. We calculated the immunopheno score and TIDE scores and tested them on the IMvigor210 and GSE78220 cohorts and found that individuals categorized in the LR group exhibited a higher likelihood of deriving therapeutic benefits from immunotherapy intervention. Additionally, we predicted that patients classified in the HR group would demonstrate enhanced sensitivity to Docetaxel using anti-tumor drugs. To summarize, we successfully developed and validated a prognostic model for LUAD by incorporating DNA replication stress as a key factor.

In silico study of HASDI (high-affinity selective DNA intercalator) as a new agent capable of highly selective recognition of the DNA sequence

Cancer as an acquired genetic disease is based on changes both in the genome itself and in transcription processes. Accordingly, it is at the DNA level that it makes sense to search for and design agents capable of effective and selective anticancer action. In this study, we used an iterative approach based on a molecular dynamics simulation to design a highly selective DNA-intercalating agent called HASDI. To confirm its selective affinity to DNA, we conducted two simulation experiments: HASDI in a complex with a DNA fragment of the EBNA1 gene (it targets 16 nucleotide pairs of this gene) and HASDI in a complex with a random DNA fragment of the KCNH2 gene. The molecular dynamics simulation was carried out in the GROMACS 2019 package. The binding energy was calculated by gmx_MMPBSA 1.5.2. The further analysis was performed using the built-in utilities of GROMACS, gmx_MMPBSA and also XMGRACE and Pymol 1.8. As a result, we determined that the EBNA1-50nt/HASDI complex was stable throughout the whole simulation trajectory. HASDI, due to the presence of a linker modified depending on a specific pair of nitrogenous bases, formed an average of 32 hydrogen bonds with a sequence of 16 nucleotide pairs. Phenazine rings were stably intercalated every 2 base pairs. The root-mean-square deviation of HASDI in such a complex fluctuated around the value of 6.5 Å and had no tendency to increase. The calculated value of the binding free energy was - 235.3 ± 7.77 kcal/mol. The KCNH2-50nt/HASDI complex, as an example of the intercalation of the designed structure into a random part of the human genome, maintained the stability of its position at a level comparable to the EBNA1-50nt/HASDI complex. The phenazine rings were constantly intercalated in their original positions, and the root-mean-square deviation fluctuated around one value, although it had a tendency to chaotic changes. At the same time, this complex was characterized by 17-19 hydrogen bonds, on average, and the binding free energy was - 193.47 ± 14.09 kcal/mol. Moreover, the DNA duplex had local single-nucleotide melting in the region of the 4th linker. According to a significant decrease in the number of hydrogen bonds, a decrease in energy gain, as well as a decrease in the stability of the DNA duplex characteristic of the KCNH2-50nt/HASDI complex compared to the target EBNA1-50nt/HASDI complex, the molecule we designed can be considered a potentially selective DNA polyintercalating agent capable of relatively accurate recognition of 16 base pairs.

Biological Activity of Newly Synthesized Benzimidazole and Benzothizole 2,5-Disubstituted Furane Derivatives

Newly designed and synthesized cyano, amidino and acrylonitrile 2,5-disubstituted furane derivatives with either benzimidazole/benzothiazole nuclei have been evaluated for antitumor and antimicrobial activity. For potential antitumor activity, the compounds were tested in 2D and 3D cell culture methods on three human lung cancer cell lines, A549, HCC827 and NCI-H358, with MTS cytotoxicity and BrdU proliferation assays in vitro. Compounds 5, 6, 8, 9 and 15 have been proven to be compounds with potential antitumor activity with high potential to stop the proliferation of cells. In general, benzothiazole derivatives were more active in comparison to benzimidazole derivatives. Antimicrobial activity was evaluated with Broth microdilution testing (according to CLSI (Clinical Laboratory Standards Institute) guidelines) on Gram-negative Escherichia coli and Gram-positive Staphylococcus aureus. Additionally, Saccharomyces cerevisiae was included in testing as a eukaryotic model organism. Compounds 5, 6, 8, 9 and 15 showed the most promising antibacterial activity. In general, the compounds showed antitumor activity, higher in 2D assays in comparison with 3D assays, on all three cell lines in both assays. In natural conditions, compounds with such an activity profile (less toxic but still effective against tumor growth) could be promising new antitumor drugs. Some of the tested compounds showed antimicrobial activity. In contrast to ctDNA, the presence of nitro group or chlorine in selected furane-benzothiazole structures did not influence the binding mode with AT-DNA. All compounds dominantly bound inside the minor groove of AT-DNA either in form of monomers or dimer and higher-order aggregates.