Glycogen phosphorylase inhibitors: a free energy perturbation analysis of glucopyranose spirohydantoin analogues

GP catalyzes the phosphorylation of glycogen to Glc-1-P. Because of its fundamental role in the metabolism of glycogen, GP has been the target for a systematic structure-assisted design of inhibitory compounds, which could be of value in the therapeutic treatment of type 2 diabetes mellitus. The most potent catalytic-site inhibitor of GP identified to date is spirohydantoin of glucopyranose (hydan). In this work, we employ MD free energy simulations to calculate the relative binding affinities for GP of hydan and two spirohydantoin analogues, methyl-hydan and n-hydan, in which a hydrogen atom is replaced by a methyl- or amino group, respectively. The results are compared with the experimental relative affinities of these ligands, estimated by kinetic measurements of the ligand inhibition constants. The calculated binding affinity for methyl-hydan (relative to hydan) is 3.75 +/- 1.4 kcal/mol, in excellent agreement with the experimental value (3.6 +/- 0.2 kcal/mol). For n-hydan, the calculated value is 1.0 +/- 1.1 kcal/mol, somewhat smaller than the experimental result (2.3 +/- 0.1 kcal/mol). A free energy decomposition analysis shows that hydan makes optimum interactions with protein residues and specific water molecules in the catalytic site. In the other two ligands, structural perturbations of the active site by the additional methyl- or amino group reduce the corresponding binding affinities. The computed binding free energies are sensitive to the preference of a specific water molecule for two well-defined positions in the catalytic site. The behavior of this water is analyzed in detail, and the free energy profile for the translocation of the water between the two positions is evaluated. The results provide insights into the role of water molecules in modulating ligand binding affinities. A comparison of the interactions between a set of ligands and their surrounding groups in X-ray structures is often used in the interpretation of binding free energy differences and in guiding the design of new ligands. For the systems in this work, such an approach fails to estimate the order of relative binding strengths, in contrast to the rigorous free energy treatment.

Cupins: the most functionally diverse protein superfamily?

The cupin superfamily of proteins, named on the basis of a conserved β-barrel fold (‘cupa’ is the Latin term for a small barrel), was originally discovered using a conserved motif found within germin and germin-like proteins from higher plants. Previous analysis of cupins had identified some 18 different functional classes that range from single-domain bacterial enzymes such as isomerases and epimerases involved in the modification of cell wall carbohydrates, through to two-domain bicupins such as the desiccation-tolerant seed storage globulins, and multidomain transcription factors including one linked to the nodulation response in legumes. Recent advances in comparative genomics, and the resolution of many more 3-D structures have now revealed that the largest subset of the cupin superfamily is the 2-oxyglutarate-Fe2+ dependent dioxygenases. The substrates for this subclass of enzyme are many and varied and in total amount to probably 50–100 different biochemical reactions, including several involved in plant growth and development. Although the majority of enzymatic cupins contain iron as an active site metal, other members contain either copper, zinc, cobalt, nickel or manganese ions as a cofactor, with each cofactor allowing a different type of chemistry to occur within the conserved tertiary structure. This review discusses the range of structures and functions found in this most diverse of superfamilies.

Proteomic analysis of arginine adducts on glyoxal-modified ribonuclease

Accumulation of advanced glycation end-products (AGEs) on proteins is associated with the development of diabetic complications. Although the overall extent of modification of protein by AGEs is limited, localization of these modifications at a few critical sites might have a significant effect on protein structure and function. In the present study, we describe the sites of modification of RNase by glyoxal under physiological conditions. Arg(39) and Arg(85), which are closest to the active site of the enzyme, were identified as the primary sites of formation of the glyoxal-derived dihydroxyimidazolidine and hydroimidazolone adducts. Lower amounts of modification were detected at Arg(10), while Arg(33) appeared to be unmodified. We conclude that dihydroxyimidazolidine adducts are the primary products of modification of protein by glyoxal, that Arg(39) and Arg(85) are the primary sites of modification of RNase by glyoxal, and that modification of arginine residues during Maillard reactions of proteins is a highly selective process.

Substrate specificity of human glutamine transaminase K as an aminotransferase and as a cysteine S-conjugate beta-lyase

Rat kidney glutamine transaminase K (GTK) exhibits broad specificity both as an aminotransferase and as a cysteine S-conjugate beta-lyase. The beta-lyase reaction products are pyruvate, ammonium and a sulfhydryl-containing fragment. We show here that recombinant human GTK (rhGTK) also exhibits broad specificity both as an aminotransferase and as a cysteine S-conjugate beta-lyase. S-(1,1,2,2-Tetrafluoroethyl)-L-CySteine is an excellent aminotransferase and beta-lyase substrate of rhGTK. Moderate aminotransferase and beta-lyase activities occur with the chemopreventive agent Se-methyl-L-selenocysteine. L-3-(2-Naphthyl)alanine, L-3-(1-naphthyl)alanine, 5-S-L-cysteinyldopamine and 5-S-L-cysteinyl-L-DOPA are measurable aminotransferase substrates, indicating that the active site can accommodate large aromatic amino acids. The alpha-keto acids generated by transamination/L-amino acid oxidase activity of the two catechol cysteine S-conjugates are unstable. A slow rhGTK-catalyzed beta-elimination reaction, as measured by pyruvate formation, was demonstrated with 5-S-L-CysteinyIdopamine, but not with 5-S-L-CySteinyl-L-DOPA. The importance of transamination, oxidation and beta-elimination reactions involving 5-S-L-cysteinyldopamine, 5-S-L-cysteinyt-L-DOPA and Se-methyl-L-selenocysteirte in human tissues and their biological relevance are discussed. (C) 2008 Elsevier Inc. All rights reserved.

Twists and turns: a tale of two shikimate-pathway enzymes

We are studying two enzymes from the shikimate pathway, dehydroquinate synthase (DHQS) and 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS). Both enzymes have been the subject of numerous studies to elucidate their reaction mechanisms. Crystal structures of DHQS and EPSPS in the presence and absence of substrates, cofactors and/or inhibitors are now available. These structures reveal movements of domains, rearrangements of loops and changes in side-chain positions necessary for the formation of a catalytically competent active site. The potential for using complementary small-angle X-ray scattering (SAXS) studies to confirm the presence of these structural differences in solution has also been explored. Comparative analysis of crystal structures, in the presence and absence of ligands, has revealed structural features critical for substrate-binding and catalysis. We have also analysed these structures by generating GRID energy maps to detect favourable binding sites. The combination of X-ray crystallography, SAXS and computational techniques provides an enhanced analysis of structural features important for the function of these complex enzymes.

Anisotropic refinement of the structure of Thermoascus aurantiacusxylanase I

The isotropic crystallographic model of the structure of xylanase I from Thermoascus aurantiacus (TAXI) has now been refined anisotropically at 1.14 Å resolution to a standard residual of R = 11.1% for all data. TAXI is amongst the five largest proteins deposited in the Protein Data Bank to have been refined with anisotropic displacement parameters (ADPs) at this level of resolution. The anisotropy analysis revealed a more isotropic distribution of anisotropy than usually observed previously. Adding ADPs resulted in high-quality electron-density maps which revealed discrepancies from the previously suggested primary sequences for this enzyme. Side-chain conformational disorder was modelled for 16 residues, including Trp275, a bulky residue at the active site. An unrestrained refinement was consistent with the protonation of the catalytic acid/base glutamate and the deprotonation of the nucleophile glutamate, as required for catalysis. The thermal stability of TAXI is reinterpreted in the light of the new refined model.

The R1441C mutation alters the folding properties of the ROC domain of LRRK2

LRRK2 is a 250 kDa multidomain protein, mutations in which cause familial Parkinson's disease. Previously, we have demonstrated that the R1441C mutation in the ROC domain decreases GTPase activity. Here we show that the R1441C alters the folding properties of the ROC domain, lowering its thermodynamic stability. Similar to small GTPases, binding of different guanosine nucleotides alters the stability of the ROC domain, suggesting that there is an alteration in conformation dependent on GDP or GTP occupying the active site. GTP/GDP bound state also alters the self-interaction of the ROC domain, accentuating the impact of the R1441C mutation on this property. These data suggest a mechanism whereby the R1441C mutation can reduce the GTPase activity of LRRK2, and highlights the possibility of targeting the stability of the ROC domain as a therapeutic avenue in LRRK2 disease.

Effects of oxidative stress on the cardiac myocyte proteome: modifications to peroxiredoxins and small heat shock proteins

Endogenous oxidative stress is a likely cause of cardiac myocyte death in vivo. We examined the early (0-2 h) changes in the proteome of isolated cardiac myocytes from neonatal rats exposed to H2O2 (0.1 mM), focussing on proteins with apparent molecular masses of between 20 and 30 kDa. Proteins were separated by two-dimensional gel electrophoresis (2DGE), located by silver-staining and identified by mass spectrometry. Incorporation of [35S]methionine or 32Pi was also studied. For selected proteins, transcript abundance was examined by reverse transcriptase-polymerase chain reaction. Of the 38 protein spots in the region, 23 were identified. Two families showed changes in 2DGE migration or abundance with H2O2 treatment: the peroxiredoxins and two small heat shock protein (Hsp) family members: heat shock 27 kDa protein 1 (Hsp25) and alphaB-crystallin. Peroxiredoxins shifted to lower pI values and this was probably attributable to 'over-oxidation' of active site Cys-residues. Hsp25 also shifted to lower pI values but this was attributable to phosphorylation. alphaB-crystallin migration was unchanged but its abundance decreased. Transcripts encoding peroxiredoxins 2 and 5 increased significantly. In addition, 10 further proteins were identified. For two (glutathione S-transferase pi, translationally-controlled tumour protein), we could not find any previous references indicating their occurrence in cardiac myocytes. We conclude that exposure of cardiac myocytes to oxidative stress causes post-translational modification in two protein families involved in cytoprotection. These changes may be potentially useful diagnostically. In the short term, oxidative stress causes few detectable changes in global protein abundance as assessed by silver-staining.

Oxidative stress and growth-regulating intracellular signaling pathways in cardiac myocytes

The toxic effects of oxidative stress on cells (including cardiac myocytes, the contractile cells of the heart) are well known. However, an increasing body of evidence has suggested that increased production of reactive oxygen species (ROS) promotes cardiac myocyte growth. Thus, ROS may be 'second messenger' molecules in their own right, and growth-promoting neurohumoral agonists might exert their effects by stimulating production of ROS. The authors review the principal growth-promoting intracellular signaling pathways that are activated by ROS in cardiac myocytes, namely the mitogen-activated protein kinase cascades (extracellular signal-regulated kinases 1/2, c-Jun N-terminal kinases, and p38-mitogen-activated protein kinases) and the phosphoinositide 3-kinase/protein kinase B (Akt) pathway. Possible mechanisms are discussed by which these pathways are activated by ROS, including the oxidation of active site cysteinyl residues of protein and lipid phosphatases with their consequent inactivation, the potential involvement of protein kinase C or the apoptosis signal-regulating kinase 1, and the current models for the activation of the guanine nucleotide binding protein Ras.

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.

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.

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.

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.