NB-506, an indolocarbazole topoisomerase I inhibitor, binds preferentially to triplex DNA

A novel competition dialysis method was used to study the structural selectivity of the nucleic acid binding of NB-506, a promising indolocarbazole anticancer agent. A pronounced preference for NB-506 binding to the DNA triplex poly [dA]:(poly[dT])(2) was observed among potential binding to 12 different nucleic acid structures and sequences. Structures included in the assay ranged from single-stranded DNA, through a variety of right-handed DNA duplexes, to multistranded triplex and tetraplex forms. RNA and left-handed Z DNA were also included in the assay. The preferential binding to triplex was confirmed by UV melting experiments. The novel and unexpected structural selectivity shown by NB-506 may arise from a complementary shape between its extended aromatic ring system and the planar triplex stack.     

Sequence and structural selectivity of nucleic acid binding ligands

The sequence and structural selectivity of 15 different DNA binding agents was explored using a novel, thermodynamically rigorous, competition dialysis procedure. In the competition dialysis method, 13 different nucleic acid structures were dialyzed against a common ligand solution. More ligand accumulated in the dialysis tube containing the structural form with the highest ligand binding affinity. DNA structural forms included in the assay ranged from single-stranded forms, through a variety of duplex forms, to multistranded triplex and tetraplex forms. Left-handed Z-DNA, RNA, and a DNA-RNA hybrid were also represented. Standard intercalators (ethidium, daunorubicin, and actinomycin D) served as control compounds and were found to show structural binding preferences fully consistent with their previously published behavior. Standard groove binding agents (DAPI, distamycin, and netropsin) showed a strong preference for AT-rich duplex DNA forms, along with apparently strong binding to the poly(dA)-[poly(dT)](2) triplex. Thermal denaturation studies revealed the apparent triplex binding to be complex, and perhaps to result from displacement of the third strand. Putative triplex (BePI, coralyne, and berberine) and tetraplex [H(2)TmPyP, 5,10,15, 20-tetrakis[4-(trimethylammonio)phenyl]-21H,23H-porphine, and N-methyl mesoporphyrin IX] selective agents showed in many cases less dramatic binding selectivity than anticipated from published reports that compared their binding to only a few structural forms. Coralyne was found to bind strongly to single-stranded poly(dA), a novel and previously unreported interaction. Finally, three compounds (berenil, chromomycin A, and pyrenemethylamine) whose structural preferences are largely unknown were examined. Pyrenemethylamine exhibited an unexpected and unprecedented preference for duplex poly(dAdT).     

Interaction of 9-O-(ω-amino) alkyl ether berberine analogs with poly(dT)·poly(dA)*poly(dT) triplex and poly(dA)·poly(dT) duplex: a comparative study

Isoquinoline alkaloids and their analogs represent an important class of molecules for their broad range of clinical and pharmacological utility. These compounds are of current interest owing to their low toxicity and excellent chemo preventive properties. These alkaloids can play important role in stabilising the nucleic acid triple helices. The present study has focused on the interaction of five 9-O-(ω-amino) alkyl ether berberine analogs with the DNA triplex poly(dT)·poly(dA)*poly(dT) and the parent duplex poly(dA)·poly(dT) studied using various biophysical techniques. Scatchard analysis of the spectral data indicated that the analogs bind both to the duplex and triplex in a non-cooperative manner in contrast to the cooperative binding of berberine to the DNA triplex. Strong intercalative binding to the DNA triplex structure was revealed from ferrocyanide quenching, fluorescence polarization and viscosity results. Thermal melting studies demonstrated higher stabilization of the Hoogsteen base paired third strand of the DNA triplex compared to the Watson–Crick strand. Circular dichroism studies suggested a stronger perturbation of the DNA triplex conformation by the alkaloid analogs compared to the duplex. The binding was entropy-driven in each case and the entropy contribution to free energy increased as the length of the alkyl side chain increased. The analogs exhibited stronger binding affinity to the triple helical structure compared to the parent double helical structure.     

Thermodynamics of nucleic acid "shape readout" by an aminosugar

Recognition of nucleic acids is important for our understanding of nucleic acid structure as well as for our understanding of nucleic acid-protein interactions. In addition to the direct readout mechanisms of nucleic acids such as H-bonding, shape recognition of nucleic acids is being increasingly recognized as playing an equally important role in DNA recognition. Competition dialysis, UV, flourescent intercalator displacement (FID), computational docking, and calorimetry studies were conducted to study the interaction of neomycin with a variety of nucleic acid conformations (shapes). At pH 5.5, the results suggest the following. (1) Neomycin binds three RNA structures [16S A site rRNA, poly(rA)·poly(rA), and poly(rA)·poly(rU)] with high affinities (K(a) ~ 10(7) M(-1)). (2) The binding of neomycin to A-form GC-rich oligomer d(A(2)G(15)C(15)T(2))(2) has an affinity comparable to those of RNA structures. (3) The binding of neomycin to DNA·RNA hybrids shows a 3-fold variance that can be attributed to their structural differences [for poly(dA)·poly(rU), K(a) = 9.4 × 10(6) M(-1), and for poly(rA)·poly(dT), K(a) = 3.1 × 10(6) M(-1)]. (4) The interaction of neomycin with DNA triplex poly(dA)·2poly(dT) yields a binding affinity (K(a)) of 2.4 × 10(5) M(-1). (5) Poly(dA-dT)(2) shows the lowest association constant for all nucleic acids studied (K(a) < 10(5)). (6) Neomycin binds to G-quadruplexes with K(a) values of ~10(4)-10(5) M(-1). (7) Computational studies show that the decrease in major groove width in the B to A transition correlates with increasing neomycin affinity. Neomycin's affinity for various nucleic acid structures can be ranked as follows: RNAs and GC-rich d(A(2)G(15)C(15)T(2))(2) structures > poly(dA)·poly(rU) > poly(rA)·poly(dT) > T·A-T triplex, G-quadruplex, B-form AT-rich, or GC-rich DNA sequences. The results illustrate the first example of a small molecule-based "shape readout" of different nucleic acid conformations.     

In silico hit optimization toward AKT inhibition: fragment-based approach,molecular docking and molecular dynamics study

Abstract

Protein kinase B also known as AKT is a cardinal node in different signaling pathways that regulates diverse cell processes. AKT has three isoforms that share high homology. Hyperactivation of each isoform is related with different types of cancer. This work describes the computational search for new inhibitors using a hit optimization process of the previously reported AKT pan inhibitor, a 2,4,6-trisubstituted pyridine. A database of new molecules was proposed using a variant of fragment-based docking methodology and previous reported considerations. Molecular docking followed by molecular dynamics studies were performed to select the best compounds and analyze their behavior. Protein–ligand complexes energy was calculated using molecular mechanics Poisson–Boltzmann surface area protocol. Further, proposed molecules were compared with the ChEMBL database of compounds assayed against AKT. Data analysis leads to determine the structural requirements necessary for a favorable interaction of the proposed ligands with the AKT pocket. Molecular dynamics data suggested that the pK_a of the ligands is important for the stability in the AKT pocket. Molecular similarity analysis shows that proposed ligands have not been previously reported. Thus, ligands with high docking scores and favorable behavior on molecular dynamics simulations are proposed as potential AKT inhibitors.     

Competition dialysis: a method for the study of structural selective nucleic acid binding

Competition dialysis is a powerful new tool for the discovery of ligands that bind to nucleic acids with structural- or sequence-selectivity. The method is based on firm thermodynamic principles and is simple to implement. In the competition dialysis experiment, an array of nucleic acid structures and sequences is dialyzed against a common test ligand solution. After equilibration, the amount of ligand bound to each structure or sequence is determined by absorbance or fluorescence measurements. Since all structures and sequences are in equilibrium with the same free ligand concentration, the amount bound is directly proportional to the ligand binding affinity. Competition dialysis thus provides a direct and quantitative measure of selectivity, and unambiguously identifies which of the samples within the array are preferred by a particular ligand. We describe here the third generation implementation of the method, in which competition dialysis was adapted for use in a 96-well plate format. In this format, we have been able to greatly expand the array of nucleic acid structures studied, and now can routinely study the interactions of a ligand of interest with 46 different structures and sequences.     

Prediction of estrogen receptor β ligands potency and selectivity by docking and MM-GBSA scoring methods using three different scaffolds

This study aimed to identify the docking and molecular mechanics-generalized born surface area (MM-GBSA) re-scoring parameters which can correlate the binding affinity and selectivity of the ligands towards oestrogen receptor β (ERβ). Three different series of ERβ ligands were used as dataset and the compounds were docked against ERβ (protein data bank (PDB) ID: 1QKM) using Glide and ArgusLab. Glide docking showed superior results when compared with ArgusLab. Docked poses were then rescored using Prime-MM-GBSA to calculate free energy binding. Correlations were made between observed activities of ERβ ligands with computationally predicted values from docking, binding energy parameters. ERβ ligands experimental binding affinity/selectivity did not correlate well with Glide and ArgusLab score. Whereas calculated Glide energy (coulomb-van der Waal interaction energy) correlated significantly with binding affinity of ERβ ligands (r²?=?0.66). MM-GBSA re-scoring showed correlation of r²?=?0.74 with selectivity of ERβ ligands. These results will aid the discovery of novel ERβ ligands with isoform selectivity.     

结合分子相似性、药效团和分子对接筛选新的HIV-1蛋白酶抑制剂

结合分子相似性、药效团和分子对接建立兼顾计算效率和预测准确度的HIV-1蛋白酶抑制剂筛选方法。首先通过对现有HIV-1蛋白酶抑制剂分子进行相似性分析,选取代表性的HIV-1蛋白酶抑制剂作为模板分子,构建和优化药效团模型,并从1万个化合物中优先筛选出500个化合物。而后采用分子对接方法进一步考察化合物与HIV-1蛋白酶结合情况,得到4个新的活性候选化合物,并进行其结合自由能计算和抗突变性分析。结果表明新候选化合物ST025723和HIV-1蛋白酶表现出较好的相互作用和抗突变性,具有深入研究的价值,同时也证明分子相似性、药效团和分子对接相结合能够快速有效地发现新颖活性候选化合物。     

Synthesis and anti-inflammatory activity of new arylidene-thiazolidine-2,4-diones as PPARγ ligands

Eight new 5-arylidene-3-benzyl-thiazolidine-2,4-diones with halide groups on their benzyl rings were synthesized and assayed in vivo to investigate their anti-inflammatory activities. These compounds showed considerable biological efficacy when compared to rosiglitazone, a potent and well-known agonist of PPARγ, which was used as a reference drug. This suggests that the substituted 5-arylidene and 3-benzylidene groups play important roles in the anti-inflammatory properties of this class of compounds. Docking studies with these compounds indicated that they exhibit specific interactions with key residues located in the site of the PPARγ structure, which corroborates the hypothesis that these molecules are potential ligands of PPARγ. In addition, competition binding assays showed that four of these compounds bound directly to the ligand-binding domain of PPARγ, with reduced affinity when compared to rosiglitazone. An important trend was observed between the docking scores and the anti-inflammatory activities of this set of molecules. The analysis of the docking results, which takes into account the hydrophilic and hydrophobic interactions between the ligands and the target, explained why the 3-(2-bromo-benzyl)-5-(4-methanesulfonyl-benzylidene)-thiazolidine-2,4-dione compound had the best activity and the best docking score. Almost all of the stronger hydrophilic interactions occurred between the substituted 5-arylidene group of this compound and the residues of the binding site.     

Ascididemin and meridine stabilise G-quadruplexes and inhibit telomerase in vitro

Ascididemin and Meridine are two marine compounds with pyridoacridine skeletons known to exhibit interesting antitumour activities. These molecules have been reported to behave like DNA intercalators. In this study, dialysis competition assay and mass spectrometry experiments were used to determine the affinity of ascididemin and meridine for DNA structures among duplexes, triplexes, quadruplexes and single-strands. Our data confirm that ascididemin and meridine interact with DNA but also recognize triplex and quadruplex structures. These molecules exhibit a significant preference for quadruplexes over duplexes or single-strands. Meridine is a stronger quadruplex ligand and therefore a stronger telomerase inhibitor than ascididemin (IC50=11 and >80 muM, respectively in a standard TRAP assay).     

Triplex hydration: nanosecond molecular dynamics simulation of the solvated triplex formed by mixed sequences

A theoretical model for the hydration pattern and motion of ions around the triple helical DNA with mixed sequences d(GACTGGTGAC)d(GTCACCAGTC)*d(GACTGGTGAC) in solution, during MD simulation, using the particle mesh Ewald sum method, is elaborated here. The AMBER 5.0 force field has been used during the simulation in solvent. The simulation studies support a dynamically stable atmosphere around the DNA triplex in solution over the entire length of the trajectory. The results have been compared with Hoogsteen triplexes and examined in the context of the observed behaviour of hydration in crystallographic data of duplexes. The dynamical organization of counterions and water molecules around the triplex formed by mixed sequences is described here. It has been observed that cations prefer to bind between two adjoining purines of the second and the third strands. The idea of localized complexes (mobile counterions in unspecific electronegative pockets around the DNA triplex with water molecules) may have important implications for understanding the specificity of the interactions of nucleic acids with proteins and other ligands.     

Binding of netropsin to DNA in complexes with polypeptides containing repetitive lysine sequences

The interaction of the antibiotic netropsin with calf thymus DNA, T4 DNA and poly(dA-dT) . poly(dA-dT) in complexes with sequential polypeptides containing repetitive lysine sequences and histone H1 was investigated using circular dichroism spectroscopy and equilibrium dialysis. Both soluble DNA-polypeptide complexes and insoluble complexes showed binding of netropsin. The possibility of displacement of polypeptides from DNA binding sites by competition with netropsin molecules was eliminated by experiments using 14C-labelled polypeptides. From the analysis of CD titration behavior as well as from the results of equilibrium dialysis studies it follows that netropsin does not compete with polypeptides for DNA binding sites, which suggests that these two ligands occupy different sites. Various explanations for minor differences in the CD behavior of the bound netropsin in the saturation region are also discussed.     

Structure-based virtual screening and biological evaluation of potent and selective ADAM12 inhibitors

We describe a series of potent and selective inhibitors of ADAM12 that were discovered using computational screening of a focused virtual library. The initial structure-based virtual screening selected 64 compounds from a 3D database of 67,062 molecules. Being evaluated by a cell-based ADAM12 activity assay, compounds 5, 11, 14, 16 were further identified as the potent and selective inhibitors of ADAM12 with low nanomolar IC₅₀ values. The mechanism underlying the potency and selectivity of a representative compound, 5, was investigated through molecular docking studies.     

Structural comparison of Mtb-DHFR and h-DHFR for design,synthesis and evaluation of selective non-pteridine analogues as antitubercular agents

Tuberculosis is an infectious disease that affects millions of population every year. Mtb-DHFR is a validated target that is vital for nucleic acids biosynthesis and therefore DNA formation and cell replication. This paper report identification and synthesis of novel compounds for selective inhibition of Mtb-DHFR and unleash the selective structural features necessary to inhibit the same. Virtual screening of databases was carried out to identify novel compounds on the basis of difference between the binding pockets of the two proteins. Consensus docking was performed to improve upon the results and best ten hits were selected. Hit 1 was subjected to analogues design and the analogues were docked against Mtb-DHFR. From the docking results 11 compounds were selected for synthesis and biological assay against H₃₇Rv. Most potent compound (IND-07) was tested for selectivity using enzymatic assay against Mtb-DHFR and h-DHFR. The compounds were found to have good inhibitory activity (25–200 µM) against H₃₇Rv and in enzyme assay against Mtb-DHFR and h-DHFR the compound was found selective towards Mtb-DHFR with selectivity index of 6.53. This work helped to identify indole moiety as novel scaffold for development of novel selective Mtb-DHFR inhibitors as antimycobacterial agents.     

Interaction of trivaline with single-stranded polyribonucleotides]

Binding of tripeptide H-Val3-(NH)2-Dns (TVP) to polyribonucleotides was studied by fluorescence methods, circular and flow linear dichroism, equilibrium dialysis and electron microscopy. It was found that TVP binds to poly(U) in monomer, dimer and tetramer forms with binding constants of about 10(3), 40, 18.10(4) M, respectively. The cooperativity parameter for peptide dimer binding is 2000. The peptide forms tetramer complexes with poly(A), poly(C), poly(G) also. The formation of a complex between the peptide tetramer and nucleic acid is accompanied by a significant increase in the fluorescence intensity. The cooperative binding of TVP dimers to poly(U), poly(A), poly(C) is accompanied by a dramatic decrease in the flexibility of polynucleotide chains. However, it has a small effect (if any) on the flexibility of the poly(G) chain. The observed similarity of thermodynamic, optical and hydrodynamic++ properties of TVP complexes with single-stranded and double-stranded nucleic acids may reflect a similarity in the geometries of peptide complexes with nucleic acids. Electron microscopy studies show that peptide binding to poly(U) and dsDNA leads to compactization of the nucleic acids caused by interaction between the peptide tetramers bound to a nucleic acid. At the first stage of the compactization process the well-organized rod-like particles are formed, each consisting of one or more single-stranded polynucleotide fibers. Increasing the peptide concentration stimulates a side-by-side association and folding of the rods with the formation of macromolecular "leech-like" structures with the thickness of 20-50 nm.     

Counterion association with native and denatured nucleic acids: an experimental approach

The melting temperature of the poly(dA) . poly(dT) double helix is exquisitely sensitive to salt concentration, and the helix-to-coil transition is sharp. Modern calorimetric instrumentation allows this transition to be detected and characterized with high precision at extremely low duplex concentrations. We have taken advantage of these properties to show that this duplex can be used as a sensitive probe to detect and to characterize the influence of other solutes on solution properties. We demonstrate how the temperature associated with poly(dA) . poly(dT) melting can be used to define the change in bulk solution cation concentration imparted by the presence of other duplex and triplex solutes, in both their native and denatured states. We use this information to critically evaluate features of counterion condensation theory, as well as to illustrate "crosstalk" between different, non-contacting solute molecules. Specifically, we probe the melting of a synthetic homopolymer, poly(dA) . poly(dT), in the presence of excess genomic salmon sperm DNA, or in the presence of one of two synthetic RNA polymers (the poly(rA) . poly(rU) duplex or the poly(rU) . poly(rA) . poly(rU) triplex). We find that these additions cause a shift in the melting temperature of poly(dA) . poly(dT), which is proportional to the concentration of the added polymer and dependent on its conformational state (B versus A, native versus denatured, and triplex versus duplex). To a first approximation, the magnitude of the observed tm shift does not depend significantly on whether the added polymer is RNA or DNA, but it does depend on the number of strands making up the helix of the added polymer. We ascribe the observed changes in melting temperature of poly(dA) . poly(dT) to the increase in ionic strength of the bulk solution brought about by the presence of the added nucleic acid and its associated counterions. We refer to this communication between non-contacting biopolymers in solution as solvent-mediated crosstalk. By comparison with a known standard curve of tm versus log[Na+] for poly(dA) . poly(dT), we estimate the magnitude of the apparent change in ionic strength resulting from the presence of the bulk nucleic acid, and we compare these results with predictions from theory. We find that current theoretical considerations correctly predict the direction of the t(m) shift (the melting temperature increases), while overestimating its magnitude. Specifically, we observe an apparent increase in ionic strength equal to 5% of the concentration of the added duplex DNA or RNA (in mol phosphate), and an additional apparent increase of about 9.5 % of the nucleic acid concentration (mol phosphate) upon denaturation of the added DNA or RNA, yielding a total apparent increase of 14.5 %. For the poly(rU) . poly(rA) . poly(rU) triplex, the total apparent increase in ionic strength corresponds to about 13.6% of the amount of added triplex (moles phosphate). The effect we observe is due to coupled equilibria between the solute molecules mediated by modulations in cation concentration induced by the presence and/or the transition of one of the solute molecules. We note that our results are general, so one can use a different solute probe sensitive to proton binding to characterize subtle changes in solution pH induced by the presence of another solute in solution. We discuss some of the broader implications of these measurements/results in terms of nucleic acid melting in multicomponent systems, in terms of probing counterion environments, and in terms of potential regulatory mechanisms.     

An In Silico Analysis of the Binding Modes and Binding Affinities of Small Molecule Modulators of PDZ-Peptide Interactions

Inhibitors of PDZ-peptide interactions have important implications in a variety of biological processes including treatment of cancer and Parkinson’s disease. Even though experimental studies have reported characterization of peptidomimetic inhibitors of PDZ-peptide interactions, the binding modes for most of them have not been characterized by structural studies. In this study we have attempted to understand the structural basis of the small molecule-PDZ interactions by in silico analysis of the binding modes and binding affinities of a set of 38 small molecules with known K_i or K_d values for PDZ2 and PDZ3 domains of PSD-95 protein. These two PDZ domains show differential selectivity for these compounds despite having a high degree of sequence similarity and almost identical peptide binding pockets. Optimum binding modes for these ligands for PDZ2 and PDZ3 domains were identified by using a novel combination of semi-flexible docking and explicit solvent molecular dynamics (MD) simulations. Analysis of the binding modes revealed most of the peptidomimectic ligands which had high K_i or K_d moved away from the peptide binding pocket, while ligands with high binding affinities remained in the peptide binding pocket. The differential specificities of the PDZ2 and PDZ3 domains primarily arise from differences in the conformation of the loop connecting βB and βC strands, because this loop interacts with the N-terminal chemical moieties of the ligands. We have also computed the MM/PBSA binding free energy values for these 38 compounds with both the PDZ domains from multiple 5 ns MD trajectories on each complex i.e. a total of 228 MD trajectories of 5 ns length each. Interestingly, computational binding free energies show good agreement with experimental binding free energies with a correlation coefficient of approximately 0.6. Thus our study demonstrates that combined use of docking and MD simulations can help in identification of potent inhibitors of PDZ-peptide complexes.     

Virtual screening of small molecules databases for discovery of novel PARP-1 inhibitors: combination of in silico and in vitro studies

Poly(ADP-ribose) polymerase-1 (PARP-1) enzyme has critical roles in DNA replication repair and recombination. Thus, PARP-1 inhibitors play an important role in the cancer therapy. In the current study, we have performed combination of in silico and in vitro studies in order to discover novel inhibitors against PARP-1 target. Structure-based virtual screening was carried out for an available small molecules database. A total of 257,951 ligands from Otava database were screened at the binding pocket of PARP-1 using high-throughput virtual screening techniques. Filtered structures based on predicted binding energy results were then used in more sophisticated molecular docking simulations (i.e. Glide/standard precision, Glide/XP, induced fit docking – IFD, and quantum mechanics polarized ligand docking – QPLD). Potential high binding affinity compounds that are predicted by molecular simulations were then tested by in vitro methods. Computationally proposed compounds as PARP-1 inhibitors (Otava Compound Codes: 7111620047 and 7119980926) were confirmed by in vitro studies. In vitro results showed that compounds 7111620047 and 7119980926 have IC₅₀ values of 0.56 and 63 μM against PARP-1 target, respectively. The molecular mechanism analysis, free energy perturbation calculations using long multiple molecular dynamics simulations for the discovered compounds which showed high binding affinity against PARP-1 enzyme, as well as structure-based pharmacophore development (E-pharmacophore) studies were also studied.     

Controlling nucleic acid secondary structure by intercalation: effects of DNA strand length on coralyne-driven duplex disproportionation

Small molecules that intercalate in DNA and RNA are powerful agents for controlling nucleic acid structural transitions. We recently demonstrated that coralyne, a small crescent-shaped molecule, can cause the complete and irreversible disproportionation of duplex poly(dA)·poly(dT) into triplex poly(dA)·poly(dT)·poly(dT) and a poly(dA) self- structure. Both DNA secondary structures that result from duplex disproportionation are stabilized by coralyne intercalation. In the present study, we show that the kinetics and thermodynamics of coralyne-driven duplex disproportionation strongly depend on oligonucleotide length. For example, disproportionation of duplex (dA)₁₆·(dT)₁₆ by coralyne reverts over the course of hours if the sample is maintained at 4°C. Coralyne-disproportioned (dA)₃₂· (dT)₃₂, on the other hand, only partially reverts to the duplex state over the course of days at the same temperature. Furthermore, the equilibrium state of a (dA)₁₆·(dT)₁₆ sample in the presence of coralyne at room temperature contains three different secondary structures [i.e. duplex, triplex and the (dA)₁₆ self-structure]. Even the well-studied process of triplex stabilization by coralyne binding is found to be a length-dependent phenomenon and more complicated than previously appreciated. Together these observations indicate that at least one secondary structure in our nucleic acid system [i.e. duplex, triplex or (dA)_n self-structure] binds coralyne in a length-dependent manner.