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

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

Crystal structure of Medicago truncatula UGT85H2--insights into the structural basis of a multifunctional (iso)flavonoid glycosyltransferase

(Iso)flavonoids are a diverse group of plant secondary metabolites with important effects on plant, animal and human health. They exist in various glycosidic forms. Glycosylation, which may determine their bioactivities and functions, is controlled by specific plant uridine diphosphate glycosyltransferases (UGTs). We describe a new multifunctional (iso)flavonoid glycosyltransferase, UGT85H2, from the model legume Medicago truncatula with activity towards a number of phenylpropanoid-derived natural products including the flavonol kaempferol, the isoflavone biochanin A, and the chalcone isoliquiritigenin. The crystal structure of UGT85H2 has been determined at 2.1 A resolution, and reveals distinct structural features that are different from those of other UGTs and related to the enzyme's functions and substrate specificities. Structural and comparative analyses revealed the putative binding sites for the donor and acceptor substrates that are located in a large cleft formed between the two domains of the enzyme, and indicated that Trp360 may undergo a conformational change after sugar donor binding to the enzyme. UGT85H2 has higher specificity for flavonol than for isoflavone. Further substrate docking combined with enzyme activity assay and kinetic analysis provided structural insights into this substrate specificity and preference.     

Crystal structure of the cofactor-binding domain of the human phase II drug-metabolism enzyme UDP-glucuronosyltransferase 2B7

Human UDP-glucuronosyltransferases (UGT) are the dominant phase II conjugative drug metabolism enzymes that also play a central role in processing a range of endobiotic compounds. UGTs catalyze the covalent addition of glucuronic acid sugar moieties to a host of therapeutics and environmental toxins, as well as to a variety of endogenous steroids and other signaling molecules. We report the 1.8-A resolution apo crystal structure of the UDP-glucuronic acid binding domain of human UGT isoform 2B7 (UGT2B7), which catalyzes the conjugative elimination of opioid, antiviral, and anticancer drugs. This is the first crystal structure of any region of a mammalian UGT drug metabolism enzyme. Designated UGT2B7 mutants at residues predicted to interact with the UDP-glucuronic acid cofactor exhibited significantly impaired catalytic activity, with maximum effects observed for amino acids closest to the glucuronic acid sugar transferred to the acceptor molecule. Homology modeling of UGT2B7 with related plant flavonoid glucosyltransferases indicates human UGTs share a common catalytic mechanism. Point mutations at predicted catalytic residues in UGT2B7 abrogated activity, strongly suggesting human UGTs also utilize a serine hydrolase-like catalytic mechanism to facilitate glucuronic acid transfer.     

植物源UDP-糖基转移酶及其分子改造*

糖基化能够增加化合物的结构多样性,有效改善水溶性、药理活性和生物利用度,对植物天然产物的药物开发至关重要。UDP-糖基转移酶(UGTs)能够催化糖基从活化的核苷酸糖供体转移到受体形成糖苷键,植物中天然产物的糖基化修饰主要通过UGTs实现。但是大多数天然UGTs的催化活性、稳定性和底物特异性较低,难以满足工业用酶的要求,限制了它们的工业化应用。近年来,通过分子改造技术改进天然UGTs的催化特性取得了突破性的研究进展。为此,概述了植物源UGTs的挖掘与表征、三维结构和催化机制,归纳了UGTs分子改造的思路和方法,包括理性设计和定向进化,重点介绍了结构域替换、序列保守分析和结构分析与定点突变的结合,总结了定向进化中的高通量筛选方法,为植物天然产物酶法糖基化的工业应用提供了参考和借鉴。     

Completion of the canonical pathway for assembly of anticancer drugs vincristine/vinblastine in Catharanthus roseus

Glycosylation is a key modification for most molecules including plant natural products, for example, flavonoids and isoflavonoids, and can enhance the bioactivity and bioavailability of the natural products. The crystal structure of plant rhamnosyltransferase UGT89C1 from Arabidopsis thaliana was determined, and the structures of UGT89C1 in complexes with UDP‐β‐l ‐rhamnose and acceptor quercetin revealed the detailed interactions between the enzyme and its substrates. Structural and mutational analysis indicated that Asp356, His357, Pro147 and Ile148 are key residues for sugar donor recognition and specificity for UDP‐β‐l ‐rhamnose. The mutant H357Q exhibited activity with both UDP‐β‐l ‐rhamnose and UDP‐glucose. Structural comparison and mutagenesis confirmed that His21 is a key residue as the catalytic base and the only catalytic residue involved in catalysis independently as UGT89C1 lacks the other catalytic Asp that is highly conserved in other reported UGTs and forms a hydrogen bond with the catalytic base His. Ser124 is located in the corresponding position of the catalytic Asp in other UGTs and is not able to form a hydrogen bond with His21. Mutagenesis further showed that Ser124 may not be important in its catalysis, suggesting that His21 and acceptor may form an acceptor‐His dyad and UGT89C1 utilizes a catalytic dyad in catalysis instead of catalytic triad. The information of structure and mutagenesis provides structural insights into rhamnosyltransferase substrate specificity and rhamnosylation mechanism.     

Conservation of structure and function among tyrosine recombinases: homology-based modeling of the lambda integrase core-binding domain

Tyrosine recombinases participate in diverse biological processes by catalyzing recombination between specific DNA sites. Although a conserved protein fold has been described for the catalytic (CAT) domains of five recombinases, structural relationships between their core-binding (CB) domains remain unclear. Despite differences in the specificity and affinity of core-type DNA recognition, a conserved binding mechanism is suggested by the shared two-domain motif in crystal structure models of the recombinases Cre, XerD and Flp. We have found additional evidence for conservation of the CB domain fold. Comparison of XerD and Cre crystal structures showed that their CB domains are closely related; the three central α-helices of these domains are superposable to within 1.44 Å. A structure-based multiple sequence alignment containing 25 diverse CB domain sequences provided evidence for widespread conservation of both structural and functional elements in this fold. Based upon the Cre and XerD crystal structures, we employed homology modeling to construct a three-dimensional structure for the λ integrase CB domain. The model provides a conceptual framework within which many previously identified, functionally important amino acid residues were investigated. In addition, the model predicts new residues that may participate in core-type DNA binding or dimerization, thereby providing hypotheses for future genetic and biochemical experiments.     

Similarities and differences among 105 members of the Int family of site-specific recombinases.

Alignments of 105 site-specific recombinases belonging to the Int family of proteins identified extended areas of similarity and three types of structural differences. In addition to the previously recognized conservation of the tetrad R-H-R-Y, located in boxes I and II, several newly identified sequence patches include charged amino acids that are highly conserved and a specific pattern of buried residues contributing to the overall protein fold. With some notable exceptions, unconserved regions correspond to loops in the crystal structures of the catalytic domains of lambda Int (Int c170) and HP1 Int (HPC) and of the recombinases XerD and Cre. Two structured regions also harbor some pronounced differences. The first comprises beta-sheets 4 and 5, alpha-helix D and the adjacent loop connecting it to alpha-helix E: two Ints of phages infecting thermophilic bacteria are missing this region altogether; the crystal structures of HPC, XerD and Cre reveal a lack of beta-sheets 4 and 5; Cre displays two additional beta-sheets following alpha-helix D; five recombinases carry large insertions. The second involves the catalytic tyrosine and is seen in a comparison of the four crystal structures. The yeast recombinases can theoretically be fitted to the Int fold, but the overall differences, involving changes in spacing as well as in motif structure, are more substantial than seen in most other proteins. The phenotypes of mutations compiled from several proteins are correlated with the available structural information and structure-function relationships are discussed. In addition, a few prokaryotic and eukaryotic enzymes with partial homology with the Int family of recombinases may be distantly related, either through divergent or convergent evolution. These include a restriction enzyme and a subgroup of eukaryotic RNA helicases (D-E-A-D proteins).     

Local Differentiation of Sugar Donor Specificity of Flavonoid Glycosyltransferase in Lamiales

Flavonoids are most commonly conjugated with various sugar moieties by UDP-sugar:glycosyltransferases (UGTs) in a lineage-specific manner. Generally, the phylogenetics and regiospecificity of flavonoid UGTs are correlated, indicating that the regiospecificity of UGT differentiated prior to speciation. By contrast, it is unclear how the sugar donor specificity of UGTs evolved. Here, we report the biochemical, homology-modeled, and phylogenetic characterization of flavonoid 7-O-glucuronosyltransferases (F7GAT), which is responsible for producing specialized metabolites in Lamiales plants. All of the Lamiales F7GATs were found to be members of the UGT88-related cluster and specifically used UDP-glucuronic acid (UDPGA). We identified an Arg residue that is specifically conserved in the PSPG box in the Lamiales F7GATs. Substitution of this Arg with Trp was sufficient to convert the sugar donor specificity of the Lamiales F7GATs from UDPGA to UDP-glucose. Homology modeling of the Lamiales F7GAT suggested that the Arg residue plays a critical role in the specific recognition of anionic carboxylate of the glucuronic acid moiety of UDPGA with its cationic guanidinium moiety. These results support the hypothesis that differentiation of sugar donor specificity of UGTs occurred locally, in specific plant lineages, after establishment of general regiospecificity for the sugar acceptor. Thus, the plasticity of sugar donor specificity explains, in part, the extraordinary structural diversification of phytochemicals.     

Structure, mechanism and engineering of plant natural product glycosyltransferases

Glycosylation is a key mechanism in determining chemical complexity and diversity of plant natural products, and influencing their chemical properties and bioactivities. Uridine diphosphate glycosyltransferases (UGTs) are the central players in these glycosylation processes for decorating natural products with sugars. Crystal structures of plant UGTs have revealed their exquisite architectures and provided the structural basis for understanding their catalytic mechanism and substrate specificity. Structure-based UGT engineering can alter substrate specificity; compromise or enhance catalytic efficiency; and confer reversibility to the glycosylation reaction. This review highlights the structural insights on plant UGTs and successes in glycosylation engineering.     

Structural relationships of homologous proteins as a fundamental principle in homology modeling

Protein structure prediction is based mainly on the modeling of proteins by homology to known structures; this knowledgebased approach is the most promising method to date. Although it is used in the whole area of protein research, no general rules concerning the quality and applicability of concepts and procedures used in homology modeling have been put forward yet. Therefore, the main goal of the present work is to provide tools for the assessment of accuracy of modeling at a given level of sequence homology. A large set of known structures from different conformational and functional classes, but various degrees of homology was selected. Pairwise structure superpositions were performed. Starting with the definition of the structurally conserved regions and determination of topologically correct sequence alignments, we correlated geometrical properties with sequence homology (defined by the 250 PAM Dayhoff Matrix) and identity. It is shown that both the topological differences of the protein backbones and the relative positions of corresponding side chains diverge with decreasing sequence identity. Below 50% identity, the deviation in regions that are structurally not conserved continually increases, thus implying that with decreasing sequence identity modeling has to take into account more and more structurally diverging loop regions that are difficult to predict. © 1993 Wiley-Liss, Inc.     

Enhanced activity of <Emphasis Type="Italic">Withania somnifera</Emphasis> family-1 glycosyltransferase (UGT73A16) via mutagenesis

This work used an approach of enzyme engineering towards the improved production of baicalin as well as alteration of acceptor and donor substrate preferences in UGT73A16. The 3D model of Withania somnifera family-1 glycosyltransferase (UGT73A16) was constructed based on the known crystal structures of plant UGTs. Structural and functional properties of UGT73A16 were investigated using docking and mutagenesis. The docking studies were performed to understand the key residues involved in substrate recognition. In the molecular model of UGT73A16, substrates binding pockets are located between N- and C-terminal domains. Modeled UGT73A16 was docked with UDP-glucose, UDP-glucuronic acid (UDPGA), kaempferol, isorhamnetin, 3-hydroxy flavones, naringenin, genistein and baicalein. The protein–ligand interactions showed that His 16, Asp 246, Lys 255, Ala 337, Gln 339, Val 340, Asn 358 and Glu 362 amino acid residues may be important for catalytic activity. The kinetic parameters indicated that mutants A337C and Q339A exhibited 2–3 fold and 6–7 fold more catalytic efficiency, respectively than wild type, and shifted the sugar donor specificity from UDP-glucose to UDPGA. The mutant Q379H displayed large loss of activity with UDP-glucose and UDPGA strongly suggested that last amino acid residue of PSPG box is important for glucuronosylation and glucosylation and highly specific to sugar binding sites. The information obtained from docking and mutational studies could be beneficial in future to engineer this biocatalyst for development of better ones.     

Sequence–structure relationships in RNA loops: establishing the basis for loop homology modeling

The specific function of RNA molecules frequently resides in their seemingly unstructured loop regions. We performed a systematic analysis of RNA loops extracted from experimentally determined three-dimensional structures of RNA molecules. A comprehensive loop-structure data set was created and organized into distinct clusters based on structural and sequence similarity. We detected clear evidence of the hallmark of homology present in the sequence–structure relationships in loops. Loops differing by <25% in sequence identity fold into very similar structures. Thus, our results support the application of homology modeling for RNA loop model building. We established a threshold that may guide the sequence divergence-based selection of template structures for RNA loop homology modeling. Of all possible sequences that are, under the assumption of isosteric relationships, theoretically compatible with actual sequences observed in RNA structures, only a small fraction is contained in the Rfam database of RNA sequences and classes implying that the actual RNA loop space may consist of a limited number of unique loop structures and conserved sequences. The loop-structure data sets are made available via an online database, RLooM. RLooM also offers functionalities for the modeling of RNA loop structures in support of RNA engineering and design efforts.     

Determination of catalytic key amino acids and UDP sugar donor specificity of the cyanohydrin glycosyltransferase UGT85B1 from Sorghum bicolor. Molecular modeling substantiated by site-specific mutagenesis and biochemical analyses

Plants produce a plethora of structurally diverse natural products. The final step in their biosynthesis is often a glycosylation step catalyzed by a family 1 glycosyltransferase (GT). In biosynthesis of the cyanogenic glucoside dhurrin in Sorghum bicolor, the UDP-glucosyltransferase UGT85B1 catalyzes the conversion of p-hydroxymandelonitrile into dhurrin. A structural model of UGT85B1 was built based on hydrophobic cluster analysis and the crystal structures of two bacterial GTs, GtfA and GtfB, which each showed approximately 15% overall amino acid sequence identity to UGT85B1. The model enabled predictions about amino acid residues important for catalysis and sugar donor specificity. p-Hydroxymandelonitrile and UDP-glucose (Glc) were predicted to be positioned within hydrogen-bonding distance to a glutamic acid residue in position 410 facilitating sugar transfer. The acceptor was packed within van der Waals distance to histidine H23. Serine S391 and arginine R201 form hydrogen bonds to the pyrophosphate part of UDP-Glc and hence stabilize binding of the sugar donor. Docking of UDP sugars predicted that UDP-Glc would serve as the sole donor sugar in UGT85B1. This was substantiated by biochemical analyses. The predictive power of the model was validated by site-directed mutagenesis of selected residues and using enzyme assays. The modeling approach has provided a tool to design GTs with new desired substrate specificities for use in biotechnological applications. The modeling identified a hypervariable loop (amino acid residues 156-188) that contained a hydrophobic patch. The involvement of this loop in mediating binding of UGT85B1 to cytochromes P450, CYP79A1, and CYP71E1 within a dhurrin metabolon is discussed.     

Structural characterization of a flavonoid glycosyltransferase from Withania somnifera

Medicinal plants are extensively utilized in traditional and herbal medicines, both in India and around the world due to the presence of diverse low molecular weight natural products such as flavonoids, alkaloids, terpenoids and sterols. Flavonoids which have health benefits for humans are the large class of phenylpropanoid-derived secondary metabolites and are mostly glycosylated by UDP-glycosyltransferases (UGTs). Although large numbers of different UGTs are known from higher plants, very few protein structures have been reported till now. In the present study, the three-dimensional model of flavonoid specific glycosyltransferases (WsFGT) from Withania somnifera was constructed based on the crystal structure of plant UGTs. The resulted model was assessed by various tools and the final refined model revealed GT-B type fold. Further, to understand the sugar donors and acceptors interactions with the active site of WsFGT, docking studies were performed. The amino acids from conserved PSPG box were interacted with sugar donor while His18, Asp110, Trp352 and Asn353 were important for catalytic function. This structural and docking information will be useful to understand the glycosylation mechanism of flavonoid glucosides.

Abbreviations

DOPE - Discrete Optimized Potential Energy, PDB - Protein Data Bank, PSPG - Plant Secondary Product Glycosyltransferase, RMSD - Root Mean Squared Deviation, UDP - Uridine diphosphate, UGT - UDP-glycosyltransferases.     

Conservation of polyproline II helices in homologous proteins: implications for structure prediction by model building.

Left-handed polyproline II (PPII) helices commonly occur in globular proteins in segments of 4-8 residues. This paper analyzes the structural conservation of PPII-helices in 3 protein families: serine proteinases, aspartic proteinases, and immunoglobulin constant domains. Calculations of the number of conserved segments based on structural alignment of homologous molecules yielded similar results for the PPII-helices, the alpha-helices, and the beta-strands. The PPII-helices are consistently conserved at the level of 100-80% in the proteins with sequence identity above 20% and RMS deviation of structure alignments below 3.0 A. The most structurally important PPII segments are conserved below this level of sequence identity. These results suggest that the PPII-helices, in addition to the other 2 secondary structure classes, should be identified as part of structurally conserved regions in proteins. This is supported by similar values for the local RMS deviations of the aligned segments for the structural classes of PPII-helices, alpha-helices, and beta-strands. The PPII-helices are shown to participate in supersecondary elements such as PPII-helix/alpha-helix. The conservation of PPII-helices depends on the conservation of a supersecondary element as a whole. PPII-helices also form links, possibly flexible, in the interdomain regions. The role of the PPII-helices in model building by homology is 2-fold; they serve as additional conserved elements in the structure allowing improvement of the accuracy of a model and provide correct chain geometry for modeling of the segments equivalenced to them in a target sequence. The improvement in model building is demonstrated in 2 test studies.     

Phylogenetic analysis of the UDP-glycosyltransferase multigene family of Arabidopsis thaliana

A class of UDP-glycosyltransferases (UGTs) defined by the presence of a C-terminal consensus sequence is found throughout the plant and animal kingdoms. Whereas mammalian enzymes use UDP-glucuronic acid, the plant enzymes typically use UDP-glucose in the transfer reactions. A diverse array of aglycones can be glucosylated by these UGTs. In plants, the aglycones include plant hormones, secondary metabolites involved in stress and defense responses, and xenobiotics such as herbicides. Glycosylation is known to regulate many properties of the aglycones such as their bioactivity, their solubility, and their transport properties within the cell and throughout the plant. As a means of providing a framework to start to understand the substrate specificities and structure-function relationships of plant UGTs, we have now applied a molecular phylogenetic analysis to the multigene family of 99 UGT sequences in Arabidopsis. We have determined the overall organization and evolutionary relationships among individual members with a surprisingly high degree of confidence. Through constructing a composite phylogenetic tree that also includes all of the additional plant UGTs with known catalytic activities, we can start to predict both the evolutionary history and substrate specificities of new sequences as they are identified. The tree already suggests that while the activities of some subgroups of the UGT family are highly conserved among different plant species, others subgroups shift substrate specificity with relative ease.     

Crystal structures of two FGF-FGFR complexes reveal the determinants of ligand-receptor specificity

To elucidate the structural determinants governing specificity in fibroblast growth factor (FGF) signaling, we have determined the crystal structures of FGF1 and FGF2 complexed with the ligand binding domains (immunoglobulin-like domains 2 [D2] and 3 [D3]) of FGF receptor 1 (FGFR1) and FGFR2, respectively. Highly conserved FGF-D2 and FGF-linker (between D2-D3) interfaces define a general binding site for all FGF-FGFR complexes. Specificity is achieved through interactions between the N-terminal and central regions of FGFs and two loop regions in D3 that are subject to alternative splicing. These structures provide a molecular basis for FGF1 as a universal FGFR ligand and for modulation of FGF-FGFR specificity through primary sequence variations and alternative splicing.     

Structural conservation in single-domain proteins: implications for homology modeling

Large-scale sequencing projects are widening the gap between the known protein universe and the fraction for which structural information has been experimentally obtained. Through the application of homology (comparative) modeling and more general structure prediction techniques, this gap can, however, be narrowed, providing indirect structural information for a considerable number of proteins. Moreover, the estimated number of existing protein folds seems to be limited and many of these yet unknown folds should be discovered by dedicated large-scale structural genomics projects. Within this perspective, homology (comparative) modeling will gain in importance, as will the use of models derived by this technique. Here we discuss how well a sequence alignment, the most common starting point for generating a model, reflects the structural conservation between homologous proteins and we show that sequence information is able to direct construction of acceptable models as far as the structural core is concerned. We also show here that the regions surrounding insertions and deletions are much less conserved than the core and discuss the implications of this observation for loop modeling.     

Homology model of surface antigen OmpC from Salmonella typhi and its functional implications

Homology based 3D structural model of the immunodominant major surface antigen OmpC from Salmonella typhi, an obligatory human pathogen, was built to understand the possible unique conformational features of its antigenic loops with respect to other immunologically cross reacting porins. The homology model was built based on the known crystal structures of the E. coli porins OmpF and PhoE. Structure based sequence alignment helped to define the structurally conserved regions (SCRs). The SCR regions of OmpC were modelled using the coordinates of corresponding regions from reference proteins. Surface exposed variable regions were modelled based on the sequence similarity and loop search in PDB. Structural refinement based on symmetry restrained energy minimization resulted in an agreeable model for the trimer of OmpC. The resulting model was compared with other porin structures, having b-barrel fold with 16 transmembrane beta-strands, and found that the variable regions are unique in terms of sequence and structure. A ranking of the loops taking into account the antigenic index, the sequence variability, the surface accessibility in the context of the trimer, and the structural variability suggests that loop 4 (151-172), loop 5 (194-218) and loop 6 (237-264) are the best ranked B-cell epitopes. The model provides possible explanations for the functional and unique immunological properties associated with the surface exposed regions and outlines the implications for structure based experimental design.