Protein-ligand docking study: diterpenes from Juniperus brevifolia as anticancer and antimicrobial agentes

REDCAT: Natural Products and related Redox Catalysts: Basic Research and Applications in Medicine and Agriculture, Aveiro, 25-27 Novembro de 2012.

Predicting the Pose of β-Casomorphin-5 and 7 in the Opioid Receptors

The opioid receptors consist of three main subtypes; μ, δ, and κ. Previous binding studies have shown that fragments of the milk protein, β-casein, known as β-casomorphins are agonists of these receptors which are selective for the μ receptor subtype. Using the crystal structures of these three receptors, computational molecular docking studies were done using the software GOLD to determine the conformation of β-casomorphin-5 and 7 when they bind to these three opioid receptors. GOLD was able to discriminate among the three receptors when docking the rigid ligands co-crystalized with the receptors. However, GOLD could not discriminate among the three receptors for either of the highly flexible β-casomorphins. A per amino acid scoring method was developed to overcome this problem. This method was used to predict the conformation of both β-casomorphin-5 and 7 in the μ receptor and determine that the two amino acid residues, Lys303 and Trp318 of the μ receptor are responsible for discriminating among the three receptor subtypes for binding of the β-casomorphin-5 and 7.

MPO Inhibitors Selected by Virtual Screening

The hemeprotein myeloperoxidase (MPO) participates in innate immune defense through its ability to generate potent microbicidal oxidants. However, these oxidants are also key mediators of the tissue damage associated with many inflammatory diseases. Thus, there is considerable interest in developing therapeutically useful MPO inhibitors. Here, we used structure-based drug design (SBDD) and ligand-based drug design (LBDD) to select for potentially new and selective MPO inhibitors. A pharmacophore model was developed based on the crystal structure of human MPO in complex with salicylhydroxamic acid (SHA), a known inhibitor of the enzyme. The pharmacophore model was used to screen the ZINC database for potential ligands, which were further filtered on the basis of their physical-chemical properties and docking score. The filtered compounds were visually inspected, and nine were purchased for experimental studies. Surprisingly, almost all of the selected compounds belonged to the aromatic hydrazide class, which had been previously described as MPO inhibitors. The compounds selected by virtual screening were shown to inhibit the chlorinating activity of MPO; the top four compounds displayed IC(50) values ranging from 1.0 to 2.8 mM. MPO inactivation by the most effective compound was shown to be irreversible. Overall, our results show that SBDD and LBDD may be useful for the rational development of new MPO inhibitors.

FERM domain interaction with myosin negatively regulates FAK in cardiomyocyte hypertrophy

Focal adhesion kinase (FAK) regulates cellular processes that affect several aspects of development and disease. The FAK N-terminal FERM (4.1 protein-ezrin-radixin-moesin homology) domain, a compact clover-leaf structure, binds partner proteins and mediates intramolecular regulatory interactions. Combined chemical cross-linking coupled to MS, small-angle X-ray scattering, computational docking and mutational analyses showed that the FAK FERM domain has a molecular cleft (similar to 998 angstrom(2)) that interacts with sarcomeric myosin, resulting in FAK inhibition. Accordingly, mutations in a unique short amino acid sequence of the FERM myosin cleft, FP-1, impaired the interaction with myosin and enhanced FAK activity in cardiomyocytes. An FP-1 decoy peptide selectively inhibited myosin interaction and increased FAK activity, promoting cardiomyocyte hypertrophy through activation of the AKT-mammalian target of rapamycin pathway. Our findings uncover an inhibitory interaction between the FAK FERM domain and sarcomeric myosin that presents potential opportunities to modulate the cardiac hypertrophic response through changes in FAK activity.

Optimization of anti-Trypanosoma cruzi oxadiazoles leads to identification of compounds with efficacy in infected mice

We recently showed that oxadiazoles have anti-Trypanosoma cruzi activity at micromolar concentrations. These compounds are easy to synthesize and show a number of clear and interpretable structure-activity relationships (SAR), features that make them attractive to pursue potency enhancement. We present here the structural design, synthesis, and anti-T. cruzi evaluation of new oxadiazoles denoted 5a-h and 6a-h. The design of these compounds was based on a previous model of computational docking of oxadiazoles on the T. cruzi protease cruzain. We tested the ability of these compounds to inhibit catalytic activity of cruzain, but we found no correlation between the enzyme inhibition and the antiparasitic activity of the compounds. However, we found reliable SAR data when we tested these compounds against the whole parasite. While none of these oxadiazoles showed toxicity for mammalian cells, oxadiazoles 6c (fluorine), 6d (chlorine), and 6e (bromine) reduced epimastigote proliferation and were cidal for trypomastigotes of T. cruzi Y strain. Oxadiazoles 6c and 6d have IC50 of 9.5 +/- 2.8 and 3.5 +/- 1.8 mu M for trypomastigotes, while Benznidazole, which is the currently used drug for Chagas disease treatment, showed an IC50 of 11.3 +/- 2.8 mu M. Compounds 6c and 6d impair trypomastigote development and invasion in macrophages, and also induce ultrastructural alterations in trypomastigotes. Finally, compound 6d given orally at 50 mg/kg substantially reduces the parasitemia in T. cruzi-infected BALB/c mice. Our drug design resulted in potency enhancement of oxadiazoles as anti-Chagas disease agents, and culminated with the identification of oxadiazole 6d, a trypanosomicidal compound in an animal model of infection. (C) 2012 Elsevier Ltd. All rights reserved.

Epigallocatechin gallate (EGCG) affects glucose metabolism and enhances fitness and life span in Drosophila melanogaster

In this study, we tested whether a standardized epigallocatechin-3-gallate (EGCG) rich green tea extract (comprising > 90% EGCG) affects fitness and lifespan as well as parameters of glucose metabolism and energy homeostasis in the fruit fly, Drosophila melanogaster. Following the application of the green tea extract a significant increase in the mean lifespan (+ 3.3 days) and the 50% survival (+ 4.3 days) as well as improved fitness was detected. These effects went along an increased expression of Spargel, the homolog of mammalian PGC1α, which has been reported to affect lifespan in flies. Intriguingly, in flies, treatment with the green tea extract decreased glucose concentrations, which were accompanied by an inhibition of α-amylase and α-glucosidase activity. Computational docking analysis proved the potential of EGCG to dock into the substrate binding pocket of α-amylase and to a greater extent into α-glucosidase. Furthermore, we demonstrate that EGCG downregulates insulin-like peptide 5 and phosphoenolpyruvate carboxykinase, major regulators of glucose metabolism, as well as the Drosophila homolog of leptin, unpaired 2. We propose that a decrease in glucose metabolism in connection with an upregulated expression of Spargel contribute to the better fitness and the extended lifespan in EGCG-treated flies.

New cytochromes P450: Tale of two variants

Cytochrome P450s, a superfamily of heme enzymes found in most living organisms. They are responsible for metabolism of many therapeutic drugs, industrial pollutants, carcinogens, and additives to foodstuffs, as well as some endogenous compounds including fatty acids and steroids. First pass drug metabolism studies represent mainly liver and small intestine elimination, and are viewed as the standard to predict therapeutic outcome. However, drug plasma levels determined after administration do not always correlate with therapeutic efficacy of the drug. Therefore, a possible explanation may come by understanding drug metabolism in extrahepatic tissues and/or at the site of drug action. Identification and characterization of novel tissue specific isoforms of P450 generated by alternative splicing of known P450 genes or as yet unidentified genes is essential to predict pharmacological outcome of drugs or the fate of a carcinogen that act at sites remote from liver. ^ Using RT-PCR, brain-specific cytochrome P450s were detected in samples of human autopsy brain. So far, we have identified two human brain variants including P450 2D7 and P450 1A1. We have shown the presence of the P450 1A1 brain specific splice variant in African Americans, Caucasians and Indians albeit different patterns of liver to brain variant ratio were seen distributed throughout each population. Interestingly, the splice variant was detected only in the brain but not in any other tissues from the same individual. Homology modeling was used to compare the variant 3D structure to the liver form structure and differences in the substrate access channels and substrate binding sites were noticed. Automated computational docking was used to predict the metabolic fate of the potent carcinogenic substrate, benzo[a]pyrene. P450 1A1 brain variant showed no binding orientations that could produce the active metabolite, whereas P450 1A1 liver form did reveal orientations capable of generating active carcinogenic product. In vitro P32 labeling studies verified the docking predictions. Therefore, the data support the hypothesis that P450 brain splice variants mediate the metabolism of xenobiotics by mechanisms distinct from the well-studied liver counterparts. ^

Structure-activity relationships for a new family of sulfonylurea herbicides

Acetohydroxyacid synthase (AHAS; EC 2.2.1.6) catalyzes the first common step in branched-chain amino acid biosynthesis. The enzyme is inhibited by several chemical classes of compounds and this inhibition is the basis of action of the sulfonylurea and imidazolinone herbicides. The commercial sulfonylureas contain a pyrimidine or a triazine ring that is substituted at both meta positions, thus obeying the initial rules proposed by Levitt. Here we assess the activity of 69 monosubstituted sulfonylurea analogs and related compounds as inhibitors of pure recombinant Arabidopsis thaliana AHAS and show that disubstitution is not absolutely essential as exemplified by our novel herbicide, monosulfuron (2-nitro-N-(4'-methyl-pyrimidin-2'-yl) phenyl-sulfonylurea), which has a pyrimidine ring with a single meta substituent. A subset of these compounds was tested for herbicidal activity and it was shown that their effect in vivo correlates well with their potency in vitro as AHAS inhibitors. Three-dimensional quantitative structure-activity relationships were developed using comparative molecular field analysis and comparative molecular similarity indices analysis. For the latter, the best result was obtained when steric, electrostatic, hydrophobic and H-bond acceptor factors were taken into consideration. The resulting fields were mapped on to the published crystal structure of the yeast enzyme and it was shown that the steric and hydrophobic fields are in good agreement with sulfonylurea-AHAS interaction geometry.

Unexpected Activity of a Novel Kunitz-Type Inhibtior: Inhibition of Cysteine Proteases but not Serine Proteases

Kunitz-type (KT) protease inhibitors are low molecular weight proteins classically defined as serine protease inhibitors. We identified a novel secreted KT inhibitor associated with the gut and parenchymal tissues of the infective juvenile stage of Fasciola hepatica, a helminth parasite of medical and veterinary importance. Unexpectedly, recombinant KT inhibitor (rFhKT1) exhibited no inhibitory activity towards serine proteases but was a potent inhibitor of the major secreted cathepsin L cysteine proteases of F. hepatica, FhCL1 and FhCL2, and of human cathepsins L and K (Ki = 0.4 nM - 27 nM). FhKT1 prevented the auto-catalytic activation of FhCL1 and FhCL2 and formed stable complexes with the mature enzymes. Pull-down experiments from adult parasite culture medium showed that rFhKT1 interacts specifically with native secreted FhCL1, FhCL2 and FhCL5. Substitution of the unusual P1 Leu15 within the exposed reactive loop of FhKT1 for the more commonly found Arg (FhKT1Leu15/Arg15) had modest adverse effects on the cysteine protease inhibition but conferred potent activity against the serine protease trypsin (Ki = 1.5 nM). Computational docking and sequence analysis provided hypotheses for the exclusive binding of FhKT1 to cysteine proteases, the importance of the Leu15 in anchoring the inhibitor into the S2 active site pocket, and the inhibitor's selectivity towards FhCL1, FhCL2 and human cathepsins L and K. FhKT1 represents a novel evolutionary adaptation of KT protease inhibitors by F. hepatica, with its prime purpose likely in the regulation of the major parasite-secreted proteases and/or cathepsin L-like proteases of its host.

Computational strategies to include protein flexibility in Ligand Docking and Virtual Screening

The dynamic character of proteins strongly influences biomolecular recognition mechanisms. With the development of the main models of ligand recognition (lock-and-key, induced fit, conformational selection theories), the role of protein plasticity has become increasingly relevant. In particular, major structural changes concerning large deviations of protein backbones, and slight movements such as side chain rotations are now carefully considered in drug discovery and development. It is of great interest to identify multiple protein conformations as preliminary step in a screening campaign. Protein flexibility has been widely investigated, in terms of both local and global motions, in two diverse biological systems. On one side, Replica Exchange Molecular Dynamics has been exploited as enhanced sampling method to collect multiple conformations of Lactate Dehydrogenase A (LDHA), an emerging anticancer target. The aim of this project was the development of an Ensemble-based Virtual Screening protocol, in order to find novel potent inhibitors. On the other side, a preliminary study concerning the local flexibility of Opioid Receptors has been carried out through ALiBERO approach, an iterative method based on Elastic Network-Normal Mode Analysis and Monte Carlo sampling. Comparison of the Virtual Screening performances by using single or multiple conformations confirmed that the inclusion of protein flexibility in screening protocols has a positive effect on the probability to early recognize novel or known active compounds.

In silico prediction of human carboxylesterase-1 (hCES1) metabolism combining docking analyses and MD simulations.

Metabolic problems lead to numerous failures during clinical trials, and much effort is now devoted in developing in silico models predicting metabolic stability and metabolites. Such models are well known for cytochromes P450 and some transferases, whereas little has been done to predict the hydrolytic activity of human hydrolases. The present study was undertaken to develop a computational approach able to predict the hydrolysis of novel esters by human carboxylesterase hCES1. The study involves both docking analyses of known substrates to develop predictive models, and molecular dynamics (MD) simulations to reveal the in situ behavior of substrates and products, with particular attention being paid to the influence of their ionization state. The results emphasize some crucial properties of the hCES1 catalytic cavity, confirming that as a trend with several exceptions, hCES1 prefers substrates with relatively smaller and somewhat polar alkyl/aryl groups and larger hydrophobic acyl moieties. The docking results underline the usefulness of the hydrophobic interaction score proposed here, which allows a robust prediction of hCES1 catalysis, while the MD simulations show the different behavior of substrates and products in the enzyme cavity, suggesting in particular that basic substrates interact with the enzyme in their unprotonated form.

Docking, virtual high throughput screening and in silico fragment-based drug design.

ABSTRACT The drug discovery process has been profoundly changed recently by the adoption of computational methods helping the design of new drug candidates more rapidly and at lower costs. In silico drug design consists of a collection of tools helping to make rational decisions at the different steps of the drug discovery process, such as the identification of a biomolecular target of therapeutical interest, the selection or the design of new lead compounds and their modification to obtain better affinities, as well as pharmacokinetic and pharmacodynamic properties. Among the different tools available, a particular emphasis is placed in this review on molecular docking, virtual high throughput screening and fragment-based ligand design.

Homology modeling and metabolism prediction of human carboxylesterase-2 using docking analyses by GriDock: a parallelized tool based on AutoDock 4.0.

Metabolic problems lead to numerous failures during clinical trials, and much effort is now devoted to developing in silico models predicting metabolic stability and metabolites. Such models are well known for cytochromes P450 and some transferases, whereas less has been done to predict the activity of human hydrolases. The present study was undertaken to develop a computational approach able to predict the hydrolysis of novel esters by human carboxylesterase hCES2. The study involved first a homology modeling of the hCES2 protein based on the model of hCES1 since the two proteins share a high degree of homology (congruent with 73%). A set of 40 known substrates of hCES2 was taken from the literature; the ligands were docked in both their neutral and ionized forms using GriDock, a parallel tool based on the AutoDock4.0 engine which can perform efficient and easy virtual screening analyses of large molecular databases exploiting multi-core architectures. Useful statistical models (e.g., r (2) = 0.91 for substrates in their unprotonated state) were calculated by correlating experimental pK(m) values with distance between the carbon atom of the substrate's ester group and the hydroxy function of Ser228. Additional parameters in the equations accounted for hydrophobic and electrostatic interactions between substrates and contributing residues. The negatively charged residues in the hCES2 cavity explained the preference of the enzyme for neutral substrates and, more generally, suggested that ligands which interact too strongly by ionic bonds (e.g., ACE inhibitors) cannot be good CES2 substrates because they are trapped in the cavity in unproductive modes and behave as inhibitors. The effects of protonation on substrate recognition and the contrasting behavior of substrates and products were finally investigated by MD simulations of some CES2 complexes.

Computational and Biological Evaluation of N-octadecyl-N'-propylsulfamide, a Selective PPARα Agonist Structurally Related to N-acylethanolamines.

To further understand the pharmacological properties of N-oleoylethanolamine (OEA), a naturally occurring lipid that activates peroxisome proliferator-activated receptor alpha (PPARα), we designed sulfamoyl analogs based on its structure. Among the compounds tested, N-octadecyl-N'-propylsulfamide (CC7) was selected for functional comparison with OEA. The performed studies include the following computational and biological approaches: 1) molecular docking analyses; 2) molecular biology studies with PPARα; 3) pharmacological studies on feeding behavior and visceral analgesia. For the docking studies, we compared OEA and CC7 data with crystallization data obtained with the reference PPARα agonist GW409544. OEA and CC7 interacted with the ligand-binding domain of PPARα in a similar manner to GW409544. Both compounds produced similar transcriptional activation by in vitro assays, including the GST pull-down assay and reporter gene analysis. In addition, CC7 and OEA induced the mRNA expression of CPT1a in HpeG2 cells through PPARα and the induction was avoided with PPARα-specific siRNA. In vivo studies in rats showed that OEA and CC7 had anorectic and antiobesity activity and induced both lipopenia and decreases in hepatic fat content. However, different effects were observed when measuring visceral pain; OEA produced visceral analgesia whereas CC7 showed no effects. These results suggest that OEA activity on the PPARα receptor (e.g., lipid metabolism and feeding behavior) may be dissociated from other actions at alternative targets (e.g., pain) because other non cannabimimetic ligands that interact with PPARα, such as CC7, do not reproduce the full spectrum of the pharmacological activity of OEA. These results provide new opportunities for the development of specific PPARα-activating drugs focused on sulfamide derivatives with a long alkyl chain for the treatment of metabolic dysfunction.

Computational drug design strategies applied to the modelling of human immunodeficiency virus-1 reverse transcriptase inhibitors

Reverse transcriptase (RT) is a multifunctional enzyme in the human immunodeficiency virus (HIV)-1 life cycle and represents a primary target for drug discovery efforts against HIV-1 infection. Two classes of RT inhibitors, the nucleoside RT inhibitors (NRTIs) and the nonnucleoside transcriptase inhibitors are prominently used in the highly active antiretroviral therapy in combination with other anti-HIV drugs. However, the rapid emergence of drug-resistant viral strains has limited the successful rate of the anti-HIV agents. Computational methods are a significant part of the drug design process and indispensable to study drug resistance. In this review, recent advances in computer-aided drug design for the rational design of new compounds against HIV-1 RT using methods such as molecular docking, molecular dynamics, free energy calculations, quantitative structure-activity relationships, pharmacophore modelling and absorption, distribution, metabolism, excretion and toxicity prediction are discussed. Successful applications of these methodologies are also highlighted.