In silico screening of quadruplex-binding ligands

Structured Waters Mediate Small Molecule Binding to G-Quadruplex Nucleic Acids

The role of G-quadruplexes in human cancers is increasingly well-defined. Accordingly, G-quadruplexes can be suitable drug targets and many small molecules have been identified to date as G-quadruplex binders, some using computer-based design methods and co-crystal structures. The role of bound water molecules in the crystal structures of G-quadruplex-small molecule complexes has been analyzed in this study, focusing on the water arrangements in several G-quadruplex ligand complexes. One is the complex between the tetrasubstituted naphthalene diimide compound MM41 and a human intramolecular telomeric DNA G-quadruplex, and the others are in substituted acridine bimolecular G-quadruplex complexes. Bridging water molecules form most of the hydrogen-bond contacts between ligands and DNA in the parallel G-quadruplex structures examined here. Clusters of structured water molecules play essential roles in mediating between ligand side chain groups/chromophore core and G-quadruplex. These clusters tend to be conserved between complex and native G-quadruplex structures, suggesting that they more generally serve as platforms for ligand binding, and should be taken into account in docking and in silico studies.

Determination of a Structural Ensemble Representing the Dynamics of a G-Quadruplex DNA

It is increasingly recognized that the structures and dynamics of G-quadruplex DNA molecules are dictated by their sequences and greatly affected by environmental factors. The core guanine tetrads (G-tetrads) coordinate cations and display a strong conformational rigidity compared with that of the connecting loops. Although long loops linking the G-tetrads are typically disfavored, when present, they provide a striking illustration of the dynamics of short, single-stranded DNA regions. In addition to their role in determining the stability of the G-quadruplex state, these loops are also interesting as potential drug targets. To characterize accurately the dynamics of this DNA state, we apply here the principles of structural ensemble determination developed in the past two decades for protein molecules to DNA molecules. We thus perform extensive molecular dynamics simulations restrained with nuclear magnetic resonance residual dipolar couplings to determine a structural ensemble of the human CEB25 minisatellite G-quadruplex, which contains a connecting loop of nine nucleotides. This structural ensemble displays a wide set of arrangements for the loop and a compact, well-defined G-quadruplex core. Our results show the importance of stacking interactions in the loop and strengthen the ability of the closing base pairs to confer a large thermodynamic stability to the G-quadruplex structure.

Tailoring a lead-like compound targeting multiple G-quadruplex structures

A focused library of analogs of a lead-like G-quadruplex (G4) targeting compound (4), sharing a furobenzoxazine naphthoquinone core and differing for the pendant groups on the N-atom of the oxazine ring, has been here analyzed with the aim of developing more potent and selective ligands. These molecules have been tested vs. topologically different G4s by the G4-CPG assay, an affinity chromatography-based method for screening putative G4 ligands. The obtained results showed that all these compounds were able to bind several G4 structures, both telomeric and extra-telomeric, thus behaving as multi-target ligands, and two of them fully discriminated G4 vs. duplex DNA. Biological assays proved that almost all the compounds produced effective DNA damage, showing marked antiproliferative effects on tumor cells in the low μM range. Combined analysis of the G4-CPG binding assays and biological data led us to focus on compound S4-5, proved to be less cytotoxic than the parent compound 4 on normal cells. An in-depth biophysical characterization of the binding of S4-5 to different G4s showed that the here identified ligand has higher affinity for the G4s and higher ability to discriminate G4 vs. duplex DNA than 4. Molecular docking studies, in agreement with the NMR data, suggest that S4-5 interacts with the accessible grooves of the target G4 structures, giving clues for its increased G4 vs. duplex selectivity.

Hit Identification of a Novel Dual Binder for h-telo/c-myc G-Quadruplex by a Combination of Pharmacophore Structure-Based Virtual Screening and Docking Refinement

It is well known that G-quadruplexes are targets of great interest for their roles in crucial biological processes, such as aging and cancer. Hence, a promising strategy for anticancer drug therapy is the stabilization of these structures by small molecules. We report a high-throughput in silico screening of commercial libraries from several different vendors by means of a combined structure-based pharmacophore model approach followed by docking simulations. The compounds selected by the virtual screening procedure were then tested for their ability to interact with human telomeric G-quadruplex folding by circular dichroism, fluorescence spectroscopy, and fluorescence intercalator displacement. Our approach resulted in the identification of a 13-[(dimethylamino)methyl]-12-hydroxy-8H-benzo[c]indolo[3,2,1-ij][1,5]naphthyridin-8-one derivative as a novel promising stabilizer of G-quadruplex structures within the human telomeric and the c-myc promoter sequences.

Quadruplex Nucleic Acids as Novel Therapeutic Targets

Quadruplex-forming sequences are widely prevalent in human and other genomes, including bacterial ones. These sequences are over-represented in eukaryotic telomeres, promoters, and 5' untranslated regions. They can form quadruplex structures, which may be transient in many situations in normal cells since they can be effectively resolved by helicase action. Mutated helicases in cancer cells are unable to unwind quadruplexes, which are impediments to transcription, translation, or replication, depending on their location within a particular gene. Small molecules that can stabilize quadruplex structures augment these effects and produce cell and proliferation growth inhibition. This article surveys the chemical biology of quadruplexes. It critically examines the major classes of quadruplex-binding small molecules that have been developed to date and the various approaches to discovering selective agents. The challenges of requiring (and achieving) small-molecule targeted selectivity for a particular quadruplex are discussed in relation to the potential of these small molecules as clinically useful therapeutic agents.

Virtual screening and optimization of Type II inhibitors of JAK2 from a natural product library

Amentoflavone has been identified as a JAK2 inhibitor by structure-based virtual screening of a natural product library. In silico optimization using the DOLPHIN model yielded analogues with enhanced potency against JAK2 activity and HCV activity in cellulo. Molecular modeling and kinetic experiments suggested that the analogues may function as Type II inhibitors of JAK2.

Searching for anthranilic acid-based thumb pocket 2 HCV NS5B polymerase inhibitors through a combination of molecular docking, 3D-QSAR and virtual screening

A combination of the following computational methods: (i) molecular docking, (ii) 3-D Quantitative Structure Activity Relationship Comparative Molecular Field Analysis (3D-QSAR CoMFA), (iii) similarity search and (iv) virtual screening using PubChem database was applied to identify new anthranilic acid-based inhibitors of hepatitis C virus (HCV) replication. A number of known inhibitors were initially docked into the "Thumb Pocket 2" allosteric site of the crystal structure of the enzyme HCV RNA-dependent RNA polymerase (NS5B GT1b). Then, the CoMFA fields were generated through a receptor-based alignment of docking poses to build a validated and stable 3D-QSAR CoMFA model. The proposed model can be first utilized to get insight into the molecular features that promote bioactivity, and then within a virtual screening procedure, it can be used to estimate the activity of novel potential bioactive compounds prior to their synthesis and biological tests.

Discovery of deoxyvasicinone derivatives as inhibitors of NEDD8-activating enzyme

NEDD8-activating enzyme (NAE) controls the specific degradation of proteins regulated by cullin-RING ubiquitin E3 ligases, and has been considered as an attractive molecular target for the development of drugs against cancer. A pharmacophore model constructed from a training set of deoxyvasicinone derivatives was used to screen 376 compounds from an analogue database. From the initial screening, the valine-linked deoxyvasicinone derivative 9 and the N-isopropyl-linked deoxyvasicinone derivative 10 emerged as the top scoring candidates. Compounds 9 and 10 showed micromolar potencies in both cell-free and cell-based systems, with selectivity for NAE over the related enzymes SAE and UAE. Molecular modelling analysis suggested that 9 and 10 may inhibit NAE by blocking the ATP-binding domain. Thus, these deoxyvasicinone derivatives could be considered as promising lead molecules for the development of more potent inhibitors of NAE.

Discovery of a small-molecule inhibitor of STAT3 by ligand-based pharmacophore screening

STAT3 modulates the transcription of a wide variety of regulatory genes involved in cell proliferation, differentiation, migration, apoptosis, and other critical cellular functions. Constitutive activation of STAT3 has been detected in a wide spectrum of human malignancies. A pharmacophore model constructed from a training set of STAT3 inhibitors binding to the SH2 domain was used to screen an in-house database of compounds, from which azepine 1 emerged as a top candidate. Compound 1 inhibited STAT3 DNA-binding activity in vitro and attenuated STAT3-directed transcription in cellulo with comparable potency to the well-known STAT3 inhibitor S3I-201. A fluorescence polarization assay revealed that compound 1 targeted the SH2 domain of STAT3. Furthermore, compound 1 inhibited STAT3 phosphorylation in cells without affecting the total expression of STAT3. This study also validates the use of pharmacophore modeling to identify inhibitors of protein-protein interactions.

Two cationic porphyrin isomers showing different multimeric G-quadruplex recognition specificity against monomeric G-quadruplexes

Ligands that can interact specifically with telomeric multimeric G-quadruplexes could be developed as promising anticancer drugs with few side effects related to other G-quadruplex-forming regions. In this paper, a new cationic porphyrin derivative, m-TMPipEOPP, was synthesized and characterized. Its multimeric G-quadruplex recognition specificity under molecular crowding conditions was compared to its isomer p-TMPipEOPP. The slight structural difference accounts for different multimeric G-quadruplex recognition specificity for the two isomers. p-TMPipEOPP can barely discriminate between multimeric and monomeric G-quadruplexes. By contrast, m-TMPipEOPP can bind with multimeric but not with monomeric G-quadruplexes. p-TMPipEOPP might bind to multimeric G-quadruplexes by two modes: sandwich-like end-stacking mode and pocket-dependent intercalative mode. Increasing the pocket size between adjacent two G-quadruplex uints is beneficial for the latter mode. m-TMPipEOPP might bind to multimeric G-quadruplexes by a side binding mode, which confers m-TMPipEOPP with higher multimeric G-quadruplex recognition specificity compared to p-TMPipEOPP. m-TMPipEOPP increases the stability of multimeric G-quadruplex under both dilute and molecular crowding conditions but its G-quadruplex-stabilizing ability is a little weaker than p-TMPipEOPP. These results provide important information for the design of highly specific multimeric G-quadruplex ligands. Another interesting finding is that pocket size is an important factor in determining the stability of multimeric G-quadruplexes.

Structure-based repurposing of FDA-approved drugs as inhibitors of NEDD8-activating enzyme

We report the discovery of an inhibitor of NEDD8-activating enzyme (NAE) by an integrated virtual screening approach. Piperacillin 1 inhibited NAE activity in cell-free and cell-based systems with high selectivity. Furthermore, piperacillin 1 was able to inhibit the degradation of the NAE downstream protein substrate p27(kip1). Our molecular modeling and kinetic studies suggested that this compound may act as a non-covalent ATP-competitive inhibitor of NAE.

Targeting C-myc G-quadruplex: dual recognition by aminosugar-bisbenzimidazoles with varying linker lengths

G-quadruplexes are therapeutically important biological targets. In this report, we present biophysical studies of neomycin-Hoechst 33258 conjugates binding to a G-quadruplex derived from the C-myc promoter sequence. Our studies indicate that conjugation of neomycin to a G-quadruplex binder, Hoechst 33258, enhances its binding. The enhancement in G-quadruplex binding of these conjugates varies with the length and composition of the linkers joining the neomycin and Hoechst 33258 units.

Mechanistic studies for the role of cellular nucleic-acid-binding protein (CNBP) in regulation of c-myc transcription

BACKGROUND

Guanine-rich sequence of c-myc nuclease hypersensitive element (NHE) III1 is known to fold in G-quadruplex and subsequently serves as a transcriptional silencer. Cellular nucleic-acid-binding protein (CNBP), a highly conserved zinc-finger protein with multiple biological functions, could bind to c-myc NHE III1 region, specifically to the single strand G-rich sequence.

METHODS

In the present study, a variety of methods, including cloning, expression and purification of protein, EMSA, CD, FRET, Ch-IP, RNA interference, luciferase reporter assay, SPR, co-immunoprecipitation, and co-transfection, were applied to investigate the mechanism for the role of CNBP in regulating c-myc transcription.

RESULTS

We found that human CNBP specifically bound to the G-rich sequence of c-myc NHE III1 region both in vitro and in cellulo, and subsequently promoted the formation of G-quadruplex. CNBP could induce a transient decrease followed by an increase in c-myc transcription in vivo. The interaction of CNBP with NM23-H2 was responsible for the increase of c-myc transcription.

CONCLUSIONS

Based on above experimental results, a new mechanism, involving G-quadruplex related CNBP/NM23-H2 interaction, for the regulation of c-myc transcription was proposed.

GENERAL SIGNIFICANCE

These findings indicated that the regulation of c-myc transcription through NHE III1 region might be governed by mechanisms involving complex protein-protein interactions, and suggested a new possibility of CNBP as a potential anti-cancer target based on CNBP's biological function in c-myc transcription.

Comparing the suitability of autodock, gold and glide for the docking and predicting the possible targets of Ru(II)-based complexes as anticancer agents

In cancer chemotherapy, metal-based complexes have been recognized as the most promising means of inhibiting cancer growth due to the successful application of cis-platin and its derivatives above many of the existing organic anticancer agents. The limitations in their rational design can be traced to the complexity of the mechanism of their operations, lack of proper knowledge of their targets and lack of force fields in docking packages to appropriately define the metal centre of the organometallic complexes. In this paper, some of the promising anticancer complexes of Ru(II) such as the rapta-based complexes formulated as [Ru(η6-p-cymene)L2(pta)] and those with unusual ligands are considered. CatB and kinases which have been experimentally confirmed as possible targets of the complexes are also predicted by the three methods as one of the most targeted receptors while TopII and HDAC7 are predicted by two and one of the methods as best targets. The interesting features of the binding of the complexes show that some of the complexes preferentially target specific macromolecules than the others, which is an indication of their specificity and possibility of their therapeutic combination without severe side effects that may come from competition for the same target. Also, introduction of unusual ligands is found to significantly improve the activities of most of the complexes studied. Strong correlations are observed for the predicted binding sites and the orientation of the complexes within the binding site by the three methods of docking. However there are disparities in the ranking of the complexes by the three method of docking, especially that of Glide.

Hit identification of IKKβ natural product inhibitor

Background

The nuclear factor-κB (NF-κB) proteins are a small group of heterodimeric transcription factors that play an important role in regulating the inflammatory, immune, and apoptotic responses. NF-κB activity is suppressed by association with the inhibitor IκB. Aberrant NF-κB signaling activity has been associated with the development of cancer, chronic inflammatory diseases and auto-immune diseases. The IKK protein complex is comprised of IKKα, IKKβ and NEMO subunits, with IKKβ thought to play the dominant role in modulating NF-κB activity. Therefore, the discovery of new IKKβ inhibitors may offer new therapeutic options for the treatment of cancer and inflammatory diseases.

Results

A structure-based molecular docking approach has been employed to discover novel IKKβ inhibitors from a natural product library of over 90,000 compounds. Preliminary screening of the 12 highest-scoring compounds using a luciferase reporter assay identified 4 promising candidates for further biological study. Among these, the benzoic acid derivative (1) showed the most promising activity at inhibiting IKKβ phosphorylation and TNF-α-induced NF-κB signaling in vitro.

Conclusions

In this study, we have successfully identified a benzoic acid derivative (1) as a novel IKKβ inhibitor via high-throughput molecular docking. Compound 1 was able to inhibit IKKβ phosphorylation activity in vitro, and block IκBα protein degradation and subsequent NF-κB activation in human cells. Further in silico optimization of the compound is currently being conducted in order to generate more potent analogues for biological tests.

Drug repositioning by structure-based virtual screening

Approved drugs have favourable or validated pharmacokinetic properties and toxicological profiles, and the repositioning of existing drugs for new indications can potentially avoid expensive costs associated with early-stage testing of the hit compounds. In recent years, technological advances in virtual screening methodologies have allowed medicinal chemists to rapidly screen drug libraries for therapeutic activity against new biomolecular targets in a cost-effective manner. This review article outlines the basic principles and recent advances in structure-based virtual screening and highlights the powerful synergy of in silico techniques in drug repositioning as demonstrated in several recent reports.

Molecular dynamics simulations of G-DNA and perspectives on the simulation of nucleic acid structures

The article reviews the application of biomolecular simulation methods to understand the structure, dynamics and interactions of nucleic acids with a focus on explicit solvent molecular dynamics simulations of guanine quadruplex (G-DNA and G-RNA) molecules. While primarily dealing with these exciting and highly relevant four-stranded systems, where recent and past simulations have provided several interesting results and novel insight into G-DNA structure, the review provides some general perspectives on the applicability of the simulation techniques to nucleic acids.