Structural, thermodynamic and functional studies of human 71 kDa heat shock cognate protein (HSPA8/hHsc70).

Human 71 kDa heat shock cognate protein (HSPA8, also known as Hsc70, Hsp70-8, Hsc71, Hsp71 or Hsp73) is a constitutively expressed chaperone that is critical for cell proteostasis. In the cytosol, HSPA8 plays a pivotal role in folding and refolding, facilitates protein trafficking across membranes and targets proteins for degradation, among other functions. Here, we report an in solution study of recombinant HSPA8 (rHSPA8) using a variety of biophysical and biochemical approaches. rHSPA8 shares several structural and functional similarities with others human Hsp70s. It has two domains with different stabilities and interacts with adenosine nucleotides with dissociation constants in the low micromolar range, which were higher in the presence of Mg2+. rHSPA8 showed lower ATPase activity than its homolog HSPA5/hGrp78/hBiP, but it was 4-fold greater than that of recombinant HSPA1A/hHsp70-1A, with which it is 86% identical. Small angle X-ray scattering indicated that rHSPA8 behaved as an elongated monomeric protein in solution with dimensions similar to those observed for HSPA1A. In addition, rHSPA8 showed structural flexibility between its compacted and extended conformations. The data also indicated that HSPA8 has capacity in preventing the aggregation of model client proteins. The present study expands the understanding of the structure and activity of this chaperone and aligns with the idea that human homologous Hsp70s have divergent functions.

Synthesis, biochemical evaluation and molecular modeling studies of nonpeptidic nitrile-based fluorinated compounds.

Aim: Compounds that block enzyme activity can kill pathogens and help develop effective and safe drugs for Chagas disease and leishmaniasis. Materials & methods: A library of nonpeptidic nitrile-based compounds was synthesized and had their inhibitory affinity tested against cruzain, Leishmania mexicana cysteine protease B and cathepsin L. Isothermal titration calorimetry experiments and molecular simulations were performed for selected compounds to obtain thermodynamic fingerprints and identify main interactions and putative modes of binding with cruzain. Results: The derivatives provided increased affinity against all enzymes compared with the lead, and thermodynamic and computational studies showed improved thermodynamic properties and a possible different mode of binding. Conclusion: Our studies culminated in 1b, a compound 60-fold more potent in cruzain than its lead that also showed entropic and enthalpic contributions favorable to Gibbs binding energy.

Synthesis and matched molecular pair analysis of covalent reversible inhibitors of the cysteine protease CPB.

Cysteine protease B (CPB) can be targeted by reversible covalent inhibitors that could serve as antileishmanial compounds. Here, sixteen dipeptidyl nitrile derivatives were synthesized, tested against CPB, and analyzed using matched molecular pairs to determine the effects of stereochemistry and p-phenyl substitution on enzyme inhibition. The compound (S)-2-(((S)-1-(4-bromophenyl)-2,2,2-trifluoroethyl)amino)-N-(1-cyanocyclopropyl)-3-phenylpropanamide (5) was the most potent CPB inhibitor (pKi = 6.82), which was also selective for human cathepsin B (pKi < 5). The inversion of the stereochemistry from S to R was more detrimental to potency when placed at the P2 position than at P3. The p-Br derivatives were more potent than the p-CH3 and p-OCH3 derivatives, probably due to intermolecular interactions with the S3 subsite.

Novel anti-prostate cancer scaffold identified by the combination of in silico and cell-based assays targeting the PI3K-AKT-mTOR pathway.

Phenotypic assays were performed in prostate cancer cell lines to describe the biological activity of PI3K-AKT-mTOR pathway inhibitors retrieved from the virtual screening initiative. These novel chemicals share in common the aminopyridine scaffold, hitting PC-3 cells in macromolar range, with selectivity index over fibroblast cell lines. Moreover, a preliminary study of the mode of action by flow cytometry assay pointed out that these compounds had a rapamycin-like response for the PI3K-AKT-mTOR pathway modulation.

Highly predictive hologram QSAR models of nitrile-containing cruzain inhibitors.

The HQSAR, molecular docking, and ROCS were applied to a data-set of 57 cruzain inhibitors. The best HQSAR model (q2 = .70, r2 = .95, [Formula: see text] = .62, [Formula: see text] = .09 and [Formula: see text] = .26), employing well-balanced, diverse training (40) and test (17) sets, was obtained using atoms (A), bonds (B), and hydrogen (H) as fragment distinctions and 6-9 as fragment sizes. This model was then used to predict the unknown potencies of 121 compounds (the V1 database), giving rise to a satisfactory predictive r2 value of .65 (external validation). By employing an extra external data-set comprising 1223 compounds (the V3 database) either retrieved from the ChEMBL or CDD databases, an overall ROC AUC score well over .70 was obtained. The contribution maps obtained with the best HQSAR model (model 3.4) are in agreement with the predicted binding mode and with the biological potencies of the studied compounds. We also screened these compounds using the ROCS method, a Gaussian-shape volume filter able to identify quickly the shapes that match a query molecule. The area under the curve (AUC) obtained with the ROC curves (ROC AUC) was .72, indicating that the method was very efficient in distinguishing between active and inactive cruzain inhibitors. These set of information guided us to propose novel cruzain inhibitors to be synthesized. Then, the best HQSAR model obtained was used to predict the pIC50 values of these new compounds. Some compounds identified using this method have shown calculated potencies higher than those which have originated them.

In silico selection and cell-based characterization of selective and bioactive compounds for androgen-dependent prostate cancer cell.

Prostate cancer is one of the most prevalent types of cancer in male population. It is a hormone driven disease, especially in its initial phase. Hence, androgen deprivation therapy (ADT) is the major chemotherapeutic effort and novel AR inhibitors with improved pharmacological profiles are needed. In this report, a novel bioactive compound was selected and investigated using in silico and cell-based assays. Neq0502 compound was selective for the testosterone stimulated AR-dependent prostate cancer cell (LNCaP, GI50=22.4µM) when compared with unstimulated LNCaP or AR-insensitive (DU145 and PC-3) cell lines. Cell cycle arrest study provided the same profile for Neq0502 and the reference drug enzalutamide. Moreover, this compound is not cytotoxic for fibroblast Balb/C 3T3 clone A31 cells up to 250µM, with a good selectivity ratio (SI>11), which could be used in compound optimization effort to a novel therapeutic alternative.

The Open Form Inducer Approach for Structure-Based Drug Design.

Many open form (OF) structures of drug targets were obtained a posteriori by analysis of co-crystals with inhibitors. Therefore, obtaining the OF structure of a drug target a priori will accelerate development of potent inhibitors. In addition to its small active site, Trypanosoma cruzi dihydroorotate dehydrogenase (TcDHODH) is fully functional in its monomeric form, making drug design approaches targeting the active site and protein-protein interactions unrealistic. Therefore, a novel a priori approach was developed to determination the TcDHODH active site in OF. This approach consists of generating an "OF inducer" (predicted in silico) to bind the target and cause steric repulsion with flexible regions proximal to the active site that force it open. We provide the first proof-of-concept of this approach by predicting and crystallizing TcDHODH in complex with an OF inducer, thereby obtaining the OF a priori with its subsequent use in designing potent and selective inhibitors. Fourteen co-crystal structures of TcDHODH with the designed inhibitors are presented herein. This approach has potential to encourage drug design against diseases where the molecular targets are such difficult proteins possessing small AS volume. This approach can be extended to study open/close conformation of proteins in general, the identification of allosteric pockets and inhibitors for other drug targets where conventional drug design approaches are not applicable, as well as the effective exploitation of the increasing number of protein structures deposited in Protein Data Bank.

Quantitative structure-retention relationships of flavonoids unraveled by immobilized artificial membrane chromatography.

The pharmacokinetic properties of flavonoids with differing degrees of lipophilicity were investigated using immobilized artificial membranes (IAMs) as the stationary phase in high performance liquid chromatography (HPLC). For each flavonoid compound, we investigated whether the type of column used affected the correlation between the retention factors and the calculated octanol/water partition (log Poct). Three-dimensional (3D) molecular descriptors were calculated from the molecular structure of each compound using i) VolSurf software, ii) the GRID method (computational procedure for determining energetically favorable binding sites in molecules of known structure using a probe for calculating the 3D molecular interaction fields, between the probe and the molecule), and iii) the relationship between partition and molecular structure, analyzed in terms of physicochemical descriptors. The VolSurf built-in Caco-2 model was used to estimate compound permeability. The extent to which the datasets obtained from different columns differ both from each other and from both the calculated log Poct and the predicted permeability in Caco-2 cells was examined by principal component analysis (PCA). The immobilized membrane partition coefficients (kIAM) were analyzed using molecular descriptors in partial least square regression (PLS) and a quantitative structure-retention relationship was generated for the chromatographic retention in the cholesterol column. The cholesterol column provided the best correlation with the permeability predicted by the Caco-2 cell model and a good fit model with great prediction power was obtained for its retention data (R(2)=0.96 and Q(2)=0.85 with four latent variables).

Evaluation of the leishmanicidal and cytotoxic effects of inhibitors for microorganism metabolic pathway enzymes.

Chemotherapy for leishmaniosis a neglected parasitic disease, is based on few drugs, which are toxic and present resistance issues. Efforts for the development of new therapies are essential for the control of leishmaniasis. Metabolic pathway enzymes are promising targets for new drugs against parasites. The search for effective drugs against key enzymes can take advantage of the similarities between metabolic pathways in different microorganisms trypanosomatids Trypanosoma cruzi and Leishmania and fungus Saccharomyces cerevisiae. In this report, leishmanicidal activity of the metabolic pathway enzymes inhibitors (IDs) of dihydroorotate dehydrogenase (DHODH), glyceraldehyde 3-phosphate dehydrogenase and cruzain-cysteine protease from T. cruzi and scitalona-desidratase, adenosine deaminase, succinate dehydrogenase complex II and hydroxynaphthalene reductase from S. cerevisiae was performed on Leishmania amazonensis extracellular promastigotes and amastigotes within macrophages. The most promising compound, ID195, which is a DHODH inhibitor was toxic against promastigotes and was selective for amastigotes over host cells.

Pharmacophore modeling for anti-Chagas drug design using the fragment molecular orbital method.

BACKGROUND: Chagas disease, caused by the parasite Trypanosoma cruzi, is a neglected tropical disease that causes severe human health problems. To develop a new chemotherapeutic agent for the treatment of Chagas disease, we predicted a pharmacophore model for T. cruzi dihydroorotate dehydrogenase (TcDHODH) by fragment molecular orbital (FMO) calculation for orotate, oxonate, and 43 orotate derivatives. METHODOLOGY/PRINCIPAL FINDINGS: Intermolecular interactions in the complexes of TcDHODH with orotate, oxonate, and 43 orotate derivatives were analyzed by FMO calculation at the MP2/6-31G level. The results indicated that the orotate moiety, which is the base fragment of these compounds, interacts with the Lys43, Asn67, and Asn194 residues of TcDHODH and the cofactor flavin mononucleotide (FMN), whereas functional groups introduced at the orotate 5-position strongly interact with the Lys214 residue. CONCLUSIONS/SIGNIFICANCE: FMO-based interaction energy analyses revealed a pharmacophore model for TcDHODH inhibitor. Hydrogen bond acceptor pharmacophores correspond to Lys43 and Lys214, hydrogen bond donor and acceptor pharmacophores correspond to Asn67 and Asn194, and the aromatic ring pharmacophore corresponds to FMN, which shows important characteristics of compounds that inhibit TcDHODH. In addition, the Lys214 residue is not conserved between TcDHODH and human DHODH. Our analysis suggests that these orotate derivatives should preferentially bind to TcDHODH, increasing their selectivity. Our results obtained by pharmacophore modeling provides insight into the structural requirements for the design of TcDHODH inhibitors and their development as new anti-Chagas drugs.

Conformational differences between the wild type and V30M mutant transthyretin modulate its binding to genistein: implications to tetramer stability and ligand-binding.

Transthyretin (TTR) is a tetrameric beta-sheet-rich transporter protein directly involved in human amyloid diseases. It was recently found that the isoflavone genistein (GEN) potently inhibits TTR amyloid fibril formation (Green et al., 2005) and is therefore a promising candidate for TTR amyloidosis treatment. Here we used structural and biophysical approaches to characterize genistein binding to the wild type (TTRwt) and to its most frequent amyloidogenic variant, the V30M mutant. In a dose-dependent manner, genistein elicited considerable increases in both mutant and TTRwt stability as demonstrated by high hydrostatic pressure (HHP) and acid-mediated dissociation/denaturation assays. TTR:GEN crystal complexes and isothermal titration calorimetry (ITC) experiments showed that the binding mechanisms of genistein to the TTRwt and to V30M are different and are dependent on apoTTR structure conformations. Furthermore, we could also identify potential allosteric movements caused by genistein binding to the wild type TTR that explains, at least in part, the frequently observed negatively cooperative process between the two sites of TTRwt when binding ligands. These findings show that TTR mutants may present different ligand recognition and therefore are of value in ligand design for inhibiting TTR amyloidosis.

Kinetic mechanism and catalysis of Trypanosoma cruzi dihydroorotate dehydrogenase enzyme evaluated by isothermal titration calorimetry.

Trypanosoma cruzi dihydroorotate dehydrogenase (TcDHODH) catalyzes the oxidation of l-dihydroorotate to orotate with concomitant reduction of fumarate to succinate in the de novo pyrimidine biosynthetic pathway. Based on the important need to characterize catalytic mechanism of TcDHODH, we have tailored a protocol to measure TcDHODH kinetic parameters based on isothermal titration calorimetry. Enzymatic assays lead to Michaelis-Menten curves that enable the Michaelis constant (K(M)) and maximum velocity (V(max)) for both of the TcDHODH substrates: dihydroorotate (K(M)=8.6+/-2.6 microM and V(max)=4.1+/-0.7 microMs(-1)) and fumarate (K(M)=120+/-9 microM and V(max)=6.71+/-0.15 microMs(-1)). TcDHODH activity was investigated using dimethyl sulfoxide (10%, v/v) and Triton X-100 (0.5%, v/v), which seem to facilitate the substrate binding process with a small decrease in K(M). Arrhenius plot analysis allowed the determination of thermodynamic parameters of activation for substrates and gave some insights into the enzyme mechanism. Activation entropy was the main contributor to the Gibbs free energy in the formation of the transition state. A factor that might contribute to the unfavorable entropy is the hindered access of substrates to the TcDHODH active site where a loop at its entrance regulates the open-close channel for substrate access.

Docking and molecular dynamics simulation of quinone compounds with trypanocidal activity.

In this work, two different docking programs were used, AutoDock and FlexX, which use different types of scoring functions and searching methods. The docking poses of all quinone compounds studied stayed in the same region in the trypanothione reductase. This region is a hydrophobic pocket near to Phe396, Pro398 and Leu399 amino acid residues. The compounds studied displays a higher affinity in trypanothione reductase (TR) than glutathione reductase (GR), since only two out of 28 quinone compounds presented more favorable docking energy in the site of human enzyme. The interaction of quinone compounds with the TR enzyme is in agreement with other studies, which showed different binding sites from the ones formed by cysteines 52 and 58. To verify the results obtained by docking, we carried out a molecular dynamics simulation with the compounds that presented the highest and lowest docking energies. The results showed that the root mean square deviation (RMSD) between the initial and final pose were very small. In addition, the hydrogen bond pattern was conserved along the simulation. In the parasite enzyme, the amino acid residues Leu399, Met400 and Lys402 are replaced in the human enzyme by Met406, Tyr407 and Ala409, respectively. In view of the fact that Leu399 is an amino acid of the Z site, this difference could be explored to design selective inhibitors of TR.

A neural networks study of quinone compounds with trypanocidal activity.

This work investigates neural network models for predicting the trypanocidal activity of 28 quinone compounds. Artificial neural networks (ANN), such as multilayer perceptrons (MLP) and Kohonen models, were employed with the aim of modeling the nonlinear relationship between quantum and molecular descriptors and trypanocidal activity. The calculated descriptors and the principal components were used as input to train neural network models to verify the behavior of the nets. The best model for both network models (MLP and Kohonen) was obtained with four descriptors as input. The descriptors were T5 (torsion angle), QTS1 (sum of absolute values of the atomic charges), VOLS2 (volume of the substituent at region B) and HOMO-1 (energy of the molecular orbital below HOMO). These descriptors provide information on the kind of interaction that occurs between the compounds and the biological receptor. Both neural network models used here can predict the trypanocidal activity of the quinone compounds with good agreement, with low errors in the testing set and a high correctness rate. Thanks to the nonlinear model obtained from the neural network models, we can conclude that electronic and structural properties are important factors in the interaction between quinone compounds that exhibit trypanocidal activity and their biological receptors. The final ANN models should be useful in the design of novel trypanocidal quinones having improved potency.

QSAR and molecular modelling studies on B-DNA recognition of minor groove binders.

Aromatic bisamidines have been proved to be efficient compounds against Leishmania spp. and Pneumocystis carinii. Although the mode of action is still not known, these molecules are supposed to be DNA minor groove binders (MGBs). This paper describes a molecular modelling study for a set of MGBs in order to rank them through their complementarity to the Dickerson Drew Dodecamer (DDD) according to their interaction energies with B-DNA. A comparative molecular field analysis (CoMFA) has shown the importance of relatively bulky positively charged groups attached to the MGB aromatic rings, and small and negatively charged substituents into the middle chain. Models were obtained for DNA denaturation related to H-bonding processes of binding modes. Validation of the model demonstrated the robustness of CoMFA in terms of independent test set of similar MGBs. GRID results allotted bioisosteric substitution of z.sbnd;Oz.sbnd; by z.sbnd;NHz.sbnd; in furan ring of furamidine and related compounds as being capable to enhance the binding to DDD.

HEPT derivatives as non-nucleoside inhibitors of HIV-1 reverse transcriptase: QSAR studies agree with the crystal structures.

The interest in the non-nucleoside inhibitors (NNIs) to the reverse transcriptase (RT) as anti-AIDS agents has grown in the last ten years. The compound 1-[(2-hydroxyethoxy)methyl]-6-(phenylthio)thymine (HEPT) is the precursor of the most studied class of NNIs, from which hundreds of derivatives have been synthesized and tested. There are at least twelve QSAR studies about the HEPT derivatives as RT inhibitors. Most of the predictions derived by these studies are related to the nature of the active site near the substituents at positions N-1 and C-5, and at the C-6 phenyl ring. The validity of these models has been checked against the 3-D structure of HIV 1 RT-HEPT complexes available. Most of these predictions were confirmed at the molecular level.