Lignin is a major component of plant cell walls that is typically underutilized in selective conversion strategies for renewable fuels and chemicals. The mechanisms by which thermal and catalytic treatments deconstruct lignin remain elusive, which is where quantum mechanical calculations can offer fundamental insights. Here, we compute homolytic bond dissociation enthalpies (BDEs) for four prevalent linkages in 69 lignin model compounds, including β-O-4, α-O-4, β-5, and biphenyl bonds, with a large range of natural and oxidized substituents. These calculations include ab initio benchmark values extrapolated to the complete basis set limit and full conformational searches for each compound. The results quantify both the relative BDEs among common lignin bonds and the effect of native and oxidized substituents on the functional groups in lignin. These data yield insights into thermal lignin deconstruction for a large range of prevalent linkages and aid in the identification of targets for catalytic cleavage.

Acid catalysis has long been used to depolymerize plant cell wall polysaccharides, and the mechanisms by which acid affects carbohydrates have been extensively studied. Lignin depolymerization, however, is not as well understood, primarily due to the heterogeneity and reactivity of lignin. We present an experimental and theoretical study of acid-catalyzed cleavage of two non-phenolic and two phenolic dimers that exhibit the β-O-4 ether linkage, the most common intermonomer bond in lignin. This work demonstrates that the rate of acid-catalyzed β-O-4 cleavage in dimers exhibiting a phenolic hydroxyl group is 2 orders of magnitude faster than in non-phenolic dimers. The experiments suggest that the major product distribution is similar for all model compounds, but a stable phenyl-dihydrobenzofuran species is observed in the acidolysis of two of the γ-carbinol containing model compounds. The presence of a methoxy substituent, commonly found in native lignin, prevents the formation of this intermediate. Reaction pathways were examined with quantum mechanical calculations, which aid in explaining the substantial differences in reactivity. Moreover, we use a radical scavenger to show that the commonly proposed homolytic cleavage pathway of phenolic β-O-4 linkages is unlikely in acidolysis conditions. Overall, this study explains the disparity between rates of β-O-4 cleavage seen in model compound experiments and acid pretreatment of biomass, and implies that depolymerization of lignin during acid-catalyzed pretreatment or fractionation will proceed via a heterolytic, unzipping mechanism wherein β-O-4 linkages are cleaved from the phenolic ends of branched, polymer chains inward toward the core of the polymer.

Lytic polysaccharide monooxygenases (LPMOs) exhibit a mononuclear copper-containing active site and use dioxygen and a reducing agent to oxidatively cleave glycosidic linkages in polysaccharides. LPMOs represent a unique paradigm in carbohydrate turnover and exhibit synergy with hydrolytic enzymes in biomass depolymerization. To date, several features of copper binding to LPMOs have been elucidated, but the identity of the reactive oxygen species and the key steps in the oxidative mechanism have not been elucidated. Here, density functional theory calculations are used with an enzyme active site model to identify the reactive oxygen species and compare two hypothesized reaction pathways in LPMOs for hydrogen abstraction and polysaccharide hydroxylation; namely, a mechanism that employs a η 1 -superoxo intermediate, which abstracts a substrate hydrogen and a hydroperoxo species is responsible for substrate hydroxylation, and a mechanism wherein a copper-oxyl radical abstracts a hydrogen and subsequently hydroxylates the substrate via an oxygen-rebound mechanism. The results predict that oxygen binds end-on (η 1 ) to copper, and that a copperoxyl-mediated, oxygen-rebound mechanism is energetically preferred. The N-terminal histidine methylation is also examined, which is thought to modify the structure and reactivity of the enzyme. Density functional theory calculations suggest that this posttranslational modification has only a minor effect on the LPMO active site structure or reactivity for the examined steps. Overall, this study suggests the steps in the LPMO mechanism for oxidative cleavage of glycosidic bonds.

Background: Lytic polysaccharide monooxygenases (LPMOs) represent a recently discovered enzymatic route to cleave carbohydrates. Results:We report the first basidiomycete LPMO structure and describe enzyme-cellulose interactions with simulation. Conclusion:We characterize the copper-containing active site and identify loops important for substrate recognition and binding. Significance: This structure is the first LPMO from a model basidiomycete fungus that contains many LPMO genes.

Research efforts in zeolite catalysis have become increasingly cognizant of the diversity in structure and function resulting from the distribution of framework aluminum atoms, through emerging reports of catalytic phenomena that fall outside those recognizable as the shape-selective ones emblematic of its earlier history. Molecular-level descriptions of how active-site distributions affect catalysis are an aspirational goal articulated frequently in experimental and theoretical research, yet they are limited by imprecise knowledge of the structure and behavior of the zeolite materials under interrogation. In experimental research, higher precision can result from more reliable control of structure during synthesis and from more robust and quantitative structural and kinetic characterization probes. In theoretical research, construction of models with specific aluminum locations and distributions seldom capture the heterogeneity inherent to the materials studied by experiment. In this Perspective, we discuss research findings that appropriately frame the challenges in developing more predictive synthesis−structure−function relations for zeolites, highlighting studies on ZSM-5 zeolites that are among the most structurally complex molecular sieve frameworks and the most widely studied because of their versatility in commercial applications. We discuss research directions to address these challenges and forge stronger connections between zeolite structure, composition, and active sites to catalytic function. Such connections promise to aid in bridging the findings of theoretical and experimental catalysis research, and transforming zeolite active site design from an empirical endeavor into a more predictable science founded on validated models.

Bond dissociation enthalpies (BDEs) of organic molecules play a fundamental role in determining chemical reactivity and selectivity. However, BDE computations at sufficiently high levels of quantum mechanical theory require substantial computing resources. In this paper, we develop a machine learning model capable of accurately predicting BDEs for organic molecules in a fraction of a second. We perform automated density functional theory (DFT) calculations at the M06-2X/def2-TZVP level of theory for 42,577 small organic molecules, resulting in 290,664 BDEs. A graph neural network trained on a subset of these results achieves a mean absolute error of 0.58 kcal mol −1 (vs DFT) for BDEs of unseen molecules. We further demonstrate the model on two applications: first, we rapidly and accurately predict major sites of hydrogen abstraction in the metabolism of drug-like molecules, and second, we determine the dominant molecular fragmentation pathways during soot formation.

Molecular dynamics simulations starting from different initial conditions are commonly used to mimic the behavior of an experimental ensemble. We show in this article that when a Langevin thermostat is used to maintain constant temperature during such simulations, extreme care must be taken when choosing the random number seeds to prevent statistical correlation among the MD trajectories. While recent studies have shown that stochastically thermostatted trajectories evolving within a single potential basin with identical random number seeds tend to synchronize, we show that there is a synchronization effect even for complex, biologically relevant systems. We demonstrate this effect in simulations of alanine trimer and pentamer and in a simulation of a temperature-jump experiment for peptide folding of a 14-residue peptide. Even in replica-exchange simulations, in which the trajectories are at different temperatures, we find partial synchronization occurring when the same random number seed is employed. We explain this by extending the recent derivation of the synchronization effect for two trajectories in a harmonic well to the case in which the trajectories are at two different temperatures. Our results suggest several ways in which mishandling selection of a pseudorandom number generator initial seed can lead to corruption of simulation data. Simulators can fall into this trap in simple situations such as neglecting to specifically indicate different random seeds in either parallel or sequential restart simulations, utilizing a simulation package with a weak pseudorandom number generator, or using an advanced simulation algorithm that has not been programmed to distribute initial seeds.

The Diels-Alder reactions of the cycloalkenes, cyclohexene through cyclopropene, with a series of dienes--1,3-dimethoxybutadiene, cyclopentadiene, 3,6-dimethyltetrazine, and 3,6-bis(trifluoromethyl)tetrazine--were studied with quantum mechanical calculations and compared with experimental values when available. The reactivities of cycloalkenes as dienophiles were found by a distortion/interaction analysis to be distortion controlled. The energies required for cycloalkenes to be distorted into the Diels-Alder transition states increase as the ring size of cycloalkenes increases from cyclopropene to cyclohexene, resulting in an increase in activation barriers. The reactivities of the dienes are controlled by both distortion and interaction energies. In normal Diels-Alder reactions with cycloalkenes, the electron-rich 1,3-dimethoxybutadiene exhibits stronger interaction energies than cyclopentadiene, but the high distortion energies required for 1,3-dimethoxybutadiene to achieve transition-state geometries overtake the favorable interaction, resulting in higher activation barriers. In inverse-electron-demand Diels-Alder reactions of 3,6-dimethyltetrazine and 3,6-bis(trifluoromethyl)tetrazine, the reactivities are mainly controlled by interaction energies.