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.

Flexibility and inhibitor binding in Cdc25 phosphatases

Cdc25 phosphatases involved in cell cycle checkpoints are now active targets for the development of anti-cancer therapies. Rational drug design would certainly benefit from detailed structural information for Cdc25s. However, only apo- or sulfate-bound crystal structures of the Cdc25 catalytic domain have been described so far. Together with previously available crystalographic data, results from molecular dynamics simulations, bioinformatic analysis, and computer-generated conformational ensembles shown here indicate that the last 30-40 residues in the C-terminus of Cdc25B are partially unfolded or disordered in solution. The effect of C-terminal flexibility upon binding of two potent small molecule inhibitors to Cdc25B is then analyzed by using three structural models with variable levels of flexibility, including an equilibrium distributed ensemble of Cdc25B backbone conformations. The three Cdc25B structural models are used in combination with flexible docking, clustering, and calculation of binding free energies by the linear interaction energy approximation to construct and validate Cdc25B-inhibitor complexes. Two binding sites are identified on top and beside the Cdc25B active site. The diversity of interaction modes found increases with receptor flexibility. Backbone flexibility allows the formation of transient cavities or compact hydrophobic units on the surface of the stable, folded protein core that are unexposed or unavailable for ligand binding in rigid and densely packed crystal structures. The present results may help to speculate on the mechanisms of small molecule complexation to partially unfolded or locally disordered proteins.

Structure-activity relationships of hypervalent organochalcogenanes as inhibitors of cysteine cathepsins V and S

A new series of organotelluranes were synthesized and investigated, and the structure-activity relationships in cysteine proteases inhibition were determinated. It was possible to identify the relevance of structural components linked to the reactivity of these compounds as inhibitors. For example, dibromo-organotelluranes showed to be more reactive than dichloro-organotelluranes towards cysteine cathepsins V and S. Besides, no remarkable enantio-selectivity was verified. In general the achiral organotelluranes were more reactive than the chiral congeners against cysteine cathepsins V and S. A reactivity order for organochalcogenanes and cysteine cathepsins was proposed after the comparison of the inhibitory potencies of organotelluranes with the related organoselenanes. (C) 2011 Elsevier Ltd. All rights reserved.

Solvation in Pure Liquids: What Can Be Learned from the Use of Pairs of Indicators?

The solvation of six solvatochromic probes in a large number of solvents (33-68) was examined at 25 degrees C. The probes employed were the following: 2,6-diphenyl-4-(2,4,6-triphenylpyridinium-1-yl) phenolate (RB); 4-[(E)2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePM; 1-methylquinolinium-8-olate, QB; 2-bromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr, 2,6-dichloro-4-(2,4,6-triphenyl pyridinium-1-yl) phenolate (WB); and 2,6-dibromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr(2), respectively. Of these, MePMBr is a novel compound. They can be grouped in three pairs, each with similar pK(a) in water but with different molecular properties, for example, lipophilicity and dipole moment. These pairs are formed by RB and MePM; QB and MePMBr; WB and MePMBr(2), respectively. Theoretical calculations were carried out in order to calculate their physicochemical properties including bond lengths, dihedral angles, dipole moments, and wavelength of absorption of the intramolecular charge-transfer band in four solvents, water, methanol, acetone, and DMSO, respectively. The data calculated were in excellent agreement with available experimental data, for example, bond length and dihedral angles. This gives credence to the use of the calculated properties in explaining the solvatochromic behaviors observed. The dependence of an empirical solvent polarity scale E(T)(probe) in kcal/mol on the physicochemical properties of the solvent (acidity, basicity, and dipolarity/polarizability) and those of the probes (pK(a), and dipole moment) was analyzed by using known multiparameter solvation equations. For each pair of probes, values of E(T)(probe) (for example, E(T)(MePM) versus E(T)(RB)) were found to be linearly correlated with correlation coefficients, r, between 0.9548 and 0.9860. For the mercyanine series, the values of E(T)(probe) also correlated linearly, with (r) of 0.9772 (MePMBr versus MePM) and 0.9919 (MePMBr(2) versus MePM). The response of each pair of probes (of similar pK(a)) to solvent acidity is the same, provided that solute-solvent hydrogen-bonding is not seriously affected by steric crowding (as in case of RB). We show, for the first time, that the response to solvent dipolarity/polarizability is linearly correlated to the dipole moment of the probes. The successive introduction of bromine atoms in MePM (to give MePMBr, then MePMBr(2)) leads to the following linear decrease: pK(a) in water, length of the phenolate oxygen-carbon bond, length of the central ethylenic bond, susceptibility to solvent acidity, and susceptibility to solvent dipolarity/polarizability. Thus studying the solvation of probes whose molecular structures are varied systematically produces a wealth of information on the effect of solute structure on its solvation. The results of solvation of the present probes were employed in order to test the goodness of fit of two independent sets of solvent solvatochromic parameters.

Complete assignment of (1)H and (13)C NMR spectra of alpha-phenylsulfinyl-N-methoxy-N-methylpropionamide and some p-substituted derivatives

The complete assignment of the (1)H and (13)C NMR spectra of the diastereomeric pairs of some alpha-arylsulfinyl-substituted N-methoxy-N-methylpropionamides with the substituents methoxy, methyl, chloro, nitro is reported. Copyright (C) 2008 John Wiley & Sons, Ltd.

Crystal and molecular structures of three 2-sulphur-substituted cyclohexanones studied by X-ray crystallography and by ab initio molecular orbital calculations

The conformational features of three 2-sulphur-substituted cyclohexanone derivatives, which differ in the number of sulphur-bound oxygen atoms, i.e. zero (I), one (II) and two (III), were investigated by single crystal X-ray crystallography and geometry optimized structures determined using Hartree-Fock method. In each of (I)-(III) an intramolecular S center dot center dot center dot O(carbonyl) interaction is found with the magnitude correlated with the oxidation state of the sulphur atom, i.e. 2.838(3) angstrom in (I) to 2.924(2) angstrom in (II) to 3.0973(18) angstrom in (III). There is an inverse relationship between the strength of this interaction and the magnitude of the carbonyl bond. The supramolecular aggregation patterns are primarily determined by C-H center dot center dot center dot O contacts and are similarly influenced by the number of oxygen atoms in the molecular structures. Thus, a supramolecular chain is found in the crystal structure of (I). With an additional oxygen atom available to participate in C-H center dot center dot center dot O interactions, as in (II), a two-dimensional array is found. Finally, a three-dimensional network is found for (III). Despite there being differences in conformations between the experimental structures and those calculated in the gas-phase, the S center dot center dot center dot O interactions persist. The presence of intermolecular C-H center dot center dot center dot O interactions involving the cyclohexanone-carbonyl group in the solid-state, disrupts the stabilising intramolecular C-H center dot center dot center dot O interaction in the energetically-favoured conformation. (I): C(12)H(13)NO(3)S, triclinic space group P (1) over bar with a = 5.392(3) angstrom b = 10.731(6) angstrom, c = 11.075(6) angstrom, alpha = 113.424(4)degrees, beta = 94.167(9)degrees, gamma = 98.444(6)degrees, V = 575.5(6) angstrom(3), Z = 2, R(1) = 0.052; (II): C(12)H(13)NO(4)S, monoclinic P2(1)/n, a = 7.3506(15) angstrom, b = 6.7814(14) angstrom, c = 23.479(5) angstrom, beta = 92.94(3)degrees, V = 1168.8(4) angstrom(3), Z = 4, R(1) = 0.046; (III): C(12)H(13)NO(5)S, monoclinic P2(1)/c, a = 5.5491(11) angstrom, b = 24.146(3) angstrom, c = 11.124(3) angstrom, beta = 114.590(10)degrees, V = 1355.3(5) angstrom(3), Z = 4, R(1) = 0.051.

Stereochemical and electronic interaction studies of some N-methoxy-N-methyl-2-[(4 `-substituted)phenylsulfinyl]propanamides

The (1)H NMR spectra of N-methoxy-N-methyl-2-[(4`-substituted)phenylsulfinyl]-propanamides [Y-Ph-S(O)CH(Me)C(O)N(OMe)Me; Y = OMe 1, Me 2, H 3. Cl 4, NO(2) 5] along with the X-ray diffraction analysis of the nitro-derivative (5). have shown the existence of two pairs of diastereomers (racemic mixture) [C(R)S(S)/C(S)S(R) (diast(1)) and C(R)S(R)/C(S)S(S) (diast(2))] in the ratio of ca. 7:3. respectively. The v(CO) IR analysis of the title compounds supported by HF and B3LYP/6-31G** calculations of 3 and of the parent N-methoxy-N-methyl-propanamide (6) by HF, have shown that diast(1) exists in an equilibrium between the two more polar and more stable quasi-cis (q-c(1) and q-c(2)) conformers and the gauche(g) conformer. The population of the g conformer in the equilibrium increases with the increase in the solvent polarity, which is attributed to a larger solvation effect on the carbonyl and sulfinyl groups. Diast(2) of compound 3 occurs in the gas phase as an equilibrium between the most stable quasi-gauche (q-g) conformer and the quasi-cis (q-c) conformer, both presenting very similar dipole moments. The former is stabilized by electrostatic and charge transfer interactions, which results in a less solvated spatial arrangement. Moreover, all conformers of both diastereomers are stabilized by several intramolecular hydrogen bonds. X-ray single crystal analysis performed for diast(1) and for diast(2) of 5 indicates that both stereoisomers assume, in the solid state, the anti-clinal (gauche) conformation. For the crystal packing, diast(1) of 5 is made up of three molecules joined through two centro-symmetric H center dot center dot center dot O hydrogen bonds. (C) 2008 Elsevier B.V. All rights reserved.

Synthesis, Spectroscopic Studies and X-ray Crystal Structures of New Pyrazoline and Pyrazole Derivatives

The synthesis and characterization of some pyrazoline compounds of 1,3-diketones with hydrazine derivatives, namely, 1-(S-benzyldithiocarbazate)-3-methyl-5-phenyl-5-hydroxypyrazoline (1); 1-(2-thiophenecarboxylic)-3-methyl-5-phenyl-5-hydroxypyrazoline (2); 1-(2-thiophenecarboxylic)-3,5-dimethyl-5-hydroxypyrazoline (3); 1-(S-benzyldithiocarbazato)-3-methyl-5-phenylpyrazole (4); 1-(2-thiophenecarboxylic)-3-methyl-5-phenylpyrazole (5) and 1-(S-benzyldithiocarbazate)-3,5-dimethylpyrazole (6) are reported. Studies by IR, ((1)H, (13)C)-NMR spectroscopies and single crystal X-ray diffraction revealed that compounds (1)(,) (2) and (3) are formed as pyrazoline, whereas (4) and (5) are formed as pyrazole derivatives only under acidic conditions. Compound (1) crystallizes in orthorhombic P2(1)2(1)2(1), a = 6.38960(10) angstrom, b = 12.9176(3) angstrom, c = 21.2552(5) angstrom, (2) crystallizes in monoclinic, P2(1)/n, a = 11.3617(2) angstrom, b = 8.4988(2) angstrom, c = 92.8900(10)angstrom and beta = 92.8900(5)degrees, (3) crystallizes in monoclinic, C2/c, a = 15.9500(5) angstrom, b = 9.3766(3) angstrom, c = 16.6910(5)angstrom and beta = 113.825(2)degrees, (4) crystallizes in monoclinic, P2(1)/c, a = 15.228(4) angstrom, b = 5.5714(13) angstrom, c = 19.956(5)angstrom and beta = 91.575(7)degrees and (6) crystallizes in orthorhombic, P2(1)2(1)2(1), a = 5.3920(2) angstrom, b = 11.2074(5) angstrom, c = 21.885(1)angstrom . The (3) derivative represents the first pyrazoline compound prepared from 2,4-pentanedione and characterized crystallographically.

Novel insights for dihydroorotate dehydrogenase class 1A inhibitors discovery

The enzyme dihydroorotate dehydrogenase (DHODH) has been suggested as a promising target for the design of trypanocidal agents. We report here the discovery of novel inhibitors of Trypanosoma cruzi DHODH identified by a combination of virtual screening and ITC methods. Monitoring of the enzymatic reaction in the presence of selected ligands together with structural information obtained from X-ray crystallography analysis have allowed the identification and validation of a novel site of interaction (S2 site). This has provided important structural insights for the rational design of T cruzi and Leishmania major DHODH inhibitors. The most potent compound (1) in the investigated series inhibits TcDHODH enzyme with K(i)(app) value of 19.28 mu M and possesses a ligand efficiency of 0.54 kcal mol(-1) per non-H atom. The compounds described in this work are promising hits for further development. (C) 2010 Elsevier Masson SAS. All rights reserved.

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. (C) 2010 Elsevier Inc. All rights reserved.

Thiosemicarbazones, semicarbazones, dithiocarbazates and hydrazide/hydrazones: Anti-Mycobacterium tuberculosis activity and cytotoxicity

The aim of this study was to identify a candidate drug for the development of anti-tuberculosis therapy from previously synthesized compounds based on the thiosemicarbazones, semicarbazones, dithio-carbazates and hydrazide/hydrazones compounds. The minimal inhibitory concentration (MIC) of these compounds against Mycobacterium tuberculosis was determined. Their in vitro cytotoxicity to J774 cells (IC(50)) was determined to establish a selectivity index (SI) (SI = IC(50)/MIC). The best compounds were the thiosemicarbazones (2, 3 and 4) and the hydrazide/hydrazones (14, 15, 16 and 18). The results are comparable to or better than those of ""first line"" or ""second line"" drugs commonly used to treat TB, suggesting these compounds as anti-TB drug candidates. (C) 2010 Elsevier Masson SAS. All rights reserved.

Thermal behavior and structural properites of plant-derived eugenyl acetate

Eugenol is an allyl chain-substituted guaiacol in the biosynthesized phenylpropanoid compound class derived from Syzygium aromaticum L. and widely used in folk medicine. Nonetheless, its pharmacological use is limited by some problems, such as instability when exposed to light and high temperature. In order to enhance stability, the eugenol molecule was structurally modified, resulting in eugenyl acetate. The eugenyl acetate`s thermal behavior and crystal structure was then characterized by differential scanning calorimetry (DSC) and X-ray diffraction (XRD) and compared to a commercial sample.

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.

On the temperature stability of extracellular hemoglobin of Glossoscolex paulistus, at different oxidation states: SAXS and DLS studies

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Structural and Biochemical Characterization of Peroxiredoxin Q beta from Xylella fastidiosa CATALYTIC MECHANISM and HIGH REACTIVITY

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)