Search

An evolutionary-based definition and classification of target evaluation units (EUs) is presented for the 14th round of the critical assessment of structure prediction (CASP14). CASP14 targets included 84 experimental models submitted by various structural groups (designated T1024-T1101). Targets were split into EUs based on the domain organization of available templates and performance of server groups. Several targets required splitting (19 out of 25 multidomain targets) due in part to observed conformation changes. All in all, 96 CASP14 EUs were defined and assigned to tertiary structure assessment categories (Topology-based FM or High Accuracy-based TBM-easy and TBM-hard) considering their evolutionary relationship to existing ECOD fold space: 24 family level, 50 distant homologs (H-group), 12 analogs (X-group), and 10 new folds. Principal component analysis and heatmap visualization of sequence and structure similarity to known templates as well as performance of servers highlighted trends in CASP14 target difficulty. The assigned evolutionary levels (i.e., H-groups) and assessment classes (i.e., FM) displayed overlapping clusters of EUs. Many viral targets diverged considerably from their template homologs and thus were more difficult for prediction than other homology-related targets. On the other hand, some targets did not have sequence-identifiable templates, but were predicted better than expected due to relatively simple arrangements of secondary structural elements. An apparent improvement in overall server performance in CASP14 further complicated traditional classification, which ultimately assigned EUs into high-accuracy modeling (27 TBM-easy and 31 TBM-hard), topology (23 FM), or both (15 FM/TBM).

This report describes the tertiary structure prediction assessment of difficult modeling targets in the 14(th) round of the Critical Assessment of Structure Prediction (CASP14). We implemented an official ranking scheme that used the same scores as the previous CASP topology-based assessment, but combined these scores with one that emphasized physically realistic models. The top performing AlphaFold2 group outperformed the rest of the prediction community on all but two of the difficult targets considered in this assessment. They provided high quality models for most of the targets (86% over GDT_TS 70), including larger targets above 150 residues, and they correctly predicted the topology of almost all the rest. AlphaFold2 performance was followed by two manual Baker methods, a Feig method that refined Zhang-server models, two notable automated Zhang server methods (QUARK and Zhang-server), and a Zhang manual group. Despite the remarkable progress in protein structure prediction of difficult targets, both the prediction community and AlphaFold2, to a lesser extent, faced challenges with flexible regions and obligate oligomeric assemblies. The official ranking of top-performing methods was supported by performance generated PCA and heatmap clusters that gave insight into target difficulties and the most successful state-of-the-art structure prediction methodologies.

The autosomal dominant disorder Schnyder corneal dystrophy (SCD) is caused by mutations in UbiA prenyltransferase domain-containing protein-1 (UBIAD1), which uses geranylgeranyl pyrophosphate (GGpp) to synthesize the vitamin K(2) subtype menaquinone-4 (MK-4). SCD is characterized by opacification of the cornea, owing to aberrant build-up of cholesterol in the tissue. We previously discovered that sterols stimulate association of UBIAD1 with ER-localized HMG-CoA reductase, which catalyzes a rate-limiting step in the synthesis of cholesterol and nonsterol isoprenoids, including GGpp. Binding to UBIAD1 inhibits sterol-accelerated ER-associated degradation (ERAD) of reductase and permits continued synthesis of GGpp in cholesterol-replete cells. GGpp disrupts UBIAD1-reductase binding and thereby allows for maximal ERAD of reductase as well as ER-to-Golgi translocation of UBIAD1. SCD-associated UBIAD1 is refractory to GGpp-mediated dissociation from reductase and remains sequestered in the ER to inhibit ERAD. Here, we report development of a biochemical assay for UBIAD1-mediated synthesis of MK-4 in isolated membranes and intact cells. Using this assay, we compared enzymatic activity of WT UBIAD1 with that of SCD-associated variants. Our studies revealed that SCD-associated UBIAD1 exhibited reduced MK-4 synthetic activity, which may result from its reduced affinity for GGpp. Sequestration in the ER protects SCD-associated UBIAD1 from autophagy and allows intracellular accumulation of the mutant protein, which amplifies the inhibitory effect on reductase ERAD. These findings have important implications not only for the understanding of SCD etiology but also for the efficacy of cholesterol-lowering statin therapy, which becomes limited, in part, because of UBIAD1-mediated inhibition of reductase ERAD.

Protein target structures for the Critical Assessment of Structure Prediction round 13 (CASP13) were split into evaluation units (EUs) based on their structural domains, the domain organization of available templates, and the performance of servers on whole targets compared to split target domains. Eighty targets were split into 112 EUs. The EUs were classified into categories suitable for assessment of high accuracy modeling (or template-based modeling [TBM]) and topology (or free modeling [FM]) based on target difficulty. Assignment into assessment categories considered the following criteria: (a) the evolutionary relationship of target domains to existing fold space as defined by the Evolutionary Classification of Protein Domains (ECOD) database; (b) the clustering of target domains using eight objective sequence, structure, and performance measures; and (c) the placement of target domains in a scatter plot of target difficulty against server performance used in the previous CASP. Generally, target domains with good server predictions had close template homologs and were classified as TBM. Alternately, targets with poor server predictions represent a mixture of fast evolving homologs, structure analogs, and new folds, and were classified as FM or FM/TBM overlap.

Adopting a proper topology is crucial for transmembrane proteins to perform their functions. We previously reported that ceramide regulates a transmembrane protein called TM4SF20 (transmembrane 4 L six family member 20) through topological inversion by altering the direction through which the protein is translocated across membranes during translation. This regulatory mechanism, denoted regulated alternative translocation (RAT), depends on a GXXXN motif present in the first transmembrane helix of TM4SF20. Here, using site-directed mutagenesis, we show that Asn-26 in the motif is crucial for RAT of TM4SF20, as it cannot be replaced even by Gln. In contrast, Gly-22 in the motif could be substituted by other small residues such as Ala and Ser without affecting RAT of TM4SF20. We further demonstrate that the GXXXN motif alone is insufficient to induce RAT of a transmembrane protein because TM4SF4, a relative of TM4SF20 that also contains the motif in the first transmembrane helix, did not undergo RAT. Using TM4SF40–TM4SF20 chimeras, we identified Pro-29 of TM4SF20 as another important element required for RAT of the protein. Substituting Pro-29 alone did not affect RAT of TM4SF20, whereas replacing Pro-29 together with either Leu-25 or Val-17 of TM4SF20 with the corresponding residues of TM4SF4 abolished RAT of TM4SF20. Because Val-17, Gly-22, Leu-25, Asn-26, and Pro-29 are predicted to reside along the same surface of the transmembrane helix, our results suggest that interactions with other proteins mediated by this surface during translocation may be critical for RAT of TM4SF20.

Upon phagocytosis into macrophages, the intracellular bacterial pathogen Legionella pneumophila secretes effector proteins that manipulate host cell components, enabling it to evade lysosomal degradation. However, the bacterial proteins involved in this evasion are incompletely characterized. Here we show that the L. pneumophila effector protein RavD targets host membrane compartments and contributes to the molecular mechanism the pathogen uses to prevent encounters with lysosomes. Protein–lipid binding assays revealed that RavD selectively binds phosphatidylinositol-3-phosphate (PI(3)P) in vitro. We further determined that a C-terminal RavD region mediates the interaction with PI(3)P and that this interaction requires Arg-292. In transiently transfected mammalian cells, mCherry-RavD colocalized with the early endosome marker EGFP-Rab5 as well as the PI(3)P biosensor EGFP-2×FYVE. However, treatment with the phosphoinositide 3-kinase inhibitor wortmannin did not disrupt localization of mCherry-RavD to endosomal compartments, suggesting that RavD's interaction with PI(3)P is not necessary to anchor RavD to endosomal membranes. Using superresolution and immunogold transmission EM, we observed that, upon translocation into macrophages, RavD was retained onto the Legionella-containing vacuole and was also present on small vesicles adjacent to the vacuole. We also report that despite no detectable effects on intracellular growth of L. pneumophila within macrophages or amebae, the lack of RavD significantly increased the number of vacuoles that accumulate the late endosome/lysosome marker LAMP-1 during macrophage infection. Together, our findings suggest that, although not required for intracellular replication of L. pneumophila, RavD is a part of the molecular mechanism that steers the Legionella-containing vacuole away from endolysosomal maturation pathways.

Vibrio parahaemolyticus is a seafood-borne pathogen that encodes two type 3 secretion systems (T3SS). The first system, T3SS1, is thought to be maintained in all strains of V. parahaemolyticus to maintain survival in the environment, whereas the second system, T3SS2, is linked to clinical isolates and disease in humans., Diverse bacterial pathogens employ effector delivery systems to disrupt vital cellular processes in the host (N. M. Alto and K. Orth, Cold Spring Harbor Perspect Biol 4:a006114, 2012, https://doi.org/10.1101/cshperspect.a006114). The type III secretion system 1 of the marine pathogen Vibrio parahaemolyticus utilizes the sequential action of four effectors to induce a rapid, proinflammatory cell death uniquely characterized by a prosurvival host transcriptional response (D. L. Burdette, M. L. Yarbrough, A Orvedahl, C. J. Gilpin, and K. Orth, Proc Natl Acad Sci USA 105:12497–12502, 2008, https://doi.org/10.1073/pnas.0802773105; N. J. De Nisco, M. Kanchwala, P. Li, J. Fernandez, C. Xing, and K. Orth, Sci Signal 10:eaa14501, 2017, https://doi.org/10.1126/scisignal.aal4501). Herein, we show that this prosurvival response is caused by the action of the channel-forming effector VopQ that targets the host V-ATPase, resulting in lysosomal deacidification and inhibition of lysosome-autophagosome fusion. Recent structural studies have shown how VopQ interacts with the V-ATPase and, while in the ER, a V-ATPase assembly intermediate can interact with VopQ, causing a disruption in membrane integrity. Additionally, we observed that VopQ-mediated disruption of the V-ATPase activates the IRE1 branch of the unfolded protein response (UPR), resulting in an IRE1-dependent activation of ERK1/2 MAPK signaling. We also find that this early VopQ-dependent induction of ERK1/2 phosphorylation is terminated by the VopS-mediated inhibitory AMPylation of Rho GTPase signaling. Since VopS dampens VopQ-induced IRE1-dependent ERK1/2 activation, we propose that IRE1 activates ERK1/2 phosphorylation at or above the level of Rho GTPases. This study illustrates how temporally induced effectors can work as in tandem as agonist/antagonist to manipulate host signaling and reveals new connections between V-ATPase function, UPR, and MAPK signaling. IMPORTANCE Vibrio parahaemolyticus is a seafood-borne pathogen that encodes two type 3 secretion systems (T3SS). The first system, T3SS1, is thought to be maintained in all strains of V. parahaemolyticus to maintain survival in the environment, whereas the second system, T3SS2, is linked to clinical isolates and disease in humans. Here, we found that first system targets evolutionarily conserved signaling systems to manipulate host cells, eventually causing a rapid, orchestrated cells death within 3 h. We have found that the T3SS1 injects virulence factors that temporally manipulate host signaling. Within the first hour of infection, the effector VopQ acts first by activating host survival signals while diminishing the host cell apoptotic machinery. Less than an hour later, another effector, VopS, reverses activation and inhibition of these signaling systems, ultimately leading to death of the host cell. This work provides example of how pathogens have evolved to manipulate the interplay between T3SS effectors to regulate host signaling pathways.

Viruses are the most abundant life form and infect practically all organisms. Consequently, these obligate parasites are a major cause of human suffering and economic loss. Rossmann-like fold is the most populated fold among α/β-folds in the Protein Data Bank and proteins containing Rossmann-like fold constitute 22% of all known proteins 3D structures. Thus, analysis of viral proteins containing Rossmann-like domains could provide an understanding of viral biology and evolution as well as could propose possible targets for antiviral therapy. We provide functional and evolutionary analysis of viral proteins containing a Rossmann-like fold found in the evolutionary classification of protein domains (ECOD) database developed in our lab. We identified 81 protein families of bacterial, archeal, and eukaryotic viruses in light of their evolution-based ECOD classification and Pfam taxonomy. We defined their functional significance using enzymatic EC number assignments as well as domain-level family annotations.

Protein domains exist by themselves or in combination with other domains to form complex multi-domain proteins. Defining domain boundaries in proteins is essential for understanding their evolution and function but is not trivial. More specifically, partitioning domains that interact by forming a single β-sheet is known to be particularly troublesome for automatic structure-based domain decomposition pipelines. Here, we study edge-to-edge β-strand interactions between domains in a protein chain, to help define the boundaries for some more difficult cases where a single β-sheet spanning over two domains gives an appearance of one. We give a number of examples where β-strands belonging to a single β-sheet do not belong to a single domain and highlight the difficulties of automatic domain parsers on these examples. This work can be used as a baseline for defining domain boundaries in homologous proteins or proteins with similar domain interactions in the future.

Cholesterol homeostasis is mediated by Scap, a polytopic endoplasmic reticulum (ER) protein that transports sterol regulatory element-binding proteins from the ER to Golgi, where they are processed to forms that activate cholesterol synthesis. Scap has eight transmembrane helices and two large luminal loops, designated Loop1 and Loop7. We earlier provided indirect evidence that Loop1 binds to Loop7, allowing Scap to bind COPII proteins for transport in coated vesicles. When ER cholesterol rises, it binds to Loop1. We hypothesized that this causes dissociation from Loop7, abrogating COPII binding. Here we demonstrate direct binding of the two loops when expressed as isolated fragments or as a fusion protein. Expressed alone, Loop1 remained intracellular and membrane-bound. When Loop7 was co-expressed, it bound to Loop1, and the soluble complex was secreted. A Loop1-Loop7 fusion protein was also secreted, and the two loops remained bound when the linker between them was cleaved by a protease. Point mutations that disrupt the Loop1-Loop7 interaction prevented secretion of the Loop1-Loop7 fusion protein. These data provide direct documentation of intramolecular Loop1-Loop7 binding, a central event in cholesterol homeostasis.

We present an overview of contact-assisted predictions in the eleventh round of critical assessment of protein structure prediction (CASP11), which included four categories: predicted contacts (Tp), correct contacts (Tc), simulated sparse NMR contacts (Ts), and cross-linking contacts (Tx). Comparison of assisted to unassisted model quality highlighted a relatively poor overall performance in CASP11 using predicted Tp and crosslinked Tx contact information. However, average model quality significantly improved in the correct Tc and simulated NMR Ts categories for most targets, where maximum improvement of unassisted models reached an impressive 70 GDT_TS. Comparison of the performance in the correct Tc category to CASP10 suggested the improvement in CASP11 model quality originated from an increased number of provided contacts per target. Group rankings based on a combination of scores used in the CASP11 free modeling (FM) assessment for each category highlight four top-performing groups, with three from the Lee lab and one from the Baker lab. We used the overall performance of these groups in each category to develop hypotheses for their relative outperformance in the correct Tc and simulated NMR Ts categories, which stemmed from the fraction of correct contacts provided (correct Tc category) and a reduced fraction of correct contacts offset by an increased coverage of the correct contacts (simulated NMR Ts category). Proteins 2016; 84(Suppl 1):164-180. © 2016 Wiley Periodicals, Inc.

Proteins and their domains evolve by a set of events commonly including the duplication and divergence of small motifs. The presence of short repetitive regions in domains has generally constituted a difficult case for structural domain classifications and their hierarchies. We developed the Evolutionary Classification Of protein Domains (ECOD) in part to implement a new schema for the classification of these types of proteins. Here we document the ways in which ECOD classifies proteins with small internal repeats, widespread functional motifs, and assemblies of small domain-like fragments in its evolutionary schema. We illustrate the ways in which the structural genomics project impacted the classification and characterization of new structural domains and sequence families over the decade.

The human genome harbors a variety of genetic variations. Single-nucleotide changes that alter amino acids in protein-coding regions are one of the major causes of human phenotypic variation and diseases. These single-amino acid variations (SAVs) are routinely found in whole genome and exome sequencing. Evaluating the functional impact of such genomic alterations is crucial for diagnosis of genetic disorders. We developed DeepSAV, a deep-learning convolutional neural network to differentiate disease-causing and benign SAVs based on a variety of protein sequence, structural and functional properties. Our method outperforms most stand-alone programs, and the version incorporating population and gene-level information (DeepSAV+PG) has similar predictive power as some of the best available. We transformed DeepSAV scores of rare SAVs in the human population into a quantity termed “mutation severity measure” for each human protein-coding gene. It reflects a gene's tolerance to deleterious missense mutations and serves as a useful tool to study gene-disease associations. Genes implicated in cancer, autism, and viral interaction are found by this measure as intolerant to mutations, while genes associated with a number of other diseases are scored as tolerant. Among known disease-associated genes, those that are mutation-intolerant are likely to function in development and signal transduction pathways, while those that are mutation-tolerant tend to encode metabolic and mitochondrial proteins., Author summary Human genetic variations in various forms are constantly found in whole genome and exome sequencing of general population and patients. It remains a challenging task to assess the functional impact of these variations. In this study, we performed comprehensive analysis of single-amino-acid variations (SAVs) in terms of their sequence, structure, and functional properties. We further developed a deep neural network-based method to predict the functional impact of SAVs. Our method is among the top performers compared to existing programs in differentiating pathogenic and benign SAVs. We designed a mutation severity measure for human protein-coding genes by aggregating the predicted scores of SAVs found in the human general population. Such a measure reflects a gene's tolerance to deleterious missense mutations and serves as a useful tool to study gene-disease associations. We found that genes implicated in cancer, autism, and viral interaction are more likely to be intolerant to mutations than genes with other diseases. Disease-associated genes with strong mutation intolerance tend to function in development and signal transduction pathways. On the other end of the mutation severity spectrum, mutation-tolerant genes often encode proteins functioning in mitochondria and metabolic pathways.

Protein target structures for the Critical Assessment of Structure Prediction round 11 (CASP11) and CASP ROLL were split into domains and classified into categories suitable for assessment of template-based modeling (TBM) and free modeling (FM) based on their evolutionary relatedness to existing structures classified by the Evolutionary Classification of Protein Domains (ECOD) database. First, target structures were divided into domain-based evaluation units. Target splits were based on the domain organization of available templates as well as the performance of servers on whole targets compared to split target domains. Second, evaluation units were classified into TBM and FM categories using a combination of measures that evaluate prediction quality and template detectability. Generally, target domains with sequence-related templates and good server prediction performance were classified as TBM, whereas targets without sequence-identifiable templates and low server performance were classified as FM. As in previous CASP experiments, the boundaries for classification were blurred due to the presence of significant insertions and deteriorations in the targets with respect to homologous templates, as well as the presence of templates with partial coverage of new folds. The FM category included 45 target domains, which represents an unprecedented number of difficult CASP targets provided for modeling. Proteins 2016; 84(Suppl 1):20-33. © 2016 Wiley Periodicals, Inc.

We present an assessment of 'template-free modeling' (FM) in CASP11and ROLL. Community-wide server performance suggested the use of automated scores similar to previous CASPs would provide a good system of evaluating performance, even in the absence of comprehensive manual assessment. The CASP11 FM category included several outstanding examples, including successful prediction by the Baker group of a 256-residue target (T0806-D1) that lacked sequence similarity to any existing template. The top server model prediction by Zhang's Quark, which was apparently selected and refined by several manual groups, encompassed the entire fold of target T0837-D1. Methods from the same two groups tended to dominate overall CASP11 FM and ROLL rankings. Comparison of top FM predictions with those from the previous CASP experiment revealed progress in the category, particularly reflected in high prediction accuracy for larger protein domains. FM prediction models for two cases were sufficient to provide functional insights that were otherwise not obtainable by traditional sequence analysis methods. Importantly, CASP11 abstracts revealed that alignment-based contact prediction methods brought about much of the CASP11 progress, producing both of the functionally relevant models as well as several of the other outstanding structure predictions. These methodological advances enabled de novo modeling of much larger domain structures than was previously possible and allowed prediction of functional sites. Proteins 2016; 84(Suppl 1):51-66. © 2015 Wiley Periodicals, Inc.

The eukaryotic protozoan parasite Trypanosoma brucei is the causative agent of human African trypanosomiasis. Polyamine biosynthesis is essential in T. brucei, and the polyamine spermidine is required for synthesis of a novel cofactor called trypanothione and for deoxyhypusine modification of eukaryotic translation initiation factor 5A (eIF5A). eIF5A promotes translation of proteins containing polyprolyl tracts in mammals and yeast. To evaluate the function of eIF5A in T. brucei, we used RNA interference (RNAi) to knock down eIF5A levels and found that it is essential for T. brucei growth. The RNAi-induced growth defect was complemented by expression of wild-type human eIF5A but not by a Lys-50 mutant that blocks modification by deoxyhypusine. Bioinformatics analysis showed that 15% of the T. brucei proteome contains 3 or more consecutive prolines and that actin-related proteins and cysteine proteases were highly enriched in the group. Steady-state protein levels of representative proteins containing 9 consecutive prolines that are involved in actin assembly (formin and CAP/Srv2p) were significantly reduced by knockdown of eIF5A. Several T. brucei polyprolyl proteins are involved in flagellar assembly. Knockdown of TbeIF5A led to abnormal cell morphologies and detached flagella, suggesting that eIF5A is important for translation of proteins needed for these processes. Potential specialized functions for eIF5A in T. brucei in translation of variable surface glycoproteins were also uncovered. Inhibitors of deoxyhypusination would be expected to cause a pleomorphic effect on multiple cell processes, suggesting that deoxyhypusine/hypusine biosynthesis could be a promising drug target in not just T. brucei but in other eukaryotic pathogens.

Chameleon sequences (ChSeqs) refer to sequence strings of identical amino acids that can adopt different conformations in protein structures. Researchers have detected and studied ChSeqs to understand the interplay between local and global interactions in protein structure formation. The different secondary structures adopted by one ChSeq challenge sequence-based secondary structure predictors. With increasing numbers of available Protein Data Bank structures, we here identify a large set of ChSeqs ranging from 6 to 10 residues in length. The homologous ChSeqs discovered highlight the structural plasticity involved in biological function. When compared with previous studies, the set of unrelated ChSeqs found represents an about 20-fold increase in the number of detected sequences, as well as an increase in the longest ChSeq length from 8 to 10 residues. We applied secondary structure predictors on our ChSeqs and found that methods based on a sequence profile outperformed methods based on a single sequence. For the unrelated ChSeqs, the evolutionary information provided by the sequence profile typically allows successful prediction of the prevailing secondary structure adopted in each protein family. Our dataset will facilitate future studies of ChSeqs, as well as interpretations of the interplay between local and nonlocal interactions. A user-friendly web interface for this ChSeq database is available at prodata.swmed.edu/chseq.

SummaryThe existence of extracellular phosphoproteins has been acknowledged for over a century. However, research in this area has been undeveloped largely because the kinases that phosphorylate secreted proteins have escaped identification. Fam20C is a kinase that phosphorylates S-x-E/pS motifs on proteins in milk and in the extracellular matrix of bones and teeth. Here, we show that Fam20C generates the majority of the extracellular phosphoproteome. Using CRISPR/Cas9 genome editing, mass spectrometry, and biochemistry, we identify more than 100 secreted phosphoproteins as genuine Fam20C substrates. Further, we show that Fam20C exhibits broader substrate specificity than previously appreciated. Functional annotations of Fam20C substrates suggest roles for the kinase beyond biomineralization, including lipid homeostasis, wound healing, and cell migration and adhesion. Our results establish Fam20C as the major secretory pathway protein kinase and serve as a foundation for new areas of investigation into the role of secreted protein phosphorylation in human biology and disease.

Summary RUVBL1 and RUVBL2 (collectively RUVBL1/2) are essential AAA+ ATPases that function as co-chaperones and have been implicated in cancer. Here we investigated the molecular and phenotypic role of RUVBL1/2 ATPase activity in non-small cell lung cancer (NSCLC). We find that RUVBL1/2 are overexpressed in NSCLC patient tumors, with high expression associated with poor survival. Utilizing a specific inhibitor of RUVBL1/2 ATPase activity, we show that RUVBL1/2 ATPase activity is necessary for the maturation or dissociation of the PAQosome, a large RUVBL1/2-dependent multiprotein complex. We also show that RUVBL1/2 have roles in DNA replication, as inhibition of its ATPase activity can cause S-phase arrest, which culminates in cancer cell death via replication catastrophe. While in vivo pharmacological inhibition of RUVBL1/2 results in modest antitumor activity, it synergizes with radiation in NSCLC, but not normal cells, an attractive property for future preclinical development.

The heme-copper oxidase (HCO) superfamily includes HCOs in aerobic respiratory chains and nitric oxide reductases (NORs) in the denitrification pathway. The HCO/NOR catalytic subunit has a core structure consisting of 12 transmembrane helices (TMHs) arranged in three-fold rotational pseudosymmetry, with six conserved histidines for heme and metal binding. Using sensitive sequence similarity searches, we detected a number of novel HCO/NOR homologs and named them HCO Homology (HCOH) proteins. Several HCOH families possess only four TMHs that exhibit the most pronounced similarity to the last four TMHs (TMHs 9–12) of HCOs/NORs. Encoded by independent genes, four-TMH HCOH proteins represent a single evolutionary unit (EU) that relates to each of the three homologous EUs of HCOs/NORs comprising TMHs 1–4, TMHs 5–8, and TMHs 9–12. Single-EU HCOH proteins could form homotrimers or heterotrimers to maintain the general structure and ligand-binding sites defined by the HCO/NOR catalytic subunit fold. The remaining HCOH families, including NnrS, have 12-TMHs and three EUs. Most three-EU HCOH proteins possess two conserved histidines and could bind a single heme. Limited experimental studies and genomic context analysis suggest that many HCOH proteins could function in the denitrification pathway and in detoxification of reactive molecules such as nitric oxide. HCO/NOR catalytic subunits exhibit remarkable structural similarity to the homotrimers of MAPEG (membrane-associated proteins in eicosanoid and glutathione metabolism) proteins. Gene duplication, fusion, and fission likely play important roles in the evolution of HCOs/NORs and HCOH proteins.

Rossmann folds are ancient, frequently diverged domains found in many biological reaction pathways where they have adapted for different functions. Consequently, discernment and classification of their homologous relations and function can be complicated. We define a minimal Rossmann-like structure motif (RLM) that corresponds for the common core of known Rossmann domains and use this motif to identify all RLM domains in the Protein Data Bank (PDB), thus finding they constitute about 20% of all known 3D structures. The Evolutionary Classification of protein structure Domains (ECOD) classifies RLM domains in a number of groups that lack evidence for homology (X-groups), which suggests that they could have evolved independently multiple times. Closely related, homologous RLM enzyme families can diverge to bind different ligands using similar binding sites and to catalyze different reactions. Conversely, non-homologous RLM domains can converge to catalyze the same reactions or to bind the same ligand with alternate binding modes. We discuss a special case of such convergent evolution that is relevant to the polypharmacology paradigm, wherein the same drug (methotrexate) binds to multiple non-homologous RLM drug targets with different topologies. Finally, assigning proteins with RLM domain to the Enzyme Commission classification suggest that RLM enzymes function mainly in metabolism (and comprise 38% of reference metabolic pathways) and are overrepresented in extant pathways that represent ancient biosynthetic routes such as nucleotide metabolism, energy metabolism, and metabolism of amino acids. In fact, RLM enzymes take part in five out of eight enzymatic reactions of the Wood-Ljungdahl metabolic pathway thought to be used by the last universal common ancestor (LUCA). The prevalence of RLM domains in this ancient metabolism might explain their wide distribution among enzymes., Author summary Protein-ligand interactions are crucial to understanding molecular-level mechanisms of cell functions. A ligand-binding site frequently acts as a protein’s functional center, where its structural scaffold orients molecules in the correct conformation for proper function such as catalysis. These binding or active sites are usually conserved within protein families and can be used to define homologous proteins. However, even homologs may have diverged to functions requiring different active sites. Here we study the structural features and evolution of active sites of proteins containing the Rossmann-like motif (RLM). RLM proteins are numerous in nature and they constitute about 20% of all known 3D structures. Taking into account the prevalence of RLMs in proteins, analysis of their ligands provides general insights into enzyme evolution as well as proposes possible domain-based strategies for drug design in light of polypharmacology. Our analysis shows that RLM enzymes function predominantly in metabolism–these proteins cover 38% of reference metabolic pathways and are involved in numerous metabolic diseases. Being one of the most ancient folds, RLM enzymes could be traced to ancient metabolic pathways used by LUCA, which are still in use in some bacteria today.

© 2016 Leija et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited., The human pathogenic parasite Trypanosoma brucei possess both de novo and salvage routes for the biosynthesis of pyrimidine nucleotides. Consequently, they do not require salvageable pyrimidines for growth. Thymidine kinase (TK) catalyzes the formation of dTMP and dUMP and is one of several salvage enzymes that appear redundant to the de novo pathway. Surprisingly, we show through analysis of TK conditional null and RNAi cells that TK is essential for growth and for infectivity in a mouse model, and that a catalytically active enzyme is required for its function. Unlike humans, T. brucei and all other kinetoplastids lack dCMP deaminase (DCTD), which provides an alternative route to dUMP formation. Ectopic expression of human DCTD resulted in full rescue of the RNAi growth phenotype and allowed for selection of viable TK null cells. Metabolite profiling by LC-MS/MS revealed a buildup of deoxypyrimidine nucleosides in TK depleted cells. Knockout of cytidine deaminase (CDA), which converts deoxycytidine to deoxyuridine led to thymidine/deoxyuridine auxotrophy. These unexpected results suggested that T. brucei encodes an unidentified 5'-nucleotidase that converts deoxypyrimidine nucleotides to their corresponding nucleosides, leading to their dead-end buildup in TK depleted cells at the expense of dTTP pools. Bioinformatics analysis identified several potential candidate genes that could encode 5'-nucleotidase activity including an HD-domain protein that we show catalyzes dephosphorylation of deoxyribonucleotide 5'-monophosphates. We conclude that TK is essential for synthesis of thymine nucleotides regardless of whether the nucleoside precursors originate from the de novo pathway or through salvage. Reliance on TK in the absence of DCTD may be a shared vulnerability among trypanosomatids and may provide a unique opportunity to selectively target a diverse group of pathogenic single-celled eukaryotes with a single drug., This work was supported by National Institutes of Health (grants AI078962 and AI034432) to MAP (https://www.niaid.nih.gov) and (grant GM007062) to CL (https://www.nigms.nih. gov), the Welch Foundation (grant I-1257) to MAP and (grant I-1686) to JJK (http://www.welch1.org), and Fundac ̧ão para a Ciência e Tecnologia (FCT, Portugal) SFRH/BD/51286/2010 (http://www.fct.pt) to FRF.

The ‘pseudokinase’ SgK196 is a protein O-mannose kinase (POMK) that catalyzes an essential phosphorylation step during biosynthesis of the laminin-binding glycan on α-dystroglycan. However, the catalytic mechanism underlying this activity remains elusive. Here we present the crystal structure of Danio rerio POMK in complex with Mg2+ ions, ADP, aluminum fluoride, and the GalNAc-β3-GlcNAc-β4-Man trisaccharide substrate, thereby providing a snapshot of the catalytic transition state of this unusual kinase. The active site of POMK is established by residues located in non-canonical positions and is stabilized by a disulfide bridge. GalNAc-β3-GlcNAc-β4-Man is recognized by a surface groove, and the GalNAc-β3-GlcNAc moiety mediates the majority of interactions with POMK. Expression of various POMK mutants in POMK knockout cells further validated the functional requirements of critical residues. Our results provide important insights into the ability of POMK to function specifically as a glycan kinase, and highlight the structural diversity of the human kinome. DOI: http://dx.doi.org/10.7554/eLife.22238.001

Catalytically inactive enzyme paralogs occur in many genomes. Some regulate their active counterparts but the structural principles of this regulation remain largely unknown. We report X-ray structures of Trypanosoma brucei S-adenosylmethionine decarboxylase alone and in functional complex with its catalytically dead paralogous partner, prozyme. We show monomeric TbAdoMetDC is inactive because of autoinhibition by its N-terminal sequence. Heterodimerization with prozyme displaces this sequence from the active site through a complex mechanism involving a cis-to-trans proline isomerization, reorganization of a β-sheet, and insertion of the N-terminal α-helix into the heterodimer interface, leading to enzyme activation. We propose that the evolution of this intricate regulatory mechanism was facilitated by the acquisition of the dimerization domain, a single step that can in principle account for the divergence of regulatory schemes in the AdoMetDC enzyme family. These studies elucidate an allosteric mechanism in an enzyme and a plausible scheme by which such complex cooperativity evolved. DOI: http://dx.doi.org/10.7554/eLife.20198.001

ATP binding cassette (ABC) transporters play critical roles in maintaining sterol balance in higher eukaryotes. The ABCG5/ABCG8 heterodimer (G5G8) mediates excretion of neutral sterols in liver and intestines1–5. Mutations disrupting G5G8 cause sitosterolaemia, a disorder characterized by sterol accumulation and premature atherosclerosis. Here we use crystallization in lipid bilayers to determine the X-ray structure of human G5G8 in a nucleotide-free state at 3.9 Å resolution, generating the first atomic model of an ABC sterol transporter. The structure reveals a new transmembrane fold that is present in a large and functionally diverse superfamily of ABC transporters. The transmembrane domains are coupled to the nucleotide-binding sites by networks of interactions that differ between the active and inactive ATPases, reflecting the catalytic asymmetry of the transporter. The G5G8 structure provides a mechanistic framework for understanding sterol transport and the disruptive effects of mutations causing sitosterolaemia.

Summary Structural backbones of iron-scavenging siderophore molecules include polyamines 1,3-diaminopropane and 1,5-diaminopentane (cadaverine). For the cadaverine-based desferroxiamine E siderophore in Streptomyces coelicolor, the corresponding biosynthetic gene cluster contains an ORF encoded by desA that was suspected of producing the cadaverine (decarboxylated lysine) backbone. However, desA encodes an l-2,4-diaminobutyrate decarboxylase (DABA DC) homologue and not any known form of lysine decarboxylase (LDC). The only known function of DABA DC is, together with l-2,4-aminobutyrate aminotransferase (DABA AT), to synthesize 1,3-diaminopropane. We show here that S. coelicolor desA encodes a novel LDC and we hypothesized that DABA DC homologues present in siderophore biosynthetic clusters in the absence of DABA AT ORFs would be novel LDCs. We confirmed this by correctly predicting the LDC activity of a DABA DC homologue from a Yersinia pestis siderophore biosynthetic pathway. The corollary was confirmed for a DABA DC homologue, adjacent to a DABA AT ORF in a siderophore pathway in the cyanobacterium Anabaena variabilis, which was shown to be a bona fide DABA DC. These findings enable prediction of whether a siderophore pathway will utilize 1,3-diaminopropane or cadaverine, and suggest that the majority of bacteria use DABA AT and DABA DC for siderophore, rather than norspermidine/polyamine biosynthesis.

Candidatus Liberibacter asiaticus (Ca. L. asiaticus) is a Gram-negative bacterium and the pathogen of Citrus Greening disease (Huanglongbing, HLB). As a parasitic bacterium, Ca. L. asiaticus harbors ABC transporters that play important roles in exchanging chemical compounds between Ca. L. asiaticus and its host. Here, we analyzed all the ABC transporter-related proteins in Ca. L. asiaticus. We identified 14 ABC transporter systems and predicted their structures and substrate specificities. In-depth sequence and structure analysis including multiple sequence alignment, phylogenetic tree reconstruction, and structure comparison further support their function predictions. Our study shows that this bacterium could use these ABC transporters to import metabolites (amino acids and phosphates) and enzyme cofactors (choline, thiamine, iron, manganese, and zinc), resist to organic solvent, heavy metal, and lipid-like drugs, maintain the composition of the outer membrane (OM), and secrete virulence factors. Although the features of most ABC systems could be deduced from the abundant experimental data on their orthologs, we reported several novel observations within ABC system proteins. Moreover, we identified seven nontransport ABC systems that are likely involved in virulence gene expression regulation, transposon excision regulation, and DNA repair. Our analysis reveals several candidates for further studies to understand and control the disease, including the type I virulence factor secretion system and its substrate that are likely related to Ca. L. asiaticus pathogenicity and the ABC transporter systems responsible for bacterial OM biosynthesis that are good drug targets.

We present an overview of the ninth round of Critical Assessment of Protein Structure Prediction (CASP9) "Template free modeling" category (FM). Prediction models were evaluated using a combination of established structural and sequence comparison measures and a novel automated method designed to mimic manual inspection by capturing both global and local structural features. These scores were compared to those assigned manually over a diverse subset of target domains. Scores were combined to compare overall performance of participating groups and to estimate rank significance. Moreover, we discuss a few examples of free modeling targets to highlight the progress and bottlenecks of current prediction methods. Notably, a server prediction model for a single target (T0581) improved significantly over the closest structure template (44% GDT increase). This accomplishment represents the "winner" of the CASP9 FM category. A number of human expert groups submitted slight variations of this model, highlighting a trend for human experts to act as "meta predictors" by correctly selecting among models produced by the top-performing automated servers. The details of evaluation are available at http://prodata.swmed.edu/CASP9/ .

The Critical assessment of protein structure prediction round 9 (CASP9) aimed to evaluate predictions for 129 experimentally determined protein structures. To assess tertiary structure predictions, these target structures were divided into domain-based evaluation units that were then classified into two assessment categories: template based modeling (TBM) and template free modeling (FM). CASP9 targets were split into domains of structurally compact evolutionary modules. For the targets with more than one defined domain, the decision to split structures into domains for evaluation was based on server performance. Target domains were categorized based on their evolutionary relatedness to existing templates as well as their difficulty levels indicated by server performance. Those target domains with sequence-related templates and high server prediction performance were classified as TMB, whereas those targets without identifiable templates and low server performance were classified as FM. However, using these generalizations for classification resulted in a blurred boundary between CASP9 assessment categories. Thus, the FM category included those domains without sequence detectable templates (25 target domains) as well as some domains with difficult to detect templates whose predictions were as poor as those without templates (five target domains). Several interesting examples are discussed, including targets with sequence related templates that exhibit unusual structural differences, targets with homologous or analogous structure templates that are not detectable by sequence, and targets with new folds. Proteins 2011; © 2011 Wiley-Liss, Inc.

Approximately 10% of human protein kinases are believed to be inactive and named pseudokinases because they lack residues required for catalysis. Here, we show that the highly conserved pseudokinase selenoprotein-O (SelO) transfers AMP from ATP to Ser, Thr, and Tyr residues on protein substrates (AMPylation), uncovering a previously unrecognized activity for a member of the protein kinase superfamily. The crystal structure of a SelO homolog reveals a protein kinase-like fold with ATP flipped in the active site, thus providing a structural basis for catalysis. SelO pseudokinases localize to the mitochondria and AMPylate proteins involved in redox homeostasis. Consequently, SelO activity is necessary for the proper cellular response to oxidative stress. Our results suggest that AMPylation may be a more widespread post-translational modification than previously appreciated and that pseudokinases should be analyzed for alternative transferase activities.

REDD1 is a conserved stress-response protein that regulates mTORCl, a critical regulator of cell growth and proliferation that is implicated in cancer. REDD1 is induced by hypoxia, and REDD1 overexpression is sufficient to inhibit mTORC1. mTORC1 is regulated by the small GTPase Rheb, which in turn is regulated by the GTPase-activating protein complex, TSC1/TSC2. REDD1 induced-mTORC1 inhibition requires the TSC1/TSC2 complex, and REDD1 has been proposed to act by directly binding to and sequestering 14-3-3 proteins away from TSC2 leading to TSC2-depedent inhibition of mTORC1. Structure/function analyses have led us to identify two segments in REDD1 that are essential for function, which act in an interdependent manner. We have determined a crystal structure of REDD1 at 2.0 Å resolution, which shows that these two segments fold together to form an intact domain with a novel fold. This domain is characterized by an α/β sandwich consisting of two antiparallel α-helices and a mixed. β-sheet encompassing an uncommon psi-loop motif. Structure-based docking and functional analyses suggest that REDD1 does not directly bind to 14-3-3 proteins. Sequence conservation mapping to the surface of the structure and mutagenesis studies demarcated a hotspot likely to interact with effector proteins that is essential for REDD1-mediated mTORC1 inhibition. © 2010 American Chemical Society.

Obesity and insulin resistance are associated with deposition of triglycerides in tissues other than adipose tissue. Previously, we showed that a missense mutation (I148M) in PNPLA3 (patatin-like phospholipase domain-containing 3 protein) is associated with increased hepatic triglyceride content in humans. Here we examined the effect of the I148M substitution on the enzymatic activity and cellular location of PNPLA3. Structural modeling predicted that the substitution of methionine for isoleucine at residue 148 would restrict access of substrate to the catalytic serine at residue 47. In vitro assays using recombinant PNPLA3 partially purified from Sf9 cells confirmed that the wild type enzyme hydrolyzes emulsified triglyceride and that the I148M substitution abolishes this activity. Expression of PNPLA3-I148M, but not wild type PNPLA3, in cultured hepatocytes or in the livers of mice increased cellular triglyceride content. Cell fractionation studies revealed that approximately 90% of wild type PNPLA3 partitioned between membranes and lipid droplets; substitution of isoleucine for methionine at position 148 did not alter the subcellular distribution of the protein. These data are consistent with PNPLA3-I148M promoting triglyceride accumulation by limiting triglyceride hydrolysis.

Iron Sensor Intracellular iron is an essential cofactor for many proteins, but can also damage macromolecules, so its levels are carefully controlled. Cellular iron homeostasis is mediated by iron regulatory proteins that regulate the expression of genes involved in iron uptake and storage. However, it is not clear how cells sense iron bioavailability (see the Perspective by Rouault ). Using different approaches, Salahudeen et al. (p. 722 , published online 17 September) and Vashisht et al. (p. 718 , published online 17 September) have identified the F-box protein FBXL5 as a human iron sensor. FBXL5 is part of an E3 ubiquitin ligase complex that regulates the degradation of iron regulatory proteins and thereby cellular iron levels. It contains a hemerythrin domain that binds iron and acts as an iron-dependent regulatory switch, causing the degradation of FBXL5 under low iron conditions. This alternative pathway for the regulation of iron homeostasis has implications for both normal cellular physiology and disease.

JAK2 (Janus kinase 2) is essential for cytokine receptor signaling, and several lines of evidence support a causal role of an activating JAK2 mutation in myeloproliferative disorders. JAK2 activity is autoinhibited by its pseudokinase domain in the basal state, and the inhibition is released by cytokine stimulation; how engagement of the cognate receptor triggers this release is unknown. From a functional screen for gain-of-function JAK2 mutations, we discovered 13 missense mutations, nine in the pseudokinase domain and four in the Src homology 2 (SH2)-pseudokinase domain linker. These mutations identified determinants for autoinhibition and inducible activation in JAK2. Two of the mutants, K539I and N622I, resulted in erythrocytosis in mice. Scanning mutagenesis of the SH2-pseudokinase domain linker indicated that its N-terminal part was essential for interaction of JAK2 with the Epo receptor, whereas certain mutations in the C-terminal region conferred constitutive activation. We further showed that substitutions for Glu(543)-Asp(544) in this linker or Leu(611), Arg(683), or Phe(694) in the hinge proximal region of the pseudokinase domain resulted in activated JAK2 mutants that could not be further stimulated by Epo. These results suggest that the SH2-pseudokinase domain linker acts as a switch that relays cytokine engagement to JAK2 activation by flexing the pseudokinase domain hinge.

The Vibrio parahaemolyticus type III effector VopS is implicated in cell rounding and the collapse of the actin cytoskeleton by inhibiting Rho guanosine triphosphatases (GTPases). We found that VopS could act to covalently modify a conserved threonine residue on Rho, Rac, and Cdc42 with adenosine 5′-monophosphate (AMP). The resulting AMPylation prevented the interaction of Rho GTPases with downstream effectors, thereby inhibiting actin assembly in the infected cell. Eukaryotic proteins were also directly modified with AMP, potentially expanding the repertoire of posttranslational modifications for molecular signaling.

Defects in Niemann-Pick, Type C-1 protein (NPC1) cause cholesterol, sphingolipids, phospholipids, and glycolipids to accumulate in lysosomes of liver, spleen, and brain. In cultured fibroblasts, NPC1 deficiency causes lysosomal retention of lipoprotein-derived cholesterol after uptake by receptor-mediated endocytosis. NPC1 contains 1278 amino acids that form 13 membrane-spanning helices and three large loops that project into the lumen of lysosomes. We showed earlier that NPC1 binds cholesterol and oxysterols. Here we localize the binding site to luminal loop-1, a 240-amino acid domain with 18 cysteines. When produced in cultured cells, luminal loop-1 was secreted as a soluble dimer. This loop bound [(3)H]cholesterol (K(d), 130 nM) and [(3)H]25-hydroxycholesterol (25-HC, K(d), 10 nM) with one sterol binding site per dimer. Binding of both sterols was competed by oxysterols (24-, 25-, and 27-HC). Unlabeled cholesterol competed strongly for binding of [(3)H]cholesterol, but weakly for [(3)H]25-HC binding. Binding of [(3)H]cholesterol but not [(3)H]25-HC was inhibited by detergents. We also studied NPC2, a soluble protein whose deficiency causes a similar disease phenotype. NPC2 bound cholesterol, but not oxysterols. Epicholesterol and cholesteryl sulfate competed for [(3)H]cholesterol binding to NPC2, but not NPC1. Glutamine 79 in luminal loop-1 of NPC-1 is important for sterol binding; a Q79A mutation abolished binding of [(3)H]cholesterol and [(3)H]25-HC to full-length NPC1. Nevertheless, the Q79A mutant restored cholesterol transport to NPC1-deficient Chinese hamster ovary cells. Thus, the sterol binding site on luminal loop-1 is not essential for NPC1 function in fibroblasts, but it may function in other cells where NPC1 deficiency produces more complicated lipid abnormalities.

X-ray crystallographic protein structures often contain disordered regions that are observed as missing electron density. Diffraction data may give little or no direct evidence as to the specific nature of disordered regions. We have developed a weighted window-based disorder predictor optimized using crystallographic data. Performance of a predictor is strongly influenced by chain termini. Optimized score adjustment values for amino- and carboxy-terminal positions demonstrate a simple, monotonic relationship between disorder and residue distance from termini. This optimized disorder predictor performs similarly to DISOPRED2 on crystallographically disordered regions. Data-optimized residue disorder propensities show strong linear correlation with experimentally determined amino acid transfer energies between water and hydrogen-bonding organic solvents, which primarily reflect residue hydrophobicity (exemplified by the Nozaki-Tanford hydrophobicity scale). Disorder propensities do not correlate as well with transfer energies between water and apolar solvents, which primarily reflect a different hydropathic property: residue hydrophilicity (also reflected by the Kyte-Doolittle hydropathy scale). Our results suggest that while hydrophobic side-chain interactions are primarily involved in determining stability of the folded conformation, hydrogen bonding, and similar polar interactions are primarily involved in conformational and interaction specificity.

Protein phosphorylation is a nearly universal post-translation modification involved in a plethora of cellular events. Even though phosphorylation of extracellular proteins had been observed, the identity of the kinases that phosphorylate secreted proteins remained a mystery until only recently. Advances in genome sequencing and genetic studies have paved the way for the discovery of a new class of kinases that localize within the endoplasmic reticulum, Golgi apparatus and the extracellular space. These novel kinases phosphorylate proteins and proteoglycans in the secretory pathway and appear to regulate various extracellular processes. Mutations in these kinases cause human disease, thus underscoring the biological importance of phosphorylation within the secretory pathway. This article is part of a Special Issue entitled: Inhibitors of Protein Kinases.

Protein phosphorylation is a fundamental mechanism regulating nearly every aspect of life. Several secreted proteins are phosphorylated, but the kinases responsible are unknown. We identified a family of atypical protein kinases which have N-terminal signal sequence and localize to the lumen of the endoplasmic reticulum and the Golgi apparatus. One of these kinases known as Fam20C is the Golgi casein kinase and it phosphorylates secretory pathway proteins within S-x-E motifs. Fam20C phosphorylates the caseins and many other secretory proteins including the small integrin-binding ligand, N-linked glycoproteins (SIBLINGs) implicated in biomineralization. Consequently, mutations in Fam20C cause an osteosclerotic bone dysplasia in humans known as Raine syndrome. Another family member, Fam20B, is a xylose kinase important in the maturation of proteoglycans. The function of other FAM 20C family members will also be addressed.

Using top-of-the-range fold-recognition methods, we have assigned a thioredoxin-like structure to a family of previously uncharacterized hypothetical proteins of bacterial origin. The DCC family, named after the conserved N-terminal DxxCxxC motif, encompasses proteins of unknown function from DUF393 (in Pfam database) and COG3011. The presence of two invariant potentially catalytic cysteine residues indicates that DCC proteins function as thiol-disulfide oxidoreductases.

We present an overview of the fifth round of Critical Assessment of Protein Structure Prediction (CASP5) fold recognition category. Prediction models were evaluated by using six different structural measures and four different alignment measures, and these scores were compared to those assigned manually over a diverse subset of target domains. Scores were combined to compare overall performance of participating groups and to estimate rank significance. The methods used by a few groups outperformed all other methods in terms of the evaluated criteria and could be considered state-of-the-art in structure prediction. We discuss a few examples of difficult fold recognition targets to highlight the progress of ab initio-type methods on difficult structure analogs and the difficulties of predicting multidomain targets and selecting prediction models. We also compared the results of manual groups to those of automatic servers evaluated in parallel by CAFASP, showing that the top performing automated server structure predictions approached those of the best manual predictors.

The carboxyl-terminal domain of phospholipase C-beta is required for its stimulation by Galpha(q) and for its Galpha(q)-specific GTPase-activating protein (GAP) activity. We subjected this domain to a combination of deletion and alanine/glycine scanning mutagenesis to detect mutations that would inhibit either responsiveness to G(q) or G(q) GAP activity. Most mutations that altered either response or GAP activity diminished both in parallel. Many of these mutations map at the interface at which the carboxyl-terminal domain was recently shown to form a dimer (Singer, A. U., et al. (2001) Nat. Struct. Biol., 9, 32-36). Most others clustered in an area that is a plausible Galpha(q) binding site. In addition, one mutation that differentially inhibited GAP activity relative to responsiveness to Galpha(q) mapped in this region at a location modeled to be in close contact with the switch II region of Galpha(q). This is the site at which RGS proteins are thought to exert their GAP activity. Last, a deletion mutation differentially inhibited the response of phospholipase C-beta1 to Galpha(q) without blocking GAP activity. Its location in the molecule suggests that moving the attachment point of the catalytic domain can disrupt its ability to be activated by Galpha(q).

Understanding the evolution of a protein, including both close and distant relationships, often reveals insight into its structure and function. Fast and easy access to such up-to-date information facilitates research. We have developed a hierarchical evolutionary classification of all proteins with experimentally determined spatial structures, and presented it as an interactive and updatable online database. ECOD (Evolutionary Classification of protein Domains) is distinct from other structural classifications in that it groups domains primarily by evolutionary relationships (homology), rather than topology (or “fold”). This distinction highlights cases of homology between domains of differing topology to aid in understanding of protein structure evolution. ECOD uniquely emphasizes distantly related homologs that are difficult to detect, and thus catalogs the largest number of evolutionary links among structural domain classifications. Placing distant homologs together underscores the ancestral similarities of these proteins and draws attention to the most important regions of sequence and structure, as well as conserved functional sites. ECOD also recognizes closer sequence-based relationships between protein domains. Currently, approximately 100,000 protein structures are classified in ECOD into 9,000 sequence families clustered into close to 2,000 evolutionary groups. The classification is assisted by an automated pipeline that quickly and consistently classifies weekly releases of PDB structures and allows for continual updates. This synchronization with PDB uniquely distinguishes ECOD among all protein classifications. Finally, we present several case studies of homologous proteins not recorded in other classifications, illustrating the potential of how ECOD can be used to further biological and evolutionary studies., Author Summary Protein structural domain databases offer a vital resource for structural bioinformatics. These databases provide functional inference for homologous structures, supply templates for structural prediction experiments, and differentiate between homologs and analogs. The rate of structure determination and deposition has increased dramatically over recent years, overwhelming the ability of current classifications to incorporate all new structures. We have developed a fast and reliable methodology for updating domain databases automatically, and created a revised hierarchy for domain classification that emphasizes evolutionary relationships. By classifying all known structures in our database with continuing automatic updates, we provide an up-to-date alternative to current resources. We illustrate several concepts that guided our classification scheme with examples of homology between domains in ECOD that are not observed in other resources.

Trypanosoma cruzi S-adenosylmethionine decarboxylase (AdoMetDC) catalyzes the pyruvoyl-dependent decarboxylation of S-adenosylmethionine (AdoMet), which is an important step in the biosynthesis of polyamines. The time course of the AdoMetDC reaction under single-turnover conditions was measured to determine the rate of the slowest catalytic step up to and including decarboxylation. Analysis of this single-turnover data yields an apparent second-order rate constant for this reaction of 3300 M(-1) s(-1) in the presence of putrescine, which corresponds to a catalytic rate of6 s(-1). This rate is minimally 100-fold faster than the steady-state rate suggesting that product release, which includes Schiff base hydrolysis, limits the overall reaction. AdoMetDC exhibits an inverse solvent isotope effect on the single-turnover kinetics, and the pH profile predicts a pK(a) of 8.9 for the basic limb. These results are consistent with a Cys residue functioning as a general acid in the rate-determining step of the single-turnover reaction. Mutation of Cys-82 to Ala reduces the rate of the single turnover reaction to 11 M(-1) s(-1) in the presence of putrescine. Further, a solvent isotope effect is not observed for the mutant enzyme. Reduction of the wild-type enzyme with cyanoborohydride traps the Schiff base between the enzyme and decarboxylated substrate, while little Schiff base species of either substrate or product was trapped with the C82A mutant. These data suggest that Cys-82 functions as a general acid/base to catalyze Schiff base formation and hydrolysis. The solvent isotope and pH effects are mirrored in single-turnover analysis of reactions without the putrescine activator, yielding an apparent second-order rate constant of 150 M(-1) s(-1). The presence of putrescine increases the single-turnover rate by 20-fold, while it has relatively little effect on the affinity of the enzyme for product. Therefore, putrescine likely activates the T. cruzi AdoMetDC enzyme by accelerating the rate of Schiff base exchange.

The gene for S-adenosylmethionine decarboxylase (AdoMetDC), a rate-limiting enzyme in the biosynthesis of polyamines, has been cloned from a Trypanosoma cruz cDNA library. The cDNA clone contains a 1.1 kb open reading frame predicted to encode a 42 kDa protein that shares 31% sequence identity to the human proenzyme. T. cruzi AdoMetDC expressed and purified from E. coli is auto-catalytically processed into two subunits of 32 kDa (alpha) and 10 kDa (beta). The catalytic activity of the purified recombinant enzyme is activated by the addition of putrescine to the reaction. To determine the effect of putrescine on the kinetics of the reaction, the velocity data collected at various substrate and putrescine concentrations were fit to the rate equation describing a non-essential activator. In the presence of fully saturating putrescine, k(cat) increases by 9-fold over the unactivated rate to 0.06 s(-1). The model derived Km for AdoMet is 0.05 mM in the absence of putrescine and the model-derived Kd for putrescine binding to free enzyme is 2.5 mM. The Km for AdoMet increases by alpha 2-fold when the enzyme is fully saturated with putrescine. Unlike human AdoMetDC, cadaverine activates the T. cruzi enzyme to a similar extent as putrescine.

Background The size of the protein sequence database has been exponentially increasing due to advances in genome sequencing. However, experimentally characterized proteins only constitute a small portion of the database, such that the majority of sequences have been annotated by computational approaches. Current automatic annotation pipelines inevitably introduce errors, making the annotations unreliable. Instead of such error-prone automatic annotations, functional interpretation should rely on annotations of ‘reference proteins’ that have been experimentally characterized or manually curated. Results The Seq2Ref server uses BLAST to detect proteins homologous to a query sequence and identifies the reference proteins among them. Seq2Ref then reports publications with experimental characterizations of the identified reference proteins that might be relevant to the query. Furthermore, a plurality-based rating system is developed to evaluate the homologous relationships and rank the reference proteins by their relevance to the query. Conclusions The reference proteins detected by our server will lend insight into proteins of unknown function and provide extensive information to develop in-depth understanding of uncharacterized proteins. Seq2Ref is available at: http://prodata.swmed.edu/seq2ref.

The covalent attachment of adenosine monophosphate (AMP) to proteins, a process called AMPylation (adenylylation), has recently emerged as a novel theme in microbial pathogenesis. Although several AMPylating enzymes have been characterized, the only known virulence protein with de-AMPylation activity is SidD from the human pathogen Legionella pneumophila. SidD de-AMPylates mammalian Rab1, a small GTPase involved in secretory vesicle transport, thereby targeting the host protein for inactivation. The molecular mechanisms underlying Rab1 recognition and de-AMPylation by SidD are unclear. Here, we report the crystal structure of the catalytic region of SidD at 1.6 Å resolution. The structure reveals a phosphatase-like fold with additional structural elements not present in generic PP2C-type phosphatases. The catalytic pocket contains a binuclear metal-binding site characteristic of hydrolytic metalloenzymes, with strong dependency on magnesium ions. Subsequent docking and molecular dynamics simulations between SidD and Rab1 revealed the interface contacts and the energetic contribution of key residues to the interaction. In conjunction with an extensive structure-based mutational analysis, we provide in vivo and in vitro evidence for a remarkable adaptation of SidD to its host cell target Rab1 which explains how this effector confers specificity to the reaction it catalyses., Author Summary The covalent attachment of adenosine monophosphate (AMP) to proteins, a process called AMPylation (adenylylation), has recently emerged as a novel theme in microbial pathogenesis. While AMPylases from various pathogenic microorganisms have recently been characterized, the only virulence protein with de-AMPylation activity known to date is the Legionella pneumophila effector SidD which catalyzes AMP removal from the host GTPase Rab1. Thus, both AMPylation and de-AMPylation constitute a novel catalytic mechanism to precisely control the function and membrane dynamics of a host Rab GTPase. In spite of this pivotal role, the molecular mechanism of AMP removal and the structural determinants for Rab1 recognition by SidD have remained largely unexplored. Here, we present the crystal structure of the de-adenylylation domain of SidD and reveal the catalytic mechanism of Rab1 de-adenylylation. Surprisingly, the structure of SidD is not related to the other known enzyme with de-AMPylation activity, the Escherichia coli GS-ATase. Instead, the catalitic domain of SidD is remarkably similar to that of the metal-dependent protein phosphatases (PPMs), however with distinctive structural features to distinguish AMPylated Rab1 from similarly modified substrates. Importantly, we provide a model for the SidD-Rab1 complex which sheds light into the specific details of substrate recognition and catalysis by this virulence factor.

Ornithine decarboxylases from Trypanosoma brucei, mouse, and Leishmania donovani share strict specificity for three basic amino acids, ornithine, lysine, and arginine. To identify residues involved in this substrate specificity and/or in the reaction chemistry, six conserved acidic resides (Asp-88, Glu-94, Asp-233, Glu-274, Asp-361, and Asp-364) were mutated to alanine in the T. brucei enzyme. Each mutation causes a substantial loss in enzyme efficiency. Most notably, mutation of Asp-361 increases the Km for ornithine by 2000-fold, with little effect on kcat, suggesting that this residue is an important substrate binding determinant. Mutation of the only strictly conserved acidic residue, Glu-274, decreases kcat 50-fold; however, substitution of N-methylpyridoxal-5'-phosphate for pyridoxal-5'-phosphate as the cofactor in the reaction restores the kcat of E274A to wild-type levels. These data demonstrate that Glu-274 interacts with the protonated pyridine nitrogen of the cofactor to enhance the electron withdrawing capability of the ring, analogous to Asp-222 in aspartate aminotransferase (Onuffer, J. J., and Kirsch, J. F. (1994) Protein Eng. 7, 413-424). Eukaryotic ornithine decarboxylase is a homodimer with two shared active sites. Residues 88, 94, 233, and 274 are contributed to each active site from the same subunit as Lys-69, while residues 361 and 364 are part of the Cys-360 subunit.