Structural evidence for evolution of the beta/alpha barrel scaffold by gene duplication and fusion

Structural, kinetic, and evolutionary peculiarities of HISN3, a plant 5'-ProFAR isomerase

Histidine biosynthesis is essential for the growth and development of plants, where it occurs within chloroplasts. The eleven reactions are catalyzed by eight enzymes, known as HISN1-8, each acting sequentially. Here, we present the crystal structures of a 5'-ProFAR isomerase (HISN3) from the model legume Medicago truncatula bound to its enzymatically synthesized substrate (ProFAR) and product (PrFAR). The active site of MtHISN3 contains a sodium cation that participates in ligand recognition, a feature not observed in bacterial and fungal structures of homologous enzymes. The steady-state kinetics of wild-type MtHISN3 revealed a slightly higher turnover rate compared to its bacterial homologs. Plant HISN3 sequences contain an unusually elongated Lys60-Ser91 fragment, while deletion of the 74-80 region resulted in a 30-fold loss in catalytic efficiency compared to the wild-type. Molecular dynamics simulations suggested that the fragment facilitates product release, thereby contributing to a higher kcat. Moreover, conservation analyses suggested a non-cyanobacterial origin for plant HISN3 enzymes, which is another instance of a non-cyanobacterial enzyme in the plant histidine biosynthetic pathway. Finally, a virtual screening campaign yielded five molecules, with the energy gains ranging between -13.6 and -13.1 kcal/mol, which provide new scaffolds for the future development of herbicides.

Mutational and structural studies of (βα)8-barrel fold methylene-tetrahydropterin reductases utilizing a common catalytic mechanism

Methylene-tetrahydropterin reductases catalyze the reduction of a methylene to a methyl group bound to a reduced pterin as C1 carrier in various one-carbon (C1) metabolisms. F420-dependent methylene-tetrahydromethanopterin (methylene-H4MPT) reductase (Mer) and the flavin-independent methylene-tetrahydrofolate (methylene-H4F) reductase (Mfr) use a ternary complex mechanism for the direct transfer of a hydride from F420H2 and NAD(P)H to the respective methylene group, whereas FAD-dependent methylene-H4F reductase (MTHFR) uses FAD as prosthetic group and a ping-pong mechanism to catalyze the reduction of methylene-H4F. A ternary complex structure and a thereof derived catalytic mechanism of MTHFR is available, while no ternary complex structures of Mfr or Mer are reported. Here, Mer from Methanocaldococcus jannaschii (jMer) was heterologously produced and the crystal structures of the enzyme with and without F420 were determined. A ternary complex of jMer was modeled on the basis of the jMer-F420 structure and the ternary complex structure of MTHFR by superimposing the polypeptide after fixing hydride-transferring atoms of the flavins on each other, and by the subsequent transfer of the methyl-tetrahydropterin from MTHFR to jMer. Mutational analysis of four functional amino acids, which are similarly positioned in the three reductase structures, indicated despite the insignificant sequence identity, a common catalytic mechanism with a 5-iminium cation of methylene-tetrahydropterin as intermediate protonated by a shared glutamate. According to structural, mutational and phylogenetic analysis, the evolution of the three reductases most likely proceeds via a convergent development although a divergent scenario requiring drastic structural changes of the common ancestor cannot be completely ruled out.

Anciently duplicated genes continuously recruited to heart expression in vertebrate evolution are associated with heart chamber increase

Although gene/genome duplications in the early stage of vertebrates have been thought to provide major resources of raw genetic materials for evolutionary innovations, it is unclear whether they continuously contribute to the evolution of morphological complexity during the course of vertebrate evolution, such as the evolution from two heart chambers (fishes) to four heart chambers (mammals and birds). We addressed this issue by our heart RNA-Seq experiments combined with published data, using 13 vertebrates and one invertebrate (sea squirt, as an outgroup). Our evolutionary transcriptome analysis showed that number of ancient paralogous genes expressed in heart tends to increase with the increase of heart chamber number along the vertebrate phylogeny, in spite that most of them were duplicated at the time near to the origin of vertebrates or even more ancient. Moreover, those paralogs expressed in heart exert considerably different functions from heart-expressed singletons: the former are functionally enriched in cardiac muscle and muscle contraction-related categories, whereas the latter play more basic functions of energy generation like aerobic respiration. These findings together support the notion that recruiting anciently paralogous genes that are expressed in heart is associated with the increase of chamber number in vertebrate evolution.

A STRP-ed definition of Structured Tandem Repeats in Proteins

Tandem Repeat Proteins (TRPs) are a class of proteins with repetitive amino acid sequences that have been studied extensively for over two decades. Different features at the level of sequence, structure, function and evolution have been attributed to them by various authors. And yet many of its salient features appear only when looking at specific subclasses of protein tandem repeats. Here, we attempt to rationalize the existing knowledge on Tandem Repeat Proteins (TRPs) by pointing out several dichotomies. The emerging picture is more nuanced than generally assumed and allows us to draw some boundaries of what is not a "proper" TRP. We conclude with an operational definition of a specific subset, which we have denominated STRPs (Structural Tandem Repeat Proteins), which separates a subclass of tandem repeats with distinctive features from several other less well-defined types of repeats. We believe that this definition will help researchers in the field to better characterize the biological meaning of this large yet largely understudied group of proteins.

Why we are made of proteins and nucleic acids: Structural biology views on extraterrestrial life

Is it a miracle that life exists on the Earth, or is it a common phenomenon in the universe? If extraterrestrial organisms exist, what are they like? To answer these questions, we must understand what kinds of molecules could evolve into life, or in other words, what properties are generally required to perform biological functions and store genetic information. This review summarizes recent findings on simple ancestral proteins, outlines the basic knowledge in textbooks, and discusses the generally required properties for biological molecules from structural biology viewpoints (e.g., restriction of shapes, and types of intra- and intermolecular interactions), leading to the conclusion that proteins and nucleic acids are at least one of the simplest (and perhaps very common) forms of catalytic and genetic biopolymers in the universe. This review article is an extended version of the Japanese article, On the Origin of Life: Coevolution between RNA and Peptide, published in SEIBUTSU BUTSURI Vol. 61, p. 232-235 (2021).

Stepwise design of pseudosymmetric protein hetero-oligomers

Pseudosymmetric hetero-oligomers with three or more unique subunits with overall structural (but not sequence) symmetry play key roles in biology, and systematic approaches for generating such proteins de novo would provide new routes to controlling cell signaling and designing complex protein materials. However, the de novo design of protein hetero-oligomers with three or more distinct chains with nearly identical structures is a challenging problem because it requires the accurate design of multiple protein-protein interfaces simultaneously. Here, we describe a divide-and-conquer approach that breaks the multiple-interface design challenge into a set of more tractable symmetric single-interface redesign problems, followed by structural recombination of the validated homo-oligomers into pseudosymmetric hetero-oligomers. Starting from de novo designed circular homo-oligomers composed of 9 or 24 tandemly repeated units, we redesigned the inter-subunit interfaces to generate 15 new homo-oligomers and recombined them to make 17 new hetero-oligomers, including ABC heterotrimers, A2B2 heterotetramers, and A3B3 and A2B2C2 heterohexamers which assemble with high structural specificity. The symmetric homo-oligomers and pseudosymmetric hetero-oligomers generated for each system share a common backbone, and hence are ideal building blocks for generating and functionalizing larger symmetric assemblies.

Rosetta's Predictive Ability for Low-Affinity Ligand Binding in Fragment-Based Drug Discovery

Fragment-based drug discovery begins with the identification of small molecules with a molecular weight of usually less than 250 Da which weakly bind to the protein of interest. This technique is challenging for computational docking methods as binding is determined by only a few specific interactions. Inaccuracies in the energy function or slight deviations in the docking pose can lead to the prediction of incorrect binding or difficulties in ranking fragments in in silico screening. Here, we test RosettaLigand by docking a series of fragments to a cysteine-depleted variant of the TIM-barrel protein, HisF (UniProtKB Q9X0C6). We compare the computational results with experimental NMR spectroscopy screens. NMR spectroscopy gives details on binding affinities of individual ligands, which allows assessment of the ligand-ranking ability using RosettaLigand and also provides feedback on the location of the binding pocket, which serves as a reliable test of RosettaLigand's ability to identify plausible binding poses. From a library screen of 3456 fragments, we identified a set of 31 ligands with intrinsic affinities to HisF with dissociation constants as low as 400 μM. The same library of fragments was blindly screened in silico. RosettaLigand was able to rank binders before non-binders with an area under the curve of the receiver operating characteristics of 0.74. The docking poses observed for binders agreed with the binding pocket identified by NMR chemical shift perturbations for all fragments. Taken together, these results provide a baseline performance of RosettaLigand in a fragment-based drug discovery setting.

The Fate of Duplicated Enzymes in Prokaryotes: The Case of Isomerases

The isomerases are a unique enzymatic class of enzymes that carry out a great diversity of chemical reactions at the intramolecular level. This class comprises about 300 members, most of which are involved in carbohydrate and terpenoid/polyketide metabolism. Along with oxidoreductases and translocases, isomerases are one of the classes with the highest ratio of paralogous enzymes. Due to its relatively small number of members, it is plausible to explore it in greater detail to identify specific cases of gene duplication. Here, we present an analysis at the level of individual isomerases and identify different members that seem to be involved in duplication events in prokaryotes. As was suggested in a previous study, there is no homogeneous distribution of paralogs, but rather they accumulate into a few subcategories, some of which differ between Archaea and Bacteria. As expected, the metabolic processes with more paralogous isomerases have to do with carbohydrate metabolism but also with RNA modification (a particular case involving an rRNA-modifying isomerase is thoroughly discussed and analyzed in detail). Overall, our findings suggest that the most common fate for paralogous enzymes is the retention of the original enzymatic function, either associated with a dosage effect or with differential expression in response to changing environments, followed by subfunctionalization and, to a much lesser degree, neofunctionalization, which is consistent with what has been reported elsewhere.

Identification and Mutation Analysis of Nonconserved Residues on the TIM-Barrel Surface of GH5_5 Cellulases for Catalytic Efficiency and Stability Improvement

Exploring the potential functions of nonconserved residues on the outer side of α-helices and systematically optimizing them are pivotal for their application in protein engineering. Based on the evolutionary structural conservation analysis of GH5_5 cellulases, a practical molecular improvement strategy was developed. Highly variable sites on the outer side of the α-helices of the GH5_5 cellulase from Aspergillus niger (AnCel5A) were screened, and 14 out of the 34 highly variable sites were confirmed to exert a positive effect on the activity. After the modular combination of the positive mutations, the catalytic efficiency of the mutants was further improved. By using CMC-Na as the substrate, the catalytic efficiency and specific activity of variant AnCel5A_N193A/T300P/D307P were approximately 2.0-fold that of AnCel5A (227 ± 21 versus 451 ± 43 ml/s/mg and 1,726 ± 19 versus 3,472 ± 42 U/mg, respectively). The half-life (t_1/2) of variant AnCel5A_N193A/T300P/D307P at 75°C was 2.36 times that of AnCel5A. The role of these sites was successfully validated in other GH5_5 cellulases. Computational analyses revealed that the flexibility of the loop 6-loop 7-loop 8 region was responsible for the increased catalytic performance. This work not only illustrated the important role of rapidly evolving positions on the outer side of the α-helices of GH5_5 cellulases but also revealed new insights into engineering the proteins that nature left as clues for us to find. IMPORTANCE A comprehensive understanding of the residues on the α-helices of the GH5_5 cellulases is important for catalytic efficiency and stability improvement. The main objective of this study was to use the evolutionary conservation and plasticity of the TIM-barrel fold to probe the relationship between nonconserved residues on the outer side of the α-helices and the catalytic efficiency of GH5_5 cellulases by conducting structure-guided protein engineering. By using a four-step nonconserved residue screening strategy, the functional role of nonconserved residues on the outer side of the α-helices was effectively identified, and a variant with superior performance and capability was constructed. Hence, this study proved the effectiveness of this strategy in engineering GH5_5 cellulases and provided a potential competitor for industrial applications. Furthermore, this study sheds new light on engineering TIM-barrel proteins.

Gene fusion as an important mechanism to generate new genes in the genus Oryza

Background

Events of gene fusion have been reported in several organisms. However, the general role of gene fusion as part of new gene origination remains unknown.

Results

We conduct genome-wide interrogations of four Oryza genomes by designing and implementing novel pipelines to detect fusion genes. Based on the phylogeny of ten plant species, we detect 310 fusion genes across four Oryza species. The estimated rate of origination of fusion genes in the Oryza genus is as high as 63 fusion genes per species per million years, which is fixed at 16 fusion genes per species per million years and much higher than that in flies. By RNA sequencing analysis, we find more than 44% of the fusion genes are expressed and 90% of gene pairs show strong signals of purifying selection. Further analysis of CRISPR/Cas9 knockout lines indicates that newly formed fusion genes regulate phenotype traits including seed germination, shoot length and root length, suggesting the functional significance of these genes.

Conclusions

We detect new fusion genes that may drive phenotype evolution in Oryza. This study provides novel insights into the genome evolution of Oryza.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02696-w.

Complex Loop Dynamics Underpin Activity, Specificity, and Evolvability in the (βα)8 Barrel Enzymes of Histidine and Tryptophan Biosynthesis

Enzymes are conformationally dynamic, and their dynamical properties play an important role in regulating their specificity and evolvability. In this context, substantial attention has been paid to the role of ligand-gated conformational changes in enzyme catalysis; however, such studies have focused on tremendously proficient enzymes such as triosephosphate isomerase and orotidine 5'-monophosphate decarboxylase, where the rapid (μs timescale) motion of a single loop dominates the transition between catalytically inactive and active conformations. In contrast, the (βα)8-barrels of tryptophan and histidine biosynthesis, such as the specialist isomerase enzymes HisA and TrpF, and the bifunctional isomerase PriA, are decorated by multiple long loops that undergo conformational transitions on the ms (or slower) timescale. Studying the interdependent motions of multiple slow loops, and their role in catalysis, poses a significant computational challenge. This work combines conventional and enhanced molecular dynamics simulations with empirical valence bond simulations to provide rich details of the conformational behavior of the catalytic loops in HisA, PriA, and TrpF, and the role of their plasticity in facilitating bifunctionality in PriA and evolved HisA variants. In addition, we demonstrate that, similar to other enzymes activated by ligand-gated conformational changes, loops 3 and 4 of HisA and PriA act as gripper loops, facilitating the isomerization of the large bulky substrate ProFAR, albeit now on much slower timescales. This hints at convergent evolution on these different (βα)8-barrel scaffolds. Finally, our work reemphasizes the potential of engineering loop dynamics as a tool to artificially manipulate the catalytic repertoire of TIM-barrel proteins.

Effect of Non-Lethal Selection on Spontaneous Revertants of Frameshift Mutations: The Escherichia coli hisF Case

Microorganisms possess the potential to adapt to fluctuations in environmental parameters, and their evolution is driven by the continuous generation of mutations. The reversion of auxotrophic mutations has been widely studied; however, little is known about the reversion of frameshift mutations resulting in amino acid auxotrophy and on the structure and functioning of the protein encoded by the revertant mutated gene. The aims of this work were to analyze the appearance of reverse mutations over time and under different selective pressures and to investigate revertant enzymes' three-dimensional structures and their correlation with a different growth ability. Escherichia coli FB182 strain, carrying the hisF892 single nucleotide deletion resulting in histidine auxotrophy, was subjected to different selective pressures, and revertant mutants were isolated and characterized. The obtained results allowed us to identify different indels of different lengths located in different positions in the hisF gene, and relations with the incubation time and the selective pressure applied were observed. Moreover, the structure of the different mutant proteins was consistent with the respective revertant ability to grow in absence of histidine, highlighting a correlation between the mutations and the catalytic activity of the mutated HisF enzyme.

Simulations reveal the key role of Arg15 in the promiscuous activity in the HisA enzyme

The HisA enzyme catalyzes the first step of histidine biosynthesis via the Amadori rearrangement of the substrate ProFAR. Since it possesses the most conserved and ancient TIM-barrel fold, it provides an ideal framework for bioengineering of a new function from ancestral enzymes. In the present study, first, the catalytic mechanism of HisA biosynthesis was elucidated using hybrid Quantum Mechanical/Molecular Mechanical calculations, and thereafter, key residues contributing towards the promiscuity for TrpF activity were revealed using several MD simulations of a wild type enzyme and its variant with the native (ProFAR) and promiscuous (PRA) substrates. Our study reveals that the two loops (βα)₁ and (βα)₅ on the catalytic site of the HisA enzyme have incredible adaptability for the native and promiscuous substrates. The conformational interplay between these two loops is substrate driven and precise bioengineering targeting these loops is key to the emergence of new functions. Furthermore, the study reveals a key role of the Arg 15 residue which is close to the catalytic center of the enzyme in the bifunctionality of the HisA enzyme by increasing the loop flexibility. Therefore, our study provides crucial information for future bioengineering work to use the HisA enzyme as a scaffold for new enzymatic activity.

Gene Fusion and Directed Evolution to Break Structural Symmetry and Boost Catalysis by an Oligomeric C-C Bond-Forming Enzyme

Gene duplication and fusion are among the primary natural processes that generate new proteins from simpler ancestors. Here we adopted this strategy to evolve a promiscuous homohexameric 4-oxalocrotonate tautomerase (4-OT) into an efficient biocatalyst for enantioselective Michael reactions. We first designed a tandem-fused 4-OT to allow independent sequence diversification of adjacent subunits by directed evolution. This fused 4-OT was then subjected to eleven rounds of directed evolution to give variant 4-OT(F11), which showed an up to 320-fold enhanced activity for the Michael addition of nitromethane to cinnamaldehydes. Crystallographic analysis revealed that 4-OT(F11) has an unusual asymmetric trimeric architecture in which one of the monomers is flipped 180° relative to the others. This gene duplication and fusion strategy to break structural symmetry is likely to become an indispensable asset of the enzyme engineering toolbox, finding wide use in engineering oligomeric proteins.

A newly introduced salt bridge cluster improves structural and biophysical properties of de novo TIM barrels

Protein stability can be fine‐tuned by modifying different structural features such as hydrogen‐bond networks, salt bridges, hydrophobic cores, or disulfide bridges. Among these, stabilization by salt bridges is a major challenge in protein design and engineering since their stabilizing effects show a high dependence on the structural environment in the protein, and therefore are difficult to predict and model. In this work, we explore the effects on structure and stability of an introduced salt bridge cluster in the context of three different de novo TIM barrels. The salt bridge variants exhibit similar thermostability in comparison with their parental designs but important differences in the conformational stability at 25°C can be observed such as a highly stabilizing effect for two of the proteins but a destabilizing effect to the third. Analysis of the formed geometries of the salt bridge cluster in the crystal structures show either highly ordered salt bridge clusters or only single salt bridges. Rosetta modeling of the salt bridge clusters results in a good prediction of the tendency on stability changes but not the geometries observed in the three‐dimensional structures. The results show that despite the similarities in protein fold, the salt bridge clusters differently influence the structural and stability properties of the de novo TIM barrel variants depending on the structural background where they are introduced.

PDB Code(s): 7OSU, 7OT7, 7OSV, 7OT8 and 7P12;

Protein Assembly by Design

Proteins are nature's primary building blocks for the construction of sophisticated molecular machines and dynamic materials, ranging from protein complexes such as photosystem II and nitrogenase that drive biogeochemical cycles to cytoskeletal assemblies and muscle fibers for motion. Such natural systems have inspired extensive efforts in the rational design of artificial protein assemblies in the last two decades. As molecular building blocks, proteins are highly complex, in terms of both their three-dimensional structures and chemical compositions. To enable control over the self-assembly of such complex molecules, scientists have devised many creative strategies by combining tools and principles of experimental and computational biophysics, supramolecular chemistry, inorganic chemistry, materials science, and polymer chemistry, among others. Owing to these innovative strategies, what started as a purely structure-building exercise two decades ago has, in short order, led to artificial protein assemblies with unprecedented structures and functions and protein-based materials with unusual properties. Our goal in this review is to give an overview of this exciting and highly interdisciplinary area of research, first outlining the design strategies and tools that have been devised for controlling protein self-assembly, then describing the diverse structures of artificial protein assemblies, and finally highlighting the emergent properties and functions of these assemblies.

Seven Amino Acid Types Suffice to Create the Core Fold of RNA Polymerase

The extant complex proteins must have evolved from ancient short and simple ancestors. The double-ψ β-barrel (DPBB) is one of the oldest protein folds and conserved in various fundamental enzymes, such as the core domain of RNA polymerase. Here, by reverse engineering a modern DPBB domain, we reconstructed its plausible evolutionary pathway started by "interlacing homodimerization" of a half-size peptide, followed by gene duplication and fusion. Furthermore, by simplifying the amino acid repertoire of the peptide, we successfully created the DPBB fold with only seven amino acid types (Ala, Asp, Glu, Gly, Lys, Arg, and Val), which can be coded by only GNN and ARR (R = A or G) codons in the modern translation system. Thus, the DPBB fold could have been materialized by the early translation system and genetic code.

Role of a high centrality residue in protein dynamics and thermal stability

Centralities determined from Residue Interaction Networks (RIN) in proteins have been used to predict aspects of their structure and dynamics. Here, we correlate the Eigenvector Centrality (Ec) with the rate constant for thermal denaturation (kden) of the HisF protein from Thermotoga maritima based on 12 single alanine substitution mutants. The molecular basis for this correlation was further explored by studying a mutant containing a replacement of a high Ec residue, Y182A, which displayed increased kden at 80 °C. The crystallographic structure of this mutant showed few changes, mostly in two flexible loops. The 1H-15N -HSQC showed only subtle changes of cross peak positions for residues located near the mutation site and scattered throughout the structure. However, the comparison of the RIN showed that Y182 is the vertex of a set of high centrality residues that spreads throughout the HisF structure, which is lacking in the mutant. Cross-correlation displacements of Cα calculated from a molecular dynamics simulation at different temperatures showed that the Y182A mutation reduced the correlated movements in the HisF structure above 70 °C. 1H-15N NMR chemical shift covariance using temperature as perturbation were consistent with these results. In conclusion the increase in temperature drives the structure of the mutant HisF-Y182A into a less connected state, richer in non-concerted motions, located predominantly in the C-terminal half of the protein where Y182 is placed. Conversely, wild-type HisF responds to increased temperature as a single unit. Hence the replacement of a high Ec residue alters the distribution of thermal energy through HisF structure.

Evolution, folding, and design of TIM barrels and related proteins

Proteins are chief actors in life that perform a myriad of exquisite functions. This diversity has been enabled through the evolution and diversification of protein folds. Analysis of sequences and structures strongly suggest that numerous protein pieces have been reused as building blocks and propagated to many modern folds. This information can be traced to understand how the protein world has diversified. In this review, we discuss the latest advances in the analysis of protein evolutionary units, and we use as a model system one of the most abundant and versatile topologies, the TIM-barrel fold, to highlight the existing common principles that interconnect protein evolution, structure, folding, function, and design.

Fungal GH25 muramidases: New family members with applications in animal nutrition and a crystal structure at 0.78Å resolution

Muramidases/lysozymes hydrolyse the peptidoglycan component of the bacterial cell wall. They are found in many of the glycoside hydrolase (GH) families. Family GH25 contains muramidases/lysozymes, known as CH type lysozymes, as they were initially discovered in the Chalaropsis species of fungus. The characterized enzymes from GH25 exhibit both β-1,4-N-acetyl- and β-1,4-N,6-O-diacetylmuramidase activities, cleaving the β-1,4-glycosidic bond between N-acetylmuramic acid (NAM) and N-acetylglucosamine (NAG) moieties in the carbohydrate backbone of bacterial peptidoglycan. Here, a set of fungal GH25 muramidases were identified from a sequence search, cloned and expressed and screened for their ability to digest bacterial peptidoglycan, to be used in a commercial application in chicken feed. The screen identified the enzyme from Acremonium alcalophilum JCM 736 as a suitable candidate for this purpose and its relevant biochemical and biophysical and properties are described. We report the crystal structure of the A. alcalophilum enzyme at atomic, 0.78 Å resolution, together with that of its homologue from Trichobolus zukalii at 1.4 Å, and compare these with the structures of homologues. GH25 enzymes offer a new solution in animal feed applications such as for processing bacterial debris in the animal gut.

Development and applications of artificial symmetrical proteins

Since the determination of the first molecular models of proteins there has been interest in creating proteins artificially, but such methods have only become widely successful in the last decade. Gradual improvements over a long period of time have now yielded numerous examples of non-natural proteins, many of which are built from repeated elements. In this review we discuss the design of such symmetrical proteins and their various applications in chemistry and medicine.

Investigating the amino acid sequences of membrane bound dihydroorotate:quinone oxidoreductases (DHOQOs): Structural and functional implications

Dihydroorotate:quinone oxidoreductases (DHOQOs) are membrane bound enzymes responsible for oxidizing dihydroorotate (DHO) to orotate with concomitant reduction of quinone to quinol. They have FMN as prosthetic group and are part of the monotopic quinone reductase superfamily. These enzymes are also members of the dihydroorotate dehydrogenases (DHODHs) family, which besides membrane bound DHOQOs, class 2, includes soluble enzymes which reduce either NAD⁺ or fumarate, class 1. As key enzymes in both the de novo pyrimidine biosynthetic pathway as well as in the energetic metabolism, inhibitors of DHOQOs have been investigated as leads for therapeutics in cancer, immunological disorders and bacterial/viral infections. This work is a thorough bioinformatic approach on the structural conservation and taxonomic distribution of DHOQOs. We explored previously established structural/functional hallmarks of these enzymes, while searching for uncharacterized common elements. We also discuss the cellular role of DHOQOs and organize the identified protein sequences within six sub-classes 2A to 2F, according to their taxonomic origin and sequence traits. We concluded that DHOQOs are present in Archaea, Eukarya and Bacteria, including the first recognition in Gram-positive organisms. DHOQOs can be the single dihydroorotate dehydrogenase encoded in the genome of a species, or they can coexist with other DHODHs, as the NAD⁺ or fumarate reducing enzymes. Furthermore, we show that the type of catalytic base present in the active site is not an absolute criterium to distinguish between class 1 and class 2 enzymes. We propose the existence of a quinone binding motif ("ExAH") adjacent to a hydrophobic cavity present in the membrane interacting N-terminal domain.

Exploring the role of the histidine biosynthetic hisF gene in cellular metabolism and in the evolution of (ancestral) genes: from LUCA to the extant (micro)organisms

Histidine biosynthesis is an ancestral pathway that was assembled before the appearance of the Last Universal Common Ancestor; afterwards, it remained unaltered in all the extant histidine-synthesizing (micro)organisms. It is a metabolic cross-road interconnecting histidine biosynthesis to nitrogen metabolism and the de novo synthesis of purines. This interconnection is due to the reaction catalyzed by the products of hisH and hisF genes. The latter gene is an excellent model to study which trajectories have been followed by primordial cells to build the first metabolic routes, since its evolution is the result of different molecular rearrangement events, i.e. gene duplication, gene fusion, gene elongation, and domain shuffling. Additionally, this review summarizes data concerning the involvement of hisF and its product in other different cellular processes, revealing that HisF very likely plays a role also in cell division control and involvement in virulence and nodule development in different bacteria. From the metabolic viewpoint, these results suggest that HisF plays a central role in cellular metabolism, highlighting the interconnections of different metabolic pathways.

The Role of Gene Elongation in the Evolution of Histidine Biosynthetic Genes

Gene elongation is a molecular mechanism consisting of an in-tandem duplication of a gene and divergence and fusion of the two copies, resulting in a gene constituted by two divergent paralogous modules. The aim of this work was to evaluate the importance of gene elongation in the evolution of histidine biosynthetic genes and to propose a possible evolutionary model for some of them. Concerning the genes hisA and hisF, which code for two homologous (β/α)₈-barrels, it has been proposed that the two extant genes could be the result of a cascade of gene elongation/domain shuffling events starting from an ancestor gene coding for just one (β/α) module. A gene elongation event has also been proposed for the evolution of hisB and hisD; structural analyses revealed the possibility of an early elongation event, resulting in the repetition of modules. Furthermore, it is quite possible that the gene elongations responsible for the evolution of the four proteins occurred before the earliest phylogenetic divergence. In conclusion, gene elongation events seem to have played a crucial role in the evolution of the histidine biosynthetic pathway, and they may have shaped the structures of many genes during the first steps of their evolution.

Ab initio folding of a trefoil-fold motif reveals structural similarity with a β-propeller blade motif

Many protein architectures exhibit evidence of internal rotational symmetry postulated to be the result of gene duplication/fusion events involving a primordial polypeptide motif. A common feature of such structures is a domain-swapped arrangement at the interface of the N- and C-termini motifs and postulated to provide cooperative interactions that promote folding and stability. De novo designed symmetric protein architectures have demonstrated an ability to accommodate circular permutation of the N- and C-termini in the overall architecture; however, the folding requirement of the primordial motif is poorly understood, and tolerance to circular permutation is essentially unknown. The β-trefoil protein fold is a threefold-symmetric architecture where the repeating ~42-mer "trefoil-fold" motif assembles via a domain-swapped arrangement. The trefoil-fold structure in isolation exposes considerable hydrophobic area that is otherwise buried in the intact β-trefoil trimeric assembly. The trefoil-fold sequence is not predicted to adopt the trefoil-fold architecture in ab initio folding studies; rather, the predicted fold is closely related to a compact "blade" motif from the β-propeller architecture. Expression of a trefoil-fold sequence and circular permutants shows that only the wild-type N-terminal motif definition yields an intact β-trefoil trimeric assembly, while permutants yield monomers. The results elucidate the folding requirements of the primordial trefoil-fold motif, and also suggest that this motif may sample a compact conformation that limits hydrophobic residue exposure, contains key trefoil-fold structural features, but is more structurally homologous to a β-propeller blade motif.

Conserved genomic neighborhood is a strong but no perfect indicator for a direct interaction of microbial gene products

BACKGROUND

The order of genes in bacterial genomes is not random; for example, the products of genes belonging to an operon work together in the same pathway. The cotranslational assembly of protein complexes is deemed to conserve genomic neighborhoods even stronger than a common function. This is why a conserved genomic neighborhood can be utilized to predict, whether gene products form protein complexes.

RESULTS

We were interested to assess the performance of a neighborhood-based classifier that analyzes a large number of genomes. Thus, we determined for the genes encoding the subunits of 494 experimentally verified hetero-dimers their local genomic context. In order to generate phylogenetically comprehensive genomic neighborhoods, we utilized the tools offered by the Enzyme Function Initiative. For each subunit, a sequence similarity network was generated and the corresponding genome neighborhood network was analyzed to deduce the most frequent gene product. This was predicted as interaction partner, if its abundance exceeded a threshold, which was the frequency giving rise to the maximal Matthews correlation coefficient. For the threshold of 16%, the true positive rate was 45%, the false positive rate 0.06%, and the precision 55%. For approximately 20% of the subunits, the interaction partner was not found in a neighborhood of ± 10 genes.

CONCLUSIONS

Our phylogenetically comprehensive analysis confirmed that complex formation is a strong evolutionary factor that conserves genome neighborhoods. On the other hand, for 55% of the cases analyzed here, classification failed. Either, the interaction partner was not present in a ± 10 gene window or was not the most frequent gene product.

More Than One HMG-CoA Lyase: The Classical Mitochondrial Enzyme Plus the Peroxisomal and the Cytosolic Ones

There are three human enzymes with HMG-CoA lyase activity that are able to synthesize ketone bodies in different subcellular compartments. The mitochondrial HMG-CoA lyase was the first to be described, and catalyzes the cleavage of 3-hydroxy-3-methylglutaryl CoA to acetoacetate and acetyl-CoA, the common final step in ketogenesis and leucine catabolism. This protein is mainly expressed in the liver and its function is metabolic, since it produces ketone bodies as energetic fuels when glucose levels are low. Another isoform is encoded by the same gene for the mitochondrial HMG-CoA lyase (HMGCL), but it is located in peroxisomes. The last HMG-CoA lyase to be described is encoded by a different gene, HMGCLL1, and is located in the cytosolic side of the endoplasmic reticulum membrane. Some activity assays and tissue distribution of this enzyme have shown the brain and lung as key tissues for studying its function. Although the roles of the peroxisomal and cytosolic HMG-CoA lyases remain unknown, recent studies highlight the role of ketone bodies in metabolic remodeling, homeostasis, and signaling, providing new insights into the molecular and cellular function of these enzymes.

Identification and Analysis of Long Repeats of Proteins at the Domain Level

Amino acid repeats play an important role in the structure and function of proteins. Analysis of long repeats in protein sequences enables one to understand their abundance, structure and function in the protein universe. In the present study, amino acid repeats of length >50 (long repeats) were identified in a non-redundant set of UniProt sequences using the RADAR program. The underlying structures and functions of these long repeats were carried out using the Gene3D for structural domains, Pfam for functional domains and enzyme and non-enzyme functional classification for catalytic and binding of the proteins. From a structural perspective, these long repeats seem to predominantly occur in certain architectures such as sandwich, bundle, barrel, and roll and within these architectures abundant in the superfolds. The lengths of the repeats within each fold are not uniform exhibiting different structures for different functions. We also observed that long repeats are in the domain regions of the family and are involved in the function of the proteins. After grouping based on enzyme and non-enzyme classes, we observed the abundant occurrence of long repeats in specific catalytic and binding of the proteins. In this study, we have analyzed the occurrence of long repeats in the protein sequence universe apart from well-characterized short tandem repeats in sequences and their structures and functions of the proteins at the domain level. The present study suggests that long repeats may play an important role in the structure and function of domains of the proteins.

Activity and Thermostability of GH5 Endoglucanase Chimeras from Mesophilic and Thermophilic Parents

Cellulases from glycoside hydrolase family 5 (GH5) are key endoglucanase enzymes in the degradation of diverse polysaccharide substrates and are used in industrial enzyme cocktails to break down biomass. The GH5 family shares a canonical (βα)8-barrel structure, where each (βα) module is essential for the enzyme's stability and activity. Despite their shared topology, the thermostability of GH5 endoglucanase enzymes can vary significantly, and highly thermostable variants are often sought for industrial applications. Based on the previously characterized thermophilic GH5 endoglucanase Egl5A from Talaromyces emersonii (TeEgl5A), which has an optimal temperature of 90°C, we created 10 hybrid enzymes with elements of the mesophilic endoglucanase Cel5 from Stegonsporium opalus (SoCel5) to determine which elements are responsible for enhanced thermostability. Five of the expressed hybrid enzymes exhibit enzyme activity. Two of these hybrids exhibited pronounced increases in the temperature optimum (10 and 20°C), the temperature at which the protein lost 50% of its activity (T50) (15 and 19°C), and the melting temperature (Tm ) (16.5 and 22.9°C) and extended half-lives (t1/2) (∼240- and 650-fold at 55°C) relative to the values for the mesophilic parent enzyme and demonstrated improved catalytic efficiency on selected substrates. The successful hybridization strategies were validated experimentally in another GH5 endoglucanase, Cel5 from Aspergillus niger (AnCel5), which demonstrated a similar increase in thermostability. Based on molecular dynamics (MD) simulations of both the SoCel5 and TeEgl5A parent enzymes and their hybrids, we hypothesize that improved hydrophobic packing of the interface between α2 and α3 is the primary mechanism by which the hybrid enzymes increase their thermostability relative to that of the mesophilic parent SoCel5.IMPORTANCE Thermal stability is an essential property of enzymes in many industrial biotechnological applications, as high temperatures improve bioreactor throughput. Many protein engineering approaches, such as rational design and directed evolution, have been employed to improve the thermal properties of mesophilic enzymes. Structure-based recombination has also been used to fuse TIM barrel fragments, and even fragments from unrelated folds, to generate new structures. However, little research has been done on GH5 endoglucanases. In this study, two GH5 endoglucanases exhibiting TIM barrel structure, SoCel5 and TeEgl5A, with different thermal properties, were hybridized to study the roles of different (βα) motifs. This work illustrates the role that structure-guided recombination can play in helping to identify sequence function relationships within GH5 enzymes by supplementing natural diversity with synthetic diversity.

Computational design of symmetrical eight-bladed β-propeller proteins

β-Propeller proteins form one of the largest families of protein structures, with a pseudo-symmetrical fold made up of subdomains called blades. They are not only abundant but are also involved in a wide variety of cellular processes, often by acting as a platform for the assembly of protein complexes. WD40 proteins are a subfamily of propeller proteins with no intrinsic enzymatic activity, but their stable, modular architecture and versatile surface have allowed evolution to adapt them to many vital roles. By computationally reverse-engineering the duplication, fusion and diversification events in the evolutionary history of a WD40 protein, a perfectly symmetrical homologue called Tako8 was made. If two or four blades of Tako8 are expressed as single polypeptides, they do not self-assemble to complete the eight-bladed architecture, which may be owing to the closely spaced negative charges inside the ring. A different computational approach was employed to redesign Tako8 to create Ika8, a fourfold-symmetrical protein in which neighbouring blades carry compensating charges. Ika2 and Ika4, carrying two or four blades per subunit, respectively, were found to assemble spontaneously into a complete eight-bladed ring in solution. These artificial eight-bladed rings may find applications in bionanotechnology and as models to study the folding and evolution of WD40 proteins.

Artificial Light Regulation of an Allosteric Bienzyme Complex by a Photosensitive Ligand

The artificial regulation of proteins by light is an emerging subdiscipline of synthetic biology. Here, we used this concept to photocontrol both catalysis and allostery within the heterodimeric enzyme complex imidazole glycerol phosphate synthase (ImGP-S). ImGP-S consists of the cyclase subunit HisF and the glutaminase subunit HisH, which is allosterically stimulated by substrate binding to HisF. We show that a light-sensitive diarylethene (1,2-dithienylethene, DTE)-based competitive inhibitor in its ring-open state binds with low micromolar affinity to the cyclase subunit and displaces its substrate from the active site. As a consequence, catalysis by HisF and allosteric stimulation of HisH are impaired. Following UV-light irradiation, the DTE ligand adopts its ring-closed state and loses affinity for HisF, restoring activity and allostery. Our approach allows for the switching of ImGP-S activity and allostery during catalysis and appears to be generally applicable for the light regulation of other multienzyme complexes.

Evolutionary convergence in the biosyntheses of the imidazole moieties of histidine and purines

Background

The imidazole group is an ubiquitous chemical motif present in several key types of biomolecules. It is a structural moiety of purines, and plays a central role in biological catalysis as part of the side-chain of histidine, the amino acid most frequently found in the catalytic site of enzymes. Histidine biosynthesis starts with both ATP and the pentose phosphoribosyl pyrophosphate (PRPP), which is also the precursor for the de novo synthesis of purines. These two anabolic pathways are also connected by the imidazole intermediate 5-aminoimidazole-4-carboxamide ribotide (AICAR), which is synthesized in both routes but used only in purine biosynthesis. Rather surprisingly, the imidazole moieties of histidine and purines are synthesized by different, non-homologous enzymes. As discussed here, this phenomenon can be understood as a case of functional molecular convergence.

Results

In this work, we analyze these polyphyletic processes and argue that the independent origin of the corresponding enzymes is best explained by the differences in the function of each of the molecules to which the imidazole moiety is attached. Since the imidazole present in histidine is a catalytic moiety, its chemical arrangement allows it to act as an acid or a base. On the contrary, the de novo biosynthesis of purines starts with an activated ribose and all the successive intermediates are ribotides, with the key β-glycosidic bondage joining the ribose and the imidazole moiety. This prevents purine ribonucleotides to exhibit any imidazole-dependent catalytic activity, and may have been the critical trait for the evolution of two separate imidazole-synthesizing-enzymes. We also suggest that, in evolutionary terms, the biosynthesis of purines predated that of histidine.

Conclusions

As reviewed here, other biosynthetic routes for imidazole molecules are also found in extant metabolism, including the autocatalytic cyclization that occurs during the formation of creatinine from creatine phosphate, as well as the internal cyclization of the Ala-Ser-Gly motif of some members of the ammonia-lyase and aminomutase families, that lead to the MIO cofactor. The diversity of imidazole-synthesizing pathways highlights the biological significance of this key chemical group, whose biosyntheses evolved independently several times.

Search for independent (β/α)4 subdomains in a (β/α)8 barrel β-glucosidase

Proteins that fold as (β/α)₈ barrels are thought to have evolved from half-barrels that underwent duplication and fusion events. The evidence is particularly clear for small barrels, which have almost identical halves. Additionally, computational calculations of the thermodynamic stability of these structures in the presence of denaturants have revealed that (β/α)₈ barrels contain two subunits or domains corresponding to half-barrels. Hence, within (β/α)₈ barrels, half-barrels are self-contained units. Here, we tested this hypothesis using β-glucosidase from the bacterium Thermotoga maritima (bglTm), which has a (β/α)₈ barrel structure. Mutations were introduced to disrupt the noncovalent contacts between its halves and reveal the presence of two domains within bglTm, thus resulting in the creation of mutants T1 (containing W12A and I217A mutations) and T2 (containing W12A, H195A, I217A and F404A mutations). Mutants T1 and T2 were properly folded, as indicated by their fluorescence spectra and enzyme kinetic parameters. T1 and wild-type bglTm were equally stable, as shown by the results of thermal inactivation, differential scanning fluorimetry and guanidine hydrochloride denaturation experiments. However, T2 showed a first-order inactivation at 80°C, a single melting temperature of 82°C and only one transition concentration (c₅₀) in 2.4 M guanidine hydrochloride. Additionally, T1 and T2 exhibited a cooperative denaturation process that followed a two-state model (m-values equal to 1.4 and 1.6 kcal/mol/M, respectively), similar to that of wild-type bglTm (1.2 kcal/mol/M). Hence, T1 and T2 each denatured as a single unit, although they contained different degrees of disruption between their halves. In conclusion, bglTm halves are equivalent in terms of their thermal and chemical stability; thus, their separate contributions to (β/α)₈ barrel unfolding cannot be disentangled.

A metagenome-derived thermostable β-glucanase with an unusual module architecture which defines the new glycoside hydrolase family GH148

The discovery of novel and robust enzymes for the breakdown of plant biomass bears tremendous potential for the development of sustainable production processes in the rapidly evolving new bioeconomy. By functional screening of a metagenomic library from a volcano soil sample a novel thermostable endo-β-glucanase (EngU) which is unusual with regard to its module architecture and cleavage specificity was identified. Various recombinant EngU variants were characterized. Assignment of EngU to an existing glycoside hydrolase (GH) family was not possible. Two regions of EngU showed weak sequence similarity to proteins of the GH clan GH-A, and acidic residues crucial for catalytic activity of EngU were identified by mutation. Unusual, a carbohydrate-binding module (CBM4) which displayed binding affinity for β-glucan, lichenin and carboxymethyl-cellulose was found as an insertion between these two regions. EngU hydrolyzed β-1,4 linkages in carboxymethyl-cellulose, but displayed its highest activity with mixed linkage (β-1,3-/β-1,4-) glucans such as barley β-glucan and lichenin, where in contrast to characterized lichenases cleavage occurred predominantly at the β-1,3 linkages of C4-substituted glucose residues. EngU and numerous related enzymes with previously unknown function represent a new GH family of biomass-degrading enzymes within the GH-A clan. The name assigned to the new GH family is GH148.

From ancestral peptides to designed proteins

The diversity of modern proteins arose through the combinatorial shuffling and differentiation of a limited number of autonomously folding domain prototypes, but the origin of these prototypes themselves has long remained poorly understood. In recent years, the proposal that they originated by repetition, accretion, and recombination from an ancestral set of peptides, which evolved as cofactors of RNA-based replication and catalysis, has gained wide acceptance, supported by the systematic identification of such ancestral peptides and the experimental recapitulation of the mechanisms by which they could have yielded the first folded proteins. Inspired by this evolutionary process, protein engineers have seized on design from pre-optimized peptide components as a powerful approach to generating proteins with novel topology and functionality.

Combinatorial Library Based on Restriction Enzyme-mediated Modular Assembly

Combinatorial approaches in directed evolution were proven to be more efficient for exploring sequence space and innovating function of protein. Here, we presented the modular assembly of secondary structures (MASS) for constructing a combinatorial library. In this approach, secondary structure elements were extracted from natural existing protein. The common linkers were flanking secondary structure elements, and then secondary structure elements were digested by Hinf I restriction endonuclease that was used in the construction of combinatorial library for the first time. The digested DNA fragments were randomly ligated in the sense orientation, then in sequence to be amplified by PCR and transformation. This approach showed that different DNA fragments without homologous sequences could be randomly assembled to create significant sequence space. With the structure analysis of recombinants, it would be beneficial to the rational design, even to the design of protein de novo, and to evolve any genetic part or circuit.

Development of a motif-based topology-independent structure comparison method to identify evolutionarily related folds

Structure conservation, functional similarities, and homologous relationships that exist across diverse protein topologies suggest that some regions of the protein fold universe are continuous. However, the current structure classification systems are based on hierarchical organizations, which cannot accommodate structural relationships that span fold definitions. Here, we describe a novel, super-secondary-structure motif-based, topology-independent structure comparison method (SmotifCOMP) that is able to quantitatively identify structural relationships between disparate topologies. The basis of SmotifCOMP is a systematically defined super-secondary-structure motif library whose representative geometries are shown to be saturated in the Protein Data Bank and exhibit a unique distribution within the known folds. SmotifCOMP offers a robust and quantitative technique to compare domains that adopt different topologies since the method does not rely on a global superposition. SmotifCOMP is used to perform an exhaustive comparison of the known folds and the identified relationships are used to produce a nonhierarchical representation of the fold space that reflects the notion of a continuous and connected fold universe. The current work offers insight into previously hypothesized evolutionary relationships between disparate folds and provides a resource for exploring novel ones. Proteins 2016; 84:1859-1874. © 2016 Wiley Periodicals, Inc.

Creation of active TIM barrel enzymes through genetic fusion of half-barrel domain constructs derived from two distantly related glycosyl hydrolases

Diverse unrelated enzymes that adopt the beta/alpha (or TIM) barrel topology display similar arrangements of beta/alpha units placed in a radial eight-fold symmetry around the barrel's axis. The TIM barrel was originally thought to be a single structural domain; however, it is now thought that TIM barrels arose from duplication and fusion of smaller half-barrels consisting of four beta/alpha units. We describe here the design, expression and purification, as well as characterization of folding, activity and stability, of chimeras of two TIM barrel glycosyl hydrolases, made by fusing different half-barrel domains derived from an endoglucanase from Clostridium cellulolyticum, CelCCA and a beta-glucosidase from Pyrococcus furiosus, CelB. We show that after refolding following purification from inclusion bodies, the two half-barrel fusion chimeras (CelCCACelB and CelBCelCCA) display catalytic activity although they assemble into large soluble oligomeric aggregated species containing chains of mixed beta and alpha structure. CelBCelCCA displays hyperthermophile-like structural stability as well as significant stability to chemical denaturation (C_m of 2.6 m guanidinium hydrochloride), whereas CelCCACelB displays mesophile-like stability (T_m of ~ 71 °C). The endoglucanase activities of both chimeras are an order of magnitude lower than those of CelB or CelCCA, whereas the beta-glucosidase activity of CelBCelCCA is about two orders of magnitude lower than that of CelB. The chimera CelCCACelB shows no beta-glucosidase activity. Our results demonstrate that half-barrel domains from unrelated sources can fold, assemble and function, with scope for improvement.

ENZYME

Pyrococcus furiosus beta-glucosidase (CelB, EC: 3.2.1.21). Clostridium cellulolyticum endoglucanase A (CelCCA, EC: 3.2.1.4).

Recurring sequence-structure motifs in (βα)8-barrel proteins and experimental optimization of a chimeric protein designed based on such motifs

An interesting way of generating novel artificial proteins is to combine sequence motifs from natural proteins, mimicking the evolutionary path suggested by natural proteins comprising recurring motifs. We analyzed the βα and αβ modules of TIM barrel proteins by structure alignment-based sequence clustering. A number of preferred motifs were identified. A chimeric TIM was designed by using recurring elements as mutually compatible interfaces. The foldability of the designed TIM protein was then significantly improved by six rounds of directed evolution. The melting temperature has been improved by more than 20°C. A variety of characteristics suggested that the resulting protein is well-folded. Our analysis provided a library of peptide motifs that is potentially useful for different protein engineering studies. The protein engineering strategy of using recurring motifs as interfaces to connect partial natural proteins may be applied to other protein folds.

A subdomain swap strategy for reengineering nonribosomal peptides

Nonribosomal peptide synthetases (NRPSs) protect microorganisms from environmental threats by producing diverse siderophores, antibiotics, and other peptide natural products. Their modular molecular structure is also attractive from the standpoint of biosynthetic engineering. Here we evaluate a methodology for swapping module specificities of these mega-enzymes that takes advantage of flavodoxin-like subdomains involved in substrate recognition. Nine subdomains encoding diverse specificities were transplanted into the Phe-specific GrsA initiation module of gramicidin S synthetase. All chimeras could be purified as soluble protein. One construct based on a Val-specific subdomain showed sizable adenylation activity and functioned as a Val-Pro diketopiperazine synthetase upon addition of the proline-specific GrsB1 module. These results suggest that subdomain swapping could be a viable alternative to previous NRPS design approaches targeting binding pockets, domains, or entire modules. The short length of the swapped sequence stretch may facilitate straightforward exploitation of the wealth of existing NRPS modules for combinatorial biosynthesis.

The TIM Barrel Architecture Facilitated the Early Evolution of Protein-Mediated Metabolism

The triosephosphate isomerase (TIM) barrel protein fold is a structurally repetitive architecture that is present in approximately 10 % of all enzymes. It is generally assumed that this ubiquity in modern proteomes reflects an essential historical role in early protein-mediated metabolism. Here, we provide quantitative and comparative analyses to support several hypotheses about the early importance of the TIM barrel architecture. An information theoretical analysis of protein structures supports the hypothesis that the TIM barrel architecture could arise more easily by duplication and recombination compared to other mixed α/β structures. We show that TIM barrel enzymes corresponding to the most taxonomically broad superfamilies also have the broadest range of functions, often aided by metal and nucleotide-derived cofactors that are thought to reflect an earlier stage of metabolic evolution. By comparison to other putatively ancient protein architectures, we find that the functional diversity of TIM barrel proteins cannot be explained simply by their antiquity. Instead, the breadth of TIM barrel functions can be explained, in part, by the incorporation of a broad range of cofactors, a trend that does not appear to be shared by proteins in general. These results support the hypothesis that the simple and functionally general TIM barrel architecture may have arisen early in the evolution of protein biosynthesis and provided an ideal scaffold to facilitate the metabolic transition from ribozymes, peptides, and geochemical catalysts to modern protein enzymes.

Evolution of a protein folding nucleus

The folding nucleus (FN) is a cryptic element within protein primary structure that enables an efficient folding pathway and is the postulated heritable element in the evolution of protein architecture; however, almost nothing is known regarding how the FN structurally changes as complex protein architecture evolves from simpler peptide motifs. We report characterization of the FN of a designed purely symmetric β-trefoil protein by ϕ-value analysis. We compare the structure and folding properties of key foldable intermediates along the evolutionary trajectory of the β-trefoil. The results show structural acquisition of the FN during gene fusion events, incorporating novel turn structure created by gene fusion. Furthermore, the FN is adjusted by circular permutation in response to destabilizing functional mutation. FN plasticity by way of circular permutation is made possible by the intrinsic C3 cyclic symmetry of the β-trefoil architecture, identifying a possible selective advantage that helps explain the prevalence of cyclic structural symmetry in the proteome.

How the Sequence of a Gene Specifies Structural Symmetry in Proteins

Internal symmetry is commonly observed in the majority of fundamental protein folds. Meanwhile, sufficient evidence suggests that nascent polypeptide chains of proteins have the potential to start the co-translational folding process and this process allows mRNA to contain additional information on protein structure. In this paper, we study the relationship between gene sequences and protein structures from the viewpoint of symmetry to explore how gene sequences code for structural symmetry in proteins. We found that, for a set of two-fold symmetric proteins from left-handed beta-helix fold, intragenic symmetry always exists in their corresponding gene sequences. Meanwhile, codon usage bias and local mRNA structure might be involved in modulating translation speed for the formation of structural symmetry: a major decrease of local codon usage bias in the middle of the codon sequence can be identified as a common feature; and major or consecutive decreases in local mRNA folding energy near the boundaries of the symmetric substructures can also be observed. The results suggest that gene duplication and fusion may be an evolutionarily conserved process for this protein fold. In addition, the usage of rare codons and the formation of higher order of secondary structure near the boundaries of symmetric substructures might have coevolved as conserved mechanisms to slow down translation elongation and to facilitate effective folding of symmetric substructures. These findings provide valuable insights into our understanding of the mechanisms of translation and its evolution, as well as the design of proteins via symmetric modules.

Case study on the evolution of hetero-oligomer interfaces based on the differences in paralogous proteins

We addressed the evolutionary trace of hetero-oligomer interfaces by comparing the structures of paralogous proteins; one of them is a monomer or homo-oligomer and the other is a hetero-oligomer. We found different trends in amino acid conservation pattern and hydrophobicity between homo-oligomer and hetero-oligomer. The degree of amino acid conservation in the interface of homo-oligomer has no obvious difference from that in the surface, whereas the degree of conservation is much higher in the interface of hetero-oligomer. The interface of homo-oligomer has a few very conserved residue positions, whereas the residue conservation in the interface of hetero-oligomer tends to be higher. In addition, the interface of hetero-oligomer has a tendency of being more hydrophobic compared with the one in homo-oligomer. We conjecture that these differences are related to the inherent symmetry in homo-oligomers that cannot exist in hetero-oligomers. Paucity of the structural data precludes statistical tests of these tendencies, yet the trend can be applied to the prediction of the interface of hetero-oligomer. We obtained putative interfaces of the subunits in CPSF (cleavage and polyadenylation specificity factor), one of the human pre-mRNA 3′-processing complexes. The locations of predicted interface residues were consistent with the known experimental data.

Tracing the structural evolution of eukaryotic ATP binding cassette transporter superfamily

The ATP binding cassette (ABC) transporters superfamily is one of the largest classes of membrane proteins. The core of the ABC transporter protein is composed of transmembrane domains (TMDs) and nucleotide binding domains (NBD). Eukaryotes ABC transporters are classified into seven main families (ABCA to ABCG) based on sequence similarity and domain organizations. With different domain number and domain organizations, eukaryote ABC transporters show diverse structures: the single structure (NBD or TMD), the ABC2 structure (NBD-NBD), the half structure (TMD-NBD or NBD-TMD) and the full structure (TMD-NBD-TMD-NBD or NBD-TMD-NBD-TMD). However, studies on how various ABC transporter gene structures evolved is still absent. Therefore, in this study, we comprehensively investigated the structural evolution of eukaryotic ABC transporters. The seven eukaryote ABC transporter families (A to G) fell into three groups: A&G group, B,C&D group and E&F group. There were at least four times the number of NBD and TMD fusion events in the origin of the half structure transporter. Two fusion modes were found in the full and ABC2 structure origination. Based on these findings, we present a putative structural evolutionary path of eukaryote ABC transporters that will increase our understanding on their origin, divergence and function.

Two independently folding units of Plasmodium profilin suggest evolution via gene fusion

Gene fusion is a common mechanism of protein evolution that has mainly been discussed in the context of multidomain or symmetric proteins. Less is known about fusion of ancestral genes to produce small single-domain proteins. Here, we show with a domain-swapped mutant Plasmodium profilin that this small, globular, apparently single-domain protein consists of two foldons. The separation of binding sites for different protein ligands in the two halves suggests evolution via an ancient gene fusion event, analogous to the formation of multidomain proteins. Finally, the two fragments can be assembled together after expression as two separate gene products. The possibility to engineer both domain-swapped dimers and half-profilins that can be assembled back to a full profilin provides perspectives for engineering of novel protein folds, e.g., with different scaffolding functions.

Biosynthesis of Histidine

The biosynthesis of histidine in Escherichia coli and Salmonella typhimurium has been an important model system for the study of relationships between the flow of intermediates through a biosynthetic pathway and the control of the genes encoding the enzymes that catalyze the steps in a pathway. This article provides a comprehensive review of the histidine biosynthetic pathway and enzymes, including regulation of the flow of intermediates through the pathway and mechanisms that regulate the amounts of the histidine biosynthetic enzymes. In addition, this article reviews the structure and regulation of the histidine (his) biosynthetic operon, including transcript processing, Rho-factor-dependent "classical" polarity, and the current model of his operon attenuation control. Emphasis is placed on areas of recent progress. Notably, most of the enzymes that catalyze histidine biosynthesis have recently been crystallized, and their structures have been determined. Many of the histidine biosynthetic intermediates are unstable, and the histidine biosynthetic enzymes catalyze some chemically unusual reactions. Therefore, these studies have led to considerable mechanistic insight into the pathway itself and have provided deep biochemical understanding of several fundamental processes, such as feedback control, allosteric interactions, and metabolite channeling. Considerable recent progress has also been made on aspects of his operon regulation, including the mechanism of pp(p)Gpp stimulation of his operon transcription, the molecular basis for transcriptional pausing by RNA polymerase, and pathway evolution. The progress in these areas will continue as sophisticated new genomic, metabolomic, proteomic, and structural approaches converge in studies of the histidine biosynthetic pathway and mechanisms of control of his biosynthetic genes in other bacterial species.

Half-Barrels Derived from a (β/α)8 Barrel β-Glycosidase Undergo an Activation Process

The evolution of (β/α)₈ barrel proteins is currently thought to have involved the fusion of two (β/α)₄ half-barrels, thereby conferring stability on the protein structure. After the formation of a whole (β/α)₈ barrel, this structure could evolve and diverge to form fully active enzymes. Interestingly, we show here that isolated (β/α)₄ half-barrels derived from the N- and C-terminal domains of the β-glucosidase Sfβgly (Sfβgly-N: residues 1 to 265; Sfβgly-C: residues 266 to 509) undergo an activation process, which renders them catalytically active. The rate constants of the activation process were calculated to be 0.029 and 0.032 h^-1 for Sfβgly-N and Sfβgly-C, respectively. Moreover, the Sfβgly-N and Sfβgly-C activation processes were simultaneous with modifications in their initial structure, which reduced the exposure of their tryptophan residues. Importantly, this activation was also coincident with an increase in the sizes of Sfβgly-N and Sfβgly-C particles. These novel observations suggest that the change in catalytic activity associated with the transition from a half to whole (β/α)₈ barrel might also have driven such an evolutionary process.