Structure-Based Virtual Screening: An Application to Human Topoisomerase II α

Investigation of reactive properties, adsorption on fullerene, DFT, molecular dynamics simulation of an anthracene derivative targeting dihydrofolate reductase and human dUTPase

Anthracenes are aromatic compounds with flexible structure and reactivity which are of great interest to theoretical and experimental chemists. Theoretical investigations of 1,4-dihydroxy-5,8-bis[2-(2-hydroxyethylamino)ethylamino]anthracene-9,10-dione (Mitoxantrone) (DDEA) based on density functional theory, molecular dynamics and adsorption on fullerene are reported in the present research. The suitable situation for adsorption with fullerene (C60) is the cyclohex-2-ene-1,4-dione ring of DDEA. Selected quantum-molecular descriptors have been calculated to predict the most reactive sites of the DDEA molecule. Interactions of DDEA with water have been studied using MD simulations. MD simulations were also used to study solubility parameter, a significant quantity for the development of pharmaceutical formulations. The affinity of DDEA on human dihydrofolate reductase and deoxyuridine triphosphatase enzymes was investigated by MD simulation of the protein-ligand complex obtained by molecular docking study.Communicated by Ramaswamy H. Sarma.

Synthesis of 7-hydroxy-6H-naphtho[2,3-c]coumarin via a TsOH-mediated tandem reaction

A concise and efficient method for the synthesis of 7-hydroxy-6H-naphtho[2,3-c]coumarin using available 1-(2-hydroxyphenyl)-2-phenylethanone and Meldrum's acid has been developed. This transformation involved a tandem aldol reaction/lactonization/Friedel-Crafts reaction to form a lactone ring and a benzene ring. It showed high atom economy with water and acetone as the byproducts. Mechanism studies demonstrated two roles of Meldrum's acid: (i) as the reagent for the tandem reaction, and (ii) as the catalyst for the Friedel-Crafts reaction. Moreover, the hydroxyl group of 7-hydroxy-6H-naphtho[2,3-c]coumarin was further functionalized efficiently by arylethynyl, aryl, and cyano groups to furnish D-π-A compounds with excellent fluorescence emissions (Φ_F = 0.14-0.78).

Analysis of plant-derived phytochemicals as anti-cancer agents targeting cyclin dependent kinase-2, human topoisomerase IIa and vascular endothelial growth factor receptor-2

Cancer is caused by a variety of pathways, involving numerous types of enzymes. Among them three enzymes i.e. Cyclin-dependent kinase-2 (CDK-2), Human topoisomerase IIα, and Vascular Endothelial Growth Factor Receptor-2 (VEGFR-2) are three of the most common enzymes that are involved in the cancer development. Although many chemical drugs are already available in the market for cancer treatment, plant sources are known to contain a wide variety of agents that are proved to possess potential anticancer activity. In this experiment, total thirty phytochemicals were analyzed against the mentioned three enzymes using different tools of bioinformatics and in silico biology like molecular docking study, drug likeness property experiment, ADME/T test, PASS prediction, and P450 site of metabolism prediction as well as DFT calculation to determine the three best ligands among them that have the capability to inhibit the mentioned enzymes. From the experiment, Epigallocatechin gallate was found to be the best ligand to inhibit CDK-2, Daidzein showed the best inhibitory activities towards the Human topoisomerase IIα, and Quercetin was predicted to be the best agent against VEGFR-2. They were also predicted to be quite safe and effective agents to treat cancer. However, more in vivo and in vitro analyses are required to finally confirm their safety and efficacy in this regard.

Inhibition of DNA Topoisomerase Type IIα (TOP2A) by Mitoxantrone and Its Halogenated Derivatives: A Combined Density Functional and Molecular Docking Study

In this study, mitoxantrone and its halogenated derivatives have been designed by density functional theory (DFT) to explore their structural and thermodynamical properties. The performance of these drugs was also evaluated to inhibit DNA topoisomerase type IIα (TOP2A) by molecular docking calculation. Noncovalent interactions play significant role in improving the performance of halogenated drugs. The combined quantum and molecular mechanics calculations revealed that CF₃ containing drug shows better preference in inhibiting the TOP2A compared to other modified drugs.

Characterisation of the DNA gyrase from the thermophilic eubacterium Thermus thermophilus

DNA gyrase is a type IIA topoisomerase found in bacteria but not in humans. The enzyme is required for bacterial DNA replication and transcription, and is an important antibacterial target that is sensitive to the widely-used fluoroquinolone drugs. Due to the emergence of fluoroquinolone resistance, the discovery of new classes of drugs that target DNA gyrase is urgent. The DNA gyrase holoenzyme is a heterodimer of subunit pairs (A2B2). The 90 kDa A subunits bind, cleave, and rejoin double stranded DNA. The enzyme introduces negative supercoils into closed circular bacterial DNA using ATP hydrolysis catalysed by the 70 kDa B subunits. Subdomains of DNA gyrase subunits have been crystallised for structural analysis and the resulting models used to improve drugs that target the DNA binding region and active site. While crystal structures are available for topoisomerase IV complexes with cleaved DNA, there is none for the complete DNA gyrase complex with substrate DNA bound. Thermophiles offer significant advantages in obtaining stable enzymes for structural and functional studies. In order to develop a capability for drug screening and structure-directed drug discovery we have reconstituted a functional and drug-sensitive DNA gyrase complex using heterologously expressed subunits from the thermophile Thermus thermophilus.

Ellagic acid and polyhydroxylated urolithins are potent catalytic inhibitors of human topoisomerase II: an in vitro study

Ellagic acid (EA), a natural polyphenol abundant in fruits and common in our diet, is under intense investigation for its chemopreventive activity resulting from multiple effects. EA inhibits topoisomerase II, but the effects on the human enzyme of urolithins, its monolactone metabolites, are not known. Therefore, the action of several synthetic urolithins toward topoisomerases II was evaluated, showing that polyhydroxylated urolithins, EA, and EA-related compounds are potent inhibitors of the α and β isoforms of human topoisomerase II at submicromolar concentrations. Competition tests demonstrate a dose-dependent relationship between ATP and the inhibition of the enzyme. Docking experiments show that the active compounds bind the ATP pocket of the human enzyme, thus supporting the hypothesis that EA and polyhydroxylated urolithins act as ATP-competitive inhibitors of human topoisomerase II.

Tyrosine 656 in topoisomerase IIβ is important for the catalytic activity of the enzyme: Identification based on artifactual +80-Da modification at this site

Topoisomerase (topo) II catalyzes topological changes in DNA. Although both human isozymes, topo IIα and β are phosphorylated, site-specific phosphorylation of topo IIβ is poorly characterized. Using LC-MS/MS analysis of topo IIβ, cleaved with trypsin, Arg C or cyanogen bromide (CNBr) plus trypsin, we detected four +80-Da modified sites: tyr656, ser1395, thr1426 and ser1545. Phosphorylation at ser1395, thr1426 and ser1545 was established based on neutral loss of H(3) PO(4) (-98 Da) in the CID spectra and on differences in 2-D-phosphopeptide maps of (32) P-labeled wild-type (WT) and S1395A or T1426A/S1545A mutant topo IIβ. However, phosphorylation at tyr656 could not be verified by 2-D-phosphopeptide mapping of (32) P-labeled WT and Y656F mutant protein or by Western blotting with phosphotyrosine-specific antibodies. Since the +80-Da modification on tyr656 was observed exclusively during cleavage with CNBr and trypsin, this modification likely represented bromination, which occurred during CNBr cleavage. Re-evaluation of the CID spectra identified +78/+80-Da fragment ions in CID spectra of two peptides containing tyr656 and tyr711, confirming bromination. Interestingly, mutation of only tyr656, but not ser1395, thr1326 or ser1545, decreased topo IIβ activity, suggesting a functional role for tyr656. These results, while identifying an important tyrosine in topo IIβ, underscore the importance of careful interpretation of modifications having the same nominal mass.

Rational design and semisynthesis of betulinic acid analogues as potent topoisomerase inhibitors

Chemical transformation studies were conducted on betulinic acid (1), a common plant-derived lupane-type triterpene. Eleven new rationally designed derivatives of 1 (2-5 and 7-13) were synthesized based on docking studies and tested for their topoisomerase I and IIalpha inhibitory activity. Semisynthetic reactions targeted C-3, C-20, and C-28 in 1. Structures of the new compounds were confirmed by spectroscopic methods (1D and 2D NMR and MS). Compound 9, 3-O-[N-(phenylsulfonyl)carbamoyl-17beta-N-(phenylsulfonyl)amide]betulinic acid, showed 1.5-fold the activity of CPT in a topoisomerase I DNA relaxation assay. Four out of 14 betulinic acid analogues (5, 9, 11, and 12) showed 1.5-fold the activity of etoposide in a topoisomerase II assay. The new analogues exhibited better cytotoxic activities against the human colon cancer cells SW948 and HCT-116 and the breast cancer cell line MDA-MB-231 compared to the parent (1). Betulinic acid (1) is a potential scaffold for the design of new topoisomerase I and IIalpha inhibitors.

Identification and characterization of topoisomerase II inhibitory peptides from soy protein hydrolysates

Topoisomerases are targets of several anticancer agents because their inhibition impedes the processes of cell proliferation and differentiation in carcinogenesis. With very limited information available on the inhibitory activities of peptides derived from dietary proteins, the objectives of this study were to employ co-immunoprecipitation to identify inhibitory peptides in soy protein hydrolysates in a single step and to investigate their molecular interactions with topoisomerase II. For this, soy protein isolates were subjected to simulated gastrointestinal digestion with pepsin and pancreatin, and the human topoisomerase II inhibitory peptides were co-immunoprecipitated and identified on a CapLC- Micromass Q-TOF Ultima API system. The inhibitory activity of these peptides from soy isolates toward topoisomerase II was confirmed using three synthetic peptides, FEITPEKNPQ, IETWNPNNKP,and VFDGEL, which have IC 50 values of 2.4, 4.0, and 7.9 mM, respectively. The molecular interactions of these peptides evaluated by molecular docking revealed interaction energies with the topoisomerase II C-terminal domain (CTD) (-186 to -398 kcal/mol) that were smaller than for the ATPase domain (-169 to -357 kcal/mol) and that correlated well with our experimental IC 50 values ( R (2) = 0.99). In conclusion, three peptides released from in vitro gastrointestinal enzyme digestion of soy proteins inhibited human topoisomerase II activity through binding to the active site of the CTD domain.

Differentiation of AmpC beta-lactamase binders vs. decoys using classification kNN QSAR modeling and application of the QSAR classifier to virtual screening

The use of inaccurate scoring functions in docking algorithms may result in the selection of compounds with high predicted binding affinity that nevertheless are known experimentally not to bind to the target receptor. Such falsely predicted binders have been termed 'binding decoys'. We posed a question as to whether true binders and decoys could be distinguished based only on their structural chemical descriptors using approaches commonly used in ligand based drug design. We have applied the k-Nearest Neighbor (kNN) classification QSAR approach to a dataset of compounds characterized as binders or binding decoys of AmpC beta-lactamase. Models were subjected to rigorous internal and external validation as part of our standard workflow and a special QSAR modeling scheme was employed that took into account the imbalanced ratio of inhibitors to non-binders (1:4) in this dataset. 342 predictive models were obtained with correct classification rate (CCR) for both training and test sets as high as 0.90 or higher. The prediction accuracy was as high as 100% (CCR = 1.00) for the external validation set composed of 10 compounds (5 true binders and 5 decoys) selected randomly from the original dataset. For an additional external set of 50 known non-binders, we have achieved the CCR of 0.87 using very conservative model applicability domain threshold. The validated binary kNN QSAR models were further employed for mining the NCGC AmpC screening dataset (69653 compounds). The consensus prediction of 64 compounds identified as screening hits in the AmpC PubChem assay disagreed with their annotation in PubChem but was in agreement with the results of secondary assays. At the same time, 15 compounds were identified as potential binders contrary to their annotation in PubChem. Five of them were tested experimentally and showed inhibitory activities in millimolar range with the highest binding constant K(i) of 135 microM. Our studies suggest that validated QSAR models could complement structure based docking and scoring approaches in identifying promising hits by virtual screening of molecular libraries.

Salvicine functions as novel topoisomerase II poison by binding to ATP pocket

Salvicine, a structurally modified diterpenoid quinone derived from Salvia prionitis, is a nonintercalative topoisomerase II (topo II) poison. The compound possesses potent in vitro and in vivo antitumor activity with a broad spectrum of anti-multidrug resistance activity and is currently in phase II clinical trials. To elucidate the distinct antitumor properties of salvicine and obtain valuable structural information of salvicine-topo II interactions, we characterized the effects of salvicine on human topo IIalpha (htopo IIalpha), including possible binding sites and molecular interactions. The enzymatic assays disclosed that salvicine mainly inhibits the catalytic activity with weak DNA cleavage action, in contrast to the classic topo II poison etoposide (VP16). Molecular modeling studies predicted that salvicine binds to the ATP pocket in the ATPase domain and superimposes on the phosphate and ribose groups. In a surface plasmon resonance binding assay, salvicine exhibited higher affinity for the ATPase domain of htopo IIalpha than ATP and ADP. Competitive inhibition tests demonstrated that ATP competitively and dose-dependently blocked the interactions between salvicine and ATPase domain of htopo IIalpha. The data illustrate that salvicine shares a common binding site with ATP and functions as an ATP competitor. To our knowledge, this is the first report to identify an ATP-binding pocket as the structural binding motif for a nonintercalative eukaryotic topo II poison. These findings collectively support the potential value of an ATP competitor of htopo IIalpha in tumor chemotherapy.

The chemistry of wine polyphenolic C-glycosidic ellagitannins targeting human topoisomerase II

Polyphenolic nonahydroxyterphenoyl-containing C-glycosidic oak ellagitannins are found in wine as a result of the aging of this beverage in oak-made barrels. Once in the slightly acidic wine (pH approximately 3-4), some of these complex natural products such as (-)-vescalagin (1), but not its C-1 epimer (-)-castalagin (2), can capture grape-derived nucleophilic entities such as ethanol, the flavanols catechin (10a) and epicatechin (10b), the anthocyanin oenin (13b), and the thiolic glutathione (16) to furnish condensation products with retention of configuration at the C-1 locus. A computer-aided rationale of this high diastereoselectivity is given. These condensation products can contribute to the modulation of organoleptic properties of the wine, as evidenced by the 23 nm bathochromic shift color absorbance observed with the novel oenin-based anthocyano-ellagitannin (15b). Hydrolysis of 1 under solvolytic conditions furnished another novel compound that we refer to as vescalene (21), in addition to the known (-)-vescalin (18). Of pharmacological importance is the fact that most of these found-in-wine water-soluble ellagitannin derivatives are much more potent than etoposide (VP-16) at inhibiting top2-mediated DNA decatenation in vitro (top2=topoisomerase II)). The known (-)-vescalin (18) and the novel vescalene (21) fully inhibited top2 at 10 microM concentration!

Dynamic Pharmacophore Model Optimization: Identification of Novel HIV-1 Integrase Inhibitors

We extended the previously described dynamic pharmacophore model studies of HIV-1 integrase (IN) by considering more key residues in the active site, including Mg2+. First, we applied a Monte Carlo sampling method to map the complementary features of the IN binding surface. Two types of dynamic pharmacophore models were generated. One considers Mg2+ as part of the IN and therefore as an excluded volume, and the other treats Mg2+ as a positively charged feature, representing a new type of pharmacophore model aimed to identify compounds potentially preventing Mg2+ binding. Second, we validated the models with 385 known active (IC50 < 20 microM) and 235 (IC50 > 100 microM) inactive IN inhibitors. Third, we used the derived models to screen our small molecule database. Twenty-two structurally novel compounds were tested in an in vitro assay specific for IN, and two of them showed IC50 < or = 10 microM for strand transfer reaction.

Virtual Screening and Scaffold Hopping Based on GRID Molecular Interaction Fields

In this study, a set of strategies for structure-based design using GRID molecular interaction fields (MIFs) to derive a pharmacophoric representation of a protein is reported. Thrombin, one of the key enzymes involved in the blood coagulation cascade, was chosen as the model system since abundant published experimental data are available related to both crystal structures and structurally diverse sets of inhibitors. First, a virtual screening methodology was developed either using a pharmacophore representation of the protein based on GRID MIFs or using GRID MIFs from the 3D structure of a set of chosen thrombin inhibitors. The search was done in a 3D multiconformation version of the Available Chemical Directory (ACD) database, which had been spiked with 262 known thrombin inhibitors (multiple conformers available per compound). The model managed to find 80% of the known thrombin inhibitors among the 74,291 conformers in the ACD by only searching 5% of the database; hence, a 15-fold enrichment of the library was achieved. Second, a scaffold hopping methodology was developed using GRID MIFs, giving the scaffold interaction pattern and the shape of the scaffold, together with the distance between the anchor points. The scaffolds reported by Dolle in the Journal of Combinatorial Chemistry summaries (2000 and 2001) and scaffolds built or derived from ligands cocomplexed with the thrombin enzyme were parameterized using a new set of descriptors and saved into a searchable database. The scaffold representation from the database was then compared to a template scaffold (from a thrombin crystal structure), and the thrombin-derived scaffolds included in the database were found among the top solutions. To validate the usefulness of the methodology to replace the template scaffold, the entire molecule was built (scaffold and side chains) and the resulting compounds were docked into the active site of thrombin. The docking solutions showed the same binding pattern as the cocomplexed compound, hence, showing that this method can be a valuable tool for medicinal chemists to select interchangeable core structures (scaffolds) in an easy manner and retaining the binding properties from the original ligand.