A thermodynamic atlas of carbon redox chemical space

Computational Prediction of Speciation Diagrams and Nucleation Mechanisms: Molecular Vanadium, Niobium, and Tantalum Oxide Nanoclusters in Solution

Understanding the aqueous speciation of molecular metal-oxo-clusters plays a key role in different fields such as catalysis, electrochemistry, nuclear waste recycling, and biochemistry. To describe the speciation accurately, it is essential to elucidate the underlying self-assembly processes. Herein, we apply a computational method to predict the speciation and formation mechanisms of polyoxovanadates, -niobates, and -tantalates. While polyoxovanadates have been widely studied, polyoxoniobates and -tantalates lack the same level of understanding. First, we propose a pentavanadate cluster ([V₅O₁₄]^3-) as a key intermediate for the formation of the decavanadate. Our computed phase speciation diagram is in particularly good agreement with the experiments. Second, we report the formation constants of the heptaniobate, [Nb₇O₂₂]^9-, decaniobate, [Nb₁₀O₂₈]^6-, and tetracosaniobate [H₉Nb₂₄O₇₂]^15-. Additionally, we compute the speciation and phase diagram of niobium, which so far was restricted to Lindqvist derivates. Finally, we predict the formation constant of the decatantalate ([Ta₁₀O₂₆]^6-) in water, even though it had only been synthesized in toluene. Furthermore, we also calculate the corresponding speciation and phase diagrams for polyoxotantalates. Overall, we show that our method can be successfully applied to different families of molecular metal oxides without any need for readjustments; therefore, it can be regarded as a trustworthy tool for exploring polyoxometalates' chemistry.

Stoichiometric ratios for biotics and xenobiotics capture effective metabolic coupling to re(de)fine biodegradation

Preserving human and environmental health requires anthropogenic pollutants to be biologically degradable. Depending on concentration, both nutrients and pollutants induce and activate metabolic capacity in the endemic bacterial consortium, which in turn aids their degradation. Knowledge on such 'acclimation' is rarely implemented in risk assessment cost-effectively. As a result, an accurate description of the mechanisms and kinetics of biodegradation remains problematic. In this study, we defined a yield 'effectivity', comprising the effectiveness at which a pollutant (substrate) enhances its own degradation by inducing (biomass) cofactors involved therein. Our architecture for calculation represents the interplay between concentration and metabolism via both stoichiometric and thermodynamic concepts. The calculus for yield 'effectivity' is biochemically intuitive, implicitly embeds co-metabolism and distinguishes 'endogenic' from 'exogenic' substances' reflecting various phenomena in biodegradation and bio-transformation studies. We combined data on half-lives of pollutants/nutrients in wastewater and surface water with transition-state rate theory to obtain also experimental values for effective yields. These quantify the state of acclimation: the portion of biodegradation kinetics attributable to (contributed by) 'natural metabolism', in view of similarity to natural substances. Calculated and experimental values showed statistically significant correspondence. Particularly, carbohydrate metabolism and nucleic acid metabolism appeared relevant for acclimation (R² = 0.11-0.42), affecting rates up to 10^4.9(±0.7) times: under steady-state acclimation, a compound stoichiometrically identical to carbohydrates or nucleic acids, is 10^3.2 to 10^4.9 times faster aerobically degraded than a compound marginally similar. Our new method, simulating (contribution by) the state of acclimation, supplements existing structure-biodegradation and kinetic models for predicting biodegradation in wastewater and surface water. The accuracy of prediction may increase when characterizing nutrients/co-metabolites in terms of, e.g., elemental analysis. We discuss strengths and limitations of our approach by comparison to empirical and mechanism-based methods.

Environmental conditions drive self-organization of reaction pathways in a prebiotic reaction network

The evolution of life from the prebiotic environment required a gradual process of chemical evolution towards greater molecular complexity. Elaborate prebiotically relevant synthetic routes to the building blocks of life have been established. However, it is still unclear how functional chemical systems evolved with direction using only the interaction between inherent molecular chemical reactivity and the abiotic environment. Here we demonstrate how complex systems of chemical reactions exhibit well-defined self-organization in response to varying environmental conditions. This self-organization allows the compositional complexity of the reaction products to be controlled as a function of factors such as feedstock and catalyst availability. We observe how Breslow's cycle contributes to the reaction composition by feeding C₂ building blocks into the network, alongside reaction pathways dominated by formaldehyde-driven chain growth. The emergence of organized systems of chemical reactions in response to changes in the environment offers a potential mechanism for a chemical evolution process that bridges the gap between prebiotic chemical building blocks and the origin of life.

Protein cost minimization promotes the emergence of coenzyme redundancy

Metabolism relies on a small class of molecules (coenzymes) that serve as universal donors and acceptors of key chemical groups and electrons. Although metabolic networks crucially depend on structurally redundant coenzymes [e.g., NAD(H) and NADP(H)] associated with different enzymes, the criteria that led to the emergence of this redundancy remain poorly understood. Our combination of modeling and structural and sequence analysis indicates that coenzyme redundancy may not be essential for metabolism but could rather constitute an evolved strategy promoting efficient usage of enzymes when biochemical reactions are near equilibrium. Our work suggests that early metabolism may have operated with fewer coenzymes and that adaptation for metabolic efficiency may have driven the rise of coenzyme diversity in living systems.

Coenzymes distribute a variety of chemical moieties throughout cellular metabolism, participating in group (e.g., phosphate and acyl) and electron transfer. For a variety of reactions requiring acceptors or donors of specific resources, there often exist degenerate sets of molecules [e.g., NAD(H) and NADP(H)] that carry out similar functions. Although the physiological roles of various coenzyme systems are well established, it is unclear what selective pressures may have driven the emergence of coenzyme redundancy. Here, we use genome-wide metabolic modeling approaches to decompose the selective pressures driving enzymatic specificity for either NAD(H) or NADP(H) in the metabolic network of Escherichia coli. We found that few enzymes are thermodynamically constrained to using a single coenzyme, and in principle a metabolic network relying on only NAD(H) is feasible. However, structural and sequence analyses revealed widespread conservation of residues that retain selectivity for either NAD(H) or NADP(H), suggesting that additional forces may shape specificity. Using a model accounting for the cost of oxidoreductase enzyme expression, we found that coenzyme redundancy universally reduces the minimal amount of protein required to catalyze coenzyme-coupled reactions, inducing individual reactions to strongly prefer one coenzyme over another when reactions are near thermodynamic equilibrium. We propose that protein minimization generically promotes coenzyme redundancy and that coenzymes typically thought to exist in a single pool (e.g., coenzyme A [CoA]) may exist in more than one form (e.g., dephospho-CoA).

Oxidative metabolisms catalyzed Earth’s oxygenation

The burial of organic carbon, which prevents its remineralization via oxygen-consuming processes, is considered one of the causes of Earth’s oxygenation. Yet, higher levels of oxygen are thought to inhibit burial. Here we propose a resolution of this conundrum, wherein Earth’s initial oxygenation is favored by oxidative metabolisms generating partially oxidized organic matter (POOM), increasing burial via interaction with minerals in sediments. First, we introduce the POOM hypothesis via a mathematical argument. Second, we reconstruct the evolutionary history of one key enzyme family, flavin-dependent Baeyer–Villiger monooxygenases, that generates POOM, and show the temporal consistency of its diversification with the Proterozoic and Phanerozoic atmospheric oxygenation. Finally, we propose that the expansion of oxidative metabolisms instigated a positive feedback, which was amplified by the chemical changes to minerals on Earth’s surface. Collectively, these results suggest that Earth’s oxygenation is an autocatalytic transition induced by a combination of biological innovations and geological changes.

How Earth’s atmosphere became oxygenated remains enigmatic. Here the authors use mathematical and phylogenetic analyses to find that Earth’s oxygenation is induced by the interactions of microbial oxidative metabolites with sediment minerals.

Unlocking Phase Diagrams for Molybdenum and Tungsten Nanoclusters and Prediction of their Formation Constants

Understanding and controlling aqueous speciation of metal oxides are key for the discovery and development of novel materials, and challenge both experimental and computational approaches. Here we present a computational method, called POMSimulator, which is able to predict speciation phase diagrams (Conc. vs pH) for multispecies chemical equilibria in solution, and which we apply to molybdenum and tungsten isopolyoxoanions (IPAs). Starting from the MO₄ monomers, and considering dimers, trimers, and larger species, the chemical reaction networks involved in the formation of [H₃₂Mo₃₆O₁₂₈]^8- and [W₁₂O₄₂]^12- are sampled in an automatic manner. This information is used for setting up ∼10⁵ speciation models, and from there, we generate the speciation phase diagrams, which show an insightful picture of the behavior of IPAs in aqueous solution. Furthermore, we predict the values of 107 formation constants for a diversity of molybdenum and tungsten molecular oxides. Among these species, we could include several pentagonal-shaped species and very reactive tungsten intermediates as well. Last but not least, the calibration employed for correcting the density functional theory (DFT) Gibbs energies is remarkably similar for both metals, which suggests that a general rule might exist for correcting computed free energies for other metals.