Fold recognition analysis of glycosyltransferase families: further members of structural superfamilies

Glycosyltransferases (GTs) are diverse enzymes organized into 65 families. X-ray crystallography and in silico studies have shown many of these to belong to two structural superfamilies: GT-A and GT-B. Through application of fold recognition and iterated sequence searches, we demonstrate that families 60, 62, and 64 may also be grouped into the GT-A fold superfamily. Analysis of conserved acidic residues suggests that catalytic sites are better conserved in superfamily GT-B than in GT-A. Although 26% and 29% of GT families may now be confidently placed in superfamilies GT-A and GT-B, respectively, the remaining 45% of families bear no discernible resemblance to either superfamily, which, given the sensitivity of modern fold recognition methods, suggests the existence of novel structural scaffolds associated with GT activity. Furthermore, bioinformatics studies indicate the apparent ease with which mechanism-inverting or retaining-may change during evolution.     

X-ray crystal structure of leukocyte type core 2 beta1,6-N-acetylglucosaminyltransferase. Evidence for a convergence of metal ion-independent glycosyltransferase mechanism

Leukocyte type core 2 beta1,6-N-acetylglucosaminyltransferase (C2GnT-L) is a key enzyme in the biosynthesis of branched O-glycans. It is an inverting, metal ion-independent family 14 glycosyltransferase that catalyzes the formation of the core 2 O-glycan (Galbeta1-3[GlcNAcbeta1-6]GalNAc-O-Ser/Thr) from its donor and acceptor substrates, UDP-GlcNAc and the core 1 O-glycan (Galbeta1-3GalNAc-O-Ser/Thr), respectively. Reported here are the x-ray crystal structures of murine C2GnT-L in the absence and presence of the acceptor substrate Galbeta1-3GalNAc at 2.0 and 2.7A resolution, respectively. C2GnT-L was found to possess the GT-A fold; however, it lacks the characteristic metal ion binding DXD motif. The Galbeta1-3GalNAc complex defines the determinants of acceptor substrate binding and shows that Glu-320 corresponds to the structurally conserved catalytic base found in other inverting GT-A fold glycosyltransferases. Comparison of the C2GnT-L structure with that of other GT-A fold glycosyltransferases further suggests that Arg-378 and Lys-401 serve to electrostatically stabilize the nucleoside diphosphate leaving group, a role normally played by metal ion in GT-A structures. The use of basic amino acid side chains in this way is strikingly similar to that seen in a number of metal ion-independent GT-B fold glycosyltransferases and suggests a convergence of catalytic mechanism shared by both GT-A and GT-B fold glycosyltransferases.     

Structures and mechanisms of glycosyltransferases

Glycosyltransferases (GTs) catalyze the transfer of a sugar moiety from an activated donor sugar onto saccharide and nonsaccharide acceptors. A sequence-based classification spreads GTs in many families thus reflecting the variety of molecules that can be used as acceptors. In contrast, this enzyme family is characterized by a more conserved three-dimensional architecture. Until recently, only two different folds (GT-A and GT-B) have been identified for solved crystal structures. The recent report of a structure for a bacterial sialyltransferase allows the definition of a new fold family. Progress in the elucidation of the structures and mechanisms of GTs are discussed in this review. To accommodate the growing number of crystal structures, we created the 3D-Glycosyltransferase database to gather structural information concerning this class of enzymes.     

A complementary bioinformatics approach to identify potential plant cell wall glycosyltransferase-encoding genes

Plant cell wall (CW) synthesizing enzymes can be divided into the glycan (i.e. cellulose and callose) synthases, which are multimembrane spanning proteins located at the plasma membrane, and the glycosyltransferases (GTs), which are Golgi localized single membrane spanning proteins, believed to participate in the synthesis of hemicellulose, pectin, mannans, and various glycoproteins. At the Carbohydrate-Active enZYmes (CAZy) database where e.g. glucoside hydrolases and GTs are classified into gene families primarily based on amino acid sequence similarities, 415 Arabidopsis GTs have been classified. Although much is known with regard to composition and fine structures of the plant CW, only a handful of CW biosynthetic GT genes-all classified in the CAZy system-have been characterized. In an effort to identify CW GTs that have not yet been classified in the CAZy database, a simple bioinformatics approach was adopted. First, the entire Arabidopsis proteome was run through the Transmembrane Hidden Markov Model 2.0 server and proteins containing one or, more rarely, two transmembrane domains within the N-terminal 150 amino acids were collected. Second, these sequences were submitted to the SUPERFAMILY prediction server, and sequences that were predicted to belong to the superfamilies NDP-sugartransferase, UDP-glycosyltransferase/glucogen-phosphorylase, carbohydrate-binding domain, Gal-binding domain, or Rossman fold were collected, yielding a total of 191 sequences. Fifty-two accessions already classified in CAZy were discarded. The resulting 139 sequences were then analyzed using the Three-Dimensional-Position-Specific Scoring Matrix and mGenTHREADER servers, and 27 sequences with similarity to either the GT-A or the GT-B fold were obtained. Proof of concept of the present approach has to some extent been provided by our recent demonstration that two members of this pool of 27 non-CAZy-classified putative GTs are xylosyltransferases involved in synthesis of pectin rhamnogalacturonan II (J. Egelund, B.L. Petersen, A. Faik, M.S. Motawia, C.E. Olsen, T. Ishii, H. Clausen, P. Ulvskov, and N. Geshi, unpublished data).     

Crystal structures of sialyltransferase from Photobacterium damselae

Sialyltransferase structures fall into either GT-A or GT-B glycosyltransferase fold. Some sialyltransferases from the Photobacterium genus have been shown to contain an additional N-terminal immunoglobulin (Ig)-like domain. Photobacterium damselae α2–6-sialyltransferase has been used efficiently in enzymatic and chemoenzymatic synthesis of α2–6-linked sialosides. Here we report three crystal structures of this enzyme. Two structures with and without a donor substrate analog CMP-3F(a)Neu5Ac contain an immunoglobulin (Ig)-like domain and adopt the GT-B sialyltransferase fold. The binary structure reveals a non-productive pre-Michaelis complex, which are caused by crystal lattice contacts that prevent the large conformational changes. The third structure lacks the Ig-domain. Comparison of the three structures reveals small inherent flexibility between the two Rossmann-like domains of the GT-B fold.     

常见致病菌糖基转移酶的结构与功能

糖基转移酶(glycosyltransferases,GTs)将糖基从活化的供体转移到糖、脂、蛋白质和核酸等受体,其参与的蛋白质糖基化是最重要的翻译后修饰(post-translational modifications,PTMs)之一。近年来越来越多的研究证明,糖基转移酶与致病菌毒力密切相关,在致病菌的黏附、免疫逃逸和定殖等生物学过程中发挥关键作用。目前,已鉴定的糖基转移酶根据其蛋白质三维结构特征分为3种类型GT-A、GT-B和GT-C,其中常见的是GT-A和GT-B型。在致病菌中发挥黏附功能的糖基转移酶,在结构上属于GT-B或GT-C型,对致病菌表面蛋白质(黏附蛋白、自转运蛋白等)进行糖基化修饰,在致病菌黏附、生物被膜的形成和毒力机制发挥具有重要作用。糖基转移酶不仅参与致病菌黏附这一感染初始过程,其中属于GT-A型的一类致病菌糖基转移酶会进入宿主细胞,通过糖基化宿主蛋白质影响宿主信号传导、蛋白翻译和免疫应答等生物学功能。本文就常见致病菌糖基转移酶的结构及其糖基化在致病机制中的作用进行综述,着重介绍了特异性糖基化高分子量(high-molecular-weight,HMW)黏附蛋白的糖基转移酶、针对富丝氨酸重复蛋白(serine-rich repeat proteins,SRRP)糖基化修饰的糖基转移酶、细菌自转运蛋白庚糖基转移酶(bacterial autotransporter heptosyltransferase,BAHT)家族、N-糖基化蛋白质系统和进入宿主细胞发挥毒力作用的大型梭菌细胞毒素、军团菌(Legionella)葡萄糖基转移酶以及肠杆菌科的效应子NleB。为揭示致病菌中糖基转移酶致病机制的系统性研究提供参考,为未来致病菌的诊断、药物设计研发以及疫苗开发等提供科学依据和思路。     

Identification of Putative Rhamnogalacturonan-II Specific Glycosyltransferases in Arabidopsis Using a Combination of Bioinformatics Approaches

Rhamnogalacturonan-II (RG-II) is a complex plant cell wall polysaccharide that is composed of an α(1,4)-linked homogalacturonan backbone substituted with four side chains. It exists in the cell wall in the form of a dimer that is cross-linked by a borate di-ester. Despite its highly complex structure, RG-II is evolutionarily conserved in the plant kingdom suggesting that this polymer has fundamental functions in the primary wall organisation. In this study, we have set up a bioinformatics strategy aimed at identifying putative glycosyltransferases (GTs) involved in RG-II biosynthesis. This strategy is based on the selection of candidate genes encoding type II membrane proteins that are tightly coexpressed in both rice and Arabidopsis with previously characterised genes encoding enzymes involved in the synthesis of RG-II and exhibiting an up-regulation upon isoxaben treatment. This study results in the final selection of 26 putative Arabidopsis GTs, including 10 sequences already classified in the CAZy database. Among these CAZy sequences, the screening protocol allowed the selection of α-galacturonosyltransferases involved in the synthesis of α4-GalA oligogalacturonides present in both homogalacturonans and RG-II, and two sialyltransferase-like sequences previously proposed to be involved in the transfer of Kdo and/or Dha on the pectic backbone of RG-II. In addition, 16 non-CAZy GT sequences were retrieved in the present study. Four of them exhibited a GT-A fold. The remaining sequences harbored a GT-B like fold and a fucosyltransferase signature. Based on homologies with glycosyltransferases of known functions, putative roles in the RG-II biosynthesis are proposed for some GT candidates.     

Sequence landmark patterns identify and characterize protein families

The three-dimensional structures of homologous proteins are usually conserved during evolution, as are critical residues in a few short sequence motifs that often constitute the active site in enzymes. The precise spatial organization of such sites depends on the lengths and positions of the secondary structural elements connecting the motifs. We show how members of protein superfamilies, such as kinesins, myosins, and G(alpha) subunits of trimeric G proteins, are identified and classed by simply counting the number of amino acid residues between important sequence motifs in their nucleotide triphosphate-hydrolyzing domains. Subfamily-specific landmark patterns (motif to motif scores) are principally due to inserts and gaps in surface loops. Unusual protein sequences and possible sequence prediction errors are detected.     

Crystal structure of mammalian alpha1,6-fucosyltransferase, FUT8

Mammalian alpha1,6-fucosyltransferase (FUT8) catalyses the transfer of a fucose residue from a donor substrate, guanosine 5'-diphosphate-beta-L-fucose to the reducing terminal N-acetylglucosamine (GlcNAc) of the core structure of an asparagine-linked oligosaccharide. Alpha1,6-fucosylation, also referred to as core fucosylation, plays an essential role in various pathophysiological events. Our group reported that FUT8 null mice showed severe growth retardation and emphysema-like lung-destruction as a result of the dysfunction of epidermal growth factor and transforming growth factor-beta receptors. To elucidate the molecular basis of FUT8 with respect to pathophysiology, the crystal structure of human FUT8 was determined at 2.6 A resolution. The overall structure of FUT8 was found to consist of three domains: an N-terminal coiled-coil domain, a catalytic domain, and a C-terminal SH3 domain. The catalytic region appears to be similar to GT-B glycosyltransferases rather than GT-A. The C-terminal part of the catalytic domain of FUT8 includes a Rossmann fold with three regions that are conserved in alpha1,6-, alpha1,2-, and protein O-fucosyltransferases. The SH3 domain of FUT8 is similar to other SH3 domain-containing proteins, although the significance of this domain remains to be elucidated. The present findings of FUT8 suggest that the conserved residues in the three conserved regions participate in the Rossmann fold and act as the donor binding site, or in catalysis, thus playing key roles in the fucose-transferring reaction.     

Phylogeny-guided characterization of glycosyltransferases for epothilone glycosylation

Glycosylation of natural products can influence their pharmacological properties, and efficient glycosyltransferases (GTs) are critical for this purpose. The polyketide epothilones are potent anti-tumour compounds, and YjiC is the only reported GT for the glycosylation of epothilone. In this study, we phylogenetically analysed 8261 GTs deposited in CAZy database and revealed that YjiC locates in a subbranch of the Macrolide I group, forming the YjiC-subbranch with 160 GT sequences. We demonstrated that the YjiC-subbranch GTs are normally efficient in epothilone glycosylation, but some showed low glycosylation activities. Sequence alignment of YjiC-subbranch showed that the 66th and 77th amino acid residues, which were close to the catalytic cavity in molecular docking model, were conserved in five high-active GTs (Q66 and P77) but changed in two low-efficient GTs. Site-directed residues swapping at the two positions in the two low-active GTs (BssGT and BamGT) and the high-active GT BsGT-1 demonstrated that the two amino acid residues played an important role in the catalytic efficiency of epothilone glycosylation. This study highlights that the potent GTs for appointed compounds are phylogenetically grouped with conserved residues for the catalytic efficiency.     

Protein folds propelled by diversity

Many proteins involved in key biological processes are modular in nature. A group of these, the beta-propeller proteins, fold by packing 4-stranded beta-sheets in a circular array. The members of this group are increasingly numerous and, although their modular building blocks all preserve the same basic conformation, they do not have similar sequences. These proteins have extreme functional and phylogenetic diversity. Here, features of the beta-propeller fold are reviewed through comparisons of available structural coordinates. Structure-based sequence alignments combined with analyses of superpositions of individual modular units reveal conserved general features such as hydrogen bonds, beta-turns and positions of hydrophobic contacts. The lack of significant sequence identity is compensated by sets of interactions which stabilise the fold differently in distinct structures. Re-occurring aspartates make contacts to exposed backbone amides in turns or peptide connections within the same sheet. The sole factor responsible for the number of sheets that assemble in the array is the size of the hydrophobic residues that pack into the cores between the sheets. Whilst there is no overall sequence conservation, it may be possible to detect new members of this fold through sequence searches that take into account the repeated nature of the modular assembly as well as the positions of hydrophobic residues and H-bonding side chains.     

植物尿苷二磷酸糖基转移酶超家族晶体结构

糖基转移酶(Glycosyltransferases,GTs)催化的糖基化反应几乎是植物中最为重要的反应。GTs家族1中的植物UGTs(UDP-dependent glycosyltransferases)成员主要运用尿苷二磷酸活化的糖作为糖基供体,因其成员众多、生物功能多样,仅仅通过序列比较和进化分析不能够精确预测其复杂的底物专一性和特有的催化机制,需要后续生化实验的进一步验证。文中主要总结了目前在蛋白结构数据库(Protein Data Bank,PDB)中报道的5种植物UGTs的晶体三维结构和定点突变功能研究进展。详细介绍了植物UGTs整体结构的特点以及蛋白与底物相互作用的细节,为更有效地生化定性UGTs以便深入理解底物专一性提供了有力的工具,从而为植物UGTs在酶工程和基因工程中的应用奠定基础。     

Functional divergence after gene duplication and sequence-structure relationship: a case study of G-protein alpha subunits

In this article, we use animal G-protein alpha subunit family as an example to illustrate a comprehensive analytical pipeline for detecting different types of functional divergence of protein families, which is phylogeny-dependent, combined with ancestral sequence inference and available protein structure information. In particular, we focus on (i) Type-I functional divergence, or site-specific rate shift, as typically exemplified by amino acid residue highly conserved in a subset of homologous genes but highly variable in a different subset of homologous genes, and (ii) Type-II functional divergence, or the shift of cluster-specific amino acid property, as exemplified by a radical shift of amino acid property between duplicate genes, which is otherwise evolutionally conserved. We utilized the software DIVERGE2 to carry out these analyses. In the case of G-protein alpha subunit gene family, we have predicted amino acid residues that are related to either Type-I or Type-II functional divergence. The inferred ancestral sequences for these sites are helpful to explore the trends of functional divergence. Finally, these predicted residues are mapped to the protein structures to test whether these residues may have 3D structure or solvent accessibility preference.     

Crystal structure of a UDP-glucose-specific glycosyltransferase from a Mycobacterium species

Glycosyltransferases (GTs) are a large and ubiquitous family of enzymes that specifically transfer sugar moieties to a range of substrates. Mycobacterium tuberculosis contains a large number of GTs, many of which are implicated in cell wall synthesis, yet the majority of these GTs remain poorly characterized. Here, we report the high resolution crystal structures of an essential GT (MAP2569c) from Mycobacterium avium subsp. paratuberculosis (a close homologue of Rv1208 from M. tuberculosis) in its apo- and ligand-bound forms. The structure adopted the GT-A fold and possessed the characteristic DXD motif that coordinated an Mn(2+) ion. Atypical of most GTs characterized to date, MAP2569c exhibited specificity toward the donor substrate, UDP-glucose. The structure of this ligated complex revealed an induced fit binding mechanism and provided a basis for this unique specificity. Collectively, the structural features suggested that MAP2569c may adopt a "retaining" enzymatic mechanism, which has implications for the classification of other GTs in this large superfamily.     

Current trends in the structure-activity relationships of sialyltransferases

Sialyltransferases (STs) represent an important group of enzymes that transfer N-acetylneuraminic acid (Neu5Ac) from cytidine monophosphate-Neu5Ac to various acceptor substrates. In higher animals, sialylated oligosaccharide structures play crucial roles in many biological processes but also in diseases, notably in microbial infection and cancer. Cell surface sialic acids have also been found in a few microorganisms, mainly pathogenic bacteria, and their presence is often associated with virulence. STs are distributed into five different families in the CAZy database (http://www.cazy.org/). On the basis of crystallographic data available for three ST families and fold recognition analysis for the two other families, STs can be grouped into two structural superfamilies that represent variations of the canonical glycosyltransferase (GT-A and GT-B) folds. These two superfamilies differ in the nature of their active site residues, notably the catalytic base (a histidine or an aspartate residue). The observed structural and functional differences strongly suggest that these two structural superfamilies have evolved independently.     

Efficient prediction of nucleic acid binding function from low-resolution protein structures

Structural genomics projects as well as ab initio protein structure prediction methods provide structures of proteins with no sequence or fold similarity to proteins with known functions. These are often low-resolution structures that may only include the positions of C alpha atoms. We present a fast and efficient method to predict DNA-binding proteins from just the amino acid sequences and low-resolution, C alpha-only protein models. The method uses the relative proportions of certain amino acids in the protein sequence, the asymmetry of the spatial distribution of certain other amino acids as well as the dipole moment of the molecule. These quantities are used in a linear formula, with coefficients derived from logistic regression performed on a training set, and DNA-binding is predicted based on whether the result is above a certain threshold. We show that the method is insensitive to errors in the atomic coordinates and provides correct predictions even on inaccurate protein models. We demonstrate that the method is capable of predicting proteins with novel binding site motifs and structures solved in an unbound state. The accuracy of our method is close to another, published method that uses all-atom structures, time-consuming calculations and information on conserved residues.     

A conserved nucleotide sequence, coding for a segment of the C-propeptide, is found at the same location in different collagen genes.

The nucleotide sequence of a segment of the chick alpha 1 type III collagen gene which codes for the C-propeptide was determined and compared with the corresponding sequence in the alpha 1 type I and alpha 2 type I collagen genes. As in the alpha 2 type I gene the coding information for the C-propeptide of the type III collagen gene is subdivided in four exons. Similarly, the amino proximal exon contains sequences for both the carboxy terminal end of the alpha-helical segment of collagen and for the beginning of the C-propeptide in both genes. Therefore, this organization of exons must have been established before these two collagen genes arose by duplication of a common ancestor. In several subsegments the deduced amino acid sequence for the C-propeptide of type III collagen shows a strong homology with the corresponding amino acid sequence in alpha 1 and alpha 2 type I. For one of these homologous amino acid sequences, however, the nucleotide sequence is much better conserved than for the others. It is possible that a mechanism of gene conversion has maintained the homogeneity of this nucleotide sequence among the interstitial collagen genes. Alternatively, the conserved nucleotide sequence may represent a regulatory signal which could function either in the DNA or in the RNA.     

Comparative analysis of structural properties of the C-type-lectin-like domain (CTLD)

The superfamily of proteins containing the C-type-lectin-like domain (CTLD) is a group of abundant extracellular metazoan proteins characterized by evolutionary flexibility and functional versatility. Several CTLDs are also found in parasitic prokaryotes and viruses. The 37 distinct currently available CTLD structures demonstrate significant structural conservation despite low or undetectable sequence similarity. Our aim in this study was to perform an extensive comparative analysis of all available CTLD structures to establish the most conserved structural features of the fold, and to test and extend the early analysis of Drickamer. By implication, these features should be those critical for maintenance of integrity of the fold. By analyzing CTLD structures superimposed by several methods, we have established groups of conserved structural positions involved in fold maintenance but not in ligand binding; these are consistent with the fold's known functional flexibility. In addition to the well-recognized disulfide bridges, groups of conserved residues are involved in hydrophobic interactions stabilizing the core of the fold and the long loop region, and in an alpha2-beta1-beta5 polar interaction. Evaluation of the conclusions of the structure comparison study compared with alignments of all available human, mouse and Caenorhabditis elegans CTLD sequences showed that conservation patterns are preserved throughout the whole CTLD sequence space. Our observations provide an improved understanding of CTLD structure, and will help in identification of new CTLDs and the mechanisms that drive and constrain the coevolution of the structure and function of the fold.     

Rheostats and Toggle Switches for Modulating Protein Function

The millions of protein sequences generated by genomics are expected to transform protein engineering and personalized medicine. To achieve these goals, tools for predicting outcomes of amino acid changes must be improved. Currently, advances are hampered by insufficient experimental data about nonconserved amino acid positions. Since the property “nonconserved” is identified using a sequence alignment, we designed experiments to recapitulate that context: Mutagenesis and functional characterization was carried out in 15 LacI/GalR homologs (rows) at 12 nonconserved positions (columns). Multiple substitutions were made at each position, to reveal how various amino acids of a nonconserved column were tolerated in each protein row. Results showed that amino acid preferences of nonconserved positions were highly context-dependent, had few correlations with physico-chemical similarities, and were not predictable from their occurrence in natural LacI/GalR sequences. Further, unlike the “toggle switch” behaviors of conserved positions, substitutions at nonconserved positions could be rank-ordered to show a “rheostatic”, progressive effect on function that spanned several orders of magnitude. Comparisons to various sequence analyses suggested that conserved and strongly co-evolving positions act as functional toggles, whereas other important, nonconserved positions serve as rheostats for modifying protein function. Both the presence of rheostat positions and the sequence analysis strategy appear to be generalizable to other protein families and should be considered when engineering protein modifications or predicting the impact of protein polymorphisms.     

Conserved key amino acid positions (CKAAPs) derived from the analysis of common substructures in proteins

An all-against-all protein structure comparison using the Combinatorial Extension (CE) algorithm applied to a representative set of PDB structures revealed a gallery of common substructures in proteins (http://cl.sdsc.edu/ce.html). These substructures represent commonly identified folds, domains, or components thereof. Most of the subsequences forming these similar substructures have no significant sequence similarity. We present a method to identify conserved amino acid positions and residue-dependent property clusters within these subsequences starting with structure alignments. Each of the subsequences is aligned to its homologues in SWALL, a nonredundant protein sequence database. The most similar sequences are purged into a common frequency matrix, and weighted homologues of each one of the subsequences are used in scoring for conserved key amino acid positions (CKAAPs). We have set the top 20% of the high-scoring positions in each substructure to be CKAAPs. It is hypothesized that CKAAPs may be responsible for the common folding patterns in either a local or global view of the protein-folding pathway. Where a significant number of structures exist, CKAAPs have also been identified in structure alignments of complete polypeptide chains from the same protein family or superfamily. Evidence to support the presence of CKAAPs comes from other computational approaches and experimental studies of mutation and protein-folding experiments, notably the Paracelsus challenge. Finally, the structural environment of CKAAPs versus non-CKAAPs is examined for solvent accessibility, hydrogen bonding, and secondary structure. The identification of CKAAPs has important implications for protein engineering, fold recognition, modeling, and structure prediction studies and is dependent on the availability of structures and an accurate structure alignment methodology. Proteins 2001;42:148-163.