Systematic structural studies of iron superoxide dismutases from human parasites and a statistical coupling analysis of metal binding specificity

Superoxide dismutases (SODs) are a crucial class of enzymes in the combat against intracellular free radical damage. They eliminate superoxide radicals by converting them into hydrogen peroxide and oxygen. In spite of their very different life cycles and infection strategies, the human parasites Plasmodium falciparum, Trypanosoma cruzi and Trypanosoma brucei are known to be sensitive to oxidative stress. Thus the parasite Fe-SODs have become attractive targets for novel drug development. Here we report the crystal structures of FeSODs from the trypanosomes T. brucei at 2.0 angstrom and T. cruzi at 1.9 angstrom resolution, and that from P. falciparum at a higher resolution (2.0 angstrom) to that previously reported. The homodimeric enzymes are compared to the related human MnSOD with particular attention to structural aspects which are relevant for drug design. Although the structures possess a very similar overall fold, differences between the enzymes at the entrance to the channel which leads to the active site could be identified. These lead to a slightly broader and more positively charged cavity in the parasite enzymes. Furthermore, a statistical coupling analysis (SCA) for the whole Fe/MnSOD family reveals different patterns of residue coupling for Mn and Fe SODs, as well as for the dimeric and tetrameric states. In both cases, the statistically coupled residues lie adjacent to the conserved core surrounding the metal center and may be expected to be responsible for its fine tuning, leading to metal ion specificity.

Molecular characterization of Arabidopsis thaliana PUF proteins - binding specificity and target candidates

PUF proteins regulate both stability and translation through sequence-specific binding to the 3` UTR of target mRNA transcripts. Binding is mediated by a conserved PUF domain, which contains eight repeats of approximately 36 amino acids each. Found in all eukaryotes, they have been related to several developmental processes. Analysis of the 25 Arabidopsis Pumilio (APUM) proteins presenting PUF repeats reveals that 12 (APUM-1 to APUM-12) have a PUF domain with 50-75% similarity to the Drosophila PUF domain. Through three-hybrid assays, we show that APUM-1 to APUM-6 can bind specifically to the Nanos response element sequence recognized by Drosophila Pumilio. Using an Arabidopsis RNA library in a three-hybrid screening, we were able to identify an APUM-binding consensus sequence. Computational analysis allowed us to identify the APUM-binding element within the 3` UTR in many Arabidopsis transcripts, even in important mRNAs related to shoot stem cell maintenance. We demonstrate that APUM-1 to APUM-6 are able to bind specifically to APUM-binding elements in the 3` UTR of WUSCHEL, CLAVATA-1, PINHEAD/ZWILLE and FASCIATA-2 transcripts. The results obtained in the present study indicate that the APUM proteins may act as regulators in Arabidopsis through an evolutionarily conserved mechanism, which may open up a new approach for investigating mRNA regulation in plants.

Crystallographic structure and substrate-binding interactions of the molybdate-binding protein of the phytopathogen Xanthomonas axonopodis pv. citri

In Xanthomonas axonopodis pv. citri (Xac or X citri), the modA gene codes for a periplasmic protein (ModA) that is capable of binding molybdate and tungstate as part of the ABC-type transporter required for the uptake of micronutrients. In this study, we report the crystallographic structure of the Xac ModA protein with bound molybdate. The Xac ModA structure is similar to orthologs with known three-dimensional structures and consists of two nearly symmetrical domains separated by a hinge region where the oxyanion-binding site lies. Phylogenetic analysis of different ModA orthologs based on sequence alignments revealed three groups of molybdate-binding proteins: bacterial phytopathogens, enterobacteria and soil bacteria. Even though the ModA orthologs are segregated into different groups, the ligand-binding hydrogen bonds are mostly conserved, except for Archaeglobus fulgidus ModA. A detailed discussion of hydrophobic interactions in the active site is presented and two new residues, Ala(38) and Ser(151), are shown to be part of the ligand-binding pocket. (c) 2007 Elsevier B.V All rights reserved.

INSECT CHYMOTRYPSINS: CHLOROMETHYL KETONE INACTIVATION AND SUBSTRATE SPECIFICITY RELATIVE TO POSSIBLE COEVOLUTIONAL ADAPTATION OF INSECTS AND PLANTS

Insect digestive chymotrypsins are present in a large variety of insect orders but their substrate specificity still remains unclear. Ewer insect chymotrypsins from 3 different insect orders (Dictyoptera, Coleoptera and two Lepidoptera) were isolated using affinity chromatography. Enzymes presented molecular masses in the range of 20 to 31 kDa and pH optima in the range of 7.5 to 10.0. Kinetic characterization. using different, colorimetric and fluorescent substrates indicated that insect chymotrypsins differ from, bovine chymotrypsin in their primary specificity toward small substrates (like N-benzoyl-L-Tyr p-nitroanilide) rather than on their preference for large substrates (exemplified by Succynil-Ala-Ala-Pro-Phe P-nitroanilide). Chloromethyl ketones (TPCK, N-alpha-tosyl-L-Phe chloromethyl ketone and Z-GGF-CK, N-carbobenzoxy-Gly-Gly-phe-CK) inactivated all chymotrypsins legated. Inactivation rates follow apparent first-order kinetics with variable second order rates (TPCK, 42 to 130 M(-1)s(-1); Z-GGF-CK, 150 to 450 M(-1)s(-1) that may be remarkably low for S. frugiperda chymotrypsin (TPCK, 6 M(-1)s(-1); Z-GGF-CK, 6.1 M(-1) s(-1)). Homology modelling and sequence alignment showed that. in lepidopteran chymotrypsins, differences in the amino acid residues in the neighborhood of the catalytic His 57 may affect its pKa, value. This is Proposed as the cause of the decrease in His 57 reactivity toward chloromethyl ketones. Such amino acid replacement in the active site is proposed. to be an adaptation to the presence of dietary ketones. (C) 2009 Wiley Periodicals, Inc.

Subsite substrate specificity of midgut insect chymotrypsins

Insect chymotrypsins are distinctively sensitive to plant protein inhibitors, suggesting that they differ in subsite architecture and hence in substrate specificities. Purified digestive chymotrypsins from insects of three different orders were assayed with internally quenched fluorescent oligopeptides with three different amino acids at P1 (Tyr, Phe, and Leu) and 13 amino acid replacements in positions P1`, P2, and P3. The binding energy (Delta G(s), calculated from Km values) and the activation energy (Delta G(T)(double dagger), determined from k(cat)/K-m values) were calculated. The hydrophobicities of each subsite were calculated from the efficiency of hydrolysis of the different amino acid replacements at that subsite. The results showed that except for S1, the other subsites (S2, S3, and S1`) vary among chymotrypsins. This result contrasts with insect trypsin data that revealed a trend along evolution, putatively associated with resistance to plant inhibitors. In spite of those differences, the data suggested that in lepidopteran chymotrypsins S2 and S1` bind the substrate ground state, whereas only S1` binds the transition state, supporting aspects of the present accepted mechanism of catalysis. 2008 Elsevier Ltd. All rights reserved.

Catalytic properties of thimet oligopeptidase H600A mutant

Thimet oligopeptidase (EC 3.4.24.15, TOP) is a metallo-oligopeptidase that participates in the intracellular metabolism of peptides. Predictions based on structurally analogous peptidases (Dcp and ACE-2) show that TOP can present a hinge-bend movement during substrate hydrolysis, what brings some residues closer to the substrate. One of these residues that in TOP crystallographic structure are far from the catalytic residues, but, moves toward the substrate considering this possible structural reorganization is His(600). In the present work, the role of His(600) of TOP was investigated by site-directed mutagenesis. TOP H600A mutant was characterized through analysis of S(1) and S(1)`, specificity, pH-activity profile and inhibition by JA-2. Results showed that TOP His(600) residue makes important interactions with the substrate, supporting the prediction that His(600) moves toward the substrate due to a hinge movement similar to the Dcp and ACE-2. Furthermore, the mutation H600A affected both K(m) and k(cat), showing the importance of His(600) for both substrate binding and/or product release from active site. Changes in the pH-profile may indicate also the participation of His(600) in TOP catalysis, transferring a proton to the newly generated NH(2)-terminus or helping Tyr(605) and/or Tyr(612) in the intermediate oxyanion stabilization. (C) 2010 Elsevier Inc. All rights reserved.

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.

Sequence and function of lysosomal and digestive cathepsin D-like proteinases of Musca domestica midgut

Musca domestica larvae display in anterior and middle midgut contents, a proteolytic activity with pH optimum of 3.0-3.5 and kinetic properties like cathepsin D. Three cDNAs coding for preprocathepsin D-like proteinases (ppCAD 1, ppCAD 2, ppCAD 3) were cloned from a M. domestica midgut cDNA library. The coded protein sequences included the signal peptide, propeptide and mature enzyme that has all conserved catalytic and substrate binding residues found in bovine lysosomal cathepsin D. Nevertheless, ppCAD 2 and ppCAD 3 lack the characteristic proline loop and glycosylation sites. A comparison among the sequences of cathepsin D-like enzymes from some vertebrates and those found in M. domestica and in the genomes of Aedes aegypti, Drosophila melanogaster, Tribolium castaneum, and Bombyx mori showed that only flies have enzymes lacking the proline loop (as defined by the motif: DxPxPx(G/A)P), thus resembling vertebrate pepsin. ppCAD 3 should correspond to the digestive cathepsin D-like proteinase (CAD) found in enzyme assays because: (1) it seems to be the most expressed CAD, based on the frequency of ESTs found. (2) The mRNA for CAD 3 is expressed only in the anterior and proximal middle midgut. (3) Recombinant procathepsin D-like proteinase (pCAD 3), after auto-activation has a pH optimum of 2.5-3.0 that is close to the luminal pH of M. domestica midgut. (4) Immunoblots of proteins from different tissues revealed with anti-pCAD 3 serum were positive only in samples of anterior and middle midgut tissue and contents. (5) CAD 3 is localized with immunogold inside secretory vesicles and around microvilli in anterior and middle midguit cells. The data support the view that on adapting to deal with a bacteria-rich food in an acid midgut region, M. domestica digestive CAD resulted from the same archetypical gene as the intracellular cathepsin D, paralleling what happened with vertebrates. The lack of the proline loop may be somehow associated with the extracellular role of both pepsin and digestive CAD 3. (C) 2009 Elsevier Ltd. All rights reserved.

Structural Aspects of the Distinct Biochemical Properties of Glutaredoxin 1 and Glutaredoxin 2 from Saccharomyces cerevisiae

Glutaredoxins (Grxs) are small (9-12 kDa) heat-stable proteins that are ubiquitously distributed. In Saccharomyces cerevisiae, seven Grx enzymes have been identified. Two of them (yGrx1 and yGrx2) are dithiolic, possessing a conserved Cys-Pro-Tyr-Cys motif. Here, we show that yGrx2 has a specific activity 15 times higher than that of yGrx1, although these two oxidoreductases share 64% identity and 85% similarity with respect to their amino acid sequences. Further characterization of the enzymatic activities through two-substrate kinetics analysis revealed that yGrx2 possesses a lower Km for glutathione and a higher turnover than yGrx1. To better comprehend these biochemical differences, the pK(a) of the N-terminal active-site cysteines (Cys27) of these two proteins and of the yGrx2-C30S mutant were determined. Since the pK(a) values of the yGrx1 and yGix2 Cys27 residues are very similar, these parameters cannot account for the difference observed between their specific activities. Therefore, crystal structures of yGrx2 in the oxidized form and with a glutathionyl mixed disulfide were determined at resolutions of 2.05 and 1.91 angstrom, respectively. Comparisons of yGrx2 structures with the recently determined structures of yGrx1 provided insights into their remarkable functional divergence. We hypothesize that the substitutions of Ser23 and Gln52 in yGrx1 by Ala23 and Glu52 in yGrx2 modify the capability of the active-site C-terminal cysteine to attack the mixed disulfide between the N-terminal active-site cysteine and the glutathione molecule. Mutagenesis studies supported this hypothesis. The observed structural and functional differences between yGrx1 and yGrx2 may reflect variations in substrate specificity. (C) 2008 Elsevier Ltd. All rights reserved.

Quantitative structure-activity relationships for a series of inhibitors of cruzain from Trypanosoma cruzi: Molecular modeling, CoMFA and CoMSIA studies

Human parasitic diseases are the foremost threat to human health and welfare around the world. Trypanosomiasis is a very serious infectious disease against which the currently available drugs are limited and not effective. Therefore, there is an urgent need for new chemotherapeutic agents. One attractive drug target is the major cysteine protease from Trypanosoma cruzi, cruzain. In the present work, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) studies were conducted on a series of thiosemicarbazone and semicarbazone derivatives as inhibitors of cruzain. Molecular modeling studies were performed in order to identify the preferred binding mode of the inhibitors into the enzyme active site, and to generate structural alignments for the three-dimensional quantitative structure-activity relationship (3D QSAR) investigations. Statistically significant models were obtained (CoMFA. r(2) = 0.96 and q(2) = 0.78; CoMSIA, r(2) = 0.91 and q(2) = 0.73), indicating their predictive ability for untested compounds. The models were externally validated employing a test set, and the predicted values were in good agreement with the experimental results. The final QSAR models and the information gathered from the 3D CoMFA and CoMSIA contour maps provided important insights into the chemical and structural basis involved in the molecular recognition process of this family of cruzain inhibitors, and should be useful for the design of new structurally related analogs with improved potency. (C) 2009 Elsevier Inc. All rights reserved.

Trypanosoma cruzi: Isolation and characterization of aspartyl proteases

Two aspartyl proteases activities were identified and isolated from Trypanosoma cruzi epimastigotes: cruzipsin-I (CZP-I) and cruzipsin-II (CZP-II). One was isolated from a soluble fraction (CZP-II) and the other was solubilized with 3-[(3-cholamidopropyl)-dimethylammonio]-1-propanesulfonate(CZP-I). The molecular mass of both proteases was estimated to be 120 kDa by HPLC gel filtration and the activity of the enzymes was detected in a doublet of bands (56 and 48 kDa) by substrate-sodium dodecyl sulphate-polyacrylamide-gelatin gel electrophoresis. Substrate specificity studies indicated that the enzymes consistently hydrolyze the cathepsin D substrate Phe-Ala-Ala-Phe (4-NO(2))-Phe-Val-Leu-O(4)MP but failed to hydrolyze serine and other protease substrates. Both proteases activities were strongly inhibited by the classic inhibitor pepstatin-A (>= 68%) and the aspartic active site labeling agent, 1,2-epoxy-3-(phenyl-nitrophenoxy) propane (>= 80%). These findings show that both proteases are novel T. cruzi acidic proteases. The physiological function of these enzymes in T. cruzi has under investigation. (c) 2009 Elsevier Inc. All rights reserved.

Electronic Structure and Spectro-Structural Correlations of Fe(III)Zn(II) Biomimetics for Purple Acid Phosphatases: Relevance to DNA Cleavage and Cytotoxic Activity

Purple acid phosphatases (PAPs) are a group of metallohydrolases that contain a dinuclear Fe(II)M(II) center (M(II) = Fe, Mn, Zn) in the active site and are able to catalyze the hydrolysis of a variety of phosphoric acid esters. The dinuclear complex [(H(2)O)Fe(III)(mu-OH)Zn(II)(L-H)](CIO(4))(2) (2) with the ligand 2-[N-bis(2-pyridylmethyl)aminomethyl]-4-methyl-6-[N-(2-pyridylmethyl)(2-hydroxybenzyl) aminomethyl]phenol (H(2)L-H) has recently been prepared and is found to closely mimic the coordination environment of the Fe(III)Zn(II) active site found in red kidney bean PAP (Neves et al. J. Am. Chem. Soc. 2007, 129, 7486). The biomimetic shows significant catalytic activity in hydrolytic reactions. By using a variety of structural, spectroscopic, and computational techniques the electronic structure of the Fe(III) center of this biomimetic complex was determined. In the solid state the electronic ground state reflects the rhombically distorted Fe(III)N(2)O(4) octahedron with a dominant tetragonal compression align ad along the mu-OH-Fe-O(phenolate) direction. To probe the role of the Fe-O(phenolate) bond, the phenolate moiety was modified to contain electron-donating or -withdrawing groups (-CH(3), -H, -Br, -NO(2)) in the 5-position. Tie effects of the substituents on the electronic properties of the biomimetic complexes were studied with a range of experimental and computational techniques. This study establishes benchmarks against accurate crystallographic struck ral information using spectroscopic techniques that are not restricted to single crystals. Kinetic studies on the hydrolysis reaction revealed that the phosphodiesterase activity increases in the order -NO(2)<- Br <- H <- CH(3) when 2,4-bis(dinitrophenyl)phosphate (2,4-bdnpp) was used as substrate, and a linear free energy relationship is found when log(k(cat)/k(0)) is plotted against the Hammett parameter a. However, nuclease activity measurements in the cleavage of double stranded DNA showed that the complexes containing the electron-withdrawing -NO(2) and electron-donating CH3 groups are the most active while the cytotoxic activity of the biomimetics on leukemia and lung tumoral cells is highest for complexes with electron-donating groups.

Role of cis-cis muconic acid in the catalysis of Pseudomonas putida chlorocatechol 1,2-dioxygenase

Chlorocatechol 1,2-dioxygenase (1,2-CCD) is a non-heme iron protein involved in the intradiol cleavage of aromatic compounds that are recalcitrant to biodegradation. In particular, 1,2-CCD catalyzes the conversion of catechol and its halogenated derivatives to cis-cis muconic acid. In this study we describe a series of experiments concerning the interaction of chlorocatechol 1,2-dioxygenase from Pseudomonas putida (Pp1,2-CCD) with cis-cis muconic acid. We used single-injection ITC to show that the reaction product inhibits enzyme kinetics. DSC and EPR measurements probed whether this was accomplished by a direct binding of the product to the enzyme active site. DSC shows that cis-cis muconic acid affects the thermal unfolding of the protein and allowed us to estimate a binding constant. Furthermore, EPR spectra of the Fe(III) center demonstrate that, upon product binding, a significant decrease in resonance intensity is observed, indicating that cis-cis muconic acid binds directly to the active site. Based on the increasing interest for understanding dioxygenases mechanism of action and, moreover, how to control such process, our data indicate that the product of the reaction does play a relevant role in the catalysis and should therefore be taken into account when one thinks about ways of regulating enzyme activity. (C) 2010 Elsevier B.V. All rights reserved.

Crystal structure of Schistosoma purine nucleoside phosphorylase complexed with a novel monocyclic inhibitor

A novel inhibitor of Schistosoma PNP was identified using an ""in silico"" approach allied to enzyme inhibition assays. The compound has a monocyclic structure which has not been previously described for PNP inhibitors The crystallographic structure of the complex was determined and used to elucidate the binding mode within the active site Furthermore, the predicted pose was very similar to that determined crystallographically, validating the methodology The compound Sm_VS1, despite its low molecular weight, possesses an IC(50) of 1 3 mu M, surprisingly low when compared with purine analogues This is presumably due to the formation of eight hydrogen bonds with key residues in the active site E203, N245 and T244. The results of this study highlight the importance of the use of multiple conformations for the target during virtual screening. Indeed the Sm_VS1 compound was only identified after flipping the N245 side chain It is expected that the structure will be of use in the development of new highly active non-purine based compounds against the Sclustosoma enzyme. (c) 2010 Elsevier B V. All rights reserved

Statistical coupling analysis of aspartic proteinases based on crystal structures of the Trichoderma reesei enzyme and its complex with pepstatin A

The crystal structures of an aspartic proteinase from Trichoderma reesei (TrAsP) and of its complex with a competitive inhibitor, pepstatin A, were solved and refined to crystallographic R-factors of 17.9% (R(free)=21.2%) at 1.70 angstrom resolution and 15.81% (R(free) = 19.2%) at 1.85 angstrom resolution, respectively. The three-dimensional structure of TrAsP is similar to structures of other members of the pepsin-like family of aspartic proteinases. Each molecule is folded in a predominantly beta-sheet bilobal structure with the N-terminal and C-terminal domains of about the same size. Structural comparison of the native structure and the TrAsP-pepstatin complex reveals that the enzyme undergoes an induced-fit, rigid-body movement upon inhibitor binding, with the N-terminal and C-terminal lobes tightly enclosing the inhibitor. Upon recognition and binding of pepstatin A, amino acid residues of the enzyme active site form a number of short hydrogen bonds to the inhibitor that may play an important role in the mechanism of catalysis and inhibition. The structures of TrAsP were used as a template for performing statistical coupling analysis of the aspartic protease family. This approach permitted, for the first time, the identification of a network of structurally linked residues putatively mediating conformational changes relevant to the function of this family of enzymes. Statistical coupling analysis reveals coevolved continuous clusters of amino acid residues that extend from the active site into the hydrophobic cores of each of the two domains and include amino acid residues from the flap regions, highlighting the importance of these parts of the protein for its enzymatic activity. (C) 2008 Elsevier Ltd. All rights reserved.