Metabolic evolution of pyranopterin-dependent biochemistry

Molybdenum (Mo)-dependent biochemistry is essential for many key metabolic pathways. However, theory and geological evidence suggests that its solubility during long intervals with low dioxygen would have limited its availability on early Earth. We developed models of metabolic evolution and found that reactions employing tungsten (W)-dependent biochemistry likely preceded Mo-dependent reactions, where Mo-usage increased dramatically after the production of dioxygen. Consistent with this finding, we analyzed genomes from over 65,000 phylogenetically diverse microbes and metagenomes from an environmental dataset, and we observed that dioxygen-utilizing prokaryotes living in aerobic niches are enriched with Mo-dependent enzymes as compared to anaerobic microbes. As an independent evaluation of this hypothesis, we combined protein language models, machine learning, and phylogenomic analysis to build a classifier for W- or Mo-pterin dependence in the DMSO reductase superfamily, and we found that W-pterin-dependent enzymes cluster near the root of the tree and that a subset of late-evolving aldehyde oxidoreductases (AORs) from aerobes are predicted to rely on Mo instead of W. Overall, our combination of metabolic modeling, phenotypic analysis, machine learning, and phylogenomic analysis suggest that Mo-pterin-dependent biochemistry likely derived from W-pterin-dependent biochemistry, and that Mo-usage increased drastically after the rise of oxygen.