Evolutionarily conserved regions and hydrophobic contacts at the superfamily level: The case of the fold-type I, pyridoxal-5'-phosphate-dependent enzymes

Structure of the Oxygen, Pyridoxal Phosphate-Dependent Capuramycin Biosynthetic Protein Cap15

Pyridoxal phosphate-dependent enzymes able to use oxygen as a co-substrate have emerged in multiple protein families. Here, we use crystallography to solve the 2.40 Å resolution crystal structure of Cap15, a nucleoside biosynthetic enzyme that catalyzes the oxidative decarboxylation of glycyl uridine. Our structural study captures the internal aldimine, pinpointing the active site lysine as K230 and showing the site of phosphate binding. Our docking studies reveal how Cap15 is able to catalyze a stereoselective deprotonation reaction, and bioinformatic analysis reveals active site residues that distinguish Cap15 from the structurally related d-glucosaminate-6-phosphate ammonia lyase and l-seryl-tRNA(Sec) selenium transferase (SelA). Our work provides the structural basis for further mechanistic investigation of a unique biosynthetic enzyme and provides a blueprint for understanding how oxygen reactivity emerged in the SelA-like protein family.

Biochemical characterization of the first step in sulfonolipid biosynthesis in Alistipes finegoldii

Sulfonolipids are unusual lipids found in the outer membranes of Gram-negative bacteria in the phylum Bacteroidetes. Sulfonolipid and its deacylated derivative, capnine, are sulfur analogs of ceramide-1-phosphate and sphingosine-1-phosphate, respectively; thus, sulfonolipid biosynthesis is postulated to be similar to the sphingolipid biosynthetic pathway. Here, we identify the first enzyme in sulfonolipid synthesis in Alistipes finegoldii, an anaerobic gut commensal bacterium, as the product of the alfi_1224 gene, cysteate acyl-acyl carrier protein (ACP) transferase (SulA). We show SulA catalyzes the condensation of acyl-ACP and cysteate (3-sulfo-alanine) to form 3-ketocapnine. Acyl-CoA is a poor substrate. We show SulA has a bound pyridoxal phosphate (PLP) cofactor that undergoes a spectral redshift in the presence of cysteate, consistent with the transition of the lysine-aldimine complex to a substrate-aldimine complex. Furthermore, the SulA crystal structure shows the same prototypical fold found in bacterial serine palmitoyltransferases (Spt), enveloping the PLP cofactor bound to Lys251. We observed the SulA and Spt active sites are identical except for Lys281 in SulA, which is an alanine in Spt. Additionally, SulA(K281A) is catalytically inactive, but binds cysteate and forms the external aldimine normally, highlighting the structural role of the Lys281 side chain in walling off the active site from bulk solvent. Finally, the electropositive groove on the protein surface adjacent to the active site entrance provides a landing pad for the electronegative acyl-ACP surface. Taken together, these data identify the substrates, products, and mechanism of SulA, the PLP-dependent condensing enzyme that catalyzes the first step in sulfonolipid synthesis in a gut commensal bacterium.

Structure and Mechanism of d-Glucosaminate-6-phosphate Ammonia-lyase: A Novel Octameric Assembly for a Pyridoxal 5'-Phosphate-Dependent Enzyme, and Unprecedented Stereochemical Inversion in the Elimination Reaction of a d-Amino Acid

d-Glucosaminate-6-phosphate ammonia-lyase (DGL) is a pyridoxal 5'-phosphate (PLP)-dependent enzyme that produces 2-keto-3-deoxygluconate 6-phosphate (KDG-6-P) in the metabolism of d-glucosaminic acid by Salmonella enterica serovar typhimurium. We have determined the crystal structure of DGL by SAD phasing with selenomethionine to a resolution of 2.58 Å. The sequence has very low identity with most other members of the aminotransferase (AT) superfamily. The structure forms an octameric assembly as a tetramer of dimers that has not been observed previously in the AT superfamily. PLP is covalently bound as a Schiff base to Lys-213 in the catalytic dimer at the interface of two monomers. The structure lacks the conserved arginine that binds the α-carboxylate of the substrate in most members of the AT superfamily. However, there is a cluster of arginines in the small domain that likely serves as a binding site for the phosphate of the substrate. The deamination reaction performed in D2O gives a KDG-6-P product stereospecifically deuterated at C3; thus, the mechanism must involve an enamine intermediate that is protonated by the enzyme before product release. Nuclear magnetic resonance (NMR) analysis demonstrates that the deuterium is located in the pro-R position in the product, showing that the elimination of water takes place with inversion of configuration at C3, which is unprecedented for a PLP-dependent dehydratase/deaminase. On the basis of the crystal structure and the NMR data, a reaction mechanism for DGL is proposed.

Large scale analyses of genotype-phenotype relationships of glycine decarboxylase mutations and neurological disease severity

Monogenetic diseases provide unique opportunity for studying complex, clinical states that underlie neurological severity. Loss of glycine decarboxylase (GLDC) can severely impact neurological development as seen in non-ketotic hyperglycinemia (NKH). NKH is a neuro-metabolic disorder lacking quantitative predictors of disease states. It is characterized by elevation of glycine, seizures and failure to thrive, but glycine reduction often fails to confer neurological benefit, suggesting need for alternate tools to distinguish severe from attenuated disease. A major challenge has been that there are 255 unique disease-causing missense mutations in GLDC, of which 206 remain entirely uncharacterized. Here we report a Multiparametric Mutation Score (MMS) developed by combining in silico predictions of stability, evolutionary conservation and protein interaction models and suitable to assess 251 of 255 mutations. In addition, we created a quantitative scale of clinical disease severity comprising of four major disease domains (seizure, cognitive failure, muscular and motor control and brain-malformation) to comprehensively score patient symptoms identified in 131 clinical reports published over the last 15 years. The resulting patient Clinical Outcomes Scores (COS) were used to optimize the MMS for biological and clinical relevance and yield a patient Weighted Multiparametric Mutation Score (WMMS) that separates severe from attenuated neurological disease (p = 1.2 e-5). Our study provides understanding for developing quantitative tools to predict clinical severity of neurological disease and a clinical scale that advances monitoring disease progression needed to evaluate new treatments for NKH.

Neurodegenerative disorders frequently have diverse, severe symptoms and health outcomes that can be difficult to predict. The rare disease non-ketotic hyperglycinemia (NKH) additionally has a wide range of disease-causing mutations in glycine decarboxylase (GLDC), a protein that breaks down glycine. But measuring glycine is not sufficient to foretell disease outcome. A method to predict whether a mutation will cause severe or more mild forms of NKH would be very helpful to both understanding the disease as well as developing treatments for it. We used computation-based approaches to develop a mutation score that comprehensively predicts how mutations decrease GLDC function. After training against clinical data, the score was able to predict whether a mutation will cause severe or attenuated disease. This study utilizes the power of computational and multidisciplinary analyses to advance understanding and treatment of genetically caused neurodegenerative diseases.

Structural and kinetic properties of serine hydroxymethyltransferase from the halophytic cyanobacterium Aphanothece halophytica provide a rationale for salt tolerance

Serine hydroxymethyltransferase (SHMT) is a pyridoxal 5'-phosphate-dependent enzyme that plays a pivotal role in cellular one‑carbon metabolism. In plants and cyanobacteria, this enzyme is also involved in photorespiration and confers salt tolerance, as in the case of SHMT from the halophilic cyanobacterium Aphanothece halophytica (AhSHMT). We have characterized the catalytic properties of AhSHMT in different salt and pH conditions. Although the kinetic properties of AhSHMT correlate with those of the mesophilic orthologue from Escherichia coli, AhSHMT appears more catalytically efficient, especially in presence of salt. Our studies also reveal substrate inhibition, previously unobserved in AhSHMT. Furthermore, addition of the osmoprotectant glycine betaine under salt conditions has a distinct positive effect on AhSHMT activity. The crystal structures of AhSHMT in three forms, as internal aldimine, as external aldimine with the l-serine substrate, and as a covalent complex with malonate, give structural insights on the possible role of specific amino acid residues implicated in the halophilic features of AhSHMT. Importantly, we observed that overexpression of the gene encoding SHMT, independently from its origin, increases the capability of E. coli to grow in high salt conditions, suggesting that the catalytic activity of this enzyme in itself plays a fundamental role in salt tolerance.

Crystal structures of aminotransferases Aro8 and Aro9 from Candida albicans and structural insights into their properties

Aminotransferases catalyze reversibly the transamination reaction by a ping-pong bi-bi mechanism with pyridoxal 5'-phosphate (PLP) as a cofactor. Various aminotransferases acting on a range of substrates have been reported. Aromatic transaminases are able to catalyze the transamination reaction with both aromatic and acidic substrates. Two aminotransferases from C. albicans, Aro8p and Aro9p, have been identified recently, exhibiting different catalytic properties. To elucidate the multiple substrate recognition of the two enzymes we determined the crystal structures of an unliganded CaAro8p, a complex of CaAro8p with the PLP cofactor bound to a substrate, forming an external aldimine, CaAro9p with PLP in the form of internal aldimine, and CaAro9p with a mixture of ligands that have been interpreted as results of the enzymatic reaction. The crystal structures of both enzymes contains in the asymmetric unit a biologically relevant dimer of 55 kDa for CaAro8 and 59 kDa for CaAro9p protein subunits. The ability of the enzymes to process multiple substrates could be related to a feature of their architecture in which the active site resides on one subunit while the substrate-binding site is formed by a long loop extending from the other subunit of the dimeric molecule. The separation of the two functions to different chemical entities could facilitate the evolution of the substrate-binding part and allow it to be flexible without destabilizing the conservative catalytic mechanism.

Computational classification of MocR transcriptional regulators into subgroups as a support for experimental and functional characterization

MocR bacterial transcriptional regulators are a subfamily within the GntR family. The MocR proteins possess an N-terminal domain containing the winged Helix-Turn-Helix (wHTH) motif and a C-terminal domain whose architecture is homologous to the fold type-I pyridoxal 5'-phosphate (PLP) dependent enzymes and whose archetypical protein is aspartate aminotransferase (AAT). The ancestor of the fold type-I PLP dependent super-family is considered one of the earliest enzymes. The members of this super-family are the product of evolution which resulted in a diversified protein population able to catalyze a set of reactions on substrates often containing amino groups. The MocR regulators are activators or repressors of gene control within many metabolic pathways often involving PLP enzymes. This diversity implies that MocR specifically responds to different classes of effector molecules. Therefore, it is of interest to compare the AAT domains of MocR from six bacteria phyla. Multi dimensional scaling and cluster analyses suggested that at least three subgroups exist within the population that reflects functional specialization rather than taxonomic origin. The AAT-domains of the three clusters display variable degree of similarity to different fold type-I PLP enzyme families. The results support the hypothesis that independent fusion events generated at least three different MocR subgroups.

Crystal structure of the C-terminal domain of Bacillus subtilis GabR reveals a closed conformation by γ-aminobutyric acid binding, inducing transcriptional activation

Bacillus subtilis GabR (BsGabR) is involved in the γ-aminobutyric acid (GABA) catabolism as a transcriptional regulator, consisting of an N-terminal helix-turn-helix DNA-binding domain and a C-terminal aminotransferase-like (AT-like) domain. Research on the C-terminal AT-like domain of BsGabR (BsGabR-CTD) has focused on the interaction with GABA as an effector, but most its functional details remain unclear. To understand the underlying mechanism, we report the crystal structure of BsGabR-CTD in complex with pyridoxal 5'-phosphate (PLP) and GABA at 2.0 Å resolution. The structure of ligand-bound BsGabR-CTD revealed two distinct monomeric states in a homodimer. One subunit is a closed-form containing the PLP-GABA adduct, and the other subunit is a PLP-bound open-form. Our structural studies provide a detailed mechanism indicating that the open-to-closed transition by the binding of GABA induces the conformational rearrangement of BsGabR-CTD, which may trigger the activation of transcription.

The aspartate aminotransferase-like domain of Firmicutes MocR transcriptional regulators

Bacterial MocR transcriptional regulators possess an N-terminal DNA-binding domain containing a conserved helix-turn-helix module and an effector-binding and/or oligomerization domain at the C-terminus, homologous to fold type-I pyridoxal 5'-phosphate (PLP) enzymes. Since a comprehensive structural analysis of the MocR regulators is still missing, a comparisons of Firmicutes MocR sequences was undertook to contribute to the understanding of the structural characteristics of the C-terminal domain of these proteins, and to shed light on the structural and functional relationship with fold type-I PLP enzymes. Results of this work suggest the presence of at least three subgroups within the MocR sequences and provide a guide for rational site-directed mutagenesis studies aimed at deciphering the structure-function relationships in this new protein family.

A subfamily of PLP-dependent enzymes specialized in handling terminal amines

The present review focuses on a subfamily of pyridoxal phosphate (PLP)-dependent enzymes, belonging to the broader fold-type I structural group and whose archetypes can be considered ornithine δ-transaminase and γ-aminobutyrate transaminase. These proteins were originally christened "subgroup-II aminotransferases" (AT-II) but are very often referred to as "class-III aminotransferases". As names suggest, the subgroup includes mainly transaminases, with just a few interesting exceptions. However, at variance with most other PLP-dependent enzymes, catalysts in this subfamily seem specialized at utilizing substrates whose amino function is not adjacent to a carboxylate group. AT-II enzymes are widespread in biology and play mostly catabolic roles. Furthermore, today several transaminases in this group are being used as bioorganic tools for the asymmetric synthesis of chiral amines. We present an overview of the biochemical and structural features of these enzymes, illustrating how they are distinctive and how they compare with those of the other fold-type I enzymes. This article is part of a Special Issue entitled: Cofactor-dependent proteins: evolution, chemical diversity and bio-applications.

The crystal structure of archaeal serine hydroxymethyltransferase reveals idiosyncratic features likely required to withstand high temperatures

Serine hydroxymethyltransferases (SHMTs) play an essential role in one-carbon unit metabolism and are used in biomimetic reactions. We determined the crystal structure of free (apo) and pyridoxal-5'-phosphate-bound (holo) SHMT from Methanocaldococcus jannaschii, the first from a hyperthermophile, from the archaea domain of life and that uses H₄MPT as a cofactor, at 2.83 and 3.0 Å resolution, respectively. Idiosyncratic features were observed that are likely to contribute to structure stabilization. At the dimer interface, the C-terminal region folds in a unique fashion with respect to SHMTs from eubacteria and eukarya. At the active site, the conserved tyrosine does not make a cation-π interaction with an arginine like that observed in all other SHMT structures, but establishes an amide-aromatic interaction with Asn257, at a different sequence position. This asparagine residue is conserved and occurs almost exclusively in (hyper)thermophile SHMTs. This led us to formulate the hypothesis that removal of frustrated interactions (such as the Arg-Tyr cation-π interaction occurring in mesophile SHMTs) is an additional strategy of adaptation to high temperature. Both peculiar features may be tested by designing enzyme variants potentially endowed with improved stability for applications in biomimetic processes.

Type I pyridoxal 5'-phosphate dependent enzymatic domains embedded within multimodular nonribosomal peptide synthetase and polyketide synthase assembly lines

BACKGROUND

Pyridoxal 5'-phosphate (PLP)-dependent enzymes of fold type I, the most studied structural class of the PLP-dependent enzyme superfamily, are known to exist as stand-alone homodimers or homotetramers. These enzymes have been found also embedded in multimodular and multidomain assembly lines involved in the biosynthesis of polyketides (PKS) and nonribosomal peptides (NRPS). The aim of this work is to provide a proteome-wide view of the distribution and characteristics of type I domains covalently integrated in these assemblies in prokaryotes.

RESULTS

An ad-hoc Hidden Markov profile was calculated using a sequence alignment derived from a multiple structural superposition of distantly related PLP-enzymes of fold type I. The profile was utilized to scan the sequence databank and to collect the proteins containing at least one type I domain linked to a component of an assembly line in bacterial genomes. The domains adjacent to a carrier protein were further investigated. Phylogenetic analysis suggested the presence of four PLP-dependent families: Aminotran_3, Beta_elim_lyase and Pyridoxal_deC, occurring mainly within mixed NRPS/PKS clusters, and Aminotran_1_2 found mainly in PKS clusters. Sequence similarity to the reference PLP enzymes with solved structures ranged from 24 to 42% identity. Homology models were built for each representative type I domain and molecular docking simulations with putative substrates were carried out. Prediction of the protein-protein interaction sites evidenced that the surface regions of the type I domains embedded within multienzyme assemblies were different from those of the self-standing enzymes; these structural features appear to be required for productive interactions with the adjacent domains in a multidomain context.

CONCLUSIONS

This work provides a systematic view of the occurrence of type I domain within NRPS and PKS assembly lines and it predicts their structural characteristics using computational methods. Comparison with the corresponding stand-alone enzymes highlighted the common and different traits related to various aspects of their structure-function relationship. Therefore, the results of this work, on one hand contribute to the understanding of the functional and structural diversity of the PLP-dependent type I enzymes and, on the other, pave the way to further studies aimed at their applications in combinatorial biosynthesis.

Type I pyridoxal 5'-phosphate dependent enzymatic domains embedded within multimodular nonribosomal peptide synthetase and polyketide synthase assembly lines

BACKGROUND

Pyridoxal 5'-phosphate (PLP)-dependent enzymes of fold type I, the most studied structural class of the PLP-dependent enzyme superfamily, are known to exist as stand-alone homodimers or homotetramers. These enzymes have been found also embedded in multimodular and multidomain assembly lines involved in the biosynthesis of polyketides (PKS) and nonribosomal peptides (NRPS). The aim of this work is to provide a proteome-wide view of the distribution and characteristics of type I domains covalently integrated in these assemblies in prokaryotes.

RESULTS

An ad-hoc Hidden Markov profile was calculated using a sequence alignment derived from a multiple structural superposition of distantly related PLP-enzymes of fold type I. The profile was utilized to scan the sequence databank and to collect the proteins containing at least one type I domain linked to a component of an assembly line in bacterial genomes. The domains adjacent to a carrier protein were further investigated. Phylogenetic analysis suggested the presence of four PLP-dependent families: Aminotran_3, Beta_elim_lyase and Pyridoxal_deC, occurring mainly within mixed NRPS/PKS clusters, and Aminotran_1_2 found mainly in PKS clusters. Sequence similarity to the reference PLP enzymes with solved structures ranged from 24 to 42% identity. Homology models were built for each representative type I domain and molecular docking simulations with putative substrates were carried out. Prediction of the protein-protein interaction sites evidenced that the surface regions of the type I domains embedded within multienzyme assemblies were different from those of the self-standing enzymes; these structural features appear to be required for productive interactions with the adjacent domains in a multidomain context.

CONCLUSIONS

This work provides a systematic view of the occurrence of type I domain within NRPS and PKS assembly lines and it predicts their structural characteristics using computational methods. Comparison with the corresponding stand-alone enzymes highlighted the common and different traits related to various aspects of their structure-function relationship. Therefore, the results of this work, on one hand contribute to the understanding of the functional and structural diversity of the PLP-dependent type I enzymes and, on the other, pave the way to further studies aimed at their applications in combinatorial biosynthesis.

Structural mimicry between SLA/LP and Rickettsia surface antigens as a driver of autoimmune hepatitis: insights from an in silico study

BACKGROUND

Autoimmune hepatitis (AIH) is a chronic, progressive liver disease, characterized by continuing hepatocellular inflammation and necrosis. A subgroup of AIH patients presents specific autoantibodies to soluble liver antigen/liver-pancreas (SLA/LP) protein, which is regarded as a highly specific diagnostic marker. Autoantigenic SLA/LP peptides are targeted by CD4+ T cells, and restricted by the allele HLA-DRB1*03:01, which confers disease susceptibility in Europeans and Americans. A positively charged residue at position 71 has been indicated as critical for AIH susceptibility in all of the HLA alleles identified to date. Though the exact molecular mechanisms underlying pathogenesis of AIH are not clear, molecular mimicry between SLA/LP and viral/bacterial antigens has been invoked.

METHODS

The immunodominant region of SLA/LP was used as query in databank searches to identify statistically significant similarities with viral/bacterial peptides. Homology modeling and docking was used to investigate the potential interaction of HLA-DRB1*03:01 with the identified peptides. By molecular mechanics means, the interactions and energy of binding at the HLA binding site was also scrutinized.

RESULTS

A statistically significant structural similarity between the immunodominant regions of SLA/LP and a region of the surface antigen PS 120 from Rickettsia spp. has been detected. The interaction of the SLA/LP autoepitope and the corresponding Rickettsia sequence with the allele HLA-DRB1*03:01 has been simulated. The obtained results predict for both peptides a similar binding mode and affinity to HLA-DRB1*03:01. A "hot spot" of interaction between HLA-DRB1*03:01 and PS 120 is located at the P4 binding pocket, and is represented by a salt bridge involving Lys at position 71 of the HLA protein, and Glu 795 of PS120 peptide.

CONCLUSIONS

These findings strongly support the notion that a molecular mimicry mechanism can trigger AIH onset. CD4+ T cells recognizing peptides of SLA/LP could indeed cross-react with foreign Rickettsia spp. antigens. Finally, the same analysis suggests a molecular explanation for the importance of position 71 in conferring the susceptibility of the allele HLA-DRB1*03:01 to AIH. The lack of a positive charge at such position could prevent HLA alleles from binding the foreign peptides and triggering the molecular mimicry event.

Alanine racemase from Tolypocladium inflatum: a key PLP-dependent enzyme in cyclosporin biosynthesis and a model of catalytic promiscuity

Cyclosporin A, a cyclic peptide produced by the fungus Tolypocladium inflatum, is a widely employed immunosuppressant drug. Its biosynthesis is strictly dependent on the action of the pyridoxal 5'-phosphate-dependent enzyme alanine racemase, which produces the d-alanine incorporated in the cyclic peptide. This enzyme has a different fold with respect to bacterial alanine racemases. The interest elicited by T. inflatum alanine racemase not only relies on its biotechnological relevance, but also on its evolutionary and structural similarity to the promiscuous enzymes serine hydroxymethyltransferase and threonine aldolase. The three enzymes represent a model of divergent evolution from an ancestral enzyme that was able to catalyse all the reactions of the modern enzymes. A protocol to express and purify with high yield recombinant T. inflatum alanine racemase was developed. The catalytic properties of the enzyme were characterized. Similarly to serine hydroxymethyltransferase and threonine aldolase, T. inflatum alanine racemase was able to catalyse retroaldol cleavage and transamination reactions. This observation corroborates the hypothesis of the common evolutionary origin of these enzymes. A three-dimensional model of T. inflatum alanine racemase was constructed on the basis of threonine aldolase crystal structure. The model helped rationalise the experimental data and explain the catalytic properties of the enzymes.

The structural biology of enzymes involved in natural product glycosylation

The glycosylation of microbial natural products often dramatically influences the biological and/or pharmacological activities of the parental metabolite. Over the past decade, crystal structures of several enzymes involved in the biosynthesis and attachment of novel sugars found appended to natural products have emerged. In many cases, these studies have paved the way to a better understanding of the corresponding enzyme mechanism of action and have served as a starting point for engineering variant enzymes to facilitate to production of differentially-glycosylated natural products. This review specifically summarizes the structural studies of bacterial enzymes involved in biosynthesis of novel sugar nucleotides.

Comprehensively surveying structure and function of RING domains from Drosophila melanogaster

Using a complete set of RING domains from Drosophila melanogaster, all the solved RING domains and cocrystal structures of RING-containing ubiquitin-ligases (RING-E3) and ubiquitin-conjugating enzyme (E2) pairs, we analyzed RING domains structures from their primary to quarternary structures. The results showed that: i) putative orthologs of RING domains between Drosophila melanogaster and the human largely occur (118/139, 84.9%); ii) of the 118 orthologous pairs from Drosophila melanogaster and the human, 117 pairs (117/118, 99.2%) were found to retain entirely uniform domain architectures, only Iap2/Diap2 experienced evolutionary expansion of domain architecture; iii) 4 evolutionary structurally conserved regions (SCRs) are responsible for homologous folding of RING domains at the superfamily level; iv) besides the conserved Cys/His chelating zinc ions, 6 equivalent residues (4 hydrophobic and 2 polar residues) in the SCRs possess good-consensus and conservation- these 4 SCRs function in the structural positioning of 6 equivalent residues as determinants for RING-E3 catalysis; v) members of these RING proteins located nucleus, multiple subcellular compartments, membrane protein and mitochondrion are respectively 42 (42/139, 30.2%), 71 (71/139, 51.1%), 22 (22/139, 15.8%) and 4 (4/139, 2.9%); vi) CG15104 (Topors) and CG1134 (Mul1) in C3HC4, and CG3929 (Deltex) in C3H2C3 seem to display broader E2s binding profiles than other RING-E3s; vii) analyzing intermolecular interfaces of E2/RING-E3 complexes indicate that residues directly interacting with E2s are all from the SCRs in RING domains. Of the 6 residues, 2 hydrophobic ones contribute to constructing the conserved hydrophobic core, while the 2 hydrophobic and 2 polar residues directly participate in E2/RING-E3 interactions. Based on sequence and structural data, SCRs, conserved equivalent residues and features of intermolecular interfaces were extracted, highlighting the presence of a nucleus for RING domain fold and formation of catalytic core in which related residues and regions exhibit preferential evolutionary conservation.

Linkage between the bacterial acid stress and stringent responses: the structure of the inducible lysine decarboxylase

The Escherichia coli inducible lysine decarboxylase, LdcI/CadA, together with the inner-membrane lysine-cadaverine antiporter, CadB, provide cells with protection against mild acidic conditions (pH∼5). To gain a better understanding of the molecular processes underlying the acid stress response, the X-ray crystal structure of LdcI was determined. The structure revealed that the protein is an oligomer of five dimers that associate to form a decamer. Surprisingly, LdcI was found to co-crystallize with the stringent response effector molecule ppGpp, also known as the alarmone, with 10 ppGpp molecules in the decamer. ppGpp is known to mediate the stringent response, which occurs in response to nutrient deprivation. The alarmone strongly inhibited LdcI enzymatic activity. This inhibition is important for modulating the consumption of lysine in cells during acid stress under nutrient limiting conditions. Hence, our data provide direct evidence for a link between the bacterial acid stress and stringent responses.

Modular prediction of protein structural classes from sequences of twilight-zone identity with predicting sequences

Background

Knowledge of structural class is used by numerous methods for identification of structural/functional characteristics of proteins and could be used for the detection of remote homologues, particularly for chains that share twilight-zone similarity. In contrast to existing sequence-based structural class predictors, which target four major classes and which are designed for high identity sequences, we predict seven classes from sequences that share twilight-zone identity with the training sequences.

Results

The proposed MODular Approach to Structural class prediction (MODAS) method is unique as it allows for selection of any subset of the classes. MODAS is also the first to utilize a novel, custom-built feature-based sequence representation that combines evolutionary profiles and predicted secondary structure. The features quantify information relevant to the definition of the classes including conservation of residues and arrangement and number of helix/strand segments. Our comprehensive design considers 8 feature selection methods and 4 classifiers to develop Support Vector Machine-based classifiers that are tailored for each of the seven classes. Tests on 5 twilight-zone and 1 high-similarity benchmark datasets and comparison with over two dozens of modern competing predictors show that MODAS provides the best overall accuracy that ranges between 80% and 96.7% (83.5% for the twilight-zone datasets), depending on the dataset. This translates into 19% and 8% error rate reduction when compared against the best performing competing method on two largest datasets. The proposed predictor provides accurate predictions at 58% accuracy for membrane proteins class, which is not considered by majority of existing methods, in spite that this class accounts for only 2% of the data. Our predictive model is analyzed to demonstrate how and why the input features are associated with the corresponding classes.

Conclusions

The improved predictions stem from the novel features that express collocation of the secondary structure segments in the protein sequence and that combine evolutionary and secondary structure information. Our work demonstrates that conservation and arrangement of the secondary structure segments predicted along the protein chain can successfully predict structural classes which are defined based on the spatial arrangement of the secondary structures. A web server is available at http://biomine.ece.ualberta.ca/MODAS/.

Protein structure classification based on conserved hydrophobic residues

Protein folding is frequently guided by local residue interactions that form clusters in the protein core. The interactions between residue clusters serve as potential nucleation sites in the folding process. Evidence postulates that the residue interactions are governed by the hydrophobic propensities that the residues possess. An array of hydrophobicity scales has been developed to determine the hydrophobic propensities of residues under different environmental conditions. In this work, we propose a graph-theory-based data mining framework to extract and isolate protein structural features that sustain invariance in evolutionary-related proteins, through the integrated analysis of five well-known hydrophobicity scales over the 3D structure of proteins. We hypothesize that proteins of the same homology contain conserved hydrophobic residues and exhibit analogous residue interaction patterns in the folded state. The results obtained demonstrate that discriminatory residue interaction patterns shared among proteins of the same family can be employed for both the structural and the functional annotation of proteins. We obtained on the average 90 percent accuracy in protein classification with a significantly small feature vector compared to previous results in the area. This work presents an elaborate study, as well as validation evidence, to illustrate the efficacy of the method and the correctness of results reported.

Accommodation of GDP-linked sugars in the active site of GDP-perosamine synthase

Perosamine (4-amino-4,6-dideoxy- d-mannose), or its N-acetylated form, is one of several dideoxy sugars found in the O-antigens of such infamous Gram-negative bacteria as Vibrio cholerae O1 and Escherichia coli O157:H7. It is added to the bacterial O-antigen via a nucleotide-linked version, namely GDP-perosamine. Three enzymes are required for the biosynthesis of GDP-perosamine starting from mannose 1-phosphate. The focus of this investigation is GDP-perosamine synthase from Caulobacter crescentus, which catalyzes the final step in GDP-perosamine synthesis, the conversion of GDP-4-keto-6-deoxymannose to GDP-perosamine. The enzyme is PLP-dependent and belongs to the aspartate aminotransferase superfamily. It contains the typically conserved active site lysine residue, which forms a Schiff base with the PLP cofactor. Two crystal structures were determined for this investigation: a site-directed mutant protein (K186A) complexed with GDP-perosamine and the wild-type enzyme complexed with an unnatural ligand, GDP-3-deoxyperosamine. These structures, determined to 1.6 and 1.7 A resolution, respectively, revealed the manner in which products, and presumably substrates, are accommodated within the active site pocket of GDP-perosamine synthase. Additional kinetic analyses using both the natural and unnatural substrates revealed that the K m for the unnatural substrate was unperturbed relative to that of the natural substrate, but the k cat was lowered by a factor of approximately 200. Taken together, these studies shed light on why GDP-perosamine synthase functions as an aminotransferase whereas another very similar PLP-dependent enzyme, GDP-4-keto-6-deoxy- d-mannose 3-dehydratase or ColD, catalyzes a dehydration reaction using the same substrate.

SCPRED: accurate prediction of protein structural class for sequences of twilight-zone similarity with predicting sequences

Background

Protein structure prediction methods provide accurate results when a homologous protein is predicted, while poorer predictions are obtained in the absence of homologous templates. However, some protein chains that share twilight-zone pairwise identity can form similar folds and thus determining structural similarity without the sequence similarity would be desirable for the structure prediction. The folding type of a protein or its domain is defined as the structural class. Current structural class prediction methods that predict the four structural classes defined in SCOP provide up to 63% accuracy for the datasets in which sequence identity of any pair of sequences belongs to the twilight-zone. We propose SCPRED method that improves prediction accuracy for sequences that share twilight-zone pairwise similarity with sequences used for the prediction.

Results

SCPRED uses a support vector machine classifier that takes several custom-designed features as its input to predict the structural classes. Based on extensive design that considers over 2300 index-, composition- and physicochemical properties-based features along with features based on the predicted secondary structure and content, the classifier's input includes 8 features based on information extracted from the secondary structure predicted with PSI-PRED and one feature computed from the sequence. Tests performed with datasets of 1673 protein chains, in which any pair of sequences shares twilight-zone similarity, show that SCPRED obtains 80.3% accuracy when predicting the four SCOP-defined structural classes, which is superior when compared with over a dozen recent competing methods that are based on support vector machine, logistic regression, and ensemble of classifiers predictors.

Conclusion

The SCPRED can accurately find similar structures for sequences that share low identity with sequence used for the prediction. The high predictive accuracy achieved by SCPRED is attributed to the design of the features, which are capable of separating the structural classes in spite of their low dimensionality. We also demonstrate that the SCPRED's predictions can be successfully used as a post-processing filter to improve performance of modern fold classification methods.

Molecular architecture of DesI: a key enzyme in the biosynthesis of desosamine

Desosamine is a 3-(dimethylamino)-3,4,6-trideoxyhexose found, for example, in such macrolide antibiotics as erthyromycin, azithromycin, and clarithromycin. The efficacies of these macrolide antibiotics are markedly reduced in the absence of desosamine. In the bacterium Streptomyces venezuelae, six enzymes are required for the production of dTDP-desosamine. The focus of this X-ray crystallographic analysis is the third enzyme in the pathway, a PLP-dependent aminotransferase referred to as DesI. The structure of DesI was solved in complex with its product, dTDP-4-amino-4,6-dideoxyglucose, to a nominal resolution of 2.1 A. Each subunit of the dimeric enzyme contains 12 alpha-helices and 14 beta-strands. Three cis-peptides are observed in each subunit, Phe 330, Pro 332, and Pro 339. The two active sites of the enzyme are located in clefts at the subunit/subunit interface. Electron density corresponding to the bound product clearly demonstrates a covalent bond between the amino group of the product and C-4' of the PLP cofactor. Interestingly, there are no hydrogen-bonding interactions between the protein and the dideoxyglucosyl group of the product (within 3.2 A). The only other sugar-modifying aminotransferase whose structure is known in the presence of product is PseC from Helicobacter pylori. This enzyme, as opposed to DesI, catalyzes amino transfer to the axial position of the sugar. A superposition of the two active sites for these proteins reveals that the major differences in ligand binding occur in the orientations of the deoxyglucosyl and phosphoryl groups. Indeed, the nearly 180 degrees difference in hexose orientation explains the equatorial versus axial amino transfer exhibited by DesI and PseC, respectively.

Molecular architecture of DesV from Streptomyces venezuelae: a PLP-dependent transaminase involved in the biosynthesis of the unusual sugar desosamine

Desosamine is a 3-(dimethylamino)-3,4,6-trideoxyhexose found in certain macrolide antibiotics such as the commonly prescribed erythromycin. Six enzymes are required for its biosynthesis in Streptomyces venezuelae. The focus of this article is DesV, which catalyzes the PLP-dependent replacement of a 3-keto group with an amino functionality in the fifth step of the pathway. For this study the three-dimensional structures of both the internal aldimine and the ketimine intermediate with glutamate were determined to 2.05 A resolution. DesV is a homodimer with each subunit containing 12 alpha-helical regions and 12 beta-strands that together form three layers of sheet. The structure of the internal aldimine demonstrates that the PLP-cofactor is held in place by residues contributed from both subunits (Asp 164 and Gln 167 from Subunit I and Tyr 221 and Asn 235 from Subunit II). When the ketimine intermediate is present in the active site, the loop defined by Gln 225 to Ser 228 from Subunit II closes down upon the active site. The structure of DesV is similar to another sugar-modifying enzyme referred to as PseC. This enzyme is involved in the biosynthesis of pseudaminic acid, which is a sialic acid-like nonulosonate found in the flagellin of Helicobacter pylori. In the case of PseC, however, the amino group is transferred to the C-4 rather than the C-3 position. Details concerning the structural analysis of DesV and a comparison of its molecular architecture to that of PseC are presented.

Crystal structure of Homo sapiens kynureninase

Kynureninase is a member of a large family of catalytically diverse but structurally homologous pyridoxal 5'-phosphate (PLP) dependent enzymes known as the aspartate aminotransferase superfamily or alpha-family. The Homo sapiens and other eukaryotic constitutive kynureninases preferentially catalyze the hydrolytic cleavage of 3-hydroxy-l-kynurenine to produce 3-hydroxyanthranilate and l-alanine, while l-kynurenine is the substrate of many prokaryotic inducible kynureninases. The human enzyme was cloned with an N-terminal hexahistidine tag, expressed, and purified from a bacterial expression system using Ni metal ion affinity chromatography. Kinetic characterization of the recombinant enzyme reveals classic Michaelis-Menten behavior, with a Km of 28.3 +/- 1.9 microM and a specific activity of 1.75 micromol min-1 mg-1 for 3-hydroxy-dl-kynurenine. Crystals of recombinant kynureninase that diffracted to 2.0 A were obtained, and the atomic structure of the PLP-bound holoenzyme was determined by molecular replacement using the Pseudomonas fluorescens kynureninase structure (PDB entry 1qz9) as the phasing model. A structural superposition with the P. fluorescens kynureninase revealed that these two structures resemble the "open" and "closed" conformations of aspartate aminotransferase. The comparison illustrates the dynamic nature of these proteins' small domains and reveals a role for Arg-434 similar to its role in other AAT alpha-family members. Docking of 3-hydroxy-l-kynurenine into the human kynureninase active site suggests that Asn-333 and His-102 are involved in substrate binding and molecular discrimination between inducible and constitutive kynureninase substrates.

A retro-evolution study of CDP-6-deoxy-D-glycero-L-threo-4-hexulose-3-dehydrase (E1) from Yersinia pseudotuberculosis: implications for C-3 deoxygenation in the biosynthesis of 3,6-dideoxyhexoses

CDP-6-deoxy-l-threo-d-glycero-4-hexulose-3-dehydrase (E1), which catalyzes C-3 deoxygenation of CDP-4-keto-6-deoxyglucose in the biosynthesis of 3,6-dideoxyhexoses, shares a modest sequence identity with other B6-dependent enzymes, albeit with two important distinctions. It is a rare example of a B6-dependent enzyme that harbors a [2Fe-2S] cluster, and a highly conserved lysine that serves as an anchor for PLP in most B6-dependent enzymes is replaced by histidine at position 220 in E1. Since alteration of His220 to a lysine residue may produce a putative progenitor of E1, the H220K mutant was constructed and tested for the ability to process the predicted substrate, CDP-4-amino-4,6-dideoxyglucose, using PLP as the coenzyme. Our data showed that H220K-E1 has no dehydrase activity, but can act as a PLP-dependent transaminase. However, the reaction is not catalytic since PLP cannot be regenerated during turnover. Reported herein are the results of this investigation and the implications for the role of His220 in the catalytic mechanism of E1.

Characterizing conserved structural contacts by pair-wise relative contacts and relative packing groups

To adequately deal with the inherent complexity of interactions between protein side-chains, we develop and describe here a novel method for characterizing protein packing within a fold family. Instead of approaching side-chain interactions absolutely from one residue to another, we instead consider the relative interactions of contacting residue pairs. The basic element, the pair-wise relative contact, is constructed from a sequence alignment and contact analysis of a set of structures and consists of a cluster of similarly oriented, interacting, side-chain pairs. To demonstrate this construct's usefulness in analyzing protein structure, we used the pair-wise relative contacts to analyze two sets of protein structures as defined by SCOP: the diverse globin-like superfamily (126 structures) and the more uniform heme binding globin family (a 94 structure subset of the globin-like superfamily). The superfamily structure set produced 1266 unique pair-wise relative contacts, whereas the family structure subset gave 1001 unique pair-wise relative contacts. For both sets, we show that these constructs can be used to accurately and automatically differentiate between fold classes. Furthermore, these pair-wise relative contacts correlate well with sequence identity and thus provide a direct relationship between changes in sequence and changes in structure. To capture the complexity of protein packing, these pair-wise relative contacts can be superimposed around a single residue to create a multi-body construct called a relative packing group. Construction of convex hulls around the individual packing groups provides a measure of the variation in packing around a residue and defines an approximate volume of space occupied by the groups interacting with a residue. We find that these relative packing groups are useful in understanding the structural quality of sequence or structure alignments. Moreover, they provide context to calculate a value for structural randomness, which is important in properly assessing the quality of a structural alignment. The results of this study provide the framework for future analysis for correlating sequence changes to specific structure changes.

Conversion of 5-aminolevulinate synthase into a more active enzyme by linking the two subunits: spectroscopic and kinetic properties

The two active sites of dimeric 5-aminolevulinate synthase (ALAS), a pyridoxal 5'-phosphate (PLP)-dependent enzyme, are located on the subunit interface with contribution of essential amino acids from each subunit. Linking the two subunits into a single polypeptide chain dimer (2XALAS) yielded an enzyme with an approximate sevenfold greater turnover number than that of wild-type ALAS. Spectroscopic and kinetic properties of 2XALAS were investigated to explore the differences in the coenzyme structure and kinetic mechanism relative to those of wild-type ALAS that confer a more active enzyme. The absorption spectra of both ALAS and 2XALAS had maxima at 410 and 330 nm, with a greater A(410)/A(330) ratio at pH approximately 7.5 for 2XALAS. The 330 nm absorption band showed an intense fluorescence at 385 nm but not at 510 nm, indicating that the 330 nm absorption species is the substituted aldamine rather than the enolimine form of the Schiff base. The 385 nm emission intensity increased with increasing pH with a single pK of approximately 8.5 for both enzymes, and thus the 410 and 330 nm absorption species were attributed to the ketoenamine and substituted aldamine, respectively. Transient kinetic analysis of the formation and decay of the quinonoid intermediate EQ(2) indicated that, although their rates were similar in ALAS and 2XALAS, accumulation of this intermediate was greater in the 2XALAS-catalyzed reaction. Collectively, these results suggest that ketoenamine is the active form of the coenzyme and forms a more prominent coenzyme structure in 2XALAS than in ALAS at pH approximately 7.5.