Search

Metal-organic cages form well-defined microenvironments that can enhance the catalytic proficiency of encapsulated transition metal complexes (TMCs). We introduce a screening protocol to efficiently identify TMCs that are promising candidates for encapsulation in the Ga 4 L 6 12- nanocage. We obtain TMCs from the Cambridge Structural Database with geometric and electronic characteristics amenable to encapsulation and mine the text of associated manuscripts to curate TMCs with documented catalytic functionality. By docking candidate TMCs inside the nanocage cavity and carrying out electronic structure calculations, we identify a subset of successfully optimized candidates (TMC-34) and observe that encapsulated guests occupy an average of 60% of the cavity volume, in line with previous observations. Notably, some guests occupy as much as 72% of the cavity as a result of linker rotation. Encapsulation has a universal effect on the electrostatic potential (ESP), systematically decreasing the ESP at the metal center of each TMC in the TMC-34 data set, while minimally altering TMC metal partial charges. Collectively these observations support geometry-based screening of potential guests and suggest that encapsulation in Ga 4 L 6 12- cages could electrostatically stabilize diverse cationic or electropositive intermediates. We highlight candidate guests with associated known reactivity and solubility most amenable for encapsulation in experimental follow-up studies.

Metal-organic frameworks (MOFs) are porous materials with applications in gas separations and catalysis, but a lack of water stability often limits their practical use given the ubiquity of water. Consequently, it is useful to predict whether a MOF is water-stable before investing time and resources into synthesis. Existing heuristics for designing water-stable MOFs lack generality and limit the diversity of explored chemistry due to narrowly defined criteria. Machine learning (ML) models offer the promise to improve the generality of predictions but require data. In an improvement on previous efforts, we enlarge the available training data for MOF water stability prediction by over 400%, adding 911 MOFs with water stability labels assigned through semiautomated manuscript analysis to curate the new data set WS24. The additional data are shown to improve ML model performance (test ROC-AUC > 0.8) over diverse chemistry for the prediction of both water stability and stability in harsher acidic conditions. We illustrate how the expanded data set and models can be used with a previously developed activation stability model in combination with genetic algorithms to quickly screen ∼10,000 MOFs from a space of hundreds of thousands for candidates with multivariate stability (upon activation, in water, and in acid). We uncover metal- and geometry-specific design rules for robust MOFs. The data set and ML models developed in this work, which we disseminate through an easy-to-use web interface, are expected to contribute toward the accelerated discovery of novel, water-stable MOFs for applications such as direct air gas capture and water treatment.

Glycan-binding proteins, or lectins, recognize distinct structural elements of polysaccharides, to mediate myriad biological functions. Targeting glycan-binding proteins involved in human disease has been challenging due to an incomplete understanding of the molecular mechanisms that govern protein-glycan interactions. Bioinformatics and structural studies of glycan-binding proteins indicate that aromatic residues with the potential for CH-π interactions are prevalent in glycan-binding sites. However, the contributions of these CH-π interactions to glycan binding and their relevance in downstream function remain unclear. An emblematic lectin, human galectin-3, recognizes lactose and N -acetyllactosamine-containing glycans by positioning the electropositive face of a galactose residue over the tryptophan 181 (W181) indole forming a CH-π interaction. We generated a suite of galectin-3 W181 variants to assess the importance of these CH-π interactions to glycan binding and function. As determined experimentally and further validated with computational modeling, variants with smaller or less electron-rich aromatic side chains (W181Y, W181F, W181H) or sterically similar but nonaromatic residues (W181M, W181R) showed poor or undetectable binding to lactose and attenuated ability to bind mucins or agglutinate red blood cells. The latter functions depend on multivalent binding, highlighting that weakened CH-π interactions cannot be overcome by avidity. Two galectin-3 variants with disrupted hydrogen bonding interactions (H158A and E184A) showed similarly impaired lactose binding. Molecular simulations demonstrate that all variants have decreased binding orientation stability relative to native galectin-3. Thus, W181 collaborates with the endogenous hydrogen bonding network to enhance binding affinity for lactose, and abrogation of these CH-π interactions is as deleterious as eliminating key hydrogen bonding interactions. These findings underscore the critical roles of CH-π interactions in carbohydrate binding and lectin function and will aid the development of novel lectin inhibitors., Competing Interests: The authors declare no competing financial interest., (© 2024 The Authors. Published by American Chemical Society.)

Metal-organic frameworks (MOFs) are promising materials for gas sensing but are often limited to single-use detection. A hybridization strategy is demonstrated synergistically deploying conductive MOFs (cMOFs) and conductive polymers (cPs) as two complementary mixed ionic-electronic conductors in high-performing stand-alone chemiresistors. This work presents significant improvement in i) sensor recovery kinetics, ii) cycling stability, and iii) dynamic range at room temperature. The effect of hybridization across well-studied cMOFs is demonstrated based on 2,3,6,7,10,11-hexahydroxytriphenylene (HHTP) and 2,3,6,7,10,11-hexaiminotriphenylene (HITP) ligands with varied metal nodes (Co, Cu, Ni). A comprehensive mechanistic study is conducted to relate energy band alignments at the heterojunctions between the MOFs and the polymer with sensing thermodynamics and binding kinetics. The findings reveal that hole enrichment of the cMOF component upon hybridization leads to selective enhancement in desorption kinetics, enabling significantly improved sensor recovery at room temperature, and thus long-term response retention. This mechanism is further supported by density functional theory calculations on sorbate-analyte interactions. It is also found that alloying cPs and cMOFs enables facile thin film co-processing and device integration, potentially unlocking the use of these hybrid conductors in diverse electronic applications., (© 2024 The Authors. Advanced Materials published by Wiley‐VCH GmbH.)

The aza Paternò-Büchi reaction is a [2+2]-cycloaddition reaction between imines and alkenes that produces azetidines, four-membered nitrogen-containing heterocycles. Currently, successful examples rely primarily on either intramolecular variants or cyclic imine equivalents. To unlock the full synthetic potential of aza Paternò-Büchi reactions, it is essential to extend the reaction to acyclic imine equivalents. Here, we report that matching of the frontier molecular orbital energies of alkenes with those of acyclic oximes enables visible light-mediated aza Paternò-Büchi reactions through triplet energy transfer catalysis. The utility of this reaction is further showcased in the synthesis of epi- penaresidin B. Density functional theory computations reveal that a competition between the desired [2+2]-cycloaddition and alkene dimerization determines the success of the reaction. Frontier orbital energy matching between the reactive components lowers transition-state energy (Δ G ǂ ) values and ultimately promotes reactivity.

Thioesters are an essential functional group in biosynthetic pathways, which has motivated their development as reactive handles in probes and peptide assembly. Thioester exchange is typically accelerated by catalysts or elevated pH. Here, we report the use of bifunctional aromatic thioesters as dynamic covalent cross-links in hydrogels, demonstrating that at physiologic pH in aqueous conditions, transthioesterification facilitates stress relaxation on the time scale of hundreds of seconds. We show that intramolecular hydrogen bonding is responsible for accelerated exchange, evident in both molecular kinetics and macromolecular stress relaxation. Drawing from concepts in the vitrimer literature, this system exemplifies how dynamic cross-links that exchange through an associative mechanism enable tunable stress relaxation without altering stiffness.

Supramolecular polymer networks contain non-covalent cross-links that enable access to broadly tunable mechanical properties and stimuli-responsive behaviors; the incorporation of multiple unique non-covalent cross-links within such materials further expands their mechanical responses and functionality. To date, however, the design of such materials has been accomplished through discrete combinations of distinct interaction types in series, limiting materials design logic. Here we introduce the concept of leveraging "nested" supramolecular crosslinks, wherein two distinct types of non-covalent interactions exist in parallel, to control bulk material functions. To demonstrate this concept, we use polymer-linked Pd 2 L 4 metal-organic cage (polyMOC) gels that form hollow metal-organic cage junctions through metal-ligand coordination and can exhibit well-defined host-guest binding within their cavity. In these "nested" supramolecular network junctions, the thermodynamics of host-guest interactions within the junctions affect the metal-ligand interactions that form those junctions, ultimately translating to substantial guest-dependent changes in bulk material properties that could not be achieved in traditional supramolecular networks with multiple interactions in series., (© 2024. The Author(s).)

Second row elements in small- and medium-rings modulate strain. Herein we report the synthesis of two novel oligosilyl-containing cycloalkynes that exhibit angle-strain, as observed by X-ray crystallography. However, the angle-strained sila-cyclooctynes are sluggish participants in cycloadditions with benzyl azide. A distortion-interaction model analysis based on density functional theory calculations was performed.

The discovery of new compounds with pharmacological properties is usually a lengthy, laborious and expensive process. Thus, there is increasing interest in developing workflows that allow for the rapid synthesis and evaluation of libraries of compounds with the aim of identifying leads for further drug development. Herein, we apply combinatorial synthesis to build a library of 90 iridium(III) complexes (81 of which are new) over two synthesise-and-test cycles, with the aim of identifying potential agents for photodynamic therapy. We demonstrate the power of this approach by identifying highly active complexes that are well-tolerated in the dark but display very low nM phototoxicity against cancer cells. To build a detailed structure-activity relationship for this class of compounds we have used density functional theory (DFT) calculations to determine some key electronic parameters and study correlations with the experimental data. Finally, we present an optimised semi-automated synthesise-and-test protocol to obtain multiplex data within 72 hours., (© 2024 The Authors. Angewandte Chemie International Edition published by Wiley-VCH GmbH.)

While computational screening with density functional theory (DFT) is frequently employed for the screening of metal-organic frameworks (MOFs) for gas separation and storage, commonly applied generalized gradient approximations (GGAs) exhibit self-interaction errors, which hinder the predictions of adsorption energies. We investigate the Hubbard U parameter to augment DFT calculations for full periodic MOFs, targeting a more precise modeling of gas molecule-MOF interactions, specifically for N2, CO2, and O2. We introduce a calibration scheme for the U parameter, which is tailored for each MOF, by leveraging higher-level calculations on the secondary building unit (SBU) of the MOF. When applied to the full periodic MOF, the U parameter calibrated against hybrid HSE06 calculations of SBUs successfully reproduces hybrid-quality calculations of the adsorption energy of the periodic MOF. The mean absolute deviation of adsorption energies reduces from 0.13 eV for a standard GGA treatment to 0.06 eV with the calibrated U, demonstrating the utility of the calibration procedure when applied to the full MOF structure. Furthermore, attempting to use coupled cluster singles and doubles with perturbative triples calculations of isolated SBUs for this calibration procedure shows varying degrees of success in predicting the experimental heat of adsorption. It improves accuracy for N2 adsorption for cases of overbinding, whereas its impact on CO2 is minimal, and ambiguities in spin state assignment hinder consistent improvements of O2 adsorption. Our findings emphasize the limitations of cluster models and advocate the use of full periodic MOF systems with a calibrated U parameter, providing a more comprehensive understanding of gas adsorption in MOFs., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).)

Polymers that release small molecules in response to mechanical force are promising candidates as next-generation on-demand delivery systems. Despite advancements in the development of mechanophores for releasing diverse payloads through careful molecular design, the availability of scaffolds capable of discharging biomedically significant cargos in substantial quantities remains scarce. In this report, we detail a nonscissile mechanophore built from an 8-thiabicyclo[3.2.1]octane 8,8-dioxide ( TBO ) motif that releases one equivalent of sulfur dioxide (SO 2 ) from each repeat unit. The TBO mechanophore exhibits high thermal stability but is activated mechanochemically using solution ultrasonication in either organic solvent or aqueous media with up to 63% efficiency, equating to 206 molecules of SO 2 released per 143.3 kDa chain. We quantified the mechanochemical reactivity of TBO by single-molecule force spectroscopy and resolved its single-event activation. The force-coupled rate constant for TBO opening reaches ∼9.0 s -1 at ∼1520 pN, and each reaction of a single TBO domain releases a stored length of ∼0.68 nm. We investigated the mechanism of TBO activation using ab initio steered molecular dynamic simulations and rationalized the observed stereoselectivity. These comprehensive studies of the TBO mechanophore provide a mechanically coupled mechanism of multi-SO 2 release from one polymer chain, facilitating the translation of polymer mechanochemistry to potential biomedical applications.

Hydrogen fluoride (HF) is a versatile reagent for material transformation, with applications in self-immolative polymers, remodeled siloxanes, and degradable polymers. The responsive in situ generation of HF in materials therefore holds promise for new classes of adaptive material systems. Here, we report the mechanochemically coupled generation of HF from alkoxy- gem -difluorocyclopropane ( g DFC) mechanophores derived from the addition of difluorocarbene to enol ethers. Production of HF involves an initial mechanochemically assisted rearrangement of g DFC mechanophore to α-fluoro allyl ether whose regiochemistry involves preferential migration of fluoride to the alkoxy-substituted carbon, and ab initio steered molecular dynamics simulations reproduce the observed selectivity and offer insights into the mechanism. When the alkoxy g DFC mechanophore is derived from poly(dihydrofuran), the α-fluoro allyl ether undergoes subsequent hydrolysis to generate 1 equiv of HF and cleave the polymer chain. The hydrolysis is accelerated via acid catalysis, leading to self-amplifying HF generation and concomitant polymer degradation. The mechanically generated HF can be used in combination with fluoride indicators to generate an optical response and to degrade polybutadiene with embedded HF-cleavable silyl ethers (11 mol %). The alkoxy- g DFC mechanophore thus provides a mechanically coupled mechanism of releasing HF for polymer remodeling pathways that complements previous thermally driven mechanisms.

We study active-site models of nonheme iron hydroxylases and their vanadium-based mimics using density functional theory to determine if vanadyl is a faithful structural mimic. We identify crucial structural and energetic differences between ferryl and vanadyl isomers owing to the differences in their ground electronic states, i.e., high spin (HS) for Fe and low spin (LS) for V. For the succinate cofactor bound to the ferryl intermediate, we predict facile interconversion between monodentate and bidentate coordination isomers for ferryl species but difficult rearrangement for vanadyl mimics. We study isomerization of the oxo intermediate between axial and equatorial positions and find the ferryl potential energy surface to be characterized by a large barrier of ca. 10 kcal/mol that is completely absent for the vanadyl mimic. This analysis reveals even starker contrasts between Fe and V in hydroxylases than those observed for this metal substitution in nonheme halogenases. Analysis of the relative bond strengths of coordinating carboxylate ligands for Fe and V reveals that all of the ligands show stronger binding to V than Fe owing to the LS ground state of V in contrast to the HS ground state of Fe, highlighting the limitations of vanadyl mimics of native nonheme iron hydroxylases.

Metalloenzymes catalyze a wide range of chemical transformations, with the active site residues playing a key role in modulating chemical reactivity and selectivity. Unlike smaller synthetic catalysts, a metalloenzyme active site is embedded in a larger protein, which makes interrogation of electronic properties and geometric features with quantum mechanical calculations challenging. Here we implement the ability to fetch crystallographic structures from the Protein Data Bank and analyze the metal binding sites in the program molSimplify. We show the usefulness of the newly created protein3D class to extract the local environment around non-heme iron enzymes containing a two histidine motif and prepare 372 structures for quantum mechanical calculations. Our implementation of protein3D serves to expand the range of systems molSimplify can be used to analyze and will enable high-throughput study of metal-containing active sites in proteins., (© 2023 The Authors. Journal of Computational Chemistry published by Wiley Periodicals LLC.)

Spin-crossover (SCO) complexes are materials that exhibit changes in the spin state in response to external stimuli, with potential applications in molecular electronics. It is challenging to know a priori how to design ligands to achieve the delicate balance of entropic and enthalpic contributions needed to tailor a transition temperature close to room temperature. We leverage the SCO complexes from the previously curated SCO-95 data set [Vennelakanti et al. J. Chem. Phys . 159, 024120 ( 2023 )] to train three machine learning (ML) models for transition temperature ( T 1/2 ) prediction using graph-based revised autocorrelations as features. We perform feature selection using random forest-ranked recursive feature addition (RF-RFA) to identify the features essential to model transferability. Of the ML models considered, the full feature set RF and recursive feature addition RF models perform best, achieving moderate correlation to experimental T 1/2 values. We then compare ML T 1/2 predictions to those from three previously identified best-performing density functional approximations (DFAs) which accurately predict SCO behavior across SCO-95, finding that the ML models predict T 1/2 more accurately than the best-performing DFAs. In addition, we study ML model predictions for a set of 18 SCO complexes for which only estimated T 1/2 values are available. Upon excluding outliers from this set, the RF-RFA RF model shows a strong correlation to estimated T 1/2 values with a Pearson's r of 0.82. In contrast, DFA-predicted T 1/2 values have large errors and show no correlation to estimated T 1/2 values over the same set of complexes. Overall, our study demonstrates slightly superior performance of ML models in comparison with some of the best-performing DFAs, and we expect ML models to improve further as larger data sets of SCO complexes are curated and become available for model training.

With the increasingly more important role of machine learning (ML) models in chemical research, the need for putting a level of confidence to the model predictions naturally arises. Several methods for obtaining uncertainty estimates have been proposed in recent years but consensus on the evaluation of these have yet to be established and different studies on uncertainties generally uses different metrics to evaluate them. We compare three of the most popular validation metrics (Spearman's rank correlation coefficient, the negative log likelihood (NLL) and the miscalibration area) to the error-based calibration introduced by Levi et al. (Sensors 2022, 22, 5540). Importantly, metrics such as the negative log likelihood (NLL) and Spearman's rank correlation coefficient bear little information in themselves. We therefore introduce reference values obtained through errors simulated directly from the uncertainty distribution. The different metrics target different properties and we show how to interpret them, but we generally find the best overall validation to be done based on the error-based calibration plot introduced by Levi et al. Finally, we illustrate the sensitivity of ranking-based methods (e.g. Spearman's rank correlation coefficient) towards test set design by using the same toy model ferent test sets and obtaining vastly different metrics (0.05 vs. 0.65)., (© 2023. The Author(s).)

The design of ion-selective materials with improved separation efficacy and efficiency is paramount, as current technologies fail to meet real-world deployment challenges. Selectivity in these materials can be informed by local ion binding in confined membrane ion channels. In this study, we utilize a data-driven approach to investigate design features in small molecular complexes coordinating ions as simplified models of ion channels. We curate a data set of 563 alkali metal coordinating molecular complexes (i.e., with Li + , Na + , or K + ) from the Cambridge Structural Database and calculate differential ion binding energies using density functional theory. Using this information, we probe when and why structures favor exchange with alternate ions. Our analysis reveals that energetic preferences are related to ion size but are largely due to chemical interactions rather than structural reorganization. We identify unique trends in the selectivity for Li + over other alkali ions, including the presence of N coordination atoms, planar coordination geometry, and small coordinating ring sizes. We use machine learning models to identify the key contributions of both geometric and electronic features in predicting selective ion binding. These physical insights offer preliminary guidance into the design of optimal membranes for ion selectivity.

Hemilabile ligands have the capacity to partially disengage from a metal center, providing a strategy to balance stability and reactivity in catalysis, but they are not straightforward to identify. We identify ligands in the Cambridge Structural Database that have been crystallized with distinct denticities and are thus identifiable as hemilabile ligands. We implement a semi-supervised learning approach using a label-spreading algorithm to augment a small negative set that is supported by heuristic rules of ligand and metal co-occurrence. We show that a heuristic based on coordinating atom identity alone is not sufficient to identify whether a ligand is hemilabile, and our trained machine-learning classification models are instead needed to predict whether a bi-, tri-, or tetradentate ligand is hemilabile with high accuracy and precision. Feature importance analysis of our models shows that the second, third, and fourth coordination spheres all play important roles in ligand hemilability.

Selective halogenation is necessary for a range of fine chemical applications, including the development of therapeutic drugs. While synthetic processes to achieve C-H halogenation require harsh conditions, enzymes such as nonheme iron halogenases carry out some types of C-H halogenation, i.e., chlorination or bromination, with ease, while others, i.e., fluorination, have never been observed in natural or engineered nonheme iron enzymes. Using density functional theory and correlated wave function theory, we investigate the differences in structural and energetic preferences of the smaller fluoride and the larger chloride or bromide intermediates throughout the catalytic cycle. Although we find that the energetics of rate-limiting hydrogen atom transfer are not strongly impacted by fluoride substitution, the higher barriers observed during the radical rebound reaction for fluoride relative to chloride and bromide contribute to the difficulty of C-H fluorination. We also investigate the possibility of isomerization playing a role in differences in reaction selectivity, and our calculations reveal crucial differences in terms of isomer energetics of the key ferryl intermediate between fluoride and chloride/bromide intermediates. While formation of monodentate isomers believed to be involved in selective catalysis is shown for chloride and bromide intermediates, we find that formation of the fluoride monodentate intermediate is not possible in our calculations, which lack additional stabilizing interactions with the greater protein environment. Furthermore, the shorter Fe-F bonds are found to increase isomerization reaction barriers, suggesting that incorporation of residues that form a halogen bond with F and elongate Fe-F bonds could make selective C-H fluorination possible in nonheme iron halogenases. Our work highlights the differences between the fluoride and chloride/bromide intermediates and suggests potential steps toward engineering nonheme iron halogenases to enable selective C-H fluorination.

Transition state search is key in chemistry for elucidating reaction mechanisms and exploring reaction networks. The search for accurate 3D transition state structures, however, requires numerous computationally intensive quantum chemistry calculations due to the complexity of potential energy surfaces. Here we developed an object-aware SE(3) equivariant diffusion model that satisfies all physical symmetries and constraints for generating sets of structures-reactant, transition state and product-in an elementary reaction. Provided reactant and product, this model generates a transition state structure in seconds instead of hours, which is typically required when performing quantum-chemistry-based optimizations. The generated transition state structures achieve a median of 0.08 Å root mean square deviation compared to the true transition state. With a confidence scoring model for uncertainty quantification, we approach an accuracy required for reaction barrier estimation (2.6 kcal mol -1 ) by only performing quantum chemistry-based optimizations on 14% of the most challenging reactions. We envision usefulness for our approach in constructing large reaction networks with unknown mechanisms., (© 2023. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Ab initio multi-reference configuration interaction (MRCI) and coupled cluster singles doubles and perturbative triples [CCSD(T)] levels of theory were used to study ground and excited electronic states of HfCO. We report potential energy curves, dissociation energies ( D e ), excitation energies, harmonic vibrational frequencies, and chemical bonding patterns of HfCO. The 3 Σ - ground state of HfCO has an 1σ 2 2σ 2 1π 2 electron configuration and a ∼30 kcal mol -1 dissociation energy with respect to its lowest-energy fragments Hf( 3 F) + CO(X 1 Σ + ). We further evaluated the D e of its isovalent HfCX (X = S, Se, Te, Po) series and observed that they increase linearly from the lighter HfCO to the heavier HfCPo with the dipole moment of the CX ligand. The same linear relationship was observed for TiCX and ZrCX. We utilized the CCSD(T) benchmark values of D e , excitation energy, and ionization energy (IE) values to evaluate density functional theory (DFT) errors with 23 exchange-correlation functionals spanning GGA, meta-GGA, global GGA hybrid, meta-GGA hybrid, range-separated hybrid, and double-hybrid functional families. The global GGA hybrid B3LYP and range-separated hybrid ωB97X performed well at representing the ground state properties of HfCO ( i.e. , D e and IE). Finally, we extended our DFT analysis to the interaction of a CO molecule with a Hf surface and observed that the surface chemisorption energy and the gas-phase molecular dissociation energy are very similar for some DFAs but not others, suggesting moderate transferability of the benchmarks on these molecules to the solid state.

Metal-organic cages/polyhedra (MOCs) are versatile building blocks for advanced polymer networks with properties that synergistically blend those of traditional polymers and crystalline frameworks. Nevertheless, constructing polyMOCs from very stable Pt(II)-based MOCs or mixtures of metal ions such as Pd(II) and Pt(II) has not, to our knowledge, been demonstrated, nor has exploration of how the dynamics of metal-ligand exchange at the MOC level may impact bulk polyMOC energy dissipation. Here, we introduce a new class of polymer metal-organic cage (polyMOC) gels featuring polyethylene glycol (PEG) strands of varied length cross-linked through bis-pyridyl-carbazole-based M 6 L 12 cubes, where M is Pd(II), Pt(II), or mixtures thereof. We show that, while polyMOCs with varied Pd(II) content have similar network structures, their average stress-relaxation rates are tunable over 3 orders of magnitude due to differences in Pd(II)- and Pt(II)-ligand exchange rates at the M 6 L 12 junction level. Moreover, mixed-metal polyMOCs display relaxation times indicative of intrajunction cooperative interactions, which stands in contrast to previous materials based on point metal junctions. Altogether, this work (1) introduces a novel MOC architecture for polyMOC design, (2) shows that polyMOCs can be prepared from mixtures of Pd(II)/Pt(II), and (3) demonstrates that polyMOCs display unique relaxation behavior due to their multivalent junctions, offering a strategy for controlling polyMOC properties independently of their polymer components.

Dithioalkylidenes are a newly-developed class of conjugate acceptors that undergo thiol exchange via an associative mechanism, enabling decoupling of key material properties for sustainability, biomedical, and sensing applications. Here, we show that the exchange rate is highly sensitive to the structure of the acceptor and tunable over four orders of magnitude in aqueous environments. Cyclic acceptors exchange rapidly, from 0.95 to 15.6 M -1 s -1 , while acyclic acceptors exchange between 3.77x10 -3 and 2.17x10 -2 M -1 s -1 . Computational, spectroscopic, and structural data suggest that cyclic acceptors are more reactive than their acyclic counterparts because of resonance stabilization of the tetrahedral exchange intermediate. We parametrize molecular reactivity with respect to computed descriptors of the electrophilic site and leverage this insight to design a compound with intermediate characteristics. Lastly, we incorporate this dynamic bond into hydrogels and demonstrate that the characteristic stress relaxation time ( τ ) is directly proportional to molecular k ex ., Competing Interests: DECLARATION OF INTERESTS The authors declare no competing interests.

The scission of the O-O bond in O 2 during respiration and the formation of the O-O bond during photosynthesis are the engines of aerobic life. Likewise, the reduction of O 2 and the oxidation of reduced oxygen species to form O 2 are indispensable components for emerging renewable technologies, including energy storage and conversion, yet discrete molecule-like systems that promote these fundamental reactions are rare. Herein, we report a square-planar tetramanganese cluster formed by self-assembly within a metal-organic framework that reversibly reduces O 2 by four electrons, facilitating the interconversion between molecular O 2 and metal-oxo species. The tetranuclear cluster spontaneously cleaves the O-O bond of O 2 at room temperature to generate a tetramanganese-bis(μ 2 -oxo) species, which, in turn, is competent for O-O bond reformation and O 2 evolution at elevated temperatures, enabled by the head-to-head orientation of two oxo species. This study demonstrates the viability of four-electron interconversion between molecular O 2 and metal-oxo species and highlights the importance of site isolation for achieving multi-electron chemistry at polynuclear metal clusters.

Spin crossover (SCO) complexes, which exhibit changes in spin state in response to external stimuli, have applications in molecular electronics and are challenging materials for computational design. We curate a dataset of 95 Fe(II) SCO complexes (SCO-95) from the Cambridge Structural Database that have available low- and high-temperature crystal structures and, in most cases, confirmed experimental spin transition temperatures (T1/2). We study these complexes using density functional theory (DFT) with 30 functionals spanning across multiple rungs of "Jacob's ladder" to understand the effect of exchange-correlation functional on electronic and Gibbs free energies associated with spin crossover. We specifically assess the effect of varying the Hartree-Fock exchange fraction (aHF) in structures and properties within the B3LYP family of functionals. We identify three best-performing functionals, a modified version of B3LYP (aHF = 0.10), M06-L, and TPSSh, that accurately predict SCO behavior for the majority of the complexes. While M06-L performs well, MN15-L, a more recently developed Minnesota functional, fails to predict SCO behavior for all complexes, which could be the result of differences in datasets used for parametrization of M06-L and MN15-L and also the increased number of parameters for MN15-L. Contrary to observations from prior studies, double-hybrids with higher aHF values are found to strongly stabilize high-spin states and therefore exhibit poor performance in predicting SCO behavior. Computationally predicted T1/2 values are consistent among the three functionals but show limited correlation to experimentally reported T1/2 values. These failures are attributed to the lack of crystal packing effects and counter-anions in the DFT calculations that would be needed to account for phenomena such as hysteresis and two-step SCO behavior. The SCO-95 set thus presents opportunities for method development, both in terms of increasing model complexity and method fidelity., (© 2023 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).)

The challenge of direct partial oxidation of methane to methanol has motivated the targeted search of metal-organic frameworks (MOFs) as a promising class of materials for this transformation because of their site-isolated metals with tunable ligand environments. Thousands of MOFs have been synthesized, yet relatively few have been screened for their promise in methane conversion. We developed a high-throughput virtual screening workflow that identifies MOFs from a diverse space of experimental MOFs that have not been studied for catalysis, yet are thermally stable, synthesizable, and have promising unsaturated metal sites for C-H activation via a terminal metal-oxo species. We carried out density functional theory calculations of the radical rebound mechanism for methane-to-methanol conversion on models of the secondary building units (SBUs) from 87 selected MOFs. While we showed that oxo formation favorability decreases with increasing 3d filling, consistent with prior work, previously observed scaling relations between oxo formation and hydrogen atom transfer (HAT) are disrupted by the greater diversity in our MOF set. Accordingly, we focused on Mn MOFs, which favor oxo intermediates without disfavoring HAT or leading to high methanol release energies─a key feature for methane hydroxylation activity. We identified three Mn MOFs comprising unsaturated Mn centers bound to weak-field carboxylate ligands in planar or bent geometries with promising methane-to-methanol kinetics and thermodynamics. The energetic spans of these MOFs are indicative of promising turnover frequencies for methane to methanol that warrant further experimental catalytic studies.

We survey more than 240 000 crystallized mononuclear transition metal complexes (TMCs) to identify trends in preferred geometric structure and metal coordination. While we observe that an increased level of d filling correlates with a lower coordination number preference, we note exceptions, and we observe undersampling of 4d/5d transition metals and 3p-coordinating ligands. For the one-third of mononuclear TMCs that are octahedral, analysis of the 67 symmetry classes of their ligand environments reveals that complexes often contain monodentate ligands that may be removable, forming an open site amenable to catalysis. Due to their use in catalysis, we analyze trends in coordination by tetradentate ligands in terms of the capacity to support multiple metals and the variability of coordination geometry. We identify promising tetradentate ligands that co-occur in crystallized complexes with labile monodentate ligands that would lead to reactive sites. Literature mining suggests that these ligands are untapped as catalysts, motivating proposal of a promising octa-functionalized porphyrin.

The cis - and trans -isomers of a silacycloheptene were selectively synthesized by the alkylation of a silyl dianion, a novel approach to strained cycloalkenes. The trans -silacycloheptene ( trans -SiCH) was significantly more strained than the cis isomer, as predicted by quantum chemical calculations and confirmed by crystallographic signatures of a twisted alkene. Each isomer exhibited distinct reactivity toward ring-opening metathesis polymerization (ROMP), where only trans -SiCH afforded high-molar-mass polymer under enthalpy-driven ROMP. Hypothesizing that the introduction of silicon might result in increased molecular compliance at large extensions, we compared poly( trans -SiCH) to organic polymers by single-molecule force spectroscopy (SMFS). Force-extension curves from SMFS showed that poly( trans -SiCH) is more easily overstretched than two carbon-based analogues, polycyclooctene and polybutadiene, with stretching constants that agree well with the results of computational simulations.

Confined fluids and electrolyte solutions in nanopores exhibit rich and surprising physics and chemistry that impact the mass transport and energy efficiency in many important natural systems and industrial applications. Existing theories often fail to predict the exotic effects observed in the narrowest of such pores, called single-digit nanopores (SDNs), which have diameters or conduit widths of less than 10 nm, and have only recently become accessible for experimental measurements. What SDNs reveal has been surprising, including a rapidly increasing number of examples such as extraordinarily fast water transport, distorted fluid-phase boundaries, strong ion-correlation and quantum effects, and dielectric anomalies that are not observed in larger pores. Exploiting these effects presents myriad opportunities in both basic and applied research that stand to impact a host of new technologies at the water-energy nexus, from new membranes for precise separations and water purification to new gas permeable materials for water electrolyzers and energy-storage devices. SDNs also present unique opportunities to achieve ultrasensitive and selective chemical sensing at the single-ion and single-molecule limit. In this review article, we summarize the progress on nanofluidics of SDNs, with a focus on the confinement effects that arise in these extremely narrow nanopores. The recent development of precision model systems, transformative experimental tools, and multiscale theories that have played enabling roles in advancing this frontier are reviewed. We also identify new knowledge gaps in our understanding of nanofluidic transport and provide an outlook for the future challenges and opportunities at this rapidly advancing frontier.

Virtual high-throughput screening (VHTS) and machine learning (ML) with density functional theory (DFT) suffer from inaccuracies from the underlying density functional approximation (DFA). Many of these inaccuracies can be traced to the lack of derivative discontinuity that leads to a curvature in the energy with electron addition or removal. Over a dataset of nearly one thousand transition metal complexes typical of VHTS applications, we computed and analyzed the average curvature ( i.e. , deviation from piecewise linearity) for 23 density functional approximations spanning multiple rungs of "Jacob's ladder". While we observe the expected dependence of the curvatures on Hartree-Fock exchange, we note limited correlation of curvature values between different rungs of "Jacob's ladder". We train ML models ( i.e. , artificial neural networks or ANNs) to predict the curvature and the associated frontier orbital energies for each of these 23 functionals and then interpret differences in curvature among the different DFAs through analysis of the ML models. Notably, we observe spin to play a much more important role in determining the curvature of range-separated and double hybrids in comparison to semi-local functionals, explaining why curvature values are weakly correlated between these and other families of functionals. Over a space of 187.2k hypothetical compounds, we use our ANNs to pinpoint DFAs for which representative transition metal complexes have near-zero curvature with low uncertainty, demonstrating an approach to accelerate screening of complexes with targeted optical gaps.

Organic-inorganic hybrid materials present new opportunities for creating low-dimensional structures with unique light-matter interaction. In this work, we report a chemically robust yellow emissive one-dimensional (1D) semiconductor, silver 2,6-difluorophenylselenolate─AgSePhF 2 (2,6), a new member of the broader class of hybrid low-dimensional semiconductors, metal-organic chalcogenolates. While silver phenylselenolate (AgSePh) crystallizes as a two-dimensional (2D) van der Waals semiconductor, introduction of fluorine atoms at the (2,6) position of the phenyl ring induces a structural transition from 2D sheets to 1D chains. Density functional theory calculations reveal that AgSePhF 2 (2,6) has strongly dispersive conduction and valence bands along the 1D crystal axis. Visible photoluminescence centered around λ p ≈ 570 nm at room temperature exhibits both prompt (110 ps) and delayed (36 ns) components. The absorption spectrum exhibits excitonic resonances characteristic of low-dimensional hybrid semiconductors, with an exciton binding energy of approximately 170 meV as determined by temperature-dependent photoluminescence. The discovery of an emissive 1D silver organoselenolate highlights the structural and compositional richness of the chalcogenolate material family and provides new insights for molecular engineering of low-dimensional hybrid organic-inorganic semiconductors.

Plants contain rapidly evolving specialized enzymes that support the biosynthesis of functionally diverse natural products. In coumarin biosynthesis, a BAHD acyltransferase-family enzyme COSY was recently discovered to accelerate coumarin formation as the only known BAHD enzyme to catalyze an intramolecular acyl transfer reaction. Here we investigate the structural and mechanistic basis for COSY's coumarin synthase activity. Our structural analyses reveal an unconventional active-site configuration adapted to COSY's specialized activity. Through mutagenesis studies and deuterium exchange experiments, we identify a unique proton exchange mechanism at the α-carbon of the o-hydroxylated trans-hydroxycinnamoyl-CoA substrates during the catalytic cycle of COSY. Quantum mechanical cluster modeling and molecular dynamics further support this key mechanism for lowering the activation energy of the rate-limiting trans-to-cis isomerization step in coumarin production. This study unveils an unconventional catalytic mechanism mediated by a BAHD-family enzyme, and sheds light on COSY's evolutionary origin and its recruitment to coumarin biosynthesis in eudicots., (© 2023. The Author(s).)

Convenient strategies for the deconstruction and reprocessing of thermosets could improve the circularity of these materials, but most approaches developed to date do not involve established, high-performance engineering materials. Here, we show that bifunctional silyl ether, i.e., R'O-SiR 2 -OR'', (BSE)-based comonomers generate covalent adaptable network analogues of the industrial thermoset polydicyclopentadiene (pDCPD) through a novel BSE exchange process facilitated by the low-cost food-safe catalyst octanoic acid. Experimental studies and density functional theory calculations suggest an exchange mechanism involving silyl ester intermediates with formation rates that strongly depend on the Si-R 2 substituents. As a result, pDCPD thermosets manufactured with BSE comonomers display temperature- and time-dependent stress relaxation as a function of their substituents. Moreover, bulk remolding of pDCPD thermosets is enabled for the first time. Altogether, this work presents a new approach toward the installation of exchangeable bonds into commercial thermosets and establishes acid-catalyzed BSE exchange as a versatile addition to the toolbox of dynamic covalent chemistry.

When a many-body wave function of a system cannot be captured by a single determinant, high-level multireference (MR) methods are required to properly explain its electronic structure. MR diagnostics to estimate the magnitude of such static correlation have been primarily developed for molecular systems and range from low in computational cost to as costly as the full MR calculation itself. We report the first application of low-cost MR diagnostics based on the fractional occupation number calculated with finite-temperature DFT to solid-state systems. To compare the behavior of the diagnostics on solids and molecules, we select metal-organic frameworks (MOFs) as model materials because their reticular nature provides an intuitive way to identify molecular derivatives. On a series of closed-shell MOFs, we demonstrate that the DFT-based MR diagnostics are equally applicable to solids as to their molecular derivatives. The magnitude and spatial distribution of the MR character of a MOF are found to have a good correlation with those of its molecular derivatives, which can be calculated much more affordably in comparison to those of the full MOF. The additivity of MR character discussed here suggests the set of molecular derivatives to be a good representation of a MOF for both MR detection and ultimately for MR corrections, facilitating accurate and efficient high-throughput screening of MOFs and other porous solids.

Prediction of the excited state properties of photoactive iridium complexes challenges ab initio methods such as time-dependent density functional theory (TDDFT) both from the perspective of accuracy and of computational cost, complicating high-throughput virtual screening (HTVS). We instead leverage low-cost machine learning (ML) models and experimental data for 1380 iridium complexes to perform these prediction tasks. We find the best-performing and most transferable models to be those trained on electronic structure features from low-cost density functional tight binding calculations. Using artificial neural network (ANN) models, we predict the mean emission energy of phosphorescence, the excited state lifetime, and the emission spectral integral for iridium complexes with accuracy competitive with or superseding that of TDDFT. We conduct feature importance analysis to determine that high cyclometalating ligand ionization potential correlates to high mean emission energy, while high ancillary ligand ionization potential correlates to low lifetime and low spectral integral. As a demonstration of how our ML models can be used for HTVS and the acceleration of chemical discovery, we curate a set of novel hypothetical iridium complexes and use uncertainty-controlled predictions to identify promising ligands for the design of new phosphors while retaining confidence in the quality of the ANN predictions., Competing Interests: The authors declare no competing financial interest., (This journal is © The Royal Society of Chemistry.)

Approximate density functional theory has become indispensable owing to its balanced cost-accuracy trade-off, including in large-scale screening. To date, however, no density functional approximation (DFA) with universal accuracy has been identified, leading to uncertainty in the quality of data generated from density functional theory. With electron density fitting and Δ-learning, we build a DFA recommender that selects the DFA with the lowest expected error with respect to the gold standard (but cost-prohibitive) coupled cluster theory in a system-specific manner. We demonstrate this recommender approach on the evaluation of vertical spin splitting energies of transition metal complexes. Our recommender predicts top-performing DFAs and yields excellent accuracy (about 2 kcal mol -1 ) for chemical discovery, outperforming both individual Δ-learning models and the best conventional single-functional approach from a set of 48 DFAs. By demonstrating transferability to diverse synthesized compounds, our recommender potentially addresses the accuracy versus scope dilemma broadly encountered in computational chemistry., (© 2022. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Silver phenylselenolate (AgSePh, also known as "mithrene") and silver phenyltellurolate (AgTePh, also known as "tethrene") are two-dimensional (2D) van der Waals semiconductors belonging to an emerging class of hybrid organic-inorganic materials called metal-organic chalcogenolates. Despite having the same crystal structure, AgSePh and AgTePh exhibit a strikingly different excitonic behavior. Whereas AgSePh exhibits narrow, fast luminescence with a minimal Stokes shift, AgTePh exhibits comparatively slow luminescence that is significantly broadened and red-shifted from its absorption minimum. Using time-resolved and temperature-dependent absorption and emission microspectroscopy, combined with subgap photoexcitation studies, we show that exciton dynamics in AgTePh films are dominated by an intrinsic self-trapping behavior, whereas dynamics in AgSePh films are dominated by the interaction of band-edge excitons with a finite number of extrinsic defect/trap states. Density functional theory calculations reveal that AgSePh has simple parabolic band edges with a direct gap at Γ, whereas AgTePh has a saddle point at Γ with a horizontal splitting along the Γ-N 1 direction. The correlation between the unique band structure of AgTePh and exciton self-trapping behavior is unclear, prompting further exploration of excitonic phenomena in this emerging class of hybrid 2D semiconductors.

Transition-metal chromophores with earth-abundant transition metals are an important design target for their applications in lighting and nontoxic bioimaging, but their design is challenged by the scarcity of complexes that simultaneously have well-defined ground states and optimal target absorption energies in the visible region. Machine learning (ML) accelerated discovery could overcome such challenges by enabling the screening of a larger space but is limited by the fidelity of the data used in ML model training, which is typically from a single approximate density functional. To address this limitation, we search for consensus in predictions among 23 density functional approximations across multiple rungs of "Jacob's ladder". To accelerate the discovery of complexes with absorption energies in the visible region while minimizing the effect of low-lying excited states, we use two-dimensional (2D)efficient global optimization to sample candidate low-spin chromophores from multimillion complex spaces. Despite the scarcity (i.e., ∼0.01%) of potential chromophores in this large chemical space, we identify candidates with high likelihood (i.e., >10%) of computational validation as the ML models improve during active learning, representing a 1000-fold acceleration in discovery. Absorption spectra of promising chromophores from time-dependent density functional theory verify that 2/3 of candidates have the desired excited-state properties. The observation that constituent ligands from our leads have demonstrated interesting optical properties in the literature exemplifies the effectiveness of our construction of a realistic design space and active learning approach., Competing Interests: The authors declare no competing financial interest., (© 2022 The Authors. Published by American Chemical Society.)

To accelerate the exploration of chemical space, it is necessary to identify the compounds that will provide the most additional information or value. A large-scale analysis of mononuclear octahedral transition metal complexes deposited in an experimental database confirms an under-representation of lower-symmetry complexes. From a set of around 1000 previously studied Fe(II) complexes, we show that the theoretical space of synthetically accessible complexes formed from the relatively small number of unique ligands is significantly (∼816k) larger. For the properties of these complexes, we validate the concept of ligand additivity by inferring heteroleptic properties from a stoichiometric combination of homoleptic complexes. An improved interpolation scheme that incorporates information about cis and trans isomer effects predicts the adiabatic spin-splitting energy to around 2 kcal/mol and the HOMO level to less than 0.2 eV. We demonstrate a multi-stage strategy to discover leads from the 816k Fe(II) complexes within a targeted property region. We carry out a coarse interpolation from homoleptic complexes that we refine over a subspace of ligands based on the likelihood of generating complexes with targeted properties. We validate our approach on nine new binary and ternary complexes predicted to be in a targeted zone of discovery, suggesting opportunities for efficient transition metal complex discovery.

Mononuclear non-heme Fe(II)- and α-ketoglutarate-dependent oxygenases (FeDOs) catalyze a site-selective C-H hydroxylation. Variants of these enzymes in which a conserved Asp/Glu residue in the Fe(II)-binding facial triad is replaced by Ala/Gly can, in some cases, bind various anionic ligands and catalyze non-native chlorination and bromination reactions. In this study, we explore the binding of different anions to an FeDO facial triad variant, SadX, and the effects of that binding on HO • vs X • rebound. We establish not only that chloride and bromide enable non-native halogenation reactions but also that all anions investigated, including azide, cyanate, formate, and fluoride, significantly accelerate and influence the site selectivity of SadX hydroxylation catalysis. Azide and cyanate also lead to the formation of products resulting from N 3 • , NCO • , and OCN • rebound. While fluoride rebound is not observed, the rate acceleration provided by this ligand leads us to calculate barriers for HO • and F • rebound from a putative Fe(III)(OH)(F) intermediate. These calculations suggest that the lack of fluorination is due to the relative barriers of the HO • and F • rebound transition states rather than an inaccessible barrier for F • rebound. Together, these results improve our understanding of the FeDO facial triad variant tolerance of different anionic ligands, their ability to promote rebound involving these ligands, and inherent rebound preferences relative to HO • that will aid efforts to develop non-native catalysis using these enzymes.

Accurate virtual high-throughput screening (VHTS) of transition metal complexes (TMCs) remains challenging due to the possibility of high multireference (MR) character that complicates property evaluation. We compute MR diagnostics for over 5,000 ligands present in previously synthesized octahedral mononuclear transition metal complexes in the Cambridge Structural Database (CSD). To accomplish this task, we introduce an iterative approach for consistent ligand charge assignment for ligands in the CSD. Across this set, we observe that the MR character correlates linearly with the inverse value of the averaged bond order over all bonds in the molecule. We then demonstrate that ligand additivity of the MR character holds in TMCs, which suggests that the TMC MR character can be inferred from the sum of the MR character of the ligands. Encouraged by this observation, we leverage ligand additivity and develop a ligand-derived machine learning representation to train neural networks to predict the MR character of TMCs from properties of the constituent ligands. This approach yields models with excellent performance and superior transferability to unseen ligand chemistry and compositions.

The immobilization of homogeneous catalysts onto supports to improve recyclability while maintaining catalytic efficiency is often a trial-and-error process limited by poor control of the local catalyst environment and few strategies to append catalysts to support materials. Here, we introduce a modular heterogenous catalysis platform that addresses these challenges. Our approach leverages the well-defined interiors of self-assembled Pd 12 L 24 metal-organic cages/polyhedra (MOCs): simple mixing of a catalyst-ligand of choice with a polymeric ligand, spacer ligands, and a Pd salt induces self-assembly of Pd 12 L 24 -cross-linked polymer gels featuring endohedrally catalyst-functionalized junctions. Semi-empirical calculations show that catalyst incorporation into the MOC junctions of these materials has minimal affect on the MOC geometry, giving rise to well-defined nanoconfined catalyst domains as confirmed experimentally using several techniques. Given the unique network topology of these freestanding gels, they are mechanically robust regardless of their endohedral catalyst composition, allowing them to be physically manipulated and transferred from one reaction to another to achieve multiple rounds of catalysis. Moreover, by decoupling the catalyst environment (interior of MOC junctions) from the physical properties of the support (the polymer matrix), this strategy enables catalysis in environments where homogeneous catalyst analogues are not viable, as demonstrated for the Au(I)-catalyzed cyclization of 4-pentynoic acid in aqueous media.

Virtual high-throughput screening (VHTS) and machine learning (ML) have greatly accelerated the design of single-site transition-metal catalysts. VHTS of catalysts, however, is often accompanied with a high calculation failure rate and wasted computational resources due to the difficulty of simultaneously converging all mechanistically relevant reactive intermediates to expected geometries and electronic states. We demonstrate a dynamic classifier approach, i.e., a convolutional neural network that monitors geometry optimizations on the fly, and exploit its good performance and transferability in identifying geometry optimization failures for catalyst design. We show that the dynamic classifier performs well on all reactive intermediates in the representative catalytic cycle of the radical rebound mechanism for the conversion of methane to methanol despite being trained on only one reactive intermediate. The dynamic classifier also generalizes to chemically distinct intermediates and metal centers absent from the training data without loss of accuracy or model confidence. We rationalize this superior model transferability as arising from the use of electronic structure and geometric information generated on-the-fly from density functional theory calculations and the convolutional layer in the dynamic classifier. When used in combination with uncertainty quantification, the dynamic classifier saves more than half of the computational resources that would have been wasted on unsuccessful calculations for all reactive intermediates being considered.

Gut microbial decarboxylation of amino acid-derived arylacetates is a chemically challenging enzymatic transformation which generates small molecules that impact host physiology. The glycyl radical enzyme (GRE) indoleacetate decarboxylase from Olsenella uli ( Ou IAD) performs the non-oxidative radical decarboxylation of indole-3-acetate (I3A) to yield skatole, a disease-associated metabolite produced in the guts of swine and ruminants. Despite the importance of IAD, our understanding of its mechanism is limited. Here, we characterize the mechanism of Ou IAD, evaluating previously proposed hypotheses of: (1) a Kolbe-type decarboxylation reaction involving an initial 1-e - oxidation of the carboxylate of I3A or (2) a hydrogen atom abstraction from the α-carbon of I3A to generate an initial carbon-centered radical. Site-directed mutagenesis, kinetic isotope effect experiments, analysis of reactions performed in D 2 O, and computational modeling are consistent with a mechanism involving initial hydrogen atom transfer. This finding expands the types of radical mechanisms employed by GRE decarboxylases and non-oxidative decarboxylases, more broadly. Elucidating the mechanism of IAD decarboxylation enhances our understanding of radical enzymes and may inform downstream efforts to modulate this disease-associated metabolism.

Machine learning (ML) has become a part of the fabric of high-throughput screening and computational discovery of materials. Despite its increasingly central role, challenges remain in fully realizing the promise of ML. This is especially true for the practical acceleration of the engineering of robust materials and the development of design strategies that surpass trial and error or high-throughput screening alone. Depending on the quantity being predicted and the experimental data available, ML can either outperform physics-based models, be used to accelerate such models, or be integrated with them to improve their performance. We cover recent advances in algorithms and in their application that are starting to make inroads toward ( a ) the discovery of new materials through large-scale enumerative screening, ( b ) the design of materials through identification of rules and principles that govern materials properties, and ( c ) the engineering of practical materials by satisfying multiple objectives. We conclude with opportunities for further advancement to realize ML as a widespread tool for practical computational materials design.

The Mo/W-containing metalloenzyme formate dehydrogenase (FDH) is an efficient and selective natural catalyst that reversibly converts CO 2 to formate under ambient conditions. In this study, we investigate the impact of the greater protein environment on the electrostatic potential (ESP) of the active site. To model the enzyme environment, we used a combination of classical molecular dynamics and multiscale quantum-mechanical (QM)/molecular-mechanical (MM) simulations. We leverage charge shift analysis to systematically construct QM regions and analyze the electronic environment of the active site by evaluating the degree of charge transfer between the core active site and the protein environment. The contribution of the terminal chalcogen ligand to the ESP of the metal center is substantial and dependent on the chalcogen identity, with similar, less negative ESPs for Se and S terminal chalcogens in comparison to O regardless of whether the metal is Mo or W. The orientation of the side chains and conformations of the cofactor also affect the ESP, highlighting the importance of sampling dynamic fluctuations in the protein. Overall, our observations suggest that the terminal chalcogen ligand identity plays an important role in the enzymatic activity of FDH, suggesting opportunities for a rational bioinspired catalyst design.

The predictive accuracy of density functional theory (DFT) is hampered by delocalization errors, especially for correlated systems such as transition-metal complexes. Two complementary strategies have been developed to reduce delocalization error: eliminating the global curvature with change in charge, and applying a linear response Hubbard U as a measure of local curvature at a metal center at fixed charge in a DFT+U framework. We investigate the relationship between the two delocalization error measures as the ligand field strength is varied with the number of strong-field ligands in a series of heteroleptic complexes or by geometrically constraining the metal-ligand bond length in homoleptic octahedral complexes. We show that across these sets of complexes an inverse relationship generally exists between global and local curvatures. We find that effects of ligand substitution on both measures of delocalization are typically additive, but the quantities seldom coincide.

Knowledge of the chemical bonding of HfO and HfB ground and low-lying electronic states provides essential insights into a range of catalysts and materials that contain Hf-O or Hf-B moieties. Here, we carry out high-level multi-reference configuration interaction theory and coupled cluster quantum chemical calculations on these systems. We compute full potential energy curves, excitation energies, ionization energies, electronic configurations, and spectroscopic parameters with large quadruple-ζ and quintuple-ζ quality correlation consistent basis sets. We also investigate equilibrium chemical bonding patterns and effects of correlating core electrons on property predictions. Differences in the ground state electron configuration of HfB(X 4 Σ - ) and HfO(X 1 Σ + ) lead to a significantly stronger bond in HfO than HfB, as judged by both dissociation energies and equilibrium bond distances. We extend our analysis to the chemical bonding patterns of the isovalent HfX (X = O, S, Se, Te, and Po) series and observe similar trends. We also note a linear trend between the decreasing value of the dissociation energy (D e ) from HfO to HfPo and the singlet-triplet energy gap (ΔE S-T ) of the molecule. Finally, we compare these benchmark results to those obtained using density functional theory (DFT) with 23 exchange-correlation functionals spanning multiple rungs of "Jacob's ladder." When comparing DFT errors to coupled cluster reference values on dissociation energies, excitation energies, and ionization energies of HfB and HfO, we observe semi-local generalized gradient approximations to significantly outperform more complex and high-cost functionals.