In Silico Identification of Piperidinyl-amine Derivatives as Novel Dual Binders of Oncogene c-myc/c-Kit G-quadruplexes

Comprehensive analysis of intramolecular G-quadruplex structures: furthering the understanding of their formalism

G-quadruplexes (G4) are helical structures found in guanine-rich DNA or RNA sequences. Generally, their formalism is based on a few dozen structures, which can produce some inconsistencies or incompleteness. Using the website ASC-G4, we analyzed the structures of 333 intramolecular G4s, of all types, which allowed us to clarify some key concepts and present new information. To each of the eight distinguishable topologies corresponds a groove-width signature and a predominant glycosidic configuration (gc) pattern governed by the directions of the strands. The relative orientations of the stacking guanines within the strands, which we quantified and related to their vertical gc successions, determine the twist and tilt of the helices. The latter impact the minimum groove widths, which represent the space available for lateral ligand binding. The G4 four helices have similar twists, even when these twists are irregular, meaning that they have various angles along the strands. Despite its importance, the vertical gc succession has no strict one-to-one relationship with the topology, which explains the discrepancy between some topologies and their corresponding circular dichroism spectra. This study allowed us to introduce the new concept of platypus G4s, which are structures with properties corresponding to several topologies.

Target-Based Small Molecule Drug Discovery for Colorectal Cancer: A Review of Molecular Pathways and In Silico Studies

Colorectal cancer is one of the most prevalent cancer types. Although there have been breakthroughs in its treatments, a better understanding of the molecular mechanisms and genetic involvement in colorectal cancer will have a substantial role in producing novel and targeted treatments with better safety profiles. In this review, the main molecular pathways and driver genes that are responsible for initiating and propagating the cascade of signaling molecules reaching carcinoma and the aggressive metastatic stages of colorectal cancer were presented. Protein kinases involved in colorectal cancer, as much as other cancers, have seen much focus and committed efforts due to their crucial role in subsidizing, inhibiting, or changing the disease course. Moreover, notable improvements in colorectal cancer treatments with in silico studies and the enhanced selectivity on specific macromolecular targets were discussed. Besides, the selective multi-target agents have been made easier by employing in silico methods in molecular de novo synthesis or target identification and drug repurposing.

Identification of Trovafloxacin, Ozanimod, and Ozenoxacin as potent c-Myc G-quadruplex stabilizers to suppress c-Myc transcription and myeloma growth

c-Myc is a major oncogene that is estimated to result in almost all human cancers and the c-Myc downregulation has become an attractive strategy for cancer treatment. For it is hard to design compounds that can directly interact with the c-Myc protein, the DNA G-quadruplex (G4) was discovered in its promoter region which was referred to as a potential drug target for controlling c-Myc expression. In this study, a combined strategy of molecular docking-based virtual screening, molecular dynamics (MD) simulation, and molecular mechanics/generalized Born surface area (MM/GBSA) free energy calculation was conducted on the existing FDA-Approved Drugs Library, eight compounds were selected for further experimental assay. Among them, five compounds exhibited dose-dependently anticancer activities against RPMI-8226 cells with IC50 values less than 18.4 μM. Further experiments showed that Trovafloxacin, Ozanimod, and Ozenoxacin decreased c-Myc mRNA level obviously and downregulated c-Myc expression significantly. In summary, compounds Trovafloxacin, Ozanimod, and Ozenoxacin might be regarded as new c-Myc G4 stabilizers for the treatment of c-Myc related cancers in the future.

Polymorphic and Higher-Order G-Quadruplexes as Possible Transcription Regulators: Novel Perspectives for Future Anticancer Therapeutic Applications

In the past two decades, significant efforts have been put into designing small molecules to target selected genomic sites where DNA conformational rearrangements control gene expression. G-rich sequences at oncogene promoters are considered good points of intervention since, under specific environmental conditions, they can fold into non-canonical tetrahelical structures known as G-quadruplexes. However, emerging evidence points to a frequent lack of correlation between small molecule targeting of G-quadruplexes at gene promoters and the expression of the associated protein, which hampers pharmaceutical applications. The wide genomic localization of G-quadruplexes along with their highly polymorphic behavior may account for this scenario, suggesting the need for more focused drug design strategies. Here, we will summarize the G4 structural features that can be considered to fulfill this goal. In particular, by comparing a telomeric sequence with the well-characterized G-rich domain of the KIT promoter, we will address how multiple secondary structures might cooperate to control genome architecture at a higher level. If this holds true, the link between drug–DNA complex formation and the associated cellular effects will need to be revisited.

Evaluating Molecular Docking Software for Small Molecule Binding to G-Quadruplex DNA

G-quadruplexes are four-stranded nucleic acid secondary structures of biological significance and have emerged as an attractive drug target. The G4 formed in the MYC promoter (MycG4) is one of the most studied small-molecule targets, and a model system for parallel structures that are prevalent in promoter DNA G4s and RNA G4s. Molecular docking has become an essential tool in structure-based drug discovery for protein targets, and is also increasingly applied to G4 DNA. However, DNA, and in particular G4, binding sites differ significantly from protein targets. Here we perform the first systematic evaluation of four commonly used docking programs (AutoDock Vina, DOCK 6, Glide, and RxDock) for G4 DNA-ligand binding pose prediction using four small molecules whose complex structures with the MycG4 have been experimentally determined in solution. The results indicate that there are considerable differences in the performance of the docking programs and that DOCK 6 with GB/SA rescoring performs better than the other programs. We found that docking accuracy is mainly limited by the scoring functions. The study shows that current docking programs should be used with caution to predict G4 DNA-small molecule binding modes.

In Silico Food-Drug Interaction: A Case Study of Eluxadoline and Fatty Meal

Food-drug interaction is an infrequently considered aspect in clinical practice. Usually, drugs are taken together with meals and what follows may adversely affect pharmacokinetic and pharmacodynamic properties, and hence, the therapeutic effects. In this study, a computational protocol was proposed to explain the different assimilations of two µ-receptors agonists, eluxadoline and loperamide, with a peculiar pharmacokinetic profile. Compared to loperamide, eluxadoline is absorbed less after the intake of a fatty meal, and the LogP values do not explain this event. Firstly, keeping in mind the different pH in the intestinal tract, the protonation states of both compounds were calculated. Then, all structures were subjected to a conformational search by using MonteCarlo and Molecular Dynamics methods, with solvation terms mimicking the water and weak polar solvent (octanol). Both computational results showed that eluxadoline has less conformational freedom in octanol, unlike loperamide, which exhibits constant behavior in both solvents. Therefore, we hypothesize that fatty meal causes the "closure" of the eluxadoline molecule to prevent the exposure of the polar groups and their interaction with water, necessary for the drug absorption. Based on our results, this work could be a reasonable "case study", useful for future investigation of the drug pharmacokinetic profile.

Natural Products Extracted from Fungal Species as New Potential Anti-Cancer Drugs: A Structure-Based Drug Repurposing Approach Targeting HDAC7

Mushrooms can be considered a valuable source of natural bioactive compounds with potential polypharmacological effects due to their proven antimicrobial, antiviral, antitumor, and antioxidant activities. In order to identify new potential anticancer compounds, an in-house chemical database of molecules extracted from both edible and non-edible fungal species was employed in a virtual screening against the isoform 7 of the Histone deacetylase (HDAC). This target is known to be implicated in different cancer processes, and in particular in both breast and ovarian tumors. In this work, we proposed the ibotenic acid as lead compound for the development of novel HDAC7 inhibitors, due to its antiproliferative activity in human breast cancer cells (MCF-7). These promising results represent the starting point for the discovery and the optimization of new HDAC7 inhibitors and highlight the interesting opportunity to apply the "drug repositioning" paradigm also to natural compounds deriving from mushrooms.

In Silico Identification and Biological Evaluation of Antioxidant Food Components Endowed with IX and XII hCA Inhibition

The tumor-associated isoenzymes hCA IX and hCA XII catalyze the hydration of carbon dioxide to bicarbonate and protons. These isoforms are highly overexpressed in many types of cancer, where they contribute to the acidification of the tumor environment, promoting tumor cell invasion and metastasis. In this work, in order to identify novel dual hCA IX and XII inhibitors, virtual screening techniques and biological assays were combined. A structure-based virtual screening towards hCA IX and XII was performed using a database of approximately 26,000 natural compounds. The best shared hits were submitted to a thermodynamic analysis and three promising best hits were identified and evaluated in terms of their hCA IX and XII inhibitor activity. In vitro biological assays were in line with the theoretical studies and revealed that syringin, lithospermic acid, and (-)-dehydrodiconiferyl alcohol behave as good hCA IX and hCA XII dual inhibitors.

Inside Perspective of the Synthetic and Computational Toolbox of JAK Inhibitors: Recent Updates

The mechanisms of inflammation and cancer are intertwined by complex networks of signaling pathways. Dysregulations in the Janus kinase/signal transducer and activator of transcription (JAK/STAT) pathway underlie several pathogenic conditions related to chronic inflammatory states, autoimmune diseases and cancer. Historically, the potential application of JAK inhibition has been thoroughly explored, thus triggering an escalation of favorable results in this field. So far, five JAK inhibitors have been approved by the Food and Drug Administration (FDA) for the treatment of different diseases. Considering the complexity of JAK-depending processes and their involvement in multiple disorders, JAK inhibitors are the perfect candidates for drug repurposing and for the assessment of multitarget strategies. Herein we reviewed the recent progress concerning JAK inhibition, including the innovations provided by the release of JAKs crystal structures and the improvement of synthetic strategies aimed to simplify of the industrial scale-up.

Multi-Targeting Bioactive Compounds Extracted from Essential Oils as Kinase Inhibitors

Essential oils (EOs) are popular in aromatherapy, a branch of alternative medicine that claims their curative effects. Moreover, several studies reported EOs as potential anti-cancer agents by inducing apoptosis in different cancer cell models. In this study, we have considered EOs as a potential resource of new kinase inhibitors with a polypharmacological profile. On the other hand, computational methods offer the possibility to predict the theoretical activity profile of ligands, discovering dangerous off-targets and/or synergistic effects due to the potential multi-target action. With this aim, we performed a Structure-Based Virtual Screening (SBVS) against X-ray models of several protein kinases selected from the Protein Data Bank (PDB) by using a chemoinformatics database of EOs. By evaluating theoretical binding affinity, 13 molecules were detected among EOs as new potential kinase inhibitors with a multi-target profile. The two compounds with higher percentages in the EOs were studied more in depth by means Induced Fit Docking (IFD) protocol, in order to better predict their binding modes taking into account also structural changes in the receptor. Finally, given its good binding affinity towards five different kinases, cinnamyl cinnamate was biologically tested on different cell lines with the aim to verify the antiproliferative activity. Thus, this work represents a starting point for the optimization of the most promising EOs structure as kinase inhibitors with multi-target features.

Understanding The Role of Reline, a Natural DES, on Temperature-Induced Conformational Changes of C-Kit G-Quadruplex DNA: A Molecular Dynamics Study

The noncanonical guanine-rich DNAs have drawn particular attention to the scientific world due to their controllable diverse and polymorphic structures. Apart from biological and medical significance, G-quadruplex DNAs are widely used in various fields such as nanotechnology, nanomachine, biosensors, and biocatalyst. So far, the applications of the G-quadruplex DNA are mainly limited in the water medium. Recently, a new generation of solvent named deep eutectic solvent (DES) has become very popular and has been widely used as a reaction medium of biocatalytic reactions and long-term storage medium for nucleic acids, even at high temperature. Hence, it is essential to understand the role of DES on temperature-induced conformational changes of a G-quadruplex DNA. In this research work, we have explored the temperature-mediated conformational dynamics of c-kit oncogene promoter G-quadruplex DNA in reline medium in the temperature range of 300-500 K, using a total of 10 μs unbiased all-atom molecular dynamics simulation. Here, from RMSD, RMSF, R_g and principal component analyses, we notice that the c-kit G-quadruplex DNA is stable up to 450 K in reline medium. However, it unfolds in water medium at 450 K. It is found that the hydrogen bonding interactions between c-kit G-quadruplex DNA and reline play a key role in the stabilization of the G-quadruplex DNA even at high temperature. Furthermore, in this work we have observed a very interesting and distinctive phenomenon of the central cation of the G-quadruplex DNA. Its position was seen to fluctuate between the two tetrad cores, that is, the region between tetrad-1 and tetrad-2 and that between tetrad-2 and tetrad-3 and vice versa at 450 and 500 K in reline medium which is absent in water medium at 450 K. Moreover, the rate of its oscillation is increased when temperature is increased.

Folding intermediate states of the parallel human telomeric G-quadruplex DNA explored using Well-Tempered Metadynamics

An increasingly comprehension of the folding intermediate states of DNA G-quadruplexes (G4s) is currently an important scientific challenge, especially for the human telomeric (h-tel) G4s-forming sequences, characterized by a highly polymorphic nature. Despite the G-triplex conformation was proposed as one of the possible folding intermediates for the antiparallel and hybrid h-tel G4s, for the parallel h-tel topology with an all-anti guanine orientation, a vertical strand-slippage involving the G-triplets was proposed in previous works through microseconds-long standard molecular dynamics simulations (MDs). Here, in order to get further insights into the vertical strand-slippage and the folding intermediate states of the parallel h-tel G4s, we have carried out a Well-Tempered Metadynamics simulation (WT-MetaD), which allowed us to retrieve an ensemble of six G4s having two/G-tetrad conformations derived by the G-triplets vertical slippage. The insights highlighted in this work are aimed at rationalizing the mechanistic characterisation of the parallel h-tel G4 folding process.

Method for Systematic Analogue Search Using the Mega SAR Matrix Database

Analogue searching is a typical requirement in hit expansion, hit-to-lead, and lead optimization projects. A new computational methodology is introduced to search for existing and virtual analogues of active compounds. The approach is based upon the SAR matrix (SARM) data structure that was originally developed for the systematic identification and structural organization of analogue series. The SARM-based analogue search algorithm further extends the capacity of current substructure-based methods by (i) simultaneously considering existing and virtual analogues that populate chemical space around query compounds, (ii) permitting not only R-group replacements but also well-defined chemical modifications in core structures to further expand the analogue space, and (iii) automatically extracting all possible analogues from large pools. In addition, as a basis for analogue searching following the SARM concept, the Mega-SARM database is introduced. Mega-SARM is derived from nearly 3.7 million compounds and contains ∼250 000 matrices with structurally related analogue series and more than 1.5 million virtual candidate compounds.

Binding of quinazolinones to c-KIT G-quadruplex; an interplay between hydrogen bonding and π-π stacking

Stabilization of G-quadruplex structures in the c-KIT promoter with the aid of ligands has become an area of great interest in potential cancer therapeutics. Understanding the binding process between ligands and G-quadruplex is essential for a discovery of selective ligands with high binding affinity to G-quadruplex. In the present work, binding mechanisms of 4-quinazolinones to c-KIT G-quadruplex were investigated theoretically by means of molecular dynamics (MD) simulations. To explore the binding affinity of ligands, binding free energy calculations were performed using the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method. We demonstrate that the key interactions in G-quadruplex-ligand complexes are π-π stacking and hydrogen bond interactions. However, neither of these two interactions alone determines the stability of the G-quadruplex-ligand complexes; rather, it is the result of an intricate interplay between the two. To further examine the nature of the binding, a free energy decomposition analysis at residue level was carried out. The results clearly demonstrate the crucial roles of two hot spot residues (DG4 and DG8) for the binding of ligands to c-KIT G-quadruplex, and highlight the importance of the planar aromatic moiety of ligands in G-quadruplex stabilization via π-π stacking interactions. Our study can assist in the design of new derivatives of 4-quinazolinone with high binding affinity for c-KIT G-quadruplex.