Docking and simulation study of the enantiospecificity of chiral epoxidation reactions catalyzed by chloroperoxidase

Chloroperoxidase (CPO), a 298-residue glycosylated protein from the fungus Caldariomyces fumago, is probably the most versatile heme enzyme yet discovered. Interest in CPO as a catalyst is based on its power to produce enantiomerically enriched products. Recent research has focused its attention on the ability of CPO to epoxidize alkenes in high regioselectivity and enantioselectivity as an efficient and environmentally benign alternative to traditional synthetic routes. There has been little work on the nature of ligand binding, which probably controls the regio- and enantiospecifity of CPO. Consequently it is here that we focus our work. We report docking calculations and computer simulations aimed at predicting the enantiospecificity of CPO-catalyzed epoxidation of three model substrates. On the basis of this work candidate mutations to improve the efficiency of CPO are predicted. In order to accomplish these aims, a simulated annealing and molecular dynamics protocol is developed to sample potentially reactive substrate/CPO complexes.

Molecular docking and geographical information systems as tools to assess the potential impact of veterinary medicines on non-target organisms and the environment

Veterinary medicines (VMs) from agricultural industry can enter the environment in a number of ways. This includes direct exposure through aquaculture, accidental spillage and disposal, and indirect entry by leaching from manure or runoff after treatment. Many compounds used in animal treatments have ecotoxic properties that may have chronic or sometimes lethal effects when they come into contact with non-target organisms. VMs enter the environment in mixtures, potentially having additive effects. Traditional ecotoxicology tests are used to determine the lethal and sometimes reproductive effects on freshwater and terrestrial organisms. However, organisms used in ecotoxicology tests can be unrepresentative of the populations that are likely to be exposed to the compound in the environment. Most often the tests are on single compound toxicity but mixture effects may be significant and should be included in ecotoxicology testing. This work investigates the use, measured environmental concentrations (MECs) and potential impact of sea lice treatments on salmon farms in Scotland. Alternative methods for ecotoxicology testing including mixture toxicity, and the use of in silico techniques to predict the chronic impact of VMs on different species of aquatic organisms were also investigated. The Scottish Environmental Protection Agency (SEPA) provided information on the use of five sea lice treatments from 2008-2011 on Scottish salmon farms. This information was combined with the recently available data on sediment MECs for the years 2009-2012 provided by SEPA using ArcGIS 10.1. In depth analysis of this data showed that from a total of 55 sites, 30 sites had a MEC higher than the maximum allowable concentration (MAC) as set out by SEPA for emamectin benzoate and 7 sites had a higher MEC than MAC for teflubenzuron. A number of sites that were up to 16 km away from the nearest salmon farm reported as using either emamectin benzoate or teflubenzuron measured positive for the two treatments. There was no relationship between current direction and the distribution of the sea lice treatments, nor was there any evidence for alternative sources of the compounds e.g. land treatments. The sites that had MECs higher than the MAC could pose a risk to non-target organisms and disrupt the species dynamics of the area. There was evidence that some marine protected sites might be at risk of exposure to these compounds. To complement this work, effects on acute mixture toxicity of the 5 sea lice treatments, plus one major metabolite 3-phenoxybenzoic acid (3PBA), were measured using an assay using the bioluminescent bacteria Aliivibrio fischeri. When exposed to the 5 sea lice treatments and 3PBA A. fischeri showed a response to 3PBA, emamectin benzoate and azamethiphos as well as combinations of the three. In order to establish any additive effect of the sea lice treatments, the efficacy of two mixture prediction equations, concentration addition (CA) and independent action ii(IA) were tested using the results from single compound dose response curves. In this instance IA was the more effective prediction method with a linear regression confidence interval of 82.6% compared with 22.6% of CA. In silico molecular docking was carried out to predict the chronic effects of 15 VMs (including the five used as sea lice control). Molecular docking has been proposed as an alternative screening method for the chronic effects of large animal treatments on non-target organisms. Oestrogen receptor alpha (ERα) of 7 non-target bony fish and the African clawed frog Xenopus laevis were modelled using SwissModel. These models were then ‘docked’ to oestradiol, the synthetic oestrogen ethinylestradiol, two known xenoestrogens dichlorodiphenyltrichloroethane (DDT) and bisphenol A (BPA), the antioestrogen breast cancer treatment tamoxifen and 15 VMs using Auto Dock 4. Based on the results of this work, four VMs were identified as being possible xenoestrogens or anti-oestrogens; these were cypermethrin, deltamethrin, fenbendazole and teflubenzuron. Further investigation, using in vitro assays, into these four VMs has been suggested as future work. A modified recombinant yeast oestrogen screen (YES) was attempted using the cDNA of the ERα of the zebrafish Danio rerio and the rainbow trout Oncorhynchus mykiss. Due to time and difficulties in cloning protocols this work was unable to be completed. Use of such in vitro assays would allow for further investigation of the highlighted VMs into their oestrogenic potential. In conclusion, VMs used as sea lice treatments, such as teflubenzuron and emamectin benzoate may be more persistent and have a wider range in the environment than previously thought. Mixtures of sea lice treatments have been found to persist together in the environment, and effects of these mixtures on the bacteria A. fischeri can be predicted using the IA equation. Finally, molecular docking may be a suitable tool to predict chronic endocrine disrupting effects and identify varying degrees of impact on the ERα of nine species of aquatic organisms.

Docking of HIV-1 Vpr to the nuclear envelope is mediated by the interaction with the nucleoporin hCG1

The HIV-1 genome contains several genes coding for auxiliary proteins, including the small Vpr protein. Vpr affects the integrity of the nuclear envelope and participates in the nuclear translocation of the preintegration complex containing the viral DNA. Here, we show by photobleaching experiments performed on living cells expressing a Vpr-green fluorescent protein fusion that the protein shuttles between the nucleus and the cytoplasm, but a significant fraction is concentrated at the nuclear envelope, supporting the hypothesis that Vpr interacts with components of the nuclear pore complex. An interaction between HIV-1 Vpr and the human nucleoporin CG1 (hCG1) was revealed in the yeast two-hybrid system, and then confirmed both in vitro and in transfected cells. This interaction does not involve the FG repeat domain of hCG1 but rather the N-terminal region of the protein. Using a nuclear import assay based on digitonin-permeabilized cells, we demonstrate that hCG1 participates in the docking of Vpr at the nuclear envelope. This association of Vpr with a component of the nuclear pore complex may contribute to the disruption of the nuclear envelope and to the nuclear import of the viral DNA.

Data of the molecular dynamics simulations of mutations in the human connexin46 docking interface

The structure of hCx26 derived from the X-ray analysis was used to generate a homology model for hCx46. Interacting connexin molecules were used as starting model for the molecular dynamics (MD) simulation using NAMD and allowed us to predict the dynamic behavior of hCx46wt and the cataract related mutant hCx46N188T as well as two artificial mutants hCx46N188Q and hCx46N188D. Within the 50 ns simulation time the docked complex composed of the mutants dissociate while hCx46wt remains stable. The data indicates that one hCx46 molecule forms 5-7 hydrogen bonds (HBs) with the counterpart connexin of the opposing connexon. These HBs appear essential for a stable docking of the connexons as shown by the simulation of an entire gap junction channel and were lost for all the tested mutants. The data described here are related to the research article entitled "The cataract related mutation N188T in human connexin46 (hCx46) revealed a critical role for residue N188 in the docking process of gap junction channels" (Schadzek et al., 2015) [1].

Low molecular weight compounds from mushrooms as potential Bcl-2 inhibitors: docking and virtual screening studies

Mushrooms have the ability to promote apoptosis in tumor cell lines, but the mechanism of action is not quite well understood. Inhibition of the interaction between Bcl-2 and pro-apoptotic proteins could be an important step that leads to apoptosis. Therefore, the discovery of compounds with the ability to inhibit Bcl-2 is an ongoing research topic in drug discovery. In this study, we started by analyzing Bcl-2 experimental structures that are currently available in Protein Data Bank database. After analysis of the more relevant Bcl-2 structures, 4 were finally selected. An analysis of the best docking methodology was then performed using a cross-docking and re-docking approach while testing 2 docking softwares: AutoDock 4 and AutoDock Vina. Autodock4 provided the best docking results and was selected to perform a virtual screening study applied to a dataset of 40 Low Molecular Weight (LMW) compounds present in mushrooms, using the selected Bcl-2 structures as target. Results suggest that steroid are the more promising family, among the analyzed compounds, and may have the ability to interact with Bcl-2 and this way promoting tumor apoptosis. The steroids that presented lowest estimated binding energy (ΔG) were: Ganodermanondiol, Cerevisterol, Ganoderic Acid X and Lucidenic Lactone; with estimated ΔG values between -8,45 and -8,23 Kcal/mol. A detailed analysis of the docked conformation of these 4 top ranked LMW compounds was also performed and illustrates a plausible interaction between the 4 top raked steroids and Bcl-2, thus substantiating the accuracy of the predicted docked poses. Therefore, tumoral apoptosis promoted by mushroom might be related to Bcl-2 inhibition mediated by steroid family of compounds.

Investigating mushroom LMW compounds as potential Bcl-2 inhibitors: docking studies using AutoDock4

The B cell CLL/lymphoma-2 (Bcl-2) family is functionally classified as either anti-apoptotic or pro-apoptotic, and the regulation of its interactions dictates survival or commitment to apoptosis. Bcl-2 family is also implicated in a wide range of diseases. In some types of cancers, including lymphomas and epithelial cancers, protein overexpression of anti-apoptotic Bcl-2 family, such as the Bcl-2 protein is indicative of cancer in an advanced stage, with a poor prognosis and resistant to chemotherapy [1]. Several reports indicate that mushrooms have the ability to promote apoptosis in tumour cell lines, but the mechanism of action is not fully understood. Inhibition of the interaction between Bcl-2 (anti-apoptotic protein) and proapoptotic proteins could be an important step in the mechanism of mushroom induced apoptosis. Therefore, the discovery of compounds with the capacity to inhibit Bcl-2 is an ongoing research topic on cancer therapy. In this work, docking studies were performed using a dataset of 40 low molecular weight (LMW) compounds present in mushrooms. The docking software AutoDock 4 was used and docking studies were performed using 5 selected Bcl-2 crystal structures as targets. Compounds with the lowest predicted binding energy (predΔG) are expected to be the more potent inhibitors. Among the tested compounds, steroids presented the lowest predΔG with several exhibiting values below -9 kcal/mol. The results are corroborated by several reports that state that steroids induce apoptosis in several tumor cells. It is thus feasible that they might act by preventing Bcl-2 from forming complexes with the respective proapoptotic protein interaction partners, namely Bak, Bax, and Bim. Moreover, previous studies on our research group demonstrated that 48 h treatment of MCF-7 cells (breast carcinoma) with Suillus collinitus methanolic extract caused a decrease in Bcl-2, highlighting the antitumor potential of this mushroom species [2]. In conclusion, the process of apoptosis promoted by mushroom extracts may be related to the inhibition of Bcl-2 by the steroid derivatives herein studied. However, further studies are needed to confirm this hypothesis.

A New Multi-Objective Approach for Molecular Docking Based on RMSD and Binding Energy

Ligand-protein docking is an optimization problem based on predicting the position of a ligand with the lowest binding energy in the active site of the receptor. Molecular docking problems are traditionally tackled with single-objective, as well as with multi-objective approaches, to minimize the binding energy. In this paper, we propose a novel multi-objective formulation that considers: the Root Mean Square Deviation (RMSD) difference in the coordinates of ligands and the binding (intermolecular) energy, as two objectives to evaluate the quality of the ligand-protein interactions. To determine the kind of Pareto front approximations that can be obtained, we have selected a set of representative multi-objective algorithms such as NSGA-II, SMPSO, GDE3, and MOEA/D. Their performances have been assessed by applying two main quality indicators intended to measure convergence and diversity of the fronts. In addition, a comparison with LGA, a reference single-objective evolutionary algorithm for molecular docking (AutoDock) is carried out. In general, SMPSO shows the best overall results in terms of energy and RMSD (value lower than 2A for successful docking results). This new multi-objective approach shows an improvement over the ligand-protein docking predictions that could be promising in in silico docking studies to select new anticancer compounds for therapeutic targets that are multidrug resistant.

New cholinesterase inhibitors for Alzheimer's disease: Structure Activity Studies (SARs) and molecular docking of isoquinolone and azepanone derivatives

A library of isoquinolinone and azepanone derivatives were screened for both acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE) activity. The strategy adopted included (a) in vitro biological assays, against eel AChE (EeAChE) and equine serum BuChE (EqBuChE) in order to determine the compounds IC50 and their dose-response activity, consolidated by (b) molecular docking studies to evaluate the docking poses and interatomic interactions in the case of the hit compounds, validated by STD-NMR studies. Compound (1f) was identified as one of these hits with an IC50 of 89.5 mu M for EeAChE and 153.8 mu M for EqBuChE, (2a) was identified as a second hit with an IC50 of 108.4 mu M (EeAChE) and 277.8 mu M (EqBuChE). In order to gain insights into the binding mode and principle active site interactions of these molecules, (R)-(1f) along with 3 other analogues (also as the R-enantiomer) were docked into both RhAChE and hBuChE models. Galantamine was used as the benchmark. The docking study was validated by performing an STD-NMR study of (1f) with EeAChE using galantamine as the benchmark.

Insights into (S)-rivastigmine inhibition of butyrylcholinesterase (BuChE): Molecular docking and saturation transfer difference NMR (STD-NMR)

Rivastigmine is a very important drug prescribed for the treatment of Alzheimer's disease (AD) symptoms. It is a dual inhibitor, in that it inhibits both acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE). For our screening program on the discovery of new rivastigmine analogue hits for human butyrylcholinesterase (hBuChE) inhibition, we investigated the interaction of this inhibitor with BuChE using the complimentary approach of the biophysical method, saturation transfer difference (STD)-NMR and molecular docking. This allowed us to obtain essential information on the key binding interactions between the inhibitor and the enzyme to be used for screening of hit compounds. The main conclusions obtained from this integrated study was that the most dominant interactions were (a) H-bonding between the carbamate carbonyl of the inhibitor and the NH group of the imidazole unit of H434, (b) stacking of the aromatic unit of the inhibitor and the W82 aromatic unit in the choline binding pocket via pi-pi interactions and (c) possible CH/pi interactions between the benzylic methyl group and the N-methyl groups of the inhibitor and W82 of the enzyme.

Synthesis And Antiproliferative Activity Of 8-hydroxyquinoline Derivatives Containing A 1,2,3-triazole Moiety.

Twelve novel 8-hydroxyquinoline derivatives were synthesized with good yields by performing copper-catalyzed Huisgen 1,3-dipolar cycloaddition (click reaction) between an 8-O-alkylated-quinoline containing a terminal alkyne and various aromatic or protected sugar azides. These compounds were evaluated in vitro for their antiproliferative activity on various cancer cell types. Protected sugar derivative 16 was the most active compound in the series, exhibiting potent antiproliferative activity and high selectivity toward ovarian cancer cells (OVCAR-03, GI50 < 0.25 μg mL(-1)); this derivative was more active than the reference drug doxorubicin (OVCAR-03, GI50 = 0.43 μg mL(-1)). In structure-activity relationship (SAR) studies, the physico-chemical parameters of the compounds were evaluated and docking calculations were performed for the α-glucosidase active site to predict the possible mechanism of action of this series of compounds.

Kinase Inhibitor Profile For Human Nek1, Nek6, And Nek7 And Analysis Of The Structural Basis For Inhibitor Specificity.

Human Neks are a conserved protein kinase family related to cell cycle progression and cell division and are considered potential drug targets for the treatment of cancer and other pathologies. We screened the activation loop mutant kinases hNek1 and hNek2, wild-type hNek7, and five hNek6 variants in different activation/phosphorylation statesand compared them against 85 compounds using thermal shift denaturation. We identified three compounds with significant Tm shifts: JNK Inhibitor II for hNek1(Δ262-1258)-(T162A), Isogranulatimide for hNek6(S206A), andGSK-3 Inhibitor XIII for hNek7wt. Each one of these compounds was also validated by reducing the kinases activity by at least 25%. The binding sites for these compounds were identified by in silico docking at the ATP-binding site of the respective hNeks. Potential inhibitors were first screened by thermal shift assays, had their efficiency tested by a kinase assay, and were finally analyzed by molecular docking. Our findings corroborate the idea of ATP-competitive inhibition for hNek1 and hNek6 and suggest a novel non-competitive inhibition for hNek7 in regard to GSK-3 Inhibitor XIII. Our results demonstrate that our approach is useful for finding promising general and specific hNekscandidate inhibitors, which may also function as scaffolds to design more potent and selective inhibitors.

Mechanism of Heparin Acceleration of Tissue Inhibitor of Metalloproteases-1 (TIMP-1) Degradation by the Human Neutrophil Elastase

Heparin has been shown to regulate human neutrophil elastase (HNE) activity. We have assessed the regulatory effect of heparin on Tissue Inhibitor of Metalloproteases-1 [TIMP-1] hydrolysis by HNE employing the recombinant form of TIMP-1 and correlated FRET-peptides comprising the TIMP-1 cleavage site. Heparin accelerates 2.5-fold TIMP-1 hydrolysis by HNE. The kinetic parameters of this reaction were monitored with the aid of a FRET-peptide substrate that mimics the TIMP-1 cleavage site in pre-steady-state conditionsby using a stopped-flow fluorescence system. The hydrolysis of the FRET-peptide substrate by HNE exhibits a pre-steady-state burst phase followed by a linear, steady-state pseudo-first-order reaction. The HNE acylation step (k(2)=21 +/- 1 s(-1)) was much higher than the HNE deacylation step (k(3)=0.57 +/- 0.05 s(-1)). The presence of heparin induces a dramatic effect in the pre-steady-state behavior of HNE. Heparin induces transient lag phase kinetics in HNE cleavage of the FRET-peptide substrate. The pre-steady-state analysis revealed that heparin affects all steps of the reaction through enhancing the ES complex concentration, increasing k(1) 2.4-fold and reducing k(-1) 3.1-fold. Heparin also promotes a 7.8-fold decrease in the k(2) value, whereas the k(3) value in the presence of heparin was increased 58-fold. These results clearly show that heparin binding accelerates deacylation and slows down acylation. Heparin shifts the HNE pH activity profile to the right, allowing HNE to be active at alkaline pH. Molecular docking and kinetic analysis suggest that heparin induces conformational changes in HNE structure. Here, we are showing for the first time that heparin is able to accelerate the hydrolysis of TIMP-1 by HNE. The degradation of TIMP-1is associated to important physiopathological states involving excessive activation of MMPs.

Sulfonyl-hydrazones of cyclic imides derivatives as potent inhibitors of the Mycobacterium tuberculosis protein tyrosine phosphatase B (PtpB)

Searching lead compounds for new antituberculosis drugs, the activity of synthetic sulfonamides and sulfonyl-hydrazones were assayed for their potential inhibitory activity towards a protein tyrosine phosphatase from Mycobacterium tuberculosis - PtpB. Four sulfonyl-hydrazones N-phenylmaleimide derivatives were active (compounds 14, 15, 19 and 21), and the inhibition of PtpB was found to be competitive with respect to the substrate p-nitrophenyl phosphate. Structure-based molecular docking simulations were performed and indicated that the new inhibitor candidates showed similar binding modes, filling the hydrophobic pocket of the protein by the establishment of van der Waals contacts, thereby contributing significantly to the complex stability.

New antitumoral agents I: In vitro anticancer activity and in vivo acute toxicity of synthetic 1,5-bis(4-hydroxy-3-methoxyphenyl)-1,4-pentadien-3-one and derivatives

This paper describes a new method for the preparation of 1,5-bis(4-hydroxy-3-methoxyphenyl)-1,4-pentadien-3-one 1 and its derivatives 2-5. This set of synthetic compounds exhibited high antitumoral activities regarding in vitro screening against several human tumor cell lines as lung carcinoma NCI-460, melanoma UACC-62, breast MCF-7, colon HT-29, renal 786-O, ovarian OVCAR-03 and ovarian expressing the resistance phenotype for adriamycin NCI-ADR/ RES, prostate PC-3, and leukemia K-562. Compounds were also tested against murine tumor cell line B16F10 melanoma and lymphocytic leukemia L1210 as well as to their effect toward normal macrophages. Specific activity against colon cancer cells HT-29 was observed for all tested compounds and suggests further studies with models of colon cancer. Compounds 1, 2, and 4 showed significant cytotoxic activity with IC(50) values <= 2.3 mu M for all human cancer cell lines. Intraperitoneal acute administration of compound 1 and 2 showed very low toxicity rate. (C) 2010 Elsevier Ltd. All rights reserved.