Thiol dioxygenases: unique families of cupin proteins

Proteins in the cupin superfamily have a wide range of biological functions in archaea, bacteria and eukaryotes. Although proteins in the cupin superfamily show very low overall sequence similarity, they all contain two short but partially conserved cupin sequence motifs separated by a less conserved intermotif region that varies both in length and amino acid sequence. Furthermore, these proteins all share a common architecture described as a six-stranded β-barrel core, and this canonical cupin or “jelly roll” β-barrel is formed with cupin motif 1, the intermotif region, and cupin motif 2 each forming two of the core six β-strands in the folded protein structure. The recently obtained crystal structures of cysteine dioxygenase (CDO), with contains conserved cupin motifs, show that it has the predicted canonical cupin β-barrel fold. Although there had been no reports of CDO activity in prokaryotes, we identified a number of bacterial cupin proteins of unknown function that share low similarity with mammalian CDO and that conserve many residues in the active-site pocket of CDO. Putative bacterial CDOs predicted to have CDO activity were shown to have similar substrate specificity and kinetic parameters as eukaryotic CDOs. Information gleaned from crystal structures of mammalian CDO along with sequence information for homologs shown to have CDO activity facilitated the identification of a CDO family fingerprint motif. One key feature of the CDO fingerprint motif is that the canonical metal-binding glutamate residue in cupin motif 1 is replaced by a cysteine (in mammalian CDOs) or by a glycine (bacterial CDOs). The recent report that some putative bacterial CDO homologs are actually 3-mercaptopropionate dioxygenases suggests that the CDO family may include proteins with specificities for other thiol substrates. A paralog of CDO in mammals was also identified and shown to be the other mammalian thiol dioxygenase, cysteamine dioxygenase (ADO). A tentative fingerprint motif for ADOs, or DUF1637 family members, is proposed. In ADOs, the conserved glutamate residue in cupin motif 1 is replaced by either glycine or valine. Both ADOs and CDOs appear to represent unique clades within the cupin superfamily.     

Crystal structure of a member of a novel family of dioxygenases (PF10014) reveals a conserved cupin fold and active site

PF10014 is a novel family of 2‐oxyglutarate‐Fe²⁺‐dependent dioxygenases that are involved in biosynthesis of antibiotics and regulation of biofilm formation, likely by catalyzing hydroxylation of free amino acids or other related ligands. The crystal structure of a PF10014 member from Methylibium petroleiphilum at 1.9 Å resolution shows strong structural similarity to cupin dioxygenases in overall fold and active site, despite very remote homology. However, one of the β‐strands of the cupin catalytic core is replaced by a loop that displays conformational isomerism that likely regulates the active site. Proteins 2014; 82:164–170. © 2013 Wiley Periodicals, Inc.     

Solution and crystal structure of BA42, a protein from the Antarctic bacterium Bizionia argentinensis comprised of a stand‐alone TPM domain

The structure of the BA42 protein belonging to the Antarctic flavobacterium Bizionia argentinensis was determined by nuclear magnetic resonance and X‐ray crystallography. This is the first structure of a member of the PF04536 family comprised of a stand‐alone TPM domain. The structure reveals a new topological variant of the four β‐strands constituting the central β‐sheet of the αβα architecture and a double metal binding site stabilizing a pair of crossing loops, not observed in previous structures of proteins belonging to this family. BA42 shows differences in structure and dynamics in the presence or absence of bound metals. The affinity for divalent metal ions is close to that observed in proteins that modulate their activity as a function of metal concentration, anticipating a possible role for BA42. Proteins 2014; 82:3062–3078. © 2014 Wiley Periodicals, Inc.     

Evolution of functional diversity in the cupin superfamily

The cupin superfamily of proteins is among the most functionally diverse of any described to date. It was named on the basis of the conserved β-barrel fold (‘cupa’ is the Latin term for a small barrel), and comprises both enzymatic and non-enzymatic members, which have either one or two cupin domains. Within the conserved tertiary structure, the variety of biochemical function is provided by minor variation of the residues in the active site and the identity of the bound metal ion. This review discusses the advantages of this particular scaffold and provides an evolutionary analysis of 18 different subclasses within the cupin superfamily.     

Solution NMR structures of the C‐domain of Tetrahymena cytoskeletal protein Tcb2 reveal distinct calcium‐induced structural rearrangements

Tcb2 is a calcium‐binding protein that localizes to the membrane‐associated skeleton of the ciliated protozoan Tetrahymena thermophila with hypothesized roles in ciliary movement, cell cortex signaling, and pronuclear exchange. Tcb2 has also been implicated in a unique calcium‐triggered, ATP‐independent type of contractility exhibited by filamentous networks isolated from the Tetrahymena cytoskeleton. To gain insight into Tcb2's structure‐function relationship and contractile properties, we determined solution NMR structures of its C‐terminal domain in the calcium‐free and calcium‐bound states. The overall architecture is similar to other calcium‐binding proteins, with paired EF‐hand calcium‐binding motifs. Comparison of the two structures reveals that Tcb2‐C's calcium‐induced conformational transition differs from the prototypical calcium sensor calmodulin, suggesting that the two proteins play distinct functional roles in Tetrahymena and likely have different mechanisms of target recognition. Future studies of the full‐length protein and the identification of Tcb2 cellular targets will help establish the molecular basis of Tcb2 function and its unique contractile properties. Proteins 2016; 84:1748–1756. © 2016 Wiley Periodicals, Inc.     

Structure of oxalate decarboxylase from Bacillus subtilis at 1.75 A resolution

Oxalate decarboxylase is a manganese-dependent enzyme that catalyzes the conversion of oxalate to formate and carbon dioxide. We have determined the structure of oxalate decarboxylase from Bacillus subtilis at 1.75 A resolution in the presence of formate. The structure reveals a hexamer with 32-point symmetry in which each monomer belongs to the cupin family of proteins. Oxalate decarboxylase is further classified as a bicupin because it contains two cupin folds, possibly resulting from gene duplication. Each oxalate decarboxylase cupin domain contains one manganese binding site. Each of the oxalate decarboxylase domains is structurally similar to oxalate oxidase, which catalyzes the manganese-dependent oxidative decarboxylation of oxalate to carbon dioxide and hydrogen peroxide. Amino acid side chains in the two metal binding sites of oxalate decarboxylase and the metal binding site of oxalate oxidase are very similar. Four manganese binding residues (three histidines and one glutamate) are conserved as well as a number of hydrophobic residues. The most notable difference is the presence of Glu333 in the metal binding site of the second cupin domain of oxalate decarboxylase. We postulate that this domain is responsible for the decarboxylase activity and that Glu333 serves as a proton donor in the production of formate. Mutation of Glu333 to alanine reduces the catalytic activity by a factor of 25. The function of the other domain in oxalate decarboxylase is not yet known.     

Correlated mutation analyses on super‐family alignments reveal functionally important residues

Correlated mutation analyses (CMA) on multiple sequence alignments are widely used for the prediction of the function of amino acids. The accuracy of CMA‐based predictions is mainly determined by the number of sequences, by their evolutionary distances, and by the quality of the alignments. These criteria are best met in structure‐based sequence alignments of large super‐families. So far, CMA‐techniques have mainly been employed to study the receptor interactions. The present work shows how a novel CMA tool, called Comulator, can be used to determine networks of functionally related residues in enzymes. These analyses provide leads for protein engineering studies that are directed towards modification of enzyme specificity or activity. As proof of concept, Comulator has been applied to four enzyme super‐families: the isocitrate lyase/phoshoenol‐pyruvate mutase super‐family, the hexokinase super‐family, the RmlC‐like cupin super‐family, and the FAD‐linked oxidases super‐family. In each of those cases networks of functionally related residue positions were discovered that upon mutation influenced enzyme specificity and/or activity as predicted. We conclude that CMA is a powerful tool for redesigning enzyme activity and selectivity. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

A hybrid NMR/SAXS‐based approach for discriminating oligomeric protein interfaces using Rosetta

Oligomeric proteins are important targets for structure determination in solution. While in most cases the fold of individual subunits can be determined experimentally, or predicted by homology‐based methods, protein–protein interfaces are challenging to determine de novo using conventional NMR structure determination protocols. Here we focus on a member of the bet‐V1 superfamily, Aha1 from Colwellia psychrerythraea. This family displays a broad range of crystallographic interfaces none of which can be reconciled with the NMR and SAXS data collected for Aha1. Unlike conventional methods relying on a dense network of experimental restraints, the sparse data are used to limit conformational search during optimization of a physically realistic energy function. This work highlights a new approach for studying minor conformational changes due to structural plasticity within a single dimeric interface in solution. Proteins 2015; 83:309–317. © 2014 Wiley Periodicals, Inc.     

Crystal structure of the Legionella pneumophila lem10 effector reveals a new member of the HD protein superfamily

Legionella pneumophila, the intracellular pathogen that can cause severe pneumonia known as Legionnaire's disease, translocates close to 300 effectors inside the host cell using Dot/Icm type IVB secretion system. The structure and function for the majority of these effector proteins remains unknown. Here, we present the crystal structure of the L. pneumophila effector Lem10. The structure reveals a multidomain organization with the largest C‐terminal domain showing strong structural similarity to the HD protein superfamily representatives. However, Lem10 lacks the catalytic His‐Asp residue pair and does not show any in vitro phosphohydrolase enzymatic activity, typical for HD proteins. While the biological function of Lem10 remains elusive, our analysis shows that similar distinct features are shared by a significant number of HD domains found in Legionella proteins, including the SidE family of effectors known to play an important role during infection. Taken together our data point to the presence of a specific group of non‐catalytic Legionella HD domains, dubbed LHDs, which are involved in pathogenesis. Proteins 2015; 83:2319–2325. © 2015 Wiley Periodicals, Inc.     

Structure and Functional Analysis of YcfD,a Novel 2-Oxoglutarate/Fe-Dependent Oxygenase Involved in Translational Regulation in Escherichia coli

The 2-oxoglutarate (2OG)/Fe^2 +-dependent oxygenases (2OG oxygenases) are a large family of proteins that share a similar overall three-dimensional structure and catalyze a diverse array of oxidation reactions. The Jumonji C (JmjC)-domain-containing proteins represent an important subclass of the 2OG oxygenase family that typically catalyze protein hydroxylation; however, recently, other reactions have been identified, such as tRNA modification. The Escherichia coli gene, ycfD, was predicted to be a JmjC-domain-containing protein of unknown function based on primary sequence. Recently, YcfD was determined to act as a ribosomal oxygenase, hydroxylating an arginine residue on the 50S ribosomal protein L-16 (RL-16). We have determined the crystal structure of YcfD at 2.7 Å resolution, revealing that YcfD is structurally similar to known JmjC proteins and possesses the characteristic double-stranded β-helix fold or cupin domain. Separate from the cupin domain, an additional globular module termed α-helical arm mediates dimerization of YcfD. We further have shown that 2OG binds to YcfD using isothermal titration calorimetry and identified key binding residues using mutagenesis that, together with the iron location and structural similarity with other cupin family members, allowed identification of the active site. Structural homology to ribosomal assembly proteins combined with GST (glutathione S-transferase)-YcfD pull-down of a ribosomal protein and docking of RL-16 to the YcfD active site support the role of YcfD in regulation of bacterial ribosome assembly. Furthermore, overexpression of YcfD is shown to inhibit cell growth signifying a toxic effect on ribosome assembly.     

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.     

Genome‐wide computational determination of the human metalloproteome

Accurate prediction of protein function in humans is important for understanding biological processes at the molecular level in biomedicine and drug design. Over a third of proteins are commonly held to bind metal, and ～10% of human proteins, to bind zinc. Therefore, an initial step in protein function prediction frequently involves predicting metal ion binding. In recent years, methods have been developed to predict a set of residues in 3D space forming the metal‐ion binding site, often with a high degree of accuracy. Here, using extensions of these methods, we provide an extensive list of human proteins and their putative metal ion binding site residues, using translated gene sequences derived from the complete, resolved human genome. Under conditions of ～90% selectivity, over 900 new human putative metal ion binding proteins are identified. A statistical analysis of resolved metal ion binding sites in the human metalloproteome is furnished and the importance of remote homology analysis is demonstrated. As an example, a novel metal‐ion binding site involving a complex of a botulinum substrate with its inhibitor is presented. On the basis of the location of the predicted site and the interactions of the contacting residues at the complex interface, we postulate that metal ion binding in this region could influence complex formation and, consequently, the functioning of the protein. Thus, this work provides testable hypotheses about novel functions of known proteins. Proteins 2015; 83:931–939. © 2015 Wiley Periodicals, Inc.     

Crystal structure of a [NiFe] hydrogenase maturation protease HybD from Thermococcus kodakarensis KOD1

A [NiFe] hydrogenase maturation protease HybD from Thermococcus kodakarensis KOD1 (TkHybD) is involved in the cleavage of the C‐terminal residues of [NiFe] hydrogenase large subunits by Ni recognition. Here, we report the crystal structure of TkHybD at 1.82 Å resolution to better understand this process. TkHybD exhibits an α/β/α sandwich fold with conserved residues responsible for the Ni recognition. Comparisons of TkHybD with homologous proteins also reveal that they share a common overall architecture, suggesting that they have similar catalytic functions. Our results including metal binding site prediction provide insight into the substrate recognition and catalysis mechanism of TkHybD. Proteins 2016; 84:1321–1327. © 2016 Wiley Periodicals, Inc.     

Discovery of binding proteins for a protein target using protein–protein docking‐based virtual screening

Target structure‐based virtual screening, which employs protein‐small molecule docking to identify potential ligands, has been widely used in small‐molecule drug discovery. In the present study, we used a protein–protein docking program to identify proteins that bind to a specific target protein. In the testing phase, an all‐to‐all protein–protein docking run on a large dataset was performed. The three‐dimensional rigid docking program SDOCK was used to examine protein–protein docking on all protein pairs in the dataset. Both the binding affinity and features of the binding energy landscape were considered in the scoring function in order to distinguish positive binding pairs from negative binding pairs. Thus, the lowest docking score, the average Z‐score, and convergency of the low‐score solutions were incorporated in the analysis. The hybrid scoring function was optimized in the all‐to‐all docking test. The docking method and the hybrid scoring function were then used to screen for proteins that bind to tumor necrosis factor‐α (TNFα), which is a well‐known therapeutic target for rheumatoid arthritis and other autoimmune diseases. A protein library containing 677 proteins was used for the screen. Proteins with scores among the top 20% were further examined. Sixteen proteins from the top‐ranking 67 proteins were selected for experimental study. Two of these proteins showed significant binding to TNFα in an in vitro binding study. The results of the present study demonstrate the power and potential application of protein–protein docking for the discovery of novel binding proteins for specific protein targets. Proteins 2014; 82:2472–2482. © 2014 Wiley Periodicals, Inc.     

Function‐based assessment of structural similarity measurements using metal co‐factor orientation

Structure comparison is widely used to quantify protein relationships. Although there are several approaches to calculate structural similarity, specifying significance thresholds for similarity metrics is difficult due to the inherent likeness of common secondary structure elements. In this study, metal co‐factor location is used to assess the biological relevance of structural alignments. The distance between the centroids of bound co‐factors adds a chemical and function‐relevant constraint to the structural superimposition of two proteins. This additional dimension can be used to define cut‐off values for discriminating valid and spurious alignments in large alignment sets. The hypothesis underlying our approach is that metal coordination sites constrain structural evolution, thus revealing functional relationships between distantly related proteins. A comparison of three related nitrogenases shows the sequence and fold constraints imposed on the protein structures up to 18 Å away from the centers of their bound metal clusters. Proteins 2014; 82:648–656. © 2013 Wiley Periodicals, Inc.     

Analyzing the effect of homogeneous frustration in protein folding

The energy landscape theory has been an invaluable theoretical framework in the understanding of biological processes such as protein folding, oligomerization, and functional transitions. According to the theory, the energy landscape of protein folding is funneled toward the native state, a conformational state that is consistent with the principle of minimal frustration. It has been accepted that real proteins are selected through natural evolution, satisfying the minimum frustration criterion. However, there is evidence that a low degree of frustration accelerates folding. We examined the interplay between topological and energetic protein frustration. We employed a C_α structure‐based model for simulations with a controlled nonspecific energetic frustration added to the potential energy function. Thermodynamics and kinetics of a group of 19 proteins are completely characterized as a function of increasing level of energetic frustration. We observed two well‐separated groups of proteins: one group where a little frustration enhances folding rates to an optimal value and another where any energetic frustration slows down folding. Protein energetic frustration regimes and their mechanisms are explained by the role of non‐native contact interactions in different folding scenarios. These findings strongly correlate with the protein free‐energy folding barrier and the absolute contact order parameters. These computational results are corroborated by principal component analysis and partial least square techniques. One simple theoretical model is proposed as a useful tool for experimentalists to predict the limits of improvements in real proteins.Proteins 2013; 81:1727–1737. © 2013 Wiley Periodicals, Inc.     

Is the bovine lysosomal phospholipase B‐like protein an amidase?

The main function of lysosomal proteins is to degrade cellular macromolecules. We purified a novel lysosomal protein to homogeneity from bovine kidneys. By gene annotation, this protein is defined as a bovine phospholipase B‐like protein 1 (bPLBD1) and, to better understand its biological function, we solved its structure at 1.9 Å resolution. We showed that bPLBD1 has uniform noncomplex‐type N‐glycosylation and that it localized to the lysosome. The first step in lysosomal protein transport, the initiation of mannose‐6‐phosphorylation by a N‐acetylglucosamine‐1‐phosphotransferase, requires recognition of at least two distinct lysines on the protein surface. We identified candidate lysines by analyzing the structural and sequentially conserved N‐glycosylation sites and lysines in bPLBD1 and in the homologous mouse PLBD2. Our model suggests that N408 is the primarily phosphorylated glycan, and K358 a key residue for N‐acetylglucosamine‐1‐phosphotransferase recognition. Two other lysines, K334 and K342, provide the required second site for N‐acetylglucosamine‐1‐phosphotransferase recognition. bPLBD1 is an N‐terminal nucleophile (Ntn) hydrolase. By comparison with other Ntn‐hydrolases, we conclude that the acyl moiety of PLBD1 substrate must be small to fit the putative binding pocket, whereas the space for the rest of the substrate is a large open cleft. Finally, as all the known substrates of Ntn‐hydrolases have amide bonds, we suggest that bPLBD1 may be an amidase or peptidase instead of lipase, explaining the difficulty in finding a good substrate for any members of the PLBD family. Proteins 2014; 82:300–311. © 2013 Wiley Periodicals, Inc.     

Proteins

Proteins continue to surprise and amaze us in the myriad of ways in which they achieve biological function. The Proteins section in this issue of Current Opinion in Structural Biology highlights several proteins in which large conformational changes and evolutionary divergence in structure and function, play essential roles in their adaptation to a variety of biological functions. In addition, fundamental advances have been made in research, spurred on by industrial interest in the use of proteins as drug targets or as catalysts. All of the reviews in this section document the fact that multiple crystal structures of a protein in different functional states, and of different members of protein families, are necessary for the composition of a complete structural picture.     

Structural features that predict real-value fluctuations of globular proteins

It is crucial to consider dynamics for understanding the biological function of proteins. We used a large number of molecular dynamics (MD) trajectories of nonhomologous proteins as references and examined static structural features of proteins that are most relevant to fluctuations. We examined correlation of individual structural features with fluctuations and further investigated effective combinations of features for predicting the real value of residue fluctuations using the support vector regression (SVR). It was found that some structural features have higher correlation than crystallographic B‐factors with fluctuations observed in MD trajectories. Moreover, SVR that uses combinations of static structural features showed accurate prediction of fluctuations with an average Pearson's correlation coefficient of 0.669 and a root mean square error of 1.04 Å. This correlation coefficient is higher than the one observed in predictions by the Gaussian network model (GNM). An advantage of the developed method over the GNMs is that the former predicts the real value of fluctuation. The results help improve our understanding of relationships between protein structure and fluctuation. Furthermore, the developed method provides a convienient practial way to predict fluctuations of proteins using easily computed static structural features of proteins. Proteins 2012; © 2012 Wiley Periodicals, Inc.