β-Propeller Blades as Ancestral Peptides in Protein Evolution

Proteins of the β-propeller fold are ubiquitous in nature and widely used as structural scaffolds for ligand binding and enzymatic activity. This fold comprises between four and twelve four-stranded β-meanders, the so called blades that are arranged circularly around a central funnel-shaped pore. Despite the large size range of β-propellers, their blades frequently show sequence similarity indicative of a common ancestry and it has been proposed that the majority of β-propellers arose divergently by amplification and diversification of an ancestral blade. Given the structural versatility of β-propellers and the hypothesis that the first folded proteins evolved from a simpler set of peptides, we investigated whether this blade may have given rise to other folds as well. Using sequence comparisons, we identified proteins of four other folds as potential homologs of β-propellers: the luminal domain of inositol-requiring enzyme 1 (IRE1-LD), type II β-prisms, β-pinwheels, and WW domains. Because, with increasing evolutionary distance and decreasing sequence length, the statistical significance of sequence comparisons becomes progressively harder to distinguish from the background of convergent similarities, we complemented our analyses with a new method that evaluates possible homology based on the correlation between sequence and structure similarity. Our results indicate a homologous relationship of IRE1-LD and type II β-prisms with β-propellers, and an analogous one for β-pinwheels and WW domains. Whereas IRE1-LD most likely originated by fold-changing mutations from a fully formed PQQ motif β-propeller, type II β-prisms originated by amplification and differentiation of a single blade, possibly also of the PQQ type. We conclude that both β-propellers and type II β-prisms arose by independent amplification of a blade-sized fragment, which represents a remnant of an ancient peptide world.     

Hemolytic Lectin CEL-III Heptamerizes via a Large Structural Transition from α-Helices to a β-Barrel during the Transmembrane Pore Formation Process

CEL-III is a hemolytic lectin isolated from the sea cucumber Cucumaria echinata. This lectin is composed of two carbohydrate-binding domains (domains 1 and 2) and one oligomerization domain (domain 3). After binding to the cell surface carbohydrate chains through domains 1 and 2, domain 3 self-associates to form transmembrane pores, leading to cell lysis or death, which resembles other pore-forming toxins of diverse organisms. To elucidate the pore formation mechanism of CEL-III, the crystal structure of the CEL-III oligomer was determined. The CEL-III oligomer has a heptameric structure with a long β-barrel as a transmembrane pore. This β-barrel is composed of 14 β-strands resulting from a large structural transition of α-helices accommodated in the interface between domains 1 and 2 and domain 3 in the monomeric structure, suggesting that the dissociation of these α-helices triggered their structural transition into a β-barrel. After heptamerization, domains 1 and 2 form a flat ring, in which all carbohydrate-binding sites remain bound to cell surface carbohydrate chains, stabilizing the transmembrane β-barrel in a position perpendicular to the plane of the lipid bilayer.     

Atomic Interaction Networks in the Core of Protein Domains and Their Native Folds

Vastly divergent sequences populate a majority of protein folds. In the quest to identify features that are conserved within protein domains belonging to the same fold, we set out to examine the entire protein universe on a fold-by-fold basis. We report that the atomic interaction network in the solvent-unexposed core of protein domains are fold-conserved, extraordinary sequence divergence notwithstanding. Further, we find that this feature, termed protein core atomic interaction network (or PCAIN) is significantly distinguishable across different folds, thus appearing to be “signature” of a domain''s native fold. As part of this study, we computed the PCAINs for 8698 representative protein domains from families across the 1018 known protein folds to construct our seed database and an automated framework was developed for PCAIN-based characterization of the protein fold universe. A test set of randomly selected domains that are not in the seed database was classified with over 97% accuracy, independent of sequence divergence. As an application of this novel fold signature, a PCAIN-based scoring scheme was developed for comparative (homology-based) structure prediction, with 1–2 angstroms (mean 1.61A) C_α RMSD generally observed between computed structures and reference crystal structures. Our results are consistent across the full spectrum of test domains including those from recent CASP experiments and most notably in the ‘twilight’ and ‘midnight’ zones wherein <30% and <10% target-template sequence identity prevails (mean twilight RMSD of 1.69A). We further demonstrate the utility of the PCAIN protocol to derive biological insight into protein structure-function relationships, by modeling the structure of the YopM effector novel E3 ligase (NEL) domain from plague-causative bacterium Yersinia Pestis and discussing its implications for host adaptive and innate immune modulation by the pathogen. Considering the several high-throughput, sequence-identity-independent applications demonstrated in this work, we suggest that the PCAIN is a fundamental fold feature that could be a valuable addition to the arsenal of protein modeling and analysis tools.     

Surface Tensiometry of Apolipoprotein B Domains at Lipid Interfaces Suggests a New Model for the Initial Steps in Triglyceride-rich Lipoprotein Assembly

Apolipoprotein B (apoB) is the principal protein component of triacylglyceride (TAG)-rich lipoproteins, including chylomicrons and very low density lipoprotein, which is the precursor to LDL (the “bad cholesterol”). TAG-rich lipoprotein assembly is initiated by the N-terminal βα1 superdomain of apoB, which co-translationally binds and remodels the luminal leaflet of the rough endoplasmic reticulum. The βα1 superdomain contains four domains and is predicted to interact directly with lipids. Using drop tensiometry, we examined the interfacial properties of the α-helical and C-sheet domains and several subdomains to establish a detailed structure-function relationship at the lipid/water interface. The adsorption, stress response, exchangeability, and pressure (Π)-area relationship were studied at both triolein/water and triolein/1-palmitoyl, 2-oleoylphosphatidylcholine/water interfaces that mimic physiological environments. The α-helical domain spontaneously adsorbed to a triolein/water interface and formed a viscoelastic surface. It was anchored to the surface by helix 6, and the other helices were ejected and/or remodeled on the surface as a function of surface pressure. The C-sheet instead formed an elastic film on a triolein/water interface and was irreversibly anchored to the lipid surface, which is consistent with the behavior of amphipathic β-strands. When both domains were adsorbed together on the surface, the C-sheet shielded a portion of the α-helical domain from the surface, which retained its globular structure. Overall, the unique secondary and tertiary structures of the N-terminal domains of apoB support the intrinsic capability of co-translational lipid recruitment. The evidence presented here allows the construction of a detailed model of the initiation of TAG-rich lipoprotein assembly.     

Resurrection of an ancient inflammatory locus reveals switch to caspase-1 specificity on a caspase-4 scaffold

Pyroptosis is a mechanism of inflammatory cell death mediated by the activation of the prolytic protein gasdermin D by caspase-1, caspase-4, and caspase-5 in human, and caspase-1 and caspase-11 in mouse. In addition, caspase-1 amplifies inflammation by proteolytic activation of cytokine interleukin-1β (IL-1β). Modern mammals of the order Carnivora lack the caspase-1 catalytic domain but express an unusual version of caspase-4 that can activate both gasdermin D and IL-1β. Seeking to understand the evolutionary origin of this caspase, we utilized the large amount of data available in public databases to perform ancestral sequence reconstruction of an inflammatory caspase of a Carnivora ancestor. We expressed the catalytic domain of this putative ancestor in Escherichia coli, purified it, and compared its substrate specificity on synthetic and protein substrates to extant caspases. We demonstrated that it activates gasdermin D but has reduced ability to activate IL-1β. Our reconstruction suggests that caspase-1 was lost in a Carnivora ancestor, perhaps upon a selective pressure for which the generation of biologically active IL-1β by caspase-1 was detrimental. We speculate that later, a Carnivora encountered selective pressures that required the production of IL-1β, and caspase-4 subsequently gained this activity. This hypothesis would explain why extant Carnivora possess an inflammatory caspase with caspase-1 catalytic function placed on a caspase-4 scaffold.     

Insight into the structure of amyloid fibrils from the analysis of globular proteins

The conversion from soluble states into cross-β fibrillar aggregates is a property shared by many different proteins and peptides and was hence conjectured to be a generic feature of polypeptide chains. Increasing evidence is now accumulating that such fibrillar assemblies are generally characterized by a parallel in-register alignment of β-strands contributed by distinct protein molecules. Here we assume a universal mechanism is responsible for β-structure formation and deduce sequence-specific interaction energies between pairs of protein fragments from a statistical analysis of the native folds of globular proteins. The derived fragment–fragment interaction was implemented within a novel algorithm, prediction of amyloid structure aggregation (PASTA), to investigate the role of sequence heterogeneity in driving specific aggregation into ordered self-propagating cross-β structures. The algorithm predicts that the parallel in-register arrangement of sequence portions that participate in the fibril cross-β core is favoured in most cases. However, the antiparallel arrangement is correctly discriminated when present in fibrils formed by short peptides. The predictions of the most aggregation-prone portions of initially unfolded polypeptide chains are also in excellent agreement with available experimental observations. These results corroborate the recent hypothesis that the amyloid structure is stabilised by the same physicochemical determinants as those operating in folded proteins. They also suggest that side chain–side chain interaction across neighbouring β-strands is a key determinant of amyloid fibril formation and of their self-propagating ability.     

Solution Structure of Tensin2 SH2 Domain and Its Phosphotyrosine-Independent Interaction with DLC-1

Background

Src homology 2 (SH2) domain is a conserved module involved in various biological processes. Tensin family member was reported to be involved in tumor suppression by interacting with DLC-1 (deleted-in-liver-cancer-1) via its SH2 domain. We explore here the important questions that what the structure of tensin2 SH2 domain is, and how it binds to DLC-1, which might reveal a novel binding mode.

Principal Findings

Tensin2 SH2 domain adopts a conserved SH2 fold that mainly consists of five β-strands flanked by two α-helices. Most SH2 domains recognize phosphorylated ligands specifically. However, tensin2 SH2 domain was identified to interact with nonphosphorylated ligand (DLC-1) as well as phosphorylated ligand.

Conclusions

We determined the solution structure of tensin2 SH2 domain using NMR spectroscopy, and revealed the interactions between tensin2 SH2 domain and its ligands in a phosphotyrosine-independent manner.     

SH3 Domains Differentially Stimulate Distinct Dynamin I Assembly Modes and G Domain Activity

Dynamin I is a highly regulated GTPase enzyme enriched in nerve terminals which mediates vesicle fission during synaptic vesicle endocytosis. One regulatory mechanism involves its interactions with proteins containing Src homology 3 (SH3) domains. At least 30 SH3 domain-containing proteins bind dynamin at its proline-rich domain (PRD). Those that stimulate dynamin activity act by promoting its oligomerisation. We undertook a systematic parallel screening of 13 glutathione-S-transferase (GST)-tagged endocytosis-related SH3 domains on dynamin binding, GTPase activity and oligomerisation. No correlation was found between dynamin binding and their potency to stimulate GTPase activity. There was limited correlation between the extent of their ability to stimulate dynamin activity and the level of oligomerisation, indicating an as yet uncharacterised allosteric coupling of the PRD and G domain. We examined the two variants, dynamin Iab and Ibb, which differ in the alternately splice middle domain α2 helix. They responded differently to the panel of SH3s, with the extent of stimulation between the splice variants varying greatly between the SH3s. This study reveals that SH3 binding can act as a heterotropic allosteric regulator of the G domain via the middle domain α2 helix, suggesting an involvement of this helix in communicating the PRD-mediated allostery. This indicates that SH3 binding both stabilises multiple conformations of the tetrameric building block of dynamin, and promotes assembly of dynamin-SH3 complexes with distinct rates of GTP hydrolysis.     

Involvement of Phospholipase Cγ1 in Mouse Egg Activation Induced by a Truncated Form of the C-kit Tyrosine Kinase Present in Spermatozoa

Microinjection of a truncated form of the c-kit tyrosine kinase present in mouse spermatozoa (tr-kit) activates mouse eggs parthenogenetically, and tr-kit– induced egg activation is inhibited by preincubation with an inhibitor of phospholipase C (PLC) (Sette, C., A. Bevilacqua, A. Bianchini, F. Mangia, R. Geremia, and P. Rossi. 1997. Development [Camb.]. 124:2267–2274). Co-injection of glutathione-S-transferase (GST) fusion proteins containing the src-homology (SH) domains of the γ1 isoform of PLC (PLCγ1) competitively inhibits tr-kit– induced egg activation. A GST fusion protein containing the SH3 domain of PLCγ1 inhibits egg activation as efficiently as the whole SH region, while a GST fusion protein containing the two SH2 domains is much less effective. A GST fusion protein containing the SH3 domain of the Grb2 adaptor protein does not inhibit tr-kit–induced egg activation, showing that the effect of the SH3 domain of PLCγ1 is specific. Tr-kit–induced egg activation is also suppressed by co-injection of antibodies raised against the PLCγ1 SH domains, but not against the PLCγ1 COOH-terminal region. In transfected COS cells, coexpression of PLCγ1 and tr-kit increases diacylglycerol and inositol phosphate production, and the phosphotyrosine content of PLCγ1 with respect to cells expressing PLCγ1 alone. These data indicate that tr-kit activates PLCγ1, and that the SH3 domain of PLCγ1 is essential for tr-kit–induced egg activation.     

Conformational Transitions of the Cross-linking Domains of Elastin during Self-assembly

Elastin is the intrinsically disordered polymeric protein imparting the exceptional properties of extension and elastic recoil to the extracellular matrix of most vertebrates. The monomeric precursor of elastin, tropoelastin, as well as polypeptides containing smaller subsets of the tropoelastin sequence, can self-assemble through a colloidal phase separation process called coacervation. Present understanding suggests that self-assembly is promoted by association of hydrophobic domains contained within the tropoelastin sequence, whereas polymerization is achieved by covalent joining of lysine side chains within distinct alanine-rich, α-helical cross-linking domains. In this study, model elastin polypeptides were used to determine the structure of cross-linking domains during the assembly process and the effect of sequence alterations in these domains on assembly and structure. CD temperature melts indicated that partial α-helical structure in cross-linking domains at lower temperatures was absent at physiological temperature. Solid-state NMR demonstrated that β-strand structure of the cross-linking domains dominated in the coacervate state, although α-helix was predominant after subsequent cross-linking of lysine side chains with genipin. Mutation of lysine residues to hydrophobic amino acids, tyrosine or alanine, leads to increased propensity for β-structure and the formation of amyloid-like fibrils, characterized by thioflavin-T binding and transmission electron microscopy. These findings indicate that cross-linking domains are structurally labile during assembly, adapting to changes in their environment and aggregated state. Furthermore, the sequence of cross-linking domains has a dramatic effect on self-assembly properties of elastin-like polypeptides, and the presence of lysine residues in these domains may serve to prevent inappropriate ordered aggregation.     

Three-dimensional structure of ribulose-1,5-bisphosphate carboxylase/oxygenase from Rhodospirillum rubrum at 2.9 A resolution

The three-dimensional structure of ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) from Rhodospirillum rubrum has been determined at 2.9 Å resolution by X-ray crystallographic methods. The MIR-electron density map was substantially improved by two-fold non-crystallographic symmetry averaging. The polypeptide chains in the dimer were traced using a graphics display system with the help of the BONES option in FRODO. The dimer has approximate dimensions of 50 x 72 x 105 Å. The enzyme subunit is a typical two-domain protein. The smaller, N-terminal domain consists of 137 amino acid residues and forms a central, mixed five-stranded β-sheet with α-helices on both sides of the sheet. The larger C-terminal domain consists of 329 amino acid residues. This domain has an eight-stranded parallel α/β barrel structure as found in triosephosphate isomerase and a number of other functionally non-related proteins. The active site in Rubisco determined by difference Fourier techniques and fitting of active site residues to the electron density map, is located at the carboxy-end of the β-strands in the α/β barrel of the C-terminal domain. There are few domain–domain interactions within the subunit. The interactions at the interface between the two subunits of the dimer are tight and extensive. There are tight contacts between the two C-terminal domains, which build up the core of the molecule. There are also interactions between the N-terminal domain of one subunit and the C-terminal domain of the second subunit, close to the active site.     

Calcium Release at Fertilization in Starfish Eggs Is Mediated by Phospholipase Cγ

Although inositol trisphosphate (IP₃) functions in releasing Ca²⁺ in eggs at fertilization, it is not known how fertilization activates the phospholipase C that produces IP₃. To distinguish between a role for PLCγ, which is activated when its two src homology-2 (SH2) domains bind to an activated tyrosine kinase, and PLCβ, which is activated by a G protein, we injected starfish eggs with a PLCγ SH2 domain fusion protein that inhibits activation of PLCγ. In these eggs, Ca²⁺ release at fertilization was delayed, or with a high concentration of protein and a low concentration of sperm, completely inhibited. The PLCγSH2 protein is a specific inhibitor of PLCγ in the egg, since it did not inhibit PLCβ activation of Ca²⁺ release initiated by the serotonin 2c receptor, or activation of Ca²⁺ release by IP₃ injection. Furthermore, injection of a PLCγ SH2 domain protein mutated at its phosphotyrosine binding site, or the SH2 domains of another protein (the phosphatase SHP2), did not inhibit Ca²⁺ release at fertilization. These results indicate that during fertilization of starfish eggs, activation of phospholipase Cγ by an SH2 domain-mediated process stimulates the production of IP₃ that causes intracellular Ca²⁺ release.     

The Natural History of Biocatalytic Mechanisms

Phylogenomic analysis of the occurrence and abundance of protein domains in proteomes has recently showed that the α/β architecture is probably the oldest fold design. This holds important implications for the origins of biochemistry. Here we explore structure-function relationships addressing the use of chemical mechanisms by ancestral enzymes. We test the hypothesis that the oldest folds used the most mechanisms. We start by tracing biocatalytic mechanisms operating in metabolic enzymes along a phylogenetic timeline of the first appearance of homologous superfamilies of protein domain structures from CATH. A total of 335 enzyme reactions were retrieved from MACiE and were mapped over fold age. We define a mechanistic step type as one of the 51 mechanistic annotations given in MACiE, and each step of each of the 335 mechanisms was described using one or more of these annotations. We find that the first two folds, the P-loop containing nucleotide triphosphate hydrolase and the NAD(P)-binding Rossmann-like homologous superfamilies, were α/β architectures responsible for introducing 35% (18/51) of the known mechanistic step types. We find that these two oldest structures in the phylogenomic analysis of protein domains introduced many mechanistic step types that were later combinatorially spread in catalytic history. The most common mechanistic step types included fundamental building blocks of enzyme chemistry: “Proton transfer,” “Bimolecular nucleophilic addition,” “Bimolecular nucleophilic substitution,” and “Unimolecular elimination by the conjugate base.” They were associated with the most ancestral fold structure typical of P-loop containing nucleotide triphosphate hydrolases. Over half of the mechanistic step types were introduced in the evolutionary timeline before the appearance of structures specific to diversified organisms, during a period of architectural diversification. The other half unfolded gradually after organismal diversification and during a period that spanned ∼2 billion years of evolutionary history.     

Differences in structural elements of Bcr-Abl oncoprotein isoforms in Chronic Myelogenous Leukemia

in silico modeling, using Psipred and ExPASy servers was employed to determine the structural elements of Bcr-Abl oncoprotein (p210^BCR-ABL) isoforms, b2a2 and b3a2, expressed in Chronic Myelogenous Leukemia (CML). Both these proteins are tyrosine kinases having masses of 210-kDa and differing only by 25 amino acids coded by the b3 exonand an amino acidsubstitution (Glu903Asp). The secondary structure elements of the two proteins show differences in five α-helices and nine β-strands which relates to differences in the SH₃, SH₂, SH1 and DNA-binding domains. These differences can result in different roles played by the two isoforms in mediating signal transduction during the course of CML.     

Identification of a New Interaction Mode between the Src Homology 2 Domain of C-terminal Src Kinase (Csk) and Csk-binding Protein/Phosphoprotein Associated with Glycosphingolipid Microdomains

Proteins with Src homology 2 (SH2) domains play major roles in tyrosine kinase signaling. Structures of many SH2 domains have been studied, and the regions involved in their interactions with ligands have been elucidated. However, these analyses have been performed using short peptides consisting of phosphotyrosine followed by a few amino acids, which are described as the canonical recognition sites. Here, we report the solution structure of the SH2 domain of C-terminal Src kinase (Csk) in complex with a longer phosphopeptide from the Csk-binding protein (Cbp). This structure, together with biochemical experiments, revealed the existence of a novel binding region in addition to the canonical phosphotyrosine 314-binding site of Cbp. Mutational analysis of this second region in cells showed that both canonical and novel binding sites are required for tumor suppression through the Cbp-Csk interaction. Furthermore, the data indicate an allosteric connection between Cbp binding and Csk activation that arises from residues in the βB/βC loop of the SH2 domain.     

Unique gating properties of C. elegans ClC anion channel splice variants are determined by altered CBS domain conformation and the R-helix linker

All eukaryotic and some prokaryotic ClC anion transport proteins have extensive cytoplasmic C-termini containing two cystathionine-β-synthase (CBS) domains. CBS domain secondary structure is highly conserved and consists of two α-helices and three β-strands arranged as β1-α1-β2-β3-α2. ClC CBS domain mutations cause muscle and bone disease and alter ClC gating. However, the precise functional roles of CBS domains and the structural bases by which they regulate ClC function are poorly understood. CLH-3a and CLH-3b are C. elegans ClC anion channel splice variants with strikingly different biophysical properties. Splice variation occurs at cytoplasmic N- and C-termini and includes several amino acids that form α2 of the second CBS domain (CBS2). We demonstrate that interchanging α2 between CLH-3a and CLH-3b interchanges their gating properties. The “R-helix” of ClC proteins forms part of the ion-conducting pore and selectivity filter and is connected to the cytoplasmic C-terminus via a short stretch of cytoplasmic amino acids termed the “R-helix linker”. C-terminus conformation changes could cause R-helix structural rearrangements via this linker. X-ray structures of three ClC protein cytoplasmic C-termini suggest that α2 of CBS2 and the R-helix linker could be closely apposed and may therefore interact. We found that mutating apposing amino acids in α2 and the R-helix linker of CLH-3b was sufficient to give rise to CLH-3a-LIKE gating. We postulate that the R-helix linker interacts with CBS2 α2, and that this putative interaction provides a pathway by which cytoplasmic C-terminus conformational changes induce conformational changes in membrane domains that in turn modulate ClC function.Key words: ClC channel, chloride channel, homology model     

Intra-domain Cross-talk Regulates Serine-arginine Protein Kinase 1-dependent Phosphorylation and Splicing Function of Transformer 2β1

Transformer 2β1 (Tra2β1) is a splicing effector protein composed of a core RNA recognition motif flanked by two arginine-serine-rich (RS) domains, RS1 and RS2. Although Tra2β1-dependent splicing is regulated by phosphorylation, very little is known about how protein kinases phosphorylate these two RS domains. We now show that the serine-arginine protein kinase-1 (SRPK1) is a regulator of Tra2β1 and promotes exon inclusion in the survival motor neuron gene 2 (SMN2). To understand how SRPK1 phosphorylates this splicing factor, we performed mass spectrometric and kinetic experiments. We found that SRPK1 specifically phosphorylates 21 serines in RS1, a process facilitated by a docking groove in the kinase domain. Although SRPK1 readily phosphorylates RS2 in a splice variant lacking the N-terminal RS domain (Tra2β3), RS1 blocks phosphorylation of these serines in the full-length Tra2β1. Thus, RS2 serves two new functions. First, RS2 positively regulates binding of the central RNA recognition motif to an exonic splicing enhancer sequence, a phenomenon reversed by SRPK1 phosphorylation on RS1. Second, RS2 enhances ligand exchange in the SRPK1 active site allowing highly efficient Tra2β1 phosphorylation. These studies demonstrate that SRPK1 is a regulator of Tra2β1 splicing function and that the individual RS domains engage in considerable cross-talk, assuming novel functions with regard to RNA binding, splicing, and SRPK1 catalysis.     

BCL::Score—Knowledge Based Energy Potentials for Ranking Protein Models Represented by Idealized Secondary Structure Elements

The topology of most experimentally determined protein domains is defined by the relative arrangement of secondary structure elements, i.e. α-helices and β-strands, which make up 50–70% of the sequence. Pairing of β-strands defines the topology of β-sheets. The packing of side chains between α-helices and β-sheets defines the majority of the protein core. Often, limited experimental datasets restrain the position of secondary structure elements while lacking detail with respect to loop or side chain conformation. At the same time the regular structure and reduced flexibility of secondary structure elements make these interactions more predictable when compared to flexible loops and side chains. To determine the topology of the protein in such settings, we introduce a tailored knowledge-based energy function that evaluates arrangement of secondary structure elements only. Based on the amino acid C_β atom coordinates within secondary structure elements, potentials for amino acid pair distance, amino acid environment, secondary structure element packing, β-strand pairing, loop length, radius of gyration, contact order and secondary structure prediction agreement are defined. Separate penalty functions exclude conformations with clashes between amino acids or secondary structure elements and loops that cannot be closed. Each individual term discriminates for native-like protein structures. The composite potential significantly enriches for native-like models in three different databases of 10,000–12,000 protein models in 80–94% of the cases. The corresponding application, “BCL::ScoreProtein,” is available at www.meilerlab.org.     

Structural Analysis of Saccharomyces cerevisiae α-Galactosidase and Its Complexes with Natural Substrates Reveals New Insights into Substrate Specificity of GH27 Glycosidases

α-Galactosidases catalyze the hydrolysis of terminal α-1,6-galactosyl units from galacto-oligosaccharides and polymeric galactomannans. The crystal structures of tetrameric Saccharomyces cerevisiae α-galactosidase and its complexes with the substrates melibiose and raffinose have been determined to 1.95, 2.40, and 2.70 Å resolution. The monomer folds into a catalytic (α/β)₈ barrel and a C-terminal β-sandwich domain with unassigned function. This pattern is conserved with other family 27 glycosidases, but this enzyme presents a unique 45-residue insertion in the β-sandwich domain that folds over the barrel protecting it from the solvent and likely explaining its high stability. The structure of the complexes and the mutational analysis show that oligomerization is a key factor in substrate binding, as the substrates are located in a deep cavity making direct interactions with the adjacent subunit. Furthermore, docking analysis suggests that the supplementary domain could be involved in binding sugar units distal from the scissile bond, therefore ascribing a role in fine-tuning substrate specificity to this domain. It may also have a role in promoting association with the polymeric substrate because of the ordered arrangement that the four domains present in one face of the tetramer. Our analysis extends to other family 27 glycosidases, where some traits regarding specificity and oligomerization can be formulated on the basis of their sequence and the structures available. These results improve our knowledge on the activity of this important family of enzymes and give a deeper insight into the structural features that rule modularity and protein-carbohydrate interactions.