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.     

Interconversion of the Specificities of Human Lysosomal Enzymes Associated with Fabry and Schindler Diseases

The human lysosomal enzymes α-galactosidase (α-GAL, EC 3.2.1.22) and α-N-acetylgalactosaminidase (α-NAGAL, EC 3.2.1.49) share 46% amino acid sequence identity and have similar folds. The active sites of the two enzymes share 11 of 13 amino acids, differing only where they interact with the 2-position of the substrates. Using a rational protein engineering approach, we interconverted the enzymatic specificity of α- GAL and α-NAGAL. The engineered α-GAL (which we call α-GAL^SA) retains the antigenicity of α-GAL but has acquired the enzymatic specificity of α-NAGAL. Conversely, the engineered α-NAGAL (which we call α-NAGAL^EL) retains the antigenicity of α-NAGAL but has acquired the enzymatic specificity of the α-GAL enzyme. Comparison of the crystal structures of the designed enzyme α-GAL^SA to the wild-type enzymes shows that active sites of α-GAL^SA and α-NAGAL superimpose well, indicating success of the rational design. The designed enzymes might be useful as non-immunogenic alternatives in enzyme replacement therapy for treatment of lysosomal storage disorders such as Fabry disease.     

High degree of conservation of the enzymes synthesizing the laminin-binding glycoepitope of α-dystroglycan

The dystroglycan (DG) complex plays a pivotal role for the stabilization of muscles in Metazoa. It is formed by two subunits, extracellular α-DG and transmembrane β-DG, originating from a unique precursor via a complex post-translational maturation process. The α-DG subunit is extensively glycosylated in sequential steps by several specific enzymes and employs such glycan scaffold to tightly bind basement membrane molecules. Mutations of several of these enzymes cause an alteration of the carbohydrate structure of α-DG, resulting in severe neuromuscular disorders collectively named dystroglycanopathies. Given the fundamental role played by DG in muscle stability, it is biochemically and clinically relevant to investigate these post-translational modifying enzymes from an evolutionary perspective. A first phylogenetic history of the thirteen enzymes involved in the fabrication of the so-called ‘M3 core’ laminin-binding epitope has been traced by an overall sequence comparison approach, and interesting details on the primordial enzyme set have emerged, as well as substantial conservation in Metazoa. The optimization along with the evolution of a well-conserved enzymatic set responsible for the glycosylation of α-DG indicate the importance of the glycosylation shell in modulating the connection between sarcolemma and surrounding basement membranes to increase skeletal muscle stability, and eventually support movement and locomotion.     

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.     

Evolutionary Genomics Reveals Lineage-Specific Gene Loss and Rapid Evolution of a Sperm-Specific Ion Channel Complex: CatSpers and CatSperβ

The mammalian CatSper ion channel family consists of four sperm-specific voltage-gated Ca²⁺ channels that are crucial for sperm hyperactivation and male fertility. All four CatSper subunits are believed to assemble into a heteromultimeric channel complex, together with an auxiliary subunit, CatSperβ. Here, we report a comprehensive comparative genomics study and evolutionary analysis of CatSpers and CatSperβ, with important correlation to physiological significance of molecular evolution of the CatSper channel complex. The development of the CatSper channel complex with four CatSpers and CatSperβ originated as early as primitive metazoans such as the Cnidarian Nematostella vectensis. Comparative genomics revealed extensive lineage-specific gene loss of all four CatSpers and CatSperβ through metazoan evolution, especially in vertebrates. The CatSper channel complex underwent rapid evolution and functional divergence, while distinct evolutionary constraints appear to have acted on different domains and specific sites of the four CatSper genes. These results reveal unique evolutionary characteristics of sperm-specific Ca²⁺ channels and their adaptation to sperm biology through metazoan evolution.     

Structural Bridges through Fold Space

Several protein structure classification schemes exist that partition the protein universe into structural units called folds. Yet these schemes do not discuss how these units sit relative to each other in a global structure space. In this paper we construct networks that describe such global relationships between folds in the form of structural bridges. We generate these networks using four different structural alignment methods across multiple score thresholds. The networks constructed using the different methods remain a similar distance apart regardless of the probability threshold defining a structural bridge. This suggests that at least some structural bridges are method specific and that any attempt to build a picture of structural space should not be reliant on a single structural superposition method. Despite these differences all representations agree on an organisation of fold space into five principal community structures: all-α, all-β sandwiches, all-β barrels, α/β and α + β. We project estimated fold ages onto the networks and find that not only are the pairings of unconnected folds associated with higher age differences than bridged folds, but this difference increases with the number of networks displaying an edge. We also examine different centrality measures for folds within the networks and how these relate to fold age. While these measures interpret the central core of fold space in varied ways they all identify the disposition of ancestral folds to fall within this core and that of the more recently evolved structures to provide the peripheral landscape. These findings suggest that evolutionary information is encoded along these structural bridges. Finally, we identify four highly central pivotal folds representing dominant topological features which act as key attractors within our landscapes.     

Phylogenomic Relationships between Amylolytic Enzymes from 85 Strains of Fungi

Fungal amylolytic enzymes, including α-amylase, gluocoamylase and α-glucosidase, have been extensively exploited in diverse industrial applications such as high fructose syrup production, paper making, food processing and ethanol production. In this paper, amylolytic genes of 85 strains of fungi from the phyla Ascomycota, Basidiomycota, Chytridiomycota and Zygomycota were annotated on the genomic scale according to the classification of glycoside hydrolase (GH) from the Carbohydrate-Active enZymes (CAZy) Database. Comparisons of gene abundance in the fungi suggested that the repertoire of amylolytic genes adapted to their respective lifestyles. Amylolytic enzymes in family GH13 were divided into four distinct clades identified as heterologous α- amylases, eukaryotic α-amylases, bacterial and fungal α-amylases and GH13 α-glucosidases. Family GH15 had two branches, one for gluocoamylases, and the other with currently unknown function. GH31 α-glucosidases showed diverse branches consisting of neutral α-glucosidases, lysosomal acid α-glucosidases and a new clade phylogenetically related to the bacterial counterparts. Distribution of starch-binding domains in above fungal amylolytic enzymes was related to the enzyme source and phylogeny. Finally, likely scenarios for the evolution of amylolytic enzymes in fungi based on phylogenetic analyses were proposed. Our results provide new insights into evolutionary relationships among subgroups of fungal amylolytic enzymes and fungal evolutionary adaptation to ecological conditions.     

Experimental evolution of viruses: Microviridae as a model system

φX174 was developed as a model system for experimental studies of evolution because of its small genome size and ease of cultivation. It has been used extensively to address statistical questions about the dynamics of adaptive evolution. Molecular changes seen during experimental evolution of φX174 under a variety of conditions were compiled from 10 experiments comprising 58 lineages, where whole genomes were sequenced. A total of 667 substitutions was seen. Parallel evolution was rampant, with over 50 per cent of substitutions occurring at sites with three or more events. Comparisons of experimentally evolved sites to variation seen among wild phage suggest that at least some of the adaptive mechanisms seen in the laboratory are relevant to adaptation in nature. Elucidation of these mechanisms is aided by the availability of capsid and pro-capsid structures for φX174 and builds on years of genetic studies of the phage life history.     

Phylogeny and evolution of chloroplast tRNAs in Adoxaceae

Chloroplasts are semiautonomous organelles found in photosynthetic plants. The major functions of chloroplasts include photosynthesis and carbon fixation, which are mainly regulated by its circular genomes. In the highly conserved chloroplast genome, the chloroplast transfer RNA genes (cp tRNA) play important roles in protein translation within chloroplasts. However, the evolution of cp tRNAs remains unclear. Thus, in the present study, we investigated the evolutionary characteristics of chloroplast tRNAs in five Adoxaceae species using 185 tRNA gene sequences. In total, 37 tRNAs encoding 28 anticodons are found in the chloroplast genome in Adoxaceae species. Some consensus sequences are found within the Ψ‐stem and anticodon loop of the tRNAs. Some putative novel structures were also identified, including a new stem located in the variable region of tRNA^Tyr in a similar manner to the anticodon stem. Furthermore, phylogenetic and evolutionary analyses indicated that synonymous tRNAs may have evolved from multiple ancestors and frequent tRNA duplications during the evolutionary process may have been primarily caused by positive selection and adaptive evolution. The transition and transversion rates are uneven among different tRNA isotypes. For all tRNAs, the transition rate is greater with a transition/transversion bias of 3.13. Phylogenetic analysis of cp tRNA suggested that the type I introns in different taxa (including eukaryote organisms and cyanobacteria) share the conserved sequences “U‐U‐x2‐C” and “U‐x‐G‐x2‐T,” thereby indicating the diverse cyanobacterial origins of organelles. This detailed study of cp tRNAs in Adoxaceae may facilitate further investigations of the evolution, phylogeny, structure, and related functions of chloroplast tRNAs.     

Ero1-α and PDIs constitute a hierarchical electron transfer network of endoplasmic reticulum oxidoreductases

Ero1-α and endoplasmic reticulum (ER) oxidoreductases of the protein disulfide isomerase (PDI) family promote the efficient introduction of disulfide bonds into nascent polypeptides in the ER. However, the hierarchy of electron transfer among these oxidoreductases is poorly understood. In this paper, Ero1-α–associated oxidoreductases were identified by proteomic analysis and further confirmed by surface plasmon resonance. Ero1-α and PDI were found to constitute a regulatory hub, whereby PDI induced conformational flexibility in an Ero1-α shuttle cysteine (Cys99) facilitated intramolecular electron transfer to the active site. In isolation, Ero1-α also oxidized ERp46, ERp57, and P5; however, kinetic measurements and redox equilibrium analysis revealed that PDI preferentially oxidized other oxidoreductases. PDI accepted electrons from the other oxidoreductases via its a′ domain, bypassing the a domain, which serves as the electron acceptor from reduced glutathione. These observations provide an integrated picture of the hierarchy of cooperative redox interactions among ER oxidoreductases in mammalian cells.     

Computational exploration of the network of sequence flow between protein structures

We investigate small sequence adjustments (of one or a few amino acids) that induce large conformational transitions between distinct and stable folds of proteins. Such transitions are intriguing from evolutionary and protein‐design perspectives. They make it possible to search for ancient protein structures or to design protein switches that flip between folds and functions. A network of sequence flow between protein folds is computed for representative structures of the Protein Data Bank. The computed network is dense, on an average each structure is connected to tens of other folds. Proteins that attract sequences from a higher than expected number of neighboring folds are more likely to be enzymes and alpha/beta fold. The large number of connections between folds may reflect the need of enzymes to adjust their structures for alternative substrates. The network of the Cro family is discussed, and we speculate that capacity is an important factor (but not the only one) that determines protein evolution. The experimentally observed flip from all alpha to alpha + beta fold is examined by the network tools. A kinetic model for the transition of sequences between the folds (with only protein stability in mind) is proposed. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Patterns of Nucleotide Substitution in Mitochondrial Protein Coding Genes of Vertebrates

Maximum likelihood methods were used to study the differences in substitution rates among the four nucleotides and among different nucleotide sites in mitochondrial protein-coding genes of vertebrates. In the 1st+2nd codon position data, the frequency of nucleotide G is negatively correlated with evolutionary rates of genes, substitution rates vary substantially among sites, and the transition/transversion rate bias (R) is two to five times larger than that expected at random. Generally, largest transition biases and greatest differences in substitution rates among sites are found in the highly conserved genes. The 3rd positions in placental mammal genes exhibit strong nucleotide composition biases and the transitional rates exceed transversional rates by one to two orders of magnitude. Tamura-Nei and Hasegawa-Kishino-Yano models with gamma distributed variable rates among sites (gamma parameter, α) adequately describe the nucleotide substitution process in 1st+2nd position data. In these data, ignoring differences in substitution rates among sites leads to largest biases while estimating substitution rates. Kimura's two-parameter model with variable-rates among sites performs satisfactorily in likelihood estimation of R, α, and overall amount of evolution for 1st+2nd position data. It can also be used to estimate pairwise distances with appropriate values of α for a majority of genes.     

Recognition of the Helical Structure of β-1,4-Galactan by a New Family of Carbohydrate-binding Modules

The microbial enzymes that depolymerize plant cell wall polysaccharides, ultimately promoting energy liberation and carbon recycling, are typically complex in their modularity and often contain carbohydrate-binding modules (CBMs). Here, through analysis of an unknown module from a Thermotoga maritima endo-β-1,4-galactanase, we identify a new family of CBMs that are most frequently found appended to proteins with β-1,4-galactanase activity. Polysaccharide microarray screening, immunofluorescence microscopy, and biochemical analysis of the isolated module demonstrate the specificity of the module, here called TmCBM61, for β-1,4-linked galactose-containing ligands, making it the founding member of family CBM61. The ultra-high resolution x-ray crystal structures of TmCBM61 (0.95 and 1.4 Å resolution) in complex with β-1,4-galactotriose reveal the molecular basis of the specificity of the CBM for β-1,4-galactan. Analysis of these structures provides insight into the recognition of an unexpected helical galactan conformation through a mode of binding that resembles the recognition of starch.     

Positive selection is the main driving force for evolution of citrus canker-causing Xanthomonas

Understanding the evolutionary history and potential of bacterial pathogens is critical to prevent the emergence of new infectious bacterial diseases. Xanthomonas axonopodis subsp. citri (Xac) (synonym X. citri subsp. citri), which causes citrus canker, is one of the hardest-fought plant bacterial pathogens in US history. Here, we sequenced 21 Xac strains (14 XacA, 3 XacA* and 4 XacA^w) with different host ranges from North America and Asia and conducted comparative genomic and evolutionary analyses. Our analyses suggest that acquisition of beneficial genes and loss of detrimental genes most likely allowed XacA to infect a broader range of hosts as compared with XacA^w and XacA*. Recombination was found to have occurred frequently on the relative ancient branches, but rarely on the young branches of the clonal genealogy. The ratio of recombination/mutation ρ/θ was 0.0790±0.0005, implying that the Xac population was clonal in structure. Positive selection has affected 14% (395 out of 2822) of core genes of the citrus canker-causing Xanthomonas. The genes affected are enriched in ‘carbohydrate transport and metabolism'' and ‘DNA replication, recombination and repair'' genes (P<0.05). Many genes related to virulence, especially genes involved in the type III secretion system and effectors, are affected by positive selection, further highlighting the contribution of positive selection to the evolution of citrus canker-causing Xanthomonas. Our results suggest that both metabolism and virulence genes provide advantages to endow XacA with higher virulence and a wider host range. Our analysis advances our understanding of the genomic basis of specialization by positive selection in bacterial evolution.     

Sequence and structural analysis of 4SNc-Tudor domain protein from Takifugu Rubripes

The fugu SN4TDR protein belongs to an evolutionarily conserved family, consisting of four repeat staphylococcal nuclease-like domains (SN1-SN4) at the N-terminus followed by Tudor and SN-like domains (TSN). Sequence analysis showed that the C-terminal TSN domain is composed of a complete SN-like domain interdigitated with a Tudor domain. In despite of low level of sequence identities, five SN-like domains have a few conserved amino acids that may play essential roles in the function of the protein. Computer modeling and secondary structural prediction of the SN-like domains revealed the presence of similar structural features of β1-β2-β3-α1-β4-β5-α2-α3, which provides a structural basis for oligonucleotides binding. The loop region L_3α for binding sites between β3 and α1 of SN-like domains are different from human p100, implying the divergence in the structures of binding sites. These results indicate that fugu SN4TDR may bind methylated ligands and/or oligonucleotides through its distant domains.     

Oxidoreductases Involved in Cell Carbon Synthesis of Methanobacterium thermoautotrophicum

Cell-free extracts of Methanobacterium thermoautotrophicum were found to contain high activities of the following oxidoreductases (at 60°C): pyruvate dehydrogenase (coenzyme A acetylating), 275 nmol/min per mg of protein; α-ketoglutarate dehydrogenase (coenzyme A acylating), 100 nmol/min per mg; fumarate reductase, 360 nmol/min per mg; malate dehydrogenase, 240 nmol/min per mg; and glyceraldehyde-3-phosphate dehydrogenase, 100 nmol/min per mg. The kinetic properties (apparent V_max and K_M values), pH optimum, temperature dependence of the rate, and specificity for electron acceptors/donors of the different oxidoreductases were examined. Pyruvate dehydrogenase and α-ketoglutarate dehydrogenase were shown to be two separate enzymes specific for factor 420 rather than for nicotinamide adenine dinucleotide (NAD), NADP, or ferredoxin as the electron acceptor. Both activities catalyzed the reduction of methyl viologen with the respective α-ketoacid and a coenzyme A-dependent exchange between the carboxyl group of the α-ketoacid and CO₂. The data indicate that the two enzymes are similar to pyruvate synthase and α-ketoglutarate synthase, respectively. Fumarate reductase was found in the soluble cell fraction. This enzyme activity coupled with reduced benzyl viologen as the electron donor, but reduced factor 420, NADH, or NADPH was not effective. The cells did not contain menaquinone, thus excluding this compound as the physiological electron donor for fumarate reduction. NAD was the preferred coenzyme for malate dehydrogenase, whereas NADP was preferred for glyceraldehyde-3-phosphate dehydrogenase. The organism also possessed a factor 420-dependent hydrogenase and a factor 420-linked NADP reductase. The involvement of the described oxidoreductases in cell carbon synthesis is discussed.     

Distinctive Patterns of Evolution of the δ-Globin Gene (HBD) in Primates

In most vertebrates, hemoglobin (Hb) is a heterotetramer composed of two dissimilar globin chains, which change during development according to the patterns of expression of α- and β-globin family members. In placental mammals, the β-globin cluster includes three early-expressed genes, ε(HBE)-γ(HBG)-ψβ(HBBP1), and the late expressed genes, δ (HBD) and β (HBB). While HBB encodes the major adult β-globin chain, HBD is weakly expressed or totally silent. Paradoxically, in human populations HBD shows high levels of conservation typical of genes under strong evolutionary constraints, possibly due to a regulatory role in the fetal-to-adult switch unique of Anthropoid primates. In this study, we have performed a comprehensive phylogenetic and comparative analysis of the two adult β-like globin genes in a set of diverse mammalian taxa, focusing on the evolution and functional divergence of HBD in primates. Our analysis revealed that anthropoids are an exception to a general pattern of concerted evolution in placental mammals, showing a high level of sequence conservation at HBD, less frequent and shorter gene conversion events. Moreover, this lineage is unique in the retention of a functional GATA-1 motif, known to be involved in the control of the developmental expression of the β-like globin genes. We further show that not only the mode but also the rate of evolution of the δ-globin gene in higher primates are strictly associated with the fetal/adult β-cluster developmental switch. To gain further insight into the possible functional constraints that have been shaping the evolutionary history of HBD in primates, we calculated dN/dS (ω) ratios under alternative models of gene evolution. Although our results indicate that HBD might have experienced different selective pressures throughout primate evolution, as shown by different ω values between apes and Old World Monkeys + New World Monkeys (0.06 versus 0.43, respectively), these estimates corroborated a constrained evolution for HBD in Anthropoid lineages, which is unlikely to be related to protein function. Collectively, these findings suggest that sequence change at the δ-globin gene has been under strong selective constraints over 65 Myr of primate evolution, likely due to a regulatory role in ontogenic switches of gene expression.     

Smooth Functional Transition along a Mutational Pathway with an Abrupt Protein Fold Switch

Recent protein design experiments have demonstrated that proteins can migrate between folds through the accumulation of substitution mutations without visiting disordered or nonfunctional points in sequence space. To explore the biophysical mechanism underlying such transitions we use a three-letter continuous protein model with seven atoms per amino acid to provide realistic sequence-structure and sequence-function mappings through explicit simulation of the folding and interaction of model sequences. We start from two 16-amino-acid sequences folding into an α-helix and a β-hairpin, respectively, each of which has a preferred binding partner with 35 amino acids. We identify a mutational pathway between the two folds, which features a sharp fold switch. By contrast, we find that the transition in function is smooth. Moreover, the switch in preferred binding partner does not coincide with the fold switch. Discovery of new folds in evolution might therefore be facilitated by following fitness slopes in sequence space underpinned by binding-induced conformational switching.