Napat Songtawee

Receptor-Based Virtual Screening of EGFR Kinase Inhibitors from the NCI Diversity Database

Epidermal growth factor receptor (EGFR) abnormalities have been associated with several types of human cancer. The crystal structures of its tyrosine kinase domain (EGFR-TK) complexed with small molecule inhibitors revealed the kinase inhibition modes, prompting us to search for novel anti-cancer drugs. A total of 1,990 compounds from the National Cancer Institute (NCI) diversity set with nonredundant structures have been tested to inhibit cancer cell lines with unknown mechanism. Cancer inhibition through EGFR-TK is one of the mechanisms of these compounds. In this work, we performed receptor-based virtual screening against the NCI diversity database. Using two different docking algorithms, AutoDock and Gold, combined with subsequent post-docking analyses, we found eight candidate compounds with high scoring functions that all bind to the ATP-competitive site of the kinase. None of these compounds belongs to the main group of the currently known EGFR-TK inhibitors. Binding mode analyses revealed that the way these compounds complexed with EGFR-TK differs from quinazoline inhibitor binding and the interaction mainly involves hydrophobic interactions. Also, the common kinase-inhibitor (NH---N and CO---HC) hydrogen bonds between the hinge region and the hit compounds are rarely observed. Our results suggest that these molecules could be developed as novel lead compounds in anti-cancer drug design.

Quinoline-based clioquinol and nitroxoline exhibit anticancer activity inducing FoxM1 inhibition in cholangiocarcinoma cells

Purpose Fork head box M1 (FoxM1) is an oncogenic transcription factor frequently elevated in numerous cancers, including cholangiocarcinoma (CCA). A growing body of evidence documents its diverse functions contributing to tumorigenesis and cancer progression. As such, discovery of agents that can target FoxM1 would be valuable for the treatment of CCA. The quinoline-based compounds, namely clioquinol (CQ) and nitroxoline (NQ), represent a new class of anticancer drug. However, their efficacy and underlying mechanisms have not been elucidated in CCA. In this study, anticancer activities and inhibitory effects of CQ and NQ on FoxM1 signaling were explored using CCA cells. Methods The effects of CQ and NQ on cell viability and proliferation were evaluated using the colorimetric 3-(4,5-dimethylthiazol-2yl)-5-(3-carboxymethoxyphenyl)-(4-sulfophenyl)-2H-tetrazolium (MTS assay). Colony formation and cell migration affected by CQ and NQ were investigated using a clonogenic and a wound healing assay, respectively. To demonstrate the agents’ effects on FoxM1 signaling, expression levels of the target genes were quantitatively determined using real-time polymerase chain reaction. Results CQ and NQ significantly inhibited cell survival of HuCCT1 and Huh28 in a dose- and a time-dependent fashion. Further investigations using the rapidly proliferating HuCCT1 cells revealed significant suppression of cell proliferation and colony formation induced by low doses of the compounds. Treatment of CQ and NQ repressed expression of cyclin D1 but enhanced expression of p21. Most importantly, upon CQ and NQ treatment, expression of oncogenic FoxM1 was markedly decreased concomitant with downregulation of various FoxM1’s downstream targets including cdc25b, CENP-B, and survivin. In addition, the compounds distinctly impaired HuCCT1 migration as well as inhibited expression of matrix metalloproteinase (MMP)-2 and MMP-9. Conclusion Collectively, this study reports for the first time the anticancer effects of CQ and NQ against CCA cells, and highlights new insights into the mechanism of actions of the quinoline-based compounds to disrupt FoxM1 signaling.

Computational study of EGFR inhibition: molecular dynamics studies on the active and inactive protein conformations

The structural diversity observed across protein kinases, resulting in subtly different active site cavities, is highly desirable in the pursuit of selective inhibitors, yet it can also be a hindrance from a structure-based design perspective. An important challenge in structure-based design is to better understand the dynamic nature of protein kinases and the underlying reasons for specific conformational preferences in the presence of different inhibitors. To investigate this issue, we performed molecular dynamics simulation on both the active and inactive wild type epidermal growth factor receptor (EGFR) protein with both type-I and type-II inhibitors. Our goal is to better understand the origin of the two distinct EGFR protein conformations, their dynamic differences, and their relative preference for Type-I inhibitors such as gefitinib and Type-II inhibitors such as lapatinib. We discuss the implications of protein dynamics from a structure-based design perspective.

Navigating the chemical space of dipeptidyl peptidase-4 inhibitors

This study represents the first large-scale study on the chemical space of inhibitors of dipeptidyl peptidase-4 (DPP4), which is a potential therapeutic protein target for the treatment of diabetes mellitus. Herein, a large set of 2,937 compounds evaluated for their ability to inhibit DPP4 was compiled from the literature. Molecular descriptors were generated from the geometrically optimized low-energy conformers of these compounds at the semiempirical AM1 level. The origins of DPP4 inhibitory activity were elucidated from computed molecular descriptors that accounted for the unique physicochemical properties inherently present in the active and inactive sets of compounds as defined by their respective half maximal inhibitory concentration values of less than 1 μM and greater than 10 μM, respectively. Decision tree analysis revealed the importance of molecular weight, total energy of a molecule, topological polar surface area, lowest unoccupied molecular orbital, and number of hydrogen-bond donors, which correspond to molecular size, energy, surface polarity, electron acceptors, and hydrogen bond donors, respectively. The prediction model was subjected to rigorous independent testing via three external sets. Scaffold and chemical fragment analysis was also performed on these active and inactive sets of compounds to shed light on the distinguishing features of the functional moieties. Docking of representative active DPP4 inhibitors was also performed to unravel key interacting residues. The results of this study are anticipated to be useful in guiding the rational design of novel and robust DPP4 inhibitors for the treatment of diabetes.

Understanding the molecular basis of EGFR kinase domain/MIG-6 peptide recognition complex using computational analyses

BackgroundEpidermal growth factor receptor (EGFR) signalling plays a major role in biological processes, including cell proliferation, differentiation and survival. Since the over-expression of EGFR causes human cancers, EGFR is an attractive drug target. A tumor suppressor endogenous protein, MIG-6, is known to suppress EGFR over-expression by binding to the C-lobe of EGFR kinase. Thus, this C-lobe of the EGFR kinase is a potential new target for EGFR kinase activity inhibition. In this study, molecular dynamics (MD) simulations and binding free energy calculations were used to investigate the protein-peptide interactions between EGFR kinase and a 27-residue peptide derived from MIG-6_s1 segment (residues 336–362).ResultsThese 27 residues of MIG-6_s1 were modeled from the published MIG-6 X-ray structure. The binding dynamics were detailed by applying the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method to predict the binding free energy. Both van der Waals interactions and non-polar solvation were favorable driving forces for binding process. Six residues of EGFR kinase and eight residues of MIG-6_s1 residues were shown to be responsible for interface binding in which we investigated per residue free energy decomposition and the results from the computational alanine scanning approach. These residues also had higher hydrogen bond occupancies than other residues at the binding interface. The results from the aforementioned calculations reasonably agreed with the previous experimental mutagenesis studies.ConclusionsMolecular dynamics simulations were used to investigate the interactions of MIG-6_s1 to EGFR kinase domain. Our study provides an insight into such interactions that is useful in guiding the design of novel anticancer therapeutics. The information on our modelled peptide interface with EGFR kinase could be a possible candidate for an EGFR dimerization inhibitor.

Insight into HIV-1 reverse transcriptase–aptamer interaction from molecular dynamics simulations

Human immunodeficiency virus-1 reverse transcriptase (HIV-1 RT) is considered to be one of the key targets for antiviral drug therapy. The emergence of the aptamers as potential inhibitors against HIV-1 reverse transcriptase has attracted the attention of the scientific community because these macromolecules can effectively inhibit HIV-1 RT with between micromolar to picomolar concentrations. However, it is not clear how aptamers interact with HIV-1 RT. We have undertaken a molecular dynamics (MD) study in order to gain a keen insight into the conformational dynamics of HIV-1 RT on the formation of a complex with an aptamer or DNA substrate. We have therefore employed three separate models: apo HIV-1 RT, HIV-1 RT with a bound RNA aptamer, and HIV-1 RT with a bound DNA substrate. The results show that HIV-1 RT complex with an aptamer was more stable than that with DNA substrate. It was found that the aptamer interacted with HIV-1 RT in a fingers-and-thumb-closed conformation, at the bound at the nucleic acid substrate binding site. We identified key residues within the HIV-1 RT-aptamer complex in order to help design, develop, and test a new aptamer based on therapies in the future.

Computational analysis of binding between malarial dihydrofolate reductases and anti-folates

BackgroundPlasmodium falciparum readily develops resistance to the anti-folates pyrimethamine and proguanil via a characteristic set of mutations in the dihydrofolate reductase (PfDHFR) gene that leads to reduced competitive drug binding at the enzymes active site. Analogous mutations can be found in the DHFR gene in isolates of Plasmodium vivax (PvDHFR) although anti-folates have not been widely used for the treatment of this infection. Here the interactions between DHFR inhibitors and modelled structures of the DHFR enzymes of Plasmodium malariae (PmDHFR) and Plasmodium ovale (PoDHFR) are described, along with an investigation of the effect of recently reported mutations within PmDHFR.MethodsDHFR models for PmDHFR and PoDHFR were constructed using the solved PfDHFR-TS and PvDHFR structures respectively as templates. The modelled structures were docked with three DHFR inhibitors as ligands and more detailed interactions were explored via simulation of molecular dynamics.ResultsHighly accurate models were obtained containing sets of residues that mediate ligand binding which are highly comparable to those mediating binding in known crystal structures. Within this set, there were differences in the relative contribution of individual residues to inhibitor binding. Modelling of PmDHFR mutant sequences revealed that PmDHFR I170M was associated with a significant reduction in binding energy to all DHFR inhibitors studied, while the other predicted resistance mutations had lesser or no effects on ligand binding.ConclusionsBinding of DHFR inhibitors to the active sites of all four Plasmodium enzymes is broadly similar, being determined by an analogous set of seven residues. PmDHFR mutations found in field isolates influenced inhibitor interactions to a varying extent. In the case of the isolated I170M mutation, the loss of interaction with pyrimethamine suggests that DHFR-inhibitor interactions in P. malariae are different to those seen for DHFRs from P. falciparum and P. vivax.

Structural and biochemical characterization of two heme binding sites on α1-microglobulin using site directed mutagenesis and molecular simulation.

BACKGROUND α1-Microglobulin (A1M) is a reductase and radical scavenger involved in physiological protection against oxidative damage. These functions were previously shown to be dependent upon cysteinyl-, C34, and lysyl side-chains, K(92, 118,130). A1M binds heme and the crystal structure suggests that C34 and H123 participate in a heme binding site. We have investigated the involvement of these five residues in the interactions with heme. METHODS Four A1M-variants were expressed: with cysteine to serine substitution in position 34, lysine to threonine substitutions in positions (92, 118, 130), histidine to serine substitution in position 123 and a wt without mutations. Heme binding was investigated by tryptophan fluorescence quenching, UV-Vis spectrophotometry, circular dichroism, SPR, electrophoretic migration shift, gel filtration, catalase-like activity and molecular simulation. RESULTS All A1M-variants bound to heme. Mutations in C34, H123 or K(92, 118, 130) resulted in significant absorbance changes, CD spectral changes, and catalase-like activity, suggesting involvement of these side-groups in coordination of the heme-iron. Molecular simulation support a model with two heme-binding sites in A1M involving the mutated residues. Binding of the first heme induces allosteric stabilization of the structure predisposing for a better fit of the second heme. CONCLUSIONS The results suggest that one heme-binding site is located in the lipocalin pocket and a second binding site between loops 1 and 4. Reactions with the hemes involve the side-groups of C34, K(92, 118, 130) and H123. GENERAL SIGNIFICANCE The model provides a structural basis for the functional activities of A1M: heme binding activity of A1M.

Molecular characterization of Galectin-8 from Nile tilapia (Oreochromis niloticus Linn.) and its response to bacterial infection.

Galectins belong to the family of galactoside-binding proteins and play a major role in the immune and inflammatory responses of vertebrates and invertebrates. The galectin family is divided into three subtypes based on molecular structure; prototypes, chimera types, and tandem-repeated types. We isolated and characterized the cDNA of galectin-8 (OnGal-8) in Nile tilapia (Oreochromis niloticus). OnGal-8 consisted of a 966 bp open reading frame (ORF) that encoded a 321 amino acid protein (43.47 kDa). Homology and phylogenetic tree analysis suggested the protein was clustered with galectin-8s from other animal species and shared at least 56.8% identity with salmon galectin-8. Structurally, the amino acid sequence included two distinct N- and C- terminus carbohydrate recognition domains (CRDs) of 135 and 133 amino acids, respectively, that were connected by a 39 amino acid polypeptide linker. The N- and C-CRDs contained two conserved WG-E-I and WG-E-T motifs, suggesting they have an important role in mediating the specific interactions between OnGal-8 and saccharide moieties such as β-galactoside. The structure of OnGal-8 was characterized by a two-fold symmetric pattern of 10-and 12-stranded antiparallel ß-sheets of both N- and C-CRDs, and the peptide linker presumably formed a random coil similar to the characteristic tandem-repeat type galectin. The expression of OnGal-8 in healthy fish was highest in the skin, intestine, and brain. Experimental challenge of Nile tilapia with S. agalactiae resulted in significant up-regulation of OnGal-8in the spleen after 5 d. Our results suggest that OnGal-8 is involved in the immune response to bacterial infection.

Molecular dynamics of the asymmetric dimers of EGFR: simulations on the active and inactive conformations of the kinase domain.

Abnormal activation of EGFR is associated with human cancer, and thus it is a key target for inhibition in cancer therapy. There is evidence suggesting that the activation mechanism of EGFR is based upon the formation of the asymmetric dimer of the kinase domains. Here, we performed MD simulations on the asymmetric dimer for both active and inactive conformations of EGFR kinase domain to investigate flexibility and intrinsic motions of the proteins. Simulations of the active conformation showed that the formation of the asymmetric dimer changes the dynamics of EGFR kinase domain by suppressing fluctuation of the protein and altering the direction of motion of the protein. In contrast, the asymmetric dimerization of the inactive conformation does not alter the overall fluctuation of the kinase domain and does not initiate destabilizing of the inactive structure. We also investigated the intermolecular interactions in the EGFR asymmetric dimers and found that in the active conformation the interactions are dominated by loop-loop contacts rather than those from the helix-helix interactions. In contrast, helix-helix interaction seemed to be more significant for the inactive kinase structure. This work helps us to better understand the conformational flexibility and dynamics of the EGFR kinase domain, as well as provides information that may be useful to develop newer classes of inhibitors that can block allosteric sites rather than the more traditional catalytic site.