The evolving world of pseudoenzymes: proteins, prejudice and zombies

Evaluation of the Structure-Function Relationship of SGNH Lipase from Streptomyces rimosus by Site-Directed Mutagenesis and Computational Approach

Streptomyces rimosus extracellular lipase (SrL) is a multifunctional hydrolase belonging to the SGNH family. Here site-directed mutagenesis (SDM) was used for the first time to investigate the functional significance of the conserved amino acid residues Ser10, Gly54, Asn82, Asn213, and His216 in the active site of SrL. The hydrolytic activity of SrL variants was determined using para-nitrophenyl (pNP) esters with C4, C8, and C16 fatty acid chains. Mutation of Ser10, Asn82, or His216, but not Gly54, to Ala abolished lipase activity for all substrates. In contrast, the Asn213Ala variant showed increased enzymatic activity for C8 and C16 pNP esters. Molecular dynamics (MD) simulations showed that the interactions between the long alkyl chain substrate (C16) and Ser10 and Asn82 were strongest in Asn213Ala SrL. In addition to Asn82, Gly54, and Ser10, several new constituents of the substrate binding site were recognized (Lys28, Ser53, Thr89, and Glu212), as well as strong electrostatic interactions between Lys28 and Glu212. In addition to the H bonds Ser10-His216 and His216-Ser214, Tyr11 interacted strongly with Ser10 and His216 in all complexes with an active enzyme form. A previously unknown strong H bond between the catalytically important Asn82 and Gly54 was uncovered, which stabilizes the substrate in an orientation suitable for the enzyme reaction.

The pseudokinase TRIB3 controls adipocyte lipid homeostasis and proliferation in vitro and in vivo

OBJECTIVE

In vivo studies in humans and mice have implicated the pseudokinase Tribbles 3 (TRIB3) in various aspects of energy metabolism. Whilst cell-based studies indicate a role for TRIB3 in adipocyte differentiation and function, it is unclear if and how these cellular functions may contribute to overall metabolic health.

METHODS

We investigated the metabolic phenotype of whole-body Trib3 knockout (Trib3KO) mice, focusing on adipocyte and adipose tissue functions. In addition, we combined lipidomics, transcriptomics, interactomics and phosphoproteomics analyses to elucidate cell-intrinsic functions of TRIB3 in pre- and mature adipocytes.

RESULTS

Trib3KO mice display increased adiposity, but their insulin sensitivity remains unaltered. Trib3KO adipocytes are smaller and display higher Proliferating Cell Nuclear Antigen (PCNA) levels, indicating potential alterations in either i) proliferation-differentiation balance, ii) impaired expansion after cell division, or iii) an altered balance between lipid storage and release, or a combination thereof. Lipidome analyses suggest TRIB3 involvement in the latter two processes, as triglyceride storage is reduced and membrane composition, which can restrain cellular expansion, is altered. Integrated interactome, phosphoproteome and transcriptome analyses support a role for TRIB3 in all three cellular processes through multiple cellular pathways, including Mitogen Activated Protein Kinase- (MAPK/ERK), Protein Kinase A (PKA)-mediated signaling and Transcription Factor 7 like 2 (TCF7L2) and Beta Catenin-mediated gene expression.

CONCLUSIONS

Our findings support TRIB3 playing multiple distinct regulatory roles in the cytoplasm, nucleus and mitochondria, ultimately controlling adipose tissue homeostasis, rather than affecting a single cellular pathway.

Enzyme function and evolution through the lens of bioinformatics

Enzymes have been shaped by evolution over billions of years to catalyse the chemical reactions that support life on earth. Dispersed in the literature, or organised in online databases, knowledge about enzymes can be structured in distinct dimensions, either related to their quality as biological macromolecules, such as their sequence and structure, or related to their chemical functions, such as the catalytic site, kinetics, mechanism, and overall reaction. The evolution of enzymes can only be understood when each of these dimensions is considered. In addition, many of the properties of enzymes only make sense in the light of evolution. We start this review by outlining the main paradigms of enzyme evolution, including gene duplication and divergence, convergent evolution, and evolution by recombination of domains. In the second part, we overview the current collective knowledge about enzymes, as organised by different types of data and collected in several databases. We also highlight some increasingly powerful computational tools that can be used to close gaps in understanding, in particular for types of data that require laborious experimental protocols. We believe that recent advances in protein structure prediction will be a powerful catalyst for the prediction of binding, mechanism, and ultimately, chemical reactions. A comprehensive mapping of enzyme function and evolution may be attainable in the near future.

Current successes and remaining challenges in protein function prediction

In recent years, improvements in protein function prediction methods have led to increased success in annotating protein sequences. However, the functions of over 30% of protein-coding genes remain unknown for many sequenced genomes. Protein functions vary widely, from catalyzing chemical reactions to binding DNA or RNA or forming structures in the cell, and some types of functions are challenging to predict due to the physical features associated with those functions. Other complications in understanding protein functions arise due to the fact that many proteins have more than one function or very small differences in sequence or structure that correspond to different functions. We will discuss some of the recent developments in predicting protein functions and some of the remaining challenges.

Natural Parameter Conditions for Singular Perturbations of Chemical and Biochemical Reaction Networks

We consider reaction networks that admit a singular perturbation reduction in a certain parameter range. The focus of this paper is on deriving "small parameters" (briefly for small perturbation parameters), to gauge the accuracy of the reduction, in a manner that is consistent, amenable to computation and permits an interpretation in chemical or biochemical terms. Our work is based on local timescale estimates via ratios of the real parts of eigenvalues of the Jacobian near critical manifolds. This approach modifies the one introduced by Segel and Slemrod and is familiar from computational singular perturbation theory. While parameters derived by this method cannot provide universal quantitative estimates for the accuracy of a reduction, they represent a critical first step toward this end. Working directly with eigenvalues is generally unfeasible, and at best cumbersome. Therefore we focus on the coefficients of the characteristic polynomial to derive parameters, and relate them to timescales. Thus, we obtain distinguished parameters for systems of arbitrary dimension, with particular emphasis on reduction to dimension one. As a first application, we discuss the Michaelis-Menten reaction mechanism system in various settings, with new and perhaps surprising results. We proceed to investigate more complex enzyme catalyzed reaction mechanisms (uncompetitive, competitive inhibition and cooperativity) of dimension three, with reductions to dimension one and two. The distinguished parameters we derive for these three-dimensional systems are new. In fact, no rigorous derivation of small parameters seems to exist in the literature so far. Numerical simulations are included to illustrate the efficacy of the parameters obtained, but also to show that certain limitations must be observed.

Dissecting the structural heterogeneity of proteins by native mass spectrometry

A single gene yields many forms of proteins via combinations of post-transcriptional/post-translational modifications. Proteins also fold into higher-order structures and interact with other molecules. The combined molecular diversity leads to the heterogeneity of proteins that manifests as distinct phenotypes. Structural biology has generated vast amounts of data, effectively enabling accurate structural prediction by computational methods. However, structures are often obtained heterologously under homogeneous states in vitro. The lack of native heterogeneity under cellular context creates challenges in precisely connecting the structural data to phenotypes. Mass spectrometry (MS) based proteomics methods can profile proteome composition of complex biological samples. Most MS methods follow the "bottom-up" approach, which denatures and digests proteins into short peptide fragments for ease of detection. Coupled with chemical biology approaches, higher-order structures can be probed via incorporation of covalent labels on native proteins that are maintained at the peptide level. Alternatively, native MS follows the "top-down" approach and directly analyzes intact proteins under nondenaturing conditions. Various tandem MS activation methods can dissect the intact proteins for in-depth structural elucidation. Herein, we review recent native MS applications for characterizing heterogeneous samples, including proteins binding to mixtures of ligands, homo/hetero-complexes with varying stoichiometry, intrinsically disordered proteins with dynamic conformations, glycoprotein complexes with mixed modification states, and active membrane protein complexes in near-native membrane environments. We summarize the benefits, challenges, and ongoing developments in native MS, with the hope to demonstrate an emerging technology that complements other tools by filling the knowledge gaps in understanding molecular heterogeneity of proteins. This article is protected by copyright. All rights reserved.

Vicinal glutamates are better phosphomimetics: Phosphorylation is required for allosteric activation of guanylyl cyclase-A

Multisite phosphorylation of guanylyl cyclase (GC)-A, also known as NPR-A or NPR1, is required for receptor activation by natriuretic peptides (NPs) because alanine substitutions for the first four GC-A phosphorylation sites produce an enzyme that cannot be stimulated by NPs. In contrast, single Glu substitutions for the first six chemically identified GC-A phosphorylation sites to mimic the negative charge of phosphate produced an enzyme that is activated by NPs but had an elevated Michaelis constant (Km), resulting in low activity. Here, we show that vicinal (double adjacent) Glu substitutions for the same sites to mimic the two negative charges of phosphate produced a near wild type (WT) enzyme with a low Km. Unlike the enzyme with single glutamate substitutions, the vicinally substituted enzyme did not require the functionally identified Ser-473-Glu substitution to achieve WT-like activity. Importantly, the negative charge associated with either phosphorylation or glutamate substitutions was required for allosteric activation of GC-A by ATP. We conclude that vicinal Glu substitutions are better phosphomimetics than single Glu substitutions and that phosphorylation is required for allosteric activation of GC-A in the absence and presence of NP. Finally, we suggest that the putative functionally identified phosphorylation sites, Ser-473 in GC-A and Ser-489 in GC-B, are not phosphorylation sites at all.

Insights into acylation mechanisms: co-expression of serine carboxypeptidase-like acyltransferases and their non-catalytic companion paralogs

Serine carboxypeptidase-like acyltransferases (SCPL-ATs) play a vital role in the diversification of plant metabolites. Galloylated flavan-3-ols highly accumulate in tea (Camellia sinensis), grape (Vitis vinifera), and persimmon (Diospyros kaki). To date, the biosynthetic mechanism of these compounds remains unknown. Herein, we report that two SCPL-AT paralogs are involved in galloylation of flavan-3-ols: CsSCPL4, which contains the conserved catalytic triad S-D-H, and CsSCPL5, which has the alternative triad T-D-Y. Integrated data from transgenic plants, recombinant enzymes, and gene mutations showed that CsSCPL4 is a catalytic acyltransferase, while CsSCPL5 is a non-catalytic companion paralog (NCCP). Co-expression of CsSCPL4 and CsSCPL5 is likely responsible for the galloylation. Furthermore, pull-down and co-immunoprecipitation assays showed that CsSCPL4 and CsSCPL5 interact, increasing protein stability and promoting post-translational processing. Moreover, phylogenetic analyses revealed that their homologs co-exist in galloylated flavan-3-ol- or hydrolyzable tannin-rich plant species. Enzymatic assays further revealed the necessity of co-expression of those homologs for acyltransferase activity. Evolution analysis revealed that the mutations of the CsSCPL5 catalytic residues may have taken place about 10 million years ago. These findings show that the co-expression of SCPL-ATs and their NCCPs contributes to the acylation of flavan-3-ols in the plant kingdom.

Protein engineering strategies to stimulate the functions of bacterial pseudokinases

Enzymes orchestrate an array of concerted functions that often culminate in the chemical conversion of substrates into products. In the bacterial kingdom, histidine kinases autophosphorylate, then transfer that phosphate to a second protein called a response regulator. Bacterial genomes can encode large numbers of histidine kinases that provide surveillance of environmental and cytosolic stresses through signal stimulation of histidine kinase activity. Pseudokinases lack these hallmark catalytic functions but often retain binding interactions and allostery. Characterization of bacterial pseudokinases then takes a fundamentally different approach than their enzymatic counterparts. Here we discuss models for how bacterial pseudokinases can utilize protein-protein interactions and allostery to serve as crucial signaling pathway regulators. Then we describe a protein engineering strategy to interrogate these models, emphasizing how signals flow within bacterial pseudokinases. This description includes design considerations, cloning strategies, and the purification of leucine zippers fused to pseudokinases. We then describe two assays to interrogate this approach. First is a C. crescentus swarm plate assay to track motility phenotypes related to a bacterial pseudokinase. Second is an in vitro coupled-enzyme assay that can be applied to test if and how a pseudokinase regulates an active kinase. Together these approaches provide a blueprint for dissecting the mechanisms of cryptic bacterial pseudokinases.

Co-expression of recombinant RIPK3:MLKL complexes using the baculovirus-insect cell system

Pseudokinase domains are found throughout the kingdoms of life and serve myriad roles in cell signaling. These domains, which resemble protein kinases but are catalytically-deficient, have been described principally as protein interaction domains. Broadly, pseudokinases have been reported to function as: allosteric regulators of conventional enzymes; scaffolds to nucleate assembly and/or localization of signaling complexes; molecular switches; or competitors of signaling complex assembly. From detailed structural and biochemical studies of individual pseudokinases, a picture of how they mediate protein interactions is beginning to emerge. Many such studies have relied on recombinant protein production in insect cells, where endogenous chaperones and modifying enzymes favor bona fide folding of pseudokinases. Here, we describe methods for co-expression of pseudokinases and their interactors in insect cells, as exemplified by the MLKL pseudokinase, which is the terminal effector in the necroptosis cell death pathway, and its upstream regulator kinase RIPK3. These methods are broadly applicable to co-expression of other pseudokinases with their interaction partners from bacmids using the baculovirus-insect cell expression system.

CRISPR deletions in cell lines for reconstitution studies of pseudokinase function

The non-catalytic cousins of protein kinases, the pseudokinases, have grown to prominence as indispensable signaling entities over the past decade, despite their lack of catalytic activity. Because their importance has only been fully embraced recently, many of the 10% of the human kinome categorized as pseudokinases are yet to be attributed biological functions. The advent of CRISPR-Cas9 editing to genetically delete pseudokinases in a cell line of interest has proven invaluable to dissecting many functions and remains the method of choice for gene knockout. Here, using the terminal effector pseudokinase in the necroptosis cell death pathway, MLKL, as an exemplar, we describe a method for genetic knockout of pseudokinases in cultured cells. This method does not retain the CRISPR guide sequence in the edited cells, which eliminates possible interference in subsequent reconstitution studies where mutant forms of the pseudokinase can be reintroduced into cells exogenously for detailed mechanistic characterization.

The Ancient Origins of Death Domains Support the 'Original Sin' Hypothesis for the Evolution of Programmed Cell Death

The discovery of caspase homologs in bacteria highlighted the relationship between programmed cell death (PCD) evolution and eukaryogenesis. However, the origin of PCD genes in prokaryotes themselves (bacteria and archaea) is poorly understood and a source of controversy. Whether archaea also contain C14 peptidase enzymes and other death domains is largely unknown because of a historical dearth of genomic data. Archaeal genomic databases have grown significantly in the last decade, which allowed us to perform a detailed comparative study of the evolutionary histories of PCD-related death domains in major archaeal phyla, including the deepest branching phyla of Candidatus Aenigmarchaeota, Candidatus Woesearchaeota, and Euryarchaeota. We identified death domains associated with executioners of PCD, like the caspase homologs of the C14 peptidase family, in 321 archaea sequences. Of these, 15.58% were metacaspase type I orthologues and 84.42% were orthocaspases. Maximum likelihood phylogenetic analyses revealed a scattered distribution of orthocaspases and metacaspases in deep-branching bacteria and archaea. The tree topology was incongruent with the prokaryote 16S phylogeny suggesting a common ancestry of PCD genes in prokaryotes and subsequent massive horizontal gene transfer coinciding with the divergence of archaea and bacteria. Previous arguments for the origin of PCD were philosophical in nature with two popular propositions being the "addiction" and 'original sin' hypotheses. Our data support the 'original sin' hypothesis, which argues for a pleiotropic origin of the PCD toolkit with pro-life and pro-death functions tracing back to the emergence of cellular life-the Last Universal Common Ancestor State.

Pseudophosphatases as Regulators of MAPK Signaling

Mitogen-activated protein kinase (MAPK) signaling pathways are highly conserved regulators of eukaryotic cell function. These enzymes regulate many biological processes, including the cell cycle, apoptosis, differentiation, protein biosynthesis, and oncogenesis; therefore, tight control of the activity of MAPK is critical. Kinases and phosphatases are well established as MAPK activators and inhibitors, respectively. Kinases phosphorylate MAPKs, initiating and controlling the amplitude of the activation. In contrast, MAPK phosphatases (MKPs) dephosphorylate MAPKs, downregulating and controlling the duration of the signal. In addition, within the past decade, pseudoenzymes of these two families, pseudokinases and pseudophosphatases, have emerged as bona fide signaling regulators. This review discusses the role of pseudophosphatases in MAPK signaling, highlighting the function of phosphoserine/threonine/tyrosine-interacting protein (STYX) and TAK1-binding protein (TAB 1) in regulating MAPKs. Finally, a new paradigm is considered for this well-studied cellular pathway, and signal transduction pathways in general.

Identification and Functional Analysis of a Pseudo-Cysteine Protease from the Midgut Transcriptome of Sphenophorus levis

The Sphenophorus levis (Coleoptera, Curculionidae) is one of the main pests of sugarcane in Brazil. Although its major digestive proteases are known, its complex digestive process still needs to be further understood. We constructed a transcriptome from the midgut of 30-day-old larvae and identified sequences similar to its major digestive protease (cysteine cathepsin Sl-CathL), however, they presented a different amino acid than cysteine in the active cleft. We identified, recombinantly produced, and characterized Sl-CathL-CS, a pseudo cysteine protease, and verified that higher gene expression levels of Sl-CathL-CS occur in the midgut of 30-day old larvae. We reverted the serine residue to cysteine and compared the activity of the mutant (Sl-CathL-mutSC) with Sl-CathL-CS. Sl-CathL-CS presented no protease activity, but Sl-CathL-mutSC hydrolyzed Z-Phe-Arg-AMC (V_max = 1017.60 ± 135.55, K_m = 10.77 mM) and was inhibited by a cysteine protease inhibitor E-64 (K_i = 38.52 ± 1.20 μM), but not by the serine protease inhibitor PMSF. Additionally, Sl-CathL-CS interacted with a sugarcane cystatin, while Sl-CathL-mutSC presented weaker interaction. Finally, protein ligand docking reinforced the differences in the catalytic sites of native and mutant proteins. These results indicate that Sl-CathL-CS is a pseudo-cysteine protease that assists protein digestion possibly by interacting with canecystatins, allowing the true proteases to work.

Inactive pseudoenzyme subunits in heterotetrameric BbsCD, a novel short-chain alcohol dehydrogenase involved in anaerobic toluene degradation

Anaerobic toluene degradation proceeds by fumarate addition to produce (R)-benzylsuccinate as first intermediate, which is further degraded via β-oxidation by five enzymes encoded in the conserved bbs operon. This study characterizes two enzymes of this pathway, (E)-benzylidenesuccinyl-CoA hydratase (BbsH), and (S,R)-2-(α-hydroxybenzyl)succinyl-CoA dehydrogenase (BbsCD) from Thauera aromatica. BbsH, a member of the enoyl-CoA hydratase family, converts (E)-benzylidenesuccinyl-CoA to 2-(α-hydroxybenzyl)succinyl-CoA and was subsequently used in a coupled enzyme assay with BbsCD, which belongs to the short-chain dehydrogenases/reductase (SDR) family. The BbsCD crystal structure shows a C2-symmetric heterotetramer consisting of BbsC₂ and BbsD₂ dimers. BbsD subunits are catalytically active and capable of binding NAD⁺ and substrate, whereas BbsC subunits represent built-in pseudoenzyme moieties lacking all motifs of the SDR family required for substrate binding or catalysis. Molecular modeling studies predict that the active site of BbsD is specific for conversion of the (S,R)-diastereomer of 2-(α-hydroxybenzyl)succinyl-CoA to (S)-2-benzoylsuccinyl-CoA by hydride transfer to the re-face of nicotinamide adenine dinucleotide (NAD)⁺ . Furthermore, BbsC subunits are not engaged in substrate binding and merely serve as scaffold for the BbsD dimer. BbsCD represents a novel clade of related enzymes within the SDR family, which adopt a heterotetrameric architecture and catalyze the β-oxidation of aromatic succinate adducts.

Promoting a more integrated approach to structure and function

The connection between structure and function is one of the fundamental tenets of biology: a biological unit's structure determines its function, and, conversely, its function depends upon its structure. Historically, important advances have been made either when understanding of structure leads to questions about function or when understanding of function raises questions about the structures involved. Consequently, considering the connections between structure and function from a broader perspective might lead to the development of novel hypotheses that move our understanding of the fundamental connections between structure and function forward. Better integration of structure and function is a key component in the broader goal of reintegrating biology within and across scales. Here, we provide examples of how integrating studies of structure and function as well as comparing structure-function relationships across biological scales can lead to scientific advances. We also emphasize the potential of integrating studies of structure and function across scales for bio-inspired design and for improving biology education.

The Roles of Pseudophosphatases in Disease

The pseudophosphatases, atypical members of the protein tyrosine phosphatase family, have emerged as bona fide signaling regulators within the past two decades. Their roles as regulators have led to a renaissance of the pseudophosphatase and pseudoenyme fields, catapulting interest from a mere curiosity to intriguing and relevant proteins to investigate. Pseudophosphatases make up approximately fourteen percent of the phosphatase family, and are conserved throughout evolution. Pseudophosphatases, along with pseudokinases, are important players in physiology and pathophysiology. These atypical members of the protein tyrosine phosphatase and protein tyrosine kinase superfamily, respectively, are rendered catalytically inactive through mutations within their catalytic active signature motif and/or other important domains required for catalysis. This new interest in the pursuit of the relevant functions of these proteins has resulted in an elucidation of their roles in signaling cascades and diseases. There is a rapid accumulation of knowledge of diseases linked to their dysregulation, such as neuropathies and various cancers. This review analyzes the involvement of pseudophosphatases in diseases, highlighting the function of various role(s) of pseudophosphatases involvement in pathologies, and thus providing a platform to strongly consider them as key therapeutic drug targets.

The expanding world of protein kinase-like families in bacteria: forty families and counting

The protein kinase-like clan/superfamily is a large group of regulatory, signaling and biosynthetic enzymes that were historically regarded as typically eukaryotic proteins, although bacterial members have also been known for a long time. In this review, we explore the diversity of bacterial protein kinase like families, and discuss functional versatility of these enzymes, both the ones acting within the bacterial cell, and those acting within eukaryotic cells as effectors during infection. We focus on novel bacterial kinase-like families discovered in the last five years. A bioinformatics perspective is held here, hence sequence and structure comparison overview is presented, and also a comparison of genomic neighbourhoods of the families. We perform a phylum-level census of the families. Also, we discuss apparent pseudokinases that turned out to perform alternative catalytic functions by repurposing their atypical kinase-like active sites. We also highlight some 'unpopular' kinase-like families that await characterisation.

The upcycled roles of pseudoenzymes in two-component signal transduction

Upon first glance at a bacterial genome, pseudoenzymes appear unremarkable due to their lack of critical motifs that facilitate catalysis. These pseudoenzymes exist within signal transduction enzymes including histidine kinases, response regulators, diguanylate cyclases, and phosphodiesterases. Here, we summarize recent studies of bacterial pseudo-histidine kinases and pseudo-response regulators that regulate cell division, capsule formation, and the circadian rhythm. These examples illuminate the mechanistic potential of catalytically dead signaling enzymes and their impact upon bacterial signal transduction. Moreover, proteins lacking characteristic catalytic features of two-component systems reveal the sophisticated underlying potential of canonical two-component systems.

Alternative splicing creates a pseudo-strictosidine β-d-glucosidase modulating alkaloid synthesis in Catharanthus roseus

Deglycosylation is a key step in the activation of specialized metabolites involved in plant defense mechanisms. This reaction is notably catalyzed by β-glucosidases of the glycosyl hydrolase 1 (GH1) family such as strictosidine β-d-glucosidase (SGD) from Catharanthus roseus. SGD catalyzes the deglycosylation of strictosidine, forming a highly reactive aglycone involved in the synthesis of cytotoxic monoterpene indole alkaloids (MIAs) and in the crosslinking of aggressor proteins. By exploring C. roseus transcriptomic resources, we identified an alternative splicing event of the SGD gene leading to the formation of a shorter isoform of this enzyme (shSGD) that lacks the last 71-residues and whose transcript ratio with SGD ranges from 1.7% up to 42.8%, depending on organs and conditions. Whereas it completely lacks β-glucosidase activity, shSGD interacts with SGD and causes the disruption of SGD multimers. Such disorganization drastically inhibits SGD activity and impacts downstream MIA synthesis. In addition, shSGD disrupts the metabolic channeling of downstream biosynthetic steps by hampering the recruitment of tetrahydroalstonine synthase in cell nuclei. shSGD thus corresponds to a pseudo-enzyme acting as a regulator of MIA biosynthesis. These data shed light on a peculiar control mechanism of β-glucosidase multimerization, an organization common to many defensive GH1 members.

Resources and Methods for Engineering "Designer" Glycan-Binding Proteins

This review provides information on available methods for engineering glycan-binding proteins (GBP). Glycans are involved in a variety of physiological functions and are found in all domains of life and viruses. Due to their wide range of functions, GBPs have been developed with diagnostic, therapeutic, and biotechnological applications. The development of GBPs has traditionally been hindered by a lack of available glycan targets and sensitive and selective protein scaffolds; however, recent advances in glycobiology have largely overcome these challenges. Here we provide information on how to approach the design of novel "designer" GBPs, starting from the protein scaffold to the mutagenesis methods, selection, and characterization of the GBPs.

Biochemical and genetic analysis of Ecm14, a conserved fungal pseudopeptidase

Background

Like most major enzyme families, the M14 family of metallocarboxypeptidases (MCPs) contains a number of pseudoenzymes predicted to lack enzyme activity and with poorly characterized molecular function. The genome of the yeast Saccharomyces cerevisiae encodes one member of the M14 MCP family, a pseudoenzyme named Ecm14 proposed to function in the extracellular matrix. In order to better understand the function of such pseudoenzymes, we studied the structure and function of Ecm14 in S. cerevisiae.

Results

A phylogenetic analysis of Ecm14 in fungi found it to be conserved throughout the ascomycete phylum, with a group of related pseudoenzymes found in basidiomycetes. To investigate the structure and function of this conserved protein, His6-tagged Ecm14 was overexpressed in Sf9 cells and purified. The prodomain of Ecm14 was cleaved in vivo and in vitro by endopeptidases, suggesting an activation mechanism; however, no activity was detectable using standard carboxypeptidase substrates. In order to determine the function of Ecm14 using an unbiased screen, we undertook a synthetic lethal assay. Upon screening approximately 27,000 yeast colonies, twenty-two putative synthetic lethal clones were identified. Further analysis showed many to be synthetic lethal with auxotrophic marker genes and requiring multiple mutations, suggesting that there are few, if any, single S. cerevisiae genes that present synthetic lethal interactions with ecm14Δ.

Conclusions

We show in this study that Ecm14, although lacking detectable enzyme activity, is a conserved carboxypeptidase-like protein that is secreted from cells and is processed to a mature form by the action of an endopeptidase. Our study and datasets from other recent large-scale screens suggest a role for Ecm14 in processes such as vesicle-mediated transport and aggregate invasion, a fungal process that has been selected against in modern laboratory strains of S. cerevisiae.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12860-020-00330-w.

Phosphoglycerate kinase: structural aspects and functions, with special emphasis on the enzyme from Kinetoplastea

Phosphoglycerate kinase (PGK) is a glycolytic enzyme that is well conserved among the three domains of life. PGK is usually a monomeric enzyme of about 45 kDa that catalyses one of the two ATP-producing reactions in the glycolytic pathway, through the conversion of 1,3-bisphosphoglycerate (1,3BPGA) to 3-phosphoglycerate (3PGA). It also participates in gluconeogenesis, catalysing the opposite reaction to produce 1,3BPGA and ADP. Like most other glycolytic enzymes, PGK has also been catalogued as a moonlighting protein, due to its involvement in different functions not associated with energy metabolism, which include pathogenesis, interaction with nucleic acids, tumorigenesis progression, cell death and viral replication. In this review, we have highlighted the overall aspects of this enzyme, such as its structure, reaction kinetics, activity regulation and possible moonlighting functions in different protistan organisms, especially both free-living and parasitic Kinetoplastea. Our analysis of the genomes of different kinetoplastids revealed the presence of open-reading frames (ORFs) for multiple PGK isoforms in several species. Some of these ORFs code for unusually large PGKs. The products appear to contain additional structural domains fused to the PGK domain. A striking aspect is that some of these PGK isoforms are predicted to be catalytically inactive enzymes or ‘dead’ enzymes. The roles of PGKs in kinetoplastid parasites are analysed, and the apparent significance of the PGK gene duplication that gave rise to the different isoforms and their expression in Trypanosoma cruzi is discussed.

Design, synthesis and biological evaluation of novel 7H-benzo [c] [1, 3] dioxolo [4, 5-f] chromen-7-one derivatives with potential anti-tumor activity

In this study, a series of novel 7H-benzo [c] [1,3] dioxolo [4, 5-f] chromen-7-one derivatives were obtained by structural modification of the lead compounds with Fissitungfine B. A total 15 compounds were designed, synthesized and evaluated as inhibitors of tumor. These target compounds have the novel chemical structures that named three six-membered rings including one lactone six-membered ring. In vitro assay, the results showed that the target compounds have a broad spectrum and strong of anti-tumor activity. Such as the target compound 4n to MCF-7 was IC₅₀ = 0.35 ± 0.01 μM, to A-549 was IC₅₀ = 0.37 ± 0.01 μM, to Hela was IC₅₀ = 0.56 ± 0.02 μM, to MDC-803 was IC₅₀ = 0.53 ± 0.02 μM and COLO-205 was IC₅₀ = 0.50 ± 0.02 μM in vitro. At the same time, in vivo anti-tumor activity assay results showed that the target compounds had a good inhibitory effect on tumor growth. Among them, the target compound 4n had the best anti-tumor activity, it could inhibit tumor growth well at a low dose. The target compound 4n could be used as a candidate drug for further research and development, in order to be used as early as application in the clinical treatment of tumors.

Cataloguing the dead: breathing new life into pseudokinase research

Pseudoenzymes are present within many, but not all, known enzyme families and lack one or more conserved canonical amino acids that help define their catalytically active counterparts. Recent findings in the pseudokinase field confirm that evolutionary repurposing of the structurally defined bilobal protein kinase fold permits distinct biological functions to emerge, many of which rely on conformational switching, as opposed to canonical catalysis. In this analysis, we evaluate progress in evaluating several members of the 'dark' pseudokinome that are pertinent to help drive this expanding field. Initially, we discuss how adaptions in erythropoietin-producing hepatocellular carcinoma (Eph) receptor tyrosine kinase domains resulted in two vertebrate pseudokinases, EphA10 and EphB6, in which co-evolving sequences generate new motifs that are likely to be important for both nucleotide binding and catalysis-independent signalling. Secondly, we discuss how conformationally flexible Tribbles pseudokinases, which have radiated in the complex vertebrates, control fundamental aspects of cell signalling that may be targetable with covalent small molecules. Finally, we show how species-level adaptions in the duplicated canonical kinase protein serine kinase histone (PSKH)1 sequence have led to the appearance of the pseudokinase PSKH2, whose physiological role remains mysterious. In conclusion, we show how the patterns we discover are selectively conserved within specific pseudokinases, and that when they are modelled alongside closely related canonical kinases, many are found to be located in functionally important regions of the conserved kinase fold. Interrogation of these patterns will be useful for future evaluation of these, and other, members of the unstudied human kinome.

Moonlighting Proteins in the Fuzzy Logic of Cellular Metabolism

The numerous interconnected biochemical pathways that make up the metabolism of a living cell comprise a fuzzy logic system because of its high level of complexity and our inability to fully understand, predict, and model the many activities, how they interact, and their regulation. Each cell contains thousands of proteins with changing levels of expression, levels of activity, and patterns of interactions. Adding more layers of complexity is the number of proteins that have multiple functions. Moonlighting proteins include a wide variety of proteins where two or more functions are performed by one polypeptide chain. In this article, we discuss examples of proteins with variable functions that contribute to the fuzziness of cellular metabolism.

Genome-wide and structural analyses of pseudokinases encoded in the genome of Arabidopsis thaliana provide functional insights

Protein Kinase-Like Non-Kinases (PKLNKs), commonly known as "pseudokinases", are homologous to eukaryotic Ser/Thr/Tyr protein kinases (PKs) but lack the crucial aspartate residue in the catalytic loop, indispensable for phosphotransferase activity. Therefore, they are predicted to be "catalytically inactive" enzyme homologs. Analysis of protein-kinase like sequences from Arabidopsis thaliana led to the identification of more than 120 pseudokinases lacking catalytic aspartate, majority of which are closely related to the plant-specific receptor-like kinase family. These pseudokinases engage in different biological processes, enabled by their diverse domain architectures and specific subcellular localizations. Structural comparison of pseudokinases with active and inactive conformations of canonical PKs, belonging to both plant and animal origin, revealed unique structural differences. The currently available crystal structures of pseudokinases show that the loop topologically equivalent to activation segment of PKs adopts a distinct-folded conformation, packing against the pseudoenzyme core, in contrast to the extended and inhibitory geometries observed for active and inactive states, respectively, of catalytic PKs. Salt-bridge between ATP-binding Lys and DFG-Asp as well as hydrophobic interactions between the conserved nonpolar residue C-terminal to the equivalent DFG motif and nonpolar residues in C-helix mediate such a conformation in pseudokinases. This results in enhanced solvent accessibility of the pseudocatalytic loop in pseudokinases that can possibly serve as an interacting surface while associating with other proteins. Specifically, our analysis identified several residues that may be involved in pseudokinase regulation and hints at the repurposing of pseudocatalytic residues to achieve mechanistic control over noncatalytic functions of pseudoenzymes.

Enzymes, pseudoenzymes, and moonlighting proteins: diversity of function in protein superfamilies

As more genome sequences are elucidated, there is an increasing need for information about the functions of the millions of proteins they encode. The function of a newly sequenced protein is often estimated by sequence alignment with the sequences of proteins with known functions. However, protein superfamilies can contain members that share significant amino acid sequence and structural homology yet catalyze different reactions or act on different substrates. Some homologous proteins differ by having a second or even third function, called moonlighting proteins. More recently, it was found that most protein superfamilies also include pseudoenzymes, a protein, or a domain within a protein, that has a three-dimensional fold that resembles a conventional catalytically active enzyme, but has no catalytic activity. In this review, we discuss several examples of protein families that contain enzymes, pseudoenzymes, and moonlighting proteins. It is becoming clear that pseudoenzymes and moonlighting proteins are widespread in the evolutionary tree, and in many protein families, and they are often very similar in sequence and structure to their monofunctional and catalytically active counterparts. A greater understanding is needed to clarify when similarities and differences in amino acid sequences and structures correspond to similarities and differences in biochemical functions and cellular roles. This information can help improve programs that identify protein functions from sequence or structure and assist in more accurate annotation of sequence and structural databases, as well as in our understanding of the broad diversity of protein functions.

Pseudophosphatase MK-STYX: the atypical member of the MAP kinase phosphatases

The regulation of the phosphorylation of mitogen-activated protein kinases (MAPKs) is essential for cellular processes such as proliferation, differentiation, survival, and death. Mutations within the MAPK signaling cascades are implicated in diseases such as cancer, neurodegenerative disorders, arthritis, obesity, and diabetes. MAPK phosphorylation is controlled by an intricate balance between MAPK kinases (enzymes that add phosphate groups) and MAPK phosphatases (MKPs) (enzymes that remove phosphate groups). MKPs are complex negative regulators of the MAPK pathway that control the amplitude and spatiotemporal regulation of MAPKs. MK-STYX (MAPK phosphoserine/threonine/tyrosine-binding protein) is a member of the MKP subfamily, which lacks the critical histidine and nucleophilic cysteine residues in the active site required for catalysis. MK-STYX does not influence the phosphorylation status of MAPK, but even so it adds to the complexity of signal transduction cascades as a signaling regulator. This review highlights the function of MK-STYX, providing insight into MK-STYX as a signal regulating molecule in the stress response, HDAC 6 dynamics, apoptosis, and neurite differentiation.

Leitmotif: protein motif scanning 2.0

MOTIVATION

Motif-HMM (mHMM) scanning has been shown to possess unique advantages over standardly used sequence-profile search methods (e.g. HMMER, PSI-BLAST) since it is particularly well-suited to discriminate proteins with variations inside conserved motifs (e.g. family subtypes) or motifs lacking essential residues (false positives, e.g. pseudoenzymes).

RESULTS

In order to make mHMM widely accessible to a broader scientific community, we developed Leitmotif, an mHMM web application with many parametrization options easily accessible through intuitive interface. Substantial improvement of performance (ROC scores) was obtained by using two novel parameters. To the best of our knowledge, Leitmotif is the only available mHMM application.

AVAILABILITY AND IMPLEMENTATION

Leitmotif is freely available at https://leitmotif.irb.hr.

CONTACT

sinisa@heuristika.hr or ivan.vujaklija@fer.hr.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Differential effects of 'resurrecting' Csp pseudoproteases during Clostridioides difficile spore germination

Clostridioides difficile is a spore-forming bacterial pathogen that is the leading cause of hospital-acquired gastroenteritis. C. difficile infections begin when its spore form germinates in the gut upon sensing bile acids. These germinants induce a proteolytic signaling cascade controlled by three members of the subtilisin-like serine protease family, CspA, CspB, and CspC. Notably, even though CspC and CspA are both pseudoproteases, they are nevertheless required to sense germinants and activate the protease, CspB. Thus, CspC and CspA are part of a growing list of pseudoenzymes that play important roles in regulating cellular processes. However, despite their importance, the structural properties of pseudoenzymes that allow them to function as regulators remain poorly understood. Our recently solved crystal structure of CspC revealed that its pseudoactive site residues align closely with the catalytic triad of CspB, suggesting that it might be possible to ‘resurrect' the ancestral protease activity of the CspC and CspA pseudoproteases. Here, we demonstrate that restoring the catalytic triad to these pseudoproteases fails to resurrect their protease activity. We further show that the pseudoactive site substitutions differentially affect the stability and function of the CspC and CspA pseudoproteases: the substitutions destabilized CspC and impaired spore germination without affecting CspA stability or function. Thus, our results surprisingly reveal that the presence of a catalytic triad does not necessarily predict protease activity. Since homologs of C. difficile CspA occasionally carry an intact catalytic triad, our results indicate that bioinformatic predictions of enzyme activity may underestimate pseudoenzymes in rare cases.

Crystal structure of DRIK1, a stress-responsive receptor-like pseudokinase, reveals the molecular basis for the absence of ATP binding

BACKGROUND

Plants reprogram metabolism and development to rapidly adapt to biotic and abiotic stress. Protein kinases play a significant role in this process by phosphorylating protein substrates that activate or inactivate signaling cascades that regulate cellular and metabolic adaptations. Despite their importance in plant biology, a notably small fraction of the plant kinomes has been studied to date.

RESULTS

In this report, we describe ZmDRIK1, a stress-responsive receptor-like pseudokinase whose expression is downregulated under water restriction. We show the structural features and molecular basis of the absence of ATP binding exhibited by ZmDRIK1. The ZmDRIK1 kinase domain lacks conserved amino acids that are essential for phosphorylation activity. The crystal structure of the ZmDRIK1 kinase domain revealed the presence of a spine formed by the side chain of the triad Leu²⁴⁰, Tyr³⁶³, and Leu³⁷⁵ that occludes the ATP binding pocket. Although ZmDRIK1 is unable to bind nucleotides, it does bind the small molecule ENMD-2076 which, in a cocrystal structure, revealed the potential to serve as a ZmDRIK1 inhibitor.

CONCLUSION

ZmDRIK1 is a novel receptor-like pseudokinase responsive to biotic and abiotic stress. The absence of ATP binding and consequently, the absence of phosphorylation activity, was proven by the crystal structure of the apo form of the protein kinase domain. The expression profiling of the gene encoding ZmDRIK1 suggests this kinase may play a role in downregulating the expression of stress responsive genes that are not necessary under normal conditions. Under biotic and abiotic stress, ZmDRIK1 is down-regulated to release the expression of these stress-responsive genes.

Structural Basis of Protein Arginine Methyltransferase Activation by a Catalytically Dead Homolog (Prozyme)

Prozymes are pseudoenzymes that stimulate the function of weakly active enzymes through complex formation. The major Trypanosoma brucei protein arginine methyltransferase, TbPRMT1 enzyme (ENZ), requires TbPRMT1 prozyme (PRO) to form an active heterotetrameric complex. Here, we present the X-ray crystal structure of the TbPRMT1 ENZ-Δ52PRO tetrameric complex with the cofactor product S-adenosyl-l-homocysteine (AdoHcy) at 2.4 Å resolution. The individual ENZ and PRO units adopt the highly-conserved PRMT domain architecture and form an antiparallel heterodimer that corresponds to the canonical homodimer observed in all previously reported PRMTs. In turn, two such heterodimers assemble into a tetramer both in the crystal and in solution with twofold rotational symmetry. ENZ is unstable in absence of PRO and incapable of forming a homodimer due to a steric clash of an ENZ-specific tyrosine within the dimerization arm, rationalizing why PRO is required to complement ENZ to form a PRMT dimer that is necessary, but not sufficient for PRMT activity. The PRO structure deviates from other, active PRMTs in that it lacks the conserved η2 3₁₀-helix within the Rossmann fold, abolishing cofactor binding. In addition to its chaperone function for ENZ, PRO substantially contributes to substrate binding. Heterotetramerization is required for catalysis, as heterodimeric ENZ-PRO mutants lack binding affinity and methyltransferase activity toward the substrate protein TbRGG1. Together, we provide a structural basis for TbPRMT1 ENZ activation by PRO heterotetramer formation, which is conserved across all kinetoplastids, and describe a chaperone function of the TbPRMT1 prozyme, which represents a novel mode of PRMT regulation.

Single-nucleotide polymorphism screening and RNA sequencing of key messenger RNAs associated with neonatal hypoxic-ischemia brain damage

A single-nucleotide polymorphism (SNP) is an alteration in one nucleotide in a certain position within a genome. SNPs are associated with disease susceptibility. However, the influences of SNPs on the pathogenesis of neonatal hypoxic-ischemic brain damage remain elusive. Seven-day-old rats were used to establish a hypoxic ischemic encephalopathy model. SNPs and expression profiles of mRNAs were analyzed in hypoxic ischemic encephalopathy model rats using RNA sequencing. Genes exhibiting SNPs associated with hypoxic ischemic encephalopathy were identified and studied by gene ontology and pathway analysis to identify their possible involvement in the disease mechanism. We identified 89 up-regulated genes containing SNPs that were mainly located on chromosome 1 and 2. Gene ontology analysis indicated that the up-regulated genes containing SNPs are mainly involved in angiogenesis, wound healing and glutamatergic synapse and biological processing of calcium-activated chloride channels. Signaling pathway analysis indicated that the differentially expressed genes play a role in glutamatergic synapses, long-term depression and oxytocin signaling. Moreover, intersection analysis of high throughput screening following PubMed retrieval and RNA sequencing for SNPs showed that CSRNP1, DUSP5 and LRRC25 were most relevant to hypoxic ischemic encephalopathy. Significant up-regulation of genes was confirmed by quantitative real-time polymerase chain reaction analysis of oxygen-glucose-deprived human fetal cortical neurons. Our results indicate that CSRNP1, DUSP5 and LRRC25, containing SNPs, may be involved in the pathogenesis of hypoxic ischemic encephalopathy. These findings indicate a novel direction for further hypoxic ischemic encephalopathy research. This animal study was approved on February 5, 2017 by the Animal Care and Use Committee of Kunming Medical University, Yunnan Province, China (approval No. kmmu2019038). Cerebral tissue collection from a human fetus was approved on September 30, 2015 by the Ethics Committee of Kunming Medical University, China (approval No. 2015-9).

Clarification of the dispensability of PDX1.2 for Arabidopsis viability using CRISPR/Cas9

BACKGROUND

PDX1.2 has recently been shown to be a regulator of vitamin B₆ biosynthesis in plants and is implicated in biotic and abiotic stress resistance. PDX1.2 expression is strongly and rapidly induced by heat stress. Interestingly, PDX1.2 is restricted to eudicota, wherein it behaves as a non-catalytic pseudoenzyme and is suggested to provide an adaptive advantage to this clade. A first report on an Arabidopsis insertion mutant claims that PDX1.2 is indispensable for viability, being essential for embryogenesis. However, a later study using an independent insertion allele suggests that knockout mutants of pdx1.2 are viable. Therefore, the essentiality of PDX1.2 for Arabidopsis viability is a matter of debate. Given the important implications of PDX1.2 in stress responses, it is imperative to clarify if it is essential for plant viability.

RESULTS

We have studied the previously reported insertion alleles of PDX1.2, one of which is claimed to be essential for embryogenesis (pdx1.2-1), whereas the other is viable (pdx1.2-2). Our study shows that pdx1.2-1 carries multiple T-DNA insertions, but the T-DNA insertion in PDX1.2 is not responsible for the loss of embryogenesis. By contrast, the pdx1.2-2 allele is an overexpressor of PDX1.2 under standard growth conditions and not a null allele as previously reported. Nonetheless, upregulation of PDX1.2 expression under heat stress is impaired in this mutant line. In wild type Arabidopsis, studies of PDX1.2-YFP fusion proteins show that the protein is enhanced under heat stress conditions. To clarify if PDX1.2 is essential for Arabidopsis viability, we generated several independent mutant lines using the CRISPR-Cas9 gene editing technology. All of these lines are viable and behave similar to wild type under standard growth conditions. Reciprocal crosses of a subset of the CRISPR lines with pdx1.2-1 recovers viability of the latter line and demonstrates that knocking out the functionality of PDX1.2 does not impair embryogenesis.

CONCLUSIONS

Gene editing reveals that PDX1.2 is dispensable for Arabidopsis viability and resolves conflicting reports in the literature on its function.

Emerging concepts in pseudoenzyme classification, evolution, and signaling

The 21st century is witnessing an explosive surge in our understanding of pseudoenzyme-driven regulatory mechanisms in biology. Pseudoenzymes are proteins that have sequence homology with enzyme families but that are proven or predicted to lack enzyme activity due to mutations in otherwise conserved catalytic amino acids. The best-studied pseudoenzymes are pseudokinases, although examples from other families are emerging at a rapid rate as experimental approaches catch up with an avalanche of freely available informatics data. Kingdom-wide analysis in prokaryotes, archaea and eukaryotes reveals that between 5 and 10% of proteins that make up enzyme families are pseudoenzymes, with notable expansions and contractions seemingly associated with specific signaling niches. Pseudoenzymes can allosterically activate canonical enzymes, act as scaffolds to control assembly of signaling complexes and their localization, serve as molecular switches, or regulate signaling networks through substrate or enzyme sequestration. Molecular analysis of pseudoenzymes is rapidly advancing knowledge of how they perform noncatalytic functions and is enabling the discovery of unexpected, and previously unappreciated, functions of their intensively studied enzyme counterparts. Notably, upon further examination, some pseudoenzymes have previously unknown enzymatic activities that could not have been predicted a priori. Pseudoenzymes can be targeted and manipulated by small molecules and therefore represent new therapeutic targets (or anti-targets, where intervention should be avoided) in various diseases. In this review, which brings together broad bioinformatics and cell signaling approaches in the field, we highlight a selection of findings relevant to a contemporary understanding of pseudoenzyme-based biology.

Nature of Long-Range Evolutionary Constraint in Enzymes: Insights from Comparison to Pseudoenzymes with Similar Structures

Enzymes are known to fine-tune their sequences to optimize catalytic function, yet quantitative evolutionary design principles of enzymes remain elusive on the proteomic scale. Recently, it was found that the catalytic site in enzymes induces long-range evolutionary constraint, where even sites distant to the catalytic site are more conserved than expected. Given that protein-fold usage is generally different between enzymes and nonenzymes, it remains an open question to what extent this long-range evolutionary constraint in enzymes is dictated, either directly or indirectly, by the special three-dimensional structure of the enzyme. To investigate this question, we have compared evolutionary properties of enzymes with those of counterpart pseudoenzymes that share the same protein fold but are catalytically inactive. We found that the long-range evolutionary constraint observed in enzymes is significantly reduced in pseudoenzyme counterparts, despite very high structural similarity (∼1.5 Å RMSD on average). Furthermore, this significant reduction in long-range evolutionary constraint is observed even in pseudoenzyme counterparts which retain the ligand-binding ability of enzymes. Finally, the distance between the site that induces the highest gradient of sequence conservation and the pseudocatalytic site in pseudoenzymes is significantly larger than the corresponding distance in enzymes. Taken together, our results suggest that the long-range evolutionary constraint in enzymes is induced mainly by the presence of the catalytic site rather than by the special three-dimensional structure of the enzyme, and that such long-range evolutionary constraint in enzymes depends mainly on the catalytic function of the active site rather than on the ligand-binding ability of the enzyme.

Using Pseudoenzymes to Probe Evolutionary Design Principles of Enzymes

Enzymes are governed by unique evolutionary design principles as their catalytic sites were shown to induce long-range evolutionary conservation gradients. We have recently used a comparative bioinformatics approach to disentangle structural determinants from other possible determinants of the evolutionary conservation gradients. The approach is based on comparing the evolutionary patterns of enzymes to those of pseudoenzymes with the same tertiary structure where the catalytic functionality is turned off. This approach provides a way to evaluate several hypotheses regarding the origin of the observed evolutionary conservation gradient in enzymes. The conclusions from such comparative analyses are important for a better understanding of the unique evolutionary design principles of enzymes, which can in turn potentially guide the design of new and improved enzymes.

Human pantothenate kinase 4 is a pseudo-pantothenate kinase

Pantothenate kinase generates 4'-phosphopantothenate in the first and rate-determining step of coenzyme A (CoA) biosynthesis. The human genome encodes three well-characterized and nearly identical pantothenate kinases (PANK1-3) plus a putative bifunctional protein (PANK4) with a predicted amino-terminal pantothenate kinase domain fused to a carboxy-terminal phosphatase domain. Structural and phylogenetic analyses show that all active, characterized PANKs contain the key catalytic residues Glu138 and Arg207 (HsPANK3 numbering). However, all amniote PANK4s, including human PANK4, encode Glu138Val and Arg207Trp substitutions which are predicted to inactivate kinase activity. Biochemical analysis corroborates bioinformatic predictions-human PANK4 lacks pantothenate kinase activity. Introducing Glu138Val and Arg207Trp substitutions to the human PANK3 and plant PANK4 abolished their robust pantothenate kinase activity. Introducing both catalytic residues back into human PANK4 restored kinase activity, but only to a low level. This result suggests that epistatic changes to the rest of the protein already reduced the kinase activity prior to mutation of the catalytic residues in the course of evolution. The PANK4 from frog, an anamniote living relative encoding the catalytically active residues, had only a low level of kinase activity, supporting the view that HsPANK4 had reduced kinase activity prior to the catalytic residue substitutions in amniotes. Together, our data show that human PANK4 is a pseudo-pantothenate kinase-a catalytically deficient variant of the catalytically active PANK4 found in plants and fungi. The Glu138Val and Arg207Trp substitutions in amniotes (HsPANK3 numbering) completely deactivated the pantothenate kinase activity that had already been reduced by prior epistatic mutations.

Tracing the origin and evolution of pseudokinases across the tree of life

Protein phosphorylation by eukaryotic protein kinases (ePKs) is a fundamental mechanism of cell signaling in all organisms. In model vertebrates, ~10% of ePKs are classified as pseudokinases, which have amino acid changes within the catalytic machinery of the kinase domain that distinguish them from their canonical kinase counterparts. However, pseudokinases still regulate various signaling pathways, usually doing so in the absence of their own catalytic output. To investigate the prevalence, evolutionary relationships, and biological diversity of these pseudoenzymes, we performed a comprehensive analysis of putative pseudokinase sequences in available eukaryotic, bacterial, and archaeal proteomes. We found that pseudokinases are present across all domains of life, and we classified nearly 30,000 eukaryotic, 1500 bacterial, and 20 archaeal pseudokinase sequences into 86 pseudokinase families, including ~30 families that were previously unknown. We uncovered a rich variety of pseudokinases with notable expansions not only in animals but also in plants, fungi, and bacteria, where pseudokinases have previously received cursory attention. These expansions are accompanied by domain shuffling, which suggests roles for pseudokinases in plant innate immunity, plant-fungal interactions, and bacterial signaling. Mechanistically, the ancestral kinase fold has diverged in many distinct ways through the enrichment of unique sequence motifs to generate new families of pseudokinases in which the kinase domain is repurposed for noncanonical nucleotide binding or to stabilize unique, inactive kinase conformations. We further provide a collection of annotated pseudokinase sequences in the Protein Kinase Ontology (ProKinO) as a new mineable resource for the signaling community.

Phylogenetic Distribution and Diversity of Bacterial Pseudo-Orthocaspases Underline Their Putative Role in Photosynthesis

Orthocaspases are prokaryotic caspase homologs - proteases, which cleave their substrates after positively charged residues using a conserved histidine - cysteine (HC) dyad situated in a catalytic p20 domain. However, in orthocaspases pseudo-variants have been identified, which instead of the catalytic HC residues contain tyrosine and serine, respectively. The presence and distribution of these presumably proteolytically inactive p20-containing enzymes has until now escaped attention. We have performed a detailed analysis of orthocaspases in all available prokaryotic genomes, focusing on pseudo-orthocaspases. Surprisingly we identified type I metacaspase homologs in filamentous cyanobacteria. While genes encoding pseudo-orthocaspases seem to be absent in Archaea, our results show conservation of these genes in organisms performing either anoxygenic photosynthesis (orders Rhizobiales, Rhodobacterales, and Rhodospirillales in Alphaproteobacteria) or oxygenic photosynthesis (all sequenced cyanobacteria, except Gloeobacter, Prochlorococcus, and Cyanobium). Contrary to earlier reports, we were able to detect pseudo-orthocaspases in all sequenced strains of the unicellular cyanobacteria Synechococcus and Synechocystis. In silico comparisons of the primary as well as tertiary structures of pseudo-p20 domains with their presumably proteolytically active homologs suggest that differences in their amino acid sequences have no influence on the overall structures. Mutations therefore affect most likely only the proteolytic activity. Our data provide an insight into diversification of pseudo-orthocaspases in Prokaryotes, their taxa-specific distribution, and allow suggestions on their taxa-specific function.

The fungal ribonuclease-like effector protein CSEP0064/BEC1054 represses plant immunity and interferes with degradation of host ribosomal RNA

The biotrophic fungal pathogen Blumeria graminis causes the powdery mildew disease of cereals and grasses. We present the first crystal structure of a B. graminis effector of pathogenicity (CSEP0064/BEC1054), demonstrating it has a ribonuclease (RNase)-like fold. This effector is part of a group of RNase-like proteins (termed RALPHs) which comprise the largest set of secreted effector candidates within the B. graminis genomes. Their exceptional abundance suggests they play crucial functions during pathogenesis. We show that transgenic expression of RALPH CSEP0064/BEC1054 increases susceptibility to infection in both monocotyledonous and dicotyledonous plants. CSEP0064/BEC1054 interacts in planta with the pathogenesis-related protein PR10. The effector protein associates with total RNA and weakly with DNA. Methyl jasmonate (MeJA) levels modulate susceptibility to aniline-induced host RNA fragmentation. In planta expression of CSEP0064/BEC1054 reduces the formation of this RNA fragment. We propose CSEP0064/BEC1054 is a pseudoenzyme that binds to host ribosomes, thereby inhibiting the action of plant ribosome-inactivating proteins (RIPs) that would otherwise lead to host cell death, an unviable interaction and demise of the fungus.

Powdery mildews are common plant diseases which affect important crop plants including cereals such as wheat and barley. The fungi that cause this disease are obligate biotrophs: they have an absolute requirement for living host cells which they penetrate with feeding structures called haustoria. These fungi must be highly effective at avoiding immune recognition which would lead to death of the host cell and the pathogen. We assume they do this by delivering effector proteins to the host. While several hundred secreted effectors have been described in cereal powdery mildews, it is unknown how they work. Here, we use X-ray crystallography and nuclear magnetic resonance (NMR) spectroscopy to determine the structure and interactions of the effector CSEP0064/BEC1054, representative of the largest class of effectors resembling fungal RNases. We find that this effector binds nucleic acids. Expression of the effector in plants increases susceptibility to infection. Moreover, transgenic CSEP0064/BEC1054 expression in wheat inhibits the degradation of host ribosomal RNA induced by ribosome-inactivating proteins (RIPs). We propose a novel mechanism of action for the RNase-like effectors in powdery mildews: they may act as pseudoenzymes to inhibit the host RIPs, known components of plant immune responses that lead to host cell death.

The demise of catalysis, but new functions arise: pseudoenzymes as the phoenixes of the protein world

Pseudoenzymes are noncatalytic homologues of enzymes and are found in most enzyme families. Although lacking catalytic activity and sometimes referred to as 'dead' enzymes, they instead resemble phoenixes because the loss of a catalytic function during evolution was associated with the development of vital new functions. They are important in regulating the activity and location of catalytically active homologues, scaffolding the assembly of signaling complexes, and regulating transcription or translation. They are key actors in cell proliferation and differentiation, proteostasis, and many other biochemical pathways and processes. They perform their functions in diverse ways, but many retain some aspects of the function of their catalytically active homologues. In some pseudoenzymes, their functions are very different from other members of their protein families, suggesting some arose from ancient moonlighting proteins during evolution. Much less is known about pseudoenzymes than their catalytically active counterparts, but a growing appreciation of their key roles in many important biochemical processes and signaling pathways has led to increased investigation in recent years. It is clear that there is still much more to learn about the structures, functions, and cellular roles of these phoenix-like proteins.

The HDAC3 enzymatic activity regulates skeletal muscle fuel metabolism

Histone deacetylase 3 (HDAC3) is a major HDAC, whose enzymatic activity is targeted by small molecule inhibitors for treating a variety of conditions. However, its enzymatic activity is largely dispensable for its function in embryonic development and hepatic lipid metabolism. HDAC3 plays a pivotal role in regulating muscle fuel metabolism and contractile function. Here, we address whether these muscular functions of HDAC3 require its enzymatic activity. By mutating the NCoR/SMRT corepressors in a knock-in mouse model named NS-DADm, we ablated the enzymatic activity of HDAC3 without affecting its protein levels. Compared to the control mice, skeletal muscles from NS-DADm mice showed lower force generation, enhanced fatigue resistance, enhanced fatty acid oxidation, reduced glucose uptake during exercise, upregulated expression of metabolic genes involved in branched-chain amino acids catabolism, and reduced muscle mass during aging, without changes in the muscle fiber-type composition or mitochondrial protein content. These muscular phenotypes are similar to those observed in the HDAC3-depleted skeletal muscles, which demonstrates that, unlike that in the liver or embryonic development, the metabolic function of HDAC3 in skeletal muscles requires its enzymatic activity. These results suggest that drugs specifically targeting HDAC3 enzyme activity could be developed and tested to modulate muscle energy metabolism and exercise performance.

The pseudokinase domains of guanylyl cyclase-A and -B allosterically increase the affinity of their catalytic domains for substrate

Natriuretic peptides regulate multiple physiologic systems by activating transmembrane receptors containing intracellular guanylyl cyclase domains, such as GC-A and GC-B, also known as Npr1 and Npr2, respectively. Both enzymes contain an intracellular, phosphorylated pseudokinase domain (PKD) critical for activation of the C-terminal cGMP-synthesizing guanylyl cyclase domain. Because ATP allosterically activates GC-A and GC-B, we investigated how ATP binding to the PKD influenced guanylyl cyclase activity. Molecular modeling indicated that all the residues of the ATP-binding site of the prototypical kinase PKA, except the catalytic aspartate, are conserved in the PKDs of GC-A and GC-B. Kinase-inactivating alanine substitutions for the invariant lysine in subdomain II or the aspartate in the DYG-loop of GC-A and GC-B failed to decrease enzyme phosphate content, consistent with the PKDs lacking kinase activity. In contrast, both mutations reduced enzyme activation by blocking the ability of ATP to decrease the Michaelis constant without affecting peptide-dependent activation. The analogous lysine-to-alanine substitution in a glutamate-substituted phosphomimetic mutant form of GC-B also reduced enzyme activity, consistent with ATP stimulating guanylyl cyclase activity through an allosteric, phosphorylation-independent mechanism. Mutations designed to rigidify the conserved regulatory or catalytic spines within the PKDs increased guanylyl cyclase activity, increased sensitivity to natriuretic peptide, or reduced the Michaelis constant in the absence of ATP, consistent with ATP binding stabilizing the PKD in a conformation analogous to that of catalytically active kinases. We conclude that allosteric mechanisms evolutionarily conserved in the PKDs promote the catalytic activation of transmembrane guanylyl cyclases.

Novel pseudo-aspartic peptidase from the midgut of the tick Rhipicephalus microplus

The characterization of Rhipicephalus microplus tick physiology can support efforts to develop and improve the efficiency of control methods. A sequence containing a domain with similarity to one derived from the aspartic peptidase family was isolated from the midgut of engorged female R. microplus. The lack of the second catalytic aspartic acid residue suggest that it may be a pseudo-aspartic peptidase, and it was named RmPAP. In this work we confirm the lack of proteolytic activity of RmPAP and investigate it’s non-proteolytic interaction with bovine hemoglobin by Surface Plasmon Resonance and phage display. Moreover we carried out RNAi interference and artificial feeding of ticks with anti-RmPAP antibodies to assess it’s possible biological role, although no changes were observed in the biological parameters evaluated. Overall, we hypothesize that RmPAP may act as a carrier of hemoglobin/heme between the tick midgut and the ovaries.

PseudoGTPase domains in p190RhoGAP proteins: a mini-review

Pseudoenzymes generally lack detectable catalytic activity despite adopting the overall protein fold of their catalytically competent counterparts, indeed 'pseudo' family members seem to be incorporated in all enzyme classes. The small GTPase enzymes are important signaling proteins, and recent studies have identified many new family members with noncanonical residues within the catalytic cleft, termed pseudoGTPases. To illustrate recent discoveries in the field, we use the p190RhoGAP proteins as an example. p190RhoGAP proteins (ARHGAP5 and ARHGAP35) are the most abundant GTPase activating proteins for the Rho family of small GTPases. These are key regulators of Rho signaling in processes such as cell migration, adhesion and cytokinesis. Structural biology has complemented and guided biochemical analyses for these proteins and has allowed discovery of two cryptic pseudoGTPase domains, and the re-classification of a third, previously identified, GTPase-fold domain as a pseudoGTPase. The three domains within p190RhoGAP proteins illustrate the diversity of this rapidly expanding pseudoGTPase group.

The N-Terminal GTPase Domain of p190RhoGAP Proteins Is a PseudoGTPase

The pseudoGTPases are a rapidly growing and important group of pseudoenzymes. p190RhoGAP proteins are critical regulators of Rho signaling and contain two previously identified pseudoGTPase domains. Here we report that p190RhoGAP proteins contain a third pseudoGTPase domain, termed N-GTPase. We find that GTP constitutively purifies with the N-GTPase domain, and a 2.8-Å crystal structure of p190RhoGAP-A co-purified with GTP reveals an unusual GTP-Mg2+ binding pocket. Six inserts in N-GTPase indicate perturbed catalytic activity and inability to bind to canonical GTPase activating proteins, guanine nucleotide exchange factors, and effector proteins. Biochemical analysis shows that N-GTPase does not detectably hydrolyze GTP, and exchanges nucleotide only under harsh Mg2+ chelation. Furthermore, mutational analysis shows that GTP and Mg2+ binding stabilizes the domain. Therefore, our results support that N-GTPase is a nucleotide binding, non-hydrolyzing, pseudoGTPase domain that may act as a protein-protein interaction domain. Thus, unique among known proteins, p190RhoGAPs contain three pseudoGTPase domains.