Molecular dynamics simulation studies for DNA sequence recognition by reactive metabolites of anticancer compounds

Unraveling Bacterial Single-Stranded Sequence Specificities: Insights from Molecular Dynamics and MMPBSA Analysis of Oligonucleotide Probes

We utilized molecular dynamics (MD) simulations and Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) free energy calculations to investigate the specificity of two oligonucleotide probes, namely probe B and probe D, in detecting single-stranded DNA (ssDNA) within three bacteria families: Enterobacteriaceae, Pasteurellaceae, and Vibrionaceae. Due to the limited understanding of molecular mechanisms in the previous research, we have extended the discussion to focus specifically on investigating the binding process of bacteria-probe DNA duplexes, with an emphasis on analyzing the binding free energy. The role of electrostatic contributions in the specificity between the oligonucleotide probes and the bacterial ssDNAs was investigated and found to be crucial. Our calculations yielded results that were highly consistent with the experimental data. Through our study, we have successfully exhibited the benefits of utilizing in-silico approaches as a powerful virtual-screening tool, particularly in research areas that demand a thorough comprehension of molecular interactions.

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.

Potential of cytochrome P450, a family of xenobiotic metabolizing enzymes, in cancer therapy

SIGNIFICANCE

Targeted cancer therapy with minimal off-target consequences has shown promise for some cancer types. Although cytochrome P450 (CYP) consists of 18 families, CYP1-4 families play key role in metabolizing xenobiotics and cancer drugs. This eventually affects the process of carcinogenesis, treatment outcome, and cancer drug resistance. Differential overexpression of CYPs in transformed cells, together with phenotypic alterations in tumors, presents a potential for therapeutic intervention.

RECENT ADVANCES

Recent advances in molecular tools and information technology have helped utilize CYPs as cancer targets. The precise expression in various tumors, X-ray crystal structures, improved understanding of the structure-activity relationship, and new approaches in the development of prodrugs have supported the ongoing efforts to develop CYPs-based drugs with a better therapeutic index.

CRITICAL ISSUES

Narrow therapeutic index, off-target effects, drug resistance, and tumor heterogeneity limit the benefits of CYP-based conventional cancer therapies. In this review, we address the CYP1-4 families as druggable targets in cancer. An emphasis is given to the CYP expression, function, and the possible mechanisms that drive expression and activity in normal and transformed tissues. The strategies that inhibit or activate CYPs for therapeutic benefits are also discussed.

FUTURE DIRECTIONS

Efforts are needed to develop more selective tools that will help comprehend molecular and metabolic alterations in tumor tissues with biological end-points in relation to CYPs. This will eventually translate to developing more specific CYP inhibitors/inducers.

Molecular Modeling Study of the Genotoxicity of the Sudan I and Sudan II Azo Dyes and Their Metabolites

Azo dyes are defined by the presence of a characteristic N=N group. Sudan I and Sudan II are synthetic azo dyes that have been used as coloring agents. Although animal toxicity studies suggest that Sudan dyes are mutagenic, their molecular mechanism of action is unknown, thus making it challenging to establish thresholds for tolerable daily intake or to understand how these molecules could be modified to ameliorate toxicity. In addition, dye metabolites, such as azobiphenyl and 4-aminobiphenyl, have been correlated with epigenetic alterations. We shed some light on the mechanisms of Sudan dye genotoxicity through a molecular modeling study of Sudan I and Sudan II dyes and two common metabolites interacting with DNA as adducts. The results suggest that all four adducts cause significant perturbations to the DNA helical conformation and structure; thus, it can be inferred that DNA repair and replication processes would be significantly impacted.

The aryl hydrocarbon receptor (AhR) as a breast cancer drug target

Breast cancer is the most common cancer in women, with more than 1.7 million diagnoses worldwide per annum. Metastatic breast cancer remains incurable, and the presence of triple-negative phenotypes makes targeted treatment impossible. The aryl hydrocarbon receptor (AhR), most commonly associated with the metabolism of xenobiotic ligands, has emerged as a promising biological target for the treatment of this deadly disease. Ligands for the AhR can be classed as exogenous or endogenous and may have agonistic or antagonistic activity. It has been well reported that agonistic ligands may have potent and selective growth inhibition activity in a number of oncogenic cell lines, and one (aminoflavone) has progressed to phase I clinical trials for breast cancer sufferers. In this study, we examine the current state of the literature in this area and elucidate the promising advances that are being made in hijacking the cytosolic-to-nuclear pathway of the AhR for the possible future treatment of breast cancer.

Molecular mechanism and binding free energy of doxorubicin intercalation in DNA

The intercalation process of binding doxorubicin (DOX) in DNA is studied by extensive MD simulations. Many molecular factors that control the binding affinity of DOX to DNA to form a stable complex are inspected and quantified by employing continuum solvation models for estimating the binding free energy. The modified MM-PB(GB)SA methodology provides a complete energetic profile of ΔGele, ΔGvDW, ΔGpolar, ΔGnon-polar, TΔStotal, ΔGdeform, ΔGcon, and ΔGion. To identify the sequence specificity of DOX, two different DNA sequences, d(CGATCG) or DNA1 and d(CGTACG) or DNA2, with one molecule (1 : 1 complex) or two molecule (2 : 1 complex) configurations of DOX were selected in this study. Our results show that the DNA deformation energy (ΔGdeform), the energy cost from translational and rotational entropic contributions (TΔStran+rot), the total electrostatic interactions (ΔGpolar-PB/GB + ΔGele) of incorporation, the intramolecular electrostatic interactions (ΔGele) and electrostatic polar solvation interactions (ΔGpolar-PB/GB) are all unfavorable to the binding of DOX to DNA. However, they are overcome by at least five favorable interactions: the van der Waals interactions (ΔGvDW), the non-polar solvation interaction (ΔGnon-polar), the vibrational entropic contribution (TΔSvib), and the standard concentration dependent free energies of DOX (ΔGcon) and the ionic solution (ΔGion). Specifically, the van der Waals interaction appears to be the major driving force to form a stable DOX-DNA complex. We also predict that DOX has stronger binding to DNA1 than DNA2. The DNA deformation penalty and entropy cost in the 2 : 1 complex are less than those in the 1 : 1 complex, thus they indicate that the 2 : 1 complex is more stable than the 1 : 1 complex. We have calculated the total binding free energy (BFE) (ΔGt-sim) using both MM-PBSA and MM-GBSA methods, which suggests a more stable DOX-DNA complex at lower ionic concentration. The calculated BFE from the modified MM-GBSA method for DOX-DNA1 and DOX-DNA2 in the 1 : 1 complex is -9.1 and -5.1 kcal mol-1 respectively. The same quantities from the modified MM-PBSA method are -12.74 and -8.35 kcal mol-1 respectively. The value of the total BFE ΔGt-sim in the 1 : 1 complex is in reasonable agreement with the experimental value of -7.7 ± 0.3 kcal mol-1.

Recent Developments and Applications of the MMPBSA Method

The Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approach has been widely applied as an efficient and reliable free energy simulation method to model molecular recognition, such as for protein-ligand binding interactions. In this review, we focus on recent developments and applications of the MMPBSA method. The methodology review covers solvation terms, the entropy term, extensions to membrane proteins and high-speed screening, and new automation toolkits. Recent applications in various important biomedical and chemical fields are also reviewed. We conclude with a few future directions aimed at making MMPBSA a more robust and efficient method.

Identification of leads for antiproliferative activity on MDA-MB-435 human breast cancer cells through pharmacophore and CYP1A1-mediated metabolism

CYP1A1 is a potential target for anticancer drug development due to its overexpression in certain cancer cells and role in cancer progression. To identify new leads for CYP1A1 mediated anticancer action, we attempted ligand based pharmacophore mapping, virtual screening of databases, molecular docking, MetaSite based filtering, and molecular dynamics simulations. Initial computational and in vitro screening identified 11 compounds from which we identified two lead compounds, ZINC33468944 and ZINC32101539, showed potential antitumor activity on MDA-MB-435 cell lines (GI50 < 0.1 μM) and CYP1A1 inhibition of 0.13 and 0.3 μM, respectively. Furthermore, the lead compounds were evaluated for CYP1A1 mediated metabolism, showing N-hydroxylated metabolites, which have potential of DNA adduct formation and cause cancerous cell death. Analysis of molecular dynamics simulations provided important guidelines for the further modification of the lead compounds. Hence, we claim the lead molecules for further development in anticancer drug discovery.

Characterization of differences in substrate specificity among CYP1A1, CYP1A2 and CYP1B1: an integrated approach employing molecular docking and molecular dynamics simulations

Recent trends in new drug discovery of anticancer drugs have made oncologists more aware of the fact that the new drug discovery must target the developing mechanism of tumorigenesis to improve the therapeutic efficacy of antineoplastic drugs. The drugs designed are expected to have high affinity towards the novel targets selectively. Current research highlights overexpression of CYP450s, particularly cytochrome P450 1A1 (CYP1A1), in tumour cells, representing a novel target for anticancer therapy. However, the CYP1 family is identified as posing significant problems in selectivity of anticancer molecules towards CYP1A1. Three members have been identified in the human CYP1 family: CYP1A1, CYP1A2 and CYP1B1. Although sequences of the three isoform have high sequence identity, they have distinct substrate specificities. The understanding of macromolecular features that govern substrate specificity is required to understand the interplay between the protein function and dynamics, design novel antitumour compounds that could be specifically metabolized by only CYP1A1 to mediate their antitumour activity and elucidate the reasons for differences in substrate specificity profile among the three proteins. In the present study, we employed a combination of computational methodologies: molecular docking and molecular dynamics simulations. We utilized eight substrates for elucidating the difference in substrate specificity of the three isoforms. Lastly, we conclude that the substrate specificity of a particular substrate depends upon the type of the active site residues, the dynamic motions in the protein structure upon ligand binding and the physico-chemical characteristics of a particular ligand. Copyright © 2016 John Wiley & Sons, Ltd.