Thioredoxin and peptide methionine sulfoxide reductase: convergence of similar structure and function in distinct structural folds.

Thioredoxin (Trx) and peptide methionine sulfoxide reductase (PMSR) are small thiol oxidoreductases implicated in antioxidant defense and redox regulation of cellular processes. Here we show that the structures of Trx and PMSR exhibit resemblance in their alphabeta core regions and that the active site cysteines in two proteins occupy equivalent positions downstream of a central beta-strand and at the N-terminus of an alpha-helix. Moreover, we identified a PMSR subfamily that contains an active site CxxC motif (two cysteines separated by two other amino acids) positioned similarly to the catalytic redox active CxxC motif in Trx. However, Trx and PMSR are characterized by distinct ancient folds that differ in both orientation of secondary structures and their patterns. Trx is a member of the Trx-fold superfamily, whereas PMSR has a unique fold not found in other proteins. The data suggest that similar structures and functions of Trx and PMSR were acquired independently during evolution and point to a general strategy of identifying new redox regulatory proteins.     

Topological variation in the evolution of new reactions in functionally diverse enzyme superfamilies

In functionally diverse enzyme superfamilies (SFs), conserved structural and active site features reflect catalytic capabilities 'hard-wired' in each SF architecture. Overlaid on this foundation, evolutionary changes in active site machinery, structural topology and other aspects of structural organization and interactions support the emergence of new reactions, mechanisms, and substrate specificity. This review connects topological with functional variation in each of the haloalkanoic acid dehalogenase (HAD) and vicinal oxygen chelate fold (VOC) SFs and a set of redox-active thioredoxin (Trx)-fold SFs to illustrate a few of the varied themes nature has used to evolve new functions from a limited set of structural scaffolds.     

A conserved cis-proline precludes metal binding by the active site thiolates in members of the thioredoxin family of proteins

Many thioredoxin-fold proteins possess a conserved cis-proline located in their C-terminal portions. This residue, as well as catalytic and resolving cysteines, is a key functional group in the active sites of these thiol-disulfide oxidoreductases. However, the specific function of the proline is poorly understood, and some thioredoxin-fold proteins lack this residue. Herein, we found that mutation of a cis-proline, Pro75, in human thioredoxin to serine, threonine, or alanine leads to the formation of an Fe2-S2 cluster in this protein. Further mutagenesis studies revealed that the first cysteine in the CxxC motif and a cysteine in the C-terminal region of the protein were responsible for metal binding. Replacement of Pro75 with arginine, a residue that occurs in place of Pro in peroxiredoxins, also led to the formation of the cluster in the thioredoxin. In addition, we found that mutation of the TxxC active site in a peroxiredoxin to the CxxC form could lead to coordination of an Fe2-S2 cluster in these proteins in vitro. Sco1, a distantly related thioredoxin-fold protein, has histidine in place of the cis-proline, and this residue binds copper. The Pro75His mutation led to increased copper binding by human thioredoxin when cells were grown in the presence of this trace element. Taken together, our data suggest that an important function of Pro75 in human thioredoxin, and likely other members of this superfamily, is to prevent metal binding by the reactive thiolate-based active site.     

An Extended Active-site Motif Controls the Reactivity of the Thioredoxin Fold

Proteins belonging to the thioredoxin (Trx) superfamily are abundant in all organisms. They share the same structural features, arranged in a seemingly simple fold, but they perform a multitude of functions in oxidative protein folding and electron transfer pathways. We use the C-terminal domain of the unique transmembrane reductant conductor DsbD as a model for an in-depth analysis of the factors controlling the reactivity of the Trx fold. We employ NMR spectroscopy, x-ray crystallography, mutagenesis, in vivo functional experiments applied to DsbD, and a comparative sequence analysis of Trx-fold proteins to determine the effect of residues in the vicinity of the active site on the ionization of the key nucleophilic cysteine of the -CXXC- motif. We show that the function and reactivity of Trx-fold proteins depend critically on the electrostatic features imposed by an extended active-site motif.     

Evidence for proton shuffling in a thioredoxin-like protein during catalysis

Proteins of the thioredoxin (Trx) superfamily catalyze disulfide-bond formation, reduction and isomerization in substrate proteins both in prokaryotic and in eukaryotic cells. All members of the Trx family with thiol-disulfide oxidoreductase activity contain the characteristic Cys-X-X-Cys motif in their active site. Here, using Poisson-Boltzmann-based protonation-state calculations based on 100-ns molecular dynamics simulations, we investigate the catalytic mechanism of DsbL, the most oxidizing Trx-like protein known to date. We observed several correlated transitions in the protonation states of the buried active-site cysteine and a neighboring lysine coupled to the exposure of the active-site thiolate. These results support the view of an internal proton shuffling mechanism during oxidation crucial for the uptake of two electrons from the substrate protein. Intramolecular disulfide-bond formation is probably steered by the conformational switch facilitating interaction with the active-site thiolate. A consistent catalytic mechanism for DsbL, probably conferrable to other proteins of the same class, is presented. Our results suggest a functional role of hydration entropy of active-site groups.     

Physical and genetic interactions of yeast Cwc21p,an ortholog of human SRm300/SRRM2, suggest a role at the catalytic center of the spliceosome

In Saccharomyces cerevisiae, Cwc21p is a protein of unknown function that is associated with the NineTeen Complex (NTC), a group of proteins involved in activating the spliceosome to promote the pre-mRNA splicing reaction. Here, we show that Cwc21p binds directly to two key splicing factors—namely, Prp8p and Snu114p—and becomes the first NTC-related protein known to dock directly to U5 snRNP proteins. Using a combination of proteomic techniques we show that the N-terminus of Prp8p contains an intramolecular fold that is a Snu114p and Cwc21p interacting domain (SCwid). Cwc21p also binds directly to the C-terminus of Snu114p. Complementary chemical cross-linking experiments reveal reciprocal protein footprints between the interacting Prp8 and Cwc21 proteins, identifying the conserved cwf21 domain in Cwc21p as a Prp8p binding site. Genetic and functional interactions between Cwc21p and Isy1p indicate that they have related functions at or prior to the first catalytic step of splicing, and suggest that Cwc21p functions at the catalytic center of the spliceosome, possibly in response to environmental or metabolic changes. We demonstrate that SRm300, the only SR-related protein known to be at the core of human catalytic spliceosomes, is a functional ortholog of Cwc21p, also interacting directly with Prp8p and Snu114p. Thus, the function of Cwc21p is likely conserved from yeast to humans.     

Structure-Based Phylogeny as a Diagnostic for Functional Characterization of Proteins with a Cupin Fold

Background

The members of cupin superfamily exhibit large variations in their sequences, functions, organization of domains, quaternary associations and the nature of bound metal ion, despite having a conserved β-barrel structural scaffold. Here, an attempt has been made to understand structure-function relationships among the members of this diverse superfamily and identify the principles governing functional diversity. The cupin superfamily also contains proteins for which the structures are available through world-wide structural genomics initiatives but characterized as “hypothetical”. We have explored the feasibility of obtaining clues to functions of such proteins by means of comparative analysis with cupins of known structure and function.

Methodology/Principal Findings

A 3-D structure-based phylogenetic approach was undertaken. Interestingly, a dendrogram generated solely on the basis of structural dissimilarity measure at the level of domain folds was found to cluster functionally similar members. This clustering also reflects an independent evolution of the two domains in bicupins. Close examination of structural superposition of members across various functional clusters reveals structural variations in regions that not only form the active site pocket but are also involved in interaction with another domain in the same polypeptide or in the oligomer.

Conclusions/Significance

Structure-based phylogeny of cupins can influence identification of functions of proteins of yet unknown function with cupin fold. This approach can be extended to other proteins with a common fold that show high evolutionary divergence. This approach is expected to have an influence on the function annotation in structural genomics initiatives.     

Anatomy of the E2 ligase fold: implications for enzymology and evolution of ubiquitin/Ub-like protein conjugation

The configuration of the active site of E2 ligases, central enzymes in the ubiquitin/ubiquitin-like protein (Ub/Ubl) conjugation systems, has long puzzled researchers. Taking advantage of the wealth of newly available structures and sequences of E2s from diverse organisms, we performed a large-scale comparative analysis of these proteins. As a result we identified a previously under-appreciated diversity in the active site of these enzymes, in particular, the spatial location of the catalytic cysteine and a constellation of associated conserved residues that potentially contributes to catalysis. We observed structural innovations of differing magnitudes occurring in various families across the E2 fold that might correlate in part with differences in target interaction. A key finding was the independent emergence on multiple occasions of a polar residue, often a histidine, in the vicinity of the catalytic cysteine in different E2 families. We propose that these convergently emerging polar residues have a common function, such as in the stabilization of oxyanion holes during Ub/Ubl transfer and spatial localization of the Ub/Ubl tails in the active site. Thus, the E2 ligases represent a rare example in enzyme evolution of high structural diversity of the active site and position of the catalytic residue despite all characterized members catalyzing a similar reaction. Our studies also indicated certain evolutionarily conserved features in all active members of the E2 superfamily that stabilize the unusual flap-like structure in the fold. These features are likely to form a critical mechanical element of the fold required for catalysis. The results presented here could aid in new experiments to understand E2 catalysis.     

The RNase H-like superfamily: new members,comparative structural analysis and evolutionary classification

Ribonuclease H-like (RNHL) superfamily, also called the retroviral integrase superfamily, groups together numerous enzymes involved in nucleic acid metabolism and implicated in many biological processes, including replication, homologous recombination, DNA repair, transposition and RNA interference. The RNHL superfamily proteins show extensive divergence of sequences and structures. We conducted database searches to identify members of the RNHL superfamily (including those previously unknown), yielding >60 000 unique domain sequences. Our analysis led to the identification of new RNHL superfamily members, such as RRXRR (PF14239), DUF460 (PF04312, COG2433), DUF3010 (PF11215), DUF429 (PF04250 and COG2410, COG4328, COG4923), DUF1092 (PF06485), COG5558, OrfB_IS605 (PF01385, COG0675) and Peptidase_A17 (PF05380). Based on the clustering analysis we grouped all identified RNHL domain sequences into 152 families. Phylogenetic studies revealed relationships between these families, and suggested a possible history of the evolution of RNHL fold and its active site. Our results revealed clear division of the RNHL superfamily into exonucleases and endonucleases. Structural analyses of features characteristic for particular groups revealed a correlation between the orientation of the C-terminal helix with the exonuclease/endonuclease function and the architecture of the active site. Our analysis provides a comprehensive picture of sequence-structure-function relationships in the RNHL superfamily that may guide functional studies of the previously uncharacterized protein families.     

A redox-active crosslinker reveals an essential and inhibitable oxidative folding network in the endoplasmic reticulum of malaria parasites

Malaria remains a major global health problem, creating a constant need for research to identify druggable weaknesses in P. falciparum biology. As important components of cellular redox biology, members of the Thioredoxin (Trx) superfamily of proteins have received interest as potential drug targets in Apicomplexans. However, the function and essentiality of endoplasmic reticulum (ER)-localized Trx-domain proteins within P. falciparum has not been investigated. We generated conditional mutants of the protein PfJ2—an ER chaperone and member of the Trx superfamily—and show that it is essential for asexual parasite survival. Using a crosslinker specific for redox-active cysteines, we identified PfJ2 substrates as PfPDI8 and PfPDI11, both members of the Trx superfamily as well, which suggests a redox-regulatory role for PfJ2. Knockdown of these PDIs in PfJ2 conditional mutants show that PfPDI11 may not be essential. However, PfPDI8 is required for asexual growth and our data suggest it may work in a complex with PfJ2 and other ER chaperones. Finally, we show that the redox interactions between these Trx-domain proteins in the parasite ER and their substrates are sensitive to small molecule inhibition. Together these data build a model for how Trx-domain proteins in the P. falciparum ER work together to assist protein folding and demonstrate the suitability of ER-localized Trx-domain proteins for antimalarial drug development.     

A genome-wide survey of human thioredoxin and glutaredoxin family pseudogenes

The thioredoxin/glutaredoxin family consists of small heat-stable proteins that have a highly conserved CXXC active site and that participate in the regulation of many redox reactions. We have searched the human genome sequence to find putative pseudogenes (non-functional copies of protein-coding genes) for all known members of this family. This survey has resulted in the identification of seven processed pseudogenes for human Trx1 and two more for human Grx1. No evidence for the presence of processed pseudogenes has been found for the remaining members of this family. In addition, we have been unable to detect any non-processed pseudogenes derived from any member of the family in the human genome. The seven thioredoxin pseudogenes can be divided into two groups: Trx1-psi2, -psi3, -psi4, -psi5 and -psi6 arose from the functional ancestor, whereas Trx1-psi1 and -psi7 originated from Trx1-psi2 and -psi6, respectively. In all cases, the pseudogenes originated after the human/rodent radiation as shown by phylogenetic analysis. This is also the case for Grx1-psi1 and Grx1-psi2, which are placed between rodent and human sequences in the phylogenetic tree. Our study provides a molecular record of the recent evolution of these two genes in the hominid lineage.     

Crystal structures of two functionally different thioredoxins in spinach chloroplasts

Thioredoxins are small ubiquitous proteins which act as general protein disulfide reductases in living cells. Chloroplasts contain two distinct thioredoxins ( f and m) with different phylogenetic origin. Both act as enzyme regulatory proteins but have different specificities towards target enzymes. Thioredoxin f (Trx f), which shares only low sequence identity with thioredoxin m (Trx m) and with all other known thioredoxins, activates enzymes of the Calvin cycle and other photosynthetic processes. Trx m shows high sequence similarity with bacterial thioredoxins and activates other chloroplast enzymes. The here described structural studies of the two chloroplast thioredoxins were carried out in order to gain insight into the structure/function relationships of these proteins. Crystal structures were determined for oxidized, recombinant thioredoxin f (Trx f-L) and at the N terminus truncated form of it (Trx f-S), as well as for oxidized and reduced thioredoxin m (at 2.1 and 2.3 A resolution, respectively). Whereas thioredoxin f crystallized as a monomer, both truncated thioredoxin f and thioredoxin m crystallized as non-covalent dimers. The structures of thioredoxins f and m exhibit the typical thioredoxin fold consisting of a central twisted five-stranded beta-sheet surrounded by four alpha-helices. Thioredoxin f contains an additional alpha-helix at the N terminus and an exposed third cysteine close to the active site. The overall three-dimensional structures of the two chloroplast thioredoxins are quite similar. However, the two proteins have a significantly different surface topology and charge distribution around the active site. An interesting feature which might significantly contribute to the specificity of thioredoxin f is an inherent flexibility of its active site, which has expressed itself crystallographically in two different crystal forms.     

A Threading-Based Method for the Prediction of DNA-Binding Proteins with Application to the Human Genome

Diverse mechanisms for DNA-protein recognition have been elucidated in numerous atomic complex structures from various protein families. These structural data provide an invaluable knowledge base not only for understanding DNA-protein interactions, but also for developing specialized methods that predict the DNA-binding function from protein structure. While such methods are useful, a major limitation is that they require an experimental structure of the target as input. To overcome this obstacle, we develop a threading-based method, DNA-Binding-Domain-Threader (DBD-Threader), for the prediction of DNA-binding domains and associated DNA-binding protein residues. Our method, which uses a template library composed of DNA-protein complex structures, requires only the target protein''s sequence. In our approach, fold similarity and DNA-binding propensity are employed as two functional discriminating properties. In benchmark tests on 179 DNA-binding and 3,797 non-DNA-binding proteins, using templates whose sequence identity is less than 30% to the target, DBD-Threader achieves a sensitivity/precision of 56%/86%. This performance is considerably better than the standard sequence comparison method PSI-BLAST and is comparable to DBD-Hunter, which requires an experimental structure as input. Moreover, for over 70% of predicted DNA-binding domains, the backbone Root Mean Square Deviations (RMSDs) of the top-ranked structural models are within 6.5 Å of their experimental structures, with their associated DNA-binding sites identified at satisfactory accuracy. Additionally, DBD-Threader correctly assigned the SCOP superfamily for most predicted domains. To demonstrate that DBD-Threader is useful for automatic function annotation on a large-scale, DBD-Threader was applied to 18,631 protein sequences from the human genome; 1,654 proteins are predicted to have DNA-binding function. Comparison with existing Gene Ontology (GO) annotations suggests that ∼30% of our predictions are new. Finally, we present some interesting predictions in detail. In particular, it is estimated that ∼20% of classic zinc finger domains play a functional role not related to direct DNA-binding.     

A new family of type III polyketide synthases in Mycobacterium tuberculosis

The Mycobacterium tuberculosis genome has revealed a remarkable array of polyketide synthases (PKSs); however, no polyketide product has been isolated thus far. Most of the PKS genes have been implicated in the biosynthesis of complex lipids. We report here the characterization of two novel type III PKSs from M. tuberculosis that are involved in the biosynthesis of long-chain alpha-pyrones. Measurement of steady-state kinetic parameters demonstrated that the catalytic efficiency of PKS18 protein was severalfold higher for long-chain acyl-coenzyme A substrates as compared with the small-chain precursors. The specificity of PKS18 and PKS11 proteins toward long-chain aliphatic acyl-coenzyme A (C12 to C20) substrates is unprecedented in the chalcone synthase (CHS) family of condensing enzymes. Based on comparative modeling studies, we propose that these proteins might have evolved by fusing the catalytic machinery of CHS and beta-ketoacyl synthases, the two evolutionarily related members with conserved thiolase fold. The mechanistic and structural importance of several active site residues, as predicted by our structural model, was investigated by performing site-directed mutagenesis. The functional identification of diverse catalytic activity in mycobacterial type III PKSs provide a fascinating example of metabolite divergence in CHS-like proteins.     

The Crystal Structure of Arabidopsis VSP1 Reveals the Plant Class C-Like Phosphatase Structure of the DDDD Superfamily of Phosphohydrolases

Arabidopsis thaliana vegetative storage proteins, VSP1 and VSP2, are acid phosphatases and belong to the haloacid dehalogenase (HAD) superfamily. In addition to their potential nutrient storage function, they were thought to be involved in plant defense and flower development. To gain insights into the architecture of the protein and obtain clues about its function, we have tested their substrate specificity and solved the structure of VSP1. The acid phosphatase activities of these two enzymes require divalent metal such as magnesium ion. Conversely, the activity of these two enzymes is inhibited by vanadate and molybdate, but is resistant to inorganic phosphate. Both VSP1 and VSP2 did not exhibit remarkable activities to any physiological substrates tested. In the current study, we presented the crystal structure of recombinant VSP1 at 1.8 Å resolution via the selenomethionine single-wavelength anomalous diffraction (SAD). Specifically, an α-helical cap domain on the top of the α/β core domain is found to be involved in dimerization. In addition, despite of the low sequence similarity between VSP1 and other HAD enzymes, the core domain of VSP1 containing conserved active site and catalytic machinery displays a classic haloacid dehalogenase fold. Furthermore, we found that VSP1 is distinguished from bacterial class C acid phosphatase P4 by several structural features. To our knowledge, this is the first study to reveal the crystal structure of plant vegetative storage proteins.     

Discarding functional residues from the substitution table improves predictions of active sites within three-dimensional structures

Substitutions of individual amino acids in proteins may be under very different evolutionary restraints depending on their structural and functional roles. The Environment Specific Substitution Table (ESST) describes the pattern of substitutions in terms of amino acid location within elements of secondary structure, solvent accessibility, and the existence of hydrogen bonds between side chains and neighbouring amino acid residues. Clearly amino acids that have very different local environments in their functional state compared to those in the protein analysed will give rise to inconsistencies in the calculation of amino acid substitution tables. Here, we describe how the calculation of ESSTs can be improved by discarding the functional residues from the calculation of substitution tables. Four categories of functions are examined in this study: protein–protein interactions, protein–nucleic acid interactions, protein–ligand interactions, and catalytic activity of enzymes. Their contributions to residue conservation are measured and investigated. We test our new ESSTs using the program CRESCENDO, designed to predict functional residues by exploiting knowledge of amino acid substitutions, and compare the benchmark results with proteins whose functions have been defined experimentally. The new methodology increases the Z-score by 98% at the active site residues and finds 16% more active sites compared with the old ESST. We also find that discarding amino acids responsible for protein–protein interactions helps in the prediction of those residues although they are not as conserved as the residues of active sites. Our methodology can make the substitution tables better reflect and describe the substitution patterns of amino acids that are under structural restraints only.     

Interaction Signatures Stabilizing the NAD(P)-Binding Rossmann Fold: A Structure Network Approach

The fidelity of the folding pathways being encoded in the amino acid sequence is met with challenge in instances where proteins with no sequence homology, performing different functions and no apparent evolutionary linkage, adopt a similar fold. The problem stated otherwise is that a limited fold space is available to a repertoire of diverse sequences. The key question is what factors lead to the formation of a fold from diverse sequences. Here, with the NAD(P)-binding Rossmann fold domains as a case study and using the concepts of network theory, we have unveiled the consensus structural features that drive the formation of this fold. We have proposed a graph theoretic formalism to capture the structural details in terms of the conserved atomic interactions in global milieu, and hence extract the essential topological features from diverse sequences. A unified mathematical representation of the different structures together with a judicious concoction of several network parameters enabled us to probe into the structural features driving the adoption of the NAD(P)-binding Rossmann fold. The atomic interactions at key positions seem to be better conserved in proteins, as compared to the residues participating in these interactions. We propose a “spatial motif” and several “fold specific hot spots” that form the signature structural blueprints of the NAD(P)-binding Rossmann fold domain. Excellent agreement of our data with previous experimental and theoretical studies validates the robustness and validity of the approach. Additionally, comparison of our results with statistical coupling analysis (SCA) provides further support. The methodology proposed here is general and can be applied to similar problems of interest.     

Alleviation of oxidative stress by potent and selective thioredoxin-mimetic peptides

One of the major enzymatic cell defenses providing protection from oxidative injury is the TrxR-Trx system. It consists of NADPH and thioredoxin reductase (TrxR), which maintain thioredoxin (Trx) in a reduced state. Perturbing the TrxR-Trx system with the selective TrxR inhibitor auranofin (AuF; 2,3,4,6-tetra-O-acetyl-1-thio-β-D-glucopyranosato-S-(triethylphosphine) gold) induces oxidative stress by keeping Trx in its oxidized state. We have prepared a family of tri- and tetra-oligopeptides derived from the canonical CxxC motif of the Trx active site and a modified CxC motif. These Trx-mimetic compounds are N- and C-terminal-blocked peptides that consist of two cysteine residues that flank the two-amino-acid CxxC motif (CB4 and CB6) or the single-amino-acid CxC motif (CB3). Catecholamine (CA) secretion in bovine chromaffin cells, which is a highly redox sensitive process, is abolished by AuF. The Trx-mimetic peptides effectively restore CA secretion, as monitored by amperometry in single cells. They also prevent the AuF-induced phosphorylation of p38 mitogen-activated protein kinase (MAPK) and c-Jun NH2-terminal kinase. In PC12 cells, the alleviation of AuF-induced ERK1/2-MAPK phosphorylation by Trx-like peptides parallels their effect of restoring CA secretion. CB3, CB4, and CB6 act intracellularly and are significantly more potent than the traditional antioxidants NAC, GSH, DTT, AD4 (NAC-amide), and ascorbic acid. Taken together, the CxxC and CxC peptides represent a new family of potent and selective redox compounds that could serve as potential candidates for prevention and treatment of oxidative-stress-related disorders.     

Protein molecular function prediction by Bayesian phylogenomics

We present a statistical graphical model to infer specific molecular function for unannotated protein sequences using homology. Based on phylogenomic principles, SIFTER (Statistical Inference of Function Through Evolutionary Relationships) accurately predicts molecular function for members of a protein family given a reconciled phylogeny and available function annotations, even when the data are sparse or noisy. Our method produced specific and consistent molecular function predictions across 100 Pfam families in comparison to the Gene Ontology annotation database, BLAST, GOtcha, and Orthostrapper. We performed a more detailed exploration of functional predictions on the adenosine-5′-monophosphate/adenosine deaminase family and the lactate/malate dehydrogenase family, in the former case comparing the predictions against a gold standard set of published functional characterizations. Given function annotations for 3% of the proteins in the deaminase family, SIFTER achieves 96% accuracy in predicting molecular function for experimentally characterized proteins as reported in the literature. The accuracy of SIFTER on this dataset is a significant improvement over other currently available methods such as BLAST (75%), GeneQuiz (64%), GOtcha (89%), and Orthostrapper (11%). We also experimentally characterized the adenosine deaminase from Plasmodium falciparum, confirming SIFTER's prediction. The results illustrate the predictive power of exploiting a statistical model of function evolution in phylogenomic problems. A software implementation of SIFTER is available from the authors.