Characterizing selective pressures on the pathway for de novo biosynthesis of pyrimidines in yeast

PEMPS: a phylogenetic software tool to model the evolution of metabolic pathways

BACKGROUND

Metabolic pathways support the enzyme flux that converts input chemicals into energy and cellular building blocks. With a constant rate of input, steady-state flux is achieved when metabolite concentrations and reaction rates remain constant over time. Individual genes undergo mutation, while selection acts on higher level functions of the pathway, such as steady-state flux where applicable. Modeling the evolution of metabolic pathways through mechanistic sets of ordinary differential equations is a piece of the genotype-phenotype map model for interpreting genetic variation and inter-specific differences. Such models can generate distinct compensatory changes and adaptive changes from directional selection, indicating single nucleotide polymorphisms and fixed differences that could affect phenotype. If used for inference, this would ultimately enable detection of selection on metabolic pathways as well as inference of ancestral states for metabolic pathway function.

RESULTS

A software tool for simulating the evolution of metabolic pathways based upon underlying biochemistry, phylogenetics, and evolutionary considerations is presented. The Python program, Phylogenetic Evolution of Metabolic Pathway Simulator (PEMPS), implements a mutation-selection framework to simulate the evolution of the pathway over a phylogeny by interfacing with COPASI to calculate the steady-state flux of the metabolic network, introducing mutations as alterations in parameter values according to a model, and calculating a fitness score and corresponding probability of fixation based on the change in steady-state flux value(s). Results from simulations are consistent with a priori expectations of fixation probabilities and systematic change in model parameters.

CONCLUSIONS

The PEMPS program simulates the evolution of a metabolic pathway with a mutation-selection modeling framework based on criteria like steady-state flux that is designed to work with SBML-formatted kinetic models, and Newick-formatted phylogenetic trees. The Python software is run on the Linux command line and is available at https://github.com/nmccloskey/PEMPS .

Computer-Based Design of a Cell Factory for High-Yield Cytidine Production

Pyrimidine ribonucleotide de novo biosynthesis pathway (PRdnBP) is an important pathway to produce pyrimidine nucleosides. We attempted to systematically investigate PRdnBP in Escherichia coli with genome-scale metabolic models and utilized the models to guide strain design. The balance of central carbon metabolism and PRdnBP affected the production of cytidine from glucose. Using Bayesian metabolic flux analysis, the effect of modified PRdnBP on the metabolic network was analyzed. The acetate overflow became coupled with PRdnBP flux, while they were originally independent under oxygen-sufficient conditions. The coupling between cytidine production and acetate secretion in the modified strain was weakened by arcA deletion, which resulted in further improving the efficient accumulation of cytidine. In total, 1.28 g/L of cytidine with a yield of 0.26 g/g glucose was produced. The yield of cytidine produced by E. coli is higher than previous reports. Our strategy provides an effective attempt to find metabolic bottlenecks in genetically engineered bacteria by using flux coupling analysis.

New Insights into the Interaction of Class II Dihydroorotate Dehydrogenases with Ubiquinone in Lipid Bilayers as a Function of Lipid Composition

The fourth enzymatic reaction in the de novo pyrimidine biosynthesis, the oxidation of dihydroorotate to orotate, is catalyzed by dihydroorotate dehydrogenase (DHODH). Enzymes belonging to the DHODH Class II are membrane-bound proteins that use ubiquinones as their electron acceptors. We have designed this study to understand the interaction of an N-terminally truncated human DHODH (HsΔ29DHODH) and the DHODH from Escherichia coli (EcDHODH) with ubiquinone (Q₁₀) in supported lipid membranes using neutron reflectometry (NR). NR has allowed us to determine in situ, under solution conditions, how the enzymes bind to lipid membranes and to unambiguously resolve the location of Q₁₀. Q₁₀ is exclusively located at the center of all of the lipid bilayers investigated, and upon binding, both of the DHODHs penetrate into the hydrophobic region of the outer lipid leaflet towards the Q₁₀. We therefore show that the interaction between the soluble enzymes and the membrane-embedded Q₁₀ is mediated by enzyme penetration. We can also show that EcDHODH binds more efficiently to the surface of simple bilayers consisting of 1-palmitoyl, 2-oleoyl phosphatidylcholine, and tetraoleoyl cardiolipin than HsΔ29DHODH, but does not penetrate into the lipids to the same degree. Our results also highlight the importance of Q₁₀, as well as lipid composition, on enzyme binding.

Protein production, kinetic and biophysical characterization of three human dihydroorotate dehydrogenase mutants associated with Miller syndrome

Miller syndrome is a rare Mendelian disorder caused by mutations in the gene encoding human dihydroorotate dehydrogenase (DHODH). Human DHODH, a Class II DHODH, is an integral protein of the inner mitochondrial membrane (IMM) catalyzing the fourth step of the de novo pyrimidine biosynthesis pathway. Here we present a summary of the state of knowledge regarding Miller syndrome in the absence of any current review on the topic. We then describe the production and characterization of three distinct DHODH missense mutations (G19E, E52G, R135C) associated with Miller syndrome by means of enzyme kinetics and biophysical techniques. These human DHODH mutants were produced both in E. coli and in insect cells using the baculovirus expression vector system. We can show that the effects of these mutations differ from each other and the wild-type enzyme with respect to decreased enzymatic activity, decreased protein stability and probably disturbance of the correct import into the IMM. In addition, our results show that the N-terminus of human DHODH is not only a structural element necessary for correct mitochondrial import and location of DHODH on the outer side of the IMM, but also influences thermal stability, enzymatic activity and affects the kinetic parameters.Supplemental data for this article is available online at https://doi.org/10.1080/15257770.2021.2023749 .

Characterizing lineage-specific evolution and the processes driving genomic diversification in chordates

Background

Understanding the origins of genome content has long been a goal of molecular evolution and comparative genomics. By examining genome evolution through the guise of lineage-specific evolution, it is possible to make inferences about the evolutionary events that have given rise to species-specific diversification. Here we characterize the evolutionary trends found in chordate species using The Adaptive Evolution Database (TAED). TAED is a database of phylogenetically indexed gene families designed to detect episodes of directional or diversifying selection across chordates. Gene families within the database have been assessed for lineage-specific estimates of dN/dS and have been reconciled to the chordate species to identify retained duplicates. Gene families have also been mapped to the functional pathways and amino acid changes which occurred on high dN/dS lineages have been mapped to protein structures.

Results

An analysis of this exhaustive database has enabled a characterization of the processes of lineage-specific diversification in chordates. A pathway level enrichment analysis of TAED determined that pathways most commonly found to have elevated rates of evolution included those involved in metabolism, immunity, and cell signaling. An analysis of protein fold presence on proteins, after normalizing for frequency in the database, found common folds such as Rossmann folds, Jelly Roll folds, and TIM barrels were overrepresented on proteins most likely to undergo directional selection. A set of gene families which experience increased numbers of duplications within short evolutionary times are associated with pathways involved in metabolism, olfactory reception, and signaling. An analysis of protein secondary structure indicated more relaxed constraint in β-sheets and stronger constraint on alpha Helices, amidst a general preference for substitutions at exposed sites. Lastly a detailed analysis of the ornithine decarboxylase gene family, a key enzyme in the pathway for polyamine synthesis, revealed lineage-specific evolution along the lineage leading to Cetacea through rapid sequence evolution in a duplicate gene with amino acid substitutions causing active site rearrangement.

Conclusion

Episodes of lineage-specific evolution are frequent throughout chordate species. Both duplication and directional selection have played large roles in the evolution of the phylum. TAED is a powerful tool for facilitating this understanding of lineage-specific evolution.

Metabolic Determinants of Enzyme Evolution in a Genome-Scale Bacterial Metabolic Network

Different genes and proteins evolve at very different rates. To identify the factors that explain these differences is an important aspect of research in molecular evolution. One such factor is the role a protein plays in a large molecular network. Here, we analyze the evolutionary rates of enzyme-coding genes in the genome-scale metabolic network of Escherichia coli to find the evolutionary constraints imposed by the structure and function of this complex metabolic system. Central and highly connected enzymes appear to evolve more slowly than less connected enzymes, but we find that they do so as a by-product of their high abundance, and not because of their position in the metabolic network. In contrast, enzymes catalyzing reactions with high metabolic flux-high substrate to product conversion rates-evolve slowly even after we account for their abundance. Moreover, enzymes catalyzing reactions that are difficult to by-pass through alternative pathways, such that they are essential in many different genetic backgrounds, also evolve more slowly. Our analyses show that an enzyme's role in the function of a metabolic network affects its evolution more than its place in the network's structure. They highlight the value of a system-level perspective for studies of molecular evolution.

The Adaptive Evolution Database (TAED): A New Release of a Database of Phylogenetically Indexed Gene Families from Chordates

With the large collections of gene and genome sequences, there is a need to generate curated comparative genomic databases that enable interpretation of results in an evolutionary context. Such resources can facilitate an understanding of the co-evolution of genes in the context of a genome mapped onto a phylogeny, of a protein structure, and of interactions within a pathway. A phylogenetically indexed gene family database, the adaptive evolution database (TAED), is presented that organizes gene families and their evolutionary histories in a species tree context. Gene families include alignments, phylogenetic trees, lineage-specific dN/dS ratios, reconciliation with the species tree to enable both the mapping and the identification of duplication events, mapping of gene families onto pathways, and mapping of amino acid substitutions onto protein structures. In addition to organization of the data, new phylogenetic visualization tools have been developed to aid in interpreting the data that are also available, including TreeThrasher and TAED Tree Viewer. A new resource of gene families organized by species and taxonomic lineage promises to be a valuable comparative genomics database for molecular biologists, evolutionary biologists, and ecologists. The new visualization tools and database framework will be of interest to both evolutionary biologists and bioinformaticians.

Characterizing the roles of changing population size and selection on the evolution of flux control in metabolic pathways

Background

Understanding the genotype-phenotype map is fundamental to our understanding of genomes. Genes do not function independently, but rather as part of networks or pathways. In the case of metabolic pathways, flux through the pathway is an important next layer of biological organization up from the individual gene or protein. Flux control in metabolic pathways, reflecting the importance of mutation to individual enzyme genes, may be evolutionarily variable due to the role of mutation-selection-drift balance. The evolutionary stability of rate limiting steps and the patterns of inter-molecular co-evolution were evaluated in a simulated pathway with a system out of equilibrium due to fluctuating selection, population size, or positive directional selection, to contrast with those under stabilizing selection.

Results

Depending upon the underlying population genetic regime, fluctuating population size was found to increase the evolutionary stability of rate limiting steps in some scenarios. This result was linked to patterns of local adaptation of the population. Further, during positive directional selection, as with more complex mutational scenarios, an increase in the observation of inter-molecular co-evolution was observed.

Conclusions

Differences in patterns of evolution when systems are in and out of equilibrium, including during positive directional selection may lead to predictable differences in observed patterns for divergent evolutionary scenarios. In particular, this result might be harnessed to detect differences between compensatory processes and directional processes at the pathway level based upon evolutionary observations in individual proteins. Detecting functional shifts in pathways reflects an important milestone in predicting when changes in genotypes result in changes in phenotypes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0962-7) contains supplementary material, which is available to authorized users.

Selection on metabolic pathway function in the presence of mutation-selection-drift balance leads to rate-limiting steps that are not evolutionarily stable

Background

While commonly assumed in the biochemistry community that the control of metabolic pathways is thought to be critical to cellular function, it is unclear if metabolic pathways generally have evolutionarily stable rate limiting (flux controlling) steps.

Results

A set of evolutionary simulations using a kinetic model of a metabolic pathway was performed under different conditions to evaluate the evolutionary stability of rate limiting steps. Simulations used combinations of selection for steady state flux, selection against the cost of molecular biosynthesis, and selection against the accumulation of high concentrations of a deleterious intermediate. Two mutational regimes were used, one with mutations that on average were neutral to molecular phenotype and a second with a preponderance of activity-destroying mutations. The evolutionary stability of rate limiting steps was low in all simulations with non-neutral mutational processes. Clustering of parameter co-evolution showed divergent inter-molecular evolutionary patterns under different evolutionary regimes.

Conclusions

This study provides a null model for pathway evolution when compensatory processes dominate with potential applications to predicting pathway functional change. This result also suggests a possible mechanism in which studies in statistical genetics that aim to associate a genotype to a phenotype assuming independent action of variants may be mis-specified through a mis-characterization of the link between individual gene function and pathway function. A better understanding of the genotype-phenotype map has potential applications in differentiating between compensatory changes and directional selection on pathways as well as detecting SNPs and fixed differences that might have phenotypic effects.

Reviewers

This article was reviewed by Arne Elofsson, David Ardell, and Shamil Sunyaev.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-016-0133-6) contains supplementary material, which is available to authorized users.

Selection on Network Dynamics Drives Differential Rates of Protein Domain Evolution

The long-held principle that functionally important proteins evolve slowly has recently been challenged by studies in mice and yeast showing that the severity of a protein knockout only weakly predicts that protein's rate of evolution. However, the relevance of these studies to evolutionary changes within proteins is unknown, because amino acid substitutions, unlike knockouts, often only slightly perturb protein activity. To quantify the phenotypic effect of small biochemical perturbations, we developed an approach to use computational systems biology models to measure the influence of individual reaction rate constants on network dynamics. We show that this dynamical influence is predictive of protein domain evolutionary rate within networks in vertebrates and yeast, even after controlling for expression level and breadth, network topology, and knockout effect. Thus, our results not only demonstrate the importance of protein domain function in determining evolutionary rate, but also the power of systems biology modeling to uncover unanticipated evolutionary forces.

Flux Control in Glycolysis Varies Across the Tree of Life

Biochemical thought posits that rate-limiting steps (defined here as points of flux control) are strongly selected as points of pathway regulation and control and are thus expected to be evolutionarily conserved. Conversely, population genetic thought based upon the concepts of mutation-selection-drift balance at the pathway level might suggest variation in flux controlling steps over evolutionary time. Glycolysis, as one of the most conserved and best characterized pathways, was studied to evaluate its evolutionary conservation. The flux controlling step in glycolysis was found to vary over the tree of life. Further, phylogenetic analysis suggested at least 60 events of gene duplication and additional events of putative positive selection that might alter pathway kinetic properties. Together, these results suggest that even with presumed largely negative selection on pathway output on glycolysis, the co-evolutionary process under the hood is dynamic.